Figures
Abstract
Many microarray and suppression subtractive hybridization (SSH) studies have analyzed the effects of environmental stress on gene transcription in marine species. However, there have been no unifying analyses of these data to identify common stress response pathways. To address this shortfall, we conducted a meta-analysis of 14 studies that investigated the effects of different environmental stressors on gene expression in oysters. The stressors tested included chemical contamination, hypoxia and infection, as well as extremes of temperature, pH and turbidity. We found that the expression of over 400 genes in a range of oyster species changed significantly after exposure to environmental stress. A repeating pattern was evident in these transcriptional responses, regardless of the type of stress applied. Many of the genes that responded to environmental stress encoded proteins involved in translation and protein processing (including molecular chaperones), the mitochondrial electron transport chain, anti-oxidant activity and the cytoskeleton. In light of these findings, we put forward a consensus model of sub-cellular stress responses in oysters.
Citation: Anderson K, Taylor DA, Thompson EL, Melwani AR, Nair SV, Raftos DA (2015) Meta-Analysis of Studies Using Suppression Subtractive Hybridization and Microarrays to Investigate the Effects of Environmental Stress on Gene Transcription in Oysters. PLoS ONE 10(3): e0118839. https://doi.org/10.1371/journal.pone.0118839
Academic Editor: John A. Craft, Glasgow Caledonian University, UNITED KINGDOM
Received: May 19, 2014; Accepted: January 15, 2015; Published: March 13, 2015
Copyright: © 2015 Anderson et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: Original data can be found in references 6-9, 13-21 and 42.
Funding: This was work supported by Macquarie University (http://www.mq.edu.au/) and by an Australian Research Council Discovery grant (DP120101946, http://www.arc.gov.au/) to DAR and others. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human activity has introduced a range of stressors into coastal marine environments around the world. These stressors exacerbate natural variability and ramify throughout all levels of biological organization, from gene expression to the broader biosphere [1]. The main sources of environmental stress in coastal areas include chemical contamination, lack of food, altered hydrological parameters (temperature, pH, dissolved oxygen, turbidity etc.) and disease. Sewage, industrial activity, agricultural runoffs and storm water outfalls represent major sources of coastal pollution. They introduce organic and inorganic chemicals, such as eutrophying nutrients, toxic metals, polycyclic aromatic hydrocarbons (PAH), pesticides and polychlorinated biphenyls, into coastal waterways [2–5]. Hydrological changes resulting from coastal engineering can exacerbate natural fluctuations in water quality, particularly by increasing the frequency of hypoxic conditions during summer months and by changing localized water temperature, pH and dissolved CO2 concentrations [5].
Numerous studies have investigated the impacts of environmental stress on gene expression in coastal marine species (e.g., [6–9]). These studies have defined a range of genes that are affected by specific types of stress. However, the data have not been integrated to test whether there is a universal environmental stress response that is common to a broad range of different stressors. Understanding whether there is such a generic reaction to stress at the sub-cellular level will be critical for the development of management practices that can withstand multifaceted anthropogenic changes to the environment.
The current review addresses this deficit by undertaking a meta-analysis of transcriptome data to identify the key molecular systems that are affected by environmental stress in edible oysters. Oysters (phylum Mollusca, class Bivalvia, order Ostreoida, family Ostreidae) were chosen for this study because of their significance in global aquaculture, their role in significantly modifying habitats (ecosystem engineers), and their frequent use in environmental biomonitoring [10,11]. Oyster farming often relies on monocultures in impacted coastal environments. Hence, the detrimental effects of environmental stress are likely to be felt severely in these aquaculture industries [12]. The potential economic impacts of environmental change on oyster farming worldwide also means that substantial effort has been put into analyzing the effects of stress on oysters.
Data mining and analysis
Data for inclusion in the meta-analysis were selected by searching GenBank’s nucleotide and expressed sequence tag (EST) databases (http://www.ncbi.nlm.nih.gov) (Fig. 1). These databases were searched in 2013 for differentially expressed genes identified by SSH or cDNA microarray analyses of oysters from the genera Crassostrea, Ostrea and Saccostrea that had been exposed to environmental stressors (as defined by their GenBank annotations). Only microarray analyses that had been effectively validated by quantitative PCR were included in the analysis. Pertinent next generation transcriptomic and shotgun proteomic studies are discussed separately.
Once the full set of differentially expressed sequences with their associated accession numbers had been identified, they were assembled into a single database and re-analyzed to confirm their identity. In all of the following analyses, we used a cut-off ≤ e-20 to represent “robust” homology. All hits that did not meet this threshold were discarded. First, accession numbers associated with the differentially expressed GenBank entries were searched using Blastn for “somewhat similar sequences” against the non-redundant “other” nucleotide database on NCBI. Robust hits (hits against homologous genes with ≤ e-20 matches) were included in our raw database. All genes that failed to meet this criterion using Blastn were searched again using Blastx (“translated nucleotide query against protein database”) with the default settings. Genes identified using Blastx searches with ≤ e-20 were also included in the raw database.
When the top blast hits in Blastn or Blastx searches had gene names that were different to that assigned to the original oyster sequence, both names were searched via Uniprot (http://www.uniprot.org/) to determine whether they were synonyms. For instance, NADH dehydrogenase (EC 1.6.99.3) was deemed to be synonymous with cytochrome c reductase, NADH oxidoreductase and 13 other gene names; see http://enzyme.expasy.org/EC/1.6.99.3). In many cases, the alternative gene name was added to our database annotation to avoid redundancy. In addition, if two or more of the entries in our database of oyster genes encoded the same gene but had been allocated different names in the literature, those redundant sequences were consolidated into a single entry. In some additional cases, genes encoded in the mitochondrial (Mt) genome were allocated accession numbers in the primary literature that corresponded to the entire Mt genome. These were re-allocated gene-specific accessions.
After this filtering process, the final consolidated database was re-analyzed to verify the putative functions of the differentially expressed genes. These functional assignments were initially checked against searches of the literature cited in the GenBank registration file for each sequence. The functional category into which the gene fell was then confirmed by reference to the UniProtKB Protein Knowledgebase (http://www.uniprot.org/uniprot/) and the primary literature. The results from these searches were cross-referenced and hierarchically preferenced so that a single biological function could be assigned to each gene sequence in the database. Preference was given to ontologies that were based on experimentally defined functions in oysters, other bivalves, and mollusks in general.
As a result of these multiple filters, only sequences from published studies with definitive gene identifications and robust functional annotations were included in subsequent analyses.
At the time of this analysis, there were 97,811 nucleotide and 223,784 EST sequences listed in GenBank under the three oyster genera investigated (Table 1). Ninety-eight percent of the CoreNucleotide sequences that we searched against were from the genus Crassostrea, 0.6% were from Saccostrea, and 1.5% were from Ostrea. More than 99% of the dbEST sequences were from the genus Crassostrea. These relative numbers of transcriptome sequences match closely with worldwide aquaculture production levels for the three genera. Over 99% of the edible oysters produced globally in 2009 were from the genus Crassostrea (primarily C. gigas; http://www.fao.org/). This demonstrates that, in the case of oysters, the depth of available sequence resources matches the worldwide economic importance of particular species.
Structure of the differentially expressed gene database
2,111 of the transcript sequences in GenBank databases were annotated as being differentially expressed in oysters that had been exposed to environmental stress. The filtering process described above, which checked the identity and putative functions of the differentially expressed genes, compacted the initial list of 2,111 sequences into a consolidated database of 617 validated non-redundant genes. The final edited list of these genes is shown in S1 Table. When genes with no known function and those not assigned a functional category (as well as equivocal entries) were removed, this list was reduced to 586 known genes with predicted functions. The differentially expressed transcripts were identified in 14 studies that used cDNA microarrays (four studies) or SSH (10 studies) to test the effects of environmental stress on oysters (Table 2). Six of the 14 studies tested gene expression in the hemolymph of oysters, four used digestive gland tissue, and the remainder used a combination of gills, mantle and/or gonads as the source of mRNA. The environmental stressors to which oysters were exposed included, hydrocarbons and other organic contaminants, metals, infectious agents, pesticides, hypoxia, extremes of pH, and thermal stress.
Six of the 14 studies included in the meta-analysis provided consistent data on the up- or down- regulation of differentially expressed genes [6, 13, 15, 18, 20, 21]. Of the 1352 differential genes identified in these studies, 59% ± 18% were upregulated in response to stress. However, a further time course analysis [9], showed that individual genes could either be up- or down- regulated depending on the duration of exposure to stress.
Ontology of differentially expressed genes
The functional ontologies of the mined genes showed that 12 discrete intracellular processes were affected by ‘stress’. The abbreviated titles for these 12 functional categories are shown in Table 3. The percentages cited in the following discussion represent the percentage of the total differential transcriptome (586 genes) comprising a particular functional category. Eight percent of the differentially expressed genes had no known function and thus could not be assigned to a functional category. Only 1% of the differentially expressed genes fell into functional categories other than the 12 listed in Table 3.
Genes associated with metabolism were the most frequently identified (Fig. 2; 28%). These included genes involved in the mitochondiral electron transport chain (10%) and in glycolysis and the tricarboxylic acid (TCA/Krebs) cycle (2%). Genes associated with translation and post translational processing were the next most abundant category, comprising 26% of the differential transcriptome. Most of these genes encoded ribosomal proteins or rRNA (18%). Stress response genes, including molecular chaperones (2%) and enzymes involved in antioxidant responses (including superoxide dismutase, SOD, and peroxiredoxin, 2%), accounted for 11% of differential transcripts. Communication, the cytoskeleton, transcription, the cell cycle and protein regulation each represented approximately 10% of the differential transcriptome. Genes involved in immune responses and membrane transport each accounted for approximately 7% of differentially expressed transcripts. Nucleic acid regulation and the extracellular matrix were the least represented functional catgories (3% or less of the differential transcriptome). The most frequently identified genes in each functional group are shown in Table 4.
The direction of change (up- or down- regulation) of genes comprising the different functional categories is shown in Fig. 3 (using data from [6, 13, 15, 18, 20, 21]). The proportion of genes that were up-regulated, relative to down-regulation, was significantly greater (t-test, p < 0.05) in the cellular stress response, cytoskeleton and electron transport chain categories.
Bars = SEM (n = 6 different studies that verifiably reported direction of change), * denotes a significant difference between up- and down-regulation for that category (t-test, p < 0.05).
The individual genes most frequently affected by stress are involved in the mitochondrial electron transport chain, translation and protein processing, subcellular stress responses and the cytoskeleton
None of the 586 individual genes identified in this analysis was found to be differentially expressed in all 10 of the stress treatments analysed (5 different types of infectious agents, two temperature extremes, two types of contaminants, and hypoxia) (Fig. 4). Actin β was the most frequently affected gene. It was differentially expressed in all treatments except “pesticide/herbicide”. Cytochrome oxidase subunit 1 (CO1) expression was effected by eight out of the 10 stressors. Commonly affected genes such as actin and CO1 were relatively rare. Only 4% of the 586 genes in this analysis were differentially expressed in more than 6 of the treatments, and the expression of 60% of genes was affected by just one of the 10 stress treatments.
Genes are listed in ranked order from most frequently affected to least frequently affected. Shaded boxes represent significantly altered gene expression in response to different stressors.
Fig. 5 shows the 10 most common differentially expressed genes (excluding ribosomes) across all stress treatments identified by name rather than accession number. These “names” represented groups of genes that encoded either different polypeptides of the same protein (for instance NADH dehydrogenase subunits), functionally related homologues in the same gene family (for example HSP70s), or synonyms for the same protein (for instance NADH dehydrogenase and NADH oxidoreductase). This pooling meant that the same “name” could occur more than once in the same stress treatment. Together, the top ten “names” accounted for 31% of the 586 differentially expressed genes identified across the 10 stress treatments. NADH dehydrogenases were by far the most commonly identified differential transcripts. By themselves, they accounted for 7% of the differential transcriptome. Their transcript abundance was significantly affected by all stressors. Among the top 10 gene names, those associated with the mitochondrial electron transport chain (NADH dehydrogenases, cytochrome c’s and ATP synthases) represented 16% of the entire differential transcriptome, whilst cytoskeletal (actin and tubulin) and stress response (metallothioneins, GSTs, HSP70s, and the antioxidant enzyme SOD) genes each comprised 6% of all differential transcripts.
Numbers in brackets represent the number of times this gene name appeared as differentially expressed in all of the different stress treatments. These data do not include genes encoding ribosomal proteins.
Effects of different types of stress on gene expression
Too few differentially expressed genes were identified in some of the 10 different stress treatments to allow robust statistical comparisons between all treatments. So, to allow us to compare the effects of different types of stressors, we collapsed the analysis to compare just 4 different general classes of stressor (infection, temperature, contamination and hypoxia). We also eliminated the least affected functional categories of genes (extracellular matrix, nucleic acid regulation and membrane transport), and combined transcription and translation into a single functional category.
The percentage of differentially expressed genes comprising the different functional categories was broadly similar between all four general classes of stressors (Fig. 6). Genes in the transcription/translation and metabolism categories were the most frequently affected by all classes of stress. However, there was significant variability (χ2 = 73.23, df = 21, p < 0.05) in the responses to different stressors. Temperature extremes affected the expression of significantly more genes involved in transcription and/or translation than the other types of stress. Similarly, hypoxia and contamination had a proportionally greater effect on genes involved in metabolism when compared to temperature extremes and infection. Despite these minor variations, the transcriptional profiles of oysters responding to the different stressors were quite similar. This suggests that there is a core set of stress response genes within oyster cells, primarily involved in metabolism, and that the transcription/translation of these genes is responsive to a broad range of environmental stressors. By corollary, the data imply that enhanced cellular energy demand is a universal response to environmental stress, and that adaptive responses to stress extend well beyond the “minimal stress proteome of cellular organisms” proposed by Kültz [22].
Genomic, next generation transcriptomic and proteomic analyses of environmental stress in oysters
A complete genome sequence is now available for Pacific oysters (C. gigas) [23], and “omics” data for edible oysters continues to expand [24]. The C. gigas genome project revealed substantial expansion, relative to other species, of gene families such as heat shock protein 70s that are involved in environmental stress responses. This was deemed to be an adaptive response to the highly stressful intertidal habitats occupied by C. gigas [23]. Other recent studies have used next generation analytical techniques to investigate stress responses in oysters; including RNAseq transcriptomic responses to salinity stress [25], virulent Vibrio strains [26] and Roseovarius Oyster Disease (ROD) [27], and shotgun proteomic analyses of heavy metal contamination [28] and elevated pCO2 [29]. These studies found differential responses to environmental stress in many of the same subcellular pathways identified in the current meta-analysis. For instance, Timmins-Schiffman et al. [29] showed that changes in the proteomes of C. gigas responding to elevated pCO2 included pathways involved in antioxidant activity, carbohydrate metabolism and transcription/ translation. Similarly, a shotgun proteomic analysis by Muralidharan et al. [28] found that proteins encoded by six of the ten most commonly affected genes identified in the current meta-analysis (Fig. 4) were also differentially regulated in S. glomerata after exposure to heavy metals. The differential proteins included NADH dehydrogenase, cytochrome c oxidase, tubulin, actin and ATP synthase. These data suggest that next generation studies investigating stress responses in oysters tend to identify the same intracellular stress response pathways as those revealed by SSH and microarrays. Based on the relatively few existing next generation analyses, these newer techniques do not seem likely to discover entirely novel mechanisms. However, their greater analytical depth will no doubt detect more components of pathways already identified by previous studies, and they may be able to discern subtle differences in responses to different types of stress. This might be particularly useful if the goal of analyses is to identify suites of biomarkers for environmental monitoring.
Conclusions and perspectives
This meta-analysis has compiled data from 14 different transcriptomic studies of oysters exposed to a range of environmental conditions. It has shown that a number of key sub-cellular systems are consistently involved in transcriptional responses to a broad range of environmental stresses. The data support a consensus model that describes a generic intracellular stress response in oysters (Fig. 7). Elements of this model, as it pertains to oysters and other organisms, have been discussed in detail by other authors [30–31]. The model highlights the importance of the mitochondrial electron transport chain, antioxidant enzymes, molecular chaperones and cytoskeletal proteins. Our meta-analysis has shown that the mitochondrial electron transport chain is the primary cellular system impacted by environmental stress and so it is central to the consensus model. Even though enhanced energy production in the mitochondria may be beneficial in terms of powering adaptive cellular processes, it is also known to elevate the production of cytotoxic reactive oxygen species (ROS), with consequent impacts on mitochondrial and cytoskeletal integrity. The failure of anti-oxidant enzymes such as SOD and molecular chaperones to limit damage caused by ROS is likely to result in cellular dysfunction and ultimately apoptotic cell death. All of these effects have been demonstrated in numerous studies of cellular function in oysters responding to stress [32–49]. Modulation of the mitochondrial electron transport chain and its downstream consequences, as depicted in this consensus model, may explain many of the impacts that environmental stress exerts at levels of biological organization from cells to ecosystems [32].
In this model, adaptive intracellular responses to stress lead to increased energy production (ATP) via the mitochondrial electron transport chain. Reactive oxygen species (ROS) are produced as by-products of ATP synthesis. Unless controlled by molecular chaperones (HSPs) and the anti-oxidant system, excessive ROS can disrupt the actin cytoskeleton and mitochondrial membranes leading to apoptotic cell death.
Supporting Information
S1 Table. The complete list of differentially expressed genes affected by environmental stress identified in this meta-analysis showing; the gene or encoded protein name, the intracellular process to which the gene was assigned, the accession number for the gene and the type of stress that significantly affected the expression of the gene (designated by X in the corresponding box).
Abbreviations: Comms, Communication; QX, QX disease/Martelia sydneyi infection, B. ost, Bonamia ostrea; V. spl, Vibrio splendidus; HV1. Ostreid Herpesvirus 1; P. mar, Perkinsus marinus; temp high, high temperature; temp low, low temperature; HC, hydorcarbons; pest/herb, pesticides/herbicides; hypox, hypoxia.
https://doi.org/10.1371/journal.pone.0118839.s002
(PDF)
Acknowledgments
D. Taylor was the recipient of an iMQRes and an Endeavour International Postgraduate Research Scholarship from Macquarie University. A Melwani is supported by a Thyne-Reid Fellowship from the Sydney Institute of Marine Science. The work was also supported in part under Australian Research Council's Discovery Projects funding scheme (project number DP120101946).
Author Contributions
Conceived and designed the experiments: DAT SVN DAR. Performed the experiments: DAT. Analyzed the data: KA DAT ARM DAR. Contributed reagents/materials/analysis tools: ELT KA. Wrote the paper: ARM KA DAT DAR. Downloaded data from NCBI: KA.
References
- 1. Peñuelas J, Sardans J, Estiarte M, Ogaya R, Carnicer J, Coll M, et al. Evidence of current impact of climate change on life: a walk from genes to the biosphere. Global Change Biol 2013;19: 2303–2338. pmid:23505157
- 2. Medeiros ID, Siebert MN, Toledo-Silva G, Rodrigues TB, Marques MRF, Bainy AC. Induced gene expression in oyster Crassostrea gigas exposed to sewage. Environ Toxicol Pharmacol. 2008;26: 362–365. pmid:21791389
- 3. Sano D, Fukushi K, Yoshida Y, Omura T.Detection of enteric viruses in municipal sewage sludge by a combination of the enzymatic virus elution method and RT-PCR. Water Res (Oxford) 2003;37: 3490–3498. pmid:12834742
- 4. Veerasingam S, Venkatachalapathy R, Sudhakar S, Raja P, Rajeswari V. Petroleum hydrocarbon concentrations in eight mollusc species along Tamilnadu coast, Bay of Bengal, India. J Environ Sci (China) 2011;23: 1129–1134. pmid:22125905
- 5. Howell ET, Chomicki KM, Kaltenecker G. Patterns in water quality on Canadian shores of Lake Ontario: Correspondence with proximity to land and level of urbanization. J. Great Lakes Res 2012;38, Supplement 4: 32–46.
- 6. Boutet I, Tanguy A, Moraga D. Response of the Pacific oyster Crassostrea gigas to hydrocarbon contamination under experimental conditions. Gene 2014;329: 147–157.
- 7. Chapman RW, Mancia A, Beal M, Veloso A, Rathburn C, Blair A, et al. The transcriptomic responses of the eastern oyster, Crassostrea virginica, to environmental conditions. Mol Ecol 2011;20: 1431–1449. pmid:21426432
- 8. Gueguen Y, Cadoret JP, Flament D, Barreau-Roumiguière C, Girardot AL, Garnier J, et al. Immune gene discovery by expressed sequence tags generated from hemocytes of the bacteria-challenged oyster, Crassostrea gigas. Gene 2003;303: 139. pmid:12559575
- 9. Meistertzheim A-l, Tanguy A, Moraga D, Thébault M-T. Identification of differentially expressed genes of the Pacific oyster Crassostrea gigas exposed to prolonged thermal stress. FEBS J 2007;274: 6392. pmid:18005253
- 10. Hedge LH, Knott NA, Johnston EL. Dredging related metal bioaccumulation in oysters. Mar Poll Bull 2009;58: 832–840. pmid:19261303
- 11. Gall ML, Poore AGB, Johnston EL. A biomonitor as a measure of an ecologically-significant fraction of metals in an industrialized harbour. J Environ Monit 2012;14: 830–838. pmid:22245937
- 12. Fitridge I. The impact and control of biofouling in marine aquaculture: A review. Biofouling 2012;28: 649–669. pmid:22775076
- 13. David E, Tanguy A, Pichavant K, Moraga D. Response of the Pacific oyster Crassostrea gigas to hypoxia exposure under experimental conditions. FEBS J 2005;272: 5635. pmid:16262701
- 14. Huvet A, Herpin A, Dégremont L, Labreuche Y, Samain J-F, Cunningham C. The identification of genes from the oyster Crassostrea gigas that are differentially expressed in progeny exhibiting opposed susceptibility to summer mortality. Gene 2004;343: 211–220. pmid:15563847
- 15. Martín-Gómez L, Villalba A, Abollo E. Identification and expression of immune genes in the flat oyster Ostrea edulis in response to bonamiosis. Gene 2012;492: 81–93. pmid:22085815
- 16. Morga B, Arzul I, Faury N, Segarra A, Chollet B, Renault T. Molecular responses of Ostrea edulis haemocytes to an in vitro infection with Bonamia ostreae. Devel Comp Immunol 2011;35: 323–333. pmid:21073892
- 17. Renault T, Faury N, Barbosa-Solomieu V, Moreau K. Suppression substractive hybridisation (SSH) and real time PCR reveal differential gene expression in the Pacific cupped oyster, Crassostrea gigas, challenged with Ostreid herpesvirus 1. Devel Comp Immunol 2011;35: 725–735. pmid:21371503
- 18. Sussarellu R, Fabioux C, Le Moullac G, Fleury E, Moraga D. Transcriptomic response of the Pacific oyster Crassostrea gigas to hypoxia. Mar Genomics 2010;3: 133–143. pmid:21798207
- 19. Tanguy A, Guo X, Ford SE. Discovery of genes expressed in response to Perkinsus marinus challenge in Eastern (Crassostrea virginica) and Pacific (C. gigas) oysters. Gene 2004;338: 121–131. pmid:15302413
- 20. Tanguy A, Boutet I, Laroche J, Moraga D. Molecular identification and expression study of differentially regulated genes in the Pacific oyster Crassostrea gigas in response to pesticide exposure. FEBS J 2005;272: 390–403. pmid:15654877
- 21. Wang S, Peatman E, Liu H, Bushek D, Ford SE, Kucuktas H, et al. Microarray analysis of gene expression in eastern oyster (Crassostrea virginica) reveals a novel combination of antimicrobial and oxidative stress host responses after dermo (Perkinsus marinus) challenge. Fish Shellfish Immunol 2010;29: 921–929. pmid:20708691
- 22. Kültz D. Molecular and evolutionary basis of the cellular stress response. Ann Rev Physiol 2005;67: 225–5723. pmid:15709958
- 23. Zhang G1, Fang X, Guo X, Li L, Luo R, Xu F, et al. Nature. 2012;490(7418): 49–54.24. pmid:22992520
- 24. Suárez-Ulloa V, Fernández-Tajes J, Manfrin C, Gerdol M, Venier P and Eirín-López JM Bivalve omics: State of the art and potential applications for the biomonitoring of harmful marine compounds. Mar Drugs 2013;11, 4370–4389. pmid:24189277
- 25. Zhao X, Yu H, Kong L, Li Q. Transcriptomic responses to salinity stress in the pacific oyster Crassostrea gigas. PLoS One 2012;7 (9)e46244. pmid:23029449
- 26. de Lorgeril J., Zenagui R., Rosa R.D., Piquemal D., Bachère E. Whole transcriptome profiling of successful immune response to Vibrio infections in the oyster Crassostrea gigas by digital gene expression analysis. PLoS One 2011;6 (8), e23142. pmid:21829707
- 27. McDowell IC, Nikapitiya C, Aguiar D, Lane CE, Istrail S, Gomez-Chiarri M Transcriptome of American oysters, Crassostrea virginica, in response to bacterial challenge: Insights into potential mechanisms of disease resistance. PLoS One 2014;9 (8) e105097. pmid:25122115
- 28. Muralidharan S, Thompson E, Raftos D, Birch G, Haynes PA. Quantitative proteomics of heavy metal stress responses in Sydney rock oysters. Proteomics 2012;12: 906–921. pmid:22539440
- 29. Timmins-Schiffman E, Coffey WD, Hua W, Nunn BL, Dickinson GH, Roberts SB. Shotgun proteomics reveals physiological response to ocean acidification in Crassostrea gigas. PeerJ PrePrints 2014;2:e388v1. pmid:24883249
- 30. Tomanek L Environmental proteomics: Changes in the proteome of marine organisms in response to environmental stress, pollutants, infection, symbiosis, and development. Ann Rev Mar Sci 2011;3:373–99. pmid:21329210
- 31. Sokolova IM Energy-limited tolerance to stress as a conceptual framework to integrate the effects of multiple stressors. Integ Comp Biol 2013;53:597–608. pmid:23615362
- 32. Lannig G, Gisela L Impact of ocean acidification on energy metabolism of oyster, Crassostrea gigas—Changes in metabolic pathways and thermal response. Mar Drugs 2010;8: 2318. pmid:20948910
- 33. Parker LM, Ross PM, O'Connor WA, Borysko L, Raftos DA, Pörtner H-O. Adult exposure influences offspring response to ocean acidification in oysters. Global Change Biol 2012;18: 82–92.
- 34. Jeong W-G, Cho S-M Long-term effect of polycyclic aromatic hydrocarbon on physiological metabolisms of the Pacific oyster, Crassostrea gigas. Aquaculture 2007;265: 343–350.
- 35. Tedesco S, Doyle H, Blasco J, Redmond G, Sheehan D Oxidative stress and toxicity of gold nanoparticles in Mytilus edulis. Aquat Toxicol 2010;100: 178–186. pmid:20382436
- 36. Bocchetti R, Regoli F. Seasonal variability of oxidative biomarkers, lysosomal parameters, metallothioneins and peroxisomal enzymes in the Mediterranean mussel Mytilus galloprovincialis from Adriatic Sea. Chemosphere 2006;65: 913–921. pmid:16678235
- 37. Donaghy L. Reactive oxygen species in unstimulated hemocytes of the pacific oyster Crassostrea gigas: A mitochondrial involvement. PloS One 2012;7(10): e46594. pmid:23056359
- 38. Kuchel RP, Raftos DA. In vitro effects of noradrenaline on Akoya pearl oyster (Pinctada imbricata) haemocytes. Fish Shellfish Immunol 2011;31: 365–372. pmid:21664977
- 39. Aladaileh S, Mohammad MG, Ferrari B, Nair SV, Raftos DA. In vitro effects of noradrenaline on Sydney rock oyster (Saccostrea glomerata) hemocytes. Comp Biochem Physiol Part A 2008;151: 691–697. pmid:18796336
- 40. Labreuche Y, Lambert C, Soudant P, Boulo V, Huvet A. Cellular and molecular hemocyte responses of the Pacific oyster, Crassostrea gigas, following bacterial infection with Vibrio aestuarianus strain 01/32. Microb Infect 2006;8: 2715–2724.
- 41. Park MS, Jo PG, An KW, Cho SH, Kil G-S. mRNA expression of antioxidant enzymes and physiological responses in the Pacific oyster, Crassostrea gigas, exposed to an hypoxic environment. Molluscan Res. 2009;29: 179–184.
- 42. Green TJ, Dixon TJ, Devic E, Adlard RD, Barnes AC. Differential expression of genes encoding anti-oxidant enzymes in Sydney rock oysters, Saccostrea glomerata (Gould) selected for disease resistance. Fish Shellfish Immunol 2009;26: 799–810. pmid:19332130
- 43. Wang Z, Wu Z, Jian J, Lu Y. Cloning and expression of heat shock protein 70 gene in the haemocytes of pearl oyster (Pinctada fucata, Gould 1850) responding to bacterial challenge. Fish Shellfish Immunol 2009;26: 639–645. pmid:19026750
- 44. Ivanina AV, Cherkasov AS, Sokolova IM. Effects of cadmium on cellular protein and glutathione synthesis and expression of stress proteins in eastern oysters, Crassostrea virginica. J Exp Biol 2008;211: 577–586. pmid:18245635
- 45. Farcy E, Voiseux C, Lebel J-M, Fievet B. Seasonal changes in mRNA encoding for cell stress markers in the oyster Crassostrea gigas exposed to radioactive discharges in their natural environment. Sci Total Environ 2007;374: 328–341. pmid:17289119
- 46. Zhang Z, Zhang Q. Molecular cloning, characterization and expression of heat shock protein 70 gene from the oyster Crassostrea hongkongensis responding to thermal stress and exposure of Cu(2) and malachite green. Gene 2012;497: 172–180. pmid:22310388
- 47. Fabbri E, Valbonesi P, Franzellitti S. HSP expression in bivalves. Invert Surv J 2008;5: 135–161.
- 48. Meng J, Zhu Q, Zhang L, Li C, Li L, She Z, et al. Genome and transcriptome analyses provide insight into the euryhaline adaptation mechanism of Crassostrea gigas. PLoS One 2013;8 (3): e58563. pmid:23554902
- 49. Lüchmann KH, Mattos JJ, Siebert MN, Dorrington TS, Toledo-Silva G, Stoco PH, et al. Suppressive subtractive hybridization libraries prepared from the digestive gland of the oyster Crassostrea brasiliana exposed to a diesel fuel water-accommodated fraction. Environ Toxicol Chem 2012;31:1249–53. pmid:22505345