Genomic islands in pathogenic and environmental microorganisms

Integrative and conjugative elements of Pasteurella multocida: Prevalence and signatures in population evolution

Pasteurella multocida (P. multocida) is a bacterial pathogen responsible for a range of infections in humans and various animal hosts, causing significant economic losses in farming. Integrative and conjugative elements (ICEs) are important horizontal gene transfer elements, potentially enabling host bacteria to enhance adaptability by acquiring multiple functional genes. However, the understanding of ICEs in P. multocida and their impact on the transmission of this pathogen remains limited. In this study, 42 poultry-sourced P. multocida genomes obtained by high-throughput sequencing together with 393 publicly available P. multocida genomes were used to analyse the horizontal transfer of ICEs. Eighty-two ICEs were identified in P. multocida, including SXT/R391 and Tn916 subtypes, as well as three subtypes of ICEHin1056 family, with the latter being widely prevalent in P. multocida and carrying multiple resistance genes. The correlations between insertion sequences and resistant genes in ICEs were also identified, and some ICEs introduced the carbapenem gene blaOXA-2 and the bleomycin gene bleO to P. multocida. Phylogenetic and collinearity analyses of these bioinformatics found that ICEs in P. multocida were transmitted vertically and horizontally and have evolved with host specialization. These findings provide insight into the transmission and evolution mode of ICEs in P. multocida and highlight the importance of understanding these elements for controlling the spread of antibiotic resistance.

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.

Disinfection and post-disinfection conditions drive bacterial and viral evolution across the environment and host

Water disinfection practices have long been established as a critical engineering intervention for controlling pathogen transmission and safeguarding individual and public health. However, recent discoveries have unveiled the significant role disinfection and post-disinfection play in accelerating the development of resistance to disinfectants and antimicrobial drugs within bacterial and viral communities in the environment. This phenomenon, in turn, may facilitate the emergence of persistent microbes and those with new genetic characteristics. These microbes may thrive in host environments with increased infectivity and resistance, posing challenges to current medical treatments and jeopardizing human health. In this perspective, we illuminate the intricate interplay between aquatic environments, microbes, and hosts and how microbial virulence evolves across the environment and host under the pressure of disinfection and post-disinfection conditions. We aim to draw attention to the previously overlooked potential risks associated with disinfection in driving the virulence evolution of bacteria and viruses, establish connections between pathogens in diverse environments and hosts within the overarching framework of the One Health concept, and ultimately provide guidelines for advancing future water disinfection technologies to effectively curb the spread of infectious diseases.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Comparative genomics of Vibrio toranzoniae strains

Vibrio toranzoniae is a marine bacterium belonging to the Splendidus clade, originally isolated from healthy clams in Galicia (NW Spain). Its isolation from different hosts and seawater indicated two lifestyles and wide geographical distribution. The aim of the present study was to determine the differences at genome level among strains, as well as to determine their phylogeny. For this purpose, whole genomes were sequenced by different technologies and the resulting sequences corrected. Genomes were annotated and compared with different online tools. Furthermore, the study of core and pan genome was examined, and the phylogeny was inferred. The content of the core genome ranged from 2,953 to 2,766 genes and that of the pangenome from 6,278 to 6,132, depending on the tool used. The comparison revealed that although the strains shared certain homology, with DDH values ranging from 77.10 to 82.30 and values of OrthoANI higher than 97%,notable differences were found related to motility, capsule synthesis, iron acquisition system or mobile genetic elements. The phylogenetic analysis of the core genome did not reveal a differentiation of the strains according to their lifestyle, but that of the pangenome pointed out certain geographical isolation in the same growing area. The study led to a reclassification of some isolates formerly described as V. toranzoniae and manifested the importance of cured deposited sequences to proper phylogenetic assignment.

A comprehensive comparative genomic analysis revealed that plant growth promoting traits are ubiquitous in strains of Stenotrophomonas

Stenotrophomonas strains, which are often described as plant growth promoting (PGP) bacteria, are ubiquitous in many environments. A total of 213 genomes of strains of Stenotrophomonas were analyzed using comparative genomics to better understand the ecological roles of these bacteria in the environment. The pan-genome of the 213 strains of Stenotrophomonas consists of 27,186 gene families, including 710 core gene families, 11,039 unique genes and 15,437 accessory genes. Nearly all strains of Stenotrophomonas harbor the genes for GH3-family cellulose degradation and GH2- and GH31-family hemicellulose hydrolase, as well as intact glycolysis and tricarboxylic acid cycle pathways. These abilities suggest that the strains of this genus can easily obtain carbon and energy from the environment. The Stenotrophomonas strains can respond to oxidative stress by synthesizing catalase, superoxide dismutase, methionine sulfoxide reductase, and disulfide isomerase, as well as managing their osmotic balance by accumulating potassium and synthesizing compatible solutes, such as betaine, trehalose, glutamate, and proline. Each Stenotrophomonas strain also contains many genes for resistance to antibiotics and heavy metals. These genes that mediate stress tolerance increase the ability of Stenotrophomonas strains to survive in extreme environments. In addition, many functional genes related to attachment and plant colonization, growth promotion and biocontrol were identified. In detail, the genes associated with flagellar assembly, motility, chemotaxis and biofilm formation enable the strains of Stenotrophomonas to effectively colonize host plants. The presence of genes for phosphate-solubilization and siderophore production and the polyamine, indole-3-acetic acid, and cytokinin biosynthetic pathways confer the ability to promote plant growth. These strains can produce antimicrobial compounds, chitinases, lipases and proteases. Each Stenotrophomonas genome contained 1-9 prophages and 17-60 genomic islands, and the genes related to antibiotic and heavy metal resistance and the biosynthesis of polyamines, indole-3-acetic acid, and cytokinin may be acquired by horizontal gene transfer. This study demonstrates that strains of Stenotrophomonas are highly adaptable for different environments and have strong potential for use as plant growth-promoting bacteria.

Elucidating doxycycline biotransformation mechanism by Chryseobacterium sp. WX1: Multi-omics insights

Doxycycline (DOX) represents a second-generation tetracycline antibiotic that persists as a challenging-to-degrade contaminant in environmental compartments. Despite its ubiquity, scant literature exists on bacteria proficient in DOX degradation. This study marked a substantial advancement in this field by isolating Chryseobacterium sp. WX1 from an activated sludge enrichment culture, showcasing its unprecedented ability to completely degrade 50 mg/L of DOX within 44 h. Throughout the degradation process, seven biotransformation products were identified, revealing a complex pathway that began with the hydroxylation of DOX, followed by a series of transformations. Employing an integrated multi-omics approach alongside in vitro heterologous expression assays, our study distinctly identified the tetX gene as a critical facilitator of DOX hydroxylation. Proteomic analyses further pinpointed the enzymes postulated to mediate the downstream modifications of DOX hydroxylation derivatives. The elucidated degradation pathway encompassed several key biological processes, such as the microbial transmembrane transport of DOX and its intermediates, the orchestration of enzyme synthesis for transformation, energy metabolism, and other gene-regulated biological directives. This study provides the first insight into the adaptive biotransformation strategies of Chryseobacterium under DOX-induced stress, highlighting the potential applications of this strain to augment DOX removal in wastewater treatment systems containing high concentrations of DOX.

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.

Maintenance of host specialisation gradients in ectomycorrhizal symbionts

Many fungi that form ectomycorrhizas exhibit a degree of host specialisation, and individual trees are frequently colonised by communities of mycorrhizal fungi comprising species that fall on a gradient of specialisation along genetic, functional and taxonomic axes of variation. By contrast, arbuscular mycorrhizal fungi exhibit little specialisation. Here, we propose that host tree root morphology is a key factor that gives host plants fine-scale control over colonisation and therefore opportunities for driving specialisation and speciation of ectomycorrhizal fungi. A gradient in host specialisation is likely driven by four proximate mechanistic 'filters' comprising partner availability, signalling recognition, competition for colonisation, and symbiotic function (trade, rewards and sanctions), and the spatially restricted colonisation seen in heterorhizic roots enables these mechanisms, especially symbiotic function, to be more effective in driving the evolution of specialisation. We encourage manipulation experiments that integrate molecular genetics and isotope tracers to test these mechanisms, alongside mathematical simulations of eco-evolutionary dynamics in mycorrhizal symbioses.

