Soil Microbial Community Proﬁling and Bacterial Metabolic Activity of Technosols as an E ﬀ ect of Soil Properties following Land Reclamation: A Case Study from the Abandoned Iron Sulphide and Uranium Mine in Rudki (South-Central Poland)

: The aims of the study were (1) to recognize the structure of bacteria diversity in Technosols developed from mine spoils containing iron (Fe) sulphides with the use of culture-independent technique, and (2) to determine microbial metabolic activities, in the context of their potential to be an adequate indicators of soil properties being the consequence of land reclamation. The study site was located in the vicinity of the abandoned Fe sulphide and uranium mine in Rudki village (Holy Cross Mts., Poland). Three soil proﬁles with di ﬀ erent chemical properties (pH, content of carbonates, soil salinity, content of total organic carbon and total nitrogen) were studied. Biodiversity was determined with the use of meta-barcoding of 16S rRNA community proﬁling analysis based on the hypervariable V3-V4 region of 16S rRNA gene (MiSeq, Illumina). The catabolic ﬁngerprinting of soil microbial communities was evaluated with the use of Biolog ® EcoPlates ™ System. It was evidenced that changes in microbial structure and their metabolic activity were the consequence of a combined e ﬀ ect of both the soil depth and soil chemical properties being the ﬁnal result of reclamation process. Consequently, microbial indicators (from phyla to genera level) indirectly testifying about success or ine ﬀ ectiveness of reclamation in technogenic soils were recommended. To our best knowledge, the present study is the ﬁrst insight into Polish Technosols biodiversity and catabolic activity.


Introduction
Technosols are soils strongly influenced by human activity or made by human [1]. They develop e.g., on mine waste disposal sites in mining areas. A specific group of Technosols are soils containing iron (Fe) sulphides. In Poland, such Technosols occur e.g., on mine spoils of abandoned Fe sulphide mines in Rudki and Wieściszowice [2].

Study Site Description and Soil Sampling
The study site was located in the vicinity of the abandoned "Staszic" mine in Rudki village (south-central Poland, 50°53′52″ N, 21°05′55″ E). Iron sulphides and uranium ores were exploited from 1920 to 1970s [2,3,5]. Three soil profiles (R1 -Rudki 1, R2 -Rudki 2, R3 -Rudki 3) developed from mine spoils containing iron sulphides and representing areas of different effects of land reclamation were investigated (Scheme 1). Detailed description of soil profiles and their properties are available elsewhere [5,39]. R1 profile was located on the former post-flotation waste tanks. These wastes were originally composed of crushed ore-bearing rock consisting mainly of dolomite and Fe sulphides, and were characterized by strong compaction. The reclamation performed in 1970s included covering the tanks with a layer of loamy material of a thickness about 40-50 cm [40]. Such a process led to the formation of soil profile with a layer of post-flotation waste (2C horizon, 45-85 cm) in the subsoil and a loamy layer in the topsoil where the following horizons were distinguished during the field work: A (1-3 cm), C1 (3-15 cm), and C2 (15-45 cm). R1 profile had high pH (6.7-7.5) throughout the profile due to carbonates presence [5] (Table S1 Supplementary material). The highest concentrations of total organic carbon (TOC) and total nitrogen (TN) were present in A and C1 horizon (4.8 and 2.7% of TOC, and 0.41 and 0.23% of TN, respectively) with decreasing trend with profile depth [5,39]. Importantly, loamy material covering post-flotation sludge contains high concentration of trace elements (e.g., Pb, As, Tl, U, Th) and consequently surface layer of R1 is enriched in these elements [2,3,5,41]. Nevertheless, the topsoil of R1 profile comprises an optimal niche for microbial colonization, whereas the subsoil (2C horizon) is not a favourable environment for development of majority of microbes [5]. The soil surface was covered with a meadow plant community dominated by Medicago sp., Rubus sp., Hypericum perforatum and Calamagrostis epigejos [5]. R1 profile was located on the former post-flotation waste tanks. These wastes were originally composed of crushed ore-bearing rock consisting mainly of dolomite and Fe sulphides, and were characterized by strong compaction. The reclamation performed in 1970s included covering the tanks with a layer of loamy material of a thickness about 40-50 cm [40]. Such a process led to the formation of soil profile with a layer of post-flotation waste (2C horizon, 45-85 cm) in the subsoil and a loamy layer in the topsoil where the following horizons were distinguished during the field work: A (1-3 cm), C1 (3-15 cm), and C2 (15-45 cm). R1 profile had high pH (6.7-7.5) throughout the profile due to carbonates presence [5] (Table S1 Supplementary material). The highest concentrations of total organic carbon (TOC) and total nitrogen (TN) were present in A and C1 horizon (4.8 and 2.7% of TOC, and 0.41 and 0.23% of TN, respectively) with decreasing trend with profile depth [5,39]. Importantly, loamy material covering post-flotation sludge contains high concentration of trace elements (e.g., Pb, As, Tl, U, Th) and consequently surface layer of R1 is enriched in these elements [2,3,5,41]. Nevertheless, the topsoil of R1 profile comprises an optimal niche for microbial colonization, whereas the subsoil (2C horizon) is not a favourable environment for development of majority of microbes [5]. The soil surface was covered with a meadow plant community dominated by Medicago sp., Rubus sp., Hypericum perforatum and Calamagrostis epigejos [5].
R2 profile was located in the area of the former "Serwis" disposal site where strongly acidic mine wastes containing Fe sulphides were deposited [5]. The reclamation works made there in 1970s included liming of waste surface and introduction of grasses and legumes [40]. R2 profile represents an area where reclamation failed as the soil was not successfully neutralized. The reaction Agronomy 2020, 10, 1795 4 of 22 throughout the profile was strongly acidic (pH H2O 3.0-3.9) and carbonates did not occur [5,39] (Table S1, Supplementary material). That profile was poor in TOC and TN, therefore it is an extreme niche for microbial colonization [5]. Surface of soil was covered by sparse meadow [2,3,5,42].
R3 profile was located nearby R2 profile, however the R3 profile represents an area where reclamation process was successful. R3 profile had pH H2O 6.7 and 7.5 in A (2-10 cm) and AC (10-20 cm) horizons, respectively [39] (Table S1, Supplementary material). It was due to the presence of carbonates added to the soil during reclamation in 1970s [5,39]. However, carbonates were absent in deeper horizons of R3 profile, i.e., C1 (20-40 cm) and C2 (40-70 cm), therefore these horizons were strongly acidic (pH H2O of 3.9) due to the presence of Fe sulphides and products of their weathering [2][3][4][5]39]. Also, TOC and TN concentrations were high only in A horizon (8% of TOC and 0.65% of TN) and they decreased downward in the profile [5]. The soil was covered with a lush meadow. Thus, R3 profile represents an area where reclamation was partially successful because the acidity was neutralized in the topsoil, (however it still remained in the subsoil) and plans were successfully introduced to the area. The soil surface was covered with a lush meadow with Alopecurus pratensis, Urtica dioica, Festuca sp. and Artemisia vulgaris as predominant plant species [5].
Profiles R1 and R3 represent areas where reclamation can be treated as successful, however the properties of soils are not perfect, because both profiles have some favourable and unfavourable properties. Profile R1 is characterized by pH~7 and presence of carbonates throughout the profile, but its unfavourable feature is the high concentration of trace elements (e.g., Pb) in the topsoil (0-45 cm) and high compaction of soil material in the subsoil (2C horizon, 45-85 cm). Profile R3 has favourable properties in the topsoil (0-20 cm) (e.g., pH~7, presence of carbonates, high TOC contents), whereas unfavourable ones in the subsoil (<20 cm) (e.g., pH < 4). Finally, R2 profile represents an area where reclamation failed and soil properties are unfavourable (pH < 4, high soil salinity) throughout the profile.
The soil samples were taken in June 2019. The samples from each soil horizon distinguished in the field were mixed in the plastic container. In total, 13 samples were taken and studied in terms of microbiological properties. Both containers and tools for soil profiles taking were washed with water and ethanol (>99% w/w) every time prior to taking the next samples [5]. Plant roots were removed prior to laboratory analysis and soil materials were sieved with 2 mm sieve. The soil samples were stored in darkness at 4 • C for a short time (no longer than 24 h) prior to analyses.

