Genome-Wide Identification and Transcriptome-Based Expression Profiling of the Sox Gene Family in the Nile Tilapia (Oreochromis niloticus)

Identification and expression analysis of Sox family genes in echinoderms

The Sox gene family, a collection of transcription factors widely distributed throughout the animal kingdom, plays a crucial role in numerous developmental processes. Echinoderms occupy a pivotal position in many research fields, such as neuroscience, sex determination and differentiation, and embryonic development. However, to date, no comprehensive study has been conducted to characterize and analyze Sox genes in echinoderms. In the present study, the evolution and expression of Sox family genes across 11 echinoderms were analyzed using bioinformatics methods. The results revealed a total of 70 Sox genes, with counts ranging from 5 to 8 across different echinoderms. Phylogenetic analysis revealed that the identified Sox genes could be categorized into seven distinct classes: the SoxB1 class, SoxB2 class, SoxC class, SoxD class, SoxE class, SoxF class and SoxH class. Notably, the SoxB1, SoxB2, and SoxF genes were ubiquitously present in all the echinoderms studied, which suggests that these genes may be conserved in echinoderms. The spatiotemporal expression patterns observed for Sox genes in the three echinoderms indicated that various Sox members perform distinct functional roles. Notably, SoxB1 is likely involved in echinoderm ovary development, while SoxH may play a crucial role in testis development in starfish and sea cucumber. In general, the present investigation provides a molecular foundation for exploring the Sox gene in echinoderms, providing a valuable resource for future phylogenetic and genomic studies.

Genome-wide identification and expression profiles of sex-related gene families in the Pacific abalone Haliotis discus hannai

In recent years, members of the Dmrt family, TGF-β superfamily and Sox family have been recognized as crucial genes for sex determination/differentiation across diverse animal species. Nevertheless, knowledge regarding the abundance and potential functions of these genes in abalone remains limited. In this study, a total of 5, 10, and 7 members of the Dmrt family, the TGF-β superfamily and the Sox family, respectively, were identified in the Pacific abalone Haliotis discus hannai. Sequence characteristics, phylogenetic relationships and spatiotemporal expression profiles of these genes were investigated. Notably, HdDmrt-04 (Dmrt1/1L-like) emerged as a potential mollusc-specific gene with a preponderance for expression in the testis. Interestingly, none of the TGF-β superfamily members exhibited specific or elevated expression in the gonads, highlighting the need for further investigation into their role in abalone sex differentiation. The Sox proteins in H. discus hannai were categorized into 7 subfamilies: B1, B2, C, D, E, F, and H. Among them, HdSox-07 (SoxH-like) was observed to play a crucial role in testis development, while HdSox-03 (SoxB1-like) and HdSox-04 (SoxC-like) probably cooperate in abalone ovary development. Taken together, the results of the present study suggested that HdDmrt-04 and HdSox-07 can be used as male-specific markers for gonad differentiation in H. discus hannai and imply conservation of their functions across invertebrates and vertebrates. Our findings provide new insights into the evolution and genetic structure of the Dmrt family, the TGF-β superfamily and the Sox family in abalone and pave the way for a deeper understanding of sex differentiation in gastropods.

Hypogonadotropic hypogonadism in male tilapia lacking a functional rln3b gene

Relaxin 3 is a neuropeptide that plays a crucial role in reproductive functions of mammals. Previous studies have confirmed that rln3a plays an important role in the male reproduction of tilapia. To further understand the significance of its paralogous gene rln3b in male fertility, we generated a homozygous mutant line of rln3b in Nile tilapia. Our findings indicated that rln3b mutation delayed spermatogenesis and led to abnormal testes structure. Knocking out rln3b gene resulted in a decrease in sperm count, sperm motility and male fish fertility. TUNEL detection revealed a small amount of apoptosis in the testes of rln3b-/- male fish at 390 days after hatching (dah). RT-qPCR analysis demonstrated that mutation of rln3b gene caused a significant downregulation of steroid synthesis-related genes such as cyp17a1, cyp11b2, germ cell marker gene, Vasa, and gonadal somatic cell marker genes of amh and amhr2. Furthermore, we found a significant down-regulation of hypothalamic-pituitary-gonadal (HPG) axis-related genes, while a significantly up-regulation of the dopamine synthetase gene in the rln3b-/- male fish. Taken together, our data strongly suggested that Rln3b played a crucial role in the fertility of XY tilapia by regulating HPG axis genes.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum)

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

Role of Sox3 in Estradiol-Induced Sex Reversal in Pelodiscus sinensis

The Chinese soft-shelled turtle Pelodiscus sinensis, an economically important species in China, exhibits significant sexual dimorphism. Males are more valuable than females owing to their wider calipash and faster growth. Estradiol (E2)-induced sex reversal is used to achieve all-male breeding of turtles; however, the mechanism of this sex reversal remains unclear. In this study, we characterized the Sox3 gene, whose expression level was high in the gonads and brain and exhibited significant sexual dimorphism in the ovary. During embryonic development, Sox3 was highly expressed at the initiation of ovarian differentiation. E2 and Sox3-RNAi treatment before sexual differentiation led to 1352, 908, 990, 1011, and 975 differentially expressed genes in five developmental stages, respectively, compared with only E2 treatment. The differentially expressed genes were clustered into 20 classes. The continuously downregulated and upregulated genes during gonadal differentiation were categorized into Class 0 (n = 271) and Class 19 (n = 606), respectively. KEGG enrichment analysis showed that Sox3 significantly affected sexual differentiation via the Wnt, TGF-β, and TNF signaling pathways and mRNA surveillance pathway. The expression of genes involved in these signaling pathways, such as Dkk4, Nog, Msi1, and Krt14, changed significantly during gonadal differentiation. In conclusion, the deletion of Sox3 may lead to significant upregulation of the mRNA surveillance pathway and TNF and Ras signaling pathways and downregulation of the Wnt and TGF-β signaling pathways, inhibiting E2-induced sex reversal. These findings suggest that Sox3 may play a certain promoting effect during E2-induced sex reversal in P. sinensis.