Hidden threats in the plastisphere: Carbapenemase-producing Enterobacterales colonizing microplastics in river water

Carbapenem resistance poses a significant burden on healthcare systems worldwide. Microplastics (MPs) have emerged as potential contributors to antibiotic resistance spread in the environment. However, the link between MPs and carbapenem resistance remains unexplored. We investigated the prevalence of carbapenem-resistant bacteria colonizing MPs placed in a river. Three replicates of a mixture of polypropylene (PP), polyethylene (PE) and polyethylene terephthalate (PET) and of PET alone were placed both upstream and downstream a wastewater treatment plant (WWTP) discharge. Carbapenem-resistant Enterobacterales (CRE) were further characterized by phenotypic tests and whole-genome sequencing. The abundance of carbapenem-resistant bacteria on MPs increased significantly downstream the WWTP. Their prevalence was higher in the MPs mixture compared to PET alone. CRE strains colonizing MPs included Klebsiella pneumoniae (n = 3), Klebsiella quasipneumoniae (n = 3), Raoultella ornithinolytica (n = 2), Enterobacter kobei (n = 1) and Citrobacter freundii (n = 1), most (n = 8) recovered after the WWTP discharge. All strains exhibited at least one of the tested virulence traits (biofilm formation at 37 °C, haemolytic activity and siderophore production), were multi-drug resistant and carried carbapenemase-encoding genes [blaKPC-3 (n = 5), blaGES-5 (n = 2) or blaKPC-3 + blaGES-5 (n = 3)]. Uncommon phenotypes of resistance to imipenem/relebactam (n = 3) and ceftazidime/avibactam (n = 2) were observed. Two blaKPC-3-positive K. pneumoniae successfully transfer this gene trough conjugation. Genome analysis predicted all strains as human pathogens. The blaKPC-3 was associated with the Tn4401d transposon on a pBK30683-like plasmid in most of the isolates (n = 7). The blaGES-5 was mostly linked to class 3 integrons. A K. pneumoniae strain belonging to the outbreak-causing high-risk clone ST15 carried both blaKPC-3 and blaCTX-M-15. Two K. quasipneumoniae isolates carried the plasmid-mediated colistin resistance gene mcr-9. Our results underscore the role of MPs as vectors for CRE dissemination, particularly following WWTPs discharges. MPs may act as carriers, facilitating the dissemination of carbapenemase-encoding genes and potentially contributing to increased CRE incidence in the environment.

The analysis of complete genome sequence and comparative genomics of Vibrio parahaemolyticus LF1113 in Hainan

Vibrio parahaemolyticus is a Gram-negative, halophilic and polymorphic coccobacillus. It is world-widely distributed and has resulted in great economic losses since its first appearance. In this study, a pathogenic strain was isolated from diseased pearl gentian grouper and identified as V. parahaemolyticus based on the sequencing results of 16S rDNA gene. In order to gain a comprehensive understanding of this isolation, the whole genome sequencing was conducted. Phylogenetic analysis of the complete genomes of 16 Vibrio species showed that LF1113, ATCC17802, ATCC33787, 2210633, FORC 004, and 160807 were the most closely related. Animal experiments demonstrated that the isolated LF1113 strain was pathogenic in a fish model. This study is the first study to describe the complete genome sequence of a V. parahaemolyticus isolate, which infected pearl gentian grouper from an outbreak in a fish factory farm in Hainan. The results will expand our understanding of genetic characteristics, pathogenesis, diagnostics and disease prevention of V. parahaemolyticus, and lay the foundation for further study.

Detection and characterization of potentially hybrid enteroaggregative Escherichia coli (EAEC) strains isolated from urinary tract infection

Uropathogenic Escherichia coli (UPEC) have the potential to receive the virulence markers of intestinal pathotypes and transform into various important hybrid pathotypes. This study aimed to investigate the frequency and characteristics of hybrid enteroaggregative E. coli (EAEC)/UPEC strains. Out of 202 UPEC strains, nine (4.5%) were detected as hybrid EAEC/UPEC. These strains carried one to four iron uptake systems. Among nine investigated pathogenicity islands (PAIs), PAI IV536, PAI II536, and PAI ICFT073 were found in 9 (100%), 3 (33.3%), and 1 (11.1%) strains, respectively. The chuA and sitA genes were detected in 5 (55.5%) and 3 (33.3%) hybrid strains, respectively. Six hybrid strains were found to be typical extraintestinal pathogenic E. coli (ExPEC) according to their virulence traits. Most of the hybrid strains belonged to the phylogenetic group E (6/9). Among the hybrid strains, seven (7/9) were able to form biofilm and adhere to cells; however, only two strains penetrated into the HeLa cells. Our findings reveal some of the virulence characteristics of hybrid strains that lead to fitness and infection in the urinary tract. These strains, with virulence factors of intestinal and non-intestinal pathotypes, may become emerging pathogens in clinical settings; therefore, further studies are needed to reveal their pathogenicity mechanisms and so that preventive measures can be taken.

The integrative and conjugative element ICECiPOL15 mediates horizontal transfer of β-lactam resistance gene in Chryseobacterium indoltheticum POL15

OBJECTIVES

The dissemination of antibiotic resistance genes (ARGs) from the environment, including agricultural sources, is of increasing concern. In this study, we examined the antibiotic resistance profile and genomic sequence of a strain of Chryseobacterium indoltheticum obtained from an agricultural location.

METHODS

The multidrug-resistant bacterial strain POL15 was isolated from the wastewater of a livestock farm in China. Whole-genome sequencing was performed followed by bioinformatics analyses to identify integrative and conjugative elements (ICEs) and ARGs. Mating assays were performed to analyse ICE transferability.

RESULTS

Whole-genome sequencing and annotation showed that the genome of POL15 encodes ARGs. Additionally, an ICE named ICECiPOL15, which carries a class C β-lactamase-encoding gene blaAQU, was identified in the POL15 genome. Genes encoding an integrase, an excisionase, a relaxase, a type IV coupling protein and conjugative transposon proteins involved in a type IV secretion system were also identified in ICECiPOL15. Sequence alignment revealed that ICECiPOL15 might have evolved from other Chryseobacterium species. The horizontal transferability of ICECiPOL15 was demonstrated by mating experiments between C. indoltheticum POL15 and Escherichia coli DL21.

CONCLUSIONS

This study represents the first characterization of a mobilizable antibiotic resistance ICE in a species of C. indoltheticum and provides evidence that C. indoltheticum strains could be important reservoirs and vehicles for ARGs on livestock farms.

SMU_1361c regulates the oxidative stress response of Streptococcus mutans

Dental caries is the most common chronic infectious disease around the world and disproportionately affects the marginalized socioeconomic group. Streptococcus mutans, considered a primary etiological agent of caries, depends on the coordinated physiological response to tolerate the oxidative stress generated by commensal species within dental plaque, which is a critical aspect of its pathogenicity. Here, we identified and characterized a novel tetracycline repressor family regulator, SMU_1361c, which appears to be acquired by the bacteria via horizontal gene transfer. Surprisingly, smu_1361c functions as a negative transcriptional regulator to regulate gene expression outside its operon and is involved in the oxidative stress response of S. mutans. The smu_1361c overexpression strain UA159/pDL278-1361c was more susceptible to oxidative stress and less competitive against hydrogen peroxide generated by commensal species Streptococcus gordonii and Streptococcus sanguinis. Transcriptomics analysis revealed that smu_1361c overexpression resulted in the significant downregulation of 22 genes, mainly belonging to three gene clusters responsible for the oxidative stress response. The conversed DNA binding motif of SMU_1361c was determined by electrophoretic mobility shift and DNase I footprinting assay with purified SMU_1361c protein; therefore, smu_1361c is directly involved in gene transcription related to the oxidative stress response. Crucially, our finding provides a new understanding of how S. mutans deals with the oxidative stress that is required for pathogenesis and will facilitate the development of new and improved therapeutic approaches for dental caries.IMPORTANCEStreptococcus mutans is the major organism associated with the development of dental caries, which globally is the most common chronic disease. To persist and survive in biofilms, S. mutans must compete with commensal species that occupy the same ecological niche. Here, we uncover a novel molecular mechanism of how tetracycline repressor family regulator smu_1361c is involved in the oxidative stress response through transcriptomics analysis, electrophoretic mobility shift assay, and DNase I footprinting assay. Furthermore, we demonstrated that smu_1361c mediates S. mutans sensitivity to oxidative stress and competitiveness with commensal streptococci. Therefore, this study has revealed a previously unknown regulation between smu_1361c and genes outside its operon and demonstrated the importance of smu_1361c in the oxidative stress response and the fitness of S. mutans within the plaque biofilms, which can be exploited as a new therapy to modulate ecological homeostasis and prevent dental caries.

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.

