A guide to null models for animal social network analysis

Null models are an important component of the social network analysis toolbox. However, their use in hypothesis testing is still not widespread. Furthermore, several different approaches for constructing null models exist, each with their relative strengths and weaknesses, and often testing different hypotheses.

In this study, I highlight why null models are important for robust hypothesis testing in studies of animal social networks. Using simulated data containing a known observation bias, I test how different statistical tests and null models perform if such a bias was unknown.

I show that permutations of the raw observational (or ‘pre‐network’) data consistently account for underlying structure in the generated social network, and thus can reduce both type I and type II error rates. However, permutations of pre‐network data remain relatively uncommon in animal social network analysis because they are challenging to implement for certain data types, particularly those from focal follows and GPS tracking.

I explain simple routines that can easily be implemented across different types of data, and supply R code that applies each type of null model to the same simulated dataset. The R code can easily be modified to test hypotheses with empirical data. Widespread use of pre‐network data permutation methods will benefit researchers by facilitating robust hypothesis testing.

Comparing species interaction networks along environmental gradients

Knowledge of species composition and their interactions, in the form of interaction networks, is required to understand processes shaping their distribution over time and space. As such, comparing ecological networks along environmental gradients represents a promising new research avenue to understand the organization of life. Variation in the position and intensity of links within networks along environmental gradients may be driven by turnover in species composition, by variation in species abundances and by abiotic influences on species interactions. While investigating changes in species composition has a long tradition, so far only a limited number of studies have examined changes in species interactions between networks, often with differing approaches. Here, we review studies investigating variation in network structures along environmental gradients, highlighting how methodological decisions about standardization can influence their conclusions. Due to their complexity, variation among ecological networks is frequently studied using properties that summarize the distribution or topology of interactions such as number of links, connectance, or modularity. These properties can either be compared directly or using a procedure of standardization. While measures of network structure can be directly related to changes along environmental gradients, standardization is frequently used to facilitate interpretation of variation in network properties by controlling for some co-variables, or via null models. Null models allow comparing the deviation of empirical networks from random expectations and are expected to provide a more mechanistic understanding of the factors shaping ecological networks when they are coupled with functional traits. As an illustration, we compare approaches to quantify the role of trait matching in driving the structure of plant-hummingbird mutualistic networks, i.e. a direct comparison, standardized by null models and hypothesis-based metaweb. Overall, our analysis warns against a comparison of studies that rely on distinct forms of standardization, as they are likely to highlight different signals. Fostering a better understanding of the analytical tools available and the signal they detect will help produce deeper insights into how and why ecological networks vary along environmental gradients.

Different but equal: the implausible assumption at the heart of neutral theory

1. The core assumption of neutral theory is that all individuals in a community have equal fitness regardless of species, and regardless of the species composition of the community. But, real communities consist of species exhibiting large trait differences; hence these differences must be subject to perfect fitness-equalizing trade-offs for neutrality to hold. 2. Here we explain that perfect equalizing trade-offs are extremely unlikely to occur in reality, because equality of fitness among species is destroyed by: (i) any deviation in the functional form of the trade-off away from the one special form that gives equal fitness; (ii) spatial or temporal variation in performance; (iii) random species differences in performance. 3. In the absence of the density-dependent processes stressed by traditional niche-based community ecology, communities featuring small amounts of (i) or (ii) rapidly lose trait variation, becoming dominated by species with similar traits, and exhibit substantially lower species richness compared to the neutral case. Communities featuring random interspecific variation in traits (iii) lose all but a few fortuitous species. 4. Thus neutrality should be viewed, a priori, as a highly improbable explanation for the long-term co-occurrence of measurably different species within ecological communities. In contrast, coexistence via niche structure and density dependence, is robust to species differences in baseline fitness, and so remains plausible. 5. We conclude that: (i) co-occurring species will typically exhibit substantial differences in baseline fitness even when (imperfect) equalizing trade-offs have been taken into account; (ii) therefore, communities must be strongly niche structured, otherwise they would lose both trait variation and species richness; (iii) nonetheless, even in strongly niche-structured communities, it is possible that the abundance of species with similar traits are at least partially free to drift.

Selective extinction drives taxonomic and functional alpha and beta diversities in island bird assemblages

Taxonomic diversity considers all species being equally different from each other and thus disregards species' different ecological functions. Exploring taxonomic and functional aspects of biodiversity simultaneously can better understand the processes of community assembly. We analysed taxonomic and functional alpha and beta diversities of breeding bird assemblages on land-bridge islands in the Thousand Island Lake, China. Given the high dispersal ability of most birds at this spatial scale (several kilometres), we predicted (i) selective extinction driving alpha and beta diversities after the creation of land-bridge islands of varying area and (ii) low taxonomic and functional beta diversities that were not correlated to spatial distance. Breeding birds were surveyed on 37 islands annually from 2007 to 2014. We decomposed beta diversity of breeding birds into spatial turnover and nestedness-resultant components, and related taxonomic and functional diversities to island area and isolation using power regression models (for alpha diversity) and multiple regression models on distance matrices (for beta diversity). We then ran simulations to assess the strength of the correlations between taxonomic and functional diversities. Results revealed that both taxonomic and functional alpha diversities increased with island area. The taxonomic nestedness-resultant and turnover components increased and decreased with difference in area, respectively, but functional counterparts did not. Isolation played a minor role in explaining alpha- and beta-diversity patterns. By partitioning beta diversity, we found low levels of overall taxonomic and functional beta diversities. The functional nestedness-resultant component dominated overall functional beta diversity, whereas taxonomic turnover was the dominant component for taxonomic beta diversity. The simulation showed that functional alpha and beta diversities were significantly correlated with taxonomic diversities, and the observed values of correlations were significantly different from null expectations of random extinction. Our assessment of island bird assemblages validated the predictions of no distance effects and low beta diversity due to pervasive dispersal events among islands and also suggested that selective extinction drives taxonomic and functional alpha and beta diversities. The contrasting turnover and nestedness-resultant components of taxonomic and functional beta diversities demonstrate the importance of considering the multifaceted nature of biodiversity when examining community assembly.

Interaction frequency, network position, and the temporal persistence of interactions in a plant-pollinator network

Ecological interactions are highly dynamic in time and space. Previous studies of plant-animal mutualistic networks have shown that the occurrence of interactions varies substantially across years. We analyzed interannual variation of a quantitative mutualistic network, in which links are weighted by interaction frequency. The network was sampled over six consecutive years, representing one of the longest time series for a community-wide mutualistic network. We estimated the interannual similarity in interactions and assessed the determinants of their persistence. The occurrence of interactions varied greatly among years, with most interactions seen in only one year (64%) and few (20%) in more than two years. This variation was associated with the frequency and position of interactions relative to the network core, so that the network consisted of a persistent core of frequent interactions and many peripheral, infrequent interactions. Null model analyses suggest that species abundances play a substantial role in generating these patterns. Our study represents an important step in the study of ecological networks, furthering our mechanistic understanding of the ecological processes driving the temporal persistence of interactions.

Relaxed Selection and the Rapid Evolution of Reproductive Genes

Evolutionary genomic studies find that reproductive protein genes, those directly involved in reproductive processes, diversify more rapidly than most other gene categories. Strong postcopulatory sexual selection acting within species is the predominant hypothesis proposed to account for the observed pattern. Recently, relaxed selection due to sex-specific gene expression has also been put forward to explain the relatively rapid diversification. We contend that relaxed selection due to sex-limited gene expression is the correct null model for tests of molecular evolution of reproductive genes and argue that it may play a more significant role in the evolutionary diversification of reproductive genes than previously recognized. We advocate for a re-evaluation of adaptive explanations for the rapid diversification of reproductive genes.

A new dynamic null model for phylogenetic community structure

