Printer Friendly

Effects of seawater acidification on gene expression: resolving broader-scale trends in sea urchins.


Since the onset of the industrial revolution, human activities have been rapidly adding carbon dioxide (C[O.sub.2]) to Earth's atmosphere (Caldeira and Wickett, 2003), and atmospheric C[O.sub.2] levels are now higher than at any point in the past 800,000 years (Luthi et al., 2008). In addition to causing the global warming responsible for rises in atmospheric and ocean temperature, anthropogenic C[O.sub.2] is also responsible for shifts in the chemical composition of the global ocean. Approximately one-quarter to one-half of atmospheric C[O.sub.2] is dissolved into our world's oceans (Sabine et al., 2004; Boyd, 2011), where it reacts with seawater to increase the concentration of hydrogen ion ([H.sup.+]) and decrease the concentration of carbonate ion (C[O.sub.3.sup.2-]) (Caldeira and Wickett, 2005). This process, termed ocean acidification, has already reduced global ocean pH by 0.1 units since pre-industrial times (IPCC, 2007), and is predicted to decrease seawater pH by 0.45 units over the next century and by 0.77 units over the next 300 years. This anticipated rate of change is faster than any experienced in the last 300 million years (Caldeira and Wickett, 2003, 2005; Honisch et al., 2012).

In the context of such rapid and widespread change, understanding the consequences of ocean acidification has emerged as a major priority in marine science research (Raven et al., 2005; National Research Council, 2010). Spanning sub-disciplines from oceanography to ecology, investigators are using a variety of scientific tools to forecast the effects of continued ocean acidification on a range of marine organisms. Genomic approaches have become an informative way to investigate responses to ocean acidification and have become more prevalent in recent years. The advent of RNA sequencing now allows genomic-based inquiry in marine organisms with little or no previous sequence information (e.g., De Wit and Palumbi, 2013) and genome sequencing of ecologically and economically important marine organisms impacted by declines in ocean pH, such as the purple sea urchin (Strongylocentrotus purpuratus Stimpson, 1857) (Pearse, 2006; Sodergren et al., 2006; Anderson et al., 2011), Pacific oyster (Crassostrea gigas) (Zhang et al., 2012), and coccolithophore Emiliana huxleyii (Benner et al., 2013) have further encouraged the use of genomics in ocean acidification research. A common application of nucleic acid sequence information involves tracking changes in gene expression in organisms exposed to near-future pC[O.sub.2] scenarios. The philosophy behind this approach is simple: marine organisms exposed to high pC[O.sub.2] seawater will induce compensatory changes in physiology that are facilitated in part by changes in gene expression (Gracey, 2003; Gracey and Cossins, 2003), and knowledge of the numbers, magnitude of change, and function of genes regulated by seawater pC[O.sub.2] informs questions surrounding the ability of these organisms to cope with ocean change (Evans and Hofmann, 2012). For example, functional analyses of genes regulated by seawater pH assists in identifying biological processes most strongly influenced by ocean acidification and provides mechanistic understanding of differential sensitivity to seawater pC[O.sub.2] among marine species or populations (Evans et al., 2013; Granados-Cifuentes et al., 2013). Analyses of gene expression may also highlight potential evolutionary responses to shifts in seawater chemistry, as the genes that contribute to short-term plastic responses to ocean acidification also represent probable targets of natural selection driving the evolution of tolerance toward more acidic oceans (Whitehead and Crawford, 2006; Crawford and Oleksiak, 2007; Whitehead, 2012).

