Phages Tune in to Host Cell Quorum Sensing.

Phages must be perfectly attuned to bacterial host cell physiology to ensure their optimal survival. Silpe and Bassler show that a Vibrio phage uses the host quorum-sensing pathway to trigger production of viral progeny at high cell density.

A Reverse Ecology Approach Based on a Biological Definition of Microbial Populations.

Delineating ecologically meaningful populations among microbes is important for identifying their roles in environmental and host-associated microbiomes. Here, we introduce a metric of recent gene flow, which when applied to co-existing microbes, identifies congruent genetic and ecological units separated by strong gene flow discontinuities from their next of kin. We then develop a pipeline to identify genome regions within these units that show differential adaptation and allow mapping of populations onto environmental variables or host associations. Using this reverse ecology approach, we show that the human commensal bacterium Ruminococcus gnavus breaks up into sharply delineated populations that show different associations with health and disease. Defining populations by recent gene flow in this way will facilitate the analysis of bacterial and archaeal genomes using ecological and evolutionary theory developed for plants and animals, thus allowing for testing unifying principles across all biology.

A Host-Produced Quorum-Sensing Autoinducer Controls a Phage Lysis-Lysogeny Decision.

Vibrio cholerae uses a quorum-sensing (QS) system composed of the autoinducer 3,5-dimethylpyrazin-2-ol (DPO) and receptor VqmA (VqmAVc), which together repress genes for virulence and biofilm formation. vqmA genes exist in Vibrio and in one vibriophage, VP882. Phage-encoded VqmA (VqmAPhage) binds to host-produced DPO, launching the phage lysis program via an antirepressor that inactivates the phage repressor by sequestration. The antirepressor interferes with repressors from related phages. Like phage VP882, these phages encode DNA-binding proteins and partner antirepressors, suggesting that they, too, integrate host-derived information into their lysis-lysogeny decisions. VqmAPhage activates the host VqmAVc regulon, whereas VqmAVc cannot induce phage-mediated lysis, suggesting an asymmetry whereby the phage influences host QS while enacting its own lytic-lysogeny program without interference. We reprogram phages to activate lysis in response to user-defined cues. Our work shows that a phage, causing bacterial infections, and V. cholerae, causing human infections, rely on the same signal molecule for pathogenesis.

A qrr noncoding RNA deploys four different regulatory mechanisms to optimize quorum-sensing dynamics.

Quorum sensing is a cell-cell communication process that bacteria use to transition between individual and social lifestyles. In vibrios, homologous small RNAs called the Qrr sRNAs function at the center of quorum-sensing pathways. The Qrr sRNAs regulate multiple mRNA targets including those encoding the quorum-sensing regulatory components luxR, luxO, luxM, and aphA. We show that a representative Qrr, Qrr3, uses four distinct mechanisms to control its particular targets: the Qrr3 sRNA represses luxR through catalytic degradation, represses luxM through coupled degradation, represses luxO through sequestration, and activates aphA by revealing the ribosome binding site while the sRNA itself is degraded. Qrr3 forms different base-pairing interactions with each mRNA target, and the particular pairing strategy determines which regulatory mechanism occurs. Combined mathematical modeling and experiments show that the specific Qrr regulatory mechanism employed governs the potency, dynamics, and competition of target mRNA regulation, which in turn, defines the overall quorum-sensing response.

Frequent nonhomologous replacement of replicative helicase loaders by viruses in Vibrionaceae.