DNA Extraction and Next Generation Sequencing (NGS)
All homogenous soil materials (250 mg in triplicate for each soil horizon) were applied for three independent DNA extractions performed with the Power Soil DNA Isolation Kit (QIAGEN, Hilden, Germany) according manufacturer's protocol. After that, the quantity and quality of isolated DNA were verified with a BioSpectrophotometer (Eppendorf, Hamburg, Germany). In the case if three DNA replicates were characterized by similar DNA content and purity, the extracted materials were pooled and mixed well in a single tube. This procedure remains in agreement with Soliman et al. [43] and Song et al. [44], who proved that pooling DNA extractions from individual soil samples increased operational taxonomic unit (OTU) richness.
Meta-barcoding of 16S rRNA community profiling analysis was performed based on the hypervariable V3-V4 region of 16S rRNA gene [9]. The following primers 341F and 785R were used for amplification of mentioned region and the library preparation [45]. The PCR reaction was carried out using a Q5 Hot Start High-Fidelity 2X Master Mix (New England Biolabs Inc., Ipswich, MA, USA) under conditions described elsewhere [9]. Next generation sequencing was performed by Genomed S.A. (Warsaw, Poland) on a MiSeq sequences (Illumina, San Diego, CA, USA) in paired-end (PE) technology, 2 × 300 nt, using Illumina v3 kit (San Diego, CA, USA).

Bioinformatic Analysis
DADA2 version 1.14 package [46] and R version 3.6.0 [47] were used to amplicon sequence variants (ASVs) analysis. The forward and reverse reads were trimmed to 250 bp based on the quality plots, and primers/adapters were removed from all analysed reads. Filtering parameters applied were as follows: maxN = 0, maxEE = 3, truncQ = 2. Sequences were also dereplicated using derepFastq with default parameters and to exact sequences were resolved using DADA [48]. Summary of the sequencing data quality obtained in the study is presented in Table S2. As a result of NGS sequencing of 13 studied soil samples 3091632 raw sequences were obtained in total. 2325317 sequences remained after the preliminary quality filtering. Next, denoised F/R quality filtering was performed, which yielded 2278140 (denoised F) and 2288207 (denoised R) sequences. The total amount of merged forward-reverse reads was 2163594 sequences. After chimera removal, 2105152 sequences remained for the analysis, meaning that ca. 32% did not meet the assumed criteria and were removed during restrictive bioinformatic data processing. Rarefaction curves were generated for all samples and did plateau for each sample studied ( Figure S1).

