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

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.

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.

Phylogenetic Analysis of Pyruvate-Ferredoxin Oxidoreductase, a Redox Enzyme Involved in the Pharmacological Activation of Nitro-Based Prodrugs in Bacteria and Protozoa

The present frontrunners in the chemotherapy of infections caused by protozoa are nitro-based prodrugs that are selectively activated by PFOR-mediated redox reactions. This study seeks to analyze the distribution of PFOR in selected protozoa and bacteria by applying comparative genomics to test the hypothesis that PFOR in eukaryotes was acquired through horizontal gene transfer (HGT) from bacteria. Furthermore, to identify other putatively acquired genes, proteome-wide and gene enrichment analyses were used. A plausible explanation for the patchy occurrence of PFOR in protozoa is based on the hypothesis that bacteria are potential sources of genes that enhance the adaptation of protozoa in hostile environments. Comparative genomics of Entamoeba histolytica and the putative gene donor, Desulfovibrio vulgaris, identified eleven candidate genes for HGT involved in intermediary metabolism. If these results can be reproduced in other PFOR-possessing protozoa, it would provide more validated evidence to support the horizontal transfer of pfor from bacteria.

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.

Genetic diversity in the metronidazole metabolism genes nitroreductases and pyruvate ferredoxin oxidoreductases in susceptible and refractory clinical samples of Giardia lamblia

The effectiveness of metronidazole against the tetraploid intestinal parasite Giardia lamblia is dependent on its activation/inactivation within the cytoplasm. There are several activating enzymes, including pyruvate ferredoxin reductase (PFOR) and nitroreductase (NR) 1 which metabolize metronidazole into toxic forms, while NR2 on the other hand inactivates it. Metronidazole treatment failures have been increasing rapidly over the last decade, indicating genetic resistance mechanisms. Analyzing genetic variation in the PFOR and NR genes in susceptible and refractory Giardia isolates may help identify potential markers of resistance. Full length PFOR1, PFOR2, NR1 and NR2 genes from clinical culturable isolates and non-cultured clinical Giardia assemblage B samples were cloned, sequenced and single nucleotide variants (SNVs) were analyzed to assess genetic diversity and alleles. A similar ratio of amino acid changing SNVs per gene length was found for the NRs; 4.2% for NR1 and 6.4% for NR2, while the PFOR1 and PFOR2 genes had less variability with a ratio of 1.1% and 1.6%, respectively. One of the samples from a refractory case had a nonsense mutation which caused a truncated NR1 gene in one out of six alleles. Further, we found three NR2 alleles with frameshift mutations, possibly causing a truncated protein in two susceptible isolates. One of these isolates was homozygous for the affected NR2 allele. Three nsSNVs with potential for affecting protein function were found in the ferredoxin domain of the PFOR2 gene. The considerable variation and discovery of mutations possibly causing dysfunctional NR proteins in clinical Giardia assemblage B isolates, reveal a potential for genetic link to metronidazole susceptibility and resistance.

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.

Dissecting the Transcriptomes of Multiple Metronidazole-Resistant and Sensitive Trichomonas vaginalis Strains Identified Distinct Genes and Pathways Associated with Drug Resistance and Cell Death

Trichomonas vaginalis is the causative agent of trichomoniasis, the most prevalent non-viral sexually transmitted infection worldwide. Metronidazole (MTZ) is the mainstay of anti-trichomonal chemotherapy; however, drug resistance has become an increasingly worrying issue. Additionally, the molecular events of MTZ-induced cell death in T. vaginalis remain elusive. To gain insight into the differential expression of genes related to MTZ resistance and cell death, we conducted RNA-sequencing of three paired MTZ-resistant (MTZ-R) and MTZ-sensitive (MTZ-S) T. vaginalis strains treated with or without MTZ. Comparative transcriptomes analysis identified that several putative drug-resistant genes were exclusively upregulated in different MTZ-R strains, such as ATP-binding cassette (ABC) transporters and multidrug resistance pumps. Additionally, several shared upregulated genes among all the MTZ-R transcriptomes were not previously identified in T. vaginalis, such as 5′-nucleotidase surE and Na⁺-driven multidrug efflux pump, which are a potential stress response protein and a multidrug and toxic compound extrusion (MATE)-like protein, respectively. Functional enrichment analysis revealed that purine and pyrimidine metabolisms were suppressed in MTZ-S parasites upon drug treatment, whereas the endoplasmic reticulum-associated degradation (ERAD) pathway, proteasome, and ubiquitin-mediated proteolysis were strikingly activated, highlighting the novel pathways responsible for drug-induced stress. Our work presents the most detailed analysis of the transcriptional changes and the regulatory networks associated with MTZ resistance and MTZ-induced signaling, providing insights into MTZ resistance and cell death mechanisms in trichomonads.

In vitro selection of Giardia duodenalis for Albendazole resistance identifies a β-tubulin mutation at amino acid E198K

Benzimidazole-2-carbamate (BZ) compounds, including Albendazole (Alb), are one of just two drug classes approved to treat the gastrointestinal protist Giardia duodenalis. Benzimidazoles bind to the tubulin dimer interface overlapping the colchicine binding site (CBS) of β-tubulin, thereby inhibiting microtubule polymerisation and disrupting microtubule networks. These BZ compounds are widely used as anthelmintic, anti-fungal and anti-giardial drugs. However, in helminths and fungi, BZ-resistance is widespread and caused by specific point mutations primarily occurring at F167, E198 and F200 in β-tubulin isoform 1. BZ-resistance in Giardia is reported clinically and readily generated in vitro, with significant implications for Giardia control. In Giardia, BZ mode of action (MOA) and resistance mechanisms are presumed but not proven, and no mutations in β-tubulin have been reported in association with Alb resistance (AlbR). Herein, we undertook detailed in vitro drug-susceptibility screens of 13 BZ compounds and 7 Alb structural analogues in isogenic G. duodenalis isolates selected for AlbR and podophyllotoxin, another β-tubulin inhibitor, as well as explored cross-resistance to structurally unrelated, metronidazole (Mtz). AlbR lines exhibited co-resistance to many structural variants in the BZ-pharmacophore, and cross-resistance to podophyllotoxin. AlbR lines were not cross-resistant to Mtz, but MtzR lines had enhanced survival in Alb. Lastly, Alb analogues with longer thioether substituents had decreased potency against our AlbR lines. In silico modelling indicated the Alb-β-tubulin interaction in Giardia partially overlaps the CBS and corresponds to residues associated with BZ-resistance in helminths and fungi (F167, E198, F200). Sequencing of Giardia β-tubulin identified a single nucleotide polymorphism resulting in a mutation from glutamic acid to lysine at amino acid 198 (E198K). To our knowledge, this is the first β-tubulin mutation reported for protistan BZ-resistance. This study provides insight into BZ mode of action and resistance in Giardia, and presents a potential avenue for a genetic test for clinically resistance isolates.

Antimicrobial resistance of the enteric protozoon Giardia duodenalis - A narrative review

Introduction

As therapy-refractory giardiasis is an emerging health issue, this review aimed at summarizing mechanisms of reduced antimicrobial susceptibility in Giardia duodenalis and strategies to overcome this problem.

Methods

A narrative review on antimicrobial resistance in G. duodenalis was based upon a selective literature research.

Results

Failed therapeutic success has been observed for all standard therapies of giardiasis comprising nitroimidazoles like metronidazole or tinidazole as first line substances but also benznidazoles like albendazole and mebendazole, the nitrofuran furazolidone, the thiazolide nitazoxanide, and the aminoglycoside paromomycin. Multicausality of the resistance phenotypes has been described, with differentiated gene expression due to epigenetic and post-translational modifications playing a considerable bigger role than mutational base exchanges in the parasite DNA. Standardized resistance testing algorithms are not available and clinical evidence for salvage therapies is scarce in spite of research efforts targeting new giardicidal drugs.

