DODGE: automated point source bacterial outbreak detection using cumulative long term genomic surveillance.

SUMMARY: The reliable and timely recognition of outbreaks is a key component of public health surveillance for foodborne diseases. Whole genome sequencing (WGS) offers high resolution typing of foodborne bacterial pathogens and facilitates the accurate detection of outbreaks. This detection relies on grouping WGS data into clusters at an appropriate genetic threshold. However, methods and tools for selecting and adjusting such thresholds according to the required resolution of surveillance and epidemiological context are lacking. Here we present DODGE (Dynamic Outbreak Detection for Genomic Epidemiology), an algorithm to dynamically select and compare these genetic thresholds. DODGE can analyse expanding datasets over time and clusters that are predicted to correspond to outbreaks (or "investigation clusters") can be named with established genomic nomenclature systems to facilitate integrated analysis across jurisdictions. DODGE was tested in two real-world Salmonella genomic surveillance datasets of different duration, 2 months from Australia and 9 years from the United Kingdom. In both cases only a minority of isolates were identified as investigation clusters. Two known outbreaks in the United Kingdom dataset were detected by DODGE and were recognized at an earlier timepoint than the outbreaks were reported. These findings demonstrated the potential of the DODGE approach to improve the effectiveness and timeliness of genomic surveillance for foodborne diseases and the effectiveness of the algorithm developed. AVAILABILITY AND IMPLEMENTATION: DODGE is freely available at https://github.com/LanLab/dodge and can easily be installed using Conda.

Emergence of Poultry-Associated Human Salmonella enterica Serovar Abortusovis Infections, New South Wales, Australia.

Salmonella enterica serovar Abortusovis is a ovine-adapted pathogen that causes spontaneous abortion. Salmonella Abortusovis was reported in poultry in 2009 and has since been reported in human infections in New South Wales, Australia. Phylogenomic analysis revealed a clade of 51 closely related isolates from Australia originating in 2004. That clade was genetically distinct from ovine-associated isolates. The clade was widespread in New South Wales poultry production facilities but was only responsible for sporadic human infections. Some known virulence factors associated with human infections were only found in the poultry-associated clade, some of which were acquired through prophages and plasmids. Furthermore, the ovine-associated clade showed signs of genome decay, but the poultry-associated clade did not. Those genomic changes most likely led to differences in host range and disease type. Surveillance using the newly identified genetic markers will be vital for tracking Salmonella Abortusovis transmission in animals and to humans and preventing future outbreaks.

A Multilocus Sequence Typing Scheme for Rapid Identification of Xanthomonas citri Based on Whole-Genome Sequencing Data.

Xanthomonas citri is a plant-pathogenic bacterium associated with a diverse range of host plant species. It has undergone substantial reclassification and currently consists of 14 different subspecies or pathovars that are responsible for a wide range of plant diseases. Whole-genome sequencing (WGS) provides a cutting-edge advantage over other diagnostic techniques in epidemiological and evolutionary studies of X. citri because it has a higher discriminatory power and is replicable across laboratories. WGS also allows for the improvement of multilocus sequence typing (MLST) schemes. In this study, we used genome sequences of Xanthomonas isolates from the NCBI RefSeq database to develop a seven-gene MLST scheme that yielded 19 sequence types (STs) that correlated with phylogenetic clades of X. citri subspecies or pathovars. Using this MLST scheme, we examined 2,911 Xanthomonas species assemblies from NCBI GenBank and identified 15 novel STs from 37 isolates that were misclassified in NCBI. In total, we identified 545 X. citri assemblies from GenBank with 95% average nucleotide identity to the X. citri type strain, and all were classified as one of the 34 STs. All MLST classifications correlated with a phylogenetic position inferred from alignments using 92 conserved genes. We observed several instances where strains from different pathovars formed closely related monophyletic clades and shared the same ST, indicating that further investigation of the validity of these pathovars is required. Our MLST scheme described here is a robust tool for rapid classification of X. citri pathovars using WGS and a powerful method for further comprehensive taxonomic revision of X. citri pathovars.

SaLTy: a novel Staphylococcus aureus Lineage Typer.

Staphylococcus aureus asymptomatically colonises 30 % of humans but can also cause a range of diseases, which can be fatal. In 2017 S. aureus was associated with 20 000 deaths in the USA alone. Dividing S. aureus isolates into smaller sub-groups can reveal the emergence of distinct sub-populations with varying potential to cause infections. Despite multiple molecular typing methods categorising such sub-groups, they do not take full advantage of S. aureus genome sequences when describing the fundamental population structure of the species. In this study, we developed Staphylococcus aureus Lineage Typing (SaLTy), which rapidly divides the species into 61 phylogenetically congruent lineages. Alleles of three core genes were identified that uniquely define the 61 lineages and were used for SaLTy typing. SaLTy was validated on 5000 genomes and 99.12 % (4956/5000) of isolates were assigned the correct lineage. We compared SaLTy lineages to previously calculated clonal complexes (CCs) from BIGSdb (n=21 173). SALTy improves on CCs by grouping isolates congruently with phylogenetic structure. SaLTy lineages were further used to describe the carriage of Staphylococcal chromosomal cassette containing mecA (SCCmec) which is carried by methicillin-resistant S. aureus (MRSA). Most lineages had isolates lacking SCCmec and the four largest lineages varied in SCCmec over time. Classifying isolates into SaLTy lineages, which were further SCCmec typed, allowed SaLTy to describe high-level MRSA epidemiology. We provide SaLTy as a simple typing method that defines phylogenetic lineages (https://github.com/LanLab/SaLTy). SaLTy is highly accurate and can quickly analyse large amounts of S. aureus genome data. SaLTy will aid the characterisation of S. aureus populations and ongoing surveillance of sub-groups that threaten human health.

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses.

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

Prevotella copri alleviates hyperglycemia and regulates gut microbiota and metabolic profiles in mice.

Prevotella copri is the dominant species of the Prevotella genus in the gut, which is genomically heterogeneous and difficult to isolate; hence, scarce research was carried out for this species. This study aimed to investigate the effect of P. copri on hyperglycemia. Thirty-nine strains were isolated from healthy individuals, and three strains (HF2123, HF1478, and HF2130) that had the highest glucose consumption were selected to evaluate the effects of P. copri supplementation on hyperglycemia. Microbiomics and non-target metabolomics were used to uncover the underlying mechanisms. Oral administration of P. copri in diabetic db/db mice increased the expression and secretion of glucagon-like peptide-1 (GLP-1), significantly improved hyperglycemia, insulin resistance, and lipid accumulation, and alleviated the pathological morphology in the pancreas, liver, and colon. P. copri changed the composition of the gut microbiota of diabetic db/db mice, which was characterized by increasing the ratio of Bacteroidetes to Firmicutes and increasing the relative abundance of genera Bacteroides, Akkermansia, and Faecalibacterium. After intervention with P. copri, fecal metabolic profiling showed that fumaric acid and homocysteine contents decreased, and glutamine contents increased. Furthermore, amino acid metabolism and cAMP/PKA signaling pathways were enriched. Our findings indicate that P. copri improved glucose metabolism abnormalities in diabetic db/db mice. Especially, one of the P. copri strains, HF2130, has shown superior performance in improving hyperglycemia, which may have the potential as a probiotic against hyperglycemia. IMPORTANCE: As a core member of the human intestinal ecosystem, Prevotelal copri has been associated with glucose metabolic homeostasis in previous studies. However, these results have often been derived from metagenomic studies, and the experimental studies have been based solely on the type of strain DSM 18205T. Therefore, more experimental evidence from additional isolates is needed to validate the results according to their high genomic heterogeneity. In this study, we isolated different branches of strains and demonstrated that P. copri could improve the metabolic profile of hyperglycemic mice by modulating microbial activity. This finding supports the causal contribution of P. copri in host glucose metabolism.

Genomic Characteristion of Opportunistic Pathogen Kluyvera Reveals a Novel CTX-M Subgroup.

A rising incidence of clinical infections has been caused by Kluyvera, a significant opportunistic pathogen. Meanwhile, Kluyvera acts as an important reservoir of blaCTX-Ms, which are the dominant genes of class A extended-spectrum ß-lactamases (ESBLs). In this work, 60 strains of Kluyvera were subjected to phylogenetic relationship reconstruction, antimicrobial susceptibility testing, and antibiotic resistance genes prediction. All mature blaCTX-Ms were gathered to perform subgroup reclassification. The findings demonstrate that Kluyvera has a large gene pool with significant genetic flexibility. Notably, 25% of strains showed simultaneous detection of ESBLs and carbapenem resistance genes. The genotypes of fourteen novel blaCTX-Ms were identified. A new subgroup classification approach for blaCTX-Ms was defined by using 20 amino acid site variants, which could split blaCTX-Ms into 10 subgroups. The results of the subgroup division were consistent with the phylogenetic clustering. More significantly, we proposed a novel blaCTX-M subgroup, KLUS, that is chromosomally encoded in K. sichuanensis and the new species put forward in this study, showing amino acid differences from the currently known sequences. Cloning and transformation tests demonstrated that the recipient bacteria had a robust phenotype of cefotaxime resistance. Closely related Kluyvera species had blaCTX-Ms in the same subgroup. Our research lays the groundwork for a deeper comprehension of Kluyvera and emphasizes how important a blaCTX-M reservoir it is. We provide an update on blaCTX-M subgroups reclassification from the aspects of phylogenetic relationship, amino acid differences, and the new subgroup KLUS, which needs to be strengthen monitored due to its strong resistance phenotype to cefotaxime.