Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish

Comparative transcriptome analysis of ovaries and testes reveals sex-biased genes and pathways in zebrafish

Zebrafish (Danio rerio) is a widely recognized and extensively studied model organism in scientific research. The regulatory mechanism of gonadal development and differentiation of this species has aroused considerable attention. Nonetheless, the major sex-biased genes and pathways associated with gonadal development remain elusive. Therefore, to comprehend this intricate process, gonadal transcriptome sequencing was carried out to identify differentially expressed genes (DEGs) between the testes and ovaries of adult zebrafish. The preliminary assessment yielded a total of 23,529,272 and 23,521,368 clean reads from the cDNA libraries of ovaries and testes. Afterward, a comparative analysis of the transcriptome revealed 3,604 upregulated and 11,371 downregulated DEGs in the ovaries compared to the testes. Of these genes, 428 were exclusively expressed in females, while 3,516 were exclusively expressed in males. Additionally, further assessments were conducted to explore the functions associated with these DEGs in various biological processes. The data revealed their involvement in sex-biased pathways, such as progesterone-mediated oocyte maturation, oocyte meiosis, cytokine-cytokine receptor interaction, and cardiac muscle contraction. Finally, the expression levels of 14 sex-biased DEGs (cdc20, ccnb1, ypel3, chn1, bmp15, rspo1, tnfsf10, egfra, acta2, cox8a, gsdf, dmrt1, star, and cyp17a1) associated with the enriched pathways were subjected to further validation through qRT-PCR. The data acquired from these investigations offer valuable resources to support further exploration of the mechanisms governing sexual dimorphism and gonadal development in zebrafish.

Effects of exogenous steroid hormones on growth, body color, and gonadal development in the Opsariichthys bidens

To investigate the effects of exogenous steroid hormones on growth, body color, and gonadal development in the Opsariichthys bidens (O. bidens), synthetic methyltestosterone (MT) and 17β-estradiol (E2) were used for 28 days' treatment of 4-month-old O. bidens before the breeding season. Our results suggested that MT had a significant growth-promoting effect (P < 0.05), whereas E2 played an inhibitory role. On the body surface, the females in the MT group showed gray stripes, and the fish in other groups showed no obvious stripes. The males with MT treatment displayed brighter blue-green stripes compared to the CK and E2 groups. The histological analysis showed that the MT significantly promoted testes development in males, blocked oocyte development, and caused massive apoptosis in females, whereas the E2 group promoted ovarian development and inhibited testes development. Based on qRT-PCR analysis, in females, the expression of igf-1, dmrt1, and cyp19a1a genes revealed that E2 treatment resulted in down-regulation of igf-1 expression and up-regulation of cyp19a1a expression. In males, igf-1 and dmrt1 were significantly up-regulated after MT treatment, and E2 treatment led to down-regulation of igf-1. Therefore, this study demonstrates that MT and E2 play an important role in reversing the morphological sex characteristics of females and males.

Identification of Sex-Specific Markers and Candidate Genes Using WGS Sequencing Reveals a ZW-Type Sex-Determination System in the Chinese Soft-Shell Turtle (Pelodiscus sinensis)

Male and female Chinese soft-shelled turtles (Pelodiscus sinensis) have sex-dimorphic growth patterns, and males have higher commercial value because of their larger size and thicker calipash. Thus, developing sex-specific markers is beneficial to studies on all-male breeding in P. sinensis. Here, we developed an accurate and efficient workflow for the screening of sex-specific sequences with ZW or XY sex determination systems. Based on this workflow, female and male P. sinensis reference genomes of 2.23 Gb and 2.26 Gb were obtained using de novo assembly. After aligning and filtering, 4.01 Mb female-specific sequences were finally identified. Subsequently, the seven developed sex-specific primer pairs were 100% accurate in preliminary, population, and embryonic validation. The presence and absence of bands for the primers of P44, P45, P66, P67, P68, and P69, as well as two and one bands for the PB1 primer, indicate that the embryos are genetically female and male, respectively. NR and functional annotations identified several sex-determining candidate genes and related pathways, including Ran, Eif4et, and Crkl genes, and the insulin signaling pathway and the cAMP signaling pathway, respectively. Collectively, our results reveal that a ZW-type sex-determination system is present in P. sinensis and provide novel insights for the screening of sex-specific markers, sex-control breeding, and the studies of the sex determination mechanism of P. sinensis.

Production of all-male non-transgenic zebrafish by conditional primordial germ cell ablation

