Inferring ecological processes from taxonomic, phylogenetic and functional trait β-diversity

The role of stochasticity in fungal community assembly: explaining apparent stochasticity with field experiments

Stochasticity is a main process in community assembly. However, experimental studies rarely target stochasticity in natural communities, and hence experimental validation of stochasticity estimates in observational studies is lacking. Here, we combine experimental and observational data to unravel the role of stochasticity in the assembly of wood-inhabiting fungi. We carried out a replicated field experiment where the natural colonization of a focal fungal species was simulated through inoculation, and the local fungal communities were monitored through DNA metabarcoding before and after the inoculations. The amount of stochasticity in fungal colonization was less pronounced than expected from the amount of unpredictability in observational data, suggesting that stochasticity may play a smaller role in fungal occurrence than previously anticipated, or that it may be a stronger influence in the dispersal and establishment phases than in colonization per se. Stochasticity was more prominent in the initial phase of community succession, with the earliest successional stage involving a higher level of stochasticity than the later stage after 2 years. We conclude that experimentally measuring the role of stochasticity in community assembly is feasible for species-rich communities under natural conditions and highlight the importance of experimentally testing the accuracy of stochasticity estimates based on observational data.

Drivers of the taxonomic and functional structuring of aquatic and terrestrial floodplain bird communities

Context

There has been a limited amount of research which comparatively examines the local and landscape scale ecological determinants of the community structure of both riparian and aquatic bird communities in floodplain ecosystems.

Objectives

Here, we quantified the contribution of local habitat structure, land cover and spatial configuration of the sampling sites to the taxonomical and functional structuring of aquatic and terrestrial bird communities in a relatively intact floodplain of the river Danube, Hungary.

Methods

We used the relative abundance of species and foraging guilds as response variables in partial redundancy analyses to determine the relative importance of each variable group.

Results

Local-scale characteristics of the water bodies proved to be less influential than land cover and spatial variables both for aquatic and terrestrial birds and both for taxonomic and foraging guild structures. Purely spatial variables were important determinants, besides purely environmental and the shared proportion of variation explained by environmental and spatial variables. The predictability of community structuring generally increased towards the lowest land cover measurement scales (i.e., 500, 250 or 125 m radius buffers). Different land cover types contributed at each scale, and their importance depended on aquatic vs terrestrial communities.

Conclusions

These results indicate the relatively strong response of floodplain bird communities to land cover and spatial configuration. They also suggest that dispersal dynamics and mass-effect mechanisms are critically important for understanding the structuring of floodplain bird communities, and should therefore be considered by conservation management strategies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10980-024-01948-3.

Arthropod community composition in urban landscapes is shaped by both environmental filtering and dispersal limitation

We assessed the relative importance of environmental filtering and dispersal limitation in structuring foliage- and ground-dwelling arthropod communities in central North Carolina. We hypothesized that differences in both local environment and dispersal distance between sites would predict differences in community composition, but that dispersal distance would be more important for ground arthropods than for foliage arthropods. We used variance partitioning to quantify the relative explanatory power of differences in the local environment and dispersal distance between sampling plots. For foliage arthropods, the full set of predictors explained 37-39% of the variance in community dissimilarity, while the full set of predictors explained 42-56% of the variance for ground arthropods. Environmental predictors alone consistently explained more variance in community dissimilarity than geographic predictors alone for ground arthropods, but the variance explained by each predictor set alone varied for foliage arthropods depending on the metric of community dissimilarity. Examining communities of multiple arthropod groups across the same spatial gradient highlights the scale-dependence of these processes and illustrates how variation in the environment can alter the relative abundance of specialist and generalist taxa.

Toward an integrative framework for microbial community coalescence

Community coalescence is defined as the mixing of intact ecological communities. From river confluences to fecal microbiota transplantation, community coalescence constitutes a common ecological occurrence affecting natural and engineered microbial systems. In this opinion article, we propose an integrative framework for microbial community coalescence to guide advances in our understanding of this important - yet underexplored - ecological phenomenon. We start by aligning community coalescence with the unified framework of biological invasion and enumerate commonalities and idiosyncrasies between these two analogous processes. Then, we discuss how organismal interactions and cohesive establishment affect coalescence outcomes with direct implications for community functioning. Last, we propose the use of ecological null modeling to study the interplay of ecological processes structuring community reassembly following coalescence.

Homogenization of Functional Diversity of Rotifer Communities in Relation to Eutrophication in an Urban River of North China

Rapid urbanization has triggered nutrient loading, which will inevitably lead to the eutrophication of water bodies and further affect the structure of aquatic populations. At present, eutrophication is a significant challenge for urban aquatic ecosystems. However, we still know little about the correlation between eutrophication in urban rivers and the composition of aquatic functional groups. The effects of urban river eutrophication on rotifer communities were investigated using an annual field survey of the Jinan section of the Xiaoqing River, a typical urban river in northern China. Using functional diversity (FD) and beta diversity, the spatiotemporal variation of the aquatic biological functional groups regime along stretches subject to different eutrophication was investigated. The functional evenness (FEve) and functional divergence (FDiv) decreased significantly with the increment of the trophic level index. Functional diversity exhibits an extremely low level across functional groups, with the richness difference (RichDiff) being an important component. The results indicate that eutrophication led to the homogenization of rotifer communities. This can be attributed to the functional homogenization of the rotifer community in the Jinan section of the Xiaoqing River. The observed homogenization may be due to widely distributed species complementing the ecological niche space. Our findings provide valuable information on the conservation of the urban river under the threat of eutrophication caused by high-intensity human activities.

Distribution and diversity of anaerobic thermophiles and putative anaerobic nickel-dependent carbon monoxide-oxidizing thermophiles in mesothermal soils and sediments

Even though thermophiles are best known from geothermal and other heated systems, numerous studies have demonstrated that they occur ubiquitously in mesothermal and permanently cold soils and sediments. Cultivation based studies of the latter have revealed that the thermophiles within them are mostly spore-forming members of the Firmicutes. Since the geographic distribution of spores is presumably unconstrained by transport through the atmosphere, similar communities (composition and diversity) of thermophiles might be expected to emerge in mesothermal habitats after they are heated. Alternatively, thermophiles might experience environmental selection before or after heating leading to divergent communities. After demonstrating the ubiquity of anaerobic thermophiles and CO uptake in a variety of mesothermal habitats and two hot springs, we used high throughput sequencing of 16S rRNA genes to assess the composition and diversity of populations that emerged after incubation at 60°C with or without headspace CO concentrations of 25%. Anaerobic Firmicutes dominated relative abundances at most sites but anaerobic thermophilic members of the Acidobacteria and Proteobacteria were also common. Nonetheless, compositions at the amplicon sequence variant (ASV) level varied among the sites with no convergence resulting from heating or CO addition as indicated by beta diversity analyses. The distinctions among thermophilic communities paralleled patterns observed for unheated "time zero" mesothermal soils and sediments. Occupancy analyses showed that the number of ASVs occupying each of n sites decreased unimodally with increasing n; no ASV occupied all 14 sites and only one each occupied 11 and 12 sites, while 69.3% of 1873 ASVs occupied just one site. Nonetheless, considerations of distances among the sites occupied by individual ASVs along with details of their distributions indicated that taxa were not dispersal limited but rather were constrained by environmental selection. This conclusion was supported by βMNTD and βNTI analyses, which showed dispersal limitation was only a minor contributor to taxon distributions.

Site and species contribution to β-diversity in terrestrial mammal communities: Evidence from multiple Neotropical forest sites

In a scenario where escalating human activities lead to several environmental changes and, consequently, affect mammal abundance and distribution, β-diversity may increase due to differences among sites. Using the ecological uniqueness approach, we analyzed β-diversity patterns of ground-dwelling mammal communities recorded through comprehensive camera trap monitoring within eight tropical forests protected areas in Mesoamerica and South America under variable landscape contexts. We aimed to investigate whether the contribution of single sites (LCBD) and single species (SCBD) to overall β-diversity could be explained by community metrics and environmental variables, and by species metrics and biological traits, respectively. Total β-diversity was also partitioned into species replacement and richness difference. We related LCBD to species richness, total relative abundance, functional indices, and environmental variables (tree basal area, protected area size, NDVI, and precipitation seasonality), and SCBD to species naïve occupancy, relative abundance, and morphoecological traits via beta regression. Our findings showed that LCBD was primarily explained by variation in species richness, rather than relative abundance and functional metrics. Protected area size and tree basal area were also important in explaining variation in LCBD. SCBD was strongly related to naïve occupancy and relative abundance, but not to biological traits, such as body mass, trophic energy level, activity cycle, and taxonomic category. Local β-diversity was a result of species replacements and to a lesser extent differences in species richness. Our approach was useful in examining and comparing the ecological uniqueness among different sites, revealing the regional scale current status of mammal diversity. High LCBD values comprised sites embedded within smaller habitat extents, hosting lower tree basal areas, and harboring low species richness. SCBD showed that relatively ubiquitous species that occur at variable abundances across sites contributed most to β-diversity.

Divergent assembly processes? A comparison of the plant and soil microbiome with plant communities in a glacier forefield

Community assembly is a result of dispersal, abiotic and biotic characteristics of the habitat as well as stochasticity. A direct comparison between the assembly of microbial and 'macrobial' organisms is hampered by the sampling of these communities in different studies, at different sites or on different scales. In a glacier forefield in the Austrian Alps, we recorded the soil and plant microbiome (bacteria and fungi) and plants that occurred in the same landscape and in close proximity in the same plots. We tested five predictions deduced from assembly processes and revealed deviating patterns of assembly in these community types. In short, microbes appeared to be less dispersal limited than plants and soil microbes, and plants strongly responded to abiotic factors whereas the leaf microbiome was plant species specific and well buffered from environmental conditions. The observed differences in community assembly processes may be attributed to the organisms' dispersal abilities, the exposure of the habitats to airborne propagules and habitat characteristics. The finding that assembly is conditional to the characteristics of the organisms, the habitat and the spatial scale under consideration is thus central for our understanding about the establishment and the maintenance of biodiversity.

High salinity gradients and intermediate spatial scales shaped similar biogeographical and co-occurrence patterns of microeukaryotes in a tropical freshwater-saltwater ecosystem

Microeukaryotes play key ecological roles in the microbial web of aquatic ecosystems. However, large knowledge gaps urgently need to be filled regarding the biogeography with associated shaping mechanisms and co-occurrence patterns of microeukaryotes under freshwater-saltwater gradients, especially true in tropical regions. Here, we investigated microeukaryotes of six mixed freshwater-saltwater regions in the Pearl River Estuary and surrounding coasts in southern China, with salinity ranging 0.1-32.0% and distances spanned up to 500 km, using molecular ecological methods. Results indicate that the biogeography of abundant and rare microeukaryotic communities was similar, both their co-occurrence patterns and biogeographical patterns were driven by deterministic and stochastic processes. The environmental factors with higher selective pressure than dispersal limitation meant that the role of deterministic process in structuring communities was more significant than that of stochastic process, and salinity played important role in structuring both microeukaryotic communities and networks. The abundant communities had stronger influence on entire microeukaryotic communities and seemed to be more sensitive to environmental changes than their rare counterparts, while rare ones had stronger interspecific relationships. Finally, the geographic scale and environmental gradients of study regions should firstly be clarified in future research on the ecological processes of microeukaryotes before conclusions are drawn.

Stochastic processes shape microeukaryotic community assembly in a subtropical river across wet and dry seasons

BACKGROUND

The deep mechanisms (deterministic and/or stochastic processes) underlying community assembly are a central challenge in microbial ecology. However, the relative importance of these processes in shaping riverine microeukaryotic biogeography is still poorly understood. Here, we compared the spatiotemporal and biogeographical patterns of microeukaryotic community using high-throughput sequencing of 18S rRNA gene and multivariate statistical analyses from a subtropical river during wet and dry seasons.

RESULTS

Our results provide the first description of biogeographical patterns of microeukaryotic communities in the Tingjiang River, the largest river in the west of Fujian province, southeastern China. The results showed that microeukaryotes from both wet and dry seasons exhibited contrasting community compositions, which might be owing to planktonic microeukaryotes having seasonal succession patterns. Further, all components of the microeukaryotic communities (including total, dominant, always rare, and conditionally rare taxa) exhibited a significant distance-decay pattern in both seasons, and these communities had a stronger distance-decay relationship during the dry season, especially for the conditionally rare taxa. Although several variables had a significant influence on the microeukaryotic communities, the environmental and spatial factors showed minor roles in shaping the communities. Importantly, these microeukaryotic communities were strongly driven by stochastic processes, with 89.9%, 88.5%, and 89.6% of the community variation explained by neutral community model during wet, dry, and both seasons, respectively. The neutral community model also explained a large fraction of the community variation across different taxonomic groups and levels. Additionally, the microeukaryotic taxa, which were above and below the neutral prediction, were ecologically and taxonomically distinct groups, which might be interactively structured by deterministic and stochastic processes.

CONCLUSIONS

This study demonstrated that stochastic processes are sufficient in shaping substantial variation in river microeukaryotic metacommunity across different hydrographic regimes, thereby providing a better understanding of spatiotemporal patterns, processes, and mechanisms of microeukaryotic community in waters.

Archaeal biogeography and interactions with microbial community across complex subtropical coastal waters

Marine Archaea are crucial in biogeochemical cycles, but their horizontal spatial variability, assembly processes, and microbial associations across complex coastal waters still lack characterizations at high coverage. Using a dense sampling strategy, we investigated horizontal variability in total archaeal, Thaumarchaeota Marine Group (MG) I, and Euryarchaeota MGII communities and associations of MGI/MGII with other microbes in surface waters with contrasting environmental characteristics across ~200 km by 16S rRNA gene amplicon sequencing. Total archaeal communities were extremely dominated by MGI and/or MGII (98.9% in average relative abundance). Niche partitioning between MGI and MGII or within each group was found across multiple environmental gradients. "Selection" was more important than "dispersal limitation" in governing biogeographic patterns of total archaeal, MGI, and MGII communities, and basic abiotic parameters (such as salinity) and inorganic/organic resources as a whole could be the main driver of "selection". While "homogenizing dispersal" also considerably governed their biogeography. MGI-Nitrospira assemblages were speculatively responsible for complete nitrification. MGI taxa commonly had negative correlations with members of Synechococcus but positive correlations with members of eukaryotic phytoplankton, suggesting that competition or synergy between MGI and phytoplankton depends on specific MGI-phytoplankton assemblages. MGII taxa showed common associations with presumed (photo)heterotrophs including members of SAR11, SAR86, SAR406, and Candidatus Actinomarina. This study sheds light on ecological processes and drivers shaping archaeal biogeography and many strong MGI/MGII-bacterial associations across complex subtropical coastal waters. Future efforts should be made on seasonality of archaeal biogeography and biological, environmental, or ecological mechanisms underlying these statistical microbial associations.

An ordination approach to explore similarities among communities

Analysis of similarities among communities can help to decipher the biogeographical, evolutionary, and ecological factors that drive local diversity. Recent indices of similarity among communities incorporate not only information on species presence and abundance but also information on how similar species are in their traits and how closely related they are in terms of taxonomy or phylogeny. Towards this aim, trait-based, taxonomic or phylogenetic similarities among species have been defined and bounded between 0 (species are maximally distinct) and 1 (species are similar). A required property for an index of similarity between two communities is that it must provide minimum similarity (0) where communities have maximally distinct species, as well as maximum similarity (1) where communities are equivalent in their trait, taxonomic or phylogenetic compositions. Here, I developed a new ordination methodology that conforms to the requirement: double similarity principal component analysis (DSPCA). DSPCA summarizes multidimensional trait-based, taxonomic or phylogenetic similarities among communities into orthogonal axes. The species that drive each similarity pattern can be identified together with their traits or with their taxonomic or phylogenetic positions. I applied this methodology to theoretical examples and to empirical data sets on bird and bat communities to illustrate key properties of DSPCA. I compared the results obtained with DSPCA with those provided by related approaches. Theoretical and empirical case studies highlight the following additional properties of DSPCA: (i) axes are orthogonal and identify independent (dis)similarity patterns between communities; (ii) the more functionally, taxonomically or phylogenetically similar communities are, the closer they are on an axis; (iii) the coordinate of a species on an axis expresses how representative the species is of the pattern identified by the axis; and (iv) a species is representative of x communities if the functional, taxonomic or phylogenetic characteristics of this species are very common within each of these x communities. DSPCA is an efficient approach to visualize functional, taxonomic and phylogenetic similarities between communities. It is also a useful alternative to recent methods dedicated to phylogenetic diversity patterns. It will be an asset for all studies that aim to compare functional, taxonomic, genetic and phylogenetic diversity.

Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes

Unraveling the relative importance of ecological processes regulating microbial community structure is a central goal in microbial ecology. Here, we used high-throughput sequencing to examine the relative contribution of selective and neutral processes in the assembly of abundant and rare subcommunities from three subtropical bays of China. We found that abundant and rare bacterial taxa were distinctly different in diversity, despite the similar biogeographic patterns and strong distance-decay relationships, but the dispersal of rare bacterial taxa was more limited than that of abundant taxa. Furthermore, the environmental (selective processes) and spatial (neutral processes) factors seemed to govern the assembly and biogeography of abundant and rare bacterial subcommunities, although both factors explained only a small fraction of variation within the rare subcommunity. More importantly, variation partitioning (based on adjusted R2 in redundancy analysis) showed that spatial factors exhibited a slightly greater influence on both abundant and rare subcommunities compared to environmental selection; however, the abundant subcommunity had a much stronger response to spatial factors (17.3% of pure variance was explained) than that shown by the rare bacteria (3.5%). These results demonstrate that environmental selection and neutral processes explained the similar biogeographic patterns of abundant and rare subcommunities, but a large proportion of unexplained variation in the rare taxa (91.1%) implies that more complex assembly mechanisms may exist to shape the rare bacterial assemblages in the three subtropical bays.