Phylogenies are increasingly applied to identify the mechanisms structuring ecological communities but progress has been hindered by a reliance on statistical null models that ignore the historical process of community assembly. Here, we address this, and develop a dynamic null model of assembly by allopatric speciation, colonisation and local extinction. Incorporating these processes fundamentally alters the structure of communities expected due to chance, with speciation leading to phylogenetic overdispersion compared to a classical statistical null model assuming equal probabilities of community membership. Applying this method to bird and primate communities in South America we show that patterns of phylogenetic overdispersion – often attributed to negative biotic interactions – are instead consistent with a species neutral model of allopatric speciation, colonisation and local extinction. Our findings provide a new null expectation for phylogenetic community patterns and highlight the importance of explicitly accounting for the dynamic history of assembly when testing the mechanisms governing community structure.

An integrated framework to improve the concept of resource specialisation

Resource specialisation, although a fundamental component of ecological theory, is employed in disparate ways. Most definitions derive from simple counts of resource species. We build on recent advances in ecophylogenetics and null model analysis to propose a concept of specialisation that comprises affinities among resources as well as their co-occurrence with consumers. In the distance-based specialisation index (DSI), specialisation is measured as relatedness (phylogenetic or otherwise) of resources, scaled by the null expectation of random use of locally available resources. Thus, specialists use significantly clustered sets of resources, whereas generalists use over-dispersed resources. Intermediate species are classed as indiscriminate consumers. The effectiveness of this approach was assessed with differentially restricted null models, applied to a data set of 168 herbivorous insect species and their hosts. Incorporation of plant relatedness and relative abundance greatly improved specialisation measures compared to taxon counts or simpler null models, which overestimate the fraction of specialists, a problem compounded by insufficient sampling effort. This framework disambiguates the concept of specialisation with an explicit measure applicable to any mode of affinity among resource classes, and is also linked to ecological and evolutionary processes. This will enable a more rigorous deployment of ecological specialisation in empirical and theoretical studies.

Neither host-specific nor random: vascular epiphytes on three tree species in a Panamanian lowland forest

BACKGROUND AND AIMS

A possible role of host tree identity in the structuring of vascular epiphyte communities has attracted scientific attention for decades. Specifically, it has been suggested that each host tree species has a specific subset of the local species pool according to its own set of properties, e.g. physicochemical characteristics of the bark, tree architecture, or leaf phenology patterns.

METHODS

A novel, quantitative approach to this question is presented, taking advantage of a complete census of the vascular epiphyte community in 0.4 ha of undisturbed lowland forest in Panama. For three locally common host-tree species (Socratea exorrhiza, Marila laxiflora, Perebea xanthochyma) null models were created of the expected epiphyte assemblages assuming that epiphyte colonization reflected random distribution of epiphytes in the forest.

KEY RESULTS

In all three tree species, abundances of the majority of epiphyte species (69-81 %) were indistinguishable from random, while the remaining species were about equally over- or under-represented compared with their occurrence in the entire forest plot. Permutations based on the number of colonized trees (reflecting observed spatial patchiness) yielded similar results. Finally, a third analysis (canonical correspondence analysis) also confirmed host-specific differences in epiphyte assemblages. In spite of pronounced preferences of some epiphytes for particular host trees, no epiphyte species was restricted to a single host.

CONCLUSIONS

The epiphytes on a given tree species are not simply a random sample of the local species pool, but there are no indications of host specificity either.

Functional and phylogenetic structure of island bird communities

Biodiversity change in anthropogenically transformed habitats is often nonrandom, yet the nature and importance of the different mechanisms shaping community structure are unclear. Here, we extend the classic Theory of Island Biogeography (TIB) to account for nonrandom processes by incorporating species traits and phylogenetic relationships into a study of faunal relaxation following habitat loss and fragmentation. Two possible mechanisms can create nonrandom community patterns on fragment islands. First, small and isolated islands might consist of similar or closely related species because they are environmentally homogeneous or select for certain shared traits, such as dispersal ability. Alternatively, communities on small islands might contain more dissimilar or distantly related species than on large islands because limited space and resource availability result in greater competitive exclusion among species with high niche overlap. Breeding birds were surveyed on 36 islands and two mainland sites annually from 2010 to 2014 in the Thousand Island Lake region, China. We assessed community structure of breeding birds on these subtropical land-bridge islands by integrating species' trait and evolutionary distances. We additionally analysed habitat heterogeneity and variance in size ratios to distinguish biotic and abiotic processes of community assembly. Results showed that functional-phylogenetic diversity increased with island area, and decreased with isolation. Bird communities on the mainland were more diverse and generally less clustered than island bird communities and not different than randomly assembled communities. Bird communities on islands tend to be functionally similar and phylogenetically clustered, especially on small and isolated islands. The nonrandom decline in species diversity and change in bird community structure with island area and isolation, along with the relatively homogeneous habitats on small islands, support the environmental filtering hypothesis. Our study demonstrates the importance of integrating multiple forms of diversity for understanding the effects of habitat loss and fragmentation, and further reveals that TIB could be extended to community measures by moving beyond assumptions of species equivalency in colonisation rates and extinction susceptibilities.

What drives the spatial distribution and dynamics of local species richness in tropical forest?

Understanding the structure and dynamics of highly diverse tropical forests is challenging. Here we investigate the factors that drive the spatio-temporal variation of local tree numbers and species richness in a tropical forest (including 1250 plots of 20 × 20 m²). To this end, we use a series of dynamic models that are built around the local spatial variation of mortality and recruitment rates, and ask which combination of processes can explain the observed spatial and temporal variation in tree and species numbers. We find that processes not included in classical neutral theory are needed to explain these fundamental patterns of the observed local forest dynamics. We identified a large spatio-temporal variability in the local number of recruits as the main missing mechanism, whereas variability of mortality rates contributed to a lesser extent. We also found that local tree numbers stabilize at typical values which can be explained by a simple analytical model. Our study emphasized the importance of spatio-temporal variability in recruitment beyond demographic stochasticity for explaining the local heterogeneity of tropical forests.

Estimating species - area relationships by modeling abundance and frequency subject to incomplete sampling

Models and data used to describe species-area relationships confound sampling with ecological process as they fail to acknowledge that estimates of species richness arise due to sampling. This compromises our ability to make ecological inferences from and about species-area relationships. We develop and illustrate hierarchical community models of abundance and frequency to estimate species richness. The models we propose separate sampling from ecological processes by explicitly accounting for the fact that sampled patches are seldom completely covered by sampling plots and that individuals present in the sampling plots are imperfectly detected. We propose a multispecies abundance model in which community assembly is treated as the summation of an ensemble of species-level Poisson processes and estimate patch-level species richness as a derived parameter. We use sampling process models appropriate for specific survey methods. We propose a multispecies frequency model that treats the number of plots in which a species occurs as a binomial process. We illustrate these models using data collected in surveys of early-successional bird species and plants in young forest plantation patches. Results indicate that only mature forest plant species deviated from the constant density hypothesis, but the null model suggested that the deviations were too small to alter the form of species-area relationships. Nevertheless, results from simulations clearly show that the aggregate pattern of individual species density-area relationships and occurrence probability-area relationships can alter the form of species-area relationships. The plant community model estimated that only half of the species present in the regional species pool were encountered during the survey. The modeling framework we propose explicitly accounts for sampling processes so that ecological processes can be examined free of sampling artefacts. Our modeling approach is extensible and could be applied to a variety of study designs and allows the inclusion of additional environmental covariates.

Weak spatiotemporal response of prey to predation risk in a freely interacting system

The extent to which prey space use actively minimizes predation risk continues to ignite controversy. Methodological reasons that have hindered consensus include inconsistent measurements of predation risk, biased spatiotemporal scales at which responses are measured and lack of robust null expectations.

