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

In vivo Validation of Hsp90 Trans-splicing in Giardia lamblia: Highlighting the Role of Cis-elements

Giardia lambliacauses giardiasis, one of the most common human infectious diseases globally. Previous studies from our lab have shown that hsp90 gene ofGiardia is split into two halves, namely hspN and hspC. The independent pre-mRNAs of these split genes join by trans-splicing, producing a full-length Hsp90 (FlHsp90) mRNA. Genetic manipulation of the participating genes is necessary to understand the mechanism and significance of such trans-splicing based expression of Hsp90. In this study, we have performed transfection based exogenous expression of hspN and/or hspC in G. lamblia. We electroporated a plasmid containing the Avi-tagged hspN component of Hsp90 and examined its fate in G. lamblia. We show that the exogenously expressed hspN RNA gets trans-spliced to endogenously expressed hspC RNA, giving rise to a hybrid-FlHsp90. We highlight the importance of cis-elements in this trans-splicing reaction through mutational analysis. The episomal plasmid carrying deletions in the intronic region of hspN, showed inhibition of the trans-splicing reaction.Additionally, exogenous hspC RNA also followed the same fate as of exogenous hspN, while upon co-transfection with episomal hspN, they underwent trans-splicing with each other. Using eGFP as a test protein, we have shown that intronic sequences of hsp90 gene can guide trans-splicing mediated repair of any associated exonic sequences. Our study provides in vivo validation of Hsp90 trans-splicing, showing crucial role of cis-elements and importantly highlights the potential of hsp90 intronic sequences to function as a minimal splicing tool.

A core UPS molecular complement implicates unique endocytic compartments at the parasite-host interface in Giardia lamblia

Unconventional protein secretion (UPS) plays important roles in cell physiology. In contrast to canonical secretory routes, UPS does not generally require secretory signal sequences and often bypasses secretory compartments such as the ER and the Golgi apparatus. Giardia lamblia is a protist parasite with reduced subcellular complexity which releases several proteins, some of them virulence factors, without canonical secretory signals. This implicates UPS at the parasite-host interface. No dedicated machinery nor mechanism(s) for UPS in Giardia are currently known, although speculations on the involvement of endocytic organelles called PV/PECs, have been put forth. To begin to address the question of whether PV/PECs are implicated in virulence-associated UPS and to define the composition of molecular machinery involved in protein release, we employed affinity purification and mass spectrometry, coupled to microscopy-based subcellular localization and signal correlation quantification to investigate the interactomes of 11 reported unconventionally secreted proteins, all predicted to be cytosolic. A subset of these are associated with PV/PECs. Extended and validated interactomes point to a core PV/PECs-associated UPS machinery, which includes uncharacterized and Giardia-specific coiled-coil proteins and NEK kinases. Finally, a subset of the alpha-giardin protein family was enriched in all PV/PECs-associated protein interactomes, highlighting a previously unappreciated role for these proteins at PV/PECs and in UPS. Taken together, our results provide the first characterization of a virulence-associated UPS protein complex in Giardia lamblia at PV/PECs, suggesting a novel link between these primarily endocytic and feeding organelles and UPS at the parasite-host interface.

Spliceosomal introns in the diplomonad parasite Giardia duodenalis revisited

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain 'WB' genome and further elucidating Giardia's gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.

Iron restriction increases the expression of a cytotoxic cysteine proteinase TvCP2 by a novel mechanism of tvcp2 mRNA alternative polyadenylation in Trichomonas vaginalis

Trichomonas vaginalis TvCP2 (TVAG_057000) is a cytotoxic cysteine proteinase (CP) expressed under iron-limited conditions. This work aimed to identify one of the mechanisms of tvcp2 gene expression regulation by iron at the posttranscriptional level. We checked tvcp2 mRNA stability under both iron-restricted (IR) and high iron (HI) conditions in the presence of actinomycin D. Greater stability of the tvcp2 mRNA under the IR than in HI conditions was observed, as expected. In silico analysis of the 3' regulatory region showed the presence of two putative polyadenylation signals in the tvcp2 transcript. By 3'-RACE assays, we demonstrated the existence of two isoforms of the tvcp2 mRNA with different 3'-UTR that resulted in more TvCP2 protein under IR than in HI conditions detected by WB assays. Additionally, we searched for homologs of the trichomonad polyadenylation machinery by an in silico analysis in the genome database, TrichDB. 16 genes that encode proteins that could be part of the trichomonad polyadenylation machinery were found. qRT-PCR assays showed that most of these genes were positively regulated by iron. Thus, our results show the presence of alternative polyadenylation as a novel iron posttranscriptional regulatory mechanism in T. vaginalis for the tvcp2 gene expression.

Stem-Loop Structures in Iron-Regulated mRNAs of Giardia duodenalis

Giardia duodenalis is a significant cause of waterborne and foodborne infections, day-care center outbreaks, and traveler's diarrhea worldwide. In protozoa such as Trichomonas vaginalis and Entamoeba histolytica, iron affects the growth, pathogenicity mechanisms, and expression of virulence genes. One of the proposed iron regulatory mechanisms is at the post-transcriptional level through an IRE/IRP-like (iron responsive element/iron regulatory protein) system. Recently, the expression of many putative giardial virulence factors in the free-iron levels has been reported in subsequent RNAseq experiments; however, the iron regulatory mechanism remains unknown. Thus, this work aimed to determine the effects of iron on the growth, gene expression, and presence of IRE-like structures in G. duodenalis. First, the parasite's growth kinetics at different iron concentrations were studied, and the cell viability was determined. It was observed that the parasite can adapt to an iron range from 7.7 to 500 µM; however, in conditions without iron, it is unable to survive in the culture medium. Additionally, the iron modulation of three genes was determined by RT-PCR assays. The results suggested that Actin, glucosamine-6-phosphate deaminase, and cytochrome b5 mRNA were down-regulated by iron. To investigate the presence of IRE-like structures, in silico analyses were performed for different mRNAs from the Giardia genome database. The Zuker mfold v2.4 web server and theoretical analysis were used to predict the secondary structures of the 91 mRNAs analyzed. Interestingly, the iron-induced downregulation of the genes analyzed corresponds to the location of the stem-loop structures found in their UTR regions. In conclusion, iron modulates the growth and expression of specific genes, likely due to the presence of IRE-like structures in G. duodenalis mRNAs.

Comprehensive characterization of Cysteine-rich protein-coding genes of Giardia lamblia and their role during antigenic variation

Giardia lamblia encodes several families of cysteine-rich proteins, including the Variant-specific Surface Proteins (VSPs) involved in the process of antigenic variation. Their characteristics, definition and relationships are still controversial. An exhaustive analysis of the Cys-rich families including organization, features, evolution and levels of expression was performed, by combining pattern searches and predictions with massive sequencing techniques. Thus a new classification for Cys-rich proteins, genes and pseudogenes that better describes their involvement in Giardia's biology is presented. Moreover, three novel characteristics exclusive to the VSP genes, comprising an Initiator element/Kozak-like sequence, an extended polyadenylation signal and a unique pattern of mutually exclusive transcript accumulation is presented as well as the finding that High Cysteine Membrane Proteins, upregulated under stress, may protect the parasite during VSP switching. These results allow better interpretation of previous reports providing the basis for further studies of the biology of this early-branching eukaryote.

Implementation of a tunable t-CRISPRi system for gene regulation in Giardia duodenalis

Giardia duodenalis, is a binuclear and microaerophilic protozoan that causes giardiasis. Up to date, several molecular approaches have been taken to understand the molecular mechanisms of diverse cellular processes in this parasitic protozoan. However, the role of many genes involved in these processes needs further analysis. The CRISPR interference (CRISPRi) system has been widely used, as a constitutive expression system for gene silencing purposes in several parasites, including Giardia. The aim of this work was to implement a tunable t-CRISPRi system in Giardia to silence abundant, moderately and low expressed genes, by constructing an optimized and inducible plasmid for the expression of both gRNA and dCas9. A doxycycline inducible pRan promoter was used to express dCas9 and each gRNA, consistently dCas9 expression and nuclear localization were confirmed by Western-blot and immunofluorescence in transfected trophozoites. The transcriptional repression was performed on α-tubulin (high expression), giardipain-1 (moderate expression) and Sir2 and Sir4 (low expression) genes. The α-tubulin gene knock-down caused by dCas9 doxycycline-induction was confirmed by a decrease in its protein expression which was of 50% and 60% at 24 and 48 h, respectively. This induced morphological alterations in flagella. The giardipain-1 knock down, showed a decrease in protein expression of 40 and 50% at 12 and 24 h, respectively, without affecting trophozoites viability, consistent with this a zymogram analysis on giardipain-1 knock down revealed a decrease in giardipain-1 protease activity. When repressing sirtuins expression, a total repression was obtained but trophozoites viability was compromised. This approach provides a molecular tool for a tailored repression to produce specific gene knockdowns.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Dual RNA Sequencing Reveals Key Events When Different Giardia Life Cycle Stages Interact With Human Intestinal Epithelial Cells In Vitro

Giardia intestinalis is a protozoan parasite causing diarrheal disease, giardiasis, after extracellular infection of humans and other mammals’ intestinal epithelial cells (IECs) of the upper small intestine. The parasite has two main life cycle stages: replicative trophozoites and transmissive cysts. Differentiating parasites (encysting cells) and trophozoites have recently been shown to be present in the same regions of the upper small intestine, whereas most mature cysts are found further down in the intestinal system. To learn more about host-parasite interactions during Giardia infections, we used an in vitro model of the parasite’s interaction with host IECs (differentiated Caco-2 cells) and Giardia WB trophozoites, early encysting cells (7 h), and cysts. Dual RNA sequencing (Dual RNAseq) was used to identify differentially expressed genes (DEGs) in both Giardia and the IECs, which might relate to establishing infection and disease induction. In the human cells, the largest gene expression changes were found in immune and MAPK signaling, transcriptional regulation, apoptosis, cholesterol metabolism and oxidative stress. The different life cycle stages of Giardia induced a core of similar DEGs but at different levels and there are many life cycle stage-specific DEGs. The metabolic protein PCK1, the transcription factors HES7, HEY1 and JUN, the peptide hormone CCK and the mucins MUC2 and MUC5A are up-regulated in the IECs by trophozoites but not cysts. Cysts specifically induce the chemokines CCL4L2, CCL5 and CXCL5, the signaling protein TRKA and the anti-bacterial protein WFDC12. The parasite, in turn, up-regulated a large number of hypothetical genes, high cysteine membrane proteins (HCMPs) and oxidative stress response genes. Early encysting cells have unique DEGs compared to trophozoites (e.g. several uniquely up-regulated HCMPs) and interaction of these cells with IECs affected the encystation process. Our data show that different life cycle stages of Giardia induce different gene expression responses in the host cells and that the IECs in turn differentially affect the gene expression in trophozoites and early encysting cells. This life cycle stage-specific host-parasite cross-talk is an important aspect to consider during further studies of Giardia’s molecular pathogenesis.

Sequences to Differences in Gene Expression: Analysis of RNA-Seq Data