Bacterial Community-Level Physiological Profiles-Biolog ® EcoPlates™
The community-level physiological profiling (CLPP) technique was performed using Biolog ® EcoPlates™ System (Biolog Inc., Hayward, CA, USA). The catabolic fingerprinting of soil microbial communities was evaluated with the use of 31 different carbon sources [51] grouped into five categories: amines and amides, amino acids, carbohydrates, carboxylic acids, polymers. 1 g of each soil sample was transferred into conical flasks holding 99 mL of sterile 0.9% NaCl, vortexed (30 min, 150 rpm, 25 • C) and cooled (30 min, 4 • C) [9]. Finally, each of the solutions were transferred into each of the wells in the EcoPlate and incubated (144 h, 28 • C).
As the most intensive metabolic activity was noted after 120 h of incubation, the results obtained in mentioned time are presented in this study. The results reading was performed every 24 h (MicroStation ID System at 590 nm). The results are shown as a percentage of utilization of individual compounds. Clustering analysis tree was prepared from all results obtained in 120 h of analysis after correction with negative control using Past 3.25 software (Oslo, Norway) [52]. Graphical presentation of the obtained results was prepared with use of Interactive Tree of Life Platform Software 5.6.2 [53]. Moreover, results of CLPP are also expressed as average well-color development (AWCD), the Shannon-Wiener (H') and Evenness (E) [54]. Graphical presentation of functional results was carried out with the use of R Studio packages (Northern Ave, Boston, MA, USA) and GraphPad Prism 8 (San Diego, CA, USA).

Biodiversity of the Studied Technosols-Phyla Taxonomic Level
Main phyla of bacteria identified by culture independent technique in the three studied soil profiles with different chemical properties (R1, R2, R3) being an effect of successful (R1 and R3 profiles) and unsuccessful (R2 profile) reclamation are shown in Figure 1. of all identified sequences) in R2-C3 soil horizon. Quite high abundance of Proteobacteria (40.49%) was noted in R2-C2 soil horizon, whereas in remaining horizons of the R2 profile oscillated at the  It was evidenced that Actinobacteria and Proteobacteria are dominants in each of the studied profiles, however their abundance was dependent on the soil depth and soil chemical properties being an effect of reclamation process. Proteobacteria were found at the highest abundance (96.94% of all Agronomy 2020, 10, 1795 7 of 22 identified sequences) in R2-C3 soil horizon. Quite high abundance of Proteobacteria (40.49%) was noted in R2-C2 soil horizon, whereas in remaining horizons of the R2 profile oscillated at the range of 14.55-25.20%, for R2-AC and R2-C1, respectively. Their lowest abundance occurred in the deepest part of R2 profile (6.13%). A predomination of Proteobacteria (49.69-61.61%) was also found in the deeper horizons of R3 profile (R3-C2 and R3-C1, respectively) ( Figure 1). In the top layers of R3 profile Proteobacteria abundance ranged from 22.83% (R3-AC) to 29.15% (R3-A). The lowest Proteobacteria share (15.44-30.18%) was found in R1 profile.Acidobacteria and Chloroflexi phyla were determined as subdominants in Technosols bacterial structure ( Figure 1). Their numbering was sensitive to soils properties (being an effect of reclamation) and the depth in the soil profile. They found the most accurate environment for colonization in the strongly acidic R2 profile, because the highest abundance of Acidobacteria (14.84-19.59%) in the surface layers of R2 was found. In R2-C2 their number amounted 6.92% whilst in the two deepest parts of R2 was strongly reduced to 0.05-1.36% ( Figure 1). In R1 profile Acidobacteria abundance was on the level of 10.48-10.60% in the two surface horizons and it decreased with depth to 0.68-6.80% in the subsoil. In R3 profile with neutral topsoil and acidic subsoil, Acidobacteria abundance reached maximal level (10.09%) in R3-A, followed by 9.63% in R3-C1 and 5.96% in R3-C2.
In contrary, the lowest (3.73%) Acidobacteria abundance occurred in R3-AC, where, in turn, a high number of Chloroflexi (18.73%) was recorded ( Figure 1). In the remaining parts of R3 profile, their numbers remained within the limits of 5.63-1.44%, displaying a clear decreasing trend with depth. The abundance of Chloroflexi in the R1 profile was opposite to that described above, as increased with depth (5.92-14.15%), reaching a maximum (20.94%) in the deepest part of R1 profile ( Figure 1). Finally, in R2 profile Chloroflexi preferred to inhabit the surface part (10.22%), whilst in the other horizons their abundance was strongly reduced to 0.27-4.75% in R2-C1 and R2-C4, respectively. The lowest Chloroflexi abundance was noted in R2-C2 and R2-C3 horizons, where amounted 0.003% and 0.026%, respectively.
An interesting data were obtained in respect to Verrucomicrobia, Bacteroidetes and Gemmatimonadetes phyla as they were not found at all in the R2 profile where reclamation failed and soil reaction was strongly acidic throughout the profile. In R1 profile, Verrucomicrobia relative abundance ranged from 0.22 to 6.60% whereas in R3 profile oscillated on the level of 0.009-4.46% ( Figure 1). In both cases, decrease of Verrucomicrobia numbers with depth was confirmed. Similarly, Bacteroidetes abundance displayed a reduction of growth with depth and ranged between 0.41-4.32% and 0.08-5.40% in R1 and R3, respectively. Decrease of Gemmatimonadetes number with the soil depth was noted only in R3 profile (0.002-4.65%), meanwhile in R1 the highest relative abundance of mentioned phylum in R1-C2 layer (6.00%) was evidenced whilst remained on a level of 2.73-3.54% in the surface parts of that profile.
The highest abundance of Plancomycetes was noted in R1 profile (0.85-6.51%), followed by R3 (0.15-2.78%). Against this background, the R2 profile stands out again, as the presence of Planctomycetes on the level of 2.83% was recorded in the surface layer (R2-AC), while in the remaining parts of the profile their numbers dropped strongly (0.006-0.63%) and in the horizon R2-C3 their presence was not noted at all. R2 profile differ from the two others also in respect to Firmicutes abundance. Firmicutes definitely dominated in (74.90%) the deepest part of R2 profile (R2-C4) whilst in other analysed horizons their abundance ranged from 0.08 to 7.75% ( Figure 1). Similar trend, showing preference to colonize the deepest parts of the soil profile by Firmicutes was confirmed in R1 and R3 profiles, where 14.15% and 1.03% of all identified sequences in layer R1-2C and R3-C2, respectively belonged to Firmicutes. The share of other bacterial phyla (e.g., Nitrosospirae, Chlorobi, Euryarcheota) illustrated in Figure 1 was negligible and usually did not exceed 2% of identified sequences.

Biodiversity of the Studied Technosols-Classes Taxonomic Level
Main classes of bacteria identified by culture independent technique in the three studied soil profiles (R1, R2 and R3) are presented in Figure 2.  At this taxonomic level, a large variation in bacterial structure was found depending on soil properties being the final effect of reclamation process. Moreover, biodiversity decreasing trend with soil depth was observed. It was the most visible in the deeper layers of R2 and R3 profiles ( Figure 2). In the case of strongly acidic R2 profile representing an area where reclamation process failed, in R2-C3 horizon Gammaproteobacteria were found as dominants (96.14%) whereas in R2-C4 Clostridia (74.35%) was a predominant taxon. In R3 profile, the deepest horizons (R3-C1 and R3-C2) were dominated by Alphaproteobacteria (55.84% and 45.44%, respectively), followed by Acidimicrobiia (14.51-24.90%), Actinobacteria (3.80-12.85%), Acidobacteria-6 (5.88-9.24%), and Gammaproteobacteria (3.80-4.95%). In R1 profile, domination of Actinobacteria (6.34-20.70%) was determined followed by Thermoleophilia (12.73-16.14%) and Alphaproteobacteria (4.65-10.54%). The bacterial structure on the class level in the deepest studied R1 horizon (R1-2C) differed most from the other layers ( Figure 2). In that horizon, Gammaproteobacteria and Ellin 6529 classes achieved the highest abundance amounting 15.07% and 12.87%, respectively.

Beta-Diversity in the Studied Technosols-Genera Taxonomic Level
A network-like Venn-diagrams, illustrating bacterial preferences for colonization of the various parts of soil profiles (surface, subsurface, subsoil) is shown in Figure S2 (Supplementary material). The division of the studied soils into these three zones was explained in Section 2.3. Not surprisingly, the majority of bacteria were found in the surface horizons and their number decreased with depth of the soil profile. The shared genera were mostly noted between surface and subsurface parts of technogenic soils meanwhile similarity in bacterial structure was reduced between surface and subsoil layers ( Figure S2).
The above observations were also confirmed by Venn diagrams (Figure 3), presenting the numbers of unique and shared core microbiome in the surface, subsurface and subsoil of studied profiles. If one considers surface horizons, it was demonstrated that the highest bacterial genera abundance occurred in the surface layer of profile (R3), followed by the surface layer of profile (R1) (both with favourable chemical properties), whereas the lowest abundance was noted in the surface layer of R2 profile with unfavourable chemical properties ( Figure 3).
A slightly different trend was noticed in relation to the subsurface layers, where by far the largest number of bacterial genera representatives (491) was found in the R1 profile, followed by R3 (55) and R2 (22) profiles. Also, in respect to subsoil layer, those from R1 and R2 profiles seemed to be the most willingly inhabited (462 and 148 of bacterial genera, respectively) whilst the deepest part of profile R3 was not an optimal niche for microorganisms to live (presence of 27 genera representatives were confirmed).

Beta-Diversity in the Studied Technosols-Genera Taxonomic Level
A network-like Venn-diagrams, illustrating bacterial preferences for colonization of the various parts of soil profiles (surface, subsurface, subsoil) is shown in Figure S2 (Supplementary material). The division of the studied soils into these three zones was explained in Section 2.3. Not surprisingly, the majority of bacteria were found in the surface horizons and their number decreased with depth of the soil profile. The shared genera were mostly noted between surface and subsurface parts of technogenic soils meanwhile similarity in bacterial structure was reduced between surface and subsoil layers ( Figure S2).
The above observations were also confirmed by Venn diagrams (Figure 3), presenting the numbers of unique and shared core microbiome in the surface, subsurface and subsoil of studied profiles. If one considers surface horizons, it was demonstrated that the highest bacterial genera abundance occurred in the surface layer of profile (R3), followed by the surface layer of profile (R1) (both with favourable chemical properties), whereas the lowest abundance was noted in the surface layer of R2 profile with unfavourable chemical properties (Figure 3). A slightly different trend was noticed in relation to the subsurface layers, where by far the largest number of bacterial genera representatives (491) was found in the R1 profile, followed by R3 (55) and R2 (22)profiles. Also, in respect to subsoil layer, those from R1 and R2 profiles seemed to be the most willingly inhabited (462 and 148 of bacterial genera, respectively) whilst the deepest part of profile R3 was not an optimal niche for microorganisms to live (presence of 27 genera representatives were confirmed).  In general, NGS resulted in identification of 7747 bacterial genera in all (13) studied soil samples (Table S3, Supplementary material). The top bacterial genera with abundance above 10% of all sequences were selected and summarized in Table 1.
Acidiphilium and Gaiella representatives were present both in the surface, subsurface and subsoil horizons of technogenic soils, however their presence depended on soil chemical properties being the final reclamation final effect (Table 1; Table S3, Supplementary material). Gaiella was found in R1 profile, whilst Acidiphiliumin R2 profile surface and subsurface part and R3 profile subsurface and subsoil zones. Acidothermus and Metallibacterium presence were confirmed in R2 profile surface and subsurface soil layers. Sulfurifustis inherence in the abundance >10% of sequences was noted in R1 profile (surface and subsoil part). The subsurface and subsoil parts of strongly acidic R2 profile seemed to be an adequate niche for Acidithiobacillus and Acidocella colonization.  In general, NGS resulted in identification of 7747 bacterial genera in all (13) studied soil samples (Table S3, Supplementary material). The top bacterial genera with abundance above 10% of all sequences were selected and summarized in Table 1.
Acidiphilium and Gaiella representatives were present both in the surface, subsurface and subsoil horizons of technogenic soils, however their presence depended on soil chemical properties being the final reclamation final effect (Table 1; Table S3, Supplementary material). Gaiella was found in R1 profile, whilst Acidiphiliumin R2 profile surface and subsurface part and R3 profile subsurface and subsoil zones. Acidothermus and Metallibacterium presence were confirmed in R2 profile surface and subsurface soil layers. Sulfurifustis inherence in the abundance >10% of sequences was noted in R1 profile (surface and subsoil part). The subsurface and subsoil parts of strongly acidic R2 profile seemed to be an adequate niche for Acidithiobacillus and Acidocella colonization. Remaining top bacterial genera with abundance above 10% of all sequences included in Table 1 were characteristic for either surface, subsurface or subsoil parts of the studied Technosols. In the surface layers of R2 profile presence of Conexibacter, Granulicella and Candidatus Udaerobacter were confirmed, meanwhile Nocardioides and Streptomyces were identified in the surface parts of R1 and R3. The subsurface layers of R2 and R1 profiles were dominated by Leptospirillum and Kribbella, respectively ( Table 1). The relatively high biodiversity occurred also in the deepest part (subsoil) of the technogenic soils. Pseudonocardia, Jatrophihabitants, Alakanibacter and Ferrithrix seemed to be characteristic genera for R3 profile (Table 1). Bacillus, Paenisporosarcina, Psychrobacter and Serratia presence were noted in the deepest part of R1 profile, whereas Desulfosporosinus and Cellulomonas were able to colonize the subsoil layer of R2 profile. The above observations were also confirmed by the calculated values of the biodiversity indices (Table 2).  The C3 horizon of R2 profile (sampleR2-C3) was characterized by the lowest diversity (H = 0.335). Simpson index of dominance (1/D) expressing 'appreciable' the number of species amounted to 0.8635-0.9927, 0.6206-0.9924 and 0.0881-0.9626 in R3, R1 and R2 profile, respectively (Table 2).
Such a trend undoubtedly indicates a direct relationship between biodiversity recovery and soil chemical properties being the effect of reclamation process. It should be pointed that in fact R3 and R1 soil profiles were characterized by a higher number of bacteria, represented in a smaller number (<10% of the sequences, Table S3, Supplementary material), whereas biodiversity in R2 was mostly restricted to the presence of the genera that were predominant (>10%), shown in Table 1.

Bacterial Community-Level Physiological Profiles-Biolog ® EcoPlates™
The Biolog ® EcoPlates™ technique was applied to characterize metabolic changes in soil communities. The calculated Shannon-Wiener index in respect to functional diversity (Table 3) confirmed that both the highest biodiversity and metabolic activity occurred in the surface parts of the soil profiles, where H amounted 3.219-3.286, 2.713-2.836 and 3.210-3.288 for R1, R2 and R3 profiles, respectively. The H decreased with soil depth to the values of: 2.736 (R3 profile), 2.636 (R1 profile) and 1.875 (R2 profile). The calculated evenness (E) index was maintained at a level of 0.983-1.248 for R1 profile, and 0.980-1.338 for R3 profile. In respect to R2, E was detected only in the surface layer (1.066), what undoubtedly testified about unfavourable soil chemical properties (strong acidity) due to failed reclamation effect.
Effect of phenotyping analysis based on results of CLPP technique was presented in Figure 4. The highest average well colour development (AWCD) index was calculated for surface horizon of R3 profile (1.497) and R1 profile (1.329), whereas the lowest was in the surface horizon of R2 profile (0.314). This testify that the highest metabolic levels occurred in the surface parts of technogenic soils profiles, where also biodiversity was the highest. AWCD displayed decreasing trend with soil depth and it was low in deepest parts of each soil profile (AWCD amounted 0.041; 0.048 and 0.048 in R1-2C, R2-C4 and R3-C2 horizons). reclamation and (2) the soil depth.
The samples with the most and the least effective carbon sources utilization were R3-A (median = 1.587, 95%CI = 1.292/1.832) and R2-C4 (median = 0.000, 95%CI = 0.000/0.012), respectively ( Figure  S3, Supplementary material). The results showed that metabolic activity determined by CLPP technique could be a good indicator of soil properties showing the reclamation process effectiveness in the studied Technosols.  It was found that carbohydrates, carboxylic acids, and amino acids were the most preferentially utilized by soil bacteria inhabiting technogenic soils. Polymers seemed to be less preferable carbon source meanwhile usage of amines and amides was the lowest (Figure 4) (Figure 4). The mean utilization of the substrates from polymer group was the highest in R2 profile (5.81-22.78%), followed by R3 profile (10.69-18.54%) and the lowest in R1 profile (10.14-17.42%). The weakest utilization was stated for amines and amides as in general their usage did not exceed 7% and in majority of samples remained on the level of 4.3-6.3%.
R3-A and R1-C1 were the samples in which the most effective substrate utilization was confirmed. Cell growth was observed in 96.66% of wells representing carbohydrates and in 100% of wells representing the remaining substrates. R2-C4 was the least effective sample which indicated cell growth in 53.33% of wells representing carbohydrates, 44.44%-carboxylic acids and amino acids, and 41.66%-polymers. Moreover, the R2-C4 sample did not indicate the cell growth in any of wells representing amines and amides.
Based on the results obtained with the Biolog ® EcoPlate™ method, a heat map was constructed ( Figure 5). As shown below, the heat map also facilitated identification of differences in the preferences of the microorganisms for the utilization of the carbon sources. The highest metabolic activity was observed in the soils from the sample R3-A, whereas a slightly lower catabolic activity was noted in the R2-C4 and R2-C1 ( Figure 5). Agronomy 2020, 10, x 14 of 22