Many fish species exhibit remarkable sexual dimorphism, with males possessing numerous advantageous traits for commercial production by aquaculture such as faster growth rate, more efficient food energy utilization for muscle development, and better breeding performance. Several studies have shown that a decrease in the number of primordial germ cells (PGCs) during early development leads predominantly to male progeny. In this study, we developed a method to obtain all-male zebrafish (Danio rerio) by targeted PGC ablation using the nitroreductase/metronidazole (NTR/Mtz) system. Embryos generated by female heterozygous Tg(nanos3:nfsB-mCherry-nanos3 3'UTR) and male wild-types (WTs) were treated with vehicle or Mtz. Compared to vehicle-treated controls, 5.0 and 10.0 mM Mtz treatment for 24 h significantly reduced the number of PGCs and yielded an exclusively male phenotype in adulthood. The gonads of offspring treated with 5.0 mM Mtz exhibited relatively normal morphology and histological characteristics. Furthermore, these males were able to chase females, spawn, and produce viable offspring, while about 20.0% of males treated with 10.0 mM Mtz were unable to produce viable offspring. The 5.0 mM Mtz treatment protocol may thus be suitable for large-scale production of fertile male offspring. Moreover, about half of these males were WT as evidenced by the absence of nfsB gene expression. It may thus be possible to breed an all-male WT fish population by Mtz-mediated PGC ablation.

Morphological and transcriptional analysis of sexual differentiation and gonadal development in a burrowing fish, the four-eyed sleeper (Bostrychus sinensis)

Four-eyed sleeper (Bostrychus sinensis) is a commercially important sea water fish, and the male individuals exhibit significant advantages in somatic growth and stress resistance, so developing sex control strategy to create all-male progeny will produce higher economic value. However, little is known about the genetic background associated with sex differentiation in this species. In this study, we investigated gonadal development and uncovered critical window stages of sexual differentiation (about 2 mph), transition from proliferation to differentiation in female germ stem cells (GSCs) (2-3 mph) and male GSCs (3-4 mph). De novo transcriptome analysis revealed candidate genes and signaling pathways associated with sexual differentiation and gonadal development in four-eyed sleeper. The results showed that sox9 and zglp1 were the earliest sex-biased transcription factors during sex differentiation. Down-regulation of chemokine, cytokines-cytokine receptors and up-regulation of cellular senescence pathway might be involved in GSC differentiation. Weighted gene correlation network analysis showed that metabolic pathway and occludin were the hub signaling and gene in ovarian development, meanwhile the MAPK signaling pathways, cellular senescence pathway and ash1l (histone H3-lysine4 N-trimethyltransferase) were the hub pathways and gene in testicular development. The present work elucidated the developmental processes of sexual differentiation and gonadal development and revealed their associated revealed genes and signaling pathways in four-eyed sleeper, providing theoretical basis for developing sex-control techniques.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.

Phosphoproteomics Reveal New Candidates in Abnormal Spermatogenesis of Pseudomales in Cynoglossus semilaevis

Phosphorylation is a post-translational modification that contributes to versatile protein functions in spermatogenesis, and the variations they generate usually results in abnormal spermatogenesis or sperm dysfunction. The sex-reversal phenomenon exists in Chinese tongue sole under certain conditions such that individuals with a ZW genotype can acquire a male phenotype and are thus called pseudomales. Pseudomale tongue sole can reach sexual maturity but produce only Z-type sperm, and the Z sperm carries paternal epigenetic information. Whether phosphorylation plays a role in the sperm abnormality of pseudomales is unknown. In this study, a phosphoproteomic analysis was performed to compare protein phosphorylation profiles between pseudomale and male testes. Altogether, we identified 14,253 phosphopeptides matching with 4843 proteins, with 1329 differentially phosphorylated peptides corresponding to 1045 differentially phosphorylated proteins (DPPs). Phosphorylation at 781 sites was upregulated and at 548 sites was downregulated. Four motifs were identified among differentially phosphorylated peptides, which were "SP", "SD", "RxxS", and "TP". Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that the cell cycle and DNA/RNA processing were significantly enriched with the genes encoding DPPs. To analyze DPP function in depth, a protein-protein interaction network was constructed, and Ran-binding protein 2 was found to play a central role in spermatogenesis by regulating several processes such as the cell cycle, eukaryotic translation, ubiquitination, and minichromosome maintenance. In kinase-associated network analyses, two "mitogen-activated protein kinase (Mapk)-centered" clusters were identified that may account for abnormal spermatogenesis in pseudomales. One cluster was centered on Mapk6, which predominantly regulated the cell cycle by interacting with several cyclin-dependent kinases, and the other was centered on the "testis-expressed kinase 1-like (Tesk1l)/Pim1l-Mapk4l- testis-expressed 14 (Tex14)" kinase cascade, which might contribute to spermatogenesis by regulating β-catenin. Taken together, these data suggested the new candidates involved in pseudomale sperm abnormalities and provided clues to discover the phosphorylated regulatory mechanism underlying tongue sole spermatogenesis.

A chromosome-level genome assembly of Hong Kong catfish (Clarias fuscus) uncovers a sex-determining region

BACKGROUND

Hong Kong catfish (Clarias fuscus) is an ecologically and economically important species that is widely distributed in freshwater regions of southern China. Hong Kong catfish has significant sexual growth dimorphism. The genome assembly of the Hong Kong catfish would facilitate study of the sex determination and evolution mechanism of the species.

RESULTS

The first high-quality chromosome-level genome of the Hong Kong catfish was constructed. The total genome was 933.4 Mb, with 416 contigs and a contig N50 length of 8.52 Mb. Using high-throughput chromosome conformation capture (Hi-C) data, the genome assembly was divided into 28 chromosomes with a scaffold N50 length of 36.68 Mb. A total of 23,345 protein-coding genes were predicted in the genome, and 94.28% of the genes were functionally annotated in public databases. Phylogenetic analysis indicated that C. fuscus and Clarias magur diverged approximately 63.7 million years ago. The comparative genome results showed that a total of 60 unique, 353 expanded and 851 contracted gene families were identified in Hong Kong catfish. A sex-linked quantitative trait locus identified in a previous study was located in a sex-determining region of 30.26 Mb (0.02 to 30.28 Mb) on chromosome 13 (Chr13), the predicted Y chromosome. This QTL region contained 785 genes, of which 18 were identified as sex-related genes.

CONCLUSIONS

This study is the first to report the chromosome-level genome assembly of Hong Kong catfish. The study provides an excellent genetic resource that will facilitate future studies of sex determination mechanisms and evolution in fish.

Inhibition of the NF-κB signaling pathway affects gonadal differentiation and leads to male bias in Paramisgurnus dabryanus

In recent years, sex-controlled breeding has emerged as an effective strategy to enhance the yields of economic animals with different growth characteristics, while increasing the economic benefits of aquaculture. It is known that the NF-κB pathway participates in gonadal differentiation and reproduction. Therefore, we used the large-scale loach as a research model for the present study and selected an effective inhibitor of the NF-κB signaling pathway (QNZ). This, to investigates the impacts of the NF-κB signaling pathway on gonadal differentiation during a critical period of gonad development and after maturation. Simultaneously, the sex ratio bias and the reproductive capacities of adult fish were analyzed. Our results indicated that the inhibition of the NF-κB signaling pathway influenced the expression of genes related to gonad development, regulated the gene expression related to the brain-gonad-liver axis of juvenile loaches, and finally impacted the gonadal differentiation of the large-scale loach and promoted a male-biased sex ratio. Meanwhile, high QNZ concentrations affected the reproductive abilities of adult loaches and inhibited the growth performance of offspring. Thus, our results deepened the exploration of sex control in fish and provided a certain research basis for the sustainable development of the aquaculture industry.

Gonadal transcriptomes reveal sex-biased expression genes associated with sex determination and differentiation in red-tail catfish (Hemibagrus wyckioides)

BACKGROUND

Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish.

RESULTS

In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways.

CONCLUSIONS

Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

Screening and characterization of sex-specific sequences through 2b-RAD sequencing in American shad (Alosa sapidissima)

American shad (Alosa sapidissima), introduced from the United States, has become one of the most expensive farmed fish in the aquatic product market of China. The shad reveals significant sexual dimorphism in growth and behaviors. For the study, five male-specific tags were identified in two-generation breeding populations of Alosa sapidissima and were verified by PCR amplification. Averages of 10,245,091 and 8,685,704 raw and enzyme reads were obtained by high-throughput sequencing of the 2b-RAD library, respectively. 301,022 unique tags were obtained from the sequences of twenty samples with sequencing depths of 0 to 500. Finally, 274,324 special tags and 29,327 SNPs were selected with a sequencing depth of 3 to 500. Eleven preliminary screening male-specific tags and three male heterogametic SNP loci were isolated. After verification by PCR amplification, five male-specific sequences of 27 bp located on chromosome 3 were screened out. Chromosome 3 could be assumed to be the sex chromosome of Alosa sapidissima. Sex-specific markers will provide invaluable and systematic animal germplasm resources to allow for the precise identification of neo-males for the all-female breeding of Alosa sapidissima in commercial aquaculture.

Comparative gonad transcriptome analysis in cobia (Rachycentron canadum)

Background: Cobia (Rachycentron canadum) is a species of fish with high commercial potential particularly due to fast growth rates. The evidence of sexual size dimorphism favoring females indicate potential benefits in having a monosex culture. However, the involvement of genetic factors responsible for sexual development and gonadal maintenance that produces phenotypic sex in cobia is largely unknown. Methods: In the present study, we performed transcriptome sequencing of cobia to identify sex-biased significantly differentially expressed genes (DEGs) in testes and ovaries. The reliability of the gonad transcriptome data was validated by qPCR analysis of eight selected significantly differential expressed sex-related candidate genes. Results: This comparative gonad transcriptomic analysis revealed that 7,120 and 4,628 DEGs are up-regulated in testes or ovaries, respectively. Further functional annotation analyses identified 76 important candidate genes involved in sex determination cascades or sex differentiation, including 42 known testis-biased DEGs (dmrt1, amh and sox9 etc.), and 34 known ovary-biased DEGs (foxl2, sox3 and cyp19a etc.). Moreover, eleven significantly enriched pathways functionally related to sex determination and sex differentiation were identified, including Wnt signaling pathway, oocyte meiosis, the TGF-beta signaling pathway and MAPK signaling pathway. Conclusion: This work represents the first comparative gonad transcriptome study in cobia. The putative sex-associated DEGs and pathways provide an important molecular basis for further investigation of cobia's sex determination, gonadal development as well as potential control breeding of monosex female populations for a possible aquaculture setting.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Sexually Dimorphic Gene Expression in X and Y Sperms Instructs Sexual Dimorphism of Embryonic Genome Activation in Yellow Catfish (Pelteobagrus fulvidraco)

Paternal factors play an important role in embryonic morphogenesis and contribute to sexual dimorphism in development. To assess the effect of paternal DNA on sexual dimorphism of embryonic genome activation, we compared X and Y sperm and different sexes of embryos before sex determination. Through transcriptome sequencing (RNA-seq) and whole-genome bisulfite sequencing (WGBS) of X and Y sperm, we found a big proportion of upregulated genes in Y sperm, supported by the observation that genome-wide DNA methylation level is slightly lower than in X sperm. Cytokine-cytokine receptor interaction, TGF-beta, and toll-like receptor pathways play important roles in spermatogenesis. Through whole-genome re-sequencing (WGRS) of parental fish and RNA-seq of five early embryonic stages, we found the low-blastocyst time point is a key to maternal transcriptome degradation and zygotic genome activation. Generally, sexual differences emerged from the bud stage. Moreover, through integrated analysis of paternal SNPs and gene expression, we evaluated the influence of paternal inheritance on sexual dimorphism of genome activation. Besides, we screened out gata6 and ddx5 as potential instructors for early sex determination and gonad development in yellow catfish. This work is meaningful for revealing the molecular mechanisms of sex determination and sexual dimorphism of fish species.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

CircRNA Identification and CircRNA-miRNA-mRNA Network in Cynoglossus semilaevis Sexual Size Dimorphism

Sexual size dimorphism (SSD), which is the sexual differences in body size, has been widely reported in various species including fishes. For Chinese tongue sole (Cynoglossus semilaevis), a flatfish exhibiting typically female-biased SSD, little is known for its epigenetic regulation mechanism, especially the role of circRNAs. Here, we identified the differently expressed abundances of circRNAs in females, males, and pseudo-males to explore the potential functions of circRNAs in Chinese tongue sole SSD. In total, 14,745 novel circRNAs were screened, among which 1461 DE circRNAs were identified from the brain, gonad, liver, and muscle in female, male, and pseudo-male individuals. The ceRNA network was subsequently constructed, including 10 circRNAs, 26 mRNAs, and 11 miRNAs. These DE mRNAs were mainly related to the mRNA surveillance pathway, metabolic pathways, and cellular senescence. Importantly, the ceRNA network has revealed that several circRNAs such as novel_circ_004374 and novel_circ_014597 may regulate homeodomain interacting protein kinase 2 (hipk2) expression by sponging miR-130-x. It is also worth exploring whether or how novel_circ_008696 regulates SET Domain Containing 2, histone lysine methyltransferase (setd2), which in turn affects the epigenetic patterns of different sexual individuals. The present study not only enriches the knowledge on the potential roles of circRNA in the physiological process, but also provides new clues for the explanation of fish SSD. In future studies, the precise function and involvement of circRNAs in female-biased SSD will require more efforts.

Potential Involvement of ewsr1-w Gene in Ovarian Development of Chinese Tongue Sole, Cynoglossus semilaevis

Simple Summary

Sexual dimorphism is a phenomenon commonly existing in animals. Chinese tongue sole Cynoglossus semilaevis is an economical marine fish with obvious female-biased size dimorphism. So, it is important to explore the molecular mechanism beyond gonadal development for sex control in aquaculture industry. RNA-binding protein Ewing Sarcoma protein-like (ewsr1) gene is important for mouse gonadal development and reproduction, however there are limited studies on this gene in teleost. In this study, two ewsr1 genes were cloned and characterized from C. semilaevis. The ewsr1-w gene, located in W chromosomes, showed female-biased expression during C. semilaevis gonadal development. In addition, knock-down effect and transcriptional regulation of Cs-ewsr1-w further suggested its essential role in ovarian development. This study broadened our understanding on ewsr1 function in teleost, and provided genetic resources for the further development of sex control breeding techniques in C. semilaevis aquaculture.

Abstract

Ewsr1 encodes a protein that acts as a multifunctional molecule in a variety of cellular processes. The full-length of Cs-ewsr1-w and Cs-ewsr1-z were cloned in Chinese tongue sole (Cynoglossus semilaevis). The open reading frame (ORF) of Cs-ewsr1-w was 1,767 bp that encoded 589 amino acids, while Cs-ewsr1-z was 1,794 bp that encoded 598 amino acids. Real-time PCR assays showed that Cs-ewsr1-w exhibited significant female-biased expression and could be hardly detected in male. It has the most abundant expression in ovaries among eight healthy tissues. Its expression in ovary increased gradually from 90 d to 3 y with C. semilaevis ovarian development and reached the peak at 3 y. After Cs-ewsr1-w knockdown with siRNA interference, several genes related to gonadal development including foxl2, sox9b and pou5f1 were down-regulated in ovarian cell line, suggesting the possible participation of Cs-ewsr1-w in C. semilaevis ovarian development. The dual-luciferase reporter assay revealed that the -733/-154 bp Cs-ewsr1-w promoter fragment exhibited strong transcription activity human embryonic kidney (HEK) 293T cell line. The mutation of a MAF BZIP Transcription Factor K (Mafk) binding site located in this fragment suggested that transcription factor Mafk might play an important role in Cs-ewsr1-w basal transcription. Our results will provide clues on the gene expression level, transcriptional regulation and knock-down effect of ewsr1 gene during ovarian development in teleost.

First account of a transient intersex in spotted scat, Scatophagus argus: a marine gonochoristic fish

This study presents the first incidence of intersex associated with testis-ova in spotted scat (Scatophagus argus) reared in a controlled environment. The testis-ova is associated with the abnormal occurrence of primary oocytes (POs) in some male testis and is referred to as ectopic primary oocytes (Ecto-PO), whiles individuals with Ecto-PO are called "Ecto-PO gonad/individuals." We investigated gonads of 129 male spotted scat aged 4-12 and 18 months after hatch (mah) by histological studies for the presence of female sexual characteristics. A total of 20 out of 88 gonads representing 22.7% of male fish aged 6-12, or 15.5% of all male fish sampled, were found to have visible Ecto-PO. At least, the Ecto-PO had an average of 7 oocytes per gonadal section, indicating high severity. The Ecto-PO appears after sex differentiation and degenerates during sexual maturation. The Ecto-PO did not significantly influence the expression pattern of male and female sex-related genes performed using qPCR. Immunofluorescence of 42sp50 specifically stained the Ecto-PO without influence from the surrounding testicular tissues. In addition, temperature did not correlate with the proliferation of the Ecto-PO, but rather gonad developmental strategy. The results show that the naturally occurring Ecto-PO might not be detrimental to testis development and could be considered a frequent-high-level incidence of natural aberration. This study highlights the intricacy of fish sex differentiation and provides a new research chapter to ascertain the mystery behind the occurrence of Ecto-PO.

Temporal variation in DNA methylation during gonadal development in a reptile with temperature-dependent sex determination

DNA methylation plays a significant role in transducing external environmental signals to a cellular response in reptiles; however, whether the methylation patterns are conserved across species remains unclear. Here, we examined the genome-wide DNA methylation differentiation between male and female hatchling gonads of the temperature-dependent sex determination (TSD) Mauremys mutica (M. mutica) using methylation-dependent restriction-site associated DNA sequencing (MethylRAD-seq) to test differentially methylated genes underlying sexual development. Several categories, including heat shock genes (HSP90A, HSP30C), histone- (KDM8) and ubiquitin-related genes (TRIM39), kinases (WNK3) and gonad differentiation or gonadal development related genes (HSD17B8, HSD17B12), were identified as candidates for future study. Additionally, we identified several regulatory pathways potentially mediating TSD thermosensitivity such as the GnRH signaling pathway and calcium signaling pathway. These findings provide evidence that sexually dimorphic DNA methylation may be associated with sex determination or sex differentiation in TSD M. mutica.

Gene duplication, conservation, and divergence of activating transcription factor 5 gene in zebrafish

Activating transcription factor 5 (Atf5) is a member of the ATF/CREB family of transcription factors and involved in diverse cellular functions and diseases in mammals. However, the function of atf5 remains largely unknown in fish. Here, we report the expression pattern and function of duplicated atf5 genes in zebrafish. The results showed that the gene structures of zebrafish atf5a and atf5b were similar to their mammalian orthologs. Zebrafish Atf5a and Atf5b shared an amino acid sequence identity of 40.7%. Zebrafish atf5a and atf5b had maternal origin with dynamic expression during embryonic development. Zebrafish atf5a mRNA is mainly enriched in olfactory epithelium, midbrain, and hindbrain, while zebrafish atf5b mRNA is mainly detected in midbrain, hindbrain, and liver during embryogenesis. The results of acute hypoxia experiment showed that atf5a mRNA was significantly upregulated in the brain, liver, and muscle, while atf5b mRNA was just increased significantly in the brain. Functional analysis showed that knockdown of atf5a affects the development of the ciliated neurons in zebrafish embryos. The effect was enhanced when atf5a MO was co-injected with atf5b MO. The development of ciliated neurons in zebrafish embryos was not affected by injection of atf5b MO alone. atf5a knockdown also affects the development of early-born olfactory neurons. The effects caused by atf5a knockdown could be rescued by atf5b mRNA. These results suggest that the duplicated atf5 genes may have evolved divergently and play redundant biological roles in the development of olfactory sensory neurons in zebrafish.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

Fish female-biased gene cyp19a1a leads to female antiviral response attenuation between sexes by autophagic degradation of MITA

From insects to mammals, both innate and adaptive immune response are usually higher in females than in males, with the sex chromosome and hormonal differences considered the main reasons. Here, we report that zebrafish cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a), an autosomal gene with female-biased expression, causes female fish to exhibit a lower antiviral response. First, we successfully constructed an infection model by intraperitoneal injection of spring viremia of carp virus (SVCV) into zebrafish (Danio rerio) and Carassius auratus herpesvirus (CaHV) in gibel carp (Carassius gibelio). Specifically, female fish were more vulnerable to viral infection than males, accompanied by a significantly weaker interferon (IFN) expression. After screening several candidates, cyp19a1a, which was highly expressed in female fish tissues, was selected for further analysis. The IFN expression and antiviral response were significantly higher in cyp19a1a-/- than in cyp19a1a+/+. Further investigation of the molecular mechanism revealed that Cyp19a1a targets mediator of IRF3 activation (MITA) for autophagic degradation. Interestingly, in the absence of MITA, Cyp19a1a alone could not elicit an autophagic response. Furthermore, the autophagy factor ATG14 (autophagy-related 14) was found interacted with Cyp19a1a to either promote or attenuate Cyp19a1a-mediated MITA degradation by either being overexpressed or knocked down, respectively. At the cellular level, both the normal and MITA-enhanced cellular antiviral responses were diminished by Cyp19a1a. These findings demonstrated a sex difference in the antiviral response based on a regulation mechanism controlled by a female-biased gene besides sex chromosome and hormonal differences, supplying the current understanding of sex differences in fish.

DNA Methylation Difference between Female and Male Ussuri Catfish (Pseudobagrus ussuriensis) in Brain and Gonad Tissues

DNA methylation has been found to be involved in sex determination and differentiation in many aquaculture species. The Ussuri catfish (Pseudobagrus ussuriensis) is a popular aquaculture fish in China with high economic value in which male-biased sex dimorphism was observed in terms of body size and body weight. In this study, DNA methylation-sensitive RAD sequencing (Methyl-RAD) was used to explore the epigenetic difference between adult male and female samples in brain and gonad tissues. In brain tissues, 5,442,496 methylated cytosine sites were found and 9.94% of these sites were from symmetric CCGG or CCWGG sites. Among these sites, 321 differential DNA methylation sites (DMSs) in 171 genes were identified, while in gonad tissues, 4,043,053 methylated cytosines sites were found in total and 11.70% of them were from CCGG or CCWGG. Among these sites, 78 differential DNA methylation sites were found which were located in 64 genes. We also found several sex-determination genes among these differential methylated genes, such as amh, gsdf and hsd11b2 in brain tissues and slco3a1, socs2 and trim47 in gonad tissues. These results provided evidence for understanding the function of DNA methylation in the sex differentiation in Pseudobagrus ussuriensis, which further deepens the relationship between gene regulation and epigenetics.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Whole-Transcriptome Analysis Identifies Gender Dimorphic Expressions of Mrnas and Non-Coding Rnas in Chinese Soft-Shell Turtle (Pelodiscus sinensis)

