- Volume 9, Issue 12, 2023

Understanding the evolution of mobile genes is important for understanding the spread of antimicrobial resistance (AMR). Many clinically important AMR genes have been mobilized by mobile genetic elements (MGEs) on the kilobase scale, such as integrons and transposons, which can integrate into both chromosomes and plasmids and lead to rapid spread of the gene through bacterial populations. Looking at the flanking regions of these mobile genes in diverse genomes can highlight common structures and reveal patterns of MGE spread. However, historically this has been a largely descriptive process, relying on gene annotation and expert knowledge. Here we describe a general method to visualize and quantify the structural diversity around genes using pangraph to find blocks of homologous sequence. We apply this method to a set of 12 clinically important beta-lactamase genes and provide interactive visualizations of their flanking regions at https://liampshaw.github.io/flanking-regions. We show that nucleotide-level variation in the mobile gene itself generally correlates with increased structural diversity in its flanking regions, demonstrating a relationship between rates of mutational evolution and rates of structural evolution, and find a bias for greater structural diversity upstream. Our framework is a starting point to investigate general rules that apply to the horizontal spread of new genes through bacterial populations.

Lactobacillus crispatus is a member of the vaginal and gastrointestinal human microbiota. Here we determined the complete genome sequence of the probiotic strain M247 combining Nanopore and Illumina technologies. The M247 genome is organized in one circular chromosome of 2 336 109 bp, with a GC content of 37.04 % and 2303 ORFs, of which 1962 could be annotated. Analysis of the M247 mobilome, which accounts for 14 % of the whole genome, revealed the presence of: (i) Tn7088, a novel 14 105 bp long integrative and mobilizable element (IME) containing 16 ORFs; (ii) ΦM247, a novel 42 510 bp long siphovirus prophage containing 52 ORFs; (iii) three clustered regularly interspaced short palindromic repeats (CRISPRs); and (iv) 226 insertion sequences (ISs) belonging to 14 different families. Tn7088 has a modular organization including a mobilization module encoding FtsK homologous proteins and a relaxase, an integration/excision module coding for an integrase and an excisionase, and an adaptation module coding for a class I bacteriocin and homologous to the listeriolysin S (lls) locus of Listeria monocytogenes . Genome-wide homology search analysis showed the presence of Tn7088-like elements in 12 out of 23 L . crispatus complete public genomes. Mobilization and integration/excision modules are essentially conserved, while the adaptation module is variable since it is the target site for the integration of different ISs. Prophage ΦM247 contains genes for phage structural proteins, DNA replication and packaging, lysogenic and lytic cycles. ΦM247-like prophages are present in seven L . crispatus complete genomes, with sequence variability mainly due to the integration of ISs. PCR and sequencing showed that the Tn7088 IME excises from the M247 chromosome producing a circular form at a concentration of 4.32×10−5 copies per chromosome, and reconstitution of the Tn7088 chromosomal target site occurred at 6.65×10−4 copies per chromosome. The ΦM247 prophage produces an excised form and a reconstituted target site at a level of 3.90×10−5 and 2.48×10−5 copies per chromosome, respectively. This study identified two novel genetic elements in L. crispatus . Tn7088 represents the first example of an IME carrying a biosynthetic gene cluster for a class I bacteriocin in L. crispatus .

This study presents the assembly and comparative genomic analysis of luminous Photobacterium strains isolated from the light organs of 12 fish species using Oxford Nanopore Technologies (ONT) sequencing. The majority of assemblies achieved chromosome-level continuity, consisting of one large (>3 Mbp) and one small (~1.5 Mbp) contig, with near complete BUSCO scores along with varying plasmid sequences. Leveraging this dataset, this study significantly expanded the available genomes for P. leiognathi and its subspecies P. ‘mandapamensis’, enabling a comparative genomic analysis between the two lineages. An analysis of the large and small chromosomes unveiled distinct patterns of core and accessory genes, with a larger fraction of the core genes residing on the large chromosome, supporting the hypothesis of secondary chromosome evolution from megaplasmids in Vibrionaceae. In addition, we discovered a proposed new species, Photobacterium acropomis sp. nov., isolated from an acropomatid host, with an average nucleotide identify (ANI) of 93 % compared to the P. leiognathi and P. ‘mandapamensis’ strains. A comparison of the P. leiognathi and P. ‘mandapamensis’ lineages revealed minimal differences in gene content, yet highlighted the former’s larger genome size and potential for horizontal gene transfer. An investigation of the lux-rib operon, responsible for light production, indicated congruence between the presence of luxF and host family, challenging its role in differentiating P. ‘mandapamensis’ from P. leiognathi . Further insights were derived from the identification of metabolic differences, such as the presence of the NADH:quinone oxidoreductase respiratory complex I in P. leiognathi as well as variations in the type II secretion system (T2S) genes between the lineages, potentially impacting protein secretion and symbiosis. In summary, this study advances our understanding of Photobacterium genome evolution, highlighting subtle differences between closely related lineages, specifically P. leiognathi and P. ‘mandapamensis’. These findings highlight the benefit of long read sequencing for bacterial genome assembly and pangenome analysis and provide a foundation for exploring early bacterial speciation processes of these facultative light organ symbionts.

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

The ICEKp258.2 genomic island (GI) has been proposed as an important factor for the emergence and success of the globally spread carbapenem-resistant Klebsiella pneumoniae sequence type (ST) 258. However, a characterization of this horizontally acquired element is lacking. Using bioinformatic and experimental approaches, we found that ICEKp258.2 is not confined to ST258 and ST512, but also carried by ST3795 strains and emergent invasive multidrug-resistant pathogens from ST1519. We also identified several ICEKp258.2-like GIs spread among different K. pneumoniae STs, other Klebsiella species and even other pathogen genera, uncovering horizontal gene transfer events between different STs and bacterial genera. Also, the comparative and phylogenetic analyses of the ICEKp258.2-like GIs revealed that the most closely related ICEKp258.2-like GIs were harboured by ST11 strains. Importantly, we found that subinhibitory concentrations of antibiotics used in treating K. pneumoniae infections can induce the excision of this GI and modulate its gene expression. Our findings provide the basis for the study of ICEKp258.2 and its role in the success of K. pneumoniae ST258. They also highlight the potential role of antibiotics in the spread of ICEKp258.2-like GIs among bacterial pathogens.

Gram-negative bacteria use type VI secretion systems (T6SSs) to antagonize neighbouring cells. Although primarily involved in bacterial competition, the T6SS is also implicated in pathogenesis, biofilm formation and ion scavenging. Enterobacter species belong to the ESKAPE pathogens, and while their antibiotic resistance has been well studied, less is known about their pathogenesis. Here, we investigated the distribution and diversity of T6SS components in isolates of two clinically relevant Enterobacter species, E. cloacae and E. bugandensis . T6SS clusters are grouped into four types (T6SSi-T6SSiv), of which type i can be further divided into six subtypes (i1, i2, i3, i4a, i4b, i5). Analysis of a curated dataset of 31 strains demonstrated that most of them encode T6SS clusters belonging to the T6SSi type. All T6SS-positive strains possessed a conserved i3 cluster, and many harboured one or two additional i2 clusters. These clusters were less conserved, and some strains displayed evidence of deletion. We focused on a pathogenic E. bugandensis clinical isolate for comprehensive in silico effector prediction, with comparative analyses across the 31 isolates. Several new effector candidates were identified, including an evolved VgrG with a metallopeptidase domain and a Tse6-like protein. Additional effectors included an anti-eukaryotic catalase (KatN), M23 peptidase, PAAR and VgrG proteins. Our findings highlight the diversity of Enterobacter T6SSs and reveal new putative effectors that may be important for the interaction of these species with neighbouring cells and their environment.

