Assembling the perfect bacterial genome using Oxford Nanopore and Illumina sequencing

Complete genome sequences of Yersinia pestis 6/69 strain isolated from a bubonic plague patient in Madagascar and its isogenic strain cured of pPCP1

We report the complete genome sequences of two valuable strains to investigate plague pathogenesis: (i) Yersinia pestis strain 6/69, which was isolated from a bubonic plague patient in Madagascar and contains pCD1, pMT1, and pPCP1 virulence plasmids, and (ii) the 6/69 strain cured of pPCP1.

Streamlining whole genome sequencing for clinical diagnostics with ONT technology

Recent advances in whole-genome sequencing (WGS) have increased the accessibility of this tool, offering substantial potential for pathogen surveillance, outbreak response, and diagnostics. However, the routine clinical adoption of WGS is hindered by factors such as high costs, technical complexity, and the requirement for bioinformatics expertise for data analysis. To address these challenges, we propose RapidONT, a workflow designed for cost-effective and accessible WGS-based pathogen analysis. RapidONT employs a mechanical shearing-based DNA extraction protocol, followed by library construction by using a multiplexing Oxford nanopore technologies (ONT) rapid barcoding kit. Flye software is used for de novo assembly without manual intervention, followed by basic assembly polishing using Medaka and Homopolish. The polished assemblies are then analyzed using the user-friendly web-based platform Pathogenwatch, which facilitates species identification, molecular typing, and antimicrobial resistance (AMR) prediction, all while requiring minimal bioinformatics expertise. The efficacy of RapidONT was evaluated using nine clinically relevant pathogens, encompassing a total of 90 gram-positive and gram-negative bacterial strains. The workflow demonstrated high accuracy in critical tasks such as multilocus sequence typing (MLST) and AMR identification, using only ONT R9.4.1 flowcell data. Notably, limitations were observed with Salmonella spp. and Neisseria gonorrhoeae. Furthermore, RapidONT enabled the generation of genomic information for 48 bacterial isolates by using a single flow cell, significantly reducing sequencing costs. This approach eliminates the need for extensive experimentation in obtaining crucial genomic information. This workflow facilitates broader WGS implementation in clinical pathogen analysis and diagnostics.

Rapid whole genome characterization of antimicrobial-resistant pathogens using long-read sequencing to identify potential healthcare transmission

OBJECTIVE

Whole genome sequencing (WGS) can help identify transmission of pathogens causing healthcare-associated infections (HAIs). However, the current gold standard of short-read, Illumina-based WGS is labor and time intensive. Given recent improvements in long-read Oxford Nanopore Technologies (ONT) sequencing, we sought to establish a low resource approach providing accurate WGS-pathogen comparison within a time frame allowing for infection prevention and control (IPC) interventions.

METHODS

WGS was prospectively performed on pathogens at increased risk of potential healthcare transmission using the ONT MinION sequencer with R10.4.1 flow cells and Dorado basecaller. Potential transmission was assessed via Ridom SeqSphere+ for core genome multilocus sequence typing and MINTyper for reference-based core genome single nucleotide polymorphisms using previously published cutoff values. The accuracy of our ONT pipeline was determined relative to Illumina.

RESULTS

Over a six-month period, 242 bacterial isolates from 216 patients were sequenced by a single operator. Compared to the Illumina gold standard, our ONT pipeline achieved a mean identity score of Q60 for assembled genomes, even with a coverage rate as low as 40×. The mean time from initiating DNA extraction to complete analysis was 2 days (IQR 2-3.25 days). We identified five potential transmission clusters comprising 21 isolates (8.7% of sequenced strains). Integrating ONT with epidemiological data, >70% (15/21) of putative transmission cluster isolates originated from patients with potential healthcare transmission links.

CONCLUSIONS

Via a stand-alone ONT pipeline, we detected potentially transmitted HAI pathogens rapidly and accurately, aligning closely with epidemiological data. Our low-resource method has the potential to assist in IPC efforts.

Sequencing Strategy to Ensure Accurate Plasmid Assembly

Despite the wide use of plasmids in research and clinical production, the need to verify plasmid sequences is a bottleneck that is too often underestimated in the manufacturing process. Although sequencing platforms continue to improve, the method and assembly pipeline chosen still influence the final plasmid assembly sequence. Furthermore, few dedicated tools exist for plasmid assembly, especially for de novo assembly. Here, we evaluated short-read, long-read, and hybrid (both short and long reads) de novo assembly pipelines across three replicates of a 24-plasmid library. Consistent with previous characterizations of each sequencing technology, short-read assemblies had issues resolving GC-rich regions, and long-read assemblies commonly had small insertions and deletions, especially in repetitive regions. The hybrid approach facilitated the most accurate, consistent assembly generation and identified mutations relative to the reference sequence. Although Sanger sequencing can be used to verify specific regions, some GC-rich and repetitive regions were difficult to resolve using any method, suggesting that easily sequenced genetic parts should be prioritized in the design of new genetic constructs.

Florfenicol administration in piglets co-selects for multiple antimicrobial resistance genes

Antimicrobial use in food-producing animals such as pigs is a significant issue due to its association with antimicrobial resistance. Florfenicol is a broad-spectrum phenicol antibiotic used in swine for various indications; however, its effect on the swine microbiome and resistome is largely unknown. This study investigated these effects in piglets treated intramuscularly with florfenicol at 1 and 7 days of age. Fecal samples were collected from treated (n = 30) and untreated (n = 30) pigs at nine different time points up until 140 days of age, and the fecal metagenomes were sequenced. The fecal microbiomes of the two groups of piglets were most dissimilar in the immediate period following florfenicol administration. These differences were driven in part by an increase in the relative abundance of Clostridium scindens, Enterococcus faecalis, and Escherichia spp. in the florfenicol-treated piglets and Fusobacterium spp., Pauljensenia hyovaginalis, and Ruminococcus gnavus in the control piglets. In addition to selecting for florfenicol resistance genes (floR, fexA, and fexB), florfenicol also selected for genes conferring resistance to the aminoglycosides, beta-lactams, or sulfonamides up until weaning at 21 days of age. Florfenicol-resistant Escherichia coli isolated from these piglets were found to carry a plasmid with floR, along with tet(A), aph(6)-Id, aph(3″)-Ib, sul2, and blaTEM-1/blaCMY-2. A plasmid carrying fexB and poxtA (phenicols and oxazolidinones) was identified in florfenicol-resistant Enterococcus avium, Enterococcus faecium, and E. faecalis isolates from the treated piglets. This study highlights the potential for co-selection and perturbation of the fecal microbial community in pre-weaned piglets administered florfenicol.IMPORTANCEAntimicrobial use remains a serious challenge in food-animal production due to its linkage with antimicrobial resistance. Antimicrobial resistance can reduce the efficacy of veterinary treatment and can potentially be transferred to humans through the food chain or direct contact with animals and their environment. In this study, early-life florfenicol treatment in piglets altered the composition of the fecal microbiome and selected for many unrelated antimicrobial resistance genes up until weaning at 21 days of age. Part of this co-selection process appeared to involve an Escherichia coli plasmid carrying a florfenicol resistance gene along with genes conferring resistance to at least four other antimicrobial classes. In addition, florfenicol selected for certain genes that provide resistance to multiple antimicrobial classes, including the oxazolidinones. These results highlight that florfenicol can co-select for multiple antimicrobial resistance genes, and their presence on mobile genetic elements suggests the potential for transfer to other bacteria.

Rapid Detection of PML::RARA Fusions in Acute Promyelocytic Leukemia: CRISPR/Cas9 Nanopore Sequencing with Adaptive Sampling

Acute promyelocytic leukemia (APL) accounts for approximately 10-15% of newly diagnosed acute myeloid leukemia cases and presents with coagulopathy and bleeding. Prompt diagnosis and treatment are required to minimize early mortality in APL as initiation of all-trans retinoic acid therapy rapidly reverses coagulopathy. The PML::RARA fusion is a hallmark of APL and its rapid identification is essential for rapid initiation of specific treatment to prevent early deaths from coagulopathy and bleeding and optimize patient outcomes. Given limitations and long turnaround time of current gene fusion diagnostic strategies, we have developed a novel amplification-free nanopore sequencing-based approach with low cost, easy setup, and fast turnaround time. We termed the approach CRISPR/Cas9-enriched nanopore sequencing with adaptive sampling (CENAS). Using CENAS, we successfully sequenced breakpoints of typical and atypical PML::RARA fusions in APL patients. Compared with the standard-of-care genetic diagnostic tests, CENAS achieved good concordance in detecting PML::RARA fusions in this study. CENAS allowed for the identification of sequence information of fusion breakpoints involved in typical and atypical PML::RARA fusions and identified additional genes (ANKFN1 and JOSD1) and genomic regions (13q14.13) involving the atypical fusions. To the best of our knowledge, involvements of the ANKFN1 gene, the JOSD1 gene, and the 13q14.13 genomic region flanking with the SIAH3 and ZC3H13 genes have not been reported in the atypical PML::RARA fusions. CENAS has great potential to develop as a point-of-care test enabling immediate, low-cost bedside diagnosis of APL patients with a PML::RARA fusion. Given the early death rate in APL patients still reaches 15%, and ~10% of APL patients are resistant to initial therapy or prone to relapse, further sequencing studies of typical and atypical PML::RARA fusion might shed light on the pathophysiology of the disease and its responsiveness to treatment. Understanding the involvement of additional genes and positional effects related to the PML and RARA genes could shed light on their role in APL and may aid in the development of novel targeted therapies.

Nanopore sequencing: flourishing in its teenage years

Over the past decade, nanopore sequencing has experienced significant advancements and changes, transitioning from an initially emerging technology to a significant instrument in the field of genomic sequencing. However, as advancements in next-generation sequencing technology persist, nanopore sequencing also improves. This paper reviews the developments, applications, and outlook on nanopore sequencing technology. Currently, nanopore sequencing supports both DNA and RNA sequencing, making it widely applicable in areas such as telomere-to-telomere (T2T) genome assembly, direct RNA sequencing (DRS), and metagenomics. The openness and versatility of nanopore sequencing have established it as a preferred option for an increasing number of research teams, signaling a transformative influence on life science research. As the nanopore sequencing technology advances, it provides a faster, more cost-effective approach with extended read lengths, demonstrating the significant potential for complex genome assembly, pathogen detection, environmental monitoring, and human disease research, offering a fresh perspective in sequencing technologies.

Vibrio cholerae lineage and pangenome diversity varies geographically across Bangladesh over one year

Cholera is a diarrhoeal disease caused by Vibrio cholerae. It remains a major public health challenge in the endemic region around the Bay of Bengal. Over decadal time scales, one lineage typically dominates the others and spreads in global pandemic waves. However, it remains unclear to what extent diverse lineages co-circulate during a single outbreak season. Defining the pool of diversity during finer time scales is important because the selective pressures that impact V. cholerae - namely antibiotics and phages - are dynamic on these time scales. To study the nationwide diversity of V. cholerae, we long-read sequenced 273 V. cholerae genomes from seven hospitals over one year (2018) in Bangladesh. Four major V. cholerae lineages were identified: known lineages BD-1, BD-2a, and BD-2b, and a novel lineage that we call BD-3. In 2022, BD-1 caused a large cholera outbreak in Dhaka, apparently outcompeting BD-2 lineages. We show that, in 2018, BD-1 was already dominant in the five northern regions, including Dhaka, consistent with an origin from India in the north. By contrast, we observed a higher diversity of lineages in the two southern regions near the coast. The four lineages differed in pangenome content, including integrative and conjugative elements (ICEs) and genes involved in resistance to bacteriophages and antibiotics. Notably, BD-2a lacked an ICE and is predicted to be more sensitive to phages and antibiotics, but nevertheless persisted throughout the year-long sampling period. Genes associated with antibiotic resistance in V. cholerae from Bangladesh in 2006 were entirely absent from all lineages in 2018-19, suggesting shifting costs and benefits of encoding these genes. Together, our results highlight the dynamic nature of the V. cholerae pangenome and the geographic structure of its lineage diversity.

Bridging Classical Methodologies in Salmonella Investigation with Modern Technologies: A Comprehensive Review

Advancements in genomics and machine learning have significantly enhanced the study of Salmonella epidemiology. Whole-genome sequencing has revolutionized bacterial genomics, allowing for detailed analysis of genetic variation and aiding in outbreak investigations and source tracking. Short-read sequencing technologies, such as those provided by Illumina, have been instrumental in generating draft genomes that facilitate serotyping and the detection of antimicrobial resistance. Long-read sequencing technologies, including those from Pacific Biosciences and Oxford Nanopore Technologies, offer the potential for more complete genome assemblies and better insights into genetic diversity. In addition to these sequencing approaches, machine learning techniques like decision trees and random forests provide powerful tools for pattern recognition and predictive modeling. Importantly, the study of bacteriophages, which interact with Salmonella, offers additional layers of understanding. Phages can impact Salmonella population dynamics and evolution, and their integration into Salmonella genomics research holds promise for novel insights into pathogen control and epidemiology. This review revisits the history of Salmonella and its pathogenesis and highlights the integration of these modern methodologies in advancing our understanding of Salmonella.

Nanopore Sequencing for T-Cell Receptor Rearrangement Analysis in Cutaneous T-Cell Lymphoma

