Stochastic and deterministic processes drive wetland community assembly across a gradient of environmental filtering

The role of deterministic and stochastic processes in community assembly is a key question in community ecology. We evaluated the effect of an abiotic filter (hydroperiod) on the partitioned diversity of three taxonomic groups (birds, vegetation, macroinvertebrates) from prairie pothole wetlands in Alberta, Canada, which naturally vary in water permanence. We observed that alpha and gamma diversity were higher in permanent than temporary wetlands (16–25% and 34–47% respectively, depending on the taxon). This suggests an influence of deterministic constraints on the number of species a wetland can support. Taxa which cannot persist in shallow, temporary wetlands are excluded by the deterministic constraints that a shortened hydroperiod imposes. In contrast, we observed that beta diversity was significantly higher (2–12%) in temporary wetlands than permanent ones, and temporary wetlands supported more unique combinations of community composition than permanent wetlands, despite having a smaller regional species pool. This observation contradicts prior mesocosm studies that found beta diversity mirrored the pattern in gamma diversity along an environmental filtering gradient. We conclude that deterministic processes are more influential in more stable permanent wetlands, whereas stochastic processes play a more important role in assembly in dynamic temporary wetlands that must disassemble and re‐establish annually. Considering three distinct taxonomic groups differing in their relative mobility, our large‐scale field study demonstrates that both stochastic and deterministic processes act together to influence the assembly of multiple communities and that the relative importance of the two processes varies consistently along a gradient of environmental filtering.     

Changes in assembly processes in soil bacterial communities following a wildfire disturbance

Although recent work has shown that both deterministic and stochastic processes are important in structuring microbial communities, the factors that affect the relative contributions of niche and neutral processes are poorly understood. The macrobiological literature indicates that ecological disturbances can influence assembly processes. Thus, we sampled bacterial communities at 4 and 16 weeks following a wildfire and used null deviation analysis to examine the role that time since disturbance has in community assembly. Fire dramatically altered bacterial community structure and diversity as well as soil chemistry for both time-points. Community structure shifted between 4 and 16 weeks for both burned and unburned communities. Community assembly in burned sites 4 weeks after fire was significantly more stochastic than in unburned sites. After 16 weeks, however, burned communities were significantly less stochastic than unburned communities. Thus, we propose a three-phase model featuring shifts in the relative importance of niche and neutral processes as a function of time since disturbance. Because neutral processes are characterized by a decoupling between environmental parameters and community structure, we hypothesize that a better understanding of community assembly may be important in determining where and when detailed studies of community composition are valuable for predicting ecosystem function.     

The Impact of Anthropogenic Disturbance on Assembly Patterns of Termite Communities

Intensification of land‐use threatens biodiversity, especially in tropical ecosystems that harbor the planet's highest species richness. This negative impact of anthropogenic disturbance on species numbers is well established, but the mechanisms underlying the community assembly processes are less well understood. Termites are of fundamental importance in tropical ecosystems where they are critical for nutrient recycling and species diversity. We tested the impact of anthropogenic disturbance on termite species diversity and assembly processes in a West African savanna applying the newest techniques of phylogenetic community analyses. Species richness dropped in areas of intensive land‐use and compositional similarity between intensive land‐use areas was high. This contrasted with a protected National Park where communities were characterized by high species richness and intermediate species turnover between sites. Slightly disturbed areas in the buffer zone surrounding the park were intermediate, they still had high species richness but similarity between sites increased. Strikingly, the assembly pattern also changed with disturbance from more phylogenetic overdispersion to more clustering (coexisting species became phylogenetically more similar), but only when the fungus‐growing termite Macrotermes bellicosus was absent. Our data suggest that the major forces structuring termite communities depend: (1) on the presence of this dominant mound‐building termite; and (2) that they change to more environmental filtering with disturbance. Anthropogenic disturbance seems to function as a filter that allows only a specific subset of species to occur. Such an effect might be widespread in ecology but it is difficult to document quantitatively. Phylogenetic community analyses can help to contribute such evidence.     

Non-Native Plants Disrupt Dual Promotion of Native Alpha and Beta Diversity

