Hierarchical Bayesian models in ecology: Reconstructing species interaction networks from non-homogeneous species abundance data

<p>The relationships among organisms and their surroundings can be of immense complexity. To describe and understand an ecosystem as a tangled bank, multiple ways of interaction and their effects have to be considered, such as predation, competition, mutualism and facilitation. Understanding the resulting interaction networks is a challenge in changing environments, e.g. to predict knock-on effects of invasive species and to understand how climate change impacts biodiversity. The elucidation of complex ecological systems with their interactions will benefit enormously from the development of new machine learning tools that aim to infer the structure of interaction networks from field data. In the present study, we propose a novel Bayesian regression and multiple changepoint model (BRAM) for reconstructing species interaction networks from observed species distributions. The model has been devised to allow robust inference in the presence of spatial autocorrelation and distributional heterogeneity. We have evaluated the model on simulated data that combines a trophic niche model with a stochastic population model on a 2-dimensional lattice, and we have compared the performance of our model with L1-penalized sparse regression (LASSO) and non-linear Bayesian networks with the BDe scoring scheme. In addition, we have applied our method to plant ground coverage data from the western shore of the Outer Hebrides with the objective to infer the ecological interactions. (C) 2012 Elsevier B.V. All rights reserved.</p>

Echinostoma caproni: Differential tegumental responses to growth in compatible and less compatible hosts

The topography of the tegument of Echinostoma caproni adults collected from high (mice) and low (rats) compatible hosts was compared by SEM. In the oral (OS) and the ventral sucker (VS) areas, a worm age-host species interaction was found with regard to the density of spines. There was a decrease in the density of spines in the adults collected from mice, whereas an increase occurred in the OS area in worms from rats over time. The tegumentary spines in adults from mice became larger and blunter. Some spines from the VS area in adults from mice at 4 wpi were multipointed. The spines of adults from rats were sharper, not covered by the tegument and no multipointed spines were observed. We detected a greater level of actin gene expression in the adults collected from rats. These facts suggest that the low compatible host induces an increased turnover of tegumentary spines. (C) 2010 Elsevier Inc. All rights reserved.

Predicting community responses to perturbations in the face of imperfect knowledge and network complexity

How best to predict the effects of perturbations to ecological communities has been a long-standing goal for both applied and basic ecology. This quest has recently been revived by new empirical data, new analysis methods, and increased computing speed, with the promise that ecologically important insights may be obtainable from a limited knowledge of community interactions. We use empirically based and simulated networks of varying size and connectance to assess two limitations to predicting perturbation responses in multispecies communities: (1) the inaccuracy by which species interaction strengths are empirically quantified and (2) the indeterminacy of species responses due to indirect effects associated with network size and structure. We find that even modest levels of species richness and connectance (similar to 25 pairwise interactions) impose high requirements for interaction strength estimates because system indeterminacy rapidly overwhelms predictive insights. Nevertheless, even poorly estimated interaction strengths provide greater average predictive certainty than an approach that uses only the sign of each interaction. Our simulations provide guidance in dealing with the trade-offs involved in maximizing the utility of network approaches for predicting dynamics in multispecies communities.

Weak interactions, omnivory and emergent food-web properties

Empirical studies have shown that, in real ecosystems, species-interaction strengths are generally skewed in their distribution towards weak interactions. Some theoretical work also suggests that weak interactions, especially in omnivorous links, are important for the local stability of a community at equilibrium. However, the majority of theoretical studies use uniform distributions of interaction strengths to generate artificial communities for study. We investigate the effects of the underlying interaction-strength distribution upon the return time, permanence and feasibility of simple Lotka-Volterra equilibrium communities. We show that a skew towards weak interactions promotes local and global stability only when omnivory is present. It is found that skewed interaction strengths are an emergent property of stable omnivorous communities, and that this skew towards weak interactions creates a dynamic constraint maintaining omnivory. Omnivory is more likely to occur when omnivorous interactions are skewed towards weak interactions. However, a skew towards weak interactions increases the return time to equilibrium, delays the recovery of ecosystems and hence decreases the stability of a community. When no skew is imposed, the set of stable omnivorous communities shows an emergent distribution of skewed interaction strengths. Our results apply to both local and global concepts of stability and are robust to the definition of a feasible community. These results are discussed in the light of empirical data and other theoretical studies, in conjunction with their broader implications for community assembly.

