Proteomic analysis was performed to compare the effects of Arsenic (As), seawater acidification (Low pH) and the combination of both stressors (Low pH + As) on Crassostrea angulata and Crassostrea gigas juveniles in the context of global environmental change. This study aimed to elucidate if two closely related Crassostrea species respond similarly to these environmental stressors, considering both single and combined exposures, to infer if the simultaneous exposure to both stressors induced a differentiated response. Identification of the most important differentially expressed proteins between conditions revealed marked differences in the response of each species towards single and combined exposures, evidencing species-related differences towards each experimental condition. Moreover, protein alterations observed in the combined exposure (Low pH + As) were substantially different from those observed in single exposures. Identified proteins and their putative biological functions revealed an array of modes of action in each condition. Among the most important, those involved in cellular structure (Actin, Atlastin, Severin, Gelsolin, Coronin) and extracellular matrix modulation (Ependymin, Tight junction ZO-1, Neprilysin) were strongly regulated, although in different exposure conditions and species. Data also revealed differences regarding metabolic modulation capacity (ATP beta, Enolase, Aconitate hydratase) and oxidative stress response (Aldehyde dehydrogenase, Lactoylglutathione, Retinal dehydrogenase) of the species, which also depended on single or combined exposures, illustrating a different response capacity of both oyster species to the presence of multiple stressors. Interestingly, alterations of piRNA abundance in C. angulata suggested genome reconfiguration in response to multiple stressors, likely an important mode of action related to adaptive evolution mechanisms previously unknown to oyster species, which requires further investigation. Our findings provide a deeper insight into the complexity of C. angulata and C. gigas responses to environmental stress at the proteome level, evidencing different capacities to endure abiotic changes, with relevance regarding the ecophysiological fitness of the species and competitive advantages in a changing environment.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2021) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2022-04-11.