Regulatory mechanisms of SoxD transcription factors and their influences on male fertility

Members of the SRY-related box (SOX) subfamily D (SoxD) of transcription factors are well conserved among vertebrate species and play important roles in different stages of male reproductive development. In mammals, the SoxD subfamily contains three members: SOX5, SOX6 and SOX13. Here, we describe their implications in testicular development and spermatogenesis, contributing to fertility. We also cover the mechanisms of action of SoxD transcription factors in gene regulation throughout male development. The specificity of activation of target genes by SoxD members depends, in part, on their post-translational modifications and interactions with other partners. Sperm production in adult males requires the coordination in the regulation of gene expression by different members of the SoxD subfamily of transcription factors in the testis. Specifically, the regulation of genes promoting adequate spermatogenesis by SoxD members is discussed in comparison between species.

Transcriptome characterization of gonadal sex differentiation in Pacific bluefin tuna, Thunnus orientalis (Temminck et Schlegel)

Tunas (genus Thunnus) are one of the most ecologically and commercially important fish worldwide. To establish a biological basis for reproduction in this globally essential species, we have recently studied crucial reproductive aspects of the Pacific bluefin tuna (T. orientalis; PBT), as a model of tuna species, based on our closed-cycle aquaculture technology. In this study, we clarified the global expression profile of the genes regulating gonadal sex differentiation in PBT, as this developmental process is vital to sexual reproduction. Based on the results of our comparative (RNA-sequencing) and temporal (qRT-PCR) transcriptome analyses using the updated genome dataset, we propose the molecular mechanisms of gonadal sex differentiation in PBT. In female gonads, foxl2 and cyp19a1a (coding aromatase) are expressed at the onset of sex differentiation. Active aromatase-mediated estrogen biosynthesis, which includes positive regulation of cyp19a1a expression by Foxl2, induces ovarian differentiation. By contrast, dmrt1 and gsdf are upregulated in differentiating male gonads lacking active estrogen synthesis. Dmrt1 and Gsdf would mainly promote testicular differentiation. Furthermore, androgen biosynthesis is upregulated in differentiating male gonad. Endogenous androgens may also be vital to testicular differentiation. This study provides the first comprehensive data clarifying the molecular basis for gonadal sex differentiation in tunas.

Genome-Wide Identification, Evolutionary and Mutational Analysis of the Buffalo Sox Gene Family

The Sox gene family constitutes transcription factors with a conserved high mobility group box (HMG) that regulate a variety of developmental processes, including sex differentiation, neural, cartilage, and early embryonic development. In this study, we systematically analyzed and characterized the 20 Sox genes from the whole buffalo genome, using comparative genomic and evolutionary analyses. All the buffalo Sox genes were divided into nine sub-groups, and each gene had a specific number of exons and introns, which contributed to different gene structures. Molecular phylogeny revealed more sequence similarity of buffalo Sox genes with those of cattle. Furthermore, evolutionary analysis revealed that the HMG domain remained conserved in the all members of the Sox gene family. Similarly, all the genes are under strong purifying selection pressure; seven segmental duplications occurred from 9.65 to 21.41 million years ago (MYA), and four potential recombination breakpoints were also predicted. Mutational analysis revealed twenty non-synonymous mutations with potential effects on physiological functions, including embryonic development and cell differentiation in the buffalo. The present study provides insights into the genetic architecture of the Sox gene family in buffalo, highlights the significance of mutations, and provides their potential utility for marker-assisted selection for targeted genetic improvement in buffalo.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Mature Gonads in Spinibarbus hollandi

Spinibarbus hollandi is an important commercial aquaculture species in southeastern China, but with long maturity period and low egg laying amount. However, there has been little study of its gonad development and reproductive regulation, which limits aquaculture production. Here, for the first time, gonadal transcriptomes of male and female S. hollandi were analyzed. A total of 167,152 unigenes were assembled, with only 48,275 annotated successfully. After comparison, a total of 21,903 differentially expressed genes were identified between male and female gonads, of which 16,395 were upregulated and 5508 were downregulated in the testis. In addition, a large number of differentially expressed genes participating in reproduction, gonad formation and differentiation, and gametogenesis were screened out and the differential expression profiles of partial genes were further validated using quantitative real-time PCR. These results will provide basic information for further research on gonad differentiation and development in S. hollandi.

Characterization of aquaporin-1ab (Aqp1ab) mRNA in mud loach (Misgurnus mizolepis) exposed to heavy metal and immunostimulant stimuli

Aquaporins (AQPs) are key proteins that regulate fluid homeostasis in cells via modulating osmotic water transport. In the present study, we identified three variants of Aqp1ab transcript (mmAQP1ab x1, mmAQP1ab x2, and mmAQP1ab x3) in mud loaches (Misgurnus mizolepis), and their expression patterns were examined in response to heavy metal and immunostimulant exposure. Mud loach Aqp1ab gene has a somewhat different organizational structure (i.e. five exons interrupted by four introns) compared to most other teleostean Aqp1ab orthologues, which have four exons. The 5'-flanking regulatory region of Aqp gene showed diverse transcription factor binding motifs, particularly those associated with stress/immune responses. Developmental expression patterns indicated that Aqp1ab mRNA was maternally inherited, presumably important for fine-tuning gene expression during embryonic and early larval developments. Expression of mud loach Aqp1ab mRNA was significantly and differentially modulated in several tissues (intestine, kidneys, spleen, and liver) in response to various heavy metal treatments. In addition, Aqp1ab gene expression was highly induced in response to immune challenge (LPS and polyI:C injections). Collectively, our results suggested that AQPs are multifunctional effectors playing diverse roles in cellular pathways relevant to immune and/or stress adaptation responses, in addition to their involvement in osmoregulation.

