Dispersal limitation and environmental heterogeneity shape scale-dependent diversity patterns in plant communities

Understanding the large-scale distribution of species diversity requires distinguishing two of the primary factors that cause compositional differences: dispersal limitation and environmental variation. In a community with a naturally discontinuous spatial structure, we asked (1) at what scale(s) nonrandom variation in species composition occurs and (2) at what scale(s) such variation is associated with spatial separation, indicative of dispersal limitation, and at what scale(s) variation is associated with environmental heterogeneity? We sampled 50 seeps (small wetlands) on five serpentine outcrops. Using a randomization model, we showed that additive beta diversity (a measure of community dissimilarity) was lower than random within seeps and higher than random among both seeps and outcrops. Using Mantel tests, we showed that plant community dissimilarity, in both the full seep assemblage as well as in a subset of seep endemics, at the two larger scales was associated with different forms of environmental heterogeneity and, at the largest scale, was also associated with geographic distance. We conclude that diversity in this system is shaped by multiple scales of heterogeneity and by dispersal limitation at the largest scale.     

Microsite-limited recruitment controls fern colonization of post-agricultural forests

Assessing the relative roles of dispersal limitation and environmental effects in population dynamics and community assembly is fundamental to understanding patterns of species distribution and diversity. In forests growing on abandoned agricultural lands, both legacies of vegetation disturbance and changes in the abiotic environment shape the diversity and composition of recovering communities. Here I specify how interactions among historical, environmental, and biological factors influence species distributions, focusing on three fern species with contrasting distributions across forests of different history in central New York, USA: Dryopteris carthusiana, Dryopteris intermedia, and Polystichum acrostichoides. Using population surveys, spore-trap and spore-bank studies, and a three-year field experiment, I compare demographic rates among species and between forest types to determine which life history stages limit colonization and which traits explain species distributions. Adult plants of all three species were larger and more likely to produce spores in post-agricultural forests than in adjacent, uncleared stands. Though lower population densities led to fewer spores in post-agricultural soils, spore availability still exceeded recruitment by four to five orders of magnitude. Sowing additional spores had relatively little effect, while microhabitat conditions had the greatest impact on establishment rates. Given similar microsites, the two forest types had equal rates of establishment, but some forest-floor features preferentially occupied by juvenile plants were less frequent in post-agricultural stands. The availability of suitable sites for establishment, created by small-scale heterogeneity on forest floors, thus limits both the growth of fern populations and the colonization of new habitats. In fact, reduced microtopographic variation in post-agricultural forests may represent a greater hindrance to plant establishment than changes in mean environmental conditions. Among the three fern species, establishment rates differed as species distributions would predict, with the strongest colonizer consistently having the highest rates and the slowest colonizer the lowest. Rather than random or trait-mediated dispersal, the different distributions of these species reflect life history traits that determine establishment rates and thus colonization ability. This case study demonstrates that ecological interactions based on the unique life histories of individual species can override dispersal in determining species distributions.     

Trait similarity patterns within grass and grasshopper communities: multitrophic community assembly at work

Trait-based community assembly theory suggests that trait variation among co-occurring species is shaped by two main processes: abiotic filtering, important in stressful environments and promoting similarity, and competition, more important in productive environments and promoting dissimilarity. Previous studies have indeed found trait similarity to decline along productivity gradients. However, these studies have always been done on single trophic levels. Here, we investigated how interactions between trophic levels affect trait similarity patterns along environmental gradients. We propose three hypotheses for the main drivers of trait similarity patterns of plants and herbivores along environmental gradients: (1) environmental control of both, (2) bottom-up control of herbivore trait variation, and (3) top-down control of grass trait variation. To test this, we collected data on the community composition and trait variation of grasses (41 species) and grasshoppers (53 species) in 50 plots in a South African savanna. Structural equation models were used to investigate how the range and spacing of within-community functional trait values of both grasses and their insect herbivores (grasshoppers; Acrididae) respond to (1) rainfall and fire frequency gradients and (2) the trait similarity patterns of the other trophic level. The analyses revealed that traits of co-occurring grasses became more similar toward lower rainfall and higher fire frequency (environmental control), while showing little evidence for top-down control. Grasshopper trait range patterns, on the other hand, were mostly directly driven by vegetation structure and grass trait range patterns (bottom-up control), while environmental factors had mostly indirect effects via plant traits. Our study shows the potential to expand trait-based community assembly theory to include trophic interactions.     

