The dominant–egalitarian transition in species-rich communities

  1. David A Kessler  Is a corresponding author
  2. Nadav M Shnerb  Is a corresponding author
  1. Department of Physics, Bar-Ilan University, Ramat-Gan, Israel

Abstract

Diverse communities of competing species are generally characterized by substantial niche overlap and strongly stochastic dynamics. Abundance fluctuations are proportional to population size, so the dynamics of rare populations is slower. Hence, once a population becomes rare, its abundance gets stuck at low values. Here, we analyze the effect of this phenomenon on community structure. We identify two distinct phases: a dominance phase, in which a tiny number of species constitute most of the community, and an egalitarian phase, where it takes a finite fraction of all species to constitute most of the community. Using empirical data from microbial, planktonic, and macroorganismal systems, we demonstrate the relevance of this transition and show how demographic stochasticity and immigration critically determine phase behavior. Our results suggest that even slight changes in noise strength or immigration rates can lead to abrupt shifts in community diversity.

Editor's evaluation

This valuable paper explores the question of when complex ecosystems will be dominated by a few species. The authors present compelling, general arguments for a phase transition from what they call a dominance phase (few species dominate biomass) to an egalitarian phase (no species dominates the biomass).

https://doi.org/10.7554/eLife.103999.sa0

Introduction

Many ecological systems contain numerous competing species, strains, or types (Hutchinson, 1961; Stomp et al., 2011; ter Steege et al., 2013; Connolly et al., 2014; Fierer et al., 2007). In such systems, it can be assumed a priori that the structure of the community essentially reflects resource competition among species (Gause, 2003; Tilman, 1982), with factors such as niche overlap and fitness differences playing a central role in the community assembly (Chesson, 2000; Chesson, 2003). Analyzing these factors and their impact is crucial for understanding the dynamics of these systems, and consequently, for our ability to intervene in these dynamics to achieve desired outcomes – from maintaining biodiversity in a changing world to successfully altering the state of the gut microbiome (David et al., 2014; Grilli, 2020; Eguíluz et al., 2019; Cooper et al., 2024; Callaghan et al., 2021).

Unfortunately, progress in this area has been quite difficult. The coexistence of many species remains largely puzzling, especially given the competitive exclusion principle (Tilman, 1982) and the strict constraints on the stability and feasibility of such complex systems (May, 1972). Moreover, quantifying the relevant parameters in diverse communities is extremely challenging, particularly considering the high level of stochasticity typically present in ecological dynamics. As a result, approaches inspired by statistical physics, which examine generic models through a few summary statistics, have gained significant popularity in recent years (Fisher and Mehta, 2014; Kessler and Shnerb, 2015; Bunin, 2017; Barbier et al., 2018; Grilli, 2020; van Nes et al., 2024).

Broadly speaking, attempts to present a generic analysis of diverse communities can be divided into two main approaches that differ in their interpretation of what determines dominance. In one class of models, the identity of the high-abundance species reflects intrinsic fitness differences (whether through growth rates, competitive advantage, or resource-use efficiency) so that the dominant species are effectively the ‘fittest’ (Fisher and Mehta, 2014; Bunin, 2017; Barbier et al., 2018; Azaele and Maritan, 2023; Marcus et al., 2022). In contrast, an alternative view holds that dominance arises from contingent stochastic dynamics: species happen to be abundant not because they are intrinsically superior, but as a result of random historical trajectories. These trajectories may reflect demographic noise, environmental stochasticity, or even deterministic chaos in systems with complex interactions.

The prototype of this second class of approaches was introduced in the neutral model proposed by Crow and Kimura, 1970 and Hubbell, 2001. In the original neutral models, demographic stochasticity – the inherent randomness in the birth–death process at the individual level – is the sole driver of abundance variations. Analytical solutions to these models are relatively easy to obtain (Volkov et al., 2003; Azaele et al., 2015) and have been quite successful in explaining the observed species abundance distributions (SADs), as well as other static patterns, in both regional and local communities (Rosindell et al., 2011). However, pure demographic stochasticity cannot account for dynamic patterns, whether evolutionary (such as the time to the most recent common ancestor) (Nee, 2005; Ricklefs, 2006) or ecological (such as the dynamics of abundance variation and similarity indices) (Leigh, 2025; Chisholm and O’Dwyer, 2014; Kalyuzhny et al., 2014b; Kalyuzhny et al., 2014a). Demographic stochasticity results in relatively slow and weak abundance variations, whereas observed variability is much stronger and occurs more rapidly (Kalyuzhny et al., 2014b; Chisholm et al., 2014).

The time-averaged neutral model (Kalyuzhny et al., 2015) addresses these limitations by relaxing the assumption of time-independent fitness. In this model, the relative fitness of each species fluctuates over time, but all species share the same mean fitness. The analysis of the time-averaged neutral model is more complex, as environmental stochasticity can facilitate coexistence through the storage effect (Chesson and Warner, 1981). However, the significance of this effect diminishes in highly diverse systems (Dean et al., 2017; Danino and Shnerb, 2018; Pande and Shnerb, 2022; Meyer et al., 2022).

Recent studies have highlighted a notable effect of environmental stochasticity in competitive communities, known as ‘stickiness’ (van Nes et al., 2024) or ‘diffusive trapping’ (Dean and Shnerb, 2020). Under environmental stochasticity, abundance variations are proportional to population size, meaning that the dynamics of rare species is slow, causing them to linger near the extinction point for extended periods. van Nes et al., 2024 suggested that this stickiness enables time-averaged neutral models to produce patterns – such as changes in abundance over time, dominant species turnover, and community egality – that closely resemble those observed in real ecological communities. Similar results were found by Mallmin et al., 2024 in a system of competing species for which the overall dynamics are chaotic. This makes sense, given that the relative fitness of a specific species depends on the abundances of its competitors, so if these abundances fluctuate strongly over time, as they do in the chaotic phase, the community will likely reach a state that resembles time-averaged neutrality.

Our main goal in this article is to identify a previously unrecognized phase transition between dominant and egalitarian communities. To describe this transition, we define the number of dominant species, S1/2, as the minimum number of species that must be grouped together at a given moment to account for more than half of the total population. The level of equality in the community is then measured by the egality ratio, S1/2/S, where S is the total species richness. We show that in the dominance phase, the growth of S1/2 with S is sublinear, so the egality parameter approaches zero as S, whereas in the egalitarian phase, S1/2 grows linearly with S, and the egality parameter converges to a finite value. We demonstrate the existence of these two phases through analysis of empirical data on patterns of commonness and rarity across a wide range of ecological systems, from tropical forest trees to the human microbiome.

During the analysis and comparisons with the empirical data, we also arrived at two important insights.

  • The basic model presented by van Nes et al., 2024 includes only growth, competition, and environmental stochasticity. We show that this dynamics alone leads to monodominance, where all populations except one decline below any finite abundance level. Coexistence requires a mechanism, such as immigration, that imposes a lower bound on species abundances. Therefore, immigration – which was considered only in the supplementary materials of van Nes et al., 2024 – is a crucial ingredient, and so even small changes in the immigration rate can shift the system between ecological phases.

  • An additional factor not included in the model of van Nes et al., 2024 is demographic stochasticity. We will show that for populations of macroorganisms, the abundance distribution of relatively rare species is significantly affected by demographic stochasticity, and to properly describe it, this mechanism must also be taken into account.

The dominant–egalitarian phase transition has practical implications. A small decrease in the immigration rate or, alternatively, a small increase in the level of stochasticity can lead to a sharp decline of the diversity within a given community. Moreover, in the section ‘Materials and methods’, we have examined cases in which neutrality is broken, including fitness differences and asymmetry in the competitive interaction matrix. The results of this study demonstrate that the phase transition between the egalitarian and dominance regimes is not an artifact of perfect neutrality, but remains robust over a finite region of parameter space near the neutral point.

Results

As we have already mentioned, without external immigration, the stickiness will eventually cause the abundance of all species, except one, to drop below any finite threshold, leading the community to a state of monodominance. To avoid disrupting the continuity of the main discussion, we detail and demonstrate this claim in ‘SI A’ (see also the supplement to Figure 1). Therefore, the basic model we study is given by

(1) dNidt=μ+Ni(1j=1SNjK)+ri(t)Nidridt=riτ+θηi(t),

where Ni, i=1,,S is the abundance of the ith species. K sets the overall carrying capacity. The growth rate of the ith species is 1+ri(t), where ri(t) is an independently fluctuating variable, with a correlation time τ, driven by the white noise ηi. The immigration rates of all species are identical and given by μ>0. Further details on the model are given in the ‘Materials and methods’ section.

Figure 1 with 1 supplement see all
Probability distribution of species abundance.

Red circles: results of a simulation of the process described in Equation 1 with parameters S=200, μ=0.6, K=2106, σ~e=0.1 and τ=8. The values of the extracted parameters for Equation 4 are c=0.0034 and σ~e2=0.0131. The solid line in blue is the expression of Equation 4 with parameters α=91.2 and β=1.56.

is an independently fluctuating variable, with a correlation time τ, driven by the white noise ηi. The immigration rates of all species are identical and given by μ>0. Further details on the model are given in the ‘Materials and methods’ section.

Stickiness and species abundance distribution

Without immigration, the carrying capacity parameter K plays no role in the dynamics. Actually, by defining Ni=Kni one may rescale Equation 1 to K=1. When immigration is introduced, the parameter K/μ expresses the ratio between the carrying capacity and the typical minimum level of population size. The larger this ratio, the stronger the stickiness effect.

On top of that, since stickiness arises from environmental stochasticity, it becomes stronger as σ~e2θ2τ2 increases. Dimensional analysis reveals that the dimensionless parameter that governs stickiness is

γKσ~e2μ.

The larger γ is, the stronger is the stickiness, and the community approaches monodominance when γ; for example, when μ0.

In order to advance in the analysis, we are interested in replacing Equation 1, which provides us with a description of an S-dimensional system where each species can affect every other species, with an effective, one-dimensional equation for a focal species, considering all others as a single rival species. In the classical neutral theory, with pure demographic stochasticity, this can be done trivially, as the species identity of a particular individual plays no role in the dynamics. In the time-averaged neutral theory, however, the situation is much more subtle.

The distinction between neutral and time-averaged neutral models – as it arises in the attempt to derive an effective one-species description – was clarified in Danino and Shnerb, 2018 and Steinmetz et al., 2020, and has to do with the difference between the time-averaged growth rate of a given species and the population-averaged growth rate of the community. In Equation 1, the linear growth rate of each species is 1+r(t). Since the r process is symmetric around zero, the time-average growth rate of each species is unity. Nevertheless, the instantaneous growth rate of the community is, on average, greater than 1. At any given moment, the fitter species are growing faster, so that on average more than 50% of the individuals belong to instantaneously beneficial species. As a result, the typical value of jNj>K, and therefore the dynamics of a single species satisfies the effective one-dimensional stochastic differential equation,

(2) dNdt=μ+(σe22c)N+σeη(t)N,

where c=E[(jNjK)/K]. The noise in the net growth rate of the ith species arises from two sources, the first being the direct effect of ηi(t), and the second being the fluctuations in the competition term, jNj/K. The two parameters c and σe may be measured, for any given values of K, S, and σ~e, from long simulations of Equation 1. Unless we are looking at the most abundant species, the correlations between the competition term and ηi are small, and the two contributions are roughly independent. Thus, the overall effect of these two contributions is

(3) σe2σ~e2+σc2

Here, σc2 is the time integral of the two-time correlation function of c and equals 2τVar[Ni/K]. For example, for σ~e=0.01, S=200, K=2106, μ=0.6, we have σc20.0031 and σe20.1031, in line with Equation 3. The Stratonovich term σe2/2 expresses the fact that for a population that sometimes grows exponentially and sometimes declines, even if its average growth rate is zero, the arithmetic mean still increases over time. Once this term is introduced, the standard Ito calculus may be applied to Equation 2.

Once the relevant parameters are calibrated, the distribution for P(N) may be extracted from Equation 2, see ‘SI 2’. The resulting distribution is

(4) P(n)=Aeα/nnβ,

where α=2μ/σe2, β=1+2c/σe2 and A is a normalization factor. Figure 1 illustrates the success of the approximation and the validity of Equation 4 in a simulation of a community of S=200 interacting species governed by Equation 1.

A similar result was presented a few months ago by Mallmin et al., 2024, who dealt with a system of competing species where the community dynamics is chaotic (but without external stochasticity). In such a case, one can consider, for each focal species, all other species as an effective external environment whose fluctuations generate stochasticity in the instantaneous growth rate of the focal species. A discussion of the similarities and differences between our case and the chaotic model will be presented below.

The egalitarian transition

The most striking implication of Equation 4 is the sharp shift in the compositional properties of the community at β=2. When β>2, abundant species are relatively rare, resulting in an ‘egalitarian’ community. Conversely, when β<2, the community composition is dominated by a few exceptionally abundant species.

To quantify the egality of the community, one may use the criterion suggested by van Nes et al., 2024, which involves comparing the total number of species, S, to the minimal number of species required to make up half of the community, S1/2. In a more egalitarian community, the fraction S1/2/S is finite, indicating that the number of dominant species is proportional to the total number of species. In contrast, in a community dominated by only a few species, the ratio S1/2/S approaches zero as S increases, meaning the number of dominant species is sublinear in S.

In Appendix C, we provide the relevant mathematical analysis, demonstrating that, as β decreases towards 2, S1/2/S monotonically decreases, and in the limit β2+,

(5) S1/2Se(ln2)/(β2).

Thus, this egality parameter vanishes in a singular manner in that limit, indicating the transition out of the egalitarian phase.

The numerical results presented in Figure 2A illustrate this phenomenon and offer several additional interesting insights. Beyond the positive indication regarding the result in the limit where S, Figure 2A shows that the value of β depends only weakly on S so it can be considered approximately constant. Moreover, as seen in Figure 2B, for β<2 the dependence of S1/2 on S follows a power law, with an exponent approaching unity from below at the phase transition point β=2.

The egaility parameter as a function of the number of species.

Left panel (A): the egality parameter S1/2/S is plotted against 1/S, the total number of species, for σ~e=0.015 (blue) 0.02 (red) 0.04 (orange) and 0.08 (green). For all these simulations μ=0.6, τ=8 and K=104. The values of the parameter β, as extracted from the simulation, are plotted (again, vs. 1/S) in the inset. As expected, for β>2 (blue) the egality parameter S1/2/S extrapolates, at S, to a finite value, meaning that the fraction of high-abundance species is finite and the community is egalitarian. This feature is also reflected in the right panel (B), where at large S, S1/20.04S for σ~e=0.02 (red) and S1/20.08S for σ~e=0.015 (blue). On the other hand, when σ~e=0.08 (green), the value of β approaches 1.8 and is definitely smaller than two. In addition, the egality parameter extrapolates to zero at large S, meaning that the community is dominated by a small number of species. The results shown in the right panel suggest S1/2S0.3. The case σ~e=0.04 (orange) represents a near marginal case, for which β extrapolates to values slightly smaller than 2, and it would appear that S1/2/S is tending toward zero. Indeed, the right panel shows a power law with exponent smaller than unity and decreasing with increasing σ~e.

Of course, there are other quantities one could employ to quantify the egality. One such possibility is the inverse participation ratio, I, defined as I=[i(Ni/K)2]1. This varies from 1 in the limit of single species dominance to S for an egalitarian community. Our measurements (not shown) indicate that, similarly to S1/2, I increases with S, though the exponents controlling the growth are slightly lower. For example, for σ~e=0.08, IS0.26, as opposed to S1/2S0.29, and for σ~e=0.4, IS0.4, while S1/2S.52.

Comparison with empirical results and the impact of demographic stochasticity

Let us now examine the SAD in several cases of diverse communities, assess their degree of alignment with the results presented above, and attempt to differentiate between dominance and egalitarian communities, linking the results to the fundamental characteristics of each system.

We analyze the community structure in four systems: the human gut microbiome (David et al., 2014), marine prokaryotes (Eguíluz et al., 2019), tropical trees (Cooper et al., 2024), and bird species (Callaghan et al., 2021). All of these communities are hyperdiverse, with thousands of species, making them reflective of the limit where 1/S is very close to zero – a limit where the distinction between egalitarian and dominance systems is clearly defined, with the relevant values of S1/2/S being those shown in Figure 2.

Figure 3a shows the excellent agreement between the empirical data for the gut microbiome and our model prediction (Equation 4). In the case of ocean prokaryotes, presented in panel (b), only a pure power law is observed, perhaps because the sampling strength is insufficient. Weak sampling shifts the distribution leftward, thus hiding the true characteristics of the SAD for small abundance populations (Maruvka et al., 2010). In both cases β<2, indicating that these systems are in the dominant phase, consistent with the very small values of S1/2/S observed in the data.

Species abundance distributions, as obtained for a few diverse communities.

The gut microbial community (OTUs of subject A from David et al., 2014) fits our formula Equation 4 quite well. Immigration is extremely weak, but the sampling power is strong enough to reveal some of the effects of the decrease in the number of species at low densities due to immigration. The corresponding distribution for the oceanic prokaryote population Eguíluz et al., 2019 is very close to a pure power law, possibly because the sampling is not deep enough (see text). In both cases, β<2, indicating that the community is dominated by only a few species. The distributions for the global bird population (Callaghan et al., 2021) and tropical trees in the Amazon Basin (Cooper et al., 2024) show a clear transition between two power-law behaviors. As explained in the text, for these macroorganisms the effect of demographic stochasticity must be taken into account. When this is done (see ‘SI D’), we obtain an excellent fit for the results using the corrected formula. Equation 6. For birds, we find β2, placing them on the margin between egalitarian and dominance communities. For tropical trees in the Amazon Basin, β is definitely larger than two, indicating that the community is indeed egalitarian. Note the similarities between the values of S1/2/S in the empirical results and the numerical experiments in Figure 2. For more information, see ‘SI E’.

In contrast, in Figure 3c and d, for birds and Amazon basin trees, respectively, the empirical distribution does not fit Equation 4. Instead, the population size distribution shows a pronounced crossover between two descending power laws, without a maximum point at a finite abundance.

To explain this, we note that trees and birds differ from microorganisms in two important ways. First, since macroorganisms tend to have larger body sizes and longer lifespans, their metabolism buffers them against environmental fluctuations. Additionally, since populations of macroorganisms often span wide geographic ranges, most environmental variations appear uncorrelated across populations. Environmental variation that affects individuals, or small groups within the population, in an uncorrelated manner contributes to demographic stochasticity rather than to environmental stochasticity, which, by definition, affects entire populations coherently. Therefore, for macroorganisms we need to account for demographic stochasticity.

In ‘SI D’, we consider the case of a time-averaged neutral community with demographic stochasticity in addition to the previously considered immigration and environmental stochasticity. The expected distribution for species abundances now takes the form

(6) P(n)=An1+2μ/σd2(σd2σe2+n)2c/σe22μ/σd2,

with σd parameterizing the strength of demographic stochasticity. Equation 6 converges to Equation 4 for values of n satisfying nσd2/σe2 and predicts a different power-law for smaller values of n. With this new formula, one can now successfully fit the empirical data in Figure 3c and d. It is noteworthy that the transition point between the two power-laws happens for values of the abundance which are quite large, reinforcing our above comments about the effect of the large spatial and temporal scales of macroorganisms on reducing σe and increasing σd.

Even with the inclusion of demographic stochasticity, the distinction between dominant and egalitarian communities depends solely on the decay of the tail and its corresponding exponent β, which, for Equation 6, remains as β=12c/σe2. Figure 3 thus illustrates the three types of behavior demonstrated in the numerical experiments shown in Figure 2: the microorganism communities are in the dominance phase, the tropical tree community is in the egalitarian phase, and the bird community appears to fall in between.

A comparison between the values of the egality parameter in Figure 2 and the values of S1/2/S in van Nes et al., 2024 (for the same range of species richness S) also shows that some communities fall within the dominance phase, while others fall within the egalitarian phase. On the other hand, the results for plankton (Ser-Giacomi et al., 2018), in which β[1...2], appear to suggest that these microorganismal communities are all in the dominance phase.

Relaxing the assumption of neutrality

So far, we have analyzed the properties of a system that is fully neutral in the deterministic limit – that is, all species have exactly the same fitness. In this section, we aim to explore, or at least provide a glimpse of, the outcomes that arise when this assumption is relaxed; in particular, we assess the robustness of our analysis in the vicinity of the neutral point.

To break perfect neutrality, we introduce heterogeneity into the interaction matrix. Specifically, we replace Equation 1 by

(7) dNidt=μ+Ni(1Ni+ji(13ϵ+ϵαi,j)NjK)+ri(t)Nidridt=riτ+θηi(t),

where αi,j are zero-mean, unit-variance Gaussian random numbers, and αi,j is independent of αj,i, so the interaction matrix is asymmetric. For ϵ=0, this reduces to our original time-averaged neutral model. The additional 3ϵ term ensures that species compete most strongly against themselves, as is commonly assumed in ecological models, since the niche overlap between individuals of the same species is maximal.

The results are presented in Figure 4, which shows S1/2 as a function of ϵ for two different values of the environmental stochasticity σ~e. As ϵ0, the system may reside either in the egalitarian phase (if σ~e is small) or in the dominance phase under strong environmental stochasticity. In contrast, for large values of ϵ, S1/2 becomes independent of σ~e.

Figure 4 with 1 supplement see all
From stochasticity-driven to chaotic dynamics.

S1/2, as obtained from simulations of Equation 7, is plotted against ϵ, for N=200, K=2106 and μ=0.6. For σ~e=0.01 (blue circles), the system is in the egalitarian phase and the main effect of the transition to chaos is an effective enhancement of the strength of stochasticity that makes the system less egalitarian, hence S1/2 decreases. The inverse effect is observed for σe=0.08 (red squares). For large values of ϵ, the effect of σe becomes negligible and the dependence of S1/2 on σe weakens significantly. The data was obtained by averaging over at least 20 different realizations of the interaction matrix αi,j.

To understand these results, note that in the absence of environmental stochasticity, the heterogeneous and asymmetric interactions lead to chaotic dynamics of the trajectories, as discussed in Arnoulx de Pirey and Bunin, 2024. Illustrations of fluctuation dynamics dominated by stochasticity versus those dominated by deterministic chaos are provided in ‘SI F’. Under chaotic dynamics, the system supports, at any instant of time, a clique of abundant species. This clique is stable in the sense that if all the species within it are placed together in a community, they will reach an equilibrium state in which all of them persist. However, the clique is invadable, meaning that at any given moment, there exist species outside the clique that are capable of invading and reaching high abundance (Kessler and Shnerb, 2025). Such a clique gives rise to a ‘shoulder’ at high abundance values (see Figure 5), since at any given moment there are several species with high abundance. In contrast, the hierarchical structure under pure environmental stochasticity leads to a simple power law, as seen earlier in Figure 1.

lnP(n) vs.lnn for ϵ={0 (solid blue line),0.05 (dashed red line)}, σ~e=0.01.

The distribution for ϵ=0 is typical of the egalitarian phase with a large slope for large n. The distribution for ϵ=0.05, however, has the slower falloff characteristic of the dominance phase. In addition, there is a shoulder for very large n, which is due to the existence of multiple sets of several transient quasi-stable species, arising from the chaotic dynamics present in the absence of environmental stochasticity.

As ϵ increases from zero to a finite value, the neutral-stochastic system discussed throughout this article transforms into a system where the dominant factor is deterministic chaos. In the chaotic phase, the dynamics still appears stochastic, but not due to fluctuations in external parameters, but rather as a result of the inherent characteristics of the deterministic dynamics.

Figure 4 demonstrates that the transition to chaos affects S1/2 in opposite ways in the two phases. In the egalitarian phase, the main effect of chaos is to increase the effective stochasticity, leading to a decline of S1/2. In the dominance phase, by contrast, S1/2 increases. This is due to the ‘shoulder’ in the abundance distribution described in Figure 5, that is, because dominance in the chaotic regime is associated with a clique rather than a single species.

Discussion

Hyperdiverse communities, like those analyzed in this article, are extremely important and frequently occur in nature. However, quantifying their specific parameters is an impossible task. Therefore, the attempt to understand the factors that dictate the community structure in these systems requires the use of models, which should preferably be as generic as possible.

Broadly speaking, there are three generic scenarios for the dynamics of multispecies systems: those involving a stable clique of resident species; those in which the composition of the group of high-abundance species changes over time due to chaotic dynamics; and those assuming neutral dynamics.

In the first case, the set of high-abundance species is determined by their relative fitness (which reflects their intrinsic growth rate and interspecific interactions), which remains fixed over time. In the second, the composition of the dominant cliques varies intermittently over time, but this turnover is governed by deterministic factors, and at any given moment, only a specific set of low-abundance species is capable of exponential growth. Moreover, the system admits a characteristic timescale, the time required for one of these rare species to invade. This timescale depends logarithmically on the migration rate μ.

In the third case, discussed here, the identity of the high-abundance species at any given moment is a random outcome of environmental stochasticity. As in the chaotic case, there is no sharp separation between core species and low-abundance ones, and the fundamental timescale of the system is determined by the stickiness phenomenon described above.

As far as can be judged from the empirical SADs of Figure 2, it appears that they cannot be explained by the first scenario. Two of its main features are a gap between species with abundance of O(K) and those with abundance of O(μ), and a truncated Gaussian-shaped SAD for the high-abundance species – both of which are not observed in our distributions.

The chaotic case yields SADs that more closely resemble those we have reviewed here: they exhibit no gap and display a power-law behavior. In fact, when attempting to derive an effective equation for a single species, one arrives, as noted above, at the very same equation. As is well known, it is often possible to identify multiple underlying models that give rise to the same abundance distribution. Therefore, matching the SAD is a necessary condition for considering a model a plausible candidate, but leaves open the possibility of alternative explanations.

To tell apart chaos from environmental stochasticity, it may be necessary to focus specifically on the dynamics of the most abundant species. If the noise is intrinsic, originating from chaotic fluctuations, then a species that has reached dominance is expected to exhibit fluctuation patterns that differ markedly from those observed when the same species is rare. In contrast, if the fluctuations stem from external factors such as weather or precipitation, we would expect the statistical properties of abundance fluctuations to be independent of the species’ current abundance level, as has been observed empirically (Kalyuzhny et al., 2014b).

Another test that may help distinguish between a chaotic and a stochastic system relates to the emergence of a fundamental timescale (Arnoulx de Pirey and Bunin, 2024) in chaotic dynamics – one associated with the invasion of rare species. Such a timescale does not appear in systems driven by external stochasticity. Like the previous approach, this too requires one to examine dynamic properties of the system, rather than relying solely on static features such as the abundance distribution.

Whatever the underlying mechanism, the net result, as we have demonstrated herein, is that in the limit of strong interspecies coupling, a transition occurs between an egalitarian phase and a dominance phase. When the immigration rate μ is reduced, total carrying capacity K increases, or when environmental stochasticity σe increases, the system can suddenly lose a significant amount of diversity at the transition point. The dependence on the total population size K is particularly interesting and has to do with the dispute about the relationships between productivity and species richness (Waide et al., 1999; Kadmon and Benjamini, 2006).

We will now discuss some potential extensions of the model described here and explore their possible implications.

Our model assumed an uncorrelated response of species to environmental variations. This treatment can be easily extended to include correlated responses, using techniques similar to those implemented by Loreau and de Mazancourt, 2008. In general, under correlated responses, the effective number of species in the community decreases.

Another interesting point relates to the interplay between the stickiness mechanism, through which environmental stochasticity causes populations to spend long periods in a state of rarity, and mechanisms such as the storage effect or relative nonlinearity (Chesson, 1982; Chesson, 2000; Ellner et al., 2016; Letten et al., 2018) that allow rare populations to invade due to environmental stochasticity. It is likely that these mechanisms weaken as the number of species increases (for storage, this has been demonstrated in several studies; Chesson and Huntly, 1989; Pande and Shnerb, 2022), and therefore, at least in diverse communities of competing species, the dominant effect will actually be that of stickiness.

The time-averaged neutral dynamic considered here assumes that the differences in average fitness between species are negligible, at least to a first approximation. This assumption is necessary in situations where niche overlap is large; otherwise, fitness differences would cause the extinction of most species. The justification for this can come from processes leading to emergent neutrality (Holt, 2006; Vergnon et al., 2012), or from the fact that environmental stochasticity itself is also a mechanism that ‘neutralizes’ fitness differences (Pande and Shnerb, 2022). Extending our model as we have done above to include weak deviations from time-averaged neutrality, perhaps using the dynamic mean-field approximation (Roy et al., 2019) , is a critical first step, which remains to be explored in more detail.

Species coexistence has long been, and remains, a theoretical puzzle of immense importance for understanding the dynamics of biological systems, with far-reaching practical implications. The research conducted in recent years has provided powerful theoretical tools that allow us to focus the discussion and understand the generic implications of community structure and the nature of species interactions on the range of possible outcomes. We believe that the parameter range we have addressed in this article – (nearly-) neutral dynamics and significant environmental stochasticity – is relevant for a wide variety of ecological systems, and we hope that our work will serve as a foundation for further studies that explain the wide range of diversity levels (for example, between a tropical forest and a tundra, or between different types of microbiome) in relation to the phase transition described here.

Materials and methods

Data compilation and analysis

Request a detailed protocol

As explained below, the phase transition we analyze in this paper manifests itself in two characteristics: the SAD and the egality parameter S1/2/S. Like any phase transition in complex systems, the distinction between the two phases becomes sharper as the system grows larger, so that S increases. To test our predictions, we used several recent large databases on communities with thousands of species: human microbiome (David et al., 2014), marine prokaryotes (Eguíluz et al., 2019), tropical forest trees (Cooper et al., 2024), and birds (worldwide) (Callaghan et al., 2021).

The (noise-free) neutral model

Request a detailed protocol

We consider a community of S populations, with differential responses to environmental variations. Ni is the abundance of the ith species, and the carrying capacity parameter is denoted by K. With no stochasticity and no immigration, Ni satisfies the equation

(8) dNidt=Ni(1j=1SNjK).

As all species admit the same growth rate and interact in a symmetric manner, the model is neutral. In particular, any solution that satisfies j=1SNj=K is a steady state of the system.

The time-averaged neutral model

Request a detailed protocol

To allow for nontrivial dynamics, one would like to add environmental stochasticity and immigration to the process described in Equation 8. The rate of immigration is denoted by μ, and the growth rate of each species fluctuates in time, undergoing an Ornstein–Uhlenbeck process with a correlation time τ. The corresponding set of stochastic differential equations, as set out in Equation 1.

(9) dNidt=μ+Ni(1j=1SNjK)+ri(t)Nidridt=riτ+θηi(t),

where ηi(t) is a white-noise process of unit strength. The index i of ηi indicates that each of the S species responds to the environmental variations in an independent manner. The correlation function of ri is

(10) ri(t)rj(t)=δi,jθ2τ2e|tt|/τ.

The equivalent white-noise strength of the environmental stochasticity, twice the diffusion constant of logNi, is then given by σ~e2=θ2τ2, the time integral of the ri correlation function.

Fokker–Planck (FP) equations

Request a detailed protocol

The main analytical results presented in this article are based on solving FP equations. To derive these equations, we assumed a relatively short correlation time τ for the environmental stochasticity, which allows the use of a single effective equation for the process described in Equation 8. The relevant considerations are detailed in ‘SI B’.

Demographic stochasticity

Request a detailed protocol

This refers to the intrinsic randomness of the birth–death process at the individual level. This stochasticity also manifests in erratic abundance variations, but the intensity of these fluctuations is weaker (compared to environmental stochasticity) when the population is large, and therefore it was not included in Equation 8 or in the analyses of van Nes et al., 2024 and Mallmin et al., 2024. However, demographic stochasticity is crucially important in small populations and in extinction processes; as we have seen, it must be considered outside the microorganism realm. To include demographic stochasticity in the treatment, one must add a term to the relevant equation that expresses noise whose amplitude scales with the square root of the population size, as opposed to environmental stochasticity, whose amplitude scales linearly with population size. The technical details can be found in ‘SI D.

Appendix 1

Supporting information

SI A: The necessity of immigration

van Nes et al., 2024 have considered various models that were proposed for a diverse community in which all species have the same mean fitness. Most cases were reviewed in the supplemental material of the aforementioned article, but the version analyzed in detail in the main text is

(SA1) dNidt=Ni(1j=1SNjK)dt+σ~eηi(t)Ni,

where η(t) is a white noise process. Ito calculus is applied, so the numerical integration procedure satisfies

(SA2) Ni,t+Δt=Ni,t+Δt[μ+Ni,t(1j=1SNj,tK)]+σ~eΔtN(0,1)Ni,t,

where N(0,1) is a random variable drawn (independently for each Δt and each species) from a normal distribution with zero mean and unit variance.

However, in the absence of an immigration term, that is, when μ=0, such a process leads to monodominance, the state where only one species survives and the abundance of all other species drops below any finite value. Figure 1—figure supplement 1 demonstrates numerically this behavior. The results presented in our main text explain this phenomenon. First, if μ=0 then α=0 in our solution, Equation 4. The resulting species abundance distribution diverges at zero abundance and cannot be normalized. The interpretation of such non-normalizable abundance distribution is that the abundance of all species will drop below any finite threshold as time goes to infinity.

One may arrive at the same conclusion from a different perspective (Dean and Shnerb, 2020). Let’s assume we have two species in a community with a fixed and rigid carrying capacity (i.e., they are playing a zero-sum game) and environmental stochasticity that determines their relative fitness.

Denoting by x the fraction of the focal species, its dynamics satisfies dx/dt=s(t)x(1x), where s(t) reflects the effect of the fluctuating environment. Therefore, the population undergoes a balanced walk in the logit space z=ln[x/(1x)]. In this type of random walk, a species that ‘gets stuck’ in the low-density region will remain there for increasing periods, as described in Section 5 of Dean and Shnerb, 2020. This is the ‘diffusive trapping’, or the stickiness of low-abundance states, that causes the system to be dominated at any given moment by a single species.

Now we can easily generalize this argument for S species: the relative fitness of each species (i.e., its fitness relative to the average fitness of the community) is a zero-mean random variable. Therefore, in the log-abundance space, it undergoes a random walk, and so all species except one will fall below any threshold of abundance over time. Thus, as seen in Figure 1—figure supplement 1, the system in the long run is dominated by a single species, and the abundance of all other species is negligible.

SI B: Derivation of Equation 4

The FP equation that corresponds to Equation 2 is

(SB1) P(n,t)n=σ222n2(n2P(n))+(cσ22)n(nP(n))μP(n)n

In steady state, P(n) is time-independent and so its time derivative vanishes, leaving one to solve the exact ordinary differential equation (after first integration),

(SB2) σ22n(n2P(n))+(cσ22)(nP(n))μP(n)=0.

Here, we have put the integration constant to zero to avoid diverging solutions that yield non-normalizable P(n) functions. Therefore,

(SB3) σ22n2P(n)n+(c+σ22)nPμP(n)=0.

That yields Equation 4 of the main text.

SI C: From egalitarian community to hyperdominant species

Let us consider a community of J individuals, whose species richness is S. By definition,

(SC1) S=0P(n)dn,

and

(SC2) J=0nP(n)dn,

From Equation 4 of the main text, we take P(n) and plug it into Equation SC1. First, we implement the condition on S to determine the normalization factor A, and find

(SC3) P(n)=Sαβ1Γ[β1]eα/nnβ

Now let us assume β>2. In that case, the integral in (S2) converges so that

(SC4) J=αβ2S.

As a metric for the egality of the community, let us take the one used in van Nes et al., 2024. First, we will check the richness S1/2 of the minimum set of species required to contain at least 50% of the entire community population. If the species within this set with the smallest abundance is represented by n1/2 individuals, this implies (implementing Equation SC4).

(SC5) n1/2nP(n)dn=αS((β1)Γ[β2,α/n1/2]Γ[β1])=J/2=α2(β2)S,

When μ0 then α0 and hence,

(SC6) Γ[β2,α/n1/2]Γ[β2]1β2(αn1/2)β2

therefore

(SC7) αn1/2=(Γ[β1]2)1/(β2),

so the ratio α/n1/2 approaches zero as β2+

The fraction of species whose abundance is equal or greater than n1/2 is obtained as

(SC8) S1/2S=1Sn1/2P(n)dn=1Γ[β1,α/n1/2]Γ[β1]=1Γ[β](Γ[β1]2)β1β2.

In particular, as β2+,

(SC9) S1/2Se(ln2)/(β2).

SI D: Demographic stochasticity

Every group of living organisms is influenced by various types of stochasticity. The literature (Kalyuzhny et al., 2014b; Lande and Engen, 2003) usually distinguishes between environmental stochasticity, which affects populations coherently, and demographic stochasticity, which affects the fitness of individuals or small groups of individuals in an incoherent manner, so that the overall reproductive success of the population reflects the sum of many random events. Weather (such as temperature or precipitation) fluctuations that are correlated over a large area and simultaneously affect entire populations contribute to environmental stochasticity, while the degree of success or failure of individuals in finding food or avoiding predators relative to its population mean contributes to demographic stochasticity.

Environmental stochasticity is represented in equations by noise that is proportional to the population size. Thus, in Equation 1 of the main text, the term ri(t), which varies stochastically, multiplies Ni. This mathematical expression reflects the situation where, in a good year, the average number of offspring for each individual increases, and in a bad year, it decreases. Demographic stochasticity, on the other hand, represents the success or failure of individuals, so that the overall result, according to the Central Limit Theorem, is fluctuations that are proportional to the square root of the population size. The corresponding term in the Langevin equation is η(t)Ni.

In an FP equation, the term that represents demographic stochasticity is

(SD1) σd222n2(nP(n))

We denote the strength of demographic stochasticity (i.e., the variance in the number of offspring per individual) by σd2, distinguishing it from σe2, which represents the strength of environmental stochasticity.

The extended version of Equation SD1 thus takes the form

(SD2) P(n,t)n=σe222n2(n2P(n))+σd222n2(nP(n))+(cσ22)n(nP(n))μP(n)n

Again, in steady state, P(n) is time-independent,

(SD3) σe22n(n2P(n))+σd22n(nP(n))+(cσe22)(nP(n))μP(n)=0,

and therefore the species abundance distribution is given by

(SD4) P(n)=An1+2μ/σd2(σd2σ2+n)2c/σ22μ/σd2

SI E: Figure 3 information

For each dataset, the data was log-binned (scaled by the inverse of the width of the bin), then the log-log histogram was fitted with the relevant expression.

Appendix 1—table 1
Summary of fitted models and parameter estimates (with 1σ uncertainties) for the four systems analyzed.
System (panel)Fitted equationFitted parameters (±1σ)
Gut microbiome (a)Fit to Equation 4 with
α=2μ/σe2,
β=1+2c/σe2, A normalization factor
A=1.77±0.060;
α=5.5×106±5×107;
β=1.31±0.01
Plankton (b)Fit to Equation 4 with
α=2μ/σe2,
β=1+2c/σe2,
A normalization factor
A=6.7±0.2;
α=0.8±0.5;
β=1.7±0.03
Birds (c)Fit to Equation 6 with
γ=2μ/σd2,
β=1+2c/σe2,
k=σd2/σe2,
A normalization factor
A=23±1;
γ=0.34±0.01;
β=2.05±0.05;
k=3.0×106±6×105
Amazon trees (d)A=14±1;
γ=2μ/σd2,
β=1+2c/σe2,
k=σd2/σe2,
A normalization factor
A=14±1;
γ=0.013±0.096;
β=2.4±0.2;
k=380±200

SI F: Stochastic and chaotic dynamics

This figure shows the abundance of 200 species over time, in the purely stochastic case (ϵ=0) and in the case dominated by chaotic dynamics (ϵ=0.05). We observe that chaotic dynamics introduce a new element into the system: a clique of species that are “dominant” at any given moment. Although there is always a species that can invade this clique and drive one or more of its members to low abundance levels, this process takes time. As a result, the number of highly abundant species is larger compared to the stochastic case, where we see a strong dominance of a single species at each moment.

Data availability

The current manuscript is a computational study, so no data have been generated for this manuscript. The simulation code can be accessed at https://github.com/kesslerda11/Sticky (copy archived at Kessler, 2024). The data analyzed was obtained from presvious published works, as indicated in the References.

The following previously published data sets were used
    1. David LA
    2. Materna AC
    3. Friedman J
    4. Campos-Baptista MI
    5. Blackburn MC
    6. Perrotta A
    7. Erdman SE
    8. Alm EJ
    (2014) European Nucleotide Archive
    ID PRJEB6518. Host lifestyle affects human microbiota on daily timescales - Subject A gut counts.
    1. Cooper D
    2. Lewis S
    3. ter Steege H
    4. Sullivan M
    5. Slik JWF
    6. Barbier N
    (2024) figshare
    Cooper consistent patterns of common species across tropical forest tree communities.
    https://doi.org/10.6084/m9.figshare.21670883.v1

References

  1. Conference
    1. Chesson PL
    (1982)
    The storage effect in stochastic competition models
    Mathematical Ecology: Proceedings, Trieste.
    1. Cooper DLM
    2. Lewis SL
    3. Sullivan MJP
    4. Prado PI
    5. Ter Steege H
    6. Barbier N
    7. Slik F
    8. Sonké B
    9. Ewango CEN
    10. Adu-Bredu S
    11. Affum-Baffoe K
    12. de Aguiar DPP
    13. Ahuite Reategui MA
    14. Aiba SI
    15. Albuquerque BW
    16. de Almeida Matos FD
    17. Alonso A
    18. Amani CA
    19. do Amaral DD
    20. do Amaral IL
    21. Andrade A
    22. de Andrade Miranda IP
    23. Angoboy IB
    24. Araujo-Murakami A
    25. Arboleda NC
    26. Arroyo L
    27. Ashton P
    28. Aymard C GA
    29. Baider C
    30. Baker TR
    31. Balinga MPB
    32. Balslev H
    33. Banin LF
    34. Bánki OS
    35. Baraloto C
    36. Barbosa EM
    37. Barbosa FR
    38. Barlow J
    39. Bastin JF
    40. Beeckman H
    41. Begne S
    42. Bengone NN
    43. Berenguer E
    44. Berry N
    45. Bitariho R
    46. Boeckx P
    47. Bogaert J
    48. Bonyoma B
    49. Boundja P
    50. Bourland N
    51. Boyemba Bosela F
    52. Brambach F
    53. Brienen R
    54. Burslem D
    55. Camargo JL
    56. Campelo W
    57. Cano A
    58. Cárdenas S
    59. Cárdenas López D
    60. de Sá Carpanedo R
    61. Carrero Márquez YA
    62. Carvalho FA
    63. Casas LF
    64. Castellanos H
    65. Castilho CV
    66. Cerón C
    67. Chapman CA
    68. Chave J
    69. Chhang P
    70. Chutipong W
    71. Chuyong GB
    72. Cintra BBL
    73. Clark CJ
    74. Coelho de Souza F
    75. Comiskey JA
    76. Coomes DA
    77. Cornejo Valverde F
    78. Correa DF
    79. Costa FRC
    80. Costa JBP
    81. Couteron P
    82. Culmsee H
    83. Cuni-Sanchez A
    84. Dallmeier F
    85. Damasco G
    86. Dauby G
    87. Dávila N
    88. Dávila Doza HP
    89. De Alban JDT
    90. de Assis RL
    91. De Canniere C
    92. De Haulleville T
    93. de Jesus Veiga Carim M
    94. Demarchi LO
    95. Dexter KG
    96. Di Fiore A
    97. Din HHM
    98. Disney MI
    99. Djiofack BY
    100. Djuikouo MNK
    101. Do TV
    102. Doucet JL
    103. Draper FC
    104. Droissart V
    105. Duivenvoorden JF
    106. Engel J
    107. Estienne V
    108. Farfan-Rios W
    109. Fauset S
    110. Feeley KJ
    111. Feitosa YO
    112. Feldpausch TR
    113. Ferreira C
    114. Ferreira J
    115. Ferreira LV
    116. Fletcher CD
    117. Flores BM
    118. Fofanah A
    119. Foli EG
    120. Fonty É
    121. Fredriksson GM
    122. Fuentes A
    123. Galbraith D
    124. Gallardo Gonzales GP
    125. Garcia-Cabrera K
    126. García-Villacorta R
    127. Gomes VHF
    128. Gómez RZ
    129. Gonzales T
    130. Gribel R
    131. Guedes MC
    132. Guevara JE
    133. Hakeem KR
    134. Hall JS
    135. Hamer KC
    136. Hamilton AC
    137. Harris DJ
    138. Harrison RD
    139. Hart TB
    140. Hector A
    141. Henkel TW
    142. Herbohn J
    143. Hockemba MBN
    144. Hoffman B
    145. Holmgren M
    146. Honorio Coronado EN
    147. Huamantupa-Chuquimaco I
    148. Hubau W
    149. Imai N
    150. Irume MV
    151. Jansen PA
    152. Jeffery KJ
    153. Jimenez EM
    154. Jucker T
    155. Junqueira AB
    156. Kalamandeen M
    157. Kamdem NG
    158. Kartawinata K
    159. Kasongo Yakusu E
    160. Katembo JM
    161. Kearsley E
    162. Kenfack D
    163. Kessler M
    164. Khaing TT
    165. Killeen TJ
    166. Kitayama K
    167. Klitgaard B
    168. Labrière N
    169. Laumonier Y
    170. Laurance SGW
    171. Laurance WF
    172. Laurent F
    173. Le TC
    174. Le TT
    175. Leal ME
    176. Leão de Moraes Novo EM
    177. Levesley A
    178. Libalah MB
    179. Licona JC
    180. de A Lima Filho D
    181. Lindsell JA
    182. Lopes A
    183. Lopes MA
    184. Lovett JC
    185. Lowe R
    186. Lozada JR
    187. Lu X
    188. Luambua NK
    189. Luize BG
    190. Maas P
    191. Magalhães JLL
    192. Magnusson WE
    193. Mahayani NPD
    194. Makana JR
    195. Malhi Y
    196. Maniguaje Rincón L
    197. Mansor A
    198. Manzatto AG
    199. Marimon BS
    200. Marimon-Junior BH
    201. Marshall AR
    202. Martins MP
    203. Mbayu FM
    204. de Medeiros MB
    205. Mesones I
    206. Metali F
    207. Mihindou V
    208. Millet J
    209. Milliken W
    210. Mogollón HF
    211. Molino JF
    212. Mohd Said MN
    213. Monteagudo Mendoza A
    214. Montero JC
    215. Moore S
    216. Mostacedo B
    217. Mozombite Pinto LF
    218. Mukul SA
    219. Munishi PKT
    220. Nagamasu H
    221. Nascimento HEM
    222. Nascimento MT
    223. Neill D
    224. Nilus R
    225. Noronha JC
    226. Nsenga L
    227. Núñez Vargas P
    228. Ojo L
    229. Oliveira AA
    230. de Oliveira EA
    231. Ondo FE
    232. Palacios Cuenca W
    233. Pansini S
    234. Pansonato MP
    235. Paredes MR
    236. Paudel E
    237. Pauletto D
    238. Pearson RG
    239. Pena JLM
    240. Pennington RT
    241. Peres CA
    242. Permana A
    243. Petronelli P
    244. Peñuela Mora MC
    245. Phillips JF
    246. Phillips OL
    247. Pickavance G
    248. Piedade MTF
    249. Pitman NCA
    250. Ploton P
    251. Popelier A
    252. Poulsen JR
    253. Prieto A
    254. Primack RB
    255. Priyadi H
    256. Qie L
    257. Quaresma AC
    258. de Queiroz HL
    259. Ramirez-Angulo H
    260. Ramos JF
    261. Reis NFC
    262. Reitsma J
    263. Revilla JDC
    264. Riutta T
    265. Rivas-Torres G
    266. Robiansyah I
    267. Rocha M
    268. Rodrigues D
    269. Rodriguez-Ronderos ME
    270. Rovero F
    271. Rozak AH
    272. Rudas A
    273. Rutishauser E
    274. Sabatier D
    275. Sagang LB
    276. Sampaio AF
    277. Samsoedin I
    278. Satdichanh M
    279. Schietti J
    280. Schöngart J
    281. Scudeller VV
    282. Seuaturien N
    283. Sheil D
    284. Sierra R
    285. Silman MR
    286. Silva TSF
    287. da Silva Guimarães JR
    288. Simo-Droissart M
    289. Simon MF
    290. Sist P
    291. Sousa TR
    292. de Sousa Farias E
    293. de Souza Coelho L
    294. Spracklen DV
    295. Stas SM
    296. Steinmetz R
    297. Stevenson PR
    298. Stropp J
    299. Sukri RS
    300. Sunderland TCH
    301. Suzuki E
    302. Swaine MD
    303. Tang J
    304. Taplin J
    305. Taylor DM
    306. Tello JS
    307. Terborgh J
    308. Texier N
    309. Theilade I
    310. Thomas DW
    311. Thomas R
    312. Thomas SC
    313. Tirado M
    314. Toirambe B
    315. de Toledo JJ
    316. Tomlinson KW
    317. Torres-Lezama A
    318. Tran HD
    319. Tshibamba Mukendi J
    320. Tumaneng RD
    321. Umaña MN
    322. Umunay PM
    323. Urrego Giraldo LE
    324. Valderrama Sandoval EH
    325. Valenzuela Gamarra L
    326. Van Andel TR
    327. van de Bult M
    328. van de Pol J
    329. van der Heijden G
    330. Vasquez R
    331. Vela CIA
    332. Venticinque EM
    333. Verbeeck H
    334. Veridiano RKA
    335. Vicentini A
    336. Vieira ICG
    337. Vilanova Torre E
    338. Villarroel D
    339. Villa Zegarra BE
    340. Vleminckx J
    341. von Hildebrand P
    342. Vos VA
    343. Vriesendorp C
    344. Webb EL
    345. White LJT
    346. Wich S
    347. Wittmann F
    348. Zagt R
    349. Zang R
    350. Zartman CE
    351. Zemagho L
    352. Zent EL
    353. Zent S
    (2024) Consistent patterns of common species across tropical tree communities
    Nature 625:728–734.
    https://doi.org/10.1038/s41586-023-06820-z
  2. Book
    1. Crow JF
    2. Kimura M
    (1970)
    An Introduction to Population Genetics Theory
    Harper & Row, Publishers.
  3. Book
    1. Gause GF
    (2003)
    The Struggle for Existence
    Courier Corporation.
  4. Book
    1. Hubbell SP
    (2001)
    The Unified Neutral Theory of Biodiversity and Biogeography, Monographs in Population Biology
    Princeton University Press.
    1. ter Steege H
    2. Pitman NCA
    3. Sabatier D
    4. Baraloto C
    5. Salomão RP
    6. Guevara JE
    7. Phillips OL
    8. Castilho CV
    9. Magnusson WE
    10. Molino J-F
    11. Monteagudo A
    12. Núñez Vargas P
    13. Montero JC
    14. Feldpausch TR
    15. Coronado ENH
    16. Killeen TJ
    17. Mostacedo B
    18. Vasquez R
    19. Assis RL
    20. Terborgh J
    21. Wittmann F
    22. Andrade A
    23. Laurance WF
    24. Laurance SGW
    25. Marimon BS
    26. Marimon B-H Jr
    27. Guimarães Vieira IC
    28. Amaral IL
    29. Brienen R
    30. Castellanos H
    31. Cárdenas López D
    32. Duivenvoorden JF
    33. Mogollón HF
    34. de A Matos FD
    35. Dávila N
    36. García-Villacorta R
    37. Stevenson Diaz PR
    38. Costa F
    39. Emilio T
    40. Levis C
    41. Schietti J
    42. Souza P
    43. Alonso A
    44. Dallmeier F
    45. Montoya AJD
    46. Fernandez Piedade MT
    47. Araujo-Murakami A
    48. Arroyo L
    49. Gribel R
    50. Fine PVA
    51. Peres CA
    52. Toledo M
    53. Aymard C GA
    54. Baker TR
    55. Cerón C
    56. Engel J
    57. Henkel TW
    58. Maas P
    59. Petronelli P
    60. Stropp J
    61. Zartman CE
    62. Daly D
    63. Neill D
    64. Silveira M
    65. Paredes MR
    66. Chave J
    67. de A Lima Filho D
    68. Jørgensen PM
    69. Fuentes A
    70. Schöngart J
    71. Cornejo Valverde F
    72. Di Fiore A
    73. Jimenez EM
    74. Peñuela Mora MC
    75. Phillips JF
    76. Rivas G
    77. van Andel TR
    78. von Hildebrand P
    79. Hoffman B
    80. Zent EL
    81. Malhi Y
    82. Prieto A
    83. Rudas A
    84. Ruschell AR
    85. Silva N
    86. Vos V
    87. Zent S
    88. Oliveira AA
    89. Schutz AC
    90. Gonzales T
    91. Trindade Nascimento M
    92. Ramirez-Angulo H
    93. Sierra R
    94. Tirado M
    95. Umaña Medina MN
    96. van der Heijden G
    97. Vela CIA
    98. Vilanova Torre E
    99. Vriesendorp C
    100. Wang O
    101. Young KR
    102. Baider C
    103. Balslev H
    104. Ferreira C
    105. Mesones I
    106. Torres-Lezama A
    107. Urrego Giraldo LE
    108. Zagt R
    109. Alexiades MN
    110. Hernandez L
    111. Huamantupa-Chuquimaco I
    112. Milliken W
    113. Palacios Cuenca W
    114. Pauletto D
    115. Valderrama Sandoval E
    116. Valenzuela Gamarra L
    117. Dexter KG
    118. Feeley K
    119. Lopez-Gonzalez G
    120. Silman MR
    (2013) Hyperdominance in the Amazonian tree flora
    Science 342:1243092.
    https://doi.org/10.1126/science.1243092

