Transcriptomic responses to ocean acidification in larval sea urchins from a naturally variable pH environment

Ocean acidification alleviated nickel toxicity to a marine copepod under multigenerational scenarios but at a cost with a loss of transcriptome plasticity during recovery

Marine ecosystem has been experiencing multiple stressors caused by anthropogenic activities, including ocean acidification (OA) and nickel (Ni) pollution. Here, we examined the individual/combined effects of OA (pCO2 1000 μatm) and Ni (6 μg/L) exposure on a marine copepod Tigriopus japonicus for six generations (F1-F6), followed by one-generation recovery (F7) in clean seawater. Ni accumulation and several important phenotypic traits were measured in each generation. To explore within-generation response and transgenerational plasticity, we analyzed the transcriptome profile for the copepods of F6 and F7. The results showed that Ni exposure compromised the development, reproduction and survival of copepods during F1-F6, but its toxicity effects were alleviated by OA. Thus, under OA and Ni combined exposure, due to their antagonistic interaction, the disruption of Ca2+ homeostasis, and the inhibition of calcium signaling pathway and oxytocin signaling pathway were not found. However, as a cost of acclimatization/adaption potential to long-term OA and Ni combined exposure, there was a loss of transcriptome plasticity during recovery, which limited the resilience of copepods to previously begin environments. Overall, our work fosters a comprehensive understanding of within- and transgenerational effects of climatic stressor and metal pollution on marine biota.

Transgenerational acclimation to acidified seawater and gene expression patterns in a sea urchin

Transgenerational responses of susceptible calcifying organisms to progressive ocean acidification are an important issue in reducing uncertainty of future predictions. In this study, a two-generation rearing experiment was conducted using mature Mesocentrotus nudus, a major edible sea urchin that occurs along the coasts of northern Japan. Morphological observations and comprehensive gene expression analysis (RNA-seq) of resulting larvae were performed to examine transgenerational acclimation to acidified seawater. Two generations of rearing experiments showed that larvae derived from parents acclimated to acidified seawater tended to have higher survival and show less reduction in body size when exposed to acidified seawater of the same pH, suggesting that a positive carry-over effect occurred. RNA-seq analysis showed that gene expression patterns of larvae originated from both acclimated and non-acclimated parents to acidified seawater tended to be different than control condition, and the gene expression pattern of larvae originated from acclimated parents was substantially different than that of larvae of non-acclimated and control parents.

Transcriptomic Responses of a Lightly Calcified Echinoderm to Experimental Seawater Acidification and Warming during Early Development

Ocean acidification (OA) and ocean warming (OW) are potential obstacles to the survival and growth of marine organisms, particularly those that rely on calcification. This study investigated the single and joint effects of OA and OW on sea cucumber Apostichopus japonicus larvae raised under combinations of two temperatures (19 °C or 22 °C) and two pCO2 levels (400 or 1000 μatm) that reflect the current and end-of-21st-century projected ocean scenarios. The investigation focused on assessing larval development and identifying differences in gene expression patterns at four crucial embryo-larval stages (blastula, gastrula, auricularia, and doliolaria) of sea cucumbers, using RNA-seq. Results showed the detrimental effect of OA on the early development and body growth of A. japonicus larvae and a reduction in the expression of genes associated with biomineralization, skeletogenesis, and ion homeostasis. This effect was particularly pronounced during the doliolaria stage, indicating the presence of bottlenecks in larval development at this transition phase between the larval and megalopa stages in response to OA. OW accelerated the larval development across four stages of A. japonicus, especially at the blastula and doliolaria stages, but resulted in a widespread upregulation of genes related to heat shock proteins, antioxidant defense, and immune response. Significantly, the negative effects of elevated pCO2 on the developmental process of larvae appeared to be mitigated when accompanied by increased temperatures at the expense of reduced immune resilience and increased system fragility. These findings suggest that alterations in gene expression within the larvae of A. japonicus provide a mechanism to adapt to stressors arising from a rapidly changing oceanic environment.

Using meta-analysis to explore the roles of global upwelling exposure and experimental design in bivalve responses to low pH

Low pH conditions, associated with ocean acidification, represent threats to many commercially and ecologically important organisms, including bivalves. However, there are knowledge gaps regarding factors explaining observed differences in biological responses to low pH in laboratory experiments. Specific sources of local adaptation such as upwelling exposure and the role of experimental design, such as carbonate chemistry parameter changes, should be considered. Linking upwelling exposure, as an individual oceanographic phenomenon, to responses measured in laboratory experiments may further our understanding of local adaptation to global change. Here, meta-analysis is used to test the hypotheses that upwelling exposure and experimental design affect outcomes of individual, laboratory-based studies that assess bivalve metabolic (clearance and respiration rate) responses to low pH. Results show that while bivalves generally decrease metabolic activity in response to low pH, upwelling exposure and experimental design can significantly impact outcomes. Bivalves from downwelling or weak upwelling areas decrease metabolic activity in response to low pH, but bivalves from strong upwelling areas increase or do not change metabolic activity in response to low pH. Furthermore, experimental temperature, exposure time and magnitude of the change in carbonate chemistry parameters all significantly affect outcomes. These results suggest that bivalves from strong upwelling areas may be less sensitive to low pH. This furthers our understanding of local adaptation to global change by demonstrating that upwelling alone can explain up to 49 % of the variability associated with bivalve metabolic responses to low pH. Furthermore, when interpreting outcomes of individual, laboratory experiments, scientists should be aware that higher temperatures, shorter exposure times and larger changes in carbonate chemistry parameters may increase the chance of suppressed metabolic activity.

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

Parentage influence on gene expression under acidification revealed through single-embryo sequencing

The dissolution of anthropogenic carbon dioxide (CO2 ) in seawater has altered its carbonate chemistry in the process of ocean acidification (OA). OA affects the viability of marine species. In particular, calcifying organisms and their early planktonic larval stages are considered vulnerable. These organisms often utilize energy reserves for metabolism rather than growth and calcification as supported by bulk RNA-sequencing (RNA-seq) experiments. Yet, transcriptomic profiling of a bulk sample reflects the average gene expression of the population, neglecting the variations between individuals, which forms the basis for natural selection. Here, we used single-embryo RNA-seq on larval sea urchin Heliocidaris crassispina, which is a commercially and ecologically valuable species in East Asia, to document gene expression changes to OA at an individual and family level. Three paternal half-sibs groups were fertilized and exposed to 3 pH conditions (ambient pH 8.0, 7.7 and 7.4) for 12 h prior to sequencing and oxygen consumption assay. The resulting transcriptomic profile of all embryos can be distinguished into four clusters, with differences in gene expressions that govern biomineralization, cell differentiation and patterning, as well as metabolism. While these responses were influenced by pH conditions, the male identities also had an effect. Specifically, a regression model and goodness of fit tests indicated a significant interaction between sire and pH on the probability of embryo membership in different clusters of gene expression. The single-embryo RNA-seq approach is promising in climate stressor research because not only does it highlight potential impacts before phenotypic changes were observed, but it also highlights variations between individuals and lineages, thus enabling a better determination of evolutionary potential.

Transcriptome Responses to Different Environments in Intertidal Zones in the Peanut Worm Sipunculus nudus

The peanut worm (Sipunculus nudus) is an important intertidal species worldwide. Species living in the same aquaculture area might suffer different environmental impacts. To increase knowledge of the molecular mechanisms underlying the response to environmental fluctuations, we performed a transcriptome analysis of S. nudus from different intertidal zones using a combination of the SMRT platform and the Illumina sequencing platform. (1) A total of 105,259 unigenes were assembled, and 23,063 unigenes were perfectly annotated. The results of the PacBio Iso-Seq and IIIumina RNA-Seq enriched the genetic database of S. nudus. (2) A total of 830 DEGs were detected in S. nudus from the different groups. In particular, 33 DEGs had differential expression in the top nine KEGG pathways related to pathogens, protein synthesis, and cellular immune response and signaling. The results indicate that S. nudus from different zones experience different environmental stresses. (3) Several DEGs (HSPA1, NFKBIA, eEF1A, etc.) in pathways related to pathogens (influenza A, legionellosis, measles, and toxoplasmosis) had higher expression in groups M and L. HSPA1 was clearly enriched in most of the pathways, followed by NFKBIA. The results show that the peanut worms from the M and L tidal flats might have suffered more severe environmental conditions. (4) Some DEGs (MKP, MRAS, and HSPB1) were upregulated in peanut worms from the H tidal flat, and these DEGs were mainly involved in the MAPK signaling pathway. These results indicate that the MAPK pathway may play a vital role in the immune response of the peanut worm to the effects of different intertidal flats. This study provides a valuable starting point for further studies to elucidate the molecular basis of the response to different environmental stresses in S. nudus.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria

The increasing concentration of CO₂ in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO₂ and reared under acidified (pH ~ 7.5, pCO₂ ~ 1200 ppm) or control (pH ~ 7.9, pCO₂ ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

Sex and gametogenesis stage are strong drivers of gene expression in Mytilus edulis exposed to environmentally relevant plasticiser levels and pH 7.7

Plastic pollution and changes in oceanic pH are both pressing environmental issues. Little emphasis, however, has been placed on the influence of sex and gametogenesis stage when investigating the effects of such stressors. Here, we examined histology and molecular biomarkers of blue mussels Mytilus edulis exposed for 7 days to a pH 7.7 scenario (- 0.4 units) in combination with environmentally relevant concentrations (0, 0.5 and 50 µg/L) of the endocrine disrupting plasticiser di-2-ethylhexyl phthalate (DEHP). Through a factorial design, we investigated the gametogenesis cycle and sex-related expression of genes involved in pH homeostasis, stress response and oestrogen receptor-like pathways after the exposure to the two environmental stressors. As expected, we found sex-related differences in the proportion of developing, mature and spawning gonads in histological sections. Male gonads also showed higher levels of the acid-base regulator CA2, but females had a higher expression of stress response-related genes (i.e. sod, cat, hsp70). We found a significant effect of DEHP on stress response-related gene expression that was dependent on the gametogenesis stage, but there was only a trend towards downregulation of CA2 in response to pH 7.7. In addition, differences in gene expression between males and females were most pronounced in experimental conditions containing DEHP and/or acidified pH but never the control, indicating that it is important to consider sex and gametogenesis stage when studying the response of mussels to diverse stressors.

A Transcriptomic Analysis of Phenotypic Plasticity in Crassostrea virginica Larvae under Experimental Acidification

Ocean acidification (OA) is a major threat to marine calcifiers, and little is known regarding acclimation to OA in bivalves. This study combined physiological assays with next-generation sequencing to assess the potential for recovery from and acclimation to OA in the eastern oyster (Crassostrea virginica) and identify molecular mechanisms associated with resilience. In a reciprocal transplant experiment, larvae transplanted from elevated pCO₂ (~1400 ppm) to ambient pCO₂ (~350 ppm) demonstrated significantly lower mortality and larger size post-transplant than oysters remaining under elevated pCO₂ and had similar mortality compared to those remaining in ambient conditions. The recovery after transplantation to ambient conditions demonstrates the ability for larvae to rebound and suggests phenotypic plasticity and acclimation. Transcriptomic analysis supported this hypothesis as genes were differentially regulated under OA stress. Transcriptomic profiles of transplanted and non-transplanted larvae terminating in the same final pCO₂ converged, further supporting the idea that acclimation underlies resilience. The functions of differentially expressed genes included cell differentiation, development, biomineralization, ion exchange, and immunity. Results suggest acclimation as a mode of resilience to OA. In addition, the identification of genes associated with resilience can serve as a valuable resource for the aquaculture industry, as these could enable marker-assisted selection of OA-resilient stocks.

A calcification-related calmodulin-like protein in the oyster Crassostrea gigas mediates the enhanced calcium deposition induced by CO2 exposure

Calcium transportation and homeostasis are essential for marine bivalves to maintain basic metabolism and build their shells. Calmodulin-like proteins (CaLPs) are important calcium sensors and buffers and can respond to ocean acidification (OA) in marine calcifiers. However, no further study of their physiological function in calcium metabolism under elevated CO₂ has been performed. Here, we identified a novel CaLP (designated CgCaLP) in the Pacific oyster Crassostrea gigas and demonstrated its participation in the calcification process: the mRNA expression level of CgCaLP peaked at the trochophore larval stage and remained high at stages when shells were shaped; the mRNA and protein of CgCaLP were more highly expressed in mantle tissue than in other tissues. Under elevated CO₂ levels, the protein expression level of CgCaLP in hemocytes increased, while in contrast, significantly decreased protein levels were detected in gill and mantle tissues. Shell dissolution caused the imbalance of calcium in hemocytes and decreased calcium absorption and transportation demand in gill and mantle tissues, inducing the molecular function allocation of CgCaLP under CO₂ exposure. Despite the decreased protein level in mantle tissue, CgCaLP was found to translocate to outer mantle epithelium (OME) cells where condensed calcium-rich deposits (CRDs) were detected. We further demonstrated that CgCaLP mRNA and protein expression levels could respond to seawater Ca²⁺ availability, suggesting that the calcium deposition capacity of oysters might be enhanced to fight against shell dissolution problems and that CgCaLP might serve as an essential participator of the process. In summary, CgCaLP might enhance calcium deposition under CO₂ exposure and thus play a significant and flexible molecular function involved in a compensation strategy of oysters to fight against the acidified ocean.

Assessment of the juvenile vulnerability of symbiont-bearing giant clams to ocean acidification

Ocean acidification (OA) severely affects marine bivalves, especially their calcification processes. However, very little is known about the fate of symbiont-bearing giant clams in the acidified oceans, which hinders our ability to develop strategies to protect this ecologically and economically important group in coral reef ecosystems. Here, we explored the integrated juvenile responses of fluted giant clam Tridacna squamosa (Lamarck, 1819) to acidified seawater at different levels of biological organization. Our results revealed that OA did not cause a significant reduction in survival and shell growth performance, indicating that T. squamosa juveniles are tolerated to moderate acidification. Yet, significantly reduced net calcification rate demonstrated the calcifying physiology sensitivity to OA, in line with significant declines in symbiont photosynthetic yield and zooxanthellae density which in turn lowered the amount of energy supply for energetically expensive calcification processes. Subsequent transcriptome sequencing and comparative analysis of differentially expressed genes revealed that the regulation of calcification processes, such as transport of calcification substrates, acid-base regulation, synthesis of organic matrix in the calcifying fluid, as well as metabolic depression were the major response to OA. Taken together, the integration of physiological and molecular responses can provide a comprehensive understanding of how the early life history stages of giant clams respond to OA and make an important leap forward in assessing their fate under future ocean conditions.

PAHs and PCBs Affect Functionally Intercorrelated Genes in the Sea Urchin Paracentrotus lividus Embryos

Polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) represent the most common pollutants in the marine sediments. Previous investigations demonstrated short-term sublethal effects of sediments polluted with both contaminants on the sea urchin Paracentrotus lividus after 2 months of exposure in mesocosms. In particular, morphological malformations observed in P. lividus embryos deriving from adults exposed to PAHs and PCBs were explained at molecular levels by de novo transcriptome assembly and real-time qPCR, leading to the identification of several differentially expressed genes involved in key physiological processes. Here, we extensively explored the genes involved in the response of the sea urchin P. lividus to PAHs and PCBs. Firstly, 25 new genes were identified and interactomic analysis revealed that they were functionally connected among them and to several genes previously defined as molecular targets of response to the two pollutants under analysis. The expression levels of these 25 genes were followed by Real Time qPCR, showing that almost all genes analyzed were affected by PAHs and PCBs. These findings represent an important further step in defining the impacts of slight concentrations of such contaminants on sea urchins and, more in general, on marine biota, increasing our knowledge of molecular targets involved in responses to environmental stressors.

