Abstract
Conventional agricultural production systems, typified by large inputs of mineral fertilizers and pesticides, reduce soil biodiversity and may negatively affect ecosystem services such as carbon fixation, nutrient cycling and disease suppressiveness. Organic soil management is thought to contribute to a more diverse and stable soil food web, but data detailing this effect are sparse and fragmented. We set out to map both the resident (rDNA) and the active (rRNA) fractions of bacterial, fungal, protozoan and metazoan communities under various soil management regimes in two distinct soil types with barley as the main crop. Contrasts between resident and active communities explained 22%, 14%, 21% and 25% of the variance within the bacterial, fungal, protozoan, and metazoan communities. As the active fractions of organismal groups define the actual ecological functioning of soils, our findings underline the relevance of characterizing both resident and active pools. All four major organismal groups were affected by soil management (pā<ā0.01), and most taxa showed both an increased presence and an enlarged activity under the organic regime. Hence, a prolonged organic soil management not only impacts the primary decomposers, bacteria and fungi, but also major representatives of the next trophic level, protists and metazoa.
Similar content being viewed by others
Introduction
Intensification of food and feed production is often accompanied by increased nutrient inputs, intense pesticide applications, frequent tillage, and irrigation management1. The negative environmental implications of these practices include eutrophication, increased salinization, soil erosion and biodiversity loss2,3. Our knowledge of the consequences of agricultural intensification on belowground biodiversity is still limited and fragmented. It has been reported that larger-sized soil biota are more negatively affected by high-input agricultural practices than soil microbes4, but is it widely acknowledged that more detailed studies are required to map the effects of soil management on soil microbiota5.
Soils harbour a quarter of the worldās biodiversity and reside among the most complex habitats on earth6,7. Soil biota plays a role in many essential soil functions such as nutrient cycling, carbon and water retention, soil texture formation and the interaction with the plant community8. Most intensive interactions between microbes and plants take place in the rhizosphere, where the plant is able to select and boost a subset of the microbial community by the release of rhizodeposits - a broad range of carbon-containing substances (e.g. root cells, mucilage, volatiles and exudates). The composition of the rhizobiome, the subset of the soil biota present in the rhizosphere, is co-determined by plant identity and age9,10,11,12,13. With the advent of affordable high throughput DNA sequencing techniques, the impact of plants on the identity and density of rhizosphere inhabitants can be mapped. Insight in this interaction could help to design soil management measures promoting a rhizobiome that would optimally support plant growth and improve crop yield14.
As in many other habitats, most soil inhabitants have to cope with unpredictable food availability15. In order to survive periods of food scarcity, various microorganisms can reversibly reduce their metabolic activity over an extended period of time16. Such a condition is referred to as a state of dormancy17. In bulk soil, typically 80% of the cells and 50% of the operational taxonomical units (OTUās) are dormant. This so-called āmicrobial seedbankā18 is alert in the sense that it can detect and respond to environmental stimuli (e.g. organic substrates) that are associated with favourable growing conditions19. Plant roots produce and release a broad spectrum of environmental stimuli and, as such, the rhizosphere is a hotspot of microbial activity20,21.
Given the typically high percentage of dormant microbiota in soil, it is relevant to discriminate between the resident and the active microbial community when considering soil ecological processes. The resident community refers here to all organisms present in a certain spatial unit of soil, whereas the active community comprises the fraction of the resident community that is metabolically active. Ribosomal (r)RNA is considered a representation of the active microbial community, while rDNA characterises the total microbial community22,23. Hence, combined profiling of community rDNA and rRNA will provide insight in both aspects of local microbial assemblies. More specifically, such a characterisation will provide information about microbial fractions, whose activity is positively or negatively affected by any kind of external influence. Although a number of soil ecological studies considered both the active and the resident microbial community24,25,26, large scale mapping of shifts in the active soil microbiome has been hampered so far by the low throughput nature and the costs of currently available kits for RNA extraction from soil. A combination of elements from various published protocols allowed us to develop a fast and affordable method for nucleic acid extraction from soil.
The aim of this study was to investigate the long-term impact of soil management regimes (including conventional and organic soil management) on resident (rDNA) and active (rRNA) microbial communities. To this end, we collected bulk and rhizosphere samples in two different growth stages of summer barley (Hordeum vulgare) grown on two distinct soil types - peaty and sandy soil - under the various types of soil management. Four major organismal groups were assessed: bacteria and fungi - representing the primary decomposers - and protists and nematodes - two major grazers on the bacterial and fungal communities. Specific variable ribosomal DNA regions were selected for the characterization of each of the four organismal groups. We hypothesise that (i) prolonged exposure to distinct soil management regimes will impact both the resident and active fractions of the primary decomposers (bacterial and fungal community), (ii) shifts in the primary decomposer community due to soil management will be translated into associated changes in the active fractions of major representatives of the next trophic level (protists and metazoa), (iii) exposure to rhizosphere will have a stronger stimulating effect on the active fractions of the primary decomposers (i.e. bacteria and fungi) than on protists and metazoa.
Results
Essentials on the characterisation of soil biota
Nucleic acids (total DNA and RNA) from 104 bulk soil and rhizosphere samples were extracted using a novel, lab-made protocol (TableĀ 1). Contrary to bacteria, fungi and protists, the 2āg soil subsamples will not provide a proper representation of the metazoan community. We included metazoa in this study because this fraction of the soil microfauna was co-extracted with the microbial community, and therefore probably in close physical contact with this community. MiSeq sequencing of organismal group-specific 16S (bacteria) or 18S (fungi, protists, and metazoa) ribosomal DNA and cDNA fragments resulted in ā 31 million reads (15.5 million forward and 15.5 million reverse), and on average ā 75,000 reads per sample. After filtering, a total of 8,297,203 sequences were retained comprehending 724 samples for all taxa together. Comprehensive sampling of the microbial community was obtained for all treatments, with average sequence coverage of 63%, 70%, 96% and 97% for respectively bacteria, protozoa, fungi and metazoan as determined by Goodās coverage estimator (Supplementary TableĀ S3).
Contrasts between resident versus active fractions of four major soil food web constituents
PERMANOVA identified Nuclei Acid (i.e. rRNA and rDNA) as the main factor responsible for the differences between samples for all four organismal groups (TableĀ 2). The factor explained 14 to 25% of the overall variance (Pā<ā0.01). The second important explanatory factor of the observed variation was Location (6 to 13%), which is assumed to be attributable mainly to differences in soil type. āVredepeelā is characterized by sandy soils, whereas peaty soil typifies āValthermondā. For the primary decomposers, a higher percentage of the overall variation was explained by soil type (11% and 13% for bacteria and fungi, respectively) when compared to the representatives of the next trophic level (8 and 6% for protists and metazoans, respectively) (for all groups Pā<ā0.01).
UniFrac, a method that uses phylogenetic distances as a measure for the comparison of microbial communities27, was used to verify the impact of individual variables. Both weighted and unweighted UniFrac demonstrated that all variables that were shown to have a significant effect on the bacterial community. This was in essence also true for fungi, although Sample Type was not significant in case of unweighted UniFrac (Pā=ā0.069). For Protozoa, only the effect of treatment was not significant while using unweighted UniFrac (Pā=ā0.097). In case of metazoa, the only non-significant variable was Time Point for the unweighted UniFrac (Pā=ā0,084) (see Supplementary TableĀ 4).
In addition to the PERMANOVA, we generated a principal coordinate analysis (PCoA) ordination of a Bray-Curtis dissimilarity matrix (Fig.Ā 1). There is a clear differentiation visible between active (rRNA) and resident communities (rDNA). The resident communities cluster (blue and light blue) and are separated from the active microbial communities (red and ochre) for all four organismal groups. This separation between clusters was most obvious for Bacteria (Fig.Ā 1A) and Protists (Fig.Ā 1C). For Fungi (Fig.Ā 1B) and Metazoans (Fig.Ā 1D) this separation between active and resident communities was visible as well, but less pronounced.
Principal coordinate analysis (PCoA) ordination of a Bray-Curtis dissimilarity matrix. Plots illustrate distances between communities (104 soil samples; for each sample both the resident (rDNA) and the active (rRNA) community were characterized) for each organismal group. (A) Bacteria; (B) Fungi; (C) Protozoa, and (D) Metazoa. Colours were used to distinguish between rRNA-bulk, rRNA-rhizosphere (ochre), and rDNA-bulk (dark blue), rDNA-rhizosphere (light blue). Locations are indicated by an ochre circle (Vredepeel, sandy soil) or a green circle (Valthermond, peaty soil).
Visualisation of the significant location effects that were demonstrated by the PERMANOVA for all four organismal groups resulted in two distinct patterns. Whereas a complete separation of clusters was observed for bacteria and protists (encircled in ochre and green in Fig.Ā 1A,C), a less clear separation was seen for fungal and metazoan communities (Fig.Ā 1B,D). It appeared to be difficult to visualize the impact of all variables for the four organismal groups with a single set of graphic settings. In Supplementary Fig.Ā S3 alternative settings were used to picture the location effect for fungi and metazoans.
Exposure to rhizosphere conditions resulted in a difference between both the resident and the active microbial community. For bacteria, the difference in community composition between bulk and rhizosphere soil is reflected at both rDNA and rRNA level (blue and light blue, and red and ochre, respectively) (Fig.Ā 1A). For protists, differences between bulk and rhizosphere soil were most pronounced at rDNA level (Fig.Ā 1C). At rRNA level, protist communities only showed differentiation in peaty soils (location Valthermond). It is noted that the PERMANOVA pinpointed significant rhizosphere effects for all four organismal groups (TableĀ 2).
Microbial taxa contributing to contrasts between resident and active communities
To pinpoint taxon-specific differences between resident and active communities and different soil management, we further analysed the samples of location āVredepeelā. Based on LEfSe with LDA thresholds for discriminative features set at ā„2 or ā¤ā2, a total of 9 bacterial, 8 fungal, 11 protozoan, and 12 metazoan orders contributed significantly to the differences between the resident (rDNA-based) and active (rRNA-based) communities (Fig.Ā 2)
LEfSe analysis of bacterial, fungal, protists and metazoan OTUs identifying taxa for which a major part of the population was active in rhizosphere (light green) or in bulk soil (light brown) (LDA score >2), or for which a major part of the population is dormant in rhizosphere (green) or in bulk soil (brown) (LDA scores <ā2).
Four bacterial orders were identified as highly active in both rhizosphere and bulk soil (Fig.Ā 2A). For barley rhizosphere the active Sphingobacteriales were distinct, whereas Bacillales, Acidimicrobiales, Solirubrobacterales and Propionibacteriales were predominately found in bulk soil. It is noted that only Propionibacteriales showed a large contrast between bulk and rhizosphere samples, the other three bacterial orders had LDA scores just below 2 in rhizosphere soils (data not shown).
Analysis of the fungal community revealed four orders that were highly active in both bulk and rhizosphere soil. The order Glomerellales was present (DNA) but barely active in both soil compartments (LDA score below -2). The bulk soil was typified by active members of the orders Filobasidiales and Endogonales. Notably, the contrast between bulk soil and rhizosphere was subtle for Filobasidiales, and more substantial for members of the Endogonales (data not shown).
Protists, being a predominantly bacterivorous group in soil, were included as major representatives of the next trophic level. Colpodida, Philasterida and Haptoria were identified as protist orders with an enhanced metabolic activity in the barley rhizosphere.
Our analyses revealed active rotifers, mites, nematodes, and insects in the rhizosphere compartment. A pronounced difference (LDA distance of 8) between the activity of the bacterivorous nematode order Rhabditida in bulk and rhizosphere soil. Rhabditida are known as extreme opportunists28. This ecological characteristic is reflected in Fig.Ā 2.
Effects of soil management regimes on community structures
Based on our analysis, compost treatment at location āValthermondā had no significant impact on any of the analysed organismal groups (see Supplementary TableĀ S5). In contrast, at location āVredepeelā the three soil management regimes (ConMin, ConSlu and Org) showed to have a significant effect on the microbial community structure. Both the PERMANOVA and a PCoA (TableĀ 3, Fig.Ā 3) indicated a distinct microbial community structure for Org fields as compared to communities found under ConMin and ConSlu management. This soil management effect was most evident for the active communities (rRNA) of Bacteria, Fungi and Protozoa. Results were verified with weighted and unweighted UniFrac. For both UniFrac variants the effect of soil management (āTreatmentā) was significant for all four organismal groups (Supplementary TableĀ S6). A similar analysis on the resident community (rDNA) revealed a comparable, though less pronounced pattern (Supplementary Fig.Ā S4). For Metazoans, no clear soil treatment effect was observed in both the active and the resident communities.
Principal coordinate analysis (PCoA) ordination of a Bray-Curtis dissimilarity matrix. Plots illustrating distances between the active fractions of communities at location Vredepeel (sandy soil) (nā=ā72) for (A) Bacteria, (B) Fungi, (C) Protozoa and (D) Metazoa. Colors were used to indicate soil management regimes: ConMin (purple), ConSlu (orange), and Organic (green).
Bacteria
In general, prolonged organic soil management on sandy soil boosted the abundance of almost all bacterial orders rDNA level (Supplementary Fig.Ā S5.1b). Out of the 38 bacterial taxa that significantly differed in abundance between ConMin and Org fields, 36 taxa were more abundant in the organic treatment. When considering the active fraction of the bacterial community, 47 taxa were significantly affected by soil management, of which 31 taxa were found to be more active in Org. Among the soil management-affected taxa, 16 showed a higher activity in ConMin fields.
As compared to two conventional soil treatments ConMin and ConSlu, prolonged organic soil management has boosted the activity (rRNA) of a range of bacterial orders. Desulfuromonadales, Clostridiales, Erysipelotrichales, Rhodocyclales, and Rhodobacterales had the highest LDA scores (LEfSe, LDA score >2, Fig.Ā 4), and their increased activity was confirmed by ANOVA (grey arrows in Supplementary Fig.Ā S5.1). The orders Bacillales, Deinococcales, Micrococcales, Acidobacteriales, Kineosporiales, and Streptomycetales were identified as indicators for conventional soil treatments. ANOVA did not confirm the status of the order Kineosporiales as indicator taxon for the ConMin treatment (red arrow 4 in Supplementary Fig.Ā S5.1). It is noted that no further attention was paid to bacterial taxa without a formal systematic name.
Discriminant active bacterial, fungal, and protozoan taxa indicated by LEfSe analysis resulting from distinct soil management types at location Vredepeel: ConMin (red), ConSlu (blue) and Org (grey). For each treatment and organismal group, six taxa with the highest LDA scores are delineated.
Fungi
Out of the 30 fungal taxa that significantly differed in OTU abundance (rDNA), 22 taxa were more abundant in organic fields compared to conventional managed fields (Supplementary Fig.Ā S5.2). When concentrating on the fungal orders of which the activity (rRNA) was affected by soil management type, 19 taxa were more active in organic soils and 10 were promoted under conventional soil management regimes (ConMin and ConSlu) (Supplementary Fig.Ā S5.2).
LEfSe analysis of fungal rRNA sequences revealed the orders Glomeromycetes (unclassified class), Cantharellales, Saccharomycetales, Trichosporonales, Agaricales, and Onygenales as indicative for the organic regime (Fig.Ā 4). In most cases, ConSlu-promoted taxa were also activated under the ConMin regime. The order Microascales was identified as a specific indicator for ConSlu-treated fields. ANOVA confirmed the indicator status of this order for the ConSlu treatment (blue arrow 1 in Supplementary Fig.Ā S5.2). The orders Paraglomerales, Eurotiales, Neocallimastigales, and Chaetothyriales were identified as indicators for conventional farming system using mineral fertilizer only.
Protozoa
The total abundance (rDNA) of 28 protist taxa was upregulated in organic treatment, while 13 taxa were promoted under conventional soil management (Supplementary Fig.Ā S5.3). The impact of organic soil management at rRNA level showed an opposite trend since only 12 out of 41 taxa showed a higher OTU abundance in the fields under organic management.
Thraustochytrida and Thecamoebida were identified as protist indicators for organic soil management. Unclassified members of the classes Labyrinthulea and Heterolobosea were characterized by higher densities and higher activity in Org fields (Fig.Ā 4., Supplementary Fig.Ā S5.3). Prasiolales, Tribonematales, Cryptofilida, Phytiales, Dermamoebida and Bicoecales were identified as indicator taxa for conventional soil management (ConMin and ConSlu) (Fig.Ā 4). Notably, ANOVA did not confirm the indicator status of Tribonematales and Bicoecales (Supplementary Fig.Ā S5.3).
Metazoa
The impact of the three soil managements regimes had little impact on the metazoan communities. Only a few orders were characterized as indicators by LEfSe analysis. Mononchida, an order of predatory nematodes, and the mollusc superfamily Nuculoidea were more actively abundant under organic farming. The nematode orders Dorylaimida and Areaolaimida were more active under the ConSlu regime, while Tylenchida and Monhysterida were indicative for ConMin. ANOVA gave non-corresponding results for a number of the aforementioned orders.
Discussion
Characterisation of both resident and active fractions of bacterial and fungal communities as well as their primary consumers, protists and metazoans, generated a holistic view on long-term effects of various soil management systems including organic farming. Analysis of 104 samples from fields with barley as main crop underlines the relevance of distinguishing active (rRNA) and resident (rDNA) communities. For all four organismal groups nucleic acid type was the most important explanatory variable. The long-term impact of soil management was significant for all four organismal groups and dozens of taxa were identified that showed an increased presence and an enlarged activity under the organic soil management regime.
A few aspects of the use of rRNA and rDNA as markers for active and resident soil biota require further attention. Ribosomal RNA is a highly abundant transcript, and rDNA is a multi-copy gene for which the number of copies varies enormously between organismal groups. Among bacteria, the number of rRNA operons is moderately diverse; typically, bacterial genes harbour 1 to 15 copies29. For protists, the number of rDNA copies is substantially higher. Focussing on a range of diatoms and dinoflagellates, individual species were shown to harbour rDNA copy number in the range of 100 to 10,00030. For fungal taxa, a recent study estimated the number of rDNA copies to range between 14 and 1,400. This variation in fungal rDNA copy number could not be linked to trophic preferences or other easily observable ecological characteristics31. Based on copy number variation, it is clear that rRNA and rDNA data from microbial communities should not be used for quantitative comparisons between organismal groups, neither should it be used for comparison of abundances or activities between high taxonomic level taxa within an organismal group. In this study, we infer rRNA to represent the active community. Nevertheless, dormant soil inhabitants may have ribosomes with functional rRNAs32 to allow organisms to resume activity as soon as condition are favourable again. However, it is reasonable to assume that the transition from a dormant to a physiologically active state will be accompanied by a substantial increase in rRNA level. Although it is appreciated that the use of rRNA/rDNA data for the characterisation of active and resident microbial communities in soil has a number of inherent constraints, within-taxon comparison of rDNA-based sequence data under various environmental conditions is likely to reveal robust and valuable information on the impact of these conditions.
Community characterisation by either of the two types of nucleic acids revealed that dormancy was a phenomenon relevant for all four organismal groups under investigation.
Bacteria
Several bacterial orders were shown to be active in both bulk soil and in the rhizosphere. Sphingobacteriales, as an exception, was indicated as specifically active in the rhizosphere (Fig.Ā 2). Pfeiffer et al.33 found a similar enrichment of this order in the rhizosphere of maize. A study by Haichar et al.34 showed Sphingobacteriales accumulation on roots of wheat, rape and barrel clover. This apparent general rhizosphere accruement is in line with the observed upregulated activity near the roots of barley. Propionibacteriales were identified as being specifically active in bulk soil. Propionibacteriales are known to contribute to both primary and secondary fermentations35. The observed activity of Propionibacteriales in bulk soil could therefore relate to the distinct fertilisation regimes.
Fungi
Both in bulk soil and in the rhizosphere representatives of the orders of Sordariales, Hypocreales, Helotiales were highly active (all belonging to the Ascomycota). This in accordance with a survey on four conventional arable fields in Austria that revealed the dominant presence of the same three fungal orders36. However, the high abundance of these Ascomycete orders is not specific for arable soils, Tedersoo et al.37 found similar orders associated with tree roots. Both these studies indicate high abundance (based on DNA data) of the three orders and our data demonstrates that these fungal orders are also highly active. Glomerellales (Ascomycota) showed a negative LDA score for bulk as well as rhizosphere, indicating dormancy. Drought has been shown to specifically decrease protein abundances of Glomerellales38. As sampling took place during a relatively dry period, this could explain the observed dormancy. Unfortunately, little is known about the ecology of Glomerellales.
Protozoa
Barley rhizosphere was characterized by active representatives of the orders of Haptoria, Colpodida and Philasterida. Haptoria are ubiquitous free-living predatory ciliates in soils39, whereas Colpodida are predominantly grazers of bacteria40. The high bacterial density in the rhizosphere could explain the accumulation of active Colpodida. Philasterida are known to occur in terrestrial habitats, but additional information about their ecology is scarce. Bulk soil was typified by active Euamoebida, a supergroup that has been indicated active in grassland and forest mineral soils41. In the same study, the paraphyletic class Variosea was found to be a dominant active class in bulk soil collected from grassland. Our study confirms the presence of active members of Variosea in bulk soil from barley fields. This class was defined only recently42, and little is known about their ecological role in soil. Physarales is suggested to play an important role in litter breakdown43 and therefore it is no surprise that they are specifically active in bulk soils.
Metazoa
As compared to bacteria, fungi and protists, metazoa showed the strongest compartment effect (R2ā=ā0.041, TableĀ 3 (Metazoa / Sample type)). This might point at a large difference in community composition between bulk soil and rhizosphere keeping in mind that the subsamples analysed in this study were too small to be representative for each of the compartments. Nevertheless, this result is in line with the highly density and activity of soil micro fauna in the immediate vicinity of plant roots that has been well reported for a range of systems44,45,46. Our analysis indicated also a number of orders as dormant (LDA <ā2). This category included the orders Haplotaxida (oligochaetes), Dorylaimida (nematodes) and Parachaela (tardigrades). For Haplotaxida and Dorylaimida the observed signals may have originated from unhatched eggs. A number of plant parasites reside within the order Dorylaimida, and these may remain unhatched until signals from a suitable host plant are perceived. Tardigrades including the Parachaela member Hypsibius dujardini are known for their ability to survive relatively dry conditions such as in the upper soil layers at the time of sample collection in a dormant state47. In case of bulk soil, the dormancy of the nematode order Rhabditida was most prominent. This is not unexpected as members of the Rhabditidae are highly opportunistic bacterivores that enter the dormant Dauer stage under unfavourable conditions (e.g. drought, food scarcity)48,49. It is noted that due the relatively large size of soil metazoa, environmental DNA could have contributed to relatively low RNA to DNA ratio which could erroneously be interpreted as a signal for dormancy50.
The effect of organic soil management on bacterial communities is relatively well documented51,52,53,54. To the best of our knowledge there are no other studies investigating the effects of distinct soil management regimes on four organismal groups simultaneously.
Bacteria
At DNA level, ANOVA revealed that 44 out of the 48 bacterial taxa were more abundant under organic management (Supplementary Fig.Ā 5.1). This observation is in line with earlier research on the effects of organic soil management55. However, at activity level (as revealed by rRNA-based analysis) āonlyā 1/3 of the 50 bacterial taxa were most active under organic soil management conditions (Supplementary Fig.Ā 5.1). As microbial activity matters in terms of soil food web functioning, these results underline the relevance of taking ā next to abundance data ā activity data into account.
The following bacterial orders contributed most to the difference between organic versus conventional soil management regimes (in terms of OTU abundance and LEfSe score): Desulfuromonadales (Ī“-proteobacteria), Clostridiales (Firmicutes), Rhodocyclales (Ī²-proteobacteria), Rhodobacterales (Ī±-proteobacteria). Desulfuromonadales, an order typifying Org fields, harbour a range of sulphate and sulphur reducing bacteria56. Their enhanced activity could relate to the slightly higher S content of the Org fields (Org:247āmgāS/kg soil, ConMin and ConSlu: 193 and 214āmgāS/kg)57,58. Clostridiales are metabolically diverse but increased abundance of members of this order has been observed upon the addition of organic matter with a high recalcitrant C content59. The Org fields studied here received crop residues as green manure, and hard-to-degrade plant remains could have promoted the Clostridiales. Rhodocyclales and Rhodobacterales may act as denitrifiers under low oxygen conditions60,61, but the background of their activation under the Org regime remains unclear.
Conventional fields were characterized by highly active Actinobacteria (Kineosporiales, Streptomycetales and Micrococcales). This was also observed in previous research on the impact of conventional and organic cropping systems on the microbial community. In a survey over 3 years, Orr et al.62 detected a similar increase in Actinobacteria, but they also showed a strong sample year effect. For the interpretation of our data this should be taken into consideration.
Fungi
In total 70% of the significantly soil management-affected fungal orders were more abundant under the organic regime. In general, this higher abundance was accompanied by higher activity. The majority of the fungal OTUs were assigned to the Ascomycota. Among the Ascomycota, there are numerous decomposers of organic substrates (such as leaf litter, wood, and manure) and more studies reported them as the major fungal phylum present in agro-ecosystems63,64. Two Basidiomycetes (Agaricales, Cantharellales) were found to be abundant in organic fields. Both orders harbour numerous wood and litter decomposer taxa65, substrates that are added to these fields in relatively large quantities. Onygenales were rarely found in the conventional fields but highly active in the organic soils. This order is associated with animal dung66,67. So, active members of the Onygenales are probably the result of the application of cattle manure in the Org fields.
The strongest fungal indicator for organic farming was an āunclassified class of Glomeromycetesā. Glomeromycetes are known to form arbuscular mycorrhizas and colonize the roots of vascular land plants including barley68. In a recent study on the same experimental farm (Vredepeel), AM fungi were also found to be more abundant under organic soil management57,58. AM fungi can stimulate the decomposition of recalcitrant organic matter, and makes nitrogen bioavailable69. Hence, the distinct type of manure used under organic management might explain the specific stimulation of Glomeromycetes.
Paraglomerales, an order of the Glomeromycetes was predominantly found in conventional systems. This finding was corroborated by Dai et al.70 who found Paraglomus to be positively associated with the conventional production of wheat. The relation between Paraglomerales and fertilization system would require further investigation.
Protozoa
In contrast to the primary decomposers, organic soil management decreased the activity of many protozoa. In parallel, we observed an increase in total abundance (rDNA) under the organic regime for a majority of the soil management-affected taxa. Increased densities of protozoa as a result of organic amendments have been reported before71. Under controlled greenhouse conditions, application of organic fertilizers increased bacterivorous and omnivorous protists, and strongly reduced the relative abundance of plant pathogenic protists72. We arenāt aware of other studies on the impact of organic amendment to soil protist activity.
Metazoa
Organic soil management stimulated the activity an order of predatory nematodes Mononchida, and members of the mollusc order Nuculoidea. Predatory Mononchida feed on other nematodes, but this does not hold for all life stages. Larval stages are too small to capture other nematodes, and they feed on bacteria73. The strongly enriched bacterial community under the organic regimes may have promoted the activity of the Mononchida. The impact of soil management on metazoa will not be discussed further, as the numbers of individuals present in the 2āg soil samples are relatively low (with some nematode taxa as an exception). Hence, sampling effects could easily obscure soil management effects.
Our results demonstrate that prolonged (>15 consecutive years) exposure to distinct soil management regimes causes shifts in the primary decomposer assemblies (bacteria and fungi) as well as changes in primary consumer communities (protists and metazoa). It is concluded that organic management practices results soil microbial communities that are demonstrable distinct from the communities under the conventional regimes. However, our fragmentary knowledge about the ecology and food preference of soil microbiota limited our ability to link most of the observed shifts to desired soil-bound ecosystem services.
Methods
Study sites
Samples were collected from barley fields at two locations in The Netherlands: (1) WUR experimental farm āVredepeelā is located in the southeast of the Netherlands (51Ā°32āN and 5Ā°51E) and is characterized by sandy soil (93,3% sand, 4.5% silt, 2.2% clay) and an organic matter (OM) content of 3ā5%. Three different soil management strategies were applied from 2001 onwards: ConMin, ConSlu and Org. ConMin fields solely received mineral fertilizer and processed organic fertilizer without organic matter (liquid mineral concentrates), and ConSlu fields were supplemented with mineral fertilizer and pig/cow slurry. In case of organic soil management, farmyard manure and cow slurry were applied, and no pesticides were used. For further details of the set up and layout of the soil management experiments see additional research reports74,75,76,77. (2) WUR experimental farm in Valthermond is situated in the northeast of the Netherlands (52Ā°50ā²N, 6Ā°55ā²E) and characterized by sandy peat soil (90% sand, 7% silt, 3% clay) and a high OM content (up to 14%). At Valthermond, the effect of the application of compost was investigated (yearly application of 15 tons (green) compost per hectare) since 2013. For this study, we made a comparison between the control and compost plots.
Soil sampling
At both experimental farms, barley (Hordeum vulgare) is one of the main crops in the crop rotation system. Due to a slight latitudinal difference, development of the barley plants in āValthermondā was one week delayed as compared to āVredepeelā. Samples were collected at two time points in spring 2017 (during the vegetative and the generative stage of the crop, see also Supplementary TableĀ S1).
At the āVredepeelā, each of the three fields was divided in 6 subfields of 540ām2 (Supplementary Fig.Ā S1). In each subfield, a bulk soil and a rhizosphere sample were collected. Rhizosphere composite samples were taken by harvesting all barley plants within a vicinity of 20āĆā20ācm. Excessive soil was removed from the root system by shaking the plants. Immediately thereafter, the plants were transferred to a field laboratory, and rhizosphere samples were collected by brushing off soil that adhered to the roots from 10 individual barley plants. For bulk soil, three cores were collected between the barley rows using an auger (Ćø1.5ācm, depth approximately: 15ācm), and thoroughly mixed in pre-labelled bags.
In total 36 composite samples (18 rhizosphere and 18 bulk) were taken at each time point. Rhizosphere soil and bulk soil samples were frozen immediately after sampling in liquid nitrogen and transported on dry ice to the laboratory and stored at ā80āĀ°C.
At the āValthermondā location, samples were taken in the first 2 meters of the subfield. In total 4 subfields of each treatment (Supplementary Fig.Ā S2) were sampled resulting in 16 samples (8 rhizosphere and 8 bulk) at each time point. Barley rhizosphere samples were collected as described above. Hence, a total of 104 soil samples (72 from āVredepeelā, 36 from āValthermondā) were used for further analysis.
DNA/RNA extraction and cDNA synthesis
Both DNA and RNA were simultaneously extracted from soil samples, using a lab-made protocol based on a combination of published protocols. This protocol uses a subsample of 2āg from a thoroughly mixed soil sample as starting point (TableĀ 1). Contrary to bacteria, fungi and protozoa, a subsample of 2āg will not give a proper representation of the metazoan community78. Metazoan DNAs were co-extracted with DNAs from other microbiota, and the resulting data represent microfauna that was presumably living in the close physical vicinity of the observed microbiota. For this reason, metazoa were included in this study.
After precipitation of humic substances from the 2āg subsample by the chemical flocculant NH4 Al (SO4)279, and the removal of proteins by a standard phenol:chloroform:isoamyl alcohol mixture, samples were physically disrupted by bead beating using a Vortex Genie 2 with tubes attached to a horizontal tube holder. After a standard isopropanol precipitation, nuclei acids were re-dissolved in a binding solution. DNA and RNA were further purified using a silica-based RP20 CommaPrep RNA extraction column80 (see TableĀ 1 for technical details). Quality and quantity of the obtained RNA and DNA was measured with a Nanodrop and Qubit. Until further processing, nucleic acid eluates were stored in ā80āĀ°C.
For synthesis of cDNA from extracted RNA the Maxima First Strand cDNA Synthesis Kit for RT-qPCR (Fermentas, Thermo Fisher Scientific Inc., USA) was used according to the manufacturerās instructions. All individual DNA and cDNA samples were diluted to 1āng/ul and 0.1āng/ul respectively, to serve as a template for PCR amplification.
PCR amplification and sequencing
For the characterisation of bacterial, fungal, protist and metazoan communities, variable regions (V) of 16S or 18S ribosomal DNA were amplified. For bacteria, the V4 region was amplified, while for protozoa, fungi and metazoa respectively the V9, V7-V8, V5-V7 regions were targeted (for details see TableĀ 4). To prepare the samples for Illumina sequencing, a two-step PCR procedure was followed. In the first PCR, locus-specific primers extended with an Illumina read area and the appropriate adapter (TableĀ 4) were used to produce primary amplicons. Three Āµl of diluted DNA or cDNA template was used with the following temperature profile: 3āmin 95āĀ°C, followed by 35āĆā(95āĀ°C, 10ās; 55āĀ°C, 20ās; 72āĀ°C, 20ās) and a final extension step of 72āĀ°C of 5āmin. This was done in triplicate for all samples and for each of the four organismal groups. The second PCR step was performed on 40x diluted amplicons of PCR step 1. This PCR 2 was conducted to attach the Illumina index and the Illumina sequencing adaptor (3āmin 95āĀ°C, followed by 10āĆā(95āĀ°C, 10ās; 60āĀ°C, 30ās; 72āĀ°C, 30ās) and a final extension step of 72āĀ°C of 5āmin). Products of PCR 1 and 2 were randomly checked on gel to ensure amplification was successful. Finally, all PCR products were pooled and sent for sequencing (Bioscience, Wageningen Research, Wageningen, The Netherlands) using the Illumina MiSeq Desktop Sequencer (2*300nt paired-end sequencing) according to the standard protocols.
Bioinformatics pipeline
The composition of microbial communities of the soil samples was analysed based on the sequencing data obtained from the Illumina MiSeq platform. Reads were first sorted into the experimental samples according to their index combination. Thereafter, they were sorted into the four organismal groups based on their locus-specific primer sequences (general run statistics can be found in Supplementary TableĀ S2).
Sequences were processed with Hydra pipeline version 1.3.381 implemented in Snakemake82. Forward and reverse reads were paired only for bacteria and fungi while single-end sequences were analysed for protozoa and metazoa. The four taxonomical groups were quality trimmed by BBDUK and then merged via VSEARCH83,84. Resulting unique sequences were then clustered at 97% similarity by using the usearch_global method implemented in VSEARCH and a representative consensus sequence per de novo OTU was determined84. The clustering algorithm also performs chimera filtering to discard likely chimeric OTUs with UCHIME algorithm in de novo mode85 implemented in VSEARCH. Sequences that passed quality filtering were then mapped to a set of representative consensus sequences to generate an OTU abundance table. Representative OTU sequences were assigned to a taxonomic classification via BLAST against the Silva database (version 12.8) for bacteria, fungi and metazoa, and against the Protist Ribosomal Reference database PR286 for protozoa using SINA87. Sequences belonging to chloroplasts, cyanobacteria and mitochondria were discarded from the bacterial dataset, and sequences not belonging to Fungi and Metazoa were removed from the 18S Fungi and Metazoa datasets. The Protozoa dataset was filtered for Streptophyta, Metazoa, fungal and unclassified Opisthokonta sequences. Low-abundance OTUs (those with abundance of <0.005% in the total data set) were discarded prior to analysis88. Samples were transformed using Hellingersā transformation for all downstream analyses.
Statistical analyses
Sampling effort was estimated by Goodās coverage89. For statistical analysis, we explored Ī² diversity patterns by performing principal coordinate analysis (PCoA) with Bray-Curtis dissimilarity using QIIME software90. Permutational multivariate analysis of variance (PERMANOVA) was used to compare the microbial community structure between soil managements taken from different sites and with different plant growth stages for active and resident community for 4 different taxa. This was performed with 999 permutations using the adonis function, based on Bray-Curtis and UniFrac (weighted and unweighted) distances using the āveganā package91 in R. To investigate the indicator taxa involved in the differences between resident and active community, a linear discriminate analysis (LDA) effect size (LEfSe) was conducted in Microbiome Analyst92 to explore the differential microbial populations at the family level for the four different taxa93. A significance level of Ī±āā¤ā0.05 was used for all biomarkers evaluated in this study.
Data Availability
The raw sequences were submitted to the NCBI Sequence Read Archive (SRA) database under study accession number BioProject ID PRJNA543417.
References
Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl. Acad. Sci. 108, 20260ā20264, https://doi.org/10.1073/pnas.1116437108 (2011).
Foley, J. A. et al. Global consequences of land use. Science 309, 570ā574, https://doi.org/10.1126/science.1111772 (2005).
Tilman, D. et al. Forecasting agriculturally driven global environmental change. Science 292, 281ā284, https://doi.org/10.1126/science.1057544 (2001).
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 91, 460ā473, https://doi.org/10.1890/09-0666.1 (2010).
El Mujtar, V., Munoz, N., Mc Cormick, B. P., Pulleman, M. & Tittonell, P. Role and management of soil biodiversity for food security and nutrition; where do we stand? Glob Food Sec. 20, 132ā144, https://doi.org/10.1016/j.gfs.2019.01.007 (2019).
Bardgett, R. The biology of soil: a community and ecosystem approach. (Oxford University Press, 2005).
Roesch, L. F. et al. Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283ā290, https://doi.org/10.1038/ismej.2007.53 (2007).
Henneron, L. et al. Fourteen years of evidence for positive effects of conservation agriculture and organic farming on soil life. Agron Sustain Dev 35, 169ā181, https://doi.org/10.1007/s13593-014-0215-8 (2015).
Chaparro, J. M., Badri, D. V. & Vivanco, J. M. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 8, 790ā803, https://doi.org/10.1038/ismej.2013.196 (2014).
Grayston, S. J., Wang, S. Q., Campbell, C. D. & Edwards, A. C. Selective influence of plant species on microbial diversity in the rhizosphere. Soil Biol Biochem 30, 369ā378, https://doi.org/10.1016/S0038-0717(97)00124-7 (1998).
Inceoglu, O., Abu Al-Soud, W., Salles, J. F., Semenov, A. V. & van Elsas, J. D. Comparative Analysis of Bacterial Communities in a Potato Field as Determined by Pyrosequencing. PLoS One 6, e23321, https://doi.org/10.1371/journal.pone.0023321 (2011).
Kuske, C. R. et al. Comparison of soil bacterial communities in rhizospheres of three plant species and the interspaces in an arid grassland. Appl Environ Microb 68, 1854ā1863, https://doi.org/10.1128/Aem.68.4.1854-1863.2002 (2002).
Schlemper, T. R. et al. Rhizobacterial community structure differences among sorghum cultivars in different growth stages and soils. FEMS Microbiol Ecol 93, 1ā11, https://doi.org/10.1093/femsec/fix096 (2017).
Knief, C. Analysis of plant microbe interactions in the era of next generation sequencing technologies. Front Plant Sci 5, 216, https://doi.org/10.3389/Fpls.2014.00216 (2014).
Morita, R. Y. Bacteria in oligotrophic environments: starvation-survival lifestyle. Vol. 1 (Chapman & Hall New York, 1997).
Blagodatskaya, E. & Kuzyakov, Y. Active microorganisms in soil: Critical review of estimation criteria and approaches. Soil Biol Biochem 67, 192ā211, https://doi.org/10.1016/j.soilbio.2013.08.024 (2013).
Stevenson, L. H. Case for Bacterial Dormancy in Aquatic Systems. Microbial Ecol 4, 127ā133, https://doi.org/10.1007/Bf02014283 (1978).
Lennon, J. T. & Jones, S. E. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol 9, 119ā130, https://doi.org/10.1038/nrmicro2504 (2011).
De Nobili, M., Contin, M., Mondini, C. & Brookes, P. C. Soil microbial biomass is triggered into activity by trace amounts of substrate. Soil Biol Biochem 33, 1163ā1170, https://doi.org/10.1016/S0038-0717(01)00020-7 (2001).
Hinsinger, P., Bengough, A. G., Vetterlein, D. & Young, I. M. Rhizosphere: biophysics, biogeochemistry and ecological relevance. Plant Soil 321, 117ā152, https://doi.org/10.1007/s11104-008-9885-9 (2009).
Reinhold-Hurek, B., Bunger, W., Burbano, C. S., Sabale, M. & Hurek, T. Roots Shaping Their Microbiome: Global Hotspots for Microbial Activity. Annu Rev Phytopathol 53, 403ā424, https://doi.org/10.1146/annurev-phyto-082712-102342 (2015).
De Vrieze, J. et al. Presence does not imply activity: DNA and RNA patterns differ in response to salt perturbation in anaerobic digestion. Biotechnol Biofuels 9, 244 (2016).
Ofek, M., Voronov-Goldman, M., Hadar, Y. & Minz, D. Host signature effect on plant root-associated microbiomes revealed through analyses of resident vs. active communities. Environ Microbiol 16, 2157ā2167, https://doi.org/10.1111/1462-2920.12228 (2014).
Baldrian, P. et al. Active and total microbial communities in forest soil are largely different and highly stratified during decomposition. ISME J. 6, 248ā258, https://doi.org/10.1038/ismej.2011.95 (2012).
Nunes, I. et al. Soil bacteria show different tolerance ranges to an unprecedented disturbance. Biol Fert Soils 54, 189ā202, https://doi.org/10.1007/s00374-017-1255-4 (2018).
Schostag, M. et al. Bacterial and protozoan dynamics upon thawing and freezing of an active layer permafrost soil. ISME J. 13, 1345ā1359 (2019).
Lozupone, C. & Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl Environ Microb 71, 8228ā8235, https://doi.org/10.1128/AEM.71.12.8228-8235.2005 (2005).
Bongers, T. The Maturity Index - an Ecological Measure of Environmental Disturbance Based on Nematode Species Composition. Oecologia 83, 14ā19, https://doi.org/10.1007/Bf00324627 (1990).
Espejo, R. T. & Plaza, N. Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA. Front Microbiol 9, 1232, https://doi.org/10.3389/fmicb.2018.01232 (2018).
Godhe, A. et al. Quantification of Diatom and Dinoflagellate Biomasses in Coastal Marine Seawater Samples by Real-Time PCR. Appl Environ Microb 74, 7174ā7182, https://doi.org/10.1128/Aem.01298-08 (2008).
Lofgren, L. A. et al. Genome-based estimates of fungal rDNA copy number variation across phylogenetic scales and ecological lifestyles. Mol Ecol 28, 721ā730, https://doi.org/10.1111/mec.14995 (2018).
Nanamiya, H. et al. Bacillus subtilis mutants harbouring a single copy of the rRNA operon exhibit severe defects in growth and sporulation. Microbiology 156, 2944ā2952, https://doi.org/10.1099/mic.0.035295-0 (2010).
Peiffer, J. A. et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. 110, 6548ā6553, https://doi.org/10.1073/pnas.1302837110 (2013).
Haichar, F. Z. et al. Plant host habitat and root exudates shape soil bacterial community structure. ISME J 2, 1221ā1230, https://doi.org/10.1038/ismej.2008.80 (2008).
Johnson, J. & Cummins, C. Cell wall composition and deoxyribonucleic acid similarities among the anaerobic coryneforms, classical propionibacteria, and strains of Arachnia propionica. J Bacteriol 109, 1047ā1066 (1972).
Klaubauf, S. et al. Molecular diversity of fungal communities in agricultural soils from Lower Austria. Fungal Divers 44, 65ā75, https://doi.org/10.1007/s13225-010-0053-1 (2010).
Tedersoo, L. et al. Ascomycetes associated with ectomycorrhizas: molecular diversity and ecology with particular reference to the Helotiales. Environ Microbiol 11, 3166ā3178, https://doi.org/10.1111/j.1462-2920.2009.02020.x (2009).
Bastida, F. et al. Differential sensitivity of total and active soil microbial communities to drought and forest management. Global Change Biol 23, 4185ā4203, https://doi.org/10.1111/gcb.13790 (2017).
VÄaÄnĆ½, P. et al. The Chaos Prevails: Molecular Phylogeny of the Haptoria (Ciliophora, Litostomatea). Protist 165, 93ā111, https://doi.org/10.1016/j.protis.2013.11.001 (2014).
VÄaÄnĆ½, P. & Foissner, W. Re-analysis of the 18S rRNA gene phylogeny of the ciliate class Colpodea. Eur J Protistol 67, 89ā105, https://doi.org/10.1016/j.ejop.2018.11.003 (2018).
Geisen, S. et al. Metatranscriptomic census of active protists in soils. ISME J. 9, 2178ā2190, https://doi.org/10.1038/ismej.2015.30 (2015).
Cavalier-Smith, T., Chao, E. E. Y. & Oates, B. Molecular phylogeny of Amoebozoa and the evolutionary significance of the unikont Phalansterium. Eur J Protistol 40, 21ā48, https://doi.org/10.1016/j.ejop.2003.10.001 (2004).
Kamono, A., Matsumoto, J., Kojima, H. & Fukui, M. Characterization of myxomycete communities in soil by reverse transcription polymerase chain reaction (RT-PCR)-based method. Soil Biol Biochem 41, 1324ā1330, https://doi.org/10.1016/j.soilbio.2009.04.001 (2009).
Chen, X., Liu, M., Hu, F., Mao, X. & Li, H. Contributions of soil micro-fauna (protozoa and nematodes) to rhizosphere ecological functions. Acta Ecologica Sinica 27, 3132ā3143, https://doi.org/10.1016/S1872-2032(07)60068-7 (2007).
Griffiths, B. S. Microbial-Feeding Nematodes and Protozoa in Soil - Their Effects on Microbial Activity and Nitrogen Mineralization in Decomposition Hotspots and the Rhizosphere. Plant Soil 164, 25ā33, https://doi.org/10.1007/Bf00010107 (1994).
Irshad, U., Villenave, C., Brauman, A. & Plassard, C. Grazing by nematodes on rhizosphere bacteria enhances nitrate and phosphorus availability to Pinus pinaster seedlings. Soil Biol Biochem 43, 2121ā2126, https://doi.org/10.1016/j.soilbio.2011.06.015 (2011).
Vasanthan, T., Nederveen, J. P. & Stone, J. Quantum-like decreased embryogenesis time with increased cold exposure time. Sci. Rep. 9, 1229, https://doi.org/10.1038/s41598-018-35396-2 (2019).
Treonis, A. M., Wall, D. H. & Virginia, R. A. The use of anhydrobiosis by soil nematodes in the Antarctic Dry Valleys. Funct Ecol 14, 460ā467, https://doi.org/10.1046/j.1365-2435.2000.00442.x (2000).
Yeates, G. W., Bongers, T., Degoede, R. G. M., Freckman, D. W. & Georgieva, S. S. Feeding-Habits in Soil Nematode Families and Genera - an Outline for Soil Ecologists. J Nematol 25, 315ā331 (1993).
Hu, S. K. et al. Protistan diversity and activity inferred from RNA and DNA at a coastal ocean site in the eastern North Pacific. FEMS Microbiol Ecol 92, fiw050, https://doi.org/10.1093/femsec/fiw050 (2016).
Hartmann, M., Frey, B., Mayer, J., Mader, P. & Widmer, F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J 9, 1177ā1194, https://doi.org/10.1038/ismej.2014.210 (2015).
Lori, M., Symnaczik, S., Mader, P., De Deyn, G. & Gattinger, A. Organic farming enhances soil microbial abundance and activity-A meta-analysis and meta-regression. PLoS One 12, https://doi.org/10.1371/journal.pone.0180442 (2017).
Lupatini, M., Korthals, G. W., de Hollander, M., Janssens, T. K. S. & Kuramae, E. E. Soil Microbiome Is More Heterogeneous in Organic Than in Conventional Farming System. Front Microbiol 7, 460ā470, https://doi.org/10.3389/fmicb.2016.02064 (2017).
Suleiman, A. K. A. et al. Recycling organic residues in agriculture impacts soil-borne microbial community structure, function and N2O emissions. Sci Total Environ 631-632, 1089ā1099, https://doi.org/10.1016/j.scitotenv.2018.03.116 (2018).
MƤder, P. et al. Soil fertility and biodiversity in organic farming. Science 296, 1694ā1697, https://doi.org/10.1126/science.1071148 (2002).
Warren, L. A., Kendra, K. E., Brady, A. L. & Slater, G. F. Sulfur Biogeochemistry of an Oil Sands Composite Tailings Deposit. Front Microbiol 6, 1533, https://doi.org/10.3389/fmicb.2015.01533 (2016).
Martinez-Garcia, L. B., Korthals, G., Brussaard, L., Jorgensen, H. B. & De Deyn, G. B. Organic management and cover crop species steer soil microbial community structure and functionality along with soil organic matter properties. Agr Ecosyst Environ 263, 7ā17, https://doi.org/10.1016/j.agee.2018.04.018 (2018).
Lupatini, M., Korthals, G. W., Roesch, L. F. W. & Kuramae, E. E. Long-term farming systems modulate multi-trophic responses. Sci Total Environ 646, 480ā490, https://doi.org/10.1016/j.scitotenv.2018.07.323 (2019).
Goldfarb, K. C. et al. Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Front Microbiol 2, 94, https://doi.org/10.3389/fmicb.2011.00094 (2011).
Saito, T., Ishii, S., Otsuka, S., Nishiyama, M. & Senoo, K. Identification of Novel Betaproteobacteria in a Succinate-Assimilating Population in Denitrifying Rice Paddy Soil by Using Stable Isotope Probing. Microbes Environ 23, 192ā200, https://doi.org/10.1264/jsme2.23.192 (2008).
Yoshida, M., Ishii, S., Otsuka, S. & Senoo, K. Temporal shifts in diversity and quantity of nirS and nirK in a rice paddy field soil. Soil Biol Biochem 41, 2044ā2051, https://doi.org/10.1016/j.soilbio.2009.07.012 (2009).
Orr, C. H., Stewart, C. J., Leifert, C., Cooper, J. M. & Cummings, S. P. Effect of crop management and sample year on abundance of soil bacterial communities in organic and conventional cropping systems. J Appl Microbiol 119, 208ā214, https://doi.org/10.1111/jam.12822 (2015).
Francioli, D. et al. Mineral vs. Organic Amendments: Microbial Community Structure, Activity and Abundance of Agriculturally Relevant Microbes Are Driven by Long-Term Fertilization Strategies. Front Microbiol 7, 1446, https://doi.org/10.3389/fmicb.2016.01446 (2016).
Lienhard, P. et al. Pyrosequencing evidences the impact of cropping on soil bacterial and fungal diversity in Laos tropical grassland. Agron Sustain Dev 34, 525ā533, https://doi.org/10.1007/s13593-013-0162-9 (2014).
Floudas, D. et al. Evolution of novel wood decay mechanisms in Agaricales revealed by the genome sequences of Fistulina hepatica and Cylindrobasidium torrendii. Fungal Genet Biol 76, 78ā92, https://doi.org/10.1016/j.fgb.2015.02.002 (2015).
Doveri, F., Pecchia, S., Vergara, M., Sarrocco, S. & Vannacci, G. A comparative study of Neogymnomyces virgineus, a new keratinolytic species from dung, and its relationships with the Onygenales. Fungal Divers 52, 13ā34, https://doi.org/10.1007/s13225-011-0120-2 (2012).
Sugiyama, M., Summerbell, R. C. & Mikawa, T. Molecular phylogeny of onygenalean fungi based on small subunit (SSU) and large subunit (LSU) ribosomal DNA sequences. Stud Mycol, 5ā23, https://doi.org/10.1007/BF02463962 (2002).
Williams, A., Manoharan, L., Rosenstock, N. P., Olsson, P. A. & Hedlund, K. Long-term agricultural fertilization alters arbuscular mycorrhizal fungal community composition and barley (Hordeum vulgare) mycorrhizal carbon and phosphorus exchange. New Phytol 213, 874ā885, https://doi.org/10.1111/nph.14196 (2017).
Hodge, A. Arbuscular mycorrhizal fungi influence decomposition of, but not plant nutrient capture from, glycine patches in soil. New Phytol 151, 725ā734, https://doi.org/10.1046/j.0028-646x.2001.00200.x (2001).
Dai, M. et al. Negative and positive contributions of arbuscular mycorrhizal fungal taxa to wheat production and nutrient uptake efficiency in organic and conventional systems in the Canadian prairie. Soil Biol Biochem 74, 156ā166, https://doi.org/10.1016/j.soilbio.2014.03.016 (2014).
Treonis, A. M. et al. Effects of organic amendment and tillage on soil microorganisms and microfauna. Appl Soil Ecol 46, 103ā110, https://doi.org/10.1016/j.apsoil.2010.06.017 (2010).
Xiong, W. et al. Soil protist communities form a dynamic hub in the soil microbiome. ISME J. 12, 634ā638, https://doi.org/10.1038/ismej.2017.171 (2018).
Yeates, G. Nematode feeding and activity: the importance of development stages. Biol Fert Soils 3, 143ā146 (1987).
de Haan, J. et al. Biologische teelt op een zuidelijke zandgrond: opbrengst, bemesting, bodemkwaliteit en stikstofverliezen. ((Wageningen University & Research, Praktijkonderzoek AGV), 2018).
de Haan, J. et al. Effect van organische stofbeheer op opbrengst, bodemkwaliteit en stikstofverliezen op een zuidelijke zandgrond. ((Wageningen University & Research, Praktijkonderzoek AGV), 2018).
Quist, C. W. et al. Organic farming practices result in compositional shifts in nematode communities that exceed crop-related changes. Appl Soil Ecol 98, 254ā260, https://doi.org/10.1016/j.apsoil.2015.10.022 (2016).
Schrama, M., de Haan, J. J., Kroonen, M., Verstegen, H. & van der Putten, W. H. Crop yield gap and stability in organic and conventional farming systems. Agr Ecosyst Environ 256, 123ā130, https://doi.org/10.1016/j.agee.2017.12.023 (2018).
Wiesel, L., Daniell, T. J., King, D. & Neilson, R. Determination of the optimal soil sample size to accurately characterise nematode communities in soil. Soil Biol Biochem 80, 89ā91, https://doi.org/10.1016/j.soilbio.2014.09.026 (2015).
Braid, M. D., Daniels, L. M. & Kitts, C. L. Removal of PCR inhibitors from soil DNA by chemical flocculation. J Microbiol Meth 52, 389ā393, https://doi.org/10.1016/S0167-7012(02)00210-5 (2003).
Miller, D. N., Bryant, J. E., Madsen, E. L. & Ghiorse, W. C. Evaluation and optimization of DNA extraction and purification procedures for soil and sediment samples. Appl Environ Microb 65, 4715ā4724 (1999).
de Hollander, M. NIOO-KNAW/hydra: Version 1.3. 3. Zenodo., https://doi.org/10.5281/zenodo.1434147 (2017).
Kƶster, J. & Rahmann, S. Snakemakeāa scalable bioinformatics workflow engine. Bioinformatics 28, 2520ā2522, https://doi.org/10.1093/bioinformatics/bts480 (2012).
Bushnell, B. BBMap: short read aligner and other bioinformatic tools v36.99, https://sourceforge.net/projects/bbmap/ (2018).
Rognes, T., Flouri, T., Nichols, B., Quince, C. & MahƩ, F. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194ā2200, https://doi.org/10.1093/bioinformatics/btr381 (2011).
Guillou, L. et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy. Nucleic Acids Res 41, D597āD604, https://doi.org/10.1093/nar/gks1160 (2013).
Pruesse, E., Peplies, J. & Glockner, F. O. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823ā1829, https://doi.org/10.1093/bioinformatics/bts252 (2012).
Bokulich, N. A. et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods 10, 57ā59, https://doi.org/10.1038/nmeth.2276 (2013).
Good, I. J. The population frequencies of species and the estimation of population parameters. Biometrika 40, 237ā264 (1953).
Caporaso, J. G. et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J 6, 1621ā1624, https://doi.org/10.1038/ismej.2012.8 (2012).
Oksanen, J. et al. Vegan: community ecology package. R package version 2.3ā0, https://CRAN.R-project.org/package=vegan (2015).
Dhariwal, A. et al. MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res 45, W180āW188, https://doi.org/10.1093/nar/gkx295 (2017).
Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol 12, R60, https://doi.org/10.1186/gb-2011-12-6-r60 (2011).
Anderson, M. A lab-made method for extracting DNA from soils. Soil Res 56, 560ā567, https://doi.org/10.1071/Sr17261 (2018).
Turner, S., Pryer, K. M., Miao, V. P. & Palmer, J. D. Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis 1. J Eukaryot Microbiol 46, 327ā338 (1999).
Verbruggen, E. et al. Testing Potential Effects of Maize Expressing the Bacillus thuringiensis Cry1Ab Endotoxin (Bt Maize) on Mycorrhizal Fungal Communities via DNA- and RNA-Based Pyrosequencing and Molecular Fingerprinting. Appl Environ Microb 78, 7384ā7392, https://doi.org/10.1128/Aem.01372-12 (2012).
Vainio, E. J. & Hantula, J. Direct analysis of wood-inhabiting fungi using denaturing gradient gel electrophoresis of amplified ribosomal DNA. Mycol Res 104, 927ā936, https://doi.org/10.1017/S0953756200002471 (2000).
Lane, D. In Nucleic Acid Techniques in Bacterial Systematics. Vol. John Wiley &Sons 115ā175 (Inc, 1991).
Medlin, L., Elwood, H. J., Stickel, S. & Sogin, M. L. The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene 71, 491ā499 (1988).
Capra, E. et al. A new primer set for DNA metabarcoding of soil Metazoa. Eur J Soil Biol 77, 53ā59, https://doi.org/10.1016/j.ejsobi.2016.10.005 (2016).
Acknowledgements
We would like to thank WUR experimental farm 't Kompas in Valthermond and WUR experimental farm Vredepeel in Vredepeel for managing and maintaining long-term field experiments on the effect of various farming systems on a broad range of agronomic parameters and the collection of all primary farm data. This research was supported by NWO Groen grants 870.15.021 and 870.15.022.
Author information
Authors and Affiliations
Contributions
P.H., A.K.A.S., J.J.d.H., E.E.K. and J.H. were responsible for the experimental design. P.H., A.K.A.S., M.H.M.H., E.E.K. and J.H. were involved in sampling the rhizosphere. S.J.J.E. developed the RNA/DNA isolation protocol. P.H. tested the RNA/DNA protocol and performed RNA/DNA isolation of the soil samples. S.J.J.E. developed the PCR two step system prior to sequencing. P.H. tested and performed the two step PCR reactions in order to prepare the sequence library. M.H.M.H. sorted the sequence data. A.K.A.S. analysed the sequence data and performed the bioinformatics and statistical analysis. P.H. and J.H. wrote the manuscript; all others co-commented on the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisherās note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleās Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleās Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Harkes, P., Suleiman, A.K.A., van den Elsen, S.J.J. et al. Conventional and organic soil management as divergent drivers of resident and active fractions of major soil food web constituents. Sci Rep 9, 13521 (2019). https://doi.org/10.1038/s41598-019-49854-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-49854-y
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
