# Deciphering the molecular adaptation of the king scallop (Pecten maximus) to heat stress using transcriptomics and proteomics

Sébastien Artigaud, Joëlle Richard, Michael As Thorne, Romain Lavaud, Jonathan Flye-Sainte-Marie, Fred Jean, Lloyd S. Peck, Melody S. Clark, Vianney Pichereau

<figure><img src="image_1.png" /></figure>

## To cite this version:

Sébastien Artigaud, Joëlle Richard, Michael As Thorne, Romain Lavaud, Jonathan Flye-Sainte-Marie, et al.. Deciphering the molecular adaptation of the king scallop (Pecten maximus) to heat stress using transcriptomics and proteomics. BMC Genomics, 2015, 16, pp.988. 10.1186/s12864-015-2132-x. hal-02552404

HAL Id: hal-02552404

https://hal.science/hal-02552404v1

Submitted on 24 Jun 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Deciphering the molecular adaptation of the king scallop (*Pecten maximus*) to heat stress using transcriptomics and proteomics

Sébastien Artigaud¹, Joëlle Richard¹, Michael AS Thorne², Romain Lavaud¹, Jonathan Flye-Sainte-Marie¹, Fred Jean¹, Lloyd S. Peck², Melody S. Clark² and Vianney Pichereau¹*

**Abstract**

**Background:** The capacity of marine species to survive chronic heat stress underpins their ability to survive warming oceans as a result of climate change. In this study RNA-Seq and 2-DE proteomics were employed to decipher the molecular response of the sub-tidal bivalve *Pecten maximus*, to elevated temperatures.

**Results:** Individuals were maintained at three different temperatures (15, 21 and 25 °C) for 56 days, representing control conditions, maximum environmental temperature and extreme warming, with individuals sampled at seven time points. The scallops thrived at 21 °C, but suffered a reduction in condition at 25 °C. RNA-Seq analyses produced 26,064 assembled contigs, of which 531 were differentially expressed, with putative annotation assigned to 177 transcripts. The proteomic approach identified 24 differentially expressed proteins, with nine identified by mass spectrometry. Network analysis of these results indicated a pivotal role for GAPDH and AP-1 signalling pathways. Data also suggested a remodelling of the cell structure, as revealed by the differential expression of genes involved in the cytoskeleton and cell membrane and a reduction in DNA repair. They also indicated the diversion of energetic metabolism towards the mobilization of lipid energy reserves to fuel the increased metabolic rate at the higher temperature.

**Conclusions:** This work provides preliminary insights into the response of *P. maximus* to chronic heat stress and provides a basis for future studies examining the tipping points and energetic trade-offs of scallop culture in warming oceans.

**Keywords:** Marine biology, Metabolism, DNA repair, Transcription regulation, Apoptosis, Energy reserves

**Background**

Near-shore marine coastal environments are experiencing considerable changes in temperature and oxygen availability due to anthropomorphic influences, which has triggered considerable modifications in biological communities [1]. Therefore, the ability of marine species and populations to cope with environmental change underpins the maintenance of current marine ecosystems. In this respect, temperature is acknowledged as one of the most important environmental threats affecting marine organisms [2, 3]. Indeed in multiple stressor

experiments investigating the interacting effects of temperature with, for example, lowered pH, the dominant effect has been identified as temperature [4–6]. Therefore understanding marine organisms’ responses to thermal acclimation and/or capacity for adaptation is a critical challenge in understanding their ecology and predicting their future distributions. For society, this is of particular interest for economically important aquaculture and fished species, such as *Pecten maximus* used in this study.

Whilst the molecular mechanisms underlying cellular responses to thermal stress have been characterised in a number of model species [7–11], the responses of marine species have largely been evaluated using a limited number of evolutionary conserved stress-responsive genes.

<footnote>* Correspondence: vianney.pichereau@univ-brest.fr
¹Laboratoire des Sciences de l'Environnement Marin, LEMAR UMR 6539 CNRS/UBO/IRD/Ifremer, Université de Brest (UBO), Institut Universitaire Européen de la Mer, Plouzané 29280, France
Full list of author information is available at the end of the article</footnote>
Prime examples of these are the heat shock proteins (HSPs), particularly the inducible 70 kDa form (HSP70). Most of these studies on bivalve molluscs have been very short-term acute thermal shocks. For example, Ivanina et al. [12] showed the overexpression of selected HSPs and metallothionein (MT) in *Crassostrea virginica* exposed for 1 h at 40 °C, and Brun et al. [13, 14] characterized the heat shock response and the acquisition of thermotolerance in selected Pectinidae (*Argopecten irradians*, *Placopecten magellanicus*), and the interplay between hypoxia and temperature stress was recently characterized in *P. maximus* [15, 16]. Lang et al. [17] used a microarray containing 1675 ESTs from *Crassostrea gigas* and *C. virginica* to characterize the transcriptomic response of different families of the Pacific oyster showing contrasting degrees of thermotolerance. Of note, they showed differential expression of genes encoding HSPs, and genes involved in lipid metabolism, protection against bacterial infections and cell structural elements (e.g. collagen) in response to an acute thermal stress. Finally, using a proteomics approach to identify differences in the thermal resilience of the mussel congeners *Mytilus galloprovincialis* and *Mytilus trossulus*, HSPs and proteins combating reactive oxygen species were identified in response to a 1 h thermal challenge [18].

Whilst studies examining short exposure to acute stress can provide interesting insights into affected pathways, they may not reflect the response to a permanent gradual shift such as the predicted increase in sea surface temperature [19]. Short- and long-term exposures can produce very different responses in genes expression profiles [20]. For example, in a recent study, the effect of both acute- (within a day) and long-term- (up to 14 days) exposure to heat stress on the gene expression of *Chlamys farreri* was studied [21]. Although, the question remains as to how long is a long-term challenge. In a study of thermal stress in *C. gigas*, a Q-PCR time course showed a graduated response in the gene expression profiles with a rapid increase at 3–7 days, a decrease at 14 days, which reduced further by days 17–24. By then, expression levels were presumed to have stabilised, and the oysters acclimated to the new temperature [22]. The only other recent study dealing with long-term adaptation to thermal challenge concerned *C. gigas* (3 months, 24 °C) [4]. This study showed the essential role of lipid mobilisation, the mTOR regulatory pathway, and ultimately the induction of apoptosis as a result of chronic elevated heat stress.

Studying organisms' adaptation to changing environments is a real challenge in the field of ecological genomics. In particular, discovery-led transcriptomic and proteomic characterizations of the responses of organisms to environmental changes offer an opportunity to understand the underlying molecular basis for

adaptation. Transcriptomic and proteomic approaches are highly complementary. NGS transcriptomic tools can provide extensive catalogues of genes, even for non-model species and essential reference data for the identification of proteins. Protein production is dependent on the efficiency of transcription and translation of a gene, with as final product, the result of a variety of post-translational modifications, such as phosphorylation. Hence the level of gene expression is not always directly correlated to that of its respective proteins, proteomics is therefore closer to phenotypes than transcriptomics. However, 2-DE based proteomics is naturally biased towards highly abundant proteins, and is often limited by the availability of genomic data, especially for non-model organisms. As a consequence, in spite of its physiological relevance, proteomic strategies classically provide much less data than transcriptomics.

In this study, we used the economically important king scallop *Pecten maximus*, as a model to understand marine species' response to thermal stress. *P. maximus* is naturally distributed along a large East Atlantic latitudinal gradient (from 31 °N to 69 °N) and lives in the subtidal zone down to a depth of 500 m (www.fao.org). It is an economically important species in the UK, France and Spain, where it is considered as a high value product [23]. In this paper, we describe the first genome-wide transcriptome analysis in the king scallop, sampled over a 56-day time course whilst subjected to a thermal challenge at three different temperatures. Coupling both transcriptomic- and proteomic- approaches, provides a complementary system biology view of heat adaptation/acclimation in this species.

## Results

### Culture conditions

These remained relatively stable throughout the 56 days of the experiment. The three temperatures were maintained at: 15.1 ± 0.2 °C as control, 21.4 ± 0.2 °C and 25.2 ± 0.9 °C (see Fig. 1). Salinity, pH and O₂ were maintained at 35.8‰ ± 0.2, 8.1 ± 0.1 and 94 ± 7 %, respectively. There were slight changes to ammonia (that increased from 4.68 µM ± 1.6 to 5.96 µM ± 1.9 and 6.5 µM ± 2.3 at 21 and 25 °C, respectively).

### Physiological response to heat stress

All the scallops used in this study came from the same 2010 production of the Tinduff hatchery, and were held under identical conditions prior the start of the experiment. Hence the animals were similar in size, age and physiological condition. Mortalities occurred in each tank during the experiment (Fig. 2). However, this mortality was not correlated to the increased temperature and more than 80 % of the animals survived in all conditions (Fig. 2b). In contrast, the condition index (CI),
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Temperatures in the three tanks during the experiment</figcaption></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Condition index (a) and survival (b) of *Pecten maximus* during a long-term exposure (56 days) to three different temperatures (15, 21 and 25 °C)</figcaption></figure>