We addressed all three challenges in a comprehensive analysis of the spatiotemporal responses of adult female elk (Cervus elaphus) to the risk of predation by wolves (Canis lupus) during winter in northern Yellowstone, USA.

We quantified spatial overlap between the winter home ranges of GPS‐collared elk and three measures of predation risk: the intensity of wolf space use, the distribution of wolf‐killed elk and vegetation openness. We also assessed whether elk varied their use of areas characterized by more or less predation risk across hours of the day, and estimated encounter rates between simultaneous elk and wolf pack trajectories. We determined whether observed values were significantly lower than expected if elk movements were random with reference to predation risk using a null model approach.

Although a small proportion of elk did show a tendency to minimize use of open vegetation at specific times of the day, overall we highlight a notable absence of spatiotemporal response by female elk to the risk of predation posed by wolves in northern Yellowstone.

Our results suggest that predator–prey interactions may not always result in strong spatiotemporal patterns of avoidance.

Statistical testing in transcriptomic-neuroimaging studies: A how-to and evaluation of methods assessing spatial and gene specificity

Multiscale integration of gene transcriptomic and neuroimaging data is becoming a widely used approach for exploring the molecular underpinnings of large‐scale brain organization in health and disease. Proper statistical evaluation of determined associations between imaging‐based phenotypic and transcriptomic data is key in these explorations, in particular to establish whether observed associations exceed “chance level” of random, nonspecific effects. Recent approaches have shown the importance of statistical models that can correct for spatial autocorrelation effects in the data to avoid inflation of reported statistics. Here, we discuss the need for examination of a second category of statistical models in transcriptomic‐neuroimaging analyses, namely those that can provide “gene specificity.” By means of a couple of simple examples of commonly performed transcriptomic‐neuroimaging analyses, we illustrate some of the potentials and challenges of transcriptomic‐imaging analyses, showing that providing gene specificity on observed transcriptomic‐neuroimaging effects is of high importance to avoid reports of nonspecific effects. Through means of simulations we show that the rate of reported nonspecific effects (i.e., effects that cannot be specifically linked to a specific gene or gene‐set) can run as high as 60%, with only less than 5% of transcriptomic‐neuroimaging associations observed through ordinary linear regression analyses showing both spatial and gene specificity. We provide a discussion, a tutorial, and an easy‐to‐use toolbox for the different options of null models in transcriptomic‐neuroimaging analyses.

Small-scale spatial variability in phylogenetic community structure during early plant succession depends on soil properties

During early plant succession, the phylogenetic structure of a community changes in response to important environmental filters and emerging species interactions. We traced the development of temperate-zone plant communities during the first 7 years of primary succession on catchment soils to explore patterns of initial species assembly. We found pronounced small-scale differences in the phylogenetic composition of neighbouring plant assemblages and a large-scale trend towards phylogenetic evenness. This small-scale variability appears to be mediated by soil properties, particularly carbonate content. Therefore, abiotic environmental conditions might counteract or even supersede the effects of interspecific competition among closely related species, which are usually predicted to exhibit patterns of phylogenetic evenness. We conclude that theories on phylogenetic community composition need to incorporate effects of small-scale variability of environmental factors.

Higher β-diversity observed for herbs over woody plants is driven by stronger habitat filtering in a tropical understory

Herbaceous plants are a key component of tropical forests. Previous work indicates that herbs contribute substantially to the species richness of tropical plant communities. However, the processes structuring tropical herb diversity, and how they contrast with woody communities, have been underexplored. Within the understory of a 50-ha forest dynamics plot in central Panama, we compared the diversity, distribution, and abundance of vascular herbaceous plants with woody seedlings (i.e., tree and lianas <1 cm DBH and ≥20 cm tall). Beta-diversity was calculated for each community using a null model approach. We then assessed the similarity in alpha and beta-diversity among herbs, tree seedlings, and liana seedlings. Strengths of habitat associations were measured using permutational ANOVA among topographic habitat-types. Variance partitioning was then used to quantify the amount of variation in species richness and composition explained by spatial and environmental variables (i.e., topography, soils, and shade) for each growth form. Species richness and diversity were highest for tree seedlings, followed by liana seedlings and then herbs. In contrast, beta-diversity was 16-127% higher for herbs compared to woody seedlings, indicating higher spatial variation in this stratum. We observed no correlation between local richness or compositional uniqueness of herbs and woody seedlings across sites, indicating that different processes control the spatial patterns of woody and herbaceous diversity and composition. Habitat associations were strongest for herbs, as indicated by greater compositional dissimilarity among habitat types. Likewise, environmental variables explained a larger proportion of the variation in species richness and composition for herbs than for woody seedlings (richness = 25%, 14%, 12%; composition = 25%, 9%, 6%, for herbs, trees, and lianas, respectively). These differences between strata did not appear to be due to differences in lifespan alone, based on data from adult trees. Our results point to contrasting assembly mechanisms for herbaceous and woody communities, with herbs showing stronger niche-derived structure. Future research on tropical herbaceous communities is likely to yield new insights into the many processes structuring diverse plant communities.

Mean reproductive traits of fungal assemblages are correlated with resource availability

Organisms have evolved a fascinating variety of strategies and organs for successful reproduction. Fruit bodies are the reproductive organ of fungi and vary considerably in size and shape among species. Our understanding of the mechanisms underlying the differences in fruit body size among species is still limited. Fruit bodies of saprotrophic fungi are smaller than those of mutualistic ectomycorrhizal fungi. If differences in fruit body size are determined by carbon acquisition, then mean reproductive traits of saprotrophic and ectomycorrhizal fungi assemblages should vary differently along gradients of resource availability as carbon acquisition seems more unpredictable and costly for saprotrophs than for ectomycorrhizal fungi. Here, we used 48 local inventories of fungal fruit bodies (plot size: 0.02 ha each) sampled along a gradient of resource availability (growing stock) across 3 years in the Bavarian Forest National Park in Germany to investigate regional and local factors that might influence the distribution of species with different reproductive traits, particularly fruit body size. As predicted, mean fruit body size of local assemblages of saprotrophic fungi was smaller than expected from the distribution of traits of the regional species pool across central and northern Europe, whereas that of ectomycorrhizal fungi did not differ from random expectation. Furthermore and also as expected, mean fruit body size of assemblages of saprotrophic fungi was significantly smaller than for assemblages of ectomycorrhizal species. However, mean fruit body sizes of not only saprotrophic species but also ectomycorrhizal species increased with resource availability, and the mean number of fruit bodies of both assemblages decreased. Our results indicate that the differences in carbon acquisition between saprotrophs and ectomycorrhizal species lead to differences in basic reproductive strategies, with implications for the breadth of their distribution. However, the differences in resource acquisition cannot explain detailed species distribution patterns at a finer, local scale based on their reproductive traits.

Phylogenetic turnover during subtropical forest succession across environmental and phylogenetic scales

Although spatial and temporal patterns of phylogenetic community structure during succession are inherently interlinked and assembly processes vary with environmental and phylogenetic scales, successional studies of community assembly have yet to integrate spatial and temporal components of community structure, while accounting for scaling issues. To gain insight into the processes that generate biodiversity after disturbance, we combine analyses of spatial and temporal phylogenetic turnover across phylogenetic scales, accounting for covariation with environmental differences. We compared phylogenetic turnover, at the species- and individual-level, within and between five successional stages, representing woody plant communities in a subtropical forest chronosequence. We decomposed turnover at different phylogenetic depths and assessed its covariation with between-plot abiotic differences. Phylogenetic turnover between stages was low relative to species turnover and was not explained by abiotic differences. However, within the late-successional stages, there was high presence-/absence-based turnover (clustering) that occurred deep in the phylogeny and covaried with environmental differentiation. Our results support a deterministic model of community assembly where (i) phylogenetic composition is constrained through successional time, but (ii) toward late succession, species sorting into preferred habitats according to niche traits that are conserved deep in phylogeny, becomes increasingly important.

