Investigating the nature of prokaryotic genomic island locations within a genome.

Horizontal gene transfer (HGT) is a powerful evolutionary force that considerably shapes the structure of prokaryotic genomes and is associated with genomic islands (GIs). A GI is a DNA segment composed of transferred genes that can be found within a prokaryotic genome, obtained through HGT. Much research has focused on detecting GIs in genomes, but here we pursue a new course, which is identifying possible preferred locations of GIs in the prokaryotic genome. Here, we identify the locations of the GIs within prokaryotic genomes to examine patterns in those locations. Prokaryotic GIs were analyzed according to the genome structure that they are located in, whether it be a circular or a linear genome. The analytical investigations employed are: (1) studying the GI locations in relation to the origin of replication (oriC); (2) exploring the distances between GIs; and (3) determining the distribution of GIs across the genomes. For each of the investigations, the analysis was performed on all of the GIs in the data set. Moreover, to void bias caused by the distribution of the genomes represented, the GIs in one genome from each species and the GIs of the most frequent species are also analyzed. Overall, the results showed that there are preferred sites for the GIs in the genome. In the linear genomes, these sites are usually located in the oriC region and terminus region, while in the circular genomes, they are located solely in the terminus region. These results also showed that the distance distribution between the GIs is almost exponential, which proves that GIs have preferred sites within genomes. The oriC and termniuns are preferred sites for the GIs and a possible natural explanation for this could be connected to the content of the oriC region. Moreover, the content of the GIs in terms of its protein families was studied and the results demonstrated that the majority of frequent protein families are close to identical in each section.

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

BACKGROUND: Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences. RESULTS: Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935 T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24. CONCLUSIONS: Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health. SIGNIFICANCE: This research reveals two Schaalia strains, NCTC 9931 T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

Genomic insights of Salmonella isolated from dry fermented sausage production chains in Spain and France.

The presence of Salmonella in dry fermented sausages is source of recalls and outbreaks. The genomic diversity of 173 Salmonella isolates from the dry fermented sausage production chains (pig carcasses, pork, and sausages) from France and Spain were investigated through their core phylogenomic relationships and accessory genome profiles. Ten different serovars and thirteen sequence type profiles were identified. The most frequent serovar from sausages was the monophasic variant of S. Typhimurium (1,4,[5],12:i:-, 72%) while S. Derby was in pig carcasses (51%). Phylogenomic clusters found in S. 1,4,[5],12:i:-, S. Derby, S. Rissen and S. Typhimurium serovars identified closely related isolates, with less than 10 alleles and 20 SNPs of difference, displaying Salmonella persistence along the pork production chain. Most of the S. 1,4,[5],12:i:- contained the Salmonella genomic island-4 (SGI-4), Tn21 and IncFIB plasmid. More than half of S. Derby strains contained the SGI-1 and Tn7. S. 1,4,[5],12:i:- genomes carried the most multidrug resistance genes (91% of the strains), whereas extended-spectrum ß-lactamase genes were found in Typhimurium and Derby serovars. Salmonella monitoring and characterization in the pork production chains, specially S. 1,4,[5],12:i:- serovar, is of special importance due to its multidrug resistance capacity and persistence in dry fermented sausages.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related.

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

The evolution of antibiotic resistance islands occurs within the framework of plasmid lineages.

Bacterial pathogens carrying multidrug resistance (MDR) plasmids are a major threat to human health. The acquisition of antibiotic resistance genes (ARGs) in plasmids is often facilitated by mobile genetic elements that copy or translocate ARGs between DNA molecules. The agglomeration of mobile elements in plasmids generates resistance islands comprising multiple ARGs. However, whether the emergence of resistance islands is restricted to specific MDR plasmid lineages remains understudied. Here we show that the agglomeration of ARGs in resistance islands is biased towards specific large plasmid lineages. Analyzing 6784 plasmids in 2441 Escherichia, Salmonella, and Klebsiella isolates, we quantify that 84% of the ARGs in MDR plasmids are found in resistance islands. We furthermore observe rapid evolution of ARG combinations in resistance islands. Most regions identified as resistance islands are shared among closely related plasmids but rarely among distantly related plasmids. Our results suggest the presence of barriers for the dissemination of ARGs between plasmid lineages, which are related to plasmid genetic properties, host range and the plasmid evolutionary history. The agglomeration of ARGs in plasmids is attributed to the workings of mobile genetic elements that operate within the framework of existing plasmid lineages.

