Search

Objectives: Escherichia coli can cause infections in the urinary tract and in normally sterile body sites leading to invasive E. coli disease (IED), including bacteraemia and sepsis, with older populations at increased risk. We aimed to estimate the theoretical coverage rate by the ExPEC4V and 9V vaccine candidates. In addition, we aimed at better understanding the diversity of E. coli isolates, including their genetic and phenotypic antimicrobial resistance (AMR), sequence types (STs), O-serotypes and the bacterial population structure., Methods: Blood and urine culture E. coli isolates (n = 304) were collected from hospitalized patients ≥60 years (n = 238) with IED during a multicentric, observational study across three continents. All isolates were tested for antimicrobial susceptibility, O-serotyped, whole-genome sequenced and bioinformatically analysed., Results: A large diversity of STs and of O-serotypes were identified across all centres, with O25b-ST131, O6-ST73 and O1-ST95 being the most prevalent types. A total of 45.4% and 64.7% of all isolates were found to have an O-serotype covered by the ExPEC4V and ExPEC9V vaccine candidates, respectively. The overall frequency of MDR was 37.4% and ST131 was predominant among MDR isolates. Low in-patient genetic variability was observed in cases where multiple isolates were collected from the same patient., Conclusions: Our results highlight the predominance of MDR O25b-ST131 E. coli isolates across diverse geographic areas. These findings provide further baseline data on the theoretical coverage of novel vaccines targeting E. coli associated with IED in older adults and their associated AMR levels., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Objectives: The novel dual-target triazaacenaphthylene, gepotidacin, recently showed promising results in its Phase III randomized controlled trial for the treatment of gonorrhoea. We investigated alterations in the gepotidacin GyrA and ParC targets in gonococci by in silico mining of publicly available global genomes (n = 33 213) and determined gepotidacin MICs in isolates with GyrA A92 alterations combined with other GyrA and/or ParC alterations., Methods: We examined gonococcal gyrA and parC alleles available at the European Nucleotide Archive. MICs were determined using the agar dilution method (gepotidacin) or Etest (four antimicrobials). Models of DNA gyrase and topoisomerase IV were obtained from AlphaFold and used to model gepotidacin in the binding site., Results: GyrA A92 alterations were identified in 0.24% of genomes: GyrA A92P/S/V + S91F + D95Y/A/N (0.208%), A92P + S91F (0.024%) and A92P (0.003%), but no A92T (previously associated with gepotidacin resistance) was found. ParC D86 alterations were found in 10.6% of genomes: ParC D86N/G (10.5%), D86N + S87I (0.051%), D86N + S88P (0.012%) and D86G + E91G (0.003%). One isolate had GyrA A92P + ParC D86N alterations, but remained susceptible to gepotidacin (MIC = 0.125 mg/L). No GyrA plus ParC alterations resulted in a gepotidacin MIC > 4 mg/L. Modelling of gepotidacin binding to GyrA A92/A92T/A92P suggested that gepotidacin resistance due to GyrA A92T might be linked to the formation of a new polar contact with DNA., Conclusions: In silico mining of 33 213 global gonococcal genomes (isolates from 1928 to 2023) showed that A92 is highly conserved in GyrA, while alterations in D86 of ParC are common. No GyrA plus ParC alterations caused gepotidacin resistance. MIC determination and genomic surveillance of potential antimicrobial resistance determinants are imperative., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Objectives: Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a major nosocomial infectious pathogen with rapidly increasing prevalence. The genomic epidemiological characteristics of CRKP nationwide, especially the evolving trends within the predominant clones, should be evaluated clearly., Methods: We collected 3415 K. pneumoniae strains from 28 hospitals across China. Antimicrobial susceptibility testing and WGS were performed. Subsequent genomic analyses, including sequence typing, K-locus (KL) identification, antimicrobial resistance gene screening, and virulence score assessment were performed. The phylogenetic relationship of clonal group 11 was determined based on core-genome analysis, and the presence of the pLVPK-like virulence plasmid in ST11 isolates was confirmed using plasmid core-gene analysis. Additionally, the trends of the ST11 lineage with different KL types on a global scale were investigated using Beast2., Results: Of the K. pneumoniae strains, 708 were identified as CRKP isolates (20.7%), of which 97.7% were MDR. ST11 was the predominant clone, and KPC-2 was the prevalent carbapenemase in China, although the prevalence of specific clones and carbapenemases varied by geographic region. Among ST11 isolates, KL47 and KL64 were the predominant KL types, and KL64 gradually replaced KL47, with a higher percentage of KL64 isolates harbouring the pLVPK-like plasmid. Global genome data showed a significant increase in the effective population size of KL64 over the last 5 years., Conclusions: The prevalence of CRKP was very high in certain regions in China. The increasing convergence of virulence and resistance, particularly in ST11-KL64 isolates, should be given more attention and further investigation., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Aims: Characterize global genomic features of 86 genomes of Salmonella Gallinarum (SG) and Pullorum (SP), which are important pathogens causing systemic infections in poultry., Methods and Results: All genomes harbored efflux pump encoding gene mdsA and gold tolerance genes golS and golT. Aminoglycoside (aac(6')-Ib, aadA5, aph(6)-Id, aph(3'')-Ib, ant(2'')-Ia), beta-lactam (blaTEM-1, blaTEM-135), efflux pump (mdsB), fosfomycin (fosA3), sulfonamide (sul1, sul2), tetracycline [tet(A)], trimethoprim (dfrA17), acid (asr), and disinfectant (qacEdelta1) resistance genes, gyrA, gyrB, and parC quinolone resistance point mutations, and mercury tolerance genes (mer) were found in different frequencies. Additionally, 310 virulence genes, pathogenicity islands (including SPI-1, 2, 3, 4, 5, 6, 9, 10, 12, 13, and 14), plasmids [IncFII(S), ColpVC, IncX1, IncN, IncX2, and IncC], and prophages (Fels-2, ST104, 500465-1, pro483, Gifsy-2, 103 203_sal5, Fels-1, RE-2010, vB_SenS-Ent2, and L-413C) were detected. MLST showed biovar-specific sequence types, and core genome MLST showed country-specific and global-related clusters., Conclusion: SG and SP global strains carry many virulence factors and important antimicrobial resistance genes. The diverse plasmids and prophages suggest genetic variability. MLST and cgMLST differentiated biovars and showed profiles occurring locally or worldwide., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Motivation: Mobile genetic elements (MGEs) are as ubiquitous in nature as they are varied in type, ranging from viral insertions to transposons to incorporated plasmids. Horizontal transfer of MGEs across bacterial species may also pose a significant threat to global health due to their capability to harbor antibiotic resistance genes. However, despite cheap and rapid whole-genome sequencing, the varied nature of MGEs makes it difficult to fully characterize them, and existing methods for detecting MGEs often do not agree on what should count. In this manuscript, we first define and argue in favor of a divergence-based characterization of mobile-genetic elements. Results: Using that paradigm, we present skandiver, a tool designed to efficiently detect MGEs from whole-genome assemblies without the need for gene annotation or markers. skandiver determines mobile elements via genome fragmentation, average nucleotide identity (ANI), and divergence time. By building on the scalable skani software for ANI computation, skandiver can query hundreds of complete assemblies against >65 000 representative genomes in a few minutes and 19 GB memory, providing scalable and efficient method for elucidating mobile element profiles in incomplete, uncharacterized genomic sequences. For isolated and integrated large plasmids (>10 kb), skandiver's recall was 48% and 47%, MobileElementFinder was 59% and 17%, and geNomad was 86% and 32%, respectively. For isolated large plasmids, skandiver's recall (48%) is lower than state-of-the-art reference-based methods geNomad (86%) and MobileElementFinder (59%). However, skandiver achieves higher recall on integrated plasmids and, unlike other methods, without comparing against a curated database, making skandiver suitable for discovery of novel MGEs., Availability and Implementation: https://github.com/YoukaiFromAccounting/skandiver., (© The Author(s) 2024. Published by Oxford University Press.)

Correct termination of transcription is essential for gene expression. In bacteria, factor-dependent termination relies on the Rho factor, that classically has three conserved domains. Some bacteria also have a functional insertion region. However, the variation in Rho structure among bacteria has not been analyzed in detail. This study determines the distribution, sequence conservation, and predicted features of Rho factors with diverse domain architectures by analyzing 2730 bacterial genomes. About half (49.8%) of the species analyzed have the typical Escherichia coli like Rho while most of the other species (39.8%) have diverse, atypical forms of Rho. Besides conservation of the main domains, we describe a duplicated RNA-binding domain present in specific species and novel variations in the bicyclomycin binding pocket. The additional regions observed in Rho proteins exhibit remarkable diversity. Commonly, however, they have exceptional amino acid compositions and are predicted to be intrinsically disordered, to undergo phase separation, or have prion-like behavior. Phase separation has recently been shown to play roles in Rho function and bacterial fitness during harsh conditions in one species and this study suggests a more widespread role. In conclusion, diverse atypical Rho factors are broadly distributed among bacteria, suggesting additional cellular roles., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The family Borreliaceae contains arthropod-borne spirochetes that cause two widespread human diseases, Lyme disease and relapsing fever. Lyme disease is a subacute, progressive illness with variable stage and tissue manifestations. Relapsing fever is an acute febrile illness with prominent bacteremia that may recur and disseminate, particularly to the nervous system. Clinical heterogeneity is a hallmark of both diseases. While human clinical manifestations are influenced by a wide variety of factors, including immune status and host genetic susceptibility, there is evidence that Borreliaceae microbial factors influence the clinical manifestations of human disease caused by this family of spirochetes. Despite these associations, the spirochete genes that influence the severity and manifestations of human disease are, for the most part, unknown. Recent work has identified lineage-specific expansions of lipoproteome-rich accessory genome elements in virulent clones of Borrelia burgdorferi. Using publicly available genome assemblies, it is shown that all Borreliaceae lineages for which sufficient sequence data are available harbor a similar pattern of strongly structured, lineage-specific expansions in their accessory genomes, particularly among lipoproteins, and that this pattern holds across phylogenetic scales including genera, species, and genotypes. The relationships among pangenome elements suggest that infrequent episodes of marked genomic change followed by clonal expansion in geographically and enzootically structured populations may account for the unique lineage structure of Borreliaceae. This analysis informs future genotype-phenotype studies among Borreliaceae and lays a foundation for studies of individual gene function guided by phylogenetic patterns of conservation, diversification, gain, and/or loss., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

High spontaneous mutation rate is crucial for obtaining ideal phenotype and exploring the relationship between genes and phenotype. How to break the genetic stability of organisms and increase the mutation frequency has become a research hotspot. Here, we present a practical and controllable evolutionary tool (oMut-Cgts) based on dual genetic level modification engineering for Corynebacterium glutamicum. Firstly, the modification engineering of transcription and replication levels based on RNA polymerase α subunit and DNA helicase Cgl0854 as the 'dock' of cytidine deaminase (pmCDA1) significantly increased the mutation rate, proving that the localization of pmCDA1 around transient ssDNA is necessary for genome mutation. Then, the combined modification and optimization of engineering at dual genetic level achieved 1.02 × 104-fold increased mutation rate. The genome sequencing revealed that the oMut-Cgts perform uniform and efficient C:G→T:A transitions on a genome-wide scale. Furthermore, oMut-Cgts-mediated rapid evolution of C. glutamicum with stress (acid, oxidative and ethanol) tolerance proved that the tool has powerful functions in multi-dimensional biological engineering (rapid phenotype evolution, gene function mining and protein evolution). The strategies for rapid genome evolution provided in this study are expected to be applicable to a variety of applications in all prokaryotic cells., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

Aims: To investigate enterococci carrying linezolid and vancomycin resistance genes from fecal samples recovered from wild boars., Methods and Results: Florfenicol- and vancomycin-resistant enterococci, isolated on selective agar plates, were screened by PCR for the presence of linezolid and vancomycin resistance genes. Five isolates carried optrA or poxtA linezolid resistance genes; one strain was resistant to vancomycin for the presence of vanA gene. All isolates were tested for their antibiotic susceptibility and subjected to Whole Genome Sequencing (WGS) analysis. In Enterococcus faecalis (E. faecalis) V1344 and V1676, the optrA was located on the new pV1344-optrA and pV1676-optrA plasmids, respectively, whereas in Enterococcus faecium (E. faecium) V1339 this gene was on a 22 354-bp chromosomal genetic context identical to the one detected in a human E. faecium isolate. In both E. faecium V1682 and E. durans V1343, poxtA was on the p1818-c plasmid previously found in a human E. faecium isolate. In E. faecium V1328, the vanA gene was on the Tn1546 transposon in turn located on a new pV1328-vanA plasmid. Only E. faecium V1682 successfully transferred the poxtA gene to an enterococcal recipient in filter mating assays., Conclusions: The occurrence of genetic elements carrying linezolid and vancomycin resistance genes in enterococci from wild boars is a matter of concern, moreover, the sharing of plasmids and transposons between isolates from wild animals, human, and environment indicates an exchange of genetic material between these settings., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Cytoplasmic incompatibility (CI), a non-Mendelian genetic phenomenon, involves the manipulation of host reproduction by Wolbachia, a maternally transmitted alphaproteobacterium. The underlying mechanism is centered around the CI Factor (CIF) system governed by two genes, cifA and cifB, where cifB induces embryonic lethality, and cifA counteracts it. Recent investigations have unveiled intriguing facets of this system, including diverse cifB variants, prophage association in specific strains, copy number variation, and rapid component divergence, hinting at a complex evolutionary history. We utilized comparative genomics to systematically classify CIF systems, analyze their locus structure and domain architectures, and reconstruct their diversification and evolutionary trajectories. Our new classification identifies ten distinct CIF types, featuring not just versions present in Wolbachia, but also other intracellular bacteria, and eukaryotic hosts. Significantly, our analysis of CIF loci reveals remarkable variability in gene composition and organization, encompassing an array of diverse endonucleases, variable toxin domains, deubiquitinating peptidases (DUBs), prophages, and transposons. We present compelling evidence that the components within the loci have been diversifying their sequences and domain architectures through extensive, independent lateral transfers and interlocus recombination involving gene conversion. The association with diverse transposons and prophages, coupled with selective pressures from host immunity, likely underpins the emergence of CIF loci as recombination hotspots. Our investigation also posits the origin of CifB-REase domains from mobile elements akin to CR (Crinkler-RHS-type) effectors and Tribolium Medea1 factor, which is linked to another non-Mendelian genetic phenomenon. This comprehensive genomic analysis offers novel insights into the molecular evolution and genomic foundations of Wolbachia-mediated host reproductive control., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Objectives: To explore the feasibility of linking data from enhanced surveillance patient questionnaires from each enteric fever case in England with genome sequencing data, including antimicrobial resistance (AMR) profiles, from the corresponding isolate of typhoidal salmonellae., Methods: After linking data we interrogated the merged dataset and assessed the utility of passive surveillance data to match and monitor antimicrobial treatment regimens in enteric fever patients with the AMR profiles of the infectious agent., Results: A high proportion of cases were given antibiotics (n = 1230/1415; 86.9%); half of the cases stated the class of antibiotic they were given (n = 630/1239) and half were prescribed cephalosporins (n = 316/630). Reported treatment with a combination of antibiotics increased with symptom severity. Nearly half of isolates (n = 644/1415; 45.5%) had mutations conferring resistance to ciprofloxacin. Based on genome-derived AMR profiles, typhoidal salmonellae isolates inferred to be susceptible to the recommended first-line antimicrobials were twice as likely to be isolated from individuals residing in the least deprived areas compared with the most deprived (n = 26/169; 15.4% versus n = 32/442; 7.2%)., Conclusions: Due to the high proportion of missing data obtained from patient interviews, we recommend a more transparent and systematic approach to recording the antibiotic prescription details by healthcare professionals in primary and secondary care. A more robust approach to data capture at this point in the care pathway would enable us to audit inconsistencies in the prescribing algorithms across England and ensure equitable treatment across all sections of society. Integrating prescribing data with the genome-derived AMR profiles of the causative agent at the individual patient level provides an opportunity to monitor the impact of treatment on clinical outcomes, and to promote best practice in real time., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Objectives: To characterize the mobile genetic elements and genetic localization of ileS2 in high-level mupirocin-resistant (Hi-MupR) methicillin-resistant Staphylococcus pseudintermedius (MRSP) and MRSA isolates recovered from canine and feline clinical samples., Methods: The identification of bacterial species and presence of mecA and ileS2 genes in MRSP and MRSA isolates were performed using MALDI-TOF MS and PCR, respectively. Antimicrobial resistance (AMR) phenotypes were determined by broth microdilution assays. The genome characteristics, ileS2-containing elements and staphylococcal cassette chromosome mec (SCCmec) were illustrated using complete circular genomes obtained from hybrid assembly of Illumina short-reads and Oxford Nanopore Technologies long-reads. These were analysed through phylogenetic and bioinformatics approaches., Results: A total of 18 MRSP clinical isolates and four MRSA clinical isolates exhibited the Hi-MupR phenotype and carried multiple AMR genes, including mecA and ileS2 genes. MRSP ST182-SCCmec V (n = 6) and ST282-ΨSCCmec57395-t10 (n = 4) contained the ileS2 transposable unit associated with IS257 on the chromosome. Three MRSA ST398-SCCmec V-t034/t4652 isolates carried ∼42 kb pSK41-like ileS2 plasmids, whereas similar ileS2 plasmids lacking tra genes were found in MRSP ST282-ΨSCCmec57395-t72/t21 isolates. Furthermore, a new group of ileS2 plasmids, carried by MRSP ST45-ΨSCCmec57395, ST433-ΨSCCmecKW21-t05 and ST2165-SCCmec IV-t06, and by one MRSA ST398-SCCmec V-t034 strain, shared the plasmid backbone with the cfr/fexA-carrying plasmid pM084526_1 in MRSA ST398., Conclusions: This study provides the first evidence of ileS2 integration into the S. pseudintermedius chromosome, which is a rare occurrence in staphylococcal species, and plasmids played a pivotal role in dissemination of ileS2 in both staphylococcal species., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Background: Stenotrophomonas maltophilia is a carbapenem-resistant Gram-negative pathogen increasingly responsible for difficult-to-treat nosocomial infections., Objectives: To describe the contemporary clinical characteristics and genome epidemiology of patients colonized or infected by S. maltophilia in a multicentre, prospective cohort., Methods: All patients with a clinical culture growing S. maltophilia were enrolled at six tertiary hospitals across Japan between April 2019 and March 2022. The clinical characteristics, outcomes, antimicrobial susceptibility and genomic epidemiology of cases with S. maltophilia were investigated., Results: In total, 78 patients were included representing 34 infection and 44 colonization cases. The median age was 72.5 years (IQR, 61-78), and males accounted for 53 cases (68%). The most common comorbidity was localized solid malignancy (39%). Nearly half of the patients (44%) were immunosuppressed, with antineoplastic chemotherapy accounting for 31%. The respiratory tract was the most common site of colonization (86%), whereas bacteraemia accounted for most infection cases (56%). The 30 day all-cause mortality rate was 21%, which was significantly higher in infection cases than colonization cases (35% versus 9%; adjusted HR, 3.81; 95% CI, 1.22-11.96). Susceptibility rates to ceftazidime, levofloxacin, minocycline and sulfamethoxazole/trimethoprim were 14%, 65%, 87% and 100%, respectively. The percentage of infection ranged from 13% in the unclassified group to 86% in genomic group 6A. The percentage of non-susceptibility to ceftazidime ranged from 33% in genomic group C to 100% in genomic groups 6 and 7 and genomic group geniculate., Conclusions: In this contemporary multicentre cohort, S. maltophilia primarily colonized the respiratory tract, whereas patients with bacteraemia had the highest the mortality from this pathogen. Sulfamethoxazole/trimethoprim remained consistently active, but susceptibility to levofloxacin was relatively low. The proportions of cases representing infection and susceptibility to ceftazidime differed significantly based on genomic groups., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Objectives: MDR and XDR Neisseria gonorrhoeae strains remain major public health concerns internationally, and quality-assured global gonococcal antimicrobial resistance (AMR) surveillance is imperative. The WHO global Gonococcal Antimicrobial Surveillance Programme (GASP) and WHO Enhanced GASP (EGASP), including metadata and WGS, are expanding internationally. We present the phenotypic, genetic and reference genome characteristics of the 2024 WHO gonococcal reference strains (n = 15) for quality assurance worldwide. All superseded WHO gonococcal reference strains (n = 14) were identically characterized., Material and Methods: The 2024 WHO reference strains include 11 of the 2016 WHO reference strains, which were further characterized, and four novel strains. The superseded WHO reference strains include 11 WHO reference strains previously unpublished. All strains were characterized phenotypically and genomically (single-molecule PacBio or Oxford Nanopore and Illumina sequencing)., Results: The 2024 WHO reference strains represent all available susceptible and resistant phenotypes and genotypes for antimicrobials currently and previously used (n = 22), or considered for future use (n = 3) in gonorrhoea treatment. The novel WHO strains include internationally spreading ceftriaxone resistance, ceftriaxone resistance due to new penA mutations, ceftriaxone plus high-level azithromycin resistance and azithromycin resistance due to mosaic MtrRCDE efflux pump. AMR, serogroup, prolyliminopeptidase, genetic AMR determinants, plasmid types, molecular epidemiological types and reference genome characteristics are presented for all strains., Conclusions: The 2024 WHO gonococcal reference strains are recommended for internal and external quality assurance in laboratory examinations, especially in the WHO GASP, EGASP and other GASPs, but also in phenotypic and molecular diagnostics, AMR prediction, pharmacodynamics, epidemiology, research and as complete reference genomes in WGS analysis., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

In Mycobacterium tuberculosis (MTB) control, whole genome sequencing-based molecular drug susceptibility testing (molDST-WGS) has emerged as a pivotal tool. However, the current reliance on a single-strain reference limits molDST-WGS's true potential. To address this, we introduce a new pan-lineage reference genome, 'MtbRf'. We assembled 'unmapped' reads from 3,614 MTB genomes (751 L1; 881 L2; 1,700 L3; and 282 L4) into 35 shared, annotated contigs (54 coding sequences [CDSs]). We constructed MtbRf through: (1) searching for contig homologues among genome database that precipitate results uniquely within Mycobacteria genus; (2) comparing genomes with H37Rv ('lift-over') to define 18 insertions; and (3) filling gaps in H37Rv with insertions. MtbRf adds 1.18% sequences to H37rv, salvaging >60% of previously unmapped reads. Transcriptomics confirmed gene expression of new CDSs. The new variants provided a moderate DST predictive value (AUROC 0.60-0.75). MtbRf thus unveils previously hidden genomic information and lays the foundation for lineage-specific molDST-WGS., (© The Author(s) 2024. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

Background: Carbapenemase-producing Enterobacterales (CPE) are challenging in healthcare, with resistance to multiple classes of antibiotics. This study describes the emergence of imipenemase (IMP)-encoding CPE among diverse Enterobacterales species between 2016 and 2019 across a London regional network., Methods: We performed a network analysis of patient pathways, using electronic health records, to identify contacts between IMP-encoding CPE-positive patients. Genomes of IMP-encoding CPE isolates were overlaid with patient contacts to imply potential transmission events., Results: Genomic analysis of 84 Enterobacterales isolates revealed diverse species (predominantly Klebsiella spp, Enterobacter spp, and Escherichia coli); 86% (72 of 84) harbored an IncHI2 plasmid carrying blaIMP and colistin resistance gene mcr-9 (68 of 72). Phylogenetic analysis of IncHI2 plasmids identified 3 lineages showing significant association with patient contacts and movements between 4 hospital sites and across medical specialties, which was missed in initial investigations., Conclusions: Combined, our patient network and plasmid analyses demonstrate an interspecies, plasmid-mediated outbreak of blaIMPCPE, which remained unidentified during standard investigations. With DNA sequencing and multimodal data incorporation, the outbreak investigation approach proposed here provides a framework for real-time identification of key factors causing pathogen spread. Plasmid-level outbreak analysis reveals that resistance spread may be wider than suspected, allowing more interventions to stop transmission within hospital networks.SummaryThis was an investigation, using integrated pathway networks and genomics methods, of the emergence of imipenemase-encoding carbapenemase-producing Enterobacterales among diverse Enterobacterales species between 2016 and 2019 in patients across a London regional hospital network, which was missed on routine investigations., Competing Interests: Potential conflicts of interest. A. C. M. is a majority shareholder in a NEXT Q company that is developing technology support systems for IPC. M. Gi. has received honoraria from Pfizer and Menarini Pharmaceuticals. J. F. T. holds shares in Oxford Nanopore Technologies. K. L. H. has received grant money from Shionogi and consulting fees from Cepheid. J. A. O. has received consulting fees from Gama Healthcare, Biointerations, Spectrum X, and Ondine. All other authors report no potential conflicts. Funding to pay the Open Access publication charges for this article was provided by Imperial College London. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Background: Streptococcus pneumoniae serotype 3 remains a problem globally. Malawi introduced 13-valent pneumococcal conjugate vaccine (PCV13) in 2011, but there has been no direct protection against serotype 3 carriage. We explored whether vaccine escape by serotype 3 is due to clonal expansion of a lineage with a competitive advantage., Methods: The distribution of serotype 3 Global Pneumococcal Sequence Clusters (GPSCs) and sequence types (STs) globally was assessed using sequences from the Global Pneumococcal Sequencing Project. Whole-genome sequences of 135 serotype 3 carriage isolates from Blantyre, Malawi (2015-2019) were analyzed. Comparative analysis of the capsule locus, entire genomes, antimicrobial resistance, and phylogenetic reconstructions were undertaken. Opsonophagocytosis was evaluated using serum samples from vaccinated adults and children., Results: Serotype 3 GPSC10-ST700 isolates were most prominent in Malawi. Compared with the prototypical serotype 3 capsular polysaccharide locus sequence, 6 genes are absent, with retention of capsule polysaccharide biosynthesis. This lineage is characterized by increased antimicrobial resistance and lower susceptibility to opsonophagocytic killing., Conclusions: A serotype 3 variant in Malawi has genotypic and phenotypic characteristics that could enhance vaccine escape and clonal expansion after post-PCV13 introduction. Genomic surveillance among high-burden populations is essential to improve the effectiveness of next-generation pneumococcal vaccines., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Staphylococcal cassette chromosome mec (SCCmec) typing is crucial for investigating methicillin-resistant Staphylococcus aureus (MRSA), relying primarily on the combination of ccr and mec gene complexes. To date, 19 ccr genes and 10 ccr gene complexes have been identified, forming 15 SCCmec types. With the vast release of bacterial genome sequences, mining the database for novel ccr gene complexes and SCC/SCCmec elements could enhance MRSA epidemiological studies. In this study, we identified 12 novel ccr genes (6 ccrA, 3 ccrB, and 3 ccrC) through mining of the National Center for Biotechnology Information (NCBI) database, forming 12 novel ccr gene complexes and 10 novel SCC elements. Overexpression of 5 groups of novel Ccr recombinases (CcrA9B3, CcrA10B1, CcrC3, CcrC4, and CcrC5) in a mutant MRSA strain lacking the ccr gene and extrachromosomal circular intermediate (ciSCC) production significantly promoted ciSCC production, demonstrating their biological activity. This discovery provides an opportunity to advance MRSA epidemiological research and develop database-based bacterial typing methods., Competing Interests: Potential conflicts of interests. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Endosymbiotic relationships have shaped eukaryotic life. As endosymbionts coevolve with their host, toward full integration as organelles, their genomes tend to shrink, with genes being completely lost or transferred to the host nucleus. Modern endosymbionts and organelles show diverse patterns of gene retention, and why some genes and not others are retained in these genomes is not fully understood. Recent bioinformatic study has explored hypothesized influences on these evolutionary processes, finding that hydrophobicity and amino acid chemistry predict patterns of gene retention, both in organelles across eukaryotes and in less mature endosymbiotic relationships. The exciting ongoing elucidation of endosymbiotic relationships affords an independent set of instances to test this theory. Here, we compare the properties of retained genes in the nitroplast, recently reported to be an integrated organelle, two related cyanobacterial endosymbionts that form "spheroid bodies" in their host cells, and a range of other endosymbionts, with free-living relatives of each. We find that in each case, the symbiont's genome encodes proteins with higher hydrophobicity and lower amino pKa than their free-living relative, supporting the data-derived model predicting the retention propensity of genes across endosymbiont and organelle genomes., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

This study delves into the genomic features of 10 Vibrio strains collected from deep-sea hydrothermal vents in the Pacific Ocean, providing insights into their evolutionary history and ecological adaptations. Through sequencing and pan-genome analysis involving 141 Vibrio species, we found that deep-sea strains exhibit larger genomes with unique gene distributions, suggesting adaptation to the vent environment. The phylogenomic reconstruction of the investigated isolates revealed the presence of 2 main clades: The first is monophyletic, consisting exclusively of Vibrio alginolyticus, while the second forms a monophyletic clade comprising both Vibrio antiquarius and Vibrio diabolicus species, which were previously isolated from deep-sea vents. All strains carry virulence and antibiotic resistance genes related to those found in human pathogenic Vibrio species which may play a wider ecological role other than host infection in these environments. In addition, functional genomic analysis identified genes potentially related to deep-sea survival and stress response, alongside candidate genes encoding for novel antimicrobial agents. Ultimately, the pan-genome we generated represents a valuable resource for future studies investigating the taxonomy, evolution, and ecology of Vibrio species., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Background: Bacteroides fragilis, an anaerobic gut bacterium and opportunistic pathogen, comprises two genetically divergent groups (or divisions) at the species level. Differences exist both in the core and accessory genomes and the beta-lactamase genes, with the cephalosporinase gene cepA represented only in division I and the carbapenemase gene cfiA only in division II., Methods: Multidrug resistance in a clinical B. fragilis strain was examined by whole-genome sequencing., Results: Strain CNM20200260 carried the antimicrobial resistance genes cepA, cfiA2, ant(6'), erm(F), mef(En2), est(T), tet(Q) and cat(A), along with 82-Phe mutation in gyrA (together with 47 amino acid changes in gyrA/B and parC/parE). bexA/B and other efflux pump genes were also observed. None of the detected insertion sequences was located upstream of cfiA2. The genome-based taxonomy coefficients (average nucleotide identity, DNA-DNA hybridization similarity and difference in genomic G + C%) with respect to genomes of the strains of B. fragilis division II and the novel species Bacteroides hominis (both cfiA-positive) met the criteria for CNM20200260 to belong to either species (>95%, >70% and <1%, respectively). No such similarity was seen with type strain NCTC 9343 or the representative genome FDAARGOS 1225 of B. fragilis (division I, cfiA-negative). Strain CNM20200260 harboured four out of nine Kyoto Encyclopedia of Genes and Genomes orthologues defined for division I and one of two defined for division II., Conclusions: This is the first description of the co-occurrence of cepA and cfiA in a Bacteroides strain, confirming the complexity of the taxonomy of this species., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Aeromonas dhakensis is reported as an emerging pathogenic species within the genus Aeromonas and is widely distributed in tropical coastal areas. This study provided a detailed description and characterization of a strain of A. dhakensis (202108B1) isolated from diseased Ancherythroculter nigrocauda in an inland region of China. Biochemical tests identified the isolate at the genus level, and the further molecular analysis of concatenated housekeeping gene sequences revealed that the strain belonged to the species A. dhakensis. The isolated A. dhakensis strain was resistant to five antibiotics, namely, penicillin, ampicillin, clindamycin, cephalexin, and imipenem, while it was susceptible to or showed intermediate resistance to most of the other 15 tested antibiotics. The isolated strain of A. dhakensis caused acute hemorrhagic septicemia and tissue damage in artificially infected A. nigrocauda, with a median lethal dose of 7.76 × 104 CFU/fish. The genome size of strain 202108B1 was 5 043 286 bp, including 1 chromosome and 4 plasmids. This is the first detailed report of the occurrence of infection caused by an A. dhakensis strain causing infection in an aquaculture system in inland China, providing important epidemiological data on this potential pathogenic species., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Motivation: ADP-ribosylation is a critical modification involved in regulating diverse cellular processes, including chromatin structure regulation, RNA transcription, and cell death. Bacterial ADP-ribosyltransferase toxins (bARTTs) serve as potent virulence factors that orchestrate the manipulation of host cell functions to facilitate bacterial pathogenesis. Despite their pivotal role, the bioinformatic identification of novel bARTTs poses a formidable challenge due to limited verified data and the inherent sequence diversity among bARTT members., Results: We proposed a deep learning-based model, ARTNet, specifically engineered to predict bARTTs from bacterial genomes. Initially, we introduced an effective data augmentation method to address the issue of data scarcity in training ARTNet. Subsequently, we employed a data optimization strategy by utilizing ART-related domain subsequences instead of the primary full sequences, thereby significantly enhancing the performance of ARTNet. ARTNet achieved a Matthew's correlation coefficient (MCC) of 0.9351 and an F1-score (macro) of 0.9666 on repeated independent test datasets, outperforming three other deep learning models and six traditional machine learning models in terms of time efficiency and accuracy. Furthermore, we empirically demonstrated the ability of ARTNet to predict novel bARTTs across domain superfamilies without sequence similarity. We anticipate that ARTNet will greatly facilitate the screening and identification of novel bARTTs from bacterial genomes., Availability and Implementation: ARTNet is publicly accessible at http://www.mgc.ac.cn/ARTNet/. The source code of ARTNet is freely available at https://github.com/zhengdd0422/ARTNet/., (© The Author(s) 2024. Published by Oxford University Press.)

Summary: The reliable and timely recognition of outbreaks is a key component of public health surveillance for foodborne diseases. Whole genome sequencing (WGS) offers high resolution typing of foodborne bacterial pathogens and facilitates the accurate detection of outbreaks. This detection relies on grouping WGS data into clusters at an appropriate genetic threshold. However, methods and tools for selecting and adjusting such thresholds according to the required resolution of surveillance and epidemiological context are lacking. Here we present DODGE (Dynamic Outbreak Detection for Genomic Epidemiology), an algorithm to dynamically select and compare these genetic thresholds. DODGE can analyse expanding datasets over time and clusters that are predicted to correspond to outbreaks (or "investigation clusters") can be named with established genomic nomenclature systems to facilitate integrated analysis across jurisdictions. DODGE was tested in two real-world Salmonella genomic surveillance datasets of different duration, 2 months from Australia and 9 years from the United Kingdom. In both cases only a minority of isolates were identified as investigation clusters. Two known outbreaks in the United Kingdom dataset were detected by DODGE and were recognized at an earlier timepoint than the outbreaks were reported. These findings demonstrated the potential of the DODGE approach to improve the effectiveness and timeliness of genomic surveillance for foodborne diseases and the effectiveness of the algorithm developed., Availability and Implementation: DODGE is freely available at https://github.com/LanLab/dodge and can easily be installed using Conda., (© The Author(s) 2024. Published by Oxford University Press.)

Background: The aac(6')-Im (aacA16) amikacin, netilmicin and tobramycin resistance gene cassette had been circulating globally undetected for many years in a sublineage of Acinetobacter baumannii global clone 2., Objectives: To identify sources for the aac(6')-Im fragment found in A. baumannii., Methods: MinION long-read sequencing and Unicycler hybrid assemblies were used to determine the genetic context of the aac(6')-Im gene. Quantitative reverse transcriptase PCR was used to measure expression., Results: Among >60 000 non-Acinetobacter draft genomes in the MRSN collection, the aac(6')-Im gene was detected in Pseudomonas putida and Enterobacter hormaechei isolates recovered from patients in Thailand between 2016 and 2019. Genomes of multiply resistant P. putida MRSN365855 and E. hormaechei MRSN791417 were completed. The class 1 integron containing the aac(6')-Im cassette was in the chromosome in MRSN365855, and in an HI2 plasmid in MRSN791417. However, MRSN791417 was amikacin susceptible and the gene was not expressed due to loss of the Pc promoter of the integron. Further examples of aac(6')-Im in plasmids from or the chromosome of various Gram-negative species were found in the GenBank nucleotide database. The aac(6')-Im context in integrons in pMRSN791417-8 and a Klebsiella plasmid pAMR200031 shared similarities with the aac(6')-Im region of AbGRI2-Im islands in A. baumannii. In other cases, the cassette array including the aac(6')-Im cassette was different., Conclusions: The aac(6')-Im gene is widespread, being found so far in several different species and in several different gene cassette arrays. The lack of amikacin resistance in E. hormaechei highlights the importance of correlating resistance gene content and antibiotic resistance phenotype., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Motivation: World Health Organization estimates that there were over 10 million cases of tuberculosis (TB) worldwide in 2019, resulting in over 1.4 million deaths, with a worrisome increasing trend yearly. The disease is caused by Mycobacterium tuberculosis (MTB) through airborne transmission. Treatment of TB is estimated to be 85% successful, however, this drops to 57% if MTB exhibits multiple antimicrobial resistance (AMR), for which fewer treatment options are available., Results: We develop a robust machine-learning classifier using both linear and nonlinear models (i.e. LASSO logistic regression (LR) and random forests (RF)) to predict the phenotypic resistance of Mycobacterium tuberculosis (MTB) for a broad range of antibiotic drugs. We use data from the CRyPTIC consortium to train our classifier, which consists of whole genome sequencing and antibiotic susceptibility testing (AST) phenotypic data for 13 different antibiotics. To train our model, we assemble the sequence data into genomic contigs, identify all unique 31-mers in the set of contigs, and build a feature matrix M, where M[i, j] is equal to the number of times the ith 31-mer occurs in the jth genome. Due to the size of this feature matrix (over 350 million unique 31-mers), we build and use a sparse matrix representation. Our method, which we refer to as MTB++, leverages compact data structures and iterative methods to allow for the screening of all the 31-mers in the development of both LASSO LR and RF. MTB++ is able to achieve high discrimination (F-1 >80%) for the first-line antibiotics. Moreover, MTB++ had the highest F-1 score in all but three classes and was the most comprehensive since it had an F-1 score >75% in all but four (rare) antibiotic drugs. We use our feature selection to contextualize the 31-mers that are used for the prediction of phenotypic resistance, leading to some insights about sequence similarity to genes in MEGARes. Lastly, we give an estimate of the amount of data that is needed in order to provide accurate predictions., Availability: The models and source code are publicly available on Github at https://github.com/M-Serajian/MTB-Pipeline., (© The Author(s) 2024. Published by Oxford University Press.)

Summary: Shotgun metagenomics allows for direct analysis of microbial community genetics, but scalable computational methods for the recovery of bacterial strain genomes from microbiomes remains a key challenge. We introduce Floria, a novel method designed for rapid and accurate recovery of strain haplotypes from short and long-read metagenome sequencing data, based on minimum error correction (MEC) read clustering and a strain-preserving network flow model. Floria can function as a standalone haplotyping method, outputting alleles and reads that co-occur on the same strain, as well as an end-to-end read-to-assembly pipeline (Floria-PL) for strain-level assembly. Benchmarking evaluations on synthetic metagenomes show that Floria is  > 3× faster and recovers 21% more strain content than base-level assembly methods (Strainberry) while being over an order of magnitude faster when only phasing is required. Applying Floria to a set of 109 deeply sequenced nanopore metagenomes took <20 min on average per sample and identified several species that have consistent strain heterogeneity. Applying Floria's short-read haplotyping to a longitudinal gut metagenomics dataset revealed a dynamic multi-strain Anaerostipes hadrus community with frequent strain loss and emergence events over 636 days. With Floria, accurate haplotyping of metagenomic datasets takes mere minutes on standard workstations, paving the way for extensive strain-level metagenomic analyses., Availability and Implementation: Floria is available at https://github.com/bluenote-1577/floria, and the Floria-PL pipeline is available at https://github.com/jsgounot/Floria&#95;analysis&#95;workflow along with code for reproducing the benchmarks., (© The Author(s) 2024. Published by Oxford University Press.)

Motivation: The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining., Results: We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux., Availability and Implementation: rhea is open source and available at: https://github.com/treangenlab/rhea., (© The Author(s) 2024. Published by Oxford University Press.)

We report that the Escherichia coli chromosome includes novel GC-rich genomic structural elements that trigger formation of post-replication gaps upon replisome passage. The two nearly perfect 222 bp repeats, designated Replication Risk Sequences or RRS, are each 650 kb from the terminus sequence dif and flank the Ter macrodomain. RRS sequence and positioning is highly conserved in enterobacteria. At least one RRS appears to be essential unless a 200 kb region encompassing one of them is amplified. The RRS contain a G-quadruplex on the lagging strand which impedes DNA polymerase extension producing lagging strand ssDNA gaps, $ \le$2000 bp long, upon replisome passage. Deletion of both RRS elements has substantial effects on global genome structure and topology. We hypothesize that RRS elements serve as topological relief valves during chromosome replication and segregation. There have been no screens for genomic sequences that trigger transient gap formation. Functional analogs of RRS could be widespread, possibly including some enigmatic G-quadruplexes in eukaryotes., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Many nonsporulating bacterial species survive prolonged resource exhaustion, by entering a state termed long-term stationary phase. Here, we performed long-term stationary phase evolutionary experiments on the bacterium Pseudomonas putida, followed by whole-genome sequencing of evolved clones. We show that P. putida is able to persist and adapt genetically under long-term stationary phase. We observed an accumulation of mutations within the evolving P. putida populations. Within each population, independently evolving lineages are established early on and persist throughout the 4-month-long experiment. Mutations accumulate in a highly convergent manner, with similar loci being mutated across independently evolving populations. Across populations, mutators emerge, that due to mutations within mismatch repair genes developed a much higher rate of mutation than other clones with which they coexisted within their respective populations. While these general dynamics of the adaptive process are quite similar to those we previously observed in the model bacterium Escherichia coli, the specific loci that are involved in adaptation only partially overlap between P. putida and E. coli., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Background: Staphylococcus pseudintermedius is a common opportunistic pathogen of companion dogs and an occasional human pathogen. Treatment is hampered by antimicrobial resistance including methicillin resistance encoded by mecA within the mobile genetic element SCCmec., Objectives: SCCmec elements are diverse, especially in non-Staphyloccocus aureus staphylococci, and novel variants are likely to be present in S. pseudintermedius. The aim was to characterize the SCCmec elements found in four canine clinical isolates of S. pseudintermedius., Material and Methods: Isolates were whole-genome sequenced and SCCmec elements were assembled, annotated and compared to known SCCmec types., Results and Discussion: Two novel SSCmec are present in these isolates. SCCmec7017-61515 is characterized by a novel combination of a Class A mec gene complex and a type 5 ccr previously only described in composite SCCmec elements. The other three isolates share a novel composite SCCmec with features of SCCmec types IV and VI., Conclusions: S. pseudintermedius is a reservoir of novel SSCmec elements that has implications for understanding antimicrobial resistant in veterinary and human medicine., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Motivation: Today, we know the function of only a small fraction of the protein sequences predicted from genomic data. This problem is even more salient for bacteria, which represent some of the most phylogenetically and metabolically diverse taxa on Earth. This low rate of bacterial gene annotation is compounded by the fact that most function prediction algorithms have focused on eukaryotes, and conventional annotation approaches rely on the presence of similar sequences in existing databases. However, often there are no such sequences for novel bacterial proteins. Thus, we need improved gene function prediction methods tailored for bacteria. Recently, transformer-based language models-adopted from the natural language processing field-have been used to obtain new representations of proteins, to replace amino acid sequences. These representations, referred to as protein embeddings, have shown promise for improving annotation of eukaryotes, but there have been only limited applications on bacterial genomes., Results: To predict gene functions in bacteria, we developed SAFPred, a novel synteny-aware gene function prediction tool based on protein embeddings from state-of-the-art protein language models. SAFpred also leverages the unique operon structure of bacteria through conserved synteny. SAFPred outperformed both conventional sequence-based annotation methods and state-of-the-art methods on multiple bacterial species, including for distant homolog detection, where the sequence similarity to the proteins in the training set was as low as 40%. Using SAFPred to identify gene functions across diverse enterococci, of which some species are major clinical threats, we identified 11 previously unrecognized putative novel toxins, with potential significance to human and animal health., Availability and Implementation: https://github.com/AbeelLab/safpred., (© The Author(s) 2024. Published by Oxford University Press.)

Aims: To investigate alternative resistance mechanisms among seven ceftazidime-avibactam (CZA)-resistant carbapenem-resistant Klebsiella pneumoniae (CRKP) strains lacking common antimicrobial resistance genes (ARGs) using whole genome sequencing., Methods and Results: ARG and virulence factors (VFs) were screened using the ARG database CARD and the VF database, respectively, and identified using genomic annotation data with BLAST+. Six strains were ST11 sequence types (STs), and one was ST2123. ST11 strains harbored more ARGs than the ST2123 strains. All seven strains carried multiple ARGs with efflux-mediated antibiotic resistance, including oqxA, oqxB, tet (A), qacEdltal, CRP, H-NS, Kpn-E, F, G, H, acrA, LptD, acrB, acrD, cpxA, mdtB, and mdtC. These efflux-mediated ARGs were identified in most strains and even all strains. Whole genome sequencing revealed that the ST11 strain carried multiple potential prophages, genomic islands, and integrative and conjugative elements, while the ST2123 strain carried an independent potential prophages and a genomic island., Conclusions: Whole genome sequencing analysis revealed that these seven CZA-resistant CRKP strains lacking common ARGs exhibited efflux-mediated antibiotic resistance-associated ARGs. The main mechanism by which CRKP resists CZA is antibiotic inactivation. Except for tet (A), no ARGs and validation experiments related to efflux were found. This study's results provide a new possibility for the resistance mechanism of CRKP to CZA, and we will verify this conclusion through experiments in the future., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Objectives: To determine the frequencies and clonal distributions of putative genetic determinants of resistance to antimicrobials applied for treatment of Clostridioides difficile infection (CDI), as documented in the genomic record., Methods: We scanned 26 557 C. difficile genome sequences publicly available from the EnteroBase platform for plasmids, point mutations and gene truncations previously reported to reduce susceptibility to vancomycin, fidaxomicin or metronidazole, respectively. We measured the antimicrobial susceptibility of 143 selected C. difficile isolates., Results: The frequency of mutations causing reduced susceptibility to vancomycin and metronidazole, respectively, increased strongly after 2000, peaking at up to 52% of all sequenced C. difficile genomes. However, both mutations declined sharply more recently, reflecting major changes in CDI epidemiology. We detected mutations associated with fidaxomicin resistance in several major genotypes, but found no evidence of international spread of resistant clones. The pCD-METRO plasmid, conferring metronidazole resistance, was detected in a single previously unreported C. difficile isolate, recovered from a hospital patient in Germany in 2008. The pX18-498 plasmid, putatively associated with decreased vancomycin susceptibility, was confined to related, recent isolates from the USA. Phenotype measurements confirmed that most of those genetic features were useful predictors of antibiotic susceptibility, even though ranges of MICs typically overlapped among isolates with and without specific mutations., Conclusions: Genomic data suggested that resistance to therapeutic antimicrobial drugs is rare in C. difficile. Public antimicrobial resistance marker databases were not equipped to detect most of the genetic determinants relevant to antibiotic therapy of CDI., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

In this prospective study, patients on home parenteral nutrition were twice as likely to be colonized with Staphylococcus aureus if their caregivers were carriers. Among S. aureus-positive patients and their caregivers, molecular analysis showed 68% genetically related strains. Despite decolonization, genetically related strains reappeared in 70% of patients., Competing Interests: Potential conflicts of interest. The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Bone and joint infections (BJIs) are difficult to treat and affect a growing number of patients, in which relapses are observed in 10-20% of case. These relapses, which call for prolonged antibiotic treatment and increase resistance emergence risk, may originate from ill-understood adaptation of the pathogen to the host. Here, we investigated 3 pairs of Escherichia coli strains from BJI cases and their relapses to unravel adaptations within patients. Whole-genome comparison presented evidence for positive selection and phenotypic characterization showed that biofilm formation remained unchanged, contrary to what is usually described in such cases. Although virulence was not modified, we identified the loss of 2 virulence factors contributing to immune system evasion in one of the studied strains. Other strategies, including global growth optimization and colicin production, likely allowed the strains to outcompete competitors. This work highlights the variety of strategies allowing in-patient adaptation in BJIs., Competing Interests: Potential conflicts of interest. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

We explored the utility of brief Mycobacterium tuberculosis whole-genome sequencing (WGS) "snapshots" at a sentinel site within Lima, Peru, for evaluating local transmission dynamics over time. Within a 17-km2 area, 15 of 70 (21%) isolates with WGS collected during 2011-2012 and 22 of 81 (27%) collected during 2020-2021 were clustered (P = .414), and additional isolates clustered with those from outside the area. Isolates from the later period were disproportionately related to large historic clusters in Lima from the earlier period. WGS snapshots at a sentinel site may not be useful for monitoring transmission, but monitoring the persistence of large transmission clusters might be., Competing Interests: Potential conflicts of interest. The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2023. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Porphyromonas gingivalis is a nonmotile, obligate anaerobic, Gram-negative bacterium known for its association with periodontal disease and its involvement in systemic diseases such as atherosclerosis, cardiovascular disease, colon cancer, and Alzheimer's disease. This bacterium produces several virulence factors, including capsules, fimbriae, lipopolysaccharides, proteolytic enzymes, and hemagglutinins. A comparative genomic analysis revealed the open pangenome of P. gingivalis and identified complete type IV secretion systems in strain KCOM2805 and almost complete type VI secretion systems in strains KCOM2798 and ATCC49417, which is a new discovery as previous studies did not find the proteins involved in secretion systems IV and VI. Conservation of some virulence factors between different strains was observed, regardless of their genetic diversity and origin. In addition, we performed for the first time a reconstruction analysis of the gene regulatory network, identifying transcription factors and proteins involved in the regulatory mechanisms of bacterial pathogenesis. In particular, QseB regulates the expression of hemagglutinin and arginine deaminase, while Rex may suppress the release of gingipain through interactions with PorV and the formatum/nitrate transporter. Our study highlights the central role of conserved virulence factors and regulatory pathways, particularly QseB and Rex, in P. gingivalis and provides insights into potential therapeutic targets., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Background: The Acinetobacter baumannii isolate called SMAL, previously used to determine the structures of capsular polysaccharide and lipooligosaccharide, was recovered in Pavia, Italy in 2002 among the collection of aminoglycoside-resistant isolates designated as SMAL type. This type was later called the Italian clone, then ST78. ST78 isolates are now widely distributed., Objectives: To establish the resistance gene complement and the location and structure of acquired resistance regions in early members of the Italian/ST78 clone., Methods: The draft genome of SMAL2002 was assembled from Illumina MiSeq reads. Contigs containing resistance genes were joined and located in the chromosome using PCR with custom primers. The resistance profile was determined using disc diffusion., Results: SMAL2002 is an ST78A isolate and includes three aminoglycoside resistance genes, aadB (gentamicin, kanamycin, tobramycin) aphA1 (kanamycin, neomycin) and aac(6')-Ian (amikacin, kanamycin, tobramycin). The aadB gene cassette is incorporated at a secondary site in a relative of the aphA1-containing, IS26-bounded pseudo-compound transposon, PTn6020. The aac(6')-Ian gene is in an adjacent IS26-bounded structure that includes sul2 (sulphonamide) and floR (florfenicol) resistance genes. The two pseudo-compound transposons overlap and are in the chromosomal hutU gene flanked by an 8 bp target site duplication. Although aac(6')-Ian was not noticed previously, the same genes and structures were found in several available draft genomes of early ST78A isolates., Conclusions: This study highlights the importance of correlating resistance profiles with resistance gene content. The location of acquired resistance genes in the SMAL2002 chromosome represents the original location in the ST78A lineage of ST78., (© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.)

Motivation: Since 2016, the number of microbial species with available reference genomes in NCBI has more than tripled. Multiple genome alignment, the process of identifying nucleotides across multiple genomes which share a common ancestor, is used as the input to numerous downstream comparative analysis methods. Parsnp is one of the few multiple genome alignment methods able to scale to the current era of genomic data; however, there has been no major release since its initial release in 2014., Results: To address this gap, we developed Parsnp v2, which significantly improves on its original release. Parsnp v2 provides users with more control over executions of the program, allowing Parsnp to be better tailored for different use-cases. We introduce a partitioning option to Parsnp, which allows the input to be broken up into multiple parallel alignment processes which are then combined into a final alignment. The partitioning option can reduce memory usage by over 4× and reduce runtime by over 2×, all while maintaining a precise core-genome alignment. The partitioning workflow is also less susceptible to complications caused by assembly artifacts and minor variation, as alignment anchors only need to be conserved within their partition and not across the entire input set. We highlight the performance on datasets involving thousands of bacterial and viral genomes., Availability and Implementation: Parsnp v2 is available at https://github.com/marbl/parsnp., (© The Author(s) 2024. Published by Oxford University Press.)

Gene duplication contributes to the evolution of expression and the origin of new genes, but the relative importance of different patterns of duplicate gene expression and mechanisms of retention remains debated and particularly poorly understood in bacteria. Here, we investigated gene expression patterns for two lab strains of the cyanobacterium Acaryochloris marina with expanding genomes that contain about 10-fold more gene duplicates compared with most bacteria. Strikingly, we observed a generally stoichiometric pattern of greater combined duplicate transcript dosage with increased gene copy number, in contrast to the prevalence of expression reduction reported for many eukaryotes. We conclude that increased transcript dosage is likely an important mechanism of initial duplicate retention in these bacteria and may persist over long periods of evolutionary time. However, we also observed that paralog expression can diverge rapidly, including possible functional partitioning, for which different copies were respectively more highly expressed in at least one condition. Divergence may be promoted by the physical separation of most Acaryochloris duplicates on different genetic elements. In addition, expression pattern for ancestrally shared duplicates could differ between strains, emphasizing that duplicate expression fate need not be deterministic. We further observed evidence for context-dependent transcript dosage, where the aggregate expression of duplicates was either greater or lower than their single-copy homolog depending on physiological state. Finally, we illustrate how these different expression patterns of duplicated genes impact Acaryochloris biology for the innovation of a novel light-harvesting apparatus and for the regulation of recA paralogs in response to environmental change., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Background: Livestock-associated MRSA (LA-MRSA) transmission/cross-contamination can occur at abattoir through colonized pigs, increasing occupational hazards and health concerns for workers. To assess this risk we used genomics to identify LA-MRSA lineages present in batches of pigs sent to slaughter and distribution of clones., Methods: WGS was performed on 85 LA-MRSA previously isolated from six abattoirs from 105 batches of pigs sent from 100 UK farms. spa typing and MLST were performed on all isolates. A mashtree tree was constructed to compare genomes of the LA-MRSA with 1281 global isolates from livestock and humans. A phylogenetic tree and pairwise SNP distance matrices were built from whole genomes of 109 isolates closest to those from abattoirs to compare evolutionary relationships and identify clones., Results: All abattoir isolates belonged to CC398 and were mainly of spa type t011, although other spa types were present. Phylogenetic analysis confirmed the abattoir isolates were most closely related to each other and to pig LA-MRSA from across Europe, indicating a common evolutionary origin with related lineages colonizing UK pigs.Comparison of genomes using SNPs suggested between one and four clones were transferring between pigs from different batches. Transmission likely occurred on farm premises, during transportation, and/or within abattoirs through contact with contaminated surfaces in lairage or post-stunning., Conclusions: Genomics forensically identified related isolates/clones circulating in pigs at slaughter, showing contamination occurs often. Results suggest that further genomic tracking will identify hotspots, and improvements in measures such as biosecurity and disinfection will help reduce risk for workers., (© Crown copyright 2024.)