Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

BMC Genomics volume 18, Article number: 613 (2017) Cite this article

Abstract

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Background

Teleost fishes exhibit remarkably diverse and plastic sexual development patterns [1], with fish sex determined either by genetics (i.e., genetic sex determination, or GSD) or by environment (i.e., environmental sex determination, or ESD). However, individuals with GSD can also be influenced by many external variables, primarily temperature—that is, GSD + EE (environmental effects) [2]. Despite the various forms of sex determination found in fishes, it is theorized that the genes involved have been well preserved throughout time. Several molecular factors were reported to be responsible for sex determination in gonochoristic species, such as Dmy in medaka [3], Amhr2 in tiger pufferfish [4], Sdy in rainbow trout [5], Amhy in Patagonian pejerrey [6], and Dmrt1 in half-smooth tongue sole [7]. These studies all detail the large diversity in sex determination. Generally speaking, sex determination and differentiation in fishes is a challenge because of the wide variety of species and much has yet to be learned about their molecular mechanisms.

The dark sleeper (Odontobutis potamophila) is a gonochoristic species widely distributed in the river systems of southeastern China [8]. It is an important commercial fish species due to its delicious taste and high nutritional value, therefore in recent years the dark sleeper has become a very promising aquaculture species in China. A previous survey on dark sleepers found that the male individual grows substantially larger and at a quicker rate than the female [9]. Thus, male individuals have more economic value than females. However, there has been little research addressing genetically selective breeding of the dark sleeper due to the dearth of genomic and transcriptomic resources available and since few of its sequences have been deposited in the NCBI GenBank. Therefore, it is important to understand the regulatory mechanisms of sex determination and gonadal differentiation in this species and enlarge the public EST database. In order to understand these mechanisms in detail at the molecular level and to provide a theoretical basis for selective breeding, the regulatory mechanisms associated with gonadal differences need to be studied at the genomic level.

RNA-seq represents an effective means of providing information about gene expression levels, information which has been widely used for analyzing the variety of expressed genes and demonstrating key genes which are expressed at a particular moment. For instance, the use of RNA-seq analysis has helped researchers to understand the molecular mechanisms underlying sex determination and differentiation and to identify differentially expressed genes in rainbow trout (Oncorhynchu smykiss) [10], Nile tilapia(Oreochromis niloticus) [11], yellow catfish (Pelteobagrus fulvidraco) [12, 13], and turbot (Scophthalmus maximus) [14]. These transcriptomic data have provided information on genes in gonads in specific conditions and identified sex-differentially expressed genes.

MicroRNAs (miRNAs) are small noncoding RNA molecules (18–22 nucleotides) which post-transcriptionally regulate gene expression [15, 16]. Most mature miRNA sequences are conserved among fish, amphibians, birds and mammals [17]. Studies have suggested that miRNAs play a crucial role in many biological processes, including metabolism, cell proliferation, differentiation, apoptosis, and developmental timing [18]. Recently, miRNA expression has been examined in the gonads of yellow catfish [19], Nile tilapia [20], and olive flounder [21] by using deep RNA-sequencing technology. These studies suggested that miRNAs are important for sex differentiation and sexual development. In addition, miRNAs are able to target numerous genes, making them a probable regulator in biological processes. Thus, gaining insight into miRNA-mRNA regulatory networks enables a greater understanding of gene expression at the post-transcriptional level [22]. Additional research related to miRNA and mRNA co-expression has also recently been conducted, focusing on the comprehensive network of post-transcriptional regulations [23,24,25].

In this study, an integration analysis was conducted in order to identify target genes of miRNAs and characterize their functional roles in gonad tissues of the dark sleeper. These findings provide us with further insight into the regulatory mechanisms associated with sex determination and gonadal differentiation. We used bioinformatics analysis and miRNA prediction algorithms in order to find miRNA-mRNA pairs. Real-time PCR was used to validate the differences between testes and ovaries for various differentially expressed genes and miRNAs. Through this research, we identified a large number of SNPs that will facilitate genetic research and selective breeding of this species, which have so far been hindered by a lack of genomic data. The findings of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Methods

Sample preparation and RNA extraction

Nine one-year-old male (35 ± 6.34 g) and nine one-year-old female (22 ± 5.65 g) dark sleepers were obtained from Nanjing Fisheries Research Institute, Jiangsu Province, China on February 16, 2016, before the breeding season. These dark sleepers were temporarily cultured in aquaria with aeration for two weeks in order to acclimatize them to laboratory conditions before they were dissected. The temperature and PH of the freshwater in the aquaria were kept at 22–24 °C and 7.2–7.4, respectively. Sequential feed restrictions were utilized for two days before the dark sleepers were used for the experiments. One-third of the culture water was replaced every day. Before excision of the gonads, genders were identified by morphological observation. Then the fish were placed in an ice bath for 3–5 min until they were lightly anesthetized. Testes and ovaries from 9 male and 9 female dark sleepers were collected for RNA extraction. Total RNA was extracted from all harvested organs using a high purity RNA fast extract reagent (BioTeke, Beijing, China). Samplings of testes (OT a, OT b, and OT c) and ovaries (OO a, OO b, and OO c) had three biological replicates, each made up of three different individual gonad tissues. The pieces of gonad from each fish were fixed in 4% PFA and sectioned at 5 μm for hematoxylin and eosin (HE) staining and observation by light microscopy to verify that the gonad sampled were in same development stage. Classification of gonad maturation state in this study refers to the Chinese black sleeper, Bostrichthys sinensis [26]. All experiments were performed according to the Guidelines for the Care and Use of Laboratory Animals in China. This study was also approved by the Ethics Committee of Experimental Animals at Nanjing Normal University. Total RNA was extracted using miRNeasy Kit (QIAGEN, USA) following the manufacturer’s protocol. The total RNA quantity and purity were analyzed with Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number > 7.0.

Transcriptome sequencing and analysis

For six cDNA library constructions, approximately 5 μg of the total RNA per sample was used for the RNA sample preparations. The library for sequencing was generated using an IlluminaTruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA). Transcriptome sequencing was paired-end sequenced with the IlluminaHiSeq 2500platform, which generated approximately 125 bp paired-end (PE) raw reads by LC Sciences (Houston, TX, USA). The clean reads from the six transcriptomes were obtained from raw data by filtering out adaptor sequences and low-quality reads, then the remaining clean reads were assembled using Trinity software for de novo transcriptome assembly without a reference genome. The quality of the assembly was critically assessed by LC Sciences before subsequent analysis. All the non-redundant sequences were annotated against the Swiss-Prot database, non-redundant protein sequences (Nr), protein family (Pfam), KEGG Ortholog database (KEGG), eukaryotic Ortholog Groups (KOG), and GeneOntology (GO). The expression level of each transcript was measured as the number of clean reads mapped to its sequence. Gene expression level was measured by RPKM with RSEM 1.2.3. We determined the FDR threshold by using DESeq. FDR < 0.05 and fold change >2 were used to identify differentially expressed genes.

