Novel sandwich immunoassay detects a shrimp AHPND-causing binary PirABVp toxin produced by Vibrio parahaemolyticus

1 Introduction

Acute hepatopancreatic necrosis disease (AHPND or early mortality syndrome) is a bacterial disease that causes significant economic losses to the global shrimp aquaculture industry (Shinn et al., 2015). Since its first outbreak in China in 2009, AHPND has been reported in several countries of the Americas and Asia (Global Seafood Alliance, 2012; Tran et al., 2013; Joshi et al., 2014; de la Peña et al., 2015; Soto-Rodriguez et al., 2015; Tun et al., 2017; Eshik et al., 2018; Dhar et al., 2019). Moreover, AHPND was reported in Korea in 2016, causing serious economic losses (Han et al., 2020). The causative agent of AHPND has been confirmed to be Vibrio parahaemolyticus strains (Vp^AHPND), which express a binary PirAB^Vp toxin homologous to the Photorhabdus insect-related toxin (Han et al., 2015; Lee et al., 2015; Sirikharin et al., 2015). The production of the PirAB^Vp toxin in Vp^AHPND is responsible for its cytotoxic effect on the epithelial cells of the shrimp hepatopancreas (HP) in the presence of a large (~70 kb) conjugative plasmid (pVA1) that contains the two toxin-subunit genes (pirA and pirB) (Yang et al., 2014; Han et al., 2015). Although the exact functions of the PirA and PirB subunits are still uncertain, the PirA subunit may play an initial stabilizing role, allowing PirB to bind with higher affinity to different glycan receptors located on the surface of the HP (Soto-Rodriguez et al., 2022).

Therefore, several DNA-based detection methods targeting toxin genes (pirA and pirB) and plasmid pVA1 have been developed, as recommended by the OIE (World Organization for Animal Health, 2017). However, those diagnostic techniques have faced difficulties due to the peculiar characteristics of the AHPND-causing pathogens as follows: (i) the emergence of mutant Vp^AHPND strains which have deletions of entire or partial pirA and/or pirB genes has been reported (Han et al., 2016); (ii) the emergence of atypical Vp^AHPND strains which contain full-length pirA and pirB but do not produce PirA^Vp toxins and fail to cause AHPND has been reported (Vicente et al., 2019). Therefore, there is an urgent need to develop a protein-based assay to detect the PirA and PirB subunits to provide more accurate and valid PCR-based assays for the detection of AHPND-causing pathogens.

Studies have reported that antibody-based immunoassays can be used to detect toxic subunits at the protein level. Monoclonal antibodies (mAbs) were previously developed by immunizing mice with concentrated culture supernatants prepared from Vp^AHPND isolates, and PirA or PirB subunits were detected by dot and western blotting analyses (Wangman et al., 2017; Wangman et al., 2018). In a recent study, an indirect enzyme-linked immunosorbent assay (iELISA) was developed with mAbs, which were generated by immunization with recombinant PirA or PirB; the indirect ELISA method detected each toxin subunit in protein lysates from several bacterial strains and AHPND-challenged shrimps (Mai et al., 2020). These antibody-based assays were successful in detecting either PirA or PirB subunit proteins; however, it has been shown that PirA and PirB interact with each other to form a heterodimeric PirAB^Vp and the complex causes significant mortality and morbidity in shrimp (Lee et al., 2015). These findings indicate that the virulence of AHPND-causing Vibrio spp. depends on heterodimeric PirAB^Vp. Thus, detecting the heterodimeric toxin complex is a challenging task for diagnosing AHPND in practice.

To the best of our knowledge, no antibody-based immunoassays are available for the detection of the heterodimeric protein complex PirAB^Vp. In this study, we developed a novel sandwich immunoassay to detect toxin complexes. Recombinant PirA and PirB subunits were produced and used as antigens for phage display panning using a naïve human single-chain variable fragment (scFv) antibody library. Selected scFvs were converted to IgG forms, and anti-PirA and anti-PirB antibodies were utilized to develop a sandwich ELISA-based immunoassay for capture and detection, respectively. The immunoassay exhibited specific detection activity for heterodimeric PirAB^Vp, but not for PirA or PirB, and successfully detected heterodimeric PirAB^Vp in protein lysates extracted from several Vibrio spp. and Vp^AHPND-infected shrimp.

2 Materials and methods

2.1 Cloning and production of recombinant PirA and PirB subunits

Full-length protein-coding sequences of pirA and pirB were obtained by PCR using the Korean Vp^AHPND strain 19–021-D1 as a DNA template (Han et al., 2020). PCR amplifications were conducted using the primer sets Vp_PirA_1F/Vp_PirA_336R and PirB_1F/PirB_1317R, according to Lee et al. (2015). After the exact nucleotide sequences were confirmed, the codons of pirA and pirB genes were optimized for expression in the Escherichia coli BL21 (DE3) expression system using the algorithm “OptimumGene” (GenScript, USA), synthesized de novo, and cloned into the pET-21b (+) expression vector (Novagen, USA). Recombinant protein expression and purification services were provided by Customs Protein Purification Service (AbClon Inc., Seoul, Korea).

2.2 Selection and screening of anti-PirA and anti-PirB scFv antibodies

For antibody selection, a naïve phage antibody display library from pooled healthy donors was constructed using the pADL-22c phagemid vector (Antibody Design Labs, USA), as described previously (Lee et al., 2018). Recombinant PirA or PirB was prepared at 10 μg/mL in phosphate buffer saline (PBS, pH 7.4) and coated on immunotube (Thermo Scientific, USA) by overnight incubation at 4°C. Immunotubes coated with each toxin subunit were incubated with phages that had been pre-blocked in PBS/2% nonfat milk (MPBS) for 1 h. Immunotubes were then washed with PBS/0.05% Tween20 (PBST) and the bound phages were eluted and propagated as previously described (Lee et al., 2018). To assess the binding activity of polyclonal scFvs, 10 μl of outputs from each round was inoculated in 5 ml 2xYT/Ampicillin (100 μg/ml) supplemented with 2% glucose and cultured until the OD600 reaches 0.4–0.5. Bacteria was pelleted by centrifugation and resuspended using 5 ml 2xYT/Ampicillin supplemented with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG), then cultured at 30°C with shaking for 16 hr. The culture was centrifuged at 4,000 rpm for 15 min and lysed with 125 μl of periplasmic extraction buffer (PEB; 30 mM Tris-HCl/20% sucrose/1 mM EDTA, pH 8.0), followed by additional periplasmic extraction using 250 μl of 5 mM MgSO₄. Two supernatant lysates extracted using PEB and MgSO₄ were mixed and used for ELISA with HRP-conjugated anti-HA antibody (Roche, Switzerland) as a secondary antibody. To screen individual scFv antibodies, single colonies picked from the output of 2^nd and 3^rd rounds were inoculated into 150 μl of 2xYT/Ampicillin/2% glucose in a 96-well U-bottom culture plate and grown for 3 hr. Bacteria was pelleted by centrifugation, resuspended with 150 μl of 2xYT/Ampicillin with 1 mM, and grown for overnight at 30°C with shaking and proper aeration. After IPTG induction, the plate was centrifuged and scFvs extracted from each pellet using 100 μl of PEB were subjected to screening by ELISA using HRP-conjugated anti-HA antibody. Clones that demonstrated at least five-fold binding over background were submitted for Sanger sequencing analysis (Bioneer, Republic of Korea) using a pelB-Forward primer (5′-AATACCTATTGCCTACGGCTG-3′) for identifying the variable heavy (VH) or variable light (VL) region sequence.

