Polypolish: Short-read polishing of long-read bacterial genome assemblies

Chromosomal integration and plasmid fusion occurring in ST20 carbapenem-resistant Klebsiella pneumoniae isolates coharboring blaNDM-1 and blaIMP-4 induce resistance transmission and fitness variation

To investigate the epidemiology of ST20 carbapenem-resistant Klebsiella pneumoniae (CRKP) in China, and further explore the genomic characteristics of blaIMP-4 and blaNDM-1 coharboring isolates and plasmid contributions to resistance and fitness. Seven ST20 CRKP isolates were collected nationwide, and antimicrobial susceptibility testing was performed. Antimicrobial resistance genes, virulence genes, and plasmid replicons were identified via whole-genome sequencing, and clonality assessed via core-genome multilocus sequence typing. Furthermore, we found four dual-metallo-β-lactamases (MBL)-harbouring isolates, the gene location was detected by Southern blotting, and plasmid location analysis showed that blaIMP-4 was located on a separate plasmid, a self-conjugative fusion plasmid, or the bacterial chromosome. These isolates were subjected to long-read sequencing, the presence of blaIMP-4 in different locations was identified by genomic comparison, and transposon units were detected via inverse PCR. We subsequently found that blaIMP-4 on the fusion plasmid and bacterial chromosome was formed via intact plasmid recombination by the IS26 and ltrA, respectively, and the circular transposon unit was related to cointegration, however, blaIMP-4 in different locations did not affect the gene stability. The blaNDM-1-harbouring plasmid contributed to the increased resistance to β-lactams and shortened survival lag time which was revealed in plasmid cured isolates. In summary, the K. pneumoniae ST20 clone is a high-risk resistant clone. With the use of ceftazidime/avibactam, MBL-positive isolates, especially dual-MBL-harbouring isolates, should be given additional attention.

OXA-48-like carbapenemases in Proteus mirabilis - novel genetic environments and a challenge for detection

OXA-48-like enzymes represent the most frequently detected carbapenemases in Enterobacterales in Western Europe, North Africa and the Middle East. In contrast to other species, the presence of OXA-48-like in Proteus mirabilis leads to an unusually susceptible phenotype with low MICs for carbapenems and piperacillin-tazobactam, which is easily missed in the diagnostic laboratory. So far, there is little data available on the genetic environments of the corresponding genes, blaOXA-48-like, in P. mirabilis. In this study susceptibility phenotypes and genomic data of 13 OXA-48-like-producing P. mirabilis were investigated (OXA-48, n=9; OXA-181, n=3; OXA-162, n=1). Ten isolates were susceptible to meropenem and ertapenem and three isolates were susceptible to piperacillin-tazobactam. The gene blaOXA-48 was chromosomally located in 7/9 isolates. Thereof, in three isolates blaOXA-48 was inserted into a P. mirabilis genomic island. Of the three isolates harbouring blaOXA-181 one was located on an IncX3 plasmid and two were located on a novel MOBF plasmid, pOXA-P12, within the new transposon Tn7713. In 5/6 isolates with plasmidic location of blaOXA-48-like, the plasmids could conjugate to E. coli recipients in vitro. Vice versa, blaOXA-48-carrying plasmids could conjugate from other Enterobacterales into a P. mirabilis recipient. These data show a high diversity of blaOXA-48-like genetic environments compared to other Enterobacterales, where genetic environments are quite homogenous. Given the difficult-to-detect phenotype of OXA-48-like-producing P. mirabilis and the location of blaOXA-48-like on mobile genetic elements, it is likely that OXA-48-like-producing P. mirabilis can disseminate, escape most surveillance systems, and contribute to a hidden spread of OXA-48-like.

Longitudinal analysis within one hospital in sub-Saharan Africa over 20 years reveals repeated replacements of dominant clones of Klebsiella pneumoniae and stresses the importance to include temporal patterns for vaccine design considerations

BACKGROUND

Infections caused by multidrug-resistant gram-negative bacteria present a severe threat to global public health. The WHO defines drug-resistant Klebsiella pneumoniae as a priority pathogen for which alternative treatments are needed given the limited treatment options and the rapid acquisition of novel resistance mechanisms by this species. Longitudinal descriptions of genomic epidemiology of Klebsiella pneumoniae can inform management strategies but data from sub-Saharan Africa are lacking.

METHODS

We present a longitudinal analysis of all invasive K. pneumoniae isolates from a single hospital in Blantyre, Malawi, southern Africa, from 1998 to 2020, combining clinical data with genome sequence analysis of the isolates.

RESULTS

We show that after a dramatic increase in the number of infections from 2016 K. pneumoniae becomes hyperendemic, driven by an increase in neonatal infections. Genomic data show repeated waves of clonal expansion of different, often ward-restricted, lineages, suggestive of hospital-associated transmission. We describe temporal trends in resistance and surface antigens, of relevance for vaccine development.

CONCLUSIONS

Our data highlight a clear need for new interventions to prevent rather than treat K. pneumoniae infections in our setting. Whilst one option may be a vaccine, the majority of cases could be avoided by an increased focus on and investment in infection prevention and control measures, which would reduce all healthcare-associated infections and not just one.

Rationally designed chromosome fusion does not prevent rapid growth of Vibrio natriegens

DNA replication is essential for the proliferation of all cells. Bacterial chromosomes are replicated bidirectionally from a single origin of replication, with replication proceeding at about 1000 bp per second. For the model organism, Escherichia coli, this translates into a replication time of about 40 min for its 4.6 Mb chromosome. Nevertheless, E. coli can propagate by overlapping replication cycles with a maximum short doubling time of 20 min. The fastest growing bacterium known, Vibrio natriegens, is able to replicate with a generation time of less than 10 min. It has a bipartite genome with chromosome sizes of 3.2 and 1.9 Mb. Is simultaneous replication from two origins a prerequisite for its rapid growth? We fused the two chromosomes of V. natriegens to create a strain carrying one chromosome with a single origin of replication. Compared to the parental, this strain showed no significant deviation in growth rate. This suggests that the split genome is not a prerequisite for rapid growth.

Vancomycin heteroresistance caused by unstable tandem amplifications of the vanM gene cluster on linear conjugative plasmids in a clinical Enterococcus faecium

Vancomycin heteroresistance is prone to missed detection and poses a risk of clinical treatment failure. We encountered one clinical Enterococcus faecium strain, SRR12, that carried a complete vanM gene cluster but was determined as susceptible to vancomycin using the broth microdilution method. However, distinct subcolonies appeared within the clear zone of inhibition in the E-test assay, one of which, named SRR12-v1, showed high-level resistance to vancomycin. SRR12 was confirmed as heteroresistant to vancomycin using population analysis profiling and displayed "revive" growth curves with a lengthy lag phase of over 13 hours when exposed to 2-32 mg/L vancomycin. The resistant subcolony SRR12-v1 was found to carry an identical vanM gene cluster to that of SRR12 but a significantly increased vanM copy number in the genome. Long-read whole genome sequencing revealed that a one-copy vanM gene cluster was located on a pELF1-like linear plasmid in SRR12. In comparison, tandem amplification of the vanM gene cluster jointed with IS1216E was seated on a linear plasmid in the genome of SRR12-v1. These amplifications of the vanM gene cluster were demonstrated as unstable and would decrease accompanied by fitness reversion after serial passaging for 50 generations under increasing vancomycin pressure or without antibiotic pressure but were relatively stable under constant vancomycin pressure. Further, vanM resistance in resistant variants was verified to be carried by conjugative plasmids with variable sizes using conjugation assays and S1-pulsed field gel electrophoresis blotting, suggesting the instability/flexibility of vanM cluster amplification in the genome and an increased risk of vanM resistance dissemination.

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.

Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing

Whole-genome reconstruction of bacterial pathogens has become an important tool for tracking transmission and antimicrobial resistance gene spread, but highly accurate and complete assemblies have largely only historically been achievable using hybrid long- and short-read sequencing. We previously found the Oxford Nanopore Technologies (ONT) R10.4/kit12 flowcell/chemistry produced improved assemblies over the R9.4.1/kit10 combination, however long-read only assemblies contained more errors compared to Illumina-ONT hybrid assemblies. ONT have since released an R10.4.1/kit14 flowcell/chemistry upgrade and recommended the use of Bovine Serum Albumin (BSA) during library preparation, both of which reportedly increase accuracy and yield. They have also released updated basecallers trained using native bacterial DNA containing methylation sites intended to fix systematic basecalling errors, including common adenosine (A) to guanine (G) and cytosine (C) to thymine (T) substitutions. To evaluate these improvements, we successfully sequenced four bacterial reference strains, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, and nine genetically diverse E. coli bloodstream infection-associated isolates from different phylogroups and sequence types, both with and without BSA. These sequences were de novo assembled and compared against Illumina-corrected reference genomes. In this small evaluation of 13 isolates we found that nanopore long-read-only R10.4.1/kit 14 assemblies with updated basecallers trained using bacterial methylated DNA produce accurate assemblies with ≥40×depth, sufficient to be cost-effective compared with hybrid ONT/Illumina sequencing in our setting.

Bacterial genome sequences of uncharacterized Chitinophaga species isolated from the International Space Station

We report four Chitinophaga sp. strains isolated from wastewater collected onboard the International Space Station. Here, we present three finished and one draft genome. Taxonomic ranks established by genome-based analysis indicate that these Chitinophaga sp. strains represent candidates for a new species.

Hybracter: Enabling Scalable, Automated, Complete and Accurate Bacterial Genome Assemblies

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants (SNVs). They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance (AMR) genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic, and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.

Complete genome sequences of Blautia hydrogenotrophica DSM 10507T isolated from human feces and Blautia coccoides DSM 935T isolated from mouse feces

We report on the closed genome sequences of the acetogen Blautia hydrogenotrophica S5a33T (DSM 10507T) and of Blautia coccoides CLC-1T (DSM 935T). The B. hydrogenotrophica S5a33T genome harbors a chromosome (3,590,609 bp) and a plasmid (7,176 bp). The B. coccoides CLC-1T genome consists of a single chromosome (6,097,890 bp).

Complete genome sequence of Escherichia coli MP1

Here, we report the complete genome sequence of Escherichia coli strain MP1, consisting of one circular chromosome and one circular plasmid. Long-read assembly was performed using a consensus approach, followed by long- and short-read polishing, and gene annotation.

Complete genome sequence of invasive Streptococcus pneumoniae serotype 11A causing meningitis in an adult patient, Italy 2022

Streptococcus pneumoniae is a major global health concern, being a common cause of meningitis in both children and adults. Here, we report the complete genome sequence of P10_PNE_LCR, a S. pneumoniae 11A strain isolated in Northern Italy from an adult patient diagnosed with meningitis.

Isolation of a PRD1-like phage uncovers the carriage of three putative conjugative plasmids in clinical Burkholderia contaminans

The Burkholderia cepacia complex (Bcc) is a group of increasingly multi-drug resistant opportunistic bacteria. This resistance is driven through a combination of intrinsic factors and the carriage of a broad range of conjugative plasmids harbouring virulence determinants. Therefore, novel treatments are required to treat and prevent further spread of these virulence determinants. In the search for phages infective for clinical Bcc isolates, CSP1 phage, a PRD1-like phage was isolated. CSP1 phage was found to require pilus machinery commonly encoded on conjugative plasmids to facilitate infection of Gram-negative bacteria genera including Escherichia and Pseudomonas. Whole genome sequencing and characterisation of one of the clinical Burkholderia isolates revealed it to be Burkholderia contaminans. B. contaminans 5080 was found to contain a genome of over 8 Mbp encoding multiple intrinsic resistance factors, such as efflux pump systems, but more interestingly, carried three novel plasmids encoding multiple putative virulence factors for increased host fitness, including antimicrobial resistance. Even though PRD1-like phages are broad host range, their use in novel antimicrobial treatments shouldn't be dismissed, as the dissemination potential of conjugative plasmids is extensive. Continued survey of clinical bacterial strains is also key to understanding the spread of antimicrobial resistance determinants and plasmid evolution.

Pathogenomic analysis and characterization of Pasteurella multocida strains recovered from human infections

Pasteurella multocida is an upper respiratory tract commensal in several mammal and bird species but can also cause severe disease in humans and in production animals such as poultry, cattle, and pigs. In this study, we performed whole-genome sequencing of P. multocida isolates recovered from a range of human infections, from the mouths of cats, and from wounds on dogs. Together with publicly available P. multocida genome sequences, we performed phylogenetic and comparative genomic analyses. While isolates from cats and dogs were spread across the phylogenetic tree, human infections were caused almost exclusively by subsp. septica strains. Most of the human isolates were capsule type A and LPS type L1 and L3; however, some strains lacked a capsule biosynthesis locus, and some strains contained a novel LPS outer-core locus, distinct from the eight LPS loci that can currently be identified using an LPS multiplex PCR. In addition, the P. multocida strains isolated from human infections contained novel mobile genetic elements. We compiled a curated database of known P. multocida virulence factor and antibiotic resistance genes (PastyVRDB) allowing for detailed characterization of isolates. The majority of human P. multocida isolates encoded a reduced range of iron receptors and contained only one filamentous hemagglutinin gene. Finally, gene-trait analysis identified a putative L-fucose uptake and utilization pathway that was over-represented in subsp. septica strains and may represent a novel host predilection mechanism in this subspecies. Together, these analyses have identified pathogenic mechanisms likely important for P. multocida zoonotic infections.IMPORTANCEPasteurella multocida can cause serious infections in humans, including skin and wound infections, pneumonia, peritonitis, meningitis, and bacteraemia. Cats and dogs are known vectors of human pasteurellosis, transmitting P. multocida via bite wounds or contact with animal saliva. The mechanisms that underpin P. multocida human predilection and pathogenesis are poorly understood. With increasing identification of antibiotic-resistant P. multocida strains, understanding these mechanisms is vital for developing novel treatments and control strategies to combat P. multocida human infection. Here, we show that a narrow range of P. multocida strains cause disease in humans, while cats and dogs, common vectors for zoonotic infections, can harbor a wide range of P. multocida strains. We also present a curated P. multocida-specific database, allowing quick and detailed characterization of newly sequenced P. multocida isolates.

Identifying functional multi-host shuttle plasmids to advance synthetic biology applications in Mesorhizobium and Bradyrhizobium

