Genomic Islands Research Articles

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 361K. 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.

Flexible genomic island conservation across freshwater and marine Methylophilaceae.

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Discovering genomic islands in unannotated bacterial genomes using sequence embedding.

Genomic islands (GEIs) are clusters of genes in bacterial genomes that are typically acquired by horizontal gene transfer. GEIs play a crucial role in the evolution of bacteria by rapidly introducing genetic diversity and thus helping them adapt to changing environments. Specifically of interest to human health, many GEIs contain pathogenicity and antimicrobial resistance genes. Detecting GEIs is, therefore, an important problem in biomedical and environmental research. There have been many previous studies for computationally identifying GEIs. Still, most of these studies rely on detecting anomalies in the unannotated nucleotide sequences or on a fixed set of known features on annotated nucleotide sequences. Here, we present TreasureIsland, which uses a new unsupervised representation of DNA sequences to predict GEIs. We developed a high-precision boundary detection method featuring an incremental fine-tuning of GEI borders, and we evaluated the accuracy of this framework using a new comprehensive reference dataset, Benbow. We show that TreasureIsland's accuracy rivals other GEI predictors, enabling efficient and faster identification of GEIs in unannotated bacterial genomes. TreasureIsland is available under an MIT license at: https://github.com/FriedbergLab/GenomicIslandPrediction.

Environmental Compatibility and Genome Flexibility of Klebsiella oxytoca Isolated from Eight Species of Aquatic Animals

Klebsiella oxytoca is an emerging pathogen that can cause life-threatening infectious diseases in humans. Recently, we firstly reported for the first time the presence of K. oxytoca in edible aquatic animals. In this study, we further investigated its bacterial environmental fitness and genome evolution signatures. The results revealed that K. oxytoca isolates (n = 8), originating from eight species of aquatic animals, were capable of growing under a broad spectrum of environmental conditions (pH 4.5–8.5, 0.5–6.5% NaCl), with different biofilm formation and swimming mobility profiles. The genome sequences of the K. oxytoca isolates were determined (5.84–6.02 Mb, 55.07–56.06% GC content). Strikingly, numerous putative mobile genetic elements (MGEs), particularly genomic islands (GIs, n = 105) and prophages (n = 24), were found in the K. oxytoca genomes, which provided the bacterium with specific adaptation traits, such as resistance, virulence, and material metabolism. Interestingly, the identified prophage-related clusters were derived from Burkholderia spp., Enterobacter spp., Klebsiella spp., Pseudomonas spp., and Haemophilus spp., suggesting phage transmission across Klebsiella and the other four genera. Many strain-specific (n = 10–447) genes were present in the K. oxytoca genomes, whereas the CRISPR-Cas protein-encoding gene was absent, indicating likely active horizontal gene transfer (HGT) and considerable genome variation in K. oxytoca evolution. Overall, the results of this study are the first to demonstrate the environmental compatibility and genome flexibility of K. oxytoca of aquatic animal origins.

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

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. 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. Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.

Genomic analysis of Oceanotoga teriensis strain UFV_LIMV02, a multidrug-resistant thermophilic bacterium isolated from an offshore oil reservoir.