2.3 scFv purification

Selected monoclonal scFvs were cultured in 100 ml 2xYT/Ampicillin/2% glucose until the absorbance at 600 nm was between 0.5–0.8. bacterial culture was centrifuged, and the pellet was resuspended in 100 ml of 2xYT/Ampicillin/1 mM IPTG and grown for overnight at 30°C with vigorous shaking. The periplasmic fraction was obtained using 2.5 ml of PEB, followed by resuspension of the remaining pellet in 5 ml of 5 mM MgSO₄. The two fractions from the extraction procedure were mixed and subjected to scFv purification using an NGC Quest 100 column (Bio-Rad Laboratories, USA) equipped with a HisTrap™ HP column (Cytiva, USA).

2.4 Reformatting, transient expression, and purification of IgGs

Using pcDNA3.4 as a backbone vector, heavy chain (HC) and light chain (LC) expression vectors were constructed by inserting HC (CH1–CH2–CH3) and LC (CL) constant regions, respectively. Restriction enzyme sites were inserted in each expression vector for cassette cloning of VH and VL. The VH or VL genes in each scFv were amplified by PCR with the appropriate restriction enzyme sites and cloned into the HC or LC expression vector. HC and LC expression vectors for producing each IgG were co-transfected into Expi293F cells (1 × 10⁸) in 50 ml Expi293F expression medium. Cells were cultured for 5 days post-transfection, and supernatants were harvested by centrifugation at 4,000 rpm for 40 min and filtered using a 0.22-μm bottle-top vacuum filter. The produced IgGs were purified using an NGC Quest 100 equipped with a MabSelect™ PrismA column (Cytiva, USA), dialyzed in PBS (pH 7.4), and analyzed by SDS-PAGE.

2.5 Indirect ELISA

To validate the binding of purified scFvs or IgGs by ELISA, recombinant PirA or PirB diluted in 100 μl of PBS (200 ng/well) was coated on a 96-well plate for overnight. The plates were washed with PBST and blocked with MPBS at room temperature for 2 h. Various concentrations of the antibodies were allowed to bind to each antigen for 1 h, followed by three washes with PBST. Bound scFvs and IgGs were detected using horseradish peroxidase (HRP)-conjugated anti-HA antibody and goat anti-human IgG (Fc-specific) (Thermo Scientific), respectively, by reacting for 30 min at room temperature. After washing with PBST three times, 100 μl of OptEIA tetramethylbenzidine (TMB) (BD Biosciences, USA) was added to each well and incubated for 5 min, followed by terminating the reaction using 100 μl of 2N H₂SO₄. Absorbance was measured at 450 nm using a SpectraMax ABS Plus plate reader (Molecular Devices, USA).

2.6 BLI Analysis

An amine-reactive 2^nd generation (AR2G) biosensor (ForteBio, USA) was used to immobilize recombinant PirA or PirB as described previously (Jang et al., 2022). The biosensors were quenched with 1 M ethanolamine–HCl (pH 8.5) for 300 s, and a zero baseline was obtained with PBS for 120 s. Kinetic analysis was performed using the Octet K2 system (ForteBio) by determining the association (K_on) and dissociation (K_off) of anti-PirA and anti-PirB antibodies for 600 s. The K_D values of each antibody were calculated based on K_on and K_off using data analysis software (HT 12.0; ForteBio).

2.7 Detection of PirA or PirB from AHPND-associated Vibrio strains

Four AHPND-pathogenic strains of V. parahaemolyticus (19-022-A1, 19-021-01, CH49, 13-028-A3), one AHPND-mutant strains of V. parahaemolyticus (13-511-A2), and one V. harveyi strain (LB4) originated from different geographical origins were used in this study (8, 10, 20) (Table 1). Vibrio spp. were isolated from the stomachs of moribund shrimp and bacterial identifications were carried out by PCR assay targeting species specific genes (toxR and TopA genes), and PCR assay targeting pirA- and pirB-like genes (Han et al., 2015). These strains were overnight cultured in TSB + (Tryptic soy broth plus 2% NaCl) at 28–29°C with gentle (100 rpm) shaking and cell pellets were obtained by centrifugation at 8,000 x g for 15 min. To extract total proteins, the pellets were lysed using 4 ml B-PER^® bacterial protein extraction reagent (Thermo Fisher Scientific, USA) supplemented with 27.3 U of DNase (Zymo Research, USA) and EDTA-free protease inhibitor (GenDEPOT, USA) in ice for 20 min, and protein lysates were harvested by centrifugation at 15,000 × g for 5 min. Lysates (250 μl/well) were coated on a 96-well plate and PirA or PirB was detected by ELISA, as described in the Indirect ELISA section.

TABLE 1

Table 1 Measurement of antibody affinity based on biolayer interferometry.

2.8 Sandwich ELISA

100 μl of anti-PirA IgG (15 μg/ml) was coated on a 96-well plate for overnight at 4°C. After blocking with MPBS for 1 h, a mixture of recombinant PirA and PirB at 1:1 ratio (each 100 nM) was serially diluted in PBS and added to each well, followed by incubation for 2 h at room temperature. Anti-PirB IgG was biotinylated for detection as previously described (Jang et al., 2022). After washing the plate with PBST three times, 100 μl of biotinylated anti-PirB IgG (15 μg/ml) was added to each well and incubated for 1 hr at room temperature. The plate was washed thrice with PBST and incubated with HRP-conjugated streptavidin (Sigma Aldrich, USA) for 1 h at room temperature. After washing thrice with PBST, a colorimetric reaction was conducted as described in the Indirect ELISA section.

2.9 Detection of PirAB^Vp in hepatopancreas lysates of Vp^AHPND-infected shrimps

17 and 15 P. vannamei (mean weight = 1 g) were stocked in 17-L tanks as infected and non-infected groups. Vp^AHPND-causing strain 13-028/A3 was grown to 1 × 10⁹ CFU/ml, mixed with shrimp feed at a 1:1 ratio, and fed to the infected group. In the non-infected group, the shrimp were fed normally without bacteria. On day 7, all shrimp were dissected under aseptic conditions to prepare HP, as described previously (Han et al., 2020). The HP samples were kept at –80°C until further use. The animal use and experimental protocols were reviewed and approved by the Animal Research Ethics Committee of the Korea Research Institute of Bioscience and Biotechnology (IACUC approval no. KRIBB-AEC-21119; April 2021). To extract total protein from shrimp HP, the stored samples were thawed at 4°C and lysed using 300 μl of 1X RIPA lysis buffer (Merck Millipore, USA) supplemented with phosphatase and protease inhibitors (GenDEPOT) in ice for 40 min. Total protein lysates were obtained by centrifugation at 10,000 rpm for 40 min, and the protein concentration was measured using the BCA assay. Protein lysates (20 μg/well) were coated and reacted to detect PirA or PirB as described in the Indirect ELISA section and PirAB^Vp as described in the Sandwich ELISA section.

2.10 Theoretical limit of detection

LoD was defined as the lowest PirAB^Vp concentration of the detected colorimetric signal, which was greater than non-specific binding. Various concentrations of recombinant PirAB^Vp or protein lysates were analyzed by sandwich ELISA, and linear regression analyses were conducted using GraphPad Prism 8.0 (GraphPad Software, USA). The standard deviation of the response (σ) and the slope (S) of the linear regression curve were used to calculate LoD using the equation described below.

$$\text{LoD }=3.3\times \text{ σ}/\text{S}$$

2.11 Antibody modeling and docking analysis

The structural model of the antibody variable fragment (Fv) was generated using RosettaAntibody on the ROSIE server using the VH and VL chain sequences of each antibody, as described previously (Weitzner et al., 2017). Docking models between 1A8 Fv and PirA (PDB:3X0T) or between 3A5 and PirB (PDB:3X0U) were generated using ZDOCK (Pierce et al., 2014). Fv-binding amino acid residues on either PirA or PirB (≤ 4 Å) were analyzed from each docking model and visualized using the PyMOL Molecular Graphics System.

2.12 Statistical analysis

The data are shown as the means ± the standard deviation (SD) of the means, and statistical data analyses were performed via unpaired two-tailed Student’s t-tests using GraphPad Prism 8.0 (GraphPad Software, USA). P values of < 0.05 indicated statistically significant differences.

3 Results

3.1 Screening of scFv antibodies binding to PirA or PirB subunit

To prepare recombinant PirA and PirB, plasmids encoding each toxin gene with a 6xHis tag were separately transformed into E. coli BL21 (DE3) cells, and the recombinant toxins were expressed by IPTG induction. Pure recombinant PirA or PirB proteins corresponding to approximately 13 or 50 kDa were purified using a Ni-NTA column and analyzed by SDS-PAGE and Coomassie Brilliant Blue staining (Figure 1A). To select scFv antibodies against the toxins, we utilized a naïve human scFv antibody library and performed three iterative panning rounds using purified recombinant PirA or PirB. After the 3^rd round of panning, polyclonal scFv antibodies from each round were expressed by IPTG induction and their binding activity to each antigen was evaluated using ELISA. As shown in Figure 1B, the binding of polyclonal anti-PirA scFv antibodies expressed from the 3^rd round output was highly increased compared to that of polyclonal scFv antibodies derived from the previous rounds or the original library, indicating that the panning rounds were successfully performed and target-binding clones were amplified in the 3^rd round output. In addition, polyclonal scFvs expressed in the 2^nd and 3^rd rounds showed similarly increased binding to PirB (Figure 1C). By screening individual scFv clones from the 3^rd round of output, we found 22 and 30 positive binders to PirA and PirB, respectively (Figure 1D), and performed DNA sequencing analysis of the clones to determine the sequences of the variable heavy (VH) and variable light (VL) chains.

FIGURE 1

Figure 1 Phage display panning and screening of scFv antibodies binding to recombinant PirA or PirB. (A) Production of recombinant PirA and PirB proteins. E. coli BL21-DE3 transformed with PirA- or PirB-expressing vector was cultured and induced by IPTG for protein expression. Each toxin subunit expressed with a 6xHis tag was purified by Ni-NTA affinity chromatography, and the molecular weight of PirA (~13 kDa) and PirB (~50 kDa) was validated using SDS-PAGE. (B, C) Phage display panning using PirA and PirB. A naïve human scFv phage library was allowed to bind to each toxin subunit coated on an immunotube, then panning rounds were performed by following washing and elution procedure. After 3rd round of panning, polyclonal scFv antibodies from each round were expressed by IPTG induction and used to evaluate their binding activity to (B) PirA, or (C) PirB, by ELISA. Data represent the mean ± SD of duplicate tests (B). ***P < 0.001. (D) Screening of anti-PirA and anti-PirB scFv antibodies. Individual TG1 colonies from the 2^nd and 3^rd round outputs were inoculated and cultured to express scFv antibodies in periplasmic region by IPTG induction. Screening of scFvs obtained from periplasmic fraction was conducted on an ELISA plate coated with each antigen. The positive ratio (OD 450 nm > 0.5) of clones tested was approximately 15% (22/144) for PirA and 31% (30/96) for PirB, respectively.

3.2 Evaluation of selected anti-PirA and anti-PirB antibodies

Based on sequencing analysis of scFv antibodies, four and two scFv antibodies, which were determined at different frequencies and sequences of the third complementarity-determining region in the variable heavy chain (CDR-H3), were identified against PirA and PirB, respectively (Figure 2A). Each scFv expressed in the periplasmic fraction after IPTG induction was purified using a Ni-NTA column, and its binding activity to each antigen was evaluated by ELISA. As shown in Figure 2B, anti-PirA 1A8 scFv exhibited the highest binding activity among the tested clones, with an apparent K_D of 7.6 nM. In addition, 3A5 scFv bound to PirB to a greater extent than did 3G7 scFv, and the apparent K_D of 3A5 was estimated at 26.3 nM (Figure 2C). Thus, we selected 1A8 and 3A5 scFv antibodies for further use.

FIGURE 2

Figure 2 Selection of a pair of anti-PirA and anti-PirB scFvs for developing an immunoassay. (A) CDR-H3 sequence and frequency of anti-PirA and anti-PirB scFv clones. Four anti-PirA scFvs (1B3, 1A8, 1C10, and 2C4) and two anti-PirB scFvs (3G7 and 3A5) were identified after sequencing analysis, and the frequency of each scFv antibody was determined. (B, C) Validation of binding activity of scFv clones. Each scFv with 6xHis tag was purified by Ni-NTA column and was titrated to assess the binding activity to (B) PirA, or (C) PirB by ELISA.

To improve the antigen-binding activity by the avidity effect and apply it to a sandwich ELISA setting, each scFv was converted into an IgG form. As depicted in Figure 3A, the VH and VL regions of each scFv gene were amplified by PCR to insert two restriction sites on the N- and C-termini of the gene and cloned into the heavy chain (HC) or light chain (LC) expression vector, which was modified by inserting constant regions of the heavy or light chain. After co-transfection of HC and LC vectors constructed for each antibody into Expi293F cells, 1A8 and 3A5 were expressed, purified, and analyzed using SDS-PAGE. IgGs were verified under non-reducing and reducing conditions and showed intact and pure production (Figure 3B). The binding activity of the IgGs was validated by titrating each antigen. As shown in Figure 3C, IgGs 1A8 and 3A5 specifically bound to PirA and PirB in a dose-dependent manner, but not vice versa. We further measured the affinities (K_D) of 1A8 and 3A5 IgG using biolayer interferometry (BLI). Recombinant PirA or PirB were immobilized on a biosensor tip, and three concentrations of each antibody were allowed to bind to each tip (Figure 3D). The association (K_on) and dissociation (K_off) values were calculated, and the K_D values of 1A8 and 3A5 were 1.68 nM and 0.29 nM, respectively (Table 1). Thus, we demonstrated that anti-PirA 1A8 and anti-PirB 3A5 are purely produced in IgG form and have a prominent binding affinity for each antigen, with promising applications in immunoassays.

FIGURE 3

Figure 3 Production and validation of 1A8 and 3A5 IgGs. (A) Schematic illustration of transition of scFv to IgG form. VH and VL chains of the selected scFvs were separately cloned into HC and LC expression vectors containing constant regions, respectively; HC and LC expression vectors constructed for each antibody were transiently co-transfected into Expi293F cells. (B) Purification of 1A8 and 3A5 IgGs. IgGs expressed in supernatants were purified by protein A column and analyzed by SDS-PAGE to evaluate their intact sizes, approximately 150 kDa in non-reducing condition, and 50 kDa (HC) and 25 kDa (LC) in reduction condition. (C) Validation of the binding activity of 1A8 and 3A5 IgGs. Dose-dependent binding of 1A8 and 3A5 IgGs was assessed using PirA (left) and PirB (right) by ELISA, respectively. (D) Affinity measurement by biolayer interferometry. PirA or PirB (200 nM) was immobilized on an AR2G biosensor and allowed to bind to 1A8 (left) or 3A5 (right) IgG (25, 75, and 100 nM). Kinetic rates and equilibrium binding constants were analyzed using global fitting analysis of the binding curves.

3.3 Development of a sandwich immunoassay for the detection of heterodimeric PirAB complex

Since the antibodies should detect the toxin subunits naturally expressed in V. parahaemolyticus strains, we investigated whether 1A8 and 3A5 exhibited binding activity towards PirA and PirB, respectively, in the protein lysates of Vibrio spp. To date, various AHPND-associated Vibro strains have been isolated from the environment or infected shrimp. We chose several V. parahaemolyticus strains in which the presence of pirA and/or pirB and the virulence causing AHPND have been confirmed (Supplementary Table 1). Total protein lysates were prepared from the strains, and the production of PirA or PirB was validated by ELISA using 1A8 or 3A5. As shown in Figure 4A, PirA was clearly detected in the lysates extracted from three strains (19-022-A1, 19-021-D1, and CH49), whereas pirA and pirB were detected by PCR, but not from V. harveyi (LB4) or other V. parahaemolyticus strains diagnosed as mutant Vp^AHPND strains. 3A5 exhibited a similar pattern of PirB detection in the lysates (Figure 4B). Next, we designed a sandwich ELISA-based immunoassay using 1A8 IgG as a capture antibody and biotinylated 3A5 IgG as a detector for PirAB^Vp, as depicted in Figure 4C. We found that the immunoassay detected not PirA or PirB but PirAB^Vp, indicating that it could only detect the heterodimeric toxin complex and not each monomeric toxin subunit (Figure 4C). We also investigated the sensitivity of the immunoassay for detecting PirAB^Vp. According to linear regression analysis after titration of recombinant heterodimeric PirAB^Vp, the LoD of the immunoassay was found to be 201.7 ng/ml (Figure 4D). We further evaluated the detectability of the immunoassay for PirAB^Vp directly produced from Vp^AHPND strains by titrating protein lysates extracted from 19-021-D1 and CH49. The linear regression analyzes showed significant coefficient of determination (R²) values in a concentration-dependent manner, and the detection sensitivity of the immunoassay was estimated as low as about 31.3 μg/ml of protein lysates extracted from both strains (Figure 4E). These data imply that the PirA and PirB toxin subunits are specifically produced by Vp^AHPND strains and natively form heterodimeric PirAB^Vp, which is detectable by the developed immunoassay.

FIGURE 4

Figure 4 Development of a sandwich immunoassay for the detection of heterodimeric PirAB^Vp complex. (A, B) Detection of PirA and PirB in protein lysates obtained from several Vibrio spp. Total proteins were extracted from culture pellets of six Vibrio strains and coated on a 96-well plate for ELISA. Anti-PirA 1A8 and anti-PirB 3A5 IgGs were used to detect (A) PirA, and (B) PirB, separately in the lysates. Data represent the mean ± SD of duplicate tests (A, B). *P < 0.05, **P < 0.01, and ***P < 0.001. (C) Detection of heterodimeric PirAB^Vp by a sandwich immunoassay. 1A8 and biotinylated 3A5 were used as a capture and a detector, respectively (left), and detectability of the sandwich immunoassay was assessed using PirA, PirB, and PirAB^Vp (right). Values represent the mean ± SD for a duplicate. *P < 0.05, **P < 0.01. (D) Sensitivity of the immunoassay in detecting PirAB. PirA and PirB were mixed together in 1:1 molar ratio, serially diluted, and detected by the immunoassay. Linear regression analysis was performed. (E) Detection of PirAB^Vp naturally produced by Vibrio spp. Protein lysates extracted from culture pellets of 19-021-D1 (left) and CH49 (right) were subject to the immunoassay.

3.4 Detection of PirAB^Vp in shrimps infected with Vp^AHPND 13-028/A3

We investigated whether the immunoassay could detect PirAB^Vp in Vp^AHPND-infected shrimp. V. parahaemolyticus 13-028/A3, a virulent Vp^AHPND strain that has been mostly used in AHPND models, was selected to infect shrimp. We first validated the expression of PirA and PirB by indirect ELISA using 1A8 and 3A5 IgGs, respectively, and determined that 13-028/A3 produced high levels of each toxin subunit (Figure 5A). In addition to detecting PirA and PirB, we analyzed the sensitivity in detecting PirAB^Vp from the lysate. As shown in Figure 5B, the LoD of the sandwich immunoassay was estimated to be 13.5 μg/ml by the linear regression, showing that the immunoassay could be sensitive enough to detect the heterodimeric PirAB^Vp from to 13-028/A3. Given the sensitivity of the immunoassay developed in this study, we examined its applicability to the diagnosis of AHPND in shrimp. We performed a challenge study by infecting shrimp with the strain and compared the PirAB^Vp detection activity between Vp^AHPND-infected and non-infected groups. We found that the immunoassay performed using hepatopancreas lysates exhibited significant detectability of PirAB^Vp in the infected group compared to the non-infected group, suggesting that 13-028/A3 highly produced heterodimeric PirAB^Vp after infecting shrimp (Figure 5C). Furthermore, we assessed the amount of PirAB^Vp in live, moribund, and dead shrimp in the infected group. The results shown in Figure 5D indicate a significantly higher level of PirAB^Vp in HP lysates from dead shrimp than in those from live shrimp post-infection. Thus, we inferred that PirAB^Vp is detectable in Vp^AHPND-infected shrimp using the developed immunoassay and plays a pivotal role in increasing shrimp mortality by inducing fatal virulence.

FIGURE 5

Figure 5 Detection of PirAB^Vp in shrimps infected by V. parahaemolyticus 13-028/A3. (A) Detection of the toxin subunits from the culture lysate of 13-028/A3 strain. Total proteins were extracted from the pellet of cultured bacteria and coated on an ELISA plate; the expression of PirA and PirB was determined by 1A8 and 3A5, respectively. CH49 and LB4 strains were used as a PirA/PirB-positive and -negative controls, respectively. Values represent the mean ± SD for a duplicate. *P < 0.05, **P < 0.01. (B) Assessment of detection sensitivity of PirAB^Vp in 13-028/A3 lysate. The culture lysate was serially diluted and the detection sensitivity (LoD) was determined estimated as 13.5 ug/mL. (C, D) Detection of PirAB^Vp complex in shrimps infected with 18-028/A3 strain. At day 7 after infection, hepatopancreas from non-infected or infected shrimps were dissected and lysed to prepare protein lysates. PirAB^Vp complex in the lysates was significantly detected in the infected group, but not in the non-infected group. The immunoassay showed a significant detection of PirAB^Vp complex from dead shrimps compared to live shrimps in the infected group.

3.5 Structural modeling reveals the binding mode of 1A8 and 3A5 to the toxin subunits

Given the amino acid sequences of the antibodies identified in this study and the previously analyzed structural information of PirA and PirB (Lee et al., 2015), we investigated the binding site of each antibody on PirA or PirB by antibody-antigen docking analysis. Structural models of variable fragments (Fv) were generated with VH and VL sequences of the 1A8 or 3A5 antibodies using RosettaAntibody, and CDR loops in the VH (CDR-H1, -H2, and -H3) or VL (CDR-L1, -L2, and -L3) were observed in each model (Figures 6A, B). Using the structures of PirA and PirB, antibody-antigen docking analysis was conducted with each Fv model to determine which amino acid residues on PirA and PirB interact with the 1A8 Fv and 3A5 Fv, respectively. As shown in Figure 6C, the CDR-H3 and CDR-L3 loops of 1A8 seemed to interact mainly with the helical region located between the fourth and fifth beta-sheets on PirA. In the case of PirB, the docking data showed that the CDR-Hs of 3A5 were mainly involved in binding to the hydrophilic region exposed as a loop on PirB (Figure 6D). Given the docking analysis data, we found that 10 amino acids on PirA and 14 amino acids on PirB might be the epitope of 1A8 and 3A5, respectively, since they were the closest residues to the CDR regions (≤4 Å). Next, we visualized two regions on each antigen: (i) the antibody-binding site expected as an epitope of each antibody, and (ii) the interaction site between PirA and PirB for heterodimerization. The results shown in Figures 6E, F indicate that the epitopes of 1A8 and 3A5 were located away from the interaction site between PirA and PirB, respectively. Taken together, these results suggest that 1A8 and 3A5 IgGs do not interfere with the interaction site on each antigen that forms the heterodimeric PirAB^Vp; therefore, the immunoassay was able to successfully detect recombinantly expressed and naturally produced PirAB^Vp.

FIGURE 6

Figure 6 Structural docking to analyze the binding mode of 1A8 and 3A5 to each toxin subunit. (A, B) Antibody variable fragment (Fv) modeling. VH and VL amino acid sequences of 1A8 (A) or 3A5 (B) were used for Fv modeling using RosettaAntibody. CDR regions in VH (red) and VL (green) are highlighted in 1A8 or 3A5 Fv model. (C, D) Antibody-antigen docking analysis. The 1A8 Fv and PirA (PDB:3X0T) (C) or the 3A5 Fv and PirB (PDB:3X0U) (D) were utilized for structural docking analysis using ZDOCK. (E) Epitope analysis in 1A8-PirA docking. Potential binding sites (within 4 Å) by CDRs in VH and VL of 1A8 were analyzed based on the docking results and highlighted on PirA. Red: residues involved in 1A8 binding (W57, G58, A59, P60, F61, M62, A63, G64, G65, and K67); Yellow: residues involved in the interaction with PirB (Y11, S12, H13, D14, W15, T16, V17, V26, D27, S28, K29, H30, G104, F105, C106, T107, I108, Y109, Y110) (Lin et al., 2017). (F) Epitope analysis in 3A5-PirB docking. Potential binding sites (within 4 Å) by CDRs in VH and VL of 3A5 were analyzed and highlighted on PirB. Red: residues involved in 3A5 binding (I317, E318, I319, H320, Y321, N368, G369, P370, E371, Q413, E414, G415, S416, and D417); Yellow: residues involved in the interaction with PirA (Y35, A36, F37, K38, A39, M40, V41, S42, F43, G44, L45, S46, N47, M247, L248, I249, W250, Q251, K252, I253, K254, E255, L256, D260, V261, F262, V263, H264, S265, N266, L267, I268, S269, Y270, P298, N299, M300, F301, G302, E303, R304, R305, P431, D432, E433, and F434) (Lin et al., 2017).

4 Discussion

It has been shown that PirA and PirB interact to form a heterodimeric PirAB^Vp and the complex causes significant mortality and morbidity in shrimp (Lee et al., 2015; Zhang et al., 2021). These results indicate that neither PirA nor PirB, but PirAB^Vp in binary form, is indispensable for causing significant toxicity after Vp^AHPND infection in hosts. Thus, it is important to detect the binary toxin to accurately determine Vp^AHPND virulence. In the present study, we demonstrated that 1A8 and 3A5 antibodies specific to PirA and PirB, respectively, can serve as a great option for the development of a novel sandwich immunoassay to detect the PirAB^Vp binary toxin. The immunoassay exhibited significant detectability of the binary toxin from not only the culture lysates of Vp^AHPND strains, but also the protein lysates of the hepatopancreas dissected from shrimp infected with the well-known Vp^AHPND strain 13-028-A3. Therefore, the novel immunoassay established in this study could specifically detect the PirAB^Vp binary toxin and was further developed for periodic surveillance of newly emerging Vibrio spp.

