The energetic physiology of juvenile mussels, Mytilus chilensis (Hupe): The prevalent role of salinity under current and predicted pCO2 scenarios

Salinity-dependent top-down effect of rotifer Brachionus plicatilis on removing harmful alga Phaeocystis globosa

Using appropriate zooplankton to transfer the primary productivity of harmful algae to higher trophic levels through food chain is an eco-friendly mode to remove harmful algae. To assess the top-down efficiency of rotifer removing Phaeocystis and the salinity effect, we adopted a series of salinities to carry out Phaeocystis-rotifer population dynamics and rotifer life-history experiments. Results showed that the time for rotifers to remove Phaeocystis population was the shortest when the salinity was ≤20 ‰. With salinity rising to above 25 ‰, although the clearance time of Phaeocystis population by rotifer was significantly prolonged, ultimately the Phaeocystis population were almost completely eliminated at all salinities. Additionally, rotifer matured and reproduced earlier at low salinity, while high salinity significantly delayed first reproductive time and decreased the total offspring. The above findings are helpful to assess the impacts of external environmental factors on the application of zooplankton to control harmful algae.

Epigenetic variation mediated by lncRNAs accounts for adaptive genomic differentiation of the endemic blue mussel Mytiluschilensis

Epigenetic variation affects gene expression without altering the underlying DNA sequence of genes controlling ecologically relevant phenotypes through different mechanisms, one of which is long non-coding RNAs (lncRNAs). This study identified and evaluated the gene expression of lncRNAs in the gill and mantle tissues of Mytilus chilensis individuals from two ecologically different sites: Cochamó (41°S) and Yaldad (43°S), southern Chile, both impacted by climatic-related conditions and by mussel farming given their use as seedbeds. Sequences identified as lncRNAs exhibited tissue-specific differences, mapping to 3.54 % of the gill transcriptome and 1.96 % of the mantle transcriptome, representing an average of 2.76 % of the whole transcriptome. Using a high fold change value (≥|100|), we identified 43 and 47 differentially expressed lncRNAs (DE-lncRNAs) in the gill and mantle tissue of individuals sampled from Cochamó and 21 and 17 in the gill and mantle tissue of individuals sampled from Yaldad. Location-specific DE-lncRNAs were also detected in Cochamó (65) and Yaldad (94) samples. Via analysis of the differential expression of neighboring protein-coding genes, we identified enriched GO terms related to metabolic, genetic, and environmental information processing and immune system functions, reflecting how the impact of local ecological conditions may influence the M. chilensis (epi)genome expression. These DE-lncRNAs represent complementary biomarkers to DNA sequence variation for maintaining adaptive differences and phenotypic plasticity to cope with natural and human-driven perturbations.

Human-induced salinity changes impact marine organisms and ecosystems

Climate change is fundamentally altering marine and coastal ecosystems on a global scale. While the effects of ocean warming and acidification on ecology and ecosystem functions and services are being comprehensively researched, less attention is directed toward understanding the impacts of human-driven ocean salinity changes. The global water cycle operates through water fluxes expressed as precipitation, evaporation, and freshwater runoff from land. Changes to these in turn modulate ocean salinity and shape the marine and coastal environment by affecting ocean currents, stratification, oxygen saturation, and sea level rise. Besides the direct impact on ocean physical processes, salinity changes impact ocean biological functions with the ecophysiological consequences are being poorly understood. This is surprising as salinity changes may impact diversity, ecosystem and habitat structure loss, and community shifts including trophic cascades. Climate model future projections (of end of the century salinity changes) indicate magnitudes that lead to modification of open ocean plankton community structure and habitat suitability of coral reef communities. Such salinity changes are also capable of affecting the diversity and metabolic capacity of coastal microorganisms and impairing the photosynthetic capacity of (coastal and open ocean) phytoplankton, macroalgae, and seagrass, with downstream ramifications on global biogeochemical cycling. The scarcity of comprehensive salinity data in dynamic coastal regions warrants additional attention. Such datasets are crucial to quantify salinity-based ecosystem function relationships and project such changes that ultimately link into carbon sequestration and freshwater as well as food availability to human populations around the globe. It is critical to integrate vigorous high-quality salinity data with interacting key environmental parameters (e.g., temperature, nutrients, oxygen) for a comprehensive understanding of anthropogenically induced marine changes and its impact on human health and the global economy.

