Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores

Genomic insights into the plasmidome of non-tuberculous mycobacteria

BACKGROUND

Non-tuberculous mycobacteria (NTM) are a diverse group of environmental bacteria that are increasingly associated with human infections and difficult to treat. Plasmids, which might carry resistance and virulence factors, remain largely unexplored in NTM.

METHODS

We used publicly available complete genome sequence data of 328 NTM isolates belonging to 125 species to study gene content, genomic diversity, and clusters of 196 annotated NTM plasmids. Furthermore, we analyzed 3755 draft genome assemblies from over 200 NTM species and 5415 short-read sequence datasets from six clinically relevant NTM species or complexes including M. abscessus, M. avium complex, M. ulcerans complex and M. kansasii complex, for the presence of these plasmids.

RESULTS

Between one and five plasmids were present in approximately one-third of the complete NTM genomes. The annotated plasmids varied widely in length (most between 10 and 400 kbp) and gene content, with many genes having an unknown function. Predicted gene functions primarily involved plasmid replication, segregation, maintenance, and mobility. Only a few plasmids contained predicted genes that are known to confer resistance to antibiotics commonly used to treat NTM infections. Out of 196 annotated plasmid sequences, 116 could be grouped into 31 clusters of closely related sequences, and about one-third were found across multiple NTM species. Among clinically relevant NTM, the presence of NTM plasmids showed significant variation between species, within (sub)species, and even among strains within (sub)lineages, such as dominant circulating clones of Mycobacterium abscessus.

CONCLUSIONS

Our analysis demonstrates that plasmids are a diverse and heterogeneously distributed feature in NTM bacteria. The frequent occurrence of closely related putative plasmid sequences across different NTM species suggests they may play a significant role in NTM evolution through horizontal gene transfer at least in some groups of NTM. However, further in vitro investigations and access to more complete genomes are necessary to validate our findings, elucidate gene functions, identify novel plasmids, and comprehensively assess the role of plasmids in NTM.

Genomic and phenotypic characterization of Chryseomicrobium imtechense from canine pyometra

Chryseomicrobium imtechense is a bacterium primarily isolated in environment samples. Here we demonstrated the virulent and antimicrobial resistant profile of the C. imtechense LBV029/19 strain, which was isolated in a pyometra infection in canines, being related as causal agent of that disease. The bacterium was recovered from purulent content of a uterus with pyometra and the identity was confirmed by both, biochemical tests followed MALDI-TOF MS and 16S-rDNA gene sequencing. Antimicrobial susceptibility test, biofilm formation assessment, genomic sequencing with Illumina MiSeq platform followed by the genome characterization and phylogenetic relationships were conducted with C. imtechense LBV029/19 strain. The phylogenetic analysis revealed a close evolutionary relationship between the C. imtechense LBV029/19 and the others publicly available 16S-rDNA genes of C. imtechense. The in vitro assays shown that C. imtechense LBV029/19 has a strong ability for biofilm formation. Additionally, the antimicrobial susceptibility tests revealed a smaller inhibition zone for novobiocin, cefuroxime, and erythromycin. To the best of our knowledge the C. imtechense LBV029/19 genome is the first genome sequenced and available of C. imtechense. Regarding the genotypic profile, resistance genes to tetracyclines, fluoroquinolones, beta-lactams, and multidrug resistance were present in this genome. Additionally, genes encoded to virulence factors, including biofilm formation, motility, and sporulation were identified. The phenotypic profile associated with the genotypic description indicates the potential of C. imtechense LBV029/19 as a primary agent of animal infectious diseases. The genomic description will contribute to the understanding of the biology of the bacterium and the pathogenesis process involved in the disease.

A novel virulent core genome multilocus sequence type CT 11424 of Listeria monocytogenes isolate causing stillbirth in Bangladesh

BACKGROUND

Listeria monocytogenes is a foodborne pathogen that can lead to severe pregnancy outcomes. This study reports the clinical and genomic characteristics of a Listeria-mediated stillbirth identified in January 2022 through the Child Health and Mortality Prevention Surveillance (CHAMPS) project in Bangladesh. The Lm-BD-CHAMPS-01 isolate was recovered from the blood and cerebrospinal fluid (CSF) of a male stillborn. Maternal history, clinical, and demographic data were collected by the CHAMPS surveillance platform. An expert panel evaluated all reports to determine the role of L. monocytogenes infection in the causal chain of stillbirth. Genomic characterization included multilocus sequence typing (MLST), core genome MLST (cgMLST), serotyping, and the presence or absence of virulence genes. Genetic divergence and phylogenetic analyses were conducted to determine the relationship with other reported isolates globally.

RESULTS

The isolate Lm-BD-CHAMPS-01 was identified as a novel cgMLST CT11424. It belonged to ST 308, Serotype 4b, Clonal Complex 1, and Phylogenetic Lineage 1. Key L. monocytogenes virulence genes facilitating the crossing of the placental barrier, including full-length inlA, LIPI-1, and LIPI-3, were detected. The isolate was closely related to clinical L. monocytogenes isolates, as determined by GrapeTree based on cgMLST. SNP-based phylogenetic analysis found Lm-BD-CHAMPS-01 to be the most distant from other CC1 isolates in the database. Possible sources of infection included the consumption of contaminated raw vegetables or exposure to pigeons.

CONCLUSIONS

This is the first genome sequence of clinical L. monocytogenes from Bangladesh, which also caused stillbirth. Rural healthcare professionals should be aware of L. monocytogenes infection risks during pregnancy. Pregnant women should be counseled on the dangers of exposure to animals or birds and consumption of potentially contaminated raw food to prevent adverse pregnancy outcomes due to L. monocytogenes infection.

Genomic analysis of antimicrobial resistant Escherichia coli isolated from manure and manured agricultural grasslands

Antimicrobial resistance (AMR) is a multifactorial issue involving an intertwining relationship between animals, humans and the environment. The environment can harbour Escherichia coli that are pathogenic or commensal. Escherichia coli is used as an indicator of environmental faecal contamination. Through culture dependent approaches this study identified 46 E. coli isolates in porcine and bovine manure, non-manured and manured soil, and manured grass. The grass isolation highlights grass as an environmental reservoir for E. coli. We also identified a diverse plasmidome with 23 different plasmid replicon types. The E. coli isolates were phenotypically antimicrobial resistant, predominantly multidrug resistant. Whole genome sequencing identified 31 antimicrobial resistance genes, and mutations in the gyrA, parC, and parE genes, conferring fluoroquinolone resistance. This study demonstrates grass as an understudied environmental niche of AMR E. coli, which directly links the environment to the grass grazing animal and vice-versa via the circular economy of manure application.

The influence of depth on the global deep-sea plasmidome

Plasmids play a crucial role in facilitating genetic exchange and enhancing the adaptability of microbial communities. Despite their importance, environmental plasmids remain understudied, particularly those in fragile and underexplored ecosystems such as the deep-sea. In this paper we implemented a bioinformatics pipeline to study the composition, diversity, and functional attributes of plasmid communities (plasmidome) in 81 deep-sea metagenomes from the Tara and Malaspina expeditions, sampled from the Pacific, Atlantic, and Indian Oceans at depths ranging from 270 to 4005 m. We observed an association between depth and plasmid traits, with the 270-1000 m range (mesopelagic samples) exhibiting the highest number of plasmids and the largest plasmid sizes. Plasmids of Alphaproteobacteria and Gammaproteobacteria were predominant across the oceans, particularly in this depth range, which also showed the highest species diversity and abundance of metabolic pathways, including aromatic compound degradation. Surprisingly, relatively few antibiotic resistance genes were found in the deep-sea ecosystem, with most being found in the mesopelagic layer. These included classes such as beta-lactamase, biocide resistance, and aminoglycosides. Our study also identified the MOBP and MOBQ relaxase families as prevalent across various taxonomic classes. This research underscores the importance of studying the plasmidome independently from the chromosomal context. Our limited understanding of the deep-sea's microbial ecology, especially its plasmidome, necessitates caution in human activities like mining. Such activities could have unforeseen impacts on this largely unexplored ecosystem.

Antimicrobial resistant Enterobacterales of clinical importance in mute swans

Urban water environments, including canals, harbours and estuaries are susceptible to contamination with antimicrobials and drug-resistant bacteria through domestic and industrial wastewater discharges and storm water overflows. There is potential for wildlife using these waters to acquire and transmit drug-resistant bacteria and antimicrobial resistance genes (ARGs) of clinical importance. This study aimed to assess clinically important drug-resistant bacteria in urban waterfowl, particularly mute swans. Faecal samples were collected from 17 mute swans in the Greater Dublin Area, Ireland during July, August, October, and November of 2022. Samples were swabbed directly onto agars to select for carbapenem resistant, Extended-spectrum Beta-lactamase (ESBL)-producing, ciprofloxacin resistant and colistin resistant bacteria. Isolates identified by MALDI-TOF as Enterobacterales were tested for susceptibility to a panel of 16 antimicrobials and real-time PCR was employed to detect cefotaximase and carbapenemase genes (CRGs). Drug-resistant isolates were characterised by Whole Genome Sequencing (WGS), including long read sequencing for carbapenemase and mobile colistin resistance (mcr) gene-producing Enterobacterales isolates. Eleven of seventeen (65 %) swan samples were positive for the resistant organism(s) (n = 35), comprising Escherichia coli (n = 32; 82 %) and other Enterobacterales (n = 3). Twenty E. coli (63 %) produced ESBL, with 16/20 (80 %) identified as positive for blaCTX-M-group 1 enzymes, comprising CTX-M-15 (n = 13), CTX-M-55 (n = 2) and CTX-M-1 (n = 1) and 4/20 (20 %) positive for blaCTX-M-group 9 enzymes CTX-M-27 (n = 2) and CTX-M-9 (n = 2). Three E. coli isolates were phenotypically ertapenem resistant, one of which was an ST4450 isolate which carried plasmid encoded blaOXA-181 and blaCMY-141 with blaCTX-M-15 identified chromosomally. One colistin resistant E. coli bore the mcr-1 gene chromosomally. Bioinformatic analysis revealed high-risk pathogenic ESBL E. coli clones including ST38 (n = 3), ST69 (n = 3), and ST131 (n = 2). The study indicates mute swans are a reservoir for drug-resistant Enterobacterales and ARGs of clinical importance and may be a useful sentinel species for antimicrobial resistance (AMR) surveillance in wildlife.