Small RNA library sequencing and analysis

To identify the small RNAs involved in gonad development and sex differentiation, six small RNA libraries were generated from the male and female gonad samples using the IlluminaTruseq Small RNA Preparation kit (Illumina, San Diego, CA, USA) according to Illumina’s TruSeq Small RNA Sample Preparation Guide. Quantity and integrity were evaluated by Aligent 2100 Bioanalyzer, and then the libraries were sequenced by Illumina Hiseq2500 50SE (single end) at the LC-BIO (Hangzhou, China) following instructions for running the instrument. The raw sequences were processed using the Illumina pipeline program. After disregarding contaminated reads, including adapter dimers, junk, low complexity, common RNA families (rRNA, tRNA, snRNA, snoRNA), and repeats, the clean reads were filtered with the software package ACGT101-miR-v4.2 (LC Sciences, Houston, Texas, USA).

The clean sequence reads were mapped with miRBase 21.0, allowing a mismatch of one or two nucleotide bases, to identify known miRNAs and novel 3p- and 5p–derived miRNAs. The differentially expressed miRNAs (DEmiRNAs) between samples were identified by using a Student’s t-test based on the experimental design. In this test, the significance threshold was set at 0.05.

miRNA target predictions

To predict the genes targeted by DE miRNAs, two computational target prediction algorithms (TargetScan 50 and miRanda 3.3a) were used to identify miRNA binding sites. Finally, the data predicted by TargetScan 50 and miRanda 3.3a were combined and the overlaps were calculated.

Quantitative real-time PCR validation of miRNAs and mRNA

In total, 43 mRNAs and 16 miRNAs were verified by using quantitative real-time PCR (qRT-PCR). For mRNA, primers for qRT-PCR are listed in Additional file 1: Table S1, qRT-PCR method with β-actin as an internal control was used to explore mRNA expression levels. qRT-PCR was performed with an SYBR Green Master kit according to the manufacturer’s protocol (Roche, Basel, Switzerland). The experiments were carried out in triplicate with a total volume of 20 μL in an ABI StepOnePlus, containing 10 μL of SYBR Green Master, 4 μL of cDNA (500 ng), and 3 μL of forward and reverse primers (2 μmol/L). The qRT-PCR was programmed at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s, and 55 °C for 1 min. The expression level was calculated by 2-△△CT method and subjected to statistical analysis. For miRNA quantification, the specific RT primers and stem-loop primers are shown in Additional file 2: Table S2. U6 was used as an internal control. The selected miRNAs were analyzed by quantifying the miRNA stem-loop. Quantitative real-time PCR (qRT-PCR) was performed on ABI StepOnePlus system (Applied Biosystems, Foster, CA, USA) using qRT-PCR reagents provided by Toyobo. The thermal cycling program was set as follows: denaturation at 95 °C for 30 s (s), and then 40 cycles of amplification including denaturation at 95 °C for 5 s, annealing at 61 °C for 30 s, and extension at 72 °C for 30 s; after 40 cycles was final extension at 72 °C for 1 min.

Detection and validation of SNPs

Potential SNPs were called using SAMtools software. SAMtools provides various utilities for manipulating alignments in the SAM format, including sorting, merging, indexing, and generating alignments in a per-position format. First, we called SNPs with the SAMtools mpileup utility. We then piped the BCF output file to SAMtools bcftools, which converted the BCF file into a VCF file. Afterwards, we piped the VCF file into vcfutils.pl with the varFilter-d100 option, which retained SNPs that had read depth higher than 100. In order to validate the accuracy of the SNPs predictions, forty SNPs were randomly selected for SNP validation. Fin tissue samples from 20 juvenile dark sleepers were collected for this SNP validation. Genomic DNA was extracted from pterygiophore tissue samples, using a cell/tissue genomic DNA extraction kit (centrifugal column type, Generay PBiotech, Shanghai). Primers were designed to amplify the flanking sequence of the selected SNPs using Primer Premier 5.0. The amplified PCR products were sequenced by utilizing an ABI3730 sequencer, and the products were analyzed using DNAMAN 8.0.

Results

Analysis of transcriptome sequencing of O. potamophila gonad tissue

Very few studies focusing on regulatory mechanisms of sex determination and gonadal differentiation in dark sleepers exist. Identifying the expression patterns of miRNAs and mRNAs in dark sleeper gonads is the first step in elucidating their gonadal development and differentiation. This experiment used mature gonads for the purpose of discovering gender specific genes and their putative regulated miRNAs in addition to finding sex-associated molecular markers.

In order to identify mRNA expression profiles in dark sleeper testes and ovaries, six cDNA libraries representing the testes (OT a, OT b, OT c) and ovaries (OO a, OO b, OO c) were constructed with total RNA and subjected to Illumina deep sequencing. In total, 52,751,768, 43,554,254, and 39,713,188 reads were obtained from the male libraries OT a, OT b, and OT c, respectively, and 51,538,466, 51,356,850, and 48,215,604 reads were obtained from the female libraries OO a, OO b, and OO c, respectively. After quality filtering, approximately 95.52%, 96.60%, and 95.27% reads from testis libraries and 97.52%, 97.64%, 97.64% reads from ovary libraries remained to be qualified for assembling. An overview of the sequencing and assembly results for these six libraries is shown in Additional file 3: Table S3. Through the Trinity de novo assembly method, we obtained 43,494 non-redundant genes, and 81,051 transcripts were achieved with an N50 of 2721 bp and an average length of 1599 bp (Additional file 4: Table S4). The length distribution of genes and transcripts larger than 200 bp is shown in Additional file 5: Figure S1.

All-unigene sequences were searched against Swiss-Prot, Nr, Pfam, KEGG, KOG, and GO databases, which returned 19,518 (44.88%), 24,291 (55.85%), 19,364 (44.52%), 13,895 (31.95%), 18,624 (42.82%), and 17,358 (39.91%) matches, respectively (Additional file 6: Table S5).

A total of 8103 male-biased and 8437 female-biased genes were identified in the gonad of the dark sleeper (Additional file 7: Figure S2), which showing significant expression difference between male and female. Clustering of the samples based on grouped males and females separately. A majority of the main 100 differentially expressed genes are female-biased genes, which were included in ZP1, ZP2 and ZP3. Among eight male-biased genes, half of them (FABP2, slc7a2, F10 and SUN1) were annotated (Fig. 1).

Fig. 1
figure 1

Heatmap of the differentially expressed genes of gonads. Heatmap of the count data of gonads for the differentially expressed genes between male and female individuals (Only the top 100 genes are included). Please see Additional file 16 for high resolution image

Full size image

