Microbial biogeography of acid mine drainage sediments at a regional scale across Southern China

Investigations of microbial biogeography in extreme environments provide unique opportunities to disentangle the roles of environment and space in microbial community assembly. Here, we reported a comprehensive microbial biogeographic survey of 90 acid mine drainage (AMD) sediment samples from 18 mining sites of various mineral types across southern China. We found that environmental selection was strong in determining the AMD habitat species pool. However, microbial alpha diversity was primarily explained by mining sites rather than environmental factors, and microbial beta diversity correlated more strongly with geographic than environmental distance at both large and small spatial scales. Particularly, the presence/absence of widespread AMD habitat generalists was only correlated with geographic distance and independent of environmental variation. These distance-decay patterns suggested that spatial processes played a more important role in determining microbial compositional variation across space; which could be explained by the reinforced impacts of dispersal limitation in less fluid, spatially structured sediment habitat with diverse pre-existing communities. In summary, our findings suggested that the deterministic assembling and spatial constraints interact to shape microbial biogeography in AMD sediments; and provided implications that spatial processes should be considered when predicting microbial dynamics in response to severe environmental change across large spatial scales.

Effects of phenotypic variation on consumer coexistence and prey community structure

A popular idea in ecology is that trait variation among individuals from the same species may promote the coexistence of competing species. However, theoretical and empirical tests of this idea have yielded inconsistent findings. We manipulated intraspecific trait diversity in a ciliate competing with a nematode for bacterial prey in experimental microcosms. We found that intraspecific trait variation inverted the original competitive hierarchy to favour the consumer with variable traits, ultimately resulting in competitive exclusion. This competitive outcome was driven by foraging traits (size, speed and directionality) that increased the ciliate's fitness ratio and niche overlap with the nematode. The interplay between consumer trait variation and competition resulted in non-additive cascading effects-mediated through prey defence traits-on prey community assembly. Our results suggest that predicting consumer competitive population dynamics and the assembly of prey communities will require understanding the complexities of trait variation within consumer species.

Climate warming promotes deterministic assembly of arbuscular mycorrhizal fungal communities

Arbuscular mycorrhizal fungi (AMF) significantly contribute to plant resource acquisition and play important roles in mediating plant interactions and soil carbon (C) dynamics. However, it remains unclear how AMF communities respond to climate change. We assessed impacts of warming and precipitation alterations (30% increase or decrease) on soil AMF communities, and examined major ecological processes shaping the AMF community assemblage in a Tibetan alpine meadow. Our results showed that warming significantly increased root biomass, and available nitrogen (N) and phosphorus (P) in soil. While precipitation alterations increased AMF abundances, they did not significantly affect the composition or diversity of AMF communities. In contrast, warming altered the composition of AMF communities and reduced their Shannon-Wiener index and Pielou's evenness. In particular, warming shifted the AMF community composition in favor of Diversisporaceae over Glomeraceae, likely through its impact on soil N and P availability. In addition, AMF communities were phylogenetically random in the unwarmed control but clustered in warming plots, implying more deterministic community assembly under climate warming. Warming enhancement of root growth, N and P availability likely reduced plant C-allocation to AMF, imposing stronger environmental filtering on AMF communities. We further proposed a conceptual framework that integrates biological and geochemical processes into a mechanistic understanding of warming and precipitation changes' effects on AMF. Taken together, these results suggest that soil AMF communities may be more sensitive to warming than expected, highlighting the need to monitor their community structure and associated functional consequences on plant communities and soil C dynamics under the future warmer climate.

The Complex Biodiversity-Ecosystem Function Relationships for the Qinghai-Tibetan Grassland Community

Despite the long history of the study of the biodiversity-ecosystem function relationship, uncertainty remains about the relationship of natural grassland ecosystems under stressful conditions. Recently, trait- and phylogenetic-based tests provide a powerful way to detect the relationship in different spaces but have seldom been applied to stressful zones on a large spatial scale. We selected Qinghai-Tibetan as the study area and collected a grassland community database involving 581 communities. We calculated biomass and species', functional, and phylogenetic diversity of each community and examined their relationships by using linear and non-linear regression models. Results showed an overall positive biodiversity-productivity relationship in species', functional and phylogenetic space. The relationship, however, was non-linear, in which biodiversity explained better the variation in community biomass when species diversity was more than a threshold, showing a weak effect of biodiversity on ecosystem function in low species diversity communities. We also found a filled triangle for the limit of the relationship between species and functional diversity, implying that functional diversity differs significantly among communities when their species diversity is low but finally converges to be a constant with increasing communities' species diversity. Our study suggests that multiple niche processes may structure the grassland communities, and their forces tend to balance in high-biodiversity communities.

Higher-order effects, continuous species interactions, and trait evolution shape microbial spatial dynamics

Persistently diverse microbial communities are one of biology’s great puzzles. Using a modeling framework that accommodates high mutation rates and a continuum of species traits, we studied microbial communities in which antagonistic interactions occur via the production of, inhibition of, and vulnerability to toxins (e.g., antibiotics). Mutation size and mobility enhanced microbial diversity and temporal persistence to extraordinarily high levels. These findings—including the discovery that the duration of the transient phase in community assembly provides a guide to equilibrial diversity—highlight the potentially critical role that antagonistic interactions play in promoting the diversity of bacterial systems. Such interactions, together with resource-driven interactions and spatial structure, may drive the enigmatic levels of biodiversity seen in microbial systems.

The assembly and maintenance of microbial diversity in natural communities, despite the abundance of toxin-based antagonistic interactions, presents major challenges for biological understanding. A common framework for investigating such antagonistic interactions involves cyclic dominance games with pairwise interactions. The incorporation of higher-order interactions in such models permits increased levels of microbial diversity, especially in communities in which antibiotic-producing, sensitive, and resistant strains coexist. However, most such models involve a small number of discrete species, assume a notion of pure cyclic dominance, and focus on low mutation rate regimes, none of which well represent the highly interlinked, quickly evolving, and continuous nature of microbial phenotypic space. Here, we present an alternative vision of spatial dynamics for microbial communities based on antagonistic interactions—one in which a large number of species interact in continuous phenotypic space, are capable of rapid mutation, and engage in both direct and higher-order interactions mediated by production of and resistance to antibiotics. Focusing on toxin production, vulnerability, and inhibition among species, we observe highly divergent patterns of diversity and spatial community dynamics. We find that species interaction constraints (rather than mobility) best predict spatiotemporal disturbance regimes, whereas community formation time, mobility, and mutation size best explain patterns of diversity. We also report an intriguing relationship among community formation time, spatial disturbance regimes, and diversity dynamics. This relationship, which suggests that both higher-order interactions and rapid evolution are critical for the origin and maintenance of microbial diversity, has broad-ranging links to the maintenance of diversity in other systems.

Divergent Abiotic Stressors Drive Grassland Community Assembly of Tibet and Mongolia Plateau

Multiple ecological processes simultaneously govern community assembly, but it remains unclear how abiotic stressors regulate the relative importance of these processes among different biogeographic regions. Therefore, we conducted a comprehensive study on the responses of community assembly to varying environmental gradients, using the mean, variance, skewness, and kurtosis of plant height (height), specific leaf area (SLA) and leaf dry matter content (LDMC) distributions on the Tibetan Plateau (TP) and the Mongolian Plateau (MP). Our results showed that the prevalence of trait convergence across all grasslands in both TP and MP seem to be the result of abiotic filtering or weaker competitive exclusion etc. These trait-convergence assembly processes decrease the functional dispersion but increase the evenness of the trait frequency distribution. The mean, variance, skewness, and kurtosis responses of grassland communities to abiotic stress varied between the TP and MP. On average, plant trait distribution was mainly driven by temperature on the TP, and low-temperature stress altered the community assembly rules. In contrast, water availability shaped plant trait frequency distributions on the MP, and drought stress mediated the balance between different assembly processes. Our results provide empirical evidence that divergent abiotic stressors regulate the grassland community assembly on the TP and MP. Together, our study speculates that different aspects of future climate change, such as climate warming and changing precipitation patterns, on community assembly are dependent on regional climatic regimes.

Dietary Fiber Influences Bacterial Community Assembly Processes in the Gut Microbiota of Durco × Bamei Crossbred Pig

Several studies have shown that dietary fiber can significantly alter the composition and structure of the gut bacterial community in humans and mammals. However, few researches have been conducted on the dynamics of the bacterial community assembly across different graded levels of dietary fiber in different gut regions. To address this, 24 Durco × Bamei crossbred pigs were randomly assigned to four experimental chows comprising graded levels of dietary fiber. Results showed that the α-and β-diversity of the bacterial community was significantly different between the cecum and the jejunum. Adding fiber to the chow significantly increased the α-diversity of the bacterial community in the jejunum and cecum, while the β-diversity decreased. The complexity of the bacterial network increased with the increase of dietary fiber in jejunal content samples, while it decreased in cecal content samples. Furthermore, we found that stochastic processes governed the bacterial community assembly of low and medium dietary fiber groups of jejunal content samples, while deterministic processes dominated the high fiber group. In addition, deterministic processes dominated all cecal content samples. Taken together, the variation of gut community composition and structure in response to dietary fiber was distinct in different gut regions, and the dynamics of bacterial community assembly across the graded levels of dietary fiber in different gut regions was also distinct. These findings enhanced our knowledge on the bacterial community assembly processes in gut ecosystems of livestock.

Biochar Addition Altered Bacterial Community and Improved Photosynthetic Rate of Seagrass: A Mesocosm Study of Seagrass Thalassia hemprichii

Seagrass meadows, as typical “blue carbon” ecosystems, play critical ecological roles in the marine ecosystem and decline every year. The application of biochar in soil has been proposed as a potential soil amendment to improve soil quality and mitigate global climate change. The effects of biochar on soil bacterial activities are integrally linked to the potential of biochar in achieving these benefits. However, biochar has been rarely applied in marine ecosystems. Whether the application of biochar could work on the seagrass ecosystem remained unknown. In this study, we investigated the responses of sediment and rhizosphere bacterial communities of seagrass Thalassia hemprichii to the biochar addition derived from maize at ratios of 5% by dry weight in the soil during a one-month incubation. Results indicated that the biochar addition significantly changed the sedimental environment with increasing pH, total phosphorus, and total kalium while total nitrogen decreased. Biochar addition significantly altered both the rhizosphere and sediment bacterial community compositions. The significant changes in rhizosphere bacterial community composition occurred after 30days of incubation, while the significant variations in sediment bacterial community composition distinctly delayed than in sediment occurred on the 14th day. Biochar application improved nitrification and denitrification, which may accelerate nitrogen cycling. As a stabilizer to communities, biochar addition decreased the importance of deterministic selection in sediment and changed the bacterial co-occurrence pattern. The biochar addition may promote seagrass photosynthesis and growth by altering the bacterial community compositions and improving nutrient circulation in the seagrass ecosystem, contributing to the seagrass health improvement. This study provided a theoretical basis for applying biochar to the seagrass ecosystem and shed light on the feasible application of biochar in the marine ecosystem.

Graphical Abstract

Why Do Antibiotics Exist?

In the struggle with antibiotic resistance, we are losing. There is now a serious threat of moving into a postantibiotic world. High levels of resistance, in terms of both frequency and strength, have evolved against all clinically approved antibiotics worldwide. The usable life span of new clinically approved antibiotics is typically less than a decade before resistance reaches frequencies so high as to require only guarded usage. However, microbes have produced antibiotics for millennia without resistance becoming an existential issue. If resistance is the inevitable consequence of antibiotic usage, as has been the human experience, why has it not become an issue for microbes as well, especially since resistance genes are as prevalent in nature as the genes responsible for antibiotic production? Here, we ask how antibiotics can exist given the almost ubiquitous presence of resistance genes in the very microbes that have produced and used antibiotics since before humans walked the planet. We find that the context of both production and usage of antibiotics by microbes may be key to understanding how resistance is managed over time, with antibiotic synthesis and resistance existing in a paired relationship, much like a cipher and key, that impacts microbial community assembly. Finally, we put forward the cohesive, ecologically based "secret society" hypothesis to explain the longevity of antibiotics in nature.

Soil Microbial Resource Limitations and Community Assembly Along a Camellia oleifera Plantation Chronosequence

Understanding soil microbial element limitation and its relation with the microbial community can help in elucidating the soil fertility status and improving nutrient management of planted forest ecosystems. The stand age of a planted forest determines the aboveground forest biomass and structure and underground microbial function and diversity. In this study, we investigated 30 plantations of Camellia oleifera distributed across the subtropical region of China that we classified into four stand ages (planted <9 years, 9–20 years, 21–60 years, and >60 years age). Enzymatic stoichiometry analysis showed that microbial metabolism in the forests was mainly limited by C and P. P limitation significantly decreased and C limitation slightly increased along the stand age gradient. The alpha diversity of the soil microbiota remained steady along stand age, while microbial communities gradually converged from scattered to clustered, which was accompanied by a decrease in network complexity. The soil bacterial community assembly shifted from stochastic to deterministic processes, which probably contributed to a decrease in soil pH along stand age. Our findings emphasize that the stand age regulated the soil microbial metabolism limitation and community assembly, which provides new insight into the improvement of C and P management in subtropical planted forest.

