Comparative analysis of the transcriptome in tissues secreting purple and white nacre in the pearl mussel Hyriopsis cumingii

HcN57, A Novel Unusual Acidic Silk-Like Matrix Protein from Hyriopsis cumingii, Participates in Framework Construction and Nacre Nucleation During Nacreous Layer Formation

In the classic molecular model of nacreous layer formation, unusual acidic matrix proteins rich in aspartic acid (Asp) residues are essential for nacre nucleation due to their great affinity for binding calcium. However, the acidic matrix proteins discovered in the nacreous layer so far have been weakly acidic with a high proportion of glutamate. In the present study, several silk-like matrix proteins, including the novel matrix protein HcN57, were identified in the ethylenediaminetetraacetic acid-soluble extracts of the nacreous layer of Hyriopsis cumingii. HcN57 is a highly repetitive protein that consists of a high proportion of alanine (Ala, 34.4%), glycine (Gly, 22.5%), and serine (Ser, 11.4%). It forms poly Ala blocks, GlynX repeats, an Ala-Gly repeat, and a Ser-Ala-rich region, exhibiting significant similarity to silk proteins found in spider species. The expression of HcN57 was specifically located in the dorsal epithelial cells of the mantle pallium and mantle center. Notably, expression of HcN57 was relatively high during nacreous layer regeneration and pearl nacre deposition, suggesting HcN57 is a silk matrix protein in the nacreous layer. Importantly, HcN57 also contains a certain content of Asp residues, making it an unusual acidic matrix protein present in the nacreous layer. These Asp residues are mainly distributed in three large hydrophilic acidic regions, which showed inhibitory activity against aragonite deposition and morphological regulation of calcite in vitro. Moreover, HcN57-dsRNA injection resulted in failure of nacre nucleation in vivo. Taken together, our results show that HcN57 is a bifunctional silk protein with poly Ala blocks and Gly-rich regions that serve as space fillers within the chitinous framework to prevent crystallization at unnecessary nucleation sites and Asp-rich regions that create a calcium ion supersaturated microenvironment for nucleation in the center of nacre tablets. These observations contribute to a better understanding of the mechanism by which silk proteins regulate framework construction and nacre nucleation during nacreous layer formation.

The lgi-miR-2d is Potentially Involved in Shell Melanin Synthesis by Targeting mitf in Manila Clam Ruditapes philippinarum

Shell color as an important economic trait is also the crucial target trait for breeding and production. MicroRNA (miRNA) is an endogenous small non-coding RNA that can post-transcriptionally regulate the expression of target genes, it plays important roles in many life activities and physiological processes, such as shell color, stress response, and disease traits. In this study, we investigated the function of lgi-miR-2d in shell melanin formation and the expression patterns of lgi-miR-2d and target gene Rpmitf in Manila clam Ruditapes philippinarum. We further explored and verified the relationship between Rpmitf and lgi-miR-2d and identified the expression level of shell color-related gene changes by RNAi and injecting the antagomir of lgi-miR-2d, respectively. Our results indicated that lgi-miR-2d antagomir affected the expression of its target gene Rpmitf. In addition, the dual-luciferase reporter assay was conducted to confirm the direct interaction between lgi-miR-2d and Rpmitf. The results showed that the expression levels of melanin-related genes such as Rpmitf and tyr were significantly decreased in the positive treatment group compared with the blank control group after the Rpmitf dsRNA injection, indicating Rpmitf plays a crucial role in the melanin synthesis pathway. Taken together, we speculated that lgi-miR-2d might be negatively modulating Rpmitf, which might regulate other shell color-related genes, thereby affecting melanin synthesis in R. philippinarum.

Freshwater pearl culture in Bangladesh: Current status and prospects

Freshwater pearl farming is an emerging sector of aquaculture in Bangladesh which plays a growing role at major jewelry markets. With some improved techniques, high quality image or designer pearls are now produced from freshwater mussels Lamellidens marginalis. Yet it is difficult to reach in conclusion as the quantities produced, culture techniques used, and the upgrading of the existing culture technique are not well documented. Furthermore, many obstacles such as proper dissemination of culture technologies among the interested peoples, optimization of the culture environment and culture methods, standardization of breeding protocol and so on need to be addressed by the scientific community. This review article reports for the first time about the status of freshwater pearl culture in Bangladesh highlighting the fundamentals of pearl production, culture techniques used in farms, challenges, and prospects for upgradation of current culture principles in Bangladesh.

A Novel Kunitz-Type Serine Protease Inhibitor (HcKuSPI) is Involved in Antibacterial Defense in Innate Immunity and Participates in Shell Formation of Hyriopsis cumingii

Serine protease inhibitors (SPIs) are abundantly reported for its inhibition against specific proteases involved in the immune responses, but SPI data related to calcareous shells are scarce. Previously, our research group has reported the proteome analysis of non-nucleated pearl powder, and a candidate matrix protein containing two Kunitz domains in the acid soluble fraction caught our attention. In the present study, the full-length cDNA sequence of HcKuSPI was obtained from Hyriopsis cumingii. HcKuSPI was specifically expressed in the mantle, with hybridization signals mainly concentrated to dorsal epithelial cells at the mantle edge and weak signals at the mantle pallium, suggesting HcKuSPI was involved in shell formation. HcKuSPI expression in the mantle was upregulated after Aeromonas hydrophila and Staphylococcus aureus challenge to extrapallial fluids (EPFs). A glutathione S transferase (GST)-HcKuSPI recombinant protein showed strong inhibitory activity against the proteases, trypsin and chymotrypsin. Moreover, HcKuSPI expression in an experimental group was significantly higher when compared with a control group during pellicle growth and crystal deposition in shell regeneration processes, while the organic shell framework of newborn prisms and nacre tablets was completely destroyed after HcKuSPI RNA interference (RNAi). Therefore, HcKuSPI secreted by the mantle may effectively neutralize excess proteases and bacterial proteases in the EPF during bacterial infection and could prevent matrix protein extracellular degradation by suppressing protease proteolytic activity, thereby ensuring a smooth shell biomineralization. In addition, GST-HcKuSPI was also crucial for crystal morphology regulation. These results have important implications for our understanding of the potential roles of SPIs during shell biomineralization.

Population Genetics of Manila Clam (Ruditapes philippinarum) in China Inferred from Microsatellite Markers

The Manila clam (Ruditapes philippinarum) is one of the most commercially important bivalves along the coast of China. With the continuous expansion of clam farming scale, it may lead to some serious problems, including loss of genetic variation, inbreeding depression, and reduced effective population size (N_e). In the present study, eleven microsatellite markers were used to investigate the genetic diversity and differentiation among 13 clam populations along the coast of China. As a result, 150 alleles were detected according to the genotyping results of eleven microsatellite loci. The observed heterozygosity (H_o) was estimated to be ranging from 0.437 to 0.678, while the expected heterozygosity (H_e) was calculated to be varying from 0.587 to 0.700. F_st values between populations ranged from 0.0046-0.1983. In particular, the Laizhou population had the highest genetic variability, which was significantly different from the others (all F_st values > 0.1). For all the clam populations, there was no significant linear regression between genetic and geographic distance, indicating that these populations do not follow a pattern of isolation by distance (IBD). Genetic structure was estimated according to NJ, principal coordinates (PCoA), and structure-based clustering. Estimates of effective population size range from dozens to thousands among different populations, based on linkage-disequilibrium and molecular coancestry methods. The results reveal the genetic diversity of clams and verify the hypothesis that clam population differentiation may be influenced by the mode of southern breeding and northern culture, providing guiding information for natural resource conservation and genetic breeding of clams.

Identification and Functional Analysis of MAPKAPK2 in Hyriopsis cumingii

MAPKAPK2 (MK2) is an important regulator of the p38 mitogen-activated protein kinase (p38 MAPK) pathway, which is involved in a plethora of cellular processes concluding the development of gamete cells in meiosis and resisting pathogenic bacterial infestation. Hyriopsis cumingii is a significant mussel resource in China and a good material for pearl breeding. To explore the role of MK2 in H. cumingii, MK2 was identified and cloned, whose full-length cDNA was 1568 bp, including 87 bp in 5′ UTR, 398 bp in 3′ UTR, and 1083 bp in the open reading frame (ORF) region, encoding 360 amino acids. The expression of MK2 was the highest in the gills. Meanwhile, there was a significant difference in the gonads. After Aeromonas hydrophila and Lipopolysaccharide (LPS) infestation, the transcript level of the MK2 was upregulated in the gills. It indicated that MK2 might be involved in the innate immune response of H. cumingii after a pathogenic attack. After quantifying H. cumingii of different ages, it was found that the expression of MK2 was highest at 1 year old. In situ hybridization (ISH) results showed that the blue-purple hybridization signal was very significant in the oocytes and egg membranes of the female gonads of H. cumingii. The expression of MK2 increased gradually at the age of 1 to 5 months and showed a downward trend at the age of 5 to 8 months. It was suggested that MK2 might play an important role in the formation of primitive germ cells in H. cumingii. To sum up, MK2 might not only be involved in the immune response against pathogenic bacterial infection but also might play an important role in the development of the gonads in H. cumingii.

Analysis of protein kinase C (HcPKC) gene expression and single-nucleotide polymorphisms related to inner shell color traits in Hyriopsis cumingii

Background

Protein kinase C (PKC) is a multifunctional serine and PKC can phosphorylate serine residues in the cytoplasmic domain of tyrosinase, thereby regulating the activity of tyrosinase. Activated PKC is bound to the melanosome membrane, and unactivated PKC is free in the cytoplasm of melanocytes. In this study, we study the role of PKC gene in the melanin synthesis pathway and its effect on the color of the nacre of H. cumingii.

Results

In this study, a HcPKC gene in H. cumingii was cloned and its effects on melanin synthesis and nacre color were studied. HcPKC was expressed in both purple and white mussels, and the level of mRNA expression was higher in the purple mussels than in white mussels. Strong and specific mRNA signals were detected in the dorsal epithelial cells of the mantle pallial layer, indicating that HcPKC may be involved in nacre formation. After SNP association with inner shell color related traits, according to the principle that 0.25 < PIC < 0.5 is medium polymorphism and PIC < 0.25 is low polymorphism, the A + 332G site on the HcPKC gene was a site of moderate polymorphism, and the other four sites were low polymorphism sex sites. There was strong linkage disequilibrium among the five loci. A haplotype was constructed and it was found that the frequency of T1 (AGGAA)in the white population was significantly higher than that in the purple population (P < 0.05).

Conclusion

The study found that HcPKC of H. cumingii can be used as a candidate gene related to inner shell color, and some of the SNP sites can be used for molecular-assisted breeding in the spinnaker mussel, providing a reference for cultivating high-quality freshwater pearls.

Identification of a coproporphyrinogen-III oxidase gene and its correlation with nacre color in Hyriopsis cumingii

Pearl color is an important factor influencing pearl value, and is affected by the nacre color of the shell in Hyriopsis cumingii. Coproporphyrinogen-III oxidase (CPOX) is a key enzyme in porphyrin synthesis, and porphyrins are involved in color formation in different organisms, including in the nacre color of mussels. In this study, a CPOX gene (HcCPOX) was identified from H. cumingii, and its amino acid sequence was found to contain a coprogen-oxidase domain. HcCPOX mRNA was expressed widely in the tissues of white and purple mussels, and the highest expression was found in the gill, followed by the fringe mantle. The expression of HcCPOX in all tissues of purple mussels (except in the middle mantle) was higher than that of white mussels. Strong hybridization signals for HcCPOX were observed in the dorsal epithelial cells of the outer fold of the mantle. The activity of CPOX in the gill, fringe mantle, and foot of purple mussels was significantly higher than that in white mussels. Moreover, the expression of HcCPOX and CPOX activity were decreased in RNA interference experiments. The findings indicate that HcCPOX might contributes to nacre color formation in H. cumingii by being involved in porphyrin synthesis.

RNA Interference by Ingested Dsrna-Expressing Bacteria to Study Porphyrin Pigmentation in Crassostrea gigas

Porphyrins are a widespread group of pigments in nature which are believed to contribute to shell colors in mollusks. Previous studies have provided candidate genes for porphyrin shell coloration, however, the linkage analysis between functional genes and porphyrin pigmentation remains unclear in mollusks. RNA interference is a powerful molecular tool for analyzing the loss of functions of genes in vivo and alter gene expression. In this study, we used unicellular alga Platymonas subcordiformis and Nitzschia closterium f. minutissima as vectors to feed oysters with Escherichia coli strain HT115 engineered to express double-stranded RNAs targeting specific genes involved in porphyrin synthesis. A strain of Crassostrea gigas with orange shell was used to target key haem pathway genes expression using the aforementioned approach. We show here that feeding the oysters with E. coli, containing dsRNA targeting pigmentation genes, can cause changes in the color of the newly deposited shell. For example, the RNAi knockdown of CgALAS and CgPBGD resulted in the loss of uroporphyrin pigmentation from the shell due to the accumulation of the pigment in the oyster's mantle. The study probed the crucial role of ALAS and PBGD genes potential functions of uroporphyrin production and shell color pigmentation in C. gigas.

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.

Cloning of a HcCreb gene and analysis of its effects on nacre color and melanin synthesis in Hyriopsis cumingii

Creb (Cyclic AMP response element binding protein) is a nuclear regulatory factor that regulates transcription through autophosphorylation. In melanocytes, cAMP's corresponding elements bind to the Creb protein to autophosphorylation and activate MITF (Microphthalmia-associated transcription factor). MITF stimulates Tyrosine(tyr) to induce melanocytes to differentiate into eumelanin and pheomelanin. In this study, a HcCreb gene in Hyriopsis cumingii was cloned and its effects on melanin synthesis and nacre color were studied. HcCreb was expressed in both purple and white mussels, and there was a significant difference in expression between adductor muscle (p<0.01) and mantle tissue (p<0.05). Other tissues did not show significant differences (except for gill tissue), and in general, the level of mRNA expression was higher in purple mussels than in white mussels. In both white and purple mussels expression levels in gill tissue was the highest, followed by the mantle. Strong and specific mRNA signals were detected in the dorsal epithelial cells of the mantle pallial layer, indicating that HcCreb may be involved in nacre formation. After arbutin treatment, the expression of HcCreb decreased significantly. By further testing the changes in mantle melanin content it was found that the melanin content after arbutin treatment decreased significantly compared to the control group (p<0.05). It is speculated that the HcCreb gene plays a role in the process of melanin synthesis and nacre color formation in H. cumingii.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Eumelanin and pheomelanin pigmentation in mollusc shells may be less common than expected: insights from mass spectrometry

BACKGROUND

The geometric patterns that adorn the shells of many phylogenetically disparate molluscan species are comprised of pigments that span the visible spectrum. Although early chemical studies implicated melanin as a commonly employed pigment, surprisingly little evidence generated with more recent and sensitive techniques exists to support these observations.

RESULTS

Here we present the first mass spectrometric investigations for the presence of eumelanin and pheomelanin in 13 different molluscan species from three conchiferan classes: Bivalvia, Cephalopoda and Gastropoda. In the bivalve Mytilus edulis we demonstrate that eumelanin mainly occurs in the outermost, non-mineralised and highly pigmented layer of the shell (often referred to as the periostracum). We also identified eumelanin in the shells of the cephalopod Nautilus pompilius and the marine gastropods Clanculus pharaonius and Steromphala adriatica. In the terrestrial gastropod Cepaea nemoralis we verify the presence of pheomelanin in a mollusc shell for the first time. Surprisingly, in a large number of brown/black coloured shells we did not find any evidence for either type of melanin.

CONCLUSIONS

We recommend methods such as high-performance liquid chromatography with mass spectrometric detection for the analysis of complex biological samples to avoid potential false-positive identification of melanin. Our results imply that many molluscan species employ as yet unidentified pigments to pattern their shells. This has implications for our understanding of how molluscs evolved the ability to pigment and pattern their shells, and for the identification of the molecular mechanisms that regulate these processes.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

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.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

'Palaeoshellomics' reveals the use of freshwater mother-of-pearl in prehistory

The extensive use of mollusc shell as a versatile raw material is testament to its importance in prehistoric times. The consistent choice of certain species for different purposes, including the making of ornaments, is a direct representation of how humans viewed and exploited their environment. The necessary taxonomic information, however, is often impossible to obtain from objects that are small, heavily worked or degraded. Here we propose a novel biogeochemical approach to track the biological origin of prehistoric mollusc shell. We conducted an in-depth study of archaeological ornaments using microstructural, geochemical and biomolecular analyses, including ‘palaeoshellomics’, the first application of palaeoproteomics to mollusc shells (and indeed to any invertebrate calcified tissue). We reveal the consistent use of locally-sourced freshwater mother-of-pearl for the standardized manufacture of ‘double-buttons’. This craft is found throughout Europe between 4200–3800 BCE, highlighting the ornament-makers’ profound knowledge of the biogeosphere and the existence of cross-cultural traditions.

Just like people do today, prehistoric humans liked to adorn themselves with beautiful objects. Shells, from creatures like clams and snails, were used to decorate clothing or worn as jewelry at least as far back as 100,000 years ago. Later people used shells as the raw materials to make beads or bracelets. Learning where the shells came from may help scientists understand why prehistoric people chose certain shells and not others. It may also offer clues about how they used natural resources and the cultural significance of these objects. But identifying the shells is difficult because they lose many of their original distinctive features when worked into ornaments.

New tools that use DNA or proteins to identify the raw materials used to craft ancient artifacts have emerged that may help. So far, scientists have mostly used these genomic and proteomic tools to identify the source of materials made from animal hide, ivory or bone – where collagen is the most abundant protein molecule. Yet it is more challenging to extract and characterize proteins or genetic material from mollusc shells. This is partly because the amount of proteins in shells is at least 300 times lower than in bone, and also because the makeup of proteins in shells is not as well-known as in collagen.

Sakalauskaite et al. have now overcome these issues by combining the analytical tools used to study the proteins and mineral content of modern shells with those of ancient protein research. They then used this approach, which they named palaeoshellomics, to extract proteins from seven “double-buttons” – pearl-like ornaments crafted by prehistoric people in Europe. The double-buttons were made between 4200 and 3800 BC and found at archeological sites in Denmark, Germany and Romania. Comparing the extracted proteins to those from various mollusc shells showed that the double-buttons were made from freshwater mussels belonging to a group known as the Unionoida.

The discovery helps settle a decade-long debate in archeology about the origin of the shells used to make double-buttons in prehistoric Europe. Ancient people often crafted ornaments from marine shells, because they were exotic and considered more prestigious. But the results on the double-buttons suggest instead that mother-of-pearl from fresh water shells was valued and used by groups throughout Europe, even those living in coastal areas. The palaeoshellomics technique used by Sakalauskaite et al. may now help identify the origins of shells from archeological and palaeontological sites.

Transcriptional changes in the Japanese scallop (Mizuhopecten yessoensis) shellinfested by Polydora provide insights into the molecular mechanism of shell formation and immunomodulation

The Japanese scallop (Mizuhopecten yessoensis) is one of the most important aquaculture species in Asian countries; however, it has suffered severe infection by Polydora in northern China in recent years, causing great economic losses. The Polydora parasitizes the shell of scallops, badly destroying the shell's structure. To investigate the molecular response mechanism of M. yessoensis to Polydora infestion, a comprehensive and niche-targeted cDNA sequence database for diseased scallops was constructed. Additionally, the transcriptional changes in the edge mantle, central mantle and hemocytes, tissues directly related to the disease, were first described in this study. The results showed that genes involved in shell formation and immunomodulation were significantly differentially expressed due to the infestation. Different transcriptional changes existed between the two mantle regions, indicating the different molecular functions likely responsible for the formation of different shell layers. The differential expression of genes for immune recognition, signal transduction and pathogen elimination presented an integrated immune response process in scallops. Moreover, neuromodulation and glycometabolism involved in the regulation process with relevant function significantly enriched. The study provides valuable information for mechanism study of shell formation and immunomodulation in scallops.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4463-x) contains supplementary material, which is available to authorized users.

Colorful seashells: Identification of haem pathway genes associated with the synthesis of porphyrin shell color in marine snails

Very little is known about the evolution of molluskan shell pigments, although Mollusca is a highly diverse, species rich, and ecologically important group of animals comprised of many brightly colored taxa. The marine snail genus Clanculus was chosen as an exceptional model for studying the evolution of shell color, first, because in Clanculus margaritarius and Clanculus pharaonius both shell and foot share similar colors and patterns; and second, because recent studies have identified the pigments, trochopuniceus (pink‐red), and trochoxouthos (yellow‐brown), both comprised of uroporphyrin I and uroporphyrin III, in both shell and colored foot tissue of these species. These unusual characteristics provide a rare opportunity to identify the genes involved in color production because, as the same pigments occur in the shell and colored foot tissue, the same color‐related genes may be simultaneously expressed in both mantle (which produces the shell) and foot tissue. In this study, the transcriptomes of these two Clanculus species along with a third species, Calliostoma zizyphinum, were sequenced to identify genes associated with the synthesis of porphyrins. Calliostoma zizyphinum was selected as a negative control as trochopuniceus and trochoxouthos were not found to occur in this species. As expected, genes necessary for the production of uroporphyrin I and III were found in all three species, but gene expression levels were consistent with synthesis of uroporphyrins in mantle and colored foot tissue only in Clanculus. These results are relevant not only to understanding the evolution of shell pigmentation in Clanculus but also to understanding the evolution of color in other species with uroporphyrin pigmentation, including (mainly marine) mollusks soft tissues and shells, annelid and platyhelminth worms, and some bird feathers.

A galectin contributes to the innate immune recognition and elimination of pathogens in the freshwater mussel Hyriopsis cumingii

Galectins are members of the lectin superfamily. They function as pattern recognition receptors in the innate immune system of vertebrates and invertebrates. A galectin homolog from the triangle sail mussel Hyriopsis cumingii (HcGal2) was cloned and characterized. HcGal2 mRNA was expressed in all tissues examined, displaying particular enrichment in mantle tissue. Interestingly, rHcGAL2 protein was only detected in the mantle, hemocytes, and gills, suggesting that post-transcriptional regulation may occur. HcGal2 expression was induced in the mantle, liver, and hemocytes after exposure to lipopolysaccharides, Gram-negative bacteria (Aeromonas hydrophila), and Gram-positive bacteria (Staphylococcus aureus). The transcript significant upregulated was also detected after implantation in the mantle, pearl sac, liver, and hemocytes. Recombinant HcGAL2 protein (rHcGAL2) agglutinated Gram-positive and Gram-negative bacteria. In addition, rHcGAL2 promoted phagocytosis by hemocytes in vivo. Our data suggest that HcGal2 functioned as a pattern recognition receptor in against the pathogenic microbes and contributed to the "non-self" recognition and elimination in H. cumingii.

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.

Design strategies and applications of nacre-based biomaterials

The field of tissue engineering and regenerative medicine relies heavily on materials capable of implantation without significant foreign body reactions and with the ability to promote tissue differentiation and regeneration. The field of bone tissue engineering in particular requires materials capable of providing enhanced mechanical properties and promoting osteogenic cell lineage commitment. While bone repair has long relied almost exclusively on inorganic, calcium phosphate ceramics such as hydroxyapatite and their composites or on non-degradable metals, the organically derived shell and pearl nacre generated by mollusks has emerged as a promising alternative. Nacre is a naturally occurring composite material composed of inorganic, calcium carbonate plates connected by a framework of organic molecules. Similar to mammalian bone, the highly organized microstructure of nacre endows the composite with superior mechanical properties while the organic phase contributes to significant bioactivity. Studies, both in vitro and in vivo, have demonstrated nacre's biocompatibility, biodegradability, and osteogenic potential, which are superior to pure inorganic minerals such as hydroxyapatite or non-degradable metals. Nacre can be used directly as a bulk implant or as part of a composite material when combined with polymers or other ceramics. While nacre has demonstrated its effectiveness in multiple cell culture and animal models, it remains a relatively underexplored biomaterial. This review introduces the formation, structure, and characteristics of nacre, and discusses the present and future uses of this biologically-derived material as a novel biomaterial for orthopedic and other tissue engineering applications.

STATEMENT OF SIGNIFICANCE

Mussel derived nacre, a biological composite composed of mineralized calcium carbonate platelets and interplatelet protein components, has recently gained interest as a potential alternative ceramic material in orthopedic biomaterials, combining the integration and mechanical capabilities of calcium phosphates with increased bioactivity derived from proteins and biomolecules; however, there is limited awareness of this material's potential. Herein, we present, to our knowledge, the first comprehensive review of nacre as a biomaterial. Nacre is a highly promising yet overlooked biomaterial for orthopedic tissue engineering with great potential in a wide variety of material systems. It is our hope that publication of this article will lead to increased community awareness of the potential of nacre as a versatile, bioactive ceramic capable of improving bone tissue regeneration and will elicit increased research effort and innovation utilizing nacre.

Hyriopsis cumingii Hic52-A novel nacreous layer matrix protein with a collagen-like structure

Nacre is a product of a precisely regulated biomineralization process and a major contributor to the luster of pearls. Nacre is composed of calcium carbonate and an organic matrix of proteins that is secreted from mollusc mantle tissue and is exclusively associated with shell formation. In this study, hic52, a novel matrix protein gene from mantle of Hyriopsis cumingii, was cloned and functionally analyzed. The full-length cDNA of hic52 encoded 542 amino acids and contained a signal peptide of 18 amino acids. Excluding the signal peptide, the theoretical molecular mass of the polypeptide was 52.2kDa. The predicted isoelectric point was 10.37, indicating a basic shell protein. The amino acid sequence of hic52 featured high proportion of Gly (28.8%) and Gln (12.4%) residues. The predicted tertiary structure was characterized as having similarities to collagen I, alpha 1 and alpha 2 in the structure. The polypeptide sequence shared no homology with collagen. The hic52 expression pattern by quantitative real-time PCR and in situ hybridization exhibits at the dorsal epithelial cells of the mantle. Expression increased during the stages of pearl sac development. The data showed that hic52 is probably a framework shell protein that mediates and controls the nacreous biomineralization process.

Co-Option and De Novo Gene Evolution Underlie Molluscan Shell Diversity

Molluscs fabricate shells of incredible diversity and complexity by localized secretions from the dorsal epithelium of the mantle. Although distantly related molluscs express remarkably different secreted gene products, it remains unclear if the evolution of shell structure and pattern is underpinned by the differential co-option of conserved genes or the integration of lineage-specific genes into the mantle regulatory program. To address this, we compare the mantle transcriptomes of 11 bivalves and gastropods of varying relatedness. We find that each species, including four Pinctada (pearl oyster) species that diverged within the last 20 Ma, expresses a unique mantle secretome. Lineage- or species-specific genes comprise a large proportion of each species' mantle secretome. A majority of these secreted proteins have unique domain architectures that include repetitive, low complexity domains (RLCDs), which evolve rapidly, and have a proclivity to expand, contract and rearrange in the genome. There are also a large number of secretome genes expressed in the mantle that arose before the origin of gastropods and bivalves. Each species expresses a unique set of these more ancient genes consistent with their independent co-option into these mantle gene regulatory networks. From this analysis, we infer lineage-specific secretomes underlie shell diversity, and include both rapidly evolving RLCD-containing proteins, and the continual recruitment and loss of both ancient and recently evolved genes into the periphery of the regulatory network controlling gene expression in the mantle epithelium.

Characterization of allograft inflammatory factor-1 in Hyriopsis cumingii and its expression in response to immune challenge and pearl sac formation

The allograft inflammatory factor-1 (AIF-1) is one of the key factors associated with inflammatory response and immune defense. In the present study, we report the identification and characterization of AIF-1 from triangle sail mussel Hyriopsis cumingii (HcAIF-1). The full-length cDNA of HcAIF-1 consisted of a 5'-terminal untranslated region (UTR) of 80 bp, a 3'-UTR of 420 bp with a poly (A) tail, and an open reading frame of 444 bp encoding a polypeptide of 147 amino acids with two conserved EF-hand Ca²⁺-binding motifs. HcAIF-1 mRNA and protein were expressed in all examined tissues and showed higher mRNA expression levels were observed in immune tissues, especially hemocytes and mantle, and the highest protein expression level was in mantle. The expression level of HcAIF-1 mRNA was significantly upregulated in hemocytes 12-48 h after lipopolysaccharide challenge. After mantle tissue implantation, the expression level of this gene in pearl sac decreased significantly at 3-48 h (P < 0.01), and then was significantly upregulated at 96 h (P < 0.05) and recovered to the control level at 21-28 d. There was significant increase HcAIF-1 transcript abundance in hemocytes 96 h (P < 0.05) after mantle tissue implantation. The phagocytosis rate was significantly enhanced in hemocytes 3-24 h (P < 0.01) after the injection of recombinant HcAIF-1 protein. These findings suggest that HcAIF-1 is important in the underlying mechanism of the innate immune responses and pearl sac formation of H. cumingii.

Construction of a high-density genetic map and QTL mapping for pearl quality-related traits in Hyriopsis cumingii

A high-density genetic map is essential for quantitative trait locus (QTL) fine mapping. In this study, 4,508 effective single nucleotide polymorphism markers (detected using specific-locus amplified fragment sequencing) and 475 microsatellites were mapped to 19 linkage groups (LGs) using a family with 157 individuals. The map spanned 2,713 cM, with an average of 259 markers and 79 loci per LG and an average inter-marker distance of 1.81 cM. To identify QTLs for pearl quality traits, 26 putatively significant QTLs were detected for 10 traits, including, three for shell width, seven for body weight, two for shell weight, two for margin mantle weight, five for inner mantle weight, and seven for shell nacre colour. Among them, five QTLs associated with shell nacre colour were mapped to LG17 and explained 19.7% to 22.8% of the trait variation; this suggests that some important genes or loci determine shell nacre colour in LG17. The linkage map and mapped QTLs for shell nacre colour would be useful for improving the quality of Hyriopsis cumingii via marker-assisted selection.

A galectin from Hyriopsis cumingii involved in the innate immune response against to pathogenic microorganism and its expression profiling during pearl sac formation

Hyriopsis cumingii is the most important freshwater pearl mussel cultured in China. The operation for implantation is one necessary technical step for pearl culture. However, implantation-induced trauma results in a series of immune responses and can enable the invasion of pathogenic microbes. Lectin proteins are found widely in nature and play important roles in innate immunity. Galectins are members of the lectin superfamily and are characterized by one or several carbohydrate recognition domains (CRDs) that produce multiple sugar binding sites on the protein. Here we cloned and characterized the H. cumingii galectin gene HcGal1, which encodes a 312 amino acid galectin protein. The HcGal1 transcript was detected in all tested H. cumingii tissues and showed higher expression specifically in immune tissues. The significant upregulation of HcGal1 expression was observed after challenging the mussel with lipopolysaccharide or Gram-negative and Gram-positive bacteria. After implantation, significant downregulation of the HcGal1 transcript was noted in the mantle, hemocytes, and pearl sac in the acute-stress stage (0-24 h) and the stage of wound healing and pearl-sac formation (24 h-7 d). In addition, significant upregulation of HcGal1 expression was observed in the liver in the stage of wound healing and pearl-sac formation. In the pearl-secretion stage (7-35 d), the HcGal1 transcript levels returned to normal in all tested tissues. We also show that recombinantly expressed and purified HcGal1 can agglutinate some Gram-negative and Gram-positive bacteria. In addition, in vivo experiments showed that the recombinant protein HcGal1 could promote phagocytosis by hemocytes. Our data suggest that HcGal1 plays a role in innate immune responses involved in pathogen recognition and wound healing.

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

The Yesso scallop Patinopecten yessoensis displays polymorphism in shell colors, which is of great interest for the scallop industry. To identify genes involved in the shell coloration, in the present study, we investigate the transcriptome differences by Illumina digital gene expression (DGE) analysis in two extreme color phenotypes, Red and White. Illumina sequencing yields a total of 62,715,364 clean sequence reads, and more than 85% reads are mapped into our previously sequenced transcriptome. There are 25 significantly differentially expressed genes between Red and White scallops. EPR (Electron paramagnetic resonance) analysis has identified EPR spectra of pheomelanin and eumelanin in the red shells, but not in the white shells. Compared to the Red scallops, the White scallops have relatively higher mRNA expression in tyrosinase genes, but lower expression in other melanogensis-associated genes. Meantime, the relatively lower tyrosinase protein and decreased tyrosinase activity in White scallops are suggested to be associated with the lack of melanin in the white shells. Our findings highlight the functional roles of melanogensis-associated genes in the melanization process of scallop shells, and shed new lights on the transcriptional and post-transcriptional mechanisms in the regulation of tyrosinase activity during the process of melanin synthesis. The present results will assist our molecular understanding of melanin synthesis underlying shell color polymorphism in scallops, as well as other bivalves, and also help the color-based breeding in shellfish aquaculture.

Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization

An external skeleton is an essential part of the body plan of many animals and is thought to be one of the key factors that enabled the great expansion in animal diversity and disparity during the Cambrian explosion. Molluscs are considered ideal to study the evolution of biomineralization because of their diversity of highly complex, robust and patterned shells. The molluscan shell forms externally at the interface of animal and environment, and involves controlled deposition of calcium carbonate within a framework of macromolecules that are secreted from the dorsal mantle epithelium. Despite its deep conservation within Mollusca, the mantle is capable of producing an incredible diversity of shell patterns, and macro- and micro-architectures. Here we review recent developments within the field of molluscan biomineralization, focusing on the genes expressed in the mantle that encode secreted proteins. The so-called mantle secretome appears to regulate shell deposition and patterning and in some cases becomes part of the shell matrix. Recent transcriptomic and proteomic studies have revealed marked differences in the mantle secretomes of even closely-related molluscs; these typically exceed expected differences based on characteristics of the external shell. All mantle secretomes surveyed to date include novel genes encoding lineage-restricted proteins and unique combinations of co-opted ancient genes. A surprisingly large proportion of both ancient and novel secreted proteins containing simple repetitive motifs or domains that are often modular in construction. These repetitive low complexity domains (RLCDs) appear to further promote the evolvability of the mantle secretome, resulting in domain shuffling, expansion and loss. RLCD families further evolve via slippage and other mechanisms associated with repetitive sequences. As analogous types of secreted proteins are expressed in biomineralizing tissues in other animals, insights into the evolution of the genes underlying molluscan shell formation may be applied more broadly to understanding the evolution of metazoan biomineralization.

Comparative Transcriptome Analysis of the Pacific Oyster Crassostrea gigas Characterized by Shell Colors: Identification of Genetic Bases Potentially Involved in Pigmentation

BACKGROUND

Shell color polymorphisms of Mollusca have contributed to development of evolutionary biology and population genetics, while the genetic bases and molecular mechanisms underlying shell pigmentation are poorly understood. The Pacific oyster (Crassostrea gigas) is one of the most important farmed oysters worldwide. Through successive family selection, four shell color variants (white, golden, black and partially pigmented) of C. gigas have been developed. To elucidate the genetic mechanisms of shell coloration in C. gigas and facilitate the selection of elite oyster lines with desired coloration patterns, differentially expressed genes (DEGs) were identified among the four shell color variants by RNA-seq.

RESULTS

Digital gene expression generated over fifteen million reads per sample, producing expression data for 28,027 genes. A total number of 2,645 DEGs were identified from pair-wise comparisons, of which 432, 91, 43 and 39 genes specially were up-regulated in white, black, golden and partially pigmented shell of C. gigas, respectively. Three genes of Abca1, Abca3 and Abcb1 which belong to the ATP-binding cassette (ABC) transporters super-families were significantly associated with white shell formation. A tyrosinase transcript (CGI_10008737) represented consistent up-regulated pattern with golden coloration. We proposed that white shell variant of C. gigas could employ "endocytosis" to down-regulate notch level and to prevent shell pigmentation.

CONCLUSION

This study discovered some potential shell coloration genes and related molecular mechanisms by the RNA-seq, which would provide foundational information to further study on shell coloration and assist in selective breeding in C. gigas.

Mantle Branch-Specific RNA Sequences of Moon Scallop Amusium pleuronectes to Identify Shell Color-Associated Genes

Amusium pleuronectes (Linnaeus) that secretes red- and white-colored valves in two branches of mantle tissues is an excellent model for shell color research. High-throughput transcriptome sequencing and profiling were applied in this project to reveal the detailed molecular mechanism of this phenotype differentiation. In this study, 50,796,780 and 54,361,178 clean reads were generated from the left branch (secreting red valve, RS) and right branch (secreting white valve, WS) using the Illumina Hiseq 2000 platform. De novo assembly generated 149,375 and 176,652 unigenes with an average length of 764 bp and 698 bp in RS and WS, respectively. Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathway analysis indicated that the differentially expressed genes were involved in 228 signaling pathways, and 43 genes were significantly enriched (P<0.01). Nineteen of 20 differentially expressed vitellogenin genes showed significantly high expression in RS, which suggested that they probably played a crucial role in organic pigment assembly and transportation of the shell. Moreover, 687 crystal formation-related (or biomineralization-related) genes were detected in A. pleuronectes, among which 144 genes exhibited significant difference between the two branches. Those genes could be classified into shell matrix framework participants, crystal nucleation and growth-related elements, upstream regulation factors, Ca level regulators, and other classifications. We also identified putative SNP and SSR markers from these samples which provided the markers for genetic diversity analysis, genetic linkage, QTL analysis. These results provide insight into the complexity of shell color differentiation in A. pleuronectes so as valuable resources for further research.

