Time-Dependent Transcriptional Changes in Axenic Giardia duodenalis Trophozoites

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.

Transcriptomic analysis of albendazole resistance in human diarrheal parasite Giardia duodenalis

Benzimidazole-2-carbamates (BZ, e.g., albendazole; ALB), which bind β-tubulin to disrupt microtubule polymerization, are one of two primary compound classes used to treat giardiasis. In most parasitic nematodes and fungi, BZ-resistance is caused by β-tubulin mutations and its molecular mode of action (MOA) is well studied. In contrast, in Giardia duodenalis BZ MOA or resistance is less well understood, may involve target-specific and broader impacts including cellular damage and oxidative stress, and its underlying cause is not clearly determined. Previously, we identified acquisition of a single nucleotide polymorphism, E198K, in β-tubulin in ALB-resistant (ALB-R) G. duodenalis WB-1B relative to ALB-sensitive (ALB-S) parental controls. E198K is linked to BZ-resistance in fungi and its allelic frequency correlated with the magnitude of BZ-resistance in G. duodenalis WB-1B. Here, we undertook detailed transcriptomic comparisons of these ALB-S and ALB-R G. duodenalis WB-1B cultures. The primary transcriptional changes with ALB-R in G. duodenalis WB-1B indicated increased protein degradation and turnover, and up-regulation of tubulin, and related genes, associated with the adhesive disc and basal bodies. These findings are consistent with previous observations noting focused disintegration of the disc and associated structures in Giardia duodenalis upon ALB exposure. We also saw transcriptional changes with ALB-R in G. duodenalis WB-1B consistent with prior observations of a shift from glycolysis to arginine metabolism for ATP production and possible changes to aspects of the vesicular trafficking system that require further investigation. Finally, we saw mixed transcriptional changes associated with DNA repair and oxidative stress responses in the G. duodenalis WB-1B line. These changes may be indicative of a role for H₂O₂ degradation in ALB-R, as has been observed in other G. duodenalis cell cultures. However, they were below the transcriptional fold-change threshold (log₂FC > 1) typically employed in transcriptomic analyses and appear to be contradicted in ALB-R G. duodenalis WB-1B by down-regulation of the NAD scavenging and conversion pathways required to support these stress pathways and up-regulation of many highly oxidation sensitive iron-sulphur (FeS) cluster based metabolic enzymes.

Gut Microbial Perturbation and Host Response Induce Redox Pathway Upregulation along the Gut-Liver Axis during Giardiasis in C57BL/6J Mouse Model

Apicomplexan infections, such as giardiasis and cryptosporidiosis, negatively impact a considerable proportion of human and commercial livestock populations. Despite this, the molecular mechanisms of disease, particularly the effect on the body beyond the gastrointestinal tract, are still poorly understood. To highlight host-parasite-microbiome biochemical interactions, we utilised integrated metabolomics-16S rRNA genomics and metabolomics-proteomics approaches in a C57BL/6J mouse model of giardiasis and compared these to Cryptosporidium and uropathogenic Escherichia coli (UPEC) infections. Comprehensive samples (faeces, blood, liver, and luminal contents from duodenum, jejunum, ileum, caecum and colon) were collected 10 days post infection and subjected to proteome and metabolome analysis by liquid and gas chromatography-mass spectrometry, respectively. Microbial populations in faeces and luminal washes were examined using 16S rRNA metagenomics. Proteome-metabolome analyses indicated that 12 and 16 key pathways were significantly altered in the gut and liver, respectively, during giardiasis with respect to other infections. Energy pathways including glycolysis and supporting pathways of glyoxylate and dicarboxylate metabolism, and the redox pathway of glutathione metabolism, were upregulated in small intestinal luminal contents and the liver during giardiasis. Metabolomics-16S rRNA genetics integration indicated that populations of three bacterial families-Autopobiaceae (Up), Desulfovibrionaceae (Up), and Akkermanasiaceae (Down)-were most significantly affected across the gut during giardiasis, causing upregulated glycolysis and short-chained fatty acid (SCFA) metabolism. In particular, the perturbed Akkermanasiaceae population seemed to cause oxidative stress responses along the gut-liver axis. Overall, the systems biology approach applied in this study highlighted that the effects of host-parasite-microbiome biochemical interactions extended beyond the gut ecosystem to the gut-liver axis. These findings form the first steps in a comprehensive comparison to ascertain the major molecular and biochemical contributors of host-parasite interactions and contribute towards the development of biomarker discovery and precision health solutions for apicomplexan infections.

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.

Risk factors for and management of metronidazole-refractory giardiasis in international travellers: A retrospective analysis

BACKGROUND

Giardia lamblia is a common cause of diarrhoea in returning travellers. Failure of the recommended first-line treatment, metronidazole, has frequently been observed. Recommendations for treatment of metronidazole-refractory giardiasis lack clarity and evidence.

METHODS

We conducted a retrospective data analysis of returned travellers with confirmed giardiasis at the Bernhard-Nocht-Clinic in Hamburg, Germany, between 2007 and 2016.

RESULTS

We identified 339 cases of giardiasis, mostly acquired in South Asia (n = 157). 308 patients received metronidazole as first-line treatment, leading to treatment failure in 93 cases. Statistical analysis suggested by far the highest risk of metronidazole treatment failure for travellers returning from South Asia (Odds Ratio 8.73). Second-line therapy consisted of various different therapy regimens. Combination therapy as second-line treatment seemed to be more effective than monotherapy. A repeat course of metronidazole proved to be futile.

CONCLUSION

This study reveals a strikingly low effectiveness of metronidazole, especially in patients returning from South Asia. Second-line treatment showed inconsistency of regimens and yielded unsatisfactory results. These findings require reconsideration of treatment strategies for giardiasis. Large prospective trials are urgently needed to assess new first-line treatment options and to help implement advice for effective, agreed second-line treatment strategies. Translational projects should be created to link the understanding of resistance mechanisms with epidemiological data and clinical outcome.

Eukaryote-conserved histone post-translational modification landscape in Giardia duodenalis revealed by mass spectrometry

Diarrheal disease caused by Giardia duodenalis is highly prevalent, causing over 200 million cases globally each year. The processes that drive parasite virulence, host immune evasion and transmission involve coordinated gene expression and have been linked to epigenetic regulation. Epigenetic regulatory systems are eukaryote-conserved, including in deep branching excavates such as Giardia, with several studies already implicating histone post-translational modifications in regulation of its pathogenesis and life cycle. However, further insights into Giardia chromatin dynamics have been hindered by a lack of site-specific knowledge of histone modifications. Using mass spectrometry, we have provided the first known molecular map of histone methylation, acetylation and phosphorylation modifications in Giardia core histones. We have identified over 50 previously unreported histone modifications including sites with established roles in epigenetic regulation, and co-occurring modifications indicative of post-translational modification crosstalk. These demonstrate conserved histone modifications in Giardia which are equivalent to many other eukaryotes, and suggest that similar epigenetic mechanisms are in place in this parasite. Further, we used sequence, domain and structural homology to annotate putative histone enzyme networks in Giardia, highlighting representative chromatin modifiers which appear sufficient for identified sites, particularly those from H3 and H4 variants. This study is to our knowledge the first and most comprehensive, complete and accurate view of Giardia histone post-translational modifications to date, and a substantial step towards understanding their associations in parasite development and virulence.

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.

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.

Differential protein expression and post-translational modifications in metronidazole-resistant Giardia duodenalis

Background

Metronidazole (Mtz) is the frontline drug treatment for multiple anaerobic pathogens, including the gastrointestinal protist, Giardia duodenalis. However, treatment failure is common and linked to in vivo drug resistance. In Giardia, in vitro drug-resistant lines allow controlled experimental interrogation of resistance mechanisms in isogenic cultures. However, resistance-associated changes are inconsistent between lines, phenotypic data are incomplete, and resistance is rarely genetically fixed, highlighted by reversion to sensitivity after drug selection ceases or via passage through the life cycle. Comprehensive quantitative approaches are required to resolve isolate variability, fully define Mtz resistance phenotypes, and explore the role of post-translational modifications therein.

Findings

We performed quantitative proteomics to describe differentially expressed proteins in 3 seminal Mtz-resistant lines compared to their isogenic, Mtz-susceptible, parental line. We also probed changes in post-translational modifications including protein acetylation, methylation, ubiquitination, and phosphorylation via immunoblotting. We quantified more than 1,000 proteins in each genotype, recording substantial genotypic variation in differentially expressed proteins between isotypes. Our data confirm substantial changes in the antioxidant network, glycolysis, and electron transport and indicate links between protein acetylation and Mtz resistance, including cross-resistance to deacetylase inhibitor trichostatin A in Mtz-resistant lines. Finally, we performed the first controlled, longitudinal study of Mtz resistance stability, monitoring lines after cessation of drug selection, revealing isolate-dependent phenotypic plasticity.

Conclusions

Our data demonstrate understanding that Mtz resistance must be broadened to post-transcriptional and post-translational responses and that Mtz resistance is polygenic, driven by isolate-dependent variation, and is correlated with changes in protein acetylation networks.

Transcriptomics Indicates Active and Passive Metronidazole Resistance Mechanisms in Three Seminal Giardia Lines

