Transcriptomic and Physiological Analysis Reveal Melanin Synthesis-Related Genes and Pathways in Pacific Oysters (Crassostrea gigas).

Shell color is one of the shell traits of molluscs, which has been regarded as an economic trait in some bivalves. Pacific oysters (Crassostrea gigas) are important aquaculture shellfish worldwide. In the past decade, several shell color strains of C. gigas were developed through selective breeding, which provides valuable materials for research on the inheritance pattern and regulation mechanisms of shell color. The inheritance patterns of different shell colors in C. gigas have been identified in certain research; however, the regulation mechanism of oyster pigmentation and shell color formation remains unclear. In this study, we performed transcriptomic and physiological analyses using black and white shell oysters to investigate the molecular mechanism of melanin synthesis in C. gigas. Several pigmentation-related pathways, such as cytochrome P450, melanogenesis, tyrosine metabolism, and the cAMP signaling pathway were found. The majority of differentially expressed genes and some signaling molecules from these pathways exhibited a higher level in the black shell oysters than in the white, especially after L-tyrosine feeding, suggesting that those differences may cause a variation of tyrosine metabolism and melanin synthesis. In addition, the in vitro assay using primary cells from mantle tissue showed that L-tyrosine incubation increased cAMP level, gene and protein expression, and melanin content. This study reveals the difference in tyrosine metabolism and melanin synthesis in black and white shell oysters and provides evidence for the potential regulatory mechanism of shell color in oysters.

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.

Heritability of chronic thermal tolerance and genetic correlations with growth traits in the Pacific oyster (Crassostrea gigas).

The increasing seawater temperature during summer months frequently results in severe mortalities in the Pacific oyster Crassostrea gigas around the world, becoming one of the most significant problems challenging the oyster farming industry. In northern China, significant recurrent summer losses of C. gigas have occurred, and its impact on oyster aquaculture has increased in recent years. Selective breeding for improved oyster resistance to high temperature could help to reduce this massive mortality, but the extent of genetic variation underlying this trait is currently unknown. In this study, we constructed 38 full-sib families using the wild C. gigas and estimated the genetic parameters by performing two month-long high-temperature challenge experiments (30 ℃). Experiment 1 was performed in March 2022 followed by experiment 2 in June 2022 (spawning season). In both challenge experiments, there were significant differences in survival among families, suggesting that C. gigas has a different ability to survive under heat stress. Notably, significantly greater mortality was observed for experiment 2, which related to reproductive status and may contribute to additional stress. Thermal tolerance was defined using both binary test survival and time of death traits. Heritability estimates for thermal tolerance were low to moderate (0.16-0.36 for experiment 1 and 0.16-0.33 for experiment 2) using both a Bayesian (MCMCglmm) and a likelihood-based (ASReml-R) approach and estimated heritability of the threshold animal model using ASReml-R (0.16) appeared to be lower compared to MCMCglmm (0.31-0.32). Notably, the genetic and phenotypic correlations for thermal tolerance between two experiments were 0.463 (BS) to 0.491 (TD) and 0.510 (family survival), respectively, which suggested a significant re-ranking of the family breeding values in different time periods. Finally, the genetic and phenotypic correlations were low between growth traits (shell height, shell length, and shell width) and thermal tolerance, suggesting that selection for these traits should be conducted separately. This study reports the first estimation of genetic parameters for chronic thermal tolerance in C. gigas and indicates that this trait is heritable and selective breeding for thermal tolerance is a feasible and promising approach to reduce summer mortality.

Multi-omic insights into the formation and evolution of a novel shell microstructure in oysters.

BACKGROUND: Molluscan shell, composed of a diverse range of architectures and microstructures, is a classic model system to study the relationships between molecular evolution and biomineralized structure formation. The shells of oysters differ from those of other molluscs by possessing a novel microstructure, chalky calcite, which facilitates adaptation to the sessile lifestyle. However, the genetic basis and evolutionary origin of this adaptive innovation remain largely unexplored. RESULTS: We report the first whole-genome assembly and shell proteomes of the Iwagaki oyster Crassostrea nippona. Multi-omic integrative analyses revealed that independently expanded and co-opted tyrosinase, peroxidase, TIMP genes may contribute to the chalky layer formation in oysters. Comparisons with other molluscan shell proteomes imply that von Willebrand factor type A and chitin-binding domains are basic members of molluscan biomineralization toolkit. Genome-wide identification and analyses of these two domains in 19 metazoans enabled us to propose that the well-known Pif may share a common origin in the last common ancestor of Bilateria. Furthermore, Pif and LamG3 genes acquire new genetic function for shell mineralization in bivalves and the chalky calcite formation in oysters likely through a combination of gene duplication and domain reorganization. CONCLUSIONS: The spatial expression of SMP genes in the mantle and molecular evolution of Pif are potentially involved in regulation of the chalky calcite deposition, thereby shaping the high plasticity of the oyster shell to adapt to a sessile lifestyle. This study further highlights neo-functionalization as a crucial mechanism for the diversification of shell mineralization and microstructures in molluscs, which may be applied more widely for studies on the evolution of metazoan biomineralization.

Parallel evolution in Crassostrea oysters along the latitudinal gradient is associated with variation in multiple genes involved in adipogenesis.

Parallel diversification provides a proper framework for studying the role of natural selection in evolution. Yet, empirical studies from ecological 'non-model' species of invertebrates are limited at the whole genome level. Here, we presented a chromosome-scale genome assembly for Crassostrea angulata and investigated the parallel genomic evolution in oysters. Specifically, we used population genomics approaches to compare two southern-northern oyster species pairs (C. angulata-C. gigas and southern-northern C. ariakensis) along the coast of China. The estimated divergence time of C. angulata and C. gigas is earlier than that of southern and northern C. ariakensis, which aligns with the overall elevated genome-wide divergence. However, the southern-northern C. ariakensis FST profile represented more extremely divergent "islands". Combined with recent reciprocal hybridization studies, we proposed that they are currently at an early stage of speciation. These two southern-northern oyster species pairs exhibited significant repeatability in patterns of genome-wide differentiation, especially in genomic regions with extremely high and low divergence. This suggested that divergent and purifying selection has contributed to the genomic parallelism between southern and northern latitudes. Top differentiated genomic regions shared in these two oyster species pairs contained candidate genes enriched for functions in energy metabolism, especially adipogenesis, which are closely related to reproductive behaviours. These genes might be good candidates for further investigation in vivo. In conclusion, our results suggest that similar divergent selection and shared genomic features could predictably transform standing genetic variation within one species pair into differences in another.

Effects of Dietary Cystine and Tyrosine Supplementation on Melanin Synthesis in the Pacific Oyster (Crassostrea gigas).

Melanogenesis is a multistep process to produce melanin for dark pigmentation in skin coloration. Previous studies in vertebrates demonstrated that cystine and tyrosine amino acids are involved in the melanin synthesis. However, very little is known about the melanogenesis in bivalve. In this study, cystine supplementation for 30 days significantly upregulated the expression of CgB-aat1, CgCbs and CgTyr and pheomelanin content in the Pacific oyster Crassostrea gigas. Transmission electron microscope (TEM) results revealed more melanosomes in the connective tissue and melanin granules were secreted in epithelium of mantle. In contrast, tyrosine supplementation had no clear effect on melanogenesis except the gene expression changes of CgB-aat1 and CgCbs. In addition, prolonged supplementation of cystine or tyrosine for 60 days had a negative impact on melanogenesis. Indeed, after 60 days, expression of most of the melanin synthesis-related genes under study was decreased, and melanin content was significantly reduced, indicating that cystine and tyrosine might inhibit production of eumelanin and pheomelanin, respectively. In addition, in vitro analysis using primary cell culture from mantle tissue indicated that incubation with cystine, tyrosine, or B-AAT1 polypeptide, CBS/TYR recombinant proteins induced the increase of CgB-aat1 and CgCbs expression in a dose-dependent manner, suggesting the presence of a regulatory network in response to cystine and tyrosine amino acids intakes in pheomelanin synthesis-related gene expression. Taken together, these data indicate that cystine-CgB-aat1-CgCbs-CgTyr axis is a potential regulator of the pheomelanin biosynthesis pathway, and thus plays an important role in the mantle pigmentation in C. gigas. This work provides a new clue for selective cultivation of oyster strains with specific shell colors in bivalve breeding.

Transcription factor CgPOU3F4-like regulates expression of pheomelanin synthesis related gene CgB-aat1 in the Pacific oyster (Crassostrea gigas).

Previous study has found that b (0, +) -type amino acid transporter 1 (CgB-aat1) plays an essential role on mantle pigmentation in the Pacific oyster Crassostrea gigas. However, the molecular regulation of CgB-aat1 gene expression remains unclear. Herein, three POU domain family members, CgPOU2F1, CgPOU3F4-like and CgPOU4F3-X1 were characterized and they all had POUs and HOX domains, respectively, which were important in transcriptional regulation. CgPOU3F4-like gene expression was the highest in mantle edge. Subsequently, the dual-luciferase reporter result showed that the core regulatory region of CgB-aat1 gene was from -632 to -350 bp of promoter. In transient co-transfection assays, the strongest activity was activated only by CgPOU3F4-like, suggesting CgPOU3F4-like was a valid transcriptional activator of CgB-aat1 gene promoter. And the structural integrity of CgPOU3F4-like was essential for its activation function. In addition, site directed mutagenesis assay was applied to detect three key binding sites between CgPOU3F4-like and core region of CgB-aat1 gene promoter, and this interaction was verified by ChIP test. Furthermore, CgPOU3F4-like knockdown by RNA interference led to obvious decreases in CgB-aat1 and cystathionine beta-synthase (CgCbs) expressions at both mRNA and protein levels. Collectively, these results indicate that CgPOU3F4-like positively regulate CgB-aat1 gene expression and it may be a critical upstream transcriptional regulation factor in pheomelanin synthesis in C. gigas.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis.

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