Using genomics to explore the epidemiology of vancomycin resistance in a sewage system

VanHAX-mediated glycopeptide resistance has been consistently high in one of the three main sewer systems in Copenhagen, Lynetten, for +20 years. To explore this for other glycopeptide resistance genes, and whether the colonization has resulted in establishment of multiple bacterial taxa, we mapped 505 shotgun metagenomic data sets from the inlet of three sewage treatment plants to 831 different glycopeptide resistance genes. Only vanHAX and vanHBX genes were differentially abundant in Lynetten. Analyses of eight contigs suggested limited variations in the flanking regions. Proximity ligation metagenomic analysis of 12 samples from Lynetten identified 441 and 5 paired reads mapping to vanHAX and vanHBX, respectively. The other end of these reads was mapped to generated metagenomic-assembled genomes and NCBI using BLAST. vanHBX could only be linked to the phylum level (Bacillota). Plasmid analysis of vanHBX Hi-C contigs showed that these were mainly located on plasmids reported found in enterococci species. Most vanHAX-linked reads could only be linked to phylum and class level, but some reads were assigned to Enterococcus faecium (7 reads), Enterococcus faecalis (4 reads), Paenibacillus apiarius (2 reads), and Paenibacillus thiaminolyticus (27 reads). Ten of the 20 Hi-C contigs-containing vanHAX were annotated as plasmid, all reported found in Enterococcus species. This study shows that while Hi-C technology is valuable for linking antimicrobial resistance genes to bacterial taxa, it suffers from challenges in reliably mapping the linked read to a genomic region with sufficient taxonomic information. Our results also suggest that over the +20 years of colonizing a sewer system, vanHAX has not become widespread across multiple taxa, remaining primarily in E. faecalis and E. faecium, with the exception of Paenibacillus.IMPORTANCELong-term colonization of microbial communities with antimicrobial-resistant bacteria is expected to result in sharing of the resistance genes between several different bacterial taxa of the communities. We investigated microbiomes from a sewer, which have been colonized with glycopeptide-resistant bacteria harboring the mobile vanHAX gene cluster for a minimum of 20 years, using metagenomics sequencing and Hi-C. We found that despite the long-term presence in the sewer, the vanHAX genes have seemingly not disseminated widely.

Description of genome sequences of arthropod-associated spirochetes of the genus Entomospira

Spirochetal bacteria isolated from arthropods of the genera Culex and Aedes are termed BR149, BR151 (Entomospira culicis), BR193 (Entomospira entomophila), and BR208 (Entomospira nematocerorum). Genome sizes assembled from Illumina MiSeq and Oxford Nanopore reads varied between 1.67 and 1.78 Mb containing three to six plasmids. GC content ranged from 38.5% to 45.76%.

The mobilome landscape of biocide-resistance in Brazilian ESKAPE isolates

The increasing frequency of antibiotic-resistant bacteria is a constant threat to global human health. Therefore, the pathogens of the ESKAPE group (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, and Enterobacter spp.) are among the most relevant causes of hospital infections responsible for millions of deaths every year. However, little has been explored about the danger of microorganisms resistant to biocides such as antiseptics and disinfectants. Widely used in domestic, industrial, and hospital environments, these substances reach the environment and can cause selective pressure for resistance genes and induce cross-resistance to antibiotics, further aggravating the problem. Therefore, it is necessary to use innovative and efficient strategies to monitor the spread of genes related to resistance to biocides. Whole genome sequencing and bioinformatics analysis aiming to search for sequences encoding resistance mechanisms are essential to help monitor and combat these pathogens. Thus, this work describes the construction of a bioinformatics tool that integrates different databases to identify gene sequences that may confer some resistance advantage about biocides. Furthermore, the tool analyzed all the genomes of Brazilian ESKAPE isolates deposited at NCBI and found a series of different genes related to resistance to benzalkonium chloride, chlorhexidine, and triclosan, which were the focus of this work. As a result, the presence of resistance genes was identified in different types of biological samples, environments, and hosts. Regarding mobile genetic elements (MGEs), around 52% of isolates containing genes related to resistance to these compounds had their genes identified in plasmids, and 48.7% in prophages. These data show that resistance to biocides can be a silent, underestimated danger spreading across different environments and, therefore, requires greater attention.

AmpC β-lactamases detected in Southeast Asian Escherichia coli and Klebsiella pneumoniae

Objectives

AmpC β-lactamases are neglected compared with ESBL as a cause of third-generation cephalosporin (3GC) resistance in Enterobacterales in low- and middle-income countries and the burden is unknown. The aim of this study was to investigate the presence of AmpC β-lactamase-producing Escherichia coli and Klebsiella pneumoniae in clinical specimens from three clinical research laboratories in Southeast Asia.

Methods

Stored clinical isolates of E. coli and K. pneumoniae resistant to ceftriaxone or ceftazidime or cefpodoxime and ESBL confirmation test negative were screened using MASTDISCS AmpC, ESBL and Carbapenemase Detection Set-D72C. Short-read WGS was performed to identify ampC genes.

Results

Of 126 isolates collected between 2010 and 2020, 31 (24.6%) and 16 (12.7%) were phenotypically AmpC and inducible AmpC positive by MASTDISCS testing, respectively. All inducible AmpC isolates were ceftriaxone susceptible and 97.7% of AmpC/inducible AmpC isolates tested against cefoxitin were resistant. Through WGS, 17 and eight different STs were detected for the AmpC/inducible AmpC E. coli and K. pneumoniae isolates, respectively. Twelve different β-lactamase resistance genes were detected, with bla CMY-2 most commonly in AmpC-positive isolates (20/31; 64.5%; 15 chromosomal, five plasmid). All inducible AmpC-positive isolates had the bla DHA-1 gene (seven chromosomal, nine plasmid).

Conclusions

Though uncommon, AmpC and inducible AmpC β-lactamases in E. coli and K. pneumoniae are an important cause of infection in Southeast Asia. With current testing methods, these infections may be going undetected, resulting in patients receiving suboptimal treatment.

Characterization and Implications of the Regionally Prevalent ST8333 Strains of Salmonella enterica Serotype 4,[5],12:i: - China, 2017-2023

What is already known about this topic?

Salmonella enterica is a significant foodborne pathogen responsible for millions of infections and thousands of deaths worldwide each year. The increasing prevalence of the ST34 Salmonella enterica serotype 4,5,12:i:-(S. 4,[5],12:i:-) poses a substantial public health threat in China.

What is added by this report?

At the end of 2023, the TraNet early-warning module, utilizing whole-genome sequencing data, identified a newly regionally prevalent sequence type of S. 4,[5],12:i:-, designated ST8333, primarily found in Sichuan Province. Analysis revealed high rates of antimicrobial resistance genes, particularly to aminoglycosides and β-lactams. Plasmid comparisons suggest that ST8333 isolates acquired multidrug-resistance plasmids of swine origin during their distribution in Sichuan Province.

What are the implications for public health practice?

The emergence of ST8333, likely associated with pork consumption, underscores the importance of food safety measures. The high prevalence of multidrug resistance poses challenges for treatment and highlights the need for vigilant surveillance and control measures. Strengthening efforts in regions like Sichuan, identified as a reservoir for multidrug-resistant ST8333 strains, is crucial to prevent further transmission and potential nationwide dissemination. Surveillance and control measures are imperative to curb the spread of multidrug-resistant ST8333 isolates, particularly in high-prevalence regions.

First Report of Carbapenem-Resistant Klebsiella michiganensis Co-Harboring bla KPC-2 and TmexCD2-ToprJ2 Isolated from Wastewater at a Tertiary Hospital in Beijing

Background

Klebsiella michiganensis is an emerging human pathogen that causes nosocomial infections. Its prevalence and spread in the environment should not be ignored. This study identified and characterized Klebsiella michiganensis co-harboring bla KPC-2 and TmexCD2-ToprJ2 in hospital wastewater samples.

Methods

Twelve K. michiganensis strains were isolated from wastewater samples collected at a tertiary hospital in Beijing, China. The genomic characteristics of K. michiganensis strains were analyzed using whole-genome sequences, providing information on the comparison between the genome of K. michiganensis strains and the reference genome, antibiotic resistance genes (ARGs), virulence genes, secretion systems, and mobile genetic elements (plasmids, insertion sequences [ISs], and prophages).

Results

Genome analysis showed that the twelve multi-drug resistant (MDR) strains carried a variety of ARGs and virulence genes, as well as four macromolecular secretion systems (T1SS, T2SS, T5aSS, T5bSS, and T4aP). The genetic environments of both the TmexCD2-ToprJ2 gene cluster and bla KPC-2 gene contained ISs. The plasmids carrying TmexCD2-ToprJ2 gene cluster of nine strains in clade 1 and two strains in clade 2 were annotated as IncR plasmid and rep_cluster_1254 type, respectively. The plasmids carrying bla KPC-2 in 10 strains in clade 1 were identified as IncU, and the plasmids carrying bla KPC-2 in the k11 and k12 strains in clade 2 were IncU and IncX6. The phylogenetic tree and heatmap revealed that the secretion system of type VI (T6SSi) existed in 10 strains in clade 1, and Type IV (T4SS) only existed in the k11 strain in clade 2. In addition, K. michiganensis strains carried 13 plasmids, 14 ISs, and 138 prophages.

Conclusion

In this study, the whole genome sequencing demonstrated the diversity of K. michiganensis genome despite 12 K. michiganensis strains from a hospital wastewater, which lays the foundation for further genetic research and drug resistance gene transmission.

Characterization of antimicrobial resistant Enterobacterales isolated from spinach and soil following zinc amendment

Antimicrobial resistant bacteria can occur in the primary food production environment. The emergence and dissemination of antimicrobial resistance (AMR) in the environment can be influenced by several factors, including the presence of heavy metals. The aim of this study was to examine the presence and characteristics of antimicrobial resistant Enterobacterales in soils and spinach grown in soils with and without zinc amendment. A total of 160 samples (92 soil and 68 spinach) were collected from two locations, in which some plots had been amended with zinc. Samples were cultured on selective agars for detection of extended-spectrum beta-lactamase-producing Enterobacterales (ESBL), carbapenem-resistant Enterobacterales and ciprofloxacin-resistant Enterobacterales. Samples were also cultured for enumeration of total Enterobacterales. Isolates were identified by MALDI-TOF. Antimicrobial susceptibility testing was carried out in accordance with EUCAST and CLSI criteria. The whole genome sequence (WGS) of selected isolates was determined. Inductively coupled plasma atomic emission spectrometry was also performed on soil samples in order to measure the concentration of zinc. In total 20 antimicrobial resistant Enterobacterales were isolated from the soil (n = 8) and spinach samples (n = 12). In both sample types, Serratia fonticola (n = 16) was the dominant species, followed by Escherichia coli (n = 1), Citrobacter freundii (n = 1) and Morganella morganii (n = 1) detected in spinach samples, and Enterobacter cloacae (n = 1) detected in a soil sample. The WGS identified genes conferring resistance to different antimicrobials in agreement with the phenotypic results; 14 S. fonticola isolates were confirmed as ESBL producers and harboured the blaFONA gene. Genes that encoded for zinc resistance and multidrug efflux pumps, transporters that can target both antimicrobials and heavy metals, were also identified. Overall, the findings of this study suggest the presence of zinc did not influence the AMR Enterobacterales in soil or spinach samples.

KleTy: integrated typing scheme for core genome and plasmids reveals repeated emergence of multi-drug resistant epidemic lineages in Klebsiella worldwide

BACKGROUND

Clinically important lineages in Klebsiella, especially those expressing multi-drug resistance (MDR), pose severe threats to public health worldwide. They arose from the co-evolution of the vertically inherited core genome and horizontal gene transfers by plasmids, which has not been systematically explored.

METHODS

We designed KleTy, which consists of dedicated typing schemes for both the core genome and plasmids in Klebsiella. We compared the performance of KleTy with many state-of-the-art pipelines using both simulated and real data.

RESULTS

Employing KleTy, we genotyped 33,272 Klebsiella genomes, categorizing them into 1773 distinct populations and predicting the presence of 87,410 plasmids from 837 clusters (PCs). Notably, Klebsiella is the center of the plasmid-exchange network within Enterobacteriaceae. Our results associated the international emergence of prevalent Klebsiella populations with only four carbapenem-resistance (CR) PCs, two hypervirulent PCs, and two hvCR-PCs encoding both carbapenemase and hypervirulence. Furthermore, we observed the ongoing international emergence of blaNDM, accompanied by the replacement of the previously dominant population, blaKPC-encoding HC1360_8 (CC258), during 2003-2018, with the emerging blaNDM-encoding HC1360_3 (CC147) thereafter. Additionally, expansions of hypervirulent carbapenem-resistant Klebsiella pneumoniae (hvCRKP) were evidenced in both populations, driven by plasmids of MDR-hypervirulence convergences.

CONCLUSIONS

The study illuminates how the global genetic landscape of Klebsiella has been shaped by the co-evolution of both the core genome and the plasmids, underscoring the importance of surveillance and control of the dissemination of plasmids for curtailing the emergence of hvCRKPs.

Genomic dynamics of high-risk carbapenem-resistant klebsiella pneumoniae clones carrying hypervirulence determinants in Egyptian clinical settings

BACKGROUND

Ongoing studies have revealed the global prevalence of severe infections caused by the hypervirulent strains of Klebsiella pneumoniae (K. pneumoniae). Meanwhile, the World Health Organization and the Centers for Disease Control declared carbapenem-resistant K. pneumoniae as an urgent public health threat, requiring swift and effective action to mitigate its spread. Low- and middle-income countries are severely impacted by such devastating infectious diseases owing to the ill implementation of antimicrobial practices and infection control policies. Having both hypervirulence and carbapenemase gene determinants, the emergence of convergent hypervirulent carbapenem-resistant K. pneumoniae is now being reported worldwide.

METHODS

In this study, we sequenced 19 carbapenemase-producing K. pneumoniae strains recovered from various clinical specimens. Additionally, we evaluated the phenotypic antimicrobial susceptibility to multiple antimicrobial classes using the VITEK2 automated system. Utilizing the sequencing data, we characterized the sequence types, serotypes, pangenome, resistance profiles, virulence profiles, and mobile genetic elements of the examined isolates. We highlighted the emergence of high-risk clones carrying hypervirulence genetic determinants among the screened isolates.

RESULTS

Our findings revealed that all carbapenem-resistant isolates exhibited either extensive- or pan-drug resistance and harbored multiple variants of resistance genes spanning nearly all the antimicrobial classes. The most prevalent carbapenemase genes detected within the isolates were blaNDM-5 and blaOXA-48. We identified high-risk clones, such as ST383-K30, ST147-K64, ST11-K15, and ST14-K2, which may have evolved into putative convergent strains by acquiring the full set of hypervirulence-associated genetic determinants (iucABCD, rmpA and/ or rmpA2, putative transporter peg-344). Additionally, this study identified ST709-K9 as a high-risk clone for the first time and uncovered that capsule types K15 and K9 carried hypervirulence genetic determinants. The most frequent Inc types found in these isolates were Col440I, IncHI1B, and Inc FII(K).

CONCLUSION

This study highlights the emergence of high-risk, extensively carbapenem-resistant K. pneumoniae strains co-carrying hypervirulence determinants in Egyptian clinical settings. This poses an imminent threat not only to Egypt but also to the global community, underscoring the urgent need for enhanced surveillance and control strategies to combat this pathogen.

Draft genome sequence of non-pathogenic Escherichia coli 5C LMG S-33222, isolated from healthy donor feces

We present the genome sequence of Escherichia coli 5C LMG S-33222, a non-pathogenic microorganism with potential probiotic features. The strain was isolated in 2021 from a fecal sample of a healthy Italian infant. The total genome size is 4,712,575 bp with a G + C content of 51%.

First report of carbapenems encoding multidrug-resistant gram-negative bacteria from a pediatric hospital in Gaza Strip, Palestine

BACKGROUND

The worldwide prevalence of multi-drug resistance (MDR) in Gram-negative bacteria (GNB), particularly related to extended-spectrum beta-lactamases (ESBLs) and carbapenemases, poses significant global public health and clinical challenges.

OBJECTIVES

To characterize ESBL-producing Gram-negative bacilli, within a pediatric hospital in Gaza using whole genome sequencing (WGS).

METHODS

A total of 158 clinical isolates of Gram-negative bacilli were collected from Al-Nasser Pediatric Hospital. These isolates were tested for ESBL production using the double disk synergy test. The antibiotic susceptibility profile was determined using the Kirby Bauer method following the Clinical and Laboratory Standard Institute guidelines. Selected 15 phenotypically MDR isolates were whole-genome sequenced and characterized for their genome-based species identity and antibiotic resistance gene profile.

RESULTS

Of the 158 isolates, 93 (58.9%) were positive for ESBL production. The frequency of Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa, Acinetobacter baumannii, Proteus mirabilis, and Serratia marcescens was 50%, 22.7%, 22.7%, 1.8%, 1.2%, and 1.2% respectively. The prevalence of ESBL among urine, pus, blood, and sputum was 64%, 44%, 23%, and 63.6%, respectively. Chloramphenicol, Imipenem, and Meropenem were the most effective antibiotics against ESBL producers. In sequenced isolates,  an average of six anti-microbial resistance (AMR) genes were noted per isolate, where one of them carried up to 13 antibiotic resistance genes. Carbapenem resistance genes such as blaKPC-2(6.6%), blaPDC-36/12 (6.6%), and blaPOM-1 (6.6%) were detected. All the sequenced E. coli isolates (n = 8) showed multiple resistance genes, mainly against β-lactamase (25.0%), aminoglycosides (37.5%), sulfonamides (37.5%), and genes conferring resistance to tetracyclines (25.0).

CONCLUSION

Our results showed a high prevalence of ESBL-producing GNB isolated from a pediatric hospital in the Gaza Strip. Various antibiotic resistance genes were identified, including those encoding ESBL and carbapenems. The results highlight the significant challenge posed by MDR in GNB and emphasize the need for effective antibiotic strategies. Given the high endemicity observed in various studies from Palestine, it is important to conduct clinical and molecular epidemiology research to identify risk factors, transmission patterns, and clinical outcomes associated with GNB strains that carry ESBL and carbapenem resistance genes.

Campylobacter coli of porcine origin exhibits an open pan-genome within a single clonal complex: insights from comparative genomic analysis

Introduction

Although Campylobacter spp., including Campylobacter coli, have emerged as important zoonotic foodborne pathogens globally, the understanding of the genomic epidemiology of C. coli of porcine origin is limited.

Methods

As pigs are an important reservoir of C. coli, we analyzed C. coli genomes that were isolated (n = 3) from pigs and sequenced (this study) them along with all other C. coli genomes for which pig intestines, pig feces, and pigs were mentioned as sources in the NCBI database up to January 6, 2023. In this paper, we report the pan-genomic features, the multi-locus sequence types, the resistome, virulome, and mobilome, and the phylogenomic analysis of these organisms that were obtained from pigs.

