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.     

Speciation and the neutral theory of biodiversity

The neutral theory of biodiversity purports that patterns in the distribution and abundance of species do not depend on adaptive differences between species (i.e. niche differentiation) but solely on random fluctuations in population size (“ecological drift”), along with dispersal and speciation. In this framework, the ultimate driver of biodiversity is speciation. However, the original neutral theory made strongly simplifying assumptions about the mechanisms of speciation, which has led to some clearly unrealistic predictions. In response, several recent studies have combined neutral community models with more elaborate speciation models. These efforts have alleviated some of the problems of the earlier approaches, while confirming the general ability of neutral theory to predict empirical patterns of biodiversity. However, the models also show that the mode of speciation can have a strong impact on relative species abundances. Future work should compare these results to diversity patterns arising from non‐neutral modes of speciation, such as adaptive radiations.     

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.     

Hubbell's fundamental biodiversity parameter and the Simpson diversity index

Central to Hubbell's neutral theory of biodiversity is a universal, dimensionless fundamental biodiversity parameter that is the product of community size and speciation rate. One of the most important discoveries of Hubbell's theory is that the species‐abundance distribution and the species–area relationship of the neutral metacommunity is completely determined by this fundamental biodiversity parameter, although the diversity patterns of the local community are collectively determined by the biodiversity parameter and migration. Using the relative abundance of species and following the concept of heterozygosity of population genetics, here we developed an analytical relationship between this biodiversity parameter and the well‐known Simpson diversity index. This relationship helps bridge the evolutionary aspect of biodiversity to the ecological and statistical aspect of the diversity. The relationship between these two parameters suggests that diversity patterns of the metacommunity can also be equally described by the Simpson index. This relationship provides an alternative approach to interpret and estimate the fundamental biodiversity parameter for the metacommunity.     

Age structure in neutral theory resolves inconsistencies related to reproductive-size threshold

Aims Neutral theory consists of a suite of models that assume ecological equivalence among individual organisms. They have been most commonly applied to tropical forest tree communities either as null models or as approximations. Neutral models typically only include reproductive adults; therefore, fitting to empirical tree community data requires defining a reproductive-size threshold, which for trees is usually set arbitrarily to a diameter at breast height (DBH) of 100 mm. The inevitable exclusion of some reproductive adults and inclusion of some saplings cause a non-random sampling bias in neutral model fits. Here, we investigate this problem and resolve it by introducing simple age structure into a neutral model.Methods We compared the performance and sensitivity of DBH threshold of three variants of a spatially explicit neutral model: the traditional model, a model incorporating random sampling and a model with two distinct age classes—reproductive adults and saplings. In the age-structured model, saplings are offspring from adults that disperse according to a Gaussian dispersal kernel around the adults. The only extra parameter is the ratio of adults to saplings, which is not a free parameter but directly measurable. We used species–area relationships (SARs) to explore the predicted effect of saplings on the species richness at different scales in our model. We then evaluated the three model variations to find the parameters required to maintain the observed level of species richness in the 50-ha plot on Barro Colorado Island (BCI). We repeated our analysis filtering the data at different minimum tree-size thresholds in order to find the effect this threshold has on our results. Lastly, we used empirical species–individual relationships (SIRs) to test the pre-existing hypothesis that environmental filtering is the primary cause of differences between the assemblage of saplings and that of adults on BCI.Important findings Our age-structured neutral model was characterized by SARs that were insensitive to the presence of saplings at large scales and highly sensitive to them at small scales. Both models without age structure were highly sensitive to the DBH threshold chosen in a way that could not be explained based on random samplings alone. The age-structured neutral model, which allowed for non-random sampling based on life stage, was consistent with species richness observations. Our analysis of empirical SIRs did not support environmental filtering as a dominant force, but it did show evidence for other differences between age classes. Age can now be easily incorporated into future studies of neutral models whenever there is a concern that a sample is not entirely composed of reproductive adult individuals. More generally, we suggest that modeling studies using tree data subject to a minimum size threshold should consider the sensitivity of their results to that threshold.     

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.     

Linking dispersal, immigration and scale in the neutral theory of biodiversity

In the classic spatially implicit formulation of Hubbell's neutral theory of biodiversity a local community receives immigrants from a metacommunity operating on a relatively slow timescale, and dispersal into the local community is governed by an immigration parameter m . A current problem with neutral theory is that m lacks a clear biological interpretation. Here, we derive analytical expressions that relate the immigration parameter m to the geometry of the plot defining the local community and the parameters of a dispersal kernel. Our results facilitate more rigorous and extensive tests of the neutral theory: we conduct a test of neutral theory by comparing estimates of m derived from fits to empirical species abundance distributions to those derived from dispersal kernels and find acceptable correspondence; and we generate a new prediction of neutral theory by investigating how the shapes of species abundance distributions change theoretically as the spatial scale of observation changes. We also discuss how our main analytical results can be used to assess the error in the mean-field approximations associated with spatially implicit formulations of neutral theory.     

