Whole genome analysis unveils genetic diversity and potential virulence determinants in Vibrio parahaemolyticus associated with disease outbreak among cultured Litopenaeus vannamei (Pacific white shrimp) in India

Genomic and Transcriptomic Analyses Reveal Multiple Strategies for Vibrio parahaemolyticus to Tolerate Sub-Lethal Concentrations of Three Antibiotics

Vibrio parahaemolyticus can cause acute gastroenteritis, wound infections, and septicemia in humans. The overuse of antibiotics in aquaculture may lead to a high incidence of the multidrug-resistant (MDR) pathogen. Nevertheless, the genome evolution of V. parahaemolyticus in aquatic animals and the mechanism of its antibiotic tolerance remain to be further deciphered. Here, we investigated the molecular basis of the antibiotic tolerance of V. parahaemolyticus isolates (n = 3) originated from shellfish and crustaceans using comparative genomic and transcriptomic analyses. The genome sequences of the V. parahaemolyticus isolates were determined (5.0-5.3 Mb), and they contained 4709-5610 predicted protein-encoding genes, of which 823-1099 genes were of unknown functions. Comparative genomic analyses revealed a number of mobile genetic elements (MGEs, n = 69), antibiotic resistance-related genes (n = 7-9), and heavy metal tolerance-related genes (n = 2-4). The V. parahaemolyticus isolates were resistant to sub-lethal concentrations (sub-LCs) of ampicillin (AMP, 512 μg/mL), kanamycin (KAN, 64 μg/mL), and streptomycin (STR, 16 μg/mL) (p < 0.05). Comparative transcriptomic analyses revealed that there were significantly altered metabolic pathways elicited by the sub-LCs of the antibiotics (p < 0.05), suggesting the existence of multiple strategies for antibiotic tolerance in V. parahaemolyticus. The results of this study enriched the V. parahaemolyticus genome database and should be useful for controlling the MDR pathogen worldwide.

Molecular docking and dynamics simulation studies uncover the host-pathogen protein-protein interactions in Penaeus vannamei and Vibrio parahaemolyticus

Shrimp aquaculture contributes significantly to global economic growth, and the whiteleg shrimp, Penaeus vannamei, is a leading species in this industry. However, Vibrio parahaemolyticus infection poses a major challenge in ensuring the success of P. vannamei aquaculture. Despite its significance in this industry, the biological knowledge of its pathogenesis remains unclear. Hence, this study was conducted to identify the interaction sites and binding affinity between several immune-related proteins of P. vannamei with V. parahaemolyticus proteins associated with virulence factors. Potential interaction sites and the binding affinity between host and pathogen proteins were identified using molecular docking and dynamics (MD) simulation. The P. vannamei-V. parahaemolyticus protein-protein interaction of Complex 1 (Ferritin-HrpE/YscL family type III secretion apparatus protein), Complex 2 (Protein kinase domain-containing protein-Chemotaxis CheY protein), and Complex 3 (GPCR-Chemotaxis CheY protein) was found to interact with -4319.76, -5271.39, and -4725.57 of the docked score and the formation of intermolecular bonds at several interacting residues. The docked scores of Complex 1, Complex 2, and Complex 3 were validated using MD simulation analysis, which revealed these complexes greatly contribute to the interactions between P. vannamei and V. parahaemolyticus proteins, with binding free energies of -22.50 kJ/mol, -30.20 kJ/mol, and -26.27 kJ/mol, respectively. This finding illustrates the capability of computational approaches to search for molecular binding sites between host and pathogen, which could increase the knowledge of Vibrio spp. infection on shrimps, which then can be used to assist in the development of effective treatment.

Genomic diversity of Vibrio spp. and metagenomic analysis of pathogens in Florida Gulf coastal waters following Hurricane Ian

IMPORTANCE

