Identification of differentially expressed proteins involved in the early larval development of the Pacific oyster Crassostrea gigas

Proteomic response of early juvenile Pacific oysters (Crassostrea gigas) to temperature

Pacific oysters (Crassostrea gigas) are a valuable aquaculture product that provides important ecosystem benefits. Among other threats, climate-driven changes in ocean temperature can impact oyster metabolism, survivorship, and immune function. We investigated how elevated temperature impacts larval oysters during settlement (19-33 days post-fertilization), using shotgun proteomics with data-independent acquisition to identify proteins present in the oysters after 2 weeks of exposure to 23 °C or 29 °C. Oysters maintained at elevated temperatures were larger and had a higher settlement rate, with 86% surviving to the end of the experiment; these oysters also had higher abundance trends of proteins related to metabolism and growth. Oysters held at 23 °C were smaller, had a decreased settlement rate, displayed 100% mortality, and had elevated abundance trends of proteins related to immune response. This novel use of proteomics was able to capture characteristic shifts in protein abundance that hint at important differences in the phenotypic response of Pacific oysters to temperature regimes. Additionally, this work has produced a robust proteomic product that will be the basis for future research on bivalve developmental processes.

Identification and Characterization of miRNAs and Their Predicted mRNAs in the Larval Development of Pearl Oyster Pinctada fucata

As an important economic shellfish, the pearl oyster, Pinctada fucata, and its larvae are an ideal model for studying molecular mechanisms of larval development in invertebrates. Larval development directly affects the quantity and quality of pearl oysters. MicroRNAs (miRNAs) may play important roles in development, but the effects of miRNA expression on P. fucata early development remain unknown. In this study, miRNA and mRNA transcriptomics of seven different P. fucata developmental stages were analyzed using Illumina RNA sequencing. A total of 329 miRNAs, including 87 known miRNAs and 242 novel miRNAs, and 33,550 unigenes, including 26,333 known genes and 7217 predicted new genes, were identified in these stages. A cluster analysis showed that the difference in the numbers of miRNAs was greatest between fertilized eggs and trochophores. In addition, the integrated mRNA transcriptome was used to predict target genes for differentially expressed miRNAs between adjacent developmental stages, and the target genes were subjected to a gene ontology enrichment analysis. Using the gene ontology annotation, 100 different expressed genes and 95 differentially expressed miRNAs were identified as part of larval development regulation. Real-time PCR was used to identify eight mRNAs and three miRNAs related to larval development. The present findings will be helpful for further clarifying the regulatory mechanisms of miRNA in invertebrate larval development.

The proliferating cell nuclear antigen (PCNA) is a potential proliferative marker in oyster Crassostrea gigas

Proliferating cell nuclear antigen (PCNA) is a crucial eukaryotic replication accessory factor in the regulation of DNA synthesis, which is always used as a proliferation marker for haematopoiesis in vertebrates. In the present study, a homologue of PCNA (named as CgPCNA) with a conserved N-terminal PCNA domain and a C-terminal PCNA domain was identified from oyster Crassostrea gigas. The deduced amino acid sequence of CgPCNA shared 85.4% and 86.6% similarities with the PCNAs identified in Mus musculus and Homo sapiens, respectively. CgPCNA was firstly clustered with PCNAs from molluscs, and then with PCNAs from arthropods to form a group falling into the invertebrate clade in the phylogenic tree. The mRNA transcripts of CgPCNA were detected in all tested tissues with higher expression level in gonad, gills and haemolymph. They were also detected in granulocytes, semi-granulocytes and agranulocytes with no significant differences, but the protein level of CgPCNA in agranulocytes was significantly higher (3.67-fold, p < 0.05) than that in granulocytes. In the haemocytes, CgPCNA was mainly distributed in the nucleus and less in the cytoplasm of haemocytes. CgPCNA protein was observed at the tubule lumen regions of gills vessels, and especially colocalized with the EdU signals. After lipopolysaccharide (LPS) and Vibrio splendidus stimulation, the expression level of CgPCNA mRNA in haemocytes was significantly (p < 0.05) up-regulated at 6 h and 12 h, which was 13.87-fold and 3.89-fold of that in control, respectively. In the oysters treated with the recombinant protein CgAstakine (rCgAstakine), the protein abundance of CgPCNA was enhanced in agranulocytes and gills, while no significant change was observed in semi-granulocytes and granulocytes. These results collectively indicated that CgPCNA was highly expressed in the newborn agranulocytes and the potential haematopoietic sites, and it might be applied as a marker for haemocytes proliferation in oysters.