Histological, elemental, and ultrastructural analysis of melanin in mantle of Pacific oyster (Crassostrea gigas).

Colorful shell of bivalve is mainly because of the biological pigments, of which melanin plays an important role in shell color formation. More and more studies focus on the genes function involved in melanin synthesis, but relatively few studies address the biochemical character and ultrastructure of melanin in bivalve from microscopic perspective. Here, we investigated the histological structure of mantle of Crassostrea gigas with orange shell color. Distribution of melanin in mantle was verified with histochemical staining. In addition, immunofluorescence technique showed that strongly positive signal of CgTYR was specific to the mantle margin, which is consistence with the location of brown granules in H&E staining. The further result of elementary composition of melanin displayed that metal Ca, Fe, and Zn were detected using scanning transmission electron microscope and energy dispersive spectroscopy mapping methods. Next, based on TEM observations, it was speculated that the series of cellular events leading to the formation and release of melanin. Melanocyte in the primary stage showed many mitochondria and rough endoplasmic reticulum as well as an extensive Golgi complex with numerous vesicles intermingled with melanosome. Subsequently, melanosome was expended and their hue gradually intensified, and Golgi complex and mitochondria were still observed in the cytoplasm. Finally, after melanosome was discharged into intercellular spaces, the disintegration of membranes in some cells, and severe cellular vacuolization. These data enrich the understanding of ultrastructural characteristic and formation of melanin in mantle of bivalve and pave the way for further investigating shell coloration at the cellular level.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation.

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

Molecular characterization of Pax7 and its role in melanin synthesis in Crassostrea gigas.

The paired-box 7 (Pax7) is a transcription factor crucial for skin color polymorphism. However, the mechanism underlying the pigmentation associated with Pax7 in mollusks have yet to be elucidated. In this study, the cDNA sequence of Pax7 in the Pacific oyster Crassostrea gigas (CgPax7) was characterized. Phylogenetically, the identity of deduced amino acid sequence was similar to that of other mollusks and contained 463 amino acids, with conserved features of paired domain (PRD), homeobox domain (HD) and octapeptide. Gene expression analysis revealed that CgPax7 was markedly increased at D-shaped larvae stage and ubiquitously expressed in six examined tissues in adult oyster. The result of whole-mount in situ hybridization (WMISH) showed a restricted pattern of CgPax7 expression on margins of shell valves at D-shaped and umbo larvae stages. Additionally, although CgPax7 silencing had no significant effect on CgMitf expression, it significantly inhibited the expressions of CgPax7, CgTyr, CgTyrp1, CgTyrp2 and CgCdk2, genes involved in Tyr-mediated melanin synthesis. Furthermore, CgPax7 knockdown obviously decreased the tyrosinase activity. Less brown-granules at mantle edge was detected by micrographic examination and melanosomes defect was observed by transmission electron microscopy. It was demonstrated that CgPax7 play a key role in melanin synthesis by regulating Tyr-pathway in C. gigas. These findings indicated the potential framework by which mollusks pigmentation.

Identification and expression of cysteine sulfinate decarboxylase, possible regulation of taurine biosynthesis in Crassostrea gigas in response to low salinity.

Taurine has been reported high amounts in marine animals to maintain osmotic balance between osmoformers and sea water. Approximately 80% of the total amino-acid content is taurine in Pacific oyster Crassostrea gigas, an intertidal and euryhaline species. In this study, we cloned the two copies of cysteine sulfinate decarboxylase (CSAD), the key enzyme in taurine biosynthesis pathway, screened in oyster genome data. Sequentially, we compared the expression patterns of CgCSAD1 and CgCSAD2 under low salinity treatment (8‰ and 15‰) using different families from two populations. There was no correlation between the expression of CSAD and the different population. Notably, CgCSAD1 increased significantly in treated groups for 24 h, but CgCSAD2 had no significant differentiation. Moreover, the results of CgCSAD1 interference provided the evidence of the positive correlation between CgCSAD1 expressions and taurine contents. The zinc finger domain showed in multi-alignment results may be the important character of CgCSAD1 as the key enzyme in taurine biosynthesis to regulate taurine pool in response to low salinity. This study provides a new evidence for the important role of taurine in adaptation to low salinity in oyster. In addition, it is a good model to discuss the function and evolution of the duplication in mollusks.