Ocean acidification research in the 'post-genomic' era: Roadmaps from the purple sea urchin Strongylocentrotus purpuratus

Dynamic responses during early development of the sea urchin Strongylocentrotus intermedius to CO2-driven ocean acidification: A microRNA-mRNA integrated analysis

To explore the dynamic molecular responses to CO2-driven ocean acidification (OA) during the early developmental stages of sea urchins, gametes of Strongylocentrotus intermedius were fertilized and developed to the four-armed larva stage in either natural seawater (as a control; pHNBS = 7.99 ± 0.01) or acidified conditions (ΔpHNBS = -0.3, -0.4, and - 0.5 units) according to the prediction for ocean pH by the end of this century. Specimens from five developmental stages (fertilization, cleavage, blastula, prism, and four-armed larva) were collected and comparative microRNA (miRNA) and mRNA transcriptome analyses were performed. The results showed that 1) a total of 22,224 differentially expressed genes (DEGs) and 51 differentially expressed miRNAs (DEMs) were identified in the OA-treated groups compared with the control group. 2) The numbers of both DEGs and DEMs were the largest at the blastula stage, indicating dramatic changes in gene expression. 3) Five "miR-1/DEG" modules were identified as potential biomarkers reflecting the response of sea urchins to OA during the early developmental period. 4) The PI3K/Akt signaling pathway was a key pathway involved in the response of S. intermedius to OA in its early developmental stages. This study deepens our understanding of the dynamic molecular regulatory mechanisms underlying sea urchin responses to CO2-driven OA.

Molecular response to CO2-driven ocean acidification in the larvae of the sea urchin Hemicentrotus pulcherrimus: Evidence from comparative transcriptome analyses

In order to explore the impact of CO2-driven ocean acidification (OA) on gene expression of sea urchins, gametes of Hemicentrotus pulcherrimus were fertilized and developed to the four-armed larvae in either seawater at current pH levels (pHNBS = 7.98) or in three laboratory-controlled OA conditions (ΔpHNBS = -0.3, -0.4, -0.5 units) based on the projections of the Intergovernmental Panel on Climate Change (IPCC) for 2100. Four-armed larval specimens were collected, and comparative transcriptome analysis was then performed. The results showed that 58 differentially expressed genes (DEGs) were identified in OA-treated groups as compared to the control. Moreover, more transition and transversion SNPs were observed in OA-treated groups than those in the control indicating a potential occurrence of adaption to OA in H. pulcherrimus larvae. Six candidate DEGs shared among OA-treated groups were identified as potential biomarkers correlated with low pH tolerance, mainly enriched in nine pathways associated with Notch signaling, dorso-ventral axis formation, oxidative phosphorylation, lysine degradation, valine, leucine and isoleucine degradation, lysosome, cell adhesion molecules, glutathione metabolism and PPAR signaling pathway. These results will not only enrich our knowledge of the impacts of OA on sea urchin larvae from the aspect of gene expression, provide a better understanding on larval forms coping with OA, but also offer more clues and biomarkers for developing protection or management strategies for sea urchins under near-future OA conditions.

Effect of elevated temperature, sea water acidification, and phenanthrene on the expression of genes involved in the shell and pearl formation of economic pearl oyster (Pinctada radiata)

Our study aims to examine the effect of some stressors on the gene expression levels of shell matrix proteins in a pearl oyster. Oysters were exposed to the different combinations of the temperature, pH, and phenanthrene concentration is currently measured in the Persian Gulf and the predicted ocean warming and acidification for 28 days. The expression of all the studied genes was significantly downregulated. Time and temperature had the greatest effects on the decreases in n19 and n16 genes expression, respectively. Aspein and msi60 genes expression were highly influenced by pH. Pearlin was affected by double interaction temperature and phenanthrene. Moreover, a correlation was observed among the expression levels of studied genes. This study represents basic data on the relationship between mRNA transcription genes involved in the shell and pearl formation and climate changes in pollutant presence conditions and acclimatizing mechanism of the oyster to the future scenario as well.

Transcriptomes reveal the involved genes in the sea urchin Mesocentrotus nudus exposed to high flow velocities

Despite the importance of flow velocity in marine ecosystems, molecular mechanisms of the water flow induced behavioral and growth changes remain largely unknown in sea urchins. The present study compared the gene expressions of the sea urchin Mesocentrotus nudus at high flow velocities (10 cm/s and 20 cm/s) and low flow velocity (2 cm/s) using transcriptomes. A total of 490 and 470 differentially expressed genes (DEGs) were discovered at 10 cm/s and 20 cm/s, respectively. There were 235 up-regulated and 255 down-regulated genes at 10 cm/s, 213 up-regulated and 257 down-regulated genes at 20 cm/s, compared with sea urchins at 2 cm/s. Further, there were 72 overlapped DEGs involved in regulation at both 10 cm/s and 20 cm/s. Gene Ontology (GO) functional annotation showed that DEGs were mainly enriched to cellular process, cell part, binding, and metabolism process. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis found that DEGs were enriched in three pathways related to amino acid metabolism and lipid metabolism. A number of genes related to growth and metabolism of sea urchins were mobilized in high flow velocity environment. We further highlighted a muscle-associated gene ankyrin-1, which is correlated with the movement of tube feet at different flow velocities. The present study provides valuable information on the molecular mechanisms of changed behaviors and growth when sea urchins are exposed to high flow velocity.

Molecular mechanisms of biomineralization in marine invertebrates

Much recent marine research has been directed towards understanding the effects of anthropogenic-induced environmental change on marine biodiversity, particularly for those animals with heavily calcified exoskeletons, such as corals, molluscs and urchins. This is because life in our oceans is becoming more challenging for these animals with changes in temperature, pH and salinity. In the future, it will be more energetically expensive to make marine skeletons and the increasingly corrosive conditions in seawater are expected to result in the dissolution of these external skeletons. However, initial predictions of wide-scale sensitivity are changing as we understand more about the mechanisms underpinning skeletal production (biomineralization). These studies demonstrate the complexity of calcification pathways and the cellular responses of animals to these altered conditions. Factors including parental conditioning, phenotypic plasticity and epigenetics can significantly impact the production of skeletons and thus future population success. This understanding is paralleled by an increase in our knowledge of the genes and proteins involved in biomineralization, particularly in some phyla, such as urchins, molluscs and corals. This Review will provide a broad overview of our current understanding of the factors affecting skeletal production in marine invertebrates. It will focus on the molecular mechanisms underpinning biomineralization and how knowledge of these processes affects experimental design and our ability to predict responses to climate change. Understanding marine biomineralization has many tangible benefits in our changing world, including improvements in conservation and aquaculture and exploitation of natural calcified structure design using biomimicry approaches that are aimed at producing novel biocomposites.

Maintenance of a Genetic Cline in the Barnacle Balanus glandula

The barnacle Balanus glandula is a broadly distributed species in the temperate northeastern Pacific that is notable for a robust genetic cline between about 36° and 40° N latitude. Prior work established the evolutionary origins of this pattern and proposed that it is maintained by environmental selection. In recent years, "climate velocity" studies in marine habitats have shown dramatic distributional shifts for many species as they track their preferred temperature range in a warming ocean. We re-sampled B. glandula across its entire geographic range to determine whether there has been any shift in this genetic distribution, a development signaling that temperature or other climate factors are maintaining this genetic cline. Additionally, we asked whether the spatially distributed mitochondrial lineages also vary in reproductive output with latitude, using location as a proxy for temperature and other coastal environmental factors. Here we show that although the distribution of the genetic cline has not appreciably changed, there is a notable association of decreased reproductive output at lower latitudes of the distribution in the "northern" lineage of B. glandula.

Transcriptomes reveal genes involved in covering and sheltering behaviors of the sea urchin Strongylocentrotus intermedius exposed to UV-B radiation

Although the potential link exists between behavioral responses to UV-B radiation and the maximization of fitness, molecular mechanisms of these UV-B induced behaviors remain poorly understood. For the first time, we investigated the transcriptomes of covered (CB), sheltered (SB) and non-protected (NA) sea urchins Strongylocentrotus intermedius exposed to UV-B radiation. A total of 330 differentially expressed genes were revealed by transcriptome comparisons. By comparing with the group NA, we found 79 up-regulated and 118 down-regulated genes in SB group, as well as 26 up-regulated and 67 down-regulated genes in group CB. There were 34 up-regulated genes and 52 down-regulated genes in group SB, compared with group CB. These differentially expressed genes failed to enrich either Gene Ontology (GO) or Kyoto Encyclopedia of Genes and Genomes (KEGG), only except an enrichment in KEGG. We highlighted TRPA1 and Opsin as key neurobiological genes involved in the molecular mechanisms of covering and sheltering behaviors of sea urchins exposed to UV-B radiation. What's more, other identified genes provide valuable resources for future investigations on the molecular basis of covering and sheltering behaviors of sea urchins.

Effect of ocean acidification on growth, calcification, and gene expression in the pearl oyster, Pinctada fucata

In this study, shell growth, shell microstructure, and expression levels of shell matrix protein genes (aspein, n16, and nacrein) that play a key role in the CaCO₃ crystal polymorphism (calcite and aragonite) of the shell were investigated in the pearl oyster Pinctada fucata at pH 8.10, 7.70, and 7.40. We found that the shell length and total weight index did not vary significantly between oysters reared at pH 8.10 and 7.70, but was significantly lower at pH 7.40. Calcium content and shell hardness were not significantly different between pH 8.10 and 7.70, but were significantly different at pH 7.40. At pH 7.40, the shell exhibited a poorly organized nacreous microstructure, and showed an apparent loss of structural integrity in the nacreous layer. The prismatic layer appeared morphologically dissimilar from the samples at pH 8.10 and 7.70. The internal layer was corroded and had dissolved. At pH 7.40, the expression levels of nacrein, aspein, and n16 decreased on day 1, and remained low between days 2 and 42. The expression levels of these genes were significantly lower at pH 7.40 than at pH 8.10 and 7.70 during days 2-42. These results suggest that ocean acidification will have a limited impact on shell growth, calcification, and associated gene expression levels at a pH of 7.70, which is projected to be reached by the end of the century. The negative effects were found on calcification and gene expression occurred at the lowest experimental pH (7.40).

De novo transcriptome sequencing and comparative analysis to discover genes involved in ovarian maturity in Strongylocentrotus nudus

Strongylocentrotus nudus is an edible sea urchin, mainly harvested in China. Correlation studies indicated that S. nudus with larger diameter have a prolonged marketing time and better palatability owing to their precocious gonads and extended maturation process. However, the molecular mechanism underlying this phenomenon is still unknown. Here, transcriptome sequencing was applied to study the ovaries of adult S. nudus with different shell diameters to explore the possible mechanism. In this study, four independent cDNA libraries were constructed, including two from the big size urchins and two from the small ones using a HiSeq™2500 platform. A total of 88,581 unigenes were acquired with a mean length of 1354bp, of which 66,331 (74.88%) unigenes could be annotated using six major publicly available databases. Comparative analysis revealed that 353 unigenes were differentially expressed (with log2(ratio)≥1, FDR≤0.001) between the two groups. Of these, 20 differentially expressed genes (DEGs) were selected to confirm the accuracy of RNA-seq data by quantitative real-time RT-PCR. Furthermore, gene ontology and KEGG pathway enrichment analyses were performed to find the putative genes and pathways related to ovarian maturity. Eight unigenes were identified as significant DEGs involved in reproduction related pathways; these included Mos, Cdc20, Rec8, YP30, cytochrome P450 2U1, ovoperoxidase, proteoliaisin, and rendezvin. Our research fills the gap in the studies on the S. nudus ovaries using transcriptome analysis.

Transcriptomic responses to seawater acidification among sea urchin populations inhabiting a natural pH mosaic

Increasing awareness of spatial and temporal variation in ocean pH suggests some marine populations may be adapted to local pH regimes and will therefore respond differently to present-day pH variation and to long-term ocean acidification. In the Northeast Pacific Ocean, differences in the strength of coastal upwelling cause latitudinal variation in prevailing pH regimes that are hypothesized to promote local adaptation and unequal pH tolerance among resident populations. In this study, responses to experimental seawater acidification were compared among embryos and larvae from six populations of purple sea urchins (Strongylocentrotus purpuratus) inhabiting areas that differ in their frequency of low pH exposure and that prior research suggests are locally adapted to seawater pH. Transcriptomic analyses demonstrate urchin populations most frequently exposed to low pH seawater responded to experimental acidification by expressing genes within major ATP-producing pathways at greater levels than populations encountering low pH less often. Multiple genes within the tricarboxylic acid cycle, electron transport chain and fatty acid beta oxidation pathways were upregulated in urchin populations experiencing low pH conditions most frequently. These same metabolic pathways were significantly over-represented among genes both expressed in a population-specific manner and putatively under selection to enhance low pH tolerance. Collectively, these data suggest natural selection is acting on metabolic gene networks to redirect ATP toward maintaining acid-base homeostasis and enhance tolerance of seawater acidification. As a trade-off, marine populations more tolerant of low pH may have less energy to put towards other aspects of fitness and to respond to additional ocean change.

