Divergent proteomic responses offer insights into resistant physiological responses of a reef-foraminifera to climate change scenarios

: Reef-dwelling calcifiers face numerous environmental stresses associated with anthropogenic carbon dioxide emissions, including ocean acidification and warming. Photosymbiont-bear-ing calcifiers, such as large benthic foraminifera, are particularly sensitive to climate change. To gain insight into their responses to near-future conditions, Amphistegina lobifera from the Gulf of Aqaba were cultured under three p CO 2 conditions (492, 963, 3182 ppm) crossed with two temperature conditions (28°C, 31°C) for two months. Differential protein abundances in host and photosymbionts were investigated alongside physiological responses and microenvironmental pH gradients assessed via proton microsensors. Over 1000 proteins were identified, of which >15% varied significantly between treatments. Thermal stress predominantly reduced protein abundances, and holobiont growth. Elevated p CO 2 caused only minor proteomic alterations and color changes. Notably, pH at the test surface decreased with increasing p CO 2 under all light/dark and temperature combinations. However, the difference between [H + ] at the test surface and [H + ] in the seawater – a measure of the organism’s mitigation of the acidified conditions – increased with light and p CO 2 . Combined stressors resulted in reduced pore sizes and increased microenvironmental pH gradients, indicating acclimative mechanisms that support calcite test production and/or preservation under climate change. Substantial proteomic variations at moderate- p CO 2 and 31°C indicate cellular modifi-cations, and putative decreases in test stability at high- p CO 2 and 31°C, contrast the LBF s’ apparently stable overall physiological performance. Our experiment shows that the effects of climate change can be missed when stressors are assessed in isolation, and that physiological responses should be assessed across organismal levels to make more meaningful inferences about the fate of reef calcifiers.


