Integrated analysis of miRNA and mRNA expression profiles identifies potential regulatory interactions during sexual development of Pacific oyster Crassostrea gigas.

One of the most important physiological process in the life cycle of most bivalves is the sexual reproduction required for species propagation. Understanding the molecular regulatory mechanisms of gender determination and gonad development have potential application in bivalve aquaculture, but there is still limited knowledge about this issue. The Pacific oyster Crassostrea gigas is a commercially important bivalve in aquaculture worldwide, as well as a representative specie to study the sexual development in bivalves for its complicated sexuality, including dioecy, sex change and rare hermaphroditism. In this study, an integrated analysis of miRNA and mRNA expression profiles of oyster gonads at different gametogenesis stages was performed to explore their possibly regulatory patterns in sexual development of C. gigas. Based on small RNA sequencing, a total of 122 miRNAs were identified, including 53 conserved and 69 novel miRNAs. Six miRNAs exhibited sex-specific pattern that expressed exclusively from either male or female gonads with increasing abundance throughout the gonadal development. Moreover, 287 putative targets of the six sex-biased miRNAs were predicted by combining previous transcriptome data. GO and KEGG enrichment analysis of target genes revealed that these sex-biased miRNAs might participate in the regulation of gametogenesis. Further, potential miRNA-mRNA interactions based on six sex-biased miRNAs were collected. The homolog of a key participator in testis development, meiotic recombination protein SPO11 in C. gigas (CgSPO11), was predicted as a target of female-biased miRNA lgi-miR-96a. This negative miRNA-mRNA expression relationship was validated in additional samples by quantitative real-time PCR (qRT-PCR) and the regulation of CgSPO11 by lgi-miR-96a was preliminarily confirmed by a dual-luciferase reporter assay. This study provides information regarding the gonad miRNAs of C. gigas and leads to the discovery of miRNA-mRNA interactions in gonad development, which will facilitate further investigations in molecular mechanisms of bivalve sex differentiation. • A total of 122 miRNAs were identified during the Pacific oyster gonad development. • Six sex-biased miRNAs might participate in the regulation of oyster gametogenesis. • The CgSPO11 gene is a target of lgi-miR-96a. [ABSTRACT FROM AUTHOR]