Within-site variation of growth rates and terminal sizes in Mytilus californianus along wave exposure and tidal gradients

Life History Strategies Drive Meso-Scale Distribution Patterns in Coastal Benthic Macroinvertebrates

The environment shapes the spatial distribution of species, but species also comprise suites of traits which may indicate their adaptability to a specific environment. This forms the basis of trait biogeography studies. We thus examined how a species distribution is not only influenced by its environment and traits, but by interactions among its traits. Trait information was collected for 150 intertidal macroinvertebrates along a 3000 km environmental and biogeographic gradient on the South African coast. This information was analysed, as functional entities (FEs) were species performing similar functions that have the same trait values and were further condensed into two trait domains (Reproduction and Lifestyle). We then defined Life History Strategies (LHS) as specific combinations of Lifestyle and Reproduction FEs. Seven combinations of Lifestyle and Reproduction formed LHS that dominated total biomass. Some of these LHS were ubiquitous, while others showed geographic patterns across our west-east environmental gradient. For Lifestyle, filter-feeders exhibited high abundances on the East (subtropical, oligotrophic) and West (cool-temperate, eutrophic) extremes of the biogeographic gradient, but differed between the two in size at reproductive maturity and larval development type. This similarity in functionality of feeding mechanism and mobility with different reproductive strategies suggests a trait trade-off (investment in one trait reduces resources for others) between the Reproduction and Lifestyle domains. Within the Reproduction domain, gonochoristic, annual planktotrophic reproduction was common across bioregions, reflecting spin-offs (investment in one trait facilitates another trait) among these traits. Gonochoristic investment in less frequent episodic reproduction is another trade-off, with investment in large size and delayed maturation being a trade-off for many reproductive cycles. Overall, although our data supports the habitat templet model (i.e., the importance of environmental drivers), it further indicates that species distribution patterns observed along the South African coast reflect strong trait interactions and biomass patterns related to their LHS.

Digestive Enzyme Activities in Mussel Mytilus californianus Endure Acute Heat Exposure in Air

The mussel Mytilus californianus is an ecosystem engineer forming beds along the coastlines of Northeastern Pacific shores. As sessile organisms, they modulate their energy balance through valve movements, feeding, and digestive functionality. A recent study observed that activity of the digestive enzyme cellulase was higher than predicted in mussels high on the shore, where temperatures are characteristically high and food availability is limited compared to low-shore habitats. In the current study, we predicted that this scavenging behavior is induced to mitigate energy losses related to heat-shock responses-that cellulase and amylase will display hyperactivity for limited recourses in the face of aerial heating. In the laboratory, we acclimated mussels to three complex diets that differed in starch and cellulose composition, followed by two acute heat shocks (+8°C) in the laboratory. Results showed no hyperactivity of amylase and cellulase in heated mussels. These results differ from previous studies that showed lowered amylase activity following heat acclimation. This difference in amylase activity across heat-stress exposure time is important when analyzing mussel bed disturbances following heat waves that compromise energy balance or cause death within adult populations.

Interpreting life-history traits, seasonal cycles, and coastal climate from an intertidal mussel species: Insights from 9000 years of synthesized stable isotope data

Understanding past coastal variability is valuable for contextualizing modern changes in coastal settings, yet existing Holocene paleoceanographic records for the North American Pacific Coast commonly originate from offshore marine sediments and may not represent the dynamic coastal environment. A potential archive of eastern Pacific Coast environmental variability is the intertidal mussel species Mytilus californianus. Archaeologists have collected copious stable isotopic (δ18O and δ13C) data from M. californianus shells to study human history at California's Channel Islands. When analyzed together, these isotopic data provide windows into 9000 years of Holocene isotopic variability and M. californianus life history. Here we synthesize over 6000 δ18O and δ13C data points from 13 published studies to investigate M. californianus shell isotopic variability across ontogenetic, geographic, seasonal, and millennial scales. Our analyses show that M. californianus may grow and record environmental information more irregularly than expected due to the competing influences of calcification, ontogeny, metabolism, and habitat. Stable isotope profiles with five or more subsamples per shell recorded environmental information ranging from seasonal to millennial scales, depending on the number of shells analyzed and the resolution of isotopic subsampling. Individual shell profiles contained seasonal cycles and an accurate inferred annual temperature range of ~ 5°C, although ontogenetic growth reduction obscured seasonal signals as organisms aged. Collectively, the mussel shell record reflected millennial-scale climate variability and an overall 0.52‰ depletion in δ18Oshell from 8800 BP to the present. The archive also revealed local-scale oceanographic variability in the form of a warmer coastal mainland δ18Oshell signal (-0.32‰) compared to a cooler offshore islands δ18Oshell signal (0.33‰). While M. californianus is a promising coastal archive, we emphasize the need for high-resolution subsampling from multiple individuals to disentangle impacts of calcification, metabolism, ontogeny, and habitat and more accurately infer environmental and biological patterns recorded by an intertidal species.

