Dynamics of neutral biodiversity

Hubbell's neutral model has become a major paradigm in ecology. Whereas the steady-state structure is well understood, results about the dynamical aspects of the model are scarce. Here we derive dynamical equations for the Simpson diversity index. Both mean and variance of the diversity are proven to satisfy stable linear system dynamics. We show that in the stationary limit we indeed recover previous results, and we supplement this with numerical simulations to validate the dynamical part of our analytical computations. These findings are especially relevant for experiments in microbial ecology, where the Simpson diversity index can be accurately measured as a function of time.     

The zero-sum assumption in neutral biodiversity theory

The neutral theory of biodiversity as put forward by Hubbell in his 2001 monograph has received much criticism for its unrealistic simplifying assumptions. These are the assumptions of functional equivalence among different species (neutrality), the assumption of point mutation speciation, and the assumption that resources are continuously saturated, such that constant resource availability implies constant community size (zero-sum assumption). Here we focus on the zero-sum assumption. We present a general theory for calculating the probability of observing a particular species-abundance distribution (sampling formula) and show that zero-sum and non-zero-sum formulations of neutral theory have exactly the same sampling formula when the community is in equilibrium. Moreover, for the non-zero-sum community the sampling formula has this same form, even out of equilibrium. Therefore, the term "zero-sum multinomial (ZSM)" to describe species abundance patterns, as coined by Hubbell [2001. The Unified Neutral Theory of Biodiversity and Biogeography, Princeton University Press, Princeton, NJ], is not really appropriate, as it also applies to non-zero-sum communities. Instead we propose the term "dispersal-limited multinomial (DLM)", thus making explicit one of the most important contributions of neutral community theory, the emphasis on dispersal limitation as a dominant factor in determining species abundances.     

Regression models for spatial prediction: their role for biodiversity and conservation

