A Novel Vascular Endothelial Growth Factor Receptor Participates in White Spot Syndrome Virus Infection in Litopenaeus vannamei

Vascular endothelial growth factor (VEGF) signaling pathway is known to play key roles in endothelial cell proliferation, migration, angiogenesis, vascular permeability, inhibition of apoptosis, and virus infection. In the present study, a novel VEGFR gene (LvVEGFR2) was identified and characterized from Litopenaeus vannamei. The deduced amino acid sequence of LvVEGFR2 possessed typical features of VEGFRs reported in other species, including six IG-like domains, a transmembrane motif, a protein kinase (PK) domain, and one tyrosine-PK active site. The transcripts of LvVEGFR2 were mainly detected in hemocytes and lymphoid organ (Oka). Subcellular localization analysis showed that LvVEGFR2 was a membrane protein. Its expression level was obviously upregulated in hemocytes and Oka of the shrimp after white spot syndrome virus (WSSV) infection. Knockdown of LvVEGFR2 gene expression by double-strand RNA mediated interference could lead to a decrease of virus copy number in WSSV-infected shrimp. The interaction between LvVEGFR2 and different LvVEGFs (LvVEGF1, LvVEGF2, and LvVEGF3) in shrimp was analyzed at the transcription level and protein level, respectively. Knockdown of LvVEGF2 or LvVEGF3 could downregulate the expression level of LvVEGFR2, and injection of the recombinant LvVEGF2 or LvVEGF3 could upregulate the expression level of LvVEGFR2. Yeast two-hybrid analysis showed that LvVEGFR2 could interact with LvVEGF2 and LvVEGF3 directly. The study improved our understanding on the VEGF signaling pathway of shrimp and its role during WSSV infection.

high-affinity binding of VEGF ligands (3,4). Different VEGFRs play various functions through binding to different VEGF ligands. VEGFR1 plays roles in recruitment of hematopoietic stem cells, migration of monocytes and macrophages as well as in cellular energy metabolism by binding to VEGF-A and VEGF-B (5,6). VEGFR2 could regulate lymphatic endothelial cells and vascular endothelial cells, promote lymphatic and vascular production, as well as regulate lymphocyte migration after activation by VEGF-A and VEGF-D (7). VEGFR3 was considered to be involved in lymphangiogenesis and maintenance of the lymphatic endothelium by interacting with VEGF-C and VEGF-D (8,9). Vascular endothelial growth factor signaling pathway is also involved in the interaction between pathogens and host cells. The Parapoxvirus Bovine papular stomatitis virus and Orf virus could encode VEGF homologs which bind to VEGFR2 of the host cells and then participate in the process of viral infection (10)(11)(12). In crustacean, several VEGF genes have been identified and reported to be involved in the host-pathogen interaction. In Eriocheir sinensis, the VEGF homologous gene EsPVF1 could be activated by Vibrio anguillarum and Pichia pastoris GS115 (13). In Marsupenaeus japonicus, two types of VEGF, MjVEGF and MjVEGF2, were also involved in innate immunity (14). In our previous study, we found that the transcription levels of the genes in VEGF signaling pathway were obviously upregulated in shrimp Fenneropenaeus chinensis under acute infection of white spot syndrome virus (WSSV) (15). Subsequently, we isolated three VEGF genes (LvVEGF1, LvVEGF2, and LvVEGF3) and a receptor gene LvVEGFR from L. vannamei. Silencing of these genes by dsRNA-mediated RNAi could reduce the in vivo copy number of WSSV and the accumulative mortality rate of the infected shrimp decreased apparently (16)(17)(18). These data suggested that the VEGF and VEGFR genes played important roles in the host immunity. However, the knowledge about the VEGF signaling pathway in crustacean is still very limited.
In the present study, a novel VEGFR gene, LvVEGFR2, from the whiteleg shrimp L. vannamei was identified. The tissue distribution of LvVEGFR2 and its response to WSSV infection were analyzed. Double-strand RNA mediated RNA interference was performed to study the function of LvVEGFR2 gene during WSSV infection. Furthermore, the study also analyzed the interaction between LvVEGFR2 and each LvVEGF. The data will not only clarify the function of LvVEGFR2 gene during WSSV infection but also improve our understanding on the VEGF signaling pathway in the crustacean.