Stochastic Community Assembly: Does It Matter in Microbial Ecology?

Understanding the mechanisms controlling community diversity, functions, succession, and biogeography is a central, but poorly understood, topic in ecology, particularly in microbial ecology. Although stochastic processes are believed to play nonnegligible roles in shaping community structure, their importance relative to deterministic processes is hotly debated. The importance of ecological stochasticity in shaping microbial community structure is far less appreciated. Some of the main reasons for such heavy debates are the difficulty in defining stochasticity and the diverse methods used for delineating stochasticity. Here, we provide a critical review and synthesis of data from the most recent studies on stochastic community assembly in microbial ecology. We then describe both stochastic and deterministic components embedded in various ecological processes, including selection, dispersal, diversification, and drift. We also describe different approaches for inferring stochasticity from observational diversity patterns and highlight experimental approaches for delineating ecological stochasticity in microbial communities. In addition, we highlight research challenges, gaps, and future directions for microbial community assembly research.

Stochastic Community Assembly: Does It Matter in Microbial Ecology?

Understanding the mechanisms controlling community diversity, functions, succession, and biogeography is a central, but poorly understood, topic in ecology, particularly in microbial ecology. Although stochastic processes are believed to play nonnegligible roles in shaping community structure, their importance relative to deterministic processes is hotly debated. The importance of ecological stochasticity in shaping microbial community structure is far less appreciated. Some of the main reasons for such heavy debates are the difficulty in defining stochasticity and the diverse methods used for delineating stochasticity. Here, we provide a critical review and synthesis of data from the most recent studies on stochastic community assembly in microbial ecology. We then describe both stochastic and deterministic components embedded in various ecological processes, including selection, dispersal, diversification, and drift. We also describe different approaches for inferring stochasticity from observational diversity patterns and highlight experimental approaches for delineating ecological stochasticity in microbial communities. In addition, we highlight research challenges, gaps, and future directions for microbial community assembly research.

Dispersal-Based Microbial Community Assembly Decreases Biogeochemical Function

Ecological mechanisms influence relationships among microbial communities, which in turn impact biogeochemistry. In particular, microbial communities are assembled by deterministic (e.g., selection) and stochastic (e.g., dispersal) processes, and the relative balance of these two process types is hypothesized to alter the influence of microbial communities over biogeochemical function. We used an ecological simulation model to evaluate this hypothesis, defining biogeochemical function generically to represent any biogeochemical reaction of interest. We assembled receiving communities under different levels of dispersal from a source community that was assembled purely by selection. The dispersal scenarios ranged from no dispersal (i.e., selection-only) to dispersal rates high enough to overwhelm selection (i.e., homogenizing dispersal). We used an aggregate measure of community fitness to infer a given community’s biogeochemical function relative to other communities. We also used ecological null models to further link the relative influence of deterministic assembly to function. We found that increasing rates of dispersal decrease biogeochemical function by increasing the proportion of maladapted taxa in a local community. Niche breadth was also a key determinant of biogeochemical function, suggesting a tradeoff between the function of generalist and specialist species. Finally, we show that microbial assembly processes exert greater influence over biogeochemical function when there is variation in the relative contributions of dispersal and selection among communities. Taken together, our results highlight the influence of spatial processes on biogeochemical function and indicate the need to account for such effects in models that aim to predict biogeochemical function under future environmental scenarios.

Effects of latitude and depth on the beta diversity of New Zealand fish communities

Marine ecosystems are difficult to sample quantitatively at increasing depth. Hence, few studies attempt to measure patterns of beta diversity for ecological communities in the deep sea. Here we (i) present and quantify large-scale gradients in fish community structure along depth and latitude gradients of the New Zealand EEZ, (ii) obtain rigorous quantitative estimates of these depth (50-1200 m) and latitudinal effects (29.15-50.91°S) and their interaction, and (iii) explicitly model how latitudinal beta diversity of fishes varies with depth. The sampling design was highly structured, replicated and stratified for latitude and depth, using data obtained from 345 standardised baited remote underwater stereo-video deployments. Results showed that gradients in fish community structure along depth and latitude were strong and interactive in New Zealand waters; latitudinal variation in fish communities progressively decreased with depth following an exponential decay (r ² = 0.96), revealing increasingly similar fish communities with increasing depth. In contrast, variation in fish community structure along the depth gradient was of a similar magnitude across all of the latitudes investigated here. We conclude that an exponential decay in beta diversity vs depth exists for fish communities present in areas shallower than the New Zealand upper continental slope.

Revisiting the role of dissimilarity of host communities in driving dissimilarity of ectoparasite assemblages: non-linear vs linear approach

We revisited the role of dissimilarity of host assemblages in shaping dissimilarity of flea assemblages using a non-linear approach. Generalized dissimilarity models (GDMs) were applied using data from regional surveys of fleas parasitic on small mammals in four biogeographical realms. We compared (1) model fit, (2) the relative effects of host compositional and phylogenetic turnover and geographic distance on flea compositional and phylogenetic turnover, and (3) the rate of flea turnover along gradients of host turnover and geographic distance with those from earlier application of a linear approach. GDMs outperformed linear models in explaining variation in flea species turnover and host dissimilarity was the best predictor of flea dissimilarity, irrespective of scale. The shape of the relationships between flea compositional turnovers along host compositional turnover was similar in all realms, whereas turnover along geographic distance differed among realms. In contrast, the rate of flea phylogenetic turnover along gradients of host phylogenetic turnover differed among realms, whereas flea phylogenetic turnover did not depend on geographic distance in any realm. We demonstrated that a non-linear approach (a) explained spatial variation in parasite community composition better than and (b) revealed patterns that were obscured by earlier linear analyses.

Making sense of metacommunities: dispelling the mythology of a metacommunity typology

Metacommunity ecology has rapidly become a dominant framework through which ecologists understand the natural world. Unfortunately, persistent misunderstandings regarding metacommunity theory and the methods for evaluating hypotheses based on the theory are common in the ecological literature. Since its beginnings, four major paradigms-species sorting, mass effects, neutrality, and patch dynamics-have been associated with metacommunity ecology. The Big 4 have been misconstrued to represent the complete set of metacommunity dynamics. As a result, many investigators attempt to evaluate community assembly processes as strictly belonging to one of the Big 4 types, rather than embracing the full scope of metacommunity theory. The Big 4 were never intended to represent the entire spectrum of metacommunity dynamics and were rather examples of historical paradigms that fit within the new framework. We argue that perpetuation of the Big 4 typology hurts community ecology and we encourage researchers to embrace the full inference space of metacommunity theory. A related, but distinct issue is that the technique of variation partitioning is often used to evaluate the dynamics of metacommunities. This methodology has produced its own set of misunderstandings, some of which are directly a product of the Big 4 typology and others which are simply the product of poor study design or statistical artefacts. However, variation partitioning is a potentially powerful technique when used appropriately and we identify several strategies for successful utilization of variation partitioning.

Groundwater-surface water mixing shifts ecological assembly processes and stimulates organic carbon turnover

Environmental transitions often result in resource mixtures that overcome limitations to microbial metabolism, resulting in biogeochemical hotspots and moments. Riverine systems, where groundwater mixes with surface water (the hyporheic zone), are spatially complex and temporally dynamic, making development of predictive models challenging. Spatial and temporal variations in hyporheic zone microbial communities are a key, but understudied, component of riverine biogeochemical function. Here, to investigate the coupling among groundwater-surface water mixing, microbial communities and biogeochemistry, we apply ecological theory, aqueous biogeochemistry, DNA sequencing and ultra-high-resolution organic carbon profiling to field samples collected across times and locations representing a broad range of mixing conditions. Our results indicate that groundwater-surface water mixing in the hyporheic zone stimulates heterotrophic respiration, alters organic carbon composition, causes ecological processes to shift from stochastic to deterministic and is associated with elevated abundances of microbial taxa that may degrade a broad suite of organic compounds.

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.

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.

The Biogeography of Adaptive Radiations and the Geographic Overlap of Sister Species

The biogeography of speciation and what can be learned about the past mode of speciation from current biogeography of sister species are recurrent problems in evolution. We used a trait- and individual-based, eco-evolutionary model to simulate adaptive radiations and recorded the geographical overlap of species during and after evolutionary branching (speciation). We compared the spatial overlap among sister species in the fully saturated community with the overlap at the speciation event. The mean geographic overlap at speciation varied continuously from complete (sympatry) to none (allopatry), depending on local and regional environmental heterogeneity and the rate of dispersal. The distribution of overlap was, however, in some cases considerably bimodal. This tendency was most expressed at large values of regional heterogeneity, corresponding to sharp environmental contrasts. The mean geographic overlap also varied during the course of a radiation, sometimes with a consistent negative trend over time. The speciations that resulted in currently observable end community sister species were therefore not an unbiased sample of all speciations throughout the radiation. Postspeciation range shifts (causing increased overlap) occurred most frequently when dispersal was high or when local habitat heterogeneity was low. Our results help us understand how the patterns of geographic mode of speciation emerge. We also show the difficulty in inferring the geographical speciation mode from phylogenies and the biogeography of extant species.

Integrating Taxonomic, Functional and Phylogenetic Beta Diversities: Interactive Effects with the Biome and Land Use across Taxa

The spatial distribution of species, functional traits and phylogenetic relationships at both the regional and local scales provide complementary approaches to study patterns of biodiversity and help to untangle the mechanisms driving community assembly. Few studies have simultaneously considered the taxonomic (TBD), functional (FBD) and phylogenetic (PBD) facets of beta diversity. Here we analyze the associations between TBD, FBD, and PBD with the biome (representing different regional species pools) and land use, and investigate whether TBD, FBD and PBD were correlated. In the study design we considered two widely used indicator taxa (birds and ants) from two contrasting biomes (subtropical forest and grassland) and land uses (tree plantations and cropfields) in the southern Neotropics. Non-metric multidimensional scaling showed that taxonomic, functional and phylogenetic distances were associated to biome and land use; study sites grouped into four groups on the bi-dimensional space (cropfields in forest and grassland, and tree plantations in forest and grassland), and that was consistent across beta diversity facets and taxa. Mantel and PERMANOVA tests showed that TBD, FBD and PBD were positively correlated for both bird and ant assemblages; in general, partial correlations were also significant. Some of the functional traits considered here were conserved along phylogeny. Our results will contribute to the development of sound land use planning and beta diversity conservation.

Estimating and mapping ecological processes influencing microbial community assembly

Ecological community assembly is governed by a combination of (i) selection resulting from among-taxa differences in performance; (ii) dispersal resulting from organismal movement; and (iii) ecological drift resulting from stochastic changes in population sizes. The relative importance and nature of these processes can vary across environments. Selection can be homogeneous or variable, and while dispersal is a rate, we conceptualize extreme dispersal rates as two categories; dispersal limitation results from limited exchange of organisms among communities, and homogenizing dispersal results from high levels of organism exchange. To estimate the influence and spatial variation of each process we extend a recently developed statistical framework, use a simulation model to evaluate the accuracy of the extended framework, and use the framework to examine subsurface microbial communities over two geologic formations. For each subsurface community we estimate the degree to which it is influenced by homogeneous selection, variable selection, dispersal limitation, and homogenizing dispersal. Our analyses revealed that the relative influences of these ecological processes vary substantially across communities even within a geologic formation. We further identify environmental and spatial features associated with each ecological process, which allowed mapping of spatial variation in ecological-process-influences. The resulting maps provide a new lens through which ecological systems can be understood; in the subsurface system investigated here they revealed that the influence of variable selection was associated with the rate at which redox conditions change with subsurface depth.

On the processes generating latitudinal richness gradients: identifying diagnostic patterns and predictions

We use a simulation model to examine four of the most common hypotheses for the latitudinal richness gradient and identify patterns that might be diagnostic of those four hypotheses. The hypotheses examined include (1) tropical niche conservatism, or the idea that the tropics are more diverse because a tropical clade origin has allowed more time for diversification in the tropics and has resulted in few species adapted to extra-tropical climates. (2) The ecological limits hypothesis suggests that species richness is limited by the amount of biologically available energy in a region. (3) The speciation rates hypothesis suggests that the latitudinal gradient arises from a gradient in speciation rates. (4) Finally, the tropical stability hypothesis argues that climatic fluctuations and glacial cycles in extratropical regions have led to greater extinction rates and less opportunity for specialization relative to the tropics. We found that tropical niche conservatism can be distinguished from the other three scenarios by phylogenies which are more balanced than expected, no relationship between mean root distance (MRD) and richness across regions, and a homogeneous rate of speciation across clades and through time. The energy gradient, speciation gradient, and disturbance gradient scenarios all produced phylogenies which were more imbalanced than expected, showed a negative relationship between MRD and richness, and diversity-dependence of speciation rate estimates through time. We found that the relationship between speciation rates and latitude could distinguish among these three scenarios, with no relation expected under the ecological limits hypothesis, a negative relationship expected under the speciation rates hypothesis, and a positive relationship expected under the tropical stability hypothesis. We emphasize the importance of considering multiple hypotheses and focusing on diagnostic predictions instead of predictions that are consistent with multiple hypotheses.

Evaluation of an integrated framework for biodiversity with a new metric for functional dispersion

Growing interest in understanding ecological patterns from phylogenetic and functional perspectives has driven the development of metrics that capture variation in evolutionary histories or ecological functions of species. Recently, an integrated framework based on Hill numbers was developed that measures three dimensions of biodiversity based on abundance, phylogeny and function of species. This framework is highly flexible, allowing comparison of those diversity dimensions, including different aspects of a single dimension and their integration into a single measure. The behavior of those metrics with regard to variation in data structure has not been explored in detail, yet is critical for ensuring an appropriate match between the concept and its measurement. We evaluated how each metric responds to particular data structures and developed a new metric for functional biodiversity. The phylogenetic metric is sensitive to variation in the topology of phylogenetic trees, including variation in the relative lengths of basal, internal and terminal branches. In contrast, the functional metric exhibited multiple shortcomings: (1) species that are functionally redundant contribute nothing to functional diversity and (2) a single highly distinct species causes functional diversity to approach the minimum possible value. We introduced an alternative, improved metric based on functional dispersion that solves both of these problems. In addition, the new metric exhibited more desirable behavior when based on multiple traits.

Distance-based functional diversity measures and their decomposition: a framework based on Hill numbers

Hill numbers (or the “effective number of species”) are increasingly used to characterize species diversity of an assemblage. This work extends Hill numbers to incorporate species pairwise functional distances calculated from species traits. We derive a parametric class of functional Hill numbers, which quantify “the effective number of equally abundant and (functionally) equally distinct species” in an assemblage. We also propose a class of mean functional diversity (per species), which quantifies the effective sum of functional distances between a fixed species to all other species. The product of the functional Hill number and the mean functional diversity thus quantifies the (total) functional diversity, i.e., the effective total distance between species of the assemblage. The three measures (functional Hill numbers, mean functional diversity and total functional diversity) quantify different aspects of species trait space, and all are based on species abundance and species pairwise functional distances. When all species are equally distinct, our functional Hill numbers reduce to ordinary Hill numbers. When species abundances are not considered or species are equally abundant, our total functional diversity reduces to the sum of all pairwise distances between species of an assemblage. The functional Hill numbers and the mean functional diversity both satisfy a replication principle, implying the total functional diversity satisfies a quadratic replication principle. When there are multiple assemblages defined by the investigator, each of the three measures of the pooled assemblage (gamma) can be multiplicatively decomposed into alpha and beta components, and the two components are independent. The resulting beta component measures pure functional differentiation among assemblages and can be further transformed to obtain several classes of normalized functional similarity (or differentiation) measures, including N-assemblage functional generalizations of the classic Jaccard, Sørensen, Horn and Morisita-Horn similarity indices. The proposed measures are applied to artificial and real data for illustration.

Taxonomic identification of algae (morphological and molecular): species concepts, methodologies, and their implications for ecological bioassessment

Algal taxonomy is a key discipline in phycology and is critical for algal genetics, physiology, ecology, applied phycology, and particularly bioassessment. Taxonomic identification is the most common analysis and hypothesis-testing endeavor in science. Errors of identification are often related to the inherent problem of small organisms with morphologies that are difficult to distinguish without research-grade microscopes and taxonomic expertise in phycology. Proposed molecular approaches for taxonomic identification from environmental samples promise rapid, potentially inexpensive, and more thorough culture-independent identification of all algal species present in a sample of interest. Molecular identification has been used in biodiversity and conservation, but it also has great potential for applications in bioassessment. Comparisons of morphological and molecular identification of benthic algal communities are improved by the identification of more taxa; however, automated identification technology does not allow for the simultaneous analysis of thousands of samples. Currently, morphological identification is used to verify molecular taxonomic identities, but with the increased number of taxa verified in algal gene libraries, molecular identification will become a universal tool in biological studies. Thus, in this report, successful application of molecular techniques related to algal bioassessment is discussed.

Ecosystem-wide morphological structure of leaf-litter ant communities along a tropical latitudinal gradient

General principles that shape community structure can be described based on a functional trait approach grounded on predictive models; increased attention has been paid to factors accounting for the functional diversity of species assemblages and its association with species richness along environmental gradients. We analyze here the interaction between leaf-litter ant species richness, the local communities' morphological structure and fundamental niche within the context of a northeast-southeast latitudinal gradient in one of the world's most species-rich ecosystems, the Atlantic Forest, representing 2,700 km of tropical rainforest along almost 20^o of latitude in eastern Brazil. Our results are consistent with an ecosystem-wide pattern in communities' structure, with relatively high species turnover but functionally analogous leaf-litter ant communities' organization. Our results suggest directional shifts in the morphological space along the environmental gradient from overdispersed to aggregated (from North to South), suggesting that primary productivity and environmental heterogeneity (altitude, temperature and precipitation in the case) determine the distribution of traits and regulate the assembly rules, shaping local leaf-litter ant communities. Contrary to the expected and most common pattern along latitudinal gradients, the Atlantic Forest leaf litter ant communities show an inverse pattern in richness, that is, richer communities in higher than in lower latitudes. The morphological specialization of communities showed more morphologically distinct communities at low latitudes and species redundancy at high latitudes. We claim that an inverse latitudinal gradient in primary productivity and environmental heterogeneity across the Atlantic forest may affect morphological diversity and species richness, enhancing species coexistence mechanisms, and producing thus the observed patterns. We suggest that a functional framework based on flexible enough traits should be pursued to allow comparisons at local, regional and global levels.

Coarse root spatial distribution determined using a ground-penetrating radar technique in a subtropical evergreen broad-leaved forest, China

Coarse roots play a critical role in forest ecosystems and both abiotic and biotic factors affect their spatial distribution. To some extent, coarse root density may reflect the quantity of root biomass and biotic competition in forests. However, using traditional methods (e.g., excavation) to study coarse roots is challenging, because those methods are time-consuming and laborious. Furthermore, these destructive methods cannot be repeated in the same forests. Therefore, the discovery of non-destructive methods for root studies will be very significant. In this study, we used a ground-penetrating radar technique to detect the coarse root density of three habitats (ridge, slope and valley) and the dominant tree species (Castanopsis eyrei and Schima superba) in a subtropical forest. We found that (i) the mean of coarse root density for these three habitats was 88.04 roots m(-2), with roots being mainly distributed at depths of 0-40 cm. Coarse root densities were lower in deeper soils and in areas far from the trunk. (ii) Coarse root densities differed significantly among the three habitats studied here with slope habitat having the lowest coarse root density. Compared with S. superba, C. eyrei had more roots distributed in deeper soils. Furthermore, coarse roots with a diameter >3 cm occurred more frequently in the valleys, compared with root densities in ridge and slope habitats, and most coarse roots occurred at soil depths of 20-40 cm. (iii) The coarse root density correlated negatively with tree species richness at soil depths of 40-60 cm. The abundances of the dominant species, such as C. eyrei, Cyclobalanopsis glauca, Pinus massoniana, had significant impacts on coarse root density. (iv) The soil depth of 0-40 cm was the "basic distribution layer" for coarse roots since the majority of coarse roots were found in this soil layer with an average root density of 84.18 roots m(-2), which had no significant linear relationships with topography, tree species richness, rarefied tree species richness and tree density. Significant relationships between coarse root density and these factors were found at the soil depth of 40-60 cm, which was the "potential distribution layer" for coarse root distribution.

Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes

Increasing evidence has emerged for non-random spatial distributions of microbes, but knowledge of the processes that cause variation in microbial assemblage among ecosystems is lacking. For instance, some studies showed that deterministic processes such as habitat specialization are important, while other studies hold that bacterial communities are assembled by stochastic forces. Here we examine the relative influence of deterministic and stochastic processes for bacterial communities from subsurface environments, stream biofilm, lake water, lake sediment and soil using pyrosequencing of the 16S ribosomal RNA gene. We show that there is a general pattern in phylogenetic signal in species ecological niches across recent evolutionary time for all studied habitats, enabling us to infer the influences of community assembly processes from patterns of phylogenetic turnover in community composition. The phylogenetic dissimilarities among-habitat types were significantly higher than within them, and the communities were clustered according to their original habitat types. For communities within-habitat types, the highest phylogenetic turnover rate through space was observed in subsurface environments, followed by stream biofilm on mountainsides, whereas the sediment assemblages across regional scales showed the lowest turnover rate. Quantifying phylogenetic turnover as the deviation from a null expectation suggested that measured environmental variables imposed strong selection on bacterial communities for nearly all sample groups. For three sample groups, spatial distance reflected unmeasured environmental variables that impose selection, as opposed to spatial isolation. Such characterization of spatial and environmental variables proved essential for proper interpretation of partial Mantel results based on observed beta diversity metrics. In summary, our results clearly indicate a dominant role of deterministic processes on bacterial assemblages and highlight that bacteria show strong habitat associations that have likely emerged through evolutionary adaptation.

Quantifying community assembly processes and identifying features that impose them

Spatial turnover in the composition of biological communities is governed by (ecological) Drift, Selection and Dispersal. Commonly applied statistical tools cannot quantitatively estimate these processes, nor identify abiotic features that impose these processes. For interrogation of subsurface microbial communities distributed across two geologically distinct formations of the unconfined aquifer underlying the Hanford Site in southeastern Washington State, we developed an analytical framework that advances ecological understanding in two primary ways. First, we quantitatively estimate influences of Drift, Selection and Dispersal. Second, ecological patterns are used to characterize measured and unmeasured abiotic variables that impose Selection or that result in low levels of Dispersal. We find that (i) Drift alone consistently governs ∼25% of spatial turnover in community composition; (ii) in deeper, finer-grained sediments, Selection is strong (governing ∼60% of turnover), being imposed by an unmeasured but spatially structured environmental variable; (iii) in shallower, coarser-grained sediments, Selection is weaker (governing ∼30% of turnover), being imposed by vertically and horizontally structured hydrological factors;(iv) low levels of Dispersal can govern nearly 30% of turnover and be caused primarily by spatial isolation resulting from limited exchange between finer and coarser-grain sediments; and (v) highly permeable sediments are associated with high levels of Dispersal that homogenize community composition and govern over 20% of turnover. We further show that our framework provides inferences that cannot be achieved using preexisting approaches, and suggest that their broad application will facilitate a unified understanding of microbial communities.

Beta diversity of demersal fish assemblages in the North-Eastern Pacific: interactions of latitude and depth

Knowledge of broad-scale global patterns in beta diversity (i.e., variation or turnover in identities of species) for marine systems is in its infancy. We analysed the beta diversity of groundfish communities along the North American Pacific coast, from trawl data spanning 32.57°N to 48.52°N and 51 m to 1200 m depth. Analyses were based on both the Jaccard measure and the probabilistic Raup-Crick measure, which accounts for variation in alpha diversity. Overall, beta diversity decreased with depth, and this effect was strongest at lower latitudes. Superimposed on this trend were peaks in beta diversity at around 400–600 m and also around 1000–1200 m, which may indicate high turnover around the edges of the oxygen minimum zone. Beta diversity was also observed to decrease with latitude, but this effect was only observed in shallower waters (<200 m); latitudinal turnover began to disappear at depths >800 m. At shallower depths (<200 m), peaks in latitudinal turnover were observed at ∼43°N, 39°N, 35°N and 31°N, which corresponded well with several classically observed oceanographic boundaries. Turnover with depth was stronger than latitudinal turnover, and is likely to reflect strong environmental filtering over relatively short distances. Patterns in beta diversity, including latitude-by-depth interactions, should be integrated with other biodiversity measures in ecosystem-based management and conservation of groundfish communities.

Beta-diversity in temperate and tropical forests reflects dissimilar mechanisms of community assembly

Site-to-site variation in species composition (β-diversity) generally increases from low- to high-diversity regions. Although biogeographical differences in community assembly mechanisms may explain this pattern, random sampling effects can create this pattern through differences in regional species pools. Here, we compared assembly mechanisms between spatially extensive networks of temperate and tropical forest plots with highly divergent species pools (46 vs. 607 species). After controlling for sampling effects, β-diversity of woody plants was similar and higher than expected by chance in both forests, reflecting strong intraspecific aggregation. However, different mechanisms appeared to explain aggregation in the two forests. In the temperate forest, aggregation reflected stronger environmental correlations, suggesting an important role for species-sorting (e.g. environmental filtering) processes, whereas in the tropics, aggregation reflected stronger spatial correlations, more likely reflecting dispersal limitation. We suggest that biogeographical differences in the relative importance of different community assembly mechanisms contribute to these striking gradients in global biodiversity.

Low functional β-diversity despite high taxonomic β-diversity among tropical estuarine fish communities

The concept of β-diversity, defined as dissimilarity among communities, has been widely used to investigate biodiversity patterns and community assembly rules. However, in ecosystems with high taxonomic β-diversity, due to marked environmental gradients, the level of functional β-diversity among communities is largely overlooked while it may reveal processes shaping community structure. Here, decomposing biodiversity indices into α (local) and γ (regional) components, we estimated taxonomic and functional β-diversity among tropical estuarine fish communities, through space and time. We found extremely low functional β-diversity values among fish communities (<1.5%) despite high dissimilarity in species composition and species dominance. Additionally, in contrast to the high α and γ taxonomic diversities, α and γ functional diversities were very close to the minimal value. These patterns were caused by two dominant functional groups which maintained a similar functional structure over space and time, despite the strong dissimilarity in taxonomic structure along environmental gradients. Our findings suggest that taxonomic and functional β-diversity deserve to be quantified simultaneously since these two facets can show contrasting patterns and the differences can in turn shed light on community assembly rules.

Dispersal, environmental niches and oceanic-scale turnover in deep-sea bivalves

Patterns of beta-diversity or distance decay at oceanic scales are completely unknown for deep-sea communities. Even when appropriate data exist, methodological problems have made it difficult to discern the relative roles of environmental filtering and dispersal limitation for generating faunal turnover patterns. Here, we combine a spatially extensive dataset on deep-sea bivalves with a model incorporating ecological dynamics and shared evolutionary history to quantify the effects of environmental filtering and dispersal limitation. Both the model and empirical data are used to relate functional, taxonomic and phylogenetic similarity between communities to environmental and spatial distances separating them for 270 sites across the Atlantic Ocean. This study represents the first ocean-wide analysis examining distance decay as a function of a broad suite of explanatory variables. We find that both strong environmental filtering and dispersal limitation drive turnover in taxonomic, functional and phylogenetic composition in deep-sea bivalves, explaining 26 per cent, 34 per cent and 9 per cent of the variation, respectively. This contrasts with previous suggestions that dispersal is not limiting in broad-scale biogeographic and biodiversity patterning in marine systems. However, rates of decay in similarity with environmental distance were eightfold to 44-fold steeper than with spatial distance. Energy availability is the most influential environmental variable evaluated, accounting for 3.9 per cent, 9.4 per cent and 22.3 per cent of the variation in functional, phylogenetic and taxonomic similarity, respectively. Comparing empirical patterns with process-based theoretical predictions provided quantitative estimates of dispersal limitation and niche breadth, indicating that 95 per cent of deep-sea bivalve propagules will be able to persist in environments that deviate from their optimum by up to 2.1 g m(-2) yr(-1) and typically disperse 749 km from their natal site.