which reflected the physiological condition of the animals, was significantly reduced at the end of the experiment for the animals kept at 25 °C (Fig. 2a).

## Differential gene expression analyses

RNA-Seq of the 22 different transcriptome libraries resulted in a total of 806 million paired end raw reads, which was reduced to 667 million paired end reads after quality clipping. These were assembled into a backbone assembly of 26,064 contigs with an average length of 1011 bp and a median length of 815 bp (fully described in [24]). As the aim of this study was to understand the molecular mechanisms underlying adaptation to long-term exposure to heat stress in *Pecten maximus*, a very stringent strategy was adopted to identify the consistently differentially expressed transcripts. This was

<figure><img src="image_7.png" /><figcaption>**Fig. 3** Representative bi-dimensional gels (pH 4–7, SDS-PAGE 12 %) for *Pecten maximus* mantle proteins maintained at 15 °C (upper gel) and 25 °C (lower gel). Differentially expressed spots are arrowed</figcaption></figure>
performed, by taking into account those displaying a ratio (treatment/ control) of at least 2 in all the last three time points.

1.67 % of the transcriptome (531 transcripts from 26,064 contigs) showed significant differential expression at the elevated temperatures with putative function for 177 transcripts, assigned by Blast sequence similarity searching against the NCBI non-redundant database, using an E-value cut off of $1e^{-10}$ (summarised in Additional file 1: Tables 1 and 2, for 25 °C and 21 °C, respectively). This annotation level of 23.6 % was similar to that obtained for other non-model organisms [25–27]. There was a considerable difference in the response of the scallops at the two temperatures, with only 95 transcripts differentially expressed in animals subjected to 21 °C, compared with 436 transcripts at 25 °C.

## Differentially expressed proteins

Using 2-DE, a mean proteome content of 956 spots was observed. Analyses revealed that 24 protein spots (2.5 % of the observable proteome) were significantly differentially expressed at 25 °C, at one or other of the time points (T4: 27 days or T7: 56 days), of which 12 were up-regulated and 12 were down-regulated (Table 1, 3). All of these 24 protein spots were subjected to trypsin digestion and Matrix-Assisted Laser Desorption Ionization Time-Of-Flight tandem mass spectrometry (MALDI TOF-TOF) analysis. This led to the identification of nine proteins using sequence similarity searching (Table 2). Of these nine putatively annotated proteins, five were identified from ESTs sequenced in the RNAseq experiments described in this study, and the other four proteins from the scallop hemocytes ESTs [28] used to build our database. This annotation level of 37 % was similar to proteomic studies in other non-model species [26, 29].

## Network analysis

All of the 177 transcripts putatively annotated using the initial Blast sequence similarity searching were used in further searches against the UniProt human database and the results (short names of genes) were entered into the STRING program. In all, 132 genes had homologs in human, and network analysis highlighted relationships between 35 of them (Fig. 4), among which only one, titin, was relied to the 21 °C condition, the other sequences were all from the 25 °C samplings. The main result showed a central node around glyceraldehyde-3-phosphate dehydrogenase (GAPDH: 10 interactants), which showed weaker links to nodes centered around the c-JUN pathway (8 interactants), calreticulin (CALR: 7 interactants), 90 kDa heat shock protein (HSP90: 5 interactants) and ankyrin (ANK1:3 interactants). Three other independent networks were also identified centering on the cell division cycle and apoptosis regulator

<table><caption>**Table 1** Values of Log₂ Fold Change for spots differentially expressed between animals maintained at 15 °C and animals maintained at 25 °C, after 27 days (T4) and 56 days (T5)</caption><thead><tr><th></th><th colspan="2">Log₂ FC (25 °C/15 °C)</th></tr></thead><tbody><tr><td>Spot</td><td>T4</td><td>T5</td></tr><tr><td>HiT-1</td><td>2.22</td><td>0.97</td></tr><tr><td>HiT-2</td><td>2.07</td><td>2.85</td></tr><tr><td>HiT-3</td><td>1.10</td><td>0.33</td></tr><tr><td>HiT-4</td><td>0.95</td><td>1.20</td></tr><tr><td>HiT-5</td><td>0.62</td><td>1.42</td></tr><tr><td>HiT-6</td><td>0.43</td><td>1.22</td></tr><tr><td>HiT-7</td><td>0.36</td><td>1.80</td></tr><tr><td>HiT-8</td><td>0.27</td><td>1.50</td></tr><tr><td>HiT-9</td><td>0.11</td><td>1.16</td></tr><tr><td>HiT-10</td><td>-0.50</td><td>1.11</td></tr><tr><td>HiT-11</td><td>-0.41</td><td>1.23</td></tr><tr><td>HiT-12</td><td>-1.14</td><td>1.42</td></tr><tr><td>LoT-1</td><td>-1.91</td><td>1.04</td></tr><tr><td>LoT-2</td><td>-2.31</td><td>0.98</td></tr><tr><td>LoT-3</td><td>-1.62</td><td>0.87</td></tr><tr><td>LoT-4</td><td>-2.14</td><td>0.80</td></tr><tr><td>LoT-5</td><td>-1.43</td><td>0.74</td></tr><tr><td>LoT-6</td><td>-1.77</td><td>0.38</td></tr><tr><td>LoT-7</td><td>-1.75</td><td>0.14</td></tr><tr><td>LoT-8</td><td>-1.50</td><td>-0.23</td></tr><tr><td>LoT-9</td><td>-1.37</td><td>-0.61</td></tr><tr><td>LoT-10</td><td>-0.14</td><td>-1.00</td></tr><tr><td>LoT-11</td><td>-0.30</td><td>-1.01</td></tr><tr><td>LoT-12</td><td>-0.33</td><td>-1.77</td></tr></tbody></table>

(CCAR1), arylformamidase (AFMID) and mucin (MUCIN/MUC5B) (for full list of gene symbols and names see Additional file 1).

## Discussion

In this paper, we used both transcriptomic and proteomic analysis of mantle tissue to characterize the acclimation of *P. maximus* to 56-days exposure to elevated thermal stress. The sea surface temperatures in the Bay of Brest vary annually between 10 and 18 °C (data from the Somlit-Brest, Coastal time-series station; 10 m depth), hence the choice of 15 °C as the control temperature, as this is the temperature in the Bay when the animals were brought into the aquaria. 21 °C was chosen as this represented the maximum temperature this species experiences in the region, with 25 °C as a more extreme, stressful temperature.

This range of temperatures is within those tested previously on *P. maximus* [15, 16, 30, 31]. Both the physiological data (condition index) and expression profiles in
<table><caption>**Table 2** List of *Pecten maximus* mantle tissue proteins identified by MS/MS whose abundance change between animals maintained at 15 °C and animals maintained at 25 °C (moderate t-test paired-comparison. fdr <0.1. absolute fold change >2)</caption><thead><tr><th>Spot</th><th>Blast hit name</th><th>Acc. Nr</th><th>EST sequence name</th><th>Sequenced peptides</th><th>Peaks score</th></tr></thead><tbody><tr><td>HiT-1</td><td>calumenin-B-like (CALUB)</td><td>XP_003704568</td><td>Contig22840</td><td>TEFMYFVHPEEGK<br>DVTVCEYTDR</td><td>83.6</td></tr><tr><td>HiT-2</td><td>calumenin-like isoform X3 (CALU)</td><td>XP_005110308</td><td>Contig_ Reproseed325</td><td>DVLVLEYTDR<br>DNAEEFEHATEQESK<br>CSWPNYVR</td><td>98.0</td></tr><tr><td>HiT-9</td><td>No hit</td><td>-</td><td>Contig_ Reproseed34795</td><td>TFVFPPADTQKPVITGIM(+15.99)R<br>GESYIHVK<br>TFVFPPADTQKPVITGIMR<br>TFQM(+15.99)FDALQYIEQGNM(+15.99)IQGR</td><td>88.0</td></tr><tr><td>LoT-2</td><td>Putative phosphoglycerate mutase (PGAM)</td><td>EKC26210</td><td>Contig_ Reproseed34491</td><td>VLISAHGNSLR<br>SYDVPPPAREDGDER<br>YAHKDASVVPR</td><td>94.3</td></tr><tr><td>LoT-3</td><td>glyceraldehyde 3-phosphate dehydrogenase (GAPDH)</td><td>XP_002434347</td><td>scallop_ rep_c46254</td><td>SSXFDANAGHAYNNNFVK<br>VPVPDVSVWDHTC(+57.02)R<br>VWSQDFC(+57.02)GDSR</td><td>98.3</td></tr><tr><td>LoT-6</td><td>gelsolin-like protein 2-like isoform X1 (GSN)</td><td>XP_005100380</td><td>Contig28884</td><td>AWDGAGQEPGTQFWR<br>QSWQVGNR<br>NSGNSGDVYXPDGGR</td><td>85.7</td></tr><tr><td>LoT-7</td><td>Kynurenine-oxoglutarate transaminase 3 (CCBL2)</td><td>EKC20610</td><td>Contig34767</td><td>NLGENFLR<br>AVNSDNANWAQYAR<br>GQVPDDGSDDPYDYK<br>IASLPGM(+15.99)WDR</td><td>97.3</td></tr><tr><td>LoT-9</td><td>peptidyl-prolyl cis-trans isomerase B precursor (PPIB)</td><td>NP_990792</td><td>Contig2622332</td><td>DFMVQGGDFSEGDGTGSK<br>VSEGMDVVR<br>KLENTEVDIENR<br>YFADENFK</td><td>98.4</td></tr><tr><td>LoT-11</td><td>Collagen alpha-6(VI) chain (COL6A6)</td><td>EKC27115</td><td>Contig35480</td><td>FNTNADQASVTEDVDDRR<br>DDAEHVGII LTDTGTTNPGR<br>REFNLDSYTKPAEVEQGLSNVR<br>NFQLELSFVK<br>ALTLLLDEGFTVR<br>VSNFDIGSDNTR</td><td>99.0</td></tr></tbody></table>

this experiment showed significant changes in animals kept at 25 °C compared with 21 °C. This is coherent with previous studies which show that *P. maximus* can rapidly acclimate (in timescales of days to weeks) to water temperatures below 23 °C [21, 30, 31]. The paucity of differential expression and stability of the condition index (Additional file 1: Table 2, Fig. 2) of the 21 °C animals indicates that these animals have successfully acclimated to this slightly elevated temperature [21]. This does not appear to be the case with the 25 °C animals where there is a reduction in condition index and this is the temperature which exhibits 82 % of the differentially

expressed transcripts. Recently published experiments have shown that prolonged exposure, in excess of a week, to 25 °C can take *P. maximus* beyond their optimal thermal window with a switch towards anaerobic metabolism [15, 16]. These analyses also showed that at 25 °C *P. maximus* reduced respiration rates, indicating that metabolism was shutting down and also significant mortalities occurred when temperature was combined with hypoxia in a multistressor experiment [15, 16]. In the experiment described here there were no significant mortalities reported at 25 °C and thus network analysis of the differentially expressed transcripts and proteins at
<figure><img src="image_8.png" /></figure>

this temperature can provide insights into cellular re-
sponses to chronic heat stress. It was possible to parti-
tion these results into different nodes and functional
categories, as described below.

GAPDH as a pivotal protein in thermal acclimation

The network analysis placed GAPDH (glyceraldehyde-
3-phosphate dehydrogenase) as a central node, linking
several other pathways. GAPDH is a key enzyme in
the glycolytic pathway, catalysing the conversion of
glyceraldehyde-3-phosphate (G3P) to 1,3-biphospho-
glycerate in the presence of NAD+ and inorganic
phosphate. Therefore the presence of this transcript
might not be unexpected due to the thermal depend-
ency of biochemical reactions. Higher temperatures
would cause an increase in the animals’ metabolic
rate with glycolysis providing the additional energy
needed, at least in the short term [32]. Indeed, evi-
dence of up-regulation of glycolytic pathways has
been identified in other molluscs thermal acclimation

studies [4, 18]. Data from our study indicate that this
scenario is unlikely, as there is a reduction in con-
dition index (Fig. 2) and putative lipid metabolism
transcripts are being up-regulated (see below). This
indicates that carbohydrate sources have been
exhausted and that stored energy reserves are being
mobilised. Thus GAPDH may be playing a role in al-
ternative biochemical pathways. Until recently, this
enzyme was purely considered as a house-keeping
gene. However recently, new extra-glycolytic roles for
GAPDH in a diverse range of cellular processes have
been discovered, such as binding nucleic acids, espe-
cially certain mRNAs and telomeric DNA [33]. It may
also function as a transcription factor interacting with
RNA pol II [34], and acts as a signaling protein, by
regulating the activity of the mTOR signaling pathway
through sequestering the GTPase Rheb [35]. Finally,
several links between GAPDH and apoptosis have also
been established. This is via the interaction of
GADPH with the ubiquitin ligase Siah1, but it also
plays a role in p53 and caspase-3 dependent apoptosis
[33, 36].

GAPDH was only identified in the proteomic ana-
lysis and was not significantly up-regulated in the
transcriptional profiling. However, it should be noted
that GAPDH acts in cells as a homotetramer, and is
known to be primarily regulated through post-
translational modifications (e.g. phosphorylation; [37]).
Whilst GAPDH was centrally positioned within the
network, there were several links to a smaller node
centered on the pleiotropic transcription factor AP1
(activating protein 1).

A major role for signaling through the AP-1 pathway ?

Numerous genes encoding nuclear proteins, or involved
in protein trafficking in the nucleus (Additional file 1:
Tables 1 and 2), were found differentially expressed in
this study, of which, AP1 (JUN) was highlighted by the
network analysis. AP-1 is known to regulate a number of
major physiological processes including cell fate. This
may either be proliferation or apoptosis, with the cellular
context critical for determining which pathway is acti-
vated [38, 39]. AP-1 is not a single protein, but a com-
plex of different dimers from several groups of the Jun
bZIP (basic region leucine zipper) domain proteins [40].
In this study, putative transcripts were identified for
other potentially related CREB/ATF bZIP transcription
factors. One of these transcripts showed high sequence
similarity to BATF (Basic leucine zipper transcription
factor, ATF-like), whilst the other matched a hypothet-
ical protein containing a bZIP domain, thus suggesting
that it belonged to the same class of transcriptional
regulators.

Like all bZIP transcription factors, proteins within the
AP-1 complex must dimerize before binding to their
DNA target sites. However, each of the proteins leading
to the AP-1 complex are differentially expressed and
regulated, so every cell contains different mixtures of
AP-1 dimers which can result in subtle differences in
function [39, 40]. This is shown through the network
analysis with links to ten other genes, known to be regu-
lated by this transcriptional factor or interacting with it,
including BATF, but also TLR4, GNRH, POLR3C, CA2,
HSP90, TRAF3 and GAPDH. Whilst the network was
built on interactions of human genes, in other species,
AP1 has been found to regulate other genes such as gel-
solin, collagen, phospholipases and DNA polymerases,
some of which are differentially expressed in this study
(additional file 1: Tables 1 and 2). Therefore it is entirely
possible that some of these other transcripts are in fact
members of the AP1 regulation pathway in *P. maximus*.

Additional links to AP1 were identified via several com-
ponents of a multi subunit autoregulatory ribonucleopro-
tein complex, potentially involved in translationally coupled

mRNA turnover [41]. This complex may associate with the
mCRD domain of c-FOS (a component of AP-1), and in-
cluded, PAIP1, HNRPD, CSDE1 and SYNCRIP. Putative
homologues of the latter two genes, SYNCRIP and CSDE1,
were found to be up-regulated at 25 °C in this study.
Interestingly, SYNCRIP was also shown to interact
with CAPRIN1 and with GAPDH in our network ana-
lysis. Of the links identified in the AP1 network, that
with TRAF3 (tumor necrosis factor receptor (TNFR)-
associated factor 3) is of particular interest. TRAF3
was one of the first identified TRAFs, initially isolated
by virtue of its binding to the cytoplasmic domain of
the TNFR-superfamily (TNFR-SF) member CD40 [42].
TRAF3 is a pivotal protein in cell signaling interact-
ing with numerous receptors, such as TLR, RLR,
CD40 etc. Interestingly, its action may require ubiqui-
tination [43], and may lead to the activation of the
JNK pathway and the negative regulation of NF-kappa
B, which is involved in the immune response and cell
atrophy [44]. AP-1 is major target of the JNK MAP
kinase pathway and it is known to regulate a number
of major physiological processes including cell fate. It
is likely in this experiment that the putative AP-1
(JUN) transcript was involved in apoptosis, as the condi-
tion index of the animals was decreasing over time.

A similar link between members of the glycolytic path-
way and apoptosis was identified in a recent transcrip-
tomic study of chronic thermal challenge in the oyster
Crassostrea gigas [4]. In that study glucose-6-phosphate
translocase was identified in the transcriptional profiles
along with several gene members involved in the mTOR
signalling pathway and apoptosis [4]. Like the study de-
scribed here, these profiles were associated with animals
with reduced condition indexes and physiological condi-
tion, which were struggling to cope with exposure to ele-
vated temperatures.

Heat shock proteins

