Search

Christina S Saghaug,1,2 Christian Klotz,3 Juha P Kallio,4 Toni Aebischer,3 Nina Langeland,1,2,5 Kurt Hanevik1 1Department of Clinical Science, University of Bergen, Bergen, Norway; 2Norwegian National Advisory Unit on Tropical Infectious Diseases, Department of Medicine, Haukeland University Hospital, Bergen, Norway; 3Department of Infectious Diseases, Unit 16 Mycotic and Parasitic Agents and Mycobacteria, Robert Koch-Institute, Berlin, Germany; 4Department of Biomedicine, University of Bergen, Bergen, Norway; 5Department of Medicine, Haraldsplass Deaconess Hospital, Bergen, NorwayCorrespondence: Christina S SaghaugDepartment of Clinical Science, University of Bergen, 8th Floor, Lab-Building, Bergen N-5021, NorwayTel +47 90 13 24 14Email christina.saghaug@uib.noPurpose: 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.Keywords: Giardia, genetic diversity, copy number variation, flavohemoprotein, oxidative stress, nitrosative stress, allele

Christina S Saghaug,1,2 Christian Klotz,3 Juha P Kallio,4 Hans-Richard Brattbakk,1,5 Tomasz Stokowy,1,5 Toni Aebischer,3 Inari Kursula,4,6 Nina Langeland,1–2,7 Kurt Hanevik1,21Department of Clinical Science, University of Bergen, Bergen, Hordaland, Norway; 2Norwegian National Advisory Unit on Tropical Infectious Diseases, Department of Medicine, Haukeland University Hospital, Bergen, Hordaland, Norway; 3Department of Infectious Diseases, Unit 16 Mycotic and Parasitic Agents and Mycobacteria, Robert Koch-Institute, Berlin, Germany; 4Department of Biomedicine, University of Bergen, Bergen, Hordaland, Norway; 5Center for Medical Genetics and Molecular Medicine, Haukeland University Hospital, Bergen, Hordaland, Norway; 6Biocenter Oulu and Faculty of Biochemistry and Molecular Medicine, University of Oulu, Oulu, Finland; 7Department of Medicine, Haraldsplass Deaconess Hospital, Bergen, Hordaland, NorwayPurpose: 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.Keywords: drug metabolism, resistance, genetic analysis, metronidazole genes, ferredoxin, genetic diversity