The wealth of species: ecological communities, complex systems and the legacy of Frank Preston

Dissimilarity in flea and host assemblages and their interaction networks along a spatial distance gradient: different patterns revealed by different network dissimilarity metrics

We investigated the distance-decay pattern (an increase in dissimilarity with increasing geographic distance) in regional assemblages of fleas and their small mammalian hosts, as well as their interaction networks, in four biogeographic realms. Dissimilarity of assemblages (βtotal) was partitioned into species richness differences (βrich) and species replacement (βrepl) components. Dissimilarity of networks was assessed using two metrics: (a) whole network dissimilarity (βWN) partitioned into species replacement (βST) and interaction rewiring (βOS) components and (b) D statistics, measuring dissimilarity in the pure structure of the networks, without using information on species identities and calculated for hosts-shared-by-fleas networks (Dh) and fleas-shared-by-hosts networks (Df). We asked whether the distance-decay pattern (a) occurs among interactor assemblages or their interaction networks; (b) depends on the network dissimilarity metric used; and (c) differs between realms. The βtotal and βrepl of flea and host assemblages increased with distance in all realms except for host assemblages in the Afrotropics. βrich for flea and host assemblages increased with distance in the Nearctic only. In networks, βWN and βST demonstrated a distance-decay pattern, whereas βOS was mainly spatially invariant except in the Neotropics. Correlations of Dh or Df and geographic distance were mostly non-significant. We conclude that investigations of dissimilarity in interaction networks should include both types of dissimilarity metrics (those that consider partner identities and those that consider the pure structure of networks). This will allow elucidating the predictability of some facets of network dissimilarity and the unpredictability of other facets.

Consistent patterns of common species across tropical tree communities

Trees structure the Earth's most biodiverse ecosystem, tropical forests. The vast number of tree species presents a formidable challenge to understanding these forests, including their response to environmental change, as very little is known about most tropical tree species. A focus on the common species may circumvent this challenge. Here we investigate abundance patterns of common tree species using inventory data on 1,003,805 trees with trunk diameters of at least 10 cm across 1,568 locations1-6 in closed-canopy, structurally intact old-growth tropical forests in Africa, Amazonia and Southeast Asia. We estimate that 2.2%, 2.2% and 2.3% of species comprise 50% of the tropical trees in these regions, respectively. Extrapolating across all closed-canopy tropical forests, we estimate that just 1,053 species comprise half of Earth's 800 billion tropical trees with trunk diameters of at least 10 cm. Despite differing biogeographic, climatic and anthropogenic histories7, we find notably consistent patterns of common species and species abundance distributions across the continents. This suggests that fundamental mechanisms of tree community assembly may apply to all tropical forests. Resampling analyses show that the most common species are likely to belong to a manageable list of known species, enabling targeted efforts to understand their ecology. Although they do not detract from the importance of rare species, our results open new opportunities to understand the world's most diverse forests, including modelling their response to environmental change, by focusing on the common species that constitute the majority of their trees.

Analysing the distance decay of community similarity in river networks using Bayesian methods

The distance decay of community similarity (DDCS) is a pattern that is widely observed in terrestrial and aquatic environments. Niche-based theories argue that species are sorted in space according to their ability to adapt to new environmental conditions. The ecological neutral theory argues that community similarity decays due to ecological drift. The continuum hypothesis provides an intermediate perspective between niche-based theories and the neutral theory, arguing that niche and neutral factors are at the opposite ends of a continuum that ranges from competitive to stochastic exclusion. We assessed the association between niche-based and neutral factors and changes in community similarity measured by Sorensen’s index in riparian plant communities. We assessed the importance of neutral processes using network distances and flow connection and of niche-based processes using Strahler order differences and precipitation differences. We used a hierarchical Bayesian approach to determine which perspective is best supported by the results. We used dataset composed of 338 vegetation censuses from eleven river basins in continental Portugal. We observed that changes in Sorensen indices were associated with network distance, flow connection, Strahler order difference and precipitation difference but to different degrees. The results suggest that community similarity changes are associated with environmental and neutral factors, supporting the continuum hypothesis.