Enterococci, especially Enterococcus faecium , are one of today’s leading causes of multidrug-resistant infections in hospital settings. The marine environment may harbour enterococci, but its role as an evolutionary niche and as a vector for the spread of enterococci is sparsely investigated. Hence, by applying enterococci in bivalves as a sentinel tool, this study aimed to describe the prevalence of enterocooci along the Norwegian coast and in addition the phylogeny of E. faecium in particular. Enterococci in batch samples of marine bivalves, harvested from 86 different locations, were quantitatively examined by a culture-dependent most probable number (MPN) method. Isolates were identified by MALDI-TOF-MS prior to antimicrobial susceptibility testing by broth microdilution. In-detail analyses of a representative selection of E. faecium isolates (n=148) were done by Illumina whole-genome sequencing, and assembled genomes were compared to closed E. faecium genomes in the public databases and to genomes from commensal and clinical isolates from Norway. Diversity among E. faecium within the same batch sample of bivalves was also explored. Enterococci were detected in 287 of the 471 examined bivalve samples, but in low concentrations with a median value of <18 MPN /100 g. From positive samples, 479 isolates of enterococci were identified belonging to ten different species, where E. faecium (n=247), Enterococcus hirae (n=114) and Enterococcus faecalis (n=66) were most frequently found. Resistance towards one or more antimicrobial agents was observed in 197 isolates (41 %), none of the isolates showed acquired resistance to vancomycin or linezolid. Phylogenetic analyses revealed high diversity among the E. faecium isolates and showed that the marine niche is dominated by strains from the non-clinical setting belonging to clade A2 (n=85) and B ( E. lactis ) (n=60). Only three isolates belonged to the hospital-associated clade A1 (ST80 and ST117). Two of these clustered with one isolate from a hospitalized patient and one from a non-hospitalized person. This study demonstrated a high prevalence, but low concentrations of enterococci in bivalves, and low levels of antimicrobial resistance. E. faecium genomes showed high population diversity and that very few E. faecium isolates in bivalves may have arisen from the human healthcare system. A systematic surveillance of target micro-organisms applying methods examining multiple isolates from the same bivalve sample provides important data to assess the enterococcal phylogeny, antimicrobial resistance and the level of faecal pollution in the marine environment.

The standalone regulator RofA is a positive regulator of the pilus locus in Streptococcus pyogenes . Found in only certain emm genotypes, RofA has been reported to regulate other virulence factors, although its role in the globally dominant emm1 S. pyogenes is unclear. Given the recent emergence of a new emm1 (M1UK) toxigenic lineage that is distinguished by three non-synonymous SNPs in rofA, we characterized the rofA regulon in six emm1 strains that are representative of the two contemporary major emm1 lineages (M1global and M1UK) using RNAseq analysis, and then determined the specific role of the M1UK-specific rofA SNPs. Deletion of rofA in three M1global strains led to altered expression of 14 genes, including six non-pilus locus genes. In M1UK strains, deletion of rofA led to altered expression of 16 genes, including nine genes that were unique to M1UK. Only the pilus locus genes were common to the RofA regulons of both lineages, while transcriptomic changes varied between strains even within the same lineage. Although introduction of the three SNPs into rofA did not impact gene expression in an M1global strain, reversal of three SNPs in an M1UK strain led to an unexpected number of transcriptomic changes that in part recapitulated transcriptomic changes seen when deleting RofA in the same strain. Computational analysis predicted that interactions with a key histidine residue in the PRD domain of RofA would differ between M1UK and M1global. RofA is a positive regulator of the pilus locus in all emm1 strains but effects on other genes are strain- and lineage-specific, with no clear, common DNA binding motif. The SNPs in rofA that characterize M1UK may impact regulation of RofA; whether they alter phosphorylation of the RofA PRD domain requires further investigation.

Between 2010 and 2015 the incidence of vancomycin-resistant Enterococcus faecium (VREfm) in Norway increased dramatically. Hence, we selected (1) a random subset of vancomycin-resistant enterococci (VRE) from the Norwegian Surveillance System for Communicable Diseases (2010–15; n=239) and (2) Norwegian vancomycin-susceptible E. faecium (VSEfm) bacteraemia isolates from the national surveillance system for antimicrobial resistance in microbes (2008 and 2014; n=261) for further analysis. Whole-genome sequences were collected for population structure, van gene cluster, mobile genetic element and virulome analysis, as well as antimicrobial susceptibility testing. Comparative genomic and phylogeographical analyses were performed with complete genomes of global E. faecium strains from the National Center for Biotechnology Information (NCBI) (1946–2022; n=272). All Norwegian VREfm and most of the VSEfm clustered with global hospital-associated sequence types (STs) in the phylogenetic subclade A1. The vanB2 subtype carried by chromosomal Tn1549 integrative conjugative elements was the dominant van type. The major Norwegian VREfm cluster types (CTs) were in accordance with concurrent European CTs. The dominant vanB-type VREfm CTs, ST192-CT3/26 and ST117-CT24, were mostly linked to a single hospital in Norway where the clones spread after independent chromosomal acquisition of Tn1549. The less prevalent vanA VRE were associated with more diverse CTs and vanA carrying Inc18 or RepA_N plasmids with toxin–antitoxin systems. Only 5 % of the Norwegian VRE were Enterococcus faecalis, all of which contained vanB. The Norwegian VREfm and VSEfm isolates harboured CT-specific virulence factor (VF) profiles supporting biofilm formation and colonization. The dominant VREfm CTs in general hosted more virulence determinants than VSEfm. The phylogenetic clade B VSEfm isolates (n=21), recently classified as Enterococcus lactis , harboured fewer VFs than E. faecium in general, and particularly subclade A1 isolates. In conclusion, the population structure of Norwegian E. faecium isolates mirrors the globally prevalent clones and particularly concurrent European VREfm/VSEfm CTs. Novel chromosomal acquisition of vanB2 on Tn1549 from the gut microbiota, however, formed a single major hospital VREfm outbreak. Dominant VREfm CTs contained more VFs than VSEfm.

