Transcriptome analysis of the gonads of olive flounder (Paralichthys olivaceus)

Global expression pattern of genes containing positively selected sites in European anchovy (Engraulis encrasicolus L.) may shed light on teleost reproduction

European anchovy is a multiple-spawning and highly fecundate pelagic fish with high economic and ecological significance. Although fecundity is influenced by nutrition, temperature and weight of spawners, high reproductive capacity is related to molecular processes in the ovary. The ovary is an essential and complex reproductive organ composed of various somatic and germ cells, which interact to facilitate the development of the ovary and functional oocytes. Revealing the ovarian transcriptome profile of highly fecundate fishes provides insights into oocyte production in teleosts. Here we use a comprehensive tissue-specific RNA sequencing which yielded 102.3 billion clean bases to analyze the transcriptional profiles of the ovary compared with other organs (liver, kidney, ovary, testis, fin, cauda and gill) and juvenile tissues of European anchovy. We conducted a comparative transcriptome and positive selection analysis of seven teleost species with varying fecundity rates to identify genes potentially involved in oogenesis and oocyte development. Of the 2,272 single copies of orthologous genes found, up to 535 genes were under positive selection in European anchovy and these genes are associated with a wide spectrum of cellular and molecular functions, with enrichments such as RNA methylation and modification, ribosome biogenesis, DNA repair, cell cycle processing and peptide/amide biosynthesis. Of the 535 positively selected genes, 55 were upregulated, and 45 were downregulated in the ovary, most of which were related to RNA and DNA transferase, developmental transcription factors, protein kinases and replication factors. Overall, our analysis of the transcriptome level in the ovarian tissue of a teleost will provide further insights into molecular processes and deepen our genetic understanding of egg production in highly fecund fish.

Transgenerational bone toxicity in F3 medaka (Oryzias latipes) induced by ancestral benzo[a]pyrene exposure: Cellular and transcriptomic insights

Benzo[a]pyrene (BaP), a ubiquitous pollutant, raises environmental health concerns due to induction of bone toxicity in the unexposed offspring. Exposure of F0 ancestor medaka (Oryzias latipes) to 1 µg/L BaP for 21 days causes reduced vertebral bone thickness in the unexposed F3 male offspring. To reveal the inherited modifications, osteoblast (OB) abundance and molecular signaling pathways of transgenerational BaP-induced bone thinning were assessed. Histomorphometric analysis showed a reduction in OB abundance. Analyses of the miRNA and mRNA transcriptomes revealed the dysregulation of Wnt signaling (frzb/ola-miR-1-3p, sfrp5/ola-miR-96-5p/miR-455-5p) and bone morphogenetic protein (Bmp) signaling (bmp3/ola-miR-96-5p/miR-181b-5p/miR-199a-5p/miR-205-5p/miR-455-5p). Both pathways are major indicators of impaired bone formation, while the altered Rank signaling in osteoclasts (c-fos/miR-205-5p) suggests a potentially augmented bone resorption. Interestingly, a typical BaP-responsive pathway, the Nrf2-mediated oxidative stress response (gst/ola-miR-181b-5p/miR-199a-5p/miR-205), was also affected. Moreover, mRNA levels of epigenetic modification enzymes (e.g., hdac6, hdac7, kdm5b) were found dysregulated. The findings indicated that epigenetic factors (e.g., miRNAs, histone modifications) may directly regulate the expression of genes associated with transgenerational BaP bone toxicity and warrants further studies. The identified candidate genes and miRNAs may serve as potential biomarkers for BaP-induced bone disease and as indicators of historic exposures in wild fish for conservation purposes.

Transcriptome sequencing reveals ovarian immune response and development during female sperm storage in viviparous black rockfish (Sebastes schlegelii)

Black rockfish (Sebastes schlegelii) is a viviparous teleost fish whose spermatozoa were transferred into the female ovary cavity and stored for up to five months and then fertilized with the matured eggs. There is no clarity about the molecular characteristics of ovarian follicles during female sperm storage in Sebastes schlegelii. In this study, histological observation, transcriptomic analysis and hormone level detection were performed in ovaries at stages of pre-mating (PRM), post-mating (POM) and pre-fertilization (PRF). Histological observation displayed that oocytes developed from the primary growth (PG) stage to the mature stage during the three stages. Furthermore, somatic cells around the oocyte were proliferated and spermatozoa were found near the layer of epithelial cells. Transcriptomic analysis showed that there were 437 and 747 differentially expressed genes (DEGs) in ovarian comparison of PRM-vs-POM and POM-vs-PRF, respectively. GO enrichment and KEGG analysis revealed that lots of DEGs from PRM-vs-POM were linked to immune-related pathways, such as antigen processing and presentation, immune response, and complement and coagulation cascade. Meanwhile, seven DEGs associated with immune response were differentially expressed after spermatozoa treatment in ovarian tissue in vitro. While the DEGs from POM-vs-PRF were mostly enriched in the pathways related to homeostasis maintenance and cellular junction and metabolism. In addition, we found increased estrogen (E2) and 11-ketotestosterone (11-KT) level and decreased testosterone level in ovarian follicles during the sperm storage period by ELISA, suggesting that sex hormones are involved in the dynamic change of ovarian follicles. In total, this study could provide new hints for understanding the immune adaption and developmental signatures of ovarian follicles post copulation in black rockfish and other viviparous fish.

Denovo RNA-Seq analysis of ovary and testis reveals potential differentially expressed transcripts associated with gonadal unsynchronization development in Onychostoma macrolepis

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered wild species. Their endangered extinction might be due to their low fertility. To further illustrate the molecular mechanism of gonad development of the male and female O. macrolepis, the present study carried out de novo testicular and ovarian transcriptome sequencing. By comparing ovary and testis, 30,869 differentially expressed unigenes (9870 in female, 20999 in male) were identified. In addition, KEGG and GO analysis suggested that the 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. RT-qPCR analysis results were consistent with transcriptome sequencing that all of gender differentiation-related genes (FOXL2, GDF9, WNT4, CYP19A1, SOX9 and GATA4), temperature-enriched genes (NOVA1, CTGF and NR4A1), clock-related genes (PER2, PER3, CRY1, CRY2, BMAL1 and CIPC) were significantly differential expression in testis compared with ovaries. Taken together, these results revealed a potential molecular mechanism that low fertility of the O. macrolepis might strong correlate with the gonadal dyssynchrony development of the male and female, which might provide theoretical basis and technical support for artificial reproduction and germplasm resource protection of the O. macrolepis.

Sex steroid profile during oocyte development and maturation in the intertidal worm Marphysa madrasi (Polychaeta: Eunicidae) from the east coast of India

Marphysa madrasi is a commercially valuable maturation diet in crustacean aquaculture. This study presents the first detailed investigation of oogenesis in the intertidal polychaete worm M. madrasi and reports the steroid profile during oocyte growth and development. Oogenesis is extraovarian type I, originating from coelomic epithelial cells, with four stages of development - primary growth, early vitellogenic, late vitellogenic, and maturation. The primary growth phase contains oogonial cells and previtellogenic oocyte clusters in the early, mid, and late stages of development form a dispersed ovary attached to blood vessels. The late previtellogenic oocytes detach from the ovary at the onset of vitellogenesis. The detached oocytes complete vitellogenesis and final maturation in the coelomic fluid as solitary free-floating cells without any connection with follicle cells. The worms display asynchronous reproduction with a heterogeneous population of developing oocytes. Steroid extracts from the polychaete homogenates in different stages of oogenesis were identified by HPLC and confirmed by LC-MS/MS. In M. madrasi, two vertebrate-type steroids, pregnenolone (P5) and 17α-hydroxyprogesterone (17-OHP) were detected and quantified. The P5 levels were low in immature worms but increased significantly by ∼ 8.3-fold in the previtellogenic stage and peaked during oocyte maturation. 17-OHP levels were low in immature worms but gradually increase as the oogenesis progress to the primary growth and early vitellogenic phase, with a significant increase (p < 0.001) during the late vitellogenic phase. Although an increase in the concentration of P5 and 17-OHP during vitellogenesis and maturation of oocytes points to a possible role in reproduction, the absence of other vertebrate-type steroids in the investigated polychaete signifies a plausible uptake of P5 and 17-OHP from the environment.

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.

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

The Role of Transposable Elements in Sexual Development

Up to 50% of most mammalian genomes are made up of transposable elements (TEs) that have the potential to mobilize around the genome. Despite this prevalence, research on TEs is only beginning to gain traction within the field of neuroscience. While TEs have long been regarded as “junk” or parasitic DNA, it has become evident that they are adaptive DNA and RNA regulatory elements. In addition to their vital role in normal development, TEs can also interact with steroid receptors, which are key elements to sexual development. In this review, we provide an overview of the involvement of TEs in processes related to sexual development- from TE activity in the germline to TE accumulation in sex chromosomes. Moreover, we highlight sex differences in TE activity and their regulation of genes related to sexual development. Finally, we speculate on the epigenetic mechanisms that may govern TEs’ role in sexual development. In this context, we emphasize the need to further the understanding of sexual development through the lens of TEs including in a variety of organs at different developmental stages, their molecular networks, and evolution.

Function of Foxl2 and Dmrt1 proteins during gonadal differentiation in the olive flounder Paralichthys olivaceus

Study on fish sex differentiation is important both from academic and practical aspects. Foxl2 and Dmrt1 are important transcription factors that should be involved in fish gonadal differentiation, but there is still no direct evidence to clarify their protein functions. Olive flounder Paralichthys olivaceus, an important mariculture fish in China, Japan, and Korea, shows sex-dimorphic growth. In this study, the Foxl2 and Dmrt1 proteins were detected in granulosa cells of the ovary and Sertoli cells of the testis, respectively, showing significant sex-dimorphic expression patterns. Then, bioactive high-purity Foxl2 and Dmrt1 recombinant proteins were obtained in vitro. Furthermore, effects of the recombinant Foxl2 and Dmrt1 during gonadal differentiation period were evaluated by intraperitoneal injection in juvenile fish. Compared with the control group, the male rate in the Dmrt1 group increased from 0 % to 82 %, showing for the first time in fish that the recombinant Dmrt1 could alter the sex phenotype. In addition, transcription levels of cyp19a and its transcription factors also changed after the recombinant Foxl2 and Dmrt1 injection. These findings reveal that Foxl2 and Dmrt1 are vital regulators for fish gonadal differentiation by regulating cyp19a expression, and also provide a new approach for sex control in fish aquaculture.

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

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

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.

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

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

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.

Identification and Characterization of lncRNAs Related to the Muscle Growth and Development of Japanese Flounder (Paralichthys olivaceus)

Long noncoding RNAs (lncRNAs) play an important role in many life activities, but the expression pattern and function of lncRNAs in Japanese flounder skeletal muscle are unclear. In this study, 751 lncRNAs were selected from skeletal muscle in different development stages of the Japanese flounder [stage A: larval 7 days post hatching (dph); stage B: juvenile about 90 dph; stage C (female) and stage D (male): adult about 24 months] using coding potential analysis methods. In total, 232, 211, 194, 28, 29, and 14 differentially expressed lncRNAs and 9549, 8673, 9181, 1821, 1080, and 557 differentially expressed mRNAs were identified in comparisons of A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. We identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and lncRNA-gene interaction networks were constructed using the Cytoscape program. In total, there were 200, 200, 200, 93, 47, and 11 cis-regulation relationships, and 29, 19, 24, 38, 8, and 47 trans-regulation relationships in the comparisons between A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. These results indicate that lncRNA may participate in the development of Japanese flounder skeletal muscle through cis- or trans-acting mechanisms, thus providing a scientific basis for further study of the biological function of lncRNA in Japanese flounder skeletal muscle. Based on these relationships, functional annotation of the related lncRNAs was performed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Differentially expressed genes associated with muscle development were enriched in multiple pairs of comparisons (e.g., differentially expressed genes LOC109640370, LOC109634180, LOC109643555, rusc1, and LOC109626999 were enriched in the actin-binding term, and differentially expressed genes LOC109640370, was, LOC109644970, LOC109643555, and itga9 were enriched in the regulation of the actin cytoskeleton pathway in the KEGG pathway analysis in the comparison of stages C and D). We predicted lncRNA-mRNA, miRNA-mRNA, and lncRNA-miRNA regulatory relationships and constructed interactive networks using Cytoscape software. Co-expression networks show that most lncRNAs interact with one or two predicted miRNAs involved in muscle growth and development. These results provide a basis for further study of the function of lncRNAs on skeletal muscle in different developmental stages of Japanese flounder.

Characteristics and sex dimorphism of 17β-hydroxysteroid dehydrogenase family genes in the olive flounder Paralichthys olivaceus

Sex steroid hormones play important roles in fish sex differentiation, gonadal development and secondary sexual characteristics. Olive flounder Paralichthys olivaceus is a valuable commercial marine fish species and has marked sexual dimorphism. However, the mechanisms of action of sex hormones in flounder sex are still unclear. In this study, a total of ten Hsd17b family genes, including Hsd17b3, -4, -7, -8, -9, -10, -12a, -12b, -14 and -15, were identified in the flounder, which encoded critical enzymes acting on sex steroid synthesis and metabolism. Hsd17b genes were distributed on eight chromosomes. Hsd17b12a and -12b were located on chromosomes 19 and 7, respectively. It was speculated that these two genes were just highly similar rather than different transcripts derived from the same gene. According to the results of domain and motif analyses, they all belonged to the SDR superfamily and contained conserved Hsd17b motifs TGxxxGxG, PGxxxT, NNAG and YxxxK. Analysis of amino acid sequences predicted that Hsd17b1, -4, -7, -12a and -14 were hydrophilic proteins. The stability of Hsd17b1, -3 and -12b proteins was predicted to be low. The various Hsd17b family genes differed in tissue expression pattern, and Hsd17b10, -12a and -12b were highly expressed in the flounder ovary. Moreover, throughout gonadal development, Hsd17b3 was highly expressed in the testis, and Hsd17b1, -12a and -12b were highly expressed in the ovary, suggesting that they might play an important role in testosterone synthesis in the testis or estrogen synthesis in the ovary. Activities of Hsd17b3 at stages I-V were all significantly higher in the testis than in the ovary (P < 0.05, P < 0.01). Transfection analysis in HEK293T cells showed that Hsd17b1 and -3 were located in both the cytoplasm and nucleus. Additionally, after challenging fish with tamoxifen, Hsd17b3 expression level in the testis decreased significantly (P < 0.01), and in the ovary no significant change was observed. Moreover, the expression of Hsd17b1 in the ovary was significantly upregulated after injection with flutamide (P < 0.05). These findings introduce the characteristics of the flounder Hsd17b in subfamily, which contribute to our understanding of the regulation of sex steroid hormone synthesis in fish gonadal development.

Gonadal transcriptome analysis of the common carp, Cyprinus carpio: Identification of differentially expressed genes and SSRs

