Shell Biosynthesis and Pigmentation as Revealed by the Expression of Tyrosinase and Tyrosinase-like Protein Genes in Pacific Oyster (Crassostrea gigas) with Different Shell Colors

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.

Mitf involved in shell pigmentation by activating tyrosinase-mediated melanin synthesis in Pacific oyster (Crassostrea gigas)

Pigmentation is frequently observed in the molluscan shells, whereas the molecular regulation about these shell pigments formation is not clear. The microphthalmia-associated transcription factor (Mitf) is an important transactivator in melanin synthesis in vertebrates. Here, the Mitf containing a highly conserved basic helix-loop-helixleucine zipper (bHLH-LZ) domain was identified in an economically important marine bivalve Pacific oyster Crassostrea gigas. The Mitf was found to widespread tissue distribution and the expression was higher in the marginal mantle than in the central mantle. Particularly, the expression level of Mitf was high in black shell color oysters compared with white shell oysters. After injecting siRNA, the expression of Mitf decreased significantly, and the efficiency of RNA interference reached 53%. Besides, knockdown Mitf obviously decreased expression of tyrosinase family genes and tyrosinase activity of mantles, indicating a potential regulatory relationship between Mitf and Tyr or Typs. Simultaneously, there was a sharply reduce in the number of the melanosomes in the outer fold of mantle by silencing of Mitf. Luciferase assays in cell culture further verified that Mitf was involved in transcriptional regulation of Typ-2 and Typ-3 genes through binding to their specific promoter regions. These data argue that Mitf is involved in shell pigmentation through activating tyrosinase-mediated melanin synthesis in C. gigas.

Proteomic analysis of shell matrix proteins from the chiton Acanthopleura loochooana

Most molluscs have mineralized shells to protect themselves. Although the remarkable mechanical properties of shells have been well-studied, the origin of shells is still elusive. Chitons are unique in molluscs because they are shelly Aculifera which diverged from Conchifera (comprising all the shell-bearing classes of molluscs) in the early pre-Cambrian. We developed a method to extract shell proteins from chiton shell plates (removing embedded soft tissues) and then compared the shell proteome to that of Conchifera groups. Sixteen shell matrix proteins from Acanthopleura loochooana were identified by proteomics, in which Nacrein-like, Pif-like proteins, and protocadherin were found. Additional evidences from shell proteome in another species Chiton densiliratus and comparative sequence alignment in five chitons supported a conserved biomineralization toolkit in chitons. Our findings shed light on the evolution of mineralization in chitons and pose a hypothesis that ancestral molluscs have already evolved biomineralization toolkits, which may facilitate the formation of mineralized shells.

Comprehensive analysis of differentially expressed mRNA, lncRNA and miRNA, and their ceRNA networks in the regulation of shell color in the Manila clam (Ruditapes philippinarum)

The regulatory mechanism of ceRNA network plays an important role in molecular function and biological processes, however, the molecular mechanism in the shell color of Ruditapes philippinarum has not yet been reported. In this study, we performed transcriptome sequencing on the mantle of R. philippinarum with different shell colors, and screened for mRNA, miRNA, and lncRNA. A total of 61 mRNAs, 3725 lncRNAs and 90 miRNAs were obtained from all the shell color comparison groups (all mRNAs, lncRNAs and miRNAs P < 0.05), and 7 mRNAs, 8 lncRNAs, and 4 miRNAs of the porphyrin pathway and melanin pathway were screened for competitive endogenous RNA (ceRNA) network construction. The results indicate that the ceRNA network composed of mRNA and lncRNA, centered around efu-miR-101, mle-bantam-3p, egr-miR-9-5p, and sma-miR-75p, may play a crucial regulatory role in shell color formation. This study reveals for the first time the mechanism of ceRNA regulatory networks in the shell color of R. philippinarum and providing important reference data for molecular breeding of shell color in R. philippinarum.

Ioxynil and diethylstilbestrol impair cardiac performance and shell growth in the mussel Mytilus coruscus

The herbicide ioxynil (IOX) and the synthetic estrogen diethylstilbestrol (DES) are environmentally relevant contaminants that act as endocrine disruptors (EDCs) and have recently been shown to be cardiovascular disruptors in vertebrates. Mussels, Mytilus coruscus, were exposed to low doses of IOX (0.37, 0.037 and 0.0037 mg/L) and DES (0.27, 0.027 and 0.0027 mg/L) via the water and the effect monitored by generating whole animal transcriptomes and measuring cardiac performance and shell growth. One day after IOX (0.37 and 0.037 mg/L) and DES (0.27 and 0.027 mg/L) exposure heart rate frequency was decreased in both groups and 0.27 mg/L DES significantly reduced heart rate frequency with increasing time of exposure (P < 0.05) and no acclimatization occurred. The functional effects were coupled to significant differential expression of genes of the serotonergic synapse pathway and cardiac-related genes at 0.027 mg/L DES, which suggests that impaired heart function may be due to interference with neuroendocrine regulation and direct cardiac effect genes. Multiple genes related to detoxifying xenobiotic substances were up regulated and genes related to immune function were down regulated in the DES group (vs. control), indicating that detoxification processes were enhanced, and the immune response was depressed. In contrast, IOX had a minor disrupting effect at a molecular level. Of note was a significant suppression (P < 0.05) by DES of shell growth in juveniles and lower doses (< 0.0027 mg/L) had a more severe effect. The shell growth depression in 0.0027 mg/L DES-treated juveniles was not accompanied by abundant differential gene expression, suggesting that the effect of 0.0027 mg/L DES on shell growth may be direct. The results obtained in the present study reveal for the first time that IOX and DES may act as neuroendocrine disrupters with a broad spectrum of effects on cardiac performance and shell growth, and that DES exposure had a much more pronounced effect than IOX in a marine bivalve.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Genomic Insights into Mollusk Terrestrialization: Parallel and Convergent Gene Family Expansions as Key Facilitators in Out-of-the-Sea Transitions

Animals abandoned their marine niche and successfully adapted to life on land multiple times throughout evolution, providing a rare opportunity to study the mechanisms driving large scale macroevolutionary convergence. However, the genomic factors underlying this process remain largely unknown. Here, we investigate the macroevolutionary dynamics of gene repertoire evolution during repeated transitions out of the sea in mollusks, a lineage that has transitioned to freshwater and terrestrial environments multiple independent times. Through phylogenomics and phylogenetic comparative methods, we examine ∼100 genomic data sets encompassing all major molluskan lineages. We introduce a conceptual framework for identifying and analyzing parallel and convergent evolution at the orthogroup level (groups of genes derived from a single ancestral gene in the species in question) and explore the extent of these mechanisms. Despite deep temporal divergences, we found that parallel expansions of ancient gene families played a major role in facilitating adaptation to nonmarine habitats, highlighting the relevance of the preexisting genomic toolkit in facilitating adaptation to new environments. The expanded functions primarily involve metabolic, osmoregulatory, and defense-related systems. We further found functionally convergent lineage-exclusive gene gains, while family contractions appear to be driven by neutral processes. Also, genomic innovations likely contributed to fuel independent habitat transitions. Overall, our study reveals that various mechanisms of gene repertoire evolution-parallelism, convergence, and innovation-can simultaneously contribute to major evolutionary transitions. Our results provide a genome-wide gene repertoire atlas of molluskan terrestrialization that paves the way toward further understanding the functional and evolutionary bases of this process.

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.

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.

Transcriptomic responses reveal impaired physiological performance of the pearl oyster following repeated exposure to marine heatwaves

Marine heatwaves are predicted to become more intense and frequent in the future, possibly threatening the survival of marine organisms and devastating their communities. While recent evidence reveals the adaptability of marine organisms to heatwaves, substantially overlooked is whether they can also adjust to repeated heatwave exposure, which can occur in nature. By analysing transcriptome, we examined the fitness and recoverability of the pearl oyster (Pinctada maxima) after two consecutive heatwaves (24 °C to 32 °C for 3 days; recovery at 24 °C for 4 days). In the first heatwave, 331 differentially expressed genes (DEGs) were found, such as AGE-RAGE, MAPK, JAK-STAT, FoxO and mTOR. Despite the recovery after the first heatwave, 2511 DEGs related to energy metabolism, body defence, cell proliferation and biomineralization were found, where 1655 of them were downregulated, suggesting a strong negative response to the second heatwave. Our findings imply that some marine organisms can indeed tolerate heatwaves by boosting energy metabolism to support molecular defence, cell proliferation and biomineralization, but this capacity can be overwhelmed by repeated exposure to heatwaves. Since recurrence of heatwaves within a short period of time is predicted to be more prevalent in the future, the functioning of marine ecosystems would be disrupted if marine organisms fail to accommodate repeated extreme thermal stress.

Regulation of Tyrosinase Gene Expression by Retinoic Acid Pathway in the Pacific Oyster Crassostrea gigas

Retinoic acid (RA) plays important roles in various biological processes in animals. RA signaling is mediated by two types of nuclear receptors, namely retinoic acid receptor (RAR) and retinoid x receptor (RXR), which regulate gene expression by binding to retinoic acid response elements (RAREs) in the promoters of target genes. Here, we explored the effect of all-trans retinoic acid (ATRA) on the Pacific oyster Crassostera gigas at the transcriptome level. A total of 586 differentially expressed genes (DEGs) were identified in C. gigas upon ATRA treatment, with 309 upregulated and 277 downregulated genes. Bioinformatic analysis revealed that ATRA affects the development, metabolism, reproduction, and immunity of C. gigas. Four tyrosinase genes, including Tyr-6 (LOC105331209), Tyr-9 (LOC105346503), Tyr-20 (LOC105330910), and Tyr-12 (LOC105320007), were upregulated by ATRA according to the transcriptome data and these results were verified by real-time quantitative polymerase chain reaction (RT-qPCR) analysis. In addition, increased expression of Tyr (a melanin-related TYR gene in C. gigas) and Tyr-2 were detected after ATRA treatment. The yeast one-hybrid assay revealed the DNA-binding activity of the RA receptors CgRAR and CgRXR, and the interaction of CgRAR with RARE present in the Tyr-2 promoter. These results provide evidence for the further studies on the role of ATRA and the mechanism of RA receptors in mollusks.