Evidence suggests warming temperatures are associated with the spread of potentially pathogenic Vibrio spp. and the emergence of human disease globally. Following Hurricane Ian, the State of Florida reported a sharp increase in the number of reported Vibrio spp. infections and deaths. Hence, monitoring of pathogens, including vibrios, and environmental parameters influencing their occurrence is critical to public health. Here, DNA sequencing was used to investigate the genomic diversity of Vibrio parahaemolyticus and Vibrio vulnificus, both potential human pathogens, in Florida coastal waters post Hurricane Ian, in October 2022. Additionally, the microbial community of water samples was profiled to detect the presence of Vibrio spp. and other microorganisms (bacteria, fungi, protists, and viruses) present in the samples. Long-term environmental data analysis showed changes in environmental parameters during and after Ian were optimal for the growth of Vibrio spp. and related pathogens. Collectively, results will be used to develop predictive risk models during climate change.

Vibrio parahaemolyticus becomes lethal to post-larvae shrimp via acquiring novel virulence factors

IMPORTANCE

As a severe emerging shrimp disease, TPD has heavily impacted the shrimp aquaculture industry and resulted in serious economic losses in China since spring 2020. This study aimed to identify the key virulent factors and related genes of the Vp TPD, for a better understanding of its pathogenicity of the novel highly lethal infectious pathogen, as well as its molecular epidemiological characteristics in China. The present study revealed that a novel protein, Vibrio high virulent protein-2 (MW >100 kDa), is responsible to the lethal virulence of V. parahaemolyticus to shrimp post-larvae. The results are essential for effectively diagnosing and monitoring novel pathogenic bacteria, like Vp TPD, in aquaculture shrimps and would be beneficial to the fisheries department in early warning of Vp TPD emergence and developing prevention strategies to reduce economic losses due to severe outbreaks of TPD. Elucidation of the key virulence genes and genomics of Vp TPD could also provide valuable information on the evolution and ecology of this emerging pathogen in aquaculture environments.

Comparative Genomic Analysis of Shrimp-Pathogenic Vibrio parahaemolyticus LC and Intraspecific Strains with Emphasis on Virulent Factors of Mobile Genetic Elements

Vibrio parahaemolyticus exhibits severe pathogenicity in humans and animals worldwide. In this study, genome sequencing and comparative analyses were conducted for in-depth characterization of the virulence factor (VF) repertoire of V. parahaemolyticus strain LC, which presented significant virulence to shrimp Litopenaeus vannamei. Strain LC, harboring two circular chromosomes and three linear plasmids, demonstrated ≥98.14% average nucleotide identities with 31 publicly available V. parahaemolyticus genomes, including 13, 11, and 7 shrimp-, human-, and non-pathogenic strains, respectively. Phylogeny analysis based on dispensable genes of pan-genome clustered 11 out of 14 shrimp-pathogenic strains and 7 out of 11 clinical strains into two distinct clades, indicating the close association between host-specific pathogenicity and accessory genes. The VFDB database revealed that 150 VFs of LC were mainly associated with the secretion system, adherence, antiphagocytosis, chemotaxis, motility, and iron uptake, whereas no homologs of the typical pathogenic genes pirA, pirB, tdh, and trh were detected. Four genes, mshB, wbfT, wbfU, and wbtI, were identified in both types of pathogenic strains but were absent in non-pathogens. Notably, a unique cluster similar to Yen-Tc, which encodes an insecticidal toxin complex, and diverse toxin-antitoxin (TA) systems, were identified on the mobile genetic elements (MGEs) of LC. Conclusively, in addition to the common VFs, various unique MGE-borne VFs, including the Yen-Tc cluster, TA components, and multiple chromosome-encoded chitinase genes, may contribute to the full spectrum of LC virulence. Moreover, V. parahaemolyticus demonstrates host-specific virulence, which potentially drives the origin and spread of pathogenic factors.

Genomic characterization and identification of virulence-related genes in Vibrio nigripulchritudo isolated from white leg shrimp Penaeus vannamei