Single-Larva RNA Sequencing Identifies Markers of Copper Toxicity and Exposure in Early Mytilus californianus Larvae

One of the challenges facing efforts to generate molecular biomarkers for toxins is distinguishing between markers that are indicative of exposure and markers that provide evidence of the effects of toxicity. Phenotypic anchoring provides an approach to help segregate markers into these categories based on some phenotypic index of toxicity. Here we leveraged the mussel embryo-larval toxicity assay in which toxicity is estimated by the fraction of larvae that exhibit an abnormal morphology, to isolate subsets of larvae that were abnormal and thus showed evidence of copper-toxicity, versus others that while exposed to copper exhibited normal morphology. Mussel larvae reared under control conditions or in the presence of increasing levels of copper (3-15 μg/L Cu²⁺) were physically sorted according to whether their morphology was normal or abnormal, and then profiled using RNAseq. Supervised differential expression analysis identified sets of genes whose differential expression was specific to the pools of abnormal larvae versus normal larvae, providing putative markers of copper toxicity versus exposure. Markers of copper exposure and copper-induced abnormality were involved in many of the same pathways, including development, shell formation, cell adhesion, and oxidative stress, yet unique markers were detected in each gene set. Markers of effect appeared to be more resolving between phenotypes at the lower copper concentration, while markers of exposure were informative at both copper concentrations.

Temporal proteomic profiling reveals insight into critical developmental processes and temperature-influenced physiological response differences in a bivalve mollusc

Background

Protein expression patterns underlie physiological processes and phenotypic differences including those occurring during early development. The Pacific oyster (Crassostrea gigas) undergoes a major phenotypic change in early development from free-swimming larval form to sessile benthic dweller while proliferating in environments with broad temperature ranges. Despite the economic and ecological importance of the species, physiological processes occurring throughout metamorphosis and the impact of temperature on these processes have not yet been mapped out.

Results

Towards this, we comprehensively characterized protein abundance patterns for 7978 proteins throughout metamorphosis in the Pacific oyster at different temperature regimes. We used a multi-statistical approach including principal component analysis, ANOVA-simultaneous component analysis, and hierarchical clustering coupled with functional enrichment analysis to characterize these data. We identified distinct sets of proteins with time-dependent abundances generally not affected by temperature. Over 12 days, adhesion and calcification related proteins acutely decreased, organogenesis and extracellular matrix related proteins gradually decreased, proteins related to signaling showed sinusoidal abundance patterns, and proteins related to metabolic and growth processes gradually increased. Contrastingly, different sets of proteins showed temperature-dependent abundance patterns with proteins related to immune response showing lower abundance and catabolic pro-growth processes showing higher abundance in animals reared at 29 °C relative to 23 °C.

Conclusion

Although time was a stronger driver than temperature of metamorphic proteome changes, temperature-induced proteome differences led to pro-growth physiology corresponding to larger oyster size at 29 °C, and to altered specific metamorphic processes and possible pathogen presence at 23 °C. These findings offer high resolution insight into why oysters may experience high mortality rates during this life transition in both field and culture settings. The proteome resource generated by this study provides data-driven guidance for future work on developmental changes in molluscs. Furthermore, the analytical approach taken here provides a foundation for effective shotgun proteomic analyses across a variety of taxa.

Proteome of larval metamorphosis induced by epinephrine in the Fujian oyster Crassostrea angulata

Background

The Fujian oyster Crassostrea angulata is an economically important species that has typical settlement and metamorphosis stages. The development of the oyster involves complex morphological and physiological changes, the molecular mechanisms of which are as yet unclear.

Results

