Non-coding RNAs (ncRNAs) are involved in several different regulatory pathways including reproduction. In teleost fish, efficacious reproduction is heavily dependent on the completion of the reproductive cycle. The presence of ncRNA, however, and their expression dynamics and putative regulatory role in mature and immature gonads have not yet been extensively explored. Therefore, the abundance of ncRNAs in mature and immature female sharpsnout seabream (Diplodus puntazzo) was investigated. The sharpsnout seabream is a rudimentary hermaphrodite which, in captivity, displays dysfunctions in the gonad maturation process. Our analyses revealed a gonad specific read length distribution with two main peaks representing miRNAs (21-26 nt) and PIWI RNA (27-34 nt). Besides, distinct expression patterns for several ncRNA biotypes including microRNAs (miRNAs), PIWI RNAs (piRNAs), and ribosomal RNAs (rRNAs) were detected. Identified miRNA accounted to 938, corresponding to ~ 13% of obtained transcripts. Among the differential expressed ncRNAs, 10 (~ 7%) were annotated as miRNA, out of which 2 were found in higher abundance in immature gonads (miR-125c and miR-24) and 8 (miR-451, miR-7a, miR-122-1, miR190a, miR129, ENSGACT00000029608, ENSGACT00000029489, and ENSGACT00000029667) were found to be higher expressed in mature gonads. Putative miRNA targets, including long non-coding RNAs (lncRNAs) and genes, are proposed. Target genes are involved in several processes of fish oocyte development, such as steroidogenesis, proteolysis, and apoptosis, and may explain hormone regulation. This study demonstrates a gonad maturation biased ncRNA profile which in turn may support the role of ncRNAs in ovarian physiology and reproductive performance of fish, stressing the specific function of each RNA biotype in oocyte development.
Keywords: Differential expression; Epigenetics; Gonad maturation; Non-coding RNA; Teleost.