Reproduction-related genes in the pearl oyster genome

Structure and function of vasa gene in gonadal gametogenesis of Pacific abalone

Pacific abalone (Haliotis discus hannai) is a marine gastropod mollusc with significant economic importance in both global fisheries and aquaculture. However, studies exploring the gonadal development and regulatory mechanisms of Haliotis discus hannai are limited. This study aimed to explore whether the vasa gene acted as a molecular marker for germ cells. Initially, the vasa gene was successfully cloned using the cDNA-end rapid amplification technique. The cloned gene had a 2478-bp-long open reading frame and encoded 825 amino acids. Then, a recombinant expression vector was constructed based on the Vasa protein, and an 87-kDa recombinant protein was prepared. Subsequently, a polyclonal antibody was prepared using the purified recombinant protein. The enzyme-linked immunosorbent assay (ELISA) confirmed the titer of the antibody to be ≥512 K. The immunohistochemical analysis revealed that Vasa was widely expressed in oogonia, Stage I oocytes, spermatogonia, and primary spermatocytes. The specific expression of Vasa in the hermaphroditic gonads of abalone was assessed using western blotting to investigate the effects of different photoperiods (12 L:12D, 24 L:0D, 18 L:6D, and 6 L:18D) on the gonadal development of abalone (P < 0.05), with higher expression levels observed in the ovarian proliferative and spermary maturing stages compared with other developmental stages (P < 0.05). Additionally, Vasa exhibited the highest expression in the spermary and ovary under a photoperiod of 18 L:6D (P < 0.05). These data demonstrated the key role of Vasa in developing germ cells in abalone. They shed light upon the molecular mechanism through which the photoperiod influenced Vasa expression and regulated gonadal development in abalone. The findings might provide theoretical references for analyzing the differentiation pattern of abalone germ cells and the genetic improvement and conservation of germplasm resources.

Molecular characterization of Vasa homolog in the pen shell Atrina pectinata: cDNA cloning and expression analysis during gonadal development

Vasa is an ATP-dependent RNA helicase of the DEAD (Asp-Glu-Ala-Asp) box family and a representative component of the germ plasm. In this study, we cloned the full-length vasa homolog in the bivalve Atrina pectinata (psvasa), and performed phylogenetic analysis, mRNA expression analysis for tissue-specific distributions, and immunostaining analysis to reveal its histological localization. The sequence of psvasa was 3587 bp in length and contained a 5' untranslated region of 150 bp, an open reading frame of 2214 bp, and a 3' untranslated region of 1223 bp. The deduced amino acid sequence of psvasa was 737 amino acids long and contained evolutionarily conserved sequences reported in other animals. The mRNA expression analysis showed the highest expression levels in the gonads. Expression was especially high in the ovaries, followed by the testes. The immunostaining analysis showed Vasa-positive cells in the developing gonads, suggesting the presence of putative germ stem cells contributing to the supply of germ cells. Furthermore, characteristic Vasa signals were observed in the basophilic nuclei of the oocytes, suggesting that psvasa plays an important role in the progression of meiosis in oocytes.

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.

Female reproduction-specific proteins, origins in marine species, and their evolution in the animal kingdom

The survival of a species largely depends on the ability of individuals to reproduce, thus perpetuating their life history. The advent of metazoans (i.e. pluricellular animals) brought about the evolution of specialized tissues and organs, which in turn led to the development of complex protein regulatory pathways. This study sought to elucidate the evolutionary relationships between female reproduction-associated proteins by analyzing the transcriptomes of representative species from a selection of marine invertebrate phyla. Our study identified more than 50 reproduction-related genes across a wide evolutionary spectrum, from Porifera to Vertebrata. Among these, a total of 19 sequences had not been previously reported in at least one phylum, particularly in Porifera. Moreover, most of the structural differences between these proteins did not appear to be determined by environmental pressures or reproductive strategies, but largely obeyed a distinguishable evolutionary pattern from sponges to mammals.

Extensive conservation of the proneuropeptide and peptide prohormone complement in mollusks

As one of the most diverse groups of invertebrate animals, mollusks represent powerful models for neurobiological and developmental studies. Neuropeptides and peptide hormones are a heterogeneous class of signalling molecules involved in chemical communication between neurons and in neuroendocrine regulation. Here we present a fine-grained view of the molluscan neuropeptide and peptide hormone toolkit. Our results expand the distribution of several peptide families (e.g., prokineticin, insulin-related peptides, prohormone-4, LFRFamide) within Lophotrochozoa and provide evidence for an early origin of others (e.g., GNXQN/prohormone-2, neuroparsin). We identified a new peptide family broadly distributed among conchiferan mollusks, the PXRX family. We found the Wnt antagonist dickkopf1/2/4 ortholog in lophotrochozoans and nematodes and reveal that the egg-laying hormone family is a DH44 homolog restricted to gastropods. Our data demonstrate that numerous peptides evolved much earlier than previously assumed and that key signalling elements are extensively conserved among extant mollusks.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

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.

Characterisation of Reproduction-Associated Genes and Peptides in the Pest Land Snail, Theba pisana