Structure classification of the proteins from Salmonella enterica pangenome revealed novel potential pathogenicity islands.

Salmonella enterica is a pathogenic bacterium known for causing severe typhoid fever in humans, making it important to study due to its potential health risks and significant impact on public health. This study provides evolutionary classification of proteins from Salmonella enterica pangenome. We classified 17,238 domains from 13,147 proteins from 79,758 Salmonella enterica strains and studied in detail domains of 272 proteins from 14 characterized Salmonella pathogenicity islands (SPIs). Among SPIs-related proteins, 90 proteins function in the secretion machinery. 41% domains of SPI proteins have no previous sequence annotation. By comparing clinical and environmental isolates, we identified 3682 proteins that are overrepresented in clinical group that we consider as potentially pathogenic. Among domains of potentially pathogenic proteins only 50% domains were annotated by sequence methods previously. Moreover, 36% (1330 out of 3682) of potentially pathogenic proteins cannot be classified into Evolutionary Classification of Protein Domains database (ECOD). Among classified domains of potentially pathogenic proteins the most populated homology groups include helix-turn-helix (HTH), Immunoglobulin-related, and P-loop domains-related. Functional analysis revealed overrepresentation of these protein in biological processes related to viral entry into host cell, antibiotic biosynthesis, DNA metabolism and conformation change, and underrepresentation in translational processes. Analysis of the potentially pathogenic proteins indicates that they form 119 clusters or novel potential pathogenicity islands (NPPIs) within the Salmonella genome, suggesting their potential contribution to the bacterium's virulence. One of the NPPIs revealed significant overrepresentation of potentially pathogenic proteins. Overall, our analysis revealed that identified potentially pathogenic proteins are poorly studied.

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

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

Investigation of genomic and pathogenicity characteristics of Streptococcus suis ST1 human strains from Guangxi Zhuang Autonomous Region (GX) between 2005 and 2020 in China.

Streptococcus suis is a significant and emerging zoonotic pathogen. ST1 and ST7 strains are the primary agents responsible for S. suis human infections in China, including the Guangxi Zhuang Autonomous Region (GX). To enhance our understanding of S. suis ST1 population characteristics, we conducted an investigation into the phylogenetic structure, genomic features, and virulence levels of 73 S. suis ST1 human strains from GX between 2005 and 2020. The ST1 GX strains were categorized into three lineages in phylogenetic analysis. Sub-lineage 3-1a exhibited a closer phylogenetic relationship with the ST7 epidemic strain SC84. The strains from lineage 3 predominantly harboured 89K-like pathogenicity islands (PAIs) which were categorized into four clades based on sequence alignment. The acquirement of 89K-like PAIs increased the antibiotic resistance and pathogenicity of corresponding transconjugants. We observed significant diversity in virulence levels among the 37 representative ST1 GX strains, that were classified as follows: epidemic (E)/highly virulent (HV) (32.4%, 12/37), virulent plus (V+) (29.7%, 11/37), virulent (V) (18.9%, 7/37), and lowly virulent (LV) (18.9%, 7/37) strains based on survival curves and mortality rates at different time points in C57BL/6 mice following infection. The E/HV strains were characterized by the overproduction of tumour necrosis factor (TNF)-α in serum and promptly established infection at the early phase of infection. Our research offers novel insights into the population structure, evolution, genomic features, and pathogenicity of ST1 strains. Our data also indicates the importance of establishing a scheme for characterizing and subtyping the virulence levels of S. suis strains.