Although several recent studies have suggested that V. parahaemolyticus is the major bacterial species that can cause AHPND in shrimp because of the acquisition of a 70-kb plasmid (pVA1) encoding the binary toxin PirAB^Vp, evidence of the transferability of pVA1-type plasmids between other Vibrio species have been found (Han et al., 2015; Lee et al., 2015; Dong et al., 2019). Moreover, the PirAB^Vp binary toxin has been identified in other Vibrio species belonging to the Harveyi clade, including V. campbellii, V. harveyi, and V. owensii (Kondo et al., 2015; Ahn et al., 2017; Dong et al., 2017; Xiao et al., 2017), and even in the Gram-positive bacterial species Micrococcus luteus (Durán-Avelar et al., 2018). PirAB^Vp binary toxins have also been reported in V. punensis, which belongs to the Orientalis clade (Restrepo et al., 2018). Whether it is still controversial that the global spread of Vp^AHPND in the shrimp industry is mainly responsible for the transferability of the pVA1-type plasmids (Fu et al., 2020), the presence of pirA and pirB genes in various bacterial species definitely poses a potential risk for the spread of emerging diseases and the emergence of mutant or atypical Vp^AHPND strains also acts as an internal hurdle to accurate diagnosis of the AHPND-causing Vibrios. Therefore, we evaluated the sensitivity of the newly developed antibody to several virulent Vp^AHPND isolates, avirulent pirA/pirB mutant Vp^AHPND isolates, and non-AHPND Vibrio isolates of the Harveyi clade. Although we were unable to confirm the potential binding activity of the newly developed antibody against the AHPND-associated virulent Vibrio spp. rather than V. parahaemolyticus because of the limited available strains, the antibody successfully detected the presence of the PirAB^Vp binary toxin in several virulent Vp^AHPND isolates from different geographical origins, including Korea, and also clearly separated the virulent and two types of avirulent pirA/pirB mutant Vp^AHPND isolates. Moreover, nonspecific binding of the newly developed antibody to the non-AHPND Vibrio isolate of the Harveyi clade was recorded. These results suggest that the newly developed antibody has strong potential for the accurate diagnosis of PirAB^Vp binary toxins of AHPND-causing Vp^AHPND in the global shrimp industry.

The established immunoassay successfully detected natural PirAB^Vp in the HP of Vp^AHPND-infected shrimp with an LoD of 13.5 ug/ml. Since the concentrations of hepatopancreatic proteins extracted from a single shrimp were on average 3 mg/ml (Supplementary Table 2), the immunoassay exhibited acceptable ability for the detection of PirAB^Vp in shrimp; therefore, achievement of higher sensitivity is not a significant consideration for the immunoassay. In addition, compared with the immunoassay and dot-blot methods previously reported (Wangman et al., 2017; Wangman et al., 2018), the immunoassay is the only approach available for determining the binary form of PirA and PirB subunits, indicating that the present technique is less time consuming and feasible to monitor the PirAB^Vp high-throughputly in a 96-well plate setting. We suggest that there is room for improvement in the immunoassay detectability. To this end, some modifications could be made by conjugating more biotin to the detector antibody 3A5 or achieving affinity maturation of capture or detector antibodies.

The crystal structures of PirA and PirB have been analyzed in previous studies, which have shown that the binary toxin is complexed heteromerically by interacting with specific amino acids in each subunit (Lee et al., 2015; Lin et al., 2017). However, the toxicity-inducing mode of action of PirAB^Vp remains unclear. The structural similarity between the PirAB^Vp binary toxin and Cry proteins has been reported, suggesting PirAB^Vp could work as a toxin that induces pore-forming activity in target cells (Lee et al., 2015; Lin et al., 2017). Cry proteins identified in Bacillus thuringiensis are insecticidal toxins composed of three active domains: i) pore-forming domain I, ii) receptor-binding domain II, and iii) sugar-binding domain III (Bravo et al., 2006; Bravo et al., 2013; Pardo-López et al., 2013; Adang et al., 2014). Based on previous investigations, the structures of the N- and C-termini of PirB are similar to those of Cry domains I and II, respectively, and the structure of PirA corresponds to Cry domain III, suggesting that the PirAB^Vp binary toxin possesses a Cry-like function (Lee et al., 2015). This suggests that the C-terminal PirB could bind to a receptor expressed on host cells and that the N-terminal PirB could exert pore-forming activity on the cellular membrane of the host, thereby causing cell death. Interestingly, the expected epitopes of the 3A5 antibody developed in this study were located in loops 14 (I317/E318/I319/H320/Y321), 18 (G369/P370/E371), and 22 (Q413/E414/G415/S416/D417) of PirB, corresponding to the potential receptor-binding region (Figure 6F). Thus, we infer that 3A5 could be utilized not only for the diagnosis of the PirAB^Vp binary toxin in Vp^AHPND-infected shrimp, but also for further therapeutic applications based on interference with the interaction of PirB and its counterpart receptor.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The animal study was approved by IACUC approval no. KRIBB-AEC-21119. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

M-YJ: Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft. JEH: Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing. DGL: Formal analysis, Investigation, Methodology, Visualization, Writing – original draft. Y-LC: Formal analysis, Investigation, Methodology, Writing – original draft. J-HJ: Investigation, Methodology, Writing – original draft. JL: Investigation, Resources, Writing – original draft. J-GP: Investigation, Resources, Writing – original draft. DHK: Investigation; Writing – review & editing. SYP: Investigation, Writing – original draft. WK: Investigation, Writing – original draft. KL: Funding acquisition, Methodology, Writing – review & editing. JHK: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing. N-KL: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Korea Research Institute of Bioscience and Biotechnology (KRIBB) Research Initiative Program (KGM5272323), by the National Research Foundation of Korea (NRF-2020R1I1A2068827) funded by the Ministry of Education, and by the Development of Technology for the Biomaterialization of Marine Fishery Byproducts of the Korea Institute of Marine Science and Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries (KIMST-20220128). This research was also supported by the National Institute of Fisheries Science (R2023004) of the Ministry of Oceans and Fisheries.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1294801/full#supplementary-material

References

Adang, M. J., Crickmore, N., Jurat-Fuentes, J. L. (2014). Diversity of Bacillus thuringiensis crystal toxins and mechanism of action. Adv. Insect Physiol. 47, 39–87. doi: 10.1016/B978-0-12-800197-4.00002-6

CrossRef Full Text | Google Scholar