Bacteria of the species Oceanotoga teriensis belong to the family Petrotogaceae, are Gram-negative bacilli, are moderately thermophilic and are included in the group of thiosulfate-reducing bacteria, being capable of significantly accelerating corrosion in metallic structures. However, no in-depth study on the genome, antibiotic resistance and mobile elements has been carried out so far. In this work, the isolation, phenotypic and genotypic characterization of the multi-resistant O. teriensis UFV_LIMV02 strain was carried out, from water samples from an offshore oil extraction platform in Rio de Janeiro (Brazil). We determined that the isolate has a genome of 2 812 778 bp in size, with 26 % GC content, organized into 34 contigs. Genomic annotation using Rapid Annotation using Subsystem Technology revealed the presence of genes related to resistance to antibiotics and heavy metals. By evaluating the antimicrobial resistance of the isolate using the disc diffusion technique, resistance was verified for the classes of antibiotics, beta-lactams, fluoroquinolones, aminoglycosides, sulfonamides, lincosamides and rifamycins, a total of 14 antibiotics. The search for genomic islands, prophages and defence systems against phage infection revealed the presence of five genomic islands in its genome, containing genes related to resistance to heavy metals and antibiotics, most of which are efflux pumps and several transposases. No prophage was found in its genome; however, nine different defence systems against phage infection were detected. When analysing the clustered regularly interspaced short palindromic repeat (CRISPR) systems, four CRISPR arrays, classified as types I-B and III-B, with 272 spacers, can provide the strain with immunity to different mobile genetic elements and bacteriophage infection. The results found in this study show that the isolate UFV_LIVM02 is an environmental bacterium, resistant to different classes of antibiotics, and that the proteins encoded by the predicted genomic islands may be associated with the development of greater resistance to antibiotics and heavy metals. They provide evidence that environmental bacteria found in offshore oil exploration residues may pose a risk for the spread of antibiotic resistance genes. More comprehensive studies on the microbial community present in oil waste are needed to assess the risks of horizontal gene transfer.

Genomic Insights into the Pathogenicity and Drug-Resistance of a Bacillus cereus Isolated from Human Teeth.

Bacillus cereus is a common bacterium found in the environment. Some strains can cause food poisoning, and very few can cause clinically severe infections, leading to death. Here, we characterized the genome sequence of B. cereus LIN78 isolated from teeth with deep caries and compared it with those of 25 other related species. Third-generation sequencing technology, bacteriological analyses, biochemistry, and mass spectrometry were applied to characterize the drug-resistance genes and virulence factors of B. cereus LIN78. The complete genome sequence of B. cereus Lin78 consists of 5647 genes distributed on a circular chromosome, a 393 kbp plasmid, and 928 pseudogenes (37.4% of whole-genome DNA). The LIN78 genome contains 14 sets of 16s, 23s, and 5s ribosomal RNA operons; 106 tRNA genes, one tmRNA, 12 genomic islands, six prophases, 64 repeats; 37 antibiotic-resistant genes; and 1119 putative virulence genes, including enterotoxins and cytolysins. The B. cereus LIN78 genome carries multiple copies of non-ribosomal polypeptide synthetase (NRPS) and post-translationally modified peptides (RiPPs). Phylogenetic analysis of the 26 B. cereus strains showed that B. cereus LIN78 is evolutionarily closely related to B. thuringiensis ATCC 10792 and B. cereus ATCC 14579. The newly isolated B. cereus carries many virulence genes, including enterotoxins and hemolysins, similar to B. anthracis, and multiple antibiotic resistance genes. These findings suggest that the strain has a potential risk of causing disease. Our studies are vital for further exploration of the evolution of B. cereus, its pathogenic mechanisms, and the control and treatment of bacterial infections.

Whole-genome comparison using complete genomes from Campylobacter fetus strains revealed single nucleotide polymorphisms on non-genomic islands for subspecies differentiation.

Bovine Genital Campylobacteriosis (BGC), caused by Campylobacter fetus subsp. venerealis, is a sexually transmitted bacterium that significantly impacts cattle reproductive performance. However, current detection methods lack consistency and reliability due to the close genetic similarity between C. fetus subsp. venerealis and C. fetus subsp. fetus. Therefore, this study aimed to utilize complete genome analysis to distinguish genetic features between C. fetus subsp. venerealis and other subspecies, thereby enhancing BGC detection for routine screening and epidemiological studies. This study reported the complete genomes of four C. fetus subsp. fetus and five C. fetus subsp. venerealis, sequenced using long-read sequencing technologies. Comparative whole-genome analyses (n = 25) were conducted, incorporating an additional 16 complete C. fetus genomes from the NCBI database, to investigate the genomic differences between these two closely related C. fetus subspecies. Pan-genomic analyses revealed a core genome consisting of 1,561 genes and an accessory pangenome of 1,064 genes between the two C. fetus subspecies. However, no unique predicted genes were identified in either subspecies. Nonetheless, whole-genome single nucleotide polymorphisms (SNPs) analysis identified 289 SNPs unique to one or the C. fetus subspecies. After the removal of SNPs located on putative genomic islands, recombination sites, and those causing synonymous amino acid changes, the remaining 184 SNPs were functionally annotated. Candidate SNPs that were annotated with the KEGG "Peptidoglycan Biosynthesis" pathway were recruited for further analysis due to their potential association with the glycine intolerance characteristic of C. fetus subsp. venerealis and its biovar variant. Verification with 58 annotated C. fetus genomes, both complete and incomplete, from RefSeq, successfully classified these seven SNPs into two groups, aligning with their phenotypic identification as CFF (Campylobacter fetus subsp. fetus) or CFV/CFVi (Campylobacter fetus subsp. venerealis and its biovar variant). Furthermore, we demonstrated the application of mraY SNPs for detecting C. fetus subspecies using a quantitative PCR assay. Our results highlighted the high genetic stability of C. fetus subspecies. Nevertheless, Campylobacter fetus subsp. venerealis and its biovar variants encoded common SNPs in genes related to glycine intolerance, which differentiates them from C. fetus subsp. fetus. This discovery highlights the potential of employing a multiple-SNP assay for the precise differentiation of C. fetus subspecies.

Expression of accessory genes in Salmonella requires the presence of the Gre factors

Genomic studies with Salmonella enterica serovar Typhimurium reveal a crucial role of horizontal gene transfer (HGT) in the acquisition of accessory cellular functions involved in host-interaction. Many virulence genes are located in genomic islands, plasmids and prophages. GreA and GreB proteins, Gre factors, interact transiently with the RNA polymerase alleviating backtracked complexes during transcription elongation. The overall effect of Gre factors depletion in Salmonella expression profile was studied. Both proteins are functionally redundant since only when both Gre factors were depleted a major effect in gene expression was detected. Remarkably, the accessory gene pool is particularly sensitive to the lack of Gre factors, with 18.6% of accessory genes stimulated by the Gre factors versus 4.4% of core genome genes. Gre factors involvement is particularly relevant for the expression of genes located in genomic islands. Our data reveal that Gre factors are required for the expression of accessory genes.

The horizontal gene transfer of perchlorate reduction genomic island in three bacteria from an ecological niche.

Three new strains of dissimilatory perchlorate-reducing bacteria (DPRB), QD19-16, QD1-5, and P3-1, were isolated from an active sludge. Phylogenetic trees based on 16S rRNA genes indicated that QD19-16, QD1-5, and P3-1 belonged to Brucella, Acidovorax, and Citrobacter, respectively, expanding the distribution of DPRB in the Proteobacteria. The three strains were gram-negative and facultative anaerobes with rod-shaped cells without flagella, which were 1.0-1.6μm long and 0.5-0.6μm wide. The three DPRB strains utilized similar broad spectrum of electron donors and acceptors and demonstrated a similar capability to reduce perchlorate within 6days. The enzyme activity of perchlorate reductase in QD19-16 toward chlorate was higher than that toward perchlorate. The high sequence similarity of the perchlorate reductase operon and chlorite dismutase genes in the perchlorate reduction genomic islands (PRI) of the three strains implied that they were monophyletic origin from a common ancestral PRI. Two transposase genes (tnp1 and tnp2) were found in the PRIs of strain QD19-16 and QD1-5, but were absent in the strain P3-1 PRI. The presence of fragments of IR sequences in the P3-1 PRI suggested that P3-1 PRI had previously contained these two tnp genes. Therefore, it is plausible to suggest that a common ancestral PRI transferred across the strains Brucella sp. QD19-16, Acidovorax sp. QD1-5, and Citrobacter sp. P3-1 through horizontal gene transfer, facilitated by transposases. These results provided a direct evidence of horizontal gene transfer of PRI that could jump across phylogenetically unrelated bacteria through transposase. KEY POINTS: • Three new DPRB strains can effectively remove high concentration of perchlorate. • The PRIs of three DPRB strains are acquired from a single ancestral PRI. • PRIs are incorporated into different bacteria genome through HGT by transposase.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59.

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Pan-genome analysis of six Paracoccus type strain genomes reveal lifestyle traits.

The genus Paracoccus capable of inhabiting a variety of different ecological niches both, marine and terrestrial, is globally distributed. In addition, Paracoccus is taxonomically, metabolically and regarding lifestyle highly diverse. Until now, little is known on how Paracoccus can adapt to such a range of different ecological niches and lifestyles. In the present study, the genus Paracoccus was phylogenomically analyzed (n = 160) and revisited, allowing species level classification of 16 so far unclassified Paracoccus sp. strains and detection of five misclassifications. Moreover, we performed pan-genome analysis of Paracoccus-type strains, isolated from a variety of ecological niches, including different soils, tidal flat sediment, host association such as the bluespotted cornetfish, Bugula plumosa, and the reef-building coral Stylophora pistillata to elucidate either i) the importance of lifestyle and adaptation potential, and ii) the role of the genomic equipment and niche adaptation potential. Six complete genomes were de novo hybrid assembled using a combination of short and long-read technologies. These Paracoccus genomes increase the number of completely closed high-quality genomes of type strains from 15 to 21. Pan-genome analysis revealed an open pan-genome composed of 13,819 genes with a minimal chromosomal core (8.84%) highlighting the genomic adaptation potential and the huge impact of extra-chromosomal elements. All genomes are shaped by the acquisition of various mobile genetic elements including genomic islands, prophages, transposases, and insertion sequences emphasizing their genomic plasticity. In terms of lifestyle, each mobile genetic elements should be evaluated separately with respect to the ecological context. Free-living genomes, in contrast to host-associated, tend to comprise (1) larger genomes, or the highest number of extra-chromosomal elements, (2) higher number of genomic islands and insertion sequence elements, and (3) a lower number of intact prophage regions. Regarding lifestyle adaptations, free-living genomes share genes linked to genetic exchange via T4SS, especially relevant for Paracoccus, known for their numerous extrachromosomal elements, enabling adaptation to dynamic environments. Conversely, host-associated genomes feature diverse genes involved in molecule transport, cell wall modification, attachment, stress protection, DNA repair, carbon, and nitrogen metabolism. Due to the vast number of adaptive genes, Paracoccus can quickly adapt to changing environmental conditions.

Comparative genomics analysis and characterization of Shiga toxin-producing Escherichia coli O157:H7 strains reveal virulence genes, resistance genes, prophages and plasmids

Escherichia coli O157:H7 is a foodborne pathogen that has been linked to global disease outbreaks. These diseases include hemorrhagic colitis and hemolytic uremic syndrome. It is vital to know the features that make this strain pathogenic to understand the development of disease outbreaks. In the current study, a comparative genomic analysis was carried out to determine the presence of structural and functional features of O157:H7 strains obtained from 115 National Center for Biotechnology Information database. These strains of interest were analysed in the following programs: BLAST Ring Image Generator, PlasmidFinder, ResFinder, VirulenceFinder, IslandViewer 4 and PHASTER. Five strains (ECP19–198, ECP19–798, F7508, F8952, H2495) demonstrated a great homology with Sakai because of a few regions missing. Five resistant genes were identified, however, Macrolide-associated resistance gene mdf(A) was commonly found in all genomes. Majority of the strains (97%) were positive for 15 of the virulent genes (espA, espB, espF, espJ, gad, chuA, eae, iss, nleA, nleB, nleC, ompT, tccP, terC and tir). The plasmid analysis demonstrated that the IncF group was the most prevalent in the strains analysed. The prophage and genomic island analysis showed a distribution of bacteriophages and genomic islands respectively. The results indicated that structural and functional features of the many O157:H7 strains differ and may be a result of obtaining mobile genetic elements via horizontal gene transfer. Understanding the evolution of O157:H7 strains pathogenicity in terms of their structural and functional features will enable the development of detection and control of transmission strategies.

