Vibrio parahaemolyticus has been identified as an emerging human pathogen worldwide with cases undergoing a global expansion over recent decades in phase with climate change. New Zealand had remained free of outbreaks until 2019, but different outbreaks have been reported consecutively since then. To provide new insights into the recent emergence of cases associated with outbreak clones over recent years, a comparative genomic study was carried out using a selection of clinical (mostly outbreak) and environmental isolates of V. parahaemolyticus obtained in New Zealand between 1973 and 2021. Among 151 isolates of clinical (n=60) and environmental (n=91) origin, 47 sequence types (STs) were identified, including 31 novel STs. The population of environmental isolates generated 30 novel STs, whereas only 1 novel ST (ST2658) was identified among the population of clinical isolates. The novel clinical ST was a single-locus variant of the pandemic ST36 strain, indicating further evolution of this pandemic strain. The environmental isolates exhibited a significant genetic heterogeneity compared to the clinical isolates. The whole-genome phylogeny separated the population of clinical isolates from their environmental counterparts, clearly indicating their distant genetic relatedness. In addition to differences in ancestral profiles and genetic relatedness, these two groups of isolates exhibited a profound difference in their virulence profiles. While the entire population of clinical isolates harboured the thermostable direct haemolysin (tdh) and/or the thermostable-related haemolysin (trh), only a few isolates of environmental origin possessed the same virulence genes. In contrast to tdh and trh, adhesin-encoding genes, vpadF and MSHA, showed a significantly (P<0.001) greater association with the environmental isolates compared to the clinical isolates. The effectors, VopQ, VPA0450 and VopS, which belong to T3SS1, were ubiquitous, being present in each isolate regardless of its origin. The effectors VopC and VopA, which belong to T3SS2, were rarely detected in any of the examined isolates. Our data indicate that the clinical and environmental isolates of V. parahaemolyticus from New Zealand differ in their population structures, ancestral profiles, genetic relatedness and virulence profiles. In addition, we identified numerous unique non-synonymous single-nucleotide polymorphisms (nsSNPs) in adhesins and effectors, exclusively associated with the clinical isolates tested, which may suggest a possible role of these mutations in the overall virulence of the clinical isolates.
Keywords: Vibrio parahaemolyticus; non-synonymous mutations; pandemic strains; population genomics; virulence; whole-genome sequencing.