Combined effects of salinity and polystyrene microplastics exposure on the Pacific oysters Crassostrea gigas: Oxidative stress and energy metabolism

Microplastics (MPs) pollution and salinity variation are two environmental stressors, but their combined effects on marine mollusks are rarely known. Oysters (Crassostrea gigas) were exposed to 1 × 10⁴ particles L^-1 spherical polystyrene MPs (PS-MPs) of different sizes (small polystyrene MPs (SPS-MPs): 6 μm, large polystyrene MPs (LPS-MPs): 50-60 μm) under three salinity levels (21, 26, and 31 psu) for 14 days. Results demonstrated that low salinity reduced PS-MPs uptake in oysters. Antagonistic interactions between PS-MPs and low salinity mainly occurred, and partial synergistic effects were mainly induced by SPS-MPs. SPS-MPs induced higher lipid peroxidation (LPO) levels than LPS-MPs. In digestive glands, low salinity decreased LPO levels and glycometabolism-related gene expression, which was related to salinity levels. Low salinity instead of MPs mainly affected metabolomics profiles of gills through energy metabolism and osmotic adjustment pathway. In conclusion, oysters can adapt to combined stressors through energy and antioxidative regulation.

Survival and physiological energetics of highly invasive mussels exposed to heatwaves

The South American mussel, Mytella strigata, is a highly invasive fouling species of great concern along intertidal shores in East and Southeast Asia, posing serious threats to native biodiversity and ecosystems. Intertidal areas, being increasingly attacked by heatwaves over the last decade, are among the most thermally challenging habitats, yet the fate of this highly invasive mussel under scenarios of hotter heatwaves remains unknown. Here, we investigated how M. strigata responded to intensifying heatwaves frequently occurring in the South China Sea. Over 97% of individuals survived the five-day-lasting heatwaves, suggesting their high ability to cope with short-term heatwaves. Virtually unaffected clearance rate and absorption efficiency throughout the course of heatwaves indicate the maintenance of energy acquisition, and significantly decreased respiration rate implies the depression of energy metabolism, generating significant decreases in the O:N ratio when heatwaves occurred. Scope for growth of heatwaves-stressed mussels significantly decreased during initial exposure and then increased over time. These findings indicate the remarkable ability of M. strigata to cope with heatwaves recorded in its invasive habitats and call the attention for the rapid spread of this highly invasive fouling species in the context of climate change.

Chromosome-Level Genome Assembly of the Blue Mussel Mytilus chilensis Reveals Molecular Signatures Facing the Marine Environment

The blue mussel Mytilus chilensis is an endemic and key socioeconomic species inhabiting the southern coast of Chile. This bivalve species supports a booming aquaculture industry, which entirely relies on artificially collected seeds from natural beds that are translocated to diverse physical-chemical ocean farming conditions. Furthermore, mussel production is threatened by a broad range of microorganisms, pollution, and environmental stressors that eventually impact its survival and growth. Herein, understanding the genomic basis of the local adaption is pivotal to developing sustainable shellfish aquaculture. We present a high-quality reference genome of M. chilensis, which is the first chromosome-level genome for a Mytilidae member in South America. The assembled genome size was 1.93 Gb, with a contig N50 of 134 Mb. Through Hi-C proximity ligation, 11,868 contigs were clustered, ordered, and assembled into 14 chromosomes in congruence with the karyological evidence. The M. chilensis genome comprises 34,530 genes and 4795 non-coding RNAs. A total of 57% of the genome contains repetitive sequences with predominancy of LTR-retrotransposons and unknown elements. Comparative genome analysis of M. chilensis and M. coruscus was conducted, revealing genic rearrangements distributed into the whole genome. Notably, transposable Steamer-like elements associated with horizontal transmissible cancer were explored in reference genomes, suggesting putative relationships at the chromosome level in Bivalvia. Genome expression analysis was also conducted, showing putative genomic differences between two ecologically different mussel populations. The evidence suggests that local genome adaptation and physiological plasticity can be analyzed to develop sustainable mussel production. The genome of M. chilensis provides pivotal molecular knowledge for the Mytilus complex.

Shedding light on variation in reproductive success through studies of population genetic structure in a Southeast Pacific Coast mussel

