A chromosome-scale reference genome for Giardia intestinalis WB

Comprehensive analysis of flavohemoprotein copy number variation in Giardia intestinalis: exploring links to metronidazole resistance

BACKGROUND

Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified.

METHODS

This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages.

RESULTS

Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages.

CONCLUSIONS

The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite.

Extreme mitochondrial reduction in a novel group of free-living metamonads

Metamonads are a diverse group of heterotrophic microbial eukaryotes adapted to living in hypoxic environments. All metamonads but one harbour metabolically altered 'mitochondrion-related organelles' (MROs) with reduced functions, however the degree of reduction varies. Here, we generate high-quality draft genomes, transcriptomes, and predicted proteomes for five recently discovered free-living metamonads. Phylogenomic analyses placed these organisms in a group we name the 'BaSk' (Barthelonids+Skoliomonads) clade, a deeply branching sister group to the Fornicata, a phylum that includes parasitic and free-living flagellates. Bioinformatic analyses of gene models shows that these organisms are predicted to have extremely reduced MRO proteomes in comparison to other free-living metamonads. Loss of the mitochondrial iron-sulfur cluster assembly system in some organisms in this group appears to be linked to the acquisition in their common ancestral lineage of a SUF-like minimal system Fe/S cluster pathway by lateral gene transfer. One of the isolates, Skoliomonas litria, appears to have lost all other known MRO pathways. No proteins were confidently assigned to the predicted MRO proteome of this organism suggesting that the organelle has been lost. The extreme mitochondrial reduction observed within this free-living anaerobic protistan clade demonstrates that mitochondrial functions may be completely lost even in free-living organisms.

Investigating Antiprotozoal Chemotherapies with Novel Proteomic Tools-Chances and Limitations: A Critical Review

Identification of drug targets and biochemical investigations on mechanisms of action are major issues in modern drug development. The present article is a critical review of the classical "one drug"-"one target" paradigm. In fact, novel methods for target deconvolution and for investigation of resistant strains based on protein mass spectrometry have shown that multiple gene products and adaptation mechanisms are involved in the responses of pathogens to xenobiotics rather than one single gene or gene product. Resistance to drugs may be linked to differential expression of other proteins than those interacting with the drug in protein binding studies and result in complex cell physiological adaptation. Consequently, the unraveling of mechanisms of action needs approaches beyond proteomics. This review is focused on protozoan pathogens. The conclusions can, however, be extended to chemotherapies against other pathogens or cancer.

Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

Harnessing the power of new genetic tools to illuminate Giardia biology and pathogenesis

Giardia is a prevalent single-celled microaerophilic intestinal parasite causing diarrheal disease and significantly impacting global health. Double diploid (essentially tetraploid) Giardia trophozoites have presented a formidable challenge to the development of molecular genetic tools to interrogate gene function. High sequence divergence and the high percentage of hypothetical proteins lacking homology to proteins in other eukaryotes have limited our understanding of Giardia protein function, slowing drug target validation and development. For more than 25 years, Giardia A and B assemblages have been readily amenable to transfection with plasmids or linear DNA templates. Here, we highlight the utility and power of genetic approaches developed to assess protein function in Giardia, with particular emphasis on the more recent clustered regularly interspaced palindromic repeats/Cas9-based methods for knockdowns and knockouts. Robust and reliable molecular genetic approaches are fundamental toward the interrogation of Giardia protein function and evaluation of druggable targets. New genetic approaches tailored for the double diploid Giardia are imperative for understanding Giardia's unique biology and pathogenesis.

Exploration of Giardia small nucleolar RNAs (snoRNAs) and their possible microRNA derivatives

Small nucleolar RNAs (snoRNAs) are short non-coding RNAs that are abundant in the nucleoli of eukaryotic cells and play a crucial role in various aspects of ribosomal RNA (rRNA) maturation, including modifications such as 2′-O-methylation or pseudouridylation. On the other hand, Giardia duodenalis is a microaerophilic, flagellated, binucleate protozoan responsible for causing giardiasis. Although numerous snoRNAs have been detected in Giardia, their investigation remains limited. Nevertheless, they have been found to play a crucial role in the rRNA precursor processing pathway and influence other cellular functions. In addition, it has been proposed that some microRNAs are generated from these snoRNAs through excision by the Giardia endoribonuclease Dicer. These microRNAs are believed to contribute to the regulation of antigenic variation, which allows the parasite to evade the host immune response. Specifically, they play a role in modulating variant-specific surface proteins (VSPs) and other cysteine-rich surface antigens (CSAs). The main objective of this study was to bring together the available data on snoRNAs in Giardia, uncovering their functions in various processes and their importance on a global scale. In addition, the research delved into potential microRNAs speculated to originate from snoRNAs, exploring their impact on cellular processes.

Giardia telomeres and telomerase