Background: Analysis of T-cell receptor (TCR) clonality is a major diagnostic tool for lymphomas, particularly for cutaneous T-cell lymphomas (CTCL) like Mycosis fungoides and Sézary syndrome. However, a fast and cost-effective workflow is needed to enable widespread use of this method. Methods: We established a procedure for TCR rearrangement analysis via Oxford Nanopore Technology (ONT) sequencing. TCR receptor rearrangements (TCR-gamma and TCR-beta chains) were analyzed in samples from 45 patients with various diagnoses: Mycosis fungoides (37/45), Sézary Syndrome (2/45), folliculotropic CTCL (1/45), and non-CTCL diagnoses as polyclonal controls (5/45). Sample types included formalin-fixed paraffin-embedded (FFPE) samples (27/45), fresh frozen samples (9/45), and CD3-isolated cells (9/45). In addition, DNA of a Jurkat cell line was used as a monoclonal control. TCR amplicons were generated employing an optimized version of the protocol from the Euro Clonality consortium. Sequencing was conducted on the ONT GridION and Illumina MiSeq platforms, followed by similar bioinformatic analysis protocols. The tumor clone frequency (TCF), a crucial prognostic factor for CTCL patients, was used for method comparison. Results: The use of an optimized amplicon protocol and adapted bioinformatic tools demonstrated a strong correlation in TCF values between both sequencing methods across all sample types (range R: 0.992-0.996; range r2: 0.984-0.991). Conclusions: In summary, ONT sequencing was able to detect TCR clonality comparable to NGS, indicating its potential as a faster and more cost-effective option for routine diagnostic use.

Matching excellence: Oxford Nanopore Technologies' rise to parity with Pacific Biosciences in genome reconstruction of non-model bacterium with high G+C content

The reconstruction of complete bacterial genomes is essential for microbial research, offering insights into genetic content, ontology and regulation. While Pacific Biosciences (PacBio) provides high-quality genomes, its cost remains a limitation. Oxford Nanopore Technologies (ONT) offers long reads at a lower cost, yet its error rate raises scepticism. Recent ONT advancements, such as new Flow cells (R10.4.1), chemistry (V14) and duplex mode, improve data quality. Our study compares ONT with PacBio and Illumina, including hybrid data. We used Propionibacterium freudenreichii, a bacterium with a genome known for being difficult to reconstruct. By combining data from ONT's Native Barcoding and a custom-developed BARSEQ method, we achieved high-quality, near-perfect genome assemblies. Our findings demonstrate, for the first time, that the combination of nanopore-only long-native with shorter PCR DNA reads (~3 kb) results in high-quality genome reconstruction, comparable to hybrid data assembly from two sequencing platforms. This endorses ONT as a cost-effective, stand-alone strategy for bacterial genome reconstruction. Additionally, we compared methylated motif detection between PacBio and ONT R10.4.1 data, showing that results comparable to PacBio are achievable using ONT, especially when utilizing the advanced Nanomotif tool.

Optimization of PCA Error Correction Conditions to Improve Efficiency of Virus Genome De Novo Synthesis

In recent years, there have been frequent global outbreaks of viral epidemics such as Zika, COVID-19, and monkeypox, which have had a huge impact on human health and society and have also spurred innovation in virus engineering technology. The rise of synthetic virus genome technology has provided researchers with a new platform to accelerate vaccine and drug development. Although DNA synthesis technology has made significant progress, the current virus genome synthesis technology still requires the assembly of short oligonucleotides of around 60 bp into kb-level lengths when constructing long segments, a process in which the commonly used polymerase chain reaction assembly (PCA) technology has high error rates and is cumbersome to operate. This study optimized the error correction conditions after PCA assembly, increasing the accuracy of synthesizing 1 kb DNA fragments from 4.2 ± 2.1% before error correction to 31.3 ± 3.1% after two rounds of correction, an improvement of over 6 times. This study provides a more efficient operational process for synthesizing virus genomes from scratch, indicating greater potential for virus engineering in epidemic prevention and control and the field of biomedicine.

Benchmarking reveals superiority of deep learning variant callers on bacterial nanopore sequence data

Variant calling is fundamental in bacterial genomics, underpinning the identification of disease transmission clusters, the construction of phylogenetic trees, and antimicrobial resistance detection. This study presents a comprehensive benchmarking of variant calling accuracy in bacterial genomes using Oxford Nanopore Technologies (ONT) sequencing data. We evaluated three ONT basecalling models and both simplex (single-strand) and duplex (dual-strand) read types across 14 diverse bacterial species. Our findings reveal that deep learning-based variant callers, particularly Clair3 and DeepVariant, significantly outperform traditional methods and even exceed the accuracy of Illumina sequencing, especially when applied to ONT's super-high accuracy model. ONT's superior performance is attributed to its ability to overcome Illumina's errors, which often arise from difficulties in aligning reads in repetitive and variant-dense genomic regions. Moreover, the use of high-performing variant callers with ONT's super-high accuracy data mitigates ONT's traditional errors in homopolymers. We also investigated the impact of read depth on variant calling, demonstrating that 10× depth of ONT super-accuracy data can achieve precision and recall comparable to, or better than, full-depth Illumina sequencing. These results underscore the potential of ONT sequencing, combined with advanced variant calling algorithms, to replace traditional short-read sequencing methods in bacterial genomics, particularly in resource-limited settings.

Using Repeated Lysis Steps Fractionates Between Heterotrophic and Cyanobacterial DNA Extracted from Xenic Cyanobacterial Cultures

Extracting DNA from cyanobacteria can be a challenge because of their diverse morphologies, challenging cellular structure, and the heterotrophic microbiome often present within cyanobacterial cultures. As such, even when our DNA yields are sufficient for sequencing, the percentage of reads coming from the cyanobacterial host can be low, leading to incomplete genomes spread across several scaffolds. In this research, we optimized a DNA isolation protocol using three iterative cell lysis steps to enrich the portion of DNA isolated coming from the cyanobacterial host rather than the heterotrophic microbiome. In order to utilize in-house nanopore sequencing, we faced a challenge in that our lysis protocol led to DNA shearing and a lower molecular weight DNA extract than is suitable for this sequencing technology. As such we used two bead-based size selection steps to remove shorter molecules of DNA before nanopore sequencing. EPI2ME analysis of the processed reads from the iterative lysis steps showed that in the first lysis the heterotrophic microbiome could make up more than half of all reads, but with each lysis the proportion of reads coming from these other species decreased. Using our iterative lysis protocol, we were able to sequence the genomes of two cyanobacteria isolated from fresh water sources around northern Mississippi, namely Leptolyngbya sp. BL-A-14 and Limnothrix sp. BL-A-16. The genomes of these isolates were assembled as closed chromosomes of 7.2 and 4.5 Mb for Leptolyngbya sp. BL-A-14 and Limnothrix sp. BL-A-16, respectively. Because some cyanobacteria have symbioses with their heterotrophic microbiome it is not always possible to prepare axenic cultures of these organisms, we hope our approach will be useful for sequencing xenic cultures of cyanobacteria, but we can also imagine applications in studying this microbiome specifically by focusing sequencing efforts on the first fraction.

Isolation and characterization of novel acetogenic Moorella strains for employment as potential thermophilic biocatalysts

Thermophilic acetogenic bacteria have attracted attention as promising candidates for biotechnological applications such as syngas fermentation, microbial electrosynthesis, and methanol conversion. Here, we aimed to isolate and characterize novel thermophilic acetogens from diverse environments. Enrichment of heterotrophic and autotrophic acetogens was monitored by 16S rRNA gene-based bacterial community analysis. Seven novel Moorella strains were isolated and characterized by genomic and physiological analyses. Two Moorella humiferrea isolates showed considerable differences during autotrophic growth. The M. humiferrea LNE isolate (DSM 117358) fermented carbon monoxide (CO) to acetate, while the M. humiferrea OCP isolate (DSM 117359) transformed CO to hydrogen and carbon dioxide (H2 + CO2), employing the water-gas shift reaction. Another carboxydotrophic hydrogenogenic Moorella strain was isolated from the covering soil of an active charcoal burning pile and proposed as the type strain (ACPsT) of the novel species Moorella carbonis (DSM 116161T and CCOS 2103T). The remaining four novel strains were affiliated with Moorella thermoacetica and showed, together with the type strain DSM 2955T, the production of small amounts of ethanol from H2 + CO2 in addition to acetate. The physiological analyses of the novel Moorella strains revealed isolate-specific differences that considerably increase the knowledge base on thermophilic acetogens for future applications.

Recovery of Vandammella animalimorsus from an immunocompetent female patient following cat bite to the lower leg

We report a case of Vandammella animalimorsus in an adult female patient following a cat scratch/bite. Animal bite/scratches may lead to zoonotic transmission of bacteria that subsequently lead to infection. Wound management and antimicrobial therapy is often necessary. The organism was initially misidentified as Neisseria animaloris/zoodegmatis and highlights the difficulty of correctly identifying some bacteria in clinical microbiology laboratories.

Borrelia burgdorferi BB0346 is an Essential, Structurally Variant LolA Homolog that is Primarily Required for Homeostatic Localization of Periplasmic Lipoproteins

In diderm bacteria, the Lol pathway canonically mediates the periplasmic transport of lipoproteins from the inner membrane (IM) to the outer membrane (OM) and therefore plays an essential role in bacterial envelope homeostasis. After extrusion of modified lipoproteins from the IM via the LolCDE complex, the periplasmic chaperone LolA carries lipoproteins through the periplasm and transfers them to the OM lipoprotein insertase LolB, itself a lipoprotein with a LolA-like fold. Yet, LolB homologs appear restricted to γ-proteobacteria and are missing from spirochetes like the tick-borne Lyme disease pathogen Borrelia burgdorferi, suggesting a different hand-off mechanism at the OM. Here, we solved the crystal structure of the B. burgdorferi LolA homolog BB0346 (LolABb) at 1.9 Å resolution. We identified multiple structural deviations in comparative analyses to other solved LolA structures, particularly a unique LolB-like protruding loop domain. LolABb failed to complement an Escherichia coli lolA knockout, even after codon optimization, signal I peptide adaptation, and a C-terminal chimerization which had allowed for complementation with an α-proteobacterial LolA. Analysis of a conditional B. burgdorferi lolA knockout strain indicated that LolABb was essential for growth. Intriguingly, protein localization assays indicated that initial depletion of LolABb led to an emerging mislocalization of both IM and periplasmic OM lipoproteins, but not surface lipoproteins. Together, these findings further support the presence of two separate primary secretion pathways for periplasmic and surface OM lipoproteins in B. burgdorferi and suggest that the distinct structural features of LolABb allow it to function in a unique LolB-deficient lipoprotein sorting system.

Long-Read Sequencing Unlocks New Insights into the Amphidinium carterae Microbiome

Dinoflagellates are one of the largest groups of marine microalgae and exhibit diverse trophic strategies. Some dinoflagellates can produce secondary metabolites that are known to be toxic, which can lead to ecologically harmful blooms. Amphidinium carterae is one species of dinoflagellate that produces toxic compounds and is used as a model for dinoflagellate studies. The impact of the microbiome on A. carterae growth and metabolite synthesis is not yet fully understood, nor is the impact of bacterial data on sequencing and assembly. An antibiotic cocktail was previously shown to eliminate 16S amplification from the dinoflagellate culture. Even with drastically reduced bacterial numbers during antibiotic treatment, bacterial sequences were still present. In this experiment, we used novel Nanopore long-read sequencing techniques on A. carterae cultures to assemble 15 full bacterial genomes ranging from 2.9 to 6.0 Mb and found that the use of antibiotics decreased the percentage of reads mapping back to bacteria. We also identified shifts in the microbiome composition and identified a potentially deleterious bacterial species arising in the absence of the antibiotic treatment. Multiple antibiotic resistance genes were identified, as well as evidence that the bacterial population does not contribute to toxic secondary metabolite synthesis.

A treasure trove of 1034 actinomycete genomes

Filamentous Actinobacteria, recently renamed Actinomycetia, are the most prolific source of microbial bioactive natural products. Studies on biosynthetic gene clusters benefit from or require chromosome-level assemblies. Here, we provide DNA sequences from >1000 isolates: 881 complete genomes and 153 near-complete genomes, representing 28 genera and 389 species, including 244 likely novel species. All genomes are from filamentous isolates of the class Actinomycetia from the NBC culture collection. The largest genus is Streptomyces with 886 genomes including 742 complete assemblies. We use this data to show that analysis of complete genomes can bring biological understanding not previously derived from more fragmented sequences or less systematic datasets. We document the central and structured location of core genes and distal location of specialized metabolite biosynthetic gene clusters and duplicate core genes on the linear Streptomyces chromosome, and analyze the content and length of the terminal inverted repeats which are characteristic for Streptomyces. We then analyze the diversity of trans-AT polyketide synthase biosynthetic gene clusters, which encodes the machinery of a biotechnologically highly interesting compound class. These insights have both ecological and biotechnological implications in understanding the importance of high quality genomic resources and the complex role synteny plays in Actinomycetia biology.

Complete genome sequences of two closely related isolates of Staphylococcus saprophyticus isolated from human fingertips

Complete genomes of two closely related isolates of Staphylococcus saprophyticus from human fingertips, SZ.YL11 and SZ.PL35w, were established through hybrid assembly. Each possesses a single circular chromosome and a circular plasmid, totaling 2,611,553 and 2,611,619 bp, respectively (with G + C 33.14% for both).

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genomic diversity and antimicrobial resistance in clinical Klebsiella pneumoniae isolates from tertiary hospitals in Southern Ghana

OBJECTIVES

