Will the emergence of core genome MLST end the role of in silico MLST?

Corynebacterium pseudotuberculosis: Whole genome sequencing reveals unforeseen and relevant genetic diversity in this pathogen

Corynebacterium pseudotuberculosis (CPS) is an important bacterial animal pathogen. CPS causes chronic, debilitating and currently incurable infectious diseases affecting a wide range of livestock and wild herbivores including camelids worldwide. Belonging to the Corynebacterium diphtheriae complex, this pathogen can also infect humans. The classical characterization of CPS is typically based on the testing of nitrate reductase activity, separating the two biovars Equi and Ovis. However, more refined resolutions are required to unravel routes of infection. This was realized in our study by generating and analyzing whole genome sequencing (WGS) data. Using newly created core genome multilocus sequence typing (cgMLST) profiles we were the first to discover isolates grouping in a cluster adjacent to clusters formed by CPS biovar Equi isolates. This novel cluster includes CPS isolates from alpacas, llamas, camels and dromedaries, which are characterized by a lack of nitrate reductase activity as encountered in biovar Ovis. This is of special interest for molecular epidemiology. Nevertheless, these isolates bear the genes of the nitrate locus, which are characteristic of biovar Equi isolates. However, sequence analysis of the genes narG and narH of the nitrate locus revealed indels leading to frameshifts and inactivity of the enzymes involved in nitrate reduction. Interestingly, one CPS isolate originating from another lama with an insertion in the MFS transporter (narT) is adjacent to a cluster formed by ovine CPS isolates biovar Equi. Based on this knowledge, the combination of biochemical and PCR based molecular biological nitrate reductase detection can be used for a fast and uncomplicated classification of isolates in routine diagnostics in order to check the origin of camelid CPS isolates. Further analysis revealed that partial sequencing of the ABC transporter substrate binding protein (CP258_RS07935) is a powerful tool to assign the biovars and the novel genomovar.

Antimicrobial resistance and genomic investigation of Salmonella isolated from retail foods in Guizhou, China

Introduction

Salmonella is a major foodborne pathogen worldwide that causes severe morbidity and mortality. It is mainly caused by consuming contaminated food, with retail food considered the primary source.

Methods

In Guizhou, China, 102 Salmonella strains isolated from 2016 to 2021 underwent phenotypic antimicrobial resistance testing and whole-genome sequencing (WGS) to understand Salmonella diversity, including serotypes, sequencing types (STs), antimicrobial genes, virulence genes, plasmid types, multi-locus sequence types (MLST), and core genome MLST (cgMLST).

Results and discussion

S.Typhimurium was the dominant serotype, and O:4(B) was the leading serogroup. The most prevalent genotype was ST40. Phenotypic antimicrobial resistance identified 66.7% of the sampled isolates as multi-drug resistant (MDR). S.Enteritidis (n = 7), S.Typhimurium (n = 1), S.Indiana (n = 1), S.Kentucky (n = 1), S.Uganda (n = 1), all of which were MDR, were resistant to Colistin. Resistance rates varied significantly across different strains and food types, particularly meat products exhibiting higher resistance. Notably, significant increases in resistance were observed from 2016 to 2021 for the following: ≥ 1 resistant (P = 0.001), MDR (P = 0.001), ampicillin (P = 0.001), tetracycline (P < 0.001), chloramphenicol (P = 0.030), and trimethoprim/sulfamethoxazole (P = 0.003). The marked escalation in drug resistance over the recent years, coupled with the varying resistance rates among food sources, underscores the growing public health concern. Our findings highlight the need for a coordinated approach to effectively monitor and respond to Salmonella infections in Guizhou, China.

Emergence and Genomic Features of a mcr-1 Escherichia coli from Duck in Hungary