Gene ontology (GO) annotation was performed to classify sex-biased genes (Fig. 2). There were 3085 and 5429 DEGs from males and females, respectively, were assigned to a GO category. The most enriched GO-terms for male-biased DEGs included “RNA-dependent DNA replication” in the biological process category, “integral to membrane” in the molecular function category and “zinc ion binding” and “ATP binding” in the cellular component category. As for female-biased DEGs, the most enriched GO-terms at biological process level was “regulation of transcription, DNA-dependent”, the most enriched GO-terms associated with molecular function were “nucleus”, and the most enriched cellular component were “zinc ion binding” and “ATP binding”.

Fig. 2
figure 2

Functional annotation of male-biased genes (a) and female-biased genes (b) based on gene ontology (GO) categorization. GO analysis was performed for three main categories: biological process, molecular function, and cellular component. Please see Additional file 17 for high resolution image

Full size image

1789 male-biased DEGs and 2450 female-biased DEGs were mapped to KEGG pathways, respectively. Enriched pathways (Fig. 3) associated with male-biased DEGs were axon guidance, ECM-receptor interaction and endocytosis etc. And enriched pathways for female-biased DEGs were axon guidance, aminoacyl-tRNA biosynthesis, endocytosis, homologous recombination, adipocytokine signaling pathway, adherens junction etc.

Fig. 3
figure 3

KEGG pathways enriched in the male-biased genes (a) and female-biased genes (b). Please see Additional file 18 for high resolution image

Full size image

Analysis of miRNA sequencing of O.potamophila gonad tissue

In order to identify miRNA differentiation between male and female individuals, we relied on the six small RNA libraries representing the testes (OT a, OT b, OT c) and ovaries (OO a, OO b, OO c) constructed with total RNA and subjected to Illumina small RNA deep sequencing. After discarding junk sequences, 8.8, 8.1, 8.7, 4.0, 4.1, and 5.9 million clean reads were generated in OT a, OT b, OT c, OO a, OO b, and OO c samples, respectively (Additional file 8: Table S6). The length distributions of miRNAs were similar among libraries in that 21 nt RNAs were the most abundant (Additional file 9: Figure S3). In total, we identified 591 conserved miRNAs belonging to 164 miRNA families, and 214 predicted novel miRNAs in the six small RNA libraries (Additional file 10: Table S7 and Additional file 11: Table S8).

Compared with groups of OT and OO, 75 miRNAs were significantly upregulated in testes, while 68 miRNAs were significantly upregulated in ovaries. (P < 0.05, Additional file 12: Table S9) We found many of the same miRNAs that were reported in the previous study, such as miR-143, miR-145, and miR-138, and detected some novel miRNAs as well. In the two libraries of differentially expressed miRNAs, hierarchical clustering results allowed the six samples to be sorted into two distinct groups (Fig. 4).

Fig. 4
figure 4

Heatmap of the differentially expressed miRNAs for gonads. Heatmap of the count data of gonads for the differentially expressed miRNAs between male and female individuals (p<0.05). Please see Additional file 19 for high resolution image

Full size image

Correlation of DE miRNAs and DE mRNAs of O. potamophila gonad tissue

In order to determine miRNA-mRNA functional interactions, DE miRNAs and their predicted target lists were investigated for cognate mRNA targets using ACGT101-CORR 1.1. There were 53,374 miRNA-mRNA pairs in the gonad tissue of the dark sleeper. There were 27,583 negative miRNA-mRNA interactions with the involvement of 143 DE miRNAs and 16,540 DE mRNAs in total. Our study described the negative interaction between 15 genes with annotation and 15 miRNAs, which were verified by using qRT-PCR (Fig. 5).

Fig. 5
figure 5

miRNA-mRNA negative correlation network. miRNAs and mRNAs were validated using real-time PCR. Please see Additional file 20 for high resolution image

Full size image

Confirmation of differential miRNAs and mRNAs by qRT-PCR

We used real-time qRT-PCR to investigate the relative expression levels of the 16 top miRNAs and 43 mRNAs related to sex determination and differentiation. The data obtained through RNA-seq and qRT-PCR was essentially identical (Table 1 and Table 2). Overall, the results of qRT-PCR validated those of RNA-seq and offered support for the reliability of the differentially expressed miRNAs and genes.

Table 1 Relative miRNA expression of top 16 DE miRNAs for comparison of the OT versus OO groups with respect to miRNA-seq and quantitative real-time PCR
Full size table
Table 2 List of important genes related to sex determination/differentiation in the dark sleeper transcriptome
Full size table

SNP detection and validation

For extended application of the RNA-Seq data, structural variations were discovered using the assembled transcriptome. A total of 1130 and 1226 SNPs were identified in males and females, respectively. A/G was the most abundant type in both males and females (Table 3). Information on all SNPs detected in males and females is presented in Additional file 13: Table S10.

Table 3 SNP type in male and female dark sleeper
Full size table

Forty SNPs in twenty dark sleepers were subjected to experimental validation. Among the 40 primer pairs designed for SNP validation, 34 could amplify target sequences. Within these amplified sequences, 19 SNPs were validated, suggesting that approximately 50% of the predicted SNPs were indeed true SNPs. The primers of validated SNPs are listed in Additional file 14: Table S11.

Discussion

As a commercial fish, the dark sleeper shows a sex dimorphism growth pattern, in that the male dark sleeper grows faster than the female. Artificial breeding and selection of dark sleepers have taken place in recent years. To further facilitate increased production of dark sleepers, the control of their sex is regarded as an important technique. Thus, understanding the detailed mechanisms involved in dark sleeper sex determination and development is necessary. This study presents mRNA and miRNA integrated analyses on the gonad tissue of the dark sleeper. A number of sex specific miRNA-mRNA pairs and SNPs were identified. Taking insight into miRNA-mRNA pairs and their corresponding SNPs can further the understanding of the regulatory mechanisms associated with sex determination and gonadal differentiation in the dark sleeper. This, then, constitutes a genetic breeding basis for this important commercial fish species.

In general, all of the RNA transcripts in the cell were included in the transcriptome data. The transcriptome reflects the mRNAs and miRNAs that are actively expressed under a particular condition. We attempted to acquire the miRNA-mRNA pairs through the use of DE miRNAs and DE mRNAs datasets and miRNA-targeting information. In the present study, a total of 53,374 miRNA-mRNA pairs were identified. In most cases, the negative correlation between miRNAs and their target mRNAs is often considered proof of miRNA targeting [27]. In our results, 27,583 miRNA-mRNA interaction sites with negative correlation were predicted with the involvement of 143 DE miRNAs and 16,540 DE mRNAs, which greatly exceeded our expectations. Moreover, this study validated 43 DE genes related to sex determination or differentiation and the top 16 DE miRNAs through the use of real-time PCR analysis.

In this study, compared with groups of OT and OO, 75 miRNAs were significantly upregulated in testes, while 68 miRNAs were significantly upregulated in ovaries. Increasing evidence suggests that miRNAs play an important role as regulators of reproduction [28, 29], and identifying the miRNA target genes is an important step for understanding their roles in gene regulatory networks.

Our results indicate that the miR-138, miR-145, miR-143, and miR-133 showed significant increased expression in testes compared with that in ovaries. miR-138 can regulate the PI3K/AKT signaling pathway by targeting PDK1 [30]. The PI3K-AKT signaling pathway is involved in many fundamental functions including testis development and spermatogenesis [31]. In yellow catfish, there are more molecules of PI3K-AKT pathway expressed in YY than in XY testis [32]. Therefore, miR-138 might be involved in testis development and spermatogenesis. miR-143 and miR-145 were proved to be co-transcribed in multipotent murine cardiac progenitors [33]. Moreover, miR-145 has been shown to play an important role in male and female differentiation, and Sox9, a male factor acting on Sertoli cells, has been tested as a direct target of miR-145 by luciferase reporter assay [34]. However, Sox9 was not differentially expressed in the gonad of the dark sleeper, which may be a result of the samples we collected were in relatively late stage. The results of our histologic section indicated that the testes were in IV stage, leading to the production of mature spermatids (Additional file 15: Figure S4). Previous studies have suggested that Sox9 overexpressed in the immature gonads of Siberian sturgeons and sablefish [35, 36]. Hence, Sox9 might play a role in an earlier stage of gonadal development in the dark sleeper. In addition, miR-143 is critical for the formation of primordial follicles and regulates ovarian development [37], and it is also a dominant miRNA detected in both ovaries and testes of Atlantic halibut [38]and Nile tilapia [20]. The predicted genes of miR-143 and miR-145 in our study include Foxl2, Fgf16, Fgf20, and Fem-1a. Foxl2 encodes a forkhead transcription factor and is involved in ovary differentiation [39], and showed that it is predominantly expressed in female gonads. However, another ovarian marker, Cyp19a1, showed that it is expressed similarly in both male and female gonads. A possible reason for this is that Cyp19a1 has been proved essential for early sex differentiation and ovarian development, and it is regarded as an early marker of ovarian differentiation [40], but Foxl2 is crucial for maintenance and functioning of adult ovaries [41]. The histologic section verified that the ovaries were in IV stage, during which mature oocytes can be produced (Additional file 14: Figure S4). Fgfs (Fibroblast growth factors) have been shown to participate in a wide variety of biological processes, including cell proliferation, migration, growth, differentiation, sex determination, and sex differentiation [42, 43]. For instance, in zebrafish, Fgf16 is required for pectoral fin bud formation [44] and Fgf20 is necessary for triggering regeneration of its fins [45]. In the XX ovaries and phase II oocytes of Nile tilapia, Fgf16 and Fgf20 showed higher expression, which suggests that Fgf16 and Fgf20 are involved in early oocyte development [46]. Moreover, the gene Fem-1 has been shown to play a role in the sex determination signaling pathway in Caenorhabditis elegans [47]. There are three conservative members of the Fem-1 family in vertebrates—namely Fem-1a, Fem-1b, and Fem-1c. Fem-1a and Fem-1c were identified as female-biased genes in our study. The production of connective tissue growth factor can be limited by miR-133 [48]. miR-133 can regulate the oocyte meiosis by down regulating the expression of cyclin B gene by using double-luciferase reporter genes assay [49].

In contrast with the testes, the most abundant expressed miRNAs in the ovaries of the dark sleeper were miR-221, miR-153, miR-222, and miR-139. Both miR-221 and miR-222 can maintain the undifferentiated state of mammalian spermatogonia [50]. ESR2 and Gata2 are among the predicted targets of miR-221 and miR-222. The ESR2 had higher expression in the testes. The expression of estrogen receptors has been shown male-biased in bluehead wrasses [51], Nile tilapia (70dah) [52], and rainbow trout (60-110dpf) [53]. During mouse gonad development, Gata2 was shown to be sex-biased. The expression of Gata2 was found to be much greater in the fetal ovaries from 11.5 days postcoitum (dpc) onwards. However, it was not detected in the fetal testes during the period studied(10.5–15.5dpc) [54]. The role of Gata2 in teleosts needs to be further studied. The role of miR-139 in fishes has yet to be studied, but an increasing number of studies have suggested that miR-139 suppresses the Wnt signaling activity by targeting several key genes such as TCF-4 and Wnt1 in this pathway [55, 56]. Wnt signaling has been known to be involved in mammalian ovary development [57], and in humans, the overexpression of Wnt4 is associated with sex reversal [58]. Rspo plays an essential role in the Wnt signaling pathway in medaka [59] and activated Rspo1 are sufficient to induce ovarian differentiation in XY medaka [60]. In our study, Rspo was shown to have male-biased expression, similar to the express pattern found in bluehead wrasses [51]. Ctnnb1 ( β-catenin) is a central molecule in the Wnt signaling pathway [61]. Ctnnb1 was highly expressed in both the ovaries and testes of the dark sleeper, but was not shown to be sex-biased. Fst is another important gene down streamed to the Wnt signaling pathway [62]. We detected some Fst genes in our study such as Fstl1, Fstl3 and Fstl5, all of which showed male-biased expression. Our results seem to verify the opinion that the female-specific function of the Wnt pathway is not a general pattern in fishes [63]. Clearly, the function of the Wnt pathway in teleost fishes needs to be further studied. miR-17-5p plays an important role in ovarian function in mice [64]. In zebrafish, miR-17a may be involved in the follicle development and maturation of oocytes [65]. Dmrt1 is one of the predicted target gene of miR-17a. Dmrt genes have been reported to have a conserved function related to testis development in vertebrates [1, 66]. We identified Dmrt1 and Dmrt3a as male-biased genes in this study. Our results suggest that Dmrt1 and Dmrt3a may both be important factors for testicular differentiation. It is worth noting that in our study, Dmrt1 and Foxl2 were both predicted to be the potential target genes of miR-17a. miR-17a in teleost fishes also warrants further study. miR-153 has been well studied in the context of human cancers, and it has been shown to be markedly down regulated in the cells that undergo EMT(epithelial-mesenchymal transition) [67]. Niu et al. [68] revealed that TGF-β2 was negatively regulated by miR-153. The signals that come from the TGF-β superfamily of proteins, such as TGF-β1, TGF-β2, and TGF-β3, play a critical role in testis development and spermatogenesis [69]. Hence, miR-153 may have a role in sex determination and differentiation of fish species.

