Abstract Banana shrimp (Fenneropenaeus merguiensis) is an economically important species in Thailand owing to the high value of globally exported frozen brine shrimps. However, the regulatory mechanisms governing spermatogenesis and testicular development in this species are poorly understood. High-throughput RNA sequencing was used to investigate the mechanisms and regulated genes involved in testis development using transcriptome profiling of juvenile and adult banana shrimp testes. Differentially expressed genes (DEGs) in these two libraries were identified and quantified to confirm gene expression. DEGs were found in 7,347 genes, with 4,465 upregulated and 2,882 downregulated. Some of these genes were designated as candidate genes, and six specific DEGs, including PRM1, SPATA20, Sry, SSRF, Sxl, and Tra-2c, were selected to confirm the reliability of the RNA-seq data using qPCR. Moreover, six non-DEGs were chosen based on testis-specific and regulatory genes that support a specific function in spermatogenesis and testis development in this species, including Dsx, Gfra2, IAG, Sox9, Sox13, and Sox14A. Furthermore, Sry, Sox14A, Sox14B and SPATA20 were identified in early stages (nauplius-postlarvae) of shrimp development to provide more information involving testes formation and development. The transcript data from this study could differentiate a group of genes required at the early and late stages of testis development and both sets of testis development. Therefore, this information would help in manipulating each stage of testicular development. Introduction Penaeid shrimp testes contain a testicular lobule where spermatogenesis occurs. Spermatogenesis is the process of haploid spermatozoa production and is divided into two stages: spermatogenesis and spermiogenesis. Spermatogenesis is the process of diploid spermatogonia and haploid spermatid production, and spermiogenesis is transforming mature spermatids into spermatozoa [[26]1]. Many studies in crustaceans have focused on spermiogenesis because the most representative process data are available for each species, including Sicyonia ingentis [[27]2], Parapenaeus longirostris [[28]3], Macrobrachium rosenbergii [[29]4], Fenneropenaeus chinensis [[30]5], Litopenaeus vannamei [[31]6], and Penaeus monodon [[32]7] However, spermatogenesis in Fenneropenaeus merguiensis is still unknown. F. merguiensis is widely distributed in tropical and subtropical waters, especially on the northern coast of Malaya, the west coast of the Philippines, Indonesia, Australia, and Thailand. This species has been observed to mature at around 6–7 months, with a slight difference between males and females [[33]8]. Owing to its high commercial importance in many areas, such as the Persian Gulf and Pakistan [[34]9], the Philippines, Thailand, Indonesia, and Australia, a better understanding of its reproductive development, especially spermatogenesis and testis development, is necessary to establish the high quality and number of spermatozoa produced. Previous transcriptome analyses using RNA-seq have been conducted in many shrimps to improve broodstock reproduction for industrial farming. Transcriptomic data from shrimp testes have been widely reported; for example, the gene expression profile of Macrobrachium nipponense testes was analyzed to identify candidate genes involved in male reproduction. Fem-1, Fem-2, and Fem-3 promote male development in the sex-determining protein complex [[35]10]. The de novo transcriptome of the Japanese mantis shrimp (Oratosquilla oratoria) was also examined. Some regulated gene expression in males has been identified, including cathepsin I and cathepsin D-like [[36]11]. The transcriptome of P. monodon was recently studied, and various organs, particularly the testes, were examined, providing helpful information for gene expression profiles involved in testis regulation [[37]12]. Several regulatory genes involved in testis differentiation and maturation have been studied; for example, cyclin A and cyclin B [[38]13], cell division cycle 2 (Cdc2) [[39]14], mitogen-activating protein kinase 1 (MAPK1) [[40]15] in P. monodon; activated protein kinase C1 (RACK1) [[41]16], cell apoptosis susceptibility (FcCAS) [[42]17] in F. chinensis; vasa-like [[43]18], spermatogonial stem-cell renewal factor (SSRF), spermatogenesis-associated proteins (SPATA2), sperm protamine P2 (PRM2), testis-specific serine proteases (TESSP) in L. vannamei [[44]19], and Doublesex and Mab-3 Related Transcription Factor 1 (Dmrt1) in Palaemon serratus [[45]20]. Currently, the regulated genes involved in testicular development and spermatogenesis in F. merguiensis are poorly understood. The transcriptomic profiles of juvenile and mature F. merguiensis testes were studied for candidate regulatory genes and involved pathways using high-throughput RNA-seq technology. The candidate genes in testis development were selected from the ≥log 2-fold change in the compared stages, and the expression of the selected candidate genes was confirmed using qPCR. The outcome of this study provides crucial information on the regulatory gene expression involved in banana shrimp testis development. Materials and methods Ethics declarations This study was approved and conducted according to the Institutional Animal Care and Use Committee (IACUC) of the Prince of Songkla University. Sample collection Juvenile and adult male shrimp were purchased from the Trang and Nakhon Si Thammarat provinces, respectively, of Thailand. The shrimp were reared from postlarva day 5 until the testis began to form an external male sex organ were classified as juvenile, approximately 6–7 cm in length and 1.5–2.0 g in weight (about postlarva day 114^th), whereas adult shrimp were wild-caught, about 16–17 cm in length, and 28.0–30.0 g in weight. Three male shrimp from each group were dissected to harvest the testes. The excised testes were divided into two pieces. One piece was fixed in Davidson’s solution for histological examination, whereas the others were immersed immediately in TRIzol reagent (GIBCO BRL, Grand Island, NY, USA) for RNA preparation. Histological analysis of testis samples Testes fixed in Davidson’s solution were processed and embedded in paraffin wax. The samples were cut into 5 μm thick sections and stained with hematoxylin and eosin to examine the histological structure of juvenile and adult shrimp testes. An Olympus BX51 microscope (Tokyo, Japan) with 40X objectives was used to capture images of all samples. RNA isolation and cDNA library construction Total RNA was isolated from each sample using TRIzol reagent (GIBCO BRL), following the manufacturer’s instructions. Potentially contaminated genomic DNA was removed using RQ1 RNase-free DNase (Promega, Madison, WI, USA). The RNA concentration and quality were determined using an Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit, Agilent, Santa Clara, CA, USA). For the construction of cDNA libraries, three RNA samples from each group were pooled. The first step in the library construction process is mRNA enrichment via oligo dT selection or rRNA depletion. RNA samples were fragmented in the second step before being reverse-transcribed to double-stranded cDNA using the N6 random primer. Before 3′-adenylation, the synthesized cDNA was repaired. Adaptors were then ligated to 3′-adenylated cDNA ends. The ligation products were purified, followed by an amplification step to enrich the purified cDNA template. The PCR products were denatured using heat, and the single-stranded DNA was cyclized using splint oligo-and DNA ligase. Each sample library was sequenced on the BGISEQ-500 platform at Beijing Genomics Institute (BGI, Shenzhen, China). The raw data for the assembled transcriptome was submitted to the Sequence Read Archive (SRA) under the accession number PRJNA961319. De novo assembly and functional annotation Following sequencing, low-quality reads (more than 20% of base qualities were less than 10), reads with adaptors, and reads with unknown bases (N bases greater than 5%) were filtered using SOAPnuke to obtain clean reads. Clean reads were assembled into unigenes using Trinity (v2.0.6) ([46]https://github.com/trinityrnaseq/trinityrnaseq/releases/tag/v2.0.6 ). Using the TGICL software (v2.0.6) ([47]https://sourceforge.net/projects/tgicl/files/tgicl/), gene family clustering was performed to obtain the final unigenes. The unigenes were then aligned to NT, NR, Clusters of Orthologous Group (COG), and SwissProt databases to perform functional annotation using BLAST (v2.2.23) and Diamond (v0.8.31) software ([48]https://github.com/bbuchfink/diamond). The Gene Ontology (GO) annotation was performed using Blast2GO (v2.5.0), providing biological functions at the molecular, cellular, and tissue system levels. Pathway enrichment analysis was performed using the online KEGG Automatic Annotation Server (KAAS) ([49]http://www.genome.jp/kegg/kaas/) to understand the high-level functions of metabolic pathways. Differential expression genes (DEGs) analysis Based on the assembly results, all of the clean reads from each sample were mapped to unigenes using Bowtie2 software (v2.2.5) and the gene expression level was calculated using RSEM (v1.2.12). The relative expression of a transcript was calculated by counting fragments per kilobase of transcript per million mapped reads (FPKM). A Poisson distribution was performed to identify differentially expressed genes (DEGs) between juvenile and adult shrimp testes, with the false discovery rate (FDR) adjusted to ≤ 0.001. A fold change ≥ 2 with the log2 ratio of DEGs between the adult testis and juvenile testis (adult testis/juvenile testis) ≥ 1 was used to judge the significance of DEGs. Quantitative real-time PCR (qRT-PCR) validation Genes involved in testis development were selected using the KEGG pathway annotation, including six specific DEGs; Sperm Protamine P1-like (PRM1), Spermatogenesis-associated protein 20 (SPATA20), Sex-determining region Y (Sry), SSRF, Sex-lethal (Sxl), and Transformer-2c (Tra-2c), as well as six non-DEGs; Doublesex (Dsx), GDNF family receptor alpha 2 (Gfra2), Insulin-like androgenic gland hormone (IAG), Sry-box transcription factor 9 (Sox9), Sry-box transcription factor 13 (Sox13), and Sry-box transcription factor 14A (Sox14A). Gene expression data from RNA-seq were validated using qRT-PCR. Briefly, total RNA from juvenile (n = 3) and adult shrimp (n = 3) testes were isolated and pooled equally, and reverse transcription was performed to synthesize cDNA using AMV reverse transcriptase (Promega) according to the manufacturer’s protocol. qRT-PCR was performed with three technical replicates in the Mx3000P qPCR system (Stratagene, San Diego, CA, USA) using the FastStart Universal SYBR Green Master (Roche, Germany). β-actin was used as an endogenous reference to eliminate sample-to-sample variations and normalize changes in specific gene expression. The amplification procedure was as follows: pre-denaturation at 94 °C for 3 min followed by 40 cycles of denaturation at 94 °C for 30 s, annealing at the required temperature for 30 s, and extension at 72 °C for 45 s. The relative expression levels of the DEGs were calculated using the 2^−ΔΔCT method, and data were expressed as mean ± SD. Analysis of gene expression involved in testes development Some genes involved in testes development, including SPATA20, SOX14A, SOX14B, and Sry were examined for expression in different stages of shrimp, including nauplius, zoea, postlarva (PL) 5, 10, 20, 30, 50, 70, 90, and adult shrimp (about 5–6 months old). RNA was extracted from whole shrimp from nauplius to PL90, while only testes from adult shrimp were collected. The cDNA was then converted using AMV Reverse Transcriptase (Promega, Madison, WI, USA) and amplified for each gene. The condition for PCR was performed as previously described in qPCR. The expression of all genes was examined using 8% polyacrylamide gel electrophoresis. Results Histological characteristics of juvenile and adult testis During the process of spermatogenesis, histological examination revealed that sperm developed into several stages, including spermatogonia (Sg), spermatocyte (Sc), spermatid (St), and spermatozoa (Sz). Spermatogonia (Sg) are cells found in the seminiferous tubules with a pale and rounded nucleus containing chromatin. They were the spermatogenic cells attached to the basal lamina of the seminiferous tubule and were engaged by the surrounding of somatic stem cell (SSC). After mitosis division, Sg was divided and transformed into Sc, a medium-sized cell with a slightly condensed spherical nucleus. Then, Sc was divided into St with a smaller and more condensed spherical nucleus after the second meiotic division. High-condensed chromatin was appeared during the maturation of St into Sz, and the nucleus becomes a distorted sickle shape and darkens. The difference in morphological appearance between juvenile and adult testes was classified. The juvenile testis ([50]Fig 1A) contained more Sg cells than the other types, whereas the adult testis ([51]Fig 1B) contained multiple stages of sperm and was primarily St and Sz. In addition, adult testis had more total cells than juvenile testis. Fig 1. Histological analysis of the juvenile and adult testes. [52]Fig 1 [53]Open in a new tab The histological analysis of the juvenile testis (A) and adult testis (B) showed various cell types in the lobules of the testes. SSC, Sg, Sc, St, and Sz indicated somatic stem cells, spermatogonia, spermatocyte, spermatid, and spermatozoa, respectively. mRNA expression profiles of the juvenile and adult testes Genes involved in the development of banana shrimp testes were identified by constructing two juvenile and adult testis libraries using the BGISEQ-500 platform. More than 77.13 and 80.63 million raw reads of juvenile and adult testis were generated. After low-quality, adaptor-polluted, and high content of unknown base (N) read, they were filtered, and the total clean reads of juveniles and adults were obtained with 69.25 and 70.18 million reads, respectively ([54]Table 1). Trinity software was used to perform de novo assembly with clean reads, and TGICL software was used to remove abundance and obtain unigenes from cluster transcripts. Table 1. Sequencing statistics analysis of juvenile and adult testes transcriptome. Sample Juvenile testis Adult testis Sequencing results Total Raw Reads (M) 77.13 80.63 Total Clean Reads (M) 69.25 70.18 Clean Reads Ratio (%) 89.78 87.04 Clean Reads Q20 (%) 96.53 96.83 Total unigene number with N50 38,216 with 2,295 40,789 with 2,610 Alignment results Total Bases 10,387,722,900 10,526,715,600 Total Reads 69,251,486 70,178,104 Total Mapped Reads 58,811,734 58,787,076 Unique Mapped Reads 36,037,600 34,150,556 Mapping Ratio (%) 84.92 83.77 Annotation results NR 24,290 (45.71%) NT 16,566 (31.17%) SwissProt 19,612 (36.91%) KEGG 20,203 (38.02%) COG 17,993 (33.86%) GO 7,324 (13.78%) InterPro 16,194 (30.47%) [55]Open in a new tab DEGs annotation and pathway analysis The expression levels of all genes were detected using PossionDis Algorithm to verify the DEGs between the juvenile and adult testes. A Scatter plot, Venn diagram, and heatmap plot were generated to show the distribution of the DEGs. In total, 7,347 genes were identified as DEGs, of which 4,465 were upregulated, and 2,882 were downregulated ([56]Fig 2A). A Venn diagram indicates the number of particular expression genes in one sample and the number expressed in both models ([57]Fig 2B). In addition, gene expression was higher in adult testes (7,007 genes) than in juvenile testes (4,029 genes). Fig 2. Summary of DEG plot. [58]Fig 2 [59]Open in a new tab A Scatter plot of DEGs (A) shows the upregulated genes with red dots, down-regulated genes with blue dots, and no-DEGs with grey dots. The vertical axis represents the log[10] of the expression level in the adult testis gene, whereas the horizontal axis represents the log[10] of the juvenile testis gene expression level. Venn diagram (B) of gene expression was analyzed between juvenile and adult testes. The result of NR annotation was investigated for species distribution to evaluate the evolutionary conservation. A total of 24,290 genes were mapped to the NR database using BLASTX. The top three species with the most abundant matched transcripts are Hyalella Azteca (25.79%), Zootermopsis nevadensis (6.32%), and Limulus Polyphemus (2.19%), all of which are members of the Arthropoda phylum, to which F. merguiensis also belongs. Among others abundant in the Penaeidae family, Litopenaeus vannamei was the most abundant shrimp (1.66%), followed by P. monodon (1.60%) and F. chinensis (0.67%) ([60]Fig 3). Fig 3. The homologous species distribution is annotated against the NR database. [61]Fig 3 [62]Open in a new tab The number showed the percentages of the best blast hits of unigenes among other species. Gene ontology (GO) analysis of DEGs was classified for functional enrichment based on three ontologies: biological process, cellular component, and molecular function ([63]Fig 4). A total of 22 GO terms belonged to the biological process category, which was mainly focused on cellular processes (505 genes), metabolic processes (451 genes), and biological regulation (229 genes). For cellular components, approximately 18 GO terms were members of "membrane" (684 genes), "membrane part" (667 genes), and "cell" (493 genes). The molecular function was the last ontology, which included 10 GO terms. The top three counts were "binding" (690 genes), "catalytic activity" (580 genes), and "transporter activity" (61 genes). Fig 4. Bar plot of enriched GO terms of DEGs. [64]Fig 4 [65]Open in a new tab The vertical axis represents the name of each Go term, whereas the horizontal axis represents the number of genes. Each color refers to the different terms of GO. In total, 20,203 genes with KEGG pathway annotations were mapped to 332 pathways. The most prominent pathway member was the metabolic pathway, which included 2,481 genes (12.28%). However, only 4,703 genes were identified as DEGs. The pathway classification of DEGs has been classified into seven categories: metabolism, genetic information processing, environmental information processing, cellular processes, organismal systems, human diseases, and drug development. Biological pathways involved in the testis development of F. merguiensis were discovered, including the GnRH signaling pathway, insulin signaling pathway, Wnt signaling pathway, Notch signaling pathway, aldosterone synthesis and secretion, TGF-beta signaling pathway, and thyroid hormone signaling pathway. DEGs involved in spermatogenesis and testis development Spermatogenesis is the process of haploid spermatozoa production. The discovery of marker genes for this process will improve our understanding of the regulatory processes that govern spermatogenesis. Some functional genes involved in late spermatogenesis were upregulated in the adult shrimp testis compared to the juvenile shrimp testis, including Argonaute 4 (AGO4), Crustacean hyperglycemic hormone (CHH), Feminization-1b (Fem-1b), Insulin receptors (Irs), Nanos (Nos), Paired Box 7 (PAX7), PRM1, SPATA20, SSRF, and Tra-2, while others were downregulated, including Sxl, Sry, and Tra-2c ([66]Table 2, [67]S1 Table). Table 2. Some functional genes involved in spermatogenesis and testis development. Gene Description log2FoldChange (adult/juvenile) P-value FDR Up/Down AGO4 Argonaute 4 1.562 0.00E+00 0 Up CHH Crustacean hyperglycemic hormone 3.665 0.00E+00 0 Up Dmrt1 Doublesex- and mab-3-related transcription factor 1 −4.644 5.65E-02 1.21E-01 - Dsx Doublesex 5.459 5.20E-01 6.49E-01 - Fem-1b Feminization-1b 1.065 3.20E-05 1.22E-04 Up IAG Insulin-like androgenic gland 2.585 2.65E-01 3.98E-01 - IGFBP Insulin-like growth factor binding protein 0.954 2.00E-304 2.29E-302 - Irs Insulin receptors 5.248 7.87E-12 5.06E-11 Up Nos Nanos 1.126 6.17E-06 2.55E-05 Up PAX7 Paired Box 7 binding protein 1 2.492 2.51E-46 4.55E-45 Up PRM1 Sperm Protamine P1-like 4.755 2.67E-15 2.05E-14 Up Sox13 SRY-related HMG-box containing transcription factor 13 1.281 4.37E-02 9.62E-02 - Sox14A SRY-related HMG-box containing transcription factor 14A 0.867 4.87E-01 6.19E-01 - Sox14B SRY-related HMG-box containing transcription factor 14B 0.347 4.17E-03 1.18E-02 - SPATA20 Spermatogenesis-associated protein 20 2.158 1.42E-07 6.80E-07 Up SSRF Spermatogonia stem-cell renewal factor 1.520 6.33E-06 2.61E-05 Up Sxl Sex-lethal −2.388 0.00E+00 0 Down Tra-2 Transformer-2 1.322 3.81E-82 1.17E-80 Up BMP7 Bone morphogenetic protein 7 −0.240 3.65E-01 5.21E-01 - Cdc2 Cell division cycle 2 −0.437 1.61E-24 1.76E-23 - GEM Gem-associated Protein 2-like Isoform X1 −0.564 8.07E-12 5.19E-11 - Gfra2 GDNF Family Receptor Alpha 2 −0.336 4.65E-02 1.01E-01 - GnRHr Gonadotropin-releasing hormone receptor −1.271 3.12E-01 4.58E-01 - HSP90 Heat shock protein 90 −0.496 0.00E+00 0.00E+00 - Pcna Proliferating cell nuclear antigen −0.755 4.76E-78 1.39E-76 - Sox5 SRY-related HMG-box containing transcription factor 5 −0.501 4.18E-04 1.39E-03 - Sox6 SRY-related HMG-box containing transcription factor 6 −0.209 4.26E-02 9.40E-02 - Sox8 SRY-related HMG-box containing transcription factor 8 −0.186 9.56E-01 9.85E-01 - Sox9 SRY-related HMG-box containing transcription factor 9 (promoter region) −7.066 2.35E-01 3.86E-01 - SPATA5 Spermatogenesis-associated protein 5 −0.923 1.99E-22 2.05E-21 - Sry Sex-determining region Y −2.445 7.19E-05 2.64E-04 Down Tra-2c Transformer-2c −1.346 8.47E-13 5.74E-12 Down VASA DEAD box family of RNA helicases −0.363 1.18E-89 3.97E-88 - [68]Open in a new tab In addition, several genes had a trend of increased expression, such as Dsx, IAG, Insulin-like growth factor binding protein (IGFBP), Sox13, Sox14A, and Sox14B. In contrast, genes had a trend to reduce expression during the testes development, including Bone morphogenetic protein 7 (BMP7), Cdc2, Gem-associated Protein 2-like Isoform X1 (GEM), Gfra2, Gonadotropin-releasing hormone receptor (GnRHr), Heat shock protein 90 (HSP90), Proliferating cell nuclear antigen (Pcna), SRY-related HMG-box containing transcription factor 5 (Sox5), SRY-related HMG-box containing transcription factor 6 (Sox6), SRY-related HMG-box containing transcription factor 8 (Soc8), Sox9, Spermatogenesis-associated protein 5 (SPATA5), and DEAD box family of RNA helicases (VASA). Although these genes were not significantly differentially expressed between juvenile and mature testes (P>0.05), the results indicated that the expression of these genes is essential in both the juvenile and mature stages of testes development. Validation of RNA-Seq data by qPCR qPCR was used to clarify the expressed genes in the transcriptome data. The male-specific gene was selected from [69]Table 2, including upregulated genes such as PRM1, SPATA20, and SSRF, as well as downregulated genes such as Sry, Sxl, and Tra-2c ([70]Fig 5A). Moreover, non-significantly different expressed genes, including Dsx, Gfra2, IAG, Sox9, Sox13, and Sox14A were also chosen for validation ([71]Fig 5B). The majority of them corresponded with transcriptome data, with the exception of Sox9, which was a +1-fold change determined by qPCR and a -7-fold change determined by transcriptome analysis. Fig 5. Validation of RNA-Seq using qPCR. [72]Fig 5 [73]Open in a new tab The expression of the selected DEGs (A) and non-DEGs (B) from RNA-Seq data was compared to qPCR in the fold change of adult/juvenile testis. The log2FC values of qPCR are shown as the means ± SD. Expression of genes involved in male sex differentiation and testes development Some selected genes such as SPATA20, Sox14A, Sox14B, and Sry were identified from each stage of shrimp from nauplius to adult shrimp to demonstrate the potential function of genes involved in male sex differentiation and testes development. The presence of these genes at an early stage, as nauplius suggested that they may play an important role in male-sex differentiation and development, however, this needs to be investigated further using other methods, such as gene knockdown or knockout, to determine the true effect of these genes’ defection on sperm synthesis. Discussion This transcriptome data analysis between juvenile and adult testis revealed mRNA expression profiles of the juvenile and adult testes that specific gene expression was required for each stage. The upregulated gene expression for spermatogenesis in adult testis, for example, SPATA20, SSRF, Nos, PRM1 ([74]Table 2, [75]S1 Table). SPATAs have been identified as testis-specific genes that regulate apoptosis during zebrafish spermatogenesis [[76]21]. The upregulation of SPATA20 was greater in adult testes than in juvenile testes in both transcriptome and qPCR analyses in this study. Furthermore, its expression was found to be consistent from larva to adult, supporting the testis development. Other SPATAs genes with testis-biased expression patterns were discovered, including SPATA5L1, SPATA6, and SPATA18, but there was no difference between juvenile and adult testes. However, SPATA2, SPATA5, and SPATA20 were found to be essential genes due to their expression was detected in the later stages of spermatogenesis in L. vannamei [[77]19]. Interestingly, the expression of SPATA20 was examined and found from the early stage at nauplius to the adult shrimp, suggesting that this gene may be involved in the process of sperm differentiation and formation, as the lack of this gene in mice resulted in sperm number reduction and abnormal of sperm morphology [[78]22]. Spermatogonia stem cells are known to play a role in self-renewal and sperm differentiation [[79]23]. Many investigations of SSRF were reported in mammals, but there have been few findings in crustaceans [[80]24–[81]26]. The high expression of SSRF allows for the renewal of the spermatogonia number in the adult testis. SSRF is named after its role in self-renewing cells in male fertility spermatogenesis, which allows stem cell populations to proliferate [[82]25]. Therefore, the expression of SSRF was increased in the adult testes rather than in juvenile testes in this study, indicating that it plays a role in testicular development. In addition, Nos genes express the conserved zinc-finger RNA-binding proteins that play roles in maintaining germline stem cell function [[83]27–[84]29]. Nos gene expression was found to be higher in the adult testes than in the juvenile testes to maintain germ cells in the adult. PRM1, which involves sperm development, also increased expression in the adult testis of this study. SSRF and PRM1 expression results relate to the histological of the adult testis ([85]Fig 1), which showed a lot of spermatogonia population and spermatozoa were found. Fem-1b was upregulated in the adult testis of F. merguiensis, which indicated a requirement for spermatogenesis. In addition, Fem-1b, a homolog of Fem-1, was found to have high expression in early embryonic development and in the testis, ovary, hepatopancreas, and muscles of E. sinensis [[86]30, [87]31]. Therefore, Fem-1b was suggested as a maternal gene expression during the early embryonic stage, and it is required during the late stages of gonadal development [[88]30]. The result of this study corresponded with the detection of Fem-1 in spermatogonia of P. vannamei [[89]32]. This study also showed that high Sxl expression was required for juvenile testis development. Sxl was studied in several decapods, such as the Redclaw crayfish: C. quadricarinatus [[90]33], E. sinensis [[91]34] and P. vannamei [[92]35]. They suggested that Sxl does not involve in sex determination, but it is required for embryonic development, sexual differentiation, and gonad development. PAX7, AGO4, CHH, and Irs also upregulated genes that supported the development of adult testes in this study. As previously reported in mice, PAX7 is present in the testes at birth and plays a role in spermatogenesis, particularly in the spermatogonial stem cell population [[93]36]. AGO4 encodes a protein that contains PAZ and PIWI domains and plays gene regulation via RNA interference and short-interfering-RNA-mediated gene silencing [[94]37]. A knockdown study of this gene in P. monodon [[95]38] using RNAi found that reducing this gene in transcript level resulted in testicular maturity and a decrease in the number of spermatogonia in the spermatophore, revealing the critical role in controlling spermatogenesis. Generally, CHH is known to elevate circulating glucose under stressful conditions and has multiple functions, such as ecdysteroidogenesis, osmoregulation, and vitellogenesis. Some CHH neuropeptides were suggested to have an extra function in reproduction at the adult stage [[96]39]. The function involved in testis development was discovered while studying this gene in M. nipponense [[97]40]. CHH was found significantly upregulated to provide energy (glucose) for adult testes development and corresponded with upregulated Irs expression for importing glucose into cells ([98]Table 2). In contrast, Tra-2c significantly down expression in the adult testis. The Tra-2 gene has been examined in various crustaceans, including P. monodon [[99]41], F. chinensis [[100]42], M. nipponense [[101]43], and C. quadricarinatus [[102]44]. Although the nucleotide sequences are similar to D. melanogaster [[103]45], the alternative splicing patterns are quite different, implying a different mechanism of sex differentiation between crustaceans and insects, which is dependent on the alternative splicing patterns of the pre-mRNA of genes involved in sex determination. Three isoforms of Tra-2 have been identified in F. chinensis: Tra-2a, Tra-2b, and Tra-2c. The transcript of Tra-2c was significantly increased at the mysis stage and showed a significantly higher expression level in females than in male. Therefore, Tra-2c was suggested to be involved in the female sex determination of F. chinensis [[104]42]. In this transcriptome, Tra-2c was more abundant in juvenile testes than in adult testes. This phenomenon may be explained by the fact that this gene may require in early testicular development, and the expression level in the larval stage or postlarval stage is interesting to be detected for clarification of Tra-2c in female sex determination in F. merguiensis. The Sox gene family comprises a group of genes that encode a Sry-like high-mobility group (HMG) box, a male-determining gene. Many transcription factors in the Sox gene family play important roles in developmental processes, such as neurogenesis, sex determination and differentiation, and testis development [[105]46]. Sry and Sox9 (non-significant DEG) decreased expression in the adult testis, which means they have a role during the early testis development. Sry has been reported to turn off ovarian genes, whereas it turned on testicular genes, especially Sox9. In addition, Sox9 turned on testicular genes and turned off ovarian genes [[106]47]. In mice, Sry and Sox9 have important roles in repressing ovarian and activating testicular differentiation genes for testis cord formation. The mechanism of the activities was reported by repressing the WNT/b-catenin transcriptional activities of ovarian-differentiating genes in a granulosa cell line [[107]48]. In this study, the detection of Sry at nauplius was expressed and extended to the adult stage, implying that this gene is involved in sex determination, which was consistent with previous findings in mice that this gene begins expressed after the genital ridges first appear and remains expressed until the first morphological signs of testis differentiation become apparent [[108]49]. Interestingly, the Sry gene was discovered with two approximate sizes in some stages of shrimp development ([109]Fig 6), requiring further investigation into the real specific function to the stage of development. Fig 6. Examining the expression of genes involved in testis development included SPATA20, Sox14A, Sox14B, and Sry, at different stages in shrimp. [110]Fig 6 [111]Open in a new tab Lane M: 100 bp DNA ladder, lanes 1–10: amplification of gene from each stage of shrimp including nauplius, zoea, postlarva 5, 10, 20, 30, 50, 70, 90, and adult, lane N: negative control containing sterile water. Sox members have been found in various crustaceans, including M. nipponense [[112]50], P. vannamei [[113]19], and P. serratus [[114]20]. Sox5, Sox14, and Sox15 are found in the testes of P. serratus, with Sox5 and Sox14 involved in male sex differentiation, whereas Sox15 has never been identified as a testis gene [[115]20]. This study found Sox5, Sox6, Sox8, Sox9, Sox13, Sox14A, and Sox14B involved in the early stage of testis development. Since the expression of Sox14A and Sox14B was detected in the nauplius and continued into adulthood, it was hypothesized that both might be involved in male-sex differentiation and sperm development, based on previous findings in S. paramamosain that high expression of Sox14 was reported in mature sperm stage than in spermatocyte and spermatid stages during the testis development [[116]51]. Additionally, the highest expression of Sox14 was detected in the fertilized egg and early embryogenesis, which had an extensive development of organs such as the optic ganglion, appendicular ganglion, and heart [[117]52, [118]53], implying that Sox14 may be involved in the embryonic development and sperm maturation and in this species. Some of the non-DEGs found in this transcriptome, including Dsx and IAG, are known to be involved in spermatogenesis and testes development [[119]54, [120]55], but after the expression was validated, upregulation of them was detected in adult testes rather than juvenile testes, indicating a potential role in testes development. In contrast, Gfra2 was discovered in the transcriptome to be non-DEG, but its expression was found to be downregulated after the expression was investigated. According to this result, Gfra2 may be involved in the early stages of spermatogenesis. This transcriptome study focused on a better understanding of the regulatory genes and factors controlling spermatogenesis and testis development in F. merguiensis due to the lack of information on this species. Most of the results showed a point in the expression of the genes involved. The effect of these genes on testis development in F. merguiensis requires further investigation. Conclusion Comparative transcriptome analysis of juvenile and adult testes from F. merguiensis was performed. Several candidate genes involved in spermatogenesis and testis development were identified using GO terms and biological pathways, with some differences in expression. The transcript data identified group the genes required at the early (Sry, Tra-2c, Gfra2, Sox9, and Sxl) and late stages (AGO4, CHH, Fem-1, Irs, Nos, PRM1, SPATA20, SSRF, Tra-2, IAG, Sox13, and Sox14A) of testis development. Therefore, this information would be useful for manipulating each stage of testicular development. In addition, some gaps were found in the knowledge about this shrimp, such as whether Dsx regulates testis development via IAG signaling and whether Tra-2c is involved in testis development. These results imply that further investigation of these genes will be useful in the regulation of spermatogenesis and testis development in this species. Supporting information S1 Table. Additional data of some functional genes involved in spermatogenesis and testis development. (DOCX) [121]Click here for additional data file.^ (22.6KB, docx) Data Availability The data are all contained within the paper. Funding Statement 1. NSRF via the Program Management Unit for Human Resources & Institutional Development, Research, and Innovation [grant number B05F630026], granted to W. Chotigeat. 2. Prince of Songkla University for funding project contract no. SCI6405032c, granted to W. Chotigeat. 3. Prince of Songkla University, for the Postdoctoral Fellowship to U. Saetan. The funders had no role in study design, data collection, analysis, the decision to publish, or manuscript preparation. References