RNA-Seq is now a routinely employed assay to measure gene expression. As the technique matured over the last decade, so have dedicated analytic tools. In this chapter, we first describe the mainstream as well as the most up-to-date protocols and their implications on downstream analysis. We then detail the steps entailing RNA-Seq analysis in three main stages: (i) preprocessing and data preparation, (ii) upstream processing, and (iii) high-level analyses. We review the most recent and relevant tools as one workflow following a stepwise order. The chapter further encompasses in-depth features of these tools. Details of the required code are made available throughout the chapter, as well as of the underlying statistics. We illustrate these steps with analysis of publicly available RNA-Seq data.

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.

The role of nuclear organization in trans-splicing based expression of heat shock protein 90 in Giardia lamblia

Hsp90 gene of G. lamblia has a split nature comprising two ORFs separated by 777 kb on chromosome 5. The ORFs of the split gene on chromosome 5 undergo transcription to generate independent pre-mRNAs that join by a unique trans-splicing reaction that remains partially understood. The canonical cis-acting nucleotide elements such as 5'SS-GU, 3'SS-AG, polypyrimidine tract and branch point adenine are present in the independent pre-mRNAs and therefore trans-splicing of Hsp90 must be assisted by spliceosomes in vivo. Using an approach of RNA-protein pull down, we show that an RNA helicase selectively interacts with HspN pre-mRNA. Our experiments involving high resolution chromosome conformation capture technology as well as DNA FISH show that the trans-spliced genes of Giardia are in three-dimensional spatial proximity in the nucleus. Altogether our study provides a glimpse into the in vivo mechanisms involving protein factors as well as chromatin structure to facilitate the unique inter-molecular post-transcriptional stitching of split genes in G. lamblia.

Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB

Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite-induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duodenalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G. duodenalis CUB formation. This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB. These results will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

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.

Analysis of fungal genomes reveals commonalities of intron gain or loss and functions in intron-poor species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

The compact genome of Giardia muris reveals important steps in the evolution of intestinal protozoan parasites

Diplomonad parasites of the genus Giardia have adapted to colonizing different hosts, most notably the intestinal tract of mammals. The human-pathogenic Giardia species, Giardia intestinalis, has been extensively studied at the genome and gene expression level, but no such information is available for other Giardia species. Comparative data would be particularly valuable for Giardia muris, which colonizes mice and is commonly used as a prototypic in vivo model for investigating host responses to intestinal parasitic infection. Here we report the draft-genome of G. muris. We discovered a highly streamlined genome, amongst the most densely encoded ever described for a nuclear eukaryotic genome. G. muris and G. intestinalis share many known or predicted virulence factors, including cysteine proteases and a large repertoire of cysteine-rich surface proteins involved in antigenic variation. Different to G. intestinalis, G. muris maintains tandem arrays of pseudogenized surface antigens at the telomeres, whereas intact surface antigens are present centrally in the chromosomes. The two classes of surface antigens engage in genetic exchange. Reconstruction of metabolic pathways from the G. muris genome suggest significant metabolic differences to G. intestinalis. Additionally, G. muris encodes proteins that might be used to modulate the prokaryotic microbiota. The responsible genes have been introduced in the Giardia genus via lateral gene transfer from prokaryotic sources. Our findings point to important evolutionary steps in the Giardia genus as it adapted to different hosts and it provides a powerful foundation for mechanistic exploration of host-pathogen interaction in the G. muris-mouse pathosystem.

Identification and evolutionary analysis of the nucleolar proteome of Giardia lamblia

Background

The nucleoli, including their proteomes, of higher eukaryotes have been extensively studied, while few studies about the nucleoli of the lower eukaryotes – protists were reported. Giardia lamblia, a protist with the controversy of whether it is an extreme primitive eukaryote or just a highly evolved parasite, might be an interesting object for carrying out the nucleolar proteome study of protists and for further examining the controversy.

Results

Using bioinformatics methods, we reconstructed G. lamblia nucleolar proteome (GiNuP) and the common nucleolar proteome of the three representative higher eukaryotes (human, Arabidopsis, yeast) (HEBNuP). Comparisons of the two proteomes revealed that: 1) GiNuP is much smaller than HEBNuP, but 78.4% of its proteins have orthologs in the latter; 2) More than 68% of the GiNuP proteins are involved in the “Ribosome related” function, and the others participate in the other functions, and these two groups of proteins are much larger and much smaller than those in HEBNuP, respectively; 3) Both GiNuP and HEBNuP have their own specific proteins, but HEBNuP has a much higher proportion of such proteins to participate in more categories of nucleolar functions.

Conclusion

For the first time the nucleolar proteome of a protist - Giardia was reconstructed. The results of comparison of it with the common proteome of three representative higher eukaryotes -- HEBNuP indicated that the simplicity of GiNuP is most probably a reflection of primitiveness but not just parasitic reduction of Giardia, and simultaneously revealed some interesting evolutionary phenomena about the nucleolus and even the eukaryotic cell, compositionally and functionally.

A new gene inventory of the ubiquitin and ubiquitin-like conjugation pathways in Giardia intestinalis

BACKGROUND

Ubiquitin (Ub) and Ub-like proteins (Ub-L) are critical regulators of complex cellular processes such as the cell cycle, DNA repair, transcription, chromatin remodeling, signal translation, and protein degradation. Giardia intestinalis possesses an experimentally proven Ub-conjugation system; however, a limited number of enzymes involved in this process were identified using basic local alignment search tool (BLAST). This is due to the limitations of BLAST’s ability to identify homologous functional regions when similarity between the sequences dips to < 30%. In addition Ub-Ls and their conjugating enzymes have not been fully elucidated in Giardia.

OBJETIVE

To identify the enzymes involved in the Ub and Ub-Ls conjugation processes using intelligent systems based on the hidden Markov models (HMMs).

METHODS

