Abstract Estradiol is an important sex steroid hormone that is involved in the regulation of crustacean ovarian development. However, the molecular regulatory mechanisms of estradiol on ovarian development are largely unknown. This study performed transcriptome sequencing of ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ of crabs after estradiol treatment (0.1μg g^-1 crab weight). A total of 23, 806 genes were annotated, and 316, 1300, 669, 142, 383 genes were expressed differently in ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ respectively. Differentially expressed gene enrichment analysis revealed several crucial pathways including protein digestion and absorption, pancreatic secretion, insect hormone biosynthesis, drug metabolism-cytochrome P450 and signal transduction pathway. Through this study, some key genes in correlation with the ovarian development and nutrition metabolism were significantly affected by estradiol, such as vitelline membrane outer layer 1-like protein, heat shock protein 70, Wnt5, JHE-like carboxylesterase 1, cytochrome P302a1, crustacean hyperglycemic hormone, neuropeptide F2, trypsin, carboxypeptidase B, pancreatic triacylglycerol lipase-like, and lipid storage droplet protein. Moreover, RT-qPCR validation demonstrated that expression of transcripts related to ovarian development (vitelline membrane outer layer 1-like protein and cytochrome P302a1) and nutrition metabolism (trypsin, glucose dehydrogenase and lipid storage droplet protein) were significantly affected by estradiol treatment. This study not only has identified relevant genes and several pathways that are involved in estradiol regulation on ovarian development of P. trituberculatus, but also provided new insight into the understanding of the molecular function mechanisms of estradiol in crustacean. Introduction In many decapod crustacean, regulation of reproductive maturation is a major process in commercial aquaculture. Reproduction in female crustaceans is characterized by ovarian development, which includes the two processes of oogenesis and vitellogenesis [[34]1, [35]2]. Crustacean ovarian development has been proved to be regulated by various hormonal factors [[36]3]. Estradiol, an important sex steroid, is an important active estrogen in most crustacean, which promotes the process of vitellogenesis and ovarian development [[37]4–[38]6]. However, the negative effects of estradiol on ovarian development and vitellogenesis were not detected on some crustaceans with a particular ovarian stage, including the immature tiger prawn, Penaeus esculentus and ridgeback shrimp Sicyonia ingentis, during the non-reproductive stage [[39]7, [40]8]. The possible reason for the contradictory results was that the effects of estradiol on vitellogenesis are ovarian stage-specific and dose-dependent for crustaceans [[41]5, [42]9, [43]10]. Furthermore, previous studies have shown that estradiol is distributed widely in several tissues such as hepatopancreas, ovary and hemolymph, and the concentration of estradiol shows a significant positive correlation with the ovarian development [[44]11–[45]13]. However, the molecular regulatory mechanisms of estradiol on ovarian development are largely unknown in crustacean. Transcriptome sequencing enables the production of high-throughput fragments of double-stranded cDNA and the rapid assembly of sequences for annotation. It facilitates gene discovery and broadens our understanding of gene networks, especially in non-model organisms with unknown genomes [[46]14]. Over the past few decades, numerous studies have shown that the endocrine regulation of crustacean ovarian development is complex with multiple hormonal factors employed to positively control ovarian development [[47]3, [48]15, [49]16]. For example, X-organ-sinus gland complex system in the eyestalk secretes many crucial neuropeptide hormones, including crustacean hyperglycemic hormone and molt inhibiting hormone that exhibit a positive influence on vitellogenin activities [[50]17–[51]19]. Brain ganglion secretes biogenic amines and gonad stimulatory hormone could also promote the vitellogenesis of crustacean, and biogenic amines can enhance the release of gonad stimulatory hormone [[52]20]. Similarly, methyl farnesoate and ecdysteroids secreted by mandibular organ and Y-organ, respectively, also exhibited gonad stimulatory function [[53]3, [54]15, [55]21, [56]22]. Remarkably, vertebrate-type sex steroids, especially estrogen and progesterone, play a vital role in crustacean vitellogenesis by stimulating related metabolic pathways initiation during vitellogenesis, such as lipogenesis [[57]15]. Those hormonal factors originating from different endocrine organs are involved in ovarian development of crustaceans individually or in synergy with one another, but data on this aspect are still inadequate [[58]3]. The swimming crab, Portunus trituberculatus, is an important marine-culture crab widely distributed in the coastal water area of East Asia, including Korea, Japan, Philippines and China [[59]23, [60]24]. Previous studies in our group have shown that estradiol widely distributed in various endocrine organs of P. trituberculatus and exogenous estradiol could promote the vitellogenesis and ovarian development [[61]10, [62]25, [63]26]. Therefore, there are two objectives in this study: (1) to obtain ten reference transcriptome assemblies for the ovary, hepatopancreas, brain ganglion, eyestalk and mandibular organ from control crabs and those exposed to an acute level of estradiol; (2) to identify estradiol-responsive gene profile differences to unravel the different mechanisms of action of estradiol in five organs. The results will provide new insight into the understanding of the molecular regulatory mechanisms of estradiol in P. trituberculatus and other homologous species. Materials and methods Experimental animals, design, and sampling Female P. trituberculatus that have just finished puberty molt (body weight = 150 ± 25 g) were collected from outdoor ponds of Qidong scientific research base, Shanghai Fisheries Research Institute, Jiangsu, China, and then acclimated for a week in the indoor circulating water system. During the experiment, these crabs were individually cultured in the culture box (L × W × D = 33 × 27.5 × 35 cm), all culture boxes were floating in two concrete tanks (Length × Width × Depth = 5.8 m × 2.4 m × 1.8 m) by foam attached to the surrounding of the boxes. There were many uniform holes on the walls of each box that allow ample water exchange with the concrete tanks. The water depth of each basket was maintained at ca. 25 cm with a layer of 5–6 cm sand provided to the bottom of each boxes for the crabs to bury [[64]27]. Prior to beginning the experiment, 20 crabs were randomly divided into control (n = 10) and treatment groups (n = 10). The crabs assigned to the treatment group received estradiol solution (Sigma-Aldrich, dissolved in absolute ethanol) injection through the arthrodial membrane at the base of the swimming-leg. The injection dose (0.1μg g^-1 crab weight) and volume (0.5μl g^-1 crab weight) were followed the previous publications on the other decapod crustacean [[65]4, [66]7, [67]28]. The control crabs received the same volume of absolute ethanol [[68]4, [69]7, [70]28]. 24 h post injection, eight crabs from each group were randomly sampled for RNA extraction and RNA-seq analysis. All samples, including the ovary (O), hepatopancreas (H), brain ganglion (BG), eyestalk (SG) and mandibular organ (MO), were frozen with liquid nitrogen and then stored at -80°C until analysis. All crabs were treated in strict accordance with the guidelines for the care and use of experimental animals established by the Administration of Affairs Concerning Experimental Animals of the State Council of the People’s Republic of China, and approved by the Committee on Experimental Animal Management of the Shanghai Ocean University. RNA extraction, cDNA library preparation and next-generation sequencing Total RNA was extracted from the tissues of P. trituberculatus using Trizol Reagent (Takara) according to the manufacturer's recommendations and genomic DNA was removed with DNase I (Takara). Subsequently, the integrity and purity were estimated by a 2100 Bioanalyzer (Agilent Technologies, USA) and NanoDrop 2000 (Thermo Fisher Scientific Inc., USA), respectively. Only high-quality RNA samples (OD260/280 ranged 1.8–2.2, RIN ≥8.0) were used to construct the sequencing library. The cDNA libraries were generated using the TruSeq RNA sample prep kit (Illumina) following manufacturer's recommendations. The synthesized cDNA libraries were checked using PicoGreen (Quantifluor^™-ST fluorometerE6090, Promega, CA, USA) and fluorospectrophotometry (Quant-iT PicoGreen dsDNA Assay Kit; Invitrogen, P7589) and quantified with Agilent 2100 (Agilent 2100 Bioanalyzer, Agilent, 2100; Agilent High Sensitivity DNA Kit, Agilent, 5067–4626). The final sequencing cDNA libraries were quantified to 4–5 pM and sequenced using the Illumina HiSeq 4000 with 150 bp pair-end reads produced (Illumina, USA). De novo assembly and annotation functional of the transcriptome Before proceeding to de novo assembly, the raw reads were first quality-filtered to remove adaptor sequences, ambiguous ‘N’ nucleotides (with an ‘N’ ratio over 10%), and low-quality sequences (with quality scores lower than 20). The clean reads were assembled into non-redundant transcripts using the Trinity software, which has been developed specifically for the de novo assembly of transcriptomes using short reads [[71]29]. The assembled unigenes were then used for sequence annotation using the databases including the NCBI non-redundant protein database (NR), Gene Ontology (GO), Evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG), Swiss-Prot protein, and Kyoto Encyclopedia of Genes and Genome (KEGG) [[72]30–[73]32]. Identification of differentially expressed genes (DEGs) The expression of all unigenes was estimated by calculating read density as ‘reads per kb per million reads’ (RPKM) using the RSEM program [[74]33]. To identify the DEGs, |log2(FoldChange)| > 1 & P-value < 0.05 were set to be the threshold for significantly different expression levels. The gene expression profiles were compared with estradiol and control group, and then all DEGs in each comparison were carried on the GO functional and KEGG pathway enrichment analysis using GO database and KEGG database. Validation of the RNA-seq profiles by quantitative real-time PCR (qPCR) Randomly selected DEGs identified by transcriptome sequencing analysis were used for RNA-seq validation by qRT-PCR. RNA from 24h samples (8 replicates from each treatment) was used for cDNA synthesis using a reverse first strand cDNA synthesis kit (RR036A, Takara Bio, Japan). The qRT-PCR assays were performed with three replicates, and the 18S gene was used as an internal control to normalize the expression level of the target genes [[75]34]. The specific primers were designed according to the unigene sequences using Primer 5.0 software ([76]Table 1). Table 1. Specific primers used to in qRT-PCR. Primer name Sequence (5′→3′) Vitelline membrane outer layer 1-like protein F: GTTCGAGCCCTACACCTTCTTG R: ACGATCTCATGTCCTGCCACTC Glucose dehydrogenase F: TATGATGCCGCATGGATTTCTG R: TCCTTGTGGACGGGTAGATGAC Lipid storage droplet protein F: GCCGTTCCGATAATCAACCTAC R: AGAGTTTGGGGCACGATCAAGT CYP302a1 F: CCAACTGTCTTGTTCTCCCCACT R: CTCTTCTGCCGAGCATGTCTCA Trypsin F: GACAGCACGGCACTTACAAACG R: TACGATGCAATCAACGCCTACA 18S F: TCCAGTTCGCAGCTTCTTCTT R: AACATCTAAGGGCATCACAGACC [77]Open in a new tab The volume of the qRT-PCR reaction system was 10 μl, containing 5 μL of 2×SYBR Master Mix (RR420A, Takara Bio, Japan), 0.2 μL of each primer and ROX Reference Dye II (RR420A, Takara Bio, Japan), 1 μL of cDNA template, and 3.4 μl of RNase Free dH[2]O. The qPCR was carried out in a FAST-7500 system (ABI-7500, ThermoFisher, Singapore) as follows: 95°C for 30 s; 40 cycles of 95°C for 5 s, 60°C for 30 s, and 72°C for 30 s. The relative quantification of qRT-PCR data was calculated using the 2^-ΔΔCt method [[78]35]. Data were presented as the mean ± standard error (SE). All statistical analyses were performed using SPSS 19.0 software (SPSS, Chicago, USA). Results Sequencing and de novo assembly A total of 441,394,570 Illumina HiSeq reads from ovary, hepatopancreas, brain ganglion, eyestalk and mandibular organ were generated. Each sample yielded more than 40 million reads of data. All raw reads and quality control statistics are presented in [79]Table 2. Table 2. Raw reads and quality control of reads for P. trituberculatus libraries. Tissue Treatment Number of raw reads GC (%) Q30 (%) Number of clean reads Reads filtered (%) Ovary estradiol (O-1) 45,496,138 56.87 93.84 44,982,290 98.87 ethanol (O-2) 47,538,614 56.64 93.29 46,931,094 98.72 Hepatopancreas estradiol (H-1) 42,206,722 57.64 93.59 41,865,022 99.19 ethanol (H-2) 43,020,694 55.87 93.95 42,674,458 99.19 Brain ganglion estradiol (BG-1) 43,396,114 53.85 93.3 43,007.618 99.1 ethanol (BG-2) 47,549,836 54.45 93.42 47.161,352 99.18 Eyestalk estradiol (SG-1) 44,105,484 55.89 93.52 43,742,076 99.17 ethanol (SG-2) 45,100,746 56.2 93.81 44,768,352 99.26 Mandibular organ estradiol (MO-1) 42,104,066 56.65 92.81 41,865,022 99.06 ethanol (MO-2) 40,876,156 57.69 93.26 40,621,404 99.37 [80]Open in a new tab De novo sequence assembly of the current dataset using the Trinity assembler generated 106,130 transcripts, with an average transcript size of 681 base pairs (bp) and N50 of 1070 bp. Finally, 84, 032 unigenes were obtained after combining the transcripts, with a total length of 49,104,718 bp. The mean unigene length was 584.36 bp, and the N50 was 789 bp. A brief summary of de novo assembly statistics was provided in [81]Table 3. The raw data were uploaded to the National Center for Biotechnology Information (NCBI, PRJNA532740). Table 3. De novo assembly statistics of the P. trituberculatus transcriptome. Transcripts Unigenes Total number 106,130 84,032 Total length 72,377,107 49,104,718 Max. length (bp) 14,331 14,331 Mean length (bp) 681.97 584.36 N50 length (bp) 1,070 789 GC% 50.32 59.36 [82]Open in a new tab Functional annotation Functional annotation of the assembly identified 84,033 unique protein-encoding genes, 23,806 (28.33%) of which could be matched to the NR database ([83]Table 4), while 9, 650 (11.48%) transcripts could be fully annotated with GO available data, illustrating the scarcity of crustacean sequences in current genomic databases. Species most represented in the BLASTx searches included Zootermopsis nevadensis (7%), Limulus polyphemus (5%) and Daphnia magna (4%), but “other” species was the largest group (74%) ([84]Fig 1). Table 4. Summary of functional annotation of P. trituberculatus transcriptome. Annotation in Database Unigene No. Percentage (%) NR 23,806 28.33 GO 9,650 11.48 KEGG 5,814 6.92 eggNOG 22,625 26.92 Swissprot 21,439 25.51 [85]Open in a new tab Fig 1. Species distribution of BLASTx hits. [86]Fig 1 [87]Open in a new tab The GO assignment programs were also used for functional categorizations of the annotated unigenes. These unigenes were classified into three groups including biological processes, cellular components, and molecular function. These unigenes were then further divided into 67 functional subgroups. Among biological processes category, metabolic process (4,079) and cellular process (4,950) were the most enriched GO terms. Within the cellular components, transcripts assigned to cell (3,762), cell part (3,734), and organelle (2,622) were the three most abundant. In molecular function category, binding-related genes (4,498) and catalytic activity (4,303) were the most enriched GO terms ([88]Fig 2). Fig 2. Gene Ontology (GO) categorization for assembled unigenes. [89]Fig 2 [90]Open in a new tab Assignments of eggNOG were further used to evaluate the completeness of the transcriptome library and the efficiency of the annotation process. These sequences were categorized into 25 categories. Of these, the “general function prediction only” (5,303) was the largest cluster and “Cell motility” (12) was the smallest one ([91]Fig 3). Fig 3. Clusters of eggNOG functional classification of P. trituberculatus transcriptome. [92]Fig 3 [93]Open in a new tab To systematically analyse the associated intracellular and intercellular metabolic pathways as well as complicated biological behaviours, the unigenes were annotated by the KEGG database. A total of 5,814 unigenes were assigned to five pathway categories, including metabolism, genetic information processing, environmental information processing, cellular processes and organismal systems, and further predicted in 33 specific pathways ([94]Fig 4). The largest pathway group was environmental information processing in the “signal transduction”, which contained 591 genes. Fig 4. Distribution of unigene numbers for the major KEGG pathway categories in the P. trituberculatus transcriptome. [95]Fig 4 [96]Open in a new tab Differentially expressed genes (DEGs) In the current study, the DEGs were selected according to the criteria of a P value ≤ 0.05 and a |log2(FoldChange)| > 1, and the clustering of transcripts according to patterns of differential expression across samples ([97]Table 5). For the ovary, 316 DEGs were detected between the estradiol and control groups, of which 76 and 240 unigenes were up-regulated and down-regulated respectively. For the hepatopancreas, 1300 DEGs were detected between the estradiol and control groups, of which 1191 and 109 unigenes were up-regulated and down-regulated respectively. For the, brain ganglion, 669 DEGs were detected between the estradiol and control groups, of which 268 and 401 unigenes were up-regulated and down-regulated respectively. For the eyestalk, 142 DEGs were detected between the estradiol and control groups, of which 21 and 121 unigenes were up-regulated and down-regulated respectively. For the mandibular organ, 383 DEGs were detected between the estradiol and control groups, of which 233 and 150 unigenes were up-regulated and down-regulated respectively. Table 5. Summary statistics of differentially expressed genes. DEG Set DEG Number Up-regulated Down-regulated O-1 vs O-2 316 76 240 H-1 vs H-2 1300 1191 109 BG-1 vs BG-2 669 268 401 SG-1 vs SG-2 142 21 121 MO-1 vs MO-2 383 233 150 [98]Open in a new tab GO and KEGG enrichment analysis of DEGs To further assign the putative functions to DEGs, we performed GO and KEGG pathway analyses. The DEGs were assigned to 205, 488, 337, 441, and 319 GO terms in the ovary, hepatopancreas, brain ganglion, eyestalk and mandibular organ, respectively. The most significantly enriched top 20 metabolic pathways were represented for each tissue between estradiol and control group ([99]Fig 5). In the ovary, the significantly enriched pathways included “Antigen processing and presentation”, “Pancreatic secretion” and “Protein digestion and absorption”. In the hepatopancreas, the significantly enriched pathways were “Purine metabolism” and “Other types of O-glycan biosynthesis”. The most significantly enriched pathways for the brain ganglion were “Ribosome-genetic information processing” and “Biosynthesis of amino acids”. The DEGs in eyestalk were significantly clustered in “Vascular smooth muscle contraction” and “Synaptic vesicle cycle”. For the mandibular organ, the significantly enriched pathways included “Oxidative phosphorylation”, “Drug metabolism-cytochrome P450” and “Metabolism of xenobiotics by cytochrome P450”. Fig 5. KEGG enrichment analyses of DEGs. [100]Fig 5 [101]Open in a new tab The enrichment factor indicates the ratio of the DEGs number to the total gene number in a certain pathway. The size and color of the dots represent the gene number and the range of P values, respectively. The 1 and 2 represent estradiol group and control group, respectively. O: ovary; H: hepatopancreas; BG: brain ganglion; SG: eyestalk; MO: mandibular organ. Key genes in response to estradiol treatment Since the proteins and lipids are the main constituents of yolk, the process of nutrients metabolism is active during the period of ovarian development. Among the genes that were found to be differentially expressed in the estradiol-injection crab compared to the control, several candidate genes that were thought to play an important role in crustacean ovarian development and nutrient metabolism were further identified ([102]Table 6). In the ovary, the vitelline membrane outer layer 1-like protein, carboxypeptidase B, trypsin-like serine proteinase and trypsin were down-regulated by estradiol injection. In the hepatopancreas, estradiol injection up-regulated the expression levels of 5 genes (pancreatic triacylglycerol lipase-like, very low-density lipoprotein receptor, fatty acid elongase protein, alcohol dehydrogenase and lipid storage droplet protein) involved in the metabolic process. For genes involved in reproductive process in hepatopancreas, the transcripts of 9 genes significantly increased in estradiol treatment, while the expression level of heat shock protein 70 significantly decreased. In brain ganglion, estradiol injection up-regulated the expression levels of cyclin B, JHE-like carboxylesterase 2 and juvenile hormone esterase-like protein, while down-regulated the transcripts of Cytochrome P 302a1. In mandibular organ, the expression levels of JHE-like carboxylesterase 1, juvenile hormone acid methyltransferase and C-type allatostatin significantly decreased in estradiol treatment. Interestingly, 4 genes (crustacean hyperglycemic hormone 1, molt-inhibiting hormone, neuropeptide F2 and SIFamide) involved in reproductive process in eyestalk were up-regulated under estradiol injection. In addition, estradiol significantly altered the expression levels of a novel transcript encoding ovary development-related protein. Table 6. Summary statistics of differentially expressed genes. Functional category Gene id Tissue Full name FC[103]^a P-value Vitellogenesis c89341_g1 Ovary Vitelline membrane outer layer 1-like protein 0.25 7.55E-03 Reproduction c112734_g3 Hepatopancreas Heat shock protein 70 0.009 3.67E-04 c111383_g1 Hepatopancreas Vasa-like protein 485.84 2.38E-03 c93926_g2 Hepatopancreas Wnt-5b Inf 2.41E-02 c120259_g1 Hepatopancreas Wnt5a ligand Inf 3.85E-02 c116484_g2 Hepatopancreas Insulin-like receptor Inf 4.42E-04 c118618_g1 Hepatopancreas Neuropeptide receptor B3 28.87 9.09E-03 c115847_g1 Hepatopancreas Putative cytoplasmic polyadenylation element-binding protein 243.78 5.73E-05 c113803_g1 Brain ganglion Cyclin B 2.90 2.07E-02 c108894_g1 Hepatopancreas Ovary development-related protein 156.86 2.69E-04 c121312_g1 Hepatopancreas JHE-like carboxylesterase 1 21.23 8.08E-03 c62034_g1 Brain ganglion JHE-like carboxylesterase 2 282.21 3.07E-09 c118386_g3 Mandibular organ JHE-like carboxylesterase 1 0.12 1.17E-04 c92476_g1 Brain ganglion Juvenile hormone esterase-like protein 2.32 1.56E-02 c108242_g2 Mandibular organ Juvenile hormone acid methyltransferase 0.33 3.03E-02 c111705_g1 Hepatopancreas Estrogen sulfotransferase 343.62 3.03E-05 c83592_g2 Brain ganglion Cytochrome P 302a1 0.48 3.87E-02 c15627_g1 Eyestalk Crustacean hyperglycemic hormone 1 210.38 1.26E-07 c87010_g1 Eyestalk Molt-inhibiting hormone 1823 3.79E-09 c82338_g1 Eyestalk Neuropeptide F2 49.75 6.54E-04 c100581_g1 Eyestalk SIFamide 24.25 4.17E-04 c100697_g1 Mandibular organ C-type allatostatin 0.08 2.85E-03 Protein metabolism c120409_g1 Ovary Carboxypeptidase B 0.12 1.07E-20 c110626_g1 Ovary Trypsin 0.15 2.64E-19 c115002_g1 Ovary Trypsin-like serine proteinase 0.16 2.80E-16 Lipid metabolism c116241_g1 Hepatopancreas Pancreatic triacylglycerol lipase-like Inf 1.06E-02 c96963_g1 Hepatopancreas Very low-density lipoprotein receptor Inf 4.71E-03 c60052_g1 Hepatopancreas Fatty acid elongase protein Inf 1.82E-02 c110929_g1 Hepatopancreas Alcohol dehydrogenase Inf 1.21E-04 c118072_g2 Hepatopancreas Lipid storage droplet protein 168.34 2.22E-04 [104]Open in a new tab ^aFold changes (Log2 ratio) in gene expression. Inf indicate gene was expressed in the estradiol group but not in the control group. Validation of transcriptomic sequencing by quantitative PCR analysis To confirm the sequencing results, 5 DEGs involved in ovarian development and nutrient metabolism were chosen for qRT-PCR analysis using the same RNA samples. These DEGs include vitelline membrane outer layer 1-like protein, trypsin, glucose dehydrogenase, lipid storage droplet protein, and cytochrome P302a1. The results showed that the expression of all genes can be detected, and the patterns of differential expression are consistent with the results from transcriptome data ([105]Fig 6). Fig 6. Comparison of gene expression patterns obtained using RNA-Seq and qRT-PCR. [106]Fig 6 [107]Open in a new tab Notes: Log-fold changes are expressed as the ratio of gene expression after normalization to 18S. The 1 and 2 represent estradiol group and control group, respectively. Discussion Estradiol-responsive gene related to ovarian development The endocrine regulation of ovarian development, involving multiple hormonal that secreted by different organs, has been studied in many crustacean species, including crayfish, shrimp, crab, lobsters [[108]3, [109]36–[110]38]. Different organs play different roles during the reproduction of crustacean, such as, ovary are the primary reproductive organs, whose normal development is crucial for reproduction [[111]3, [112]39]. As the central organ for the absorption and storage of nutrients, hepatopancreas is also an important site for the synthesis of vitellogenin and metabolism of sex steroid hormones [[113]34, [114]40–[115]43]. Brain ganglion, eyestalk (X-organ-sinus gland complex) and mandibular organ are known to be the major endocrine organ that secrete many hormones involved in crustacean reproduction [[116]3, [117]44]. This initial study describes the transcriptomic response of P. trituberculatus to short-term estradiol injection, and these unigenes could be used to reveal new insights into the estradiol regulation of ovarian development. In the present study, a large number of key genes that may participate in the ovarian development were detected in the transcriptomic sequencing data. Vg is considered to be a biological marker of ovarian development, which is widely used in the evaluation of reproduction and endocrine disruption for both vertebrates and invertebrates [[118]45–[119]47]. However, estradiol injection did not significantly increase the expression of Vg in the transcriptome. This may be attributed to the time lag between the hormone induction and Vg expression. Similar results were also found in Scylla paramamosain [[120]48]. Another gene associated with vitellogenesis, the transcript of vitelline membrane outer layer 1-like protein was regulated by estradiol injection, which is an important protein for developing oocytes and its major function is to avoid mixing of yolk and albumen [[121]49]. Although estradiol did not significantly upregulate the gene expression related to vitellogenesis, the expression level of genes related to ovarian development were significantly altered by estradiol injection in various tissues of P. trituberculatus. Heat shock proteins are very conserved and involved in protein folding, degradation, and transportation, and their expression levels were regulated by sex steroid hormones [[122]50, [123]51]. It has been proved that Heat shock protein 70 could be involved in estrogen nuclear-initiated steroid signaling as a molecular chaperone to negatively regulate Vg gene expression in Metapenaeus ensis [[124]52, [125]53]. Therefore, the expression level of Hsp70 was decreased may be related with the promotion of vitellogenesis by estradiol. Vasa is one of the important regulatory factors that determine the development of the reproductive system [[126]54, [127]55]. In crustacean, vasa cDNAs have been reported in a variety of crustaceans, including Litopenaeus vannamei [[128]56], Fenneropenaeus chinensis [[129]57], Macrobtachium rosenbergii [[130]58], Parhyale hawaiensis [[131]59], Scylla paramamosain [[132]60] and Eriocheir sinensis [[133]61]. It has been proved that vasa are specifically expressed in gonads, and exogenous estradiol promotes the ovarian maturation by up-regulated the expression level of vasa in Sparus aurata [[134]62]. Therefore, the up-regulation of vasa expression level in this study may be related to the positive regulation of estradiol on ovarian development. The wnt signaling pathway is a conservative signaling network, which takes part in embryonic development, cell differentiation and proliferation, and the process of growth regulation [[135]63]. As a signaling factor, Wnt5 plays various roles in vertebrate development [[136]64–[137]67]. Importantly, Wnt5 could coordinate with Wnt4 to initiate the meiosis for ovarian follicular growth in mammals [[138]68]. The expression level of Wnt5a ligand and Wnt5b significantly increased after estradiol injection, indicating estradiol might regulate the ovarian development of P. trituberculatus by affecting the expression of Wnt signal pathway related genes. Further research is required to elucidate the roles of these genes in the reproductive process of crab. Progesterone is a hormone controlling the reproductive development, which not only stimulates yolk protein synthesis and ovarian maturation in Penaeus vannamei and Scylla paramamosain, but also stimulates pawning in Metapenaeus ensis [[139]69–[140]71]. It has also been reported that progesterone induces oocyte maturation by the activation of cyclin B [[141]72]. In this study, two transcripts, encoding cytoplasmic polyadenylation element binding protein and Cyclin B, involved in regulation of progesterone on oocyte maturation pathway were found to be upregulated by estradiol injection. This indicated that an interaction between estradiol and progesterone-mediated oocyte maturation pathway might exist in crab. Ovarian development is modulated by a variety of hormones, Thus, the expression of thirteen genes encoding hormones is also affected by estradiol in this study. Methyl farnesoate, a unepoxidated form of insect juvenile hormone III, is secreted by the mandibular organ in crustacean, and has obvious stimulatory function on Vg synthesis in various decapod species, including Libinia emarginata, Cancer pagurus and Procambarus clarkii [[142]4, [143]21, [144]73, [145]74]. The biosynthetic pathway of methyl farnesoate has been well established in arthropods, which is divided into two steps [[146]75]. In the first step, farnesyl pyrophosphate (FPP) is produced by the classical mevalonate pathway [[147]76]. Then, the FPP is hydrolyzed to farnesol and oxidized to farnesoic acid that will be converted to methyl farnesoate via farnesoic acid O-methyl transferase [[148]77]. However, little is known about its degradation in crustaceans. In insects, the juvenile hormone esterase, a carboxylesterase, is responsible for JH inactivation [[149]78]. Previous studies showed that methyl farnesoate is metabolized to farnesoic acid in vitro by esterases that present in crustacean tissues [[150]79]. In this study, four transcripts related to methyl farnesoate metabolism pathway were identified, including two JHE-like carboxylesterase, Juvenile hormone esterase-like protein and Juvenile hormone acid methyltransferase. Estradiol up-regulated the expression level of JHE-like carboxylesterase 1, JHE-like carboxylesterase 2 and Juvenile hormone esterase-like protein in hepatopancreas and brain ganglion, indicating that estradiol promoted the methyl farnesoate inactivation in two tissues. However, estradiol reduced the transcript of juvenile hormone acid methyltransferase and JHE-like carboxylesterase 1 in mandibular organ, which may indicate that the effect of estradiol on methyl farnesoate metabolism in mandibular organ is different from hepatopancreas and brain ganglion. As we know, Cytochrome P 302a1(CYP302a1) plays a critical role during ecdysteroid, another important reproductive hormone, biosynthesis process, that catalyzes the conversion of 5β-ketodiol to 20-hydroxyecdysone [[151]80]. In the present study, estradiol reduced the transcript of CYP302a1, indicating that an increase in concentration of estradiol may lead to a decrease of ecdysone synthesis. Moreover, estrogen sulfotransferase (EST), the enzyme responsible for the sulfonation and inactivation of estrogens, plays an important role in estrogen homeostasis [[152]81]. Previous studies have shown that the crustacean hepatopancreas is an important site for steroid hormone catabolism [[153]72, [154]82]. The expression level of EST was up-regulated in hepatopancreas indicating that the crustaceans may have a self-regulation mechanism to metabolize and excrete the excess exogenous estradiol. Similar result was also found in other crab [[155]6]. In crustaceans, crustacean hyperglycemic hormone family, a very important reproductive neuropeptides that mainly synthesized in the X-organ in the eyestalk, possess five types peptides, including crustacean hyperglycemic hormone, ion transport peptides, molting inhibitory hormone, vitellogenesis-inhibiting hormone, and mandibular organ-inhibiting hormone, which exhibit diversified physiological functions [[156]83]. Previous studies showed that crustacean hyperglycemic hormone and molting inhibitory hormone promoted the process of vitellogenesis and ovarian development in decapod crustacean [[157]17, [158]19]. Besides crustacean hyperglycemic hormone family, SIFamide, Neuropeptide F2 and C-type allatostatin were also secreted by eyestalk and cerebral ganglia. It has been reported that Neuropeptides may involve in stimulating food intake to supply energy for ovarian maturation in Scylla paramamosain, while allatostatins reduce the methyl farnesoate synthesis and hence repress ovarian maturation [[159]84]. In this present study, estradiol injection significantly up-regulated the expression of crustacean hyperglycemic hormone 1, molting inhibitory hormone, SIFamide and Neuropeptide F2, while the expression of C-type allatostatin was significantly down-regulated. It can be speculated that estradiol regulates ovarian development by affecting the balance of neuropeptides in crustaceans. Estradiol-responsive gene related to nutrition metabolism Based on previous studies, estradiol is involved in regulating metabolic processes [[160]85, [161]86]. Therefore, further understanding the changes in genes associated with protein and lipid metabolism is essential in the transcriptome of P. trituberculatus to investigate its role on reproduction of crustacean. In this present study, three protease genes including trypsin, trypsin-like serine proteinase and carboxypeptidase B, were down-regulated by estradiol injection in the ovary of P. trituberculatus. Trypsin and carboxypeptidase are enzymes that cleave the peptide bond of different amino acid residue, while trypsin is reported to be involved in regulating the yolk degradation and yolk degradome activation [[162]87, [163]88]. Proteins are one of the main constituents of yolk [[164]2, [165]89]. It was speculated that estradiol indirectly promotes the process of yolk formation by reducing the expression of protein catabolism genes. Six candidate genes for lipid metabolism were up-regulated by estradiol in the hepatopancreas of P. trituberculatus, including pancreatic triacylglycerol lipase-like, very low-density lipoprotein receptor, fatty acid elongase protein, lipid storage droplet protein and alcohol dehydrogenase. In crustacean, lipid accumulating oocytes provides fuel for the biosynthetic processes of oogenesis and vitellogenesis and is later used by developing larvae [[166]90]. Lipid storage droplets (LSDs) are organelles that accumulate lipid for both long- and short-term storage occurring in many animal cell types [[167]91]. A recent study in Penaeus monodon indicated that the expression level of LSD significantly up-regulated during ovarian maturation [[168]92, [169]93]. In this study, estradiol promoted the expression of LSD in hepatopancreas, indicating that estradiol can provide energy for ovarian development of P. trituberculatus by promoting the accumulation of LSD in the hepatopancreas. Although above results have shown that exogenous estradiol had positive effects on genes relative to ovarian development of P. trituberculatus, the regulation mechanism of estradiol is largely unknown [[170]5, [171]38, [172]94]. In vertebrate animals, it is well-known that cellular estrogen signaling was mediated primarily via the estrogen receptor (ER), a family of nuclear hormone receptor-type transcription factors [[173]95]. Unfortunately, to date, the ER gene has not been functionally identified for any crustacean species, which may have been lost during the evolution of arthropods [[174]96]. However, the estrogen related receptor (ERR) with high sequence homology with the vertebrate ER has been found in P. trituberculatus [[175]97]. In this study, estradiol injection altered the transcription of ERR in five tissues of P. trituberculatus ([176]S1 Table). Previous study has shown that ERR involved in regulating the gonadal development and cellular energy balance in vertebrate [[177]98, [178]99]. Our previous results showed that the expression level of ERR was also altered when long-term injection of estradiol into P. trituberculatus [[179]26]. In future studies, it is necessary to in-depth study the role of ERR in the regulation of estradiol in crustaceans. Conclusion This study investigates the transcriptome sequencing of ovary, hepatopancreas, brain ganglion, eyestalk and mandibular organ of a swimming crab within a relatively short period after estradiol injection (24 h, 0.1μg g^-1 crab weight). Through this study, some key genes in correlation with the reproduction and nutrition metabolism were significantly affected by estradiol. These results will serve as important resources for future experiments that further investigate the role and regulation of estradiol in P. trituberculatus. However, the function of estradiol-responsive gene in ovarian development and their cross-talking warrants further investigation. Supporting information S1 Table. The expression level of receptor and steroidogenic-related genes in five tissues. (XLSX) [180]Click here for additional data file.^ (10.3KB, xlsx) Data Availability All relevant data are within the paper and its Supporting Information files. Funding Statement This study was funded by the two projects (No. 41276158 and No. 41606169) from the Natural Science Foundation of China and a special research project (No. SZ-LYG2017019) for North Jiangsu area from Science and Technology Department of Jiangsu Province. Infrastructure costs were partially supported by the research project (No. A1-2801-18-1003) for high-level university in Shanghai from Shanghai Education Commission and Collaborative Innovation Project for Mari-culture industry in East China Sea from Ningbo University. References