SKESA: strategic k-mer extension for scrupulous assemblies

A survey of k-mer methods and applications in bioinformatics

The rapid progression of genomics and proteomics has been driven by the advent of advanced sequencing technologies, large, diverse, and readily available omics datasets, and the evolution of computational data processing capabilities. The vast amount of data generated by these advancements necessitates efficient algorithms to extract meaningful information. K-mers serve as a valuable tool when working with large sequencing datasets, offering several advantages in computational speed and memory efficiency and carrying the potential for intrinsic biological functionality. This review provides an overview of the methods, applications, and significance of k-mers in genomic and proteomic data analyses, as well as the utility of absent sequences, including nullomers and nullpeptides, in disease detection, vaccine development, therapeutics, and forensic science. Therefore, the review highlights the pivotal role of k-mers in addressing current genomic and proteomic problems and underscores their potential for future breakthroughs in research.

Characterisation of glucose-induced protein fragments among the order Enterobacterales using matrix-assisted laser desorption ionization-time of flight mass spectrometry

To characterise the glucose-induced protein fragments by MALDI-TOF MS analysis, we compared data for samples from Escherichia coli cultured in media with or without glucose. Characteristic peaks were observed in the presence of glucose, and MS/MS revealed Asr-specific fragments. The amino acid sequences of the fragments suggested sequence-specific proteolysis. Blast-analysis revealed that numerous Enterobacterales harbored genes encoding Asr as well as E. coli. Here, we analysed 32 strains from 20 genera and 25 species of seven Enterobacterales families. We did not detect changes in the mass spectra of four strains of Morganellaceae lacking asr, whereas peaks of Asr-specific fragments were detected in the other 28 strains. We therefore concluded that the induction of Asr production in the presence of glucose is common among the Enterobacterales, except for certain Morganellaceae species. In members of family Budviciaceae, unfragmented Asr was detected. Molecular genetic information suggested that the amino acid sequences of Asr homologs are diverse, with fragments varying in number and size, indicating that Asr may serve as a discriminative biomarker for identifying Enterobacterales species.

A multicenter study on accuracy and reproducibility of nanopore sequencing-based genotyping of bacterial pathogens

Nanopore sequencing has shown the potential to democratize genomic pathogen surveillance due to its ease of use and low entry cost. However, recent genotyping studies showed discrepant results compared to gold-standard short-read sequencing. Furthermore, although essential for widespread application, the reproducibility of nanopore-only genotyping remains largely unresolved. In our multicenter performance study involving five laboratories, four public health-relevant bacterial species were sequenced with the latest R10.4.1 flow cells and V14 chemistry. Core genome MLST analysis of over 500 data sets revealed highly strain-specific typing errors in all species in each laboratory. Investigation of the methylation-related errors revealed consistent DNA motifs at error-prone sites across participants at read level. Depending on the frequency of incorrect target reads, this either leads to correct or incorrect typing, whereby only minimal frequency deviations can randomly determine the final result. PCR preamplification, recent basecalling model updates and an optimized polishing strategy notably diminished the non-reproducible typing. Our study highlights the potential for new errors to appear with each newly sequenced strain and lays the foundation for computational approaches to reduce such typing errors. In conclusion, our multicenter study shows the necessity for a new validation concept for nanopore sequencing-based, standardized bacterial typing, where single nucleotide accuracy is critical.

Draft genome sequences of three poultry Salmonella Shamba isolates from South Africa

Nontyphoidal Salmonella enterica serovars are foodborne pathogens commonly transmitted through poultry products. Draft genome sequences of three Salmonella enterica subsp. enterica serovar Shamba isolates which were obtained from poultry house dust in South Africa are reported herein.

Genomic epidemiology of mecC-carrying Staphylococcus aureus isolates from human clinical cases in New Zealand

In 2011, a novel methicillin resistance gene, mecC, was described in human and bovine Staphylococcus aureus isolates. mecC-positive S. aureus is most commonly associated with livestock and wildlife populations across Europe and is particularly prevalent in hedgehogs, but only occasionally causes human infections. In this study, we characterize and investigate the origin of two human S. aureus isolates containing mecC genes from New Zealand. The two isolates were identified from patients with severe invasion infections as part of an S. aureus bacteraemia study. Whole-genome sequencing was used to characterize staphylococcal cassette chromosome mec (SCCmec) elements and perform phylogenetic comparisons with publicly available strains from mecC-associated clonal complexes, including isolates from hedgehogs from New Zealand and Europe/United Kingdom (UK), and livestock, wildlife and human isolates from Europe/UK. The two isolates from our study have almost identical SCCmec type XI elements containing a mecC gene. However, this gene contains a premature stop codon, consistent with the methicillin-susceptible phenotype observed for these isolates. Core genome SNP analyses showed that the two isolates are 234 SNPs apart and are most closely related to an isolate obtained from a New Zealand hedgehog. However, there are considerable differences in the mecC mobile element between the human and hedgehog isolates, indicating the presence of an as-yet-unknown reservoir of mecC S. aureus in the New Zealand environment.

River waters in Greece: A reservoir for clinically relevant extended-spectrum-β-lactamases-producing Escherichia coli

In the current study, the genotypic characteristics such as antimicrobial resistance and virulence genes, and plasmid replicons and phenotypic characteristics such as biofilm formation and antimicrobial resistance of 87 extended-spectrum beta-lactamase (ESBL)-producing E. coli (ESBL-Ec) isolated from 7 water bodies in northern Greece were investigated. Our data show a high prevalence (60.0 %) of ESBL-Ec in surface waters that exhibit high genetic diversity, suggesting multiple sources of their transmission into the aquatic environment. When evaluating the antimicrobial resistance of isolates, wide variation in their resistance profiles has been detected, with all isolates being multi-drug resistant (MDR). Regarding biofilm formation capacity and phylogenetic groups, the majority (54.0 %, 47/87) of ESBL-Ec were classified as no biofilm producers mainly assigned to phylogroup A (35.6 %; 31/87), followed by B2 (26.5 %; 23/87). PCR screening showed that a high proportion of the isolates tested positive for the blaCTX-M-1 group genes (69 %, 60/87), followed by blaTEM (55.2 %, 48/87), blaOXA (25.3 %, 22/87) and blaCTX-M-9 (17.2 %, 15/87). A subset of 28 ESBL-Ec strains was further investigated by applying whole genome sequencing (WGS), and among them, certain clinically significant sequence types were identified, such as ST131 and ST10. The corresponding in silico analysis predicted all these isolates as human pathogens, while a significant proportion of WGS-ESBL-Ec were assigned to extraintestinal pathogenic E. coli (ExPEC; 32.1 %), and urinary pathogenic E. coli (UPEC; 28.6 %) pathotypes. Comparative phylogenetic analysis, showed that the genomes of the ST131-O25:H4-H30 isolates are genetically linked to the human clinical strains. Here, we report for the first time the detection of a plasmid-mediated mobile colistin resistance gene in ESBL-Ec in Greece isolated from an environmental source. Overall, this study underlines the role of surface waters as a reservoir for antibiotic resistance genes and for presumptive pathogenic ESBL-Ec.

Transmission of Brucella canis in a canine kennel following introduction of an infected dog

Brucella canis is a zoonotic pathogen and the main causative agent of canine brucellosis. In the Netherlands, B. canis had previously only been detected in individual cases of imported dogs. However, an outbreak of B. canis occurred for the first time in a cohort of autochthonous dogs in a breeding kennel in 2019. The outbreak began with a positive serological test result of an imported intact male dog showing clinical symptoms of brucellosis. Consequently, urine and blood samples were collected and tested positive for B. canis by culture, matrix-assisted laser desorption/ionization - time of flight mass spectrometry (MALDI-TOF MS) and whole-genome-sequencing (WGS). Screening of the contact dogs in the kennel where the index case was kept, revealed that antibodies against B. canis could be detected in 23 out of 69 dogs (34 %) by serum agglutination test (SAT). Of the 23 seropositive dogs, B. canis could be cultured from the urine and/or heparin samples of 19 dogs (83 %). This outbreak represents the first documented case of transmission of B. canis to autochthonous contact dogs in the Netherlands. WGS revealed all B. canis isolates belonged to the same cluster, which means the transmission of B. canis in the breeding kennel was most likely caused by the introduction of one infected dog. Comparing this cluster with data from other B. canis isolates, it also appears that characteristic clusters of B. canis are present in several endemic countries. These clusters seem to remain stable over time and may help in locating the origin of new isolates found. This outbreak showed that the international movement of dogs from endemic countries poses a threat to the canine population, while serological screening and WGS proved to be valuable tools for respectively screening and the epidemiological investigation.