Discussion
We previously proved that an area of the abandoned mine in Rudki seems to be a perfect place to examine the relationship between soil chemical characteristics, microbiological activity and culturable bacteria diversity in Technosols after land reclamation conducted in 1970s [5]. Herein, we extended our microbiological analyses and determined the biodiversity of the bacteria using a culture independent technique to which we added a picture of their metabolic activity using the community level physiological profiling (CLPP) technique. The present study gives a comprehensive insight into Technosols biodiversity and catabolic activity and is a part of contemporary study trend regarding microbiological properties of Technosols [55][56][57][58]. Physical and chemical features of the studied soil profiles were described in detail in our previous papers [2][3][4][5]39], however a short summary of soil properties was shown in "Materials and methods" section and in Table S1 (Supplementary material) in order to allow for a better understanding of further discussion of soil biodiversity.
In the present study, we hypothesized that both culture-independent bacterial biodiversity and metabolic activity is related with soil properties and soil depth, as well as could be an adequate indirect indicators of reclamation process effectiveness in technogenic soils. The obtained results seem to confirm the validity of this hypothesis. As previous reports on Technosols microbiological properties were limited to determination of their enzymatic activities, microbial biomass [16][17][18]56,58] and culturable autochthonous microbiota structure [5,19,20,39], herein we attempted to determine the biodiversity of these soils with the use of culture-independent technique. By using NGS tools, biodiversity of bacteria at taxonomic levels of phyla, classes and genera in three different profiles of Technosols was recognized. Hafeez et al. [57] indicated that at the phyla and classes taxonomic levels microbial community structure in constructed Technosols remained in very close similarity to "natural" soils. Starting from the highest taxonomic level, the dominance of Proteobacteria and Actinobacteria was demonstrated in the studied technogenic soils (Figure 1). Other study reported that Proteobacteria dominance in Technosols were supported by: Bacteroidetes, Firmicutes, Chloroflexi, and Actinobacteria presence [57]. Proteobacteria, Acidobacteria and Bacteroidetes were found as dominant phyla in a Technosols historically contaminated by the Nevertheless, these differences depended on the type of metabolized substrate and the effect of reclamation. The i-Erythritol, alpha-Ketobutyric Acid, Itaconic Acid, Glycyl-L-Glutamic Acid, D-Galactonic Acid gamma-Lactone, L-Serine and alpha-Cyclodextrin were characterized by the highest utilization rate in the soil of R3-A. Also the beta-Methyl-DGlucoside, Glucose-1-Phosphate, D-Xylose, Phenylethylamine, D-Malic Acid and Glycogen were characterized by the highest utilization rate in the soil of R1-A ( Figure 5). We proved that preferences for particular carbon sources utilization by the soil bacterial communities depended on (1) the soil chemical properties being an effect of reclamation and (2) the soil depth.
The samples with the most and the least effective carbon sources utilization were R3-A (median = 1.587, 95%CI = 1.292/1.832) and R2-C4 (median = 0.000, 95%CI = 0.000/0.012), respectively ( Figure S3, Supplementary material). The results showed that metabolic activity determined by CLPP technique could be a good indicator of soil properties showing the reclamation process effectiveness in the studied Technosols.

Discussion
We previously proved that an area of the abandoned mine in Rudki seems to be a perfect place to examine the relationship between soil chemical characteristics, microbiological activity and culturable bacteria diversity in Technosols after land reclamation conducted in 1970s [5]. Herein, we extended our microbiological analyses and determined the biodiversity of the bacteria using a culture independent technique to which we added a picture of their metabolic activity using the community level physiological profiling (CLPP) technique. The present study gives a comprehensive insight into Technosols biodiversity and catabolic activity and is a part of contemporary study trend regarding microbiological properties of Technosols [55][56][57][58]. Physical and chemical features of the studied soil profiles were described in detail in our previous papers [2][3][4][5]39], however a short summary of soil properties was shown in "Materials and methods" section and in Table S1 (Supplementary material) in order to allow for a better understanding of further discussion of soil biodiversity.
In the present study, we hypothesized that both culture-independent bacterial biodiversity and metabolic activity is related with soil properties and soil depth, as well as could be an adequate indirect indicators of reclamation process effectiveness in technogenic soils. The obtained results seem to confirm the validity of this hypothesis. As previous reports on Technosols microbiological properties were limited to determination of their enzymatic activities, microbial biomass [16][17][18]56,58] and culturable autochthonous microbiota structure [5,19,20,39], herein we attempted to determine the biodiversity of these soils with the use of culture-independent technique. By using NGS tools, biodiversity of bacteria at taxonomic levels of phyla, classes and genera in three different profiles of Technosols was recognized. Hafeez et al. [57] indicated that at the phyla and classes taxonomic levels microbial community structure in constructed Technosols remained in very close similarity to "natural" soils. Starting from the highest taxonomic level, the dominance of Proteobacteria and Actinobacteria was demonstrated in the studied technogenic soils (Figure 1). Other study reported that Proteobacteria dominance in Technosols were supported by: Bacteroidetes, Firmicutes, Chloroflexi, and Actinobacteria presence [57]. Proteobacteria, Acidobacteria and Bacteroidetes were found as dominant phyla in a Technosols historically contaminated by the burning of chemical ammunitions [56]. The fact that Proteobacteria dominated in the soil environment is well known [6,8,9,56,57], however current results evidenced that Proteobacteria are able to colonize not only the topsoil but also the subsoil, what remains in agreement with other studies [59,60]. However, the fact that these bacteria possess ability for adaptation to unfavourable environmental conditions (strong acidity), as in the whole R2 profile or in deeper part of R3 profile ( Figure 1) is a new finding. On the other hand, deeper parts of R2 profile (C3 and C4 horizons) with strong acidity and high salinity are microbiologically interested because domination of both Proteobacteria and Planctomycetes was stated (Figure 1). This shows the high resistance of these phyla to unfavourable environmental conditions and the ability to live in soil with low pH, high salinity, and low concentration of soil organic matter. Members of the bacterial phylum Planctomycetes are known to be able to inhabit a wide range of aquatic and terrestrial environments with diverse environmental conditions, i.e., polar desert soils at the northern land limit of the arctic polar region [61] and low-temperature environments (i.e., boreal wetlands) [62,63]. Also it was suggested that Planctomycetes participate in degradation of plant-derived polymers, exoskeletons of peat-inhabiting arthropods, as well as exopolysaccharides produced by other bacteria [63]. However, till now there was no information about their abundant occurrence in strongly acidic Technosols developed from mine wastes containing Fe sulphides. That finding suggests that Planctomyctes ( Figure 1) could inhabit soil having extreme properties even after reclamation process failure. In contrary, Actinobacteria abundance could be recommended as potential indicator of effectiveness of reclamation, as they were present in higher abundance in R1 and R3 profiles rather than in R2, where their presence was limited to surface horizons. Very interesting findings were noted in respect to Verrucomicrobia, Bacteroidetes and Gemmatimonadetes phyla (Figure 1) whose presence was not demonstrated at all in strongly acidic R2 profile. This proves that mentioned phyla can be treated as indirect indicators of the reclamation process success, whilst the recorded lack of their presence in the profile representing an area where reclamation failed confirms the ineffectiveness of the reclamation attempt.
On the classes taxonomic level, Gammaproteobacteria and Clostridia could be undoubtedly treated as indicators of reclamation ineffectiveness in Technosols developed from mine wastes containing Fe sulphides, as both group dominated in the deepest parts of R2 profile (Figure 2). The ability for colonizing acidic soils by Gammaproteobacteria and Clostridia were reported previously [64,65], however their great abundance was generally noted in the surface soil layers [56,64,65], whereas in the present study we found that this bacteria representatives are also able to colonize the subsoil parts. Moreover, Gammaproteobacteria dominance was confirmed in Technosols constructed for remediation of a contaminated industrial wasteland [57] as well as in Technosols contaminated with arsenic [56].
Our study at genera taxonomic level showed that, if one considers surface soil horizons, the highest bacterial abundance occurred in surface horizons of R3 and R1 profiles where reclamation processes improved soil properties, whereas the lowest abundance was noted in the surface horizons of R2 profile representing an area where reclamation process failed. This trend was confirmed by biodiversity indices (Table 2), especially by H index that reached the highest values in R3 profile, then in surface part of R1 profile, whilst the lowest H value was noted in R2 profile. In general, higher biodiversity represented by genera at lower abundance <10% was typical of R3 and R1 profiles having favourable and partly favourable properties and representing successfully reclaimed areas, whereas in strongly acidic R2 profile where reclamation failed, mostly genera of large numbers (>10% sequences, Table 1) were present. All of noted genera in R2 profile (Conexibacter, Acidothermus, Granulicella, Metallibacterium, Acidithiobacillus, Acidocella, Leptospirillum, Desulfosporosinus, Cellulomonas) could be treated as potential indicators testifying about inefficiency of reclamation process in Technosols developed from mine wastes containing Fe sulphides. In majority, they are represented by genera preferring acidic conditions as their optimal niche for colonization and are classified as sulphur oxidizing bacteria [66,67]. Therefore, their presence in such extreme environment as strongly acidic technogenic soils containing Fe sulphides is justified. Leptospirillum presence was also noted in Technosols contaminated by the burning of chemical ammunitions [56]. The success of reclamation processes can be proved by the abundance (>10%) of genera: Nocardioides, Streptomyces, Pseudonocardia, Jatrophihabitans, Alkanibacter, Ferrithrix (Table 1). Actinobacteria is a phylum and class of Gram-positive bacteria that are ubiquitous in nature, also in the soil environments [68]. Their presence in soil is usually beneficial as they can: (1) be utilized as biofertilizers for sustainable agriculture, (2) enhance plant growth and soil health though different plant growth promoting attributes, i.e., solubilization of phosphorus, potassium and zinc, production of Fe-chelating compounds and phytohormones, as well as by biological nitrogen fixation [68]. In addition to these important ecological functions, recommendation of Actinobacteria representatives (Nocardioides, Streptomyces, Pseudonocardia, Jatrophihabitans) as indicators of reclamation success in areas where Fe sulphides occur in anthropogenic soil substrates is a new approach. Alkanibacter was noted to be present in soil contaminated with petroleum substances [69] or tetracycline [70], whereas Ferritrix of Firmicutes phylum is mainly known as acidophilic, autotrophic and Fe/S oxidizing bacteria with ability to facilitate bioremediation processes [71] or as microbial methanol-utilizer [72].
Community-level physiological profiling (CLPP) analyses from very diverse environments are frequently used in order to characterize the metabolic versatility of the whole environmental bacterial communities [37]. Despite advantages of CLPP that include the simplicity of the protocol and the largely reduced cost of analysis, many limitations in the use of this approach for complex environmental samples have been previously reported (i.e., preference of fast-growing bacteria in the assay, the data analysis and the interpretation of the CLPP results [73][74][75][76]. However, CLPP is still commonly used to select microorganisms that are able to utilize various carbon compounds and that can be subsequently identified using molecular methods [37,76]. Due to the fact that every bacterial community has a characteristic reaction pattern with different optical density values for different carbon compounds, a specific pattern called a 'metabolic fingerprint' is received as an effect of CLPP assay [76,77]. Herein, we applied CLPP, together with NGS, in order to prove that CLPP could be extremely effective in recognizing the influence of soil properties being an effect of land reclamation on microbiological properties of Technosols containing Fe sulphides. Utilization of all substrates in Biolog ® EcoPlates™ analysis in studied samples ( Figure S3, Supplementary material) clearly demonstrated the differences in metabolic activity of microorganisms which coincided with soil chemical properties (Table S2, Supplementary material) being an effect of reclamation processes, as well as with soil depth. Most of the carbon substrates were more intensively utilized in R3 profile, followed by the surface layers of R1, where soil chemical properties (pH, content of carbonates, TOC, and TN) were favourable due to successful reclamation. We found that carbohydrates and carboxylic acids were the most preferentially utilized by soil bacteria inhabiting technogenic soils. Our results may be supported by the findings of Kenarova et al. [78] and Martinez-Toledo et al. [79] who indicated that heavy metals contaminated soils generate environmental stress, demanding a high energy generation from microorganisms. Consequently, carbon sources with a high energy input, i.e., carbohydrates and carboxylic acids are firstly metabolized [78,79] in contaminated environments. Moreover, it was evidenced that D-Galactonic Acid gamma-Lactone, found in the current study as substrate with the highest utilization in R3, belongs to carbon source that is the most used to induce resistance to heavy metals contamination in microbial communities [78]. Studies on engine-oil contaminated soils have shown that the use of different remediation methods increases the activity of microorganisms and functional richness [80]. It has also been demonstrated that the addition of heavy metals to the sludge affects the decrease of biological activity in the sample [81]. Very low catabolic activity was typical of strongly acidic R2 profile representing an area where reclamation failed. Mentioned trends were reflected by calculated AWCD values, reaching the highest levels in the surface parts of technogenic soils (R3 > R1 > R2). A negative correlation between the AWCD and heavy metals presence in soil environment was noted in other studies [82,83]. Research on Technosols in comparison to forest soil [84] displayed also lower AWCD. The fact that CLPP is a suitable assay for estimation of reclamation results was also confirmed by calculation of Shannon E index (Table 3), focused on the evenness of values across all utilized substrates [76] which remained uncountable (not detectable) in majority of horizons of R2 profile. However, it is worth mentioning that CLPP does not reflect the functional abilities of the entire soil microbial community but that of limited subset of microbial fraction capable to grow on the carbon sources provided in the Biolog ® EcoPlates™ [85]. Therefore, in order to receive a comprehensive insight into microbial structure and activity in the heterogeneous soil environment such as Technosols the combination of two techniques (CLPP and NGS) is justified [9,32,86].

Conclusions
The combination of two techniques: NGS and CLPP resulted in comprehensive recognition of bacterial biodiversity and their catabolic activity in three Technosols developed from mine wastes containing Fe sulphides in Poland. All soils had different chemical properties being an effect of reclamation processes conducted in 1970s. It was evidenced that changes in microbial structure were the consequence of a combined effect of (1) the soil chemical properties and (2) the soil depth.
The following taxa (from higher to lower taxonomic level) were found in soils representing areas where soil chemical properties were improved by reclamation process: Verrucomicrobia, Bacteroidetes, Gemmatimonadetes, Alphaproteobacteria, Acidimicrobia, Nocardioides, Strepromyces, Pseudonocardia, Jatrophihabitants, Alkanibacter and Ferrithrix. Therefore they can be treated as an indirect microbial indicators of reclamation process success in technogenic soils developed from mine wastes containing Fe sulphides. On the contrary, reclamation ineffectiveness manifested by strong acidity of Technosols of the study area could be confirmed by (1) the lack of the following phyla: Verrucomicrobia, Bacteroidetes and Gemmatimonadetes, as well as (2) the high abundance of the following bacterial representatives: Gammaproteobacteria, Clostridia, Conexibacter, Acidothermus, Granulicella, Metallibacterium, Acidithiobacillus, Acidocella, Leptospirillum, Desulfosporosinus and Cellulomonas.
The CLPP analysis allowed for recognition of microbial preferences for the intensity of carbon substrate utilization by obtaining catabolic fingerprinting. CLPP (analogically like NGS) evidenced that metabolic activity depended on (1) the chemical properties of soils being the result of reclamation process and (2) the soil depth. Most of the carbon substrates were more intensively utilized in topsoil horizons of the successfully reclaimed profile R3 and R1, whereas very low catabolic activity was characteristic for strongly acidic R2 profile where reclamation failed. The most easily utilized group of compounds were as follows: carbohydrates > carboxylic acids > amino acids > polymers > amines and amides.
Our findings indicate that both soil microbial community profiling and bacterial metabolic activity in technogenic soils developed from mine wastes containing Fe sulphides can be a suitable indirect indicators of reclamation process effectiveness of mine waste spoils from iron sulphide mines.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/11/1795/s1, Figure S1: Rarefaction curves for the studied soil samples based on the number of unique sequences (ASVs). Figure S2: Network-like Venn diagrams of the bacterial genera presence in the surface, subsurface and subsoil of