Outcome of Different Sequencing and Assembly Approaches on the Detection of Plasmids and Localization of Antimicrobial Resistance Genes in Commensal Escherichia coli

Single Laboratory Evaluation of the Q20+ Nanopore Sequencing Kit for Bacterial Outbreak Investigations

Leafy greens are a significant source of produce-related Shiga toxin-producing Escherichia coli (STEC) outbreaks in the United States, with agricultural water often implicated as a potential source. Current FDA outbreak detection protocols are time-consuming and rely on sequencing methods performed in costly equipment. This study evaluated the potential of Oxford Nanopore Technologies (ONT) with Q20+ chemistry as a cost-effective, rapid, and accurate method for identifying and clustering foodborne pathogens. The study focuses on assessing whether ONT Q20+ technology could facilitate near real-time pathogen identification, including SNP differences, serotypes, and antimicrobial resistance genes. This pilot study evaluated different combinations of two DNA extraction methods (Maxwell RSC Cultured Cell DNA kit and Monarch high molecular weight extraction kits) and two ONT library preparation protocols (ligation and the rapid barcoding sequencing kit) using five well-characterized strains representing diverse foodborne pathogens. High-quality, closed bacterial genomes were obtained from all combinations of extraction and sequencing kits. However, variations in assembly length and genome completeness were observed, indicating the need for further optimization. In silico analyses demonstrated that Q20+ nanopore sequencing chemistry accurately identified species, genotype, and virulence factors, with comparable results to Illumina sequencing. Phylogenomic clustering showed that ONT assemblies clustered with reference genomes, though some indels and SNP differences were observed, likely due to sequencing and analysis methodologies rather than inherent genetic variation. Additionally, the study evaluated the impact of a change in the sampling rates from 4 kHz (260 bases pair second) to 5 kHz (400 bases pair second), finding no significant difference in sequencing accuracy. This evaluation workflow offers a framework for evaluating novel technologies for use in surveillance and foodborne outbreak investigations. Overall, the evaluation demonstrated the potential of ONT Q20+ nanopore sequencing chemistry to assist in identifying the correct strain during outbreak investigations. However, further research, validation studies, and optimization efforts are needed to address the observed limitations and fully realize the technology's potential for improving public health outcomes and enabling more efficient responses to foodborne disease threats.

Widespread prevalence of plasmid-mediated blaCTX-M type extended-spectrum beta-lactamase Escherichia coli in backyard broiler production systems in the United States

Extended-spectrum beta-lactamase (ESBL) Escherichia coli (E. coli) is an emerging pathogen of high concern given its resistance to extended-spectrum cephalosporins. Broiler chicken, which is the number one consumed meat in the United States and worldwide, can be a reservoir of ESBL E. coli. Backyard poultry ownership is on the rise in the United States, yet there is little research investigating prevalence of ESBL E. coli in this setting. This study aims to identify the prevalence and antimicrobial resistance profiles (phenotypically and genotypically) of ESBL E. coli in some backyard and commercial broiler farms in the U.S. For this study ten backyard and ten commercial farms were visited at three time-points across flock production. Fecal (n = 10), litter/compost (n = 5), soil (n = 5), and swabs of feeders and waterers (n = 6) were collected at each visit and processed for E. coli. Assessment of ESBL phenotype was determined through using disk diffusion with 3rd generation cephalosporins, cefotaxime and ceftazidime, and that with clavulanic acid. Broth microdilution and whole genome sequencing were used to investigate both phenotypic and genotypic resistance profiles, respectively. ESBL E. coli was more prevalent in backyard farms with 12.95% of samples testing positive whereas 0.77% of commercial farm samples were positive. All isolates contained a blaCTX-M gene, the dominant variant being blaCTX-M-1, and its presence was entirely due to plasmids. Our study confirms concerns of growing resistance to fourth generation cephalosporin, cefepime, as roughly half (51.4%) of all isolates were found to be susceptible dose-dependent and few were resistant. Resistance to non-beta lactams, gentamicin and ciprofloxacin, was also detected in our samples. Our study identifies prevalence of blaCTX-M type ESBL E. coli in U.S. backyard broiler farms, emphasizing the need for interventions for food and production safety.

Sensitivity and consistency of long- and short-read metagenomics and epicPCR for the detection of antibiotic resistance genes and their bacterial hosts in wastewater

Wastewater surveillance is a powerful tool to assess the risks associated with antibiotic resistance in communities. One challenge is selecting which analytical tool to deploy to measure risk indicators, such as antibiotic resistance genes (ARGs) and their respective bacterial hosts. Although metagenomics is frequently used for analyzing ARGs, few studies have compared the performance of long-read and short-read metagenomics in identifying which bacteria harbor ARGs in wastewater. Furthermore, for ARG host detection, untargeted metagenomics has not been compared to targeted methods such as epicPCR. Here, we 1) evaluated long-read and short-read metagenomics as well as epicPCR for detecting ARG hosts in wastewater, and 2) investigated the host range of ARGs across the wastewater treatment plant (WWTP) to evaluate host proliferation. Results highlighted long-read revealed a wider range of ARG hosts compared to short-read metagenomics. Nonetheless, the ARG host range detected by long-read metagenomics only represented a subset of the hosts detected by epicPCR. The ARG-host linkages across the influent and effluent of the WWTP were characterized. Results showed the ARG-host phylum linkages were relatively consistent across the WWTP, whereas new ARG-host species linkages appeared in the WWTP effluent. The ARG-host linkages of several clinically relevant species found in the effluent were identified.

Do we still need Illumina sequencing data? Evaluating Oxford Nanopore Technologies R10.4.1 flow cells and the Rapid v14 library prep kit for Gram negative bacteria whole genome assemblies

The best whole genome assemblies are currently built from a combination of highly accurate short-read sequencing data and long-read sequencing data that can bridge repetitive and problematic regions. Oxford Nanopore Technologies (ONT) produce long-read sequencing platforms and they are continually improving their technology to obtain higher quality read data that is approaching the quality obtained from short-read platforms such as Illumina. As these innovations continue, we evaluated how much ONT read coverage produced by the Rapid Barcoding Kit v14 (SQK-RBK114) is necessary to generate high-quality hybrid and long-read-only genome assemblies for a panel of carbapenemase-producing Enterobacterales bacterial isolates. We found that 30× long-read coverage is sufficient if Illumina data are available, and that more (at least 100× long-read coverage is recommended for long-read-only assemblies. Illumina polishing is still improving single nucleotide variants (SNVs) and INDELs in long-read-only assemblies. We also examined if antimicrobial resistance genes could be accurately identified in long-read-only data, and found that Flye assemblies regardless of ONT coverage detected >96% of resistance genes at 100% identity and length. Overall, the Rapid Barcoding Kit v14 and long-read-only assemblies can be an optimal sequencing strategy (i.e., plasmid characterization and AMR detection) but finer-scale analyses (i.e., SNV) still benefit from short-read data.

Hybrid Illumina-Nanopore assembly improves identification of multilocus sequence types and antimicrobial resistance genes of Staphylococcus aureus isolated from Vermont dairy farms: comparison to Illumina-only and R9.4.1 nanopore-only assemblies

Antimicrobial resistance (AMR) in Staphylococcus aureus is a pressing public health challenge with significant implications for the dairy industry, encompassing bovine mastitis concerns and potential zoonotic threats. To delve deeper into the resistance mechanisms of S. aureus, this study employed a hybrid whole genome assembly approach that synergized the precision of Illumina with the continuity of Oxford Nanopore. A total of 62 isolates, collected from multiple sources from Vermont dairy farms, were sequenced using the GridION Oxford Nanopore R9.4.1 platform and the Illumina platform, and subsequently processed through our long-read first bioinformatics pipeline. Our analyses showcased the hybrid-assembled genome's superior completeness compared to Oxford Nanopore (R9.4.1)-only or Illumina-only assembled genomes. Furthermore, the hybrid assembly accurately determined multilocus sequence typing (MLST) strain types across all isolates. The comprehensive probe for antibiotic resistance genes (ARGs) using databases like CARD, Resfinder, and MEGARES 2.0 characterized AMR in S. aureus isolates from Vermont dairy farms, and revealed the presence of notable resistance genes, including beta-lactam genes blaZ, blaI, and blaR. In conclusion, the hybrid assembly approach emerged as a tool for uncovering the genomic nuances of S. aureus isolates collected from multiple sources on dairy farms. Our findings offer a pathway for detecting AMR gene prevalence and shaping AMR management strategies crucial for safeguarding human and animal health.

Application of phylodynamics to identify spread of antimicrobial-resistant Escherichia coli between humans and canines in an urban environment

The transmission of antimicrobial resistant bacteria in the urban environment is poorly understood. We utilized genomic sequencing and phylogenetics to characterize the transmission dynamics of antimicrobial resistant Escherichia coli (AMR-Ec) cultured from putative canine (caninep) and human feces present on urban sidewalks in San Francisco, California. We isolated a total of fifty-six AMR-Ec isolates from human (n = 20) and caninep (n = 36) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antimicrobial resistance (AMR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and caninesp from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Our results provide evidence for multiple sharing events of AMR-Ec between humans and caninesp. In particular, we found one instance of likely transmission from caninesp to humans as well as an additional local outbreak cluster consisting of one caninep and one human sample. Based on this analysis, it appears that non-human feces act as an important reservoir of clinically relevant AMR-Ec within the urban environment for this study population. This work showcases the utility of genomic epidemiology to reconstruct potential pathways by which antimicrobial resistance spreads.

Plasmid-encoded gene duplications of extended-spectrum β-lactamases in clinical bacterial isolates

Introduction

The emergence of extended-spectrum β-lactamase (ESBL)-producing Enterobacteriaceae is an urgent and alarming One Health problem. This study aimed to investigate duplications of plasmid-encoded ESBL genes and their impact on antimicrobial resistance (AMR) phenotypes in clinical and screening isolates.

Methods

Multi-drug-resistant bacteria from hospitalized patients were collected during routine clinical surveillance from January 2022 to June 2023, and their antimicrobial susceptibility patterns were determined. Genotypes were extracted from long-read whole-genome sequencing data. Furthermore, plasmids and other mobile genetic elements associated with ESBL genes were characterized, and the ESBL genes were correlated to ceftazidime minimal inhibitory concentration (MIC).

Results

In total, we identified four cases of plasmid-encoded ESBL gene duplications that match four genetically similar plasmids during the 18-month surveillance period: five Escherichia coli and three Klebsiella pneumoniae isolates. As the ESBL genes were part of transposable elements, the surrounding sequence regions were duplicated as well. In-depth analysis revealed insertion sequence (IS)-mediated transposition mechanisms. Isolates with duplicated ESBL genes exhibited a higher MIC for ceftazidime in comparison to isolates with a single gene copy (3-256 vs. 1.5-32 mg/L, respectively).

Conclusion

ESBL gene duplications led to an increased phenotypic resistance against ceftazidime. Our data suggest that ESBL gene duplications by an IS-mediated transposition are a relevant mechanism for how AMR develops in the clinical setting and is part of the microevolution of plasmids.

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.

Exploring Nocardia's ecological spectrum and novel therapeutic frontiers through whole-genome sequencing: unraveling drug resistance and virulence factors

Nocardia farcinica is the leading pathogen responsible for nocardiosis, a life-threatening infection primarily affecting immunocompromised patients. In this study, the genomic sequence of a clinically isolated N. farcinica sample was sequenced. Subsequently, the assembled genome was annotated to identify antimicrobial resistance and virulence genes, as well as plasmid and prophages. The analysis of the entire genome size was 6,021,225 bp, with a GC content of 70.78% and consists of 103 contigs and N50 values of 292,531 bp. The genome analysis revealed the presence of several antimicrobial resistance genes, including RbpA, mtrA, FAR-1, blaFAR-1, blaFAR-1_1, and rox. In addition, virulence genes such as relA, icl, and mbtH were also detected. The present study signifies that N. farcinica genome is pivotal for the understanding of antimicrobial resistance and virulence genes is crucial for comprehending resistance mechanism, and developing effective strategies to combat bacterial infections effectively, especially adhesins and toxins. This study aids in identifying crucial drug targets for combating multidrug-resistant N. farcinica in the future.

Three Rounds of Read Correction Significantly Improve Eukaryotic Protein Detection in ONT Reads

BACKGROUND

Eukaryotes' whole-genome sequencing is crucial for species identification, gene detection, and protein annotation. Oxford Nanopore Technology (ONT) is an affordable and rapid platform for sequencing eukaryotes; however, the relatively higher error rates require computational and bioinformatic efforts to produce more accurate genome assemblies. Here, we evaluated the effect of read correction tools on eukaryote genome completeness, gene detection and protein annotation.

METHODS

Reads generated by ONT of four eukaryotes, C. albicans, C. gattii, S. cerevisiae, and P. falciparum, were assembled using minimap2 and underwent three rounds of read correction using flye, medaka and racon. The generates consensus FASTA files were compared for total length (bp), genome completeness, gene detection, and protein-annotation by QUAST, BUSCO, BRAKER1 and InterProScan, respectively.

RESULTS

Genome completeness was dependent on the assembly method rather than on the read correction tool; however, medaka performed better than flye and racon. Racon significantly performed better than flye and medaka in gene detection, while both racon and medaka significantly performed better than flye in protein-annotation.

CONCLUSION

We show that three rounds of read correction significantly affect gene detection and protein annotation, which are dependent on assembly quality in preference to assembly completeness.

Role of natural transformation in the evolution of small cryptic plasmids in Synechocystis sp. PCC 6803

Small cryptic plasmids have no clear effect on the host fitness and their functional repertoire remains obscure. The naturally competent cyanobacterium Synechocystis sp. PCC 6803 harbours several small cryptic plasmids; whether their evolution with this species is supported by horizontal transfer remains understudied. Here, we show that the small cryptic plasmid DNA is transferred in the population exclusively by natural transformation, where the transfer frequency of plasmid-encoded genes is similar to that of chromosome-encoded genes. Establishing a system to follow gene transfer, we compared the transfer frequency of genes encoded in cryptic plasmids pCA2.4 (2378 bp) and pCB2.4 (2345 bp) within and between populations of two Synechocystis sp. PCC 6803 labtypes (termed Kiel and Sevilla). Our results reveal that plasmid gene transfer frequency depends on the recipient labtype. Furthermore, gene transfer via whole plasmid uptake in the Sevilla labtype ranged among the lowest detected transfer rates in our experiments. Our study indicates that horizontal DNA transfer via natural transformation is frequent in the evolution of small cryptic plasmids that reside in naturally competent organisms. Furthermore, we suggest that the contribution of natural transformation to cryptic plasmid persistence in Synechocystis is limited.

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.

Evaluation of metagenomic assembly methods for the detection and characterization of antimicrobial resistance determinants and associated mobilizable elements

Antimicrobial resistance genes (ARGs) can be transferred between members of a bacterial population by mobile genetic elements (MGE). Understanding the risk of these transfer events is important in monitoring and predicting antimicrobial resistance (AMR), especially in the context of a One Health Continuum. However, there is no universally accepted method for detection of ARGs and MGEs, and especially for determining their linkages. This study used publicly available shotgun metagenomic DNA short-read (Illumina, 100 bp paired-end) sequence data from samples across the One Health Continuum (including beef cattle composite feces from feedlots, catch basin water at feedlots, agricultural soil from feedlot manured surrounding fields, and urban/municipal sewage influent from two municipal wastewater treatment plants) to develop a workflow to identify and associate ARGs and MGEs. ARG- and MGE-based targeted-assemblies with available short-read data were unable to meet this analysis goal. In contrast, de novo assembly of contigs provided enough sequence context to associate ARGs and MGEs, without compromising discovery rate. However, to estimate the relative abundance of these elements, unassembled sequence data must still be used.

The impact of applying various de novo assembly and correction tools on the identification of genome characterization, drug resistance, and virulence factors of clinical isolates using ONT sequencing

Oxford Nanopore sequencing technology (ONT) is currently widely used due to its affordability, simplicity, and reliability. Despite the advantage ONT has over next-generation sequencing in detecting resistance genes in mobile genetic elements, its relatively high error rate (10-15%) is still a deterrent. Several bioinformatic tools are freely available for raw data processing and obtaining complete and more accurate genome assemblies. In this study, we evaluated the impact of using mix-and-matched read assembly (Flye, Canu, Wtdbg2, and NECAT) and read correction (Medaka, NextPolish, and Racon) tools in generating complete and accurate genome assemblies, and downstream genomic analysis of nine clinical Escherichia coli isolates. Flye and Canu assemblers were the most robust in genome assembly, and Medaka and Racon correction tools significantly improved assembly parameters. Flye functioned well in pan-genome analysis, while Medaka increased the number of core genes detected. Flye, Canu, and NECAT assembler functioned well in detecting antimicrobial resistance genes (AMR), while Wtdbg2 required correction tools for better detection. Flye was the best assembler for detecting and locating both virulence and AMR genes (i.e., chromosomal vs. plasmid). This study provides insight into the performance of several read assembly and read correction tools for analyzing ONT sequencing reads for clinical isolates.

Illustration on phenotypic and genotypic characteristics of typical multi-antibiotic resistant bacteria in aquatic environments through complete genomes and comparative genomics

Antibiotic-resistant bacteria, especially multi-antibiotic-resistant bacteria (MARBs), greatly threaten environmental safety and human health. However, studies on the phenotypic resistance and complete genotypic characterization of MARB in aquatic environments are lacking. In this study, a multi-resistant superbug (TR3) was screened by the selective pressure of multi-antibiotics from the activated sludge of the aeration tanks of urban wastewater treatment plants (WWTPs) in 5 different regions of China. Based on the 16 S rDNA sequence alignment it was found that the sequence similarity between strain TR3 and Aeromonas was as high as 99.50%. The genome-wide sequence showed that the base content of the chromosome of strain TR3 is 4,521,851 bp. It contains a plasmid with a length of 9182 bp. All antibiotic resistance genes (ARGs) of strain TR3 are located on the chromosome, which means that it has passage stability. There are multiple types of resistance genes in the genome and plasmid of strain TR3, enduing it with resistance to 5 antibiotics (ciprofloxacin, tetracycline, ampicillin, clarithromycin, and kanamycin), accompanied by the strongest resistance to kanamycin (aminoglycosides) and the worst resistance to clarithromycin (quinolones). From the perspective of gene expression, we show the resistance mechanism of strain TR3 to different types of antibiotics. In addition, the potential pathogenicity of strain TR3 is also discussed. Chlorine and ultraviolet (UV) sterilization on strain TR3 showed that UV is ineffective at low intensity, and it is easy to be revived by light. A low concentration of hypochlorous acid is effective for sterilization, but it can cause the release of DNA, becoming a potential source of ARGs discharged from WWTPs to environmental water bodies.

Phylodynamics Uncovers the Transmission of Antibiotic-Resistant Escherichia coli between Canines and Humans in an Urban Environment

The role of canines in transmitting antibiotic resistant bacteria to humans in the urban environment is poorly understood. To elucidate this role, we utilized genomic sequencing and phylogenetics to characterize the burden and transmission dynamics of antibiotic resistant Escherichia coli (ABR-Ec) cultured from canine and human feces present on urban sidewalks in San Francisco, California. We collected a total of fifty-nine ABR-Ec from human (n=12) and canine (n=47) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antibiotic resistance (ABR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and canines from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Overall, we found human and canine samples to carry similar amounts and profiles of ABR genes. Our results provide evidence for multiple transmission events of ABR-Ec between humans and canines. In particular, we found one instance of likely transmission from canines to humans as well as an additional local outbreak cluster consisting of one canine and one human sample. Based on this analysis, it appears that canine feces act as an important reservoir of clinically relevant ABR-Ec within the urban environment. Our findings support that public health measures should continue to emphasize proper canine feces disposal practices, access to public toilets and sidewalk and street cleaning. Importance: Antibiotic resistance in E. coli is a growing public health concern with global attributable deaths projected to reach millions annually. Current research has focused heavily on clinical routes of antibiotic resistance transmission to design interventions while the role of alternative reservoirs such as domesticated animals remain less well understood. Our results suggest canines are part of the transmission network that disseminates high-risk multidrug resistance in E. coli within the urban San Francisco community. As such, this study highlights the need to consider canines, and potentially domesticated animals more broadly, when designing interventions to reduce the prevalence of antibiotic resistance in the community. Additionally, it showcases the utility of genomic epidemiology to reconstruct the pathways by which antimicrobial resistance spreads.

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.

Phylogenetic Analysis of Escherichia coli Isolated from Australian Feedlot Cattle in Comparison to Pig Faecal and Poultry/Human Extraintestinal Isolates

The similarity of commensal Escherichia coli isolated from healthy cattle to antimicrobial-resistant bacteria causing extraintestinal infections in humans is not fully understood. In this study, we used a bioinformatics approach based on whole genome sequencing data to determine the genetic characteristics and phylogenetic relationships among faecal Escherichia coli isolates from beef cattle (n = 37) from a single feedlot in comparison to previously analysed pig faecal (n = 45), poultry extraintestinal (n = 19), and human extraintestinal E. coli isolates (n = 40) from three previous Australian studies. Most beef cattle and pig isolates belonged to E. coli phylogroups A and B1, whereas most avian and human isolates belonged to B2 and D, although a single human extraintestinal isolate belonged to phylogenetic group A and sequence type (ST) 10. The most common E. coli sequence types (STs) included ST10 for beef cattle, ST361 for pig, ST117 for poultry, and ST73 for human isolates. Extended-spectrum and AmpC β-lactamase genes were identified in seven out of thirty-seven (18.9%) beef cattle isolates. The most common plasmid replicons identified were IncFIB (AP001918), followed by IncFII, Col156, and IncX1. The results confirm that feedlot cattle isolates examined in this study represent a reduced risk to human and environmental health with regard to being a source of antimicrobial-resistant E. coli of clinical importance.

Nanopore Electrochemistry for Pathogen Detection

Pathogen infections have seriously threatened human health, and there is an urgent demand for rapid and efficient pathogen identification to provide instructions in clinical diagnosis and therapeutic intervention. Recently, nanopore technology, a rapidly maturing technology which delivers ultrasensitive sensing and high throughput in real-time and at low cost, has achieved success in pathogen detection. Furthermore, the remarkable development of nanopore sequencing, for example, the MinION sequencer from Oxford Nanopore Technologies (ONT) as a competitive sequencing technology, has facilitated the rapid analysis of disease-related microbiomes at the whole-genome level and on a large scale. Here, we highlighted recent advances in nanopore approaches for pathogen detection at the single-molecule level. We also overviewed the applications of nanopore sequencing in pathogenic bacteria identification and diagnosis. In the end, we discussed the challenges and future developments of nanopore technology as promising tools for the management of infections, which may be helpful to aid understanding as well as decision-making.

Combination of Whole Genome Sequencing and Metagenomics for Microbiological Diagnostics

Whole genome sequencing (WGS) provides the highest resolution for genome-based species identification and can provide insight into the antimicrobial resistance and virulence potential of a single microbiological isolate during the diagnostic process. In contrast, metagenomic sequencing allows the analysis of DNA segments from multiple microorganisms within a community, either using an amplicon- or shotgun-based approach. However, WGS and shotgun metagenomic data are rarely combined, although such an approach may generate additive or synergistic information, critical for, e.g., patient management, infection control, and pathogen surveillance. To produce a combined workflow with actionable outputs, we need to understand the pre-to-post analytical process of both technologies. This will require specific databases storing interlinked sequencing and metadata, and also involves customized bioinformatic analytical pipelines. This review article will provide an overview of the critical steps and potential clinical application of combining WGS and metagenomics together for microbiological diagnosis.

Genomic analysis reveals the role of integrative and conjugative elements in plant pathogenic bacteria

Background

ICEs are mobile genetic elements found integrated into bacterial chromosomes that can excise and be transferred to a new cell. They play an important role in horizontal gene transmission and carry accessory genes that may provide interesting phenotypes for the bacteria. Here, we seek to research the presence and the role of ICEs in 300 genomes of phytopathogenic bacteria with the greatest scientific and economic impact.

Results

Seventy-eight ICEs (45 distinct elements) were identified and characterized in chromosomes of Agrobacterium tumefaciens, Dickeya dadantii, and D. solani, Pectobacterium carotovorum and P. atrosepticum, Pseudomonas syringae, Ralstonia solanacearum Species Complex, and Xanthomonas campestris. Intriguingly, the co-occurrence of four ICEs was observed in some P. syringae strains. Moreover, we identified 31 novel elements, carrying 396 accessory genes with potential influence on virulence and fitness, such as genes coding for functions related to T3SS, cell wall degradation and resistance to heavy metals. We also present the analysis of previously reported data on the expression of cargo genes related to the virulence of P. atrosepticum ICEs, which evidences the role of these genes in the infection process of tobacco plants.

Conclusions

Altogether, this paper has highlighted the potential of ICEs to affect the pathogenicity and lifestyle of these phytopathogens and direct the spread of significant putative virulence genes in phytopathogenic bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00275-1.

Characterization of qnrB-carrying plasmids from ESBL- and non-ESBL-producing Escherichia coli

BACKGROUND

Escherichia coli carrying clinically important antimicrobial resistances [i.e., against extended-spectrum-beta-lactamases (ESBL)] are of high concern for human health and are increasingly detected worldwide. Worryingly, they are often identified as multidrug-resistant (MDR) isolates, frequently including resistances against quinolones/fluoroquinolones.

RESULTS

Here, the occurrence and genetic basis of the fluoroquinolone resistance enhancing determinant qnrB in ESBL-/non-ESBL-producing E. coli was investigated. Overall, 33 qnrB-carrying isolates out of the annual German antimicrobial resistance (AMR) monitoring on commensal E. coli (incl. ESBL-/AmpC-producing E. coli) recovered from food and livestock between 2013 and 2018 were analysed in detail. Whole-genome sequencing, bioinformatics analyses and transferability evaluation was conducted to characterise the prevailing qnrB-associated plasmids. Furthermore, predominant qnrB-carrying plasmid-types were subjected to in silico genome reconstruction analysis. In general, the qnrB-carrying E. coli were found to be highly heterogenic in their multilocus sequence types (STs) and their phenotypic resistance profiles. Most of them appeared to be MDR and exhibited resistances against up to ten antimicrobials of different classes. With respect to qnrB-carrying plasmids, we found qnrB19 located on small Col440I plasmids to be most widespread among ESBL-producing E. coli from German livestock and food. This Col440I plasmid-type was found to be highly conserved by exhibiting qnrB19, a pspF operon and different genes of unassigned function. Furthermore, we detected plasmids of the incompatibility groups IncN and IncH as carriers of qnrB. All qnrB-carrying plasmids also exhibited virulence factors and various insertion sequences (IS). The majority of the qnrB-carrying plasmids were determined to be self-transmissible, indicating their possible contribution to the spread of resistances against (fluoro)quinolones and other antimicrobials.

CONCLUSION

In this study, a diversity of different plasmid types carrying qnrB alone or in combination with other resistance determinants (i.e., beta-lactamase genes) were found. The spread of these plasmids, especially those carrying antimicrobial resistance genes against highest priority critically important antimicrobial agents, is highly unfavourable and can pose a threat for public health. Therefore, the dissemination pathways and evolution of these plasmids need to be further monitored.

Exploiting genomics to mitigate the public health impact of antimicrobial resistance

Antimicrobial resistance (AMR) is a major global public health threat, which has been largely driven by the excessive use of antimicrobials. Control measures are urgently needed to slow the trajectory of AMR but are hampered by an incomplete understanding of the interplay between pathogens, AMR encoding genes, and mobile genetic elements at a microbial level. These factors, combined with the human, animal, and environmental interactions that underlie AMR dissemination at a population level, make for a highly complex landscape. Whole-genome sequencing (WGS) and, more recently, metagenomic analyses have greatly enhanced our understanding of these processes, and these approaches are informing mitigation strategies for how we better understand and control AMR. This review explores how WGS techniques have advanced global, national, and local AMR surveillance, and how this improved understanding is being applied to inform solutions, such as novel diagnostic methods that allow antimicrobial use to be optimised and vaccination strategies for better controlling AMR. We highlight some future opportunities for AMR control informed by genomic sequencing, along with the remaining challenges that must be overcome to fully realise the potential of WGS approaches for international AMR control.

Ultra-deep long-read sequencing detects IS-mediated gene duplications as a potential trigger to generate arrays of resistance genes and a mechanism to induce novel gene variants such as blaCTX-M-243

BACKGROUND

Extended-spectrum β-lactamases (ESBLs) are enzymes that can render their hosts resistant to various β-lactam antibiotics. CTX-M-type enzymes are the most prevalent ESBLs and the main cause of resistance to third-generation cephalosporins in Enterobacteriaceae. The number of described CTX-M types is continuously rising, currently comprising over 240 variants. During routine screening we identified a novel blaCTX-M gene.

OBJECTIVES

To characterize a novel blaCTX-M variant harboured by a multidrug-resistant Escherichia coli isolate of sequence type ST354.

METHODS

Antibiotic susceptibilities were determined using broth microdilution. Genome and plasmid sequences were reconstructed using short- and long-read sequencing. The novel blaCTX-M locus was analysed using long-read and Sanger sequencing. Plasmid polymorphisms were determined in silico on a single plasmid molecule level.

RESULTS

The novel blaCTX-M-243 allele was discovered alongside a nearly identical blaCTX-M-104-containing gene array on a 219 kbp IncHI2A plasmid. CTX-M-243 differed from CTX-M-104 by only one amino acid substitution (N109S). Ultra-deep (2300-fold coverage) long-read sequencing revealed dynamic scaling of the blaCTX-M genetic contexts from one to five copies. Further antibiotic resistance genes such as blaTEM-1 also exhibited sequence heterogeneity but were stable in copy number.

CONCLUSIONS

We identified the novel ESBL gene blaCTX-M-243 and illustrate a dynamic system of varying blaCTX-M copy numbers. Our results highlight the constant emergence of new CTX-M family enzymes and demonstrate a potential evolutionary platform to generate novel ESBL variants and possibly other antibiotic resistance genes.

Application and Challenge of 3rd Generation Sequencing for Clinical Bacterial Studies

Over the past 25 years, the powerful combination of genome sequencing and bioinformatics analysis has played a crucial role in interpreting information encoded in bacterial genomes. High-throughput sequencing technologies have paved the way towards understanding an increasingly wide range of biological questions. This revolution has enabled advances in areas ranging from genome composition to how proteins interact with nucleic acids. This has created unprecedented opportunities through the integration of genomic data into clinics for the diagnosis of genetic traits associated with disease. Since then, these technologies have continued to evolve, and recently, long-read sequencing has overcome previous limitations in terms of accuracy, thus expanding its applications in genomics, transcriptomics and metagenomics. In this review, we describe a brief history of the bacterial genome sequencing revolution and its application in public health and molecular epidemiology. We present a chronology that encompasses the various technological developments: whole-genome shotgun sequencing, high-throughput sequencing, long-read sequencing. We mainly discuss the application of next-generation sequencing to decipher bacterial genomes. Secondly, we highlight how long-read sequencing technologies go beyond the limitations of traditional short-read sequencing. We intend to provide a description of the guiding principles of the 3rd generation sequencing applications and ongoing improvements in the field of microbial medical research.

Comparison of Reference-Based Assembly and De Novo Assembly for Bacterial Plasmid Reconstruction and AMR Gene Localization in Salmonella enterica Serovar Schwarzengrund Isolates

It is well established that plasmids carrying multiple antimicrobial resistance (AMR) genes can be easily transferred among bacterial isolates by horizontal gene transfer. Previous studies have shown that a combination of short- and long-read approaches is effective in reconstructing accurate plasmids. However, high-quality Illumina short reads mapped onto the long reads in the context of an AMR hybrid monitoring strategy have not yet been explored. Hence, this study aimed to improve the reconstruction of plasmids, including the localization of AMR genes, using the above-described parameters on whole-genome sequencing (WGS) results. To the best of our knowledge, this study is the first to use S1 nuclease pulsed-field gel electrophoresis (S1-PFGE) to confirm the number and sizes of plasmids detected by in silico-based predictions in Salmonella strains. Our results showed that de novo assembly did not detect the number of bacterial plasmids more accurately than reference-based assembly did. As this new hybrid mapping strategy surpassed de novo assembly in bacterial reconstruction, it was further used to identify the presence and genomic location of AMR genes among three Salmonella enterica serovar Schwarzengrund isolates. The AMR genes identified in the bacterial chromosome among the three Salmonella enterica serovar Schwarzengrund isolates included: AAC(3)-IV, AAC(6')-Iy, aadA2, APH(4)-Ia, cmlA1, golS, mdsA, mdsB, mdsC, mdtK, qacH, sdiA, sul2, sul3, and TEM-1 genes. Moreover, the presence of TEM-1, AAC(3)-IV, aadA2, APH(4)-Ia, cmlA1, dfrA12, floR, sul1, sul3, and tet(A) genes found within three IncFIB plasmids and one IncX1 plasmid highlight their possible transmission into the environment, which is a public health risk. In conclusion, the generated data using this new hybrid mapping strategy will contribute to the improvement of AMR monitoring and support the risk assessment of AMR dissemination.

Mobile Genetic Elements Drive Antimicrobial Resistance Gene Spread in Pasteurellaceae Species

Mobile genetic elements (MGEs) and antimicrobial resistance (AMR) drive important ecological relationships in microbial communities and pathogen-host interaction. In this study, we investigated the resistome-associated mobilome in 345 publicly available Pasteurellaceae genomes, a large family of Gram-negative bacteria including major human and animal pathogens. We generated a comprehensive dataset of the mobilome integrated into genomes, including 10,820 insertion sequences, 2,939 prophages, and 43 integrative and conjugative elements. Also, we assessed plasmid sequences of Pasteurellaceae. Our findings greatly expand the diversity of MGEs for the family, including a description of novel elements. We discovered that MGEs are comparable and dispersed across species and that they also co-occur in genomes, contributing to the family's ecology via gene transfer. In addition, we investigated the impact of these elements in the dissemination and shaping of AMR genes. A total of 55 different AMR genes were mapped to 721 locations in the dataset. MGEs are linked with 77.6% of AMR genes discovered, indicating their important involvement in the acquisition and transmission of such genes. This study provides an uncharted view of the Pasteurellaceae by demonstrating the global distribution of resistance genes linked with MGEs.

Hybrid Assembly Provides Improved Resolution of Plasmids, Antimicrobial Resistance Genes, and Virulence Factors in Escherichia coli and Klebsiella pneumoniae Clinical Isolates

Emerging new sequencing technologies have provided researchers with a unique opportunity to study factors related to microbial pathogenicity, such as antimicrobial resistance (AMR) genes and virulence factors. However, the use of whole-genome sequence (WGS) data requires good knowledge of the bioinformatics involved, as well as the necessary techniques. In this study, a total of nine Escherichia coli and Klebsiella pneumoniae isolates from Norwegian clinical samples were sequenced using both MinION and Illumina platforms. Three out of nine samples were sequenced directly from blood culture, and one sample was sequenced from a mixed-blood culture. For genome assembly, several long-read, (Canu, Flye, Unicycler, and Miniasm), short-read (ABySS, Unicycler and SPAdes) and hybrid assemblers (Unicycler, hybridSPAdes, and MaSurCa) were tested. Assembled genomes from the best-performing assemblers (according to quality checks using QUAST and BUSCO) were subjected to downstream analyses. Flye and Unicycler assemblers performed best for the assembly of long and short reads, respectively. For hybrid assembly, Unicycler was the top-performing assembler and produced more circularized and complete genome assemblies. Hybrid assembled genomes performed substantially better in downstream analyses to predict putative plasmids, AMR genes and β-lactamase gene variants, compared to MinION and Illumina assemblies. Thus, hybrid assembly has the potential to reveal factors related to microbial pathogenicity in clinical and mixed samples.

Complete genome for Actinobacillus pleuropneumoniae serovar 8 reference strain 405: comparative analysis with draft genomes for different laboratory stock cultures indicates little genetic variation

We report here the complete genome sequence of the widely studied Actinobacillus pleuropneumoniae serovar 8 reference strain 405, generated using the Pacific Biosciences (PacBio) RS II platform. Furthermore, we compared draft sequences generated by Illumina sequencing of six stocks of this strain, including the same original stock used to generate the PacBio sequence, held in different countries and found little genetic variation, with only three SNPs identified, all within the degS gene. However, sequences of two small plasmids, pARD3079 and p405tetH, detected by Illumina sequencing of the draft genomes were not identified in the PacBio sequence of the reference strain.

Interplay between Phenotypic Resistance to Relevant Antibiotics in Gram-Negative Urinary Pathogens: A Data-Driven Analysis of 10 Years' Worth of Antibiogram Data

The global emergence of antimicrobial resistance (AMR) has become a critical issue for clinicians, as it puts the decades of developments in the medical field in jeopardy, by severely limiting the useful therapeutic arsenal of drugs, both in nosocomial and community-acquired infections. In the present study, a secondary analysis of taxonomic and resistance data was performed, corresponding to urinary tract infections (UTIs) caused by Gram-negative bacteria, detected between 1 January 2008 to 31 December 2017 at the Albert Szent-Györgyi Health Center, University of Szeged. The following were identifiable from the data collected: year of isolation; outpatient (OP)/inpatient (IP) origin of the isolate; taxonomy; and susceptibility/resistance to selected indicator antibiotics. Principal component analysis (PCA) and a correlation matrix were used to determine the association between the presences of resistance against indicator antibiotics in each taxonomic group. Overall, data from n = 16,240 outpatient and n = 13,964 inpatient Gram-negative UTI isolates were included in the data analyses. In E. coli, strong positive correlations were seen between resistance to ciprofloxacin (CIP) and gentamicin (GEN) resistance (OP: r = 0.6342, p = 0.049; IP: r = 0.9602, p < 0.001), whereas strong negative correlations were shown for fosfomycin (FOS) and nitrofurantoin (NIT) resistance (OP: r = -0.7183, p = 0.019; IP: r = -0.7437; p = 0.014). For Klebsiella spp. isolates, CIP resistance showed strong positive correlation with resistance to third-generation cephalosporins (3GC) and GEN (r = 0.7976, p = 0.006 and r = 0.7428, p = 0.014, respectively) in OP isolates, and with resistance to trimethoprim-sulfamethoxazole (SXT) and FOS (r = 0.8144, p = 0.004 and r = 0.7758, p < 0.001, respectively) in IP isolates. For members of the Citrobacter-Enterobacter-Serratia group, the resistance among indicator antibiotics showed a strong positive correlation, with the exception of FOS resistance. In the Proteus-Providencia-Morganella group, the strongest association was noted between CIP and SXT resistance (OP: r = 0.9251, p < 0.001; IP: r = 0.8007; p = 0.005). In the case of OP Acinetobacter spp., CIP showed strong and significant positive correlations with most indicator antibiotics, whereas for IP isolates, strong negative correlations arose among imipenem (IMI) resistance and resistance to other drugs. For Pseudomonas spp., strong and positive correlations were noted among resistance to β-lactam antibiotics and aminoglycosides, with the exception of ceftazidime (CEFT), showing strong, but negative correlations. Though molecular tests and sequencing-based platforms are now considered as the gold-standard for AMR surveillance, standardized collection of phenotypic resistance data and the introduction of Big Data analytic methods may be a viable alternative for molecular surveillance, especially in low-resource settings.

Antibiotic-resistant bacteria, antibiotic resistance genes, and antibiotic residues in wastewater from a poultry slaughterhouse after conventional and advanced treatments

Slaughterhouse wastewater is considered a reservoir for antibiotic-resistant bacteria and antibiotic residues, which are not sufficiently removed by conventional treatment processes. This study focuses on the occurrence of ESKAPE bacteria (Enterococcus spp., S. aureus, K. pneumoniae, A. baumannii, P. aeruginosa, Enterobacter spp.), ESBL (extended-spectrum β-lactamase)-producing E. coli, antibiotic resistance genes (ARGs) and antibiotic residues in wastewater from a poultry slaughterhouse. The efficacy of conventional and advanced treatments (i.e., ozonation) of the in-house wastewater treatment plant regarding their removal was also evaluated. Target culturable bacteria were detected only in the influent and effluent after conventional treatment. High abundances of genes (e.g., bla_TEM, bla_CTX-M-15, bla_CTX-M-32, bla_OXA-48, bla_CMY and mcr-1) of up to 1.48 × 10⁶ copies/100 mL were detected in raw influent. All of them were already significantly reduced by 1–4.2 log units after conventional treatment. Following ozonation, mcr-1 and bla_CTX-M-32 were further reduced below the limit of detection. Antibiotic residues were detected in 55.6% (n = 10/18) of the wastewater samples. Despite the significant reduction through conventional and advanced treatments, effluents still exhibited high concentrations of some ARGs (e.g., sul1, ermB and bla_OXA-48), ranging from 1.75 × 10² to 3.44 × 10³ copies/100 mL. Thus, a combination of oxidative, adsorptive and membrane-based technologies should be considered.

Phenotypic and Genotypic Properties of Fluoroquinolone-Resistant, qnr-Carrying Escherichia coli Isolated from the German Food Chain in 2017

Fluoroquinolones are the highest priority, critically important antimicrobial agents. Resistance development can occur via different mechanisms, with plasmid-mediated quinolone resistance (PMQR) being prevalent in the livestock and food area. Especially, qnr genes, commonly located on mobile genetic elements, are major drivers for the spread of resistance determinants against fluoroquinolones. We investigated the prevalence and characteristics of qnr-positive Escherichia (E.) coli obtained from different monitoring programs in Germany in 2017. Furthermore, we aimed to evaluate commonalities of qnr-carrying plasmids in E. coli. We found qnr to be broadly spread over different livestock and food matrices, and to be present in various sequence types. The qnr-positive isolates were predominantly detected within selectively isolated ESBL (extended spectrum beta-lactamase)-producing E. coli, leading to a frequent association with other resistance genes, especially cephalosporin determinants. Furthermore, we found that qnr correlates with the presence of genes involved in resistance development against quaternary ammonium compounds (qac). The detection of additional point mutations in many isolates within the chromosomal QRDR region led to even higher MIC values against fluoroquinolones for the investigated E. coli. All of these attributes should be carefully taken into account in the risk assessment of qnr-carrying E. coli from livestock and food.