Competition among eucalyptus trees depends on genetic variation and resource supply

Genetic variation and environmental heterogeneity fundamentally shape the interactions between plants of the same species. According to the resource partitioning hypothesis, competition between neighbors intensifies as their similarity increases. Such competition may change in response to increasing supplies of limiting resources. We tested the resource partitioning hypothesis in stands of genetically identical (clone-origin) and genetically diverse (seed-origin) Eucalyptus trees with different water and nutrient supplies, using individual-based tree growth models. We found that genetic variation greatly reduced competitive interactions between neighboring trees, supporting the resource partitioning hypothesis. The importance of genetic variation for Eucalyptus growth patterns depended strongly on local stand structure and focal tree size. This suggests that spatial and temporal variation in the strength of species interactions leads to reversals in the growth rank of seed-origin and clone-origin trees. This study is one of the first to experimentally test the resource partitioning hypothesis for intergenotypic vs. intragenotypic interactions in trees. We provide evidence that variation at the level of genes, and not just species, is functionally important for driving individual and community-level processes in forested ecosystems.     

Effects of howler monkey reintroduction on ecological interactions and processes

Rewilding has been an increasingly popular tool to restore plant–animal interactions and ecological processes impaired by defaunation. However, the reestablishment of such processes has seldom been assessed. We investigated the restoration of ecological interactions following the reintroduction of the brown howler monkey (Alouatta guariba) to a defaunated Atlantic forest site. We expected the reintroduction to restore plant–animal interactions and interactions between howlers and dung beetles, which promote secondary seed dispersal. We estimated the number of interactions expected to be restored by the reintroduction to provide the baseline interaction richness that could be restored. We followed the reintroduced howler monkeys twice a week for 24 months (337 hours total) to assess their diet. We used howler monkey dung in secondary seed dispersal experiments with 2484 seed mimics to estimate the removal rates by dung beetles and collected the beetles to assess community attributes. We compared the potential future contribution of howler monkeys and other frugivores to seed dispersal based on the seed sizes they disperse in other areas where they occur. In 2 years, howler monkeys consumed 60 animal-dispersed plant species out of the 330 estimated. Twenty-one dung beetle species were attracted to experimentally provided dung; most of them were tunnelers, nocturnal, and large-sized (>10 mm). On average 30% (range 0–100%) of the large seed mimics (14 mm) were moved by dung beetles. About 91% of the species consumed by howlers (size range 0.3–34.3 mm) overlapped in seed size with those removed by dung beetles. In our study area, howler monkeys may consume more large-seeded fruit species than most other frugivores, highlighting their potential to affect forest regeneration. Our results show reintroductions may effectively restore ecological links and enhance ecological processes.     

Evolving metacommunities: toward an evolutionary perspective on metacommunities

The metacommunity framework predicts that local coexistence depends on the outcome of local species interactions and regional migration. In analogous fashion, spatial structure among populations can shape species interactions through evolutionary mechanisms. Yet, most metacommunity theories assume that populations do not evolve. Here, we evaluate how evolution shapes local species coexistence and exclusion within the multiscale and multispecies context embodied by the metacommunity framework. In general, coexistence in joint ecological-evolutionary models requires low to intermediate dispersal rates that can promote maintenance of both regional species and genetic diversity. These conditions support a set of key mechanisms that modify patterns of species coexistence including local adaptation, gene storage effects, genetic rescue effects, spatial genetic subsidies, and metacommunity evolution. Multispecies extensions indicate that correlated selection can further alter the outcome of interspecific interactions depending on the magnitude and direction of correlations and shape of fitness trade-offs. We suggest that an evolving metacommunity perspective has the potential to generate novel predictions about community structure and function by incorporating the genetic and species diversity that characterize natural communities. In adopting such a perspective, we seek to facilitate understanding about the interactions between evolutionary and metacommunity dynamics.     

Direct male-male competition can facilitate invasion of new colour types in Lake Victoria cichlids

The possibility that disruptive sexual selection alone can cause sympatric speciation is currently a subject of much debate. The initial difficulty for new and rare ornament phenotypes to invade a population, and the stabilisation of the resulting polymorphism in trait and preference make this hypothesis problematic. Recent theoretical work indicates that the invasion is facilitated if males with the new phenotype have an initial advantage in male-male competition. We studied a pair of sympatric incipient species of cichlids from Lake Victoria, in which the red (Pundamilia nyererei) and blue males (P. pundamilia) vigorously defend territories. Other studies suggested that red phenotypes may have repeatedly invaded blue populations in independent episodes of speciation. We hypothesised that red coloration confers an advantage in male-male competition, assisting red phenotypes to invade. To test this hypothesis, we staged contests between red and blue males from a population where the phenotypes are interbreeding morphs or incipient species. We staged contests under both white and green light condition. Green light effectively masks the difference between red and blue coloration. Red males dominated blue males under white light, but their competitive advantage was significantly diminished under green light. Contests were shorter when colour differences were visible. Experience of blue males with red males did not affect the outcome of a contest. The advantage of red over blue in combats may assist the red phenotype to invade blue populations. The apparently stable co-existence of red and blue incipient species in many populations of Lake Victoria cichlids is discussed.Communicated by C. St. Mary     

Determining the drivers of population structure in a highly urbanized landscape to inform conservation planning

Understanding the environmental contributors to population structure is of paramount importance for conservation in urbanized environments. We used spatially explicit models to determine genetic population structure under current and future environmental conditions across a highly fragmented, human‐dominated environment in Southern California to assess the effects of natural ecological variation and urbanization. We focused on 7 common species with diverse habitat requirements, home‐range sizes, and dispersal abilities. We quantified the relative roles of potential barriers, including natural environmental characteristics and an anthropogenic barrier created by a major highway, in shaping genetic variation. The ability to predict genetic variation in our models differed among species: 11–81% of intraspecific genetic variation was explained by environmental variables. Although an anthropogenically induced barrier (a major highway) severely restricted gene flow and movement at broad scales for some species, genetic variation seemed to be primarily driven by natural environmental heterogeneity at a local level. Our results show how assessing environmentally associated variation for multiple species under current and future climate conditions can help identify priority regions for maximizing population persistence under environmental change in urbanized regions.     

Interactive effects of disturbance and dispersal directionality on species richness and composition in metacommunities

Dispersal among ecological communities is usually assumed to be random in direction, or to vary in distance or frequency among species. However, a variety of natural systems and types of organisms may experience dispersal that is biased by directional currents or by gravity on hillslopes. We developed a general model for competing species in metacommunities to evaluate the role of directionally biased dispersal on species diversity, abundance, and traits. In parallel, we tested the role of directionally biased dispersal on communities in a microcosm experiment with protists and rotifers. Both the model and experiment independently demonstrated that diversity in local communities was reduced by directionally biased dispersal, especially dispersal that was biased away from disturbed patches. Abundance of species (and composition) in local communities was a product of disturbance intensity but not dispersal directionality. High disturbance selected for species with high intrinsic growth rates and low competitive abilities. Overall, our conclusions about the key role of dispersal directionality in (meta)communities seem robust and general, since they were supported both by the model, which was set in a general framework and not parameterized to fit to a specific system, and by a specific experimental test with microcosms.     

Productivity alters the scale dependence of the diversity-invasibility relationship

At small scales, areas with high native diversity are often resistant to invasion, while at large scales, areas with more native species harbor more exotic species, suggesting that different processes control the relationship between native and exotic species diversity at different spatial scales. Although the small-scale negative relationship between native and exotic diversity has a satisfactory explanation, we lack a mechanistic explanation for the change in relationship to positive at large scales. We investigated the native-exotic diversity relationship at three scales (range: 1-4000 km2) in California serpentine, a system with a wide range in the productivity of sites from harsh to lush. Native and exotic diversity were positively correlated at all three scales; it is rarer to detect a positive relationship at the small scales within which interactions between individuals occur. However, although positively correlated on average, the small-scale relationship between native and exotic diversity was positive at low-productivity sites and negative at high-productivity sites. Thus, the change in the relationship between native and exotic diversity does not depend on spatial scale per se, but occurs whenever environmental conditions change to promote species coexistence rather than competitive exclusion. This occurred within a single spatial scale when the environment shifted from being locally unproductive to productive.     

Neutral theory and the evolution of ecological equivalence

Since the publication of the unified neutral theory in 2001, there has been much discussion of the theory, pro and con. The hypothesis of ecological equivalence is the fundamental yet controversial idea behind neutral theory. Assuming trophically similar species are demographically alike (symmetric) on a per capita basis is only an approximation, but it is equivalent to asking: How many of the patterns of ecological communities are the result of species similarities, rather than of species differences? The strategy behind neutral theory is to see how far one can get with the simplification of assuming ecological equivalence before introducing more complexity. In another paper, I review the empirical evidence that led me to hypothesize ecological equivalence among many of the tree species in the species-rich tropical forest on Barro Colorado Island (BCI). In this paper, I develop a simple model for the evolution of ecological equivalence or niche convergence, using as an example evolution of the suite of life history traits characteristic of shade tolerant tropical tree species. Although the model is simple, the conclusions from it seem likely to be robust. I conclude that ecological equivalence for resource use are likely to evolve easily and often, especially in species-rich communities that are dispersal and recruitment limited. In the case of the BCI forest, tree species are strongly dispersal- and recruitment-limited, not only because of restricted seed dispersal, but also because of low recruitment success due to heavy losses of the seedling stages to predators and pathogens and other abiotic stresses such as drought. These factors and the high species richness of the community strongly reduce the potential for competitive exclusion of functionally equivalent or nearly equivalent species.     

Dispersal, competition, and shifting patterns of diversity in a degraded oak savanna

Diversity is a balance between processes that add and limit species (e.g., dispersal vs. competition), but reconciling their contributions remains a challenge. Recruit-ment-based models predict that dispersal barriers are most limiting for diversity, while competition-based models predict that dispersal matters only when competition is minimized. Testing these models is difficult because their influence varies with scale and site productivity. In a degraded oak savanna, we used plot-level (seed additions, burning) and site-level (proportions of regional functional groups found locally) analyses in areas with variable soil depth to examine how dispersal and competition influence diversity. At the plot level, added species persisted where they were formerly absent, but few established naturally despite fire-induced resource enrichment and nearby populations, revealing the importance of dispersal limitation for diversity. This result did not vary with soil depth or standing crop. Although competition could not prevent establishment in unburned plots, it significantly lowered survival, indicating that resource limitations exacerbate dispersal inefficiencies. At the site level, the concordance between regional and local diversity for native species was associated with soil depth heterogeneity, not dispersal or competition. This suggests that persistence is determined primarily by the influence of the environment on population demographics. Given that the formation of new populations is unlikely, those remaining appear to be confined to optimal habitat where they resist competitive or stochastic displacement, possibly explaining why species loss is rare despite substantial habitat loss and invasion. For exotics, there was no relationship between diversity and soil depth heterogeneity. Annuals with presumed dispersal capabilities were significantly overrepresented in all sites while perennial forbs, the largest regional functional group, were significantly underrepresented. We interpret the native-exotic discrepancies as reflecting the recent arrival of exotics (150 years ago), suggesting that local establishment occurs slowly even for species with regional prevalence. The accumulation lag may be explained by the need for founder populations to be demographically stable; otherwise persistence requires continual immigration favoring overrepresentation by dispersers. Our findings support the view that dispersal limitation restricts diversity within plant communities, but suggests that the impacts of environment on demographic performance ultimately determine the pattern and rate of community assembly.     

Mass effects mediate coexistence in competing shrews

Recent developments in metacommunity theory have raised awareness that processes occurring at regional scales might interfere with local dynamics and affect conditions for the local coexistence of competing species. Four main paradigms are recognized in this context (namely, neutral, patch-dynamics, species-sorting, and mass-effect), which differ according to the role assigned to ecological or life-history differences among competing species, as well as to the relative time scale of regional vs. local dynamics. We investigated the patterns of regional and local coexistence of two species of shrews (Crocidura russula and Sorex coronatus) sharing a similar diet (generalist insectivores) over four generations, in a spatially structured habitat at the altitudinal limit of their distributions. Local populations were small, and regional dynamics were strong, with high rates of extinction and recolonization. Niche analysis revealed significant habitat differentiation on a few important variables, including temperature and availability of winter resting sites. In sites suitable for both species, we found instances of local coexistence with no evidence of competitive exclusion. Patterns of temporal succession did not differ from random, with no suggestion of a colonization-competition trade-off. Altogether, our data provide support for the mass-effect paradigm, where regional coexistence is mediated by specialization on different habitat types, and local coexistence by rescue effects from source sites. The strong regional dynamics and demographic stochasticity, together with high dispersal rates, presumably contributed to mass effects by overriding local differences in specific competitive abilities.     

The community effects of phenotypic and genetic variation within a predator population

Natural populations are heterogeneous mixtures of individuals differing in physiology, morphology, and behavior. Despite the ubiquity of phenotypic variation within natural populations, its effects on the dynamics of ecological communities are not well understood. Here, we use a quantitative genetics framework to examine how phenotypic variation in a predator affects the outcome of apparent competition between its two prey species. Classical apparent competition theory predicts that prey have reciprocally negative effects on each other. The addition of phenotypic trait variation in predation can marginalize these negative effects, mediate coexistence, or generate positive indirect effects between the prey species. Long-term coexistence or facilitation, however, can be preceded by long transients of extinction risk whenever the heritability of phenotypic variation is low. Greater heritability can circumvent these ecological transients but also can generate oscillatory and chaotic dynamics. These dramatic changes in ecological outcomes, in the sign of indirect effects, and in stability suggest that studies which ignore intraspecific trait variation may reach fundamentally incorrect conclusions regarding ecological dynamics.     

Change in avian functional fingerprints of a Neotropical montane forest over 100 years as an indicator of ecosystem integrity

Ecologically relevant traits of organisms in an assemblage determine an ecosystem's functional fingerprint (i.e., the shape, size, and position of multidimensional trait space). Quantifying changes in functional fingerprints can therefore provide information about the effects of diversity loss or gain through time on ecosystem condition and is a promising approach to monitoring ecological integrity. This, however, is seldom possible owing to limitations in historical surveys and a lack of data on organismal traits, particularly in diverse tropical regions. Using data from detailed bird surveys from 4 periods across more than a century, and morphological and ecological traits of 233 species, we quantified changes in the avian functional fingerprint of a tropical montane forest in the Andes of Colombia. We found that 78% of the variation in functional space, regardless of period, was described by 3 major axes summarizing body size, dispersal ability (indexed by wing shape), and habitat breadth. Changes in species composition significantly altered the functional fingerprint of the assemblage and functional richness and dispersion decreased 35–60%. Owing to species extirpations and to novel additions to the assemblage, functional space decreased over time, but at least 11% of its volume in the 2010s extended to areas of functional space that were unoccupied in the 1910s. The assemblage now includes fewer large-sized species, more species with greater dispersal ability, and fewer habitat specialists. Extirpated species had high functional uniqueness and distinctiveness, resulting in large reductions in functional richness and dispersion after their loss, which implies important consequences for ecosystem integrity. Conservation efforts aimed at maintaining ecosystem function must move beyond seeking to sustain species numbers to designing complementary strategies for the maintenance of ecological function by identifying and conserving species with traits conferring high vulnerability such as large body size, poor dispersal ability, and greater habitat specialization. Article impact statement: Changes in functional fingerprints provide a means to quantify the integrity of ecological assemblages affected by diversity loss or gain.     

The reef environment and competitive success in the Corallimorpharia

Competitive success within coral reef communities is controlled by various factors. In addition to competitive abilities in direct interactions with a contestant, external influences such as disturbance caused by nutrient input may determine the outcome of antagonistic interactions. We examined the competitive success of corallimorpharians on coral reefs by investigating their distribution patterns within reefs and how well they perform in interference competition with staghorn corals in different environments. Substrate composition and corallimorpharian growth were examined on three reefs in Tanzania under different disturbance regimes using the line-intercept transect and point techniques. A transplant experiment was conducted in which staghorn corals (Acropora formosa) were exposed to the polyps of Rhodactis rhodostoma to establish how competition between corals and corallimorpharians affects their respective distributions. Within reefs corallimorpharians seemed to be more competitive in shallow waters. This could be due to both environmental factors as well as varied competitive abilities depending on surrounding benthos that changed with depth. Reef environment also seemed to influence corallimorpharian growth among reefs as they had the highest densities in the areas with the highest nutrient loads. The transplant experiment revealed that the corallimorpharians had a competitive advantage over the corals, and in comparisons of reefs influenced by different degrees of disturbance, corallimorpharians were most competitive in the area with the highest nutrient content. Hence, stress on coral reefs in the form of raised nutrient loads may favour the competitive success of corallimorpharians.Communicated by P.W. Sammarco, Chauvin     

Size-mediated non-trophic interactions and stochastic predation drive assembly and dynamics in a seabird community