In this study, changes in proteins were investigated during larval settlement and metamorphosis of Crassostrea angulata using epinephrine induction. Protein abundance and identity were characterized using label-free quantitative proteomics, tandem mass spectrometry (MS/ MS), and Mascot methods. The results showed that more than 50% (764 out of 1471) of the quantified proteins were characterized as differentially expressed. Notably, more than two-thirds of the differentially expressed proteins were down-regulated in epinephrine-induced larvae. The results showed that “metabolic process” was closely related to the development of settlement and metamorphosis; 5 × 10^− 4 M epinephrine induced direct metamorphosis of larvae and was non-toxic. Calmodulin and MAPK pathways were involved in the regulation of settlement of the oyster. Expression levels of immune-related proteins increased during metamorphosis. Hepatic lectin-like proteins, cadherins, calmodulin, calreticulin, and cytoskeletal proteins were involved in metamorphosis. The nervous system may be remodeled in larval metamorphosis induced by epinephrine. Expression levels of proteins that were enriched in the epinephrine signaling pathway may reflect the developmental stage of the larvae, that may reflect whether or not larvae were directly involved in metamorphosis when the larvae were treated with epinephrine.

Conclusion

The study provides insight into proteins that function in energy metabolism, immune responses, settlement and metamorphosis, and shell formation in C. angulata. The results contribute valuable information for further research on larval settlement and metamorphosis.

Graphical abstract

Transcriptome analysis and histopathological observations of Geloina erosa gills upon Cr(VI) exposure

The heavy metal contamination like Cr(VI) has been increased by human activities and that threats the ecosystem health of mangrove areas. Bioindicator is an emerging tool in the environmental contamination assessment. The objective of this study was to investigate the Geloina erosa response mechanisms and sensitivities of several biomarkers in the Cr(VI) exposure and identify the G. erosa capability of being used as heavy metals bioindicator. In this study, G. erosa was exposed to 100 μmol·L^-1 Cr(VI) for 48 h. After transcriptome sequencing, a total of 134,817 unigenes were obtained, including 12,555 up-regulated and 18,829 down-regulated differentially expressed genes and were validated through quantitative real-time PCR. In addition, a total of 12,185 SSRs and 1,428,214 candidate SNPs were identified from all the G. erosa transcriptome libraries. Histopathology of the gill indicated the Cr(VI) exposure induced damage of the organ leading to its immunization, detoxification or apoptosis reactions. Among eight genes of the selected biomarkers, Calm, HSP70, CYP450, ATG5, TLR2, MYD88 and CASP8 were up-regulated, while TLR4 was down-regulated in response to the Cr(VI) exposure.

Identification of three cell populations from the shell gland of a bivalve mollusc

The molluscan larval shell formation is a complicated process. There is evidence that the mantle of the primary larva (trochophore) contains functionally different cell populations with distinct gene expression profiles. However, it remains unclear how these cells are specified. In the present study, we identified three cell populations from the shell gland in earlier stages (gastrula) from the bivalve mollusc Crassostrea gigas. These cell populations were determined by analyzing the co-expression relationships among six potential shell formation (pSF) genes using two-color hybridization. The three cell populations, which we designated as SGCPs (shell gland cell populations), formed a concentric-circle pattern from outside to inside of the shell gland. SGCP I was located in the outer edge of the shell gland and the cells expressed pax2/5/8, gata2/3, and bmp2/4. SGCP II was located more internally and the cells expressed two engrailed genes. The last population, SGCP III, was located in the central region of the shell gland and the cells expressed lox4. Determination of the gene expression profiles of SGCPs would help trace their origins and fates and elucidate how these cell populations are specified. Moreover, potential roles of the SGCPs, e.g., development of sensory cells and shell biogenesis, are suggested. Our results reveal the internal organization of the embryonic shell gland at the molecular level and add to the knowledge of larval shell formation.

Dynamic response in the larval geoduck (Panopea generosa) proteome to elevated pCO2

Pacific geoducks (Panopea generosa) are clams found along the northeast Pacific coast where they are important components of coastal and estuarine ecosystems and a major aquaculture product. The Pacific coastline, however, is also experiencing rapidly changing ocean habitat, including significant reductions in pH. To better understand the physiological impact of ocean acidification on geoduck clams, we characterized for the first time the proteomic profile of this bivalve during larval development and compared it to that of larvae exposed to low pH conditions. Geoduck larvae were reared at pH 7.5 (ambient) or pH 7.1 in a commercial shellfish hatchery from day 6 to day 19 postfertilization and sampled at six time points for an in-depth proteomics analysis using high-resolution data-dependent analysis. Larvae reared at low pH were smaller than those reared at ambient pH, especially in the prodissoconch II phase of development, and displayed a delay in their competency for settlement. Proteomic profiles revealed that metabolic, cell cycle, and protein turnover pathways differed between the two pH and suggested that differing phenotypic outcomes between pH 7.5 and 7.1 are likely due to environmental disruptions to the timing of physiological events. In summary, ocean acidification results in elevated energetic demand on geoduck larvae, resulting in delayed development and disruptions to normal molecular developmental pathways, such as carbohydrate metabolism, cell growth, and protein synthesis.

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.

Gene expression correlated with delay in shell formation in larval Pacific oysters (Crassostrea gigas) exposed to experimental ocean acidification provides insights into shell formation mechanisms

Background

Despite recent work to characterize gene expression changes associated with larval development in oysters, the mechanism by which the larval shell is first formed is still largely unknown. In Crassostrea gigas, this shell forms within the first 24 h post fertilization, and it has been demonstrated that changes in water chemistry can cause delays in shell formation, shell deformations and higher mortality rates. In this study, we use the delay in shell formation associated with exposure to CO₂-acidified seawater to identify genes correlated with initial shell deposition.

Results

By fitting linear models to gene expression data in ambient and low aragonite saturation treatments, we are able to isolate 37 annotated genes correlated with initial larval shell formation, which can be categorized into 1) ion transporters, 2) shell matrix proteins and 3) protease inhibitors. Clustering of the gene expression data into co-expression networks further supports the result of the linear models, and also implies an important role of dynein motor proteins as transporters of cellular components during the initial shell formation process.

Conclusions

Using an RNA-Seq approach with high temporal resolution allows us to identify a conceptual model for how oyster larval calcification is initiated. This work provides a foundation for further studies on how genetic variation in these identified genes could affect fitness of oyster populations subjected to future environmental changes, such as ocean acidification.

Electronic supplementary material

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

Proteomic analysis of trochophore and veliger larvae development in the small abalone Haliotis diversicolor

Background

Haliotis diversicolor is commercially important species. The trochophore and veliger are distinct larval stages in gastropod development. Their development involves complex morphological and physiological changes. We studied protein changes during the embryonic development of H. diversicolor using two dimensional electrophoresis (2-DE) and label-free methods, tandem mass spectrometry (MS/ MS), and Mascot for protein identification.

Results

A total of 150 2-DE gel spots were identified. Protein spots showed upregulation of 15 proteins and downregulation of 28 proteins as H. diversicolor developed from trochophore to veliger larvae. Trochophore and veliger larvae were compared using a label-free quantitative proteomic approach. A total of 526 proteins were identified from both samples, and 104 proteins were differentially expressed (> 1.5 fold). Compared with trochophore larvae, veliger larvae had 55 proteins upregulated and 49 proteins downregulated. These differentially expressed proteins were involved in shell formation, energy metabolism, cellular and stress response processes, protein synthesis and folding, cell cycle, and cell fate determination. Compared with the 5 protein (fructose-bisphosphate aldolase, 14–3-3ε, profilin, actin-depolymerizing factor (ADF)/cofilin) and calreticulin) expression patterns, the mRNA expression exhibited similar patterns except gene of fructose-bisphosphate aldolase.

Conclusion

Our results provide insight into novel aspects of protein function in shell formation, torsion, and nervous system development, and muscle system differentiation in H. diversicolor larvae. “Quality control” proteins were identified to be involved in abalone larval development.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4203-7) contains supplementary material, which is available to authorized users.

A comparative proteomic analysis reveals important proteins for the fertilization and early embryonic development of the oyster Crassostrea gigas

Molluscan development involves important features that are important to understanding not only molluscan ontogeny but also animal evolution. To gain insight into the gamete proteome and protein function in fertilization and early development, we analyzed the proteomes of unfertilized oocytes and early embryos (2/4-cell stage) of the Pacific oyster, Crassostrea gigas. An oocyte reference map containing 116 protein spots, of which 69 were identified, revealed a high abundance of vitellogenin-derived protein spots. The differentially regulated protein spots during fertilization were screened using comparative proteomic approaches. In total, 18 differentially regulated protein spots were screened, and 15 of these were identified and divided into three groups. The proteins belonging to the first group function in energy supply and antioxidation and are proposed to ensure successful fertilization by regulating the levels of adenosine triphosphate, resisting oxidative stress, and preventing polyspermy. The proteins of the second group are associated with protein synthesis and modification, reflecting active protein synthesis after fertilization. The three proteins belonging to the final group are hypothesized to function in the regulation of embryonic development through the establishment of cell polarity and modulation of methylation reactions in nuclei. These results will enhance our knowledge of molluscan fertilization and development.