In aquaculture, the Chinese soft-shelled turtle (Pelodiscus sinensis) is an economically important species with remarkable gender dimorphism in its growth patterns. However, the underlying molecular mechanisms of this phenomenon have not been elucidated well. Here, we conducted a whole-transcriptome analysis of the female and male gonads of P. sinensis. Overall, 7833 DE mRNAs, 619 DE lncRNAs, 231 DE circRNAs, and 520 DE miRNAs were identified. Some "star genes" associated with sex differentiation containing dmrt1, sox9, and foxl2 were identified. Additionally, some potential genes linked to sex differentiation, such as bmp2, ran, and sox3, were also isolated in P. sinensis. Functional analysis showed that the DE miRNAs and DE ncRNAs were enriched in the pathways related to sex differentiation, including ovarian steroidogenesis, the hippo signaling pathway, and the calcium signaling pathway. Remarkably, a lncRNA/circRNA-miRNA-mRNA interaction network was constructed, containing the key genes associated with sex differentiation, including fgf9, foxl3, and dmrta2. Collectively, we constructed a gender dimorphism profile of the female and male gonads of P. sinensis, profoundly contributing to the exploration of the major genes and potential ncRNAs involved in the sex differentiation of P. sinensis. More importantly, we highlighted the potential functions of ncRNAs for gene regulation during sex differentiation in P. sinensis as well as in other turtles.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Sexual dimorphism in Odontobutis sinensis brain-pituitary-gonad axis and liver highlighted by histological and transcriptomic approach

In this study, sexual dimorphism in Chinese dark sleeper (Odontobutis sinensis) brain-pituitary-gonad axis and liver was highlighted by histological and transcriptomic approach. The results showed that there were two significant differences between males and females. Firstly, males grew larger and faster than females. Transcriptomic analysis and qPCR validation indicated that two key growth genes, insulin-like growth factor (igf) and 25-hydroxyvitamin D-1 alpha hydroxylase (cyp27b), were more highly detected in male liver than that in female liver. Secondly, histological analysis displayed that the liver in males showed an obvious ivory fatty phenomenon with more fat vacuoles and lipid droplet aggregation compared to that in females. Transcriptomic analysis indicated that the transcript level of vitellogenin (vtg) in male liver were significantly lower than that in female liver. After 17β-estradiol (E2) treatment of primary cultured hepatocytes, the vtg mRNA expression was induced significantly, while dihydrotestosterone (DHT) treatment had little effect on it. Generally, this study will provide some ideas for further exploring the mechanism of sexual dimorphism in Odontobutis sinensis.

Sex-Specific Genomic Region Identification and Molecular Sex Marker Development of Rock Bream (Oplegnathus fasciatus)

Rock bream (Oplegnathus fasciatus) is a valuable commercial marine teleost species, which exhibits sexual dimorphism in growth performance. However, the absence of a rapid and cost-effective sex identification method based on sex-specific genetic marker has impeded study on sex determination mechanisms and breeding applications. In the present study, we firstly developed the PCR method for identifying potential sex-specific sequences in Oplegnathus fasciatus with the next-generation sequencing. Sex-specific genomic regions/loci for sex determination were discovered on Chr2 and Chr6 by genome-wide association analysis, sequencing depth, and heterozygosity comparison between females and males. Candidate sex-determining genes (CCDC63, ITR, WNT4) were furtherly detected in transcriptome data of testes and ovaries. Taken together, a male-specific 34-bp deletion on the Chr2 was identified and developed into molecular marker of sex for O. fasciatus. After validation in individuals with known phenotypic sexes, the accuracy was 100%. This study gives an insight into the mechanism of sex determination in O. fasciatus, and the gender marker is crucial both for future genomic research and for development of efficient and sustainable aquaculture practice.

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.

A High-Density Genetic Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Red Swamp Crayfish (Procambarus clarkii)