Genomic insights on carotenoid synthesis by extremely halophilic archaea Haloarcula rubripromontorii BS2, Haloferax lucentense BBK2 and Halogeometricum borinquense E3 isolated from the solar salterns of India

Haloarchaeal cultures were isolated from solar salterns of Goa and Tamil Nadu and designated as BS2, BBK2 and E3. These isolates grew with a characteristic bright orange to pink pigmentation and were capable of growing in media containing upto 25% (w/vol) NaCl. Whole genome sequencing (WGS) of the three haloarchaeal strains BS2, BBK2 and E3 indicated an assembled genomic size of 4.1 Mb, 3.8 Mb and 4 Mb with G + C content of 61.8, 65.6 and 59.8% respectively. Phylogenetic analysis based on the 16S rRNA gene sequence revealed that the archaeal isolates belong to Haloarcula, Haloferax and Halogeometricum genera. Haloarcula rubripromontorii BS2  was predicted to have 4292 genes with 4242 CDS regions, 46 tRNAs, 6 rRNAs and 3 misc_RNAs. In case of Haloferax lucentense  BBK2,, 3840 genes with 3780 CDS regions were detected along with 52 tRNAs, 5 rRNAs and 3 misc_RNAs. Halogeometricum borinquense  E3 contained 4101 genes, 4043 CDS regions, 52 tRNAs, 4 rRNAs, and 2 misc_RNAs. The functional annotation and curation of the haloarchaeal genome, revealed C50 carotenoid biosynthetic genes like phytoene desaturase/carotenoid 3' -4' desaturase (crtI), lycopene elongase (ubiA/lyeJ) and carotenoid biosynthesis membrane protein (cruF) in the three isolates. Whereas crtD (C-3',4' desaturase), crtY (lycopene cyclase) and brp/blh (β-carotene dioxygenase) genes were identified only in BS2.

Draft genome of Thermoanaerobacter thermohydrosulfuricus strain AK152, a novel thermophilic and anaerobic bacterium isolated from a hot spring in Iceland

We present the draft genome of the bacterium Thermoanaerobacter thermohydrosulfuricus strain AK152, a thermophilic, endospore-spore-forming, anaerobe isolated from a hot spring in Grensdalur, in Southwestern Iceland. This assembled genome will lay the foundation for identifying the carboxylic and amino acid fermentation pathways, suggesting biotechnological applications for this strain.

Draft genome sequence of the fungal biocontrol agent, Bacillus velezensis Kos

Here, we report the draft genome sequence of Bacillus velezensis strain Kos, isolated from casing soil used during Agaricus bisporus cultivation in Dublin, Ireland. B. velezensis Kos exhibits a suppressive ability toward Cladobotryum mycophilum, Trichoderma aggressivum, and Lecanicillium fungicola, which are common threats to A. bisporus production, cultivation, and quality.

Evolutionary genomic analyses of canine E. coli infections identify a relic capsular locus associated with resistance to multiple classes of antimicrobials

Infections caused by antimicrobial-resistant Escherichia coli are the leading cause of death attributed to antimicrobial resistance (AMR) worldwide, and the known AMR mechanisms involve a range of functional proteins. Here, we employed a pan-genome wide association study (GWAS) approach on over 1,000 E. coli isolates from sick dogs collected across the US and Canada and identified a strong statistical association (empirical P < 0.01) of AMR, involving a range of antibiotics to a group 1 capsular (CPS) gene cluster. This cluster included genes under relaxed selection pressure, had several loci missing, and had pseudogenes for other key loci. Furthermore, this cluster is widespread in E. coli and Klebsiella clinical isolates across multiple host species. Earlier studies demonstrated that the octameric CPS polysaccharide export protein Wza can transmit macrolide antibiotics into the E. coli periplasm. We suggest that the CPS in question, and its highly divergent Wza, functions as an antibiotic trap, preventing antimicrobial penetration. We also highlight the high diversity of lineages circulating in dogs across all regions studied, the overlap with human lineages, and regional prevalence of resistance to multiple antimicrobial classes.

IMPORTANCE

Much of the human genomic epidemiology data available for E. coli mechanism discovery studies has been heavily biased toward shiga-toxin producing strains from humans and livestock. E. coli occupies many niches and produces a wide variety of other significant pathotypes, including some implicated in chronic disease. We hypothesized that since dogs tend to share similar strains with their owners and are treated with similar antibiotics, their pathogenic isolates will harbor unexplored AMR mechanisms of importance to humans as well as animals. By comparing over 1,000 genomes with in vitro antimicrobial susceptibility data from sick dogs across the US and Canada, we identified a strong multidrug resistance association with an operon that appears to have once conferred a type 1 capsule production system.

Genomic diversity of Campylobacter jejuni and Campylobacter coli isolated from the Ethiopian dairy supply chain

Campylobacteriosis outbreaks have previously been linked to dairy foods. While the genetic diversity of Campylobacter is well understood in high-income countries, it is largely unknown in low-income countries, such as Ethiopia. This study therefore aimed to conduct the first genomic characterization of Campylobacter isolates from the Ethiopian dairy supply chain to aid in future epidemiological studies. Fourteen C. jejuni and four C. coli isolates were whole genome sequenced using an Illumina platform. Sequences were analyzed using the bioinformatics tools in the GalaxyTrakr platform to identify MLST types, and single nucleotide polymorphisms, and infer phylogenetic relationships among the studied isolates. Assembled genomes were further screened to detect antimicrobial resistance and virulence gene sequences. Among 14 C. jejuni, ST 2084 and ST 51, which belong to the clonal complexes ST-353 and ST-443, respectively, were identified. Among the 4 sequenced C. coli isolates, two isolates belonged to ST 1628 and two to ST 830 from the clonal complex ST-828. The isolates of C. jejuni ST 2084 and ST 51 carried β-lactam resistance gene blaOXA-605, a fluoroquinolone resistance-associated mutation T86I in the gryA gene, and a macrolide resistance-associated mutation A103V in 50S L22. Only ST 2084 isolates carried the tetracycline resistance gene tetO. Conversely, all four C. coli ST 830 and ST 1628 isolates carried tetO, but only ST 1628 isolates also carried blaOXA-605. Lastly, C. jejuni ST 2084 isolates carried a total of 89 virulence genes, and ST 51 isolates carried up to 88 virulence genes. Among C. coli, ST 830 isolates carried 71 genes involved in virulence, whereas two ST 1628 isolates carried up to 82 genes involved in virulence. Isolates from all identified STs have previously been isolated from human clinical cases, demonstrating a potential food safety concern. This finding warrants further monitoring of Campylobacter in dairy foods in Ethiopia to better understand and manage the risks associated with Campylobacter contamination and transmission.

Development and implementation of a core genome multilocus sequence typing scheme for Yersinia enterocolitica: a tool for surveillance and outbreak detection

Yersinia enterocolitica (Y. enterocolitica) is the most frequent etiological agent of yersiniosis and has been responsible for several national outbreaks in Norway and elsewhere. A standardized high-resolution method, such as core genome Multilocus Sequence Typing (cgMLST), is needed for pathogen traceability at the national and international levels. In this study, we developed and implemented a cgMLST scheme for Y. enterocolitica. We designed a cgMLST scheme in SeqSphere + using high-quality genomes from different Y. enterocolitica biotype sublineages. The scheme was validated if more than 95% of targets were found across all tested Y. enterocolitica: 563 Norwegian genomes collected between 2012 and 2022 and 327 genomes from public data sets. We applied the scheme to known outbreaks to establish a threshold for identifying major complex types (CTs) based on the number of allelic differences. The final cgMLST scheme included 2,582 genes with a median of 97.9% (interquartile range 97.6%-98.8%) targets found across all tested genomes. Analysis of outbreaks identified all outbreak strains using single linkage clustering at four allelic differences. This threshold identified 311 unique CTs in Norway, of which CT18, CT12, and CT5 were identified as the most frequently associated with outbreaks. The cgMLST scheme showed a very good performance in typing Y. enterocolitica using diverse data sources and was able to identify outbreak clusters. We recommend the implementation of this scheme nationally and internationally to facilitate Y. enterocolitica surveillance and improve outbreak response in national and cross-border outbreaks.

Escherichia coli in urban marine sediments: interpreting virulence, biofilm formation, halotolerance, and antibiotic resistance to infer contamination or naturalization

Marine sediments have been suggested as a reservoir for pathogenic bacteria, including Escherichia coli. The origins, and properties promoting survival of E. coli in marine sediments (including osmotolerance, biofilm formation capacity, and antibiotic resistance), have not been well-characterized. Phenotypes and genotypes of 37 E. coli isolates from coastal marine sediments were characterized. The isolates were diverse: 30 sequence types were identified that have been previously documented in humans, livestock, and other animals. Virulence genes were found in all isolates, with more virulence genes found in isolates sampled from sediment closer to the effluent discharge point of a wastewater treatment plant. Antibiotic resistance was demonstrated phenotypically for one isolate, which also carried tetracycline resistance genes on a plasmid. Biofilm formation capacity varied for the different isolates, with most biofilm formed by phylogroup B1 isolates. All isolates were halotolerant, growing at 3.5% NaCl. This suggests that the properties of some isolates may facilitate survival in marine environments and can explain in part how marine sediments can be a reservoir for pathogenic E. coli. As disturbance of sediment could resuspend bacteria, this should be considered as a potential contributor to compromised bathing water quality at nearby beaches.

Rapid Reversal of Carbapenemase-Producing Pseudomonas aeruginosa Epidemiology from blaVIM- to blaNDM-harbouring Isolates in a Greek Tertiary Care Hospital

Carbapenemase-producing Pseudomonas aeruginosa strains present a specific geographical distribution regarding the type of carbapenemase-encoding genes that they harbor. For more than twenty years, VIM-type enzymes were the only major carbapenemases that were detected among P. aeruginosa isolates in Greece until the emergence of NDM-1-encoding P. aeruginosa in early 2023. In the present study, we present the rapid reversal of the carbapenemase-producing P. aeruginosa epidemiology from blaVIM- to blaNDM-harbouring isolates that occurred in our hospital since then. Between January 2023 and February 2024, 139 isolates tested positive for carbapenemase production with the NG-Test CARBA 5 immunochromatographic assay. Eight isolates were processed with the Hybrispot antimicrobial resistance direct flow chip molecular assay, and the first NDM-producing isolate was further analyzed through whole genome sequencing and bioinformatics analysis. Multiple resistance genes were detected by molecular techniques in accordance with the extensively drug-resistant phenotype. The isolate that was subjected to whole-genome sequencing belonged to the P. aeruginosa high-risk clone ST308, and the blaNDM was located in the chromosome in accordance with previously reported data. During the study period, NDM-producing isolates were increasingly detected, and only five months after their emergence, they overcame VIM producers. Our results indicate the potential of this new clone to spread rapidly and predominate within healthcare institutions, further restricting the already limited treatment options.

Methylomonas rivi sp. nov., Methylomonas rosea sp. nov., Methylomonas aurea sp. nov. and Methylomonas subterranea sp. nov., type I methane-oxidizing bacteria isolated from a freshwater creek and the deep terrestrial subsurface

Four methane-oxidizing bacteria, designated as strains WSC-6T, WSC-7T, SURF-1T, and SURF-2T, were isolated from Saddle Mountain Creek in southwestern Oklahoma, USA, and the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The strains were Gram-negative, motile, short rods that possessed intracytoplasmic membranes characteristic of type I methanotrophs. All four strains were oxidase-negative and weakly catalase-positive. Colonies ranged from pale pink to orange in colour. Methane and methanol were the only compounds that could serve as carbon and energy sources for growth. Strains WSC-6T and WSC-7T grew optimally at lower temperatures (25 and 20 °C, respectively) compared to strains SURF-1T and SURF-2T (40 °C). Strains WSC-6T and SURF-2T were neutrophilic (optimal pH of 7.5 and 7.3, respectively), while strains WSC-7T and SURF-1T were slightly alkaliphilic, with an optimal pH of 8.8. The strains grew best in media amended with ≤0.5% NaCl. The major cellular fatty acids were C14 : 0, C16 : 1  ω8c, C16 : 1  ω7c, and C16 : 1  ω5c. The DNA G+C content ranged from 51.5 to 56.0 mol%. Phylogenetic analyses indicated that the strains belonged to the genus Methylomonas, with each exhibiting 98.6-99.6% 16S rRNA gene sequence similarity to closely related strains. Genome-wide estimates of relatedness (84.5-88.4% average nucleotide identity, 85.8-92.4% average amino acid identity and 27.4-35.0% digital DNA-DNA hybridization) fell below established thresholds for species delineation. Based on these combined results, we propose to classify these strains as representing novel species of the genus Methylomonas, for which the names Methylomonas rivi (type strain WSC-6T=ATCC TSD-251T=DSM 112293T), Methylomonas rosea (type strain WSC-7T=ATCC TSD-252T=DSM 112281T), Methylomonas aurea (type strain SURF-1T=ATCC TSD-253T=DSM 112282T), and Methylomonas subterranea (type strain SURF-2T=ATCC TSD-254T=DSM 112283T) are proposed.