Dual proteomics of infected macrophages reveal bacterial and host players involved in the Francisella intracellular life cycle and cell to cell dissemination by merocytophagy.

Bacterial pathogens adapt and replicate within host cells, while host cells develop mechanisms to eliminate them. Using a dual proteomic approach, we characterized the intra-macrophage proteome of the facultative intracellular pathogen, Francisella novicida. More than 900 Francisella proteins were identified in infected macrophages after a 10-h infection. Biotin biosynthesis-related proteins were upregulated, emphasizing the role of biotin-associated genes in Francisella replication. Conversely, proteins encoded by the Francisella pathogenicity island (FPI) were downregulated, supporting the importance of the F. tularensis Type VI Secretion System for vacuole escape, not cytosolic replication. In the host cell, over 300 proteins showed differential expression among the 6200 identified during infection. The most upregulated host protein was cis-aconitate decarboxylase IRG1, known for itaconate production with antimicrobial properties in Francisella. Surprisingly, disrupting IRG1 expression did not impact Francisella's intracellular life cycle, suggesting redundancy with other immune proteins or inclusion in larger complexes. Over-representation analysis highlighted cell-cell contact and actin polymerization in macrophage deregulated proteins. Using flow cytometry and live cell imaging, we demonstrated that merocytophagy involves diverse cell-to-cell contacts and actin polymerization-dependent processes. These findings lay the groundwork for further exploration of merocytophagy and its molecular mechanisms in future research.Data are available via ProteomeXchange with identifier PXD035145.

Using the pan-genomic framework for the discovery of genomic islands in the haloarchaeon Halorubrum ezzemoulense.

In this study, we use pan-genomics to characterize the genomic variability of the widely dispersed halophilic archaeal species Halorubrum ezzemoulense (Hez). We include a multi-regional sampling of newly sequenced, high-quality draft genomes. The pan-genome graph of the species reveals 50 genomic islands that represent rare accessory genetic capabilities available to members. Most notably, we observe rearrangements that have led to the insertion/recombination/replacement of mutually exclusive genomic islands in equivalent genome positions ("homeocassettes"). These conflicting islands encode for similar functions, but homologs from islands located between the same core genes exhibit high divergence on the amino acid level, while the neighboring core genes are nearly identical. Both islands of a homeocassette often coexist in the same geographic location, suggesting that either island may be beyond the reach of selective sweeps and that these loci of divergence between Hez members are maintained and persist long term. This implies that subsections of the population have different niche preferences and rare metabolic capabilities. After an evaluation of the gene content in the homeocassettes, we speculate that these islands may play a role in the speciation, niche adaptability, and group selection dynamics in Hez. Though homeocassettes are first described in this study, similar replacements and divergence of genes on genomic islands have been previously reported in other Haloarchaea and distantly related Archaea, suggesting that homeocassettes may be a feature in a wide range of organisms outside of Hez.IMPORTANCEThis study catalogs the rare genes discovered in strains of the species Halorubrum ezzemoulense (Hez), an obligate halophilic archaeon, through the perspective of its pan-genome. These rare genes are often found to be arranged on islands that confer metabolic and transport functions and contain genes that have eluded previous studies. The discovery of divergent, but homologous islands occupying equivalent genome positions ("homeocassettes") in different genomes, reveals significant new information on genome evolution in Hez. Homeocassette pairs encode for similar functions, but their dissimilarity and distribution imply high rates of recombination, different specializations, and niche preferences in Hez. The coexistence of both islands of a homeocassette pair in multiple environments demonstrates that both islands are beyond the reach of selective sweeps and that these genome content differences between strains persist long term. The switch between islands through recombination under different environmental conditions may lead to a greater range of niche adaptability in Hez.

