Comparative transcriptome analysis of ovary and testis reveals potential sex-related genes and pathways in spotted knifejaw Oplegnathus punctatus

Exploitation and Application of a New Genetic Sex Marker Based on Intron Insertion Variation of erc2 Gene in Oplegnathus punctatus

Spotted knifejaw (Oplegnathus punctatus), one of the most valuable mariculture species, grows with significant sexual dimorphism, with males growing significantly faster than females. O. punctatus not only has excellent growth characteristics and high food value, but also shows high economic value in aquaculture, which has become a hotspot in the field of aquaculture. The current insufficiency of sex marker identification in O. punctatus restricts the process of its unisexual breeding. Rapid identification of sex will help to study the mechanisms of sex determination and accelerate the development of sex-controlled breeding. With the completion of the sequencing of the male and female genomes of O. punctatus, the efficient and precise development of genetic sex markers has been made possible. In this study, we used genome-wide information combined with molecular biology techniques from marker sequences to further establish a rapid method for DNA insertion variant detection in the intron of O. punctatus erc2 gene, which can be used to rapidly, accurately, and efficiently identify whether DNA insertion occurs in the intron of O. punctatus erc2 gene to be detected, and to identify the sex of O. punctatus to be detected. It could also be distinguished by agarose gel electrophoresis, which would shorten the time for accurate identification and improves the detection efficiency. Homozygous comparison of male and female individuals showed that the length of the DNA fragment of the erc2 gene was 239 bp on chromosome X1 and 1173 bp on chromosome Y. It can therefore be inferred that a 934 bp insertion fragment exists on the Y chromosome. The PCR amplification results showed that two DNA fragments of 1173 bp and 239 bp could be amplified in male O. punctatus, and the 1173 bp fragment was a marker fragment specific to the variant intron erc2 gene, while only a single DNA fragment of 239 bp was amplified in female O. punctatus. It has important significance and application value in the study of neurotransmitter transmission and environmental adaptability of female and male fish based on erc2 gene, as well as the identification of male and female sex, the preparation of high male fry, and family breeding.

De Novo Assembly, Characterization, and Comparative Transcriptome Analysis of Mature Male and Female Gonads of Rabbitfish (Siganus oramin) (Bloch & Schneider, 1801)

The rabbitfish, Siganus oramin, is a commercially important table fish in southeastern China. However, there have been few studies on its gonad development and reproduction regulation. Comparative transcriptome analysis was first performed on adult male and female gonads of S. oramin. In total, 47,070 unigenes were successfully assembled and 22,737 unigenes were successfully annotated. Through comparative transcriptome analysis of male and female gonads, a total of 6722 differentially expressed genes were successfully identified, with 3528 upregulated genes and 3154 downregulated genes in the testes. In addition, 39 differentially expressed reproduction-related genes were identified. Finally, quantitative real-time PCR was used to validate the expression levels of several differentially expressed genes. These results provide important data for further studying the function of reproduction-related genes and the molecular mechanism regulating gonad development and reproduction in S. oramin.

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

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

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.

Expression pattern and cellular localization of pepsinogen in early development and induced by different diets in the spotted knifejaw (Oplegnathus punctatus)

To solve the high mortality rate of early-stage larval feed conversion during aquaculture in Oplegnathus punctatus, the investigation of the structural and functional characteristics of the gastric tissue was conducted. Histological results showed that the gastric gland rudiment appeared at 17 dph. The basic structure of the stomach was fully developed between 26 and 35 dph. Two pepsinogen genes, named OpPGA1 and OpPGA2, were identified in the spotted knifejaw genome. qPCR results of developmental period showed that the two genes were low in expression during early development (5 and 15 dph). At 20 dph, the two genes started to show trace expression, and at 30 dph the mRNA expression levels of OpPGA1 and OpPGA2 reached the highest levels. Results of pepsin activity detection during the development period showed lower activity was detected 22 dph, followed by a peak at 30 dph. Under different feeding inductions, OpPGA1 showed the highest expression in the basic diet group and hard-shell group, while the expression level in the phytophagous group remained consistently low. The mRNA expression level of OpPGA2 in the phytophagous group was significantly higher than in other groups. Enzyme activity determination under different feeding inductions showed slightly higher enzyme activity in the basic diet group and crustacean group. The results of in situ hybridization showed that the mRNA of both OpPGA1 and OpPGA2 genes was both expressed in gastric gland cells. These information can contribute to the development of practical feeding methods in terms of digestive physiology for the development of larvae.

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

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

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

Comparative Studies on Duplicated foxl2 Paralogs in Spotted Knifejaw Oplegnathus punctatus Show Functional Diversification

