Molecular cloning and gene expression of Cg-Foxl2 during the development and the adult gametogenetic cycle in the oyster Crassostrea gigas

Expression and functional analysis of Fushi Tarazu transcription factor 1 (FTZ-F1) in the regulation of steroid hormones during the gonad development of Fujian Oyster, Crassostrea angulata

Crassostrea angulata, a major shellfish cultivated in Southern China, has experienced a notable surge in commercial value in recent years. Understanding the molecular mechanisms governing their reproductive processes holds significant implications for advancing aquaculture practices. In this study, we cloned the orphan nuclear receptor gene, Fushi Tarazu transcription factor 1 (FTZ-F1), of C. angulata and investigated its functional role in the gonadal development. The full-length cDNA of FTZ-F1 spans 2357 bp and encodes a protein sequence of 530 amino acids. Notably, the amino acid sequence of FTZ-F1 in C. angulata shares remarkable similarity with its homologues in other species, particularly in the DNA-binding region (>90%) and ligand-binding region (>44%). In C. angulata, the highest expression level of FTZ-F1 was observed in the ovary, exhibiting more than a 200-fold increase during the maturation stage compared to the initiation stage (P < 0.001). Specifically, FTZ-F1 was mainly expressed in the follicular cells surrounding the oocytes of C. angulata. Upon inhibiting FTZ-F1 gene expression in C. angulata through RNA interference (RNAi), a substantial reduction in the expression of genes involved in the synthesis of sex steroids in the gonads, including 3β-HSD, Cyp17, and follistatin, was observed. In addition, estradiol (E2) and testosterone (T) levels also showed a decrease upon FTZ-F1 silencing, resulting in a delayed gonadal development. These results indicate that FTZ-F1 acts as a steroidogenic factor, participating in the synthesis and regulation of steroid hormones and thus playing an important role in the reproductive and endocrine systems within oysters.

Characterization and expression patterns of the Fox gene family under heat and cold stress in Manila clam Ruditapes philippinarum based on genome-wide identification

In this study, the Fox gene family of Ruditapes philippinarum was identified by bioinformatics analysis and genome data. The results showed that a total of 21 Fox genes were identified in R. philippinarum, which were divided into 16 subfamilies, including two members of Foxa subfamily (Foxa1, Foxa2), three members of Foxl subfamily (Foxl1b, Foxl1a, FOXL2), three members of Foxn subfamily (FOXN3, FOX4A, Foxn4b) and one member of other families. The chromosome distribution, domains, conserved motifs, introns, exons and protein tertiary structures of these 21 Fox genes were predicted. By analyzing the RNA-seq data of R. philippinarum, it was found that the Fox gene family was differentially expressed in different tissues, different developmental stages and under heat and cold stress. Most of Fox genes were highly expressed in four tissues: labial palp, gonad, gill and foot. Most of the Fox genes were highly expressed in blastula stage. Most of the Fox genes were highly expressed in high temperature group of two populations, and Foxo, FOXG1 were highly expressed in low temperature group. In addition, qPCR showed that the expression levels of Foxo and Foxj1b genes increased significantly under acute cold stress. Therefore, we speculate that Fox genes may play important roles in embryo development and the temperature stress of R. philippinarum, and this study provides a basis for further exploring the molecular mechanism of low temperature tolerance mediated by Fox.

Transcriptomic and metabolomic analyses reveal sex-related differences in the gonads of Pinctada fucata martensii

Pinctada fucata martensii is an economically important bivalve mollusk, as this species makes a major contribution to seawater pearl production. Pearl production efficiency varies between the sexes of P. f. martensii, but many aspects of the molecular mechanisms underlying sex determination and sex differentiation in P. f. martensii remain unclear. Here, transcriptomic and metabonomic analyses were conducted to identify the major genes and metabolic changes associated with sex determination and gametogenesis. We identified a total of 3426 differentially expressed genes (DEGs) between females and males. These included Fem-1c and Foxl2, which are involved in sex determination and sex differentiation, and SOHLH2, Nanos1 and TSSK4, which are involved in gametogenesis. We also identified a total of 5231 significant differential metabolites (SDMs) between females and males. These DEGs were enriched in 47 metabolic pathways, including "ABC transporters," "purine metabolism," and "glycerophospholipid metabolism." Our findings provide new insights into the molecular mechanisms underlying sex determination, sex differentiation, and gametogenesis and will aid future studies of P. f. martensii.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Genome-wide identification of STATs and analysis of their role in sex determination in Pacific oysters (Crassostrea gigas)

STAT (signal transducer and activator of the transcription) proteins, are a group of highly conserved transcription factors and fundamental components of the JAK-STAT signaling pathway. They play crucial roles in a variety of biological processes, such as immunity, proliferation, differentiation, and growth. However, little information is known regarding their role in gonad development and sex determination in mollusks. In this study, we identified 3 STAT genes in Pacific Oyster Crassostrea gigas. Phylogenetic analysis showed that STATs from mollusks were highly conserved, and most of them had four identical motif regions, except for the STAT1 and STAT3 predicted sequences from Crassostrea hongkongensis. Tissue expression analysis indicated CgSTAT1 had a high expression level in most tissues, while CgSTAT3 had a low expression level in most tissues. Expression analysis of early developmental stages showed CgSTAT1 had a higher expression level from egg to D shaped larva and a lower expression level in subsequent stages. In contrast CgSTAT1, CgSTAT2 had a reverse expression pattern. Expression analysis of different developmental stages of diploid gonads indicated that CgSTAT1 had a higher expression level at the S1 and S3 stages relative to the S2 stage in females, while in males the S3 stage had a higher expression than than the S2 stage. The expression level of CgSTAT1 between diploids and triploids in females differed significantly, but there were no significant differences in males. Expression of CgSTAT2 differed significantly between diploid and triploid males. These data suggest an important role for STATs in sex differentiation in diploid and triploid oysters. Our study is the first to explore the role of STATs in sex differentiation and gonadal development in oysters, and will help us better understand the molecular mechanisms of sex differentiation in shellfish.

Effects of Benzo[a]pyrene on Food Metabolism and Reproductive Endocrine and Ovarian Development in Female Scallop Chlamys farreri at Different Reproductive Stages

Benzo[a]pyrene (B[a]P), a polycyclic aromatic hydrocarbon (PAH) with the most carcinogenic effects of all the PAHs, has multiple toxic effects on marine bivalves. We investigated the interference mechanism of B[a]P on food metabolism (sugars, proteins, and sugars), and on reproductive endocrine and ovarian development in female scallops (Chlamys farreri). Scallops were exposed to different concentrations of B[a]P concentrations of 0, 0.38, 3.8, and 38 μg/L throughout gonadal development. Total cholesterol and triglyceride contents in the digestive glands were increased, and their synthesis genes were upregulated. The plasma glucose contents decreased with the inhibition of glycogen synthesis genes and the induction of glycolysis genes in the digestive gland. The results showed that B[a]P had endocrine-disrupting effects on scallops, that it negatively affected genes related to ovarian cell proliferation, sex differentiation, and egg development, and that it caused damage to ovarian tissue. Our findings supplement the information on B[a]P disruption in gonadal development of marine bivalves. Environ Toxicol Chem 2024;00:1-14. © 2023 SETAC.

Starvation-induced changes in sex ratio involve alterations in sex-related gene expression and methylation in Pacific oyster Crassostrea gigas

Aquatic animals are subject to varying degrees of starvation stress in their natural habitats due to food limitations. Consequently, starvation is a crucial environmental factor for sex determination in many species; however, limited research has been conducted on the effects of starvation on sex determination in shellfish. Here, four full sibling families of Pacific oysters were established and subjected to starvation stress. The results demonstrated that starvation caused the sex ratio (female to male) to change from 1:0.78 to 1:1.44 and resulted in a delay in gonadal development. Further studies revealed that the expression levels of DNA methylation-related genes Dnmt1 (DNA methyltransferase 1), Dnmt3a/b (DNA methyltransferase 3a/b) and Tet3 (tet methylcytosine dioxygenase 3) decreased under starvation stress. Conversely, the upregulation of Dmrt1 (doublesex and mab-3 related transcription factor 1), a gene typically associated with males, in females suggests that the increased proportion of males may be linked to starvation-induced high expression of this particular gene. In addition, the gene Dgkd (diacylglycerol kinase delta), which is involved in the regulation of second messenger protein kinase C, was differentially methylated between males and females, with the methylation level of this gene gradually increasing with male development, while the methylation level of this gene decreased under starvation stress. Correlation analysis of Dgkd methylation levels with expression levels showed a negative correlation between DNA methylation and gene expression. Finally dual fluorescence reporter experiments confirmed that DNA methylation suppressed Dgkd expression in vitro. In summary, the results suggest that starvation may regulate Dgkd gene expression through DNA methylation and thus affect Dmrt1 expression, thereby causing sex reversal in the Pacific oyster. The outcomes resolved how environmental factors are involved in sex reversal from an epigenetic perspective and provided a theoretical basis for sex control breeding in the Pacific oyster.

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.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

Identification and functional analysis of foxl2 and nodal in sea cucumber, Apostichopus japonicus

Sea cucumber (Apostichopus japonicus) is an important mariculture species in China. To date, the mechanisms of sex determination and differentiation in sea cucumber remain unclear. Identifying sex-specific molecular markers is an effective method for revealing the genetic basis of sex determination and sex differentiation. In this study, foxl2 and nodal homologous genes were identified in A. japonicus. Foxl2 exhibited dynamic and sexually dimorphic expression patterns in the gonads, with prominent expression in the ovaries and minimal expression in the testis according to real-time quantitative PCR (RT-qPCR) study. As nodal was specifically expressed in the ovary, it could serve as an ovary-specific marker in sea cucumber. Additionally, knockdown of foxl2 or nodal using RNA interference (RNAi) led to the down-regulation of piwi, germ cell-less, and dmrt1, suggesting that foxl2 and nodal may play important roles in gonad maintenance of sea cucumber. Overall, this study adds to our understanding of the roles of foxl2 and nodal in the gonadal development of A. japonicus, which provides further insight into the mechanisms of sea cucumber sex determination and differentiation.

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.

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.

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.

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.

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.

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.

A natural antisense transcript of the fem-1 gene was found expressed in female gonads during the characterization, expression profile, and cellular localization of the fem-1 gene in Pacific white shrimp Penaeus vannamei

The fem-1 gene in Caenorhabditis elegans is involved in sex differentiation; it is specifically required for all aspects of male development. In this study, the full-length cDNA of the fem-1 (Pvfem-1) gene was isolated from the Pacific whiteleg shrimp Penaeus vannamei. The Pvfem-1 transcript is 3778 nt long and encodes a putative protein (PvFEM-1) of 638 amino acids that presented eight ankyrin repeats. The translated protein showed a significant (P < 0.05) structural similitude by superposition with C. elegans FEM-1 protein. Pvfem-1 expression was evaluated by qPCR and in situ hybridization (ISH) during embryogenesis, larval development, and gonads of both genders in subadult and adult life stages. Pvfem-1 was found expressed in brain, intestine, hepatopancreas, and in the gonads of both genders in subadults and adults when quantified by RT-qPCR. A significant finding was the discovery of a natural antisense transcript (NAT) of Pvfem-1 by ISH. It was present in the oocyte nucleus of subadult female shrimp gonads but was not seen within oocytes from adult females, although it was detected in follicular cells, suggesting a possible post-transcriptional regulation of Pvfem-1 in female gonad. Conversely, in males, no NAT was observed, and Pvfem-1 was found expressed in spermatogonia of both, subadult and adult shrimps indicating a function in male sexual differentiation and gametes generation. This study represents the first step for future functional analysis that is expected to contribute to clarifying the role of Pvfem-1 in sex differentiation and determination.

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.

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.

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.

Effects of a parental exposure to diuron on Pacific oyster spat methylome