Multi-omics reveals largely distinct transcript- and protein-level responses to the environment in an intertidal mussel

Organismal responses to stressful environments are influenced by numerous transcript- and protein-level mechanisms, and the relationships between expression changes at these levels are not always straightforward. Here, we used paired transcriptomic and proteomic datasets from two previous studies from gill of the California mussel, Mytilus californianus, to explore how simultaneous transcript and protein abundance patterns may diverge under different environmental scenarios. Field-acclimatized mussels were sampled from two disparate intertidal sites; individuals from one site were subjected to three further treatments (common garden, low-intertidal or high-intertidal outplant) that vary in temperature and feeding time. Assessing 1519 genes shared between the two datasets revealed that both transcript and protein expression patterns differentiated the treatments at a global level, despite numerous underlying discrepancies. There were far more instances of differential expression between treatments in transcript only (1451) or protein only (226) than of the two levels shifting expression concordantly (68 instances). Upregulated expression of cilium-associated transcripts (likely related to feeding) was associated with relatively benign field treatments. In the most stressful treatment, transcripts, but not proteins, for several molecular chaperones (including heat shock proteins and endoplasmic reticulum chaperones) were more abundant, consistent with a threshold model for induction of translation of constitutively available mRNAs. Overall, these results suggest that the relative importance of transcript- and protein-level regulation (translation and/or turnover) differs among cellular functions and across specific microhabitats or environmental contexts. Furthermore, the degree of concordance between transcript and protein expression can vary across benign versus acutely stressful environmental conditions.

Elevated aerial temperature modulates digestive enzyme activities in Mytilus californianus

The marine intertidal mussel Mytilus californianus aggregates to form beds along the Pacific shores of North America. As a sessile organism it must cope with fluctuations in temperature during low-tide aerial exposure, which elevates maintenance costs and negatively affects its overall energy budget. The function of its digestive gland is to release enzymes that break apart ingested polymers for subsequent nutrient absorption. The effects of elevated aerial warming acclimation on the functioning of digestive gland enzymes are not well studied. In this study we asked whether digestive gland carbohydases and proteases could be overstimulated in warm condition to possibly mitigate the costs related to the heat-shock response. We compared mussels acclimated to a + 9 °C heat-shock during daily low-tide aerial exposure to mussels acclimated to isothermal tidal conditions in a simulated intertidal system. The results showed fairly consistent activities of cellulase, trypsin, and amino-peptidase across tidal variation and between thermal treatments; however, amylase activity was lower in warmed versus cool mussels across low and high-tide. We also observed the expression of heat-shock genes in gill tissue during warm tidal conditions, suggestive that moderate temperatures during aerial exposure can induce a stress response.

Effects of Intertidal Position on Metabolism and Behavior in the Acorn Barnacle, Balanus glandula