Interaction strength, food web topology and the relative importance of species in food webs

P&gt;1. We established complex marine communities, consisting of over 100 species, in large subtidal experimental mesocosms. We measured the strength of direct interactions and the net strength of direct and indirect interactions between the species in those communities, using a combination of theoretical and empirical approaches.

The human micro-vascular enclothelial cells in vitro interaction with atomic-nitrogen-doped diamond-like carbon thin films

This paper reports the initial response of atomic nitrogen doped diamond like carbon (DLC) to endothelial cells in vitro. The introduction of nitrogen atoms/molecules to the diamond like carbon structures leads to an atomic structural change favorable to the attachment of human micro-vascular enclothelial cells. Whilst the semi-conductivity induced by nitrogen in DLC is thought to play a part, the increase in the inion-bonded N atoms and N-2 molecules in the atomic doped species (with the exclusion of the charged species) seems to contribute to the improved attachment of human microvascular endothelial cells. The increased endothelial attachment is associated with a lower work function and slightly higher water contact angle in the atomic doped films, where the heavy charged particles are excluded. The films used in the study were synthesized by the RF PECVD technique followed by post deposition doping with nitrogen, and afterwards the films were characterized by XPS, Raman spectroscopy, SIMS and Kelvin probe. The water contact angles were measured, and the counts of the adherent endothelial cells on the samples were carried out. This study is relevant and contributory to improving biocompatibility of surgical implants and prostheses.

Invader-invader interactions in relation to environmental heterogeneity leads to zonation of two invasive amphipods, Dikerogammarus villosus (Sowinsky) and Gammarus tigrinus Sexton: amphipod pilot species project (AMPIS) report 6

As biological invasions continue, interactions occur not only between invaders and natives, but increasingly new invaders come into contact with previous invaders. Whilst this can lead to species replacements, co-existence may occur, but we lack knowledge of processes driving such patterns. Since environmental heterogeneity can determine species richness and co-existence, the present study examines habitat use and its mediation of the predatory interaction between invasive aquatic amphipods, the Ponto-Caspian Dikerogammarus villosus and the N. American Gammarus tigrinus. In the Dutch Lake IJsselmeer, we found broad segregation of D. villosus and G. tigrinus by habitat type, the former predominating in the boulder zone and the latter in the soft sediment. However, the two species co-exist in the boulder zone, both on the short and longer terms. We used an experimental simulation of habitat heterogeneity and show that both species utilize crevices, different sized holes in a plastic grid, non-randomly. These amphipods appear to optimise the use of holes with respect to their 'C-shape' body size. When placed together, D. villosus adults preyed on G. tigrinus adults and juveniles, while G. tigrinus adults preyed on D. villosus juveniles. Juveniles were also predators and both species were cannibalistic. However, the impact on G. tigrinus of the superior intraguild predator, D. villosus, was significantly reduced where experimental grids were present as compared to absent. This mitigation of intraguild predation between the two species in complex habitats may explain the co-existence of these two invasive species.

Role of behaviour in biological invasions and species distributions; lessons from interactions between the invasive Gammarus pulex and the native G-duebeni (Crustacea : Amphipoda)

In attempting to understand the distributions of both introduced species and the native species on which they impact, there is a growing trend to integrate studies of behaviour with more traditional life history/ecological approaches. The question of what mechanisms drive the displacement of the freshwater amphipod Gammarus duebeni by the often introduced G pulex is presented as a case study Patterns of displacement are well documented throughout Europe, but the speed and direction of displacement between these species can be varied. From early studies proposing interspecific competition as causal in these patterns, I review research progress to date. I show there has been no evidence for interspecific competition operating, other than the field patterns themselves, a somewhat tautological argument. Rather, the increased recognition of behavioural attributes with respect to the cannibalistic and predatory nature of these species gave rise to a series of studies unravelling the processes driving field patterns. Both species engage in 'intraguild predation' (IGP), with moulting females particularly vulnerable to predation by congeneric males. G pulex is more able both to engage in and avoid this interaction with G duebeni. However, several factors mediate the strength and asymmetry of this IGP, some biotic (e.g. parasitism) and others abiotic (e.g. water chemistry). Further, a number of alternative hypotheses that may account for the displacement (hybridization; parasite transmission) have been tested and rejected. While interspecific competition has been modelled mathematically and found to be a weak interaction relative to IGP, mechanisms of competition between these Gammarus species remain largely untested empirically. Since IGP may be finely balanced in some circumstances, I conclude that the challenge to detect interspecific competition remains and we require assessment of its role, if any, in the interaction between these species. Appreciation of behavioural attributes and their mediation should allow us to more fully understand, and perhaps predict, species introductions and resultant distributions.

Loss of functionally unique species may gradually undermine ecosystems

Functionally unique species contribute to the functional diversity of natural systems, often enhancing ecosystem functioning. An abundance of weakly interacting species increases stability in natural systems, suggesting that loss of weakly linked species may reduce stability. Any link between the functional uniqueness of a species and the strength of its interactions in a food web could therefore have simultaneous effects on ecosystem functioning and stability. Here, we analyse patterns in 213 real food webs and show that highly unique species consistently tend to have the weakest mean interaction strength per unit biomass in the system. This relationship is not a simple consequence of the interdependence of both measures on body size and appears to be driven by the empirical pattern of size structuring in aquatic systems and the trophic position of each species in the web. Food web resolution also has an important effect, with aggregation of species into higher taxonomic groups producing a much weaker relationship. Food webs with fewer unique and less weakly interacting species also show significantly greater variability in their levels of primary production. Thus, the loss of highly unique, weakly interacting species may eventually lead to dramatic state changes and unpredictable levels of ecosystem functioning.

Manipulating Interaction Strengths and the Consequences for Trivariate Patterns in a Marine Food Web

We are experiencing a global extinction crisis as a result of climate change and human-induced alteration of natural habitats, with large predators at high trophic levels in food webs being particularly vulnerable. Unfortunately, there is a scarcity of food web data that can be used to assess how species extinctions alter the structure and stability of temporally and spatially replicated networks. We established a series of large experimental mesocosms in a shallow subtidal benthic marine system and constructed food webs for each replicate. After 6 months of community assembly, we removed large predators from the core communities of 20 experimental food webs, based on the strength of their trophic interactions, and monitored the changes in the networks' structure and stability over an 8-month period. Our analyses revealed the importance of allometric relationships and size-structuring in natural communities as a means of preserving food web structure and sustainability, despite significant changes in the diversity, stability and productivity of the system.

Measurement of interaction strength innature

Understanding and predicting the dynamics of multispecies systems generally require estimates of interaction strength among species. Measuring interaction strength is difficult because of the large number of interactions in any natural system, long-term feedback, multiple pathways of effects between species pairs, and possible nonlinearities in interaction-strength functions. Presently, the few studies that extensively estimate interaction strength suggest that distributions of interaction strength tend to be skewed toward few strong and many weak interactions. Modeling studies indicate that such skewed patterns tend to promote system stability and arise during assembly of persistent communities. Methods for estimating interaction strength efficiently from traits of organisms, such as allometric relationships, show some promise. Methods for estimating community response to environmental perturbations without an estimate of interaction strength may also be of use. Spatial and temporal scale may affect patterns of interaction strength, but these effects require further investigation and new multispecies modeling frameworks. Future progress will be aided by development of long-term multispecies time series of natural communities, by experimental tests of different methods for estimating interaction strength, and by increased understanding of nonlinear functional forms.

Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs

When it is difficult to develop selective ligands within a family of related G-protein-coupled receptors (GPCRs), chemically engineered receptors activated solely by synthetic ligands (RASSLs) are useful alternatives for probing receptor function. In the present work, we explored whether a RASSL of the free fatty acid receptor 2 (FFA2) could be developed on the basis of pharmacological variation between species orthologs. For this, bovine FFA2 was characterized, revealing distinct ligand selectivity compared with human FFA2. Homology modeling and mutational analysis demonstrated a single mutation in human FFA2 of C4.57G resulted in a human FFA2 receptor with ligand selectivity similar to the bovine receptor. This was exploited to generate human FFA2-RASSL by the addition of a second mutation at a known orthosteric ligand interaction site, H6.55Q. The resulting FFA2-RASSL displayed a &gt;100-fold loss of activity to endogenous ligands, while responding to the distinct ligand sorbic acid with pEC(50) values for inhibition of cAMP, 5.83 ± 0.11; Ca(2+) mobilization, 4.63 ± 0.05; ERK phosphorylation, 5.61 ± 0.06; and dynamic mass redistribution, 5.35 ± 0.06. This FFA2-RASSL will be useful in future studies on this receptor and demonstrates that exploitation of pharmacological variation between species orthologs is a powerful method to generate novel chemically engineered GPCRs.-Hudson, B. D., Christiansen, E., Tikhonova, I. G., Grundmann, M., Kostenis, E., Adams, D. R., Ulven, T., Milligan, G. Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs.

Possible interactions between environmental factors in determining species optima

<p>Questions: 1. Indicator values, such as those of Ellenberg, for different environmental factors are seen as independent. We tested for the presence of interactions between environmental factors ( soil moisture and reaction) to see if this assumption is simplistic. 2. How close are Ellenberg indicator values (IVs) related to the observed optima of species response curves in an area peripheral to those where they have been previously employed and 3. Can the inclusion of bryophytes add to the utility of IVs?</p><p>Location: South Uist, Outer Hebrides, Scotland, UK.</p><p>Methods: Two grids (ca. 2000 m x 2000 m) were sampled at 50-m intervals across the transition from machair to upland communities covering an orthogonal gradient of both soil pH ( reaction) and soil moisture content. Percentage cover data for vascular plants, bryophytes and lichens were recorded, along with pH and moisture content of the underlying sand/soil/peat. Reaction optima, derived from species response curves calculated using HOF models, were compared between wet and dry sites, and moisture optima between acidic and basic samples. Optima for the whole data set were compared to Ellenberg IVs to assess their performance in this area, with and without the inclusion of bryophytes.</p><p>Results: A number of species showed substantially different pH optima at high and low soil moisture contents (18% of those tested) and different soil moisture optima at high and low pH (49%). For a number of species the IVs were poor predictors of their actual distribution across the sampled area. Bryophytes were poor at explaining local variation in the environmental factors and also their inclusion with vascular plants negatively affected the strength of relationships.</p><p>Conclusions: A substantial number of species showed an interaction between soil moisture and reaction in determining their optima on the two respective gradients. It should be borne in mind that IVs such as Ellenberg's may not be independent of one another.</p>

Large-scale spatial patterns in the distribution of Collembola (Hexapoda) species in Antarctic terrestrial ecosystems

<p>We tested whether the distribution of three common springtail species (Gressittacantha terranova, Gomphiocephalus hodgsoni and Friesea grisea) in Victoria Land (Antarctica) could be modelled as a function of latitude, longitude, altitude and distance from the sea.</p><p>Victoria Land, Ross Dependency, Antarctica.</p><p>Generalized linear models were constructed using species presence/absence data relative to geographical features (latitude, longitude, altitude, distance from sea) across the species' entire ranges. Model results were then integrated with the known phylogeography of each species and hypotheses were generated on the role of climate as a major driver of Antarctic springtail distribution.</p><p>Based on model selection using Akaike's information criterion, the species' distributions were: hump-shaped relative to longitude and monotonic with altitude for Gressittacantha terranova; hump-shaped relative to latitude and monotonic with altitude for Gomphiocephalus hodgsoni; and hump-shaped relative to longitude and monotonic with latitude, altitude and distance from the sea for Friesea grisea.</p><p>No single distributional pattern was shared by the three species. While distributions were partially a response to climatic spatial clines, the patterns observed strongly suggest that past geological events have influenced the observed distributions. Accordingly, present-day spatial patterns are likely to have arisen from the interaction of historical and environmental drivers. Future studies will need to integrate a range of spatial and temporal scales to further quantify their respective roles.</p>