Comprehensive data on the genomic epidemiology of hospital-associated Klebsiella pneumoniae in Ghana are scarce. This study investigated the genomic diversity, antimicrobial resistance patterns, and clonal relationships of 103 clinical K. pneumoniae isolates from five tertiary hospitals in Southern Ghana-predominantly from paediatric patients aged under 5 years (67/103; 65%), with the majority collected from urine (32/103; 31%) and blood (25/103; 24%) cultures.

METHODS

We generated hybrid Nanopore-Illumina assemblies and employed Pathogenwatch for genotyping via Kaptive [capsular (K) locus and lipopolysaccharide (O) antigens] and Kleborate (antimicrobial resistance and hypervirulence) and determined clonal relationships using core-genome MLST (cgMLST).

RESULTS

Of 44 distinct STs detected, ST133 was the most common, comprising 23% of isolates (n = 23/103). KL116 (28/103; 27%) and O1 (66/103; 64%) were the most prevalent K-locus and O-antigen types. Single-linkage clustering highlighted the global spread of MDR clones such as ST15, ST307, ST17, ST11, ST101 and ST48, with minimal allele differences (1-5) from publicly available genomes worldwide. Conversely, 17 isolates constituted novel clonal groups and lacked close relatives among publicly available genomes, displaying unique genetic diversity within our study population. A significant proportion of isolates (88/103; 85%) carried resistance genes for ≥3 antibiotic classes, with the blaCTX-M-15 gene present in 78% (n = 80/103). Carbapenem resistance, predominantly due to blaOXA-181 and blaNDM-1 genes, was found in 10% (n = 10/103) of the isolates.

CONCLUSIONS

Our findings reveal a complex genomic landscape of K. pneumoniae in Southern Ghana, underscoring the critical need for ongoing genomic surveillance to manage the substantial burden of antimicrobial resistance.

Phage-Mediated Digestive Decolonization in a Gut-On-A-Chip Model: A Tale of Gut-Specific Bacterial Prosperity

Infections due to antimicrobial-resistant bacteria have become a major threat to global health. Some patients may carry resistant bacteria in their gut microbiota. Specific risk factors may trigger the conversion of these carriages into infections in hospitalized patients. Preventively eradicating these carriages has been postulated as a promising preventive intervention. However, previous attempts at such eradication using oral antibiotics or probiotics have led to discouraging results. Phage therapy, the therapeutic use of bacteriophage viruses, might represent a worthy alternative in this context. Taking inspiration from this clinical challenge, we built Gut-On-A-Chip (GOAC) models, which are tridimensional cell culture models mimicking a simplified gut section. These were used to better understand bacterial dynamics under phage pressure using two relevant species: Pseudomonas aeruginosa and Escherichia coli. Model mucus secretion was documented by ELISA assays. Bacterial dynamics assays were performed in GOAC triplicates monitored for 72 h under numerous conditions, such as pre-, per-, or post-bacterial timing of phage introduction, punctual versus continuous phage administration, and phage expression of mucus-binding properties. The potential genomic basis of bacterial phage resistance acquired in the model was investigated by variant sequencing. The bacterial "escape growth" rates under phage pressure were compared to static in vitro conditions. Our results suggest that there is specific bacterial prosperity in this model compared to other in vitro conditions. In E. coli assays, the introduction of a phage harboring unique mucus-binding properties could not shift this balance of power, contradicting previous findings in an in vivo mouse model and highlighting the key differences between these models. Genomic modifications were correlated with bacterial phage resistance acquisition in some but not all instances, suggesting that alternate ways are needed to evade phage predation, which warrants further investigation.

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.

Sequencing Strategy to Ensure Accurate Plasmid Assembly

Despite the wide use of plasmids in research and clinical production, the need to verify plasmid sequences is a bottleneck that is too often underestimated in the manufacturing process. Although sequencing platforms continue to improve, the method and assembly pipeline chosen still influence the final plasmid assembly sequence. Furthermore, few dedicated tools exist for plasmid assembly, especially for de novo assembly. Here, we evaluated short-read, long-read, and hybrid (both short and long reads) de novo assembly pipelines across three replicates of a 24-plasmid library. Consistent with previous characterizations of each sequencing technology, short-read assemblies had issues resolving GC-rich regions, and long-read assemblies commonly had small insertions and deletions, especially in repetitive regions. The hybrid approach facilitated the most accurate, consistent assembly generation and identified mutations relative to the reference sequence. Although Sanger sequencing can be used to verify specific regions, some GC-rich and repetitive regions were difficult to resolve using any method, suggesting that easily sequenced genetic parts should be prioritized in the design of new genetic constructs.

How low can you go? Short-read polishing of Oxford Nanopore bacterial genome assemblies

It is now possible to assemble near-perfect bacterial genomes using Oxford Nanopore Technologies (ONT) long reads, but short-read polishing is usually required for perfection. However, the effect of short-read depth on polishing performance is not well understood. Here, we introduce Pypolca (with default and careful parameters) and Polypolish v0.6.0 (with a new careful parameter). We then show that: (1) all polishers other than Pypolca-careful, Polypolish-default and Polypolish-careful commonly introduce false-positive errors at low read depth; (2) most of the benefit of short-read polishing occurs by 25× depth; (3) Polypolish-careful almost never introduces false-positive errors at any depth; and (4) Pypolca-careful is the single most effective polisher. Overall, we recommend the following polishing strategies: Polypolish-careful alone when depth is very low (<5×), Polypolish-careful and Pypolca-careful when depth is low (5-25×), and Polypolish-default and Pypolca-careful when depth is sufficient (>25×).

Comprehensive genomic and plasmid characterization of multidrug-resistant bacterial strains by R10.4.1 nanopore sequencing

The escalating prevalence of multidrug-resistant (MDR) bacteria pose a significant public health threat. Understanding the genomic features and deciphering the antibiotic resistance profiles of these pathogens is crucial for development of effective surveillance and treatment strategies. In this study, we employed the R10.4.1 nanopore sequencing technology, specifically through the use of the MinION platform, to analyze eight MDR bacterial strains originating from clinical, ecological and food sources. A single 72-hour sequencing run could yield approximately 12 million reads which covered a total of 34 gigabases (Gbp). The nanopore R10.4.1 data was processed using the Flye assembler, successfully assembling the genomes of eight bacterial strains and their 18 plasmids. Notably, the assemblies generated solely from R10.4.1 nanopore data closely matched those from next-generation sequencing data. Diverse antibiotic resistance patterns and specific resistance genes in the test strains were identified. Hospital strains that exhibited multidrug resistance were found to harbor various resistance genes that encode efflux pumps and extended-spectrum β-lactamases. Environmental and food sources were found to display resistance profiles in a species-specific manner. The composition of structurally complex plasmids in the test strains could also be revealed by analysis of nanopore long reads, which also suggested evidence of horizontal transfer of plasmids between different bacterial species. These findings provide valuable insights into the genetic characteristics of MDR bacteria and demonstrating the practicality of nanopore sequencing technology for detecting of resistance elements in bacterial pathogens.

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.

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.

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. They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance 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.

