Comparative Genomics Studies on the dmrt Gene Family in Fish

Transforming growth factor-β (TGF-β): A master signal pathway in teleost sex determination

Sex determination and differentiation in fish has always been a hot topic in genetic breeding of aquatic animals. With the advances in next-generation sequencing (NGS) in recent years, sex chromosomes and sex determining genes can be efficiently identified in teleosts. To date, master sex determination genes have been elucidated in 114 species, of which 72 species have sex determination genes belonging to TGF-β superfamily. TGF-β is the only signaling pathway that the largest proportion of components, which including ligands (amhy, gsdfy, gdf6), receptors (amhr, bmpr), and regulator (id2bby), have opportunity recognized as a sex determination gene. In this review, we focus on the recent studies about teleost sex-determination genes within TGF-β superfamily and propose several hypotheses on how these genes regulate sex determination process. Differing from other reviews, our review specifically devotes significant attention to all members of the TGF-β signal pathway, not solely the sex determination genes within the TGF-β superfamily. However, the functions of the paralogous genes of TGF superfamily are still needed ongoing research. Further studies are required to more accurately interpret the molecular mechanism of TGF-β superfamily sex determination genes.

Genome-wide identification and expressional profile of the Dmrt gene family in the swimming crab (Portunus trituberculatus)

The swimming crab, Portunus trituberculatus is one of crucial aquaculture crabs with significant differences in growth and economic performance between male and female swimming crabs. Consequently, the culture of female populations presents higher economic value. The doublesex and mab-3 related transcription factor (Dmrt) gene family are known to play crucial role in gonad differentiation and development. However, there is limited information about this gene family in Portunus trituberculatus. In this study, we identified seven members of the Dmrt gene family in P. trituberculatus based on the published transcriptome and genome data and designated as Ptdmrt-1, Ptdoublesex (Ptdsx), Ptidmrt-1, Ptdmrt-11E, Ptidmrt-2, Ptdmrt-99B, and Ptdmrt-3 based on the homology analysis results, respectively. These Ptdmrt genes distributed across 6 chromosomes and were predicted to encode 283 aa, 288 aa, 529 aa, 436 aa, 523 aa, 224 aa, and 435 aa protein precursors, respectively. The expression patterns of these dmrt genes were characterized by qRT-PCR and gonad transcriptome data. The results showed that five members (Ptdmrt-99B, Ptdsx, Ptdmrt-1, Ptdmrt-3, and Ptdmrt-11E) were differentially expressed between the testis and ovary. In addition, their expression patterns from zoea 2 to juvenile 1 were characterized by published transcriptome data and the results showed that they were lowly expressed and did not exhibit notable difference except for Ptdsx during early development. Overall, majority of Ptdmrt genes were involved in gonad differentiation in the swimming crab. Current findings provide a solid foundation for further exploration of the roles of these genes in gonad development and differentiation in P. trituberculatus.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

A delayed and unsynchronized ovary development as revealed by transcriptome of brain and pituitary of Coilia nasus

Coilia nasus is an anadromous fish that has been successfully domesticated in the last decade due to its high economic value. The fish exhibits a delayed ovary development during the reproductive season, despite breeding and selection for five to six offspring. The molecular mechanism of the delayed ovary development is still unknown, so the obstacles have not been removed in the large-scale breeding program. This study aims to investigate the key genes regulating ovarian development by comparing the transcriptomes of ovarian-stage IV and stage II brain/pituitary of Coilia nasus. Ovarian stages were validated by histological sections. A total of 75,097,641 and 66,735,592 high-quality reads were obtained from brain and pituitary transcriptomes, respectively, and alternatively spliced transcripts associated with gonadal development were detected. Compared to ovarian Ⅱ- brain, 515 differentially expressed genes (DEGs) were upregulated and 535 DEGs were downregulated in ovarian Ⅳ- brain, whereas 470 DEGs were upregulated and 483 DEGs were downregulated in ovarian Ⅳ- pituitary compared to ovarian Ⅱ- pituitary. DEGs involved in hormone synthesis and secretion and in the GnRH signaling pathway were screened. Weighted gene co-expression network analysis identified gene co-expression modules that were positively correlated with ovarian phenotypic traits. The hub genes Smad4 and TRPC4 in the modules were co-expressed with DEGs including Kiss1 receptor and JUNB, suggesting that ovarian development is controlled by a hypothalamic-pituitary-gonadal axis. Our results have provided new insights that advance our understanding of the molecular mechanism of C. nasus reproductive functions and will be useful for future breeding.

Towards a comprehensive understanding of the muscle proteome in Schizothorax labiatus: Insights from seasonal variations, metabolic responses, and reproductive signatures in the River Jhelum

