Changes in the Abundance and Community Complexity of Soil Nematodes in Two Rice Cultivars Under Elevated Ozone

The atmospheric ozone concentrations have substantially increased in the surface layer over the past decades, and consequently exhibited a strong influence on soil microbial communities and functions. However, the effect of elevated ozone (eO₃) on the abundance, diversity, and structural complexity of soil nematode communities are elusive under different rice (Oryza sativa L.) cultivars. Here, the soil nematode community was investigated in two rice cultivars (Hybrid, Shanyou 63 vs. Japonica, Wuyujing 3) under open-top chambers (OTC) with control and eO₃ conditions. The results showed that the abundance of soil nematode community was altered by eO₃, but the responses were dependent on crop cultivars. The eO₃ decreased the total abundance and simplified the network complexity of the soil nematode community for both cultivars. However, eO₃ increased the abundance of c-p 4 in Shanyou 63, rather than Wuyujing 3, indicating that the hybrid rice cultivar could tradeoff the adverse impacts of eO₃ on the functional group of soil nematodes. Similarly, bacterivores belonging to K-strategy (c-p 4) increased under eO₃ in Shanyou 63, suggesting that the soil food web formed a bacteria-dominated channel under eO₃ for the hybrid rice cultivar. This study shed new light on the critical importance of rice cultivars in shaping the impacts of eO₃ on the soil micro-food web. Therefore, breeding and biotechnological approaches may become valuable pathways to improve soil health by shaping the community structures of the soil micro-food web in response to climate change in the future.

Introduction

Ozone is one of the most important atmospheric pollutants in the surface layer, seriously threatening the health of agricultural ecosystems (Tang et al., 2013). Since the industrial revolution, the concentration of ozone in the surface layer has been increasing and attracting extensive interest (Phillips et al., 2009). Previous studies have revealed that elevated ozone (eO₃) can damage plant leaf tissue and show detrimental effects on food production (Ashmore, 2005). For example, Zhu et al. (2011) suggested that a 25% increase in the ozone concentration caused a significant decrease of 20% in the grain yield of winter wheat in a fully open-air field condition. A meta-analysis aiming at investigating the impacts of eO₃ on crop yields showed that higher ozone concentrations (ca 31–50 ppb) resulted in the 8.9, 9.7, and 17.5% of declines in the yield production of barley, wheat, and rice, respectively (Feng and Kobayashi, 2009). Furthermore, eO₃ has been demonstrated to inhibit plant photosynthesis and negatively influence carbon allocation from plant to soil (Andersen, 2003; Ainsworth, 2017). Therefore, eO₃ can not only affect plant growth and crop production but also threaten the belowground microbial community and biodiversity (Agathokleous et al., 2020).

The soil micro-food web is indispensable for providing soil functions and services, such as food security, nutrient dynamic, biodiversity, and soil health (De Vries et al., 2013; Zhu et al., 2020). The soil nematode community serving as an essential component of the soil micro-food web is one of the most widely distributed invertebrates in global terrestrial ecosystems (Bongers and Ferris, 1999; Van Den Hoogen et al., 2019). Soil nematodes can simultaneously occupy different trophic levels and act as critical bio-indicators for predicting climate change in agroecosystems (Wang et al., 2021). However, large uncertainty exists regarding the influence of eO₃ on the soil nematode community. Because of the toxic effect of ozone, eO₃ shows an adverse impact on the abundance and diversity of soil nematodes in cropland fields (Bao et al., 2014; Zhang et al., 2021). A previous study found that the functional diversity of the soil nematode community was sensitive to eO₃, and the diversity index decreased whereas the dominance index increased under eO₃ in the open-top chambers (OTCs) for soybean plants (Bao et al., 2014). Li et al. (2015) showed that the legacy effect of eO₃ led to a significant decrease in the number of bacterivores and fungivores, while the total abundance of soil nematodes was increased by eO₃ in winter wheat fields (Li et al., 2015).

Plant cultivars can greatly contribute to the variable effects of eO₃ on soil nematode community, mainly due to the differences in the quantity and quality of plant resources (Li et al., 2012; Zhang et al., 2021). Several studies have verified that crop cultivars exhibited different sensitivity to eO₃ stress (Shi et al., 2009; Ainsworth, 2017). In a rice-wheat rotation system, wheat cultivars alleviated the residual effects of eO₃ on soil fungivores (Li et al., 2016). Li et al. (2012) found that the impact of eO₃ on soil nematode community was more sensitive in ozone-tolerant cultivars than in ozone-sensitive cultivars in winter wheat fields. Rice has been considered vulnerable to eO₃ compared with other major crop species (Feng and Kobayashi, 2009). The hybrid Indica cultivar, which can express the Bt on gene to resist the pests and disease invasion, has been widely cultivated to increase the yield (Tu et al., 2000). However, greater reductions in the leaf photosynthesis and grain yield through oxidative damage to cells have been found in the hybrid Indica cultivar exposed to eO₃ when compared with conventional cultivars (Xu et al., 2007; Pang et al., 2009). Moreover, the tropospheric ozone concentration could rise 40–60% by the end of the 21st century, resulting from a dramatic increase in the emissions of O₃ precursors (e.g., VOCs and NOx) (Meehl et al., 2007). To date, the response and underlying mechanism of soil nematode community to eO₃ remain unclear under different rice cultivars.