Many purported pseudogenes in bacterial genomes are bona fide genes

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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).

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.

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.

Deciphering the genomes of motility-deficient mutants of Vibrio alginolyticus 138-2

The motility of Vibrio species plays a pivotal role in their survival and adaptation to diverse environments and is intricately associated with pathogenicity in both humans and aquatic animals. Numerous mutant strains of Vibrio alginolyticus have been generated using UV or EMS mutagenesis to probe flagellar motility using molecular genetic approaches. Identifying these mutations promises to yield valuable insights into motility at the protein structural physiology level. In this study, we determined the complete genomic structure of 4 reference specimens of laboratory V. alginolyticus strains: a precursor strain, V. alginolyticus 138-2, two strains showing defects in the lateral flagellum (VIO5 and YM4), and one strain showing defects in the polar flagellum (YM19). Subsequently, we meticulously ascertained the specific mutation sites within the 18 motility-deficient strains related to the polar flagellum (they fall into three categories: flagellar-deficient, multi-flagellar, and chemotaxis-deficient strains) by whole genome sequencing and mapping to the complete genome of parental strains VIO5 or YM4. The mutant strains had an average of 20.6 (±12.7) mutations, most of which were randomly distributed throughout the genome. However, at least two or more different mutations in six flagellar-related genes were detected in 18 mutants specifically selected as chemotaxis-deficient mutants. Genomic analysis using a large number of mutant strains is a very effective tool to comprehensively identify genes associated with specific phenotypes using forward genetics.

The diagnosis of central nervous system infections in resource-limited settings and the use of novel and molecular diagnostic platforms to improve diagnosis

INTRODUCTION

Central nervous system infections (CNSI) disproportionately affect individuals in low-resource settings where diagnosis is challenging; large proportions of patients never receive a confirmed microbiological diagnosis resulting in inadequate management and high mortality. The epidemiology of CNSI varies globally and conventional diagnostics deployed in resource-limited settings have significant limitations, with an urgent need for improved diagnostic strategies.

AREAS COVERED

This review describes molecular platforms and other novel diagnostics used in the diagnosis of CNSI that are applicable to resource-limited settings. An extensive literature search of Medline and PubMed was performed. The emphasis is on investigations targeting infections of relevance to resource-limited settings either due to variation in regional CNSI epidemiology or due to increased prevalence in patients with immunosuppression. This includes commercially available multiplex PCR platforms, mycobacterial PCR platforms, and rapid diagnostics tests. To offer a framework for the optimal implementation in clinical settings, existing evidence highlighting the advantages and limitations of available platforms is reviewed.

EXPERT OPINION

The implementation of molecular platforms and other novel diagnostics has the potential to transform CNSI diagnosis in resource-limited settings, with several examples of successful rollout of novel diagnostics such as Xpert MTB/RIF Ultra and cryptococcal antigen testing.

A recently collected Xanthomonas translucens isolate encodes TAL effectors distinct from older, less virulent isolates

Xanthomonas translucens, the causal agent of bacterial leaf streak disease (BLS) in cereals, is a re-emerging pathogen that is becoming increasingly destructive across the world. While BLS has caused yield losses in the past, there is anecdotal evidence that newer isolates may be more virulent. We observed that two X. translucens isolates collected from two sites in Colorado, USA, are more aggressive on current wheat and barley varieties compared to older isolates, and we hypothesize that genetic changes between recent and older isolates contribute to the differences in isolate aggressiveness. To test this, we phenotyped and genetically characterized two X. translucens isolates collected from Colorado in 2018, which we designated CO236 (from barley) and CO237 (from wheat). Using pathovar-specific phenotyping and PCR primers, we determined that CO236 belongs to pathovar translucens (Xtt) and CO237 belongs to pathovar undulosa (Xtu). We sequenced the full genomes of the isolates using Oxford Nanopore long-read sequencing, and compared their whole genomes against published X. translucens genomes. This analysis confirmed our pathovar designations for Xtt CO236 and Xtu CO237, and showed that, at the whole-genome level, there were no obvious genomic structural changes between Xtt CO236 and Xtu CO237 and other respective published pathovar genomes. Focusing on pathovar undulosa (Xtu CO237), we then compared putative type III effectors among all available Xtu isolate genomes and found that they were highly conserved. However, there were striking differences in the presence and sequence of various transcription activator-like effectors between Xtu CO237 and published undulosa genomes, which correlate with isolate virulence. Here, we explore the potential implications of the differences in these virulence factors, and provide possible explanations for the increased virulence of recently emerged isolates.

Characterization of Rosculus vilicus sp. nov., a rhizarian amoeba interacting with Mycobacterium avium subsp. paratuberculosis

Free-living amoebae are described as potential reservoirs for pathogenic bacteria in the environment. It has been hypothesized that this might be the case for Mycobacterium avium subsp. paratuberculosis, the bacterium responsible for paratuberculosis. In a previous work, we isolated an amoeba from a water sample in the environment of infected cattle and showed that this amoeba was associated with Mycobacterium avium subsp. paratuberculosis. While a partial 18S rRNA gene has allowed us to suggest that this amoeba was Rosculus-like, at that time we were not able to sub-cultivate it. In the present study, we succeeded in cultivating this strain at 20-25°C. This amoeba is among the smallest (5-7 μm) described. The sequencing of the whole genome allowed us to extract the full 18S rRNA gene and propose this strain as a new species of the Rosculus genus, i.e., R. vilicus. Of note, the mitochondrial genome is particularly large (184,954 bp). Finally, we showed that this amoeba was able to phagocyte Mycobacterium avium subsp. paratuberculosis and that the bacterium was still observed within amoebae after at least 3 days. In conclusion, we characterized a new environmental amoeba species at the cellular and genome level that was able to interact with Mycobacterium avium subsp. paratuberculosis. As a result, R. vilicus is a potential candidate as environmental reservoir for Mycobacterium avium subsp. paratuberculosis but further experiments are needed to test this hypothesis.

Genome mining for macrolactam-encoding gene clusters allowed for the network-guided isolation of β-amino acid-containing cyclic derivatives and heterologous production of ciromicin A