APH(3')-Ie, an aminoglycoside-modifying enzyme discovered in a rabbit-derived Citrobacter gillenii isolate

Background

Aminoglycoside-modifying enzymes (AMEs) play an essential role in bacterial resistance to aminoglycoside antimicrobials. With the development of sequencing techniques, more bacterial genomes have been sequenced, which has aided in the discovery of an increasing number of novel resistance mechanisms.

Methods

The bacterial species was identified by 16S rRNA gene homology and average nucleotide identity (ANI) analyses. The minimum inhibitory concentration (MIC) of each antimicrobial was determined by the agar dilution method. The protein was expressed with the pCold I vector in E. coli BL21, and enzyme kinetic parameters were examined. The whole-genome sequence of the bacterium was obtained via the Illumina and PacBio sequencing platforms. Reconstruction of the phylogenetic tree, identification of conserved functional residues, and gene context analysis were performed using the corresponding bioinformatic techniques.

Results

A novel aminoglycoside resistance gene, designated aph(3')-Ie, which confers resistance to ribostamycin, kanamycin, sisomicin and paromomycin, was identified in the chromosome of the animal bacterium Citrobacter gillenii DW61, which exhibited a multidrug resistance phenotype. APH(3')-Ie showed the highest amino acid identity of 74.90% with the functionally characterized enzyme APH(3')-Ia. Enzyme kinetics analysis demonstrated that it had phosphorylation activity toward four aminoglycoside substrates, exhibiting the highest affinity (K m, 4.22 ± 0.88 µM) and the highest catalytic efficiency [k cat/K m, (32.27 ± 8.14) × 104] for ribomycin. Similar to the other APH(3') proteins, APH(3')-Ie contained all the conserved functional sites of the APH family. The aph(3')-Ie homologous genes were present in C. gillenii isolates from different sources, including some of clinical significance.

Conclusion

In this work, a novel chromosomal aminoglycoside resistance gene, designated aph(3')-Ie, conferring resistance to aminoglycoside antimicrobials, was identified in a rabbit isolate C. gillenii DW61. The elucidation of the novel resistance mechanism will aid in the effective treatment of infections caused by pathogens carrying such resistance genes.

Bacterial and Genetic Features of Raw Retail Pork Meat: Integrative Analysis of Antibiotic Susceptibility, Whole-Genome Sequencing, and Metagenomics

The global antibiotic resistance crisis, driven by overuse and misuse of antibiotics, is multifaceted. This study aimed to assess the microbiological and genetic characteristics of raw retail pork meat through various methods, including the isolation, antibiotic susceptibility testing (AST), whole-genome sequencing (WGS) of selected indicator bacteria, antibiotic residue testing, and metagenomic sequencing. Samples were purchased from 10 pre-selected retail stores in Gauteng, South Africa. The samples were aseptically separated, with portions sent to an external laboratory for isolating indicator bacteria and testing for antibiotic residues. Identification of the isolated bacteria was reconfirmed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). AST was performed using the Microscan Walkaway system (Beckman Coulter, Brea, CA, USA). WGS and metagenomic sequencing were performed using the Illumina NextSeq 550 instrument (San Diego, CA, USA). The isolated E. coli and E. faecalis exhibited minimal phenotypic resistance, with WGS revealing the presence of tetracycline resistance genes. Both the isolated bacteria and meat samples harboured tetracycline resistance genes and the antibiotic residue concentrations were within acceptable limits for human consumption. In the metagenomic context, most identified bacteria were of food/meat spoilage and environmental origin. The resistome analysis primarily indicated beta-lactam, tetracycline and multidrug resistance genes. Further research is needed to understand the broader implications of these findings on environmental health and antibiotic resistance.

Average nucleotide identity-based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital- and community-acquired infections in humans worldwide. Due to the high incidence of infection, S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole-genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of average nucleotide identity revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10%-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location. Notably, strain-diffuse genes were associated with prophages; strain-concentrated genes were associated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic resistance genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core, and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

IMPORTANCE

We analyzed the genomic diversity of Staphylococcus aureus, a globally prevalent bacterial species that causes serious infections in humans. Our goal was to build a genetic picture of the different strains of S. aureus and which genes may be associated with them. We reprocessed >84,000 genomes and subsampled to remove redundancy. We found that individual samples sharing >99.5% of their genome could be grouped into strains. We also showed that a portion of genes that are present in intermediate frequency in the species are strongly associated with some strains but completely absent from others, suggesting a role in strain specificity. This work lays the foundation for understanding individual gene histories of the S. aureus species and also outlines strategies for processing large bacterial genomic data sets.

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

A multicenter genomic epidemiological investigation in Brazil, Chile, and Mexico reveals the diversity and persistence of Salmonella populations in surface waters

This study examined the diversity and persistence of Salmonella in the surface waters of agricultural regions of Brazil, Chile, and Mexico. Research groups (three in 2019-2020 and five in 2021-2022) conducted a long-term survey of surface water across 5-8 months annually (n = 30 monthly). On-site, each team filtered 10-L water samples with modified Moore Swabs to capture Salmonella, which were then isolated and identified using conventional microbiological techniques. Salmonella isolates were sequenced on Illumina platforms. Salmonella was present in 1,493/3,291 water samples (45.8%), with varying isolation rates across countries and years. Newport, Infantis, and Typhimurium were the most frequent among the 128 different serovars. Notably, 22 serovars were found in all three countries, representing almost half of the 1,911 different isolates collected. The resistome comprised 72 antimicrobial resistance (AMR) genes and six point mutations in three genes. At least one AMR determinant was observed in 33.8% (646/1,911) of the isolates, of which 47.4% (306/646) were potentially multidrug resistant. Phylogeny based on core genome multilocus sequence typing (cgMLST) showed that most isolates clustered according to sequence type and country of origin. Only 14 cgMLST multi-country clusters were detected among the 275 clusters. However, further analysis confirmed that close genetic relatedness occurred mostly among isolates from the same country, with three exceptions. Interestingly, isolates closely related phylogenetically were recovered over multiple years within the same country, indicating the persistence of certain Salmonella in those areas. In conclusion, surface waters in these regions are consistently contaminated with diverse Salmonella, including strains that persist over time.IMPORTANCESalmonella is a leading foodborne pathogen responsible for millions of illnesses, hospitalizations, and deaths annually. Although Salmonella-contaminated water has now been recognized as an important contamination source in the agrifood chain, there is a lack of knowledge on the global occurrence and diversity of Salmonella in surface water. Moreover, there has been insufficient research on Salmonella in surface waters from Latin American countries that are major producers and exporters of agricultural products. Incorporating genetic profiling of Salmonella isolates from underrepresented regions, such as Latin America, enhances our understanding of the pathogen's ecology, evolution, antimicrobial resistance, and pathogenicity. Moreover, leveraging genomic data derived from pathogens isolated from diverse geographical areas is critical for assessing the potential public health risk posed by the pathogen and expediting investigations of foodborne outbreaks. Ultimately, global efforts contribute significantly to reducing the incidence of foodborne infections.

Genetic Homogeneity of Francisella tularensis subsp. mediasiatica Strains in Kazakhstan

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.

Overcoming Microbial Inhibition of S. Sonnei Through the Exploitation of Genomically Predicted Antibiotic Resistance Profiles for the Development of Food Enrichment Media

Linking outbreaks of Shigella spp. to specific foods is challenging due to poor selectivity of current enrichment media. We have previously shown that enrichment media, tailored to the genomically-predicted antimicrobial resistance (AMR) of Shiga toxigenic E. coli strains, enhances their isolation from foods. This study investigates the application of this approach for Shigella isolation. The AMR gene profiles of 21,908 published S. sonnei genomes indicated a high prevalence of genes conferring resistance to streptomycin (aadA, aph(3″)-Ib, aph(6)-Id, 92.8%), sulfonamides (sul1, sul2, 74.8%), and/or trimethoprim (dfrA, 96.2%). Genomic analysis and antibiotic susceptibility testing conducted with a panel of 17 outbreak-associated S. sonnei strains confirmed the correlation of AMR gene detection with resistance phenotypes. Supplementation of Shigella Broth (SB) with up to 400 µg/mL of trimethoprim or sulfadiazine did not suppress the growth of sensitive strains, whereas 100 µg/mL of streptomycin increased the selectivity of this broth. All three antibiotics increased the selectivity of modified Tryptone Soya Broth (mTSB). Based on these results, supplemented media formulations were developed and assessed by measuring the relative growth of S. sonnei in cultures coinoculated with a strain of bacteriocin-producing E. coli that is inhibitory to Shigella growth. S. sonnei was not recovered from cocultures grown in SB or mTSB without antibiotics. In contrast, media supplemented with streptomycin at 50 and 100 µg/mL, trimethoprim at 25 and 50 µg/mL, and sulfadiazine at 100 µg/mL increased the relative proportion of S. sonnei in postenrichment cultures. The enhanced recovery of resistant S. sonnei strains achieved in this study indicates that, in cases where genomic data are available for clinical S. sonnei isolates, customization of selective enrichment media based on AMR gene detection could be a valuable tool for supporting the investigation of foodborne shigellosis outbreaks.

Lacrimispora brassicae sp. nov. isolated from fermented cabbage, and proposal of Clostridium indicum Gundawar et al. 2019 and Clostridium methoxybenzovorans Mechichi et al. 1999 as heterotypic synonyms of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Lacrimispora indolis (McClung and McCoy 1957) Haas and Blanchard 2020, respectively

A Gram-stain-negative, endospore-forming, rod-shaped, indole-producing bacterial strain, designated YZC6T, was isolated from fermented cabbage. Strain YZC6T grew at 10-37  °C, pH 5.5-8.5, and with up to 2  % (w/v) NaCl. The major cellular fatty acids were C16 : 0 and C18 : 1 cis 11 dimethyl acetal. Phylogenetic analysis of the 16S rRNA gene revealed that strain YZC6T belonged to the genus Lacrimispora and was closely related to Lacrimispora aerotolerans DSM 5434T (98.3  % sequence similarity), Lacrimispora saccharolytica WM1T (98.1  %), and Lacrimispora algidixylanolytica SPL73T (98.1  %). The average nucleotide identity based on blast (below 87.8  %) and digital DNA-DNA hybridization (below 36.1 %) values between the novel isolate and its corresponding relatives showed that strain YZC6T could be readily distinguished from its closely related species. Based on genotypic, phenotypic, and chemotaxonomic data, a novel Lacrimispora species, Lacrimispora brassicae sp. nov., was proposed, with YZC6T as the type strain (=MAFF 212518T=JCM 32810T=DSM 112100T). This study also proposed Clostridium indicum Gundawar et al. 2019 as a later heterotypic synonym of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Clostridium methoxybenzovorans Mechichi et al. 1999 as a later heterotypic synonym of Lacrimispora indolis (McClung and McCpy 1957) Haas and Blanchard 2020.

What's on a prophage: analysis of Salmonella spp. prophages identifies a diverse range of cargo with multiple virulence- and metabolism-associated functions

The gain of mobile elements, such as prophages, can introduce cargo to the recipient bacterium that could facilitate its persistence in or expansion to a new environment, such as a host. While previous studies have focused on identifying and characterizing the genetic diversity of prophages, analyses characterizing the cargo that prophages carry have not been extensively explored. We characterized prophage regions from 303 Salmonella spp. genomes (representing 254 unique serovars) to assess the distribution of prophages in diverse Salmonella. On average, prophages accounted for 3.7% (0.1%-8.8%) of the total genomic content of each isolate. Prophage regions annotated as Gifsy 1 and Salmon Fels 1 were the most commonly identified intact prophages, suggesting that they are common throughout the Salmonella genus. Among 21,687 total coding sequences (CDSs) from intact prophage regions in subsp. enterica genomes, 7.5% (median; range: 1.1%-47.6%) were categorized as having a function not related to prophage integration or phage structure, some of which could potentially provide a functional attribute to the host Salmonella cell. These predicted functions could be broadly categorized into CDSs involved in: (i) modification of cell surface structures (i.e., glycosyltransferases); (ii) modulation of host responses (e.g., SodC/SodA, SopE, ArtAB, and typhoid toxin); (iii) conferring resistance to heavy metals and antimicrobials; (iv) metabolism of carbohydrates, amino acids, and nucleotides; and (v) DNA replication, repair, and regulation. Overall, our systematic analysis of prophage cargo highlights a broader role for prophage cargo in influencing the metabolic, virulence, and resistance characteristics of Salmonella.

IMPORTANCE

Lysogenic bacteriophages (phages) can integrate their genome into a bacterial host's genome, potentially introducing genetic elements that can affect the fitness of the host bacterium. The functions of prophage-encoded genes are important to understand as these genes could be mobilized and transferred to a new host. Using a large genomic dataset representing >300 isolates from all known subspecies and species of Salmonella, our study contributes important new findings on the distribution of prophages and the types of cargo that diverse Salmonella prophages carry. We identified a number of coding sequences (CDSs) annotated as having cell surface-modifying attributes, suggesting that prophages may have played an important role in shaping Salmonella's diverse surface antigen repertoire. Furthermore, our characterization of prophages suggests that they play a broader role in facilitating the acquisition and transfer of CDSs associated with metabolism, DNA replication and repair, virulence factors, and to a lesser extent, antimicrobial resistance.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

Core genome multilocus sequence typing (cgMLST) applicable to the monophyletic Klebsiella oxytoca species complex

The environmental bacterium Klebsiella oxytoca displays an alarming increase of antibiotic-resistant strains that frequently cause outbreaks in intensive care units. Due to its prevalence in the environment and opportunistic presence in humans, molecular surveillance (including resistance marker screening) and high-resolution cluster analysis are of high relevance. Furthermore, K. oxytoca previously described in studies is rather a species complex (KoSC) than a single species comprising at least six closely related species that are not easily differentiated by standard typing methods. To reach a discriminatory power high enough to identify and resolve clusters within these species, whole genome sequencing is necessary. The resolution is achievable with core genome multilocus sequence typing (cgMLST) extending typing of a few housekeeping genes to thousands of core genome genes. CgMLST is highly standardized and provides a nomenclature enabling cross laboratory reproducibility and data exchange for routine diagnostics. Here, we established a cgMLST scheme not only capable of resolving the KoSC species but also producing reliable and consistent results for published outbreaks. Our cgMLST scheme consists of 2,536 core genome and 2,693 accessory genome targets, with a percentage of good cgMLST targets of 98.31% in 880 KoSC genomes downloaded from the National Center for Biotechnology Information (NCBI). We also validated resistance markers against known resistance gene patterns and successfully linked genetic results to phenotypically confirmed toxic strains carrying the til gene cluster. In conclusion, our novel cgMLST enables highly reproducible typing of four different clinically relevant species of the KoSC and thus facilitates molecular surveillance and cluster investigations.

Complete genome sequence of Yokenella regensburgei isolated from a patient with urinary tract infection in India

Yokenella regensburgei, an environmental organism, is an emerging pathogen in patients chiefly with immune suppression. We report the draft genome of Y. regensburgei, strain UU2206353, isolated from the urinary tract of an immunocompetent individual. The assembled genome consisted of 4,669,536 bp distributed over 20 contigs with 4,283 protein-coding genes.

The GEA pipeline for characterizing Escherichia coli and Salmonella genomes

Salmonella enterica and Escherichia coli are major food-borne human pathogens, and their genomes are routinely sequenced for clinical surveillance. Computational pipelines designed for analyzing pathogen genomes should both utilize the most current information from annotation databases and increase the coverage of these databases over time. We report the development of the GEA pipeline to analyze large batches of E. coli and S. enterica genomes. The GEA pipeline takes as input paired Illumina raw reads files which are then assembled followed by annotation. Alternatively, assemblies can be provided as input and directly annotated. The pipeline provides predictive genome annotations for E. coli and S. enterica with a focus on the Center for Genomic Epidemiology tools. Annotation results are provided as a tab delimited text file. The GEA pipeline is designed for large-scale E. coli and S. enterica genome assembly and characterization using the Center for Genomic Epidemiology command-line tools and high-performance computing. Large scale annotation is demonstrated by an analysis of more than 14,000 Salmonella genome assemblies. Testing the GEA pipeline on E. coli raw reads demonstrates reproducibility across multiple compute environments and computational usage is optimized on high performance computers.

Comparative transcriptomic analysis of Staphylococcus epidermidis associated with periprosthetic joint infection under in vivo and in vitro conditions

Staphylococcus epidermidis is part of the commensal microbiota of the skin and mucous membranes, though it can also act as a pathogen in certain scenarios, causing a range of infections, including periprosthetic joint infection (PJI). Transcriptomic profiling may provide insights into mechanisms by which S. epidermidis adapts while in a pathogenic compared to a commensal state. Here, a total RNA-sequencing approach was used to profile and compare the transcriptomes of 19 paired PJI-associated S. epidermidis samples from an in vivo clinical source and grown in in vitro laboratory culture. Genomic comparison of PJI-associated and publicly available commensal-state isolates were also compared. Of the 1919 total transcripts found, 145 were from differentially expressed genes (DEGs) when comparing in vivo or in vitro samples. Forty-two transcripts were upregulated and 103 downregulated in in vivo samples. Of note, metal sequestration-associated genes, specifically those related to staphylopine activity (cntA, cntK, cntL, and cntM), were upregulated in a subset of clinical in vivo compared to laboratory grown in vitro samples. About 70% of the total transcripts and almost 50% of the DEGs identified have not yet been annotated. There were no significant genomic differences between known commensal and PJI-associated S. epidermidis isolates, suggesting that differential genomics may not play a role in S. epidermidis pathogenicity. In conclusion, this study provides insights into phenotypic alterations employed by S epidermidis to adapt to infective and non-infected microenvironments, potentially informing future therapeutic targets for related infections.

Phenotypic and genomic characterization of Methanothermobacter wolfeii strain BSEL, a CO2-capturing archaeon with minimal nutrient requirements

A new variant of Methanothermobacter wolfeii was isolated from an anaerobic digester using enrichment cultivation in anaerobic conditions. The new isolate was taxonomically identified via 16S rRNA gene sequencing and tagged as M. wolfeii BSEL. The whole genome of the new variant was sequenced and de novo assembled. Genomic variations between the BSEL strain and the type strain were discovered, suggesting evolutionary adaptations of the BSEL strain that conferred advantages while growing under a low concentration of nutrients. M. wolfeii BSEL displayed the highest specific growth rate ever reported for the wolfeii species (0.27 ± 0.03 h-1) using carbon dioxide (CO2) as unique carbon source and hydrogen (H2) as electron donor. M. wolfeii BSEL grew at this rate in an environment with ammonium (NH4+) as sole nitrogen source. The minerals content required to cultivate the BSEL strain was relatively low and resembled the ionic background of tap water without mineral supplements. Optimum growth rate for the new isolate was observed at 64°C and pH 8.3. In this work, it was shown that wastewater from a wastewater treatment facility can be used as a low-cost alternative medium to cultivate M. wolfeii BSEL. Continuous gas fermentation fed with a synthetic biogas mimic along with H2 in a bubble column bioreactor using M. wolfeii BSEL as biocatalyst resulted in a CO2 conversion efficiency of 97% and a final methane (CH4) titer of 98.5%v, demonstrating the ability of the new strain for upgrading biogas to renewable natural gas.IMPORTANCEAs a methanogenic archaeon, Methanothermobacter wolfeii uses CO2 as electron acceptor, producing CH4 as final product. The metabolism of M. wolfeii can be harnessed to capture CO2 from industrial emissions, besides producing a drop-in renewable biofuel to substitute fossil natural gas. If used as biocatalyst in new-generation CO2 sequestration processes, M. wolfeii has the potential to accelerate the decarbonization of the energy generation sector, which is the biggest contributor of CO2 emissions worldwide. Nonetheless, the development of CO2 sequestration archaeal-based biotechnology is still limited by an uncertainty in the requirements to cultivate methanogenic archaea and the unknown longevity of archaeal cultures. In this study, we report the adaptation, isolation, and phenotypic characterization of a novel variant of M. wolfeii, which is capable of maximum growth with minimal nutrients input. Our findings demonstrate the potential of this variant for the production of renewable natural gas, paving the way for the development of more efficient and sustainable CO2 sequestration processes.

AMRViz enables seamless genomics analysis and visualization of antimicrobial resistance

We have developed AMRViz, a toolkit for analyzing, visualizing, and managing bacterial genomics samples. The toolkit is bundled with the current best practice analysis pipeline allowing researchers to perform comprehensive analysis of a collection of samples directly from raw sequencing data with a single command line. The analysis results in a report showing the genome structure, genome annotations, antibiotic resistance and virulence profile for each sample. The pan-genome of all samples of the collection is analyzed to identify core- and accessory-genes. Phylogenies of the whole genome as well as all gene clusters are also generated. The toolkit provides a web-based visualization dashboard allowing researchers to interactively examine various aspects of the analysis results. Availability: AMRViz is implemented in Python and NodeJS, and is publicly available under open source MIT license at https://github.com/amromics/amrviz .

Gut-associated functions are favored during microbiome assembly across a major part of C. elegans life

The microbiome expresses a variety of functions that influence host biology. The range of functions depends on the microbiome's composition, which can change during the host's lifetime due to neutral assembly processes, host-mediated selection, and environmental conditions. To date, the exact dynamics of microbiome assembly, the underlying determinants, and the effects on host-associated functions remain poorly understood. Here, we used the nematode Caenorhabditis elegans and a defined community of fully sequenced, naturally associated bacteria to study microbiome dynamics and functions across a major part of the worm's lifetime of hosts under controlled experimental conditions. Bacterial community composition initially shows strongly declining levels of stochasticity, which increases during later time points, suggesting selective effects in younger animals as opposed to more random processes in older animals. The adult microbiome is enriched in genera Ochrobactrum and Enterobacter compared to the direct substrate and a host-free control environment. Using pathway analysis, metabolic, and ecological modeling, we further find that the lifetime assembly dynamics increase competitive strategies and gut-associated functions in the host-associated microbiome, indicating that the colonizing bacteria benefit the worm. Overall, our study introduces a framework for studying microbiome assembly dynamics based on stochastic, ecological, and metabolic models, yielding new insights into the processes that determine host-associated microbiome composition and function.

IMPORTANCE

The microbiome plays a crucial role in host biology. Its functions depend on the microbiome composition that can change during a host's lifetime. To date, the dynamics of microbiome assembly and the resulting functions still need to be better understood. This study introduces a new approach to characterize the functional consequences of microbiome assembly by modeling both the relevance of stochastic processes and metabolic characteristics of microbial community changes. The approach was applied to experimental time-series data obtained for the microbiome of the nematode Caenorhabditis elegans across the major part of its lifetime. Stochastic processes played a minor role, whereas beneficial bacteria as well as gut-associated functions enriched in hosts. This indicates that the host might actively shape the composition of its microbiome. Overall, this study provides a framework for studying microbiome assembly dynamics and yields new insights into C. elegans microbiome functions.

Prevalence and genetic diversity of azole-resistant Malassezia pachydermatis isolates from canine otitis and dermatitis: A 2-year study

Despite previous reports on the emergence of Malassezia pachydermatis strains with decreased susceptibility to azoles, there is limited information on the actual prevalence and genetic diversity of azole-resistant isolates of this yeast species. We assessed the prevalence of azole resistance in M. pachydermatis isolates from cases of dog otitis or skin disease attended in a veterinary teaching hospital during a 2-year period and analyzed the ERG11 (encoding a lanosterol 14-α demethylase, the primary target of azoles) and whole genome sequence diversity of a group of isolates that displayed reduced azole susceptibility. Susceptibility testing of 89 M. pachydermatis isolates from 54 clinical episodes (1-6 isolates/episode) revealed low minimum inhibitory concentrations (MICs) to most azoles and other antifungals, but 11 isolates from six different episodes (i.e., 12.4% of isolates and 11.1% of episodes) had decreased susceptibility to multiple azoles (fluconazole, itraconazole, ketoconazole, posaconazole, ravuconazole, and/or voriconazole). ERG11 sequencing of these 11 azole-resistant isolates identified eight DNA sequence profiles, most of which contained amino acid substitutions also found in some azole-susceptible isolates. Analysis of whole genome sequencing (WGS) results revealed that the azole-resistant isolates from the same episode of otitis, or even different episodes affecting the same animal, were more genetically related to each other than to isolates from other dogs. In conclusion, our results confirmed the remarkable ERG11 sequence variability in M. pachydermatis isolates of animal origin observed in previous studies and demonstrated the value of WGS for disentangling the epidemiology of this yeast species.

Genomic characterization of antimicrobial resistance in 61 aquatic bacterial isolates

Antimicrobial resistance (AMR) in pathogens important to aquatic animal health is of increasing concern but vastly understudied. Antimicrobial therapy is used to both treat and prevent bacterial disease in fish and is critical for a viable aquaculture industry and for maintenance of wild fish populations. Unfortunately, phenotypic antimicrobial susceptibility testing is technically difficult for bacteria recovered from aquatic animal hosts resulting in challenges in resistance monitoring using traditional methods. Whole-genome sequencing provides an appealing methodology for investigation of putative resistance. As part of the ongoing efforts of the FDA CVM Vet-LIRN to monitor AMR, source laboratories cultured and preliminarily identified pathogenic bacteria isolated from various fish species collected in 2019 from across the United States. Sixty-one bacterial isolates were evaluated using whole-genome sequencing. We present here the assembled draft genomes, AMR genes, predicted resistance phenotypes, and virulence factors of the 61 isolates and discuss concurrence of the identifications made by source laboratories using matrix-assisted laser desorption/time-of-flight mass spectrometry.

Characterization of MLST-99 Salmonella Typhimurium and the monophasic variant I:4,[5],12:i:- isolated from Canadian Atlantic coast shellfish

Salmonella enterica subsp. enterica Typhimurium and its monophasic variant I 1;4,[5],12:i:- (MVST) are responsible for thousands of reported cases of salmonellosis each year in Canada, and countries worldwide. We investigated S. Typhimurium and MVST isolates recovered from raw shellfish harvested in Atlantic Canada by the Canadian Food Inspection Agency (CFIA) over the past decade, to assess the potential impact of these isolates on human illness and to explore possible routes of shellfish contamination. Whole-genome sequence analysis was performed on 210 isolates of S. Typhimurium and MVST recovered from various food sources, including shellfish. The objective was to identify genetic markers linked to ST-99, a sequence type specifically associated with shellfish, which could explain their high prevalence in shellfish. We also investigated the genetic similarity amongst CFIA ST-99 isolates recovered in different years and geographical locations. Finally, the study aimed to enhance the molecular serotyping of ST-99 isolates, as they are serologically classified as MVST but are frequently misidentified as S. Typhimurium through sequence analysis. To ensure recovery of ST-99 from shellfish was not due to favourable growth kinetics, we measured the growth rates of these isolates relative to other Salmonella and determined that ST-99 did not have a faster growth rate and/or shorter lag phase than other Salmonella evaluated. The CFIA ST-99 isolates from shellfish were highly clonal, with up to 81 high-quality single nucleotide variants amongst isolates. ST-99 isolates both within the CFIA collection and those isolated globally carried numerous unique deletions, insertions and mutations in genes, including some considered important for virulence, such as gene deletions in the type VI secretion system. Interestingly, several of these genetic characteristics appear to be unique to North America. Most notably was a large genomic region showing a high prevalence in genomes from Canadian isolates compared to those from the USA. Although the functions of the majority of the proteins encoded within this region remain unknown, the genes umuC and umuD, known to be protective against UV light damage, were present. While this study did not specifically examine the effects of mutations and insertions, results indicate that these isolates may be adapted to survive in specific environments, such as ocean water, where wild birds and/or animals serve as the natural hosts. Our hypothesis is reinforced by a global phylogenetic analysis, which indicates that isolates obtained from North American shellfish and wild birds are infrequently connected to isolates from human sources. These findings suggest a distinct ecological niche for ST-99, potentially indicating their specialization and adaptation to non-human hosts and environments, such as oceanic habitats.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

From Investigating a Case of Cellulitis to Exploring Nosocomial Infection Control of ST1 Legionella pneumophila Using Genomic Approaches

Legionella pneumophila can cause a large panel of symptoms besides the classic pneumonia presentation. Here we present a case of fatal nosocomial cellulitis in an immunocompromised patient followed, a year later, by a second case of Legionnaires' disease in the same ward. While the first case was easily assumed as nosocomial based on the date of symptom onset, the second case required clear typing results to be assigned either as nosocomial and related to the same environmental source as the first case, or community acquired. To untangle this specific question, we applied core-genome multilocus typing (MLST), whole-genome single nucleotide polymorphism and whole-genome MLST methods to a collection of 36 Belgian and 41 international sequence-type 1 (ST1) isolates using both thresholds recommended in the literature and tailored threshold based on local epidemiological data. Based on the thresholds applied to cluster isolates together, the three methods gave different results and no firm conclusion about the nosocomial setting of the second case could been drawn. Our data highlight that despite promising results in the study of outbreaks and for large-scale epidemiological investigations, next-generation sequencing typing methods applied to ST1 outbreak investigation still need standardization regarding both wet-lab protocols and bioinformatics. A deeper evaluation of the L. pneumophila evolutionary clock is also required to increase our understanding of genomic differences between isolates sampled during a clinical infection and in the environment.

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.

Many purported pseudogenes in bacterial genomes are bona fide genes

BACKGROUND

Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly.

RESULTS

Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors.

CONCLUSIONS

Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

Comprehensive whole genome analysis of Staphylococcus aureus isolates from dairy cows with subclinical mastitis

Staphylococcus species are the primary cause of mastitis in dairy cows across the world. Staphylococcus aureus has recently become a pathogen that is zoonotic and multidrug resistant. This study aimed to sequence whole genomes of 38 S. aureus isolates from 55 subclinical mastitis dairy cows of 7 small-scale farmers in the Free State Province, South Africa and document and their antimicrobial and virulence genes. The 38 isolates were grouped by the in silico multi-locus sequencing types (MLST) into seven sequence types (STs), that is (ST 97, 352, 152, 243) and three new STs (ST8495, ST8500, and ST8501). Thirty-three S. aureus isolates were divided into 7 core single-nucleotide polymorphism (SNP) clusters. Among the 9 distinct spa-types that were detected, Spa-types t2883 accounted for the majority of isolates at 12 (31.57%), followed by t416 with 11 (28.94%) and t2844 with 5 (13.15%). The data also revealed the identification of four (4) plasmids, with Rep_N (rep20) accounting for the majority of isolates with 17 (44.73%), followed by Inc18 (repUS5) with 2 (5.26%). These isolates included 11 distinct antimicrobial resistance genes and 23 genes linked to bacterial virulence. Surprisingly, no methicillin resistance associated genes were detected in these isolates. Genome data of the current study will contribute to understanding epidemiology S. aureus genotypes and ultimately aid in developing treatment and control plans to stop the spread of mastitis in the Free State province and South Africa as a whole.

Antimicrobial Resistance in Selected Bacteria from Food Animals in New Zealand 2018-2022

Antimicrobial resistance (AMR) presents a significant threat to human health worldwide. One important source of antimicrobial-resistant infections in humans is exposure to animals or animal products. In a phased survey, we investigated AMR in 300 Escherichia coli isolates and 300 enterococci (Enterococcus faecalis and E. faecium) isolates each from the carcasses of poultry, pigs, very young calves, and dairy cattle (food animals); all Salmonella isolates from poultry, very young calves, and dairy cattle; and 300 Campylobacter (Campylobacter jejuni and C. coli) isolates from poultry. The highest resistance levels in E. coli were found for sulfamethoxazole, tetracycline, and streptomycin, for all food animals. Cefotaxime-resistant E. coli were not found and low resistance to ciprofloxacin, colistin, and gentamicin was observed. The majority of enterococci isolates from all food animals were bacitracin-resistant. Erythromycin- and/or tetracycline-resistant enterococci isolates were found in varying proportions from all food animals. Ampicillin- or vancomycin-resistant enterococci isolates were not identified, and ciprofloxacin-resistant E. faecalis were not found. Salmonella isolates were only recovered from very young calves and all eight isolates were susceptible to all tested antimicrobials. Most Campylobacter isolates were susceptible to all tested antimicrobials, although 16.6% of C. jejuni were resistant to quinolones and tetracycline. Results suggest that AMR in E. coli, enterococci, Salmonella, and Campylobacter isolates from food animals in New Zealand is low, and currently, AMR in food animals poses a limited public health risk. Despite the low prevalence of AMR in this survey, ongoing monitoring of antimicrobial susceptibility in bacteria from food animals is recommended, to ensure timely detection of AMR with potential impacts on animal and human health.

Prevention of legionella infections from toilet flushing cisterns

INTRODUCTION

Immunocompromised patients are at an increased risk of severe legionella infections. We present the results of an outbreak investigation initiated following a fatal case of hospital-acquired legionellosis linked to contaminated water from a toilet-flushing cistern. Additionally, we provide experimental data on the growth of Legionella spp. in flushing cisterns and propose a straightforward protocol for prevention.

METHODS

We monitored the growth of Legionella spp. in the building's hot- and cold-water systems using quantitative bacterial culture on selective agar. Molecular typing of Legionella pneumophila isolates from the infected patient and the water system was conducted through core-genome multi-locus sequence typing (cgMLST).

RESULTS

Legionella contamination in the hospital building's cold-water system was significantly higher than in the hot-water system and significantly higher in toilet flushing cistern's water compared with cold water from bathroom sinks and showers. Isolates from the patient and from the flushing cistern of the patient's bathroom were identical by cgMLST. In an experimental setting, daily toilet flushing for a period of 21 days resulted in a 67% reduction in the growth of Legionella spp. in the water of toilet flushing cisterns. Moreover, a one-time disinfection of cisterns with peracetic acid, followed by daily flushing, decreased legionella growth to less than 1% over a period of at least seven weeks in these setting.

CONCLUSIONS

One-time disinfection of highly contaminated cisterns with peracetic acid and daily toilet flushing as short-term measure can significantly reduce legionella contamination in flushing cisterns. These measures may aid in preventing legionella infection among immunocompromised patients.

Emergence of OXA-48-producing Klebsiella pneumoniae in Lithuania, 2023: a multi-cluster, multi-hospital outbreak

In 2023, an increase of OXA-48-producing Klebsiella pneumoniae was noticed by the Lithuanian National Public Health Surveillance Laboratory. Whole genome sequencing (WGS) of 106 OXA-48-producing K. pneumoniae isolates revealed three distinct clusters of carbapenemase-producing K. pneumoniae high-risk clones, including sequence type (ST) 45 (n = 35 isolates), ST392 (n = 32) and ST395 (n = 28), involving six, six and nine hospitals in different regions, respectively. These results enabled targeted investigation and control, and underscore the value of national WGS-based surveillance for antimicrobial resistance.

Analysis of Emerging Variants of Turkey Reovirus using Machine Learning

Avian reoviruses continue to cause disease in turkeys with varied pathogenicity and tissue tropism. Turkey enteric reovirus has been identified as a causative agent of enteritis or inapparent infections in turkeys. The new emerging variants of turkey reovirus, tentatively named turkey arthritis reovirus (TARV) and turkey hepatitis reovirus (THRV), are linked to tenosynovitis/arthritis and hepatitis, respectively. Turkey arthritis and hepatitis reoviruses are causing significant economic losses to the turkey industry. These infections can lead to poor weight gain, uneven growth, poor feed conversion, increased morbidity and mortality and reduced marketability of commercial turkeys. To combat these issues, detecting and classifying the types of reoviruses in turkey populations is essential. This research aims to employ clustering methods, specifically K-means and Hierarchical clustering, to differentiate three types of turkey reoviruses and identify novel emerging variants. Additionally, it focuses on classifying variants of turkey reoviruses by leveraging various machine learning algorithms such as Support Vector Machines, Naive Bayes, Random Forest, Decision Tree, and deep learning algorithms, including convolutional neural networks (CNNs). The experiments use real turkey reovirus sequence data, allowing for robust analysis and evaluation of the proposed methods. The results indicate that machine learning methods achieve an average accuracy of 92%, F1-Macro of 93% and F1-Weighted of 92% scores in classifying reovirus types. In contrast, the CNN model demonstrates an average accuracy of 85%, F1-Macro of 71% and F1-Weighted of 84% scores in the same classification task. The superior performance of the machine learning classifiers provides valuable insights into reovirus evolution and mutation, aiding in detecting emerging variants of pathogenic TARVs and THRVs.

Antimicrobial resistance and heterogeneity of Neisseria gonorrhoeae isolated from patients attending sexually transmitted infection clinics in Lusaka, Zambia

BACKGROUND

Antimicrobial resistance (AMR) of Neisseria gonorrhoeae is a threat to public health as strains have developed resistance to antimicrobials available for the treatment of gonorrhea. Whole genome sequencing (WGS) can detect and predict antimicrobial resistance to enhance the control and prevention of gonorrhea. Data on the molecular epidemiology of N. gonorrhoeae is sparse in Zambia. This study aimed to determine the genetic diversity of N. gonorrhoeae isolated from patients attending sexually transmitted infection (STI) clinics in Lusaka, Zambia.

METHODS

A cross-sectional study that sequenced 38 N. gonorrhoeae isolated from 122 patients with gonorrhea from 2019 to 2020 was conducted. The AMR profiles were determined by the E-test, and the DNA was extracted using the NucliSens easyMaG magnetic device. Whole genome sequencing was performed on the Illumina NextSeq550 platform. The Bacterial analysis pipeline (BAP) that is readily available at: https://cge.cbs.dtu.dk/services/CGEpipeline-1.1 was used for the identification of the species, assembling the genome, multi-locus sequence typing (MLST), detection of plasmids and AMR genes. Phylogeny by single nucleotide polymorphisms (SNPs) was determined with the CCphylo dataset.

RESULTS

The most frequent STs with 18.4% of isolates each were ST7363, ST1921 and ST1582, followed by ST1583 (13%), novel ST17026 (7.9%), ST1588 (7.9%), ST1596 (5.3%), ST11181 (5.3%), ST11750 (2.6/%) and ST11241 (2.6%) among the 38 genotyped isolates. The blaTeM-1B and tetM (55%) was the most prevalent combination of AMR genes, followed by blaTeM-1B (18.4%), tetM (15.8%), and the combination of blaTeM-1B, ermT, and tetL was 2.6% of the isolates. The AMR phenotypes were predicted in ciprofloxacin, penicillin, tetracycline, azithromycin, and cefixime. The combination of mutations 23.7% was gryA (S91F), parC (E91G), ponA (L421) and rpsJ (V57M), followed by 18.4% in gyrA (S91F), ponA (L421P), rpsJ (V57M), and 18.4% in gyrA (D95G, S91F), ponA (L421P), and rpsJ (V57M). The combinations in gyrA (D95G, S91F) and rpsJ (V57M), and gyrA (D95G, S91F), parC (E91F), ponA (L421P) and rpsJ (V57M) were 13.2% each of the isolates. Plasmid TEM-1 (84.2%), tetM (15.8%), and gonococcal genetic island (GGI) was detected in all isolates.

CONCLUSION

This study revealed remarkable heterogeneity of N. gonorrhoeae with blaTEM-1, tetM, ponA, gyrA, and parC genes associated with high resistance to penicillin, tetracycline, and ciprofloxacin demanding revision of the standard treatment guidelines and improved antimicrobial stewardship in Zambia.

Draft genome sequences of Escherichia spp. isolates from New Zealand environmental sources

Escherichia coli is often used as a fecal indicator bacterium for water quality monitoring. We report the draft genome sequences of 500 Escherichia isolates including newly described Escherichia species, namely Escherichia marmotae, Escherichia ruysiae, and Escherichia whittamii, obtained from diverse environmental sources to assist with improved public health risk assessments.

Chiral Nanoclusters as Alternative Therapeutic Strategies to Confront the Health Threat from Antibiotic-Resistant Pathogens

Pseudomonas aeruginosa (P. aeruginosa), a drug-resistant Gram-negative pathogen, is listed among the "critical" group of pathogens by the World Health Organization urgently needing efficacious antibiotics in the clinics. Nanomaterials especially silver nanoparticles (AgNPs) due to the broad-spectrum antimicrobial activity are tested in antimicrobial therapeutic applications. Pathogens rapidly develop resistance to AgNPs; however, the health threat from antibiotic-resistant pathogens remains challenging. Here we present a strategy to prevent bacterial resistance to silver nanomaterials through imparting chirality to silver nanoclusters (AgNCs). Nonchiral AgNCs with high efficacy against P. aeruginosa causes heritable resistance, as indicated by a 5.4-fold increase in the minimum inhibitory concentration (MIC) after 9 repeated passages. Whole-genome sequencing identifies a Rhs mutation related to the wall of Gram-negative bacteria that possibly causes morphology changes in resistance compared to susceptible P. aeruginosa. Nevertheless, AgNCs with laevorotary chirality (l-AgNCs) induce negligible resistance even after 40 repeated passages and maintain a superior antibacterial efficiency at the MIC. l-AgNCs also show high cytocompatibility; negligible cytotoxicity to mammalian cells including JB6, H460, HEK293, and RAW264.7 is observed even at 30-fold MIC. l-AgNCs thus are examined as an alternative to levofloxacin in vivo, healing wound infections of P. aeruginosa efficaciously. This work provides a potential opportunity to confront the rising threat of antimicrobial resistance by developing chiral nanoclusters.

VirulenceFinder for Enterococcus faecium and Enterococcus lactis: an enhanced database for detection of putative virulence markers by using whole-genome sequencing data

Enterococcus faecium (Efm) is a leading cause of hospital-associated (HA) infections, often enriched in putative virulence markers (PVMs). Recently, the Efm clade B was assigned as Enterococcus lactis (Elts), which usually lack HA-Efm infection markers. Available databases for extracting PVM are incomplete and/or present an intermix of genes from Efm and Enterococcus faecalis, with distinct virulence profiles. In this study, we constructed a new database containing 27 PVMs [acm, scm, sgrA, ecbA, fnm, sagA, hylEfm, ptsD, orf1481, fms15, fms21-fms20 (pili gene cluster 1, PGC-1), fms14-fms17-fms13 (PGC-2), empA-empB-empC (PGC-3), fms11-fms19-fms16 (PGC-4), ccpA, bepA, gls20-glsB1, and gls33-glsB] from nine reference genomes (seven Efm + two Elts). The database was validated against these reference genomes and further evaluated using a collection of well-characterized Efm (n = 43) and Elts (n = 7) control strains, by assessing PVM presence/absence and its variants together with a genomic phylogeny constructed as single-nucleotide polymorphisms. We found a high concordance between the phylogeny and in silico findings of the PVM, with Elts clustering separately and mostly carrying Elts-specific PVM gene variants. Based on our validation results, we recommend using the database with raw reads instead of assemblies to avoid missing gene variants. This newly constructed database of 27 PVMs will enable a more comprehensive characterization of Efm and Elts based on WGS data. The developed database exhibits scalability and boasts a range of applications in public health, including diagnostics, outbreak investigations, and epidemiological studies. It can be further used in risk assessment for distinguishing between safe and unsafe enterococci.IMPORTANCEThe newly constructed database, consisting of 27 putative virulence markers, is highly scalable and serves as a valuable resource for the comprehensive characterization of these closely related species using WGS data. It holds significant potential for various public health applications, including hospital outbreak investigations, surveillance, and risk assessment for probiotics and feed additives.

Genomic epidemiology of hypervirulent Listeria monocytogenes CC619: Population structure, phylodynamics and virulence

Listeria monocytogenes is a ubiquitous foodborne pathogen causing human and animal listeriosis with high mortality. Neurological and maternal-neonatal listeriosis outbreaks in Europe and the United States were frequently associated with clonal complexes CC1, CC2 and CC6 harboring Listeria Pathogenicity Island-1 (LIPI-1), as well as CC4 carrying both LIPI-1 and LIPI-4. However, human listeriosis in China was predominantly linked to CC87 and CC619 from serotype 1/2b. To understand the genetic evolution and distribution patterns of CC619, we characterized the epidemic history, population structure, and transmission feature of CC619 strains through analysis of 49,421 L. monocytogenes genomes globally. We found that CC619 was uniquely distributed in China, and closely related with perinatal infection. As CC619 strains were being mainly isolated from livestock and poultry products, we hypothesized that pigs and live chicken were the reservoirs of CC619. Importantly, all CC619 strains not only harbored the intact LIPI-1 and LIPI-4, but these also carried LIPI-3 that could facilitate host colonization and invasion. The deficiency of LIPI-3 or LIPI-4 markedly decreased L. monocytogenes colonization capacity in a model of intragastric infection in the mouse. Altogether, our findings suggest that the hypervirulent CC619 harboring three pathogenicity islands LIPI-1, LIPI-3 and LIPI-4 is a putatively persistent population in various foods, environment, and human population, warranting the further research for deciphering its pathogenicity and strengthening epidemiological surveillance.