Proteomics is a very advanced technique used for defining correlations, compositions and activities of hundreds of proteins from organisms as well as effectively used in identifying particular proteins with varying peptide lengths and amino acid counts. In the present study, an endeavour has been put forth to create muscle proteome expression of snow trout, Schizothorax labiatus. Liquid chromatography-mass spectrometry (LC-MS) using label free quantification (LFQ) technique has extensively been carried out to explore changes in protein metabolism and its composition to discriminate across species, clarify functions and pinpoint protein biomarkers from organisms. In LFQ technique, the abundances of proteins are determined based on the signal intensities of their corresponding peptides in mass spectrometry. The main benefit of using this method is that it doesn't require pre-labelling proteins with isotopic tags, which streamlines the experimental procedure and gets rid of any bias that might have been caused by the labelling process. LFQ techniques frequently offer a wider dynamic range, making it possible to detect and quantify proteins over a broad range of abundances obtained from the complex biological materials including fish muscle. The results of proteomic analysis could provide an insight in understanding about how various proteins are expressed in response to environmental challenges. For proteomic study, two different weight groups of S. labiatus were taken from River Jhelum based on biological, physiological and logistical factors. These groups corresponded to different life stages, such as younger size and adults/brooders in order to capture potential variations in the muscle proteome related to growth and development. The proteomic analysis of S. labiatus depicted that an overall of 220 proteins in male and 228 in female fish of group 1 were noted. However, when male and female S. labiatus were examined based on spectral count and peptide abundance using ProteinLynx Global Software, a total of 10 downregulated and 32 upregulated proteins were found. In group 2 of S. labiatus, a total of 249 proteins in male and 301 in female fish were documented. When the two genders of S. labiatus were likened to one another by LFQ technique, a total of 41 downregulated and 06 upregulated proteins were identified. The variability in the protein numbers between two fish weight groups reflected biological differences, influenced by factors such as age, developmental stages, physiological condition and reproductive activities. During the study, it was observed that S. labiatus exhibited downregulated levels of proteins that were involved in feeding and growth. The contributing factors to this manifestation could be explained by lower feeding and metabolic activity of fish and decreased food availability during winter in River Jhelum. Contrarily, the fish immune response proteins were found to be significantly over-expressed in S. labiatus, indicating that the environment was more likely to undergo increased microbial infection, pollution load and anthropogenic activities. In addition, it was also discovered that there was an upregulated expression of the reproductive proteins in S. labiatus, which could be linked to the fish's pre-spawning time as the fish used in this study was collected in the winter season which is the pre-spawning period of the fish. Therefore, the present study would be useful in obtaining new insights regarding the molecular makeup of species, methods of adaptation and reactions to environmental stresses. This information contributes to our understanding of basic science and may have applications in environmental monitoring, conservation and preservation of fish species.

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

Genome-wide identification and expression analysis of Dmrt genes in bivalves

In recent years, some common themes in the development of sex-specific traits in different animal lineages have started to emerge since the discovery of the Dmrt (doublesex-mab3-related transcription factor gene) genes. Bivalves are characterized by a diversity of sexual systems, including simultaneous hermaphroditism, sequential hermaphroditism, and strict gonochorism. However, to date, no research has focused on the genome-wide characterization and analysis of Dmrt genes in bivalves. In this study, the identification and analysis of Dmrt genes in 15 bivalves were performed using bioinformatics methods. A total of 55 Dmrt genes were retrieved in the studied bivalve genomes. The number of Dmrt genes in different species ranged from 3 to 5. The phylogenetic tree showed that Dmrt genes in bivalves can be subdivided into 5 classes: the Dmrt2-like class, Dmrt3-like class, Dmrt4/5-like class, Dsx-like class, and scallop-specific Dmrt class. The Ka/Ks ratios suggested that all Dmrt classes underwent purifying selection pressure. Furthermore, the spatiotemporal expression of Dmrt genes in four bivalve species suggested that different Dmrt genes may have different functions, and scallop-specific Dmrt genes may play a key role in sex determination/differentiation. In general, this study provides a molecular basis for in-depth examination of the functions of Dmrt genes and phylogenomic analyses in bivalves.

Genomic organization and gene evolution of two warm temperature acclimation proteins (Wap65s) of Micropterus salmoides and their responses to temperature and bacterial/viral infections

Warm temperature acclimation-related 65-kDa proteins (Wap65s) are fish plasma acute-phase glycoproteins homologous to hemopexin with high affinity and clearance for heme. The study characterized Mswap65-1 and Mswap65-2 genes in Micropterus salmoides. Structural analysis showed MsWap65s contained conserved heme-binding sites. MsWap65-1 had a chloride-binding site similar to hemopexin, while MsWap65-2 had an additional calcium-binding site. Phylogenetic and Ka/Ks analysis showed that fish Wap65s were evolutionarily conserved and underwent strong purifying selection. Functional divergence analysis indicated that fish Wap65-2 retained the putative function of ancestral Wap65, while Wap65-1 underwent neofunctional differentiation. QPCR showed Mswap65s were predominantly expressed in liver, but prolonged hyperthermy inhibited Mswap65-2 expression. Mswap65-2 expression was up-regulated in liver and spleen after Nocardia seriolae infection, while Mswap65-1 was down-regulated. MsWap65-2 may be associated with pathogenesis and play potential role in pathogen resistance. LMBV infection resulted in both significant downregulation of Mswap65s were both significantly down-regulated, with differences observed between sexes. We speculated the immune system might suppress expression after viral infection. Exogenous rMsWap65s were prepared, and injection of rMsWap65s alleviated phenylhydrazine-induced hemolysis and inhibited increases in heme, complement C3 and inflammatory symptoms. Our results contribute to an advanced understanding of the functions and mechanisms of MsWap65s in stress resistance.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Gonadal Transcriptome Analysis and Sequence Characterization of Sex-Related Genes in Cranoglanis bouderius

In China, the Cranoglanis bouderius is classified as a national class II-protected animal. The development of C. bouderius populations has been affected by a variety of factors over the past few decades, with severe declines occurring. Considering the likelihood of continued population declines of the C. bouderius in the future, it is critical to investigate the currently unknown characteristics of gonadal differentiation and sex-related genes for C. bouderius conservation. In this study, the Illumina sequencing platform was used to sequence the gonadal transcriptome of the C. bouderius to identify the pathways and genes related to gonadal development and analyze the expression differences in the gonads. A total of 12,002 DEGs were identified, with 7220 being significantly expressed in the ovary and 4782 being significantly expressed in the testis. According to the functional enrichment results, the cell cycle, RNA transport, apoptosis, Wnt signaling pathway, p53 signaling pathway, and prolactin signaling pathway play important roles in sex development in the C. bouderius. Furthermore, the sequence characterization and evolutionary analysis revealed that AMH, DAX1, NANOS1, and AR of the C. bouderius are highly conserved. Specifically, the qRT-PCR results from various tissues showed significant differences in AMH, DAX1, NANOS1, and AR expression levels in the gonads of both sexes of C. bouderius. These analyses indicated that AMH, DAX1, NANOS1, and AR may play important roles in the differentiation and development of C. bouderius gonads. To our best knowledge, this study is the first to analyze the C. bouderius gonadal transcriptome and identify the structures of sex-related genes, laying the foundation for future research.

Full-length gonad transcriptome analysis of Amur sturgeon Dmrt family genes: identification, characterization, and expression patterns during gonadal differentiation

The regulatory mechanisms that govern sex differentiation in sturgeon are still poorly understood. The doublesex and Mab-3-related transcription factor (Dmrt) gene family is known for its extensive roles in sex determination and differentiation across vertebrates. This study aimed to identify new members of sturgeon Dmrt family genes and core actors in the gonadal differentiation of Amur sturgeon. A full-length gonad transcriptome database was exploited to identify Dmrt gene orthologs. Analyses of phylogenetic relationships and selection pressure were performed, and tissue expression profiles and spatiotemporal expression patterns in gonads were then analyzed using real-time PCR. In total, five Dmrt family genes were identified from the full-length gonad transcriptome, including Dmrt2, DmrtA1, DmrtA2, DmrtB1a, and DmrtB1b. Phylogenetic analysis showed that these genes were clustered into clades corresponding to the doublesex/Mav-3 (DM) genes of vertebrates. Furthermore, the analysis of evolutionary selective pressure indicated that DmrtB1a and DmrtB1b were subject to positive selection, suggesting the existence of adaptive evolution in sturgeon. The extensive tissue expression profiling of each Dmrt family gene revealed typical characteristics. Remarkably, according to a spatiotemporal expression pattern analysis, in later stages, DmrtB1b expression increased in testes and was significantly higher in testes than in ovaries at 24 months after hatching (M) and 36 M. This study provides a genetic resource of full-length Dmrt family genes and increases the understanding of Dmrt functions in sex differentiation in sturgeon.

Chromosome-Scale, Haplotype-Resolved Genome Assembly of Non-Sex-Reversal Females of Swamp Eel Using High-Fidelity Long Reads and Hi-C Data

The Asian swamp eel (Monopterus albus) is an excellent model species for studying sex change and chromosome evolution. M. albus is also widely reared in East Asia and South-East Asia because of its great nutritional value. The low fecundity of this species (about 300 eggs per fish) greatly hinders fries production and breeding programs. Interestingly, about 3–5% of the eels could remain as females for 3 years and lay more than 3,000 eggs per fish, which are referred to as non-sex-reversal (NSR) females. Here, we presented a new chromosome-level genome assembly of such NSR females using Illumina, HiFi, and Hi-C sequencing technologies. The new assembly (Mal.V2_NSR) is 838.39 Mb in length, and the N50 of the contigs is 49.8 Mb. Compared with the previous assembly obtained using the continuous long-read sequencing technology (Mal.V1_CLR), we found a remarkable increase of continuity in the new assembly Mal.V2_NSR with a 20-times longer contig N50. Chromosomes 2 and 12 were assembled into a single contig, respectively. Meanwhile, two highly contiguous haplotype assemblies were also obtained, with contig N50 being 14.54 and 12.13 Mb, respectively. BUSCO and Merqury analyses indicate completeness and high accuracy of these three assemblies. A comparative genomic analysis revealed substantial structural variations (SVs) between Mal.V2_NSR and Mal.V1_CLR and two phased haplotype assemblies, as well as whole chromosome fusion events when compared with the zig-zag eel. Additionally, our newly obtained assembly provides a genomic view of sex-related genes and a complete landscape of the MHC genes. Therefore, these high-quality genome assemblies would provide great help for future breeding works of the swamp eel, and it is a valuable new reference for genetic and genomic studies of this species.

Genome-Wide Identification, Phylogeny, and Expression Profile of the Dmrt (Doublesex and Mab-3 Related Transcription Factor) Gene Family in Channel Catfish (Ictalurus punctatus)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM domain related to sex determination and differentiation, which has been systematically described in various teleost fish, but less in channel catfish (Ictalurus punctatus), an important global aquaculture species in the US and China. In this study, seven Dmrt genes from channel catfish genome were identified and analyzed using bioinformatics methods. Seven IpDmrt genes were distributed unevenly across five chromosomes. Synteny analysis revealed that Dmrt1, Dmrt2a, Dmrt2b, Dmrt3, Dmrt4, and Dmrt5 were relatively conserved in teleost fish. Tissue distribution analysis showed that IpDmrt1, IpDmrt2b, IpDmrt5, and IpDmrt6 exhibited sexually dimorphic expression patterns and, among them, IpDmrt1 and IpDmrt6 had high expression levels in the testes, while IpDmrt2b and IpDmrt5 had more significant expression levels in the ovaries than in other tissues. After 17β-estradiol treatment, IpDmrt2b and IpDmrt5 were significantly up regulated, while the expression of IpDmrt1 and IpDmrt6 was significantly repressed in XY channel catfish ovaries compared with XX channel catfish ovaries. The present study provides a comprehensive insight into the Dmrt gene family of channel catfish. The results suggest that IpDmrt1 and IpDmrt6 may play an important role in testis differentiation/development, while IpDmrt2b and IpDmrt5 are critical in ovary development in this species.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Discovery of the Dmrt gene family members based on transcriptome analysis in mud crab Scylla paramamosain

Doublesex and mab-3 related transcription factors (Dmrts) play crucial roles in sex determination/differentiation and gonad development. The information on Dmrts and their functions are still scarce in mud crab Scylla paramamosain. In this study, 12 published transcriptome data of S. paramamosain were retrieved, pooled, and assembled. From the assembly, 7 Dmrt gene family members were identified and consisted of Spdmrt-like, Spdmrt-1a, Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdoublesex (Spdsx), and Spidmrt-2. These dmrt genes were predicted to encode 224 aa, 465 aa, 435 aa, 276 aa, 520 aa, 552 aa, and 266 aa protein precursors, respectively. The expression patterns of the dmrt genes were characterized by semi-quantitative PCR. The Spdmrt-like and Spdmrt-1a were exclusively detected in gonads, of which both expression levels in the testis were higher than that in the ovary. The Spdmrt-3, Spdmrt-11E, Spidmrt-1, Spdsx, and Spidmrt-2 were observed in various tissues; all these genes were sexually dimorphic except for dmrt-11E. Specifically, the expression level of Spdmrt-3 and Spidmrt-2 were higher in the testis than that in the ovary. On the contrary, the Spdsx and Spidmrt-1 expression level were higher in ovary than that in testis. The present study's findings provided a fundamental understanding of Dmrt gene family members involving sex determination/differentiation and gonad development in the S. paramamosain.