Ammonia availability has a crucial role in agriculture as it ensures healthy plant growth and increased crop yields. Since diazotrophs are the only organisms capable of reducing dinitrogen to ammonia, they have great ecological importance and potential to mitigate the environmental and economic costs of synthetic fertilizer use. Rhizobia are especially valuable being that they can engage in nitrogen-fixing symbiotic relationships with legumes, and they demonstrate great diversity and plasticity in genomic and phenotypic traits. However, few rhizobial species have sufficient genetic tractability for synthetic biology applications. This study established a basic genetic toolbox with antibiotic resistance markers, multi-host shuttle plasmids and a streamlined protocol for biparental conjugation with Mesorhizobium and Bradyrhizobium species. We identified two repABC origins of replication from Sinorhizobium meliloti (pSymB) and Rhizobium etli (p42d) that were stable across all three strains of interest. Furthermore, the NZP2235 genome was sequenced and phylogenetic analysis determined its reclassification to Mesorhizobium huakuii. These tools will enable the use of plasmid-based strategies for more advanced genetic engineering projects and ultimately contribute towards the development of more sustainable agriculture practices by means of novel nitrogen-fixing organelles, elite bioinoculants, or symbiotic association with nonlegumes.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

Snakemake workflows for long-read bacterial genome assembly and evaluation

With the advancement of long-read sequencing technologies and their increasing use for bacterial genomics, several methods for generating genome assemblies from error-prone long reads have been developed. These are complemented by various tools for assembly polishing using either long reads, short reads, or reference genomes. End users are therefore left with a plethora of possible combinations of programs for obtaining a final trusted assembly. Hence, there is also a need to measure the completeness and accuracy of such assemblies, for which, again, several evaluation methods implemented in various programs are available. In order to automatically run multiple genome assembly and evaluation programs at once, I developed two workflows for the workflow management system Snakemake, which provide end users with an easy-to-run solution for testing various genome assemblies from their sequencing data. Both workflows use the conda packaging system, so there is no need for manual installation of each program.

Availability & Implementation

The workflows are available as open source software under the MIT license at github.com/pmenzel/ont-assembly-snake and github.com/pmenzel/score-assemblies.

Klebsiella pneumoniae exhibiting a phenotypic hyper-splitting phenomenon including the formation of small colony variants

In this study, we characterized a Klebsiella pneumoniae strain in a patient with shrapnel hip injury, which resulted in multiple phenotypic changes, including the formation of a small colony variant (SCV) phenotype. Although already described since the 1960s, there is little knowledge about SCV phenotypes in Enterobacteriaceae. The formation of SCVs has been recognized as a bacterial strategy to evade host immune responses and compromise the efficacy of antimicrobial therapies, leading to persistent and recurrent courses of infections. In this case, 14 isolates with different resisto- and morpho-types were distinguished from the patient's urine and tissue samples. Whole genome sequencing revealed that all isolates were clonally identical belonging to the K. pneumoniae high-risk sequence type 147. Subculturing the SCV colonies consistently resulted in the reappearance of the initial SCV phenotype and three stable normal-sized phenotypes with distinct morphological characteristics. Additionally, an increase in resistance was observed over time in isolates that shared the same colony appearance. Our findings highlight the complexity of bacterial behavior by revealing a case of phenotypic "hyper-splitting" in a K. pneumoniae SCV and its potential clinical significance.

Comparative genomics of an extensively drug resistant strain Klebsiella pneumoniae IITR008 with international high-risk clonal lineage ST147 isolated from river water

Carbapenem resistant Klebsiella pneumoniae causing severe infection resulting in morbidity and mortality have become a global health concern. K. pneumoniae with sequence type ST147 is an international high-risk clonal lineage, genomic studies have been done on K. pneumoniae ST147 isolated from clinical origin but genomic data for environmental K. pneumoniae ST147 is very scarce. Herein, K. pneumoniae IITR008, an extensively drug resistant and potentially hypervirulent bacterium, was isolated from Triveni Sangam, the confluence of three rivers where religious congregations are organized. Phenotypic, genomic and comparative genomic analysis of strain IITR008 was performed. Antibiotic susceptibility profiling revealed resistance to 9 different classes of antibiotics including ß-lactams, ß-lactam combination agents, carbapenem, aminoglycoside, macrolide, quinolones, cephams, phenicol, and folate pathway antagonists and was found to be susceptible to only tetracycline. The strain IITR008 possesses hypervirulence genes namely, iutA and iroN in addition to numerous virulence factors coding for adherence, regulation, iron uptake, secretion system and toxin. Both the IITR008 chromosome and plasmid pIITR008_75 possess a plethora of clinically relevant antibiotic-resistant genes (ARGs) including blaCTX-M-15, blaTEM-1, and blaSHV-11, corroborating the phenotypic resistance. Comparative genomic analysis with other ST147 K. pneumoniae provided insights on the phylogenetic clustering of IITR008 with a clinical strain isolated from a patient in Czech with recent travel history in India and other clinical strains isolated from India and Pakistan. According to the 'One Health' perspective, surveillance of antibiotic resistance in the environment is crucial to impede its accelerated development in diverse ecological niches.

Complete genome sequence of new Microbacterium sp. strain ET2, isolated from roots of leafless orchid

Whole-genome sequence of ET2 strain, isolated from the roots of leafless orchid, constitutes a single circular chromosome of 3,604,840 bp (69.44% G + C content). BLAST+-based average nucleotide identity (ANIb) and digital DNA-DNA hybridization values indicate that ET2 may be a novel Microbacterium species. Genes putatively involved in plant-microbial interactions were predicted.

Complete genome sequences of a Canadian strain of enteroaggregative Escherichia coli (EAEC) with multiple metals and antimicrobial resistance genes isolated from municipal waste-activated sludge

Enteroaggregative Escherichia coli (EAEC) is an emerging food-borne pathogen causing acute or persistent diarrhea in humans. Here, we report the complete genome sequence of a strain of EAEC with multiple metals and antimicrobial resistance genes isolated from a waste-activated sludge collected from a Canadian municipal wastewater treatment plant.

Carbapenem-resistant hypervirulent ST23 Klebsiella pneumoniae with a highly transmissible dual-carbapenemase plasmid in Chile

BACKGROUND

The convergence of hypervirulence and carbapenem resistance in the bacterial pathogen Klebsiella pneumoniae represents a critical global health concern. Hypervirulent K. pneumoniae (hvKp) strains, frequently from sequence type 23 (ST23) and having a K1 capsule, have been associated with severe community-acquired invasive infections. Although hvKp were initially restricted to Southeast Asia and primarily antibiotic-sensitive, carbapenem-resistant hvKp infections are reported worldwide. Here, within the carbapenemase production Enterobacterales surveillance system headed by the Chilean Public Health Institute, we describe the isolation in Chile of a high-risk ST23 dual-carbapenemase-producing hvKp strain, which carbapenemase genes are encoded in a single conjugative plasmid.

RESULTS