DICEP: An integrative approach to augmenting genomic island detection.

Mobilization of clusters of genes called genomic islands (GIs) across bacterial lineages facilitates dissemination of traits, such as, resistance against antibiotics, virulence or hypervirulence, and versatile metabolic capabilities. Robust delineation of GIs is critical to understanding bacterial evolution that has a vast impact on different life forms. Methods for identification of GIs exploit different evolutionary features or signals encoded within the genomes of bacteria, however, the current state-of-the-art in GI detection still leaves much to be desired. Here, we have taken a combinatorial approach that accounted for GI specific features such as compositional bias, aberrant phyletic pattern, and marker gene enrichment within an integrative framework to delineate GIs in bacterial genomes. Our GI prediction tool, DICEP, was assessed on simulated genomes and well-characterized bacterial genomes. DICEP compared favorably with current GI detection tools on real and synthetic datasets.

Identification of a multidrug resistance genomic island harboring a nonfunctional optrA gene in Campylobacter coli of chicken origin.

Campylobacter spp., such as Campylobacter jejuni and Campylobacter coli, are important zoonotic Gram-negative pathogens that cause acute intestinal diseases in humans. In this study, a retrospective analysis was conducted on previously collected Campylobacter isolates from antimicrobial resistance surveillance. A total of 29 optrA-positive C. coli strains were identified and subjected to second-generation sequencing. Multilocus sequence typing and single nucleotide polymorphism analyses demonstrated that the 29 optrA-positive isolates were genetically homogeneous. Notably, among the 29 isolated strains, the ΔoptrA variants exhibit a nonsense mutation at position 979 where the base C is substituted by T, leading to the formation of a premature termination codon. The alignment of sequences and genetic environmental characteristics suggested that ΔoptrA located on a chromosomally carried multidrug-resistant genomic island. There are other resistant genes on the multidrug resistance genomic island, such as aph(2'')-If, aph(3')-III, aadE, tet(O), tet(L), cat, erm(A), optrA and blaOXA-61. As a result, the 29 ΔoptrA-positive strains displayed susceptibility to both florfenicol and linezolid. The ΔoptrA gene is linked to the erm(A) gene, resulting in the formation of translocatable unit (TU) that are encompassed by two copies of IS1216 mobile elements. Multiple occurrences of similar TUs have been documented in numerous C. coli and provided evidence for the significance of TUs in facilitating the transfer of drug resistance genes in C. coli.

A novel variant in Salmonella genomic island 1 of multidrug-resistant Salmonella enterica serovar Kentucky ST198.

Salmonella enterica serovar Kentucky ST198 is a major health threat due to its resistance to ciprofloxacin and several other drugs, including third-generation cephalosporins. Many drug-resistant genes have been identified in the Salmonella genomic island 1 variant K (SGI1-K). In this study, we investigated the antimicrobial resistance (AMR) profile and genotypic relatedness of two isolates of ciprofloxacin-resistant (CIPR) S. Kentucky ST198 from poultry in Northeastern Thailand. We successfully assembled the complete genomes of both isolates, namely SSSE-01 and SSSE-03, using hybrid de novo assembly of both short- and long-read sequence data. The complete genomes revealed their highly similar genomic structures and a novel variant of SGI1-K underlying multidrug-resistant (MDR) patterns, including the presence of blaTEM-1b, which confers resistance to beta-lactams, including cephalosporins and lnu(F) which confers resistance to lincomycin and other lincosamides. In addition, the chromosomal mutations in the quinolone resistance-determining region (QRDR) were found at positions 83 (Ser83Phe) and 87 (Asp87Asn) of GyrA and at positions 57 (Thr57Ser) and 80 (Ser80Ile) of ParC suggesting high resistance to ciprofloxacin. We also compared SSSE-01 and SSSE-03 with publicly available complete genome data and revealed significant variations in SGI1-K genetic structures and variable relationships to antibiotic resistance. In comparison to the other isolates, SGI1-K of SSSE-01 and SSSE-03 had a relatively large deletion in the backbone, spanning from S011 (traG∆) to S027 (resG), and the inversion of the IS26-S044∆-yidY segment. Their MDR region was characterized by the inversion of a large segment, including the mer operon and the relocation of IntI1 and several resistance genes downstream of the IS26-S044∆-yidY segment. These structural changes were likely mediated by the recombination of IS26. The findings broaden our understanding of the possible evolution pathway of SGI1-K in fostering drug resistance, which may provide opportunities to control these MDR strains.IMPORTANCEThe emergence of ciprofloxacin-resistant (CIPR) Salmonella Kentucky ST198 globally has raised significant concerns. This study focuses on two poultry isolates from Thailand, revealing a distinct Salmonella genomic island 1 variant K (SGI1-K) genetic structure. Remarkably, multiple antibiotic resistance genes (ARGs) were identified within the SGI1-K as well as other locations in the chromosome, but not in plasmids. Comparing the SGI1-K genetic structures among global and even within-country isolates unveiled substantial variations. Intriguingly, certain isolates lacked ARGs within the SGI1-K, while others had ARGs relocated outside. The presence of chromosomal extended-spectrum ß-lactamase (ESBL) genes and lincosamide resistance, lnu(F), gene, could potentially inform the choices of the treatment of CIPRS. Kentucky ST198 infections in humans. This study highlights the importance of understanding the diverse genetic structures of SGI1-K and emphasizes the role of animals and humans in the emergence of antimicrobial resistance.

German coasts harbor non-O1/non-O139 Vibrio cholerae with clinical virulence gene profiles.

Non-O1/non-O139 Vibrio cholerae (NOVC) are ubiquitous in aquatic ecosystems. In rare cases, they can cause intestinal and extra-intestinal infections in human. This ability is associated with various virulence factors. The presence of NOVC in German North Sea and Baltic Sea was observed in previous studies. However, data on virulence characteristics are still scarce. Therefore, this work aimed to investigating the virulence potential of NOVC isolated in these two regions. In total, 31 NOVC strains were collected and subjected to whole genome sequencing. In silico analysis of the pathogenic potential was performed based on the detection of genes involved in colonization and virulence. Phenotypic assays, including biofilm formation, mobility and human serum resistance assays were applied for validation. Associated toxin genes (hlyA, rtxA, chxA and stn), pathogenicity islands (Vibrio pathogenicity island 2 (VPI-II) and Vibrio seventh pathogenicity island 2 (VSP-II)) and secretion systems (Type II, III and VI secretion system) were observed. A maximum likelihood analysis from shared core genes revealed a close relationship between clinical NOVCs published in NCBI and environmental strains from this study. NOVC strains are more mobile at 37 °C than at 25 °C, and 68% of the NOVC strains could form strong biofilms at both temperatures. All tested strains were able to lyse erythrocytes from both human and sheep blood. Additionally, one strain could survive up to 60% and seven strains up to 40% human serum at 37 °C. Overall, the genetic virulence profile as well as the phenotypic virulence characteristics of the investigated NOVC from the German North Sea and Baltic Sea suggest potential human pathogenicity.

Hypervirulent Klebsiella pneumoniae employs genomic island encoded toxins against bacterial competitors in the gut.

The hypervirulent lineages of Klebsiella pneumoniae (HvKp) cause invasive infections such as Klebsiella-liver abscess. Invasive infection often occurs after initial colonization of the host gastrointestinal tract by HvKp. Over 80% of HvKp isolates belong to the clonal group 23 sublineage I that has acquired genomic islands (GIs) GIE492 and ICEKp10. Our analysis of 12 361 K. pneumoniae genomes revealed that GIs GIE492 and ICEKp10 are co-associated with the CG23-I and CG10118 HvKp lineages. GIE492 and ICEKp10 enable HvKp to make a functional bacteriocin microcin E492 (mccE492) and the genotoxin colibactin, respectively. We discovered that GIE492 and ICEKp10 play cooperative roles and enhance gastrointestinal colonization by HvKp. Colibactin is the primary driver of this effect, modifying gut microbiome diversity. Our in vitro assays demonstrate that colibactin and mccE492 kill or inhibit a range of Gram-negative Klebsiella species and Escherichia coli strains, including Gram-positive bacteria, sometimes cooperatively. Moreover, mccE492 and colibactin kill human anaerobic gut commensals that are similar to the taxa found altered by colibactin in the mouse intestines. Our findings suggest that GIs GIE492 and ICEKp10 enable HvKp to kill several commensal bacterial taxa during interspecies interactions in the gut. Thus, acquisition of GIE492 and ICEKp10 could enable better carriage in host populations and explain the dominance of the CG23-I HvKp lineage.

Aminoglycoside resistance genes in early members of the Acinetobacter baumannii ST78A (SMAL, Italian clone) reside in an IS26-bounded island in the chromosome.

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.

VIM-type metallo-ß-lactamase (MBL)-encoding genomic islands in Pseudomonas spp. in Poland: predominance of clc-like integrative and conjugative elements (ICEs).

OBJECTIVES: To characterize VIM-type metallo-ß-lactamase (MBL)-encoding genomic islands (GIs) in Pseudomonas aeruginosa and P. putida group isolates from Polish hospitals from 2001-2015/16. METHODS: Twelve P. aeruginosa and 20 P. putida group isolates producing VIM-like MBLs were selected from a large collection of these based on epidemiological and typing data. The organisms represented all major epidemic genotypes of these species spread in Poland with chromosomally located blaVIM gene-carrying integrons. The previously determined short-read sequences were complemented by long-read sequencing in this study. The comparative structural analysis of the GIs used a variety of bioinformatic tools. RESULTS: Thirty different GIs with blaVIM integrons were identified in the 32 isolates, of which 24 GIs from 26 isolates were integrative and conjugative elements (ICEs) of the clc family. These in turn were dominated by 21 variants of the GI2/ICE6441 subfamily with a total of 19 VIM integrons, each inserted in the same position within the ICE's Tn21-like transposon Tn4380. The three other ICEs formed a novel ICE6705 subfamily, lacking Tn4380 and having different VIM integrons located in another site of the elements. The remaining six non-ICE GIs represented miscellaneous structures. The presence of various integrons in the same ICE sublineage, and of the same integron in different GIs, indicated circulation and recombination of the integron-carrying genetic platforms across Pseudomonas species/genotypes. CONCLUSIONS: Despite the general diversity of the blaVIM-carrying GIs in Pseudomonas spp. in Poland, a clear predominance of broadly spread and rapidly evolving clc-type ICEs was documented, confirming their significant role in antimicrobial resistance epidemiology.

Combining CRISPR/Cas9 and natural excision for the precise and complete removal of mobile genetic elements in bacteria.

Horizontal gene transfer, facilitated by mobile genetic elements (MGEs), is an adaptive evolutionary process that contributes to the evolution of bacterial populations and infectious diseases. A variety of MGEs not only can integrate into the bacterial genome but also can survive or even replicate like plasmids in the cytoplasm, thus requiring precise and complete removal for studying their strategies in benefiting host cells. Existing methods for MGE removal, such as homologous recombination-based deletion and excisionase-based methods, have limitations in effectively eliminating certain MGEs. To overcome these limitations, we developed the Cas9-NE method, which combines the CRISPR/Cas9 system with the natural excision of MGEs. In this approach, a specialized single guide RNA (sgRNA) element is designed with a 20-nucleotide region that pairs with the MGE sequence. This sgRNA is expressed from a plasmid that also carries the Cas9 gene. By utilizing the Cas9-NE method, both the integrative and circular forms of MGEs can be precisely and completely eliminated through Cas9 cleavage, generating MGE-removed cells. We have successfully applied the Cas9-NE method to remove four representative MGEs, including plasmids, prophages, and genomic islands, from Vibrio strains. This new approach not only enables various investigations on MGEs but also has significant implications for the rapid generation of strains for commercial purposes.IMPORTANCEMobile genetic elements (MGEs) are of utmost importance for bacterial adaptation and pathogenicity, existing in various forms and multiple copies within bacterial cells. Integrated MGEs play dual roles in bacterial hosts, enhancing the fitness of the host by delivering cargo genes and potentially modifying the bacterial genome through the integration/excision process. This process can lead to alterations in promoters or coding sequences or even gene disruptions at integration sites, influencing the physiological functions of host bacteria. Here, we developed a new approach called Cas9-NE, allowing them to maintain the natural sequence changes associated with MGE excision. Cas9-NE allows the one-step removal of integrated and circular MGEs, addressing the challenge of eliminating various MGE forms efficiently. This approach simplifies MGE elimination in bacteria, expediting research on MGEs.

Genomic islands of speciation harbor genes underlying coloration differences in a pair of Neotropical seedeaters.

Incomplete speciation can be leveraged to associate phenotypes with genotypes, thus providing insights into the traits relevant to the reproductive isolation of diverging taxa. We investigate the genetic underpinnings of the phenotypic differences between Sporophila plumbea and Sporophila beltoni. Sporophila beltoni has only recently been described based, most notably, based on differences in bill coloration (yellow vs. black in S. plumbea). Both species are indistinguishable through mtDNA or reduced-representation genomic data, and even whole-genome sequencing revealed low genetic differentiation. Demographic reconstructions attribute this genetic homogeneity to gene flow, despite divergence in the order of millions of generations. We found a narrow hybrid zone in southern Brazil where genetically, yet not phenotypically, admixed individuals appear to be prevalent. Despite the overall low genetic differentiation, we identified 3 narrow peaks along the genome with highly differentiated SNPs. These regions harbor 6 genes, one of which is involved in pigmentation (EDN3) and is a candidate for controlling bill color. Within the outlier peaks, we found signatures of resistance to gene flow, as expected for islands of speciation. Our study shows how genes related to coloration traits are likely involved in generating prezygotic isolation and establishing species boundaries early in speciation.

Speeding genomic island discovery through systematic design of reference database composition.

BACKGROUND: Genomic islands (GIs) are mobile genetic elements that integrate site-specifically into bacterial chromosomes, bearing genes that affect phenotypes such as pathogenicity and metabolism. GIs typically occur sporadically among related bacterial strains, enabling comparative genomic approaches to GI identification. For a candidate GI in a query genome, the number of reference genomes with a precise deletion of the GI serves as a support value for the GI. Our comparative software for GI identification was slowed by our original use of large reference genome databases (DBs). Here we explore smaller species-focused DBs. RESULTS: With increasing DB size, recovery of our reliable prophage GI calls reached a plateau, while recovery of less reliable GI calls (FPs) increased rapidly as DB sizes exceeded ~500 genomes; i.e., overlarge DBs can increase FP rates. Paradoxically, relative to prophages, FPs were both more frequently supported only by genomes outside the species and more frequently supported only by genomes inside the species; this may be due to their generally lower support values. Setting a DB size limit for our SMAll Ranked Tailored (SMART) DB design speeded runtime ~65-fold. Strictly intra-species DBs would tend to lower yields of prophages for small species (with few genomes available); simulations with large species showed that this could be partially overcome by reaching outside the species to closely related taxa, without an FP burden. Employing such taxonomic outreach in DB design generated redundancy in the DB set; as few as 2984 DBs were needed to cover all 47894 prokaryotic species. CONCLUSIONS: Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.