sox19 regulates ovarian steroidogenesis in common carp

In teleost, ovarian steroidogenesis governed by the neuroendocrine system is also regulated by several transcription factors of gonadal origin. Investigating the synchronized interactions between the transcriptional and the hormonal factors is vital to comprehend the mechanisms that lead to gonadal differentiation. This study signifies the role of sry-related box (sox) 19 in ovarian steroidogenesis regulation of the common carp, Cyprinus carpio. Analysis of tissue distribution displayed higher sox19 expression in brain and ovary, and gonadal ontogeny showed higher expression of sox19 at 80 days post hatch (dph). Higher sox19 mRNA expression during spawning and increase of sox19 post human chorionic gonadotropin induction substantiate gonadotropin dependency. Estradiol-17β treatment but not 17α-methyl-di-hydroxy-testosterone to 50 dph common carp for inducing mono-sex, elevated sox19 expression substantially. Sox19 protein was observed in granulosa cells of the follicular layer in common carp ovary. Higher sox19 expression was detected in isolated granulosa and theca cells, in vitro. Transient gene silencing with sox19-siRNA caused downregulation of various ovary-related genes including those specific to activator protein-1 factors, fibroblast growth factors, wnt-signaling, steroidogenic genes along with certain transcription factors. Serum 17α, 20β-dihydroxy-4-pregnen-3-one and estradiol-17β reduced significantly post sox19 silencing, in vivo. Concomitantly, a decrease in aromatase activity was detected post sox19-siRNA treatment, in vivo. This study demonstrates the impact of sox19 in the regulation of common carp ovarian growth and steroidogenesis.

Genome-wide identification and transcriptome-based expression profiling of the Sox gene family in the spinyhead croaker (Collichthys lucidus)

Sox genes encode transcription factors with a high-mobility group (HMG) box, playing critical roles in the initiation and maintenance of a variety of developmental processes, such as sex determination and differentiation. In the present study, we identified 26 Sox genes in the genome of spinyhead croaker Collichthys lucidus (Richardson, 1844) with homology-based analysis of the HMG box. The transcriptome-based expression profiles revealed that the expression of the Sox gene in gonads began to differ between sexes when the body length was 2.74 ± 0.24 cm. At that time, three Sox genes (Sox11b, Sox8a and Sox19) were significantly upregulated, accompanied by the downregulation of 12 Sox genes in the ovary, and six Sox genes were temporarily significantly upregulated in the testis. Afterwards, the expression profile of Sox genes changed only with a small amplitude in both the ovary and testis. For adult tissues, huge differences were observed in the expression profiles of Sox genes between ovaries and testes, as well as small differences in somatic tissues between sexes. These results provide clues to further decipher the role of Sox genes in the processes of sex determination and differentiation in spinyhead croaker and other teleosts.

The transcription factor Sox30 is involved in Nile tilapia spermatogenesis

Spermatogenesis is a complex process in which spermatogonial stem cells differentiate and develop into mature spermatozoa. The transcriptional regulatory network involved in fish spermatogenesis remains poorly understood. Here, we demonstrate in Nile tilapia that the Sox transcription factor family member Sox30 is specifically expressed in the testes and mainly localizes to spermatocytes and spermatids. CRISPR/Cas9-mediated sox30 mutation results in abnormal spermiogenesis, reduction of sperm motility, and male subfertility. Comparative transcriptome analysis shows that sox30 mutation alters the expression of genes involved in spermatogenesis. Further chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), ChIP-PCR, and luciferase reporter assays reveal that Sox30 positively regulates the transcription of ift140 and ptprb, two genes involved in spermiogenesis, by directly binding to their promoters. Taken together, our data indicate that Sox30 plays essential roles in Nile tilapia spermatogenesis by directly regulating the transcription of the spermiogenesis-related genes ift140 and ptprb.

A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal

The yellowfin seabream Acanthopagrus latus is the economically most important Sparidae fish in the northern South China Sea. As euryhaline fish, they are perfect model for investigating osmoregulatory mechanisms in teleosts. Moreover, the reproductive biology of hermaphrodites has long been intriguing; however, little information is known about the molecular pathways underlying their sex change. Here, we report a chromosome level reference genome of A. latus generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The draft genome of yellowfin seabream was 806 Mb, with 732 Mb scaffolds anchored on 24 chromosomes. The contig N50 and scaffold N50 were 2.6 Mb and 30.17 Mb, respectively. The assembly is of high integrity and includes 92.23% universal single-copy orthologues based on benchmarking universal single-copy orthologs (BUSCO) analysis. A total of 19,631 protein-coding genes were functionally annotated in the reference genome. Moreover, ARRDC3 and GSTA gene families which related to osmoregulation underwent an extensive expansion in two euryhaline sparids fish genomes compared to other teleost genomes. Moreover, integrating sex-specific transcriptome analyses, several genes related to the transforming growth factor beta (TGF-β) signalling pathway involved in sex differentiation and development. This genomic resource will not only be valuable for studying the osmoregulatory mechanisms in estuarine fish and sex determination in hermaphrodite vertebrate species, but also provide useful genomic tools for facilitating breeding of the yellowfin seabream.