Genomic islands and their role in fitness traits of two key sepsis-causing bacterial pathogens

To survive and establish a niche for themselves, bacteria constantly evolve. Toward that, they not only insert point mutations and promote illegitimate recombinations within their genomes but also insert pieces of 'foreign' deoxyribonucleic acid, which are commonly referred to as 'genomic islands' (GEIs). The GEIs come in several forms, structures and types, often providing a fitness advantage to the harboring bacterium. In pathogenic bacteria, some GEIs may enhance virulence, thus altering disease burden, morbidity and mortality. Hence, delineating (i) the GEIs framework, (ii) their encoded functions, (iii) the triggers that help them move, (iv) the mechanisms they exploit to move among bacteria and (v) identification of their natural reservoirs will aid in superior tackling of several bacterial diseases, including sepsis. Given the vast array of comparative genomics data, in this short review, we provide an overview of the GEIs, their types and the compositions therein, especially highlighting GEIs harbored by two important pathogens, viz. Acinetobacter baumannii and Klebsiella pneumoniae, which prominently trigger sepsis in low- and middle-income countries. Our efforts help shed some light on the challenges these pathogens pose when equipped with GEIs. We hope that this review will provoke intense research into understanding GEIs, the cues that drive their mobility across bacteria and the ways and means to prevent their transfer, especially across pathogenic bacteria.

Toxinome-the bacterial protein toxin database

IMPORTANCE

Microbes use protein toxins as important tools to attack neighboring cells, microbial or eukaryotic, and for self-killing when attacked by viruses. These toxins work through different mechanisms to inhibit cell growth or kill cells. Microbes also use antitoxin proteins to neutralize the toxin activities. Here, we developed a comprehensive database called Toxinome of nearly two million toxins and antitoxins that are encoded in 59,475 bacterial genomes. We described the distribution of bacterial toxins and identified that they are depleted by bacteria that live in hot and cold temperatures. We found 5,161 cases in which toxins and antitoxins are densely clustered in bacterial genomes and termed these areas "Toxin Islands." The Toxinome database is a useful resource for anyone interested in toxin biology and evolution, and it can guide the discovery of new toxins.

Adaptive evolution of Sphingopyxis sp. MC4 conferred degradation potential for persistent β- and δ-Hexachlorocyclohexane (HCH) isomers

Hexachlorocyclohexane (HCH), an organochlorine pesticide imposes several harmful impacts on the ecosystem. β- and δ-isomers of HCH are highly toxic, persistent, and recalcitrant to biodegradation, slow and incomplete degradation of β- and δ- isomers have been reported in a few strains. We have isolated a strain designated as Sphingopyxis strain MC4 that can tolerate and degrade high concentrations of α-, β-, γ- and δ-HCH isomers. To date, no other Sphingopyxis strain has been reported to degrade β- and δ-isomers. To understand the underlying genetic makeup contributing to adaptations, the whole genome of strain MC4 was sequenced. Comparative genome analysis showed that strain MC4 harbors the complete pathway (lin genes) required for HCH degradation. Genetic footprints such as presence of lin genes on genomic islands, IS6100 elements in close proximity of lin genes, and synteny in lin flanking regions with other strains reflects the horizontal gene transfer in strain MC4. Positive selection and HGT drive the adaptive evolution of strain MC4 under the pressure of HCH contamination that it experienced in its surrounding niche. In silico analyses showed efficient binding of β- and δ-isomers with enzymes leading to rapid degradation that need further validation by cloning and biochemical experiments.

Genomic analysis reveals the population structure and antimicrobial resistance of avian Pasteurella multocida in China

OBJECTIVES

To investigate the population structure and antimicrobial resistance (AMR) of avian Pasteurella multocida in China.

METHODS

Utilizing WGS analysis, we explored the phylogeny using a dataset of 546 genomes, comprising avian P. multocida isolates from China (n = 121), the USA (n = 165), Australia(n = 153), Bangladesh (n = 3) and isolates of other hosts from China (n = 104). We examined the integrative and conjugative element (ICE) structures and the distribution of their components carrying resistance genes, and reconstructed the evolutionary history of A:L1:ST129 (n = 110).

RESULTS

The population structure of avian P. multocida in China was dominated by the A:L1:ST129 clone with limited genetic diversity. A:L1:ST129 isolates possessed a broader spectrum of resistance genes at comparatively higher frequencies than those from other hosts and countries. The novel putative ICEs harboured complex resistant clusters that were prevalent in A:L1:ST129. Bayesian analysis predicted that the A:L1:ST129 clone emerged around 1923, and evolved slowly.

CONCLUSIONS

A:L1:ST129 appears to possess a host predilection towards avian species in China, posing a potential health threat to other animals. The complex AMR determinants coupled with high frequencies may strengthen the population dominance of A:L1:ST129. The extensive antimicrobial utilization in poultry farming and the mixed rearing practices could have accelerated AMR accumulation in A:L1:ST129. ICEs, together with their resistant clusters, significantly contribute to resistance gene transfer and facilitate the adaptation of A:L1:ST129 to ecological niches. Despite the genetic stability and slow evolution rate, A:L1:ST129 deserves continued monitoring due to its propensity to retain resistance genes, warranting global attention to preclude substantial economic losses.

Mobilome-driven partitions of the resistome in Salmonella

IMPORTANCE

Antimicrobial resistance (AMR) has become a significant global challenge, with an estimated 10 million deaths annually by 2050. The emergence of AMR is mainly attributed to mobile genetic elements (MGEs or mobilomes), which accelerate wide dissemination among pathogens. The interaction between mobilomes and AMR genes (or resistomes) in Salmonella, a primary cause of diarrheal diseases that results in over 90 million cases annually, remains poorly understood. The available fragmented or incomplete genomes remain a significant limitation in investigating the relationship between AMR and MGEs. Here, we collected the most extensive closed Salmonella genomes (n = 1,817) from various sources across 58 countries. Notably, our results demonstrate that resistome transmission between Salmonella lineages follows a specific pattern of MGEs and is influenced by external drivers, including certain socioeconomic factors. Therefore, targeted interventions are urgently needed to mitigate the catastrophic consequences of Salmonella AMR.

Complex interactions between diverse mobile genetic elements drive the evolution of metal-resistant bacterial genomes

In this study, we compared the genomes of three metal-resistant bacteria isolated from mercury-contaminated soil. We identified diverse and novel MGEs with evidence of multiple LGT events shaping their genomic structure and heavy metal resistance. Among the three metal-resistant strains, Sphingobium sp SA2 and Sphingopyxis sp SE2 were resistant to multiple metals including mercury, cadmium, copper, zinc and lead. Pseudoxanthomonas sp SE1 showed resistance to mercury only. Whole genome sequencing by Illumina and Oxford Nanopore technologies was undertaken to obtain comprehensive genomic data. The Sphingobium and Sphingopyxis strains contained multiple chromosomes and plasmids, whereas the Pseudoxanthomonas strain contained one circular chromosome. Consistent with their metal resistance profiles, the strains of Sphingobium and Sphingopyxis contained a higher quantity of diverse metal resistance genes across their chromosomes and plasmids compared to the single-metal resistant Pseudoxanthomonas SE1. In all three strains, metal resistance genes were principally associated with various novel MGEs including genomic islands (GIs), integrative conjugative elements (ICEs), transposons, insertion sequences (IS), recombinase in trio (RIT) elements and group II introns, indicating their importance in facilitating metal resistance adaptation in a contaminated environment. In the Pseudoxanthomonas strain, metal resistance regions were largely situated on a GI. The chromosomes of the strains of Sphingobium and Sphingopyxis contained multiple metal resistance regions, which were likely acquired by several GIs, ICEs, numerous IS elements, several Tn3 family transposons and RIT elements. Two of the plasmids of Sphingobium were impacted by Tn3 family transposons and ISs likely integrating metal resistance genes. The two plasmids of Sphingopyxis harboured transposons, IS elements, an RIT element and a group II intron. This study provides a comprehensive annotation of complex genomic regions of metal resistance associated with novel MGEs. It highlights the critical importance of LGT in the evolution of metal resistance of bacteria in contaminated environments.

Comprehensive probiogenomics analysis of the commensal Escherichia coli CEC15 as a potential probiotic strain

BACKGROUND

Probiotics have gained attention for their potential maintaining gut and immune homeostasis. They have been found to confer protection against pathogen colonization, possess immunomodulatory effects, enhance gut barrier functionality, and mitigate inflammation. However, a thorough understanding of the unique mechanisms of effects triggered by individual strains is necessary to optimize their therapeutic efficacy. Probiogenomics, involving high-throughput techniques, can help identify uncharacterized strains and aid in the rational selection of new probiotics. This study evaluates the potential of the Escherichia coli CEC15 strain as a probiotic through in silico, in vitro, and in vivo analyses, comparing it to the well-known probiotic reference E. coli Nissle 1917. Genomic analysis was conducted to identify traits with potential beneficial activity and to assess the safety of each strain (genomic islands, bacteriocin production, antibiotic resistance, production of proteins involved in host homeostasis, and proteins with adhesive properties). In vitro studies assessed survival in gastrointestinal simulated conditions and adhesion to cultured human intestinal cells. Safety was evaluated in BALB/c mice, monitoring the impact of E. coli consumption on clinical signs, intestinal architecture, intestinal permeability, and fecal microbiota. Additionally, the protective effects of both strains were assessed in a murine model of 5-FU-induced mucositis.

RESULTS

CEC15 mitigates inflammation, reinforces intestinal barrier, and modulates intestinal microbiota. In silico analysis revealed fewer pathogenicity-related traits in CEC15, when compared to Nissle 1917, with fewer toxin-associated genes and no gene suggesting the production of colibactin (a genotoxic agent). Most predicted antibiotic-resistance genes were neither associated with actual resistance, nor with transposable elements. The genome of CEC15 strain encodes proteins related to stress tolerance and to adhesion, in line with its better survival during digestion and higher adhesion to intestinal cells, when compared to Nissle 1917. Moreover, CEC15 exhibited beneficial effects on mice and their intestinal microbiota, both in healthy animals and against 5FU-induced intestinal mucositis.

CONCLUSIONS

These findings suggest that the CEC15 strain holds promise as a probiotic, as it could modulate the intestinal microbiota, providing immunomodulatory and anti-inflammatory effects, and reinforcing the intestinal barrier. These findings may have implications for the treatment of gastrointestinal disorders, particularly some forms of diarrhea.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Complete genome sequence and comparative analysis of a Vibrio vulnificus strain isolated from a clinical patient

Vibrio vulnificus is an opportunistic, global pathogen that naturally inhabits sea water and is responsible for most vibriosis-related deaths. We investigated the genetic characteristics of V. vulnificus isolated from the clinical blood culture specimen of a patient with hepatitis B virus cirrhosis in 2018 (named as V. vulnificus VV2018) by whole genome sequencing (WGS). VV2018 belonged to a novel sequencing type 620 (ST620) and comprised two circular chromosomes, containing 4,389 potential coding sequences (CDSs) and 152 RNA genes. The phylogenetic tree of single nucleotide polymorphisms (SNPs) using 26 representative genomes revealed that VV2108 grouped with two other V. vulnificus strains isolated from humans. The pan-genome of V. vulnificus was constructed using 26 representative genomes to elucidate their genetic diversity, evolutionary characteristics, and virulence and antibiotic resistance profiles. The pan-genome analysis revealed that VV2018 shared a total of 3,016 core genes (≥99% presence), including 115 core virulence factors (VFs) and 5 core antibiotic resistance-related genes, and 309 soft core genes (≥95 and <99% presence) with 25 other V. vulnificus strains. The varG gene might account for the cefazolin resistance, and comparative analysis of the genetic context of varG revealed that two genes upstream and downstream of varG were conserved. The glycosylation (pgl) like genes were found in VV2018 compared with Pgl-related proteins in Neisseria that might affect the adherence of the strain in hosts. The comparative analysis of VV2018 would contribute to a better understanding of the virulence and antibiotic resistance profiles of V. vulnificus. Meanwhile much work remains to be done to better understand the function of pgl-like genes in V. vulnificus.

Recursive genome engineering decodes the evolutionary origin of an essential thymidylate kinase activity in Pseudomonas putida KT2440

IMPORTANCE

Investigating fundamental aspects of metabolism is vital for advancing our understanding of the diverse biochemical capabilities and biotechnological applications of bacteria. The origin of the essential thymidylate kinase function in the model bacterium Pseudomonas putida KT2440, seemingly interrupted due to the presence of a large genomic island that disrupts the cognate gene, eluded a satisfactory explanation thus far. This is a first-case example of an essential metabolic function, likely acquired by horizontal gene transfer, which "landed" in a locus encoding the same activity. As such, foreign DNA encoding an essential dNMPK could immediately adjust to the recipient host-instead of long-term accommodation and adaptation. Understanding how these functions evolve is a major biological question, and the work presented here is a decisive step toward this direction. Furthermore, identifying essential and accessory genes facilitates removing those deemed irrelevant in industrial settings-yielding genome-reduced cell factories with enhanced properties and genetic stability.

The Genomic-Driven Discovery of Glutarimide-Containing Derivatives from Burkholderia gladioli

Glutarimide-containing polyketides exhibiting potent antitumor and antimicrobial activities were encoded via conserved module blocks in various strains that favor the genomic mining of these family compounds. The bioinformatic analysis of the genome of Burkholderia gladioli ATCC 10248 showed a silent trans-AT PKS biosynthetic gene cluster (BGC) on chromosome 2 (Chr2C8), which was predicted to produce new glutarimide-containing derivatives. Then, the silent polyketide synthase gene cluster was successfully activated via in situ promoter insertion and heterologous expression. As a result, seven glutarimide-containing analogs, including five new ones, gladiofungins D-H (3-7), and two known gladiofungin A/gladiostatin (1) and 2 (named gladiofungin C), were isolated from the fermentation of the activated mutant. Their structures were elucidated through the analysis of HR-ESI-MS and NMR spectroscopy. The structural diversities of gladiofungins may be due to the degradation of the butenolide group in gladiofungin A (1) during the fermentation and extraction process. Bioactivity screening showed that 2 and 4 had moderate anti-inflammatory activities. Thus, genome mining combined with promoter engineering and heterologous expression were proved to be effective strategies for the pathway-specific activation of the silent BGCs for the directional discovery of new natural products.

That's no moon, it's a Starship: Giant transposons driving fungal horizontal gene transfer

To date, most reports of horizontal gene transfer (HGT) in fungi rely on genome sequence data and are therefore an indirect measure of HGT after the event has occurred. However, a novel group of class II-like transposons known as Starships may soon alter this status quo. Starships are giant transposable elements that carry dozens of genes, some of which are host-beneficial, and are linked to many recent HGT events in the fungal kingdom. These transposons remain active and mobile in many fungal genomes and their transposition has recently been shown to be driven by a conserved tyrosine-recombinase called 'Captain'. This perspective explores some of the remaining unanswered questions about how these Starship transposons move, both within a genome and between different species. We seek to outline several experimental approaches that can be used to identify the genes essential for Starship-mediated HGT and draw links to other recently discovered giant transposons outside of the fungal kingdom.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

Multi-omic profiling of a novel activated sludge strain Sphingobacterium sp. WM1 reveals the mechanism of tetracycline biodegradation and its merits of potential application