A novel genealogical approach to neutral biodiversity theory

Current neutral theory in community ecology views local biodiversity as a result of the interplay between speciation, extinction and immigration. Simulations and a mean‐field approximation have been used to study this neutral theory. As simulations have limitations of convergence and the mean‐field approximation ignores dependencies between species’ abundances when applied to species‐abundance data, there is still no final conclusion whether the neutral theory or the traditional lognormal model describes community structure best. We present a novel analytical framework, based on the genealogy of individuals in the local community, to overcome the problems of previous approaches, and show, using Bayesian statistics, that the lognormal model provides a slightly better fit to the species‐abundance distribution of a much‐discussed tropical tree community. A key feature of our approach is that it shows the tight link between genetic and species diversity, which creates important perspectives to future integration of evolutionary and community ecological theory.     

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.     

Community structure of vascular epiphytes: a neutral perspective

Vascular epiphytes form a diverse group of almost 30 000 species, yet theory concerning their community structure is still largely lacking. We therefore employed the simplest models of biodiversity, (near-)neutral models, to generate hypotheses concerning their community structure. With recently developed tools for (near-)neutral models we analyzed species abundance data from many samples in Central and South America which we divided into four metacommunities (Mesoamerica, Central America, Amazonia and Paraná), where for each metacommunity we considered two subsets differing in dispersal syndrome: an animal-dispersed guild and a wind-dispersed guild. We considered three models differing in the underlying speciation mode. Across all metacommunities, we found observed patterns to be indistinguishable from patterns generated by neutral or near-neutral processes. Furthermore, we found that subdivision in different dispersal guilds was often supported, with recruitment limitation being stronger for animal-dispersed species than for wind-dispersed species. This is the first time that (near-)neutral theory has been applied to epiphyte communities. Future efforts with additional data sets and more refined models are expected to further improve our understanding of community structure in epiphytes and will have to test the generality of our findings.     

Inferring the parameters of the neutral theory of biodiversity using phylogenetic information and implications for tropical forests

We develop a statistical method to infer the parameters of Hubbell's neutral model of biodiversity using data on local species abundances and their phylogenetic relatedness. This method uses the approximate Bayesian computation (ABC) approach, where the data are summarized into a small number of informative summary statistics. We used three statistics: the number of species in the sample, Shannon H index of evenness and Shao and Sokal's B ₁ index of phylogenetic tree imbalance. Our approach was found to outperform previous methods, illustrating the potential of ABC methods in ecology. Applying it to four large tropical forest tree data sets, the best-fit immigration rates m were found to be two orders of magnitude smaller and regional diversities θ larger than previously reported for the same data. This implies that neutral-compatible regional pools of tropical trees should extend over continental scales, and that m measures, in this context, mostly the frequency of long-distance dispersal events.     

Applying ecological models to communities of genetic elements: the case of neutral theory

A promising recent development in molecular biology involves viewing the genome as a mini‐ecosystem, where genetic elements are compared to organisms and the surrounding cellular and genomic structures are regarded as the local environment. Here, we critically evaluate the prospects of ecological neutral theory (ENT), a popular model in ecology, as it applies at the genomic level. This assessment requires an overview of the controversy surrounding neutral models in community ecology. In particular, we discuss the limitations of using ENT both as an explanation of community dynamics and as a null hypothesis. We then analyse a case study in which ENT has been applied to genomic data. Our central finding is that genetic elements do not conform to the requirements of ENT once its assumptions and limitations are made explicit. We further compare this genome‐level application of ENT to two other, more familiar approaches in genomics that rely on neutral mechanisms: Kimura's molecular neutral theory and Lynch's mutational‐hazard model. Interestingly, this comparison reveals that there are two distinct concepts of neutrality associated with these models, which we dub ‘fitness neutrality’ and ‘competitive neutrality’. This distinction helps to clarify the various roles for neutral models in genomics, for example in explaining the evolution of genome size.     

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.     

A neutral theory with environmental stochasticity explains static and dynamic properties of ecological communities

Understanding the forces shaping ecological communities is crucial to basic science and conservation. Neutral theory has made considerable progress in explaining static properties of communities, like species abundance distributions (SADs), with a simple and generic model, but was criticised for making unrealistic predictions of fundamental dynamic patterns and for being sensitive to interspecific differences in fitness. Here, we show that a generalised neutral theory incorporating environmental stochasticity may resolve these limitations. We apply the theory to real data (the tropical forest of Barro Colorado Island) and demonstrate that it much better explains the properties of short‐term population fluctuations and the decay of compositional similarity with time, while retaining the ability to explain SADs. Furthermore, the predictions are considerably more robust to interspecific fitness differences. Our results suggest that this integration of niches and stochasticity may serve as a minimalistic framework explaining fundamental static and dynamic characteristics of ecological communities.