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

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

Soil acidification mediates changes in soil bacterial community assembly processes in response to agricultural intensification

Agricultural intensification is known to alter the assembly of soil microbial communities, which regulate several critical ecosystem processes. However, the underlying ecological processes driving changes in microbial community assembly, particularly at the regional scale, remain poorly understood. Using 16S rDNA sequencing, we characterized soil bacterial community assembly in three land-use types with increasing land-use intensity: open fields cultivated with main crops (CF) or vegetables (VF), and greenhouses cultivated with vegetables (VG). Compared with CF, VF and VG altered bacterial community composition and decreased spatial turnover rates of edaphic variables and bacterial communities. Bacterial community assembly was primarily governed by deterministic processes; however, bacterial communities in VF and VG were phylogenetically less clustered and more influenced by variable selection and less by dispersal limitation. Soil pH was the most important edaphic variable mediating the changes in bacterial community assembly processes induced by agricultural intensification. Specifically, decreasing soil pH led to stochastic assembly of bacterial community. Soil pH was lower in more intensively managed lands, especially in case of VG (pH range: 5.86–7.42). Overall, agricultural intensification altered soil bacterial community assembly processes, which was associated with soil acidification. These findings may have implications for improving soil quality and agroecosystem sustainability.     

The balance between deterministic and stochastic processes in structuring lake bacterioplankton community over time

One major goal in microbial ecology is to establish the importance of deterministic and stochastic processes for community assembly. This is relevant to explain and predict how diversity changes at different temporal scales. However, understanding of the relative quantitative contribution of these processes and particularly of how they may change over time is limited. Here, we assessed the importance of deterministic and stochastic processes based on the analysis of the bacterial microbiome in one alpine oligotrophic and in one subalpine mesotrophic lake, which were sampled over two consecutive years at different time scales. We found that in both lakes, homogeneous selection (i.e., a deterministic process) was the main assembly process at the annual scale and explained 66.7% of the bacterial community turnover, despite differences in diversity and temporal variability patterns between ecosystems. However, in the alpine lake, homogenizing dispersal (i.e., a stochastic process) was the most important assembly process at the short‐term (daily and weekly) sampling scale and explained 55% of the community turnover. Alpha diversity differed between lakes, and seasonal stability of the bacterial community was more evident in the oligotrophic lake than in the mesotrophic one. Our results demonstrate how important forces that govern temporal changes in bacterial communities act at different time scales. Overall, our study validates on a quantitative basis, the importance and dominance of deterministic processes in structuring bacterial communities in freshwater environments over long time scales.     

环境选择和扩散限制驱动温带森林土壤细菌群落的构建

环境选择和扩散限制是生态系统中生物群落构建的两个基本过程,而两者相对作用的大小因研究尺度、群落属性和类型等有所不同.目前对温带亚高山森林土壤微生物群落构建的驱动因子和机制尚缺乏了解.本文利用PCR-DGGE技术研究庞泉沟自然保护区内5种典型森林包括华北落叶松林、青杄林、白杄林、油松林以及桦树林的6个森林土壤细菌群落(Lp MC1、Lp MC2、Pw MC、Pm MC、Pt MC、BMC)的结构特征及其影响因素,分析细菌群落结构与环境因子的相关性,以及土壤因子、植被和空间因素对细菌群落结构的影响.结果表明:研究区各样地土壤细菌群落的结构和生物多样性具有显著差异,低海拔落叶松和油松土壤细菌群落多样性较高(20条带),白杄林土壤细菌群落(13条带)多样性最低,高海拔落叶松土壤细菌群落多样性最高;土壤环境因子,如pH、土壤含水量、总碳、总氮、土壤有机质、速效磷以及土壤酶活性与土壤细菌群落多样性和结构显著相关;样地土壤细菌群落的beta多样性与群落的空间距离呈显著相关,表明扩散限制对群落结构具有一定的影响;方差分解分析结果显示,6个样地细菌群落结构的驱动因素大小依次为土壤因子(0.27)、空间因素(0.19)和植被(0.15);将区域土壤微生物作为"源群落",微宇宙试验结果显示,土壤因子是细菌群落结构形成的主要驱动力(0.35),同时源群落丰富的物种多样性对微宇宙土壤细菌群落结构具有显著影响.总之,在局域尺度下,环境选择对温带森林土壤细菌群落结构动态和多样性发挥主导作用,地理距离对群落结构具有显著影响,即确定性过程和随机过程共同决定局域森林土壤细菌群落结构,前者占主导地位.对于土壤细菌群落而言,扩散群落的组成和结构受到源群落的多样性特征和环境因子的双重影响.     

Environmental and biotic drivers of soil microbial β‐diversity across spatial and phylogenetic scales

Soil microbial communities play a key role in ecosystem functioning but still little is known about the processes that determine their turnover (β‐diversity) along ecological gradients. Here, we characterize soil microbial β‐diversity at two spatial scales and at multiple phylogenetic grains to ask how archaeal, bacterial and fungal communities are shaped by abiotic processes and biotic interactions with plants. We characterized microbial and plant communities using DNA metabarcoding of soil samples distributed across and within eighteen plots along an elevation gradient in the French Alps. The recovered taxa were placed onto phylogenies to estimate microbial and plant β‐diversity at different phylogenetic grains (i.e. resolution). We then modeled microbial β‐diversities with respect to plant β‐diversities and environmental dissimilarities across plots (landscape scale) and with respect to plant β‐diversities and spatial distances within plots (plot scale). At the landscape scale, fungal and archaeal β‐diversities were mostly related to plant β‐diversity, while bacterial β‐diversities were mostly related to environmental dissimilarities. At the plot scale, we detected a modest covariation of bacterial and fungal β‐diversities with plant β‐diversity; as well as a distance–decay relationship that suggested the influence of ecological drift on microbial communities. In addition, the covariation between fungal and plant β‐diversity at the plot scale was highest at fine or intermediate phylogenetic grains hinting that biotic interactions between those clades depends on early‐evolved traits. Altogether, we show how multiple ecological processes determine soil microbial community assembly at different spatial scales and how the strength of these processes change among microbial clades. In addition, we emphasized the imprint of microbial and plant evolutionary history on today's microbial community structure.     

Proportion of mycorrhiza-associated trees mediates community assemblages of soil fungi but not of bacteria

Recent studies have shown that mycorrhizal trees can greatly influence soil microbial communities, which in turn play important roles in the function offorest ecosystems. However, there is lack of understanding how the composition of trees with different mycorrhizal types affects soil microbial communities. Here, we collected 1606 soil samples from a 25-ha subtropical forest plot to investigate how the proportion of arbuscular mycorrhizal (AM) versus ectomycorrhizal (EcM) trees mediated soil microbial assemblages. Results showed the alpha diversities of both soil fungal and bacterial communities were significantly positively correlated with the ratio of AM/EcM trees. The AM/EcM tree ratio was important to the fungal community assembly, whereas soil pH was key to the bacterial communities. The increase in the AM/EcM tree ratio decreased the importance of stochastic forces in assembling fungal communities, while it had no significant effect on the bacterial communities. The differential importance of the AM/EcM tree ratio to fungal and bacterial communities highlights the role of mycorrhiza-associated tree composition in regulating soil microbial communities. This finding suggests that forests with different AM/EcM tree ratios would have different soil microbial communities, potentially leading to differences in soil nutrient cycling and in return different tree diversity and forest productivity.     

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

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

Environmental degradation results in contrasting changes in the assembly processes of stream bacterial and fungal communities

Environmental degradation may have strong effects on community assembly processes. We examined the assembly of bacterial and fungal communities in anthropogenically altered and near‐pristine streams. Using pyrosequencing of bacterial and fungal DNA from decomposed alder Alnus incana leaves, we specifically examined if environmental degradation deterministically decreases or increases the compositional turnover of bacterial and fungal communities. Our results showed that near‐pristine streams and anthropogenically altered streams supported distinct fungal and bacterial communities. The mechanisms assembling these communities were different in near‐pristine and altered environments. Environmental disturbance homogenized bacterial communities, whereas fungal communities were more dissimilar in disturbed sites than in near‐pristine sites. Compositional variation of both bacteria and fungi was related to water chemistry variables in disturbed sites, further implying the influence of environmental degradation on community assembly. Bacterial and fungal communities in near‐pristine streams were weakly controlled by environmental factors, suggesting that the relative importance of niche‐based versus neutral processes in assembling microbial communities may strongly depend on the spatial scale and local environmental context. Our results thus suggest that environmental degradation may strongly affect the composition and β‐diversity of stream microbial communities colonizing leaf litter, and that the direction of the change can be different between bacteria and fungi. A better understanding of the environmental tolerances of microbes and the mechanisms assembling microbial communities in natural environmental settings is needed to predict how environmental alteration is likely to affect microbial communities.     

Body size determines soil community assembly in a tropical forest

Tropical forests shelter an unparalleled biological diversity. The relative influence of environmental selection (i.e., abiotic conditions, biotic interactions) and stochastic–distance‐dependent neutral processes (i.e., demography, dispersal) in shaping communities has been extensively studied for various organisms, but has rarely been explored across a large range of body sizes, in particular in soil environments. We built a detailed census of the whole soil biota in a 12‐ha tropical forest plot using soil DNA metabarcoding. We show that the distribution of 19 taxonomic groups (ranging from microbes to mesofauna) is primarily stochastic, suggesting that neutral processes are prominent drivers of the assembly of these communities at this scale. We also identify aluminium, topography and plant species identity as weak, yet significant drivers of soil richness and community composition of bacteria, protists and to a lesser extent fungi. Finally, we show that body size, which determines the scale at which an organism perceives its environment, predicted the community assembly across taxonomic groups, with soil mesofauna assemblages being more stochastic than microbial ones. These results suggest that the relative contribution of neutral processes and environmental selection to community assembly directly depends on body size. Body size is hence an important determinant of community assembly rules at the scale of the ecological community in tropical soils and should be accounted for in spatial models of tropical soil food webs.     

Biogeographic responses and niche occupancy of microbial communities following long-term land-use change

Understanding the effects of forest-to-agriculture conversion on microbial diversity has been a major goal in soil ecological studies. However, linking community assembly to the ruling ecological processes at local and regional scales remains challenging. Here, we evaluated bacterial community assembly patterns and the ecological processes governing niche specialization in a gradient of geography, seasonality, and land-use change, totaling 324 soil samples, 43 habitat characteristics (abiotic factors), and 16 metabolic and co-occurrence patterns (biotic factors), in the Brazilian Atlantic Rainforest, a subtropical biome recognized as one the world’s largest and most threatened hotspots of biodiversity. Pairwise beta diversities were lower in pastures than in forest and no-till soils. Pasture communities showed a predominantly neutral model, regarding stochastic processes, with moderate dispersion, leading to biotic homogenization. Most no-till and forest microbial communities followed a niche-based model, with low rates of dispersal and weak homogenizing selection, indicating niche specialization or variable selection. Historical and evolutionary contingencies, as represented by soil type, season, and dispersal limitation were the main drivers of microbial assembly and processes at the local scale, markedly correlated with the occurrence of endemic microbes. Our results indicate that the patterns of assembly and their governing processes are dependent on the niche occupancy of the taxa evaluated (generalists or specialists). They are also more correlated with historical and evolutionary contingencies and the interactions among taxa (i.e., co-occurrence patterns) than the land-use change itself.

     

Environmental changes affect the assembly of soil bacterial community primarily by mediating stochastic processes

