Skip to main content

Screening and identification of cyprinid herpesvirus 2 (CyHV-2) ORF55-interacting proteins by phage display



Cyprinid herpesvirus 2 (CyHV-2) is a pathogenic fish virus belonging to family Alloherpesviridae. The CyHV-2 gene encoding thymidine kinase (TK) is an important virulence-associated factor. Therefore, we aimed to investigate the biological function of open reading frame 55 (ORF55) in viral replication.


Purified CyHV-2 ORF55 protein was obtained by prokaryotic expression, and the interacting peptide was screened out using phage display. Host interacting proteins were then predicted and validated.


ORF55 was efficiently expressed in the prokaryotic expression system. Protein and peptide interaction prediction and dot-blot overlay assay confirmed that peptides identified by phage display could interact with the ORF55 protein. Comparing the peptides to the National Center for Biotechnology Information database revealed four potential interacting proteins. Reverse transcription quantitative PCR results demonstrated high expression of an actin-binding Rho-activating protein in the latter stages of virus-infected cells, and molecular docking, cell transfection and coimmunoprecipitation experiments confirmed that it interacted with the ORF55 protein.


During viral infection, the ORF55 protein exerts its biological function through interactions with host proteins. The specific mechanisms remain to be further explored.


Herpesviral haematopoietic necrosis (HVHN), a highly lethal disease of crucian carp (Carassius carassius), is caused by cyprinid herpesvirus 2 (CyHV-2). This disease was first detected in goldfish (Carassius auratus) cultured in Aichi and Nara prefectures in Japan in 1992, with a mortality rate of 100% [1]. Since 2009, CyHV-2 has been detected in Gibel Carp (Carassius auratus gibelio) cultured in some areas of Jiangsu Province in China, and the virus rapidly became endemic throughout the main carp farming areas in China [2]. CyHV-2 has since dramatically spread across the world [3,4,5,6]. CyHV-2 is a member of the Alloherpesviridae family (genus Cyprinivirus), which also comprises cyprinid herpesvirus 3 (CyHV-3) and carp pox virus (CyHV-1) [7]. CyHV-2 is a DNA virus with a vesicular membrane and a large linear double-stranded DNA genome. The virions are oval, with a diameter of approximately 170 − 200 nm [8].

Developing the molecular epidemiological characterisation of CyHV-2 is an important basis for research on the infection mechanism, and for prevention and control of the virus. Some genes of CyHV-2 have been used for research on virus detection, preparation of antibody, and genotype distinction and diversity analysis [9,10,11]. Additionally, the CyHV-2 open reading frame 55 (ORF55) gene was predicted to encode thymidine kinase (TK), and enzyme participating in nucleic acid metabolism [12]. In previous studies, the function of TK was investigated in herpes simplex virus 1 (HSV-1) [13,14,15], duck plague virus (DPV) [16] and Epstein-Barr (EB) virus [17]. The TK gene has since been studied in Pseudorabies virus (PRV) [18]. The TK gene was not essential for the growth of PRV, but it was essential for viral proliferation, and it was also the most important gene determining PRV critical virulence. Additionally, deletion of the TK gene could significantly decrease the virulence of PRV, while also greatly reducing the replication, transmission and lethal ability of the virus in nerve tissues, but immunogenicity was not affected [19].

Thus, TK plays an important role during viral infection. Although the experimental evidence mentioned above linked the TK gene with virulence, its overall specific mechanism of action remains elusive. From previous experimental results, we deduced that ORF55 also has a virulence effect, and is an important virulence gene in CyHV-2. Therefore, further exploration of the ORF55 gene is important.

Phage display technology (PDT) is a powerful approach first reported in 1985 [20]. PDT is a well-established method for selecting peptides and proteins with novel binding functions [21]. A foreign protein is first expressed in fusion with one of the phage capsid proteins by cloning the foreign DNA sequence in frame with the gene encoding the chosen capsid protein [22]. PDT allows a direct link and identification between a large number of random peptides and DNA-coding sequences. It is used extensively in medicine, biology and chemistry [23, 24]. However, how the CyHV-2 ORF55 protein performs its biological functions in goldfish cells remains poorly understood. Screening and identifying peptides or proteins that may interact with the ORF55 protein could help to elucidate its mechanism of action.

In this study, we cloned the ORF55 gene into a prokaryotic expression vector using seamless cloning, and used PDT to screen interacting partners. Bioinformatics prediction identified four potential interacting proteins. This study lays a foundation for further exploration of the biological functions of the ORF55 protein during viral infection.

Materials and methods

Bacterial strains, cloning vectors, cells and virus

Escherichia coli DH5α and E. coli BL21 (DE3) strains were obtained from TransGen Biotech (Beijing, China). The prokaryotic expression vector pMAL-c5X was purchased from New England Biolabs (NEB; Beijing, China). The CyHV-2 YC-01 strain (GenBank accession No. MN593216; 275,367 bp) was isolated and purified from Crucian Carp showing obvious symptoms of the early stages gill haemorrhage in the laboratory at Shanghai Ocean University, and preserved at − 80 °C. A cell line from the caudal fin of C. auratus gibelio (GiCF) [25] was established in previous experiments and cultured in M199 (Genomcell, Zhejiang, China) medium supplemented with 12% fetal bovine serum (FBS; Every Green, Zhejiang, China) at 27 °C. This cell line could be stably infected with CyHV-2 and used for in vitro virus replication. HEK-293 T cells were cultured in Dulbecco’s modified Eagle’s medium with high glucose containing 10% FBS in a temperature-controlled incubator at 37 °C with 5% CO2.

Sequence-based prediction

The National Center for Biotechnology Information (NCBI) database of the National Library of Medicine was used to obtain the sequence of the TK gene from the CyHV-2 YC-01 strain. The ExPASy ProtParam tool (, a server that calculates theoretical protein parameters including molecular weight, amino acids, isoelectric point (pI) and instability index, was used assess the protein physical and chemical properties [26]. In order to explore the functions of the ORF55 protein, we applied different in silico approaches using Cell-PLoc (, PSORTb (, lncLocator (, iLoc-LncRNA ( and Euk-mPLOC ( to predict localisation and other characteristics [27,28,29,30,31]. The secondary structure of the ORF55 protein sequence was predicted by various web servers including SOPMA (, NovoPro (, GOR ( and PSSpred ( [32,33,34,35]. Also, we downloaded the TK gene sequences of other species and performed evolutionary tree analysis using MEGA version 7.0. For homology analysis as well as sequence alignment, TK genes of all published CyHV-2 strains and representative CyHV-1 and CyHV-3 strains were selected for comparison. Protein tertiary structure was predicted using I-TASSER (

Seamless cloning

The ORF55 gene was amplified using PCR with sense and antisense primers 5′-GCGATATCGTCGACGGATCCATGGCTTTTCTGGAGTTGGTGCTGG-3′ and 5′-CCTGCAGGGAATTCGGATCCTCATACAAAAGAACAAGGGGCATCC-3′. Restrictions sites for BamHI were inserted in both primes. Thermal cycling included an initial denaturation step at 98 °C for 5 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 58 °C for 30 s, and extension at 72 °C for 30 s, and a final extension at 72 °C for 5 min. The subsequent ligation reaction was carried out at 50 °C according to the pEASY-Basic Seamless Cloning and Assembly Kit (TransGen Biotech), after which ligation products were transformed into E. coli DH5α. Positive colonies were screened and plasmid was isolated and sent to Sangon (Shanghai, China) for DNA sequencing.

Protein purification

Transformed E. coli BL-21 (DE3) cells were grown in Luria Bertani (LB) broth media containing 100 µg/mL ampicillin, and isopropylthiogalactoside (IPTG) was added to induce protein expression. Cells were collected by centrifugation and resuspended in 40 mL of 30 mM Tris–HCl, 20% sucrose, 1 mM EDTA, and mixed with 40 mL MgSO4. Following complete ultrasonic lysis and clarification by centrifugation, the bacterial supernatant was incubated with Amylose Resin (Sangon) on ice and transferred to an affinity chromatography column. After column washing, protein was eluted in elution buffer (20 mM Tris–HCl, 0.2 mol/L NaCl and 1 mM EDTA) containing 10 mM maltose.

Western blotting

The bacterial solution before and after induction, as well as the purified protein, were separated by 12% SDS-PAGE at 80 V for 30 min followed by 120 V for 60 min. Gels were stained with Coomassie Brilliant Blue. For western blotting, purified protein was subjected to electroblotting on a polyvinylidene difluoride (PVDF) membrane, which was subsequently blocked with 5% nonfat milk. The membrane was incubated with anti-MBP-Tag monoclonal antibody as primary antibody (Proteintech, USA) then with goat anti-mouse IgG (H + L)-horseradish peroxidase (HRP) conjugate (Bioworld Technology, USA) as secondary antibody. Results were detected using an enhanced chemiluminescent (ECL) detection system.

Phage display

The ORF55 protein concentration was determined from a standard curve prepared using bovine serum albumin (Sangon) and a Bicinchoninic Acid (BCA) Protein Quantification Kit (Beijing Solarbio Science & Technology Co. Ltd., Beijing, China). Protein was seeded in a 12-well plate at a concentration of 100 μg/mL, with three replicates for each experimental group. The phage library (NEB) was bound to the ORF55 protein according to the instruction manual of the Random Ph.D.-12 Phage Display Peptide Library Kit. Following three rounds of panning, the plate with plaques was sent to Genewiz Biotechnology Co. Ltd. (Tianjin, China) for sequencing and identification. The plaque sequencing were compared with NCBI crucian carp and goldfish databases.

Interaction validation

As described above, we obtained polypeptides that may interact with the ORF55 protein through PDT. The molecular docking tool from Home for researchers was used to predict interactions between the ORF55 protein and polypeptides, and to confirm their specific interaction sites. Dot blot experiments were also performed to verify interactions. Serially diluted peptides were blotted onto a PVDF membrane, soaked in ORF55 protein solution, and successively incubated with primary and secondary antibodies as described above for western blotting.

Protein prediction

Sequences of the four potential interacting proteins were downloaded from NCBI databases, and their sequences in other hosts were obtained for alignment, tertiary structure prediction, evolutionary tree analysis, and prediction of other physicochemical properties.

Preliminary in vitro validation

GiCF cells were infected with the supernatant of CyHV-2 YC-01 and sampled at various timepoints to observe whether the virus solution could alter expression of the above genes. RNA was extracted from sampled cells using TRIzol reagent and reverse-transcribed to cDNA as a template for gene amplification. Gene expression was measured using PrimeScript RT Master Mix (Takara, Beijing, China) from amplified cDNA. A real-time system instrument was used to quantify gene expression levels. Reactions contained 6.5 µL of TB Green Premix Ex Taq II (2 × ; Takara), 0.5 µL forward/reverse primers, and 500 ng cDNA. Reactions were repeated in triplicate.

Subcellular localisation analysis

Based on the above results, two proteins most likely to interact with the ORF55 protein were selected for further validation. Eukaryotic expression plasmids were constructed using seamless cloning technology as described above. Two luciferase reporter genes were selected as vectors. Unlike prokaryotic expression, a free endotoxin plasmid extraction kit was required for amplification. Recombinant plasmids were co-transfected into HEK-293 T cells. After 40 h, cell nuclei were stained using DAPI, and photographed using a confocal microscope.

Coimmunoprecipitation analysis

Samples of cells co-transfected as described above were collected to determine whether proteins were successfully expressed. Cells were lysed using NP40 lysate, and green fluorescent protein (GFP) fusion-tagged proteins and their interacting partners were selected by rabbit anti-GFP immunomagnetic beads. Western blotting was performed to probe interactions between proteins using anti-RFP Tag mouse monoclonal antibody as primary antibody as described above.

Statistical analysis

Statistical analysis was performed using GraphPad Prism ( and IBM SPASS statistics (SPASS 19.0). Results are presented as mean ± standard error (SE). Statistical significance was calculated using one-way analysis of variance (ANOVA). Values were considered significant (*) for p-values of 0.01 to 0.05, very significant (**) for p-values of 0.001 to 0.01, and extremely significant (***) for p-values 0 to 0.001.


Physiochemical parameters

The physicochemical properties of the CyHV-2 YC-01 ORF55 protein were predicted using the Protparam web tool (Table 1). The predicted pI of ORF55 was < 7, suggesting it is acidic. The instability index value indicated that it is an unstable protein in CyHV-2 and CyHV-3, but stable in CyHV-1. The calculated aliphatic index indicated that it was thermally over a wide temperature range. ORF55 was classified as polar based on the grand average hydropathy (GRAVY) value. Knowing where the ORF55 protein is localised is critical to predicting its functions. All different algorithm-based servers predicted that the ORF55 protein is in the cytoplasm (Table 2). Almost all servers predicted that the ORF55 protein is rich in the alpha-helix secondary structure, which is likely important for protein structure and function (Table 3).

Table 1 Physiochemical profile of ORF55 protein in different virus strain type
Table 2 Subcellular localization prediction of ORF55 protein
Table 3 Secondary structure profile of Rv1636 protein by the different webserver

Bioinformatics analysis of the ORF55 protein

A phylogenetic tree of TKs from different viruses was constructed using the neighbour-joining method in MEGA 7.0 software, revealing the evolutionary journey and relatedness of TK genes (Fig. 1A). Homology analysis in MegAlign 7.1 showed that all CyHV-2 strain types shared 100% sequence identity and > 50% identity with CyHV-1 and CyHV-3 (Additional file 2: Fig. 1G). All CyHV-2 strain types shared 98.6% similarity (Fig. 1B). Among them, ST-J1 and YC-01 strains were completely identical, with nine more bases (GGACGAGGA) than other strains. The results indicate that the ORF55 gene is highly conserved between CyHV-2 strains. The protein tertiary structure analysis (Fig. 1C) showed that there was a Pfam: TK domain located at 2–176, which laid a foundation for subsequent in-depth analysis of the reciprocal binding sites.

Fig. 1
figure 1

Bioinformatics analysis of the CyHV-2 ORF55 gene. A Evolutionary tree. Trees were generated with the neighbour-joining method using MEGA version 7. TK genes from different viruses were selected for simultaneous evolutionary tree generation, and TK from CyHV-2 is marked with a black square. Bootstrap values (1000 replicates) below 70% are hidden. B Amino acid sequences of TKs CyHV-1, CyHV-2 and CyHV-3. C Tertiary structure prediction. The tertiary structure of the ORF55 protein was generated by I-TASSER

Validation of protein expression

PCR results showed that the ORF55 gene was 645 bp in length, and 904 bp after connecting the prokaryotic expression vector, as expected (Fig. 2A). Expression of the ORF55 recombinant plasmid resulted in protein expression in E. coli BL-21 (DE3) following induction with 0.5 mM IPTG (Fig. 2B). SDS-PAGE results revealed a single band corresponding to MBP-tagged recombinant ORF55 protein, with molecular weights of 45 kDa for MBP and 70 kDa for ORF55 (Fig. 2C). Western blotting displayed a single band for purified recombinant ORF55 protein (Fig. 2D).

Fig. 2
figure 2

PCR amplification, prokaryotic expression, and validation of the ORF55 gene. A DNA extracted from the virus supernatant of infected cells was used as the amplification template, and the ORF55 gene was amplified by PCR with specific primers. M, DNA Marker 2000; 1, MBP gene (partial); 2, ORF55 gene; 3, MBP-ORF55 gene; 4, negative control. B Induced expression analysis of the ORF55 protein. The bacterial solution before and after IPTG induction was collected by centrifugation and analysed by SDS-PAGE. M, Protein markers; 1, uninduced; 2, induced. C SDS-PAGE analysis. The recombinant ORF55 protein was purified and the collected fractions were subjected to SDS-PAGE. The target protein was purified to a high concentration, and molecular weights of 45 kDa and 70 kDa were observed for bands corresponding to the MBP tag and ORF55 recombinant protein, respectively. M, Protein markers; 1, MBP protein; 2, ORF55 recombinant protein. D Western blotting analysis of purified ORF55 recombinant protein. A single band was observed. M, Protein markers; 1, MBP protein; 2, ORF55 recombinant protein

Screening of interacting polypeptides

A standard curve was established using BSA and the ORF55 protein concentration was measured at 562 nm using a spectrophotometer (Additional file 2: Fig. 2). During the three rounds of phage screening, the screening conditions were optimised by gradually reducing the amount of coated protein, increasing the concentration of Tween-20 and washing with detergent, and reducing the incubation duration with phage to obtain positive clones with high affinity (Table 4A). To determine phage enrichment for specific binding to the ORF55 protein in each round of screening requires a quantifiable indicator, in addition to the visual counting of phage plaques obtained in each round of phage panning (Fig. 3), hence phage titres after each round of screening were measured. The results showed that the correlation between the number of selection rounds and the titres of phages bound to ORF55 protein was positive (Table 4B). The recovery yield of phages was increased approximately 84-fold after the third round of panning compared with the first round. This indicates an obvious enrichment in the specific binding capability of phages.

Table 4 Enrichment, titer analysis and polypeptide screening of phages
Fig. 3
figure 3

Phage plaque plating. Determining the phage enrichment for specific binding to the ORF55 protein for each round of screening requires a quantifiable indicator, hence visual counting and observation of phage plaques were performed for all three rounds of phage panning. Each panning was repeated three times. A First round of screening. B Second round of screening. C Third round of screening

Prediction of interacting proteins

After the third round of screening, 50 plaque samples were randomly selected and sent to Genewiz for DNA sequencing. Based on the resulting data, nucleotide sequences were translated into amino acid sequences using DNAMAN (version 6). The results showed that the polypeptide sequence LSPGANSHVSRH had the highest homology, with 38 repetitions (Table 4C), suggesting it binds strongly to the ORF55 protein.

Validation of polypeptide-protein interactions

Through analysis of the phage display results described above, we obtained a polypeptide sequence of 12 amino acids. In this experiment, we used the dot blot method for validation. When the concentration of the ORF55 protein was > 0.5 mg/mL, the dot blot results showed a significant band, with the MBP protein tag alone serving as a negative control (Fig. 4A).

Fig. 4
figure 4

Prediction and verification of the interaction between ORF55 protein and polypeptide. A Dot blot results. Dot blot overlay analysis of ORF55 interaction with polypeptide in vitro. B Prediction of interactions. Chimera was used to construct 3D models, and Autodock 4.2 was used to perform protein-peptide docking

In order to further understand the role of the ICG (ligand) in nanostructure formation, ORF55 was used as a macromolecule and the polypeptide was used as a ligand for in silico analysis of the interaction between them during NP formation. The docking results revealed active binding between the ORF55 protein and the polypeptide. Specifically, binding between the two was due to interactions between a lysine residue of the ORF55 protein and a sulphonate group of the polypeptide (Fig. 4B). It was found that the docked ligand is close to the lysine and lies within hydrogen bonding distance of the sulphonate group. The distance from the lysine residue of PLL to the adjacent oxygen atom of ICG is in the range of 1.8 Å to 2.2 Å. Docking studies revealed that the docked ICG forms two hydrogen bonds with the lysine, with a bond length < 2.2 Å, suggesting possible aggregation of the ICG molecule within PLL/salt aggregates.

Bioinformatics analysis of potential interacting proteins

Amino acid sequences were analysed by alignment to the NCBI crucian carp and goldfish databases, and the ORF55 protein was predicted to interact with four proteins (Table 4D). They were gastrula zinc finger protein XlCGF8.2DB-like (ZFP, GenBank accession NO. XP_026094447.1), WD repeat-containing protein 7 isoform X1 (WDR7, XP_026052722.1), actin-binding Rho-activating protein (ABRA, XP_026140919.1) and zinc finger protein 516-like Isoform X1 (ZFP516, XP_026087732.1). A brief analysis of the respective physical and chemical parameters of the above proteins was performed (Additional file 1: Table S1). The Swiss-model ( was used to predict the tertiary structures (Fig. 5) and the evolutionary tree (Additional file 2: Fig. 3G) of the four proteins. The results showed that all four shared at least 70% sequence similarity between homologs in different species, indicating high conservation between their respective gene families. The results of conserved domain analysis showed that all four proteins had at least one domain or superfamily structure (Additional file 1: Table S2) potentially playing a specific role in the process of viral infection.

Fig. 5
figure 5

Tertiary structure analysis. Tertiary structures for five potential interacting proteins were predicted

RT-qPCR verification

GiCF cells were infected with the supernatant of CyHV-2 YC-01 and sampled at different timepoints (0 h, 12 h, 24 h, 36 h and 48 h). A significant CPE was observed under the microscope ~ 48 h after viral infection, and all cells had died ~ 144 h after infection (Fig. 6A). According to the nucleic acid sequences published by NCBI, RT-qPCR primers were designed (Table 5). The results showed that the relative expression of ZFP, ABRA and ZFP516 was upregulated in the latter stages of viral infection, without a significant change in WDR7 expression at any of the timepoints (Fig. 6B − E). Melting curves of RT-qPCR experiments are shown in Additional file 2: Fig. 4. Additionally, viral titres were tested, and the results showed a significant increase with increasing time, as predicted (Fig. 6F).

Fig. 6
figure 6

Images of virus-infected cells and RT-qPCR results. A Viral infected cells. GiCF cells were infected with the supernatant of strain YC-01 stored in our laboratory, and sampled at different timepoints (0 h, 12 h, 24 h, 36 h, 48 h, 54 h, 60 h, 66 h, 72 h, 78 h, 84 h, 90 h, 96 h, 108 h and 120 h). Over time, CPE of cells gradually became obvious until all cells died. BE RT-qPCR results. RNA was extracted by Trizol reagent and reverse-transcribed to cDNA to be used as template. The above genes were quantified using a real-time instrument. Error bars represent standard deviation of the mean for experiments performed in triplicate (*p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001). F Virus titre analysis. Titre values after viral infection of cells were determined at different timepoints

Table 5 Oligonucleotide primers and conditions for Real-time PCR

Subcellular localisation

Through the results of the relative expression of four genes described above, we inferred that ABRA might interact with the ORF55 protein, and attempted to validate this. Both pEGFP-N1 and pDsRed-Express-N1 were chosen to connect ORF55 and ABRA genes, using primers listed in Table 6. Both were cotransfected into HEK-293 T cells, and RFP empty vector served as a control. After 40 h, cells were transferred to a confocal microscope, observed, and photographed. The results showed that the ORF55 protein colocalised with ABRA in the nuclei of 293 T cells, but it did not colocalise with RFP, indicating that the RFP tag and ORF55 protein were located in different organelles, respectively. (Fig. 7).

Table 6 Primer sequences for Dual-luciferase reporter gene assay
Fig. 7
figure 7

Subcellular localisation. HEK 293 T cells were transfected with pEGFP-N1-ORF55, pDsRed-Express-N1-ZFP and pDsRed-Express-N1-ABRA for 40 h. RFP was used as a negative control


Cell samples after successful transfection as described above were analysed by western blotting. The experiment was divided into two parts; the first set of samples contained RFP and GFP-ORF55 protein, and the second set of samples contained RFP-ABRA and GFP-ORF55 protein. Both were recognised by GFP antibodies in the input (IP) group, while the RFP antibody only recognised 26 kDa and 43 kDa proteins, respectively. The IP group confirmed that the GFP magnetic beads had successfully adsorbed the GFP-ORF 55 protein. The final coimmunoprecipitation results revealed a clear band at 43 kDa in the second set of samples, implying an interaction between the GFP-ORF55 protein and the RFP-ABRA protein, but no interaction with the RFP protein (Fig. 8A).

Fig. 8
figure 8

Verification of the interaction between ORF55 protein and ABRA. A Molecular docking of ORF55 protein and ABRA. Home for researchers was used to analyse the specific interaction binding sites of the two proteins. B Results of coimmunoprecipitation. After taking cell samples for validation analysis, immunoprecipitation magnetic beads were used to purify and separate interacting proteins

Furthermore, we used molecular docking to verify the interaction between the ORF55 protein and ABRA, and to make specific binding site predictions. ZDOCK Score values and their best pose interactions were calculated (Table 7). PROK2 forms hydrogen bonds with amino acids including ARG167-ASP109 and ARG167-ASP112. Comprehensive analysis revealed that proteins PROK2 and CSF3 formed stable protein docking models (Fig. 8B).

Table 7 Results of molecular docking


Outbreaks of CyHV-2 with a mortality rate of > 90% and no effective treatment measures have been reported in many countries and regions. There are three virus types in the Alloherpesviridae family, among which CyHV-1 is the least harmful, hence it is rarely reported. CyHV-3, also known as Koi herpesvirus (KHV) [36], causes a large number of koi deaths worldwide, hence it has received international attention. There have been some studies on virulence genes of KHV [37], but research on CyHV-2 protein is limited, and the main focus has been on the identification and functions of major membrane proteins and capsid proteins, as well as immunological detection methods [38]. The function and virulence mechanism of the CyHV-2 virulence remain unknown.

In the field of virus research, the E. coli prokaryotic expression system has many advantages including relatively low cost and simple operation, hence it is widely employed, as demonstrated for Duck Tembusu virus (DTMUV) non-structural protein 2A (NS2A) [39]. The CyHV-2 protein received some attention [40]. Therefore, we wished to set up a prokaryotic expression system and purify recombinant ORF55 in the present work. According to the preliminary analysis of physicochemical parameters, the ORF55 protein was unstable and not highly hydrophilic (Table 1), indicating that the final purified protein may not be highly soluble and/or easily degraded. Based on these analyses, the pMAL-c5X vector was selected for gene cloning because MBP is known to be highly soluble [41] and can enhance the biological activity of fused target proteins [42]. We successfully cloned the ORF55 gene into the pMAL-c5X vector, and expressed and purified the MBP-ORF55 fusion protein.

We then explored the reference lists of relevant studies on potential interacting proteins (Table 4D). Zinc finger proteins (ZFPs) are a class of proteins that contains short, stable, self-folding ‘finger’ structures that bind zinc ions. Studies have shown that zinc-finger proteins are natural host endogenous immune factors that regulate expression of cytokines such as interleukins and interferons, and directly regulate lymphocyte growth and differentiation [43]. The human specificity protein 1 (SP1) transcription factor belongs to the C2H2-type zinc finger family, and has been widely studied in HeLa cells (human cervical cancer cells) [44]. This protein regulates gene expression by interacting with other proteins, acting as a negative or positive regulator of gene expression [45]. In addition, a class of zinc-finger antiviral proteins has been widely studied in land animal viruses, including Small Ruminant Morbillivirus (SRMV) [46], Hepatitis E virus (HEV) [47] and Human immunodeficiency virus [48]. Members of this class of proteins can have specific inhibitory effects on viral replication.

Actin-binding Rho-activating Protein (ABRA, also known as Stars) is a muscle-specific actin-binding protein, a novel and evolutionarily conserved actin-binding protein (ABP) that binds to I-Band in sarcomeres and actin filaments in transfected cells to activate Rho signalling events [49]. ABRA may be involved in regulating actin cytoskeleton functions and actin polymerisation in various cells. Actin is a globular multifunctional protein that forms microfilaments, present in essentially all eukaryotic cells, and is the most abundant cytoskeletal protein in eukaryotic cells. Actin plays a regulatory role in viral infection, and is closely related to virus adsorption and entry, virus replication, virus assembly, virus release, and virus transmission. Studies have shown that actins including Profilins [50], Cofilin [51, 52] and Drebrin [53] can play a key role in viral infection [54]. Rho family proteins were the first cloned proteins in the Ras superfamily; they have GTP hydrolysis activity and play an important role in the regulation of cytoskeletal recombination. Rho is highly expressed in various malignancies and is closely related to tumorigenesis, invasion and metastasis [55].

We also investigated WD repeat-containing protein 7 (WDR7). F-box and WD repeat domain containing 7 (FBXW7) is a member of the F-box protein family that functions as the substrate recognition component of the SCF E3 ubiquitin ligase. FBXW7 is a critical tumour suppressor and is able to control proteasome-mediated degradation of oncoproteins [56]. Herein, the WDR7 gene did not show an obvious change in expression, indicating that it might not be altered during viral infection of cells.

In this experiment, PDT was used to explore a polypeptides interacting with the ORF55 protein, and predict potential interacting host proteins. Putative peptides can then be tested in cells to assess inhibitory or promoting effects on the virus, and further explore its mechanism. Through molecular docking and in vitro experiments, we determined that ABRA could interact with the ORF55 protein. In future work, the mechanism of interaction will be further investigated. For example, RNA interference (RNAi) and overexpression can be used to determine the nature of interactions, and signalling pathways can be studied by transcriptome analysis. For the other three identified interacting proteins, other methods can employed to probe interactions with the ORF55 protein. This screening results provide a basis for exploring the specific biological functions of the ORF55 protein in the mechanism of virus infection.


In summary, we successfully constructed a prokaryotic expression system for the CyHV-2 YC-01 TK gene, and applied the resulting purified ORF55 protein in phage display experiments. A polypeptide sequence interacting with the ORF55 protein was successfully screened, and four potential interacting proteins were predicted using the NCBI database. Following validation by RT-qPCR and coimmunoprecipitation experiments, we determined that the ORF55 protein could interact with ABRA. The findings lay a foundation for the further investigation of the infection mechanism of ORF55.

Availability of data and materials

Not applicable.


  1. Jung SJ, Miyazaki T. Herpesviral haematopoietic necrosis of goldfish, Carassius auratus (L.). J Fish Diseases. 1995;18(3):211–20.

    Article  Google Scholar 

  2. Wu T, Ding Z, Ren M, An L, Xiao Z, Liu P, Gu W, Meng Q, Wang W. The histo- and ultra-pathological studies on a fatal disease of Prussian carp (Carassius gibelio) in mainland China associated with cyprinid herpesvirus 2 (CyHV-2). Aquaculture. 2013;412–413:8–13.

    Article  Google Scholar 

  3. Goodwin AE, Merry GE, Sadler J. Detection of the herpesviral hematopoietic necrosis disease agent (Cyprinid herpesvirus 2) in moribund and healthy goldfish: validation of a quantitative PCR diagnostic method. Dis Aquat Org. 2006;69:137–43.

    Article  CAS  Google Scholar 

  4. Jeffery KR, Bateman K, Bayley A, Feist SW, Hulland J, Longshaw C, Stone D, Woolford G, Way K. Isolation of a cyprinid herpesvirus 2 from goldfish, Carassius auratus (L.), in the UK. J Fish Dis. 2007;30:649–56.

    Article  CAS  PubMed  Google Scholar 

  5. Boitard PM, Baud M, Labrut S, de Boisseson C, Jamin M, Bigarre L. First detection of cyprinid herpesvirus 2 (CyHV-2) in goldfish (Carassius auratus) in France. J Fish Dis. 2016;39:673–80.

    Article  CAS  PubMed  Google Scholar 

  6. Giovannini S, Bergmann SM, Keeling C, Lany C, Schütze H, Schmidt-Posthaus H. Herpesviral hematopoietic necrosis in goldfish in Switzerland: early lesions in clinically normal goldfish (Carassius auratus). Vet Pathol. 2016;53:847–52.

    Article  CAS  PubMed  Google Scholar 

  7. Tang R, Lu L, Wang B, Yu J, Wang H. Identification of the immediate-early genes of cyprinid herpesvirus 2. Viruses-Basel. 2020.

    Article  Google Scholar 

  8. Xu J, Zeng L, Zhang H, Zhou Y, Ma J, Fan Y. Cyprinid herpesvirus 2 infection emerged in cultured gibel carp, Carassius auratus gibelio in China. Vet Microbiol. 2013;166:138–44.

    Article  PubMed  Google Scholar 

  9. Marsh M, Helenius A. Virus entry: open sesame. Cell. 2006;124:729–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kong SY, Jiang YS, Wang Q, Lu JF, Xu D, Lu LQ. Detection methods of Cyprinid herpesvirus 2 infection in silver crucian carp (Carassius auratus gibelio) via a pORF72 monoclonal antibody. J Fish Dis. 2017;40:1791–8.

    Article  CAS  PubMed  Google Scholar 

  11. Zhou Y, Jiang N, Ma J, Fan Y, Zhang L, Xu J, Zeng L. Protective immunity in gibel carp, Carassius gibelio of the truncated proteins of cyprinid herpesvirus 2 expressed in Pichia pastoris. Fish Shellfish Immunol. 2015;47:1024–31.

    Article  CAS  PubMed  Google Scholar 

  12. Davison AJ, Kurobe T, Gatherer D, Cunningham C, Korf I, Fukuda H, Hedrick RP, Waltzek TB. Comparative genomics of carp herpesviruses. J Virol. 2013;87:2908–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Graessmann A, Sandberg G, Guhl E, Graessmann M. Methylation of single sites within the herpes simplex virus tk coding region and the simian virus 40 T-antigen intron causes gene inactivation. Mol Cell Biol. 1994;14:2004–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Hay KA, Gaydos A, Tenser RB. The role of herpes simplex thymidine kinase expression in neurovirulence and latency in newborn vs adult mice. J Neuroimmunol. 1995;61(1):41–52.

    Article  CAS  PubMed  Google Scholar 

  15. Yi QY, Bai ZS, Cai B, Chen N, Chen LS, Yuan T, Mao JH. HSV-TK/GCV can induce cytotoxicity of retinoblastoma cells through autophagy inhibition by activating MAPK/ERK. Oncol Rep. 2018;40:682–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wen Y, Cheng A, Wang M, Ge H, Shen C, Liu S, Xiang J, Jia R, Zhu D, Chen X, et al. A Thymidine Kinase recombinant protein-based ELISA for detecting antibodies to duck plague virus. Virol J. 2010.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Gill MB, Kutok JL, Fingeroth JD. Epstein-Barr virus thymidine kinase is a centrosomal resident precisely localized to the periphery of centrioles. J Virol. 2007;81:6523–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Zhang L, Ruan K, Sang G, Xu Z, Tong W, Yu H, Shan T, Gao F, Li L, Kong N, et al. Tk-deleted pseudorabies virus retains high pathogenicity in rats. J Vet Res. 2021;65:401–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Zhao Y, Wang LQ, Zheng HH, Yang YR, Liu F, Zheng LL, Jin Y, Chen HY. Construction and immunogenicity of a gE/gI/TK-deleted PRV based on porcine pseudorabies virus variant. Mol Cell Probes. 2020;53:101605.

    Article  CAS  PubMed  Google Scholar 

  20. Smith GP. Filamentous fusion phage: novel expression vectors that display cloned antigens on the virion surface. Science. 1985;228:1315–7.

    Article  CAS  PubMed  Google Scholar 

  21. Dias RL, Fasan R, Moehle K, Renard A, Obrecht D, Robinson JA. Protein ligand design: from phage display to synthetic protein epitope mimetics in human antibody Fc-binding peptidomimetics. J Am Chem Soc. 2006;128:2726–32.

    Article  CAS  PubMed  Google Scholar 

  22. Mandrup OA, Friis NA, Lykkemark S, Just J, Kristensen P. A novel heavy domain antibody library with functionally optimized complementarity determining regions. PLoS ONE. 2013;8:e76834.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lee KJ, Lee JH, Chung HK, Choi J, Park J, Park SS, Ju EJ, Park J, Shin SH, Park HJ, et al. Novel peptides functionally targeting in vivo human lung cancer discovered by in vivo peptide displayed phage screening. Amino Acids. 2015;47:281–9.

    Article  CAS  PubMed  Google Scholar 

  24. Wang G, Yin P, Wang J, Ma P, Wang Y, Cai Y, Qi H, Liu A. Specific heptapeptide screened from pIII phage display library for sensitive enzyme-linked chemiluminescence immunoassay of vascular endothelial growth factor. Sensors Actuators B-Chem. 2021.

    Article  Google Scholar 

  25. Lu J, Xu D, Lu L. A novel cell line established from caudal fin tissue of Carassius auratus gibelio is susceptible to cyprinid herpesvirus 2 infection with the induction of apoptosis. Virus Res. 2018;258:19–27.

    Article  CAS  PubMed  Google Scholar 

  26. Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, Hochstrasser DF. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999;112:531–52.

    Article  CAS  PubMed  Google Scholar 

  27. Chou KC, Shen HB. Cell-PLoc: a package of Web servers for predicting subcellular localization of proteins in various organisms. Nat Protoc. 2008;3:153–62.

    Article  CAS  PubMed  Google Scholar 

  28. Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, Dao P, Sahinalp SC, Ester M, Foster LJ, et al. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010;26:1608–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Cao Z, Pan X, Yang Y, Huang Y, Shen HB. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018;34:2185–94.

    Article  CAS  PubMed  Google Scholar 

  30. Su ZD, Huang Y, Zhang ZY, Zhao YW, Wang D, Chen W, Chou KC, Lin H. iLoc-lncRNA: predict the subcellular location of lncRNAs by incorporating octamer composition into general PseKNC. Bioinformatics. 2018;34:4196–204.

    Article  CAS  PubMed  Google Scholar 

  31. Chou KC, Shen HB. A new method for predicting the subcellular localization of eukaryotic proteins with both single and multiple sites: Euk-mPLoc 20. PLoS ONE. 2010;5:e9931.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Geourjon C, Deléage G. SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput Appl Biosci. 1995;11:681–4.

    Article  CAS  PubMed  Google Scholar 

  33. Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292:195–202.

    Article  CAS  PubMed  Google Scholar 

  34. Sen TZ, Jernigan RL, Garnier J, Kloczkowski A. GOR V server for protein secondary structure prediction. Bioinformatics. 2005;21:2787–8.

    Article  CAS  PubMed  Google Scholar 

  35. Yan R, Xu D, Yang J, Walker S, Zhang Y. A comparative assessment and analysis of 20 representative sequence alignment methods for protein structure prediction. Sci Rep. 2013;3:2619.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Hedrick RP, Gilad O, Yun S, Spangenberg JV, Marty GD, Nordhausen RW, Kebus MJ, Bercovier H, Eldar A. A herpesvirus associated with mass mortality of Juvenile and Adult Koi, a strain of common carp. J Aquat Anim Health. 2000;12:44–57.;2.

    Article  CAS  PubMed  Google Scholar 

  37. Boutier M, Gao Y, Vancsok C, Suárez NM, Davison AJ, Vanderplasschen A. Identification of an essential virulence gene of cyprinid herpesvirus 3. Antiviral Res. 2017;145:60–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Shen Z, Jiang Y, Lu J, Sano M, Xu D, Lu L. Application of a monoclonal antibody specific for the ORF92 capsid protein of Cyprinid herpesvirus 2. J Virol Methods. 2018;261:22–7.

    Article  CAS  PubMed  Google Scholar 

  39. Jiang B, Zhang W, Liu P, Yang C, Wang M, Jia R, Zhu D, Liu M, Yang Q, Wu Y, et al. The prokaryotic expression, polyclonal antibody preparation, and subcellular localization of the transmembrane protein NS2A of the duck Tembusu virus. Acta Virol. 2020;64:380–5.

    Article  CAS  PubMed  Google Scholar 

  40. Haq WY, Kang SK, Lee SB, Kang HC, Choi YJ, Lee CN, Kim YS. High-level soluble expression of bioactive porcine myostatin propeptide in E coli. Appl Microbiol Biotechnol. 2013;97:8517–27.

    Article  CAS  PubMed  Google Scholar 

  41. Waugh DS. The remarkable solubility-enhancing power of Escherichia coli maltose-binding protein. Postepy Biochem. 2016;62:377–82.

    Article  CAS  PubMed  Google Scholar 

  42. Lee SB, Park SK, Kim YS. Maltose binding protein-fusion enhances the bioactivity of truncated forms of pig myostatin propeptide produced in E coli. PLoS ONE. 2017;12:e0174956.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Laity JH, Lee BM, Wright PE. Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol. 2001;11:39–46.

    Article  CAS  PubMed  Google Scholar 

  44. Li L, He S, Sun JM, Davie JR. Gene regulation by Sp1 and Sp3. Biochem Cell Biol. 2004;82:460–71.

    Article  CAS  PubMed  Google Scholar 

  45. Liu H, He J, Koh SP, Zhong Y, Liu Z, Wang Z, Zhang Y, Li Z, Tam BT, Lin P, et al. Reprogrammed marrow adipocytes contribute to myeloma-induced bone disease. Sci Transl Med. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Tang A, Tang J, Miao Q, Zhu J, Guo H, Liu C, Meng C, Li C, Chen Z, Liu G. Zinc finger antiviral protein (ZAP) inhibits small ruminant morbillivirus replication in vitro. Vet Microbiol. 2021;260:109163.

    Article  CAS  PubMed  Google Scholar 

  47. Yu W, Ji H, Long F, Chen S, He Q, Xia Y, Cong C, Yang C, Wei D, Huang F. Inhibition of hepatitis E virus replication by zinc-finger antiviral Protein synergizes with IFN-β. J Viral Hepat. 2021;28:1219–29.

    Article  CAS  PubMed  Google Scholar 

  48. Ficarelli M, Antzin-Anduetza I, Hugh-White R, Firth AE, Sertkaya H, Wilson H, Neil SJD, Schulz R, Swanson CM. CpG dinucleotides inhibit HIV-1 replication through zinc finger antiviral protein (ZAP)-dependent and -independent mechanisms. J Virol. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Arai A, Spencer JA, Olson EN. STARS, a striated muscle activator of Rho signaling and serum response factor-dependent transcription. J Biol Chem. 2002;277:24453–9.

    Article  CAS  PubMed  Google Scholar 

  50. Ohta K, Matsumoto Y, Nishio M. Profilin2 is required for filamentous actin formation induced by human parainfluenza virus type 2. Virology. 2019;533:108–14.

    Article  CAS  PubMed  Google Scholar 

  51. Lv X, Li Z, Guan J, Hu S, Zhang J, Lan Y, Zhao K, Lu H, Song D, He H, et al. Porcine hemagglutinating encephalomyelitis virus activation of the integrin α5β1-FAK-cofilin pathway causes cytoskeletal rearrangement to promote its invasion of N2a cells. J Virol. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Zheng K, Xiang Y, Wang X, Wang Q, Zhong M, Wang S, Wang X, Fan J, Kitazato K, Wang Y Epidermal growth factor receptor-PI3K signaling controls cofilin activity to facilitate herpes simplex virus 1 entry into neuronal cells. mBio 2014, 5, e00958–00913, .

  53. Li B, Ding S, Feng N, Mooney N, Ooi YS, Ren L, Diep J, Kelly MR, Yasukawa LL, Patton JT, et al. Drebrin restricts rotavirus entry by inhibiting dynamin-mediated endocytosis. Proc Natl Acad Sci U S A. 2017;114:E3642-e3651.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Borovac J, Bosch M, Okamoto K. Regulation of actin dynamics during structural plasticity of dendritic spines: Signaling messengers and actin-binding proteins. Mol Cell Neurosci. 2018;91:122–30.

    Article  CAS  PubMed  Google Scholar 

  55. Malliri A, Collard JG. Role of Rho-family proteins in cell adhesion and cancer. Curr Opin Cell Biol. 2003;15:583–9.

    Article  CAS  PubMed  Google Scholar 

  56. Yeh CH, Bellon M, Nicot C. FBXW7: a critical tumor suppressor of human cancers. Mol Cancer. 2018;17:115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


This work was supported by grants from the Key Research and Development Plan of Jiangsu Province (No. BE2021369) and the earmarked fund for CARS (CARS-45-16).

Author information

Authors and Affiliations



Conceptualization, H.W. and L. L.; Methodology, M.Q.; Software, M.Q. and S. X.; Validation, M.Q, S. X. and J.W.; Formal Analysis, M.Q.; Investigation, M.Q.; Data Curation, M.Q. and Y. Y; Writing – Original Draft Preparation, M.Q. and H.W.; Writing – Review & Editing, H.W.; Supervision, H.W.; Project Administration, H.W.; Funding Acquisition, H.W. and F. Y. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hao Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent to publication

Not applicable.

Competing interests

The authors declare that they have no competing of interest.

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.

Analysis of physical and chemical parameters and conserved domain of interactive protein.

Additional file 2.

Bioinformatic analysis and melting curves.

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Qian, M., Xiao, S., Yang, Y. et al. Screening and identification of cyprinid herpesvirus 2 (CyHV-2) ORF55-interacting proteins by phage display. Virol J 20, 66 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: