Next Article in Journal
Effect of Organic Amendments on the Productivity of Rainfed Lowland Rice in the Kilombero Floodplain of Tanzania
Next Article in Special Issue
Soil Properties after Eight Years of the Use of Strip-Till One-Pass Technology
Previous Article in Journal
Soybean Response to Weather and Climate Conditions in the Krasnodar and Primorye Territories of Russia over the Past Decades
Previous Article in Special Issue
Comparison of the Effects of Different Crop Production Systems on Soil Physico-Chemical Properties and Microbial Activity under Winter Wheat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Changes in Glomalin-Related Soil Proteins (GRSP) Content, Microbial Diversity and Physical Properties Depending on the Type of Soil as the Important Biotic Determinants of Soil Quality

by
Anna Gałązka
1,*,
Jacek Niedźwiecki
2,
Jarosław Grządziel
1 and
Karolina Gawryjołek
1
1
Department of Agriculture Microbiology, Institute of Soil Science and Plant Cultivation—State Research Institute, Czartoryskich 8, 24-100 Pulawy, Poland
2
Department of Soil Science Erosion and Land Conservation, Institute of Soil Science and Plant Cultivation–State Research Institute, Czartoryskich 8, 24-100 Pulawy, Poland
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(9), 1279; https://doi.org/10.3390/agronomy10091279
Submission received: 20 August 2020 / Revised: 24 August 2020 / Accepted: 26 August 2020 / Published: 29 August 2020
(This article belongs to the Special Issue Effects of Agricultural Management on Soil Properties and Health)

Abstract

:
The aim of the study was to evaluate the changes in glomalin-related soil proteins (GRSP) content, microbial diversity and soil physical quality depending on the type of soil measures of soil improvement and changes in soil health. The study was based on a 100-year stationary field microplot experiment where the soil profiles were collected with preserving the natural soil horizons. The microplot experiment was carried out on eight different soil types: Brunic Arenosol (Dystric I), Rendzic Leptosol, Fluvic Cambisol, Haplic Cambisol (Eutric), Gleyic Phaeozem, Brunic Arenosol (Dystric II), Haplic Cambisol (Eutric II) and Haplic Cambisol (Dystric). These soils are the most common types of agricultural soils in Poland. Relatively significant correlations with the soil quality, physical parameters and the glomalin-related soil proteins have been found. The study determined the total GRSP (T-GRSP) and easily extractable GRSP (EE-GRSP) levels in soils as well as the soil physical quality index and soil’s microbial biodiversity. The GRSP depended on the type of soil and correlated with S-Index and also was responsible for the unique chemical and physical properties of soils. Soils characterized by the highest T-GRSP content belonged to the group of very good and good soil physical quality characterized also by high biological activity, for which there were strong correlations with such parameters as dehydrogenase activity (DHA), microbial biomass content (MBC), microbial nitrogen content (MBN) and total bacteria number (B). The highest T-GRSP content and higher microbial diversity were found in Gleyic Phaeozem, Rendzic Leptosol and Fluvic Cambisol. The T-GRSP and EE-GRSP content were additionally correlated with the number of AMF spores. Very poor and poor soil physical quality according to S-Index characterized Brunic Arenosol (Dystric I) and Haplic Cambisol (Dystric). This research indicates that a specific edaphone of soil microorganisms and GRSP content may be of great importance when assessing a soil’s quality and improvements in soil health. The abundance of glomalin-producing fungi significantly affects the quality of the soil. This effect is particularly important for agricultural soils are threatened by ongoing land degradation.

1. Introduction

Organic matter (OM) and the products of its biochemical transformation have a favourable effect on a complex of soil properties that determine its fertility and productivity [1,2]. OM plays a major role in the formation of soil structure and exerts a strong impact on water infiltration, root development and resistance to erosion [3]. Organic matter serves as a store of nutrients and determines cation-exchange capacity. The quality of Polish agricultural soils are rated as the poorest in Europe [3,4]. Production potential of an average hectare of Polish soil equal to 0.6 ha of arable land in the EU. That is why it is so important to properly manage soil resources in the context of the proper functioning of ecosystems and the special protection of high-performance soils [3]. This is mainly due to the type of parent material of soils. Over 70% of Poland’s soils were mainly formed from Pleistocene sands, highly blurred and sorted by glacial waters [4]. Also, some soils were formed from boulder clays. In Poland, such soil types as Brunic Arenosol (Dystric I), Rendzic Leptosol, Fluvic Cambisol, Haplic Cambisol (Eutric), Gleyic Phaeozem, Brunic Arenosol (Dystric II), Haplic Cambisol (Eutric II) and Haplic Cambisol (Dystric) constitute about 75% of the land. These soils are widely used in agriculture for the cultivation of cereals [3,4]. The fundamental assumption of sustainable agriculture is to protect the natural environment and ensure biodiversity in agrocenoses. The abundance of soil microorganisms and higher enzymatic activity are sensitive indicators of the soil biological activity and an appropriate complex of properties determining soil fertility. The well-functioning soil is determined by enhanced biological activity and also a healthy soil also providing a more suitable ecosystem for soil biota. Any changes in soil properties may induce changes in the abundance and activity of soil microorganisms in terms of their species composition and function [2,5].
The quality of soils is connected not only with their microbiological properties but also with its physical and chemical properties [2,6]. Water availability and pH are also important factors affecting the number and activity of microorganisms and ecological relationships between them [7]. In soils with properties favourable for plant production, e.g., texture, soil aggregates, higher OM content and higher pH, intensive development of bacteria is observed [7]. On the other hand, the growth of fungi over bacteria is favoured by drier soil conditions [8]. Soil type is also of great importance for environmental biodiversity. The study of Grządziel et al. [9,10] indicated that differences in bacterial and fungal composition are affected directly by the soil type [5,9,10,11]. Mycorrhizal fungi also play a large role in maintaining proper biocenosis of the soil environment and, at the same time, soil stability [12,13].
Arbuscular mycorrhizal fungi (AMF) are one of the most important biotic determinants of soil environment quality [12,13]. They play a very important role in agrarian ecosystems and, through the production of glomalin related soil proteins, can exert a positive effect on plant growth and health depending on the type of soil [13]. GRSP are glicoproteins that are very characteristic in terms of their structure and physicochemical properties [14]. Glomalin related soil proteins in literature is reportedly a stable and persistent protein produced in copious quantities by mycorrhizal fungi and may be an important pool of organic N in soil [13,14]. GRSP is only operationally defined by its extraction method as the total glomalin related soil proteins (T-GRSP) and easily extracted glomalin related soil proteins (EE-GRSP) [12]. Gillespie et al. [14] conducted research aimed at evaluation of glomalin structure and their characterization. On the basis of modern research methods such as synchrotron-based X-ray absorption near-edge structure (XANES) spectroscopy and pyrolysis field-ionization mass spectrometry (Py-FIMS) authors revealed that GRSP contains a consortium of proteins along with many impurities [14]. On the basis of employing proteomic techniques, Gillespie et al. [14] found that glomalin related soil proteins itself may be a thioredoxin-containing chaperone. The authors did not find no homologies with proteins or DNA of mycorrhizal origin were detected. In addition on the basis of proteomics techniques the authors have proved that GRSP contains large amounts of soil-related heat-stable proteins and proteins of non-mycorrhizal origin [14]. Therefore, based on the Gillespie et al. research [14], it should be considered that the extraction procedure that defines GRSP as a mixture of compounds. Hence, due to the very complicated nature of this substance the chemical nature of glomalin has yet to be conclusively determined. However, it is extremely interesting how GRSP itself can affect soil quality and its physicochemical and biological properties.
GRSP (T-GRSP and EE-GRSP) are stable, water-insoluble and degradation-resistant molecules stabilizing soil aggregates and protecting them from disintegration [15,16]. The hydrophobic domains of the molecules are responsible for the difficulty in their extraction and their low solubility. On the one hand, such properties make GRSP very stable compounds, ideal for protecting soil aggregates from degradation. On the other hand, they are difficult to identify and establish the exact structure of these molecules [16].
Similarly, the physicochemical properties of soils may have a significant impact on the glomalin content in soil [12]. The brief literature review presented above clearly indicates that GRSP not only play an important role in the functioning of the symbiotic system between plant roots and AMF but they are also a very important fraction of soil OM determining the structural properties of soils and improving soil stability and quality [14].
Soil quality is concerned with physical, chemical and biological aspects and defines the capacity of agricultural soil to deliver multiple functions (productivity, nutrient cycling, water filtration, carbon sequestration and microbial diversity) [16,17,18]. This capacity can be depleted, resulting in degradation of the system by compaction, erosion or loss of soil OM. An index as described by Dexter [16] using S-theory was identified as a key metric of overall soil quality. Based upon the S-theory, the S-term offers a singular value which is considered reflective of overall soil physical quality, and the S-term can be considered suggestive of not only physical quality but also both chemical and biological quality. The S-theory proposes that soil physical properties and behaviour are essentially controlled by soil structure, which is expressed as pore size distribution. Thus, the S-Index offers a simple scale that has the same physical meaning regardless of soil type and can therefore be used to compare soil physical quality across soil types and spatial scales [18].
The aim of the study was to evaluate the changes in glomalin-related soil proteins (GRSP) content, microbial diversity and soil physical quality depending on the type of soil as the important determinants of soil quality. The connections and relations between the soil quality, physical parameters including S Index and the glomalin-related soil proteins in context to maintenance of soil health and quality will be looking for.

2. Materials and Methods

2.1. Field Experiment and Soil Samples