Non-native species can alter patterns of species diversity at multiple spatial scales, but the processes that underlie multi-scale effects remain unclear. Here we show that non-native species reduce native diversity at multiple scales through simultaneous disruption of two processes of native community assembly: species immigration, which enhances alpha diversity, and community divergence, which enhances beta diversity. Community divergence refers to the process in which local communities diverge over time in species composition because the history of species immigration and, consequently, the way species affect one another within communities are variable among communities. Continuous experimental removal of species over four years of floodplain succession revealed that, when non-native species were excluded, stochastic variation in the timing of a dominant native species’ arrival allowed local communities to diverge, thereby enhancing beta diversity, without compromising promotion of alpha diversity by species immigration. In contrast, when non-native species were allowed to enter experimental plots, they not only reduced native alpha diversity by limiting immigration, but also diminished the dominant native species’ role in enhancing native beta diversity. Our results highlight the importance of community assembly and succession for understanding multi-scale effects of non-native species.     

A conceptual framework for the spatial analysis of functional trait diversity

Functional trait diversity is a popular tool in modern ecology, mainly used to infer assembly processes and ecosystem functioning. Patterns of functional trait diversity are shaped by ecological processes such as environmental filtering, species interactions and dispersal that are inherently spatial, and different processes may operate at different spatial scales. Adding a spatial dimension to the analysis of functional trait diversity may thus increase our ability to infer community assembly processes and to predict change in assembly processes following disturbance or land‐use change. Richness, evenness and divergence of functional traits are commonly used indices of functional trait diversity that are known to respond differently to large‐scale filters related to environmental heterogeneity and dispersal and fine‐scale filters related to species interactions (competition). Recent developments in spatial statistics make it possible to separately quantify large‐scale patterns (variation in local means) and fine‐scale patterns (variation around local means) by decomposing overall spatial autocorrelation quantified by Moran's coefficient into its positive and negative components using Moran eigenvector maps (MEM). We thus propose to identify the spatial signature of multiple ecological processes that are potentially acting at different spatial scales by contrasting positive and negative components of spatial autocorrelation for each of the three indices of functional trait diversity. We illustrate this approach with a case study from riparian plant communities, where we test the effects of disturbance on spatial patterns of functional trait diversity. The fine‐scale pattern of all three indices was increased in the disturbed versus control habitat, suggesting an increase in local scale competition and an overall increase in unexplained variance in the post‐disturbance versus control community. Further research using simulation modeling should focus on establishing the proposed link between community assembly rules and spatial patterns of functional trait diversity to maximize our ability to infer multiple processes from spatial community structure.     

Deterministic processes dominate nematode community structure in the Fynbos Mediterranean heathland of South Africa

In many ecosystems, the factors that determine landscape-scale community structure of soil nematodes are poorly understood. We were interested in discovering whether deterministic or stochastic factors dominate nematode community variation. We used a novel metagenetic approach to investigate variation in nematode community structure in the Fynbos vegetation of South Africa. We compared 23 samples of soil nematode communities from five different Fynbos landscapes. Nematode DNA was 454-pyrosequenced for the 18S rRNA gene. We investigated the community structure, diversity, and the relative role of both deterministic (niche-based) and neutral processes may play in delimiting the nematode phylogenetic community structure in different Fynbos vegetation types. Nematode diversity showed no relationship to any measured soil parameter. The phylogenetic signal showed that more closely related types of nematodes in the Fynbos tended to co-occur more often than would be expected by chance, demonstrating that (a) closely related lineages occupy similar niches spatially and (b) community variation is influenced more by determinism than stochasticity. The standardised beta mean nearest taxon distance (ses.βMNTD) index showed no association with vegetation type. Both ses.NTI (nearest taxon index) and ses.βMNTD deviated significantly from null models, indicating that deterministic processes were important in the assembly of nematode communities. Furthermore, at local scale, the ses.NTI was significantly higher than null expectations, indicating that co-occurrence of related nematode lineages is determined by the differences in environmental conditions across the sites. We conclude that in the Fynbos there is niche overlap between closely related types of nematodes, that nematode speciation tends to occur conservatively into closely related niches, and that the phylogenetic community structure reveals that deterministic (rather than stochastic) processes are more important in delimiting the community assembly.     

The scale dependence of fungal community distribution in paddy soil driven by stochastic and deterministic processes

Fungal communities play important roles in terrestrial ecosystem functioning. Unraveling the relative importance of stochastic versus deterministic processes in shaping biogeographic patterns of fungal communities has long been a challenge in microbial ecology, owing to high biodiversity and difficulties in identifying fungal taxa. Using a unique anthropogenic system of geographically isolated paddy ‘islands’, we collected 198 soil samples with a spatially explicit design to examine how ecological processes shaped fungal biogeographic patterns. Fungal community structure showed scale-dependent distance-decay relationships. Stochastic processes (dispersal and drift) contributed more to community assembly than deterministic processes (selection) at the local scale, which was largely attributed to drift. In contrast, deterministic processes contributed more to community assembly than stochastic processes at the regional scale, with soil dissolved organic carbon being the most important measured factor. Collectively, scale dependence of fungal biogeographical patterns in paddy soils is influenced by differential contribution of deterministic and stochastic processes.     

Predators temper the relative importance of stochastic processes in the assembly of prey metacommunities

Communities assemble through a combination of stochastic processes, which can make environmentally similar communities divergent (high β-diversity), and deterministic processes, which can make environmentally similar communities convergent (low β-diversity). Top predators can influence both stochasticity (e.g. colonization and extinction events) and determinism (e.g. size of the realized species pool), in community assembly, and thus their net effect is unknown. We investigated how predatory fish influenced the scaling of prey diversity in ponds at local and regional spatial scales. While fish reduced both local and regional richness, their effects were markedly more intense at the regional scale. Underlying this result was that the presence of fish made localities within metacommunities more similar in their community composition (lower β-diversity), suggesting that fish enhance the deterministic, relative to the stochastic, components of community assembly. Thus, the presence of predators can alter fundamental mechanisms of community assembly and the scaling of diversity within metacommunities.     

Species abundance distributions as a proxy for the niche-neutrality continuum

Aims Species abundance distributions (SADs) are often used to verify mechanistic theories underlying community assembly. However, it is now accepted that SADs alone are not sufficient to reveal biological mechanisms. Recent attention focuses on the relative importance of stochastic dispersal processes versus deterministic processes such as interspecific competition and environmental filtering. Here, we combine a study of the commonness and rarity of species (i.e. the SAD) with mechanistic processes underlying community composition. By comparing the occurrence frequencies of each and every species with its abundance, we quantify the relative contributions of common and rare species to the maintenance of community structure. Essentially, we relate the continuum between commonness and rarity with that of niches and neutrality.Methods An individual-based, spatially explicit model was used to simulate local communities in niche spaces with the same parameters. We generated sets of assemblages from which species were eliminated in opposing sequences: from common to rare and from rare to common, and investigated the relationship between the abundance and frequency of species. We tested the predictions of our model with empirical data from a field experiment in the environmentally homogeneous alpine meadows of the Qinghai–Tibetan plateau.Important findings Our simulations support the widespread notion that common species maintain community structure, while rare species maintain species diversity, in both local and regional communities. Our results, both from theoretical simulations and from empirical observations, revealed positive correlations between the abundance of a particular species and its occurrence frequency. SAD curves describe a continuum between commonness and rarity. Removing species from the 'rare' end of this continuum has little effect on the similarity of communities, but removing species from the 'common' end of the continuum causes significant increases in beta diversity, or species turnover, between communities. In local communities distributed in a homogenous habitat, species located at the 'common' end of the continuum should be selected by environmental filtering, with niche space partitioning governed by interspecific competition. Conversely, species located at the 'rare' end of the continuum are most likely subject to stochastic dispersal processes. Species situated at intermediate locations on this continuum are therefore determined by niche and neutral processes acting together. Our results suggest that, in homogeneous habitats, SAD curves describing the common-rare continuum may also be used to describe the continuum between niches and neutrality.     

Neutral processes and species sorting in benthic microalgal community assembly: effects of tidal resuspension

Benthic microalgae (BMA) provide vital food resources for heterotrophs and stabilize sediments with their extracellular secretions. A central goal in ecology is to understand how processes such as species interactions and dispersal, contribute to observed patterns of species abundance and distribution. Our objectives were to assess the effects of sediment resuspension on microalgal community structure. We tested whether taxa‐abundance distributions could be predicted using neutral community models (NCMs) and also specific hypotheses about passive migration: (i) As migration decreases in sediment patches, BMA α‐diversity will decrease, and (ii) As migration decreases, BMA community dissimilarity (β‐diversity) will increase. Co‐occurrence indices (checkerboard score and variance ratio) were also computed to test for deterministic factors, such as competition and niche differentiation, in shaping communities. Two intertidal sites (mudflat and sand bar) differing in resuspension regime were sampled throughout the tidal cycle. Fluorometry and denaturing gradient gel electrophoresis were utilized to investigate diatom community structure. Observed taxa‐abundances fit those predicted from NCMs reasonably well (R² of 0.68–0.93), although comparisons of observed local communities to artificial randomly assembled communities rejected the null hypothesis that diatom communities were assembled solely by stochastic processes. No co‐occurrence tests indicated a significant role for competitive exclusion or niche partitioning in microalgal community assembly. In general, predictions about relationships between migration and species diversity were supported for local community dynamics. BMA at low tide (lowest migration) exhibited reduced α‐diversity as compared to periods of immersion at both mudflat and sand bar sites. β‐diversity was higher during low tide emersion on the mudflat, but did not differ temporally at the sand bar site. In between‐site metacommunity comparisons, low‐ and high‐resuspension sites exhibited distinct community compositions while the low‐energy mudflats contained higher microalgal biomass and greater α‐diversity. To our knowledge this is the first study to test the relevance of neutral processes in structuring marine microalgal communities. Our results demonstrate a prominent role for stochastic factors in structuring local BMA community assembly, although unidentified nonrandom processes also appear to play some role. High passive migration, in particular, appears to help maintain species diversity and structure communities in both sand and muddy habitats.     

Relative Roles of Deterministic and Stochastic Processes in Driving the Vertical Distribution of Bacterial Communities in a Permafrost Core from the Qinghai-Tibet Plateau,China

Understanding the processes that influence the structure of biotic communities is one of the major ecological topics, and both stochastic and deterministic processes are expected to be at work simultaneously in most communities. Here, we investigated the vertical distribution patterns of bacterial communities in a 10-m-long soil core taken within permafrost of the Qinghai-Tibet Plateau. To get a better understanding of the forces that govern these patterns, we examined the diversity and structure of bacterial communities, and the change in community composition along the vertical distance (spatial turnover) from both taxonomic and phylogenetic perspectives. Measures of taxonomic and phylogenetic beta diversity revealed that bacterial community composition changed continuously along the soil core, and showed a vertical distance-decay relationship. Multiple stepwise regression analysis suggested that bacterial alpha diversity and phylogenetic structure were strongly correlated with soil conductivity and pH but weakly correlated with depth. There was evidence that deterministic and stochastic processes collectively drived bacterial vertically-structured pattern. Bacterial communities in five soil horizons (two originated from the active layer and three from permafrost) of the permafrost core were phylogenetically random, indicator of stochastic processes. However, we found a stronger effect of deterministic processes related to soil pH, conductivity, and organic carbon content that were structuring the bacterial communities. We therefore conclude that the vertical distribution of bacterial communities was governed primarily by deterministic ecological selection, although stochastic processes were also at work. Furthermore, the strong impact of environmental conditions (for example, soil physicochemical parameters and seasonal freeze-thaw cycles) on these communities underlines the sensitivity of permafrost microorganisms to climate change and potentially subsequent permafrost thaw.     

Rare taxa maintain the stability of crop mycobiomes and ecosystem functions

Plants harbour highly diverse mycobiomes which sustain essential functions for host health and productivity. However, ecological processes that govern the plant–mycobiome assembly, interactions and their impact on ecosystem functions remain poorly known. Here we characterized the ecological role and community assembly of both abundant and rare fungal taxa along the soil–plant continuums (rhizosphere, phyllosphere and endosphere) in the maize–wheat/barley rotation system under different fertilization practices at two contrasting sites. Our results indicate that mycobiome assembly is shaped predominantly by compartment niche and host species rather than by environmental factors. Moreover, crop-associated fungal communities are dominated by few abundant taxa mainly belonging to Sordariomycetes and Dothideomycetes, while the majority of diversity within mycobiomes are represented by rare taxa. For plant compartments, the abundant sub-community is mainly determined by stochastic processes. In contrast, the rare sub-community is more sensitive to host selection and mainly governed by deterministic processes. Furthermore, our results demonstrate that rare taxa play an important role in fungal co-occurrence network and ecosystem functioning like crop yield and soil enzyme activities. These results significantly advance our understanding of crop mycobiome assembly and highlight the key role of rare taxa in sustaining the stability of crop mycobiomes and ecosystem functions.     

Phylogenetic distance–decay patterns are not explained by local community assembly processes in freshwater lake microbial communities