Several other genes related to sex determination and differentiation were also identified in this study. Dnali1 is a type of axonemal dynein, a large motor protein complex which enables the movement of eukaryotic cilia and flagella. Research has shown that male infertility or lateralization defects can occur due to a disruption of axonemal dynein function [70]. Dnali1 was found to be a male-biased gene in the olive flounder as well [71]. Tdrd7 has been shown to function during spermatogenesis [72] and is regarded as the marker of the testis. In this study, we found Tdrd7 to be a male-biased gene showing a similar expression pattern with the Asian seabass [73]. Piwi proteins have been shown to act as functional partners with Tdrds, working together to regulate spermiogenesis [74]. Piwil was reported to be a male-biased gene in the yellow catfish [75], in contrast to the findings of our study. There is evidence that Piwi might be involved in mammalian oocytes and early embryos [76]. However, further investigation is needed to gain a better understanding of the functions of Piwil. Finally, Zp genes, found on the extracellular matrix of oocytes, function as receptors for the binding of sperm in mammals [77, 78]. In our study, Zona pellucida proteins Zp1, Zp2, and Zp3 were identified as female-biased genes.

Conclusion

The present study describes the first integrated analysis of mRNA and miRNA found in gonad tissue of dark sleeper. There were 27,583 negative miRNA-mRNA interactions with the involvement of 143 DE miRNAs and 16,540 DE mRNAs in total. A significant number of common mRNAs and miRNAs involved in sex determination and differentiation were obtained by comparing the mRNA and miRNA expression profiles in the gonad tissue of the dark sleeper. It is a good case for analyzing mRNA and miRNA co-expression in order to uncover the sex differences between males and females of non-model fish species by using NGST.

Abbreviations

Cyp19a :

Cytochrome P450, family 19, subfamily A

Dmrt1 :

Doublesex and mab-3-related transcription factor 1

Dmrt3a :

Doublesex and mab-3-related transcription factor 3a

Dnali1 :

Dynein axonemal light intermediate chain 1

ESR2 :

Estrogen receptor 2

Fem-1 :

Feminization-1

Fgf16 :

Fibroblast growth factor 16

Fgf20 :

Fibroblast growth factor 20

Foxl2 :

Forkhead box protein L2

Fst :

Follistatin-like protein

RPKM:

Reads per kilobase per million

Rspo :

R-spondin

Tdrd7 :

Tudor domain-containing protein 7

TGF-β:

Transforming growth factor beta

Zp :

Zona pellucida proteins