The objective of this study is therefore to identify the effects of eO₃ on the abundance, diversity, and structural complexity of soil nematode community under different rice cultivars through an OTCs experiment. The changes in the ecological environment and agricultural production in response to eO₃ were explored simultaneously. We hypothesized that (1) eO₃ would reduce the abundance and diversity of soil nematodes and simplify the structural complexity of the soil micro-food web, due to decreased plant biomass and soil food resources; and (2) soil nematode community would be more susceptible to eO₃ under the conventional cultivar compared with hybrid Indica cultivar, owing to the relatively low crop biomass or plant C resources input to the soil associated with the conventional cultivar.

Materials and Methods

Site Description

The experimental site is located at Jiangdu, in Jiangsu Province (119°43′E, 32°25′N), China. The site belongs to a subtropical monsoon anmate, with the average precipitation and temperature of 1131.3 mm and 16.2°C, respectively, from 2009 to 2018, respectively. The soil is classified as Gleyic Stagnic Anthrosol. The physicochemical properties of topsoil (0–15 cm) were as follows: soil pH 7.05, total carbon (TC) 6.94 g kg⁻¹, and total nitrogen (TN) 1.05 g kg⁻¹.

Experimental Design

This experiment was performed in six open-top chambers (OTCs) (i.e., three control conditions, CK vs. three elevated ozone, eO₃). The octagonal OTCs of 2.3 m in height and 4.8 m in diameter had aluminum alloy frames covered with walls of transparent tempered glass. In this study, an ozone generator (HY003, Chuangcheng Co., Jinan, China) was used to produce ozone, and then the produced gas was mixed and exchanged by a high-power fan (2200 w, CX125, Quanfeng Co., Shanghai, China). An ozone analyzer (Model 49i, Thermo Scientific, USA) was used to monitor and record the ozone concentrations at the plant canopy in the chambers. The flow of compressed oxygen and ozone concentrations was controlled by a mass flow meter according to actual and target ozone concentrations in the OTCs (Shang et al., 2021). The ozone fumigation began on 25 July 2020, with a daily duration between 8:00 a.m. and 6:00 p.m. The average ozone concentrations were 39.6 ± 1.9 ppb and 87.3 ± 0.9 ppb for CK and eO₃ chambers, respectively.

Rice cultivars of Shanyou 63 (hybrid Indica cultivar) and Wuyujing 3 (conventional japonica cultivar) have commonly been planted in this area. In this study, rice seeds were sown on 20th May 2020. Rice seedlings were transplanted into round pots (31 cm in height 22 cm in diameter, ca 15 kg soil) on 2nd July 2020, each pot had 3 hills, and then all pots were transferred to the chambers. In total, each treatment had three OTCs replicates, and each OTCs had 10 round pots, i.e., 5 round pots per cultivar (Shang et al., 2021).

Sample Collection and Analysis

After rice harvest, five-round pot soils were collected and then mixed as a sample. A total of 12 soil samples (2 ozone levels × 2 cultivars × 3 replicates) were collected. After the removal of all roots and plant residues, all samples were stored at 4°C until further analysis.

Soil microbial biomass carbon (MBC) and nitrogen (MBN) were extracted with 0.5 mol L⁻¹ K₂SO₄ after fumigation with ethanol-free chloroform and determined by a multi N/C Analyzer (TOC Analyzer, Germany), the differences in fumigated samples were used to calculate soil microbial biomass with coefficients of 0.38 (MBC) and 0.45 and (MBN), respectively. Soil dissolved organic carbon (DOC) was analyzed by a multi-N/C Analyzer (Jena TOC Analyzer, Germany) (Wu et al., 1990). Soil dissolved organic nitrogen (DON) was calculated based on soil contents of NH${}_4^{+}$-N, NO${}_3^{-}$-N, and TN. Soil NH${}_4^{+}$-N and NO${}_3^{-}$-N concentrations were extracted with 2 mol L⁻¹ KCL, and then were analyzed through a flow injection auto-analyzer (Skalar, The Netherlands). Soil pH was determined by 1/2.5 (soil/water ratio) with Mettler-Toledo pH (Lu, 2000).

Soil Nematode Analysis

Soil nematodes were extracted using a modified Baermann funnel as previously described (Wang et al., 2019). The abundances of soil nematodes were quantified based on a Motic microscope (40 × and 400 ×). The genera of soil nematodes were identified using 150 randomly selected individuals per soil sample by a Motic microscope. The trophic groups of the nematode community were assigned to bacterivores (BF), fungivores (FF), herbivores (PF), and omnivores-predators (OP), and values of 1–5 c-p referred to http://Nemaplex.ucdavis.edu. The abundance of soil nematode was converted to the individuals per 100 g of dry soil. The functional groups of soil nematodes were assigned to colonizers (c) and persisters (p) based on the c-p scale ranging from 1 to 5 (Bongers, 1990). The nematodes of c-p 1 resemble the r-strategy and occur in highly disturbed conditions. On the other hand, the nematodes of c-p 5 resemble the K-strategy, which is featured by high sensitivity and low fecundity (Neher and Darby, 2009). Soil nematodes were assigned to different functional groups combining the c-p class and the trophic groups of the soil nematode genus (Ferris et al., 2001). The Shannon Diversity Index (H') was used to calculate the diversity of soil nematode communities.

Statistical Analysis

All data were collated in Microsoft Excel 2016 and expressed as the mean of the three replicates plus/minus the standard deviation. All statistics were performed by R and IBM SPSS 22.0 software (SPSS Inc., Chicago, United States). The difference was considered statistically significant at the level of P < 0.05. General linear model analysis was carried out to explore the main influences of eO₃ and rice cultivars on the soil nematode community. Pearson correlation analysis was used to investigate the relationship between environmental factors and soil nematodes.

To determine the relationship between soil microbial community and soil environmental parameters, redundancy analysis (RDA) was performed by the ‘vegan' package in R. The non-metric multidimensional scaling analysis (NMDS) was carried out to reveal the distribution pattern and factors influencing soil nematode community under different treatments.

Structural equation modeling analysis (SEM) was carried out to explore the effects of eO₃ and rice cultivars on soil nematodes using AMOS (Amos Development Corporation, Chicago, IL, USA, Version 24.0). Before the SEM, we reduced the dimensionality of variables for soil abiotic factors (pH, NO${}_3^{-}$-N, NH${}_4^{+}$-N, DOC, and DON), biotic factors (plant, MBC, and MBN), and the dynamics of soil nematodes at different trophic groups through principal component analysis (PCA) (Chen et al., 2013).

The possible pairwise Spearman's rank correlations of nematode community were calculated by the ‘corr.test' function in the ‘psych' package in R. The connections represented significant (P ≤ 0.05) and strong (R² ≥ 0.6) correlations. A Gephi v.0.9.2 was carried out to visualize the correlation network of soil nematode community (Bastian et al., 2009).

Results

Variations of Edaphic Properties

Elevated ozone reduced soil DOC and DON by 17.5 and 40.7%, respectively, since there was no significant interaction between ozone and cultivar. However, eO₃ and rice cultivars did not significantly affect soil pH, NH${}_4^{+}$-N, or NO${}_3^{-}$-N content (Table 1). The significant interactive effect of eO₃ and rice cultivars was observed for MBN (P = 0.027): eO₃ only decreased MBN for Shanyou 63.

TABLE 1

Table 1. Changes in soil variables under elevated ozone conditions.

The Abundance and Diversity of Soil Nematodes

There were 38 genera of soil nematodes across all of the treatments. The herbivore nematodes (Helicotylenchus and Meloidogyne) were the dominant genera among all treatments (Table S1). Elevated ozone reduced the total nematode abundance by 28.4% (P = 0.052) (Table 2, Table S2), but had no effects on that of different trophic groups or the diversity of soil nematode community (Table S2, Figure S1). By contrast, eO₃ significantly decreased the abundances of c-p 1+2 groups (P = 0.067), rather than that of c-p 3 and c-p 5 in both cultivars (Table S2). Shanyou 63 reversed the negative impacts of eO₃ on the abundance of c-p 4, as indicated by the significant interaction between eO₃ and cultivars (P = 0.034). Similarly, eO₃ increased the abundance of bacterivores belonging to c-p 4 (i.e., BF4) for Shanyou 63, rather than Wuyujing 3 (Table 2). Elevated ozone had no significant impact on the relative abundance of nematode functional groups for both cultivars, while slightly increasing that of c-p 4 for Shanyou 63 (Figure S2).

TABLE 2

Table 2. Changes in the abundance of soil nematodes (individuals per 100 g dry soil) under elevated ozone conditions.

Soil Nematode Community

The effects of eO₃ and rice cultivars were evaluated by NMDS (Figure 1). Soil nematode community in CK was clearly separated from those under eO₃ in Shanyou 63, but not for Wuyujing 3, indicating that the ozone effect was largely dependent on rice cultivars (Figure 1). Network analysis was used to visualize the structural complexity of the soil nematode community, which showed that the networks were simplified by eO₃ (Figure 2). In the network, most nodes were derived from BF and PF, which showed decreasing node numbers under eO₃. Total nodes in the network were minimized by eO₃, and the connectivity in the eO₃ treatment was 21.1% lower than CK.

FIGURE 1

Figure 1. Non-multidimensional scaling (NMDS) analysis of the nematode community under elevated ozone.

FIGURE 2

Figure 2. The complexity and interactions of soil nematode community under elevated ozone. The connections represented strong (R² > 0.6) and significant (P < 0.05) correlations. The size of each node was proportional to the number of connections (degrees). Red lines represented significantly positive and purple lines represent significantly negative correlations. Different colors of nodes represented nematode trophic groups.

Relationships Between Soil Nematodes and Soil Properties

Soil pH was the most important parameter contributing to the changes in soil nematode community, followed by soil MBC and MBN (Figure 3). Correlation analysis revealed that the abundance of FF was positively correlated with soil pH, DON, and grain yield (Figure 4). The abundance of OP was positively correlated with soil NH${}_4^{+}$-N but negatively related to soil MBC. There was a significant positive correlation between the c-p 1+ 2 groups and MBN as well as grain yield. The abundance of the c-p 5 groups was significantly and negatively correlated with MBN.

FIGURE 3

Figure 3. Redundancy analysis (RDA) of the relationships between soil properties and soil nematode community under elevated ozone.

FIGURE 4

Figure 4. Correlations between soil properties, grain yield and soil nematode community under elevated ozone. DOC, dissolved organic carbon; DON, dissolved organic nitrogen; MBC, microbial biomass carbon; MBN, microbial biomass nitrogen; TN, total nematodes; BF, bacterivores; FF, fungivores; PF, plant parasites; OP, omnivorous-predatory; H′, Shannon-Wiener diversity index. * The asterisk represents a significant correlation at the 0.05 probability level.

The Chi-square (χ²) test of SEM fit (χ² = 4.635, P = 0.704, df = 7), CFI (comparative fit index) = 1.000, and RMSEA (root square mean error of approximation) < 0.001, indicated that the SEM was of excellent fit (Figure 5). Elevated ozone positively correlated with soil abiotic factors (covariance coefficient = 0.573), but negatively correlated with biotic factors (covariance coefficient = −0.671). Rice cultivars only negatively correlated with biotic factors (covariance coefficient = −0.663). Therefore, eO₃ and cultivars explained 40.5 and 56.7% of the total variances in abiotic factors, and biotic factors, respectively. Elevated ozone, cultivars, soil abiotic factors, and biotic factors pathways explained 26.3, 44.5, and 47.6% of the total variances in bacterivores, fungivores, and herbivores, respectively. Interestingly, herbivores were directly related to eO₃ (covariance coefficient = −0.908) and statistically associated with soil abiotic factors (covariance coefficient = 0.554). The model explained 53.5% of the variance in omnivores-predators, which positively correlated with bacterivores (covariance coefficient = 0.521).

FIGURE 5

Figure 5. Structural equation models (SEM) analysis of elevated ozone and rice cultivars on soil nematodes. The results of the model fitting: Chi-square (χ²) = 4.635, P = 0.704, df = 7, CFI (comparative fit index) = 1, RMSEA (root square mean error of approximation) < 0.001. Red and purple arrows showed positive and negative effects, while solid arrows and dashed arrows showed the directions and effects that were significant (P < 0.05) and non-significant (P > 0.05), respectively.

Discussion

Elevated Ozone Reduced the Abundance and Community Complexity of Soil Nematodes

We found that eO₃ reduced the total abundance of soil nematodes (Table 2). This result was consistent with several previous studies, which found that the abundance of soil nematodes was significantly decreased by eO₃ (Bao et al., 2014; Li et al., 2015). This could be explained by the fact that ozone is an effective biocide, and eO₃ can inhibit the biomass and reproduction of soil bacteria and fungi (Kim et al., 1999; Kanerva et al., 2008). As a consequence, the reduction in bacterial and fungal biomass would affect the soil food web and the abundance of soil nematodes (Zhang et al., 2021).

The SEMs results showed that eO₃ directly and negatively correlated with herbivores (Figure 5), indicating that eO₃ directly inhibited the abundance of soil nematodes. Consistent with Qiu et al. (2009), it was found that the root-knot nematodes were reduced by 24 and 68%, and the free-living nematodes were decreased by 19 and 52% at the ozone concentration of 50 and 250 kg hm⁻², respectively. Previous studies have confirmed that ozone is an active oxidant with biocidal characteristics, which could directly kill soil nematodes (Mitsugi et al., 2017; Veronico et al., 2017). Therefore, the survival rate of natural nematodes can be directly minimized as the application rate of ozone increases (Mitsugi et al., 2017).

