Abstract
Vibrio parahaemolyticus is a bacterial pathogen that becomes lethal to Penaeus shrimps when acquiring the pVA1-type plasmid carrying the PirABvp genes, causing acute hepatopancreatic necrosis disease (AHPND). This disease causes significant losses across the world, with outbreaks reported in Southeast Asia, Mexico, and South America. Virulence level and mortality differences have been reported in isolates from different locations, and whether this phenomenon is caused by plasmid-related elements or genomic-related elements from the bacteria remains unclear. Here, nine genomes of South American AHPND-causing V. parahaemolyticus (VPAHPND) isolates were assembled and analyzed using a comparative genomics approach at (i) whole-genome, (ii) secretion system, and (iii) plasmid level, and then included for a phylogenomic analysis with another 86 strains. Two main results were obtained from our analyses. First, all isolates contained pVA1-type plasmids harboring the toxin coding genes, and with high similarity with the prototypical sequence of Mexican-like origin, while phylogenomic analysis showed some level of heterogeneity with discrete clusters and wide diversity compared to other available genomes. Second, although a high genomic similarity was observed, variation in virulence genes and clusters was observed, which might be relevant in the expression of the disease. Overall, our results suggest that South American pathogenic isolates are derived from various genetic lineages which appear to have acquired the plasmid through horizontal gene transfer. Furthermore, pathogenicity seems to be a multifactorial trait where the degree of virulence could be altered by the presence or variations of several virulence factors. IMPORTANCE AHPND have caused losses of over $2.6 billion to the aquaculture industry around the world due to its high mortality rate in shrimp farming. The most common etiological agent is V. parahaemolyticus strains possessing the pVA1-type plasmid carrying the PirABvp toxin. Nevertheless, complete understanding of the role of genetic elements and their impact in the virulence of this pathogen remains unclear. In this work, we analyzed nine South American AHPND-causing V. parahaemolyticus isolates at a genomic level, and assessed their evolutionary relationship with other 86 strains. We found that all our isolates were highly similar and possessed the Mexican-type plasmid, but their genomic content did not cluster with other Mexican strains, but instead were spread across all isolates. These results suggest that South American VPAHPND have different genetic backgrounds, and probably proceed from diverse geographical locations, and acquire the pVA1-type plasmid via horizontal gene transfer at different times.