Evolutionary trends of carbapenem-resistant and susceptible Acinetobacter baumannii isolates in a major tertiary care setting from North India

Emergence of carbapenem-resistant A. baumannii (CRAB) is a global, ongoing healthcare concern. CRAB is among the topmost priority pathogens, with various studies focusing on its global population structure and resistant allelic profiles. However, carbapenem-susceptible A. baumannii (CSAB) isolates are often overlooked due to their sensitivity to beta-lactams, which can provide important insights into origin of CRAB lineages and isolates. In the present study, we report genomic investigation of CRAB and CSAB coexisting in Indian hospital setting. MLST based population structure and phylogenomics suggest they mainly follow distinct evolutionary routes forming two phylogroups. PG-I exclusively for a successful clone (ST2) of CRAB and PG-II comprises diversified CSAB isolates except PG3373, which is CRAB. Additionally, there are few CRAB isolates not belonging to PG-I and sharing clonal relationship with CSAB isolates indicating role of genome plasticity towards extensive drug resistance in the nosocomial environment. Further, genealogical analysis depicts prominent role of recombination in emergence and evolution of a major CRAB lineage. Further, CRAB isolates are enriched in resistomes as compared to CSAB isolates, which were encoded on the genomic island. Such comparative genomic insights will aid in our understanding and localized management of rapidly evolving pandrug resistant nosocomial pathogens.

Clonality and Diversity in the Soft Rot Dickeya solani Phytopathogen.

Bacterial diversity analyses often suffer from a bias due to sampling only from a limited number of hosts or narrow geographic locations. This was the case for the phytopathogenic species Dickeya solani, whose members were mainly isolated from a few hosts-potato and ornamentals-and from the same geographical area-Europe and Israel, which are connected by seed trade. Most D. solani members were clonal with the notable exception of the potato isolate RNS05.1.2A and two related strains that are clearly distinct from other D. solani genomes. To investigate if D. solani genomic diversity might be broadened by analysis of strains isolated from other environments, we analysed new strains isolated from ornamentals and from river water as well as strain CFBP 5647 isolated from tomato in the Caribbean island Guadeloupe. While water strains were clonal to RNS05.1.2A, the Caribbean tomato strain formed a third clade. The genomes of the three clades are highly syntenic; they shared almost 3900 protein families, and clade-specific genes were mainly included in genomic islands of extrachromosomal origin. Our study thus revealed both broader D. solani diversity with the characterisation of a third clade isolated in Latin America and a very high genomic conservation between clade members.

In silico, in vitro and in vivo safety assessment of three Limosilactobacillus mucosae strains for probiotic candidates

Previous studies indicated that Limosilactobacillus mucosae had a good ability to regulate lipid metabolism, however, there was limited research on its safety. In this study, we assessed the safety of three L. mucosae strains including genotypes and phenotype, subacute oral toxicity test, and their influences on immunosuppressed mice. The results showed that there were no virulence genes or genetic islands with clear pathogenicity in the genomes of L. mucosae DCC1HL5, FZJTZ26M3, and GD16M9. All three strains were resistant to chloramphenicol, additionally, DCC1HL5 was resistant to kanamycin and tetracycline, but the possibility of horizontal transfer of their resistance genes was relatively low. L. mucosae DCC1HL5, FZJTZ26M3, and GD16M9 did not produce biogenic amine, and α-hemolysis was found. The results of the 28-day subacute oral toxicity test showed that 109 CFU of each strain did not cause significant abnormalities in various aspects such as body weight, organ indexes, histopathology, hematological parameters, fasting blood glucose, serum lipids, liver parameters, and intestinal inflammation in both healthy and immunosuppressed mice. Therefore, L. mucosae DCC1HL5, FZJTZ26M3, and GD16M9 have the potential for safe application in functional foods.

Subinhibitory antibiotic concentrations promote the excision of a genomic island carried by the globally spread carbapenem-resistant Klebsiella pneumoniae sequence type 258.

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.

Remodulation of bacterial transcriptome after acquisition of foreign DNA: the case of irp-HPI high-pathogenicity island in Vibrio anguillarum.

The high-pathogenicity island irp-HPI is widespread in Vibrionaceae and encodes the siderophore piscibactin, as well as the regulator PbtA that is essential for its expression. In this work, we aim to study whether PbtA directly interacts with irp-HPI promoters. Furthermore, we hypothesize that PbtA, and thereby the acquisition of irp-HPI island, may also influence the expression of other genes elsewhere in the bacterial genome. To address this question, an RNAseq analysis was conducted to identify differentially expressed genes after pbtA deletion in Vibrio anguillarum RV22 genetic background. The results showed that PbtA not only modulates the irp-HPI genes but also modulates the expression of a plethora of V. anguillarum core genome genes, inducing nitrate, arginine, and sulfate metabolism, T6SS1, and quorum sensing, while repressing lipopolysaccharide (LPS) production, MARTX toxin, and major porins such as OmpV and ChiP. The direct binding of the C-terminal domain of PbtA to piscibactin promoters (PfrpA and PfrpC), quorum sensing (vanT), LPS transporter wza, and T6SS structure- and effector-encoding genes was demonstrated by electrophoretic mobility shift assay (EMSA). The results provide valuable insights into the regulatory mechanisms underlying the expression of irp-HPI island and its impact on Vibrios transcriptome, with implications in pathogenesis.IMPORTANCEHorizontal gene transfer enables bacteria to acquire traits, such as virulence factors, thereby increasing the risk of the emergence of new pathogens. irp-HPI genomic island has a broad dissemination in Vibrionaceae and is present in numerous potentially pathogenic marine bacteria, some of which can infect humans. Previous works showed that certain V. anguillarum strains exhibit an expanded host range plasticity and heightened virulence, a phenomenon linked to the acquisition of the irp-HPI genomic island. The present work shows that this adaptive capability is likely achieved through comprehensive changes in the transcriptome of the bacteria and that these changes are mediated by the master regulator PbtA encoded within the irp-HPI element. Our results shed light on the broad implications of horizontal gene transfer in bacterial evolution, showing that the acquired DNA can directly mediate changes in the expression of the core genome, with profounds implications in pathogenesis.

The effect of divergent and parallel selection on the genomic landscape of divergence.

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Unveiling highly efficient degradation mechanism of strobilurin fungicides by strain Hyphomicrobium sp. DY-1

Strobilurin fungicides are a type of agricultural chemical commonly used to control fungal diseases in crops. However, the abuse of strobilurin fungicides leads to heavy residue in the environment, posing serious threats to ecosystems and human health. In this regard, the strains Hyphomicrobium sp. DY-1 and Hyphomicrobium sp. B1 could hydrolyze the strobilurin fungicides effectively (trifloxystrobin, azoxystrobin, pyraclostrobin, and picoxystrobin), generating the corresponding parent acid. Therefore, the present study sought to explore the molecular degradation mechanism of strobilurin fungicides by strain DY-1. The results showed that strain DY-1 exhibited a faster degradation rate for strobilurin fungicides than strain B1. Comparative genomic analysis and functional verification of both strains revealed that strain DY-1 possesses six strobilurin hydrolases (TriH, StrH, StrH1, StrH2, StrH3, and StrH4), whereas strain B1 possesses only three strobilurin hydrolases (TriH, StrH, and StrH1). All StrHs (StrH1, StrH2, StrH3, and StrH4) shared 50%–77% identity with StrH and contained conserved catalytic motif Ser-Glu-His. Further investigation into the distribution of strobilurin fungicide hydrolases revealed that strH and strH1 were aligned together and surrounded by transposases in the strains DY-1 and B1 and a genomic island containing strH and strH1 derived from the Hyphomicrobium genus. Overall, the function, diversity, and distribution of strobilurin hydrolases will provide insights into the microbial degradation of strobilurin fungicides.