Patterns and Drivers of nirK-Type and nirS-Type Denitrifier Community Assembly along an Elevation Gradient

nirK-type and nirS-type denitrifier communities mediate the conversion of nitrite to nitric oxide, which is the key step in denitrification. Results of previous studies have indicated that nirK-type and nirS-type denitrifiers may occupy different niches; however, the mechanisms and drivers of their responses to environmental changes within community assembly are poorly understood. In this study, we evaluated the distribution and assembly of nirK-type and nirS-type denitrifier communities along an elevation gradient from 1,800 to 4,100 m at Mount Gongga, China. Results showed that elevational patterns of alpha diversity in nirK-type and nirS-type denitrifier communities followed hump-backed patterns along the elevation gradient. However, nirK-type denitrifier communities formed two distinct clusters that were primarily separated by elevation, whereas nirS-type denitrifier communities formed three distinct clusters that were primarily separated by forest type along the elevation gradient. Moreover, deterministic processes were dominant in governing the assemblages of nirK-type and nirS-type denitrifiers. Soil pH was a key factor influencing the alpha and beta diversity of the nirK-type denitrifier communities, whereas plant richness was a primary variable influencing nirS-type denitrifiers. Additionally, our work revealed that soil denitrification potential was mainly explained by the variation in the beta diversity of denitrifier communities rather than the alpha diversity of denitrifier communities or denitrifier abundances over a large elevation gradient, and nirK-type denitrifiers played more important roles in soil denitrification. These results may contribute to predicting the consequences of global changes on denitrifier communities and their ecological services. IMPORTANCE Mount Gongga is the highest peak in the Hengduan Mountain region and is located at the southeastern fringe of the Tibetan Plateau, Sichuan Province, southwest China. As a transitional zone between the Tibetan Plateau and Sichuan Basin, Gongga Mountain features particularly diverse topography, geology, climate, and biodiversity and is a globally significant hot spot of biodiversity. In this contribution, we comprehensively describe the diversity and assembly of denitrifier communities along an elevation gradient on Gongga Mountain. Our findings established for the first time that the distribution patterns of beta diversity and driving factors differed between nirK-type and nirS-type denitrifier communities, and deterministic processes were dominant in shaping communities of denitrifiers. Moreover, the beta diversity of denitrifier communities rather than alpha diversity or denitrifier abundance played an important role in explaining denitrification potential, and the beta diversity of nirK-type denitrifier communities was more important than nirS-type denitrifier communities in soil denitrification. This work provides crucial insights into the spatial distribution of denitrifier communities and their ecological function and increases our understanding of the mechanisms underlying spatial distribution of community assembly along large elevation gradients.

Contrasting Community Assembly Forces Drive Microbial Structural and Potential Functional Responses to Precipitation in an Incipient Soil System

Microbial communities in incipient soil systems serve as the only biotic force shaping landscape evolution. However, the underlying ecological forces shaping microbial community structure and function are inadequately understood. We used amplicon sequencing to determine microbial taxonomic assembly and metagenome sequencing to evaluate microbial functional assembly in incipient basaltic soil subjected to precipitation. Community composition was stratified with soil depth in the pre-precipitation samples, with surficial communities maintaining their distinct structure and diversity after precipitation, while the deeper soil samples appeared to become more uniform. The structural community assembly remained deterministic in pre- and post-precipitation periods, with homogenous selection being dominant. Metagenome analysis revealed that carbon and nitrogen functional potential was assembled stochastically. Sub-populations putatively involved in the nitrogen cycle and carbon fixation experienced counteracting assembly pressures at the deepest depths, suggesting the communities may functionally assemble to respond to short-term environmental fluctuations and impact the landscape-scale response to perturbations. We propose that contrasting assembly forces impact microbial structure and potential function in an incipient landscape; in situ landscape characteristics (here homogenous parent material) drive community structure assembly, while short-term environmental fluctuations (here precipitation) shape environmental variations that are random in the soil depth profile and drive stochastic sub-population functional dynamics.

Functional recovery of secondary tropical forests

One-third of all Neotropical forests are secondary forests that regrow naturally after agricultural use through secondary succession. We need to understand better how and why succession varies across environmental gradients and broad geographic scales. Here, we analyze functional recovery using community data on seven plant characteristics (traits) of 1,016 forest plots from 30 chronosequence sites across the Neotropics. By analyzing communities in terms of their traits, we enhance understanding of the mechanisms of succession, assess ecosystem recovery, and use these insights to propose successful forest restoration strategies. Wet and dry forests diverged markedly for several traits that increase growth rate in wet forests but come at the expense of reduced drought tolerance, delay, or avoidance, which is important in seasonally dry forests. Dry and wet forests showed different successional pathways for several traits. In dry forests, species turnover is driven by drought tolerance traits that are important early in succession and in wet forests by shade tolerance traits that are important later in succession. In both forests, deciduous and compound-leaved trees decreased with forest age, probably because microclimatic conditions became less hot and dry. Our results suggest that climatic water availability drives functional recovery by influencing the start and trajectory of succession, resulting in a convergence of community trait values with forest age when vegetation cover builds up. Within plots, the range in functional trait values increased with age. Based on the observed successional trait changes, we indicate the consequences for carbon and nutrient cycling and propose an ecologically sound strategy to improve forest restoration success.

Ectomycorrhizal fungal decay traits along a soil nitrogen gradient

The extent to which ectomycorrhizal (ECM) fungi decay soil organic matter (SOM) has implications for accurately predicting forest ecosystem response to climate change. Investigating the distribution of gene traits associated with SOM decay among ectomycorrhizal fungal communities could improve understanding of SOM dynamics and plant nutrition. We hypothesized that soil inorganic nitrogen (N) availability structures the distribution of ECM fungal genes associated with SOM decay and, specifically, that ECM fungal communities occurring in inorganic N-poor soils have greater SOM decay potential. To test this hypothesis, we paired amplicon and shotgun metagenomic sequencing of 60 ECM fungal communities associating with Quercus rubra along a natural soil inorganic N gradient. Ectomycorrhizal fungal communities occurring in low inorganic N soils were enriched in gene families involved in the decay of lignin, cellulose, and chitin. Ectomycorrhizal fungal community composition was the strongest driver of shifts in metagenomic estimates of fungal decay potential. Our study simultaneously illuminates the identity of key ECM fungal taxa and gene families potentially involved in the decay of SOM, and we link rhizomorphic and medium-distance hyphal morphologies with enhanced SOM decay potential. Coupled shifts in ECM fungal community composition and community-level decay gene frequencies are consistent with outcomes of trait-mediated community assembly processes.

The evolutionary assembly of forest communities along environmental gradients: recent diversification or sorting of pre-adapted clades?

Recent studies have demonstrated that ecological processes that shape community structure and dynamics change along environmental gradients. However, much less is known about how the emergence of the gradients themselves shape the evolution of species that underlie community assembly. In this study, we address how the creation of novel environments leads to community assembly via two nonmutually exclusive processes: immigration and ecological sorting of pre-adapted clades (ISPC), and recent adaptive diversification (RAD). We study these processes in the context of the elevational gradient created by the uplift of the Central Andes. We develop a novel approach and method based on the decomposition of species turnover into within- and among-clade components, where clades correspond to lineages that originated before mountain uplift. Effects of ISPC and RAD can be inferred from how components of turnover change with elevation. We test our approach using data from over 500 Andean forest plots. We found that species turnover between communities at different elevations is dominated by the replacement of clades that originated before the uplift of the Central Andes. Our results suggest that immigration and sorting of clades pre-adapted to montane habitats is the primary mechanism shaping tree communities across elevations.

Opposing community assembly patterns for dominant and nondominant plant species in herbaceous ecosystems globally

Biotic and abiotic factors interact with dominant plants-the locally most frequent or with the largest coverage-and nondominant plants differently, partially because dominant plants modify the environment where nondominant plants grow. For instance, if dominant plants compete strongly, they will deplete most resources, forcing nondominant plants into a narrower niche space. Conversely, if dominant plants are constrained by the environment, they might not exhaust available resources but instead may ameliorate environmental stressors that usually limit nondominants. Hence, the nature of interactions among nondominant species could be modified by dominant species. Furthermore, these differences could translate into a disparity in the phylogenetic relatedness among dominants compared to the relatedness among nondominants. By estimating phylogenetic dispersion in 78 grasslands across five continents, we found that dominant species were clustered (e.g., co-dominant grasses), suggesting dominant species are likely organized by environmental filtering, and that nondominant species were either randomly assembled or overdispersed. Traits showed similar trends for those sites (<50%) with sufficient trait data. Furthermore, several lineages scattered in the phylogeny had more nondominant species than expected at random, suggesting that traits common in nondominants are phylogenetically conserved and have evolved multiple times. We also explored environmental drivers of the dominant/nondominant disparity. We found different assembly patterns for dominants and nondominants, consistent with asymmetries in assembly mechanisms. Among the different postulated mechanisms, our results suggest two complementary hypotheses seldom explored: (1) Nondominant species include lineages adapted to thrive in the environment generated by dominant species. (2) Even when dominant species reduce resources to nondominant ones, dominant species could have a stronger positive effect on some nondominants by ameliorating environmental stressors affecting them, than by depleting resources and increasing the environmental stress to those nondominants. These results show that the dominant/nondominant asymmetry has ecological and evolutionary consequences fundamental to understand plant communities.

Patch size drives colonization by aquatic insects, with minor priority effects of a cohabitant

Patch size is one of the most important factors affecting the distribution and abundance of species, and recent research has shown that patch size is an important niche dimension affecting community structure in aquatic insects. Building on this result, we examined the impact of patch size in conjunction with presence of larval anurans on colonization by aquatic insects. Hyla chrysoscelis (Cope's gray treefrog) larvae are abundant and early colonists in fishless lentic habitats, and these larvae can fill multiple ecological roles. By establishing larvae in mesocosms prior to colonization, we were able to assess whether H. chrysoscelis larvae have priority effects on aquatic insect assemblages. We conducted a series of three experiments in naturally colonized experimental landscapes to test whether (1) H. chrysoscelis larval density affects insect colonization, (2) variation in patch size affects insect colonization, and (3) the presence and larval density of H. chrysoscelis shift colonization of insects between patches of different size. Larval density independently had almost no effect on colonization, while patch size had species-specific effects consistent with prior work. When larvae and patch size were tested in conjunction, patch size had numerous, often strong, species-specific effects on colonization; larval density had effects largely limited to the assemblages of colonizing beetles and water bugs, with few effects on individual species. Higher larval densities in large mesocosms shifted some insect colonization to smaller patches, resulting in higher beta diversity among small patches in proximity to high density large mesocosms. This indicates establishing H. chrysoscelis larvae prior to insect colonization can likely create priority effects that slightly shape insect communities. Our results support the importance of patch size in studying species abundances and distributions and also indicate that colonization order plays an important role in determining the communities found within habitat patches.

Diversity Distribution, Driving Factors and Assembly Mechanisms of Free-Living and Particle-Associated Bacterial Communities at a Subtropical Marginal Sea

Free-living (FL) and particle-associated (PA) bacterioplankton communities play critical roles in biogeochemical cycles in the ocean. However, their community composition, assembly process and functions in the continental shelf and slope regions are poorly understood. Based on 16S rRNA gene amplicon sequencing, we investigated bacterial communities’ driving factors, assembly processes and functional potentials at a subtropical marginal sea. The bacterioplankton community showed specific distribution patterns with respect to lifestyle (free living vs. particle associated), habitat (slope vs. shelf) and depth (surface vs. DCM and Bottom). Salinity and water temperature were the key factors modulating turnover in the FL community, whereas nitrite, silicate and phosphate were the key factors for the PA community. Model analyses revealed that stochastic processes outweighed deterministic processes and had stronger influences on PA than FL. Homogeneous selection (Hos) was more responsible for the assembly and turnover of FL, while drift and dispersal limitation contributed more to the assembly of PA. Importantly, the primary contributor to Hos in PA was Gammaproteobacteria:Others, whereas that in FL was Cyanobacteria:Bin6. Finally, the PICRUSt2 analysis indicated that the potential metabolisms of carbohydrates, cofactors, amino acids, terpenoids, polyketides, lipids and antibiotic resistance were markedly enriched in PA than FL.

Decay of similarity across tropical forest communities: integrating spatial distance with soil nutrients

Understanding the mechanisms that drive the change of biotic assemblages over space and time is the main quest of community ecology. Assessing the relative importance of dispersal and environmental species selection in a range of organismic sizes and motilities has been a fruitful strategy. A consensus for whether spatial and environmental distances operate similarly across spatial scales and taxa, however, has yet to emerge. We used censuses of four major groups of organisms (soil bacteria, fungi, ground insects, and trees) at two observation scales (1-m² sampling point vs. 2,500-m² plots) in a topographically standardized sampling design replicated in two tropical rainforests with contrasting relationships between spatial distance and nutrient availability. We modeled the decay of assemblage similarity for each taxon set and site to assess the relative contributions of spatial distance and nutrient availability distance. Then, we evaluated the potentially structuring effect of tree composition over all other taxa. The similarity of nutrient content in the litter and topsoil had a stronger and more consistent selective effect than did dispersal limitation, particularly for bacteria, fungi, and trees at the plot level. Ground insects, the only group assessed with the capacity of active dispersal, had the highest species turnover and the flattest nonsignificant distance-decay relationship, suggesting that neither dispersal limitation nor nutrient availability were fundamental drivers of their community assembly at this scale of analysis. Only the fungal communities at one of our study sites were clearly coordinated with tree composition. The spatial distance at the smallest scale was more important than nutrient selection for the bacteria, fungi, and insects. The lower initial similarity and the moderate variation in composition identified by these distance-decay models, however, suggested that the effects of stochastic sampling were important at this smaller spatial scale. Our results highlight the importance of nutrients as one of the main environmental drivers of rainforest communities irrespective of organismic or propagule size and how the overriding effect of the analytical scale influences the interpretation, leading to the perception of greater importance of dispersal limitation and ecological drift over selection associated with environmental niches at decreasing observation scales.

Community Assembly and Co-occurrence Patterns Underlying the Core and Satellite Bacterial Sub-communities in the Tibetan Lakes

Microbial communities normally comprise a few core species and large numbers of satellite species. These two sub-communities have different ecological and functional roles in natural environments, but knowledge on the assembly processes and co-occurrence patterns of the core and satellite species in Tibetan lakes is still sparse. Here, we investigated the ecological processes and co-occurrence relationships of the core and satellite bacterial sub-communities in the Tibetan lakes via 454 sequencing of 16S rRNA gene. Our studies indicated that the core and satellite bacterial sub-communities have similar dominant phyla (Proteobacteria, Bacteroidetes, and Actinobacteria). But the core sub-communities were less diverse and exhibited a stronger distance-decay relationship than the satellite sub-communities. In addition, topological properties of nodes in the network demonstrated that the core sub-communities had more complex and stable co-occurrence associations and were primarily driven by stochastic processes (58.19%). By contrast, the satellite sub-communities were mainly governed by deterministic processes (62.17%). Overall, this study demonstrated the differences in the core and satellite sub-community assembly and network stability, suggesting the importance of considering species traits to understand the biogeographic distribution of bacterial communities in high-altitude lakes.

Integrating Stochastic and Deterministic Process in the Biogeography of N2-Fixing Cyanobacterium Candidatus Atelocyanobacterium Thalassa

UCYN-A is one of the most widespread and important marine diazotrophs. Its unusual distribution in both cold/warm and coastal/oceanic waters challenges current understanding about what drives the biogeography of diazotrophs. This study assessed the community assembly processes of the nitrogen-fixing cyanobacterium UCYN-A, developing a framework of assembly processes underpinning the microbial biogeography and diversity. High-throughput sequencing and a qPCR approach targeting the nifH gene were used to investigate three tropical seas: the Bay of Bengal, the Western Pacific Ocean, and the South China Sea. Based on the neutral community model and two types of null models calculating the β-nearest taxon index and the normalized stochasticity ratio, we found that stochastic assembly processes could explain 66-92% of the community assembly; thus, they exert overwhelming influence on UCYN-A biogeography and diversity. Among the deterministic processes, temperature and coastal/oceanic position appeared to be the principal environmental factors driving UCYN-A diversity. In addition, a close linkage between assembly processes and UCYN-A abundance/diversity/drivers can provide clues for the unusual global distribution of UCYN-A.

Species Distribution Pattern and Their Contribution in Plant Community Assembly in Response to Ecological Gradients of the Ecotonal Zone in the Himalayan Region

The ecotonal zones support populations that are acclimated to changing, fluctuating, and unstable conditions, and as a result, these populations are better equipped to adjust to expected change. In this context, a hypothesis was tested that there must be vegetation dominated by unique indicator plant species under the influence of ecological gradients in the ecotonal zone of Manoor Valley (northwestern Himalaya), Pakistan. Keeping the aforementioned hypothesis in mind, detailed field studies were conducted during different seasons in 2015-18. Line transect sampling and phytosociological characteristics (density, frequency, cover, and their relative values and Importance Value) were implemented as ecological methods. This investigation documented 97 plant species recorded from seven sampling sites. The community distribution modelling revealed that the ecological variables separate the seven sampling sites into two major plant communities (Indigofera-Parrotiopsis-Bistorta and Ziziphus-Leptopus-Quercus) recognized by TWINSPAN. The IBP communities showed a positive and significant correlation with altitude (1789.6-1896.3 m), sandy soil texture with a slightly acidic pH (6.4-6.5), and higher phosphorous (9-13 mg kg^-1). In contrast with this, the ZLQ community was recognized on the southern slope under the strong influence of high electrical conductivity (2.82-5.4 dsm^-1), organic matter (1.08-1.25%), calcium carbonate (5.8-7.6 mg kg^-1), potassium (202-220 mg kg^-1), and temperature (28.8-31.8 °C). Hence, both communities were found on opposite axes with clear differences based on the ecological gradients. NMDS clustered different species with similar habitats and different stands with common species, showing that plant species and stands were in a linear combination with ecological gradients. The IPB community has the maximum number of plant species (87 species), Shannon value (H' = 4), Simpson value (0.98), and Pielou's evenness value (0.96). Thus, the multivariate approaches revealed unique vegetation with sharp boundaries between communities which might be due to abrupt environmental changes.

The Role Transition of Dietary Species Richness in Modulating the Gut Microbial Assembly and Postweaning Performance of a Generalist Herbivore

When facing a food shortage, generalist herbivores can respond by expanding their dietary species richness (DSR) to maximize energy collection, regardless of whether forages are preferred or not. Higher DSR usually indicates higher nutrient adequacy and better health. However, the high-DSR diet containing a large proportion of preferred species or a large proportion of less-preferred species means different things to an animal. It is still unknown how different shift patterns in DSR would affect distinctly the performance of animals via altering gut microbiota. We examined the gut microbial composition, diversity, community assembly processes, and performance of a generalist herbivore, Lasiopodomys brandtii, in a feeding experiment with increased levels of simulated DSR shifting from preferred plant species to less preferred ones. We found the survival rate and body growth of Brandt's voles showed a dome-shaped association with DSR: species performance increased initially with the increase of preferred plant species but declined with the increase of less-preferred food items. Several microbial taxa and functions closely related to the metabolism of amino acids and short-chain fatty acids also showed a dome-shaped association with DSR, which is consistent with the observation of performance change. However, the alpha diversities of gut microbiota increased linearly with DSR. The null model and phylogenetic analysis suggested that stochastic processes dominate at low DSR diets, whereas deterministic processes prevail at high DSR diets. These results suggest that the role of DSR in regulating animal performance by gut microbiota depends on the number of preferred forage items. IMPORTANCE The plant species diversity varies greatly under the influence of both climate change and human disturbance, which may negatively affect the productivity as well as the variability of organisms (e.g., small herbivores) at the next trophic level. It is still unknown how gut microbiota of small herbivores respond to such changes in dietary species richness. Our manipulative food experiment revealed that dietary species richness can affect the composition, functions, and community assembly of gut microbiota of Brandt's vole in a nonlinear way. Given the fast-growing interest in therapeutic diets to treat dysbiosis and to improve health conditions, our study highlights the need to consider not just the variety of consumed food but also the principles of rational nutrition.

High historical variability weakens the effects of current climate differentiation on microbial community dissimilarity and assembly

Understanding the influences of global climate change on soil microbial communities is essential in evaluating the terrestrial biosphere's feedback to this alarming anthropogenic disturbance. However, little is known about how intra-site historical climate variability can mediate the influences of current climate differences on community dissimilarity and assembly. To fill this gap, we examined and disentangled the interactive effects of historical climate variability and current climate differences on the soil bacterial community dissimilarity and stochasticity of community assembly among 143 sites from 28 forests across eastern China. We hypothesize that the relative importance of stochasticity and community dissimilarity are related to historical climate variability and that an increasing sum of intra-site historical variability enhances stochasticity while reduces dissimilarity between two communities. To test our hypothesis, we statistically controlled for covariates between sites including differences in soil chemistry, plant diversity, spatial distance, and seasonal climate variations at annual timescales. We observed that an increase in inter-site current climate differences led to a reduced impact of stochasticity in community assembly and a pronounced divergence between communities. In stark contrast, when communities were subjected to a high level of intra-site historical climate fluctuation, the observed impact incurred from current climate differences was substantially weakened. Moreover, the influence of increased historical variability was consistent along the gradient of current temperature differences between sites. However, effects induced by historical fluctuation in precipitation were disproportional and only evident when small inter-site differences were observed. Consequently, if the prior climate variability is ignored, especially regarding environmental factors like temperature, we assert that the influence current climate differentiation has on regulating community dissimilarity and assembly stochasticity will be underestimated. Together, our findings highlight the importance and need of explicitly controlling the mean and the historical variability of climate factors for the next "generation" of climate change experiments to come.