Pertussis remains a public health concern in South Africa, with an increase in reported cases and outbreaks in recent years. Whole genome sequencing was performed on 32 Bordetella pertussis isolates sourced from three different surveillance programmes in South Africa between 2015 and 2019. Genome sequences were characterized using multilocus sequence typing, vaccine antigen genes (ptxP, ptxA, ptxB, prn and fimH) and overall genome structure. All isolates were sequence type 2 and harboured the pertussis toxin promoter allele ptxP3. The dominant genotype was ptxP3-ptxA1-ptxB2-prn2-fimH2 (31/32, 96.9 %), with no pertactin-deficient or other mutations in vaccine antigen genes identified. Amongst 21 isolates yielding closed genome assemblies, eight distinct genome structures were detected, with 61.9 % (13/21) of the isolates exhibiting three predominant structures. Increases in case numbers are probably not due to evolutionary changes in the genome but possibly due to other factors such as the cyclical nature of B. pertussis disease, waning immunity due to the use of acellular vaccines and/or population immunity gaps.

Thailand experienced five waves of coronavirus disease 2019 (COVID-19) between 2020 and 2022, with the Bangkok Metropolitan Region (BMR) being at the centre of all outbreaks. The molecular evolution of the causative agent of the disease, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has previously been characterized in Thailand, but a detailed spatiotemporal analysis is still lacking. In this study, we comprehensively reviewed the development and timelines of the five COVID-19 outbreaks in Thailand and the public health responses, and also conducted a phylogenetic analysis of 27 913 SARS-CoV-2 genomes from Thailand, together with 7330 global references, to investigate the virus’s spatiotemporal evolution during 2020 and 2022, with a particular focus on the BMR. Limited cross-border transmission was observed during the first four waves in 2020 and 2021, but was common in 2022, aligning well with the timeline of change in the international travel restrictions. Within the country, viruses were mostly restricted to the BMR during the first two waves in 2020, but subsequent waves in 2021 and 2022 saw extensive nationwide transmission of the virus, consistent with the timeline of relaxation of disease control measures employed within the country. Our results also suggest frequent epidemiological connections between Thailand and neighbouring countries during 2020 and 2021 despite relatively stringent international travel controls. The overall sequencing rate of the viruses circulating in the BMR was ~0.525 %, meeting the recommended benchmark, and our analysis supports that this is sufficient for tracking of the trend of the virus burden and genetic diversity. Our findings reveal insights into the local transmission dynamics of SARS-CoV-2 in Thailand, and provide a valuable reference for planning responses to future outbreaks.

In 2010 a single isolate of a trimethoprim-resistant multilocus sequence type 5, Panton–Valentine leucocidin-positive, community-associated methicillin-resistant   Staphylococcus aureus   (PVL-positive ST5 CA-MRSA), colloquially named WA121, was identified in northern Western Australia (WA). WA121 now accounts for ~14 % of all WA MRSA infections. To gain an understanding of the genetic composition and phylogenomic structure of WA121 isolates we sequenced the genomes of 155 WA121 isolates collected 2010–2021 and present a detailed genomic description. WA121 was revealed to be a single clonally expanding lineage clearly distinct from sequenced ST5 strains reported outside Australia. WA121 strains were typified by the presence of the distinct PVL phage φSa2wa-st5, the recently described methicillin resistance element SCCmecIVo carrying the trimethoprim resistance (dfrG) transposon Tn4791, the novel β-lactamase transposon Tn7702 and the epidermal cell differentiation inhibitor (EDIN-A) plasmid p2010-15611-2. We present evidence that SCCmecIVo together with Tn4791 has horizontally transferred to   Staphylococcus argenteus   and evidence of intragenomic movement of both Tn4791 and Tn7702. We experimentally demonstrate that p2010-15611-2 is capable of horizontal transfer by conjugative mobilization from one of several WA121 isolates also harbouring a pWBG749-like conjugative plasmid. In summary, WA121 is a distinct and clonally expanding Australian PVL-positive CA-MRSA lineage that is increasingly responsible for infections in indigenous communities in northern and western Australia. WA121 harbours a unique complement of mobile genetic elements and is capable of transferring antimicrobial resistance and virulence determinants to other staphylococci.