The intertidal zone is characterized by persistent, tidally-driven fluctuations in both abiotic (e.g., temperature, oxygen, and salinity) and biotic (e.g., food availability and predation) factors, which make this a physiologically challenging habitat for resident organisms. The relative magnitude and degree of variability of environmental stress differ between intertidal zones, with the most extreme physiological stress often being experienced by organisms in the high intertidal. Given that so many of the constantly shifting parameters in this habitat are primary drivers of metabolic rate (e.g., temperature, [O₂], and food availability), we hypothesized that sessile conspecifics residing in different tidal zones would exhibit distinct “metabolic phenotypes,” a term we use to collectively describe the organisms’ baseline metabolic performance and capacity. To investigate this hypothesis, we collected acorn barnacles (Balanus glandula) from low, mid, and high intertidal positions in San Luis Obispo Bay, CA, and measured a suite of biochemical (whole-animal citrate synthase (CS) and lactate dehydrogenase (LDH) activity, and aerial [D-lactate]), physiological (O₂ consumption rates), morphological (body size), and behavioral (e.g., cirri beat frequency and percentage of time operculum open) indices of metabolism. We found tidal zone-dependent differences in B. glandula metabolism that primarily related to anaerobic capacity, cirral activity patterns, and body size. Barnacles from the low intertidal tended to have a greater capacity for anaerobic metabolism (i.e., increased LDH activity and increased baseline [D-lactate]), have reduced cirral beating activity—and presumably reduced feeding—when submerged, and be smaller in size compared to conspecifics in the high intertidal. We did not, however, see any D-lactate accumulation in barnacles from any tidal height throughout 96 h of air exposure. This trend indicates that the enhanced capacity of low intertidal barnacles for anaerobic metabolism may have evolved to support metabolism during more prolonged episodes of emersion or during events other than emersion (e.g., coastal hypoxia and predation). There were also no significant differences in CS activity or baseline O₂ consumption rates (in air or seawater at 14°C) across tidal heights, which imply that aerobic metabolic capacity may not be as sensitive to tidal position as anaerobic processes. Understanding how individuals occupying different shore heights differ in their metabolic capacity becomes increasingly interesting in the context of global climate change, given that the intertidal zone is predicted to experience even greater extremes in abiotic stress.

Inter-individual physiological variation in responses to environmental variation and environmental change: Integrating across traits and time

Greater understanding of physiological responses to climate change demands deeper comprehension of the causes and consequences of physiological variation. Increasingly, population trait means are being deconstructed into variable signals at the level of individuals. We advocate for greater consideration of such inter-individual physiological variation and how it both depends on and interacts with environmental variability. First, we review several studies on the intertidal mussel Mytilus californianus to illustrate how the magnitude of inter-individual variation may depend on the environmental context analyzed (i.e., is the mean condition benign or stressful?) and/or on the specific physiological metric investigated. Stressful conditions may reveal or mask variation in disparate ways at different levels of analysis (e.g., transcriptome vs. proteome), but we often lack crucial information regarding the relationships among these different physiological metrics and their consequences for fitness. We then reanalyze several published datasets to ask whether individuals employ divergent strategies over time in response to acute heat stress; such time-dependence would further complicate interpretation of physiological variation. However, definitive conclusions are precluded by limited sample sizes and short timescales in extant datasets. A key remaining challenge is to extend these analytical frameworks to longer periods over which individuals in a population experience repeated, but spatially variable, episodic stress events. We conclude that variation at multiple levels of analysis should be investigated over longer periods and, where possible, within individuals (or genotypes) experiencing repeated environmental challenges. Although difficult in practice, such studies will facilitate improved understanding of potential population-level physiological responses to climate change.

Responses of the mussel Brachidontes rodriguezii to aerial exposure: implications on growth and physiological condition

The length of aerial exposure (i.e., environmental conditions) of an organism, due to daily tides, induces physiological responses. A mark–recapture field experiment was conducted in two intertidal zones (low tide and high tide) using the stain calcein AM to determine growth, as measured by shell length, of the mussel Brachidontes rodriguezii (d’Orbigny, 1842) along its vertical distribution off the north Argentinean coast. In the high intertidal zone, B. rodriguezii exhibited slower growth in shell length because of the physiological stress resulting from aerial exposure during low tides. In the low intertidal zone and during spring, B. rodriguezii exhibited faster growth in shell length. It is suggested that growth of B. rodriguezii relies on several environmental factors (temperature, photoperiod, and aerial exposure) that fluctuate according to tidal height.