Plasmids carrying high-risk resistance mechanisms in pathogenic E. coli have gained particular attention in veterinary medicine, especially since the discovery of the colistin resistance gene, mcr-1. Here, we provide the first evidence of its emergence and describe the complete mcr-1 plasmid sequence of a multi-resistant avian pathogenic E. coli (APEC) strain from waterfowl in Hungary. Whole-genome sequencing analysis and core-genome MLST were performed to characterize the genome structure of the mcr-1 plasmid and to reveal the phylogenetic relation between the Hungarian duck strain Ec45-2020 and the internationally circulating mcr-1-positive E. coli strains from poultry and humans. Results showed that plasmid pEc45-2020-33kb displayed a high level of genome identity with mcr-1 plasmids of IncX4 type widespread among human, animal and food reservoirs of enteric bacteria of public health. The mcr-1-positive E. coli strain Ec45-2020 belongs to the ST162 genotype, considered as one of the globally disseminated zoonotic genotypes of MDR E. coli. In accordance with international findings, our results underline the importance of continuous surveillance of enteric bacteria with high-risk antimicrobial resistance genotypes, including neglected animals, such as waterfowls, as possible reservoirs for the colistin resistance gene mcr-1.

Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria

Core genome multilocus sequence typing (cgMLST) employs a strategy where the set of orthologous genes common to all members of a group of organisms are used for phylogenetic analysis of the group members. The Bacillus cereus group consists of species with pathogenicity towards insect species as well as warm-blooded animals including humans. While B. cereus is an opportunistic pathogen linked to a range of human disease conditions, including emesis and diarrhoea, Bacillus thuringiensis is an entomopathogenic species with toxicity toward insect larvae, and therefore used as a biological pesticide worldwide. Bacillus anthracis is a classical obligate pathogen causing anthrax, an acute lethal condition in herbivores as well as humans, and which is endemic in many parts of the world. The group also includes a range of additional species, and B. cereus group bacteria have been subject to analysis with a wide variety of phylogenetic typing systems. Here we present, based on analyses of 173 complete genomes from B. cereus group species available in public databases, the identification of a set of 1568 core genes which were used to create a core genome multilocus typing scheme for the group which is implemented in the PubMLST system as an open online database freely available to the community. The new cgMLST system provides unprecedented resolution over existing phylogenetic analysis schemes covering the B. cereus group.

Genomes-based MLST, cgMLST, wgMLST and SNP analysis of Salmonella Typhimurium from animals and humans

Salmonella Typhimurium (S. Typhimurium) is an important food-borne and zoonotic pathogen that causes salmonellosis. With the development of whole genome sequencing (WGS), genome-based typing has been widely applied to bacteriology. In this study, we investigated genotyping and phylogenetic clusters of S. Typhimurium isolates from humans and animals in different provinces (including Beijing, Shandong, Guangxi, Shaanxi, Henan, and Shanghai) of China during 2009-2018 using multi locus sequence typing (MLST), core genome MLST (cgMLST), whole genome MLST (wgMLST) and single nucleotide polymorphism (SNP) based on WGS. 29 S. Typhimurium isolates from chicken (n = 22), sick pigeon (n = 2), patients (n = 4) and sick swine (n = 1) were tested. MLST analysis showed S. Typhimurium strains were divided into four STs, namely ST19 (n = 14), ST34 (n = 12), ST128 (n = 2) and ST1544 (n = 1). cgMLST and wgMLST divided 29 strains into 27 cgSTs and 29 wgST, respectively. Phylogenetic clustering showed that all isolates were divided into 4 clusters and 4 singletons. SNP analysis was used to examine MLST, cgMLST, wgMLST analysis. Finally, comparisons of MLST, cgMLST, wgMLST, and SNP were analyzed and the results showed their precision increased in order. In summary, genomic typing and phylogenetic relationships of 29 S. Typhimurium strains from different sources in China were analyzed. These findings were beneficial to investigate molecular pathogenesis, bacterial diversity, and traceability analysis of Salmonella.

Application of a core genome sequence typing (cgMLST) pipeline for surveillance of Clostridioides difficile in China

Introduction

Clostridioides difficile (C. difficile) is a nosocomial bacterial pathogen that causes antibiotic-associated diarrhea mediated by cellular exotoxins secreted into the intestine during bacterial growth. Multilocus sequence typing (MLST) and PCR ribotyping are the main molecular typing for C. difficile. Whole genome sequencing (WGS) core genome multilocus sequence typing (cgMLST) was developed for genetic evolution and outbreak investigation of C. difficile with higher precision and accuracy.

Methods

A total of 699 whole (complete and draft) genome sequences of distinct C. difficile strains were used in this study to identify core gene set (2469 core genes) and the cgMLST scheme for the phylogeny analysis of C. difficile. This cgMLST pipeline was then carried the Chinese Pathogen Identification Net (China PIN) for surveillance of C. difficile in China. Within the China PIN, 195 WGS of C. difficile and an outbreak of CDI with 12 WGS of C. difficile were used to evaluate the cgMLST pipeline.

