Both the upwelling site HBH (21˚56'18.01"N, 120˚44'45.54"E) and the non-upwelling site HWN (22˚01'23.30"N, 120˚41'18.29"E) have been well characterized with respect to their oceanography [13] and coral reef ecology (e.g., [22] [23] ). Although the difference in temperature variability between study sites was the sole focus of prior environmental physiology works on seriatoporid corals from HBH and HWN [12] , other oceanographic variables were shown by Liu et al. [22] to differ between them, and these additional seawater quality parameters could likewise influence the physiology of the resident corals. We therefore undertook a meta-analysis of seawater quality at each site by pooling our own published data ( [12] [14] [15] ) with those of Liu et al. [22] . We now briefly review how these data were collected.

The temperature regimes of both the upwelling site HBH and the non-upwelling, Taiwan Strait, control site HWN (Figure 1) were characterized in detail in 2010 at 7 - 8 m depth (the depth of coral collection [12] ); specifically, temperature was measured at hourly intervals with HOBO® Pendant data loggers (Onset) for one year. Although the mean annual temperature did not differ between the two sites, the variability did (Table 1); the mean monthly temperature range was 2-fold higher at HBH than HWN due to spring-tide upwelling events that occur during the boreal summer at the former site only [12] .

As mentioned above, prior works on seriatoporid corals from HBH and HWN have focused mainly on these temperature differences as being responsible for the physiological heterogeneity documented across the sampled coral colonies (e.g., [12] [14] [15] [16] [17] ). However, we considered additional seawater quality parameters herein as being potential drivers of physiological variation between corals of the two study sites by incorporating the data from Liu et al. [22] , in which stations 2 and 7 - 8 correspond to HWN and HBH, respectively. Specifically, Liu et al. [22] measured salinity, pH, dissolved oxygen (DO) content (%), biochemical oxygen demand (BOD₅; mg/L), suspended solid levels (mg/L), turbidity (NTU), and concentrations of nitrite (μg/L), nitrate (μg/L), ammonia (μg/L), phosphate (μg/L), silicate (μg/L), and chlorophyll a (chl-a; μg/L) at monthly intervals between 2001 and 2008, and their respective averages across this timespan have been included in Table 1. Photosynthetically active radiation (PAR) was instead measured at the depth of coral collection (7 - 8 m) at hourly intervals by Mayfield et al. [12] , with the mean hourly daytime PAR calculated across several days in May 2010 (the month of coral collection) presented in Table 1. Finally, coral cover was estimated at each site as in Tkachenco et al. [24] , with modifications made by Liu et al. [22] , and the mean average live coral cover (ALCC) values of the latter work have been included in Table 1 for both study sites. Data were compared statistically across sites as described below.

Sampling of six S. hystrix colonies (see insets of Figure 1) was undertaken at each site in May 2010 (a month in which upwelling events were frequent [12] ) during a period of stable (i.e., non-upwelling) temperatures (~26˚C) as in Mayfield et al. [14] . Only visibly healthy corals were sampled, and there were no evident signs of stress (e.g., excessive release of mucus, bacterial infection, bleaching, recent tissue loss, etc.) in any of the 12 colonies. Once at the surface, three small (~50 mg) biopsies taken from each colony were submerged in ~2 ml of TRIzol™ (Life Technologies), and the 36 biopsies were later homogenized with a mortar and pestle at the lab as in Mayfield et al. [14] . Additional tissue biopsies were taken immediately upon removal of the colonies from the ocean in order to determined in situ chl-a content and Symbiodinium density (described in detail in Mayfield et al. [14] ). The colonies from which these “in situ biopsies” were taken were transported to the laboratory and fragmented into nubbins for use in a laboratory-based reciprocal transplant [12] . Physiological- and molecular-scale

^aunequal variance; ^blog-transformed data; ^csquare root-transformed data.

data from the same coral colonies as those analyzed herein ( [12] [14] [15] [16] [17] ) have been summarized in Table 1 in a similar manner as for the previously acquired seawater quality data in order to more thoroughly assess environmental differences in the physiology of S. hystrix in Southern Taiwan. This meta analysis was also carried out to better understand the effect of experimental husbandry on coral physiology (sensu Mayfield et al. [14] ), since Symbiodinium densities and chl-a concentrations, as well as the maximum dark-adapted yield of photosystem II (Fv/Fm), were documented both in situ and after four weeks of aquarium husbandry for corals of both study sites.

Proteins from one of the three technical replicates from each of three randomly chosen colonies from each of the two sites were purified as in Mayfield et al. [16] . These six samples represent those biopsies preserved immediately upon removal from the source colonies from the ocean and were therefore meant to be representative of “in situ protein concentrations.” Proteins (n = 6) were extracted as in Mayfield et al. [25] , precipitated in acetone, washed, dried, dissolved in ~150 μL of rehydration buffer (9.5 M urea, 2% CHAPS, 0.5% carrier ampholytes, and 65 mM dithiothreitol), and 20 μL of solubilized protein were quantified with the 2D Quant kit (Amersham Biosciences). The transcriptomes of all 12 colonies, including the 6 whose proteomes were analyzed herein, were sequenced previously [17] , albeit after three weeks of acclimation and one week of experimentation (stable vs. variable temperature regimes [12] ).

Proteins (~130 μg) were electrophoresed across 2D as in Mayfield et al. [16] ; however, as this book chapter is not freely available except on ABM’s personal website (http://coralreefdiagnostics.com/), certain key details, which can also be found in Mayfield et al. [26] , been reiterated below. Briefly, upon electrophoresing the protein samples (n = 3 for each of two sites of origin: HBH vs. HWN) across 2D (isoelectric focusing for isoelectric point [pI] determination and SDS-PAGE for molecular weight [in kDa] determination) as described in Mayfield et al. [26] , gels were fixed and stained with SYPRO® Ruby (Life Technologies). They were then imaged by a Typhoon Trio™ scanner (GE Healthcare), and ImageQuant TL software (GE Healthcare) was used to identify USPs, or, when USPs were not uncovered, DCPs, between each pair of HBH and HWN gels. Two gels were run simultaneously on each of three days: one HWN sample and one HBH sample. DCP/USP pIs and molecular weights were recorded from the 1^st and 2^nd HBH vs. HWN gel pairs. When these same protein spots were also differentially concentrated between gels, or uniquely synthesized by corals of one of the two sites, in the third pair of biological replicates, they were removed from the gel in which their concentration was higher (referred to herein as the “representative” gel), digested with trypsin, and prepared for MS as in Mayfield et al. [16] [26] .

In general, only protein spots that were uniquely synthesized by coral samples of one site and not the other in all three pairs of gels were removed from the final, representative gel and sequenced, though two protein spots that were instead more concentrated by HBH samples in all three gel pairs (i.e., DCPs, rather than USPs) were removed from the final HBH gel and analyzed by MS. In total three HBH > HWN and six HWN > HBH spots were in-gel digested with trypsin and purified as in Mayfield et al. [16] [26] . Then, 2 μL of purified, trypsin-digested peptides were injected into a nano-liquid chromatography system and detected by an LTQ Orbitrap “Discovery Hybrid Fourier Transform” mass spectrometer (Thermo-Fisher) at a resolution of 30,000 coupled with a nanospray source that was executed in positive ion mode. The nano-UPLC system (“nanoACQUITY”), desalting (Symmetry C18; 5 μm × 180 μm × 20 mm), and analytical (BEH C18; 1.7 μm × 75 μm × 150 mm) columns were all purchased from Waters. The peptide eluate from the column was directed to the nanospray source, and the MS was operated in data-dependent mode.

When oceanographic (e.g., temperature), ecological (ALCC), and molecular physiological (i.e., the 23 response variables analyzed previously in the sampled colonies; Table 1) data were normally distributed and of homogeneous variance across the two study sites, their means were compared across sites with student’s t-tests. Wilcoxon tests were used when log- or square root-transformations did not lead to normally distributed (Shapiro-Wilk test, P < 0.05) or homogeneously variable (Levene’s test, P < 0.05) datasets between sites. In certain cases, both in situ data (Symbiodinium density, Fv/Fm, and chl-a content) and data following four weeks of experimental husbandry (described in Mayfield et al. [14] ) were available. In these cases, repeated-measures ANOVAs were instead used to determine whether husbandry had a differential effect on corals of the two study sites. All univariate statistical analyses were carried out with JMP® (ver. 12.0.1). Finally, a multivariate approach aimed at modeling site differences in coral physiology (sensu Mayfield [27] ) was carried out with PRIMER (ver. 6) with all 23 response variables (standardized prior to analysis), testing site of origin as the fixed factor. PERMANOVA is based on similarity (Bray-Curtis) between samples and so does not require that data are normally distributed and of homogenous variance (as does MANOVA, which is sensitive to significant deviations from normality and cannot be used in cases such as this study in which there are more response variables [n = 23] than samples [n = 12]).

.MGF data files from the mass spectrometer (n = 6) were directly uploaded into the MS-SCAN program featured on the S. hystrix-Symbiodinium transcriptome server (http://symbiont.iis.sinica.edu.tw/s_hystrix/static/html/#mscan), and all default conditions of the MS-GF+ script [28] upon which MS-SCAN is based were used (discussed on the open access website housing the script: https://omics.pnl.gov/software/ms-gf); this included up to two missed cleavages allowed. However, as 1) Chiva et al. [29] found that sequence datasets featuring such missed cleavages are not inherently biased with respect to quantification, and 2) Mayfield et al. [16] found the S. hystrix-Symbiodinium proteome to be lysine- and arginine-rich, we generally included sequences containing over two missed cleavages provided that either 15 amino acids (AA) were sequenced at minimum or, alternatively, two peptides mapped to the same reference protein whose collective length was ≥15 AA. A decoy database (sensu [30] ) was not queried to calculate a false discovery rate since there is no fully sequenced proteome for a coral or Symbiodinium. Instead, the aforementioned peptide length criteria were increased above the default minimum of 6 AA of MS-GF+ to ensure that only peptides that could be identified with confidence were considered in the analysis. This undoubtedly reduced the total number of DCPs identified. Additional details of MS-GF+, which, unlike Mascot (Matrix Sciences), is free, can be found in Kim and Pevzner [28] .

When peptides could be identified with confidence (≥15 AA mapping to a translated contig in the S. hystrix-Symbiodinium transcriptome) with MS-SCAN, they were assigned a compartment of origin (host, Symbiodinium, intermediate [either host coral or Symbiodinium], or unknown), as well as a protein identity and functional category (e.g., metabolism) when the top hit contig (mRNA) hypothetically encoding the sequenced peptide aligned significantly (e < 10⁻⁵) to a functionally characterized protein in the NCBI nr database. The full mRNA sequence of the top hit contig derived from MS-SCAN analysis of the spectral data, rather than the peptide sequence itself, was used as the query (BLASTx) of the NCBI database; this is because trypsin-digested peptides are generally short (mean length = 22 ± 9 [std. dev.] AA herein; see Table S1 located at the end of the manuscript.), and such short sequences may not align significantly to any homolog on a public sequence repository like NCBI. However, when >30 contiguous AA were sequenced from a single protein, the peptide sequence itself was additionally BLASTed (BLASTp) against the NCBI database to attempt to corroborate the BLASTx analysis of the respective mRNA.

The compartmental breakdowns of the differentially concentrated proteomes were compared to the S. hystrix: Symbiodinium mRNA ratio of 1.8 [17] with z-tests. This approach aimed to determine whether one member of the endosymbiosis contributed relatively more USPs/DCPs than the other for 1) all DCPs, 2) the HWN > HBH DCPs, and 3) the HBH > HWN DCPs. Two-sample proportion tests were used to determine whether 1) certain functional categories differed in proportional abundance between the HWN > HBH and HBH > HWN proteomes and 2) certain functional categories were significantly over-represented relative to the stable vs. variable differentially concentrated proteome of Mayfield et al. [17] . For all such statistical analyses, an alpha level of 0.05 was established a priori.

