- Research article
- Open access
- Published:
Soil microbial communities following 20 years of fertilization and crop rotation practices in the Czech Republic
Environmental Microbiome volume 17, Article number: 13 (2022)
Abstract
Background
Although fertilization and crop rotation practices are commonly used worldwide in agriculture to maximize crop yields, their long-term effect on the structures of soil microorganisms is still poorly understood. This study investigated the long-term impact of fertilization and crop rotation on soil microbial diversity and the microbial community structure in four different locations with three soil types. Since 1996, manure (MF; 330 kg N/ha), sewage sludge (SF; 330 and SF3x; 990 kg N/ha), and NPK (NPK; 330 kg N/ha) fertilizers were periodically applied to the soils classified as chernozem, luvisol and cambisol, which are among the most abundant or fertile soils used for agricultural purposes in the world. In these soils, potato (Solanum tuberosum L.), winter wheat (Triticum aestivum L.), and spring barley (Hordeum vulgare L.) were rotated every three years.
Results
Soil chemistry, which was significantly associated with location, fertilization, crop rotation, and the interaction of fertilization and location, was the dominant driver of soil microbial communities, both prokaryotic and fungal. A direct effect of long-term crop rotation and fertilization on the structure of their communities was confirmed, although there was no evidence of their influence on microbial diversity. Fungal and bacterial communities responded differently to fertilization treatments; prokaryotic communities were only significantly different from the control soil (CF) in soils treated with MF and SF3x, while fungal communities differed across all treatments. Indicator genera were identified for different treatments. These taxa were either specific for their decomposition activities or fungal plant pathogens. Sequential rotation of the three crops restricted the growth of several of the indicator plant pathogens.
Conclusions
Long-term fertilization and crop rotation significantly altered microbial community structure in the soil. While fertilization affected soil microorganisms mainly through changes in nutrient profile, crop rotations lead to the attraction and repulsion of specific plant pathogens. Such changes in soil microbial communities need to be considered when planning soil management.
Background
The soil microorganisms represent an integral component determining the soil quality and fertility [1]. They play a pivotal role in biogeochemical cycling [2], soil health [3], plant production [4], or disease suppression efficacy [5]. The soil microbial communities are prone to changes in physical, biological, and chemical soil properties [6,7,8] and have a profound influence on plant development [9]. Any manipulation with these properties directly shapes the structure of microbial community and has a crucial impact on soil functioning and crop production as a consequence [10]. Therefore, the microbial community in soil reflects the impact of soil management on the overall ecosystem and should be considered when planning soil management.
Fertilization aims to maximize crop production by creating more favorable conditions for plant cultivation. The market supply includes a wide range of different kinds of fertilizers, from chemical (mineral or synthetic), biological (microbial inocultants) to organic (sewage sludge or manure), or their combinations. Chemical fertilizers are applied on the soil to address the specific needs of treated soil and cultivated crop [11]. Common chemical fertilizers aim to increase the availability primarily of nitrogen (N), and then also of phosphate (P), potassium (K), and sulphur (S). The provided nutrients are simple compounds so they can be quicky utilized by plants (e.g., ammonium nitrate, ammonium sulphate, potassium chloride, calcium nitrate, elemental sulphur, phosphate rock, superphosphate, or urea) [12]. Contrary to chemical fertilizers, organic fertilizers are complex organic materials that are slowly decomposed in soils [13]. They are commonly produced from organic wastes, such as animal manure, sewage sludge, food processing wastes, municipal solid waste, and food waste [12]. For instance, the annual production of sewage sludge is more than 10 million tDM (tonnes of dry matter), and approximately 50% is used in agriculture in Europe [14, 15]. It is because recycling of such wastes fulfills the concept of “circular economy” that aims to eliminate the waste by their conversion to a valuable source of nutrients [16, 17].
The impact of fertilization on soil properties varies depending on fertilizer type and its composition [18]. Simple compounds in chemical fertilizers leads to significant increase of crop yields, but also might increase the risk of eutrophication, soil acidification, and decrease microbial biomass, biomass carbon, and microbial respiration [19, 20]. Contrary to that, the nutrients in organic fertilizers are released by microbial activity over a longer period [21]. As a result, the microbial biomass and diversity increases [22, 23] and bacterial and fungal community structures are changed [24, 25]. Furthermore, organic fertilizers augments the soil with vast numbers of microbes [26, 27]. Given the importance of microorganism in agriculture, it is especially important to monitor and evaluate the long-term effects of the fertilizers on soil microbes, and compare their long-term impact within each other.
Crop rotation is another typical agricultural practice that is used to avoid nutrient losses, enhance crop yield and mitigate the impact of soil-borne pathogens on crop production [28, 29]. Similarly to the application of fertilizer, crop rotation practices have been found to affect microbial diversity [30, 31] and microbial structures in soil [32]. The diversity and biomass have been described to increase with an increasing number of rotated crops [33,34,35]. Monoculture or crops growing in short rotations (two or three crops rotated on the same field) mostly lead to a decline in crop yield, which might be attributed to the enrichment of plant pathogens that can persist in the soil for the whole rotation cycle [36]. Long rotations, on the other hand, contribute to the increase of disease suppression [37]. The insertion of non-host crops into long rotations decreases the number of pathogens in the soil by not providing the host materials [36]. Unfortunately, this approach does not ensure the total elimination of soil-borne pathogens from the soil.
Our previous studies have shown how long-term fertilization influences soil enzymatic activity [38] and endophytic communities in potatoes [39]. This study focused on the effects of long-term crop rotations, and organic and chemical fertilizer application on soil bacterial and fungal communities in various soil types. The experimental fields were established at four sites in the Czech Republic that differed in their climate characteristics and soil characteristics. Chernozem, luvisol, and cambisol belong among the most widespread or the most fertile soil groups. Chernozems cover approximately 230 million ha, are among the most productive soil types in the world, and are the tenth most abundant soil group in the EU [40]. Luvisols cover approximately 600 million ha, and cambisols cover 1.500 million ha [41], being the second most abundant and the most abundant, respectively, soil groups in the EU [40]. This work aimed to: (1) determine the influence of long-term fertilization and crop rotation on microbial diversity in soils, (2) discover whether the bacterial and fungal community structures are affected, (3) describe the specific alterations of the communities if there were any, and (4) draw conclusions in the context of the soil chemistry. Since the experimental sites were established in different environments with various soil types, we expect the microbial communities to be associated primarily with these sites. However, we hypothesised that long-term fertilization regimes would influence the soil chemistry and, therefore, become the dominant driver of bacterial and fungal community structures within each experimental site. We also hypothesised that communities in organically-treated soils would differ from the chemically-treated ones because of more complex composition of manure and sewage sludge. Furthermore, the crop rotation is expected to repress specific plant pathogens and, consequently, be significantly associated with microbial community structure.
Material and methods
Experimental design
Four experimental field plots were established in the Czech Republic in 1996 and since then were regularly fertilized with organic and chemical fertilizers. Due to the geographically distinct locations of the field plots, different environmental conditions were achieved. Soils at the field plots differed in their soil type (determined using Casagrande’s areometric method [42]), soil chemistry, and nutrient and trace element contents (Table 1, Additional file 1: Table S1).
Experimental field plots were divided into three sections, Section A, Section B, and Section C. In each of these sections, potato (Solanum tuberosum L.), winter wheat (Triticum aestivum L.) and spring barley (Hordeum vulgare L.) were rotated in that order every three years, but the sections differed in a current grown crop. Another practice was the application of four different fertilizers: (1) sewage sludge (330 kg N/ha; SF), (2) sewage sludge (990 kg N/ha; SF3x), (3) cow manure (330 kg N/ha; MF), and (4) NPK (330–90–300 N-P-K kg/ha; NPK). NPK fertilizer was applied in the form of calcium ammonium nitrate (source of N), triple superphosphate (source of P) and potassium salt (source of K). Non-treated soil was used as a control (CF). Following typical agricultural practice, organic fertilizers (SF, SF3x, MF) were applied every third year in the autumn before plowing, i.e., only before potato planting. Chemical fertilizer (NPK) was applied regularly throughout the three-year rotation cycle. The application rate of fertilizers was always determined based on their nitrogen content (determined by the Kjeldahl method) for the whole three-year period. The sewage sludge used at all experimental sites always originated from a wastewater treatment plant in Prague and was stabilized anaerobically at 55 °C [43]. Anaerobic stabilization of sewage sludge took place in two stages. And the average total retention time (both stages together) was around 30 days. In contrast, different cattle manures provided by local farmers were applied to each plot. Manure originated mainly from deep cattle litter. After removal from the stable, manure was properly stored for 6–14 months in the field storage (depending on the possibilities of the experimental station). Immediately after weighing the required dose of manure at the field storage, manure is applied and plowed after application. No chemical fertilizer was applied to the organically fertilized plots. An illustration of the experimental field plots is displayed in Fig. 1. The spatial distribution of the 15 sub-plots (60 m2 each), each with a unique combination of fertilization treatment and crop rotation, varied at each location.
Schematic illustration of experimental design. In practice, each location had 15 sub-plots (60 m2) with various spatial distribution of the sub-plots. Each sub-plot also included a protection area to prevent fertiliser spreading to other sub-fileds
Three composite samples of bulk soil were taken from each sub-plot. Each sample consisted of six individual sub-samples distributed over a 30-cm radius collected with a soil sampler probe (0.75-in. diameter). The sub-samples were taken from the topsoil layer to a depth of 30 cm and pooled together. The total amount of soil in the composite sample, approximately 60 g, was furhter sieved through a 2 mm mesh and thoroughly mixed. In total, 180 pooled soil samples were collected from all locations in March of 2016. At the sampling time, two sections of each field plot were not vegetated: (1) Section A after potato harvesting, (2) Section C after barley harvesting. The third part, Section B, was vegetated with winter wheat. Thus, the last time the organic fertilization was applied was 0.5 years ago for Section A, 1.5 years ago for Section B, and 2.5 years ago for Section C.
For soil chemistry analyses, the three composite samples, previously sieved through a 2 mm mesh and thoroughly mixed, were compiled, further homogenized, and divided into three technical replicates. For analysis of microbial communities, an aliquot of 0.5 g of each composite sample (three representative samples per subplot) was used in downstream analyses.
Soil chemistry analysis
Soil samples were air-dried and sieved through a 2-mm mesh. Soil pH was determined in a 0.2 mol/L KCl (2:5 w/v) solution [44]. The cation exchange capacity (CEC) was calculated as the sum of Ca, Mg, K, Na, Fe, Mn, and Al extractable in 0.1 mol/L BaCl2 (w/v = 1:20) for 2 h [45]. For the determination of total carbon, sulfur, and nitrogen in soils, a CHNS Vario MACRO cube (Elementar Analysensysteme GmbH, Germany) analyzer was used. In this instrument, roughly 25 mg of the soil was burned in a catalytic furnace, and C and N content was determined by using a thermal conductivity detector. Inorganic–N forms (N–NH4+ and N–NO3−) were determined via a SKALAR SAN PLUS SYSTEM continuous flow segmented analyzer (Skalar, Netherlands).
To determine the available macro- and micronutrient content, as well as those of potentially toxic elements in the soils, the Mehlich III extraction procedure was used. One gram of soil was mixed with 10 ml of the extraction mixture (mixed for 10 min) of the following composition: 0.2 mol/L of acetic acid, 0.25 mol/L of ammonium nitrate, 0.013 mol/L of nitric acid, 0.015 mol/L of ammonium fluoride and 0.001 mol/L of ethylenediaminetetraacetic acid (EDTA) [46]. Inductively coupled plasma-atomic emission spectrometry (ICP-OES) using an Agilent 720 (Agilent Technologies Inc., USA) equipped with a two-channel peristaltic pump, a Struman-Masters spray chamber, and a V-groove pneumatic nebulizer made of inert material was used to determine the content of Cd, Cu, Fe, Mn, Pb, Zn, P, and S in soil extracts. Spectrometry parameters were as follows: 1.2 kW; plasma flow: 15.0 L/min; auxiliary flow: 0.75 L/min; nebulizer flow: 0.9 L/min. Flame atomic absorption spectrometry (F-AAS, Varian 280FS, Varian, Australia; airflow of 13.5 L/min, acetylene flow of 2.2 L/min, burner height of 13.5 cm, and a nebulizer uptake rate of 5 mL/ min) was used to determine the content of Ca, Mg, and K in the extracts. The monitored soil chemical parameters are summarized in Additional file 1: Table S1.
DNA isolation
Metagenomic DNA was extracted from 0.5 g of each sieved soil sample using a FastDNA™ Spin Kit for Soil (MP Biomedicals, USA) and further purified with Genomic DNA Clean and Concentrator™ (ZYMO Research, USA) according to the manufacturers’ instructions. DNA concentration and purity were determined spectrophotometrically using a NanoDrop ND-1000 (NanoDrop Technologies, USA).
16S rRNA gene and ITS region amplification and sequencing
Amplicons were prepared using two sequential polymerase chain reactions (PCRs) with specific primers. For 16S rRNA gene amplicons, 515 forward (5′-GTGYCAGCMGCNGCGG-3′) and 926 reverse (5′-CCGYCAATTYMTTTRAGTTT-3′) primers were used targeting the hypervariable regions V4-V5 [47]. The first 15-µl reaction consisted of 0.02 U/µl KAPA HiFi HotStart ReadyMix (Kapa Biosystems, USA), metagenomic DNA (~ 10 ng), 0.3 µM of each primer (Sigma-Aldrich, USA), and water for molecular biology (Sigma-Aldrich, USA). The temperature program for the first PCR was set as follows: 5 min at 95 °C, 20 s at 98 °C, 28—30 cycles of 15 s at 56 °C, 15 s at 72 °C, and a final extension of 5 min at 72 °C [48]. For the second PCR, which consisted of 8 – 10 cycles, the same primers as in the first PCR modified with adaptor tags and internal barcodes were used. The 25-µl reaction contained 0.02 U/µl KAPA HiFi HotStart ReadyMix (Kapa Biosystems, USA), 1 µM of each primer (Sigma-Aldrich, USA), 0.5 µl of the previous unpurified PCR product used as a DNA template, and water for molecular biology (Sigma-Aldrich, USA). The same temperature program as for the first run of PCR was used for the second PCR, except that the annealing temperature was 56 °C. The amplification of the ITS region followed the same protocol with a few modifications. Primers 5.8S-Fun forward (5′-AACTTTYRRCAAYGGATCWCT-3′) and ITS4-Fun reverse (5′-AGCCTCCGCTTATTGATATGCTTAART-3′) [49] were used and for the second PCR modified with the same adaptor tags and barcodes as the primers targeted for the 16S rRNA gene. The reactions for the first and second PCRs were prepared according to the same procedure as for 16S rRNA gene amplification. The temperature program for the first PCR was following: 5 min at 95 °C, 20 s at 98 °C, 28 – 30 cycles of 15 s at 50 °C, 15 s at 72 °C, and a final extension of 5 min at 72 °C. For the second PCR, the temperature program was the same as for 16S rRNA amplicons.
The resultant amplicons of 16S rRNA genes and ITS regions were purified using SPRIselect magnetic beads (Beckman Coulter, USA), eluted into water for molecular biology and sent on ice packs to the Core Facility for Nucleic Acid Analysis at the University of Alaska Fairbanks (AK, USA), where the sequencing analysis was performed using an Illumina Miseq platform as described earlier [48]. Briefly, the DNA concentration was normalized to 1 – 2 ng/µl using a SequalPrep Kit (Thermo Fisher Scientific, USA). With the 16S rRNA gene amplicons, the amplicons were mixed, and non-diluted, 1.5-fold diluted, and threefold diluted technical replicates were prepared and subjected to sequencing. With the ITS region amplicons, no technical replicates from sequencing were conducted, only non-diluted samples were used after the normalization of DNA concentration.
In order to maintain the proper parameters for sequence data processing, amplicons of mock communities were prepared by mixing the genomic DNA of selected strains and processed along with the amplicons of soil samples. The preparation procedure of the 16S rRNA gene mock community is described in a previous study [50], with the exception that the mock community consisted of genomic DNA of the following 12 strains: Arthrobacter chlorophenolicus A6, Achromobacter xylosoxidans A8, Bacillus pumillus SAFR-032, Micrococcus luteus NCTC 2665, Methylobacterium radiotolerans JCM 2831, Pseudomonas alcaliphila JAB1, Pseudomonas veronii 20a2, Rhodobacter jostii RHA1, Burkholderia xenovorans LB400, Cupriavidus necator H850, Pandoraea pnomenusa B-356, and Rhizobium radiobacter C58. The fungal mock community was prepared from the genomic DNA of Saccharomyces cerevisiae DBM 2101, Rhodotorula mucilaginosa DBM 19, Fusarium culmorum DBM 4044, Aspergillus niger DBM 4054, Alternaria sp. DBM 4245, Candida intermedia DSK 46, and Penicillium chrysogenum DBM 4062. The genomic DNA from both bacterial and fungal cultures was isolated with a DNeasy Plant Mini Kit (Qiagen, Germany), and DNA concentration was measured using a Qubit fluorimeter and Qubit® dsDNA HS Assay Kit (both Life Technologies, USA). The final mock community was prepared by mixing the genomic DNA from the cultures in equal concentrations. As a negative control 0.5 ml of water for molecular biology (Sigma-Aldrich, St Louis, MO, USA) was used and proccesed by the same manipulation as soil samples (DNA isolation and purification, PCR amplification and storage). The negative controles were done in four replicates.
Sequence data processing
Amplicon sequence variants (ASVs) were obtained from raw sequence data using the package DADA2 [51] in the statistical program R [52], with minor modifications to the recommended DADA2 pipeline [51] based on the results of an analysis of mock communities. Briefly, after the primers’ sequences were trimmed off, the 16S rRNA gene sequences were filtered and trimmed using the following parameters: truncLen = c(247,174), maxN = 0, maxEE = c(2,2), truncQ = 2. The ITS region sequences were filtered and trimmed using the following parameters: trimLeft = c(0, 0), maxN = 0, maxEE = c(2,2), truncQ = 10. Chimeric sequences were detected and removed using the “consensus” method. To further reduce sequencing errors, ASVs of 16S rRNA genes that differed by one base [47, 48] and ASVs of ITS region that differed by up to two bases were merged while keeping the most abundant sequence as valid. With the 16S rRNA gene sequence dataset, technical replicates from sequencing were merged while omitting all sequences that were only present in one of the replicates. The taxonomy assignment was performed with the assignTaxonomy function (minimal bootstrap value 50) using the rdp_train_set_16 [53] database for 16S rRNA gene sequences and the Warcup v2 database [54] for ITS region sequences. To avoid pseudoreplication, ASV abundances of samples taken from each sub-plot with a unique combination of fertilization and crop rotation were summed together. All obtained MiSeq reads were deposited in the NCBI Short Read Archive under the BioProject accession number PRJNA587449.
Statistical analyses of soil chemistry
Highly correlated continuous chemical variables were identified using Pearson’s correlation coefficient (coefficient threshold set to |0.7|). The correlations were graphically visualized using the psych package [55] in R. All chemical variables were standardized by subtracting their mean value and dividing by their standard deviation. The influence of field location, fertilization, crop rotation, and interaction of these factors on soil chemistry was tested with permutational multivariate analysis of variance (PERMANOVA) based on Euclidean distance. In the final model, the variables were ordered by decreasing R2 values. The soil chemistry did not meet assumptions of normality (according to the Shapiro–Wilk test); hence further analyses were performed with a non-parametric Kruskal–Wallis test with a false discovery rate (FDR) correction of p-values [56]. A pairwise Wilcoxon rank-sum test was used to compare the pairs of five fertilization treatments with each other, again with an FDR correction of p-values. The same procedure was followed to test the influence of crop rotation on the soil chemistry.
Multivariate statistical analyses of microbial community data
Further analyses of microbial data were conducted in the phyloseq [57] and vegan [58] packages in R. ASVs with no assigned taxonomy at the phylum level were removed from analyses (3.99% of all 16S rRNA gene reads; 5.67% of all ITS region reads). The analyzed datasets were rarefied to the smallest sample size (1700 for 16S rRNA ASVs and 24,000 for ITS ASVs). For testing null hypotheses, samples representing sub-plots with different combinations of fertilizer and crop rotation were permuted within the localities. Alpha-diversity was assessed by calculating Shannon and Simpson diversity indices [59]. The Kruskal Wallis rank-sum test was used to determine statistically significant differences among the diversity indices of individual treatments. The analysis of microbial diversity changes associated with different fertilization and crop rotation treatments was at first performed on the whole community dataset and then also on the datasets of each location separately.
The ASVs of prokaryotic and fungal communities were merged at the genus level and genera with ten or fewer sequences were removed. A non-metric multidimensional scaling (NMDS) based on the Bray Curtis distance was used as an ordination projection. To analyze the underlying trends affecting microbial community data, a significant correlation of environmental variables with the ordination configuration was conducted using the envfit function (vegan package). The variables that were found to significantly (p < 0.05) correlate with the sample distribution within the ordination space were fitted to the NMDS ordination. The length of the fitted arrows, representing the gradient direction of variables, were scaled by their correlation (square root of R2) so that weaker predictors had shorter arrows.
The hellinger transformation was applied to the final dataset (i.e., abundance values were divided by the total abundance, and then the square root of the result was taken) [60]. Permutational multivariate analysis of variance (PERMANOVA) was used to assess the effect of fertilization and/or crop rotation on microbial community structures. PERMANOVA was based on Bray Curtis distance, and the factors of location, soil type, fertilization, crop rotation, and their interaction were used to design the final model as follows: the test was first performed for each variable separately, then the variables were ordered according to their decreasing R2 values in the final model. The statistical significance of all factors was tested in this model, but the permutations never involved blending together samples from different locations. Pairwise PERMANOVA was used to identify significant differences between pairs of fertilization and crop rotation levels. Bonferroni correction [61] was used to calculate the corrected p-value, and permutations were again restricted within the levels of location factor. Redundancy analysis (RDA) was used to test the hypothesis that the tested factors influence hellinger-transformed microbial community data. The factors analyzed with RDA included location and soil type, fertilization regime, crop rotation, and the interaction between fertilization and crop rotation. While each of the factors was tested separately, the other factors were used as covariates to exclude their effect from the tests. An ANOVA-like permutation test was performed for RDA to assess the significance of the factors. Only those factors with p < 0.05 were applied in a variation partitioning analysis to assess the extent to which they explained the community data. FUNguild was used to assign identified fungal taxa to functional groups [62] followed only with plant pathogens investigation.
Indicator genera analysis
An indicator species-like analysis (multipatt, indicspecies packages) was used to identify indicator genera for each type of treatment [63, 64]. Indicator genera are taxa that are highly associated with a specific environmental condition such as experimental treatments. The measure of association is calculated based on the probability of taxon occurrence given a sample group and probability of sample group assignment upon taxon detection [65]. The false discovery rate (FDR) correction method was applied to p values, and genera with padj < 0.1 were classified as the indicators. The indicator genera were confirmed to be present in a minimum of four samples from the treatment to which they were assigned.
Results
Chemical parameters of soils
The chemical parameters of the soils examined, including soil elements (Ca, Mg, Mn, Fe, C, Cu, S, Zn, N, Cd, Pb, P, K) and pH values, are summarized in Additional file 1: Table S1. The following positive correlations were found among the chemical parameters (Additional file 1: Fig. S1; Pearson’s correlation, coefficient threshold set to |0.7|): (1) pH with Ca and Mg, (2) Ca further correlated with C and Mn, (3) Cu with Cd and Zn, (4) N with C, (5) Fe with P. Negative correlation was observed for Fe with Mn and Ca.
The location of the experimental field plots (R2 = 0.54), fertilization (R2 = 0.27), crop rotation (R2 = 0.03) and the interaction of fertilization with location (R2 = 0.06) were all found to significantly influence the soil chemistry (p < 0.05, PERMANOVA). While pH significantly differed across the locations and soil types (Kruskal–Wallis test, padj < 0.05), no significant association was found between fertilization treatment and soil pH. Similarly, the levels of almost all soil elements significantly differed across the locations and soil types (Kruskal–Wallis test, padj < 0.05), and the fertilization regime was only significantly associated with 7 out of the 13 soil elements determined (Kruskal–Wallis test, padj < 0.05) (for details, see Additional file 1: Table S2). Compared to the control soil (CF), significantly higher concentrations of Cu, S, Zn, Fe, P, and N were found in SF3x treatment; Zn, Fe, P were significantly increased in SF treatment, and P with K were significantly higher in MF treatment (pairwise Wilcoxon test, padj < 0.05). In total, six elements were significantly increased in multiple treated soils, and the highest concentration of all these elements was measured in SF3x. Only one element (K) was significantly influenced by one treatment (MF), in which it reached the highest concentration in contrast to other treatments. No significant difference of soil chemistry between CF and NPK treatments was observed. Crop rotation was only significantly associated with the concentration of Cd (Kruskal–Wallis test, padj < 0.05).
Microbial diversity
In total, 6,405 unique 16S rRNA gene ASVs and 5,878 unique ITS ASVs were obtained from all samples. Prokaryotic and fungal Shannon and Simpson alpha-diversity indices (Fig. 2) were not significantly associated with fertilization, crop rotation or location (Kruskal–Wallis test, p > 0.05).
Shannon and Simpson diversity indices calculated from prokaryotic (a) and fungal (b) sequence data according to different fertilization and crop rotation treatments; (fertilizer): control (CF), cow manure (MF, 330 kg N/ha), NPK (NPK, 330–90-330 kg/ha), sewage sludge (SF, 330 kg N/ha), sewage sludge (SF3x, 990 kg N/ha); (crop): Section A (after potato), Section B (wheat), Section C (after barley)
Influence of soil chemistry on microbial communities
NMDS ordinations of prokaryotic (stress 0.084, Fig. 3) and fungal communities (stress 0.095, Fig. 4) demonstrated that samples clustered according to the location and soil type. Both of these parameters were significantly correlated with the configuration of the ordination (location: R2 = 0.90 for both communities, p < 0.001; soil type: R2 = 0.93 and 0.87 for prokaryotes and fungi respectively, p < 0.001). Almost all soil chemistry parameters significantly correlated with the sample distribution as well (Figs. 3, 4). The strongest predictors for both prokaryotic and fungal communities were soil pH (R2 = 0.88 for prokaryotes and R2 = 0.85 for fungi, p < 0.001) and Ca (R2 = 0.88 for prokaryotes and R2 = 0.87 for fungi, p < 0.001). C, Fe, Mg, and Mn were other strong predictors for both prokaryotic and fungal communities (R2 from 0.43 to 0.72, p < 0.001). The only difference between predictors of the prokaryotic and fungal community was that Pb was significantly associated with prokaryotic community structure (R2 = 0.22, p < 0.001) and K only with fungal community structure (R2 = 0.19, p < 0.001). P was associated (p > 0.1) with neither prokaryotic nor fungal soil community structures.
Non-metric multidimensional scaling (NMDS) ordination (stress = 0.084) of samples based on composition of soil prokaryotic community (ASVs). Samples are coded by the location (color) and soil type (symbols) they originated from. Arrows represent soil chemical parameters that significantly (p < 0.05) correlated with the ordination configuration
Non-metric multidimensional scaling (NMDS) ordination (stress = 0.095) of samples based on composition of soil fungal community (ASVs). Samples are coded by the location (color) and soil type (symbols) they originated from. Arrows represent soil chemical parameters that significantly (p < 0.05) correlated with the ordination configuration
Influence of fertilization and crop rotation on microbial communities
The influence of fertilization and crop rotation on microbial communities was determined on clustered data of 212 prokaryotic and 402 fungal genera. Both prokaryotic and fungal community structures were significantly associated with crop rotation and fertilization regimes, after partialling out the influence of location and soil type (PERMANOVA, p < 0.001). The interaction effect of crop rotation and fertilization was not statistically significant on either community structure. Compared to control soil (CF), the prokaryotic community was significantly different in soils treated with SF3x and MF (Pairwise PERMANOVA, padj < 0.001), and the fungal community was significantly different (Pairwise PERMANOVA, padj < 0.05) in all fertilized soils (SF, SF3x, MF and NPK). These results correlate with RDA ordination (Fig. 5a) showing that prokaryotic communities were the most different in MF- and SF3x-treated soils from the CF soil, while fungal community structures differed in all fertilized soils. All crops included in rotations were significantly associated with both bacterial and fungal communities (Pairwise PERMANOVA, padj < 0.05) which is also shown in RDA ordination (Fig. 5b). Variation partitioning analysis showed that the influence of fertilization accounted for 5% of the variation in prokaryotes and 10% in fungi (RDA, p < 0.01), whereas the impact of crop rotations was 3% for both prokaryotic and fungal communities (RDA, p < 0.01). Location and soil type were also included in the variation partitioning analysis to determine their influence. Their significant impact (RDA, p < 0.001) on prokaryotic and fungal communities was 48% and 30%, respectively. Using FUNGuild, 6% of identified fungal taxa were assigned to the plant pathogen functional group. Despite location was the main driver of the structure of fungal plant pathogens in bulk soil (R2 = 57%, padj < 0.01, PERMANOVA), fertilizer application and crop rotation had also a significant effect on the structure (R2 = 4%, padj < 0.01). Specifically, a significant difference in the plant pathogen community structure was observed between CF and SF3x fertilization treatments (padj < 0.05, Pairwise PERMANOVA) and all pair comparisons of crop rotation systems (padj < 0.05).
Redundancy analysis (RDA) ordinations based on Bray–Curtis distance. Subfigures represents: a prokaryotic and fungal community structure in soils collected from different fertilization regimes: control (CF), manure (MF; 330 kg N/ha), NPK (NPK; 330–90-330 kg/ha), sewage sludge (SF; 330 kg N/ha), sewage sludge (SF3x; 990 kg N/ha); b prokaryotic and fungal community structure in soils after different crop rotation: Section A (samples were collected after potato harvesting), Section B (samples were collected from field vegetated with wheat), Section C (samples were collected after barley harvesting)
Analysis of indicator genera
Indicator genera associated with one or more types of fertilization regimes or crop rotation treatment are listed in Table 2. Four prokaryotic and three fungal genera were determined to be indicators of the manure amendment (MF), which included Thermoflavimicrobium, Halocella, Clostridium cluster XlVa (Clostridium coccoides group), Erythrobacter, Mycosymbioces, Remersonia, and Mycothermus. Three prokaryotic genera, Rhodanobacter, Micromonospora, and Candidimonas, were found to be indicators of SF3x, whereas no indicator genera were found for SF. However, two prokaryotic (Pseudoxanthomonas, Kaistia) and two fungal genera (Cephaliophora, Leucosporidium) were found to be indicators of both SF and SF3x. One prokaryotic (Romboutsia) and one fungal (Trichophyton) genus were identified as indicators of all organic fertilization regimes (i.e., MF, SF, and SF3x). Pseudaleuria was found to be an indicator of all types of fertilization regimes (MF, SF, SF3x, and NPK).
Ten fungal genera and no prokaryotic genera were found to be indicators of different phases of three-year crop rotation (Table 3).
Discussion
This study is focused on changes in soil microbial communities after 20 years of regular fertilization and crop rotation practices. The major merit of this study lies in evaluating the influence of both agricultural practices on three common soil types used for agricultural purposes and different locations. The experimental field plots were settled over 100 km apart, so large-scale gradients of environmental conditions were achieved. Temperature, cation exchange capacity, oxidizable carbon and clay content differed among the field plots (Table 1), and these factors were previously found to shape the soil microbial community [66, 67]. The soil chemistry also significantly differed across the experimental field plots, with pH being the strongest predictor of both fungal and prokaryotic community structures (Figs. 3, 4). In fact, the pH gradient was almost parallel with the first axis of NMDS ordination, reaching the highest values in chernozem soil. Soil pH was previously reported to be one of the key edaphic determinants of soil microbial community structure [68,69,70]. It plays, among others, a crucial role in nutrient availability for plants and the solubility of soil elements [71], and is therefore often used as a soil quality indicator [72]. Neutral soil pH, which was found in luvisol and chernozem, is mostly reported to be the optimum for nutrient availability [71].
Since location and soil type, hence soil chemistry and environmental conditions, were the major determinants of microbial community structure, they were always controlled for in further analyses when determining the influence of fertilization and crop rotation on microbial communities. Separating out the major determinants of the soil community helped us to obtain insights into the influence of agricultural practices, regardless of what other environmental factors shaped the community. Although only three soil types were included in the experiment, they are among the most widespread or fertile soils in the world [41], which makes this research globally relevant for agricultural systems.
Broadly used alpha-diversity indices, such as Shannon or Simpson, have been reported to be significantly affected after long-term fertilization [23, 73,74,75]. In our study, there was no evidence of a significant influence of location, soil type, long-term fertilization or crop rotation on microbial diversity, which is in agreement with other studies that found no significant effect of fertilization on either prokaryotic [76, 77] or fungal diversity [24].
However, the absence of significant changes in microbial diversity does not necessarily reflect the stability of the microbial community structure [78]. In our study, both fertilization and crop rotation were found to be significantly associated with the microbial community structures, both prokaryotic and fungal (Fig. 5). Since the soil chemistry was significantly altered by fertilization, we assume that the impact of fertilization on microbial succession was more likely indirect: first affecting soil nutrients, which subsequently affect the communities. This suggestion is in agreement with the results of other study [79]. The significant changes in soil chemistry were found in fertilized soils, especially after the application of manure (MF) or higher doses of sewage sludge (SF3x). While the manure increased concentrations of K and P, the macronutrients essential for plant growth [80], sewage sludge (SF3x) was found to be also associated with higher levels of N, P, S and heavy metals and micronutrients, such as Cu, Zn, Fe. These metals were previously detected in sewage sludge [81], and their presence raised concerns about using sewage sludge as a soil amendment [82]. It is important to note, though, that in our study, the concentration of these elements was below the risk limits [83]. Even though these heavy metals might pose a health risk at high concentrations, they are essential to plant health and growth [84].
Although the soil chemistry was significantly associated with fertilization treatments, pH was not one of the altered parameters. This finding is in contrast with our expectations and with the results of several other studies [7, 85], which showed that long-term fertilization resulted in an altered soil pH. It is important to note, though, that organic fertilizers (MF, SF, SF3x) were only applied once every three years, and the sample collection was performed at least half a year after the application of the last organic fertilizer. Therefore, we assume that temporal pH changes may occur directly upon fertilization, but did not last in the long-term. Compared to organic fertilizers, NPK application was performed periodically every year, however, this treatment was not found to significantly influence any of the monitored soil chemical parameters. Bearing that in mind, we assume that NPK doses were too low to impact the soil chemistry in the long-term.
Organic fertilizers are also a source of other substances, such as humic and fulvic acids [86] and pollutants [43, 87] when used as fertilizers. Bearing that in mind, organic fertilizers may stimulate the growth of specific indigenous populations [22] involved in the degradation of these compounds [88]. The direct influence of, in particular, organic fertilizers on soil microbial communities cannot be ruled out either; organic fertilizers contain their own microbiota [27] which are directly transferred to the soil with the application of the fertilizer [89].
Our results indicated that microbial community structure was significantly associated with the fertilization treatment, despite the fact that the final fertilizer application was performed half a year before sampling. These findings suggest that long-term fertilization has a persistent influence on both prokaryotic and fungal communities, although the impact may be different on these two communities (Fig. 5). In fact, while prokaryotic communities only significantly differed from the control soil (CF) with SF3x and MF treatments, fungal communities were significantly different from all fertilization treatments tested. Because SF3x and MF were the treatments in which the soil chemistry significantly differed the most from the control soil (CF), we conclude that the prokaryotic communities, rather than fungi, are more likely to be indirectly influenced by the fertilizers and reflect the current nutrient composition of soils. This finding is in agreement with another study [90] suggesting that nutrients, independently of other edaphic factors, strongly influence the bacterial community structure in the soil. The response of fungi, which generally have tenfold lower growth rates than bacteria [91], to soil chemistry changes is then much slower [92]. Hence, the fungal community structure does not reflect the current nutrient composition in the soil as the bacterial community.
Long-term fertilization in this study was significantly associated with the presence of specific taxa (Table 2), the majority of which are involved in nutrient cycling in soils. The fungus Pseudaleuria, which was found to be an indicator of all fertilization regimes, was previously associated with healthy and disease-suppressive soils [93, 94], thus emphasizing the role of fertilization regimes in the improvement of soil fitness. Thermoflavimicrobium, Halocella, and Clostridium cluster XIVa, all of which were found to be indicators of MF, are typical bacterial fermenters involved in the anaerobic digestion of sugars, lipids, and proteins [95], which are abundant in manure. The presence of the genus Halocella in manured soil was already described in other studies [96, 97], in which this genus was found to be the most dominant, together with bacilli. Erythrobacter was found in the later phases of the composting process of manure, suggesting that it can survive the higher temperatures that occur during compost maturation [98]. Similarly, the fungal genera Mycosymbioces, Remersonia, and Mycothermus, which were also found to be significantly associated with manure application, produce hemicellulosic hydrolytic and related enzymes [99, 100].
Sewage sludge was previously reported to contain several pollutants, such as polychlorinated biphenyls, polybrominated diphenyl ethers, polychlorinated dibenzodioxins, dibenzofurans, polycyclic aromatic hydrocarbons (PAHs), polyhalogenated organic compounds (hexachlorohexane, hexabromocyclododecane), heavy metals [43, 87, 101] and micropollutants, such as antibiotics [102, 103], pharmaceuticals [104] or synthetic musks [105]. These chemicals may be deposited in sewage sludge from the disposal of household chemicals, human metabolic waste, or urban runoff [106]. Nevertheless, the application of sewage sludge introduces a lot of important nutrients into arable soil, together with strains capable of degrading the pollutants [107]. In this respect, taxa commonly associated with the degradation of a wide range of organic pollutants were found to be significantly enriched in soils amended with sewage sludge. The enriched taxa included Rhodanobacter, which has the ability to degrade halogenated pollutants [108, 109] and high-molecular-weight PAHs [110]; Micromonospora, Pseudoxanthomonas and Kaistia, whose biodegradation roles are mostly linked with the degradation of aromatic [111] and polyaromatic hydrocarbons [112,113,114] as well as nitrogen-containing organic pollutants [115]; or Micromonospora, which is involved in the degradation of petroleum hydrocarbons [113] and other chemicals [116]. Additionaly, Micromonospora is also beneficial to plant health and growth. It can stimulate plant immunity by enhancing the jasmonate-regulated self-defense system [117]. Despite the fact that around 50% of sewage sludge is used as a fertilizer in Europe [14], European legislation sets the limits for metals in soil [83], while the limits for organic contaminants are not included, and the Member States have issued their own norms [14, 118]. Importantly, the concentrations of organic pollutants in our experimental soils were not significantly associated with the fertilization regime [119].
Crop rotation is an agricultural practice often used to sustain soil fertility, since the continuous planting of one crop decreases the nutrient content and promotes the build-up of host-specific plant pathogens [36]. The incorporation of other crops with different nutrient demands into the rotations helps restore the soil chemistry and decreases the number of soil-borne pathogens [36, 120]. Although microbial diversity has been reported to depend on the currently planted crop and to alter according to the order of crops rotated [30], no significant differences in diversity indices were detected in our study. As the same plants were rotated in the same order for a continuous period of 20 years, we concluded that a consistent history of long-term crop rotations contributes to the stabilization of microbial diversity. The other factor which could infuenced the unchanged diversity is the sampling the bulk soil.
Although no significant changes in alpha-microbial diversity were detected, the community structure was significantly associated with current crop rotation. Plants are able to alter the soil microbial community composition via root exudates and other rhizodeposits [8, 32, 121] which differ among plant species, cultivars or even their developmental stage [122]. Our finding indicates that the plant secondary metabolites may persist in the soil even after crop-harvesting, and together with plant residues influence the soil microbial community structure. However, further studies are needed to reveal a potential persistent influence of rhizodeposits on soil microbial communities.
A significant association between crop rotation and microbial communities was also shown by identifying indicator microorganisms that were typical of each phase of rotation. These indicators consisted only of fungi; no specific bacterial genera were found to be significantly associated with rotations (Table 3). Not surprisingly, the overwhelming majority of the observed indicator fungi were plant pathogens, except for the arbuscular mycorrhizal fungus Glomus [123], which was found in soil after potato harvesting. The genus Bipolaris, which was an indicator also found in soil after potato harvesting, includes barley and wheat hemibiotrophs [124] and is able to create dormant forms persisting in soil [125]. Blumeria, which was found to be associated with wheat growth, is an obligate biotrophic pathogen causing a powdery mildew of cereals, wheat included [126], and is able to produce its survival form, chasmothecia [127]. Zymoseptoria, detected during the wheat phase, is a biotrophic wheat pathogen [128], and Gaeumannomyces, detected after barley harvesting, is a soil-borne pathogen that causes the „take-all “ disease of wheat [129]. These results showed that the sequential rotation of the three crops restricted the growth of several pathogens, such as Zymoseptoria and Blumeria, over long periods, and very likely hinders outbreaks of them. However, it seems that Bipolaris and Gaeumannomyces were able to persist in soil even during the crop rotations. Since all crop rotational systems significantly shaped the profile of fungal plant pathogens in soil, it can be assumed that it is the plant secondary metabolites and plant residues that attract or repeal plant pathogens. Prolonging of the rotation cycle by the insertion of more plants might, therefore, be needed to suppress more specific plant pathogens.
Conclusion
Our study is one of the few that investigates the effect of regular fertilization and crop rotation on soil microbial structure and diversity over the long term. The fact that the experiment was performed in multiple geographic locations with different soil types increases its general insights into how these common agricultural practices influence the soil microbiome under diverse environmental conditions. Soil chemistry that differed among the experimental filed plots was the major determinant of microbial community structures, with pH being the strongest predictor. Within the experimental field plots, the soil chemistry was further influenced by fertilization and crop rotation, which consequently modified the bacterial and fungal communities. However, the response of prokaryotic and fungal communities to long-term fertilization treatments differed. Prokaryotes significantly differed from CF only in SF3x and MF treatments while fungal communities were significantly shifted from CF by each fertilization treatment. The crop rotations significantly altered both prokaryotic and fungal communities, with the main effect on the relative abundance of plant-pathogens. In summary, our study provides deeper insights into plant-fertilizer-microbe interactions in soil, and our findings can potentially serve to improve our understanding of soil quality and globally contribute to sustainable agriculture.
Availability of data and materials
Raw data of 16S rRNA gene amplicon sequences supporting the findings of the present study are available in the Sequence Read Archive of NCBI under BioProject accession PRJNA587449.
References
Reeves DW. The role of soil organic matter in maintaining soil quality in continuous cropping systems. Soil Till Res. 1997;43:131–67.
Falkowski PG, Fenchel T, Delong EF. The microbial engines that drive earth’s biogeochemical cycles. Science. 2008;320:1034–9. https://doi.org/10.1126/science.1153213.
Chaparro JM, Sheflin AM, Manter DK, Vivanco JM. Manipulating the soil microbiome to increase soil health and plant fertility. Biol Fertil Soils. 2012;48:489–99. https://doi.org/10.1007/s00374-012-0691-4.
van der Heijden MGA, Bardgett RD, van Straalen NM. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett. 2008;11:296–310. https://doi.org/10.1111/j.1461-0248.2007.01139.x.
Hu J, Wei Z, Friman V-P, Gu S, Wang X, Eisenhauer N, Yang T, Ma J, Shen Q, Xu Y, et al. Probiotic diversity enhances rhizosphere microbiome function and plant disease suppression. MBio. 2016. https://doi.org/10.1128/mBio.01790-16.
Cleveland CC, Liptzin D. C:N: P stoichiometry in soil: is there a “redfield ratio” for the microbial biomass? Biogeochemistry. 2007;85:235–52. https://doi.org/10.1007/s10533-007-9132-0.
Zhang Y, Shen H, He X, Thomas BW, Lupwayi NZ, Hao X, Thomas MC, Shi X. Fertilization shapes bacterial community structure by alteration of soil PH. Front Microbiol. 2017. https://doi.org/10.3389/fmicb.2017.01325.
Dennis PG, Miller AJ, Hirsch PR. Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities? Root exudates and rhizosphere bacteria. FEMS Microbiol Ecol. 2010;72:313–27. https://doi.org/10.1111/j.1574-6941.2010.00860.x.
Gaiero JR, McCall CA, Thompson KA, Day NJ, Best AS, Dunfield KE. Inside the root microbiome: bacterial root endophytes and plant growth promotion. Am J Bot. 2013;100:1738–50. https://doi.org/10.3732/ajb.1200572.
Seghers D, Wittebolle L, Top EM, Verstraete W, Siciliano SD. Impact of agricultural practices on the Zea mays L. endophytic community. Appl Environ Microbiol. 2004;70:1475–82. https://doi.org/10.1128/AEM.70.3.1475-1482.2004.
Hermida L, Agustian J. Slow release urea fertilizer synthesized through recrystallization of urea incorporating natural bentonite using various binders. Environ Technol Innov. 2019;13:113–21. https://doi.org/10.1016/j.eti.2018.11.005.
Bünemann EK, Schwenke GD, Van Zwieten L. Impact of agricultural inputs on soil organisms—a review. Soil Res. 2006;44:379. https://doi.org/10.1071/SR05125.
Golabi MH, Denney MJ, Iyekar C. Value of composted organic wastes as an alternative to synthetic fertilizers for soil quality improvement and increased yield. Compost Sci Util. 2007;15:267–71. https://doi.org/10.1080/1065657X.2007.10702343.
Collivignarelli A, Frattarola CM, Padovani K. Torretta legislation for the reuse of biosolids on agricultural land in Europe: overview. Sustainability. 2019;11:6015. https://doi.org/10.3390/su11216015.
Eurostat (Overview, 2019) Sewage sludge production and disposal from urban wastewater. https://Ec.Europa.Eu/Eurostat/Web/Environment/Water.
Chew C, Yen N, Ho Y. Show transformation of biomass waste into sustainable organic fertilizers. Sustainability. 2019;11:2266. https://doi.org/10.3390/su11082266.
EUR-Lex Directive EU/2018/851 EUR-Lex Directive EU/2018/851 of the European Parliament and of the Council of 30 May 2018 Amending Directive 2008/98/EC on Waste. Off. J. Eur. Communities 2018, 150, 109–140
Blanchet G, Gavazov K, Bragazza L, Sinaj S. Responses of soil properties and crop yields to different inorganic and organic amendments in a swiss conventional farming system. Agr Ecosyst Environ. 2016;230:116–26. https://doi.org/10.1016/j.agee.2016.05.032.
Minakova EA, Shlichkov AP, Arinina AV. Approaches to management of anthropogenic eutrophication caused by loading from mineral fertilizers. IOP Conf Ser Earth Environ Sci. 2019;272:032006. https://doi.org/10.1088/1755-1315/272/3/032006.
Zhang T, Chen HYH, Ruan H. Global negative effects of nitrogen deposition on soil microbes. ISME J. 2018;12:1817–25. https://doi.org/10.1038/s41396-018-0096-y.
Nyamangara J. Plant and environment interactions. J Environ Qual. 2003;32:8.
Chu H, Lin X, Fujii T, Morimoto S, Yagi K, Hu J, Zhang J. Soil microbial biomass, dehydrogenase activity, bacterial community structure in response to long-term fertilizer management. Soil Biol Biochem. 2007;39:2971–6. https://doi.org/10.1016/j.soilbio.2007.05.031.
Ge Y, Zhang J, Zhang L, Yang M, He J. Long-term fertilization regimes affect bacterial community structure and diversity of an agricultural soil in Northern China. J Soils Sediments. 2008;8:43–50. https://doi.org/10.1065/jss2008.01.270.
Francioli D, Schulz E, Lentendu G, Wubet T, Buscot F, Reitz T. Mineral vs. organic amendments: microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front Microbiol. 2016. https://doi.org/10.3389/fmicb.2016.01446.
Lentendu G, Wubet T, Chatzinotas A, Wilhelm C, Buscot F, Schlegel M. Effects of long-term differential fertilization on eukaryotic microbial communities in an arable soil: a multiple barcoding approach. Mol Ecol. 2014;23:3341–55. https://doi.org/10.1111/mec.12819.
Gomez E, Ferreras L, Toresani S. Soil bacterial functional diversity as influenced by organic amendment application. Biores Technol. 2006;97:1484–9. https://doi.org/10.1016/j.biortech.2005.06.021.
Stiborova H, Wolfram J, Demnerova K, Macek T, Uhlik O. Bacterial community structure in treated sewage sludge with mesophilic and thermophilic anaerobic digestion. Folia Microbiol. 2015;60:531–9. https://doi.org/10.1007/s12223-015-0396-9.
Berzsenyi Z, Győrffy B, Lap D. Effect of crop rotation and fertilisation on maize and wheat yields and yield stability in a long-term experiment. Eur J Agron. 2000;13:225–44. https://doi.org/10.1016/S1161-0301(00)00076-9.
Dixon GR, Tilston EL, editors. Soil microbiology and sustainable crop production. Dordrecht: Springer; 2010.
Larkin RP, Honeycutt CW. Effects of different 3-year cropping systems on soil microbial communities and rhizoctonia diseases of potato. Phytopathology. 2006;96:68–79. https://doi.org/10.1094/PHYTO-96-0068.
Xuan DT, Guong VT, Rosling A, Alström S, Chai B, Högberg N. Different crop rotation systems as drivers of change in soil bacterial community structure and yield of rice, Oryza sativa. Biol Fertil Soils. 2012;48:217–25. https://doi.org/10.1007/s00374-011-0618-5.
Soman C, Li D, Wander MM, Kent AD. Long-term fertilizer and crop-rotation treatments differentially affect soil bacterial community structure. Plant Soil. 2017;413:145–59. https://doi.org/10.1007/s11104-016-3083-y.
McDaniel MD, Grandy AS, Tiemann LK, Weintraub MN. Crop rotation complexity regulates the decomposition of high and low quality residues. Soil Biol Biochem. 2014;78:243–54. https://doi.org/10.1016/j.soilbio.2014.07.027.
Tiemann LK, Grandy AS, Atkinson EE, Marin-Spiotta E, McDaniel MD. Crop rotational diversity enhances belowground communities and functions in an agroecosystem. Ecol Lett. 2015;18:761–71. https://doi.org/10.1111/ele.12453.
Venter ZS, Jacobs K, Hawkins H-J. The impact of crop rotation on soil microbial diversity: a meta-analysis. Pedobiologia. 2016;59:215–23. https://doi.org/10.1016/j.pedobi.2016.04.001.
Bennett AJ, Bending GD, Chandler D, Hilton S, Mills P. Meeting the demand for crop production: the challenge of yield decline in crops grown in short rotations. Biol Rev. 2012;87:52–71. https://doi.org/10.1111/j.1469-185X.2011.00184.x.
Peralta AL, Sun Y, McDaniel MD, Lennon JT. Crop rotational diversity increases disease suppressive capacity of soil microbiomes. Ecosphere. 2018;9:16.
Kracmarova M, Kratochvilova H, Uhlik O, Strejcek M, Szakova J, Cerny J, Tlustos P, Balik J, Demnerova K, Stiborova H. Response of soil microbes and soil enzymatic activity to 20 years of fertilization. Agronomy. 2020;10:1542. https://doi.org/10.3390/agronomy10101542.
Kracmarova M, Karpiskova J, Uhlik O, Strejcek M, Szakova J, Balik J, Demnerova K, Stiborova H. Microbial communities in soils and endosphere of Solanum tuberosum L. and their response to long-term fertilization. Microorganisms. 2020;8:1377. https://doi.org/10.3390/microorganisms8091377.
Tóth, G.; Montanarella, L.; Stolbovoy, V.; Máté, F.; Bódis, K.; Jones, A.; Panagos, P.; Van Liedekerke, M.; European Commission; Joint Research Centre; et al. Soils of the European Union.; Publications Office: Luxembourg, 2008; ISBN 978-92-79-09530-6.
Status of the World’s Soil Resources: Main Report.; Food and Agriculture Organization of the United Nationsand Intergovernmental Technical Panel on Soils.; FAO and ITPS: Rome, Italy, 2015; ISBN 978-92-5-109004-6.
Rząsa S, Owczarzak W. Methods for the granulometric analysis of soil for science and practice. Pol J Soil Sci. 2013;46:50.
Stiborova H, Kolar M, Vrkoslavova J, Pulkrabova J, Hajslova J, Demnerova K, Uhlik O. Linking toxicity profiles to pollutants in sludge and sediments. J Hazard Mater. 2017;321:672–80. https://doi.org/10.1016/j.jhazmat.2016.09.051.
Zbíral J. Comparison of methods for soil PH determination. Rostl Výr. 2001;47:463–7.
Standard of Soil Quality – Determination of Effective Cation Exchange Capacity and Base Saturation Level Using Barium Chloride Solution; 1994.
Mehlich A. Mehlich 3 soil test extractant: a modification of Mehlich 2 extractant. Commun Soil Sci Plant Anal. 1984;15:1409–16.
Fraraccio S, Strejcek M, Dolinova I, Macek T, Uhlik O. Secondary compound hypothesis revisited: selected plant secondary metabolites promote bacterial degradation of Cis-1,2-dichloroethylene (CDCE). Sci Rep. 2017. https://doi.org/10.1038/s41598-017-07760-1.
Lopez-Echartea E, Strejcek M, Mateju V, Vosahlova S, Kyclt R, Demnerova K, Uhlik O. Bioremediation of chlorophenol-contaminated sawmill soil using pilot-scale bioreactors under consecutive anaerobic-aerobic conditions. Chemosphere. 2019;227:670–80. https://doi.org/10.1016/j.chemosphere.2019.04.036.
Taylor DL, Walters WA, Lennon NJ, Bochicchio J, Krohn A, Caporaso JG, Pennanen T. Accurate estimation of fungal diversity and abundance through improved lineage-specific primers optimized for illumina amplicon sequencing. Appl Environ Microbiol. 2016;82:7217–26. https://doi.org/10.1128/AEM.02576-16.
Uhlik O, Wald J, Strejcek M, Musilova L, Ridl J, Hroudova M, Vlcek C, Cardenas E, Mackova M, Macek T. Identification of bacteria utilizing biphenyl, benzoate, and naphthalene in long-term contaminated soil. PLoS ONE. 2012;7:e40653. https://doi.org/10.1371/journal.pone.0040653.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from illumina amplicon data. Nat Methods. 2016;13:581–3. https://doi.org/10.1038/nmeth.3869.
R Core Team. R: a language and environment for statistical computing in R foundation for statistical computing. Vienna, Austria, 2017.
Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM. Ribosomal database project: data and tools for high throughput RRNA analysis. Nucleic Acids Res. 2014;42:D633–42. https://doi.org/10.1093/nar/gkt1244.
Deshpande V, Wang Q, Greenfield P, Charleston M, Porras-Alfaro A, Kuske CR, Cole JR, Midgley DJ, Tran-Dinh N. Fungal identification using a Bayesian classifier and the warcup training set of internal transcribed spacer sequences. Mycologia. 2016;108:1–5. https://doi.org/10.3852/14-293.
Revelle, W.R. Psych: procedures for personality and psychological research. Software. Photographer 80; 2017.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodol). 1995;57:289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217. https://doi.org/10.1371/journal.pone.0061217.
Oksanen J, Blanchet FG, Kindt R, Legendre P, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan: community ecology package. R Package Version 2.5-6; 2019.
Callahan BJ, Sankaran K, Fukuyama JA, McMurdie PJ, Holmes SP. Bioconductor workflow for microbiome data analysis: from raw reads to community analyses. F1000Res. 2016;5:1492. https://doi.org/10.12688/f1000research.8986.1.
Legendre P, Gallagher EDG. Ecologically meaningful transformations for ordination of species data. Oecologia. 2001;129:271–80.
Cabin RT, Mitchell RT. To bonferroni or not to bonferroni: when and how are the questions. Bull Ecol Soc Am. 2000;81:246–8.
Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, Schilling JS, Kennedy PG. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016;20:241–8. https://doi.org/10.1016/j.funeco.2015.06.006.
Cáceres M, Jansen F. Package “indicspecies”: relationship between species and groups of sites (Version 1.7.6). 2016. https://cran.r-project.org/web/packages/indicspecies/indicspecies.pdf
Cáceres MD, Legendre P, Moretti M. Improving indicator species analysis by combining groups of sites. Oikos. 2010;119:1674–84. https://doi.org/10.1111/j.1600-0706.2010.18334.x.
Cáceres MD, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566–74. https://doi.org/10.1890/08-1823.1.
Lladó S, López-Mondéjar R, Baldrian P. Drivers of microbial community structure in forest soils. Appl Microbiol Biotechnol. 2018;102:4331–8. https://doi.org/10.1007/s00253-018-8950-4.
Xue P-P, Carrillo Y, Pino V, Minasny B, McBratney AB. Soil properties drive microbial community structure in a large scale transect in South Eastern Australia. Sci Rep. 2018;8:11725. https://doi.org/10.1038/s41598-018-30005-8.
Wakelin SA, Macdonald LM, Rogers SL, Gregg AL, Bolger TP, Baldock JA. Habitat selective factors influencing the structural composition and functional capacity of microbial communities in agricultural soils. Soil Biol Biochem. 2008;40:803–13. https://doi.org/10.1016/j.soilbio.2007.10.015.
Waldrop MP, Holloway JM, Smith DB, Goldhaber MB, Drenovsky RE, Scow KM, Dick R, Howard D, Wylie B, Grace JB. The interacting roles of climate, soils, and plant production on soil microbial communities at a continental scale. Ecology. 2017;98:1957–67. https://doi.org/10.1002/ecy.1883.
Xu Z, Zhang T, Wang S, Wang Z. Soil PH and C/N ratio determines spatial variations in soil microbial communities and enzymatic activities of the agricultural ecosystems in Northeast China: Jilin Province Case. Appl Soil Ecol. 2020;155:103629. https://doi.org/10.1016/j.apsoil.2020.103629.
Penn CJ, Camberato JJ. A critical review on soil chemical processes that control how soil PH affects phosphorus availability to plants. Agriculture. 2019;9:120. https://doi.org/10.3390/agriculture9060120.
Schloter M, Dilly O, Munch JC. Indicators for evaluating soil quality. Agric Ecosyst Environ. 2003;98:255–62. https://doi.org/10.1016/S0167-8809(03)00085-9.
Hartmann M, Frey B, Mayer J, Mäder P, Widmer F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 2015;9:1177–94.
Luo G, Rensing C, Chen H, Liu M, Wang M, Guo S, Ling N, Shen Q. Deciphering the associations between soil microbial diversity and ecosystem multifunctionality driven by long-term fertilization management. Funct Ecol. 2018;32:1103–16. https://doi.org/10.1111/1365-2435.13039.
Zhao J, Ni T, Li Y, Xiong W, Ran W, Shen B, Shen Q, Zhang R. Responses of bacterial communities in arable soils in a rice-wheat cropping system to different fertilizer regimes and sampling times. PLoS ONE. 2014;9:e85301. https://doi.org/10.1371/journal.pone.0085301.
Fierer N, Lauber CL, Ramirez KS, Zaneveld J, Bradford MA, Knight R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 2012;6:1007–17. https://doi.org/10.1038/ismej.2011.159.
Riber L, Poulsen PHB, Al-Soud WA, Skov Hansen LB, Bergmark L, Brejnrod A, Norman A, Hansen LH, Magid J, Sørensen SJ. Exploring the immediate and long-term impact on bacterial communities in soil amended with animal and urban organic waste fertilizers using pyrosequencing and screening for horizontal transfer of antibiotic resistance. FEMS Microbiol Ecol. 2014;90:206–24. https://doi.org/10.1111/1574-6941.12403.
Hartmann M, Widmer F. Community structure analyses are more sensitive to differences in soil bacterial communities than anonymous diversity indices. Appl Environ Microbiol. 2006;72:7804–12. https://doi.org/10.1128/AEM.01464-06.
Marschner P. Structure and function of the soil microbial community in a long-term fertilizer experiment. Soil Biol Biochem. 2003;35:453–61. https://doi.org/10.1016/S0038-0717(02)00297-3.
Kulcheski FR, Côrrea R, Gomes IA, de Lima JC, Margis R. NPK macronutrients and MicroRNA homeostasis. Front Plant Sci. 2015. https://doi.org/10.3389/fpls.2015.00451.
Baldwin A, Brown TA, Beckett BHT, Elliott GEP. The forms of combination of Cu and Zn in digested sewage sludge. Water Res. 1983;17:1935–44. https://doi.org/10.1016/0043-1354(83)90217-8.
Saha S, Saha BN, Pati S, Pal B, Hazra GC. Agricultural use of sewage sludge in India: benefits and potential risk of heavy metals contamination and possible remediation options—a review. Int J Environ Technol Manag. 2017;20:183–99.
Council Directive 86/278/EEC Council Directive 86/278/EEC on the protection of the environment, and in particular of the soil, when sewage sludge is used in agriculture, OJ L 181, 4.7.1986, pp. 6–12.
Rengel Z. Heavy metals as essential nutrients. In: Prsad MNV, Hagemyer J, editors. Heavy metal stress in plants. Berlin: Springer; 1999. p. 231–51.
Liu J, Liu M, Wu M, Jiang C, Chen X, Cai Z, Wang B, Zhang J, Zhang T, Li Z. Soil PH rather than nutrients drive changes in microbial community following long-term fertilization in acidic ultisols of Southern China. J Soils Sediments. 2018;18:1853–64. https://doi.org/10.1007/s11368-018-1934-2.
Iakimenko O, Otabbong E, Sadovnikova L, Persson J, Nilsson I, Orlov D, Ammosova Y. Dynamic transformation of sewage sludge and farmyard manure components. 1. Content of humic substances and mineralisation of organic carbon and nitrogen in incubated soils. Agric Ecosyst Environ. 1996;58:121–6. https://doi.org/10.1016/0167-8809(95)01006-84.
Clarke RM, Cummins E. Evaluation of “classic” and emerging contaminants resulting from the application of biosolids to agricultural lands: a review. Hum Ecol Risk Assess Int J. 2015;21:492–513. https://doi.org/10.1080/10807039.2014.930295.
Ventorino V, Pascale A, Adamo P, Rocco C, Fiorentino N, Mori M, Faraco V, Pepe O, Fagnano M. Comparative assessment of autochthonous bacterial and fungal communities and microbial biomarkers of polluted agricultural soils of the Terra Dei Fuochi. Sci Rep. 2018;8:14281. https://doi.org/10.1038/s41598-018-32688-5.
Manyi-Loh C, Mamphweli S, Meyer E, Makaka G, Simon M, Okoh A. An overview of the control of bacterial pathogens in cattle manure. IJERPH. 2016;13:843. https://doi.org/10.3390/ijerph13090843.
Knelman JE, Schmidt SK, Lynch RC, Darcy JL, Castle SC, Cleveland CC, Nemergut DR. Nutrient addition dramatically accelerates microbial community succession. PLoS ONE. 2014;9:e102609. https://doi.org/10.1371/journal.pone.0102609.
Rousk J, Baath E. Fungal biomass production and turnover in soil estimated using the acetate-in-ergosterol technique. Soil Biol. 2007;5:2173–7.
Sun S, Li S, Avera BN, Strahm BD, Badgley BD. Soil bacterial and fungal communities show distinct recovery patterns during forest ecosystem restoration. Appl Environ Microbiol. 2017;83:e00966-17. https://doi.org/10.1128/AEM.00966-17.
Gao L, Wang R, Gao J, Li F, Huang G, Huo G, Liu Z, Tang W, Shen G. Analysis of the structure of bacterial and fungal communities in disease suppressive and disease conducive tobacco-planting soils in China. Soil Res. 2020;58:35. https://doi.org/10.1071/SR19204.
Xu L, Ravnskov S, Larsen J, Nilsson RH, Nicolaisen M. Soil Fungal community structure along a soil health gradient in pea fields examined using deep amplicon sequencing. Soil Biol Biochem. 2012;46:26–32. https://doi.org/10.1016/j.soilbio.2011.11.010.
Zhang J, Loh K-C, Lee J, Wang C-H, Dai Y, Wah Tong Y. Three-stage anaerobic co-digestion of food waste and horse manure. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-01408-w.
Toyota K, Kuninaga S. Comparison of soil microbial community between soils amended with or without farmyard manure. Appl Soil Ecol. 2006;33:39–48. https://doi.org/10.1016/j.apsoil.2005.09.002.
Holman DB, Hao X, Topp E, Yang HE, Alexander TW. Effect of co-composting cattle manure with construction and demolition waste on the archaeal, bacterial, and fungal microbiota, and on antimicrobial resistance determinants. PLoS ONE. 2016;11:e0157539. https://doi.org/10.1371/journal.pone.0157539.
Selvam A, Xu D, Zhao Z, Wong JWC. Fate of tetracycline, sulfonamide and fluoroquinolone resistance genes and the changes in bacterial diversity during composting of swine manure. Biores Technol. 2012;126:383–90. https://doi.org/10.1016/j.biortech.2012.03.045.
McPhillips K, Waters DM, Parlet C, Walsh DJ, Arendt EK, Murray PG. Purification and characterisation of a β-1,4-xylanase from Remersonia thermophila CBS 540.69 and its application in bread making. Appl Biochem Biotechnol. 2014;172:1747–62. https://doi.org/10.1007/s12010-013-0640-1.
Wang K, Yin X, Mao H, Chu C, Tian Y. Changes in structure and function of fungal community in cow manure composting. Biores Technol. 2018;255:123–30. https://doi.org/10.1016/j.biortech.2018.01.064.
Stiborova H, Vrkoslavova J, Lovecka P, Pulkrabova J, Hradkova P, Hajslova J, Demnerova K. Aerobic biodegradation of selected polybrominated diphenyl ethers (PBDEs) in wastewater sewage sludge. Chemosphere. 2015;118:315–21. https://doi.org/10.1016/j.chemosphere.2014.09.048.
Chen Q, An X, Li H, Su J, Ma Y, Zhu Y-G. Long-term field application of sewage sludge increases the abundance of antibiotic resistance genes in soil. Environ Int. 2016;92–93:1–10. https://doi.org/10.1016/j.envint.2016.03.026.
Huang K, Tang J, Zhang X-X, Xu K, Ren H. A comprehensive insight into tetracycline resistant bacteria and antibiotic resistance genes in activated sludge using next-generation sequencing. Int J Mol Sci. 2014;15:10083–100. https://doi.org/10.3390/ijms150610083.
Jelic A, Gros M, Ginebreda A, Cespedes-Sánchez R, Ventura F, Petrovic M, Barcelo D. Occurrence, partition and removal of pharmaceuticals in sewage water and sludge during wastewater treatment. Water Res. 2011;45:1165–76. https://doi.org/10.1016/j.watres.2010.11.010.
Stevens JL, Northcott GL, Stern GA, Tomy GT, Jones KC. PAHs, PCBs, PCNs, organochlorine pesticides, synthetic musks, and polychlorinated n-alkanes in UK. Sewage sludge: survey results and implications. Environ Sci Technol. 2003;37:462–7. https://doi.org/10.1021/es020161y.
Pérez S, Guillamón M, Barceló D. Quantitative analysis of polycyclic aromatic hydrocarbons in sewage sludge from wastewater treatment plants. J Chromatogr A. 2001;938:57–65. https://doi.org/10.1016/S0021-9673(01)01338-3.
Schloter M, Nannipieri P, Sørensen SJ, van Elsas JD. Microbial indicators for soil quality. Biol Fertil Soils. 2018;54:1–10. https://doi.org/10.1007/s00374-017-1248-3.
Nalin R, Simonet P, Vogel TM, Normand P. Rhodanobacter lindaniclasticus Gen. Nov., Sp. Nov., a lindane-degrading bacterium. Int J Syst Bacteriol. 1999;49:19–23. https://doi.org/10.1099/00207713-49-1-19.
Uhlik O, Musilova L, Ridl J, Hroudova M, Vlcek C, Koubek J, Holeckova M, Mackova M, Macek T. Plant secondary metabolite-induced shifts in bacterial community structure and degradative ability in contaminated soil. Appl Microbiol Biotechnol. 2013;97:9245–56. https://doi.org/10.1007/s00253-012-4627-6.
Lafortune I, Juteau P, Déziel E, Lépine F, Beaudet R, Villemur R. bacterial diversity of a consortium degrading high-molecular-weight polycyclic aromatic hydrocarbons in a two-liquid phase biosystem. Microb Ecol. 2009;57:455–68. https://doi.org/10.1007/s00248-008-9417-4.
Kim JM, Le NT, Chung BS, Park JH, Bae J-W, Madsen EL, Jeon CO. Influence of soil components on the biodegradation of benzene, toluene, ethylbenzene, and o-, m-, and p-Xylenes by the newly isolated bacterium Pseudoxanthomonas spadix BD-A59. Appl Environ Microbiol. 2008;74:7313–20. https://doi.org/10.1128/AEM.01695-08.
Cui Z, Lai Q, Dong C, Shao Z. Biodiversity of polycyclic aromatic hydrocarbon-degrading bacteria from deep sea sediments of the Middle Atlantic Ridge. Environ Microbiol. 2008;10:2138–49. https://doi.org/10.1111/j.1462-2920.2008.01637.x.
Gojgic-Cvijovic GD, Milic JS, Solevic TM, Beskoski VP, Ilic MV, Djokic LS, Narancic TM, Vrvic MM. Biodegradation of petroleum sludge and petroleum polluted soil by a bacterial consortium: a laboratory study. Biodegradation. 2012;23:1–14. https://doi.org/10.1007/s10532-011-9481-1.
Klankeo P, Nopcharoenkul W, Pinyakong O. Two novel pyrene-degrading Diaphorobacter Sp. and Pseudoxanthomonas Sp. isolated from soil. J Biosci Bioeng. 2009;108:488–95. https://doi.org/10.1016/j.jbiosc.2009.05.016.
Zhang X, Qu Y, You S, Ma Q, Zhou H, Zhang L, Zhang L, Jing J, Liu L. Bioremediation of nitrogen-containing organic pollutants using phenol-stimulated activated sludge: performance and microbial community analysis: phenol-stimulated sludge for n-heterocyclic bioremediation. J Chem Technol Biotechnol. 2018;93:3199–207. https://doi.org/10.1002/jctb.5676.
Vidal CM, Vitale AA, Viale AA. Degradation of naphthalene-2-sulfonate by strains of micromonospora. Rev Argent Microbiol. 1999;31(1):42–8.
Ciancio A, Pieterse CMJ, Mercado-Blanco J. Editorial: harnessing useful rhizosphere microorganisms for pathogen and pest biocontrol. Front Microbiol. 2016. https://doi.org/10.3389/fmicb.2016.01620.
Smith SR. Organic contaminants in sewage sludge (biosolids) and their significance for agricultural recycling. Proc R Soc A. 2009;367:4005–41. https://doi.org/10.1098/rsta.2009.0154.
Pulkrabová J, Černý J, Száková J, Švarcová A, Gramblička T, Hajšlová J, Balík J, Tlustoš P. Is the long-term application of sewage sludge turning soil into a sink for organic pollutants? Evidence from field studies in the Czech Republic. J Soils Sediments. 2019;19:2445–58. https://doi.org/10.1007/s11368-019-02265-y.
Raaijmakers JM, Paulitz TC, Steinberg C, Alabouvette C, Moënne-Loccoz Y. The rhizosphere: a playground and battlefield for soilborne pathogens and beneficial microorganisms. Plant Soil. 2009;321:341–61. https://doi.org/10.1007/s11104-008-9568-6.
Ridl J, Kolar M, Strejcek M, Strnad H, Stursa P, Paces J, Macek T, Uhlik O. Plants rather than mineral fertilization shape microbial community structure and functional potential in legacy contaminated soil. Front Microbiol. 2016. https://doi.org/10.3389/fmicb.2016.00995.
Papik J, Folkmanova M, Polivkova-Majorova M, Suman J, Uhlik O. The invisible life inside plants: deciphering the riddles of endophytic bacterial diversity. Biotechnol Adv. 2020;44:107614. https://doi.org/10.1016/j.biotechadv.2020.107614.
Sommermann L, Geistlinger J, Wibberg D, Deubel A, Zwanzig J, Babin D, Schlüter A, Schellenberg I. Fungal community profiles in agricultural soils of a long-term field trial under different tillage, fertilization and crop rotation conditions analyzed by high-throughput ITS-amplicon sequencing. PLoS ONE. 2018;13:e0195345. https://doi.org/10.1371/journal.pone.0195345.
McDonald MC, Ahren D, Simpfendorfer S, Milgate A, Solomon PS. The discovery of the virulence gene ToxA in the wheat and barley pathogen bipolaris sorokiniana. Mol Plant Pathol. 2018;19:432–9. https://doi.org/10.1111/mpp.12535.
Gupta PK, Chand R, Vasistha NK, Pandey SP, Kumar U, Mishra VK, Joshi AK. Spot blotch disease of wheat: the current status of research on genetics and breeding. Plant Pathol. 2018;67:508–31. https://doi.org/10.1111/ppa.12781.
Zhu M, Riederer M, Hildebrandt U. Very-long-chain aldehydes induce appressorium formation in ascospores of the wheat powdery mildew fungus Blumeria graminis. Fungal Biol. 2017;121:716–28. https://doi.org/10.1016/j.funbio.2017.05.003.
Dean R, Van Kan JAL, Pretorius ZA, Hammond-Kosack KE, Di Pietro A, Spanu PD, Rudd JJ, Dickman M, Kahmann R, Ellis J, et al. The top 10 fungal pathogens in molecular plant pathology: top 10 fungal pathogens. Mol Plant Pathol. 2012;13:414–30. https://doi.org/10.1111/j.1364-3703.2011.00783.x.
Fones HN, Eyles CJ, Kay W, Cowper J, Gurr SJ. A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by Zymoseptoria tritici. Fungal Genet Biol. 2017;106:51–60. https://doi.org/10.1016/j.fgb.2017.07.002.
Durán P, Jorquera M, Viscardi S, Carrion VJ, de la Mora ML, Pozo MJ. Screening and characterization of potentially suppressive soils against Gaeumannomyces graminis under extensive wheat cropping by Chilean indigenous communities. Front Microbiol. 2017. https://doi.org/10.3389/fmicb.2017.01552.
Acknowledgements
The authors wish to thank Benjamin John Watson-Jones MEng for English language correction.
Funding
We wish to thank the European Regional Development Fund (project No. CZ.02.1.01/0.0/0.0/16_019/0000845) for founding the sample collection and the Czech Science Foundation (project no. 19-02836S) for funding this research, within which the methodology applied for data analysis was developed.
Author information
Authors and Affiliations
Contributions
Conceptualization, KD, HS and OU; methodology, HS, OU and JB; sample preparation, JC; formal analysis, MK and MS; investigation, MK; data curation, MK, PK, and MS; writing—original draft preparation, MK; writing—review and editing, OU and HS; visualization, MK; supervision, HS, JS, JB, and KD; project administration, HS and JS; funding acquisition, KD, JB, JT, and JS. All authors read and approved the final manuscript.
Corresponding authors
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Table S1. Table of soil physicochemical parameters. Fig. S1. Scatter plot of each pair of physicochemical parameters showing Pearson’s correlation coefficients. Table S2. The influence of location, soil type, fertilization and crop rotation on soil chemical properties analysed with Kruskal-Wallis test.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kracmarova, M., Uhlik, O., Strejcek, M. et al. Soil microbial communities following 20 years of fertilization and crop rotation practices in the Czech Republic. Environmental Microbiome 17, 13 (2022). https://doi.org/10.1186/s40793-022-00406-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40793-022-00406-4