Empirical abundance distributions are more uneven than expected given their statistical baseline

Exploring and accounting for the emergent properties of ecosystems as complex systems is a promising horizon in the search for general processes to explain common ecological patterns. For example the ubiquitous hollow-curve form of the species abundance distribution is frequently assumed to reflect ecological processes structuring communities, but can also emerge as a statistical phenomenon from the mathematical definition of an abundance distribution. Although the hollow curve may be a statistical artefact, ecological processes may induce subtle deviations between empirical species abundance distributions and their statistically most probable forms. These deviations may reflect biological processes operating on top of mathematical constraints and provide new avenues for advancing ecological theory. Examining ~22,000 communities, we found that empirical SADs are highly uneven and dominated by rare species compared to their statistical baselines. Efforts to detect deviations may be less informative in small communities-those with few species or individuals-because these communities have poorly resolved statistical baselines. The uneven nature of many empirical SADs demonstrates a path forward for leveraging complexity to understand ecological processes governing the distribution of abundance, while the issues posed by small communities illustrate the limitations of using this approach to study ecological patterns in small samples.

Niche dissociated assembly drives insular lizard community organization

Interspecific competition for resources leading to niche partitioning is considered as one of the major drivers of community assembly. Competitive niche partitioning is diagnosed from species co-occurrence, species abundance distributions (SADs), and body size distributions of species. For several decades, studies have explored these patterns for the relative significance of interspecific competition in shaping communities. We explored these patterns in a finite assemblage of insectivorous lizards in the Andaman & Nicobar Islands, both at the level of archipelago and individual islands. Negative geographic co-occurrences occurred only between species pairs in islands separated by deep ocean channels. Ecologically similar species did not show positive co-occurrence in guild co-occurrence analyses, indicating that the negative geographical co-occurrences between species in islands were due to historical allopatry. Species abundance distribution was best explained by a Pareto distribution in both metacommunity and local communities. There was no predictable spacing of body sizes among co-existing species in local communities. The empirical data on insular lizard community on species co-occurrence, SADs, and body size ratios does not lend support to assortment of species in islands caused by niche subdivision. Such niche-dissociated assembly of species in islands might be an important factor in formation of biological communities, regardless of geographic scale.

Emergent niche structuring leads to increased differences from neutrality in species abundance distributions

Species abundance distributions must reflect the dynamic processes involved in community assembly, but whether and when specific processes lead to distinguishable signals is not well understood. Biodiversity and species abundances may be shaped by a variety of influences, but particular attention has been paid to competition, which can involve neutral dynamics, where competitor abundances are governed only by demographic stochasticity and immigration, and dynamics driven by trait differences that enable stable coexistence through the formation of niches. Key recent studies of the species abundance patterns of communities with niches employ simple models with pre-imposed niche structure. These studies suggest that species abundance distributions are insensitive to the relative contributions of niche and neutral processes, especially when diversity is much higher than the number of niches. Here we analyze results from a stochastic population model with competition driven by trait differences. With this model, niche structure emerges as clumps of species that persist along the trait axis, and leads to more substantial differences from neutral species abundance distributions than have been previously shown. We show that heterogeneity in "between-niche" interaction strength (i.e., in the strength of competition between species in different niches) plays the dominant role in shaping the species abundances along the trait axis, acting as a biotic filter favoring species at the centers of niches. Furthermore, we show that heterogeneity in "within-niche" interactions (i.e., in the competition between species in the same niche) counteracts the influence of heterogeneity in "between-niche" interactions on the SAD to some degree. Our results suggest that competitive interactions that produce niches can also influence the shapes of SADs.

Academic inequality through the lens of community ecology: a meta-analysis

Ecological assemblages are generally characterized by a few dominant species and numerous others. Such unequal distributions of dominance also emerge in human society, including in scientific communities. Here, based on formal community ecological analyses, we show the temporal trends in the number of scientific publication in the discipline of "ecology." Based on this, we infer possible factors causing the imbalance of reputation and dominance among countries. We relied on 454 ecological meta-analysis papers published from 1998 to 2014, which sourced over 29,000 original publications. Formal meta-analyses are essential for synthesizing findings from individual studies and are critical for assessing issues and informing policy. We found that, despite the rapid expansion of outlets for ecology papers (analogous to an increase in carrying capacity, in ecological systems), country diversity as determined from first author affiliations (analogous to species diversity) did not increase. Furthermore, a country identity was more powerful than the popularity of the scientific topic and affected the chance of publication in high-profile journals, independent of the potential novelty of findings and arguments of the papers, suggesting possible academic injustice. Consequently, a rank order and hierarchy has been gradually formed among countries. Notably, this country-dominance rank is not only specific to this scientific domain but also universal across different societal situations including sports and economics, further emphasizing that inequality and hierarchical structure exist even in modern human society. Our study demonstrates a need for having robust frameworks to facilitate equality and diversity in the scientific domain in order to better inform society and policy.

Separating macroecological pattern and process: comparing ecological, economic, and geological systems

Theories of biodiversity rest on several macroecological patterns describing the relationship between species abundance and diversity. A central problem is that all theories make similar predictions for these patterns despite disparate assumptions. A troubling implication is that these patterns may not reflect anything unique about organizational principles of biology or the functioning of ecological systems. To test this, we analyze five datasets from ecological, economic, and geological systems that describe the distribution of objects across categories in the United States. At the level of functional form ('first-order effects'), these patterns are not unique to ecological systems, indicating they may reveal little about biological process. However, we show that mechanism can be better revealed in the scale-dependency of first-order patterns ('second-order effects'). These results provide a roadmap for biodiversity theory to move beyond traditional patterns, and also suggest ways in which macroecological theory can constrain the dynamics of economic systems.

The control of rank-abundance distributions by a competitive despotic species

Accounting for differences in abundances among species remains a high priority for community ecology. While there has been more than 80 years of work on trying to explain the characteristic S shape of rank-abundance distributions (RADs), there has been recent conjecture that the form may not depend on ecological processes per se but may be a general phenomenon arising in many unrelated disciplines. We show that the RAD shape can be influenced by an ecological process, namely, interference competition. The noisy miner (Manorina melanocephala) is a hyperaggressive, 'despotic' bird that occurs over much of eastern Australia (>10(6) km(2)). We compiled data for bird communities from 350 locations within its range, which were collected using standard avian survey methods. We used hierarchical Bayesian models to show that the RAD shape was much altered when the abundance of the strong interactor exceeded a threshold density; RADs consistently were steeper when the density of the noisy miner ≥2.5 birds ha(-1). The structure of bird communities at sites where the noisy miner exceeded this density was very different from that at sites where the densities fell below the threshold: species richness and Shannon diversity were much reduced, but mean abundances and mean avian biomass per site did not differ substantially.

How species richness and total abundance constrain the distribution of abundance

The species abundance distribution (SAD) is one of the most intensively studied distributions in ecology and its hollow-curve shape is one of ecology's most general patterns. We examine the SAD in the context of all possible forms having the same richness (S) and total abundance (N), i.e. the feasible set. We find that feasible sets are dominated by similarly shaped hollow curves, most of which are highly correlated with empirical SADs (most R(2) values > 75%), revealing a strong influence of N and S on the form of the SAD and an a priori explanation for the ubiquitous hollow curve. Empirical SADs are often more hollow and less variable than the majority of the feasible set, revealing exceptional unevenness and relatively low natural variability among ecological communities. We discuss the importance of the feasible set in understanding how general constraints determine observable variation and influence the forms of predicted and empirical patterns.