Results

The result displayed that mostly tested C. difficile isolates could be successfully divided into 5 classic clades and the outbreak event was also successfully identified.

Discussion

The results are meaningful and offer a practicable pipeline for a national-wide surveillance of C. difficile in China.

FastMLST: A Multi-core Tool for Multilocus Sequence Typing of Draft Genome Assemblies

Multilocus Sequence Typing (MLST) is a precise microbial typing approach at the intra-species level for epidemiologic and evolutionary purposes. It operates by assigning a sequence type (ST) identifier to each specimen, based on a combination of alleles of multiple housekeeping genes included in a defined scheme. The use of MLST has multiplied due to the availability of large numbers of genomic sequences and epidemiologic data in public repositories. However, data processing speed has become problematic due to the massive size of modern datasets. Here, we present FastMLST, a tool that is designed to perform PubMLST searches using BLASTn and a divide-and-conquer approach that processes each genome assembly in parallel. The output offered by FastMLST includes a table with the ST, allelic profile, and clonal complex or clade (when available), detected for a query, as well as a multi-FASTA file or a series of FASTA files with the concatenated or single allele sequences detected, respectively. FastMLST was validated with 91 different species, with a wide range of guanine-cytosine content (%GC), genome sizes, and fragmentation levels, and a speed test was performed on 3 datasets with varying genome sizes. Compared with other tools such as mlst, CGE/MLST, MLSTar, and PubMLST, FastMLST takes advantage of multiple processors to simultaneously type up to 28 000 genomes in less than 10 minutes, reducing processing times by at least 3-fold with 100% concordance to PubMLST, if contaminated genomes are excluded from the analysis. The source code, installation instructions, and documentation of FastMLST are available at https://github.com/EnzoAndree/FastMLST

Serotyping, MLST, and Core Genome MLST Analysis of Salmonella enterica From Different Sources in China During 2004-2019

Salmonella enterica (S. enterica) is an important foodborne pathogen, causing food poisoning and human infection, and critically threatening food safety and public health. Salmonella typing is essential for bacterial identification, tracing, epidemiological investigation, and monitoring. Serotyping and multilocus sequence typing (MLST) analysis are standard bacterial typing methods despite the low resolution. Core genome MLST (cgMLST) is a high-resolution molecular typing method based on whole genomic sequencing for accurate bacterial tracing. We investigated 250 S. enterica isolates from poultry, livestock, food, and human sources in nine provinces of China from 2004 to 2019 using serotyping, MLST, and cgMLST analysis. All S. enterica isolates were divided into 36 serovars using slide agglutination. The major serovars in order were Enteritidis (31 isolates), Typhimurium (29 isolates), Mbandaka (23 isolates), and Indiana (22 isolates). All strains were assigned into 43 sequence types (STs) by MLST. Among them, ST11 (31 isolates) was the primary ST. Besides this, a novel ST, ST8016, was identified, and it was different from ST40 by position 317 C → T in dnaN. Furthermore, these 250 isolates were grouped into 185 cgMLST sequence types (cgSTs) by cgMLST. The major cgST was cgST235530 (11 isolates), and only three cgSTs contained isolates from human and other sources, indicating a possibility of cross-species infection. Phylogenetic analysis indicated that most of the same serovar strains were putatively homologous except Saintpaul and Derby due to their multilineage characteristics. In addition, serovar I 4,[5],12:i:- and Typhimurium isolates have similar genomic relatedness on the phylogenetic tree. In conclusion, we sorted out the phenotyping and genotyping diversity of S. enterica isolates in China during 2004-2019 and clarified the temporal and spatial distribution characteristics of Salmonella from different hosts in China in the recent 16 years. These results greatly supplement Salmonella strain resources, genetic information, and traceability typing data; facilitate the typing, traceability, identification, and genetic evolution analysis of Salmonella; and therefore, improve the level of analysis, monitoring, and controlling of foodborne microorganisms in China.

Historical, current, and emerging tools for identification and serotyping of Shigella

The Shigella genus includes serious foodborne disease etiologic agents, with 4 species and 54 serotypes. Identification at species and serotype levels is a crucial task in microbiological laboratories. Nevertheless, the genetic similarity between Shigella spp. and Escherichia coli challenges the correct identification and serotyping of Shigella spp., with subsequent negative repercussions on surveillance, epidemiological investigations, and selection of appropriate treatments. For this purpose, multiple techniques have been developed historically ranging from phenotype-based methods and single or multilocus molecular techniques to whole-genome sequencing (WGS). To facilitate the selection of the most relevant method, we herein provide a global overview of historical and emerging identification and serotyping techniques with a particular focus on the WGS-based approaches. This review highlights the excellent discriminatory power of WGS to more accurately elucidate the epidemiology of Shigella spp., disclose novel promising genomic targets for surveillance methods, and validate previous well-established methods.

Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan

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.

Multilocus Sequence Typing of Staphylococcus aureus

Multilocus sequence typing (MLST) has been successfully used to differentiate and trace the bacterial species and pathogens that cause outbreaks or epidemics of infectious diseases. MLST provides a powerful solution for molecular epidemiological characterization of bacterial strains, including Staphylococcus aureus, by using the sequences of the internal region of seven housekeeping genes. In previous studies, we utilized MLST to analyze the genotypes of S. aureus isolates from pediatric patients with cystic fibrosis and revealed three prevalent ST types, including ST5, ST30, and ST8 in these isolates. In this chapter, we describe a detailed procedure of MLST for genotyping S. aureus.

Molecular epidemiology in current times

Motivated to find options for prevention or intervention, molecular epidemiology aims to identify the host and microbial factors that determine the transmission, manifestation and progression of infectious disease. The genotyping of cultivatable bacterial strains is performed by either anonymous fingerprinting techniques or sequence-based exploration of variable genomic sites. Multilocus sequence typing of housekeeping genes and allele profiling of the core genome have become standard techniques of bacterial strain typing that may be supplemented by whole genome sequencing to explore all single nucleotide variants and/or the composition of the accessory genome. Next, novel protocols to investigate host and microbiome based upon smart third generation sequencing technologies are being developed for an effective surveillance, rapid diagnosis and real-time tracking of infectious diseases.

A Core Genome Multilocus Sequence Typing Scheme for Streptococcus mutans

Streptococcus mutans is one of the primary pathogens responsible for the development of dental caries. Recent whole-genome sequencing (WGS)-based core genome multilocus sequence typing (cgMLST) approaches have been employed in epidemiological studies of specific human pathogens. However, this approach has not been reported in studies of S. mutans Here, we therefore developed a cgMLST scheme for S. mutans We surveyed 199 available S. mutans genomes as a means of identifying cgMLST targets, developing a scheme that incorporated 594 targets from the S. mutans UA159 reference genome. Sixty-eight sequence types (STs) were identified in this cgMLST scheme (cgSTs) in 80 S. mutans isolates from 40 children that were sequenced in this study, compared to 35 STs identified by multilocus sequence typing (MLST). Fifty-six cgSTs (82.35%) were associated with a single isolate based on our cgMLST scheme, which is significantly higher than in the MLST scheme (11.43%). In addition, 58.06% of all MLST profiles with ≥2 isolates were further differentiated by our cgMLST scheme. Topological analyses of the maximum likelihood phylogenetic trees revealed that our cgMLST scheme was more reliable than the MLST scheme. A minimum spanning tree of 145 S. mutans isolates from 10 countries developed based upon the cgMLST scheme highlighted the diverse population structure of S. mutans This cgMLST scheme thus offers a new molecular typing method suitable for evaluating the epidemiological distribution of this pathogen and has the potential to serve as a benchmark for future global studies of the epidemiological nature of dental caries.IMPORTANCE Streptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries.

A Core Genome Multilocus Sequence Typing Scheme for Pseudomonas aeruginosa

