Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes

Multidrug-resistant Shigella flexneri outbreak affecting humans and non-human primates in New Mexico, USA

Shigellosis is a gastrointestinal infection caused by species of Shigella. A large outbreak of Shigella flexneri serotype 2a occurred in Albuquerque, New Mexico between May 2021 and November 2023 that involved humans and non-human primates (NHP) from a local zoo. We analyzed the genomes of 202 New Mexican isolates as well as 15 closely related isolates from other states, and four from NHP. The outbreak was initially detected within men who have sex with men but then predominantly affected people experiencing homelessness. Nearly 70% of cases were hospitalized and there was one human death. The outbreak extended into Albuquerque's BioPark Zoo, causing high morbidity and six deaths in NHPs. All isolates were multidrug-resistant, including towards fluoroquinolones, a first line treatment option which led to treatment failures in human and NHP populations. We show the circulation of the same S. flexneri strain in humans and NHPs, causing fatalities in both populations. This study demonstrates the threat of antimicrobial resistant organisms to vulnerable human and NHP populations and emphasizes the value of genomic surveillance within a One Health framework.

Genomic analysis of community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) causing infections in children-a Spanish multicenter study

Objectives

Staphylococcus aureus is one of the most common human pathogens causing skin and soft tissue infections (SSTIs) among children. This study investigated the molecular traits of community-associated methicillin-resistant S. aureus (CA-MRSA) isolates causing infections in children in Spain.

Methods

Antibiotic susceptibility testing and whole-genome sequencing were performed in 98 CA-MRSA isolates (4.2 median age, 52% males). The phylogenetic relationship, antibiotic resistance, virulence, and plasmid replicon genes content were investigated.

Results

Resistance rates were found as follows: Erythromycin, 42.9%, which could be explained due to the presence of erm(C), mph(C), and msr(A) genes; tobramycin, 27.5%, which could be explained due to the presence of aac(6')-Ie/aph(2″)-Ia and aadD1 genes; tetracycline, 25.5%, which could be explained mainly due to the presence of tet(K) genes; levofloxacin and moxifloxacin, 19.4%, which could be explained primarily due to the mutations in gyrA and parC genes; and gentamicin, 15.3%, which could be explained due to the presence of aac(6')-Ie/aph(2″)-Ia gene. The most prevalent lineage was ST8-IVc and t008. Most isolates were genetically diverse, except for three groups of isolates from the same hospital and one group of isolates from different hospitals. These had less than or equal to 5 allele differences by core-genome multilocus sequence typing (cgMLST) analysis or 0-6 core single-nucleotide polymorphisms (SNPs) by core-genome SNP-based analysis. Phage-encoded Panton-Valentine leukocidin (PVL) genes were found in 75.5% of the isolates. Other common virulence genes were related to adhesion (capA and capP), lipid degradation (geh), hemolysis (hlb, hld, hlgABC, and hly/hla), and tissue destruction (sspAB).

Conclusion

This study observed a high genetic diversity among CA-MRSA isolates causing community-acquired infections in children in Spain, with ST8-IVc as the most prevalent lineage. Nevertheless, genetic relatedness of some isolates from the same as well as different hospitals suggests the dissemination of CA-MRSA among children by contact.

Temporality and Genetic Relatedness of Salmonella in a Pork Processing Facility

The goal of this study was to investigate the prevalence and genetic relatedness of Salmonella enterica in meat and contact surfaces from two processing lines at a pork processing plant over a commercial production schedule. Across 192 samples, there was no significant difference in Salmonella prevalence between Bootjack Trim (BJ) and Boston Butt Trim (BBT) meat (11.5% vs. 11.5%, P = 1.0), though prevalence was higher in meat than on contact surfaces for both the BJ (11.5% vs. 0%, P = 0.01) and BBT (11.5% vs. 3.1%, P = 0.08) processing lines. Both Salmonella prevalence and identified serotypes clustered within four distinct processing windows that spanned multiple dates and processing lines. Phylogenetic analysis using core single nucleotide polymorphisms (SNPs) identified a highly related Salmonella I4,[5],12:i:- strain (N = 33, 0-2 SNPs difference across all isolates) in both the BJ and BBT lines, persisting over consecutive days within one processing window. Similarly, a highly related Salmonella London strain (N = 18, 0-1 SNPs) was found across both processing lines on three processing dates that spanned 28 days. Additional highly related strains of Salmonella Typhimurium (N = 8, 0-1 SNPs) and Salmonella Agona (N = 7, 0-3 SNPs) were also detected across multiple dates. Strains of S. I4,[5],12:i:- and S. London were genetically distinct (>30 SNPs) from publicly available genomes from isolates obtained from other pork processing plants located in the Upper Midwest. Overall, findings suggested that Salmonella prevalence varies across processing lines and production schedules. However, the high phylogenetic relatedness among the Salmonella serotypes suggests a common source may have been present prior to each primal cut being processed into subprimal cuts.

Identification of Novel Staphylococcus aureus Core and Accessory Virulence Patterns in Chronic Rhinosinusitis

Staphylococcus aureus (S. aureus) colonizes the nasal cavities of both healthy individuals and patients with chronic rhinosinusitis (CRS) with (CRSwNP) and without (CRSsNP) nasal polyps. Treatment-resistant S. aureus biofilms and intracellular persistence are common in CRS patients, requiring the expression of specific virulence factor genes to transition into these forms. We hypothesized that S. aureus isolates from non-diseased controls, CRSsNP patients, and CRSwNP patients would exhibit distinct virulence factor patterns contributing to persistence and intracellular survival in CRS patients. Nasal swabs from seventy-seven individuals yielded S. aureus cultures in eight non-diseased controls, eight CRSsNP patients, and five CRSwNP patients. Whole-genome sequencing analyzed stress, antimicrobial resistance, and virulence genes, including plasmids and prophages. Four virulence factor gene patterns emerged: a core set (hlgA, icaC, hlgB, hlgC, hld, and aur) present in all isolates, and accessory sets, including the enterotoxin gene cluster (seo, sem, seu, sei, and sen) and a partial/complete invasive virulence factor set (splE, splA, splB, lukE, and lukD) (p = 0.001). CRSwNP isolates exhibited incomplete carriage of the core set, with frequent loss of scn, icaC, and hlgA (p < 0.05). These findings suggest that S. aureus has clusters of virulence factors that may act in concert to support the survival and persistence of the bacteria, resulting in enhanced pathogenicity. This may manifest clinically with resistant disease and refractoriness to antibiotics.

Emergence of invasive Streptococcus dysgalactiae subsp. equisimilis in Spain (2012-2022): genomic insights and clinical correlations

OBJECTIVES

An increase in Streptococcus dysgalactiae subsp. equisimilis (SDSE) infections has been documented worldwide. This study aims to analyze invasive disease caused by SDSE (iSDSE) in adults over an 11-year period in Spain.

METHODS

We conducted a retrospective, laboratory-based study of iSDSE detected at Hospital Universitari de Bellvitge (HUB) from 2012 to 2022 (n = 89) and isolates collected in three Spanish hospitals in 2018 (n = 22). Clinical data from HUB were collected. Isolates were tested for antimicrobial susceptibility (European Committee on Antimicrobial Susceptibility Testing 2023), subjected to whole genome sequencing and analyzed for mobile genetic elements (MGEs). A mouse model was used to analyze virulence.

RESULTS

iSDSE episodes at HUB occurred predominantly in older patients with comorbidities (particularly, diabetes, chronic heart disease, and malignancies). Whole genome sequencing revealed a high genetic diversity, with the most common lineages being CC15, CC17, and CC20. Various virulence factors, including the superantigen speG, were identified. Macrolides, lincosamides, and tetracyclines exhibited the highest resistance rates (>27%) and changed over time, linked to multiple MGEs. The mouse model highlighted the virulence of the CC20-stG62647 lineage, but these results were discordant with clinical data.

CONCLUSION

iSDSE incidence is increasing and associated with older patients with comorbidities. Genetically, SDSE is diverse with a high capacity to integrate MGEs carrying resistance determinants. Mouse model studies showed the enhanced virulence of the CC20-stG62647 lineage. These findings underscore the need for ongoing surveillance of this emerging pathogen.