As an emerging pollutant, the antibiotic tetracycline (TC) has been consistently detected in wastewater and activated sludge. Biodegradation represents a potentially crucial pathway to dissipate TC contamination. However, few efficient TC-degrading bacteria have been isolated and a comprehensive understanding of the molecular mechanisms underlying TC degradation is still lacking. In this study, a novel TC-degrading bacterium, designated as Sphingobacterium sp. WM1, was successfully isolated from activated sludge. Strain WM1 exhibited a remarkable performance in degrading 50 mg/L TC within 1 day under co-metabolic conditions. Genomic analysis of the strain WM1 unveiled the presence of three functional tetX genes. Unraveling the complex molecular mechanisms, transcriptome analysis highlighted the role of upregulated transmembrane transport and accelerated electron transport in facilitating TC degradation. Proteomics confirmed the up-regulation of proteins involved in cellular biosynthesis/metabolism and ribosomal processes. Crucially, the tetX gene-encoding protein showed a significant upregulation, indicating its role in TC degradation. Heterologous expression of the tetX gene resulted in TC dissipation from an initial 51.9 mg/L to 4.2 mg/L within 24 h. The degradation pathway encompassed TC hydroxylation, transforming into TP461 and subsequent metabolites, which effectively depleted TC's inhibitory activity. Notably, the tetX genes in strain WM1 showed limited potential for horizontal gene transfer. Collectively, strain WM1's potent TC degradation capacity signals a promise for enhancing TC clean-up strategies.

Genomic and pathogenicity islands of Listeria monocytogenes-overview of selected aspects

Listeria monocytogenes causes listeriosis, a disease characterized by a high mortality rate (up to 30%). Since the pathogen is highly tolerant to changing conditions (high and low temperature, wide pH range, low availability of nutrients), it is widespread in the environment, e.g., water, soil, or food. L. monocytogenes possess a number of genes that determine its high virulence potential, i.e., genes involved in the intracellular cycle (e.g., prfA, hly, plcA, plcB, inlA, inlB), response to stress conditions (e.g., sigB, gadA, caspD, clpB, lmo1138), biofilm formation (e.g., agr, luxS), or resistance to disinfectants (e.g., emrELm, bcrABC, mdrL). Some genes are organized into genomic and pathogenicity islands. The islands LIPI-1 and LIPI-3 contain genes related to the infectious life cycle and survival in the food processing environment, while LGI-1 and LGI-2 potentially ensure survival and durability in the production environment. Researchers constantly have been searching for new genes determining the virulence of L. monocytogenes. Understanding the virulence potential of L. monocytogenes is an important element of public health protection, as highly pathogenic strains may be associated with outbreaks and the severity of listeriosis. This review summarizes the selected aspects of L. monocytogenes genomic and pathogenicity islands, and the importance of whole genome sequencing for epidemiological purposes.

Comparative genomics of two Vietnamese Helicobacter pylori strains, CHC155 from a non-cardia gastric cancer patient and VN1291 from a duodenal ulcer patient

Helicobacter pylori is involved in the etiology and severity of several gastroduodenal diseases; however, plasticity of the H. pylori genome makes complete genome assembly difficult. We report here the full genomes of H. pylori strains CHC155 and VN1291 isolated from a non-cardia gastric cancer patient and a duodenal ulcer patient, respectively, and their virulence demonstrated by in vitro infection. Whole-genome sequences were obtained by combining long- and short-reads with a hybrid-assembly approach. Both CHC155 and VN1291 genome possessed four kinds of genomic island: a cag pathogenicity island (cagPAI), two type 4 secretion system islands within an integrative and conjugative element (tfs ICE), and prophage. CHC155 and VN1291 carried East Asian-type cagA and vacA s1m1, and outer membrane protein genes, including two copies of oipA. Corresponded to genetic determinants of antibiotic resistance, chromosomal mutations were identified in CHC155 (rdxA, gyrA, and 23S rRNA) and VN1291 (rdxA, 23S rRNA, and pbp1A). In vitro infection of AGS cells by both strains induced the cell scattering phenotype, tyrosine phosphorylation of CagA, and promoted high levels of IL8 secretion, indicating fully intact phenotypes of the cagPAI. Virulence genes in CHC155 and VN1291 genomes are crucial for H. pylori pathogenesis and are risk factors in the development of gastric cancer and duodenal ulcer. Our in vitro studies indicate that the strains CHC155 and VN1291 carry the pathogenic potential.

Cascade capture, oxidization and inactivation for removing multi-species pollutants, antimicrobial resistance and pathogenicity from hospital wastewater

As reservoirs of pathogens, antimicrobial resistant microorganisms and a wide variety of pollutants, hospital wastewaters (HWWs) need to be effectively treated before discharge. This study employed the functionalized colloidal microbubble technology as one-step fast HWW treatment. Inorganic coagulant (monomeric Fe(III)-coagulant or polymeric Al(III)-coagulant) and ozone were used as surface-decorator and gaseous core modifier, respectively. The Fe(III)- or Al(III)-modified colloidal gas (or, ozone) microbubbles (Fe(III)-CCGMBs, Fe(III)-CCOMBs, Al(III)-CCGMBs and Al(III)-CCOMBs) were constructed. Within 3 min, CCOMBs decreased COD_Cr and fecal coliform concentration to the levels meeting the national discharge standard for medical organization. Regrowth of bacteria was inhibited and biodegradability of organics was increased after the simultaneous oxidation and cell-inactivation process. The metagenomics analysis further reveals that Al(III)-CCOMBs performed best in capturing the virulence genes, antibiotic resistance genes and their potential hosts. The horizontal transfer of those harmful genes could be effectively hampered thanks to the removal of mobile genetic elements. Interestingly, the virulence factors of adherence, micronutrient uptake/acquisition and phase invasion could facilitate the interface-dominated capture. Featured as cascade processes of capture, oxidation and inactivation in the one-step operation, the robust Al(III)-CCOMB treatment is recommended for the HWW treatment and the protection of downstream aquatic environment.

Potential mobile units drive the horizontal transfer of phytoplasma effector phyllogen genes

Phytoplasmas are obligate intracellular plant pathogenic bacteria that can induce phyllody, which is a type of abnormal floral organ development. Phytoplasmas possess phyllogens, which are effector proteins that cause phyllody in plants. Phylogenetic comparisons of phyllogen and 16S rRNA genes have suggested that phyllogen genes undergo horizontal transfer between phytoplasma species and strains. However, the mechanisms and evolutionary implications of this horizontal gene transfer are unclear. Here, we analyzed synteny in phyllogen flanking genomic regions from 17 phytoplasma strains that were related to six 'Candidatus' species, including three strains newly sequenced in this study. Many of the phyllogens were flanked by multicopy genes within potential mobile units (PMUs), which are putative transposable elements found in phytoplasmas. The multicopy genes exhibited two distinct patterns of synteny that correlated with the linked phyllogens. The low level of sequence identities and partial truncations found among these phyllogen flanking genes indicate that the PMU sequences are deteriorating, whereas the highly conserved sequences and functions (e.g., inducing phyllody) of the phyllogens suggest that the latter are important for phytoplasma fitness. Furthermore, although their phyllogens were similar, PMUs in strains related to 'Ca. P. asteris' were often located in different regions of the genome. These findings strongly indicate that PMUs drive the horizontal transfer of phyllogens among phytoplasma species and strains. These insights improve our understanding of how symptom-determinant genes have been shared among phytoplasmas.

Comprehensive genomic analysis of hypocholesterolemic probiotic Enterococcus faecium LR13 reveals unique proteins involved in cholesterol-assimilation

Hypercholesterolemia is a major risk factor for cardiovascular diseases (CVDs). Chemotherapeutic agents for CVDs exhibit several side effects. Specific probiotics with hypocholesterolemic effects can be safe and effective alternatives to chemotherapeutics. Here, we have analyzed and compared the genome of a novel rhizospheric Enterococcus faecium LR13 cholesterol-assimilating probiotic with other probiotic/pathogenic E. faecium strains to discern genetic factors underlying probiotic efficacy and cholesterol-assimilation. Genomic analyses of E. faecium probiotic strains revealed that LR13 and WEFA23 (cholesterol-assimilating probiotics) harbored 21 unique proteins absent in non-cholesterol-assimilating probiotics. Of these, 14 proteins could directly help in cholesterol-assimilation by producing short chain fatty acids, lipid (sterol) transport and membrane stabilization, and bile salt hydrolase activity. This suggests that cholesterol-assimilation is an intrinsic, strain-specific trait exhibited by probiotics with a specific genetic constitution. Moreover, the unique proteins identified in this study can serve as biomarkers for discerning/characterizing cholesterol-assimilating probiotics as novel biotherapeutics against CVDs.

Silent gene clusters encode magnetic organelle biosynthesis in a non-magnetotactic phototrophic bacterium

Horizontal gene transfer is a powerful source of innovations in prokaryotes that can affect almost any cellular system, including microbial organelles. The formation of magnetosomes, one of the most sophisticated microbial mineral-containing organelles synthesized by magnetotactic bacteria for magnetic navigation in the environment, was also shown to be a horizontally transferrable trait. However, the mechanisms determining the fate of such genes in new hosts are not well understood, since non-adaptive gene acquisitions are typically rapidly lost and become unavailable for observation. This likely explains why gene clusters encoding magnetosome biosynthesis have never been observed in non-magnetotactic bacteria. Here, we report the first discovery of a horizontally inherited dormant gene clusters encoding biosynthesis of magnetosomes in a non-magnetotactic phototrophic bacterium Rhodovastum atsumiense. We show that these clusters were inactivated through transcriptional silencing and antisense RNA regulation, but retain functionality, as several genes were able to complement the orthologous deletions in a remotely related magnetotactic bacterium. The laboratory transfer of foreign magnetosome genes to R. atsumiense was found to endow the strain with magnetosome biosynthesis, but strong negative selection led to rapid loss of this trait upon subcultivation, highlighting the trait instability in this organism. Our results provide insight into the horizontal dissemination of gene clusters encoding complex prokaryotic organelles and illuminate the potential mechanisms of their genomic preservation in a dormant state.

Vibrio parahaemolyticus Epidemiology and Pathogenesis: Novel Insights on an Emerging Foodborne Pathogen

The epidemiological dynamics of V. parahaemolyticus´ infections have been characterized by the abrupt appearance of outbreaks in remote areas where these diseases had not been previously detected, without knowing the routes of entry of the pathogens in the new area. However, there are recent studies that show the link between the appearance of epidemic outbreaks of Vibrio and environmental factors such as oceanic transport of warm waters, which has provided a possible mechanism for the dispersion of Vibrio diseases globally. Despite this evidence, there is little information on the possible routes of entry and transport of infectious agents from endemic countries to the entire world. In this sense, the recent advances in genomic sequencing tools are making it possible to infer possible biogeographical patterns of diverse pathogens with relevance in public health like V. parahaemolyticus. In this chapter, we will address several general aspects about V. parahaemolyticus, including their microbiological and genetic detection, main virulence factors, and the epidemiology of genotypes involved in foodborne outbreaks globally.

Genomic Island-Encoded Histidine Kinase and Response Regulator Coordinate Mannose Utilization with Virulence in Enterohemorrhagic Escherichia coli

Enterohemorrhagic Escherichia coli (EHEC) is a highly adaptive pathogen and has acquired diverse genetic elements, such as genomic islands and prophages, via horizontal gene transfer to promote fitness in vivo. Two-component signaling systems (TCSs) allow bacteria to sense, respond to, and adapt to various environments. This study identified a putative two-component signaling system composed of the histidine kinase EDL5436 (renamed LmvK) and the response regulator EDL5428 (renamed LmvR) in EHEC. lmvK and lmvR along with EDL5429 to EDL5434 (EDL5429-5434) between them constitute the OI167 genomic island and are highly associated with the EHEC pathotype. EDL5429-5434 encode transporters and metabolic enzymes that contribute to growth on mannose and are directly upregulated by LmvK/LmvR in the presence of mannose, as revealed by quantitative PCR (qPCR) and DNase I footprint assays. Moreover, LmvR directly activates the expression of the type III secretion system in response to mannose and promotes the formation of attaching and effacing lesions on HeLa cells. Using human colonoid and mouse infection models, we show that lmvK and lmvR contributed greatly to adherence and microcolony (MC) formation ex vivo and colonization in vivo. Finally, RNA sequencing and chromatin immunoprecipitation coupled with sequencing analyses identified additional direct targets of LmvR, most of which are involved in metabolism. Given that mannose is a mucus-derived sugar that induces virulence and is preferentially used by EHEC during infection, our data revealed a previously unknown mechanism by which EHEC recognizes the host metabolic landscape and regulates virulence expression accordingly. Our findings provide insights into how pathogenic bacteria evolve by acquiring genetic elements horizontally to adapt to host environments. IMPORTANCE The gastrointestinal tract represents a complex and challenging environment for enterohemorrhagic Escherichia coli (EHEC). However, EHEC is a highly adaptable pathogen, requiring only 10 to 100 CFUs to cause infection. This ability was achieved partially by acquiring mobile genetic elements, such as genomic islands, that promote overall fitness. Mannose is an intestinal mucus-derived sugar that stimulates virulence and is preferentially used by EHEC during infection. Here, we characterize the OI167 genomic island of EHEC, which encodes a novel two-component signaling system (TCS) and transporters and metabolic enzymes (EDL5429-5434) involved in mannose utilization. The TCS directly upregulates EDL5429-5434 and genes encoding the type III secretion system in the presence of mannose. Moreover, the TCS contributes greatly to EHEC virulence ex vivo and in vivo. Our data demonstrate an elegant example in which EHEC strains evolve by acquiring genetic elements horizontally to recognize the host metabolic landscape and regulate virulence expression accordingly, leading to successful infections.

Insights into the genome of Methylobacterium sp. NMS14P, a novel bacterium for growth promotion of maize, chili, and sugarcane

A novel methylotrophic bacterium designated as NMS14P was isolated from the root of an organic coffee plant (Coffea arabica) in Thailand. The 16S rRNA sequence analysis revealed that this new isolate belongs to the genus Methylobacterium, and its novelty was clarified by genomic and comparative genomic analyses, in which NMS14P exhibited low levels of relatedness with other Methylobacterium-type strains. NMS14P genome consists of a 6,268,579 bp chromosome, accompanied by a 542,519 bp megaplasmid and a 66,590 bp plasmid, namely pNMS14P1 and pNMS14P2, respectively. Several genes conferring plant growth promotion are aggregated on both chromosome and plasmids, including phosphate solubilization, indole-3-acetic acid (IAA) biosynthesis, cytokinins (CKs) production, 1-aminocyclopropane-1-carboxylate (ACC) deaminase activity, sulfur-oxidizing activity, trehalose synthesis, and urea metabolism. Furthermore, pangenome analysis showed that NMS14P possessed the highest number of strain-specific genes accounting for 1408 genes, particularly those that are essential for colonization and survival in a wide array of host environments, such as ABC transporter, chemotaxis, quorum sensing, biofilm formation, and biosynthesis of secondary metabolites. In vivo tests have supported that NMS14P significantly promoted the growth and development of maize, chili, and sugarcane. Collectively, NMS14P is proposed as a novel plant growth-promoting Methylobacterium that could potentially be applied to a broad range of host plants as Methylobacterium-based biofertilizers to reduce and ultimately substitute the use of synthetic agrochemicals for sustainable agriculture.

Transmission of antibiotic resistance genes through mobile genetic elements in Acinetobacter baumannii and gene-transfer prevention

Antibiotic resistance is a major global public health concern. Acinetobacter baumannii is a nosocomial pathogen that has emerged as a global threat because of its high levels of resistance to many antibiotics, particularly those considered as last-resort antibiotics, such as carbapenems. Mobile genetic elements (MGEs) play an important role in the dissemination and expression of antibiotic resistance genes (ARGs), including the mobilization of ARGs within and between species. We conducted an in-depth, systematic investigation of the occurrence and dissemination of ARGs associated with MGEs in A. baumannii. We focused on a cross-sectoral approach that integrates humans, animals, and environments. Four strategies for the prevention of ARG dissemination through MGEs have been discussed: prevention of airborne transmission of ARGs using semi-permeable membrane-covered thermophilic composting; application of nanomaterials for the removal of emerging pollutants (antibiotics) and pathogens; tertiary treatment technologies for controlling ARGs and MGEs in wastewater treatment plants; and the removal of ARGs by advanced oxidation techniques. This review contemplates and evaluates the major drivers involved in the transmission of ARGs from the cross-sectoral perspective and ARG-transfer prevention processes.