The study was based on a long-term stationary field microplot experiment established in 1880–1881 [19]. The experiment is located in areas belonging to the Institute of Soil Science and Plant Cultivation—State Research Institute in Pulawy, Poland (geographical coordinates of the experiment 51°24′52.9″ N 21°57′34.5″ E). Plots of 12 m2 and 1 m deep were separated with concrete walls and filled with soil. The soil was placed in microplots in accordance with the natural profile systems. The soil profiles came from the surrounding area of Pulawy, Poland [19]. Soil horizons were placed in accordance their natural arrangement, to specially prepared concreted micro-parcels. Microplots do not have a concrete bottom in order to ensure contact of the soil with the natural substrate. The soil profiles were collected without destroying the soil structure, preserving the natural soil horizons. From the beginning of the experiment, all soils were treated equally including the same cultivated plant, agriculture management (hand tillage; hand-picked weeds, digging soil with a hoe), irrigation and fertilization. The number of microplots were presented in Figure 1A. An example of a micro plot was presented in Figure 1B.
Potassium and phosphorus fertilization was applied before sowing in the amount of 100 kg K2O and 80 kg P2O5/ha. Nitrogen fertilization 90 kg N/ha was used on two dates: 50 kg N/ha before sowing and 40 kg N/ha during the time of firing. Soil samples were taken every year in April when there were no plants. In this publication, the authors decided to analyze the soil itself without plants. The same soil samples had previously been tested in the context of the biodiversity measure [8,10].
Samples of bulk soil were collected in April 2016–2018, annually as three biological replicates from each plot at a depth of 0–15 cm (a total of 15 sampling occasions with an each microplot; five for each biological replicates). Figure 2 shows an exemplary scheme of sampling from the M1 plot. Subsequently, soil samples were, sieved through a 2 mm sieve, and stored at −20 °C prior to DNA extraction and at 4 °C until analysis of biological activity. Moreover, the soils always were under the same weather conditions and hand tillage management.
The microplot experiment represents the most typical and common soil types of agricultural soils found in Poland [6]. Three soils were characterized by high organic carbon, pH 7.3–7.5 and good quality for plant cultivation: Rendzic Leptosol, Gleyic Phaeozem and Fluvic Cambisol. Another two soils, Haplic Cambisol (Eutric II) and Haplic Cambisol (Dystric), exhibited low organic carbon, pH 5.4–5.6 and relatively good conditions for cultivation of certain plants. The last three soils were acidic (pH 4.0–4.7), characterized by low organic carbon, and poor quality for cultivation: Brunic Arenosol (Dystric I), Brunic Arenosol (Dystric II) and Haplic Cambisol (Eutric I) (Table 1). The full description of the experiment along with engravings can be found in the supplement for the paper Grządziel and Gałązka [11]; https://www.frontiersin.org/articles/10.3389/fmicb.2019.00006/full#supplementary-material. The main advantage of this microplot experiment is that all soils are collected in the very small area, and are subjected to the same conditions (geographical location, weather, the same climate etc.). Since 1979, the plots have been planted mostly with cereals as the main crop with mineral fertilizers, and mainly mustard, phacelia or leguminous plants as the second crop. In 1984, the plots were fertilized with compost (80 t/ha) under potatoes [8]. In 2016–2018, the same plant species (winter wheat) was always grown in all plots at the same time and the soil received the same NPK fertilization.
The microplot experiment has been described in detail in the papers by Czaban et al. [8], Grządziel and Gałązka [10,11], Grządziel et al. [9] and Siebielec et al. [20].
Soil types were determined in accordance with the World Reference Base (WRB), whereas soil texture classes were determined on the basis of the United States Department of Agriculture (USDA) Soil Texture Classes System (Table 1). Based on the S Index, these soils can be divided into a few groups: very poor soil physical quality (S ≥ 0.025, Brunic Arenosol Distric I), poor soil physical quality (S ≤ 0.02, Brunic Arenosol Distric II), good soil physical quality [0.050 > S ≥ 0.035, Fluvic Cambisol, Rendzic Leptosol, Haplic Cambisol (Dystric), Haplic Cambisol (Eutric II), Haplic Cambisol (Eutric I)] and very good soil physical quality (S ≥ 0.050, Gleyic Phaeozem). It was closely related to the total glomalin related soil protein content. This division on the basis of the S index value was closely related to the GRSP content (Table 1).
The meteorological conditions during the growing seasons (2016–2018) was presented in Table 2. The average annual April temperature reaches +8.7 °C. The annual April rainfall is 42.9 mm. The snow cover lasts 80–85 days and the vegetation period lasts 210–220 days a year.

2.2. Clay Dispersion Measurements

The method used for determination of readily dispersible clay (RDC) in normalized turbidity units (NTU/(g/L)), is that described by Czyż and Dexter [21]. The amounts of total clay (TC) and RDC were determined by a turbidimetric method using a standard procedure. The stability of soils was calculated as 1-(RDC/TC) NTU/(g/L).

2.3. Soil Water Retention

Soil water retention curves were measured with undisturbed soil samples using standard sand boxes (Eijkelkamp, Giesbeek, The Netherlands) for low levels of water potential (10, 20, 40 and 100 hPa) and pressure plates (Soil Moisture Inc., Santa Barbara, CA, USA) for high levels of water potential (250, 500, 1000, 2000, 4000, 8000 and 15,000 hPa). The mean values of the water content corresponding to each level of water potential were fitted to the van Genuchten equation [22].

2.4. Soil Physical Quality S-Index

Soil physical quality S-Index, as the slope the of the water retention curve at its inflection point, was calculated based on the equation described by Dexter [16] and Dexter and Czyż [17]. S-Index is a measure of soil microstructure which controls many of the soil physical properties and it is hoped that it may be useful for the assessment of soil quality. Based on experience in a wide range of soil properties and on experimental data of various soils, Dexter presented categories of soil physical quality ‘S’ [16]. The S-Index is derived from the relationship between the gravimetric soil water content and the natural log of matric tension. This is calculated as the slope of the soil water retention curve (SWRC).

2.5. Microbial Communities and Enzymatic Analysis

Microbiological counts were expressed as the number of colony-forming units (CFU) per gram of dry mass of soil (d. m.). The total number of bacteria (B) was determined by the dilution method on agarized soil extract. The total number of fungi was determined on Martin’s medium [23]. Ammonifying bacteria (AM) were determined according to Rodina’s method [24] and Azotobacter spp. according to Fenglerowa’s method [25]. Microbial biomass was determined by the chloroform fumigation-extraction method according to Ghani et al. [26]. Enzymatic activity was determined spectrophotometrically: soil dehydrogenase activity using the TTC method (Polish Standard. PN-EN ISO 23753-1) [27], phosphatases activity by the p-NPP method [28].

2.6. Glomalin Related Soil Proteins Content Analysis

Glomalin related soil proteins content was determined as proposed by Wright et al. [29] with modifications described by Gałązka et al. [7]. GRSP concentrations were determined spectrophotometrically as both T-GRSP and EE-GRSP [24]. TG-GRSP was obtained by repeated extraction from 1 g of ground dry-sieved soil with 8 mL of 50 mM citrate, pH 8.0 at 121 °C for 60 min. The EE-GRSP protein was extracted from 1 g of ground dry-sieved soil with 8 mL of 20 mM citrate, pH 7.0 at 121 °C for 30 min. After each autoclaving cycle, the supernatant was removed by centrifugation at 4000 rpm for 15 min and stored. Extracts from each cycle were pooled, centrifuged at 9000 rpm for 5 min to remove soil particles and then analysed. Extracts from each replicate were pooled and analysed, and the protein in the supernatant was determined by the Bradford dye-binding assay with bovine serum albumin as the standard on a 96-plate reader (Victor Perkin Elmer, Waltham, MA, USA).

2.7. Determination of Number of AMF Spores in Soils

The total number of AMF spores was evaluated according to the Phillips and Hayman method [30]. Root mycorrhizal structures were stained with 0.05% (w/v) trypan blue. The number of AMF spores was determined in averaged samples of 100 g of air-dried soil (2 × 50 g), while AMF spores separated on filters with mesh sizes of 150, 100 and 50 μm were counted, which, when added together, give the number of spores in a soil sample.

2.8. Total DNA Extraction from Soil

Fresh soil samples were weighed; similar amounts (300–350 mg) were collected into 1.5 mL tubes and used to extract DNA with a FastDNA™ SPIN Kit for Soil (MP Biomedical) according to the manufacturer’s instructions. The concentration and purity of extracted DNA were measured in a NanoDrop™ 1000 Spectrophotometer (Thermo Scientific, Waltham, MA, USA). As a result, approximately 60 ng/µL DNA was obtained from each sample which corresponds to 17.5 µg of DNA from 1 g of soil sample, with an A260/280 ratio of 1.4–1.6. For further experiments, each DNA sample was diluted to a final concentration of about 10 ng/µL.

2.9. PCR Amplification

Fragments of 16S rRNA gene were amplified to perform genetic differentiation of the bacteria found in soils. To increase the specificity of the reaction, a two-step PCR method was performed. In the first step, a 1492 bp fragment of the 16S rRNA gene was amplified with primers 27F and 1492R. The PCR mixture (total 20 µL volume) contained 4 µL Silver Hot Start PCR MIX LOAD (Syngen Biotech, Poland), 0.15 µM of each primer, 20 ng of template DNA, and ultrapure water. Each primer set can be found in Table 3. The amplification products were used as templates for PCR-DGGE with two sets of primers. GC-338F and 518R were used to amplify the hypervariable V3 region of about 197 bp present in most of the Bacteria. Amplification of a 280 bp product of Actinobacteria was performed with Act235F and Eub518R primers. The PCR mixture (total volume 50 µL) contained 25 µL of DreamTaq Green PCR MasterMix (Thermo Scientific), 0.4 µM of each primer, 2 µL of diluted 16S rDNA PCR product, and ultrapure water. Detailed PCR parameters are presented in Table 4. Amplicons were confirmed by electrophoresis on 1% agarose gel stained with ethidium bromide in 1× TAE buffer using a GeneRuler 100 bp DNA Ladder (Thermo Scientific).

2.10. Denaturing Gradient Gel Electrophoresis (DGGE)

Amplicons were separated using the DCode™ Universal Mutation Detection System (Bio-Rad Laboratories, Hercules, CA, USA) with modifications to the method described by Muyzer et al. [32]. Equal quantities of DNA (650 ng) were loaded into 8% (w/v) PLA gels (40% acrylamide/bis solution, 37.5 : 1) in 1× TAE buffer. The denaturing gradient for Bacteria-V3 amplicons ranged from 40% to 60% and for Actinobacteria from 50% to 70%, assuming that 100% is 7 M urea and 40% (v/v) deionized formamide. Electrophoresis was performed at 55 V for 16 h (Table 3). A constant temperature of 60 °C was maintained. Gels were stained for 30 min using 10,000 times-diluted SYBR® Gold Nucleic Acid Gel Stain (Thermo Scientific) and visualized on a UV transilluminator. Images were then analysed using Quantity One 4.6.9 software (Bio-Rad).

2.11. Community-Level Physiological Profiling (CLPP) Analysis Using Biolog EcoPlates

The CLPP was evaluated using Biolog EcoPlates (Biolog Inc., Hayward, CA, USA) with 31 different carbon sources [34]. The soil suspension for the inoculation of wells in microplates was prepared as follows. One gram of soil was weighed and transferred to a conical flask holding 99 cm3 sterile 0.9% NaCl, vortexed for 30 min at 150 rpm and at 25 °C, after which the samples were cooled for 30 min to 4 °C. After that, 120 mm3 was transferred to each well of an EcoPlate and incubated in the dark at 25 °C for 120 h. The experiment included three replications. The results were read on a MicroStation ID system (Biolog Inc., Hayward, CA, USA). The extent to which carbon sources were used was determined through the reduction of colourless tetrazolium chloride (TTC) to red triphenyl-formasane (TPF) [34,35]. The intensity of colour development was recorded at λ = 590 nm for a period of 120 h at 24 h intervals. The activities of soil microorganisms are based on all carbon sources and on grouped sources defined as amines and amides, amino acids, carbohydrate, carboxylic acids and polymers. The results were expressed as average well-colour development (AWCD) and Shannon–Wiener (H’) index.

2.12. Statistical Analysis

Statistical analyses were performed using the packet STATISTICA.PL (10) (StatSoft. Inc., Tulsa, OK, USA). The variables had a normal distribution. Collected data were subjected to analysis of variance (ANOVA) for comparison of means. Significant differences were calculated according to Tukey’s post hoc test at a significance level of p < 0.05. Cluster analysis, including grouping of the treatments and features, was performed on standardized data (95 carbon sources in triplicate) from the average absorbance values at 120 h (Biolog EcoPlate). The dendrogram was prepared with scaled bond distances on the axis (Ward’s method), and the boundary was marked according to Sneath’s criteria. AWCD was evaluated according to the formula AWCD = Σ(C − R)/95, where C is the absorbency in each well and R is the absorbency in the control well. The Shannon–Wiener (H’) index was evaluated in accordance with the formula H’ = −Σpi (lnpi), where pi is the ratio of the absorbance of each well to the absorbency of all wells. The results were also submitted to principal component analysis (PCA). All soil quality parameters tested in this study were used as variables in the PCA analysis. The PCA analysis was performed using the packet STATISTICA.PL (10), the data were standardized.

3. Results

3.1. Physical and Chemical Properties of Soil

The physical and chemical properties of soils are presented in Table 1. The carbon content of soil samples ranged from 0.522% (Brunic Arenosol Dystric II) to 3.101% (Gleyic Phaeozem). The nitrogen content in soil samples ranged from 0.046% (Brunic Arenosol Dystric II) to 0.225% (Gleyic Phaeozem). The RDC/TC ratio was highest in Brunic Arenosol Dystric I (0.313) and lowest in Fluvic Cambisol (0.054).
Based on the parameters of van Genuchten, the S-Index was calculated (Table 5). Gleyic Phaeozem was characterized by the highest value of the S-Index (0.076)—very good SPQ (Soil Physical Quality). A higher S-Index value also corresponded to a higher glomalin content (Table 1). The highest GRSP content was found in Gleyic Phaeozem (6.33 mg/g d.m. soil), while the lowest was found in Brunic Arenosol (Dystric I)—2.87 mg/g d.m. of soil; this value corresponded to the lowest S-Index for this soil (0.018)—very poor SPQ.

3.2. GRSP Content

Based on an ANOVA analysis, it was found that both the type of soil and the year of sampling were significant in assessing the soil’s biological activity (Table 6). As variables in the analysis, the following activities were chosen: T-GRSP and EE-GRSP content, enzymatic activity, total number of bacteria, total number of fungi and indicators of soil functional metabolism evaluation on the example of Biolog EcoPlate analysis. Detailed information is presented in Table 6.
The concentrations of T-GRSP and EE-GRSP depended primarily on the type of soil (Table 7). The highest concentrations of T-GRSP were found in Gleyic Phaeozem which was also characterized by a higher organic carbon (Corg) content (3.1%). In addition, a high glomalin content was found in the following soils: Rendzic Leptosol, Fluvic Cambisol and Haplic Cambisol (Eutric II) (4.43, 4.59 and 4.68 mg/g d.m. of soil, respectively) (Table 5). Soils with a low T-GRSP content were: Brunic Arenosol (Dystric I), Brunic Arenosol (Dystric II), and Haplic Cambisol (Dystric) (2.87, 3.07 and 3.79 mg/g d.m. of soil, respectively). The T-GRSP and EE-GRSP content was also strongly associated with the type of soil. The highest EEG-GRSP content was found in Haplic Cambisol (Eurtic I), Gleyic Phaeozem and Haplic Cambisol (Eutric II), (2.31, 2.38 and 2.44 mg/g d.m. of soil, respectively) (Table 7).
The highest number of AMF spores was found in Haplic Cambisol (Eutric II), Haplic Cambisol (Dystric) and Fluvic Cambisol (Table 8). A different AMF spore content was also found depending on the diameter of the spores. The total number of 50 µm diameter spores was highest in Haplic Cambisol (Dystric), but the number of 150 µm diameter spores was highest in Haplic Cambisol (Eutric II). The lowest number of AMF spores was found in Brunic Arenosol (Distric I).

3.3. Microbial Diversity

The microbiological diversity differed depending on the type of soil. As a result of PCR-DGGE analysis, two separate fingerprint images were obtained (Figure 3A,B).
Despite the high background signal generated by the large amount of DNA applied to a gel, there can be seen clearly distinct bands, which are characteristic for certain soils. Although both sets of primers generated similar results, dividing samples into two (Figure 4A) or three (Figure 4B) main groups, based on the results from the clustering, it can be observed that the use of Actinobacteria-specific primers resulted in more discriminating results. Each of these methods put together good-quality soils, while acidic, poor-quality soils were grouped into a separate cluster, providing the background for subsequent analyses.
The results of our study show that not only are biochemical differences detectable in different soil types, but also the soil microbiome differs significantly among microplots. The result is not only a wide genetic diversity between samples using universal primers, but also considerable variation within the phylum Actinobacteria. These observations are essential, notably proving the soil buffering capacity. In such a complex biological material as soil, it is important to adjust the amplicon quantity loaded into a gel, at the same time remembering that the more PCR products there are, the more the band smearing effect is observed. This issue is caused by many factors, for example interactions between molecules, formation of single-stranded DNA during PCR, and electrophoresis [36].

3.4. Biological Activity of Soil

The basic enzymatic activity (dehydrogenase activity, acid and alkaline phosphatase activities) are presented on Figure 5. The enzymatic activity was strictly dependent on the type of soil and year. The highest dehydrogenase activity was found in Gleyic Phaeozem and the lowest in Brunic Arenosol (Dystric I and II). The highest alkaline phosphatase activity was observed in Gleyic Phaeozem. Very slight changes of enzymatic activities were observed also in years. In general, soils maintained enzymatic activity in the years 2016–2018 at a similar level, which indicates a stabilized microbiome of soils.
The highest content of microbial biomass of carbon (MBC) was observed in Gleyic Phaeozem and the lowest in Haplic Cambisol (Eurtic II) (Figure 6). The highest content of microbial biomass of nitrogen (MBN) was observed in Gleyic Phaeozem, but equally was also high in Rendzic Laptosol and Fluvic Cambisol.
The total number of bacteria and fungi also varied depending on the type of soil and year (Figure 7). In the case of the total number of bacteria and total number of fungi, much greater variability was observed in years.
Table 9 presents a summary of means from 2016–2018 as changes in soil biological activity. Enzyme activity and C and N biomass as well as the total number of bacteria and fungi are discussed above. On the other hand, in the case of the total number of ammonification bacteria, the highest number of these bacteria was found in soil samples Rendzic Leptosol, Gleyic Phaeozem and Haplic Cambisol (Dystric), (Table 8). A higher total number of Azotobacter spp. was observed only in soil Rendzic Leptosol, Fluvic Cambisol and Gleyic Phaeozem (Table 9). The presence of Azotobacter bacteria depended on the pH of soil. Azotobacter spp. prefers neutral soils (pH 7.0–7.5).
Based on the parameters of microbiological activity, soil types were grouped on two main groups. Based on Ward’s cluster analysis, two main groups were distinguished (Figure 8). The first group included soil types with high biological activity: Gleyic Phaeozem, Fluvic Cambisol and Rendzic Leptosol, while the second group included soils with lower soil biological activity: Haplic Cambisol (Distric), Haplic Cambisol (Eutric II), Haplic Cambisol (Eutric I), Brunic Arenosol (Distric I) and Brunic Arenosol (Distric II). In addition, the first group distinguished two subgroups, the first containing Gleyic Phaeozem and the second Fluvic Cambisol and Rendzic Leptosol. In turn, in the second group, three subgroups were distinguished, the first containing Haplic Cambisol (Distric), the second Haplic Cambisol (Eutric II), Haplic Cambisol (Eutric I) and the third), Brunic Arenosol (Distric I) and Brunic Arenosol (Distric II), (Figure 8).
The Pearson’s correlation coefficient of selected microbiological, chemical and biochemical parameters and GRSP content are presented in Table 10. Statistically significant correlations were found between the TG-GRSP content and the enzymatic activity of soils and the total number of bacteria as well as the content of carbon (Corg) and nitrogen in the biomass of microorganisms (Table 10). The S-Index value was closely correlated with the glomalin content (T-GRSP and EE-GRSP), dehydrogenase activity (DHA) and the total number of bacteria. The glomalin content was additionally correlated with the number of AMF spores with a diameter of 150 µm.
The correlation coefficients for selected parameters of soil biological activity and GRSP content between the two principal components PC1 and PC2 in the PCA are presented in Table 9. The principal components of the PCA analysis showed strong correlations between the parameters of soil types and microbial indicators. Selected indicators of soil microbial diversity explained 71.73% of biological variability in soils (Table 11).
The bi-plot PCA of selected microbiological and biochemical parameters, glomalin content, biodiversity indices and Biolog EcoPlate data from soil samples incubated for 120 h was presented in Figure 9. The Gleyic Phaeozem, Fluvic Cambisol and Rendzic Leptosol were characterized by high biological activity, for which there were strong correlations with such parameters as DHA, MBC, MBN, S-Index and total bacteria number (Figure 9). The Haplic Cambisol (Eutric I), Haplic Cambisol (Eutric II) and Haplic Cambisol (Distric) were characterized by biological activity and there were strong correlations with total number of fungi and acid phosphatase. The Brunic Arenosol (Distric I) and Brunic Arenosol (Distric II) were characterized by very low biological activity of soil (Figure 9).

4. Discussion

Investigations carried out in different countries have shown that GRSP is commonly recorded in soils [7,8]. It has been detected in agricultural soils, forest and meadow soils, and agricultural wasteland [13]. The highest concentration of GRSP was recorded in tropical forest soils (over 100 mg/g soil), whereas desert soils were characterized by the lowest glomalin concentrations (below 0.5 mg/g soil) [36,37]. In Poland, there have been no studies on the content of this important fraction of OM in different types of soils so far. Therefore, the aim of the study was to evaluate the changes in glomalin-related soil proteins (GRSP) content, microbial diversity and soil physical quality depending on the type of soil measures of soil improvement. Thanks to this, it will be possible to better understand how some microbiological and physicochemical properties affect the physical quality of soils, defined by the S index.
The highest T-GRSP and EE-GRSP content and higher microbial diversity were found in Gleyic Phaeozem and Rendzic Leptosol. The soils characterized by a low T-GRSP content were Brunic Arenosol (Dystric II). Soils with the lowest T-GRSP and EE-GRSP content were also characterized by the lowest soil stability measured as the proportion between TC content and RDC content. The glomalin related soil proteins content (T-GRSP and EE-GRSP) depended on the type of soil and correlated with S-Index. According to the S index, soils with very poor and poor physical quality were characterized by the lowest content of T-GRSP, as well as with a low content of organic matter, very acidic reaction and sandy. In turn, soils characterized by the highest T-GRSP content belonged to the group of very good and good soil physical quality. Soils with a high T-GRSP content were more stable because the clay contained in them is not easily dispersible. This is all the more important because all the tested soil types are used in Poland for agriculture [3,4]. This indicates the need to care for the quality of the soil, in particular by maintaining or increasing the content of organic matter and maintaining the optimal soil pH for crops. Hence increas the soil quality and health will be crucial for improving plant growth and yield.
Small-scale environmental experiments are often linked with some concerns. One of them is that if microplots are located close to each other, they will eventually mix and unify the composition of the soil microbiota. These results demonstrate, however, that even if the effect occurred in the last 130 years, each soil retains its individuality in terms of microbiological composition. It also suggests that this negative effect may be only superficial, and a depth of 0–15 cm seems to be more or less resistant to such fluctuations. Soils from the original locations for experiment were placed in micro plots while maintaining the natural system of soil horizons and have undergone over 130 years of unchanged cereal production.
GRSP plays a significant role mainly in the formation of the soil aggregate stability, as it facilitates adhesion of mineral particles [14,38]. This protein covers the surface of soil particles, creating a characteristic, protective coating. Hence, due to the production of glomalin, AMF have aroused growing interest as a potential factor supporting plant growth and health and improving the soil structure [15]. Increasing soil inoculation by AMF can have an important impact on plant yield and soil health.
The role of AMF in soil aggregation and, hence, in enhanced water infiltration, protection of OM from decomposition, and enhanced soil resistance to erosion seems equally important or perhaps even more important than their role in nutrient uptake by plants. Therefore, control of the growth of mycorrhizal fungi and glomalin production may be useful for monitoring desertification and degradation of soil. Unfavourable conditions for the growth of endomycorrhizal fungi, i.e., soil contamination, climatic change or cultivation techniques, may be the reason for increased glomalin production aimed at protection of the mycelium [38].
Similarly, the physicochemical properties of soils may have a significant impact on the glomalin content in soil [13]. Soil microorganisms are the most important and at the same time the most diverse group of organisms that are part of the Earth’s biosphere, creating multi-species communities and a network of relationships between individual physiological groups [39,40]. Therefore, the processes of soil formation and soil degradation carried out by communities of microorganisms should be seen as the sum of functions for which microbial communities are responsible, and not only individual species. Studies on the activity of microorganisms in communities are necessary to learn about the ecology of microorganisms in biocenoses and should be analysed in relation to existing environmental conditions and time [40].
In the paper by Czaban et al. [8], (research conducted on the same plots) the authors showed that soil microorganisms differing in their water requirements are associated with different soil particulars. The total number of bacteria is closely related to small particles and small pores usually filled with water, but non-xerophilic fungi are associated with medium particles and medium pores which are usually the region of near the air–water interface. Moreover, the number of colony forming units (CFU) of xerotolerant/xerophilic fungi and particularly formulated indices of xerotolerance of soil fungal communities (in contrast to CFU of bacteria) are negatively correlated with water content, soil OM content, soil pH and an index of soil quality [5]. Our results also show that the microbiological activity of soils and the GRSP content (T-GRSP and EE-GRSP) are closely related to the physical quality of soils and the S-Index. A strong correlation between S-Index and T-GRSP was demonstrated (Pearson’s correlation coefficient = 0.867). In addition, a close relationship between TG-GRSP and organic carbon (Pearson’s correlation coefficient = 0.879) and nitrogen content (Pearson’s correlation coefficient = 0.937) was demonstrated. It is very important in the context of preserving the original physical properties of these soils. Siebielec et al. [20] by conducting research on the same field experiment, determined that for over 20 years, soils have maintained at a similar level parameters such as pH, organic carbon content, the amount of bioavailable potassium, phosphorus or magnesium. Our research shows that soils of good physical quality are characterized by a high content of T-GRSP, EE-GRSP [41,42].
Also, a negative impact of RDC content on physical soil quality was presented by Gate et al. [43] and Dexter and Czyż [17]. Czyż and Dexter [21] showed that clay dispersion, measured by the turbidity of the suspension, is a sensitive indicator of soil stability, which is important for protecting the soil environment. Basically, the soil is of poor physical quality when it shows low resistance to degradation caused by a high content of RDC. Which may be the reason for the low water infiltration rate, increased surface runoff, poor coherence, low aeration and root density as well as mechanization difficulties [16,44,45].
In this study a significant number of AMF spores were found in soils. AMF spores may remain after harvesting plants from the previous year. And despite the fact that soil samples were taken before sowing (the plants were not yet in the field), so AMF spores were found in the soil. In our previous paper on the evaluation of the structural diversity of fungi on the same plots, only a small presence of AMF was found (<1%, using next generation sequencing (NGS) [11]. Our results shows that the number of AMF spores was from 2.85 to 15.32 per 1 g of soil, i.e., 100 g is ~300 to 500 pieces of AMF spores. Calculated for the amount used for DNA isolation (about 350 mg of soil), in our previous paper with NGS, there are about 2 spores in the soil sample. Therefore, with such a small amount taken for testing in the sequencing of the ITS variable region for fungi in the NGS method, it could happen that this material has not been amplified. Similar results in AMF spores were obtained by authors of other papers [46,47,48]. Therefore, the results obtained by us are all the more valuable because both a modern research method (NGS in paper Grządziel and Gałązka [11]) and traditional methods of counting AMF spores were used in the present study. These methods will complement each other very well.
The amount of OM in soils is the basic indicator of soil quality assessment in terms of their physicochemical properties, such as sorption and buffer abilities, as well as their biological transformation processes, important from the point of view of habitat functioning and referred to as biological activity. A high organic carbon content in soils is a factor stabilizing their structure, reducing susceptibility to compaction and degradation due to water and wind erosion, and positively affecting the physical quality of soils [1]. Our results also confirm this, indicating that the highest physical quality of soils was characterized by the soil with the highest organic matter content (Gleyic Phaeozem). This is all the more important as organic matter is a fundamental indicator of soil health. In addition, soils rich in organic matter are also of great importance in carbon sequestration [44,45].
The European Parliament resolution (EP Resolution 2007) lists among the main threats to the proper functioning of soils a decrease in the content of OM, because its loss can cause many adverse phenomena associated with the disturbance of chemical, physical and biological processes that directly affect the quality of soil and its fertility [46,47,48,49,50]. Among the 8 types of analyzed soils, Gleyic Phaeozem, as a soil rich in organic matter, is the most sensitive to. Their share in the area of mineral soils in Poland does not exceed 2% [3]. That is why it is so important to properly manage these soils so as not to lose organic matter in them. The preservation of soil organic matter is important not only for maintaining the productive functions of soils, but also for the role of soils in the sequestration (binding) of carbon dioxide from the atmosphere, contributing to the reduction of the greenhouse effect [4]. The soil with the highest agricultural productivity in Poland constitutes only about 26% of all arable land, which covers about 60% of the country’s area. This group also includes Gleyic Phaeozems, which in Poland are predominantly over 90% agricultural land. Haplic Cambisols—developed from loess, some Fluvisols and Rendzic Leptosols are also highly productive. Hence the greatest care to maintain these soils of the best quality. These are the most productive soils, where the highest yields are obtained. Moreover, these are soils characterized by high water retention and high resistance to chemical degradation. For these reasons, in Poland, first of all, soils with the highest productivity should be protected against non-agricultural use and against degradation [3]. The remaining soil types are low productivity soils and account for more than 70% of the agricultural land area. The majority of them are sandy soils with low water retention and low organic carbon content [3,4].
Due to the diversity and activity of soil microorganisms, different types of soil aggregates are formed, in which the cells are bacterial and fungal hyphae, as well as polysaccharides synthesized by microorganisms, ensuring stability and giving appropriate structure to the soil [51,52,53]. The benefits of improving soil structure directly related to maintaining biodiversity, especially in agricultural environments, include: reduced risk of erosion, better root penetration and better access to soil moisture and nutrients, easier germination due to reduced surface crusting, and greater infiltration, retention and availability water due to improved porosity.
In the paper by Sarapatka et al. [53], the authors showed the importance of glomalin in soil degradation. The authors studied the sensitivity and changes of glomalin in the soil environment as a result of soil degradation, how it correlates with other soil characteristics, and how useful it is for indicating changes in soil properties. Statistically significant differences were recorded between erosional and depositional sites, the latter showing higher concentrations of both glomalin and characteristics relating to soil OM and content of nutrients.
Therefore, it is extremely important to search for and develop effective low-cost environmentally friendly solutions and techniques that can be effective tools for improving the quality of the soil environment, in particular in light of increasing the OM content in soil sustainable agriculture [54]. Undoubtedly, increased presence of glomalin-producing mycorrhizal fungi in the soil can be one of such solutions. Furthermore, investigations of the impact of soil type on the microbiological quality of the soil environment indicated by selected parameters of soil biological quality allowed a profound analysis of not only quantitative but also qualitative changes occurring in soil OM and microbial communities. They also allowed assessment of these relationships in a broader aspect of soil productivity protection. The obtained research results may be useful in planning the soil sustainable management because they indicate the role of microbiological properties on the physical quality of soils, which is one of the elements determining soil health. However, study by Lobry de Bruyn and Andrews [55], revealed there may be some problems when it comes to using soil health research as a decision-making tool, and all of these guide best practice. To make and maintain soil sustainable management, farmers want to be confident that the benefits to their soil are real. Soil assessment and monitoring can benefit us all stakeholders, including farmers if undertaken with clear guidelines and within supportive structures.
Undoubtedly, increased presence of glomalin-producing mycorrhizal fungi in the soil can be one of such solutions. Furthermore, investigations of the impact of soil type on the microbiological quality of the soil environment indicated by selected parameters of soil biological quality allowed a profound analysis of not only quantitative but also qualitative changes occurring in soil OM and microbial communities. They also allowed assessment of these relationships in a broader aspect of soil productivity protection.

5. Conclusions

  • The content of glomalin related soil proteins depended not only on the soil type but also on the biological and physical activities and correlated with S-Index and was the important biotic determinant of soil quality.
  • Significant correlations were found between the total carbon content (%C) and the concentration of T-GRSP and EE-GRSP. Soils characterized by highest GRSP content belonged to the group of very good and good SPQ characterized also by high biological activity.
  • Soil biological activity depended strictly on the type of soil and remained stable over the years examined. The Gleyic Phaeozem, Fluvic Cambisol and Rendzic Leptosol were characterized by high biological activity, for which there were strong correlations with such parameters as DHA, MBC, MBN, S-Index and total bacteria number. The Haplic Cambisol (Eutric I), Haplic Cambisol (Eutric II) and Haplic Cambisol (Distric) were characterized by moderate biological activity and there were strong correlations with total number of fungi and acid phosphatase. The Brunic Arenosol (Distric I) and Brunic Arenosol (Distric II) were characterized by very low biological activity of soil.
  • This research indicates that a specific edaphone of soil microorganisms may be of great importance when assessing potential soil resistance to degradation and improving quality and health. Our research indicates a significant role of soil organic carbon content on soil quality, which may also have important implications and connections to carbon sequestration, preventing the negative effects of climate change and improving soil health and biodiversity of soil microorganisms.

Author Contributions

Conceptualization, A.G.; methodology, A.G., J.N., J.G. and K.G.; software, A.G. and J.G.; investigation, A.G. and J.N.; writing—original draft preparation, A.G., J.N. and J.G.; visualization, A.G., J.N. and J.G.; supervision, A.G.; funding acquisition, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research was conducted within the statutory activity of IUNG-PIB (1.20) and the frames of Task 1.2 “Agricultural soils evaluation taking into account proper functioning of soil ecosystems and indication of actions preventing degradation processes” and 1.4. “Evaluation and formation of biodiversity of soil and microbial activity of soil with regard to habitat conditions and management system”, Multi—Annual Programme IUNG—PIB (2016–2020).

Acknowledgments

The authors thank technical staff Jarosław Ciepiel, Olga Ciotucha and Andrzej Księżniak for their help in laboratory experiments. The authors thank the Department of Cereal Crop Production for providing microplots for research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ukalska-Jaruga, A.; Klimkowicz-Pawlas, A.; Smreczak, B. Characterization of organic matter fractions in the top layer of soils under different land uses in Central-Eastern Europe. Soil Use Manag. 2019, 35, 595–606. [Google Scholar] [CrossRef]
  2. Kibblewhite, M.G.; Ritz, K.; Swift, M.J. Soil health in agricultural systems. Philos. Trans. R. Soc. B 2008, 363, 685–701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Skłodowski, P.; Bielska, A. Properties and fertility of soils in Poland–a basis for the formation of agro-environmental relations. Water Environ. Rural Areas 2009, 9, 203–2014. (In Polish) [Google Scholar]
  4. Krasowicz, S.; Oleszek, W.; Horabik, J.; Dębicki, R.; Jankowiak, J.; Stuczyński, T.; Jadczyszyn, J. Rational managenent of the soil environment in Poland. Pol. J. Agron. 2011, 7, 43–58. (In Polish) [Google Scholar]
  5. Grządziel, J. Structural and Functional Diversity of Microorganisms in Selected Soil Types Characteristic of Poland. Ph.D. Thesis, Insitute of Soil Science and Plant Cultivation, State Research Institute, Puławy, Poland, 23 June 2020. [Google Scholar]
  6. Czyż, E.A.; Dexter, A.R.; Niedźwiecki, J. Soil water retention of long-term permanent plot experiments. Soil Sci. Ann. 2003, 4, 27–37. (In Polish) [Google Scholar]
  7. Gałązka, A.; Gawryjołek, K.; Grządziel, J.; Księżak, A. Effect of different agricultural management practices on soil biological parameters including glomalin fraction. Plant Soil Environ. 2017, 63, 300–306. [Google Scholar]
  8. Czaban, J.; Wróblewska, B.; Niedźwiecki, J.; Sułek, A. Relationships between numbers of microbial communities in Polish agricultural soils and properties of these soils, paying special attention to xerophilic/xerotolerant fungi. Pol. J. Environ. Stud. 2010, 19, 1171–1183. [Google Scholar]
  9. Grządziel, J.; Furtak, K.; Gałązka, A. Community-level physiological profiles of microorganisms from different types of soil that are characteristic to Poland—A long-term microplot experiment. Sustainability 2019, 11, 56. [Google Scholar] [CrossRef] [Green Version]
  10. Grządziel, J.; Gałązka, A. Microplot long-term experiment reveals strong soil type influence on bacteria composition and its functional diversity. Appl. Soil Ecol. 2018, 124, 117–123. [Google Scholar] [CrossRef]
  11. Grządziel, J.; Gałązka, A. Fungal biodiversity of the most common types of Polish soil in a long-term microplot experiment. Front. Microbiol. 2019, 10, 6. [Google Scholar] [CrossRef] [Green Version]
  12. Jamiołkowska, A.; Księżniak, A.; Hetman, B.; Kopacki, M.; Skwaryło-Bednarz, B.; Gałązka, A.; Thanoon, A.H. Interaction of arbuscular mycorrhizal fungi with plants and soil microflora. Acta Sci. Pol. Hortorum Cultus 2017, 16, 89–95. [Google Scholar] [CrossRef]
  13. Gałązka, A.; Gawryjołek, K. Glomalin-soil glicoprotein produced by arbuscular mycorhizal fungus. Postępy Mikrobiol. 2015, 54, 331–343. [Google Scholar]
  14. Gillespie, A.W.; Farrell, R.E.; Walley, F.L.; Ross, A.R.S.; Leinweber, P.; Eckhardt, K.U.; Tom, Z.; Regier, T.Z.; Blyth, R.I.R. Glomalin-related soil protein contains non-mycorrhizal-related heat-stable proteins, lipids and humic materials. Soil Biol. Biochem. 2011, 43, 766–777. [Google Scholar] [CrossRef]
  15. Wang, W.; Zhong, Z.; Wang, Q.; Wang, H.; Fu, Y.; He, X. Glomalin contributed more to carbon, nutrients in deeper soils, and differently associated with climates and soil properties in vertical profiles. Sci. Rep. 2017, 7, 13003. [Google Scholar] [CrossRef] [Green Version]
  16. Dexter, A.R. Soil physical quality: Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma 2004, 120, 201–214. [Google Scholar] [CrossRef]
  17. Dexter, A.R.; Czyz, E.A. Applications of S-theory in the study of soil physical degradation and its consequences. Land Degrad. Dev. 2007, 18, 369–381. [Google Scholar] [CrossRef]
  18. Gate, O.P.; Czyz, E.A.; Dexter, A.R. Soil physical quality, S, as a basis for relationships between some key physical properties of arable soils. Adv. Geoecology 2006, 38, 102–109. [Google Scholar]
  19. Strzemski, M. The History of Polish Soil Science from the Dawn of Polish Literature to the Establishment of the Second Polish Republic; PWRiL: Warszawa, Poland, 1980. (In polish) [Google Scholar]
  20. Siebielec, G.; Siebielec, S.; Podolska, G. Comparison of microbial and chemical characteristics of soil types after over 100 years of cereal production. Pol. J. Agron. 2015, 23, 88–100. [Google Scholar]
  21. Czyż, E.; Dexter, A.R. Mechanical dispersion of clay from soil into water: Readily-dispersed and spontaneously-dispersed clay. Int. Agrophys. 2015, 29, 31–37. [Google Scholar] [CrossRef]
  22. Genuchten, M.T.V. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  23. Martin, J.P. Use of acid, rose bengal and streptomycin in the plate method for estimating soil fungi. Soil Sci. 1950, 69, 215–232. [Google Scholar] [CrossRef]
  24. Rodina, A. Microbiological Methods for the Study of Water; PWRiL: Warszawa, Poland, 1968. (In Polish) [Google Scholar]
  25. Fenglerowa, W. Simple method for counting azotobacter in soil samples. Acta Microbiol. Pol. 1965, 14, 203–206. [Google Scholar]
  26. Ghani, A.; Dexter, M.; Perrott, K.W. Hot-water extractable carbon in soils: A sensitive measurement for determining impacts of fertilization, grazing and cultivation. Soil Biol. Biochem. 2003, 35, 1231–1243. [Google Scholar] [CrossRef]
  27. Determination of Dehydrogenase Activity in Soil Using 2,3,5 Triphenyltetrazolium Chloride (TTC); Polish Standard: Warsaw, Poland, 2011; PN-EN ISO 23753-1.
  28. Tabatabai, M.A. Soil enzymes. In Methods of Soil Analysis, Part 2; Page, A.L., Miller, R.H., Keeney, D.R., Eds.; American Society of Agronomy and Soil Science Society of America: Madison, WI, USA, 1982. [Google Scholar]
  29. Wright, S.F.; Franke-Synder, M.; Morton, J.B.; Upadhyaya, A. Time course study and partial characterization of a protein on hyphae of arbuscular mycorrhizal fungi during active colonization of roots. Plant Soil 1996, 181, 193–203. [Google Scholar] [CrossRef]
  30. Phillips, J.M.; Hayman, D.S. Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection. Trans. Br. Mycol. Soc. 1970, 55, 158–161. [Google Scholar] [CrossRef]
  31. Weisburg, W.G.; Barns, S.M.; Pelletier, D.A.; Lane, D.J. 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 1991, 173, 697–703. [Google Scholar] [CrossRef] [Green Version]
  32. Muyzer, G.; Waal, E.C.D.E.; Uitierlinden, A.G. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 1993, 59, 695–700. [Google Scholar] [CrossRef] [Green Version]
  33. Stach, J.E.M.; Maldonado, L.A.; Ward, A.C.; Goodfellow, M.; Bull, A.T. New primers for the class actinobacteria: Application to marine and terrestrial environments. Environ. Microbiol. 2003, 5, 828–841. [Google Scholar] [CrossRef] [Green Version]
  34. Gałązka, A.; Gawryjołek, K.; Grządziel, J.; Frąc, M.; Księżak, J. Microbial community diversity and the interaction of soil under maize growth in different cultivation techniques. Plant Soil Environ. 2017, 63, 264–270. [Google Scholar] [CrossRef] [Green Version]
  35. Chau, J.F.; Bagtzoglou, A.C.; Willig, M.R. The effect of soil texture on richness and diversity of bacterial communities. Environ. Forensics 2011, 12, 333–341. [Google Scholar] [CrossRef]
  36. Nikolausz, M.; Sipos, R.; Révész, S.; Székely, A.; Márialigeti, K. Observation of bias associated with re-amplification of DNA isolated from denaturing gradient gels. FEMS Microbiol. Lett. 2005, 244, 385–390. [Google Scholar] [CrossRef] [PubMed]
  37. Furtak, K.; Gawryjołek, K.; Gajda, A.; Gałązka, A. Effects of maize and winter wheat grown under different cultivation techniques on biological activity of soil. Plant Soil Environ. 2017, 63, 449–454. [Google Scholar]
  38. Gałązka, A.; Gawryjołek, K.; Gajda, A.; Furtak, K.; Księżniak, A.; Jończyk, K. Assessment of the glomalins content in the soil under winter wheat in different crop production systems. Plant Soil Environ. 2018, 64, 32–37. [Google Scholar]
  39. Latour, X.; Philippot, L.; Corberand, T.; Lemanceau, P. The establishment of an introduced community of fluorescent pseudomonads in the soil and in the rhizosphere is affected by the soil type. FEMS Microbiol. Ecol. 1999, 30, 163–170. [Google Scholar] [CrossRef] [PubMed]
  40. Schreiter, S.; Ding, G.-C.; Heuer, H.; Neumann, G.; Sandmann, M.; Grosch, R.; Kropf, S.; Smalla, K. Effect of the soil type on the microbiome in the rhizosphere of field-grown lettuce. Front. Microbiol. 2014, 5, 1–13. [Google Scholar] [CrossRef] [PubMed]
  41. Sessitsch, A.; Weilharter, A.; Gerzabek, M.H.; Kirchmann, H.; Kandeler, E. Microbial population structures in soil particle size fractions of a long-term fertilizer field experiment. Appl. Environ. Microbiol. 2001, 67, 4215–4224. [Google Scholar] [CrossRef] [Green Version]
  42. Lehman, R.M.; Cambardella, C.A.; Stott, D.E.; Acosta-Martinez, V.; Manter, D.K.; Buyer, J.S.; Maul, J.E.; Smith, J.L.; Collins, H.P.; Halvorson, J.J.; et al. Understanding and enhancing soil biological health: The solution for reversing soil degradation. Sustainability 2015, 7, 988–1027. [Google Scholar] [CrossRef] [Green Version]
  43. Gaţe, O.P.; Czyż, E.A.; Dexter, A.R. Effects of readily-dispersible clay on soil quality and root growth. In Plant Growth in Relation to Soil Physical Conditions; Lipiec, J., Walczak, R., Józefaciuk, G., Eds.; Institute of Agrophysics, Polish Academy of Science: Lublin, Poland, 2004. [Google Scholar]
  44. Lal, R. Soil health and carbon management. Food Energy Secur. 2016, 5, 212–222. [Google Scholar] [CrossRef]
  45. Norris, C.E.; Bean, G.M.; Cappellazzi, S.B.; Cope, M.; Greub, K.L.H.; Liptzin, D.; Rieke, E.L.; Tracy, P.W.; Morgan, C.L.S.; Honeycutt, C.W. Introducing the North American project to evaluate soil health. Meas. Agron. J. 2020, 1–21. [Google Scholar] [CrossRef]
  46. He, F.; Chen, H.; Tang, M. Arbuscular mycorrhizal fungal communities are influenced by host tree species on the loess plateau, northwest China. Forests 2019, 10, 930. [Google Scholar] [CrossRef] [Green Version]
  47. Verzeaux, J.; Nivelle, E.; Roger, D.; Hirel, B.; Dubois, F.; Tetu, T. Spore density of arbuscular mycorrhizal fungi is fostered by six years of a No-Till system and is correlated with environmental parameters in a silty loam soil. Agronomy 2017, 7, 38. [Google Scholar] [CrossRef]
  48. Melo, C.D.; Walker, C.; Krüger, C.; Borges, P.A.V.; Luna, S.; Mendonça, D.; Fonseca, H.M.A.C.; Machado, A.C. Environmental factors driving arbuscular mycorrhizal fungal communities associated with endemic woody plant Picconia azorica on native forest of Azores. Ann. Microbiol. 2019, 69, 1309–1327. [Google Scholar] [CrossRef]
  49. Onet, A.; Dinca, P.; Laslo, V.; Teusdea, A.C.; Vasile, D.L.; Enescu, R.E.; Crisan, V.E. Biological indicators for evaluating soil quality improvement in a soil degraded by erosion processes. J. Soils Sediments 2019, 19, 2393–2404. [Google Scholar] [CrossRef]
  50. Nannipieri, P.; Ascher, J.; Ceccherini, M.T.; Landi, L.; Pietramellara, G.; Renella, G. Microbial diversity and soil functions. Eur. J. Soil Sci. 2003, 54, 655–670. [Google Scholar] [CrossRef]
  51. Girvan, M.S.; Bullimore, J.; Pretty, J.N.; Osborn, M.; Ball, A.S. Soil type is the primary determinant of the composition of the total and active bacterial communities in arable soils. Appl. Environ. Microbiol. 2003, 69, 1800–1809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Wolińska, A.; Kuźniar, A.; Zielenkiewicz, U.; Banach, A.; Błaszczyk, M. Indicators of arable soils fatigue–bacterial families and genera: A metagenomic approach. Ecol. Ind. 2018, 93, 490–500. [Google Scholar] [CrossRef]
  53. Sarapatka, B.; Alvarado-Solano, D.P.; Čižmár, D. Can glomalin content be used as an indicator for erosion damage to soil and related changes in organic matter characteristics and nutrients? Catena 2019, 181, 104078. [Google Scholar] [CrossRef]
  54. Seaton, F.M.; George, P.B.L.; Lebron, I.; Jones, D.L.J.; Creer, S.; Robinson, D.A. Soil textural heterogeneity impacts bacterial but not fungal diversity. Soil Biol. Biochem. 2020, 144, 107766. [Google Scholar] [CrossRef]
  55. Lobry de Bruyn, L.; Andrews, S. Are Australian and United States farmers using soil information for soil health management? Sustainability 2016, 8, 304. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The microplots experiment (A). The microplots are presented in the period before sowing, in early spring, April (B). (photo by J. Grządziel, materials from Ph.D. thesis), [5].
Figure 1. The microplots experiment (A). The microplots are presented in the period before sowing, in early spring, April (B). (photo by J. Grządziel, materials from Ph.D. thesis), [5].
Agronomy 10 01279 g001
Figure 2. The scheme of soil sampling from the M1 plot.
Figure 2. The scheme of soil sampling from the M1 plot.
Agronomy 10 01279 g002
Figure 3. PCR-DGGE fingerprinting of 16S rRNA gene: (A) Bacteria V3 regions of 16S rRNA gene; (B) Actinobacteria regions of 16S RNA gene. SL—unrelated standard lanes. Bands that are unique for an indicated lane are marked with an asterisk. Two examples of bands that are specific only to good-quality soils are marked with an arrow. The soil symbols were given in Table 1.
Figure 3. PCR-DGGE fingerprinting of 16S rRNA gene: (A) Bacteria V3 regions of 16S rRNA gene; (B) Actinobacteria regions of 16S RNA gene. SL—unrelated standard lanes. Bands that are unique for an indicated lane are marked with an asterisk. Two examples of bands that are specific only to good-quality soils are marked with an arrow. The soil symbols were given in Table 1.
Agronomy 10 01279 g003
Figure 4. Similarity matrix of PCR-DGGE fingerprint with cluster analysis. Scale bar represents similarity in percentages: from ~60% to 100% (A) and from ~30% to 100% (B).
Figure 4. Similarity matrix of PCR-DGGE fingerprint with cluster analysis. Scale bar represents similarity in percentages: from ~60% to 100% (A) and from ~30% to 100% (B).
Agronomy 10 01279 g004
Figure 5. Enzymatic activities in soil, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Dehydrogenase activity; (B) Acid phosphatase activity; (C) Alkaline phosphatase activity.
Figure 5. Enzymatic activities in soil, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Dehydrogenase activity; (B) Acid phosphatase activity; (C) Alkaline phosphatase activity.
Agronomy 10 01279 g005
Figure 6. Microbial biomass analysis, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Microbial Biomass of Carbon (MBC); (B) Microbial Biomass of Nitrogen (MBN).
Figure 6. Microbial biomass analysis, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Microbial Biomass of Carbon (MBC); (B) Microbial Biomass of Nitrogen (MBN).
Agronomy 10 01279 g006
Figure 7. Total number of bacteria and fungi analysis, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Total number of bacteria; (B) Total number of fungi.
Figure 7. Total number of bacteria and fungi analysis, ANOVA analysis (the soils were organised in this case from highest to lowest soil quality): (A) Total number of bacteria; (B) Total number of fungi.
Agronomy 10 01279 g007
Figure 8. Dendrogram of the bond distances between of selected microbiological, chemical and biochemical parameters and GRSP content (see Table 9).
Figure 8. Dendrogram of the bond distances between of selected microbiological, chemical and biochemical parameters and GRSP content (see Table 9).
Agronomy 10 01279 g008
Figure 9. Principal component analysis of microbial parameters, biodiversity indices and Biolog EcoPlate data from soil samples incubated for 120 h. See Table 11.
Figure 9. Principal component analysis of microbial parameters, biodiversity indices and Biolog EcoPlate data from soil samples incubated for 120 h. See Table 11.
Agronomy 10 01279 g009
Table 1. Classification and basic characteristic of soils used in the experiment (the soils were grouped according to the S Index value from the lowest value to the highest).
Table 1. Classification and basic characteristic of soils used in the experiment (the soils were grouped according to the S Index value from the lowest value to the highest).
Soil Type 1SymbolPlot NumberSoil Texture Class 2Particle Size Distribution [%]S-IndexT-GRSPC [%]N [%]C/NRDC/TCpH
Sand
(2.0–0.05 mm)
Silt
(0.05–0.002 mm)
Clay
(<0.002 mm)
Brunic Arenosol (Dystric I)BA (I)M1loamy sand782110.0182.870.6540.06310.380.3134.0
Brunic Arenosol (Dystric II)BA (II)M6loamy sand851410.0253.070.5220.04611.350.2884.5
Fluvic CambisolFCM3sandy loam5234140.0464.591.3640.10612.870.0547.5
Rendzic LeptosolRLM2sandy loam603280.0504.431.9670.11517.100.0847.3
Haplic Cambisol (Dystric)HC (D)M8fine sandy loam712720.0563.790.9160.09010.180.1265.6
Haplic Cambisol (Eutric II)HC (EII)M7medium sandy loam583840.0584.680.8580.0919.430.0985.4
Haplic Cambisol (Eutric I)—developed from loessHC (EI)M4fine sandy loam494830.0604.040.8420.0899.460.1794.7
Gleyic PhaeozemGPhM5fine sandy loam544330.0766.333.1010.22513.780.0957.4
1 Determined in accordance with the World Reference Base for Soil Resources (WRB, 2008). 2 Determined on the basis of the USDA (United States Department of Agriculture) Natural Resources Conservation Service; Soil Taxonomy A Basic System of Soil Classification for Making and Interpreting Soil Surveys. Particle size distribution [%]—Casagrande method modified by Prószyński [PN—R 04032:1998] T-GRSP—total glomalin related soil proteins [mg/g d.m. of soil]. C [%], N [%]-Determinations of carbon and nitrogen by C-N Analyzer RDC/TC—readily dispersible clay (RDC)/total clay (TC). pH H2O-potentiometric method-[PN-ISO 10390:1997]. S-Index classification: S ≥ 0.050 very good soil physical quality; 0.050 > S ≥ 0.035 good soil physical quality; 0.035 > S ≥ 0.020 poor soil physical quality; 0.020 > S very poor soil physical quality.
Table 2. Meteorological conditions during the growing seasons (2016–2018).
Table 2. Meteorological conditions during the growing seasons (2016–2018).
Average Air Temperature (°C)
IIIIIIIVVVIVIIVIIIIXXXIXIIYear
20161.31.05.18.613.517.920.422.615.37.05.24.010.2
2017−3.33.74.39.615.619.820.118.715.67.73.20.89.7
2018−4.6−1.16.27.914.718.919.420.114.19.54.62.59.4
Sum of Precipitation (mm)
IIIIIIIVVVIVIIVIIIIXXXIXIIYear
201655.39.945.828.5111.632.055.43.6126.030.047.425.2570.7
201733.064.553.038.472.227.986.641.920.6100.244.964.6647.8
201810.742.431.631.948.235.4150.376.8105.294.953.720.5731.6
Table 3. Primers used in the experiment.
Table 3. Primers used in the experiment.
Primer SetSequence 5′–3′TargetReference
Bacteria-specific
27F
1492R
>AGAGTTTGATCCTGGCTCAG
>GGTTACCTTGTTACGACTT
16S rRNA gene[31]
GC-338F
Eub518R
>GC-clamp-ACTCCTACGGGAGGCAGCAG
>ATTACCGCGGCTGCTGG
V3 hypervariable region[32]
Actinobacteria-specific
Actino235F
Eub518R
>GC-clamp-CGCGGCCTATCAGCTTGTTG
>ATTACCGCGGCTGCTGG
16S rRNA gene[33]
GC-clamp sequence
> CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGG (40 nt)[32]
Table 4. PCR and PCR-DGGE amplification parameters.
Table 4. PCR and PCR-DGGE amplification parameters.
PCR Step16S rRNA PCRV3 PCR-DGGEActino PCR-DGGE
Initial denaturation95 °C, 15 min95 °C, 3 min95 °C, 3 min
Denaturation95 °C, 30 s95 °C, 30 s95 °C, 30 s
Annealing52 °C, 60 s55 °C, 30 s55 °C, 30 s
Extension72 °C, 90 s72 °C, 30 s72 °C, 30 s
Final extension72 °C, 7 min72 °C, 5 min72 °C, 5 min
Number of cycles353232
Table 5. Mean values of parameters of the van Genuchten equation [20] and of the S-Index (the soils were grouped according to the S Index value from the lowest value to the highest).
Table 5. Mean values of parameters of the van Genuchten equation [20] and of the S-Index (the soils were grouped according to the S Index value from the lowest value to the highest).
Soil Typeθs
(kg/kg)
θr
(kg/kg)
α
(hPa−1)
n
(-)
m
(-)
S
(-)
Brunic Arenosol (Dystric I)0.25280.01180.01011.68890.4079−0.0185
Brunic Arenosol (Dystric II)0.23600.02020.01161.57550.3652−0.0250
Fluvic Cambisol0.25280.03490.00151.68890.4079−0.0456
Rendzic Leptosol0.28040.05040.00771.87520.4667−0.0501
Haplic Cambisol (Dystric)0.24920.03490.01011.66460.3992−0.0562
Haplic Cambisol (Eutric II)0.24560.03360.00671.64030.3903−0.0581
Haplic Cambisol (Eutric I)0.25400.02520.00671.69700.4107−0.0603
Gleyic Phaeozem0.26840.08070.01011.79420.4426−0.0764
Explanation: θs—water content at saturation; θr—residual water content; α—adjustable scaling factor; m, n—adjustable shape factors; S—soil physical quality index.
Table 6. Results of an ANOVA analysis of variance (p ≤ 0.05) for the entire data set. Multidimensional significance tests. Parameterization with sigma-restrictions-decomposition of effective hypotheses.
Table 6. Results of an ANOVA analysis of variance (p ≤ 0.05) for the entire data set. Multidimensional significance tests. Parameterization with sigma-restrictions-decomposition of effective hypotheses.
Wilks ValueF Valuep
Type of soil0.00035.340.000
Year0.1526.970.000
Type of soil_Year0.0052.530.000
Table 7. Total and easily extractable glomalin related soil proteins content in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest).
Table 7. Total and easily extractable glomalin related soil proteins content in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest).
Total Glomalin Related Soil Proteins (T-GRSP), (mg/g d.m. of Soil)
Year201620172018Mean
Brunic Arenosol (Dystric I)3.273.232.122.87 c
Brunic Arenosol (Dystric II)2.583.503.143.07 c
Fluvic Cambisol4.564.494.714.59 b
Rendzic Leptosol5.134.563.604.43 b
Haplic Cambisol (Dystric)3.924.033.453.79 b,c
Haplic Cambisol (Eutric II)4.684.824.564.68 b
Haplic Cambisol (Eutric I)4.234.253.654.04 b
Gleyic Phaeozem6.436.595.986.33 a
Easily Extractable Glomalin Related Soil Proteins (EE-GRSP), (mg/g d.m. of Soil)
Brunic Arenosol (Dystric I)1.721.541.621.63 c
Brunic Arenosol (Dystric II)1.721.771.541.67 c
Fluvic Cambisol2.002.001.911.97 c
Rendzic Leptosol2.462.212.062.24 a
Haplic Cambisol (Dystric)1.952.101.972.00 b
Haplic Cambisol (Eutric II)2.672.462.182.44 a
Haplic Cambisol (Eutric I)2.432.292.222.31 a
Gleyic Phaeozem2.322.582.252.38 a
Treatment means separated by different letters are significantly different (Tukey’s mean separation test, p < 0.05); average for years 2016–2018 (n = 9)
Table 8. Number of AMF spores in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest). Treatment means separated by different letters are significantly different (Tukey’s mean separation test, p < 0.05); average for years 2016–2018 (n = 9).
Table 8. Number of AMF spores in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest). Treatment means separated by different letters are significantly different (Tukey’s mean separation test, p < 0.05); average for years 2016–2018 (n = 9).
Soil TypeTotal Number of All SporesTotal Number by Spore Diameter
50 µm100 µm150 µm
Brunic Arenosol (Dystric I)188 d126 c54 c8 c
Brunic Arenosol (Dystric II)406 c296 b92 c18 b
Fluvic Cambisol534 a312 b184 a38 a
Rendzic Leptosol374 c284 b76 c14 b
Haplic Cambisol (Dystric)592 a396 a172 a24 b
Haplic Cambisol (Eutric II)586 a324 b204 a58 a
Haplic Cambisol (Eutric I)464 b360 b88 c16 c
Gleyic Phaeozem468 b308 b116 b44 a
Table 9. Microbial diversity in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest). Treatment means separated by different letters are significantly different (Tukey’s mean separation test, p < 0.05); average for years 2016–2018 (n = 9).
Table 9. Microbial diversity in different soil types (the soils were grouped according to the S Index value from the lowest value to the highest). Treatment means separated by different letters are significantly different (Tukey’s mean separation test, p < 0.05); average for years 2016–2018 (n = 9).
Soil TypeDHAAlPAcPMBCMBNBacteriaFungiAmoAzo
Brunic Arenosol (Dystric I)5.524 d8.525 d34.628 b53.867 c4.074 c12.949 c16.906 b6.475 b0.000 c
Brunic Arenosol (Dystric II)5.401 d8.648 d41.257 b83.378 b7.476 b0.352 d21.795 a2.109 c0.000 a
Fluvic Cambisol70.478 b77.268 b30.418 b59.378 c17.418 a49.772 a6.079 c7.219 b199.467 a
Rendzic Leptosol93.348 a90.357 b37.644 b150.222 a19.268 a39.089 b11.114 b11.497 a173.985 a
Haplic Cambisol (Dystric)24.968 c17.038 c57.903 a67.556 b5.150 c26.862 b27.936 a12.177 a0.000 c
Haplic Cambisol (Eutric II)40.641 c28.975 c54.069 a45.689 c7.990 b20.862 b8.196 c7.078 b1.863 b
Haplic Cambisol (Eutric I)44.232 c29.256 c56.882 a39.911 c8.953 b54.072 a7.184 c6.806 a1.513 b
Gleyic Phaeozem115.036 a128.880 a43.044 a138.578 a20.185 a65.961 a8.294 c9.874 b187.218 a
Explanation: DHA—dehydrogenases (μg formazan/g d.m. (dry matter) of soil/24 h); AlP—alkaline phosphatase (μg p-nitrophenol/g d.m. of soil/h); AcP—acid phosphatase (μg p-nitrophenol/g d.m. of soil/h); MBC—microbial biomass C (μg/g d.m. of soil); MBN—microbial biomass N (μg/g d.m. of soil); bacteria (108 CFU (colony-forming units)/g d.m. of soil); fungi (104 CFU/g d.m. of soil); Am—ammonification bacteria (106 CFU/g d.m. of soil); Azo—Azotobacter spp. (101 CFU/g d.m. of soil).
Table 10. Pearson’s correlation coefficient (p ≤ 0.05) of selected microbiological, chemical and biochemical parameters and glomalin content.
Table 10. Pearson’s correlation coefficient (p ≤ 0.05) of selected microbiological, chemical and biochemical parameters and glomalin content.
AlPAcPMBCMBNBacteriaFungiAmoAzoN [%]C [%]RDC/TGS-IndexT-GRSPEE-GRSPTotal Number of the Spores (µm)All AMF Spores
15010050
DHA0.979 *−0.2000.7030.945 *0.836 *−0.6220.5510.878 *0.882 *0.934 *0.734 *0.725 *0.897 *0.6790.3820.0910.2250.224
AlP −0.3430.744 *0.960 *0.789 *−0.5730.4750.925 *0.899 *0.962 *0.6990.6250.876 *0.5390.3610.0440.1120.133
AcP−0.343 −0.328−0.4400.0250.2250.213−0.592−0.044−0.222−0.4500.4960.0490.4820.1920.3030.6680.554
MBC0.744 *−0.328 0.731*0.315−0.0530.4400.6690.6220.796 *0.3110.2870.4970.226−0.051−0.319−0.034−0.162
MBN0.960 *−0.4400.731 * 0.723 *−0.6110.3610.959 *0.749 *0.861 *0.822 *0.5120.771 *0.4750.2860.0260.1310.127
Bacteria0.789 *0.0250.3150.723 * −0.6470.5020.6730.799 *0.745 *0.6940.785 *0.795 *0.6610.2750.1000.3510.289
Fungi−0.5730.225−0.053−0.611−0.647 0.035−0.530−0.443−0.422−0.775 *−0.417−0.593−0.599−0.430−0.0950.036−0.075
Amo0.4750.2130.4400.3610.5020.035 0.3910.5230.5390.2290.5880.4590.4750.0820.1860.3040.272
Azo0.925 *−0.5920.6690.959 *0.673−0.5300.391 0.6950.813 *0.817 *0.3790.6860.2950.2450.0700.0150.071
N [%]0.899 *−0.0440.6220.749 *0.799 *−0.4430.5230.695 0.960 *0.3940.775 *0.937 *0.6170.4630.1020.1890.217
C [%]0.962−0.2220.796 *0.861 *0.745 *−0.4220.5390.813 *0.960 * 0.4930.6590.879 *0.5330.328−0.0290.0940.087
RDC/TG0.699−0.4500.3110.822 *0.694−0.775 *0.2290.817 *0.3940.493 0.3310.4950.4000.1450.0630.0950.103
S-Index0.6250.4960.2870.5120.785 *−0.4170.5880.3790.775 *0.6590.331 0.867 *0.905 *0.5880.4240.6820.666
T-GRSP0.876 *0.0490.4970.771 *0.795 *−0.5930.4590.6860.937 *0.879 *0.4950.867 * 0.784 *0.790 *0.3530.3600.656
EE-GRSP0.5390.4820.2260.4750.661−0.5990.4750.2950.6170.5330.4000.905 *0.784 * 0.6000.3650.5460.560
150 µm0.3610.192−0.0510.2860.275−0.4300.0820.2450.4630.3280.1450.5880.790 *0.600 0.809 *0.3740.698
100 µm0.0440.303−0.3190.0260.100−0.0950.1860.0700.102−0.0290.0630.4240.3530.3650.809 * 0.5820.878 *
50 µm0.1120.668−0.0340.1310.3510.0360.3040.0150.1890.0940.0950.6820.3600.5460.3740.582 0.895 *
All spores0.1330.554−0.1620.1270.289−0.0750.2720.0710.2170.0870.1030.6660.6560.5600.6980.878 *0.895 *
* Indicates statistically significant (p ≤ 0.05); average of 2016–2018 (n = 9) Corg—organic carbon content (g/kg of soil); N (g/kg of soil); see Table 9.
Table 11. Results of principal component analysis.
Table 11. Results of principal component analysis.
VariablePC1 (48.35%)PC2 (23.38%)
DHA−0.9740.164
AlP−0.9310.243
AcP0.010−0.783
MBC−0.6010.149
MBN−0.8630.405
Bacteria−0.886−0.165
Fungi0.578−0.383
Amo−0.634−0.340
Azo−0.7910.487
N [%]−0.908−0.098
C [%]−0.9070.051
S-Index−0.851−0.366
RDC/TC0.827−0.160
T-GRSP−0.9380.032
EE-GRSP−0.768−0.201
Amines and amides−0.197−0.873
Carboxylic and acetic acids−0.218−0.893
Carbohydrates−0.214−0.596
Amino acids0.719−0.392
Polymers−0.393−0.513
AWCD−0.024−0.889
Shannon index−0.477−0.624
Correlation coefficients for principal component (PC) factors; average for years 2016–2018 (n = 9). For abbreviations, see Table 9.

