Comparative Analysis of the Biochemical Composition, Amino Acid, and Fatty Acid Contents of Diploid, Triploid, and Tetraploid Crassostrea gigas.

Tetraploid oysters are artificially produced oysters that do not exist in nature. The successful breeding of 100% triploid oysters resolved the difficulties of traditional drug-induced triploids, such as the presence of drug residues and a low triploid induction rate. However, little is known concerning the biochemical composition and nutrient contents of such tetraploids. Therefore, we investigated compositional differences among diploid, triploid, and tetraploid Crassostrea gigas as well as between males and females of diploids and tetraploids. The findings indicated that glycogen, EPA, ∑PUFA, and omega-3 contents were significantly higher in triploid oysters than in diploids or tetraploids; tetraploid oysters had a significantly higher protein content, C14:0, essential amino acid, and flavor-presenting amino acid contents than diploids or triploids. For both diploid and tetraploids, females had significantly higher levels of glutamate, methionine, and phenylalanine than males but lower levels of glycine and alanine. In addition, female oysters had significantly more EPA, DHA, omega-3, and total fatty acids, a result that may be due to the fact that gonadal development in male oysters requires more energy to sustain growth, consumes greater amounts of nutrients, and accumulates more proteins. With these results, important information is provided on the production of C. gigas, as well as on the basis and backing for the genetic breeding of oysters.

Comparative transcriptome analysis reveals molecular damage associated with cryopreservation in Crassostrea angulata D-larvae rather than to cryoprotectant exposure.

BACKGROUND: The Portuguese oyster Crassostrea angulata, a bivalve of significant economic and ecological importance, has faced a decline in both production and natural populations due to pathologies, climate change, and anthropogenic factors. To safeguard its genetic diversity and improve reproductive management, cryopreservation emerges as a valuable strategy. However, the cryopreservation methodologies lead to some damage in structures and functions of the cells and tissues that can affect post-thaw quality. Transcriptomics may help to understand the molecular consequences related to cryopreservation steps and therefore to identify different freezability biomarkers. This study investigates the molecular damage induced by cryopreservation in C. angulata D-larvae, focusing on two critical steps: exposure to cryoprotectant solution and the freezing/thawing process. RESULTS: Expression analysis revealed 3 differentially expressed genes between larvae exposed to cryoprotectant solution and fresh larvae and 611 differentially expressed genes in cryopreserved larvae against fresh larvae. The most significantly enriched gene ontology terms were "carbohydrate metabolic process", "integral component of membrane" and "chitin binding" for biological processes, cellular components and molecular functions, respectively. Kyoto Encyclopedia of Genes and Genomes enrichment analysis identified the "neuroactive ligand receptor interaction", "endocytosis" and "spliceosome" as the most enriched pathways. RNA sequencing results were validate by quantitative RT-PCR, once both techniques presented the same gene expression tendency and a group of 11 genes were considered important molecular biomarkers to be used in further studies for the evaluation of cryodamage. CONCLUSIONS: The current work provided valuable insights into the molecular repercussions of cryopreservation on D-larvae of Crassostrea angulata, revealing that the freezing process had a more pronounced impact on larval quality compared to any potential cryoprotectant-induced toxicity. Additionally, was identify 11 genes serving as biomarkers of freezability for D-larvae quality assessment. This research contributes to the development of more effective cryopreservation protocols and detection methods for cryodamage in this species.

Temperature effects on plasmalogen profile and quality characteristics in Pacific oyster (Crassostrea gigas) during depuration.

The quality of Pacific oyster (Crassostrea gigas) can be affected by many factors during depuration, in which temperature is the major element. In this study, we aim to determine the quality and plasmalogen changes in C. gigas depurated at different temperatures. The quality was significantly affected by temperature, represented by varying survival rate, glycogen content, total antioxidant capacity, alkaline phosphatase activity between control and stressed groups. Targeted MS analysis demonstrated that plasmalogen profile was significantly changed during depuration with PUFA-containing plasmalogen species being most affected by temperature. Proteomics analysis and gene expression assay further verified that plasmalogen metabolism is regulated by temperature, specifically, the plasmalogen synthesis enzyme EPT1 was significantly downregulated by high temperature and four plasmalogen-related genes (GPDH, PEDS, Pex11, and PLD1) were transcriptionally regulated. The positive correlations between the plasmalogen level and quality characteristics suggested plasmalogen could be regarded as a quality indicator of oysters during depuration.

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

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

CgADAR1 involved in regulating the synthesis of interferon-like protein in Crassostrea gigas.

Adenosine deaminases acting on RNA 1 (ADAR1) is a dsRNA adenosine (A)-to-inosine (I) editing enzyme that regulates the innate immune response against virus invasion. In the present study, a novel CgADAR1 was identified from the oyster Crassostrea gigas. The open reading frame (ORF) of CgADAR1 was of 3444 bp encoding a peptide of 1147 amino acid residues with two Zα domains, one dsRNA binding motif (DSRM) and one RNA adenosine deaminase domain (ADEAMc). The mRNA transcripts of CgADAR1 were detected in all the examined tissues, with higher expression levels in mantle and gill, which were 7.11-fold and 4.90-fold (p < 0.05) of that in labial palp, respectively. The mRNA transcripts of CgADAR1 in haemocytes were significantly induced at 24 h and 36 h after Poly (A: U) stimulation, which were 6.03-fold (p < 0.01) and 1.37-fold (p < 0.001) of that in control group, respectively. At 48 h after Poly (A:U) stimulation, the mRNA expression of CgRIG-Ⅰ, CgIRF8 and CgIFNLP significantly increased, which were 4.36-fold (p < 0.001), 1.82-fold (p < 0.05) and 1.92-fold (p < 0.05) of that in control group. After CgADAR1 expression was inhibited by RNA interference (RNAi), the mRNA expression levels of CgMDA5, CgRIG-Ⅰ, CgTBK1, CgIRF8 and CgIFNLP were significantly increased, which were 11.88-fold, 11.51-fold, 2.22-fold, 2.85-fold and 2.52-fold of that in control group (p < 0.001), and the phosphorylation level of CgTBK1 was also significantly increased. These results suggested that CgADAR1 played a regulation role in the early stages of viral infection by inhibiting the synthesis of interferon-like protein.

In silico structural features of the CgNR5A: CgDAX complex and its role in regulating gene expression of CYP target genes in Crassostrea gigas.

Contamination of aquatic environments has been steadily increasing due to human activities. The Pacific oyster Crassostrea gigas has been used as a key species in studies assessing the impacts of contaminants on human health and the aquatic biome. In this context, cytochrome P450 (CYPs) play a crucial role in xenobiotic metabolism. In vertebrates many of these CYPs are regulated by nuclear receptors (NRs) and little is known about the NRs role in C. gigas. Particularly, the CgNR5A represents a homologue of SF1 and LRH-1 found in vertebrates. Members of this group can regulate genes of CYPs involved in lipid/steroid metabolism, with their activity regulated by other NR, called as DAX-1, generating a NR complex on DNA response elements (REs). As C. gigas does not exhibit steroid biosynthesis pathways, CgNR5A may play other physiological roles. To clarify this issue, we conducted an in silico investigation of the interaction between CgNR5A and DNA to identify potential C. gigas CYP target genes. Using molecular docking and dynamics simulations of the CgNR5A on DNA molecules, we identified a monomeric interaction with extended REs. This RE was found in the promoter region of 30 CYP genes and also the NR CgDAX. When the upstream regulatory region was analyzed, CYP2C39, CYP3A11, CYP4C21, CYP7A1, CYP17A1, and CYP27C1 were mapped as the main genes regulated by CgNR5A. These identified CYPs belong to families known for their involvement in xenobiotic and lipid/steroid metabolism. Furthermore, we reconstructed a trimeric complex, previously proposed for vertebrates, with CgNR5A:CgDAX and subjected it to molecular dynamics simulations analysis. Heterotrimeric complex remained stable during the simulations, suggesting that CgDAX may modulate CgNR5A transcriptional activity. This study provides insights into the potential physiological processes involving these NRs in the regulation of CYPs associated with xenobiotic and steroid/lipid metabolism.

Morphological and Molecular Analysis Identified a Subspecies of Crassostrea ariakensis (Fujita, 1913) along the Coast of Asia.

Crassostrea ariakensis (Fujita, 1913) is one of the most important economic and ecological oysters that is naturally distributed along the coast of Asia, separated by the Yangtze River estuary. They are usually compared as different populations, while there is no consensus on whether C. ariakensis in northern and southern areas should be considered as two species or subspecies. Here, we analyzed morphological characteristics, COI, 16s rRNA, mitogenome sequences, and species delimitation analysis (ASAP and PTP) to resolve the intraspecific taxonomic status of the C. ariakensis. Phylogenetic and ASAP analysis highlight that C. ariakensis was divided into N-type and S-type. PTP was unable to differentiate between the two types of C. ariakensis. The divergence time of N-type and S-type C. ariakinsis is estimated to be 1.6 Mya, using the relaxed uncorrelated lognormal clock method. Additionally, significant morphological differences exist between the two groups in terms of the adductor muscle scar color. Despite these differences, the COI (0.6%) and 16S rRNA (0.6%) genetic distance differences between N-type and S-type C. ariakensis has not yet reached the interspecific level. These results suggest that N-type and S-type C. ariakensis should be treated as different subspecies and renamed as C. ariakensis ariakensis subsp. nov and C. ariakensis meridioyangtzensis subsp. nov.

The c.503A>G polymorphism in ZIP1-II of Pacific oyster Crassostrea gigas associated with zinc content.

The Pacific oyster Crassostrea gigas is renowned for its high zinc content, but the significant variation among individuals diminishes its value as a reliable source of zinc supplementation. The Zrt/Irt-like protein 1 (ZIP1), a pivotal zinc transporter that facilitates zinc uptake in various organisms, plays crucial roles in regulating zinc content. In the present study, polymorphisms of a ZIP1 gene in C. gigas (CgZIP1-II) were investigated, and their association with zinc content was evaluated through preliminary association analysis in 41 oysters and verification analysis in another 200 oysters. A total of 17 single nucleotide polymorphisms (SNPs) were identified in the exonic region of CgZIP1-II gene, with c.503A>G significantly associated with zinc content. Protein sequence and structure prediction showed that c.503A>G caused a p.Met110Val nonsynonymous mutation located in the metal-binding region of CgZIP1-II, which could influence its affinity for zinc ions, thereby modulating its zinc transport functionality. These results indicate the potential influence of CgZIP1-II polymorphisms on zinc content and provide candidate markers for selecting C. gigas with high zinc content.

CgSHIP2 negatively regulates the mRNA expressions of CgIL-17s in response to Vibrio splendidus stimulation in Crassostrea gigas.

SH2 domain containing inositol polyphosphate5-phosphatase-2 (SHIP2) is a member of the 5-phosphatase family, acting as a vital negative regulator of immune response in vertebrates. In the present study, a SHIP2 homologue (designed as CgSHIP2) was identified from Pacific oyster, Crassostrea gigas. There was a SH2 domain, an IPPc domain and a SAM domain in CgSHIP2. The mRNA transcripts of CgSHIP2 were widely expressed in all the tested tissues with the highest expression in haemolymph. The mRNA expressions of CgSHIP2 in haemocytes increased significantly at 6, 12, 48 and 72 h after Vibrio splendidus stimulation. The positive green signals of CgSHIP2 protein were mainly located in cytoplasm of haemocytes. After the expression of CgSHIP2 was inhibited by RNA interference, the mRNA transcripts of interleukin 17s (CgIL-17-1, CgIL-17-2, CgIL-17-3 and CgIL-17-6) in the haemocytes increased significantly at 24 h after V. splendidus stimulation, which were 8.15-fold (p < 0.001), 3.44-fold (p < 0.05), 2.15-fold (p < 0.01) and 4.63-fold (p < 0.05) compared with that in NC-RNAi group, respectively. Obvious branchial swelling and cilium shedding in gills were observed in CgSHIP2-RNAi group at 24 h after V. splendidus stimulation. The results suggested that CgSHIP2 played an important role in controlling inflammatory response induced by bacteria in oysters.

miRNA-seq provides novel insight into the response to hyper- and hypo- salinity acclimation in Crassostrea hongkongensis.

The Hong Kong oyster, Crassostrea hongkongensis, is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis.

Genetic divergence accessed with microsatellite markers reflects the time of Crassostrea gigas genetic breeding in Brazil.

The Pacific Oyster was introduced on Santa Catarina Island in 1987, experiencing processes of selection and genetic breeding since then. Such procedures may have led to the establishment of specific strains, given the saltier and warmer conditions of the Atlantic Ocean. This study employed microsatellite markers to compare allelic patterns of oysters cultivated in Santa Catarina, the USA, and Asia. Specific allelic patterns were revealed in the Santa Catarina samples, reflecting the time of selection/breeding of the oyster in this region. This result supports the effectiveness of the selection/breeding procedures and the demand for protection of this commercially important genetic resource.

Transcriptomic investigation and biomarker discovery for zinc response in oysters Crassostrea gasar.

In an era of unprecedented industrial and agricultural growth, metal contamination in marine environments is a pressing concern. Sentinel organisms such as the mangrove oyster Crassostrea gasar provide valuable insights into these environments' health. However, a comprehensive understanding of the molecular mechanisms underlying their response to metal exposure remains elusive. To address this gap, we reanalyzed the 454-sequencing data of C. gasar, utilizing an array of bioinformatics workflow of CDTA (Combined De Novo Transcriptome Assembly) to generate a more representative assembly. In parallel, C. gasar individuals were exposed to two concentrations of zinc (850 and 4500 µg L-1 Zn) for 48 h to understand their molecular responses. We utilized Trinotate workflow for the 11,684-CDTA unigenes annotation, with most transcripts aligning with the genus Crassostrea. Our analysis indicated that 67.3% of transcript sequences showed homology with Pfam, while 51.4% and 54.5%, respectively had GO and KO terms annotated. We identified potential metal pollution biomarkers, focusing on metal-related genes, such as those related to the GSH biosynthesis (CHAC1 and GCLC-like), to zinc transporters (ZNT2-like), and metallothionein (MT-like). The evolutionary conservation of these genes within the Crassostrea genus was assessed through phylogenetic analysis. Further, these genes were evaluated by qPCR in the laboratory exposed oysters. All target genes exhibited significant upregulation upon exposure to Zn at both 850 and 4500 µg L-1, except for GCLC-like, which showed upregulation only at the higher concentration of 4500 µg L-1. This result suggests distinct activation thresholds and complex interactions among these genes in response to varying Zn concentrations. Our study provides insights into the molecular responses of C. gasar to Zn, adding valuable tools for monitoring metal pollution in marine ecosystems using the mangrove oyster as a sentinel organism.

Exploration of cell-cell interactions and the notch signaling pathway in the gonadal niche of Crassostrea gigas.

The Notch signaling pathway plays a pivotal role in governing cell fate determinations within the gonadal niche. This study provides an extensive elucidation of the male and female gonadal niches within Crassostrea gigas. Examination via transmission electron microscopy revealed the presence of desmosome-like connection not only between germ cells and niche cells but also among adjacent niche cells within the oyster gonad. Transcriptomic analysis identified several putative Notch pathway components, including CgJAG1, CgNOTCH1, CgSuh, and CgHey1. Phylogenetic analysis indicated a close evolutionary relationship between CgJAG1, CgNOTCH1, and CgHey1 and Notch members present in Drosophila. Expression profiling results indicated a notable abundance of CgHey1 in the gonads, while CgJAG1 and CgNOTCH1 displayed distinct expression patterns associated with sexual dimorphism. In situ hybridization findings corroborated the predominant expression of CgJAG1 in male niche cells, while CgNOTCH1 was expressed in both male and female germ cells, as well as female niche cells. These findings demonstrate the important role of the Notch signaling pathway in the gonadal niche of oysters.

Molecular characterization of transcription factor CREB3L2 and CREB3L3 and their role in melanogenesis in Pacific oysters (Crassostrea gigas).