Regional adaptation defines sensitivity to future ocean acidification

Physiological responses to temperature are known to be a major determinant of species distributions and can dictate the sensitivity of populations to global warming. In contrast, little is known about how other major global change drivers, such as ocean acidification (OA), will shape species distributions in the future. Here, by integrating population genetics with experimental data for growth and mineralization, physiology and metabolomics, we demonstrate that the sensitivity of populations of the gastropod Littorina littorea to future OA is shaped by regional adaptation. Individuals from populations towards the edges of the natural latitudinal range in the Northeast Atlantic exhibit greater shell dissolution and the inability to upregulate their metabolism when exposed to low pH, thus appearing most sensitive to low seawater pH. Our results suggest that future levels of OA could mediate temperature-driven shifts in species distributions, thereby influencing future biogeography and the functioning of marine ecosystems.

Global warming is expected to lead to shifts in species' geographic ranges to track preferred temperatures. Here, the authors show that populations of the common periwinkle vary in their sensitivity to ocean acidification, another major global change driver, which could further restrict range shifts caused by warming.

Quantitative analysis of oyster larval proteome provides new insights into the effects of multiple climate change stressors

The metamorphosis of planktonic larvae of the Pacific oyster (Crassostrea gigas) underpins their complex life-history strategy by switching on the molecular machinery required for sessile life and building calcite shells. Metamorphosis becomes a survival bottleneck, which will be pressured by different anthropogenically induced climate change-related variables. Therefore, it is important to understand how metamorphosing larvae interact with emerging climate change stressors. To predict how larvae might be affected in a future ocean, we examined changes in the proteome of metamorphosing larvae under multiple stressors: decreased pH (pH 7.4), increased temperature (30 °C), and reduced salinity (15 psu). Quantitative protein expression profiling using iTRAQ-LC-MS/MS identified more than 1300 proteins. Decreased pH had a negative effect on metamorphosis by down-regulating several proteins involved in energy production, metabolism, and protein synthesis. However, warming switched on these down-regulated pathways at pH 7.4. Under multiple stressors, cell signaling, energy production, growth, and developmental pathways were up-regulated, although metamorphosis was still reduced. Despite the lack of lethal effects, significant physiological responses to both individual and interacting climate change related stressors were observed at proteome level. The metamorphosing larvae of the C. gigas population in the Yellow Sea appear to have adequate phenotypic plasticity at the proteome level to survive in future coastal oceans, but with developmental and physiological costs.

Sea Star Wasting Disease in the Keystone Predator Pisaster ochraceus in Oregon: Insights into Differential Population Impacts, Recovery, Predation Rate, and Temperature Effects from Long-Term Research

Sea star wasting disease (SSWD) first appeared in Oregon in April 2014, and by June had spread to most of the coast. Although delayed compared to areas to the north and south, SSWD was initially most intense in north and central Oregon and spread southward. Up to 90% of individuals showed signs of disease from June-August 2014. In rocky intertidal habitats, populations of the dominant sea star Pisaster ochraceus were rapidly depleted, with magnitudes of decline in density among sites ranging from -2x to -9x (59 to 84%) and of biomass from -2.6x to -15.8x (60 to 90%) by September 2014. The frequency of symptomatic individuals declined over winter and persisted at a low rate through the spring and summer 2015 (~5-15%, at most sites) and into fall 2015. Disease expression included six symptoms: initially with twisting arms, then deflation and/or lesions, lost arms, losing grip on substrate, and final disintegration. SSWD was disproportionally higher in orange individuals, and higher in tidepools. Although historically P. ochraceus recruitment has been low, from fall 2014 to spring 2015 an unprecedented surge of sea star recruitment occurred at all sites, ranging from ~7x to 300x greater than in 2014. The loss of adult and juvenile individuals in 2014 led to a dramatic decline in predation rate on mussels compared to the previous two decades. A proximate cause of wasting was likely the "Sea Star associated Densovirus" (SSaDV), but the ultimate factors triggering the epidemic, if any, remain unclear. Although warm temperature has been proposed as a possible trigger, SSWD in Oregon populations increased with cool temperatures. Since P. ochraceus is a keystone predator that can strongly influence the biodiversity and community structure of the intertidal community, major community-level responses to the disease are expected. However, predicting the specific impacts and time course of change across west coast meta-communities is difficult, suggesting the need for detailed coast-wide investigation of the effects of this outbreak.