References

  1. Adolfi MC, Carreira AC, Jesus LW, Bogerd J, Funes RM, Schartl M, Sogayar MC, Borella MI. Molecular cloning and expression analysis of dmrt1 and sox9 during gonad development and male reproductive cycle in the lambari fish. Astyanax altiparanae Reprod Biol Endocrinol. 2015;13:2.

    Article  PubMed  Google Scholar 

  2. Ospina-Alvarez N, Piferrer F. Temperature-dependent sex determination in fish revisited: prevalence, a single sex ratio response pattern, and possible effects of climate change. PLoS One. 2008;3(7):e2837.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, Morrey CE, Shibata N, Asakawa S, Shimizu N, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. Nature. 2002;417(6888):559–63.

    Article  CAS  PubMed  Google Scholar 

  4. Kamiya T, Kai W, Tasumi S, Oka A, Matsunaga T, Mizuno N, Fujita M, Suetake H, Suzuki S, Hosoya S, et al. A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu). PLoS Genet. 2012;8(7):e1002798.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yano A, Guyomard R, Nicol B, Jouanno E, Quillet E, Klopp C, Cabau C, Bouchez O, Fostier A, Guiguen Y. An immune-related gene evolved into the master sex-determining gene in rainbow trout. Oncorhynchus mykiss Curr Biol. 2012;22(15):1423–8.

    Article  CAS  PubMed  Google Scholar 

  6. Hattori RS, Murai Y, Oura M, Masuda S, Majhi SK, Sakamoto T, Fernandino JI, Somoza GM, Yokota M, Strussmann CA. A Y-linked anti-Mullerian hormone duplication takes over a critical role in sex determination. Proc Natl Acad Sci U S A. 2012;109(8):2955–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Chen S, Zhang G, Shao C, Huang Q, Liu G, Zhang P, Song W, An N, Chalopin D, Volff JN, et al. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nat Genet. 2014;46(3):253–60.

    Article  CAS  PubMed  Google Scholar 

  8. Hou X, Zhu F, Yin S, Zhang L, Hu Y, Wang Y, Jia Y, Zhang G, Li L. Genetic diversity of Odontobutis Potamophila from different geographic populations inferred from mtDNA control region. Mitochondrial DNA. 2014;25(5):400–6.

    Article  CAS  PubMed  Google Scholar 

  9. Mei J, Gui JF. Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish. Sci China Life Sci. 2015;58(2):124–36.

    Article  CAS  PubMed  Google Scholar 

  10. Salem M, Rexroad CE 3rd, Wang J, Thorgaard GH, Yao J. Characterization of the rainbow trout transcriptome using sanger and 454-pyrosequencing approaches. BMC Genomics. 2010;11:564.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Tao W, Yuan J, Zhou L, Sun L, Sun Y, Yang S, Li M, Zeng S, Huang B, Wang D. Characterization of gonadal transcriptomes from Nile tilapia (Oreochromis niloticus) reveals differentially expressed genes. PLoS One. 2013;8(5):e63604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Chen X, Mei J, Wu J, Jing J, Ma W, Zhang J, Dan C, Wang W, Gui JF. A comprehensive transcriptome provides candidate genes for sex determination/differentiation and SSR/SNP markers in yellow catfish. Mar Biotechnol. 2015;17(2):190–8.

    Article  CAS  PubMed  Google Scholar 

  13. Lu J, Zheng M, Zheng J, Liu J, Liu Y, Peng L, Wang P, Zhang X, Wang Q, Luan P, et al. Transcriptomic analyses reveal novel genes with sexually dimorphic expression in yellow catfish (Pelteobagrus fulvidraco) brain. Mar Biotechnol. 2015;17(5):613–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ribas L, Pardo BG, Fernandez C, Alvarez-Dios JA, Gomez-Tato A, Quiroga MI, Planas JV, Sitja-Bobadilla A, Martinez P, Piferrer F. A combined strategy involving sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus). BMC Genomics. 2013;14:180.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.

    Article  CAS  PubMed  Google Scholar 

  16. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Kloosterman WP, Plasterk RH. The diverse functions of microRNAs in animal development and disease. Dev Cell. 2006;11(4):441–50.

    Article  CAS  PubMed  Google Scholar 

  18. Gong G, Sha Z, Chen S, Li C, Yan H, Chen Y, Wang T. Expression profiling analysis of the microRNA response of Cynoglossus semilaevis to Vibrio anguillarum and other stimuli. Mar Biotechnol. 2015;17(3):338–52.

    Article  CAS  PubMed  Google Scholar 

  19. Jing J, Wu J, Liu W, Xiong S, Ma W, Zhang J, Wang W, Gui JF, Mei J. Sex-biased miRNAs in gonad and their potential roles for testis development in yellow catfish. PLoS One. 2014;9(9):e107946.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Xiao J, Zhong H, Zhou Y, Yu F, Gao Y, Luo Y, Tang Z, Guo Z, Guo E, Gan X, et al. Identification and characterization of microRNAs in ovary and testis of Nile tilapia (Oreochromis niloticus) by using solexa sequencing technology. PLoS One. 2014;9(1):e86821.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Gu Y, Zhang L, Chen X. Differential expression analysis of Paralichthys olivaceus microRNAs in adult ovary and testis by deep sequencing. Gen Comp Endocrinol. 2014;204:181–4.

    Article  CAS  PubMed  Google Scholar 

  22. Liu W, Ling S, Sun W, Liu T, Li Y, Zhong G, Zhao D, Zhang P, Song J, Jin X, et al. Circulating microRNAs correlated with the level of coronary artery calcification in symptomatic patients. Sci Rep. 2015;5:16099.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Fuller-Carter PI, Carter KW, Anderson D, Harvey AR, Giles KM, Rodger J. Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration. BMC Genomics. 2015;16:591.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Tang Z, Yang Y, Wang Z, Zhao S, Mu Y, Li K. Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs. Sci Rep. 2015;5:15544.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Zhuang X, Li Z, Lin H, Gu L, Lin Q, Lu Z, Tzeng CM. Integrated miRNA and mRNA expression profiling to identify mRNA targets of dysregulated miRNAs in non-obstructive azoospermia. Sci Rep. 2015;5:7922.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Zhang YT, Liu DT, Zhu Y, Chen SX, Hong WS. Cloning and olfactory expression of progestin receptors in the Chinese black sleeper Bostrichthys sinensis. Gen Comp Endocrinol. 2016;230–231:87–102.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Alshalalfa M. MicroRNA response elements-mediated miRNA-miRNA interactions in prostate cancer. Adv Bioinforma. 2012;2012:839837.

    Google Scholar 

  28. Murchison EP, Stein P, Xuan Z, Pan H, Zhang MQ, Schultz RM, Hannon GJ. Critical roles for dicer in the female germline. Genes Dev. 2007;21(6):682–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Lei L, Jin S, Gonzalez G, Behringer RR, Woodruff TK. The regulatory role of dicer in folliculogenesis in mice. Mol Cell Endocrinol. 2010;315(1–2):63–73.

    Article  CAS  PubMed  Google Scholar 

  30. Liu Y, Yang K, Sun X, Fang P, Shi H, Xu J, Xie M, Li M. MiR-138 suppresses airway smooth muscle cell proliferation through the PI3K/AKT signaling pathway by targeting PDK1. Exp Lung Res. 2015;41(7):363–9.

    Article  PubMed  Google Scholar 

  31. Ciraolo E, Morello F, Hobbs RM, Wolf F, Marone R, Iezzi M, Lu X, Mengozzi G, Altruda F, Sorba G, et al. Essential role of the p110beta subunit of phosphoinositide 3-OH kinase in male fertility. Mol Biol Cell. 2010;21(5):704–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wu J, Xiong S, Jing J, Chen X, Wang W, Gui JF, Mei J. Comparative Transcriptome analysis of differentially expressed genes and signaling pathways between XY and YY testis in yellow catfish. PLoS One. 2015;10(8):e0134626.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Cordes KR, Sheehy NT, White MP, Berry EC, Morton SU, Muth AN, Lee TH, Miano JM, Ivey KN, Srivastava D. miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature. 2009;460(7256):705–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Yang B, Guo H, Zhang Y, Chen L, Ying D, Dong S. MicroRNA-145 regulates chondrogenic differentiation of mesenchymal stem cells by targeting Sox9. PLoS One. 2011;6(7):e21679.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Berbejillo J, Martinez-Bengochea A, Bedo G, Vizziano-Cantonnet D. Expression of dmrt1 and sox9 during gonadal development in the Siberian sturgeon (Acipenser baerii). Fish Physiol Biochem. 2013;39(1):91–4.

    Article  CAS  PubMed  Google Scholar 

  36. Smith EK, Guzman JM, Luckenbach JA. Molecular cloning, characterization, and sexually dimorphic expression of five major sex differentiation-related genes in a Scorpaeniform fish, sablefish (Anoplopoma fimbria). Comp Biochem Physiol B Biochem Mol Biol. 2013;165(2):125–37.

    Article  CAS  PubMed  Google Scholar 

  37. Zhang J, Ji X, Zhou D, Li Y, Lin J, Liu J, Luo H, Cui S. miR-143 is critical for the formation of primordial follicles in mice. Front Biosci. 2013;18:588–97.

    Article  CAS  Google Scholar 

  38. Bizuayehu TT, Babiak J, Norberg B, Fernandes JM, Johansen SD, Babiak I. Sex-biased miRNA expression in Atlantic halibut (Hippoglossus hippoglossus) brain and gonads. Sex Dev. 2012;6(5):257–66.

    Article  CAS  PubMed  Google Scholar 

  39. Benayoun BA, Dipietromaria A, Bazin C, Veitia RA. FOXL2: at the crossroads of female sex determination and ovarian function. Adv Exp Med Biol. 2009;665:207–26.

    Article  CAS  PubMed  Google Scholar 

  40. Guiguen Y, Fostier A, Piferrer F, Chang CF. Ovarian aromatase and estrogens: a pivotal role for gonadal sex differentiation and sex change in fish. Gen Comp Endocrinol. 2010;165(3):352–66.

    Article  CAS  PubMed  Google Scholar 

  41. Uhlenhaut NH, Jakob S, Anlag K, Eisenberger T, Sekido R, Kress J, Treier AC, Klugmann C, Klasen C, Holter NI, et al. Somatic sex reprogramming of adult ovaries to testes by FOXL2 ablation. Cell. 2009;139(6):1130–42.

    Article  CAS  PubMed  Google Scholar 

  42. Boilly B, Vercoutter-Edouart AS, Hondermarck H, Nurcombe V, Le Bourhis X. FGF signals for cell proliferation and migration through different pathways. Cytokine Growth Factor Rev. 2000;11(4):295–302.

    Article  CAS  PubMed  Google Scholar 

  43. DeVore DL, Horvitz HR, Stern MJ. An FGF Receptor signaling pathway is required for the normal cell migrations of the sex myoblasts in C. elegans hermaphrodites. Cell 1995;83(4):611–620.

  44. Nomura R, Kamei E, Hotta Y, Konishi M, Miyake A, Itoh N. Fgf16 is essential for pectoral fin bud formation in zebrafish. Biochem Biophys Res Commun. 2006;347(1):340–6.

    Article  CAS  PubMed  Google Scholar 

  45. Whitehead GG, Makino S, Lien CL, Keating MT. fgf20 is essential for initiating zebrafish fin regeneration. Science. 2005;310(5756):1957–60.

    Article  CAS  PubMed  Google Scholar 

  46. Sun YL, Zeng S, Ye K, Yang C, Li MH, Huang BF, Sun LN, Zhou LY, Wang DS. Involvement of FGF9/16/20 subfamily in female germ cell development of the Nile tilapia. Oreochromis niloticus Fish Physiol Biochem. 2012;38(5):1427–39.

    Article  CAS  PubMed  Google Scholar 

  47. Doniach T, Hodgkin J. A sex-determining gene, fem-1, required for both male and hermaphrodite development in Caenorhabditis Elegans. Dev Biol. 1984;106(1):223–35.

    Article  CAS  PubMed  Google Scholar 

  48. Duisters RF, Tijsen AJ, Schroen B, Leenders JJ, Lentink V, van der Made I, Herias V, van Leeuwen RE, Schellings MW, Barenbrug P, et al. miR-133 and miR-30 regulate connective tissue growth factor: implications for a role of microRNAs in myocardial matrix remodeling. Circ Res. 2009;104(2):170–8.

    Article  CAS  PubMed  Google Scholar 

  49. Song YN, Shi LL, Liu ZQ, Qiu GF. Global analysis of the ovarian microRNA transcriptome: implication for miR-2 and miR-133 regulation of oocyte meiosis in the Chinese mitten crab, Eriocheir sinensis (Crustacea:Decapoda). BMC Genomics. 2014;15:547.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Yang QE, Racicot KE, Kaucher AV, Oatley MJ, Oatley JM. MicroRNAs 221 and 222 regulate the undifferentiated state in mammalian male germ cells. Development. 2013;140(2):280–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Liu H, Lamm MS, Rutherford K, Black MA, Godwin JR, Gemmell NJ. Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish. Biol Sex Differ. 2015;6:26.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Ijiri S, Kaneko H, Kobayashi T, Wang DS, Sakai F, Paul-Prasanth B, Nakamura M, Nagahama Y. Sexual dimorphic expression of genes in gonads during early differentiation of a teleost fish, the Nile tilapia Oreochromis Niloticus. Biol Reprod. 2008;78(2):333–41.

    Article  CAS  PubMed  Google Scholar 

  53. Baron D, Houlgatte R, Fostier A, Guiguen Y. Expression profiling of candidate genes during ovary-to-testis trans-differentiation in rainbow trout masculinized by androgens. Gen Comp Endocrinol. 2008;156(2):369–78.

    Article  CAS  PubMed  Google Scholar 

  54. Siggers P, Smith L, Greenfield A. Sexually dimorphic expression of Gata-2 during mouse gonad development. Mech Dev. 2002;111(1–2):159–62.

    Article  CAS  PubMed  Google Scholar 

  55. Gu W, Li X, Wang J. miR-139 regulates the proliferation and invasion of hepatocellular carcinoma through the WNT/TCF-4 pathway. Oncol Rep. 2014;31(1):397–404.

    Article  CAS  PubMed  Google Scholar 

  56. Mi L, Li Y, Zhang Q, Zhao C, Peng Y, Yang G, Zheng X. MicroRNA-139-5p regulates C2C12 cell myogenesis through blocking Wnt/beta-catenin signaling pathway. Biochem Cell Biol. 2015;93(1):8–15.

    Article  CAS  PubMed  Google Scholar 

  57. Kopczynski J, Kowalik A, Chlopek M, Wang ZF, Gozdz S, Lasota J, Miettinen M. Oncogenic activation of the Wnt/beta-catenin signaling pathway in signet ring stromal cell tumor of the ovary. Appl Immunohistochem Mol Morphol. 2016;24(5):e28–33.

    Article  CAS  PubMed  Google Scholar 

  58. Jordan BK, Mohammed M, Ching ST, Delot E, Chen XN, Dewing P, Swain A, Rao PN, Elejalde BR, Vilain E. Up-regulation of WNT-4 signaling and dosage-sensitive sex reversal in humans. Am J Hum Genet. 2001;68(5):1102–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Zhou L, Charkraborty T, Yu X, Wu L, Liu G, Mohapatra S, Wang D, Nagahama Y. R-spondins are involved in the ovarian differentiation in a teleost, medaka (Oryzias latipes). BMC Dev Biol. 2012;12:36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Zhou L, Charkraborty T, Zhou Q, Mohapatra S, Nagahama Y, Zhang Y. Rspo1-activated signalling molecules are sufficient to induce ovarian differentiation in XY medaka (Oryzias latipes). Sci Rep. 2016;6:19543.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Chassot AA, Gregoire EP, Lavery R, Taketo MM, de Rooij DG, Adams IR, Chaboissier MC. RSPO1/beta-catenin signaling pathway regulates oogonia differentiation and entry into meiosis in the mouse fetal ovary. PLoS One. 2011;6(10):e25641.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Yao HH, Matzuk MM, Jorgez CJ, Menke DB, Page DC, Swain A, Capel B. Follistatin operates downstream of Wnt4 in mammalian ovary organogenesis. Dev Dyn. 2004;230(2):210–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Nicol B, Guerin A, Fostier A, Guiguen Y. Ovary-predominant wnt4 expression during gonadal differentiation is not conserved in the rainbow trout (Oncorhynchus mykiss). Mol Reprod Dev. 2012;79(1):51–63.

    Article  CAS  PubMed  Google Scholar 

  64. Otsuka M, Zheng M, Hayashi M, Lee JD, Yoshino O, Lin S, Han J. Impaired microRNA processing causes corpus luteum insufficiency and infertility in mice. J Clin Invest. 2008;118(5):1944–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Abramov R, Fu G, Zhang Y, Peng C. Expression and regulation of miR-17a and miR-430b in zebrafish ovarian follicles. Gen Comp Endocrinol. 2013;188:309–15.

    Article  CAS  PubMed  Google Scholar 

  66. Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, Sinclair AH. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009;461(7261):267–71.

    Article  CAS  PubMed  Google Scholar 

  67. Xu Q, Sun Q, Zhang J, Yu J, Chen W, Zhang Z. Downregulation of miR-153 contributes to epithelial-mesenchymal transition and tumor metastasis in human epithelial cancer. Carcinogenesis. 2013;34(3):539–49.

    Article  CAS  PubMed  Google Scholar 

  68. Niu G, Li B, Sun L, An C. MicroRNA-153 inhibits osteosarcoma cells proliferation and invasion by targeting TGF-beta2. PLoS One. 2015;10(3):e0119225.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Itman C, Mendis S, Barakat B, Loveland KL. All in the family: TGF-beta family action in testis development. Reproduction. 2006;132(2):233–46.

    Article  CAS  PubMed  Google Scholar 

  70. Rashid S, Breckle R, Hupe M, Geisler S, Doerwald N, Neesen J. The murine Dnali1 gene encodes a flagellar protein that interacts with the cytoplasmic dynein heavy chain 1. Mol Reprod Dev. 2006;73(6):784–94.

    Article  CAS  PubMed  Google Scholar 

  71. Fan Z, You F, Wang L, Weng S, Wu Z, Hu J, Zou Y, Tan X, Zhang P. Gonadal transcriptome analysis of male and female olive flounder (Paralichthys olivaceus). Biomed Res Int. 2014;291067

  72. Lachke SA, Alkuraya FS, Kneeland SC, Ohn T, Aboukhalil A, Howell GR, Saadi I, Cavallesco R, Yue Y, Tsai AC, et al. Mutations in the RNA granule component TDRD7 cause cataract and glaucoma. Science. 2011;331(6024):1571–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Ravi P, Jiang J, Liew WC, Orban L. Small-scale transcriptomics reveals differences among gonadal stages in Asian seabass (Lates calcarifer). Reprod Biol Endocrinol. 2014;12:5.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Kojima K, Kuramochi-Miyagawa S, Chuma S, Tanaka T, Nakatsuji N, Kimura T, Nakano T. Associations between PIWI proteins and TDRD1/MTR-1 are critical for integrated subcellular localization in murine male germ cells. Genes Cells. 2009;14(10):1155–65.

    Article  CAS  PubMed  Google Scholar 

  75. Sun F, Liu S, Gao X, Jiang Y, Perera D, Wang X, Li C, Sun L, Zhang J, Kaltenboeck L, et al. Male-biased genes in catfish as revealed by RNA-Seq analysis of the testis transcriptome. PLoS One. 2013;8(7):e68452.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Rosenkranz D, Han CT, Roovers EF, Zischler H, Ketting RF. Piwi proteins and piRNAs in mammalian oocytes and early embryos: from sample to sequence. Genom data. 2015;5:309–13.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Hamazaki TS, Nagahama Y, Iuchi I, Yamagami K. A glycoprotein from the liver constitutes the inner layer of the egg envelope (zona pellucida interna) of the fish. Oryzias latipes Dev Biol. 1989;133(1):101–10.

    Article  CAS  PubMed  Google Scholar 

  78. Wassarman PM. Mammalian fertilization: molecular aspects of gamete adhesion, exocytosis, and fusion. Cell. 1999;96(2):175–83.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study was supported by the Jiangsu Province Science and Technology Support Program (BE2013441), Nanjing Science and Technology Plan Projects (201505059), Project Foundation of the Academic Program Development of Jiangsu Higher Education Institution (PAPD), and The Innovation of Graduate Student Training Project of Jiangsu Province (KYLX16-1281)

Availability of data and materials

The dataset generated by deep-sequencing platforms was deposited in the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/Traces/sra). The NCBI Sequence Read Archive accession number is GSE93019.

Author information

Authors and Affiliations

  1. College of Life Sciences, Key Laboratory of Biodiversity and Biotechnology of Jiangsu Province, Nanjing Normal University, Nanjing, Jiangsu, 210023, China

    Cheng Zhao, Guosong Zhang, Shaowu Yin, Zecheng Li & Qintao Wang

  2. Co-Innovation Center for Marine Bio-Industry Technology of Jiangsu Province, Lianyungang, Jiangsu, 222005, China

    Cheng Zhao, Guosong Zhang, Shaowu Yin, Zecheng Li & Qintao Wang

  3. Nanjing Institute of Fisheries Science, Nanjing, Jiangsu, 210036, China

    Shuqiao Chen & Guoqin Zhou

Authors
  1. Cheng Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Guosong Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Shaowu Yin
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zecheng Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qintao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Shuqiao Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Guoqin Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YSW ZC and ZGS conceived this study, and designed and supervised the experiments; ZGS and ZC performed the experiments; ZGS and LZC conducted the data analysis and prepared figures and tables. ZGS, ZC, CSQ and ZGQ wrote the manuscript. WQT and ZC conducted the SNP validation. All of the authors reviewed and approved the manuscript.