Several microbial genomes lack textbook-defined essential genes. If an essential gene is absent from a genome, then an evolutionarily independent gene of unknown function complements its function. Here, we identified frequent nonhomologous replacement of an essential component of DNA replication initiation, a replicative helicase loader gene, in Vibrionaceae. Our analysis of Vibrionaceae genomes revealed two genes with unknown function, named vdhL1 and vdhL2, that were substantially enriched in genomes without the known helicase-loader genes. These genes showed no sequence similarities to genes with known function but encoded proteins structurally similar with a viral helicase loader. Analyses of genomic syntenies and coevolution with helicase genes suggested that vdhL1/2 encodes a helicase loader. The in vitro assay showed that Vibrio harveyi VdhL1 and Vibrio ezurae VdhL2 promote the helicase activity of DnaB. Furthermore, molecular phylogenetics suggested that vdhL1/2 were derived from phages and replaced an intrinsic helicase loader gene of Vibrionaceae over 20 times. This high replacement frequency implies the host's advantage in acquiring a viral helicase loader gene.

Shrimp shapes a resistance trait against vibriosis by memorizing the colonization resistance of intestinal microbiota.

Vibriosis is one of the most serious diseases that commonly occurs in aquatic animals, thus, shaping a steady inherited resistance trait in organisms has received the highest priority in aquaculture. Whereas, the mechanisms underlying the development of such a resistance trait are mostly elusive. In this study, we constructed vibriosis-resistant and susceptible families of the Pacific white shrimp Litopenaeus vannamei after four generations of artificial selection. Microbiome sequencing indicated that shrimp can successfully develop a colonization resistance trait against Vibrio infections. This trait was characterized by a microbial community structure with specific enrichment of a single probiotic species (namely Shewanella algae), and notably, its formation was inheritable and might be memorized by host epigenetic remodeling. Regardless of the infection status, a group of genes was specifically activated in the resistant family through disruption of complete methylation. Specifically, hypo-methylation and hyper-expression of genes related to lactate dehydrogenase (LDH) and iron homeostasis might provide rich sources of specific carbon (lactate) and ions for the colonization of S. algae, which directly results in the reduction of Vibrio load in shrimp. Lactate feeding increased the survival of shrimp, while knockdown of LDH gene decreased the survival when shrimp was infected by Vibrio pathogens. In addition, treatment of shrimp with the methyltransferase inhibitor 5-azacytidine resulted in upregulations of LDH and some protein processing genes, significant enrichment of S. algae, and simultaneous reduction of Vibrio in shrimp. Our results suggest that the colonization resistance can be memorized as epigenetic information by the host, which has played a pivotal role in vibriosis resistance. The findings of this study will aid in disease control and the selection of superior lines of shrimp with high disease resistance.

An inhibitory mechanism of AasS, an exogenous fatty acid scavenger: Implications for re-sensitization of FAS II antimicrobials.

Antimicrobial resistance is an ongoing "one health" challenge of global concern. The acyl-ACP synthetase (termed AasS) of the zoonotic pathogen Vibrio harveyi recycles exogenous fatty acid (eFA), bypassing the requirement of type II fatty acid synthesis (FAS II), a druggable pathway. A growing body of bacterial AasS-type isoenzymes compromises the clinical efficacy of FAS II-directed antimicrobials, like cerulenin. Very recently, an acyl adenylate mimic, C10-AMS, was proposed as a lead compound against AasS activity. However, the underlying mechanism remains poorly understood. Here we present two high-resolution cryo-EM structures of AasS liganded with C10-AMS inhibitor (2.33 Å) and C10-AMP intermediate (2.19 Å) in addition to its apo form (2.53 Å). Apart from our measurements for C10-AMS' Ki value of around 0.6 µM, structural and functional analyses explained how this inhibitor interacts with AasS enzyme. Unlike an open state of AasS, ready for C10-AMP formation, a closed conformation is trapped by the C10-AMS inhibitor. Tight binding of C10-AMS blocks fatty acyl substrate entry, and therefore inhibits AasS action. Additionally, this intermediate analog C10-AMS appears to be a mixed-type AasS inhibitor. In summary, our results provide the proof of principle that inhibiting salvage of eFA by AasS reverses the FAS II bypass. This facilitates the development of next-generation anti-bacterial therapeutics, esp. the dual therapy consisting of C10-AMS scaffold derivatives combined with certain FAS II inhibitors.

YTHDC1/CRM1 Facilitates m6A-Modified circRNA388 Nuclear Export to Induce Coelomocyte Autophagy via the miR-2008/ULK Axis in Apostichopus japonicus.

N 6-methyladenosine (m6A), the most prevalent internal modification in eukaryotic RNA, was able to mediate circular RNA (circRNA) function in many immune processes. Nevertheless, the functional role of m6A-modified circRNAs in innate immunity of invertebrates remained unclear. In this study, we identified m6A-modified circRNA388 from cultured sea cucumber (Apostichopus japonicus) coelomocytes, which was mainly detected in cytoplasm after Vibrio splendidus infection. A knockdown assay indicated that cytoplasm circRNA388 promoted coelomocyte autophagy and decreased the number of intracellular V. splendidus. Mechanistically, the circRNA388 in the cytoplasm directly sponged miR-2008 to block its interaction with Unc-51-like kinase 1 from A. japonicus (AjULK) and further promoted autophagy to resist V. splendidus infection. More importantly, we found that m6A modification was vital to circRNA388 nuclear export with YTH domain-containing protein 1 from A. japonicus (AjYTHDC1) as the reader. AjYTHDC1 facilitated the nuclear export of m6A-modified circRNA388 via interaction with exportin-1 (chromosomal maintenance 1) from A. japonicus (AjCRM1). Knockdown of AjCRM1 could significantly decrease the content of cytoplasm circRNA388. Overall, our results provide the first evidence that nuclear export of m6A-modified circRNA388 is dependent on the novel AjCRM1 to our knowledge, which was further promoted coelomocyte autophagy by miR-2008/AjULK axis to clear intracellular V. splendidus.

Chemical Proteomics Reveals Nε-Fatty-Acylation of Septins by Rho Inactivation Domain (RID) of the Vibrio MARTX Toxin to Alter Septin Localization and Organization.

Vibrio species, the Gram-negative bacterial pathogens causing cholera and sepsis, produce multiple secreted virulence factors for infection and pathogenesis. Among these is the multifunctional-autoprocessing repeats-in-toxin (MARTX) toxin that releases several critical effector domains with distinct functions inside eukaryotic host cells. One such effector domain, the Rho inactivation domain (RID), has been discovered to catalyze long-chain Nε-fatty-acylation on lysine residues of Rho GTPases, causing inactivation of Rho GTPases and disruption of the host actin cytoskeleton. However, whether RID modifies other host proteins to exert additional functions remains to be determined. Herein, we describe the integration of bioorthogonal chemical labeling and quantitative proteomics to globally profile the target proteins modified by RID in living cells. More than 246 proteins are identified as new RID substrates, including many involved in GTPase regulation, cytoskeletal organization, and cell division. We demonstrate that RID extensively Nε-fatty-acylates septin proteins, the fourth cytoskeletal component of mammalian cells with important roles in diverse cellular processes. While affinity purification and mass spectrometry analysis show that RID-mediated Nε-fatty-acylation does not affect septin-septin interactions, this modification increases the membrane association of septins and confers localization to detergent-resistant membrane rafts. As a result, the filamentous assembly and organization of septins are disrupted by RID-mediated Nε-fatty-acylation, further contributing to cytoskeletal and mitotic defects that phenocopy the effects of septin depletion. Overall, our work greatly expands the substrate scope and function of RID and demonstrates the role of RID-mediated Nε-fatty-acylation in manipulating septin localization and organization.

Identification of promoter activity in gene-less cassettes from Vibrionaceae superintegrons.

