A Dynamic Energy Budget (DEB) model for the keystone predator Pisaster ochraceus

Resource and seasonality drive interspecific variability in simulations from a dynamic energy budget model

Animals show a vast array of phenotypic traits in time and space. Such variation patterns have traditionally been described as ecogeographical rules; for example, the tendency of size and clutch size to increase with latitude (Bergmann's and Lack's rules, respectively). Despite considerable research into these variation patterns and their consequences for biodiversity and conservation, the processes behind trait variation remain controversial. Here, we show how food variability, largely set by climate and weather, can drive interspecific trait variation by determining individual energy input and allocation trade-offs. Using a dynamic energy budget (DEB) model, we simulated different food environments, as well as interspecific variability in the parameters for energy assimilation, mobilization and allocation to soma. We found that interspecific variability is greater when the resource is non-limiting in both constant and seasonal environments. Our findings further show that individuals can reach larger biomass and greater reproductive output in a seasonal environment than in a constant environment of equal average resource due to the peaks of food surplus. Our results agree with the classical patterns of interspecific trait variation and provide a mechanistic understanding that supports recent hypotheses which explain them: the resource and the eNPP (net primary production during the growing season) rules. Due to the current alterations to ecosystems and communities, disentangling the mechanisms of trait variation is increasingly important to understand and predict biodiversity dynamics under climate change, as well as to improve conservation strategies.

The role of Dynamic Energy Budgets in conservation physiology

The contribution of knowledge, concepts and perspectives from physiological ecology to conservation decision-making has become critical for understanding and acting upon threats to the persistence of sensitive species. Here we review applications of dynamic energy budget (DEB) theory to conservation issues and discuss how this theory for metabolic organization of all life on earth (from bacteria to whales) is well equipped to support current and future investigations in conservation research. DEB theory was first invented in 1979 in an applied institution for environmental quality assessment and mitigation. The theory has since undergone extensive development and applications. An increasing number of studies using DEB modelling have provided valuable insights and predictions in areas that pertain to conservation such as species distribution, evolutionary biology, toxicological impacts and ecosystem management. We discuss why DEB theory, through its mechanistic nature, its universality and the wide range of outcomes it can provide represents a valuable tool to tackle some of the current and future challenges linked to maintaining biodiversity, ensuring species survival, ecotoxicology, setting water and soil quality standards and restoring ecosystem structure and functioning in a changing environment under the pressure of anthropogenic driven changes.

Effects of triclosan exposure on the energy budget of Ruditapes philippinarum and R. decussatus under climate change scenarios

We built a simulation model based on Dynamic Energy Budget theory (DEB) to assess the growth and reproductive potential of the native European clam Ruditapes decussatus and the introduced Manila clam Ruditapes philippinarum under current temperature and pH conditions in a Portuguese estuary and under those forecasted for the end of the 21st c. The climate change scenario RCP8.5 predicts temperature increase of 3 °C and a pH decrease of 0.4 units. The model was run under additional conditions of exposure to the emerging contaminant triclosan (TCS) and in the absence of this compound. The parameters of the DEB model were calibrated with the results of laboratory experiments complemented with data from the literature available for these two important commercial shellfish resources. For each species and experimental condition (eight combinations), we used data from the experiments to produce estimates for the key parameters controlling food intake flux, assimilation flux, somatic maintenance flux and energy at the initial simulation time. The results showed that the growth and reproductive potential of both species would be compromised under future climate conditions, but the effect of TCS exposure had a higher impact on the energy budget than forecasted temperature and pH variations. The egg production of R. philippinarum was projected to suffer a more marked reduction with exposure to TCS, regardless of the climatic factor, while the native R. decussatus appeared more resilient to environmental causes of stress. The results suggest a likely decrease in the rates of expansion of the introduced R. philippinarum in European waters, and negative effects on fisheries and aquaculture production of exposure to emerging contaminants (e.g., TCS) and climate change.

Dynamic Energy Budget model suggests feeding constraints and physiological stress in black-lip pearl oysters, 5 years post mass-mortality event