In addition, the impact of eO₃ on plant growth could indirectly inhibit the activity and reproduction of the soil nematode community (Wang et al., 2021; Zhang et al., 2021). Ozone as a phytotoxic compound, can inhibit plant growth and decrease root exudates and carbon allocation (Hu et al., 2018), resulting in reduced resource availability for the soil microbial community. It has been acknowledged that ozone can decrease the transportation of carbohydrates to rhizodeposition, which inhibited the reproduction of arbuscular mycorrhizal fungi, and thus provided fewer nutrients for soil nematodes (Wang et al., 2011; Zhang et al., 2021). Li et al. (2012) found that DOC content was decreased by eO₃, and the effect of eO₃ on the soil food web was associated with soil carbon and nitrogen contents. In the present study, we found that eO₃ decreased the contents of DOC and DON (Table 1), which correlated to the abundances of FF and PF (Figure 3). Furthermore, the crop yields under eO₃ were observed to be reduced by 37.8% in the present study (Figure S3). The detrimental effect of eO₃ on resource partitioning between above- and below-ground may lead to an indirect influence on the stability of the soil nematode community.

Consistent with our first hypothesis, the structural complexity of the soil nematode community was simplified by eO₃ (Figure 2), indicating that eO₃ reduced the complexity and interactions of the soil food web. As mentioned above, the abundance of soil nematode was decreased under eO₃, which was partially ascribed to the reduction in crop production (Figure S3). The variation in plant inputs to the soil under eO₃ can propagate through carbon resources to alter the complexity and interactions of the soil food web (Chung et al., 2006). Moreover, the declines in root biomass induced by eO₃ meant fewer resources available for the soil nematode communities (Bao et al., 2014). Accordingly, eO₃ had a direct and negative impact on soil nutrients such as DOC and DON (Table 1) and subsequently propagated into soil nematode communities, especially herbivores (Figure 5). Given that the network complexity of the feeding links was significantly correlated with the stability and stress tolerance of the soil food web (Beckerman et al., 2006), our findings thus highlighted that eO₃ went against the stability and interactions of the soil micro-food web.

Rice Cultivars Altered the Effects of eO₃ on Functional Groups of Soil Nematodes

Soil nematode functional groups were sensitive to eO₃ and cultivar effects (Li et al., 2012). Interestingly, the abundance of bacterivores belonging to K-strategies (BF4) was increased under eO₃ (Table 2). The increase in the abundance of bacterivores within K-strategies (BF4) groups may be caused by the stability of the bacterial feeding channel under the stress conditions (De Vries et al., 2013). It has been recognized that the bacterial feeding channel (bacteria and bacterivores) is more resilient than other trophic groups of the slowly growing biotas under climate change (De Vries et al., 2013). Our result showed that bacterivores were directly and statistically associated with omnivores-predators (Figure 5), indicating that the soil micro-food web formed a bacterial-dominated channel under eO₃ in paddy soil.

Our results suggested that rice cultivars altered the adverse effects of eO₃ on soil nematodes. In favor of our hypothesis, the cultivar of Shanyou 63 reversed the negative impact of eO₃ on the abundance of soil nematodes belonging to c-p 4 (Table 2). This is consistent with a previous study by Li et al. (2012), who observed that eO₃ reduced the abundance of bacterivores in the ozone-sensitive wheat but increased it in the ozone-tolerant wheat. Rice has been considered vulnerable to eO₃ (Feng and Kobayashi, 2009). The response of rice yield to eO₃ showed varietal difference (Pang et al., 2009), which indirectly altered the effect of eO₃ on soil nematodes through changes in the input of nutrient resources. Therefore, the negative effect of eO₃ on the soil nematode community could be weakened under the hybrid Indica cultivar in paddy soil.

Crop cultivars played an essential role in determining the response of soil nematode community to eO₃ (Li et al., 2012; Feng et al., 2017). The community composition of soil nematodes responded differently to specific plants under eO₃ (De Deyn et al., 2004). In this study, rice cultivars significantly influenced the nematode abundance of c-p 1+2. In particular, the nematode abundance of c-p 1+2 groups in Wuyujing 3 was decreased by 44.7% compared with Shanyou 63 (Table 2). Plants can influence soil biota through changes in the inputs of root exudates and plant litter (Zhang et al., 2019). Previous studies have revealed that the aboveground biomass, yield, and harvest index of Shanyou 63 were higher than that of Wuyujin 3 (Xu et al., 2007; Pang et al., 2009). Data from the same experiment showed that the grain yield of Shanyou 63 was 31.9% higher than that of Wuyujing 3 (Figure S3). The changes in specific soil nematode groups could therefore depend on plant inputs of soil food resources and energy, and the drop in crop biomass could directly result in the decreased abundance of soil nematodes under climate change.

In this study, eO₃ increased the abundances of c-p 4 groups and bacterivores belonging to K-strategies in Shanyou 63, rather than Wuyujing 3, indicating that the hybrid Indica cultivar mitigated the negative impact of eO₃ on the functional groups of soil nematode community. This finding provided important insights for the adaptation of soil food web and ecosystem functions to climate change in the agricultural ecosystem. Considering that the stress derived from eO₃ on soil nematodes and food production continues to increase in near future (Ainsworth, 2017). It is an essential task to develop new rice cultivars and benefit from the positive aspects of cultivars with the capacity for restraining the adverse effects of climate change (i.e., eO₃) on soil health and ecosystem functions (Ortiz et al., 2008; Wang et al., 2019).

