1
|
Shevtsov A, Izbanova U, Amirgazin A, Kairzhanova A, Dauletov A, Kiyan V, Vergnaud G. Genetic Homogeneity of Francisella tularensis subsp. mediasiatica Strains in Kazakhstan. Pathogens 2024; 13:581. [PMID: 39057808 PMCID: PMC11279412 DOI: 10.3390/pathogens13070581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2024] [Revised: 07/05/2024] [Accepted: 07/06/2024] [Indexed: 07/28/2024] Open
Abstract
Tularemia is an acute febrile disease caused by the Gram-negative bacillus Francisella tularensis. Based on genetic and phenotypic characteristics, three subspecies are distinguished: tularensis, holarctica, and mediasiatica. F. tularensis subsp. mediasiatica remains the least studied subspecies. Over the past decade, new foci of distribution of F. tularensis subsp. mediasiatica have been discovered in Russia (Siberia), expanding the possible distribution area by thousands of kilometers. This article provides whole genome single nucleotide polymorphism (wgSNP) and polymorphic tandem repeats (MLVA) analyses of 28 mediasiatica strains isolated between 1965 and 2004 in Kazakhstan. Despite high genetic homogeneity, MLVA with eleven loci (MLVA11) demonstrates a high discriminatory ability (diversity index, 0.9497). The topological structure of the trees based on wgSNP and MLVA is not comparable; however, clustering remains congruent for most outbreaks, with the exception of two strains from one outbreak that are identical in terms of wgSNP but differ at three tandem repeat loci. Based on wgSNP, the strains are assigned to one of the three currently known mediasiatica sublineages, lineage M.I, together with other historical strains maintained in collections in Russia and Sweden. wgSNP shows limited previously unknown genetic diversity, with the M.I lineage size being only 118 SNPs. The wgSNP genotype is not strongly correlated with year and place of isolation.
Collapse
Affiliation(s)
- Alexandr Shevtsov
- National Center for Biotechnology, Astana 010000, Kazakhstan; (A.A.); (A.K.); (A.D.); (V.K.)
| | - Uinkul Izbanova
- Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty 050000, Kazakhstan
| | - Asylulan Amirgazin
- National Center for Biotechnology, Astana 010000, Kazakhstan; (A.A.); (A.K.); (A.D.); (V.K.)
| | - Alma Kairzhanova
- National Center for Biotechnology, Astana 010000, Kazakhstan; (A.A.); (A.K.); (A.D.); (V.K.)
| | - Ayan Dauletov
- National Center for Biotechnology, Astana 010000, Kazakhstan; (A.A.); (A.K.); (A.D.); (V.K.)
| | - Vladimir Kiyan
- National Center for Biotechnology, Astana 010000, Kazakhstan; (A.A.); (A.K.); (A.D.); (V.K.)
| | - Gilles Vergnaud
- Institute for Integrative Biology of the Cell (I2BC), Université Paris-Saclay, CEA, CNRS, 91198 Gif-sur-Yvette, France
| |
Collapse
|
2
|
Izbanova U, Lukhnova L, Sadovskaya V, Zhumadilova Z, Meka-Mechenko T, Shevtsov A, Baitursyn B, Turebekov N, Tukhanova N. Characterization of tularemia foci in the Republic of Kazakhstan from 2000 to 2020. FRONTIERS IN EPIDEMIOLOGY 2024; 4:1291690. [PMID: 38456077 PMCID: PMC10910903 DOI: 10.3389/fepid.2024.1291690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/10/2023] [Accepted: 02/06/2024] [Indexed: 03/09/2024]
Abstract
The wide distribution of tularemia in the territory of Kazakhstan is associated with landscape and geographical characteristics. This is explained by a combination of natural factors: the presence of certain types of rodents-reservoirs and sources, ectoparasites-carriers of the causative agent of tularemia. The study of the current spatial and temporal characterization of tularemia in Kazakhstan from 2000 to 2020 will determine the epidemiological status of tularemia and improve the monitoring system in Kazakhstan. In this work we demonstrated the results of a retrospective survey of natural foci of tularemia: analysis of vector, small mammal and human data. The spatial and temporal characteristics of tularemia from 2000 to 2020 in the territory of Kazakhstan were studied in comparison with historical data, including the description of tularemia outbreaks, the clinical picture, and the source of infection, transmission factors, and geographical coordinates of outbreak registration. Sampling was carried out by trapping rodents on snap traps and collecting ticks by rodent combing and by "flagging" methods. For the last 20 years, 85 human cases of tularemia have been reported. During the period from 2000 to 2020, more than 600 strains of F. tularensis were isolated from field rodents and ticks in the natural foci of tularemia. MLVA typing of F. tularensis strains isolated from natural foci of tularemia in Kazakhstan over the past 20 years. The results of retrospective monitoring indicate that currently active foci of tularemia include the Aktobe, West Kazakhstan, Almaty, East Kazakhstan, and Pavlodar regions. Low-activity natural foci are located in the territory of the Akmola, Karaganda, North Kazakhstan, Kostanay, Atyrau, Zhambyl, and Kyzylorda regions. There are no active natural foci of tularemia in the Mangystau and Turkestan regions. The widespread occurrence of tularemia in the country is associated with landscape and geographical features that contribute to the circulation of the pathogen in the natural focus. An analysis of natural foci of tularemia showed that it is necessary to continue monitoring studies of carriers and vectors for the presence of the causative agent of the F. tularensis, in order to prevent mass cases of human disease.
Collapse
Affiliation(s)
- U. Izbanova
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - L. Lukhnova
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - V. Sadovskaya
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - Z. Zhumadilova
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - T. Meka-Mechenko
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - A. Shevtsov
- National Center for Biotechnology, Astana, Kazakhstan
| | - B. Baitursyn
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - N. Turebekov
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| | - N. Tukhanova
- M.Aikimbayev’s National Scientific Center for Especially Dangerous Infections, Almaty, Kazakhstan
| |
Collapse
|
3
|
Isidro J, Escudero R, Luque-Larena JJ, Pinto M, Borges V, González-Martín-Niño R, Duarte S, Vieira L, Mougeot F, Vidal D, Herrera-Rodríguez D, Rodríguez-Pastor R, Herrero-Cófreces S, Jubete-Tazo F, Gomes JP, Lopes de Carvalho I. Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples. Front Microbiol 2024; 14:1277468. [PMID: 38249473 PMCID: PMC10797068 DOI: 10.3389/fmicb.2023.1277468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 11/23/2023] [Indexed: 01/23/2024] Open
Abstract
Introduction Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited. Methods To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity. Results We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries. Discussion Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.
Collapse
Affiliation(s)
- Joana Isidro
- Genomics and Bioinformatics Unit, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
| | - Raquel Escudero
- Reference and Research Laboratory on Special Pathogens, National Centre for Microbiology (CNM), Carlos II Health Institute (ISCIII), Madrid, Spain
| | - Juan José Luque-Larena
- Departamento de Ciencias Agroforestales, Instituto Universitario de Investigación en Gestión Forestal Sostenible (iuFOR), E.T.S. Ingenierías Agrarias, Universidad de Valladolid, Palencia, Spain
| | - Miguel Pinto
- Genomics and Bioinformatics Unit, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
| | - Vítor Borges
- Genomics and Bioinformatics Unit, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
| | - Rosa González-Martín-Niño
- Reference and Research Laboratory on Special Pathogens, National Centre for Microbiology (CNM), Carlos II Health Institute (ISCIII), Madrid, Spain
| | - Sílvia Duarte
- Technology and Innovation Unit, Department of Human Genetics, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
| | - Luís Vieira
- Technology and Innovation Unit, Department of Human Genetics, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
| | - François Mougeot
- Instituto de Investigación en Recursos Cinegéticos (IREC-CSIC, UCLM, JCCM), Ciudad Real, Spain
| | - Dolors Vidal
- Área de Microbiología, Facultad de Medicina, Universidad de Catilla-La Mancha (UCLM), Ciudad Real, Spain
| | - Daniel Herrera-Rodríguez
- Instituto de Investigación en Recursos Cinegéticos (IREC-CSIC, UCLM, JCCM), Ciudad Real, Spain
- Área de Microbiología, Facultad de Medicina, Universidad de Catilla-La Mancha (UCLM), Ciudad Real, Spain
| | - Ruth Rodríguez-Pastor
- Department of Parasitology, Veterinary Faculty, University of Zaragoza, Zaragoza, Spain, Zaragoza, Spain
- Departamento de Parasitología, Facultad de Veterinaria, Universidad de Zaragoza, Instituto Agroalimentario de Aragón-IA2 (Universidad de Zaragoza-CITA), Zaragoza, Spain
| | - Silvia Herrero-Cófreces
- Departamento de Ciencias Agroforestales, Instituto Universitario de Investigación en Gestión Forestal Sostenible (iuFOR), E.T.S. Ingenierías Agrarias, Universidad de Valladolid, Palencia, Spain
| | - Fernando Jubete-Tazo
- Departamento de Ciencias Agroforestales, Instituto Universitario de Investigación en Gestión Forestal Sostenible (iuFOR), E.T.S. Ingenierías Agrarias, Universidad de Valladolid, Palencia, Spain
| | - João Paulo Gomes
- Genomics and Bioinformatics Unit, National Institute of Health Doutor Ricardo Jorge (INSA), Lisbon, Portugal
- Veterinary and Animal Research Center (CECAV), Faculty of Veterinary Medicine, Lusófona University, Lisbon, Portugal
| | - Isabel Lopes de Carvalho
- Emergency and Biopreparedness Unit, National Institute of Health Doutor Ricardo Jorge, Lisbon, Portugal
| |
Collapse
|
4
|
Köppen K, Rydzewski K, Doellinger J, Myrtennäs K, Forsman M, Appelt S, Scholz H, Heuner K. Phenotypic and genotypic discrimination of Francisella tularensis ssp. holarctica clades. Int J Med Microbiol 2023; 313:151583. [PMID: 37331050 DOI: 10.1016/j.ijmm.2023.151583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 06/12/2023] [Accepted: 06/13/2023] [Indexed: 06/20/2023] Open
Abstract
Francisella tularensis is the causative agent of tularemia, a zoonotic disease with a wide host range. F. tularensis ssp. holarctica (Fth) is of clinical relevance for European countries, including Germany. Whole genome sequencing methods, including canonical Single Nucleotide Polymorphism (canSNP) typing and whole genome SNP typing, have revealed that European Fth strains belong to a few monophyletic populations. The majority of German Fth isolates belong to two basal phylogenetic clades B.6 (biovar I) and B.12 (biovar II). Strains of B.6 and B.12 seem to differ in their pathogenicity, and it has been shown that strains of biovar II are resistant against erythromycin. In this study, we present data corroborating our previous data demonstrating that basal clade B.12 can be divided into clades B.71 and B.72. By applying phylogenetic whole genome analysis as well as proteome analysis, we could verify that strains of these two clades are distinct from one another. This was confirmed by measuring the intensity of backscatter light on bacteria grown in liquid media. Strains belonging to clades B.6, B.71 or B.72 showed clade-specific backscatter growth curves. Furthermore, we present the whole genome sequence of strain A-1341, as a reference genome of clade B.71, and whole proteomes comparison of Fth strains belonging to clades B.6, B.71 and B.72. Further research is necessary to investigate phenotypes and putative differences in pathogenicity of the investigated different clades of Fth to better understand the relationship between observed phenotypes, pathogenicity and distribution of Fth strains.
Collapse
Affiliation(s)
- Kristin Köppen
- Working group: Cellular Interactions of Bacterial Pathogens, Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany
| | - Kerstin Rydzewski
- Working group: Cellular Interactions of Bacterial Pathogens, Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany
| | - Joerg Doellinger
- Centre for Biological Threats and Special Pathogens, Proteomics and Spectroscopy (ZBS 6), Robert Koch Institute, Berlin, Germany
| | - Kerstin Myrtennäs
- Division of CBRN Defence and Security, Swedish Defence Research Agency (FOI), Umeå, Sweden
| | - Mats Forsman
- Division of CBRN Defence and Security, Swedish Defence Research Agency (FOI), Umeå, Sweden
| | - Sandra Appelt
- Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany
| | - Holger Scholz
- Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany
| | - Klaus Heuner
- Working group: Cellular Interactions of Bacterial Pathogens, Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany; Centre for Biological Threats and Special Pathogens, Highly Pathogenic Microorganisms (ZBS 2), Robert Koch Institute, Berlin, Germany.
| |
Collapse
|
5
|
Rijks JM, Tulen AD, Notermans DW, Reubsaet FAG, de Vries MC, Koene MGJ, Swaan CM, Maas M. Tularemia Transmission to Humans, the Netherlands, 2011–2021. Emerg Infect Dis 2022; 28:883-885. [PMID: 35318921 PMCID: PMC8962900 DOI: 10.3201/eid2804.211913] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
We used national registry data on human cases of Francisella tularensis subspecies holarctica infection to assess transmission modes among all 26 autochthonous cases in the Netherlands since 2011. The results indicate predominance of terrestrial over aquatic animal transmission sources. We recommend targeting disease-risk communication toward hunters, recreationists, and outdoor professionals.
Collapse
|
6
|
Shevtsov V, Kairzhanova A, Shevtsov A, Shustov A, Kalendar R, Abdrakhmanov S, Lukhnova L, Izbanova U, Ramankulov Y, Vergnaud G. Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan. PLoS Negl Trop Dis 2021; 15:e0009419. [PMID: 33999916 PMCID: PMC8158875 DOI: 10.1371/journal.pntd.0009419] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 05/27/2021] [Accepted: 04/29/2021] [Indexed: 11/29/2022] Open
Abstract
Tularemia is a highly dangerous zoonotic infection due to the bacteria Francisella tularensis. Low genetic diversity promoted the use of polymorphic tandem repeats (MLVA) as first-line assay for genetic description. Whole genome sequencing (WGS) is becoming increasingly accessible, opening the perspective of a time when WGS might become the universal genotyping assay. The main goal of this study was to describe F. tularensis strains circulating in Kazakhstan based on WGS data and develop a MLVA assay compatible with in vitro and in silico analysis. In vitro MLVA genotyping and WGS were performed for the vaccine strain and for 38 strains isolated in Kazakhstan from natural water bodies, ticks, rodents, carnivores, and from one migratory bird, an Isabellina wheatear captured in a rodent burrow. The two genotyping approaches were congruent and allowed to attribute all strains to two F. tularensis holarctica lineages, B.4 and B.12. The seven tandem repeats polymorphic in the investigated strain collection could be typed in a single multiplex PCR assay. Identical MLVA genotypes were produced by in vitro and in silico analysis, demonstrating full compatibility between the two approaches. The strains from Kazakhstan were compared to all publicly available WGS data of worldwide origin by whole genome SNP (wgSNP) analysis. Genotypes differing at a single SNP position were collected within a time interval of more than fifty years, from locations separated from each other by more than one thousand kilometers, supporting a role for migratory birds in the worldwide spread of the bacteria. Genotyping of Francisella tularensis has become a routine practice in epidemiology. Despite rapidly accumulating knowledge, the phylogeography of the pathogen is still poorly understood and discussions about geographic and temporal origins continue. One important reason is the poor characterization of the pathogen in many tularemia-endemic countries. This article describes the genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan using tandem repeat polymorphisms as well as whole genome sequencing. Thirty-nine strains were analyzed and two lineages were identified, namely B.4 and B.12. The study demonstrates a wider distribution of genotype B.4 in Asia, and identified a more basal branching point in this subclade. The obtained data support the Asian origin hypothesis for F. tularensis. The finding of identical genotypes in strains separated in time by decades and a thousand-kilometers geographic distance, confirms the ability of the bacteria for long-term preservation and fast long distances spread. The isolation of F. tularensis subsp. holarctica from the bird species Isabellina wheatear allows speculating about a major contribution of birds to the phylogeography of the pathogen. A genotyping protocol was developed utilizing seven polymorphic tandem repeats, two of which were identified within the framework of this work. The in vitro and in silico results are identical when using sequencing reads of 300 base-pairs or more.
Collapse
Affiliation(s)
| | - Alma Kairzhanova
- National Center for Biotechnology, Nur Sultan, Kazakhstan
- S. Seifullin Kazakh Agrotechnical University, Nur Sultan, Kazakhstan
| | - Alexandr Shevtsov
- National Center for Biotechnology, Nur Sultan, Kazakhstan
- * E-mail: (AS); (GV)
| | | | | | | | - Larissa Lukhnova
- National Scientific Center for Especially Dangerous Infections named by Masgut Aykimbayev, Almaty, Kazakhstan
| | - Uinkul Izbanova
- National Scientific Center for Especially Dangerous Infections named by Masgut Aykimbayev, Almaty, Kazakhstan
| | - Yerlan Ramankulov
- National Center for Biotechnology, Nur Sultan, Kazakhstan
- School of Science and Technology Nazarbayev University, Nur Sultan, Kazakhstan
| | - Gilles Vergnaud
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell, Gif-sur-Yvette, France
- * E-mail: (AS); (GV)
| |
Collapse
|
7
|
Glinšek Biškup U, Kogoj R, Korva M, Knap N, Cerar T, Knapič T, Petrovec M, Avšič-Županc T. Characterization of Tularemia Cases in Slovenia with Multiple-Locus Variable-Number Tandem Repeat Analysis. Vector Borne Zoonotic Dis 2021; 21:351-357. [PMID: 33601964 PMCID: PMC8086403 DOI: 10.1089/vbz.2020.2711] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Francisella tularensis is the etiologic agent of tularemia, a bacterial zoonotic disease. The genome of F. tularensis shows a recent evolutionary change, especially in reservoirs. Variable number of tandem repeats (VNTR) is described as a high-speed molecular clock and can thus be used as a high-resolution typing system. The main objective of our study was to investigate the molecular diversity of F. tularensis strains and reveal possible sources of infection. Using real-time PCR targeting the ISFtu2 region, we successfully amplified targeted DNA in 13/31 Slovenian patients with a clinical diagnosis of tularemia, and with PCR targeting the fopA gene, we obtained 11/13 PCR products. Sequencing revealed that all samples were identified as F. tularensis subsp. holarctica. We successfully obtained one F. tularensis isolate from a lymph node aspirate by culture on chocolate agar. Our isolate was clustered into major clade B12 (subclade B43). We optimized VNTR typing to be used directly on clinical samples. Multiple-locus VNTR analysis (MLVA) revealed five unique MLVA types; 45.5% samples had the same MLVA type, another 27.3% shared a different MLVA type, and each of the remaining had a unique MLVA type. Most samples differed at only two VNTR markers (Ft-M03 and Ft-M06). Additionally, we investigated samples from small mammals (n = 532) and Ixodes ricinus ticks (n = 232) captured in the same geographical area in which patients with tularemia were found. No F. tularensis DNA was detected in samples of small mammals or I. ricinus ticks. The diversity of MLVA types in Slovenia was high, despite the small region, but most of the samples from the same region shared the same MLVA type. Our results suggest that MLVA is a useful tool for quick molecular characterization of F. tularensis directly from patient samples, especially when investigating geographically localized outbreaks.
Collapse
Affiliation(s)
- Urška Glinšek Biškup
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Rok Kogoj
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Miša Korva
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Nataša Knap
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Tjaša Cerar
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Tea Knapič
- Slovenian Museum of Natural History, Ljubljana, Slovenia
| | - Miroslav Petrovec
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| | - Tatjana Avšič-Županc
- Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia
| |
Collapse
|
8
|
Kevin M, Girault G, Caspar Y, Cherfa MA, Mendy C, Maurin M, Ponsart C, Madani N. High-resolution melting PCR assay as a powerful tool for the epidemiological surveillance of tularemia in Western Europe. INFECTION GENETICS AND EVOLUTION 2021; 90:104741. [PMID: 33556556 DOI: 10.1016/j.meegid.2021.104741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Revised: 12/14/2020] [Accepted: 01/26/2021] [Indexed: 11/19/2022]
Abstract
In Europe, tularemia is caused by Francisella tularensis subsp. holarctica and is a sporadic disease affecting mainly wildlife animals and humans. Classification of this species relies on canonical single nucleotide polymorphisms (canSNPs). Four main clades have been described for F. tularensis subsp. holarctica: B.4, B.6, B.12 and B.16. Phylogeographic studies have shown that clade B.6 is predominant in Western Europe and B.12 in Eastern and Central Europe. Based on this global phylogeny, we aimed to design a molecular typing assay for all genetic subclades of subclade B.11, which is the predominant subclade in clade B.6. We designed high-resolution melting (HRM) primers for the screening of 109 canSNPs divided in seven orders of discrimination for the molecular epidemiology analysis and tracking of Francisella tularensis subsp. holarctica in Western Europe.
Collapse
Affiliation(s)
- Maëllys Kevin
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France
| | - Guillaume Girault
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France
| | - Yvan Caspar
- Laboratoire de Bactériologie-Hygiène hospitalière, Centre National de Référence des Francisella, CHU-Grenoble-Alpes, Grenoble, France
| | - Moulay Ali Cherfa
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France
| | - Christiane Mendy
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France
| | - Max Maurin
- Laboratoire de Bactériologie-Hygiène hospitalière, Centre National de Référence des Francisella, CHU-Grenoble-Alpes, Grenoble, France
| | - Claire Ponsart
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France
| | - Nora Madani
- Paris-Est University/ANSES, Animal Health Laboratory, Maisons-Alfort, France.
| |
Collapse
|
9
|
Linde J, Homeier-Bachmann T, Dangel A, Riehm JM, Sundell D, Öhrman C, Forsman M, Tomaso H. Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany. Microorganisms 2020; 8:microorganisms8121932. [PMID: 33291395 PMCID: PMC7761992 DOI: 10.3390/microorganisms8121932] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Revised: 11/26/2020] [Accepted: 12/02/2020] [Indexed: 02/04/2023] Open
Abstract
Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates-obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)-were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.
Collapse
Affiliation(s)
- Jörg Linde
- Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, 07743 Jena, Germany;
- Correspondence: ; Tel.: +49-3641-804-2320
| | - Timo Homeier-Bachmann
- Institute of Epidemiology, Friedrich-Loeffler-Institut, 17493 Greifswald-Insel Riems, Germany;
| | - Alexandra Dangel
- Bayerisches Landesamt für Gesundheit und Lebensmittelsicherheit, 85764 Oberschleißheim, Germany; (A.D.); (J.M.R.)
| | - Julia M. Riehm
- Bayerisches Landesamt für Gesundheit und Lebensmittelsicherheit, 85764 Oberschleißheim, Germany; (A.D.); (J.M.R.)
| | - David Sundell
- CBRN Defence and Security, Swedish Defence Research Agency (FOI), SE-901 82 Umeå, Sweden; (D.S.); (C.Ö.); (M.F.)
| | - Caroline Öhrman
- CBRN Defence and Security, Swedish Defence Research Agency (FOI), SE-901 82 Umeå, Sweden; (D.S.); (C.Ö.); (M.F.)
| | - Mats Forsman
- CBRN Defence and Security, Swedish Defence Research Agency (FOI), SE-901 82 Umeå, Sweden; (D.S.); (C.Ö.); (M.F.)
| | - Herbert Tomaso
- Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, 07743 Jena, Germany;
| |
Collapse
|
10
|
Busch A, Homeier-Bachmann T, Abdel-Glil MY, Hackbart A, Hotzel H, Tomaso H. Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis. PLoS Negl Trop Dis 2020; 14:e0008018. [PMID: 32991594 PMCID: PMC7523947 DOI: 10.1371/journal.pntd.0008018] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Accepted: 12/27/2019] [Indexed: 12/31/2022] Open
Abstract
By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.
Collapse
Affiliation(s)
- Anne Busch
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany
- * E-mail:
| | - Timo Homeier-Bachmann
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Epidemiology, Friedrich-Loeffler-Institut, Greifswald-Insel Riems, Germany
| | - Mostafa Y. Abdel-Glil
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany
| | - Anja Hackbart
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany
| | - Helmut Hotzel
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany
| | - Herbert Tomaso
- Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany
| |
Collapse
|
11
|
Mihelčić M, Marečić V, Ožanič M, Kelava I, Knežević M, Šantić M. Epidemiologic and Epizootic Data of Tularemia in the Past and in the Recent History in Croatia. Microorganisms 2020; 8:microorganisms8050721. [PMID: 32408663 PMCID: PMC7284554 DOI: 10.3390/microorganisms8050721] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2020] [Revised: 05/08/2020] [Accepted: 05/09/2020] [Indexed: 11/20/2022] Open
Abstract
Tularemia is a zoonotic disease caused by Francisella tularensis. A large number of recent studies have provided an update on the disease characteristics and the distribution across Europe. In Croatia, most of the clinical cases, as well as the reports of the disease in animals, date from the 20th century. In that period, epidemic and epizootic research had given detailed information about endemic regions and their characteristics, including suspected animal hosts and vectors. The region along the middle course of the Sava River, called Middle Posavina, is described as an endemic region, i.e., a “natural focus” of tularemia, in Croatia. In the 21st century, cases of human tularemia are being reported sporadically, with ulceloglandular, oropharyngeal and typhoid forms of disease. A majority of the described cases are linked with the consumption of contaminated food or water. The disease outbreaks still occur in areas along the course of the river Sava and in northwest Croatia. In this review article, we have summarized epidemiologic and epizootic data of tularemia in the past and in recent Croatian history.
Collapse
|
12
|
Appelt S, Köppen K, Radonić A, Drechsel O, Jacob D, Grunow R, Heuner K. Genetic Diversity and Spatial Segregation of Francisella tularensis Subspecies holarctica in Germany. Front Cell Infect Microbiol 2019; 9:376. [PMID: 31781515 PMCID: PMC6851236 DOI: 10.3389/fcimb.2019.00376] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 10/17/2019] [Indexed: 01/28/2023] Open
Abstract
Francisella tularensis is an intracellular pleomorphic bacterium and the causative agent of tularemia, a zoonotic disease with a wide host range. Among the F. tularensis subspecies, especially F. tularensis subsp. holarctica is of clinical relevance for European countries. The study presented herein focuses namely on genetic diversity and spatial segregation of F. tularensis subsp. holarctica in Germany, as still limited information is available. The investigation is based on the analysis of 34 F. tularensis subsp. holarctica isolates and one draft genome from an outbreak strain. The isolates were cultured from sample material being that of primarily human patients (n = 25) and free-living animals (n = 9). For six of 25 human isolates, epidemiological links between disease onset and tick bites could be established, confirming the importance of arthropod linked transmission of tularemia in Germany. The strains were assigned to three of four major F. tularensis subsp. holarctica clades: B.4, B.6, and B.12. Thereby, B.6 and B.12 clade members were predominantly found; only one human isolate was assigned to clade B.4. Also, it turned out that eight isolates which caused pneumonia in patients clustered into the B.6 clade. Altogether, eight different final subclades were assigned to clade B.6 (biovar I, erythromycin sensitive) and six to B.12 (biovar II, erythromycin resistant) in addition to one new final B.12 subclade. Moreover, for 13 human and 3 animal isolates, final subclade subdivisions were not assigned (B.12 subdivisions B.33 and B.34, and B.6 subdivision B.45) because official nomenclatures are not available yet. This gives credit to the genetic variability of F. tularensis subsp. holarctica strains in Germany. The results clearly point out that the given genetic diversity in Germany seems to be comparably high to that found in other European countries including Scandinavian regions. A spatial segregation of B.6 and B.12 strains was found and statistically confirmed, and B.12 clade members were predominantly found in eastern parts and B.6 members more in western to southern parts of Germany. The portion of B.12 clade members in northeastern parts of Germany was 78.5% and in southwestern parts 1.9%.
Collapse
Affiliation(s)
- Sandra Appelt
- Centre for Biological Threats and Special Pathogens (ZBS2), Robert Koch Institute, Berlin, Germany
| | - Kristin Köppen
- Working Group Cellular Interactions of Bacterial Pathogens, ZBS2, Robert Koch Institute, Berlin, Germany
| | - Aleksandar Radonić
- Methodology and Research Infrastructure Genome Sequencing (MF2), Robert Koch Institute, Berlin, Germany
| | - Oliver Drechsel
- Bioinformatics (MF1), Robert Koch Institute, Berlin, Germany
| | - Daniela Jacob
- Centre for Biological Threats and Special Pathogens (ZBS2), Robert Koch Institute, Berlin, Germany
| | - Roland Grunow
- Centre for Biological Threats and Special Pathogens (ZBS2), Robert Koch Institute, Berlin, Germany
| | - Klaus Heuner
- Working Group Cellular Interactions of Bacterial Pathogens, ZBS2, Robert Koch Institute, Berlin, Germany
| |
Collapse
|