As a member of the forkhead box L gene family, foxl2 plays a significant role in gonadal development and the regulation of reproduction. During the evolution of deuterostome, whole genome duplication (WGD)-enriched lineage diversifications and regulation mechanisms occurs. However, only limited research exists on foxl2 duplication in teleost or other vertebrate species. In this study, two foxl2 paralogs, foxl2 and foxl2l, were identified in the transcriptome of spotted knifejaw (Oplegnathus punctatus), which had varying expressions in the gonads. The foxl2 was expressed higher in the ovary, while foxl2l was expressed higher in the testis. Phylogenetic reconstruction, synteny analysis, and the molecular evolution test confirmed that foxl2 and foxl2l likely originated from the first two WGD. The expression patterns test using qRT-PCR and ISH as well as motif scan analysis revealed evidence of potentially functional divergence between the foxl2 and foxl2l paralogs in spotted knifejaw. Our results indicate that foxl2 and foxl2l may originate from the first two WGD, be active in transcription, and have undergone functional divergence. These results shed new light on the evolutionary trajectories of foxl2 and foxl2l and highlights the need for further detailed functional analysis of these two duplicated paralogs.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Tea polyphenols, astaxanthin, and melittin can significantly enhance the immune response of juvenile spotted knifejaw (Oplegnathus punctatus)

The frequent occurrence of diseases seriously hampers the sustainable development of the spotted knifejaw (Oplegnathus punctatus) breeding industry. Our previous genome-wide scan and cross-species comparative genomic analysis revealed that the immune gene family (Toll-like receptors, TLR) members of O. punctatus underwent a significant contraction event (tlr1, tlr2, tlr14, tlr5, and tlr23). To address immune genetic contraction may result in reduced immunity, we investigated whether adding different doses (0, 200, 400, 600, and 800 mg/kg) of immune enhancers (tea polyphenols, astaxanthin, and melittin) to the bait after 30 days of continuous feeding could stimulate the immune response of O. punctatus. We found that the expression of tlr1, tlr14, tlr23 genes in immune organs (spleen and head kidney) was stimulated when tea polyphenols were added at 600 mg/kg. The tlr2 (400 mg/kg), tlr14 (200 mg/kg), tlr5 (200 mg/kg), and tlr23 (200 mg/kg) genes expression of intestine were elevated in the tea polyphenol group. When the addition of astaxanthin is 600 mg/kg, it can effectively stimulate the expression of tlr14 gene in immune organs (liver, spleen and head kidney). In the astaxanthin group, the expression of the genes tlr1 (400 mg/kg), tlr14 (600 mg/kg), tlr5 (400 mg/kg) and tlr23 (400 mg/kg) reached their highest expression in the intestine. Besides, the addition of 400 mg/kg of melittin can effectively induce the expression of tlr genes in the liver, spleen and head kidney, except the tlr5 gene. The tlr-related genes expression in the intestine was not significantly elevated in the melittin group. We hypothesize that the immune enhancers could enhance the immunity of O. punctatus by increasing the expression of tlr genes, and thereby leading to increased resistance to diseases. Meanwhile, our findings further demonstrated that significant increases in weight gain rate (WGR), visceral index (VSI), and feed conversion rate (FCR) were observed at 400 mg/kg, 200 mg/kg and 200 mg/kg of tea polyphenols, astaxanthin and melittin in the diet, respectively. Overall, our study provided valuable insights for future immunity enhancement and viral infection prevention in O. punctatus, as well as offered guidance for the healthy development of the O. punctatus breeding industry.

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

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

Comparative transcriptome sequencing analysis of female and male Decapterus macrosoma

Sexual growth dimorphism is a common phenomenon in teleost fish and has led to many reproductive strategies. Growth- and sex-related gene research in teleost fish would broaden our understanding of the process. In this study, transcriptome sequencing of shortfin scad Decapterus macrosoma was performed for the first time, and a high-quality reference transcriptome was constructed. After identification and assembly, a total of 58,475 nonredundant unigenes were obtained with an N50 length of 2,266 bp, and 28,174 unigenes were successfully annotated with multiple public databases. BUSCO analysis determined a level of 92.9% completeness for the assembled transcriptome. Gene expression analysis revealed 2,345 differentially expressed genes (DEGs) in the female and male D. macrosoma, 1,150 of which were female-biased DEGs, and 1,195 unigenes were male-biased DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the DEGs were mainly involved in biological processes including protein synthesis, growth, rhythmic processes, immune defense, and vitellogenesis. Then, we identified many growth- and sex-related genes, including Igf, Fabps, EF-hand family genes, Zp3, Zp4 and Vg. In addition, a total of 19,573 simple sequence repeats (SSRs) were screened and identified from the transcriptome sequences. The results of this study can provide valuable information on growth- and sex-related genes and facilitate further exploration of the molecular mechanism of sexual growth dimorphism.

