Skip to main content

Bacteriostatic effects of benzyl isothiocyanate on Vibrio parahaemolyticus: Transcriptomic analysis and morphological verification



Foodborne illness caused by Vibrio parahaemolyticus (V. parahaemolyticus) is generally associated with the consumption of seafood. Fish and other seafood can be contaminated with V. parahaemolyticus, natural inhabitants of the marine, estuarine, and freshwater environment. In this study, the antibacterial activities of benzyl isothiocyanate (BITC) against V. parahaemolyticus were investigated by both transcriptomic analysis and morphological verification.


Treatment with 1/8 minimum inhibitory concentration (1/8 MIC) BITC resulted in 234 upregulated genes and 273 downregulated genes. The results validated by quantitative real-time polymerase chain reaction (qRT-PCR) revealed that the relative expression levels of the six genes VP0820, VP0548, VP2233, VPA2362, fliA and fliG were only 31.0%, 31.1%, 55.8%, 57.0%, 75.3%, and 79.9% of the control group, respectively. Among them, genes VP2233, fliA and fliG are related to flagella and VP2362 can regulate a protein relevant to biofilm formation. Morphologically, we verified that the swimming diffusion diameter of V. parahaemolyticus was significantly reduced by 14.9% by bacterial swimming ability, and biofilm formation was significantly inhibited by treatment with 1/8 MIC BITC by crystal violet quantification assay.


These results indicated that 1/8 MIC BITC had antibacterial effect on V. parahaemolyticus by inhibiting virulence gene expression related to flagella and biofilm.


Vibrio parahaemolyticus is a moderately halophilic gram-negative bacterium mainly in the form of rods, arcs, etc., without capsules and spores, and has become the main reason of aquatic product poisoning worldwide since it was identified in 1950 [1]. The colonies of V. parahaemolyticus CGMCC 1.1614, known as ATCC 33,845, are round, smooth, complete and cream colored. It can grow at temperature of 5 to 44 °C, and survive in the pH range from 4.8 to 11.0. It is common in coastal estuaries and marine environments and is usually associated with seawater and marine organisms such as zooplankton, plankton, mollusks and shellfish, fish and crabs [1]. Bacterial gastroenteritis related to seafood is mainly caused by direct or cross-contamination by V. parahaemolyticus [1]. Poisoning usually manifests as diarrhea, headache, vomiting, nausea, abdominal cramps and other symptoms. Therefore, it is necessary to take various measures to reduce the contamination of V. parahaemolyticus.

Antibiotics have long been applied to prevent and control marine pollution, but their long-term use has led to bacterial resistance to antibiotics and even the induction of aplastic anemia in human [2]. The hydrolysis of glucosinolate can produce isothiocyanates (ITCs), which are organic sulfides in cruciferous plants. ITCs have shown various beneficial effects, including antibacterial, antiviral, anti-inflammatory, anticancer, neuroprotective, chemical prevention and anti-parasitic properties. Benzyl isothiocyanate (BITC) is a type of ITC strong antibacterial capacity. Studies have found that BITC indicates antifungal activity and can effectively inhibit various molds [3] and bacteria, such as Campylobacter jejuni [4], Salmonella typhimurium [5], Staphylococcus aureus [6]. Studies have indicated the antibacterial mechanism of BITC. For example, BITC can affect the biofilm integrity, bacterial morphology, and membrane potential of Pseudomonas aeruginosa, E. coli and S. aureus [7]. However, knowledge about the antibacterial effects of BITC on V. parahaemolyticus is limited.

Scholars have conducted studies on the antibacterial mechanisms of natural products in V. parahaemolyticus. The major constituent of black seed volatile oil could affect quorum sensing, biofilm, and virulence-associated genes to reduce the virulence of V. parahaemolyticus [8]. Dihydromyricetin is the main bioactive component of Ampelopsis grossedentata, which can inhibit V. parahaemolyticus by decreasing the activity of proline dehydrogenase and inducing an increase in cell injury, proline content and cell surface hydrophobicity [9]. Banu et al. [10] studied the effect of essential oil from Cinnamomum tamala on V. parahaemolyticus and found that the virulence factors were regulated by polysaccharides, cytotoxins, flagella, lipopolysaccharides and biofilms. Song et al. [11] used the transcriptome to study the antibacterial effect of 1/4 MIC BITC on V. parahaemolyticus. However, no study has reported the bacteriostatic effect of BITC on V. parahaemolyticus by combining transcriptional level and morphological characteristics. Therefore, the antibacterial mechanism of BITC against V. parahaemolyticus can be studied from the genetic level and the destruction of the bacterial membrane.

In the current study, we determined the antibacterial effect of BITC at a 1/8 MIC subinhibitory concentration on V. parahaemolyticus. The differentially expressed genes (DEGs) were detected by RNA-seq, from which we screened out genes related to virulence and verified their expression by qRT-PCR. Morphologically, the influence of BITC on V. parahaemolyticus motility and biofilm formation was also verified.