Efforts to characterize ocean-acidification-induced changes in gene expression have been implemented using both small- and large-scale approaches across a variety of marine taxa. From candidate gene approaches that target a small number of genes with well-described functions, to high-throughput approaches that track expression across the entire transcriptome or proteome, expression-based ocean acidification studies have been performed in corals (e.g., Moya et al., 2012; Vidal-Dupiol et al., 2013), crabs (e.g., Fehsenfeld et al., 2011), coccolithophores (e.g., Jones et al., 2013), oysters (e.g., Tomanek et al., 2011; Dineshram et al., 2012; Liu et al., 2012), mussels (e.g., Huning et al., 2013), barnacles (e.g., Wong et al., 2011), and other important marine species (e.g., Hauton et al., 2009; Crawfurd et al., 2011). Sea urchins have emerged as an ideal system for monitoring changes in gene expression caused by ocean acidification. Urchin embryos, larvae, and adults all construct calcium carbonate skeletons and are therefore predicted to be susceptible to the decreases in carbonate ion concentration that accompany ocean acidification (Dupont et al., 2010; Hofmann et al., 2010; Kroeker et al., 2010; Byrne et al., 2011, 2013). The importance of the purple sea urchin as a laboratory model used around the world expedited the sequencing of its genome and stimulated the use of urchins in genetic experiments that study the effects of ocean acidification. This combination of ecological relevance and experimental power has led to the publication of eight manuscripts that analyzed gene expression during exposure to acidified seawater in four species of urchins (Strongylocentrotus purpuratus, Lytechnius pictus, Hemicentrotus pulcherrimus, and Paracentrotus lividus). While these studies are undoubtedly informative, gene expression data are sensitive to differences in experimental methodology, which prevents the genetic response to ocean acidification from being directly compared even between studies using the same species. For example, pronounced differences in treatment regimes between studies almost certainly influence the number and type of genes deemed responsive to seawater pC[O.sub.2]. Martin et al. (2011) detected increases in the expression of three key genes involved in calcification in Mediterranean urchin (P. lividus) larvae cultured in low pH seawater, a result that was not replicated in studies of purple sea urchins (S. purpuratus) that used lower pC[O.sub.2] treatment levels (Todgham and Hofmann, 2009; Evans et al., 2013). Discrepancies such as these ultimately hinder the ability of researchers to resolve common or divergent responses to ocean acidification between species or populations. Other experimental parameters, such as method of detecting differential expression (e.g., quantitative polymerase chain reaction (qPCR) vs. microarray vs. RNA sequencing), statistical approaches (supervised approaches including Student's t-tests and analyses of variance vs. unsupervised approaches such as principal components analysis), and life stage (early stage vs. adult) also influence gene expression, with these differences again limiting broader-scale conclusions about the effects of ocean acidification on gene expression.

One way to overcome methodological differences between studies and reveal common patterns across datasets is to perform a meta-analysis. Meta-analyses have proven effective in finding predominant trends within ocean acidification data, despite individual studies reporting conflicting results (Ries et al., 2009: Hendriks et al., 2010; Kroeker et al., 2010, 2013). Although a sufficient mass of gene expression data now exists to enable comparative analyses of the genetic response to ocean acidification in sea urchins, we are unaware of any efforts to do so. In this study, we perform meta-analyses of genes deemed significantly differentially expressed in sea urchins exposed to simulated ocean acidification as a means of identifying major outcomes across these studies. We use the compiled gene expression data to address five prevailing hypotheses regarding the effects of ocean acidification on sea urchins and marine organisms in general: exposure to high pC[O.sub.2] seawater (i) depresses metabolism, (ii) impairs calcification, (iii) disrupts pH homeostasis, (iv) delays development; and finally, (v) future pC[O.sub.2] regimes will exceed organism tolerance limits. Our results largely support these five hypotheses and reassert the vulnerability of sea urchins in future oceans (Wittmann and Portner 2013). We also aim to identify a core response to seawater acidification in urchins that could be used to increase our understanding of species-specific responses to ocean acidification and develop more accurate biomarkers of exposure to low pH seawater. We resolved a set of 26 genes common to multiple studies that further emphasizes the impact of ocean acidification on metabolism, calcification, and pH regulation.

Materials and Methods

Data selection

We searched the biological literature for studies that tracked changes in gene expression in sea urchins exposed to acidified seawater. Literature searches were conducted using the IS I Web of Science database and different combinations of relevant keywords that included ocean acidification, seawater acidification, hypercapnia, sea urchin, urchin, echinoderm, gene expression, gene, genomics, transcriptomics, and proteomics. Literature searches were performed until 21 September 2013. We retained studies regardless of the method used to measure gene expression, the pC[O.sub.2]/pH exposure level used, or the life stage or species investigated. Database searches using these criteria returned eight manuscripts that were used in subsequent analyses (see Table 1).

Compilation of differentially expressed gene lists

Each of the eight manuscripts were parsed for genes that the authors considered significantly differentially expressed in response to seawater acidification. Differentially expressed genes from Evans et al. (2013) and Padilla-Gamino et al. (2013) were extracted from lists contained in the online supplementary information accompanying those manuscripts. Differentially expressed genes from the remaining six manuscripts were extracted from figures, tables, or other sources present in each manuscript. Using this gene expression data, we complied four gene sets that were used in downstream analyses: (i) the complete list of all differentially expressed genes, (ii) genes down-regulated by high pC[O.sub.2] seawater, (iii) genes up-regulated by high pC[O.sub.2] seawater, and (iv) genes differentially expressed in two or more studies.

To generate the list of genes differentially expressed in at least two studies, we first assigned all differentially expressed genes to GLEAN gene models (e.g., SPU_X-XXXXX) that were used to annotate the S. purpuratus genome (as described in Oliver et al., 2010; e-scores <0.001). This conversion step allowed us to use a common set of identifiers throughout our analyses and was necessary because the software used to determine intersection of gene sets requires an exact character match in the chosen identifier (Oliveros, 2007). Gene names could not be used to determine intersection because annotations of the same gene

often use slightly different nomenclature (e.g., sodium/potassium transporting ATPase alpha subunit v.v. [Na.sup.+]/[K.sup.+]-ATPase subunit alpha). GLEAN gene model identiners were listed in Todgham and Hofmann (2009), Hammond and Hoffman (2012), Evans et al. (2013), and Padilla-Gamino et al. (2013). GLEAN identifiers were also listed in the transcriptomic study of Lytechnius pictus, as a microarray containing S. purpuratus sequences was used (O'Donnell et al., 2010). For the remaining three studies, GenBank sequence identifiers listed in Kurihara et al. (2012), Martin et al. (2011), and Stumpp et al. (2011a) were used in a Basic Local Alignment Search Tool (BLAST) search to annotate genes differentially expressed in Hemicentrotus pulcherrimus, Paracentrotus lividus, or .S', purpuratus with GLEAN identifiers, respectively. We used the online program Venny (Oliveros, 2007) to identify overlapping GLEAN identifiers across the eight studies. Only genes with identical GLEAN gene model identifiers were considered the same gene.

Gene functional annotation and gene set enrichment analysis

Gene functional annotation and gene set enrichment analyses (GSEA) were performed on each of the four lists to resolve the biological significance of genes regulated by seawater pC[O.sub.2] in sea urchins. Gene functional annotation serves to organize lists of genes into functionally related groups called ontologies, which assists in defining the range of molecular, cellular and physiological processes regulated by seawater acidification. To do so, each differentially expressed gene was assigned a UniProt identifier using GLEAN models to retrieve corresponding RNA sequences from the S. purpuratus transcriptome, ver. 3.0 (Cameron et ai, 2009). These sequences were then used in a BLAST search of the Uniprot database (UniProt Consortium, 2014) to retrieve corresponding UniProt identifiers for each gene ( Finally UniProt identifiers were used to assign each gene to one or more ontologies using the PANTHER Database, ver. 8.0 (Thomas et al., 2003, 2006). Three different ontology databases within PANTHER were used: (i) Gene Ontology: Biological Process, (ii) Gene Ontology: Molecular Function, and (iii) PANTHER Protein Class. Integrating information across these three databases provided a method to validate results between databases and more accurately identify processes strongly affected by seawater pC[O.sub.2].

UniProt identifiers and the PANTHER Database were also used to perform GSEA (Thomas et al., 2003; Thomas et al., 2006; Mi and Thomas, 2009). GSEA is a statistical tool that uses gene functional information to identify categories of genes significantly over- or under-represented within a user-defined list. Significance using the PANTHER database is determined by the binomial statistic; that is, the probability that the number of genes observed from a given ontology occurred by chance relative to the total number of genes from this ontology present in the entire S. purpuratus transcriptome. The same three PANTHER ontology databases were again used to summarize results. Following PANTHER recommendations, we considered ontologies with Bonferroni false-discovery rate corrected P-values < 0.05 to be significant (Thomas et al., 2006). A separate GSEA was performed on gene sets that excluded five genes identified as differentially expressed using qPCR (but not also differentially expressed in one of the four microarray studies), as these genes were not selected randomly and may bias results. Results of the GSEA excluding the qPCR genes are available in the supplementary material (http://www. for comparison. However, interpretation of results focuses on analyses performed using all genes.

An alternative approach was taken to investigate the impact of seawater pC[O.sub.2] on the expression of biomineralization genes. Genes involved in urchin biomineralization are very well-characterized within the scientific literature (Oliveri et al., 2002; Livingston et al., 2006), but annotate poorly in functional databases like PANTHER due to the absence of orthologous proteins in organisms used to populate the databases. As a result, functional analyses using these databases are not appropriate for determining whether a taxa-specific process such as biomineralization is strongly affected by seawater acidification. To circumvent this issue, we used an urchin-specific list of 345 biomineralization genes used previously by Pespeni et al. (2013a) and searched for overlap between this list and the list of differentially expressed, down-regulated, up-regulated, and overlapping genes defined in this study. A Fisher's exact test was then used to determine whether the proportion of genes involved in biomineralization from each list was higher than expected given the total number of biomineralization genes encoded in the S. purpuratus genome.


All genes regulated by seawater acidification in urchins

Our database search yielded eight manuscripts reporting expression changes in response to elevated seawater pC[O.sub.2] among four different species of sea urchin (Table 1). A total of 602 genes were considered significantly differentially expressed across these eight studies (Fig. SI) (http://www. The vast majority of differentially expressed genes originated from the four studies that used high-throughput microarray-based approaches (n = 594), while the remaining eight were obtained from studies using candidate gene approaches and qPCR. During exposure to acidified seawater, 483 genes were down-regulated (80%) while only 134 (22%) were up-regulated. Interestingly, 15 differentially expressed genes exhibited conflicting directions of change between different studies (Fig. S1). Forty-three of the 602 genes differentially expressed during ocean acidification in urchins were also found within the urchin-specific list of biomineralization genes, including 10 spicule matrix proteins and three isoforms of carbonic anhydrase (Table 2).

We next used gene functional annotation to determine the complete range of biological processes regulated by seawater acidification in sea urchins. We were able to functionally annotate 591 of the 602 differentially expressed genes to ontologies associated with 17 major biological processes (with most genes annotating to more than one functional category) (Fig. 1). While genes with a variety of functions were differentially expressed, genes involved in metabolism were by far the most numerous, with 251 genes (42% of total) annotating to the "metabolic process" ontology within the PANTHER Biological Process database. Given the very large number of genes contributing to this ontology, we also examined the number of genes annotating to finer-scale ontologies within the "metabolic process" parent ontology (Fig. 2A). Genes involved protein (n = 104: e.g., dipeptidase 1, 26S protease regulatory subunit 6A, matrix metalloproteinase-17); nucleic acid (n = 91: e.g., ribonuclease T2, DNA helicase INO80); carbohydrate (n = 44: e.g., lactate dehydrogenase, pyruvate dehydrogenase, citrate synthase, succinate dehydrogenase, isocitrate dehydrogenase); lipid (n = 30; e.g., elongation of very long chain fatty acids protein 6, long-chain-fatty-acid-CoA ligase 1); and amino acid (n = 25: e.g., L-asparaginase, cysteine dioxygenase type 1, cysteine/glutamate transporter) metabolism were all differentially expressed in large numbers. The next most frequent ontology was "cellular process" containing 149 genes (25% of total) (Fig. 2B). The cellular process ontology includes genes involved in cell communication (n = 102: e.g., MAP kinase-activated protein kinase 5); cell adhesion (n = 39: e.g., integrin alpha-4, protocaherin-12); cell motion (n = 29: e.g., alpha and beta tubulin); and the cell cycle (n = 47: e.g., cyclin-dependent kinase 1). The ontology containing the third largest number of genes was "transport" containing 112 genes (19% of total) involved in the transport of ions (n = 54: e.g., sodium/potassium-transporting ATPase subunit alpha); amino acids (n = 11 : e.g., proton-coupled amino acid transporter 1); lipids (n = 6: e.g., apolipoprotein D); nucleotides (n = 6: e.g., ADP/ ATP translocase 3); vitamins (n = 9; e.g., sodium-dependent vitamin C transporter); and other molecules (Fig 2C).

While functional annotation is an appropriate method to identify the full range of biological processes affected by seawater acidification, these analyses do not consider the relative proportions of genes changing within specific ontologies. For example, the S. purpuratus genome encodes for many more genes annotating to metabolic ontologies than to other ontologies; consequently, it may be expected that a large number of metabolic genes will respond to seawater acidification simply due to their abundance in the genome. To resolve this issue and identify processes most strongly impacted by seawater acidification in urchins, we performed GSEA, which identifies ontologies where the number of contributing genes are in significantly higher proportion relative to their total number in the S. purpuratus genome. Using GSEA we identified two biological processes most strongly influenced by pC[O.sub.2] in sea urchins: ion transport and metabolism (Table 3). Ontologies associated with ion transport were the most numerous and had the lowest P-values in the GSEA. Significantly over-represented ontologies relating to ion transport included "ion transport," "cation transport," and "transport" from the Biological Processes database; "transporter activity," "transmembrane transporter activity," and "cation transmembrane transporter activity" from the Molecular Function database; and "transporter" and "cation transporter" from the Protein Class database. A large number of ontologies clearly linked to metabolic changes were also significantly over-represented (Table 3). These data indicate that the metabolic genes present in the functional annotation were not simply the result of a larger number of metabolic genes in the S. purpuratus genome. Significantly over-represented metabolic ontologies from the Biological Process database included "generation of precursor metabolites and energy," "carbohydrate metabolic process," "oxidative phosphorylation," "respiratory electron transport chain," and "tricarboxylic acid cycle." A significantly higher proportion of genes involved in biomineralization were differentially expressed in urchins exposed to high pC[O.sub.2] seawater than would be expected based upon the number of biomineralization genes encoded in the S. purpuratus genome (Fisher's exact test P < 0.001; Fig. S2; mental/).

Genes down-regulated by seawater acidification in urchins

We performed similar analyses on the set of genes down-regulated by pC[O.sub.2] to identify biological processes repressed by seawater acidification. Since input gene lists for down-regulated genes and all differentially expressed genes were largely the same (483 of 602 genes, or 80% similarity), functional annotation of down-regulated genes mirrored that of the complete list, with genes annotating most frequently to "metabolic process" (n = 194), "cellular process" (n = 114), and "transport" (n = 100) ontologies (Fig. 1). The same trend was apparent in the GSEA performed on down-regulated genes, with the majority of significantly overrepresented ontologies functioning in ion transport and metabolism (Figs. S3, S4; supplemental/). Thirty-eight down-regulated genes are involved in urchin biomineralization, including seven spicule matrix proteins and two carbonic anhydrases (Fig. S2). This is again a significantly higher proportion than would be expected based upon the number of biomineralization genes encoded in the S. purpuratus genome (Fisher's exact test P < 0.001).

Genes up-regulated by seawater acidification in urchins

We next performed analyses on the 134 up-regulated genes as a means to resolve processes activated in response to seawater acidification in urchins. These genes functionally annotated to 16 different ontologies with the largest number of up-regulated genes annotating to "metabolic process" (n = 65), followed by "cellular process" (n = 39), "cellular communication" (n = 23), and "transport" (n = 22) (Fig. 1). We also performed GSEA on the set of upregulated genes to resolve biological processes activated in response to seawater acidification. Far fewer ontologies were significantly over-represented among genes up-regulated as compared to genes down-regulated by seawater acidification, a function of the much lower number of input genes (134 up-regulated genes compared to 483 down-regulated genes). However, four ontologies were significantly over-represented among the set of up-regulated genes: "cellular calcium ion homeostasis" and "homeostatic process" from the Biological Process database, "cation transmembrane transporter activity" from the Molecular Function database, and "cation transporter" from the Protein Class database (Fig. S3). Despite the induction of many genes involved in metabolism, ontologies with apparent metabolic functions were not significantly over-represented among up-regulated genes. Nine of the 134 genes up-regulated during ocean acidification in urchins were involved in biomineralization (Fig. S2), which again was a significantly higher proportion than expected by chance (Fisher's exact test P <0.001).

Genes regulated by seawater acidification in more than one study

A total of 26 genes were common to at least two studies, five genes were differentially expressed in at least three studies, and one gene, [Na.sup.+], [K.sup.+]-ATPase alpha subunit, was differentially expressed in four studies. No genes were found to intersect in more than four studies (Fig. SI). Strikingly, 25 of the 26 overlapping genes were down-regulated (94%), while only spicule matrix protein SM30 was induced by seawater acidification in multiple studies. Genes annotating to metabolic and ion transport ontologies were most frequent among the 26 overlapping genes: 10 genes annotated to "transport," including the [Na.sup.+], [K.sup.+]-ATPase alpha subunit and H(+)-transporting ATPase beta subunit; and seven genes annotated to "metabolic process," including succinyl-CoA synthetase alpha subunit and cytochrome c oxidase subunit IV precursor (Fig. 1). GSEA revealed that ontologies associated with ion transport and metabolism were significantly over-represented, and several ontologies also over-represented in the complete list of differentially expressed genes were present, including "oxidative phosphorylation," "generation of precursor metabolites and energy," "respiratory electron transport chain," "cation transport," and "cation transmembrane transporter activity" (Table 4). "Cellular calcium ion homeostasis" was the only ontology significantly over-represented in both the set of overlapping genes and the set of up-regulated genes. We did not perform separate GSEA on subsets of up- or down-regulated genes from within the 26 overlapping genes given the low number of up-regulated genes in this set (n = 1) and the high degree of similarity between the set of down-regulated genes and all overlapping genes (25 of 26 genes). Seven of the 26 genes (39%) common to at least two studies intersected with the list of urchin biomineralization genes, including H(+-transporting ATPase beta subunit and msp130 protein (i.e., spicule matrix protein 130) (Fig. S2). Once again, this represents a significantly higher proportion than would be expected by chance (Fisher's exact test P < 0.001).

Genes not strongly affected by seawater acidification

We were also interested in determining the effect of seawater acidification on genes involved in development and the cell stress response. Exposure to elevated pC[O.sub.2] negatively affects development in marine embryos and larvae, including echinoderms (Dupont and Thomdyke, 2009; Dupont et al., 2010). We hypothesized that impaired development caused by exposure to seawater acidification in urchins would alter the expression of developmentally regulated genes, and in turn, these genes would be detected as differentially expressed when compared to expression levels in urchins raised under control conditions. Seventy-nine genes regulated by seawater acidification in urchins functionally annotated to the "developmental process" ontology (Fig. 1). However, GSEA showed that ontologies relating to development were not significantly over-represented among any of the four gene sets analyzed (Fig. S3).

As organisms approach environmental tolerance limits, a conserved set of genes functioning to repair and prevent cellular damage are induced (Kultz, 2003, 2005; Petrak et ai, 2008; Wang et al., 2009a; Evans and Hofmann, 2012). Little evidence for cellular damage or the induction of a cellular stress response was apparent in the urchin expression data. Most notably, 80% of genes affected by seawater acidification were down-regulated, a pattern that does not support the induction of genes involved in the cellular stress response. The set of up-regulated genes was not significantly over-represented for ontologies with apparent functions in the cellular stress response, and there was very limited evidence that cellular function was severely compromised. Four genes indicative of macromolecular damage were induced: heat shock protein 40, DNAJb11 (molecular chaperones involved in re-folding damaged proteins; Feder and Hofmann, 1999), ubiquitin (which targets irreparable proteins for destruction; Dutton and Hofmann 2008; Hofmann and Somero 1995), and caspase activation and recruitment domain 15 (part of the signaling cascade that ultimately causes apoptosis; Festjens et al., 2006) (Fig. SI).


Inherent problems with comparing changes in gene expression caused by ocean acidification

As ecologically and economically important marine calcifiers, sea urchins have emerged as ideal models for studying the effects of ocean acidification. Many researchers have exploited nucleic acid sequence resources provided by the Strongylocentrotus purpuratus genome to undertake studies that monitor shifts in gene expression triggered by exposure to acidified seawater. However, much like studies of ocean acidification in general, differences in experimental methods, target species or populations, and statistical approaches make direct comparison of results problematic. Comparisons of gene expression data between studies are a necessity to properly relate new results to that of previous work, but most of these comparisons are made in the context of major caveats that reduce their impact. For example, four studies analyzing changes in gene expression caused by seawater acidification in urchins used qPCR to evaluate shifts in mRNA abundance for a small number of genes (Table 1). Enhanced sensitivity of qPCR compared to alternative approaches like microarrays could lead to genes being reported as differentially expressed in one study, but unchanged in another (Wang et al., 2006). The trade-off associated with qPCR-based inquiry is that while sensitive, expression can only be tracked across a small number of pre-selected genes, which prevents discovery of novel genes with potentially important roles in the response to seawater acidification. Four other studies used higher-throughput transcriptomic approaches to examine genes regulated by seawater pC[O.sub.2] in urchins (Table 1). These studies eliminate biases associated with gene selection and are capable of monitoring expression across the transcriptome, but are themselves associated with methodological discrepancies that make comparing gene lists challenging. For example, O'Donnell et al. (2010) used a 5. purpuratus microarray platform to track acidification-induced changes in gene expression in the painted urchin, Lytechinus pictus. This approach, termed heterologous hybridization, limits the number of genes in which expression data can be acquired because L. pictus mRNAs will not bind strongly to divergent S. purpuratus sequences spotted on the array (Renn et al., 2004; Buckley, 2007). Statistical approaches to analyzing high-throughput and candidate gene expression also differ dramatically. Whereas a small number of genes can be assessed using standard statistical approaches, analyzing high-throughput expression data requires more specialized techniques. In most cases, arbitrary filtering criteria are applied to high-throughput datasets to reduce the number of genes used in downstream statistical tests (Tseng et al., 2001; Hackstadt and Hess, 2009), with the consequence of ignoring genes with low levels of expression but potentially important roles in the response to ocean acidification. For example, both microarray and RNA-sequencing datasets are often filtered by removing genes whose fluorescence intensity or sequence count, respectively, is below a certain threshold in all samples. Consequently, genes are often eliminated from analyses due to low expression in a single replicate. Importantly, low expression in one sample could reflect minor technical discrepancies in sample preparation or processing just as easily as it could a biological response or lack thereof. The emergence of RNA-sequencing technology provides enhanced sensitivity compared to other transcriptomic techniques such as microarrays (Wang, Z., et al., 2009), but seems to result in the identification of a far greater number of significantly differentially expressed genes compared to microarray experiments (Sirbu et al., 2012). This difference in sensitivity is a potentially important issue considering that the number of differentially expressed genes is often used as an indicator of how strongly an organism is affected by environmental change.

Meta-analyses have been used successfully to resolve broader-scale trends within the sometimes contradictory scientific literature of ocean acidification (Ries et al., 2009; Kroeker et al., 2010, 2013; Hendriks et al., 2010), and here we performed a meta-analysis of eight studies monitoring gene expression in sea urchins in order to resolve key questions surrounding the effects of ocean acidification. A central goal was to determine whether gene expression data collected to date support five prevailing hypotheses about the effects of ocean acidification on echinoderms. Gene expression profiles in urchins strongly support that exposure to high pC[O.sub.2] seawater will depress metabolism, impact biomineralization, and disrupt pH homeostasis. Expression patterns also provide inconclusive evidence that ocean acidification slows developmental rate, and little evidence that short-term exposure to future ocean pH regimes will exceed organism tolerance limits.

The pervasive effects of ocean acidification

An important conclusion from this study is that the impact of ocean acidification on urchins extends beyond biomineralization and pH regulation. Functional annotation of 591 genes regulated by seawater pC[O.sub.2] indicates that ocean acidification is a pervasive environmental variable that affects a wide range of biological phenomenon not typically associated with it, from reproduction to cell adhesion to the immune response (Fig. 1). Mechanistic bases for susceptibility or resilience toward future acidification could lie within these largely overlooked processes, and urchin expression data underscore the need to evaluate a wider range of biological processes in future studies. Trends in the scientific literature imply that the ocean acidification research community has already reached a similar conclusion and is actively exploring the consequences of seawater pC[O.sub.2] on processes outside calcification and pH regulation. For example, a link between ocean acidification and immune function has been demonstrated repeatedly (Hernroth et al., 2011, 2012; Matozzo et al., 2012), including a recent study in the green sea urchin, S. droebachiensis, that suggests a direct link between extracellular pH and phagocyte numbers (Dupont and Thomdyke, 2012). Byssal threads used to anchor mussels to rocks within the intertidal are weakened and compromise attachment to substratum when secreted in more acidic seawater (O'Donnell et al., 2013), data that provide at least some indication as to why the expression of genes encoding adhesive proteins may be modified by ocean acidification in sea urchins. However, mechanistic understanding of these newly discovered responses is sparse, and questions surrounding the nature of pC[O.sub.2]-regulated changes in, for example, cell communication, are at this point speculative. These unanswered questions illustrate that searching for impacts beyond calcification and pH regulation will continue to be an important step in understanding the consequences ocean acidification.

Biological processes most strongly affected by ocean acidification

The expression of a disproportionate number of genes involved in metabolism, ion transport, and biomineralization were modified by seawater acidification in sea urchins, providing clear evidence that these three processes will be most strongly affected by changes in ocean pC[O.sub.2]. This conclusion is congruent with well-established outcomes regarding the physiological impact of ocean acidification on marine organisms--specifically, that exposure to high pC[O.sub.2] seawater depresses metabolism (Pdrtner, 2008), triggers compensatory movements of ions across cell membranes (Stumpp et al., 2012), and influences biomineralization (Kroeker et al., 2010; Byrne et al., 2013). These conclusions validate gene expression as an accurate method to assess responses to shifting environmental conditions and endorse the use of expression-based approaches in future ocean acidification research.

Metabolic depression. Eighty percent of the genes differentially expressed in sea urchins exposed to acidified seawater were down-regulated, suggesting that most biological processes were repressed. By far the largest number of down-regulated genes annotated to the "metabolic process" ontology, including genes participating in the metabolism of carbohydrates, lipids, proteins, and nucleic acids (Figs. 1, 2A). These data strongly suggest that future ocean acidification will have an inhibitory effect on the metabolism of sea urchins. Results of the GSEA suggest that elevated seawater pC[O.sub.2] triggers a decrease in the production of ATP by down-regulating the expression of genes involved in the tricarboxylic acid cycle, electron transport chain, and oxidative phosphorylation. Significantly over-represented ontologies among the set of down-regulated genes included "respiratory electron transport chain," "tricarboxylic acid cycle," and "oxidative phosphorylation" (Fig. S3; http:// The importance of modulating these ATP-producing pathways is highlighted by the differential expression of succinyl CoA synthetase, which catalyzes the conversion of succinyl CoA to succinate as part of the tricarboxylic acid cycle (Li et al., 2013); and cytochrome c oxidase subunit IV, the last enzyme in the electron transport chain (Arnold, 2012). These two genes were among only 26 others significantly differentially expressed in two studies. Furthermore, these metabolic pathways are closely connected, as the tricarboxylic acid cycle supplies precursor metabolites to the electron transport chain, and the electron transport chain is a fundamental component of ATP production via oxidative phosphorylation (Leninger, 1999).

ATP production is tightly coupled to organism energy demands, thus a decline in ATP production is an indication that urchins either exhibit a reduced physiological demand for ATP in high pC[O.sub.2] seawater or strategically suppress metabolism as a means to extend the duration of tolerance until more favorable conditions return (Todgham and Hofmann, 2009). We hypothesize that the down-regulation of genes involved in ATP production reflects metabolic depression rather than declining ATP requirements in more acidic seawater. Adaptive responses to the environment are energetically expensive endeavors that require large amounts of ATP to fuel processes such as protein synthesis (Kiiltz, 2003, 2005). The ATP-dependent processes involved in acclimating to changes in seawater pC[O.sub.2] are unlikely to be mediated by down-regulating genes involved in the cell's primary ATP production pathways (Le., oxidative phosphorylation). Alternatively, metabolic depression occurs through the coordinated suppression of ATP-generating pathways (Guppy and Withers, 1999) and has been observed in diverse organisms exposed to acidified seawater, including corals (Nakamura et al., 2011), sipunculids (Reipschlager and Portner, 1996), mussels (Michaelidis et ai, 2005), crabs (Carter et al., 2013), and scallops (Schalkhausser et al., 2012). Metabolic depression is considered an adaptive strategy that allows organisms to temporarily extend the duration of tolerance by inhibiting energetically expensive processes (Fabry et al., 2008). In purple sea urchins, the species of choice for five of the eight studies analyzed here, temporary metabolic depression may promote survival during episodic upwelling events that expose nearshore environments in the northeast Pacific to transient declines in seawater pH that can reach as low as 7.4 (Evans et al., 2013). However, metabolic depression provides only short-term relief from the negative effects of seawater acidification, and the use of this strategy to cope with permanently high pC[O.sub.2] seawater is unlikely. Vital metabolic functions cannot be repressed over the long-term without negative consequences for growth and development, and the obvious trend toward metabolic depression in urchins raises uncertainty as to whether these organisms have the capacity to acclimatize to oceans permanently high in pC[O.sub.2]- The potential vulnerability of urchins to conditions expected in future oceans is supported by their reduced abundance in the low pH marine ecosystems surrounding volcanic C[O.sub.2] vents (Hall-Spencer et al., 2008; Kroeker et al., 2011) and by the fact that urchins are considered to be among the poorest acid-base regulators (Calosi et al., 2013). However, caution should be exercised in extrapolating trends observed in short-term experiments, like those evaluated here, to chronically high pC[O.sub.2] oceans. The difficulty in implementing long-term ocean acidification experiments has resulted in a paucity of data on the capacity of marine populations to acclimate to long-term increases in seawater acidity. The potential for misstep in using brief exposures to infer long-term survival in future oceans is illustrated in a study of the presumably vulnerable cold-water coral Lophelia pertusa, which exhibited a decline in calcification during a one-week exposure to low pH seawater, but was actually able to increase calcification rates over the course of a 6-month exposure (Form and Riebesell, 2012). The physiological response to seawater acidification also differed temporally over the course of a 2-month exposure in the mussel Mytilus edulis (Thomsen and Melzner, 2010). More long-term studies will be necessary to determine whether urchins can reverse metabolic depression and acclimatize to permanently high pC[O.sub.2] conditions. Discretion is also needed when using contemporary genotypes to assess fitness in future oceans. Accumulating evidence suggests that purple sea urchins will be capable of rapidly adapting to more acidic oceans (Kelly et al., 2013; Pespeni et al., 2013b). In this case evolution, rather than phenotypic plasticity, may allow urchin populations to persist in a more acidic ocean.

Ion transport. Ion transport was the second most strongly affected biological process in sea urchins exposed to low pH seawater. One quarter of the genes that could be functionally annotated were associated with ion transport (Figs. 1, 2C), and both the up-regulated and down-regulated gene sets were significantly over-represented for ontologies relating to the movement of ions across cell membranes (Table 3, Fig. S3). Differential expression of ion transporters is likely a function of pH regulation. Although urchins are considered poor acid-base regulators (e.g., Calosi et al., 2013), accumulating evidence suggests that they possess some capacity to compensate for decreases in extracellular pH. Larval green sea urchins (congeners of S. purpuratus) compensate for acid-base imbalance with a bicarbonate buffering system that uses energy derived from transmembrane sodium gradients to export hydrogen ions (Stumpp et al., 2012). A total of 15 genes involved in the transport of sodium were differentially expressed in response to high pC[O.sub.2] seawater in urchins, including a [K.sup.+]-dependent [Na.sup.+], [Ca.sup.+2] exchanger, subunits of the [Na.sup.+]/[K.sup.+]-transporting ATPase, and two [Na.sup.+] channels (Fig. S1; http://www.biolbull. org/supplemental/). An important trend regarding the expression of transport proteins is that despite the very strong signature of metabolic depression, urchins up-regulated several ATP-dependent ion transporters. Both the alpha-3 and beta-1 subunits of the sodium/potassium-transporting ATPase and the H(+)-transporting ATPase beta subunit increased in expression during exposure to acidified seawater. In fact, the alpha subunit of the sodium/potassium-transporting ATPase was the only gene deemed significantly differentially expressed in four studies (Fig. SI). The transmembrane movement of ions is metabolically expensive and can consume up to 77% of larval metabolism in sea urchins (Leong and Manahan, 1997). Induction of energy-demanding ion transporters within an overall framework of metabolic depression underscores the importance of ion transport in maintaining homeostasis in high pC[O.sub.2] seawater. Although sodium gradients may fuel acid-base regulation, the actual transmembrane movement of bicarbonate and hydrogen ions that will ultimately help restore pH homeostasis primarily depends on two genes: the [Na.sup.+]/[H.sup.+]and [Cl.sup.-]/HC[O.sup.-.sub.3] exchangers (Wheatly and Henry, 1992). The expressions of these transport proteins did not change in any of the eight studies analyzed here, indicating that although modifying the abundance of sodium transporters is potentially important to pH regulation, the actual movement of protons and bicarbonate required to restore acid-base homeostasis may rely solely on changes in the activity of these transporters rather than on increasing abundance.

Biomineralization. Biomineralization is highly sensitive to the changes in seawater chemistry expected with continued ocean acidification (Kroeker et al., 2010), and increases in seawater pC[O.sub.2] have been shown to reduce calcification rates in a wide range of calcifying organisms, including urchins (Shirayama and Thornton, 2005; Ries et al., 2009; Byrne et al., 2011; Yu et al., 2011, 2013; Courtney et al., 2013). Larval urchins secrete fragile but functionally important skeletons, and consequently, the impact of ocean acidification on skeletogenesis in early-life-stage urchins has received considerable attention. Indeed, all eight of the studies analyzed here chose to investigate early life stages (Table 1). Stunted growth and impaired biomineralization appears to be a nearly universal response to future levels of seawater acidification in early-life-stage urchins: significant reductions in the length of the larval arms or supporting skeletal rods have been documented across 13 species from tropical, temperate, and polar habitats (Byrne et al., 2013). An advantage to using gene expression in assessing biomineralization in urchins is that gene networks underlying skeletogenesis are extremely well characterized (Oliveri et al., 2002; Livingston et al., 2006). The ubiquitous effect of pC[O.sub.2] on biomineralization in larval urchins is consistent with patterns of gene expression resolved in this study. The complete set of differentially expressed genes (43 biomineralization genes of 602 differentially expressed), the set of up-regulated (9 biomineralization genes of 134 up-regulated genes), the set of down-regulated (38 biomineralization genes of 483 down-regulated genes), and the set of overlapping genes (7 biomineralization genes of 26 overlapping genes) all contained significantly higher proportions of genes involved in biomineralization relative to the total number of biomineralization genes encoded by the S. purpuratus genome (Fig. S2; mental/). Thirty-eight of the 43 biomineralization genes differentially expressed in response to acidified seawater in urchins were down-regulated, a result that supports skeletogenesis being impaired.

Mechanistically, reduced biomineralization may be linked to insufficient abundance of spicule matrix proteins, critical components of the urchin endoskeleton that participate exclusively in calcification (Wilt, 2002) and function to bind and sequester [Ca.sup.+2] ions before deposition into the endoskeleton (Farach-Carson et al., 1989). Seven of the 10 differentially expressed spicule matrix proteins were downregulated, including members of the spicule matrix protein 16, 29, 30, 50, and 130 gene families (Table 2). We hypothesize that impaired calcification, down-regulation of spicule matrix proteins, and decreased expression of other genes involved in biomineralization are occurring as part of an overall metabolic depression in urchins. This hypothesis is congruent with the predominant effect of pC[O.sub.2] on the expression of metabolic genes, and several studies have suggested that impaired urchin skeletogenesis in acidified seawater is primarily an energetic problem rather than a function of reduced carbonate saturation in seawater. Intracellular pH regulation is an energy-demanding process and therefore consumes resources that would normally be directed toward skeletogenesis, with negative consequences for skeletal development (Stumpp et al., 201 la; Byrne et al., 2013).

A reduced ability to calcify in more acidic seawater reinforces the vulnerability of urchins in future oceans (Wittmann and Portner, 2013). Larval skeletons are functionally important for feeding and swimming, and while the degree to which biomineralization is affected between species and habitat types varies (Byrne et al., 2013), even seemingly small impacts on skeletogenesis can reduce the fitness of urchin embryos and larvae. Food capture in larval urchins depends upon the surface area of ciliary bands that surround the larval arms. A 7% reduction in band length translates to a 16% reduction in feeding rates in the echinoid sand dollar Dendraster excentricus (Strathmann, 1971). The ability of urchin larvae to disperse could also be affected by small changes in skeletal morphology that negatively affect swimming performance (Pennington and Strathmann, 1990; Strathmann and Grunbuam, 2006). Reduced skeletogenesis may also influence ecosystem-level processes, as the inability of larval urchins to construct their skeleton may leave them at a greater risk for predation (Allen, 2008) and may even disturb energy transfer across trophic levels if decreased calcification translates into an overall decrease in larval size (Sheridan and Bickford, 2011).

The set of 125 up-regulated genes was not significantly over-represented for metabolic ontologies, further supporting the inhibitory effect of increasing seawater pC[O.sub.2] on urchin metabolism. Alternatively, genes up-regulated in response to acidified seawater were significantly over-represented for ontologies relating to calcium homeostasis, including several calcium transporters (e.g., similar to voltage-dependent T-type calcium channel) and binding proteins (e.g., calmodulin-like protein 5). Shifts in the expression of genes that modify the availability of calcium may represent an adaptive response to promote calcification in a carbonate-limited environment (Evans et al., 2013). Seawater is the ultimate source of calcium incorporated into urchin skeletons, and large amounts of calcium must first be imported via calcium transporters before being deposited into the skeleton (Wilt, 2002). Evidence for this potential function is strengthened by similar transcriptomic responses to high pC[O.sub.2] seawater in the corals Acropora millepora and Pocillopora damicornis (Moya et al., 2012; Vidal-Dupiol et al., 2013). Both studies reported increases in the abundance of calcium transporters, and both predicted that these expression changes were related to calcification, either directly via an increase in the amount available for calcification (Vidal-Dupiol et al., 2013) or indirectly through the alteration of calcium-dependent signaling pathways (Moya et al., 2012). Although the connection between calcium transport and calcification is well-supported, the up-regulation of calcium transporters and binding proteins in the context of so many other biomineralization genes being down-regulated raises several questions. For example, why would urchins choose to up-regulate genes involved in calcium homeostasis rather than spicule matrix proteins or other genes involved in biomineralization? Also why, despite changes in the expression of calcium-related genes, is skeletogenesis still impaired in the majority of urchins tested? (Byrne et al., 2013). The purpose of inducing genes involved in calcium homeostasis is also confounding considering the clear evidence for metabolic depression in urchins exposed to high pC[O.sub.2] seawater and that the transport of calcium, like the transport of sodium and other ions, is ATP-dependent. Uptake of calcium from seawater is temporally variable during early urchin development (Tellis et al., 2013), and differences in the life-stage sampled between each study could plausibly lead to the differential expression of genes controlling calcium import. An alternative speculation is that fluxes in calcium are related to intracellular signaling events of which calcium is a major player (Moya et al., 2012). Given that calcium transport appears to be regulated by seawater pC[O.sub.2] in both urchins and corals, experiments addressing these uncertainties may be worth pursuing.

Processes not strongly affected by ocean acidification

Our analyses of gene expression were informative both in terms of identifying biological processes most commonly regulated by seawater pC[O.sub.2] in sea urchins, and also in defining processes not strongly affected by seawater acidification. Functional annotation showed that ocean acidification is a pervasive stressor that affects a range of biological processes; however, genes involved in two processes commonly associated with the response to ocean acidification and environmental perturbation in general did not appear to be widely affected: development and the cellular stress response.

Development. Numerous studies have examined the developmental consequences of ocean acidification, with developmental delay reported to occur in a wide range of species including crustaceans, molluscs, and echinoderms (Dupont and Thorndyke, 2009). While development in high pC[O.sub.2] seawater has been reported to cause developmental delay in urchins (Martin et al., 2011 ; Stumpp et al., 201 la, b), this response is not unequivocal (Byrne et al., 2013), and other studies have failed to detect an effect on embryonic (cell cycle progression, first mitotic division, onset of DNA synthesis, and mitotic spindle formation; Place and Smith, 2012) or larval development (prism and early pluteus stages; Padilla-Gamino et al., 2013). The genetic underpinnings of sea urchin development have been extensively studied and involve the coordinated expression of genes within well-described developmental programs (Wei et al., 2006). Should ocean acidification slow development, the expression of developmentally regulated genes should also be delayed, causing these genes to be detected as significantly differentially expressed relative to those of urchins developing at normal rates in lower pC[O.sub.2] seawater. We observed limited evidence for this response in the gene expression patterns analyzed here. Ontologies relating to development were not significantly over-represented among genes differentially expressed, suggesting that development is not among the processes strongly affected by seawater acidification. Seventy-nine genes regulated by seawater acidification in urchins annotated to the "developmental process" ontology, a result that does not refute a pC[O.sub.2]-mediated developmental delay in urchins. However, these shifts in gene expression could also be driven by compensatory responses within specific developmental pathways rather than by an overall developmental delay. Todgham and Hofmann (2009) suggest that the differential expression of key transcription factors involved in early developmental programs reflects an effect on gene networks regulating skeletogenesis rather than an overall developmental delay, and a similar mechanism was proposed by Evans et al (2013) as a means for promoting calcification in carbonate-limited environments. Disentangling an overall developmental delay from effects on specific developmental pathways is likely not feasible using gene expression exclusively. More directed research will be required to determine which of these two possibilities, if either, is most prevalent in urchins.

Cellular stress response. Environmental stress levels that approach organism tolerance limits induce a specific set of proteins that repair and prevent further macromolecular damage; these include chaperones that re-fold denatured proteins, members of protein degradation pathways, and cell cycle regulators that arrest cell division and initiate programmed cell death (Kultz, 2003, 2005; Petrak et al., 2008; Wang, P., et al., 2009). While we acknowledge that gene expression is an indirect measure of short-term organism tolerance limits, genes comprising the cellular stress response are nonetheless often used for exactly this purpose (?e.g., Lockwood and Somero, 2010; Tomanek, 2012). Expression profiles in urchins suggest that short-term exposure to near-future pC[O.sub.2] regimes do not approach survivable limits and that the negative impacts of ocean acidification on urchins will more likely manifest through reductions in fitness (e.g., impaired calcification, fertilization, metabolic deficits). This conclusion is largely consistent with recent meta-analyses on the effects of near-future ocean acidification on urchins (Dupont et al., 2010; Byrne et al., 2013). However, experiments conducted in urchins to date present only a snapshot of life in future oceans and provide little information as to how these animals will cope with chronic exposure to high pC[O.sub.2] seawater. Some urchin species are extremely long-lived and may experience a century worth of ocean acidification in a single lifetime (Ebert and Southon, 2003).


The authors thank Dr. Melissa Pespeni (University of Vermont) for sharing the list of urchin-specific biomineralization genes. This research was supported by start-up funds provided to TGE by California State University East Bay.

Literature Cited

Allen, J. D. 2008. Size-specific predation on marine invertebrate larvae. Biol. Bull. 214: 42-49.

Anderson, S. C., J. Milis Fleming, R. Watson, and H. K. Lotze. 2011. Rapid global expansion of invertebrate fisheries: trends, drivers and ecosystem effects. PLoS One 6(3): e14735.

Arnold, S. 2012. The power of life--Cytochrome c oxidase takes center stage in metabolic control, cell signalling and survival. Mitochondrion 12: 46-56.

Benner, L, R. E. Diner, S. C. Lefebvre, D. Li, T. Komada, E. J. Carpenter, and J. H. Stillman. 2013. Emiliana huxleyi increases calcification but not expression of calcification-related genes in long-term exposure to elevated temperature and pC[O.sub.2]. Philos. Trans. R. Soc. B 368: doi:org/10.1098/rstb.2013.0049.

Boyd, P. W. 2011. Beyond ocean acidification. Nat. Geosci. 4: 273-274.

Buckley, B. A. 2007. Comparative environmental genomics in nonmodel species: using heterologous hybridization to DNA-based microarrays. J. Exp. Biol. 210: 1602-1606.

Byrne, M., M. Ho, E. Wong, M. A Soars, P. Selvakumaraswamy, H. Shepard-Brennard, S. A. Dworjanyn, and A. R. Davis. 2011. Unshelled abalone and corrupted urchins: development of marine calcifiers in a changing ocean. Proc. R. Soc. B. 278: 2376-2383.

Byrne, M., M. Lamare, D. Winter, S. A. Dworjanyan, and S. Uthicke. 2013. The stunting effect of a high C[O.sub.2] ocean on calcification and development in sea urchin larvae, a synthesis from the tropics to the poles. Proc. R. Soc. B. 368: doi:org/10.1098/rstb.2012.0439.

Caldeira, K., and M. E. Wickett. 2003. Oceanography: anthropogenic carbon and ocean pH. Nature 425: 365.

Caldeira, K., and M. E. Wickett. 2005. Ocean model predictions of chemistry changes from carbon dioxide emissions to the atmosphere and ocean. J. Geophys. Res. 110: 1-12.

Calosi, P., S. P. S. Rastrick, M. Graziano, S. C. Thomas, C. Baggini, H. A. Cater, J. M. Hall-Spencer, M. Milazzo, and J. Spicer. 2013. Distribution of sea urchins living near shallow water C[O.sub.2] vents is dependent upon species acid-base and ion-regulatory abilities. Mar. Pollut. Bull. 73: 470-484.

Carter, H. A., L. Ceballos-Osuna, N. A. Miller, and J. A. Stillman. 2013. Impact of ocean acidification on metabolism and energetics during early life stages of the intertidal porcelain crab Petrolisthes cinctipes. J. Exp. Biol. 216: 1412-1422.

Courtney, T., I. Westfield, and J. B. Ries. 2013. C[O.sub.2]-induced ocean acidification impairs calcification in the tropical urchin Echinometra viridis. J. Exp. Mar. Biol. Ecol. 440: 169-175.

Crawford, D. L., and M. F. Oleksiak. 2007. The biological importance of measuring individual variation. J. Exp. Biol. 210: 1613-1621.

Crawfurd, K. J., J. A. Raven, G. L. Wheeler, E. J. Baxter, and I. Joint. 2011. The response of Thalassiosira pseudonana to long-term exposure to increased C[O.sub.2] and decreased pH. PLoS One 6: e26695.

De Wit, P., and S. R. Palumbi. 2013. Transcriptome-wide polymorphisms of red abalone (Haliotis rufescens) reveal patterns of gene flow and local adaptation. Mol. Ecol. 22: 2884-2897.

Dineshram, R., K. K. Wong, S. Xiao, Z. Yu, P. Y. Qian, and V. Thiyagarajan. 2012. Analysis of Pacific oyster larval proteome and its response to high-C[O.sub.2]. Mar. Pollut. Bull. 64: 2160-2167.

Dupont, S., and M. C. Thorndyke. 2009. Impact of C[O.sub.2]-driven ocean acidification on invertebrates early life-history--What we know, what we need to know and what we can do. Biogeosciences 6: 3109-3131.

Dupont, S., and M. Thorndyke. 2012. Relationship between C[O.sub.2]driven changes in extracellular acid-base balance and cellular immune response in two polar echinoderm species. J. Exp. Mar. Biol. Ecol. 424: 32-37.

Dupont, S., O. Ortega-Martinez, and M. Thorndyke. 2010. Impact of near future acidification on echinoderms. Ecotoxicology 19: 449-462.

Dutton, J. M., and G. E. Hofmann. 2008. Spatial and temporal variation in distribution and protein ubiquitination for Mytilus congeners in the California hybrid zone. Mar. Biol. 154: 1067-1075.

Ebert, T. A., and J. R. Southon. 2003. Red urchins (Strongylocentrotus francisanus) can live over 100 years: confirmation with A-bomb [sup.14]carbon. Fish. Bull. 101: 915-922.

Evans, T. G., and G. E. Hofmann. 2012. Defining the limits of physiological plasticity: how gene expression can assess and predict the consequences of ocean change. Philos. Trans. R. Soc. B 367: 1733-1745.

Evans, T. G., F. Chan, B. A. Menge, and G. E. Hofmann. 2013. Transcriptomic responses to ocean acidification in larval sea urchins from a naturally low pH environment. Mol. Ecol. 22: 1609-1625.

Fabry, V. J., B. A. Seibel, R. A. Feely, and J. C. Orr. 2008. Impacts of ocean acidification on marine fauna and ecosystem processes. ICES J. Mar. Sci. 65: 414-432.

Farach-Carson, M. C., D. D Carson, J. L. Collier, W. J. Lennarz, H. R. Park, and G. C. Wright. 1989. A calcium-binding, asparagine-linked oligosaccharide is involved in skeleton formation in the sea urchin embryo. J. Cell Biol. 109: 1289-1299.

Feder, M. E., and G. E. Hofmann. 1999. Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu. Rev. Physiol. 61: 243-282.

Fehsenfeld, S., R. Kiko, Y. Appelhaus, D. W. Towle, M. Zimmer, and F. Melzner. 2011. Effects of elevated seawater pC[O.sub.2] on gene expression patterns in the gills of the green crab, Carcinus maenas. BMC Genomics 12: 488.

Festjens, N., S. Cornelis, M. Lamkanfi, and P. Vandenabeele. 2006. Caspase-containing complexes in the regulation of cell death and inflammation. Biol. Chem. 387: 1005-1016.

Form, A. U., and U. Riebesell. 2012. Acclimation to ocean acidification during long-term C[O.sub.2] exposure in the cold-water coral Lophelia pertusa. Glob. Change Biol. 18: 843-853.

Gracey, A. Y. 2003. Interpreting physiological responses to environmental change through gene expression profiling. J. Exp. Biol. 210: 1584-1592.

Gracey, A. Y., and A. R. Cossins. 2003. Application of microarray technology in environmental and comparative physiology. Annu. Rev. Physiol. 65: 231-259.

Granados-Cifuentes, C., A. J. Bellantuono, T. Ridgway, O. Hoegh-Guldberg, and M. Rodriguez-Lanetty. 2013. High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity. BMC Genomics 14: 228.

Guppy, M., and P. Withers, 1999. Metabolic depression in animals: physiological perspectives and biochemical generalizations. Biol. Rev. 74: 1-40.

Hackstadt, A. J., and A. M. Hess. 2009. Filtering for increased power for microarray data analysis. BMC Bioinform. 10: 11.

Hall-Spencer, J. M., R. Rodolfo-Metalpa, S. Martin, E. Ransome, M. Fine, S. M. Turner, S. J. Rowley, D. Tedesco, and M. C. Buia. 2008.

Volcanic carbon dioxide vents show ecosystem effects of ocean acidification. Nature 454: 96-99.

Hammond, L. T. M., and G. E. Hofmann. 2012. Early developmental gene regulation in Strongylocentrotus purpuratus embryos in response to elevated C[O.sub.2] seawater conditions. J. Exp. Biol. 215: 2445-2454.

Hauton, C., T. Tyrrell, and J. Williams. 2009. The subtle effects of sea water acidification on the amphipod Gammarus locusta. Biogeosciences 6: 1479-1489.

Hendriks, L E., C. M. Duarte, and M. Alvarez. 2010. Vulnerability of marine biodiversity to ocean acidification: a meta-analysis. Estuar. Coast. Shelf Sci. 86: 157-164.

Hernroth, B., S. Baden, M. Thorndyke, and S. Dupont. 2011. Immune suppression of the echinoderm Asterias rubens (L.) following long-term ocean acidification. Aquat. Toxicol. 103: 222-224.

Hernroth, B., H. N. Skold, K. Wiklander, F. Jutfelt, and S. Baden. 2012. Simulated climate change causes immune suppression and protein damage in the crustacean Nephrops norvegicus. Fish Shellfish Immunol. 33: 1095-1101.

Hofmann, G. E., and G. N. Somero. 1995. Evidence for protein damage at environmental temperatures-seasonal cycles in levels of ubiquitin conjugates and Hsp70 in the intertidal mussel Mytilus trossulus. J. Exp. Biol. 198: 1509-1518.

Hofmann, G. E., J. P. Barry, P. J. Edmunds, R. D. Gates, II. A. Hutchins, T. Klinger, and M. A. Sewell. 2010. The effect of ocean acidification on calcifying organisms in marine ecosystems: an organism-to-ecosystem perspective. Annu. Rev. Ecol. Evol. Syst. 41: 127-147.

Honisch, B., A. Ridgwell, D. N. Schmidt, E. Thomas, S. J. Gibbs, A. Slujis, R. Zeebe, L. Kump, R. C. Martindale, S. E. Greene el al. 2012. The geological record of ocean acidification. Science 335: 1058-1063.

Huning, A. K., F. Melzner, J. Thomsen, M. A. Gutowska, L. Kramer, S. Frickenhaus, P. Rossenstiel, H. O. Portner, E. E. R. Philipp, and M. Lucassen. 2013. Impacts of seawater acidification on mantle gene expression patterns of the Baltic Sea blue mussel: implications for shell formation and energy metabolism. Mar. Biol. 160: 1845-1861.

IPCC (Intergovernmental Panel on Climate Change). 2007. Climate Change 2007: Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom.

Jones, B. M., D. M. Iglesias-Rodriguez, P. J. Skipp, R. J. Edwards, M. J. Greaves, J. R. Young, H. Elderfield, and D. C. O'Connor. 2013. Responses of the Emiliana huxleyi proteome to ocean acidification. PLoS One 8: e61868.

Kelly, M. W., J. L. Padilla-Gamino, and G. E. Hofmann. 2013. Natural variation and the capacity to adapt to ocean acidification in the keystone sea urchin Strongylocentrotus purpuratus. Glob. Change Biol. 19: 2536-2546.

Kroeker. K. J., R. L. Kordas, R. N. Crim, and G. S. Singh. 2010. Meta-analysis reveals negative yet variable effects of ocean acidification on marine organisms. Ecol. Lett. 13: 1419-1434.

Kroeker, K. J., F. Micheli, M. C. Gambi, and T. R. Martz. 2011. Divergent ecosystem responses within a benthic marine community to ocean acidification. Proc. Nat. Acad. Sci. USA 108: 14515-14520.

Kroeker, K. J., R. L. Kordas, R. Crim, 1. E. Hendriks, L. Ramajo, G. S. Singh, C. M. Duarte, and J. P. Gatusso. 2013. Impacts of ocean acidification on marine organisms: quantifying sensitivities and interaction with warming. Glob. Change Biol. 19: 1884-1896.

Kiiltz, D. 2003. Evolution of the cellular stress proteome: from monophyletic origin to ubiquitous function. J. Exp. Biol. 206: 3119-3124.

Kultz, D. 2005. Molecular and evolutionary basis of the cellular stress response. Annu. Rev. Physiol. 67: 225-257.

Kurihara, H., Y. Takano, D. Kurokawa, and K. Akasaka. 2012. Ocean acidification reduces biomineralization-related gene expression in the sea urchin, Hemicentrotus pulcherrimus. Mar. Biol. 159: 2819-2826.

Leninger, A. L. 1999. Principles of Biochemistry, 3rd ed., D. L. Nelson and M. M. Cox, eds. Worth, New York. Pp. 620-635, 708-721.

Leong, P. K. K., and D. T. Manahan. 1997. Metabolic importance of [Na.sup.+]/[K.sup.+]-ATPase activity during sea urchin development. J. Exp. Biol. 200: 2881-2892.

Li, X., F. Wu, and D. A. Beard. 2013. Identification of the kinetic mechanism of succinyl-CoA synthetase. BioSci. Rep. 33: 145-163.

Liu, W., Z. Huang, J. Lin, and M. He. 2012. Seawater acidification and elevated temperature affect gene expression patterns of the pearly oyster Pinctada fucata. PLoS One 7: e33679.

Livingston, B., C. Killian. F. Wilt, A. Cameron. M. Landrum. O.

Ermolaeva, V. Sapojnikov, D. Maglott, A. Buchanan, and C. Ettensohn. 2006. A genome-wide analysis of biomineralization-related proteins in the sea urchin Strongylocentrotus purpuratus. Dev. Biol. 300: 335-348.

Lockwood, B. L., and G. N. Somero. 2010. Transcriptomic responses to heat stress in invasive and native blue mussels (genus Mytilus): molecular correlates of invasive success. J. Exp. Biol. 213: 3548-3558.

Luthi, D., M. Le Floch, B. Bereiter, T. Blunier, J. M. Barnola, U. Siegenthaler, D. Raynaud, J. Jouzel, H. Fischer, K. Kawamura, and T. F. Stocker. 2008. High-resolution carbon dioxide concentration record 650,000-800,000 years before present. Nature 453: 379-82.

Martin, S., S. Richier, M. L. Pedrotti, S. Dupont, C. Castejon, Y. Gerakis, M. E. Kerros, F. Oberhansli, J. L. Teyssie, R. Jeffree, and

J. P. Gattusso. 2011. Early development and molecular plasticity in the Mediterranean sea urchin Paracentrotus lividis exposed to C[O.sub.2]-driven acidification. J. Exp. Biol. 214: 1357-1368.

Matozzo, V., A. Chinellato, M. Munari, L. Finos, M. Bressan, and M. G. Marin. 2012. First evidence of immunomodulation in bivalves under seawater acidification and increased temperature. PLoS One 7: e33820.

Mi, H., and P. Thomas. 2009. PANTHER Pathway: an ontology- based pathway database coupled with data analysis tools. Methods Mol. Biol. 563: 123-140.

Michaelidis, B., C. Ouzounis, A. Paleras, and H. O. Portner. 2005. Effects of long-term moderate hypercapnia on acid-base balance and growth rate in marine mussels Mytilus galloprovinciallis. Mar. Ecol. Prog. Ser. 293: 109-118.

Moya, A., L. Huisman, E. E. Ball, L. C. Grasso, C. M. Chua, H. N. Woo, J. P. Gatusso, S. Foret, and D. J. Miller. 2012. Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to C[O.sub.2]-driven acidification during the initiation of calcification. Mol. Ecol. 21: 2440-2454.

Nakamura, M., S. Ohki, A. Suzuki, and K. Sakai. 2011. Coral larvae under ocean acidification: survival, metabolism, and metamorphosis. PLoS One 6: e14521.

National Research Council. 2010. America's Climate Choices: Panel on Adapting to the Impacts of Climate Change. The National Academies Press, Washington DC. 325 pp.

O'Donnell, M. J., A. E. Todgham, M. A. Sewell, L. T. M. Hammond, K. Ruggiero, N. A. Fangue, M. L. Zippay, and G. E. Hofmann. 2010. Ocean acidification alters skeletogenesis and gene expression in larval sea urchins. Mar. Ecol. Prog. Ser. 398: 157-171.

O'Donnell, M. J., M. N. George, and E. Carrington. 2013. Mussel byssus attachment weakened by ocean acidification. Nat. Clim. Change. 3: 587-590.

Oliver, T. A., D. A. Garfield, M. K. Manier, R. Haygood, G. A. Wray,

and S. R. Palumbi. 2010. Whole-genome positive selection and habitat-driven evolution in a shallow and deep-sea urchin. Genome Biol. Evol. 2: 800-814.

Oliveri, P., D. M. Carrick, and E. H. Davidson. 2002. A regulatory gene network that directs micromere specification in the sea urchin embryo. Dev. Biol. 246: 209-228.

Oliveros, J. C. 2007. VENNY. An interactive tool for comparing lists with Venn Diagrams. [Online] Available: tools/venny/index.html [2014, April 28].

Padilla-Gamino, J. L., M. W. Kelly, T. G. Evans, and G. E. Hofmann. 2013. Temperature and C[O.sub.2] additively regulate physiology, morphology, and genomic responses of larval sea urchins, Strongylocentrotus purpuratus. Proc. R. Soc. B 280: 20130155.

Pearse, J. S. 2006. Ecological role of purple sea urchins. Science 314: 940-941.

Pennington, J. T., and R. R. Strathmann. 1990. Consequences of the calcite skeletons of plantonic echinoderm larvae for orientation, swimming, and shape. Biol. Bull. 179: 121-133.

Pespeni, M. H., B. T. Barney, and S. R. Palumbi. 2013a. Differences in the regulation of growth and biomineralization genes revealed through long-term common-garden acclimation and experimental genomics in the purple sea urchin. Evolution 67: 1901-1914.

Pespeni, M. H., E. Sanford, B. Gaylord, T. M. Hill, J. D. Hosfelt, H. K. Jaris, M. LaVigne, E. A. Lenz, A. D. Russell, M. K. Young, and S. R. Palumbi. 2013b. Evolutionary change during experimental ocean acidification. Proc. Natl. Acad. Sci. USA 110: 6937-6942.

Petrak, J., R. Ivanek, O. Toman, R. Cmejla, J. Cmejlova, D. Vyoral, J. Zivny, and C. D. Vulpe. 2008. Deja vu in proteomics: a hit parade of repeatedly identified differentially expressed proteins. Proteomics 8: 1744-1749.

Place, S. P., and B. W. Smith. 2012. Effects of seawater acidification on cell cycle control mechanisms in Strongylocentrotus purpuratus embryos. PLoS One 7: e34068.

Portner, H. 0.2008. Ecosystem effects of ocean acidification in times of ocean warming: a physiologist's view. Mar. Ecol. Prog. Ser. 373: 203-217.

Raven, J., K. Caldeira, H. Elderfield, O. Hoegh-Guldberg, P. Liss, U. Riebesell, J. Shepherd, C. Turley, and A. Watson. 2005. Ocean Acidification Due to Increasing Atmospheric Carbon Dioxide, Royal Society Policy Document. The Royal Society, London.

Reipschlager, A., and H. O. Portner. 1996. Metabolic depression during environmental stress: the role of extra- versus intracellular pH in Sipunculus nudus. J. Exp. Biol. 199: 1801-1807.

Renn, S. C. P., N. Aubin-Horth, and H. A. Hofmann. 2004. Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray. BMC Genomics 5: 42.

Ries, J. B., A. L. Cohen, and D. C. McCorkle. 2009. Marine calcifiers exhibit mixed responses to C[O.sub.2]-induced ocean acidification. Geology 37: 1131-1134.

Sabine, C. L., R. A. Feely, N. Gruber, R. M. Key, K. Lee, J. L. Bullistcr, R. Wanninkof, C. S. Wong, D. W. R. Wallace, B. Tilbrook et al. 2004. The oceanic sink for anthropogenic C[O.sub.2]. Science 305: 367-371.

Schalkhausser, B., C. Bock, K. Stemmer, T. Brey, H. O. Portner and G. Laning. 2012. Impact of ocean acidification on escape performance of the king scallop, Pecten maximus, from Norway. Mar. Biol. 160: 1995-2006.

Sheridan, J. A., and D. Bickford. 2011. Shrinking body size as an ecological response to climate change. Nat. Clim. Change 1: 401-406.

Shirayama, Y., and H. Thornton. 2005. Effect of increased atmospheric C[O.sub.2] on shallow water marine benthos. J. Geophys. Res. 110: C09S08.

Sirbu, A., G. Kerr, M. Crane, and H. J. Ruskin. 2012. RNA-Seq vs dual- and single-channel microarray data: sensitivity analysis for differential expression and clustering. PLoS One 7: e50986.

Sodergren, E., G. M. Weinstock, E. H. Davidson, A. R. Cameron, R. A. Gibbs, R. C. Angerer, L. M. Angerer, M. I. Arnone, D. R. Burgess, R. D. Burke et al., 2006. The genome of the sea urchin Strongylocentrotus purpuratus. Science 314: 941-952.

Strathmann, R. R. 1971. The feeding behavior of plantotrophic echinoderm larvae mechanisms: regulation and rates of suspension feeding. J. Exp. Mar. Biol. Ecol. 6: 109-160.

Strathmann, R. R., and D. Grunbaum. 2006. Good eaters, poor swimmers: compromises in larval form. Integr. Comp. Biol. 46: 312-322.

Stumpp, M., S. Dupont, M. C. Thorndyke and F. Melzner. 2011a. C[O.sub.2] induced seawater acidification impacts sea urchin larval development. II. Gene expression patterns in pluteus larvae. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 160: 320-330.

Stumpp, M., J. Wren, F. Melzner, M. C. Thorndyke, and S. T. Dupont. 2011b. C[O.sub.2] induced seawater acidification impacts sea urchin larval development. I. Elevated metabolic rates decrease scope for growth and induce developmental delay. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 160: 331-340.

Stumpp, M., M. Y. Hu, F. Melzner, M. A. Gutowska, N. Dorey, N. Himmerkus, W. C. Holtmann, S. T. Dupont, M. C. Thorndyke, and M. Bleich. 2012. Acidified seawater impacts sea urchin larvae pH regulatory systems relevant for calcification. Proc. Natl. Acad. Sci. USA 109: 18192-18197.

Tellis, M. S., M. M. Lauer, S. Nadella, A. Bianchini, and C. M. Wood. 2013. Ionic status, calcium uptake, and [Ca.sup.2+]-ATPase activity during early development in the purple sea urchin (Strongylocentrotus purpuratus). Comp. Biochem. Physiol. A Mol. Integr. Physiol. 166: 212-211.

Thomas, P. D., A. Kejarlwal, M. J. Campbell, H. Mi, D. N. Guo, I. Istvan, B. Ulitsky-Lazareva, A. Muruganujan, S. Rabkin, J. A. Vanfergriff et al. 2003. PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 31: 334-341.

Thomas, P. D., A. Kejariwal, N. Guo, N. H. Mi, M. J. Campbell, A. Muruganujan, and B. Ulitsky-Lazareva. 2006. Applications for protein sequence-function evolution data: mRNA/protein expression analysis and coding SNP scoring tools. Nucleic Acids Res. 34: W645-W650.

Thomsen, J., and F. Melzner. 2010. Moderate seawater acidification does not elicit long-term metabolic depression in the blue mussel Mytilus edulis. Mar Biol. 157: 2667-2676.

Todgham, A. E., and G. E. Hofmann. 2009. Transcriptomic response of sea urchin larvae Strongylocentrotus purpuratus to C[O.sub.2]-driven seawater acidification. J. Exp. Biot. 212: 2579-2594.

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

Tomanek, L., M. J. Zuzow, A. V. Ivanina, E. Beniash, and I. M. Sokolova. 2011. Proteomic response to elevated pC[O.sub.2] in eastern oysters, Crassostrea virginica: evidence for oxidative stress. J. Exp. Biol. 214: 1836-1844.

Tseng, G. C., M. K. Oh, L. Rohlin, and W. H. Wong. 2001. Issues in cDNA microarray analysis: quality filtering, channel normalization, models of variations and assessment of gene effects. Nucleic Acids Res. 29: 2549-2557.

UniProt Consoritum. 2014. Activities at the Universal Protein Resource (UniProt). Nucl. Acids. Res. 42: D191-D198.

Vidal-Dupiol, J., D. Zoccola, E. Tambutte, C. Grunau, C. Cosseau, K. M. Smith, M. Freltag, N. M. Dheilly, D. Allemand, and S. Tambutte. 2013. Genes related to ion-transport and energy production are upregulated in response to C[O.sub.2]-driven pH decrease in corals: new insights from transcriptome analysis. PLoS One 8: e58652.

Wang, P., F. G. Bouwman, and E. C. M. Mariman. 2009. Generally detected proteins in comparative proteomics--a matter of cellular stress response. Proteomics 9: 2955-2966.

Wang, Y., C. Barbacioru, F. Hyland, W. Xiao, K. L. Hunkapiller, J. Blake, F. Chan, C. Gonzalez, L. Zhang, and R. R. Samasha. 2006. Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC Genomics 7: 59.

Wang, Z., M. Gerstein, and M. Snyder. 2009. RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 10: 57-63.

Wei, Z., R. C. Angerer, and L. M. Angerer. 2006. A database of mRNA expression patterns for the sea urchin embryo. Dev. Biol. 300: 476-484.

Wheatly, M. G., and R. P. Henry. 1992. Extracellular and intracellular acid-base regulation in crustaceans. J. Exp. Zool. 263: 127-142.

Whitehead, A. 2012. Comparative genomics in ecological physiology: toward a more nuanced understanding of acclimation and adaptation. J. Exp. Biol. 215: 884-891.

Whitehead, A., and D. L. Crawford. 2006. Variation within and among species in gene expression: raw material for evolution. Mol. Ecol. 15: 1197-1211.

Wilt, F. H. 2002. Biomineralization of the spicules of sea urchin embryos. Zoo/. Sci. 19: 253-261.

Wittmann, A. C., and H. O. Portner. 2013. Sensitivities of extant animal taxa to ocean acidification. Nat. Clim. Change 3: 995-1001.

Wong, K. K. W., A. C. Lane, P. T. Y. Leung, and L. V. Thiyagarajan. 2011. Response of larval barnacle proteome to C[O.sub.2]-driven seawater acidification. Comp. Biochem. Physiol. D Genomics Proteomics. 6: 310-321.

Yu, P. C., P. G. Matson, T. R. Martz, and G. E. Hofmann. 2011. The ocean acidification seascape and its relationship to the performance of calcifying marine invertebrates: laboratory experiments on the development of urchin larvae framed by environmentally-relevant pC[O.sub.2]/pH. J. Exp. Mar. Biol. Ecol. 400: 288-295.

Yu, P. C., M. A. Sewell, P. G. Matson, E. B. Rivest, L. Kapsenberg, and G. E. Hofmann. 2013. Growth attenuation with developmental schedule progression in embryos and early larvae of Sterechinus neumayeri raised under elevated C[O.sub.2]. PLoS One 8: e52448.

Zhang, G., X. Fang, X. Guo, L. Li, F. Xu, P. Yang, L. Zhang, X. Wang, H. Qi, Z. Xiong et al. 2012. The oyster genome reveals stress adaptation and complexity of shell formation. Nature 490: 4954.


Department of Biological Sciences, California State University East Bay, Hayward, California 94542

Received 5 November 2013; accepted 27 March 2014.

* To whom correspondence should be addressed. E-mail: tyler.evans@

Tabic 1
List of publications used in analyses

    Publication             Species            Life stage

Evans et al. (2013)    Strongylocentrotus   gastrula,
                       purpuratus (2)       echinoplutei

Hammond and Hofmann    Strongylocentrotus   hatching, primary
(2012)                 purpuratus (3)       invagination,

Kurihara et al.        Hemicentrotus        gastrula, prism,
(2012)                 pulcherrimus         pluteus

Martin et al. (2011)   Paracentrotus        1-, 2- and 3-day
                       lividus              embryos

O'Donnell et al.       Lytechinus pictus    echinoplutei

Padilla-Gamino et      Strongylocentrotus   gastrula,
al. (2013)             purpuratus (3)       echinoplutei

Stumpp et al.          Strongylocentrotus   2, 4, 7 days past
(2011a)                purpuratus (4)       fertilization

Todgham and Hofmann    Strongylocentrotus   early prism
(2009)                 purpuratus (3)

                          pC[O.sub.2]       Method of detection
    Publication           ([micro]atm)              (1)

Evans et al. (2013)    435, 813, 1255       microarray (28,036

Hammond and Hofmann    380, 1000, 1350      qPCR (7 genes)

Kurihara et al.        380, 1000, 2000      qPCR (3 genes)

Martin et al. (2011)   Experiment 1: 397,   qPCR (6 genes)
                       743, 1188, 1962,
                       3555, 6632

                       Experiment 2: 412,
                       704, 1161, 1996,
                       3624, 6590

O'Donnell et al.       380, 540, 970        microarray (1057
(2010)                                      genes)

Padilla-Gamino et      400, 1100            microarray (28,036
al. (2013)                                  genes)

Stumpp et al.          399, 1318            qPCR (27 genes)

Todgham and Hofmann    380, 540, 1020       microarray (1057
(2009)                                      genes)

(1) Number of genes assayed in parenthesis.

(2) Collected from Fogarty Creek, Oregon.

(3) Collected from Santa Barbara, California.

(4) Collected from Newport Beach, California.

Table 2
All differentially expressed genes involved in urchin

GLEAN Gene Model                      Gene Name

SPU_000062         "Tubulin beta chain (Beta tubulin)"

SPU_000142         "PREDICTED: similar to collagen XVIII
                   [Strongylocentrotus purpuratus]"

SPU_000310         "PREDICTED: similar to Fit-1 tyrosine kinase

SPU_000826         "PREDICTED: similar to SM30-alpha protein
                   [Strongylocentrotus purpuratus]"

SPU_000827         "SM30-alpha protein [Strongylocentrotus

SPU_002788         "tubulin, beta 4 [Xenopus tropicalisr

SPU_003768         "3 alpha procollagen [Strongylocentrotus

SPU_004135         "Unknown (protein for MGC:139465) [Bos taurus]"

SPU_004551         "winged helix transcription factor Forkhead-1
                   [Strongylocentrotus purpuratus]"

SPU_004867         "PREDICTED: similar to Spicule matrix protein
                   SM30 [Strongylocentrotus purpuratus]"

SPU_004876         "PREDICTED: similar to ferritin
                   [Strongylocentrotus purpuratus]"

SPU_004885         "PREDICTED: similar to apolipoprotein D; apoD
                   [Strongylocentrotus purpuratus]"

SPU_005296         "PREDICTED: similar to H(+)-transporting ATPase
                   beta subunit"

SPU_005990         "spicule matrix protein SM29 [Strongylocentrotus

SPU_006387         "PREDICTED: similar to mspl30 protein

                   [Strongylocentrotus purpuratus]"

SPU_007231         "PREDICTED: similar to tetraspanin family

SPU_007484         "cyclophilin 1 [Strongylocentrotus purpuratus]"

SPU_007882         "C-type lectin domain protein
                   [Strongylocentrotus purpuratus]"

SPU_008305         "PREDICTED: similar to cyclophilin
                   [Strongylocentrotus purpuratus]"

SPU_009076         "PREDICTED: similar to 6a fibrillar collagen
                   [Strongylocentrotus purpuratus]"

SPU_009549         "PREDICTED: similar to Sdcbp-prov protein
                   [Strongylocentrotus purpuratus]"

SPU_012518         "PREDICTED: similar to carbonic anhydrase
                   [Strongylocentrotus purpuratus]"

SPU_012548         "PREDICTED: similar to secreted protein acidic
                   and rich in cysteine"

SPU_013821         "PREDICTED: similar to Mesenchyme-specific cell
                   surface glycoprotein"

SPU_013823         "PREDICTED: similar to mspl30 protein
                   [Strongylocentrotus purpuratus]"

SPU_014492         "PREDICTED: similar to primary mesenchyme
                   specific protein"

SPU_015276         "PREDICTED: similar to ubiquitin
                   [Strongylocentrotus purpuratus]"

SPU_015708         "alpha2(IV)-like collagen [Strongylocentrotus

SPU_017585         "PREDICTED: similar to extracellular
                   calcium-sensing receptor"

SPU_017605         "mitochondrial ATP synthase alpha subunit

SPU_017641         "similar to transmembrane receptor [Oikopleura

SPU_018405         "similar to roundabout homolog 4"

SPU_018408         "biomineralization protein SpP16
                   [Strongylocentrotus purpuratus]"

SPU_018811         "spicule matrix protein SM50 [Strongylocentrotus

SPU_022116         "PREDICTED: similar to ENSANGP00000020852
                   [Strongylocentrotus purpuratus]"

SPU_022346         "PREDICTED: similar to CA XIV, partial
                   [Strongylocentrotus purpuratus]"

SPU_024809         "PREDICTED: similar to CA XIV, partial
                   [Strongylocentrotus purpuratus]"

SPU_025235         "PREDICTED: similar to Egf16-prov protein
                   [Strongylocentrotus purpuratus]"

SPU 026483         "PREDICTED: similar to Carbonic anhydrase 8
                   [Strongylocentrotus purpuratus]"

SPU_026949         "PREDICTED: similar to melanotransferrin/EOS47
                   [Strongylocentrotus purpuratus]"

SPU_027906         "Secreted trypsin-like serine protease [Hahella
                   chejuensis KCTC 2396]"

SPU_028275         "osteonectin [Strongylocentrotus purpuratus]"

SPU_030147         "Primary mesenchyme-specific protein PM27"

Table 3
Significantly over-represented metabolism and
ion transport ontologies among the set of all
differentially expressed genes

                                        Number of       Number
                                        Genes in       of Genes
                                        Reference   Differentially
            Gene Oncology                 List        Expressed

Biological Process
  ion transport                            365            54
  cation transport                         313            46
  transport                               1424           114
  generation of precursor
    metabolites and energy                 203            28
  carbohydrate metabolic process           460            44
  extracellular transport                  69             13
  oxidative phosphorylation                47             10
  respiratory electron transport
    chain                                  182            22
  amino acid transport                     58             11
  metabolic process                       4486           268
  tricarboxylic acid cycle                 24             7
Molecular Function
  transporter activity                     554            77
  transmembrane transporter activity       520            70
  cation transmembrane transporter
    activity                               197            29
  oxidoreductase activity                  444            47
  hydrogen ion transmembrane
    transporter activity                   29             8
  carbohydrate transmembrane
    transporter activity                   23             10
Protein Class
  transporter                              607            77
  cation transporter                       135            28
  oxidoreductase                           443            47
  ATP synthase                             29             8
  carbohydrate transporter                 23             7

            Gene Oncology                P-Value       Rank

Biological Process
  ion transport                          1.08E-09    2 of 15
  cation transport                       5.53E-08    3 of 15
  transport                              1.07E-04    5 of 15
  generation of precursor
    metabolites and energy               5.99E-04    6 of 15
  carbohydrate metabolic process         1.23E-02    8 of 15
  extracellular transport                1.35E-02    9 of 15
  oxidative phosphorylation              3.43E-02    11 of 15
  respiratory electron transport
    chain                                4.24E-02    12 of 15
  amino acid transport                   4.46E-02    13 of 15
  metabolic process                      4.46E-02    14 of 15
  tricarboxylic acid cycle               4.91E-02    15 of 15
Molecular Function
  transporter activity                   4.02E-13    1 of 15
  transmembrane transporter activity     4.07E-11    2 of 15
  cation transmembrane transporter
    activity                             9.84E-05    4 of 15
  oxidoreductase activity                4.69E-04    5 of 15
  hydrogen ion transmembrane
    transporter activity                 2.33E-02    7 of 15
  carbohydrate transmembrane
    transporter activity                 3.32E-02    8 of 15
Protein Class
  transporter                            3.13E-11    1 of 15
  cation transporter                     1.67E-07    3 of 15
  oxidoreductase                         5.34E-04    4 of 15
  ATP synthase                           2.81E-02    7 of 15
  carbohydrate transporter               4.01E-02    8 of 15

Table 4
Significantly over-represented ontologies among
the set of genes differentially expressed in at
least two studies

                            Number of      Number of
                            Genes in         Genes        Bonferroni
                            Reference   Differentially    Corrected
      Gene Ontology           List         Expressed       P-Value

Biological Process
  homeostatic process          60              5           1.45E-04
    phosphorylation            47              4           1.87E-03
  generation of
    metabolites and
    energy                     203             6           3.91E-03
  cellular calcium ion
    homeostasis                23              3           7.51E-03
  respiratory electron
    transport chain            182             5           2.91E-02
  cation transport             313             6           4.23E-02
Molecular Function
  transporter activity         554             9           2.62E-03
    transporter activity       520             8           1.29E-02
  cation transmembrane
    transporter activity       197             5           3.66E-02
Protein Class
  cation transporter           135             5           7.56E-03
  transporter                  601             8           4.25E-02
COPYRIGHT 2014 University of Chicago Press
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Evans, Tyler G.; Watson-Wynn, Priscilla
Publication:The Biological Bulletin
Date:Jun 1, 2014
Previous Article:The skeleton of postmetamorphic echinoderms in a changing world.
Next Article:Effects of ocean acidification on population dynamics and community structure of crustose coralline algae.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters