- Plant Conservation and Population Biology, Department of Biology, KU Leuven, Leuven, Belgium
Conversion of natural ecosystems into agricultural land may strongly affect the soil microbiome and the functioning of the soil ecosystem. Alternative farming systems, such as organic farming, have therefore been advocated to reduce this impact, yet the outcomes of different agricultural management regimes often remain ambiguous and their evaluations mostly lack a proper more natural benchmark. We used high-throughput amplicon sequencing, linear models, redundancy analyses, and co-occurrence network analyses to investigate the effect of organic and integrated pest management (IPM) on soil fungal and bacterial communities in both the crop and drive rows of apple orchards in Belgium, and we included semi-natural grasslands as a benchmark. Fungi were strongly influenced by agricultural management, with lower diversity indices and distinct communities in IPM compared to organic orchards, whereas IPM orchards had a higher AMF abundance and the most complex and connected fungal communities. Bacterial diversity indices, community composition, and functional groups were less affected by management, with only a higher network connectivity and abundance of keystone taxa in organic drive rows. On the other hand, none of the agricultural soil microbiomes matched the complexity and connectedness of our semi-natural benchmark, demonstrating that even more nature-friendly agricultural management practices strongly affect the soil microbiome and highlighting the essential role of (semi-)natural systems as a harbor of robust and functionally diverse fungal and bacterial communities.
Introduction
Soils and their associated biodiversity play a crucial role in providing a wide range of ecosystem services in terrestrial ecosystems (Barrios, 2007; De Vries et al., 2013; van der Heijden and Wagg, 2013; Bardgett and Van Der Putten, 2014). Bacteria and fungi are essential parts of soil biodiversity, mediating fundamental ecological functions, such as nutrient cycling and carbon sequestration (Graham et al., 2016; Paul, 2016; Steffan and Dharampal, 2019), and directly influencing plant health through symbiotic and pathogenic interactions (Van Der Heijden et al., 2008; De Coninck et al., 2015). Ongoing global conversion of natural ecosystems into agricultural land poses a continuing threat to soil biodiversity and soil functioning (Orgiazzi et al., 2016), and due to the still increasing demand for food, feed, fiber, and bioenergy crops, millions of square kilometers of natural ecosystems are threatened by the expansion of agricultural land (Tilman et al., 2017; Williams et al., 2020). Soil biodiversity in agricultural lands is known to be continuously exposed to a variety of threats, including pesticide and fertilizer applications and frequent soil disturbance (Orgiazzi et al., 2016). Evidence is increasing that this results in an agricultural soil microbiome that is significantly less complex and highly different from that of more natural land uses (Creamer et al., 2016; Lammel et al., 2021).
In an effort to reduce the impact of agriculture on the environment and biodiversity, alternative farming systems, such as organic farming, are commonly advocated. Organic practices exclude the use of synthetic pesticides and fertilizers and rely more on natural processes for biological pest control and nutrient supply (Reganold and Wachter, 2016). Globally, organic farming is a booming sector with 72.3 million ha of organic farmland as of 2019, increasing more than 5-fold since 1999 (Schlatter et al., 2021). This not only indicates the need for a thorough comparison of soil microbiomes between organic and conventional farming systems, but also between farming systems and a benchmark of natural or semi-natural land. Organic farms generally harbor distinct microbial communities compared to conventional farms (Hartmann et al., 2015; Wang et al., 2016; Lupatini et al., 2017; Jurburg et al., 2020), although some studies did not find differences in bacterial community composition (Pershina et al., 2015; Hendgen et al., 2018). Soil microbiomes in organic farming systems have also been shown to have a higher total microbial biomass and activity compared to conventional systems (Esperschütz et al., 2007; Reilly et al., 2013; Lori et al., 2017; Martínez-García et al., 2018). Additionally, a higher microbial diversity has been found under organic management (Tuck et al., 2014; Hartmann et al., 2015). On the other hand, Jurburg et al. (2020) did not find differences in microbial richness between conventional and organic farms, while others found inconsistent responses of fungal and bacterial diversity (Hendgen et al., 2018; Lupatini et al., 2018; Finn et al., 2021).
Overall, the response of soil communities to agricultural management is complicated and cannot be captured by simple biodiversity indices (Hartmann et al., 2015). Due to their essential importance for ecosystem functioning, changes in the community composition and complexity of soil microbiomes can be expected to be more important than in species diversity alone (Wagg et al., 2019). The composition of functional groups is known to strongly impact soil functioning through the regulation of key ecosystem processes (Barrios, 2007). Arbuscular mycorrhizal fungi (AMF) for instance form a mutualistic relationship with a wide range of plant species and allocate nutrients to the plant in exchange for photosynthates (Van Der Heijden et al., 2008; Powell and Rillig, 2018). In addition, AMF can alter soil abiotic conditions and other plant biotic interactions, further influencing plant production and ecosystem functioning (Chaudhary, 2016; Powell and Rillig, 2018). Microbiomes are highly structured communities and form complex networks where microorganisms associate either directly or indirectly, in turn affecting functioning (Fuhrman, 2009; Banerjee et al., 2016). Creating co-occurrence networks can help elucidate these complex interactions and, together with network analysis, it allows both the quantification of association shifts in different environmental settings (Barberán et al., 2012; Layeghifard et al., 2017; de Vries et al., 2018) and the identification of keystone taxa (Herren and McMahon, 2018). A higher network complexity is expected to result in more resilient microbial communities, due to complementation and redundancy effects, and has been found to positively influence nutrient cycling (Wagg et al., 2019). Agricultural intensification has been shown to lower root fungal network complexity and the number of keystone taxa (Banerjee et al., 2019), whereas abandoned arable fields harbor more connected networks as ecological restoration progresses, improving carbon uptake (Morriën et al., 2017). Throughout the succession toward more natural ecosystems, soil microbiomes generally become more complex and stable and contain less soil pathogens and more mycorrhizal fungi (Kardol et al., 2006; Walker et al., 2010; Wagg et al., 2020).
So far, most research investigating effects of agricultural practices on soil microbial communities has been performed in arable systems, whereas perennial crop systems, such as orchards, have received less attention. Studies in orchards have mostly focused on a limited number of taxa, often AMF, and to our knowledge, the microbial network structure under different orchard management practices has not yet been investigated (Radić et al., 2014; Meyer et al., 2015; Meng et al., 2016; Van Geel et al., 2017). Furthermore, in order to be able to fully appreciate the effects of different agricultural management strategies on the soil microbiome and to quantify to what extent microbial communities have been degraded as a result of agricultural practices, it is highly informative to include a natural or semi-natural benchmark in the comparison (Katayama et al., 2019).
In this study, we used high-throughput amplicon sequencing to characterize soil fungal and bacterial communities in 30 Belgian apple orchards under conventional farming following Integrated Pest Management (IPM), and in 30 apple orchards under certified organic management. Furthermore, we added 30 semi-natural grasslands as a natural benchmark. Our objectives were to evaluate the effects of the agricultural management system on (i) taxon richness and diversity; (ii) phylogenetic diversity, (iii) community composition, and (iv) functional group abundance, all in comparison with the semi-natural reference system. Additionally, we aimed to (v) compare co-occurrence networks and their complexity between the management and land use types and to identify key microbial taxa within these networks.
Materials and Methods
Study Sites and Sampling
The sampled locations were located in Flanders, the northern region of Belgium, from 50°44’52.6” to 51°16’12.9“N and 3°39’14.2” to 5°24’47.5“E (Supplementary Figure S1). The region has a temperate oceanic climate with an average temperature of 9.8°C and an average annual precipitation of 925 mm. Soil samples were taken in October 2018 in 30 orchards managed under IPM and in 30 organically managed orchards. IPM orchards were fertilized using synthetic fertilizers at different time points during the year, whereas pests were managed using a combination of biological, mechanical, and physical methods, such as mating disruption, with chemical pesticides applied only as a last resort. All organic orchards were in agreement with European legislation [Council Regulation (EC) No 834/2007], meaning that they applied organic fertilizers and barred synthetic pesticides. Organic orchards instead also relied on non-chemical methods, such as the propagation of natural pest predators, complemented with organic pesticides, such as sulfur when necessary. All orchards contained Jonagold apple trees (Malus domestica “Jonagold”) at production age. In each orchard, three plots were randomly established at least 10 m from the orchard edge. In each plot, five soil samples (10 cm depth) were taken using a soil auger (2 cm diameter) from two 1 m2 sub-plots, one in the crop row, and one in the adjacent drive row, and pooled into one bulk sample per sub-plot. Crop rows are either sprayed with herbicides (IPM) or mechanically weeded (organic) to remove weeds, while drive rows are mown to maintain herbaceous vegetation. In addition, 30 reference semi-natural grasslands with similar soil type, hydrology, and slope were sampled in three 1 m2 plots in their center in a similar way. These unfertilized grasslands are managed by mowing or grazing for nature conservation purposes. This resulted in 90 soil samples from crop rows and 90 samples from drive rows in both the organic and IPM orchards (360 samples in total) and 90 samples from semi-natural grasslands. After mixing, the soil samples were divided in two parts, one part was stored at 5°C for chemical analyses, and the other part was stored at −20°C for DNA extraction.
Soil Chemical Analyses
Soil water content was determined by drying soil at 50°C for 24 h or until a stable weight was achieved. To quantify organic matter content, the dried soil was subsequently put in a muffle furnace at 600°C to combust organic compounds. Soil pH was measured using a glass electrode in a 1:10 soil/water mixture. We determined plant-available nitrogen content in the soil, under the form of ammonium and nitrate availability, by shaking 10 g of soil in 200 ml of 1 M potassium chloride solution for 1 h. Subsequent colorimetric analysis was performed using a segmented flow auto analyzer (Skalar, Breda, the Netherlands). As a measure of plant-available phosphorus in the soil, Olson P was quantified by shaking 2 g of dried soil for 30 min with a 0.5 M sodium bicarbonate solution at pH 8.5, followed by colorimetric analysis using the molybdenum blue method (Lajtha et al., 1999).
DNA Extraction, PCR Amplification, and High-Throughput Amplicon Sequencing
Genomic DNA was extracted from 250 mg of soil using the soil DNA Isolation Plus Kit (Norgen Biotek Crop., Thorold, ON, Canada) following the manufacturer’s instructions. The DNA extracts were diluted 10 times before PCR amplification. Fungal communities were profiled through PCR amplification targeting the ITS2 region of the rRNA gene using the sample-specific barcode-labeled versions of the primers ITS86F and ITS4 (Waud et al., 2014; dual-index sequencing strategy, Kozich et al., 2013). Bacterial communities were characterized targeting the hypervariable V4 region of the 16S rRNA gene with the 515F/806R primer pair (Caporaso et al., 2011). PCR was performed in a 25 μl volume, containing 0.25 mM of each dNTP, 0.5 μM of each primer, ×1 ALLin HiFi Buffer, 1 U ALLin HiFi DNA Polymerase (HighQu, Kraichtal, Germany), and 1 μl genomic DNA. DNA samples were denatured at 95°C for 1 min, followed by 35 cycles of 15 s at 95°C, 15 s at 60°C, and 15 s at 72°C, and finally 3 min at 72°C. Amplicons were purified using the Agencourt AMPure XP kit (Beckman Coulter Life Sciences, Indianapolis, IN, United States). Next, purified double-stranded DNA (dsDNA) amplicons were quantified using the Quant-iT PicoGreen dsDNA Assay Kit and Qubit fluorometer (Invitrogen), and pooled in equimolar quantities in four amplicon libraries. These amplicon libraries were loaded on an agarose gel, and the amplicon of the expected size was excised and purified using the QIAquick Gel Extraction Kit (Qiagen). After dilution to 2 nM, the amplicons were sequenced at the Genomics Core UZ Leuven on an Illumina Miseq platform using 2 × 250 paired-end Illumina chemistry with v2 500 cycle reagent kit (Illumina, San Diego, CA, United States).
Bioinformatics
Quality-filtered sequences were clustered into operational taxonomical units (OTUs) following the recommended pipeline in USEARCH (v. 11; Edgar, 2013). First, paired-end reads from the obtained demultiplexed FASTQ files were merged to create consensus sequences (fastq_mergepairs). Next, reads were quality filtered (fastq_filter), allowing a maximum expected error of 0.5 for the individual sequences. Sequences were then dereplicated and sorted by abundance (fastx_uniques). To improve the accuracy of diversity estimates, singletons were removed before clustering (Brown et al., 2015). Finally, sequences were clustered into OTUs defined at 97% sequence similarity (cluster_otus), during which chimeric sequences were also removed (Edgar, 2013).
Data Analyses
All analyses were performed using R (R Core Team, 2021, v. 4.0.0) and separately for the fungal and the bacterial datasets. More than 20,100 fungal and 36,000 bacterial effective sequences were obtained on average for each sample. Sufficient sequencing depth was assessed and confirmed for all samples through rarefaction using the vegan package (Oksanen et al., 2020). To avoid bias due to differences in sequencing depth, reads were randomly resampled from each sample to an equal number of reads (10.000 for fungi and 20.000 for bacteria). Statistical analyses were done in two steps. First, all data from the orchards were compiled in order to assess the effects of orchard management (IPM vs. organic) and location (crop row vs. drive row), further referred to as management models. Second, to model the effects of land use (IPM vs. organic. vs. semi-natural grassland) on the microbiome, the data from the herbaceous drive rows of both orchard types and the data from the reference semi-natural grasslands were used. We omitted the crop rows from this analysis in order to obtain a more conservative estimate of the effects of agricultural management as compared to the semi-natural benchmark. These models are further referred to as land use models. A set of covariates, i.e., soil characteristics (water content, organic matter, pH, ammonium, nitrate, and P) and spatial predictors (distance-based Moran’s eigenvector maps, dbMEMs) were added to all models. dbMEMs were included to account for possible effects from the spatial structure of the study sites (Legendre and Legendre, 2012).
Operational taxonomical unit taxonomic diversity was quantified by applying the Hill number framework, using both diversity of order 0 (OTU richness) and order 1 (1D; Hill, 1973), and phylogenetic diversity (PD) was calculated using Faith’s PD (Faith, 1992). A generalized linear mixed model (GLMM) following a Poisson distribution and log link function was used to test effects of land use and management and the variables described above on OTU richness, while linear mixed models (LMM) were used to test effects on 1D and PD. To assess effects on functional groups, we assigned functionality using the FUNGuild for fungal OTUs (Nguyen et al., 2016) and METAGENASSIST for bacterial OTUs (Arndt et al., 2012). The read count of fungal saprotrophs, AMF, and plant pathogens as well as the four most dominant metabolic groups of bacteria (sulfite reducer, ammonia oxidizer, dehalogenation, and xylan degrader) were used as response variables in individual GLMMs following a Poisson distribution and log link function. In all management models, sub-plot was nested within plot as a random factor, whereas plot was included as a random factor in all land use models. Model selection was based on the Akaike’s information criterion (AIC) with a backward selection procedure. Variance inflation factors (VIFs) were calculated to detect multicollinearity between the explanatory variables and were low (all <3). Next, we explored variation in fungal and bacterial community composition through canonical redundancy analyses (RDA), using the vegan package (Oksanen et al., 2020) with management and location (management models) or land use (land use models), soil characteristics, and dbMEM variables as explanatory variables. Forward model selection was done with 999 permutations.
To identify specific OTUs that were associated with a specific land use and location, we performed correlation-based association analyses using the multipatt function of the indicspecies package (De Caceres et al., 2016). For visualization, all OTUs with significant associations (p < 0.05) were displayed as target nodes in a network with the five land use and location types as source nodes. Networks were visualized using the spring-embedded layout algorithm in Cytoscape (Shannon et al., 2003). We then investigated co-occurrence patterns by constructing networks for each land use type separately. Interactions within the networks consisted of positive and significant Spearman’s correlations (ρ > 0.6 and p < 0.05). All network topology parameters were calculated using the igraph package in R (Csardi and Nepusz, 2006), whereas the networks were visualized in Gephi (Bastian et al., 2009). The nodes, corresponding with OTUs, are the fundamental units in these networks, while the edges form the connections between these nodes. The degree of a node then represents the total number of edges associated with that node. Network diameter is the largest distance between two nodes in a network, and average path length is the mean of the shortest distance between all pairs of network nodes in the network. Closeness centrality quantifies the inverse distance to all other nodes in the network and reveals the importance of nodes for the dissemination of information. Betweenness centrality on the other hand indicates the level in which nodes act as a bridge connecting different clusters of the graph. Modularity expresses the strength of division within a network into clusters. The OTUs with the highest degree and highest closeness centrality and the lowest betweenness centrality were identified as keystone taxa (Berry and Widder, 2014). OTUs with degree > 5, closeness centrality > 0.005, and betweenness centrality < 40 were selected as keystone taxa for the fungal networks, and OTUs with degree > 15, closeness centrality > 0.03, and betweenness centrality < 50 for the bacterial networks.
Results
Diversity Indices
In the management model, fungal OTU richness was significantly affected by management (p = 0.0004), location (p < 0.0001), and by their interaction (p < 0.0001), with a higher richness in organic orchards and in drive rows. The difference between crop row and drive row was smaller in organic orchards (Figure 1A; Supplementary Table S1). In contrast, fungal 1D was only affected by location in the management model (p < 0.0001), with a higher 1D in drive rows than in crop rows (Figure 1A; Supplementary Table S1). Fungal PD was significantly higher in organic orchards (p = 0.0084) and in drive rows (p = 0.0001) in the management model. Bacterial OTU richness was not significantly affected by management or location in the management model, but bacterial 1D and PD were significantly higher in crop rows compared to drive rows (both p < 0.0001; Figure 1B; Supplementary Table S1). In the land use models, fungal OTU richness was significantly lower in grasslands compared to organic orchards (p = 0.015), whereas grasslands and organic orchards both did not significantly differ in fungal OTU richness compared to IPM orchards (Figure 1A; Supplementary Table S1). Fungal 1D and PD were unaffected by land use (Figure 1A; Supplementary Table S1). All bacterial diversity indices were significantly lower in grasslands compared to both orchard types (p < 0.0001, p = 0.0003, and p = 0.0004 for OTU richness, 1D, and PD, respectively; Figure 1B; Supplementary Table S1). Overall, the significant differences that were found between management types for the studied diversity indices in the management models disappeared in the land use models.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/830668/fmicb-13-830668-HTML/image_m/fmicb-13-830668-g001.jpg)
Figure 1. Diversity indices in integrated pest management (IPM) and organic orchards with a semi-natural grassland reference. Modeled mean values ±1SE and raw data points of Hill numbers of order 0 [operational taxonomical unit (OTU) richness] and order 1 (1D) and phylogenetic diversity (PD) for fungal (A) and bacterial (B) communities across land use types. When land use or location was not selected during the model selection procedure, only raw data points are shown. Left panels: crop rows (dark) and drive rows (light) of IPM and organic orchards. Right panels: drive rows of IPM and organic orchards with a semi-natural grassland reference. Same lowercase letter indicates no statistically significant (p < 0.05) difference between land uses.
Community Composition
The management RDA selected management, location, and their interaction as significantly explaining variation in fungal orchard communities (Figure 2A; Supplementary Table S2). Selected covariates included phosphorus, organic matter, and soil pH, whereas spatial vectors also significantly explained community variation. The final model explained 17% of the variation. The land use RDA selected land use as a factor significantly explaining fungal community composition together with covariates phosphorus, water content, soil pH, and spatial vectors (Figure 2A; Supplementary Table S2). The final model explained 16% of the variation. Variation in bacterial community composition in orchards (land use model) was significantly affected by the factor location while phosphorus, water content, organic matter, and soil pH and spatial vectors were selected as covariates (Figure 2B; Supplementary Table S2). The final model explained 16% of the variation. In the land use RDA, land use was selected as a factor explaining variation in bacterial community composition together with covariates soil pH, phosphorus, water content, organic matter, nitrate, and spatial vectors (Figure 2B; Supplementary Table S2). The final model explained 30% of the variation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/830668/fmicb-13-830668-HTML/image_m/fmicb-13-830668-g002.jpg)
Figure 2. Redundancy analysis (RDA) ordination biplots of communities in IPM and organic orchards with a semi-natural grassland reference. RDA biplots generated of fungal (A) and bacterial (B) communities in crop rows and drive rows of IPM and organic orchards (upper plots) and for drive rows of IPM and organic orchards with a semi-natural grassland reference (lower plots). Arrows indicate variables selected in the model selection procedure (Supplementary Table S2).
Functional Groups
The number of saprotrophic fungal reads in the management model was not influenced by management or location (Figure 3A; Supplementary Table S3). The abundance of AMF was influenced by management (p = 0.045), location (p < 0.0001), and their interaction (p < 0.0001), with a higher abundance in IPM and in drive rows, though the difference between drive rows and crop rows was larger in organic orchards. The pathogen abundance was significantly affected by location (p < 0.0001) and the interaction between management and location (p < 0.0001; Figure 3; Supplementary Table S3), with more pathogen OTUs in crop rows compared to drive rows in IPM orchards. In the land use model, a significantly higher number of saprotroph reads was found in grasslands compared to both orchard types (p = 0.0002; Figure 3A; Supplementary Table S3). The number of AMF reads was significantly lower in organic orchards compared to IPM and grasslands (p = 0.024; Figure 3B; Supplementary Table S3). The number of plant pathogen reads was lower in grasslands compared to both orchard types (p = 0.013; Figure 3C; Supplementary Table S3). Bacterial functional groups were less affected by management and land use than fungal functional groups and only the number of ammonia oxidizer reads in the land use model was significantly affected by land use (p < 0.0001; Supplementary Figure S2; Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/830668/fmicb-13-830668-HTML/image_m/fmicb-13-830668-g003.jpg)
Figure 3. Read count of fungal functional groups in IPM and organic orchards with a semi-natural grassland reference. Modeled mean values ±1SE and raw data points of saprotrophic fungi, arbuscular mycorrhizal fungi (AMF), and plant pathogens across land use types. When land use or location was not selected during the model selection procedure, only raw data points are shown. Left panels: crop rows (dark) and drive rows (light) of IPM and organic orchards. Right panels: drive rows of IPM and organic orchards with a semi-natural grassland reference. Same lowercase letter indicates no statistically significant (p < 0.05) difference between land uses.
Indicator Taxa
In the fungal network, 258 OTUs significantly associated with one or more of the three land use types. Of these OTUs, 63.9% were associated with just one land use type and almost half of these taxa (30.6%) were associated with IPM crop rows. These taxa included a higher number of Glomeromycota, in addition to mainly Ascomycota and Basidiomycota. Both IPM orchard locations combined shared more OTUs with semi-natural grasslands (17.1%) than did organic orchards (3.1%; Figure 4A). In the bacterial networks, a total of 870 OTUs were significantly associated with at least one land use type, with 42.3% of these OTUs associated with one particular type and a fairly even distribution of the number of habitat specialists and their assigned phyla across land use types. However, we found more Acidobacteria habitat specialists in organic and more Planctomycetes habitat specialists in IPM orchards. There was more shared indicator bacterial OTUs between organic orchards and grasslands (19.4%) than between IPM orchards and grasslands (7.7%; Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/830668/fmicb-13-830668-HTML/image_m/fmicb-13-830668-g004.jpg)
Figure 4. Correlation-based association networks across land use types. Networks are constructed at the OTU level of the fungal (A) and bacterial (B) communities across IPM and organic crop rows and drive rows with a semi-natural grassland reference. Only statistically significant OTU associations are displayed as source nodes (p < 0.05).
Co-occurrence Networks
We separately evaluated co-occurrence networks for all land use types involved in this study, based on significant correlations, revealing notable differences in the topology between networks. Both for fungal and bacterial networks, the grassland network was the most densely populated, complex and connected compared to all orchard networks, and it contained the highest number of keystone taxa (Figure 5; Supplementary Table S5). The fungal grassland network contained 31 keystone taxa, mainly belonging to the phyla Ascomycota and Basidiomycota. Interestingly, the fungal IPM crop row network was the most complex of all fungal orchard networks (Figure 5; Supplementary Table S5). It contained eight keystone taxa belonging to the phyla Ascomycota and Glomeromycota. The other fungal networks were less complex and connected and contained less keystone taxa; the organic drive row network had three keystone taxa, while the organic crop row network had one, and the IPM drive row network contained no keystone taxa at all (Figure 5; Supplementary Table S5). Bacterial networks were more complex and contained a higher number of nodes than fungal networks. There was a high abundance of keystone taxa in the bacterial grassland network, with a total of 110 taxa belonging to a variety of phyla—mainly Chloroflexi, Planctomycetes, and Proteobacteria (Figure 5). Of the bacterial orchard networks, the organic drive row network seems to be best connected, with a degree equal to the grassland network, a relatively high closeness centrality, high betweenness centrality, and low average path length (Figure 5; Supplementary Table S5). It contained 22 keystone taxa, most belonging to Proteobacteria, Actinobacteria, and Chloroflexi. The other orchard networks showed similar topology values (Figure 5; Supplementary Table S5). The organic crop row and IPM drive row each had one keystone taxon of the phylum Proteobacteria. The bacterial IPM crop row network did not contain any keystone taxa.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/830668/fmicb-13-830668-HTML/image_m/fmicb-13-830668-g005.jpg)
Figure 5. Co-occurrence networks across land use types. Networks are generated for fungal (upper plots) and bacterial (lower plots) communities across IPM and organic crop rows and drive rows with a semi-natural grassland reference. Only interactions with positive and significant Spearman’s correlations are displayed (p > 0.6 and p < 0.05), with OTUs as nodes and interactions as edges. Large nodes indicate the keystone taxa.
Discussion
To our knowledge, this is the first study to evaluate the effects of IPM and organic management in a perennial agricultural system on the soil microbiome in reference to a semi-natural benchmark. We generally found significant effects of orchard management on fungal diversity indices, community composition, the abundance of functional groups, and the complexity of co-occurrence networks. In contrast, whereas the complexity of bacterial networks was affected by orchard management, bacterial diversity, community composition, and the abundance of functional groups were not. On the other hand, differences in microbial diversity between the drive rows of both orchard management systems were very small and disappeared when compared to a semi-natural benchmark. None of the agricultural soil microbiomes matched the complexity and connectedness of the semi-natural benchmark.
Management and Location Effects on Soil Fungi
Fungal communities were significantly different and OTU richness and phylogenetic diversity were higher in organic orchards compared to IPM orchards, and in drive rows compared to crop rows. However, the higher OTU richness found in organic orchards compared to IPM orchards was more noticeable in crop rows and less so in drive rows. Previous work in perennial agricultural systems did not find management effects on fungal OTU richness, yet they did find management effects on fungal community composition (Hendgen et al., 2018; Jurburg et al., 2020). An important driver may be the use of herbicides in the IPM crop rows, which prevents the emergence of herbaceous plant species and might affect plant species diversity in the adjacent drive row through spray drift as well. On the other hand, mechanical weeding in organic crop rows allows the emergence of some herbs during part of the growing season. As plant species richness has been shown to positively correlate with soil fungal diversity, this could explain the lower fungal richness, phylogenetic diversity, and altered community composition in the IPM orchards (Vukicevich et al., 2016). All investigated fungal diversity indices were higher in drive rows compared to crop rows, confirming the findings of Hendgen et al. (2018) in vineyards. The herbaceous drive rows act as a permanent vegetation cover where the rhizosphere delivers a suitable habitat for soil microbes (Berg and Smalla, 2009). The root exudates supply nutrient sources, such as organic carbon and amino acids, and anti-pathogenic compounds, as well as take part in plant-microorganism interactions through the production of signaling molecules, e.g., acting as chemoattractant to microbes or microbial growth promotors (Haichar et al., 2014). Additionally, frequent mowing and litter decomposition in the drive rows increase carbon inputs to the soil, acting as a nutrient source for microbes. Indeed, we found higher organic matter content in drive rows compared to their adjacent crop row (Supplementary Table S6).
General fungal richness aside, understanding the management response of key fungal functional groups within the soil microbiome is critical as they strongly impact soil functioning (Barrios, 2007; Lüneberg et al., 2019; Lekberg et al., 2021). Both AMF abundance and the prevalence of Glomeromycota in the indicator species analysis and as keystone taxa in the co-occurrence networks were higher in the IPM than in the organic orchards. This is in line with previous studies finding either lower (Di Giacinto et al., 2020) or similar AMF occurrences (Turrini et al., 2017) in organic compared to IPM orchards. Whereas the response of AMF to herbicide applications has been reported to be highly variable, AMF may negatively respond to mechanical weeding in organic crop rows, likely explaining the increased occurrence of AMF in IPM orchards (Turrini et al., 2017; Mandl et al., 2018; Hage-Ahmed et al., 2019). On the other hand, we found no management effects on the number of plant pathogen reads, aligning with recent work (Hannula et al., 2021). While AMF were more abundant in drive rows, pathogenic fungi were more abundant in crop rows. This is in line with other work that has shown that cover crop diversity increased AMF root colonization, likely due to the increased availability of root carbon sources, whereas it reduced the amount of soil-borne pathogens, possibly because of the antagonistic role of AMF (Haichar et al., 2014; Vukicevich et al., 2016; Bowles et al., 2017).
Our fungal networks were relatively similar in both management types and across locations in terms of complexity and connectivity, with a notable exception for the more complex and connected IPM crop row network. It is known that the complexity of microbial co-occurrence networks in the soil is reduced with increasing disturbance and compaction (Hartmann et al., 2014; Karimi et al., 2017; Xia et al., 2022). However, the soil in IPM crop rows is not mechanically disturbed, as opposed to organic crop rows, which are mechanically weeded, or drive rows, which are frequently compacted by heavy machinery. The increased complexity of our IPM crop row network could thus be attributed to the relatively undisturbed soil under IPM crop rows. This is contrary to its relatively low observed fungal richness and phylogenetic diversity, highlighting the limited usefulness of simple diversity indices when quantifying management impacts (Hartmann et al., 2015; Karimi et al., 2017).
Management and Location Effects on Soil Bacteria
Unlike fungi, orchard management did not affect bacterial diversity indices, nor bacterial community composition. This lack of management effects has previously also been observed for bacterial diversity indices in coffee plantations (Jurburg et al., 2020) and for bacterial community composition in vineyards (Hendgen et al., 2018). Bacteria have indeed been found to be less sensitive to land use, as they respond more directly to soil properties (Petersen et al., 2019). As soil characteristics were mostly accounted for in our analyses through integrating them in the models as covariates, this might explain the lack of a management effect on bacterial diversity and community composition. It is for instance known that soil bacterial diversity and community composition are heavily influenced by soil pH levels and organic matter content, as observed here (Fierer and Jackson, 2006; Rousk et al., 2010; Griffiths et al., 2011; Yashiro et al., 2016; Karimi et al., 2017). In addition, the spatial autocorrelation identified here has also been detected in other systems, indicating the importance of spatial processes in the distribution of bacterial communities (Constancias et al., 2015; Shi et al., 2015). Fungi are also known to be more sensitive to plant diversity loss than bacteria due to their often more specific interactions with plant species, both through root exudates and litter inputs (Eisenhauer et al., 2017; Jurburg et al., 2020). On the other hand, bacterial diversity and phylogenetic diversity were higher in crop rows, which also harbored distinct bacterial communities compared to drive rows. Bacteria have a relatively high resistance to herbicides and physical soil disturbances, which could have led to the higher diversity indices found in crop rows (Schlatter et al., 2017; Chou et al., 2018; Mandl et al., 2018).
The bacterial co-occurrence networks were generally relatively similar across both management types, which were also found in an arable long-term trial cultivating wheat and rapeseed (Finn et al., 2021). However, the organic drive row network was more complex and contained more keystone taxa than the other orchard networks. The increased complexity in the organic drive rows might be due to the increased plant species diversity observed in its permanent herbaceous vegetation crop (Vukicevich et al., 2016). Additionally, bacterial networks have been found to have an increased number of positive links under a cover crop, indicating decreased competition and increased cooperation, potentially due to greater niche separation (Zheng et al., 2018). Again, we observed a disparity between the diversity indices responses and the network analyses. Whereas nutrient additions and soil disturbances in crop rows homogenize microhabitats in the soil, potentially resulting in less complex networks, they increased the overall nutrient levels, possibly leading to higher diversity indices (Gossner et al., 2016; Szoboszlay et al., 2017).
Land Use Effects on the Soil Microbiome
When comparing the effects of orchard management on microbial diversity indices in the drive rows to the semi-natural grassland benchmark, the latter had a lower fungal OTU richness and a lower value for all bacterial diversity indices compared to both orchard types. There was also a clear difference in the fungal community composition among all three land use types, though there was a higher similarity between both orchard systems. Additionally, bacterial community composition did not differ between orchard types, while there was a clear separation from the grassland communities. Other studies have reported both an increased fungal and bacterial diversity and a community shift in agricultural systems compared to more natural systems (Szoboszlay et al., 2017; Terrat et al., 2017; George et al., 2019; Petersen et al., 2019; Lammel et al., 2021). This has been attributed to the higher fertilizer inputs, increased soil compaction, and more frequent soil disturbance transferring nutrients into the soil and increasing microbial diversity indices (Hartmann et al., 2014; Szoboszlay et al., 2017; Longepierre et al., 2021). Additionally, disturbance might induce dormancy in some microbial taxa, or reduce dominant OTUs, opening opportunities for a range of weaker competitors and thus increasing diversity (Wang et al., 2014; Szoboszlay et al., 2017; George et al., 2019). These land use effects were more apparent for bacteria compared to fungi, which is in line with similar findings in tropical agricultural systems (Petersen et al., 2019; Lammel et al., 2021).
We found functionally enhanced communities in semi-natural grasslands compared to the herbaceous drive rows in both orchard types. The higher saprotroph read counts found in the semi-natural grasslands compared to the drive rows in both orchard types might be attributed to the much higher plant diversity found in these semi-natural grasslands, which changes the quality and diversity of litter inputs to the soil, affecting these primary decomposers (Millard and Singh, 2010; Francioli et al., 2021). Surprisingly, the number of AMF reads in organic drive rows was lower than in IPM drive rows, while we found no difference in AMF abundance between IPM and the benchmark grasslands. This is in contrast to reduced AMF abundances found under synthetic fertilization compared to a non-fertilized control and reduced AMF richness in grasslands compared to arable systems (Hannula et al., 2021; Zhu et al., 2021). Instead, other factors may influence AMF abundance in drive rows. It has for instance been suggested that compaction is higher in organic drive rows, potentially due to increased traffic as a result of higher labor inputs and the frequent need for fungicide applications under wet conditions, which may reduce mycorrhizal formation (Entry et al., 2002; Vogeler et al., 2006).
Both the fungal and bacterial semi-natural grassland co-occurrence networks were distinctly more complex, more connected and contained more keystone taxa than any of the orchard networks. In addition, the bacterial grassland network contained nodes from a much wider range of phyla in comparison with orchard networks. The increased abundance and diversity of keystone taxa in the semi-natural grassland networks further amplified network complexity (Herren and McMahon, 2018). As we only included positive correlations, links between network nodes are associated with potential co-dependent relationships, such as biotic interactions and niche sharing (Karimi et al., 2017; Xue et al., 2022). The tightened positive interactions found in both grassland networks could be related to improved trophic interactions and cooperation, as increased nutrient levels may displace the need for robust trophic interactions between different microbial functional groups (Holtkamp et al., 2011; Wu et al., 2021). Additionally, network complexity in orchards might be reduced by increased disturbance, compaction, and fertilization due to the homogenization of available soil niches (Fraterrigo et al., 2005).
Agricultural management substantially affects soil microbiomes and the effects of conventional and organic management have frequently been evaluated. However, to adequately quantify the importance of these effects, it is important to include a (semi-)natural benchmark. This study found effects of both IPM and organic agricultural management on the soil microbiome in the crop and drive rows of apple orchards. Our results indicate the importance of management practices that minimize both soil disturbance and compaction and favor a more species-rich permanent vegetation cover for maintaining more complex and connected microbiomes in perennial agricultural systems. However, the agricultural soil microbiomes in the drive rows were no match for the complexity and connectedness found in the semi-natural benchmark, demonstrating that even more nature-friendly agricultural management practices strongly affect the soil microbiome and highlighting the need for the conservation of (semi-)natural systems to harbor robust and functionally diverse fungal and bacterial communities.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA786545.
Author Contributions
EH: conceptualization and design, field work, data analysis and interpretation, drafting the manuscript, revision, and editing. GP: laboratory analysis, manuscript revision, and editing. OH: conceptualization and design, manuscript revision and editing, and supervising the research. All authors contributed to the article and approved the submitted version.
Funding
This project was funded by Research Foundation—Flanders (71635) and by the KU Leuven Special Research Fund (BOF; project C24/18/034).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We wish to thank all apple growers and land managers involved for access to their land as well as Proefcentrum Fruitteelt vzw for helping us contact them. We are also grateful to Kasper Van Acker and Michiel Stas for assistance during field work and Dries Grauwels and Elvira Vassilieva for performing soil analyses. We also thank Maarten Van Geel and Kenny Helsen for discussions regarding data analyses.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.830668/full#supplementary-material
References
Arndt, D., Xia, J., Liu, Y., Zhou, Y., Guo, A. C., Cruz, J. A., et al. (2012). METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 40, 88–95. doi: 10.1093/nar/gks497
Banerjee, S., Kirkby, C. A., Schmutter, D., Bissett, A., Kirkegaard, J. A., and Richardson, A. E. (2016). Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biol. Biochem. 97, 188–198. doi: 10.1016/j.soilbio.2016.03.017
Banerjee, S., Walder, F., Büchi, L., Meyer, M., Held, A. Y., Gattinger, A., et al. (2019). Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. ISME J. 13, 1722–1736. doi: 10.1038/s41396-019-0383-2
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Bardgett, R. D., and Van Der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi: 10.1038/nature13855
Barrios, E. (2007). Soil biota, ecosystem services and land productivity. Ecol. Econ. 64, 269–285. doi: 10.1016/j.ecolecon.2007.03.004
Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: an open source software for exploring and manipulating networks.” in International AAAI Conference on Weblogs and Social Media. Association for the Advancement of Artificial Intelligence. 361–362.
Berg, G., and Smalla, K. (2009). Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 68, 1–13. doi: 10.1111/j.1574-6941.2009.00654.x
Berry, D., and Widder, S. (2014). Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. 5:219. doi: 10.3389/fmicb.2014.00219
Bowles, T. M., Jackson, L. E., Loeher, M., and Cavagnaro, T. R. (2017). Ecological intensification and arbuscular mycorrhizas: a meta-analysis of tillage and cover crop effects. J. Appl. Ecol. 54, 1785–1793. doi: 10.1111/1365-2664.12815
Brown, S. P., Veach, A. M., Rigdon-Huss, A. R., Grond, K., Lickteig, S. K., Lothamer, K., et al. (2015). Scraping the bottom of the barrel: are rare high throughput sequences artifacts? Fungal Ecol. 13, 221–225. doi: 10.1016/j.funeco.2014.08.006
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U. S. A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Chaudhary, V. B. (2016). Data descriptor: MycoDB, a global database of plant response to mycorrhizal fungi. Nature 3:160028. doi: 10.1038/sdata.2016.28
Chou, M. Y., Vanden Heuvel, J., Bell, T. H., Panke-Buisse, K., and Kao-Kniffin, J. (2018). Vineyard under-vine floor management alters soil microbial composition, while the fruit microbiome shows no corresponding shifts. Sci. Rep. 8, 11039–11039. doi: 10.1038/s41598-018-29346-1
Constancias, F., Saby, N. P. A., Terrat, S., Dequiedt, S., Horrigue, W., Nowak, V., et al. (2015). Contrasting spatial patterns and ecological attributes of soil bacterial and archaeal taxa across a landscape. Microbiology 4, 518–531. doi: 10.1002/mbo3.256
Creamer, R. E., Hannula, S. E., Leeuwen, J. P. V., Stone, D., Rutgers, M., Schmelz, R. M., et al. (2016). Ecological network analysis reveals the inter-connection between soil biodiversity and ecosystem function as affected by land use across Europe. Appl. Soil Ecol. 97, 112–124. doi: 10.1016/j.apsoil.2015.08.006
Csardi, G., and Nepusz, T. (2006). The igraph software package for complex network research. Int. J. Complex Syst. 1695, 1–9.
De Coninck, B., Timmermans, P., Vos, C., Cammue, B. P. A., and Kazan, K. (2015). What lies beneath: belowground defense strategies in plants. Trends Plant Sci. 20, 91–101. doi: 10.1016/j.tplants.2014.09.007
de Vries, F. T., Griffiths, R. I., Bailey, M., Craig, H., Girlanda, M., Gweon, H. S., et al. (2018). Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9:3033. doi: 10.1038/s41467-018-05516-7
De Vries, F. T., Thébault, E., Liiri, M., Birkhofer, K., Tsiafouli, M. A., Bjørnlund, L., et al. (2013). Soil food web properties explain ecosystem services across European land use systems. Proc. Natl. Acad. Sci. U. S. A. 110, 14296–14301. doi: 10.1073/pnas.1305198110
Di Giacinto, S., Friedel, M., Poll, C., Döring, J., Kunz, R., and Kauer, R. (2020). Vineyard management system affects soil microbiological properties. Oeno One 54, 131–143. doi: 10.20870/oeno-one.2020.54.1.2578
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Eisenhauer, N., Lanoue, A., Strecker, T., Scheu, S., Steinauer, K., Thakur, M. P., et al. (2017). Root biomass and exudates link plant diversity with soil bacterial and fungal biomass. Sci. Rep. 7, 1–8. doi: 10.1038/srep44641
Entry, J. A., Rygiewicz, P. T., Watrud, L. S., and Donnelly, P. K. (2002). “Arbuscular mycorrhizal response to adverse soil conditions,” Arbuscular Mycorrhizae: Interactions in Plants, Rhizosphere and Soils. eds. Sharma, A., and Johri, B. (Science Publishers), 135–158.
Esperschütz, J., Gattinger, A., Mäder, P., Schloter, M., and Fließbach, A. (2007). Response of soil microbial biomass and community structures to conventional and organic farming systems under identical crop rotations. FEMS Microbiol. Ecol. 61, 26–37. doi: 10.1111/j.1574-6941.2007.00318.x
Faith, D. P. (1992). Conservation evaluation and phylogenetic diversity. Biol. Conserv. 61, 1–10. doi: 10.1016/0006-3207(92)91201-3
Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 103, 626–631. doi: 10.1073/pnas.0507535103
Finn, D. R., Lee, S., Lanzén, A., Bertrand, M., Nicol, G. W., and Hazard, C. (2021). Cropping systems impact changes in soil fungal, but not prokaryote, alpha-diversity and community composition stability over a growing season in a long-term field trial. FEMS Microbiol. Ecol. 97, 1–17. doi: 10.1093/femsec/fiab136
Francioli, D., van Rijssel, S. Q., van Ruijven, J., Termorshuizen, A. J., Cotton, T. E. A., Dumbrell, A. J., et al. (2021). Plant functional group drives the community structure of saprophytic fungi in a grassland biodiversity experiment. Plant Soil 461, 91–105. doi: 10.1007/s11104-020-04454-y
Fraterrigo, J. M., Turner, M. G., Pearson, S. M., and Dixon, P. (2005). Effects of past land use on spatial heterogeneity of soil nutrients in southern Appalachian forests. Ecol. Monogr. 75, 215–230. doi: 10.1890/03-0475
Fuhrman, J. A. (2009). Microbial community structure and its functional implications. Nature 459, 193–199. doi: 10.1038/nature08058
George, P. B. L., Lallias, D., Creer, S., Seaton, F. M., Kenny, J. G., Eccles, R. M., et al. (2019). Divergent national-scale trends of microbial and animal biodiversity revealed across diverse temperate soil ecosystems. Nat. Commun. 10, 1107–1111. doi: 10.1038/s41467-019-09031-1
Gossner, M. M., Lewinsohn, T. M., Kahl, T., Grassein, F., Boch, S., Prati, D., et al. (2016). Land-use intensification causes multitrophic homogenization of grassland communities. Nature 540, 266–269. doi: 10.1038/nature20575
Graham, E. B., Knelman, J. E., Schindlbacher, A., Siciliano, S., Breulmann, M., Yannarell, A., et al. (2016). Microbes as engines of ecosystem function: when does community structure enhance predictions of ecosystem processes? Front. Microbiol. 7:214. doi: 10.3389/fmicb.2016.00214
Griffiths, R. I., Thomson, B. C., James, P., Bell, T., Bailey, M., and Whiteley, A. S. (2011). The bacterial biogeography of British soils. Environ. Microbiol. 13, 1642–1654. doi: 10.1111/j.1462-2920.2011.02480.x
Hage-Ahmed, K., Rosner, K., and Steinkellner, S. (2019). Arbuscular mycorrhizal fungi and their response to pesticides. Pest Manag. Sci. 75, 583–590. doi: 10.1002/ps.5220
Haichar, F. E. Z., Santaella, C., Heulin, T., and Achouak, W. (2014). Root exudates mediated interactions belowground. Soil Biol. Biochem. 77, 69–80. doi: 10.1016/j.soilbio.2014.06.017
Hannula, S. E., Di Lonardo, D. P., Christensen, B. T., Crotty, F. V., Elsen, A., van Erp, P. J., et al. (2021). Inconsistent effects of agricultural practices on soil fungal communities across 12 European long-term experiments. Eur. J. Soil Sci. 72, 1902–1923. doi: 10.1111/ejss.13090
Hartmann, M., Frey, B., Mayer, J., Mäder, P., and Widmer, F. (2015). Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 9, 1177–1194. doi: 10.1038/ismej.2014.210
Hartmann, M., Niklaus, P. A., Zimmermann, S., Schmutz, S., Kremer, J., Abarenkov, K., et al. (2014). Resistance and resilience of the forest soil microbiome to logging-associated compaction. ISME J. 8, 226–244. doi: 10.1038/ismej.2013.141
Hendgen, M., Hoppe, B., Döring, J., Friedel, M., Kauer, R., Frisch, M., et al. (2018). Effects of different management regimes on microbial biodiversity in vineyard soils. Sci. Rep. 8, 9393–9314. doi: 10.1038/s41598-018-27743-0
Herren, C. M., and McMahon, K. D. (2018). Keystone taxa predict compositional change in microbial communities. Environ. Microbiol. 20, 2207–2217. doi: 10.1111/1462-2920.14257
Hill, M. O. (1973). Diversity and evenness: a unifying notation and its consequences. Ecology 54, 427–432. doi: 10.2307/1934352
Holtkamp, R., van der Wal, A., Kardol, P., van der Putten, W. H., de Ruiter, P. C., and Dekker, S. C. (2011). Modelling C and N mineralisation in soil food webs during secondary succession on ex-arable land. Soil Biol. Biochem. 43, 251–260. doi: 10.1016/j.soilbio.2010.10.004
Jurburg, S. D., Shek, K. L., and McGuire, K. (2020). Soil microbial composition varies in response to coffee agroecosystem management. FEMS Microbiol. Ecol. 96, 1–11. doi: 10.1093/femsec/fiaa164
Kardol, P., Martijn Bezemer, T., and Van Der Putten, W. H. (2006). Temporal variation in plant-soil feedback controls succession. Ecol. Lett. 9, 1080–1088. doi: 10.1111/j.1461-0248.2006.00953.x
Karimi, B., Maron, P. A., Chemidlin-Prevost Boure, N., Bernard, N., Gilbert, D., and Ranjard, L. (2017). Microbial diversity and ecological networks as indicators of environmental quality. Environ. Chem. Lett. 15, 265–281. doi: 10.1007/s10311-017-0614-6
Katayama, N., Bouam, I., Koshida, C., and Baba, Y. G. (2019). Biodiversity and yield under different land-use types in orchard/vineyard landscapes: a meta-analysis. Biol. Conserv. 229, 125–133. doi: 10.1016/j.biocon.2018.11.020
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Lajtha, K., Driscoll, C. T., Jarrell, W. M., and Elliot, E. T. (1999). “Soil phosphorus. standard methods in long term ecological research element analysis,” in Standard Soil Methods Long-Term Ecology Research. New York: (Oxford University Press), 115–142.
Lammel, D. R., Nüsslein, K., Cerri, C. E. P., Veresoglou, S. D., and Rillig, M. C. (2021). Soil biota shift with land use change from pristine rainforest and Savannah (Cerrado) to agriculture in southern Amazonia. Mol. Ecol. 30, 4899–4912. doi: 10.1111/mec.16090
Layeghifard, M., Hwang, D. M., and Guttman, D. S. (2017). Disentangling interactions in the microbiome: a network perspective. Trends Microbiol. 25, 217–228. doi: 10.1016/j.tim.2016.11.008
Lekberg, Y., Arnillas, C. A., Borer, E. T., Bullington, L. S., Fierer, N., Kennedy, P. G., et al. (2021). Nitrogen and phosphorus fertilization consistently favor pathogenic over mutualistic fungi in grassland soils. Nat. Commun. 12, 3484–3488. doi: 10.1038/s41467-021-23605-y
Longepierre, M., Widmer, F., Keller, T., Weisskopf, P., Colombi, T., Six, J., et al. (2021). Limited resilience of the soil microbiome to mechanical compaction within four growing seasons of agricultural management. ISME Commun. 1, 1–13. doi: 10.1038/s43705-021-00046-8
Lori, M., Symnaczik, S., Mäder, P., De Deyn, G., and Gattinger, A. (2017). Organic farming enhances soil microbial abundance and activity—a meta-analysis and meta-regression. PLoS One 12:e0180442. doi: 10.1371/journal.pone.0180442
Lüneberg, K., Schneider, D., Brinkmann, N., Siebe, C., and Daniel, R. (2019). Land use change and water quality use for irrigation alters drylands soil fungal community in the Mezquital Valley, Mexico. Front. Microbiol. 10:1220. doi: 10.3389/fmicb.2019.01220
Lupatini, M., Korthals, G. W., de Hollander, M., Janssens, T. K. S., and Kuramae, E. E. (2017). Soil microbiome is more heterogeneous in organic than in conventional farming system. Front. Microbiol. 7:2064. doi: 10.3389/fmicb.2016.02064
Lupatini, M., Korthals, G. W., Roesch, L. F. W., and Kuramae, E. E. (2018). Long-term farming systems modulate multi-trophic responses. Sci. Total Environ. 646, 480–490. doi: 10.1016/j.scitotenv.2018.07.323
Mandl, K., Cantelmo, C., Gruber, E., Faber, F., Friedrich, B., and Zaller, J. G. (2018). Effects of glyphosate-, glufosinate- and flazasulfuron-based herbicides on soil microorganisms in a vineyard. Bull. Environ. Contam. Toxicol. 101, 562–569. doi: 10.1007/s00128-018-2438-x
Martínez-García, L. B., Korthals, G., Brussaard, L., Jørgensen, H. B., and De Deyn, G. B. (2018). Organic management and cover crop species steer soil microbial community structure and functionality along with soil organic matter properties. Agric. Ecosyst. Environ. 263, 7–17. doi: 10.1016/j.agee.2018.04.018
Meng, J., Li, L., Liu, H., Li, Y., Li, C., Wu, G., et al. (2016). Biodiversity management of organic orchard enhances both ecological and economic profitability. PeerJ 4, e2137–e2122. doi: 10.7717/PEERJ.2137
Meyer, A. H., Wooldridge, J., and Dames, J. F. (2015). Effect of conventional and organic orchard floor management practices on arbuscular mycorrhizal fungi in a ‘Cripp’s pink’/M7 apple orchard soil. Agric. Ecosyst. Environ. 213, 114–120. doi: 10.1016/j.agee.2015.07.026
Millard, P., and Singh, B. K. (2010). Does grassland vegetation drive soil microbial diversity? Nutr. Cycl. Agroecosyst. 88, 147–158. doi: 10.1007/s10705-009-9314-3
Morriën, E., Hannula, S. E., Snoek, L. B., Helmsing, N. R., Zweers, H., De Hollander, M., et al. (2017). Soil networks become more connected and take up more carbon as nature restoration progresses. Nat. Commun. 8:14349. doi: 10.1038/ncomms14349
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., and McGlinn, D. (2020). Vegan: Community Ecology Package. R package version 2.5-7.
Orgiazzi, A., Panagos, P., Yigini, Y., Dunbar, M. B., Gardi, C., Montanarella, L., et al. (2016). A knowledge-based approach to estimating the magnitude and spatial patterns of potential threats to soil biodiversity. Sci. Total Environ. 545-546, 11–20. doi: 10.1016/j.scitotenv.2015.12.092
Paul, E. A. (2016). The nature and dynamics of soil organic matter: plant inputs, microbial transformations, and organic matter stabilization. Soil Biol. Biochem. 98, 109–126. doi: 10.1016/j.soilbio.2016.04.001
Pershina, E., Valkonen, J., Kurki, P., Ivanova, E., Chirak, E., Korvigo, I., et al. (2015). Comparative analysis of prokaryotic communities associated with organic and conventional farming systems. PLoS One 10:e0145072. doi: 10.1371/journal.pone.0145072
Petersen, I. A. B., Meyer, K. M., and Bohannan, B. J. M. (2019). Meta-analysis reveals consistent bacterial responses to land use change across the tropics. Front. Ecol. Evol. 7:391. doi: 10.3389/fevo.2019.00391
Powell, J. R., and Rillig, M. C. (2018). Biodiversity of arbuscular mycorrhizal fungi and ecosystem function. New Phytol. 220, 1059–1075. doi: 10.1111/nph.15119
Radić, T., Likar, M., Hančević, K., Bogdanović, I., and Pasković, I. (2014). Occurrence of root endophytic fungi in organic versus conventional vineyards on the Croatian coast. Agric. Ecosyst. Environ. 192, 115–121. doi: 10.1016/j.agee.2014.04.008
Reganold, J. P., and Wachter, J. M. (2016). Organic agriculture in the twenty-first century. Nat. plants 2:15221. doi: 10.1038/nplants.2015.221
Reilly, K., Cullen, E., Lola-Luz, T., Stone, D., Valverde, J., Gaffney, M., et al. (2013). Effect of organic, conventional and mixed cultivation practices on soil microbial community structure and nematode abundance in a cultivated onion crop. J. Sci. Food Agric. 93, 3700–3709. doi: 10.1002/jsfa.6206
Rousk, J., Bååth, E., Brookes, P. C., Lauber, C. L., Lozupone, C., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Schlatter, B., Trávníček, J., Meier, C., Keller, O., and Willer, H. (eds.) (2021). “Current statistics on organic agriculture worldwide: area, operators and market,” in The World of Organic Agriculture. Statistics and Emerging Trends 2021. (Bonn: Research Institute of Organic Agriculture FiBL, Frick, and IFOAM – Organics International), 32–34.
Schlatter, D. C., Yin, C., Hulbert, S., Burke, I., and Paulitz, T. (2017). Impacts of repeated glyphosate use on wheat-associated bacteria are small and depend on glyphosate use history. Appl. Environ. Microbiol. 83, e01354–e01317. doi: 10.1128/AEM.01354-17
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models. Genome Res. 13:426. doi: 10.1101/gr.1239303.metabolite
Shi, Y., Grogan, P., Sun, H., Xiong, J., Yang, Y., Zhou, J., et al. (2015). Multi-scale variability analysis reveals the importance of spatial distance in shaping Arctic soil microbial functional communities. Soil Biol. Biochem. 86, 126–134. doi: 10.1016/j.soilbio.2015.03.028
Steffan, S. A., and Dharampal, P. S. (2019). Undead food-webs: integrating microbes into the food-chain. Food Webs 18:e00111. doi: 10.1016/j.fooweb.2018.e00111
Szoboszlay, M., Dohrmann, A. B., Poeplau, C., Don, A., and Tebbe, C. C. (2017). Impact of land-use change and soil organic carbon quality on microbial diversity in soils across Europe. FEMS Microbiol. Ecol. 93, 1–12. doi: 10.1093/femsec/fix146
Terrat, S., Horrigue, W., Dequietd, S., Saby, N. P. A., Lelièvre, M., Nowak, V., et al. (2017). Mapping and predictive variations of soil bacterial richness across France. PLoS One 12, e0186766–e0186768. doi: 10.1371/journal.pone.0186766
Tilman, D., Clark, M., Williams, D. R., Kimmel, K., Polasky, S., and Packer, C. (2017). Future threats to biodiversity and pathways to their prevention. Nature 546, 73–81. doi: 10.1038/nature22900
Tuck, S. L., Winqvist, C., Mota, F., Ahnström, J., Turnbull, L. A., and Bengtsson, J. (2014). Land-use intensity and the effects of organic farming on biodiversity: a hierarchical meta-analysis. J. Appl. Ecol. 51, 746–755. doi: 10.1111/1365-2664.12219
Turrini, A., Agnolucci, M., Palla, M., Tomé, E., Tagliavini, M., Scandellari, F., et al. (2017). Species diversity and community composition of native arbuscular mycorrhizal fungi in apple roots are affected by site and orchard management. Appl. Soil Ecol. 116, 42–54. doi: 10.1016/j.apsoil.2017.03.016
Van Der Heijden, M. G. A., Bardgett, R. D., and Van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
van der Heijden, M. G. A., and Wagg, C. (2013). Soil microbial diversity and agro-ecosystem functioning. Plant Soil 363, 1–5. doi: 10.1007/s11104-012-1545-4
Van Geel, M., Verbruggen, E., De Beenhouwer, M., Van Rennes, G., Lievens, B., and Honnay, O. (2017). High soil phosphorus levels overrule the potential benefits of organic farming on arbuscular mycorrhizal diversity in northern vineyards. Agric. Ecosyst. Environ. 248, 144–152. doi: 10.1016/j.agee.2017.07.017
Vogeler, I., Cichota, R., Sivakumaran, S., Deurer, M., and McIvor, I. (2006). Soil assessment of apple orchards under conventional and organic management. Aust. J. Soil Res. 44, 745–752. doi: 10.1071/SR06096
Vukicevich, E., Lowery, T., Bowen, P., Úrbez-Torres, J. R., and Hart, M. (2016). Cover crops to increase soil microbial diversity and mitigate decline in perennial agriculture. A review. Agron. Sustain. Dev. 36:48. doi: 10.1007/s13593-016-0385-7
Wagg, C., Hautier, Y., Pellkofer, S., Banerjee, S., Schmid, B., and van der Heijden, M. G. A. (2020). Diversity and asynchrony in soil microbial communities stabilizes ecosystem functioning. bioRxiv [Preprint]. doi: 10.1101/2020.10.02.324012
Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E. E., and van der Heijden, M. G. A. (2019). Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat. Commun. 10, 4841–4810. doi: 10.1038/s41467-019-12798-y
Walker, L. R., Wardle, D. A., Bardgett, R. D., and Clarkson, B. D. (2010). The use of chronosequences in studies of ecological succession and soil development. J. Ecol. 98, 725–736. doi: 10.1111/j.1365-2745.2010.01664.x
Wang, G., Mayes, M. A., Gu, L., and Schadt, C. W. (2014). Representation of dormant and active microbial dynamics for ecosystem modeling. PLoS One 9:e89252. doi: 10.1371/journal.pone.0089252
Wang, W., Wang, H., Feng, Y., Wang, L., Xiao, X., Xi, Y., et al. (2016). Consistent responses of the microbial community structure to organic farming along the middle and lower reaches of the Yangtze River. Sci. Rep. 6, 1–11. doi: 10.1038/srep35046
Waud, M., Busschaert, P., Ruyters, S., Jacquemyn, H., and Lievens, B. (2014). Impact of primer choice on characterization of orchid mycorrhizal communities using 454 pyrosequencing. Mol. Ecol. Resour. 14, 679–699. doi: 10.1111/1755-0998.12229
Williams, D. R., Clark, M., Buchanan, G. M., Ficetola, G. F., Rondinini, C., and Tilman, D. (2020). Proactive conservation to prevent habitat losses to agricultural expansion. Nat. Sustain. 4, 314–322. doi: 10.1038/s41893-020-00656-5
Wu, X., Yang, J., Ruan, H., Wang, S., Yang, Y., Naeem, I., et al. (2021). The diversity and co-occurrence network of soil bacterial and fungal communities and their implications for a new indicator of grassland degradation. Ecol. Indic. 129:107989. doi: 10.1016/j.ecolind.2021.107989
Xia, Q., Zheng, N., Heitman, J. L., and Shi, W. (2022). Soil pore size distribution shaped not only compositions but also networks of the soil microbial community. Appl. Soil Ecol. 170:104273. doi: 10.1016/j.apsoil.2021.104273
Xue, P., Minasny, B., and McBratney, A. B. (2022). Land-use affects soil microbial co-occurrence networks and their putative functions. Appl. Soil Ecol. 169:104184. doi: 10.1016/j.apsoil.2021.104184
Yashiro, E., Pinto-Figueroa, E., Buri, A., Spangenberg, J. E., Adatte, T., Niculita-Hirzel, H., et al. (2016). Local environmental factors drive divergent grassland soil bacterial communities in the western Swiss Alps. Appl. Environ. Microbiol. 82, 6303–6316. doi: 10.1128/AEM.01170-16
Zheng, W., Zhao, Z., Gong, Q., Zhai, B., and Li, Z. (2018). Effects of cover crop in an apple orchard on microbial community composition, networks, and potential genes involved with degradation of crop residues in soil. Biol. Fertil. Soils 54, 743–759. doi: 10.1007/s00374-018-1298-1
Keywords: apple orchards, soil microbiome, agricultural management, semi-natural reference, organic agriculture, integrated pest management
Citation: Hulsmans E, Peeters G and Honnay O (2022) Soil Microbiomes in Apple Orchards Are Influenced by the Type of Agricultural Management but Never Match the Complexity and Connectivity of a Semi-natural Benchmark. Front. Microbiol. 13:830668. doi: 10.3389/fmicb.2022.830668
Edited by:
Yurong Liu, Huazhong Agricultural University, ChinaReviewed by:
Xiangnan Li, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences (CAS), ChinaDi Liu, China Agricultural University, China
Qingfeng Wang, Shanghai Academy of Agricultural Sciences, China
Copyright © 2022 Hulsmans, Peeters and Honnay. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eva Hulsmans, eva.hulsmans@kuleuven.be