Bacterial strain and culture

The V. parahaemolyticus CGMCC 1.1614 (tdh + , tlh + , trh -) was purchased from China General Microbiological Culture Collection Center. The strain was stored in physiological saline with 10% glycerol at − 80 °C. 100 µL of frozen bacteria solution was added into 10 mL tryptone soy broth medium with 3% sodium chloride (3% NaCl-TSB) liquid medium and cultured at 37 °C overnight. The cultured bacterial cells were streaked onto tryptone soy agar medium with 3% sodium chloride (3% NaCl-TSA) medium and incubated at 37 °C. Next, the activated single colony was inoculated into 3% NaCl-TSB liquid medium and cultured with shaking at 37 °C.

Antimicrobial tests

BITC was purchased from Sigma-Aldrich (CAS: 622–78-6), and its Flavor and Extract Manufacturing Association (FEMA) number is 4428. The MIC of BITC of V. parahaemolyticus was determined to be 9.54 µmol/L by the broth microdilution method. 1/8 MIC BITC was added to 100 mL 3% NaCl-TSB liquid culture medium with V. parahaemolyticus in log phase. Bacterial solution without BITC was incubated at 37 °C on a shaker (150 r/min) for 6 h. Bacterial solution without BITC was used as a control.

Enrichment and sequencing of RNA

V. parahaemolyticus was treated with 1/8 MIC BITC (E_BITC) or without BITC (C_BITC) for 6 h. Total RNA was extracted by the RNAprep Pure Cell/Bacterial Kit (Tiangen Biotech, Beijing, China) as recommended by the manufacturer. The preparation of sequencing library used the Illumina's NEBNext UltraTM Directed RNA Library Preparation Kit (NEB, Ipswich, MA, USA), and then added the index codes. All the samples were evaluated for product and library quality using an Agilent Bioanalyzer 2100 system (G2939B; Agilent Technologies, Palo Alto, CA, USA). The recount data obtained in the gene expression level analysis by DEGseq 2 software were used to analyze and screen the DEGs, wherein the screening standard was p < 0.05. DEGs is analyzed using volcanic plots, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The statistical enrichment of DEGs in the KEGG pathway was analyzed by KOBAS software.

qRT-PCR verification of the RNA-seq results

We used qRT-PCR to investigate the differential gene expression between different groups with 16S rRNA as the reference gene. The annealing temperature was set at 58 °C and the obtainment of melting curve was 60 to 95 °C. Additionally, the relative expression was calculated by 2−∆∆Ct method. The specific primer sequences for qRT-PCR are listed in Additional file 1. The genomic DNA was removed from the total RNA as recommended by instructions. The cDNA templates were reverse transcribed using the PrimeScript™ RT Reagent kit with gDNA Eraser (TaKaRa, Dalian, China). qRT-PCR was performed by using SYBR®Premix Ex Taq™II (TliRNaseH Plus) (Takara, Dalian, China) as recommended by the manufacturer.

Mobility measurement

The mobility of V. parahaemolyticus was determined based on the method of Butler et al. [12]. One milliliter of the log phase of V. parahaemolyticus was added to 99 mL of fresh 3% NaCl-TSB liquid medium, followed by incubation for 6 h. Next, 15 mL of swimming medium was sterilized and cooled to approximately 45 °C, then added 500 μL of 1/8 MIC BITC stock solution or physiological salin as the experimental group or control group, respectively. Thereafter, 3 μL of bacterial solution was added dropwise to the plate, and then the diameter of the inhibition zone was counted after incubation at 37 °C for 12 h. The average diameters of inhibition zone of the control group and experimental group were obtained to show the antibacterial effect of BITC on V. parahaemolyticus mobility. The data comparison of each group was analyzed by student’s t test. P < 0.05 indicated a statistically significant difference.

Biofilm measurement

The biofilm formation of V. parahaemolyticus affected by BITC was investigated according to the method of Pratt and Kolter [13]. Three-percent NaCl-TSB liquid medium, bacterial solution and 1/8 MIC BITC were added to a 96-well plate at ratio of 8:1:1, followed by incubating at 25 °C to form the biofilm. 3% NaCl-TSB liquid medium with and without bacterial culture was used as the controls. After treatment with phosphate buffer saline (PBS), anhydrous methanol, crystal violet and glacial acetic acid, the absorbance was measured at 590 nm using a microplate reader (Spectra Max M2, Molecular Devices, CA, USA). The experiment was repeated three times, and each group included three parallel determinations. The data comparison of two groups was analyzed by student’s t test. P < 0.05 indicated a statistically significant difference.


Transcriptome results

