Next Article in Journal
Biodiversity and Phylogeny of North Atlantic Euphrosinidae (Annelida)
Next Article in Special Issue
Soil Invertebrate Communities as Indicator of Ecological Conservation Status of Some Fertilised Grasslands from Romania
Previous Article in Journal
Ice Cod Arctogadus glacialis (Peters, 1874) in Northeast Greenland—A First Sketch of Spatial Occurrence and Abundance
Previous Article in Special Issue
Advancing Mitochondrial Metagenomics: A New Assembly Strategy and Validating the Power of Seed-Based Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Assemblages of Two Groups of Collembola (Strong Furca and Weak Furca) under Different Agricultural Management Systems, Northeastern China

1
State Environmental Protection Key Laboratory of Wetland Ecology and Vegetation Restoration, School of Environment, Northeast Normal University, Changchun 130117, China
2
Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China
3
Key Laboratory of Vegetation Ecology, Ministry of Education, Northeast Normal University, Changchun 130024, China
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(11), 994; https://doi.org/10.3390/d14110994
Submission received: 16 October 2022 / Revised: 12 November 2022 / Accepted: 15 November 2022 / Published: 17 November 2022
(This article belongs to the Special Issue Soil Fauna Diversity under Global Change)

Abstract

:
Agriculture itself has been considered one of the leading reasons for biodiversity loss because of the huge quantity of land dedicated to just this activity, and agricultural intensification has impacted soil organisms at several taxonomic levels. Soil biota is a vital component of the agricultural system, providing essential ecosystem services while also having synergistic impacts on crop yield. Preservation of their diversity becomes a major element of an agricultural sustainability strategy. Many studies focused on agricultural activities’ effects on soil organisms, but few of them have focused on their effects on the co-occurrence patterns of their communities. Collembola communities are frequently employed as a substitute for soil organisms; thus, as a surface-dwelling arthropods representative, we investigated assemblages of soil Collembola in reaction to the arrangement of 6 treatments varying in crop rotation (MC: monoculture (corn) vs. CS: corn–soybean rotation) and tillage types (MP: mould ploughing; RT: ridge tillage; and NT: no-tillage). We hypothesized that Collembola communities with strong furca would respond well to the agricultural practices than those with weak furca, and there would be strong co-occurrence between species of Collembola communities belonging to the treatments with less intensity of soil disturbance and more variation in crop rotation. Our study found no clear evidence of a beneficial effect of crop rotation on Collembola communities. Although Collembola with strong furca shows higher abundances in plots with mold plowing, weak furca abundances were not showing any difference (differences in abilities to move fast from harsh habitats could be the reason for different responses of these two groups). Network analysis revealed that Collembola assemblages seem to occur more responsive to tillage intensity than crop rotation. Network graphs of treatments with ridge tillage are significantly more clustered than all others. For the first time, we can show that assemblages of springtails in agriculture were distinguished by a pattern of co-occurrence alongside agricultural practices (crop rotation, soil tillage), showing variations in the disturbance of soil and soil nutrients. Our results, contrary to our expectations, demonstrated that the effects of agricultural activities on Collembola abundance and diversity could be weak after long-term application of the same treatment, but still, they will clearly affect the bonds between Collembola species by affecting their co-occurrence pattern in Collembola communities.

1. Introduction

The primary objective of ecology is determining how communities change alongside environmental circumstances, and it is a necessity for estimating ecological reactions to global changes [1,2]. Agriculture itself has been considered one of the leading causes of biodiversity loss because of the huge quantity of land dedicated to just this activity [3,4]. Agricultural intensification has impacted soil organisms at several taxonomic levels. [5,6]. Agricultural practices like grazing, tillage, crop rotation, drainage, and exhaustive use of fertilizers and pesticides have major impacts on the soil flora and fauna [3,7]. However, the effect on soil organisms, especially Collembola, is extremely variable and highly dependent on the character and frequency of agricultural techniques applied, such as crop rotation, soil tillage, pesticide use, etc. [8,9]. Decreased or no-tillage techniques, on the other hand, can be beneficial for preserving populations of native species [3]. Global demands for food manufacturing are rising, as is the requirement for sustainable agricultural intensification [10]. Knowing ways to maximize yields while minimizing the loss of structure and function of soil is vital for agriculture’s long-term viability [11]. It is essential and helpful to understand the impacts of agricultural management intensity on communities of soil fauna in order to sustain, re-establish, and enhance the roles of biodiversity in the soil function of agricultural areas [5]. Goals of conservation agriculture are to build a sustainable agricultural system that maximizes crop productivity while significantly decreasing adverse environmental effects [12,13]. To do this, agronomical techniques should incorporate ecological processes to depend on more efficiency in the self-regulation of agroecosystems [14,15]. Usually, sustainable systems are centered on minimizing the disturbance of soil through the reduced intensity of tillage, reduced nitrogen fertilization, and increased crop residue retention to replenish organic matter content in soil [16,17].
The soil biota performs a variety of services that are critical in defining the physical, chemical, and biological characteristics of soils, particularly agricultural soils [18,19,20]. Soil communities of mesofauna are made up of populations of various species (for example, Collembola) that interact with one another and are possibly impacted by environmental factors and the management of soil [21]. Collembola performs a critical role in ecological processes and soil performances. For instance, they have a substantial impact on nutrient cycling and microbial community control by consuming soil microbes and dead organic debris, therefore influencing plant production [22,23,24,25,26,27]. Many writers have suggested the study of springtails as an indicator of ecological conditions in this situation [28,29,30]. Significant changes have been observed in the quantity and diversity of soil fauna, including both pests and useful organisms, by variations in tillage, rotation, and residue methods. Tillage alters the soil’s chemical and physical qualities, especially affecting the matrices that enable the growth of microbial populations [31,32,33]. The composition and structure of a community are widely acknowledged as indicators of ecological status [34,35,36], as well as the utilization of biological data to measure ecological health is a hot topic of research nowadays [37]. Despite the fact that diversity can be utilized to distinguish the structure of ecosystems, another significant element of the system is the variation within the abundance of its constituents [38].
Generally, beneficial impacts of modern conservation agriculture strategies on communities of springtails have been noticed [39,40,41,42,43], but these effects might vary depending on the soil depth investigated [42]. Wardle [44] concluded as a result of 106 pieces of research that while switching from conventional to no-tillage treatment, Collembola is often mildly stimulated, most likely as a result of relieving stress. Additionally, it has been demonstrated that returning crop remains (organic material) has a significant effect on microclimatic conditions as well as increases trophic resources for decomposers [45,46], as well as springtails, also see [47]. Nevertheless, various research indicated contradicting results, with conservation approaches having no or a negative influence on springtails when contrasted to conventional treatment [5,17,48,49,50,51,52]. Thus, we remain unclear about the agricultural strategies that are beneficial to Collembola populations. The impacts of crop rotation and monoculture on soil mesofauna are unclear in a sum of studies. Akkerhuis, et al. [53] discovered a considerably higher number of Mites in crops grown on a rotation of six-year than in crops grown on a rotation of three-year, but no significant differences in Collembola between treatments were found. According to Rebek, et al. [54], lower input cropping systems (grasslands) are more favorable to epitaphic Collembola than higher management cropping systems (monoculture). Gruss and Twardowski [55] discovered that the collembolan assemblage’s diversity was not significantly differed between crop rotation and monoculture, and the structure of species was typically similar for each treatment. Additionally, numerous studies compared only one agricultural technique, such as residue management, tillage type, or nitrogen fertilization rate. The impacts of merging these different techniques have attracted little attention. In the aboveground systems, both physical disruption and resource availability are assumed to be major attributes of biological variety [48,56].
In this study, we did an experiment to investigate the responses of Collembola communities of two groups of Collembola differed based on the strength of furca (strong furca and weak furca) under different treatments varying according to tillage (mould plow, ridge tillage, no-tillage) and crop rotation (continuous corn, corn–soybean) to compare their different reactions. The main objective of this study is to estimate the differences between the responses of two groups of Collembola (strong furca and weak furca) to agricultural practices, so we can choose the group that responds well for further research and explore the co-occurrences of these Collembola assemblages by using Network analysis for the first time. In order to mitigate the objectives of this study, we hypothesized that: (1) Collembola communities with strong furca will give good response to the agricultural practices than Collembola communities with weak furca because of their fast movement. (2) There will be strong co-occurrence between species of Collembola communities belonging to the treatments with less intensity of soil disturbance and more variation in crop rotation. This study will be able to enrich our understanding of the different reactions of these two groups of Collembola and select one of them in future analysis for a better understanding of the responses of Collembola communities to agricultural techniques.

