Functional responses of an invasive mud crab across a salinity gradient

Environmental gradients may alter the ecological impacts of invasive alien species. In marine systems such as the Baltic Sea, current salinity is variable and seawater freshening is projected in future, potentially facilitating novel keystone predators. Here, we examine the influence of salinity variation in the western Baltic Sea (i.e. ambient 10, then 7 and 4 ppt) on the functional response (FR) of the Harris mud crab Rhithropanopeus harrisii towards benthic macroinvertebrate prey at different densities. Rhithropanopeus harrisii displayed a Type II FR across salinities towards larval chironomids, due to a consistently high resource consumption rate at low prey densities. Feeding rates were significantly reduced at 4 ppt (mean 6 chironomid prey killed day^-1) compared to 10 ppt and 7 ppt (9 killed day^-1). Search efficiencies tended to be greatest at 10 ppt, whereas handling times were shortest - and maximum feeding rate highest - at the intermediate 7 ppt. These results suggest a slight reduction in predatory impact by R. harrisii at lower salinities. Nevertheless, across most prey densities, FRs were not significantly different, indicating sustained interaction strength across a range of salinity regimes.

Density-dependent ecosystem service delivery under shifting temperatures by dung beetles

Increases in the frequency and magnitude of suboptimal temperatures as a result of climate change are subjecting insects to unprecedented stresses. This may negatively affect their fitness and the efficiency of their ecosystem service provision. Dung beetles are ecosystem service providers: through feeding on and burying dung, they facilitate nutrient recycling, secondary seed dispersal, parasite control, soil bioturbation and dung decomposition. As such, prediction of how dung beetles respond to multiple anthropogenic environmental changes is critical for the conservation of ecosystem services. Here, we quantified ecosystem services via dung utilisation and dung ball production in three telecoprid species: Allogymnopleurus indigaceous, Scarabaeus zambezianus and Khepher prodigiosus. We examined ecosystem service efficiency factorially under different beetle densities towards different dung masses and under three temperature treatments (21 °C, 28 °C and 35 °C). Khepher prodigiosus, exhibited greatest dung utilisation efficiency overall across dung masses, compared to both S. zambezianus and A. indigaceous. Dung removal was exhibited under all the tested temperatures by all tested species, and therefore the sub-optimal temperatures employed here did not fully inhibit ecosystem service delivery. However, emergent effects among temperatures, beetle species and beetle density further affected removal efficiency: S. zambezianus and A. indigaceous utilisation increased with both warming and beetle density, whereas K. prodigiosus performance was less temperature- and density-dependent. Beetles also tended to exhibit positive density-dependence as dung supply increased. The numbers of dung balls produced differed across species, and increased with temperature and densities, with S. zambezianus producing significantly most balls overall. Our study provides novel evidence for differential density-dependent ecosystem service delivery among species across stressful temperature regimes and emergent effects for dung mass utilisation. This information is essential for biodiversity-ecosystem-function and is critical for the conservation of functionally efficacious species, with implications for natural capital conservation policy in rapidly changing environments.

Attack by a common parasitoid stabilizes population dynamics of multi-host communities

We model the population dynamics of two host species attacked by a common parasitoid using a discrete-time formalism that captures their population densities from year to year. It is well known starting from the seminal work of Nicholson and Bailey that a constant parasitoid attack rate leads to an unstable host-parasitoid interaction. However, a Type III functional response, where the parasitoid attack rate accelerates with increasing host density stabilizes the population dynamics. We first consider a scenario where both host species are attacked by a parasitoid with the same Type III functional response. Our results show that sufficient fast acceleration of the parasitoid attack rate stabilizes the population dynamics of all three species. For two symmetric host species, the extent of acceleration needed to stabilize the three-species equilibrium is exactly the same as that needed for a single host-parasitoid interaction. However, asymmetry can lead to scenarios where the removal of a host species from a stable interaction destabilizes the interaction between the remaining host species and the parasitoid. Next, we consider a situation where one of the host species is attacked at a constant rate (i.e., Type I functional response), and the other species is attacked via a Type III functional response. We identify parameter regimes where a Type III functional response to just one of the host species stabilizes the three species interaction. In summary, our results show that a generalist parasitoid with a Type III functional response to one or many host species can play a key role in stabilizing population dynamics of host-parasitoid communities in apparent competition.

Breathing space: deoxygenation of aquatic environments can drive differential ecological impacts across biological invasion stages