Comparative transcriptome provides insights into differentially expressed genes between testis and ovary of Onychostoma macrolepis in reproduction period

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered fishery species inhabiting the river of Qinling Mountains and some flowing freshwaters in China. The declining population of O. macrolepis caused by asynchrony of male and female development prompted us to focus on genetic regulation of its reproduction. In this study, high-throughput RNA-sequencing technology was applied to assemble and annotate the transcriptome of O. macrolepis testis and ovary. The results showed that a number of 338089335 (ovary:163216500, testis:174872835) raw sequences were obtained. After non-redundant analysis, a number of 207826065 (ovary:102334008, testis:105492057) high quality reads were obtained and predicted as unigenes, in which 201,038,682 unigenes were annotated with multiple databases. Taking the ovarian transcriptome as a control, comparative transcriptome analysis showed that 9918 differentially expressed genes (DEGs) up-regulated in the testis and 13,095 DEGs down-regulated. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, gonadal development, steroid biosynthesis pathways, MAPK signaling pathway and Wnt signaling pathway. Finally, the expression patterns of 19 unigenes were validated by using quantitative real-time polymerase chain reaction (qRT-PCR). This study illustrates a potential molecular mechanism on the unsynchronized male and female development of the O. macrolepis during the reproduction period in June and provides a theoretical basis for future artificial reproduction.

Small RNA sequencing reveals sex-related miRNAs in Collichthys lucidus

Collichthys lucidus (C. lucidus) is an economically important fish species, exhibiting sexual dimorphism in its growth rate. However, there is a lack of research on its underlying sex-related mechanisms. Therefore, small RNA sequencing was performed to better comprehend these sex-related molecular mechanisms. In total, 171 differentially expressed miRNAs (DE-miRNAs) were identified between the ovaries and testes. Functional enrichment analysis revealed that the target genes of DE-miRNAs were considerably enriched in the p53 signaling, PI3K–Akt signaling, and TGF-beta signaling pathways. In addition, sex-related miRNAs were identified, and the expression of miR-430c-3p and miR-430f-3p was specifically observed in the gonads compared with other organs and their expression was markedly upregulated in the testes relative to the ovaries. Bmp15 was a target of miR-430c-3p and was greatly expressed in the ovaries compared with the testes. Importantly, miR-430c-3p and bmp15 co-expressed in the ovaries and testes. This research provides the first detailed miRNA profiles for C. lucidus concerning sex, likely laying the basis for further studies on sex differentiation in C. lucidus.

Identification of Male-Specific Molecular Markers by Recombination of RhoGEF10 Gene in Spotted Knifejaw (Oplegnathus punctatus)

The spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high ecological value, food value, and fishing value, and its growth has obvious sex dimorphism. The rapid identification of its sex is beneficial to the development of sex determination and breeding. In this study, the method of comparative genomics and PCR amplification was used to further establish a rapid detection method for the recombinant RhoGEF10 gene in O. punctatus, which can quickly, accurately, and efficiently identify the sex of the O. punctatus to be tested. The homologous comparison results of male and female individuals showed that the DNA fragment length of the RhoGEF10 gene on the X1 chromosome was 326 bp, and the DNA fragment length on the Y chromosome was 879 bp. Therefore, it can be concluded that there is an insert fragment of 553 bp on the Y chromosome. PCR amplification results showed that the two DNA fragments of 879 bp and 326 bp were amplified in the Y chromosome and X1 chromosome of the male O. punctatus (X1X2Y), respectively, and the 879 bp fragment was a unique marker fragment of the recombinant RhoGEF10 gene; The female O. punctatus (X1X1X2X2) only a single DNA fragment of 326 bp was amplified. At the same time, the inserted fragment of the male individual resulted in partial inactivation of the RhoGEF10 protein, which in turn resulted in a slowing of peripheral nerve conduction velocity and thinning of the myelin sheath in male O. punctatus. The method shortens the time for accurate identification of the O. punctatus RhoGEF10 gene recombination and improves the detection efficiency. It is of great significance and application value in the research of nerve conduction and myelin development, male and female sex identification, the preparation of high male seedlings, and family selection based on the RhoGEF10 gene in the O. punctatus.

Comparative Transcriptome Analyses Characterize Expression Signatures Among Males, Females, Neo-Males, and Gynogenetic Females in the Yellow Drum (Nibea albiflora)

The yellow drum (Nibea albiflora) is one of the most important marine economic fish in China, and its sexually dimorphic growth makes it preferable for mono-sex culture. Although gynogenesis and neo-male induction techniques have been established, the molecular pathways and regulatory mechanisms of sex determination and maintenance in gynogenetic females and neo-males remains far from fully understood. In this study, the gene expression profiles were investigated in the gonads and brains of wild-type male, wild-type female, neo-male, and gynogenetic female yellow drum using comparative transcriptome analyses. Generally, a total of 52,999 novel transcripts were obtained in RNA-seq, of which 45,651 were isoforms of known protein-coding genes, 1,358 novel protein-coding genes, and 5,990 long non-coding RNAs. We found that the differences between wild-type males and neo-males and between wild-type females and gynogenetic females were relatively small at both the histological and transcriptomic levels, indicating that artificial gynogenesis or hormonal sex reversal may have minimal effects on normal female or male life function, respectively. In the brain, pathways such as “Oocyte meiosis”, “Cell cycle”, and “Riboflavin metabolism” were found to be significantly enriched. In the gonads, pathways such as “Prolactin signaling pathway”, “PPAR signaling pathway”, “Cholesterol metabolism”, and “Jak-STAT signaling pathway” were found to play important roles in maintaining the regular proliferation and differentiation of females and males in yellow drum. In particular, we found that zp4 might be an effective molecular marker to differentiate between gynogenetic and normal females owing to its unique expression pattern. The results of this study may help to elucidate the molecular mechanisms involved in sex maintenance in the gonads and brain and provide basic data for genetic breeding of the yellow drum.

Transcription Analysis for Core Networks of lncRNAs–mRNAs: Implication for Potential Role in Sterility of Crassostrea gigas

Simple Summary

This study reveals the expression profiles of lncRNA in the gonads of the Pacific oyster Crassostrea gigas. The potential function of lncRNAs was predicted in the case of antisense and cis-regulatory mechanisms based on their physical positions and their coexpression relationships in the case of trans regulation. Sterility-related DEGs and DELs were chosen for subsequent analysis, demonstrating that trans-regulatory lncRNAs might play a vital role in the gametogenesis of C. gigas. We constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites based on pathway results, in which 28 lncRNAs and their 54 trans-regulatory genes were detected. Among 28 sterility-specific lncRNAs, MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites play the most significant role.

Abstract

Long noncoding RNA (lncRNA), a type of non-protein-coding transcript, is emerging as a crucial regulator of gene expression. However, few roles of lncRNA in the reproductive process of the Pacific oyster (Crassostrea gigas) have been defined, especially in the regulatory mechanism of sterile triploids gametogenesis. To uncover the potential role of lncRNA, the gonads of diploids, sterile triploids, and partially sterile triploids underwent RNA sequencing. A total of 9618 reliable lncRNAs were identified. The target relationship between lncRNA and mRNA was predicted based on cis, trans, and antisense regulation with bioinformatic software. We chose differentially expressed lncRNAs and mRNAs when sterile triploids were compared to partially sterile triploids and diploids for subsequent functional enrichment analysis. Findings revealed that trans-regulatory lncRNAs might play a significant role in the gametogenesis of C. gigas. Combining pathway results, we constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites. Fifty-four genes related to cell division, germline-cell maintenance, and glycogen metabolism were found to be associated with sterility. A total of 28 candidate lncRNAs were predicted to trans-regulate these genes. We speculated that MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites were highly important as they were predicted to regulate more sterility-specific genes than others. Our work collectively identified sterility-related lncRNAs and implicated the potential mechanism of lncRNA-mediated regulation in the gametogenesis of sterile triploid oysters.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

Comparative transcriptomic analysis of gonadal development and renewal in the ovoviviparous black rockfish (Sebastes schlegelii)

BACKGROUND

The black rockfish (Sebastes schlegelii) has an ovoviviparous reproductive pattern and long-term sperm storage, resulting in asynchronous gonadal development between the sexes. However, the comprehensive understanding of gonadal development in black rockfish has not yet been achieved. Here, we studied gonadal development and germ cell renewal using histology and RNA-seq.

RESULTS

In this study, RNA-seq was performed on testes and ovaries to characterize key pathways and genes that are active during development and gamete maturation in black rockfish. Differentially expressed genes (DEGs) were identified and annotated in 4 comparisons (F_III vs. F_IV, F_IV vs. F_V, M_III vs. M_IV and M_IV vs. M_V). Based on analysis of DEGs enriched in the testis, 11 and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped to the M_III vs. M_IV group and the M_IV vs. M_V group, respectively. DEGs in ovarian development were also classified into 10 groups according to their biological functions. The expression patterns of the selected genes determined by qPCR were significantly correlated with the RNA-Seq results, supporting the reliability and accuracy of the RNA-Seq analysis. E2 levels showed down regulation from previtellogenesis to mature stage in female and T level showed down regulation from spermatogenesis to regressed stage in the male.