While water and sediment microbial communities exhibit pronounced spatio-temporal patterns in freshwater lakes, the underlying drivers are yet poorly understood. Here, we evaluated the importance of spatial and temporal variation in abiotic environmental factors for bacterial and microeukaryotic community assembly and distance–decay relationships in water and sediment niches in Hongze Lake. By sampling across the whole lake during both Autumn and Spring sampling time points, we show that only bacterial sediment communities were governed by deterministic community assembly processes due to abiotic environmental drivers. Nevertheless, consistent distance–decay relationships were found with both bacterial and microeukaryotic communities, which were relatively stable with both sampling time points. Our results suggest that spatio-temporal variation in environmental factors was important in explaining mainly bacterial community assembly in the sediment, possibly due lesser disturbance. However, clear distance–decay patterns emerged also when the community assembly was stochastic. Together, these results suggest that abiotic environmental factors do not clearly drive the spatial structuring of lake microbial communities, highlighting the need to understand the role of other potential drivers, such as spatial heterogeneity and biotic species interactions.     

Deterministic tropical tree community turnover: evidence from patterns of functional beta diversity along an elevational gradient

Explaining the mechanisms that produce the enormous diversity within and between tropical tree communities is a pressing challenge for plant community ecologists. Mechanistic hypotheses range from niche-based deterministic to dispersal-based stochastic models. Strong tests of these hypotheses require detailed information regarding the functional strategies of species. A few tropical studies to date have examined trait dispersion within individual forest plots using species trait means in order to ask whether coexisting species tend to be more or less functionally similar than expected given a null model. The present work takes an alternative approach by: (i) explicitly incorporating population-level trait variability; and (ii) quantifying the functional beta diversity in a series of 15 tropical forest plots arrayed along an elevational gradient. The results show a strong pattern of decay in community functional similarity with elevation. These observed patterns of functional beta diversity are shown to be highly non-random and support a deterministic model of tropical tree community assembly and turnover.     

Quantifying the role of deterministic assembly and stochastic drift in a natural community of Arctic mosses

The interpretation of natural plant communities frequently invokes species‐sorting controlled by niche differences along spatial environmental gradients. This process of niche structuring can be explained by reference to functional traits, which provide a mechanistic explanation for community structure. In contrast, models explaining species coexistence obviate the limiting effect of niche difference, by invoking processes which cause species‐level drift, e.g. demographic stochasticity. This paper investigates a simple habitat with strong gradients (moss communities in a patterned arctic wetland) to identify signature‐patterns under‐pinning the relative importance of deterministic assembly and stochastic drift in a natural community. First, ordination analysis was used to confirm community composition structured by a range of nine carefully selected functional traits. Second, to determine whether traits explaining community composition might also explain species richness, local species richness (sR) was compared to (1) observed trait diversity and (2) expected trait diversity based on permutation tests, which are used to simulate null community assembly for different values of sR. Traits explaining species composition, consistent with deterministic niche structuring, do not appear to maintain sR. This surprising result was explained by decomposing the community into individual pair‐wise comparisons, i.e. species niche‐differences and association (χ²). Results support deterministic processes via the sorting of species with similar and contrasting niches, at opposite ends of a composite environmental gradient. Nevertheless, stochastic drift is apparent in the random structure of a majority of pair‐wise associations; in addition, a species’ abundance was in general not related to environmental distance from response‐optima. We suggest therefore that spatial pattern in the moss community is a balance between deterministic forces with respect to species traits and controlling environmental gradients, and stochastic drift, which weakens this deterministic structure.     

Compositional divergence and convergence in local communities and spatially structured landscapes

Community structure depends on both deterministic and stochastic processes. However, patterns of community dissimilarity (e.g. difference in species composition) are difficult to interpret in terms of the relative roles of these processes. Local communities can be more dissimilar (divergence) than, less dissimilar (convergence) than, or as dissimilar as a hypothetical control based on either null or neutral models. However, several mechanisms may result in the same pattern, or act concurrently to generate a pattern, and much research has recently been focusing on unravelling these mechanisms and their relative contributions. Using a simulation approach, we addressed the effect of a complex but realistic spatial structure in the distribution of the niche axis and we analysed patterns of species co-occurrence and beta diversity as measured by dissimilarity indices (e.g. Jaccard index) using either expectations under a null model or neutral dynamics (i.e., based on switching off the niche effect). The strength of niche processes, dispersal, and environmental noise strongly interacted so that niche-driven dynamics may result in local communities that either diverge or converge depending on the combination of these factors. Thus, a fundamental result is that, in real systems, interacting processes of community assembly can be disentangled only by measuring traits such as niche breadth and dispersal. The ability to detect the signal of the niche was also dependent on the spatial resolution of the sampling strategy, which must account for the multiple scale spatial patterns in the niche axis. Notably, some of the patterns we observed correspond to patterns of community dissimilarities previously observed in the field and suggest mechanistic explanations for them or the data required to solve them. Our framework offers a synthesis of the patterns of community dissimilarity produced by the interaction of deterministic and stochastic determinants of community assembly in a spatially explicit and complex context.     

Dispersal and neutral sampling mediate contingent effects of disturbance on plant beta‐diversity: a meta‐analysis

A major challenge in ecology, conservation and global‐change biology is to understand why biodiversity responds differently to similar environmental changes. Contingent biodiversity responses may depend on how disturbance and dispersal interact to alter variation in community composition (β‐diversity) and assembly mechanisms. However, quantitative syntheses of these patterns and processes across studies are lacking. Using null‐models and meta‐analyses of 22 factorial experiments in herbaceous plant communities across Europe and North America, we show that disturbance diversifies communities when dispersal is limited, but homogenises communities when combined with increased immigration from the species pool. In contrast to the hypothesis that disturbance and dispersal mediate the strength of niche assembly, both processes altered β‐diversity through neutral‐sampling effects on numbers of individuals and species in communities. Our synthesis suggests that stochastic effects of disturbance and dispersal on community assembly play an important, but underappreciated, role in mediating biotic homogenisation and biodiversity responses to environmental change.     

Stochastic processes play more important roles in driving the dynamics of rarer species

Aims One major goal of modern community ecology is to understand how deterministic and stochastic processes combine to drive community assembly. However, little empirical knowledge is known about how their relative importance varies between common and rare species.Methods We exploited two 30-year data sets of plant communities in a temperate steppe using two different methods. One is a null model method, and the other is a recently developed direct-calculation method.Important findings We found that stochastic processes tended to be more important in influencing rare than common species. This finding suggests that stochastic forces may play a more important role in structuring communities with more rare species, providing a possible solution to the debate on the varied importance of deterministic and stochastic processes among different communities.     

Interactive effects of disturbance and dispersal on community assembly

The traditional debate on alternative community states has been over whether or not they exist. Studies of community assembly have examined the role of assembly history in driving community divergence, but the context in which assembly history becomes important is a continued topic of interest. In this study, we created communities of bacterivorous ciliated protists in laboratory microcosms and manipulated assembly history, disturbance frequency, and the presence of dispersal among local communities to investigate the mechanisms behind community divergence. Specifically, we sought to understand how the role of assembly history changed in response to disturbance, dispersal, and the combination of the two. Assembly history influenced the identity of the dominant species through priority effects, and dispersal and disturbance showed interactive effects on both alpha and beta diversity. Dispersal substantially increased alpha diversity, but only in the absence of disturbance, and it reduced beta diversity, but not in the presence of low or mixed disturbance. These results demonstrate that the role of assembly history and the strength of priority effects depend on community context, suggesting that understanding the interactions between various factors shaping community assembly is important for understanding how ecological communities are structured.     

Deterministic processes have limited impacts on foliar fungal endophyte communities along a savanna-forest successional gradient

Patterns and drivers of succession provide insight into the mechanisms that govern community assembly, but remain poorly understood for microbial communities. We assess whether successional trends of trees are mirrored by foliar endophyte communities of three tree species across a deterministic woody successional gradient. Additionally, we test the relative contribution of abiotic predictors, biotic factors, and spatial distance between sites in predicting composition and richness of endophyte communities. Unlike the tree community, endophyte communities showed no consistent evidence of deterministic succession. Host identity was the most important factor structuring endophyte community composition; within hosts, spatial distance from the indigenous forest and between samples was important, while environmental predictors had small and inconsistent effects. Much variation in endophyte composition remained unexplained. In contrast, endophyte richness was well-explained by predictor variables. Host identity was most important in predicting endophyte richness, while the effect of other predictors on richness differed between host species. We conclude that deterministic succession in trees did not result in deterministic succession in endophyte communities; instead community assembly was most strongly influenced by host identity; while within hosts, neutral processes may be more important for endophyte assembly than deterministic factors.