Pseudomonas aeruginosa is a ubiquitous microorganism and an important opportunistic pathogen responsible for a broad spectrum of infections mainly in immunosuppressed and critically ill patients. Molecular investigations traditionally rely on pulsed field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). In this work we propose a core genome multilocus sequence typing (cgMLST) scheme for P. aeruginosa, a methodology that combines traditional MLST principles with whole genome sequencing data. All publicly available complete P. aeruginosa genomes, representing the diversity of this species, were used to establish a cgMLST scheme targeting 2,653 genes. The scheme was then tested using genomes available at contig, chromosome and scaffold levels. The proposed cgMLST scheme for P. aeruginosa typed over 99% (2,314/2,325) of the genomes available for this study considering at least 95% of the cgMLST target genes present. The absence of a certain number gene targets at the threshold considered for both the creation and validation steps due to low genome sequence quality is possibly the main reason for this result. The cgMLST scheme was compared with previously published whole genome single nucleotide polymorphism analysis for the characterization of the population structure of the epidemic clone ST235 and results were highly similar. In order to evaluate the typing resolution of the proposed scheme, collections of isolates belonging to two important STs associated with cystic fibrosis, ST146 and ST274, were typed using this scheme, and ST235 isolates associated with an outbreak were evaluated. Besides confirming the relatedness of all the isolates, earlier determined by MLST, the higher resolution of cgMLST denotes that it may be suitable for surveillance programs, overcoming possible shortcomings of classical MLST. The proposed scheme is publicly available at: https://github.com/BioinformaticsHIAEMolecularMicrobiology/cgMLST-Pseudomonas-aeruginosa.

Clostridium difficile ribotype 017 - characterization, evolution and epidemiology of the dominant strain in Asia

Clostridium difficile ribotype (RT) 017 is an important toxigenic C. difficile RT which, due to a deletion in the repetitive region of the tcdA gene, only produces functional toxin B. Strains belonging to this RT were initially dismissed as nonpathogenic and circulated largely undetected for almost two decades until they rose to prominence following a series of outbreaks in the early 2000s. Despite lacking a functional toxin A, C. difficile RT 017 strains have been shown subsequently to be capable of causing disease as severe as that caused by strains producing both toxins A and B. While C. difficile RT 017 strains can be found in almost every continent today, epidemiological studies suggest that the RT is endemic in Asia and that the global spread of this MLST clade 4 lineage member is a relatively recent event. C. difficile RT 017 transmission appears to be mostly from human to human with only a handful of reports of isolations from animals. An important feature of C. difficile RT 017 strains is their resistance to several antimicrobials and this has been documented as a possible factor driving multiple outbreaks in different parts of the world. This review summarizes what is currently known regarding the emergence and evolution of strains belonging to C. difficile RT 017 as well as features that have allowed it to become an RT of global importance.

Fourier-Transform InfraRed Spectroscopy Can Quickly Type Gram-Negative Bacilli Responsible for Hospital Outbreaks

The typing of epidemic bacterial pathogens in hospitals relies on DNA-based, expensive, and time-consuming techniques, that are often limited to retrospective studies. However, the quick identification of epidemic pathogens in the routine of the microbiology laboratories would expedite infection control procedures that limit the contamination of new patients. IR Biotyper (Bruker Daltonics GmbH) is a new typing machine based on Fourier-transform infrared (FTIR) spectroscopy which generates spectra, aiming at typing the micro-organisms within 3 h. This technique discriminates the isolates by exploring the differences of the surface cell polysaccharides. In this work, we evaluated the ability of the FTIR spectroscopy to recognize Gram-negative bacilli clones responsible for hospital outbreaks. Isolates of Pseudomonas aeruginosa (n = 100), Klebsiella pneumoniae (n = 16), Enterobacter cloacae (n = 23), and Acinetobacter baumannii (n = 20) were typed by the reference methods Multi-Locus Sequence Typing (defining sequence types – STs) along with or without pulsed field gel electrophoresis (PFGE) (defining pulsotypes), and by FTIR spectroscopy. The congruence of FTIR spectroscopy clustering was compared to those of MLST and PFGE by Adjusted Rand index and Adjusted Wallace coefficient. We found that FTIR spectroscopy accurately clustered P. aeruginosa, K. pneumoniae, and E. cloacae isolates belonging to the same ST. The performance of the FTIR spectroscopy was slightly lower for A. baumannii. Furthermore, FTIR spectroscopy also correctly clustered P. aeruginosa isolates having a similar pulsotype. Overall, the IR Biotyper can quickly (in less than 3 h) detect the spread of clones of P. aeruginosa, K. pneumoniae, E. cloacae, and A. baumannii. The use of this technique by clinical microbiology laboratories may help to tackle the spread of epidemic clones by the quick implementation of infection control measures.