Common carp (Cyprinus carpio) is a world-wide freshwater fish of eutrophic waters. C. carpio, have various reproductive traits, including early sexual maturity, that may make them excellent, large, realistic, aquaculture model species. In the present work, de novo assembly of gonadal (testicular and ovarian) transcriptomes from juvenile common carp was performed to identify genes involved in gonadal development. A total of 81,757 and 43,257 transcripts with average lengths of 769 and 856 bp, were obtained from the immature testicular and ovarian transcriptomes, respectively. About 84,367 unigenes were constructed after removing redundancy involving representation of transcripts in both gonadal transcriptomes. Gene ontology (39,171 unigenes), clusters of orthologous group's analysis (6651 unigenes) and Kyoto encyclopedia of genes, and genomes automatic annotation server analysis (4783 unigenes) were performed to identify potential genes along with their functions. Furthermore, 18,342 (testis) and 8693 (ovary) simple sequence repeats were identified. About 298 differentially expressed genes were identified, of which 171 and 127 genes were up-regulated in testis and ovary, respectively. Quantitative real-time reverse transcription PCR was performed to validate differential expression of selected genes in testis and ovary. Nearly 809 genes related to reproduction were identified, sex-wise expression pattern of genes related to steroid synthesis, endocrine regulation, germ cell maintenance and others factors related to gonadal differentiation was observed, and expression analysis of nanos, ad4bp/sf-1, and gdf9 was performed. The present study identified certain important genes/factors involved in the gonadal development of C. carpio which may provide insights into the understanding of sex-differentiation and gonadal development processes.

Rapid evolution of piRNA pathway and its transposon targets in Japanese flounder (Paralichthys olivaceus)

Piwi-interacting RNA (piRNA) pathway is essential for germline specification, gametogenesis, and genome integrity as defense against transposable elements (TEs). This pathway has been suggested to have undergone rapid adaptive evolution in spite of its conserved role in TE silencing. However, with diverse sexual development patterns, piRNA pathway evolution and its adaptation to transposon activity in teleost lineages remain less known. This study illustrated the evolutionary significance of piRNA pathway via a systematic comparative analysis on diverse teleosts, including flatfish lineages. Molecular evolution of piRNA pathway and microRNA/small interfering RNA pathway genes indicated a faster evolution of piRNA pathway in teleosts than in mammals. Positive selection was detected at the PAZ (Piwi-Argonaute-Zwille) domain involved in Piwi-piRNA interaction, thereby suggesting that the amino acid substitutions were adaptive to their functions in teleost piRNA pathway. Notably, Piwil1 evolved faster in Japanese flounder than in other teleosts, and the piRNA pathway genes expressed higher in testes than in ovaries. In addition, gonadal transcriptomic analysis revealed male under-represented TE families mainly from DNA transposons, which were the potential targets of the complex formed by male-biased Piwi genes and male over-represented piRNAs in Japanese flounder testes. The potential piRNA-TE regulatory relationships suggested that the rapidly evolved piRNA pathway in Japanese flounder was likely involved in the regulation of transposon activity in germlines and could play important roles in Japanese flounder gonadal development and spermatogenesis.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

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.

Gene expression profiles alteration after infection of virus, bacteria, and parasite in the Olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is one of economically valuable fish species in the East Asia. In comparison with its economic importance, available genomic information of the olive flounder is very limited. The mass mortality caused by variety of pathogens (virus, bacteria and parasites) is main problem in aquaculture industry, including in olive flounder culture. In this study, we carried out transcriptome analysis using the olive flounder gill tissues after infection of three types of pathogens (Virus; Viral hemorrhagic septicemia virus, Bacteria; Streptococcus parauberis, and Parasite; Miamiensis avidus), respectively. As a result, we identified total 12,415 differentially expressed genes (DEG) from viral infection, 1,754 from bacterial infection, and 795 from parasite infection, respectively. To investigate the effects of pathogenic infection on immune response, we analyzed Gene ontology (GO) enrichment analysis with DEGs and sorted immune-related GO terms per three pathogen groups. Especially, we verified various GO terms, and genes in these terms showed down-regulated expression pattern. In addition, we identified 67 common genes (10 up-regulated and 57 down-regulated) present in three pathogen infection groups. Our goals are to provide plenty of genomic knowledge about olive flounder transcripts for further research and report genes, which were changed in their expression after specific pathogen infection.

GATA4 is a transcriptional regulator of R-spondin1 in Japanese flounder (Paralichthys olivaceus)