CONCLUSIONS

The categories "intercellular interaction and cytoskeleton", "molecule amplification" and "repair in the cell cycle" were revealed to be crucial in testis development and spermatogenesis, as was the biosynthesis of a series of metabolites. Our results provide comprehensive insight into black rockfish gonadal development and provide a basis for further study of reproductive physiology and molecular biology in ovoviviparity teleosts.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

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

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

Genome-Wide DNA Methylation and Transcriptome Analyses Reveal Epigenetic and Genetic Mechanisms Underlying Sex Maintenance of Adult Chinese Alligator

The sexes of Chinese alligators are determined during embryonic development and remain fixed thereafter. In this study, we investigated the genetic and epigenetic mechanisms underlying sex maintenance in Chinese alligators through RNA sequencing and bisulfite sequencing data analyses of the adult gonads. We identified the genes and pathways (e. g., DMRT1-SOX9-AMH pathway for males and oocyte meiotic maturation pathway for females) involved in male and female sex maintenance and gonadal development of adult Chinese alligators. In contrast to their expression patterns in the embryo, both DMRT1 and the steroid hormone biosynthesis related genes showed a male-biased expression in adult gonads. The overall DNA methylation density and level were higher in testes than in ovaries. Hypermethylation in the gene bodies enhanced the expression of male-biased genes (such as DMRT1-SOX9-AMH and steroid hormone biosynthesis related genes) in the testis, as opposed to the normalization of gene expression. Our results provide insights into the genetic and epigenetic mechanisms underlying sex maintenance in adult Chinese alligators, and are expected to contribute to the development of scientific programs for the successful conservation of this endangered species.

Application of quantitative transcriptomics in evaluating the ex vivo effects of per- and polyfluoroalkyl substances on Atlantic cod (Gadus morhua) ovarian physiology

Because of their global consumption and persistence, per- and polyfluoroalkyl substances (PFASs), are ubiquitously distributed in the environment, as well as in wildlife and humans. In the present study, we have employed an ex vivo organ culture technique, based on the floating agarose method, of Atlantic cod ovarian tissue to investigate the effects of three different concentrations of PFOS, PFOA (1, 5 and 25 μM) and PFNA (0.5, 5 and 50 μM), used singly and in also in combination (1×, 20× and 100×). In the 1× exposure mixture, concentrations were decided based on their proportional levels (in molar equivalents) relative to PFOS, which is the most abundant PFAS in cod liver from a 2013 screening project. To investigate the detailed underlying mechanisms and biological processes, transcriptome sequencing was performed on exposed ovarian tissue. The number of differentially expressed genes (DEGs) having at least 0.75 log₂-fold change was elevated in high, compared to low and medium concentration exposures. The highest PFNA, PFOA and PFOS concentrations, and the highest (100×) mixture exposure, showed 40, 68, 1295, and 802 DEGs, respectively. The latter two exposure groups shared a maximum of 438 DEGs. In addition, they both shared the majority of functionally enriched pathways belonging to biological processes such as cellular signaling, cell adhesion, lipid metabolism, immunological responses, cancer, reproduction and metabolism. Shortlisted DEGs that were specifically annotated to reproduction associated gene ontology (GO) terms were observed only in the highest PFOS and mixture exposure groups. These transcripts contributed to ovarian key events such as steroidogenesis (star, cyp19a1a), oocyte growth (amh), maturation (igfbp5b, tgfβ2, tgfβ3), and ovulation (pgr, mmp2). Contrary to other PFAS congeners, the highest PFOS concentration showed almost similar transcript expression patterns compared to the highest mixture exposure group. This indicates that PFOS is the active component of the mixture that significantly altered the normal functioning of female gonads, and possibly leading to serious reproductive consequences in teleosts.

Comparative transcriptomic analysis for identification of candidate sex-related genes and pathways in Crimson seabream (Parargyrops edita)

Teleost fishes display the largest array of sex-determining systems among animals, resulting in various reproductive strategies. Research on sex-related genes in teleosts will broaden our understanding of the process, and provide important insight into the plasticity of the sex determination process in vertebrates in general. Crimson seabream (Parargyrops edita Tanaka, 1916) is one of the most valuable and abundant fish resources throughout Asia. However, little genomic information on P. edita is available. In the present study, the transcriptomes of male and female P. edita were sequenced with RNA-seq technology. A total of 388,683,472 reads were generated from the libraries. After filtering and assembling, a total of 79,775 non redundant unigenes were obtained with an N50 of 2,921 bp. The unigenes were annotated with multiple public databases, including NT (53,556, 67.13%), NR (54,092, 67.81%), Swiss-Prot (45,265, 56.74%), KOG (41,274, 51.74%), KEGG (46,302, 58.04%), and GO (11,056, 13.86%) databases. Comparison of the unigenes of different sexes of P. edita revealed that 11,676 unigenes (9,335 in females, 2,341 in males) were differentially expressed between males and females. Of these, 5,463 were specifically expressed in females, and 1,134 were specifically expressed in males. In addition, the expression levels of ten unigenes were confirmed to validate the transcriptomic data by qRT-PCR. Moreover, 34,473 simple sequence repeats (SSRs) were identified in SSR-containing sequences, and 50 loci were randomly selected for primer development. Of these, 36 loci were successfully amplified, and 19 loci were polymorphic. Finally, our comparative analysis identified many sex-related genes (zps, amh, gsdf, sox4, cyp19a, etc.) and pathways (MAPK signaling pathway, p53 signaling pathway, etc.) of P. edita. This informative transcriptomic analysis provides valuable data to increase genomic resources of P. edita. The results will be useful for clarifying the molecular mechanism of sex determination and for future functional analyses of sex-associated genes.

Antibacterial functions of a novel fish-egg lectin from spotted knifejaw (Oplegnathus punctatus) during host defense immune responses

Fish-egg lectins (FELs) have been identified in several teleost species and have been proved to play important roles in innate immune system against pathogen infection. In this study a novel fish-egg lectin (OppFEL) was identified from spotted knifejaw (Oplegnathus punctatus), and the expression patterns against bacterial infection was characterized. The amino acid sequence is highly homologous to other teleost FELs, containing five repeats of the conserved TECPR domain. Expression of OppFEL was widely observed in examined tissues, with the most abundant transcripts observed in gill, showing a pattern of tissue specific expression. The OppFEL expression was significantly up-regulated following a Gram-negative bacterium (Vibrio anguillarum) challenge in vivo, suggesting participation in host antibacterial immune responses. Recombinant OppFEL protein (rOppFEL) possessed calcium dependent binding capacities and agglutination to four Gram-negative bacterium and two Gram-positive bacterium. Sugar binding assay revealed that rOppFEL specifically bound to insoluble lipopolysaccharide and peptidoglycan. In addition, rOppFEL was also proved to have hemagglutinating activity against erythrocytes from Mus musculus, O. punctatus, Sebastes schlegelii and Paralichthys olivaceus. Dual-luciferase analysis showed that overexpression of OppFEL could suppress the activity of NF-κB in a dose dependent manner. Taken together, these results suggest that OppFEL is a unique fish-egg lectin that possesses apparent immunomodulating property and is involved in host defense against pathogens invasion.

Isolation of a Male-Specific Molecular Marker and Development of a Genetic Sex Identification Technique in Spotted Knifejaw (Oplegnathus punctatus)

Sex-specific DNA markers are very helpful for identifying genetic sex and studying sex determination mechanisms in fish. To identify the sex-specific markers of spotted knifejaw (Oplegnathus punctatus), we performed a comparative analysis of the female and male genomes. In this study, an 18 bp insertion was identified in the male genome after verification by sequencing depth and PCR. An effective and rapid method based on PCR was then developed to identify the genetic sex. A male-female-shared primer pair and a male-specific primer were designed for PCR amplification to avoid false-negative phenomena. To examine the primers in practice, we utilized hundreds of spotted knifejaw fish from different groups to identify their genetic sex, and the results were consistent with their phenotypic sex. The male-specific DNA marker would be helpful for artificial breeding, Y chromosome assembly and further study of the sex determination mechanism. This study is the first to identify an effective sex-specific marker in spotted knifejaw.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

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

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

Transcriptome Sequencing Reveals the Traits of Spermatogenesis and Testicular Development in Large Yellow Croaker (Larimichthys crocea)

Larimichthys crocea is an economically important marine fish in China. To date, the molecular mechanisms underlying testicular development and spermatogenesis in L. crocea have not been thoroughly elucidated. In this study, we conducted a comparative transcriptome analysis between testes (TES) and pooled multiple tissues (PMT) (liver, spleen, heart, and kidney) from six male individuals. More than 54 million clean reads were yielded from TES and PMT libraries. After mapping to the draft genome of L. crocea, we acquired 25,787 genes from the transcriptome dataset. Expression analyses identified a total of 3853 differentially expressed genes (DEGs), including 2194 testes-biased genes (highly expressed in the TES) and 1659 somatic-biased genes (highly expressed in the PMT). The dataset was further annotated by blasting with multi-databases. Functional genes and enrichment pathways involved in spermatogenesis and testicular development were analyzed, such as the neuroactive ligand-receptor interaction pathway, gonadotropin-releasing hormone (GnRH) and mitogen-activated protein kinase (MAPK) signaling pathways, cell cycle pathway, and dynein, kinesin, myosin, actin, heat shock protein (hsp), synaptonemal complex protein 2 (sycp2), doublesex- and mab-3-related transcription factor 1 (dmrt1), spermatogenesis-associated genes (spata), DEAD-Box Helicases (ddx), tudor domain-containing protein (tdrd), and piwi genes. The candidate genes identified by this study lay the foundation for further studies into the molecular mechanisms underlying testicular development and spermatogenesis in L. crocea.

A novel C-type lectin from spotted knifejaw, Oplegnathus punctatus possesses antibacterial and anti-inflammatory activity

C-type lectin is a type of carbohydrate-binding protein and plays significant roles in innate immune response against pathogen infection. To date, thousands of C-type lectin had been identified in teleost. In the present study, we isolated a novel isoform of C-type lectin (OppCTL) from spotted knifejaw (Oplegnathus punctatus). The OppCTL encoded a typical Ca²⁺-dependent carbohydrate-binding protein, and was mainly expressed in liver in a tissue specific fashion. The expression of OppCTL was significantly up-regulated following Vibrio anguillarum infection in vivo, suggesting involvement in immune response. Hemagglutination analysis showed that the recombinant OppCTL (rOppCTL) could agglutinate erythrocyte from Mus musculus, Oplegnathus punctatus, Sebastes schlegelii and Paralichthys olivaceus. The rOppCTL could bind and agglutinate all tested bacteria. The rOppCTL possessed capacities of calcium-dependent agglutination to all tested bacteria. Sugar binding assay revealed that rOppCTL could also bind to the glycoconjugates of the bacterial surface, including lipopolysaccharide and peptidoglycan. Interestingly, Dual-luciferase analysis revealed that OppCTL could inhibit the activity of NF-κB in HEK-293T cells after OppCTL overexpression. Taken together, these results indicate that OppCTL has immune activity capable of defending invading pathogens and possesses potential immunoregulatory activity, enriching our understanding of the function of C-type lectin.

Characterization of the European Sea Bass (Dicentrarchus labrax) Gonadal Transcriptome During Sexual Development

The European sea bass is one of the most important cultured fish in Europe and has a marked sexual growth dimorphism in favor of females. It is a gonochoristic species with polygenic sex determination, where a combination between still undifferentiated genetic factors and environmental temperature determines sex ratios. The molecular mechanisms responsible for gonadal sex differentiation are still unknown. Here, we sampled fish during the gonadal developmental period (110 to 350 days post fertilization, dpf), and performed a comprehensive transcriptomic study by using a species-specific microarray. This analysis uncovered sex-specific gonadal transcriptomic profiles at each stage of development, identifying larger number of differentially expressed genes in ovaries when compared to testis. The expression patterns of 54 reproduction-related genes were analyzed. We found that hsd17β10 is a reliable marker of early ovarian differentiation. Further, three genes, pdgfb, snx1, and nfy, not previously related to fish sex differentiation, were tightly associated with testis development in the sea bass. Regarding signaling pathways, lysine degradation, bladder cancer, and NOD-like receptor signaling were enriched for ovarian development while eight pathways including basal transcription factors and steroid biosynthesis were enriched for testis development. Analysis of the transcription factor abundance showed an earlier increase in females than in males. Our results show that, although many players in the sex differentiation pathways are conserved among species, there are peculiarities in gene expression worth exploring. The genes identified in this study illustrate the diversity of players involved in fish sex differentiation and can become potential biomarkers for the management of sex ratios in the European sea bass and perhaps other cultured species.

Transcriptome sequencing and comparative analysis of adult ovary and testis identify potential gonadal maintenance-related genes in Mauremys reevesii with temperature-dependent sex determination

Mauremys reevesii is a classical organism with temperature-dependent sex determination (TSD). Gonad development in early life has recently received considerable attention but gonadal maintenance after sex differentiation in turtles with TSD remains a mystery. In this study, we sequenced the transcriptomes for the adult testis and ovary using RNA-seq, and 36,221 transcripts were identified. In total, 1,594 differentially expressed genes (DEGs) were identified where 756 DEGs were upregulated in the testis and 838 DEGs were upregulated in the ovary. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis suggested that the TGF-beta signaling pathway and Hedgehog signaling pathway have important roles in testis maintenance and spermatogenesis, whereas the Hippo signaling pathway and Wnt signaling pathway are likely to participate in ovary maintenance. We determined the existence of antagonistic networks containing significant specific-expressed genes and pathways related to gonadal maintenance and gametogenesis in the adult gonads of M. reevesii. The candidate gene Fibronectin type 3 and ankyrin repeat domains 1 (FANK1) might be involved with the regulation of testis spermatogenesis.

