Biodiversity plays pivotal roles in ecosystem functioning and provision of ecosystem services that are crucial to human well-being. These services include providing food and water; managing floods, pests, and diseases; and supporting photosynthesis, nutrient cycling, soil formation, and crop pollination that sustain all other services (Millennium Ecosystem Assessment (MEA), 2003). Unfortunately, modern human civilization occurs at the expense of biodiversity. Land transformation is the principal driving force for biodiversity loss. Human activity has extensively transformed the land surface by agricultural intensification and urbanization (Vitousek et al., 1997). Urbanization and agricultural practices such as burning, tillage, fertilizer applications, and mono-cultural cropping practices affect below-ground biodiversity and its functions including decomposition, nutrient cycling, bioremediation, and pest and disease regulation (Giller et al., 1997). Despite its diverse benefits, biodiversity in soils is understudied compared to above-ground biodiversity.
Soil is a dynamic system in which organisms interact with each other and form complex food webs (Hunt and Wall, 2002). Nematodes are at the central place in the soil food web because they represent multiple trophic levels including primary, secondary, and tertiary consumer levels (Yeates et al., 1993). The structure of a nematode community provides good information on the condition of the soil food web since nematodes are specific in their food sources and are most abundant in all habitats where decomposition occurs (Bongers and Bongers, 1998). Yeates et al. (1993) assigned nematodes to different trophic groups such as bacterivores, fungivores, herbivores, predators, and omnivores based on their feeding habits. Bacterivores, fungivores, and herbivores are considered as nematode trophic groups in the lower hierarchy of the soil food web and predators and omnivores are considered as nematode trophic groups in the higher hierarchy of the soil food web (Yodzis, 2001). Nematode trophic interactions contribute to regulating nutrient dynamics in soil. Bacterivores and fungivores promote N and C mineralization by feeding on decomposing bacterial and fungal biomass. Nematode trophic groups in the higher hierarchy of the soil food web maintain ecological balance between decomposition and mineralization by regulating bacterivores and fungivores (Ingham et al., 1985). In addition, predators act as biocontrol agents by feeding on plant feeding nematodes (Bilgrami and Brey, 2005). Bongers (1990) developed a colonizer-persister (c-p) scale for nematodes by allocating the nematode taxa to one of five c-p groups ranging from colonizers (c) with a c-p value 1 to persisters (p) with a c-p value 5 through intermediate values based on their life history characteristics and survival strategies. Nematodes with small size, short life span, and high fecundity are assigned to c-p 1, whereas those with the longest-lived nematodes, low fecundity, and slow in development are placed in c-p 5 (Bongers, 1990). Many useful indices for nematode faunal analysis have been developed based on trophic groups and c-p scale. Consequently, nematodes can be used as indicators of structure and function of soil food webs and overall ecosystem conditions (Ferris et al., 2001).
A plethora of published literature exists on how different ecosystems affect the abundance (number of nematodes) and richness (number of taxa) of nematodes. However, there is no single consensus about the pattern of nematode abundance and richness in different ecosystems across the published literature. Some authors have reported that richness is high in forest ecosystems and abundance is high in agricultural ecosystems (Yeates and Bongers, 1999; Ferris et al., 2001; Yeates, 2007; Cardoso et al., 2015) but others have stated the converse (Neher et al., 2005; Briar et al., 2007; Darby et al., 2007; Kimenju et al., 2009). The existence of a large body of literature with diverse results creates the need to synthesize quantitative summaries in order to draw general conclusions across studies and test key hypotheses regarding patterns and processes governing soil biodiversity. Meta-analysis is a tractable and powerful statistical tool developed to generate a quantitative summary of all the published literature and draw conclusions across multiple studies (Arnqvist and Wooster, 1995). Therefore, meta-analysis was chosen to address this issue.
The specific objective of this study was to assess the influence of agricultural intensification and urbanization on nematode richness and abundance compared to forest and grassland ecosystems through meta-analysis of published literature on a global scale. The richness and abundance of nematodes were compared using different moderator levels or explanatory variables. We hypothesized that overall richness, overall abundance, and richness and abundance of nematodes of each trophic group and c-p guild are greater in forest and natural grassland (NGL) ecosystems compared to urban, agriculture and disturbed grassland (DGL) ecosystems.
The Web of Science core database was systematically searched for relevant publications on October 7, 2016, with the following combination of search terms: (‘nematode communities’ or ‘soil nematodes’ or ‘nematode diversity’ or ‘nematode abundance’ or ‘nematode biodiversity’) and (‘grassland’ or ‘forest’ or ‘agriculture’ or ‘prairie’ or ‘urban’), which resulted in 1,613 articles. Criteria for including an article in the analysis were: studies were conducted in forest, grassland, urban, or agriculture ecosystems; studies identified nematodes to family or genus level; studies reported mean abundance or richness expressed per grams or cm3 of soil; soil samples were collected from natural conditions; and studies reported sample size. Criteria for excluding an article were: studies conducted in controlled conditions like microcosms, mesocosms, pots, or greenhouses; studies expressing abundance of nematodes as relative abundance instead of absolute abundance; and studies reporting data for total free-living nematodes instead of each trophic group. Among the 1,613 articles, 598 relevant articles that contained data on richness and abundance of nematodes in different ecosystems were selected by examining titles and abstracts. Among the 598 articles, 111 articles (Supplementary Material) met the inclusion criteria and were selected for data extraction. Among the 111 articles, 91 expressed data in grams and 20 expressed data in cm3. The first 200 articles from a Google Scholar search were examined using the above search terms, which did not produce additional articles. A spreadsheet was constructed by extracting data from each article on authors, title, year of publication, unit of soil, richness and abundance of nematodes of each trophic group and each c-p guild, overall richness and overall abundance of nematodes, treatment, sample size, and type of ecosystem. Overall richness and overall abundance of nematodes were calculated by adding the number of genera/families and abundance of nematodes of either all trophic groups or all c-p guilds, respectively. Richness and abundance of nematodes under each trophic group and each c-p guild were calculated by adding the number of genera/families and abundance of nematodes corresponding to each guild and each trophic group, respectively. If there was more than one treatment in an article, they were considered as distinct studies in the meta-analysis. For example, there were two treatments, conventional-conservation tillage and organic-conservation tillage in Sánchez-Moreno et al. (2009), these two treatments were considered as two distinct studies. Based on these criteria, a total of 667 studies were subjected for meta-analysis of which 449 studies conducted in agriculture, 28 conducted in DGL, 74 conducted in forest, 36 conducted in NGL, and 80 conducted in urban ecosystems. Soil units in nematode studies are typically expressed as grams (Briar et al., 2007) or in cm3 (Wang et al., 2006). Therefore, the richness and abundance of nematodes expressed per 100 g of soil and 100 cm3 of soil were analyzed separately. Richness and abundance of nematodes per 100 g of soil were compared across all five ecosystems. However, the data expressed per 100 cm3 of soil were compared across only four ecosystems as no urban ecosystem studies using 100 cm3 were available. Abundance of nematodes that was not expressed per 100 g or cm3 of soil was converted to 100 g or cm3 of soil. However, richness of nematodes was not converted because increase in richness cannot be assessed with increase in the quantity of soil.
Effect size typically represents the strength of the relationship between two variables or two groups (treatment and control) but can also refer to the estimate of a single group or value such as richness or abundance of each study (Borenstein et al., 2009). Summary effect size is defined as weighted mean of richness or abundance of all studies in each ecosystem. Meta-analyses were conducted to compare the summary effect sizes of overall richness and overall abundance of nematodes and nematodes of each trophic group and each c-p guild per soil weight and volume basis among different ecosystems such as forest, NGL, DGL, agriculture, and urban ecosystems. Overall richness and overall abundance of nematodes per grams and per cm3 of soil were considered as four main effect sizes; richness and abundance of nematodes per grams and per cm3 of soil in each trophic group and each c-p guild were considered as subgroup effect sizes.
The types of ecosystems, forest, NGL, DGL, agriculture, and urban, were considered as moderator levels. These five ecosystems were assumed to have different regimes of disturbance where forest and NGL are considered less disturbed, whereas agriculture and urban ecosystems are considered highly disturbed from continuous human intervention. The moderator was chosen to determine the influence of disturbance on soil health.
The procedures and terminology of Borenstein et al. (2009) were followed in this analysis. Comprehensive meta-analysis (CMA) software was used to estimate effects of different levels of moderator on nematodes based on their confidence intervals,
A sensitivity analysis was conducted to assess the stability and consistency of the summary effects. The summary effect was recalculated by removing one study at a time. This measures how sensitive the results are to any one study. The potential presence of publication bias was tested using the Begg and Mazumdar rank (Kendall) correlation test and graphically by examining summary effect sizes vs their standard errors in funnel plots (Begg and Mazumdar, 1994; Borenstein et al., 2009).
A total of 44 summary effect sizes were tested in the meta-analysis performed, of which 40 summary effect sizes were significantly heterogeneous (
Heterogeneity statistics for the summary effect sizes per 100 g and per 100 cm3 of soil.
100 g | 100 cm3 | |||||
---|---|---|---|---|---|---|
Summary effect | ||||||
Overall richness | 740.37 | 0.000 | 23.37 | 49.75 | 0.103 | 12.42 |
Overall abundance | 525.42 | 0.007 | 2.67 | 320.94 | 0.000 | 10.28 |
Richness of c-p 1 | 347.88 | 0.000 | 8.50 | 129.66 | 0.000 | 66.24 |
Richness of c-p 2 | 486.97 | 0.000 | 15.43 | 79.06 | 0.000 | 34.16 |
Richness of c-p 3 | 453.61 | 0.000 | 32.25 | 147.73 | 0.000 | 73.73 |
Richness of c-p 4 | 520.05 | 0.000 | 29.18 | 42.39 | 0.357 | 7.62 |
Richness of c-p 5 | 390.74 | 0.001 | 4.56 | 54.84 | 0.025 | 17.00 |
Abundance of c-p 1 | 553.15 | 0.009 | 2.43 | 330.58 | 0.000 | 6.07 |
Abundance of c-p 2 | 422.77 | 0.028 | 2.57 | 186.77 | 0.000 | 27.30 |
Abundance of c-p 3 | 1,299.77 | 0.000 | 2.07 | 224.18 | 0.000 | 43.86 |
Abundance of c-p 4 | 609.75 | 0.000 | 13.95 | 70.48 | 0.088 | 9.27 |
Abundance of c-p 5 | 730.70 | 0.000 | 9.33 | 159.32 | 0.000 | 14.05 |
Richness of bacterivores | 584.92 | 0.000 | 17.00 | 76.29 | 0.000 | 29.57 |
Richness of fungivores | 392.01 | 0.000 | 18.47 | 105.16 | 0.000 | 57.06 |
Richness of herbivores | 358.48 | 0.000 | 15.42 | 69.50 | 0.000 | 37.84 |
Richness of predators | 267.55 | 0.000 | 18.34 | 50.88 | 0.061 | 14.51 |
Richness of omnivores | 446.01 | 0.000 | 18.48 | 48.12 | 0.135 | 11.56 |
Abundance of bacterivores | 519.91 | 0.001 | 3.80 | 396.91 | 0.000 | 9.81 |
Abundance of fungivores | 645.08 | 0.034 | 1.61 | 357.16 | 0.000 | 17.18 |
Abundance of herbivores | 762.77 | 0.015 | 1.62 | 430.30 | 0.001 | 3.92 |
Abundance of predators | 768.10 | 0.000 | 6.72 | 144.93 | 0.000 | 18.25 |
Abundance of omnivores | 747.91 | 0.000 | 11.09 | 344.69 | 0.000 | 12.77 |
Sensitivity analysis indicates the contribution of each study to the summary effect, which is measured by the change in the summary effect in its absence. The summary effect size of overall abundance per 100 g of soil was most affected by the removal of treatment B4 at Bohemia in the study conducted by Čermák et al. (2011). This study reduced the summary effect size from 1,208.00 to 1,186.23 (Supplementary Material, Table 1). Similarly, the summary effect size of overall richness per 100 g of soil was most influenced by the removal of Renčo and Baležentiené (2015), grassland (control) treatment, reducing the summary effect size from 27.35 to 27.21 (Supplementary Material, Table 2). The summary effect size of overall abundance per 100 cm3 was most affected by the removal of the Bulluck et al. (2002), cotton-gin trash (harvest) treatment. This study reduced the summary effect size from 649.22 to 634.56 (Supplementary Material, Table 3). The summary effect size of overall richness per 100 cm3 soil was most influenced by the removal of the control treatment from Kapagianni et al. (2010) from 28.97 to 28.70 (Supplementary Material, Table 4). These results indicated that no single study changed any of the summary effect sizes to any important degree. Funnel plots did not show any observable patterns between standard errors and point estimate values, indicating no publication bias in this meta-analysis. In addition, the Begg and Mazumdar rank correlation test gave absolute Kendall tau values for all four summary effect sizes of less than 0.22, suggesting no publication bias.
Overall nematode richness expressed per 100 g soil was highest in forest compared to NGL, DGL, urban, and agriculture (
Effect of ecosystem on genus-level nematode richness. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing overall richness of nematodes per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode richness. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing overall richness of nematodes per 100 cm3 of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode richness of each c-p guild. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing richness of nematodes at c-p guilds 1 to 5 per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode richness of each c-p guild. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing richness of nematodes at c-p guilds 1 to 5 per 100 cm3 of soil in different ecosystems. Letter
The richness of bacterivores, fungivores, and predators per 100 g of soil was higher in forest ecosystems than in the other ecosystems and the richness of omnivores per 100 g of soil was higher in forest ecosystems than in disturbed ecosystems. The richness of herbivores per 100 g of soil was higher in forest ecosystems than in agricultural ecosystems (
Effect of ecosystem on genus-level nematode richness of each trophic group. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing richness of nematodes of each trophic group per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode richness of each trophic group. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing richness of nematodes of each trophic group per 100 cm3 of soil in different ecosystems. Letter
The overall abundance of nematodes per 100 g of soil was similar in forest, NGL, DGL, and agricultural ecosystems. However, overall abundance of nematodes per 100 g of soil was higher in agricultural ecosystems than in urban ecosystems (
Effect of ecosystem on genus-level nematode abundance. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing overall abundance of nematodes per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode abundance. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing overall abundance of nematodes per 100 cm3 in different ecosystems. Letter
The abundance of c-p 1 guild per 100 g of soil was higher in agriculture ecosystems than in NGL ecosystems; abundance of c-p 2 guild per 100 g of soil was higher in DGL and agriculture ecosystems than in urban ecosystems and abundance of c-p 3 guild per 100 g of soil was higher in agriculture ecosystems than in NGL and urban ecosystems In contrast, the abundance of c-p 4 and c-p 5 guilds per 100 g of soil was higher in undisturbed ecosystems than disturbed ecosystems (
Effect of ecosystem on genus-level nematode abundance of each c-p guild. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing abundance of nematodes at c-p guilds 1 to 5 per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode abundance of each c-p guild. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing abundance of nematodes at c-p guilds 1 to 5 per 100 cm3 of soil in different ecosystems. Letter
The abundance of bacterivores per 100 g of soil was higher in agriculture than in NGL and urban ecosystems; abundance of fungivores per 100 g of soil was higher in agriculture than in urban ecosystems and abundance of herbivores per 100 g of soil was higher in DGL ecosystems than in urban ecosystems. On the other hand, abundance of predators and omnivores per 100 g of soil was higher in undisturbed ecosystems than in disturbed ecosystems (
Effect of ecosystem on genus-level nematode abundance of each trophic group. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing abundance of nematodes of each trophic group per 100 g of soil in different ecosystems. Letter
Effect of ecosystem on genus-level nematode abundance of each trophic group. Mean values are the weighted summary effect sizes and the bars represent standard error for comparing abundance of nematodes of each trophic group per 100 cm3 of soil in different ecosystems. Letter
Soil nematode assemblages can serve as ecological indicators since different nematode taxa vary in their sensitivity to disturbances in a terrestrial ecosystem (Bongers, 1990; Neher et al., 2005). Extensive research has been conducted on abundance and richness of nematode assemblages in different ecosystems but very few studies have been conducted to compare the impact of disturbances on nematode abundance and richness among two or more ecosystems (Neher et al., 2005; Briar et al., 2007; McSorley and Wang, 2009; Cardoso et al., 2015). Recently, meta-analysis was conducted using the literature published on soil nematodes to analyze soil energy pathways in different ecosystems (Zhao and Neher, 2014) and the effect of organic and inorganic fertilizers on soil nematodes in croplands (Liu et al., 2016). Meta-analysis was conducted to study the collective impact of anthropogenic disturbances on nematode assemblages by comparing five ecosystems with a gradient of human disturbance. Disturbances that are considered anthropogenic include physical disturbances such as burning, tillage, soil solarization, and harvesting; chemical disturbances such as addition of organic amendments and inorganic fertilizers in agriculture ecosystems; heavy metal pollution; building and road construction in urban settings; and seeding, tillage, harvesting, fertilizer application, and grazing rate in DGL were considered as anthropogenic disturbances. Forests and NGL with little to no direct human intervention were considered as undisturbed ecosystems.
The results from data expressed per 100 g of soil show that the overall richness of nematodes was highest in forest ecosystems compared to NGL, DGL, agriculture, and urban ecosystems. These results supported the hypothesis that the richness of nematodes is higher in undisturbed ecosystems than in human-disturbed ecosystems (Wasilewska, 1979; Bongers and Bongers, 1998; Briar et al., 2007; Darby et al., 2007). These results were congruent with the general statement that ecosystems with less or no disturbance support greater richness of soil biota (Hooper et al., 2005) consistent with the results of Hanel (1993), Ivezic et al. (2000), Neher et al. (2005), Brmez et al. (2007), Yeates (2007), Jiao et al. (2008), Cardoso et al. (2012, 2015). High richness in forest points to the stability of the ecosystem.
The richness of nematodes of all c-p guilds was higher in forest ecosystems due to little or no disturbance but the richness of c-p 1 was higher in agricultural ecosystems along with forest and NGL ecosystems. Nematodes in the c-p 1 guild are considered enrichment opportunists as most are bacterial feeders, which are most active in the presence of abundant resources (De Goede et al., 1993). The high richness of c-p 1 taxa in agricultural ecosystems may be due to continuous addition and incorporation of fertilizers and organic matter. After addition of nutrients or organic matter incorporation into the soil, c-p 1 guild nematodes respond immediately and flourish in number due to increased microbial activity, resulting from the newly available nutrients (Ettema and Bongers, 1993). Richness of nematodes in c-p 3, c-p 4 and c-p 5 guilds, which are sensitive to disturbance, was higher in forest ecosystem due to little or no disturbance. Nematodes of higher c-p guilds were found to be sensitive to disturbances (Park et al., 2010; Cardoso et al., 2015). High richness of higher c-p guilds indicates a mature and stable ecosystem (Bongers, 1990, 1999).
The richness of nematodes of all trophic groups except herbivores was highest in forest ecosystems. This result is consistent with the reports of Briar et al. (2007), Jiao et al. (2008), and Kimenju et al. (2009). Forests typically support a greater richness of organisms including nematodes due to the absence of human intervention such as tillage, monocultures, cultivated lawns, and application of fertilizers and amendments. Nematode trophic groups in the higher hierarchy of the soil food web such as omnivores and predators are particularly sensitive to disturbances (Korthals et al., 1996) and therefore are rich in undisturbed forest ecosystems. The presence of these nematodes maintains ecological balance by regulating nematode trophic groups in the lower hierarchy of the soil food web including plant feeding nematodes (Bilgrami and Brey, 2005).
Overall nematode abundance was similar in all ecosystems except urban ecosystems. Although high nematode abundance in an ecosystem represents high productivity of the ecosystem (Ritz and Trudgill, 1999), the high abundance in DGL and agriculture ecosystems was mostly attributed to high abundance of c-p 2, an indication of more stressful soil food web populated by recalcitrant bacterivores (Ferris et al., 2001). The higher abundance in forest and NGL ecosystems could be contributed by the higher abundance of predators and omnivores, most of which belong to c-p 4 and c-p 5 guilds.
The nematodes of c-p 1, c-p 2, and c-p 3 guilds were similar in all ecosystems, whereas the abundance of nematodes of c-p 4 and c-p 5 guilds was highest in forest and NGL ecosystems. The similar abundance of lower c-p guilds in disturbed ecosystems along with undisturbed ecosystems may be attributed to the incorporation of plant material and fertilizers, which favor microbial activity; thus, microbivorous colonizers with a high reproduction rate dominate these disturbed ecosystems (Bongers, 1990; Freckman and Ettema, 1993; Brmež et al., 2006; Brmez et al., 2007). Moreover, nematodes of lower c-p guilds are tolerant to disturbance (Bongers, 1990). On the other hand, the abundance of nematodes of higher c-p guilds, which are sensitive to disturbances, was highest in undisturbed ecosystems, which might be due to the absence of anthropogenic intervention such as tillage and fertilizer applications (Wasilewska, 1995; Grewal et al., 2011). High abundance of higher c-p guilds indicates mature soil food webs in an ecosystem (Neher, 1999; Yeates and Bongers, 1999).
The abundance of bacterivores, fungivores and herbivores was highest in DGL and agriculture ecosystems, whereas the abundance of predators and omnivores was highest in forest and NGL ecosystems. These results are consistent with the findings of Ivezic et al. (2000), Hanel (1993) and Háněl (2010). The abundance of nematode trophic groups in the lower hierarchy of soil food web is highest in disturbed ecosystems because bacterivores and fungivores with c-p 2 are tolerant and responding to more stressful soil environment (Bongers, 1990). High abundance of herbivores in disturbed ecosystems may be due to lack of omnivores and predators that potentially feed on herbivores. On the other hand, the high abundance of predators and omnivores in forest and NGL ecosystems may be due to lack of human intervention (Ferris and Ferris, 1974; Wasilewska, 1979, 1995; Hanel, 1993; Cardoso et al., 2012). Perturbations in an ecosystem may increase the abundance of trophic groups in the lower hierarchy of soil food web (bacterivores, fungivores, and herbivores) but decrease the abundance of nematode trophic groups in the higher hierarchy of the soil food web (predators and omnivores), which play a crucial role in regulating the lower groups including herbivores. Therefore, losing these regulators may be detrimental to nutrient cycling dynamics and agricultural management.
Overall richness, overall abundance, and richness and abundance of each c-p guild and each trophic group per 100 cm3 of soil in all four ecosystems were analyzed as no urban ecosystem studies using 100 cm3 were available. Summary effect sizes of overall richness, richness and abundance of c-p 4, predator and omnivore richness were not significantly different (Table 1). The overall abundance, abundance of nematodes of all c-p guilds, and abundance of nematodes of all trophic groups expressed per 100 cm3 of soil followed a somewhat similar pattern as that of 100 g of soil. However, overall richness, richness of all c-p guilds, and richness of all trophic groups expressed per 100 cm3 differed from those for 100 g of soil. This ambiguity may be due to the lower numbers of studies that used abundance per 100 cm3 than that expressed in per g soil, low statistical power, or the variation in the quantity of soil depending on its compactness, bulk density, and soil moisture.
Comprehensive meta-analyses of distinct ecosystems with different schemes of human intervention from 111 publications, using random effects model and non-parametric variance, confirmed that nematode richness was higher in less disturbed ecosystems (forest and NGL) compared to more disturbed ecosystems (agriculture, DGL, and urban ecosystems), nematode abundance of trophic groups in the lower hierarchy of the soil food web was higher in more disturbed ecosystems and nematode abundance of trophic groups in the higher hierarchy of the soil food web was higher in less disturbed ecosystems, consistent with general findings from previous works in the field of nematode ecology.