Giardia duodenalis, the protozoan responsible for giardiasis, is a significant contributor to millions of diarrheal diseases worldwide. Despite the availability of treatments for this parasitic infection, therapeutic failures are alarmingly frequent. Thus, there is a clear need to identify new therapeutic targets. Giardia telomeres were previously identified, but our understanding of these structures and the critical role played by Giardia telomerase in maintaining genomic stability and its influence on cellular processes remains limited. In this regard, it is known that all Giardia chromosomes are capped by small telomeres, organized and protected by specific proteins that regulate their functions. To counteract natural telomere shortening and maintain high proliferation, Giardia exhibits constant telomerase activity and employs additional mechanisms, such as the formation of G-quadruplex structures and the involvement of transposable elements linked to telomeric repeats. Thus, this study aims to address the existing knowledge gap by compiling the available information (until 2023) about Giardia telomeres and telomerase, focusing on highlighting the distinctive features within this parasite. Furthermore, the potential feasibility of targeting Giardia telomeres and/or telomerase as an innovative therapeutic strategy is discussed.

The Giardia lamblia ribosome structure reveals divergence in several biological pathways and the mode of emetine function

Giardia lamblia is a deeply branching protist and a human pathogen. Its unusual biology presents the opportunity to explore conserved and fundamental molecular mechanisms. We determined the structure of the G. lamblia 80S ribosome bound to tRNA, mRNA, and the antibiotic emetine by cryo-electron microscopy, to an overall resolution of 2.49 Å. The structure reveals rapidly evolving protein and nucleotide regions, differences in the peptide exit tunnel, and likely altered ribosome quality control pathways. Examination of translation initiation factor binding sites suggests these interactions are conserved despite a divergent initiation mechanism. Highlighting the potential of G. lamblia to resolve conserved biological principles; our structure reveals the interactions of the translation inhibitor emetine with the ribosome and mRNA, thus providing insight into the mechanism of action for this widely used antibiotic. Our work defines key questions in G. lamblia and motivates future experiments to explore the diversity of eukaryotic gene regulation.

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.

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.1 50 years ago, "herpes-type" viral particles2 were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes,3,4 but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ∼300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

Genomic survey maps differences in the molecular complement of vesicle formation machinery between Giardia intestinalis assemblages

Giardia intestinalis is a globally important microbial pathogen with considerable public health, agricultural, and economic burden. Genome sequencing and comparative analyses have elucidated G. intestinalis to be a taxonomically diverse species consisting of at least eight different sub-types (assemblages A-H) that can infect a great variety of animal hosts, including humans. The best studied of these are assemblages A and B which have a broad host range and have zoonotic transmissibility towards humans where clinical Giardiasis can range from asymptomatic to diarrheal disease. Epidemiological surveys as well as previous molecular investigations have pointed towards critical genomic level differences within numerous molecular pathways and families of parasite virulence factors within assemblage A and B isolates. In this study, we explored the necessary machinery for the formation of vesicles and cargo transport in 89 Canadian isolates of assemblage A and B G. intestinalis. Considerable variability within the molecular complement of the endolysosomal ESCRT protein machinery, adaptor coat protein complexes, and ARF regulatory system have previously been reported. Here, we confirm inter-assemblage, but find no intra-assemblage variation within the trafficking systems examined. This variation includes losses of subunits belonging to the ESCRTIII as well as novel lineage specific duplications in components of the COPII machinery, ARF1, and ARFGEF families (BIG and CYTH). Since differences in disease manifestation between assemblages A and B have been controversially reported, our findings may well have clinical implications and even taxonomic, as the membrane trafficking system underpin parasite survival, pathogenesis, and propagation.

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.

Adaptation of the late ISC pathway in the anaerobic mitochondrial organelles of Giardia intestinalis

Mitochondrial metabolism is entirely dependent on the biosynthesis of the [4Fe-4S] clusters, which are part of the subunits of the respiratory chain. The mitochondrial late ISC pathway mediates the formation of these clusters from simpler [2Fe-2S] molecules and transfers them to client proteins. Here, we characterized the late ISC pathway in one of the simplest mitochondria, mitosomes, of the anaerobic protist Giardia intestinalis that lost the respiratory chain and other hallmarks of mitochondria. In addition to IscA2, Nfu1 and Grx5 we identified a novel BolA1 homologue in G. intestinalis mitosomes. It specifically interacts with Grx5 and according to the high-affinity pulldown also with other core mitosomal components. Using CRISPR/Cas9 we were able to establish full bolA1 knock out, the first cell line lacking a mitosomal protein. Despite the ISC pathway being the only metabolic role of the mitosome no significant changes in the mitosome biology could be observed as neither the number of the mitosomes or their capability to form [2Fe-2S] clusters in vitro was affected. We failed to identify natural client proteins that would require the [2Fe-2S] or [4Fe-4S] cluster within the mitosomes, with the exception of [2Fe-2S] ferredoxin, which is itself part of the ISC pathway. The overall uptake of iron into the cellular proteins remained unchanged as also observed for the grx5 knock out cell line. The pull-downs of all late ISC components were used to build the interactome of the pathway showing specific position of IscA2 due to its interaction with the outer mitosomal membrane proteins. Finally, the comparative analysis across Metamonada species suggested that the adaptation of the late ISC pathway identified in G. intestinalis occurred early in the evolution of this supergroup of eukaryotes.

The nucleolus of Giardia and its ribosomal biogenesis

Giardia duodenalis is a protozoan intestinal parasite that causes a significant number of infections worldwide each year, particularly in low-income and developing countries. Despite the availability of treatments for this parasitic infection, treatment failures are alarmingly common. As a result, new therapeutic strategies are urgently needed to effectively combat this disease. On the other hand, within the eukaryotic nucleus, the nucleolus stands out as the most prominent structure. It plays a crucial role in coordinating ribosome biogenesis and is involved in vital processes such as maintaining genome stability, regulating cell cycle progression, controlling cell senescence, and responding to stress. Given its significance, the nucleolus presents itself as a valuable target for selectively inducing cell death in undesirable cells, making it a potential avenue for anti-Giardia treatments. Despite its potential importance, the Giardia nucleolus remains poorly studied and often overlooked. In light of this, the objective of this study is to provide a detailed molecular description of the structure and function of the Giardia nucleolus, with a primary focus on its involvement in ribosomal biogenesis. Likewise, it discusses the targeting of the Giardia nucleolus as a therapeutic strategy, its feasibility, and the challenges involved.

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.

Genome-Wide Libraries for Protozoan Pathogen Drug Target Screening Using Yeast Surface Display

The lack of genetic tools to manipulate protozoan pathogens has limited the use of genome-wide approaches to identify drug or vaccine targets and understand these organisms' biology. We have developed an efficient method to construct genome-wide libraries for yeast surface display (YSD) and developed a YSD fitness screen (YSD-FS) to identify drug targets. We show the efficacy of our method by generating genome-wide libraries for Trypanosoma brucei, Trypanosoma cruzi, and Giardia lamblia parasites. Each library has a diversity of ∼105 to 106 clones, representing ∼6- to 30-fold of the parasite's genome. Nanopore sequencing confirmed the libraries' genome coverage with multiple clones for each parasite gene. Western blot and imaging analysis confirmed surface expression of the G. lamblia library proteins in yeast. Using the YSD-FS assay, we identified bonafide interactors of metronidazole, a drug used to treat protozoan and bacterial infections. We also found enrichment in nucleotide-binding domain sequences associated with yeast increased fitness to metronidazole, indicating that this drug might target multiple enzymes containing nucleotide-binding domains. The libraries are valuable biological resources for discovering drug or vaccine targets, ligand receptors, protein-protein interactions, and pathogen-host interactions. The library assembly approach can be applied to other organisms or expression systems, and the YSD-FS assay might help identify new drug targets in protozoan pathogens.

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.

Genomic comparisons confirm Giardia duodenalis sub-assemblage AII as a unique species

Giardia duodenalis is a parasitic flagellated protozoan which infects a wide range of mammalian hosts, including humans, and is subdivided into at least eight genetic assemblages commonly thought to represent cryptic species. Molecular studies have shown that G. duodenalis assemblage A, which parasitizes humans and animals, contains several phylogenetically distinct groupings known as sub-assemblages. Molecular studies employing poor phylogenetic-resolution markers routinely recover these sub-assemblages, implying that they represent evolutionarily distinct clades and possibly cryptic species, a hypothesis which is supported by epidemiologic trends. Here, we further tested this hypothesis by using available data from 41 whole genomes to characterize sub-assemblages and coalescent techniques for statistical estimation of species boundaries coupled to functional gene content analysis, thereby assessing the stability and distinctiveness of clades. Our analysis revealed two new sub-assemblage clades as well as novel signatures of gene content geared toward differential host adaptation and population structuring via vertical inheritance rather than recombination or panmixia. We formally propose sub-assemblage AII as a new species, Giardia hominis, while preserving the name Giardia duodenalis for sub-assemblage AI. Additionally, our bioinformatic methods broadly address the challenges of identifying cryptic microbial species to advance our understanding of emerging disease epidemiology, which should be broadly applicable to other lower eukaryotic taxa of interest. Giardia hominis n. sp. Zoobank LSID: urn:lsid: zoobank.org:pub:4298F3E1-E3EF-4977-B9DD-5CC59378C80E.

A chromosome-scale reference genome for Spironucleus salmonicida

Spironucleus salmonicida is a diplomonad causing systemic infection in salmon. The first S. salmonicida genome assembly was published 2014 and has been a valuable reference genome in protist research. However, the genome assembly is fragmented without assignment of the sequences to chromosomes. In our previous Giardia genome study, we have shown how a fragmented genome assembly can be improved with long-read sequencing technology complemented with optical maps. Combining Pacbio long-read sequencing technology and optical maps, we are presenting here this new S. salmonicida genome assembly in nine near-complete chromosomes with only three internal gaps at long repeats. This new genome assembly is not only more complete sequence-wise but also more complete at annotation level, providing more details into gene families, gene organizations and chromosomal structure. This near-complete reference genome will aid comparative genomics at chromosomal level, and serve as a valuable resource for the diplomonad community and protist research.

Measurement(s)	genomic_DNA • sequence_assembly • sequence feature annotation
Technology Type(s)	SMRT Sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Spironucleus salmonicida

Ribosomal RNA Transcription Machineries in Intestinal Protozoan Parasites: A Bioinformatic Analysis

PURPOSE

Ribosome biogenesis is a key process in all living organisms, energetically expensive and tightly regulated. Currently, little is known about the components of the ribosomal RNA (rRNA) transcription machinery that are present in intestinal parasites, such as Giardia duodenalis, Cryptosporidium parvum, and Entamoeba histolytica. Thus, in the present work, an analysis was carried out looking for the components of the rRNA transcription machinery that are conserved in intestinal parasites and if these could be used to design new treatment strategies.

METHODS

The different components of the rRNA transcription machinery were searched in the studied parasites with the NCBI BLAST tool in the EuPathDB Bioinformatics Resource Center database. The sequences of the RRN3 and POLR1F orthologs were aligned and important regions identified. Subsequently, three-dimensional models were built with different bioinformatic tools and a structural analysis was performed.

RESULTS

Among the protozoa examined, C. parvum is the parasite with the fewest identifiable components of the rRNA transcription machinery. TBP, RRN3, POLR1A, POLR1B, POLR1C, POLR1D, POLR1F, POLR1H, POLR2E, POLR2F and POLR2H subunits were identified in all species studied. Furthermore, the interaction regions between RRN3 and POLR1F were found to be conserved and could be used to design drugs that inhibit rRNA transcription in the parasites studied.

CONCLUSION

The inhibition of the rRNA transcription machinery in parasites might be a new therapeutic strategy against these microorganisms.

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.

Nucleosome Structures Built from Highly Divergent Histones: Parasites and Giant DNA Viruses

In eukaryotes, genomic DNA is bound with histone proteins and packaged into chromatin. The nucleosome, a fundamental unit of chromatin, regulates the accessibility of DNA to enzymes involved in gene regulation. During the past few years, structural analyses of chromatin architectures have been limited to evolutionarily related organisms. The amino acid sequences of histone proteins are highly conserved from humans to yeasts, but are divergent in the deeply branching protozoan groups, including human parasites that are directly related to human health. Certain large DNA viruses, as well as archaeal organisms, contain distant homologs of eukaryotic histone proteins. The divergent sequences give rise to unique and distinct nucleosome architectures, although the fundamental principles of histone folding and DNA contact are highly conserved. In this article, we review the structures and biophysical properties of nucleosomes containing histones from the human parasites Giardia lamblia and Leishmania major, and histone-like proteins from the Marseilleviridae amoeba virus family. The presented data confirm the sharing of the overall DNA compaction system among evolutionally distant species and clarify the deviations from the species-specific nature of the nucleosome.

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.

Efficient CRISPR/Cas9-mediated gene disruption in the tetraploid protist Giardia intestinalis

CRISPR/Cas9-mediated genome editing has become an extremely powerful technique used to modify gene expression in many organisms, including parasitic protists. Giardia intestinalis, a protist parasite that infects approximately 280 million people around the world each year, has been eluding the use of CRISPR/Cas9 to generate knockout cell lines due to its tetraploid genome. In this work, we show the ability of the in vitro assembled CRISPR/Cas9 components to successfully edit the genome of G. intestinalis. The cell line that stably expresses Cas9 in both nuclei of G. intestinalis showed effective recombination of the cassette containing the transcription units for the gRNA and the resistance marker. This highly efficient process led to the removal of all gene copies at once for three independent experimental genes, mem, cwp1 and mlf1. The method was also applicable to incomplete disruption of the essential gene, as evidenced by significantly reduced expression of tom40. Finally, testing the efficiency of Cas9-induced recombination revealed that homologous arms as short as 150 bp can be sufficient to establish a complete knockout cell line in G. intestinalis.

The controversies surrounding Giardia intestinalis assemblages A and B

Giardia intestinalis continues to be one of the most encountered parasitic diseases around the world. Although more frequently detected in developing countries, Giardia infections nonetheless pose significant public health problems in developed countries as well. Molecular characterisation of Giardia isolates from humans and animals reveals that there are two genetically different assemblages (known as assemblage A and B) that cause human infections. However, the current molecular assays used to genotype G. intestinalis isolates are quite controversial. This is in part due to a complex phenomenon where assemblages are incorrectly typed and underreported depending on which targeted locus is sequenced. In this review, we outline current knowledge based on molecular epidemiological studies and raise questions as to the reliability of current genotyping assays and a lack of a globally accepted method. Additionally, we discuss the clinical symptoms caused by G. intestinalis infection and how these symptoms vary depending on the assemblage infecting an individual. We also introduce the host-parasite factors that play a role in the subsequent clinical presentation of an infected person, and explore which assemblages are most seen globally.

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

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

A hybrid sequencing and assembly strategy for generating culture free Giardia genomes

Giardia duodenalis is a pathogenic intestinal protozoan parasite of humans and many other animals. Giardia duodenalis is found throughout the world, and infection is known to have adverse health consequences for human and other mammalian hosts. Yet, many aspects of the biology of this ubiquitous parasite remain unresolved. Whole genome sequencing and comparative genomics can provide insight into the biology of G. duodenalis by helping to reveal traits that are shared by all G. duodenalis assemblages or unique to an individual assemblage or strain. However, these types of analyses are currently hindered by the lack of available G. duodenalis genomes, due, in part, to the difficulty in obtaining the genetic material needed to perform whole genome sequencing. In this study, a novel approach using a multistep cleaning procedure coupled with a hybrid sequencing and assembly strategy was assessed for use in producing high quality G. duodenalis genomes directly from cysts obtained from feces of two naturally infected hosts, a cat and dog infected with assemblage A and D, respectively. Cysts were cleaned and concentrated using cesium chloride gradient centrifugation followed by immunomagnetic separation. Whole genome sequencing was performed using both Illumina MiSeq and Oxford Nanopore MinION platforms. A hybrid assembly strategy was found to produce higher quality genomes than assemblies from either platform alone. The hybrid G. duodenalis genomes obtained from fecal isolates (cysts) in this study compare favorably for quality and completeness against reference genomes of G. duodenalis from cultured isolates. The whole genome assembly for assemblage D is the most contiguous genome available for this assemblage and is an important reference genome for future comparative studies. The data presented here support a hybrid sequencing and assembly strategy as a suitable method to produce whole genome sequences from DNA obtained from G. duodenalis cysts which can be used to produce novel reference genomes necessary to perform comparative genomics studies of this parasite.

•

Assemblage A and D genomes were generated directly from cysts isolated from feces.

•

Genomes were sequenced using Illumina and Oxford Nanopore sequencing platforms.

•

A hybrid sequencing/assembly strategy was used to generate G. duodenalis genomes.

•

A hybrid strategy yields reference quality genomes from fecal isolates.

•

These methods have generated the most contiguous Assemblage D genome to date.

Not all Is SET for Methylation: Evolution of Eukaryotic Protein Methyltransferases

Dynamic posttranslational modifications to canonical histones that constitute the nucleosome (H2A, H2B, H3, and H4) control all aspects of enzymatic transactions with DNA. Histone methylation has been studied heavily for the past 20 years, and our mechanistic understanding of the control and function of individual methylation events on specific histone arginine and lysine residues has been greatly improved over the past decade, driven by excellent new tools and methods. Here, we will summarize what is known about the distribution and some of the functions of protein methyltransferases from all major eukaryotic supergroups. The main conclusion is that protein, and specifically histone, methylation is an ancient process. Many taxa in all supergroups have lost some subfamilies of both protein arginine methyltransferases (PRMT) and the heavily studied SET domain lysine methyltransferases (KMT). Over time, novel subfamilies, especially of SET domain proteins, arose. We use the interactions between H3K27 and H3K36 methylation as one example for the complex circuitry of histone modifications that make up the "histone code," and we discuss one recent example (Paramecium Ezl1) for how extant enzymes that may resemble more ancient SET domain KMTs are able to modify two lysine residues that have divergent functions in plants, fungi, and animals. Complexity of SET domain KMT function in the well-studied plant and animal lineages arose not only by gene duplication but also acquisition of novel DNA- and histone-binding domains in certain subfamilies.

High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis

Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.

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.

Unexpected organellar locations of ESCRT machinery in Giardia intestinalis and complex evolutionary dynamics spanning the transition to parasitism in the lineage Fornicata

BACKGROUND

Comparing a parasitic lineage to its free-living relatives is a powerful way to understand how that evolutionary transition to parasitism occurred. Giardia intestinalis (Fornicata) is a leading cause of gastrointestinal disease world-wide and is famous for its unusual complement of cellular compartments, such as having peripheral vacuoles instead of typical endosomal compartments. Endocytosis plays an important role in Giardia's pathogenesis. Endosomal sorting complexes required for transport (ESCRT) are membrane-deforming proteins associated with the late endosome/multivesicular body (MVB). MVBs are ill-defined in G. intestinalis, and roles for identified ESCRT-related proteins are not fully understood in the context of its unique endocytic system. Furthermore, components thought to be required for full ESCRT functionality have not yet been documented in this species.

RESULTS

We used genomic and transcriptomic data from several Fornicata species to clarify the evolutionary genome streamlining observed in Giardia, as well as to detect any divergent orthologs of the Fornicata ESCRT subunits. We observed differences in the ESCRT machinery complement between Giardia strains. Microscopy-based investigations of key components of ESCRT machinery such as GiVPS36 and GiVPS25 link them to peripheral vacuoles, highlighting these organelles as simplified MVB equivalents. Unexpectedly, we show ESCRT components associated with the endoplasmic reticulum and, for the first time, mitosomes. Finally, we identified the rare ESCRT component CHMP7 in several fornicate representatives, including Giardia and show that contrary to current understanding, CHMP7 evolved from a gene fusion of VPS25 and SNF7 domains, prior to the last eukaryotic common ancestor, over 1.5 billion years ago.