Conclusion

This study demonstrated that eO₃ significantly reduced the total abundance and simplified the structural complexity of the soil nematode community. However, the hybrid Indica cultivar could tradeoff the negative impact of eO₃ on the abundance of bacterivores belonging to K-strategies, suggesting that the hybrid Indica cultivar may alleviate the threat of eO₃ to the soil micro–food web. This study revealed that the soil nematode community was altered by eO₃, but the variations of nematode functional groups were dependent on rice cultivars. Our study highlighted that the breeding and biotechnological approaches could be valuable approaches to tradeoff the adverse impacts of eO₃ on the soil nematode community, and improve the potential feedback of soil micro-food web in agricultural ecosystems.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

JW and YT analyzed the data. YT and YS performed the laboratory work. JW and XS drafted the manuscript. XS and GZ contributed ideas to the study and revised the first draft. All authors contributed to the article and approved the submitted version.

Funding

The research was supported by the National Natural Science Foundation of China (Grant Nos. 42077209, 31901165, 32071631, and 41907022), the Sino-German Mobility Programme (M-0105), the Open Project of Jiangsu Key Laboratory of Crop Genetics and Physiology (YCSL202004), and Natural Science Foundation of Fujian Province, China (Grant No. 2020J01186), and the Open Fund of Key Laboratory of Agrometeorology of Jiangsu Province (JKLAM2001). The design of the study and collection of data was supported by the Startup Foundation for Introducing Talent of Nanjing University of Information Science and Technology (003035), Nanjing, China.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.916875/full#supplementary-material

References

Agathokleous, E., Feng, Z., Oksanen, E., Sicard, P., Wang, Q., Saitanis, C. J., et al. (2020). Ozone affects plant, insect, and soil microbial communities: A threat to terrestrial ecosystems and biodiversity. Sci. Adv. 6, eabc1176. doi: 10.1126./sciadv.abc1176

PubMed Abstract | CrossRef Full Text | Google Scholar

Ainsworth, E. A. (2017). Understanding and improving global crop response to ozone pollution. Plant J. 90, 886–897. doi: 10.1111/tpj.13298

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersen, C. P. (2003). Source–sink balance and carbon allocation below ground in plants exposed to ozone. New Phytol. 157, 213–228. doi: 10.1046/j.1469-8137.2003.00674.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashmore, M. R. (2005). Assessing the future global impacts of ozone on vegetation. Plant Cell Environ. 28, 949–964. doi: 10.1111/j.1365-3040.2005.01341.x

CrossRef Full Text | Google Scholar

Bao, X., Li, Q., Hua, J., Zhao, T., and Liang, W. (2014). Interactive effects of elevated ozone and UV-B radiation on soil nematode diversity. Ecotoxicology. 23, 11–20. doi: 10.1007/s10646-013-1146-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: an open source software for exploring and manipulating networks,” in Third International AAAI Conference on Weblogs and Social Media.

Google Scholar

Beckerman, A. P., Petchey, O. L., and Warren, P. H. (2006). Foraging biology predicts food web complexity. Proc. Natl. Acad. Sci. 103, 13745. doi: 10.1073/pnas.0603039103

PubMed Abstract | CrossRef Full Text | Google Scholar

Bongers, T. (1990). The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 83, 14e19. doi: 10.1007/BF00324627

PubMed Abstract | CrossRef Full Text | Google Scholar

Bongers, T., and Ferris, H. (1999). Nematode community structureas a bioindicator in environmental monitoring. Trends Ecol. Evol. 14, 224–228. doi: 10.1016/S0169-5347(98)01583-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D., Lan, Z., Bai, X., Grace, J. B., Bai, Y., and van der Heijden, M. (2013). Evidence that acidification-induced declines in plant diversity and productivity are mediated by changes in below-ground communities and soil properties in a semi-arid steppe. J. Ecol. 101, 1322–1334. doi: 10.1111/1365-2745.12119

CrossRef Full Text | Google Scholar

Chung, H., Zak, D. R., and Lilleskov, E. A. (2006). Fungal community composition and metabolism under elevated CO₂ and O₃. Oecologia 147, 143–154. doi: 10.1007/s00442-005-0249-3

PubMed Abstract | CrossRef Full Text | Google Scholar

De Deyn, G. B., Raaijmakers, C. E., Van Ruijven, J., Berendse, F., and Van Der Putten, W. H. (2004). Plant species identity and diversity effects on different trophic levels of nematodes in the soil food web. Oikos 106, 576–586. doi: 10.1111/j.0030-2004.13265.x

CrossRef Full Text | Google Scholar

De Vries, F. T., Thebault, E., Liiri, M., Birkhofer, K., Tsiafouli, M. A., Bjornlund, L., et al. (2013). Soil food web properties explain ecosystem services across European land use systems. Proc. Natl. Acad. Sci. 110, 14296–14301. doi: 10.1073/pnas.1305198110

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, Z., and Kobayashi, K. (2009). Assessing the impacts of current and future concentrations of surface ozone on crop yield with meta-analysis. Atmos. Environ. 43, 1510–1519. doi: 10.1016/j.atmosenv.11.033

CrossRef Full Text | Google Scholar

Feng, Z., Tang, H., and Kobayashi, K. (2017). “Effects of ozone on crops in China,” in Air pollution impacts on plants in East Asia (Tokyo: Springer), 175–194. doi: 10.1007/978-4-431-56438-6_12

CrossRef Full Text | Google Scholar

Ferris, H., Bongers, T., and De Goede, R. (2001). A framework for soil food web diagnostics: extension of the nematode faunal analysis concept. Appl. Soil Ecol. 18, 13e29. doi: 10.1016/S0929-1393(01)00152-4

CrossRef Full Text | Google Scholar

Hu, E., Yuan, Z., Zhang, H., Zhang, W., Wang, X., Jones, S. B., et al. (2018). Impact of elevated tropospheric ozone on soil C, N and microbial dynamics of winter wheat. Agric. Ecosyst. Environ. 253, 166–176. doi: 10.1016/j.agee.11.010

CrossRef Full Text | Google Scholar

Kanerva, T., Palojärvi, A., Rämö, K., and Manninen, S. (2008). Changes in soil microbial community structure under elevated tropospheric O₃ and CO₂. Soil Biol. Biochem. 40, 2502–2510. doi: 10.1016/j.soilbio.06.007

CrossRef Full Text | Google Scholar

Kim, J. G., Yousef, A. E., and Dave, S. (1999). Application of ozone for enhancing the microbiological safety and quality of foods: a review. J. Food Protect. 62, 1071–1087. doi: 10.4315/0362-028X-62.9.1071

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Bao, X., Lu, C., Zhang, X., Zhu, J., Jiang, Y., et al. (2012). Soil microbial food web responses to free-air ozone enrichment can depend on the ozone-tolerance of wheat cultivars. Soil Biol. Biochem. 47, 27–35. doi: 10.1016/j.soilbio.12.012

CrossRef Full Text | Google Scholar

Li, Q., Yang, Y., Bao, X., Liu, F., Liang, W., and Zhu, J. (2015), Legacy effects of elevated ozone on soil biota plant growth. Soil Biol. Biochem. 91, 50–57. 10.1016/j.soilbio.08.029.

Google Scholar

Li, Q., Yang, Y., Bao, X., Zhu, J., Liang, W., and Bezemer, T. M. (2016). Cultivar specific plant-soil feedback overrules soil legacy effects of elevated ozone in a rice-wheat rotation system. Agric. Ecosyst. Environ. 232, 85–92. doi: 10.1016/j.agee.07.025

CrossRef Full Text | Google Scholar

Lu, R. (2000). Methods of Soil and Agro-Chemistry Analysis. Beijing: China Agricultural Science and Technology Press, 62–141.

PubMed Abstract | Google Scholar

Meehl, G. A., Stocker, T. F., Collins, W. D., Friedlingstein, P., Gaye, A. T., Gregory, J. M., et al. (2007). Global Climate Projections. Cambridge, UK/New York, NY: Cambridge University Press.

Google Scholar

Mitsugi, F., Abiru, T., Ikegami, T., Ebihara, K., and Nagahama, K. (2017). Treatment of nematode in soil using surface barrier discharge ozone generator. IEEE T. Plasma Sci. 45, 3076–3081. doi: 10.1109/TPS.2017.2708706

CrossRef Full Text | Google Scholar

Neher, D. A., and Darby, B. J. (2009). “General community indices that can be used for analysis of nematode assemblages,” in Nematodes as Environmental Indicators, 107e123. doi: 10.1079./9781845933852.0107

CrossRef Full Text | Google Scholar

Ortiz, R., Sayre, K. D., Govaerts, B., Gupta, R., Subbarao, G., Ban, T., et al. (2008). Climate change: can wheat beatthe heat? Agric. Ecosyst. Environ. 126, 46–58. doi: 10.1016/j.agee.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, J., Kobayashi, K., and Zhu, J. (2009). Yield and photosynthetic characteristics of flag leaves in Chinese rice (Oryza sativa L.) varieties subjected to free-air release of ozone. Agric. Ecosyst. Environ. 132, 203–211. doi: 10.1016/j.agee.03.012

CrossRef Full Text | Google Scholar

Phillips, D. L., Johnson, M. G., Tingey, D. T., and Storm, M. J. (2009). Elevated CO₂ and O₃ effects on fine-root survivorship in ponderosa pine mesocosms. Oecologia 160, 827–837. doi: 10.1007/s00442-009-1339-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, J. J., Westerdahl, B. B., and Pryor, A. (2009). Reduction of root-knot nematode, Meloidogyne javanica, and ozone mass transfer in soil treated with ozone. J. Nematol. 41, 241.