GATA4 is a well-known transcription factor of the GATA family implicated in regulation of sex determination and gonadal development in mammals. In this study, we cloned the full-length cDNA of Paralichthys olivaceus gata4 (Po-gata4). Phylogenetic, gene structure, and synteny analysis showed that Po-GATA4 is homologous to GATA4 of teleost and tetrapod. Po-gata4 transcripts were detected in Sertoli cells, spermatogonia, oogonia and oocytes, with higher transcript levels overall in the testis than the ovary. The promoter region of P. olivaceus R-spondin1was found to contain a GATA4-binding motif. Results of CBA (cleaved amplified polymorphic sequence-based binding assay) indicated that GATA4 could indeed bind to the promoter sequence of R-spondin1. Moreover, human GATA4 recombinant protein could upregulate R-spondin1 in P. olivaceus ovary cells and FBCs (flounder brain cell line). In FBCs, overexpression of Po-gata4 resulted in elevated transcript levels of R-spondin1. Taken together, our results indicate that Po-GATA4 is involved in gonadal development by regulating R-spondin1 expression.

Different expression patterns of sperm motility-related genes in testis of diploid and tetraploid cyprinid fish†

Sperm motility is an important standard to measure the fertility of male. In our previous study, we found that the diploid spermatozoa from allotetraploid hybrid (4nAT) had longer durations of rapid and slow progressive motility than haploid spermatozoa from common carp (COC). In this study, to explore sperm motility-related molecular mechanisms, we compared the testis tissues transcriptomes from 2-year-old male COC and 4nAT. The RNA-seq data revealed that 2985 genes were differentially expressed between COC and 4nAT, including 2216 upregulated and 769 downregulated genes in 4nAT. Some differentially expressed genes, such as tubulin genes, dynein, axonemal, heavy chain(dnah) genes, mitogen-activated protein kinase(mapk) genes, tektin 4, FOX transcription factors, proteasome genes, and ubiquitin carboxyl-terminal hydrolase(uchl) genes, are involved in the regulation of cell division, flagellar and ciliary motility, gene transcription, cytoskeleton, energy metabolism, and the ubiquitin–proteasome system, suggesting that these genes were related to sperm motility of the 4nAT. We confirmed the differential expression of 12 such genes in 4nAT by quantitative PCR. By western blotting, we also confirmed increased expression of Uchl3 in 4nAT testis. In addition, we identified 1915 and 2551 predicted long noncoding RNA (lncRNA) transcripts from testis tissue transcriptomes of COC and 4nAT, respectively. Of these, 1575 lncRNAs were specifically expressed in 4nAT and 939 were specifically expressed in COC. This study provides insights into the transcriptome profile of testis tissues from diploid and tetraploid, which are useful for research on regulatory mechanisms behind sperm motility in male polyploidy.

Transcriptome profiling provides gene resources for understanding gill immune responses in Japanese flounder (Paralichthys olivaceus) challenged with Edwardsiella tarda

Marine organisms are commonly under threats from various pathogens. Edwardsiella tarda is one of the fish pathogens that seriously infect cultured and wild fish species. Bacteremia caused by E. tarda can be a fatal disease in humans. Fish gill is a mucosa-associated lymphoid tissue that directly contacted with sea water. Generating gill transcriptomic resources that challenged by E. tarda is crucial for understanding the molecular mechanisms underlying gill immune responses. In this study, we performed transcriptome profiling of gene expression in Japanese flounder gills (Paralichthys olivaceus) challenged by E. tarda with different stress duration. An average of 40 million clean reads per library were obtained, of which approximately 83.2% were successfully mapped to the reference genome. 456 and 1037 differential expressed genes (DEGs) were identified at 8 h and 48 h post-injection, respectively. Gene annotation analysis and protein-protein interaction networks were conducted to obtain the key interaction relationships of immune-related DEGs during pathogens infection. 24 hub genes with multiple protein-protein interaction relationships or involved in multiple KEGG signaling pathways were discovered and validated by qRT-PCR. These hub genes mainly participated in Leukocyte transendothelial migration signaling pathway, B cell receptor signaling pathway, Wnt signaling pathway and Apoptosis signaling pathway. This study represents the first gill transcriptomic analysis based on protein-protein interaction networks in fish and provides valuable gene resources for understanding the fish gill immunity, which can pave the way to understand the molecular mechanisms of immune responses with E. tarda infection.