MaTerials anD MeThODs experimental animals
Healthy adult Pacific whiteleg shrimp cultured in our lab, with a body length of 16.1 ± 2.0 cm and a body weight of 12.1 ± 1.30 g, were used for tissue distribution analysis. Shrimp with a body length of 7.6 ± 0.5 cm and a body weight of 3.5 ± 1.0 g were used for WSSV infection, RNA interference, and recombinant protein injection experiments. All shrimp were acclimated in air-pumped circulating sea water at 25°C before experiment.

Tissue collection
For tissue distribution analysis, hemolymph of 12 healthy adult shrimp were obtained by using a sterile syringe preloaded with equal volume of modified anticoagulant Alsever solution and then centrifuged immediately at 800 g at 4°C for 10 min (19). The hemocytes were collected and preserved in liquid nitrogen. Then, the hepatopancreas, heart, muscle, lymphoid organ (Oka), intestine, stomach, gill, epidermis, eyestalk, testis, and ovary were dissected and kept in liquid nitrogen.
For WSSV challenge experiment, virus particles were prepared according to the method described by Sun et al. (20). The extracted WSSV was diluted to 500 copies/μl with sterile phosphate-buffered saline (PBS). Experimental animals were divided into 2 groups including WSSV group and PBS group, with 64 individuals in each group. 20 µl WSSV solution containing 10 4 copies of WSSV was injected into each shrimp in the WSSV group and an equal volume of PBS was injected into each shrimp in the PBS group. All shrimp were cultured in air-pumped circulating sea water with a temperature of 25 ± 1°C. Hemocytes and lymphoid organ (Oka) of 12 shrimp in each group were collected separately at 1, 3, 6, 12, and 24 hour (h) post WSSV infection (hpi).

Total rna extraction and cDna synthesis
Total RNA of each tissue was extracted using RNAiso Plus reagent (Takara, Beijing) following the manufacturer's protocol. The RNA concentration was assessed by Nanodrop 2000 (Thermo Fisher Scientific, USA) and the RNA quality was assessed by electrophoresis on 1% agarose gel. All cDNA samples were synthesized from 1 µg total RNA with PrimeScript RT Reagent Kit (Takara, Beijing). First, genomic DNA (gDNA) was removed with gDNA Eraser. And then, the first-strand cDNA was synthesized by PrimeScript RT Enzyme following the manufacturer's instructions with random primers.
cloning of the Open reading Frame (OrF) sequence of lvVegFr2 The ORF sequence of LvVEGFR2 cDNA was obtained from an Illumina-based transcriptome sequencing database of L. vannamei built in our lab (21). Primers LvVEGFR2-F1/LvVEGFR2-R1 and LvVEGFR2-F2/LvVEGFR2-R2 (Table 1) were designed to validate the sequence. The PCR program for the amplification of LvVEGFR2 was as follows: 1 cycle of denaturation at 94°C for 5 min, 35 cycles of denaturation at 94°C for 45 s, annealing at 57°C for 45 s, and extension at 72°C for 45 s, followed by an extension at 72°C for 7 min. The specific products were assessed by electrophoresis on 1% agarose gel. Then, the amplified products were purified using TIANgel Midi purification kit (Tiangen, China) and cloned into pMD19-T vector (Takara, Beijing) and then transformed into Trans 5α competent cell for sequencing.

sequence analysis
The complete ORF region and amino acid sequence of LvVEGFR2 was deduced using ORF finder. 1 The ORF sequence and deduced protein sequence of LvVEGFR2 were analyzed with BLAST   4 Eight protein sequences of VEGFR family members from different species were used to perform multiple alignments using ClustalW2 5 and a phylogenic tree was constructed by the neighbor-joining (NJ) algorithm using the MEGA 4 software. The reliability of the tree was tested by bootstrapping using 1,000 replications. The information of VEGFRs used in the present study was shown in Table 2.
Quantitative real-time Pcr Detection of lvVegFr2 mrna expression SYBR Green-based quantitative real-time PCR (qPCR) was performed to detect the gene expression level of LvVEGFR2.
A pair of primers LvVEGFR2-qF and LvVEGFR2-qR ( construction of pegFP-n1 Plasmid with the Transmembrane Motif of lvVegFr2 and Transfection into the Mammalian 293T cells A plasmid containing the predicted transmembrane motif (TM) of LvVEGFR2 was constructed to study the subcellular localization of LvVEGFR2 protein in mammalian 293T cells. The plasmid, designated as pEGFP-VR2, was constructed based on the expression plasmid pEGFP-N1 (Clontech, USA). The nucleotide sequence between the restriction sites Age I and Not I on pEGFP-N1 was replaced by a designed nucleotide sequence which encoded a IgK_secretion tag (METDTLLLWVLLLWVPGSTGD), the full length of enhanced green fluorescent protein (EGFP), and the predicted TM of LvVEGFR2 and its flanking sequence (from Val 644 to Ser 702 ), and flanked with the restriction sites Age I and Not I. The designed nucleotide sequence was synthesized by Sangon Biotech company (Shanghai, China) and sub-cloned into pEGFP-N1 using the restriction sites Age I and Not I. Three micrograms of the plasmids pEGFP-VR2 and pEGFP-N1 were transfected into the mammalian 293T cells with Lipofectamine 3000 Reagent (Thermo Fisher Scientific, USA) following the manufacture's instruction, respectively. The transfected cells were cultured for 24 h, fixed with 4% paraformaldehyde, and stained with 100 ng/ml DAPI (4′,6-diamidino-2-phenylindole) solution. The green and blue fluorescence signals were detected and merged.

Preparation and Optimization of Double-strand rna
To better understand the function of LvVEGFR2, we detected the effects of interference of endogenous LvVEGFR2 on the propagation of WSSV in L. vannamei. One pair of primers with T7 promoter sequences, LvVEGFR2-dsF and LvVEGFR2-dsR ( Table 1) were designed to amplify a 701-bp cDNA fragments of LvVEGFR2 as the template for dsRNA synthesis. Primers of dsEGFP-F and dsEGFP-R with T7 promoter sequences ( In order to optimize the silencing efficiency of dsRNA of LvVEGFR2, 30 healthy shrimp were divided into six groups. Shrimp in each group were injected with negative control or dsRNA with different dosages, including 1, 2, and 4 µg for each individual. At 48 h after interference, cephalothoraxes of four individuals in each group were isolated and frozen in liquid nitrogen for total RNA extraction. The transcription level of LvVEGFR2 was detected by qPCR with primers LvVEGFR2-qF and LvVEGFR2-qR as described in Section "Quantitative Realtime PCR Detection of LvVEGFR2 mRNA Expression. "

rna interference and WssV infection
In order to explore the effect of LvVEGFR2 silencing on WSSV propagation, 40 shrimp were divided into 2 groups including dsEGFP group and dsVR2 group. One microgram of dsRNA was injected into each shrimp after optimization of RNA interference effect. One microgram of dsEGFP and 1 µg dsLvVEGFR2 dissolved in 20-µl PBS was injected into each shrimp in dsEGFP group and dsVR2 group, respectively. At 48 h after dsRNA injection, 10 4 copies of WSSV dissolved in 10 µl PBS was injected into each individual. The pleopods of 15 individuals from each group were collected at 24 and 36 h after WSSV injection. The samples collected from each group at each time were divided into three subgroups and frozen in liquid nitrogen for DNA extraction and WSSV copy number detection.
DNA was extracted from the frozen pleopods using Genomic DNA Kit (Tiangen, China) following the manufacturer's instructions. The WSSV copy number in each DNA samples was quantified by qPCR with primers VP28-qF and VP28-qR ( Table 1) according to the method described by Sun et al. (20). Briefly, the plasmid DNA containing a 281 bp fragment of VP28 gene from WSSV was constructed and transfected into Escherichia coli. The plasmid was then extracted and quantified, and the copy number was calculated. Standard curve was constructed using 10-fold dilutions of the plasmid DNA ranging from 10 8 to 10 3 . The qPCR was performed with the diluted plasmid DNA and the pleopods DNA under the following conditions: denaturation at 94°C for 2 min; 40 cycles of 94°C for 20 s, 55°C for 20 s, and 72°C for 20 s. The WSSV copy number per nanogram of pleopods DNA was calculated based on the standard curve.

analysis of the effects of lvVegFs on the expression regulation of lvVegFr2
Detection on the Transcriptional Level of LvVEGFR2 after Silencing of LvVEGFs Previously, we identified three VEGF genes in L. vannamei (17,18). In order to explore the effect of LvVEGF genes on LvVEGFR2, the samples used in the previous studies, including shrimp at 48 h after RNA interference by dsRNA of LvVEGF1, LvVEGF2, and LvVEGF3 were used to detect expression of LvVEGFR2. QPCR with primers LvVEGFR2-qF and LvVEGFR2-qR ( Table 1) were performed as described in Section "Quantitative Real-time PCR Detection of LvVEGFR2 mRNA Expression. "

Detection on the Transcription Level of LvVEGFR2 after Injection of Recombinant LvVEGFs
Primers LvVEGF1-pF/LvVEGF1-pR, LvVEGF2-pF/LvVEGF1-pR, and LvVEGF3-pF/LvVEGF3-pR were designed to amply the ORF excluding sequence encoding signal peptide of LvVEGF1, LvVEGF2, and LvVEGF3. The PCR fragments were cloned into the pCzn 1 vector (Zoonbio Biotechnology, China) to construct recombinant plasmids. The recombinant plasmid pCzn 1-VF1, A Novel VEGFR Participates in WSSV Infection Frontiers in Immunology | www.frontiersin.org November 2017 | Volume 8 | Article 1457 pCzn 1-VF2, and pCzn 1-VF3 were transformed and successfully expressed in competent cells E. coli BL21 (DE3). Recombinant LvVEGFs protein pVF1, pVF2, and pVF3 were examined by SDSpolyacrylamide gel electrophoresis (SDS-PAGE), and purified by Ni-IDA-Sepharose CL-6B column (Bioz, USA). The purified pVF1, pVF2, and pVF3 protein were quantified and diluted to the concentration of 1.8 µg/µl. Twenty experimental shrimp were divided into four groups. Each group was injected with 36 µg pVF1, pVF2, pVF3, or pET30A dissolved in 20 µl PBS, respectively. In addition, 20 µl PBS was injected into shrimp as blank group. At 3 h after injection of pVF1, pVF2, or pVF3 protein, cephalothoraxes of four individuals in each group were collected and frozen in liquid nitrogen for total RNA extraction and cDNA synthesis as described in Section "Total RNA Extraction and cDNA Synthesis. " The transcriptional levels of LvVEGFR2 in these samples were detected through qPCR as described in Section "Quantitative Real-time PCR Detection of LvVEGFR2 mRNA Expression. " Yeast Two-hybrid assay

Yeast Transformation and Galactosidase Assays
The prey plasmid pGAD-VF1, pGAD-VF2, or pGAD-VF3 was co-transformed into yeast strain Y2H Gold with the bait plasmid pGBK-VR2(1-4) by the lithium acetate transformation procedure according to CLONTECH Matchmaker protocol manual (Clontech, USA). Transformants were then plated directly onto selective growth media, SD/-Leu/-Trp (DDO) plate. DDO is a commercial synthetic defined medium lacking both of leucine and tryptophan for selecting the bait and prey plasmids (Clontech, USA). Clones which contained both bait and prey plasmids could grew on the DDO plate. In addition, prey plasmids co-transformed with pGBKT7 and pGBK-VR2(1-4) co-transformed with pGADT7 were used to detect autoactivation. The co-transformed control plasmids pGBK-p53 and pGAD-T antigen were used as positive control, and pGBK-Lam and pGAD-T antigen were used as negative control, respectively. Transformants were allowed to grow at 30°C, usually for 3-5 days, until colonies were large enough (diameter 2-3 mm usually) for galactosidase activity assay.
After colonies grown, several fast-growing colonies were picked up and plated onto SD/-Leu/-Trp/X-α-Gal/AbA (DDO/ X/A) plate for primary screening. The plates were allowed to grow at 30°C, usually for 3-5 days. After 3-5 days, colonies were then picked and plated onto SD/-Ade/-His/-Leu/-Trp/X-α-Gal/AbA (QDO/X/A) plate. QDO is also a commercial synthetic defined medium lacking Adenine, histidine, leucine, and tryptophan. X-α-Gal and Aureobasid in A were added into QDO to make up QDO/X/A. Yeast colonies that express Mel1 turn blue in the presence of the chromogenic substrate X-α-Gal (Clontech, USA). This medium was used at the end of the two-hybrid screen to confirm protein interactions.

Detection on the expression levels of lvFaK and lvPi3K
In order to know whether LvVEGFR2 could regulate downstream signaling pathways, we detected the expression levels of LvFAK and LvPI3K after silencing of LvVEGFR2 gene. Primers LvFAK-qF/LvFAK-qR and LvPI3K-qF/LvPI3K-qR ( Table 1) were used to detect the expression level LvFAK and LvPI3K genes. The program was as follows: denaturation at 94°C for 2 min; 40 cycles of 94°C for 15 s, annealing temperature for 20 s, and 72°C for 20 s. The PCR product was denatured to produce melting-curve to check the specificity of the PCR product.

Data analysis
All the assays described above were biologically repeated for three times. For quantitative real-time PCR, four replicates were set for each sample. The relative transcription level of LvVEGFR2 were obtained using 2 −ΔΔCt method (22) and the WSSV copy number per nanogram of DNA was obtained according to the standard curve. In addition, the numerical data from each experiment were analyzed to calculate the mean and standard deviation of triplicate assays. An independent sample t-test was used to analyze the difference between two groups by SPSS 16.0. P < 0.05 was considered statistically significant.

The Full-length cDnas of lvVegFr2 and sequence analysis
The full-length ORF of LvVEGFR2 cDNA was 3,594 bp encoding 1,197 amino acid residues (accession number: MF417824). As shown in Figure 1, the deduced amino acids sequence of LvVEGFR2 contained five immunoglobulin subtype domains (Ser 48 to Phe 141 , Phe 156 to Lys 230 , Ala 243 to Lys 331 , Thr 349 to Leu 430 , and Asp 555 to Val 646 ) and one immunoglobulin-like domain (Pro 413 to Gln 552 ) were predicted in the extracellular region. TM of 21 aa was located from Asn 649 to Gly 671 followed by a protein kinase (PK) domain from Ile 726 to Leu 1007 . Moreover, a tyrosine-PK active site (Val 971 to Leu 983 ) existed in the PK domain. Phylogenetic analysis showed that LvVEGFR2 was first clustered together with LvVEGFR1 and then with VEGFR members from insects, while VEGFRs from vertebrates were clustered into another branch (Figure 2). Sequence alignment showed a 55% similarity between LvVEGFR1 and LvVEGFR2 amino acid sequences (data not shown).

Tissue Distribution and subcellular localization of lvVegFr2
The expression levels of LvVEGFR2 in different tissues were shown in Figure 3. The transcript of LvVEGFR2 was mainly detected in  Table 2).  hemocytes, followed by Oka and heart. The expression level of LvVEGFR2 was relatively lower in other tissues. Subcellular localization analysis showed that the fluorescence signals expressed by the constructed plasmid pEGFP-VR2 that contained the predicted transmembrane motif of LvVEGFR2 were mainly detected on the cell membrane of 293T cells, while the fluorescence signals expressed the control plasmid pEGFP-N1 were mainly detected in the cytoplasm of 293T cells (Figure 4).

silencing of lvVegFr2 affected In Vivo Virus Propagation
RNA interference based on dsRNA was performed to study the function of LvVEGFR2 gene. First, the silencing efficiency of LvVEGFR2 was detected under different dsRNA dosages, including 1, 2, and 4 µg for each shrimp. After optimization to the dsRNA dosage for interference, the dose of 4 µg dsRNA per shrimp, with 52.5% inhibition efficiency (Figure 6), was used for further RNAi experiment. The WSSV copy number was detected in shrimp after dsRNA and WSSV injection. At 48 and 72 hpi, the WSSV copy number in shrimp from dsEGFP group was about 4 × 10 5 copies/ng DNA, while it was about 1 × 10 5 copies/ng DNA in shrimp from dsVR2 group, which was significantly lower (P < 0.05) than that of dsEGFP group (Figure 7).

silencing of lvVegFs reduced the Transcription level of lvVegFr2
In order to know whether LvVEGFs could regulate the transcription of LvVEGFR2, LvVEGF1, LvVEGF2, and LvVEGF3 were silenced by dsRNA and then the transcriptional level of LvVEGFR2 was detected. After silencing of LvVEGF1 by 48.6%, the transcriptional level of LvVEGFR2 had no obvious change ( Figure 8A). After silencing of LvVEGF2 by 87.0%, the transcriptional level of LvVEGFR2 was downregulated by 55.3% (Figure 8B). After silencing of LvVEGF3 by 86.6%, the transcriptional level of LvVEGFR2 was downregulated by 59.3% ( Figure 8C).    Injection of the recombinant protein pVF1, pVF2, and pVF3 could also influence the transcriptional level of LvVEGFR2. As shown in Figure 9, injection of pVF1 did not upregulate the transcriptional level of LvVEGFR2 when compared to pET30A injection and PBS control. After injection of pVF2 or pVF3, the transcriptional level of LvVEGFR2 was obviously upregulated by 44 and 106% compared to pET30A injection and by 104 and 192% when compared with PBS control, respectively.  Detection on the interaction between lvVegFs and lvVegFr2 by Yeast Two-hybrid system In order to further investigate the interaction between LvVEGFs and LvVEGFR2, the yeast two-hybrid system was used. Expression of reporter genes in the positive control, which was co-transformed with pGBK-p53 and pGAD-T antigen plasmids, could generate blue color in the yeast cells (Figures 10A-C, zone 1). When pGAD-VF2 and pGBK-VR2(1-4) plasmids or pGAD-VF3 and pGBK-VR2(1-4) plasmids were co-expressed in the yeast cells, the blue color was also detected (Figures 10B,C, zone 2), while co-transformation of pGAD-VF1 and pGBK-VR2(1-4) plasmids into the yeast cells did not generate blue color ( Figure 10A, zone 2). These data indicated that LvVEGF2 and LvVEGF3 could interact with the Ig domains region of LvVEGFR2. Self-activation and negative controls were set and no blue signal was detected (Figures 10A-C, zones 3-5).
silencing of lvVegFr2 Decreased expression levels of lvFaK and lvPi3K genes In order to know whether LvVEGFR2 had impact on other signaling pathways, we detected the expression levels of LvFAK and LvPI3K in shrimp after LvVEGFR2 silencing. After silencing of LvVEGFR2 gene, the expression level of LvFAK and LvPI3K was downregulated by about 38 and 32%, respectively (Figure 11).

DiscUssiOn
The VEGF signaling pathway in mammals is considered essential in many physiological processes in angiogenesis, hematopoiesis and lymphogenesis (23,24). In crustacean, VEGF signaling pathway participates in innate immunity and is responsive to pathogen infection (13,14). Previously, we characterized a VEGFR gene (LvVEGFR1) and three VEGF genes in L. vannamei and data showed that they were all involved in WSSV infection (16)(17)(18). Here, we identified a new VEGFR gene, LvVEGFR2, in L. vannamei and characterized its function during WSSV infection. VEGFR is a kind of tyrosine kinase receptor which carries several immunoglobulin-like domains, a single transmembrane domain, a juxtamembrane domain, and a kinase domain split by a kinase insert and a carboxyl terminus (24). Sequence analysis A Novel VEGFR Participates in WSSV Infection Frontiers in Immunology | www.frontiersin.org November 2017 | Volume 8 | Article 1457  showed that LvVEGFR2 shared similar intracellular domains with VEGFR family members. In the extracellular region, most VEGFRs such as LvVEGFR1, HsVEGFR1, HsVEGFR2, and DmVEGF-R had seven IG-like domains, while the extracellular region of LvVEGFR2 was composed of six IG-like domains which were similar to HsVEGFR3 and DmVEGFR (1,2,25). Although phylogenetic analysis showed that LvVEGFR2 was first clustered with LvVEGFR1, the sequence similarity between them was only 55%. These data indicated that LvVEGFR2 was a novel VEGFR gene which was distinguished from LvVEGFR1 in L. vannamei. In order to know whether LvVEGFR2 was involved in immune response in shrimp, we first detected its spatial distribution and its temporal expression during WSSV infection. LvVEGFR2 showed the highest expression levels in hemocytes and Oka. Hemocytes play key roles in crustacean innate immunity such as phagocytosis (26,27), storage and release of the proPO (28), and so on. Oka is one of the most specific and effective organ for clearance of bacteria (29) and virus (30). At 24 h post WSSV infection, the expression level of LvVEGFR2 gene in hemocytes and Oka was significantly upregulated. Different from LvVEGFR2, the previously reported LvVEGFR1 gene was upregulated at 1 h post WSSV infection in hemocytes (16). The data indicated that LvVEGFR1 might be an early response gene to WSSV infection, while LvVEGFR2 might be a late response gene. It should be also noted that there was a significant decrease of LvVEGFR2 expression at 1 hpi in Oka, which was similar to the decrease of LvVEGFR1 at 1 hpi in intestine (16). A possible reason for the decreases might be that WSSV infection recruited VEGFs from other tissues into hemolymph since WSSV would first enter the hemolymph after muscle injection.
Interference of LvVEGFR1 gene by double-strand RNA could inhibit the in vivo propagation of WSSV in shrimp (16). Silencing of other members in shrimp VEGF signaling pathway, including LvVEGF1, LvVEGF2, and LvVEGF3 genes, could also reduce WSSV copy number in WSSV-infected shrimp (17,18). Similarly, silencing of LvVEGFR2 could obviously reduce the virus copy number in WSSV-infected shrimp. These data collectively supported the previous conclusion that inhibition of the VEGF signaling pathway could restrict the propagation of WSSV in shrimp.
The VEGF signaling pathway is directly activated by binding VEGFs to their receptors VEGFRs. The extracellular IG domains of VEGFRs are considered as the binding site of VEGFs (3,4). In shrimp L. vannamei, silencing of LvVEGF1 and LvVEGF2 genes with dsRNA reduced the expression level of LvVEGFR1 gene (17), and LvVEGF3 could interact with the second to fifth Ig domains of LvVEGFR1 under yeast two-hybrid system (18), which indicated that LvVEGFR1 might be the receptor for LvVEGF1, LvVEGF2, and LvVEGF3. In the present study, the transcription level of LvVEGFR2 was downregulated after silencing of LvVEGF2 and LvVEGF3, and upregulated by injection of recombinant LvVEGF2 and LvVEGF3 proteins. However, silencing of LvVEGF1 and injection of recombinant LvVEGF1 did not influence the expression level of LvVEGFR2 gene. Meanwhile, LvVEGF2 and LvVEGF3, rather than LvVEGF1, could interact with the extracellular Ig domains of LvVEGFR2. These data indicated that the LvVEGFR2 was also the receptor of LvVEGF2 and LvVEGF3.
After activation, VEGFRs usually function through regulating many genes in cellular signaling pathways such as STAT (31), PI3K/Akt (32), FAK (33), ERK (34), and so on. These genes or related signaling pathways play important roles during WSSV infection in crustaceans. WSSV infection could activate STAT (35) and then use it to promote WSSV ie1 gene expression (36). FAK could not only promote WSSV infection but also participate in immune defense (37). Suppression of PI3K-Akt-mTOR pathway inhibited in vivo WSSV propagation (38). Presently, the expression level of LvFAK and LvPI3K was downregulated after LvVEGFR2 silencing, which was similar to our previous reports on the LvVEGFR1 and LvVEGFs genes (16)(17)(18). Therefore, we inferred that the VEGF signaling pathway in shrimp affected WSSV propagation might through regulating downstream signaling pathways.
In conclusion, a novel VEGFR gene (LvVEGFR2) was identified in shrimp L. vannamei in the present study. LvVEGFR2 was responsive to WSSV infection and could regulate the in vivo propagation of WSSV. LvVEGFR2 interacted with extracellular ligands LvVEGF2 and LvVEGF3, and they could regulate expression of downstream genes including LvFAK and LvPI3K. These results, together with our previous report on LvVEGFR1 and LvVEGFs, collectively indicated that the VEGF signaling pathway in shrimp might participate in WSSV infection through regulating its downstream signaling pathways.
aUThOr cOnTriBUTiOns SL, FL, and JX conceived and designed the project. SL and ZW performed all the experiments and prepared the figures. KY prepared for the experimental animals. SL, ZW, and FL wrote the manuscript. All authors reviewed the manuscript.