In our transcriptome data genes encoding both the
70 kDa (HSP70) and 90 kDa (HSP90) genes and a small
HSP40 were up-regulated at 25 °C, along with several re-
lated genes, including an activator of HSP90 (AHSA1) and
other molecular chaperones e.g. PPIB, PDIA, CALR (Add-
itional file 1: Table 1). Proteins involved in proteins recyc-
ling (linked to ubiquitin dependent protein recycling)
were also found. This indicates a cellular requirement for
enhanced protein folding capacities and removal of mis-
folded proteins, which are classic indicators of cell stress.
To date, the only universal response to heat stress docu-
mented is the production of heat shock proteins (HSPs),
often called the heat shock response (HSR) [45] and these
genes are often monitored as a measure of environmental
stress in molluscs [4, 13, 14, 17, 18, 22]. The majority of
the experiments in scallops to date have used HSPs to test
thermal tolerance, but these genes have largely been identified in acute exposures [13, 14, 21]. What is particularly interesting with this *P. maximus* experiment, is that putative HSPs and chaperone transcripts were still expressed after a chronic exposure and production of these proteins is known to be energetically costly [46].

However, the heat shock response is complex and dependent on the duration, type and intensity of the stress [47] and these proteins are not just involved in response to heat shock. HSPs have a wider role in relation to apoptosis and the inflammatory process [48]. Experimental induction of HSPs has been shown, in a variety of mammalian tissues, to inhibit apoptosis and considerably enhance cell survival [49, 50]. Thus, although many molluscs show expansions of heat shock gene family members [51, 52], these may not all be required for response to thermal stress. For example, the Pacific oyster may have up to 88 HSP70 genes in the genome, a study which exposed this species to elevated temperatures for 3 months, only identified a single HSP7012B and HSP60 gene in the expression data [4]. The poor numbers of HSPs identified as up-regulated in the oyster study and the continued expression of HSPs in our study might reflect the functional diversity of HSPs in molluscs.

## Long term-exposure to heat stress and apoptosis

The previous results, detailed above, on the activation of the AP-1 signalling pathway and up-regulation of HSPs, suggest a tension between cell renewal and apoptosis in the mantle tissue of 25 °C animals. The AP-1 complex can either promote cell proliferation or apoptosis and HSPs can have an anti-apoptotic function. In the transcriptome data a number of differentially expressed genes encoded apoptosis-related proteins. In particular, 4 transcripts were putatively annotated as inhibitors of apoptosis (3 were annotated as BIRC7, and 1 as PIAP) (Additional file 1: Table 1). The BIRC7 proteins repress apoptosis through the inhibition of CASP3, CASP7 and CASP9, and through their E3 ubiquitin-protein ligase activity. In addition, the major regulator of apoptosis (CCAR1) was also found up-regulated during heat stress. Network analysis showed interactions with CROP, which encodes a cisplatin resistance associated over-expressed protein (cisplatin induces cell death through induction of apoptosis [53], CSTF3 (Cleavage stimulation factor 77 kDa), and SMC1A (Structural maintenance of chromosomes protein 1A). A number of genes related to the ubiquitination and protein degradation were also up-regulated, which would be required if mis-folded proteins were accumulating. So the fact that after 56 days of being cultured at 25 °C HSPs and anti-apoptotic transcripts were significantly up-regulated in *P. maximus* indicates that the scallops were probably subject to chronic stress and having to prioritise the unfolded protein response either/and regulate apoptosis. The latter of

which would certainly agree with the results showing a reduction in body mass, the decreased condition index and the mobilisation of energy reserves, as described below. This agrees with previous long-term studies in molluscs. In these, if the temperature increase was within the acclimation capacity of the animals, little response was seen in terms of gene expression profiles [21, 22], whereas lack of acclimation was demonstrated by a measurable cellular stress response [4, 17, 54].

## Heat stress and mobilisation of stored energy reserves

The transcription profiles suggest that the animals at 25 °C were struggling to consume and metabolise sufficient food to support life at the higher temperature. This agrees with the condition index (CI) data, which show a steady decrease in the CI of the 25 °C animals with time. When animals are subjected to increased temperatures, their metabolic rate also increases [55]. Our results suggest that they boost energy production to deal with this and also the increased energetic demand associated with mounting a robust stress response. Numerous studies have revealed the up-regulation of metabolism genes in hyperthermally stressed organisms such as the yeast *Saccharomyces cerevisiae* [56], the goby *Gillichthys mirabilis* [57], the mussels *Mytilus californianus* [58] and *M. trosulus* [18]. In these studies, most of genes related to energy metabolism, which were up-regulated, were related to glycolysis and the subsequent respiratory metabolism. This was not the case with the transcriptomic data described here. Two proteins related to the glycolysis (GAPDH and PGAM) were found to have their abundance decreased at 25 °C, while several genes associated with lipid metabolism were up-regulated, i.e. APMAP, PLA2G16 (2 contigs), AMACR, and ACOT4, suggesting a switch towards the mobilization of stored lipids. A similar observation was found in one of the few studies dealing with the adaptation to prolonged exposure to heat stress and low pH in a mollusc bivalves, the Pacific oyster [4].

## Remodelling of cell membranes and cytoskeleton

Changes to the basic cell structure were indicated by the differential expression of several transcripts, at both 21 °C and 25 °C (Additional file 1: Tables 1 and 2). It is well known that temperature affects the lipid constitution of cell membranes [59] as it is essential to maintain the membrane lipids in the liquid crystalline phase to maintain function. This is mainly achieved by altering the relative amounts of unsaturated fatty acids (UFAs) within the membranes, a process termed homeoviscous adaptation (HVA) [60]. In animals, the main enzyme involved in fatty acid desaturation is the enzyme Δ9-acyl-CoA desaturase. In our study, we showed that the gene encoding this enzyme (SCD) was strongly down-regulated in the mantle
cells of scallops subjected to 25 °C. This would most probably result in a decreased level of unsaturation in membrane fatty acids, thus reducing the fluidity of membranes at the higher temperatures.

Cytoskeleton transcripts also showed differential expression at both 21 °C and 25 °C (Additional file 1: Tables 1 and 2). At the proteomic level several proteins involved in structure were also identified as decreased in abundance at 25 °C (GSN, COL6A6, Table 2). Beside from their structural role, cytoskeleton proteins have been shown to have a protective role in response to reactive oxygen species produced as a result of cell stress [61]. However a wider range of cytoskeletal proteins, such as tubulin, were implicated in studies investigating thermal stress in mussels, with a suggestion that this was due to the changing requirements for cell proliferation [18, 54, 62].

## DNA damage and repair

Our transcriptomic study revealed the differential regulation of numerous genes encoding proteins putatively involved in DNA repair and replication, both at 21 and 25 °C. In particular, we observed the down-regulation of not less than five different contigs annotated as Uracil DNA glycosylase (UDG). Numerous studies have indicated that environmental stress can affect genome stability in eukaryotes [63]. For example, in mammalian cells, hypoxia and starvation can suppress error-free DNA repair pathways (e.g. mismatch repair and homologous recombination) and cause an increase in mutagenesis [64–67] and Shor et al. [68] recently discovered that the yeast environmental stress response regulates mutagenesis. Base excision repair (BER) is one of the most important mechanism protecting cells from the deleterious effects of endogenous DNA damage induced by hydrolysis, reactive oxygen species or other intracellular metabolites that modify the DNA base structure. The pivotal enzymes involved in BER are DNA glycosylases, which remove different types of modified or damaged bases by cleavage of the N-glycosidic bond between the base and the 2-deoxyribose moieties of the nucleotide residues. Different DNA glycosylases, including UDG, remove different kinds of damage, and the specificity of the repair pathway is determined by the type of glycosylase involved. Uracil appears in DNA as a result of cytosine deamination, one of the most common DNA spontaneous modifications (that naturally occur in genomes at a rate of around 100 lesions per cell per day). Since this lesion is directly mutagenic, producing C→T transition, the production of UDG is essential for all living organisms. Therefore, the down-regulation of not less than five different contigs annotated as UDG in this study indicates that these cells may experience increased

mutation rates, which could ultimately lead to increased cell death.

## Acute stress versus time limited thermal tolerance versus long term acclimation

Laboratory-based thermal tolerance trials can provide valuable data on species' capacities to respond to changing environmental conditions, particularly their ability to acclimate [69]. If these studies include molecular analyses, then the cellular responses underlying sub-lethal exposures can identify which biochemical pathways are up regulated and therefore help predict energetic trade-offs and long-term survival [70]. To date, the majority of such multidisciplinary thermal trials on scallop species have involved acute exposures with the consequential up-regulation of heat shock proteins [13, 14]. The one longer term study of 14 days showed an immediate “classical” stress response, but by 14 days the animals appeared to have fully acclimated to the new temperature with gene expression levels returned to those of the controls [21]. This expression profile is typical of temperate species where the thermal stress applied is within the normal range experienced by that species and not extreme [22]. In this study we were primarily interested in long term effects and therefore did not examine molecular responses at acute timescales. The transcriptional profiles reflect this, with a lack of up-regulation of antioxidants which often accompany acute responses [5, 18, 21, 22]. The lack of response and stable condition index of the 21 °C animals mirrors the 14 days study [21] on the Zhikong scallop *Chlamys fareri* which demonstrated successful acclimation. It should be noted that adult acclimation does not necessarily guarantee stable reproductive capacity and successful recruitment of juveniles, longer term, generational studies would be required to test this. In contrast, the 25 °C animals showed many similarities in their expression profiles with a chronic 3 month exposure to elevated temperature in *C. gigas*, which in spite of a lack of animal mortality, showed a decrease in condition index and induction of apoptosis and mobilisation of stored energy reserves [4]. If both the experiment described here and the *C. gigas* study had been extended by several months, we would predict significant mortalities, as the animals are merely resisting the higher temperatures and are exhibiting time limited thermal tolerance.

## Conclusions

The aim of this study was to assess the molecular mechanisms underlying responses of the king scallop to prolonged heat stress. It was clear from both the condition indices and the molecular results that the scallops were able to support long-term culture at 21 °C, but 25 °C
resulted in chronic cell stress, a reduction in animal mass and potential viability. Combining RNAseq based transcriptomics and proteomics approaches we identified differential expression in a number of different pathways. Whatever the temperature, the classical heat shock response appeared reduced. Rather, our main results suggested a remodelling of the cell structure as revealed by the differential expression of many genes involved in the cytoskeleton. Interestingly, several genes involved in DNA repair appeared down-regulated, suggesting that cells turned towards a mutator state in response to stress. We also showed, in animals maintained at 25 °C, modifications of cell membrane properties (through the modulation of fatty acids unsaturation) and an apparent mobilization of lipid energy reserves. Finally, our results also suggested a central role for the AP-1 signaling pathway and a potential tension with regard to apoptosis, and led us to propose a pivotal role for GAPDH, a common glycolytic enzyme which has recently been shown to play a critical role in cellular functioning. These are clearly key pathways to investigate in future studies into the thermal resilience of the king scallop *P. maximus*.

# Methods

## Biological material

Juveniles of *Pecten maximus* (1 year old; average length 34.0 mm ± 4.1 mm S.D.) used in this study came from the same cohort spawned at the Tinduff hatchery (Bay of Brest, France) in Spring 2010. They were then raised in small cages at Sainte-Anne du Portzic in the Bay of Brest for 1 year before being transferred to tanks at the Argenton Shellfish Laboratory (Brittany, France) of the French Research Institute for Exploitation of the Sea (IFREMER). Thus, all animals were similar in size, age and physiological condition.

## Experimental culture conditions

The scallops were randomly divided into three batches of 440 individuals placed into 250 L tanks. In each tank, the 440 individuals were sub-divided into two groups and placed in plastic cages without sediment. Each tank was supplied with 1 µm filtered, UV sterilized seawater at a controlled temperature of 15 °C (the temperature in the Bay of Brest at the time) with a flow rate of 80 L.h⁻¹. A 12:12 photoperiod regime was used. Complete water changes were performed once a week. The animals were continuously fed with a cultured mix of microalgae (50 % *Isochrysis affinis* galbana (Tahitian strain *T.iso*) and 50 % *Chaetoceros calcitrans*) via a peristaltic pump delivering directly into the inflow water supply and mixed in the tanks by aquarium pumps. The supply of microalgae to each tank was regulated twice a day by counting the microalgae at the inlet and outlet of the

tank with a Coulter counter to ensure that the outflow water contained 1.5 × 10⁶ µm³.mL⁻¹ of microalgae [71]. This ensured that the scallops were fed *ad libitum* irrespective of the temperature and the number of individuals in the tank. The scallops were initially left at the ambient temperature of 15.1 ± 0.2 °C for 16 days to allow the animals to acclimatise to the rearing conditions and to identify any damaged or unhealthy individuals prior to experimentation. Then, one tank was maintained as a control at approximately 15 °C (15.1 ± 0.2 °C), while the temperature was increased at 1 °C/day to reach 21 °C (21.4 ± 0.2 °C) and 25 °C (25.2 ± 0.9 °C) in the two other tanks. During the experiment, the temperature was recorded every 15 min in the three tanks using button temperature data loggers (Signatrol SL52T Button Data Logger). Additional environmental parameters were also monitored: salinity, pH, O₂ and ammonia. The condition index (CI) was estimated according to the Lucas and Beninger [72] method. For this, 15 to 30 animals were sampled every week. Soft tissues were dried at 75 °C during 48 h and shells were air-dried. Tissues and shells were weighted, and the condition index (CI) was calculated following the formula:

$$CI = (Soft Tissues Dry Weight/Shell Dry Weight) \times 100$$

## Molecular analyses

All molecular analyses were carried out using the mantle tissue. The choice of mantle as the experimental tissue was influenced by several factors. A number of transcriptome studies have been published for this organ, which enhanced the annotation potential of transcripts [25, 27, 73]. Also, this organ has multiple functions including shell formation, secretion of the ligament and sensorial activities and hence is very responsive to external stimuli [74] and transcription profiles of the mantle can act as an effective proxy of early whole-animal response [4]. At the beginning of the experiment (T0), 10 (four for transcriptomic and six for proteomic) individuals were sampled and then, 10 individuals were sampled in each temperature treatment after 3 days (T1), 14 days (T2), 21 days (T3), 27 days (T4) and 56 days (T5). The scallops were quickly dissected and mantle tissue was flash frozen in liquid nitrogen and stored at -80 °C until further analysis.

## RNA extraction and sequencing

Total RNA was extracted from the mantle tissue of 4 individuals per condition at each time point using TRI Reagent® Solution (Life Technologies) according to manufacturer's instructions (22 libraries in total). RNA from the four individuals was then pooled for each time point separately. RNA quality and concentration were
determined using an Agilent 2100 RNA Nanochip (Agilent, Santa Clara, CA, USA) and by NanoDrop ND-1000 Spectrophotometer (260/280 nm, NanoDrop Technologies, Wilmington, DE, USA), respectively. Samples were then treated with a DNA-free Kit (Ambion, Austin, TX) to remove genomic DNA, and reverse transcription was carried out using 500 ng of total RNA, random hexamers, and MMLV reverse transcriptase (Promega) according to manufacturer's protocols. The production of Illumina libraries for mRNA-seq and the transcriptome sequencing using Illumina HiSeq™ 2500 sequencing machine (HiSeq 100 pair-ends) was conducted by the Genome Analysis Centre (Norwich, UK), as described in [24].

### RNA-Seq data sets

The RNA libraries yielded a total of 806 million paired end reads. Raw reads were filtered and trimmed using the FASTX-toolkit (Version 0.0.13 from Assaf Gordon Hannon lab) and rRNA contamination removed using ribopicker [75] and cutadapt (Version 1.1; [76]), with a final quality check performed using fastQC (Version 0.10.0; http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/). The contigs were assembled using SOAPdenovo [77] and a kmer size of 89 was used to construct the initial de novo transcriptome assembly, resulting in 1,311,367 contigs, which were further assembled with CAP3 [78]. Redundancy was determined by self-Blasting and the use of CD-HIT (95 % similarity [79]). *De novo* assembly produced a total number of 26,064 contigs [24]. The transcripts were then processed through the Blast2GO pipeline to produce putative annotations and functional classifications based on Blastx results against the GenBank NR database release 190. To enable differential genes expression analyses, normalization was carried out by dividing counts by library size using TPM (Transcripts per million). Ratios were calculated for each time point with 25 °C versus 15 °C treatment and 21 °C versus 15 °C. Only the transcripts with a ratio higher than two successfully at T3, T4 and T5 were kept for further analyses.

### Protein extraction and separation

Mantle tissues from six individuals from the 25 and 15 °C treatment at the two last time points (T4 and T5) were crushed with a mixer mill (MM400, RETSCH, Haan, Germany). 100 mg of this powder was homogenized in 100 mM Tris-HCl (pH 6.8) with 1 % of Protease inhibitor mix (GE Healthcare), centrifuged (4 °C, 50,000 g, 5 min) and supernatants removed to new tubes. Nucleic acids were then removed (nuclease mix, GE Healthcare, following manufacturer's instructions). Samples were precipitated at 4 °C using TCA 20 % (1/

1:v/v, overnight). After centrifugation (4 °C, 20,000 g, 30 min), pellets were washed with acetone 70 % and re-suspended in urea/thiourea buffer (2 M thiourea, 7 M urea, 4 % CHAPS, 1 % DTT) containing 1 % IPG (pH 4–7, GE Healthcare). Protein concentrations were determined using a modified Bradford assay [80], and all samples were adjusted to 400 µg of proteins in 250 µl. Prior to isoelectric focusing, IPG strips (pH 4–7, 13 cm, GE Healthcare) were passively rehydrated with 250 µl protein solution in wells for 14 h. Isoelectric focusing was conducted using the following protocol: 250 V for 15 min, 500 V for 2 h, gradient voltage increase to 1000 V for 1 h, gradient voltage increase to 8000 V for 2 h 30, 8000 V for 3 h, and reduced to 500 V (Ettan IPGphor3, GE Healthcare). To prepare for second dimension SDS-PAGE electrophoresis, strips were incubated in equilibration buffer (50 mM Tris-HCl pH 8.8, 6 M urea, 30 % glycerol, 2 % SDS and 0.002 % Bromophenol Blue) for two 15-min intervals, first with 1 g.l⁻¹ dithiothreitol and then with 48 g.l⁻¹ iodoacetamide. IPG strips were then placed on top of 12 % polyacrylamide gels, which were run in 10 °C thermo-regulated device (SE 600 Ruby, Amersham Biosciences) at 10 mA per gel for 1 h and then 30 mA per gel until complete migration. Gels were subsequently stained with Coomassie Blue (PhastGel, GE Healthcare) and destained in 30 % methanol and 7 % acetic acid. The resulting gels were scanned with a transparency scanner (Epson Perfection V700) in gray scale with 16-bit depth and a resolution of 400 dpi.

### Differential proteomics analysis

Images were aligned and spots were detected and quantified using the Progenesis SameSpots software (version 3.3, Nonlinear Dynamics) applying the automated algorithm. All detected spots were manually carefully checked and artifact spots were removed. Data were exported as volume raw values and statistical analyzes were conducted in R [81] using the packages prot2D [82] and Limma [83]. Data were normalized (quantile normalization) and the samples from the 25 °C treatment were paired-compared using moderated *t*-test to the samples from 15 °C treatment (control) for each sampling date (with 6 replicates per group). A global correction by false discovery rate (fdr; [84]) was used, in order to take into account multiple comparisons issues and paired-comparison correction. Spots with an FDR threshold lower than 0.1 and an absolute fold change superior to two were considered as differentially expressed.

### Mass spectrometry and protein identification

Proteins for which the abundance was significantly changed between 25 °C and 15 °C were excised from
gels and prepared for mass spectrometry (MS) analysis, essentially as described in [85]. Briefly, gel pieces were first washed in 50 mM ammonium bicarbonate (BICAM), dehydrated in 100 % acetonitrile (ACN), vacuum-dried, and then rehydrated with BICAM containing 0.5 µg sequencing grade porcine trypsin (Promega). After overnight incubation at 37 °C, peptides were extracted from the gels by alternatively washing with 50 mM BICAM and ACN, and with 5 % formic acid and ACN. Between each step, the supernatants were pooled, and finally concentrated by evaporation using a centrifugal evaporator (SpeedVac). Samples were resuspended in trifluoroacetic acid (TFA; 0.1 % in water), and mixed with the α-cyano-4-hydroxycinnamic acid (HCCA, 10 mg·ml⁻¹ of a ACN/TFA/water (60/4/36:v/v/v) solution), and spotted on a polished steel target using the dried droplet method. Peptides were then analyzed by Matrix-Assisted Laser Desorption Ionization Time-Of-Flight tandem mass spectrometry (MALDI TOF-TOF) in positive ion reflector mode, using an Autoflex III (Bruker Daltonics) mass spectrometer. The FlexControl software (v3.0, Bruker Daltonics) was set up to acquire successively PMF spectra and MS/MS from the dominant peaks. Mass spectra were analyzed with FlexAnalysis (v 3.0; Bruker Daltonics) by applying the following conditions: TopHat algorithm for baseline subtraction, Savitzky-Golay analysis for smoothing (0.2 m/z; number of cycles: 1) and SNAP algorithm for peaks detection (signal-to-noise ratio: 6 for MS and 1.5 for MS/MS). The charge state of the peptides was assumed to be +1. Fragments of porcine trypsin were used for internal mass calibration. Proteins were subsequently identified with PEAKS (v.5.2, Bioinformatics Solutions) using MS/MS-based identification and *de novo* sequencing. The search parameters against a custom-made database were set as follows: carbamidomethylation of cysteine was set as a fixed modification, oxidation of methionine and phosphorylation of serine, threonine and tyrosine were set as variable modifications, one missing cleavage during trypsin digestion was allowed and the tolerance for precursor-ion mass tolerance was set to 1 Da. The database was constructed by combining *P. maximus* sequences from two sources: sequences acquired during this study [24] and sequences from a recent transcriptional analysis of *P. maximus* hemocytes cells [28]. Overall, the database included a total of 252,888 *P. maximus* expressed sequence tags (ESTs). Protein identification was considered as unambiguous when a minimum of two peptides matched with a minimum score of 20. False discovery rates were also estimated using a reverse database as decoy. EST database sequences were annotated by sequence similarity searches against a non-redundant database using the Blast algorithm from NCBI (http://www.ncbi.nlm.nih.gov/BLAST).

## Network analysis

To better understand the relationships between the proteins and genes identified in this study,

network analyses using the String 9.05 algorithm [86] (http://string-db.org) were performed. STRING is a database of known and predicted protein interactions, which includes both direct (physical) and indirect (functional) associations. STRING quantitatively integrates interaction data and calculate a score for each protein interaction based on the genomic context, high-throughput experiments, co-expression and literature. This kind of network analysis is strongly dependent on the quality of the used database, and is thus much more efficient using databases from the most studied model organisms. Hence, the putative human homologues of the genes and proteins up- or down- regulated in this study were identified and used as the input data for the STRING program.

## Availability of supporting data

The sequence data related to this work has been deposited in the GenBank SRA, accession number: SRP040427. The contigs, and annotation are available from http://ramada.nerc-bas.ac.uk/repository/entry/show/Polar+Data+Centre/NERC-BAS+Datasets/Genomics/.

## Additional file

**Additional file 1: Table S1.** List of genes found deregulated at 25 °C. E-values, gene names and accession numbers are those given by Blast2GO for the best hit. Short names are, when possible, those of the best human homolog of the best hit (these names were used for String 9.05 analyses). Induction factors are given as Log2. **Table S2.** List of genes found deregulated at 21 °C. E-values, gene names and accession numbers are those given by Blast2GO for the best hit. Short names are, when possible, those of the best human homolog of the best hit (these names were used for String 9.05 analyses). Induction factors are given as Log2. (DOCX 251 kb)

## Abbreviations

MT: Metallothionein; HSP: Heat shock protein; 2-DE: Two-dimensional protein electrophoresis; Q-PCR: Quantitative PCR; NGS: Next-Generation Sequencing; MALDI TOF-TOF: Matrix-assisted laser desorption ionization time-of-flight tandem mass spectrometry; EST: Expressed sequence tag;

GAPDH: Glyceraldehyde 3 phosphate dehydrogenase; CALR: Calreticulin;

ANK1: Ankyrin; CCAR1: Cell division cycle and apoptosis regulator;

AFMID: Arylformamidase; G3P: Glyceraldehyde-3-phosphate; bZIP: Basic

region leucine zipper; BATF: Basic leucine zipper transcription factor, ATF like;

CREB: cAMP response element-binding; TLR4: Toll-like receptor 4;

GNRH: Gonadotropin-releasing hormone; CA2: Carbonic anhydrase 2;

TRAF: Tumor necrosis factor receptor associated protein;

PAIP1: Polyadenylate binding protein-interacting protein 1;

HNRPD: Heterogeneous nuclear ribonucleoprotein D; CSDE1: Cold shock

domain containing E1; SYNCRIP: Heterogenous nuclear ribonucleoprotein;

JNK: c-JUN N-terminal kinase; CASP: Caspase; UDG: Uracil DNA glycosylase;

CAN: Acetonitrile; SDS: Sodium dodecyl sulfate; DTT: Dithiothreitol;

BICAM: Ammonium bicarbonate; FDR: False discovery rate;

TFA: Trifluoroacetic acid; HCCA: α-cyano-4-hydroxycinnamic acid;

PMF: Peptide mass fingerprinting; MS/MS: Tandem mass spectrometry.

## Competing interests

The authors declare that they have no competing interests.
## Authors' contributions

SA, JR, JFSM, FJ, MSC, VP conceived the study and designed the experiments. SA, JR, RL, JFSM, FJ performed the experiments and analyzed the physiological data. SA, MAST, JR, MSC, VP analyzed the transcriptome data. SA, JR, MSC, VP analyzed the proteome data. JR, LSP and MSC revised the manuscript. SA, VP and MSC wrote the manuscript. All authors approved the final version.

## Acknowledgements

This research was granted by the Région Bretagne, i.e. the Pemadapt project (ref. 6368) and a doctoral fellowship to S.A. (Protmar project, ref. 6197). MSC, MAST and LSP were financed by core NERC funding to the British Antarctic Survey Polar Sciences for Planet Earth programme. COMANCHE (ANR-2010- STRA-010) and LabexMER (ANR-10-LABX-19) programs from the French National Research Agency (ANR), and a NERC Theme Leader Award to LSP (NE/G018), also supported our research. For temperatures in the bay of Brest, we used data provided by the "Service d'Observation en Milieu Littoral", INSU-CNRS, SOMLIT-Brest, IUEM/UBO.

## Author details

¹Laboratoire des Sciences de l'Environnement Marin, LEMAR UMR 6539 CNRS/UBO/IRD/Ifremer, Université de Brest (UBO), Institut Universitaire Européen de la Mer, Plouzané 29280, France. ²British Antarctic Survey, Natural Environment Research Council, High Cross, Madingley Road, Cambridge CB3 0ET, UK.

Received: 20 June 2015 Accepted: 22 October 2015

Published online: 23 November 2015

## References

1. Harley CDG, Randall Hughes A, Hultgren KM, Miner BG, Sorte CJB, Thomber CS, et al. The impacts of climate change in coastal marine systems. Ecol Lett. 2006;9:228–41.

2. Perry AL, Low PJ, Ellis JR, Reynolds JD. Climate change and distribution shifts in marine fishes. Science. 2005;308:1912–5.

3. Rosenzweig C, Karoly D, Vicarelli M, Neofotis P, Wu Q, Casassa G, et al. Attributing physical and biological impacts to anthropogenic climate change. Nature. 2008;453:353–7.

4. Clark MS, Thorne MAS, Amaral A, Vieira F, Batista FM, Reis J, et al. Identification of molecular and physiological responses to chronic environmental challenge in an invasive species: the Pacific oyster, *Crassostrea gigas*. Ecol Evol. 2013;3:3283–97.

5. Matoo OB, Ivanina AV, Ullstad C, Beniash E, Sokolova IM. Interactive effects of elevated temperature and CO(2) levels on metabolism and oxidative stress in two common marine bivalves (*Crassostrea virginica* and *Mercenaria mercenaria*). Comp Biochem Physiol A Mol Integr Physiol. 2013;164:545–53.

6. Gazeau F, Alliouane S, Bock C, Bramanti L, López Correa M, Gentile M, et al. Impact of ocean acidification and warming on the Mediterranean mussel (*Mytilus galloprovincialis*). Front. Mar. Sci. 2014;1:62. doi:10.3389/fmar.2014.00062

7. Murtha J. Characterization of the heat shock response in mature zebrafish (*Danio rerio*). Exp Gerontol. 2003;38:683–91.

8. Lenski RE, Bennett AF. Evolutionary Response of *Escherichia coli* to Thermal Stress. Am Nat. 1993;142:547.

9. Liu XD, Liu PC, Santoro N, Thiele DJ. Conservation of a stress response: human heat shock transcription factors functionally substitute for yeast HSF. EMBO J. 1997;16:6466–77.

10. Konkel M, Kim B, Klena J, Young C, Ziprin R. Characterization of the thermal stress response of *Campylobacter jejuni*. Infect Immun. 1998;66:3666–72.

11. Wang W, Vinocur B, Shoseyov O, Altman A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004;9:244–52.

12. Ivanina AV, Taylor AC, Sokolova IM. Effects of elevated temperature and cadmium exposure on stress protein response in eastern oysters *Crassostrea virginica* (Gmelin). Aquat Toxicol. 2009;91:245–54.

13. Brun NT, Bricelj VM, MacRae TH, Ross NW. Heat shock protein responses in thermally stressed bay scallops, *Argopecten irradians*, and sea scallops, *Placopecten magellanicus*. J Exp Mar Bio Ecol. 2008;358:151–62.

14. Brun NT, Bricelj VM, MacRae TH, Ross NW. Acquisition of thermotolerance in bay scallops, *Argopecten irradians* irradians, via differential induction of heat shock proteins. J Exp Mar Bio Ecol. 2009;371:77–83.

15. Artigaud S, Lacroix C, Pichereau V, Flye-Sainte-Marie J. Respiratory response to combined heat and hypoxia in the marine bivalves *Pecten maximus* and *Mytilus* spp. Comp Biochem Physiol A Mol Integr Physiol. 2014;175:135–40.

16. Artigaud S, Lacroix C, Richard J, Flye-Sainte-Marie J, Bargelloni L, Pichereau V. Proteomic responses to hypoxia at different temperatures in the great scallop (*Pecten maximus*). PeerJ. 2015;3:e871.

17. Lang RP, Bayne C, Camara M, Cunningham C, Jenny M, Langdon C. Transcriptome Profiling of Selectively Bred Pacific Oyster *Crassostrea gigas* Families that Differ in Tolerance of Heat Shock. Mar Biotechnol. 2009;11:650–68.

18. Lockwood BL, Sanders JG, Somero GN. Transcriptomic responses to heat stress in invasive and native blue mussels (genus *Mytilus*): molecular correlates of invasive success. J Exp Biol. 2010;213:3548–58.

19. 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.

20. Clark MS, Peck LS. HSP70 heat shock proteins and environmental stress in Antarctic marine organisms: A mini-review. Mar Genomics. 2009;2:11–8.

21. Fu X, Sun Y, Wang J, Xing Q, Zou J, Li R, et al. Sequencing-based gene network analysis provides a core set of gene resource for understanding thermal adaptation in Zhikong scallop *Chlamys farreri*. Mol Ecol Resour. 2014;14:184–98.

22. 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–402.

23. Brand A: The European scallop fisheries for *Pecten maximus*, *Aequipecten opercularis* and *Mimachlamys varia*. In Developments in aquaculture and fisheries science. Volume 35. Elsevier Science, Boston. Edited by Shumway SE, Parsons GJ; 2006:991–1058.

24. Artigaud S, Thorne MAS, Richard J, Lavaud R, Jean F, Flye-Sainte-Marie J, et al. Deep sequencing of the mantle transcriptome of the great scallop *Pecten maximus*. Mar Genomics. 2014;15:3–4.

25. Jackson DJ, McDougall C, Woodcroft B, Moase P, Rose RA, Kube M, et al. Parallel Evolution of Nacre Building Gene Sets in Molluscs. Mol Biol Evol. 2010;27:591–608.

26. Forné I, Abián J, Cerdà J. Fish proteome analysis: model organisms and non-sequenced species. Proteomics. 2010;10:858–72.

27. Kinoshita S, Wang N, Inoue H, Maeyama K, Okamoto K, Nagai K, et al. Deep sequencing of ESTs from nacreous and prismatic layer producing tissues and a screen for novel shell formation-related genes in the pearl oyster. PLoS One. 2011;6:e21238.

28. Pauletto M, Milan M, Moreira R, Novoa B, Figueras A, Babbucci M, et al. Deep transcriptome sequencing of *Pecten maximus* hemocytes: a genomic resource for bivalve immunology. Fish Shellfish Immunol. 2014;37:154–65.

29. Dowd WW. Challenges for biological interpretation of environmental proteomics data in non-model organisms. Integr Comp Biol. 2012;52:705–20.

30. Dickie L. Effects of high temperature on survival of the giant scallop. J Fish Res board Canada. 1958;15:1189–211.

31. Laing I. Effect of temperature and ration on growth and condition of king scallop (*Pecten maximus*) spat. Aquaculture. 2000;183:325–34.

32. Schmidt-Nielsen K. Animal Physiology: Adaptation and Environment. 4th ed. Cambridge, New York, Melbourne: Cambridge University Press; 1990.

33. Nicholls C, Li H, Liu J-P. GAPDH: A common enzyme with uncommon functions. Clin Exp Pharmacol Physiol. 2012;39:674–9.

34. Mitsuzawa H, Kimura M, Kanda E, Ishihama A. Glyceraldehyde-3-phosphate dehydrogenase and actin associate with RNA polymerase II and interact with its Rpb7 subunit. FEBS Lett. 2005;579:48–52.

35. Lee MN, Ha SH, Kim J, Koh A, Lee CS, Kim JH, et al. Glycolytic flux signals to mTOR through glyceraldehyde-3-phosphate dehydrogenase-mediated regulation of Rheb. Mol Cell Biol. 2009;29:3991–4001.

36. Sirover MA. New nuclear functions of the glycolytic protein, glyceraldehyde-3-phosphate dehydrogenase, in mammalian cells. J Cell Biochem. 2005;95:45–52.

37. Tisdale EJ. Glyceraldehyde-3-phosphate dehydrogenase is phosphorylated by protein kinase C/λ and plays a role in microtubule dynamics in the early secretory pathway. J Biol Chem. 2002;277:3334–41.

38. Leppä S, Bohmann D. Diverse functions of JNK signaling and c-Jun in stress response and apoptosis. Oncogene. 1999;18:6158–62.

39. Ayemar M, Wisniewska M, Weitzman JB. A role for AP-1 in apoptosis: the case for and against. Biochimie. 2003;85:747–52.
40. Hess J, Angel P, Schorpp-Kistner M. AP-1 subunits: quarrel and harmony among siblings. J Cell Sci. 2004;117:5965-73.

41. Grosset C, Chen C-YA XN, Sonenberg N, Jacquemin-Sablon H, Shyu A-B. A Mechanism for Translationally Coupled mRNA Turnover: Interaction between the Poly(A) Tail and a c-fos RNA Coding Determinant via a Protein Complex. Cell. 2000;103:29-40.

42. Hu HM, O'Rourke K, Boguski MS, Dixit VM. A novel RING finger protein interacts with the cytoplasmic domain of CD40. J Biol Chem. 1994;269:30069-72.

43. Hildebrand JM, Yi Z, Buchta CM, Poovassery J, Stunz LL, Bishop GA. Roles of tumor necrosis factor receptor associated factor 3 (TRAF3) and TRAF5 in immune cell functions. Immunol Rev. 2011;244:55-74.

44. Salminen A, Huuskonen J, Ojala J, Kauppinen A, Kaarniranta K, Suuronen T. Activation of innate immunity system during aging: NF-kB signaling is the molecular culprit of inflammm-aging. Ageing Res Rev. 2008;7:83-105.

45. Gross M. Emergency Services: A Birds Eye Perspective on the Many Different Functions of Stress Proteins. Curr Protein Pept Sci. 2004;5:213-23.

46. Sørensen JG, Loeschcke V. Studying stress responses in the post-genomic era: its ecological and evolutionary role. J Biosci. 2007;32:447-56.

47. Clark MS, Peck LS. Triggers of the HSP70 stress response: environmental responses and laboratory manipulation in an Antarctic marine invertebrate (*Nacella concinna*). Cell Stress Chaperones. 2009;14:649-60.

48. Roberts RJ, Agius C, Saliba C, Bossier P, Sung YY. Heat shock proteins (chaperones) in fish and shellfish and their potential role in relation to fish health: a review. J Fish Dis. 2010;33:789-801.

49. Landry J, Chrétien P, Lambert H, Hickey E, Weber LA. Heat shock resistance conferred by expression of the human HSP27 gene in rodent cells. J Cell Biol. 1989;109(1):7-15.

50. Musch MW, Sugi K, Straus D, Chang EB. Heat-shock protein 72 protects against oxidant-induced injury of barrier function of human colonic epithelial Caco2/bbe cells. Gastroenterology. 1999;117:115-22.

51. Zhang G, Fang X, Guo X, Li L, Luo R, Xu F, et al. The oyster genome reveals stress adaptation and complexity of shell formation. Nature. 2012;490:49-54.

52. Luchmann KH, Clark MS, Bainy ACD, Gilbert JA, Craft JA, Chipman JK, et al. Key metabolic pathways involved in xenobiotic biotransformation and stress responses revealed by transcriptomics of the mangrove oyster *Crassostrea brasiliana*. Aquat Toxicol. 2015;166:10-20.

53. Henkels KM, Turchi JJ. Cisplatin-induced apoptosis proceeds by caspase-3-dependent and-independent pathways in cisplatin-resistant and-sensitive human ovarian cancer cell lines. Cancer Res. 1999;59:3077-83.

54. Tomanek L, Zuzow MJ, Ivanina AV, Beniash E, Sokolova IM. Proteomic response to elevated PCO₂ level in eastern oysters, *Crassostrea virginica*: evidence for oxidative stress. J Exp Biol. 2011;21(Pt 11):1836-44.

55. Clarke A, Johnston NM. Scaling of metabolic rate with body mass and temperature in teleost fish. J Anim Ecol. 1999;68:893-905.

56. Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, et al. Genomic expression programs in the response of yeast cells to environmental changes. Sci Signal. 2000;11:4241.

57. Buckley BA, Gracey AY, Somero GN. The cellular response to heat stress in the goby *Gillichthys mirabilis*: a cDNA microarray and protein-level analysis. J Exp Biol. 2006;209:2660-77.

58. Gracey AY, Chaney ML, Boomhower JP, Tyburczy WR, Connor K, Somero GN. Rhythms of Gene Expression in a Fluctuating Intertidal Environment. Curr Biol. 2008;18:1501-7.

59. Cossins AR. Temperature Adaptation of Biological Membranes. London: Chapel Hill, NC: Portland Press; 1994. c1994.

60. Hazel JR. Thermal adaptation in biological membranes: is homeoviscous adaptation the explanation? Annu Rev Physiol. 1995;57:19-42.

61. Farah ME, Sirotkin V, Haarer B, Kakhnashvili D, Amberg DC. Diverse protective roles of the actin cytoskeleton during oxidative stress. Cytoskeleton (Hoboken). 2011;68:340-54.

62. Tomanek L. Environmental proteomics of the mussel *Mytilus*: implications for tolerance to stress and change in limits of biogeographic ranges in response to climate change. Integr Comp Biol. 2012;52:648-64.

63. Galhardo RS, Hastings PJ, Rosenberg SM. Mutation as a stress response and the regulation of evolvability. Crit Rev Biochem Mol Biol. 2007;42:399-435.

64. Goncharova EI, Nádas A, Rossman TG. Serum deprivation, but not inhibition of growth per se, induces a hypermutable state in Chinese hamster G12 cells. Cancer Res. 1996;56:752-6.

65. Koshiji M, To KK-W, Hammer S, Kumamoto K, Harris AL, Modrich P, et al. HIF-1α induces genetic instability by transcriptionally downregulating MutSa expression. Mol Cell. 2005;17:793-803.

66. Mihaylova VT, Bindra RS, Yuan J, Campisi D, Narayanan L, Jensen R, et al. Decreased expression of the DNA mismatch repair gene Mlh1 under hypoxic stress in mammalian cells. Mol Cell Biol. 2003;23:3265-73.

67. Yuan J, Narayanan L, Rockwell S, Glazer PM. Diminished DNA repair and elevated mutagenesis in mammalian cells exposed to hypoxia and low pH. Cancer Res. 2000;60:4372-6.

68. Shor E, Fox CA, Broach JR. The Yeast Environmental Stress Response Regulates Mutagenesis Induced by Proteotoxic Stress. PLoS Genet. 2013;9:e1003680.

69. Somero GN. The physiology of climate change: how potentials for acclimatization and genetic adaptation will determine “winners” and “losers”. J Exp Biol. 2010;213:912-20.

70. Truebano M, Burns G, Thorne MAS, Hillyard G, Peck LS, Skibinski DOF, et al. Transcriptional response to heat stress in the Antarctic bivalve *Laternula elliptica*. J Exp Mar Bio Ecol. 2010;391:65-72.

71. Utting S, Millican P. The role of diet in hatchery conditioning of *Pecten maximus* L.: a review. Aquaculture. 1998;165:167-78.

72. Lucas A, Beninger PG. The use of physiological condition indices in marine bivalve aquaculture. Aquaculture. 1985;44:187-200.

73. Clark MS, Thorne MA, Vieira FA, Cardoso JC, Power DM, Peck LS. Insights into shell deposition in the Antarctic bivalve *Laternula elliptica*: gene discovery in the mantle transcriptome using 454 pyrosequencing. BMC Genomics. 2010;11:362.

74. Joubert C, Piquemal D, Marie B, Manchon L, Pierrat F, Zanella-Cléon I, et al. Transcriptome and proteome analysis of *Pinctada margaritifera* calcifying mantle and shell: focus on biomineralization. BMC Genomics. 2010;11:613.

75. Schmieder R, Lim YW, Edwards R. Identification and removal of ribosomal RNA sequences from metatranscriptomes. Bioinformatics. 2012;28:433-5.

76. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10-2.

77. Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigsci. 2012;1:18.

78. Huang X. CAP3: A DNA Sequence Assembly Program. Genome Res. 1999;9:868-77.

79. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658-9.

80. Link AJ, Ramagli LS. Quantifying Protein in 2-D PAGE Solubilization Buffers BT - 2-D Proteome Analysis Protocols. In: Totowa NJ, editor. 2-D Proteome Analysis Protocols, vol. 112. 1998. p. 99-104. Methods in Molecular Biology.

81. R Core Team R version 3. 0.2: R: A language and environment for statistical computing. *R Found Stat Comput Vienna, Austria* 2013: http://www.R-project.org

82. Artigaud S, Gauthier O, Pichereau V. Identifying differentially expressed proteins in two-dimensional electrophoresis experiments: inputs from transcriptomics statistical tools. Bioinformatics. 2013;29:2729-34.

83. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.

84. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289-300.

85. Artigaud S, Lavaud R, Thébault J, Jean F, Strand Ø, Strohmeier T, et al. Proteomic-based comparison between populations of the Great Scallop, *Pecten maximus*. J Proteomics. 2014;105:164-73.

86. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41 (Database issue):D808-815.