Environmental epigenetic is an emerging field that studies the cause-effect relationship between environmental factors and heritable trait via an alteration in epigenetic marks. This field has received much attentions since the impact of environmental factors on different epigenetic marks have been shown to be associated with a broad range of phenotypic disorders in natural ecosystems. Chemical pollutants have been shown to affect immediate epigenetic information carriers of several aquatic species but the heritability of the chromatin marks and the consequences for long term adaptation remain open questions. In this work, we investigated the impact of the diuron herbicide on the DNA methylation pattern of spat from exposed Crassotrea gigas genitors. This oyster is one of the most important mollusk species produced worldwide and a key coastal economic resource in France. The whole genome bisulfite sequencing (WGBS, BS-Seq) was applied to obtain a methylome at single nucleotide resolution on DNA extracted from spat issued from diuron exposed genitors comparatively to control spat. We showed that the parental diuron exposure has an impact on the DNA methylation pattern of its progeny. Most of the differentially methylated regions occurred within coding sequences and we showed that this change in methylation level correlates with RNA level only in a very small group of genes. Although the DNA methylation profile is variable between individuals, we showed conserved DNA methylation patterns in response to parental diuron exposure. This relevant result opens perspectives for the setting of new markers based on epimutations as early indicators of marine pollutions.

Transcriptome Sequencing and Comparative Analysis of Ovary and Testis Identifies Potential Key Sex-Related Genes and Pathways in Scallop Patinopecten yessoensis

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

Foxl2 and Its Relatives Are Evolutionary Conserved Players in Gonadal Sex Differentiation

Foxl2 is a member of the large family of Forkhead Box (Fox) domain transcription factors. It emerged during the last 15 years as a key player in ovarian differentiation and oogenesis in vertebrates and especially mammals. This review focuses on Foxl2 genes in light of recent findings on their evolution, expression, and implication in sex differentiation in animals in general. Homologs of Foxl2 and its paralog Foxl3 are found in all metazoans, but their gene evolution is complex, with multiple gains and losses following successive whole genome duplication events in vertebrates. This review aims to decipher the evolutionary forces that drove Foxl2/3 gene specialization through sub- and neo-functionalization during evolution. Expression data in metazoans suggests that Foxl2/3 progressively acquired a role in both somatic and germ cell gonad differentiation and that a certain degree of sub-functionalization occurred after its duplication in vertebrates. This generated a scenario where Foxl2 is predominantly expressed in ovarian somatic cells and Foxl3 in male germ cells. To support this hypothesis, we provide original results showing that in the pea aphid (insects) foxl2/3 is predominantly expressed in sexual females and showing that in bovine ovaries FOXL2 is specifically expressed in granulosa cells. Overall, current results suggest that Foxl2 and Foxl3 are evolutionarily conserved players involved in somatic and germinal differentiation of gonadal sex.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

Transcriptomics Analysis of Crassostrea hongkongensis for the Discovery of Reproduction-Related Genes

BACKGROUND

The reproductive mechanisms of mollusk species have been interesting targets in biological research because of the diverse reproductive strategies observed in this phylum. These species have also been studied for the development of fishery technologies in molluscan aquaculture. Although the molecular mechanisms underlying the reproductive process have been well studied in animal models, the relevant information from mollusks remains limited, particularly in species of great commercial interest. Crassostrea hongkongensis is the dominant oyster species that is distributed along the coast of the South China Sea and little genomic information on this species is available. Currently, high-throughput sequencing techniques have been widely used for investigating the basis of physiological processes and facilitating the establishment of adequate genetic selection programs.

RESULTS

The C.hongkongensis transcriptome included a total of 1,595,855 reads, which were generated by 454 sequencing and were assembled into 41,472 contigs using de novo methods. Contigs were clustered into 33,920 isotigs and further grouped into 22,829 isogroups. Approximately 77.6% of the isogroups were successfully annotated by the Nr database. More than 1,910 genes were identified as being related to reproduction. Some key genes involved in germline development, sex determination and differentiation were identified for the first time in C.hongkongensis (nanos, piwi, ATRX, FoxL2, β-catenin, etc.). Gene expression analysis indicated that vasa, nanos, piwi, ATRX, FoxL2, β-catenin and SRD5A1 were highly or specifically expressed in C.hongkongensis gonads. Additionally, 94,056 single nucleotide polymorphisms (SNPs) and 1,699 simple sequence repeats (SSRs) were compiled.

CONCLUSIONS

Our study significantly increased C.hongkongensis genomic information based on transcriptomics analysis. The group of reproduction-related genes identified in the present study constitutes a new tool for research on bivalve reproduction processes. The large group of molecular markers discovered in this study will be useful for population screening and marker assisted selection programs in C.hongkongensis aquaculture.