Genomic diversity, pathogenicity and antimicrobial resistance of Escherichia coli isolated from poultry in the southern United States

Escherichia coli (E. coli) are typically present as commensal bacteria in the gastro-intestinal tract of most animals including poultry species, but some avian pathogenic E. coli (APEC) strains can cause localized and even systematic infections in domestic poultry. Emergence and re-emergence of antimicrobial resistant isolates (AMR) constrain antibiotics usage in poultry production, and development of an effective vaccination program remains one of the primary options in E. coli disease prevention and control for domestic poultry. Thus, understanding genetic and pathogenic diversity of the enzootic E. coli isolates, particularly APEC, in poultry farms is the key to designing an optimal vaccine candidate and to developing an effective vaccination program. This study explored the genomic and pathogenic diversity among E. coli isolates in southern United States poultry. A total of nine isolates were recovered from sick broilers from Mississippi, and one from Georgia, with epidemiological variations among clinical signs, type of housing, and bird age. The genomes of these isolates were sequenced by using both Illumina short-reads and Oxford Nanopore long-reads, and our comparative analyses suggested data from both platforms were highly consistent. The 16 s rRNA based phylogenetic analyses showed that the 10 bacteria strains are genetically closer to each other than those in the public database. However, whole genome analyses showed that these 10 isolates encoded a diverse set of reported virulence and AMR genes, belonging to at least nine O:H serotypes, and are genetically clustered with at least five different groups of E. coli isolates reported by other states in the United States. Despite the small sample size, this study suggested that there was a large extent of genomic and serological diversity among E. coli isolates in southern United States poultry. A large-scale comprehensive study is needed to understand the overall genomic diversity and the associated virulence, and such a study will be important to develop a broadly protective E. coli vaccine.

Genomic features, antimicrobial susceptibility, and epidemiological insights into Burkholderia cenocepacia clonal complex 31 isolates from bloodstream infections in India

Introduction

Burkholderia cepacia complex (Bcc) clonal complex (CC) 31, the predominant lineage causing devastating outbreaks globally, has been a growing concern of infections in non-cystic fibrosis (NCF) patients in India. B. cenocepacia is very challenging to treat owing to its virulence determinants and antibiotic resistance. Improving the management of these infections requires a better knowledge of their resistance patterns and mechanisms.

Methods

Whole-genome sequences of 35 CC31 isolates obtained from patient samples, were analyzed against available 210 CC31 genomes in the NCBI database to glean details of resistance, virulence, mobile elements, and phylogenetic markers to study genomic diversity and evolution of CC31 lineage in India.

Results

Genomic analysis revealed that 35 isolates belonging to CC31 were categorized into 11 sequence types (ST), of which five STs were reported exclusively from India. Phylogenetic analysis classified 245 CC31 isolates into eight distinct clades (I-VIII) and unveiled that NCF isolates are evolving independently from the global cystic fibrosis (CF) isolates forming a distinct clade. The detection rate of seven classes of antibiotic-related genes in 35 isolates was 35 (100%) for tetracyclines, aminoglycosides, and fluoroquinolones; 26 (74.2%) for sulphonamides and phenicols; 7 (20%) for beta-lactamases; and 1 (2.8%) for trimethoprim resistance genes. Additionally, 3 (8.5%) NCF isolates were resistant to disinfecting agents and antiseptics. Antimicrobial susceptibility testing revealed that majority of NCF isolates were resistant to chloramphenicol (77%) and levofloxacin (34%). NCF isolates have a comparable number of virulence genes to CF isolates. A well-studied pathogenicity island of B. cenocepacia, GI11 is present in ST628 and ST709 isolates from the Indian Bcc population. In contrast, genomic island GI15 (highly similar to the island found in B. pseudomallei strain EY1) is exclusively reported in ST839 and ST824 isolates from two different locations in India. Horizontal acquisition of lytic phage ST79 of pathogenic B. pseudomallei is demonstrated in ST628 isolates Bcc1463, Bcc29163, and BccR4654 amongst CC31 lineage.

Discussion

The study reveals a high diversity of CC31 lineages among B. cenocepacia isolates from India. The extensive information from this study will facilitate the development of rapid diagnostic and novel therapeutic approaches to manage B. cenocepacia infections.

SourceFinder: a Machine-Learning-Based Tool for Identification of Chromosomal, Plasmid, and Bacteriophage Sequences from Assemblies

