Crispin M. Mutshinda

What drives community dynamics

The search for general mechanisms of community assembly is a major focus of community ecology. The common practice so far has been to examine alternative assembly theories using dichotomist approaches of the form neutrality versus niche, or compensatory dynamics versus environmental forcing. In reality, all these mechanisms will be operating, albeit with different strengths. While there have been different approaches to community structure and dynamics, including neutrality and niche differentiation, less work has gone into separating out the temporal variation in species abundances into relative contributions from different components. Here we use a refined statistical machinery to decompose temporal fluctuations in species abundances into contributions from environmental stochasticity and inter-/intraspecific interactions, to see which ones dominate. We apply the methodology to community data from a range of taxa. Our results show that communities are largely driven by environmental fluctuations, and that member populations are, to different extents, regulated through intraspecific interactions, the effects of interspecific interactions remaining broadly minor. By decomposing the temporal variation in this way, we have been able to show directly what has been previously inferred indirectly: compensatory dynamics are in fact largely outweighed by environmental forcing, and the latter tends to synchronize the population dynamics.

A multispecies perspective on ecological impacts of climatic forcing

1. In the prevailing context of concerns over climate change and its potential impacts on ecosystems, evaluating ecological consequences of climatic forcing has become a critical issue. 2. Historical data on the abundance of organisms have been extensively used to characterize the ecological effects of climatic forcing through specific weather and/or climatic variables, with most of the studies confined to single population models. 3. However, population responses to environmental fluctuations typically depend upon positive and negative feedbacks induced by interactions with other species. It is therefore important to integrate the insights gained from single population approaches into a multispecies perspective. 4. Here we combine the hierarchical Bayesian modelling approach with the state-space formulation to extend the scope of previously proposed models of population dynamics under climatic forcing to multi-species systems. 5. We use our model to analyse long-term macro-moth (Lepidoptera) community data from the Rothamsted Insect Survey network in the UK, using winter rainfall and winter temperature as environmental covariates. 6. The effects of the two weather variables were consistent across species, being negative for winter rainfall and positive for winter temperature. The two weather variables jointly explained 15-40% of the total environmental variation affecting the dynamics of individual species, and could explain up to 90% of covariances in species dynamics. 7. The contribution of interspecific interactions to community-level variation was found to be weak compared to the contributions of environmental forcing and intraspecific interactions.

Extended Bayesian LASSO for Multiple Quantitative Trait Loci Mapping and Unobserved Phenotype Prediction

The Bayesian LASSO (BL) has been pointed out to be an effective approach to sparse model representation and successfully applied to quantitative trait loci (QTL) mapping and genomic breeding value (GBV) estimation using genome-wide dense sets of markers. However, the BL relies on a single parameter known as the regularization parameter to simultaneously control the overall model sparsity and the shrinkage of individual covariate effects. This may be idealistic when dealing with a large number of predictors whose effect sizes may differ by orders of magnitude. Here we propose the extended Bayesian LASSO (EBL) for QTL mapping and unobserved phenotype prediction, which introduces an additional level to the hierarchical specification of the BL to explicitly separate out these two model features. Compared to the adaptiveness of the BL, the EBL is “doubly adaptive” and thus, more robust to tuning. In simulations, the EBL outperformed the BL in regard to the accuracy of both effect size estimates and phenotypic value predictions, with comparable computational time. Moreover, the EBL proved to be less sensitive to tuning than the related Bayesian adaptive LASSO (BAL), which introduces locus-specific regularization parameters as well, but involves no mechanism for distinguishing between model sparsity and parameter shrinkage. Consequently, the EBL seems to point to a new direction for QTL mapping, phenotype prediction, and GBV estimation.

Integrating the niche and neutral perspectives on community structure and dynamics.

Elucidating the mechanisms underlying the assembly and dynamics of ecological communities is a fundamental goal of ecology. Two conceptual approaches have emerged in this respect: the niche-assembly view and the neutral perspective. The debate as to which approach best explains the biodiversity patterns observed in nature is becoming outdated, as ecologists increasingly agree on the existence of a niche-neutral continuum of community dynamical behaviors. However, attempts to make the continuum idea operational and measurable remain sparse. Here, we propose a model-based approach to achieving this. The proposed methodology consists of separating out fluctuations in species abundances into niche-mediated and stochastic factors, linking the niche configuration to community dynamics through competition, and adding demographic stochasticity. This results in a comprehensive framework including neutrality and strict niche segregation as extreme cases. We develop an index of departure from neutral drift as a surrogate for community position on the niche-neutral continuum. We evaluate the performance of our modeling approach with simulated data, and subsequently use the model to analyze rodent web-trapping data from a real-world system. The model fitting is carried out with a Bayesian approach using Markov chain Monte Carlo simulation methods.

Environmental control of the dominant phytoplankton in the Cariaco basin: a hierarchical Bayesian approach

Abstract We develop a hierarchical Bayesian model linking the abundance of individual phytoplankton species with over a decade (1995–2011) of environmental data from the Cariaco Ocean Time Series Program in the Cariaco Basin, Venezuela, to characterize how phytoplankton respond to environmental forcing. Temperature, salinity, irradiance, and macronutrient concentrations account for 39% of the variation in log cell abundance across 67 species. Individual phytoplankton taxa varied widely in their response to these environmental variables. A principal component analysis of the environmental response profiles clearly distinguishes the responses of diatoms and dinoflagellates to environmental forcing. Phytoplankton abundance primarily varied with temperature, pH, and irradiance, with salinity and macronutrient concentrations acting as secondary drivers. In the aggregate, our results demonstrate that environmental changes, whether short-term or a result of climate change, should be expected to have dramatic consequences on the taxonomic composition of phytoplankton communities.

Ecological equivalence of species within phytoplankton functional groups

