scholarly journals Population Structure of Mycobacterium bovis in Germany: a Long-Term Study Using Whole-Genome Sequencing Combined with Conventional Molecular Typing Methods

2020 ◽  
Vol 58 (11) ◽  
Author(s):  
Thomas A. Kohl ◽  
Katharina Kranzer ◽  
Sönke Andres ◽  
Thierry Wirth ◽  
Stefan Niemann ◽  
...  

ABSTRACT Mycobacterium bovis is the primary cause of bovine tuberculosis (bTB) and infects a wide range of domestic animal and wildlife species and humans. In Germany, bTB still emerges sporadically in cattle herds, free-ranging wildlife, diverse captive animal species, and humans. In order to understand the underlying population structure and estimate the population size fluctuation through time, we analyzed 131 M. bovis strains from animals (n = 38) and humans (n = 93) in Germany from 1999 to 2017 by whole-genome sequencing (WGS), mycobacterial interspersed repetitive-unit–variable-number tandem-repeat (MIRU-VNTR) typing, and spoligotyping. Based on WGS data analysis, 122 out of the 131 M. bovis strains were classified into 13 major clades, of which 6 contained strains from both human and animal cases and 7 only strains from human cases. Bayesian analyses suggest that the M. bovis population went through two sharp anticlimaxes, one in the middle of the 18th century and another one in the 1950s. WGS-based cluster analysis grouped 46 strains into 13 clusters ranging in size from 2 to 11 members and involving strains from distinct host types, e.g., only cattle and also mixed hosts. Animal strains of four clusters were obtained over a 9-year span, pointing toward autochthonous persistent bTB infection cycles. As expected, WGS had a higher discriminatory power than spoligotyping and MIRU-VNTR typing. In conclusion, our data confirm that WGS and suitable bioinformatics constitute the method of choice to implement prospective molecular epidemiological surveillance of M. bovis. The population of M. bovis in Germany is diverse, with subtle, but existing, interactions between different host groups.

2021 ◽  
Vol 9 (8) ◽  
pp. 1585
Author(s):  
Ana C. Reis ◽  
Liliana C. M. Salvador ◽  
Suelee Robbe-Austerman ◽  
Rogério Tenreiro ◽  
Ana Botelho ◽  
...  

Classical molecular analyses of Mycobacterium bovis based on spoligotyping and Variable Number Tandem Repeat (MIRU-VNTR) brought the first insights into the epidemiology of animal tuberculosis (TB) in Portugal, showing high genotypic diversity of circulating strains that mostly cluster within the European 2 clonal complex. Previous surveillance provided valuable information on the prevalence and spatial occurrence of TB and highlighted prevalent genotypes in areas where livestock and wild ungulates are sympatric. However, links at the wildlife–livestock interfaces were established mainly via classical genotype associations. Here, we apply whole genome sequencing (WGS) to cattle, red deer and wild boar isolates to reconstruct the M. bovis population structure in a multi-host, multi-region disease system and to explore links at a fine genomic scale between M. bovis from wildlife hosts and cattle. Whole genome sequences of 44 representative M. bovis isolates, obtained between 2003 and 2015 from three TB hotspots, were compared through single nucleotide polymorphism (SNP) variant calling analyses. Consistent with previous results combining classical genotyping with Bayesian population admixture modelling, SNP-based phylogenies support the branching of this M. bovis population into five genetic clades, three with apparent geographic specificities, as well as the establishment of an SNP catalogue specific to each clade, which may be explored in the future as phylogenetic markers. The core genome alignment of SNPs was integrated within a spatiotemporal metadata framework to further structure this M. bovis population by host species and TB hotspots, providing a baseline for network analyses in different epidemiological and disease control contexts. WGS of M. bovis isolates from Portugal is reported for the first time in this pilot study, refining the spatiotemporal context of TB at the wildlife–livestock interface and providing further support to the key role of red deer and wild boar on disease maintenance. The SNP diversity observed within this dataset supports the natural circulation of M. bovis for a long time period, as well as multiple introduction events of the pathogen in this Iberian multi-host system.


2021 ◽  
Vol 12 ◽  
Author(s):  
Ana Pelerito ◽  
Alexandra Nunes ◽  
Teresa Grilo ◽  
Joana Isidro ◽  
Catarina Silva ◽  
...  

Brucellosis is an important zoonosis that is emerging in some regions of the world, gaining increased relevance with the inclusion of the causing agent Brucella spp. in the class B bioterrorism group. Until now, multi-locus VNTR Analysis (MLVA) based on 16 loci has been considered as the gold standard for Brucella typing. However, this methodology is laborious, and, with the rampant release of Brucella genomes, the transition from the traditional MLVA to whole genome sequencing (WGS)-based typing is on course. Nevertheless, in order to avoid a disruptive transition with the loss of massive genetic data obtained throughout the last decade and considering that the transition timings will vary considerably among different countries, it is important to determine WGS-based MLVA alleles of the nowadays sequenced genomes. On this regard, we aimed to evaluate the performance of a Python script that had been previously developed for the rapid in silico extraction of the MLVA alleles, by comparing it to the PCR-based MLVA procedure over 83 strains from different Brucella species. The WGS-based MLVA approach detected 95.3% of all possible 1,328 hits (83 strains×16 loci) and showed an agreement rate with the PCR-based MLVA procedure of 96.4% for MLVA-16. According to our dataset, we suggest the use of a minimal depth of coverage of ~50x and a maximum number of ~200 contigs as guiding “boundaries” for the future application of the script. In conclusion, the evaluated script seems to be a very useful and robust tool for the in silico determination of MLVA profiles of Brucella strains, allowing retrospective and prospective molecular epidemiological studies, which are important for maintaining an active epidemiological surveillance of brucellosis.


2021 ◽  
Author(s):  
Francisco M Ortuno ◽  
Carlos Loucera ◽  
Carlos S Casimiro-Soriguer ◽  
Jose A Lepe ◽  
Pedro Camacho Martinez ◽  
...  

The current SARS-CoV-2 pandemic has emphasized the utility of viral whole genome sequencing in the surveillance and control of the pathogen. An unprecedented ongoing global initiative is increasingly producing hundreds of thousands of sequences worldwide. However, the complex circumstances in which viruses are sequenced, along with the demand of urgent results, causes a high rate of incomplete and therefore useless, sequences. However, viral sequences evolve in the context of a complex phylogeny and therefore different positions along the genome are in linkage disequilibrium. Therefore, an imputation method would be able to predict missing positions from the available sequencing data. We developed impuSARS, an application that includes Minimac, the most widely used strategy for genomic data imputation and, taking advantage of the enormous amount of SARS-CoV-2 whole genome sequences available, a reference panel containing 239,301 sequences was built. The impuSARS application was tested in a wide range of conditions (continuous fragments, amplicons or sparse individual positions missing) showing great fidelity when reconstructing the original sequences. The impuSARS application is also able to impute whole genomes from commercial kits covering less than 20% of the genome or only from the Spike protein with a precision of 0.96. It also recovers the lineage with a 100% precision for almost all the lineages, even in very poorly covered genomes (< 20%). Imputation can improve the pace of SARS-CoV-2 sequencing production by recovering many incomplete or low-quality sequences that would be otherwise discarded. impuSARS can be incorporated in any primary data processing pipeline for SARS-CoV-2 whole genome sequencing.


2018 ◽  
Vol 56 (9) ◽  
Author(s):  
Thomas A. Kohl ◽  
Christian Utpatel ◽  
Stefan Niemann ◽  
Irmgard Moser

ABSTRACTBovine tuberculosis (bTB) caused byMycobacterium bovisis a transmissible disease notifiable to the World Organization for Animal Health and to the European Union, with ongoing efforts of surveillance and eradication in every EU member state. In Germany, a country which has been declared officially free from bovine tuberculosis since 1997 by the EU,M. bovisinfections still occur sporadically in cattle and other mammals, including humans. Here, the transmission routes of a bTB outbreak in a wildlife park in Germany affecting different cervid species, bison, lynx, and pot-bellied pigs were followed by employing whole-genome sequencing (WGS) combined with spoligotyping and mycobacterial interspersed repetitive-unit–variable-number tandem-repeat (MIRU-VNTR) typing. One singleM. bovisstrain persisted from 2002 to 2015, and transmission between the park and a distantly located captive cervid farm was verified. The spoligotyping patterns remained identical, while MIRU-VNTR typing of 24 loci of the standardized panel and locus 2163a as an additional locus revealed one change at locus 2165 in a strain from a fallow deer and one at locus 2461 in isolates from red deer over the whole time period. WGS analysis confirmed the close relatedness of the isolates, with a maximum of 12 single nucleotide polymorphisms (SNPs) detected between any two sequenced isolates. In conclusion, our data confirm a longitudinal outbreak ofM. bovisin a German wildlife park and provide the first insights into the dynamics of different genotyping markers inM. bovis.


