Roles of Crustacean Female Sex Hormone 1a in a Protandric Simultaneous Hermaphrodite Shrimp

Crustacean female sex hormone (CFSH) plays a pivotal role in the development of secondary sex characteristics in dioecious crustaceans. However, until now the knowledge concerning its functions in hermaphroditic species is scanty. Herein, we explored the function of CFSH (Lvit-CFSH1a) in the peppermint shrimp Lysmata vittata, a species characterized by a rare reproductive system of protandric simultaneous hermaphroditism (PSH). Lvit-CFSH1a cDNA was 1,220-bp in length with a 720-bp ORF encoded a polypeptide of 239-aa. RT-PCR showed that Lvit-CFSH1a was exclusively expressed in the eyestalk ganglion. For female physiology, it was found that Lvit-CFSH1a was indispensable for the development of female gonopores, but it might not involve vitellogenesis of the species. For male physiology, Lvit-CFSH1a suppressed Lvit-IAG2 expression in short-term silencing experiment and recombinant protein injection experiment, but did not affect male sexual differentiation in long-term silencing experiment. In addition, silencing the Lvit-CFSH1a gene impeded individual growth in L. vittata.


INTRODUCTION
Gonochorism is the most common reproductive strategy in decapod crustacean species (Juchault, 1999). In dioecious decapod crustaceans, androgenic gland hormone (AGH) and insulin-like androgenic gland hormone (IAG) secreted mainly by androgenic gland (AG), play critical roles in sexual differentiation (Manor et al., 2007). Implantation of AG into females induced masculinization in the red swamp crayfish Procambarus clarkii (Taketomi and Nishikawa, 1996), the giant freshwater prawn Macrobrachium rosenbergii (Nagamine et al., 1980;Malecha et al., 1992) and the red claw crayfish Cherax quadricarinatus (Khalaila et al., 2001;Barki et al., 2003), and vice versa (Barki et al., 2006). Likewise, silencing of IAG genes in males blocked male sexual differentiation while stimulating feminization. For instance, knockdown of IAG led to the arrest of testicular spermatogenesis and incomplete development of secondary sexual characteristics (appendices masculinae) in young M. rosenbergii males (Ventura et al., 2009). Moreover, silencing IAG gene could also feminize male-related phenotypes in the intersex C. quadricarinatus (Rosen et al., 2010) and the male Chinese mitten crab Eriocheir sinensis (Fu et al., 2020), and even induce fully sex reversal in young M. rosenbergii males (Ventura et al., 2009). Therefore, by virtue of its universal role as a master regulator of crustacean male development, IAG was also termed the sexual "IAG-switch" (Levy and Sagi, 2020). It was presumed that females arose as the absence of AG or IAG (Ventura et al., 2009).
Interestingly, a review of the literature shows that Caridean shrimps exhibit several other protandry sexual systems apart from gonochorism (Bauer, 2000). To date, however, there are rare reports elaborating sexual differentiation mechanism in protandric crustaceans. In the strictly sequential protandric hermaphroditism (SPH) shrimp Pandalus platyceros, Pnp-IAG knockdown elevated the expression of vitellogenin in the hepatopancreas and promoted transformation of the gonad from ovotestis to ovary . In the protandric simultaneous hermaphroditism (PSH) shrimp Lysmata wurdemanni, it was suggested that IAG was possibly responsible for the maintenance of the male reproductive activity in euhermaphrodite phase (Zhang et al., 2017). More detailed studies were performed in another PSH shrimp L. vittata. It was demonstrated that Lvit-IAG1 and Lvit-IAG2 jointly regulated male sexual differentiation (Liu et al., 2021a,b). Meanwhile, Lvit-IAG1 was also suggested to regulate the ovarian development by inhibiting Lvit-GIHs expression (Liu et al., 2021b). Recently, two CFSH transcripts (Lvit-CFSH1a and Lvit-CFSH1b) have been identified from L. vittata by transcriptomic analysis (Bao et al., 2020), however, no further in-depth studies are conducted.
In this study, we explored the function of Lvit-CFSH1a in L. vittata, a species that displays the unique PSH sexual system, whereby individuals first mature as males; and then, with increasing age and size, acquire reproductive functions of both males and females (Bauer, 2000;Alves et al., 2019). Considering that CFSH acts as an inhibitor of IAG (Liu et al., 2018), we hypothesized that Lvit-CFSH1a not only regulated female sexual differentiation, but also involved male sexual differentiation via inhibiting Lvit-IAGs expression of the PSH species. To validate this hypothesis, we performed both short-term and long-term gene knockdown via RNA interference (RNAi). In addition, we expressed recombinant Lvit-CFSH1a mature peptide (rLvit-CFSH1a) using a prokaryotic expression system, and carried out in vivo experiment to examine the effect of Lvit-CFSH1a on the expression of Lvit-IAG1, Lvit-IAG2, Lvit-Vg, and Lvit-VgR.

