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.