Conclusion

In case of therapeutic failure of first line nitroimidazoles, salvage strategies including various options for combination therapy exist in spite of limited evidence and lacking routine diagnostic-compatible assays for antimicrobial susceptibility testing in G. duodenalis. Sufficiently powered clinical and diagnostic studies are needed to overcome both the lacking evidence regarding salvage therapy and the diagnostic neglect of antimicrobial resistance.

Giardiasis treatment: an update with a focus on refractory disease

PURPOSE OF REVIEW

Giardiasis remains a common cause of diarrhea and intestinal enteropathy globally. Here we give an overview of clinical treatment studies and discuss potential mechanisms and molecular targets for in-vitro testing of drug resistance.

RECENT FINDINGS

Giardia is a cause of disease both in diarrheal and nondiarrheal cases. The prevalence of treatment refractory giardiasis is increasing. Recent studies reveal 5-nitroimidazole refractory infection occurs in up to 50% of cases. Mechanisms of drug resistance are not known. Placebo controlled studies of drug efficacy, taking the self-limiting course of giardiasis into account, has not been reported. No randomized controlled trials of treatment of refractory infection have been performed the last 25 years. Based on the clinical studies reported, combination treatment with a 5-nitroimidazole and a benzimidazole is more effective than repeated courses of 5-nitroimidazole or monotherapies in refractory cases. Quinacrine is effective in refractory cases, but potentially severe side effects limit its use.

SUMMARY

A combination of a 5-nitroimidazole and albendazole or mebendazole, and quinacrine monotherapy, are rational choices in nitroimidazole refractory infections, but randomized controlled studies are needed. Further research into more recent clinical isolates is necessary to uncover mechanisms for the increase in metronidazole refractory giardiasis observed during the last decade.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis-a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

Genetic Diversity of the Flavohemoprotein Gene of Giardia lamblia: Evidence for High Allelic Heterozygosity and Copy Number Variation

Purpose

The flavohemoprotein (gFlHb) in Giardia plays an important role in managing nitrosative and oxidative stress, and potentially also in virulence and nitroimidazole drug tolerance. The aim of this study was to analyze the genetic diversity of gFlHb in Giardia assemblages A and B clinical isolates.

Methods

gFlHb genes from 20 cultured clinical Giardia isolates were subjected to PCR amplification and cloning, followed by Sanger sequencing. Sequences of all cloned PCR fragments from each isolate were analyzed for single nucleotide variants (SNVs) and compared to genomic Illumina sequence data. Identical clone sequences were sorted into alleles, and diversity was further analyzed. The number of gFlHb gene copies was assessed by mining PacBio de novo assembled genomes in eight isolates. Homology models for assessment of SNV's potential impact on protein function were created using Phyre2.

Results

A variable copy number of the gFlHb gene, between two and six copies, depending on isolate, was found. A total of 37 distinct sequences, representing different alleles of the gFlHb gene, were identified in AII isolates, and 41 were identified in B isolates. In some isolates, up to 12 different alleles were found. The total allelic diversity was high for both assemblages (>0.9) and was coupled with a nucleotide diversity of <0.01. The genetic variation (SNVs per CDS length) was 4.8% in sub-assemblage AII and 5.4% in assemblage B. The number of non-synonymous (ns) SNVs was high in gFIHb of both assemblages, 1.6% in A and 3.0% in B, respectively. Some of the identified nsSNV are predicted to alter protein structure and possibly function.

Conclusion

In this study, we present evidence that gFlHb, a putative protective enzyme against oxidative and nitrosative stress in Giardia, is a variable copy number gene with high allelic diversity. The genetic variability of gFlHb may contribute metabolic adaptability against metronidazole toxicity.

In vitro-induction of metronidazole-resistant Giardia duodenalis is not associated with nucleotide alterations in the genes involved in pro-drug activation

Giardiasis is an infectious disease caused by Giardia duodenalis. The pro-drug metronidazole (MTZ) is the first-line treatment for giardiasis. Parasite’s proteins as pyruvate:ferredoxin oxidoreductase (PFOR), ferredoxin (Fd), nitroreductase-1 (NR-1) and thioredoxin reductase (TrxR) participate in MTZ activation. Here, we showed Giardia trophozoites long-term exposed to MTZ presented higher IC₅₀ than controls, showing the drug influenced the parasite survival. That reduction in MTZ’s susceptibility does not seem to be related to mutations in the genes pfor, fd, nr-1 or trxr. It points that different mechanism as alterations in other metabolic pathways can account for Giardia resistance to MTZ therapy.

Whole-genome sequencing of dog-specific assemblages C and D of Giardia duodenalis from single and pooled cysts indicates host-associated genes

Giardia duodenalis (syn. Giardia intestinalis or Giardia lamblia) infSAects over 280 million people each year and numerous animals. G. duodenalis can be subdivided into eight assemblages with different host specificity. Unculturable assemblages have so far resisted genome sequencing efforts. In this study, we isolated single and pooled cysts of assemblages C and D from dog faeces by FACS, and sequenced them using multiple displacement amplification and Illumina paired-end sequencing. The genomes of assemblages C and D were compared with genomes of assemblages A and B from humans and assemblage E from ruminants and pigs. The genomes obtained from the pooled cysts and from the single cysts were considered complete (>99 % marker genes observed) and the allelic sequence heterozygosity (ASH) values of assemblages C and D were 0.89 and 0.74 %, respectively. These ASH values were slightly higher than for assemblage B (>0.43 %) and much higher than for assemblages A and E, which ranged from 0.002 to 0.037 %. The flavohaemoglobin and 4Fe-4S binding domain family encoding genes involved in O₂ and NO detoxification were only present in assemblages A, B and E. Cathepsin B orthologs were found in all genomes. Six clades of cathepsin B orthologs contained one gene of each genome, while in three clades not all assemblages were represented. We conclude that whole-genome sequencing from a single Giardia cyst results in complete draft genomes, making the genomes of unculturable Giardia assemblages accessible. Observed differences between the genomes of assemblages C and D on one hand and the assemblages A, B and E on the other hand are possibly associated with host specificity.

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.

Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy

Giardia lamblia, a causative agent of persistent diarrhea in humans, domestic animals, and cattle, is usually treated with nitro compounds. Consequently, enzymes involved in anaerobic nitro reduction have been investigated in detail as potential targets. Their role within the normal metabolic context is, however, not understood. Using ¹H high-resolution magic angle spinning (HR-MAS) NMR spectroscopy, we analyzed the metabolomes of G. lamblia trophozoites overexpressing three nitroreductases (NR1–NR3) and thioredoxin reductase (TrxR), most likely a scavenger of reactive oxygen species, as suggested by the results published in this study. We compared the patterns to convenient controls and to the situation in the nitro drug resistant strain C4 where NR1 is downregulated. We identified 27 metabolites in G. lamblia trophozoites. Excluding metabolites of high variability among different wildtype populations, only trophozoites overexpressing NR1 presented a distinct pattern of nine metabolites, in particular arginine catabolites, differing from the respective controls. This pattern matched a differential pattern between wildtype and strain C4. This suggests that NR1 interferes with arginine and thus energy metabolism. The exact metabolic function of NR1 (and the other nitroreductases) remains to be elucidated.

Drug resistance in Giardia: Mechanisms and alternative treatments for Giardiasis