Mass-mortality events of marine species can disturb the structure of communities. While identifying the causes of mass-mortality events is crucial for implementing recovery strategies, monitoring is challenging in remote locations. Black-lip pearl oysters (Pinctada margaritifera) are farmed for producing black pearls within remote atolls of French Polynesia. Previous mass-mortality events have resulted in the collapse of oysters and other species; however, the causes and conditions that favour recovery are unclear. We investigated the potential for oyster population recovery 5 years after a mortality event at Takaroa Atoll (Tuamotu Archipelago). Temperature, food availability (total chlorophyll-a), growth and reproduction were monitored. Growth was also simulated using a Dynamic Energy Budget model. Despite favourable conditions, reduced growth and reproduction signalled an energetic deficit. The model overpredicted growth, and supported the hypotheses that individuals are unable to profit from the phytoplankton available and maintenance costs are high in Takaroa, ultimately explaining their poor physiological condition.

Seascape connectivity of European anchovy in the Central Mediterranean Sea revealed by weighted Lagrangian backtracking and bio-energetic modelling

Ecological connectivity is one of the most important processes that shape marine populations and ecosystems, determining their distribution, persistence, and productivity. Here we use the synergy of Lagrangian back-trajectories, otolith-derived ages of larvae, and satellite-based chlorophyll-a to identify spawning areas of European anchovy from ichthyoplanktonic data, collected in the Strait of Sicily (Central Mediterranean Sea), i.e., the crucial channel in between the European and African continents. We obtain new evidence of ecosystem connectivity between North Africa and recruitment regions off the southern European coasts. We assess this result by using bio-energetic modeling, which predicts species-specific responses to environmental changes by producing quantitative information on functional traits. Our work gives support to a collaborative and harmonized use of Geographical Sub-Areas, currently identified by the General Fisheries Commission for the Mediterranean. It also confirms the need to incorporate climate and environmental variability effects into future marine resources management plans, strategies, and directives.

Range-extending coral reef fishes trade-off growth for maintenance of body condition in cooler waters

As ocean waters warm due to climate change, tropical species are shifting their ranges poleward to remain within their preferred thermal niches. As a result, novel communities are emerging in which tropical species interact with local temperate species, competing for similar resources, such as food and habitat. To understand how range-extending coral reef fish species perform along their leading edges when invading temperate ecosystems, we studied proxies of their fitness, including somatic growth (length increase), feeding rates, and body condition, along a 730-km latitudinal gradient situated in one of the global warming hotspots. We also studied co-occurring temperate species to assess how their fitness is affected along their trailing edges under ocean warming. We predicted that tropical fishes would experience reduced performance as they enter novel communities with suboptimal environmental conditions. Our study shows that although tropical fish maintain their body condition (based on three proxies) and stomach fullness across all invaded temperate latitudes, they exhibit decreased in situ growth rates, activity levels, and feeding rates in their novel temperate environment, likely a result of lower metabolic rates in cooler waters. We posit that tropical fishes face a growth-maintenance trade-off under the initial phases of ocean warming (i.e. at their leading edges), allowing them to maintain their body condition in cooler temperate waters but at the cost of slower growth. Temperate fish exhibited no distinct patterns in body condition and performance along the natural temperature gradient studied. However, in the face of future climate change, when metabolism is no longer stymied by low water temperatures, tropical range-extending species are likely to approach their native-range growth rates along their leading edges, ultimately leading to increased competitive interactions with local species in temperate ecosystems.

Disease epidemic and a marine heat wave are associated with the continental-scale collapse of a pivotal predator (Pycnopodia helianthoides)

Multihost infectious disease outbreaks have endangered wildlife, causing extinction of frogs and endemic birds, and widespread declines of bats, corals, and abalone. Since 2013, a sea star wasting disease has affected >20 sea star species from Mexico to Alaska. The common, predatory sunflower star (Pycnopodia helianthoides), shown to be highly susceptible to sea star wasting disease, has been extirpated across most of its range. Diver surveys conducted in shallow nearshore waters (n = 10,956; 2006-2017) from California to Alaska and deep offshore (55 to 1280 m) trawl surveys from California to Washington (n = 8968; 2004-2016) reveal 80 to 100% declines across a ~3000-km range. Furthermore, timing of peak declines in nearshore waters coincided with anomalously warm sea surface temperatures. The rapid, widespread decline of this pivotal subtidal predator threatens its persistence and may have large ecosystem-level consequences.

Applicability of Dynamic Energy Budget (DEB) models across steep environmental gradients