Share and Cite

MDPI and ACS Style

Gałązka, A.; Niedźwiecki, J.; Grządziel, J.; Gawryjołek, K. Evaluation of Changes in Glomalin-Related Soil Proteins (GRSP) Content, Microbial Diversity and Physical Properties Depending on the Type of Soil as the Important Biotic Determinants of Soil Quality. Agronomy 2020, 10, 1279. https://doi.org/10.3390/agronomy10091279

AMA Style

Gałązka A, Niedźwiecki J, Grządziel J, Gawryjołek K. Evaluation of Changes in Glomalin-Related Soil Proteins (GRSP) Content, Microbial Diversity and Physical Properties Depending on the Type of Soil as the Important Biotic Determinants of Soil Quality. Agronomy. 2020; 10(9):1279. https://doi.org/10.3390/agronomy10091279

Chicago/Turabian Style

Gałązka, Anna, Jacek Niedźwiecki, Jarosław Grządziel, and Karolina Gawryjołek. 2020. "Evaluation of Changes in Glomalin-Related Soil Proteins (GRSP) Content, Microbial Diversity and Physical Properties Depending on the Type of Soil as the Important Biotic Determinants of Soil Quality" Agronomy 10, no. 9: 1279. https://doi.org/10.3390/agronomy10091279

APA Style

Gałązka, A., Niedźwiecki, J., Grządziel, J., & Gawryjołek, K. (2020). Evaluation of Changes in Glomalin-Related Soil Proteins (GRSP) Content, Microbial Diversity and Physical Properties Depending on the Type of Soil as the Important Biotic Determinants of Soil Quality. Agronomy, 10(9), 1279. https://doi.org/10.3390/agronomy10091279

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

Article Metrics

Back to TopTop