The influence of climate change on the ecological impacts of invasive alien species (IAS) remains understudied, with deoxygenation of aquatic environments often-overlooked as a consequence of climate change. Here, we therefore assessed how oxygen saturation affects the ecological impact of a predatory invasive fish, the Ponto-Caspian round goby (Neogobius melanostomus), relative to a co-occurring endangered European native analogue, the bullhead (Cottus gobio) experiencing decline in the presence of the IAS. In individual trials and mesocosms, we assessed the effect of high, medium and low (90%, 60% and 30%) oxygen saturation on: (1) functional responses (FRs) of the IAS and native, i.e. per capita feeding rates; (2) the impact on prey populations exerted; and (3) how combined impacts of both fishes change over invasion stages (Pre-invasion, Arrival, Replacement, Proliferation). Both species showed Type II potentially destabilising FRs, but at low oxygen saturation, the invader had a significantly higher feeding rate than the native. Relative Impact Potential, combining fish per capita effects and population abundances, revealed that low oxygen saturation exacerbates the high relative impact of the invader. The Relative Total Impact Potential (RTIP), modelling both consumer species’ impacts on prey populations in a system, was consistently higher at low oxygen saturation and especially high during invader Proliferation. In the mesocosm experiment, low oxygen lowered RTIP where both species were present, but again the IAS retained high relative impact during Replacement and Proliferation stages at low oxygen. We also found evidence of multiple predator effects, principally antagonism. We highlight the threat posed to native communities by IAS alongside climate-related stressors, but note that solutions may be available to remedy hypoxia and potentially mitigate impacts across invasion stages.

Estimation for the bivariate quantile varying coefficient model with application to diffusion tensor imaging data analysis

Despite interest in the joint modeling of multiple functional responses such as diffusion properties in neuroimaging, robust statistical methods appropriate for this task are lacking. To address this need, we propose a varying coefficient quantile regression model able to handle bivariate functional responses. Our work supports innovative insights into biomedical data by modeling the joint distribution of functional variables over their domains and across clinical covariates. We propose an estimation procedure based on the alternating direction method of multipliers and propagation separation algorithms to estimate varying coefficients using a B-spline basis and an $L_2$ smoothness penalty that encourages interpretability. A simulation study and an application to a real-world neurodevelopmental data set demonstrates the performance of our model and the insights provided by modeling functional fractional anisotropy and mean diffusivity jointly and their association with gestational age and sex.

Responses of Mytilus galloprovincialis hemocytes to environmental strains of Vibrio parahaemolyticus, Vibrio alginolyticus, Vibrio vulnificus

Marine bivalves are exposed to different types of bacteria in the surrounding waters, in particular of the Vibrio genus. In the hemocytes of the mussel Mytilus spp. immune responses to different vibrios have been largely characterized. However, little information is available on the hemocyte responses to human pathogenic vibrios commonly detected in coastal waters and bivalve tissues that are involved in seafood-borne diseases. In this work, functional parameters of the hemocytes from the Mediterranean mussel M. galloprovincialis were evaluated in response to in vitro challenge with different vibrios isolated from environmental samples of the Adriatic sea (Italy): V. parahaemolyticus Conero, V. alginolyticus 1513 and V. vulnificus 509. V. parahaemolyticus ATCC 43996 was used for comparison. At the 50:1 bacteria hemocyte ratio, only V. parahaemolyticus strains induced significant lysosomal membrane destabilisation. Stimulation of extracellular lysozyme release, total ROS, O₂^- and NO production were observed, although to different extents and with distinct time courses for different vibrios, V. vulnificus 509 in particular. Further comparisons between V. parahaemolyticus Conero and V. vulnificus 509 showed that only the latter induced dysregulation of the phosphorylation state of p38 MAP Kinase and apoptotic processes. The results indicate that mussel hemocytes can mount an efficient immune response towards V. parahaemolyticus and V. alginolyticus strains, whereas V. vulnificus 509 may affect the hemocyte function. This is the first report on immune responses of mussels to local environmental isolates of human pathogenic vibrios. These data reinforce the hypothesis that Mytilus hemocytes show specific responses to different vibrio species and strains.

The effect of nanoparticles on plankton dynamics: a mathematical model

A simple modification of the Rosenzweig-MacArthur predator (zooplankton)-prey (phytoplankton) model with the interference of the predators by adding the effect of nanoparticles is proposed and analyzed. It is assumed that the effect of these particles has a potential to reduce the maximum physiological per-capita growth rate of the prey. The dynamics of nanoparticles is assumed to follow a simple Lotka-Volterra uptake term. Our study suggests that nanoparticle induce growth suppression of phytoplankton population can destabilize the system which leads to limit cycle oscillation. We also observe that if the contact rate of nanoparticles and phytoplankton increases, then the equilibrium densities of phytoplankton as well as zooplankton decrease. Furthermore, we observe that the depletion/removal of nanoparticles from the aquatic system plays a crucial role for the stable coexistence of both populations. Our investigation with various types of functional response suggests that Beddington functional response is the most appropriate representation of the interaction of phytoplankton-nanoparticles in comparison to other widely used functional responses.