This paper is an introduction to a Special Issue on regression models for spatial predictions published in Biodiversity and Conservation following an international workshop held in Switzerland in 2001 (http://leba.unige.ch/workshop). This introduction describes how the exponential growth in computing power has improved our ability to reach spatially explicit assessment of biodiversity and to develop cost-effective conservation management. New questions arising from these modern approaches are listed, while papers presenting examples of applications are briefly introduced.     

The influence of spatial inhomogeneities on neutral models of geographical variation : I. Formulation

The influence of spatial variation in the carrying capacity and migration rate of a geographical barrier on the one-dimensional stepping-stone model is studied. The monoecious, diploid population is subdivided into an infinite linear array of panmictic colonies that exchange gametes. In each deme, the rate of self-fertilization is equal to the reciprocal of the number of individuals in that deme. Generations are discrete and nonoverlapping; the analysis is restricted to a single locus in the absence of selection; every allele mutates to new alleles at the same rate. In the diffusion approximation, a partial differential equation that incorporates spatial (and temporal) variation in the carrying capacity and migration rate is derived for the probability of identity. Transition conditions that simultaneously take into account discontinuities in the carrying capacity and migration rate are established: the probability of identity is continuous, but its partial derivatives are not, their ratio being a simple function of the carrying capacities and migrational variance on the two sides of the inhomogeneity. The partial derivatives of the probability of identity are continuous across a geographical barrier, whereas the probability of identity itself has a discontinuity proportional to the partial derivative at the barrier, the constant of proportionality being a measure of the difficulty of crossing the barrier.     

The unified neutral theory of biodiversity and biogeography at age ten

A decade has now passed since Hubbell published The Unified Neutral Theory of Biodiversity and Biogeography. Neutral theory highlights the importance of dispersal limitation, speciation and ecological drift in the natural world and provides quantitative null models for assessing the role of adaptation and natural selection. Significant advances have been made in providing methods for understanding neutral predictions and comparing them with empirical data. In this review, we describe the current state-of-the-art techniques and ideas in neutral theory and how these are of relevance to ecology. The future of neutral theory is promising, but its concepts must be applied more broadly beyond the current focus on species-abundance distributions.     

A new sampling formula for neutral biodiversity

The neutral model of biodiversity, proposed by Hubbell (The Unified Neutral Theory of Biodiversity and Biogeography, Princeton University Press, Princeton, NJ, 2001) to explain the diversity of functionally equivalent species, has been subject of hot debate in community ecology. Whereas Hubbell studied the model mostly by simulations, recently analytical treatments have yielded expressions of the expected number of species of a particular abundance in a local community with dispersal limitation. Moreover, a formula has been offered for the joint likelihood of observing a given species‐abundance dataset in a local community with dispersal limitation, but this formula is too complicated to allow practical applications. Here, I present a much simplified expression that can be regarded as an enhanced version of the famous Ewens sampling formula. It can be used in maximum likelihood methods for quick estimation of the model parameters, using all information in the data, and for model comparison. I also show how to rapidly generate examples of species‐abundance distributions for a given set of model parameters and how to calculate Simpson's diversity index.     

Random sampling does not exclude spatial dependence: The importance of neutral models for ecological hypothesis testing

Lájer (2007) notes that, to investigate phytosociological and ecological relationships, many authors apply traditional inferential tests to sets of relevés obtained by non-random methods. Unfortunately, this procedure does not provide reliable support for hypothesis testing because non-random sampling violates the assumptions of independence required by many parametric inferential tests. Instead, a random sampling scheme is recommended. Nonetheless, random sampling will not eliminate spatial autocorrelation. For instance, a classical law of geography holds that everything in a piece of (biotic) space is interrelated, but near objects are more related than distant ones. Because most ecological processes that shape community structure and species coexistence are spatially explicit, spatial autocorrelation is a vital part of almost all ecological data. This means that, independently from the underlying sampling design, ecological data are generally spatially autocorrelated, violating the assumption of independence that is generally required by traditional inferential tests. To overcome this drawback, randomization tests may be used. Such tests evaluate statistical significance based on empirical distributions generated from the sample and do not necessarily require data independence. However, as concerns hypothesis testing, randomization tests are not the universal remedy for ecologists, because the choice of inadequate null models can have significant effects on the ecological hypotheses tested. In this paper, I emphasize the need of developing null models for which the statistical assumptions match the underlying biological mechanisms.     

The neutral theory is dead. The current significance and standing of neutral and nearly neutral theories

Comparative studies of DNA sequences provide opportunities for testing the neutral and the selection theories of molecular evolution. In particular, the separate estimation of the numbers of synonymous and nonsynonymous substitutions is a powerful tool for detecting selection of the latter. The difference in the patterns of these two types of substitutions of mammalian genes turned out to be in accord with the slightly deleterious or nearly neutral mutation theory for nonsynonymous changes. Interaction systems at the amino acid level were suggested to be responsible for such nearly neutral, or very weak, selection. Synonymous substitutions are not strictly neutral, but because of their minute effect, random drift predominates such that the rate of substitution is only slightly less than the completely neutral prediction. It was concluded that the strictly neutral theory has not held up as well as the nearly neutral theory, yet remains invaluable as a null hypothesis for detecting selection. On the other hand, the main difference between the nearly neutral and the traditional selection theories is that the former predicts rapid evolution in small populations, whereas the latter predicts rapid evolution in large populations.     

Demographic analysis of Hubbell's neutral theory of biodiversity

Hubbell's neutral model is increasingly applied in both theoretical and empirical studies but so far little attention has been paid to the ecological mechanisms that determine species diversity in neutral communities. In this contribution we use a stochastic individual-based Markovian model to provide an explicit derivation of Hubbell's local community model from the fundamental processes of reproduction, mortality, and immigration, and show that such derivation provides important insights on the mechanisms regulating species diversity that cannot be obtained from the original model and its previous extensions. One important insight is that the basic parameters of Hubbell's model, community size (J) and the probability that a dying individual will be replaced by an immigrant (m), cannot be considered independent and that their interdependency leads to a counterintuitive trade-off between community size and species diversity. We further demonstrate that Hubbell's treatment of community size as a free parameter hides fundamental mechanisms that influence species diversity through their effect on the size of the community. For example, while in Hubbell's model immigration can only increase species diversity by promoting colonization rates, the demographic derivation shows that immigration can also promote species diversity by reducing extinction rates. Our demographic derivation also unifies previous contrasting predictions about the effect of reproduction on species diversity by showing that both positive and negative effects are possible, and that the balance between the two effects depends on the size of the community. The demographic derivation also reconciles an apparent contradiction between Hubbell's theory and patch occupancy theory, and integrates three previously proposed mechanisms of species diversity, the More Individuals Hypothesis, the rescue effect, and the dilution effect, within a single, unified framework.     

Relaxing the zero-sum assumption in neutral biodiversity theory

The zero-sum assumption is one of the ingredients of the standard neutral model of biodiversity by Hubbell. It states that the community is saturated all the time, which in this model means that the total number of individuals in the community is constant over time, and therefore introduces a coupling between species abundances. It was shown recently that a neutral model with independent species, and thus without any coupling between species abundances, has the same sampling formula (given a fixed number of individuals in the sample) as the standard model [Etienne, R.S., Alonso, D., McKane, A.J., 2007. The zero-sum assumption in neutral biodiversity theory. J. Theor. Biol. 248, 522-536]. The equilibria of both models are therefore equivalent from a practical point of view. Here we show that this equivalence can be extended to a class of neutral models with density-dependence on the community-level. This result can be interpreted as robustness of the model, i.e. insensitivity of the model to the precise interaction of the species in a neutral community. It can also be interpreted as a lack of resolution, as different mechanisms of interactions between neutral species cannot be distinguished using only a single snapshot of species abundance data.     

The robustness of neutral models of geographical variation

The approximation of diploid migration by gametic dispersion is studied. The monoecious, diploid population is subdivided into panmictic colonies that exchange migrants. Generations are discrete and nonoverlapping; the analysis is restricted to a single locus in the absence of selection; every allele mutates to a new allele at the same rate u. Diploid-migration models without self-fertilization and with selfing at the “random” rate (equal to the reciprocal of the deme size in each deme) are investigated; in the gametic-dispersion models, selfing occurs at the random rate. It is shown for the unbounded stepping-stone model in one and two dimensions, the circular stepping-stone model, and the island model that the probabilitities of identity in state at equilibrium for diploid migration are close to those for gametic dispersion if the mutation rate is small or the deme size is large. Explicit error bounds are presented in all the above cases. It is also proved that if the number of demes is finite and the migration matrix is arbitrary but time independent and ergodic, then in the strong-migration approximation the equilibrium and the ultimate rate and pattern of convergence of both diploid-dispersion models are close to the corresponding gametic-dispersion formulae. For the strong-migration approximation at equilibrium, migration must dominate both mutation and random drift; for the convergence results, it suffices that migration dominate random drift. All the results apply to a dioecious population if the migration pattern and mutation rate are sex independent.     

Modes of speciation and the neutral theory of biodiversity

Hubbell's neutral theory of biodiversity has generated much debate over the need for niches to explain biodiversity patterns. Discussion of the theory has focused on its neutrality assumption, i.e. the functional equivalence of species in competition and dispersal. Almost no attention has been paid to another critical aspect of the theory, the assumptions on the nature of the speciation process. In the standard version of the neutral theory each individual has a fixed probability to speciate. Hence, the speciation rate of a species is directly proportional to its abundance in the metacommunity. We argue that this assumption is not realistic for most speciation modes because speciation is an emergent property of complex processes at larger spatial and temporal scales and, consequently, speciation rate can either increase or decrease with abundance. Accordingly, the assumption that speciation rate is independent of abundance (each species has a fixed probability to speciate) is a more natural starting point in a neutral theory of biodiversity. Here we present a neutral model based on this assumption and we confront this new model to 20 large data sets of tree communities, expecting the new model to fit the data better than Hubbell's original model. We find, however, that the data sets are much better fitted by Hubbell's original model. This implies that species abundance data can discriminate between different modes of speciation, or, stated otherwise, that the mode of speciation has a large impact on the species abundance distribution. Our model analysis points out new ways to study how biodiversity patterns are shaped by the interplay between evolutionary processes (speciation, extinction) and ecological processes (competition, dispersal).     

A neutral metapopulation model of biodiversity in river networks

In this paper, we develop a stochastic, discrete, structured metapopulation model to explore the dynamics and patterns of biodiversity of riparian vegetation. In the model, individual plants spread along a branched network via directional dispersal and undergo neutral ecological drift. Simulation results suggest that in comparison to 2-D landscapes with non-directional dispersal, river networks with directional dispersal have lower local (alpha) and overall (gamma) diversities, but higher between-community (beta) diversity, implying that riparian species are distributed in a more localized pattern and more vulnerable to local extinction. The relative abundance patterns also change, such that higher percentages of species are in low-abundance, or rare, classes, accompanied by concave rank-abundance curves. In contrast to existing theories, the results suggest that in river networks, increased directional dispersal reduces alpha diversity. These altered patterns and trends result from the combined effects of directionality of dispersal and river network structure, whose relative importance is in need of continuing study. In addition, riparian communities obeying neutral dynamics seem to exhibit abrupt changes where large tributaries confluence; this pattern may provide a signature to identify types of interspecific dynamics in river networks.     

Understanding how biodiversity unfolds through time under neutral theory

Theoretical predictions for biodiversity patterns are typically derived under the assumption that ecological systems have reached a dynamic equilibrium. Yet, there is increasing evidence that various aspects of ecological systems, including (but not limited to) species richness, are not at equilibrium. Here, we use simulations to analyse how biodiversity patterns unfold through time. In particular, we focus on the relative time required for various biodiversity patterns (macroecological or phylogenetic) to reach equilibrium. We simulate spatially explicit metacommunities according to the Neutral Theory of Biodiversity (NTB) under three modes of speciation, which differ in how evenly a parent species is split between its two daughter species. We find that species richness stabilizes first, followed by species area relationships (SAR) and finally species abundance distributions (SAD). The difference in timing of equilibrium between these different macroecological patterns is the largest when the split of individuals between sibling species at speciation is the most uneven. Phylogenetic patterns of biodiversity take even longer to stabilize (tens to hundreds of times longer than species richness) so that equilibrium predictions from neutral theory for these patterns are unlikely to be relevant. Our results suggest that it may be unwise to assume that biodiversity patterns are at equilibrium and provide a first step in studying how these patterns unfold through time.     

Time-dependent solutions of the spatially implicit neutral model of biodiversity

Previous research into the neutral theory of biodiversity has focused mainly on equilibrium solutions rather than time-dependent solutions. Understanding the time-dependent solutions is essential for applying neutral theory to ecosystems in which time-dependent processes, such as succession and invasion, are driving the dynamics. Time-dependent solutions also facilitate tests against data that are stronger than those based on static equilibrium patterns. Here I investigate the time-dependent solutions of the classic spatially implicit neutral model, in which a small local community is coupled to a much larger metacommunity through immigration. I present explicit general formulas for the eigenvalues, left eigenvectors and right eigenvectors of the models’s transition matrix. The time-dependent solutions can then be expressed in terms of these eigenvalues and eigenvectors. Some of these results are translated directly from existing results for the classic Moran model of population genetics (the Moran model is equivalent to the spatially implicit neutral model after a reparameterization); others of the results are new. I demonstrate that the asymptotic time-dependent solution corresponding to just these first two eigenvectors can be a good approximation to the full time-dependent solution. I also demonstrate the feasibility of a partial eigendecomposition of the transition matrix, which facilitates direct application of the results to a biologically relevant example in which a newly invading species is initially present in the metacommunity but absent from the local community.     

Testing the standard neutral model of biodiversity in lake communities

Hubbell's (2001) neutral model describes how local communities are structured if population dynamics are statistically identical among species in a constant, possibly patchy, environment with random speciation. Tests of this model have been restricted largely to terrestrial communities. Here we tested the fit of this neutral model to fish, zooplankton and phytoplankton species–abundance distributions from 30 well-studied lake communities varying widely in lake size and productivity. We measured the fit of the communities to the neutral model in three ways. All but two zooplankton (7 of 9) and all but three fish (9 of 12) communities were consistent with all three measures of fit. However, all nine phytoplankton communities did not fit the neutral model by at least one measure. This result for phytoplankton communities represents to date the most consistent failure of the standard neutral model to predict the shape of species-abundance distributions.     

The current limitations on tissue banking from an academic perspective

For research on human physiology and pathologies the most relevant results come from human tissue, necessitating the creation of more tissue banks. This need is acknowledged by academics, clinical researchers and the pharmaceutical industry. For academics, the major obstacles to establishing tissue banks are the somewhat cumbersome ethical procedures, a perceived lack of demand for human tissue and insufficient knowledge about supply and its demographic differences. The causes are inter-related: confusing and time-consuming ethics applications cause some researchers to avoid human tissue work and expend research efforts on animal studies, leading to a false presumption of a lower level of demand for human tissue. Lack of knowledge about why rates of donation are low, and why there are differences in donation for different organs, leads to an uncertainty about supply. This too poses a problem for tissue bank establishment, and further research into this area is required.