β-Amino acid-containing macrolactams represent a structurally diverse group of bioactive natural products derived from polyketides; however we are currently lacking a comprehensive overview about their abundance across bacterial families and the underlying biosynthetic diversity. In this study, we employed a targeted β-amino acid-specific homology-based multi-query search to identify potential bacterial macrolactam producers. Here we demonstrate that approximately 10% of each of the identified actinobacterial genera harbor a biosynthetic gene cluster (BGC) encoding macrolactam production. Based on our comparative study, we propose that mutations occurring in specific regions of polyketide synthases (PKS) are the primary drivers behind the variation in macrolactam ring sizes. We successfully validated two producers of ciromicin A from the genus Amycolatopsis, revised the composition of the biosynthetic gene cluster region mte of macrotermycins, and confirmed the ciromicin biosynthetic pathway through heterologous expression. Additionally, network-based metabolomic analysis uncovered three previously unreported macrotermycin congeners from Amycolatopsis sp. M39. The combination of targeted mining and network-based analysis serves as a powerful tool for identifying macrolactam producers and our studies will catalyze the future discovery of yet unreported macrolactams.

Complete genome sequence of Yersinia pseudotuberculosis strain SP-1303 from lineage 8, associated with Far East scarlet-like fever

We report the complete genome sequence of Yersinia pseudotuberculosis strain SP-1303, identified as part of lineage 8 and associated with Far East scarlet-like fever. The genome includes the chromosome, the Yersinia-virulence plasmid (pYV) encoding a type III secretion system essential for virulence, the pVM82 plasmid, and two cryptic plasmids.

The intra-host evolutionary landscape and pathoadaptation of persistent Staphylococcus aureus in chronic rhinosinusitis

Chronic rhinosinusitis (CRS) is a common chronic sinonasal mucosal inflammation associated with Staphylococcus aureus biofilm and relapsing infections. This study aimed to determine rates of S. aureus persistence and pathoadaptation in CRS patients by investigating the genomic relatedness and antibiotic resistance/tolerance in longitudinally collected S. aureus clinical isolates. A total of 68 S. aureus paired isolates (34 pairs) were sourced from 34 CRS patients at least 6 months apart. Isolates were grown into 48 h biofilms and tested for tolerance to antibiotics. A hybrid sequencing strategy was used to obtain high-quality reference-grade assemblies of all isolates. Single nucleotide variants (SNV) divergence in the core genome and sequence type clustering were used to analyse the relatedness of the isolate pairs. Single nucleotide and structural genome variations, plasmid similarity, and plasmid copy numbers between pairs were examined. Our analysis revealed that 41 % (14/34 pairs) of S. aureus isolates were persistent, while 59 % (20/34 pairs) were non-persistent. Persistent isolates showed episode-specific mutational changes over time with a bias towards events in genes involved in adhesion to the host and mobile genetic elements such as plasmids, prophages, and insertion sequences. Furthermore, a significant increase in the copy number of conserved plasmids of persistent strains was observed. This was accompanied by a significant increase in biofilm tolerance against all tested antibiotics, which was linked to a significant increase in biofilm biomass over time, indicating a potential biofilm pathoadaptive process in persistent isolates. In conclusion, our study provides important insights into the mutational changes during S. aureus persistence in CRS patients highlighting potential pathoadaptive mechanisms in S. aureus persistent isolates culminating in increased biofilm biomass.

Genomic Analysis of Two Cold-Active Pseudoalteromonas Phages Isolated from the Continental Shelf in the Arctic Ocean

Cold-active bacteriophages are bacterial viruses that infect and replicate at low temperatures (≤4 °C). Understanding remains limited of how cold-active phage-host systems sustain high viral abundance despite the persistently low temperatures in pelagic sediments in polar seas. In this study, two Pseudoalteromonas phages, ACA1 and ACA2, were isolated from sediment core samples of the continental shelf in the western Arctic Ocean. These phages exhibited successful propagation at a low temperature of 1 °C and displayed typical myovirus morphology with isometric icosahedral heads and contractile tails. The complete genome sequences of phages ACA1 and ACA2 were 36,825 bp and 36,826 bp in size, respectively, sharing almost the same gene content. These are temperate phages encoding lysogeny-related proteins such as anti-repressor, immunity repressor and integrase. The absence of cross-infection between the host strains, which were genomically distinct Pseudoalteromonas species, can likely be attributed to heavy divergence in the anti-receptor apparently mediated by an associated diversity-generating retroelement. HHpred searching identified genes for all of the structural components of a P2-like phage (family Peduoviridae), although the whole of the Peduoviridae family appeared to be divided between two anciently diverged tail modules. In contrast, Blast matching and whole genome tree analysis are dominated by a nonstructural gene module sharing high similarity with Pseudoalteromonas phage C5a (founder of genus Catalunyavirus). This study expands the knowledge of diversity of P2-like phages known to inhabit Peudoalteromonas and demonstrates their presence in the Arctic niche.

De novo genome assembly resolving repetitive structures enables genomic analysis of 35 European Mycoplasmopsis bovis strains

Mycoplasmopsis (M.) bovis, the agent of mastitis, pneumonia, and arthritis in cattle, harbors a small genome of approximately 1 Mbp. Combining data from Illumina and Nanopore technologies, we sequenced and assembled the genomes of 35 European strains and isolate DL422_88 from Cuba. While the high proportion of repetitive structures in M. bovis genomes represent a particular challenge, implementation of our own pipeline Mycovista (available on GitHub www.github.com/sandraTriebel/mycovista ) in a hybrid approach enabled contiguous assembly of the genomes and, consequently, improved annotation rates considerably. To put our European strain panel in a global context, we analyzed the new genome sequences together with 175 genome assemblies from public databases. Construction of a phylogenetic tree based on core genes of these 219 strains revealed a clustering pattern according to geographical origin, with European isolates positioned on clades 4 and 5. Genomic data allowing assignment of strains to tissue specificity or certain disease manifestations could not be identified. Seven strains isolated from cattle with systemic circular condition (SCC), still a largely unknown manifestation of M. bovis disease, were located on both clades 4 and 5. Pairwise association analysis revealed 108 genomic elements associated with a particular clade of the phylogenetic tree. Further analyzing these hits, 25 genes are functionally annotated and could be linked to a M. bovis protein, e.g. various proteases and nucleases, as well as ten variable surface lipoproteins (Vsps) and other surface proteins. These clade-specific genes could serve as useful markers in epidemiological and clinical surveys.

IS1-related large-scale deletion of chromosomal regions harbouring the oxygen-insensitive nitroreductase gene nfsB causes nitrofurantoin heteroresistance in Escherichia coli