Red swamp crayfish (Procambarus clarkii) is a commercially important species in global aquaculture and most successfully invasive freshwater shrimp in China. In order to determine the genetic basis of growth- and sex-related traits, a high-density genetic linkage map was constructed using 2b-RAD sequencing technology in a full-sib family. The consensus map contains 4,878 SNP markers assigned to 94 linkage groups (LGs) and spanned 6,157.737 cM with an average marker interval of 1.26 cM and 96.93% genome coverage. The quantitative trait locus (QTL) mapping for growth and sex traits was performed for the first time. QTL mapping uncovers 28 QTLs for growth-related traits in nine LGs, explaining 7.9-14.4% of the phenotypic variation, and identifies some potential candidate growth-related genes such as mih, lamr, golgb1, nurf301, and tbcd1 within the QTL intervals. A single major locus for sex determination was revealed in LG20 that explains 59.3-63.7% of the phenotypic variations. Some candidate sex-related genes, such as vps4bl, ssrf, and acot1, were identified in the QTL intervals and found to be differentially expressed in the muscle tissues between the females and the males. Furthermore, the identified SNPs were revealed to be female heterozygotes, suggesting that red swamp crayfish might have the female heterogametic ZZ/ZW sex determination system. The present study provides a valuable resource for marker-assisted selection and genetic improvement and for further genetic and genomic research in red swamp crayfish.

Tissue-Specific Expression Pattern in Ancherythroculter nigrocauda, a Sexually Size Dimorphic Fish

Certain members of the Actinopterygii class are known to exhibit sexual dimorphism (SD) that results in major phenotypic differences between male and female fishes of a species. One of the most common differences between the two sexes is in body weight, a factor with a high economic value in aquaculture. In this study, we used RNA sequencing (RNA-seq) to study the liver and brain transcriptomes of Ancherythroculter nigrocauda, a fish exhibiting SD. Females attain about fourfold body weight of males at sexual maturity. Sample clustering showed that both sexes were grouped well with their sex phenotypes. In addition, 2,395 and 457 differentially expressed genes (DEGs) were identified in the liver and brain tissues, respectively. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses predicted the association of PPAR signaling, cytochrome P450, and steroid hormone biosynthesis to the differences in sexual size. In addition, weighted gene co-expression network analyses (WGCNA) were conducted, and the green module was identified to be significantly correlated with sexual size dimorphism (SSD). Altogether, these results improve our understanding of the molecular mechanism underlying SSD in A. nigrocauda.

Integration of Transcriptome and Methylome Highlights the Roles of Cell Cycle and Hippo Signaling Pathway in Flatfish Sexual Size Dimorphism

Sexual size dimorphism (SSD) is the difference in segments or body size between sexes prevalent in various species. Understanding the genetic architecture of SSD has remained a significant challenge owing to the complexity of growth mechanisms and the sexual influences among species. The Chinese tongue sole (Cynoglossus semilaevis), which exhibits a female-biased SSD and sex reversal from female to pseudomale, is an ideal model for exploring SSD mechanism at the molecular level. The present study aimed to integrate transcriptome and methylome analysis to unravel the genetic and epigenetic changes in female, male, and pseudomale C. semilaevis. The somatotropic and reproductive tissues (brain, liver, gonad, and muscle) transcriptomes were characterized by RNA-seq technology. Transcriptomic analysis unravelled numerous differentially expressed genes (DEGs) involved in cell growth and death-related pathways. The gonad and muscle methylomes were further employed for screening differentially methylated genes (DMGs). Relatively higher DNA methylation levels were observed in the male and pseudomale individuals. In detail, hypermethylation of the chromosome W was pronounced in the pseudomale group than in the female group. Furthermore, weighted gene co-expression network analysis showed that turquoise and brown modules positively and negatively correlated with the female-biased SSD, respectively. A combined analysis of the module genes and DMGs revealed the female-biased mRNA transcripts and hypomethylated levels in the upstream and downstream regions across the cell cycle-related genes. Moreover, the male and pseudomale-biased gene expression in the hippo signaling pathway were positively correlated with their hypermethylation levels in the gene body. These findings implied that the activation of the cell cycle and the inhibition of the hippo signaling pathway were implicated in C. semilaevis female-biased SSD. In addition, the dynamic expression pattern of the epigenetic regulatory factors, including dnmt1, dnmt3a, dnmt3b, and uhrf1, among the different sexes correspond with their distinct DNA methylation levels. Herein, we provide valuable clues for understanding female-biased SSD in C. semilaevis.

ASER: Animal Sex Reversal Database

Sex reversal, representing extraordinary sexual plasticity during the life cycle, not only triggers reproduction in animals but also affects reproductive and endocrine system-related diseases and cancers in humans. Sex reversal has been broadly reported in animals; however, an integrated resource hub of sex reversal information is still lacking. Here, we constructed a comprehensive database named ASER (Animal Sex Reversal) by integrating sex reversal-related data of 18 species from teleostei to mammalia. We systematically collected 40,018 published papers and mined the sex reversal-associated genes (SRGs), including their regulatory networks, from 1611 core papers. We annotated homologous genes and computed conservation scores for whole genomes across the 18 species. Furthermore, we collected available RNA-seq datasets and investigated the expression dynamics of SRGs during sex reversal or sex determination processes. In addition, we manually annotated 550 in situ hybridization (ISH), fluorescence in situ hybridization (FISH), and immunohistochemistry (IHC) images of SRGs from the literature and described their spatial expression in the gonads. Collectively, ASER provides a unique and integrated resource for researchers to query and reuse organized data to explore the mechanisms and applications of SRGs in animal breeding and human health. The ASER database is publicly available at http://aser.ihb.ac.cn/.

Screening and identification of female-specific DNA sequences in octaploid sturgeon using comparative genomics with high-throughput sequencing

In this study, six candidate female-specific DNA sequences of octaploid Amur sturgeon (Acipenser schrenckii) were identified using comparative genomic approaches with high-throughput sequencing data. Their specificity was confirmed by traditional PCR. Two of these sex-specific sequences were also validated as female-specific in other eight sturgeon species and two hybrid sturgeons. The identified female-specific DNA fragments suggest that the family Acipenseridae has a ZZ/ZW sex-determining system. However, one of the two DNA sequences has been deleted in some sturgeons such as Sterlet sturgeon (Acipenser ruthenus), Beluga (Huso huso) and Kaluga (H. dauricus). The difference of sex-specific sequences among sturgeons indicates that there are different sex-specific regions among species of sturgeon. This study not only provided the sex-specific DNA sequences for management, conservation and studies of sex-determination mechanisms in sturgeons, but also confirmed the capability of the workflow to identify sex-specific DNA sequences in the polyploid species with complex genomes.