Identification and characterization of miRNAs transcriptome in the South African abalone, Haliotis midae

Aquatic animal diseases are one of the most important limitations to the growth of aquaculture. miRNAs represent an important class of small ncRNAs able to modulate host immune and stress responses. In Mollusca, a large phylum of invertebrates, miRNAs have been identified in several species. The current preliminary study identified known miRNAs from the South African abalone, Haliotis midae. The economic and ecological importance of abalone makes this species a suitable model for studying and understanding stress response in marine gastropods. Furthermore, the identification of miRNA, represents an alternative and powerful tool to combat infectious disease.

Assessment of housekeeping genes as internal references in quantitative expression analysis during early development of oyster

The early development of mollusks exhibits important characteristics from the developmental and evolutionary perspective. With the increasing number of genome-wide studies, accurate analyses of quantitative gene expression during development are impeded by the lack of validated reference genes. To improve the situation, in this study, we analyzed the expression stability of seven candidate housekeeping genes during early development of the Pacific oyster Crassostrea gigas: actin, glyceraldehyde-3-phosphate dehydrogenase (gapdh), α subunit of elongation factor 1 (elf1α), adp-ribosylation factor 1 (arf1), heterogeneous nuclear ribonucleoprotein q, ubiquitin-conjugating enzyme e2d2 and ribosomal protein s18. We focused on 11 stages from oocyte to D-veliger, which include crucial developmental processes such as axis determination, gastrulation and shell formation. Gene expression stabilities were assessed with the three commonly used programs geNorm, NormFinder and BestKeeper. Although the results obtained with the three programs varied to some extent, in general, arf1, elf1α and gapdh were highly ranked and actin was poorly ranked. This analysis also indicated that multiple genes should be used for normalization, and we concluded that arf1-elf1α-gapdh should be used as internal references. The findings of this study will help researchers to obtain accurate results in future quantitative gene expression analysis of development in bivalve mollusks.

Characterization and expression of a novel caspase gene: Evidence of the expansion of caspases in Crassostrea gigas

Caspases are a group of cysteine-aspartate proteases involved in apoptosis and a variety of non-apoptotic processes. In this study, a novel caspase gene was cloned and its potential role in apoptosis was investigated. The caspase gene (CgCasp 3/7) has an open reading frame of 1626bp encoding 541 amino acids containing the conserved functional domains and motifs of effector caspases. Its amino acid sequence shows low identity with the other effector caspases of Crassostrea gigas and contains a unique long intersubunit linker (IL). The CgCasp 3/7 mRNA was expressed highly in oocytes and then decreased gradually after fertilization, indicating CgCasp 3/7 could function in oocyte apoptosis. In adult tissues, it is located primarily in the gills and hepatopancreas. We examined the mRNA expression of CgCasp 3/7 in gills of oysters immersed in ambient (17°C) or heated (27°C) seawater. The thermal stress stimulated mRNA expression of CgCasp 3/7 by 2.5- and 4.1-fold at 2h and 6h post-treatment, respectively, indicating CgCasp3/7 was involved in the early response to thermal stress. To examine the function of the IL, CgCasp 3/7 and CgCasp 3/7-T (with a truncated IL) were expressed using an in vitro translation system and their DEVDase activity was measured. Both proteins showed a significantly higher level of DEVDase activity than control, but CgCasp3/7-T had lower DEVDase activity than CgCasp3/7, indicating CgCasp3/7 had DEVDase activity and the IL was required for maximal DEVDase activity. Our study adds to the complexity of caspases in C. gigas.

Developmental Transcriptome Analysis and Identification of Genes Involved in Larval Metamorphosis of the Razor Clam, Sinonovacula constricta