As the transcriptomes of nubbins generated from the colonies analyzed herein were sequenced previously (6 nubbins from each site of origin sequenced after one month of experimental husbandry [17] ), the expression data from all mRNAs encoding the 53 DCPs uncovered herein were acquired from the S. hystrix-Symbiodinium transcriptome server (http://symbiont.iis.sinica.edu.tw/s_hystrix/static/html/#stat), and 2-way ANOVAs were performed to determine the effects of site of origin (HWN vs. HBH), temperature treatment (stable vs. variable; see Mayfield et al. [12] [17] for details.), and their interaction. When a site of origin difference was statistically significant (P < 0.05) and matched that observed at the protein level documented herein, congruency between mRNA expression and protein concentration was said to have occurred. For instance, if the 2-way ANOVA revealed that a gene was expressed at 2-fold-higher levels in corals of HBH, and the protein spot was more highly concentrated (or uniquely identified) in the final, representative HBH protein gel, molecular congruency was deemed to have been verified. Two-sample proportion tests were used to determine whether congruency differed significantly across 1) compartments (host coral vs. Symbiodinium), 2) site of origin (HBH vs. HWN), and 3) experiments (this study compared to a variable temperature study performed with these same samples [17] ).

Although most seawater quality parameters were similar between HWN and HBH, several differed significantly upon pooling data over an eight-year period (Table 1); notably, silicate concentration, suspended solid levels, and turbidity were all significantly higher at HBH (P < 0.05), which abuts a marina, than HWN, which is adjacent to Taiwan’s National Museum of Marine Biology and Aquarium (NMMBA; Figure 1; where all laboratory work was undertaken). As reported previously [12] , the mean monthly temperature range at HBH of ~6˚C was approximately double that of HWN (~3˚C; Table 1) in the year of coral sampling (2010).

In addition to the 2D + MS results generated herein, we aimed to first provide a brief overview of previous findings obtained from these same coral colonies. We first summarize the influence of a one-month experimental husbandry on coral nubbins from the two study sites that were fragmented from the same source colonies from which biopsies were taken herein for 2D + MS analysis (“Effect of experimental husbandry”); these findings are described in detail in Mayfield et al. [14] . Then, we briefly summarize the findings of two prior works (Mayfield et al. [12] [15] ) in which 23 molecular-physiological response variables were measured in coral nubbins generated from the same source colonies as those analyzed herein, albeit exposed to stable or variable temperatures for one week (following three weeks of acclimation; “Coral molecular physiology differences between upwelling (HBH) and non-upwelling (HWN) study sites”). Upon presenting the 2D + MS data produced herein (“2D + MS”), we then compare the results obtained to those of another study with these same samples [17] that instead sought to model temperature-related, rather than site of origin-associated, differences in protein concentrations (“Comparison with a stable vs. variable temperature regime study carried out with S. hystrix”). It should be mentioned that all 12 colonies were of the identical genotype (determined by analysis of microsatellites [15] ) and hosted Symbiodinium of clade C exclusively [15] .

As discussed in more detail in Mayfield et al. [14] , some coral response variables did change in response to four weeks of experimental husbandry (three weeks of 26˚C-acclimation and one week of experimentation [stable vs. variable temperature regime]; see Mayfield et al. [12] for details.), and husbandry differentially affected corals of the two sites; Symbiodinium density increased more after four weeks of husbandry for corals of HWN than for those of HBH, whereas chl-a concentration and Fv/Fm increased more over this period for corals of HBH (Table 1).

Of the 23 response variables measured in corals of the two study sites (Table 1; see Mayfield et al. [12] [15] for details.), only the protein/DNA ratio and Symbiodinium photosystem I (subunit III; psI) mRNA expression differed significantly between corals of the upwelling and non-upwelling sites (Table 1). Regarding the former, S. hystrix colonies of HWN were characterized by 1.4-fold higher protein/DNA ratios than colonies of the same genotype and Symbiodinium assemblages of HBH. It should be mentioned that, despite this difference, equal quantities of protein were loaded into all 2D gels. In contrast, psI mRNA expression was 1.5-fold higher in Symbiodinium (clade C [15] ) populations within corals of HBH. Although PERMANOVA did not detect a multivariate effect of site of origin across the 23 response variables (Table 1), it is possible that a hypothesis-neutral, proteomics-based approach could nevertheless uncover site-related differences in coral proteo-biology; such was indeed the case, and these findings are discussed in detail below.

Of the 6 and 3 protein spots concentrated at higher levels by samples of HWN (Figure 2(a)) and HBH (Figure 2(b)), respectively, 6 and 1, respectively, were uniquely synthesized by corals of one site of origin and not the other; two DCPs were additionally isolated from the representative HBH gel (Figure 2(b)). From the 6 HWN > HBH and 3 HBH > HWN protein spots, 38 (Table 2) and 15 (Table 3) proteins, respectively, were identified with MS-SCAN using the S. hystrix-Symbiodinium transcriptome as a reference database (Figure 3). Regarding the compartmental breakdown of all 53 proteins identified (Figure 3(a)), 55 and 38% were of host coral and Symbiodinium origin, respectively. This 1.4:1 ratio did not differ significantly from the 1.8:1 mRNA ratio of this coral [17] (z-test, P > 0.05), nor did it differ from the 1.6:1 host/Symbiodinium DCP ratio of another proteomic study of S. hystrix [17] (2-sample proportion test, P > 0.05).

With respect to the functional breakdown of all 49 proteins that could be assigned a compartment of origin (Figure 3(d)), nearly half aligned to proteins that had not been characterized. Of the remaining 28 proteins (57%) that did

align to characterized proteins, the dominant functional categories were cytoskeleton, stress response, transport, and transcription. Of the 14 HBH > HWN DCPs (Figure 3(e)), cytoskeleton and the stress response were the most represented categories, and proteins involved in the stress response were significantly more likely to be documented at higher levels by corals of HBH (22% of the differentially expressed proteome) than conspecifics of HWN (Figure 3(f); 3%; 2-sample proportion test, P < 0.05). For the 35 HWN > HBH DCPs (Figure 3(f)), proteins involved in transcription and transport were instead of higher proportional abundance.

When looking only at the 29 host coral DCPs (Figure 3(g)), proteins involved in the stress response, cytoskeleton, and transcription were most abundant. As when the host and Symbiodinium data were pooled (Figure 3(d)-(f); discussed above), the former process was proportionally more abundant in the HBH > HWN proteome (Figure 3(h)); 22% of the 9 HBH > HWN and 0% of the 20 HWN > HWN host DCPs were involved in the stress response (Figure 3(i); 2-sample proportion test, P < 0.05). Instead, transport and transcription were

the most numerically dominant functional categories in the HWN > HBH differentially concentrated host coral proteome (Figure 3(i)).

^aMaintained at lower concentrations in corals exposed to variable temperature regimes [17] ; ^bClosely related to a protein documented at lower concentrations in corals exposed to variable temperature regimes [17] .

^aMaintained at lower concentrations in corals exposed to variable temperature regimes [17] .

In general, the same functional categories differed in proportional abundance between sites of origin for Symbiodinium (Figure 3(j)): stress response, transcription, and transport. However, the Symbiodinium differentially concentrated proteome also featured proteins involved in metabolism. When compared to another study performed with Taiwanese S. hystrix samples [16] [17] , the Symbiodinium proteome sequenced herein featured significantly higher proportions of proteins involved in the stress response, transport, and transcription (2-sample proportion tests, P < 0.05). In the stable vs. variable temperature experiment, on the other hand, proteins associated with lipid bodies and mRNA processing were more likely to be maintained at different levels by Symbiodinium populations (2-sample proportion tests, P < 0.05).

Only five Symbiodinium proteins were identified across the three HBH > HWN protein spots, and only one could be identified with confidence (Figure 3(k)): a stress-associated HSPB1-associated protein 1-like protein was concentrated at higher levels in Symbiodinium from corals of HBH (Table 3). In contrast, the majority of the 15 HWN > HBH Symbiodinium DCPs could be identified with confidence (Figure 3(l)), and these proteins were involved in metabolism, transport, and transcription. No functional category differed in proportional abundance between host and Symbiodinium for the HBH > HWN or HWN > HBH differentially expressed proteomes.

The congruency between mRNA and protein expression did not differ between the two compartments of the S. hystrix-Symbiodinium endosymbiosis (2-sample proportion tests, P > 0.05). Specifically, 2 of the 29 host coral molecules (7%) and 1 of the 20 Symbiodinium molecules (5%) demonstrated congruency between mRNA expression and protein concentration. These three molecules included 1) a host O-aminophenol oxidase involved in the oxidative stress response (the lone HBH > HWN molecule demonstrating congruency between mRNA and protein expression; 7%; Table 3), 2) one host coral molecule of unknown function (HWN > HBH; Table 2), and 3) one Symbiodinium molecule of unknown function (HWN > HBH; Table 2). In addition, one HWN > HBH USP of unknown cellular origin (Table 2) also demonstrated congruency between mRNA expression and protein concentration. The overall mRNA vs. protein congruency of 7.5% (4/53) was significantly higher than that of another study performed with Taiwanese S. hystrix samples (2-sample proportion test, P < 0.05); Mayfield et al. [17] found only 2 molecules out of 167 total (<2%) that demonstrated congruency between mRNA expression and protein concentration. A more detailed comparison with this only other proteomic assessment of S. hystrix [17] can be found below.

Besides a host coral beta-gamma crystallin protein, only three additional proteins were found to be differentially concentrated between sites of origin herein and between stable and variable temperature treatments in S. hystrix nubbins made from these same colonies by Mayfield et al. [17] . All three proteins were of host coral origin, though only two could be identified with confidence: a serine-arginine repetitive matrix protein (RMP) involved in mRNA processing and a translation factor known as eukaryotic translation initiation factor 3 subunit A-like. The serine-arginine RMP was one of only four DCPs uncovered herein involved in mRNA processing and editing (7.5%), the most temperature-sensitive cellular process documented previously in S. hystrix [17] (16%). In contrast to the variable temperature study, proteins whose concentrations differed between the two study sites herein were more likely to be involved in the stress response, transcription, transport, and the cytoskeleton.

Cell and molecular biology-driven approaches have aided in developing our understanding of both the fundamental biology of anthozoan-dinoflagellate endosymbioses [31] [32] [33] [34] , as well as their environmental physiology [35] [36] [37] . Herein we utilized a differential proteomics approach to uncover proteins whose concentrations differed between corals of an upwelling site (HBH) and a non-upwelling site (HWN). In addition to the difference in the mean monthly temperature range between these two sites uncovered in prior works (e.g., [12] ), a meta-analysis of previously published data conducted herein also revealed significant differences in silicate concentrations, suspended solid levels, and turbidity between them; specifically, these parameters were all significantly higher at HBH. As the PAR levels reaching the sampled corals were identical at each site, the elevated turbidity and suspended solid levels at HBH did not evidently result in a decrease in PAR reaching the sampled S. hystrix colonies. However, it is possible that such suspended solids affected other wavelengths of light that were not measured; for instance, high levels of suspended particulate matter lead to the attenuation of ultraviolet radiation (UVR) in wastewater [38] .

If diminished UVR levels were reaching the HBH S. hystrix colonies, then they could be hypothesized to be at lower risk of high temperature + UVR-induced bleaching [39] . In fact, relatively more proteins in the HBH>HWN differentially expressed host + Symbiodinium proteome (22%) were involved in the stress response compared to the 35-protein HWN > HBH proteome (3%). These three, presumably stress-targeted HBH > HWN proteins included two host coral proteins, beta-gamma crystallin and O-aminophenol oxidase, as well as the Symbiodinium protein known as HSPB1-associated protein 1-like. Beta-gamma crystallin, which has only ever been hypothesized to be involved in the stress response (as well as calcium binding [40] ), was actually found to be down-regulated in S. hystrix specimens exposed to a variable temperature regime for one week [17] and represents one of only four proteins found to differ in concentration between sites of origin herein and across temperature treatments in Mayfield et al. [17] ; as such, and in reiterating the recommendation of Meyer et al. [41] , we advocate that the role of beta-gamma crystallin in coral thermal adaptation be more thoroughly characterized in future works.

Given the fact that mean monthly temperature range, silicate concentrations, and suspended solid levels and turbidity were all higher at HBH (the upwelling site), it is tempting to speculate that the relatively higher number of stress response-associated proteins within the HBH > HWN proteome is indicative of sub-lethal levels of cellular stress in corals of HBH. However, the physiological performance of these corals did not differ from those of HWN ( [12] [15] and Table 1). This may mean, in contrast, that such constitutively elevated stress protein levels instead represent “front-loading” [25] [42] , whereby stress proteins are maintained at high intracellular concentrations such that the molecular machinery requisite for dealing with large shifts in the abiotic environment (e.g., an upwelling event in the case of corals of HBH) is engaged at any given time. This strategy is relatively uncommon in nature given the high energetic expense of being constitutively “stressed” at all times [43] , though it does characterize the cellular biology of some intertidal invertebrates [44] . By exposing corals from both study sites to variable temperature regimes and elevated suspended solid levels in the laboratory, it may ultimately be determined whether the protein expression signatures of corals of HBH documented herein are resultant of cellular stress due to, for instance, direct impacts of sediments with coral tissues (sensu [45] [46] ), or, alternatively, represent a protective response to counter future environmental change. If healthy corals indeed constitutively maintain high cellular concentrations of stress-associated proteins, as has also been documented in the most-remote reaches of French Polynesia [35] , then this may complicate the interpretation of data derived from molecular biomarker panels aimed at assessing coral health; in other words, are high concentrations of molecular chaperones and other stress proteins indicative of healthy corals or stressed ones?

Osmoregulation has been hypothesized to be the crux of the coral stress and bleaching response [47] . This theory stems from the fact that high temperature- and light-induced photoinhibition [48] would presumably result in a reduction in osmolyte flux from Symbiodinium to host. This would lead to a drop in osmotic pressure in the endosymbiotic gastrodermal cell [49] , which would manifest in changes in cytoskeletal architecture [50] . This may explain why a large number of host cytoskeleton-associated proteins were differentially concentrated between corals of the two sites. Specifically, actin and gelsolin (which is involved in actin assembly/disassembly) were synthesized at higher levels by corals of HBH. The fact that these corals are more likely to undergo osmotic pressure fluctuations as a result of upwelling-induced temperature changes may explain why they constitutively synthesize larger quantities of gelsolin, in particular, as this protein acts in cytoskeleton rebuilding. Not surprisingly then, genes encoding proteins involved in the cytoskeleton have found to be differentially expressed in corals exposed to elevated temperatures [51] . Given these findings, as well as the fact that osmoregulation is the most energetically expensive task a cell undertakes, the role of osmoregulation in the coral thermal stress response should be more thoroughly characterized in future works. It should be noted here that, due to the rigid cell walls of Symbiodinium, the osmotic stress associated with temperature + UVR stress-derived photoinhibition mentioned above is not hypothesized to dramatically affect Symbiodinium cell volume [47] ; therefore, it is unsurprising that no Symbiodinium cytoskeleton proteins were uncovered herein.

Unfortunately, it has become commonplace in the coral biology field to make conjectures about protein behavior based on mRNA expression data alone (sensu [52] [53] ), despite the fact that mRNA vs. protein congruency has been found to be as low as 0% in Symbiodinium populations within S. hystrix nubbins of another study [17] . Likewise, the degree of congruency between mRNA expression and protein concentration was markedly low herein; only 4 of the 53 DCPs uncovered were associated with an mRNA whose expression also differed significantly between sites of origin (7.5%). One such molecule, the host coral O-aminophenol oxidase, was discussed above in the context of the coral stress/environmental acclimation response. Given the low congruency between mRNA expression and protein concentration in both this study (7.5%) and others (2% in Mayfield et al. [17] and 10.5% in Mayfield et al. [26] ), we advise that those researchers looking to model the response of anthozoan-dinoflagellate endosymbioses to environmental change instead exploit proteomics-based approaches in their experiments; unlike the mRNAs that encode them, proteins actually enact physiological changes in cells and are likely to play key roles in thermal acclimation in both coral hosts [54] and their in hospite Symbiodinium populations.