Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes

Cloning, characterization, and spatio-temporal expression patterns of HdhSPARC and its responses to multiple stressors

SPARC is an extracellular Ca2+-binding, secreted glycoprotein that plays a dynamic role in the growth and development of organisms. This study aimed to describe the isolation, characterization, and expression analysis of HdhSPARC in Pacific abalone (Haliotis discus hannai) to infer its potential functional role. The isolated HdhSPARC was 1633 bp long, encoding a polypeptide of 284 amino acid residues. Structurally, the SPARC protein in abalone is comprised of three biological domains. However, the structure of this protein varied between vertebrates and invertebrates, as suggested by their distinct clustering patterns in phylogenetic analysis. In early development, HdhSPARC was variably expressed, and higher expression was found in veliger larvae. Moreover, HdhSPARC was highly expressed in juvenile abalone with rapid growth compared to their slower-growing counterparts. Among the testicular development stages, the growth stage exhibited higher HdhSPARC expression. HdhSPARC was also upregulated during muscle remodeling and shell biomineralization, as well as in response to different stressors such as heat shock, LPS, and H2O2 exposure. However, this gene was downregulated in Cd-exposed abalone. The present study first comprehensively characterized the HdhSPARC gene, and its spatio-temporal expressions were analyzed along with its responses to various stressors.

Evolution of the expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

Spatial-temporal expression analysis of lineage-restricted shell matrix proteins reveals shell field regionalization and distinct cell populations in the slipper snail Crepidula atrasolea

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea , to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata -restricted SMPs contradicts the generally-held notion that an animal’s biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea . We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

RNA-Seq comparative study reveals molecular effectors linked to the resistance of Pinna nobilis to Haplosporidium pinnae parasite

With the intensification of maritime traffic, recently emerged infectious diseases have become major drivers in the decline and extinction of species. Since 2016, mass mortality events have decimated the endemic Mediterranean Sea bivalve Pinna nobilis, affecting ca. 100% of individuals. These events have largely been driven by Haplosporidium pinnae's infection, an invasive species which was likely introduced by shipping. While monitoring wild populations of P. nobilis, we observed individuals that survived such a mass mortality event during the summer of 2018 (France). We considered these individuals resistant, as they did not show any symptoms of the disease, while the rest of the population in the area was devastated. Furthermore, the parasite was not detected when we conducted a PCR amplification of a species-specific fragment of the small subunit ribosomal DNA. In parallel, the transcriptomic analysis showed evidence of some parasite RNA indicating that the resistant individuals had been exposed to the parasite without proliferating. To understand the underlying mechanisms of resistance in these individuals, we compared their gene expression with that of susceptible individuals. We performed de novo transcriptome assembly and annotated the expressed genes. A comparison of the transcriptomes in resistant and susceptible individuals highlighted a gene expression signature of the resistant phenotype. We found significant differential expressions of genes involved in immunity and cell architecture. This data provides the first insights into how individuals escape the pathogenicity associated with infection.

Genetic differentiation of a critically endangered population of the limpet Patella candei candei d’Orbigny, 1840, in the Canary Islands

The adoption of measures to protect the viability of threatened populations should be supported by empirical data identifying appropriate conservation units and management strategies. The global population of the majorera limpet, P. candei candei d’Orbigny, 1840, is restricted to the Macaronesian islands in the NE Atlantic, including near-to-extinct and healthy populations in Fuerteventura and Selvagens, respectively. The taxonomic position, genetic diversity and intra- and interspecific relationships of these populations are unclear, which is hindering the implementation of a recovery plan for the overexploited majorera limpet on Fuerteventura. In this study, ddRAD-based genome scanning was used to overcome the limitations of mitochondrial DNA-based analysis. As a result, P. candei candei was genetically differentiated from the closely related P. candei crenata for the first time. Moreover, genetic differentiation was detected between P. candei candei samples from Selvagens and Fuerteventura, indicating that translocations from the healthy Selvagens source population are inadvisable. In conclusion, the majorera limpet requires population-specific management focused on the preservation of exceptional genetic diversity with which to face future environmental challenges.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Characterization of Thrombospondin Type 1 Repeat in Haliotis diversicolor and Its Possible Role in Osteoinduction

Thrombospondin repeats (TSR) are important peptide domains present in the sequences of many extracellular and transmembrane proteins with which a variety of ligands interact. In this study, we characterized HdTSR domains in the ADAMTS3 protein of Thai abalone, Haliotis diversicolor, based on the transcriptomic analysis of its mantle tissues. PCR amplification and localization studies demonstrated the existence of HdTSR transcript and protein in H. diversicolor tissues, particularly in both the inner and outer mantle epithelial folds. We, therefore, generated a short recombinant protein, termed HdTSR1/2, based on the existence of the WxxWxxW or WxxxxW motif (which binds to TGF-β, a known signaling in bone formation/repair) in HdTSR1 and HdTSR2 sequences and used it to test the osteoinduction function in the pre-osteoblastic cell line, MC3T3-E1. This recombinant protein demonstrated the ability to induce the differentiation of MC3T3-E1 cells by the concentration- and time-dependent upregulation of many known osteogenic markers, including RUNX2, COL1A1, OCN, and OPN. We also demonstrated the upregulation of the SMAD2 gene after cell treatment with HdTSR1/2 proteinindicating its possible interaction through TGF-β, which thus activates its downstream signaling cascade and triggers the biomineralization process in the differentiated osteoblastic cells. Together, HdTSR domains existed in an extracellular ADAMTS3 protein in the mantle epithelium of H. diversicolor and played a role in osteoinduction as similar to the other nacreous proteins, opening up its possibility to be developed as an inducing agent of bone repair.