Comparison of Gonadal Transcriptomes Uncovers Reproduction-Related Genes with Sexually Dimorphic Expression Patterns in Diodon hystrix

Diodon hystrix is a new and emerging aquaculture species in south China. However, due to the lack of understanding of reproductive regulation, the management of breeding and reproduction under captivity remains a barrier for the commercial aquaculture of D. hystrix. More genetic information is needed to identify genes critical for gonadal development. Here, the first gonadal transcriptomes of D. hystrix were analyzed and 151.89 million clean reads were generated. All reads were assembled into 57,077 unigenes, and 24,574 could be annotated. By comparing the gonad transcriptomes, 11,487 differentially expressed genes were obtained, of which 4599 were upregulated and 6888 were downregulated in the ovaries. Using enrichment analyses, many functional pathways were found to be associated with reproduction regulation. A set of sex-biased genes putatively involved in gonad development and gametogenesis were identified and their sexually dimorphic expression patterns were characterized. The detailed transcriptomic data provide a useful resource for further research on D. hystrix reproductive manipulation.

Genome wide identification, phylogeny, and synteny analysis of sox gene family in common carp (Cyprinus carpio)

•

27 SOX (high-mobility group HMG-box) genes were identified in the C. carp genome.

•

SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative protein series from 307 to 509 amino acids.

•

Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %.

•

Chromosomal location and synteny analysis display all SOX gene are located on different chromosomes.

Common carp (Cyprinus carpio) is a commercial fish species valuable for nutritious components and plays a vital role in human healthy nutrition. The SOX (SRY-related genes systematically characterized by a high-mobility group HMG-box) encoded important gene regulatory proteins, a family of transcription factors found in a broad range of animal taxa and extensively known for its contribution in multiple developmental processes including contribution in sex determination across phyla. In our current study, we initially accomplished a genome-wide analysis to report the SOX gene family in common carp fish based on available genomic sequences of zebrafish retrieved from gene repository databases, we focused on the global identification of the Sox gene family in Common carp among wide range of vertebrates and teleosts based on bioinformatics tools and techniques and explore the evolutionary relationships. In our results, a total of 27 SOX (high-mobility group HMG-box) domain genes were identified in the C. carp genome. The full length sequences of SOX genes ranging from 3496 (SOX6) to 924bp (SOX17b) which coded with putative proteins series from 307 to 509 amino acids and all gene having exon number expect SOX9 and SOX13. All the SOX proteins contained at least one conserved DNA-binding HMG-box domain and two (SOX7 and SOX18) were found C terminal. The Gene ontology revealed SOX proteins maximum involvement is in metabolic process 49.796 %, average in biological regulation 45.188 %, biosynthetic process (19.992 %), regulation of cellular process 39.68, 45.508 % organic substance metabolic process, multicellular organismal process 23.23 %,developmental process 21.74 %, system development 16.59 %, gene expression 16.05 % and 14.337 % of RNA metabolic process. Chromosomal location and syntanic analysis show all SOX gene are located on different chromosomes and apparently does not fallow the unique pattern. The maximum linkage of chromosome is (2) on Unplaced Scaffold region. Finally, our results provide important genomic suggestion for upcoming studies of biochemical, physiological, and phylogenetic understanding on SOX genes among teleost.

Systematic identification and expression analysis of the Sox gene family in spotted sea bass (Lateolabrax maculatus)

The Sox gene family encodes a set of transcription factors characterized by a conserved Sry-related high mobility group (HMG)-box domain, which performs a series of essential biological functions in diverse tissues and developmental processes. In this study, the Sox gene family was systematically characterized in spotted sea bass (Lateolabrax maculatus). A total of 26 Sox genes were identified and classified into eight subfamilies, namely, SoxB1, SoxB2, SoxC, SoxD, SoxE, SoxF, SoxH and SoxK. The phylogenetic relationship, exon-intron and domain structure analyses supported their annotation and classification. Comparison of gene copy numbers and chromosome locations among different species indicated that except tandem duplicated paralogs of Sox17/Sox32, duplicated Sox genes in spotted sea bass were generated from teleost-specific whole genome duplication during evolution. In addition, qRT-PCR was performed to detect the expression profiles of Sox genes during development and adulthood. The results showed that the expression of 16 out of 26 Sox genes was induced dramatically at different starting points after the multicellular stage, which is consistent with embryogenesis. At the early stage of sex differentiation, 9 Sox genes exhibited sexually dimorphic expression patterns, among which Sox3, Sox19 and Sox6b showed the most significant ovary-biased expression. Moreover, the distinct expression pattern of Sox genes was observed in different adult tissues. Our results provide a fundamental resource for further investigating the functions of Sox genes in embryonic processes, sex determination and differentiation as well as controlling the homeostasis of adult tissues in spotted sea bass.

A comparative genomic database of skeletogenesis genes: from fish to mammals