Robust ecological forecasting requires accurate predictions of physiological responses to environmental drivers. Energy budget models facilitate this by mechanistically linking biology to abiotic drivers, but are usually ground-truthed under relatively stable physical conditions, omitting temporal/spatial environmental variability. Dynamic Energy Budget (DEB) theory is a powerful framework capable of linking individual fitness to environmental drivers and we tested its ability to accommodate variability by examining model predictions across the rocky shore, a steep ecotone characterized by wide fluctuations in temperature and food availability. We parameterized DEB models for co-existing mid/high-shore (Mytilus galloprovincialis) and mid/low-shore (Perna perna) mussels on the south coast of South Africa. First, we assumed permanently submerged conditions, and then incorporated metabolic depression under low tide conditions, using detailed data of tidal cycles, body temperature and variability in food over 12 months at three sites. Models provided good estimates of shell length for both species across the shore, but predictions of gonadosomatic index were consistently lower than observed. Model disagreement could reflect the effects of details of biology and/or difficulties in capturing environmental variability, emphasising the need to incorporate both. Our approach provides guidelines for incorporating environmental variability and long-term change into mechanistic models to improve ecological predictions.

Decimation by sea star wasting disease and rapid genetic change in a keystone species, Pisaster ochraceus

Standing genetic variation enables or restricts a population's capacity to respond to changing conditions, including the extreme disturbances expected to increase in frequency and intensity with continuing anthropogenic climate change. However, we know little about how populations might respond to extreme events with rapid genetic shifts, or how population dynamics may influence and be influenced by population genomic change. We use a range-wide epizootic, sea star wasting disease, that onset in mid-2013 and caused mass mortality in Pisaster ochraceus to explore how a keystone marine species responded to an extreme perturbation. We integrated field surveys with restriction site-associated DNA sequencing data to (i) describe the population dynamics of mortality and recovery, and (ii) compare allele frequencies in mature P. ochraceus before the disease outbreak with allele frequencies in adults and new juveniles after the outbreak, to identify whether selection may have occurred. We found P. ochraceus suffered 81% mortality in the study region between 2012 and 2015, and experienced a concurrent 74-fold increase in recruitment beginning in late 2013. Comparison of pre- and postoutbreak adults revealed significant allele frequency changes at three loci, which showed consistent changes across the large majority of locations. Allele frequency shifts in juvenile P. ochraceus (spawned from premortality adults) were consistent with those seen in adult survivors. Such parallel shifts suggest detectable signals of selection and highlight the potential for persistence of this change in subsequent generations, which may influence the resilience of this keystone species to future outbreaks.

The evolution of gonad expenditure and gonadosomatic index (GSI) in male and female broadcast-spawning invertebrates

Sedentary broadcast-spawning marine invertebrates, which release both eggs and sperm into the water for fertilization, are of special interest for sexual selection studies. They provide unique insight into the early stages of the evolutionary succession leading to the often-intense operation of both pre- and post-mating sexual selection in mobile gonochorists. Since they are sessile or only weakly mobile, adults can interact only to a limited extent with other adults and with their own fertilized offspring. They are consequently subject mainly to selection on gamete production and gamete success, and so high gonad expenditure is expected in both sexes. We review literature on gonadosomatic index (GSI; the proportion of body tissue devoted to gamete production) of gonochoristic broadcast spawners, which we use as a proxy for gonad expenditure. We show that such taxa most often have a high GSI that is approximately equal in both sexes. When GSI is asymmetric, female GSI usually exceeds male GSI, at least in echinoderms (the majority of species recorded). Intriguingly, though, higher male GSI also occurs in some species and appears more common than female-biased GSI in certain orders of gastropod molluscs. Our limited data also suggest that higher male GSI may be the prevalent pattern in sperm casters (where only males release gametes). We explore how selection might have shaped these patterns using game theoretic models for gonad expenditure that consider possible trade-offs with (i) somatic maintenance or (ii) growth, while also considering sperm competition, sperm limitation, and polyspermy. Our models of the trade-off between somatic tissue (which increases survival) and gonad (which increases reproductive success) predict that GSI should be equal for the two sexes when sperm competition is intense, as is probably common in broadcast spawners due to synchronous spawning in aggregations. Higher female GSI occurs under low sperm competition. Sperm limitation appears unlikely to alter these conclusions qualitatively, but can also act as a force to keep male GSI high, and close to that of females. Polyspermy can act to reduce male GSI. Higher male than female GSI is predicted to be less common (as observed in the data), but can occur when ova/ovaries are sufficiently more resource-intensive to produce than sperm/testes, for which some evidence exists. We also show that sex-specific trade-offs between gonads and growth can generate different life-history strategies for males and females, with males beginning reproduction earlier. This could lead to apparently higher male GSI in empirical studies if immature females are included in calculations of mean GSI. The existence of higher male GSI nonetheless remains somewhat problematic and requires further investigation. When sperm limitation is low, we suggest that the natural logarithm of the male/female GSI ratio may be a suitable index for sperm competition level in broadcast spawners, and that this may also be considered as an index for internally fertilizing taxa.