Molecular Characterization of Carbonic Anhydrase II (CA II) and Its Potential Involvement in Regulating Shell Formation in the Pacific Abalone, Haliotis discus hannai

Carbonic anhydrases (CAs) are a family of metalloenzymes that can catalyze the reversible interconversion of CO₂/HCO₃^–, ubiquitously present in both prokaryotes and eukaryotes. In the present study, a CA II (designated as HdhCA II) was sequenced and characterized from the mantle tissue of the Pacific abalone. The complete sequence of HdhCA II was 1,169 bp, encoding a polypeptide of 349 amino acids with a NH₂-terminal signal peptide and a CA architectural domain. The predicted protein shared 98.57% and 68.59% sequence identities with CA II of Haliotis gigantea and Haliotis tuberculata, respectively. Two putative N-linked glycosylation motifs and two cysteine residues could potentially form intramolecular disulfide bond present in HdhCA II. The phylogenetic analysis indicated that HdhCA II was placed in a gastropod clade and robustly clustered with CA II of H. gigantea and H. tuberculata. The highest level of HdhCA II mRNA expression was detected in the shell forming mantle tissue. During ontogenesis, the mRNA of HdhCA II was detected in all stages, with larval shell formation stage showing the highest expression level. The in situ hybridization results detected the HdhCA II mRNA expression in the epithelial cells of the dorsal mantle pallial, an area known to express genes involved in the formation of a nacreous layer in the shell. This is the first report of HdhCA II in the Pacific abalone, and the results of this study indicate that this gene might play a role in the shell formation of abalone.

Characterization of Insulin-Like Growth Factor Binding Protein 7 (Igfbp7) and Its Potential Involvement in Shell Formation and Metamorphosis of Pacific Abalone, Haliotis discus hannai

Insulin-like growth factor binding proteins (IGFBPs) are secreted proteins that play an important role in IGF regulation of growth and development of vertebrate and invertebrates. In this study, the IGFBP7 gene was cloned and characterized from mantle tissues of H. discus hannai, and designated as Hdh IGFBP7. The full-length cDNA sequence transcribed from the Hdh IGFBP7 gene was 1519-bp long with an open reading frame of 720-bp corresponding to a putative polypeptide of 239 amino acids. The molecular mass of its mature protein was approximately 23.44 KDa with an estimated isoelectric point (pI) of 5.35, and it shared significant homology with IGFBP7 gene of H. madaka. Hdh IGFBP7 has a characteristic IGFBP N-terminal domain (22–89 aa), a kazal-type serine proteinase inhibitor domain (77–128), and an immunoglobulin-like C2 domain (144–223). Furthermore, twelve cysteine residues and a signature motif of IGFBPs (XCGCCXXC) were found in its N-terminal domain. Phylogenetic analysis revealed that Hdh IGFBP7 was aligned with IGFBP7 of H. madaka. Tissue distribution analysis showed that the mRNA of Hdh IGFBP7 was expressed in all examined tissues, with the highest expression level observed in the mantle and gill tissues. The expression level of Hdh IGFBP7 mRNA was relatively higher at the juvenile stage during its metamorphosis period. In situ hybridization showed that Hdh IGFBP7 transcript was expressed in epithelial cells of the dorsal mantle pallial and mucus cells of the branchial epithelium in gill. These results provide basic information for future studies on the role of IGFBP7 in IGF regulation of shell growth, development and metamorphosis of abalone.

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.

The calcitonin-like system is an ancient regulatory system of biomineralization

Biomineralization is the process by which living organisms acquired the capacity to accumulate minerals in tissues. Shells are the biomineralized exoskeleton of marine molluscs produced by the mantle but factors that regulate mantle shell building are still enigmatic. This study sought to identify candidate regulatory factors of molluscan shell mineralization and targeted family B G-protein coupled receptors (GPCRs) and ligands that include calcium regulatory factors in vertebrates, such as calcitonin (CALC). In molluscs, CALC receptor (CALCR) number was variable and arose through lineage and species-specific duplications. The Mediterranean mussel (Mytilus galloprovincialis) mantle transcriptome expresses six CALCR-like and two CALC-precursors encoding four putative mature peptides. Mussel CALCR-like are activated in vitro by vertebrate CALC but only receptor CALCRIIc is activated by the mussel CALCIIa peptide (EC50 = 2.6 ×10-5 M). Ex-vivo incubations of mantle edge tissue and mantle cells with CALCIIa revealed they accumulated significantly more calcium than untreated tissue and cells. Mussel CALCIIa also significantly decreased mantle acid phosphatase activity, which is associated with shell remodelling. Our data indicate the CALC-like system as candidate regulatory factors of shell mineralization. The identification of the CALC system from molluscs to vertebrates suggests it is an ancient and conserved calcium regulatory system of mineralization.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

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.

