Cg-SoxE and Cg-β-catenin, two new potential actors of the sex-determining pathway in a hermaphrodite lophotrochozoan, the Pacific oyster Crassostrea gigas

Polystyrene microplastics photo-aged under simulated sunlight influences gonadal development in the Pacific oyster

Microplastics (MPs) aging in natural ecosystems are caused by solar irradiation. Photo-aged MPs in aquatic systems are a major threat to molluscs. In this study, polystyrene (PS) photo-aging was simulated using a sunlight simulator. After exposure of Crassostrea gigas to photo-aged PS, a decreased gonadosomatic index, coupled with histological alterations, suggested an inhibitory effect on the gonadal development of bivalves. As the concentration of aged PS increased, the inhibitory effects on gonadal development became more severe. The sex hormone (testosterone and estradiol) and energy metabolism (glycogen, lipid, and protein content) differences between C. gigas males and females suggested a disruption of sex hormonal homeostasis and a shift in energy allocation strategy, which may have affected reproduction, especially female oysters. In addition, the substantial downregulation of SOX-8, SOX-E, Piwi1, and TGF-β genes may be contributing factors causing the inhibitory effect of aged PS on the gonadal development of C. gigas. This study provides an essential reference for evaluating the reproductive health risks posed by aged MPs and offers novel insights and perspectives for exploring the impact of MPs under natural conditions.

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.

Cytogenetic characterization of the golden mussel (Limnoperna fortunei) reveals the absence of sex heteromorphic chromosomes

The golden mussel (Limnoperna fortunei) is an aggressive invasive species in South America, where it endangers native species and freshwater ecosystems, in addition to causing extensive economic losses, mainly to the hydroelectric sector. Currently, there's no efficient control method available and the invasion has progressed across the continent. Its high reproduction rate is one of the key factors of the golden mussel's high invasive potential and, recently, efforts have been done in order to understand the reproduction and the sexual features of this species. However, its cytogenetics characterization is incipient and the possible occurrence of sex-specific cytogenetic features was never investigated. In this study, we aimed to characterize the chromosomal morphometry, the distribution profile of heterochromatin, and to detect possible sex-related epigenetic marks in the golden mussel. Results revealed that the karyotypic structure is similar in both sexes and no chromosome heteromorphism was observed between males and females specimens. The data increment the cytogenetic characterization of Limnoperna fortunei and contribute for future studies that aim to further investigate its reproduction and underlying sex determination processes.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

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.

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.

Integrated mRNA and miRNA expression profile analysis of female and male gonads in Hyriopsis cumingii

Hyriopsis cumingii is an important species for freshwater pearl cultivation in China. In terms of pearl production, males have larger pearls and better glossiness than females, but there are few reports focusing on the sex of H. cumingii. In this study, six mRNA and six microRNA (miRNA) libraries were prepared from ovaries and testes. Additionally, 28,502 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Compared with testis, 14,360 mRNAs and 20 miRNAs were up-regulated in ovary, 14,142 mRNAs and 12 miRNAs were down-regulated. In DEGs, the known genes related to sex determinism and/or differentiation were also identified, such as DMRT1, SOX9, SF1 for males, FOXL2 for females, and other potentially significant candidate genes. Three sex-related pathways have also been identified, which are Wnt, Notch, and TGF-beta. In 32 DEMs, the three miRNAs (miR-9-5p, miR-92, miR-184) were paid more attention, they predicted 28 target genes, which may also be candidates for sex-related miRNAs and genes. Differential miRNAs target genes analysis reveals the pathway associated with oocyte meiosis and spermatogenesis. Overall, the findings of the study provide significant insights to enhance our understanding of sex differentiation and/or sex determination mechanisms for H. cumingii.

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.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Identification of Hc-β-catenin in freshwater mussel Hyriopsis cumingii and its involvement in innate immunity and sex determination

β-catenin is a multifunctional protein that participates in a variety of physiological activities, including immune regulation, sex determination, nervous system development and, cell differentiation. However, the function of β-catenin in freshwater mussel Hyriopsis cumingii remains unclear. Herein, the gene encoding β-catenin from H. cumingii (Hc-β-catenin) was cloned and characterised. The full-length 5544 bp gene includes an open reading frame (ORF) of 2463 bp encoding a putative protein of 820 amino acids residues containing 12 armadillo (ARM) repeats. After injecting H. cumingii with Aeromonas hydrophila or lipopolysaccharides, Hc-β-catenin transcription was induced in hemocytes and gills, and the greatest responses occurred at 24 h after bacterial challenge, confirming an important role in immune responses. Quantitative real-time PCR analysis showed that Hc-β-catenin mRNA was distributed in the gill, foot, liver, kidney, mantle, adductor muscle and gonad of male and female mussels. In gonad, Hc-β-catenin expression was markedly higher in females than males. During the embryonic period, Hc-β-catenin expression was highest at 3 day. In 1-, 2- and 3-year-old mature mussels, Hc-β-catenin expression in female gonad tissue was notably higher than in males. In situ hybridisation revealed a significant hybridisation signal in female gonads, indicating that Hc-β-catenin is a pro-ovarian, anti-testis gene. Our findings demonstrate that Hc-β-catenin is important in immune regulation and sex determination in freshwater mussel.

SOHLHs Might Be Gametogenesis-Specific bHLH Transcriptional Regulation Factors in Crassostrea gigas

The self-renewal and differentiation of germ cells are essential for gametogenesis and reproduction. In mammals, the transcription factors SOHLH1 and SOHLH2, two members of the bHLH family, are specifically expressed in the gonads, and play an important role in spermatocyte and oocyte differentiation. In our previous study, we performed a phylogenetic analysis of the Lophotrochozoa bHLH genes, and two Sohlh were identified in the Pacific oyster Crassostrea gigas. Based on the genomes of other species that have complete genomic information, we further analyzed the phylogenetics of the Sohlh in this study. The results indicate that the Sohlh are ancient genes that were lost in many species during evolution, including in some invertebrates, and lower vertebrates. The phylogenetic tree shows that Sohlh1 and Sohlh2 are located in different scaffolds and that they have low similarity, suggesting early separation in invertebrates. We used RNA-seq and RT-PCR to examine the mRNA expression of the Sohlh in C. gigas (termed Cg-Sohlh), we found that Cg-Sohlh1, and Cg-Sohlh2 are specifically expressed in the gonads. During gonadal development, the mRNA expression levels of both genes increased from the proliferative stage and reached the highest level at the growth stage (P < 0.05). Then, the expression level decreased until the resting stage. In addition, immunohistochemistry was used to determine that the Cg-SOHLH1 protein was specifically expressed in the spermatogonia and spermatocytes. Cg-Sohlh2 mRNA was expressed in both the male and female gonads, while Cg-Sohlh1 mRNA was highly expressed in the female gonads at all developmental stages except for the resting stage. These data indicate that Cg-SOHLH might be gonad-specific regulatory factors, similar to mammalian SOHLH, and that Cg-SOHLH1 might be involved in spermatogonial differentiation. This study lays the foundation to further determine the functional role of SOHLH in mollusk gametogenesis and provides a foundation to better understand the regulatory mechanism of gametogenesis in invertebrates.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

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

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

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

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

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

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

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

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.

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.

Analyses of Sox-B and Sox-E Family Genes in the Cephalopod Sepia officinalis: Revealing the Conserved and the Unusual

Cephalopods provide an unprecedented opportunity for comparative studies of the developmental genetics of organ systems that are convergent with analogous vertebrate structures. The Sox-family of transcription factors is an important class of DNA-binding proteins that are known to be involved in many aspects of differentiation, but have been largely unstudied in lophotrochozoan systems. Using a degenerate primer strategy we have isolated coding sequence for three members of the Sox family of transcription factors from a cephalopod mollusk, the European cuttlefish Sepia officinalis: Sof-SoxE, Sof-SoxB1, and Sof-SoxB2. Analyses of their expression patterns during organogenesis reveals distinct spatial and temporal expression domains. Sof-SoxB1 shows early ectodermal expression throughout the developing epithelium, which is gradually restricted to presumptive sensory epithelia. Expression within the nervous system appears by mid-embryogenesis. Sof-SoxB2 expression is similar to Sof-SoxB1 within the developing epithelia in early embryogenesis, however appears in largely non-overlapping expression domains within the central nervous system and is not expressed in the maturing sensory epithelium. In contrast, Sof-SoxE is expressed throughout the presumptive mesodermal territories at the onset of organogenesis. As development proceeds, Sof-SoxE expression is elevated throughout the developing peripheral circulatory system. This expression disappears as the circulatory system matures, but expression is maintained within undifferentiated connective tissues throughout the animal, and appears within the nervous system near the end of embryogenesis. SoxB proteins are widely known for their role in neural specification in numerous phylogenetic lineages. Our data suggests that Sof-SoxB genes play similar roles in cephalopods. In contrast, Sof-SoxE appears to be involved in the early stages of vasculogenesis of the cephalopod closed circulatory system, a novel role for a member of this gene family.

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.

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.