Phylogeography often focuses on the spatial dimension of genetic diversity, rarely including the temporal dynamics occurring interannually among local populations, which can provide insight into past variations in reproductive success. Currently, there is an intense aquaculture industry of Mytilus spp. on the Southeast Pacific Coast which depends entirely on the spat released by natural populations forming a relevant and sensitive social-ecological system. Temporal and spatial spat variability from natural mussel beds could be related to interannual reproductive dynamics with variable reproductive success and recruitment, which leave genetic signatures. Temporal and spatial genetic structure was evaluated in six natural beds in the Southeast Pacific (from 39°25'S to 43°07'S) on the most abundant and widespread Mytilus lineage detected, Mytilus cf. chilensis, in 4 consecutive years. Analyses included data from >180 individuals per year, with a total of 751 (mitochondrial COI) and 747 (nuclear H1) individuals, respectively. Overall, both markers showed high haplotype diversity and low spatial and temporal genetic differentiation. Likely, the high dispersal capacity of Mytilus cf. chilensis maintains population homogeneity and prevents diversity erosion. The slight differences in genetic variance of COI were better explained by differences among sites (space), and conversely, the H1 genetic variance was better explained by interannual (temporal) comparisons, which could explain temporal variability in spat availability. This study highlights the important insights achieved with the evaluation of both temporal and spatial population genetic structures in marine species with high reproductive output, which can condition the success and sustainability of the relevant social-ecological system.

Metabolic rate allometry in intertidal mussels across environmental gradients: The role of coastal carbonate system parameters in mediating the effects of latitude and temperature

We assess the role of direct and indirect effects of coastal environmental drivers (including the parameters of the carbonate system) on energy expenditure (MR) and body mass (M) of the intertidal mussel, Perumytilus purpuratus, across 10 populations distributed over 2800 km along the Southern Eastern Pacific (SEP) coast. We find biogeographic and local variation in carbonate system variables mediates the effects of latitude and temperature on metabolic rate allometry along the SEP coast. Also, the fitted Piecewise Structural Equation models (PSEM) have greater predictive ability (conditional R² = 0.95) relative to the allometric scaling model (R² = 0.35). The largest standardized coefficients for MR and M were determined by the influence of temperature and latitude, followed by pCO₂, pH, total alkalinity, and salinity. Thus, physiological diversity of P. purpuratus along the SEP coast emerges as the result of direct and indirect effects of biogeographic and local environmental variables.

Ocean acidification alters anti-predator responses in a competitive dominant intertidal mussel

Widespread intertidal mussels are exposed to a variety of natural and anthropogenic stressors. Even so, our understanding of the combined influence of stressors such as predation risk and ocean acidification (OA) on these species remains limited. This study examined the response of the purple mussel (Perumytilus purpuratus), a species distributed along Pacific southeastern rocky shores, to the effects of predation risk and OA. Using a laboratory 2 × 2 cross design, purple mussels were either devoid or exposed to predator cues from the muricid snail Acanthina monodon, while simultaneously exposing them to current (500 ppm) or projected OA conditions (1500 ppm). The response of purple mussels to these factors was assessed using growth, calcification, clearance, and metabolic rates, in addition to byssus production. After 60 d, the presence of predator cues reduced mussel growth in width and length, and in the latter case, OA enhanced this response making the effects of predator cues more severe. Calcification rates were driven by the interaction between the two stressors, whereas clearance rates increased only in response to OA, likely explaining some of the growth results. Mussel byssus production also increased with pCO₂ but interacted with predation risk: in the absence of predator cues, byssus production increased with OA. These results suggest that projected levels of OA may alter and in some cases prevail over the natural response of purple mussels to predation risk. Considering the role played by this mussel as a dominant competitor and ecosystem engineer in rocky shores, these results have community-wide implications.

Morphological, physiological and behavioral responses of an intertidal snail, Acanthina monodon (Pallas), to projected ocean acidification and cooling water conditions in upwelling ecosystems