Evolution and diversity of alpha-carbonic anhydrases in the mantle of the Mediterranean mussel (Mytilus galloprovincialis)

The α-carbonic anhydrases (α-CAs) are a large and ancient group of metazoan-specific enzymes. They generate bicarbonate from metabolic carbon dioxide and through calcium carbonate crystal formation play a key role in the regulation of mineralized structures. To better understand how α-CAs contribute to shell mineralization in the marine Mediterranean mussel (Mytilus galloprovincialis) we characterized them in the mantle. Phylogenetic analysis revealed that mollusc α-CA evolution was affected by lineage and species-specific events. Ten α-CAs were found in the Mediterranean mussel mantle and the most abundant form was named, MgNACR, as it grouped with oyster nacreins (NACR). Exposure of the Mediterranean mussel to reduced water salinity (18 vs 37 ppt), caused a significant reduction (p < 0.05) in mantle esterase activity and MgNACR transcript abundance (p < 0.05). Protonograms revealed multiple proteins in the mantle with α-CA hydratase activity and mapped to a protein with a similar size to that deduced for monomeric MgNACR. Our data indicate that MgNACR is a major α-CA enzyme in mantle and that by homology with oyster nacreins likely regulates mussel shell production. We propose that species-dependent α-CA evolution may contribute to explain the diversity of bivalve shell structures and their vulnerability to environmental changes.

Insights into the Evolution of Shells and Love Darts of Land Snails Revealed from Their Matrix Proteins

Over the past decade, many skeletal matrix proteins that are possibly related to calcification have been reported in various calcifying animals. Molluscs are among the most diverse calcifying animals and some gastropods have adapted to terrestrial ecological niches. Although many shell matrix proteins (SMPs) have already been reported in molluscs, most reports have focused on marine molluscs, and the SMPs of terrestrial snails remain unclear. In addition, some terrestrial stylommatophoran snails have evolved an additional unique calcified character, called a "love dart," used for mating behavior. We identified 54 SMPs in the terrestrial snail Euhadra quaesita, and found that they contain specific domains that are widely conserved in molluscan SMPs. However, our results also suggest that some of them possibly have evolved independently by domain shuffling, domain recruitment, or gene co-option. We then identified four dart matrix proteins, and found that two of them are the same proteins as those identified as SMPs. Our results suggest that some dart matrix proteins possibly have evolved by independent gene co-option from SMPs during dart evolution events. These results provide a new perspective on the evolution of SMPs and "love darts" in land snails.

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.

In-depth proteomic analyses of Haliotis laevigata (greenlip abalone) nacre and prismatic organic shell matrix

Background

The shells of various Haliotis species have served as models of invertebrate biomineralization and physical shell properties for more than 20 years. A focus of this research has been the nacreous inner layer of the shell with its conspicuous arrangement of aragonite platelets, resembling in cross-section a brick-and-mortar wall. In comparison, the outer, less stable, calcitic prismatic layer has received much less attention. One of the first molluscan shell proteins to be characterized at the molecular level was Lustrin A, a component of the nacreous organic matrix of Haliotis rufescens. This was soon followed by the C-type lectin perlucin and the growth factor-binding perlustrin, both isolated from H. laevigata nacre, and the crystal growth-modulating AP7 and AP24, isolated from H. rufescens nacre. Mass spectrometry-based proteomics was subsequently applied to to Haliotis biomineralization research with the analysis of the H. asinina shell matrix and yielded 14 different shell-associated proteins. That study was the most comprehensive for a Haliotis species to date.

Methods

The shell proteomes of nacre and prismatic layer of the marine gastropod Haliotis laevigata were analyzed combining mass spectrometry-based proteomics and next generation sequencing.

Results

We identified 297 proteins from the nacreous shell layer and 350 proteins from the prismatic shell layer from the green lip abalone H. laevigata. Considering the overlap between the two sets we identified a total of 448 proteins. Fifty-one nacre proteins and 43 prismatic layer proteins were defined as major proteins based on their abundance at more than 0.2% of the total. The remaining proteins occurred at low abundance and may not play any significant role in shell fabrication. The overlap of major proteins between the two shell layers was 17, amounting to a total of 77 major proteins.

Conclusions

The H. laevigata shell proteome shares moderate sequence similarity at the protein level with other gastropod, bivalve and more distantly related invertebrate biomineralising proteomes. Features conserved in H. laevigata and other molluscan shell proteomes include short repetitive sequences of low complexity predicted to lack intrinsic three-dimensional structure, and domains such as tyrosinase, chitin-binding, and carbonic anhydrase. This catalogue of H. laevigata shell proteins represents the most comprehensive for a haliotid and should support future efforts to elucidate the molecular mechanisms of shell assembly.

Electronic supplementary material

The online version of this article (10.1186/s12953-018-0139-3) contains supplementary material, which is available to authorized users.

Proteomics Studies on the three Larval Stages of Development and Metamorphosis of Babylonia areolata

The ivory shell, Babylonia areolata, is a commercially important aquaculture species in the southeast coast of mainland China. The middle veliger stage, later veliger stage, and juvenile stage are distinct larval stages in B. areolata development. In this study, we used label-free quantification proteomics analysis of the three developmental stages of B. areolata. We identified a total of 5,583 proteins, of which 1,419 proteins expression level showed significant differential expression. The results of gene ontology enrichment analysis showed that the number of proteins involved in metabolic and cellular processes were the most abundant. Those proteins mostly had functions such as binding, catalytic activity and transporter activity. The results of Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that the number of proteins involved in the ribosome, carbon metabolism, and lysosome pathways were the most abundant, indicating that protein synthesis and the immune response were active during the three stages of development. This is the first study to use proteomics and real-time PCR to study the early developmental stages of B. areolata, which could provide relevant data on gastropod development. Our results provide insights into the novel aspects of protein function in shell formation, body torsion, changes in feeding habits, attachment and metamorphosis, immune-related activities in B. areolata larvae.

Functional Involvement of Carbonic Anhydrase in the Lysosomal Response to Cadmium Exposure in Mytilus galloprovincialis Digestive Gland

Carbonic anhydrase (CA) is a ubiquitous metalloenzyme, whose functions in animals span from respiration to pH homeostasis, electrolyte transport, calcification, and biosynthetic reactions. CA is sensitive to trace metals in a number of species. In mussels, a previous study demonstrated CA activity and protein expression to be enhanced in digestive gland by cadmium exposure. The aim of the present work was to investigate the functional meaning, if any, of this response. To this end the study addressed the possible involvement of CA in the lysosomal system response of digestive gland cells to metal exposure. The in vivo exposure to acetazolamide, specific CA inhibitor, significantly inhibited the acidification of the lysosomal compartment in the digestive gland cells charged with the acidotropic probe LysoSensor Green D-189, demonstrating in vivo the physiological contribution of CA to the acidification of the lysosomes. Under CdCl₂ exposure, CA activity significantly increased in parallel to the increase of the fluorescence of LysoSensor Green charged cells, which is in turn indicative of proliferation and/or increase in size of lysosomes. Acetazolamide exposure was able to completely inhibit the cadmium induced Lysosensor fluorescence increase in digestive gland cells. In conclusion, our results demonstrated the functional role of CA in the lysosomal acidification of Mytilus galloprovincialis digestive gland and its involvement in the lysosomal activation following cadmium exposure. CA induction could physiologically respond to a prolonged increased requirement of H⁺ for supporting lysosomal acidification during lysosomal activation.

Matricellular Proteins: Functional Insights From Non-mammalian Animal Models

The extracellular matrix (ECM) has central roles in tissue integrity and remodeling throughout the life span of animals. While collagens are the most abundant structural components of ECM in most tissues, tissue-specific molecular complexity is contributed by ECM glycoproteins. The matricellular glycoproteins are categorized primarily according to functional criteria and represented predominantly by the thrombospondin, tenascin, SPARC/osteonectin, and CCN families. These proteins do not self-assemble into ECM fibrils; nevertheless, they shape ECM properties through interactions with structural ECM proteins, growth factors, and cells. Matricellular proteins also promote cell migration or morphological changes through adhesion-modulating or counter-adhesive actions on cell-ECM adhesions, intracellular signaling, and the actin cytoskeleton. Typically, matricellular proteins are most highly expressed during embryonic development. In adult tissues, expression is more limited unless activated by cues for dynamic tissue remodeling and cell motility, such as occur during inflammatory response and wound repair. Many insights in the complex roles of matricellular proteins have been obtained from studies of gene knockout mice. However, with the exception of chordate-specific tenascins, these are highly conserved proteins that are encoded in many animal phyla. This review will consider the increasing body of research on matricellular proteins in nonmammalian animal models. These models provide better access to the very earliest stages of embryonic development and opportunities to study biological processes such as limb and organ regeneration. In aggregate, this research is expanding concepts of the functions and mechanisms of action of matricellular proteins.

De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea

Crassadoma gigantea is an important commercial marine bivalve species in Baja California and Mexico. In this study, we have applied RNA-Seq technology to profile the transcriptome of the C. gigantea for the first time. A total of 80,832,518 raw reads were produced from a Illumina HiSeq4000 platform, and 77,306,198 (95.64%) clean reads were generated after trimming the adaptor sequences. The transcriptome assembled into 158,855 transcripts with an N50 size of 1995 bp and an average size of 1008 bp. A number of DNA repair related genes, such as MSH3, EGF, TGF, IGF, FGF, encoding different groups of growth factors were found in the transcriptome data of C. gigantean. In addition, immune related genes Toll-like receptor (TLR) including TLR1, TLR2, TLR3, TLR4, TLR5, TLR6, TLR7, TLR8, and TLR9 was also observed in C. gigantean. A set of 12 polymorphic microsatellite loci were firstly developed and characterized in C. gigantea. The results show that the number of alleles and expected heterozygosity ranged from 3 to 9 and from 0.254 to 0.820, respectively. The average polymorphic information content was 0.790. These microsatellite loci will facilitate future studies of population structure and conservation genetics in this species.

Treatment of Patellar Lower Pole Fracture with Modified Titanium Cable Tension Band Plus Patellar Tibial Tunnel Steel "8" Reduction Band

Purpose: To determine the efficacy of modified titanium tension band plus patellar tendon tunnel steel 8 "reduction band" versus titanium cable tension band fixation for the treatment of patellar lower pole fracture. Materials and Methods: 58 patients with lower patella fracture were enrolled in this study, including 30 patients treated with modified titanium cable tension band plus patellar tibial tunnel wire "8" tension band internal fixation (modified group), and 28 patients with titanium cable tension band fixation. All patients were followed up for 9∼15 months with an average of 11.6 months. Results: Knee flexion was significantly improved in the modified group than in the titanium cable tension band group (111.33 ± 13 degrees versus 98.21 ± 21.70 degrees, P = 0.004). The fracture healing time showed no significant difference. At the end of the follow-up, the improvement excellent rate was 93.33% in the modified group, and 82.14% in the titanium cable tension band group. Titanium cable tension band internal fixation loosening was found in 2 cases, including 1 case of treatment by two surgeries without loose internal fixation. Conclusions: The modified titanium cable tension band with "8" tension band fixation showed better efficacy for lower patella fractures than titanium cable tension band fixation.

Toxic Effects and Molecular Mechanism of Different Types of Silver Nanoparticles to the Aquatic Crustacean Daphnia magna

Silver nanoparticles (AgNPs) have been assessed to have a high exposure risk for humans and aquatic organisms. Toxicity varies considerably between different types of AgNPs. This study aimed to investigate the toxic effects of AgNPs with different particle sizes (40 and 110 nm) and different surface coatings (sodium citrate and polyvinylpyrrolidone, PVP) on Daphnia magna and their mechanisms of action. The results revealed that the citrate-coated AgNPs were more toxic than PVP-coated AgNPs and that the 40 nm AgNPs were more toxic than the 110 nm AgNPs. Transcriptome analysis further revealed that the toxic effects of AgNPs on D. magna were related to the mechanisms of ion binding and several metabolic pathways, such as the "RNA polymerase" pathway and the "protein digestion and absorption" pathway. Moreover, the principal component analysis (PAC) results found that surface coating was the major factor that determines the toxicities compared to particle size. These results could help us better understand the possible mechanism of AgNP toxicity in aquatic invertebrates at the transcriptome level and establish an important foundation for revealing the broad impacts of nanoparticles on aquatic environments.

Hemolymph and gill carbonic anhydrase are more sensitive to aquatic contamination than mantle carbonic anhydrase in the mangrove oyster Crassostrea rhizophorae

Carbonic anhydrase (CA) is a ubiquitous metalloenzyme of great importance in several physiological processes. Due to its physiological importance and sensitivity to various pollutants, CA activity has been used as biomarker of aquatic contamination. Considering that in bivalves the sensitivity of CA to pollutants seems to be tissue-specific, we proposed here to analyze CA activity of hemolymph, gill and mantle of Crassostrea rhizophorae collected in two tropical Brazilian estuaries with different levels of anthropogenic impact, in dry and rainy season. We found increased carbonic anhydrase activity in hemolymph, gill and mantle of oysters collected in the Paraíba Estuary (a site of high anthropogenic impact) when compared to oysters from Mamanguape Estuary (inserted in an area of environmental preservation), especially in the rainy season. CA of hemolymph and gill were more sensitive than mantle CA to aquatic contamination. This study enhances the suitability of carbonic anhydrase activity for field biomarker applications with bivalves and brings new and relevant information on hemolymph carbonic anhydrase activity as biomarker of aquatic contamination.

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.

De novo assembly, characterization and annotation for the transcriptome of Sarcocheilichthys sinensis

The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

RESULTS