Transgenerational plasticity responses of oysters to ocean acidification differ with habitat

Transgenerational plasticity (TGP) has been identified as a critical mechanism of acclimation that may buffer marine organisms against climate change, yet whether the TGP response of marine organisms is altered depending on their habitat is unknown. Many marine organisms are found in intertidal zones where they experience episodes of emersion (air exposure) daily as the tide rises and recedes. During episodes of emersion, the accumulation of metabolic carbon dioxide (CO2) leads to hypercapnia for many species. How this metabolic hypercapnia impacts the TGP response of marine organisms to climate change is unknown as all previous transgenerational studies have been done under subtidal conditions, where parents are constantly immersed. Here, we assess the capacity of the ecologically and economically important oyster, Saccostrea glomerata, to acclimate to elevated CO2 dependent on habitat, across its vertical distribution, from the subtidal to intertidal zone. Tidal habitat altered both the existing tolerance and transgenerational response of S. glomerata to elevated CO2. Overall, larvae from parents conditioned in an intertidal habitat had a greater existing tolerance to elevated CO2 than larvae from parents conditioned in a subtidal habitat, but had a lower capacity for beneficial TGP following parental exposure to elevated CO2. Our results suggest that the TGP responses of marine species will not be uniform across their distribution and highlights the need to consider the habitat of a species when assessing TGP responses to climate change stressors.

Unique Genomic and Phenotypic Responses to Extreme and Variable pH Conditions in Purple Urchin Larvae

Environmental variation experienced by a species across space and time can promote the maintenance of genetic diversity that may be adaptive in future global change conditions. Selection experiments have shown that purple sea urchin, Strongylocentrotus purpuratus, populations have adaptive genetic variation for surviving pH conditions at the "edge" (pH 7.5) of conditions experienced in nature. However, little is known about whether populations have genetic variation for surviving low-pH events beyond those currently experienced in nature or how variation in pH conditions affects organismal and genetic responses. Here, we quantified survival, growth, and allele frequency shifts in experimentally selected developing purple sea urchin larvae in static and variable conditions at three pH levels: pH 8.1 (control), pH 7.5 (edge-of-range), and pH 7.0 (extreme). Variable treatments recovered body size relative to static treatments, but resulted in higher mortality, suggesting a potential tradeoff between survival and growth under pH stress. However, within each pH level, allele frequency changes were overlapping between static and variable conditions, suggesting a shared genetic basis underlying survival to mean pH regardless of variability. In contrast, genetic responses to pH 7.5 (edge) versus pH 7.0 (extreme) conditions were distinct, indicating a unique genetic basis of survival. In addition, loci under selection were more likely to be in exonic regions than regulatory, indicating that selection targeted protein-coding variation. Loci under selection in variable pH 7.5 conditions, more similar to conditions periodically experienced in nature, performed functions related to lipid biosynthesis and metabolism, while loci under selection in static pH 7.0 conditions performed functions related to transmembrane and mitochondrial processes. While these results are promising in that purple sea urchin populations possess genetic variation for surviving extreme pH conditions not currently experienced in nature, they caution that increased acidification does not result in a linear response but elicits unique physiological stresses and survival mechanisms.

Effect of Fermentation on Nutritional Quality, Growth and Hematological Parameters of Rats Fed Sorghum-Soybean-Orange flesh Sweet Potato Complementary Diet

The protein quality of complementary foods developed from fermented and unfermented sorghum, soybeans, and orange-fleshed sweet potato (OFSP) flour blends was evaluated using rat model. The test diet was as follows: UF2: unfermented sorghum (56%), soybean (17%), and OFSP (27%); UF3: unfermented sorghum (59%), soybean (31%), and OFSP (10%); F2: fermented sorghum (56%), soybean (17%), and OFSP (27%); and F3: fermented sorghum (59%), soybean (31%), and OFSP (10%), while cerelac served as positive control, corn starch (basal diet), and ogi (negative control). Forty-nine Wistar albino rats were grouped and fed with diets for 28 days. The growth, hematological, serum parameters of animals, protein quality, and proximate composition of developed diet were determined. Fermentation significantly improved the protein content and nutritional indices of experimental animals. Moisture content ranged from 2.5% to 9.24%, protein (7.09%-25.29%), ash (1.09%-3.71%), fat (10.28%-15.24%), and fiber (0.85%-3.17%). The biological values (BV) ranged from 75.11% to 78.44%. The weight gained in rat fed the formulated diet ranged from 46.0 g to 77.3 g and was highest in F3. The packed cell volume (PCV), hemoglobin concentration (HBC), red blood cell (RBC), and lymphocytes were highest in F3. Urea nitrogen and creatinine of the rats fed with formulated diets ranged from 3.58 to 15.32 mg/dl and 1.56 to 6.15 mg/dl, respectively. Sample F3 is a protein-rich complementary food that is comparable to ogi and suitable to manage malnutrition and support growth in children. However, clinical trials on the formulated diet are needed to further substantiate its nutritional potentials.

Gene Expression Changes after Parental Exposure to Metals in the Sea Urchin Affect Timing of Genetic Programme of Embryo Development

It is widely accepted that phenotypic traits can be modulated at the epigenetic level so that some conditions can affect the progeny of exposed individuals. To assess if the exposure of adult animals could result in effects on the offspring, the Mediterranean sea urchin and its well-characterized gene regulatory networks (GRNs) was chosen as a model. Adult animals were exposed to known concentrations of zinc and cadmium (both individually and in combination) for 10 days, and the resulting embryos were followed during the development. The oxidative stress occurring in parental gonads, embryo phenotypes and mortality, and the expression level of a set of selected genes, including members of the skeletogenic and endodermal GRNs, were evaluated. Increased oxidative stress at F₀, high rates of developmental aberration with impaired gastrulation, in association to deregulation of genes involved in skeletogenesis (dri, hex, sm50, p16, p19, msp130), endodermal specification (foxa, hox11/13b, wnt8) and epigenetic regulation (kat2A, hdac1, ehmt2, phf8 and UBE2a) occurred either at 24 or 48 hpf. Results strongly indicate that exposure to environmental pollutants can affect not only directly challenged animals but also their progeny (at least F₁), influencing optimal timing of genetic programme of embryo development, resulting in an overall impairment of developmental success.

Gene expression patterns of red sea urchins (Mesocentrotus franciscanus) exposed to different combinations of temperature and pCO2 during early development

BACKGROUND

The red sea urchin Mesocentrotus franciscanus is an ecologically important kelp forest herbivore and an economically valuable wild fishery species. To examine how M. franciscanus responds to its environment on a molecular level, differences in gene expression patterns were observed in embryos raised under combinations of two temperatures (13 °C or 17 °C) and two pCO2 levels (475 μatm or 1050 μatm). These combinations mimic various present-day conditions measured during and between upwelling events in the highly dynamic California Current System with the exception of the 17 °C and 1050 μatm combination, which does not currently occur. However, as ocean warming and acidification continues, warmer temperatures and higher pCO2 conditions are expected to increase in frequency and to occur simultaneously. The transcriptomic responses of the embryos were assessed at two developmental stages (gastrula and prism) in light of previously described plasticity in body size and thermotolerance under these temperature and pCO2 treatments.

RESULTS

Although transcriptomic patterns primarily varied by developmental stage, there were pronounced differences in gene expression as a result of the treatment conditions. Temperature and pCO2 treatments led to the differential expression of genes related to the cellular stress response, transmembrane transport, metabolic processes, and the regulation of gene expression. At each developmental stage, temperature contributed significantly to the observed variance in gene expression, which was also correlated to the phenotypic attributes of the embryos. On the other hand, the transcriptomic response to pCO2 was relatively muted, particularly at the prism stage.

CONCLUSIONS

M. franciscanus exhibited transcriptomic plasticity under different temperatures, indicating their capacity for a molecular-level response that may facilitate red sea urchins facing ocean warming as climate change continues. In contrast, the lack of a robust transcriptomic response, in combination with observations of decreased body size, under elevated pCO2 levels suggest that this species may be negatively affected by ocean acidification. High present-day pCO2 conditions that occur due to coastal upwelling may already be influencing populations of M. franciscanus.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

Transcriptional changes of Pacific oyster Crassostrea gigas reveal essential role of calcium signal pathway in response to CO2-driven acidification

There is increasing evidence that ocean acidification (OA) has a significant impact on marine organisms. However, the ability of most marine organisms to acclimate to OA and the underlying mechanisms are still not well understood. In the present study, whole transcriptome analysis was performed to compare the impacts of short- (7 days, named as short group) and long- (60 days, named as long group) term CO₂ exposure (pH 7.50) on Pacific oyster Crassostrea gigas. The responses of C. gigas to short- and long-term CO₂ exposure shared common mechanisms in metabolism, membrane-associated transportation and binding processes. Long-term CO₂ exposure induced significant expression of genes involved in DNA or RNA binding, indicating the activated transcription after long-term CO₂ exposure. Oysters in the short-term group underwent significant intracellular calcium variation and oxidative stress. In contrast, the intracellular calcium, ROS level in hemocytes and H₂O₂ in serum recovered to normal levels after long-term CO₂ exposure, suggesting the compensation of physiological status and mutual interplay between calcium and oxidative level. The compensation was supported by the up-regulation of a series of calcium binding proteins (CBPs) and calmodulins (CaMs) related signal pathway. The results provided valuable information to understand the molecular mechanism underlying the responses of Pacific oyster to the acidified ocean and might have implications for predicting the possible effects of global climate changes on oyster aquaculture.

Biogeography of ocean acidification: Differential field performance of transplanted mussels to upwelling-driven variation in carbonate chemistry

Ocean acidification (OA) represents a serious challenge to marine ecosystems. Laboratory studies addressing OA indicate broadly negative effects for marine organisms, particularly those relying on calcification processes. Growing evidence also suggests OA combined with other environmental stressors may be even more deleterious. Scaling these laboratory studies to ecological performance in the field, where environmental heterogeneity may mediate responses, is a critical next step toward understanding OA impacts on natural communities. We leveraged an upwelling-driven pH mosaic along the California Current System to deconstruct the relative influences of pH, ocean temperature, and food availability on seasonal growth, condition and shell thickness of the ecologically dominant intertidal mussel Mytilus californianus. In 2011 and 2012, ecological performance of adult mussels from local and commonly sourced populations was measured at 8 rocky intertidal sites between central Oregon and southern California. Sites coincided with a large-scale network of intertidal pH sensors, allowing comparisons among pH and other environmental stressors. Adult California mussel growth and size varied latitudinally among sites and inter-annually, and mean shell thickness index and shell weight growth were reduced with low pH. Surprisingly, shell length growth and the ratio of tissue to shell weight were enhanced, not diminished as expected, by low pH. In contrast, and as expected, shell weight growth and shell thickness were both diminished by low pH, consistent with the idea that OA exposure can compromise shell-dependent defenses against predators or wave forces. We also found that adult mussel shell weight growth and relative tissue mass were negatively associated with increased pH variability. Including local pH conditions with previously documented influences of ocean temperature, food availability, aerial exposure, and origin site enhanced the explanatory power of models describing observed performance differences. Responses of local mussel populations differed from those of a common source population suggesting mussel performance partially depended on genetic or persistent phenotypic differences. In light of prior research showing deleterious effects of low pH on larval mussels, our results suggest a life history transition leading to greater resilience in at least some performance metrics to ocean acidification by adult California mussels. Our data also demonstrate “hot” (more extreme) and “cold” (less extreme) spots in both mussel responses and environmental conditions, a pattern that may enable mitigation approaches in response to future changes in climate.

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.

Molecular adaptation of molluscan biomineralisation to high-CO2 oceans - The known and the unknown

High-CO₂ induced ocean acidification (OA) reduces the calcium carbonate (CaCO₃) saturation level (Ω) and the pH of oceans. Consequently, OA is causing a serious threat to several ecologically and economically important biomineralising molluscs. Biomineralisation is a highly controlled biochemical process by which molluscs deposit their calcareous structures. In this process, shell matrix proteins aid the nucleation, growth and assemblage of the CaCO₃ crystals in the shell. These molluscan shell proteins (MSPs) are, ultimately, responsible for determination of the diverse shell microstructures and mechanical strength. Recent studies have attempted to integrate gene and protein expression data of MSPs with shell structure and mechanical properties. These advances made in understanding the molecular mechanism of biomineralisation suggest that molluscs either succumb or adapt to OA stress. In this review, we discuss the fate of biomineralisation process in future high-CO₂ oceans and its ultimate impact on the mineralised shell's structure and mechanical properties from the perspectives of limited substrate availability theory, proton flux limitation model and the omega myth theory. Furthermore, studying the interplay of energy availability and differential gene expression is an essential first step towards understanding adaptation of molluscan biomineralisation to OA, because if there is a need to change gene expression under stressors, any living system would require more energy than usual. To conclude, we have listed, four important future research directions for molecular adaptation of molluscan biomineralisation in high-CO₂ oceans: 1) Including an energy budgeting factor while understanding differential gene expression of MSPs and ion transporters under OA. 2) Unraveling the genetic or epigenetic changes related to biomineralisation under stressors to help solving a bigger picture about future evolution of molluscs, and 3) Understanding Post Translational Modifications of MSPs with and without stressors. 4) Understanding carbon uptake mechanisms across taxa with and without OA to clarify the OA theories on Ω.

Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey

For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.

Ocean Acidification and Coastal Marine Invertebrates: Tracking CO2 Effects from Seawater to the Cell

In the last few decades, numerous studies have investigated the impacts of simulated ocean acidification on marine species and communities, particularly those inhabiting dynamic coastal systems. Despite these research efforts, there are many gaps in our understanding, particularly with respect to physiological mechanisms that lead to pathologies. In this review, we trace how carbonate system disturbances propagate from the coastal environment into marine invertebrates and highlight mechanistic links between these disturbances and organism function. We also point toward several processes related to basic invertebrate biology that are severely understudied and prevent an accurate understanding of how carbonate system dynamics influence organismic homeostasis and fitness-related traits. We recommend that significant research effort be directed to studying cellular phenotypes of invertebrates acclimated or adapted to elevated seawater pCO₂ using biochemical and physiological methods.

The flagellar protein Enkurin is required for mouse sperm motility and for transport through the female reproductive tract

Enkurin was identified initially in mouse sperm where it was suggested to act as an intracellular adaptor protein linking membrane calcium influx to intracellular signaling pathways. In order to examine the function of this protein, a targeted mutation was introduced into the mouse Enkurin gene. Males that were homozygous for this mutated allele were subfertile. This was associated with lower rates of sperm transport in the female reproductive tract, including reduced entry into the oviduct and slower migration to the site of fertilization in the distal oviduct, and with poor progressive motility in vitro. Flagella from wild-type animals exhibited symmetrical bending and progressive motility in culture medium, and demembranated flagella exhibited the "curlicue" response to Ca2+ in vitro. In contrast, flagella of mice homozygous for the mutated allele displayed only asymmetric bending, nonprogressive motility, and a loss of Ca2+-responsiveness following demembrantion. We propose that Enkurin is part of a flagellar Ca2+-sensor that regulates bending and that the motility defects following mutation of the locus are the proximate cause of subfertility.

Intraspecific variation in the response of the estuarine European isopod Cyathura carinata (Krøyer, 1847) to ocean acidification

In the present study the model isopod, Cyathura carinata were exposed to four pH_NIST treatments (control: 7.9; 7.5, 7, 6.5) in order to determine the tolerance and pH threshold value this estuarine species withstand under future acidification scenarios. Seawater acidification significantly affected the lifespan of C. carinata, where population density was remarkably reduced at the lowest pH treatment. The longevity, survivorship and swimming activity (related to the acquisition of energy) of these isopods decreased with decreasing pH. Furthermore, to determine the possible metabolic plasticity of this species, the swimming activity, the Na+/K + -ATPase activity (relevant for osmoregulation process), and the RNA:DNA ratio (an indicator of fitness) were measure from two populations of C. carinata, one inhabiting a stable environment (pH_NIST 7.5-8.0) and one inhabiting a fluctuating pCO₂ regimes (pH 3.3-8.5) subjected to three pH treatments (7.9, 7.0 and 6.5). The population from high fluctuating pCO₂ conditions showed capacity to withstand to pH 6.5, as well as higher longevity and metabolic plasticity, when compared with the population from the habitat with slight pCO₂ variation. These results indicate that Cyathura population from stable environments could be vulnerable to ocean acidification because it could trigger detrimental effects on its survival energy budget, and growth. However, ocean acidification has limited effect on the energy budget and survival of C. carinata population from highly variable habitats, suggesting that they are able to cope with the elevated energy demand. The difference showed between populations is likely an indication of genetic differentiation in tolerance to ocean acidification, possibly attributable to local adaptations, which could provide the raw material necessary for adaptation to future conditions. In addition, our results suggest that when assessing marine crustacean responses to changing environments on a global scale, variability in population and metabolic responses need to be considered.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Little evidence of adaptation potential to ocean acidification in sea urchins living in "Future Ocean" conditions at a CO2 vent

Ocean acidification (OA) can be detrimental to calcifying marine organisms, with stunting of invertebrate larval development one of the most consistent responses. Effects are usually measured by short-term, within-generation exposure, an approach that does not consider the potential for adaptation. We examined the genetic response to OA of larvae of the tropical sea urchin Echinometra sp. C. raised on coral reefs that were either influenced by CO2 vents (pH ~ 7.9, future OA condition) or nonvent control reefs (pH 8.2). We assembled a high quality de novo transcriptome of Echinometra embryos (8 hr) and pluteus larvae (48 hr) and identified 68,056 SNPs. We tested for outlier SNPs and functional enrichment in embryos and larvae raised from adults from the control or vent sites. Generally, highest F ST values in embryos were observed between sites (intrinsic adaptation, most representative of the gene pool in the spawned populations). This comparison also had the highest number of outlier loci (40). In the other comparisons, classical adaptation (comparing larvae with adults from the control transplanted to either the control or vent conditions) and reverse adaptation (larvae from the vent site returned to the vent or explanted at the control), we only observed modest numbers of outlier SNPs (6-19) and only enrichment in two functional pathways. Most of the outliers detected were silent substitutions without adaptive potential. We conclude that there is little evidence of realized adaptation potential during early development, while some potential (albeit relatively low) exists in the intrinsic gene pool after more than one generation of exposure.

Rare genetic variation and balanced polymorphisms are important for survival in global change conditions

Standing genetic variation is important for population persistence in extreme environmental conditions. While some species may have the capacity to adapt to predicted average future global change conditions, the ability to survive extreme events is largely unknown. We used single-generation selection experiments on hundreds of thousands of Strongylocentrotus purpuratus sea urchin larvae generated from wild-caught adults to identify adaptive genetic variation responsive to moderate (pH 8.0) and extreme (pH 7.5) low-pH conditions. Sequencing genomic DNA from pools of larvae, we identified consistent changes in allele frequencies across replicate cultures for each pH condition and observed increased linkage disequilibrium around selected loci, revealing selection on recombined standing genetic variation. We found that loci responding uniquely to either selection regime were at low starting allele frequencies while variants that responded to both pH conditions (11.6% of selected variants) started at high frequencies. Loci under selection performed functions related to energetics, pH tolerance, cell growth and actin/cytoskeleton dynamics. These results highlight that persistence in future conditions will require two classes of genetic variation: common, pH-responsive variants maintained by balancing selection in a heterogeneous environment, and rare variants, particularly for extreme conditions, that must be maintained by large population sizes.

Expression of calcification-related ion transporters during blue mussel larval development

The physiological processes driving the rapid rates of calcification in larval bivalves are poorly understood. Here, we use a calcification substrate-limited approach (low dissolved inorganic carbon, C T) and mRNA sequencing to identify proteins involved in bicarbonate acquisition during shell formation. As a secondary approach, we examined expression of ion transport and shell matrix proteins (SMPs) over the course of larval development and shell formation. We reared four families of Mytilus edulis under ambient (ca. 1865 µmol/kg) and low C T (ca. 941 µmol/kg) conditions and compared expression patterns at six developmental time points. Larvae reared under low C T exhibited a developmental delay, and a small subset of contigs was differentially regulated between ambient and low C T conditions. Of particular note was the identification of one contig encoding an anion transporter (SLC26) which was strongly upregulated (2.3-2.9 fold) under low C T conditions. By analyzing gene expression profiles over the course of larval development, we are able to isolate sequences encoding ion transport and SMPs to enhance our understanding of cellular pathways underlying larval calcification processes. In particular, we observe the differential expression of contigs encoding SLC4 family members (sodium bicarbonate cotransporters, anion exchangers), calcium-transporting ATPases, sodium/calcium exchangers, and SMPs such as nacrein, tyrosinase, and transcripts related to chitin production. With a range of candidate genes, this work identifies ion transport pathways in bivalve larvae and by applying comparative genomics to investigate temporal expression patterns, provides a foundation for further studies to functionally characterize the proteins involved in larval calcification.

Differential responses to ocean acidification between populations of Balanophyllia elegans corals from high and low upwelling environments

Ocean acidification (OA), the global decrease in surface water pH from absorption of anthropogenic CO₂ , may put many marine taxa at risk. However, populations that experience extreme localized conditions, and are adapted to these conditions predicted in the global ocean in 2,100, may be more tolerant to future OA. By identifying locally adapted populations, researchers can examine the mechanisms used to cope with decreasing pH. One oceanographic process that influences pH is wind-driven upwelling. Here we compare two Californian populations of the coral Balanophyllia elegans from distinct upwelling regimes, and test their physiological and transcriptomic responses to experimental seawater acidification. We measured respiration rates, protein and lipid content, and gene expression in corals from both populations exposed to pH levels of 7.8 and 7.4 for 29 days. Corals from the population that experiences lower pH due to high upwelling maintained the same respiration rate throughout the exposure. In contrast, corals from the low upwelling site had reduced respiration rates, protein content and lipid-class content at low pH exposure, suggesting they have depleted their energy reserves. Using RNA-Seq, we found that corals from the high upwelling site upregulated genes involved in calcium ion binding and ion transport, most likely related to pH homeostasis and calcification. In contrast, corals from the low upwelling site downregulated stress response genes at low pH exposure. Divergent population responses to low pH observed in B. elegans highlight the importance of multi-population studies for predicting a species' response to future OA.

Aberrant gene expression profiles in Mediterranean sea urchin reproductive tissues after metal exposures

Marine organisms are simultaneously exposed to numerous pollutants, among which metals probably represent the most abundant in marine environments. In order to evaluate the effects of metal exposure at molecular level in reproductive tissues, we profiled the sea urchin transcriptional response after non-lethal exposures using pathway-focused mRNA expression analyses. Herein, we show that exposures to relatively high concentrations of both essential and toxic metals hugely affected the gonadic expression of several genes involved in stress-response, detoxification, transcriptional and post-transcriptional regulation, without significant changes in gonadosomatic indices. Even though treatments did not result in reproductive tissues visible alterations, metal exposures negatively affected the main mechanisms of stress-response, detoxification and survival of adult P. lividus. Additionally, transcriptional changes observed in P. lividus gonads may cause altered gametogenesis and maintenance of heritable aberrant epigenetic effects. This study leads to the conclusion that exposures to metals, as usually occurs in polluted coastal areas, may affect sea urchin gametogenesis, thus supporting the hypothesis that parental exposure to environmental stressors affects the phenotype of the offspring.

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.

Biotic and environmental stress induces nitration and changes in structure and function of the sea urchin major yolk protein toposome

The major yolk protein toposome plays crucial roles during gametogenesis and development of sea urchins. We previously found that nitration of toposome increases in the gonads of a Paracentrotus lividus population living in a marine protected area affected by toxic blooms of Ostreospsis cf. ovata, compared to control populations. This modification is associated with ovatoxin accumulation, high levels of nitric oxide in the gonads, and a remarkable impairment of progeny development. However, nothing is known about the environmental-mediated-regulation of the structure and biological function of toposome. Here, we characterize through wide-ranging biochemical and structural analyses the nitrated toposome of sea urchins exposed to the bloom, and subsequently detoxified. The increased number of nitrated tyrosines in toposome of sea urchins collected during algal bloom induced structural changes and improvement of the Ca²⁺-binding affinity of the protein. After 3 months’ detoxification, ovatoxin was undetectable, and the number of nitric oxide-modified tyrosines was reduced. However, the nitration of specific residues was irreversible and occurred also in embryos treated with metals, used as a proxy of environmental pollutants. The structural and functional changes of toposome caused by nitration under adverse environmental conditions may be related to the defective development of sea urchins’ progeny.

Effects of ocean acidification on copepods

Ocean acidification (OA) leads to significant changes in seawater carbon chemistry, broadly affects marine organisms, and considered as a global threat to the fitness of marine ecosystems. Due to the crucial role of copepods in marine food webs of transferring energy from primary producers to higher trophic levels, numerous studies have been conducted to examine the impacts of OA on biological traits of copepods such as growth and reproduction. Under OA stress, the copepods demonstrated species-specific and stage-dependent responses. Notably, different populations of the same copepod species demonstrated different sensitivities to the increased pCO₂. In copepods, the deleterious effects of OA are also reinforced by other naturally occurring co-stressors (e.g., thermal stress, food deprivation, and metal pollution). Given that most OA stress studies have focused on the effects of short-term exposure (shorter than a single generation), experiments using adults might have underestimated the damaging effects of OA and the long-term multigenerational exposure to multiple stressors (e.g., increased pCO₂ and food shortage) will be required. Particularly, omics-based technologies (e.g., genomics, proteomics, and metabolomics) will be helpful to better understand the underlying processes behind biological responses (e.g., survival, development, and offspring production) at the mechanistic level which will improve our predictions of the responses of copepods to climate change stressors including OA.

Transcriptomic response of the Antarctic pteropod Limacina helicina antarctica to ocean acidification

Background

Ocean acidification (OA), a change in ocean chemistry due to the absorption of atmospheric CO₂ into surface oceans, challenges biogenic calcification in many marine organisms. Ocean acidification is expected to rapidly progress in polar seas, with regions of the Southern Ocean expected to experience severe OA within decades. Biologically, the consequences of OA challenge calcification processes and impose an energetic cost.

Results

In order to better characterize the response of a polar calcifier to conditions of OA, we assessed differential gene expression in the Antarctic pteropod, Limacina helicina antarctica. Experimental levels of pCO₂ were chosen to create both contemporary pH conditions, and to mimic future pH expected in OA scenarios. Significant changes in the transcriptome were observed when juvenile L. h. antarctica were acclimated for 21 days to low-pH (7.71), mid-pH (7.9) or high-pH (8.13) conditions. Differential gene expression analysis of individuals maintained in the low-pH treatment identified down-regulation of genes involved in cytoskeletal structure, lipid transport, and metabolism. High pH exposure led to increased expression and enrichment for genes involved in shell formation, calcium ion binding, and DNA binding. Significant differential gene expression was observed in four major cellular and physiological processes: shell formation, the cellular stress response, metabolism, and neural function. Across these functional groups, exposure to conditions that mimic ocean acidification led to rapid suppression of gene expression.

Conclusions

Results of this study demonstrated that the transcriptome of the juvenile pteropod, L. h. antarctica, was dynamic and changed in response to different levels of pCO₂. In a global change context, exposure of L. h. antarctica to the low pH, high pCO₂ OA conditions resulted in a suppression of transcripts for genes involved in key physiological processes: calcification, metabolism, and the cellular stress response. The transcriptomic response at both acute and longer-term acclimation time frames indicated that contemporary L. h. antarctica may not have the physiological plasticity necessary for adaptation to OA conditions expected in future decades. Lastly, the differential gene expression results further support the role of shelled pteropods such as L. h. antarctica as sentinel organisms for the impacts of ocean acidification.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4161-0) contains supplementary material, which is available to authorized users.

Regulation of gene expression is associated with tolerance of the Arctic copepod Calanus glacialis to CO2-acidified sea water

Ocean acidification is the increase in seawater pCO ₂ due to the uptake of atmospheric anthropogenic CO ₂, with the largest changes predicted to occur in the Arctic seas. For some marine organisms, this change in pCO ₂, and associated decrease in pH, represents a climate change-related stressor. In this study, we investigated the gene expression patterns of nauplii of the Arctic copepod Calanus glacialis cultured at low pH levels. We have previously shown that organismal-level performance (development, growth, respiration) of C. glacialis nauplii is unaffected by low pH. Here, we investigated the molecular-level response to lowered pH in order to elucidate the physiological processes involved in this tolerance. Nauplii from wild-caught C. glacialis were cultured at four pH levels (8.05, 7.9, 7.7, 7.5). At stage N6, mRNA was extracted and sequenced using RNA-seq. The physiological functionality of the proteins identified was categorized using Gene Ontology and KEGG pathways. We found that the expression of 151 contigs varied significantly with pH on a continuous scale (93% downregulated with decreasing pH). Gene set enrichment analysis revealed that, of the processes downregulated, many were components of the universal cellular stress response, including DNA repair, redox regulation, protein folding, and proteolysis. Sodium:proton antiporters were among the processes significantly upregulated, indicating that these ion pumps were involved in maintaining cellular pH homeostasis. C. glacialis significantly alters its gene expression at low pH, although they maintain normal larval development. Understanding what confers tolerance to some species will support our ability to predict the effects of future ocean acidification on marine organisms.

Genomewide transcriptional reprogramming in the seagrass Cymodocea nodosa under experimental ocean acidification

Here, we report the first use of massive-scale RNA-sequencing to explore seagrass response to CO₂ -driven ocean acidification (OA). Large-scale gene expression changes in the seagrass Cymodocea nodosa occurred at CO₂ levels projected by the end of the century. C. nodosa transcriptome was obtained using Illumina RNA-Seq technology and de novo assembly, and differential gene expression was explored in plants exposed to short-term high CO₂ /low pH conditions. At high pCO₂ , there was a significant increased expression of transcripts associated with photosynthesis, including light reaction functions and CO₂ fixation, and also to respiratory pathways, specifically for enzymes involved in glycolysis, in the tricarboxylic acid cycle and in the energy metabolism of the mitochondrial electron transport. The upregulation of respiratory metabolism is probably supported by the increased availability of photosynthates and increased energy demand for biosynthesis and stress-related processes under elevated CO₂ and low pH. The upregulation of several chaperones resembling heat stress-induced changes in gene expression highlighted the positive role these proteins play in tolerance to intracellular acid stress in seagrasses. OA further modifies C. nodosa secondary metabolism inducing the transcription of enzymes related to biosynthesis of carbon-based secondary compounds, in particular the synthesis of polyphenols and isoprenoid compounds that have a variety of biological functions including plant defence. By demonstrating which physiological processes are most sensitive to OA, this research provides a major advance in the understanding of seagrass metabolism in the context of altered seawater chemistry from global climate change.

Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification

Ocean acidification (OA) is increasing due to anthropogenic CO₂ emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

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.

The development of contemporary European sea bass larvae (Dicentrarchus labrax) is not affected by projected ocean acidification scenarios

Ocean acidification is a recognized consequence of anthropogenic carbon dioxide (CO₂) emission in the atmosphere. Despite its threat to marine ecosystems, little is presently known about the capacity for fish to respond efficiently to this acidification. In adult fish, acid-base regulatory capacities are believed to be relatively competent to respond to hypercapnic conditions. However, fish in early life stage could be particularly sensitive to environmental factors as organs and important physiological functions become progressively operational during this period. In this study, the response of European sea bass (Dicentrarchus labrax) larvae reared under three ocean acidification scenarios, i.e., control (present condition, [Formula: see text] = 590 µatm, pH total = 7.9), low acidification (intermediate IPCC scenario, [Formula: see text] = 980 µatm, pH total = 7.7), and high acidification (most severe IPCC scenario, [Formula: see text] = 1520 µatm, pH total = 7.5) were compared across multiple levels of biological organizations. From 2 to 45 days-post-hatching, the chronic exposure to the different scenarios had limited influence on the survival and growth of the larvae (in the low acidification condition only) and had no apparent effect on the digestive developmental processes. The high acidification condition induced both faster mineralization and reduction in skeletal deformities. Global (microarray) and targeted (qPCR) analysis of transcript levels in whole larvae did not reveal any significant changes in gene expression across tested acidification conditions. Overall, this study suggests that contemporary sea bass larvae are already capable of coping with projected acidification conditions without having to mobilize specific defense mechanisms.

Differential protein expression using proteomics from a crustacean brine shrimp (Artemia sinica) under CO2-driven seawater acidification

Gradually increasing atmospheric CO₂ partial pressure (pCO₂) has caused an imbalance in carbonate chemistry and resulted in decreased seawater pH in marine ecosystems, termed seawater acidification. Anthropogenic seawater acidification is postulated to affect the physiology of many marine calcifying organisms. To understand the possible effects of seawater acidification on the proteomic responses of a marine crustacean brine shrimp (Artemia sinica) three groups of cysts were hatched and further raised in seawater at different pH levels (8.2 as control and 7.8 and 7.6 as acidification stress levels according to the predicted levels at the end of this century and next century, respectively) for 1, 7 and 14 days followed by examination of the protein expression changes via two-dimensional gel electrophoresis. Searches of protein databases revealed that 67 differential protein spots were altered due to lower pH level (7.6 and 7.8) stress in comparison to control groups (pH 8.2) by mass spectrometry. Generally, these differentially expressed proteins included the following: 1) metabolic process-related proteins involved in glycolysis and glucogenesis, nucleotide/amino acid/fatty acid metabolism, protein biosynthesis, DNA replication and apoptosis; 2) stress response-related proteins, such as peroxiredoxin, thioredoxin peroxidase, 70-kDa heat shock protein, Na/K ATPase, and ubiquinol-cytochrome c reductase; 3) immune defence-related proteins, such as prophenoloxidase and ferritin; 4) cytoskeletal-related proteins, such as myosin light chain, TCP1 subunit 2, tropomyosin and tubulin alpha chain; and 5) signal transduction-related proteins, such as phospholipase C-like protein, 14-3-3 zeta, translationally controlled tumour protein and RNA binding motif protein. Taken together, these data support the idea that CO₂-driven seawater acidification may affect protein expression in the crustacean A. sinica and possibly also in other species that feed on brine shrimp in the ecosystem, particularly marine food webs.

Biological responses to environmental heterogeneity under future ocean conditions

Organisms are projected to face unprecedented rates of change in future ocean conditions due to anthropogenic climate-change. At present, marine life encounters a wide range of environmental heterogeneity from natural fluctuations to mean climate change. Manipulation studies suggest that biota from more variable marine environments have more phenotypic plasticity to tolerate environmental heterogeneity. Here, we consider current strategies employed by a range of representative organisms across various habitats - from short-lived phytoplankton to long-lived corals - in response to environmental heterogeneity. We then discuss how, if and when organismal responses (acclimate/migrate/adapt) may be altered by shifts in the magnitude of the mean climate-change signal relative to that for natural fluctuations projected for coming decades. The findings from both novel climate-change modelling simulations and prior biological manipulation studies, in which natural fluctuations are superimposed on those of mean change, provide valuable insights into organismal responses to environmental heterogeneity. Manipulations reveal that different experimental outcomes are evident between climate-change treatments which include natural fluctuations vs. those which do not. Modelling simulations project that the magnitude of climate variability, along with mean climate change, will increase in coming decades, and hence environmental heterogeneity will increase, illustrating the need for more realistic biological manipulation experiments that include natural fluctuations. However, simulations also strongly suggest that the timescales over which the mean climate-change signature will become dominant, relative to natural fluctuations, will vary for individual properties, being most rapid for CO2 (~10 years from present day) to 4 decades for nutrients. We conclude that the strategies used by biota to respond to shifts in environmental heterogeneity may be complex, as they will have to physiologically straddle wide-ranging timescales in the alteration of ocean conditions, including the need to adapt to rapidly rising CO2 and also acclimate to environmental heterogeneity in more slowly changing properties such as warming.

De novo assembly of a tadpole shrimp (Triops newberryi) transcriptome and preliminary differential gene expression analysis

Next-generation sequencing techniques, such as RNA sequencing, have provided a wealth of genomic information for nonmodel species. Transcriptomic information can be used to quantify the patterns of gene expression, which can identify how environmental differences invoke organismal stress responses and provide a gauge in predicting species adaptability. In our study, we used RNA sequencing to characterize the first transcriptome from a naupliar tadpole shrimp (Triops newberryi) to identify the genes expressed during the early life history stages and which could be important for future genomic studies. RNA was extracted from naupliar T. newberryi that were reared in a laboratory-controlled setting and in two different water types, a native and a non-native condition. A total of six replicates, three per condition, were sequenced with the Illumina Hi-Seq 2000 achieving 365 M 50-nt reads. High-quality reads were produced and de novo assembly was used to construct a T. newberryi transcriptome that was approximately 24.8 M base pairs. More than 10 000 peptides were predicted from the assembly, and genes were sorted into gene ontology categories. The use of different water conditions allowed for a preliminary differential gene expression analysis in order to compare the changes in gene expression between conditions. There were 299 differentially expressed genes between water conditions that might serve as a focal point for future genomic studies of Triops acclimation to different environments. The Triops transcriptome could serve as vital genomic information for additional studies on Branchiopod crustaceans.

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.

Interactive Effects of Seawater Acidification and Elevated Temperature on the Transcriptome and Biomineralization in the Pearl Oyster Pinctada fucata

Interactive effects of ocean acidification and ocean warming on marine calcifiers vary among species, but little is known about the underlying mechanisms. The present study investigated the combined effects of seawater acidification and elevated temperature (ambient condition: pH 8.1 × 23 °C, stress conditions: pH 7.8 × 23 °C, pH 8.1 × 28 °C, and pH 7.8 × 28 °C, exposure time: two months) on the transcriptome and biomineralization of the pearl oyster Pinctada fucata, which is an important marine calcifier. Transcriptome analyses indicated that P. fucata implemented a compensatory acid-base mechanism, metabolic depression and positive physiological responses to mitigate the effects of seawater acidification alone. These responses were energy-expensive processes, leading to decreases in the net calcification rate, shell surface calcium and carbon content, and changes in the shell ultrastructure. Elevated temperature (28 °C) within the thermal window of P. fucata did not induce significant enrichment of the sequenced genes and conversely facilitated calcification, which was detected to alleviate the negative effects of seawater acidification on biomineralization and the shell ultrastructure. Overall, this study will help elucidate the mechanisms by which pearl oysters respond to changing seawater conditions and predict the effects of global climate change on pearl aquaculture.