Clinical characteristics and genomic changes of recurrent Methicillin-Resistant Staphylococcus aureus bacteremia

BACKGROUND

Recurrent or persistent methicillin-resistant Staphylococcus aureus (MRSA) bacteremia presents significant clinical challenges. Comprehensive genomic-scale studies on the genetic changes in MRSA that correspond to refractory bacteremia are lacking.

METHOD

From 2011 to 2019, MRSA blood isolates were collected from patients with persistent or recurrent bacteremia at a teaching hospital in southern Taiwan. Whole-genome sequencing (WGS) captured the genomic changes in strains responsible for refractory bacteremia, and the altered susceptibilities to specific antimicrobial agents were assessed through measurements of minimal inhibitory concentrations (MICs).

RESULT

A total of 35 MRSA blood isolates from 15 patients with recurrent or persistent bacteremia were analyzed. Reduced susceptibilities to at least one anti-MRSA agent developed in strains from seven (46.7 %) patients. Of them, a non-synonymous mutation on a global regulator mgrA was associated with reduced daptomycin susceptibility, while an increase in vancomycin MIC was linked to mutations in genes encoding LCP family protein. A 16-fold increase in MIC to fusidic acid was connected to a mutation in the elongation factor G. These recurrent strains commonly exhibited a loss or acquisition of adhesion genes that were involved in biofilm formation, including fnbA, fnbB, and sdrD, and easG series genes of type VII secretion system.

CONCLUSION

Changes in the susceptibility of successive strains to common anti-MRSA agents were frequently observed in recurrent MRSA bacteremia. These changes were linked to modifications in genes of regulatory cascade, peptidoglycan binding, adhesion, and type VII secretion system.

Comparison of pharyngeal and invasive isolates of Streptococcus pyogenes by whole-genome sequencing in Toronto, Canada

Invasive Group A streptococcal (iGAS) infections are rising in Canada and wordwide. The 2022-2023 Ontario iGAS season was among the highest recorded, a trend continuing in 2023-2024. We sequenced 38 invasive (blood) and 117 non-invasive (pharyngeal) Streptococcus pyogenes clinical isolates from Toronto (January-May 2023) to compare between the two cohorts and against published sequences to determine if any genomic changes accounted for the trend. Results demonstrated limited clustering with one small totally invasive cluster (emm49) with both invasive and non-invasive isolates represented across a diverse set of lineages. Non-invasive isolates were predominantly emm12 (70.1%), whereas invasive isolates included emm12 (26.32%), emm49 (23.68%), and emm1 (13.16%) with most emm1 strains containing the 27 SNPs that define the hypervirulent M1UK clone (58.33%). Although there were no differences in the presence of overall virulence factors/adhesin genes between cohorts, there were statistically more superantigen and DNase genes in non-invasive isolates and a rare phage gene was significantly associated with invasiveness across three emm-types. The prevalence of individual virulence factor/adhesin genes also differed between our cohorts, including a higher likelihood of speA, enn, mrp, ideS/Mac, fbaA, and fbaB in invasive isolates. There were also no significant differences across the 11 antimicrobial resistance genes identified. Finally, pharyngeal isolates had larger hydrolysis and hemolysis zones, and covS deletions were observed in only seven invasive strains. Despite there being no genetic signature that differentiated our isolates, we observed several features that were predominant in invasive strains which provides further insights into the factors that contribute to GAS invasiveness.IMPORTANCEIncreasing rates of invasive Group A streptococcal (iGAS) infections are being seen both in Canada and worldwide, which is leading to a greater disease burden caused by this pathogen. Leveraging whole-genome sequencing gives us an opportunity to better understand the underlying genetic mechanisms of streptococcal disease. By utilizing this technique, we shed light on the circulating invasive and non-invasive strains of Streptococcus pyogenes in the largest urban area in Canada from January to May 2023. GAS strains causing non-invasive disease were found to have a higher abundance of superantigen and DNase genes, whereas invasive isolates were more likely to contain M-like protein genes, the superantigen speA, the protease ideS/Mac, and/or the fibronectin-binding proteins fbaA and fbaB. This work provides valuable insights into iGAS disease which will help with surveillance, epidemiology as well as developing treatment and preventative modalities to help curb the disease burden caused by this globally important pathogen.

Influence of Sequencing Technology on Pangenome-Level Analysis and Detection of Antimicrobial Resistance Genes in ESKAPE Pathogens

As sequencing costs decrease, short-read and long-read technologies are indispensable tools for uncovering the genetic drivers behind bacterial pathogen resistance. This study explores the differences between the use of short-read (Illumina) and long-read (Oxford Nanopore Technologies [ONT]) sequencing in detecting antimicrobial resistance (AMR) genes in ESKAPE pathogens (ie, Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter cloacae). Utilizing a dataset of 1385 whole genome sequences and applying commonly used bioinformatic methods in bacterial genomics, we assessed the differences in genomic completeness, pangenome structure, and AMR gene and point mutation identification. Illumina presented higher genome completeness, while ONT identified a broader pangenome. Hybrid assembly outperformed both Illumina and ONT at identifying key AMR genetic determinants, presented results closer to Illumina's completeness, and revealed ONT-like pangenomic content. Notably, Illumina consistently detected more AMR-related point mutations than its counterparts. This highlights the importance of method selection based on research goals, particularly when using publicly available data ranging a wide timespan. Differences were also observed for specific gene classes and bacterial species, underscoring the need for a nuanced understanding of technology limitations. Overall, this study reveals the strengths and limitations of each approach, advocating for the use of Illumina for common AMR analysis, ONT for studying complex genomes and novel species, and hybrid assembly for a more comprehensive characterization, leveraging the benefits of both technologies.

Genomic investigation of MRSA bacteremia relapse reveals diverse genomic profiles but convergence in bacteremia-associated genes

Background

Recurrence of methicillin-resistant Staphylococcus aureus (MRSA) bacteremia is a high risk complication for patients. Distinguishing persistent lineages from new infections is not standardized across clinical studies.

Methods

We investigated factors contributing to recurrence of MRSA bacteremia among subjects in Philadelphia, Pennsylvania. Subject demographics and clinical history were collected and paired with whole-genome sequences of infection isolates. Recurrent bacteremia episodes were recorded and defined as relapse infections (same lineage) or new infections by genomic criteria, where a relapse contains isolates <=25 single nucleotide polymorphisms (SNP) different, and by clinical criteria. All isolates were assessed for pairwise SNP distances, common mutations, and signatures of within-host adaptation using the McDonald-Kreitman test. Clusters of transmission between relapse-associated isolates and other subject lineages were identified.

Results