1.There are tens of thousands of species of phytoplankton found throughout the tree of life. Despite this diversity, phytoplankton are often aggregated into a few functional groups according to metabolic traits or biogeochemical role. We investigate the extent to which phytoplankton species dynamics are neutral within functional groups. 2.Seasonal dynamics in many regions of the ocean are known to affect phytoplankton at the functional group level leading to largely predictable patterns of seasonal succession. It is much more difficult to make general statements about the dynamics of individual species. 3.We use a 7 year time-series at station L4 in the Western English Channel with 57 diatom and 17 dinoflagellate species enumerated weekly to test if the abundance of diatom and dinoflagellate species vary randomly within their functional group envelope or if each species is driven uniquely by external factors. 4.We show that the total biomass of the diatom and dinoflagellate functional groups is well predicted by irradiance and temperature and quantify trait values governing the growth rate of both functional groups. The biomass dynamics of the functional groups are not neutral and each has their own distinct responses to environmental forcing. Compared to dinoflagellates, diatoms have faster growth rates, and grow faster under lower irradiance, cooler temperatures, and higher nutrient conditions. 5.The biomass of most species vary randomly within their functional group biomass envelope, most of the time. As a consequence, modelers will find it difficult to predict the biomass of most individual species. Our analysis supports the approach of using a single set of traits for a functional group and suggests that it should be possible to determine these traits from natural communities.

Which UGC features drive web purchase intent? A spike-and-slab Bayesian Variable Selection Approach

– The purpose of this paper is to identify user-generated content (UGC) features that determine web purchase decision making. , – The authors embed a spike-and-slab Bayesian variable selection mechanism into a logistic regression model to identify the UGC features that are critical to web purchase intent. This enables us to make a highly reliable analysis of survey data. , – The results indicate that the web purchase decision is driven by the relevance, up-to-dateness and credibility of the UGC information content. , – The results show that the characteristics of UGC are seen as positive and the medium enables consumers to sort information and concentrate on aspects of the message that are similar to traditional word-of-mouth (WOM). One important implication is the relative importance of credibility which has been previously hypothesized to be lower in the electronic word-of-mouth (e-WOM) context. The results show that consumers consider credibility important as the improved technology provides more possibilities to find out about that factor. A limitation is that the data are not fully representative of the general population but our Bayesian method gives us high analytical quality. , – The study shows that UGC impacts consumer online purchase intentions. Marketers should understand the wide range of media that provide UGC and they should concentrate on the relevance, up-to-dateness and credibility of product information that they provide. , – The analytical quality of the spike- and- slab Bayesian method suggests a new way of understanding the impact of aspects of UGC on consumers.

Swift block-updating EM and pseudo-EM procedures for Bayesian shrinkage analysis of quantitative trait loci.

IntroductionVirtually all existing expectation-maximization (EM) algorithms for quantitative trait locus (QTL) mapping overlook the covariance structure of genetic effects, even though this information can help enhance the robustness of model-based inferences.ResultsHere, we propose fast EM and pseudo-EM-based procedures for Bayesian shrinkage analysis of QTLs, designed to accommodate the posterior covariance structure of genetic effects through a block-updating scheme. That is, updating all genetic effects simultaneously through many cycles of iterations.ConclusionSimulation results based on computer-generated and real-world marker data demonstrated the ability of our method to swiftly produce sensible results regarding the phenotype-to-genotype association. Our new method provides a robust and remarkably fast alternative to full Bayesian estimation in high-dimensional models where the computational burden associated with Markov chain Monte Carlo simulation is often unwieldy. The R code used to fit the model to the data is provided in the online supplementary material.

A Decision Rule for Quantitative Trait Locus Detection Under the Extended Bayesian LASSO Model

Bayesian shrinkage analysis is arguably the state-of-the-art technique for large-scale multiple quantitative trait locus (QTL) mapping. However, when the shrinkage model does not involve indicator variables for marker inclusion, QTL detection remains heavily dependent on significance thresholds derived from phenotype permutation under the null hypothesis of no phenotype-to-genotype association. This approach is computationally intensive and more importantly, the hypothetical data generation at the heart of the permutation-based method violates the Bayesian philosophy. Here we propose a fully Bayesian decision rule for QTL detection under the recently introduced extended Bayesian LASSO for QTL mapping. Our new decision rule is free of any hypothetical data generation and relies on the well-established Bayes factors for evaluating the evidence for QTL presence at any locus. Simulation results demonstrate the remarkable performance of our decision rule. An application to real-world data is considered as well.

Bayesian shrinkage analysis of QTLs under shape-adaptive shrinkage priors, and accurate re-estimation of genetic effects

The successful implementation of Bayesian shrinkage analysis of high-dimensional regression models, as often encountered in quantitative trait locus (QTL) mapping, is contingent upon the choice of suitable sparsity-inducing priors. In practice, the shape (that is, the rate of tail decay) of such priors is typically preset, with no regard for the range of plausible alternatives and the fact that the most appropriate shape may depend on the data at hand. This study is presumably the first attempt to tackle this oversight through the shape-adaptive shrinkage prior (SASP) approach, with a focus on the mapping of QTLs in experimental crosses. Simulation results showed that the separation between genuine QTL effects and spurious ones can be made clearer using the SASP-based approach as compared with existing competitors. This feature makes our new method a promising approach to QTL mapping, where good separation is the ultimate goal. We also discuss a re-estimation procedure intended to improve the accuracy of the estimated genetic effects of detected QTLs with regard to shrinkage-induced bias, which may be particularly important in large-scale models with collinear predictors. The re-estimation procedure is relevant to any shrinkage method, and is potentially valuable for many scientific disciplines such as bioinformatics and quantitative genetics, where oversaturated models are booming.