Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes

Establishment and Characterization of Bisexually Fertile Triploid Dwarf Surf Clam Mulinia lateralis

Triploids are widely used to rapidly achieve genetic improvements of organisms due to their fast growth and enhanced environmental adaptability. Artificially induced triploids are generally considered to be infertile owing to the obvious inhibition of gonadal development. Recently, some fertile individuals with reduced advantages have been found in triploid bivalves, which is a notable deviation from the original intention of artificially inducing triploids. This study utilized dwarf surf clams (Mulinia lateralis), a promising model organism of bivalves, to develop a model for exploring the potential mechanism of triploid reproduction. The results showed that the optimal induction condition for triploid M. lateralis, determined by orthogonal experiments, was 0.5 mg/L cytochalasin B (CB) to inhibit PB2 for 20 min, resulting in a triploidy rate of 95.57% and a hatching rate of 60.25%. By tracking the development of M. lateralis, we found that the induced triploids could develop normally to maturity and exhibited significant growth and survival advantages post-metamorphosis. Although the triploidy rate exhibited a slight decline overtime, it remained high, with a ratio of 90.63% at 120 dpf. Histological observation confirmed that the gonadal development pattern of triploid M. laterali was similar to that of diploids, but it also showed characteristics such as developmental retardation, few mature gametes, and gamete gigantism. The dynamic expression of genes related to gonadal development provided further molecular evidence for this phenomenon. Additionally, 82.6% of triploid M. laterali exhibited normal spawning behavior, produced fewer but larger viable gametes, and could generate offspring with full developmental potential. Flow cytometry analysis revealed that sperm of triploid M. laterali was aneuploid, with a DNA content of about 1.5 times that of diploid sperm, and the ploidy levels of mating offspring were 2N (DD, diploid female × diploid male), 2.5N (DT, diploid female × triploid male), 2.5N (TD, triploid female × diploid male), and 3N (TT, triploid female × triploid male), respectively. Overall, the artificially induced triploid M. laterali has been confirmed to be bisexually fertile, which will provide a unique model for exploring the underlying mechanisms of advantageous trait formation and fertility regulation in triploids, and offer a valuable platform for the study of ploidy control and polyploidization in bivalves. Please check and confirm that the authors and their respective affiliations have been correctly identified and amend if necessary. Yes, i have checked and it is OK.

Fem-1 Gene of Chinese White Pine Beetle (Dendroctonus armandi): Function and Response to Environmental Treatments

Dendroctonus armandi (Tsai and Li) (Coleoptera: Curculionidae: Scolytinae) is regarded as the most destructive forest pest in the Qinling and Bashan Mountains of China. The sex determination of Dendroctonus armandi plays a significant role in the reproduction of its population. In recent years, the role of the fem-1 gene in sex determination in other insects has been reported. However, the function and expression of the fem-1 gene in Dendroctonus armandi remain uncertain. In this study, three fem-1 genes were cloned and characterized. These were named Dafem-1A, Dafem-1B, and Dafem-1C, respectively. The expression levels of these three Dafem-1 genes vary at different stages of development and between the sexes. In response to different environmental treatments, including temperature, nutrients, terpenoids, and feeding duration, significant differences were observed between the three Dafem-1 genes at different developmental stages and between males and females. Furthermore, injection of double-stranded RNA (dsRNA) targeting the expressions of the Dafem-1A, Dafem-1B, and Dafem-1C genes resulted in increased mortality, deformity, and decreased emergence rates, as well as an imbalance in the sex ratio. Following the interference with Dafem-1A and Dafem-1C, no notable difference was observed in the expression of the Dafem-1B gene. Similarly, after the interference with the Dafem-1B gene, no significant difference was evident in the expression levels of the Dafem-1A and Dafem-1C genes. However, the interference of either the Dafem-1A or Dafem-1C gene results in the downregulation of the other gene. The aforementioned results demonstrate that the Dafem-1A, Dafem-1B, and Dafem-1C genes play a pivotal role in the regulation of life development and sex determination. Furthermore, it can be concluded that external factors such as temperature, nutrition, terpenoids, and feeding have a significant impact on the expression levels of the Dafem-1A, Dafem-1B, and Dafem-1C genes. This provides a crucial theoretical foundation for further elucidating the sex determination mechanism of Dendroctonus armandi.

Genome-wide identification and expression profiles of sex-related gene families in the Pacific abalone Haliotis discus hannai

In recent years, members of the Dmrt family, TGF-β superfamily and Sox family have been recognized as crucial genes for sex determination/differentiation across diverse animal species. Nevertheless, knowledge regarding the abundance and potential functions of these genes in abalone remains limited. In this study, a total of 5, 10, and 7 members of the Dmrt family, the TGF-β superfamily and the Sox family, respectively, were identified in the Pacific abalone Haliotis discus hannai. Sequence characteristics, phylogenetic relationships and spatiotemporal expression profiles of these genes were investigated. Notably, HdDmrt-04 (Dmrt1/1L-like) emerged as a potential mollusc-specific gene with a preponderance for expression in the testis. Interestingly, none of the TGF-β superfamily members exhibited specific or elevated expression in the gonads, highlighting the need for further investigation into their role in abalone sex differentiation. The Sox proteins in H. discus hannai were categorized into 7 subfamilies: B1, B2, C, D, E, F, and H. Among them, HdSox-07 (SoxH-like) was observed to play a crucial role in testis development, while HdSox-03 (SoxB1-like) and HdSox-04 (SoxC-like) probably cooperate in abalone ovary development. Taken together, the results of the present study suggested that HdDmrt-04 and HdSox-07 can be used as male-specific markers for gonad differentiation in H. discus hannai and imply conservation of their functions across invertebrates and vertebrates. Our findings provide new insights into the evolution and genetic structure of the Dmrt family, the TGF-β superfamily and the Sox family in abalone and pave the way for a deeper understanding of sex differentiation in gastropods.

Gonadal Transcriptome Sequencing Analysis Reveals the Candidate Sex-Related Genes and Signaling Pathways in the East Asian Common Octopus, Octopus sinensis

The East Asian common octopus (Octopus sinensis) is an economically important species among cephalopods. This species exhibits a strict dioecious and allogamous reproductive strategy, along with a phenotypic sexual dimorphism, where the third right arm differentiates into hectocotylus in males. However, our understanding of the molecular mechanisms that underlie sex determination and differentiation in this species remains limited. In the present study, we surveyed gene-expression profiles in the immature male and female gonads of O. sinensis based on the RNA-seq, and a total of 47.83 Gb of high-quality data were generated. Compared with the testis, we identified 8302 differentially expressed genes (DEGs) in the ovary, of which 4459 genes were up-regulated and 3843 genes were down-regulated. Based on the GO enrichment, many GO terms related to sex differentiation were identified, such as sex differentiation (GO: 0007548), sexual reproduction (GO: 0019953) and male sex differentiation (GO: 0046661). A KEGG classification analysis identified three conserved signaling pathways that related to sex differentiation, including the Wnt signaling pathway, TGF-β signaling pathway and Notch signaling pathway. Additionally, 21 sex-related DEGs were selected, of which 13 DEGs were male-biased, including Dmrt1, Foxn5, Foxj1, Sox30, etc., and 8 DEGs were female-biased, including Sox14, Nanos3, β-tubulin, Suh, etc. Ten DEGs were used to verify the expression patterns in the testis and ovary using the RT-qPCR method, and the results showed that the expression level shown by RT-qPCR was consistent with that from the RNA-seq, which confirmed the reliability of the transcriptome data. The results presented in this study will not only contribute to our understanding of sex-formation mechanisms in O. sinensis but also provide the foundational information for further investigating the molecular mechanisms that underline its gonadal development and facilitate the sustainable development of octopus artificial breeding.

A Testis-Specific DMRT1 (Double Sex and Mab-3-Related Transcription Factor 1) Plays a Role in Spermatogenesis and Gonadal Development in the Hermaphrodite Boring Giant Clam Tridacna crocea

The testis-specific double sex and mab-3-related transcription factor 1 (DMRT1) has long been recognized as a crucial player in sex determination across vertebrates, and its essential role in gonadal development and the regulation of spermatogenesis is well established. Here, we report the cloning of the key spermatogenesis-related DMRT1 cDNA, named Tc-DMRT1, from the gonads of Tridacna crocea (T. crocea), with a molecular weight of 41.93 kDa and an isoelectric point of 7.83 (pI). Our hypothesis is that DMRT1 machinery governs spermatogenesis and regulates gonadogenesis. RNAi-mediated Tc-DMRT1 knockdown revealed its critical role in hindering spermatogenesis and reducing expression levels in boring giant clams. A histological analysis showed structural changes, with normal sperm cell counts in the control group (ds-EGFP) but significantly lower concentrations of sperm cells in the experimental group (ds-DMRT1). DMRT1 transcripts during embryogenesis exhibited a significantly high expression pattern (p < 0.05) during the early zygote stage, and whole-embryo in-situ hybridization confirmed its expression pattern throughout embryogenesis. A qRT-PCR analysis of various reproductive stages revealed an abundant expression of Tc-DMRT1 in the gonads during the male reproductive stage. In-situ hybridization showed tissue-specific expression of DMRT1, with a positive signal detected in male-stage gonadal tissues comprising sperm cells, while no signal was detected in other stages. Our study findings provide an initial understanding of the DMRT1 molecular machinery controlling spermatogenesis and its specificity in male-stage gonads of the key bivalve species, Tridacna crocea, and suggest that DMRT1 predominantly functions as a key regulator of spermatogenesis in giant clams.

Gonadal transcriptome analysis of paradise fish Macropodus opercularis to reveal sex-related genes

Macropodus opercularis is an ornamental fish species endemic to China, with obvious sexual dimorphism in phenotype. To obtain the gene expression profile of the gonads of M. opercularis and explore its sex-related genes, six cDNA libraries were constructed from the sexually mature M. opercularis, and RNA-seq analysis was performed. The sequenced clean data were assembled by de novo splicing to generate 171,415 unigenes, and differentially expressed genes (DEGs) screening revealed that there were 41,638 DEGs in the gonads of M. opercularis. By comparing those DEGS in the ovary with the testis, we found 29,870 DEGs were upregulated and 11,768 DEGs were downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis showed that GO terms related to cell cycle and gamete formation were enriched, and pathway signals related to sex differences, such as FoxO signalling pathway and PI3K-Akt signalling pathway, were also detected. Reverse transcript fluorescence quantitative PCR (RT-qPCR) validation of 14 DEGs associated with sex differences showed that the RT-qPCR results were consistent with RNA-Seq analysis, and five genes, foxl2, sox3, foxo, zar1, cyp19a1, were significantly expressed in the ovaries. dmrt1, cyp11b, amh, sf1, sox9, gdf6, dmrt3, fstl1 and hsd11b2, a total of nine genes were significantly expressed in the testis. The results of this study provide a basis for the study of gonadal differentiation, developmental mechanisms and related functional genes in M. opercularis.

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

Reproductive condition of the black-lip pearl oyster Pinctada margaritifera during the lunar phase

This study analyzed the relationship between the lunar phase and the reproductive cycle of Pinctada margaritifera inhabiting Weno Island, Chuuk Lagoon, Micronesia. We measured indicators of maturity (gonadosomatic index [GSI] and sexual maturation-related genes) and investigated changes in the gonadal maturity stages (GMS) of P. margaritifera over lunar cycle. GSI was higher around the full moon. GMS of P. margaritifera were classified as the early gametogenesis stage, ripe and spawning stage, and spent and degenerating stage. A large percentage of oysters was observed in the ripe and spawning stage at the first quarter moon in female and the full moon in male as well as in the spent and degenerating stages at the third quarter moon in both sexes. In addition, the expression of doublesex- and mab-3-related transcription factor 2 (DMRT2) in the male P. margaritifera black-lip pearl oyster was the highest during the full and third quarter moon phases, whereas no difference in expression was observed with the lunar phase in females. In contrast, the expression of vitellogenin (VTG) was the highest in female P. margaritifera during the first and third quarters. No difference in expression was observed according to the lunar phase in males. The results suggest that the lunar phase directly affects the expression of sexually mature gonads in P. margaritifera black-lip pearl oyster.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.

DNA methylation differences between male and female gonads of the oyster reveal the role of epigenetics in sex determination

DNA methylation involved in sex determination mechanism by regulating gene expression related to sex determination networks are common in vertebrates. However, the mechanism linking epigenetics in invertebrates and sex determination has remained elusive. Here, methylome of the male and female gonads in the oyster Crassostrea gigas were conducted to explore the role of epigenetics in invertebrate sex determination. Comparative analysis of gonadal DNA methylation of females and males revealed that male gonads displayed a higher level of DNA methylation and a greater number of hypermethylated genes. Luxury genes presented hypomethylation, while housekeeping genes got hypermethylation. Genes in the conserved signaling pathways, rather than the key master genes in the sex determination pathway, were the major targets of substantial DNA methylation modification. The negative correlation of expression and promoter methylation in the diacylglycerol kinase delta gene (Dgkd) - a ubiquitously expressed gene - indicated DNA methylation may fine turn the expression of Dgkd and be involved in the process of sex determination. Dgkd can be used as an epigenetic marker to distinguish male C. gigas based on the different methylation regions in the promoter region. The results suggest that DNA methylation mechanisms played potential functional impacts in the sex determination in oysters, which is helpful to deepen the understanding of sex determination in invertebrate.

Full-Length Transcriptome Sequences Provide Insight Into Hermaphroditism of Freshwater Pearl Mussel Hyriopsis schlegelii

The freshwater mussel Hyriopsis schlegelii is a cultured bivalve in China, and the quality of the pearls produced is affected by the type of gonads. However, because of the lack of a published genome and the complexity of sex determination, research on sex reversal and development of this species is limited. In this study, Illumina RNA-seq and PacBio Isoform Sequencing (Iso-Seq) were combined to analyze the gonads of H. schlegelii. A total of 201,481 high-quality transcripts were generated. The study identified 7,922 differentially expressed genes in three comparison group (females versus males, hermaphrodites versus females, and hermaphrodites versus males). Twenty-four genes were identified as potential sex-related genes, including sox9 and wnt4 involved in sex determination, and vtg, cyp17a1 and 17β-hsd2 involved in gonadal development. We also speculated a possible pathways for the formation of hermaphroditism in H. schlegelii. The data provide a clear view of the transcriptome for H. schlegelii gonads and will be valuable in elucidating the mechanisms of gonad development.

Proteome and Transcriptome Analysis of Gonads Reveals Intersex in Gigantidas haimaensis

Sex has proven to be one of the most intriguing areas of research across evolution, development, and ecology. Intersex or sex change occurs frequently in molluscs. The deep-sea mussel Gigantidas haimaensis often dominates within Haima cold seep ecosystems, but details of their reproduction remain unknown. Herein, we conducted a combined proteomic and transcriptomic analysis of G. haimaensis gonads to provide a systematic understanding of sexual development in deep-sea bivalves. A total of 2,452 out of 42,238 genes (5.81%) and 288 out of 7,089 proteins (4.06%) were significantly differentially expressed between ovaries and testes with a false discovery rate (FDR) <0.05. Candidate genes involved in sexual development were identified; among 12 differentially expressed genes between sexes, four ovary-biased genes (β-catenin, fem-1, forkhead box L2 and membrane progestin receptor α) were expressed significantly higher in males than females. Combining histological characteristics, we speculate that the males maybe intersex undergoing sex change, and implied that these genes may be involved in the process of male testis converting into female gonads in G. haimaensis. The results suggest that this adaptation may be based on local environmental factors, sedentary lifestyles, and patchy distribution, and sex change may facilitate adaptation to a changing environment and expansion of the population. The findings provide a valuable genetic resource to better understand the mechanisms of sex change and survival strategies in deep-sea bivalves.

