Bivalve mollusks have fascinatingly diverse modes of reproduction. However, research investigating sex determination and reproductive regulation in this group of animals is still in its infancy. In this study, transcriptomes of three ovaries and three testes of Yesso scallop were sequenced and analyzed. Transcriptome comparison revealed that 4394 genes were significantly different between ovaries and testes, of which 1973 were ovary-biased (upregulated in the ovaries) and 2421 were testis-biased. Crucial sex-determining genes that were previously reported in vertebrates and putatively present in bivalves, namely FOXL2, DMRT, SOXH, and SOXE, were investigated. The genes all possessed conserved functional domains and were detected in the gonads. Except for PySOXE, the other three genes were significantly differentially expressed between the ovaries and testes. PyFOXL2 was ovary-biased, and PyDMRT and PySOXH were testis-biased, suggesting that these three genes are likely to be key candidates for scallop sex determination/differentiation. Furthermore, GO and KEGG enrichment analyses were conducted for both ovary- and testis-biased genes. Interestingly, both neurotransmitter transporters and GABAergic synapse genes were overrepresented in the ovary-biased genes, suggesting that neurotransmitters, such as GABA and glycine, are likely to participate in scallop ovary development. Our study will assist in better understanding of the molecular mechanisms underlying bivalve sex determination and reproductive regulation.

Background: Reverse transcription quantitative PCR (RT-qPCR) is widely used for gene expression analysis in various organisms. Its accuracy largely relies on the stability of reference genes, making reference gene selection a vital step in RT-qPCR experiments. However, previous studies in mollusks only focused on the reference genes widely used in vertebrates. Results: In this study, we conducted the transcriptome-wide identification of reference genes in the bivalve mollusk Mizuhopecten yessoensis based on 60 transcriptomes covering early development, adult tissues and gonadal development. A total of 964, 1210 and 2097 candidate reference genes were identified, respectively, resulting in a core set of 568 genes. Functional enrichment analysis showed that these genes are significantly overrepresented in Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ribosomes, energy production, etc. Six genes (RS23, EF1A, NDUS4, SELR1, EIF3F, and OLA1) were selected from the candidate genes for RT-qPCR validation, together with 6 commonly used reference genes (ACT, CYTC, HEL, EF1B, GAPDH and RPL16). Stability analyses using geNorm, NormFinder and the comparative delta-Ct method revealed that the new candidate reference genes are more stable than the traditionally used genes, and ACT and CYTC are not recommended under either of the three circumstances. There was a significant correlation between the Ct of RT-qPCR and the log 2 (TPM) of RNA-Seq data (Ct = − 0.94 log 2 (TPM) + 29.67, R 2 = 0.73), making it easy to estimate the Ct values from transcriptome data prior to RT-qPCR experiments. Conclusion: Our study represents the first transcriptome-wide identification of reference genes for early development, adult tissues, and gonadal development in the Yesso scallop and will benefit gene expression studies in other bivalve mollusks.

Neuropeptides play essential roles in regulation of reproduction and growth in marine molluscs. But their function in marine bivalves – a group of animals of commercial importance – is largely unexplored due to the lack of systematic identification of these molecules. In this study, we sequenced and analyzed the transcriptome of nerve ganglia of Yesso scallop Patinopecten yessoensis, from which 63 neuropeptide genes were identified based on BLAST and de novo prediction approaches, and 31 were confirmed by proteomic analysis using the liquid chromatography-tandem mass spectrometry (LC-MS/MS). Fifty genes encode known neuropeptide precursors, of which 20 commonly exist in bilaterians and 30 are protostome specific. Three neuropeptides that have not yet been reported in bivalves were identified, including calcitonin/DH31, lymnokinin and pleurin. Characterization of glycoprotein hormones, insulin-like peptides, allatostatins, RFamides, and some reproduction, cardioactivity or feeding related neuropeptides reveals scallop neuropeptides have conserved molluscan neuropeptide domains, but some (e.g., GPB5, APGWamide and ELH) are characterized with bivalve-specific features. Thirteen potentially novel neuropeptides were identified, including 10 that may also exist in other protostomes, and 3 (GNamide, LRYamide, and Vamide) that may be scallop specific. In addition, we found neuropeptides potentially related to scallop shell growth and eye functioning. This study represents the first comprehensive identification of neuropeptides in scallop, and would contribute to a complete understanding on the roles of various neuropeptides in endocrine regulation in bivalve molluscs.

Sex determination and differentiation have long been a research hotspot in metazoans. However, little is known about when and how sex differentiation occurs in most mollusks. In this study, we conducted a combined morphological and molecular study on sex differentiation in the Yesso scallop Patinopecten yessoensis. Histological examination on gonads from 5- to 13-month-old juveniles revealed that the morphological sex differentiation occurred at 10 months of age. To determine the onset of molecular sex differentiation, molecular markers were screened for early identification of sex. The gonadal expression profiles of eight candidate genes for sex determination or differentiation showed that only two genes displayed sexually dimorphic expression, with FOXL2 being abundant in ovaries and DMRT1L in testes. In situ hybridization revealed that both of them were detected in germ cells and follicle cells. We therefore developed LOG10(DMRT1L/FOXL2) for scallop sex identification and confirmed its feasibility in differentiated individuals. By tracing its changes in 5- to 13-month-old juveniles, molecular sex differentiation time was determined: some scallops differentiate early in September when they are 7 months old, and some do late in December when they are 10 months old. Two kinds of coexpression patterns were found between FOXL2 and DMRT1L: expected antagonism after differentiation and unexpected coordination before differentiation. Our results revealed that scallop sex differentiation co-occurs with the formation of follicles, and molecular sex differentiation is established prior to morphological sex differentiation. Our study will assist in a better understanding of the molecular mechanism underlying bivalve sex differentiation.

No abstract