Is there an ecological basis for species abundance distributions?

Community ecologists have attempted to explain species abundance distribution (SAD) shape for more than 80 years, but usually without relating SAD shape explicitly to ecological variables. We explored whether the scale (total assemblage abundance) and shape (assemblage evenness) of avifaunal SADs were related to ecological covariates. We used data on avifaunas, in-site habitat structure and landscape context that were assembled from previous studies; this amounted to 197 transects distributed across 16,000 km(2) of the box-ironbark forests of southeastern Australia. We used Bayesian conditional autoregressive models to link SAD scale and shape to these ecological covariates. Variation in SAD scale was relatable to some ecological covariates, especially to landscape vegetation cover and to tree height. We could not find any relationships between SAD shape and ecological covariates. SAD shape, the core component in SAD theory, may hold little information about how assemblages are governed ecologically and may result from statistical processes, which, if general, would indicate that SAD shape is not useful for distinguishing among theories of assemblage structure.

Species-area relationships are controlled by species traits

The species-area relationship (SAR) is one of the most thoroughly investigated empirical relationships in ecology. Two theories have been proposed to explain SARs: classical island biogeography theory and niche theory. Classical island biogeography theory considers the processes of persistence, extinction, and colonization, whereas niche theory focuses on species requirements, such as habitat and resource use. Recent studies have called for the unification of these two theories to better explain the underlying mechanisms that generates SARs. In this context, species traits that can be related to each theory seem promising. Here we analyzed the SARs of butterfly and moth assemblages on islands differing in size and isolation. We tested whether species traits modify the SAR and the response to isolation. In addition to the expected overall effects on the area, traits related to each of the two theories increased the model fit, from 69% up to 90%. Steeper slopes have been shown to have a particularly higher sensitivity to area, which was indicated by species with restricted range (slope  = 0.82), narrow dietary niche (slope  = 0.59), low abundance (slope  = 0.52), and low reproductive potential (slope  = 0.51). We concluded that considering species traits by analyzing SARs yields considerable potential for unifying island biogeography theory and niche theory, and that the systematic and predictable effects observed when considering traits can help to guide conservation and management actions.

Properties of similarity indices under niche-based and dispersal-based processes in communities

Little is known about the combined impact of habitat filtering and dispersal limitation on species turnover patterns. To gain new insights, we constructed a spatially explicit community model wherein we controlled dispersal distances, the strength of habitat filtering, and the grain of habitat heterogeneity to study the distance decay of several (dis)similarity indices. The impact of habitat filtering is dependent on the ratio between the grain of habitats and the mean dispersal distance. The behavior of (dis)similarity indices varies. First, incidence-based measures of species overlap are less affected by habitat filtering than are abundance-based indices. Second, species identity-based indices, derived from population genetics' F(ST), show interesting capacities to infer dispersal processes under neutrality but are also highly sensitive to habitat filtering. All indices except F(ST)-related indices under neutrality are very sensitive to overall species richness. Hence, community patterns showing contrasted diversity levels should be compared with caution. Partitioning similarity indices within and between habitats appears to be an efficient approach to assess the strength of habitat filtering, and we show that a torus-translation test is powerful for this purpose. We finally highlight the need for further analytical studies to achieve theoretical expectations of similarity decay under dispersal and niche processes.

Universal ecological patterns in college basketball communities

The rank abundance of common and rare species within ecological communities is remarkably consistent from the tropics to the tundra. This invariant patterning provides one of ecology's most enduring and unified tenets: most species rare and a few very common. Increasingly, attention is focused upon elucidating biological mechanisms that explain these species abundance distributions (SADs), but these evaluations remain controversial. We show that college basketball wins generate SADs just like those observed in ecological communities. Whereas college basketball wins are structured by competitive interactions, the result produces a SAD pattern indistinguishable from random wins. We also show that species abundance data for tropical trees exhibits a significant-digit pattern consistent with data derived from complex structuring forces. These results cast doubt upon the ability of SAD analysis to resolve ecological mechanism, and their patterning may reflect statistical artifact as much as biological processes.

Similarity in ectoparasite faunas of Palaearctic rodents as a function of host phylogenetic, geographic or environmental distances: which matters the most?

Different host species harbour parasite faunas that are anywhere from very similar to very different in species composition. A priori, the similarity in the parasite faunas of any two host species should decrease with increases in either the phylogenetic distance, the distinctness of the environments occupied or the geographical distance between these hosts. We tested these predictions using extensive data on the faunas of fleas (Insecta: Siphonaptera) and gamasid mites (Acari: Parasitiformes) parasitic on rodents across the Palaearctic. For each pair of host species, we computed the similarity in parasite faunas based on both species composition as well as the phylogenetic and/or taxonomic distinctness of parasite species. Phylogenetic distances between hosts were based on patristic distances through a rodent phylogeny, geographic distances were computed from geographic range data, and environmental dissimilarity was measured from the average climatic and vegetation scores of each host range. Using multiple regressions on distance matrices to assess the separate explanatory power of each of the three dependent variables, environmental dissimilarity between the ranges of host species emerged as the best predictor of dissimilarity between parasite faunas, especially for fleas; in the case of mites, phylogenetic distance between host species was also important. A closer look at the data indicates that the flea and mite faunas of two hosts inhabiting different environments are always different, whilst hosts living in similar environments can have either very similar or dissimilar parasite faunas. Additional tests showed that dissimilarity in flea or mite faunas between host geographic ranges was best explained by dissimilarity in vegetation, followed by dissimilarity in climatic conditions. Thus, external environmental factors may play greater roles than commonly thought in the evolution of host-parasite associations.

Species abundance distribution results from a spatial analogy of central limit theorem

The frequency distribution of species abundances [the species abundance distribution (SAD)] is considered to be a fundamental characteristic of community structure. It is almost invariably strongly right-skewed, with most species being rare. There has been much debate as to its exact properties and the processes from which it results. Here, we contend that an SAD for a study plot must be viewed as spliced from the SADs of many smaller nonoverlapping subplots covering that plot. We show that this splicing, if applied repeatedly to produce subplots of progressively larger size, leads to the observed shape of the SAD for the whole plot regardless of that of the SADs of those subplots. The widely reported shape of an SAD is thus likely to be driven by a spatial parallel of the central limit theorem, a statistically convergent process through which the SAD arises from small to large scales. Exact properties of the SAD are driven by species spatial turnover and the spatial autocorrelation of abundances, and can be predicted using this information. The theory therefore provides a direct link between SADs and the spatial correlation structure of species distributions, and thus between several fundamental descriptors of community structure. Moreover, the statistical process described may lie behind similar frequency distributions observed in many other scientific fields.

Coupled microbial and human systems: evidence for a relationship between infectious disease and gross national product

We provide evidence that maternal metabolic energy is diverted to increased birth rates in nations experiencing high infectious disease risk. The "economic stoichiometry" of such situations limits the availability and distribution of metabolic resources available for national production. Lowering disease risk, and thus the metabolic energy required for replacement human biomass production, makes energy available for national production during the demographic transition, and increases the national GDP.

Zero sum, the niche, and metacommunities: long-term dynamics of community assembly

Recent models of community assembly, structure, and dynamics have incorporated, to varying degrees, three mechanistic processes: resource limitation and interspecific competition, niche requirements of species, and exchanges between a local community and a regional species pool. Synthesizing 30 years of data from an intensively studied desert rodent community, we show that all of these processes, separately and in combination, have influenced the structural organization of this community and affected its dynamical response to both natural environmental changes and experimental perturbations. In addition, our analyses suggest that zero-sum constraints, niche differences, and metacommunity processes are inextricably linked in the ways that they affect the structure and dynamics of this system. Explicit consideration of the interaction of these processes should yield a deeper understanding of the assembly and dynamics of other ecological communities. This synthesis highlights the role that long-term data, especially when coupled with experimental manipulations, can play in assessing the fundamental processes that govern the structure and function of ecological communities.

