Integration of Next Generation Sequencing and EPR Analysis to Uncover Molecular Mechanism Underlying Shell Color Variation in Scallops

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.

Outer fold is sole effective tissue among three mantle folds with regard to oyster shell colour

Molluscs constitute the second largest phylum of animals in the world, and shell colour is one of their most important phenotypic characteristics. In this study, we found among three folds on the mantle edge of oyster, only the outer fold had the same colour as the shell. Transcriptome and mantle cutting experiment indicated that the outer fold may be mainly reflected in chitin framework formation and biomineralisation. There were obvious differences in SEM structure and protein composition between the black and white shell periostraca. The black shell periostraca had more proteins related to melanin biosynthesis and chitin binding. Additionally, we identified an uncharacterized protein gene (named as CgCBP) ultra-highly expressed only in the black outer fold and confirmed its function of chitin-binding and CaCO₃ precipitation promoting. RNAi also indicated that CgCBP knockdown could change the structure of shell periostracum and reduce shell pigmentation. All these results suggest that the mantle outer fold plays multiple key roles in shell periostraca bioprocessing, and shell periostracum structure affected by chitin-binding protein is functionally correlated with shell pigmentation. The investigation of oyster shell periostracum structure and shell colour will provide a better understanding in pigmentation during biological mineralisation in molluscs.

Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing

BACKGROUND

Mollusca, a phylum of highly rich species, possess vivid shell colours, but the underlying molecular mechanism remains to be elucidated. DNA methylation, one of the most common epigenetic modifications in eukaryotes, is believed to play a vital role in various biological processes. However, analysis of the effects of DNA methylation on shell colouration has rarely been performed in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation.

RESULTS

In the present study, to reveal the role of epigenetic regulation in shell colouration, WGBS, the "gold standard" of DNA methylation analysis, was first performed on the mantle tissues of Yesso scallops (Patinopecten yessoensis) with different shell colours (brown and white), and DNA methylomes at single-base resolution were generated. About 3% of cytosines were methylated in the genome of the Yesso scallop. A slight increase in ^mCG percentage and methylation level was found in brown scallops. Sequence preference of nearby methylated cytosines differed between high and low methylation level sites and between the brown- and white-shelled scallops. DNA methylation levels varied among the different genomic regions; all the detected regions in the brown group exhibited higher methylation levels than the white group. A total of 41,175 DMRs (differentially methylated regions) were detected between brown and white scallops. GO functions and pathways associated with the biosynthesis of melanin and porphyrins were significantly enriched for DMRs, among which several key shell colour-related genes were identified. Further, different correlations between mRNA expression levels and DNA methylation status were found in these genes, suggesting that DNA methylation regulates shell colouration in the Yesso scallop.

CONCLUSIONS

This study provides genome-wide DNA methylation landscapes of Yesso scallops with different shell colours, offering new insights into the epigenetic regulatory mechanism underlying shell colour.

Integrated Analysis of Coding Genes and Non-coding RNAs Associated with Shell Color in the Pacific Oyster (Crassostrea gigas)

Molluscan shell color polymorphism is important in genetic breeding, while the molecular information mechanism for shell coloring is unclear. Here, high-throughput RNA sequencing was used to compare expression profiles of coding and non-coding RNAs (ncRNAs) from Pacific oyster Crassostrea gigas with orange and black shell, which were from an F2 family constructed by crossing an orange shell male with a black shell female. First, 458, 13, and 8 differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) were identified, respectively. Functional analysis suggested that the DEGs were significantly enriched in 9 pathways including tyrosine metabolism and oxidative phosphorylation pathways. Several genes related to melanin synthesis and biomineralization expressed higher whereas genes associated with carotenoid pigmentation or metabolism expressed lower in orange shell oyster. Then, based on the ncRNA analysis, 163 and 20 genes were targeted by 13 and 8 differentially expressed lncRNAs (DELs) and miRNAs (DEMs), severally. Potential DELs-DEMs-DEGs interactions were also examined. Seven DEMs-DEGs pairs were detected, in which tyrosinase-like protein 1 was targeted by lgi-miR-133-3p and lgi-miR-252a and cytochrome P450 was targeted by dme-miRNA-1-3p. These results revealed that melanin synthesis-related genes and miRNAs-mRNA interactions functioned on orange shell coloration, which shed light on the molecular regulation of shell coloration in marine shellfish.

Genome-wide identification, characterisation and expression analysis of the ALAS gene in the Yesso scallop (Patinopecten yessoensis) with different shell colours