Skeletogenesis is a complex process that requires a rigorous control at multiple levels during osteogenesis, such as signaling pathways and transcription factors. The skeleton among vertebrates is a highly conserved organ system, but teleost fish and mammals have evolved unique traits or have lost particular skeletal elements in each lineage. In present study, we constructed a skeletogenesis database containing 4101, 3715, 2996, 3300, 3719 and 3737 genes in Danio rerio, Oryzias latipes, Gallus gallus, Xenopus tropicalis, Mus musculus and Homo sapiens genome, respectively. Then, we found over 55% of the genes are conserved in the six species. Notably, there are 181 specific-genes in the human genome without orthologues in the other five genomes, such as the ZNF family (ZNF100, ZNF101, ZNF14, CALML6, CCL4L2, ZIM2, HSPA6, etc); and 31 genes are identified explicitly in fish species, which are mainly involved in TGF-beta, Wnt, MAPK, Calcium signaling pathways, such as bmp16, bmpr2a, eif4e1c, wnt2ba, etc. Particularly, there are 20 zebrafish-specific genes (calm3a, si:dkey-25li10, drd1a, drd7, etc) and one medaka-specific gene (c-myc17) that may alter skeletogenesis formation in the corresponding species. The database provides the new systematic genomic insights into skeletal development from teleosts to mammals, which may help to explain some of the complexities of skeletal phenotypes among different vertebrates and provide a reference for the treatment of skeletal diseases as well as for applications in the aquaculture industry.

Role of sox30 in regulating testicular steroidogenesis of common carp

Expression of transcription factors is crucial for the regulation of steroidogenesis and gonadal development in fish. SRY-related box (SOX) proteins regulate gene expression of various events related to vertebrate reproduction. This study reports the role of sox30 and its influence on sox9a/b in regulating testicular steroidogenesis of the common carp, Cyprinus carpio. Tissue distribution showed predominant expression of sox30 in gonads, while gonadal ontogeny indicated significant dimorphic expression of sox30 from 120 days post hatch. Higher sox30 transcripts during the spawning season, an elevation of sox30 after human chorionic gonadotropin induction, and 11-ketotestosterone (11-KT) treatment authenticate gonadotropin dependency. Treatment of 17α-methyl-di-hydroxy-testosterone to juvenile common carp for mono-sex induction, vis-à-vis elevated sox30 expression. Sox30 protein was detected abundantly in spermatocytes and spermatid/sperm of carp testis. Transient silencing of sox30 using small interfering RNAs decreased sox9a/b expression, lead to downregulation of certain molecule/factor, transcription factor, germ/stem cell marker, and steroidogenesis-related enzyme genes. Serum testosterone and 11-KT decreased significantly upon transient silencing of sox30, in vivo. Concomitantly, a reduction in testicular microsomal 11-β hydroxysteroid dehydrogenase activity was observed. These results demonstrate the influence of sox30 as well as sox9a/b in the regulation of testicular steroidogenesis in common carp.

Dll4 Blockade Promotes Angiogenesis in Nonhealing Wounds of Sox7-Deficient Mice

Objective: This study aimed to elucidate the role of the proangiogenic transcription factors Sox7 and Sox17 in the wound healing process and investigate the therapeutic potential of Dll4 blockade, which is an upstream regulator of Sox17, for the treatment of nonhealing wounds.

Approach: After generating a full-thickness skin defect wound model of endothelial Sox7- and/or Sox17-deficient mice, we measured the wound healing rates and performed histological analysis. The effects of an anti-Dll4 antibody on wound angiogenesis in Sox7-deficient mice and db/db diabetic mice were assessed.

Results: Sox7 and/or Sox17 deletion delayed wound healing. Moreover, the loss of Sox7 and Sox17 inhibited wound angiogenesis, without affecting the expression of the other. Of interest, after anti-Dll4 antibody treatment, Sox17 levels were increased and the suppression of angiogenesis was alleviated in Sox7-deficient mice and db/db diabetic mice. Consequently, Dll4 blockade effectively recovered the observed delay in wound healing.

Innovation: The proangiogenic role of Sox7 and Sox17 in wound angiogenesis was addressed and effective treatment of nonhealing wounds by Dll4 blockade was suggested.

Conclusion: This study revealed the proangiogenic role of the transcription factors Sox7 and Sox17 in wound angiogenesis. Furthermore, we suggest a novel method for treating nonhealing wounds by particularly targeting the Dll4–Sox17 axis.

Transcriptome analysis of blood for the discovery of sex-related genes in ricefield eel Monopterus albus

The blood acts as a transfer channel for a variety of factors in the whole body. The ricefield eel (Monopterus albus) is a protogynous hermaphrodite vertebrate. Until now, no research has reported an analysis of the blood transcriptome during the process of sexual development in the ricefield eel. In this study, the transcriptome sequencing of blood samples from male and female ricefield eels was completed with a total of 34.70 Gb clean data. The clean data of each sample all reached 5.23 GB, and the percent of the Q30 basic group was 88.62% and above. A total of 106,369 unigenes were obtained after assembly, including 13,296 unigenes with a length of more than 1 kb. Further functional annotation analysis showed that there are 28,522 unigenes that can be annotated. The annotations of genes with differential expression revealed that there were 563 genes with significant differential expression in the blood of male and female ricefield eels, including 91 upregulated genes and 472 downregulated genes. Among which, 14 genes may be closely related to sex differentiation, the qPCR was used to confirmed the expression pattern of those genes and result shown that 11 genes were downregulated and 3 genes were upregulated, consistent with the results of our RNA-Seq analysis. This blood transcript dataset will open future research avenues on ricefield eel sex development and differentiation.

Sox Gene Family Revealed Genetic Variations in Autotetraploid Carassius auratus

The Sox gene family encoded transcription factors that played key roles in developmental processes in vertebrates. To further understand the evolutionary fate of the Sox gene family in teleosts, the Sox genes were comprehensively characterized in fish of different ploidy levels, including blunt snout bream (2n = 48, Megalobrama amblycephala, BSB), goldfish (2n = 100, Carassius auratus red var., 2nRCC), and autotetraploid C. auratus (4n = 200, 4nRCC). The 4nRCC, which derived from the whole genome duplication (WGD) of 2nRCC, were obtained through the distant hybridization of 2nRCC (♀) × BSB (♂). Compared with the 26 Sox genes in zebrafish (2n = 50, Danio rerio), 26, 47, and 92 putative Sox genes were identified in the BSB, 2nRCC, and 4nRCC genomes, respectively, and classified into seven subfamilies (B1, B2, C, D, E, F, and K). Comparative analyses showed that 89.36% (42/47) of Sox genes were duplicated in 2nRCC compared with those in BSB, while 97.83% (90/92) of Sox genes were duplicated in 4nRCC compared with those in 2nRCC, meaning the Sox gene family had undergone an expansion in BSB, 2nRCC, and 4nRCC, respectively, following polyploidization events. In addition, potential gene loss, genetic variations, and paternal parent SNP locus insertion occurred during the polyploidization events. Our data provided new insights into the evolution of the Sox gene family in polyploid vertebrates after several rounds of WGD events.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Transcriptome Profiling and Analysis of Genes Associated with High Temperature-Induced Masculinization in Sex-Undifferentiated Nile Tilapia Gonad

Artificially high temperatures during critical thermosensitive periods (TSPs) can induce the sex reversal of Nile tilapia (Oreochromis niloticus) females into pseudomales; Nile tilapia is a GSD + TE (genotypic plus temperature effects) fish species. Previous studies have shown that water temperature affects the expression levels of many genes in the gonad or brain in various teleost species. However, few studies on the effect of temperature at the whole-gonad transcriptomic level in the early stage of sex differentiation have been reported in fish species exhibiting GSD + TE. In this study, RNA-Seq was performed to characterize the transcriptomic profile and identify genes exhibiting temperature- and sex-biased expressions in the Nile tilapia gonad at 21 dpf. A total of 42 genes were found to be associated with both high-temperature treatment and sex development, as the expression levels of these genes differed in both FC (female control) vs MC (male control) and FC vs FT (high temperature-treated females in the TSP). Among these genes, the transcriptional alterations of many male sex determination and differentiation genes, such as Dmrt1, Gsdf, and the DNA damage-inducible protein GADD45 alpha, suggested that the male pathway is initiated after high-temperature treatment and that its initiation may play a role in high temperature-induced masculinization in Nile tilapia. The qRT-PCR validation results for thirteen differentially expressed genes showed that the Pearson's correlation of the log₁₀ fold change values between the qPCR and RNA-Seq results was 0.70 (p < 0.001), indicating the accuracy and reliability of the RNA-Seq results. Our study provides insights into how high-temperature treatment induces the sex reversal of Nile tilapia females.

Identification and functional analysis of SOX transcription factors in the genome of the Chinese soft-shell turtle (Pelodiscus sinensis)

SOX transcription factors play an irreplaceable role in biological developmental processes. Sox genes have been identified in a wide variety of species; however, their identification and functional analysis in the genome of the Chinese soft-shell turtle (Pelodiscus sinensis) have not been performed. In the present study, the Chinese soft-shell turtle genome was found to contain 17 Sox genes, which were categorized into seven groups according to their phylogenetic relationships. Gene structure and protein motif analysis of the Sox genes showed that within the same phylogenetic group, their exon-intron number and motif structure of the Sox family were relatively conserved, but diverged in the comparison between different groups. Sexual dimorphism expression analysis for the Sox genes displayed that Sox8 and Sox9 were upregulated in the testis, while Sox3, Sox7, Sox11, and Sox13 were upregulated in the ovary. A correlation network analysis of SOX transcription factors with their target genes analysis showed that Sox3 correlated negatively with Sox9 and gata4. Sox11 and Sox7 correlated negatively with gata4. Sox8 and Sox9 correlated positively with gata4. Therefore, the genome-wide identification and functional analysis of the Sox gene family will be useful to further reveal the functions of Sox genes in the Chinese soft-shell turtle.

Beta-Hydroxysteroid Dehydrogenase Genes in Orange-Spotted Grouper (Epinephelus coioides): Genome-Wide Identification and Expression Analysis During Sex Reversal

Beta-hydroxysteroid dehydrogenases (β-HSDs) are a group of steroidogenic enzymes that are involved in steroid biosynthesis and metabolism, and play a crucial role in mammalian physiology and development, including sex determination and differentiation. In the present study, a genome-wide analysis identified the numbers of β-hsd genes in orange-spotted grouper (Epinephelus coioides) (19), human (Homo sapiens) (22), mouse (Mus musculus) (24), chicken (Gallus gallus) (16), xenopus (Xenopus tropicalis) (24), coelacanth (Latimeria chalumnae) (17), spotted gar (Lepisosteus oculatus) (14), zebrafish (Danio rerio) (19), fugu (Takifugu rubripes) (19), tilapia (Oreochromis niloticus) (19), medaka (Oryzias latipes) (19), stickleback (Gasterosteus aculeatus) (17) and common carp (Cyprinus carpio) (27) samples. A comparative analysis revealed that the number of β-hsd genes in teleost fish was no greater than in tetrapods due to gene loss followed by a teleost-specific whole-genome duplication event. Based on transcriptome data from grouper brain and gonad samples during sex reversal, six β-hsd genes had relatively high expression levels in the brain, indicating that these genes may be required for neurogenesis or the maintenance of specific biological processes in the brain. In the gonad, two and eight β-hsd genes were up- and downregulated, respectively, indicating their important roles in sex reversal. Our results demonstrated that β-hsd genes may be involved in the sex reversal of grouper by regulating the synthesis and metabolism of sex steroid hormones.

Genome-Wide Identification, Evolution and Expression of the Complete Set of Cytoplasmic Ribosomal Protein Genes in Nile Tilapia