Nitrofurantoin is a broad-spectrum first-line antimicrobial used for managing uncomplicated urinary tract infection (UTI). Loss-of-function mutations in chromosomal genes nfsA, nfsB and ribE of Escherichia coli are known to reduce nitrofurantoin susceptibility. Here, we report the discovery of nitrofurantoin heteroresistance in E. coli clinical isolates and a novel genetic mechanism associated with this phenomenon. Subpopulations with lower nitrofurantoin susceptibility than major populations (hereafter, nitrofurantoin-resistant subpopulations) in two E. coli blood isolates (previously whole-genome sequenced) were identified using population analysis profiling. Each isolate was known to have a loss-of-function mutation in nfsA. From each isolate, four nitrofurantoin-resistant isolates were derived at a nitrofurantoin concentration of 32 mg l-1, and a comparator isolate was obtained without any nitrofurantoin exposure. Genomes of derived isolates were sequenced on Illumina and Nanopore MinION systems. Genetic variation between isolates was determined based on genome assemblies and read mapping. Nitrofurantoin minimum inhibitory concentrations (MICs) of both blood isolates were 64 mg l-1, with MICs of major nitrofurantoin-susceptible populations varying from 4 to 8 mg l-1. Two to 99 c.f.u. per million demonstrated growth at the nitrofurantoin concentration of 32 mg l-1, which is distinct from that of a homogeneously susceptible or resistant isolate. Derived nitrofurantoin-resistant isolates had 11-66 kb deletions in chromosomal regions harbouring nfsB, and all deletions were immediately adjacent to IS1-family insertion sequences. Our findings demonstrate that the IS1-associated large-scale genetic deletion is a hitherto unrecognized mechanism of nitrofurantoin heteroresistance and could compromise UTI management. Further, frequencies of resistant subpopulations from nitrofurantoin-heteroresistant isolates may challenge conventional nitrofurantoin susceptibility testing in clinical settings.

Sequencing and characterization of Helcococcus ovis: a comprehensive comparative genomic analysis of virulence

BACKGROUND

Helcococcus ovis (H. ovis) is an emerging bacterial pathogen that commonly causes opportunistic respiratory, mammary, and uterine infections across mammalian hosts. This study applied long- and short-read whole genome sequencing technologies to identify virulence factors in five H. ovis isolates with low, medium, and high virulence phenotypes.

RESULTS

The resulting assemblies contained one circular chromosome ranging from 1,744,566 to 1,850,083 bp in length and had a mean GC content of 27.6%. Phylogenetic and nucleotide identity analyses found low virulence strain KG38 to be part of a clade that forms an outgroup apart from the rest of the H. ovis taxon. Assembling the first complete genomes of the species revealed major genomic rearrangements in KG38. One to six prophage regions were identified in each genome. A novel pathogenicity island was found exclusively in the two high virulence strains (KG37 and KG104), along with two hypothetical transmembrane proteins designated as putative VFs. Finally, three zinc ABC transporters and three Type-II/IV secretion systems were identified as possible virulence determinants in this species. The low virulence strain KG38 has fewer intact paralogs of these operons in its genome compared to the higher virulence isolates, which strongly suggests a role in virulence. This strain is also missing four putative virulence factors (VFs) found in other isolates associated with adherence (collagen adhesin precursor), immune evasion (choline-binding protein A and a PspA-like hypothetical protein) and cell wall synthesis (glycerol-3-phosphate cytidylyltransferase).

CONCLUSIONS

In this study, we assembled reference-quality complete genomes for five H. ovis strains to identify putative virulence factors. Phylogenetic analyses of H. ovis isolates revealed the presence of a clade representing a potentially novel species within the genus Helcococcus. A novel pathogenicity island and two hypothetical transmembrane proteins were found exclusively in high-virulence strains. The identification of Zinc ABC transporters and Type-II/IV secretion systems as possible virulence determinants, along with the differences in operon content between the low and high virulence isolates, strongly suggests they also play a role in the bacterium's pathogenicity. Taken together, these findings are a valuable first step toward deciphering the pathogenesis of H. ovis infections.

Genomic dissection of endemic carbapenem resistance reveals metallo-beta-lactamase dissemination through clonal, plasmid and integron transfer

Infections caused by metallo-beta-lactamase-producing organisms (MBLs) are a global health threat. Our understanding of transmission dynamics and how MBLs establish endemicity remains limited. We analysed two decades of blaIMP-4 evolution in a hospital using sequence data from 270 clinical and environmental isolates (including 169 completed genomes) and identified the blaIMP-4 gene across 7 Gram-negative genera, 68 bacterial strains and 7 distinct plasmid types. We showed how an initial multi-species outbreak of conserved IncC plasmids (95 genomes across 37 strains) allowed endemicity to be established through the ability of blaIMP-4 to disseminate in successful strain-genetic setting pairs we termed propagators, in particular Serratia marcescens and Enterobacter hormaechei. From this reservoir, blaIMP-4 persisted through diversification of genetic settings that resulted from transfer of blaIMP-4 plasmids between bacterial hosts and of the integron carrying blaIMP-4 between plasmids. Our findings provide a framework for understanding endemicity and spread of MBLs and may have broader applicability to other carbapenemase-producing organisms.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

A hybrid and poly-polish workflow for the complete and accurate assembly of phage genomes: a case study of ten przondoviruses

Bacteriophages (phages) within the genus Przondovirus are T7-like podoviruses belonging to the subfamily Studiervirinae, within the family Autographiviridae, and have a highly conserved genome organisation. The genomes of these phages range from 37 to 42 kb in size, encode 50-60 genes and are characterised by the presence of direct terminal repeats (DTRs) flanking the linear chromosome. These DTRs are often deleted during short-read-only and hybrid assemblies. Moreover, long-read-only assemblies are often littered with sequencing and/or assembly errors and require additional curation. Here, we present the isolation and characterisation of ten novel przondoviruses targeting Klebsiella spp. We describe HYPPA, a HYbrid and Poly-polish Phage Assembly workflow, which utilises long-read assemblies in combination with short-read sequencing to resolve phage DTRs and correcting errors, negating the need for laborious primer walking and Sanger sequencing validation. Our assembly workflow utilised Oxford Nanopore Technologies for long-read sequencing for its accessibility, making it the more relevant long-read sequencing technology at this time, and Illumina DNA Prep for short-read sequencing, representing the most commonly used technologies globally. Our data demonstrate the importance of careful curation of phage assemblies before publication, and prior to using them for comparative genomics.

Plassembler: an automated bacterial plasmid assembly tool

SUMMARY

With recent advances in sequencing technologies, it is now possible to obtain near-perfect complete bacterial chromosome assemblies cheaply and efficiently by combining a long-read-first assembly approach with short-read polishing. However, existing methods for assembling bacterial plasmids from long-read-first assemblies often misassemble or even miss bacterial plasmids entirely and accordingly require manual curation. Plassembler was developed to provide a tool that automatically assembles and outputs bacterial plasmids using a hybrid assembly approach. It achieves increased accuracy and computational efficiency compared to the existing gold standard tool Unicycler by removing chromosomal reads from the input read sets using a mapping approach.

AVAILABILITY AND IMPLEMENTATION

Plassembler is implemented in Python and is installable as a bioconda package using 'conda install -c bioconda plassembler'. The source code is available on GitHub at https://github.com/gbouras13/plassembler. The full benchmarking pipeline can be found at https://github.com/gbouras13/plassembler_simulation_benchmarking, while the benchmarking input FASTQ and output files can be found at https://doi.org/10.5281/zenodo.7996690.