Introduction
Anthropogenic carbon dioxide (CO2) emissions affect ecosystems worldwide [1]. Marine environments are impacted by both ocean warming (OW) and ocean acidification (OA) that result from rising atmospheric CO2. The associated rise in seawater temperature poses a major threat to coral reef ecosystems [2]. Alongside the altered carbonate chemis-try and saturation states [3][4][5], warming represents a substantial challenge for marine calcifiers to overcome [6] as it leads to the breakdown of the symbiosis between host organisms, like corals and foraminifera, and their photosymbionts. OA has various impacts on biomineralization and skeleton construction [7,8], with photosymbiotic reef organisms, such as scleractinian corals, often showing non-linear and even parabolic responses [9][10][11]. The responses of photosymbiotic organisms can be difficult to predict due to their differential yet potentially complementary modification of carbonate chemistry within their cytoplasm, but also in their microenvironment, during photosynthesis and calcification [12][13][14]. Moreover, global stressors are acting in concert, potentially compounding the physiological responses of shallow-water marine calcifiers [15], which adds to the complexity of predicting the fate of coral reef organisms under future global change scenarios.
Large benthic foraminifera (LBF) are protistan eukaryotes that build calcareous shells, known as tests, and host endosymbiotic microalgae. The proliferation of these and other calcifying tropical/subtropical reef organisms, including corals, are augmented by photosynthesizing endosymbiotic microalgae. The energetic advantages of this association arise from the photosynthate-derived carbohydrates that are released by the symbionts to the host [16,17], thus facilitating LBF to be major calcium carbonate producers within reef ecosystems and contribute significantly to the global carbon cycle [18,19]. The taxonomic family Amphisteginidae has occurred in the fossil record for the past 50-60 million years, and is one of the most scientifically recognized and widespread modern families of LBF [20]. Amphistegina spp. appear more resilient to thermal stress than many other LBF [21][22][23], which could be an adaptation from ancestors of the Paleogene, when atmospheric CO2 fluctuated substantially and was generally 2-to 3-fold higher than Neogene levels [24]. During the middle Eocene, Amphistegina spp. and other hyaline LBF were extremely diverse and dominant in shallow tropical marine carbonate environments [25], when atmospheric CO2 was similar to or higher than today [26]. Recently, 18S SSU rDNA markers have shown how rich diversification and broad geographic distribution have contributed to the withstanding of Amphistegina spp. to extinction events throughout evolutionary history [27], however, it remains unclear which adaptive traits exactly facilitated these photosymbiotic calcifiers to thrive throughout such variable environmental conditions and maintain their long-lasting evolutionary success. In modern settings, Amphistegina spp. are prevalent throughout the tropical/subtropical coastal environments [28,29], and are known to occur in such high densities that they can locally represent >1500 g per m 2 CaCO3 of marine sediments [30], contributing substantially to sand production and beach nourishment.
Generally, LBF seem to be more tolerant of and adaptive to ocean warming and acidification compared to corals, but the underlying mechanisms of variations in sensitivity to global change are not fully understood [31]. Ocean warming and high light intensities are known to cause LBF bleaching [32,33]. The reductions in photosymbiont densities or pigmentations are associated with impaired host performance, including reduced photosynthesis and respiration, motility, growth and antioxidant production (e.g. [23,34,35]). Such effects are often accompanied by increased proportions of malformed or broken tests, infections, and reproductive dysfunction [36]. Yet, there are LBF with remarkably high thermal thresholds, such as A. lobifera from the Gulf of Aqaba [22]-similar to many other marine species from the northern Red Sea [37,38]. The results of laboratory and fieldbased experiments on the effects of OA, however, range from predicting extinction of LBF by 2100 [39], and other negative impacts like test dissolution, reduced growth, lower chl a concentration and photosynthesis (790 -1100 µatm [21,[40][41][42][43][44]), to no significant effects on photosynthesis, chl a and net production/growth (up to 1925 µatm [41,45]) or positive effects such as increased symbiont performance or higher growth rates (580 -770 µatm or 1169 -1662 μatm [40,41,[44][45][46][47]).
The effects of the combination of these stressors on LBF physiology has received little study thus far [21,41,48]. Each compartment of the holobiont (photosymbiont or host) could react depending on the environmental conditions, functional group and life cycle of the organisms, and the test mineralogy [19], thus making the interactive effects of future global change on LBF hard to predict. Mass spectrometry-based bottom-up proteomics was recently proven to be a powerful tool for detecting cellular responses among Amphistegina spp. [49], especially as this approach allows quantification of changes in algal symbiont and/or the host foraminifera simultaneously. Application of proteomics approaches may thus reveal the bio-functional mechanisms involved in maintaining homeostasis and symbiosis under rapidly changing environmental conditions.
To test the tolerance of a globally distributed shallow-water LBF to the independent and combined effects of OA and OW, a crossed-factor laboratory experiment was set up. Amphistegina lobifera from the northern Red Sea were cultured under elevated temperature (T, control = 28°C and high = 31°C) and three pCO2 scenarios (ambient / control = 492 ppm, year 2100 / moderate = 963 ppm, and year 2500 / high = 3182 ppm [1]) for two months. Effects on its physiological performance, diffusive boundary layer (DBL) pH, and host and photosymbiont proteomes were investigated by measuring proton concentrations on the test surface and in the surrounding seawater during dark and light (via microsensors), holobiont mortality (as frequency of empty tests), growth (as changes in cross-sectional area area), and color (on the L*a*b color scale as proxy for photo-physiological condition), chlorophyll a concentration (as proxy for photosymbiont biomass), test pore size (on newest chamber), as well as changes in relative abundance of proteins associated with host foraminifera and/or their photosymbiotic diatoms (via liquid chromatography-tandem mass spectrometry [LC-MS/MS] proteomics). This combined approach focuses on detecting molecular alterations at the proteomic level, representing mechanisms that contribute to holobiont functioning, as well as the downstream functional performance development of LBF to future climate change. Elucidating the cellular mechanisms responsible for acclimatization and resilience of photosymbiotic calcifiers will ultimately contribute to our understanding of how coral reefs evolve, aiding in their conservation and restoration.

Materials and Methods
Amphistegina lobifera harbor diatoms as endosymbionts [22,[50][51][52]. Adult specimens of A. lobifera reach sizes of >1 mm in diameter [53,54], the spherical perforate test is built of intermediate-Mg calcite, ca. 4-6 mol% MgCO3 [55]. At depths <10 m A. lobifera is the shallowest dwelling foraminifera of the genus and is commonly found attached to coral rubble in the upper reef slope, reef flat, and back reef zones [56,57]. This genus has a broad geographic distribution and specimens from the Red Sea, along with populations from the western Indian Ocean and the Mediterranean Sea, belong to lineage Ia [27].

Sample collection
Collection procedures followed established protocols [58]. Initially, coral rubble was collected from the narrow and patchy fringing reef (29°30.619 N, 34°55.247 E) in the northernmost extension of the Red Sea, the Gulf of Aqaba (GoA) [28,59]. Samples were obtained by snorkeling in 2−4 m deep water in front of the Inter-University Institute for Marine Sciences (IUI) in Eilat, Israel, in March 2016. LBF specimens were removed from the rubble and transported to the Leibniz Centre for Tropical Marine Research (ZMT), Bremen, in thermally insulated containers for <12 h. Temperature loggers (HOBO®, USA) were placed in the containers at a temperature of 23.3°C (standard deviation (SD) = 1.3), which is within the normal daily temperature range for the species. Upon arrival at the marine experimental facility (MAREE) at ZMT, the samples were cultured under low-light intensities of 15−20 μmol photons m -1 s -1 (LICOR LI-250) to recover from transport. Since the water temperature during collection at springtime (ca. 23°C) was relatively low compared to the target-temperatures of the experiment (28, 31°C), the specimens were slowly acclimated to the experimental target temperatures by exposing them to a temperature of 24−25.5°C during the first week. After this, the foraminifera were distributed amongst the experimental aquaria initially maintained at 26°C, after which the temperature was ramped up by 1°C per day until the final temperature levels of 28°C or 31°C were reached. Although mean annual temperatures at the collection site in the GoA rarely reach 28°C, that temperature was selected as the control because a prior study showed that growth and photosynthetic performance of A. lobifera from the GoA were not impacted by a threeweek exposure to temperatures of up to 30°C [22]. As their thermal threshold appears to be slightly above this value, and because global sea-surface temperatures are expected to rise by an average of 3°C (RCP 8.5 [1]), 31°C was used as the high-temperature treatment. These experimental temperatures also correspond to the summertime thermal range at the sampling location (Fine et al., 2013). Temperatures were maintained at 28°C (SD 0.2) and 31°C (SD 0.3) using 125-watt aquarium heaters (EHEIM).

Experimental design
A crossed-factor laboratory experiment with six treatments, two T levels (control = 28°C, high = 31°C), and three pCO2 scenarios (ambient / control = 492 ppm; year 2100 / moderate = 963 ppm; year 2500 / high = 3182 ppm [1]) were replicated four times (n = 24 tanks in total), was set up. The culturing array consisted of six 284-L seawater reservoirs, which were comprised of two seawater basins and four 10-L replicate tanks. Specimens were reared for eight weeks in each of the 10-L aquaria (≥ 18 specimens per replicate tank, > 432 specimens in total). This study was conducted in parallel with scleractinian coral specimens in the same experimental array. Corals held in the same experimental array [15,60] were fed live Artemia sp. hatched from ~40 mg of eggs within 24 hrs every second day, which also supplied sufficient organic matter and nutrients to the foraminifera maintained in the same tanks, especially since A. lobifera has previously been shown to depend very little on heterotrophic feeding [61,62]. Natural oligotrophic seawater collected from Svalbard, Norway, was used.
Foraminifera were housed in 48 custom-made transparent plastic containers cut open on both sides (approximately 1 × 2 × 5 cm) and covered with 400 µm nylon mesh fixed by coral glue (EcoTech Marine, USA). This arrangement allowed water to circulate through the submerged containers. Two of these containers were placed in each of the 24 10-L replicate tanks (4 replicate tanks per treatment). One of the two containers held eight specimens that served for assessing growth and mortality. The other container held at least ten specimens, one to two of which underwent pH microelectrode measurements, chl a extractions, and scanning electron microscopy (SEM) imaging, while the remaining eight specimens were used for proteomic analysis.
The tanks were illuminated for 12 h per day with compact fluorescent lights (ultraactinic and white; 96 W, 1000 K) combined with standard white fluorescent lights (32W, T8, 6500 K), with a maximum photosynthetically active radiation (PAR) of ca. 250 µmol photons m -2 s -1 to meet the light intensity needed for the zooxanthellate corals cultured in the same experiment. To account for the needs of A. lobifera, the irradiance was reduced by a factor of ten (to 25−35 μmol photons m -1 s -1 ) by shading with the black nylon mesh. Light measurements were conducted using a quantum light meter (LI-COR LI-250A) with a submersible spherical micro-quantum sensor (Waltz, Germany).
The seawater in all of the experimental tanks was equilibrated with mixed gases formulated at three pCO2 levels. The control-pCO2 treatment was set at 492 ppm (SD 80) to emulate control (C) present-day conditions, the moderate-pCO2 treatment (M) was set at 963 ppm (SD 134) to emulate predicted year-2100 conditions, and the high-pCO2 treatment (H) was set at 3182 ppm (SD 817) to emulate predicted year-2500 conditions [1,5]. To achieve the target pCO2 conditions, CO2-free air and pure CO2 gas were mixed with solenoid-valve-controlled mass flow controllers and bubbled via flexible microporous air stones into the two large seawater reservoirs of each experimental treatment group.
Seawater was continuously pumped from the lower large reservoir to the four replicate tanks contained within the upper large reservoir of each treatment group using a 30watt aquarium pump (output = 2700 L hr -1 , Aqua Medic Eco Runner 2700). Seawater was distributed to each replicate tank via two horizontal pipes with holes drilled along their long-axes and placed at both ends of the tank, one at the top of the tank and one at the bottom, so that the corals and foraminifera specimens received laminar, circulating flow. Seawater overflowed the replicate tanks into the large upper reservoir then into the lower reservoir where it was filtered by activated charcoal (Deltec brand protein skimmer) and a filter medium pad. The seawater was continuously recirculated from the lower reservoir to the replicate tanks in the upper reservoirs using an aquarium pump. Seawater from the reserve tank was added to each experimental treatment group at a rate of 0.6 L/hour, which allows seawater to mix at near-equilibrium within each tank, thereby maintaining target carbonate chemistry of the experimental seawater's dissolved constituents via calcification and photosynthesis (e.g., Ca 2+ , alkalinity, nutrients). Each treatment group was covered with a transparent Plexiglas lid to prevent in-/outgassing from the seawater and to minimize evaporation.
Following the methods of Guillermic et al. [15], salinity, pH, and temperature were measured three times per week in all of the replicate tanks using a combination probe kit (WTW Multi 3430 Set K). The pH and salinity sensors (Thermo Orion TM ) were calibrated at the appropriate tank temperature (28°C or 31°C ) before each measurement using NISTtraceable pH 4 and 7 NBS buffers and Certified Seawater Reference Material (Dickson standard batch #154) of salinity 33.347, respectively. Total alkalinity (TA) and dissolved inorganic carbon (DIC) samples were collected once per week at midday, thereby minimizing disparities in DIC due to diurnal variations in the photosynthetic activity and respiration of the tank inhabitants. DIC samples were collected in sealed 25 ml borosilicate glass vials and TA samples were collected in 50 ml centrifuge tubes that were twice flushed with seawater prior to use. Subsequently, DIC and TA samples were preserved with 10 µm and 20 µm, respectively, of saturated mercuric chloride solution and refrigerated until analysis. TA was measured by open-cell potentiometric Gran titration (accuracy 20 mol kg -1 ) and DIC was measured with a Shimadzu DOC analyzer converted for DIC analysis (accuracy 10 mol kg -1 ). Every seventh sample, the TA and DIC instruments were calibrated with Certified Seawater Reference Material (Dickson Batch #154). Concentration of carbonate ions, bicarbonate ions, aqueous CO2, aragonite (ΩA), and calcite saturation states (ΩC) were calculated from measured TA, DIC, salinity, and temperature with CO2SYS software [63] using Roy et al. [64] values for the K1 and K2 carbonic acid constants, the Mucci [65] value for the stoichiometric aragonite and calcite solubility products and with an atmospheric pressure of 1.015 atm.

Microelectrode measurements
Microsensors are commonly used to characterize the DBL that exists between the surface of an organism and the surrounding seawater (e.g. [7,12,66,67]). Among foraminifera, oxygen and pH gradients of the DBL reflect photosynthesis, respiration, and biomineralization activity [12][13][14]. To investigate the DBL of LBF, liquid ion exchange (LIX) pH microelectrodes were fabricated in the de Beer laboratory (Max Planck Institute for Marine Microbiology, Bremen). Microsensors were produced following de Beer et al. [68], but were filled with hydrogen ionophore II cocktails A and B (Sigma Aldrich), with cocktail B having an epoxy added to it to prevent leakage of the ionophore out of the sensor. Microelectrode measurements were performed in two reservoirs of seawater, sourced from the corresponding experimental treatment tanks, located adjacent to the experimental tanks and the microelectrode system. These reservoirs were bubbled with the corresponding treatment gases and maintained at the corresponding treatment temperature using 125-watt aquarium heaters (EHEIM, 2015), mimicking the respective treatment conditions in the experimental aquaria. The seawater was circulated between the two reservoirs through two 5.4 L flow-through chambers (30 x 12 x 15 cm). All pH microelectrode measurements were performed within these smaller flow-through chambers, with gentle water movement maintained with a small powerhead (12V Tarwig). Following the methods of Glas et al. [12], single specimens were carefully removed from the experimental culturing setup and placed on the bottom of aquaria. LIX membrane microelectrodes (8-20 μm tip diameter, ≤20 s response time) connected to a high-impedance mV meter were mounted on a motorized micromanipulator (Micos VT-80). Prior to measurements, the microelectrodes were prepared and calibrated in pH buffer solutions (Mettler Toledo, pHNBS 7.0 and 9.21) at treatment temperature, following de Beer et al. [68]. Microelectrode tips were positioned on the LBF test surface. Exact positioning of the specimen and the microsensor tip was achieved using a Stemi SV6 binocular (Zeiss).
To determine the effects of photosynthesis and respiration on the pH microenvironment under the different pCO2 and temperature treatments, the gradient in pH away from the test surface was measured in three specimens per treatment. The foraminifera were allowed to acclimate to the conditions and recover from handling for five minutes before initiating pH measurements. Measurements were first conducted under light until a stable signal was reached, after which the light was turned off. When pH at the test surface reached a constant level, the light was turned on again and measurement was continued until a stable measurement was reached. This was followed by moving the microelectrode along a 1000-µm vertical profile away from the LBF test surface. A step size of 1 µm was used for the first 100 µm in order to generate high-spatial-resolution pH profiles within the DBL at the surface of the test. After this, step size was increased to 20 µm in order to characterize pH of the seawater adjacent to the DBL. Between measurements, a fine brush was gently applied to the LBF to prevent it from articulating and moving away from the sensor. The difference in [H + ] between the DBL (immediately adjacent to test surface) and the seawater ('∆[H + ]') was calculated as an indicator of the specimen's mitigation of the acidified conditions.

Holobiont mortality and growth
At the end of the 8-week experiment, the proportions of live-to-dead specimens in each treatment were recorded. Empty tests, determined by lack of cytoplasmic color and reticulopodial activity [69], indicate either stress-related death, natural death (old age), or reproduction in which the host's cytoplasm undergoes either asexual or sexual, leading to the release of hundreds of juvenile cells into the environment. Based on test size homogeneity, it is assumed all specimens used were from the same cohort. Additionally, there was no evidence of reproduction throughout the experimental period (i.e., no juvenile cells). Therefore, the occurrence of empty tests is assumed most likely to be due to a response to stress resulting in mortality. Mortality rate was calculated as a percentage of all specimens (n = 195) from each treatment.
A subset of LBF from each treatment (n = 4 replicates per treatment) was photographed with a digital light microscope (VHX-5000, Keyence) using fixed settings immediately before and after the experiment (as described in [23]). The surface areas of all specimens were determined with the image analysis software Fiji [70]. Growth was calculated as the average daily increase in cross-sectional surface area of all specimens from each sample over the eight-week experimental period following ter Kuile & Erez [54].

Color and chlorophyll a
Holobiont color, as a proxy for photo-physiological condition, was determined using the International Commission on Illumination (CIE) L*a*b* color space values of each specimen following Hosono et al. (2012). In each image, holobiont color and background color were transformed into CIE L*a*b* color space by using the color space converter in Fiji [71], as described in Stuhr et al. [23]. The resulting color values represent the three coordinates within the CIE color space: L* indicates whiteness (brightness) of the color (0 = black, 100 = white) and was shown to correlate to photosynthetic performance as measured via oxygen production rate [72], while a* indicates the position between green (negative values) and magenta (positive values), and b* indicates the position between blue (negative values) and yellow (positive values), serving as indicators for photopigment content. A total of 175 specimens (i.e., those that were alive at the end of the experiment) were used for color measurements.
Chl a concentration was analyzed as a proxy of photosymbiont biomass, adapted from established protocols [32]. Single specimens from three randomly chosen tanks from each treatment, which were also used for microelectrode measurements (n = 3 per treatment), were each placed into 320 µl of 90% ethanol to extract pigment, heated to 78°C for 5 min, vortexed, and extracted in the dark for 24 h at 4°C. Then 300 µl of the extraction liquid was measured at absorbances of 630, 663, and 750 nm using a microplate reader (Tecan infinite M200 PRO). After extraction, specimens were dried at 40°C for 24 h and weighed with an analytical balance (Mettler Toledo MX5). Chl a concentrations were calculated following the method of Jeffrey and Humphrey (1975) and normalized to dried test weight [23]. The same procedure was repeated with the remaining specimens (n = 1 -4 per replicate tank / 7 -10 per treatment) to ensure that the single specimens selected for microelectrode measurements represent the average photosymbiont condition of their respective treatment. The results of prior total chl a contents were added to the respective replicates of these measurements and normalized by their collective dry weight to calculate average chl a concentrations.

Test surface, pore size and breakage
Images of all remaining LBF 'tests' were made by scanning electron microscopy (SEM) to detect dissolution or other alteration of the test surface. The dried samples were mounted on conductive carbon discs attached to specimen stubs and sputter-coated with gold for 35 sec under vacuum. Images of the surface morphology of the mounted specimens were made by secondary electron (SE) microscopy (Tescan VEGA3 XMU), initially at the highest magnification that maintained the entire test in the SEM's field of view, and then at ca. 2250-fold magnification. As Amphistegina add a new layer of calcite to their entire test surface during chamber addition (growth), these high magnification images revealed the newest growth of each specimen. Dissolution of each specimen's test surface was assessed visually by comparison with tests formed under the control condition and also by comparison with images of LBF reference samples that were cultured under control conditions. Representative sets of images were assembled for each treatment. Changes in test microstructure were quantitatively assessed by comparing average pore sizes determined in SE images (n = 72) using the function 'Analyze particles' in ImageJ (Version 1.52f; method modified from [73]). The occurrence of partial test breakage visible in SE images, possibly induced during handling, was counted and its percentage in each treatment was calculated.

Statistics
Permutational multivariate analyses of variance (PERMANOVA) were conducted on all ∆[H + ] data, using PAST 4.04. Two-way analysis was used for the fixed effects of pCO2 and T, and a one-way analysis for the fixed effect of light (vs. dark). This was followed by pairwise comparison of ∆[H + ] data amongst treatments, with p-values adjusted for multiple testing following Benjamini & Hochberg [74] at a false-discovery-rate (FDR) q-threshold of 0.1. For statistical analyses of physiological parameters, residual and normality plots indicated that assumptions of equal group variance and normality were not violated for chl a concentrations, L*a*b* values, growth rates, or pore sizes, and the Shapiro-Wilk test was used to confirm normal distribution of these data. Changes in these parameters under the respective fixed factors T and pCO2, or their interaction, were analyzed using linear models for two-factorial analysis of variance (ANOVA). In cases where the linear model indicated significant effects of either of the factors, or a significant interaction between the factors, a posthoc Tukeys' Honest Significant Difference (HSD) test was conducted. For the binary mortality (dead) and breakage data, a nominal logistic model was used with the fixed factors T and pCO2 or their interaction. All analyses of physiological parameters were conducted in JMP Pro 14.3.0.

Identification of photosymbionts
LBF are known to harbor a wide range of microalgae. To identify the dominant photosymbiont species in the A. lobifera population from Eilat, a 442-base pairs (bp) fragment of the 3' end of the 18S rDNA was extracted from six specimens, amplified, and sequenced as described in Stuhr et al. [35]. The resulting diatom sequences were deposited at European Nucleotide Archive (ENA) under the study accession number PRJEB33832. Sanger sequences were trimmed and assembled using Geneious 8.1.9. All resulting consensus sequences were identified via BLASTn [75] searches using the default settings for highly similar sequences (megablast) against the NCBI database (accessed December 2018) and several representative sequences (with an e-value below 1e-166) were manually selected from the results. These sequences were complemented with foraminifera symbiont sequences from [22,35], and other closely related algae formerly reported as photosymbionts in LBF (e.g. in Lee et al. [50]). The collected sequences were aligned using MAFFT v7.310 [76] with the default options and trimmed to equal length at the 3' and 5' prime end in BioEdit [77]. Subsequently, identical sequences were removed using Alter [78] and the collapse haplotypes command. The final 407-bp alignment was used to calculate a Maximum Likelihood phylogenetic tree based on the GTR model with gamma distribution and 1000 bootstrap replicates in Mega 6.06 [79]. The final tree is shown with a 50% bootstrap cut-off ( Figure S1).

Label-free peptide quantification and protein identification
To gain insight into the underlying cellular responses of the photosymbionts and the host foraminifera, a label-free quantitative proteomics approach based on liquid chromatography coupled to tandem mass-spectrometry was applied, as described in [49] and provided in appendix A. At the end of the experiment, eight or nine pooled specimens from three replicate tanks for each treatment (i.e., >144 specimens in total) were frozen at -80°C and used for the proteomics analysis. As three out of the initial four biological replicate samples from the control treatment (control pCO2 and 28°C = 'C28') were damaged during sample preparation due to defective laboratory equipment, this treatment had to be removed from further statistical analysis. However, technical replicates of the remaining C28 sample, as well as biological replicates from samples that were collected before the temperature ramping of the experiment (t0), were measured. These data are provided along with all final treatment data in the ProteomeXchange Consortium via the PRIDE partner repository [80] (dataset identifier PXD009890).
A homology-based search approach allowed peptide sequences to be assigned to either the host or photosymbiont component and identified at the protein level. An updated and enlarged reference dataset for foraminifera, containing the recently published genome sequences of the calcareous benthic foraminifera Globobulimina [81] was used, combined with the publicly available genome data of three diatom species used previously [49]. Clusters of homologous protein sequences were generated from amino acid sequences of all proteins identified in the database search with a sequence identity cut-off of 0.7. Normalized abundance values of the peptides were used to calculate protein cluster abundances. Subsequently, protein descriptions and gene ontology (GO) annotation were obtained via eggNOG-mapper v2 [82,83] using default settings (minimum hit e-value = 0.001, minimum hit bit-score = 60, minimum % of query coverage = 20).

Differential protein abundance and functional annotations
To identify proteins with significantly different abundances between the treatments, multivariate analyses with post hoc identifications were conducted in R 4.0.0. [84], adapting procedures of proteomics studies on anemones [85,86]. To aid in parametric modelfitting, data were log-transformed, and Bray−Curtis dissimilarities were calculated for each pairwise sample comparison to identify the differentially abundant protein-clusters. Analysis of variance by permutation of dissimilarities was conducted using the "adonis" function in the R package "vegan" version 2.5-6 [87]. First, a PERMANOVA, with fixed effects of pCO2, T, and their interaction, was conducted on all proteins, only host-associated proteins, and only photosymbiont-associated proteins. This was followed by pairwise comparisons for significant differences using the "pairwise.perm.manova" function in the R package "RVAideMemoire" [88], with p-values adjusted for multiple testing at a FDR of 0.05 [74]. All PERMANOVAs used 9999 permutations. Then, a generalized linear model (GLM) with fixed effect treatments was fitted and compared against a null model (intercept only) using an analysis of deviance  2 test with α = 0.05. Differentially abundant proteins (DAPs) were identified with the FDR procedure [74] at a q-threshold of 0.05, correcting for multiple testing. The corresponding log2 fold-changes (log2FC) were derived from average normalized abundance values between treatments.
Further descriptions of DAPs were extracted via OmicsBox 1.3.11 (BioBam Bioinformatics, 2019) by searching via the blast program blastp-fast against all organisms deposited in the public NCBI non-redundant protein sequence database on July 9 th 2020 (evalue-hit-filter: 0.001, annotation cut-off: 55, HSP length cut-off: 33, number of hits: 20).
Any proteins without a match were designated as hypothetical proteins. For proteins without description in the non-redundant NCBI database, those given by the Swiss-Prot database or eggNOG annotations were used, if available. DAPs were grouped by major functional categories inferred from GO term annotations describing molecular functions, biological processes, and cellular components, though it should be noted that most of them can function in multiple ways and are not restricted to these categories.
A Fisher's Exact test was conducted to determine whether DAP frequencies differed between host and photosymbionts, and multivariate analyses of log-transformed abundance values of DAPs in JMP Pro 14.3.0. Hierarchical cluster analyses of proteins and samples or treatments, using average distance linkage, as well as principal component analysis (PCA) on row-wise correlations, were performed. Further PCAs were run on total normalized abundance values of each functional category of host-or photosymbiont-associated DAPs. To further focus on those most 'regulated' proteins that varied significantly in abundance between treatments and also showed pronounced fold-changes, only DAPs with log2FC ≥|1| between those treatments with Tukeys' p ≤0.05 were extracted and highlighted.

Physicochemical parameters of seawater treatments
A summary of the measured aquarium seawater parameters T, DIC, TA, calculated seawater pH on seawater-scale (pHSW), calcite (ΩC) and aragonite saturation states (ΩA) of the treatment seawater, and pCO2 of gas in equilibrium with seawater are provided in Table 1. Salinity and carbonate ion concentrations are found in the supplementary material (Table S1). In all treatments, temperature was maintained at either 28 ± 0.2°C (± indicates SD) or 30.9 ± 0.3°C, salinity stayed at 35.6 ± 0.6, and seawater remained supersaturated with respect to aragonite (ΩA >1) and calcite (ΩC >1). Calculated pHSW of the experimental treatments was 8.16 ± 0.11, 7.97 ± 0.07, and 7.50 ± 0.11 in the control (491 ppm), moderate (963 ppm), and high (3192 ppm) pCO2 treatments, respectively. Table 1. Mean physical and chemical conditions ± SD in the different treatments of elevated pCO2 (C = 492 ppm, M = 963 ppm, H = 3182 ppm) crossed with 28°C and 31°C. Temperature = T, dissolved organic carbon = DIC, total alkalinity = TA, pH on seawater scale = pHSW, pCO2 of gas in equilibrium with seawater = pCO2, calcite saturation state = ΩC, and aragonite saturation state = ΩA. Additional data on water conditions are provided in Table S1.

Identification of photosymbionts
The photosymbiont sequences recovered from these A. lobifera represent previously non-sequenced taxa, but phylogenetic analysis as well as BLAST searches confirm that the photosymbionts are a member of the same diatom assemblage as described previously from A. lobifera and other amphisteginids [22,35,52,89]. Molecular fingerprinting of symbionts in the present study revealed that the directly sequenced A. lobifera specimens (n = 6) contained six slightly variable haplotypes ( Figure S1) clustering closely with sequences from A. lobifera photosymbionts, collected from the Red and the Mediterranean Sea [22], within the order Fragilariales.

Microelectrode measurements
Microsensor measurements showed that pH at the test surface (pHts) decreased with increasing pCO2 under all light/temperature conditions (Figure 1b) with light and temperature. The most-elevated pH at the test surface in comparison to the surrounding water was consequently found at high-pCO231°C under light (Figure 1b).

Color and chlorophyll a
There was no interactive effect of pCO2 and T on L*, a*, or b* values, nor on chl a concentration ( Figure 2, Table S4). However, there was a significant effect of pCO2 on b* (ANOVA p = 0.040, Figure 2d). Post-hoc tests suggest a marginally significant (Tukeys' p = 0.053) color shift from yellower to bluer tones at high-pCO2 compared to control-pCO2 regardless of temperature.

Holobiont mortality and growth
Holobiont mortality remained low throughout the experiment ranging between 2.9% in the control (28°C, 492 ppm) and 15.2% at control-pCO231°C ( Figure 3c). There was no significant effect of pCO2 or T on mortality rates throughout the experiment (p ≥0.136; Table S2).
There was no interactive effect of pCO2 and T (p = 0.174) on test growth, however, there was a significant reduction in growth under elevated T (p = 0.038, Table S4, Figure  3a).  Table S2). However, the p-value for the interaction term is marginally significant (defined as 0.05 < p < 0.07), suggesting a possible interaction between pCO2 and temperature on breakage frequency.
Mean pore size was significantly impacted by the interaction of pCO2T (p = 0.015, Figure 3b; Figure 2). Tukeys' HSD post-hoc analysis shows that at 31°C, under moderate and high pCO2, pore size was significantly reduced when compared to control-pCO2 conditions (Tukeys' p = 0.030 and p = 0.040, respectively; Table S4).  (28°C, a-c; 31°C, d-f). Entire specimens shown via light microscopy (upper rows) and low-magnification SE imaging (middle rows), with the terminal chamber shown via highmagnification SE imaging. Scale bars in the full-specimen images represent 500 µm, while bars in high-magnification images represent 50 µm.

Label-free peptide quantification and protein identification
A total of 1103 protein clusters (i.e., groups of proteins with ≥70% sequence identity) were detected and identified from A. lobifera and its photosymbiotic diatoms (Table S5). Fifteen of these clusters were comprised of specific protein sequences belonging to tubulins, actin, heat shock proteins (HSPs), and peptidyl-prolyl cis-trans isomerase, from both the host and photosymbionts. Approximately 40% of the protein clusters (hereafter referred to as 'proteins') were assigned to the photosymbionts, whereas 60% were associated with the host foraminifera (Figure 5a). The only successfully recovered sample from the control treatment (control-pCO228°C) was removed from further analyses due to the lack of replication. Permutational multivariate analyses of variance conducted on all proteins, as well as only those associated with either host or photosymbionts, all indicated significant single effects of pCO2 or T (2-way PERMANOVA, p < 0.030 and 0.018, respectively, Table 2) and no interactive effects. However, adjusted p-values of post-hoc pairwise comparisons were non-significant. Nevertheless, the effect of T was overall more significant than the effect of pCO2, especially for the photosymbionts. We employed a GLM using the five remaining 'treatment' combinations as single factors to detect significant differences in protein abundance. In total, 173 proteins were significantly differently abundant between treatments (q <0.05, Table S6). Of these DAPs, 55 were from the symbiont, 113 were from the host, and 5 were assigned to both ( Figure  5a). Hence, the proportion of differentially abundant host proteins was significantly higher than that of photosymbiont proteins (2-tailed Fisher's exact test, p = 0.049). Most significant differences were found between moderate-pCO231°C in contrast to all other treatments in both compartments, and least at control-pCO231°C (Figure 5b, Table S7). Second highest percentage of DAPs in the host was found in response to high-pCO228°C, while it was comparably low in the photosymbionts, where high-pCO231°C caused the second most DAPs.
Multivariate analyses of all DAPs confirmed this distinction of the moderate-pCO231°C treatment, and further illustrated high similarities among protein abundance patterns of control-pCO231°C and high-pCO231°C treatments, as well as moderate-pCO228°C and high-pCO228°C treatments (Figure 5c). The strong overlap between samples of these treatments led to partial clustering of alike samples ( Figure S2). However, the distinction of moderate-pCO231°C from all other treatments and the otherwise main grouping by temperature were evident throughout, as well as when host and photosymbiont DAPs were analyzed separately ( Figures S3 & S4). Principal component analysis (PCA) showed that 35.9% of overall variations were explained by the differentiation along the x-axis, which mainly differentiated moderate-pCO231°C from all other treatments, while 24.1% were explained by the separation of temperature treatment groups along the y-axis (Figure 5d). In moderate-pCO231°C, most proteins were less abundant, but those clustering close to this treatment were highly abundant and mostly host-associated. The differentiation between 31°C and 28°C was associated with overall reduced protein abundances in the host and, particularly, the photosymbionts, as most proteins clustered in the upper half of the diagram (Figures 5d & S6).
The origin of the most abundant proteins, beta tubulin and polyubiquitin, could not be determined due to their strong homology among organisms. Also, three other highly conserved and abundant proteins [actin, heat shock protein (HSP) 90, and ADP-ribosylation factor] could not be clearly associated to either host or photosymbionts. Polyubiquitin increased with pCO2, but not at 31°C. In contrast, the other four proteins were most abundant in control-and high-pCO231°C, but not at moderate-pCO231°C. HSP90 and ADPribosylation factor, which are responsible for protein folding and transport, also showed increases at high-compared to moderate-pCO2 at 28°C and 31°C, respectively ( Figure S2, Table S6). All proteins of unclear origin were removed from downstream analyses. Separate analyses of host and photosymbiont-associated DAPs showed overall alike patterns of protein variations (Figures S3-S6). Differences were mainly found for moderate-pCO228°C, which was overall most similar to the other 28°C treatment in the host, but highly similar to control-pCO231°C in the photosymbionts. Generally, there was a further separation along the y-axis in the photosymbiont, showing more distinct proteome variations of moderate-pCO228°C and high-pCO231°C in comparison to the host.

Functional categories of photosymbiont-and host-associated DAPs
The gene ontology analysis of DAPs revealed the biological processes, molecular functions, and cellular components these proteins are involved in (Table S6, Figures S3,  S4). The 97 most strongly regulated proteins (i.e. DAPs that also showed log2FC ≥|1| between treatments with Tukeys' p≤ 0.05) were grouped by major functional categories based on GO annotations and summarized (Tables 3 and 4). PCAs of total DAP abundances per functional category revealed the underlying contribution of these groups driving the proteomic differences between treatments ( Figure 6). . Arrows indicate the direction and intensity of functional categories driving differences between the treatments, with triangles connecting the triplicate samples of each treatment (legend as in Figure 5).
In the photosymbionts, a separation along the x-axis explained 48.6% of the variation, mostly separating 31°C from 28°C treatments (Figure 6a). High T (31°C) dominantly caused decreases of proteins involved in light-harvesting and photosynthesis, proton transport or cell redox homeostasis, metabolic and (fatty acid or other) biosynthetic processes, but also cytoskeleton (development) and protein transport decreased, while some proteins involved in translation, protein folding, biosynthetic processes, and proton transmembrane transport showed increases. High-pCO231°C led to a decrease of enzymes involved in fatty-acid and other biosynthetic processes, whereas proteins responsible for protein folding increased. At moderate-pCO231°C, stronger decreases were found, specifically of proteins involved in photosynthetic light-harvesting, (carbohydrate) metabolism, biosynthetic, and transport processes. At moderate-pCO228°C, some proteins of photosystem I, or involved in metabolism (e.g., the citric acid cycle) and cytoskeleton development-related processes were relatively elevated (Table 3). Table 3. Regulated proteins of photosymbionts of Amphistegina lobifera in response to elevated pCO2 ('C' = 492 ppm, 'M' = 963 ppm, 'H' = 3182 ppm) at 28°C or 31°C, without the C28 treatment (GLM, q≤ 0.05, Tukeys' p ≤0.05, and log2FC ≥|1|). Proteins are arranged by their major functions and processes according to GO terms, however, these roles are not exclusive as many proteins can have various functions. Directions and intensities of log2FC are illustrated by color, significant differences in average protein abundance between treatments are denoted in bold. Respective annotated gene ontologies, accession numbers, and remaining DAPs are provided in Table S6. In the host, the impact of these functional protein groups on treatment differentiation was even stronger (Figure 6b). Elevated temperature (31°C) caused general decreases of ribosomal and other translation-related proteins, and various metabolic (catalytic) enzymes (e.g., those involved in the citric acid cycle, pentose-phosphate pathway, and electron transport). Likewise, some proteasomal and proteins involved in the regulation of protein folding or ubiquitination, cytoskeleton organization, or microtubule-based movement, decreased. In contrast, some proteins responsible for protein (re-)folding and degradation (proteasome subunits), and microtubule-based transport processes, increased, especially at moderate-pCO231°C. In comparison, proteins responsible for (transmembrane) ion transport or signal transduction were particularly increased at control-and high-pCO231°C. Several proteasomal, metabolic and biosynthetic proteins were increased with pCO2, especially at high-pCO228°C, whereas proton transporters seemed to decrease under that treatment (Table 4).

Discussion
LBF are important CaCO3 producers and ecological indicators integral to tropical marine habitats [19]. Increasing global CO2 emissions lead to simultaneous ocean acidification and warming, representing two major environmental threats to calcifying marine organisms. Understanding thresholds of physiological tolerance and possible acclimatization mechanisms to changing physicochemical conditions will help to better understand resilience mechanisms and the potential for persistence of coral reef organisms under near-future climate change. The geologic record of LBF throughout periods with high and strongly fluctuating atmospheric pCO2 and temperatures [24,25] has spurred interest in the abilities of photosymbiont-bearing foraminifera to withstand or adapt to rapidly changing ocean conditions. In our experiment, A. lobifera persists under the combined stressors of elevated pCO2 and T over an eight-week interval, showing only few effects induced by single stressor treatments or their interaction at the organism-level, but significant effects of temperature or pCO2 on the host and photosymbiont proteome.

Bleaching and Mortality
Consistent with some previous studies on A. lobifera [12,43], we did not see significant impacts of elevated temperature or CO2 (or their combined effects) on photosymbiont physiology as expressed by holobiont color and chl a (Figure 2). While ocean warming is currently considered the most harmful factor of climate change on coral reef organisms, our study documented a resistance to 31°C and elevated pCO2. Importantly, we did not observe any significant effects on A. lobifera mortality in response to either of these global stressors, nor to their combined effects (Figure 3c). However, similar studies have shown that elevated pCO2, in combination with high temperatures [21,41] or elevated copper concentrations [46], induces bleaching and lowered chl a concentration in other species of LBF. This suggests a species-specific and perhaps a study-specific response to significant interactions between multiple stressors on symbiont physiology. While elevated temperature alone has been documented to disrupt photosymbiosis as evidenced by reductions in photosymbiont density, pigment density, rate of photosynthesis, or coloration [23,32,35,90], it is important to note the thermal history of specimens used in other studies likely differed from that of the present study-in which LBF were collected from a site previously shown to harbor highly thermally resistant populations of LBF and scleractinian corals [22,37]. Ultimately, the present study is consistent with previous experiments demonstrating no or low mortality in various LBF species exposed to OA [41,46,47,91].
The proteomic response of the photosymbiont compartment offers interesting mechanistic insights into responses to OW/OA, showing these were primarily driven by elevated T and only secondarily by elevated pCO2. At 31°C, particularly at moderate pCO231°C, proteins involved in photosynthesis (e.g., fucoxanthin chlorophyll a/c protein, thylakoid lumenal protein 15kDa protein 2, PSI P700 apoproteins), several enzymes involved in (carbohydrate) metabolic processes (e.g., fructose-bisphosphate aldolase, alanine dehydrogenase, inorganic pyrophosphate), and various biosynthetic processes (e.g., uroporphyrinogen decarboxylase, 3-isoproylmalate dehydrogenase) were decreased. In contrast, proteins associated with unfolded protein response (e.g., HSP90, HSP70-type, ATP-dependent chaperone ClpB) were elevated at 31°C (Figures 5, 6a, S3, S5-S6). Thus, we document here a similar but less intense proteomic symbiont stress response to OW compared to former studies on A. gibbosa [49]. Concurrently, the only changes observed between pCO2 treatments at 28°C in diatom-associated proteins included slightly elevated chaperones, which are responsible for protein folding and repair, signifying a minor response of the photosymbionts to pCO2 levels approaching 3000 ppm (Tables 3 & S6). Minimal gene expression responses to OA were correspondingly found for the in hospite photosymbiotic dinoflagellates of corals [92,93], which suggest that host-derived tissue may act as a buffer, or lower susceptibility of endosymbiotic microalgae to OA. This is also in line with the significantly lower proportion of DAPs in the photosymbionts compared to the host found here amongst A. lobifera from the GoA. As responses suggesting cellular stress or acclimatization mechanisms are only detectable at the proteome level, and not reflected in bleaching or mortality data, our results confirm the comparably high stability of the photosymbiosis of A. lobifera originating from the GoA [22], and emphasize corresponding proteomic responses to ocean warming to lessen the physiological impact of ocean warming and acidification conditions expected for years 2100 and 2500 (RCP 8.5 [1]). The high resilience of A. lobifera is further evident from an exposure experiment showing that seven species of LBF (including three Amphistegina species) surviving a five-day incubation at temperatures up to 40°C and at pH between 5.9 and 7.4 [21].

Growth
While the LBF in our study did not exhibit significant mortality and bleaching, mean growth rates (2D planar extension) were generally lower under 31°C than under 28°C ( Figure 3). Such reduced growth is consistent with prior temperature-stress studies on LBF (e.g. [23,32,35,41]). This is confirmed by the significant decrease in host protein abundances at 31°C (Figures 5, 6b, S4-S6) in the present study. This caused a depletion of mitochondrial components (e.g., metabolic enzymes involved in the tricarboxylic acid cycle) and components of the electron transport chain (e.g., putative glyoxalase I, UDP-glucose 4-epimerase, isocitrate dehydrogenase [NADP], cytochrome c oxidase 1), but also of ribosomal and others proteins involved in translation (e.g., 60S and 40S ribosomal proteins, ribonuclease P protein subunit p25-like). Downregulation of reactive oxygen species (ROS)-producing reactions is a common proteomic response to reduce oxidative stress in marine organisms exposed to heat stress [94], while pathways that generate NADPH as a reducing equivalent for scavenging ROS may increase. This also emphasizes a potential decrease in energy production, which may affect resource allocation to the host or other cellular processes. Moreover, proteasome regulatory components (e.g., 26S proteasome non-ATPase regulatory subunit 1 homologue) decreased at 31°C, while proteasome subunits responsible for protein catabolic processes increased (Tables 4 & S6). Simultaneously, increases in abundance were detected for a few proteins involved in protein transport and folding, or stress response (e.g., endoplasmin, chaperone proteins), and var-ious proteins responsible for signal transduction, which may partly play roles in exocytosis, phagocytosis, or programmed cell death (Rab GTPase, Ras-related C3 botulinum toxin substrate 1, Ras-related protein RABB1c). All of these functional groups are likewise major traits of differential gene expression in the heat and bleaching responses of corals [95]. Hence, the dominant decreases in proteins from the 31°C treatments suggest cellular damage, reductions in energy production and cell-cycle related processes, and mild enhancement of unfolded protein response and repair processes.
The result of these proteomic changes suggests that fewer energy resources were available as a result of being reallocated to stress response mechanisms, potentially compromising growth. In comparison, one-month exposure of A. gibbosa to thermal stress (32°C) likewise lead to reduced growth, but also to putative inhibition of reproduction, and to visible bleaching along with impacts on photophysiological parameters [23]. Accordingly, the proteomic response of A. gibbosa was more prominent, e.g., by stronger increased abundance of protein-folding and stress-related proteins [49]. These disparities indicate a higher resistance of A. lobifera to ocean warming, which may be related to the differing photosymbiont responses. While A. gibbosa hosts only one lineage of photosymbiotic diatoms [35], A. lobifera in contrast hosts a much higher diversity of photosymbionts [22,52], such that the higher flexibility in the photosymbiotic association may facilitate increased overall resilience of the holobiont. This is supported by the proteomic data, which show that despite reductions in photosynthesis-related, metabolic, and biosynthetic proteins in the photosymbiont, overall bleaching and mortality was not induced by elevated temperature.
An effect of elevated pCO2 alone on LBF growth was not found in the present study, which is similar to growth results of the closely related LBF species A. gibbosa that was not impacted by pCO2 of up to ca. 2000 ppm [43], suggesting that energetic reserves are relatively stable under OA conditions. Conversely, Prazeres et al. [47] detected significantly reduced buoyant weight growth rates in A. lessonii in response to seawater pH levels of 7.7 and lower, while no significant impacts on growth rates expressed by changes in surface area were found for the same individuals. Fujita et al. [40] likewise found heavier tests at intermediate pCO2 (580 and 770 µatm), but lower test weight at 970 µatm in two other hyaline diatom-bearing foraminifera species. Recent work on the impacts of OA on coral calcification have documented a reduction in skeleton density, but not linear extension, in response to increasing pCO2 [96]. In our study, although not statistically significant, the lowest growth rates were found in the high-pCO231°C treatment (Figure 3a). A marginally significant increase in the frequency of partial breakage was also observed in this treatment (Figure 3d), which may have arisen from underlying changes in test structure and density associated with the reduced growth rates under these conditions. Future studies into the effects of ocean acidification on LBF should therefore use methods that assess changes in test density, weight, and calcification rate.
Another impact of OA documented in the present study is a significant decrease in test pore size under high-pCO231°C compared to control-pCO231°C (Figure 3b), which suggests that the structure of LBF tests may change under elevated pCO2, at least when combined with thermal stress. Such reduction in LBF pore size could limit the gas exchange between the environment and the LBF's inner cell. A 14 CO2-experiment demonstrated that the pores of amphisteginids allow the exchange of photosynthesis-and respiration-related gases, but not larger molecules [97]. Leutenegger & Hansen [97] also showed very slow CO2-uptake in A. lobifera compared to other species, and suggested that the organic pore elements may influence transport rates. Hence, the foraminiferal endoplasm and intercellular photosymbionts may confer some degree of protection from acidified environmental conditions.

Buffering of OA in the LBF microenvironment
Contrary to the majority of results documenting strong negative responses to CO2induced OA, such as dissolution of biogenic CaCO3 [98,99] and reducing overall fitness [43], the present study showed a generally mild response of LBF to OA. This is further documented through SEM images of test surfaces from living A. lobifera in treatments, which show no obvious signs of dissolution ( Figure 4). Some prior studies on the response of amphisteginid and other LBF tests to environmental stress [21,40,43] suggest that the foraminiferal organic matrix, its ectoplasm, and/or reticulopodial network help protect the test from dissolution [100]. These processes are likely assisted by increasing pH-as observed in the present study-and, thus, calcite saturation state of the DBL around LBF tests.
Our study notably shows that pH at the test surface (pHts) during dark is similar to surrounding seawater pH (pHSW), while under light conditions pHts at the LBF surface is elevated relative to pHSW (Figure 1b)-apparently due to the drawdown of DIC via symbiont photosynthesis under light conditions. Because DIC is a limiting substrate for photosynthesis in LBF [101], elevated pCO2 could therefore raise the photosynthetic output of the photosymbionts. Light-induced photosynthetic drawdown of DIC [67], as well as CO2uptake for fixation in the Calvin cycle in the dark, outweighing respired CO2 [12] and/or active regulation of CO2 in the foraminiferal microenvironment have been previously proposed for LBF [13,102], all of which would effectively minimize test dissolution by elevating pHts. Under light, the difference in proton concentration between seawater and test surface (∆[H + ]) became significantly larger at high-pCO231°C (Figure 1a), such that the resulting pHts (7.84 ± 0.21) was approaching pHts at moderate-pCO231°C (7.96 ± 0.04). Such elevated pH at the LBFs' test surface under light was similarly observed in A. radiata at 1141 and 2151 µatm [12], potentially indicating a shielding effect of maintaining a strong proton gradient in the DBL due to enhanced photosynthetic uptake of aqueous CO2 under conditions of elevated pCO2. Ecologically, photosynthesis and DBL thickness of LBF depend on flow rates [67], which they can control to some degree through their own motility, potentially allowing them to control DBL thickness and its protective function.
The mechanisms primarily considered responsible for the observed persistence of calcification in LBF under OA are active elevation of pH at the site of calcification via proton pumping (e.g., by V-type H + -ATPase, Ca 2+ -ATPase and Mg-ATPase [47,103]) and/or antiporter exchange, and/or drawdown of aqueous CO2 via photosynthesis [67], which compensate for the CO2-induced decrease in seawater pH and maintain pH homeostasis for cellular function [104]. Yet, proton gradients (∆[H + ]) in the LBF microenvironment weakened significantly under dark conditions (Figure 1a), with pH at the test surface becoming nearly indistinguishable from seawater pH) under all treatments ( Figure  1b)-suggesting that light is essential for maintaining the chemically protective pH gradients in the DBL of this species.
Notably, the A. lobifera maintained a much higher pHts and stronger ∆[H + ] in the highest pCO2 treatment under 31°C compared with 28°C, suggesting that this species is able to mitigate the impacts of extreme OA under the higher temperature treatment. This trend is not borne out in the growth data (as measured by change in cross-sectional area of test), which under the highest pCO2 treatment showed higher growth under 28°C versus 31°C. Interestingly, members of several electron transport complexes (e.g., cytochrome b-c1 complex subunit Rieske, cytochrome c oxidase subunit 5b-1, NADH-ubiquinone oxidoreductase, NADH dehydrogenase [ubiquinone] flavoprotein 1) increased with elevated pCO2, particularly at high-pCO228°C (H28), but not at 31°C (Figures S4-S6, Tables 4 & S6), suggesting that these may contribute to cellular homeostasis of this species under OA in the absence of thermal stress, but that these processes are not reflected in DBL chemistry.
Although the PERMANOVA indicated significant effects of pCO2 on global protein abundances in both compartments, the generally low number of DAPs and lack of proteomic data showing consistent regulation of common proton pumps (i.e. V-type proton ATPase subunits and vacuolar membrane proton pumps) under the elevated pCO2 treatments suggest that host cellular pH homeostasis was maintained without substantial regulatory response. Only a slight increase of these in control-and high-pCO231°C hint at active regulation of CO2 in the LBFs' microenvironment. In contrast, direct comparison of proteomic responses to high-and moderate-pCO228°C (Tables 4 & S6) suggests slightly reduced ion transport and protein phosphorylation (pyrophosphate-energized vacuolar membrane proton pump 1-like, V-type proton ATPase subunits, CAMK/CDPK protein kinase) under high compared to moderate pCO2. The few increased host-proteins under high pCO2 were involved in carbohydrate metabolic processes and proteolysis (e.g., UDPglucose 4-epimerase GalE, proteasome subunit alpha type-4-3). Similarly, a surprising lack of changes in abundance of ATP-driven proton pumps or even lower levels of proton pumps were evident in similarly designed coral OA experiments [105,106]. Nevertheless, though not statistically significant, ∆H + was also apparently enhanced at high-pCO231°C in the dark (Figure 1a), suggesting that other mechanisms facilitating proton regulation in the DBL may have been at work. These mechanisms could include higher activity of known ion pumps, or enhancement of so far unrecognized proton-regulating mechanisms and/or proteins.

Non-linear responses to multiple climate change stressors
The overall fitness in A. lobifera was largely unaffected by the climate change scenarios to which they were exposed in the present study. However, contrasting responses to the interaction of OA and OW were identified. Previous studies have highlighted taxonspecific and non-linear responses ( [8,11]) to OA, which were observed in the present study in the high amounts of DAPs, indicating cellular acclimatization mechanisms that facilitate the physiological functions of LBF, in response to moderate-pCO231°C, but not to high-pCO231°C,.
The proteomic response to moderate-pCO231°C is clearly anomalous in intensity and direction of variation (Figures 5c & 6), although this treatment did not cause detectable differences in other aspects of LBF physiology. The abundances of most DAPs decreased in the moderate-pCO231°C treatment, with only a small proportion increasing. The decreased proteins are involved in light harvesting, metabolic processes and protontransmembrane transport (e.g., fucoxanthin chlorophyll a/c proteins, divinyl chlorophyllide a 8-vinyl-reductase, thylakoid luminal 15 kDa protein 2), lipid trafficking (translocation protein sec 62), or cytoplasm organization (tubulin-specific chaperone, cell division protein FtsZ), biosynthetic processes and cell redox homeostasis (e.g., uroporphyrinogen decarboxylase, thioredoxin). Proteins involved in chlorophyll biosynthesis (e.g., geranylgeranyl reductase) increased, probably to offset the decrease in light-harvesting antenna proteins. Simultaneously, other proteins involved in biosynthetic and metabolic processes or cell redox homeostasis (e.g., aspartate-ammonia ligase, 40S ribosomal protein, dihydrolipoyl dehydrogenase), and protein folding and degradation (e.g., Clp-protease ATP binding subunit, ATP-dependent chaperone ClpB) increased. Overall, this hints at a stress-induced response of the photosymbionts under the moderate-pCO231°C treatment, probably related to oxidative stress. Although not detectable in chl a (Figure 2a), the photo-physiological condition of the symbionts may have been compromised. Hence, despite a presumably fertilizing effect of elevated DIC availability on photosynthesis, the ∆H + under light in this treatment stayed at the same level as under the control-pCO231°C treatment ( Figure 1). Thus, the proteomic response observed at moderate-pCO231°C was probably caused by oxidative and cell-redox stress offsetting the benefit of increased DIC available for symbiont photosynthesis under the year-2100 scenario.
Such cellular responses in the photosymbionts may affect holobiont performance by reducing transfer of resources. In the host proteome, in contrast, various proteins were enriched, and respiratory activity was highest in this treatment, as indicated by the lowered pHts in the dark. Some proteins involved in fatty-acid oxidation and mitochondrial electron transport (e.g., fatty-acid CoA synthase, phytanoyl-CoA dioxygenase) were increased, while other enzymes involved in fatty-acid and glucose metabolism (e.g., longchain-fatty-acid-CoA ligase, fructose-bisphosphate aldolase) were depleted. Altered fattyacid metabolism was previously identified as part of the core transcriptomic response of corals to OA [107]. The downregulation of fatty-acid catabolism and simultaneous upregulation of a fatty-acid synthesis indicate increased lipid storage, possibly a common acclimatization mechanism of photosymbiotic calcifiers to long-term acidification. However, several other metabolic proteins decreased, along with cell cycle-related processes (e.g., cullin-associated NEDD8-dissociated protein 1, eukaryotic initiation factor 4E), and signal transduction and regulation (e.g. GTPase HRas, Ras-related protein Rab-21), while proteins involved in translational, protein folding and degradation (e.g., 60S ribosomal protein L5, molecular chaperone Dank/HSP70, peroxiredoxin, autophagy-related protein 3), cytoskeleton or microtubule movement (dynein intermediate and heavy chain, clathrin light chain, villin-2), and protein phosphorylation or transport (e.g., putative vesicle-fusing ATPase, sorting nexin 1), increased. These processes likely reflect an oxidative stress response, including repair mechanisms and cytoskeleton restructuring, as observed previously [49,85].
LBF from the high-pCO231°C treatment exhibited an unexpectedly low number of DAPs and showed high similarity to LBF from the control-pCO231°C treatment ( Figures  5, 6, S2-S4, Table S7). Yet, pHts was significantly elevated to ca. 7.84 in the light (Figure 1) under high-pCO231°C. This indicates enhanced photosynthesis at ca. high-pCO231°C, although no significant change in chl a, color, or photosynthesis-related proteins were detected ( Figure 2, Tables 3 & S6). Thus, higher chloroplast efficiency under elevated pCO2 possibly enhanced photosynthetic productivity [12] without further proteomic adjustments. Decreases of a few proteins related to biosynthetic processes, oxidative stress, and cell redox homeostasis (e.g., acetyl-coA carboxylase, taurine catabolism dioxygenase TauD, glutamate synthase 2, L-ascorbate peroxidase,) and increases of DAPs, related to protein folding and proton transport (e.g., T-complex protein 1 subunit eta, Vtype proton ATPase proteolipid subunit), imply minor cellular effects that could indicate an acclimative mechanism that can potentially maintain cell redox homeostasis. This mild response underscores the high fitness of photosymbionts in the high-pCO231°C treatment, potentially facilitated by the fertilizing effect of OA on diatoms [8] and enhanced enzyme activities at high T. This further means the response of LBF photosymbionts to OA is non-linear when combined with warming. Hence, the impacts of 31°C and putative cell redox stress observed in moderate-pCO231°C are neutralized by the fertilizing effects of elevated pCO2 at a threshold between ca. 1000 and 3000 ppm. This seems to contradict the observation that DIC uptake of A. lobifera is close to saturation at ambient pCO2 [101], as higher temperatures facilitate higher enzymatic activities, which were presumably a limiting factor for CO2-fixation via symbiont photosynthesis in previous studies.
In the host proteome, only a very minor response was detectable at high-pCO231°C, therefore, no major proteomic responses were required, probably due to fully functional photosymbionts and/or the chemical protection conferred from the strong proton gradient maintained in their DBL, or limitations in gas exchange through reductions in pore size. This underscores how well-adjusted these prolific calcifiers are to low pH environments, as long as the activity of their photosymbionts is unaffected.

Conclusions
Climate change will have diverse and complex ecological consequences for marine ecosystems. Overall, the effects of two-months exposure to 31°C elicited a more adverse response than exposure to CO2-induced acidification at pCO2 conditions up to ca. 3000 ppm, which indicates the relatively high tolerance of A. lobifera to OA, potentially a result of a conserved adaptive trait in this and many other photosymbiotic Red Sea species [22,37,38,108,109]. However, both elevated pCO2 and temperature, but not their interaction, had significant effects on foraminifera and diatom proteomes. The frequency of these differentially abundant proteins was significantly higher in the foraminifera, indicating a stronger proteomic response in the host compared to its endosymbiotic diatoms. A mild heat-stress response was shown by the predominant decrease in protein abundances and reduced holobiont growth, indicating reduced holobiont performance through chronic sublethal cellular stress. OA showed even fewer effects on LBF. All pCO2 scenarios appear to be within the adaptive range of A. lobifera from the GoA. This sets their CO2 threshold above those previously proposed for Amphistegina spp. (e.g., for Great Barrier Reef populations: 1141 µatm pCO2 / pHNBS 7.85 [12]) and other hyaline foraminifera (e.g., from the northwest Pacific ~1000 µatm / pHSW 7.7 [40]). Yet at 31°C, elevated pCO2 induced cellular and morphological alterations (e.g., strongest proteome variations at moderate-pCO2 [year-2100]), indicating non-linear responses to combined stressors, and adaptive mechanism. Microenvironmental proton gradients [∆H + ] increasing with elevated pCO2 and light denote the organism's ability to mitigate the acidified conditions, especially during photosynthetic activity of the endosymbiotic diatoms. These responses may come at energetic costs that ultimately limit the host's resistance to climate change, as previous studies suggest [40,47], resulting in the reduced test growth and marginally increased frequency in test breakage under higher-pCO2 conditions observed in the present study. In general, the effects of OA (and its combination with other stressors) on LBF are highly variable amongst species and studies [19], possibly owing to differences in test structure (hyaline vs. porcellaneous), mineral polymorph solubility (aragonite vs. low-Mg calcite vs. high-Mg calcite [110], calcification mechanism [111], environmental preconditioning, compounding environmental factors (e.g., light, flow, nutrient availability, pesticides), life stage, flexibility in feeding strategies, and/or photosymbiont community.
Ancestral Amphisteginidae and other hyaline LBF, such as nummulitids, were diverse and dominant in shallow tropical marine carbonate environments during the middle Eocene [24,25], when atmospheric CO2 levels were similar to or markedly higher than today and underwent strong fluctuations [26,112]. The survival and, in some cases, proliferation of foraminifera throughout this interval invited the hypothesis that LBF possess important adaptations to changing pCO2. The ability of Amphistegina to continue calcifying under a wide range of pCO2 conditions, as observed in the present study, may have allowed this family of LBF to thrive during the middle Eocene and may facilitate their persistence in the future. The relative tolerance of these LBF to near-future climate change suggests a growing dominance in inter-reef and reef flat foraminifera assemblages [57] and further range expansions [22,113], along with a potential increased importance as carbonate producers. It also confirms the low susceptibility of coral reef organisms from the Gulf of Aqaba [22,37], which suggests that the northern Red Sea is a potentially important marine refuge from climate change [114], requiring strengthened protection and study.
The presented exploration of LBF holobiont response to the combined and singular effects of OA and OW at the level of the proteome yields valuable information about cellular stress responses, and the potential for physiological plasticity and resilience in the face of global climate change. Whether the putative acclimatization mechanisms can uphold overall fitness over longer time scales, and at what costs (e.g., energetic costs and/or physiological tradeoffs), should be assessed in future studies. With regard to the lack of signs of stress in the evaluated physiological parameters, but apparent proteomic responses, this study highlights the importance of considering impacts at several functional levels of the organism, as sublethal effects that can be (temporarily) overcome by acclimatization mechanisms at the molecular level, may not be detectable otherwise. Our experiment also shows that the effects of climate change can be missed when related stressors are assessed in isolation, and that physiological responses should be assessed across a range of organismal levels in order to make more realistic projections about the fate of reef calcifiers.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Figure S1: Phylogenetic assignment of photosymbionts; Figure S2: Heatmap of differentially abundant proteins per sample; Figure S3: Heatmap of differentially abundant proteins in photosymbionts; Figure  S4: Heatmap of differentially abundant proteins of the host foraminifera; Figure S5: Principal component analysis of differentially abundant proteins in (a) photosymbionts and (b) foraminiferal host including GO term summaries; Figure S6: Principal component analysis of differentially abundant proteins in (a) photosymbionts and (b) foraminiferal host including protein descriptions; Figure S7: Number of up-and down-regulated proteins comparing all treatments; Figure S8: Bar plots showing the number of gene ontology term hits of (a) pCO2-and (b) temperature-regulated proteins; Table  S1: Additional data for mean chemical conditions; Table S2: Results of nominal logistic models for the frequencies of test breakage and mortality; GuHCl, and 2 mM CaCl2 in 50 mM NH4HCO3 (pH 7.8) was added to each sample. After incubation at 37°C for 14 h, digested peptides were recovered by centrifugation followed by consecutive washing steps with 50 µl of 50 mM NH4HCO3 and 50 µl of ultra-pure water. Samples were acidified by addition of 20 µl of 10% TFA, desalted using SPEC C18, 4 mg sorbent (Agilent) as per manufacturer's instructions, and quality controlled as described previously [120].
A TIC based normalization was performed by measurement of aliquots of all samples on a nano LC-MS system to account for systematic errors. Samples (each 1 µg after normalization) were analysed using an Ultimate 3000 nano RSLC system coupled to a Q Exactive HF mass spectrometer (both Thermo Scientific). Preconcentration on a 100 µm × 2 cm C18 trapping column for 10 min using 0.1% TFA with a flow rate of 20 µl/min was followed by separation on a 75 µm x 50 cm C18 main column (both PepMap RSLC, Thermo Scientific) with a 120 min LC gradient ranging from 3-35% of buffer B: 84% ACN, 0.1% FA at a flow rate of 250 nl/min. The Q Exactive HF was operated in data-dependent acquisition mode and MS survey scans were acquired from m/z 300 to 1500 at a resolution of 60000 using the polysiloxane ion at m/z 371.101236 as lock mass [121]. Quadrupole isolation was set to 0.4 m/z and the fifteen most intense signals (Top15) were subjected to higher energy collisional dissociation with a normalized collision energy of 27% at a resolution of 15000 and a dynamic exclusion of 12 s. Automated gain control target values were set to 3 × 106 for MS and 5 × 104 MS/MS. Maximum injection times were 120 ms and 250 ms, respectively. Precursor ions with charge states of +1, >+5 or unassigned were excluded from MS/MS analysis.
Progenesis LC-MS software (version 3.0.6039.34628) from Nonlinear Dynamics (Newcastle upon Tyne, U.K.) and the statistical software R (version 3.3.162) [84] were used for label-free quantification. Raw MS data was imported into Progenesis and aligned to one automatically selected LC-MS reference file. After peak picking, MS/MS spectra were exported .mgf files. Identification of proteins and peptides was performed by the PEAKS Studio software suite, version 7.5, de novo to SPIDER [122]. Searches were performed in a decoy-fusion manner against a concatenated database comprising an updated and enlarged reference dataset for foraminifera, containing the recently published genome sequences of the calcareous benthic foraminifera Globobulimina [81], combined with the publicly available genome data of three diatom species (Bacillariophyta) from the NCBI database on 2016/03/15. Protein sequence databases were used as such, while nucleotide sequences were translated in all six reading frames prior to concatenation. Precursor mass tolerance was set to 10 ppm and fragment ion tolerance to 0.02 Da. Full tryptic enzyme specificity was selected, allowing a maximum of two missed cleavages. For de novo and DB search, carbamidomethylation of cysteines was defined as fixed modification and oxidation of methionine as variable modification. A maximum number of two variable modifications per peptide was allowed for de novo and DB search, while for the PTM and SPI-DER search the default settings were used. PSM-and peptide-level FDR was limited to 1% and proteins had to be identified with at least one unique peptide in order to be reported. Identifications meeting these criteria were re-imported into Progenesis to calculate the normalized abundances on the peptide-level.
Amino acid sequences of all proteins identified in the database search were uploaded in .fasta file format to the CD-HIT suite web server of the Weizhong Li Lab [123] in order to generate cluster of homologous protein sequences. Sequence identity cut-off was set to 0.7 and minimum alignment coverage for the longer sequence was set to 0.0. All other parameters were left as default. Normalized abundance-values of the peptides were used to calculate protein cluster abundance. Only peptides unique to a given protein cluster were used for quantification, hence, peptides that could not be clearly assigned to either of the compartments (host or symbionts) were removed from further analysis.