CONCLUSIONS

Our findings show that ESCRT machinery in G. intestinalis is far more varied and complete than previously thought, associates to multiple cellular locations, and presents changes in ESCRT complement which pre-date adoption of a parasitic lifestyle.

Giardia duodenalis: Biology and Pathogenesis

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

Real-Time PCR for Molecular Detection of Zoonotic and Non-Zoonotic Giardia spp. in Wild Rodents

Giardiasis in humans is a gastrointestinal disease transmitted by the potentially zoonotic Giardia duodenalis genotypes (assemblages) A and B. Small wild rodents such as mice and voles are discussed as potential reservoirs for G. duodenalis but are predominantly populated by the two rodent species Giardia microti and Giardia muris. Currently, the detection of zoonotic and non-zoonotic Giardia species and genotypes in these animals relies on cumbersome PCR and sequencing approaches of genetic marker genes. This hampers the risk assessment of potential zoonotic Giardia transmissions by these animals. Here, we provide a workflow based on newly developed real-time PCR schemes targeting the small ribosomal RNA multi-copy gene locus to distinguish G. muris, G. microti and G. duodenalis infections. For the identification of potentially zoonotic G. duodenalis assemblage types A and B, an established protocol targeting the single-copy gene 4E1-HP was used. The assays were specific for the distinct Giardia species or genotypes and revealed an analytical sensitivity of approximately one or below genome equivalent for the multi-copy gene and of about 10 genome equivalents for the single-copy gene. Retesting a biobank of small rodent samples confirmed the specificity. It further identified the underlying Giardia species in four out of 11 samples that could not be typed before by PCR and sequencing. The newly developed workflow has the potential to facilitate the detection of potentially zoonotic and non-zoonotic Giardia species in wild rodents.

Re-Discovery of Giardiavirus: Genomic and Functional Analysis of Viruses from Giardia duodenalis Isolates

Giardiasis, caused by the protozoan parasite Giardia duodenalis, is an intestinal diarrheal disease affecting almost one billion people worldwide. A small endosymbiotic dsRNA viruses, G. lamblia virus (GLV), genus Giardiavirus, family Totiviridae, might inhabit human and animal isolates of G. duodenalis. Three GLV genomes have been sequenced so far, and only one was intensively studied; moreover, a positive correlation between GLV and parasite virulence is yet to be proved. To understand the biological significance of GLV infection in Giardia, the characterization of several GLV strains from naturally infected G. duodenalis isolates is necessary. Here we report high-throughput sequencing of four GLVs strains, from Giardia isolates of human and animal origin. We also report on a new, unclassified viral sequence (designed GdRV-2), unrelated to Giardiavirus, encoding and expressing for a single large protein with an RdRp domain homologous to Totiviridae and Botybirnaviridae. The result of our sequencing and proteomic analyses challenge the current knowledge on GLV and strongly suggest that viral capsid protein translation unusually starts with a proline and that translation of the RNA-dependent RNA polymerase (RdRp) occurs via a +1/-2 ribosomal frameshift mechanism. Nucleotide polymorphism, confirmed by mass-spectrometry analysis, was also observed among and between GLV strains. Phylogenetic analysis indicated the occurrence of at least two GLV subtypes which display different phenotypes and transmissibility in experimental infections of a GLV naïve Giardia isolate.

Ribosomal DNA in the protozoan parasite Giardia duodenalis has a differential chromatin distribution and epigenetic markings across the subunits

Giardia duodenalis is a parasite that causes a large number of diarrheal diseases around the world. It is noteworthy that in a large number of processes, Giardia requires fewer components than other eukaryotes, even without some organelles such as mitochondria and peroxisomes. Despite this, core histones are known to exist in Giardia and epigenetic marks have been found on them, suggesting that they somehow control the expression of certain genes. The regulation of the expression of ribosomal DNA (rDNA) is essential, since it is required to maintain adequate levels of ribosomes and, given the nature of tandem repeat, it is a feasible area to create genomic instability. In Giardia, it is not known how this process occurs, but as in other eukaryotes, it is suggested through various epigenetic mechanisms. Thus, in the present work we seek to identify how chromatin is distributed through the Giardia rDNA and if there were histone marks that could control its expression.

The reduced ARF regulatory system in Giardia intestinalis pre-dates the transition to parasitism in the lineage Fornicata

Giardia intestinalis is an enteric pathogen with an extremely modified membrane trafficking system, lacking canonical compartments such as the Golgi, endosomes, and intermediate vesicle carriers. By comparison the fornicate relatives of Giardia possess greater endomembrane system complexity. In eukaryotes, the ADP ribosylation factor (ARF) GTPase regulatory system proteins, which consist of the small GTPase ARF1, and its guanine exchange nucleotide factors (GEFs) and GTPase activating proteins (GAPs), coordinate temporal and directional trafficking of cargo vesicles by recognizing and interacting with heterotetrameric coat complexes at pre-Golgi and post-Golgi interfaces. To understand the evolution of this regulatory system across the fornicate lineage, we have performed comparative genomic and phylogenetic analyses of the ARF GTPases, and their regulatory GAPs and GEFs in fornicate genomes and transcriptomes. Prior to our analysis of the fornicates, we first establish that the ARF GAP sub-family ArfGAP with dual PH domains (ADAP) is sparsely distributed but present in at least four eukaryotic supergroups and thus was likely present in the Last Eukaryotic Common Ancestor (LECA). Next, our collective comparative genomic and phylogenetic investigations into the ARF regulatory proteins in fornicates identify a duplication of ARF1 GTPase yielding two paralogues of ARF1F proteins, ancestral to all fornicates and present in all examined isolates of Giardia. However, the ARF GEF and ARF GAP complement is reduced compared with the LECA. This investigation shows that the system was significantly streamlined prior to the fornicate ancestor but was not further reduced concurrent with a transition into a parasitic lifestyle.

Suitability of current typing procedures to identify epidemiologically linked human Giardia duodenalis isolates

Background

Giardia duodenalis is a leading cause of gastroenteritis worldwide. Humans are mainly infected by two different subtypes, i.e., assemblage A and B. Genotyping is hampered by allelic sequence heterozygosity (ASH) mainly in assemblage B, and by occurrence of mixed infections. Here we assessed the suitability of current genotyping protocols of G. duodenalis for epidemiological applications such as molecular tracing of transmission chains.

Methodology/Principal findings

Two G. duodenalis isolate collections, from an outpatient tropical medicine clinic and from several primary care laboratories, were characterized by assemblage-specific qPCR (TIF, CATH gene loci) and a common multi locus sequence typing (MLST; TPI, BG, GDH gene loci). Assemblage A isolates were further typed at additional loci (HCMP22547, CID1, RHP26, HCMP6372, DIS3, NEK15411).

Of 175/202 (86.6%) patients the G. duodenalis assemblage could be identified: Assemblages A 25/175 (14.3%), B 115/175 (65.7%) and A+B mixed 35/175 (20.0%). By incorporating allelic sequence heterozygosity in the analysis, the three marker MLST correctly identified 6/9 (66,7%) and 4/5 (80.0%) consecutive samples from chronic assemblage B infections in the two collections, respectively, and identified a cluster of five independent patients carrying assemblage B parasites of identical MLST type. Extended MLST for assemblage A altogether identified 5/6 (83,3%) consecutive samples from chronic assemblage A infections and 15 novel genotypes. Based on the observed A+B mixed infections it is estimated that only 75% and 50% of assemblage A or B only cases represent single strain infections, respectively. We demonstrate that typing results are consistent with this prediction.

Conclusions/Significance

Typing of assemblage A and B isolates with resolution for epidemiological applications is possible but requires separate genotyping protocols. The high frequency of multiple infections and their impact on typing results are findings with immediate consequences for result interpretation in this field.

Giardia duodenalis is a leading cause of gastroenteritis worldwide. Humans are mainly infected by the two different genetic subtypes, assemblage A and B. Molecular typing tools for epidemiological applications such as tracking transmission, attribution to a source and outbreak investigations have been developed and are highly desirable. However, to what degree the tetraploid genome with allelic sequence heterogeneity (ASH), and the frequent occurrence of mixed, assemblage A and B infections hamper performance is unclear. Here, we assessed the suitability of current genotyping protocols for deciphering the molecular epidemiology of G. duodenalis. Against a common reporting bias, we incorporated ASH in the analysis and we show that typing with resolution for epidemiological applications is possible for both, assemblage A and B isolates, but requires separate protocols. We also demonstrate how the high frequency of multiple infections overall impacts on typing results, which has immediate consequences for result interpretation in this field.

Sirtuin GdSir2.4 participates in the regulation of rRNA transcription in the Giardia duodenalis parasite

Giardia duodenalis is a parasite of great medical interest due to the number of infections it causes worldwide each year. Although research on epigenetic mechanisms in this protist has only begun recently, epigenetic regulation has already been shown to have important roles in encystation, antigenic variation, and resistance to antibiotics in Giardia. In this work, we show that a Giardia ortholog of Sir2, GdSir2.4, is involved in the silencing of rRNA expression. Our results demonstrate that GdSir2.4 localizes to the nucleolus, and its binding to the intergenic spacer region of the rDNA is associated with the deacetylation of the chromatin in this region. Given the importance of the regulation of rRNA expression to maintain adequate levels of ribosomes and genomic stability within the cells, GdSir2.4 can be considered a target to create new therapeutic agents against this parasite.

Global Lysine Acetylation and 2-Hydroxyisobutyrylation Profiling Reveals the Metabolism Conversion Mechanism in Giardia lamblia