Vibrio nigripulchritudo causes vibriosis in penaeid shrimps. Here, we used Illumina and Nanopore sequencing technologies to sequence the genomes of three of its strains (TUMSAT-V. nig1, TUMSAT-V. nig2, and TUMSAT-V. nig3) to explore opportunities for disease management. Putative virulence factors and mobile genetic elements were detected while evaluating the phylogenetic relationship of each isolated strain. The genomes consisted of two circular chromosomes (I and II) plus one or two plasmids. The size of chromosome I ranged from 4.02 to 4.07 Mb with an average GC content of 46%, while the number of predicted CDSs ranged from 3563 to 3644. The size of chromosome II ranged from 2.16 to 2.18 Mb, with an average GC content of 45.5%, and the number of predicted CDSs ranged from 1970 to 1987. Numerous virulence genes were identified related to adherence, antiphagocytosis, chemotaxis, motility, iron uptake, quorum sensing, secretion systems, and toxins in all three genomes. Higher numbers of prophages and genomic islands found in TUMSAT-V. nig1 suggest that the strain has experienced numerous horizontal gene transfer events. The presence of antimicrobial resistance genes suggests that the strains have multidrug resistance. Comparative genomic analysis showed that all three strains belonged to the same clade.

Taurine metabolism is modulated in Vibrio-infected Penaeus vannamei to shape shrimp antibacterial response and survival

BACKGROUND

Numerous microorganisms are found in aquaculture ponds, including several pathogenic bacteria. Infection of cultured animals by these pathogens results in diseases and metabolic dysregulation. However, changes in the metabolic profiles that occur at different infection stages in the same ponds and how these metabolic changes can be modulated by exogenous metabolites in Penaeus vannamei remain unknown.

RESULTS

Here, we collected gastrointestinal tract (GIT) samples from healthy, diseased, and moribund P. vannamei in the same aquaculture pond for histological, metabolic, and transcriptome profiling. We found that diseased and moribund shrimp with empty GITs and atrophied hepatopancreas were mainly infected with Vibrio parahaemolyticus and Vibrio harveyi. Although significant dysregulation of crucial metabolites and their enzymes were observed in diseased and moribund shrimps, diseased shrimp expressed high levels of taurine and taurine metabolism-related enzymes, while moribund shrimp expressed high levels of hypoxanthine and related metabolism enzymes. Moreover, a strong negative correlation was observed between taurine levels and the relative abundance of V. parahaemolyticus and V. harveyi. Besides, exogenous taurine enhanced shrimp survival against V. parahaemolyticus challenge by increasing the expression of key taurine metabolism enzymes, mainly, cysteine dioxygenase (CDO) and cysteine sulfinic acid decarboxylase (CSD).

CONCLUSIONS

Our study revealed that taurine metabolism could be modulated by exogenous supplementation to improve crustacean immune response against pathogenic microbes. Video Abstract.

Genomic and Transcriptomic Analysis Reveal Multiple Strategies for the Cadmium Tolerance in Vibrio parahaemolyticus N10-18 Isolated from Aquatic Animal Ostrea gigas Thunberg

The waterborne Vibrio parahaemolyticus can cause acute gastroenteritis, wound infection, and septicemia in humans. Pollution of heavy metals in aquatic environments is proposed to link high incidence of the multidrug-resistant (MDR) pathogen. Nevertheless, the genome evolution and heavy metal tolerance mechanism of V. parahaemolyticus in aquatic animals remain to be largely unveiled. Here, we overcome the limitation by characterizing an MDR V. parahaemolyticus N10-18 isolate with high cadmium (Cd) tolerance using genomic and transcriptomic techniques. The draft genome sequence (4,910,080 bp) of V. parahaemolyticus N10-18 recovered from Ostrea gigas Thunberg was determined, and 722 of 4653 predicted genes had unknown function. Comparative genomic analysis revealed mobile genetic elements (n = 11) and heavy metal and antibiotic-resistance genes (n = 38 and 7). The bacterium significantly changed cell membrane structure to resist the Cd2+ (50 μg/mL) stress (p < 0.05). Comparative transcriptomic analysis revealed seven significantly altered metabolic pathways elicited by the stress. The zinc/Cd/mercury/lead transportation and efflux and the zinc ATP-binding cassette (ABC) transportation were greatly enhanced; metal and iron ABC transportation and thiamine metabolism were also up-regulated; conversely, propanoate metabolism and ribose and maltose ABC transportation were inhibited (p < 0.05). The results of this study demonstrate multiple strategies for the Cd tolerance in V. parahaemolyticus.

Comparative genome analyses of five Vibrio penaeicida strains provide insights into their virulence-related factors