2. Materials and Methods

2.1. Study Site

The Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Experimental Station (44°12′ N, 125°33′ E) was used for this study since it is located in a continental temperate monsoon zone. The mean yearly temperature is 4.4 degrees Celsius, and the mean yearly precipitation is 520.3 mm, with most of the precipitation (70%) between June and August. Prior to the experiment, the experimental location had already been homogenized using maize (corn) as monoculture and conventional tillage treatment for more than 30 years. The experiment was started in the fall of 2012 on a clay loam soil classed as Black soil (average soil texture was 360 g kg−1 clay, 245 g kg−1 silt, and 395 g kg−1 sand) [57].

2.2. Experimental Design

A randomized block design was used to conduct the experiment. Diverse tillage methods such as mould plow (MP), ridge tillage (RT), and no-tillage (NT), as well as two crop rotation systems, monoculture (corn) and corn–soybean rotation, were used in different combinations to differentiate the 6 treatments being used in this experiment. The MP system included residue return after harvest, drop mouldboard plow (20 cm), and secondary seedbed preparations in the spring, which included disking (7.5–10 cm), ridge-building, and harrowing. In RT, ridges were generated using a scrubber and a moderated lister and kept maintained with a cultivator, and the residue was restored in June of every year. For the NT, without disturbing the soil unless when planting with a no-till planter. Soybean remains were restored straight to the surface of the soil, while corn residue was chopped into roughly 30 cm parts and placed on the surface soil with 30–35 cm of standing stubble after harvesting. For corn, starter fertilizer was sprayed using a planter at the rate of 89-kg N ha–1, 51-kg P ha–1, and 51-kg K ha–1, while for soybean, the rate was 40-kg N ha–1, 49-kg Pha–1, and 53-kg K ha–1. The rates of application of N, P, and K stayed similar across all treatments of tillage. Four replicates of all six treatments were prepared in this design (24 plots in total). Each four-replication treatment area measured 5.2 m × 30 m, with 5 m buffers rows across treatments. The six treatments observed in this experiment are:
  • MP-MC: Mould plowing among monoculture (corn)
  • MP-CS: Mould plowing among corn–soybean rotation
  • RT-MC: Ridge tillage among monoculture (corn)
  • RT-CS: Ridge tillage among corn–soybean rotation
  • NT-MC: No-tillage among monoculture (corn)
  • NT-CS: No-tillage among corn–soybean rotation

2.3. Sampling Design

The current study sampled the corn phase of all six treatments. At the study site, we collected surface living Collembola using pitfall traps [58,59]. Pitfall traps expose the captured species’ movement activities and thus can be used to estimate their activity density [60]. A total of three samples were chosen from each of the 24 plots (northeast, central and southwest) (in total, 72 sampling points). In September 2020, we set a pitfall trap (10 cm depth, 5.6 cm diameter) at each sampling point [59]. The kind of trap used for sampling was with alcohol (100%) in a cup for each of 3 replicates of samples. They were spaced approximately 2 m apart along a transect running through the interrow in autumn, when the majority of Collembola species are active, allowing for the capture of most surface-active species. After a seven-day exposure period (15–22 September 2020), the traps were removed. All captured invertebrate specimens were preserved in pure alcohol (100%) for subsequent analysis [61]. The springtails were placed on slides one by one with 100 percent alcohol. They were recognized and then counted to the species level using an optical microscope [62,63,64,65]. Individuals (springtails) captured in alcohol were numbered and recognized utilizing classification keys [65], and we compared the sequences in Gen Bank and also identified traits of Collembola because our identification mainly relied on the traits of Collembola and our gene sequences were mainly auxiliary functions, so we did not get all the data. We have provided the names of all the species (Supplementary Data S2). After that, we collected three specimens of each specie and performed DNA extraction using a DNA extraction Kit. The extracted DNA can be immediately carried out for the next experiment or stored at −20 ℃. We use two types of primers (LCO1490 and HC02198) to prepare extracted DNA for further PCR analysis. After PCR, we send it to the company to get COI gene sequence data to differentiate between species (For the identification of all the species, please see Supplementary Data S1). Analyses were carried out at the system level of soil management. The total of individuals (abundance) and a total of species (specie richness) was then computed [66]. Later we divided the data into two parts (1) Species with strong furca and (2) Species with weak furca. Furca is an important functional trait, and strong/weak furca represents the development degree of different furca. Significant differences in their living habits and range of activities; therefore, they play different functions in the farmland ecosystem. Collembola with strong furca must have a retinaculum to manage the furca when they do not jump or move, which are tucked under the abdomen. Collembola with weak furca seldom use their retinaculum; this is our criteria. We did not upload the gene sequence to GenBank, and our identification work mainly relied on morphological identification. Thus, the data were collected as; Strong furca abundance (SA), Strong furca species richness (SSR), Weak furca abundance (WA), and Weak furca species richness (WSR).

2.4. Statistical Analysis

