Genome sequencing analysis and validation of infestation-related functional genes of Vibrio parahaemolyticus LG2206 isolated from the hepatopancreas of diseased mud crab (Scylla paramamosain) in South China

Vibrio parahaemolyticus (V. parahaemolyticus) is a major bacterial pathogen found in brackish environments, leading to disease outbreaks and great economic losses in the mud crab industry. This study investigated the molecular mechanism of V. parahaemolyticus infecting mud crabs through genome sequencing analysis, survival experiments, and the expression patterns of related functional genes. A strain of V. parahaemolyticus with high pathogenicity and lethality was isolated from diseased mud crab in South China. The genome sequencing results showed that the genome size of V. parahaemolyticus was a circular chromosome of 3,357,271 bp, with a GC content of 45 %, containing 2985 protein-coding genes, denoted as V. parahaemolyticus LG2206. Genome analysis data revealed that a total of 113 adherence coding genes were obtained, including 120 virulence factor coding genes, 37 type III secretion system (T3SS) coding genes, and 277 sequences of T3SS effectors. Survival experiments showed that the mortality was 20 % within 96 h in the 1 × 104 CFU/mL infection group, 90 % in the 3.2 × 105 CFU/mL treatment group, and 100 % in the 1 × 106 CFU/mL treatment group. The LD50 of V. parahaemolyticus LG2206 was determined as 4.6 × 104 CFU/mL. Six genes of znuA and fliD (flagellin encoding genes), yscE and yscR (T3SS encoding genes), and nfuA and htpX (virulence factor encoding genes) were selected and validated by quantitative real-time PCR analysis after infection with 4.6 × 104 CFU/mL of V. parahaemolyticus LG2206 for 96 h. The expression of the six genes exhibited a significant up-regulation trend at all tested time points. The results indicated that the infestation-related genes screened in the experiment play important roles in the infestation process. This study provides timely and effective information to further analyze the molecular mechanism of V. parahaemolyticus infection and develop comprehensive measures for disease prevention and control. •A high pathogenicity and lethality of V. parahaemolyticus was isolated, and name as LG2206.•The genome size of LG2206 was 3,357,271 bp, containing 2985 protein-coding genes.•The LG2206 codes 120 virulence factor, 37 type III secretion system (T3SS), and 277 sequences of T3SS effectors.•The LD50 of LG2206 was determined as 4.6 × 104 CFU/mL.