In our previous studies, Song et al. [11] determined the sensitivity of V. parahaemolyticus CGMCC 1.1614 to BITC, and the MIC was 9.54 μmol/L. Transcriptome analysis can reveal the underlying molecular mechanisms by which food additives act on pathogens. We performed transcriptome sequencing by the Illumina Hi Seq TM2500 high-throughput sequencing platform, generating two transcriptome databases for the control group without BITC (C_BITC) and the experimental group with 1/8 MIC BITC treatment (E_BITC) to understand the bacteriostatic effect of BITC on V. parahaemolyticus. The RNA-seq results were uploaded to the Gene Expression Omnibus (GEO) database, numbered GSE152671. The total quantity of sequencing data we obtained was 9.52 Gb, the error rate of the single-base position of all samples was less than 1%, and Q20 was greater than 80%, which were up to the standards of sequencing quality control (see Additional file 2).

The volcano map shows the distribution of DEGs. According to Fig. 1, E_BITC obtained 507 DEGs compared with C_BITC, among which 234 were upregulated and 273 were downregulated.

Fig. 1
figure 1

Volcano plot of differentially expressed genes (DEGs) in Vibrio parahaemolyticus treated with 1/8 MIC BITC

The top 10% genes with the most significantly differential expression were shown in the Table 1, of which 19 were up-regulated and 31 were down-regulated. It was shown that BITC regulated 14 enzyme-related genes, including cryptic beta-D-galactosidase subunit alpha, membrane-bound lytic murein transglycosylase D, orotate phosphoribosyltransferase, long-chain-fatty-acid-CoA ligase, manganese-dependent inorganic pyrophosphatase, oxaloacetate decarboxylase subunit gamma, carbamoyl phosphate synthase small subunit, GTP cyclohydrolase I, ribonuclease PH, quinolinate phosphoribosyltransferase, 5'-deoxynucleotidase, short chain dehydrogenase, glycerol kinase and cytidylate kinase. Oxaloacetate decarboxylase α/β/γ complex is a membrane-bound enzyme complex in Vibrio cholerae, and its assembly is affected the oxaloacetate decarboxylase subunit gamma regulated by VP2545 [14]. The influence of BITC on the activities of these enzyme-related genes should be discussed in the further study. BITC also showed effects on the other genes related to bacterial virulence. For examples, BITC upregulated molecular chaperone DnaK, which is related to protein folding [15]. BITC also affected the expression of transfer protein genes such as uhpT and VP2869, which regulate sugar phosphate antiporter and sodium/solute symporter, respectively (Table 1). Membrane associated protein genes such as VP1286 and VP1091, which regulate integral membrane protein and transmembrane protein, respectively, affecting cell membrane permeability, biofilm formation and septum formation (Table 1). BITC also regulated some functional proteins, such as sodium/solute symporter, carbon starvation protein A, integral membrane protein, heavy metal membrane efflux protein, RhlE protein, sugar phosphate antiporter, transmembrane protein affecting septum formation and cell membrane permeability, lipoprotein, heat shock protein 90, phosphate ABC transporter ATP-binding protein, integral membrane protein transporter, SpoOM-like protein, 30S ribosomal protein S21, thiamine biosynthesis protein ThiC, NadC family protein and nucleoid occlusion protein. Among them VP0821-regulated heat shock protein 90 is important for the virulence and spread of protozoan parasites [16].

Table 1 Top 10% genes with significant differential expression identified in E_BITC and C_BITC from RNA sequencing

Besides, rpsU-encoded ribosomal protein S21 can affect the motility and biofilm formation of Bacillus subtilis [17]. The down-regulated genes VP1892 and FliG are related to bacterial chemotaxis (see Additional file 3). FliG regulates flagellar motor switch protein G (Table 2) and VP1892 regulates methyl-accepting chemotaxis protein. VP2826 and VPA1006 were down-regulated, regulating transporter and LysR family transcriptional regulator, respectively. According to the objective of present study, we focused on the antibacterial mechanism of BITC by screening differentially expressed genes related to virulence, as shown in Table 2 (VP0820, VP0548, VP2233, VPA2362, fliA and fliG), in which VP0820, VP0548 can regulate ToxR, thereby affecting the virulence of V. parahaemolyticus; VP2233, fliA and fliG genes are related to flagella, thereby affecting motility, and VP2362-encoded outer membrane protein that is related to the biofilm formation.

Table 2 Differentially expressed genes relevant to virulence from RNA sequencing