The numbers of Collembola have been utilized to determine density prior to statistical data analysis (ind m-2). The data set of abundance and species richness have been subjected to a normality and homogeneity test. Using the R software, the main effects of analysis of variance (ANOVA) were operated to identify the considerable variances between the Collembola abundances and community compositions of different experimental treatment. In order to figure out what the impacts of the agricultural activities on soil Collembola assemblages are, operating R software, version 4.0.4, network analysis was used to visualize a correlation matrix of Spearman rank based on springtail biomass data [67,68]. The co-occurrence pattern of collembolan taxa was investigated for each treatment using network analysis. Species showing relative abundances higher than 0.05 percent were chosen. If Spearman’s correlation coefficient (r) is >0.6 and the p-value is <0.01, Spearman’s correlation among the two species was declared statistically robust [69]. Every robust correlation detected through pairwise comparisons of abundances of genera creates a correlation network, with every single node representing a single genus and every single edge representing a significant and strong association among the nodes. To characterize the produced networks topology, a set of parameters (numbers of nodes and edges, clustering coefficient, average path length, graph density, average degree, modularity, and network diameter) were measured in the R environment utilizing the igraph package [70], and an interactive platform called Gephi was used to visualize the networks [71,72,73]. However, 10,000 Erdős-Rényi randomized networks were produced to correlate the original network’s topology, with every edge holding the equal possibility of giving to any node [74]. Unless otherwise stated, all the statistical analyses were contained using the R environment (http:// www.r-project.org, accessed on 25 September 2021).

3. Results

3.1. Community Structure Indices

In total, 9098 surface-active Collembola were recorded from 9 taxonomic families and 35 species. Out of 35 species, specie 1, 2, 3, and 8 were the dominant species because they appeared in all the treatments with high abundances almost (>10%), and species 17, 28, 30, and 35 were the rare species that appeared only once in one treatment out of 6 treatments. Collembola having strong furca were counted 3562 individuals, and Collembola having weak furca were counted 5666 individuals out of a total number of Collembola. Out of 35 species, 16 species have strong furca, and 19 have weak furca. From our six treatments, the lowest total number of Collembola was found in NT-MC (793 individuals), and the highest number of total Collembola was found in NT-CS (2854 individuals). The species richness of all the treatments varied between 9 to 28, although the overall species richness did not change significantly among treatments at the 5% error level (p < 0.05). The Hypogastruridae family was the most abundant, accounting for 54.55% of total individuals. The Entomobryidae family was the second most abundant, accounting for 36.97% of total individuals. In this experiment, four dominant species that show the highest abundances (>10%) in almost all the six treatments belong to Entomobryidae and Hypogastruridae families. (Supplementary Data S2) Therefore, these two families are the most dominant families found in our research experiment.

3.2. Different Responses of Species Richness and Abundances of Two Collembola Groups (Strong Furca and Weak Furca) to Different Tillage Types and Rotation Types

The ANOVA-based analysis of both data sets of two Collembola groups displayed those abundances of Collembola belonging to strong furca species differed significantly among tillage types (MP > RT > NT) but did not show any variations based on rotation type. The capitalized letters (A, B, AB) indicate the significant differences between tillage types (mold plowing, ridge tillage, and no-tillage). Different letters mean they have significantly different abundances, A is significantly higher than B, and AB means it does not have significant differences with both A and B, and the lowercase letter (a) indicates no significant differences between the two rotation types (Table 1; Figure 1A). However, the specie richness of strong furca species remains unchanged by both factors (Table 1; Figure 1B). Two-way ANOVA revealed that no agricultural management element (soil tillage type, plant rotation type) had a significant influence on specie richness and abundances of Collembola having weak furca (Table 1; Figure 1C, D). There were no significant changes discovered between the diversity (Shannon index), relative abundance (Simpson index), and evenness (pielou) of the six types of treatments observed in this research (Table 2; Figure 2A–C).

3.3. Co-Occurrence Pattern of Collembola Communities under Six Different Treatments

Collembola networks for each community consist of different numbers of nodes (genera) and edges (Table 3). The nodes in the networks were assigned to nine Collembola families (Figure 3). Significant topological characteristics were evaluated to explain the intricate structure of inter-relationships between nodes (taxa) of each treatment under the experiment and to compare these networks with each other (Table 3). This network analysis shows clear differences between the network graphs of different types of soil tillage, but there are no clear differences between networks belonging to different crop rotation types (Figure 3). Here is the pattern for the types of soil tillage showing the strong networks among all the treatments (RT > NT > MP). For the RT-MC network with 25 nodes and 72 edges (with 2.88 edges/node mean), the average path length (APL) was 2.091 edges with a diameter of 4 edges. Graph density was 0.240 as well as the average degree was 2.880. In terms of clustering and modularity, our results came in at 0.629 (CC) and 0.341 (MD), respectively; values > 0.4 indicate that this network comprises a modular structure (Table 3). The structural characteristics of this network indicated that this network graph of Collembola was constructed of highly associated genera and presented a “small world” topology (it is, the RT-MC network graph is more significantly clustered than the other five graphs) (Table 3, Figure 3a).

4. Discussion

The focus of this research is to differentiate the responses of abundances of two groups of Collembola species that differed on the basis of the strength of Collembola furca and to examine the effects of agricultural activities on springtail co-occurrences. Due to the fact that all of the studied plots had previously been managed with various treatments for at least 12 years using various types of tillage (MP, RT, and NT) and crop rotations (continuous corn, corn–soybean), we were capable of assessing the long-term effects on soil Collembola assemblages due to soil tillage and rotation type in particular. Our results were at odds with previous studies [75,76,77,78] That suggested that the tillage influences on the soil communities would be alleviated by the establishment of plant roots throughout the growth period. These variations may vary depending on the conservation tillage duration, which was much longer in our study than in the studies cited (long-term, > 10 years, and 10 years short-term). It may not be enough to apply tillage techniques for less than 10 years to develop long benefits in soils [79,80]. The results showed that the abundance of Collembola with strong furca increased knowingly under conventional tillage (MP) compared to reduced tillage (RT, NT) in both rotation types. Meanwhile, Collembola with strong furca decreased under zero tillage with both rotations (continuous corn and corn–soybean rotation) compared with conventional tillage. On the other hand, tillage and rotation management did not significantly influence the abundance of Collembola with weak furca, corroborating Hypothesis 1 that Collembola with strong furca will respond well to these agricultural practices. Certain species are known to achieve their highest abundance in the spring, whereas others are highest in the autumn [81]. Loring, et al. [82] recommended that certain collembolan morphological types could be targeted by distinct tillage practices [83]. In a study population of Pratylenchus thornei with other plant-parasitic and non-plant parasitic nematodes have been assessed in another research. It shows that continuous maize had a greater P. thornei population than corn–wheat rotation, while conventional tillage had a higher P. thornei population than No-tillage [84]. Few significant impacts on total Collembola abundance were found by Sabatini, et al. [85], which can be attributed to differences in long-term conventional plowing and reduced tillage. However, some species, such as Protaphorura armata and Isotoma notabilis, were favored by reduced tillage, whereas the others were favored by conventional plowing. MooRE, et al. [86] found that no-tillage had a detrimental effect on Collembola (springtails) in comparison to conventional plowing. On the other side, the differences in agricultural practices between plots did not result in a growth in collembolan species diversity of both Collembola groups, corroborating with certain prior research that failed to demonstrate an influence of agricultural practices on soil animal communities. This might be because (1) the tillage and rotation types highly appropriate for Collembola were underreported, (2) apparently little Collembola diversity in all experimental plots excludes a growth in overall species richness under different treatments, and (3) Collembola communities learn to adjust and normally work under harsh condition after long-term continuous application of harmful agricultural practices as well as (4) soil characteristics, particularly texture (particle size distribution), influenced Collembola distribution despite our best efforts (randomized seeding) to account for this natural result in our experimental design. In addition, a possible reason can be the effects of these agricultural activities on soil; the Collembola community structure may be weakened after a 12-year time period. For example, in an experiment, although treatment shows differences in abundances, the plots under conventional tillage were very similar to the plots with ECO tillage treatments in terms of species diversity and trends [83]. Another study found no significant difference in the average species richness between management systems [87]. Moreover, several factors, e.g., environmental circumstances, soil characteristics (soil pH, soil type, organic matter content, soil structure, etc.), species community (soil fauna, plants, and microbes), and period of sampling can affect the density, composition and species activity [8,88,89,90,91]. This change could be explained by regional or local differences in Collembola abundance [92], but there were no differences found in the Shannonindex, Simpsonindex, and pielou of the total number of Collembola between six different treatments studied in our experiment. According to a study, there was very little difference between NT and CLI systems in terms of species richness, biodiversity (Margalef), evenness (Pielou), and diversity (Shannon-Wiener) [87]. In another study, it was discovered that the soil management technique tested had no effect on the total number of individuals or their richness [87].
Coulibaly et al. [6] result demonstrated an overall rise in springtails species richness and density through decreasing tillage intensity, which is consistent with several past pieces of research [42,44,83]. However, Capelle et al. [51] found contradictory results that showed no change in Collembola abundance between reduced and conventional tillage and observed a greater number of collembolan in conventional tillage in comparison to the systems with no-tillage, in a German data-based meta-analysis. Capelle et al. [51] showed that tillage systems with reduced-tillage intensity enhanced the total earthworm population, species diversity, and biomass considerably. At the same time, Nakamoto, et al. [93] discovered a similar effect on nematode counts. House and Parmelee, on the other hand, discovered that under no-tillage, all main microarthropod taxa, including Collembola, had higher abundances than under conventional tillage. However, according to another study, Collembola species richness was positively correlated with higher tillage intensities. Under grapevines, the inter-row Collembola activity density was positively linked with inter-row plowing and herbicide applications [61]. However, in our research, there were no noticeable impacts of soil tillage or rotation type on specie richness, with each treatment having a similar average number of species. Widespread, Common and eurytopic species dominated the assemblages. Out of the thirty-five species obtained, four (species 1, 2, 3, 8) were more abundant in almost all the treatments. Entomobryidae and Hypogastruridae families were the most abundant families accounting for 91.5% of the total number of Collembola collected during this experiment. The latter was apparently attributable to these families’ ability to reproduce quickly and invade new areas. In a study, it was also discovered that species similar to both systems (conventional and conservation) had greater densities and a higher difference in their densities than those unique species [87]. It is essential to notice that only one -time sampling was done during this research (in September), and it is well known that the communities of Collembola alter greatly during the year. The management techniques had a restricted influence on the average abundance and species richness of Collembola. Nevertheless, springtail community co-occurrence differed according to soil management systems, corroborating Hypothesis 2 that agricultural practices will alter species co-occurrence and community structure. There are clear differences between the network graphs of different communities of Collembola belonging to different tillage types that can be seen visually. However, we must acknowledge that crop rotation type had no influence on Collembola populations. Although it has been repeatedly demonstrated that plant species richness and diversity encourage the growth, densities, and diversification of soil fauna [94,95], possibly by increasing the quality or diversity of underground resources, numerous studies have found no significant, or even contradictory, effects of plant variety on soil organisms [95], emphasizing the lack of confidence in the obtained results. Differences in the influence of plant diversity on Collembola species richness seen in several research suggest that these influences are minor or disguised by interaction processes, all of which confirm our findings. Crop diversity’s effect may be influenced by additional factors, for example, farming systems (e.g., organic vs. conventional farming), the surrounding landscape [96,97], and soil characteristics [98]. Besides that, site conditions might well influence how soil organisms respond to changes in land use. Soil organisms living in the topsoil (0–10 cm) more accurately reflect present management practices than soil biota living deeper in the soil, which may more accurately reflect management history as well as properties of soil [99]. Only monoculture planting in experimental plots can also be the reason for not significant effects of rotation type on Collembola assemblages.
According to our research, the ANOVA and network analysis results show that tillage type is the primary factor responsible for the fate abundances of strong furca and assemblages of springtails regarding co-occurrence and species composition. Certainly, Collembola seemed to react more to the tillage intensities than to some other management intensities [6]. Other studies, however, suggest that springtails are not always harmed by the mechanical disturbance caused by repeated tillage, at least in comparison to alternate tillage on rather compacted soils [59,100]. However, to understand the underlying processes, more comprehensive studies with additional sampling for every season are required. It’s also worth noting that we just examined epigeic Collembola species, and other life forms of Collembola could respond differently or more accurately because techniques for calculating the ecomorphological index show that euedaphic Collembola species are more responsive to many mechanical disturbances than epigeic Collembola species [101]. Still, there is a gap in the knowledge about the response of springtails present in the agricultural ecosystem. Therefore, further studies are needed to have a deeper understanding of this kind of production system, apart from improving our ability to better understand the relationships among springtails, environmental factors, and certain soil properties, as well as to identify the cause and effects of springtails and different responses of different groups of Collembola species to understand which group can give a better understanding of this relationship of Collembola and their agricultural environment. Secondly, some of the mysterious differences could be ascribed to unmeasured crop as well as soil characteristics, for instance, pore spaces factors (tortuosity, size distribution) and roots of the crop, so further research must be done in this field to understand the impacts of tillage and crop rotation on soil collembolan assemblages including other important factors more deeply.

5. Conclusions

This study, along with the abundance and richness of Collembola, allowed us to investigate the co-occurrence pattern of collembolans under six distinct experimental treatments on one site. Our study established that the abundances of Collembola with strong furca varied considerably across the tillage types, but the abundances of Collembola with weak furca did not respond to tillage or rotation type. Although, tillage had a greater impact on the co-occurrence pattern of Collembola assemblages than crop rotation. This reveals that collembolan assemblages respond to mechanical disturbances. Network graphs of treatments with ridge tillage were more significantly clustered than other graphs. These results indicate that Collembola with strong furca responds better to these disturbances than Collembola with weak furca, but it will affect the community structure of total Collembola by affecting co-occurrence between species. In further studies, seasonal sampling is required to predict the dynamic pattern of Collembola assemblage. Secondly, investigating all Collembola life forms could react more clearly to these agricultural activities, and their co-occurrence patterns will give a better understanding of this relationship.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d14110994/s1, Supplementary Data S1: Collembola species information that was identified during the research; Supplementary Data S2: Table S1. The dominance of each species except species 20 and 31 (because of 0% in all treatments) in 6 types of treatments used in our experiment is shown in this table by (***) for more than 10% of the total number of Collembola collected from the treatment type mentioned in the column, (**) for 1-10% of the total number of Collembola collected from the treatment type mentioned in the column and (*) for 0-1% of the total number of Collembola collected from the treatment type mentioned in the column.

Author Contributions

Conceptualization, T.H. and D.W.; Data curation, T.H., D.W. and L.C.; Funding acquisition, D.W.; Methodology, T.H. and D.W.; Writing—original draft, T.H.; Writing—review and editing, T.H., D.W., S.L., Y.W. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA28020201), the National Natural Sciences Foundation of China (No. 42071059), the Program of Introducing Talents of Discipline to Universities (No. B16011), the National Science & Technology Fundamental Resources Investigation Program of China (2018FY100300).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We greatly appreciate our laboratory colleagues, such as Zhuoma Wan, Pinting Guan, Qiaoqiao Ji, Xue Pan, and Jingjing Yang, who provided us with field sampling support and whom we could always address in order to solve any problem during preparing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fukami, T.; Wardle, D.A. Long-term ecological dynamics: Reciprocal insights from natural and anthropogenic gradients. Proc. R. Soc. B Biol. Sci. 2005, 272, 2105–2115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Laliberté, E.; Kardol, P.; Didham, R.K.; Teste, F.P.; Turner, B.L.; Wardle, D.A. Soil fertility shapes belowground food webs across a regional climate gradient. Ecol. Lett. 2017, 20, 1273–1284. [Google Scholar] [CrossRef] [PubMed]
  3. McLaughlin, A.; Mineau, P. The impact of agricultural practices on biodiversity. Agric. Ecosyst. Environ. 1995, 55, 201–212. [Google Scholar] [CrossRef]
  4. González de Molina, M.; Soto Fernández, D.; Guzmán Casado, G.; Infante-Amate, J.; Aguilera Fernández, E.; Vila Traver, J.; García Ruiz, R. The Social Metabolism of Spanish Agriculture, 1900–2008: The Mediterranean Way Towards Industrialization; Springer Nature: Berlin, Germany, 2020. [Google Scholar]
  5. Postma-Blaauw, M.B.; de Goede, R.G.M.; Bloem, J.; Faber, J.H.; Brussaard, L. Soil biota community structure and abundance under agricultural intensification and extensification. Ecology 2010, 91, 460–473. [Google Scholar] [CrossRef] [PubMed]
  6. Coulibaly, S.; Coudrain, V.; Hedde, M.; Brunet, N.; Mary, B.; Recous, S.; Chauvat, M. Effect of different crop management practices on soil Collembola assemblages: A 4-year follow-up. Appl. Soil Ecol. 2017, 119, 354–366. [Google Scholar] [CrossRef]
  7. Ghosh, A.; Misra, S.; Bhattacharyya, R.; Sarkar, A.; Singh, A.K.; Tyagi, V.C.; Kumar, R.V.; Meena, V.S. Agriculture, dairy and fishery farming practices and greenhouse gas emission footprint: A strategic appraisal for mitigation. Environ. Sci. Pollut. Res. 2020, 27, 10160–10184. [Google Scholar] [CrossRef]
  8. Cluzeau, D.; Guernion, M.; Chaussod, R.; Martin-Laurent, F.; Villenave, C.; Cortet, J.; Ruiz-Camacho, N.; Pernin, C.; Mateille, T.; Philippot, L. Integration of biodiversity in soil quality monitoring: Baselines for microbial and soil fauna parameters for different land-use types. Eur. J. Soil Biol. 2012, 49, 63–72. [Google Scholar] [CrossRef]
  9. Mahfouz, M.; Mohamed, M. Towards optimization of entomopathogenic nematodes for more service in the biological control. J. Nematol. 2019, 51, 77. [Google Scholar]
  10. Garnett, T.; Appleby, M.C.; Balmford, A.; Bateman, I.J.; Benton, T.G.; Bloomer, P.; Burlingame, B.; Dawkins, M.; Dolan, L.; Fraser, D. Sustainable intensification in agriculture: Premises and policies. Science 2013, 341, 33–34. [Google Scholar] [CrossRef]
  11. Ball, B.; Bingham, I.; Rees, R.; Watson, C.; Litterick, A. The role of crop rotations in determining soil structure and crop growth conditions. Can. J. Soil Sci. 2005, 85, 557–577. [Google Scholar] [CrossRef] [Green Version]
  12. Bommarco, R.; Kleijn, D.; Potts, S.G. Ecological intensification: Harnessing ecosystem services for food security. Trends Ecol. Evol. 2013, 28, 230–238. [Google Scholar] [CrossRef] [PubMed]
  13. Hobbs, P.R. Conservation agriculture: What is it and why is it important for future sustainable food production? J. Agric. Sci. Camb. 2007, 145, 127–137. [Google Scholar] [CrossRef] [Green Version]
  14. Altieri, M.A. The ecological role of biodiversity in agroecosystems. In Invertebrate Biodiversity as Bioindicators of Sustainable Landscapes; Elsevier: Amsterdam, The Netherlands, 1999; pp. 19–31. [Google Scholar]
  15. Wezel, A.; Casagrande, M.; Celette, F.; Vian, J.-F.; Ferrer, A.; Peigné, J. Agroecological practices for sustainable agriculture. A review. Agron. Sustain. Dev. 2014, 34, 1–20. [Google Scholar] [CrossRef] [Green Version]
  16. Michels, K.; Sivakumar, M.; Allison, B. Wind erosion control using crop residue I. Effects on soil flux and soil properties. Field Crops Res. 1995, 40, 101–110. [Google Scholar] [CrossRef] [Green Version]
  17. Scheunemann, N.; Maraun, M.; Scheu, S.; Butenschoen, O. The role of shoot residues vs. crop species for soil arthropod diversity and abundance of arable systems. Soil Biol. Biochem. 2015, 81, 81–88. [Google Scholar] [CrossRef]
  18. Brevik, E.; Cerdà, A.; Mataix-Solera, J.; Pereg, L.; Quinton, J.; Six, J.; Van Oost, K. The interdisciplinary nature of SOIL. Soil 2015, 1, 117–129. [Google Scholar] [CrossRef] [Green Version]
  19. Brown, S.; Fu, G.; White, L.W. Complementary teaching approaches facilitating interdisciplinary soil science education. Can. J. Soil Sci. 2021, 102, 433–443. [Google Scholar] [CrossRef]
  20. Hou, D.; Bolan, N.S.; Tsang, D.C.; Kirkham, M.B.; O’Connor, D. Sustainable soil use and management: An interdisciplinary and systematic approach. Sci. Total Environ. 2020, 729, 138961. [Google Scholar] [CrossRef]
  21. Wahl, J.; Theron, P.; Maboeta, M. Soil mesofauna as bioindicators to assess environmental disturbance at a platinum mine in South Africa. Ecotoxicol. Environ. Saf. 2012, 86, 250–260. [Google Scholar] [CrossRef]
  22. Endlweber, K.; Scheu, S. Interactions between mycorrhizal fungi and Collembola: Effects on root structure of competing plant species. Biol. Fertil. Soils 2007, 43, 741–749. [Google Scholar] [CrossRef]
  23. Menta, C. Soil fauna diversity-function, soil degradation, biological indices, soil restoration. Biodivers. Conserv. Util. A Divers. World 2012. [Google Scholar] [CrossRef]
  24. Rzeszowski, K.; Zadrożny, P.; Nicia, P. The effect of soil nutrient gradients on Collembola communities inhabiting typical urban green spaces. Pedobiologia 2017, 64, 15–24. [Google Scholar] [CrossRef]
  25. Kitching, R.L.; Dahlsjö, C.A.; Eggleton, P. Invertebrates and the complexity of tropical ecosystems. Biotropica 2020, 52, 207–214. [Google Scholar] [CrossRef]
  26. Coulibaly, S.F.; Winck, B.R.; Akpa-Vinceslas, M.; Mignot, L.; Legras, M.; Forey, E.; Chauvat, M.J.A.S.E. Functional assemblages of Collembola determine soil microbial communities and associated functions. Front. Environ. Sci. 2019, 7, 52. [Google Scholar] [CrossRef] [Green Version]
  27. Winck, B.R.; Chauvat, M.; Coulibaly, S.F.; Santonja, M.; Saccol de Sá, E.L.; Forey, E. Functional collembolan assemblages induce different plant responses in Lolium perenne. Plant Soil 2020, 452, 347–358. [Google Scholar] [CrossRef]
  28. Frampton, G.K. The potential of Collembola as indicators of pesticide usage: Evidence and methods from the UK arable ecosystem. Pedobiologia 1997, 41, 179–184. [Google Scholar]
  29. Kopeszki, H. An active bioindication method for the diagnosis of soil properties using Collembola. Pedobiologia 1997, 41, 159–166. [Google Scholar]
  30. Van Straalen, N.M.; Verhoef, H.A. The development of a bioindicator system for soil acidity based on arthropod pH preferences. J. Appl. Ecol. 1997, 34, 217–232. [Google Scholar] [CrossRef]
  31. Andersen, A. Plant protection in spring cereal production with reduced tillage. II. Pests and beneficial insects. Crop Prot. 1999, 18, 651–657. [Google Scholar] [CrossRef]
  32. Bockus, W.; Shroyer, J. The impact of reduced tillage on soilborne plant pathogens. Annu. Rev. Phytopathol. 1998, 36, 485–500. [Google Scholar] [CrossRef]
  33. Vanhée, B.; Devigne, C. Differences in collembola species assemblages (Arthropoda) between spoil tips and surrounding environments are dependent on vegetation development. Sci. Rep. 2018, 8, 18067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Cairns, J.; Pratt, J.R. A history of biological monitoring using benthic macroinvertebrates. Freshw. Biomonitoring Benthic Macroinvertebrates 1993, 10, 27. [Google Scholar]
  35. Dickens, C.; Graham, P. Biomonitoring for effective management of wastewater discharges and the health of the river environment. Aquat. Ecosyst. Health Manag. 1998, 1, 199–217. [Google Scholar] [CrossRef]
  36. Carlisle, D.M.; Meador, M.R.; Moulton II, S.R.; Ruhl, P.M. Estimation and application of indicator values for common macroinvertebrate genera and families of the United States. Ecol. Indic. 2007, 7, 22–33. [Google Scholar] [CrossRef]
  37. Sandler, R.V.; Falco, L.B.; Di Ciocco, C.A.; Castro-Huerta, R.; Coviella, C.E. The Degree of Change of Collembolan Community Structure Related to Anthropic Soil Disturbance; PeerJ PrePrints: London, UK, 2014; pp. 2167–9843. [Google Scholar]
  38. da Fonseca, J.C.; Sarkar, S. Soil microarthropods in two different managed ecological systems (Tripura, India). Appl. Soil Ecol. 1998, 9, 105–107. [Google Scholar] [CrossRef]
  39. Cortet, J.; Ronce, D.; Poinsot-Balaguer, N.; Beaufreton, C.; Chabert, A.; Viaux, P.; de Fonseca, J.P.C. Impacts of different agricultural practices on the biodiversity of microarthropod communities in arable crop systems. Eur. J. Soil Biol. 2002, 38, 239–244. [Google Scholar] [CrossRef]
  40. Dittmer, S.; Schrader, S. Longterm effects of soil compaction and tillage on Collembola and straw decomposition in arable soil. Pedobiologia 2000, 44, 527–538. [Google Scholar] [CrossRef]
  41. Winter, J.; Voroney, R.; Ainsworth, D. Soil microarthropods in long-term no-tillage and conventional tillage corn production. Can. J. Soil Sci. 1990, 70, 641–653. [Google Scholar] [CrossRef] [Green Version]
  42. Petersen, H. Effects of non-inverting deep tillage vs. conventional ploughing on collembolan populations in an organic wheat field. Eur. J. Soil Biol. 2002, 38, 177–180. [Google Scholar] [CrossRef]
  43. Gupta, V.; Roper, M.; Thompson, J.; Pratley, J.; Kirkegaard, J. Harnessing the benefits of soil biology in conservation agriculture. In Australian Agriculture in 2020: From Conservation to Automation; Pratley, J., Kirkegaard, J., Eds.; Agronomy Australia and Charles Sturt University: Wagga Wagga, NSW, Australia; pp. 237–253.
  44. Wardle, D.A. Impacts of disturbance on detritus food webs in agro-ecosystems of contrasting tillage and weed management practices. In Advances in Ecological Research; Elsevier: Amsterdam, The Netherlands, 1995; Volume 26, pp. 105–185. [Google Scholar]
  45. Axelsen, J.A.; Kristensen, K.T. Collembola and mites in plots fertilised with different types of green manure. Pedobiologia 2000, 44, 556–566. [Google Scholar] [CrossRef]
  46. Kladivko, E.J. Tillage systems and soil ecology. Soil Tillage Res. 2001, 61, 61–76. [Google Scholar] [CrossRef]
  47. Ngosong, C.; Raupp, J.; Scheu, S.; Ruess, L. Low importance for a fungal based food web in arable soils under mineral and organic fertilization indicated by Collembola grazers. Soil Biol. Biochem. 2009, 41, 2308–2317. [Google Scholar] [CrossRef]
  48. Cole, L.; Buckland, S.; Bardgett, R.D. Influence of disturbance and nitrogen addition on plant and soil animal diversity in grassland. Soil Biol. Biochem. 2008, 40, 505–514. [Google Scholar] [CrossRef]
  49. Holland, J.M. The environmental consequences of adopting conservation tillage in Europe: Reviewing the evidence. Agric. Ecosyst. Environ. 2004, 103, 1–25. [Google Scholar] [CrossRef]
  50. Tsiafouli, M.A.; Thébault, E.; Sgardelis, S.P.; De Ruiter, P.C.; Van Der Putten, W.H.; Birkhofer, K.; Hemerik, L.; De Vries, F.T.; Bardgett, R.D.; Brady, M.V. Intensive agriculture reduces soil biodiversity across Europe. Glob. Chang. Biol. 2015, 21, 973–985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. van Capelle, C.; Schrader, S.; Brunotte, J. Tillage-induced changes in the functional diversity of soil biota–A review with a focus on German data. Eur. J. Soil Biol. 2012, 50, 165–181. [Google Scholar] [CrossRef]
  52. Moos, J.H.; Schrader, S.; Paulsen, H.M. Minor changes in collembolan communities under different organic crop rotations and tillage regimes. Landbauforsch. J. Sustain. Org. Agric. Syst. 2020, 70, 113–128. [Google Scholar]
  53. Akkerhuis, G.J.O.; De Ley, F.; Zwetsloot, H.; Ponge, J.-F.; Brussaard, L. Soil microarthropods (Acari and Collembola) in two crop rotations on a heavy marine clay soil. Rev. D’ecologie Et De Biol. Du Sol 1988, 25, 175–202. [Google Scholar]
  54. Rebek, E.; Hogg, D.; Young, D. Effect of four cropping systems on the abundance and diversity of epedaphic springtails (Hexapoda: Parainsecta: Collembola) in southern Wisconsin. Environ. Entomol. 2002, 31, 37–46. [Google Scholar] [CrossRef] [Green Version]
  55. Gruss, I.; Twardowski, J. The assemblages of soil-dwelling springtails (Collembola) in winter rye under long-term monoculture and crop rotation. Zemdirb. Agric. 2016, 103, 159–166. [Google Scholar] [CrossRef] [Green Version]
  56. Lawton, J.H. Community Ecology in a Changing World; Oldendorf/Luhe (Germany) Ecology Inst.: Oldendorf/Luhe, Germany, 2000. [Google Scholar]
  57. Baillie, I. Soil Survey Staff 1999, Soil Taxonomy: A Basic System of Soil Classification for Making and Interpreting Soil Surveys, Agricultural Handbook 436, Natural Resources Conservation Service, USDA, Washington DC, USA; Wiley Online Library: Hoboken, NJ, USA, 2001; p. 869. [Google Scholar]
  58. Thomson, L.; Neville, P.; Hoffmann, A. Effective trapping methods for assessing invertebrates in vineyards. Aust. J. Exp. Agric. 2004, 44, 947–953. [Google Scholar] [CrossRef]
  59. Buchholz, J.; Querner, P.; Paredes, D.; Bauer, T.; Strauss, P.; Guernion, M.; Scimia, J.; Cluzeau, D.; Burel, F.; Kratschmer, S. Soil biota in vineyards are more influenced by plants and soil quality than by tillage intensity or the surrounding landscape. Sci. Rep. 2017, 7, 17445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Melbourne, B.A. Bias in the effect of habitat structure on pitfall traps: An experimental evaluation. Aust. J. Ecol. 1999, 24, 228–239. [Google Scholar] [CrossRef] [Green Version]
  61. Fiera, C.; Ulrich, W.; Popescu, D.; Buchholz, J.; Querner, P.; Bunea, C.-I.; Strauss, P.; Bauer, T.; Kratschmer, S.; Winter, S. Tillage intensity and herbicide application influence surface-active springtail (Collembola) communities in Romanian vineyards. Agric. Ecosyst. Environ. 2020, 300, 107006. [Google Scholar] [CrossRef]
  62. Betsch, J.-M. Éléments pour une monographie des collemboles symphyplèones (hexapodes, aptérygotes). Éditions Du Muséum 1980, 116, 1–227. [Google Scholar]
  63. Bretfeld, G. Synopses on Palaearctic Collembola: Symphypleona. Abh. Und Ber. Des Nat. Gorlitz. 1999, 71, 1–318. [Google Scholar]
  64. Potapov, M. Synopses on Palaearctic Collembola: Isotomidae. Abh. Und Ber. Des Nat. Gorlitz. 2001, 73, 1–603. [Google Scholar]
  65. Bellinger, P.F.; Christiansen, K.A.; Janssens, F. Checklist of the Collembola of the World. Available online: http://www.collembola.org (accessed on 11 May 2021).
  66. Vandewalle, M.; De Bello, F.; Berg, M.P.; Bolger, T.; Doledec, S.; Dubs, F.; Feld, C.K.; Harrington, R.; Harrison, P.A.; Lavorel, S. Functional traits as indicators of biodiversity response to land use changes across ecosystems and organisms. Biodivers. Conserv. 2010, 19, 2921–2947. [Google Scholar] [CrossRef] [Green Version]
  67. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  68. Morriën, E.; Hannula, S.E.; Snoek, L.B.; Helmsing, N.R.; Zweers, H.; De Hollander, M.; Soto, R.L.; Bouffaud, M.-L.; Buée, M.; Dimmers, W. Soil networks become more connected and take up more carbon as nature restoration progresses. Nat. Commun. 2017, 8, 14349. [Google Scholar] [CrossRef] [Green Version]
  69. Barberán, A.; Bates, S.T.; Casamayor, E.O.; Fierer, N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012, 6, 343–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Csardi, G.; Nepusz, T. The igraph software package for complex network research. InterJournal Complex Syst. 2006, 1695, 1–9. [Google Scholar]
  71. Newman, M.E. The structure and function of complex networks. SIAM Rev. 2003, 45, 167–256. [Google Scholar] [CrossRef] [Green Version]
  72. Newman, M.E. Modularity and community structure in networks. Proc. Natl. Acad. Sci. USA 2006, 103, 8577–8582. [Google Scholar] [CrossRef] [Green Version]
  73. Bastian, M.; Heymann, S.; Jacomy, M. Gephi: An open Source Software for Exploring and Manipulating Networks. In Proceedings of the Third International AAAI Conference on Weblogs and Social Media, San Jose, CA, USA, 17–20 May 2009. [Google Scholar]
  74. Erdos, P.; Rényi, A. On the evolution of random graphs. Publ. Math. Inst. Hung. Acad. Sci 1960, 5, 17–60. [Google Scholar]
  75. Degrune, F.; Theodorakopoulos, N.; Colinet, G.; Hiel, M.-P.; Bodson, B.; Taminiau, B.; Daube, G.; Vandenbol, M.; Hartmann, M. Temporal dynamics of soil microbial communities below the seedbed under two contrasting tillage regimes. Front. Microbiol. 2017, 8, 1127. [Google Scholar] [CrossRef]
  76. Pausch, J.; Hofmann, S.; Scharroba, A.; Kuzyakov, Y.; Ruess, L. Fluxes of root-derived carbon into the nematode micro-food web of an arable soil. Food Webs 2016, 9, 32–38. [Google Scholar] [CrossRef]
  77. Spedding, T.; Hamel, C.; Mehuys, G.; Madramootoo, C. Soil microbial dynamics in maize-growing soil under different tillage and residue management systems. Soil Biol. Biochem. 2004, 36, 499–512. [Google Scholar] [CrossRef]
  78. Zhang, B.; He, H.; Ding, X.; Zhang, X.; Zhang, X.; Yang, X.; Filley, T.R. Soil microbial community dynamics over a maize (Zea mays L.) growing season under conventional-and no-tillage practices in a rainfed agroecosystem. Soil Tillage Res. 2012, 124, 153–160. [Google Scholar] [CrossRef]
  79. Six, J.; Ogle, S.M.; Jay Breidt, F.; Conant, R.T.; Mosier, A.R.; Paustian, K. The potential to mitigate global warming with no-tillage management is only realized when practised in the long term. Glob. Chang. Biol. 2004, 10, 155–160. [Google Scholar] [CrossRef] [Green Version]
  80. Zhang, S.; Chen, X.; Jia, S.; Liang, A.; Zhang, X.; Yang, X.; Wei, S.; Sun, B.; Huang, D.; Zhou, G. The potential mechanism of long-term conservation tillage effects on maize yield in the black soil of Northeast China. Soil Tillage Res. 2015, 154, 84–90. [Google Scholar] [CrossRef]
  81. Cole, L.J.; McCracken, D.I.; Foster, G.N.; Aitken, M.N. Using Collembola to assess the risks of applying metal-rich sewage sludge to agricultural land in western Scotland. Agric. Ecosyst. Environ. 2001, 83, 177–189. [Google Scholar] [CrossRef]
  82. Loring, S.; Snider, R.J.; Robertson, L.S. The effects of three tillage practices on Collembola and Acarina populations. Pedobiologia 1981, 22, 172–184. [Google Scholar]
  83. Brennan, A.; Fortune, T.; Bolger, T. Collembola abundances and assemblage structures in conventionally tilled and conservation tillage arable systems. Pedobiologia 2006, 50, 135–145. [Google Scholar] [CrossRef]
  84. Govaerts, B.; Mezzalama, M.; Sayre, K.D.; Crossa, J.; Nicol, J.M.; Deckers, J. Long-term consequences of tillage, residue management, and crop rotation on maize/wheat root rot and nematode populations in subtropical highlands. Appl. Soil Ecol. 2006, 32, 305–315. [Google Scholar] [CrossRef]
  85. Sabatini, M.A.; Rebecchi, L.; Cappi, C.; Bertolani, R.; Fratello, B. Long-term effects of three different continuous tillage practices on Collembola populations. Pedobiologia 1997, 41, 185–193. [Google Scholar]
  86. MooRE, J.C.; Snider, R.J.; Robertson, L.S. Effects of different management practices on Collembola and Acarina in corn production systems. I: The effects of no-tillage and atrazine. Pedobiologia 1984, 26, 143–152. [Google Scholar]
  87. de Oliveira Filho, L.C.I.; Zeppelini, D.; Sousa, J.P.; Baretta, D.; Klauberg-Filho, O. Collembola community structure under different land management in subtropical Brazil. Ann. Appl. Biol. 2020, 177, 294–307. [Google Scholar] [CrossRef]
  88. Baretta, D.; Ferreira, C.S.; Sousa, J.P.; Cardoso, E.J.B.N. Colêmbolos (Hexapoda: Collembola) como bioindicadores de qualidade do solo em áreas com Araucaria angustifolia. Rev. Bras. De Ciência Do Solo 2008, 32, 2693–2699. [Google Scholar] [CrossRef]
  89. Pflug, A.; Wolters, V. Collembola communities along a European transect. Eur. J. Soil Biol. 2002, 38, 301–304. [Google Scholar] [CrossRef]
  90. Rożen, A.; Sobczyk, Ł.; Liszka, K.; Weiner, J. Soil faunal activity as measured by the bait-lamina test in monocultures of 14 tree species in the Siemianice common-garden experiment, Poland. Appl. Soil Ecol. 2010, 45, 160–167. [Google Scholar] [CrossRef]
  91. Schröder, P. Mesofauna. In Perspectives for Agroecosystem Management; Elsevier: Amsterdam, The Netherlands, 2008; pp. 293–306. [Google Scholar]
  92. Culik, M.P.; Martins, D.d.S.; Ventura, J.A. Collembola (Arthropoda: Hexapoda) communities in the soil of papaya orchards managed with conventional and integrated production in Espírito Santo, Brazil. Biota Neotrop. 2006, 6, No. 3. [Google Scholar] [CrossRef] [Green Version]
  93. Nakamoto, T.; Yamagishi, J.; Miura, F. Effect of reduced tillage on weeds and soil organisms in winter wheat and summer maize cropping on Humic Andosols in Central Japan. Soil Tillage Res. 2006, 85, 94–106. [Google Scholar] [CrossRef]
  94. Eisenhauer, N.; Milcu, A.; Sabais, A.C.; Bessler, H.; Brenner, J.; Engels, C.; Klarner, B.; Maraun, M.; Partsch, S.; Roscher, C. Plant diversity surpasses plant functional groups and plant productivity as driver of soil biota in the long term. PLoS ONE 2011, 6, e16055. [Google Scholar] [CrossRef]
  95. Milcu, A.; Partsch, S.; Langel, R.; Scheu, S. The response of decomposers (earthworms, springtails and microorganisms) to variations in species and functional group diversity of plants. Oikos 2006, 112, 513–524. [Google Scholar] [CrossRef]
  96. Bengtsson, J.; Ahnström, J.; WEIBULL, A.C. The effects of organic agriculture on biodiversity and abundance: A meta-analysis. J. Appl. Ecol. 2005, 42, 261–269. [Google Scholar] [CrossRef]
  97. Jeanneret, P.; Schüpbach, B.; Luka, H. Quantifying the impact of landscape and habitat features on biodiversity in cultivated landscapes. Agric. Ecosyst. Environ. 2003, 98, 311–320. [Google Scholar] [CrossRef]
  98. Heiniger, C.; Barot, S.; Ponge, J.-F.; Salmon, S.; Meriguet, J.; Carmignac, D.; Suillerot, M.; Dubs, F. Collembolan preferences for soil and microclimate in forest and pasture communities. Soil Biol. Biochem. 2015, 86, 181–192. [Google Scholar] [CrossRef]
  99. Querner, P.; Bruckner, A.; Drapela, T.; Moser, D.; Zaller, J.; Frank, T. Landscape and site effects on Collembola diversity and abundance in winter oilseed rape fields in eastern Austria. Agric. Ecosyst. Environ. 2013, 164, 145–154. [Google Scholar] [CrossRef]
  100. Pfingstmann, A.; Paredes, D.; Buchholz, J.; Querner, P.; Bauer, T.; Strauss, P.; Kratschmer, S.; Winter, S.; Zaller, J. Contrasting effects of tillage and landscape structure on spiders and springtails in vineyards. Sustainability 2019, 11, 2095. [Google Scholar] [CrossRef] [Green Version]
  101. Joimel, S.; Schwartz, C.; Hedde, M.; Kiyota, S.; Krogh, P.H.; Nahmani, J.; Pérès, G.; Vergnes, A.; Cortet, J. Urban and industrial land uses have a higher soil biological quality than expected from physicochemical quality. Sci. Total Environ. 2017, 584, 614–621. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Analysis of variance (ANOVA) barplot of the treatments data showing (A) abundance of strong furca (B) specie richness of strong furca (C) abundance of weak furca (D) specie richness of weak furca across three tillage types are differentiated by capitalized letters (A, B, AB) and across corn, and corn–soybean rotation is differentiated by lower case letters (a, b) Superscripts indicate significant differences between factors at the 0.05 level.