Thermal history and gape of individual Mytilus californianus correlate with oxidative damage and thermoprotective osmolytes

The ability of animals to cope with environmental stress depends – in part – on past experience, yet knowledge of the factors influencing an individual's physiology in nature remains underdeveloped. We used an individual monitoring system to record body temperature and valve gaping behavior of rocky intertidal zone mussels (Mytilus californianus). Thirty individuals were selected from two mussel beds (wave-exposed and wave-protected) that differ in thermal regime. Instrumented mussels were deployed at two intertidal heights (near the lower and upper edges of the mussel zone) and in a continuously submerged tidepool. Following a 23-day monitoring period, measures of oxidative damage to DNA and lipids, antioxidant capacities (catalase activity and peroxyl radical scavenging) and tissue contents of organic osmolytes were obtained from gill tissue of each individual. Univariate and multivariate analyses indicated that inter-individual variation in cumulative thermal stress is a predominant driver of physiological variation. Thermal history over the outplant period was positively correlated with oxidative DNA damage. Thermal history was also positively correlated with tissue contents of taurine, a thermoprotectant osmolyte, and with activity of the antioxidant enzyme catalase. Origin site differences, possibly indicative of developmental plasticity, were only significant for catalase activity. Gaping behavior was positively correlated with tissue contents of two osmolytes. Overall, these results are some of the first to clearly demonstrate relationships between inter-individual variation in recent experience in the field and inter-individual physiological variation, in this case within mussel beds. Such micro-scale, environmentally mediated physiological differences should be considered in attempts to forecast biological responses to a changing environment.

From classical to nonparametric growth models: Towards comprehensive modelling of mussel growth patterns

Understanding biological processes, such as growth, is crucial to development management and sustainability plans for bivalve populations. Von Bertalanffy and Gompertz models have been commonly used to fit bivalve growth. These models assume that individual growth is only determined by size, overlooking the effects of environmental and intrinsic conditions on growth patterns. The comparison between classical models and nonparametric GAM (generalized additive models) fits conducted in this work shows that the latter provide a more realistic approach of mussel growth measured in terms of shell length, and dry weight of hard and soft tissues. GAM fits detected a reduction in growth during the cold season, under unfavourable nutritional conditions. These fits also captured the decoupling between hard and soft tissue growth, widely addressed in the literature but not incorporated in growth models. In addition a GAM fit of condition index allowed us to explain annual changes in resources allocation, identifying the asymptotic growth of shell and the effects of the reproductive cycle on soft tissue fluctuations.

Modulation of digestive physiology and biochemistry in Mytilus californianus in response to feeding level acclimation and microhabitat

The intertidal mussel Mytilus californianus is a critical foundation species that is exposed to fluctuations in the environment along tidal- and wave-exposure gradients. We investigated feeding and digestion in mussels under laboratory conditions and across environmental gradients in the field. We assessed whether mussels adopt a rate-maximization (higher ingestion and lower assimilation) or a yield-maximization acquisition (lower ingestion and higher assimilation) strategy under laboratory conditions by measuring feeding physiology and digestive enzyme activities. We used digestive enzyme activity to define resource acquisition strategies in laboratory studies, then measured digestive enzyme activities in three microhabitats at the extreme ends of the tidal- and wave-exposure gradients within a stretch of shore (<20 m) projected sea-ward. Our laboratory results indicated that mussels benefit from a high assimilation efficiency when food concentration is low and have a low assimilation efficiency when food concentration is high. Additionally, enzyme activities of carbohydrases amylase, laminarinase and cellulase were elevated when food concentration was high. The protease trypsin, however, did not increase with increasing food concentration. In field conditions, low-shore mussels surprisingly did not have high enzyme activities. Rather, high-shore mussels exhibited higher cellulase activities than low-shore mussels. Similarly, trypsin activity in the high-shore-wave-sheltered microhabitat was higher than that in high-shore-wave-exposed. As expected, mussels experienced increasing thermal stress as a function of reduced submergence from low to high shore and shelter from wave-splash. Our findings suggest that mussels compensate for limited feeding opportunities and thermal stress by modulating digestive enzyme activities.