ConFindr: rapid detection of intraspecies and cross-species contamination in bacterial whole-genome sequence data

Molecular characterization of multidrug-resistant E. coli recovered from diarrheagenic children under 5 years from Mukuru Informal Settlement, Nairobi, Kenya, based on whole-genome sequencing analysis

High genomic plasticity within Escherichia coli enables it to acquire and accumulate genetic material through horizontal gene transfer. In this study, we sought to investigate the virulence genes, phylogroups, antibiotic resistance genes, plasmid replicons, multilocus sequence types (MLST), and core genome MLST of multidrug-resistant E. coli recovered from diarrheagenic children under 5 years from Mukuru Informal Settlement in Nairobi, Kenya. A total of 39 multidrug-resistant (MDR) strains had their DNA extracted, and whole-genome sequencing was done using the Illumina HiSeq 2000 platform. Twenty-six E. coli assemblies were analyzed using web-based bioinformatics tools available at the Centre for Genomic Epidemiology and EnteroBase. The isolates were categorized into four main phylogroups, where 10/26 (38.5%) belonged to the B2 phylogroup, 4/26 (15.4%) belonged to D, 3/26 (11.5%) belonged to A, 1/26 (3.8%) belonged to B1, while 8/26 (30.8%) were not determined. FimH30 was predominantly found in the most frequent phylogroup B2 and sequence type (ST) 131. The most common beta-lactam resistance genes were bla TEM-1B and blaCTXM 15, followed by three fluoroquinolone resistance genes [qnrS1 6/26 (23.1%), qnrB4 2/26 (7.7%), and aac(6')-Ib-cr, 8/26 (30.8%)]. Of 26 isolates, 15 had at least one amino acid substitution in the housekeeping genes gyrA (p.S83L), gyrA (p.D87N), parC (p.S80I), parC (p.E84V), parC (p.S57T), and parE (p.I529L), associated with resistance to fluoroquinolones. A total of 40 diverse virulence genes were detected among the isolates. Thirteen different STs were isolated from the E. coli genomes, which included ST 131, ST 3036, ST 38, ST 10, ST 12569, ST 15271, ST 2076, ST 311, ST 3572, ST 394, ST 453, ST 46, and ST 1722. Only two isolates (2/26, 7.7%) from the Municipal City Council clinic were genetically related. Additionally, the most abundant plasmid replicon identified belonged to the IncF family, IncFII(pRSB107), in particular, followed by the Col family. The study highlighted the first E. coli ST46 to harbor the bla NDM5 gene encoded in Col(BS512), IncFII(pRSB107), and IncFIB(AP001918) plasmid replicons in Kenya. We further demonstrated the diversity of MDR E. coli associated with diarrhea in an endemic setting in Kenya.

IMPORTANCE

This study investigated the molecular characterization of multidrug-resistant Escherichia coli isolated from diarrheagenic children under 5 years of age in Mukuru Informal Settlement in Nairobi, Kenya. This is an important addition to the genomic analysis data of multi-drug resistant diarrheal Escherichia coli in Kenya. The use of whole-genome sequencing to identify and characterize these isolates is valuable and provides valuable insights into the molecular epidemiology of E. coli in the region.

Whole-genome sequencing of soil- and foodborne Bacillus cereus sensu lato indicates no clear association between their virulence repertoire, genomic diversity and food matrix

Bacillus cereus sensu lato is frequently involved in foodborne toxico-infections and is found in various foodstuff. It is unclear whether certain strains have a higher affinity for specific food matrices, which can be of interest for risk assessment. This study reports the characterization by whole-genome sequencing of 169 B. cereus isolates, isolated from 12 food types and soil over two decades. Any potential links between the food matrix of isolation, the isolate's genetic lineage and/or their (putative) virulence gene reservoir were investigated. More than 20 % of the strains contained the genes for the main potential enterotoxins (nheABC, hblCDA and cytK_2). Cereulide biosynthesis genes and genes encoding hemolysins and phospholipases, were detected in multiple isolates. Strain typing revealed a high diversity, as illustrated by 84 distinct sequence types, including 26 not previously described. This diversity was also reflected in the detection of all seven panC types and 71 unique virulence gene profiles. Core-genome MLST was used for phylogenomic investigation of the entire collection and SNP-based clustering was performed on the four most abundant sequence types, which did not reveal a clear affinity for specific B. cereus lineages or (putative) virulence genes for certain food matrices. Additionally, minimal genetic overlap was observed between soil and foodborne isolates. Clusters of closely-related isolates with common epidemiological metadata were detected. However, some isolates from different food matrices or collected several years apart were found to be genetically identical. This study provides elements that can be used for risk assessment of B. cereus in food.

Multi-country and intersectoral assessment of cluster congruence between pipelines for genomics surveillance of foodborne pathogens

Different laboratories employ different Whole-Genome Sequencing (WGS) pipelines for Food and Waterborne disease (FWD) surveillance, casting doubt on the comparability of their results and hindering optimal communication at intersectoral and international levels. Through a collaborative effort involving eleven European institutes spanning the food, animal, and human health sectors, we aimed to assess the inter-pipeline clustering congruence across all resolution levels and perform an in-depth comparative analysis of cluster composition at outbreak level for four important foodborne pathogens: Listeria monocytogenes, Salmonella enterica, Escherichia coli, and Campylobacter jejuni. We found a general concordance between allele-based pipelines for all species, except for C. jejuni, where the different resolution power of allele-based schemas led to marked discrepancies. Still, we identified non-negligible differences in outbreak detection and demonstrated how a threshold flexibilization favors the detection of similar outbreak signals by different laboratories. These results, together with the observation that different traditional typing groups (e.g., serotypes) exhibit a remarkably different genetic diversity, represent valuable information for future outbreak case-definitions and WGS-based nomenclature design. This study reinforces the need, while demonstrating the feasibility, of conducting continuous pipeline comparability assessments, and opens good perspectives for a smoother international and intersectoral cooperation towards an efficient One Health FWD surveillance.

Prevalence of Shiga Toxin-Producing Escherichia coli (STEC) and Risk Characterization Based on Virulence Genes in Retail Raw Ground Meat of Beef, Veal, and Lamb in Canada

Shiga toxin-producing Escherichia coli (STEC) are potentially pathogenic E. coli that may cause mild to severe gastrointestinal illnesses. STEC-contaminated foods of animal origin have been the most frequently implicated sources of foodborne outbreaks. A multiyear (2016 to 2021) targeted survey was conducted to investigate the prevalence of STEC in retail ground meats (beef, veal, and lamb). Samples were screened for the presence of Shiga toxin genes (stx) to identify presumptive samples, followed by culture and molecular confirmation of isolates to confirm the presence of stx genes and subsequent characterization by whole-genome-sequencing (WGS) for O serogroup and virulence genes (e.g., stx, eae, aggR). A total of 175 STEC strains were isolated from a total of 148 samples where the presence of viable STEC was confirmed out of 2,398 ground meat samples. This represented 1.2% (7/589 positive, 8 unique strains) of the beef samples, 4.7% (58/1,241 positive, 67 unique strains) of the veal samples, and 14.6% (83/568 positive, 100 unique strains) of the lamb samples. The intimin virulence gene, eae, was identified in the STEC strains of veal origin (9/67, 13.4%) only and were classified as belonging to risk level 1 (1/67), level 3 (2/67), and level 4 (6/67) according to the FAO/WHO risk categories. Risk level 2 STEC strains were of beef (2/8, 25.0%), veal (8/67, 11.9%), and lamb (1/100, 1.0%) origin. The majority of the STEC strains, 75.0% (6/8) of the beef, 67.2% (45/67) of the veal, and 94.0% (94/100) of the lamb STEC strains were classified as risk level 5 (lowest level) of the FAO/WHO risk categories. This study's findings indicate that the current food safety control measures implemented for ground meats in Canada are effective at maintaining an acceptable level of possible contamination with STEC strains associated with severe clinical outcomes. Continued application of effective control measures and safe food handling practices by consumers will minimize the potential risk of foodborne infections from raw ground meats.

Profiling the resistome and virulome of Bacillus strains used for probiotic-based sanitation: a multicenter WGS analysis

BACKGROUND

Healthcare-associated infections (HAIs) caused by microbes that acquire antimicrobial resistance (AMR) represent an increasing threat to human health worldwide. The high use of chemical disinfectants aimed at reducing the presence of pathogens in the hospital environment can simultaneously favor the selection of resistant strains, potentially worsening AMR concerns. In the search for sustainable ways to control bioburden without affecting this aspect, probiotic-based sanitation (PBS) using Bacillus spp. was proposed to achieve stable reduction of pathogens, AMR, and associated HAIs. Although Bacillus probiotics are classified as nonpathogenic, comprehensive data about the potential genetic alterations of these probiotics following prolonged contact with surrounding pathogens are not yet available. This study aimed to assess in depth the genetic content of PBS-Bacillus isolates to evaluate any eventual variations that occurred during their usage.

RESULTS

WGS analysis was used for the precise identification of PBS-Bacillus species and detailed profiling of their SNPs, resistome, virulome, and mobilome. Analyses were conducted on both the original PBS detergent and 172 environmental isolates from eight hospitals sanitized with PBS over a 30-month period. The two species B. subtilis and B. velezensis were identified in both the original product and the hospital environment, and SNP analysis revealed the presence of two clusters in each species. No virulence/resistance genes or mobile conjugative plasmids were detected in either the original PBS-Bacillus strain or any of the analyzed environmental isolates, confirming their high genetic stability and their low/no tendency to be involved in horizontal gene transfer events.

CONCLUSIONS

The data obtained by metagenomic analysis revealed the absence of genetic sequences associated with PBS-Bacillus and the lack of alterations in all the environmental isolates analyzed, despite their continuous contact with surrounding pathogens. These results support the safety of the Bacillus species analyzed. Further metagenomic studies aimed at profiling the whole genomes of these and other species of Bacillus, possibly during longer periods and under stress conditions, would be of interest since they may provide further confirmation of their stability and safety.

Comparative genomics of Mycobacterium avium subsp. hominissuis strains within a group of captive lowland tapirs

Within a group of three captive lowland tapirs (Tapirus terrestris) suffering from clinically apparent mycobacteriosis, non-tuberculous Mycobacterium avium subsp. hominissuis (MAH) strains were isolated from the animals and the tapir's enclosure. Based on MIRU-VNTR findings, which identified two closely related INMV profiles (124 and 246), a micro-evolutionary event was assumed, and four available MAH strains were submitted to whole genome sequencing (short- and long-read technologies). Surprisingly, the differences based on single nucleotide polymorphisms (SNPs) were exceptionally high between the four strains, i.e., between 841 and 11,166 bases, due to a strong impact of homologous recombination. Thus, an ad hoc core genome multilocus sequence typing (cgMLST) scheme was created and pangenome analysis was conducted for determining the genomic similarity between the strains. The INMV246 isolate obtained from sputum on the enclosure floor and one INMV124 isolate of tapir #2 showed the highest congruence, suggesting that both originated from a shared source. The other two INMV124 isolates were genomically distinct from these strains. Nevertheless, in all four strains two plasmids were detected, which were highly conserved between the strains. The study showed that the genomic variability between MAH strains isolated from the same site within a short period of time can be exceptionally high and the influence of homologous recombination needs to be considered when determining MAH strain relationships, particularly via SNP analyses.

Comprehensive molecular epidemiology of Acinetobacter baumannii from diverse sources in Nigeria

BACKGROUND

Acinetobacter baumannii, a Gram-negative bacterium, is a public health threat due to its role in nosocomial infections and increasing antibiotic resistance. In Nigeria, data on the molecular epidemiology of A. baumannii is scarce. This study investigates the genetic diversity and the presence of antimicrobial resistance determinants and virulence-related genes in whole-genome sequencing data of 189 Nigerian A. baumannii isolates deposited in public repositories. Genotypes were determined in-silico by multilocus sequence typing (MLST) and core genome MLST (cgMLST). Further, antimicrobial resistance (AMR) and virulence-related genes were analyzed.

RESULTS

Most isolates (57.67%) originated from South-west Nigeria. Isolates of human origin accounted for 33.86%, while environmental sources comprised 6.87%, and 59.27% lacked information on the source of isolation. The cgMLST analysis revealed a multitude of genomic lineages circulating in Nigeria. The MLST Oxford scheme identified 44 sequence types (STs) in 62.96% of strains, with ST1089 being the most prevalent. The MLST Pasteur could assign 95.77% of strains to 49 STs, with ST2(IC2) and ST85(IC9) being the most dominant. Antimicrobial resistance analysis detected 168 genes encoding resistance to 12 antibiotic classes, with cephalosporin, carbapenem, and aminoglycoside resistance genes being the most prevalent. Notably, blaADC-79 (23.81%), blaOXA-23 (30.69%), and aph(3″)-Ib (30%) were frequent variants encountered. Seventeen multi-efflux system genes conferring resistance to multiple antibiotic classes were identified. Virulence gene analysis revealed 137 genes encoding six mechanisms, with genes for nutritional factors, effector delivery systems, and biofilm production being the most prevalent.

CONCLUSION

This study highlights the diversity in AMR and virulence genes of A. baumannii in Nigeria, emphasizing the need for ongoing genomic surveillance to inform infection control and develop antibiotic resistance management strategies.

Core genome multilocus sequence typing schemes for epidemiological investigation of Taylorella equigenitalis and Taylorella asinigenitalis

Taylorella equigenitalis is the causative agent of contagious equine metritis, an internationally regulated sexually-transmitted infection in horses, which is of great concern as it usually results in temporary infertility. Taylorella asinigenitalis, the second member of the genus, is mainly found in donkeys and is considered non-pathogenic, although a first natural outbreak was reported in 2019 in the United Arab Emirates. Multilocus sequence typing (MLST) is currently used to study the epidemiology of Taylorella spp. but, while highly transposable and reproducible, it only focuses on < 0.5 % of the genome (seven genes). We therefore aimed to develop a robust core genome MLST (cgMLST) based on the analysis of 370 T. equigenitalis and 68 T. asinigenitalis genomes belonging to 46 and 18 sequence types (STs), respectively. Typing results based on 1333 loci (84.0 % of the genome) from T. equigenitalis genomes and 1255 loci (80.3 %) from T. asinigenitalis genomes showed that the discriminatory power of both species-specific cgMLSTs was greater than that of MLST, with 368 and 68 distinct core genome STs (cgSTs), respectively. Clustering was congruent between the cgMLST and MLST methods, with few inconsistencies for T. equigenitalis. Maximum allelic distance between epidemiologically-related strains was used to define cgMLST clustering thresholds, set at ≤ 15 and 20 allelic distances for T. equigenitalis and T. asinigenitalis, respectively. These parameters grouped the cgSTs into 47 and 11 clonal groups (CGs), respectively. Overall, the cgMLST method outperformed conventional MLST in distinguishing clonal strains from epidemiologically-linked strains, supporting the hypothesis that typing based on a few housekeeping genes does not always accurately reflect genomic relatedness between strains, and making cgMLST more suitable for outbreak investigations.

Galaxy @Sciensano: a comprehensive bioinformatics portal for genomics-based microbial typing, characterization, and outbreak detection

The influx of whole genome sequencing (WGS) data in the public health and clinical diagnostic sectors has created a need for data analysis methods and bioinformatics expertise, which can be a bottleneck for many laboratories. At Sciensano, the Belgian national public health institute, an intuitive and user-friendly bioinformatics tool portal was implemented using Galaxy, an open-source platform for data analysis and workflow creation. The Galaxy @Sciensano instance is available to both internal and external scientists and offers a wide range of tools provided by the community, complemented by over 50 custom tools and pipelines developed in-house. The tool selection is currently focused primarily on the analysis of WGS data generated using Illumina sequencing for microbial pathogen typing, characterization and outbreak detection, but it also addresses specific use cases for other data types. Our Galaxy instance includes several custom-developed 'push-button' pipelines, which are user-friendly and intuitive stand-alone tools that perform complete characterization of bacterial isolates based on WGS data and generate interactive HTML output reports with key findings. These pipelines include quality control, de novo assembly, sequence typing, antimicrobial resistance prediction and several relevant species-specific assays. They are tailored for pathogens with active genomic surveillance programs, and clinical relevance, such as Escherichia coli, Listeria monocytogenes, Salmonella spp. and Mycobacterium tuberculosis. These tools and pipelines utilize internationally recognized databases such as PubMLST, EnteroBase, and the NCBI National Database of Antibiotic Resistant Organisms, which are automatically synchronized on a regular basis to ensure up-to-date results. Many of these pipelines are part of the routine activities of Belgian national reference centers and laboratories, some of which use them under ISO accreditation. This resource is publicly available for noncommercial use at https://galaxy.sciensano.be/ and can help other laboratories establish reliable, traceable and reproducible bioinformatics analyses for pathogens encountered in public health settings.

Development and evaluation of a sensitive approach for detection and recovery of third-generation cephalosporin- and carbapenem-resistant Enterobacterales from ready-to-eat frozen stone fruit

Antimicrobial resistance (AMR) is a global public health threat, but the role of foods in its dissemination is poorly understood. We examined the incidence of foodborne bacteria carrying AMR genes considered high-priority research targets by the World Health Organization. Frozen, ready-to-eat, avocado, coconut, mango, and peach (n = 161) were tested for bacteria encoding extended-spectrum β-lactamases (ESBLs) and carbapenemases. Over 600 presumptive-positive isolates were recovered and analyzed with a pooled sequencing (Pool-seq) strategy. Coconut samples exhibited the highest bacterial loads and prevalence/diversity of AMR genes. Isolates harbouring the β-lactamase genes blactx-m, blatem, and blashv, identified in 14 coconut and 2 mango samples, were further characterized by whole-genome sequencing and antimicrobial susceptibility testing. The most common gene was blactx-m-15, detected in 20 unique strains. Two carbapenemase-producing strains were isolated from coconut: Enterobacter roggenkampii encoding blandm-1 and Escherichia coli encoding blandm-5. Subsequent quantitative PCR (qPCR) analysis of enrichments for blactx-m/blandm indicated a potentially higher prevalence of these genes than observed by colony screening. This study presents a practical method for recovering ESBL- and carbapenemase-producing bacteria from foods. Mapping their distribution in food products is crucial to assessing the role of foods in the global spread of AMR and developing effective public health interventions.

Epidemiology of Staphylococcus aureus food isolates: Comparison of conventional methods with whole genome sequencing typing methods

A variety of methods exists for typing bacteria. However, guidelines for the application and interpretation of typing tools in epidemiologic investigations of Staphylococcus aureus are lacking. This study aimed to identify appropriate typing methods for S. aureus population studies and outbreak investigation. We compared pulsed-field gel electrophoresis (PFGE), seven loci multi-locus sequence typing (MLST), core genome MLST (cgMLST), core single nucleotide polymorphism (cSNP), and enterotoxin (se/SE) profiles on 351 S. aureus isolates. The discriminatory power, concordance, and congruence of typing results were assessed. cgMLST, cSNP, and PFGE yielded the highest discrimination value, followed by se/SE typing and MLST. The best concordance of results was found between cgMLST and cSNP, while the best congruence was observed for cgMLST and cSNP with all methods, followed by PFGE with MLST. The strengths and weaknesses of each method are highlighted. For population structure, cgMLST and cSNP performed better than PFGE and MLST in terms of resolution of clusters and in phylogenetic inference. Enterotoxin profiles matched with MLST groups, suggesting the use of se/SE typing to predict MLST results. For the retrospective analysis of 31 outbreaks, all methods performed almost equally to discriminate epidemiologically related strains and can be used to unambiguously distinguish outbreak strains.

Genomic characterization of antimicrobial-resistance and virulence factors in Salmonella isolates obtained from pig farms in Antioquia, Colombia

BACKGROUND

Occurrence of antimicrobial-resistant Salmonella strains has been reported worldwide, because of inappropriate use of antimicrobial products in either humans or animals. The presence of multidrug resistant Salmonella in pig production systems had been reported in Antioquia, Colombia.

AIM

To identify antimicrobial resistance genes (ARG) in different Salmonella spp. strains isolated from pig productions in Antioquia, Colombia. Methods: Samples were received at the Diagnostic Unit of the Faculty of Agrarian Sciences at the University of Antioquia, from January 1, 2019, to January 2021. A total of 28 isolates of Salmonella spp. were included, which presented phenotypic resistance to more than one antibiotic used in pig farms. Whole genome sequencing (WGS) was performed in the Unit of Genomic of Agrosavia using an automated pipeline from the GHRU- Sanger Institute, employing the Illumina MiSeq platform.

RESULTS

WGS revealed 34 ARGs among these isolates. In 25 isolates (89%) more than two ARGs were found. Genes encoding resistance were found for 10 different groups of antibiotics (beta-lactam, aminoglycosides, chloramphenicol, rifampicins, lincosamides, fluoroquinolones, tetracyclines, sulfonamides and trimethoprim). The most frequently observed MDR profile in Typhimurium isolates was AMP-CEX-CEP-CEF-EFT-CEQ-FLU-ENR-TE-FFC-SXT.

CONCLUSION

The presence of multi-drug resistant Salmonella strains in pigs destined for human consumption in Antioquia, Colombia was determined. This research emphasizes the utmost importance of epidemiological tools to understand the presence and spreading of antimicrobial resistance genes in pig farms. Additionally, it highlights the critical need for developing educational programs and public policies to help reduce the spread of antimicrobial resistance in production systems.

Genomic diversity of Salmonella Typhimurium and its monophasic variant in pig and pork production in France

Salmonella Typhimurium and its monophasic variant (Salmonella 4,[5],12:i:-) are among the most prevalent serovars worldwide. Even though these serovars have been the focus of many studies, their spread has not yet been investigated in French pig herds and slaughterhouses at a regional scale. Here, we characterized the genomic diversity of 188 Salmonella strains belonging to sequence type (ST) 19 and 34. These strains were isolated from pigs in metropolitan France between 2014 and 2019. Samples were collected from 10 regions, three of which together represent 75% of French pig production in 2020. To contextualize the French Salmonella genomes at a worldwide level, 193 ST 34 genomes from three continents and 14 countries were also included. This study revealed little diversity in ST 34 strains circulating in France, suggesting that one or two clones are spreading within pig herds and slaughterhouses. In silico virulence and antimicrobial resistance genes were investigated to understand the prevalence of these strains among farmed pigs and in the slaughterhouse environment. A comparison with ST 34 isolates from other countries highlighted the genomic specificity of the ST 34 monophasic variants in France, with some exceptions concerning isolates from bordering countries. This work provides new insights into the dynamics of S. Typhimurium and its monophasic variant sampled in French pig herds and slaughterhouses.

IMPORTANCE

Salmonellosis is a leading cause of bacterial infection in humans and animals around the world. This study provides a snapshot of the genomic diversity of one of the most prevalent Salmonella serovars (Salmonella Typhimurium and Salmonella 4,[5],12:i:-) circulating on French pig farms between 2013 and 2021. We investigated the link between geographical and genomic diversity. The analyses revealed little diversity of the strains, suggesting that one or two clones are spreading within French pig herds. We also in silico screened genetic elements that could explain the prevalence of these strains among farmed pigs and in the slaughterhouse environment. Finally, the comparison with isolates from other countries highlighted the genomic specificity of these two French sequence type 34 clones. This work provides new insights into the dynamics of S. Typhimurium and S. 4,[5],12:i:- sampled from pig herds and slaughterhouses in France, thus laying the foundations for future analyses.

Collagen binding adhesin restricts Staphylococcus aureus skin infection

Staphylococcus aureus causes approximately 80% of skin and soft tissue infections (SSTIs). Collagen is the most abundant human extracellular matrix protein with critical roles in wound healing, and S. aureus encodes a collagen binding adhesin (Cna). The role of this protein during skin infections is unknown. Here we report that inability to bind collagen results in worsened pathology of intradermal Δcna S. aureus infection. WT/Cna+ S. aureus showed reduced infection severity, aggregate formation, and significantly improved clearance of bacteria. Cna binds to the collagen-like domain of serum C1q protein to reduce its opsonophagocytic functions. We demonstrate that infection of C1qKO mice with WT bacteria show results similar to the Δcna group. Conversely, inability to bind collagen resulted in an amplified inflammatory response caused in part by macrophage and neutrophil small molecule mediators released at the infection site (MMP-9, MMP-12, LTB4), resulting in increased immune cell infiltration and death.

Molecular Characterization of Pseudomonas aeruginosa Clinical Isolates Through Whole-Genome Sequencing: A Comprehensive Analysis of Biotypes, Sequence Types, and Antimicrobial and Virulence Genes

Introduction Pseudomonas aeruginosa (PA) is a bacterial species commonly isolated from human clinical specimens. Despite being present in the environment as a saprophyte, PA possesses the ability to cause human infections, especially among debilitated patients. It is therefore essential to understand the genomic imprints of antimicrobial resistance (AMR) and virulence genes associated with PA isolated from patient samples. Methods The study carried out next-generation sequencing (NGS) or whole-genome sequencing (WGS) of nine PA strains isolated from various human clinical specimens from patients at Prathima Institute of Medical Sciences, Karimnagar, India. All the isolates were identified by conventional microbiological methods and confirmed by automated systems. Antimicrobial susceptibility patterns of the isolates were carried out using the Kirby-Bauer disc diffusion method. Additionally, NGS/WGS was done to evaluate the carriage of AMR and virulence genes associated with each PA strain. Sequence type was identified through multi-locus sequence typing (MLST). Results The genotype and phenotypic antimicrobial susceptibility patterns revealed the same (11.11% resistance) results with carbapenems and fluoroquinolone antibiotics. However, discordant antimicrobial susceptibility patterns were noticed with trimethoprim-sulfamethoxazole (66.66% resistance phenotype vs. 100% sensitive genotype), aminoglycosides (100% sensitive phenotype vs. 100% resistant genotype), and beta-lactamase/extended-spectrum beta-lactamase (ESBL) (44.44% sensitive phenotype vs. 100% resistant genotype) antibiotics. All (100%, 9/9) the PA isolates included in the study demonstrated the presence of multiple antibiotic resistance and virulence genes. The antibiotic resistance genes identified included aph, aad, aac, bla PDC, bla OXA, bla VIM, catB7, fosA, qnrVC1, and crpP. All (100%, 9/9) isolates demonstrated the presence of class C beta-lactamase bla PDC and class B metallo-beta-lactamase bla OXA. Only one (11.11%, 1/9) isolate showed the presence of subclass B1 metallo-beta-lactamase bla VIM. Among the virulence genes identified were toxA, fih, xcp, wzz, pvc, pvd, and many others. This study showed the presence of ST244, a high-risk PA strain with global significance. Conclusions PA is an opportunistic pathogen, and its isolation among hospitalized patients should be carefully evaluated. Tracking PA for the presence of high-risk sequence types and the prevalence of resistance and virulence genes could improve the understanding of the organism. Molecular data obtained from this study demonstrated that the PA isolates carried multiple antibiotic-resistant and virulence genes that could potentially enable them to cause invasive infections and treatment failures. The data obtained from this study could be applied to devise treatment and management strategies favorable to the hospital or healthcare institution.

John A, Tierce R, Zheng X, Zhou C, Noval MG, O'Keeffe A, Podkowik M, Gonzales S, Inglima K, Desvignes L, Hochman SE, Stapleford KA, Thorpe LE, Pironti A, Shopsin B, Cadwell K, Dittmann M, Torres VJ. SARS-CoV-2 infection predisposes patients to coinfection with Staphylococcus aureus

Severe COVID-19 has been associated with coinfections with bacterial and fungal pathogens. Notably, patients with COVID-19 who develop Staphylococcus aureus bacteremia exhibit higher rates of mortality than those infected with either pathogen alone. To understand this clinical scenario, we collected and examined S. aureus blood and respiratory isolates from a hospital in New York City during the early phase of the pandemic from both SARS-CoV-2+ and SARS-CoV-2- patients. Whole genome sequencing of these S. aureus isolates revealed broad phylogenetic diversity in both patient groups, suggesting that SARS-CoV-2 coinfection was not associated with a particular S. aureus lineage. Phenotypic characterization of the contemporary collection of S. aureus isolates from SARS-CoV-2+ and SARS-CoV-2- patients revealed no notable differences in several virulence traits examined. However, we noted a trend toward overrepresentation of S. aureus bloodstream strains with low cytotoxicity in the SARS-CoV-2+ group. We observed that patients coinfected with SARS-CoV-2 and S. aureus were more likely to die during the acute phase of infection when the coinfecting S. aureus strain exhibited high or low cytotoxicity. To further investigate the relationship between SARS-CoV-2 and S. aureus infections, we developed a murine coinfection model. These studies revealed that infection with SARS-CoV-2 renders mice susceptible to subsequent superinfection with low cytotoxicity S. aureus. Thus, SARS-CoV-2 infection sensitizes the host to coinfections, including S. aureus isolates with low intrinsic virulence.

IMPORTANCE

The COVID-19 pandemic has had an enormous impact on healthcare across the globe. Patients who were severely infected with SARS-CoV-2, the virus causing COVID-19, sometimes became infected with other pathogens, which is termed coinfection. If the coinfecting pathogen is the bacterium Staphylococcus aureus, there is an increased risk of patient death. We collected S. aureus strains that coinfected patients with SARS-CoV-2 to study the disease outcome caused by the interaction of these two important pathogens. We found that both in patients and in mice, coinfection with an S. aureus strain lacking toxicity resulted in more severe disease during the early phase of infection, compared with infection with either pathogen alone. Thus, SARS-CoV-2 infection can directly increase the severity of S. aureus infection.

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.

NanoMGT: Marker gene typing of low complexity mono-species metagenomic samples using noisy long reads

Rapid advancements in sequencing technologies have led to significant progress in microbial genomics, yet challenges persist in accurately identifying microbial strain diversity in metagenomic samples, especially when working with noisy long-read data from platforms like Oxford Nanopore Technologies (ONT). In this article, we introduce NanoMGT, a tool designed to enhance marker gene typing in low-complexity mono-species samples, leveraging the unique properties of long reads. NanoMGT excels in its ability to accurately identify mutations amidst high error rates, ensuring the reliable detection of multiple strain-specific marker genes. Our tool implements a novel scoring system that rewards mutations co-occurring across different reads and penalizes densely grouped, likely erroneous variants, thereby achieving a good balance between sensitivity and precision. A comparative evaluation of NanoMGT, using a simulated multi-strain sample of seven bacterial species, demonstrated superior performance relative to existing tools and the advantages of using a threshold-based filtering approach to calling minority variants in ONT's sequencing data. NanoMGT's potential as a post-binning tool in metagenomic pipelines is particularly notable, enabling researchers to more accurately determine specific alleles and understand strain diversity in microbial communities. Our findings have significant implications for clinical diagnostics, environmental microbiology, and the broader field of genomics. The findings offer a reliable and efficient approach to marker gene typing in complex metagenomic samples.

Complete genome sequence of a Salmonella enterica strain of serovar Agona associated with an outbreak in France

We report here the complete genome of one Salmonella Agona strain isolated in 2017 from a dried milk powder in France.

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.

Genomic study of European Clostridioides difficile ribotype 002/sequence type 8

Clostridioides difficile has significant clinical importance as a leading cause of healthcare-associated infections, with symptoms ranging from mild diarrhoea to severe colitis, and possible life-threatening complications. C. difficile ribotype (RT) 002, mainly associated with MLST sequence type (ST) 8, is one of the most common RTs found in humans. This study aimed at investigating the genetic characteristics of 537 C. difficile genomes of ST8/RT002. To this end, we sequenced 298 C. difficile strains representing a new European genome collection, with strains from Germany, Denmark, France and Portugal. These sequences were analysed against a global dataset consisting of 1,437 ST8 genomes available through Enterobase. Our results showed close genetic relatedness among the studied ST8 genomes, a diverse array of antimicrobial resistance (AMR) genes and the presence of multiple mobile elements. Notably, the pangenome analysis revealed an open genomic structure. ST8 shows relatively low overall variation. Thus, clonal isolates were found across different One Health sectors (humans, animals, environment and food), time periods, and geographical locations, suggesting the lineage's stability and a universal environmental source. Importantly, this stability did not hinder the acquisition of AMR genes, emphasizing the adaptability of this bacterium to different selective pressures. Although only 2.4 % (41/1,735) of the studied genomes originated from non-human sources, such as animals, food, or the environment, we identified 9 cross-sectoral core genome multilocus sequence typing (cgMLST) clusters. Our study highlights the importance of ST8 as a prominent lineage of C. difficile with critical implications in the context of One Health. In addition, these findings strongly support the need for continued surveillance and investigation of non-human samples to gain a more comprehensive understanding of the epidemiology of C. difficile.

A retrospective investigation of the population structure and geospatial distribution of Salmonella Paratyphi A in Kathmandu, Nepal

Salmonella Paratyphi A, one of the major etiologic agents of enteric fever, has increased in prevalence in recent decades in certain endemic regions in comparison to S. Typhi, the most prevalent cause of enteric fever. Despite this increase, data on the prevalence and molecular epidemiology of S. Paratyphi A remain generally scarce. Here, we analysed the whole genome sequences of 216 S. Paratyphi A isolates originating from Kathmandu, Nepal between 2005 and 2014, of which 200 were from patients with acute enteric fever and 16 from the gallbladder of people with suspected chronic carriage. By exploiting the recently developed genotyping framework for S. Paratyphi A (Paratype), we identified several genotypes circulating in Kathmandu. Notably, we observed an unusual clonal expansion of genotype 2.4.3 over a four-year period that spread geographically and systematically replaced other genotypes. This rapid genotype replacement is hypothesised to have been driven by both reduced susceptibility to fluoroquinolones and genetic changes to virulence factors, such as functional and structural genes encoding the type 3 secretion systems. Finally, we show that person-to-person is likely the most common mode of transmission and chronic carriers seem to play a limited role in maintaining disease circulation.

Genomic insights of Salmonella isolated from dry fermented sausage production chains in Spain and France

The presence of Salmonella in dry fermented sausages is source of recalls and outbreaks. The genomic diversity of 173 Salmonella isolates from the dry fermented sausage production chains (pig carcasses, pork, and sausages) from France and Spain were investigated through their core phylogenomic relationships and accessory genome profiles. Ten different serovars and thirteen sequence type profiles were identified. The most frequent serovar from sausages was the monophasic variant of S. Typhimurium (1,4,[5],12:i:-, 72%) while S. Derby was in pig carcasses (51%). Phylogenomic clusters found in S. 1,4,[5],12:i:-, S. Derby, S. Rissen and S. Typhimurium serovars identified closely related isolates, with less than 10 alleles and 20 SNPs of difference, displaying Salmonella persistence along the pork production chain. Most of the S. 1,4,[5],12:i:- contained the Salmonella genomic island-4 (SGI-4), Tn21 and IncFIB plasmid. More than half of S. Derby strains contained the SGI-1 and Tn7. S. 1,4,[5],12:i:- genomes carried the most multidrug resistance genes (91% of the strains), whereas extended-spectrum β-lactamase genes were found in Typhimurium and Derby serovars. Salmonella monitoring and characterization in the pork production chains, specially S. 1,4,[5],12:i:- serovar, is of special importance due to its multidrug resistance capacity and persistence in dry fermented sausages.

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.

Molecular Characterization of Klebsiella pneumoniae Clinical Isolates Through Whole-Genome Sequencing: A Comprehensive Analysis of Serotypes, Sequence Types, and Antimicrobial and Virulence Genes

Introduction Antimicrobial resistance (AMR) has become a menace, spreading among bacterial species globally. AMR is now recognized as a silent pandemic responsible for treatment failures. Therefore, an effective surveillance mechanism is warranted to understand the bacterial species isolated from human clinical specimens. The present study employed next-generation sequencing (NGS) or whole-genome sequencing (WGS) to identify the resistance and virulence genes, sequence type, and serotypes. Methods This study included 18 multidrug-resistant (MDR) Klebsiella pneumoniae (K. pneumoniae) isolates obtained from patients suffering from different infections attending the Prathima Institute of Medical Sciences, Karimnagar, India. All isolates were identified, and antimicrobial susceptibility profiles were determined through conventional microbiological techniques and confirmed by automated systems. All the isolates were investigated using NGS or WGS to identify the genes coding for resistance, such as extended-spectrum beta-lactamases (ESBLs), metallo-beta-lactamases, and virulence genes. Multilocus sequence typing (MLST) was conducted to identify the sequence types, and Kleborate analysis was performed to confirm the species, genes for AMR, and virulence and evaluate the capsular polysaccharide (KL) and cell wall/lipopolysaccharide (O) serotypes carried by the isolates. Results The mean age of the patients was 46.11±20.35 years. Among the patients included, 12 (66.66%) were males and 6 (33.33%) were females. A high percentage (>50%) of hypervirulent K. pneumoniae (hvKp) strains that had genes coding for AMR and plasmids having the potential to carry blaNDM and resistance genes were observed. Among the isolates, 16 (88.88%) revealed the presence of multiple antibiotic-resistant genes with evidence of at least one gene coding for beta-lactamase resistance. There was a high prevalence of blaSHV (17/18; 94.44%) and blaCTX-M-15 (16/18; 88.88%) AMR genes. Other AMR genes identified included blaTEM (83.33%; 15/18) and blaOXA (14/18; 77.77%). Two (11.11%) strains each showed the presence of blaNDM-1 and blaNDM-5 genes. The virulence genes identified included gapA, infB, mdh, pgi, phoE, rpoB, tonB, and ybt. The most frequent K. pneumoniae serotypes found were KL51:O1v2 (3/18, 16.66%), KL17:O1v1 (3/18, 16.66%), and KL64:O2v1 (3/18, 16.66%). KL64 (4/18; 22.22%) was the most common capsular serotype identified among the isolates. The most frequent MLST-based sequence type (ST) identified included ST-147 (5/18, 27.77%), followed by ST-231 (3/18, 16.66%) and ST-101 (2/18, 11.11%). Conclusions The molecular analysis of K. pneumoniae isolates revealed multiple AMR, plasmid, and virulence genes. Additionally, many global STs were noticed by MLST. The results noted a high prevalence of hvKp strains. Molecular characterization of bacterial strains using NGS/WGS is important to understand the epidemiology of bacterial strains and the antibiotic resistance and virulence genes they are potentially carrying. The data obtained from this study may be utilized to devise careful antibiotic-prescribing approaches and improve patient management practices.

Molecular Characterization of Escherichia coli Causing Urinary Tract Infections Through Next-Generation Sequencing: A Comprehensive Analysis of Serotypes, Sequence Types, and Antimicrobial and Virulence Genes

Introduction An enormous increase in antimicrobial resistance (AMR) among bacteria isolated from human clinical specimens contributed to treatment failures. Increased surveillance through next-generation sequencing (NGS) or whole genome sequencing (WGS) could facilitate the study of the epidemiology of drug-resistant bacterial strains, resistance genes, and other virulence determinants they are potentially carrying. Methods This study included 30 Escherichia coli (E. coli) isolates obtained from patients suffering from urinary tract infections (UTIs) attending Prathima Institute of Medical Sciences, Karimnagar, India. All bacterial isolates were identified, and antimicrobial susceptibility patterns were determined through conventional microbiological techniques and confirmed by automated systems. All the isolates were investigated using NGS to identify genes coding for resistance, such as extended-spectrum beta-lactamases (ESBLs), metallo-beta-lactamases, and virulence genes. Multilocus sequence typing (MLST) was used to understand the prevalent strain types, and serotyping was carried out to evaluate the type of O (cell wall antigen) and H (flagellar antigen) serotypes carried by the isolates. Results The conventional antimicrobial susceptibility testing revealed that 15 (50%) isolates were resistant to imipenem (IPM), 10 (33.33%) were resistant to amikacin (AK), 13 (43.33%) were resistant to piperacillin-tazobactam (PTZ), 17 (56.66%) were resistant to cephalosporins, and 14 (46.66%) were resistant to nitrofurantoin (NIT). Among the isolates, 26 (86.66%) had revealed the presence of multiple antibiotic-resistant genes with evidence of at least one gene coding for beta-lactamase resistance. There was a high prevalence of blaCTX-M (19/30, 63.33%) genes, followed by blaTEM and blaOXA-1. The blaNDM-5 gene was found in three isolates (3/30, 10%). The virulence genes identified in the present study were iutA, sat, iss, and papC, among others. The E. coli serotype found predominantly belonged to O25:H4 (5, 16.66%), followed by O102:H6 (4, 13.33%). A total of 16 MLST variants were identified among the examined samples. Of the MLST-based sequence types (STs) identified, ST-131 (7, 23.33%) was the predominant one, followed by ST-167 (3, 10%) and ST-12 (3, 10%). Conclusions The study results demonstrated that the E. coli strains isolated from patients suffering from UTIs potentially carried antimicrobial resistance and virulence genes and belonged to different strain types based on MLST. Careful evaluation of bacterial strains using molecular analyses such as NGS could facilitate an improved understanding of bacterial antibiotic resistance and its virulence potential. This could enable physicians to choose appropriate antimicrobial agents and contribute to better patient management, thereby preventing the emergence and spread of drug-resistant bacteria.

ContScout: sensitive detection and removal of contamination from annotated genomes

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Methods for Genomic Epidemiology of Bacterial Pathogens: Example Salmonella

Genomics has revolutionized how we characterize and monitor infectious diseases for public health. The surveillance and characterization of Salmonella has improved drastically within the past decade. In this chapter, we discuss the prerequisites for good bacterial genomics studies and make note of advantages and disadvantages of this research approach. We discuss methods for outbreak detection and the evolutionary and epidemiological characterization of Salmonella spp. We provide an outline for determining the sequence type and serotype of isolates, building a core genome phylogenetic tree, and detecting antimicrobial resistance genes, virulence factors, and mobile genetic elements. These methods can be used to study other pathogenic bacterial species.

Three novel sequencing types from seventeen Staphylococcus aureus genomes isolated from dairy cows milk in the Free State Province of South Africa

Staphylococcus aureus is one of the major pathogens causing bovine mastitis, which results in huge economic losses in the dairy industry worldwide. Here, we report genome sequences of 17 S. aureus strains, with three novel sequencing types (ST8495, ST8500, and ST8501) isolated from the milk of dairy cows with subclinical mastitis.

Identification of plasmids in avian-associated Escherichia coli using nanopore and illumina sequencing

BACKGROUND

Avian pathogenic Escherichia coli (APEC) are the causative agents of colibacillosis in chickens, a disease which has significant economic impact on the poultry industry. Large plasmids detected in APEC are known to contribute to strain diversity for pathogenicity and antimicrobial resistance, but there could be other plasmids that are missed in standard analysis. In this study, we determined the impact of sequencing and assembly factors for the detection of plasmids in an E. coli whole genome sequencing project.

RESULTS

Hybrid assembly (Illumina and Nanopore) combined with plasmid DNA extractions allowed for detection of the greatest number of plasmids in E. coli, as detected by MOB-suite software. In total, 79 plasmids were identified in 19 E. coli isolates. Hybrid assemblies were robust and consistent in quality regardless of sequencing kit used or if long reads were filtered or not. In contrast, long read only assemblies were more variable and influenced by sequencing and assembly parameters. Plasmid DNA extractions allowed for the detection of physically smaller plasmids, but when averaged over 19 isolates did not significantly change the overall number of plasmids detected.

CONCLUSIONS

Hybrid assembly can be reliably used to detect plasmids in E. coli, especially if researchers are focused on large plasmids containing antimicrobial resistance genes and virulence factors. If the goal is comprehensive detection of all plasmids, particularly if smaller sized vectors are desired for biotechnology applications, the addition of plasmid DNA extractions to hybrid assemblies is prudent. Long read sequencing is sufficient to detect many plasmids in E. coli, however, it is more prone to errors when expanded to analyze a large number of isolates.

Cont-ID: detection of sample cross-contamination in viral metagenomic data

BACKGROUND

High-throughput sequencing (HTS) technologies completed by the bioinformatic analysis of the generated data are becoming an important detection technique for virus diagnostics. They have the potential to replace or complement the current PCR-based methods thanks to their improved inclusivity and analytical sensitivity, as well as their overall good repeatability and reproducibility. Cross-contamination is a well-known phenomenon in molecular diagnostics and corresponds to the exchange of genetic material between samples. Cross-contamination management was a key drawback during the development of PCR-based detection and is now adequately monitored in routine diagnostics. HTS technologies are facing similar difficulties due to their very high analytical sensitivity. As a single viral read could be detected in millions of sequencing reads, it is mandatory to fix a detection threshold that will be informed by estimated cross-contamination. Cross-contamination monitoring should therefore be a priority when detecting viruses by HTS technologies.

RESULTS

We present Cont-ID, a bioinformatic tool designed to check for cross-contamination by analysing the relative abundance of virus sequencing reads identified in sequence metagenomic datasets and their duplication between samples. It can be applied when the samples in a sequencing batch have been processed in parallel in the laboratory and with at least one specific external control called Alien control. Using 273 real datasets, including 68 virus species from different hosts (fruit tree, plant, human) and several library preparation protocols (Ribodepleted total RNA, small RNA and double-stranded RNA), we demonstrated that Cont-ID classifies with high accuracy (91%) viral species detection into (true) infection or (cross) contamination. This classification raises confidence in the detection and facilitates the downstream interpretation and confirmation of the results by prioritising the virus detections that should be confirmed.

CONCLUSIONS

Cross-contamination between samples when detecting viruses using HTS (Illumina technology) can be monitored and highlighted by Cont-ID (provided an alien control is present). Cont-ID is based on a flexible methodology relying on the output of bioinformatics analyses of the sequencing reads and considering the contamination pattern specific to each batch of samples. The Cont-ID method is adaptable so that each laboratory can optimise it before its validation and routine use.

Whole genome sequencing of increased number of azithromycin-resistant Shigella flexneri 1b isolates in Ontario

Azithromycin (AZM) resistance among Shigella is a major public health concern. Here, we investigated the epidemiology of Shigella flexneri serotype 1b recovered during 2016-2018 in Ontario, to describe the prevalence and spread of AZM resistance. We found that 72.3% (47/65) of cases were AZM-resistant (AZMR), of which 95.7% (45/47) were among males (P < 0.001). Whole-genome based phylogenetic analysis showed three major clusters, and 56.9% of isolates grouped within a single closely-related cluster (0-10 ∆SNP). A single AZMR clonal population was persistent over 3 years and involved 67.9% (36/53) of all male cases, and none reported international travel. In 2018, a different AZMR cluster appeared among adult males not reporting travel. A proportion of isolates (10.7%) with reduced susceptibility to ciprofloxacin (CIP) due to S83L mutation in gyrA were AZM susceptible, and 71.4% reported international travel. Resistance to AZM was due to the acquisition of mph gene-bearing incFII plasmids having > 95% nucleotide similarity to pKSR100. Plasmid-borne resistance limiting treatment options to AZM, ceftriaxone (CRO) and CIP was noted in a single isolate. We characterized AZMR isolates circulating locally among males and found that genomic analysis can support targeted prevention and mitigation strategies against antimicrobial-resistance.

Transmission dynamics variability of lineage 2 Mycobacterium tuberculosis strains in Kobe, Japan, determined using population-based whole-genome sequencing analysis

Currently, tuberculosis (TB) in Japan is highly prevalent among elderly patients who were born during a time when TB was highly prevalent. Mycobacterium tuberculosis (Mtb) lineage 2 (L2) is the predominant strain in the country. Moreover, the proportion of foreign-born patients with TB has been increasing. This epidemiological situation in Japan motivated us to explore the heterogeneity in transmission dynamics among the sublineages of Mtb L2 within this aging population. For this purpose, we conducted a population-based whole genome sequencing analysis of 550 Mtb strains in Kobe, Japan, and employed pairwise single nucleotide polymorphism (SNP) distance clustering and terminal branch length (TBL) distribution analysis to assess Mtb transmission. The genomic clustering rate with a threshold of ≤5 SNPs was significantly lower in elderly patients aged 70 years or higher than in non-elderly patients. The elderly patient group showed significantly longer TBL than the non-elderly group. These results supported the notion that reactivation of distant infection is a major driving force for the high incidence of TB in elderly individuals. The age group distribution and frequency of lineages/sublineages were found to significantly differ between foreign-born and Japan-born patients. The increased proportion of foreign-born patients might have resulted in more strain diversity in Japan. The L2.2.A sublineage demonstrated a significant association with elderly patients and exhibited lower transmission rates, which indicate to be prone to reactivate from long-term latency. In contrast, L2.2.Modern, showed a strong association with younger and foreign-born patients. This sublineage showed a high genomic cluster rate, suggesting its high transmissibility. The other three major sublineages, namely L2.2.AA2, L2.2.AA3.1, and L2.2.AA3.2, exhibited a consistent increase in cluster rates across varying SNP thresholds, indicating their relatively recent emergence as endemic sublineages in Japan. In conclusion, this study highlights distinct differences in the transmission dynamics of L2 sublineages within an aging society.

Harmonization of supervised machine learning practices for efficient source attribution of Listeria monocytogenes based on genomic data

BACKGROUND

Genomic data-based machine learning tools are promising for real-time surveillance activities performing source attribution of foodborne bacteria such as Listeria monocytogenes. Given the heterogeneity of machine learning practices, our aim was to identify those influencing the source prediction performance of the usual holdout method combined with the repeated k-fold cross-validation method.

METHODS

A large collection of 1 100 L. monocytogenes genomes with known sources was built according to several genomic metrics to ensure authenticity and completeness of genomic profiles. Based on these genomic profiles (i.e. 7-locus alleles, core alleles, accessory genes, core SNPs and pan kmers), we developed a versatile workflow assessing prediction performance of different combinations of training dataset splitting (i.e. 50, 60, 70, 80 and 90%), data preprocessing (i.e. with or without near-zero variance removal), and learning models (i.e. BLR, ERT, RF, SGB, SVM and XGB). The performance metrics included accuracy, Cohen's kappa, F1-score, area under the curves from receiver operating characteristic curve, precision recall curve or precision recall gain curve, and execution time.

RESULTS

The testing average accuracies from accessory genes and pan kmers were significantly higher than accuracies from core alleles or SNPs. While the accuracies from 70 and 80% of training dataset splitting were not significantly different, those from 80% were significantly higher than the other tested proportions. The near-zero variance removal did not allow to produce results for 7-locus alleles, did not impact significantly the accuracy for core alleles, accessory genes and pan kmers, and decreased significantly accuracy for core SNPs. The SVM and XGB models did not present significant differences in accuracy between each other and reached significantly higher accuracies than BLR, SGB, ERT and RF, in this order of magnitude. However, the SVM model required more computing power than the XGB model, especially for high amount of descriptors such like core SNPs and pan kmers.

CONCLUSIONS

In addition to recommendations about machine learning practices for L. monocytogenes source attribution based on genomic data, the present study also provides a freely available workflow to solve other balanced or unbalanced multiclass phenotypes from binary and categorical genomic profiles of other microorganisms without source code modifications.

Antimicrobial overproduction sustains intestinal inflammation by inhibiting Enterococcus colonization

Loss of antimicrobial proteins such as REG3 family members compromises the integrity of the intestinal barrier. Here, we demonstrate that overproduction of REG3 proteins can also be detrimental by reducing a protective species in the microbiota. Patients with inflammatory bowel disease (IBD) experiencing flares displayed heightened levels of secreted REG3 proteins that mediated depletion of Enterococcus faecium (Efm) from the gut microbiota. Efm inoculation of mice ameliorated intestinal inflammation through activation of the innate immune receptor NOD2, which was associated with the bacterial DL-endopeptidase SagA that generates NOD2-stimulating muropeptides. NOD2 activation in myeloid cells induced interleukin-1β (IL-1β) secretion to increase the proportion of IL-22-producing CD4+ T helper cells and innate lymphoid cells that promote tissue repair. Finally, Efm was unable to protect mice carrying a NOD2 gene variant commonly found in IBD patients. Our findings demonstrate that inflammation self-perpetuates by causing aberrant antimicrobial activity that disrupts symbiotic relationships with gut microbes.

Whole-Genome Sequencing-Based Screening of MRSA in Patients and Healthcare Workers in Public Hospitals in Benin

Methicillin-resistant Staphylococcus aureus (MRSA) constitutes a serious public health concern, with a considerable impact on patients' health, and substantial healthcare costs. In this study, patients and healthcare workers (HCWs) from six public hospitals in Benin were screened for MRSA. Strains were identified as MRSA using conventional microbiological methods in Benin, and confirmed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry in Belgium. Whole-genome sequencing (WGS) was used on the confirmed MRSA isolates, to characterize their genomic content and study their relatedness. Amongst the 305 isolates (304 wound swabs and 61 nasal swabs) that were collected from patients and HCWs, we detected 32 and 15 cases of MRSA, respectively. From this collection, 27 high-quality WGS datasets were obtained, which carried numerous genes and mutations associated with antimicrobial resistance. The mecA gene was detected in all the sequenced isolates. These isolates were assigned to five sequence types (STs), with ST8 (55.56%, n = 15/27), ST152 (18.52%, n = 5/27), and ST121 (18.52%, n = 5/27) being the most common. These 27 isolates carried multiple virulence genes, including the genes encoding the Panton-Valentine leukocidin toxin (48.15%, n = 13/27), and the tst gene (29.63%, n = 8/27), associated with toxic shock syndrome. This study highlights the need to implement a multimodal strategy for reducing the risk of the cross-transmission of MRSA in hospitals.

The Good, The Bad, and The Misleading: How to Improve the Quality of 'Genome Announcements'?

When comparing the requirements of diverse journals to publish microbial 'Genome Reports,' we noticed that some mostly focus on benchmarking universal single-copy orthologs scores as a quality measure, while the exclusion of possible contaminating sequences from genomic resources and the possible misidentification of the target microbes receive less attention. To deal with these quality issues, we suggest that DNA barcodes that are widely accepted for the identification of the target microbe species should be extracted from newly reported genome resources and included in phylogenetic analyses to confirm the identity of the sequenced microorganisms before Genome Reports are published. This approach, applied, for example, by the journal IMA Fungus, largely prevents the misidentification of the microbes that are targeted for whole-genome sequencing (WGS). In addition, contig similarity values, including GC content, remapping coverage of WGS reads, and BLASTN searches against the National Center for Biotechnology Information nucleotide database, would also reveal contamination issues. The values of these two recommendations to improve the publication criteria for microbial Genome Reports in diverse journals are demonstrated here through analyses of a draft genome published in Molecular Plant-Microbe Interactions and then retracted due to contaminations. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Role of Mobile Genetic Elements in Virulence Factor Carriage from Symptomatic and Asymptomatic Cases of Escherichia coli Bacteriuria

Uropathogenic Escherichia coli (UPEC) is extremely diverse genotypically and phenotypically. Individual strains can variably carry diverse virulence factors, making it challenging to define a molecular signature for this pathotype. For many bacterial pathogens, mobile genetic elements (MGEs) constitute a major mechanism of virulence factor acquisition. For urinary E. coli, the total distribution of MGEs and their role in the acquisition of virulence factors is not well defined, including in the context of symptomatic infection versus asymptomatic bacteriuria (ASB). In this work, we characterized 151 isolates of E. coli, derived from patients with either urinary tract infection (UTI) or ASB. For both sets of E. coli, we catalogued the presence of plasmids, prophage, and transposons. We analyzed MGE sequences for the presence of virulence factors and antimicrobial resistance genes. These MGEs were associated with only ~4% of total virulence associated genes, while plasmids contributed to ~15% of antimicrobial resistance genes under consideration. Our analyses suggests that, across strains of E. coli, MGEs are not a prominent driver of urinary tract pathogenesis and symptomatic infection. IMPORTANCE Escherichia coli is the most common etiological agent of urinary tract infections (UTIs), with UTI-associated strains designated "uropathogenic" E. coli or UPEC. Across urinary strains of E. coli, the global landscape of MGEs and its relationship to virulence factor carriage and clinical symptomatology require greater clarity. Here, we demonstrate that many of the putative virulence factors of UPEC are not associated with acquisition due to MGEs. The current work enhances our understanding of the strain-to-strain variability and pathogenic potential of urine-associated E. coli and points toward more subtle genomic differences distinguishing ASB from UTI isolates.

Enteric fever cluster identification in South Africa using genomic surveillance of Salmonella enterica serovar Typhi

The National Institute for Communicable Diseases in South Africa participates in national laboratory-based surveillance for human isolates of Salmonella species. Laboratory analysis includes whole-genome sequencing (WGS) of isolates. We report on WGS-based surveillance of Salmonella enterica serovar Typhi (Salmonella Typhi) in South Africa from 2020 through 2021. We describe how WGS analysis identified clusters of enteric fever in the Western Cape Province of South Africa and describe the epidemiological investigations associated with these clusters. A total of 206 Salmonella Typhi isolates were received for analysis. Genomic DNA was isolated from bacteria and WGS was performed using Illumina NextSeq technology. WGS data were investigated using multiple bioinformatics tools, including those available at the Centre for Genomic Epidemiology, EnteroBase and Pathogenwatch. Core-genome multilocus sequence typing was used to investigate the phylogeny of isolates and identify clusters. Three major clusters of enteric fever were identified in the Western Cape Province; cluster one (n=11 isolates), cluster two (n=13 isolates), and cluster three (n=14 isolates). To date, no likely source has been identified for any of the clusters. All isolates associated with the clusters, showed the same genotype (4.3.1.1.EA1) and resistome (antimicrobial resistance genes: bla _TEM-1B, catA1, sul1, sul2, dfrA7). The implementation of genomic surveillance of Salmonella Typhi in South Africa has enabled rapid detection of clusters indicative of possible outbreaks. Cluster identification allows for targeted epidemiological investigations and a timely, coordinated public health response.

A panel of diverse Klebsiella pneumoniae clinical isolates for research and development

Klebsiella pneumoniae are a leading cause of healthcare-associated infections worldwide. In particular, strains expressing extended-spectrum β-lactamases (ESBLs) and carbapenemases pose serious treatment challenges, leading the World Health Organization (WHO) to designate ESBL and carbapenem-resistant Enterobacteriaceae as 'critical' threats to human health. Research efforts to combat these pathogens can be supported by accessibility to diverse and clinically relevant isolates for testing novel therapeutics. Here, we describe a panel of 100 diverse K. pneumoniae isolates that are publicly available to assist the research community in this endeavour. Whole-genome sequencing (WGS) was performed on 3878 K. pneumoniae clinical isolates housed at the Multidrug-Resistant Organism Repository and Surveillance Network. The isolates were cultured from 63 facilities in 19 countries between 2001 and 2020. Core-genome multilocus sequence typing and high-resolution single-nucleotide polymorphism-based phylogenetic analyses captured the genetic diversity of the collection and were used to select the final panel of 100 isolates. In addition to known multidrug-resistant (MDR) pandemic lineages, the final panel includes hypervirulent lineages and isolates with specific and diverse resistance genes and virulence biomarkers. A broad range of antibiotic susceptibilities, ranging from pan-sensitive to extensively drug-resistant isolates, are described. The panel collection, and all associated metadata and genome sequences, are available at no additional cost and will be an important resource for the research community and for the design and development of novel antimicrobial agents and diagnostics against this important pathogen.

Capsular Polysaccharide Is Essential for the Virulence of the Antimicrobial-Resistant Pathogen Enterobacter hormaechei

Nosocomial infections caused by multidrug-resistant (MDR) Enterobacter cloacae complex (ECC) pathogens are on the rise. However, the virulence strategies employed by these pathogens remain elusive. Here, we study the interaction of ECC clinical isolates with human serum to define how this pathogen evades the antimicrobial action of complement, one of the first lines of host-mediated immune defense. We identified a small number of serum-sensitive strains, including Enterobacter hormaechei strain NR3055, which we exploited for the in vitro selection of serum-resistant clones. Comparative genomics between the serum-sensitive NR3055 strain and the isolated serum-resistant clones revealed a premature stop codon in the wzy gene of the capsular polysaccharide biosynthesis locus of NR3055. The complementation of wzy conferred serum resistance to NR3055, prevented the deposition of complement proteins on the bacterial surface, inhibited phagocytosis by human neutrophils, and rendered the bacteria virulent in a mouse model of peritonitis. Mice exposed to a nonlethal dose of encapsulated NR3055 were protected from subsequent lethal infections by encapsulated NR3055, whereas mice that were previously exposed to unencapsulated NR3055 succumbed to infection. Thus, capsule is a key immune evasion determinant for E. hormaechei, and it is a potential target for prophylactics and therapeutics to combat these increasingly MDR human pathogens. IMPORTANCE Infections caused by antimicrobial resistant bacteria are of increasing concern, especially those due to carbapenem-resistant Enterobacteriaceae pathogens. Included in this group are species of the Enterobacter cloacae complex, regarding which there is a paucity of knowledge on the infection biology of the pathogens, despite their clinical relevance. In this study, we combine techniques in comparative genomics, bacterial genetics, and diverse models of infection to establish capsule as an important mechanism of Enterobacter pathogens to resist the antibacterial activity of serum, a first line of host defense against bacterial infections. We also show that immune memory targeting the Enterobacter capsule protects against lethal infection. The further characterization of Enterobacter infection biology and the immune response to infection are needed for the development of therapies and preventative interventions targeting these highly antibiotic resistant pathogens.

Genomic Characterization of Salmonella enterica serovar Weltevreden Associated with Human Diarrhea

Salmonella Weltevreden is an emerging pathogen associated with human diarrhea, and knowledge of the genomics and epidemiology of this serovar is still limited. In this study, we performed whole-genome sequencing of 96 S. Weltevreden isolates recovered from diarrheal patients and 62 isolates from food animals in China between 2006 and 2017. Together, with an additional 199 genome sequences of S. Weltevreden published in NCBI, we performed an analysis on all 357 S. Weltevreden genome sequences. Our results demonstrated that the majority of S. Weltevreden from diarrheal patients from China (97.92%, 94/96) and the other regions in the world (94.97%, 189/199) identified in this study were sequence type (ST) 365. The remaining types were ST3771 (n = 3), ST22 (n = 1), ST155 (n = 1), and ST684 (n = 1). In addition, ST365 was also widely recovered from animals, food, and environmental samples in different regions of the world. Phylogenetic analysis and pulsed-field gel electrophoresis (PFGE) revealed that S. Weltevreden from diarrheal patients was closely related to those recovered from food and environmental specimens. We also showed that S. Weltevreden did not exhibit severe antimicrobial resistance profiles, suggesting administering antibiotics is still effective for controlling the agent. Interestingly, we found that S. Weltevreden strains carried a number of virulence factor genes, and a 100.03-kb IncFII(S) type plasmid was widely distributed in S. Weltevreden strains. Elimination of this plasmid decreased the bacterial capacity to infect both Caco-2 cells and C57BL/6 mice, suggesting the importance of this plasmid for bacterial virulence. Our results contribute to the understanding of the epidemiology and virulence of S. Weltevreden. IMPORTANCE Salmonella Weltevreden is a pathogen associated with human diarrheal diseases found across the globe. However, knowledge of the genomics and epidemiology of this pathogen is still limited. In this study, we found S. Weltevreden sequence type (ST) 365 is commonly recovered from diarrheal patients in China and many other regions of the world, and there is no major difference between the Chinese isolates and the global isolates at the phylogenetic level. We also demonstrated that ST365 was widely recovered from animal, food, and environmental samples collected in different, global regions. Importantly, we discovered an IncFII(S) type plasmid commonly carried by S. Weltevreden strains of human, animal, and food origins, and this plasmid is likely to contribute to the bacterial pathogenesis. These findings enhance our understanding of the emergence of S. Weltevreden involved in diarrheal outbreaks and the global spread of S. Weltevreden strains.

Tell me if you prefer bovine or poultry sectors and I'll tell you who you are: Characterization of Salmonella enterica subsp. enterica serovar Mbandaka in France

Introduction

In north-western France, Salmonella enterica susp. enterica serovar Mbandaka (S. Mbandaka) is most frequently isolated from bovine and dairy samples. While this serovar most often results in asymptomatic carriage, for a number of years it has caused episodes of abortions, which have serious economic consequences for the sector. Interestingly, this serovar is also isolated from Gallus gallus in the same geographic zone. Despite its prevalence in bovines in north-western France, S. Mbandaka has not been broadly studied at the genomic level, and its prevalence and host adaptation are still not fully understood.

Methods

In this study, we analyzed the genomic diversity of 304 strains of S. Mbandaka isolated from the bovine and poultry sectors in this area over a period of 5 years. A phylogenetic analysis was carried out and two approaches were followed to identify conserved genes and mutations related to host associations. The first approach targeted the genes compiled in the MEGARESv2, Resfinder, VFDB and SPI databases. Plasmid and phage contents were also investigated. The second approach refers to an in-house algorithm developed for this study that computes sensitivity, specificity, and accuracy of accessory genes and core variants according to predefined genomes groups.

Results and discussion

All the analyzed strains belong to the multi-locus sequence type profile ST413, and the phylogenomic analysis revealed main clustering by host (bovine and poultry), emphasizing the circulation of 12 different major clones, of which seven circulate in poultry and five in the bovine sector in France and a likely food production chain adaptation of these clones. All strains present resistance determinants including heavy metals and biocides that could explain the ability of this serovar to survive and persist in the environment, within herds, and in food processing plants. To explore the wild animal contribution to the spread of this serovar in north-western France, we retrieved S. Mbandaka genomes isolated from wild birds from EnteroBase and included them in the phylogenomic analysis together with our collection. Lastly, screening of accessory genes and major variants allowed us to identify conserved specific mutations characteristic of each major cluster. These mutations could be used to design useful probes for food safety surveillance.

Non-random genetic alterations in the cyanobacterium Nostoc sp. exposed to space conditions

Understanding the impact of long-term exposure of microorganisms to space is critical in understanding how these exposures impact the evolution and adaptation of microbial life under space conditions. In this work we subjected Nostoc sp. CCCryo 231-06, a cyanobacterium capable of living under many different ecological conditions, and also surviving in extreme ones, to a 23-month stay at the International Space Station (the Biology and Mars Experiment, BIOMEX, on the EXPOSE-R2 platform) and returned it to Earth for single-cell genome analysis. We used microfluidic technology and single cell sequencing to identify the changes that occurred in the whole genome of single Nostoc cells. The variant profile showed that biofilm and photosystem associated loci were the most altered, with an increased variant rate of synonymous base pair substitutions. The cause(s) of these non-random alterations and their implications to the evolutionary potential of single bacterial cells under long-term cosmic exposure warrants further investigation.

Background Filtering of Clinical Metagenomic Sequencing with a Library Concentration-Normalized Model

Metagenomic next-generation sequencing (mNGS) can accurately detect pathogens in clinical samples. However, wet-lab contamination constrains mNGS analysis and may result in erroneous interpretation of results. Many existing methods rely on large-scale observational microbiome studies and may not be applicable to clinical mNGS tests. By generation of a pretrained profile of common laboratory contaminants, we developed an mNGS noise-filtering model based on the inverse linear relationship between microbial sequencing reads and sample library concentration, named the background elimination and correction by library concentration-normalized (BECLEAN) model. Its efficacy was evaluated with bacteria- and yeast-spiked samples and 28 cerebrospinal fluid (CSF) specimens. The diagnostic accuracy, precision, sensitivity, and specificity of BECLEAN with reference to conventional methods and diagnosis were 92.9%, 86.7%, 100%, and 86.7%, respectively. BECLEAN led to a dramatic reduction of background noise without affecting the true-positive rate and thus can provide a time-saving and convenient tool in various clinical settings. IMPORTANCE Most of the existing methods to remove wet-lab contamination rely on large-scale observational microbiome studies and may not be applicable to clinical mNGS testing in individual cases. In clinical settings, only a handful of samples might be sequenced in a run. The lab-specific microbiome can complicate existing statistical approaches for removing contamination from small-scale clinical metagenomic sequencing data sets; thus, use of a preliminary lab-specific training set is necessary. Our study provides a rapid and accurate background-filtering tool for clinical metagenomic sequencing by generation of a pretrained profile of common laboratory contaminants. Notably, our work demonstrates that the inverse linear relationship between microbial sequencing reads and library concentration can serve to identify true contaminants and evaluate the relative abundance of a taxon in samples by comparing the observed microbial reads to the model-predicted value. Our findings extend the previously published research and demonstrate confirmatory results in clinical settings.

Brucella ceti and Brucella pinnipedialis genome characterization unveils genetic features that highlight their zoonotic potential

The Gram-negative bacteria Brucella ceti and Brucella pinnipedialis circulate in marine environments primarily infecting marine mammals, where they cause an often-fatal disease named brucellosis. The increase of brucellosis among several species of cetaceans and pinnipeds, together with the report of sporadic human infections, raises concerns about the zoonotic potential of these pathogens on a large scale and may pose a threat to coastal communities worldwide. Therefore, the characterization of the B. ceti and B. pinnipedialis genetic features is a priority to better understand the pathological factors that may impact global health. Moreover, an in-depth functional analysis of the B. ceti and B. pinnipedialis genome in the context of virulence and pathogenesis was not undertaken so far. Within this picture, here we present the comparative whole-genome characterization of all B. ceti and B. pinnipedialis genomes available in public resources, uncovering a collection of genetic tools possessed by these aquatic bacterial species compared to their zoonotic terrestrial relatives. We show that B. ceti and B. pinnipedialis genomes display a wide host-range infection capability and a polyphyletic phylogeny within the genus, showing a genomic structure that fits the canonical definition of closeness. Functional genome annotation led to identifying genes related to several pathways involved in mechanisms of infection, others conferring pan-susceptibility to antimicrobials and a set of virulence genes that highlight the similarity of B. ceti and B. pinnipedialis genotypes to those of Brucella spp. displaying human-infecting phenotypes.

Genomic elements located in the accessory repertoire drive the adaptation to biocides in Listeria monocytogenes strains from different ecological niches

In response to the massive use of biocides for controlling Listeria monocytogenes (hereafter Lm) contaminations along the food chain, strains showing biocide tolerance emerged. Here, accessory genomic elements were associated with biocide tolerance through pangenome-wide associations performed on 197 Lm strains from different lineages, ecological, geographical and temporal origins. Mobile elements, including prophage-related loci, the Tn6188_qacH transposon and pLMST6_emrC plasmid, were widespread across lineage I and II food strains and associated with tolerance to benzalkonium-chloride (BC), a quaternary ammonium compound (QAC) widely used in food processing. The pLMST6_emrC was also associated with tolerance to another QAC, the didecyldimethylammonium-chloride, displaying a pleiotropic effect. While no associations were detected for chemically reactive biocides (alcohols and chlorines), genes encoding for cell-surface proteins were associated with BC or polymeric biguanide tolerance. The latter was restricted to lineage I strains from animal and the environment. In conclusion, different genetic markers, with polygenic nature or not, appear to have driven the Lm adaptation to biocide, especially in food strains but also from animal and the environment. These markers could aid to monitor and predict the spread of biocide tolerant Lm genotypes across different ecological niches, finally reducing the risk of such strains in food industrial settings.

Core Genome Multilocus Sequence Typing Scheme for Improved Characterization and Epidemiological Surveillance of Pathogenic Brucella

Brucellosis poses a significant burden to human and animal health worldwide. Robust and harmonized molecular epidemiological approaches and population studies that include routine disease screening are needed to efficiently track the origin and spread of Brucella strains. Core genome multilocus sequence typing (cgMLST) is a powerful genotyping system commonly used to delineate pathogen transmission routes for disease surveillance and control. Except for Brucella melitensis, cgMLST schemes for Brucella species are currently not established. Here, we describe a novel cgMLST scheme that covers multiple Brucella species. We first determined the phylogenetic breadth of the genus using 612 Brucella genomes. We selected 1,764 genes that were particularly well conserved and typeable in at least 98% of these genomes. We tested the new scheme on 600 genomes and found high agreement with the whole-genome-based single nucleotide polymorphism (SNP) analysis. Next, we applied the scheme to reanalyze the genome of Brucella strains from epidemiologically linked outbreaks. We demonstrated the applicability of the new scheme for high-resolution typing required in outbreak investigations as previously reported with whole-genome SNP methods. We also used the novel scheme to define the global population structure of the genus using 1,322 Brucella genomes. Finally, we demonstrated the possibility of tracing distribution of Brucella strains by performing cluster analysis of cgMLST profiles and found nearly identical cgMLST profiles in different countries. Our results show that sequencing depth of more than 40-fold is optimal for allele calling with this scheme. In summary, this study describes a novel Brucella-wide cgMLST scheme that is applicable in Brucella molecular epidemiology and helps in accurately tracking and thus controlling the sources of infection. The scheme is publicly accessible and should represent a valuable resource for laboratories with limited computational resources and bioinformatics expertise.

A retrospective and regional approach assessing the genomic diversity of Salmonella Dublin

From a historically rare serotype, Salmonella enterica subsp. enterica Dublin slowly became one of the most prevalent Salmonella in cattle and raw milk cheese in some regions of France. We present a retrospective genomic analysis of 480 S. Dublin isolates to address the context, evolutionary dynamics, local diversity and the genesis processes of regional S. Dublin outbreaks events between 2015 and 2017. Samples were clustered and assessed for correlation against metadata including isolation date, isolation matrices, geographical origin and epidemiological hypotheses. Significant findings can be drawn from this work. We found that the geographical distance was a major factor explaining genetic groups in the early stages of the cheese production processes (animals, farms) while down-the-line transformation steps were more likely to host genomic diversity. This supports the hypothesis of a generalised local persistence of strains from animal to finished products, with occasional migration. We also observed that the bacterial surveillance is representative of diversity, while targeted investigations without genomics evidence often included unrelated isolates. Combining both approaches in phylogeography methods allows a better representation of the dynamics, of outbreaks.