GO analysis was performed on the DEGs in V. parahaemolyticus treated with 1/8 MIC BITC. We found that these genes were annotated into 1679 GO terms, and the enrichment was selected from each group. Thirty significant GO terms were identified, with three main domains, molecular functions (12 subclasses), biological processes (10 subclasses) and cellular components (8 subclasses). Among these, 891 were annotated to biological processes, 227 were annotated to cellular components, and 561 were annotated to molecular function. Figure 2a is the most significantly enriched composition selected from each domain. Among them, the macromolecular complex in the cellular component had 23 genes upregulated and 27 genes downregulated. Six upregulated genes rplC, rplD, rplW, rplB, rpsS and rplV, and 11 downregulated genes rpsR, rpsU, rpsF, rpsI, rplM, rplN, RplU, VP1210, rpmA, rpmE and rpmG all regulate ribosome-related proteins (see Additional file 4). Studies have shown that ribosomal proteins can affect bacterial motility and biofilm formation [17]. Downregulated genes VP1892, VP1904, VP2629 and VPA1000 all regulate the methyl-accepting chemotaxis protein (see Additional file 4), it is related to flagellar movement [18]. Both cilium and motile cilium had significant changes in 5 genes, the upregulated gene VP1392 and downregulated genes VP0246, VP0417, alaS and VP2629 respectively affected the synthesis of cilia-related ClpA/B-type protease, hypothetical protein, hypothetical protein, alanyl-tRNA synthetase, and methyl-accepting chemotaxis protein to influence motility (see Additional file 4). Six genes were significant changed in the groups of intrinsic component of organelle membrane and integral component of organelle membrane. The upregulated gene metF regulates 5, 10-methylenetetrahydrofolate reductase, and downregulated genes VP0388, VP0246, VP0470, VP0939, and VPA1370 regulate type I restriction enzyme M protein, hypothetical protein, carbamoyl phosphate synthase small subunit, hypothetical protein and hypothetical protein, respectively (see Additional file 4). Three genes VP0295, VPA1128, and VPA1735 were upregulated, while eight genes VP1092, VP1256, VP1741, nhaB, VP2351, VP2545, VP2778, and VP2826 were downregulated in sodium ion transport of biological processes. They regulated multiple proteins related to sodium ion transport, such as sodium/sulfate symporter, acyl-CoA carboxylase alpha chain, acridine efflux pump, NadC family protein, sodium/alanine symporter, sodium/proton antiporter, Na( +)-translocating NADH-quinone reductase subunit A, oxaloacetate decarboxylase subunit gamma, FKBP-type peptidylprolyl isomerase and transporter (see Additional file 4). In the molecular function, oxidoreductase activity had 32 upregulated genes and 44 downregulated genes, transferase activity, transferring one-carbon groups had 4 upregulated genes and 19 downregulated genes. Multiple genes related to oxidoreductase such as VP0068, VP0235, VP1017, VPA0278, phhA, VP0442, sdhC, VP1710, VP2014, and VPA0566, which regulate glutathione reductase, epimerase/dehydratase, arginyl-tRNA-protein transferase, isopentenyl pyrophosphate isomerase, phenylalanine 4-monooxygenase, ubiquinol-cytochrome c reductase, cytochrome b, succinate dehydrogenase cytochrome b556 large membrane subunit, glucose-6-phosphate 1-dehydrogenase, tetrathionate reductase subunit A and alcohol dehydrogenase, respectively (see Additional file 4). Among the DEGs in the term of transferase activity, transferring one-carbon groups, there are multiple methyltransferase-related genes such as VPA0046, ubiE, VP0594, VP0954, yebU, VP1933, VP2477, and gidB, which regulate methylated-DNA–protein -cysteine S-methyltransferase, ubiquinone/menaquinone biosynthesis methyltransferase, RNA methyltransferase, 16S rRNA (cytosine(1407)-C(5))-methyltransferase RsmF, 3-demethylubiquinone-9 3-methyltransferase, 16S rRNA methyltransferase and 16S rRNA methyltransferase GidB (see Additional file 4). Some studies have confirmed that methyltransferase can affect bacterial motility and biofilm formation [19]. These terms might be related to the motility and biofilm formation in V. parahaemolyticus.

Fig. 2
figure 2

GO analysis (a) and KEGG analysis (b) of differentially expressed genes (DEGs). a Thirty significantly enriched GO terms are shown. *Means significantly enriched GO terms. b The 20 special KEGG pathways above include 507 single genes. The closer the q-val is to zero, the higher is the degree of gene enrichment

The analysis of KEGG enrichment was conducted following GO analysis. In total, 507 DEGs were enriched into 78 pathways compared with the control group. Bacterial chemotaxis pathways can affect the direction of bacterial movement by adjusting the rotation direction of flagella movement. Bacterial chemotaxis pathways included 7 differentially expressed genes, of which VP1628, VP1892, VP1904, VP2629, and VPA1000 regulate methyl-accepting chemotaxis protein, while VP2230 and fliG regulate chemotaxis protein CheZ and flagellar motor switch protein G, respectively (see Additional file 3). Flagellar assembly pathways included three DEGs flgF, fliI, and fliG, which are related to flagellar basal body rod protein FlgF, flagellar-specific ATP synthase and flagellar motor switch protein G, respectively (see Additional file 3). Figure 2b shows 20 pathways with the most enrichment, such as ABC transporters, biosynthesis of secondary metabolites, carbon metabolism, microbial metabolism in diverse environments and metabolic pathways.

qRT-PCR verification of the RNA-seq results

qRT-PCR is a common and important technical method to study the gene expression levels. The DEGs related to virulence were screened by RNA-seq, and the reliability of the data was verified by qRT-PCR. According to RNA-seq, six genes related to virulence with obviously decreased expression were screened (Table 2). The expression of six virulence-related genes was investigated in V. parahaemolyticus treated with 1/8 MIC BITC by qRT-PCR. The virulence-related genes VP0820, VP0548, VP2233, VPA2362, fliA, and fliG were downregulated (Fig. 3). The relative expression of these genes was significantly reduced with values of 31.0%, 31.1%, 55.8%, 57.0%, 75.3%, and 79.9% of the control group, respectively. These results verified that BITC at 1/8 MIC could effectively inhibit virulence-related gene expression of V. parahaemolyticus.

Fig. 3
figure 3

Inhibition of the gene expressions related to virulrnce in V. parahaemolyticus by 1/8 MIC BITC. VP0820 and VP0548 encode the ToxR protein gene, VP2233, fliA and fliG are related genes of flagellin synthesis, and VPA2362 encodes the outer membrane protein gene. The data in the figure were derived from the average of three parallel experiments. *p < 0.05 indicates significant differences

Mobility measurement

The motility of V. parahaemolyticus is affected by the flagella of bacterium, which are connected with the virulence of bacteria. According to the transcriptome results, the expression of flagella-related genes VP2233, fliA, and fliG was downregulated. By detecting the movement of V. parahaemolyticus, we verified the antibacterial effect of BITC on V. parahaemolyticus morphologically. The motility of V. parahaemolyticus is related to the flagella of bacterium. As shown in Fig. 4, BITC at 1/8 MIC significantly suppressed the mobility of V. parahaemolyticus. BITC significantly reduced the diameter of swimming diffusion by 14.9% (p < 0.05) compared with that of control group. This finding agreed with the downregulation of VP2233, fliA, and fliG gene expression in the transcriptome results (Table 2). Therefore, swimming is related to the bacteriostatic effect of BITC in V. parahaemolyticus.

Fig. 4
figure 4

Inhibitory influences of BITC on the swimming ability of V. parahaemolyticus. The data in the figure were derived from the average of three parallel

Biofilm measurement

The formation of biofilm is related to the virulence of bacteria. In most cases, bacteria adhere to the surface of object during formation of biofilms to maintain survival, causing food safety problems. The mechanism of BITC in food-borne pathogenic bacteria can be investigated by exploring the biofilm changes of V. parahaemolyticus after BITC treatment. VP2362 encodes an outer membrane protein, which is related to the biofilm formation of V. parahaemolyticus. According to the transcriptome results, the relative expression of this gene was downregulated by treatment with BITC. We expected to verify the transcriptome results by measuring the biofilm formation of V. parahaemolyticus and understand the mechanism of BITC effects. The formation of biofilm of V. parahaemolyticus was significantly inhibited by 1/8 MIC BITC treatment (p < 0.05), and BITC functioned to scavenge biofilm (Fig. 5). These results were consistent with the downregulation of VP2362 gene expression (see Additional file 4).

Fig. 5
figure 5

Inhibitory influences of BITC on the biofilm formation of V. parahaemolyticus. The data in the figure were derived from the average of three parallel


Transcriptome sequencing technology is an important technical approach to screen DEGs and determine gene functions. It has advantages in terms of time, data yield, cost, coverage, and accuracy of data [20]. It is also possible to obtain gene-related biological function data by analyzing multilevel gene regulation expression. Studies based on prokaryotic transcriptome sequencing have included E. coli [21], S. aureus [22], and Listeria monocytogenes [23]. Moreover, the difference in gene expression of V. parahaemolyticus in seawater has been studied under different oxygen environments [24]. Therefore, RNA sequencing was applied in this study to understand the corresponding changes of V. parahaemolyticus CGMCC 1.1614 under the action of 1/8 MIC BITC. According to the research purpose, we screened the DEGs related to virulence at the molecular level and then verified the gene expression using qRT-PCR. From a molecular biology perspective, the impact of BITC on both the virulence and morphology of V. parahaemolyticus was also investigated. Studies have shown that the motility of bacteria will be inhibited when sodium ion transport is disturbed [25]. Oxidoreductase and transferase are very important in the formation of biofilms [26, 27]. RNA sequencing results of Salmonella Enteritidis under acid stress have shown that the terms of macromolecular complexes, transferase activity, and transferring one-carbon groups are the most highly represented [28]. Combined with the GO analysis results, multiple pathways related to motility and biofilm were significantly enriched (sodium ion transport, oxidoreductase activity, transferase activity, transferring one-carbon groups and macromolecular complex). Therefore, it is necessary to investigate the influence of BITC on V. parahaemolyticus of motility and biofilm formation.

Recent findings have confirmed that flagella are adhesive and invasive as a potential class of virulence factors [29]. Bacterial motility is closely related to flagella, and V. parahaemolyticus has two different types of flagella to accommodate life in different situations [24]. Studies have shown that bacterial motility is also inhibited when flagella are affected. When histone-like nucleoid structuring protein (H-NS) acted on the flagellin lateral flagellar A (lafA), the swarming motilities of Vibrio hemolyticus [30] was inhibited. Quinazoline-2, 4-diamino analogs inhibit the aggregation of polar flagella and movement of V. cholerae [31], and phenethyl isothiocyanate (PEITC) can reduce the migration capacity of E. coli [32]. In this study, the transcriptome results indicated that the flagella-related gene expression of VP2233, VP2232, and VP2248 were significantly inhibited by 1/8 MIC BITC treatment. Additionally, the morphological verification results of Fig. 4 were similar to the transcriptome results, further proving that the flagella and bacterial motility are related to the bacteriostatic effects of BITC on V. parahaemolyticus.

The antimicrobial effects of ITCs are thought to be associated with damage to bacterial membrane integrity [33]. ITCs, including BITC, prevent the formation of biofilms in various bacteria such as P. aeruginosa, L. monocytogenes, S. aureus, and E. coli [32]. VP2362 is a gene encoding an outer membrane protein, which is connected with the biofilm formation of V. parahaemolyticus. By BITC treatment, gene expression level of VP2362 was significantly downregulated compared with that of control group. Previous findings have confirmed that the cell surface integrity, neutralizing host defense mechanisms, cell adhesion and invasion, and inhibition of the complement system are all affected by outer membrane proteins [34]. Therefore, the outer membrane protein is considered to be a virulence factor. The expression of membrane genes was significantly downregulated after BITC treatment at 1/8 MIC (p < 0.05). We found that the formation of V. parahaemolyticus of biofilm was inhibited by BITC treatment, which was consistent with the gene expression results.

In this paper, the bacteriostatic effect of BITC on V. parahaemolyticus was reported by RNA-seq technology from the gene level, showing multiple pathways of secondary metabolites. Various virulence genes were inhibited that were further verified morphologically. Further studies will focus on other important DEGs and the corresponding proteomics studies to obtain more comprehensive mechanism for bacteriostatic effects of BITC.


In this paper, the bacteriostatic effect of BITC on V. parahaemolyticus was reported by RNA-seq technology from the gene level, showing multiple pathways of secondary metabolites. Various virulence genes were inhibited that were further verified morphologically. Further studies will focus on other important DEGs and the corresponding proteomics studies to obtain more comprehensive mechanism for bacteriostatic effects of BITC.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.



Benzyl isothiocyanate


Experimental group


Control group


Minimum inhibitory concentration


Quantitative real-time PCR




Differentially expressed genes


Gene ontology


Kyoto Encyclopedia of Genes and Genomes


Phosphate buffer saline


Gene expression omnibus


  1. 1.

    Su YC, Liu C. Vibrio parahaemolyticus: a concern of seafood safety. Food Microbiol. 2007;24(6):549–58.

    Article  Google Scholar 

  2. 2.

    Zhu M, Li W, Dong X, Chen Y, Lu Y, Lin B, et al. Benzyl-isothiocyanate induces apoptosis and inhibits migration and invasion of hepatocellular carcinoma cells in vitro. J Cancer. 2016;8(2):240–8.

    Article  Google Scholar 

  3. 3.

    Federica S, Juan MQ, Fernando BL, Jordi M, Mónica FF, Giuseppe M. Shelf life improvement of the loaf bread using allyl, phenyl and benzyl isothiocyanates against Aspergillus parasiticus. Lwt Food Sci Technol. 2017;78:208–14.

    Article  Google Scholar 

  4. 4.

    Virginie D, Bachar A, Gwennola E, Marion T, Albert R, Odile T, et al. Antimicrobial activities of isothiocyanates against Campylobacter jejuni Isolates. Front Cell Infect Microbiol. 2012;2(53):53.

    Google Scholar 

  5. 5.

    Niu TX, Wang XN, Wu HY, Bi JR, Zhang GL. Transcriptomic analysis, motility and biofilm formation characteristics of Salmonella typhimurium exposed to benzyl isothiocyanate treatment. Int J Mol Sci. 2020;21(3):1025.

    CAS  Article  Google Scholar 

  6. 6.

    Wang X, Wu H, Niu T, Bi J, Zhang G. Downregulated expression of virulence factors induced by benzyl isothiocyanate in Staphylococcus Aureus: a transcriptomic analysis. Int J Mol Sci. 2019;20(21):5441.

    CAS  Article  Google Scholar 

  7. 7.

    Anabela B, Ana CA, Carla F, Maria JS, Lúcia CS, Manuel S. Antibacterial activity and mode of action of selected glucosinolate hydrolysis products against bacterial pathogens. J Food Sci Technol. 2015;52:4737–48.

    Article  Google Scholar 

  8. 8.

    Guo D, Yang Z, Zheng X, Kang S, Xia X. Thymoquinone inhibits biofilm formation and attachment-invasion in host cells of Vibrio parahaemolyticus. Foodborne Pathog Dis. 2019;16(10):671–8.

    CAS  Article  Google Scholar 

  9. 9.

    Ding L, Xiao S, Dan L, Pang W. Effect of dihydromyricetin on proline metabolism of Vibrio parahaemolyticus: Inhibitory mechanism and interaction with molecular docking simulation. J Food Biochem. 2018; 42(1):e12463.

  10. 10.

    Banu SF, Rubini D, Murugan R, Vadivel V, Gowrishankar S, Pandian SK, et al. Exploring the antivirulent and sea food preservation efficacy of essential oil combined with DNase on Vibrio parahaemolyticus. LWT- Food Sci Technol. 2018;95:107–15.

    Article  Google Scholar 

  11. 11.

    Song J, Hou HM, Wu HY, Li KX, Wang Y, Zhou QQ, et al. Transcriptomic analysis of vibrio parahaemolyticus reveals different virulence gene expression in response to benzyl isothiocyanate. Molecules. 2019;24(4):761.

    Article  Google Scholar 

  12. 12.

    Butler MT, Wang Q, Gadagkar HR. Cell density and mobility protect swarming bacteria against antibiotics. Proc Natl Acad Sci USA. 2010;107(8):3776–81.

    CAS  Article  Google Scholar 

  13. 13.

    Pratt LA, Kolter R. Genetic analysis of Escherichia coli biofilm formation: roles of flagella, motility, chemotaxis and type I pili. Mol Microbiol. 2010;30(2):285–93.

    Article  Google Scholar 

  14. 14.

    Dahinden P, Pos KM, Dimroth P. Identification of a domain in the α-subunit of the oxaloacetate decarboxylase Na+ pump that accomplishes complex formation with the γ-subunit. FEBS J. 2010;272(3):846–55.

    Article  Google Scholar 

  15. 15.

    Langer T, Lu C, Echols H, Flanagan J, Hayer MK, Hartl FU. Successive action of DnaK, DnaJ and GroEL along the pathway of chaperone-mediated protein folding. Nature. 1992;356(6371):683–9.

    CAS  Article  Google Scholar 

  16. 16.

    Meetali S, Shalini S, Alok B, Utpal T. Heat shock protein 90 regulates encystation in Entamoeba. Front Microbiol. 2015;6:1125.

    Google Scholar 

  17. 17.

    Takada H, Morita M, Shiwa Y, Sugimoto R, Yoshikawa H. Cell motility and biofilm formation in Bacillus subtilis are affected by the ribosomal proteins, S11 and S21. Biosci Biotechnol Biochem. 2014;78(5):898–907.

    CAS  Article  Google Scholar 

  18. 18.

    Engström P, Hazelbauer GL. Methyl-accepting chemotaxis proteins are distributed in the membrane independently from basal ends of bacterial flagella. BBA Biomembranes. 1982;686(1):19–26.

    Article  Google Scholar 

  19. 19.

    Marczak M, Żebracki K, Koper P, Turska-Szewczuk A, Skorupska A. Mgl2 is a hypothetical methyltransferase involved in exopolysaccharide production, biofilm formation, and motility in rhizobium leguminosarum bv. trifolii. Mol Plant Microbe Interact. 2019;32(7):899–911.

    CAS  Article  Google Scholar 

  20. 20.

    Bhardwaj J, Chauhan R, Swarnkar MK, Chahota RK, Singh AK, Shankar R, et al. Comprehensive transcriptomic study on horse gram (Macrotyloma uniflorum): De novo assembly, functional characterization and comparative analysis in relation to drought stress. BMC Genom. 2013;14:647.

    CAS  Article  Google Scholar 

  21. 21.

    Tao H, Christoph B, Craig R, Frederick RB, Tyrrell C. Functional genomics: expression analysis of Escherichia coli growing on minimal and rich media. J Bacteriol. 1999;181(20):6425.

    CAS  Article  Google Scholar 

  22. 22.

    Tan X, Nan Q, Wu C, Sheng J, Rui Y, Zheng B, et al. Transcriptome analysis of the biofilm formed by methicillin-susceptible. Sci Rep. 2015;5:11997.

    Article  Google Scholar 

  23. 23.

    Pieta L, Escudero F, Jacobus AP, Cheiran KP, Gross J, Moya M, et al. Comparative transcriptomic analysis of Listeria monocytogenes reveals upregulation of stress genes and downregulation of virulence genes in response to essential oil extracted from Baccharis psiadioides. Ann Microbiol. 2017;67:479–90.

    CAS  Article  Google Scholar 

  24. 24.

    Feng L, Zhang K, Gao M, Shi C, Han J. Inactivation of Vibrio parahaemolyticus by aqueous ozone. J Microbiol Biotechnol. 2018;28(8):1233–46.

    CAS  Article  Google Scholar 

  25. 25.

    Jaques S, Kim YK, Mccarter LL. Mutations conferring resistance to phenamil and amiloride, inhibitors of sodium-driven motility of Vibrio parahaemolyticus. Proc Natl Acad Sci USA. 1999;96(10):5740–5.

    CAS  Article  Google Scholar 

  26. 26.

    Naeem A, Fazle RS, Ute R, Mikael R, Michael H. Modulation of biofilm-formation in Salmonella enterica Serovar Typhimurium by the periplasmic DsbA/DsbB oxidoreductase system requires the GGDEF-EAL domain protein STM3615. PLoS ONE. 2014;9(8):e106095.

    Article  Google Scholar 

  27. 27.

    Siala W, Kucharíková S, Braem A, Vleugels J, Tulkens PM, Mingeot-Leclercq M-P, et al. The antifungal caspofungin increases fluoroquinolone activity against Staphylococcus aureus biofilms by inhibiting N-acetylglucosamine transferase. Nat Commun. 2016;7(1):13286.

    CAS  Article  Google Scholar 

  28. 28.

    Jia K, Wang G, Liang L, Wang M, Wang H, Xu X. Preliminary transcriptome analysis of mature biofilm and planktonic cells of Salmonella enteritidis exposure to acid stress. Front Microbiol. 2017;8:1861.

    Article  Google Scholar 

  29. 29.

    Haiko J, Westerlund-Wikström B. The role of the bacterial flagellum in adhesion and virulence. Biology. 2013;2(4):1242–67.

    Article  Google Scholar 

  30. 30.

    Wang Y, Yiquan Z, Yin Z, Wang J, Zhu Y, Peng H, et al. H-NS represses transcription of the flagellin gene lafA of lateral flagella in Vibrio parahaemolyticus. Can J Microbiol. 2017;64:69–74.

    CAS  Article  Google Scholar 

  31. 31.

    Rasmussen L, White EL, Pathak A, Ayala JC, Wang H, Wu JH, et al. A high-throughput screening assay for inhibitors of bacterial motility identifies a novel inhibitor of the Na+-driven flagellar motor and virulence gene expression in Vibrio cholerae. Antimicrob Agents Chemother. 2011;55(9):4134–43.

    CAS  Article  Google Scholar 

  32. 32.

    Abreu AC, Borges A, Mergulhão F, Simões M. Use of phenyl isothiocyanate for biofilm prevention and control. Int Biodeterior Biodegradation. 2014;86:34–41.

    CAS  Article  Google Scholar 

  33. 33.

    Abreu AC, Borges A, Saavedra MJ, Simões M. Antibacterial activity of phenyl isothiocyanate on Escherichia coli and Staphylococcus aureus. Med Chem. 2013;9:756–61.

    CAS  Article  Google Scholar 

  34. 34.

    Koebnik R, Locher KP, Gelder PV. Structure and function of bacterial outer membrane proteins: barrels in a nutshell. Mol Microbiol. 2010;37(2):239–53.

    Article  Google Scholar 

Download references


Not applicable.


This work was supported by the National Key R and D Program of China (No. 2019YFC1605902), the Liaoning Provincial Natural Science Foundation of China (No. 2019-MS-021), and the Program for Innovative Talents of Higher Learning Institutions of Liaoning (No. LR2019009).

Author information




GZ conceived and designed the project. GZ and JL drafted the manuscript. JL, KZ and JS performed the experiments. HW and JB analyzed data. HH modified the manuscript. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Gongliang Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Sequences of specific primers for qRT–PCR (doc).

Additional file 2.

Data of RNA sequencing (doc).

Additional file 3.

Differentially expressed genes related to bacterial motility from KEGG pathways (doc).

Additional file 4.

Differentially expressed genes from GO terms (doc).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, J., Zhang, K., Song, J. et al. Bacteriostatic effects of benzyl isothiocyanate on Vibrio parahaemolyticus: Transcriptomic analysis and morphological verification. BMC Biotechnol 21, 56 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Benzyl isothiocyanate
  • Vibrio parahaemolyticus
  • RNA-seq
  • Motility
  • Biofilm