The razor clam Sinonovacula constricta is an important commercial species. The deficiency of developmental transcriptomic data is becoming the bottleneck of further researches on the mechanisms underlying settlement and metamorphosis in early development. In this study, de novo transcriptome sequencing was performed for S. constricta at different early developmental stages by using Illumina HiSeq 2000 paired-end (PE) sequencing technology. A total of 112,209,077 PE clean reads were generated. De novo assembly generated 249,795 contigs with an average length of 585 bp. Gene annotation resulted in the identification of 22,870 unigene hits against the NCBI database. Eight unique sequences related to metamorphosis were identified and analyzed using real-time PCR. The razor clam reference transcriptome would provide useful information on early developmental and metamorphosis mechanisms and could be used in the genetic breeding of shellfish.

Genome-wide identification and characterization of long intergenic noncoding RNAs and their potential association with larval development in the Pacific oyster

An increasing amount of evidence suggests that long intergenic noncoding RNAs (lincRNAs) may play diverse roles in many cellular processes. However, little is known about lincRNAs in marine invertebrates. Here, we presented the first identification and characterization of lincRNAs in the Pacific oyster (Crassostrea gigas). We developed a pipeline and identified 11,668 lincRNAs in C. gigas based on RNA-Seq resources available. These lincRNAs exhibited many common characteristics with vertebrate lincRNAs: relatively short length, low exon numbers, low expression, and low sequence conservation. 1,175 lincRNAs were expressed in a tissue-specific manner, with 35.2% preferentially expressed in male gonad. 776 lincRNAs were specifically expressed in juvenile during different developmental stages. In addition, 47 lincRNAs were found to be potentially related to oyster settlement and metamorphosis. Such diverse temporal and spatial patterns of expression suggest that these lincRNAs might function in cell differentiation during early development, as well as sex differentiation and reproduction. Based on a co-expression network analysis, five lincRNAs were detected that have an expression correlation with key hub genes in four modules significantly correlated with larval development. Our study provides the first large-scale identification of lincRNAs in molluscs and offers new insights into potential functions of lincRNAs in marine invertebrates.

Non-additive effects of ocean acidification in combination with warming on the larval proteome of the Pacific oyster, Crassostrea gigas

Increasing atmospheric carbon dioxide results in ocean acidification and warming, significantly impacting marine invertebrate larvae development. We investigated how ocean acidification in combination with warming affected D-veliger larvae of the Pacific oyster Crassostrea gigas. Larvae were reared for 40h under either control (pH8.1, 20 °C), acidified (pH7.9, 20 °C), warm (pH8.1, 22 °C) or warm acidified (pH7.9, 22 °C) conditions. Larvae in acidified conditions were significantly smaller than in the control, but warm acidified conditions mitigated negative effects on size, and increased calcification. A proteomic approach employing two-dimensional electrophoresis (2-DE) was used to quantify proteins and relate their abundance to phenotypic traits. In total 12 differentially abundant spots were identified by nano-liquid chromatography-tandem mass spectrometry. These proteins had roles in metabolism, intra- and extra-cellular matrix formations, stress response, and as molecular chaperones. Seven spots responded to reduced pH, four to increased temperature, and six to acidification and warming. Reduced abundance of proteins such as ATP synthase and GAPDH, and increased abundance of superoxide dismutase, occurred when both pH and temperature changes were imposed, suggesting altered metabolism and enhanced oxidative stress. These results identify key proteins that may be involved in the acclimation of C. gigas larvae to ocean acidification and warming.

SIGNIFICANCE

Increasing atmospheric CO2 raises sea surface temperatures and results in ocean acidification, two climatic variables known to impact marine organisms. Larvae of calcifying species may be particularly at risk to such changing environmental conditions. The Pacific oyster Crassostrea gigas is ecologically and commercially important, and understanding its ability to acclimate to climate change will help to predict how aquaculture of this species is likely to be impacted. Modest, yet realistic changes in pH and/or temperature may be more informative of how populations will respond to contemporary climate change. We showed that concurrent acidification and warming mitigates the negative effects of pH alone on size of larvae, but proteomic analysis reveals altered patterns of metabolism and an increase in oxidative stress suggesting non-additive effects of the interaction between pH and temperature on protein abundance. Thus, even small changes in climate may influence development, with potential consequences later in life.

A Label-Free Proteomic Analysis on Competent Larvae and Juveniles of the Pacific Oyster Crassostrea gigas

Current understandings on the molecular mechanisms underlying bivalve metamorphosis are still fragmentary, and a comprehensive description is required. In this study, using a large-scale label-free proteomic approach, we described and compared the proteomes of competent larvae (CL) and juveniles (JU) of the Pacific oyster, Crassostrea gigas. A total of 788 proteins were identified: 392 in the CL proteome and 636 in the JU proteome. Gene Ontology analysis of the proteome from each sample revealed active metabolic processes in both stages. Further quantitative analyses revealed 117 proteins that were differentially expressed between the two samples. These proteins were divided into eight groups: cytoskeleton and cell adhesion, protein synthesis and degradation, immunity and stress response, development of particular tissues, signal regulation, metabolism and energy supply, transport, and other proteins. A certification experiment using real-time PCR assay confirmed 20 of 30 examined genes exhibited the same trends at the mRNA and protein levels. The differentially expressed proteins may play roles in tissue remodeling, signal transduction, and organ development during and after metamorphosis. Novel roles were proposed for some differentially expressed proteins, such as chymotrypsin. The results of this work provide an overview of metamorphosis and post-metamorphosis development of C. gigas at the protein level. Future studies on the functions of the differentially expressed proteins will help to obtain a more in-depth understanding of bivalve metamorphosis.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Identification of conserved and novel microRNAs in the Pacific oyster Crassostrea gigas by deep sequencing

MicroRNAs (miRNAs) play important roles in regulatory processes in various organisms. To date many studies have been performed in the investigation of miRNAs of numerous bilaterians, but limited numbers of miRNAs have been identified in the few species belonging to the clade Lophotrochozoa. In the current study, deep sequencing was conducted to identify the miRNAs of Crassostrea gigas (Lophotrochozoa) at a genomic scale, using 21 libraries that included different developmental stages and adult organs. A total of 100 hairpin precursor loci were predicted to encode miRNAs. Of these, 19 precursors (pre-miRNA) were novel in the oyster. As many as 53 (53%) miRNAs were distributed in clusters and 49 (49%) precursors were intragenic, which suggests two important biogenetic sources of miRNAs. Different developmental stages were characterized with specific miRNA expression patterns that highlighted regulatory variation along a temporal axis. Conserved miRNAs were expressed universally throughout different stages and organs, whereas novel miRNAs tended to be more specific and may be related to the determination of the novel body plan. Furthermore, we developed an index named the miRNA profile age index (miRPAI) to integrate the evolutionary age and expression levels of miRNAs during a particular developmental stage. We found that the swimming stages were characterized by the youngest miRPAIs. Indeed, the large-scale expression of novel miRNAs indicated the importance of these stages during development, particularly from organogenetic and evolutionary perspectives. Some potentially important miRNAs were identified for further study through significant changes between expression patterns in different developmental events, such as metamorphosis. This study broadened the knowledge of miRNAs in animals and indicated the presence of sophisticated miRNA regulatory networks related to the biological processes in lophotrochozoans.

The identification and characteristics of immune-related microRNAs in haemocytes of oyster Crassostrea gigas

Background

MicroRNAs (miRNAs) represent a class of small ncRNAs that repress gene expression on the post-transcriptional level by the degradation or translation inhibition of target mRNA.

Methodology

Three small RNA libraries from oyster haemocytes were sequenced on the Illumina platform to investigate the latent immunomodulation of miRNAs after bacteria challenge and heat stress. Totally, 10,498,663, 8,588,606 and 9,679,663 high-quality reads were obtained in the control, bacteria and bacteria+heat library, respectively, from which 199 oyster miRNAs including 71 known and 128 novel ones were identified. Among these miRNAs, 6 known and 23 novel ones were predicted to possess more than one precursor-coding region, and cgi-miR-10a, cgi-miR-184b, cgi-miR-100, cgi-miR-1984 and cgi-miR-67a were observed to be the most abundant miRNAs in the control library. The expression levels of 22 miRNAs in the bacteria library were significantly higher than those in the control library, while there were another 33 miRNAs whose expression levels were significantly lower than that in the control library. Meanwhile, the expression levels of 65 miRNAs in the bacteria+heat library changed significantly compared to those in the bacteria library. The target genes of these differentially expressed miRNAs were annotated, and they fell in immune and stress-related GO terms including antioxidant, cell killing, death, immune system process, and response to stimulus. Furthermore, there were 42 differentially expressed miRNAs detected in both control/bacteria and bacteria/bacteria+heat comparisons, among which 9 miRNAs displayed the identical pattern in the two comparisons, and the expression alterations of 8 miRNAs were confirmed using quantitative RT-PCR.

Conclusions

These results indicated collectively that immune challenge could induce the expression of immune-related miRNAs, which might modulate the immune response such as redox reaction, phagocytosis and apoptosis, and the expression of some immune-related miRNAs could be also regulated by heat stress to improve the environmental adaption of oyster.

Identification of a tyrosinase gene potentially involved in early larval shell biogenesis of the Pacific oyster Crassostrea gigas

The larval shell emerges early in embryogenesis of mollusks, but the detailed mechanisms of its biogenesis remain to be determined. In this study, we cloned a tyrosinase gene (cgi-tyr1) that potentially functioned in larval shell biogenesis from the Pacific oyster Crassostrea gigas, a worldwide bivalve species. Sequence analysis of cgi-tyr1 revealed that it had typical copper-binding domains and a signal peptide. Through whole mount in situ hybridization and an electron scanning microscopic observation, we detected the expression of cgi-tyr1 firstly in the saddle-shaped shell field in trochophores, indicating that cgi-tyr1 might participate in the biogenesis of the initial non-calcified shell of trochophores. In the following development to early D-veliger, cells in the central region of shell field exhibited no detectable cgi-tyr1 expression, and cgi-tyr1 expression was sustained only in the edge of the shell field and the hinge region, indicating that cgi-tyr1 might function fundamentally in shell growth from trochophore to early D-veliger. Unexpectedly, cgi-tyr1 expression was not detected after the D-veliger stage. This indicated that other molecules might function in later shell development. Our results suggested a role for a tyrosinase gene that specifically functioned in the initial phase of the larval shell biogenesis of C. gigas. This work would shed light on future studies on larval shell development and might be helpful to understand how the molluscan shell emerged during evolution.

Identification of differential expressed proteins and characterization their mRNA expression in thermally stressed Apostichopus japonicus

In this study, we present a comparative proteomic analysis of the global protein expression changes in sea cucumber after 7 days exposure at 25°C. Using two-dimensional electrophoresis followed by MALDI-TOF MS/MS, 27 protein spots with significant differences in abundance were identified and characterized. The identified proteins belonged primarily to the following four functional categories: cytoskeletal, material and energy metabolism, calcium homeostasis and extracellular matrix. The mRNA expression levels of 7 differentially expressed proteins were further assessed by qRT-PCR. The expression levels of 6 genes, including collagen, ATP synthase, major yolk protein, ferritin, nectin and protein disulfide isomerase showed significant differences under thermal stress, and among them, only two genes-ATP synthase and major yolk protein-showed consistent levels of protein and mRNA expression. Our results offer insight into the complex changes in protein turnover during higher temperature exposure in sea cucumber.

Shotgun proteomics as a viable approach for biological discovery in the Pacific oyster

Shotgun proteomics offers an efficient means to characterize proteins in a complex mixture, particularly when sufficient genomic resources are available. In order to assess the practical application of shotgun proteomics in the Pacific oyster (Crassostrea gigas), liquid chromatography coupled with tandem mass spectrometry was used to characterize the gill proteome. Using information from the recently published Pacific oyster genome, 1043 proteins were identified. Biological samples (n = 4) and corresponding technical replicates (three) were similar in both specific proteins identified and expression, as determined by normalized spectral abundance factor. A majority of the proteins identified (703) were present in all biological samples. Functional analysis of the protein repertoire illustrates that these proteins represent a wide range of biological processes, supporting the dynamic function of the gill. These insights are important for understanding environmental influences on the oyster, because the gill tissue acts as the interface between the oyster and its environment. In silico analysis indicated that this sequencing effort identified a large proportion of the complete gill proteome. Together, these data demonstrate that shotgun sequencing is a viable approach for biological discovery and will play an important role in future studies of oyster physiology.