Expression of tyrosinase-like protein genes and their functional analysis in melanin synthesis of Pacific oyster (Crassostrea gigas)

Color polymorphism in Mollusca is of great interest for consumer preference. Although the heritability of shell color variation has been conducted by experimental crossing, little is known about molecular basis involved in these patterns. Tyrosinase-like proteins are important enzymes which are members of the type-3 copper protein superfamily. In this research, two tyrosinase-like protein genes including CgTyp-1 and CgTyp-3 were identified in the Pacific oyster Crassostrea gigas. Tissue expression analysis showed that CgTyp-1 and CgTyp-3 were dominantly expressed in the mantle. Particularly, they were expressed significantly higher in the edge mantle than that in the central mantle whether on the left or right mantles. Additionally, expressions of CgTyp-1 and CgTyp-3 were mainly found in the black shell color oysters, with relative lower levels in the white shell color oysters. In situ hybridization showed that positive signals for CgTyp-1 and CgTyp-3 were both detected within the outer epithelium of the outer fold either in the black or white shell color oysters. After interference, the expression levels of CgTyp-1 and CgTyp-3 mRNA were significantly attenuated, and the efficiency of RNAi reached 84.72% and 71.58%, respectively. Besides, knockdown CgTyp-1 or CgTyp-3, obviously decreased the tyrosinase activity of mantles. Furthermore, the number of the melanosomes within epithelium of the outer fold was sharply reduced by silencing of each Typ. These findings argue that CgTyp-1 and CgTyp-3 may be involved in the melanin synthesis, which lends insight into regulation mechanism of shell pigmentation in C. gigas.

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.

Heme-Peroxidase 2 Modulated by POU2F1 and SOX5 is Involved in Pigmentation in Pacific Oyster (Crassostrea gigas)

Color polymorphism is frequently observed in molluscan shellfish, while the molecular regulation of shell pigmentation is not well understood. Peroxidase is a key enzyme involved in melanogenesis. Here, we identified a heme-peroxidase 2 gene (CgHPX2), and characterized the expression patterns and transcriptional regulation of CgHPX2 in the Pacific oyster Crassostrea gigas. Tissues expression analysis showed that CgHPX2 was a mantle-specific gene and primarily expressed in the edge mantle in black shell color oyster compared with white shell oyster. In situ hybridization showed that strong signals for CgHPX2 were detected in the both inner and outer surface of the outer fold of mantle in the black shell color oyster, whereas positive signals in white shell oyster were mainly localized in the outer surface of the outer fold of mantle. In the embryos and larvae, a high expression level of CgHPX2 was detected in the trochophore stage in both black and white shell color oysters. The temporal localization of CgHPX2 was mainly detected in the shell gland and edge mantle of trochophore and calcified shell larvae, respectively. In addition, a 2227 bp of 5' flanking region sequence of CgHPX2 was cloned, which contained a presumed core promoter region and many potential transcription factor binding sites. Further luciferase assay experiment confirmed that POU domain, class 2, transcription factor 1 (POU2F1), and SRY-box transcription factor 5 (SOX5) were involved in transcriptional regulation of CgHPX2 gene through binding to its specific promoter region. After CgPOU2F1 and CgSOX5 RNA interference, the CgHPX2 gene expression was significantly decreased. These results suggested that CgPOU2F1 and CgSOX5 might be two important transcription factors that positively regulated the expression of CgHPX2 gene, improving our understanding of the transcriptional regulation of molluscan shell pigmentation.

MiRNA-mRNA Integration Analysis Reveals the Regulatory Roles of MiRNAs in Shell Pigmentation of the Manila clam (Ruditapes philippinarum)

The shell color of the Manila clam (Ruditapes philippinarum) is an economically important trait. We used high-throughput sequencing and transcriptome analysis to study the molecular mechanisms that underlie shell color formation and regulation in this species. We constructed small RNA libraries from mantle tissues from four shell color strains of Manila clam, subjected them to high-throughput sequencing. Notably, the results suggested that a number of pigment-associated genes including Mitf, HERC2, were negatively regulated by nvi-miR-2a, tgu-miR-133-3p, respectively. They might be involved in melanin formation via the activation of the melanogenesis pathway. And aae-miR-71-5p and dme-miR-7-5p linked to shell formation-related genes such as Calmodulin and IMSP3 were considered to participate in the calcium signaling pathway. We then used quantitative PCR to verify the candidate miRNAs and target genes in different shell color groups. Our results indicated that miR-7, miR-71, and miR-133 may regulate target mRNAs to participate in shell color pigmentation. These results provide the foundation to further characterize miRNA effects on the regulation of shell color and have significant implications for the breeding of new varieties of clams.