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

Sex-specific mRNA alternative splicing patterns and Dmrt1 isoforms contribute to sex determination and differentiation of oyster

Alternative splicing (AS) of pre-mRNA is a crucial mechanism that regulates the expression of genes involved in sex determination and differentiation. Despite its importance, AS has been rarely characterized in molluscs. In this study, PacBio Iso-Seq was employed to obtain full-length transcriptome and unveil AS patterns of gonads in the Pacific oyster Crassostrea gigas. A total of 24,783 AS events were identified across 6259 genes, with many enriched in phosphorylation-related processes. Splicing factors were found to drive a high frequency of AS events in gonads. Significant sex-based differences in isoform abundance and the incidence of AS events were observed. Comparative analysis of mature female and male gonads revealed a subset of overlapping differential alternative splicing genes and differentially expressed genes enriched in processes related to microtubule function and cell motility. In addition, the expression levels of sex-biased genes were found correlated with their isoform number in both female and male gonads. A novel isoform of Dmrt1 was identified with male specific expression in mature gonads. This study provides the first comprehensive understanding of full-length transcriptome and AS patterns in molluscan gonads, shedding light on the post-transcriptional regulatory mechanisms underlying sex determination and differentiation in molluscs and potentially across other animals.

Vertebrate endocrine disruptors induce sex-reversal in blue mussels

Mollusks are the second most diverse animal phylum, yet little is known about their endocrinology or how they respond to endocrine disrupting compound (EDC) pollution. Characteristic effects of endocrine disruption are reproductive impairment, skewed sex ratios, development of opposite sex characteristics, and population decline. However, whether classical vertebrate EDCs, such as steroid hormone-like chemicals and inhibitors of steroidogenesis, exert effects on mollusks is controversial. In the blue mussel, Mytilus edulis, EDC exposure is correlated with feminized sex ratios in wild and laboratory mussels, but sex reversal has not been confirmed. Here, we describe a non-destructive qPCR assay to identify the sex of M. edulis allowing identification of males and females prior to experimentation. We exposed male mussels to 17α-ethinylestradiol and female mussels to ketoconazole, EDCs that mimic vertebrate steroid hormones or inhibit their biosynthesis. Both chemicals changed the sex of individual mussels, interfered with gonadal development, and disrupted gene expression of the sex differentiation pathway. Impacts from ketoconazole treatment, including changes in steroid levels, confirmed a role for steroidogenesis and steroid-like hormones in mollusk endocrinology. The present study expands the possibilities for laboratory and field monitoring of mollusk species and provides key insights into endocrine disruption and sexual differentiation in bivalves.

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.

Phenotypic sorting of individual male and female intersex Cherax quadricarinatus and analysis of molecular differences in the gonadal transcriptome

Cherax quadricarinatus exhibit sexual dimorphism, with males outpacing females in size specification and growth rate. However, there is limited understanding of the molecular mechanisms underlying sex determination and sex differentiation in crustaceans. To study the differences between intersex individuals and normal individuals, this study counted the proportion of intersex individuals in the natural population, collected the proportion of 7 different phenotypes in 200 intersex individuals, and observed the differences in tissue sections. RNA-seq was used to study the different changes in the transcriptome of normal and intersex gonads. The results showed that: the percentage of intersex in the natural population was 1.5 %, and the percentage of different types of intersex ranged from 0.5 % to 22.5 %; the sections revealed that the development of normal ovaries was stagnant at the primary oocyte stage when intersex individuals with ovaries were present; We screened for pathways and genes that may be associated with gonadal development and sex, including ovarian steroid synthesis, estrogen signaling pathway, oocyte meiosis, progesterone-mediated oocyte maturation, etc. Relevant genes including tra2a, dmrta2, ccnb2, foxl2, and smad4. This study provides an important molecular basis for sex determination, sex-controlled breeding, and unisex breeding in red crayfish.

Gamma irradiation-induced offspring masculinization is associated with epigenetic changes in female zebrafish

