Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Seawater carbonate chemistry and biological processes during experiments with four species of marine diatoms

The effect of pCO2 on carbon acquisition and intracellular assimilation was investigated in the three bloom-forming diatom species, Eucampia zodiacus (Ehrenberg), Skeletonema costatum (Greville) Cleve, Thalassionema nitzschioides (Grunow) Mereschkowsky and the non-bloom-forming Thalassiosira pseudonana (Hust.) Hasle and Heimdal. In vivo activities of carbonic anhydrase (CA), photosynthetic O2 evolution, CO2 and HCO3? uptake rates were measured by membrane-inlet mass spectrometry (MIMS) in cells acclimated to pCO2 levels of 370 and 800 ?atm. To investigate whether the cells operate a C4-like pathway, activities of ribulose-1,5-bisphosphate carboxylase (RubisCO) and phosphoenolpyruvate carboxylase (PEPC) were measured at the mentioned pCO2 levels and a lower pCO2 level of 50 ?atm. In the bloom-forming species, extracellular CA activities strongly increased with decreasing CO2 supply while constantly low activities were obtained for T. pseudonana. Half-saturation concentrations (K1/2) for photosynthetic O2 evolution decreased with decreasing CO2 supply in the two bloom-forming species S. costatum and T. nitzschioides, but not in T. pseudonana and E. zodiacus. With the exception of S. costatum, maximum rates (Vmax) of photosynthesis remained constant in all investigated diatom species. Independent of the pCO2 level, PEPC activities were significantly lower than those for RubisCO, averaging generally less than 3%. All examined diatom species operate highly efficient CCMs under ambient and high pCO2, but differ strongly in the degree of regulation of individual components of the CCM such as Ci uptake kinetics and extracellular CA activities. The present data do not suggest C4 metabolism in the investigated species.

Seawater carbonate chemistry and sea urchin Paracentrotus lividus biological processes during experiments, 2011

Ocean acidification is predicted to have significant effects on benthic calcifying invertebrates, in particular on their early developmental stages. Echinoderm larvae could be particularly vulnerable to decreased pH, with major consequences for adult populations. The objective of this study was to understand how ocean acidification would affect the initial life stages of the sea urchin Paracentrotus lividus, a common species that is widely distributed in the Mediterranean Sea and the NE Atlantic. The effects of decreased pH (elevated PCO2) were investigated through physiological and molecular analyses on both embryonic and larval stages. Eggs and larvae were reared in Mediterranean seawater at six pH levels, i.e. pHT 8.1, 7.9, 7.7, 7.5, 7.25 and 7.0. Fertilization success, survival, growth and calcification rates were monitored over a 3 day period. The expression of genes coding for key proteins involved in development and biomineralization was also monitored. Paracentrotus lividus appears to be extremely resistant to low pH, with no effect on fertilization success or larval survival. Larval growth was slowed when exposed to low pH but with no direct impact on relative larval morphology or calcification down to pHT 7.25. Consequently, at a given time, larvae exposed to low pH were present at a normal but delayed larval stage. More surprisingly, candidate genes involved in development and biomineralization were upregulated by factors of up to 26 at low pH. Our results revealed plasticity at the gene expression level that allows a normal, but delayed, development under low pH conditions.

Treatment for Complicated Grief : State of the Science and Ways Forward

Peer reviewed

Design, Synthesis, and Biological Characterization of Largazole Analogues

Histone deacetylases (HDACs) have been shown to play key roles in tumorigenesis, and

have been validated as effective enzyme target for cancer treatment. Largazole, a marine natural

product isolated from the cyanobacterium Symploca, is an extremely potent HDAC inhibitor that

has been shown to possess high differential cytotoxicity towards cancer cells along with excellent

HDAC class-selectivity. However, improvements can be made in the isoform-selectivity and

pharmacokinetic properties of largazole.

In attempts to make these improvements and furnish a more efficient biochemical probe

as well as a potential therapeutic, several largazole analogues have been designed, synthesized,

and tested for their biological activity. Three different types of analogues were prepared. First,

different chemical functionalities were introduced at the C2 position to probe the class Iselectivity profile of largazole. Additionally, docking studies led to the design of a potential

HDAC8-selective analogue. Secondly, the thiol moiety in largazole was replaced with a wide

variety of othe zinc-binding group in order to probe the effect of Zn2+ affinity on HDAC

inhibition. Lastly, three disulfide analogues of largazole were prepared in order to utilize a

different prodrug strategy to modulate the pharmacokinetic properties of largazole.

Through these analogues it was shown that C2 position can be modified significantly

without a major loss in activity while also eliciting minimal changes in isoform-selectivity. While

the Zn2+-binding group plays a major role in HDAC inhibition, it was also shown that the thiol

can be replaced by other functionalities while still retaining inhibitory activity. Lastly, the use of

a disulfide prodrug strategy was shown to affect pharmacokinetic properties resulting in varying

functional responses in vitro and in vivo.

v

Largazole is already an impressive HDAC inhibitor that shows incredible promise.

However, in order to further develop this natural product into an anti-cancer therapeutic as well as

a chemical probe, improvements in the areas of pharmacokinetics as well as isoform-selectivity

are required. Through these studies we plan on building upon existing structure–activity

relationships to further our understanding of largazole’s mechanism of inhibition so that we may

improve these properties and ultimately develop largazole into an efficient HDAC inhibitor that

may be used as an anti-cancer therapeutic as well as a chemical probe for the studying of

biochemical systems.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

Seawater carbonate chemistry and biological processes during experiments with postlarvae of barnacle of Semibalanus balanoides and Elminius modestus, 2010

Ocean acidification and global warming are occurring concomitantly, yet few studies have investigated how organisms will respond to increases in both temperature and CO2. Intertidal microcosms were used to examine growth, shell mineralogy and survival of two intertidal barnacle post-larvae, Semibalanus balanoides and Elminius modestus, at two temperatures (14 and 19°C) and two CO2 concentrations (380 and 1,000 ppm), fed with a mixed diatom-flagellate diet at 15,000 cells ml-1 with flow rate of 10 ml-1 min-1. Control growth rates, using operculum diameter, were 14 ± 8 µm day-1 and 6 ± 2 µm day-1 for S. balanoides and E. modestus, respectively. Subtle, but significant decreases in E. modestus growth rate were observed in high CO2 but there were no impacts on shell calcium content and survival by either elevated temperature or CO2. S. balanoides exhibited no clear alterations in growth rate but did show a large reduction in shell calcium content and survival under elevated temperature and CO2. These results suggest that a decrease by 0.4 pH(NBS) units alone would not be sufficient to directly impact the survival of barnacles during the first month post-settlement. However, in conjunction with a 4-5°C increase in temperature, it appears that significant changes to the biology of these organisms will ensue.

Seawater carbonate chemistry, pigments and biological processes during experiments with coralline alga Corallina sessilis

Previous studies have shown that increasing atmospheric CO2 concentrations affect calcification in some planktonic and macroalgal calcifiers due to the changed carbonate chemistry of seawater. However, little is known regarding how calcifying algae respond to solar UV radiation (UVR, UVA+UVB, 280-400 nm). UVR may act synergistically, antagonistically or independently with ocean acidification (high CO2/low pH of seawater) to affect their calcification processes. We cultured the articulated coralline alga Corallina sessilis Yendo at 380 ppmv (low) and 1000 ppmv (high) CO2 levels while exposing the alga to solar radiation treatments with or without UVR. The presence of UVR inhibited the growth, photosynthetic O2evolution and calcification rates by13%, 6% and 3% in the low and by 47%, 20% and 8% in the high CO2 concentrations, respectively, reflecting a synergistic effect of CO2 enrichment with UVR. UVR induced significant decline of pH in the CO2-enriched cultures. The contents of key photosynthetic pigments, chlorophyll a and phycobiliproteins decreased, while UV-absorptivity increased under the highpCO2/low pH condition. Nevertheless, UV-induced inhibition of photosynthesis increased when the ratio of particulate inorganic carbon/particulate organic carbon decreased under the influence of CO2-acidified seawater, suggesting that the calcified layer played a UV-protective role. Both UVA and UVB negatively impacted photosynthesis and calcification, but the inhibition caused by UVB was about 2.5-2.6 times that caused by UVA. The results imply that coralline algae suffer from more damage caused by UVB as they calcify less and less with progressing ocean acidification.