High-throughput genome sequencing technologies enable the investigation of complex genetic interactions, including the horizontal gene transfer of plasmids and bacteriophages. However, identifying these elements from assembled reads remains challenging due to genome sequence plasticity and the difficulty in assembling complete sequences. In this study, we developed a classifier, using random forest, to identify whether sequences originated from bacterial chromosomes, plasmids, or bacteriophages. The classifier was trained on a diverse collection of 23,211 chromosomal, plasmid, and bacteriophage sequences from hundreds of bacterial species. In order to adapt the classifier to incomplete sequences, each complete sequence was subsampled into 5,000 nucleotide fragments and further subdivided into k-mers. This three-class classifier succeeded in identifying chromosomes, plasmids, and bacteriophages using k-mer distributions of complete and partial genome sequences, including simulated metagenomic scaffolds with minimum performance of 0.939 area under the receiver operating characteristic curve (AUC). This classifier, implemented as SourceFinder, has been made available as an online web service to help the community with predicting the chromosomal, plasmid, and bacteriophage sources of assembled bacterial sequence data (https://cge.food.dtu.dk/services/SourceFinder/). IMPORTANCE Extra-chromosomal genes encoding antimicrobial resistance, metal resistance, and virulence provide selective advantages for bacterial survival under stress conditions and pose serious threats to human and animal health. These accessory genes can impact the composition of microbiomes by providing selective advantages to their hosts. Accurately identifying extra-chromosomal elements in genome sequence data are critical for understanding gene dissemination trajectories and taking preventative measures. Therefore, in this study, we developed a random forest classifier for identifying the source of bacterial chromosomal, plasmid, and bacteriophage sequences.

Genome-Based Analysis of Virulence Factors and Biofilm Formation in Novel P. aeruginosa Strains Isolated from Household Appliances

In household washing machines, opportunistic pathogens such as Pseudomonas aeruginosa are present, which represent the household as a possible reservoir for clinical pathogens. Here, four novel P. aeruginosa strains, isolated from different sites of household appliances, were investigated regarding their biofilm formation. Only two isolates showed strong surface-adhered biofilm formation. In consequence of these phenotypic differences, we performed whole genome sequencing using Oxford Nanopore Technology together with Illumina MiSeq. Whole genome data were screened for the prevalence of 285 virulence- and biofilm-associated genes as well as for prophages. Linking biofilm phenotypes and parallelly appearing gene compositions, we assume a relevancy of the las quorum sensing system and the phage-encoded bacteriophage control infection gene bci, which was found on integrated phi297 DNA in all biofilm-forming isolates. Additionally, only the isolates revealing strong biofilm formation harbored the ϕCTX-like prophage Dobby, implicating a role of this prophage on biofilm formation. Investigations on clinically relevant pathogens within household appliances emphasize their adaptability to harsh environments, with high concentrations of detergents, providing greater insights into pathogenicity and underlying mechanisms. This in turn opens the possibility to map and characterize potentially relevant strains even before they appear as pathogens in society.

Fluxes of the Amazon River plume nutrients and microbes into marine sponges

Sponges have co-evolved with microbes for over 400 myr. Previous studies have demonstrated that sponges can be classified according to the abundance of microbes in their tissues as Low Microbial Abundance (LMA) and High Microbial Abundance (HMA). While LMA sponges rely mainly on water column microbes, HMA appear to rely much more on symbiotic fermentative and autotrophic microbes maintained in their tissues. However, it is unclear if this pattern holds when comparing different species of tropical sponges under extreme nutrient conditions and sediment loads in the water column, such as the Great Amazon Reef System (GARS), which covers an area of ~56,000 km² off the Amazon River mouth. Sponges are the major GARS benthic components. However, these sponges' microbiome across the GARS is still unknown. Here, we investigated water quality, isotopic values (δ¹³C and δ¹⁵N), metagenomic and lipidomic profiles of sponges obtained from different sectors throughout the GARS. >180 million shotgun metagenomic reads were annotated, covering 22 sponge species. Isotopic and lipidomic analyses suggested LMA sponges rely on the Amazon River Plume for nutrition. HMA sponges (N = 15) had higher Roseiflexus and Nitrospira abundance, whereas LMA sponges (N = 7) had higher Prochlorococcus and Pelagibacter abundance. Functional data revealed that the LMA sponge microbiomes had greater number of sequences related to phages and prophages as well as electron transport and photophosphorylation which may be related to photosynthetic processes associated with the Prochlorococcus and Synechococcus found in the LMA. The higher phages abundance in LMA sponges could be related to these holobionts' reduced defense towards phage infection. Meanwhile, HMA sponge microbiomes had higher Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR abundance, which may be involved in defense against phage infection. This study sheds light on the nutrient fluxes and microbes from the Amazon River plume into the sponge holobionts.

The Origin, Function, Distribution, Quantification, and Research Advances of Extracellular DNA

In nature, DNA is ubiquitous, existing not only inside but also outside of the cells of organisms. Intracellular DNA (iDNA) plays an essential role in different stages of biological growth, and it is defined as the carrier of genetic information. In addition, extracellular DNA (eDNA) is not enclosed in living cells, accounting for a large proportion of total DNA in the environment. Both the lysis-dependent and lysis-independent pathways are involved in eDNA release, and the released DNA has diverse environmental functions. This review provides an insight into the origin as well as the multiple ecological functions of eDNA. Furthermore, the main research advancements of eDNA in the various ecological environments and the various model microorganisms are summarized. Furthermore, the major methods for eDNA extraction and quantification are evaluated.

Comparative genomic analyses reveal genetic characteristics and pathogenic factors of Bacillus pumilus HM-7

Bacillus pumilus plays an important role in industrial application and biocontrol activities, as well as causing humans and plants disease, leading to economic losses and biosafety concerns. However, until now, the pathogenesis and underlying mechanisms of B. pumilus strains remain unclear. In our previous study, one representative isolate of B. pumilus named HM-7 has been recovered and proved to be the causal agent of fruit rot on muskmelon (Cucumis melo). Herein, we present a complete and annotated genome sequence of HM-7 that contains 4,111 coding genes in a single 3,951,520 bp chromosome with 41.04% GC content. A total of 3,481 genes were functionally annotated with the GO, COG, and KEGG databases. Pan-core genome analysis of HM-7 and 20 representative B. pumilus strains, as well as six closely related Bacillus species, discovered 740 core genes and 15,205 genes in the pan-genome of 21 B. pumilus strains, in which 485 specific-genes were identified in HM-7 genome. The average nucleotide identity (ANI), and whole-genome-based phylogenetic analysis revealed that HM-7 was most closely related to the C4, GR8, MTCC-B6033, TUAT1 and SH-B11 strains, but evolutionarily distinct from other strains in B. pumilus. Collinearity analysis of the six similar B. pumilus strains showed high levels of synteny but also several divergent regions for each strains. In the HM-7 genome, we identified 484 genes in the carbohydrate-active enzymes (CAZyme) class, 650 genes encoding virulence factors, and 1,115 genes associated with pathogen-host interactions. Moreover, three HM-7-specific regions were determined, which contained 424 protein-coding genes. Further investigation of these genes showed that 19 pathogenesis-related genes were mainly associated with flagella formation and secretion of toxic products, which might be involved in the virulence of strain HM-7. Our results provided detailed genomic and taxonomic information for the HM-7 strain, and discovered its potential pathogenic mechanism, which lay a foundation for developing effective prevention and control strategies against this pathogen in the future.

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

Identification of Novel Antimicrobial Resistance Genes Using Machine Learning, Homology Modeling, and Molecular Docking

Antimicrobial resistance (AMR) threatens the healthcare system worldwide with the rise of emerging drug resistant infectious agents. AMR may render the current therapeutics ineffective or diminish their efficacy, and its rapid dissemination can have unmitigated health and socioeconomic consequences. Just like with many other health problems, recent computational advances including developments in machine learning or artificial intelligence hold a prodigious promise in deciphering genetic factors underlying emergence and dissemination of AMR and in aiding development of therapeutics for more efficient AMR solutions. Current machine learning frameworks focus mainly on known AMR genes and are, therefore, prone to missing genes that have not been implicated in resistance yet, including many uncharacterized genes whose functions have not yet been elucidated. Furthermore, new resistance traits may evolve from these genes leading to the rise of superbugs, and therefore, these genes need to be characterized. To infer novel resistance genes, we used complete gene sets of several bacterial strains known to be susceptible or resistant to specific drugs and associated phenotypic information within a machine learning framework that enabled prioritizing genes potentially involved in resistance. Further, homology modeling of proteins encoded by prioritized genes and subsequent molecular docking studies indicated stable interactions between these proteins and the antimicrobials that the strains containing these proteins are known to be resistant to. Our study highlights the capability of a machine learning framework to uncover novel genes that have not yet been implicated in resistance to any antimicrobials and thus could spur further studies targeted at neutralizing AMR.

Insights on the bacterial composition of Parmigiano Reggiano Natural Whey Starter by a culture-dependent and 16S rRNA metabarcoding portrait

Natural whey starters (NWS) are undefined bacterial communities produced daily from whey of the previous cheese-making round, by application of high temperature. As a result, in any dairy plant, NWS are continuously evolving, undefined mixtures of several strains and/or species of lactic acid bacteria, whose composition and performance strongly depend on the selective pressure acting during incubation. While NWS is critical to assure consistency to cheese-making process, little is known about the composition, functional features, and plant-to-plant fluctuations. Here, we integrated 16S rRNA metabarcoding and culture-dependent methods to profile bacterial communities of 10 NWS sampled in the production area of Parmigiano Reggiano cheese. 16S rRNA metabarcoding analysis revealed two main NWS community types, namely NWS type-H and NWS type-D. Lactobacillus helveticus was more abundant in NWS type-H, whilst Lactobacillus delbrueckii/St. thermophilus in NWS type-D, respectively. Based on the prediction of metagenome functions, NWS type-H samples were enriched in functional pathways related to galactose catabolism and purine metabolism, while NWS type-D in pathways related to aromatic and branched chain amino acid biosynthesis, which are flavor compound precursors. Culture-dependent approaches revealed low cultivability of individual colonies as axenic cultures and high genetic diversity in the pool of cultivable survivors. Co-culturing experiments showed that fermentative performance decreases by reducing the bacterial complexity of inoculum, suggesting that biotic interactions and cross-feeding relationships could take place in NWS communities, assuring phenotypic robustness. Even though our data cannot directly predict these ecological interactions, this study provides the basis for experiments targeted at understanding how selective regime affects composition, bacterial interaction, and fermentative performance in NWS.

Comparative analysis of Vibrio cholerae isolates from Ghana reveals variations in genome architecture and adaptation of outbreak and environmental strains

Recurrent epidemics of cholera denote robust adaptive mechanisms of Vibrio cholerae for ecological shifting and persistence despite variable stress conditions. Tracking the evolution of pathobiological traits requires comparative genomic studies of isolates from endemic areas. Here, we investigated the genetic differentiation among V. cholerae clinical and environmental isolates by highlighting the genomic divergence associated with gene decay, genome plasticity, and the acquisition of virulence and adaptive traits. The clinical isolates showed high phylogenetic relatedness due to a higher frequency of shared orthologs and fewer gene variants in contrast to the evolutionarily divergent environmental strains. Divergence of the environmental isolates is linked to extensive genomic rearrangements in regions containing mobile genetic elements resulting in numerous breakpoints, relocations, and insertions coupled with the loss of virulence determinants acf, zot, tcp, and ctx in the genomic islands. Also, four isolates possessed the CRISPR-Cas systems with spacers specific for Vibrio phages and plasmids. Genome synteny and homology analysis of the CRISPR-Cas systems suggest horizontal acquisition. The marked differences in the distribution of other phage and plasmid defense systems such as Zorya, DdmABC, DdmDE, and type-I Restriction Modification systems among the isolates indicated a higher propensity for plasmid or phage disseminated traits in the environmental isolates. Our results reveal that V. cholerae strains undergo extensive genomic rearrangements coupled with gene acquisition, reflecting their adaptation during ecological shifts and pathogenicity.