Sex-determining Gene Research Articles

Gene Variant of Barrier to Autointegration Factor 2 (Banf2w) Is Concordant with Female Determination in Cichlids.

Oreochromis fishes exhibit variability of sex-determination (SD) genes whose characterization contributes to understanding of the sex differentiation network, and to effective tilapia farming, which requires all-male culture. However, O. niloticus (On) amh is the only master-key regulator (MKR) of SD that has been mapped (XY/XX SD-system on LG23). In O. aureus (Oa), LG3 controls a WZ/ZZ SD-system that has recently been delimited to 9.2 Mbp, with an embedded interval rich with female-specific variation, harboring two paics genes and banf2. Developing genetic markers within this interval and using a hybrid Oa stock that demonstrates no recombination repression in LG3, we mapped the critical SD region to 235 Kbp on the orthologous On physical map (p < 1.5 × 10−26). DNA-seq assembly and peak-proportion analysis of variation based on Sanger chromatograms allowed the characterization of copy-number variation (CNV) of banf2. Oa males had three exons capable of encoding 90-amino-acid polypeptides, yet in Oa females, we found an extra copy with an 89-amino-acid polypeptide and three non-conservative amino acid substitutions, designated as banf2w. CNV analysis suggested the existence of two to five copies of banf2 in diploidic Cichlidae. Disrupting the Hardy–Weinberg equilibrium (p < 4.2 × 10−3), banf2w was concordant with female determination in Oa and in three cichlids with LG3 WZ/ZZ SD-systems (O. tanganicae, O. hornorum and Pelmatolapia mariae). Furthermore, exclusive RNA-seq expression in Oa females strengthened the candidacy of banf2w as the long-sought LG3 SD MKR. As banf genes mediate nuclear assembly, chromatin organization, gene expression and gonad development, banf2w may play a fundamental role inducing female nucleus formation that is essential for WZ/ZZ SD.

Doublesex mediates species-, sex-, environment- and trait-specific exaggeration of size and shape.

Context-dependent trait exaggeration is a major contributor to phenotypic diversity. However, the genetic modifiers instructing development across multiple contexts remain largely unknown. We use the arthropod tibia, a hotspot for segmental differentiation, as a paradigm to assess the developmental mechanisms underlying the context-dependent structural exaggeration of size and shape through nutritional plasticity, sexual dimorphism and segmental differentiation. Using an RNAseq approach in the sexually dimorphic and male-polyphenic dung beetle Digitonthophagus gazella, we find that only a small portion (3.7%) of all transcripts covary positively in expression level with trait size across contexts. However, RNAi-mediated knockdown of the conserved sex-determination gene doublesex suggests that it functions as a context-dependent master mediator of trait exaggeration in D. gazella as well as the closely related dung beetle Onthophagus taurus. Taken together, our findings suggest (i) that the gene networks associated with trait exaggeration are highly dependent on the precise developmental context, (ii) that doublesex differentially shapes morphological exaggeration depending on developmental contexts and (iii) that this context-specificity of dsx-mediated trait exaggeration may diversify rapidly. This mechanism may contribute to the resolution of conflict arising from environment-dependent antagonistic selection among sexes and divergent developmental contexts in a wide range of animals.

A role for sex-determination genes in life history evolution? Doublesex mediates sexual size dimorphism in the gazelle dung beetle.

An organism's fitness depends strongly on its age and size at maturation. Although the evolutionary forces acting on these critical life history traits have been heavily scrutinized, the developmental mechanisms underpinning intraspecific variation in adult size and development time remain much less well-understood. Using RNA interference, I here show that the highly conserved sex-determination gene doublesex (dsx) mediates sexual size dimorphism (SSD) in the gazelle dung beetle Digitonthophagus gazella. Because doublesex undergoes sex-specific splicing and sex-limited isoforms regulate different target genes, this suggests that dsx contributes to the resolution of intralocus sexual conflict in body size. However, these results contrast with previous studies demonstrating that dsx does not affect body size or SSD in Drosophila. This indicates that intraspecific body size variation is underlain by contrasting developmental mechanisms in different insect lineages. Furthermore, although male D.gazella have a longer development time than females, sexual bimaturism was not affected by dsx expression knockdown. In addition, and in contrast to secondary sexual morphology, dsx did not significantly affect nutritional plasticity in life history. Taken together, these findings indicate that dsx signalling contributes to intraspecific life history variation but that dsx's function in mediating sexual dimorphism in life history differs among traits and species. More generally, these findings suggest that genes ancestrally tasked with sex determination have been co-opted into the developmental regulation of life history traits and may represent an underappreciated mechanism of life history evolution.

