Three Copies of zbed1 Specific in Chromosome W Are Essential for Female-Biased Sexual Size Dimorphism in Cynoglossus semilaevis.

The sex chromosome, especially specific in one sex, generally determines sexual size dimorphism (SSD), a phenomenon with dimorphic sexual difference in the body size. For Cynoglossus semilaevis, a flatfish in China, although the importance of chromosome W and its specific gene zbed1 in female-biased SSD have been suggested, its family members and regulation information are still unknown. At present, three zbed1 copies gene were identified on chromosome W, with no gametologs. Phylogenetic analysis for the ZBED family revealed an existence of ZBED9 in the fish. Nine members were uncovered from C. semilaevis, clustering into three kinds, ZBED1, ZBED4 and ZBEDX, which is less than the eleven kinds of ZBED members in mammals. The predominant expression of zbed1 in the female brain and pituitary tissues was further verified by qPCR. Transcription factor c/ebpα could significantly enhance the transcriptional activity of zbed1 promoter, which is opposite to its effect on the male determinant factor-dmrt1. When zbed1 was interfered with, piwil1, esr2 and wnt7b were up-regulated, while cell-cycle-related genes, including cdk4 and ccng1, were down-regulated. Thus, zbed1 is involved in cell proliferation by regulating esr2, piwil1, cell cycle and the Wnt pathway. Further research on their interactions would be helpful to understand fish SSD.

Identification and Functional Analysis of Cynoglossus semilaevis Z-Linked E3 Ubiquitin Ligase rnf34.

The high proportion of males in C. semilaevis hinders their industrial development. The genetic ZW individual can become a pseudomale by sex reversal. And the pseudomale can produce Z-sperm (with epigenetic information to cause sex reversal) while W-sperm is absent, which leads to an even higher male proportion in offspring. Recently, with the development of transcriptomic technologies, research on spermatogenesis in C. semilaevis has been focused on the ubiquitination pathway. In this study, we analyzed the function of the ubiquitin ligase rnf34 gene on the Z chromosome. A qPCR experiment showed that its expression level in the gonad was the highest among different tissues. In the ovary, the expression gradually increased with development from 40 days post-hatching (dph) to 1.5 years post-hatching (yph). In the testis, rnf34 showed increased expression from 40 dph to 6 months post-hatching (mpf) and stabilized up until 1.5 ypf. In situ hybridization showed that the mRNA of rnf34 was mainly distributed in the germ cells of the testis and the ovary. In vivo siRNA-mediated knockdown of the rnf34 gene in male fish affected the expression of a series of genes related to sex differentiation and spermatogenesis. These results provide genetic data on the molecular mechanisms of gonadal development and spermatogenesis in C. semilaevis.

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.

Identification and Functional Analysis of foxo Genes in Chinese Tongue Sole (Cynoglossus semilaevis).

The Chinese tongue sole (Cynoglossus semilaevis) is a traditional, precious fish in China. Due to the large growth difference between males and females, the investigation of their sex determination and differentiation mechanisms receives a great deal of attention. Forkhead Box O (FoxO) plays versatile roles in the regulation of sex differentiation and reproduction. Our recent transcriptomic analysis has shown that foxo genes may participate in the male differentiation and spermatogenesis of Chinese tongue sole. In this study, six Csfoxo members (Csfoxo1a, Csfoxo3a, Csfoxo3b, Csfoxo4, Csfoxo6-like, and Csfoxo1a-like) were identified. Phylogenetic analysis indicated that these six members were clustered into four groups corresponding to their denomination. The expression patterns of the gonads at different developmental stages were further analyzed. All members showed high levels of expression in the early stages (before 6 months post-hatching), and this expression was male-biased. In addition, promoter analysis found that the addition of C/EBPα and c-Jun transcription factors enhanced the transcriptional activities of Csfoxo1a, Csfoxo3a, Csfoxo3b, and Csfoxo4. The siRNA-mediated knockdown of the Csfoxo1a, Csfoxo3a, and Csfoxo3b genes in the testicular cell line of Chinese tongue sole affected the expression of genes related to sex differentiation and spermatogenesis. These results have broadened the understanding of foxo's function and provide valuable data for studying the male differentiation of tongue sole.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis).

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

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

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.

Identification and Expression Pattern of cyp26b1 Gene in Gonad of the Chinese Tongue Sole (Cynoglossus semilaevis).

As an RA-metabolizing enzyme, cyp26b1 has a substantial impact on RA-signaling pathways. The cyp26b1 gene from the Chinese tongue sole was cloned and identified in this investigation. The cyp26b1 ORF was 1536 bp in length and encoded a 512 amino acid protein. A quantitative real-time PCR (qPCR) indicated that the cyp26b1 expression is no significant sexual dimorphism in the gonads at the 80 days post-hatching (dph) stages. After 4 months post-hatching (mph), the expression of cyp26b1 showed sexual dimorphism and lower level of expression in the ovaries than in the testes. An in situ hybridization demonstrated that cyp26b1 mRNA was primarily located in the testis. Interestingly, the cyp26b1 mRNA probe was also detected in the ovaries. These results suggested that cyp26b1 participates in the sex-differentiation and gonadal development of the Chinese tongue sole.

A patched1 gene homologue participates in female differentiation of Cynoglossus semilaevis.

Patched (Ptch) is a receptor in the hedgehog signaling pathway, essential for animal development. Our previous study showed that ptch1 gene participates in the maintenance of the male germline and spermatogenesis in Cynoglossus semilaevis (csptch1). In this study, we identified a patched1 gene homolog (csptch1 x1). The csptch1 x1 gene is 5761 bp long, with a 4638 bp coding sequence that encodes 1545 amino acids. The Csptch1 x1 protein has 12 transmembrane regions and sterol-sensing domains and is highly homologous to the csptch1 (91 amino acids difference). Expression pattern analysis showed that csptch1 x1 is expressed in eight different tissues of adult tongue sole, and the expression is significantly higher in tissues of female than that in male tissues. The expression pattern in developmental stages was also analyzed. csptch1 x1 could be detected at the 1-cell stage and was highly expressed at the blastocyst, somite, and blastopore closing stages, implying that it participates in cell differentiation. In ovarian development, the expression of csptch1 x1 was initiated at 20 days after hatching (dah) and was significantly high at 35-50 and 95-150 dah. In situ hybridization showed that csptch1 x1 was predominantly expressed in primordial germ cells, oocytes, and follicular cells, but the expression of the gene was lower in the testis. These results suggest that csptch1 x1 may be mainly involved in female differentiation and ovarian development, different from the role of csptch1 in spermatogenesis.

Identification of stress-related genes by co-expression network analysis based on the improved turbot genome.

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties.

Transcriptomic analysis reveals the gene expression profiles in the spleen of spotted knifejaw (Oplegnathus punctatus) infected by Vibrio harveyi.

As one of the most valuable maricultured species, spotted knifejaw (Oplegnathus punctatus) has high popularity in eastern Asia. In recent years, diseases caused by Vibrio harveyi have brought huge economic losses in spotted knifejaw industry. To better understand the molecular mechanisms of immune response about V. harveyi resistance in spotted knifejaw, a comparative transcriptome analysis was performed on spleen tissues at five different time points post-infection (0, 12, 24, 48 and 72 hpi). A total of 4279 differentially expressed genes (DEGs) were identified. KEGG pathways analysis showed that multiple immune-related pathways were significant regulated, including Toll-like receptor signaling pathway, ECM-receptor interaction pathway, cytokine-cytokine receptor interaction pathway and hematopoietic cell lineage pathway. Weighted gene co-expression network analysis showed that several immune-related pathways of the highest correlation with 12 hpi (cor = 0.89, P = 7e-06) were significantly enriched. In addition, 12 hpi was a turning point for 7 gene clusters out of 9 that were divided according to gene expression patterns. Therefore, we speculated that 12 hpi might be a very critical time point for spotted knifejaw against V. harveyi infection. Additionally, qRT-PCR was carried out to validate the expressions of 12 DEGs. This study provided the first systematical transcriptome analysis of spotted knifejaw against V. harveyi. The results could help us better understand the dynamic immune responses of spotted knifejaw against bacterial infection, and provide useful information for antibacterial defense in spotted knifejaw industry as well.

pitpß_w Encoding Phosphatidylinositol Transfer Protein Is Involved in Female Differentiation of Chinese Tongue Sole, Cynoglossus semilaevis.

Phosphatidylinositol transfer protein (pitp) plays an important role in phospholipid transfer in animals. A pitp variant (pitpß_w) in Chinese tongue sole was identified by transcriptomic analysis for its female-biased expression. The coding sequence of pitpß_w was 816 bp, encoding a 371-amino-acid protein. pitpß_w showed female-biased expression and was relatively high in brain, muscle, and ovary tissues. In different developmental stages of the ovary, pitpß_w could be detected from 40 days until 3 years post hatching, and the highest expression was observed at 90 days. In situ hybridization revealed that pitpß_w was predominantly localized in early-stage oocytes (I-III stages). After siRNA-mediated knockdown of pitpß_w in an ovarian cell line, the expression of sox9a was reduced, while that of figla_tv1 and sox9b was significantly increased. Our findings suggest that pitpß_w might be involved in female differentiation and early oogenesis.

Genome-wide identification, immune response profile and functional characterization of IL-10 from spotted knifejaw (Oplegnathus punctatus) during host defense against bacterial and viral infection.

Interleukin 10 (IL-10), a pleiotropic cytokine, plays an essential role in multiple immunity responses. In the current study, the sequences of IL-10 family were identified from spotted knifejaw (Oplegnathus punctatus) whole genome, and O. punctatus IL-10 (OpIL-10) was cloned and characterized. OpIL-10 encodes 187 amino acids with a typical IL-10 family signature motif and predicted α-helices. It shared high identities with Notolabrus celidotus IL-10 and Epinephelus Lanceolatus IL-10. OpIL-10 was widely detected in healthy tissues, with the abundant expression in liver and skin. It was significantly up-regulated in the six immune-related tissues (liver, spleen, kidney, intestine, gill and skin) after infection against Vibrio harveyi and spotted knifejaw iridovirus (SKIV). Dual-luciferase analysis showed that OpIL-10 overexpression could suppress the activity of NF-κB. Meanwhile, OpIL-10 knockdown caused the down-regulation of five immune-related genes in JAK2/STAT3 signaling pathway and NF-κB signaling pathway, including IL-10R2, TYK2, STAT3, NOD2, and IκB. In addition, LPS and poly I:C stimulated expression of pro-inflammatory cytokines, including IL-6, IL-1ß, IL-8, and IL-12, were lower with recombinant OpIL-10 (rOp IL-10) than the control group, indicating the anti-inflammatory roles of rOpIL-10. Taken together, these results indicated OpIL-10 as a negative regulator in the inflammatory responses of spotted knifejaw against bacterial and viral infection, which would help us better understand the role of IL-10 in teleost immunity.

Sustainable use of CRISPR/Cas in fish aquaculture: the biosafety perspective.

Aquaculture is becoming the primary source of seafood for human diets, and farmed fish aquaculture is one of its fastest growing sectors. The industry currently faces several challenges including infectious and parasitic diseases, reduced viability, fertility reduction, slow growth, escapee fish and environmental pollution. The commercialization of the growth-enhanced AquAdvantage salmon and the CRISPR/Cas9-developed tilapia (Oreochromis niloticus) proffers genetic engineering and genome editing tools, e.g. CRISPR/Cas, as potential solutions to these challenges. Future traits being developed in different fish species include disease resistance, sterility, and enhanced growth. Despite these notable advances, off-target effect and non-clarification of trait-related genes among other technical challenges hinder full realization of CRISPR/Cas potentials in fish breeding. In addition, current regulatory and risk assessment frameworks are not fit-for purpose regarding the challenges of CRISPR/Cas notwithstanding that public and regulatory acceptance are key to commercialization of products of the new technology. In this study, we discuss how CRISPR/Cas can be used to overcome some of these limitations focusing on diseases and environmental release in farmed fish aquaculture. We further present technical limitations, regulatory and risk assessment challenges of the use of CRISPR/Cas, and proffer research strategies that will provide much-needed data for regulatory decisions, risk assessments, increased public awareness and sustainable applications of CRISPR/Cas in fish aquaculture with emphasis on Atlantic salmon (Salmo salar) breeding.

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.

A Z-Linked E3 Ubiquitin Ligase Cs-rchy1 Is Involved in Gametogenesis in Chinese Tongue Sole, Cynoglossus semilaevis.

Ubiquitin ligase (E3) plays a versatile role in gonadal development and spermatogenesis in mammals, while its function in fish is little reported. In this study, a Z-chromosome linked ubiquitin ligase rchy1 in C. semilaevis (Cs-rchy1) was cloned and characterized. The full-length cDNA was composed of 1962 bp, including 551 bp 5'UTR, 736 bp 3'UTR, and 675 bp ORF encoding a 224-amino-acid (aa) protein. Cs-rchy1 was examined among seven different tissues and found to be predominantly expressed in gonads. In testis, Cs-rchy1 could be detected from 40 days post hatching (dph) until 3 years post hatching (yph), but there was a significant increase at 6 months post hatching (mph). In comparison, the expression levels in ovary were rather stable among different developmental stages. In situ hybridization showed that Cs-rchy1 was mainly localized in germ cells, that is, spermatid and spermatozoa in testis and stage I, II and III oocytes in ovary. In vitro RNA interference found that Cs-rchy1 knockdown resulted in the decline of sox9 and igf1 in ovarian cell line and down-regulation of cyp19a in the testicular cell line. These data suggested that Cs-rchy1 might participate in gonadal differentiation and gametogenesis, via regulating steroid hormone synthesis.

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.

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.

Cloning, expression and functional analysis of the desert hedgehog (dhh) gene in Chinese tongue sole (Cynoglossus semilaevis).

Desert hedgehog (dhh) is a gene that is crucial for spermatogenesis and Leydig cell differentiation, but little is known regarding its influence on gonadal differentiation and development in fish. To understand its function, we cloned and characterized the dhh gene from Cynoglossus semilaevis (csdhh). The full length csdhh cDNA was 2473 bp, including a 1386 bp open reading frame (ORF), a 475 bp 5'-UTR, and a 612 bp 3'-UTR, encoding a predicted protein of 461 amino acid residues. Phylogenetic analysis showed that the putative protein belongs to the hedgehog (HH) family, and contains typical HH-N and HH-C domains. Amino acid sequence analysis revealed that CsDhh shares many features with Dhh analogues in other teleost species. Real-time quantitative PCR showed that csdhh was detected in eight different tissues in male and female tongue sole. During early embryonic development, the relative expression of the csdhh was significantly higher in the neural stage than in other embryonic developmental stages (P < 0.05). csdhh was detected at 20 days after hatching (dah) and at the critical period of male gonadal differentiation (80-95 dah), the relative expression of the csdhh was significantly higher in the male gonads than the female gonads. In 5, 8, and 12 month old gonads, the relative expression of the csdhh was significantly higher in male and pseudo-male than in female fish. The in situ hybridization (ISH) results showed that the hybridization signal was strongly expressed in primary and secondary spermatocytes, spermatids, and sertoli cells of the 1-year-old fish testis, with only weak signal expression in the corresponding ovarian tissue. These results suggest that csdhh is highly conserved in evolution and plays an important role in spermatogenesis in males and pseudo-males.

Molecular characterization and expression pattern of inhibin α and ßb in Chinese tongue sole (Cynoglossus semilaevis).

Inhibin plays important roles in vertebrate reproduction and development. In this study, we have cloned two genes encoding inhibin subunits, inhα and ihnßb, in Chinese tongue sole. inhα consists of 1032 bp, encoding a 343 amino-acid protein. inhßb is composed of 1275 bp, encoding a 424 amino-acid protein. Phylogenetic tree analysis indicated that INHα and INHßB were independently evolved. qPCR showed that inhα expression of in male testis was higher than that in ovary and pseudomale testis, while the expression of inhßb in ovary was higher than that in male and pseudomale testis. During gonadal developmental stages, inhα expression reached highest at 120 days post hatching (dph) both in ovary and testis, then showed decline in ovary but it was first decreased and then increased in the testis. Similarly, inhßb expression in ovary was low at 50-80 dph. At 120 dph, its expression was significantly increased to the peak level, and then gradually decreased. inhßb expression in testis maintained at a low level. During the embryonic developmental stages, inhα displayed the highest expression at 32-cell stage, whereas inhßb reached the highest expression at blastula stages. In situ hybridization data showed that both of inhα and inhßb were detected in oocytes of all stages. In male testis, inhα and inhßb was localized in spermatogonia, spermatocytes, spermatozoa, sertoli and leydig cells. In pseudomale testis, inhα showed the similar pattern in male testis, while the inhßb was detected in spermatocytes and spermatozoa. These data suggested that inhα may participate the spermatogenesis and oogenesis of Chinese tongue sole, while inhßb might predominantly function in oogenesis.