Integrons are genetic platforms that acquire new genes encoded in integron cassettes (ICs), building arrays of adaptive functions. ICs generally encode promoterless genes, whose expression relies on the platform-associated Pc promoter, with the cassette array functioning as an operon-like structure regulated by the distance to the Pc. This is relevant in large sedentary chromosomal integrons (SCIs) carrying hundreds of ICs, like those in Vibrio species. We selected 29 gene-less cassettes in four Vibrio SCIs, and explored whether their function could be related to the transcription regulation of adjacent ICs. We show that most gene-less cassettes have promoter activity on the sense strand, enhancing the expression of downstream cassettes. Additionally, we identified the transcription start sites of gene-less ICs through 5'-RACE. Accordingly, we found that most of the superintegron in Vibrio cholerae is not silent. These promoter cassettes can trigger the expression of a silent dfrB9 cassette downstream, increasing trimethoprim resistance >512-fold in V. cholerae and Escherichia coli. Furthermore, one cassette with an antisense promoter can reduce trimethoprim resistance when cloned downstream. Our findings highlight the regulatory role of gene-less cassettes in the expression of adjacent cassettes, emphasizing their significance in SCIs and their clinical importance if captured by mobile integrons.

Scavenger receptor B2, a type III membrane pattern recognition receptor, senses LPS and activates the IMD pathway in crustaceans.

The immune deficiency (IMD) pathway is critical for elevating host immunity in both insects and crustaceans. The IMD pathway activation in insects is mediated by peptidoglycan recognition proteins, which do not exist in crustaceans, suggesting a previously unidentified mechanism involved in crustacean IMD pathway activation. In this study, we identified a Marsupenaeus japonicus B class type III scavenger receptor, SRB2, as a receptor for activation of the IMD pathway. SRB2 is up-regulated upon bacterial challenge, while its depletion exacerbates bacterial proliferation and shrimp mortality via abolishing the expression of antimicrobial peptides. The extracellular domain of SRB2 recognizes bacterial lipopolysaccharide (LPS), while its C-terminal intracellular region containing a cryptic RHIM-like motif interacts with IMD, and activates the pathway by promoting nuclear translocation of RELISH. Overexpressing shrimp SRB2 in Drosophila melanogaster S2 cells potentiates LPS-induced IMD pathway activation and diptericin expression. These results unveil a previously unrecognized SRB2-IMD axis responsible for antimicrobial peptide induction and restriction of bacterial infection in crustaceans and provide evidence of biological diversity of IMD signaling in animals. A better understanding of the innate immunity of crustaceans will permit the optimization of prevention and treatment strategies against the arising shrimp diseases.

NTR 2.0: a rationally engineered prodrug-converting enzyme with substantially enhanced efficacy for targeted cell ablation.

Transgenic expression of bacterial nitroreductase (NTR) enzymes sensitizes eukaryotic cells to prodrugs such as metronidazole (MTZ), enabling selective cell-ablation paradigms that have expanded studies of cell function and regeneration in vertebrates. However, first-generation NTRs required confoundingly toxic prodrug treatments to achieve effective cell ablation, and some cell types have proven resistant. Here we used rational engineering and cross-species screening to develop an NTR variant, NTR 2.0, which exhibits ~100-fold improvement in MTZ-mediated cell-specific ablation efficacy, eliminating the need for near-toxic prodrug treatment regimens. NTR 2.0 therefore enables sustained cell-loss paradigms and ablation of previously resistant cell types. These properties permit enhanced interrogations of cell function, extended challenges to the regenerative capacities of discrete stem cell niches, and novel modeling of chronic degenerative diseases. Accordingly, we have created a series of bipartite transgenic reporter/effector resources to facilitate dissemination of NTR 2.0 to the research community.

Two novel genes identified by large-scale transcriptomic analysis are essential for biofilm and rugose colony development of Vibrio vulnificus.

Many pathogenic bacteria form biofilms to survive under environmental stresses and host immune defenses. Differential expression (DE) analysis of the genes in biofilm and planktonic cells under a single condition, however, has limitations to identify the genes essential for biofilm formation. Independent component analysis (ICA), a machine learning algorithm, was adopted to comprehensively identify the biofilm genes of Vibrio vulnificus, a fulminating human pathogen, in this study. ICA analyzed the large-scale transcriptome data of V. vulnificus cells under various biofilm and planktonic conditions and then identified a total of 72 sets of independently co-regulated genes, iModulons. Among the three iModulons specifically activated in biofilm cells, BrpT-iModulon mainly consisted of known genes of the regulon of BrpT, a transcriptional regulator controlling biofilm formation of V. vulnificus. Interestingly, the BrpT-iModulon additionally contained two novel genes, VV1_3061 and VV2_1694, designated as cabH and brpN, respectively. cabH and brpN were shared in other Vibrio species and not yet identified by DE analyses. Genetic and biochemical analyses revealed that cabH and brpN are directly up-regulated by BrpT. The deletion of cabH and brpN impaired the robust biofilm and rugose colony formation. CabH, structurally similar to the previously known calcium-binding matrix protein CabA, was essential for attachment to the surface. BrpN, carrying an acyltransferase-3 domain as observed in BrpL, played an important role in exopolysaccharide production. Altogether, ICA identified two novel genes, cabH and brpN, which are regulated by BrpT and essential for the development of robust biofilms and rugose colonies of V. vulnificus.

Novel secreted STPKLRR from Vibrio splendidus AJ01 promotes pathogen internalization via mediating tropomodulin phosphorylation dependent cytoskeleton rearrangement.

We previously demonstrated that the flagellin of intracellular Vibrio splendidus AJ01 could be specifically identified by tropomodulin (Tmod) and further mediate p53-dependent coelomocyte apoptosis in the sea cucumber Apostichopus japonicus. In higher animals, Tmod serves as a regulator in stabilizing the actin cytoskeleton. However, the mechanism on how AJ01 breaks the AjTmod-stabilized cytoskeleton for internalization remains unclear. Here, we identified a novel AJ01 Type III secretion system (T3SS) effector of leucine-rich repeat-containing serine/threonine-protein kinase (STPKLRR) with five LRR domains and a serine/threonine kinase (STYKc) domain, which could specifically interact with tropomodulin domain of AjTmod. Furthermore, we found that STPKLRR directly phosphorylated AjTmod at serine 52 (S52) to reduce the binding stability between AjTmod and actin. After AjTmod dissociated from actin, the F-actin/G-actin ratio decreased to induce cytoskeletal rearrangement, which in turn promoted the internalization of AJ01. The STPKLRR knocked out strain could not phosphorylated AjTmod and displayed lower internalization capacity and pathogenic effect compared to AJ01. Overall, we demonstrated for the first time that the T3SS effector STPKLRR with kinase activity was a novel virulence factor in Vibrio and mediated self-internalization by targeting host AjTmod phosphorylation dependent cytoskeleton rearrangement, which provided a candidate target to control AJ01 infection in practice.

Eukaryotic catecholamine hormones influence the chemotactic control of Vibrio campbellii by binding to the coupling protein CheW.

SignificanceHost-emitted stress hormones significantly influence the growth and behavior of various bacterial species; however, their cellular targets have so far remained elusive. Here, we used customized probes and quantitative proteomics to identify the target of epinephrine and the α-adrenoceptor agonist phenylephrine in live cells of the aquatic pathogen Vibrio campbellii. Consequently, we have discovered the coupling protein CheW, which is in the center of the chemotaxis signaling network, as a target of both molecules. We not only demonstrate direct ligand binding to CheW but also elucidate how this affects chemotactic control. These findings are pivotal for further research on hormone-specific effects on bacterial behavior.

Gas-Selective Catalytic Regulation by a Newly Identified Globin-Coupled Sensor Phosphodiesterase Containing an HD-GYP Domain from the Human Pathogen Vibrio fluvialis.

Globin-coupled sensors constitute an important family of heme-based gas sensors, an emerging class of heme proteins. In this study, we have identified and characterized a globin-coupled sensor phosphodiesterase containing an HD-GYP domain (GCS-HD-GYP) from the human pathogen Vibrio fluvialis, which is an emerging foodborne pathogen of increasing public health concern. The amino acid sequence encoded by the AL536_01530 gene from V. fluvialis indicated the presence of an N-terminal globin domain and a C-terminal HD-GYP domain, with HD-GYP domains shown previously to display phosphodiesterase activity toward bis(3',5')-cyclic dimeric guanosine monophosphate (c-di-GMP), a bacterial second messenger that regulates numerous important physiological functions in bacteria, including in bacterial pathogens. Optical absorption spectral properties of GCS-HD-GYP were found to be similar to those of myoglobin and hemoglobin and of other bacterial globin-coupled sensors. The binding of O2 to the Fe(II) heme iron complex of GCS-HD-GYP promoted the catalysis of the hydrolysis of c-di-GMP to its linearized product, 5'-phosphoguanylyl-(3',5')-guanosine (pGpG), whereas CO and NO binding did not enhance the catalysis, indicating a strict discrimination of these gaseous ligands. These results shed new light on the molecular mechanism of gas-selective catalytic regulation by globin-coupled sensors, with these advances apt to lead to a better understanding of the family of globin-coupled sensors, a still growing family of heme-based gas sensors. In addition, given the importance of c-di-GMP in infection and virulence, our results suggested that GCS-HD-GYP could play an important role in the ability of V. fluvialis to sense O2 and NO in the context of host-pathogen interactions.

Using Vibrio natriegens for High-Yield Production of Challenging Expression Targets and for Protein Perdeuteration.

Production of soluble proteins is essential for structure/function studies; however, this usually requires milligram amounts of protein, which can be difficult to obtain with traditional expression systems. Recently, the Gram-negative bacterium Vibrio natriegens emerged as a novel and alternative host platform for production of proteins in high yields. Here, we used a commercial strain derived from V. natriegens (Vmax X2) to produce soluble bacterial and fungal proteins in milligram scale, which we struggled to achieve in Escherichia coli. These proteins include the cholera toxin (CT) and N-acetyl glucosamine-binding protein A (GbpA) from Vibrio cholerae, the heat-labile enterotoxin (LT) from E. coli and the fungal nematotoxin CCTX2 from Coprinopsis cinerea. CT, GbpA, and LT are secreted by the Type II secretion system in their natural hosts. When these three proteins were produced in Vmax, they were also secreted and could be recovered from the growth media. This simplified the downstream purification procedure and resulted in considerably higher protein yields compared to production in E. coli (6- to 26-fold increase). We also tested Vmax for protein perdeuteration using deuterated minimal media with deuterium oxide as solvent and achieved a 3-fold increase in yield compared to the equivalent protocol in E. coli. This is good news, since isotopic labeling is expensive and often ineffective but represents a necessary prerequisite for some structural biology techniques. Thus, Vmax represents a promising host for production of challenging expression targets and for protein perdeuteration in amounts suitable for structural biology studies.

High sensitivity and low-cost flavin luciferase (FLUXVc)-based reporter gene for mammalian cell expression.

Luciferase-based gene reporters generating bioluminescence signals are important tools for biomedical research. Amongst the luciferases, flavin-dependent enzymes use the most economical chemicals. However, their applications in mammalian cells are limited due to their low signals compared to other systems. Here, we constructed Flavin Luciferase from Vibrio campbellii (Vc) for Mammalian Cell Expression (FLUXVc) by engineering luciferase from V. campbellii (the most thermostable bacterial luciferase reported to date) and optimizing its expression and reporter assays in mammalian cells which can improve the bioluminescence light output by >400-fold as compared to the nonengineered version. We found that the FLUXVc reporter gene can be overexpressed in various cell lines and showed outstanding signal-to-background in HepG2 cells, significantly higher than that of firefly luciferase (Fluc). The combined use of FLUXVc/Fluc as target/control vectors gave the most stable signals, better than the standard set of Fluc(target)/Rluc(control). We also demonstrated that FLUXVc can be used for testing inhibitors of the NF-κB signaling pathway. Collectively, our results provide an optimized method for using the more economical flavin-dependent luciferase in mammalian cells.

Recovering high-quality bacterial genomes from cross-contaminated cultures: a case study of marine Vibrio campbellii.

BACKGROUND: Environmental monitoring of bacterial pathogens is critical for disease control in coastal marine ecosystems to maintain animal welfare and ecosystem function and to prevent significant economic losses. This requires accurate taxonomic identification of environmental bacterial pathogens, which often cannot be achieved by commonly used genetic markers (e.g., 16S rRNA gene), and an understanding of their pathogenic potential based on the information encoded in their genomes. The decreasing costs of whole genome sequencing (WGS), combined with newly developed bioinformatics tools, now make it possible to unravel the full potential of environmental pathogens, beyond traditional microbiological approaches. However, obtaining a high-quality bacterial genome, requires initial cultivation in an axenic culture, which is a bottleneck in environmental microbiology due to cross-contamination in the laboratory or isolation of non-axenic strains. RESULTS: We applied WGS to determine the pathogenic potential of two Vibrio isolates from coastal seawater. During the analysis, we identified cross-contamination of one of the isolates and decided to use this dataset to evaluate the possibility of bioinformatic contaminant removal and recovery of bacterial genomes from a contaminated culture. Despite the contamination, using an appropriate bioinformatics workflow, we were able to obtain high quality and highly identical genomes (Average Nucleotide Identity value 99.98%) of one of the Vibrio isolates from both the axenic and the contaminated culture. Using the assembled genome, we were able to determine that this isolate belongs to a sub-lineage of Vibrio campbellii associated with several diseases in marine organisms. We also found that the genome of the isolate contains a novel Vibrio plasmid associated with bacterial defense mechanisms and horizontal gene transfer, which may offer a competitive advantage to this putative pathogen. CONCLUSIONS: Our study shows that, using state-of-the-art bioinformatics tools and a sufficient sequencing effort, it is possible to obtain high quality genomes of the bacteria of interest and perform in-depth genomic analyses even in the case of a contaminated culture. With the new isolate and its complete genome, we are providing new insights into the genomic characteristics and functional potential of this sub-lineage of V. campbellii. The approach described here also highlights the possibility of recovering complete bacterial genomes in the case of non-axenic cultures or obligatory co-cultures.

Prophages carrying Zot toxins on different Vibrio genomes: A comprehensive assessment using multilayer networks.

Vibrios, a group of bacteria that are among the most abundant in marine environments, include several species such as Vibrio cholerae and Vibrio parahaemolyticus, which can be pathogenic to humans. Some species of Vibrio contain prophages within their genomes. These prophages can carry genes that code for toxins, such as the zonula occludens toxin (Zot), which contribute to bacterial virulence. Understanding the association between different Vibrio species, prophages and Zot genes can provide insights into their ecological interactions. In this study, we evaluated 4619 Vibrio genomes from 127 species to detect the presence of prophages carrying the Zot toxin. We found 2030 potential prophages with zot-like genes in 43 Vibrio species, showing a non-random association within a primarily modular interaction network. Some prophages, such as CTX or Vf33, were associated with specific species. In contrast, prophages phiVCY and VfO3K6 were found in 28 and 20 Vibrio species, respectively. We also identified six clusters of Zot-like sequences in prophages, with the ZOT2 cluster being the most frequent, present in 34 Vibrio species. This analysis helps to understand the distribution patterns of zot-containing prophages across Vibrio genomes and the potential routes of Zot-like toxin dissemination.