Comparative Transcriptome Analysis Reveals Sex-Based Differences during the Development of the Adult Parasitic Wasp Cotesia vestalis (Hymenoptera: Braconidae).

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.

In Silico Prediction of Transcription Factor Collaborations Underlying Phenotypic Sexual Dimorphism in Zebrafish (Danio rerio).

The transcriptional regulation of gene expression in higher organisms is essential for different cellular and biological processes. These processes are controlled by transcription factors and their combinatorial interplay, which are crucial for complex genetic programs and transcriptional machinery. The regulation of sex-biased gene expression plays a major role in phenotypic sexual dimorphism in many species, causing dimorphic gene expression patterns between two different sexes. The role of transcription factor (TF) in gene regulatory mechanisms so far has not been studied for sex determination and sex-associated colour patterning in zebrafish with respect to phenotypic sexual dimorphism. To address this open biological issue, we applied bioinformatics approaches for identifying the predicted TF pairs based on their binding sites for sex and colour genes in zebrafish. In this study, we identified 25 (e.g., STAT6-GATA4; JUN-GATA4; SOX9-JUN) and 14 (e.g., IRF-STAT6; SOX9-JUN; STAT6-GATA4) potentially cooperating TFs based on their binding patterns in promoter regions for sex determination and colour pattern genes in zebrafish, respectively. The comparison between identified TFs for sex and colour genes revealed several predicted TF pairs (e.g., STAT6-GATA4; JUN-SOX9) are common for both phenotypes, which may play a pivotal role in phenotypic sexual dimorphism in zebrafish.

Efficient TALEN-mediated gene knockin at the bovine Y chromosome and generation of a sex-reversal bovine

Functional elucidation of bovine Y-chromosome genes requires available genome editing technologies. Meanwhile, it has yet to be proven whether the bovine Sry gene is the main or single factor involved in the development of the male phenotype in bovine. Here, we efficiently knocked out four Y-linked genes (Sry, ZFY, DDX3Y, and EIF2S3Y) in bovine fetal fibroblasts (BFFs) with transcription activator-like effector nucleases (TALENs) individually. Furthermore, we used TALEN-mediated gene knockin at the Sry gene and generated a sex-reversal bovine by somatic cell nuclear transfer (SCNT). The resulting bovine had only one ovary and was sterile. We demonstrate, for the first time, that the Sry gene is an important sex-determining gene in bovine. Our method lays a solid foundation for detecting the biology of the bovine Y chromosome, as it may provide an alternative biological model system for the study of mammalian sex determination, and new methods for the practical application in agricultural, especially for sex predetermination.

Possible Epigenetic Origin of a Recurrent Gynandromorph Pattern in Megachile Wild Bees.

Simple SummaryGynandromorphs, i.e., individuals with a mix of male and female body parts, are known for many species of insects and other animals with separate sexes. This anomaly is generally regarded as the result of localized genetic mutations in sex-determining genes. We analyzed the specific mix of male and female characters in naturally occurring gynandromorphs of 21 species of the wild bee genus Megachile and found a recurrent pattern. Based on the regularity of this pattern, and the current knowledge on sex determination and sex differentiation in the relatively closely-related honey bee, we argue that the origin of these composite phenotypes is possibly epigenetic, rather than genetic, i.e., produced by some defects in the maintenance of the regulatory signals that control sex differentiation at the level of single cell lineages, rather than triggered by genetic mutations.Gynandromorphs, i.e., individuals with a mix of male and female traits, are common in the wild bees of the genus Megachile (Hymenoptera, Apoidea). We described new transverse gynandromorphs in Megachile pilidens Alfkeen, 1924 and analyze the spatial distribution of body parts with male vs. female phenotype hitherto recorded in the transverse gynandromorphs of the genus Megachile. We identified 10 different arrangements, nine of which are minor variants of a very general pattern, with a combination of male and female traits largely shared by the gynandromorphs recorded in 20 out of 21 Megachile species in our dataset. Based on the recurrence of the same gynandromorph pattern, the current knowledge on sex determination and sex differentiation in the honey bee, and the results of recent gene-knockdown experiments in these insects, we suggest that these composite phenotypes are possibly epigenetic, rather than genetic, mosaics, with individual body parts of either male or female phenotype according to the locally expressed product of the alternative splicing of sex-determining gene transcripts.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy.

Simple Summarydmy is considered to be the male-determining gene in Japanese medaka (Oryzias latipes) and Hainan medaka (Oryzias. curvinotus), both of which have the XX/XY sex-determination system. Here, we found a group of medaka in the Sanya River (named SY-medaka) and confirmed that SY-medaka belongs to O. curvinotus by morphological characteristics and mitochondrial phylogenetic analysis. Through genetic sex identification, genome re-sequencing and gonadal transcriptome analysis, it was preliminary confirmed that SY-medaka did not contain dmy. Our results provide a basis for further studies of the mechanism underlying sex determination in Oryzias and functional genomics and reproduction biology in O. curvinotus.Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

Transcriptome analysis provides insights into long noncoding RNAs in medaka gonads

Long non-coding RNAs (lncRNAs) are gradually regarded as regulators in sex determination and gonad development of various animals. Medaka (Oryzias latipes) is an excellent reproductive research model with sex-determining genes. However, the regulation of gonadal lncRNAs on medaka reproductive development remains unknown. Here, 5317 lncRNAs were obtained from medaka ovary and testis by Illumina HiSeq4000, among which 177 lncRNAs were up-regulated and 120 lncRNAs were down-regulated in the testis compared to the ovary. In addition, 6904 cis-regulated target genes were predicted from 3099 lncRNAs. GO and KEGG enrichment analysis showed that these target genes were mainly involved in phosphorylation, metabolic, metabolism of xenobiotics by cytochrome P450, insulin secretion, and GnRH signaling pathways. Furthermore, six highly expressed lncRNAs were randomly selected to verify the sequencing data by quantitative real time PCR (qRT-PCR). Next, in situ hybridization revealed that one of the sex-biased lncRNA MSTRG.14827.1 was highly expressed in immature germ cells, indicating MSTRG.14827.1 may play a key role in gametogenesis. Taken together, this study provides emerging lncRNA libraries and opens new avenues for future investigation of lncRNAs in medaka.

Toxic effects of octocrylene on zebrafish larvae and liver cell line (ZFL).

Octocrylene (OC) is a broad-spectrum ultraviolet-absorbing chemical used in sunscreen and other personal care products. Its health effects are a concern because it has been detected in water, fish, humans, and food chains. In vivo and in vitro investigations were performed in zebrafish (Danio rerio) larvae and a zebrafish liver cell line (ZFL), respectively, to understand the potential risks and molecular mechanisms of OC toxicity. The 96-h median lethal concentration (LC50) of OC was determined to be 251.8μM in larvae and 5.5μM in ZFL cells. Quantitative real-time PCR (qRT-PCR) showed that OC induced the expression of genes for CYPs (CYP1A, CYP3A65), estrogen receptors (ERα, ERβ1, GPER), vitellogenin (VTG1), and sex determination (BRCA2, CYP19A, DMRT1, SOX9A), both in vitro and in vivo. A whole-transcriptome sequencing method was used to evaluate the gene expression profile of larvae exposed to OC. OC was found to mediate the biosynthesis of estrogens (such as estriol) and affect the antioxidant pathway (glutathione transferases and peroxisome). These findings clarify the toxic effects and molecular mechanisms of OC and support banning its use in cosmetics.

Caste-biased Expression of fem and Amdsx Genes in Apis mellifera ligustica (Hymenoptera: Apidae)

Sex determination and caste differentiation are two crucial processes for morphology building in honeybees. It is unclear whether there is an interaction between these two processes. Here, we investigated the expression of fem and Amdsx genes between female castes of honeybees. We found that the expression of fem and Amdsx is significantly higher in queens than in workers, and this expression was positively regulated by juvenile hormone (JH). Our results suggest that sex-determining genes fem and Amdsx are also involved in honeybee caste differentiation.

Chromosome level reference of Atlantic halibut Hippoglossus hippoglossus provides insight into the evolution of sexual determination systems.

Changes in the genetic mechanisms that control sexual determination have occurred independently across the tree of life, and with exceptional frequency in teleost fishes. To investigate the genomic changes underlying the evolution of sexual determination, we sequenced a chromosome-level genome, multitissue transcriptomes, and reduced representation population data for the Atlantic halibut (Hippoglossus hippoglossus), which has an XY/XX sex determination mechanism and has recently diverged (0.9-3.8Ma) from the Pacific halibut (Hippoglossus stenolepis), which has a ZZ/ZW system. We used frequency and coverage-based population approaches to identify a putative sex-determining factor, GSDF. We characterized regions with elevated heterozygosity and linkage disequilibrium indicating suppression of recombination across a nascent sex chromosome. We detected testis-specific expression of GSDF, the sequence of which is highly conserved across flatfishes. Based on evidence from genome-wide association, coverage, linkage disequilibrium, testis and brain transcriptomes, and sequence conservation with other flatfishes, we propose a mechanism for the recent evolution of an XY sex-determination mechanism in Atlantic halibut. Changes to the ancestral sex-determining gene DMRT1 in regulating the downstream gene GSDF probably coincided with GSDF, or a proximal regulatory element of it, becoming the primary sex-determining factor. Our results suggest changes to a small number of elements can have drastic repercussions for the genomic substrate available to sex-specific evolutionary forces, providing insight into how certain elements repeatedly evolve to control sex across taxa. Our chromosome-level assembly, multitissue transcriptomes, and population genomic data provide a valuable resource and understanding of the evolution of sexual systems in fishes.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration.

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Evolution of a Multiple Sex-Chromosome System by Three-Sequential Translocations among Potential Sex-Chromosomes in the Taiwanese Frog Odorrana swinhoana.

Translocation between sex-chromosomes and autosomes generates multiple sex-chromosome systems. It happens unexpectedly, and therefore, the evolutionary meaning is not clear. The current study shows a multiple sex chromosome system comprising three different chromosome pairs in a Taiwanese brown frog (Odorrana swinhoana). The male-specific three translocations created a system of six sex-chromosomes, ♂X1Y1X2Y2X3Y3-♀X1X1X2X2X3X3. It is unique in that the translocations occurred among three out of the six members of potential sex-determining chromosomes, which are known to be involved in sex-chromosome turnover in frogs, and the two out of three include orthologs of the sex-determining genes in mammals, birds and fishes. This rare case suggests sex-specific, nonrandom translocations and thus provides a new viewpoint for the evolutionary meaning of the multiple sex chromosome system.

Reconstruction of the Origin of a Neo-Y Sex Chromosome and Its Evolution in the Spotted Knifejaw, Oplegnathus punctatus.

Sex chromosomes are a peculiar constituent of the genome because the evolutionary forces that fix the primary sex-determining gene cause genic degeneration and accumulation of junk DNA in the heterogametic partner. One of the most spectacular phenomena in sex chromosome evolution is the occurrence of neo-Y chromosomes, which lead to X1X2Y sex-determining systems. Such neo-sex chromosomes are critical for understanding the processes of sex chromosome evolution because they rejuvenate their total gene content. We assembled the male and female genomes at the chromosome level of the spotted knifejaw (Oplegnathus punctatus), which has a cytogenetically recognized neo-Y chromosome. The full assembly and annotation of all three sex chromosomes allowed us to reconstruct their evolutionary history. Contrary to other neo-Y chromosomes, the fusion to X2 is quite ancient, estimated at 48 Ma. Despite its old age and being even older in the X1 homologous region which carries a huge inversion that occurred as early as 55–48 Ma, genetic degeneration of the neo-Y appears to be only moderate. Transcriptomic analysis showed that sex chromosomes harbor 87 genes, which may serve important functions in the testis. The accumulation of such male-beneficial genes, a large inversion on the X1 homologous region and fusion to X2 appear to be the main drivers of neo-Y evolution in the spotted knifejaw. The availability of high-quality assemblies of the neo-Y and both X chromosomes make this fish an ideal model for a better understanding of the variability of sex determination mechanisms and of sex chromosome evolution.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing.

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

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.