Comparative phylotranscriptomics reveals putative sex differentiating genes across eight diverse bivalve species

Mollusks, especially bivalves, exhibit a great diversity of sex determining mechanisms, including both genetic and environmental sex determination. Some bivalve species can be gonochoristic (separate sexes), while others are hermaphroditic (sequential or simultaneous). Several models have been proposed for specific bivalve species, utilizing information gained from gene expression data, as well as limited RAD-seq data (e.g., from Crassostrea gigas). However, these mechanisms are not as well studied as those in model organisms (e.g., Mus musculus, Drosophila melanogaster, Caenorhabditis elegans) and many genes involved in sex differentiation are not well characterized. We used phylotranscriptomics to better understand which possible sex differentiating genes are in bivalves and how these genes relate to similar genes in diverse phyla. We collected RNAseq data from eight phylogenetically diverse bivalve species: Argopecten irradians, Ensis directus, Geukensia demissa, Macoma tenta, Mercenaria mercenaria, Mya arenaria, Mytilus edulis, and Solemya velum. Using these data, we assembled representative transcriptomes for each species. We then searched for candidate sex differentiating genes using BLAST and confirmed the identity of nine genes using phylogenetics analyses from nine phyla. To increase the confidence of identification, we included ten bivalve genomes in our analyses. From the analysis of doublesex and mab-3 related transcription factor (DMRT) genes, we confirmed the identify of a Mollusk-specific sex determining DMRT gene: DMRT1L. Based on gene expression data from M. edulis and previous research, DMRT1L and FoxL2 are key genes for male and female development, respectively.

Gonadal transcriptomes associated with sex phenotypes provide potential male and female candidate genes of sex determination or early differentiation in Crassostrea gigas, a sequential hermaphrodite mollusc

Background

In the animal kingdom, mollusca is an important phylum of the Lophotrochozoa. However, few studies have investigated the molecular cascade of sex determination/early gonadal differentiation within this phylum. The oyster Crassostrea gigas is a sequential irregular hermaphrodite mollusc of economic, physiological and phylogenetic importance. Although some studies identified genes of its sex-determining/−differentiating pathway, this particular topic remains to be further deepened, in particular with regard to the expression patterns. Indeed, these patterns need to cover the entire period of sex lability and have to be associated to future sex phenotypes, usually impossible to establish in this sequential hermaphrodite. This is why we performed a gonadal RNA-Seq analysis of diploid male and female oysters that have not changed sex for 4 years, sampled during the entire time-window of sex determination/early sex differentiation (stages 0 and 3 of the gametogenetic cycle). This individual long-term monitoring gave us the opportunity to explain the molecular expression patterns in the light of the most statistically likely future sex of each oyster.

Results

The differential gene expression analysis of gonadal transcriptomes revealed that 9723 genes were differentially expressed between gametogenetic stages, and 141 between sexes (98 and 43 genes highly expressed in females and males, respectively). Eighty-four genes were both stage- and sex-specific, 57 of them being highly expressed at the time of sex determination/early sex differentiation. These 4 novel genes including Trophoblast glycoprotein-like, Protein PML-like, Protein singed-like and PREDICTED: paramyosin, while being supported by RT-qPCR, displayed sexually dimorphic gene expression patterns.

Conclusions

This gonadal transcriptome analysis, the first one associated with sex phenotypes in C. gigas, revealed 57 genes highly expressed in stage 0 or 3 of gametogenesis and which could be linked to the future sex of the individuals. While further study will be needed to suggest a role for these factors, some could certainly be original potential actors involved in sex determination/early sex differentiation, like paramyosin and could be used to predict the future sex of oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07838-1.

RNA-seq analysis of the mantle transcriptome from Mytilus edulis during a seasonal spawning event in deep and shallow water culture sites on the northeast coast of Newfoundland, Canada

The blue mussel (Mytilus edulis) has global commercial and ecological importance both in wild and cultured conditions. However there is a qualitative and quantitative lack of knowledge of the molecular mechanisms associated with its reproductive physiology, especially with reference to environmental interactions. Here we initiated a transcriptomic analysis (RNA-sequencing (RNA-seq)) of the mantle from both sexes sampled during a seasonal spawning event and from two culture depths (shallow-5 m; deep- 15 m). Mantle libraries were produced from 3 males and 3 females sampled from each of two shallow sites and two deep sites for a total of 12 replicate male and 12 replicate female libraries (24 total libraries). Overall a total of 2.3 billion raw 100 base reads with an average of 96.5 million reads/library were obtained and assembled into 296,118 transcripts with an average length of 568 bp. Overall, 315 transcripts from male libraries and 25 from female libraries were found to be upregulated in deep water as compared to shallow (edgeR adjusted p value ≤ 0.05). Conversely, 126 transcripts from male libraries and 135 from female libraries were found to be significantly downregulated at the same depth. Thirteen transcripts were selected for qPCR validation based on importance in reproduction, antimicrobial defense and metabolism. Of these, 9 RNA-seq identified transcripts were shown by qPCR to be differentially expressed between groups: 2 were upregulated in deep compared with shallow water (dhx38, mt-co1), 2 were upregulated for female compared with male mantle (pias2, mapkap1) and 6 genes (fndc3a, acbd3, klhl10, ccnb3, armc4, mt-co1) showed to be upregulated in males compared to females. The majority of qPCR studied transcripts were identified as involved in gamete development based on the UniProt database. This study further characterizes the importance of the mantle transcriptome during reproductive activities of M. edulis.

Sexual Development of the Hermaphroditic Scallop Argopecten irradians Revealed by Morphological, Endocrine and Molecular Analysis

Simultaneous or functional hermaphrodites possessing both ovary and testis at the same time are good materials for studying sexual development. However, previous research on sex determination and differentiation was mainly conducted in gonochoristic species and studies on simultaneous hermaphrodites are still limited. In this study, we conducted a combined morphological, endocrine and molecular study on the gonadal development of a hermaphroditic scallop Argopecten irradians aged 2–10 month old. Morphological analysis showed that sex differentiation occurred at 6 months of age. By examining the dynamic changes of progesterone, testosterone and estradiol, we found testosterone and estradiol were significantly different between the ovaries and testes almost throughout the whole process, suggesting the two hormones may be involved in scallop sex differentiation. In addition, we identified two critical sex-related genes FoxL2 and Dmrt1L, and investigated their spatiotemporal expression patterns. Results showed that FoxL2 and Dmrt1L were female- and male-biased, respectively, and mainly localized in the germ cells and follicular cells, indicating their feasibility as molecular markers for early identification of sex. Further analysis on the changes of FoxL2 and Dmrt1L expression in juveniles showed that significant sexual dimorphic expression of FoxL2 occurred at 2 months of age, earlier than that of Dmrt1L. Moreover, FoxL2 expression was significantly correlated with estradiol/testosterone ratio (E₂/T). All these results indicated that molecular sex differentiation occurs earlier than morphological sex differentiation, and FoxL2 may be a key driver that functions through regulating sex steroid hormones in the scallop. This study will deepen our understanding of the molecular mechanism underlying sex differentiation and development in spiralians.

Comparative transcriptome analysis of three gonadal development stages reveals potential genes involved in gametogenesis of the fluted giant clam (Tridacna squamosa)

BACKGROUND

Gonad development and differentiation is an essential function for all sexually reproducing species, and many aspects of these developmental processes are highly conserved among the metazoa. However, the mechanisms underlying gonad development and gametogenesis remain unclear in Tridacna squamosa, a large-size bivalve of great ecological value. They are protandrous simultaneous hermaphrodites, with the male gonad maturing first, eventually followed by the female gonads. In this study, nine gonad libraries representing resting, male and hermaphrodite stages in T. squamosa were performed to identify the molecular mechanisms.

RESULTS

Sixteen thousand four hundred ninety-one unigenes were annotated in the NCBI non-redundant protein database. Among the annotated unigenes, 5091 and 7328 unigenes were assigned to Gene Ontology categories and the Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway database, respectively. A total of 4763 differentially expressed genes (DEGs) were identified by comparing male to resting gonads, consisting of 3499 which were comparatively upregulated in males and 1264 which were downregulated in males. Six hundred-ninteen DEGs between male and hermaphroditic gonads were identified, with 518 DEGs more strongly expressed in hermaphrodites and 101 more strongly expressed in males. GO (Gene Ontology) and KEGG pathway analyses revealed that various biological functions and processes, including functions related to the endocrine system, oocyte meiosis, carbon metabolism, and the cell cycle, were involved in regulating gonadal development and gametogenesis in T. squamosa. Testis-specific serine/threonine kinases 1 (TSSK1), TSSK4, TSSK5, Doublesex- and mab-3-related transcription factor 1 (DMRT1), SOX, Sperm surface protein 17 (SP17) and other genes were involved in male gonadal development in Tridacna squamosal. Both spermatogenesis- (TSSK4, spermatogenesis-associated protein 17, spermatogenesis-associated protein 8, sperm motility kinase X, SP17) and oogenesis-related genes (zona pellucida protein, Forkhead Box L2, Vitellogenin, Vitellogenin receptor, 5-hydroxytryptamine, 5-hydroxytryptamine receptor) were simultaneously highly expressed in the hermaphroditic gonad to maintain the hermaphroditism of T. squamosa.

CONCLUSION

All these results from our study will facilitate better understanding of the molecular mechanisms underlying giant clam gonad development and gametogenesis, which can provided a base on obtaining excellent gametes during the seed production process for giant clams.

Characterizing transcriptome in female scallop Chlamys farreri provides new insights into the molecular mechanisms of reproductive regulation during ovarian development and spawn

Bivalve mollusks are descendants of an early-Cambrian lineage and have successfully evolved unique strategies for reproduction. Nonetheless, the molecular mechanisms underlying reproductive regulation in mollusks remain to be elucidated. In this study, transcriptomes of ovary at four reproductive stages in female Chlamys farreri were characterized by RNA-Seq. Regarding signaling pathways, ECM-receptor interaction pathway, mTOR signaling pathway, Fanconi anemia pathway, FoxO signaling pathway, Wnt signaling pathway and Hedgehog signaling pathway were enriched during ovarian development processes. In addition, pathways related to energy metabolism such as Nitrogen metabolism and Arachidonic acid metabolism were enriched at spawn stage. Interestingly, Neuroactive ligand-receptor interaction was significantly enriched involved in ovarian development and spawn, and indicated the potential functions of nervous system on reproductive regulation in C. farreri. What's more, this study identified and characterized fourteen genes involved in "sex hormones synthesis and regulation", "ovarian development and spawn" and "maternal immunity" during the four reproductive stages in C. farreri. We determined that CYP17 uniquely affected gamete release by influencing the physiological balance among the steroid hormones and showed that receptors of the 5-HT and GABA neurotransmitters were tightly associated with ovarian maturation. Furthermore, to the best of our knowledge, this is the first study to report the maternal effect gene Zar1 in bivalve mollusks, likewise the maternal immunity genes displayed coordinated and cooperative expression during reproductive periods, which strengthened the environmental adaptation mechanisms of bivalves. Taken together, this study provides the first dynamic transcriptomic analysis of C. farreri at four key reproductive stages, which will assist in revealing the molecular mechanisms underlying bivalves on reproductive regulation in ovarian development and spawn.

De novo transcriptome assembly of four organs of Collichthys lucidus and identification of genes involved in sex determination and reproduction

The spinyhead croaker (Collichthys lucidus) is a commercially important fish species, which is mainly distributed in the coastal regions of China. However, little is known about the molecular regulatory mechanism underlying reproduction in C. lucidus. A de novo transcriptome assembly in brain, liver, ovary and testis tissues of C. lucidus was performed. Illumina sequencing generated 60,322,004, 57,044,284, 60,867,978 and 57,087,688 clean reads from brain, liver, ovary and testis tissues of C. lucidus, respectively. Totally, 131,168 unigenes with an average length of 644 bp and an N50 value of 1033 bp were assembled. In addition, 1288 genes were differentially expressed between ovary and testis, including 442 up-regulated and 846 down-regulated in ovary. Functional analysis revealed that the differentially expressed genes between ovary and testis were mainly involved in the function of sexual reproduction, sex differentiation, development of primary male sexual characteristics, female gamete generation, and male sex differentiation. A number of genes which might be involved in the regulation of reproduction and sex determination were found, including HYAL and SYCP3 and BMP15. Furthermore, 35,476 simple sequence repeats (SSRs) were identified in this transcriptome dataset, which would contribute to further genetic and mechanism researches. De novo transcriptome sequencing analysis of four organs of C. lucidus provides rich resources for understanding the mechanism of reproductive development of C. lucidus and further investigation of the molecular regulation of sex determination and reproduction of C. lucidus.

Identification and expression profiles of Fox transcription factors in the Yesso scallop (Patinopecten yessoensis)

The forkhead box (Fox) gene family is a family of transcription factors that play important roles in a variety of biological processes in vertebrates, including early development and cell proliferation and differentiation. However, at present, studies on the mollusk Fox family are relatively lacking. In the present study, the Fox gene family of the Yesso scallop (Patinopecten yessoensis) was systematically identified. In addition, the expression profiles of the Fox gene family in early development and adult tissues were analyzed. The results showed that there were 26 Fox genes in P. yessoensis. Of the 26 genes, 24 belonged to 20 subfamilies. The Fox genes belonging to the I, Q1, R and S subfamilies were absent in P. yessoensis. The other 2 genes formed 2 independent clades with the Fox genes of other mollusks and protostomes. They might be new members of the Fox family and were named FoxY and FoxZ. P. yessoensis contained a FoxC-FoxL1 gene cluster similar in structure to that of Branchiostoma floridae, suggesting that the cluster might already exist in the ancestors of bilaterally symmetrical animals. The gene expression analysis of Fox showed that most of the genes were continuously expressed in multiple stages of early development, suggesting that Fox genes might be widely involved in the regulation of embryo and larval development of P. yessoensis. Nine Fox genes were specifically expressed in certain tissues, such as the nerve ganglia, foot, ovary, testis, and gills. For the 9 genes that were differentially expressed between the testis and ovary, their expression levels were analyzed during the 4 developmental stages of gonads. The results showed that FoxL2, FoxE and FoxY were highly expressed in the ovary during all developmental stages, while FoxZ was highly expressed in the testis during all developmental stages. The results suggested that these genes might play an important role in sex maintenance or gametogenesis. The present study could provide a reference for evolutionary and functional studies of the Fox family in metazoans.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

De novo gonad transcriptome analysis of the common littoral shrimp Palaemon serratus: novel insights into sex-related genes

BACKGROUND

The common littoral shrimp Palaemon serratus is an economically important decapod resource in some European communities. Aquaculture practices prevent the genetic deterioration of wild stocks caused by overfishing and at the same time enhance the production. The biotechnological manipulation of sex-related genes has the proved potential to improve the aquaculture production but the scarcity of genomic data about P. serratus hinders these applications. RNA-Seq analysis has been performed on ovary and testis samples to generate a reference gonadal transcriptome. Differential expression analyses were conducted between three ovary and three testis samples sequenced by Illumina HiSeq 4000 PE100 to reveal sex-related genes with sex-biased or sex-specific expression patterns.

RESULTS

A total of 224.5 and 281.1 million paired-end reads were produced from ovary and testis samples, respectively. De novo assembly of ovary and testis trimmed reads yielded a transcriptome with 39,186 transcripts. The 29.57% of the transcriptome retrieved at least one annotation and 11,087 differentially expressed genes (DEGs) were detected between ovary and testis replicates. Six thousand two hundred seven genes were up-regulated in ovaries meanwhile 4880 genes were up-regulated in testes. Candidate genes to be involved in sexual development and gonadal development processes were retrieved from the transcriptome. These sex-related genes were discussed taking into account whether they were up-regulated in ovary, up-regulated in testis or not differentially expressed between gonads and in the framework of previous findings in other crustacean species.

CONCLUSIONS

This is the first transcriptome analysis of P. serratus gonads using RNA-Seq technology. Interesting findings about sex-related genes from an evolutionary perspective (such as Dmrt1) and for putative future aquaculture applications (Iag or vitellogenesis genes) are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of this shrimp.

Hepatopancreas and ovarian transcriptome response to different dietary soybean lecithin levels in Portunus trituberculatus

Ovaries (O) are specialized tissues that play critical roles in producing oocytes and hormones. The crustacean hepatopancreas (H) is a metabolic organ that plays important functions including absorption, storage of nutrients and vitellogenesis during growth and ovarian development. However, genetic information on the biological functions of the crustacean ovaries and hepatopancreas are limited. This study compared the transcriptome in the ovary and the hepatopancreas of female P. trituberculatus fed two different diets containing 0% (SL0) and 4% soybean lecithin (SL4), respectively during the growth and ovarian maturation stages by Illumina HiSeq4000 sequencer. The differences between ovary and hepatopancreas of P. trituberculatus were also compared at transcriptional level. A total of 55,667 unigenes were obtained with mean length of 962 bps across the four treatment groups (SL0_O, SL4_O, SL0_H and SL4_H). In ovary, there were 257 differentially expressed genes (DEGs) between SL0_O and SL4_O, with 145 down- and 112 up-regulated genes in the SL4_O group. Candidate genes involved in ovarian development were detected in SL4_H group. In hepatopancreas, 146 DEGs were found between SL0_H and SL4_H, including 43 down- and 103 up-regulated genes in the SL4_H group. The specific DEGs were mainly involved with lipid related metabolism pathways, including fat digestion and absorption, PPAR signaling pathway and insulin resistance. 14,725 DEGs were found in the comparison between SL0_O and SL4_H, including 7250 up- and 7475 down-regulated genes in the SL4_H group. The specific DEGs were mainly involved with lipid (fat digestion and absorption, linoleic acid metabolism), hormone (steroid hormone biosynthesis, ovarian steroidogenesis, etc), and amino acid (phenylalanine metabolism, arginine biosynthesis, tyrosine) related metabolism pathways. Crabs fed the SL4 diet exhibited higher gene expression of cryptocyanin 1 (cc1), cryptocyanin 2 (cc2) and neuroparsin 1 (np1) in hepatopancreas and ovarian than those fed the SL0 diet, however, crab fed SL4 diet showed higher gene expression of fatty acid-binding protein 1 (fabp1), vitellogenin (vtg) and Delta-6 desaturase-like protein (fadsd6) in hepatopancreas than those fed the SL0 diet. Moreover, crabs fed the SL0 diet had lower gene expression of vtg, extracellular copper‑zinc superoxide dismutase (cuznsod) and estrogen sulfotransferase (ests) in ovary compared to those fed the diet containing 4% soybean lecithin. These results might provide important clues with respect to elucidating the molecular mechanisms underlying the regulation of phospholipid on the gonadal development and lipid metabolism of P. trituberculatus.

Transcriptomics analysis revealing candidate genes and networks for sex differentiation of yesso scallop (Patinopecten yessoensis)

Background

The Yesso scallop, Patinopecten (Mizuhopecten) yessoensis, is a commercially important bivalve in the coastal countries of Northeast Asia. It has complex modes of sex differentiation, but knowledge of the mechanisms underlying this sex determination and differentiation is limited.

Results

In this study, the gonad tissues from females and males at three developmental stages were used to investigate candidate genes and networks for sex differentiation via RNA-Req. A total of 901,980,606 high quality clean reads were obtained from 18 libraries, of which 417 expressed male-specific genes and 754 expressed female-specific genes. Totally, 10,074 genes differentially expressed in females and males were identified. Weighted gene co-expression network analysis (WGCNA) revealed that turquoise and green gene modules were significantly positively correlated with male gonads, while coral1 and black modules were significantly associated with female gonads. The most important gene for sex determination and differentiation was Pydmrt 1, which was the only gene discovered that determined the male sex phenotype during early gonadal differentiation. Enrichment analyses of GO terms and KEGG pathways revealed that genes involved in metabolism, genetic and environmental information processes or pathways are sex-biased. Forty-nine genes in the five modules involved in sex differentiation or determination were identified and selected to construct a gene co-expression network and a hypothesized sex differentiation pathway.

Conclusions

The current study focused on screening genes of sex differentiation in Yesso scallop, highlighting the potential regulatory mechanisms of gonadal development in P. yessoensis. Our data suggested that WCGNA can facilitate identification of key genes for sex differentiation and determination. Using this method, a hypothesized P. yessoensis sex determination and differentiation pathway was constructed. In this pathway, Pydmrt 1 may have a leading function.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6021-6) contains supplementary material, which is available to authorized users.

Molecular Cloning and Sexually Dimorphic Expression Analysis of nanos2 in the Sea Urchin, Mesocentrotus nudus

Sea urchin (Mesocentrotus nudus) is an economically important mariculture species in China and the gonads are the solely edible parts to human. The molecular mechanisms of gonad development have attracted increasing attention in recent years. Although the nanos2 gene has been identified as a germ cell marker in several invertebrates, little is known about nanos2 in adult sea urchins. Hereinto, we report the characterization of Mnnano2, an M. nudus nanos2 homology gene. Mnnanos2 is a maternal factor and can be detected continuously during embryogenesis and early ontogeny. Real-time quantitative PCR (RT-qPCR) and section in situ hybridization (ISH) analysis revealed a dynamic and sexually dimorphic expression pattern of Mnnano2 in the gonads. Its expression reached the maximal level at Stage 2 along with the gonad development in both ovary and testis. In the ovary, Mnnanos2 is specifically expressed in germ cells. In contrast, Mnnanos2 is expressed in both nutritive phagocytes (NP) cells and male germ cells in testis. Moreover, knocking down of Mnnanos2 by means of RNA interference (RNAi) reduced nanos2 and boule expression but conversely increased the expression of foxl2. Therefore, our data suggest that Mnnanos2 may serve as a female germ cell marker during gametogenesis and provide chances to uncover its function in adult sea urchin.

Phenotypic Stability of Sex and Expression of Sex Identification Markers in the Adult Yesso Scallop Mizuhopecten yessoensis throughout the Reproductive Cycle

Simple Summary

Bivalve sex is thought to fluctuate depending on environmental conditions. So far, there has been no investigation on the phenotypic stability of sex in the commercially important Yesso scallop Mizuhopecten yessoensis. The present study revealed that the sex of the Yesso scallop is stable after initial sex differentiation and that this species maintains a sex-stable maturation system throughout its life. In addition, gonad differentiation for each sex was precisely characterized by using molecular markers throughout the maturational cycle.

Abstract

The objective of the present study was to analyze the phenotypic stability of sex after sex differentiation in the Yesso scallop, which is a gonochoristic species that has been described as protandrous. So far, no study has investigated in detail the sexual fate of the scallop after completion of sex differentiation, although bivalve species often show annual sex change. In the present study, we performed a tracking experiment to analyze the phenotypic stability of sex in scallops between one and two years of age. We also conducted molecular marker analyses to describe sex differentiation and gonad development. The results of the tracking experiment revealed that all scallops maintained their initial sex phenotype, as identified in the last reproductive period. Using molecular analyses, we characterized my-dmrt2 and my-foxl2 as sex identification markers for the testis and ovary, respectively. We conclude by proposing that the Yesso scallop is a sex-stable bivalve after its initial sex differentiation and that it maintains a sex-stable maturation system throughout its life. The sex-specific molecular markers identified in this study are useful tools to assess the reproductive status of the Yesso scallop.

Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis

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₂(TPM) of RNA-Seq data (Ct = − 0.94 log₂(TPM) + 29.67, R² = 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.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5661-x) contains supplementary material, which is available to authorized users.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

Transcriptome analysis for identifying possible causes of post-reproductive death of Sepia esculenta based on brain tissue

BACKGROUND

The subpeduncle lobe/olfactory lobe-optic gland axis is called the endocrine regulation center of cephalopods. However, little is known about the mechanism of the subpeduncle lobe/olfactory lobe-optic gland axis regulate the sexual maturation and post-reproductive death of Sepia esculenta Hoyle.

OBJECTIVES

The primary objective of this study was to provide basic information for revealing the mechanism of the subpeduncle lobe/olfactory lobe-optic axis regulating the rapid post-reproductive death of S. esculenta.

METHODS

In this paper, Illumina sequencing based transcriptome analysis was performed on the brain tissue of female S. esculenta in the three key developmental stages: growth stage (BG), spawning stage (BS), and post-reproductive death stage (BA).

RESULTS

A total of 66.19 Gb Illumina sequencing data were obtained. A comparative analysis of the three stages showed 2609, 3333, and 170 differentially expressed genes (DEGs) in BG-vs-BA, BG-vs-BA, and BS-vs-BA, respectively. The Gene Ontology (GO) enrichment analysis of DEGs revealed that the regulation of cyclin-dependent protein serine/threonine kinase activity, oxidative phosphorylation, and respiratory chain were significantly enriched. The significant enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway identified pathways associated with the regulation of death, such as the mammalian target of rapamycin (mTOR) signaling pathway, AMPK signaling pathway, oxidative phosphorylation, and cell cycle.

CONCLUSION

The post-reproductive death of S. esculenta was found to be a complex energy steady-state regulation network system. The mTOR acted as an energy receptor and had a key role in regulating energy homeostasis.

Transcriptome analysis of air-breathing land slug, Incilaria fruhstorferi reveals functional insights into growth, immunity, and reproduction

Background

Incilaria (= Meghimatium) fruhstorferi is an air-breathing land slug found in restricted habitats of Japan, Taiwan and selected provinces of South Korea (Jeju, Chuncheon, Busan, and Deokjeokdo). The species is on a decline due to depletion of forest cover, predation by natural enemies, and collection. To facilitate the conservation of the species, it is important to decide on a number of traits related to growth, immunity and reproduction addressing fitness advantage of the species.

Results

The visceral mass transcriptome of I. fruhstorferi was enabled using the Illumina HiSeq 4000 sequencing platform. According to BUSCO (Benchmarking Universal Single-Copy Orthologs) method, the transcriptome was considered complete with 91.8% of ortholog genes present (Single: 70.7%; Duplicated: 21.1%). A total of 96.79% of the raw read sequences were processed as clean reads. TransDecoder identified 197,271 contigs that contained candidate-coding regions. Of a total of 50,230 unigenes, 34,470 (68.62% of the total unigenes) annotated to homologous proteins in the Protostome database (PANM-DB). The GO term and KEGG pathway analysis indicated genes involved in metabolism, phosphatidylinositol signalling system, aminobenzoate degradation, and T-cell receptor signalling pathway. Many genes associated with molluscan innate immunity were categorized under pathogen recognition receptor, TLR signalling pathway, MyD88 dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related. The reproduction-associated unigenes showed homology to protein fem-1, spermatogenesis-associated protein, sperm associated antigen, and testis expressed sequences, among others. In addition, we identified key growth-related genes categorized under somatotrophic axis, muscle growth, chitinases and collagens. A total of 4822 Simple Sequence Repeats (SSRs) were also identified from the unigene sequences of I. fruhstorferi.

Conclusions

This is the first available genomic information for non-model land slug, I. fruhstorferi focusing on genes related to growth, immunity, and reproduction, with additional focus on microsatellites and repeating elements. The transcriptome provides access to greater number of traits of unknown relevance in the species that could be exploited for in-depth analyses of evolutionary plasticity and making informed choices during conservation planning. This would be appropriate for understanding the dynamics of the species on a priority basis considering the ecological, health, and social benefits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5526-3) contains supplementary material, which is available to authorized users.

Whole transcriptome sequencing and biomineralization gene architecture associated with cultured pearl quality traits in the pearl oyster, Pinctada margaritifera

Background

Cultured pearls are unique gems produced by living organisms, mainly molluscs of the Pinctada genus, through the biomineralization properties of pearl sac tissue. Improvement of P. margaritifera pearl quality is one of the biggest challenges that Polynesian research has faced to date. To achieve this goal, a better understanding of the complex mechanisms related to nacre and pearl formation is essential and can now be approached through the use of massive parallel sequencing technologies. The aim of this study was to use RNA-seq to compare whole transcriptome expression of pearl sacs that had producing pearls with high and low quality. For this purpose, a comprehensive reference transcriptome of P. margaritifera was built based on multi-tissue sampling (mantle, gonad, whole animal), including different living stages (juvenile, adults) and phenotypes (colour morphotypes, sex).

Results

Strikingly, few genes were found to be up-regulated for high quality pearls (n = 16) compared to the up-regulated genes in low quality pearls (n = 246). Biomineralization genes up-regulated in low quality pearls were specific to prismatic and prism-nacre layers. Alternative splicing was further identified in several key biomineralization genes based on a recent P. margaritifera draft genome.

Conclusion

This study lifts the veil on the multi-level regulation of biomineralization genes associated with pearl quality determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5443-5) contains supplementary material, which is available to authorized users.

SOX2 participates in spermatogenesis of Zhikong scallop Chlamys farreri

As an important transcription factor, SOX2 involves in embryogenesis, maintenance of stem cells and proliferation of primordial germ cell (PGC). However, little was known about its function in mature gonads. Herein, we investigated the SOX2 gene profiles in testis of scallop, Chlamys farreri. The level of C. farreri SOX2 (Cf-SOX2) mRNA increased gradually along with gonadal development and reached the peak at mature stage, and was located in all germ cells, including spermatogonia, spermatocytes, spermatids and spermatozoa. Knockdown of Cf-SOX2 using RNAi leaded to a mass of germ cells lost, and only a few spermatogonia retained in the nearly empty testicular acini after 21 days. TUNEL assay showed that apoptosis occurred in spermatocytes. Furthermore, transcriptome profiles of the testes were compared between Cf-SOX2 knockdown and normal scallops, 131,340 unigenes were obtained and 2,067 differential expression genes (DEGs) were identified. GO and KEGG analysis showed that most DEGs were related to cell apoptosis (casp2, casp3, casp8), cell proliferation (samd9, crebzf, iqsec1) and spermatogenesis (htt, tusc3, zmynd10, nipbl, mfge8), and enriched in p53, TNF and apoptosis pathways. Our study revealed Cf-SOX2 is essential in spermatogenesis and testis development of C. farreri and provided important clues for better understanding of SOX2 regulatory mechanisms in bivalve testis.

Transcriptomic analysis of the red and green light responses in Columba livia domestica