The use of chemotherapeutic drugs is the main resource against clinical giardiasis due to the lack of approved vaccines. Resistance of G. duodenalis to the most used drugs to treat giardiasis, metronidazole and albendazole, is a clinical issue of growing concern and yet unknown impact, respectively. In the search of new drugs, the completion of the Giardia genome project and the use of biochemical, molecular and bioinformatics tools allowed the identification of ligands/inhibitors for about one tenth of ≈150 potential drug targets in this parasite. Further, the synthesis of second generation nitroimidazoles and benzimidazoles along with high-throughput technologies have allowed not only to define overall mechanisms of resistance to metronidazole but to screen libraries of repurposed drugs and new pharmacophores, thereby increasing the known arsenal of anti-giardial compounds to some hundreds, with most demonstrating activity against metronidazole or albendazole-resistant Giardia. In particular, cysteine-modifying agents which include omeprazole, disulfiram, allicin and auranofin outstand due to their pleiotropic activity based on the extensive repertoire of thiol-containing proteins and the microaerophilic metabolism of this parasite. Other promising agents derived from higher organisms including phytochemicals, lactoferrin and propolis as well as probiotic bacteria/fungi have also demonstrated significant potential for therapeutic and prophylactic purposes in giardiasis. In this context the present chapter offers a comprehensive review of the current knowledge, including commonly prescribed drugs, causes of therapeutic failures, drug resistance mechanisms, strategies for the discovery of new agents and alternative drug therapies.

On the molecular and cellular effects of omeprazole to further support its effectiveness as an antigiardial drug

Research on Giardia lamblia has accumulated large information about its molecular cell biology and infection biology. However, giardiasis is still one of the commonest parasitic diarrheal diseases affecting humans. Additionally, an alarming increase in cases refractory to conventional treatment has been reported in low prevalence settings. Consequently, efforts directed toward supporting the efficient use of alternative drugs, and the study of their molecular targets appears promising. Repurposing of proton pump inhibitors is effective in vitro against the parasite and the toxic activity is associated with the inhibition of the G. lamblia triosephosphate isomerase (GlTIM) via the formation of covalent adducts with cysteine residue at position 222. Herein, we evaluate the effectiveness of omeprazole in vitro and in situ on GlTIM mutants lacking the most superficial cysteines. We studied the influence on the glycolysis of Giardia trophozoites treated with omeprazole and characterized, for the first time, the morphological effect caused by this drug on the parasite. Our results support the effectiveness of omeprazole against GlTIM despite of the possibility to mutate the druggable amino acid targets as an adaptive response. Also, we further characterized the effect of omeprazole on trophozoites and discuss the possible mechanism involved in its antigiardial effect.

Genetic variation in metronidazole metabolism and oxidative stress pathways in clinical Giardia lamblia assemblage A and B isolates

Purpose: Treatment-refractory Giardia cases have increased rapidly within the last decade. No markers of resistance nor a standardized susceptibility test have been established yet, but several enzymes and their pathways have been associated with metronidazole (MTZ) resistant Giardia. Very limited data are available regarding genetic variation in these pathways. We aimed to investigate genetic variation in metabolic pathway genes proposed to be involved in MTZ resistance in recently acquired, cultured clinical isolates.

Methods: Whole genome sequencing of 12 assemblage A2 and 8 assemblage B isolates was done, to decipher genomic variation in Giardia. Twenty-nine genes were identified in a literature search and investigated for their single nucleotide variants (SNVs) in the coding/non-coding regions of the genes, either as amino acid changing (non-synonymous SNVs) or non-changing SNVs (synonymous).

Results: In Giardia assemblage B, several genes involved in MTZ activation or oxidative stress management were found to have higher numbers of non-synonymous SNVs (thioredoxin peroxidase, nitroreductase 1, ferredoxin 2, NADH oxidase, nitroreductase 2, alcohol dehydrogenase, ferredoxin 4 and ferredoxin 1) than the average variation. For Giardia assemblage A2, the highest genetic variability was found in the ferredoxin 2, ferredoxin 6 and in nicotinamide adenine dinucleotide phosphate (NADPH) oxidoreductase putative genes. SNVs found in the ferredoxins and nitroreductases were analyzed further by alignment and homology modeling. SNVs close to the iron-sulfur cluster binding sites in nitroreductase-1 and 2 and ferredoxin 2 and 4 could potentially affect protein function. Flavohemoprotein seems to be a variable-copy gene, due to higher, but variable coverage compared to other genes investigated.

Conclusion: In clinical Giardia isolates, genetic variability is common in important genes in the MTZ metabolizing pathway and in the management of oxidative and nitrosative stress and includes high numbers of non-synonymous SNVs. Some of the identified amino acid changes could potentially affect the respective proteins important in the MTZ metabolism.

Resistance formation to nitro drugs in Giardia lamblia: No common markers identified by comparative proteomics

In order to elucidate the question whether resistance to nitro drugs in G. lamblia is due to common resistance markers, trophozoites of three resistant G. lamblia strains, namely C4, 1062ID10, and 713M3 were grown in the presence of the two nitro drugs metronidazole and nitazoxanide and compared to their corresponding wild-types WBC6, 106, and 713 by mass spectometry shotgun analysis of their proteomes. Depending on the strain and the nitro drug, more than 200 to 500 differentially expressed proteins were identified, but there were no common patterns across strains and drugs. All resistant strains underwent antigenic variation with distinct surface antigens like variant surface proteins or cysteine rich proteins depending on strain and nitro compound. A closer look on enzymes involved in nitroreduction and detoxification of nitro radicals, NO or O₂ suggested the existence of distinct strategies for each drug and each strain. Therefore, we conclude that resistance to nitro drugs in G. lamblia is not correlated with a specific pattern of differentially expressed proteins and therefore seems not to be the result of a directed process.

•

Is resistance to nitro drugs in G. lamblia due to common resistance markers?

•

Three resistant strains were grown in the presence of two nitro drugs separately and compared to wild-types by MS shotgun analysis.

•

More than 200 to 500 differentially expressed proteins identified depending on strain and drug.

•

No common patterns across strains and drugs.

•

Strain specific antigenic variation and strategies linked to nitro reduction.

Oxygen induces the expression of invasion and stress response genes in the anaerobic salmon parasite Spironucleus salmonicida

Background

Spironucleus salmonicida is an anaerobic parasite that can cause systemic infections in Atlantic salmon. Unlike other diplomonad parasites, such as the human pathogen Giardia intestinalis, Spironucleus species can infiltrate the blood stream of their hosts eventually colonizing organs, skin and gills. How this presumed anaerobe can persist and invade oxygenated tissues, despite having a strictly anaerobic metabolism, remains elusive.

Results

To investigate how S. salmonicida response to oxygen stress, we performed RNAseq transcriptomic analyses of cells grown in the presence of oxygen or antioxidant-free medium. We found that over 20% of the transcriptome is differentially regulated in oxygen (1705 genes) and antioxidant-depleted (2280 genes) conditions. These differentially regulated transcripts encode proteins related to anaerobic metabolism, cysteine and Fe-S cluster biosynthesis, as well as a large number of proteins of unknown function. S. salmonicida does not encode genes involved in the classical elements of oxygen metabolism (e.g., catalases, superoxide dismutase, glutathione biosynthesis, oxidative phosphorylation). Instead, we found that genes encoding bacterial-like oxidoreductases were upregulated in response to oxygen stress. Phylogenetic analysis revealed some of these oxygen-responsive genes (e.g., nadh oxidase, rubrerythrin, superoxide reductase) are rare in eukaryotes and likely derived from lateral gene transfer (LGT) events into diplomonads from prokaryotes. Unexpectedly, we observed that many host evasion- and invasion-related genes were also upregulated under oxidative stress suggesting that oxygen might be an important signal for pathogenesis.

Conclusion

While oxygen is toxic for related organisms, such as G. intestinalis, we find that oxygen is likely a gene induction signal for host invasion- and evasion-related pathways in S. salmonicida. These data provide the first molecular evidence for how S. salmonicida could tolerate oxic host environments and demonstrate how LGT can have a profound impact on the biology of anaerobic parasites.

Electronic supplementary material

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

Annotation of the Giardia proteome through structure-based homology and machine learning

Background

Large-scale computational prediction of protein structures represents a cost-effective alternative to empirical structure determination with particular promise for non-model organisms and neglected pathogens. Conventional sequence-based tools are insufficient to annotate the genomes of such divergent biological systems. Conversely, protein structure tolerates substantial variation in primary amino acid sequence and is thus a robust indicator of biochemical function. Structural proteomics is poised to become a standard part of pathogen genomics research; however, informatic methods are now required to assign confidence in large volumes of predicted structures.

Aims

Our aim was to predict the proteome of a neglected human pathogen, Giardia duodenalis, and stratify predicted structures into high- and lower-confidence categories using a variety of metrics in isolation and combination.

Methods

We used the I-TASSER suite to predict structural models for ∼5,000 proteins encoded in G. duodenalis and identify their closest empirically-determined structural homologues in the Protein Data Bank. Models were assigned to high- or lower-confidence categories depending on the presence of matching protein family (Pfam) domains in query and reference peptides. Metrics output from the suite and derived metrics were assessed for their ability to predict the high-confidence category individually, and in combination through development of a random forest classifier.

Results

We identified 1,095 high-confidence models including 212 hypothetical proteins. Amino acid identity between query and reference peptides was the greatest individual predictor of high-confidence status; however, the random forest classifier outperformed any metric in isolation (area under the receiver operating characteristic curve = 0.976) and identified a subset of 305 high-confidence-like models, corresponding to false-positive predictions. High-confidence models exhibited greater transcriptional abundance, and the classifier generalized across species, indicating the broad utility of this approach for automatically stratifying predicted structures. Additional structure-based clustering was used to cross-check confidence predictions in an expanded family of Nek kinases. Several high-confidence-like proteins yielded substantial new insight into mechanisms of redox balance in G. duodenalis-a system central to the efficacy of limited anti-giardial drugs.

Conclusion

Structural proteomics combined with machine learning can aid genome annotation for genetically divergent organisms, including human pathogens, and stratify predicted structures to promote efficient allocation of limited resources for experimental investigation.

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.

Treatment-refractory giardiasis: challenges and solutions

Giardia is the commonest parasitic diarrheal pathogen affecting humans and a frequent cause of waterborne/foodborne parasitic diseases worldwide. Prevalence of giardiasis is higher in children, living in poor, low hygiene settings in developing countries, and in travelers returning from highly endemic areas. The clinical picture of giardiasis is heterogeneous, with high variability in severity of clinical disease. It can become chronic or be followed by post-infectious sequelae. An alarming increase in cases refractory to the conventional treatment with nitroimidazoles (ie, metronidazole) has been reported in low prevalence settings, such as European Union countries, especially in patients returning from Asia. In view of its relevance, we aim in this review to recapitulate present clinical knowledge about Giardia, with a special focus on the challenge of treatment-refractory giardiasis. We propose a working definition of clinically drug-resistant giardiasis, summarize knowledge regarding resistance mechanisms, and discuss its clinical management according to research-based evidence and medical practice. Advances in development and identification of novel drugs and potential non-pharmacological alternatives are also reviewed with the overall aim to define knowledge gaps and suggest future directions for research.

Physiological aspects of nitro drug resistance in Giardia lamblia

For over 50 years, metronidazole and other nitro compounds such as nitazoxanide have been used as a therapy of choice against giardiasis and more and more frequently, resistance formation has been observed. Model systems allowing studies on biochemical aspects of resistance formation to nitro drugs are, however, scarce since resistant strains are often unstable in culture. In order to fill this gap, we have generated a stable metronidazole- and nitazoxanide-resistant Giardia lamblia WBC6 clone, the strain C4. Previous studies on strain C4 and the corresponding wild-type strain WBC6 revealed marked differences in the transcriptomes of both strains. Here, we present a physiological comparison between trophozoites of both strains with respect to their ultrastructure, whole cell activities such as oxygen consumption and resazurin reduction assays, key enzyme activities, and several metabolic key parameters such as NAD(P)⁺/NAD(P)H and ADP/ATP ratios and FAD contents. We show that nitro compound-resistant C4 trophozoites exhibit lower nitroreductase activities, lower oxygen consumption and resazurin reduction rates, lower ornithine-carbamyl-transferase activity, reduced FAD and NADP(H) pool sizes and higher ADP/ATP ratios than wildtype trophozoites. The present results suggest that resistance formation against nitro compounds is correlated with metabolic adaptations resulting in a reduction of the activities of FAD-dependent oxidoreductases.

Genetic Indicators of Drug Resistance in the Highly Repetitive Genome of Trichomonas vaginalis

Trichomonas vaginalis, the most common nonviral sexually transmitted parasite, causes ∼283 million trichomoniasis infections annually and is associated with pregnancy complications and increased risk of HIV-1 acquisition. The antimicrobial drug metronidazole is used for treatment, but in a fraction of clinical cases, the parasites can become resistant to this drug. We undertook sequencing of multiple clinical isolates and lab derived lines to identify genetic markers and mechanisms of metronidazole resistance. Reduced representation genome sequencing of ∼100 T. vaginalis clinical isolates identified 3,923 SNP markers and presence of a bipartite population structure. Linkage disequilibrium was found to decay rapidly, suggesting genome-wide recombination and the feasibility of genetic association studies in the parasite. We identified 72 SNPs associated with metronidazole resistance, and a comparison of SNPs within several lab-derived resistant lines revealed an overlap with the clinically resistant isolates. We identified SNPs in genes for which no function has yet been assigned, as well as in functionally-characterized genes relevant to drug resistance (e.g., pyruvate:ferredoxin oxidoreductase). Transcription profiles of resistant strains showed common changes in genes involved in drug activation (e.g., flavin reductase), accumulation (e.g., multidrug resistance pump), and detoxification (e.g., nitroreductase). Finally, we identified convergent genetic changes in lab-derived resistant lines of Tritrichomonas foetus, a distantly related species that causes venereal disease in cattle. Shared genetic changes within and between T. vaginalis and Tr. foetus parasites suggest conservation of the pathways through which adaptation has occurred. These findings extend our knowledge of drug resistance in the parasite, providing a panel of markers that can be used as a diagnostic tool.

A review on metronidazole: an old warhorse in antimicrobial chemotherapy

The 5-nitroimidazole drug metronidazole has remained the drug of choice in the treatment of anaerobic infections, parasitic as well as bacterial, ever since its development in 1959. In contrast to most other antimicrobials, it has a pleiotropic mode of action and reacts with a large number of molecules. Importantly, metronidazole, which is strictly speaking a prodrug, needs to be reduced at its nitro group in order to become toxic. Reduction of metronidazole, however, only takes place under very low concentrations of oxygen, explaining why metronidazole is exclusively toxic to microaerophilic and anaerobic microorganisms. In general, resistance rates amongst the pathogens treated with metronidazole have remained low until the present day. Nevertheless, metronidazole resistance does occur, and for the treatment of some pathogens, especially Helicobacter pylori, metronidazole has become almost useless in some parts of the world. This review will give an account on the current status of research on metronidazole's mode of action, metronidazole resistance in eukaryotes and prokaryotes, and on other 5-nitroimidazoles in use.

Nitroimidazoles: Molecular Fireworks That Combat a Broad Spectrum of Infectious Diseases

Infectious diseases claim millions of lives every year, but with the advent of drug resistance, therapeutic options to treat infections are inadequate. There is now an urgent need to develop new and effective treatments. Nitroimidazoles are a class of antimicrobial drugs that have remarkable broad spectrum activity against parasites, mycobacteria, and anaerobic Gram-positive and Gram-negative bacteria. While nitroimidazoles were discovered in the 1950s, there has been renewed interest in their therapeutic potential, particularly for the treatment of parasitic infections and tuberculosis. In this review, we summarize different classes of nitroimidazoles that have been described in the literature in the past five years, from approved drugs and clinical candidates to examples undergoing preclinical or early stage development. The relatively "nonspecific" mode of action and resistance mechanisms of nitromidazoles are discussed, and contemporary strategies to facilitate nitroimidazole drug development are highlighted.