A Novel Matrix Protein Hic31 from the Prismatic Layer of Hyriopsis Cumingii Displays a Collagen-Like Structure

In this study, we clone and characterize a novel matrix protein, hic31, from the mantle of Hyriopsis cumingii. The amino acid composition of hic31 consists of a high proportion of Glycine residues (26.67%). Tissue expression detection by RT-PCR indicates that hic31 is expressed specifically at the mantle edge. In situ hybridization results reveals strong signals from the dorsal epithelial cells of the outer fold at the mantle edge, and weak signals from inner epithelial cells of the same fold, indicating that hic31 is a prismatic-layer matrix protein. Although BLASTP results identify no shared homology with other shell-matrix proteins or any other known proteins, the hic31 tertiary structure is similar to that of collagen I, alpha 1 and alpha 2. It has been well proved that collagen forms the basic organic frameworks in way of collagen fibrils and minerals present within or outside of these fibrils. Therefore, hic31 might be a framework-matrix protein involved in the prismatic-layer biomineralization. Besides, the gene expression of hic31 increase in the early stages of pearl sac development, indicating that hic31 may play important roles in biomineralization of the pearl prismatic layer.

Characterization of the mantle transcriptome of yesso scallop (Patinopecten yessoensis): identification of genes potentially involved in biomineralization and pigmentation

The Yesso scallop Patinopecten yessoensis is an economically important marine bivalve species in aquaculture and fishery in Asian countries. However, limited genomic resources are available for this scallop, which hampers investigations into molecular mechanisms underlying their unique biological characteristics, such as shell formation and pigmentation. Mantle is the special tissue of P. yessoensis that secretes biomineralization proteins inducing shell deposition as well as pigmentation on the shells. However, a current deficiency of transcriptome information limits insight into mechanisms of shell formation and pigmentation in this species. In this study, the transcriptome of the mantle of P. yessoensis was deeply sequenced and characterized using Illumina RNA-seq technology. A total of 86,521 unique transcripts are assembled from 55,884,122 reads that passed quality filters, and annotated, using Gene Ontology classification. A total of 259 pathways are identified in the mantle transcriptome, including the calcium signaling and melanogenesis pathways. A total of 237 unigenes that are homologous to 102 reported biomineralization genes are identified, and 121 unigenes that are homologous to 93 known proteins related to melanin biosynthesis are found. Twenty-three annotated unigenes, which are mainly homologous to calmodulin and related proteins, Ca2+/calmodulin-dependent protein kinase, adenylate/guanylate cyclase, and tyrosinase family are potentially involved in both biomineralization and melanin biosynthesis. It is suggested that these genes are probably not limited in function to induce shell deposition by calcium metabolism, but may also be involved in pigmentation of the shells of the scallop. This potentially supports the idea that there might be a link between calcium metabolism and melanin biosynthesis, which was previously found in vertebrates. The findings presented here will notably advance the understanding of the sophisticated processes of shell formation as well as shell pigmentation in P. yessoensis and other bivalve species, and also provide new evidence on gene expression for the understanding of pigmentation and biomineralization not only in invertebrates but also probably in vertebrates.

Transcriptomics of developing embryos and organs: A raising tool for evo-devo

Comparative transcriptomics has become an important tool for revisiting many evo-devo questions and exploring new ones, and its importance is likely to increase in the near future, partly because RNA-seq data open many new possibilities. The aim of this opinion piece is twofold. In the first section, we discuss the particularities of transcriptomic studies in evo-devo, focusing mainly on RNA-seq data. The preliminary processing steps (getting coding sequences as well as expression levels) are challenging, because many studied species do not have a sequenced genome. The next step (interpreting expression differences) is also challenging, due to several issues with interpreting expression levels in complex tissues, managing developmental stages and species heterochronies, and the problem of conceptualizing expression differences. In the second section, we discuss some past and possible future applications of transcriptomic approaches (using microarray or RNA-seq) to three major themes in evo-devo: the evolution of the developmental toolkit, the genetic and developmental basis for phenotypic changes, and the general rules of the evolution of development. We believe that conceptual and technical tools are necessary in order to fully exploit the richness of multispecies transcriptomic time-series data.

Evolution of the tyrosinase gene family in bivalve molluscs: independent expansion of the mantle gene repertoire

Tyrosinase is a copper-containing enzyme that mediates the hydroxylation of monophenols and oxidation of o-diphenols to o-quinones. This enzyme is involved in a variety of biological processes, including pigment production, innate immunity, wound healing, and exoskeleton fabrication and hardening (e.g. arthropod skeleton and mollusc shell). Here we show that the tyrosinase gene family has undergone large expansions in pearl oysters (Pinctada spp.) and the Pacific oyster (Crassostrea gigas). Phylogenetic analysis reveals that pearl oysters possess at least four tyrosinase genes that are not present in the Pacific oyster. Likewise, C. gigas has multiple tyrosinase genes that are not orthologous to the Pinctada genes, indicating that this gene family has expanded independently in these bivalve lineages. Many of the tyrosinase genes in these bivalves are expressed at relatively high levels in the mantle, the organ responsible for shell fabrication. Detailed comparisons of tyrosinase gene expression in different regions of the mantle in two closely related pearl oysters, P. maxima and P. margaritifera, reveals that recently evolved orthologous tyrosinase genes can have markedly different expression profiles. The expansion of tyrosinase genes in these oysters and their co-option into the mantle's gene regulatory network is consistent with mollusc shell formation being underpinned by a rapidly evolving transcriptome.

Molecular biology approaches in bioadhesion research