We performed an HMM search of functional Pfam domains found in the key enzymes of these pathways in Giardia’s proteome. Each open reading frame identified was analysed by sequence homology, domain architecture, and transcription levels.

FINDINGS

We identified 118 genes, 106 of which corresponded to the ubiquitination process (Ub, E1, E2, E3, and DUB enzymes). The E3 ligase group was the largest group with 82 members; 71 of which harbored a characteristic RING domain. Four Ub-Ls were identified and the conjugation enzymes for NEDD8 and URM1 were described for first time. The 3D model for Ub-Ls displayed the β-grasp fold typical. Furthermore, our sequence analysis for the corresponding activating enzymes detected the essential motifs required for conjugation.

MAIN CONCLUSIONS

Our findings highlight the complexity of Giardia’s Ub-conjugation system, which is drastically different from that previously reported, and provides evidence for the presence of NEDDylation and URMylation enzymes in the genome and transcriptome of G. intestinalis.

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.

mRNA Polyadenylation Machineries in Intestinal Protozoan Parasites

In humans, mRNA polyadenylation involves the participation of about 20 factors in four main complexes that recognize specific RNA sequences. Notably, CFIm25, CPSF73, and PAP have essential roles for poly(A) site selection, mRNA cleavage, and adenosine residues polymerization. Besides the relevance of polyadenylation for gene expression, information is scarce in intestinal protozoan parasites that threaten human health. To better understand polyadenylation in Entamoeba histolytica, Giardia lamblia, and Cryptosporidium parvum, which represent leading causes of diarrhea worldwide, genomes were screened for orthologs of human factors. Results showed that Entamoeba histolytica and C. parvum have 16 and 12 proteins out of the 19 human proteins used as queries, respectively, while G. lamblia seems to have the smallest polyadenylation machinery with only six factors. Remarkably, CPSF30, CPSF73, CstF77, PABP2, and PAP, which were found in all parasites, could represent the core polyadenylation machinery. Multiple genes were detected for several proteins in Entamoeba, while gene redundancy is lower in Giardia and Cryptosporidium. Congruently with their relevance in the polyadenylation process, CPSF73 and PAP are present in all parasites, and CFIm25 is only missing in Giardia. They conserve the functional domains and predicted folding of human proteins, suggesting they may have the same roles in polyadenylation.

The Oxymonad Genome Displays Canonical Eukaryotic Complexity in the Absence of a Mitochondrion

The discovery that the protist Monocercomonoides exilis completely lacks mitochondria demonstrates that these organelles are not absolutely essential to eukaryotic cells. However, the degree to which the metabolism and cellular systems of this organism have adapted to the loss of mitochondria is unknown. Here, we report an extensive analysis of the M. exilis genome to address this question. Unexpectedly, we find that M. exilis genome structure and content is similar in complexity to other eukaryotes and less "reduced" than genomes of some other protists from the Metamonada group to which it belongs. Furthermore, the predicted cytoskeletal systems, the organization of endomembrane systems, and biosynthetic pathways also display canonical eukaryotic complexity. The only apparent preadaptation that permitted the loss of mitochondria was the acquisition of the SUF system for Fe-S cluster assembly and the loss of glycine cleavage system. Changes in other systems, including in amino acid metabolism and oxidative stress response, were coincident with the loss of mitochondria but are likely adaptations to the microaerophilic and endobiotic niche rather than the mitochondrial loss per se. Apart from the lack of mitochondria and peroxisomes, we show that M. exilis is a fully elaborated eukaryotic cell that is a promising model system in which eukaryotic cell biology can be investigated in the absence of mitochondria.

Giardia intestinalis cystatin is a potent inhibitor of papain, parasite cysteine proteases and, to a lesser extent, human cathepsin B

Cystatins are important regulators of papain-like cysteine proteases. In the protozoan parasite Giardia intestinalis, papain-like cysteine proteases play an essential role in the parasite's biology and pathogenicity. Here, we characterized a cysteine protease inhibitor of G. intestinalis that belongs to type-I-cystatins. The parasite cystatin is shown to be a strong inhibitor of papain (K_i  ≈ 0.3 nm) and three parasite cysteine proteases (CP14019, CP16160 and CP16779, K_i  ≈ 0.9-5.8 nm), but a weaker inhibitor of human cathepsin B (K_i  ≈ 79.9 nm). The protein localizes mainly in the cytoplasm. Together, these data suggest that cystatin of G. intestinalis plays a role in the regulation of cysteine protease activities in the parasite and, possibly, in the interaction with the host.

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.

Characterization of mRNA polyadenylation in the apicomplexa

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

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.

Host specificity in the Giardia duodenalis species complex

Giardia duodenalis is a unicellular flagellated parasite that infects the gastrointestinal tract of a wide range of mammalian species, including humans. Investigations of protein and DNA polymorphisms revealed that G. duodenalis should be considered as a species complex, whose members, despite being morphologically indistinguishable, can be classified into eight groups, or Assemblages, separated by large genetic distances. Assemblages display various degree of host specificity, with Assemblages A and B occurring in humans and many other hosts, Assemblage C and D in canids, Assemblage E in hoofed animals, Assemblage F in cats, Assemblage G in rodents, and Assemblage H in pinnipeds. The factors determining host specificity are only partially understood, and clearly involve both the host and the parasite. Here, we review the results of in vitro and in vivo experiments, and clinical observations to highlight relevant biological and genetic differences between Assemblages, with a focus on human infection.

Giardia duodenalis Infections in Humans and Other Animals in China

Giardia duodenalis is an important zoonotic pathogen in both public and veterinary health, and has been genotyped into at least eight assemblages (A-H), each with a distinct host range. In recent years, this intestinal protozoan parasite has been identified widely in humans and various other animals, and has even been recorded in environmental contaminants. Along with whole genome sequencing of G. duodenalis, multilocus sequence typing is increasingly being used to characterize G. duodenalis isolates. Here, we review the epidemiology, genotyping, and subtyping of G. duodenalis from humans and a wide range of other animals, as well as from wastewater, in China.