Giardia lamblia (G. lamblia) is the cause of giardiasis, a common infection that affects the general population of the world. Despite the constant possibility of damage because of their own metabolism, G. lamblia has survived and evolved to adapt to various environments. However, research on energy-metabolism conversion in G. lamblia is limited. This study aimed to reveal the dynamic metabolism conversion mechanism in G. lamblia under sugar starvation by detecting global lysine acetylation (Kac) and 2-hydroxyisobutyrylation (Khib) sites combined with quantitative proteome analyses. A total of 2999 acetylation sites on 956 proteins and 8877 2-hydroxyisobutyryl sites on 1546 proteins were quantified under sugar starvation. Integrated Kac and Khib data revealed that modified proteins were associated with arginine biosynthesis, glycolysis/gluconeogenesis, and alanine, aspartate, and glutamate metabolisms. These findings suggest that Kac and Khib were ubiquitous and provide deep insight into the metabolism conversion mechanism in G. lamblia under sugar starvation. Overall, these results can help delineate the biology of G. lamblia infections and reveal the evolutionary rule from prokaryote to eukaryote.

Transfection With Plasmid Causing Stable Expression of a Foreign Gene Affects General Proteome Pattern in Giardia lamblia Trophozoites

Giardia lamblia is an important causative agent of persistent diarrhea in humans, domestic animals, and cattle. Basic research is usually performed with the strain WBC6 and includes genetic manipulations such as transfections. Here, we investigate how transfection with a plasmid causing stable expression of a foreign gene affects the whole proteome pattern. Using shotgun mass spectrometry, we compare the proteomes of untransfected trophozoites to trophozoites transfected with Escherichia coli glucuronidase A (GusA). Besides GusA, which is detected in the transfected trophozoites only, the proteomes of untransfected and transfected trophozoites differ by 132 differentially expressed proteins. In particular, transfection induces antigenic variation. Since transfection causing stable expression affects the proteome pattern, transfection experiments should take into account this effect. Due to a unique peptide panel, GusA is an example for a suitable internal standard for experiments involving transfected cells. Data are available via ProteomeXchange with identifier PXD022565.

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.

High Cysteine Membrane Proteins (HCMPs) Are Up-Regulated During Giardia-Host Cell Interactions

Giardia intestinalis colonizes the upper small intestine of humans and animals, causing the diarrheal disease giardiasis. This unicellular eukaryotic parasite is not invasive but it attaches to the surface of small intestinal epithelial cells (IECs), disrupting the epithelial barrier. Here, we used an in vitro model of the parasite’s interaction with host IECs (differentiated Caco-2 cells) and RNA sequencing (RNAseq) to identify differentially expressed genes (DEGs) in Giardia, which might relate to the establishment of infection and disease induction. Giardia trophozoites interacted with differentiated Caco-2 cells for 1.5, 3, and 4.5 h and at each time point, 61, 89, and 148 parasite genes were up-regulated more than twofold, whereas 209, 265, and 313 parasite genes were down-regulated more than twofold. The most abundant DEGs encode hypothetical proteins and members of the High Cysteine Membrane Protein (HCMP) family. Among the up-regulated genes we also observed proteins associated with proteolysis, cellular redox balance, as well as lipid and nucleic acid metabolic pathways. In contrast, genes encoding kinases, regulators of the cell cycle and arginine metabolism and cytoskeletal proteins were down-regulated. Immunofluorescence imaging of selected, up-regulated HCMPs, using C-terminal HA-tagging, showed localization to the plasma membrane and peripheral vesicles (PVs). The expression of the HCMPs was affected by histone acetylation and free iron-levels. In fact, the latter was shown to regulate the expression of many putative giardial virulence factors in subsequent RNAseq experiments. We suggest that the plasma membrane localized and differentially expressed HCMPs play important roles during Giardia-host cell interactions.

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.

Comparative proteomics of three Giardia lamblia strains: investigation of antigenic variation in the post-genomic era

Giardia lamblia is a causative agent of persistent diarrhoea widespread in regions with low hygienic standards. Laboratory research is based on cloned lines issuing from various patient isolates typed in the late 1980s and 90s using restriction analysis and serology. In the present study, we compared the well-characterized strain WBC6 with another clone of the parent WB isolate termed WBA1 and with a clone from another isolate, GS/M-83-H7, using shotgun mass spectrometry proteomics. We identified 398 proteins differentially expressed between the GS and both WB isolates and 97 proteins differentially expressed between the two WB isolates. We investigated the expression levels of the predominant variant-specific surface proteins (VSPs) in each clone and matched the previously described major VSPs of each strain to the corresponding open reading frame sequences identified by whole-genome sequencing efforts. Furthermore, since the original WB isolate comes from a patient treated with metronidazole, we compared the susceptibilities of the strains to nitro compounds, as well the expression levels of enzymes involved in nitro reduction and on the corresponding enzyme activities and found distinct differences between the three strains.