Ribosomal proteins (RPs) are indispensable in ribosome biogenesis and protein synthesis, and play a crucial role in diverse developmental processes. In the present study, we carried out a comprehensive analysis of RPs in chordates and examined the expression profiles of the complete set of 92 cytoplasmic RP genes in Nile tilapia. The RP genes were randomly distributed throughout the tilapia genome. Phylogenetic and syntenic analyses revealed the existence of duplicated RP genes from 2R (RPL3, RPL7, RPL22 and RPS27) and 3R (RPL5, RPL19, RPL22, RPL41, RPLP2, RPS17, RPS19 and RPS27) in tilapia and even more from 4R in common carp and Atlantic salmon. The RP genes were found to be expressed in all tissues examined, but their expression levels differed among different tissues. Gonadal transcriptome analysis revealed that almost all RP genes were highly expressed, and their expression levels were highly variable between ovaries and testes at different developmental stages in tilapia. No sex- and stage-specific RP genes were found. Eleven RP genes displayed sexually dimorphic expression with nine higher in XY gonad and two higher in XX gonad at all stages examined, which were proved to be phenotypic sex dependent. Quantitative real-time PCR and immunohistochemistry ofRPL5b and RPL24 were performed to validate the transcriptome data. The genomic resources and expression data obtained in this study will contribute to a better understanding of RPs evolution and functions in chordates.

Genome-wide analysis of the Chinese sturgeon sox gene family: identification, characterisation and expression profiles of different tissues

The sox family is assumed to be responsible for a number of developmental systems. Genome sequencing technology makes it possible to scan sox genes and conduct characteristic analyses of different species. In fish, full characterisation of sox genes at the genome-wide level has been reported for pufferfish Takifugu rubripes, medaka Oryzias latipes, tilapia Oreochromis niloticus and channel catfish Ictalurus punctatus. However, no systematic investigation of the sox family in sturgeons (Acipenseridae) has been reported to date. This study conducted genome-wide identification of the sox genes in the Chinese sturgeon Acipenser sinensis and profiled their tissue distribution between male and female individuals. In total, 19 sox genes were identified, including soxb1, b2, c, d, e, f and h, in the Chinese sturgeon. Genomic structure analysis indicated relatively conserved exon-intron structures in each sox group and phylogenetic analysis supported the previous classification of the sox family. Most of the sox genes showed a tissue-specific expression pattern, indicating the possible involvement of Chinese sturgeon sox genes at different developmental processes such as cardiac and gonadal development. This study provides a comprehensive resource of Chinese sturgeon sox genes and enables a better understanding of the evolution and function of the sox family.

Comparative transcriptome analysis of male and female gonads reveals sex-biased genes in spotted scat (Scatophagus argus)

Scatophagus argus is a new emerging aquaculture fish in East and Southeast Asia. To date, research on reproductive development and regulation in S. argus is lacking. Additionally, genetic and genomic information about reproduction, such as gonadal transcriptome data, is also lacking. Herein, we report the first gonadal transcriptomes of S. argus and identify genes potentially involved in reproduction and gonadal development. A total of 136,561 unigenes were obtained by sequencing of testes (n = 3) and ovaries (n = 3) at stage III. Genes upregulated in males and females known to be involved in gonadal development and gametogenesis were identified, including male-biased dmrt1, amh, gsdf, wt1a, sox9b, and nanos2, and female-biased foxl2, gdf9, bmp15, sox3, zar1, and figla. Serum estradiol-17β and 11-ketotestosterone levels were biased in female and male fish, respectively. Sexual dimorphism of serum steroid hormone levels were interpreted after expression analysis of 20 steroidogenesis-related genes, including cyp19a1a and cyp11b2. This gonadal transcript dataset will help investigate functional genes related to reproduction in S. argus.

Transcription of the Sox30 Gene Is Positively Regulated by Dmrt1 in Nile Tilapia

The Sox family member Sox30 is highly expressed in the testis of several vertebrate species and has been shown to play key roles in spermiogenesis. However, its transcription regulation remains unclear. Here, we analyzed the Sox30 promoter from the teleost fish Nile tilapia (Oreochromis niloticus) and predicted a putative cis-regulatory element (CRE) for doublesex and mab-3 related transcription factor 1 (Dmrt1), a male-specific transcription factor involved in male sex differentiation. Transcriptional profiling revealed that Sox30 and Dmrt1 similarly exhibited a high expression in tilapia testes from 90 days after hatching (dah) to 300 dah, and the transcription of the Sox30 gene was reduced about one-fold in the testes of male tilapia with Dmrt1 knockdown. Further dual-luciferase reporter assay confirmed that Dmrt1 overexpression significantly promoted transcriptional activity of the Sox30 promoter and this promotion was decreased following the mutation of putative CRE for Dmrt1 within the Sox30 promoter. Chromatin immunoprecipitation-based PCR (ChIP-PCR) and electrophoretic mobility shift assay (EMSA) demonstrated that Dmrt1 directly binds to putative CRE within the Sox30 promoter. These results together indicate that Dmrt1 positively regulates the transcription of the tilapia Sox30 gene by directly binding to specific CRE within the Sox30 promoter.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Phenotypic Stability of Sex and Expression of Sex Identification Markers in the Adult Yesso Scallop Mizuhopecten yessoensis throughout the Reproductive Cycle

Simple Summary

Bivalve sex is thought to fluctuate depending on environmental conditions. So far, there has been no investigation on the phenotypic stability of sex in the commercially important Yesso scallop Mizuhopecten yessoensis. The present study revealed that the sex of the Yesso scallop is stable after initial sex differentiation and that this species maintains a sex-stable maturation system throughout its life. In addition, gonad differentiation for each sex was precisely characterized by using molecular markers throughout the maturational cycle.

Abstract

The objective of the present study was to analyze the phenotypic stability of sex after sex differentiation in the Yesso scallop, which is a gonochoristic species that has been described as protandrous. So far, no study has investigated in detail the sexual fate of the scallop after completion of sex differentiation, although bivalve species often show annual sex change. In the present study, we performed a tracking experiment to analyze the phenotypic stability of sex in scallops between one and two years of age. We also conducted molecular marker analyses to describe sex differentiation and gonad development. The results of the tracking experiment revealed that all scallops maintained their initial sex phenotype, as identified in the last reproductive period. Using molecular analyses, we characterized my-dmrt2 and my-foxl2 as sex identification markers for the testis and ovary, respectively. We conclude by proposing that the Yesso scallop is a sex-stable bivalve after its initial sex differentiation and that it maintains a sex-stable maturation system throughout its life. The sex-specific molecular markers identified in this study are useful tools to assess the reproductive status of the Yesso scallop.

Profiling expression changes of genes associated with temperature and sex during high temperature-induced masculinization in the Nile tilapia brain

Fish sex-determining mechanisms can be classified as genotypic (GSD), temperature (TSD), or genotypic plus temperature effects (GSD+TE). Previous studies have shown that culturing water temperature during thermosensitive periods (TSP) could affect the expression of many genes in the gonad in some fish. However, few studies have focused on gene expression changes in the brain after temperature treatment during TSP in fish species. In this study, three families were developed by crossing XX neomales with XX females and one of them was used for transcriptome analysis. The results showed that a total of 105, 3164 and 4666 DEGs were respectively obtained in FC (female control) vs. FT (high temperature-treated females at TSP), FC vs. MC (male control), and MC vs. FT comparison groups. By profiling analysis, we show that the mRNA expression levels of 16 differentially expressed genes (DEGs) exhibited significant downregulation or upregulation after high temperature treatment and reached a similar level as that in MC. Among the 16 DEGs, LOC100699848 (lysine specific demethylase 6A) and Jarid2 contained JmjC domain, showing the possible important role of JmjC domain in response to temperature treatment in Nile tilapia. Kdm6b (lysine demethylase 6B) and Jarid2 have been shown to play important roles in reptile TSD, showing the relative conservation of underlying regulation mechanisms between TSD in reptile and TSD or GSD+TE in fish species. Finally, the transcriptome profiling was validated by quantitative real-time PCR in nine selected genes. These results provide a direction for investigating the GSD+TE molecular mechanism in fish species.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.

The evolution and possible role of two Sox8 genes during sex differentiation in Japanese flounder (Paralichthys olivaceus)

Sox8 genes, as members of the Sox family, have been studied widely in mammals. However, regulation of sox8 genes in teleosts has rarely been studied, and functional analysis of these genes in teleosts has rarely been performed. Here, two duplicates of sox8 genes were identified in Japanese flounder, Posox8a and Posox8b. The analysis of expression showed that Posox8a and Posox8b were expressed in Sertoli cells of the testis, indicating that they play important roles in development and functional maintenance of the testis. Positive selection and phylogenetic analysis found that both Posox8a and Posox8b underwent the purification selection during evolutionary and that sox8 was most likely to be the ancestor sox8a. These results suggested that both Posox8a and Posox8b had important biological functions after generation from three rounds of whole-genome duplication in Japanese flounder. The functional differentiation of Posox8a and Posox8b was verified using cell transfection and dual-luciferase reporter assays; Posox8a overexpression-promoted 3β-hydroxysteroid dehydrogenase expression and Posox8b overexpression-promoted cytochrome P450 aromatase (cyp19a1; P450arom) expression. Finally, combined with Posox8a and Posox8b expression analysis from 30 to 100 days after hatch, we speculated that Posox8a and Posox8b might participate in the process of sex differentiation and gonadogenesis by regulating sex hormone biosynthesis in the Japanese flounder. Our study is the first to demonstrate the possible mechanism of Posox8a and Posox8b in Japanese flounder sex differentiation and gonadogenesis, laying a solid foundation for functional studies of sox8 genes in teleosts.

Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

DNA methyltransferases (dnmts) are responsible for DNA methylation and play important roles in organism development. In this study, seven dnmts genes (dnmt1, dnmt2, dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1, dnmt3bb.2) were identified in Nile tilapia. Comprehensive analyses of dnmts were performed using available genome databases from representative animal species. Phylogenetic analysis revealed that the dnmts family were highly conserved in teleosts. Based on transcriptome data from eight adult tilapia tissues, the dnmts were found to be dominantly expressed in the head kidney, testis and ovary. Analyses of the gonadal transcriptome data in different developmental stages revealed that all dnmts were expressed in both ovary and testis, and four de novo dnmts (dnmt3aa, dnmt3ab, dnmt3bb.1, dnmt3bb.2) showed higher expression in the testis than in the ovary. Furthermore, during sex reversal induced by Fadrozole, the expression of these four de novo dnmts increased significantly in treated group compared to female control group. By in situ hybridization, the seven dnmts were found to be expressed mainly in phase I and II oocytes of the ovary and spermatocytes of the testis. When gonads were incubated with a methyltransferase inhibitor (5-AzaCdR) in vitro, the expression of dnmts genes were down-regulated significantly, while the expression of cyp19a1a (a key gene in female pathway) and dmrt1 (a key gene in male pathway) increased significantly. Our results revealed the conservation of dnmts during evolution and indicated a potential role of dnmts in epigenetic regulation of gonadal development.