Animals
The experimental animals (L. vittata) were artificial-bred at the Fisheries Research Institute of Fujian Province, in Xiamen, China. After transport to the laboratory, shrimps were acclimated in seawater aquaria at temperature of 26 ± 1 • C and salinity of 32 ± 1 PSU for 2 days. During that period, they were fed with a commercially formulated shrimp diet daily. Two developmental phases covering four gonadal development stages were defined and described for L. vittata according to the previous research (Chen et al., 2019).

cDNA Cloning of Lvit-CFSH1a
Total RNA was extracted from the eyestalk ganglion of shrimps at gonadal development stage III using TRIzol R reagent (Invitrogen) according to the manufacturer's instructions. The first-strand cDNA for fragment and 3 -untranslated region (3 UTR) cloning was generated with 1 µg total RNA using RevertAid First Strand cDNA Synthesis Kit (Fermentas). Fragment of Lvit-CFSH1a was obtained from a de novo transcriptomic library of L. vittata and polymerase chain reaction (PCR) by primer pair (CFSH1aF/CFSH1aR) was performed to verify its accuracy. Seminested PCR was performed for 3 UTR cloning. The 5 UTR of Lvit-CFSH1a was obtained by a method of rapid amplification of cDNA ends (RACE) with SMART TM RACE cDNA Amplification Kit (Clontech) according to the manufacturer's protocol. Primers used in cDNA cloning were list in Table 1.

Genomic DNA Amplification of Lvit-CFSH1a
EasyPure R Marine Animal Genomic DNA Kit (TransGen) was used to extract the genomic DNA from the eyestalk ganglion of L. vittata at gonadal development stage III. Specific primers were designed for amplification of Lvit-CFSHs genomic DNA (gDNA) ( Table 1). The PCR reaction was performed with LA-Taq polymerase (TaKaRa) under the following conditions: 95 • C for 5 min; 35 cycles of 95 • C for 30 s, 60 • C for 30 s, and 72 • C for 2 min, followed by 72 • C for 10 min final extension.

The Quantitative Real-Time PCR Assays
Primers used for quantitative real-time PCR (qRT-PCR) were from previous studies (Liu et al., 2021a,b). Amplification efficiency of each primer pair was determined before used for qRT-PCR assays. The cDNA was diluted fourfolds using RNase-free water before it was utilized in qRT-PCR detection. Components including, 10 µl TB Green Premix Ex Taq II (2X) (TaKaRa), 2 µl diluted cDNA, 0.5 µl forward/reverse primer (10 µM), as well as and 7 µl RNase-free water, were used for a 20 µl qRT-PCR reaction system. The reaction was performed using in 7500 Real-Time PCR (Applied Biosystems): 95 • C for 30 s, followed by 40 cycles of 95 • C for 15 s, 58.5 • C for 15 s and 72 • C for 30 s. The result was calculated using 2 − Ct method, whereby Lvit-βactin (GenBank accession number: MT114194) was utilized as the reference gene. Tissue Expression of Lvit-CFSH1a in L. vittata The total RNA was extracted from various tissues (eyestalk ganglion, brain, thoracic ganglion, abdominal ganglion, ovary, testis, androgenic gland, hepatopancreas, stomach, intestine, heart, gill, and muscle) as described in section "cDNA Cloning of Lvit-CFSH1a." The first-strand cDNA was produced from 1 µg total RNA using PrimeScript TM RT reagent Kit with gDNA Eraser (Perfect Real Time) (TaKaRa). Tissue expression profile was determined by RT-PCR under the following conditions: 95 • C for 5 min; 35 cycles of 95 • C for 30 s, 58.5 • C for 30 s, and 72 • C for 30 s, followed by 72 • C for 5 min final extension. Meanwhile, Lvitβ-actin (GenBank accession number: MT114194) was amplified as a positive control. RT-PCR products were examined using 1.5% agarose gel imaged by UV detector (Geldoc, Thermo Fisher Scientific).

Expression Profile of Lvit-CFSH1a During Gonadal Development
In order to examine expression profile of Lvit-CFSH1a during gonadal development, total RNA was extracted from eyestalk ganglion of L. vittata at different gonadal development stages (n = 5). It was followed by the synthesis of cDNA as described in section "Tissue Expression of Lvit-CFSH1a in L. vittata" and qRT-PCR detection as described in section "The Quantitative Real-Time PCR Assays, " respectively.

Short-Term Silencing Experiment in vivo
To evaluate the efficacy of gene knockdown via RNA interference (RNAi), a short-term silencing experiment was carried out with L. vittata at gonadal development stage I. A total of 15 shrimps (carapace length 2.99 ± 0.10 mm, body weight 42.79 ± 4.28 mg) were randomly and equally assigned to following 3 treatment groups (n = 5): dsRNA Lvit-CFSH1a-injected, dsRNA GFPinjected and PBS-injected. dsRNA (2 µg/g) (Liu et al., 2021a,b) was delivered via intramuscular injection in the abdominal segment of shrimp. Meanwhile, the PBS-injected group received an equivalent volume of PBS. Sampling was performed 24 h after injection (Shi et al., 2020). Following anesthesia on ice for 5 min, eyestalk ganglion (EG) and androgenic gland (AG) were collected to test the effect of dsRNA-mediated silencing on the specific genes by qRT-PCR. While RNA extraction and qRT-PCR were performed as described above in sections " cDNA Cloning of Lvit-CFSH1a" and "The Quantitative Real-Time PCR Assays, " the first-strand cDNA was generated with 200 ng total RNA by TransScript R II One-Step gDNA Removal and cDNA short SuperMix Kit (TransGen).

Long-Term Silencing Experiment in vivo
A long-term silencing experiment was conducted to determine the potential roles played by Lvit-CFSH1a in sexual differentiation and gonadal development. Shrimps (carapace length 3.12 ± 0.17 mm, body weight 47.80 ± 6.42 mg) at gonadal development stage I were randomly divided into three groups (n = 13) as described in section "Short-Term Silencing Experiment in vivo." Similar dose of dsRNA (2 µg/g) or equal volume of PBS was injected into the abdominal segment of shrimp once every 4 days for a total of 8 injections in 29-days duration during which shrimps were kept in seawater aquaria under the following conditions: temperature, 26 ± 1 • C; salinity, 32 ± 1 PSU; 12L:12D photoperiod. On day 30 (24 h after the 8th injection), all of the shrimp were sampled after anesthetization on ice (Shi et al., 2020). Measurements of carapace length and body weight were recorded. Changes in external sexual features and gonadal development were assessed and recorded as described in our laboratory (Liu et al., 2021a,b). Samples of eyestalk ganglion, androgenic gland, ovarian region of the gonad and hepatopancreas were collected to test the effect of long-term silencing on Lvit-CFSH1b (GenBank accession number: MT114198), Lvit-IAG1 (GenBank accession number: MT114196), Lvit-IAG2 (GenBank accession number: MT114197), Lvit-Vg (GenBank accession number: MT113122) and Lvit-VgR (GenBank accession number: MT114195) by qRT-PCR. RNA extraction, the first-strand cDNA synthesis, qRT-PCR were also performed as described in section "Short-Term Silencing Experiment in vivo." The remaining gonad tissue was fixed in modified Bouin's fixative (25 ml 37-40% formaldehyde, 75 ml saturation picric acid and 5 ml glacial acetic acid) at 4 • C for 24 h, followed by gradient alcohol dehydration, paraffin embedding, preparation of 6 µm sections and staining with hematoxylin and eosin (H & E) for histological observation.
In vivo Effect of rLvit-CFSH1a on Gene Expression rLvit-CFSH1a was expressed using a prokaryotic expression system and purified by immobilized metal-affinity chromatography (IMAC) (Bornhorst and Falke, 2000). The fragment encoding the mature peptide of Lvit-CFSH1a was cloned into pET-His vector with restriction enzyme sites (EcoRI and NheI). The generated constructs (pET-His-CFSH1a) was transformed into E. coli TransB (DE3) and induced at 16 • C after adding isopropyl-beta-D-thiogalactopyranoside (IPTG, 0.2 mM final concentration). After 20 h, bacterial cells were harvested by centrifugation. Because rLvit-CFSH1a expressed as inclusion bodies (Supplementary Figure 1A), purification was performed under denaturing conditions (8M urea) with Ni Sepharose TM 6 Fast Flow (GE Healthcare) according to the manufacturer's instructions (Supplementary Figure 1B). Purified rLvit-CFSH1a was renatured by graded urea dialysis and confirmed by Western blot analysis. Samples were separated by 12% SDS-PAGE gel electrophoresis. The electrophoresed proteins were transferred to a PVDF membrane and blocked with 5% bull serum albumin (BSA)-PBS for 1 h at room temperature. Following blocking, the membrane was washed three times with PBS. After that, it was incubated in ProteinFind R Anti-His Mouse Monoclonal Antibody (1:2,000, TransGen) for 1 h at 37 • C. Horseradish peroxidase activity was detected with Western Blotting Mouse IgG DAB Chromogenic Reagent Kit (incubated with Goat Anti-Mouse IgG/HRP, 1:4,000, 1 h, 37 • C; Solarbio) according to the manufacturer's instructions. A total of 10 shrimps (carapace length 5.17 ± 0.22 mm, body weight 179.02 ± 29.41 mg) at gonadal development stage II were randomly and equally assigned to following 2 treatment groups (n = 5): rLvit-CFSH1a-injected and PBS-injected. rLvit-CFSH1a (2 µg/g) was delivered via intramuscular injection in the abdominal segment of shrimp. Meanwhile, the PBS-injected group received an equivalent volume of PBS. Sampling was performed 24 h after injection. After shrimps were anesthetized on ice for 5 min, AG, ovarian region and hepatopancreas were obtained to test the effect of rLvit-CFSH1a-injection on the specific genes by qRT-PCR. RNA extraction, the first-strand cDNA synthesis, qRT-PCR were also performed as described in section "Short-Term Silencing Experiment in vivo."

Bioinformatics Analyses
The primers used for cDNA cloning, gDNA cloning, dsRNA preparation and prokaryotic expression were designed with Primer 5.0. The open reading frame (ORF) was predicted by ORF Finder software 1 . We adopted the SignalP-5.0 Server 2 to predict the signal peptides, whereas cysteine residues and putative disulfide bonds were predicted via DiANNA 1.1 web server 3 . Further, sequence alignment of deduced amino acid sequences with reported sequences was performed using the Clustal Omega website 4 . N-glycosylation motif was predicted by NetNGlyc 1.0 Server 5 .
The Maximum Likelihood method with 1000 bootstrap replicates based on the JTT matrix-based model in MEGA7 was applied to generate a phylogenetic tree entailing the deduced amino acid sequence of decapoda CFSH mature peptides. Most of CFSH sequences were borrowed from previous works (Thongbuakaew et al., 2019); other sequences were shown in     Table 2. The eight conserved cysteine residues were boxed in red and other cysteine residues were shown in green. The light blue bar showed conserved IL-17domain. The N-glycosylation motif was shown in gray background. (B) Phylogenetic tree of CFSHs in decapod crustaceans. Sequences used in the alignment were the same as those used in multiple sequence alignment. Phylogenetic analysis was conducted by Maximum Likelihood method based on the JTT matrix-based model in MEGA7.
The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) was shown next to the branches. Lvit-CFSH1a was indicated with blue box.

Statistical Analyses
Normality of data was established by the Kolmogorov-Smirnov test. All the data were presented in a normal distribution and tested for variances homogeneity by the Levene's test. All statistical analyses were performed using the SPSS 18.0 software. Statistical significance (p < 0.05) of data regarding Lvit-CFSH1a expression profile and silencing experiments was determined using one-way ANOVA followed by Tukey's multiple range tests. The t-test was used to analyze the data of rLvit-CFSH1a injection experiment. All data were presented as mean ± SD.

Homology and Phylogenetic Analysis
Multiple sequence alignment of deduced mature peptide of CFSHs was shown in Figure 2A. Eight conserved cysteine residues and an interleukin 17 (IL-17) domain were found in Lvit-CFSH1a, and that was faultlessly aligned with other CFSHs (Figure 2A). According to the former studies Tsutsui et al., 2018), the CFSHs were grouped into two subtypes: type I and type II. The type I CFSHs possessed a single conserved N-glycosylation motif while the type II CFSHs contained additional one or two cysteine residues (Figure 2A). Lvit-CFSH1a mature peptide shared the highest identity with Can-CFSH1a (84.15%). Phylogenetic analysis demonstrated that the CFSHs in the decapod crustaceans formed two major clades: type I and type II CFSH. Lvit-CFSH1a was classified into type I CFSH ( Figure 2B).

Expression Profiles of Lvit-CFSH1a
RT-PCR was performed on L. vittata at the gonadal development stage III to determine the spatial distribution profile of Lvit-CFSH1a. The findings were that Lvit-CFSH1a was exclusively expressed in the eyestalk ganglion ( Figure 3A). The relative expression of Lvit-CFSH1a in the eyestalk ganglion during gonadal development was also assessed by qRT-PCR. No significant difference was observed in Lvit-CFSH1a expression levels during gonadal development (F 3 , 16 = 0.265, p = 0.850) (Figure 3B).

Effects of Lvit-CFSH1a Silencing on Growth and Development of Sexual Characteristics
At the end of the 30-day long-term trial, we recorded the average carapace length and bodyweight of the shrimps. Compared with individuals from the PBS (4.94 ± 0.06 mm, 140.36 ± 1.37 mg) or dsRNA GFP (4.82 ± 0.09 mm, 142.90 ± 2.05 mg) treatments, shrimps in dsRNA Lvit-CFSH1a (4.09 ± 0.05 mm, 101.90 ± 3.84 mg) treatment were significantly smaller [carapace length: F (2 , 15) = 5.441, p < 0.05; body weight: F (2 , 15) = 9.124, p < 0.05] (Figures 6A,B). Moreover, changes in female and male sexual characteristics were also documented through photography (Figure 7). On the one hand, Lvit-CFSH1a gene knockdown led to retardation of female sexual characteristics. In the PBS and dsRNA GFP treatment, gonophores bulged out like a frustum surrounded by lush feathery setae (Figures 7A,B). Female gonophores in dsRNA Lvit-CFSH1a treatment were hypogenesis ( Figure 7C). Contrary to the controls, gonophores were less visually evident or completely disappeared. The feathery setae surrounding the female gonopores located at the base of the third pair of pereiopods were substantially sparse ( Figure 7C). On the other hand, no significant difference in male characteristics (cincinnuli, AM and male gonopore) was observed (Figures 6C, 7D-L).

Effects of Lvit-CFSH1a Silencing on Gonadal Development
Morphological and histological characteristics of gonads following Lvit-CFSH1a knockdown were recorded. Knockdown of Lvit-CFSH1a led to retardation of ovarian region. In dsRNA Lvit-CFSH1a treatment, the ovarian region was thinner ( Figure 8C). Additionally, the oocyte size significantly decreased following injection of dsRNA Lvit-CFSH1a (Figures 6D, 8D-F). For the PBS and dsRNA GFP treatment, the average oocyte diameter was 55.57 ± 0.51 µm (n = 7) and 53.87 ± 0.28 µm (n = 6), respectively. Following knockdown of Lvit-CFSH1a, the oocyte diameter decreased to 48.81 ± 0.46 µm (n = 5). However, Lvit-CFSH1a knockdown had no significant effect on testicular development. Testicular regions of the three treatments were cloudy white (Figures 8A-C). Gonadal histology further showed that similar compositions of germ cell types were observed among the three treatments. Abundant of spermatocytes I (Sc I), spermatid (Sd), and spermatozoa (Sz) were found in the testicular region of the three treatments (Figures 8G-I).

DISCUSSION
CFSH is a crucial hormone involved in development of femalerelated phenotypes in dioecious crustaceans (Zmora and Chung, 2014;Jiang et al., 2020). However, the biological functions of CFSH in PSH crustaceans have been rarely investigated.
Our findings revealed that RNAi induced a specific knockdown of the Lvit-CFSH1a transcripts level by 83.7% in a short-term experiment, and by 81.7% in a long-term experiment. Furthermore, we assessed the influence of Lvit-CFSH1a on the development of female features by comparing The expression levels of Lvit-IAG1, Lvit-IAG2, Lvit-Vg, and Lvit-VgR were detected following in vivo injection with PBS or rLvit-CFSH1a. The gene expression levels were standardized by Lvit-β-actin expression levels and represented as mean ± SD (t-test, with **p < 0.01; n = 5).
female external features and ovarian development status. In the present study, we found that Lvit-CFSH1a gene knockdown led to retardation of female sexual characteristics. Numbers of feathery setae surrounding the female gonopores were substantially sparse. This result is similar to those described in C. sapidus and S. paramamosain (Zmora and Chung, 2014;Jiang et al., 2020). However, results varied when it came to the influences on ovarian development. In fact, to date the effect of CFSH on ovarian development remains obscure. In C. sapidus, knockdown of Cs-CFSH showed no significant effect on the ovarian development (Zmora and Chung, 2014). Study in M. japonicus suggested that CFSH might participate in some reproductive process other than vitellogenesis (Tsutsui et al., 2018). Research in M. rosenbergii also demonstrated that both type I and type II CFSHs had unknown effect on ovarian development (Thongbuakaew et al., 2019). In the present study, a repression phenomenon was observed in the ovarian region of L. vittata in long-term silencing experiment. Contrary to the controls, ovarian region in the Lvit-CFSH1a silencing treatment was less developed with significantly smaller oocytes. Moreover, qRT-PCR data lent support to histomorphological results. It is known that both the Vg gene expression levels in hepatopancreas and VgR gene expression levels in ovary are approved indexes of ovarian development (Warrier and Subramoniam, 2002;Subramoniam, 2011;Jia et al., 2013;Urtgam et al., 2015). Following Lvit-CFSH1a knockdown, Lvit-Vg expression in the hepatopancreas and Lvit-VgR expression in the ovary were found significantly down-regulated. It seemed plausible that Lvit-CFSH1a stimulated ovarian development via promoting vitellogenesis. We also purified Lvit-CFSH1a mature peptide and conducted a further in vivo injection experiment in the study. However, results of rLvit-CFSH1a administration demonstrated that rLvit-CFSH1a affect neither vitellogenin synthesis in the hepatopancreas nor vitellogenesis in the gonad. Moreover, no significant difference was observed in Lvit-CFSH1a expression levels during gonadal development (Figure 3). These results jointly suggested that Lvit-CFSH1a might regulate ovarian development via some unknown process other than vitellogenesis. It is worth noting that, as shown in Figures 6A,B, Lvit-CFSH1a silencing induced significantly slower growth in the treated shrimps. It is known that growth and reproduction are closely related, gonadal development and fecundity usually increase with body size (Heino and Kaitala, 2001;Michalakis et al., 2013). The same phenomenon has also been observed in L. vittata (Chen et al., 2019). Thus, another possibility is that long-term Lvit-CFSH1a silencing caused significantly slower growth in the treated shrimps and this indirectly hindered the ovarian development, along with smaller oocytes and down-regulation of Lvit-Vg and Lvit-VgR expression.
Simultaneously, we evaluated the effects of Lvit-CFSH1a on male development in the study. In S. paramamosain, it was reported that CFSH acts as an inhibitor of IAG (Liu et al., 2018). In dioecious decapod crustaceans, there is usually one IAG gene in a species (Li et al., 2012;Chung, 2014;Huang et al., 2017). Transcripts generated by alternative splicing of the same IAG gene have different functions in different organs (Li et al., 2012;Chung, 2014;Huang et al., 2017). Nevertheless, two IAG genes were identified in the PSH shrimp L. vittata, and they cooperatively modulated male sexual differentiation of the species (Liu et al., 2021a,b). Specifically, Lvit-IAG1 was suggested to be closely related to the development of both AM and male gonopores, and participated in primary-tosecondary spermatocyte transition (Liu et al., 2021b); whereas Lvit-IAG2 was somewhat related to the development of AM, and participated in secondary spermatocyte-to-spermatid transition (Liu et al., 2021a). In the present study, results of short-term silencing experiment and rLvit-CFSH1a injection experiment demonstrated that Lvit-CFSH1a acted as an inhibitor of Lvit-IAG2 rather than Lvit-IAG1. Thus, we speculated that knockdown of Lvit-CFSH1a might slightly stimulate male sexual differentiation in long-term silencing experiment, such as relatively longer AM and more spermatid/spermatozoa in the testicular regions. However, at the end of the 30-day longterm trial, no significant promotion in male sexual differentiation was observed. Long-term Lvit-CFSH1a knockdown affected neither testicular maturation nor development of male-related phenotypes. Notably, Lvit-IAG2 expression was significantly suppressed, which was contradictory to the above hypothesis. The following may explain these apparently contradictory results. Previous studies suggested that Lvit-IAG2 is more than a sexual differentiation regulator, and it also stimulates the growth of the PSH species (Liu et al., 2021a). In the present study, silencing the Lvit-CFSH1a gene impeded individual growth of the species. Moreover, no significant difference was observed in Lvit-CFSH1a expression levels during the life cycle of the PSH species (Figure 3). Based on the above results, we proposed that Lvit-CFSH1a might involve some unknown biological processes and ultimately influence individual growth. In the short-term trial (Figures 4, 9), Lvit-IAG2 expression was significantly suppressed by Lvit-CFSH1a, but in the long-term experiment the adverse effects of Lvit-CFSH1a silencing on growth became more obvious (Figure 6), and it in turn affected expression of growth-related gene (e.g., Lvit-IAG2). With the decrease of Lvit-IAG2, the rate of male differentiation in dsRNA Lvit-CFSH1a treatment was slower than the PBS or dsRNA GFP treatment, which eventually led to similar development of male-related phenotypes among the three treatments.
In summary, we characterized a CFSH gene, Lvit-CFSH1a, from the eyestalk ganglion of the PSH species L. vittata. This study showed that Lvit-CFSH1a regulated female-related phenotypes, but didn't evidently affect the male sexual differentiation. In addition, Lvit-CFSH1a might participate in the regulation of individual growth.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
FL contributed to conceptualization, methodology, software, validation, formal analysis, investigation, data curation, visualization, and writing-original draft preparation of the study. HY contributed to conceptualization, methodology, validation, data curation, writing-review and editing, supervision, project administration, and funding acquisition. WS and LH contributed to investigation. GW contributed to funding acquisition. ZZ contributed to funding acquisition and provided resources. All authors contributed to manuscript revision, read, and approved the submitted version.

FUNDING
This work was supported by the Special Fund of Marine and Fishery Structure Adjustment in Fujian (2020HYJG01 and 2020HYJG08).