Development and validation of a predictive model for the growth of Vibrio parahaemolyticus in post-harvest shellstock oysters

Information is limited about the growth and survival of naturally-occurring Vibrio parahaemolyticus in live oysters under commercially relevant storage conditions harvested from different regions and in different oyster species. This study produced a predictive model for the growth of naturally-occurring V. parahaemolyticus in live Eastern oysters (Crassostrea virginica) harvested from the Chesapeake Bay, MD, USA and stored at 5–30 °C until oysters gapped. The model was validated with model-independent data collected from Eastern oysters harvested from the Chesapeake Bay and Mobile Bay, AL, USA and Asian (C. ariakensis) oysters from the Chesapeake Bay, VA, USA. The effect of harvest season, region and water condition on growth rate (GR) was also tested. At each time interval, two samples consisting of six oysters each were analyzed by a direct-plating method for total V. parahaemolyticus. The Baranyi D-model was fitted to the total V. parahaemolyticus growth and survival data. A secondary model was produced using the square root model. V. parahaemolyticus slowly inactivated at 5 and 10 °C with average rates of − 0.002 and − 0.001 log cfu/h, respectively. The average GRs at 15, 20, 25, and 30 °C were 0.038, 0.082, 0.228, and 0.219 log cfu/h, respectively. The bias and accuracy factors of the secondary model for model-independent data were 1.36 and 1.46 for Eastern oysters from Mobile Bay and the Chesapeake Bay, respectively. V. parahaemolyticus GRs were markedly lower in Asian oysters. Harvest temperature, salinity, region and season had no effect on GRs. The observed GRs were less than those predicted by the U.S. Food and Drug Administration's V. parahaemolyticus quantitative risk assessment.

► A predictive model for the growth of naturally-occurring V. parahaemolyticus in live oysters was developed.
► V. parahaemolyticus was slowly inactivated at 5 and 10 °C.
► The average growth rates at 15, 20, 25, and 30 °C were 0.038, 0.082, 0.228, and 0.219 log cfu/h, respectively.
► Harvest temperature, salinity, region and season had no effect on growth rates.