The use of molecular biology tools in the field of bioadhesion is still in its infancy. For new research groups who are considering taking a molecular approach, the techniques presented here are essential to unravelling the sequence of a gene, its expression and its biological function. Here we provide an outline for addressing adhesion-related genes in diverse organisms. We show how to gradually narrow down the number of candidate transcripts that are involved in adhesion by (1) generating a transcriptome and a differentially expressed cDNA list enriched for adhesion-related transcripts, (2) setting up a BLAST search facility, (3) perform an in situ hybridization screen, and (4) functional analyses of selected genes by using RNA interference knock-down. Furthermore, latest developments in genome-editing are presented as new tools to study gene function. By using this iterative multi-technologies approach, the identification, isolation, expression and function of adhesion-related genes can be studied in most organisms. These tools will improve our understanding of the diversity of molecules used for adhesion in different organisms and these findings will help to develop innovative bio-inspired adhesives.

Transcriptomic analysis of candidate osmoregulatory genes in the eastern oyster Crassostrea virginica

BACKGROUND

The eastern oyster, Crassostrea virginica, is a euryhaline species that can thrive across a wide range of salinities (5-35). As with all estuarine species, individual oysters must be able to regulate their osmotic balance in response to constant temporal variation in salinity. At the population level, recurrent viability selection may be an additional mechanism shaping adaptive osmoregulatory phenotypes at the margins of oyster salinity tolerance. To identify candidate genes for osmoregulation, we sequenced, assembled, and annotated the transcriptome of wild juvenile eastern oysters from 'high' and 'low' salinity regimes. Annotations and candidates were mostly based on the Pacific oyster (Crassostrea gigas) genome sequence so osmoregulatory relevance in C. virginica was explored by testing functional enrichment of genes showing spatially discrete patterns of expression and by quantifying coding sequence divergence.

RESULTS

The assembly of sequence reads and permissive clustering of potentially oversplit alleles resulted in 98,729 reftigs (contigs and singletons). Of these, 50,736 were annotated with 9,307 belonging to a set of candidate osmoregulatory genes identified from the C. gigas genome. A total of 218,777 SNPs (0.0185 SNPs/bp) were identified in annotated reftigs of C. virginica. Amino acid divergence between translations of C. virginica annotated reftigs and C. gigas coding sequence averaged 23.2 % with an average dN/dS ratio of 0.074, suggesting purifying selection on protein sequences. The high and low salinity source oysters each expressed a subset of genes unique to that group, and the functions for these annotated genes were consistent with known molecular mechanisms for osmotic regulation in molluscs.

CONCLUSIONS

Most of the osmoregulatory gene candidates experimentally identified in C. gigas are present in this C. virginica transcriptome. In general these congeners show coding sequence divergence too high to make the C. gigas genome a useful reference for C. virginica bioinformatics. However, strong purifying selection is characteristic of the osmoregulatory candidates so functional annotations are likely to correspond. An initial examination of C. virginica presence/absence expression patterns across the salinity gradient in a single estuary suggests that many of these candidates have expression patterns that co-vary with salinity, consistent with osmoregulatory function in C. virginica.

Investigating protein function in biomineralized tissues using molecular biology techniques

We describe modern molecular biology methods currently used in the study of biomineralization. We focus our descriptions on two areas of biomineralization research in which these methods have been particularly powerful. The first area is the use of modern molecular methods to identify and characterize the so-called occluded matrix proteins present in mineralized tissues. More specifically, we describe the use of RNA-seq and the next generation of DNA sequencers and the use of direct protein sequencing and mass spectrometers as ways of identifying proteins present in mineralized tissues. The second area is the use of molecular methods to examine the function of proteins in biomineralization. RNA interference (RNAi), morpholino antisense, and other methods are described and discussed as ways of elucidating protein function.

Identification of housekeeping genes suitable for gene expression analysis in the pearl mussel, Hyriopsis cumingii, during biomineralization

Quantitative real-time polymerase chain reaction is a sensitive technique for quantifying gene expression levels. One or more appropriate reference genes must be selected to accurately compare mRNA transcript levels across different samples and tissues. The freshwater pearl, Hyriopsis cumingii (Lea), is an important economic species cultured in China. To date, no reference genes for gene expression analysis in this species have been validated. This study aimed to compare the relative expression of seven housekeeping genes across different tissue types and in the mantle or pearl sac during three biomineralization processes: seasonal shell growth, shell healing and pearl-sac formation in H. cumingii. Three programs evaluated the expression stabilities of the seven genes: BestKeeper, geNorm and NormFinder. The beta actin gene (ACTB), commonly used as a housekeeping gene in many studies, was the least stable. The expressions of Ubiquitin (Ubi) and Ribosomal protein L18 (Rpl18) and Elongation factor 1-alpha (EF1α) were more stable than the remaining four genes. Therefore, we suggest that Ubi, Rpl18 and EF1α are suitable reference genes. The three selected reference genes are expected to facilitate analysis of gene expressions during shell or pearl formation in H. cumingii.

Cloning, differential tissue expression of a novel hcApo gene, and its correlation with total carotenoid content in purple and white inner-shell color pearl mussel Hyriopsis cumingii

As a molecular carrier and storage protein, apolipoprotein (Apo) mediates the intracellular uptake of lipids, proteins, vitamins and carotenoids. In this study, we identified a novel Apo gene, designated hcApo, from the freshwater pearl mussel Hyriopsis cumingii. The complete hcApo cDNA consists of 4104 nucleotides with an open reading frame encoding 1155 amino acid residues. The hcApo protein contains a conserved lipoprotein N-terminal domain (LPD-N) that is a characteristic of the large lipid transfer protein (LLTP) superfamily. The hcApo mRNA is constitutively expressed in a wide range of tissues with the highest expression level in the liver. Moreover, differential expression analysis revealed that the hcApo gene is more highly expressed in the liver, kidney, mantle and gill of purple line mussels compared to white line mussels. In situ hybridization investigations of the precise expression site of hcApo mRNA in the mantle showed that hcApo mRNA is specifically expressed in the outer epithelial cells of the middle fold and the inner epithelial cells of the outer fold of the mantle, as well as throughout the outer epithelium of the outer fold and ventral mantle. Another very important finding is that significantly positive correlation existed between the hcApo gene expression level and the total carotenoid content in purple line mussels. These findings may provide a better understanding of the roles of hcApo in the molecular mechanisms of shell formation and coloring of H. cumingii.