The experimental shrimp were cultured in the laboratory of the Institute of Oceanology, Chinese Academy of Sciences (Qingdao, Shandong, China), at a temperature of 25 ± 1°C, salinity of 30%, and pH of 7.5 ± 0.1, the photoperiod was maintained at 12L:12D. The aquaculture seawater was filtered, sterilized, and continuously oxygenated. All shrimp were fed three times per day at 9:00 a.m., 2:00 p.m., and 7:00 p.m. with equal weights of commercial food pellets (Dale Feed Company, Yantai, China). The average weight of the shrimp was 4.0 ± 0.8g (7.5 ± 0.5cm). The animal study was reviewed and approved by the ethics committee of the Institute of Oceanology, Chinese Academy of Sciences. We declare that all animal experiments in this study were conducted in accordance with the guidelines of UK Animals Act, 1986 and EU Directive 2010/63/EU. In these experiments, no any endangered or protected species were used.

We collected all genes annotated as Ras superfamily or small G protein genes from the L. vannamei genome database (http://www.shrimpbase.net). At the same time, we screened all Ras superfamily genes in the previous RNA-Seq data from different developmental stages, molting stages, different adult tissues, and WSSV infection of the shrimp (Wei et al., 2014; Gao et al., 2015; Wang et al., 2019; Zhang et al., 2019). After comparing all the obtained sequences, eliminating redundant sequences, merging overlapping fragments and connecting broken genes, all non-redundant candidate sequences were initially identified and compared by blastx (https://blast.ncbi.nlm.nih.gov/Blast.cgi) and SMART (https://smart.embl.de/) to confirm the gene members of Ras superfamily.

In order to determine evolutionary relationships of the Ras superfamily members in L. vannamei, all identified Ras GTPase homologous sequences of this species and a number of representative sequences of other animals, fungi, protists and plants were aligned by ClustalW in MEGA X (www.megasoftware.net/), and phylogenetic trees were constructed using the maximum likelihood (ML) method of MEGA X with 1000 bootstrap repeats (Whelan and Goldman, 2001; Kumar et al., 2018). The phylogenetic trees were then visualized using iTOL (Letunic and Bork, 2007).

In order to illustrate the gene structure of the Ras GTPases of L. vannamei, the gene location, gene length, open reading frame (ORF), exon number and deduced amino acid number of all the obtained genes were analyzed in detail. ExPASy (http://web.expasy.org/protparam/) was used for protein molecular weight, isoelectric point analysis. Using the gene structure display server (GSDS) program (Hu et al., 2015), the coding sequence (CDS) of each gene was compared with the genome sequence of L. vannamei, and the intron and exon arrangement diagram of each gene was obtained. The gene structure of Ras GTPase proteins was drawn by TBools software (Chen et al., 2020), and the conservative domains analysis was carried out according to the amino acid sequences. In order to better analyze the conserved motifs in Ras superfamily, a conserved motif map using the Multiple EM for Motif Elicitation (MEME, https://meme-suite.org/meme/tools/meme) program was constructed to provide more detailed evidence for clarifying the structural characteristics between different categories of Ras superfamily numbers.

According to the previous RNA-Seq data at different developmental stages, molting stages, different adult tissues and WSSV infection of L. vannamei (Wei et al., 2014; Gao et al., 2015; Zhang et al., 2019), the RPKM (Reads per kilo base per million mapped reads) value of each unigene or transcript were obtained. The expression heatmaps of Ras GTPase genes in different transcriptomes were constructed and a part of representative genes were verified and analyzed using RT-qPCR subsequently.

The healthy WSSV-free shrimp L. vannamei (9 ± 1g) were collected from laboratory culture tanks for WSSV infection study, each shrimp was injected into 1000 copies of live WSSV particles suspended in 10μl sterile phosphate-buffered saline (PBS) in vivo WSSV challenge group, the control group shrimps were injected into the same volume of PBS. A total of 15 individuals were randomly assigned in each group and equally divided into three parallel subgroups as biological replicates. At 6 hpi, same size shrimp were picked, sacrificed and dissected. A total of 3 tissues, hemocyte, lymphoid (Oka) organ, hepatopancreas were sampled and frozen in liquid nitrogen and stored at -80 °C for total RNA extraction.

About 4.0 ± 0.8g (7.5 ± 0.5cm) untreated shrimp were picked, sacrificed and dissected. A total of 12 tissues, hemocyte, muscle, intestines, stomach, lymphoid organ, gill, hepatopancreas, eye stalk, brain, ventral nerve, epidermis and heart, were sampled and frozen in liquid nitrogen and stored at -80 °C for total RNA extraction. According to the manufacturer’s instructions, RNAiso Plusreagent (TaKaRa, Japan) was used to isolate total RNA from these tissues. Then, the quality and the concentration of RNA were detected using 1% agarose gel electrophoresis and Nanodrop 2000 (Thermo Fisher Scientific, United States). The first-strand cDNA was synthesized with the PrimeScript First Strand cDNA Synthesis Kit (TaKaRa, Japan) using 1.5 µg RNA as template. The specific cDNA synthesis were as follows: the first step was to remove genomic DNA, 5×genomic DNA eraser buffer was added to the template at 42°C for 5 min, and then moving on to step 2 immediately, at 37°C for 1 h and 85°C for 5 s. Finally, the cDNA was stored at -80°C until use.

The 18S rRNA gene was selected as the internal reference gene. Primer3Plus (http://www.primer3plus) was used to design thirteen pairs of primers of the selected Ras superfamily genes. All primer sequences used in this study are shown in Supplementary Table S1. Then, an Eppendorf Mastercycler ep realplex (Eppendorf, Hamburg, Germany) was used to perform the RT-qPCR. The SuperReal PreMix Plus (SYBR Green) (TIANGEN, Beijing, China), template, primers, and DEPC-treated water were mixed in a certain proportion (Supplementary Table S2); each sample includes four technical replicates. After qPCR, melting temperature (Tm) analysis showed a single peak and a single PCR band was identified, indicating that both primers were suitable for further experiments. The qPCR steps were as follows: 94°C for 2 min, 40 cycles of 94°C for 20 s, 62°C or 55°C for 20 s (the annealing temperature of rtRas superfamily genes-F/R and rt18S-F/R were 62 and 55°C, respectively), and 72°C for 20 s. Eventually, relative expression levels were analyzed by the 2^-△△Ct method (Livak and Schmittgen, 2001).

The different groups were subjected to one-way ANOVA tests using SPSS (https://www.ibm.com/cn-zh/analytics/spss-statistics-software) (version 20).

In this study, a total of 108 Ras superfamily genes of L. vannamei were screened and identified from genome and transcriptome data. The identified genes were classified and analyzed according to gene homology, conserved domain, gene annotation and evolutionary relationship. Using the SMART tool, the domains of 50 sequences were identified as the classical domains of the Ras superfamily: Ras, Rab, Arf (Sar), Rho respectively, and 49 sequences contained a small-GTPase domain, which have similar GTP binding conserved motifs to the classical Ras superfamily domains (Table 1), all them shared a set of conserved G box GDP/GTP-binding motif elements (Figure 1). Among the small-GTPase domain containing sequences, five have large molecular weight, they have not only a small-GTPase domain, but also other domains, such as BTB, FF, EFh, and RPT1 (Table 1). The remaining 9 sequences do not have the classical Ras superfamily domains or the small-GTPase domain, but they all have GTPase-related domains, and all the residues are perfectly conserved with Ras GTPases, mainly G1 and G3 motif (Figure 2A, corresponding to motif symbol 1 and 2).

From the perspective of gene location, these Ras superfamily genes were distributed in 98 scaffolds, Ras, Rho, Rab and Arf families each had 1-2 two-gene clusters, and Arfs also had a three-gene cluster (LvArf4, LvArf14, and LvArf15) (Table 1).

In order to distinguish these Ras superfamily members, we compared identity as assayed by smartblast and blastp, and found that most of them were classical Ras GTPases, such as Ras, Rab, Arf, Rho, Ran, named LvRas1-24, LvRab1-35, LvArf1-20, LvRho1-10, LvRan. A small number (18) are unconventional or newly discovered members of the Ras superfamily, named LvRbj, LvRGK1-3, LvGPN1-3, LvREM, LvIFT1-2, LvRGBP1-2, LvOBG1-2, Lvseptin, LvEFCAB4B, LvSRP, LvGPo (Table 1).

In order to clarify the phylogenetic relationship among Ras superfamily members of L. vannamei, we constructed phylogenetic trees using Ras superfamily genes from L. vannamei and those identified from other species (Table 1 and Supplementary Table S3). Phylogenetic analysis showed that the Ras GTPase proteins were clustered in the clades of each subfamily, although they come from different species (Figure 2A). The result suggested that the Ras superfamily members of L. vannamei could be divided into two large monophyletic groups. The Ras family and Rho family could be clustered into a clade, the number of Rab family is the largest (35 members) and form a clade with Ran family, then clustered with Ras family and Rho family into a larger clade. Arf family members were clustered together and far away from other four families, other unconventional Ras superfamily numbers were clustered together, and then clustered with the Arf family clade formed another large clade (Figure 2A).

According to the analyses of gene structure and conserved motifs, most Ras superfamily genes share a set of conserved G box GDP/GTP-binding motif elements: GXXXXGKS/T, T, DXXGQ/TE, NXXD and SXK (X stands for any amino acid), a few genes contain only G1 box, G3 box, and G4 box (Figures 1, 2B). However, different Ras GTPases have many specific motifs, resulting in different domains (Figure 2), and the differences of their amino acid sequences mainly occur in the amino terminal and carboxyl terminal, which are considered as major protein modification sites (Figure 1).

In order to clarify the expression patterns of Ras superfamily genes of L. vannamei, we summarized four different transcriptional profiles of 108 Ras GTPase genes into heatmaps of different tissues, development stages, molting stages and WSSV infection status (Figures 3–6; Supplementary Figures S1-6).

In different tissues of adult shrimp, most members of Ras superfamily showed low expression, a few genes were highly expressed in their respective families, such as LvRab18, LvRab20, LvRho3, LvRho5, LvRas4, LvRas10, LvArf1, LvArf10, and LvRan. (Figure 3). LvRan had the highest expression level in all the examined tissues, and its expression level is more than 3-5 times that of other members (Figure 3; Supplementary Figures S1, S2A). Other highly expressed Ras superfamily genes usually presented in some specific tissues: LvRas4 was highly expressed in hemocyte and intestines; LvRab18 was highly expressed in muscle; LvArf1 was highly expressed in brain, thoracic ganglion, ventral nerve and intestines; and LvRho5 was highly expressed in antenna, hepatopancreas and intestines. Some unconventional Ras superfamily members were highly expressed in testis and ovary tissues, such as LvGPN1, LvGPN3, LvIFT2 and LvOBG1.

At early development, most members of Ras superfamily were expressed at a higher level. LvRas1, LvRas2, LvRas5 and LvRas10 are the main expressed Ras family genes, in addition to LvRas10 at gastrula (gast) and limb bud embryo (Lbe) stages, the expression of Ras genes was not very obvious in early development as a whole. The Rab family showed two expression patterns, some genes were highly expressed mainly during the zygote to gast stages, while others were highly expressed after limb bud embryo stages. Arf12, Rho1 and Rho5 were the most expressed genes in their respective families, and their expression trends remained relatively stable during different stages (Figure 4; Supplementary Figure S3). Ran was the gene with the highest expression of all Ras superfamily numbers in the early developmental stages of L. vannamei, but with little fluctuations (Supplementary Figure S4A). Interestingly, most of the unconventional Ras superfamily genes had high expression levels at early development stages. LvSRP and LvGPo were highly expressed at the whole stages, showing that play an important role in early development. LvREM, LvGPN2, LvIFT2, and LvRGBP1 had similar expression patterns, they were all highly expressed after gastrula stage (gast). In contrast, LvRbj and LvGPN1 were highly expressed before gastrula stage (gast) (Supplementary Figures S4B-H).

In different molting stages of L. vannamei, most members of Ras superfamily showed lower expression (Figure 5; Supplementary Figure S5), and only a few genes were highly expressed in respective different families, such as LvRab18, LvRho1, LvArf12, LvRan, LvRas5. Similarly, a few unconventional members were highly expressed at different molting stages, such as LvRGBP1, LvSRP and LvGPo (Supplementary Figures S6B-H).

The hemocyte, hepatopancreas, lymphoid (Oka) organs are the three main immune related tissues. The expression level of many Ras superfamily genes was up-regulated in different organs after WSSV infection. In lymphoid organs, 19 Rab genes, 8 Rho genes, 11 Ras genes, 10 Arf genes, LvRGBP2, LvSRP, LvEFCAB4B and LvGPo were significantly up-regulated. In hepatopancreas, 11 Rab genes, 5 Rho genes, 4 Ras genes, 10 Arf genes, LvIFT1 and LvGPN2 genes were significantly up-regulated. In hematocyte, 6 Rab genes, 1 Rho gene, 7 Ras genes and 8 Arf genes were up-regulated, most unconventional Ras superfamily genes were also significantly up-regulated (except for LvRGKs, LvSRP and LvRGBP2) (Figure 6; Supplementary Table S2). The types of up-regulated Ras superfamily genes were different in different tissues, and it is possible that these genes play different roles in the immune process. For example, LvRan was up-regulated in hemocyte and hepatopancreas, but down-regulated in lymphoid organs after WSSV infection. In lymphoid organs, most of the Ras superfamily genes were up-regulated and a few genes were down-regulated, while was opposite in hepatopancreas. Among them, some Ras and Arf genes and LvRGK showed no changes in their expression levels in the three organs after WSSV infection, indicating that these genes are not affected by immune response.

In order to verify the accuracy of the expression level of different tissues, the expression of twelve Ras superfamily genes were verified by qRT-PCR in twelve different adult tissues (Figure 7). LvRas5 and LvRas10 were highly expressed in hepatopancreas and ventral nerve respectively, and LvRas6 was highly expressed in ventral nerve, lymphoid organ and brain. LvRab18, LvRab20, LvRho1 and LvRho5 were relatively expressed in most tissues except eyestalk, epidermis and muscle, among the expression tissues, these genes were mainly highly expressed in ventral nerve, brain, lymphoid organ and hepatopancreas, which is similar to the expressed pattern of verified Ras genes above. LvArf1 and LvArf12 were expressed in most tissues, and the three tissues with the highest expression were brain, ventral nerve and gill. LvRan was highly expressed in ventral nerve, and is about 2-4 times higher than other tissues. In the unconventional Ras superfamily members, LvGPN1 and LvRGBP1 had similar expression patterns with the classical Ras superfamily members, which were mainly highly expressed in hepatopancreas, brain, heart and other tissues. In sum, the most expression levels of these genes were consistent with respective RNA-Seq results, and some differences might be caused by individual differences. We found that almost all verified members of Ras superfamily were expressed in almost all detected tissues, and showed the lowest expression level in eyestalk, epidermis and muscle tissues. It is possible that Ras GTPases play a weak role in these tissues.

Ras family is the most typical and common class of Ras superfamily. In L. vannamei, multiple Ras family members such as Ras, Rap, Ral, Rheb, RHEs and RIT were found. Gene expression analysis showed that these Ras genes were expressed in different tissues, indicating their respective functions. Most members of the Ras family in mammals are called oncogenes and well-studied. Ras proteins receive signals from cell surface receptors, these signals are transmitted among proteins through different pathways, and finally affect a variety of biological functions, such as development, proliferation, differentiation, and survival (Goitre et al., 2014). If the Ras gene was mutated, these signal pathways will be destroyed, which will lead to a variety of tumors and cancers (Bos, 1989; Goodsell, 1999; Murugan et al., 2019). In addition, some Ras family genes are considered to be involved in animal growth and reproduction. For example, the Ras gene is related to the size of catfish head (Geng et al., 2016), it is involved in the regulation of insulin pathway during oocyte vitellogenesis in female oysters (Jouaux et al., 2012). During vitellogenesis of the marine flounder Solea senegalensis, the expression of Ras homologous gene was up-regulated (Tingaud-Sequeira et al., 2009). In L. vannamei, Ras family genes also were highly expressed in ovary, such as LvRas4 and LvRas10. In previous study, we conducted growth candidate gene association analysis in two independent populations of L. vannamei, the results showed that the SNP of Ras related protein gene Rap-2a (LvRas1) was significantly correlated with growth traits (Yu et al., 2019). The above research proves that Ras family may play an important role in growth and reproduction of shrimp.

There are also other classic Ras family genes in L. vannamei, such as M-Ras and Ras related protein R-Ras2. M-Ras alone constitutes a small Ras subfamily in mammals, in L. vannamei, M-Ras gene was named as LvRas2, and showed high relative expression under most conditions, indicating that M-Ras plays an important role in shrimp. Classical Ras binds to Raf and activates ERK pathway (Endo, 2020), while M-Ras mediated cell transformation related to the weak activation of RAF/MEK/ERK pathway, other downstream effectors are also involved in the pathway, which can induce neuronal differentiation (Castro et al., 2012). In mice, both classical Ras and M-Ras are highly expressed in the central nervous system (Sun et al., 2006). In addition, neuronal differentiation of p12 cells in rats also requires the induction of classical Ras and M-Ras (Sassone-Corsi et al., 1989). R-Ras2 has also been shown to be essential for correct axonal myelination and accurate neurotransmission in mouse (Gutierrez-Erlandsson et al., 2013; Sanz-Rodriguez et al., 2018). In L. vannamei, both M-Ras (LvRas2) and R-Ras2 (LvRas6) gene were highly expressed in the thoracic ganglion, which further showed that Ras family genes may play an important role in nervous system development and signal transduction.

At present, there are few studies on Ras family of L. vannamei, but as these typical Ras superfamily proteins may play an important role in metabolism, growth, reproduction and neural development, the related mechanisms need to be further studied.

Rab is a kind of regulatory small molecule GTPase protein, found on eukaryotic cell and organelle membrane (Seabra et al., 2002), and it also is a group with the largest number of genes of the Ras superfamily in L. vannamei. For example, there are 57 Rab genes in Arabidopsis, 30 in D. melanogaster, and 61 in Mus musculus (Rojas et al., 2012). The main members of Rab family include Rab1, Rab3A and Rab5c.

Rab, called Ypt in yeast, plays an important role in vesicle transport (Takai et al., 2001). Rab8 and Rab5 in human have the same function as YPT in yeast, and they share high sequence similarity (Molendijk et al., 2004). It has also been reported that Rab protein plays a role in vesicle transport in plants. For example, Rab1, Rab2, Rab4, Rab5 and Rab6 in Arabidopsis have been proven to play an important role in regulating the morphology of endoplasmic reticulum, Golgi apparatus and plasma membrane vesicles (Stenmark and Olkkonen, 2001). In the squid Loligo pealei, it was verified that Myo5a and Rab3A directly bind and interact with synaptic vesicles (SVs) and participate in the transport of neuronal vesicles (Wöllert et al., 2011). In L. vannamei, studies showed that Rab might be related to the resistance mechanism of shrimp induced by environmental stress (Wang et al., 2015). In addition, Rab gene also plays an important role in shrimp immunity, Rab like protein had been proved to interact with immune related membrane proteins to regulate the phagocytosis of shrimp hemolymph cells against WSSV, and Rab27 mutation leads to immune deficiency (Wu et al., 2008; Han et al., 2011; Chen et al., 2021). Similarly, Rab gene (Rab6a) in M. japonicus was up-regulated during WSSV infection (Wu and Zhang, 2007). In this study of WSSV-infected L. vannamei, several Rab family genes were significantly up-regulated in immune organs such as Oka organ, hepatopancreas and hemocyte. Therefore, Rab family may affect vesicle morphology, vesicle transport and immune regulation in shrimp.

In 1985, Madaule found a Ras superfamily member Rho (RAS homolog) in the marine gastropod molluscs Aplysia (Kawasaki et al., 2004). The study confirmed that Ras has highly sequence homology with Rho, with 35% consistency in amino acid sequence. Moreover, they have the same C-terminal required for membrane attachment. In L. vannamei, Rho family were found in the same clade of the phylogenetic tree with Ras family. These results suggest that Rho family is closely related to Ras family. In vertebrates, the Rho family has undergone considerable expansion and differentiated into more than 10 subfamilies, including Rho subtypes (A, B, C, D, G, E), Rac subtypes (1, 2, 3), Cdc42, Rnd (1, 2, 3), TCL, Rho H/TTF, Chp, Wrch-1, Rif, Rho BTB1, Rho BTB2, Miro-1 and Miro-2 (Burridge and Wennerberg, 2004; Wennerberg and Der, 2004). Five Rho subfamilies, including Rhol, Cdc42, Rac2, MIG-2, Rho BTB1, were identified in L. vannamei.

As a representative and well-studied member of Rho family, Cdc42 plays a significant role in a variety of cellular processes that are dependent on the actin cytoskeleton, such as cytokinesis, cell migration, phagocytosis, morphogenesis, axon myelination, intracellular trafficking, and tumor occurrence (Etienne-Manneville and Hall, 2002; Sahai and Marshall, 2002; Moon and Zheng, 2003). Rac is another member of the Rho family with more research. In mammals, Rac is mainly involved in promoting the malignant proliferation and migration of tumor cells (Chan et al., 2007). Compared with normal striated muscle tissue, the expression of Rac1 and Cdc42 was significantly high (P < 0.05) in rhabdomyosarcoma (RMS) tissue (Li et al., 2021). In addition, Rac also plays a role in immunity, the expression of Rac2 gene in the large yellow croaker Pseudosciaena crocea was significantly up-regulated after challenge by Vibrio parahaemolyticus (Liu et al., 2017). Injection of Vibrio alginolyticus into L. vannamei induced up-regulated expression of LvRac1 (LvRho1) in hepatopancreas, then after knocking down LvRac1 and stimulating V. alginolyticus, the mortality of L. vannamei was significantly increased relative to that of the control group (Cha et al., 2015). Rho1 is another important Rho GTPase protein, and previous research has shown a role for this in immunity, cytoskeleton dynamics and embryonic development. For example, Rho1 affected the development of eggs by regulating the hormone level of the braconid wasp Microplitis mediator (Magie and Parkhurst, 2005). Rho1 also regulated the rearrangement of cytoskeleton, and then affected the cellular immunity of the cotton bollworm Helicoverpa armigera (Li et al., 2010). In the purple sea urchin Strongylocentrotus purpuratus, Rho could affect SpROCK expression by the Rho dependent signal pathway, which is essential for early embryonic development (Aguirre-Armenta et al., 2011). In this study, we suggest Rho family may play an important role in development, cytoskeleton dynamics and immunity in L. vannamei.

Although the Ran is small family in number, it plays an important role in many species. In H. armigera, it was reported that Ran participated in the 20-hydroxyecdysone (20E) signal transduction pathway by regulating the location of ecdysone receptor-B1 (EcR-B1) (He et al., 2010). In the brown planthopper Nilaparvata lugens, NlRan knockdown significantly delayed development and affected reproduction (Liao et al., 2019). Only one Ran gene exists in L. vannamei, and its expression is higher than that of other Ras superfamily genes. It was previously reported that the highest expression of Ran is the black tiger shrimp, Penaeus monodon is in the ovary (Zhou et al., 2012). In the expression profile of testis maturation stages of scallop, Ran expression increased dramatically during meiosis and spermatogenesis (Hino et al., 2012). In addition, there is a similar situation in mammals. Ran had a high level of expression from the late pachytene spermatocytes to early round spermatocytes in mice (López-Casas et al., 2003). Moreover, the cellular localization of Ran also changed during spermatogenesis (Kierszenbaum et al., 2002). On the other hand, Ran may also involve in the immune process. In the Kuruma shrimp M. japonicus, Ran played a vital role in antiviral immunity (Han and Zhang, 2007). Another study found that Ran interacts with myosin in M. japonicus, which can regulate blood cell phagocytosis. RNAi knockdown led to a significant increase in virus copy number in M. japonicus, and overexpression of Ran resulted in a significant decrease in virus copy number (Liu et al., 2009). When IL-4 and lysophosphatidylcholine were respectively injected into shrimp, the results indicated that the two molecules could enhance the Ran GTPase activity and improve hemocytic phagocytosis against WSSV (Zhao et al., 2011). In our study, LvRan was significantly up-regulated in hemocytes and hepatopancreas, but down-regulated in Oka after WSSV infection. The above research suggests that Ran may participate in the 20E signaling pathway, regulate the gonad development and antiviral immunity, and then affect the growth and reproduction in shrimp.

The Arf family is mainly divided into three subfamilies: Arf, Arf like proteins (ARLs) and Sar1. Relevant studies have confirmed that Arf and the components that promote Arf function played an important role in mediating the transport of endoplasmic reticulum to the Golgi (Balch et al., 1992; Dong et al., 2010). For example, Arf1 regulated vesicle formation, Golgi assembly and promoting vesicle division, and Arf6 promoted membrane invagination on the cell surface during endocytosis (D’Souza-Schorey and Chavrier, 2006). Arf and Rab have similar functions in vesicle formation and transportation, indicating that the connection between GTPase-mediated signaling pathways requires different Ras superfamily proteins to fulfill a common task by a cooperation (Mitin et al., 2005). In D. melanogaster, Arf6 was only involved in spermatogenesis (Lambaerts et al., 2009), it was not required for early development in mice (Doherty and McMahon, 2009). In L. vannamei, Arf4 was almost not expressed in the early development stages, however, it had a higher expression in adults. There is almost no research on Arf in crustaceans, we found most Arf genes highly expressed in brain and ventral nerve of L. vannamei, so it is speculated that Arf should have a certain effect on the nervous system in shrimp.

The researches of the five classical Ras GTPase families have been very in-depth. Furthermore, many genes have regions predicted as or similar to the small GTPase domain, but they’re unconventional, for example, RJL family, RGK family, GPN family, IFT family, OBG family and Septin family. They have the same motif as classical Ras superfamily members, and the prediction region acts as a signal converter or molecular switch, called GTP-binding protein and has GTPase activity. Unlike the other members of the classical Ras GTPase, most these non-classical members showed extremely low levels of expression in adult tissues of L. vannamei, but most of them showed high levels of expression during early development and WSSV infection, this indicates that their functions are relatively specific and may be related to early development and immunity.