Figure 1. Analysis of variance (ANOVA) barplot of the treatments data showing (A) abundance of strong furca (B) specie richness of strong furca (C) abundance of weak furca (D) specie richness of weak furca across three tillage types are differentiated by capitalized letters (A, B, AB) and across corn, and corn–soybean rotation is differentiated by lower case letters (a, b) Superscripts indicate significant differences between factors at the 0.05 level.
Diversity 14 00994 g001aDiversity 14 00994 g001b
Figure 2. Analysis of variance (ANOVA) barplot of the treatments data showing (A) Shannon wiener index, (B) Simpson index, and (C) Pielou index of total Collembola across three tillage types (mould plow, ridge tillage, and No-tillage) in corn and corn–soybean rotation. Superscripts indicate significant differences between factors at the 0.05 level.
Figure 2. Analysis of variance (ANOVA) barplot of the treatments data showing (A) Shannon wiener index, (B) Simpson index, and (C) Pielou index of total Collembola across three tillage types (mould plow, ridge tillage, and No-tillage) in corn and corn–soybean rotation. Superscripts indicate significant differences between factors at the 0.05 level.
Diversity 14 00994 g002aDiversity 14 00994 g002b
Figure 3. Network of co-occurring Collembola species based on correlation analysis. Communities of six treatments, RT-MC (a), RT-CS (b), MP-MC (c), MP-CS (d), NT-MC (e), and NT-CS (f), included in this experiment to identify the different correlations between species. A strong (Spearman’s r = 0.6) and significant (p = 0.01) correlation is defined as a connection. The width of every connection connecting two nodes (edge) is proportional to the value of Spearman’s correlation coefficients; the size of every node is proportional to the relative abundance. 9 Collembola families color the co-occurring network.
Figure 3. Network of co-occurring Collembola species based on correlation analysis. Communities of six treatments, RT-MC (a), RT-CS (b), MP-MC (c), MP-CS (d), NT-MC (e), and NT-CS (f), included in this experiment to identify the different correlations between species. A strong (Spearman’s r = 0.6) and significant (p = 0.01) correlation is defined as a connection. The width of every connection connecting two nodes (edge) is proportional to the value of Spearman’s correlation coefficients; the size of every node is proportional to the relative abundance. 9 Collembola families color the co-occurring network.
Diversity 14 00994 g003aDiversity 14 00994 g003b
Table 1. Results of analysis of variance (ANOVA) of the effects of tillage type and crop rotation on the abundance and species richness of Collembola with strong furca and Collembola with weak furca.
Table 1. Results of analysis of variance (ANOVA) of the effects of tillage type and crop rotation on the abundance and species richness of Collembola with strong furca and Collembola with weak furca.
FactorsAbundanceRichness
Weak FurcaStrong FurcaWeak FurcaStrong Furca
FPFpFpFp
Til0.0240.9763.8910.0390.3240.7270.6560.531
Mag4.1980.0550.1640.691.3620.2580.0570.813
Til × Mag1.0620.3662.5650.1050.6760.5210.2480.783
Table 2. Results of analysis of variance (ANOVA) of the effects of tillage type and crop rotation on the Shanon index, Simpson index, and Pielou index of Collembola communities.
Table 2. Results of analysis of variance (ANOVA) of the effects of tillage type and crop rotation on the Shanon index, Simpson index, and Pielou index of Collembola communities.
FactorsShanonSimpsonPielou
FpFPFP
Til0.2620.7730.3640.70.090.914
Mag0.3550.5580.3850.5430.8870.359
Til × Mag0.4020.6750.6790.521.9680.169
Table 3. Network data for the communities of six treatments (RT-MC, RT-CS, MP-MC, MP-CS, NT-MC, NT-CS) were included in this experiment to identify the correlation between species.
Table 3. Network data for the communities of six treatments (RT-MC, RT-CS, MP-MC, MP-CS, NT-MC, NT-CS) were included in this experiment to identify the correlation between species.
TreatmentsNodes NumberEdges NumberModularity (MD)Clustering Coefficient (cc)Average Path Length (APL)Network Diameter (NA)Graph Density (GD)Average Degree (AD)
MT-MC-SR17220.4590.5561.530.1621.294
MT-CS-SR19200.6090.5512.41860.1171.053
RT-MC-SR25720.3410.6292.09140.242.88
RT-CS-SR22490.3970.542.51960.2122.227
NT-MC-SR19320.4820.5492.6250.1871.684
NT-CS-SR25460.4240.6222.80260.1531.84
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Habib, T.; Liu, S.; Chang, L.; Wu, Y.; Hao, C.; Wu, D. Investigating the Assemblages of Two Groups of Collembola (Strong Furca and Weak Furca) under Different Agricultural Management Systems, Northeastern China. Diversity 2022, 14, 994. https://doi.org/10.3390/d14110994

AMA Style

Habib T, Liu S, Chang L, Wu Y, Hao C, Wu D. Investigating the Assemblages of Two Groups of Collembola (Strong Furca and Weak Furca) under Different Agricultural Management Systems, Northeastern China. Diversity. 2022; 14(11):994. https://doi.org/10.3390/d14110994

Chicago/Turabian Style

Habib, Tayyiba, Shuchen Liu, Liang Chang, Yunga Wu, Cao Hao, and Donghui Wu. 2022. "Investigating the Assemblages of Two Groups of Collembola (Strong Furca and Weak Furca) under Different Agricultural Management Systems, Northeastern China" Diversity 14, no. 11: 994. https://doi.org/10.3390/d14110994

APA Style

Habib, T., Liu, S., Chang, L., Wu, Y., Hao, C., & Wu, D. (2022). Investigating the Assemblages of Two Groups of Collembola (Strong Furca and Weak Furca) under Different Agricultural Management Systems, Northeastern China. Diversity, 14(11), 994. https://doi.org/10.3390/d14110994

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop