Comparative genomic analysis unravels the transmission pattern and intra-species divergence of acute hepatopancreatic necrosis disease (AHPND)-causing Vibrio parahaemolyticus strains

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.

Biological Function of Prophage-Related Gene Cluster ΔVpaChn25_RS25055~ΔVpaChn25_0714 of Vibrio parahaemolyticus CHN25

Vibrio parahaemolyticus is the primary foodborne pathogen known to cause gastrointestinal infections in humans. Nevertheless, the molecular mechanisms of V. parahaemolyticus pathogenicity are not fully understood. Prophages carry virulence and antibiotic resistance genes commonly found in Vibrio populations, and they facilitate the spread of virulence and the emergence of pathogenic Vibrio strains. In this study, we characterized three such genes, VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055, within the largest prophage gene cluster in V. parahaemolyticus CHN25. The deletion mutants ΔVpaChn25_RS25055, ΔVpaChn25_0713, ΔVpaChn25_0714, and ΔVpaChn25_RS25055-0713-0714 were derived with homologous recombination, and the complementary mutants ΔVpaChn25_0713-com, ΔVpaChn25_0714-com, ΔVpaChn25_RS25055-com, ΔVpaChn25_RS25055-0713-0714-com were also constructed. In the absence of the VpaChn25_RS25055, VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055-0713-0714 genes, the mutants showed significant reductions in low-temperature survivability and biofilm formation (p < 0.001). The ΔVpaChn25_0713, ΔVpaChn25_RS25055, and ΔVpaChn25_RS25055-0713-0714 mutants were also significantly defective in swimming motility (p < 0.001). In the Caco-2 model, the above four mutants attenuated the cytotoxic effects of V. parahaemolyticus CHN25 on human intestinal epithelial cells (p < 0.01), especially the ΔVpaChn25_RS25055 and ΔVpaChn25_RS25055-0713-0714 mutants. Transcriptomic analysis showed that 15, 14, 8, and 11 metabolic pathways were changed in the ΔVpaChn25_RS25055, ΔVpaChn25_0713, ΔVpaChn25_0714, and ΔVpaChn25_RS25055-0713-0714 mutants, respectively. We labeled the VpaChn25_RS25055 gene with superfolder green fluorescent protein (sfGFP) and found it localized at both poles of the bacteria cell. In addition, we analyzed the evolutionary origins of the above genes. In summary, the prophage genes VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055 enhance V. parahaemolyticus CHN25's survival in the environment and host. Our work improves the comprehension of the synergy between prophage-associated genes and the evolutionary process of V. parahaemolyticus.

Comparative analysis unravels genetic recombination events of Vibrio parahaemolyticus recA gene

Vibrio parahaemolyticus is a gram-negative bacterium capable of causing diseases in humans and aquatic animals. The global relationships among V. parahaemolyticus genomes have been studied using multilocus sequence typing (MLST). Recently, the MLST gene recA has shown difficulties in amplification and/or a larger PCR fragment for some V. parahaemolyticus genomes due to genetic recombination. We aimed to investigate these recombination events of recA gene by analyzing 500 publicly available whole genomes from the NCBI database. The genomes with untypable recA genes were separated using BIGSdb and CGEMLST 2.0 servers, followed by annotation with RAST and NCBI pipelines. Moreover, the variable nature of V. parahaemolyticus was investigated by wgMLST analysis. The hypothetical proteins in recombinant regions were analyzed with VCIMPred tool. In the results, 3 genomes were detected with recA gene recombination, in which 2 were associated with phages and 1 to an AHPND causing strain. All 3 recombinant regions had a G + C content of 39%-40% with 15-30 ORFs, including a newly incorporated recA gene. These acquired recA genes were closely related to 3 different genera namely Aliivibrio, Photobacterium, and Vibrio. The wgMLST analysis indicated genetic recombination events occur independently among V. parahaemolyticus on a global scale. The in silico analysis revealed 4 hypothetical proteins associated with virulence factors in recombinant regions. The present study confirms, recombination events of V. parahaemolyticus recA gene, are diverse and may have an impact on the evolutionary process. Moreover, understanding these genetic recombination events of the recA gene is necessary to determine their STs and, therefore assessing epidemiological relationships.

New Insights into the Mechanism of Action of PirAB from Vibrio Parahaemolyticus

PirAB toxins secreted by Vibrio parahaemolyticus (Vp) harbor the pVA1 virulence plasmid, which causes acute hepatopancreatic necrosis disease (AHPND), an emerging disease in Penaeid shrimp that can cause 70–100% mortality and that has resulted in great economic losses since its first appearance. The cytotoxic effect of PirAB^Vp on the epithelial cells of the shrimp hepatopancreas (Hp) has been extensively documented. New insights into the biological role of the PirB^Vp subunit show that it has lectin-like activity and recognizes mucin-like O-glycosidic structures in the shrimp Hp. The search for toxin receptors can lead to a better understanding of the infection mechanisms of the pathogen and the prevention of the host disease by blocking toxin–receptor interactions using a mimetic antagonist. There is also evidence that Vp AHPND changes the community structure of the microbiota in the surrounding water, resulting in a significant reduction of several bacterial taxa, especially Neptuniibacter spp. Considering these findings, the PirAB^vp toxin could exhibit a dual role of damaging the shrimp Hp while killing the surrounding bacteria.

Whole-genome sequence analysis through online web interfaces: a review

The recent development of whole-genome sequencing technologies paved the way for understanding the genomes of microorganisms. Every whole-genome sequencing (WGS) project requires a considerable cost and a massive effort to address the questions at hand. The final step of WGS is data analysis. The analysis of whole-genome sequence is dependent on highly sophisticated bioinformatics tools that the research personal have to buy. However, many laboratories and research institutions do not have the bioinformatics capabilities to analyze the genomic data and therefore, are unable to take maximum advantage of whole-genome sequencing. In this aspect, this study provides a guide for research personals on a set of bioinformatics tools available online that can be used to analyze whole-genome sequence data of bacterial genomes. The web interfaces described here have many advantages and, in most cases exempting the need for costly analysis tools and intensive computing resources.

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

Vibrio parahaemolyticus has caused widespread mortality in Indian shrimp aquaculture in recent years. However, there are insufficient genome data for the isolates from Indian shrimp vibriosis to analyze genetic diversity and track the acquisition of genetic features that could be involved in virulence and fitness. In this study, we have performed genome analysis of V. parahaemolyticus isolated from moribund shrimps collected from shrimp farms along coastal Karnataka, India, for better understanding of their diversity and virulence. Five newly sequenced genomes of V. parahaemolyticus along with 40 genomes retrieved from NCBI were subjected to comparative genome analysis. The sequenced genomes had an overall genome size of 5.2 Mb. MLST analysis and core genome phylogenomic analysis revealed considerable genetic diversity among the isolates obtained from the moribund shrimps. Interestingly, none of the V. parahaemolyticus isolates possessed the classical features (PirAB) of the strains associated with Acute Hepatopancreatic Necrosis Disease (AHPND). This study also revealed the presence of multiple virulence attributes, including ZOT, ACE and RTX toxins, secretion systems, and mobile genetic elements. The findings of this study provide insights into the possible transition of an environmental V. parahaemolyticus to emerge as pathogens of aquaculture species by increasing its virulence and host adaptation. Future studies focusing on continuous genomic surveillance of V. parahaemolyticus are required to study the evolution and transmission of new variants in shrimp aquaculture, as well as to design and implement biosecurity programs to prevent disease outbreaks.

Characterization of vB_VpaP_MGD2, a newly isolated bacteriophage with biocontrol potential against multidrug-resistant Vibrio parahaemolyticus

Vibrio parahaemolyticus is a major foodborne pathogen and is also pathogenic to shrimp. Due to the emergence of multidrug-resistant V. parahaemolyticus strains, bacteriophages have shown promise as antimicrobial agents that could be used for controlling antibiotic-resistant strains. Here, a V. parahaemolyticus phage, vB_VpaP_MGD2, was isolated from a clam (Meretrix meretrix) and further characterized to evaluate its potential capability for biocontrol. Podophage vB_VpaP_MGD2 had a wide host range and was able to lyse 27 antibiotic-resistant V. parahaemolyticus strains. A one-step growth curve showed that vB_VpaP_MGD2 has a short latent period of 10 min and a large burst size of 244 phages per cell. Phage vB_VpaP_MGD2 was able to tolerate a wide range of temperature (30 °C-50 °C) and pH (pH 3-pH 10). Two multidrug-resistant strains (SH06 and SA411) were suppressed by treatment with phage vB_VpaP_MGD2 at a multiplicity of infection of 100 for 24 h without apparent regrowth of bacterial populations. The frequency of mutations causing bacteriophage resistance was relatively low (3.1 × 10^-6). Phage vB_VpaP_MGD2 has a double-stranded DNA with a genome size of 45,105 bp. Among the 48 open reading frames annotated in the genome, no lysogenic genes or virulence genes were detected. Sequence comparisons suggested that vB_VpaP_MGD2 is a member of a new species in the genus Zindervirus within the subfamily Autographivirinae. This is the first report of a member of the genus Zindervirus that can infect V. parahaemolyticus. These findings suggest that vB_VpaP_MGD2 may be a candidate biocontrol agent against early mortality syndrome/acute hepatopancreatic necrosis disease (EMS/AHPND) caused by multidrug-resistant V. parahaemolyticus in shrimp production.

Coordination of primary metabolism and virulence factors expression mediates the virulence of Vibrio parahaemolyticus towards cultured shrimp (Penaeus vannamei)

AIMS

Acute hepatopancreatic necrosis disease (AHPND) caused by Vibrio parahaemolyticus has emerged as a severe bacterial disease of cultured shrimp. To identify the key virulence factors, two AHPND-causing V. parahaemolyticus (Vp_AHPND ) strains (123 and 137) and two non-Vp_AHPND strains (HZ56 and ATCC 17082) were selected.

METHODS AND RESULTS

Challenge tests showed that the four strains exhibited different virulence towards shrimp with cumulative mortalities at 48 h postinfection (hpi) ranging from 10 to 92%. The expression of pirAB^VP in strain 123 and 137 was not significantly different. Genomic analysis revealed that the two Vp_AHPND strains contain a plasmid with the PirAB^VP toxins (pirAB^VP ) flanked by the insertion sequence (ISVal1) that has been identified in various locations of chromosomes in Vp_AHPND strains. The two Vp_AHPND strains possessed almost identical virulence factors, while ISVal1 disrupted three genes related to flagellar motility in strain 137. Phenotype assay showed that strain 123 possessed the highest growth rate and swimming motility, followed by strain 137, suggesting that the disruption of essential genes mediated by ISVal1 significantly affected the virulence level. Transcriptome analysis of two Vp_AHPND strains (123 and 137) further suggested that virulence genes related to the capsule, flagella and primary metabolism were highly expressed in strain 123.

CONCLUSIONS

Here for the first time, it is demonstrated that the virulence of Vp_AHPND is not only determined by the expression of pirAB^VP , but also is mediated by ISVal1 which affects the genes involved in flagellar motility and primary metabolism.

SIGNIFICANCE AND IMPACT OF THE STUDY

The genomic and transcriptomic analysis of Vp_AHPND strains provides valuable information on the virulence factors affecting the pathogenicity of Vp_AHPND.

Delineating the key virulence factors and intraspecies divergence of Vibrio harveyi via whole-genome sequencing

Vibrio harveyi is one of the major pathogens in aquaculture. To identify the key virulence factors affecting pathogenesis of V. harveyi towards fish, we conducted a field investigation for three representative fish farms infected with V. harveyi. Multilocus sequence typing (MLST) and whole-genome sequencing were conducted to delineate the phylogenetic relationship and genetic divergence of V. harveyi. A total of 25 V. harveyi strains were isolated from the diseased fish and groundwater and were subtyped into 12 sequence types by MLST. Five virulence genes, mshB, pilA, hutR, ureB, and ureG, were variably presented in the sequenced strains. The virulence gene profiles strongly correlated with the distinct pathogenicity of V. harveyi strains, with a strain harboring all five genes exhibiting the highest virulence towards fish. Phenotype assay confirmed that reduced virulence correlated with decreased motility and biofilm formation ability. Additionally, three types of type VI secretion system, namely T6SS1, T6SS2, and T6SS3, were identified in V. harveyi strains, which can be classified into six, four, and 12 subtypes, respectively. In conclusion, the results indicated that the virulence level of V. harveyi is mainly determined by the above virulence genes, which may play vital roles in environmental adaptation for V. harveyi.

Horizontal Plasmid Transfer Promotes the Dissemination of Asian Acute Hepatopancreatic Necrosis Disease and Provides a Novel Mechanism for Genetic Exchange and Environmental Adaptation

Global outbreaks of shrimp acute hepatopancreatic necrosis disease (AHPND) caused by V. parahaemolyticus represent an urgent issue for the shrimp industry. This study revealed that the transmission mode of AHPND consists of two steps, the transregional dissemination of V. parahaemolyticus and the horizontal transfer of an AHPND-associated plasmid. Surprisingly, the introduction of the AHPND-associated plasmid also offers a novel mechanism of genetic exchange mediated by insertion sequences, and it improved the fitness of V. parahaemolyticus in a harsh environment. The results presented herein suggest that current shrimp farming practices promote genetic mixture between endemic and oceanic V. parahaemolyticus populations, which introduced the plasmid and accelerated bacterial adaptation by the acquisition of ecologically important functions. This entails a risk of the emergence of new virulent populations both for shrimp and humans. This study improves our understanding of the global dissemination of the AHPND-associated plasmid and highlights the urgent need to improve biosecurity for shrimp farming.

Vibrio parahaemolyticus is an important foodborne pathogen and has recently gained particular notoriety because it causes acute hepatopancreatic necrosis disease (AHPND) in shrimp, which has caused significant economic loss in the shrimp industry. Here, we report a whole-genome analysis of 233 V. parahaemolyticus strains isolated from humans, diseased shrimp, and environmental samples collected between 2008 and 2017, providing unprecedented insight into the historical spread of AHPND. The results show that V. parahaemolyticus is genetically diverse and can be divided into 84 sequence types (STs). However, genomic analysis of three STs of V. parahaemolyticus identified seven transmission routes in Asia since 1996, which promoted the transfer of an AHPND-associated plasmid. Notably, the insertion sequence (ISVal1) from the plasmid subsequently mediated the genetic exchange among V. parahaemolyticus STs and resulted in the deletion of an 11-kb region regulating cell mobility and the production of capsular polysaccharides. Phenotype assays confirmed that this deletion enhanced biofilm formation, providing a novel mechanism for environmental adaptation. We conclude that the transmission mode of AHPND consists of two steps, the transmission of V. parahaemolyticus and the subsequent horizontal transfer of the AHPND-associated plasmid. This plasmid allows ISVal1 to mediate genetic exchange and improve pathogen fitness in shrimp ponds. Current shrimp farming practices promoted such genetic exchanges, which highlighted a risk of the emergence of new virulent populations, with potentially devastating consequences for both aquaculture and human health. This study addressed the basic questions regarding the transmission mechanism of AHPND and provided novel insights into shrimp and human disease management.

IMPORTANCE Global outbreaks of shrimp acute hepatopancreatic necrosis disease (AHPND) caused by V. parahaemolyticus represent an urgent issue for the shrimp industry. This study revealed that the transmission mode of AHPND consists of two steps, the transregional dissemination of V. parahaemolyticus and the horizontal transfer of an AHPND-associated plasmid. Surprisingly, the introduction of the AHPND-associated plasmid also offers a novel mechanism of genetic exchange mediated by insertion sequences, and it improved the fitness of V. parahaemolyticus in a harsh environment. The results presented herein suggest that current shrimp farming practices promote genetic mixture between endemic and oceanic V. parahaemolyticus populations, which introduced the plasmid and accelerated bacterial adaptation by the acquisition of ecologically important functions. This entails a risk of the emergence of new virulent populations both for shrimp and humans. This study improves our understanding of the global dissemination of the AHPND-associated plasmid and highlights the urgent need to improve biosecurity for shrimp farming.

The Complete Genome of an Endogenous Nimavirus (Nimav-1_LVa) From the Pacific Whiteleg Shrimp Penaeus (Litopenaeus) Vannamei

White spot syndrome virus (WSSV), the lone virus of the genus Whispovirus under the family Nimaviridae, is one of the most devastating viruses affecting the shrimp farming industry. Knowledge about this virus, in particular, its evolution history, has been limited, partly due to its large genome and the lack of other closely related free-living viruses for comparative studies. In this study, we reconstructed a full-length endogenous nimavirus consensus genome, Nimav-1_LVa (279,905 bp), in the genome sequence of Penaeus (Litopenaeus) vannamei breed Kehai No. 1 (ASM378908v1). This endogenous virus seemed to insert exclusively into the telomeric pentanucleotide microsatellite (TAACC/GGTTA)n. It encoded 117 putative genes, with some containing introns, such as g012 (inhibitor of apoptosis, IAP), g046 (crustacean hyperglycemic hormone, CHH), g155 (innexin), g158 (Bax inhibitor 1 like). More than a dozen Nimav-1_LVa genes are involved in the pathogen-host interactions. We hypothesized that g046, g155, g158, and g227 (semaphorin 1A like) were recruited host genes for their roles in immune regulation. Sequence analysis indicated that a total of 43 WSSV genes belonged to the ancestral/core nimavirus gene set, including four genes reported in this study: wsv112 (dUTPase), wsv206, wsv226, and wsv308 (nucleocapsid protein). The availability of the Nimav-1_LVa sequence would help understand the genetic diversity, epidemiology, evolution, and virulence of WSSV.