A general framework for the distance-decay of similarity in ecological communities

Species spatial turnover, or beta-diversity, induces a decay of community similarity with geographic distance known as the distance-decay relationship. Although this relationship is central to biodiversity and biogeography, its theoretical underpinnings remain poorly understood. Here, we develop a general framework to describe how the distance-decay relationship is influenced by population aggregation and the landscape-scale species-abundance distribution. We utilize this general framework and data from three tropical forests to show that rare species have a weak influence on distance-decay curves, and that overall similarity and rates of decay are primarily influenced by species abundances and population aggregation respectively. We illustrate the utility of the framework by deriving an exact analytical expression of the distance-decay relationship when population aggregation is characterized by the Poisson Cluster Process. Our study provides a foundation for understanding the distance-decay relationship, and for predicting and testing patterns of beta-diversity under competing theories in ecology.

Species abundance distributions: moving beyond single prediction theories to integration within an ecological framework

Species abundance distributions (SADs) follow one of ecology's oldest and most universal laws--every community shows a hollow curve or hyperbolic shape on a histogram with many rare species and just a few common species. Here, we review theoretical, empirical and statistical developments in the study of SADs. Several key points emerge. (i) Literally dozens of models have been proposed to explain the hollow curve. Unfortunately, very few models are ever rejected, primarily because few theories make any predictions beyond the hollow-curve SAD itself. (ii) Interesting work has been performed both empirically and theoretically, which goes beyond the hollow-curve prediction to provide a rich variety of information about how SADs behave. These include the study of SADs along environmental gradients and theories that integrate SADs with other biodiversity patterns. Central to this body of work is an effort to move beyond treating the SAD in isolation and to integrate the SAD into its ecological context to enable making many predictions. (iii) Moving forward will entail understanding how sampling and scale affect SADs and developing statistical tools for describing and comparing SADs. We are optimistic that SADs can provide significant insights into basic and applied ecological science.

Global patterns of city size distributions and their fundamental drivers

Urban areas and their voracious appetites are increasingly dominating the flows of energy and materials around the globe. Understanding the size distribution and dynamics of urban areas is vital if we are to manage their growth and mitigate their negative impacts on global ecosystems. For over 50 years, city size distributions have been assumed to universally follow a power function, and many theories have been put forth to explain what has become known as Zipf's law (the instance where the exponent of the power function equals unity). Most previous studies, however, only include the largest cities that comprise the tail of the distribution. Here we show that national, regional and continental city size distributions, whether based on census data or inferred from cluster areas of remotely-sensed nighttime lights, are in fact lognormally distributed through the majority of cities and only approach power functions for the largest cities in the distribution tails. To explore generating processes, we use a simple model incorporating only two basic human dynamics, migration and reproduction, that nonetheless generates distributions very similar to those found empirically. Our results suggest that macroscopic patterns of human settlements may be far more constrained by fundamental ecological principles than more fine-scale socioeconomic factors.

The maximum entropy formalism and the idiosyncratic theory of biodiversity

Why does the neutral theory, which is based on unrealistic assumptions, predict diversity patterns so accurately? Answering questions like this requires a radical change in the way we tackle them. The large number of degrees of freedom of ecosystems pose a fundamental obstacle to mechanistic modelling. However, there are tools of statistical physics, such as the maximum entropy formalism (MaxEnt), that allow transcending particular models to simultaneously work with immense families of models with different rules and parameters, sharing only well-established features. We applied MaxEnt allowing species to be ecologically idiosyncratic, instead of constraining them to be equivalent as the neutral theory does. The answer we found is that neutral models are just a subset of the majority of plausible models that lead to the same patterns. Small variations in these patterns naturally lead to the main classical species abundance distributions, which are thus unified in a single framework.