Ahn, Y. S., Piamsomboon, P., Tang, K. F. J., Han, J. E., Kim, J. H. (2017). Complete Genome Sequence of Acute Hepatopancreatic Necrosis Disease-Causing Vibrio campbellii LA16-V1, Isolated from Penaeus vannamei Cultured in a Latin American Country. Genome Announc. 5, e01011–e01017. doi: 10.1128/genomeA.01011-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo, A., Gill, S. S., Soberón, M. (2006). Mode of action of Bacillus thuringiensis Cry and Cyt toxins and their potential for insect control. Toxicon. 49, 423–435. doi: 10.1016/j.toxicon.2006.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo, A., Gómez, I., Porta, H., García-Gómez, B. I., Rodriguez-Almazan, C., Pardo, L., et al. (2013). Evolution of Bacillus thuringiensis Cry toxins insecticidal activity. Microb. Biotechnol. 6, 17–26. doi: 10.1111/j.1751-7915.2012.00342.x

PubMed Abstract | CrossRef Full Text | Google Scholar

de la Peña, L. D., Cabillon, N. A., Catedral, D. D., Amar, E. C., Usero, R. C., Monotilla, W. D., et al. (2015). Acute hepatopancreatic necrosis disease (AHPND) outbreaks in Penaeus vannamei and P. monodon cultured in the Philippines. Dis. Aquat. Organ. 116, 251–254. doi: 10.3354/dao02919

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhar, A. K., Piamsomboon, P., Aranguren Caro, L. F., Kanrar, S., Adami, R., Jr., Juan, Y. S. (2019). First report of acute hepatopancreatic necrosis disease (AHPND) occurring in the USA. Dis. Aquat. Organ. 132, 241–247. doi: 10.3354/dao03330

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, X., Bi, D., Wang, H., Zou, P., Xie, G., Wan, X., et al. (2017). pirAB vp-bearing Vibrio parahaemolyticus and Vibrio campbellii pathogens isolated from the same AHPND-affected pond possess highly similar pathogenic plasmids. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.01859

CrossRef Full Text | Google Scholar

Dong, X., Song, J., Chen, J., Bi, D., Wang, W., Ren, Y., et al. (2019). Conjugative Transfer of the pVA1-Type Plasmid Carrying the pirAB(vp) Genes Results in the Formation of New AHPND-Causing Vibrio. Front. Cell. Infect. Microbiol. 9. doi: 10.3389/fcimb.2019.00195

CrossRef Full Text | Google Scholar

Durán-Avelar, M. J., Vázquez-Reyes, A., González-Mercado, A. L., Zambrano-Zaragoza, J. F., Ayón-Pérez, M. F., Agraz-Cibrián, J. M., et al. (2018). pirA-and pirB-like gene identification in Micrococcus luteus strains in Mexico. J. Fish Dis. 41, 1667–1673. doi: 10.1111/jfd.1287

PubMed Abstract | CrossRef Full Text | Google Scholar

Eshik, M. M. E., Abedin, M. M., Punom, N. J., Begum, M. K., Rahman, M. S. (2018). Molecular identification of AHPND positive Vibrio parahaemolyticus causing an outbreak in south-west shrimp farming regions of Bangladesh. J. Bangladesh Acad. Sci. 41, 127–135. doi: 10.3329/jbas.v41i2.35492

CrossRef Full Text | Google Scholar

Fu, S., Wei, D., Yang, Q., Xie, G., Pang, B., Wang, Y., et al. (2020). Horizontal plasmid transfer promotes the dissemination of asian acute hepatopancreatic necrosis disease and provides a novel mechanism for genetic exchange and environmental adaptation. mSystems. 5, e00799–e00719. doi: 10.1128/mSystems.00799-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Global Seafood Alliance (2012) Early mortality syndrome affects shrimp in Asia. Available at: https://www.globalseafood.org/advocate/early-mortality-syndrome-affects-shrimp-in-asia/.

Google Scholar

Han, J. E., Choi, S.-K., Han, S.-H., Lee, S. C., Jeon, H. J., Lee, C., et al. (2020). Genomic and histopathological characteristics of Vibrio parahaemolyticus isolated from an acute hepatopancreatic necrosis disease outbreak in Pacific white shrimp (Penaeus vannamei) cultured in Korea. Aquaculture. 524, 735284. doi: 10.1016/j.aquaculture.2020.735284

CrossRef Full Text | Google Scholar

Han, J., Tang, K., Aranguren, F., Piamsomboon, P. (2016). Characterization and pathogenicity of acute hepatopancreatic necrosis disease natural mutants, PirABVp (-) V. parahaemolyticus, and PirAB^Vp (+) V. campbellii strains. Aquaculture. 470, 84–90. doi: 10.1016/j.aquaculture.2016.12.022

CrossRef Full Text | Google Scholar

Han, J. E., Tang, K. F., Tran, L. H., Lightner, D. V. (2015). Photorhabdus insect-related (Pir) toxin-like genes in a plasmid of Vibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis disease (AHPND) of shrimp. Dis. Aquat. Organ. 113, 33–40. doi: 10.3354/dao02830

PubMed Abstract | CrossRef Full Text | Google Scholar

Jang, J. H., Kim, S., Kim, S. G., Lee, J., Lee, D. G., Jang, J., et al. (2022). A sensitive immunodetection assay using antibodies specific to staphylococcal enterotoxin B produced by baculovirus expression. Biosensors (Basel). 12, 787. doi: 10.3390/bios12100787

PubMed Abstract | CrossRef Full Text | Google Scholar

Joshi, J., Srisala, J., Truong, V. H., Chen, I.-T., Nuangsaeng, B., Suthienkul, O., et al. (2014). Variation in Vibrio parahaemolyticus isolates from a single Thai shrimp farm experiencing an outbreak of acute hepatopancreatic necrosis disease (AHPND). Aquaculture. 428, 297–302. doi: 10.1016/j.aquaculture.2014.03.030

CrossRef Full Text | Google Scholar

Kondo, H., Van, P. T., Dang, L. T., Hirono, I. (2015). Draft genome sequence of non-vibrio parahaemolyticus acute hepatopancreatic necrosis disease strain KC13.17.5, isolated from diseased shrimp in Vietnam. Genome Announc. 3, e00978–e00915. doi: 10.1128/genomeA.00978-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, N. K., Bidlingmaier, S., Su, Y., Liu, B. (2018). Modular construction of large non-immune human antibody phage-display libraries from variable heavy and light chain gene cassettes. Methods Mol. Biol. 1701, 61–82. doi: 10.1007/978-1-4939-7447-4_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C. T., Chen, I. T., Yang, Y. T., Ko, T. P., Huang, Y. T., Huang, J. Y., et al. (2015). The opportunistic marine pathogen Vibrio parahaemolyticus becomes virulent by acquiring a plasmid that expresses a deadly toxin. Proc. Natl. Acad. Sci. U. S. A. 112, 10798–10803. doi: 10.1073/pnas.1503129112

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, S.-J., Hsu, K.-C., Wang, H.-C. (2017). Structural insights into the cytotoxic mechanism of Vibrio parahaemolyticus PirA vp and PirB vp toxins. Mar. Drugs 15, 373. doi: 10.3390/md15120373

PubMed Abstract | CrossRef Full Text | Google Scholar

Mai, H. N., Cruz-Flores, R., Dhar, A. K. (2020). Development of an indirect Enzyme Linked Immunoassay (iELISA) using monoclonal antibodies against Photorhabdus insect related toxins, PirA(Vp) and PirB(Vp) released from Vibrio spp. J. Microbiol. Methods 176, 106002. doi: 10.1016/j.mimet.2020.106002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardo-López, L., Soberón, M., Bravo, A. (2013). Bacillus thuringiensis insecticidal three-domain Cry toxins: mode of action, insect resistance and consequences for crop protection. FEMS Microbiol. Rev. 37, 3–22. doi: 10.1111/j.1574-6976.2012.00341.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pierce, B. G., Wiehe, K., Hwang, H., Kim, B. H., Vreven, T., Weng, Z. (2014). ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinformatics. 30, 1771–1773. doi: 10.1093/bioinformatics/btu097

PubMed Abstract | CrossRef Full Text | Google Scholar

Restrepo, L., Bayot, B., Arciniegas, S., Bajana, L., Betancourt, I., Panchana, F., et al. (2018). PirVP genes causing AHPND identified in a new Vibrio species (Vibrio punensis) within the commensal Orientalis clade. Sci. Rep. 8, 13080. doi: 10.1038/s41598-018-30903-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinn, A. P., Pratoomyot, J., Bron, J. E., Paladini, G., Brooker, E. E., Brooker, A. J. (2015). Economic costs of protistan and metazoan parasites to global mariculture. Parasitology. 142, 196–270. doi: 10.1017/S0031182014001437

PubMed Abstract | CrossRef Full Text | Google Scholar

Sirikharin, R., Taengchaiyaphum, S., Sanguanrut, P., Chi, T. D., Mavichak, R., Proespraiwong, P., et al. (2015). Characterization and PCR detection of binary, pir-like toxins from vibrio parahaemolyticus isolates that cause acute hepatopancreatic necrosis disease (AHPND) in shrimp. PloS One 10, e0126987. doi: 10.1371/journal.pone.0126987

PubMed Abstract | CrossRef Full Text | Google Scholar

Soto-Rodriguez, S. A., Gomez-Gil, B., Lozano-Olvera, R., Betancourt-Lozano, M., Morales-Covarrubias, M. S. (2015). Field and experimental evidence of Vibrio parahaemolyticus as the causative agent of acute hepatopancreatic necrosis disease of cultured shrimp (Litopenaeus vannamei) in Northwestern Mexico. Appl. Environ. Microbiol. 81, 1689–1699. doi: 10.1128/AEM.03610-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Soto-Rodriguez, S. A., Lozano-Olvera, R., Ramos-Clamont Montfort, G., Zenteno, E., Sanchez-Salgado, J. L., Vibanco-Perez, N., et al. (2022). New insights into the mechanism of action of pirAB from vibrio parahaemolyticus. Toxins (Basel). 14, 243. doi: 10.3390/toxins14040243

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, L., Nunan, L., Redman, R. M., Mohney, L. L., Pantoja, C. R., Fitzsimmons, K., et al. (2013). Determination of the infectious nature of the agent of acute hepatopancreatic necrosis syndrome affecting penaeid shrimp. Dis. Aquat. Organ. 105, 45–55. doi: 10.3354/dao02621

PubMed Abstract | CrossRef Full Text | Google Scholar

Tun, K., Kanrar, S., Fitzsimmons, K., McLain, J., Aranguren, L., Oo, M., et al. (2017). Acute hepatopancreatic necrosis disease (AHPND) in black tiger shrimp (Penaeus monodon), Pacific white shrimp Penaeus vannamei), and fresh water shrimp (Macrobrachium rosenbergii) (Meeting abstract, Kuala Lumpur, Malaysia: Asia Pacific Aquaculture).

Google Scholar

Vicente, A., Taengphu, S., Hung, A. L., Mora, C. M., Dong, H. T., Senapin, S. (2019). Detection of Vibrio campbellii and V. parahaemolyticus carrying full-length PirAB^Vp but only V. campbellii produces PirVp toxins. Aquaculture. 519, 734708. doi: 10.1016/j.aquaculture.2019.734708

CrossRef Full Text | Google Scholar

Wangman, P., Chaivisuthangkura, P., Sritunyalucksana, K., Taengchaiyaphum, S., Senapin, S., Pengsuk, C., et al. (2017). Development of monoclonal antibodies specific to ToxA and ToxB of Vibrio parahaemolyticus that cause acute hepatopancreatic necrosis disease (AHPND). Aquaculture. 474, 75–81. doi: 10.1016/J.AQUACULTURE.2017.03.039

CrossRef Full Text | Google Scholar

Wangman, P., Longyant, S., Taengchaiyaphum, S., Senapin, S., Sithigorngul, P., Chaivisuthangkura, P. (2018). PirA & B toxins discovered in archived shrimp pathogenic Vibrio campbellii isolated long before EMS/AHPND outbreaks. Aquaculture. 497, 494–502. doi: 10.1016/j.aquaculture.2018.08.025

CrossRef Full Text | Google Scholar

Weitzner, B. D., Jeliazkov, J. R., Lyskov, S., Marze, N., Kuroda, D., Frick, R., et al. (2017). Modeling and docking of antibody structures with Rosetta. Nat. Protoc. 12, 401–416. doi: 10.1038/nprot.2016.180

PubMed Abstract | CrossRef Full Text | Google Scholar

World Organization for Animal Health (2017) Aquatic animal health code. Available at: https://rr-europe.woah.org/wp-content/uploads/2020/08/oie-aqua-code_2019_en.pdf.

Google Scholar

Xiao, J., Liu, L., Ke, Y., Li, X., Liu, Y., Pan, Y., et al. (2017). Shrimp AHPND-causing plasmids encoding the PirAB toxins as mediated by pirAB-Tn903 are prevalent in various Vibrio species. Sci. Rep. 7, 42177. doi: 10.1038/srep42177

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y. T., Chen, I. T., Lee, C. T., Chen, C. Y., Lin, S. S., Hor, L. I., et al. (2014). Draft genome sequences of four strains of vibrio parahaemolyticus, three of which cause early mortality syndrome/acute hepatopancreatic necrosis disease in shrimp in China and Thailand. Genome Announc. 2, e00816–e00814. doi: 10.1128/genomeA.00816-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Bader, T. J., Lange, M. D., Shoemaker, C. A., Beck, B. H. (2021). Toxicity of recombinant PirA and PirB derived from Vibrio parahaemolyticus in shrimp. Microb. Pathog. 155, 104886. doi: 10.1016/j.micpath.2021.104886

PubMed Abstract | CrossRef Full Text | Google Scholar