Porphyrins, one of the most common shell pigments, are by-products of the haem pathway. 5-Aminolaevulinate synthase (ALAS) is the first and rate-limiting enzyme in this pathway and has been well studied in vertebrate species. However, the function of ALAS in shell colouration has been poorly studied in molluscs, which are renowned for their colourful shells. In the present study, an ALAS gene, named PyALAS, was identified through whole-genome scanning in the Yesso scallop (Patinopecten yessoensis), an economically and evolutionarily important bivalve species in which the shell colour represents polymorphism. Two conserved domains were detected in the PyALAS protein sequence, including a Preseq-ALAS domain and a 5-ALAS domain, confirming the identification of PyALAS. Phylogenetic analysis of the ALAS proteins among various invertebrate and vertebrate species revealed a high consistency between the molecular evolution of ALAS and the species taxonomy. PyALAS was ubiquitously expressed in most adult tissues of the Yesso scallop. The left mantle expressed a significantly higher level of PyALAS than the right side in brown scallops, whereas there was no significant difference in white scallops. Significantly different expression levels of PyALAS was also detected between the two different shell colour strains. These data indicate that PyALAS plays an important role in shell colouration in Yesso scallops and the present study provides new insights into the molecular mechanism of shell colouration in molluscs.

A Novel Tyrosinase Gene Plays a Potential Role in Modification the Shell Organic Matrix of the Triangle Mussel Hyriopsis cumingii

Although tyrosinases have been speculated to participate in the shell formation of mollusks, there is still a lack of experimental evidence to support this assumption. In this study, a novel tyrosinase designated HcTyr2 was isolated and characterized from the freshwater mussel Hyriopsis cumingii. The change in HcTyr2 mRNA expression during the process of embryonic development was detected by real-time quantitative PCR. The result showed that the expression of HcTyr2 mRNA was significantly upregulated at the stages of gastrulae and unmatured glochidia (P < 0.05), suggesting that this gene might fundamentally participate in the biogenesis and growth of the initial shell. Meanwhile, the upregulation of HcTyr2 mRNA at the stages of shell regeneration 24 h and 9 days after shell notching in the mantle edge (P < 0.05) implied that it might play an important role in shell periostracum and nacre formation by mediating the cross-linking of quinoproteins to promote the maturity of organic matrix. Additionally, the knockdown of HcTyr2 mRNA by RNA interference resulted in not only the suppression of periostracum growth but also structural disorder of nacre aragonite tablets, as detected by scanning electron microscopy. These results suggested that HcTyr2 might regulate the growth of shell by its oxidative ability to transform soluble matrix proteins into insoluble matrix proteins, then promoting the maturity of the shell organic framework in H. cumingii. In general, our results suggested the importance of HcTyr2 in the shell formation and regeneration of H. cumingii.

Histological and Expression Differences Among Different Mantle Regions of the Yesso Scallop (Patinopecten yessoensis) Provide Insights into the Molecular Mechanisms of Biomineralization and Pigmentation

The molecular mechanisms of shell formation and pigmentation are issues of great interest in molluscan studies due to the unique physical and biological properties of shells. The Yesso scallop, Patinopecten yessoensis, is one of the most important maricultural bivalves in Asian countries, and its shell color shows polymorphism. To gain more information about the underlying mechanisms of shell formation and pigmentation, this study presents the first analyses of histological and transcriptional differences between different mantle regions of the Yesso scallop, which are thought to be responsible for the formation of different shell layers. The results showed major microstructural differences between the edge and central mantles, which were closely associated with their functions. Different biomineralization-related GO functions, which might participate in the formation of different shell layers, were significantly enriched in the different mantle regions, indicating the different molecular functions of the two mantle regions in shell formation. The melanogenesis pathway, which controls melanin biosynthesis, was the most significantly enriched pathway in the DEGs between the two mantle regions, indicating its important role in shell pigmentation. Tyr, the key and rate-limiting gene in melanogenesis, was expressed at a remarkably high level in the central mantle, while the upstream regulatory genes included in melanogenesis were mainly upregulated in the edge mantle, suggesting the different molecular functions of the two mantle regions in shell pigmentation.

Genome-wide identification, characterization and expression analysis of the MITF gene in Yesso scallops (Patinopecten yessoensis) with different shell colors

The microphthalmia-associated transcription factor (MITF) is the center of the regulator network of melanin synthesis in vertebrates. However, the role of MITF in shell color formation is poorly studied in mollusks. In the present study, an MITF gene, PyMITF, was first identified at the whole-genome level in Yesso scallop (Patinopecten yessoensis), an evolutionarily and economically important species, the shell color of which shows polymorphism. The PyMITF is a large gene spanning ~37 kb in the genome with 7 introns and 8 exons. A basic helix-loop-helix leucine zipper (bHLH-LZ) domain was detected in the PyMITF protein sequence, which can bind the canonical E-box sequence in the promoter region of the downstream genes. Phylogenetic analysis of the MITFs among vertebrates and invertebrates revealed that the molecular evolution of MITFs was consistent with the species taxonomy. Different expression levels of PyMITF were detected among different shell color strains, indicating the important role of PyMITF involved in shell pigmentation. Besides, PyMITF was expressed at a significantly higher level in the central mantle than that in the edge mantle, proving the participation of the central mantle in shell color formation in molecular level for the first time. The work provides valuable information for the molecular mechanism study of shell color formation.

Differences between fast and slow muscles in scallops revealed through proteomics and transcriptomics

BACKGROUND

Scallops possess striated and catch adductor muscles, which have different structure and contractile properties. The striated muscle contracts very quickly for swimming, whereas the smooth catch muscle can keep the shells closed for long periods with little expenditure of energy. In this study, we performed proteomic and transcriptomic analyses of differences between the striated (fast) and catch (slow) adductor muscles in Yesso scallop Patinopecten yessoensis.

RESULTS

Transcriptomic analysis reveals 1316 upregulated and 8239 downregulated genes in slow compared to fast adductor muscle. For the same comparison, iTRAQ-based proteomics reveals 474 differentially expressed proteins (DEPs), 198 up- and 276 downregulated. These DEPs mainly comprise muscle-specific proteins of the sarcoplasmic reticulum, extracellular matrix, and metabolic pathways. A group of conventional muscle proteins-myosin heavy chain, myosin regulatory light chain, myosin essential light chain, and troponin-are enriched in fast muscle. In contrast, paramyosin, twitchin, and catchin are preferentially expressed in slow muscle. The association analysis of proteomic and transcriptomic data provides the evidences of regulatory events at the transcriptional and posttranscriptional levels in fast and slow muscles. Among 1236 differentially expressed unigenes, 22.7% show a similar regulation of mRNA levels and protein abundances. In contrast, more unigenes (53.2%) exhibit striking differences between gene expression and protein abundances in the two muscles, which indicates the existence of fiber-type specific, posttranscriptional regulatory events in most of myofibrillar proteins, such as myosin heavy chain, titin, troponin, and twitchin.

CONCLUSIONS

This first, global view of protein and mRNA expression levels in scallop fast and slow muscles reveal that regulatory mechanisms at the transcriptional and posttranscriptional levels are essential in the maintenance of muscle structure and function. The existence of fiber-type specific, posttranscriptional regulatory mechanisms in myofibrillar proteins will greatly improve our understanding of the molecular basis of muscle contraction and its regulation in non-model invertebrates.

Transcriptional profiling of long non-coding RNAs in mantle of Crassostrea gigas and their association with shell pigmentation

Long non-coding RNAs (lncRNAs) play crucial roles in diverse biological processes and have drawn extensive attention in the past few years. However, lncRNAs remain poorly understood about expression and roles in Crassostrea gigas, a potential model organism for marine molluscan studies. Here, we systematically identified lncRNAs in the mantles of C. gigas from four full-sib families characterized by white, black, golden, and partially pigmented shell. Using poly(A)-independent and strand-specific RNA-seq, a total of 441,205,852 clean reads and 12,243 lncRNA transcripts were obtained. LncRNA transcripts were relatively short with few exons and low levels of expression in comparison to protein coding mRNA transcripts. A total of 427 lncRNAs and 349 mRNAs were identified to differentially express among six pairwise groups, mainly involving in biomineralization and pigmentation through functional enrichment. Furthermore, a total of 6 mRNAs and their cis-acting lncRNAs were predicted to involve in synthesis of melanin, carotenoid, tetrapyrrole, or ommochrome. Of them, chorion peroxidase and its cis-acting lincRNA TCONS_00951105 are implicated in playing an essential role in the melanin synthetic pathway. Our studies provided the first systematic characterization of lncRNAs catalog expressed in oyster mantle, which may facilitate understanding the molecular regulation of shell colour diversity and provide new insights into future selective breeding of C. gigas for aquaculture.

A Genome-Wide Association Study Identifies the Genomic Region Associated with Shell Color in Yesso Scallop, Patinopecten yessoensis

The shell color polymorphism widely exists in economic shellfish, which not only results in a better visual perception but also shows great value as an economic trait for breeding. Small numbers of reddish-orange shell Yesso scallops, Patinopecten yessoensis, were found in cultured populations compared to the brown majority. In this study, a genome-wide association study was conducted to understand the genetic basis of shell color. Sixty-six 2b-RAD libraries with equal numbers of reddish-orange and brown shell individuals were constructed and sequenced using the Illumina HiSeq 2000 platform. A total of 322,332,684 high-quality reads were obtained, and the average sequencing depth was 18.4×. One genomic region on chromosome 11 that included 239 single-nucleotide polymorphisms (SNPs) was identified as significantly associated with shell color. After verification by high-resolution melting in another population, two SNPs were selected as specific loci for reddish-orange shell color. These two SNPs could be used to improve the selective breeding progress of true-breeding strains with complete reddish-orange scallops. In addition, within the significantly associated genomic region, candidate genes were identified using marker sequences to search the draft genome of Yesso scallop. Three genes (LDLR, FRIS, and FRIY) with known functions in carotenoid metabolism were identified. Further study using high-performance liquid chromatography proved that the relative level of carotenoids in the reddish-orange shells was 40 times higher than that in the brown shells. These results suggested that the accumulation of carotenoids contributes to the formation of reddish-orange shells.