Next-generation sequencing was used to retrieve transcriptomes of representatives of seven out of the eight recent class-level taxa of mollusks. Similarity searches, phylogenetic inferences, and a detailed manual curation were used to identify and confirm the orthology of numerous molluscan Hox and ParaHox genes, which resulted in a comprehensive catalog that highlights the evolution of these genes in Mollusca and other metazoans. The identification of a specific molluscan motif in the Hox paralog group 5 and a lophotrochozoan ParaHox motif in the Gsx gene is described. Functional analyses using KEGG and GO tools enabled a detailed description of key developmental genes expressed in important pathways such as Hedgehog, Wnt, and Notch during development of the respective species. The KEGG analysis revealed Wnt8, Wnt11, and Wnt16 as Wnt genes hitherto not reported for mollusks, thereby enlarging the known Wnt complement of the phylum. In addition, novel Hedgehog (Hh)-related genes were identified in the gastropod Lottia cf. kogamogai, demonstrating a more complex gene content in this species than in other mollusks.

CONCLUSIONS

The use of de novo transcriptome assembly and well-designed in silico protocols proved to be a robust approach for surveying and mining large sequence data in a wide range of non-model mollusks. The data presented herein constitute only a small fraction of the information retrieved from the analysed molluscan transcriptomes, which can be promptly employed in the identification of novel genes and gene families, phylogenetic inferences, and other studies using molecular tools. As such, our study provides an important framework for understanding some of the underlying molecular mechanisms involved in molluscan body plan diversification and hints towards functions of key developmental genes in molluscan morphogenesis.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

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.

Characterization of the mantle transcriptome in bivalves: Pecten maximus, Mytilus edulis and Crassostrea gigas

The calcareous shells secreted by bivalve molluscs display diverse and species specific structural compositions, which indicates possible divergent biomineralization processes. Thus, studying multiple mollusc species will provide a more comprehensive understanding of shell formation. Here, the transcriptomes of the mantle tissues responsible for shell deposition were characterized in three commercially relevant bivalve species. Using high-throughput sequencing and bioinformatics tools, de novo transcriptome assemblies of mantle tissues were generated for the mussel Mytilus edulis, the oyster Crassostrea gigas and the scallop Pecten maximus. These transcriptomes were annotated, and contigs with similarity to proteins known to have shell formation roles in other species were identified. Comparison of the shell formation specific proteins in the three bivalves indicates the possibility of species specific shell proteins.

Characterisation of the mantle transcriptome and biomineralisation genes in the blunt-gaper clam, Mya truncata

Members of the Myidae family are ecologically and economically important, but there is currently very little molecular data on these species. The present study sequenced and assembled the mantle transcriptome of Mya truncata from the North West coast of Scotland and identified candidate biomineralisation genes. RNA-Seq reads were assembled to create 20,106 contigs in a de novo transciptome, 18.81% of which were assigned putative functions using BLAST sequence similarity searching (cuttoff E-value 1E-10). The most highly expressed genes were compared to the Antarctic clam (Laternula elliptica) and showed that many of the dominant biological functions (muscle contraction, energy production, biomineralisation) in the mantle were conserved. There were however, differences in the constitutive expression of heat shock proteins, which were possibly due to the M. truncata sampling location being at a relatively low latitude, and hence relatively warm, in terms of the global distribution of the species. Phylogenetic analyses of the Tyrosinase proteins from M. truncata showed a gene expansion which was absent in L. elliptica. The tissue distribution expression patterns of putative biomineralisation genes were investigated using quantitative PCR, all genes showed a mantle specific expression pattern supporting their hypothesised role in shell secretion. The present study provides some preliminary insights into how clams from different environments - temperate versus polar - build their shells. In addition, the transcriptome data provides a valuable resource for future comparative studies investigating biomineralisation.

Allatostatin-type A, kisspeptin and galanin GPCRs and putative ligands as candidate regulatory factors of mantle function

Allatostatin-type A (AST-A), kisspeptin (KISS) and galanin (GAL) G-protein coupled receptor (GPCR) systems share a common ancestral origin in arthropods and the vertebrates where they regulate metabolism and reproduction. The molluscs are the second most diverse phylum in the animal kingdom, they occupy an important phylogenetic position, and their genome is more similar to deuterostomes than the arthropods and nematodes and thus they are good models for studies of gene family evolution and function. This mini-review intends to extend the current knowledge about AST-A, KISS and GAL GPCR system evolution and their putative function in the mollusc mantle. Comparative evolutionary analysis of the target GPCR systems was established by identifying homologues in genomes and tissue transcriptome datasets available for molluscs and comparing them to those of other metazoan systems. Studies in arthropods have revealed the existence of the AST-A system but the loss of homologues of the KISS and GAL systems. Homologues of the insect AST-AR and vertebrate KISSR genes were found in molluscs but putative GALR genes were absent. Receptor gene number suggested that members of this family have suffered lineage specific evolution during the molluscan radiation. In molluscs, orthologues of the insect AST-A peptides were not identified but buccalin peptides that are structurally related were identified and are putative receptor agonists. The identification of AST-AR and KISSR genes in molluscs strengthens the hypotheses that in metazoans members of the AST-AR subfamily share evolutionary proximity with KISSRs. The variable number of receptors and large repertoire of buccalin peptides may be indicative of the functional diversity of the AST-AR/KISSR systems in molluscs. The identification of AST-A and KISS receptors and ligands in the mantle transcriptome indicates that in molluscs they may have acquired a novel function and may play a role in shell development or sensory detection in the mantle.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

Integration of Transcriptomic and Proteomic Approaches Provides a Core Set of Genes for Understanding of Scallop Attachment

Attachment is an essential physiological process in life histories of many marine organisms. Using a combination of transcriptomic and proteomic approach, scallop byssal proteins (Sbps) and their associated regulatory network genes were investigated for the first time. We built the first scallop foot transcriptome library, and 75 foot-specific genes were identified. Through integration of transcriptomic-proteomic approach, seven unique Sbps were identified. Of them, three showed significant amino acid sequence homology to known proteins. In contrast, the rest did not show significant protein matches, indicating they are possibly novel proteins. Our transcriptomic and proteomic analyses also suggest that post-translational modification may be one of the significant features for Sbps as well. Taken together, our study provides the first multidimensional collection of a core set of genes that may be potentially involved in scallop byssal attachment.

Draft genome assemblies and predicted microRNA complements of the intertidal lophotrochozoans Patella vulgata (Mollusca, Patellogastropoda) and Spirobranchus (Pomatoceros) lamarcki (Annelida, Serpulida)

MicroRNAs (miRNA) are small non-coding RNAs that act post-transcriptionally to regulate gene expression levels. Some studies have indicated that microRNAs may have low homoplasy, and as a consequence the phylogenetic distribution of microRNA families has been used to study animal evolutionary relationships. Limited levels of lineage sampling, however, may distort such analyses. Lophotrochozoa is an under-sampled taxon that includes molluscs, annelids and nemerteans, among other phyla. Here, we present two novel draft genomes, those of the limpet Patella vulgata and polychaete Spirobranchus (Pomatoceros) lamarcki. Surveying these genomes for known microRNAs identifies numerous potential orthologues, including a number that have been considered to be confined to other lineages. RT-PCR demonstrates that some of these (miR-1285, miR-1287, miR-1957, miR-1983 and miR-3533), previously thought to be found only in vertebrates, are expressed. This study provides genomic resources for two lophotrochozoans and reveals patterns of microRNA evolution that could be hidden by more restricted sampling.

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.

Lectin-like molecules in transcriptome of Littorina littorea hemocytes

The common periwinkle Littorina littorea was introduced in the list of models for comparative immunobiology as a representative of phylogenetically important taxon Caenogastropoda. Using Illumina sequencing technology, we de novo assembled the transcriptome of Littorina littorea hemocytes from 182 million mRNA-Seq pair-end 100 bp reads into a total of 15,526 contigs clustered in 4472 unigenes. The transcriptome profile was analyzed for presence of carbohydrate-binding molecules in a variety of architectural contexts. Hemocytes' repertoire of lectin-like proteins bearing conserved carbohydrate-recognition domains (CRDs) is highly diversified, including 11 of 15 lectin families earlier described in animals, as well as the novel members of lectin family found for the first time in mollusc species. The new molluscan lineage-specific domain combinations were confirmed by cloning and sequencing, including the fuco-lectin related molecules (FLReMs) composed of N-terminal region with no sequence homology to any known protein, a middle Fucolectin Tachylectin-4 Pentaxrin (FTP) domain, and a C-terminal epidermal growth factor (EGF) repeat region. The repertoire of lectin-like molecules is discussed in terms of their potential participation in the receptor phase of immune response. In total, immune-associated functions may be attributed to 70 transcripts belonging to 6 lectin families. These lectin-like genes show low overlap between species of invertebrates, suggesting relatively rapid evolution of immune-associated genes in the group. The repertoire provides valuable candidates for further characterization of the gene functions in mollusc immunity.

De novo transcriptome sequencing of the snail Echinolittorina malaccana: identification of genes responsive to thermal stress and development of genetic markers for population studies

Echinolittorina snails inhabit the upper intertidal rocky shore and face strong selection pressures from thermal extremes and fluctuations. Revealing the molecular processes of adaptive significance is greatly obstructed by the scarcity of genomic resource for these taxa. Here, we reported the first comprehensive transcriptome dataset for the genus Echinolittorina. Using Illumina HiSeq 2000 platform, about 52 M and 54 M paired-end clean reads were, respectively, generated for the control and heat-stressed libraries. Totally, 115,211 unique transcript fragments (unigenes) were assembled, with an average length of 453 bp and a N50 size of 492 bp. Approximately one third of the unigenes could be annotated according to their homology matches against the Nr, Swiss-Prot, COG, or KEGG databases, and they were found to represent 23,098 non-redundant genes. Gene expression comparison revealed that 1,267 and 6,663 annotated genes were, respectively, up- and downregulated with at least twofold changes upon heat stress. Gene Ontology and KEGG pathway analyses indicated that there were overrepresented amount of genes enriched in a broad spectrum of biological processes and pathways, including those associated with cytoskeleton organization, developmental regulation, signaling transduction, infection, and cardiac function. In addition, a transcriptome-wide search for polymorphic loci yielded a total of 11,228 simple sequence repeats (SSRs) from 9,938 unigenes and 138,631 single nucleotide polymorphism (SNP) and insertion/deletion (INDEL) sites among 22,770 unigenes. The large number of transcript sequences acquired, the biological pathways identified, and the candidate microsatellite and SNP/INDEL loci discovered in the study will serve as valuable resources for further investigations of genetic differentiation and thermal adaptation among populations.

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.