Molecular signatures discriminating the male and the female sexual pathways in the pearl oyster Pinctada margaritifera

The genomics of economically important marine bivalves is studied to provide better understanding of the molecular mechanisms underlying their different reproductive strategies. The recently available gonad transcriptome of the black-lip pearl oyster Pinctada margaritifera is a novel and powerful resource to study these mechanisms in marine mollusks displaying hermaphroditic features. In this study, RNAseq quantification data of the P. margaritifera gonad transcriptome were analyzed to identify candidate genes in histologically-characterized gonad samples to provide molecular signatures of the female and male sexual pathway in this pearl oyster. Based on the RNAseq data set, stringent expression analysis identified 1,937 contigs that were differentially expressed between the gonad histological categories. From the hierarchical clustering analysis, a new reproduction model is proposed, based on a dual histo-molecular analytical approach. Nine candidate genes were identified as markers of the sexual pathway: 7 for the female pathway and 2 for the male one. Their mRNA levels were assayed by real-time PCR on a new set of gonadic samples. A clustering method revealed four principal expression patterns based on the relative gene expression ratio. A multivariate regression tree realized on these new samples and validated on the previously analyzed RNAseq samples showed that the sexual pathway of P. margaritifera can be predicted by a 3-gene-pair expression ratio model of 4 different genes: pmarg-43476, pmarg-foxl2, pmarg-54338 and pmarg-fem1-like. This 3-gene-pair expression ratio model strongly suggests only the implication of pmarg-foxl2 and pmarg-fem1-like in the sex inversion of P. margaritifera. This work provides the first histo-molecular model of P. margaritifera reproduction and a gene expression signature of its sexual pathway discriminating the male and female pathways. These represent useful tools for understanding and studying sex inversion, sex differentiation and sex determinism in this species and other related species for aquaculture purposes such as genetic selection programs.

Transcriptomic profiling of gametogenesis in triploid Pacific Oysters Crassostrea gigas: towards an understanding of partial sterility associated with triploidy

BACKGROUND

Triploidy can occur in many animal species but is often lethal. Among invertebrates, amphibians and fishes, triploids are viable although often sterile or infertile. Most triploids of the Pacific oyster Crassostrea gigas are almost sterile (named "3nβ") yet a low but significant proportion show an advanced gametogenesis (named "3nα"). These oysters thus constitute an interesting model to study the effect of triploidy on germ cell development. We used microarrays to compare the gonad transcriptomes of diploid 2n and the abovementioned triploid 3nβ and 3nα male and female oysters throughout gametogenesis.

RESULTS

All triploids displayed an upregulation of genes related to DNA repair and apoptosis and a downregulation of genes associated with cell division. The comparison of 3nα and 3nβ transcriptomes with 2n revealed the likely involvement of a cell cycle checkpoint during mitosis in the successful but delayed development of gonads in 3nα individuals. In contrast, a disruption of sex differentiation mechanisms may explain the sterility of 3nβ individuals with 3nβ females expressing male-specific genes and 3nβ males expressing female-specific genes.

CONCLUSIONS

The disruption of sex differentiation and mitosis may be responsible for the impaired gametogenesis of triploid Pacific oysters. The function of the numerous candidate genes identified in our study should now be studied in detail in order to elucidate their role in sex determination, mitosis/meiosis control, pachytene cell cycle checkpoint, and the control of DNA repair/apoptosis.

Genomic analysis of the Pacific oyster (Crassostrea gigas) reveals possible conservation of vertebrate sex determination in a mollusc