Seawater carbonate chemistry and Gadus morhua length, weight and sperm biological processes, 2010

Ocean acidification, as a consequence of increasing marine pCO2, may have severe effects on the physiology of marine organisms. However, experimental studies remain scarce, in particular concerning fish. While adults will most likely remain relatively unaffected by changes in seawater pH, early life-history stages are potentially more sensitive - particularly the critical stage of fertilization, in which sperm motility plays a central role. In this study, the effects of ocean acidification (decrease of pHT to 7.55) on sperm motility of Baltic cod, Gadus morhua, were assessed. We found no significant effect of decreased pH on sperm speed, rate of change of direction or percent motility for the population of cod analyzed. We predict that future ocean acidification will probably not pose a problem for sperm behavior, and hence fertilization success, of Baltic cod.

Seawater carbonate chemistry and biological processes of Amphisorus hemprichii, Amphisorus kudakajimensis and Calcarina gaudichaudii during experiments, 2011

Ocean acidification, which like global warming is an outcome of anthropogenic CO2emissions, severely impacts marine calcifying organisms, especially those living in coral reef ecosystems. However, knowledge about the responses of reef calcifiers to ocean acidification is quite limited, although coral responses are known to be generally negative. In a culture experiment with two algal symbiont-bearing, reef-dwelling foraminifers, Amphisorus kudakajimensis and Calcarina gaudichaudii, in seawater under five different pCO2 conditions, 245, 375, 588, 763 and 907 µatm, maintained with a precise pCO2-controlling technique, net calcification of A. kudakajimensis was reduced under higher pCO2, whereas calcification of C. gaudichaudii generally increased with increased pCO2. In another culture experiment conducted in seawater in which bicarbonate ion concentrations were varied under a constant carbonate ion concentration, calcification was not significantly different between treatments in Amphisorus hemprichii, a species closely related to A. kudakajimensis, or in C. gaudichaudii. From these results, we concluded that carbonate ion and CO2 were the carbonate species that most affected growth ofAmphisorus and Calcarina, respectively. The opposite responses of these two foraminifer genera probably reflect different sensitivities to these carbonate species, which may be due to their different symbiotic algae.

Seawater carbonate chemistry and biological processes during experiments with oyster Crassostrea gigas, 2009

An increasing number of studies are now reporting the effects of ocean acidification on a broad range of marine species, processes and systems. Many of these are investigating the sensitive early life-history stages that several major reviews have highlighted as being potentially most susceptible to ocean acidification. Nonetheless there remain few investigations of the effects of ocean acidification on the very earliest, and critical, process of fertilization, and still fewer that have investigated levels of ocean acidification relevant for the coming century. Here we report the effects of near-future levels of ocean acidification (?0.35 pH unit change) on sperm swimming speed, sperm motility, and fertilization kinetics in a population of the Pacific oyster Crassostrea gigas from western Sweden. We found no significant effect of ocean acidification - a result that was well-supported by power analysis. Similar findings from Japan suggest that this may be a globally robust result, and we emphasise the need for experiments on multiple populations from throughout a species' range. We also discuss the importance of sound experimental design and power analysis in meaningful interpretation of non-significant results.

Seawater carbonate chemistry and biological processes of Emiliania huxleyi (strains RCC1256 and NZEH) during experiments, 2011

With respect to their sensitivity to ocean acidification, calcifiers such as the coccolithophore Emiliania huxleyi have received special attention, as the process of calcification seems to be particularly sensitive to changes in the marine carbonate system. For E. huxleyi, apparently conflicting results regarding its sensitivity to ocean acidification have been published (Iglesias-Rodriguez et al., 2008a; Riebesell et al., 2000). As possible causes for discrepancies, intra-specific variability and different effects of CO2 manipulation methods, i.e. the manipulation of total alkalinity (TA) or total dissolved inorganic carbon (DIC), have been discussed. While Langer et al. (2009) demonstrate a high degree of intra-specific variability between strains of E. huxleyi, the question whether different CO2 manipulation methods influence the cellular responses has not been resolved yet. In this study, closed TA as well as open and closed DIC manipulation methods were compared with respect to E. huxleyi's CO2-dependence in growth rate, POC- and PIC-production. The differences in the carbonate chemistry between TA and DIC manipulations were shown not to cause any differences in response patterns, while the latter differed between open and closed DIC manipulation. The two strains investigated showed different sensitivities to acidification of seawater, RCC1256 being more negatively affected in growth rates and PIC production than NZEH.