Transcriptome profiling based on protein-protein interaction networks provides a core set of genes for understanding blood immune response mechanisms against Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

Marine organisms are commonly under threat from various pathogens. Edwardsiella tarda is one of the fish pathogens that can infect both cultured and wild fish species. E. tarda can also infect other vertebrates, including amphibians, reptiles, and mammals. Bacteremia caused by E. tarda can be a fatal disease in humans. Blood acts as a pipeline for the fish immune system. Generating blood transcriptomic resources from fish challenged by E. tarda is crucial for understanding molecular mechanisms underlying blood immune response process. In this study, we performed transcriptome-wide gene expression profiling of Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h E. tarda stress. An average of 37 million clean reads per library was obtained, and approximately 85.6% of these reads were successfully mapped to the reference genome. In addition, 808 and 1265 differential expression genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to search immune-related DEGs. A protein-protein interaction network was constructed to obtain the interaction relationship of immune genes during pathogens stress. Based on KEGG and protein association networks analysis, 30 hub genes were discovered and validated by quantitative RT-PCR. This study represents the first transcriptome analysis based on protein-protein interaction networks in fish and provides us with valuable gene resources for the research of fish blood immunity, which can significantly assist us to further understand the molecular mechanisms of humans and other vertebrates against E. tarda.

Sequencing-based network analysis provides a core set of gene resource for understanding kidney immune response against Edwardsiella tarda infection in Japanese flounder

Marine organisms are under a frequent threat from various pathogens. Edwardsiella tarda is one of the major fish pathogens infecting both cultured and wild fish species. It can also infect a variety of other vertebrates, including amphibians, reptiles, and mammals, and bacteremia caused by E. tarda can be fatal in humans. The kidney is the largest lymphoid organ in fish, and generating kidney transcriptomic information under different stresses is crucial for understanding molecular mechanisms underlying the immune responses in the kidneys. In this study, we performed transcriptome-wide gene expression profiling of the Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h of E. tarda infection. An average of 40 million clean reads per library was obtained, and approximately 81.6% of these reads were successfully mapped to the reference genome. In addition, 1319 and 4439 differentially expressed genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to search immune-related DEGs. A protein-protein interaction network was constructed to ascertain the relationship between interacting immune genes during pathogen-induced stress. Based on the KEGG and protein association networks analysis, 24 hub genes were discovered and validated by qRT-PCR. To our knowledge, this study is the first to represent the kidney transcriptome analysis based on protein-protein interaction networks in fish. Our results provide valuable gene resources for further research on kidney immune response in fish, which can significantly improve our understanding of the molecular mechanisms underlying the immune response to E. tarda in humans and other vertebrates.

Sexually dimorphic expression in developing and adult gonads shows an important role of gonadal soma-derived factor during sex differentiation in olive flounder (Paralichthys olivaceus)

Gonadal soma-derived factor (gsdf) is a new member of transforming growth factor beta (TGF-β) superfamily. As a teleost- and gonad-specific growth factor, gsdf has been indicated to play an important role in early germ cell development. However, little is known about its significance in germ cell development of olive flounder (Paralichthys olivaceus). In the present study, a 1338 bp gsdf gene was isolated from P. olivaceus for the first time. Bioinformatic analysis revealed that the genomic structure and synteny relationship of gsdf in teleosts were conserved. Quantitative real-time PCR (qRT-PCR) showed that gsdf expressed before sex gonadal differentiation, and the expression level increased rapidly after initiation of sex differentiation in males. In adult individuals, the expression of gsdf was higher in testis than that in ovary (P<0.01). In situ hybridization (ISH) indicated that gsdf mRNA was detected in the somatic cells of both males and females, and also in the cytoplasm of oocytes. These results suggested that gsdf might play an important role as initial switches to promote testis differentiation and participate in early germ cell development, such as proliferation and differentiation of spermatogonia and oogonia in P. olivaceus.