Transcriptomic analysis of diplomonad parasites reveals a trans-spliced intron in a helicase gene in Giardia

Background

The mechanisms by which DNA sequences are expressed is the central preoccupation of molecular genetics. Recently, ourselves and others reported that in the diplomonad protist Giardia lamblia, the coding regions of several mRNAs are produced by ligation of independent RNA species expressed from distinct genomic loci. Such trans-splicing of introns was found to affect nearly as many genes in this organism as does classical cis-splicing of introns. These findings raised questions about the incidence of intron trans-splicing both across the G. lambliatranscriptome and across diplomonad diversity in general, however a dearth of transcriptomic data at the time prohibited systematic study of these questions.

Methods

I leverage newly available transcriptomic data from G. lamblia and the related diplomonad Spironucleus salmonicidato search for trans-spliced introns. My computational pipeline recovers all four previously reported trans-spliced introns in G. lamblia, suggesting good sensitivity.

Results

Scrutiny of thousands of potential cases revealed only a single additional trans-spliced intron in G. lamblia, in the p68 helicase gene, and no cases in S. salmonicida. The p68 intron differs from the previously reported trans-spliced introns in its high degree of streamlining: the core features of G. lamblia trans-spliced introns are closely packed together, revealing striking economy in the implementation of a seemingly inherently uneconomical molecular mechanism.

Discussion

These results serve to circumscribe the role of trans-splicing in diplomonads both in terms of the number of genes effected and taxonomically. Future work should focus on the molecular mechanisms, evolutionary origins and phenotypic implications of this intriguing phenomenon.

Cyst-Wall-Protein-1 is fundamental for Golgi-like organelle neogenesis and cyst-wall biosynthesis in Giardia lamblia

The genome of the protozoan parasite Giardia lamblia is organized in two diploid nuclei, which has so far precluded complete analysis of gene function. Here we use a previously developed Cre/loxP-based knock-out and selection marker salvage strategy in the human-derived isolate WB-C6 to eliminate all four copies of the Cyst-Wall-Protein-1 locus (CWP1). Because these loci are silenced in proliferating trophozoites and highly expressed only in encysting cells, CWP1 ablation allows functional characterization of a conditional phenotype in parasites induced to encyst. We show that encysting Δcwp1 cells are unable to establish the stage-regulated trafficking machinery with Golgi-like encystation-specific vesicles required for cyst-wall formation but show morphological hallmarks of cyst development and karyokinesis. This ‘pseudocyst' phenotype is rescued by transfection of Δcwp1 cells with an episomally maintained CWP1 expression vector. Genome editing in genera Giardia and Trypanosoma are the only reported examples addressing questions on pathogen transmission within the Excavata supergroup.

Giardia lamblia is a human protozoan parasite with two diploid nuclei, which makes complete knock-out of a gene of interest challenging. Here the authors use a Cre/loxP-based approach to knock-out cyst-wall protein 1 (cwp1) and show that CWP1 is essential for cyst-wall biosynthesis.

Inhibición parcial de dos genes que codifican para proteínas spliceosomales en Giardia intestinalis

<p><strong>Introducción.</strong> <em>Giardia intestinalis</em> es un organismo tempranamente divergente en el que recientemente se demostró la presencia de intrones. La maquinaria responsable de la remoción de intrones en eucariotes superiores es el spliceosoma, conformado por 5 ribonucleoproteínas. Cada una tiene un ARN pequeño nuclear, un set de 7 proteínas Sm (B, D1, D2, D3, E, F y G) y varias proteínas específicas. En <em>G. intestinalis</em>, se han identificado los genes de algunas proteínas del spliceosoma por bioinformática. Aunque se asume que el spliceosoma es el responsable del splicing en el parásito, su caracterización bioquímica no ha sido realizada.</p><p><strong>Objetivo.</strong> Inhibir dos genes que codifican para proteínas del spliceosoma de <em>G. intestinalis</em> con el fin de determinar si esta inhibición afecta el crecimiento o la enquistación del parásito.</p><p><strong>Materiales y métodos.</strong> Se clonaron en un vector específico para <em>G. intestinalis</em> secuencias antisentido de los genes que codifican para las proteínas spliceosomales SmB y SmD3 del parásito. Posteriormente se transfectó <em>G. intestinalis</em> con los vectores recombinantes y se seleccionaron aquellos parásitos que lo incorporaron. Se confirmó la disminución del mensajero por PCR en tiempo real y se evaluó el crecimiento y la enquistación en parásitos silvestres y transfectados.</p><p><strong>Resultados.</strong> Se observó una disminución del 40% y 70% en el mARN de SmB y SmD3, respectivamente. El crecimiento y la enquistación no se vieron afectados en estos parásitos.</p><p><strong>Conclusión. </strong>La disminución de SmB y SmD3 no afectan el parásito, indicando que el spliceosoma sigue siendo funcional o que el splicing no es una función vital del parásito.</p>

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.