Seawater carbonate chemistry and biological processes of sand dollars Dendraster excentricus during experiments, 2011

Reduction in global ocean pH due to the uptake of increased atmospheric CO2 is expected to negatively affect calcifying organisms, including the planktonic larval stages of many marine invertebrates. Planktonic larvae play crucial roles in the benthic-pelagic life cycle of marine organisms by connecting and sustaining existing populations and colonizing new habitats. Calcified larvae are typically denser than seawater and rely on swimming to navigate vertically structured water columns. Larval sand dollars Dendraster excentricus have calcified skeletal rods supporting their bodies, and propel themselves with ciliated bands looped around projections called arms. Ciliated bands are also used in food capture, and filtration rate is correlated with band length. As a result, swimming and feeding performance are highly sensitive to morphological changes. When reared at an elevated PCO2 level (1000 ppm), larval sand dollars developed significantly narrower bodies at four and six-arm stages. Morphological changes also varied between four observed maternal lineages, suggesting within-population variation in sensitivity to changes in PCO2 level. Despite these morphological changes, PCO2 concentration alone had no significant effect on swimming speeds. However, acidified larvae had significantly smaller larval stomachs and bodies, suggesting reduced feeding performance. Adjustments to larval morphologies in response to ocean acidification may prioritize swimming over feeding, implying that negative consequences of ocean acidification are carried over to later developmental stages.

Seawater carbonate chemistry and biological processed during experiments with corals Porites sp. & Stylophora pistillata, 2010

Uptake of anthropogenic CO2 by the oceans is altering seawater chemistry with potentially serious consequences for coral reef ecosystems due to the reduction of seawater pH and aragonite saturation state (omega arag). The objectives of this long-term study were to investigate the viability of two ecologically important reef-building coral species, massive Porites sp. and Stylophora pistilata, exposed to high pCO2(or low pH) conditions and to observe possible changes in physiologically related parameters as well as skeletal isotopic composition. Fragments of Porites sp. and S. pistilata were kept for 6-14 months under controlled aquarium conditions characterized by normal and elevated pCO2 conditions, corresponding to pHTvalues of 8.09, 7.49, and 7.19, respectively. In contrast with shorter, and therefore more transient experiments, the long experimental timescale achieved in this study ensures complete equilibration and steady state with the experimental environment and guarantees that the data provide insights into viable and stably growing corals. During the experiments, all coral fragments survived and added new skeleton, even at seawater omega arag <1, implying that the coral skeleton is formed by mechanisms under strong biological control. Measurements of boron (B), carbon (C) and oxygen (O) isotopic composition of skeleton, C isotopic composition of coral tissue and symbiont zooxanthellae, along with physiological data (such as skeletal growth, tissue biomass, zooxanthellae cell density and chlorophyll concentration) allow for a direct comparison with corals living under normal conditions and sampled simultaneously. Skeletal growth and zooxanthellae density were found to decrease, whereas coral tissue biomass (measured as protein concentration) and zooxanthellae chlorophyll concentrations increased under high pCO2 (low pH) conditions. Both species showed similar trends of delta11B depletion and delta18O enrichment under reduced pH, whereas the delta13C results imply species-specific metabolic response to high pCO2 conditions. The skeletal delta11B values plot above seawater delta11B vs. pH borate fractionation curves calculated using either the theoretically derived deltaB value of 1.0194 (Kakihana et al., Bull. Chem. Soc. Jpn. 50(1977), 158) or the empirical deltaB value of 1.0272 (Klochko et al., EPSL 248 (2006), 261). However, the effective deltaB must be greater than 1.0200 in order to yield calculated coral skeletal delta11B values for pH conditions where omega arag >1. The delta11B vs. pH offset from the literature seawater delta11B vs. pH fractionation curves suggests a change in the ratio of skeletal material laid down during dark and light calcification and/or an internal pH regulation, presumably controlled by ion-transport enzymes. Finally, seawater pH significantly influences skeletal delta13C and delta18O. This must be taken into consideration when reconstructing paleo-environmental conditions from coral skeleton