Vibrio penaeicida (family Vibrionaceae) is an important bacterial pathogen that affects Japanese shrimp aquaculture. Only two whole-genome sequences of V. penaeicida are publicly available, which has hampered our understanding of the pathogenesis of shrimp vibriosis caused by this bacterium. To gain insight into the genetic features, evolution and pathogenicity of V. penaeicida, we sequenced five V. penaeicida strains (IFO 15640^T, IFO 15641, IFO 15642, TUMSAT-OK1 and TUMSAT-OK2) and performed comparative genomic analyses. Virulence factors and mobile genetic elements were detected. Furthermore, average nucleotide identities (ANIs), clusters of orthologous groups and phylogenetic relationships were evaluated. The V. penaeicida genome consists of two circular chromosomes. Chromosome I sizes ranged from 4.1 to 4.3 Mb, the GC content ranged from 43.9 to 44.1 %, and the number of predicted protein-coding sequences (CDSs) ranged from 3620 to 3782. Chromosome II sizes ranged from 2.2 to 2.4 Mb, the GC content ranged from 43.5 to 43.8 %, and the number of predicted CDSs ranged from 1992 to 2273. All strains except IFO 15641 harboured one plasmid, having sizes that ranged from 150 to 285 kb. All five genomes had typical virulence factors, including adherence, anti-phagocytosis, flagella-related proteins and toxins (repeats-in-toxin and thermolabile haemolysin). The genomes also contained factors responsible for iron uptake and the type II, IV and VI secretion systems. The genome of strain TUMSAT-OK2 tended to encode more prophage regions than the other strains, whereas the genome of strain IFO 15640^T had the highest number of regions encoding genomic islands. For comparative genome analysis, we used V. penaeicida (strain CAIM 285^T) as a reference strain. ANIs between strain CAIM 285^T and the five V. penaeicida strains were >95 %, which indicated that these strains belong to the same species. Orthology cluster analysis showed that strains TUMSAT-OK1 and TUMSAT-OK2 had the greatest number of shared gene clusters, followed by strains CAIM 285^T and IFO 15640^T. These strains were also the most closely related to each other in a phylogenetic analysis. This study presents the first comparative genome analysis of V. penaeicida and these results will be useful for understanding the pathogenesis of this bacterium.

Characterization of GefA, a GGEEF domain-containing protein that modulates Vibrio parahaemolyticus motility, biofilm formation, and virulence

Vibrio parahaemolyticus is a significant food-borne pathogen that causes economic and public health problems worldwide and has the high capacity to adapt to diverse environments and hosts. The second messenger cyclic diguanylate monophosphate (c-di-GMP) allows bacteria to shift from a planktonic form to a communal multicellular lifestyle and plays an important role in bacterial survival and transmission. Here we characterized single-domain c-di-GMP synthetases in V. parahaemolyticus and identified a novel GGEEF domain-containing protein designated GefA that modulates bacterial swarming motility, biofilm formation, and virulence. GefA inhibits swarming motility by regulating the expression of lateral flagella, while it enhances biofilm formation by controlling exopolysaccharide biosynthesis. Under high-c-di-GMP conditions caused by scrABC knock-out, we found that GefA is bifunctional, as it has no effect on swarming motility but retains the ability to regulate biofilm formation. Subsequent studies suggested that GefA regulates the expression of type III secretion system 1 (T3SS1), which is an important virulence factor in V. parahaemolyticus. Here, we also revealed that the flagella participate in the infection of V. parahaemolyticus. We found that both the T3SS1 and flagella contribute to the GefA-mediated virulence of V. parahaemolyticus in the zebrafish model. Our results expand the knowledge of the V. parahaemolyticus c-di-GMP synthetases and their roles in social behaviors and pathogenicity. Importance The c-di-GMP metabolic enzymes constitute one of the largest clusters of potential orthologues in V. parahaemolyticus. However, the specific roles that these individual c-di-GMP metabolic enzymes play are largely unknown. Here, we identified a GGEEF domain-containing protein designated GefA that regulates bacterial behaviors and virulence. We also demonstrated that flagella participate in the infection of this bacterium, through which GefA regulates the bacterial virulence. To our knowledge, the roles that c-di-GMP and flagella play in V. parahaemolyticus virulence have never been revealed before. Our findings contribute to a better understanding of the function of c-di-GMP and its synthetases in V. parahaemolyticus.