Phenotypic and molecular tests of this strain revealed an extensive resistance to at least 15 antibiotic classes and the production of KPC-2 and VIM-1 carbapenemases. Unexpectedly, this isolate lacked hypermucoviscosity, challenging this commonly used hvKp identification criteria. Complete genome sequencing and analysis confirmed the K1 capsular type, the KpVP-1 virulence plasmid, and the GIE492 and ICEKp10 genomic islands carrying virulence factors strongly associated with hvKp. Although this isolate belonged to the globally disseminated hvKp clonal group CG23-I, it is unique, as it formed a clade apart from a previously reported Chilean ST23 hvKp isolate and acquired an IncN KPC-2 plasmid highly disseminated in South America (absent in other hvKp genomes), but now including a class-I integron carrying blaVIM-1 and other resistance genes. Notably, this isolate was able to conjugate the double carbapenemase plasmid to an E. coli recipient, conferring resistance to 1st -5th generation cephalosporins (including combinations with beta-lactamase inhibitors), penicillins, monobactams, and carbapenems.

CONCLUSIONS

We reported the isolation in Chile of high-risk carbapenem-resistant hvKp carrying a highly transmissible conjugative plasmid encoding KPC-2 and VIM-1 carbapenemases, conferring resistance to most beta-lactams. Furthermore, the lack of hypermucoviscosity argues against this trait as a reliable hvKp marker. These findings highlight the rapid evolution towards multi-drug resistance of hvKp in Chile and globally, as well as the importance of conjugative plasmids and other mobile genetic elements in this convergence. In this regard, genomic approaches provide valuable support to monitor and obtain essential information on these priority pathogens and mobile elements.

A synthetic biology approach to assemble and reboot clinically relevant Pseudomonas aeruginosa tailed phages

The rise in the frequency of antibiotic resistance has made bacterial infections, specifically Pseudomonas aeruginosa, a cause for greater concern. Phage therapy is a promising solution that uses naturally isolated phages to treat bacterial infections. Ecological limitations, which stipulate a discrete host range and the inevitable evolution of resistance, may be overcome through a better understanding of phage biology and the utilization of engineered phages. In this study, we developed a synthetic biology approach to construct tailed phages that naturally target clinically relevant strains of Pseudomonas aeruginosa. As proof of concept, we successfully cloned and assembled the JG024 and DMS3 phage genomes in yeast using transformation-associated recombination cloning and rebooted these two phage genomes in two different strains of P. aeruginosa. We identified factors that affected phage reboot efficiency like the phage species or the presence of antiviral defense systems in the bacterial strain. We have successfully extended this method to two other phage species and observed that the method enables the reboot of phages that are naturally unable to infect the strain used for reboot. This research represents a critical step toward the construction of clinically relevant, engineered P. aeruginosa phages.IMPORTANCEPseudomonas aeruginosa is a bacterium responsible for severe infections and a common major complication in cystic fibrosis. The use of antibiotics to treat bacterial infections has become increasingly difficult as antibiotic resistance has become more prevalent. Phage therapy is an alternative solution that is already being used in some European countries, but its use is limited by the narrow host range due to the phage receptor specificity, the presence of antiviral defense systems in the bacterial strain, and the possible emergence of phage resistance. In this study, we demonstrate the use of a synthetic biology approach to construct and reboot clinically relevant P. aeruginosa tailed phages. This method enables a significant expansion of possibilities through the construction of engineered phages for therapy applications.

Closing the gap: Oxford Nanopore Technologies R10 sequencing allows comparable results to Illumina sequencing for SNP-based outbreak investigation of bacterial pathogens

Whole-genome sequencing has become the method of choice for bacterial outbreak investigation, with most clinical and public health laboratories currently routinely using short-read Illumina sequencing. Recently, long-read Oxford Nanopore Technologies (ONT) sequencing has gained prominence and may offer advantages over short-read sequencing, particularly with the recent introduction of the R10 chemistry, which promises much lower error rates than the R9 chemistry. However, limited information is available on its performance for bacterial single-nucleotide polymorphism (SNP)-based outbreak investigation. We present an open-source workflow, Prokaryotic Awesome variant Calling Utility (PACU) (https://github.com/BioinformaticsPlatformWIV-ISP/PACU), for constructing SNP phylogenies using Illumina and/or ONT R9/R10 sequencing data. The workflow was evaluated using outbreak data sets of Shiga toxin-producing Escherichia coli and Listeria monocytogenes by comparing ONT R9 and R10 with Illumina data. The performance of each sequencing technology was evaluated not only separately but also by integrating samples sequenced by different technologies/chemistries into the same phylogenomic analysis. Additionally, the minimum sequencing time required to obtain accurate phylogenetic results using nanopore sequencing was evaluated. PACU allowed accurate identification of outbreak clusters for both species using all technologies/chemistries, but ONT R9 results deviated slightly more from the Illumina results. ONT R10 results showed trends very similar to Illumina, and we found that integrating data sets sequenced by either Illumina or ONT R10 for different isolates into the same analysis produced stable and highly accurate phylogenomic results. The resulting phylogenies for these two outbreaks stabilized after ~20 hours of sequencing for ONT R9 and ~8 hours for ONT R10. This study provides a proof of concept for using ONT R10, either in isolation or in combination with Illumina, for rapid and accurate bacterial SNP-based outbreak investigation.

Draft genome of Brasenia schreberi, a worldwide distributed and endangered aquatic plant

OBJECTIVES

Brasenia is a monotypic genus in the family of Cabombaceae. The only species, B. schreberi, is a macrophyte distributed worldwide. Because it requires good water quality, it is endangered in China and other countries due to the deterioration of aquatic habitats. The young leaves and stems of B. schreberi are covered by thick mucilage, which has high medical value. As an allelopathic aquatic plant, it can also be used in the management of aquatic weeds. Here, we present its assembled and annotated genome to help shed light on medial and allelopathic substrates and facilitate their conservation.

DATA DESCRIPTION

Genomic DNA and RNA extracted from B. schreberi leaf tissues were used for whole genome and RNA sequencing using a Nanopore and/or MGI sequencer. The assembly was 1,055,148,839 bp in length, with 92 contigs and an N50 of 22,379,495 bp. The repetitive elements in the assembly were 555,442,205 bp. A completeness assessment of the assembly with BUSCO and compleasm indicated 88.4 and 90.9% completeness in the Eudicots database and 95.4 and 96.6% completeness in the Embryphyta database. Gene annotation revealed 67,747 genes that coded for 73,344 proteins.

Unveiling genome plasticity and a novel phage in Mycoplasma felis: Genomic investigations of four feline isolates

Mycoplasma felis has been isolated from diseased cats and horses, but to date only a single fully assembled genome of this species, of an isolate from a horse, has been characterized. This study aimed to characterize and compare the completely assembled genomes of four clinical isolates of M. felis from three domestic cats, assembled with the aid of short- and long-read sequencing methods. The completed genomes encoded a median of 759 ORFs (range 743-777) and had a median average nucleotide identity of 98.2 % with the genome of the available equid origin reference strain. Comparative genomic analysis revealed the occurrence of multiple horizontal gene transfer events and significant genome reassortment. This had resulted in the acquisition or loss of numerous genes within the Australian felid isolate genomes, encoding putative proteins involved in DNA transfer, metabolism, DNA replication, host cell interaction and restriction modification systems. Additionally, a novel mycoplasma phage was detected in one Australian felid M. felis isolate by genomic analysis and visualized using cryo-transmission electron microscopy. This study has highlighted the complex genomic dynamics in different host environments. Furthermore, the sequences obtained in this work will enable the development of new diagnostic tools, and identification of future infection control and treatment options for the respiratory disease complex in cats.

Xanthomonas rydalmerensis sp. nov., a non-pathogenic member of Group 1 Xanthomonas

Five bacterial isolates were isolated from Fragaria × ananassa in 1976 in Rydalmere, Australia, during routine biosecurity surveillance. Initially, the results of biochemical characterisation indicated that these isolates represented members of the genus Xanthomonas. To determine their species, further analysis was conducted using both phenotypic and genotypic approaches. Phenotypic analysis involved using MALDI-TOF MS and BIOLOG GEN III microplates, which confirmed that the isolates represented members of the genus Xanthomonas but did not allow them to be classified with respect to species. Genome relatedness indices and the results of extensive phylogenetic analysis confirmed that the isolates were members of the genus Xanthomonas and represented a novel species. On the basis the minimal presence of virulence-associated factors typically found in genomes of members of the genus Xanthomonas, we suggest that these isolates are non-pathogenic. This conclusion was supported by the results of a pathogenicity assay. On the basis of these findings, we propose the name Xanthomonas rydalmerensis, with DAR 34855T = ICMP 24941 as the type strain.

Genomic analysis of diverse environmental Acinetobacter isolates identifies plasmids, antibiotic resistance genes, and capsular polysaccharides shared with clinical strains

Acinetobacter baumannii, an important pathogen known for its widespread antibiotic resistance, has been the focus of extensive research within its genus, primarily involving clinical isolates. Consequently, data on environmental A. baumannii and other Acinetobacter species remain limited. Here, we utilized Illumina and Nanopore sequencing to analyze the genomes of 10 Acinetobacter isolates representing 6 different species sourced from aquatic environments in South Australia. All 10 isolates were phylogenetically distinct compared to clinical and other non-clinical Acinetobacter strains, often tens of thousands of single-nucleotide polymorphisms from their nearest neighbors. Despite the genetic divergence, we identified pdif modules (sections of mobilized DNA) carrying clinically important antimicrobial resistance genes in species other than A. baumannii, including carbapenemase oxa58, tetracycline resistance gene tet(39), and macrolide resistance genes msr(E)-mph(E). These pdif modules were located on plasmids with high sequence identity to those circulating in globally distributed A. baumannii ST1 and ST2 clones. The environmental A. baumannii isolate characterized here (SAAb472; ST350) did not possess any native plasmids; however, it could capture two clinically important plasmids (pRAY and pACICU2) with high transfer frequencies. Furthermore, A. baumannii SAAb472 possessed virulence genes and a capsular polysaccharide type analogous to clinical strains. Our findings highlight the potential for environmental Acinetobacter species to acquire and disseminate clinically important antimicrobial resistance genes, underscoring the need for further research into the ecology and evolution of this important genus.IMPORTANCEAntimicrobial resistance (AMR) is a global threat to human, animal, and environmental health. Studying AMR in environmental bacteria is crucial to understand the emergence and dissemination of resistance genes and pathogens, and to identify potential reservoirs and transmission routes. This study provides novel insights into the genomic diversity and AMR potential of environmental Acinetobacter species. By comparing the genomes of aquatic Acinetobacter isolates with clinical and non-clinical strains, we revealed that they are highly divergent yet carry pdif modules that encode resistance to antibiotics commonly used in clinical settings. We also demonstrated that an environmental A. baumannii isolate can acquire clinically relevant plasmids and carries virulence factors similar to those of hospital-associated strains. These findings suggest that environmental Acinetobacter species may serve as reservoirs and vectors of clinically important genes. Consequently, further research is warranted to comprehensively understand the ecology and evolution of this genus.

Ecology- and genome-based identification of the Bifidobacterium adolescentis prototype of the healthy human gut microbiota

Bifidobacteria are among the first microbial colonizers of the human gut, being frequently associated with human health-promoting activities. In the current study, an in silico methodology based on an ecological and phylogenomic-driven approach allowed the selection of a Bifidobacterium adolescentis prototype strain, i.e., B. adolescentis PRL2023, which best represents the overall genetic content and functional features of the B. adolescentis taxon. Such features were confirmed by in vitro experiments aimed at evaluating the ability of this strain to survive in the gastrointestinal tract of the host and its ability to interact with human intestinal cells and other microbial gut commensals. In this context, co-cultivation of B. adolescentis PRL2023 and several gut commensals revealed various microbe-microbe interactions and indicated co-metabolism of particular plant-derived glycans, such as xylan.IMPORTANCEThe use of appropriate bacterial strains in experimental research becomes imperative in order to investigate bacterial behavior while mimicking the natural environment. In the current study, through in silico and in vitro methodologies, we were able to identify the most representative strain of the Bifidobacterium adolescentis species. The ability of this strain, B. adolescentis PRL2023, to cope with the environmental challenges imposed by the gastrointestinal tract, together with its ability to switch its carbohydrate metabolism to compete with other gut microorganisms, makes it an ideal choice as a B. adolescentis prototype and a member of the healthy microbiota of adults. This strain possesses a genetic blueprint appropriate for its exploitation as a candidate for next-generation probiotics.

Genome sequencing unveils blaKPC-2-harboring plasmids as drivers of enhanced resistance and virulence in nosocomial Klebsiella pneumoniae

The threat posed by Klebsiella pneumoniae in healthcare settings has worsened due to the evolutionary advantages conferred by blaKPC-2-harboring plasmids (pKPC-2). However, the specific evolutionary pathway of nosocomial K. pneumoniae carrying pKPC-2 and its transmission between patients and healthcare environments are not yet well understood. Between 1 August and 31 December 2019, 237 ST11 KPC-2-producing-carbapenem-resistant K. pneumoniae (CRKP) (KPC-2-CRKP) were collected from patient or ward environments in an intensive care unit and subjected to Illumina sequencing, of which 32 strains were additionally selected for Nanopore sequencing to obtain complete plasmid sequences. Bioinformatics analysis, conjugation experiments, antimicrobial susceptibility tests, and virulence assays were performed to identify the evolutionary characteristics of pKPC-2. The pKPC-2 plasmids were divided into three subgroups with distinct evolutionary events, including Tn3-mediated plasmid homologous recombination, IS26-mediated horizontal gene transfer, and dynamic duplications of antibiotic resistance genes (ARGs). Surprisingly, the incidence rates of multicopy blaKPC-2, blaSHV-12, and blaCTX-M-65 were quite high (ranging from 27.43% to 67.01%), and strains negative for extended-spectrum β-lactamase tended to develop multicopy blaKPC-2. Notably, the presence of multicopy blaSHV-12 reduced sensitivity to ceftazidime/avibactam (CZA), and the relative expression level of blaSHV-12 in the CZA-resistant group was 6.12 times higher than that in the sensitive group. Furthermore, a novel hybrid pKPC-2 was identified, presenting enhanced virulence levels and decreased susceptibility to CZA. This study emphasizes the notable prevalence of multicopy ARGs and provides a comprehensive insight into the intricate and diverse evolutionary pathways of resistant plasmids that disseminate among patients and healthcare environments.IMPORTANCEThis study is based on a CRKP screening program between patients and ward environments in an intensive care unit, describing the pKPC-2 (blaKPC-2-harboring plasmids) population structure and evolutionary characteristics in clinical settings. Long-read sequencing was performed in genetically closely related strains, enabling the high-resolution analysis of evolution pathway between or within pKPC-2 subgroups. We revealed the extremely high rates of multicopy antibiotic resistance genes (ARGs) in clinical settings and its effect on resistance profile toward novel β-lactam/β-lactamase inhibitor combinations, which belongs to the last line treatment choices toward CRKP infection. A novel hybrid pKPC-2 carrying CRKP with enhanced resistance and virulence level was captured during its clonal spread between patients and ward environment. These evidences highlight the threat of pKPC-2 to CRKP treatment and control. Thus, surveillance and timely disinfection in clinical settings should be practiced to prevent transmission of CRKP carrying threatful pKPC-2. And rational use of antibiotics should be called for to prevent inducing of pKPC-2 evolution, especially the multicopy ARGs.

A new chromosome-level genome assembly and annotation of Cryptosporidium meleagridis

Cryptosporidium spp. are medically and scientifically relevant protozoan parasites that cause severe diarrheal illness in infants and immunosuppressed populations as well as animals. Although most human Cryptosporidium infections are caused by C. parvum and C. hominis, there are several other human-infecting species including C. meleagridis, which is commonly observed in developing countries. Here, we polished and annotated a long-read genome sequence assembly for C. meleagridis TU1867, a species which infects birds and humans. The genome sequence was generated using a combination of whole genome amplification (WGA) and long-read Oxford Nanopore Technologies sequencing. The assembly was then polished with Illumina data. The chromosome-level genome assembly is 9.2 Mbp with a contig N50 of 1.1 Mb. Annotation revealed 3,923 protein-coding genes. A BUSCO analysis indicates a completeness of 96.6% (n=446), including 430 (96.4%) single-copy and 1 (0.224%) duplicated apicomplexan conserved gene(s). The new C. meleagridis genome assembly is nearly gap-free and provides a valuable new resource for the Cryptosporidium community and future studies on evolution and host-specificity.

Whole-genome comparison of two same-genotype macrolide-resistant Bordetella pertussis isolates collected in Japan

The emergence of macrolide-resistant Bordetella pertussis (MRBP) is a significant problem because it reduces treatment options for pertussis and exacerbates the severity and spread of the disease. MRBP has been widely prevalent in mainland China since the 2010s and has been sporadically detected in other Asian countries. In Japan, two MRBP clinical strains were first isolated in Tokyo and Osaka between June and July 2018. The isolates BP616 in Osaka and BP625 in Tokyo harbored the same virulence-associated allelic genes (including ptxP1, ptxA1, prn1, fim3A, and fhaB3) and MT195 genotype and exhibited similar antimicrobial susceptibility profiles. However, despite their simultaneous occurrence, a distinguishable epidemiological link between these isolates could not be established. To gain further insight into the genetic relationship between these isolates in this study, we performed whole-genome analyses. Phylogenetic analysis based on genome-wide single-nucleotide polymorphisms revealed that the isolates belonged to one of the three clades of Chinese MRBP isolates, but there were 11 single-nucleotide polymorphism differences between BP616 and BP625. Genome structure analysis revealed two large inversions (202 and 523 kbp) and one small transposition (3.8 kbp) between the genomes. These findings indicate that the two Japanese MRBP isolates are closely related to Chinese MRBP isolates but are genomically distinct, suggesting that they were introduced into Japan from mainland China through different transmission routes.

The first complete genome sequence of Staphylococcus aureus ST5477

This study presents the first complete genome of Staphylococcus aureus ST5477, one of the most common sequence types (ST) from bovine in eastern Africa. The genome consists of a 2,723,132-bp circular chromosome and a 3,044-bp plasmid. This strain was collected in 2017 from cow milk in Tanzania.

Complete genome sequence of Pseudomonas entomophila strain TVIN-A01

The complete genome sequence of Pseudomonas entomophila strain TVIN-A01 is reported in this manuscript. It is a commonly used laboratory strain of Pseudomonas entomophila known to infect insects and, hence, often studied in host-pathogen interactions. Oxford Nanopore sequencing and Illumina sequencing were performed to assemble the genome completely.

Complete genome sequences of 11 Streptococcus pneumoniae strains isolated from acute infectious purpura fulminans, sepsis, and pneumonia

The complete genome sequences of 11 Japanese Streptococcus pneumoniae isolates were determined by hybrid assembly of long and short reads, including two strains isolated from patients with acute infectious purpura fulminans, six strains from patients with sepsis, and three strains from patients with pneumonia.

Genomic and ecological approaches to identify the Bifidobacterium breve prototype of the healthy human gut microbiota

Members of the genus Bifidobacterium are among the first microorganisms colonizing the human gut. Among these species, strains of Bifidobacterium breve are known to be commonly transmitted from mother to her newborn, while this species has also been linked with activities supporting human wellbeing. In the current study, an in silico approach, guided by ecology- and phylogenome-based analyses, was employed to identify a representative strain of B. breve to be exploited as a novel health-promoting candidate. The selected strain, i.e., B. breve PRL2012, was found to well represent the genetic content and functional genomic features of the B. breve taxon. We evaluated the ability of PRL2012 to survive in the gastrointestinal tract and to interact with other human gut commensal microbes. When co-cultivated with various human gut commensals, B. breve PRL2012 revealed an enhancement of its metabolic activity coupled with the activation of cellular defense mechanisms to apparently improve its survivability in a simulated ecosystem resembling the human microbiome.

The role of hypermutation and collateral sensitivity in antimicrobial resistance diversity of Pseudomonas aeruginosa populations in cystic fibrosis lung infection

Pseudomonas aeruginosa is an opportunistic pathogen which causes chronic, drug-resistant lung infections in cystic fibrosis (CF) patients. In this study, we explore the role of genomic diversification and evolutionary trade-offs in antimicrobial resistance (AMR) diversity within P. aeruginosa populations sourced from CF lung infections. We analyzed 300 clinical isolates from four CF patients (75 per patient) and found that genomic diversity is not a consistent indicator of phenotypic AMR diversity. Remarkably, some genetically less diverse populations showed AMR diversity comparable to those with significantly more genetic variation. We also observed that hypermutator strains frequently exhibited increased sensitivity to antimicrobials, contradicting expectations from their treatment histories. Investigating potential evolutionary trade-offs, we found no substantial evidence of collateral sensitivity among aminoglycoside, beta-lactam, or fluoroquinolone antibiotics, nor did we observe trade-offs between AMR and growth in conditions mimicking CF sputum. Our findings suggest that (i) genomic diversity is not a prerequisite for phenotypic AMR diversity, (ii) hypermutator populations may develop increased antimicrobial sensitivity under selection pressure, (iii) collateral sensitivity is not a prominent feature in CF strains, and (iv) resistance to a single antibiotic does not necessarily lead to significant fitness costs. These insights challenge prevailing assumptions about AMR evolution in chronic infections, emphasizing the complexity of bacterial adaptation during infection.IMPORTANCEUpon infection in the cystic fibrosis (CF) lung, Pseudomonas aeruginosa rapidly acquires genetic mutations, especially in genes involved in antimicrobial resistance (AMR), often resulting in diverse, treatment-resistant populations. However, the role of bacterial population diversity within the context of chronic infection is still poorly understood. In this study, we found that hypermutator strains of P. aeruginosa in the CF lung undergoing treatment with tobramycin evolved increased sensitivity to tobramycin relative to non-hypermutators within the same population. This finding suggests that antimicrobial treatment may only exert weak selection pressure on P. aeruginosa populations in the CF lung. We further found no evidence for collateral sensitivity in these clinical populations, suggesting that collateral sensitivity may not be a robust, naturally occurring phenomenon for this microbe.

Historic methicillin-resistant Staphylococcus aureus: expanding current knowledge using molecular epidemiological characterization of a Swiss legacy collection

BACKGROUND

Few methicillin-resistant Staphylococcus aureus (MRSA) from the early years of its global emergence have been sequenced. Knowledge about evolutionary factors promoting the success of specific MRSA multi-locus sequence types (MLSTs) remains scarce. We aimed to characterize a legacy MRSA collection isolated from 1965 to 1987 and compare it against publicly available international and local genomes.

METHODS

We accessed 451 historic (1965-1987) MRSA isolates stored in the Culture Collection of Switzerland, mostly collected from the Zurich region. We determined phenotypic antimicrobial resistance (AMR) and performed whole genome sequencing (WGS) using Illumina short-read sequencing on all isolates and long-read sequencing on a selection with Oxford Nanopore Technology. For context, we included 103 publicly available international assemblies from 1960 to 1992 and sequenced 1207 modern Swiss MRSA isolates from 2007 to 2022. We analyzed the core genome (cg)MLST and predicted SCCmec cassette types, AMR, and virulence genes.

RESULTS

Among the 451 historic Swiss MRSA isolates, we found 17 sequence types (STs) of which 11 have been previously described. Two STs were novel combinations of known loci and six isolates carried previously unsubmitted MLST alleles, representing five new STs (ST7843, ST7844, ST7837, ST7839, and ST7842). Most isolates (83% 376/451) represented ST247-MRSA-I isolated in the 1960s, followed by ST7844 (6% 25/451), a novel single locus variant (SLV) of ST239. Analysis by cgMLST indicated that isolates belonging to ST7844-MRSA-III cluster within the diversity of ST239-MRSA-III. Early MRSA were predominantly from clonal complex (CC)8. From 1980 to the end of the twentieth century, we observed that CC22 and CC5 as well as CC8 were present, both locally and internationally.

CONCLUSIONS

The combined analysis of 1761 historic and contemporary MRSA isolates across more than 50 years uncovered novel STs and allowed us a glimpse into the lineage flux between Swiss-German and international MRSA across time.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Phenotypes of a Pseudomonas aeruginosa hypermutator lineage that emerged during prolonged mechanical ventilation in a patient without cystic fibrosis

Hypermutator lineages of Pseudomonas aeruginosa arise frequently during the years of airway infection experienced by patients with cystic fibrosis and bronchiectasis but are rare in the absence of chronic infection and structural lung disease. Since the onset of the COVID-19 pandemic, large numbers of patients have remained mechanically ventilated for extended periods of time. These patients are prone to acquire bacterial pathogens that persist for many weeks and have the opportunity to evolve within the pulmonary environment. However, little is known about what types of adaptations occur in these bacteria and whether these adaptations mimic those observed in chronic infections. We describe a COVID-19 patient with a secondary P. aeruginosa lung infection in whom the causative bacterium persisted for >50 days. Over the course of this infection, a hypermutator lineage of P. aeruginosa emerged and co-existed with a non-hypermutator lineage. Compared to the parental lineage, the hypermutator lineage evolved to be less cytotoxic and less virulent. Genomic analyses of the hypermutator lineage identified numerous mutations, including in the mismatch repair gene mutL and other genes frequently mutated in individuals with cystic fibrosis. Together, these findings demonstrate that hypermutator lineages can emerge when P. aeruginosa persists following acute infections such as ventilator-associated pneumonia and that these lineages have the potential to affect patient outcomes.IMPORTANCEPseudomonas aeruginosa may evolve to accumulate large numbers of mutations in the context of chronic infections such as those that occur in individuals with cystic fibrosis. However, these "hypermutator" lineages are rare following acute infections. Here, we describe a non-cystic fibrosis patient with COVID-19 pneumonia who remained mechanically ventilated for months. The patient became infected with a strain of P. aeruginosa that evolved to become a hypermutator. We demonstrate that hypermutation led to changes in cytotoxicity and virulence. These findings are important because they demonstrate that P. aeruginosa hypermutators can emerge following acute infections and that they have the potential to affect patient outcomes in this setting.

Development of a SNP-based strain-identified method for Streptococcus thermophilus CICC 6038 and Lactobacillus delbrueckii ssp. bulgaricus CICC 6047 using pan-genomics analysis

The health benefits conferred by probiotics is specific to individual probiotic strains, highlighting the importance of identifying specific strains for research and production purposes. Streptococcus thermophilus CICC 6038 and Lactobacillus delbrueckii ssp. bulgaricus CICC 6047 are exceedingly valuable for commercial use with an excellent mixed-culture fermentation. To differentiate these 2 strains from other S. thermophilus and L. delbrueckii ssp. bulgaricus, a specific, sensitive, accurate, rapid, convenient, and cost-effective method is required. In this study, we conducted a pan-genome analysis of S. thermophilus and L. delbrueckii ssp. bulgaricus to identify species-specific core genes, along with strain-specific single-nucleotide polymorphisms (SNPs). These genes were used to develop suitable PCR primers, and the conformity of sequence length and unique SNPs was confirmed by sequencing for qualitative identification at the strain level. The results demonstrated that SNPs analysis of PCR products derived from these primers could distinguish CICC 6038 and CICC 6047 accurately and reproducibly from the other strains of S. thermophilus and L. delbrueckii ssp. bulgaricus, respectively. The strain-specific PCR method based on SNPs herein is universally applicable for probiotics identification. It offers valuable insights into identifying probiotics at the strain level that is fit-for-purpose in quality control and compliance assessment of commercial dairy products.

Twelve new metagenome-assembled genomes from non-axenic culture of Griffithsia monilis (Rhodophyta)

We report 12 metagenome-assembled genomes from a non-axenic culture of the red alga Griffithsia monilis Harvey, some of which are distantly related to publicly available genomes.

Complete genome sequence of Bordetella parapertussis strain 400431-b, isolated from a protracted course of whooping cough in Austria, 2023

Here, we report the complete genome of Bordetella parapertussis strain 400431-b isolated from a nasopharyngeal swab from a 4-year-old patient presenting with a protracted course of whooping cough, vaccinated with three doses of diphtheria-tetanus-acellular pertussis-hepatitis B-poliomyelitis-Haemophilus influenzae type b vaccine.

Complete genome sequence of eight LEE-negative Shiga toxin-producing Escherichia coli strains isolated from patients with hemolytic-uremic syndrome

Major serotypes of Shiga toxin-producing Escherichia coli (STEC) carry a locus of enterocyte effacement (LEE), which is required for attaching and effacing lesion formation. Genome information of LEE-negative STEC is scarce despite their virulence potential. We present the complete genome sequences of eight LEE-negative STEC isolates from hemolytic-uremic syndrome patients.

Colicins and T6SS-based competition systems enhance enterotoxigenic E. coli (ETEC) competitiveness

Diarrheal diseases are still a significant problem for humankind, causing approximately half a million deaths annually. To cause diarrhea, enteric bacterial pathogens must first colonize the gut, which is a niche occupied by the normal bacterial microbiota. Therefore, the ability of pathogenic bacteria to inhibit the growth of other bacteria can facilitate the colonization process. Although enterotoxigenic Escherichia coli (ETEC) is one of the major causative agents of diarrheal diseases, little is known about the competition systems found in and used by ETEC and how they contribute to the ability of ETEC to colonize a host. Here, we collected a set of 94 fully assembled ETEC genomes by performing whole-genome sequencing and mining the NCBI RefSeq database. Using this set, we performed a comprehensive search for delivered bacterial toxins and investigated how these toxins contribute to ETEC competitiveness in vitro. We found that type VI secretion systems (T6SS) were widespread among ETEC (n = 47). In addition, several closely related ETEC strains were found to encode Colicin Ia and T6SS (n = 8). These toxins provide ETEC competitive advantages during in vitro competition against other E. coli, suggesting that the role of T6SS as well as colicins in ETEC biology has until now been underappreciated.

KombOver: Efficient k-core and K-truss based characterization of perturbations within the human gut microbiome

The microbes present in the human gastrointestinal tract are regularly linked to human health and disease outcomes. Thanks to technological and methodological advances in recent years, metagenomic sequencing data, and computational methods designed to analyze metagenomic data, have contributed to improved understanding of the link between the human gut microbiome and disease. However, while numerous methods have been recently developed to extract quantitative and qualitative results from host-associated microbiome data, improved computational tools are still needed to track microbiome dynamics with short-read sequencing data. Previously we have proposed KOMB as a de novo tool for identifying copy number variations in metagenomes for characterizing microbial genome dynamics in response to perturbations. In this work, we present KombOver (KO), which includes four key contributions with respect to our previous work: (i) it scales to large microbiome study cohorts, (ii) it includes both k-core and K-truss based analysis, (iii) we provide the foundation of a theoretical understanding of the relation between various graph-based metagenome representations, and (iv) we provide an improved user experience with easier-to-run code and more descriptive outputs/results. To highlight the aforementioned benefits, we applied KO to nearly 1000 human microbiome samples, requiring less than 10 minutes and 10 GB RAM per sample to process these data. Furthermore, we highlight how graph-based approaches such as k-core and K-truss can be informative for pinpointing microbial community dynamics within a myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) cohort. KO is open source and available for download/use at: https://github.com/treangenlab/komb.

Genomic characterization of Mycoplasma edwardii isolated from a dog bite induced cat wound reveals multiple horizontal gene transfer events and loss of the CRISPR/Cas system

A domestic short hair cat (Felis catus) suffering from a purulent wound infection resulting from a dog bite was sampled for bacterial culture and isolation as the wound had been unresponsive to prolonged antimicrobial treatment. A mycoplasma was isolated from the wound. Whole genome sequencing of the isolate was performed using short-read Illumina and long-read Oxford Nanopore chemistry, and the organism was identified as Mycoplasma edwardii. Comparison of the genome sequence of the isolate to a reference M. edwardii genome sequence (canid isolate) identified the loss of several key bacterial factors involved in genome editing, as well the insertion of several novel ORFs most closely related to those found in other canine mycoplasmas, specifically Mycoplasma canis, M. cynos, M. molare and M. maculosa. This is only the second known report of disease caused by M. edwardii in a non-canid species, and the first report of it infecting and causing clinical disease in a cat.

Complete genomes of Clostridium botulinum type B(F) isolates associated with a 1995 foodborne botulism outbreak from commercial pâté reveals a recombination event disrupting the ntnh gene

Foodborne botulism is a neuroparalytic disease caused by ingestion of foods contaminated with botulinum neurotoxin (BoNT), produced by Clostridium botulinum. In 1995 a husband and wife from Québec, Canada, were hospitalized for several months with prolonged muscle paralysis after ingesting a commercial pâté de campagne. Examination of faecal samples from both patients and the pâté produced viable Group I (proteolytic) C. botulinum type B from each of the three samples. Whole genome sequencing revealed that all three isolates contain identical bont/B5 and bont/F2 genes encoded on a plasmid. Both faecal isolate genomes were identical in chromosome and plasmid length, as well as gene content. The genome of the pâté isolate was nearly identical to that of the faecal isolates with the notable difference of a missing 13-gene insertion on the bont/B5 cluster disrupting the ntnh gene. Examination of the insertion revealed several mobile genetic elements that participate in recombination.