Despite the prevalence of sex in animal kingdom, we have only limited understanding of how sex is determined and evolved in many taxa. The mollusc Pacific oyster Crassostrea gigas exhibits complex modes of sexual reproduction that consists of protandric dioecy, sex change, and occasional hermaphroditism. This complex system is controlled by both environmental and genetic factors through unknown molecular mechanisms. In this study, we investigated genes related to sex-determining pathways in C. gigas through transcriptome sequencing and analysis of female and male gonads. Our analysis identified or confirmed novel homologs in the oyster of key sex-determining genes (SoxH or Sry-like and FoxL2) that were thought to be vertebrate-specific. Their expression profile in C. gigas is consistent with conserved roles in sex determination, under a proposed model where a novel testis-determining CgSoxH may serve as a primary regulator, directly or indirectly interacting with a testis-promoting CgDsx and an ovary-promoting CgFoxL2. Our findings plus previous results suggest that key vertebrate sex-determining genes such as Sry and FoxL2 may not be inventions of vertebrates. The presence of such genes in a mollusc with expression profiles consistent with expected roles in sex determination suggest that sex determination may be deeply conserved in animals, despite rapid evolution of the regulatory pathways that in C. gigas may involve both genetic and environmental factors.

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

BACKGROUND

Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis.

RESULTS

The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera.

CONCLUSIONS

The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

Vitellogenin of Fujian oyster, Crassostrea angulata: Synthesized in the ovary and controlled by estradiol-17β

In this study, we cloned a full-length cDNA encoding vitellogenin (Vg) in the Fujian oyster Crassostrea angulata. The complete Vg cDNA consists of 5160 nucleotides with a long open reading frame encoding 1641 amino acid residues. The deduced amino acid sequence shared high similarity with the Vgs of other mollusc, fish, nematode and arthropod species, particularly in the N-terminal region. We analyzed the spatiotemporal expression of caVg transcripts by Real-time Quantitative PCR. In common with other mollusc Vgs, the caVg gene was expressed primarily in the ovary, and the levels were 348 and 177 times higher in maturation and ripeness stages (P<0.01), respectively, than in the partially spent stage. There was negligible expression in male oysters. In situ hybridization analysis further localized caVg mRNA to the follicle cells (also named auxiliary cells) surrounding the oocytes in the ovary. Moreover, in vivo waterborne exposure experiments in early gametogenesis oysters showed that estradiol-17β (E2) administration resulted in a significant increase in caVg mRNA expression. We conclude that caVg is synthesized in the follicle cell surrounding the vitellogenic oocyte in C. angulata, and directly passed to oocytes through the extracellular space without mediation through hemolymph. Also, we hypothesize that this process is mediated by E2 in a dose dependent.

Reproduction-related genes in the pearl oyster genome

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

Cloning and expression analysis of Foxl2 during the reproductive cycle in Korean rockfish, Sebastes schlegeli

Foxl2 is a member of the winged helix/forkhead family of transcription factors and is known to regulate ovarian aromatase, which plays a crucial role in ovarian differentiation. To address the role of Foxl2 in gonads and brain during gonadal development, we isolated the full-length cDNA of Foxl2 and analyzed its spatiotemporal expression patterns in the viviparous teleost Korean rockfish, Sebastes schlegeli. Tissue distribution pattern revealed that the Foxl2 was detected in the liver, fat, gill, brain, and ovary, but could hardly be found in the testis. Reverse transcriptase PCR suggested that Foxl2 in Korean rockfish may involve in ovary development in the study of expression level during gonads development. It also revealed that the stage of highest expression level for Foxl2 was almost much earlier than cyp19a1a and cyp19a1b during the gonadal development stage in gonads and brain except for cyp19a1a in brain. Furthermore, the expression pattern of Foxl2 as well as aromatases may imply the role of Foxl2 in the up-regulation of aromatases not only in the female fish but also in male.

Cloning, expression and identification of two glutathione S-transferase isoenzymes from Perna viridis

Glutathione S-transferases (GSTs; EC 2.5.1.18) are phase II enzymes involved in major detoxification reactions of xenobiotic in many organisms. In the present study, two classes of GSTs (PvGST1 and PvGST2) were cloned from P. viridis by rapid amplification of cDNA ends method. Sequence alignments and phylogenetic analysis together supported that PvGST1 and PvGST2 belonged to the pi and omega classes, respectively. The PvGST1 cDNA was 1214 nucleotides (nt) in length and contained a 618 nt open reading frame (ORF) encoding 206 amino acid residues, and had 46 nt of 5'-untranslated region (UTR) and a 3' UTR of 550 nt including a tailing signal (AATAAA) and a poly (A) tail. The molecular mass of the predicted PvGST1 was 23.815kDa, with the calculated isoelectric point being 5.39. PvGST2 was 1093bp, consisting of a 5' UTR of 13bp, a 3' UTR of 246bp and an ORF of 834bp. The deduced protein was composed of 278 amino acids, with an estimated molecular mass of 32.476kDa and isoelectric point of 8.88. Tissue distribution analysis of the PvGST1 and PvGST2 mRNA revealed that the GST expression level was higher in digestive gland and gonad, while lower in gill and mantle in both genders. Molecular modeling analysis of two GSTs implicated their various functions account for their different enzymatic features.

Oyster sex determination is influenced by temperature - first clues in spat during first gonadic differentiation and gametogenesis

The sex-determining system of Crassostrea gigas is still poorly known, especially regarding the potential influence of temperature. In order to address this question, mRNA expressions of actors of the molecular cascade (Cg-DMl, Cg-SoxE, Cg-β-catenin, Cg-Foxl2/Cg-Foxl2os) and of Oyvlg, a germ cell marker, were investigated by real-time PCR in spat grown at different temperatures (18, 22, 25 and 28°C). In parallel, gonadic differentiation, gametogenesis and sex ratios were assessed by histology at each of these temperatures. Whatever the temperature, Cg-DMl, Cg-SoxE, Cg-β-catenin and Oyvlg expressions peaked at the same developmental stage, always after Cg-Foxl2/Cg-Foxl2os (around 40-44dpf for spat grown at 18°C). Temperatures increased the kinetics of first gonadic differentiation and gametogenesis. At 25°C a significant switch occurred in sex ratio towards males and in the balance of expression between male and female genes, in favor of males. A slight gametogenesis disturbance was also observed. These results strengthen the hypotheses about the sex-determining time window and molecular cascade governing the development of C. gigas, with notably the involvement of Cg-Foxl2/Cg-Foxl2os in the very early steps. They also suggest an influence of temperature on the oyster's sex determination which, associated to genetic control, would induce a mixed sex determination system (GSD+TSD).

Cloning and characterization of the dopamine like receptor in the oyster Crassostrea angulata: expression during the ovarian cycle

We cloned and characterized a complete cDNA encoding a dopamine receptor (DAR) named Ca-DA1R from Fujian oyster, Crassostrea angulata. The 2843 bp long cDNA sequence includes a 916-bp 5'-UTR, the 1197 bp ORF which encodes a putative protein of 399 amino acids, and a 729 bp 3'-UTR. The Ca-DA1R sequence possesses typical characteristics of a D1 receptor: two main features being a short third intracellular loop and a long inner COOH-terminal tail domain. Using a real-time PCR approach, expression profiles of Ca-DA1R were analyzed in adult tissues and during the four stages of ovarian development. Ca-DA1R was expressed ubiquitously, although transcript levels varied between tissues, with higher mRNA levels detected in the ovary, labial palps and mantle. During the four stages of ovarian development, Ca-DA1R mRNA expression level was higher in the proliferation stage than in the other three stages during the ovary cycle. In situ hybridization results reveal that the Ca-DA1R mRNA is mainly expressed in the epithelium of the gonoducts. These observations suggest that Ca-DA1R binding of DA probably plays an important role in early ovarian development and via regulating oocyte locomotion cooperates with the 5-HT receptor system during the ovarian cycle in C. angulata.

Cloning and characterization of the follistatin gene from Crassostrea angulata and its expression during the reproductive cycle

Follistatin is an activin-binding protein that prevents activin from binding to its receptor and neutralizes its activity. Follistatin plays a key role in regulating folliculogenesis and the development of ovary. However, limited information on follistatin genes from molluscs is available until now. By using Race, real-time PCR, in situ hybridization and in silico analysis, a full-length cDNA of follistatin of the Portuguese oyster Crassostrea angulata was acquired. The full-length (1297 bp) cDNA of Ca-follistatin encodes a peptide of 241 amino acids. The similarity of its deduced amino acid sequence to these of other invertebrate species was about 60%. Ca-follistatin mRNA transcript was most abundantly expressed in ovary (p<0.05), and it was also expressed in testis, adductor muscle, mantle, gill and visceral mass. In situ hybridization revealed that the expression and distribution of Ca-follistatin gene were expressed exclusively in granulosa cells, neither in cumulus oophorus nor in oocytes. During the reproductive cycle of female oyster (initiation stage, maturation stage, ripeness stage and partially spent stage), the expression of Ca-follistatin in the ovary continuously increased from initiation to ripeness stages attaining its highest value (p<0.05), then the expression level decreased sharply to the lowest point in the partially spent stage (p<0.05), whereas the Ca-follistatin mRNA transcript of male oyster in the testis maintained a relatively stable low level during the first three stages, and also noticeably decreased thereafter (p<0.05). These findings suggest that follistatin is likely to play an important role in the ovary development of oysters by autocrine signaling.