Transcriptomic analysis revealed gene expression profiles during the sex differentiation of Chinese tongue sole (Cynoglossus semilaevis)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process.

Characterization of Chinese tongue sole (Cynoglossus semilaevis) 24-dehydrocholesterol reductase: Expression profile, epigenetic modification, and its knock-down effect

The sexual size dimorphism of the Chinese tongue sole (Cynoglossus semilaevis) has greatly obstructed its sustainable development; however, the underlying mechanism remains unclear. Based on C. semilaevis transcriptomic information, 24-dehydrocholesterol reductase (dhcr24) was identified in steroid biosynthesis, showing female-liver-biased expression. Dhcr24 has been reported to participate in various processes, such as cholesterol synthesis, oxidative stress response, neuroprotection, and cell survival. The present study assessed its role in the sexual size dimorphism in fish. First, detailed expression pattern analysis showed that dhcr24 mRNAs were extensively expressed in tissues and the highest levels were found in the liver and gonads of females. Analysis of the dhcr24 promoter region demonstrated different DNA methylation statuses in female, male, and pseudomale gonads with higher epigenetic modification in males. The confirmation of transcription activity of the dhcr24 promoter and putative transcription factors (e.g., ER, AR, SREBP, and POU1F1a) provides the foundation for studying its regulatory mechanism. Finally, dhcr24-siRNA mediated knock-down assay using C. semilaevis liver cells showed that steroid biosynthesis related genes (e.g., ebp, dhcr7, and sc5d), core component of PI3K/Akt pathway (e.g., pi3k), and igf1r exhibited different expression patterns. Further investigation on the interplay between steroid hormones, dhcr24, PI3K/Akt, and IGF-1 systems will be valuable to better understand the mechanism underlying the sexual size dimorphism in C. semilaevis.

The Snakeskin Gourami (Trichopodus pectoralis) Tends to Exhibit XX/XY Sex Determination

The snakeskin gourami (Trichopodus pectoralis) has a high meat yield and is one of the top five aquaculture freshwater fishes in Thailand. The species is not externally sexually dimorphic, and its sex determination system is unknown. Understanding the sex determination system of this species will contribute to its full-scale commercialization. In this study, a cytogenetic analysis did not reveal any between-sex differences in chromosomal patterns. However, we used genotyping-by-sequencing to identify 4 male-linked loci and 1 female-linked locus, indicating that the snakeskin gourami tends to exhibit an XX/XY sex determination system. However, we did not find any male-specific loci after filtering the loci for a ratio of 100:0 ratio of males:females. This suggests that the putative Y chromosome is young and that the sex determination region is cryptic. This approach provides solid information that can help identify the sex determination mechanism and potential sex determination regions in the snakeskin gourami, allowing further investigation of genetic improvements in the species.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Emergence of breeding tubercles and puberty onset in male zebrafish: evidence for a dependence on body growth

The presence of breeding tubercles (BTs) on the pectoral fins has been investigated as a typical male secondary sexual characteristic (SSC) that distinguish males from females in adult zebrafish. Nonetheless, the earliest occurrence of these tubercles and its association with puberty onset and body growth remain unclear. In this study, using morphological, histological and statistical analyses, the authors examined the first appearance of BTs and puberty onset in male zebrafish, with particular emphasis on the potential impact of body growth on them. The results of this study revealed that BTs distributed along the first five branched pectoral fin rays were the earliest manifestation of male SSCs, which is significantly strongly correlated with body weight (R²  = 0.9609, P < 0.001), and could be used as a "gold standard" for the earliest sex distinction (<0.1 g in weight). Using the first appearance of BTs (<0.20 mm² ) as a metric, the authors established that male puberty commenced at a body weight of c. 0.056 ± 0.015 g or a standard length of 10.99 ± 1.051 mm (mean ± S.D.). In this study, the authors thus established a simple method that can be used to sex live zebrafish at the pubertal stage and provides the first evidence for the relationship of BTs and male puberty initiation with body growth. These findings will accordingly lay a foundation for exploring mechanisms of the SSCs and male puberty onset in zebrafish and other teleost fish.

Characterization, expression and CpG methylation analysis of Dmrt1 and its response to steroid hormone in blotched snakehead (Channa maculata)

Dmrt1 is an important transcriptional regulator that plays critical role in male gonadogenesis, testicular differentiation and development. In this study, Dmrt1 was cloned from blotched snakehead (Channa maculata), which is designated as CmDmrt1. CmDmrt1 encoded a putative protein with 293 amino acids and presented an extremely conserved DM domain. It was nearly expressed in the gonads, and the expression was more than 15 times higher in the testis than in the ovary. 1851 bp promoter sequence of CmDmrt1 was characterized and the methylation levels of the CpG sites were analyzed to detect sex-related differences. A significant negative correlation between CmDmrt1 expression and CpG methylation level of its promoter was found in the testis and ovary. During gonadal development, CmDmrt1 transcription displayed strong male-biased expression patterns, increased with the maturation of testis and reached the peak at 195 days after hatching (dah), which indicates a significant role of Dmrt1 in spermatogenesis. Steroid treatment could influence CmDmrt1 expression, and long-term 17β-estradiol (E2) treatment could induce the male-to-female secondary sex reversal (SSR), which resulted in the differentiated testis transformed to ovary or ovotestis. Meanwhile, CmDmrt1 expression was down-regulated to fairly low level in the ovary of the SSR XY fish, which was similar to that in normal XX females ovary. Our research illustrates that Dmrt1 is linked to testis differentiation and spermatogenesis in blotched snakehead, providing information for functional studies on sex differentiation and gonadal development of C. maculata, and scientific basis for the production practice of all-male snakehead breeding.

Identification of Sex-Specific Markers Through 2b-RAD Sequencing in the Sea Urchin (Mesocentrotus nudus)