2015 ◽  
Vol 54 (2) ◽  
pp. 333-342 ◽  
Author(s):  
Jason C. Kwong ◽  
Karolina Mercoulia ◽  
Takehiro Tomita ◽  
Marion Easton ◽  
Hua Y. Li ◽  
...  

Whole-genome sequencing (WGS) has emerged as a powerful tool for comparing bacterial isolates in outbreak detection and investigation. Here we demonstrate that WGS performed prospectively for national epidemiologic surveillance ofListeria monocytogeneshas the capacity to be superior to our current approaches using pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), multilocus variable-number tandem-repeat analysis (MLVA), binary typing, and serotyping. Initially 423L. monocytogenesisolates underwent WGS, and comparisons uncovered a diverse genetic population structure derived from three distinct lineages. MLST, binary typing, and serotyping results inferredin silicofrom the WGS data were highly concordant (>99%) with laboratory typing performed in parallel. However, WGS was able to identify distinct nested clusters within groups of isolates that were otherwise indistinguishable using our current typing methods. Routine WGS was then used for prospective epidemiologic surveillance on a further 97L. monocytogenesisolates over a 12-month period, which provided a greater level of discrimination than that of conventional typing for inferring linkage to point source outbreaks. A risk-based alert system based on WGS similarity was used to inform epidemiologists required to act on the data. Our experience shows that WGS can be adopted for prospectiveL. monocytogenessurveillance and investigated for other pathogens relevant to public health.


2014 ◽  
Vol 53 (1) ◽  
pp. 212-218 ◽  
Author(s):  
Xiangyu Deng ◽  
Nikki Shariat ◽  
Elizabeth M. Driebe ◽  
Chandler C. Roe ◽  
Beth Tolar ◽  
...  

A retrospective investigation was performed to evaluate whole-genome sequencing as a benchmark for comparing molecular subtyping methods forSalmonella entericaserotype Enteritidis and survey the population structure of commonly encounteredS. entericaserotype Enteritidis outbreak isolates in the United States. A total of 52S. entericaserotype Enteritidis isolates representing 16 major outbreaks and three sporadic cases collected between 2001 and 2012 were sequenced and subjected to subtyping by four different methods: (i) whole-genome single-nucleotide-polymorphism typing (WGST), (ii) multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA), (iii) clustered regularly interspaced short palindromic repeats combined with multi-virulence-locus sequence typing (CRISPR-MVLST), and (iv) pulsed-field gel electrophoresis (PFGE). WGST resolved all outbreak clusters and provided useful robust phylogenetic inference results with high epidemiological correlation. While both MLVA and CRISPR-MVLST yielded higher discriminatory power than PFGE, MLVA outperformed the other methods in delineating outbreak clusters whereas CRISPR-MVLST showed the potential to trace major lineages and ecological origins ofS. entericaserotype Enteritidis. Our results suggested that whole-genome sequencing makes a viable platform for the evaluation and benchmarking of molecular subtyping methods.


2020 ◽  
Author(s):  
Miguel Pinto ◽  
Vítor Borges ◽  
Joana Isidro ◽  
João Carlos Rodrigues ◽  
Luís Vieira ◽  
...  