Ocean acidification (OA) is expected to rise towards the end of the 21st century altering the life history traits in marine organisms. Upwelling systems will not escape OA, but unlike other areas of the ocean, cooling effects are expected to intensify in these systems. Regardless, studies evaluating the combined effects of OA and cooling remain scarce. We addressed this gap using a mesocosm system, where we exposed juveniles of the intertidal muricid snail Acanthina monodon to current and projected pCO₂ (500 vs. 1500 ppm) and temperature (15 vs. 10 °C) from the southeast Pacific upwelling system. After 9 weeks of experimental exposure to those conditions, we conducted three estimations of growth (wet weight, shell length and shell peristomal length), in addition to measuring calcification, metabolic and feeding rates and the ability of these organisms to return to the normal upright position after being overturned (self-righting). Growth, feeding and calcification rates increased in projected cooling conditions (10 °C) but were unaffected by pCO₂ or the interaction between pCO₂ and temperature. Instead, metabolic rates were driven by pCO₂, but a significant interaction with temperature suggests that in cooler conditions, metabolic rates will increase when associated with high pCO₂ levels. Snail self-righting times were not affected across treatments. These results suggest that colder temperatures projected for this area would drive this species growth, feeding and calcification, and consequently, some of its population biology and productivity. However, the snails may need to compensate for the increase in metabolic rates under the effects of ocean acidification. Although A. monodon ability to adjust to individual or combined stressors will likely account for some of the changes described here, our results point to a complex dynamic to take place in intertidal habitats associated with upwelling systems.

Plasticity in organic composition maintains biomechanical performance in shells of juvenile scallops exposed to altered temperature and pH conditions

The exposure to environmental variations in pH and temperature has proven impacts on benthic ectotherms calcifiers, as evidenced by tradeoffs between physiological processes. However, how these stressors affect structure and functionality of mollusk shells has received less attention. Episodic events of upwelling of deep cold and low pH waters are well documented in eastern boundary systems and may be stressful to mollusks, impairing both physiological and biomechanical performance. These events are projected to become more intense, and extensive in time with ongoing global warming. In this study, we evaluate the independent and interactive effects of temperature and pH on the biomineral and biomechanical properties of Argopecten purpuratus scallop shells. Total organic matter in the shell mineral increased under reduced pH (~ 7.7) and control conditions (pH ~ 8.0). The periostracum layer coating the outer shell surface showed increased protein content under low pH conditions but decreasing sulfate and polysaccharides content. Reduced pH negatively impacts shell density and increases the disorder in the orientation of calcite crystals. At elevated temperatures (18 °C), shell microhardness increased. Other biomechanical properties were not affected by pH/temperature treatments. Thus, under a reduction of 0.3 pH units and low temperature, the response of A. purpuratus was a tradeoff among organic compounds (biopolymer plasticity), density, and crystal organization (mineral plasticity) to maintain shell biomechanical performance, while increased temperature ameliorated the impacts on shell hardness. Biopolymer plasticity was associated with ecophysiological performance, indicating that, under the influence of natural fluctuations in pH and temperature, energetic constraints might be critical in modulating the long-term sustainability of this compensatory mechanism.

Variable Environments in an Upwelling System Trigger Differential Thermal Sensitivity in a Low Intertidal Chiton

Environmental variability in coastal oceans associated with upwelling dynamics probably is one of the most pervasive forces affecting the physiological performance of marine life. As the environmental temperature is the abiotic factor with major incidence in the physiology and ecology of marine ectotherms, the abrupt temperature changes in upwelling systems could generate important variations in these organisms’ functional processes. The relationship between ambient temperature and physiological performance can be described through a thermal performance curve (TPC). The parameters of this curve usually show geographic variation usually is in accordance with the predictions of the climate variability hypothesis (CVH), which states that organisms inhabiting more variable environments should have broader ranges of environmental tolerance in order to cope with the fluctuating environmental conditions they experience. Here we study the effect generated by the environmental variability in an active upwelling zone on the physiological performance of the marine ectotherm Achanthopleura echinata. In particular, we compared the parameters of the TPC and the metabolic rate of two populations of A. echinata, one found in high semi-permanent upwelling (Talcaruca), while the other is situated in an adjacent area with seasonal upwelling (Los Molles) and therefore more stable environmental conditions. Our results show that: (1) oxygen consumption increases with body size and this effect is more significant in individuals from the Talcaruca population, (2) optimal temperature, thermal breadth, upper critical limit and maximum performance were higher in the population located in the area of high environmental heterogeneity and (3) individuals from Talcaruca showed greater variance in optimal temperature, thermal breadth, upper critical limit but not in maximum performance. Although it is clear that a variable environment affects the thermal physiology of organisms, expanding their tolerance ranges and generating energy costs in the performance of individuals, it is relevant to note that upwelling systems are multifactorial phenomena where the rise of water masses modifies not only temperature, but also decreases O2, pH, and increases pCO2 which in turn could modify metabolism and TPC.