Molecular cloning of Foxl2 gene and the effects of endocrine-disrupting chemicals on its mRNA level in rare minnow, Gobiocypris rarus

Endocrine-disrupting chemicals (EDCs) can affect normal sexual differentiation in fish. Foxl2, one forkhead transcription factor, plays an important role in ovarian differentiation in the early development of the female gonad in mammals and fish. How EDCs affect Foxl2 expression is little known. In this study, we isolated a Foxl2 cDNA from the ovary of rare minnow Gobiocypris rarus and examined its expression during early development stages and in different adult tissues. Then, we analyzed Foxl2 expression in G. rarus juvenile following 3-day exposure to 17α- ethinylestradiol (EE2), 4-n-nonylphenol (NP), and bisphenol A (BPA). Alignment of known Foxl2 sequences among vertebrates showed high identity in forkhead domain and C-terminal region with other vertebrate proteins. Quantitative RT-PCR analysis showed that Foxl2 expression was linear decrease and cyp19a1a, the downstream target gene of Foxl2, had no correlation with Foxl2 from 18 to 50 days post fertilization (dpf). Among different adult tissues, Foxl2 is mainly expressed in ovary, brain, gill, eye, and male spleen. In the 3-day exposure, the juvenile fish to EDCs, 0.1 nM EE2, and 1 nM BPA significantly up-regulated the expression of Foxl2 gene, while NP had no effect on Foxl2 expression. Altogether, these results provide basic data for further study on how Foxl2 mediates EDCs impact on the sexual differentiation in G. rarus.

Gametogenesis in the Pacific oyster Crassostrea gigas: a microarrays-based analysis identifies sex and stage specific genes

Background

The Pacific oyster Crassostrea gigas (Mollusca, Lophotrochozoa) is an alternative and irregular protandrous hermaphrodite: most individuals mature first as males and then change sex several times. Little is known about genetic and phenotypic basis of sex differentiation in oysters, and little more about the molecular pathways regulating reproduction. We have recently developed and validated a microarray containing 31,918 oligomers (Dheilly et al., 2011) representing the oyster transcriptome. The application of this microarray to the study of mollusk gametogenesis should provide a better understanding of the key factors involved in sex differentiation and the regulation of oyster reproduction.

Methodology/Principal Findings

Gene expression was studied in gonads of oysters cultured over a yearly reproductive cycle. Principal component analysis and hierarchical clustering showed a significant divergence in gene expression patterns of males and females coinciding with the start of gonial mitosis. ANOVA analysis of the data revealed 2,482 genes differentially expressed during the course of males and/or females gametogenesis. The expression of 434 genes could be localized in either germ cells or somatic cells of the gonad by comparing the transcriptome of female gonads to the transcriptome of stripped oocytes and somatic tissues. Analysis of the annotated genes revealed conserved molecular mechanisms between mollusks and mammals: genes involved in chromatin condensation, DNA replication and repair, mitosis and meiosis regulation, transcription, translation and apoptosis were expressed in both male and female gonads. Most interestingly, early expressed male-specific genes included bindin and a dpy-30 homolog and female-specific genes included foxL2, nanos homolog 3, a pancreatic lipase related protein, cd63 and vitellogenin. Further functional analyses are now required in order to investigate their role in sex differentiation in oysters.

Conclusions/Significance

This study allowed us to identify potential markers of early sex differentiation in the oyster C. gigas, an alternative hermaphrodite mollusk. We also provided new highly valuable information on genes specifically expressed by mature spermatozoids and mature oocytes.