Corresponding author

Correspondence to Shaowu Yin.

Ethics declarations

Ethics approval and consent to participate

All experiments were performed according to the Guidelines for the Care and Use of Laboratory Animals in China. This study was also approved by the Ethics Committee for Experimental Animals at Nanjing Normal University (Permit 91, No. SYXK2015–0028).

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict interests of this article.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

RT-qPCR primers for mRNAs. (DOCX 18 kb)

Additional file 2: Table S2.

RT-qPCR primers for miRNAs. (DOCX 13 kb)

Additional file 3: Table S3.

Summary of sequence data generated for dark sleeper transcriptome, and quality filtering. (DOCX 13 kb)

Additional file 4: Table S4.

Assembly statistics of reads. (DOCX 12 kb)

Additional file 5: Figure S1.

Distribution of assembled genes and transcript length. (DOCX 170 kb)

Additional file 6: Table S5.

Blast analysis of non-redundant unigenes against public databases. (DOCX 12 kb)

Additional file 7: Figure S2.

Differentially expressed genes in the two sexes for gonads of dark sleeper. (DOCX 126 kb)

Additional file 8: Table S6.

Overview of reads for sRNA-seq from raw data to high quality reads, and quality filtering. (DOCX 19 kb)

Additional file 9: Figure S3.

Length distribution of small RNA in six miRNA libraries. (DOCX 109 kb)

Additional file 10: Table S7.

List of miRNA member in each family in dark sleeper. (DOCX 23 kb)

Additional file 11: Table S8.

List of know-miRNA and novel-miRNA of dark sleeper. (DOCX 57 kb)

Additional file 12: Table S9.

List of differentially expressed miRNAs of dark sleeper in testis and ovary. (DOCX 34 kb)

Additional file 13: Table S10.

All SNP detected in gonad tissue of dark sleeper. (DOCX 262 kb)

Additional file 14: Table S11.

Information for validated SNPs. (DOCX 17 kb)

Additional file 15: Figure S4.

Histological section analysis of testis and ovary structure in dark sleeper with hematoxylin and eosin staining. (DOCX 1183 kb)

Additional file 16:

High resolution image of Figure 1. (TIF 816 kb)

Additional file 17:

High resolution image of Figure 2. (TIF 1.16 mb)

Additional file 18:

High resolution image of Figure 3. (TIF 1.79 mb)

Additional file 19:

High resolution image of Figure 4. (TIF 1.04 mb)

Additional file 20:

High resolution image of Figure 5. (TIF 371 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, C., Zhang, G., Yin, S. et al. Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila) . BMC Genomics 18, 613 (2017). https://doi.org/10.1186/s12864-017-3995-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-3995-9

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us