A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9ppt which showed best growth performance. Total sequence data generated was 467.8Mbp, consisting of 4,116,424 reads with an average length of 112bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

Biomineral proteins from Mytilus edulis mantle tissue transcriptome

The common blue mussel, Mytilus edulis, has a bimineralic shell composed of approximately equal proportions of the two major polymorphs of calcium carbonate: calcite and aragonite. The exquisite biological control of polymorph production is the focus of research interest in terms of understanding the details of biomineralisation and the proteins involved in the process of complex shell formation. Recent advances in ease and availability of pyrosequencing and assembly have resulted in a sharp increase in transcriptome data for invertebrate biominerals. We have applied Roche 454 pyrosequencing technology to profile the transcriptome for the mantle tissue of the bivalve M. edulis. A comparison was made between the results of several assembly programs: Roche Newbler assembler versions 2.3, 2.5.2 and 2.6 and MIRA 3.2.1 and 3.4.0. The Newbler and MIRA assemblies were subsequently merged using the CAP3 assembler to give a higher consensus in alignments and a more accurate estimate of the true size of the M. edulis transcriptome. Comparison sequence searches show that the mantle transcripts for M. edulis encode putative proteins exhibiting sequence similarities with previously characterised shell proteins of other species of Mytilus, the Bivalvia Pinctada and haliotid gastropods. Importantly, this enhanced transcriptome has detected several transcripts that encode proteins with sequence similarity with previously described shell biomineral proteins including Shematrins and lysine-rich matrix proteins (KRMPs) not previously found in Mytilus.

Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica)

As a benthic filter-feeder of estuaries, the eastern oyster, Crassostrea virginica, faces tremendous exposure to microbial pathogens. How eastern oysters without adaptive immunity survive in pathogen-rich environments is of fundamental interest, but studies on its immune system are hindered by the lack of genomic resources. We sequenced the transcriptome of an adult oyster with short Illumina reads and assembled 66,229 contigs with a N50 length of 1,503 bp. The assembly covered 89.4 % of published ESTs and 97.9 % of mitochondrial genes demonstrating its quality. A set of 39,978 contigs and unigenes (>300 bp) were identified and annotated by searching public databases. Analysis of the gene set yielded a diverse set of 657 genes related to innate immunity, including many pertaining to pattern recognition, effectors, signal transduction, cytokines, and apoptosis. Gene families encoding C1q domain containing proteins, CTLD, IAPs, Ig_I-set, and TRAFs expanded in C. virginica and Crassostrea gigas. Many key genes of the apoptosis system including IAP, BAX, BAC-2, caspase, FADD, and TNFR were identified, suggesting C. virginica possess advanced apoptosis and apoptosis-regulating systems. Our results show that short Illumina reads can produce transcriptomes of highly polymorphic genomes with coverage and integrity comparable to that from longer 454 reads. The expansion and high diversity in gene families related to innate immunity, point to a complex defense system in the lophotrochozoan C. virginica, probably in adaptation to a pathogen-rich environment.

Characterization of the Zhikong scallop (Chlamys farreri) mantle transcriptome and identification of biomineralization-related genes

Chlamys farreri is a significant species in aquaculture and fishery in East Asia. A deep understanding of its shell formation by studying the transcriptome of the mantle, a key organ in shell formation, could provide important guidance for its culture. Thus, we sequenced and analyzed the mantle transcriptome of C. farreri. The 77,975 unigenes were generated after Illumina sequencing and de novo assembly. The unigenes were annotated using authoritative databases (non-redundant (NR), COG, Gene Ontology (GO), and KEGG) to obtain functional information. BLASTX alignment was performed between unigenes and reported proteins related to biomineralization. The results identified 53 homologous genes representing 17 matrix proteins, most of which are involved in calcite formation, and 171 homologies with 26 proteins related to general processes of biomineralization. The discovery and unusually high expression of MSP-1 suggested its importance in scallops. Homologous unigenes with aragonite-formation-related matrix proteins were much fewer compared with those related to calcite formation. The results implied that, in C. farreri, the number and proportion of matrix proteins related to aragonite formation is much lower than those related to calcite formation, which was consistent with the proportions of aragonite and calcite in C. farreri shells. Thus, the formation of different polymorphs of calcium carbonate (calcite and aragonite) in molluskan shells is regulated by different groups of proteins. Moreover, 17 candidate unigenes, which are probably involved in biomineralization, were predicted by screening for gene products with secreted domains and tandem-arranged repeat units. Our results contribute to the understanding of biomineralization processes and the evolution of shell formation.