Theoretical and empirical evidence suggests that body size is a major life-history trait impacting on the structure and functioning of complex food webs. However, long-term analyses of size-dependent interactions within simpler network modules, for instance, competitive guilds, are scant. Here, we model the assembly dynamics of the largest breeding seabird community in the Mediterranean basin during the last 30 years. This unique data set allowed us to test, through a "natural experiment," whether body size drove the assembly and dynamics of an ecological guild growing from very low numbers after habitat protection. Although environmental stochasticity accounted for most of community variability, the population variance explained by interspecific interactions, albeit small, decreased sharply with increasing body size. Since we found a demographic gradient along a body size continuum, in which population density and stability increase with increasing body size, the numerical effects of interspecific interactions were proportionally higher on smaller species than on larger ones. Moreover, we found that the per capita interaction coefficients were larger the higher the size ratio among competing species, but only for the set of interactions in which the species exerting the effect was greater. This provides empirical evidence for long-term asymmetric interspecific competition, which ultimately prompted the local extinction of two small species during the study period. During the assembly process stochastic predation by generalist carnivores further triggered community reorganizations and global decays in population synchrony, which disrupted the pattern of interspecific interactions. These results suggest that the major patterns detected in complex food webs can hold as well for simpler sub-modules of these networks involving non-trophic interactions, and highlight the shifting ecological processes impacting on assembling vs. asymptotic communities.     

Quantifying the complexities of Saccharomyces cerevisiae's ecosystem engineering via fermentation

The theory of niche construction suggests that organisms may engineer environments via their activities. Despite the potential of this phenomenon being realized by Darwin, the capability of niche construction to generally unite ecological and evolutionary biology has never been empirically quantified. Here I quantify the fitness effects of Saccharomyces cerevisiae's ecosystem engineering in a natural ferment in order to understand the interaction between ecological and evolutionary processes. I show that S. cerevisiae eventually dominates in fruit niches, where it is naturally initially rare, by modifying the environment through fermentation (the Crabtree effect) in ways which extend beyond just considering ethanol production. These data show that an additional cause of S. cerevisiae's competitive advantage over the other yeasts in the community is due to the production of heat via fermentation. Even though fermentation is less energetically efficient than respiration, it seems that this trait has been selected for because its net effect provides roughly a 7% fitness advantage over the other members of the community. These data provide an elegant example of niche construction because this trait clearly modifies the environment and therefore the selection pressures to which S. cerevisiae, and other organisms that access the fruit resource, including humans, are exposed to.     

Measuring impacts on species with models and metrics of varying ecological and computational complexity

Approaches to assess the impacts of landscape disturbance scenarios on species range from metrics based on patterns of occurrence or habitat to comprehensive models that explicitly include ecological processes. The choice of metrics and models affects how impacts are interpreted and conservation decisions. We explored the impacts of 3 realistic disturbance scenarios on 4 species with different ecological and taxonomic traits. We used progressively more complex models and metrics to evaluate relative impact and rank of scenarios on the species. Models ranged from species distribution models that relied on implicit assumptions about environmental factors and species presence to highly parameterized spatially explicit population models that explicitly included ecological processes and stochasticity. Metrics performed consistently in ranking different scenarios in order of severity primarily when variation in impact was driven by habitat amount. However, they differed in rank for cases where dispersal dynamics were critical in influencing metapopulation persistence. Impacts of scenarios on species with low dispersal ability were better characterized using models that explicitly captured these processes. Metapopulation capacity provided rank orders that most consistently correlated with those from highly parameterized and data-rich models and incorporated information about dispersal with little additional computational and data cost. Our results highlight the importance of explicitly considering species’ ecology, spatial configuration of habitat, and disturbance when choosing indicators of species persistence. We suggest using hybrid approaches that are a mixture of simple and complex models to improve multispecies assessments.     

The effect of positive interactions on community structure in a multi-species metacommunity model along an environmental gradient

Positive interactions are widely recognized as playing a major role in the organization of community structure and diversity. As such, recent theoretical and empirical works have revealed the significant contribution of positive interactions in shaping species’ geographical distributions, particularly in harsh abiotic conditions. In this report, we explore the joint influence of local dispersal and an environmental gradient on the spatial distribution, structure and function of communities containing positive interactions. While most previous theoretical efforts were limited to modelling the dynamics of single pairs of associated species being mutualist or competitor, here we employ a spatially explicit multi-species metacommunity model covering a rich range of interspecific interactions (mutualism, competition and exploitation) along an environmental gradient. We find that mutualistic interactions dominate in communities with low diversity characterized by limited species dispersal and poor habitat quality. On the other hand, the fraction of mutualistic interactions decreases at the expense of exploitation and competition with the increase in diversity caused by higher dispersal and/or habitat quality. Our multi-species model exemplifies the ubiquitous presence of mutualistic interactions and the role of mutualistic species as facilitators for the further establishment of species during ecosystem assembly. We therefore argue that mutualism is an essential component driving the origination of complex and diverse communities.