Identification and Characterization of microRNAs in the Gonads of Litopenaeus vannamei Using High-Throughput Sequencing

Although the expression of miRNAs has been widely applied to investigate on gonads, the role of miRNAs in the gonadal development of white Pacific shrimp (Litopenaeus vannamei) remains unknown. In this study, we performed high-throughput sequencing to identify the sex-related microRNAs (miRNAs) that elucidated the regulatory mechanisms on the gonadal differentiation of L. vannamei. We obtained a total of 29,671,557 and 28,526,942 raw reads from the ovaries and testes library, respectively. We then mapped 26,365,828 (92.73%) of the ovarian clean sequences and 23,694,294 (85.65%) of the testicular clean sequences for a transcriptome reference sequence of L. vannamei. After blasting the miRNA sequences against the miRBase database, we identified 153 significantly differentially expressed miRNAs between the ovaries and testes. To confirm the high-throughput sequencing results, we used a reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) to verify the expression patterns of the seven most differentially expressed miRNAs (i.e., novel_mir23, miR-92b-3p_3, miR-12-5p_2, novel_mir67, miR-279_1, let-7-5p_6, miR-263a-5p_1). According to the results of RT-qPCR, most of the miRNAs were expressed consistently with the high-throughput sequencing results. In addition, the target genes significantly enriched several Kyoto Encyclopedia of Genes and Genome (KEGG) pathways that were closely related to gonadal differentiation and development, including extracellular matrix-receptor interaction, Hedgehog signaling pathway, protein digestion and absorption and cell adhesion molecules (CAMs). This study revealed the first miRNAs sequencing of L. vannamei gonads. We identified sex-related differentially expressed miRNAs and KEGG pathways, which will be helpful to facilitate future research into the regulatory mechanism on the gonadal differentiation of L. vannamei.