PubMed Abstract | Google Scholar

Shang, B., Fu, R., Agathokleous, E., Dai, L, Zhang, G, Wu, R., et al. (2021). Ethylenediurea offers moderate protection against ozone-induced rice yield loss under high ozone pollution. Sci. The Total Environ. 806, 151341. doi: 10.1016/j.scitotenv.2021.151341

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, G., Yang, L., Wang, Y., Kobayashi, K., Zhu, J., Tang, H., et al. (2009). Impact of elevated ozone concentration on yield of four Chinese rice cultivars under fully open-air field conditions. Agric. Ecosyst. Environ. 131, 178–184. doi: 10.1016/j.agee.01.009

CrossRef Full Text | Google Scholar

Tang, H., Takigawa, M., Liu, G., Zhu, J., and Kobayashi, K. (2013). A projection of ozone-induced wheat production loss in China and India for the years 2000 and 2020 with exposure-based and flux-based approaches. Glob. Change Biol. 19, 2739–2752. doi: 10.1111/gcb.12252

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, J., Zhang, G., Datta, K., Xu, C., He, Y., Zhang, Q., et al. (2000). Field performance of transgenic elite commercial hybrid rice expressing Bacillus thuringiensis δ-endotoxin. Nat. Biotechnol. 18, 1101–1104. doi: 10.1038/80310

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Den Hoogen, J., Geisen, S., Routh, D., Ferris, H., Traunspurger, W., Wardle, D. A., et al. (2019). Soil nematode abundance and functional group composition at a global scale. Nature 572, 194–198. doi: 10.1038/s41586-019-1418-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Veronico, P., Paciolla, C., Sasanelli, N., De Leonardis, S., and Melillo, M. T. (2017). Ozonated water reduces susceptibility in tomato plants to Meloidogyne incognita by the modulation of the antioxidant system. Mol. Plant Pathol. 18, 529–539. doi: 10.1111/mpp.12413

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Li, M., Zhang, X., Liu, X., Li, L., Shi, X., et al. (2019). Changes in soil nematode abundance and composition under elevated [CO₂] and canopy warming in a rice paddy field. Plant Soil 445, 425–437. doi: 10.1007/s11104-019-04330-4

CrossRef Full Text | Google Scholar

Wang, J., Shi, X., Li, L., and Zhang, X. (2021). Changes in soil nematodes in rhizosphere and non-rhizosphere soils following combined elevated [CO₂] and canopy warming in a winter wheat field. Geoderma. 386, 114907. doi: 10.1016/j.geoderma.2020.114907

CrossRef Full Text | Google Scholar

Wang, S., Feng, Z., Wang, X., and Gong, W. (2011). Arbuscular mycorrhizal fungi alter the response of growth and nutrient uptake of snap bean (Phaseolus vulgaris L.) to O₃. J. Environ. Sci. 23, 968–974. doi: 10.1016/S1001-0742(10)60503-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Joergensen, R. G., Pommerening, B., Chaussod, R., and Brookes, P. C. (1990). Measurement of soil microbial biomass C by fumigation-extraction—an automated procedure. Soil Biol. Biochem. 22, 1167–1169. doi: 10.1016/0038-0717(90)90046-3

CrossRef Full Text | Google Scholar

Xu, G., Zhang, J., Lam, H. M., Wang, Z., and Yang, J. (2007). Hormonal changes are related to the poor grain filling in the inferior spikelets of rice cultivated under non-flooded and mulched condition. Field Crops Res. 101, 53–61. doi: 10.1016/j.fcr.09.008

CrossRef Full Text | Google Scholar

Zhang, P., Li, B., Wu, J., and Hu, S. (2019). Invasive plants differentially affect soil biota through litter and rhizosphere pathways: a meta-analysis. Ecol. Lett. 22, 200–210. doi: 10.1111/ele.13181

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wang, H., Wang, Y., Zhang, X., Zhao, T., and Mahamood, M. (2021). Organic input practice alleviates the negative impacts of elevated ozone on soil microfood-web. J. Clean. Prod. 290, 125773. doi: 10.1016/j.jclepro.2020.125773

CrossRef Full Text | Google Scholar

Zhu, D., Ding, J., Yin, Y., Ke, X., O'Connor, P., and Zhu, Y. G. (2020). Effects of earthworms on the microbiomes and antibiotic resistomes of detritus fauna and phyllospheres. Environ. Sci. Technol. 54, 6000–6008. doi: 10.1021/acs.est.9b04500

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Feng, Z., Sun, T., Liu, X., Tang, H., Zhu, J., et al. (2011). Effects of elevated ozone concentration on yield of four Chinese cultivars of winter wheat under fully open-air field conditions. Glob. Change Biol. 17, 2697–2706. doi: 10.1111/j.1365-2011.02400.x

CrossRef Full Text | Google Scholar