Wolf in sheep's clothing: Model misspecification undermines tests of the neutral theory for life histories

Understanding the processes behind change in reproductive state along life‐history trajectories is a salient research program in evolutionary ecology. Two processes, state dependence and heterogeneity, can drive the dynamics of change among states. Both processes can operate simultaneously, begging the difficult question of how to tease them apart in practice. The Neutral Theory for Life Histories (NTLH) holds that the bulk of variations in life‐history trajectories is due to state dependence and is hence neutral: Once previous (breeding) state is taken into account, variations are mostly random. Lifetime reproductive success (LRS), the number of descendants produced over an individual's reproductive life span, has been used to infer support for NTLH in natura. Support stemmed from accurate prediction of the population‐level distribution of LRS with parameters estimated from a state dependence model. We show with Monte Carlo simulations that the current reliance of NTLH on LRS prediction in a null hypothesis framework easily leads to selecting a misspecified model, biased estimates and flawed inferences. Support for the NTLH can be spurious because of a systematic positive bias in estimated state dependence when heterogeneity is present in the data but ignored in the analysis. This bias can lead to spurious positive covariance between fitness components when there is in fact an underlying trade‐off. Furthermore, neutrality implied by NTLH needs a clarification because of a probable disjunction between its common understanding by evolutionary ecologists and its translation into statistical models of life‐history trajectories. Irrespective of what neutrality entails, testing hypotheses about the dynamics of change among states in life histories requires a multimodel framework because state dependence and heterogeneity can easily be mistaken for each other.

Assembly Processes and Co-occurrence Patterns of Abundant and Rare Bacterial Community in the Eastern Indian Ocean

Microbial communities are composed of many rare species and a few abundant species. Considering the disproportionate importance of rare species for ecosystem functioning, it is important to understand the mechanisms structuring the rare and abundant components of a diverse community in response to environmental changes. Here, we used a 16S ribosomal RNA gene sequencing approach to investigate the bacterial community diversity in the Eastern Indian Ocean (EIO) during the monsoon and intermonsoon. We employed a phylogenetic null model and network analysis to evaluate the assembly processes and co-occurrence pattern of the microbial community. We found that higher bacterial diversity was detected in the intermonsoon with high temperature and low Chlorophyll a concentrations and N/P ratios. The balance between ecological deterministic processes and stochastic processes varied with seasons in the EIO. Meanwhile, conditionally rare taxa (CRT) were more likely modulated by variable selection processes than always rare taxa (ART) and abundant taxa (AT) (CRT > ART > AT). By linking assembly process and species co-occurrence, we demonstrated that the microbial co-occurrence associations tended to be higher when deterministic processes (mainly variable selection) were weaker. This negative trend was observed in rare species rather than abundant species. The linkage could enhance our understanding of the underlying mechanisms underpinning the generation and maintenance of microbial community diversity.

Environmental filtering and dispersal limitation jointly shaped the taxonomic and phylogenetic beta diversity of natural forests in southern China

AIM

The mechanisms underlying the maintenance of biodiversity remain to be elucidated. Taxonomic diversity alone remains an unresolved issue, especially in terms of the mechanisms of species co-existence. We hypothesized that phylogenetic information could help to elucidate the mechanism of community assembly and the services and functions of ecosystems. The aim of this study was to explore the mechanisms driving floral diversity in subtropical forests and evaluate the relative effects of these mechanisms on diversity variation, by combining taxonomic and phylogenetic information.

LOCATION

We examined 35 1-ha tree stem-mapped plots across eight national nature reserves in Guangxi Zhuang Autonomous Region, China.

TAXON

Trees.

METHODS

We quantified the taxonomic and phylogenetic β-diversity between each pair of plots using the (abundance-based) Rao's quadratic entropy and the (incidence-based) Sørensen dissimilarity indices. Using a null model approach, we compared the observed β-diversity with the expected diversity at random and calculated the standard effect size of the observed β-diversity deviation. Furthermore, we used distance-based redundancy analysis (dbRDA) to partition the variations in taxonomic and phylogenetic observed β-diversity and β-deviation into four parts to assess the environmental and spatial effects.

RESULTS

The taxonomic β-deviation was related to and higher than the phylogenetic β-deviation (r = .74). This indicated that the species turnover between pairwise plots was mainly the turnover of closely related species. Higher taxonomic and phylogenetic β-deviation were mainly concentrated in the pairwise karst and nonkarst forest plots, indicating that the species in karst forests and nonkarst forests were predominantly from distantly related clades. A large proportions of the variation in taxonomic and phylogenetic β-deviation were explained by the joint effect of environmental and spatial variables, while the contribution of environmental variables was greater than that of spatial variables, probably owing to the influence of the sampling scale dependence, integrality of sampling size and species pool, and the unique climatic and geomorphic characteristics.

MAIN CONCLUSIONS

Our study highlights the importance of phylogeny in biodiversity research. The incorporation of taxonomic and phylogenetic information provides a perspective to explore potential underlying mechanisms that have shaped species assemblages and phylogenetic patterns in biodiversity hotspots.

Spatial scale changes the relationship between beta diversity, species richness and latitude

The relationship between β-diversity and latitude still remains to be a core question in ecology because of the lack of consensus between studies. One hypothesis for the lack of consensus between studies is that spatial scale changes the relationship between latitude and β-diversity. Here, we test this hypothesis using tree data from 15 large-scale forest plots (greater than or equal to 15 ha, diameter at breast height ≥ 1 cm) across a latitudinal gradient (3-30^o) in the Asia-Pacific region. We found that the observed β-diversity decreased with increasing latitude when sampling local tree communities at small spatial scale (grain size ≤0.1 ha), but the observed β-diversity did not change with latitude when sampling at large spatial scales (greater than or equal to 0.25 ha). Differences in latitudinal β-diversity gradients across spatial scales were caused by pooled species richness (γ-diversity), which influenced observed β-diversity values at small spatial scales, but not at large spatial scales. Therefore, spatial scale changes the relationship between β-diversity, γ-diversity and latitude, and improving sample representativeness avoids the γ-dependence of β-diversity.

When do we have the power to detect biological interactions in spatial point patterns?

Uncovering the roles of biotic interactions in assembling and maintaining species-rich communities remains a major challenge in ecology. In plant communities, interactions between individuals of different species are expected to generate positive or negative spatial interspecific associations over short distances. Recent studies using individual-based point pattern datasets have concluded that (a) detectable interspecific interactions are generally rare, but (b) are most common in communities with fewer species; and (c) the most abundant species tend to have the highest frequency of interactions. However, it is unclear how the detection of spatial interactions may change with the abundances of each species, or the scale and intensity of interactions. We ask if statistical power is sufficient to explain all three key results.We use a simple two-species model, assuming no habitat associations, and where the abundances, scale and intensity of interactions are controlled to simulate point pattern data. In combination with an approximation to the variance of the spatial summary statistics that we sample, we investigate the power of current spatial point pattern methods to correctly reject the null model of pairwise species independence.We show the power to detect interactions is positively related to both the abundances of the species tested, and the intensity and scale of interactions, but negatively related to imbalance in abundances. Differences in detection power in combination with the abundance distributions found in natural communities are sufficient to explain all the three key empirical results, even if all pairwise interactions are identical. Critically, many hundreds of individuals of both species may be required to detect even intense interactions, implying current abundance thresholds for including species in the analyses are too low. Sy n thesis. The widespread failure to reject the null model of spatial interspecific independence could be due to low power of the tests rather than any key biological process. Since we do not model habitat associations, our results represent a first step in quantifying sample sizes required to make strong statements about the role of biotic interactions in diverse plant communities. However, power should be factored into analyses and considered when designing empirical studies.