Sex ratio variation is a key topic in ecology, because of its direct effects on population dynamics and thus, on animal conservation strategies. Among factors affecting sex ratio, types of sex determination systems have a central role, since some species could have a sex determined by genetic factors, environmental factors or a mix of those two. Yet, most studies on the factors affecting sex determination have focused on temperature or endocrine-disrupting chemicals (EDCs), and much less is known regarding other factors. Exposure to gamma irradiation was found to trigger offspring masculinization in zebrafish. Here we aimed at deciphering the potential mechanisms involved, by focusing on stress (i.e. cortisol) and epigenetic regulation of key genes involved in sex differentiation in fish. Cortisol levels in exposed and control (F0) zebrafish females' gonads were similar. However, irradiation increased the DNA methylation level of foxl2a and cyp19a1a in females of the F0 and F1 generation, respectively, while no effects were detected in testis. Overall, our results suggest that parental exposure could alter offspring sex ratio, at least in part by inducing methylation changes in ovaries.

Genome-wide identification and expression analysis of Dmrt gene family and their role in gonad development of Pacific oyster (Crassostrea gigas)

Doublesex and Mab-3-related transcription factor (Dmrt) is a type of transcription factor with a zinc-finger DM structural domain, which plays a significant role in sex determination and differentiation in animals. Although Dmrt has been reported in many vertebrates and invertebrates, it has rarely been studied in bivalves. In this study, a total of three members of the Dmrt gene family were identified and characterized in Crassostrea gigas, and all these CgDmrt genes contained a conserved DM domain. Analysis of the phylogenetic tree and gene structure revealed that Dmrt genes clustered on one branch may have similar functions in bivalves. Expression profiling of CgDmrt mRNA in different tissues and stages of gonad development indicated that CgDmrt3 exhibited sexually dimorphic expression and played an important role in the development of the male gonad in C. gigas. Furthermore, analysis of CgDmrt mRNA expression between fertile triploids and sterile triploids showed that CgDmrt3 may be involved in sperm production. Collectively, the systematic analysis of the CgDmrt genes will provide potential insights into the function of these genes in gonadal development.

Transcriptional regulation of IAG by dsx and foxl-2 in mud crab (Scylla paramamosain)

Scylla paramamosain is an important cultured crab species on the southeast coast of China. However, the molecular regulation mechanism of its gonadal development still has not been thoroughly studied. Dsx (doublesex) and foxl-2 (forkhead transcription factor gene 2) are important transcription factors involved in gonadal development. So far, studies on the functions of dsx and foxl-2 in crustaceans are very limited. Insulin-like androgenic gland hormone (IAG) is an effector molecule that regulates the differentiation, development and sex maintenance of testes in crustaceans. In this study, the promoter region of Sp-IAG was predicted, and several potential binding sites of dsx and foxl-2 were found. Site-directed mutagenesis was performed on the predicted potential binding sites, and their promoter activity was analyzed. The results showed that there was a dsx and a foxl-2 binding site, respectively, that could regulate the expression of Sp-IAG. In order to verify the regulatory effect of these two transcription factors on Sp-IAG, we constructed the expression plasmids of dsx and foxl-2 and co-transfected them into HEK293T cell lines with the promoter of Sp-IAG, respectively. The results showed that dsx could significantly promote the expression of Sp-IAG, while foxl-2 could inhibit its expression substantially. Then we carried out in vivo RNA interference experiment on mud crabs. The expression of dsx and foxl-2 in crabs was interfered respectively. The results of qRT-PCR showed that the expression of Sp-IAG was significantly inhibited after interfering with dsx, while significantly increased after interfering with foxl-2, which was consistent with the cell experiment. In conclusion, dsx and foxl-2 transcription factors play opposite roles in regulating the expression of Sp-IAG.

Toward invasive mussel genetic biocontrol: Approaches, challenges, and perspectives

Invasive freshwater mussels, such as the zebra (Dreissena polymorpha), quagga (Dreissena rostriformis bugensis), and golden (Limnoperna fortunei) mussel have spread outside their native ranges throughout many regions of the North American, South American, and European continents in recent decades, damaging infrastructure and the environment. This review describes ongoing efforts by multiple groups to develop genetic biocontrol methods for invasive mussels. First, we provide an overview of genetic biocontrol strategies that have been applied in other invasive or pest species. Next, we summarize physical and chemical methods that are currently in use for invasive mussel control. We then describe the multidisciplinary approaches our groups are employing to develop genetic biocontrol tools for invasive mussels. Finally, we discuss the challenges and limitations of applying genetic biocontrol tools to invasive mussels. Collectively, we aim to openly share information and combine expertise to develop practical tools to enable the management of invasive freshwater mussels.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

Genome-wide identification and expression analysis of Dmrt genes in bivalves

In recent years, some common themes in the development of sex-specific traits in different animal lineages have started to emerge since the discovery of the Dmrt (doublesex-mab3-related transcription factor gene) genes. Bivalves are characterized by a diversity of sexual systems, including simultaneous hermaphroditism, sequential hermaphroditism, and strict gonochorism. However, to date, no research has focused on the genome-wide characterization and analysis of Dmrt genes in bivalves. In this study, the identification and analysis of Dmrt genes in 15 bivalves were performed using bioinformatics methods. A total of 55 Dmrt genes were retrieved in the studied bivalve genomes. The number of Dmrt genes in different species ranged from 3 to 5. The phylogenetic tree showed that Dmrt genes in bivalves can be subdivided into 5 classes: the Dmrt2-like class, Dmrt3-like class, Dmrt4/5-like class, Dsx-like class, and scallop-specific Dmrt class. The Ka/Ks ratios suggested that all Dmrt classes underwent purifying selection pressure. Furthermore, the spatiotemporal expression of Dmrt genes in four bivalve species suggested that different Dmrt genes may have different functions, and scallop-specific Dmrt genes may play a key role in sex determination/differentiation. In general, this study provides a molecular basis for in-depth examination of the functions of Dmrt genes and phylogenomic analyses in bivalves.

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log₁₀ (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log₁₀ (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log₁₀ (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log₁₀ (Foxl2) and Log₁₀ (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

A chromosome-level assembly supports genome-wide investigation of the DMRT gene family in the golden mussel (Limnoperna fortunei)

BACKGROUND

The golden mussel (Limnoperna fortunei) is a highly invasive species that causes environmental and socioeconomic losses in invaded areas. Reference genomes have proven to be a valuable resource for studying the biology of invasive species. While the current golden mussel genome has been useful for identifying new genes, its high fragmentation hinders some applications.

FINDINGS

In this study, we provide the first chromosome-level reference genome for the golden mussel. The genome was built using PacBio HiFi, 10X, and Hi-C sequencing data. The final assembly contains 99.4% of its total length assembled to the 15 chromosomes of the species and a scaffold N50 of 97.05 Mb. A total of 34,862 protein-coding genes were predicted, of which 84.7% were functionally annotated. A significant (6.48%) proportion of the genome was found to be in a hemizygous state. Using the new genome, we have performed a genome-wide characterization of the Doublesex and Mab-3 related transcription factor gene family, which has been proposed as a target for population control strategies in other species.

CONCLUSIONS

From the applied research perspective, a higher-quality genome will support genome editing with the aim of developing biotechnology-based solutions to control invasion. From the basic research perspective, the new genome is a high-quality reference for molecular evolutionary studies of Mytilida and other Lophotrochozoa, and it may be used as a reference for future resequencing studies to assess genomic variation among different golden mussel populations, unveiling potential routes of dispersion and helping to establish better control policies.

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.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

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.

Discovery of Nanos1 and Nanos2/3 as Germ Cell Markers During Scallop Gonadal Development

Nanos are conserved genes involved in germline cell specification and differentiation. However, little is known about the role of different members of Nanos family in germ cell development in mollusks. In the present study, we conducted genome-wide identification of Nanos family in an economically important scallop Patinopecten yessoensis, and detected their expression in adult tissues and during early development. Two Nanos genes (PyNanos1, PyNanos2/3) were identified, both of which have the N-terminal NOT1-interacting motif and C-terminal (CCHC)₂ zinc finger domain. Expression profiles showed that PyNanos1 and PyNanos2/3 were primarily expressed in the gonads, with PyNanos1 being localized in the oogonia, oocytes, and spermatogonia, while PyNanos2/3 being specifically in spermatogonia. The results suggest that PyNanos are germ cell specific and may play crucial roles in gametogenesis in the scallop. PyNanos1 is a maternal gene, which is distributed uniformly at early cleavage, and restricted to 2-3 cell clusters from blastulae to trochophore larvae, suggesting its potential role in the formation of PGCs. Zygotically expressed PyNanos2/3 displayed a similar signal with PyNanos1 in the trochophore larvae, suggesting it may also participate in the formation and/or maintenance of PGCs. This study will benefit germplasm exploitation and conservation in bivalves, and facilitate a better understanding of the evolution of Nanos family and the role of different Nanos in germ cell development in mollusks.

Identification and Characterization of MicroRNAs Involving in Initial Sex Differentiation of Chlamys farreri Gonads

Simple Summary

Sex formation of gonads encompasses two ancient and highly conserved biological processes, sex determination and sex differentiation. The processes are strictly regulated by a complex of gene networks. There is increasing evidence that miRNAs play key roles in many biological processes. however, information is limited in their contribution to sex differentiation in animals. In the present study, we identified the novel miRNAs involved in sex-related genes regulation and explored the miRNA–mRNA networks underlying the posttranscriptional regulation during the initial sex differentiation in Zhikong scallop, Chlamys farreri. Our findings provide an important basis for studying the sex differentiation mechanisms, as well as developing sex control techniques in bivalves.

Abstract

Research on expressional regulation of genes at the initial sex differentiation of gonads will help to elucidate the mechanisms of sex determination and differentiation in animals. However, information on initial sex differentiation of gonads is limited in bivalves. MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that can regulate the target gene expression at the posttranscription level by degrading the mRNA or repressing the mRNA translation. In the present study, we investigated the small RNAs transcriptome using the testes and ovaries of Zhikong scallop Chlamys farreri juveniles with a shell height of 5.0 mm, a critical stage of initial sex differentiation of gonads. A total of 75 known mature miRNAs and 103 novel miRNAs were identified. By comparing the expression of miRNAs between the ovary and testis, 11 miRNAs were determined to be differentially expressed. GO annotations and KEGG analyses indicated that many putative target genes that matched to these differentially expressed miRNAs participated in the regulation of sex differentiation. Furthermore, two selected miRNAs, cfa-novel_miR65 and cfa-miR-87a-3p_1, were confirmed to downregulate expressions of Foxl2 (a female-critical gene) and Klf4 (a male-critical gene), respectively, using a dual-luciferase reporter analysis. Our findings provided new insights into the initial sex differentiation of gonads regulated by miRNAs in 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.

The Effect of Temperature on Gonadal Sex Differentiation of Yesso Scallop Patinopecten yessoensis

Many marine organisms are generally poikilotherms, making seawater temperature one of the most important environmental factors affecting gonadal sex differentiation. Mollusca is the second-largest animal phylum with diverse reproductive systems, but studies on the impact of temperature on sex differentiation are limited to a few sequential hermaphrodites. By combining morphological and molecular analyses, we investigated the effect of temperature on gonadal sex differentiation of a commercially important gonochoristic scallop Patinopecten yessoensis in the field and under laboratory conditions. Based on the relative expression of FoxL2 and Dmrt1L in the gonads of 6- to 12 month-old scallops, we found the scallops start to differentiate at 7 months old in September when the seawater temperature was 21°C. To eliminate the effect of factors other than temperature on sex differentiation, we compared the gonadal development of juvenile scallops at different temperatures (21, 16 and 11°C) under laboratory conditions. After 50 days of treatment, the 11°C group contain more germ cell types, and have higher sex differentiation rates than the 21°C group. But no obvious sex bias was observed. These results suggest that high temperature (21°C) inhibits sex differentiation, whereas low temperature (11°C) accelerates sex differentiation by 2 months for this cold-water species. It also supports juvenile P. yessoensis is gonochoristic rather than protandrous hermaphroditic. Our study addresses for the first time an environmental influence associated with genetic controls on scallop sex differentiation. It will facilitate a better understanding of how environmental factors affect gonadal development in poikilotherms, especially in the less studied molluscs.

Screening and Identification of Transcription Factors Potentially Regulating Foxl2 Expression in Chlamys farreri Ovary

Foxl2 is an evolutionarily conserved female sex gene, which is specifically expressed in the ovary and mainly involved in oogenesis and ovarian function maintenance. However, little is known about the mechanism that regulates Foxl2 specific expression during the ovary development. In the present study, we constructed the gonadal yeast one-hybrid (Y1H) library of Chlamysfarreri with ovaries and testes at different developmental stages using the Gateway technology. The library capacity was more than 1.36 × 10⁷ CFU, and the length of the inserted fragment was 0.75 Kb~2 Kb, which fully met the demand of yeast library screening. The highly transcriptional activity promoter sequence of C. farreri Foxl2 (Cf-Foxl2) was determined at -1000~-616 bp by dual-luciferase reporter (DLR) assay and was used as bait to screen possible transcription factors from the Y1H library. Eleven candidate factors, including five unannotated factors, were selected based on Y1H as well as their expressional differences between ovaries and testes and were verified for the first time to be involved in the transcriptional regulation of Cf-Foxl2 by RT-qPCR and DLR. Our findings provided valuable data for further studying the specific regulation mechanism of Foxl2 in the ovary.

Tissue distribution and seasonal accumulation of carotenoids in Yesso scallop (Mizuhopecten yessoensis) with orange adductor muscle

Carotenoids are colored compounds with important physiological functions. The Haida golden scallop, which has an orange adductor muscle, is a carotenoid-enriched variety of scallop Mizuhopecten yessoensis, an important aquaculture shellfish. In this study, we investigated the tissue distribution of the carotenoids, pectenolone and pectenoxanthin, in both Haida golden scallop and normal Yesso scallop. Both carotenoids were detected in all the sampled tissues of the two scallops, except in the adductor muscle of normal scallop. There were significantly more carotenoids in Haida golden scallop than in normal scallop, in the tissues of the mantle, female gonad, kidney, and adductor muscle. Increased carotenoid concentrations were detected in Haida golden scallop adductor muscle during the spring spawning season, indicating the effects of reproduction on muscle carotenoids accumulation. This study was the first systematic investigation of carotenoid distribution in Yesso scallop tissues and will benefit future research on carotenoid accumulation and function in scallops.

Characterization of the foxl2 gene involved in the vtg expression in mud crab (Scylla paramamosain)

Forkhead box protein L2 (Foxl2) is involved in multiple physiological processes, such as ovarian development, granulosa cell differentiation, ovarian follicle development, and oocyte growth. In this study, a Spfoxl2 gene encoded 530 amino acid protein with characteristic forkhead (FH) domain was identified from transcriptome data of mud crab Scylla paramamosain and validated the accuracy by PCR technology. Meanwhile, the orthologues of the Spfoxl2 gene in other 14 crustacean species were identified with the same method. Further multiple sequence alignment analysis revealed the Foxl2 was highly conserved, especially in the FH domain, even completely identical in several species. Besides, the semi-quantitative PCR (Sq-PCR) result showed Spfoxl2 gene was mainly expressed in the gonad (testis and ovary). Further quantitative real-time PCR (qRT-PCR) result demonstrated its expression level in the testis was significantly higher than that in the ovary (p < 0.01). In addition, the qRT-PCR result showed that in zoea V, megalopa, and larval I, the expression level of Spfoxl2 in megalopa is the highest. In addition, a putative Foxl2 binding site was identified on the promoter region of Spvtg, and knockdown of Spfoxl2 mediated by RNAi technology increased the expression of Spvtg in the ovary, suggesting Spfoxl2 might be the upstream negative regulator of Spvtg. Overall, this study provided new insights into the role of Spfoxl2 in ovary development through regulating Spvtg expression in S. paramamosain.

Identification and functional analysis of a sex-biased transcriptional factor Foxl2 in the bay scallop Argopecten irradians irradians

Transcription factor Foxl2 is an evolutionarily conserved gene playing pivotal roles in regulation of early ovarian differentiation and maintenance in animals. However, the Foxl2 gene has not been thoroughly studied in hermaphroditic scallops. In this study, we cloned and characterized a Foxl2 (designated as AiFoxl2) from the bay scallop Argopecten irradians irradians. The open reading frame of AiFoxl2 was 1122 bp encoding 373 amino acids residues and contained a conserved forkhead box domain. Quantitative real-time PCR showed that AiFoxl2 was mainly expressed in the ovary. Moreover, the highest expression of AiFoxl2 in the ovary was detected at proliferative stage and growing stage, while the lowest level was found at resting stage. During the embryonic and larval development, expression of AiFoxl2 was found first in fertilized eggs, increased significantly at the blastula stage, and reached peak value at the D-larvae stage. When AiFoxl2 was knocked down, testis development-related genes (Dmrt1, Sox7 and Sox9) were up-regulated significantly while the ovary development-related genes (Vg, HSD14, and GATA-1) were down-regulated manifestly. These findings suggested that AiFoxl2 was a female-related gene in A. i. irradians and may be involved in regulation of ovarian development and differentiation.

Transgenerational epigenetic sex determination: Environment experienced by female fish affects offspring sex ratio

Sex determination is a complex process that can be influenced by environment in various taxa. Disturbed environments can affect population sex ratios and thus threaten their viability. Emerging evidences support a role of epigenetic mechanisms, notably DNA methylation, in environmental sex determination (ESD). In this work, using zebrafish as model and a transgenerational experiment comprising 4 successive generations, we report a strength link between the promotor methylation level of three genes in female gonads and population sex ratio. One generation of zebrafish was exposed throughout its lifetime to cadmium (Cd), a non-essential metal, at an environmentally relevant concentration. The subsequent generations were not exposed. At the first and the third generation a subset of individuals was exposed to an elevated temperature, a well-known masculinizing factor in zebrafish. While heat was associated to an increase in the methylation level of cyp19a1a gene and population masculinization, foxl2a/dmrt1 methylation levels appeared to be influenced by Cd and fish density leading to offspring feminization. Ancestral Cd exposure indeed led to a progressive feminization of the population over generations and affected the sex plastic response of zebrafish in response to heat. The effect of Cd on the methylation level of foxl2a was observed until the third generation, supporting potential transgenerational inheritance. Our results support (i) a key role of cyp19a1a methylation in SD in zebrafish in response to environmental cues and (ii) the fact that the environment experienced by parents, namely mothers in the present case, can affect their offspring sex ratio via environment-induced DNA methylation changes in gonads.

Variance in expression and localization of sex-related genes CgDsx, CgBHMG1 and CgFoxl2 during diploid and triploid Pacific oyster Crassostrea gigas gonad differentiation

Several evolutionarily conserved classes of transcriptional regulators were involved in diverse sex determination and differentiation pathways across taxa, whereas their roles in most mollusks is still limited. The Pacific oyster Crassostrea gigas, a dioecious bivalve with sex reversal, could be an ideal model for this issue because of its complex sexuality and potential disruption of sex differentiation in triploid individuals. Here, two mRNA splicing isoforms of a DM domain gene CgDsx and two isoforms of a novel sex-related CgBHMG1 (ortholog of BHMG1 in mammals) were identified in C. gigas. Real time PCR showed that two isoforms of CgDsx and one isoform of CgBHMG1 displayed male-specific expression in diploid oysters, opposite with the female-specific CgFoxl2 (a potential factor of female gonadic differentiation). Interestingly, the four sex-specific transcripts in diploid oyster were expressed in triploid oysters with opposite sex, triploid hermaphrodites and individuals at stage I that sex could not be determined. Subsequent in situ hybridization analysis on gonads of diploid oysters revealed predominant expression of CgDsx in spermatogonia of testes, CgBHMG1 in spermatocytes of testes and follicle cells of ovaries, and CgFoxl2 in follicle cells of ovaries and some male germ cells in testes. And aberrant co-expression of the three genes in triploid oysters was localized in gonadal tubules of gonads at stage I, ovarian follicle cells and undetermined gonial cells in nontypical hermaphroditic gonads with rare female materials. From the above, temporal and spatial expression of sex-related genes in diploid and triploid gonads indicated that CgDsx and CgFoxl2 might mainly function in C. gigas sex differentiation, and CgBHMG1 appeared as a factor involved in meiosis. This work will help to illuminate the gene network of sex differentiation in bivalves and provides new sight on this issue from comparison between diploid and triploid individuals.

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.

Potential GnRH and steroidogenesis pathways in the scallop Patinopecten yessoensis

Gonadotropin-releasing hormone (GnRH) controls synthesis of sex steroid hormones through hypothalamic-pituitary-gonadal (HPG) axis in vertebrates. But in mollusks, research on GnRH and steroidogenesis pathways is still limited. In this study, we first identified two gonadotropin receptor like genes (LGR and LGR5L) and four steroidogenesis-related genes (CYP17A, HSD17B12, HSD3B1 and HSD3B2) in the scallop Patinopecten yessoensis. By examining the expression of 11 genes in the ganglia and/or gonad as well as the concentration of progesterone, testosterone and estradiol in the gonad, we postulate that a potential GnRH signaling pathway (GnRH-GnRHR-GPB5-LGR/LGR5L) in the cerebral and pedal ganglia (CPG) and steroidogenesis pathway (CYP17A, HSD17B12 and HSD3B1) in the gonad are involved in regulating sex steroid hormones. E₂/T index that indicates aromatase activity is higher in the ovary than testis and is positively correlated with the expression of FOXL2 in the gonad, implying the presence of aromatase in the scallop. In addition, we confirmed that expression of most of the downstream genes in the two pathways was significantly elevated after injection of mature py-GnRH peptide. This study would contribute to a new understanding of the molecular basis underlying reproduction regulation by GnRH in mollusks.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Chromosome-level analysis of the Crassostrea hongkongensis genome reveals extensive duplication of immune-related genes in bivalves

Crassostrea hongkongensis is a popular and important native oyster species that is cultured mainly along the coast of the South China Sea. However, the absence of a reference genome has restricted genetic studies and the development of molecular breeding schemes for this species. Here, we combined PacBio and 10 × Genomics technologies to create a C. hongkongensis genome assembly, which has a size of 610 Mb, and is close to that estimated by flow cytometry (~650 Mb). Contig and scaffold N50 are 2.57 and 4.99 Mb, respectively, and BUSCO analysis indicates that 95.8% of metazoan conserved genes are completely represented. Using a high-density linkage map of its closest related species, C. gigas, a total of 521 Mb (85.4%) was anchored to 10 haploid chromosomes. Comparative genomic analyses with other molluscs reveal that several immune- or stress response-related genes extensively expanded in bivalves by tandem duplication, including C1q, Toll-like receptors and Hsp70, which are associated with their adaptation to filter-feeding and sessile lifestyles in shallow sea and/or deep-sea ecosystems. Through transcriptome sequencing, potential genes and pathways related to sex determination and gonad development were identified. The genome and transcriptome of C. hongkongensis provide valuable resources for future molecular studies, genetic improvement and genome-assisted breeding of oysters.

Expression profiling of the Kdm genes in scallop Patinopecten yessoensis suggests involvement of histone demethylation in regulation of early development and gametogenesis

Histone demethylation modification is an important means of gene expression regulation and is widely involved in biological processes such as animal reproduction and development. Histone lysine demethylases (Kdm) plays an important role in the demethylation of histones. To understand the role of histone demethylation in scallops, we identified the Kdm gene family of the Yesso scallop Patinopecten yessoensis, and analyzed its expression during the gonad and early development. The results showed that the P. yessoensis has a complete Kdm family including seventeen members that belong to sixteen subfamilies (Hif1an, Hspbap1, Jarid2, Jmjd4, Jmjd6, Jmjd7, Jmjd8, Kdm1, Kdm2, Kdm3, Kdm4, Kdm5, Kdm6, Kdm7, Kdm8 and Kdm9). The Kdm genes showed five different expression patterns in the early development of scallop, with some of them (e.g. Jmjd7, Jmjd8 and Kdm8) being highly expressed in only one or two stage and the others (e.g. Kdm1A, Kdm9, Jmjd4 and Jmjd6) in several consecutive stages. During gonadal development, the Kdm genes also display various expression patterns. Some genes (e.g. Kdm2, Kdm4 and Jmjd7) display preferential expression in the testis, and the others have no obvious sex bias but show stage preference (resting, proliferative, growing or maturation stage). These results suggest that various histone demethylation modifications (e.g. H3K4, H3K9 and H3K27) may participate in the regulation of gametogenesis and early development of Yesso scallop. It will facilitate a better understanding of the epigenetic contributions to mollusk development.

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.

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.

Dynamics of DNA Methylation and DNMT Expression During Gametogenesis and Early Development of Scallop Patinopecten yessoensis

DNA methylation reprograms during gametogenesis and embryo development, which is essential for germ cell specification and genomic imprinting in mammals. Corresponding process remains poorly investigated in molluscs. Here, we examined global DNA methylation level in the gonads of scallop Patinopecten yessoensis during gametogenesis and in embryos/larvae at different stages. DNA methylation level fluctuates during gametogenesis and early development, peaking at proliferative stage of ovary, growing stage of testis, and in blastulae. To understand the mechanisms underlying these changes, we conducted genome-wide characterization of DNMT family and investigated their expression profiles based on transcriptomes and in situ hybridization. Three genes were identified, namely PyDNMT1, PyDNMT2, and PyDNMT3. Expression of PyDnmt3 agrees with DNA methylation level during oogenesis and early development, suggesting PyDNMT3 may participate in de novo DNA methylation that occurs mainly at proliferative stage of ovary and testis, and in blastulae and gastrulae. PyDnmt1 expression is positively correlated with DNA methylation level during spermatogenesis, and is higher at maturation stage of ovary and in 2-8 cell embryos than other stages, implying possible involvement of PyDNMT1 in DNA methylation maintenance during meiosis and embryonic development. This study will facilitate better understanding of the developmental epigenetic reprogramming in bivalve molluscs.