Host nutritional immunity utilizes metal deprivation to help prevent microbial infection. To investigate bacterial adaptation to such restrictive conditions, we conducted experimental evolution with two metal sequestering agents. Ethylenediaminetetraacetic acid (EDTA) and diethylenetriamine pentamethylene phosphonic acid (DTPMP) were selected as ligands because they differentially affect cellular levels of iron, manganese and zinc in Escherichia coli . Mutants of E. coli strain BW25113 were isolated after cultivation at sub-minimum inhibitory concentration (MIC) chelant levels and genetic changes potentially responsible for tolerance were identified by whole-genome sequencing. In EDTA-selected strains, mutations in the promoter region of yeiR resulted in elevated gene expression. The yeiR product, a zinc-specific metallochaperone, was confirmed to be primarily responsible for EDTA resistance. Similarly, in two of the DTPMP-selected strains, a promoter mutation increased expression of the fepA-entD operon, which encodes components of the ferric-enterobactin uptake pathway. However, in this case improved DTPMP tolerance was only detectable following overexpression of FepA or EntD in trans. Additional mutations in the cadC gene product, an acid-response regulator, preserved the neutrality of the growth medium by constitutively activating expression of the cadAB regulon. This study uncovers specific resistance mechanisms for zinc and iron starvation that could emerge by selection against host nutritional immunity or competition with heterologous metallophores. It also provides insight into the specific metals affected by these two widely used chelators critical for their antibacterial mode of action.

Bacteria cooperate by working collaboratively to defend their colonies, share nutrients, and resist antibiotics. Nevertheless, our understanding of these remarkable behaviours primarily comes from studying a few well-characterized species. Consequently, there is a significant gap in our understanding of microbial social traits, particularly in natural environments. To address this gap, we can use bioinformatic tools to identify genes that control cooperative or otherwise social traits. Existing tools address this challenge through two approaches. One approach is to identify genes that encode extracellular proteins, which can provide benefits to neighbouring cells. An alternative approach is to predict gene function using annotation tools. However, these tools have several limitations. Not all extracellular proteins are cooperative, and not all cooperative behaviours are controlled by extracellular proteins. Furthermore, existing functional annotation methods frequently miss known cooperative genes. We introduce SOCfinder as a new tool to find bacterial genes that control cooperative or otherwise social traits. SOCfinder combines information from several methods, considering if a gene is likely to [ 1 ] code for an extracellular protein [ 2 ], have a cooperative functional annotation, or [ 3 ] be part of the biosynthesis of a cooperative secondary metabolite. We use data on two extensively-studied species ( P. aeruginosa and B. subtilis ) to show that SOCfinder is better at finding known cooperative genes than existing tools. We also use theory from population genetics to identify a signature of kin selection in SOCfinder cooperative genes, which is lacking in genes identified by existing tools. SOCfinder opens up a number of exciting directions for future research, and is available to download from https://github.com/lauriebelch/SOCfinder.

In this study, we characterized 54 clinical isolates of Shigella collected in North Lebanon between 2009 and 2017 through phenotypic and genomic analyses. The most prevalent serogroup was S. sonnei, accounting for 46.3 % (25/54) of the isolates, followed by S. flexneri (27.8 %, 15/54), S. boydii (18.5 %, 10/54) and S. dysenteriae (7.4 %, 4/54). Only three isolates were pan-susceptible, and 87 % (47/54) of the isolates had multidrug resistance phenotypes. Notably, 27.8 % (15/54) of the isolates were resistant to third-generation cephalosporins (3GCs) and 77.8 % (42/54) were resistant to nalidixic acid. 3GC resistance was mediated by the extended-spectrum beta-lactamase genes bla CTX-M-15 and bla CTX-M-3, which were present on various plasmids. Quinolone resistance was conferred by single point mutations in the gyrA DNA gyrase gene, leading to GyrA S83L, GyrA D87Y or GyrA S83A amino acid substitutions. This is the first study, to our knowledge, to provide genomic insights into the serotypes of Shigella circulating in Lebanon and the various antimicrobial resistance determinants carried by these strains.