Whole-genome-based phylogenetic analyses provide new insights into the evolution of springtails (Hexapoda: Collembola).

Springtails (Collembola) stand as one of the most abundant, widespread, and ancient terrestrial arthropods on earth. However, their evolutionary history and deep phylogenetic relationships remain elusive. In this study, we employed phylogenomic approaches to elucidate the basal relationships among Collembola. We sampled whole-genome data representing all major collembolan lineages in proportion to their known diversity. To account for potential phylogenomic biases, we implemented various data extraction, locus sampling, and signal filtering strategies to generate matrices. Subsequently, we applied a diverse array of tree-searching and rate-modelling methods to reconstruct the phylogeny. Our analyses, utilizing different matrices and methods, converged on the same unrooted relationships among collembolan ingroups, supporting the current ordinal classification and challenging the monophyly of Arthropleona and Symphypleona s.l. However, discrepancies across analyses existed in the root of Collembola. Among various root positions, those based on more informative matrices and biologically realistic models, favoring a basal topology of Entomobryomorpha + (Symphypleona s.s. + (Neelipleona + Poduromorpha)), were supported by subsequent methodological assessment, topology tests, and rooting analyses. This optimal topology suggests multiple independent reduction of the pronotum in non-poduromorph orders and aligns with the plesiomorphic status of neuroendocrine organs and epicuticular structure of Entomobryomorpha. Fossil-calibrated dating analyses based on the optimal topology indicated late-Paleozoic to mid-Mesozoic origins of the crown Collembola and four orders. In addition, our results questioned the monophyly of Isotomidae and Neanuridae, underscoring the need for further attention to the systematics of these families. Overall, this study provides novel insights into the phylogenetic backbone of Collembola, which will inform future studies on the systematics, ecology, and evolution of this significant arthropod lineage.

[Research advances in trait-based approaches in soil animal community ecology]. / åŸºäºŽåŠŸèƒ½æ€§çŠ¶çš„åœŸå£¤åŠ¨ç‰©ç¾¤è½ç”Ÿæ€å­¦ç ”ç©¶è¿›å±•.

Functional traits are indicators of the responses and adaptation of organisms to environmental changes and cascade to a series of ecosystem functions. The functional traits of soil animals are sensitive to environmental factors and may characterize and predict the changes of ecosystem functions. Multiple dimensions of biodiversity that combing species, phylogenetic, and functional diversity improves the understanding of distribution patterns, community assembly mechanisms and ecosystem functions of soil animals. In this review, we listed the categories of soil animal functional traits and their ecological significance, and summarized current researches on the responses of soil animal communities to environmental changes and the community assembly processes based on trait-based approaches. We proposed to strengthen the study on the impacts of eco-evolution processes of biotic interactions to soil animal functional traits, establish the database of soil animal functional traits, and apply trait-based approaches in the ecological restoration in the future, which would benefit soil biodiversity conservation and sustainability of soil ecosystems.

Nutrient-induced acidification modulates soil biodiversity-function relationships.

Nutrient enrichment is a major global change component that often disrupts the relationship between aboveground biodiversity and ecosystem functions by promoting species dominance, altering trophic interactions, and reducing ecosystem stability. Emerging evidence indicates that nutrient enrichment also reduces soil biodiversity and weakens the relationship between belowground biodiversity and ecosystem functions, but the underlying mechanisms remain largely unclear. Here, we explore the effects of nutrient enrichment on soil properties, soil biodiversity, and multiple ecosystem functions through a 13-year field experiment. We show that soil acidification induced by nutrient enrichment, rather than changes in mineral nutrient and carbon (C) availability, is the primary factor negatively affecting the relationship between soil diversity and ecosystem multifunctionality. Nitrogen and phosphorus additions significantly reduce soil pH, diversity of bacteria, fungi and nematodes, as well as an array of ecosystem functions related to C and nutrient cycling. Effects of nutrient enrichment on microbial diversity also have negative consequences at higher trophic levels on the diversity of microbivorous nematodes. These results indicate that nutrient-induced acidification can cascade up its impacts along the soil food webs and influence ecosystem functioning, providing novel insight into the mechanisms through which nutrient enrichment influences soil community and ecosystem properties.

Competition-Induced Macroscopic Superlubricity of Ionic Liquid Analogues by Hydroxyl Ligands Revealed by in Situ Raman.

High load-bearing capacity is one of the crucial indicators for liquid superlubricants to move toward practicality. However, some of the current emerging systems not only have low contact pressures but also are highly susceptible to further degradation due to water adsorption and even superlubricity failure. Herein, a novel choline chloride-based ionic liquid analogues (ILAs) of a superlubricant with triethanolamine (TEOA) as the H-bond donor is reported for the first time; it obtains an ultralow coefficient of friction (0.005) and high load-bearing capacity (360 MPa, more than 2 times that of similar systems) due to adsorption of a small amount of water (<5 wt %) from the air. In situ Raman combined with 1H NMR and FTIR techniques reveals that adsorbed water competes with the hydroxyl group of TEOA for coordination with Cl-, leading to the conversion of some strong H-bonds to weak H-bonds in ILAs; the localized strong H-bonds and weak H-bonds endow the ILAs with high load-bearing capacity and the formation of ultralow shear-resistance sliding interfaces, respectively, under the shear motion. This study proposes a strategy to modulate the interactions between liquid species using adsorbed water from air as a competing ligand, which provides new insights into the design of ILA-based macroscopic liquid superlubricants with a high load-bearing capacity.

Active microbial population dynamics and life strategies drive the enhanced carbon use efficiency in high-organic matter soils.

Microbial carbon use efficiency (CUE) is a critical parameter that controls carbon storage in soil, but many uncertainties remain concerning adaptations of microbial communities to long-term fertilization that impact CUE. Based on H218O quantitative stable isotope probing coupled with metagenomic sequencing, we disentangled the roles of active microbial population dynamics and life strategies for CUE in soils after a long-term (35 years) mineral or organic fertilization. We found that the soils rich in organic matter supported high microbial CUE, indicating a more efficient microbial biomass formation and a greater carbon sequestration potential. Organic fertilizers supported active microbial communities characterized by high diversity and a relative increase in net growth rate, as well as an anabolic-biased carbon cycling, which likely explains the observed enhanced CUE. Overall, these results highlight the role of population dynamics and life strategies in understanding and predicting microbial CUE and sequestration in soil.IMPORTANCEMicrobial CUE is a major determinant of global soil organic carbon storage. Understanding the microbial processes underlying CUE can help to maintain soil sustainable productivity and mitigate climate change. Our findings indicated that active microbial communities, adapted to long-term organic fertilization, exhibited a relative increase in net growth rate and a preference for anabolic carbon cycling when compared to those subjected to chemical fertilization. These shifts in population dynamics and life strategies led the active microbes to allocate more carbon to biomass production rather than cellular respiration. Consequently, the more fertile soils may harbor a greater microbially mediated carbon sequestration potential. This finding is of great importance for manipulating microorganisms to increase soil C sequestration.

Long-term organic amendments increase the vulnerability of microbial respiration to environmental changes: Evidence from field and laboratory studies.

Organic amendments can improve soil fertility and microbial diversity, making agroecosystems more resilient to stress. However, it is uncertain whether organic amendments will enhance the functional capacity of soil microbial communities, thereby mitigating fluctuations in microbial respiration caused by environmental changes. Here, we examined the impacts of long-term organic amendments on the dynamics of microbial catabolic capacity (characterized by enzyme activities and carbon source utilization) and microbial respiration, as well as their interrelationships during a period with fluctuating temperature and rainfall in the field. We then subjected the field soil samples to laboratory heating disturbances to further evaluate the importance of microbial catabolic capacity in explaining patterns of microbial respiration. In both field and laboratory experiments, organic amendments tended to increase the stability of microbial catabolic capacity, but significantly increased the vulnerability of microbial respiration to environmental changes. However, the direction and driving factors of microbial respiration affected by environmental changes differed between the field and laboratory experiments. Environmental changes in the field suppressed the promotional effects of organic amendments on microbial respiration mainly through reducing microbial catabolic capacity, while laboratory heating further enhanced microbial respiration mainly due to increased soil resource availability. Together, these findings suggest that increased microbial respiration variations under organic amendments may potentially increase the uncertainty in predicting soil carbon emissions in the scenario of ongoing climate/anthropogenic changes, and highlight the necessity of linking laboratory studies on environmental changes to field conditions.

Highly Concentrated Electrolyte Superlubricants Enhanced by Interfacial Water Competition Around Chemically Active MgO Additives.