Colorful shells in mollusks are commonly attributable to the presence of biological pigments. In Pacific oysters, the inheritance patterns of several shell colors have been investigated, but little is known about the molecular mechanisms of melanogenesis and pigmentation. cAMP-response element binding proteins (CREB) are important transcription factors in the cAMP-mediated melanogenesis pathway. In this study, we characterized two CREB genes (CREB3L2 and CREB3L3) from Pacific oysters. Both of them contained a conserved DNA-binding and dimerization domain (a basic-leucine zipper domain). CREB3L2 and CREB3L3 were expressed highly in the mantle tissues and exhibited higher expression levels in the black-shell oyster than in the white. Masson-Fontana melanin staining and immunofluorescence analysis showed that the location of CREB3L2 protein was generally consistent with the distribution of melanin in oyster edge mantle. Dual-luciferase reporter assays revealed that CREB3L2 and CREB3L3 could activate the microphthalmia-associated transcription factor (MITF) promoter and this process was regulated by the level of cAMP. Additionally, we found that cAMP regulated melanogenic gene expression through the CREB-MITF-TYR axis. These results implied that CREB3L2 and CREB3L3 play important roles in melanin synthesis and pigmentation in Pacific oysters.

Polymorphisms in the cysteine dioxygenase gene and their association with taurine content in the Pacific oyster Crassostrea gigas.

The Pacific oyster Crassostrea gigas is rich in taurine, which is crucial for its adaptation to the fluctuating intertidal environment and presents significant potential in improving taurine nutrition and boosting immunity in humans. Cysteine dioxygenase (CDO) is a key enzyme involved in the initial step of taurine biosynthesis and plays a crucial role in regulating taurine content in the body. In the present study, polymorphisms of CDO gene in C. gigas (CgCDO) and their association with taurine content were evaluated in 198 individuals. A total of 24 single nucleotide polymorphism (SNP) loci were identified in the exonic region of CgCDO gene by direct sequencing. Among these SNPs, c.279G>A and c.287C>A were found to be significantly associated with taurine content, with the GG and AA genotype at the two loci exhibiting enhanced taurine accumulation (p < 0.05). Haplotype analysis revealed that the 279GG/287AA haplotype had the highest taurine content of 29.24 mg/g, while the 279AA/287CC haplotype showed the lowest taurine content of 21.19 mg/g. These results indicated that the SNPs of CgCDO gene could influence the taurine content in C. gigas and have potential applications in the selective breeding of high-taurine varieties.

Determination, expression and characterization of an UDP-N-acetylglucosamine:α-1,3-D-mannoside ß-1,2-N-acetylglucosaminyltransferase I (GnT-I) from the Pacific oyster, Crassostrea gigas.

Molluscs are intermediate hosts for several parasites. The recognition processes, required to evade the host's immune response, depend on carbohydrates. Therefore, the investigation of mollusc glycosylation capacities is of high relevance to understand the interaction of parasites with their host. UDP-N-acetylglucosamine:α-1,3-D-mannoside ß-1,2-N-acetylglucosaminyltransferase I (GnT-I) is the key enzyme for the biosynthesis of hybrid and complex type N-glycans catalysing the transfer of N-acetylglucosamine from UDP-N-acetylglucosamine to the α-1,3 Man antenna of Man5GlcNAc2. Thereby, the enzyme produces a suitable substrate for further enzymes, such as α-mannosidase II, GlcNAc-transferase II, galactosyltransferases or fucosyltransferases. The sequence of GnT- I from the Pacific oyster, Crassostrea gigas, was obtained by homology search using the corresponding human enzyme as the template. The obtained gene codes for a 445 amino acids long type II transmembrane glycoprotein and shared typical structural elements with enzymes from other species. The enzyme was expressed in insect cells and purified by immunoprecipitation using protein A/G-plus agarose beads linked to monoclonal His-tag antibodies. GnT-I activity was determined towards the substrates Man5-PA, MM-PA and GnM-PA. The enzyme displayed highest activity at pH 7.0 and 30 °C, using Man5-PA as the substrate. Divalent cations were indispensable for the enzyme, with highest activity at 40 mM Mn2+, while the addition of EDTA or Cu2+ abolished the activity completely. The activity was also reduced by the addition of UDP, UTP or galactose. In this study we present the identification, expression and biochemical characterization of the first molluscan UDP-N-acetylglucosamine:α-1,3-D-mannoside ß-1,2-N-acetylglucosaminyltransferase I, GnT-I, from the Pacific oyster Crassostrea gigas.

Molecular confirmation and functional study of signal transducer and activator of transcription genes in the Pacific oyster Crassostrea gigas.

The JAK (Janus kinase)-STAT (Signal transducer and activator of transcription) is a well-known functional signaling pathway that plays a key role in several important biological activities such as apoptosis, cell proliferation, differentiation, and immunity. However, limited studies have explored the functions of STAT genes in invertebrates. In the present study, the gene sequences of two STAT genes from the Pacific oyster (Crassostrea gigas), termed CgSTAT-Like-1 (CgSTAT-L1) and CgSTAT-Like-2 (CgSTAT-L2), were obtained using polymerase chain reaction (PCR) amplification and cloning. Multiple sequence comparisons revealed that the sequences of crucial domains of these proteins were conserved, and the similarity with the protein sequence of other molluscan STAT is close to 90 %. The phylogenetic analyses indicated that CgSTAT-L1 and CgSTAT-L2 are novel members of the mollusk STAT family. Quantitative real-time PCR results implied that CgSTAT-L1 and CgSTAT-L2 mRNA expression was found in all tissues, and significantly induced after challenge with lipopolysaccharide (LPS), peptidoglycan (PGN), or poly(I:C). After that, dual-luciferase reporter assays denoted that overexpression of CgSTAT-L1 and CgSTAT-L2 significantly activated the NF-κB signaling, and, interestingly, the overexpressed CgSTAT proteins potentiated LPS-induced NF-κB activation. These results contributed a preliminary analysis of the immune-related function of STAT genes in oysters, laying the foundation for deeper understanding of the function of invertebrate STAT genes.

The expression patterns of exosomal miRNAs in the Pacific oyster after high-temperature stress or Vibrio stimulation.

The exosomal miRNA plays a crucial role in the intercellular communication response to environmental stress and pathogenic stimulation. In the present study, the expression of exosomal miRNAs in the Pacific oyster Crassostrea gigas after high-temperature stress or Vibrio splendidus stimulation was investigated through high-throughput sequencing. The exosomes were identified to be teardrop-like vesicles with the average size of 81.7 nm by transmission electron microscopy. There were 66 known miRNAs and 33 novel miRNAs identified, of which 10 miRNAs were differentially expressed after both high-temperature stress and Vibrio stimulation compared to the control group. A total of 1868 genes were predicted as the putative targets of miRNAs, of which threonine aspartase 1-like was targeted by the highest number of related miRNAs. The robustness and reliability of miRNA expression from the sRNA sequencing data were verified by employing eight miRNAs for qPCR. GO and KEGG clustering analyses revealed that apoptosis was significantly enriched by the target genes of differentially expressed exosomal miRNAs after high-temperature stress, and autophagy and cytokine activity were significantly enriched after Vibrio stimulation. Energy metabolism was found to be significantly shared in the target gene enrichments after both high-temperature stress and Vibrio stimulation. These findings would improve our understanding of the regulatory mechanisms of exosomal miRNAs in C. gigas after high-temperature stress or Vibrio stimulation.

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.

Transcriptional and post-translational regulation of MITF mediated by bHLH domain during the melanogenesis and melanocyte proliferation in Crassostrea gigas.

Melanocyte differentiation is orchestrated by the master regulator transcription factor MITF. However, its ability to discern distinct binding sites linked to effective gene regulation remains poorly understood. This study aims to assess how co-activator acetyltransferase interacts with MITF to modulate their related lysine action, thereby mediating downstream gene regulation, including DNA affinity, stability, transcriptional activity, particularly in the process of shell pigmentation. Here, we have demonstrated that the CgMITF protein can be acetylated, further enabling selective amplification of the melanocyte maturation program. Collaboration with transcriptional co-regulator p300 advances MITF dynamically interplay with downstream targeted gene promoters. We have established that MITF activation was partially dependent on the bHLH domain, which was well conserved across species. The bHLH domain contained conserved lysine residues, including K6 and K43, which interacted with the E-box motif of downstream targeted-genes. Mutations at K6 and K43 lead to a decrease in the binding affinity of the E-box motif. CgMITF protein bound to the E-box motif within the promoter regions of the tyrosinase-related genes, contributing to melanogenesis, and also interacted with the E-box motif within the TBX2 promoter regions, associated with melanocyte proliferation. We elucidated how the bHLH domain links the transcriptional regulation and acetylation modifications in the melanocyte development in C. gigas.