Sex-specific markers play an important role in revealing sex-determination mechanism. Sea urchin (Mesocentrotus nudus) is an economically important mariculture species in several Asian countries and its gonads are the sole edible parts for people. However, growth rate and immunocompetence differ by sex in this species, sex-specific markers have not been identified, and the sex-determination mechanism of sea urchin remains undetermined. In this study, type IIB endonuclease restriction-site associated DNA sequencing (2b-RAD-seq) and a genome survey of M. nudus were performed, and three female-specific markers and three female heterogametic single nucleotide polymorphism (SNP) loci were identified. We validated these sex-specific markers via PCR amplification in a large number of individuals, including wild and artificially bred populations. Several open reading frames (ORFs) were predicted, although there are no potential genes known for sex determination and sex differentiation within the scaffold in which the sex-specific markers are located. Importantly, the female-specific sequences and female heterozygous SNP loci indicate that a female heterogametic and male homogametic ZW/ZZ sex-determination system should exist in M. nudus. The results provide a solid basis for revealing the sex-determination mechanism of this species, and open up new possibilities for developing sex-control breeding in sea urchin.

Testis-specific expression pattern of dmrt1 and its putative regulatory region in the sea urchin (Mesocentrotus nudus)

Sea urchin (Mseocentrotus nudus) is an economically important mariculture species in several Asian countries. The growth rate and immunocompetence differ by sex in this species. However, the mechanisms of sex determination in M. nudus have remain unclear. In the present study, we focus on the dmrt1 gene of M. nudus (Mndmrt1) to investigate its dynamic expression pattern during different developmental stages. Real-time quantitative PCR (RT-qPCR) revealed that Mndmrt1 exhibits testis-specific expression and undetectable during the whole embryogenesis. With the development of ontogenetic, Mndmrt1 transcripts are first detected at 9 months post-fertilization (mpf). In addition, both the transcripts and protein of Mndmrt1 gene were specifically expressed in spermatogonia and spermatocytes, indicating that it might be a male germ cells marker in sea urchin. Significantly, the 1441 bp promoter sequence of Mndmrt1 gene was obtained by DNA walking, and one positive regulatory region at -1197/ -968 in the promoter, as well as one negative regulatory region at -1441/ -1198 have been identified by promoter activity analysis. Moreover, two regulatory regions contain multiple putative binding sites for transcription factors, including Sp1, Egr1, Sox5, CEBP, GATA and SRY. These findings suggest that Mndmrt1 may be related to testis differentiation and spermatogenesis in sea urchin and will provide an insight into understanding the regulatory mechanism of the dmrt1 gene.

A sex-linked SNP mutation in amhr2 is responsible for male differentiation in obscure puffer (Takifugu obscurus)

BACKGROUND

Anti-Mullerian hormone receptor type II (Amhr2) is a key receptor of Amh signaling in regulating gonad development. The amhr2 gene has been identified in numerous species, including a few teleost fishes. However, the roles of Amhr2 in Amh signaling in fish are poorly studied.

METHODS AND RESULTS

In this study, an amhr2 homolog from obscure puffer (Takifugu obscurus) was identified, and its molecular characteristics were systematically analyzed. Expression analysis revealed that amhr2 was highly expressed in the gonads of adult pufferfish and significantly upregulated during sex differentiation. Significantly, a sex-linked SNP site was verified in obscure puffer amhr2. Females exhibited a homozygous genotype (C/C), while males possessed a heterozygous genotype (C/G), resulting in an amino acid variation (His/Asp384) in the kinase domain of Amhr2. Then, the functions of the different Amhr2 genotypes were further investigated. The male genotype protein (Amhr2^D384) showed an enhanced ability to interact with the type I receptor (Bmpr1a) compared to the female genotype (Amhr2^H384). The phosphorylation levels of Smads and activity of the target gene (id3) induced by the male genotype were also much higher than those induced by the female genotype. These results confirmed that the male genotype had an enhanced effect on the Amh signaling pathway compared with the female genotype.

CONCLUSIONS

This study provides direct experimental evidence for the roles of different Amhr2 genotypes in pufferfish and suggests that amhr2 is responsible for male sex differentiation in obscure puffer.

The application of genome editing technology in fish

The advent and development of genome editing technology has opened up the possibility of directly targeting and modifying genomic sequences in the field of life sciences with rapid developments occurring in the last decade. As a powerful tool to decipher genome data at the molecular biology level, genome editing technology has made important contributions to elucidating many biological problems. Currently, the three most widely used genome editing technologies include: zinc finger nucleases (ZFN), transcription activator like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats (CRISPR). Researchers are still striving to create simpler, more efficient, and accurate techniques, such as engineered base editors and new CRISPR/Cas systems, to improve editing efficiency and reduce off-target rate, as well as a near-PAMless SpCas9 variants to expand the scope of genome editing. As one of the important animal protein sources, fish has significant economic value in aquaculture. In addition, fish is indispensable for research as it serves as the evolutionary link between invertebrates and higher vertebrates. Consequently, genome editing technologies were applied extensively in various fish species for basic functional studies as well as applied research in aquaculture. In this review, we focus on the application of genome editing technologies in fish species detailing growth, gender, and pigmentation traits. In addition, we have focused on the construction of a zebrafish (Danio rerio) disease model and high-throughput screening of functional genes. Finally, we provide some of the future perspectives of this technology.

Evidence of Sex Differentiation Based on Morphological Traits During the Early Development Stage of Mud Crab Scylla paramamosain

In order to uncover the sexual difference in morphology and how early they appear during the development stage of mud crab Scylla paramamosain, we measured, observed, and biostatistically analyzed morphological traits related to sex. For unveiling the morphological differences between sexes, morphological traits involving abdomen width (AW), carapace length (CL), and carapace width (CW) were first measured during the crablet development stage of S. paramamosain in the present study. The correlation analyses and path analyses exhibited that sexual dimorphism in the third abdomen width (AW3) and fourth abdomen width (AW4) could be used for sex identification from stage C VI (stage VI of crablet). Based on the stepwise discriminant analysis and standardized traits, a sex discriminant equation was constructed, which is capable for sex identification in crablets from stage C VI. Observations for secondary sexual traits and abdomen morphology (shape and pleopods) using a dissecting microscope or scanning electron microscope indicated that sexes are easily identified at stage C VIII according to the abdomen shape; meanwhile, at stage C II based on pleopod difference, and at stage C I by the presence or absence of gonopores. The findings in this study contribute greatly to the accuracy of sex identification of S. paramamosain during the early development stage, which promotes the understanding of the morphological differentiation mechanism of sex.