Results and discussion

Our analysis revealed that, in addition to having an open pan-genome, majority (63%) of the typeable isolates of C. coli of pig origin belonged to a single clonal complex, ST-828. The resistome of these C. coli isolates was predominated by the genes tetO (53%), blaOXA-193 (49%), and APH (3')-IIIa (21%); however, the virulome analysis revealed a core set of 37 virulence genes. Analysis of the mobile genetic elements in the genomes revealed wide diversity of the plasmids and bacteriophages, while 30 transposons were common to all genomes of C. coli of porcine origin. Phylogenomic analysis showed two discernible clusters comprising isolates originating from Japan and another set of isolates comprising mostly copies of a type strain stored in three different culture collections.

A single amplified genome catalog reveals the dynamics of mobilome and resistome in the human microbiome

BACKGROUND

The increase in metagenome-assembled genomes (MAGs) has advanced our understanding of the functional characterization and taxonomic assignment within the human microbiome. However, MAGs, as population consensus genomes, often aggregate heterogeneity among species and strains, thereby obfuscating the precise relationships between microbial hosts and mobile genetic elements (MGEs). In contrast, single amplified genomes (SAGs) derived via single-cell genome sequencing can capture individual genomic content, including MGEs.

RESULTS

We introduce the first substantial SAG dataset (bbsag20) from the human oral and gut microbiome, comprising 17,202 SAGs above medium-quality without co-assembly. This collection unveils a diversity of bacterial lineages across 312 oral and 647 gut species, demonstrating different taxonomic compositions from MAGs. Moreover, the SAGs showed cellular-level evidence of the translocation of oral bacteria to the gut. We also identified broad-host-range MGEs harboring antibiotic resistance genes (ARGs), which were not detected in the MAGs.

CONCLUSIONS

The difference in taxonomic composition between SAGs and MAGs indicates that combining both methods would be effective in expanding the genome catalog. By connecting mobilomes and resistomes in individual samples, SAGs could meticulously chart a dynamic network of ARGs on MGEs, pinpointing potential ARG reservoirs and their spreading patterns in the microbial community. Video Abstract.

Bifidobacterium fermentum sp. nov. and Bifidobacterium aquikefiricola sp. nov., isolated from water kefir

Four strains, representing two novel Bifidobacterium species, were isolated from water kefir, a fermented beverage. 16S rRNA gene analysis suggested that the novel species share high identities (98.82-98.89%) with Bifidobacterium aquikefiri LMG 28769T. Complete genomes were assembled with a short- and long-read hybrid sequencing approach. In agreement with the 16S rRNA gene analysis, phylogenetics with 117 marker genes places the novel species closest to B. aquikefiri LMG 28769T as well. The isolates have average nucleotide identity (ANI) scores ranging from 81.46 to 84.84% and digital DNA-DNA hybridization (dDDH) scores from 23.9 to 38.5% with the closest related species, as well as ANI scores between the proposed new species of 80.50%, indicating that the isolates represent two novel species. Matrix-assisted laser desorption/ionization-time of flight chemotaxonomic analysis supported the gene-based taxonomic placement. We propose the names Bifidobacterium fermentum sp. nov. and Bifidobacterium aquikefiricola sp. nov. for these novel species within the Bifidobacterium genus. The proposed type strain B. fermentum WK012_4_13T (= LMG 33104T = DSM 116073T; GenBank accession number GCF_041080835.1) has a genome size of 2.43 Mbp, with a G+C content of 56.00 mol%. The proposed type strain for B. aquikefiricola WK041_4_12T (= LMG 33105T = DSM 116074T; GenBank accession number GCF_041080795.1) has a genome size of 2.36 Mbp and a G+C content of 53.94 mol%. B. fermentum cells are Gram-positive staining, non-motile, non-spore-forming, fructose-6-phosphate phosphoketolase (F6PPK)-positive, catalase- and oxidase-negative and bacillary club shaped. B. aquikefiricola cells are Gram-positive staining, non-motile, non-spore-forming, F6PPK-positive, catalase- and oxidase-negative and square rod shaped.

Transmission of dominant strains of Campylobacter jejuni and Campylobacter coli between farms and retail stores in Ecuador: Genetic diversity and antimicrobial resistance

Thermotolerant Campylobacter is an important zoonotic pathogen known for causing gastroenteritis in humans, with poultry as its primary reservoir. A total of 468 samples were collected, of which 335 were chicken carcass samples (representing the food component), and 133 were chicken caeca samples (representing the animal component). These samples underwent culture, with colonies examined under a microscope. Species identification was achieved through multiplex PCR. Additionally, antimicrobial susceptibility profiles were determined using the Kirby-Bauer method, testing for sensitivity to gentamicin, ciprofloxacin, tetracycline, and erythromycin. Additionally, 55 C. jejuni (62.5%) and 33 C. coli (37.5%) isolates were selected for whole genome sequencing (WGS). A High prevalence of Campylobacter was observed, with rates of 95.5% (n = 127, CI95%: 92.5% - 98.5%) in the animal component and 72.5% (n = 243, CI95%: 69.9% - 75.1%) in the food component. Specifically, C. jejuni was detected in 33.1% (n = 42) of poultry farms and 38.3% (n = 93) of chicken carcasses, while C. coli was found in 64.6% (n = 82) of poultry farms and 60.5% (n = 147) of chicken carcasses. Antimicrobials with the highest rates of resistance (67%-100%) were ciprofloxacin and tetracycline, in both animal and food component isolates. Erythromycin resistance was notable, ranging from 22% to 33%, with only two C. jejuni isolates from retail were resistant to gentamicin. Furthermore, multidrug resistance was identified in 23% (20 isolates) of the Campylobacter isolates. Genetic analysis revealed the presence of fourteen resistance genes in both C. jejuni and C. coli isolates, including tet(O), blaOXA-460, blaOXA-184, blaOXA-489, blaOXA-193, blaOXA-784, blaOXA-603, aph(3')-IIIa, aad9, aph(2'')-If, aadE-Cc, sat4, and ant(6)-Ia. Additionally, twenty-five plasmids were detected in the 88 Campylobacter isolates examined. Interestingly, most isolates also harbored genes encoding putative virulence factors associated with pathogenicity, invasion, adherence, and production of cytolethal distending toxin (cdt): cheV, cheA, cheW, cheY, flaA, flgR, flaC, flaD, flgB, flgC, ciaB, ciaC. The WGS analysis showed the presence of several cgSTs in both animal and food components, with nine of them widely disseminated between components. Moreover, C. coli and C. jejuni isolates from different sources presented less than 11 single nucleotide polymorphisms (SNPs), suggesting clonality (16 isolates). Further analysis using SNP tree demonstrated widespread distribution of certain C. jejuni and C. coli clones across multiple farms and retail stores. This study presents, for the first-time, insights into the clonality, plasmid diversity, virulence, and antimicrobial resistance (AMR) of thermotolerant Campylobacter strains originating from the Ecuadorian poultry industry. The identification of AMR genes associated with the main antibiotics used in the treatment of campylobacteriosis in humans, highlights the importance of the prudent use of antimicrobials in the poultry industry. Additionally, this research remarks the need for regional studies to understand the epidemiology of this pathogen.

Insights into the ecology of the infant gut plasmidome

Plasmids are small DNA molecules that enable bacteria to share beneficial traits, influencing microbial communities. However, their role within the human gut microbiome remains largely unknown. In this study, we investigate the gut microbiomes of 34 mother-child cohorts, employing a plasmid analysis workflow to understand the impact of plasmids on the gut microbiome. We create a plasmid phylogenetic tree, devise a method for assigning plasmid hosts, and examine potential plasmid transfer networks. Our research discovers a wide variety of previously unidentified plasmid sequences, indicating that current databases do not fully represent the gut plasmidome. Interestingly, infants display greater plasmid diversity compared to mothers and other healthy adults. We find that Bacteroidota, a major bacterial phylum, serves as the primary host for gut plasmids and plays a dominant role in gut plasmid transfer events. Additionally, plasmids broaden the genetic capabilities of bacteria, with their influence on bacterial function becoming more apparent as children's gut microbiomes develop. This study sheds light on the role of plasmids in the infant gut microbiome, making a significant contribution to our understanding of plasmid biology.

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Diversity and phage sensitivity to phages of porcine enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli (ETEC) is a diverse and poorly characterized E. coli pathotype that causes diarrhea in humans and animals. Phages have been proposed for the veterinary biocontrol of ETEC, but effective solutions require understanding of porcine ETEC diversity that affects phage infection. Here, we sequenced and analyzed the genomes of the PHAGEBio ETEC collection, gathering 79 diverse ETEC strains isolated from European pigs with post-weaning diarrhea (PWD). We identified the virulence factors characterizing the pathotype and several antibiotic resistance genes on plasmids, while phage resistance genes and other virulence factors were mostly chromosome encoded. We experienced that ETEC strains were highly resistant to Enterobacteriaceae phage infection. It was only by enrichment of numerous diverse samples with different media and conditions, using the 41 ETEC strains of our collection as hosts, that we could isolate two lytic phages that could infect a large part of our diverse ETEC collection: vB_EcoP_ETEP21B and vB_EcoS_ETEP102. Based on genome and host range analyses, we discussed the infection strategies of the two phages and identified components of lipopolysaccharides ( LPS) as receptors for the two phages. Our detailed computational structural analysis highlights several loops and pockets in the tail fibers that may allow recognition and binding of ETEC strains, also in the presence of O-antigens. Despite the importance of receptor recognition, the diversity of the ETEC strains remains a significant challenge for isolating ETEC phages and developing sustainable phage-based products to address ETEC-induced PWD.IMPORTANCEEnterotoxigenic Escherichia coli (ETEC)-induced post-weaning diarrhea is a severe disease in piglets that leads to weight loss and potentially death, with high economic and animal welfare costs worldwide. Phage-based approaches have been proposed, but available data are insufficient to ensure efficacy. Genome analysis of an extensive collection of ETEC strains revealed that phage defense mechanisms were mostly chromosome encoded, suggesting a lower chance of spread and selection by phage exposure. The difficulty in isolating lytic phages and the molecular and structural analyses of two ETEC phages point toward a multifactorial resistance of ETEC to phage infection and the importance of extensive phage screenings specifically against clinically relevant strains. The PHAGEBio ETEC collection and these two phages are valuable tools for the scientific community to expand our knowledge on the most studied, but still enigmatic, bacterial species-E. coli.

Complete genome sequence of Listeria seeligeri strain 43-1 isolated from a Danish forest freshwater swamp

Here, we report the genome sequence of Listeria seeligeri 43-1 isolated from a Danish freshwater swamp using Oxford Nanopore sequencing. The isolate shared a high genomic similarity to two other L. seeligeri isolates from soil and water.

Benchmarking short and long read polishing tools for nanopore assemblies: achieving near-perfect genomes for outbreak isolates

BACKGROUND

Oxford Nanopore provides high throughput sequencing platforms able to reconstruct complete bacterial genomes with 99.95% accuracy. However, even small levels of error can obscure the phylogenetic relationships between closely related isolates. Polishing tools have been developed to correct these errors, but it is uncertain if they obtain the accuracy needed for the high-resolution source tracking of foodborne illness outbreaks.

RESULTS

We tested 132 combinations of assembly and short- and long-read polishing tools to assess their accuracy for reconstructing the genome sequences of 15 highly similar Salmonella enterica serovar Newport isolates from a 2020 onion outbreak. While long-read polishing alone improved accuracy, near perfect accuracy (99.9999% accuracy or ~ 5 nucleotide errors across the 4.8 Mbp genome, excluding low confidence regions) was only obtained by pipelines that combined both long- and short-read polishing tools. Notably, medaka was a more accurate and efficient long-read polisher than Racon. Among short-read polishers, NextPolish showed the highest accuracy, but Pilon, Polypolish, and POLCA performed similarly. Among the 5 best performing pipelines, polishing with medaka followed by NextPolish was the most common combination. Importantly, the order of polishing tools mattered i.e., using less accurate tools after more accurate ones introduced errors. Indels in homopolymers and repetitive regions, where the short reads could not be uniquely mapped, remained the most challenging errors to correct.

CONCLUSIONS

Short reads are still needed to correct errors in nanopore sequenced assemblies to obtain the accuracy required for source tracking investigations. Our granular assessment of the performance of the polishing pipelines allowed us to suggest best practices for tool users and areas for improvement for tool developers.

Phylogenomic and genomic analysis reveals unique and shared genetic signatures of Mycobacterium kansasii complex species

Species belonging to the Mycobacterium kansasii complex (MKC) are frequently isolated from humans and the environment and can cause serious diseases. The most common MKC infections are caused by the species M. kansasii (sensu stricto), leading to tuberculosis-like disease. However, a broad spectrum of virulence, antimicrobial resistance and pathogenicity of these non-tuberculous mycobacteria (NTM) are observed across the MKC. Many genomic aspects of the MKC that relate to these broad phenotypes are not well elucidated. Here, we performed genomic analyses from a collection of 665 MKC strains, isolated from environmental, animal and human sources. We inferred the MKC pangenome, mobilome, resistome, virulome and defence systems and show that the MKC species harbours unique and shared genomic signatures. High frequency of presence of prophages and different types of defence systems were observed. We found that the M. kansasii species splits into four lineages, of which three are lowly represented and mainly in Brazil, while one lineage is dominant and globally spread. Moreover, we show that four sub-lineages of this most distributed M. kansasii lineage emerged during the twentieth century. Further analysis of the M. kansasii genomes revealed almost 300 regions of difference contributing to genomic diversity, as well as fixed mutations that may explain the M. kansasii's increased virulence and drug resistance.

The detailed analysis of the microbiome and resistome of artisanal blue-veined cheeses provides evidence on sources and patterns of succession linked with quality and safety traits

BACKGROUND

Artisanal cheeses usually contain a highly diverse microbial community which can significantly impact their quality and safety. Here, we describe a detailed longitudinal study assessing the impact of ripening in three natural caves on the microbiome and resistome succession across three different producers of Cabrales blue-veined cheese.

RESULTS

Both the producer and cave in which cheeses were ripened significantly influenced the cheese microbiome. Lactococcus and the former Lactobacillus genus, among other taxa, showed high abundance in cheeses at initial stages of ripening, either coming from the raw material, starter culture used, and/or the environment of processing plants. Along cheese ripening in caves, these taxa were displaced by other bacteria, such as Tetragenococcus, Corynebacterium, Brevibacterium, Yaniella, and Staphylococcus, predominantly originating from cave environments (mainly food contact surfaces), as demonstrated by source-tracking analysis, strain analysis at read level, and the characterization of 613 metagenome-assembled genomes. The high abundance of Tetragenococcus koreensis and Tetragenococcus halophilus detected in cheese has not been found previously in cheese metagenomes. Furthermore, Tetragenococcus showed a high level of horizontal gene transfer with other members of the cheese microbiome, mainly with Lactococcus and Staphylococcus, involving genes related to carbohydrate metabolism functions. The resistome analysis revealed that raw milk and the associated processing environments are a rich reservoir of antimicrobial resistance determinants, mainly associated with resistance to aminoglycosides, tetracyclines, and β-lactam antibiotics and harbored by aerobic gram-negative bacteria of high relevance from a safety point of view, such as Escherichia coli, Salmonella enterica, Acinetobacter, and Klebsiella pneumoniae, and that the displacement of most raw milk-associated taxa by cave-associated taxa during ripening gave rise to a significant decrease in the load of ARGs and, therefore, to a safer end product.

CONCLUSION

Overall, the cave environments represented an important source of non-starter microorganisms which may play a relevant role in the quality and safety of the end products. Among them, we have identified novel taxa and taxa not previously regarded as being dominant components of the cheese microbiome (Tetragenococcus spp.), providing very valuable information for the authentication of this protected designation of origin artisanal cheese. Video Abstract.

Tracing clinically-relevant antimicrobial resistances in Acinetobacter baumannii-calcoaceticus complex across diverse environments: A study spanning clinical, livestock, and wastewater treatment settings

Acinetobacter baumannii has become a prominent nosocomial pathogen, primarily owing to its remarkable ability to rapidly acquire resistance to a wide range of antimicrobial agents and its ability to persist in diverse environments. However, there is a lack of data on the molecular epidemiology and its potential implications for public health of A. baumannii strains exhibiting clinically significant resistances that originate from non-clinical environments. Therefore, the genetic characteristics and resistance mechanisms of 80 A. baumannii-calcoaceticus (ABC) complex isolates, sourced from environments associated with poultry and pig production, municipal wastewater treatment plants (WWTPs), and clinical settings, were investigated. In total, our study classified 54 isolates into 29 previously described sequence types (STs), while 26 isolates exhibited as-yet-unassigned STs. We identified a broad range of A. baumannii STs originating from poultry and pig production environments (e.g., ST10, ST238, ST240, ST267, ST345, ST370, ST372, ST1112 according to Pasteur scheme). These STs have also been documented in clinical settings worldwide, highlighting their clinical significance. These findings also raise concerns about the potential zoonotic transmission of certain STs associated with livestock environments. Furthermore, we observed that clinical isolates exhibited the highest diversity of antimicrobial resistance genes (ARGs). In contrast to non-clinical isolates, clinical isolates typically carried a significantly higher number of ARGs, ranging from 10 to 15. They were also the exclusive carriers of biocide resistance genes and acquired carbapenemases (blaOXA-23, blaOXA-58, blaOXA-72, blaGIM-1, blaNDM-1). Additionally, we observed that clinical strains displayed an increased capacity for carrying plasmids and undergoing genetic transformation. This heightened capability could be linked to the intense selective pressures commonly found within clinical settings. Our study provides comprehensive insights into essential aspects of ABC isolates originating from livestock-associated environments and clinical settings. We explored their resistance mechanisms and potential implications for public health, providing valuable knowledge for addressing these critical issues.

Airborne antibiotic resistome and microbiome in pharmaceutical factories

Antimicrobial resistance is considered to be one of the biggest public health problems, and airborne transmission is an important but under-appreciated pathway for the spread of antibiotic resistance genes (ARGs) in the environment. Previous research has shown pharmaceutical factories to be a major source of ARGs and antibiotic resistant bacteria (ARB) in the surrounding receiving water and soil environments. Pharmaceutical factories are hotspots of antibiotic resistance, but the atmospheric transmission and its environmental risk remain more concerns. Here, we conducted a metagenomic investigation into the airborne microbiome and resistome in three pharmaceutical factories in China. Soil (average: 38.45%) and wastewater (average: 28.53%) were major contributors of airborne resistome. ARGs (vanR/vanS, blaOXA, and CfxA) conferring resistance to critically important clinically used antibiotics were identified in the air samples. The wastewater treatment area had significantly higher relative abundances of ARGs (average: 0.64 copies/16S rRNA). Approximately 28.2% of the detected airborne ARGs were found to be associated with plasmids, and this increased to about 50% in the wastewater treatment area. We have compiled a list of high-risk airborne ARGs found in pharmaceutical factories. Moreover, A total of 1,043 viral operational taxonomic units were identified and linked to 47 family-group taxa. Different CRISPR-Cas immune systems have been identified in bacterial hosts in response to phage infection. Similarly, higher phage abundance (average: 2451.70 PPM) was found in the air of the wastewater treatment area. Our data provide insights into the antibiotic resistance gene profiles and microbiome (bacterial and non-bacterial) in pharmaceutical factories and reveal the potential role of horizontal transfer in the spread of airborne ARGs, with implications for human and animal health.

Complete genome sequences of 17 Salmonella enterica serovar Schwarzengrund isolates carrying an IncFIB-IncFIC (FII) fusion plasmid

Seventeen Salmonella enterica serovar Schwarzengrund isolates from chicken (n = 9) and clinical samples including stool (n = 6), urine (n = 1), and gallbladder (n = 1) were sequenced and found to carry an IncFIB-IncFIC (FII) fusion plasmid of approximately 145 Kb. This information provides reference genomic data for comparative studies of S. Schwarzengrund pathogenicity and plasmid genetics.

Identification of knowledge gaps in whole-genome sequence analysis of multi-resistant thermotolerant Campylobacter spp

BACKGROUND

Campylobacter spp. is the most frequent cause of bacterial food-borne gastroenteritis and a high priority antibiotic resistant bacterium according to the World Health Organization (WHO). European monitoring of thermotolerant Campylobacter spp. does not reflect the global burden of resistances already circulating within the bacterial population worldwide.

METHODS

We systematically compared whole genome sequencing with comprehensive phenotypic antimicrobial susceptibility, analyzing 494 thermotolerant Campylobacter poultry isolates from Vietnam and Germany. Any discrepancy was checked by repeating the wet lab and improving the dry lab part. Selected isolates were additionally analyzed via long-read Oxford Nanopore technology, leading to closed chromosomes and plasmids.

RESULTS

Overall, 22 different resistance genes and gene variants (e. g. erm(B), aph(3')-IIIa, aph(2'')-If, catA, lnu(C), blaOXA, sat4) and point mutations in three distinct genes (gyrA, 23S rRNA, rpsL) associated with AMR were present in the Campylobacter isolates. Two AMR genes were missing in the database and one falsely associated with resistance. Bioinformatic analysis based on short-read data partly failed to identify tet(O) and aadE, when the genes were present as duplicate or homologous gene variants. Intriguingly, isolates also contained different determinants, redundantly conferring resistance to chloramphenicol, gentamicin, kanamycin, lincomycin and streptomycin. We found a novel tet(W) in tetracycline sensitive strains, harboring point mutations. Furthermore, analysis based on assemblies from short-read data was impaired to identify full length phase variable aad9, due to variations of the poly-C tract within the gene. The genetic determinant responsible for gentamicin resistance of one isolate from Germany could not be identified. GyrT86I, presenting the main determinant for (fluoro-)quinolone resistance led to a rare atypical phenotype of ciprofloxacin resistance but nalidixic acid sensitivity. Long-read sequencing predicted AMR genes were mainly located on the chromosome, and rarely on plasmids. Predictions from long- and short-read sequencing, respectively, often differed. AMR genes were often organized in multidrug resistance islands (MDRI) and partially located in proximity to transposase genes, suggesting main mobilization of resistance determinants is via natural transformation and transposition in Campylobacter.

CONCLUSIONS

The results of this study suggest that there is frequent resistance gene duplication, mosaicism, and mutation leading to gene variation and truncation in Campylobacter strains that have not been reported in previous studies and are missing from databases. Furthermore, there is a need for deciphering yet unknown resistance mechanisms and resistance spread in thermotolerant Campylobacter spp. that may pose a challenge to global food safety.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Genomic dissection of Escherichia marmotae provides insights into diversity and pathogenic potential

Anthropogenic activities enhance the interconnection of human, animal, and environmental habitats and drive the evolution and inter-niche transmission of bacteria. Clear identification of emerging bacteria and pathogen control is therefore a public health priority. In 2015, the novel Escherichia species Escherichia marmotae was assigned, but due to the lack of appropriate detection and typing technologies, the One Health impact of this species is still being unraveled. E. marmotae represents a missing link in the impact of Escherichia spp. Here, we report 25 E. marmotae identified by next-generation sequencing that were previously phenotypically characterized as Escherichia coli during national zoonosis monitoring of food-producing animals. Applying fastANI to 153 738 published Escherichia spp. genome assemblies, we identified further 124 E. marmotae, originally classified as E. coli. Phylogenomics of all 149 isolates reveals an undefined population structure that is independent of the ecological niche. We highlight the phenotypic, genomic, and plasmid diversity of E. marmotae and provide evidence for gene flow across the species. The latter is illustrated by the acquisition of antibiotic resistance plasmids and pathogenicity islands, such as the type III secretion system. Thus, our comprehensive genomic overview of an emerging potential opportunistic pathogen underlines the importance of improved detection and characterization.

Differential responses of the gut microbiome and resistome to antibiotic exposures in infants and adults

Despite their crucial importance for human health, there is still relatively limited knowledge on how the gut resistome changes or responds to antibiotic treatment across ages, especially in the latter case. Here, we use fecal metagenomic data from 662 Danish infants and 217 young adults to fill this gap. The gut resistomes are characterized by a bimodal distribution driven by E. coli composition. The typical profile of the gut resistome differs significantly between adults and infants, with the latter distinguished by higher gene and plasmid abundances. However, the predominant antibiotic resistance genes (ARGs) are the same. Antibiotic treatment reduces bacterial diversity and increased ARG and plasmid abundances in both cohorts, especially core ARGs. The effects of antibiotic treatments on the gut microbiome last longer in adults than in infants, and different antibiotics are associated with distinct impacts. Overall, this study broadens our current understanding of gut resistome dynamics and the impact of antibiotic treatment across age groups.

Complete genome sequence of Methylomonas sp. UP202 isolated from an urban waterway sediment

We report the complete genome sequence of Methylomonas sp. UP202 isolated from an urban waterway sediment in Singapore. The genome contains genes involved in methane, methanol, formaldehyde, and formate oxidation. It also contains genes utilizing various nitrogen sources such as nitrogen, nitrate, nitrite, urea, and ammonium.

Genomic epidemiological investigation of an outbreak of Serratia marcescens neurosurgical site infections associated with contaminated haircutting toolkits in a hospital barber shop

BACKGROUND

Nine surgical site infections caused by Serratia marcescens were diagnosed in neurosurgical patients in a 3500-bed hospital between 2nd February and 6th April 2022.

OBJECTIVE

To trace the source of infections caused by S. marcescens to expedite termination of the outbreak and prevent future epidemics.

METHODS

A review of all surgical procedures and cultures yielding S. marcescens since February 2022 was conducted. Samples were collected from patients and environmental sources. S. marcescens isolates were characterized by antibiotic susceptibility testing. Whole-genome sequencing (WGS) was used to investigate genetic relationships. Resistance genes, virulence genes and plasmid replicons were identified.

RESULTS

S. marcescens was isolated from patients' puncture fluid, cerebrospinal fluid and other secretions, and was also cultured from the barbers' haircutting tools, including leather knives, slicker scrapers and razors. In total, 15 isolates were obtained from patients and eight isolates were obtained from haircutting tools. All isolates exhibited identical antibiotic resistance patterns. WGS revealed close clustering among the 23 isolates which differed significantly from previous strains. Three resistance genes and nine virulence-associated genes were detected in all isolates, and 19 of 23 isolates harboured an MOBP-type plasmid. The results confirmed an outbreak of S. marcescens, which was traced to contaminated haircutting tools in the hospital barber shop. The outbreak ended after extensive reinforcement of infection control procedures and re-education of the barbers.

CONCLUSIONS

These results highlight the risk of postoperative infections related to pre-operative skin preparation, and demonstrate the value of next-generation sequencing tools to expedite outbreak investigations.

Genetic characterization of a multidrug-resistant Salmonella enterica serovar Agona isolated from a dietary supplement in Germany

Salmonella enterica subsp. enterica serovar Agona has a history of causing food-borne outbreaks and any emergence of multidrug-resistant (MDR) isolates in novel food products is of concern. Particularly, in food products frequently consumed without sufficient heating prior to consumption. Here, we report about the MDR isolate, 18-SA00377, which had been isolated from a dietary supplement in Germany in 2018 and submitted to the German National Reference Laboratory for Salmonella. WGS-based comparative genetic analyses were conducted to find a potential reservoir of the isolate itself or mobile genetic elements associated with MDR. As a phylogenetic analysis did not yield any closely related S. Agona isolates, either globally or from Germany, a detailed analysis of the largest plasmid (295,499 bp) was performed as it is the main carrier of resistances. A combined approach of long-read and short-read sequencing enabled the assembly of the isolate's chromosome and its four plasmids. Their characterization revealed the presence of 23 different antibiotic resistance genes (ARGs), conferring resistance to 12 different antibiotic drug classes, as well as genes conferring resistance to six different heavy metals. The largest plasmid, pSE18-SA00377-1, belongs to the IncHI2 plasmid family and carries 16 ARGs, that are organized as two distinct clusters, with each ARG associated with putative composite transposons. Through a two-pronged approach, highly similar plasmids to pSE18-SA00377-1 were identified in the NCBI database and a search for Salmonella isolates with a highly similar ARG resistance profile was conducted. Mapping and structural comparisons between pSE18-SA00377-1 and these plasmids and Salmonella isolates showed that both the plasmid backbone and identical or similar ARG clusters can be found not only in Salmonella isolates, originating mostly from a wide variety of livestock, but also in a diverse range of bacterial genera of varying geographical origins and isolation sources. Thus, it can be speculated that the host range of pSE18-SA00377-1 is not restricted to Salmonella and its spread already occurred in different bacterial populations. Overall, this hints at a complex history for pSE18-SA00377-1 and highlights the importance of surveilling multidrug-resistant S. enterica isolates, especially in novel food items that are not yet heavily regulated.

Large-scale genomic survey and characterization of mcr genes carried by foodborne Cronobacter isolates

IMPORTANCE

Cronobacter is an emerging foodborne opportunistic pathogen, which can cause neonatal meningitis, bacteremia, and NEC by contaminating food. However, the entire picture of foodborne Cronobacter carriage of the mcr genes is not known. Here, we investigated the mcr genes of Cronobacter isolates by whole-genome sequencing and found 133 previously undescribed Cronobacter isolates carrying mcr genes. Further genomic analysis revealed that these mcr genes mainly belonged to the mcr-9 and mcr-10. Genomic analysis of the flanking structures of mcr genes revealed that two core flanking structures were prevalent in foodborne Cronobacter isolates, and the flanking structure carrying IS1R was found for the first time in this study.

MetaCC allows scalable and integrative analyses of both long-read and short-read metagenomic Hi-C data

Metagenomic Hi-C (metaHi-C) can identify contig-to-contig relationships with respect to their proximity within the same physical cell. Shotgun libraries in metaHi-C experiments can be constructed by next-generation sequencing (short-read metaHi-C) or more recent third-generation sequencing (long-read metaHi-C). However, all existing metaHi-C analysis methods are developed and benchmarked on short-read metaHi-C datasets and there exists much room for improvement in terms of more scalable and stable analyses, especially for long-read metaHi-C data. Here we report MetaCC, an efficient and integrative framework for analyzing both short-read and long-read metaHi-C datasets. MetaCC outperforms existing methods on normalization and binning. In particular, the MetaCC normalization module, named NormCC, is more than 3000 times faster than the current state-of-the-art method HiCzin on a complex wastewater dataset. When applied to one sheep gut long-read metaHi-C dataset, MetaCC binning module can retrieve 709 high-quality genomes with the largest species diversity using one single sample, including an expansion of five uncultured members from the order Erysipelotrichales, and is the only binner that can recover the genome of one important species Bacteroides vulgatus. Further plasmid analyses reveal that MetaCC binning is able to capture multi-copy plasmids.

plASgraph2: using graph neural networks to detect plasmid contigs from an assembly graph

Identification of plasmids from sequencing data is an important and challenging problem related to antimicrobial resistance spread and other One-Health issues. We provide a new architecture for identifying plasmid contigs in fragmented genome assemblies built from short-read data. We employ graph neural networks (GNNs) and the assembly graph to propagate the information from nearby nodes, which leads to more accurate classification, especially for short contigs that are difficult to classify based on sequence features or database searches alone. We trained plASgraph2 on a data set of samples from the ESKAPEE group of pathogens. plASgraph2 either outperforms or performs on par with a wide range of state-of-the-art methods on testing sets of independent ESKAPEE samples and samples from related pathogens. On one hand, our study provides a new accurate and easy to use tool for contig classification in bacterial isolates; on the other hand, it serves as a proof-of-concept for the use of GNNs in genomics. Our software is available at https://github.com/cchauve/plasgraph2 and the training and testing data sets are available at https://github.com/fmfi-compbio/plasgraph2-datasets.

Scalable and versatile container-based pipelines for de novo genome assembly and bacterial annotation

Background: Advancements in DNA sequencing technology have transformed the field of bacterial genomics, allowing for faster and more cost effective chromosome level assemblies compared to a decade ago. However, transforming raw reads into a complete genome model is a significant computational challenge due to the varying quality and quantity of data obtained from different sequencing instruments, as well as intrinsic characteristics of the genome and desired analyses. To address this issue, we have developed a set of container-based pipelines using Nextflow, offering both common workflows for inexperienced users and high levels of customization for experienced ones. Their processing strategies are adaptable based on the sequencing data type, and their modularity enables the incorporation of new components to address the community's evolving needs. Methods: These pipelines consist of three parts: quality control, de novo genome assembly, and bacterial genome annotation. In particular, the genome annotation pipeline provides a comprehensive overview of the genome, including standard gene prediction and functional inference, as well as predictions relevant to clinical applications such as virulence and resistance gene annotation, secondary metabolite detection, prophage and plasmid prediction, and more. Results: The annotation results are presented in reports, genome browsers, and a web-based application that enables users to explore and interact with the genome annotation results. Conclusions: Overall, our user-friendly pipelines offer a seamless integration of computational tools to facilitate routine bacterial genomics research. The effectiveness of these is illustrated by examining the sequencing data of a clinical sample of Klebsiella pneumoniae.

Diverse bacteria colonizing leaves and the rhizosphere of lettuce degrade azoxystrobin

Concerns about the possible effects of pesticide residues on both the environment and human health have increased worldwide. Bioremediation by the use of microorganisms to degrade or remove these residues has emerged as a powerful technology. However, the knowledge about the potential of different microorganisms for pesticide degradation is limited. This study focused on the isolation and characterisation of bacterial strains with the potential to degrade the active fungicide ingredient azoxystrobin. Potential degrading bacteria were tested in vitro and in the greenhouse, and the genomes of the best degrading strains were sequenced and analysed. We identified and characterised 59 unique bacterial strains, which were further tested in vitro and in greenhouse trials for their degradation activity. The best degraders from a foliar application trial in the greenhouse were identified as Bacillus subtilis strain MK101, Pseudomonas kermanshahensis strain MK113 and Rhodococcus fascians strain MK144 and analysed by whole genome sequencing. Genome analysis revealed that these three bacterial strains encode several genes predicted to be involved in the degradation of pesticides e.g., benC, pcaG, pcaH, however we could not find any specific gene previously reported to be involved in azoxystrobin degradation e.g., strH. Genome analysis pinpointed to some potential activities involved in plant growth promotion.

Comparative genomic analysis of Fusobacterium nucleatum reveals high intra-species diversity and cgmlst marker construction

BACKGROUND

Fusobacterium nucleatum is a one of the most important anaerobic opportunistic pathogens in the oral and intestinal tracts of human and animals. It can cause various diseases such as infections, Lemierre's syndrome, oral cancer and colorectal cancer. The comparative genomic studies on the population genome level, have not been reported.

RESULTS

We analyzed all publicly available Fusobacterium nucleatums' genomic data for a comparative genomic study, focusing on the pan-genomic features, virulence genes, plasmid genomes and developed cgmlst molecular markers. We found the pan-genome shows a clear open tendency and most of plasmids in Fusobacterium nucleatum are mainly transmitted intraspecifically.

CONCLUSIONS

Our comparative analysis of Fusobacterium nucleatum systematically revealed the open pan-genomic features and phylogenetic tree based on cgmlst molecular markers. What's more, we also identified common plasmid typing among genomes. We hope that our study will provide a theoretical basis for subsequent functional studies.

Intensive care unit sinks are persistently colonized with multidrug resistant bacteria and mobilizable, resistance-conferring plasmids

Contamination of hospital sinks with microbial pathogens presents a serious potential threat to patients, but our understanding of sink colonization dynamics is largely based on infection outbreaks. Here, we investigate the colonization patterns of multidrug-resistant organisms (MDROs) in intensive care unit sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. Using culture-based methods, we recovered 822 bacterial isolates representing 104 unique species and genomospecies. Genomic analyses revealed long-term colonization by Pseudomonas spp. and Serratia marcescens strains across multiple rooms. Nanopore sequencing uncovered examples of long-term persistence of resistance-conferring plasmids in unrelated hosts. These data indicate that antibiotic resistance (AR) in Pseudomonas spp. is maintained both by strain colonization and horizontal gene transfer (HGT), while HGT maintains AR within Acinetobacter spp. and Enterobacterales, independent of colonization. These results emphasize the importance of proactive, genomic-focused surveillance of built environments to mitigate MDRO spread. IMPORTANCE Hospital sinks are frequently linked to outbreaks of antibiotic-resistant bacteria. Here, we used whole-genome sequencing to track the long-term colonization patterns in intensive care unit (ICU) sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. We analyzed 822 bacterial genomes, representing over 100 different species. We identified long-term contamination by opportunistic pathogens, as well as transient appearance of other common pathogens. We found that bacteria recovered from the ICU had more antibiotic resistance genes (ARGs) in their genomes compared to matched community spaces. We also found that many of these ARGs are harbored on mobilizable plasmids, which were found shared in the genomes of unrelated bacteria. Overall, this study provides an in-depth view of contamination patterns for common nosocomial pathogens and identifies specific targets for surveillance.

Psychrotrophic Bacteria Equipped with Virulence and Colonization Traits Populate the Ice Cream Manufacturing Environment

Several microbial taxa have been associated with food processing facilities, and they might resist by attaching on tools and equipment even after sanitation procedures, producing biofilms that adhere to the surfaces and might embed other microorganisms, including spoilers and pathogens. There is increasing evidence that these communities can be transferred to the final product. To explore the microbial contamination routes in a facility producing ice creams, we collected foods and environmental swabs from industrial surfaces of equipment and tools and performed taxonomic and functional analyses of the microbial DNA extracted from the environmental samples. Our results suggest that complex communities dominated by psychrotrophic bacteria (e.g., Pseudomonas and Acinetobacter spp.) inhabit the food processing environment, and we demonstrate that these communities might be transferred from the surfaces to the products. Functional analysis performed on environmental samples highlighted the presence of several genes linked to antimicrobial resistance and adherence on abiotic surfaces; such genes were more abundant on food contact (FC) than on other surfaces. Metagenome-assembled genomes (MAGs) of Pseudomonas stutzeri showed genes linked with biofilm formation and motility, which are surely linked to colonizing capabilities in the processing lines. The study highlights clear potential advantages of applying microbiome mapping in the food industry for source tracking of microbial contamination and for planning appropriate ad hoc sanitization strategies. IMPORTANCE Several microbial species might permanently establish in food processing facilities, thus contributing to food loss. In fact, food contact surfaces might transfer microorganisms to intermediates and products, potentially representing a hazard to human health. In this work, we provide evidence of the existence of complex microbial communities overcoming sanitation in an ice cream-producing facility. These communities harbored several genes that could potentially lead to attachment to surfaces and antimicrobial resistance. Also, prediction of routes of contamination showed that several potential spoilage taxa might end up in the final product. Importantly, in this work, we show that mapping the environmental microbiome is a high-resolution technique that might help food business operators ensure food quality and safety through detection of potentially hazardous microorganisms.

Whole genome sequence-based analysis of Staphylococcus aureus isolated from bovine mastitis in Thuringia, Germany

Background

Bovine mastitis is a common disease of dairy cattle causing major economic losses due to reduced yield and poor quality of milk worldwide. The current investigation aimed to gain insight into the genetic diversity, antimicrobial resistance profiles and virulence associated factors of Staphylococcus (S.) aureus isolated from clinical bovine mastitis in dairy farms in Thuringia, Germany.

Methods

Forty Staphylococcus aureus isolates collected from clinical bovine mastitis cases from 17 Thuringian dairy farms were phenotyped and genetically characterized using whole genome sequencing.

Results

Out of 40 S. aureus, 30 (75%) were confirmed as methicillin resistant isolates. The isolates showed elevated antimicrobial resistance against penicillin, tetracycline and oxacillin, i.e., 77.5, 77.5, and 75%, respectively. Lower resistance rates were found against moxifloxacin, ciprofloxacin, gentamicin and trimethoprim/sulfamethoxazole, i.e., 35, 35, 30, and 22.5%, respectively. While resistance against clindamycin and erythromycin was rarely found (5 and 2.5%, respectively). All isolates were susceptible to linezolid, teicoplanin, vancomycin, tigecycline, fosfomycin, fusidic acid and rifampicin. These isolates were further allocated into five different sequence types: ST398 (n = 31), ST1074 (n = 4), ST504 (n = 3), ST582 (CC15) (n = 1) and ST479 (n = 1). These isolates were also assigned to seven clusters with up to 100 SNP which has facilitated geographical mapping and epidemiological distribution in Thuringia. Strains belonging to ST398 were classified into clusters 1, 2, 3, 4 and 7. The isolates of ST504 were of cluster 5, those of ST1074 were belonging to cluster 6. Resistance genes blaZ, blaI and blaR associated with penicillin resistance were found in 32 (80%) strains, all except one were belonging to ST398. Methicillin resistance associated mecA was identified in 30 (96.8%) isolates of ST398. All tetracycline and erythromycin resistant isolates were of ST398, and all harbored both tetM and ermA. About 90.3% of tetracycline resistant isolates assigned to ST398 were also carrying tetK gene. The point mutations parC_S80F, gyrA_S84L and parC_S80Y in gyrA and parC associated with quinolone resistance were found in all phenotypically resistant isolates to ciprofloxacin and moxifloxacin (n = 14). Sixty-eight virulence genes were identified among isolates. Both lukD/E and lukM/F-PV-P83 were identified in 22.5% of isolates, all were non-ST398.

Conclusion

In this study, ST398 had the highest potential to cause disease and had a massive prevalence in bovine mastitis cases. Five different sequence types and seven clusters were identified in the federal state of Thuringia. The circulation of some clusters in the same region over several years shows the persistence of cluster-associated infection despite the intensive veterinary care. On the other hand, some regions had different clusters at the same year or in different consecutive years. Different sequence types and associated different clusters of S. aureus were geographically widely distributed among dairy farms in Thuringia. The findings of this study show that various clusters have the potential to spread over a large geographical scale. The detection of LA-MRSA on dairy farms, which is known for cabapility to widely spread among different groups of animals, humans and their environment urges for the implementation of national wide strategic programs. The identification of CA-MRSA among the isolates such as ST398 poses a significant risk for the transmission of such strains between animals and humans on dairy farms.

Insights into the physiological and genomic characterization of three bacterial isolates from a highly alkaline, terrestrial serpentinizing system

The terrestrial serpentinite-hosted ecosystem known as "The Cedars" is home to a diverse microbial community persisting under highly alkaline (pH ~ 12) and reducing (Eh < -550 mV) conditions. This extreme environment presents particular difficulties for microbial life, and efforts to isolate microorganisms from The Cedars over the past decade have remained challenging. Herein, we report the initial physiological assessment and/or full genomic characterization of three isolates: Paenibacillus sp. Cedars ('Paeni-Cedars'), Alishewanella sp. BS5-314 ('Ali-BS5-314'), and Anaerobacillus sp. CMMVII ('Anaero-CMMVII'). Paeni-Cedars is a Gram-positive, rod-shaped, mesophilic facultative anaerobe that grows between pH 7-10 (minimum pH tested was 7), temperatures 20-40°C, and 0-3% NaCl concentration. The addition of 10-20 mM CaCl₂ enhanced growth, and iron reduction was observed in the following order, 2-line ferrihydrite > magnetite > serpentinite ~ chromite ~ hematite. Genome analysis identified genes for flavin-mediated iron reduction and synthesis of a bacillibactin-like, catechol-type siderophore. Ali-BS5-314 is a Gram-negative, rod-shaped, mesophilic, facultative anaerobic alkaliphile that grows between pH 10-12 and temperatures 10-40°C, with limited growth observed 1-5% NaCl. Nitrate is used as a terminal electron acceptor under anaerobic conditions, which was corroborated by genome analysis. The Ali-BS5-314 genome also includes genes for benzoate-like compound metabolism. Anaero-CMMVII remained difficult to cultivate for physiological studies; however, growth was observed between pH 9-12, with the addition of 0.01-1% yeast extract. Anaero-CMMVII is a probable oxygen-tolerant anaerobic alkaliphile with hydrogenotrophic respiration coupled with nitrate reduction, as determined by genome analysis. Based on single-copy genes, ANI, AAI and dDDH analyses, Paeni-Cedars and Ali-BS5-314 are related to other species (P. glucanolyticus and A. aestuarii, respectively), and Anaero-CMMVII represents a new species. The characterization of these three isolates demonstrate the range of ecophysiological adaptations and metabolisms present in serpentinite-hosted ecosystems, including mineral reduction, alkaliphily, and siderophore production.

Plasmer: an Accurate and Sensitive Bacterial Plasmid Prediction Tool Based on Machine Learning of Shared k-mers and Genomic Features

Identification of plasmids in bacterial genomes is critical for many factors, including horizontal gene transfer, antibiotic resistance genes, host-microbe interactions, cloning vectors, and industrial production. There are several in silico methods to predict plasmid sequences in assembled genomes. However, existing methods have evident shortcomings, such as unbalance in sensitivity and specificity, dependency on species-specific models, and performance reduction in sequences shorter than 10 kb, which has limited their scope of applicability. In this work, we proposed Plasmer, a novel plasmid predictor based on machine-learning of shared k-mers and genomic features. Unlike existing k-mer or genomic-feature based methods, Plasmer employs the random forest algorithm to make predictions using the percent of shared k-mers with plasmid and chromosome databases combined with other genomic features, including alignment E value and replicon distribution scores (RDS). Plasmer can predict on multiple species and has achieved an average the area under the curve (AUC) of 0.996 with accuracy of 98.4%. Compared to existing methods, tests of both sliding sequences and simulated and de novo assemblies have consistently shown that Plasmer has outperforming accuracy and stable performance across long and short contigs above 500 bp, demonstrating its applicability for fragmented assemblies. Plasmer also has excellent and balanced performance on both sensitivity and specificity (both >0.95 above 500 bp) with the highest F1-score, which has eliminated the bias on sensitivity or specificity that was common in existing methods. Plasmer also provides taxonomy classification to help identify the origin of plasmids. IMPORTANCE In this study, we proposed a novel plasmid prediction tool named Plasmer. Technically, unlike existing k-mer or genomic features-based methods, Plasmer is the first tool to combine the advantages of the percent of shared k-mers and the alignment score of genomic features. This has given Plasmer (i) evident improvement in performance compared to other methods, with the best F1-score and accuracy on sliding sequences, simulated contigs, and de novo assemblies; (ii) applicability for contigs above 500 bp with highest accuracy, enabling plasmid prediction in fragmented short-read assemblies; (iii) excellent and balanced performance between sensitivity and specificity (both >0.95 above 500 bp) with the highest F1-score, which eliminated the bias on sensitivity or specificity that commonly existed in other methods; and (iv) no dependency of species-specific training models. We believe that Plasmer provides a more reliable alternative for plasmid prediction in bacterial genome assemblies.

Oxford nanopore technologies-a valuable tool to generate whole-genome sequencing data for in silico serotyping and the detection of genetic markers in Salmonella

Bacteria of the genus Salmonella pose a major risk to livestock, the food economy, and public health. Salmonella infections are one of the leading causes of food poisoning. The identification of serovars of Salmonella achieved by their diverse surface antigens is essential to gain information on their epidemiological context. Traditionally, slide agglutination has been used for serotyping. In recent years, whole-genome sequencing (WGS) followed by in silico serotyping has been established as an alternative method for serotyping and the detection of genetic markers for Salmonella. Until now, WGS data generated with Illumina sequencing are used to validate in silico serotyping methods. Oxford Nanopore Technologies (ONT) opens the possibility to sequence ultra-long reads and has frequently been used for bacterial sequencing. In this study, ONT sequencing data of 28 Salmonella strains of different serovars with epidemiological relevance in humans, food, and animals were taken to investigate the performance of the in silico serotyping tools SISTR and SeqSero2 compared to traditional slide agglutination tests. Moreover, the detection of genetic markers for resistance against antimicrobial agents, virulence, and plasmids was studied by comparing WGS data based on ONT with WGS data based on Illumina. Based on the ONT data from flow cell version R9.4.1, in silico serotyping achieved an accuracy of 96.4 and 92% for the tools SISTR and SeqSero2, respectively. Highly similar sets of genetic markers comparing both sequencing technologies were identified. Taking the ongoing improvement of basecalling and flow cells into account, ONT data can be used for Salmonella in silico serotyping and genetic marker detection.