In this study, 108 paired White King pigeons, randomly divided into three compartments were exposed to green light, red light, and white light followed by 15 h of light exposure, for a 6-month period. Three female birds from each group were selected and ovarian stromal tissue was collected. Pigeon reproductive data were also recorded every day. We performed transcriptome assembly on several tissue samples using Illumina Hiseq 2000 and analyzed differentially expressed genes involving follicle development mechanisms. Reproductive data confirmed that exposure to red and green lights improved pigeon reproduction. In total, approximately 158,080 unigenes with an average length of 753 bp were obtained using the Trinity program. Gene ontology, clusters of orthologous groups, and the Kyoto encyclopedia of genes were used to annotate and classify these unigenes. Large numbers of differentially expressed genes were discovered through pairwise comparisons between groups treated with monochromatic light versus white light. Some of these genes are associated with steroid hormone biosynthesis, cell cycle and circadian rhythm. Furthermore, qRT-PCR was used to detect the relative expression levels of randomly selected genes. A total of 17,419 potential simple sequence repeats were also identified. Our study provides insights into potential molecular mechanisms and genes that regulate pigeon reproduction in response to monochromatic light exposure. Our results and data will facilitate a further investigation into the molecular mechanisms behind the effects of red and green lights on follicle development and reproduction in the pigeon.

Sex Determination, Sexual Development, and Sex Change in Slipper Snails

Sex determination and sexual development are highly diverse and controlled by mechanisms that are extremely labile. While dioecy (separate male and female functions) is the norm for most animals, hermaphroditism (both male and female functions within a single body) is phylogenetically widespread. Much of our current understanding of sexual development comes from a small number of model systems, limiting our ability to make broader conclusions about the evolution of sexual diversity. We present the calyptraeid gastropods as a model for the study of the evolution of sex determination in a sequentially hermaphroditic system. Calyptraeid gastropods, a group of sedentary, filter-feeding marine snails, are sequential hermaphrodites that change sex from male to female during their life span (protandry). This transition includes resorption of the penis and the elaboration of female genitalia, in addition to shifting from production of spermatocytes to oocytes. This transition is typically under environmental control and frequently mediated by social interactions. Males in contact with females delay sex change to transition at larger sizes, while isolated males transition more rapidly and at smaller sizes. This phenomenon has been known for over a century; however, the mechanisms that control the switch from male to female are poorly understood. We review here our current understanding of sexual development and sex determination in the calyptraeid gastropods and other molluscs, highlighting our current understanding of factors implicated in the timing of sex change and the potential mechanisms. We also consider the embryonic origins and earliest expression of the germ line and the effects of environmental contaminants on sexual development.

De novo male gonad transcriptome draft for the marine mussel Perumytilus purpuratus with a focus on its reproductive-related proteins

Perumytilus purpuratus is a marine mussel considered a bioengineer species with a broad distribution in the Pacific and Atlantic coast of South America. Studies have shown two geographically and genetically differentiated subpopulations at molecular level and in sperm morphological traits. To open avenues for molecular research on P. purpuratus, a global de novo transcriptome from gonadal tissue of mature males was sequenced using the Illumina platform. From a total of 126.38 million reads, 37,765 transcripts were successfully annotated. BUSCO analysis determined a level of 89% completeness for the assembled transcriptome. The functional gene ontology (GO) annotation indicated that, in terms of abundance, the transcripts related with molecular function were the most represented, followed by those related with biological process and cellular components. Additionally, a subset of GO annotations generated using the "sperm" term resulted in a total of 1,294 sequences where the biological process category was the more represented, with transcripts strongly associated to sperm-processes required for fertilization, and with processes where the sperm-egg interaction could be implicated. Our work will contribute to the evolutionary understanding of the molecular mechanisms related to tissue-specific functions. This work reports the first male gonad transcriptome for the mussel P. purpuratus, generating a useful transcriptomic resource for this species and other closely related mytilids.

Dimorphism of sex and gonad-development-related genes in male and female lined seahorse, Hippocampus erectus, based on transcriptome analyses

Seahorse is characterized by its male pregnancy and sex-role reversal. To better understand the sexual dimorphism of male and female seahorses based on essential genes, we performed systematic transcriptome studies for both genders. A total of 157,834,590 cleaned reads were obtained and assembled into 129,268 transcripts and 31,764 could be annotated. Results showed that 176 up-regulated and 391 down-regulated transcripts were identified in the male seahorses compared with those in females. Genes involved in sex differentiation, such as dmrt1, sox9, fem1 and vasa, were identified and characterized. Moreover, the essential genes involved in reproductive molecular pathway were identified and analyzed in seahorses. In conclusion, the present study provides an archive for the future systematic research on seahorse sex differentiation.

FOXL2 and DMRT1L Are Yin and Yang Genes for Determining Timing of Sex Differentiation in the Bivalve Mollusk Patinopecten yessoensis

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 LOG₁₀(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.

Proteome and Transcriptome Analysis of Ovary, Intersex Gonads, and Testis Reveals Potential Key Sex Reversal/Differentiation Genes and Mechanism in Scallop Chlamys nobilis

Bivalve mollusks exhibit hermaphroditism and sex reversal/differentiation. Studies generally focus on transcriptional profiling and specific genes related to sex determination and differentiation. Few studies on sex reversal/differentiation have been reported. A combination analysis of gonad proteomics and transcriptomics was conducted on Chlamys nobilis to provide a systematic understanding of sex reversal/differentiation in bivalves. We obtained 4258 unique peptides and 93,731 unigenes with good correlation between messenger RNA and protein levels. Candidate genes in sex reversal/differentiation were found: 15 genes differentially expressed between sexes were identified and 12 had obvious sexual functions. Three novel genes (foxl2, β-catenin, and sry) were expressed highly in intersex individuals and were likely involved in the control of gonadal sex in C. nobilis. High expression of foxl2 or β-catenin may inhibit sry and activate 5-HT receptor and vitellogenin to maintain female development. High expression of sry may inhibit foxl2 and β-catenin and activate dmrt2, fem-1, sfp2, sa6, Amy-1, APCP4, and PLK to maintain male function. High expression of sry, foxl2, and β-catenin in C. nobilis may be involved in promoting and maintaining sex reversal/differentiation. The downstream regulator may not be dimorphic expressed genes, but genes expressed in intersex individuals, males and females. Different expression patterns of sex-related genes and gonadal histological characteristics suggested that C. nobilis may change its sex from male to female. These findings suggest highly conserved sex reversal/differentiation with diverged regulatory pathways during C. nobilis evolution. This study provides valuable genetic resources for understanding sex reversal/differentiation (intersex) mechanisms and pathways underlying bivalve reproductive regulation.

Deciphering the Link between Doubly Uniparental Inheritance of mtDNA and Sex Determination in Bivalves: Clues from Comparative Transcriptomics

Bivalves exhibit an astonishing diversity of sexual systems and sex-determining mechanisms. They can be gonochoric, hermaphroditic or androgenetic, with both genetic and environmental factors known to determine or influence sex. One unique sex-determining system involving the mitochondrial genome has also been hypothesized to exist in bivalves with doubly uniparental inheritance (DUI) of mtDNA. However, the link between DUI and sex determination remains obscure. In this study, we performed a comparative gonad transcriptomics analysis for two DUI-possessing freshwater mussel species to better understand the mechanisms underlying sex determination and DUI in these bivalves. We used a BLAST reciprocal analysis to identify orthologs between Venustaconcha ellipsiformis and Utterbackia peninsularis and compared our results with previously published sex-specific bivalve transcriptomes to identify conserved sex-determining genes. We also compared our data with other DUI species to identify candidate genes possibly involved in the regulation of DUI. A total of ∼12,000 orthologous relationships were found, with 2,583 genes differentially expressed in both species. Among these genes, key sex-determining factors previously reported in vertebrates and in bivalves (e.g., Sry, Dmrt1, Foxl2) were identified, suggesting that some steps of the sex-determination pathway may be deeply conserved in metazoans. Our results also support the hypothesis that a modified ubiquitination mechanism could be responsible for the retention of the paternal mtDNA in male bivalves, and revealed that DNA methylation could also be involved in the regulation of DUI. Globally, our results suggest that sets of genes associated with sex determination and DUI are similar in distantly-related DUI species.

Sex determination and differentiation genes in a functional hermaphrodite scallop, Nodipecten subnodosus

The lion-paw, Nodipecten subnodosus is one of three scallop species commercially exploited on the west coast of the Peninsula of Baja California. Because nothing is known about sex determination and sexual differentiation in hermaphrodite scallops, in the present work, a global transcriptomic analysis was performed in two early developmental stages, settling eyed-larvae and spat, as well as in three tissues (undifferentiated gonad, digestive gland, and adductor muscle). Over 27 million Illumina paired-end reads were obtained through the MiSeq platform. After processing the reads a total of 243,774 transcripts were assembled with an N50 of 980 and an average length of 775nt. A total of 43,252 proteins were inferred and 36,103 transcripts had at least one homolog in the SwissProt database according to a blastx search. After differential expression analyses and GO annotations it was possible to identify several sex-related genes in the scallop, including one known to be involved in the sex determination pathway of the hermaphrodite model organism Caenorhabditis elegans, N. subnodosus-sex1 (Ns-sex1). Other interesting sex determination and differentiation genes were Ns-dmrta2, Ns-sox9, Ns-wnt4, Ns-doa, Ns-ovo, Ns-vir, among others. Most of these genes were mainly expressed in the testis region, suggesting their participation in male gonad region sex differentiation. These results represent the first available information on the genetics of sex determination and differentiation in a functional hermaphrodite scallop.

The Whole-Genome and Transcriptome of the Manila Clam (Ruditapes philippinarum)

The manila clam, Ruditapes philippinarum, is an important bivalve species in worldwide aquaculture including Korea. The aquaculture production of R. philippinarum is under threat from diverse environmental factors including viruses, microorganisms, parasites, and water conditions with subsequently declining production. In spite of its importance as a marine resource, the reference genome of R. philippinarum for comprehensive genetic studies is largely unexplored. Here, we report the de novo whole-genome and transcriptome assembly of R. philippinarum across three different tissues (foot, gill, and adductor muscle), and provide the basic data for advanced studies in selective breeding and disease control in order to obtain successful aquaculture systems. An approximately 2.56 Gb high quality whole-genome was assembled with various library construction methods. A total of 108,034 protein coding gene models were predicted and repetitive elements including simple sequence repeats and noncoding RNAs were identified to further understanding of the genetic background of R. philippinarum for genomics-assisted breeding. Comparative analysis with the bivalve marine invertebrates uncover that the gene family related to complement C1q was enriched. Furthermore, we performed transcriptome analysis with three different tissues in order to support genome annotation and then identified 41,275 transcripts which were annotated. The R. philippinarum genome resource will markedly advance a wide range of potential genetic studies, a reference genome for comparative analysis of bivalve species and unraveling mechanisms of biological processes in molluscs. We believe that the R. philippinarum genome will serve as an initial platform for breeding better-quality clams using a genomic approach.

Transcriptome analysis of sex-related genes in the blood clam Tegillarca granosa

BACKGROUND

Blood clams (Tegillarca granosa) are one of the most commercial shellfish in China and South Asia with wide distribution in Indo-Pacific tropical to temperate estuaries. However, recent data indicate a decline in the germplasm of this species. Furthermore, the molecular mechanisms underpinning reproductive regulation remain unclear and information regarding genetic diversity is limited. Understanding the reproductive biology of shellfish is important in interpreting their embryology development, reproduction and population structure. Transcriptome sequencing (RNA-seq) rapidly obtains genetic sequence information from almost all transcripts of a particular tissue and currently represents the most prevalent and effective method for constructing genetic expression profiles.

RESULTS

Non-reference RNA-seq, an Illumina HiSeq2500 Solexa system, and de novo assembly were used to construct a gonadal expression profile of the blood clam. A total of 63.75 Gb of clean data, with at least 89.46% of Quality30 (Q30), were generated which was then combined into 214,440 transcripts and 125,673 unigenes with a mean length of 1,122.63 and 781.30 base pairs (bp). In total, 27,325 genes were annotated by comparison with public databases. Of these, 2,140 and 2,070 differentially expressed genes (DEGs) were obtained (T05 T08 vs T01 T02 T04, T06 T07 vs T01 T02 T04; in which T01-T04 and T05-T08 represent biological replicates of individual female and male clams, respectively) and classified into two groups according to the evaluation of biological replicates. Then 35 DEGs and 5 sex-related unigenes, in other similar species, were investigated using qRT-PCR, the results of which were confirmed to data arising from RNA-seq. Among the DEGs, sex-related genes were identified, including forkhead box L2 (Foxl2), sex determining region Y-box (Sox), beta-catenin (β-catenin), chromobox homolog (CBX) and Sex-lethal (Sxl). In addition, 6,283 simple sequence repeats (SSRs) and 614,710 single nucleotide polymorphisms (SNPs) were identified from the RNA-seq results.

CONCLUSIONS

This study provided the first complete gonadal transcriptome data for the blood clam and allowed us to search many aspects of gene sequence information, not limited to gender. This data will improve our understanding of the transcriptomics and reproductive biology of the blood clam. Furthermore, molecular markers such as SSRs and SNPs will be useful in the analysis of genetic evolution, bulked segregant analysis (BSA) and genome-wide association studies (GWAS). Our transcriptome data will therefore provide important genetic information for the breeding and conservation of germplasm.

Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis discus hannai

In order to characterize the female or male transcriptome of the Pacific abalone and further increase genomic resources, we sequenced the mRNA of full-length complementary DNA (cDNA) libraries derived from pooled tissues of female and male Haliotis discus hannai by employing the Iso-Seq protocol of the PacBio RSII platform. We successfully assembled whole full-length cDNA sequences and constructed a transcriptome database that included isoform information. After clustering, a total of 15,110 and 12,145 genes that coded for proteins were identified in female and male abalones, respectively. A total of 13,057 putative orthologs were retained from each transcriptome in abalones. Overall Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyzed in each database showed a similar composition between sexes. In addition, a total of 519 and 391 isoforms were genome-widely identified with at least two isoforms from female and male transcriptome databases. We found that the number of isoforms and their alternatively spliced patterns are variable and sex-dependent. This information represents the first significant contribution to sex-preferential genomic resources of the Pacific abalone. The availability of whole female and male transcriptome database and their isoform information will be useful to improve our understanding of molecular responses and also for the analysis of population dynamics in the Pacific abalone.

De novo transcriptome sequencing and analysis of male, pseudo-male and female yellow perch, Perca flavescens

Transcriptome sequencing could facilitate discovery of sex-biased genes, biological pathways and molecular markers, which could help clarify the molecular mechanism of sex determination and sexual dimorphism, and assist with selective breeding in aquaculture. Yellow perch has unique gonad system and sexual dimorphism and is an alternative model to study mechanism of sex determination, sexual dimorphism and sexual selection. In this study, we performed the de novo assembly of yellow perch gonads and muscle transcriptomes by high throughput Illumina sequencing. A total of 212,180 contigs were obtained, ranging from 127 to 64,876 bp, and N50 of 1,066 bp. The assembly RNA-Seq contigs (≥200bp) were then used for subsequent analyses, including annotation, pathway analysis, and microsatellites discovery. No female- and pseudo-male-biased genes were involved in any pathways while male-biased genes were involved in 29 pathways, and neuroactive ligand receptor interaction and enzyme of trypsin (enzyme code, EC: 3.4.21.4) was highly involved. Pyruvate kinase (enzyme code, EC: 2.7.1.40), which plays important roles in cell proliferation, was highly expressed in muscles. In addition, a total of 183,939 SNPs, 11,286 InDels and 41,479 microsatellites were identified. This study is the first report on transcriptome information in Percids, and provides rich resources for conducting further studies on understanding the molecular basis of sex determinations, sexual dimorphism, and sexual selection in fish, and for population studies and marker-assisted selection in Percids.