The low concentration of water-based lubricants and the high chemical inertness of the additives involved are often regarded as basic norms in the design of liquid lubricants. Herein, a novel liquid superlubricant of an aqueous solution containing a relatively high concentration of salt, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), is reported for the first time, and the superlubricity stability and load-bearing capacity of the optimized system (MgO0.10/LiTFSI10) are effectively strengthened by the addition of only trace (0.10 wt %) water-chemically active MgO additives. It demonstrates higher applicable loads, lower COF (∼0.004), and stability relative to the base solution. Only a trace amount of MgO additive is needed for the superlubricity, which makes up for the cost and environmental deficiencies of LiTFSI10. The weak interaction region between free water and the outer-layer water of Li+ hydration shells becomes a possible ultralow shear resistance sliding interface; the Mg(OH)2 layer, generated by the reaction of MgO with water, further creates additional weakly interacting interfaces, leading to the formation of an asymmetric contact between the clusters/particles within the hydrodynamic film by moderating the competition between interfacial water and free water, thus achieving high load-bearing macroscopic superlubricity. This study deepens the contribution of electrolyte concentration to ionic hydration and superlubricity due to the low shear slip layer formed by interfacial water competition with water-activated solid additives, providing new insights into the next generation of high load-bearing water-based liquid superlubricity systems.

Linking nematodes and ecosystem function: a trait-based framework.

Trait-based approaches are being increasingly adopted to understand species' ecological strategies and how organisms influence ecosystem function. Trait-based research on soil organisms, however, remains poorly developed compared with that for plants. The abundant and diverse soil nematodes are prime candidates to advance trait-based approaches belowground, but a unified trait framework to describe nematode ecological strategies and assess their linkages with ecosystem function is lacking. We categorized nematode traits as morphological, physiological, life history, and community clusters, and proposed the nematode economics spectrum (NES) to better understand nematode ecological strategies and their association with ecosystem function. We argue that bridging the NES and the plant economics spectrum will facilitate a more holistic understanding of ecosystem carbon and nutrient cycling under global change.

Unravelling High-Load Superlubricity of Ionic Liquid Analogues by In Situ Raman: Incomplete Hydration Induced by Competitive Exchange of External Water with Crystalline Water.

A high load-carrying capacity is the key to the practicality of liquid superlubricity, but it is difficult to achieve high load and low friction simultaneously by relying solely on a liquid film. Herein, a choline chloride-based ionic liquid analogue (ILA) macroscale superlubricant is first reported by tuning down strong hydrogen bonding in the ILA via introducing 2-10 wt % water, with a high load of 160 MPa and a low coefficient of friction of 0.006-0.008. In situ Raman reveals that competitive exchange between external water and crystalline water induces weak H-bond-dominated incomplete hydration, conferring a low-shear interface and considerable load-carrying capacity inside the lubricant. It is a hydrodynamic lubrication film rather than a tribochemical/physical adsorption film, allowing it to be applied to friction pairs of various materials. This study unveils the principle of water mediation of high-viscosity ILAs and also provides new insights into the design of practical ILA-based superlubrication materials with high load-carrying capacity.

Long-term manure inputs induce a deep selection on agroecosystem soil antibiotic resistome.

Although the enrichment of antibiotic resistance genes (ARGs) in diverse organic soils have been explored, understanding of the ecological processes governing the composition of ARGs in long-term organically fertilized soils still remains limited across typical agricultural regions. Thus, the distribution and assembly of ARG profile in three typical agricultural soils (black soil, fluvo-aquic soil, and red soil) under long-term contrasting fertilization regimes (chemical-only vs organic-only) were investigated using high-throughput qPCR (HT-qPCR). The application of organic manure significantly increased the abundance and number of ARGs across soils, as compared to those with chemical fertilizer. Organic manure application enriched the abundance of mobile genetic elements (MGEs), which were positively associated with ARGs. In addition, it is long-term organic fertilizer that enriched the number and abundance of opportunist and specialist ARGs in the fluvo-aquic and red soils, but not black soils. The number and abundance of most generalist ARGs did not change significantly among different fertilization or soil types. The assembly process of the ARG profiles tends to be more deterministic in organically fertilized soils than in chemically fertilized soils. These results suggest that long-term organic fertilizer application may contribute to the persistence and health risk of the soil antibiotic resistomes (especially specialist ARGs).

Management effects on soil nematode abundance differ among functional groups and land-use types at a global scale.

Anthropogenic land use is threatening global biodiversity. As one of the most abundant animals on Earth, nematodes occupy several key positions in belowground food webs and contribute to many ecosystem functions and services. However, the effects of land use on nematode abundance and its determinants remain poorly understood at a global scale. To characterize nematodes' responses to land use across trophic groups, we used a dataset of 6,825 soil samples globally to assess how nematode abundance varies among regional land-use types (i.e. primary vegetation, secondary vegetation, pasture, cropland and urban) and local land-use intensities (i.e. human-managed or not). We also quantified the interactive effects of land use and environmental predictors (i.e. mean annual temperature, annual precipitation, soil organic carbon, soil pH, global vegetation biomass and global vegetation productivity) on nematode abundance. We found that total nematode abundance and the abundance of bacterivores, fungivores, herbivores, omnivores and predators generally increased or were not affected under management across land-use types. Specifically, the most numerically abundant bacterivores were higher in managed than in unmanaged secondary vegetation habitats and urban areas, and herbivores were more abundant in managed than in unmanaged primary and secondary vegetation habitats. Furthermore, the numbers of significant environmental predictors of nematode abundance were reduced and the magnitude and the direction of the predictors were changed under management. We also found that nematode abundance was more variable and less determined by environmental factors in urban than in other land-use types. These findings challenge the view that human land use decreases animal abundance across trophic groups, but highlight that land use is altering the trophic composition of soil nematodes and its relationships with the environment at the global scale.

Synergies between Heat Disturbance and Inoculum Size Promote the Invasion Potential of a Bacterial Pathogen in Soil.

Inoculum size contributes to the invasion potential of pathogens in the soil. However, the role of inoculum size in determining the fate of pathogens in disturbed soils remains unclear. Herein, we investigated the survival rates of a bacterial pathogen, Ralstonia solanacearum, in soils subjected to heat as a simulated disturbance. Our results revealed that heating increased soil resource availability but reduced resource differentiation between R. solanacearum and indigenous bacterial communities. In both non-heated and heated soils, invader abundances increased with inoculum size, with a greater magnitude in heated soils. Inoculum size and heat-induced increases in soil-available carbon and nitrogen best predicted invasion success. Altogether, our findings suggested that the invasion by soil pathogens could be predicted by synergies between heat perturbation and inoculum size.

Phylogenomics of Elongate-Bodied Springtails Reveals Independent Transitions from Aboveground to Belowground Habitats in Deep Time.

Soil has become a major hotspot of biodiversity studies, yet the pattern and timing of the evolution of soil organisms are poorly known because of the scarcity of paleontological data. To overcome this limitation, we conducted a genome-based macroevolutionary study of an ancient, diversified, and widespread lineage of soil fauna, the elongate-bodied springtails (class Collembola, order Entomobryomorpha). To build the first robust backbone phylogeny of this previously refractory group, we sampled representatives of major higher taxa (6 out of 8 families, 11 out of 16 subfamilies) of the order with an emphasis on the most problematic superfamily Tomoceroidea, applied whole-genome sequencing methods, and compared the performance of different combinations of data sets (universal single-copy orthologs [USCO] vs. ultraconserved elements]) and modeling schemes. The fossil-calibrated timetree was used to reconstruct the evolution of body size, sensory organs, and pigmentation to establish a time frame of the ecomorphological divergences. The resultant trees based on different analyses were congruent in most nodes. Several discordant nodes were carefully evaluated by considering method fitness, morphological information, and topology test. The evaluation favored the well-resolved topology from analyses using USCO amino acid matrices and complex site-heterogeneous models (CAT$+$GTR and LG$+$PMSF (C60)). The preferred topology supports the monophyletic superfamily Tomoceroidea as an early-diverging lineage and a sister relationship between Entomobryoidea and Isotomoidea. The family Tomoceridae was recovered as monophyletic, whereas Oncopoduridae was recovered as paraphyletic, with Harlomillsia as a sister to Tomoceridae and hence deserving a separate family status as Harlomillsiidae Yu and Zhang fam. n. Ancestral Entomobryomorpha were reconstructed as surface-living, supporting independent origins of soil-living groups across the Paleozoic-Mesozoic, and highlighting the ancient evolutionary interaction between aboveground and belowground fauna. [Collembola; phylogenomics; soil-living adaptation; whole-genome sequencing.].

Dinotefuran alters Collembola-fungi-bacteria interactions that control mineralization of maize and soil organic carbon.

Rare studies investigated influence of neonicotinoid insecticides on the whole soil biota including non-target invertebrates and microorganisms. And less is known about the consequent intervention on soil C processes. This study aimed to decipher Collembola-fungi-bacteria interactive effects on pathways of maize C translocation, combining isotopic tracer analysis of relevant compartments with high-throughput sequencing for bacterial and fungal genetic profiles. Dinotefuran was applied at 0 or 100 µg kg-1 (a simulating residual dosage) to microcosms containing soils, Collembola and 13C labelled maize. Dinotefuran drastically reduced the density and maize-derived biomass C of Collembola, while intensifying antagonistic associations between soil organisms, with flourishing growth of Ascomycota and Actinobacteria, e.g., Streptomyces. This led to higher soil organic C (SOC) mineralization (elevated by 9.8-10.5%) across soils, attributing to the shift in microbial taxonomic and functional guild, e.g., with the increased abundance of genes aligned to cytochrome P450. Maize decomposition was controlled by Collembola that primarily fed on maize, via grazing behavior that facilitated labile maize C preferred decomposers, e.g., Xanthomonadaceae. These findings elucidate the influence of minute dinotefuran on intra-linkages between biomes (Collembola, fungi and bacteria), and highlight such legacy effects on maize and SOC mineralization.

Home-field advantage in soil respiration and its resilience to drying and rewetting cycles.

Climate change is expected to increase extreme weather events, such as more extreme drought and rainfall incidences, with consequences for ecosystem carbon (C) cycling. An understanding of how drying and rewetting (DRW) events affect microbe-mediated soil processes is therefore critical to the predictions of future climate. Here, a reciprocal-transplant experiment was conducted using two soils originated from distinct climate and agricultural managements to evaluate how soil biotic and abiotic properties regulate soil respiration and its resilience to simulated DRW cycles. We found that regardless of the DRW intensity, the effects of microbial community on soil respiration and its resilience to DRW cycles were dependent on soil type. Soil microbial communities yielded higher respiration rates and resilience in native than foreign soils under both one and four DRW cycles, supporting the "home-field advantage" hypothesis. Structural equation modeling demonstrated that soil pH and total C directly influenced soil respiration, but effects of soil abiotic properties on respiration resilience were mediated by microbial communities. Among microbial drivers, the microbial C utilization capacity (as characterized by community-level physiological profile, C-acquisition enzyme activities and microbial metabolic quotients) was the best predictor of respiration resilience to DRW cycles, followed by microbial biomass carbon/nitrogen ratio and microbial community composition. Our study suggests that soil microbial communities may have adapted to their historical conditions, which facilitates the resilience of soil respiration to changing environments, but this adaptation may accelerate C loss from soils facing increasingly variable climate.

Molecular phylogeny and trait evolution in an ancient terrestrial arthropod lineage: Systematic revision and implications for ecological divergence (Collembola, Tomocerinae).

Phylogenetic assessments of functional traits are important for mechanistically understanding the interactions between organisms and environments, but such practices are strongly limited by the availability of phylogenetic frameworks. The tomocerin springtails are an ancient, widespread and ecologically important group of terrestrial arthropods, whereas their phylogeny and trait evolution remained unaddressed. In the present study, we conducted the first comprehensive phylogenetic reconstruction of Tomocerinae, based on a multi-loci molecular dataset covering all major lineages within the subfamily, using Bayesian inference (BI), maximum-likelihood (ML) and maximum-parsimony (MP) approaches. Divergence time was estimated and ancestral character state reconstruction (ACSR) was performed to trace the evolutionary history of five ecomorphological traits correlated with sensory and locomotory functions. Our results support the monophyly of Tomocerinae, and indicate that current classification of Tomocerinae only partially reflects evolutionary relationships, notably the commonest and speciose genus Tomocerus is polyphyletic. The subfamily probably originated in Early Cretaceous and diversified in two Cretaceous and one Eocene radiation events. As indicated by the evolutionary patterns of functional traits, multiple ecological divergences took place during the diversification of Tomocerinae. The study suggests a potential underestimation of ecological divergence and functional diversity in terrestrial arthropods, calls for an update of present trait databases, and demonstrates the value of macroevolutionary knowledge for improving the trait-based ecology. In addition, Tomocerus, Tomocerina and Tritomurus are redefined, a new genus Yoshiicerusgen. n. and new subgenera Coloratomurussubgen. n., Ciliatomurussubgen. n., Striatomurussubgen. n. and Ocreatomurussubgen. n. are described in the appendix.

Trophic interactions as determinants of the arbuscular mycorrhizal fungal community with cascading plant-promoting consequences.

BACKGROUND: The soil mycobiome is composed of a complex and diverse fungal community, which includes functionally diverse species ranging from plant pathogens to mutualists. Among the latter are arbuscular mycorrhizal fungi (AMF) that provide phosphorous (P) to plants. While plant hosts and abiotic parameters are known to structure AMF communities, it remains largely unknown how higher trophic level organisms, including protists and nematodes, affect AMF abundance and community composition. RESULTS: Here, we explored the connections between AMF, fungivorous protists and nematodes that could partly reflect trophic interactions, and linked those to rhizosphere P dynamics and plant performance in a long-term manure application setting. Our results revealed that manure addition increased AMF biomass and the density of fungivorous nematodes, and tailored the community structures of AMF, fungivorous protists, and nematodes. We detected a higher abundance of AMF digested by the dominant fungivorous nematodes Aphelenchoides and Aphelenchus in high manure treatments compared to no manure and low manure treatments. Structural equation modeling combined with network analysis suggested that predation by fungivorous protists and nematodes stimulated AMF biomass and modified the AMF community composition. The mycorrhizal-fungivore interactions catalyzed AMF colonization and expression levels of the P transporter gene ZMPht1;6 in maize roots, which resulted in enhanced plant productivity. CONCLUSIONS: Our study highlights the importance of predation as a key element in shaping the composition and enhancing the biomass of AMF, leading to increased plant performance. As such, we clarify novel biological mechanism of the complex interactions between AMF, fungivorous protists, and nematodes in driving P absorption and plant performance. Video Abstract.

Effects of historical legacies on soil nematode communities are mediated by contemporary environmental conditions.

Both contemporary and historical factors are documented to be crucial in regulating species diversity and distribution. Soil fauna contribute substantially to global biodiversity and ecosystem functioning, while it is unclear whether and to what extent historical factors shape their diversity patterns. Here, we used soil nematodes as a model organism to test historical effects on soil fauna and to investigate the relative importance of climatic, soil, and historical factors. Based on nematode distribution data in 16 natural sites at a large scale (ranging from 22 to 40°N) in mainland China, we conducted elastic net regression model to test the effects of climatic (e.g., mean and seasonality of temperature/precipitation), soil (e.g., soil carbon, nitrogen, and pH), and historical (e.g., temperature/precipitation anomaly and the velocity of the change since the Last Glacial Maximum) variables on nematode genus richness and Shannon's diversity. Additionally, variation partitioning was used to determine the contribution of the three predictor sets to the explanation of both Jaccard and Bray-Curtis community dissimilarity. We found that climate generally explained more variations in both diversity and composition than soil and historical predictors in our samples. We also showed that although historical factors (e.g., temperature change velocity) were correlated with nematode diversity and composition, the pure effects of these historical factors were negligible. In other words, the historical effects were commonly represented by their interactions with current climatic and soil factors within our selected sites. Our results indicated that contemporary factors, especially climate, may outperform historical factors in regulating soil nematode diversity patterns at large scales.

Evaluating the effects of transgenic Bt rice cultivation on soil stability.

Insecticidal crystal (Cry) proteins produced by genetically modified rice that enter the soil via pollen dispersal, plant residues, and root exudation may disturb soil health. In the present study, we assessed the influences of transgenic Bt rice (i.e., HH1 with Cry1Ab/Cry1Ac) cultivation on the dynamics of soil carbon and nutrients under field conditions during 2013-2016. Transgenic treatments (transgenic Bt rice vs. its parental line (i.e., MH63) of non-Bt rice) have no consistently significant effects on soil property, including available nitrogen, available phosphorus, available potassium, total nitrogen, and total phosphorus, while apparent seasonal changes were observed. Besides, the variations of soil nutrients in the paddy field of transgenic Bt rice did not exceed their resistance capacities, except total organic carbon (TOC; RS (resistance) = 1.51) and total potassium (TK; RS = 2.62) in 2013 and TK (RS = 1.94) in 2014. However, the TOC and soil nutrient of TK in the paddy field of transgenic Bt rice have recovered to the pre-perturbation status after harvest (RL (resilience) = 1.01, F = 0.01, P = 0.91; RL = 0.98, F = 0.34, P = 0.58; RL = 0.99, F = 1.26, P = 0.29). Moreover, the paddy yield of transgenic Bt rice was consistently higher than that of its parental line of non-Bt rice. These results suggested that the cultivation of transgenic Bt rice has no adverse impact on soil stability in terms of soil carbon and nutrients and paddy yield.