Giardia duodenalis is an intestinal parasite that causes 200-300 million episodes of diarrhoea annually. Metronidazole (Mtz) is a front-line anti-giardial, but treatment failure is common and clinical resistance has been demonstrated. Mtz is thought to be activated within the parasite by oxidoreductase enzymes, and to kill by causing oxidative damage. In G. duodenalis, Mtz resistance involves active and passive mechanisms. Relatively low activity of iron-sulfur binding proteins, namely pyruvate:ferredoxin oxidoreductase (PFOR), ferredoxins, and nitroreductase-1, enable resistant cells to passively avoid Mtz activation. Additionally, low expression of oxygen-detoxification enzymes can allow passive (non-enzymatic) Mtz detoxification via futile redox cycling. In contrast, active resistance mechanisms include complete enzymatic detoxification of the pro-drug by nitroreductase-2 and enhanced repair of oxidized biomolecules via thioredoxin-dependent antioxidant enzymes. Molecular resistance mechanisms may be largely founded on reversible transcriptional changes, as some resistant lines revert to drug sensitivity during drug-free culture in vitro, or passage through the life cycle. To comprehensively characterize these changes, we undertook strand-specific RNA sequencing of three laboratory-derived Mtz-resistant lines, 106-2ID10, 713-M3, and WB-M3, and compared transcription relative to their susceptible parents. Common up-regulated genes encoded variant-specific surface proteins (VSPs), a high cysteine membrane protein, calcium and zinc channels, a Mad-2 cell cycle regulator and a putative fatty acid α-oxidase. Down-regulated genes included nitroreductase-1, putative chromate and quinone reductases, and numerous genes that act proximal to PFOR. Transcriptional changes in 106-2ID10 diverged from those in 713-r and WB-r (r ≤ 0.2), which were more similar to each other (r = 0.47). In 106-2ID10, a nonsense mutation in nitroreductase-1 transcripts could enhance passive resistance whereas increased transcription of nitroreductase-2, and a MATE transmembrane pump system, suggest active drug detoxification and efflux, respectively. By contrast, transcriptional changes in 713-M3 and WB-M3 indicated a higher oxidative stress load, attributed to Mtz- and oxygen-derived radicals, respectively. Quantitative comparisons of orthologous gene transcription between Mtz-resistant G. duodenalis and Trichomonas vaginalis, a closely related parasite, revealed changes in transcripts encoding peroxidases, heat shock proteins, and FMN-binding oxidoreductases, as prominent correlates of resistance. This work provides deep insight into Mtz-resistant G. duodenalis, and illuminates resistance-associated features across parasitic species.

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.

Stem-Loop RT-qPCR as an Efficient Tool for the Detection and Quantification of Small RNAs in Giardia lamblia

Stem-loop quantitative reverse transcription PCR (RT-qPCR) is a molecular technique used for identification and quantification of individual small RNAs in cells. In this work, we used a Universal ProbeLibrary (UPL)-based design to detect-in a rapid, sensitive, specific, and reproducible way-the small nucleolar RNA (snoRNA) GlsR17 and its derived miRNA (miR2) of Giardia lamblia using a stem-loop RT-qPCR approach. Both small RNAs could be isolated from both total RNA and small RNA samples. Identification of the two small RNAs was carried out by sequencing the PCR-amplified small RNA products upon ligation into the pJET1.2/blunt vector. GlsR17 is constitutively expressed during the 72 h cultures of trophozoites, while the mature miR2 is present in 2-fold higher abundance during the first 48 h than at 72 h. Because it has been suggested that miRNAs in G. lamblia have an important role in the regulation of gene expression, the use of the stem-loop RT-qPCR method could be valuable for the study of miRNAs of G. lamblia. This methodology will be a powerful tool for studying gene regulation in G. lamblia, and will help to better understand the features and functions of these regulatory molecules and how they work within the RNA interference (RNAi) pathway in G. lamblia.

Divergent Transcriptional Responses to Physiological and Xenobiotic Stress in Giardia duodenalis

Understanding how parasites respond to stress can help to identify essential biological processes. Giardia duodenalis is a parasitic protist that infects the human gastrointestinal tract and causes 200 to 300 million cases of diarrhea annually. Metronidazole, a major antigiardial drug, is thought to cause oxidative damage within the infective trophozoite form. However, treatment efficacy is suboptimal, due partly to metronidazole-resistant infections. To elucidate conserved and stress-specific responses, we calibrated sublethal metronidazole, hydrogen peroxide, and thermal stresses to exert approximately equal pressure on trophozoite growth and compared transcriptional responses after 24 h of exposure. We identified 252 genes that were differentially transcribed in response to all three stressors, including glycolytic and DNA repair enzymes, a mitogen-activated protein (MAP) kinase, high-cysteine membrane proteins, flavin adenine dinucleotide (FAD) synthetase, and histone modification enzymes. Transcriptional responses appeared to diverge according to physiological or xenobiotic stress. Downregulation of the antioxidant system and α-giardins was observed only under metronidazole-induced stress, whereas upregulation of GARP-like transcription factors and their subordinate genes was observed in response to hydrogen peroxide and thermal stressors. Limited evidence was found in support of stress-specific response elements upstream of differentially transcribed genes; however, antisense derepression and differential regulation of RNA interference machinery suggest multiple epigenetic mechanisms of transcriptional control.