Neisseria gonorrhoeae , the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managing N. gonorrhoeae infections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to perform N. gonorrhoeae clustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspected N. gonorrhoeae genomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with available N. gonorrhoeae typing methods [i.e. MLST (multilocus sequence typing), NG-MAST ( N. gonorrhoeae multi-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis of N. gonorrhoeae genetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.


mSphere ◽  
2019 ◽  
Vol 4 (2) ◽  
Author(s):  
Yachen Hu ◽  
Zhenyu Wang ◽  
Bin Qiang ◽  
Yaohui Xu ◽  
Xiang Chen ◽  
...  

ABSTRACTSalmonella entericasubspeciesentericaserovar Gallinarum biovar Pullorum (S. Pullorum) is the etiological agent of pullorum disease, causing white diarrhea with high mortality in chickens. There are many unsolved issues surrounding the epidemiology ofS. Pullorum, including its origin and transmission history as well as the discordance between its phenotypic heterogeneity and genetic monomorphism. In this paper, we report the results of whole-genome sequencing of a panel of 97S. Pullorum strains isolated between 1962 and 2014 from four countries across three continents. We utilized 6,795 core genome single nucleotide polymorphisms (SNPs) to reconstruct a phylogenetic tree within a spatiotemporal Bayesian framework, estimating that the most recent common ancestor ofS. Pullorum emerged in ∼914 CE (95% confidence interval [95%CI], 565 to 1273 CE). The extantS. Pullorum strains can be divided into four distinct lineages, each of which is significantly associated with geographical distribution. The intercontinental transmissions of lineages III and IV can be traced to the mid-19th century and are probably related to the “Hen Fever” prevalent at that time. Further genomic analysis indicated that the loss or pseudogenization of functional genes involved in metabolism and virulence inS. Pullorum has been ongoing since before and after divergence from the ancestor. In contrast, multiple prophages and plasmids have been acquired byS. Pullorum, and these have endowed it with new characteristics, especially the multidrug resistance conferred by two large plasmids in lineage I. The results of this study provide insight into the evolution ofS. Pullorum and prove the efficiency of whole-genome sequencing in epidemiological surveillance of pullorum disease.IMPORTANCEPullorum disease, an acute poultry septicemia caused bySalmonellaGallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history ofS. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution ofS. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.


Author(s):  
M.I. Terekhova ◽  
◽  
E.V. Rogacheva ◽  
I.A. Derevyanchenko ◽  
L.A. Kraeva ◽  
...  

The increasing number of antibiotic-resistant isolates of L. monocytogenes is required to establish a genotypic resistance profile to ensure appropriate antibiotic therapy of listeriosis. In this study, whole-genome sequencing and de novo assembly was performed on L. monocytogenes strains from St. Petersburg and the Vologda region. We obtained the MLST ST, phylogenetic lineage and PCR-serogroups in silico for isolates under the study, revealed genes and mutations associated with antibiotic resistance. In general, the genetic composition was similar between the strains from different regions and included a wide range of antibiotic resistance mechanisms. Listeria strains possessed genes that code for resistance to β-lactam antibiotics, fluoroquinolones, tetracyclines and macrolides, — classes that are commonly used in the treatment of listeria infection. The present study is important in the sanitary and epidemiological surveillance of listeriosis in Russia.


2015 ◽  
Vol 53 (11) ◽  
pp. 3565-3573 ◽  
Author(s):  
Anne Holmes ◽  
Lesley Allison ◽  
Melissa Ward ◽  
Timothy J. Dallman ◽  
Richard Clark ◽  
...  

Detailed laboratory characterization ofEscherichia coliO157 is essential to inform epidemiological investigations. This study assessed the utility of whole-genome sequencing (WGS) for outbreak detection and epidemiological surveillance ofE. coliO157, and the data were used to identify discernible associations between genotypes and clinical outcomes. One hundred fiveE. coliO157 strains isolated over a 5-year period from human fecal samples in Lothian, Scotland, were sequenced with the Ion Torrent Personal Genome Machine. A total of 8,721 variable sites in the core genome were identified among the 105 isolates; 47% of the single nucleotide polymorphisms (SNPs) were attributable to six “atypical”E. coliO157 strains and included recombinant regions. Phylogenetic analyses showed that WGS correlated well with the epidemiological data. Epidemiological links existed between cases whose isolates differed by three or fewer SNPs. WGS also correlated well with multilocus variable-number tandem repeat analysis (MLVA) typing data, with only three discordant results observed, all among isolates from cases not known to be epidemiologically related. WGS produced a better-supported, higher-resolution phylogeny than MLVA, confirming that the method is more suitable for epidemiological surveillance ofE. coliO157. A combination ofinsilicoanalyses (VirulenceFinder, ResFinder, and local BLAST searches) were used to determinestxsubtypes, multilocus sequence types (15 loci), and the presence of virulence and acquired antimicrobial resistance genes. There was a high level of correlation between the WGS data and our routine typing methods, although some discordant results were observed, mostly related to the limitation of short sequence read assembly. The data were used to identify sublineages and clades ofE. coliO157, and when they were correlated with the clinical outcome data, they showed that one clade, Ic3, was significantly associated with severe disease. Together, the results show that WGS data can provide higher resolution of the relationships betweenE. coliO157 isolates than that provided by MLVA. The method has the potential to streamline the laboratory workflow and provide detailed information for the clinical management of patients and public health interventions.


Sign in / Sign up

Export Citation Format

Share Document