Time-Dependent Transcriptional Changes in Axenic Giardia duodenalis Trophozoites

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. Although G. duodenalis can be cultured axenically, significant gaps exist in our understanding of the molecular biology and metabolism of this pathogen. The present study employed RNA sequencing to characterize the mRNA transcriptome of G. duodenalis trophozoites in axenic culture, at log (48 h of growth), stationary (60 h), and declining (96 h) growth phases. Using ~400-times coverage of the transcriptome, we identified 754 differentially transcribed genes (DTGs), mainly representing two large DTG groups: 438 that were down-regulated in the declining phase relative to log and stationary phases, and 281 that were up-regulated. Differential transcription of prominent antioxidant and glycolytic enzymes implicated oxygen tension as a key factor influencing the transcriptional program of axenic trophozoites. Systematic bioinformatic characterization of numerous DTGs encoding hypothetical proteins of unknown function was achieved using structural homology searching. This powerful approach greatly informed the differential transcription analysis and revealed putative novel antioxidant-coding genes, and the presence of a near-complete two-component-like signaling system that may link cytosolic redox or metabolite sensing to the observed transcriptional changes. Motif searching applied to promoter regions of the two large DTG groups identified different putative transcription factor-binding motifs that may underpin global transcriptional regulation. This study provides new insights into the drivers and potential mediators of transcriptional variation in axenic G. duodenalis and provides context for static transcriptional studies.

Giardia is the most common gastrointestinal protozoan parasite of humans. This parasite causes diarrheal disease and is correlated with post-infectious conditions such as irritable bowel syndrome. In the absence of a vaccine, treatment is limited to drugs such as metronidazole, against which clinical resistance is reported. Effective control of Giardia requires a detailed understanding of its biology, and in turn, complete characterization of the standard in vitro culture system. Using RNA sequencing assisted by informatics to functionally annotate hypothetical proteins, we investigated transcriptional changes in axenic Giardia trophozoites at three growth phases over 96 hours. We found two large groups of differentially transcribed genes that indicate changes in the antioxidant system and central carbon metabolism over time. A putative novel signaling pathway may act together with putative transcription factor-binding motifs to regulate these transcriptional changes. Our results suggest that dissolved oxygen in Giardia culture medium may cause oxidative stress early during in vitro growth and that oxygen depletion may limit the efficiency of glycolysis in the declining phase. This work enhances our understanding of the transcriptional flexibility and metabolism of Giardia in vitro.

Comparative genomic analyses of freshly isolated Giardia intestinalis assemblage A isolates

Background

The diarrhea-causing protozoan Giardia intestinalis makes up a species complex of eight different assemblages (A-H), where assemblage A and B infect humans. Comparative whole-genome analyses of three of these assemblages have shown that there is significant divergence at the inter-assemblage level, however little is currently known regarding variation at the intra-assemblage level. We have performed whole genome sequencing of two sub-assemblage AII isolates, recently axenized from symptomatic human patients, to study the biological and genetic diversity within assemblage A isolates.

Results

Several biological differences between the new and earlier characterized assemblage A isolates were identified, including a difference in growth medium preference. The two AII isolates were of different sub-assemblage types (AII-1 [AS175] and AII-2 [AS98]) and showed size differences in the smallest chromosomes. The amount of genetic diversity was characterized in relation to the genome of the Giardia reference isolate WB, an assemblage AI isolate. Our analyses indicate that the divergence between AI and AII is approximately 1 %, represented by ~100,000 single nucleotide polymorphisms (SNP) distributed over the chromosomes with enrichment in variable genomic regions containing surface antigens. The level of allelic sequence heterozygosity (ASH) in the two AII isolates was found to be 0.25–0.35 %, which is 25–30 fold higher than in the WB isolate and 10 fold higher than the assemblage AII isolate DH (0.037 %). 35 protein-encoding genes, not found in the WB genome, were identified in the two AII genomes. The large gene families of variant-specific surface proteins (VSPs) and high cysteine membrane proteins (HCMPs) showed isolate-specific divergences of the gene repertoires. Certain genes, often in small gene families with 2 to 8 members, localize to the variable regions of the genomes and show high sequence diversity between the assemblage A isolates. One of the families, Bactericidal/Permeability Increasing-like protein (BPIL), with eight members was characterized further and the proteins were shown to localize to the ER in trophozoites.

Conclusions

Giardia genomes are modular with highly conserved core regions mixed up by variable regions containing high levels of ASH, SNPs and variable surface antigens. There are significant genomic variations in assemblage A isolates, in terms of chromosome size, gene content, surface protein repertoire and gene polymorphisms and these differences mainly localize to the variable regions of the genomes. The large genetic differences within one assemblage of G. intestinalis strengthen the argument that the assemblages represent different Giardia species.

Electronic supplementary material

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

Giardia fatty acyl-CoA synthetases as potential drug targets

Giardiasis caused by Giardia intestinalis (syn. G. lamblia, G. duodenalis) is one of the leading causes of diarrheal parasitic diseases worldwide. Although limited drugs to treat giardiasis are available, there are concerns regarding toxicity in some patients and the emerging drug resistance. By data-mining genome sequences, we observed that G. intestinalis is incapable of synthesizing fatty acids (FA) de novo. However, this parasite has five long-chain fatty acyl-CoA synthetases (GiACS1 to GiACS5) to activate FA scavenged from the host. ACS is an essential enzyme because FA need to be activated to form acyl-CoA thioesters before they can enter subsequent metabolism. In the present study, we performed experiments to explore whether some GiACS enzymes could serve as drug targets in Giardia. Based on the high-throughput datasets and protein modeling analyses, we initially studied the GiACS1 and GiACS2, because genes encoding these two enzymes were found to be more consistently expressed in varied parasite life cycle stages and when interacting with host cells based on previously reported transcriptome data. These two proteins were cloned and expressed as recombinant proteins. Biochemical analysis revealed that both had apparent substrate preference toward palmitic acid (C16:0) and myristic acid (C14:0), and allosteric or Michaelis–Menten kinetics on palmitic acid or ATP. The ACS inhibitor triacsin C inhibited the activity of both enzymes (IC₅₀ = 1.56 μM, K_i = 0.18 μM for GiACS1, and IC₅₀ = 2.28 μM, K_i = 0.23 μM for GiACS2, respectively) and the growth of G. intestinalis in vitro (IC₅₀ = 0.8 μM). As expected from giardial evolutionary characteristics, both GiACSs displayed differences in overall folding structure as compared with their human counterparts. These observations support the notion that some of the GiACS enzymes may be explored as drug targets in this parasite.

Drug resistance in Giardia duodenalis

Giardia duodenalis is a microaerophilic parasite of the human gastrointestinal tract and a major contributor to diarrheal and post-infectious chronic gastrointestinal disease world-wide. Treatment of G. duodenalis infection currently relies on a small number of drug classes. Nitroheterocyclics, in particular metronidazole, have represented the front line treatment for the last 40 years. Nitroheterocyclic-resistant G. duodenalis have been isolated from patients and created in vitro, prompting considerable research into the biomolecular mechanisms of resistance. These compounds are redox-active and are believed to damage proteins and DNA after being activated by oxidoreductase enzymes in metabolically active cells. In this review, we explore the molecular phenotypes of nitroheterocyclic-resistant G. duodenalis described to date in the context of the protist's unusual glycolytic and antioxidant systems. We propose that resistance mechanisms are likely to extend well beyond currently described resistance-associated enzymes (i.e., pyruvate ferredoxin oxidoreductases and nitroreductases), to include NAD(P)H- and flavin-generating pathways, and possibly redox-sensitive epigenetic regulation. Mechanisms that allow G. duodenalis to tolerate oxidative stress may lead to resistance against both oxygen and nitroheterocyclics, with implications for clinical control. The present review highlights the potential for systems biology tools and advanced bioinformatics to further investigate the multifaceted mechanisms of nitroheterocyclic resistance in this important pathogen.

Natural genetic variation impacts expression levels of coding, non-coding, and antisense transcripts in fission yeast

Our current understanding of how natural genetic variation affects gene expression beyond well-annotated coding genes is still limited. The use of deep sequencing technologies for the study of expression quantitative trait loci (eQTLs) has the potential to close this gap. Here, we generated the first recombinant strain library for fission yeast and conducted an RNA-seq-based QTL study of the coding, non-coding, and antisense transcriptomes. We show that the frequency of distal effects (trans-eQTLs) greatly exceeds the number of local effects (cis-eQTLs) and that non-coding RNAs are as likely to be affected by eQTLs as protein-coding RNAs. We identified a genetic variation of swc5 that modifies the levels of 871 RNAs, with effects on both sense and antisense transcription, and show that this effect most likely goes through a compromised deposition of the histone variant H2A.Z. The strains, methods, and datasets generated here provide a rich resource for future studies.

Both endo-siRNAs and tRNA-derived small RNAs are involved in the differentiation of primitive eukaryote Giardia lamblia

Small RNAs (sRNAs), including microRNAs and endogenous siRNAs (endo-siRNAs), regulate most important biologic processes in eukaryotes, such as cell division and differentiation. Although sRNAs have been extensively studied in various eukaryotes, the role of sRNAs in the early emergence of eukaryotes is unclear. To address these questions, we deep sequenced the sRNA transcriptome of four different stages in the differentiation of Giardia lamblia, one of the most primitive eukaryotes. We identified a large number of endo-siRNAs in this fascinating parasitic protozoan and found that they were produced from live telomeric retrotransposons and three genomic regions (i.e., endo-siRNA generating regions [eSGRs]). eSGR-derived endo-siRNAs were proven to target mRNAs in trans. Gradual up-regulation of endo-siRNAs in the differentiation of Giardia suggested that they might be involved in the regulation of this process. This hypothesis was supported by the impairment of the differentiation ability of Giardia when GLDICER, essential for the biogenesis of endo-siRNAs, was knocked down. Endo-siRNAs are not the only sRNA regulators in Giardia differentiation, because a great number of tRNAs-derived sRNAs showed more dramatic expression changes than endo-siRNAs in this process. We totally identified five novel kinds of tRNAs-derived sRNAs and found that the biogenesis in four of them might be correlated with that of stress-induced tRNA-derived RNA (sitRNA), which was discovered in our previous studies. Our studies reveal an unexpected complex panorama of sRNA in G. lamblia and shed light on the origin and functional evolution of eukaryotic sRNAs.

Orthology detection combining clustering and synteny for very large datasets

The elucidation of orthology relationships is an important step both in gene function prediction as well as towards understanding patterns of sequence evolution. Orthology assignments are usually derived directly from sequence similarities for large data because more exact approaches exhibit too high computational costs. Here we present PoFF, an extension for the standalone tool Proteinortho, which enhances orthology detection by combining clustering, sequence similarity, and synteny. In the course of this work, FFAdj-MCS, a heuristic that assesses pairwise gene order using adjacencies (a similarity measure related to the breakpoint distance) was adapted to support multiple linear chromosomes and extended to detect duplicated regions. PoFF largely reduces the number of false positives and enables more fine-grained predictions than purely similarity-based approaches. The extension maintains the low memory requirements and the efficient concurrency options of its basis Proteinortho, making the software applicable to very large datasets.

Secondary loss of a cis-spliced intron during the divergence of Giardia intestinalis assemblages

BACKGROUND

Giardia intestinalis is a parasitic unicellular eukaryote with a highly reduced genome, in which only six cis-spliced and four trans-spliced introns have been discovered. However, we anticipate that more cis- and trans-spliced introns likely remain unidentified in genes encoding hypothetical proteins that occupy ca. 2/3 of all of the open reading frames (ORFs) in the Giardia genome. Consequently, comprehensive surveys of introns in ORFs for hypothetical proteins are critical for better understanding of the intron evolution in this organism.

RESULTS