Article and author information

Author details

  1. David A Kessler

    Department of Physics, Bar-Ilan University, Ramat-Gan, Ramat-Gan, Israel
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    David.Kessler@biu.ac.il
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5279-1655
  2. Nadav M Shnerb

    Department of Physics, Bar-Ilan University, Ramat-Gan, Ramat-Gan, Israel
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    nadav.shnerb@gmail.com
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4418-6284

Funding

Israeli Ministry of Innovation, Science and Technology (Israel-Italy Cooperation Program)

  • Nadav M Shnerb

Israel Science Foundation (1614/21)

  • Nadav M Shnerb

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Emil Mallmin and Egbert van Nes for the interesting exchange of ideas that contributed to the shaping of this work. NMS acknowledges funding from the Israeli Ministry of Innovation, Science and Technology (MOST), under the framework of the Israel–Italy Cooperation Program. DAK acknowledges funding from the Israel Science Foundation, Grant 1614/21.

Copyright

© 2025, Kessler and Shnerb

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 765
    views
  • 34
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. David A Kessler
  2. Nadav M Shnerb
(2025)
The dominant–egalitarian transition in species-rich communities
eLife 14:e103999.
https://doi.org/10.7554/eLife.103999

Share this article

https://doi.org/10.7554/eLife.103999