Cell Metabolism Control Through O-GlcNAcylation of STAT5: A Full or Empty Fuel Tank Makes a Big Difference for Cancer Cell Growth and Survival

O-GlcNAcylation is a post-translational modification that influences tyrosine phosphorylation in healthy and malignant cells. O-GlcNAc is a product of the hexosamine biosynthetic pathway, a side pathway of glucose metabolism. It is essential for cell survival and proper gene regulation, mirroring the metabolic status of a cell. STAT3 and STAT5 proteins are essential transcription factors that can act in a mutational context-dependent manner as oncogenes or tumor suppressors. They regulate gene expression for vital processes such as cell differentiation, survival, or growth, and are also critically involved in metabolic control. The role of STAT3/5 proteins in metabolic processes is partly independent of their transcriptional regulatory role, but is still poorly understood. Interestingly, STAT3 and STAT5 are modified by O-GlcNAc in response to the metabolic status of the cell. Here, we discuss and summarize evidence of O-GlcNAcylation-regulating STAT function, focusing in particular on hyperactive STAT5A transplant studies in the hematopoietic system. We emphasize that a single O-GlcNAc modification is essential to promote development of neoplastic cell growth through enhancing STAT5A tyrosine phosphorylation. Inhibition of O-GlcNAcylation of STAT5A on threonine 92 lowers tyrosine phosphorylation of oncogenic STAT5A and ablates malignant transformation. We conclude on strategies for new therapeutic options to block O-GlcNAcylation in combination with tyrosine kinase inhibitors to target neoplastic cancer cell growth and survival.

Transcriptome Analysis of Male and Female Mature Gonads of Silver Sillago (Sillago sihama)

Silver sillago (Sillago sihama) is an emerging commercial marine aquaculture species in China. To date, fundamental information on S. sihama, such as genomic information, is lacking, and no data are available on the gonad transcriptome of S. sihama. Here, the first gonadal transcriptomes of S. sihama have been constructed and genes potentially involved in gonadal development and reproduction identified. Illumina sequencing generated 60.18 million clean reads for the testis and 59.10 million for the ovary. All reads were assembled into 74,038 unigenes with a mean length of 1,004 bp and N50 value of 2,190 bp. Among all the predictable unigenes, a total of 34,104 unigenes (46%) were searched against multiple databases, including 33,244 unigenes annotated in the RefSeq Non- Redundant database at NCBI, and 28,924 in Swiss-Prot. By comparing the ovary and testis, 35,367 unigenes were identified as being differentially expressed between males and females, of which 29,127 were upregulated in the testis and 6,240 were upregulated in the ovary. Numerous differentially expressed genes (DEGs) known to be involved in gonadal development and gametogenesis were identified, including amh, dmrt1, gsdf, cyp19a1a, gnrhr, and zps. Using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, the top 20 KEGG pathways with highest number of DEGs were found to be involved in regulating gonadal development and gametogenesis in S. sihama. Moreover, 22,666 simple sequence repeats (SSRs) were identified in 14,577 SSR-containing sequences. The findings provide a valuable dataset for future functional analyses of sex-associated genes and molecular marker assisted selection in S. sihama.

Discovery and functional characterization of microRNAs and their potential roles for gonadal development in spotted knifejaw, Oplegnathus punctatus

The spotted knifejaw (Oplegnathus punctatus) is a newly emerging economical fishery species in China. Studies focused on the regulation of gonadal development and gametogenesis of spotted knifejaw are still insufficient. As a key post-transcriptional regulator, miRNAs have been shown to play important roles in development and reproduction systems. In this study, small RNA deep sequencing in ovary and testis of spotted knifejaw were performed to screen miRNA expression patterns. After sequencing and bioinformatics analysis, a total of 247 conserved known miRNAs and 41 novel miRNAs were identified in spotted knifejaw gonads for the first time. In addition, 36 miRNAs were differentially expressed between testis and ovary. The putative target genes of differentially expressed (DE) miRNAs were significantly enriched in several pathways related to sexual differentiation and gonadal development, such as steroid hormone biosynthesis. Sequencing data was validated through qRT-PCR analysis of selected DE miRNAs. Dual-luciferase reporter analyses of filtered miRNA-target gene pairs confirmed that opu-miR-27b-3p targeted in piwi2 and mov10l1 3' UTRs and down-regulated their expressions in spotted knifejaw. The notion that mov10l1 and piwi2 enhance germ cells proliferation and regulate gonadal development and gametogenesis suggests that opu-miR-27b-3p may attenuated this process in the gonads of spotted knifejaw. These findings provided insights into regulatory roles of gonadal miRNAs and supplied fundamental resources for further studies on miRNA-mediated post-transcriptional regulation in reproductive system of spotted knifejaw.