In this study, we identified two novel cis-spliced introns in the draft genome data of G. intestinalis strain WB, by surveying the conserved sequence motifs shared amongst the previously known introns. G. intestinalis strains can be divided into phylogenetically distinct assemblages A-H, and all the introns identified in past studies are shared among the published genome data from strains WB, DH, GS, and P15 representing assemblages A1, A2, B, and E, respectively. Nevertheless one of the two novel introns identified in this study was found to be absent in strain P15.

CONCLUSION

By considering the organismal relationship among G. intestinalis assemblages A1, A2, B, and E, one of the two introns identified in this study has highly likely been lost after the divergence of the assemblages. On the basis of a sequence comparison between the intron-bearing loci in WB, DH, and GS genomes and the homologous but intron-free locus in P15 genome, we propose that the loss of this particular intron was mediated by integration of the DNA fragment reverse-transcribed from mature mRNAs.

The microRNAs in an ancient protist repress the variant-specific surface protein expression by targeting the entire coding sequence

microRNAs (miRNA) have been detected in the deeply branched protist, Giardia lamblia, and shown to repress expression of the family of variant-specific surface proteins (VSPs), only one of which is expressed in Giardia trophozoite at a given time. Three next-generation sequencing libraries of Giardia Argonaute-associated small RNAs were constructed and analyzed. Analysis of the libraries identified a total of 99 new putative miRNAs with a size primarily in the 26 nt range similar to the size previously predicted by the Giardia Dicer crystal structure and identified by our own studies. Bioinformatic analysis identified multiple putative miRNA target sites in the mRNAs of all 73 VSPs. The effect of miRNA target sites within a defined 3′-region were tested on two vsp mRNAs. All the miRNAs showed partial repression of the corresponding vsp expression and were additive when the targeting sites were separately located. But the combined repression still falls short of 100%. Two other relatively short vsp mRNAs with 15 and 11 putative miRNA target sites identified throughout their ORFs were tested with their corresponding miRNAs. The results indicate that; (1) near 100% repression of vsp mRNA expression can be achieved through the combined action of multiple miRNAs on target sites located throughout the ORF; (2) the miRNA machinery could be instrumental in repressing the expression of vsp genes in Giardia; (3) this is the first time that all the miRNA target sites in the entire ORF of a mRNA have been tested and shown to be functional.

Giardia lamblia is a protozoan parasite causing the diarrheal disease giardiasis. Variant-specific surface proteins (VSP) in Giardia are likely involved in its evasion of host immune response. Their expression is regulated by microRNAs (miRNA). To determine the full complement of miRNAs in Giardia, three cDNA libraries of Giardia Argonaute associated small RNAs were constructed and analyzed to identify a total of 105 miRNAs. Bioinformatic target identification showed that 102 of the 105 miRNAs find their putative target sites in vsp mRNAs. When only the target sites within the 3′ region,100 nts upstream of the stop codon, were tested against their corresponding miRNAs, however, only partial repression of VSP expression was observed. When all the miRNA target sites in the open reading frames of vsp mRNAs were examined, however, they all turned out to be functional. A saturation of them with the corresponding miRNAs resulted in a full repression of VSP expression, suggesting that this is the mechanism of miRNA repression of VSP expression in Giardia. The ability of miRNAs to regulate target sites throughout the entire open reading frame also provides the first indication that all the miRNA target sites in an mRNA are functional.

The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments

Spironucleus salmonicida causes systemic infections in salmonid fish. It belongs to the group diplomonads, binucleated heterotrophic flagellates adapted to micro-aerobic environments. Recently we identified energy-producing hydrogenosomes in S. salmonicida. Here we present a genome analysis of the fish parasite with a focus on the comparison to the more studied diplomonad Giardia intestinalis. We annotated 8067 protein coding genes in the ∼12.9 Mbp S. salmonicida genome. Unlike G. intestinalis, promoter-like motifs were found upstream of genes which are correlated with gene expression, suggesting a more elaborate transcriptional regulation. S. salmonicida can utilise more carbohydrates as energy sources, has an extended amino acid and sulfur metabolism, and more enzymes involved in scavenging of reactive oxygen species compared to G. intestinalis. Both genomes have large families of cysteine-rich membrane proteins. A cluster analysis indicated large divergence of these families in the two diplomonads. Nevertheless, one of S. salmonicida cysteine-rich proteins was localised to the plasma membrane similar to G. intestinalis variant-surface proteins. We identified S. salmonicida homologs to cyst wall proteins and showed that one of these is functional when expressed in Giardia. This suggests that the fish parasite is transmitted as a cyst between hosts. The extended metabolic repertoire and more extensive gene regulation compared to G. intestinalis suggest that the fish parasite is more adapted to cope with environmental fluctuations. Our genome analyses indicate that S. salmonicida is a well-adapted pathogen that can colonize different sites in the host.

Studies of model organisms are very powerful. However, to appreciate the enormous diversity of genetic and cell biological processes we need to extend the number of available model organisms. For example, there are very few model organisms for diverse microbial eukaryotes, a group of organisms which indeed represents the vast majority of the eukaryotic diversity. To this end, we have developed a system to do genetic modification on the Atlantic salmon pathogen Spironucleus salmonicida. Using this system we could show that the organism is capable of producing hydrogen within specialised compartments. Here we present the genome sequence of S. salmonicida together with a thorough annotation. We compare the results with the closest available model organism, the human intestinal parasite Giardia intestinalis. The fish parasite has a more elaborate system for regulation of gene expression, as well as a larger metabolic capacity. This indicates that S. salmonicida is a well-adapted pathogen that can deal with fluctuating environments, an important trait to be able to establish systemic infections in the host. The development of S. salmonicida into a model system will benefit the studies of fish infections, as well as cell biological processes.

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”.