Adaptive Differences in Gene Expression in Farm-Impacted Seedbeds of the Native Blue Mussel Mytilus chilensis

The study of adaptive population differences is relevant for evolutionary biology, as it evidences the power of selective local forces relative to gene flow in maintaining adaptive phenotypes and their underlying genetic determinants. However, human-mediated hybridization through habitat translocations, a common and recurrent aquaculture practice where hybrids could eventually replace local genotypes, risk populations' ability to cope with perturbations. The endemic marine mussel Mytilus chilensis supports a booming farming industry in the inner sea of Chiloé Island, southern Chile, which entirely relies on artificially collected seeds from natural beds that are translocated to ecologically different fattening centers. A matter of concern is how farm-impacted seedbeds will potentially cope with environmental shifts and anthropogenic perturbations. This study provides the first de novo transcriptome of M. chilensis; assembled from tissue samples of mantles and gills of individuals collected in ecologically different farm-impacted seedbeds, Cochamó (41°S) and Yaldad (43°S). Both locations and tissue samples differentially expressed transcripts (DETs) in candidate adaptive genes controlling multiple fitness traits, involved with metabolism, genetic and environmental information processing, and cellular processes. From 189,743 consensus contigs assembled: 1,716 (Bonferroni p value ≤ 0.05) were DETs detected in different tissues of samples from different locations, 210 of them (fold change ≥ | 100|) in the same tissue of samples from a different location, and 665 (fold change ≥ | 4|) regardless of the tissue in samples from a different location. Site-specific DETs in Cochamó (169) and Yaldad (150) in candidate genes controlling tolerance to temperature and salinity shifts, and biomineralization exhibit a high number of nucleotide genetic variants with regular occurrence (frequency > 99%). This novel M. chilensis transcriptome should help assessing and monitoring the impact of translocations in wild and farm-impacted mussel beds in Chiloé Island. At the same time, it would help designing effective managing practices for conservation, and translocation traceability.

Environmental refuges increase performance of juvenile mussels Mytilus chilensis: Implications for mussel seedling and farming strategies

Estuarine ecosystems are characterized by a wide physical-chemical variation that in the context of global change scenarios may be exacerbated in the future. The fitness of resident organisms is expected to be influenced by such variation and, hence, its study is a priority. Some of that variation relates to water vertical stratification, which may create "environmental refuges" or distinct layers of water with conditions favoring the fitness of some individuals and species. This study explored the performance of juvenile mussels (M. chilensis) settled in two distinctive water depths (1 m and 4 m) of the Reloncaví fjord (southern Chile) by conducting a reciprocal transplants experiment. Salinity, saturation state and the contents of CO₃ in seawater were among the factors that best explained the differences between the two layers. In such environmental conditions, the mussel traits that responded to such variation were growth and calcification rates, with significantly higher values at 4 m deep, whereas the opposite, increased metabolic stress, was higher in mussels raised and transplanted to the surface waters (1 m). Such differences support the notion of an environmental refuge, where species like mussels can find better growth conditions and achieve higher performance levels. These results are relevant considering the importance of M. chilensis as a shellfish resource for aquaculture and a habitat forming species. In addition, these results shed light on the variable responses exhibited by estuarine organisms to small-scale changes in the characteristics of the water column, which in turn will help to better understand the responses of the organisms to the projected scenarios of climate global change.

The combined effects of salinity and pH on shell biomineralization of the edible mussel Mytilus chilensis

Coastal ecosystems influenced by river discharges are subjected to important environmental changes. Understanding how marine biota cope with its environment is relevant in predicting the responses to future conditions imposed by climate change. To date, a large number of studies have addressed the role of pH on shell and biomineralization properties on multiple calcifying species; however the role of salinity in combination with other stressors has been poorly studied. In particular, the edible mussel Mytilus chilensis, an important marine resource of the Chilean coasts, inhabits estuarine areas which show high natural variability in terms of pH and salinity. Here, we studied how M. chilensis shell periostracum, shell organic matrix and crystal orientation are affected by different pH (8.1 and 7.7) and salinity conditions (30, 25 and 20 psu), isolated and in combination, at different time intervals. Our results show differences in the plasticity of the different biomineralogical properties studied during the experiment under the different pH and salinity treatments. While the periostracum thickness and the total shell organic matter were not affected by pH and salinity, the periostracum organic composition did. Higher amounts of polysaccharides were observed under low pH conditions after 20 days of experiment, while after 60 days, low salinity was responsible for the decrease of the polysaccharides and proteins in the periostracum. Low salinity also produced a major disorder in crystal organization at the outer shell surface. Finally, total shell weight was only affected by low pH conditions under lower salinity conditions (20 psu). From the results, in the majority of the shell properties observed we did not observe any combined effect of pH and salinity. Also, we detected that the magnitude of the impacts of salinity and pH are variable and time-dependent. This would be suggesting some level of acclimatization of M. chilensis to lower pH and salinity conditions.

Plastic response of the oyster Ostrea chilensis to temperature and pCO2 within the present natural range of variability

Estuaries are characterized by high fluctuation of their environmental conditions. Environmental parameters measured show that the seawater properties of the Quempillén estuary (i.e. temperature, salinity, pCO₂, pH and ΩCaCO₃) were highly fluctuating and related with season and tide. We test the effects of increasing temperature and pCO₂ in the seawater on the physiological energetics of the bivalve Ostrea chilensis. Juvenile oysters were exposed to an orthogonal combination of three temperatures (10, 15, and 20°C) and two pCO₂ levels (~400 and ~1000 μatm) for a period of 60 days to evaluate the temporal effect (i.e. 10, 20, 30, 60 days) on the physiological rates of the oysters. Results indicated a significant effect of temperature and time of exposure on the clearance rate, while pCO₂ and the interaction between pCO₂ and the other factors studied did not show significant effects. Significant effects of temperature and time of exposure were also observed on the absorption rate, but not the pCO₂ nor its interaction with other factors studied. Oxygen consumption was significantly affected by pCO₂, temperature and time. Scope for growth was only significantly affected by time; despite this, the highest values were observed for individuals subject to to 20°C and to ~1000 μatm pCO₂. In this study, Ostrea chilensis showed high phenotypic plasticity to respond to the high levels of temperature and pCO₂ experienced in its habitat as no negative physiological effects were observed. Thus, the highly variable conditions of this organism’s environment could select for individuals that are more resistant to future scenarios of climate change, mainly to warming and acidification.

Assessing the impact of elevated pCO2 within and across generations in a highly invasive fouling mussel (Musculista senhousia)

Marine biofouling by the swiftly spreading invasive mussel (Musculista senhousia) has caused serious ecological and economic consequences in the global coastal waters. However, the fate of this highly invasive fouling species in a rapidly acidifying ocean remains unknown. Here, we demonstrated the impacts of ocean acidification within and across generations, to understand whether M. senhousia has the capacity to acclimate to changing ocean conditions. During the gonadal development, exposure of mussels to elevated pCO₂ caused significant decreases of survival, growth performance and condition index, and shifted the whole-organism energy budget by inflating energy expenses to fuel compensatory processes, eventually impairing the success of spawning. Yet, rapid transgenerational acclimation occurred during the early life history stage and persisted into adulthood. Eggs spawned from CO₂-exposed mussels were significantly bigger compared with those from non-CO₂-exposed mussels, indicating increased maternal provisioning into eggs and hence conferring larvae resilience under harsh conditions. Larvae with a prior history of transgenerational exposure to elevated pCO₂ developed faster and had a higher survival than those with no prior history of CO₂ exposure. Transgenerational exposure significantly increased the number of larvae completing metamorphosis. While significant differences in shell growth were no longer observed during juvenile nursery and adult grow-out, transgenerationally exposed mussels displayed improved survival in comparison to non-transgenerationally exposed mussels. Metabolic plasticity arose following transgenerational acclimation, generating more energy available for fitness-related functions. Overall, the present study demonstrates the remarkable ability of M. senhousia to respond plastically and acclimate rapidly to changing ocean conditions.