Sexual dimorphism: Neural circuit switches in the Drosophila brain

The sex-determining genes Double Sex and Fruitless are expressed in sexually differentiated neurons of the Drosophila brain. A tiny cluster of neurons, aDN cells, serves as a key circuit switch with sexually dimorphic properties: those of female flies respond to visual signals in males, while those of male flies respond to smell and humidity in females, supporting effective courtship and communal egg laying behaviors, respectively.

A long reads-based de-novo assembly of the genome of the Arlee homozygous line reveals chromosomal rearrangements in rainbow trout.

Currently, there is still a need to improve the contiguity of the rainbow trout reference genome and to use multiple genetic backgrounds that will represent the genetic diversity of this species. The Arlee doubled haploid line was originated from a domesticated hatchery strain that was originally collected from the northern California coast. The Canu pipeline was used to generate the Arlee line genome de-novo assembly from high coverage PacBio long-reads sequence data. The assembly was further improved with Bionano optical maps and Hi-C proximity ligation sequence data to generate 32 major scaffolds corresponding to the karyotype of the Arlee line (2 N = 64). It is composed of 938 scaffolds with N50 of 39.16 Mb and a total length of 2.33 Gb, of which ∼95% was in 32 chromosome sequences with only 438 gaps between contigs and scaffolds. In rainbow trout the haploid chromosome number can vary from 29 to 32. In the Arlee karyotype the haploid chromosome number is 32 because chromosomes Omy04, 14 and 25 are divided into six acrocentric chromosomes. Additional structural variations that were identified in the Arlee genome included the major inversions on chromosomes Omy05 and Omy20 and additional 15 smaller inversions that will require further validation. This is also the first rainbow trout genome assembly that includes a scaffold with the sex-determination gene (sdY) in the chromosome Y sequence. The utility of this genome assembly is shown through the improved annotation of the duplicated genome loci that harbor the IGH genes on chromosomes Omy12 and Omy13.

Differential expression analysis and identification of sex-related genes by gonad transcriptome sequencing in estradiol-treated and non-treated Ussuri catfish Pseudobagrus ussuriensis.

The Ussuri catfish (Pseudobagrus ussuriensis) has an XX/XY sex determination system but its sex determination gene(s) remain unknown. To better understand the molecular sex determination mechanism, transcriptome analysis was conducted to obtain sex-related gene expression profiles. Transcriptome analyses were made of male and female developing/differentiating gonads by high-throughput RNA sequencing, including gonads from fish given an estradiol-induced sex reversal treatment. A total of 81,569 unigenes were assembled and 39,904 were significantly matched to known unique proteins by comparison with public databases. Twenty specifically expressed and 142 differentially expressed sex-related genes were extracted from annotated data by comparing the treatment groups. These genes are involved in spermatogenesis (e.g., Dnali1, nectin3, klhl10, mybl1, Katnal1, Eno4, Mns1, Spag6, Tsga10, Septin7), oogenesis (e.g., Lagr5, Fmn2, Npm2, zar1, Fbxo5, Fbxo43, Prdx4, Nrip1, Lfng, Atrip), gonadal development/differentiation (e.g., Cxcr4b, Hmgb2, Cftr, Ch25h, brip1, Prdm9, Tdrd1, Star, dmrt1, Tut4, Hsd17b12a, gdf9, dnd, arf1, Spata22), and estradiol response (e.g., Mmp14, Lhcgr, vtg1, vtg2, esr2b, Piwil1, Aifm1, Hsf1, gdf9). Dmrt1 and gdf9 may play an essential role in sex determination in P. ussuriensis. The expression patterns of six random genes were validated by quantitative real-time PCR, which confirmed the reliability and accuracy of the RNA-seq results. These data provide a valuable resource for future studies of gene expression and for understanding the molecular mechanism of sex determination/differentiation and gonadal development/differentiation (including hormone-induced sexual reversal) in Ussuri catfish. This has the potential to assist in producing monosex Ussuri catfish to increase aquacultural productivity.