Increased understanding of the molecular components involved in reproduction may assist in understanding the evolutionary adaptations used by animals, including hermaphrodites, to produce offspring and retain a continuation of their lineage. In this study, we focus on the Mediterranean snail, Theba pisana, a hermaphroditic land snail that has become a highly invasive pest species within agricultural areas throughout the world. Our analysis of T. pisana CNS tissue has revealed gene transcripts encoding molluscan reproduction-associated proteins including APGWamide, gonadotropin-releasing hormone (GnRH) and an egg-laying hormone (ELH). ELH isoform 1 (ELH1) is known to be a potent reproductive peptide hormone involved in ovulation and egg-laying in some aquatic molluscs. Two other non-CNS ELH isoforms were also present in T. pisana (Tpi-ELH2 and Tpi-ELH3) within the snail dart sac and mucous glands. Bioactivity of a synthetic ELH1 on sexually mature T. pisana was confirmed through bioassay, with snails showing ELH1-induced egg-laying behaviours, including soil burrowing and oviposition. In summary, this study presents a detailed molecular analysis of reproductive neuropeptide genes in a land snail and provides a foundation for understanding ELH function.

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.

Bivalve-specific gene expansion in the pearl oyster genome: implications of adaptation to a sessile lifestyle

INTRODUCTION

Bivalve molluscs have flourished in marine environments, and many species constitute important aquatic resources. Recently, whole genome sequences from two bivalves, the pearl oyster, Pinctada fucata, and the Pacific oyster, Crassostrea gigas, have been decoded, making it possible to compare genomic sequences among molluscs, and to explore general and lineage-specific genetic features and trends in bivalves. In order to improve the quality of sequence data for these purposes, we have updated the entire P. fucata genome assembly.

RESULTS

We present a new genome assembly of the pearl oyster, Pinctada fucata (version 2.0). To update the assembly, we conducted additional sequencing, obtaining accumulated sequence data amounting to 193× the P. fucata genome. Sequence redundancy in contigs that was caused by heterozygosity was removed in silico, which significantly improved subsequent scaffolding. Gene model version 2.0 was generated with the aid of manual gene annotations supplied by the P. fucata research community. Comparison of mollusc and other bilaterian genomes shows that gene arrangements of Hox, ParaHox, and Wnt clusters in the P. fucata genome are similar to those of other molluscs. Like the Pacific oyster, P. fucata possesses many genes involved in environmental responses and in immune defense. Phylogenetic analyses of heat shock protein70 and C1q domain-containing protein families indicate that extensive expansion of genes occurred independently in each lineage. Several gene duplication events prior to the split between the pearl oyster and the Pacific oyster are also evident. In addition, a number of tandem duplications of genes that encode shell matrix proteins are also well characterized in the P. fucata genome.

CONCLUSIONS

Both the Pinctada and Crassostrea lineages have expanded specific gene families in a lineage-specific manner. Frequent duplication of genes responsible for shell formation in the P. fucata genome explains the diversity of mollusc shell structures. These duplications reveal dynamic genome evolution to forge the complex physiology that enables bivalves to employ a sessile lifestyle in the intertidal zone.

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.

Neuropeptides encoded by the genomes of the Akoya pearl oyster Pinctata fucata and Pacific oyster Crassostrea gigas: a bioinformatic and peptidomic survey

BACKGROUND

Oysters impart significant socio-ecological benefits from primary production of food supply, to estuarine ecosystems via reduction of water column nutrients, plankton and seston biomass. Little though is known at the molecular level of what genes are responsible for how oysters reproduce, filter nutrients, survive stressful physiological events and form reef communities. Neuropeptides represent a diverse class of chemical messengers, instrumental in orchestrating these complex physiological events in other species.

RESULTS

By a combination of in silico data mining and peptide analysis of ganglia, 74 putative neuropeptide genes were identified from genome and transcriptome databases of the Akoya pearl oyster, Pinctata fucata and the Pacific oyster, Crassostrea gigas, encoding precursors for over 300 predicted bioactive peptide products, including three newly identified neuropeptide precursors PFGx8amide, RxIamide and Wx3Yamide. Our findings also include a gene for the gonadotropin-releasing hormone (GnRH) and two egg-laying hormones (ELH) which were identified from both oysters. Multiple sequence alignments and phylogenetic analysis supports similar global organization of these mature peptides. Computer-based peptide modeling of the molecular tertiary structures of ELH highlights the structural homologies within ELH family, which may facilitate ELH activity leading to the release of gametes.

CONCLUSION

Our analysis demonstrates that oysters possess conserved molluscan neuropeptide domains and overall precursor organization whilst highlighting many previously unrecognized bivalve idiosyncrasies. This genomic analysis provides a solid foundation from which further studies aimed at the functional characterization of these molluscan neuropeptides can be conducted to further stimulate advances in understanding the ecology and cultivation of oysters.

Molecular characterization and analysis of a putative 5-HT receptor involved in reproduction process of the pearl oyster Pinctada fucata

5-HT (5-hydroxytryptamine; serotonin) has been linked to a variety of biological roles including gonad maturation and sequential spawning in bivalve molluscs. To gain a better understanding of the effects of 5-HT on developmental regulation in the pearl oyster Pinctada fucata, the isolation, cloning, and expression of the 5-HT receptor was investigated in this study. A full-length cDNA (2541 bp) encoding a putative 5-HT receptor (5-HTpf) of 471 amino acids was isolated from the ovary of the pearl oyster. It shared 71% and 51% homology, respectively, with the Crassostrea gigas 5-HT receptor and the Aplysia californica 5-HT1ap. The 5-HTpf sequence possessed the typical characteristics of seven transmembrane domains and a long third inner loop. Phylogenetic analysis also indicated that 5-HTpf was classified into the 5-HT1 subtype together with other invertebrate 5-HT1 receptors. Quantitative RT-PCR showed that 5-HTpf is widely expressed in all tissues tested, is involved in the gametogenesis cycle, embryonic and larval development stages, and expression is induced by E2 in ovarian tissues. These results suggest that 5-HTpf is involved in the reproductive process, specifically in the induction of oocyte maturation and spawning of P. fucata.

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.