Comparison of wild-type KT2440 and genome-reduced EM42 Pseudomonas putida strains for muconate production from aromatic compounds and glucose

Pseudomonas putida KT2440 is a robust, aromatic catabolic bacterium that has been widely engineered to convert bio-based and waste-based feedstocks to target products. Towards industrial domestication of P. putida KT2440, rational genome reduction has been previously conducted, resulting in P. putida strain EM42, which exhibited characteristics that could be advantageous for production strains. Here, we compared P. putida KT2440- and EM42-derived strains for cis,cis-muconic acid production from an aromatic compound, p-coumarate, and in separate strains, from glucose. To our surprise, the EM42-derived strains did not outperform the KT2440-derived strains in muconate production from either substrate. In bioreactor cultivations, KT2440- and EM42-derived strains produced muconate from p-coumarate at titers of 45 g/L and 37 g/L, respectively, and from glucose at 20 g/L and 13 g/L, respectively. To provide additional insights about the differences in the parent strains, we analyzed growth profiles of KT2440 and EM42 on aromatic compounds as the sole carbon and energy sources. In general, the EM42 strain exhibited reduced growth rates but shorter growth lags than KT2440. We also observed that EM42-derived strains resulted in higher growth rates on glucose compared to KT2440-derived strains, but only at the lowest glucose concentrations tested. Transcriptomics revealed that genome reduction in EM42 had global effects on transcript levels and showed that the EM42-derived strains that produce muconate from glucose exhibit reduced modulation of gene expression in response to changes in glucose concentrations. Overall, our results highlight that additional studies are warranted to understand the effects of genome reduction on microbial metabolism and physiology, especially when intended for use in production strains.

Detection of chromosome-mediated blaNDM-1-carrying Aeromonas spp. in the intestinal contents of fresh water river fish in Ho Chi Minh City, Vietnam

The spread of antibiotic-resistant bacteria is a global problem that should be addressed through the perspective of the "one health" concept. The purpose of this study was to determine the contamination rate of antibiotic-resistant Aeromonas spp. in fresh water river fish purchased from a fish market in Vietnam. We then defined the pattern of antibiotic resistance to assess antibiotic-resistant contamination. Antibiotic-resistant Aeromonas spp. were detected in the intestinal contents of 32 of 80 fish. blaNDM-1 was detected in seven strains. Extended-spectrum β-lactamase and AmpC β-lactamase-related genes were detected in 28 strains, including blaCTX-M-55, blaCTX-M-15, blaCTX-M-1, and blaDHA,blaFOX, and blaMOX. The blaNDM-1 detected in the seven Aeromonas spp. strains were found chromosomally. This finding suggests that the blaNDM gene is stable in the natural environment and may spread widely into animals and humans via Aeromonas spp. with a transposon. Our results suggest the importance of continuing to monitor carbapenemase genes in Aeromonas spp. to evaluate the possibility that they may spread in other Enterobacterales, and to elucidate the mechanism of spread.