Among 411 sequential subjects with MRSA bacteremia, 32 experienced recurrent bacteremia episodes, with 24 subjects having exclusively relapse infections, six with infections exclusively from a new strain, and two patients with both relapse and new infections. No concordance between a genomic and a clinical definition of relapse was evident (Cohen's Kappa = 0.18, CI: -0.41). Recurrence-associated lineages exhibited signatures of positive selection(G-test:<0.01). Genes with SNPs occurring in multiple relapse lineages have roles in antibiotic resistance and virulence, including 5 lineages with mutations in mprF and 3 lineages with mutations in rpoB, which corresponded with evolved phenotypic changes in daptomycin and rifampin resistance.

Conclusions

Recurrent infections have a diverse strain background. Relapses can be readily distinguished from newly acquired infections using genomic sequencing but not clinical criteria.

Comparative genomic analysis of emerging non-typeable Haemophilus influenzae (NTHi) causing emerging septic arthritis in Atlanta

Background

Haemophilus influenzae is a Gram-negative bacterium that can exist as a commensal organism or cause a range of diseases, from ear infections to invasive conditions like meningitis. While encapsulated H. influenzae strains have historically been linked to severe diseases, non-typeable Haemophilus influenzae (NTHi) strains, lacking an intact capsule locus, have emerged as the leading cause of invasive H. influenzae infections, particularly following the widespread use of the H. influenzae serotype b (Hib) vaccine.

Methods

In response to a significant increase in invasive NTHi infections among persons living with HIV in metropolitan Atlanta during 2017-2018, we conducted a comparative genomic analysis of two predominant NTHi clones, C1 and C2, identified during this period. These clones correspond to multilocus sequence types ST164 and ST1714, respectively. We analyzed the genomic characteristics of C1 and C2 using whole genome sequencing data and compared them to a broader pangenome of H. influenzae strains to identify potential virulence factors and genetic adaptations.

Results

Both C1 and C2 isolates were highly related within their clusters, with C1 showing a maximum of 132 SNPs and C2 showing 149 SNPs within their respective core genomes. Genomic analysis revealed significant deletions in known virulence genes, surprisingly suggesting possible attenuation of virulence. No unique accessory genes were identified that distinguished C1 and C2 from other H. influenzae strains, although both clusters exhibited a consistent loss of the pxpB gene (encoding 5-oxoprolinase subunit), replaced by a mobile cassette containing genes potentially involved in sugar metabolism. All C1 and C2 isolates showed potential enrichment in accessory genes associated with systemic infections.

Conclusions

Our study suggests that while C1 and C2 clones possess some genetic markers potentially linked to systemic infections, there are no definitive unique genetic factors that distinguish these clones as more virulent than other H. influenzae strains. The expansion of these clones in a vulnerable population may reflect both chance introduction and potential adaptations to the host environment. Further research is needed to understand the implications of these genetic findings on the clinical management and prevention of invasive NTHi infections.

RSYD-BASIC: a bioinformatic pipeline for routine sequence analysis and data processing of bacterial isolates for clinical microbiology

Background. Whole-genome sequencing of bacterial isolates is increasingly becoming routine in clinical microbiology; however, subsequent analysis often needs to be started by a bioinformatician even for comprehensive pipelines. To increase the robustness of our workflow and free up bioinformatician work hours for development and advanced analysis, we aimed to produce a robust, customizable bioinformatic pipeline for bacterial genome assembly and routine analysis results that could be initiated by non-bioinformaticians. Results. We introduce the RSYD-BASIC pipeline for bacterial isolate sequence analysis and provide a demonstration of its functionality with two datasets composed of publicly available sequences, in which comparable results are obtained in most cases. In some instances, the pipeline provided additional information, corresponding to in vitro results where these could be obtained. In routine use at our department, the pipeline has already yielded clinically relevant results, allowing us to type a variety of bacterial pathogens isolated in our clinical laboratory. We also demonstrate how RSYD-BASIC results aided in disproving a potential outbreak. Conclusion. With the RSYD-BASIC pipeline, we present a configurable reads-to-results analysis pipeline operated by non-expert users that greatly eases the investigation of potential outbreaks by expert end users. Results obtained with publicly available sequences show comparable performance to the original methods, while underlining the importance of standardized methods.

Emergence of ST3390: A non-pigmented HA-MRSA clone with enhanced virulence

Background

One of the most successful and widely-distributed hospital-associated lineages of MRSA is CC5. These strains are known from widespread antibiotic resistance, but less severe disease than CA-MRSA counterparts. Recently, CC5 descendant lineages have appeared globally with hypervirulent properties.

Methods

Herein we use genomic analyses to study the epidemiology of a rare CC5 MRSA sequence type, ST3390, circulating within Tampa General Hospital (TGH). We employ genetic tools alongside in vitro and in vivo models of virulence to study the pathogenic capabilities of strains.

Results

To date, there have only been 50 recorded instances of infection caused by ST3390 globally, with 36 of those occurring at TGH. Genomic analysis of strains identified numerous spa-types, with a t010 cluster found only at TGH. Exploration of AMR genes detected the presence of unique hybrid SCCmec types, with ~90% of TGH strains possessing components of SCCmecIa, SCCmecIIa and/or SCCmecVIII. Phenotypically, all ST3390 strains lack the staphyloxanthin pigment, which is mediated by a conserved 6 aa in frame deletion within the staphyloxanthin biosynthesis protein CrtN. ST3390 strains display high levels of cytotoxicity towards human neutrophils compared to other CC5 lineages, with several isolates displaying hypervirulence in animal models of infection.

Conclusions

This is the first study to characterize the pathogenicity and genomic architecture of the rare MRSA lineage ST3390. Our work provides a deeper understanding of the clonal expansion of CC5, and the wider diversification of S. aureus isolates within patient populations.

The association between undetected heteroresistance and antibiotic treatment failure in complicated urinary tract infection

Background

Antibiotic resistance is a worsening public health threat. One poorly understood aspect of this problem is unexpected antibiotic treatment failure; when an infecting isolate is deemed susceptible to a given antibiotic, yet treatment with that drug fails. It has been proposed that heteroresistance may be an explanation for at least some unexplained treatment failures. Heteroresistance occurs when a bacterial isolate harbors a minor subpopulation of resistant cells which coexists with a majority susceptible population. The clinical relevance of heteroresistance is not clear.

Methods

We obtained 291 index isolates from 288 unique patients in the piperacillin/tazobactam arm of the ALLIUM phase 3 clinical trial for the treatment of Gram-negative pathogens causing complicated urinary tract infections. The MIC for all isolates was below the piperacillin/tazobactam resistance breakpoint according to standard antimicrobial susceptibility testing. We performed population analysis profiles on these isolates to detect piperacillin/tazobactam heteroresistance and conducted a post hoc analysis to examine the impact of heteroresistance on clinical outcomes.

Findings

We observed that 33/288 (11.5%) of the patients were infected with isolates that were heteroresistant to piperacillin/tazobactam and that patients infected with heteroresistant isolates had an increased rate of treatment failure when compared to patients infected with a non-heteroresistant isolate (odds ratio [OR] 2.13, 95% CI 1.02, 4.41; adjusted OR 1.74, 95% CI 0.82, 3.71). Further, patients without a removable catheter were at particular risk of treatment failure from infection with heteroresistant isolates.

Interpretation

These data demonstrate that patients infected with a piperacillin/tazobactam heteroresistant isolate are at an increased risk for piperacillin/tazobactam treatment failure. The results help contextualize commonly observed unexpected antibiotic treatment failure and highlight heteroresistance as a potential cause.

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

Epidemiologic and Genomic Surveillance of Vibrio cholerae and Effectiveness of Single-Dose Oral Cholera Vaccine, Democratic Republic of the Congo

We conducted 4 years of epidemiologic and genomic surveillance of single-dose effectiveness of a killed whole-cell oral cholera vaccine (kOCV) and Vibrio cholerae transmission in the Democratic Republic of the Congo. We enrolled 1,154 patients with diarrhea; 342 of those had culture-confirmed cholera. We performed whole-genome sequencing on clinical and water V. cholerae isolates from 200 patient households, which showed annual bimodal peaks of V. cholerae clade AFR10e infections. A large clonal cholera outbreak occurred 14 months after a kOCV campaign of >1 million doses, likely because of low (9%) vaccine coverage in informal settlements. Clinical and water isolates collected in the same household were closely related, suggesting person-to-person and water-to-person transmission. Single-dose kOCV vaccine effectiveness 24 months after vaccination was 59.8% (95% CI 19.7%-79.9%), suggesting modest single-dose kOCV protection. kOCV campaigns combined with water, sanitation, and hygiene programs should be used to reduce cholera in disease-endemic settings worldwide.

Genomic epidemiology and genetic characteristics of clinical Campylobacter species cocirculating in West Bengal, India, 2019, using whole genome analysis

Campylobacter species are the most common pathogens responsible for foodborne gastroenteritis worldwide. India is a region with frequent diarrheal infections and a high level of Campylobacter infection incidence, but the detailed genomic information is limited. This study aimed to characterize 112 isolates of Campylobacter from diarrhea patients at two hospitals in Kolkata, West Bengal, by whole genome analysis. The Campylobacter isolates consisted of 90 C. jejuni, 20 C. coli, and 2 C. lari isolates. Multilocus sequence typing analysis revealed that the largest sequence type (ST) populations were ST-2131 in C. jejuni and ST-830 in C. coli and seven novel STs were found in C. jejuni and one in C. coli. Notably, ST-2131, which is rarely seen elsewhere, was positive for a sialylated LOS-related gene (wlaN +neuA + cstIII) associated with Guillain-Barré syndrome. Antibiotic resistance factors predicted from the genome sequence included blaOXA variants (58.9%), tet(O) (54.5%), tet(W) (0.9%), ant(6)-Ia (0.9%), mutation in GyrA (T86I, T86I+D90N, T86I+P104S, T86I+D90N+P104S) (79.5%), and mutation in 23S rRNA (A2075G) (12.5%). In addition to the high drug resistance of Campylobacter in Kolkata, Campylobacter pathogens were circulating that may be associated with Guillain-Barré syndrome. This study indicates the importance of genomic analysis in the surveillance of pathogens, which provides genomic information on genetic diversity, virulence mechanisms, and determinants of antimicrobial resistance.

Persistent cross-species transmission systems dominate Shiga toxin-producing Escherichia coli O157:H7 epidemiology in a high incidence region: A genomic epidemiology study

Several areas of the world suffer a notably high incidence of Shiga toxin-producing Escherichia coli. To assess the impact of persistent cross-species transmission systems on the epidemiology of E. coli O157:H7 in Alberta, Canada, we sequenced and assembled E. coli O157:H7 isolates originating from collocated cattle and human populations, 2007-2015. We constructed a timed phylogeny using BEAST2 using a structured coalescent model. We then extended the tree with human isolates through 2019 to assess the long-term disease impact of locally persistent lineages. During 2007-2015, we estimated that 88.5% of human lineages arose from cattle lineages. We identified 11 persistent lineages local to Alberta, which were associated with 38.0% (95% CI 29.3%, 47.3%) of human isolates. During the later period, six locally persistent lineages continued to be associated with human illness, including 74.7% (95% CI 68.3%, 80.3%) of reported cases in 2018 and 2019. Our study identified multiple locally evolving lineages transmitted between cattle and humans persistently associated with E. coli O157:H7 illnesses for up to 13 y. Locally persistent lineages may be a principal cause of the high incidence of E. coli O157:H7 in locations such as Alberta and provide opportunities for focused control efforts.

Circulation of a Unique Klebsiella pneumoniae Clone, ST147 NDM-1/OXA-48, in Two Diverse Hospitals in Calabria (Italy)

Background/Objectives: Carbapenem-resistant Klebsiella pneumoniae has become endemic in Europe, including in Italy, where its prevalence has risen dramatically, primarily due to epidemic clones harboring metallo-enzymes. This study aims to investigate the dissemination of K. pneumoniae strains co-producing OXA-48 and NDM-1 between two hospitals in southern Italy using molecular analyses. Methods: A total of 49 K. pneumoniae strains, predominantly co-producing OXA-48 and NDM-1, were collected between March and December 2023. Antibiotic susceptibility testing was conducted following EUCAST guidelines. Whole-genome sequencing (Illumina MiSeq) and bioinformatics tools (CARD, CLC Genomics Workbench) were used to identify resistance and virulence genes, capsule loci, and phylogenetic relationships. Results: All isolates exhibited multidrug-resistant or extensively drug-resistant profiles, including resistance to ceftazidime/avibactam and meropenem/vaborbactam. Genomic analysis revealed diverse resistance genes such as blaOXA-48, blaNDM-1, blaCTX-M-15, and blaSHV variants. Virulence genes associated with capsules, fimbriae, and siderophores were widespread. Most strains were classified as ST147 by MLST and contained various plasmids known to carry antimicrobial resistance. Phylogenetic analysis confirmed their clonal relatedness, highlighting the intra-hospital dissemination of high-risk clones. Conclusions: High-risk K. pneumoniae clones, particularly ST147, pose significant challenges in healthcare settings due to the extensive antimicrobial resistance driven by plasmid-borne resistance genes, including those that co-produce carbapenemases, like blaNDM-1 and blaOXA-48. Molecular monitoring of these clones is essential for improving targeted infection control strategies, mitigating the spread of multidrug-resistant pathogens, and managing their clinical impact effectively.

Solu: a cloud platform for real-time genomic pathogen surveillance

BACKGROUND

Genomic surveillance is extensively used for tracking public health outbreaks and healthcare-associated pathogens. Despite advancements in bioinformatics pipelines, there are still significant challenges in terms of infrastructure, expertise, and security when it comes to continuous surveillance. The existing pipelines often require the user to set up and manage their own infrastructure and are not designed for continuous surveillance that demands integration of new and regularly generated sequencing data with previous analyses. Additionally, academic projects often do not meet the privacy requirements of healthcare providers.

RESULTS

We present Solu, a cloud-based platform that integrates genomic data into a real-time, privacy-focused surveillance system.

EVALUATION

Solu's accuracy for taxonomy assignment, antimicrobial resistance genes, and phylogenetics was comparable to established pathogen surveillance pipelines. In some cases, Solu identified antimicrobial resistance genes that were previously undetected. Together, these findings demonstrate the efficacy of our platform.

CONCLUSIONS

By enabling reliable, user-friendly, and privacy-focused genomic surveillance, Solu has the potential to bridge the gap between cutting-edge research and practical, widespread application in healthcare settings. The platform is available for free academic use at https://platform.solugenomics.com .

Sequencing Analysis of Invasive Carbapenem-Resistant Klebsiella pneumoniae Isolates Secondary to Gastrointestinal Colonization

Klebsiella pneumoniae represent a common invasive infection etiological agent, whose potential carbapenem-resistance and hypermucoviscosity complicate the patient's management. Infection development often derives from gastrointestinal colonization; thus, it is fundamental to monitor asymptomatic K. pneumoniae colonization through surveillance protocols, especially for intensive care and immunocompromised patients. We described a six-month routine screening protocol from the Policlinico of Catania (Italy), while blood samples were collected from the same patients only in cases of a systemic infection suspicion. All the patients who had dissemination episodes were furtherly investigated through next-generation sequencing, analyzing both colonizing and disseminating strains. This study documents emerging invasive sequence types such as ST101, ST307, and ST395, mainly revealing blaNDM or blaKPC genes, along with siderophores and hyperproduction capsule markers as virulence factors. Most of the detected factors are presumably related to a specific plasmid content, which are extremely varied and rich. In conclusion, active surveillance through sequencing is essential to enhance awareness of local epidemiology within high-risk multi-drug resistance areas. A random sequencing analysis on the most warning microorganisms could enhance sequence typing (ST) awareness within specific settings, allowing for better prevention control strategies on their eventual persistence or diffusion.

Genetic Concordance of Staphylococcus aureus From Oropharyngeal and Sputum Cultures in People With Cystic Fibrosis

BACKGROUND

People with cystic fibrosis (CF) may not expectorate sputum at young ages or after they receive CFTR modulators. While oropharyngeal swabs are commonly used to test for lower airway pathogens, it is unknown whether Staphylococcus aureus from the oropharynx matches the strain(s) infecting the lungs. Our goal was to determine whether oropharyngeal and sputum isolates of S. aureus are genetically distinct in a cohort of patients with CF.

METHODS

We obtained historical S. aureus isolates from patients who intermittently expectorated sputum in 2018, and we prospectively cultured S. aureus from oropharyngeal swabs and sputum from subjects with CF between August 2020 and February 2022. We performed short-read whole genome sequencing, determined sequence type, and performed phylogenetic analysis using S. aureus core genome single nucleotide polymorphisms (SNPs). We assigned isolates from a patient to the same strain if they had the same sequence type and differed by ≤ 60 SNPs or the isolates were not disturbed by clade breaker analysis.

RESULTS

36 subjects had S. aureus in ≥ 1 oropharyngeal swab and ≥ 1 sputum in 2018. In the prospective collection, 31 subjects had synchronous oropharyngeal swab and sputum collections. Although polyclonal infections were detected, sputum and oropharyngeal isolates of S. aureus typically matched the same strain within study subjects, both over the span of 2018 (31/36 patients) and when collected simultaneously from 2020 to 2022 (29/31 patients).

CONCLUSIONS

In patients with CF who intermittently produce sputum, oropharyngeal swabs identify S. aureus with genetic and phenotypic similarity to those cultured from sputum.

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

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

Recent Antimicrobial Resistance Situation and Mechanisms of Resistance to Key Antimicrobials in Enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli (ETEC) is a major cause of diarrhea in developing countries and is regularly imported into developed countries as a major cause of traveler's diarrhea. ETEC is usually self-limiting and not necessarily treated with antimicrobials, although antimicrobial treatment is recommended in malnourished children, severe cases, and traveler's diarrhea. However, resistant strains to representative therapeutic agents such as ciprofloxacin and azithromycin have been reported in recent years, and multidrug-resistant ETEC has also emerged. This review discusses the recent antimicrobial resistance surveillance in ETEC and the mechanisms of resistance to major antimicrobials.

Insights into Within-Host Evolution and Dynamics of Oral and Intestinal Streptococci Unveil Niche Adaptation

The oral-gut axis is a complex system linking the oral cavity and gastrointestinal tract, impacting host health and microbial composition. This study investigates genetic changes and adaptive mechanisms employed by streptococci-one of the few genera capable of colonizing oral and intestinal niches-within the same individual. We conducted whole-genome sequencing (WGS) on 218 streptococcal isolates from saliva and fecal samples of 14 inflammatory bowel disease (IBD) patients and 12 healthy controls. Our analysis identified 16 streptococcal species, with Streptococcus infantis, S. mitis, S. parasanguinis, S. australis, and S. salivarius being the most prevalent. S. infantis dominated the oral niche in both IBD patients (33%) and healthy controls (26%). It was also the primary species in fecal samples from IBD patients and the second most prevalent in those from healthy controls. S. parasanguinis was more prevalent in the gut than in the oral cavity in both groups. Comparative genomics demonstrated a within-host microevolution of streptococci, showing adaptations via recombination and acquisition of mobile genetic elements (MGEs). Intestinal streptococcal genomes exhibited a higher proportion of intact phages and a significantly greater acquisition of the tetA gene, which confers tetracycline resistance compared to oral genomes. Core-genome single-nucleotide polymorphisms (SNPs) analysis showed significant genetic divergence between oral and intestinal streptococcal genomes within the same individual. Our findings also unveil distinct niche-specific mutation signatures within intestinal genomes, indicating the emergence of distinct clonal lineages within each niche and suggesting that within-host streptococcal evolution is individual-dependent, initiated in the oral cavity.

Effectiveness of a single dose of oral cholera vaccine: findings from epidemiological and genomic surveillance of Vibrio Cholerae in the Democratic Republic of the Congo (PICHA7 Program)

This study investigated whole-cell oral cholera vaccine (kOCV) single-dose effectiveness and transmission dynamics of Vibrio cholerae through 4 years of epidemiological and genomic surveillance in Democratic Republic of the Congo (DRC). Whole genome sequencing was performed on clinical and water V. cholerae strains from 200 patient households and found annual bimodal peaks of V. cholerae clade AFR10e. 1154 diarrhea patients were enrolled with 342 culture confirmed cholera patients. A large clonal cholera outbreak occurred 18 months after a kOCV campaign of >1 million doses of Euvichol-Plus, likely because of low vaccine coverage in informal settlements (9%). Clinical and water V. cholerae strains in the same household were more closely related than different households suggesting both person-to-person and water-to-person transmission. Single-dose kOCV vaccine effectiveness in the first 24 month after vaccination was 56.9% (95% CI: 18.6%-77.2%), suggesting a single-dose provided modest protection against medically attended cholera during the 24 months post-vaccination.

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

Objectives

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

Methods

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

Results

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

Conclusions

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

Eco-Evolutionary Drivers of Vibrio parahaemolyticus Sequence Type 3 Expansion: Retrospective Machine Learning Approach

BACKGROUND

Environmentally sensitive pathogens exhibit ecological and evolutionary responses to climate change that result in the emergence and global expansion of well-adapted variants. It is imperative to understand the mechanisms that facilitate pathogen emergence and expansion, as well as the drivers behind the mechanisms, to understand and prepare for future pandemic expansions.

OBJECTIVE

The unique, rapid, global expansion of a clonal complex of Vibrio parahaemolyticus (a marine bacterium causing gastroenteritis infections) named Vibrio parahaemolyticus sequence type 3 (VpST3) provides an opportunity to explore the eco-evolutionary drivers of pathogen expansion.

METHODS

The global expansion of VpST3 was reconstructed using VpST3 genomes, which were then classified into metrics characterizing the stages of this expansion process, indicative of the stages of emergence and establishment. We used machine learning, specifically a random forest classifier, to test a range of ecological and evolutionary drivers for their potential in predicting VpST3 expansion dynamics.

RESULTS

We identified a range of evolutionary features, including mutations in the core genome and accessory gene presence, associated with expansion dynamics. A range of random forest classifier approaches were tested to predict expansion classification metrics for each genome. The highest predictive accuracies (ranging from 0.722 to 0.967) were achieved for models using a combined eco-evolutionary approach. While population structure and the difference between introduced and established isolates could be predicted to a high accuracy, our model reported multiple false positives when predicting the success of an introduced isolate, suggesting potential limiting factors not represented in our eco-evolutionary features. Regional models produced for 2 countries reporting the most VpST3 genomes had varying success, reflecting the impacts of class imbalance.

CONCLUSIONS

These novel insights into evolutionary features and ecological conditions related to the stages of VpST3 expansion showcase the potential of machine learning models using genomic data and will contribute to the future understanding of the eco-evolutionary pathways of climate-sensitive pathogens.

Xenorhabdus bharatensis sp. nov., Xenorhabdus entomophaga sp. nov., Xenorhabdus siamensis sp. nov., and Xenorhabdus thailandensis sp. nov. Isolated from Steinernema Entomopathogenic Nematodes

Four Gram-stain-negative bacterial strains, CS20T, AUT15.5T, XENO-11T, and CCN3.3T, isolated from Steinernema entomopathogenic nematodes, were found to represent novel species within the genus Xenorhabdus (Gammaproteobacteria, Morganellaceae). In this study, we described these new species using whole-genome phylogenomic reconstructions, sequence identity values from core genome sequences, and phenotypic characterization. Phylogenetic reconstructions based on 16S rRNA gene sequences showed that: (i) strain CS20T is closely related to X. stockiae DSM 17904T, (ii) strain AUT15.5T is closely related to X. budapestensis DSM 16342T, (iii) strain XENO-11T is closely related to X. khoisanae DSM 25463T, and (iv) strain CCN3.3T is closely related to X. griffiniae DSM 17911T. The 16S rRNA gene sequence similarity value between strain CS20T and X. stockiae DSM 17904T is 97.8%, between strain AUT15.5T and X. budapestensis DSM 16342T is 98.1%, between strain XENO-11T and X. khoisanae DSM 25463T is 97.8%, and between strain CCN3.3T and X. griffiniae DSM 17911T is 98.6%. Phylogenomic reconstructions using whole-genome sequences showed that: (i) strain CS20T is closely related to X. stockiae DSM 17904T and X. innexi DSM 16336T, (ii) strain AUT15.5T is closely related to X. indica DSM 17382T, (iii) strain XENO-11T is closely related to X. khoisanae DSM 25463T, and (iv) strain CCN3.3T is closely related to X. griffiniae DSM 17911T. The digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values between strains CS20T, AUT15.5T, XENO-11T, and CCN3.3T and the type strains of their more closely related species are below the 70% and the 95-96% divergence thresholds, respectively, used for prokaryotic species delineation. Hence, we propose the following four new species: Xenorhabdus bharatensis sp. nov. (the type strain is CS20T=CCM 9320T=CCOS 2070T), X. entomophaga sp. nov. (the type strain is XENO-11T=CCM 9389T=CCOS 2111T), X. siamensis sp. nov. (the type strain is AUT15.5T=CCM 9405T=CCOS 2116T), and X. thailandensis sp. nov. (the type strain is CCN3.3T=CCM 9406T=CCOS 2115T). The following biochemical tests may be useful for differentiating the novel species from their more closely related taxa: acetoin production, arginine dihydrolase, citrate utilization, gelatinase, glucose oxidation, indole production, and tryptophan deaminase. Our study contributes to a better understanding of the biodiversity and phylogenetic relationships of entomopathogenic bacteria associated with insect parasitic nematodes.

Intraspecific Diversity of Staphylococcus aureus Populations Isolated from Cystic Fibrosis Respiratory Infections

Chronic bacterial infections are often polymicrobial, comprising multiple bacterial species or variants of the same species. Because chronic infections may last for decades, they have the potential to generate high levels of intraspecific variation through within-host diversification over time, and the potential for superinfections to occur through the introduction of multiple pathogen populations to the ongoing infection. Traditional methods for identifying infective agents generally involve isolating one single colony from a given sample, usually after selecting for a specific pathogen or antibiotic resistance profile. Isolating a recognized virulent or difficult to treat pathogen is an important part of informing clinical treatment and correlative research; however, these reductive methods alone, do not provide researchers or healthcare providers with the potentially important perspective on the true pathogen population structure and dynamics over time. To begin to address this limitation, in this study, we compare findings on Staphylococcus aureus single colonies versus and pools of colonies taken from fresh sputum samples from three patients with cystic fibrosis to isolates collected from the same sputum samples and processed by the clinical microbiology laboratory. Phenotypic and genotypic analysis of isolated S. aureus populations revealed coexisting lineages in two of three sputum samples as well as population structures that were not reflected in the single colony isolates. Altogether, our observations presented here demonstrate that clinically relevant diversity can be missed with standard sampling methods when assessing chronic infections. More broadly, this work outlines the potential impact that comprehensive population-level sampling may have for both research efforts and more effective treatment practices.

Identifying gene-level mechanisms of successful dispersal of Vibrio parahaemolyticus during El Niño events

El Niño events, the warm phase of the El Niño Southern Oscillation, facilitate the movement of warm surface waters eastwards across the Pacific Ocean. Marine organisms transported by these waters can act as biological corridors for water-borne bacteria with attachment abilities. El Niño events have been hypothesized as driving the recent emergence of Vibrio parahaemolyticus (Vp) variants, marine bacterium causing gastroenteritis, in South America, but the lack of a robust methodological framework limited any further exploration. Here, we introduce two new analysis approaches to explore Vp dynamics in South America, which will be central to uncovering Vp dynamics in the future. Distributed non-linear lag models found that strong El Niño events increase the relative probability of Vp detection in Peru, with a 3-4-month lag time. Machine learning found that the presence of a specific gene (vopZ) involved in attachment to plankton in a pandemic Vp clone in South America was temporally associated with strong El Niño events, offering a possible strategy for survival over long-range dispersal, such as that offered by El Niño events. Robust surveillance of marine pathogens and methodological development are necessary to produce resolute conclusions on the effect of El Niño events on water-borne diseases.

Photorhabdus viridis sp. nov. Isolated from Heterorhabditis zealandica Entomopathogenic Nematodes

A novel bacterial species, Photorhabdus viridis sp. nov., represented by strain GreenT, isolated from Heterorhabditis zealandica MJ2C entomopathogenic nematodes, is described. Phylogenetic reconstructions using 16S rRNA gene sequences show that strain GreenT is closely related to P. thracensis DSM 15199 T. The 16rRNA gene sequences of these two strains are 98.8% identical. Phylogenetic reconstructions using whole-genome sequences show that strain GreenT is closely related to P. tasmaniensis DSM 22387 T, P. thracensis DSM 15199 T, and P. temperata DSM 14550 T. Digital DNA-DNA hybridization (dDDH) values between strain GreenT and its three more close relative species, P. tasmaniensis DSM 22387 T, P. thracensis DSM 15199 T, and P. temperata DSM 14550 T, are 49%, 59%, and 59%, respectively. In addition, average nucleotide identity (ANI) values between GreenT and P. tasmaniensis DSM 22387 T, P. thracensis DSM 15199 T, and P. temperata DSM 14550 T are 92.4%, 94.4%, and 94.6%, respectively. The novel species also differs in their biochemical capacities from the biochemical capacities of their more closely related taxa. The following biochemical tests may be particularly useful in this context: Arginine dihydrolase, gelatinase, and glucose and mannitol oxidation. Given the clear phylogenetic separation, the sequence divergence values, and the phenotypic differences, we conclude that strain GreenT represents a novel bacterial species, for which we propose the name Photorhabdus viridis sp. nov. with GreenT (= CCM 9407 T = CCOS 2117 T = MJ2CT) as the type strain. Our study contributes to a better understanding of the taxonomy and biodiversity of an important bacterial group with great biotechnological and agricultural potential.

Genomic characterization of MRSA recovered from people with cystic fibrosis during two Spanish multicentre studies (2013 and 2021)

Background

Chronic bronchopulmonary infection due to MRSA in people with cystic fibrosis (pwCF) has been associated with accelerated decline in lung function, increased hospitalizations and increased mortality.

Material and methods

We studied microbiological and genomic characteristics of MRSA isolates recovered from pwCF in two Spanish multicentre studies (2013, 2021). Antimicrobial susceptibility was performed. WGS was carried out to determine population structure [MLST, spa-typing, staphylococcal cassette chromosome mec (SCCmec)], resistome and virulome. Clinical charts of MRSA-infected and MRSA-non-infected pwCF were also reviewed.

Results

MRSA infection prevalence decreased between 2013 (29/341, 8.5%) and 2021 (21/326, 6.4%) (P = 0.378). Differences in lung function were observed between infected and non-infected patients (P < 0.005). A higher prevalence of hospital-acquired (HA) clones was found compared with community-acquired (CA) clones (2013: 67% versus 33%; and 2021: 71% versus 29%). Overall, we noted clustering of isolates based on year of sampling, type of acquisition and clonal complex (CC). HA-MRSA population was dominated by CC5, with ST125-MRSA-IVc-t067 the most prevalent lineage (37%). A higher clonal diversity was detected among CA-MRSA. One Panton-Valentine leucocidin (PVL)-positive strain (ST8-MRSA-IV) and three strains of porcine origin (two ST398-MRSA-V-t011, one ST398-MRSA-V-t8567) were found. Additionally, acquired resistance genes (n = 24) were detected, including the cfr gene conferring linezolid resistance. A higher gentamicin resistance was found in 2021 (42%) compared with 2013 (7%) (P = 0.046), associated with the aac(6')-aph(2″) gene.

Conclusions

Despite a decrease in MRSA prevalence, we showed its potential impact on CF severity and progression. Moreover, we observed great genotypic and phenotypic diversity in MRSA isolates from pwCF as well as an MDR trait.

Molecular diversity in fusidic acid-resistant Methicillin Susceptible Staphylococcus aureus

Objectives

The recent emergence of fusidic acid (FA)-resistant Staphylococcus aureus has underscored the importance of active surveillance in isolating these strains. The molecular basis of fusidic acid resistance and the carriage of virulence factors in four borderline oxacillin-resistant Staphylococcus aureus (BORSA) clinical strains was assessed through phenotypical and genotypical methods.

Methods

All S. aureus clinical strains were obtained from various hospital units in Sicily. In vitro antibiotic susceptibility testing was conducted. WGS was performed using the Illumina MiSeq Platform, and data analysis was carried out to determine ST, resistome and virulome profiles.

Results

Genotypic characterization revealed that the strains belong to four STs: ST630, ST8, ST15, and ST1. FA resistance was associated with mutations in the fusA gene or fusB and fusC genes. Additionally, one case exhibited resistance to mupirocin, related to the presence of the mupA gene. Borderline MIC values were observed for cefoxitin in three out of four cases, leading to their categorization as BORSA. Virulence gene content was complex and diversified, with one testing positive for the lukS/F genes, coding for PVL toxin.

Conclusions

Resistance to FA is multifactorial, involving point mutations in chromosomal genes or association with mobile genetic elements. Monitoring the resistance to these antibiotics might help to manage and eradicate mupirocin- and FA-resistant S. aureus strains, which are also known to be important carriers of virulence determinants.

Host and antibiotic jointly select for greater virulence in Staphylococcus aureus

Widespread antibiotic usage has resulted in the rapid evolution of drug-resistant bacterial pathogens and poses significant threats to public health. Resolving how pathogens respond to antibiotics under different contexts is critical for understanding disease emergence and evolution going forward. The impact of antibiotics has been demonstrated most directly through in vitro pathogen passaging experiments. Independent from antibiotic selection, interactions with hosts have also altered the evolutionary trajectories and fitness landscapes of pathogens, shaping infectious disease outcomes. However, it is unclear how interactions between hosts and antibiotics impact the evolution of pathogen virulence. Here, we evolved and re-sequenced Staphylococcus aureus, a major bacterial pathogen, varying exposure to host and antibiotics to tease apart the contributions of these selective pressures on pathogen adaptation. After 12 passages, S. aureus evolving in Caenorhabditis elegans nematodes exposed to a sub-minimum inhibitory concentration of antibiotic (oxacillin) became highly virulent, regardless of whether the ancestral pathogen was methicillin-resistant (MRSA) or methicillin-sensitive (MSSA). Host and antibiotic exposure selected for reduced drug susceptibility in MSSA lineages while increasing MRSA total growth outside hosts. We identified mutations in genes involved in complex regulatory networks linking virulence and metabolism, including codY , agr , and gdpP , suggesting that rapid adaptation to infect hosts may have pleiotropic effects. In particular, MSSA populations under selection from host and antibiotic accumulated mutations in the global regulator gene codY , which controls biofilm formation in S. aureus. These populations had indeed evolved more robust biofilms-a trait linked to both virulence and antibiotic resistance-suggesting evolution of one trait can confer multiple adaptive benefits. Despite evolving in similar environments, MRSA and MSSA populations proceeded on divergent evolutionary paths, with MSSA populations exhibiting more similarities across replicate populations. Our results underscore the importance of considering multiple and concurrent selective pressures as drivers of pervasive pathogen traits.

Evolutionary dynamics of the successful expansion of pandemic Vibrio parahaemolyticus ST3 in Latin America

The underlying evolutionary mechanisms driving global expansions of pathogen strains are poorly understood. Vibrio parahaemolyticus is one of only two marine pathogens where variants have emerged in distinct climates globally. The success of a Vibrio parahaemolyticus clone (VpST3) in Latin America- the first spread identified outside its endemic region of tropical Asia- provided an invaluable opportunity to investigate mechanisms of VpST3 expansion into a distinct marine climate. A global collection of VpST3 isolates and novel Latin American isolates were used for evolutionary population genomics, pangenome analysis and combined with oceanic climate data. We found a VpST3 population (LatAm-VpST3) introduced in Latin America well before the emergence of this clone in India, previously considered the onset of the VpST3 epidemic. LatAm-VpST3 underwent successful adaptation to local conditions over its evolutionary divergence from Asian VpST3 isolates, to become dominant in Latin America. Selection signatures were found in genes providing resilience to the distinct marine climate. Core genome mutations and accessory gene presences that promoted survival over long dispersals or increased environmental fitness were associated with environmental conditions. These results provide novel insights into the global expansion of this successful V. parahaemolyticus clone into regions with different climate scenarios.

A hotspot of diversity: novel Shewanella species isolated from Baltic Sea sediments delineate a sympatric species complex

Two bacterial strains, SP1S1-4T and SP2S1-2T, were isolated from sediment samples collected in the Stockholm archipelago in November 2021. Following whole-genome sequencing, these strains were identified as tentatively belonging to two novel Shewanella genospecies, based on digital DNA-DNA hybridization, as implemented in the Type Strain Genome Server. Shewanella septentrionalis, Shewanella baltica and Shewanella hafniensis were, in this order and within a narrow genomic relatedness range, their closest genotypic relatives. Additional sampling and sequencing efforts led to the retrieval of distinct isolates that were monophyletic with SP1S1-4T and SP2S1-2T, respectively, based on phylogenomic analysis of whole-genome sequences. Comparative analyses of genome sequence data, which included blast-based average nucleotide identity, core genome-based and core proteome-based phylogenomics, in addition to MALDI-TOF MS-based protein profiling, confirmed the distinctness of the putative novel genospecies with respect to their closest genotypic relatives. A comprehensive phenotypic characterisation of SP1S1-4T and SP2S1-2T revealed only minor differences with respect to the type strains of S. septentrionalis, S. baltica and S. hafniensis. Based on the collective phylogenomic, proteomic, and phenotypic evidence presented here, we describe two novel genospecies within the genus Shewanella, for which the names Shewanella scandinavica sp. nov. and Shewanella vaxholmensis sp. nov. are proposed. The type strains are, respectively, SP2S1-2T (=CCUG 76457T=CECT 30688T), with a draft genome sequence of 5 041 805 bp and a G+C content of 46.3 mol%, and SP1S1-4T (=CCUG 76453T=CECT 30684T), with a draft genome sequence of 4 920147 bp and a G+C content of 46.0 mol%. Our findings suggest the existence of a species complex formed by the species S. baltica, S. septentrionalis, S. scandinavica sp. nov., and S. vaxholmensis sp. nov., with S. hafniensis falling in the periphery, where distinct genomic species clusters could be identified. However, this does not exclude the possibility of a continuum of genomic diversity within this sedimental ecosystem, as discussed herein with additional sequenced isolates.

Average nucleotide identity-based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital- and community-acquired infections in humans worldwide. Due to the high incidence of infection, S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole-genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of average nucleotide identity revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10%-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location. Notably, strain-diffuse genes were associated with prophages; strain-concentrated genes were associated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic resistance genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core, and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

IMPORTANCE

We analyzed the genomic diversity of Staphylococcus aureus, a globally prevalent bacterial species that causes serious infections in humans. Our goal was to build a genetic picture of the different strains of S. aureus and which genes may be associated with them. We reprocessed >84,000 genomes and subsampled to remove redundancy. We found that individual samples sharing >99.5% of their genome could be grouped into strains. We also showed that a portion of genes that are present in intermediate frequency in the species are strongly associated with some strains but completely absent from others, suggesting a role in strain specificity. This work lays the foundation for understanding individual gene histories of the S. aureus species and also outlines strategies for processing large bacterial genomic data sets.

zDB: bacterial comparative genomics made easy

The analysis and comparison of genomes rely on different tools for tasks such as annotation, orthology prediction, and phylogenetic inference. Most tools are specialized for a single task, and additional efforts are necessary to integrate and visualize the results. To fill this gap, we developed zDB, an application integrating a Nextflow analysis pipeline and a Python visualization platform built on the Django framework. The application is available on GitHub (https://github.com/metagenlab/zDB) and from the bioconda channel. Starting from annotated Genbank files, zDB identifies orthologs and infers a phylogeny for each orthogroup. A species phylogeny is also constructed from shared single-copy orthologs. The results can be enriched with Pfam protein domain prediction, Cluster of Orthologs Genes and Kyoto Encyclopedia of Genes and Genomes annotations, and Swissprot homologs. The web application allows searching for specific genes or annotations, running Blast queries, and comparing genomic regions and whole genomes. The metabolic capacities of organisms can be compared at either the module or pathway levels. Finally, users can run queries to examine the conservation of specific genes or annotations across a chosen subset of genomes and display the results as a list of genes, Venn diagram, or heatmaps. Those features make zDB useful for both bioinformaticians and researchers more accustomed to laboratory research.IMPORTANCEGenome comparison and analysis rely on many independent tools, leaving to scientists the burden to integrate and visualize their results for interpretation. To alleviate this burden, we have built zDB, a comparative genomics tool that includes both an analysis pipeline and a visualization platform. The analysis pipeline automates gene annotation, orthology prediction, and phylogenetic inference, while the visualization platform allows scientists to easily explore the results in a web browser. Among other features, the interface allows users to visually compare whole genomes and targeted regions, assess the conservation of genes or metabolic pathways, perform Blast searches, or look for specific annotations. Altogether, this tool will be useful for a broad range of applications in comparative studies between two and hundred genomes. Furthermore, it is designed to allow sharing of data sets easily at a local or international scale, thereby supporting exploratory analyses for non-bioinformaticians on the genome of their favorite organisms.

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

Dissemination of meticillin-resistant Staphylococcus aureus sequence type 8 (USA300) in Taiwan

BACKGROUND

In Taiwan, sequence type (ST) 239 and ST59 were two major clones among meticillin-resistant Staphylococcus aureus (MRSA) clinical isolates in the past two decades. USA300 (ST8) prevailed in the Americas but not in outside areas. Recently USA300 (ST8) emerged and was increasingly identified in Taiwan; we thus conducted an island-wide study to explore the role of USA300 among MRSA isolates.

METHODS

One hundred MRSA bloodstream isolates identified in 2020 from each of the six participating hospitals in Taiwan were collected and characterized. The first 10 ST8 isolates from each hospital were further analysed by whole-genome sequencing.

RESULTS

Of the 590 confirmed MRSA isolates, a total of 22 pulsotypes and 21 STs were identified. The strain of pulsotype AI/ST8 was the most common lineage identified, accounting for 187 isolates (31.7%) and dominating in five of six hospitals, followed by pulsotype A/ST239 (14.7%), pulsotype C/ST59 (13.9%) and pulsotype D/ST59 (9.2%). Of the 187 pulsotype AI/ST8 isolates, 184 isolates were characterized as USA300 and clustered in three major sub-pulsotypes, accounting for 78%. Ninety per cent of the 60 ST8 isolates for whole-genome sequencing were clustered in three major clades.

CONCLUSIONS

In 2020, USA300 became the most common clone of MRSA in Taiwan, accounting for >30% of MRSA bloodstream isolates island wide. Most of USA300 isolates circulating in Taiwan might have been imported on multiple occasions and evolved into at least three successful local clades. MRSA USA300 has successfully established its role in Taiwan, an area outside of the Americas.

Photorhabdus africana sp. nov. isolated from Heterorhabditis entomopathogenic nematodes

One Gram-negative, rod-shaped bacterial strain, isolated from an undescribed Heterorhabditis entomopathogenic nematode species was characterized to determine its taxonomic position. The 16S rRNA gene sequences indicate that it belongs to the class Gammaproteobacteria, to the family Morganellaceae, to the genus Photorhabdus, and likely represents a novel bacterial species. This strain, designated here as CRI-LCT, was therefore molecularly, biochemically, and morphologically characterized to describe the novel bacterial species. Phylogenetic reconstructions using 16S rRNA gene sequences show that CRI-LCT is closely related to P. laumondii subsp. laumondii TT01T and to P. laumondii subsp. clarkei BOJ-47T. The 16rRNA gene sequences between CRI-LCT and P. laumondii subsp. laumondii TT01T are 99.1% identical, and between CRI-LCT and P. laumondii subsp. clarkei BOJ-47T are 99.2% identical. Phylogenetic reconstructions using whole genome sequences show that CRI-LCT is closely related to P. laumondii subsp. laumondii TT01T and to P. laumondii subsp. clarkei BOJ-47T. Moreover, digital DNA-DNA hybridization (dDDH) values between CRI-LCT and its two relative species P. laumondii subsp. laumondii TT01T and P. laumondii subsp. clarkei BOJ-47T are 65% and 63%, respectively. In addition, we observed that average nucleotide identity (ANI) values between CRI-LCT and its two relative species P. laumondii subsp. laumondii TT01T and P. laumondii subsp. clarkei BOJ-47T are 95.8% and 95.5%, respectively. These values are below the 70% dDDH and the 95-96% ANI divergence thresholds that delimits prokaryotic species. Based on these genomic divergence values, and the phylogenomic separation, we conclude that CRI-LCT represents a novel bacterial species, for which we propose the name Photorhabdus africana sp. nov. with CRI-LCT (= CCM 9390T = CCOS 2112T) as the type strain. The following biochemical tests allow to differentiate P. africana sp. nov. CRI-LCT from other species of the genus, including its more closely related taxa: β-Galactosidase, citrate utilization, urease and tryptophan deaminase activities, indole and acetoin production, and glucose and inositol oxidation. Our study contributes to a better understanding of the taxonomy and biodiversity of this important bacterial group with great biotechnological and agricultural potential.

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

Strain Differences in Bloodstream and Skin Infection: Methicillin-Resistant Staphylococcus aureus Isolated in 2018-2021 in a Single Health System

Staphylococcus aureus is a common cause of skin and soft-tissue infections (SSTIs) and has become the most common cause of bloodstream infections (BSIs) in recent years, but whether the strains causing these two clinical syndromes overlap has not been studied adequately. USA300/500 (clonal complex [CC] 8-sequence type [ST] 8) and USA100 (CC5-ST5) have dominated among methicillin-resistant S aureus (MRSA) strains in the United States since the early 2000s. We compared the genomes of unselected MRSA isolates from 131 SSTIs with those from 145 BSIs at a single US center in overlapping periods in 2018-2021. CC8 MRSA was more common among SSTIs, and CC5 was more common among BSIs, consistent with prior literature. Based on clustering genomes with a threshold of 15 single-nucleotide polymorphisms, we identified clusters limited to patients with SSTI and separate clusters exclusively comprising patients with BSIs. However, we also identified eight clusters that included at least one SSTI and one BSI isolate. This suggests that virulent MRSA strains are transmitted from person to person locally in the healthcare setting or the community and that single lineages are often capable of causing both SSTIs and BSIs.

Genome characterisation and comparative analysis of Schaalia dentiphila sp. nov. and its subspecies, S. dentiphila subsp. denticola subsp. nov., from the human oral cavity

BACKGROUND

Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences.

RESULTS

Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935 T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24.

CONCLUSIONS

Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health.

SIGNIFICANCE

This research reveals two Schaalia strains, NCTC 9931 T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

AMRViz enables seamless genomics analysis and visualization of antimicrobial resistance

We have developed AMRViz, a toolkit for analyzing, visualizing, and managing bacterial genomics samples. The toolkit is bundled with the current best practice analysis pipeline allowing researchers to perform comprehensive analysis of a collection of samples directly from raw sequencing data with a single command line. The analysis results in a report showing the genome structure, genome annotations, antibiotic resistance and virulence profile for each sample. The pan-genome of all samples of the collection is analyzed to identify core- and accessory-genes. Phylogenies of the whole genome as well as all gene clusters are also generated. The toolkit provides a web-based visualization dashboard allowing researchers to interactively examine various aspects of the analysis results. Availability: AMRViz is implemented in Python and NodeJS, and is publicly available under open source MIT license at https://github.com/amromics/amrviz .

Clonal dissemination of Acinetobacter radioresistens among Humboldt penguins (Spheniscus humboldti) inhabiting a barren northern Peruvian island

Acinetobacter spp. are often isolated from natural sources, but knowledge about their presence in wild animals is fragmented and uncomplete. The present study aimed to characterize a series of Acinetobacter radioresistens isolated from Humboldt penguins (Spheniscus humboldti). Fifteen Humboldt penguins from an inhabited northern Peruvian island were sampled. Microorganisms were identified by MALDI-TOF MS. Antibiotic susceptibility to 12 antimicrobial agents was established, and clonal relationships were determined. A representative isolate was selected for whole genome sequencing (WGS). A. radioresistens were isolated from the feces of 12 (80%) Humboldt penguins, being susceptible to all the antimicrobial agents tested, except eight cefotaxime-intermediate isolates. All A. radioresistens were clonally related. WGS showed that the isolate belonged to ST1972, the presence of two chromosomal encoded carbapenemases (blaOXA-23 and a putative subclass B3 metallo-β-lactamase), and a series of point mutations in antibiotic-resistance related chromosomal genes, which were considered as polymorphisms. In addition, a few virulence factors, including a capsule-encoding operon, superoxide dismutases, catalases, phospholipases and a siderophore receptor were identified. The present results suggest that A. radioresistens may be a common member of the gut microbiota of Humboldt penguins, but further studies in other geographical areas are needed to establish this finding.

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.

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.

Draft genome sequence of Bacillus velezensis strain 3TSA-3, a potential probiotic for Pacific white shrimp Penaeus vannamei postlarvae isolated from commercial hatchery tanks

We report the draft genome of Bacillus velezensis strain 3TSA-3, isolated from Pacific white shrimp Penaeus vannamei postlarvae collected from a hatchery tank with high survival despite the presence of pathogenic Vibrio. The strain possesses genes encoding bacteriocins and lacks virulence factor genes, characteristics for a potential aquaculture probiotic.

Parallel evolution of alternate morphotypes of Chryseobacterium gleum during experimental evolution with Caenorhabditis elegans

Microbial evolution within polymicrobial communities is a complex process. Here, we report within-species diversification within multispecies microbial communities during experimental evolution with the nematode Caenorhabditis elegans. We describe morphological diversity in the target species Chryseobacterium gleum, which developed a novel colony morphotype in a small number of replicate communities. Alternate morphotypes coexisted with original morphotypes in communities, as well as in single-species experiments using evolved isolates. We found that the original and alternate morphotypes differed in motility and in spatial expansion in the presence of C. elegans. This study provides insight into the emergence and maintenance of intraspecies diversity in the context of microbial communities.