Both ‘species fitness difference’‐based deterministic processes, such as competitive exclusion and environmental filtering, and ‘species fitness difference’‐independent stochastic processes, such as birth/death and dispersal/colonization, can influence the assembly of soil microbial communities. However, how both types of processes are mediated by anthropogenic environmental changes has rarely been explored. Here we report a novel and general pattern that almost all anthropogenic environmental changes that took place in a grassland ecosystem affected soil bacterial community assembly primarily through promoting or restraining stochastic processes. We performed four experiments mimicking 16 types of environmental changes and separated the compositional variation of soil bacterial communities caused by each environmental change into deterministic and stochastic components, with a recently developed method. Briefly, because the difference between control and treatment communities is primarily caused by deterministic processes, the deterministic change was quantified as (mean compositional variation between treatment and control) – (mean compositional variation within control). The difference among replicate treatment communities is primarily caused by stochastic processes, so the stochastic change was estimated as (mean compositional variation within treatment) – (mean compositional variation within control). The absolute of the stochastic change was greater than that of the deterministic change across almost all environmental changes, which was robust for both taxonomic and functional‐based criterion. Although the deterministic change may become more important as environmental changes last longer, our findings showed that changes usually occurred through mediating stochastic processes over 5 years, challenging the traditional determinism‐dominated view.     

Environmental filtering affects soil fungal community composition more than dispersal limitation at regional scales

The relative importance of dispersal limitation versus environmental filtering for community assembly has received much attention for macroorganisms. These processes have only recently been examined in microbial communities. Instead, microbial dispersal has mostly been measured as community composition change over space (i.e., distance decay). Here we directly examined fungal composition in airborne wind currents and soil fungal communities across a 40 000 km² regional landscape to determine if dispersal limitation or abiotic factors were structuring soil fungal communities. Over this landscape, neither airborne nor soil fungal communities exhibited compositional differences due to geographic distance. Airborne fungal communities shifted temporally while soil fungal communities were correlated with abiotic parameters. These patterns suggest that environmental filtering may have the largest influence on fungal regional community assembly in soils, especially for aerially dispersed fungal taxa. Furthermore, we found evidence that dispersal of fungal spores differs between fungal taxa and can be both a stochastic and deterministic process. The spatial range of soil fungal taxa was correlated with their average regional abundance across all sites, which may imply stochastic dispersal mechanisms. Nevertheless, spore volume was also negatively correlated with spatial range for some species. Smaller volume spores may be adapted to long-range dispersal, or establishment, suggesting that deterministic fungal traits may also influence fungal distributions. Fungal life-history traits may influence their distributions as well. Hypogeous fungal taxa exhibited high local abundance, but small spatial ranges, while epigeous fungal taxa had lower local abundance, but larger spatial ranges. This study is the first, to our knowledge, to directly sample air dispersal and soil fungal communities simultaneously across a regional landscape. We provide some of the first evidence that soil fungal communities are mostly assembled through environmental filtering and experience little dispersal limitation.     

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

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

Contrasting primary successional trajectories of fungi and bacteria in retreating glacier soils

Early community assembly of soil microbial communities is essential for pedogenesis and development of organic legacies. We examined fungal and bacterial successions along a well‐established temperate glacier forefront chronosequence representing ~70 years of deglaciation to determine community assembly. As microbial communities may be heavily structured by establishing vegetation, we included nonvegetated soils as well as soils from underneath four plant species with differing mycorrhizal ecologies (Abies lasiocarpa, ectomycorrhizal; Luetkea pectinata, arbuscular mycorrhizal; Phyllodoce empetriformis, ericoid mycorrhizal; Saxifraga ferruginea, nonmycorrhizal). Our main objectives were to contrast fungal and bacterial successional dynamics and community assembly as well as to decouple the effects of plant establishment and time since deglaciation on microbial trajectories using high‐throughput sequencing. Our data indicate that distance from glacier terminus has large effects on biomass accumulation, community membership, and distribution for both fungi and bacteria. Surprisingly, presence of plants rather than their identity was more important in structuring bacterial communities along the chronosequence and played only a very minor role in structuring the fungal communities. Further, our analyses suggest that bacterial communities may converge during assembly supporting determinism, whereas fungal communities show no such patterns. Although fungal communities provided little evidence of convergence in community structure, many taxa were nonrandomly distributed across the glacier foreland; similar taxon‐level responses were observed in bacterial communities. Overall, our data highlight differing drivers for fungal and bacterial trajectories during early primary succession in recently deglaciated soils.     

Comparison of soil microbial communities inhabiting vineyards and native sclerophyllous forests in central Chile

Natural ecosystems provide services to agriculture such as pest control, soil nutrients, and key microbial components. These services and others in turn provide essential elements that fuel biomass productivity. Responsible agricultural management and conservation of natural habitats can enhance these ecosystem services. Vineyards are currently driving land‐use changes in many Mediterranean ecosystems. These land‐use changes could have important effects on the supporting ecosystems services related to the soil properties and the microbial communities associated with forests and vineyard soils. Here, we explore soil bacterial and fungal communities present in sclerophyllous forests and organic vineyards from three different wine growing areas in central Chile. We employed terminal restriction fragment length polymorphisms (T‐RFLP) to describe the soil microbial communities inhabiting native forests and vineyards in central Chile. We found that the bacterial community changed between the sampled growing areas; however, the fungal community did not differ. At the local scale, our findings show that fungal communities differed between habitats because fungi species might be more sensitive to land‐use change compared to bacterial species, as bacterial communities did not change between forests and vineyards. We discuss these findings based on the sensitivity of microbial communities to soil properties and land‐use change. Finally, we focus our conclusions on the importance of naturally derived ecosystem services to vineyards.     

Abundant fungi adapt to broader environmental gradients than rare fungi in agricultural fields

Soil communities are intricately linked to ecosystem functioning, and a predictive understanding of how communities assemble in response to environmental change is of great ecological importance. Little is known about the assembly processes governing abundant and rare fungal communities across agro‐ecosystems, particularly with regard to their environmental adaptation. By considering abundant and rare taxa, we tested the environmental thresholds and phylogenetic signals for ecological preferences of fungal communities across complex environmental gradients to reflect their environmental adaptation, and explored the factors influencing their assembly based on the large‐scale soil survey in agricultural fields across eastern China. We found that the abundant taxa exhibited remarkably broader response thresholds and stronger phylogenetic signals for the ecological preferences across environmental gradients compared to the rare taxa. Neutral processes played a key role in shaping the abundant subcommunity compared to the rare subcommunity. Null model analysis revealed that the abundant subcommunity was less clustered phylogenetically and governed primarily by dispersal limitation, while homogeneous selection was the major assembly process in the rare subcommunity. Soil available sulfur was the major factor mediating the balance between stochastic and deterministic processes of both the abundant and rare subcommunities, as indicated by an increase in stochasticity with higher available sulfur concentration. Based on macroecological spatial scale datasets, our study revealed the potential broader environmental adaptation of abundant fungal taxa compared to rare fungal taxa, and identified the factors mediating their distinct community assembly processes in agricultural fields. These results contribute to our understanding of the mechanisms underlying the generation and maintenance of fungal diversity in response to global environmental change.     

Distinct patterns of soil bacterial and fungal community assemblages in subtropical forest ecosystems under warming

Climate change globally affects soil microbial community assembly across ecosystems. However, little is known about the impact of warming on the structure of soil microbial communities or underlying mechanisms that shape microbial community composition in subtropical forest ecosystems. To address this gap, we utilized natural variation in temperature via an altitudinal gradient to simulate ecosystem warming. After 6 years, microbial co-occurrence network complexity increased with warming, and changes in their taxonomic composition were asynchronous, likely due to contrasting community assembly processes. We found that while stochastic processes were drivers of bacterial community composition, warming led to a shift from stochastic to deterministic drivers in dry season. Structural equation modelling highlighted that soil temperature and water content positively influenced soil microbial communities during dry season and negatively during wet season. These results facilitate our understanding of the response of soil microbial communities to climate warming and may improve predictions of ecosystem function of soil microbes in subtropical forests.     

Similar Processes but Different Environmental Filters for Soil Bacterial and Fungal Community Composition Turnover on a Broad Spatial Scale

Spatial scaling of microorganisms has been demonstrated over the last decade. However, the processes and environmental filters shaping soil microbial community structure on a broad spatial scale still need to be refined and ranked. Here, we compared bacterial and fungal community composition turnovers through a biogeographical approach on the same soil sampling design at a broad spatial scale (area range: 13300 to 31000 km²): i) to examine their spatial structuring; ii) to investigate the relative importance of environmental selection and spatial autocorrelation in determining their community composition turnover; and iii) to identify and rank the relevant environmental filters and scales involved in their spatial variations. Molecular fingerprinting of soil bacterial and fungal communities was performed on 413 soils from four French regions of contrasting environmental heterogeneity (Landes<Burgundy≤Brittany<<South-East) using the systematic grid of French Soil Quality Monitoring Network to evaluate the communities’ composition turnovers. The relative importance of processes and filters was assessed by distance-based redundancy analysis. This study demonstrates significant community composition turnover rates for soil bacteria and fungi, which were dependent on the region. Bacterial and fungal community composition turnovers were mainly driven by environmental selection explaining from 10% to 20% of community composition variations, but spatial variables also explained 3% to 9% of total variance. These variables highlighted significant spatial autocorrelation of both communities unexplained by the environmental variables measured and could partly be explained by dispersal limitations. Although the identified filters and their hierarchy were dependent on the region and organism, selection was systematically based on a common group of environmental variables: pH, trophic resources, texture and land use. Spatial autocorrelation was also important at coarse (80 to 120 km radius) and/or medium (40 to 65 km radius) spatial scales, suggesting dispersal limitations at these scales.     

Fungal infection alters the selection,dispersal and drift processes structuring the amphibian skin microbiome

Symbiotic microbial communities are important for host health, but the processes shaping these communities are poorly understood. Understanding how community assembly processes jointly affect microbial community composition is limited because inflexible community models rely on rejecting dispersal and drift before considering selection. We developed a flexible community assembly model based on neutral theory to ask: How do dispersal, drift and selection concurrently affect the microbiome across environmental gradients? We applied this approach to examine how a fungal pathogen affected the assembly processes structuring the amphibian skin microbiome. We found that the rejection of neutrality for the amphibian microbiome across a fungal gradient was not strictly due to selection processes, but was also a result of species‐specific changes in dispersal and drift. Our modelling framework brings the qualitative recognition that niche and neutral processes jointly structure microbiomes into quantitative focus, allowing for improved predictions of microbial community turnover across environmental gradients.     

Microplastics reduce soil microbial network complexity and ecological deterministic selection

Microplastics have been proposed as emerging threats for terrestrial systems as they may potentially alter the physicochemical/biophysical soil environments. Due to the variety of properties of microplastics and soils, the microplastic-induced effects in soil ecosystems are greatly manifold. Here, we studied effects of three polymer microplastics (polyamide-6, polyethylene, and polyethylene terephthalate) on soil properties with four different soil types. The success patterns, interaction relationships, and assembly processes of soil bacterial communities were also studied. Microplastics have the potential to promote CO₂ emissions and enhance the soil humification. Even though microplastics did not significantly alter the diversity and composition of the soil microbial community, the application of microplastics decreased the network complexity and stability, including network size, connectivity, and the number of module and keystone species. The bacterial community assembly was governed by deterministic selection (77.3%–90.9%) in all treatments, while microplastics increased the contribution of stochastic processes from 9.1% in control to 13.6%–22.7%. The neutral model results also indicated most of the bacterial taxa were present in the predicted neutral region (approximately 98%), suggesting the importance of stochastic processes. These findings provided a fundamental insight in understanding the effects of microplastics on soil ecosystems.