Physics of metabolic organization

We review the most comprehensive metabolic theory of life existing to date. A special focus is given to the thermodynamic roots of this theory and to implications that the laws of physics-such as the conservation of mass and energy-have on all life. Both the theoretical foundations and biological applications are covered. Hitherto, the foundations were more accessible to physicists or mathematicians, and the applications to biologists, causing a dichotomy in what always should have been a single body of work. To bridge the gap between the two aspects of the same theory, we (i) adhere to the theoretical formalism, (ii) try to minimize the amount of information that a reader needs to process, but also (iii) invoke examples from biology to motivate the introduction of new concepts and to justify the assumptions made, and (iv) show how the careful formalism of the general theory enables modular, self-consistent extensions that capture important features of the species and the problem in question. Perhaps the most difficult among the introduced concepts, the utilization (or mobilization) energy flow, is given particular attention in the form of an original and considerably simplified derivation. Specific examples illustrate a range of possible applications-from energy budgets of individual organisms, to population dynamics, to ecotoxicology.

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

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

Ochre star mortality during the 2014 wasting disease epizootic: role of population size structure and temperature

Over 20 species of asteroids were devastated by a sea star wasting disease (SSWD) epizootic, linked to a densovirus, from Mexico to Alaska in 2013 and 2014. For Pisaster ochraceus from the San Juan Islands, South Puget Sound and Washington outer coast, time-series monitoring showed rapid disease spread, high mortality rates in 2014, and continuing levels of wasting in the survivors in 2015. Peak prevalence of disease at 16 sites ranged to 100%, with an overall mean of 61%. Analysis of longitudinal data showed disease risk was correlated with both size and temperature and resulted in shifts in population size structure; adult populations fell to one quarter of pre-outbreak abundances. In laboratory experiments, time between development of disease signs and death was influenced by temperature in adults but not juveniles and adult mortality was 18% higher in the 19 °C treatment compared to the lower temperature treatments. While larger ochre stars developed disease signs sooner than juveniles, diseased juveniles died more quickly than diseased adults. Unusual 2-3 °C warm temperature anomalies were coincident with the summer 2014 mortalities. We suggest these warm waters could have increased the disease progression and mortality rates of SSWD in Washington State.

Parameter Estimations of Dynamic Energy Budget (DEB) Model over the Life History of a Key Antarctic Species: The Antarctic Sea Star Odontaster validus Koehler, 1906

Marine organisms in Antarctica are adapted to an extreme ecosystem including extremely stable temperatures and strong seasonality due to changes in day length. It is now largely accepted that Southern Ocean organisms are particularly vulnerable to global warming with some regions already being challenged by a rapid increase of temperature. Climate change affects both the physical and biotic components of marine ecosystems and will have an impact on the distribution and population dynamics of Antarctic marine organisms. To predict and assess the effect of climate change on marine ecosystems a more comprehensive knowledge of the life history and physiology of key species is urgently needed. In this study we estimate the Dynamic Energy Budget (DEB) model parameters for key benthic Antarctic species the sea star Odontaster validus using available information from literature and experiments. The DEB theory is unique in capturing the metabolic processes of an organism through its entire life cycle as a function of temperature and food availability. The DEB model allows for the inclusion of the different life history stages, and thus, becomes a tool that can be used to model lifetime feeding, growth, reproduction, and their responses to changes in biotic and abiotic conditions. The DEB model presented here includes the estimation of reproduction handling rules for the development of simultaneous oocyte cohorts within the gonad. Additionally it links the DEB model reserves to the pyloric caeca an organ whose function has long been ascribed to energy storage. Model parameters described a slowed down metabolism of long living animals that mature slowly. O. validus has a large reserve that—matching low maintenance costs- allow withstanding long periods of starvation. Gonad development is continuous and individual cohorts developed within the gonads grow in biomass following a power function of the age of the cohort. The DEB model developed here for O. validus allowed us to increase our knowledge on the ecophysiology of this species, providing new insights on the role of food availability and temperature on its life cycle and reproduction strategy.