Structural Shifts in the Soil Prokaryotic Communities Marking the Podzol-Forming Process on Sand Dumps

: This work describes the microbial community structure of the continuously revegetated chronosequence of a former sand quarry, which demonstrates a unique example of nearly complete soil restoration in less than 100 years. Samples were collected at five time points (0, 3, 30, 70 years and mature soil) from the entire set of soil horizons, revealing the history of pedogenesis. Real-time PCR was applied to quantitatively describe the bacterial and archaeal communities. High-throughput sequencing of the bacterial and archaeal V4 variable region of the 16S rRNA gene was used to identify abundant microbial taxa. A beta-diversity analysis revealed that the prokaryotic community structure responded strongly to the processes of organic matter accumulation and the corresponding evolution of the soil into discrete horizons. Changes in soil microbiota in the course of soil profile evolution revealed three groups of prokaryotes, which tended to accumulate in the specific soil horizons and might be associated with the certain soil-forming processes, including plant roots growth. This research showed the heuristic potential of soil horizon profiling in microbiological studies as opposed to the formal depth-dependent separation of the soil layers. The results allowed us to trace the relationship between the structure of the soil prokaryotic community and the peculiarities of the evolution of the podzolic soil profile as well as to identify the microbial indicators and drivers of primary pedogenesis.


Introduction
Soil evolution is defined as the process of soil development from the initial stage-known as the "zero stage" of soil formation [1]-through the separation of well-defined horizons of zonal soil types (the so-called "embryonic" soil profile) to the formation of the mature soilthe final stage in the ecological succession [2,3].Thus, the soil profile-the unique packet of recognizable soil horizons-is the result of complex interactions between four soil-forming factors (biota, climate, parent material and relief) over time.The combination of five soilforming factors (time, climate, topography, parent material, living organisms) is unique to the particular ecosystem, so the rate and specific characteristics of soil development will depend on the vegetation, parent material, and climate [4][5][6][7].Although the soil microbiome is an essential link between all these factors, its relevance in soil development is still poorly Soil Syst.2024, 8, 9 2 of 25 understood.There are many examples of research where the biodiversity (mainly fungi, bacteria and protozoa and to a lesser extent archaea) was studied in different soil types as well as in various soil horizons [8][9][10][11].Many of them lack the complete description of the soil profile, dividing the soil layers according to the depth [12] and saying nothing of the dynamic changes of microbial diversity on the successive stages of soil development.
It is essential to study soil as a continuously evolving natural body that demonstrates vertical stratification features.Studies of the former soil iterations help us to match certain microorganisms with particular soil-forming processes and predict their potential impacts on soil fertility and function.However, it is often difficult to conduct such studies due to the paucity of suitable natural sites where different stages of soil succession over time can be observed.
Currently, we can only obtain glimpses of soil evolution from two types of studies.The first is devoted to the description of the microbial diversity in the soil horizons of mature soils [13][14][15][16][17][18][19][20][21].These studies have demonstrated a measurable decrease in the number of bacteria from the upper to lower soil layers [13,14] and have recorded maintenance of the microbial counts with depth [20].The archaea vs. bacteria ratio has been found to increase with depth, whereas fungal vs. bacterial abundance has the reverse tendency.Will et al. [14] provided an exhaustive description of the taxonomic diversity of bacteria inhabiting different soil horizons, finding Actinobacteria, Bacteroidetes, Cyanobacteria, Fibrobacteres, Firmicutes, Spirochaetes, Verrucomicrobia, Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria among the dominant population in the A horizon, whereas Acidobacteria, Chloroflexi, Gemmatimonadetes, Nitrospira, TM7, and WS3 were dominant in the B horizon.The diversification in prokaryotic communities with depth has also been demonstrated with the use of RNA analysis [22].
The second type of soil evolution study describes microbial succession in disturbed soil environments, including former mining sites [23][24][25][26], areas of deglaciation [27][28][29], post-fire ecosystems [30,31], and volcanogenic soils [32].The basic methodology of these studies has been to convert time into a spatial gradient, resulting in the ability to visualize continuous stages of soil evolution, which are known as chronosequences [27,33].Soil chronosequences represent successive stages of soil formation, starting from the same parent material and evolving in a specific climatic zone [21].
Mining activity has an enormous impact on natural ecosystems [34,35].An analysis of the studies on microbial successions from former mining sites revealed that many were devoted to complex and biologically unavailable substrates, such as uranium or oil deposits [36][37][38] and coal mining sites [23,39].The majority of these studies used conventional microbiological approaches rather than metagenomics methods.Anthropogenic disturbance in this case leads to complete disturbance of the land cover, and soil formation starts from the so-called "zero stage" [1].In the case of natural disturbances, for example, windthrow successions, the soil layer is not completely destroyed, and this model of chronoseries is much less informative than the first one.
Generally, in the initial stages of biological succession, soil is inhabited by pioneer prokaryotic communities that are involved in erosion of the rock parent material and the accumulation of organic substances [27,40].Consequently, they form substrates suitable for plant growth [40].Plants, in turn, promote the second stage of soil microbiome formation and the increase in the diversity of heterotrophic microorganisms [41,42].Simultaneously, the functional diversity of the soil microorganisms increases in this stage, providing stability to the entire ecosystem [43].
Studies on primary soil successions have been described in numerous publications [40,41,44,45], whereas the long-term development of soil was less investigated and usually restricted to the upper soil layers [21,46].Among the physicochemical properties of soil that influence microbial diversity, increased concentrations of iron and aluminum oxides have been noted with increasing age of the soil [47,48].Prokaryotic communities undergo substantial changes in their structure during the initial stages of succession with the rate of change decreasing thereafter [21,46,47,49].Some authors have reported a decrease in the number of bacteria as soil ages due to transformation of the organic matter from biologically available to inert forms [39].This decrease in bacteria counts was accompanied by the active growth of fungal communities in the final stages of succession [21,26].The number of archaea increased in mature soils, which can also be associated with the accumulation of specific substances and the partitioning of ecological niches [21].The maintenance of such specific ecological niches seems to be particular to lower soil horizons where organic matter forms complexes with the mineral matrix of the soil [47,50,51].
During the initial stages of soil development, the microbial community is enriched with the bacterial taxa Acidobacteria, Bacteroidetes, Verrucomicrobia, Comamonadaceae, and Betaproteobacteria [27].Terrimonas, Nocardioides, Solirubrobacter, and Microvirga dominate in mature soils [20].The rate of succession decreases in the final stages of soil development [52].The decrease in microbial biomass and increase in the number of oligotrophic groups of bacteria can largely be explained by modification of the organic matter and its transformation into inert forms [52][53][54].
The changes observed in microbial community compositions over time [27,28,45,54] are entirely dependent on changes occurring in the soil physicochemical properties, particularly pH, and the amount of organic carbon (C), phosphorus, and nitrogen (N) [45,46,55].The diversity of archaea is affected by N and plant community composition [45].After 30-40 years of succession, the values of the physicochemical parameters stabilize and are maintained until the final stage [39].
The main disadvantage of all the aforementioned studies is the lack of a comprehensive description of the soil evolutionary process, including separation of the early stages of soil genesis from the intermediate stages, and stages of soil maturation.The main obstacle to such studies is identifying sites with relatively high rates of soil evolution.In this respect, noteworthy sites are where pine forests have formed in the sand quarries located in the Leningrad region of Russia, which preserve soil chronosequences of 100 to 200 years.This paper aims to characterize the soil prokaryotic communities in the entire soil profile using metagenomics data from different stages of Podzol soil chronosequences that formed on the former spoil heaps of a sandy quarry.The objectives of this study were to (1) evaluate the main indices of soil development in the chronosequence, including the soil morphology and basic chemical and physical properties; (2) characterize the soil microbial community composition of differently aged sites using 16S rRNA high-throughput sequencing; (3) analyze the levels of microbial diversity in the different ages of soil and multiple soil horizons; and (4) identify the horizon-and stage-specific microorganisms that mark the steps of soil profile development.

Study Sites
The chronosequence of sandy soils in the Maluksa sand quarry (Kirovsky District, Leningrad Region, Russia) is one of the foremost examples of Podzol chronosequences in northwestern Russia (Figure 1).According to previous studies, this chronosequence includes soils ranging from the initial sod-sandy soils with no features of the podzolization process [3] to embryonic Podzols with evident features of Spodic profile differentiation [56].Soil organic matter dynamics, as well as the composition of humic substances, have previously been analyzed for differently aged soils in the same area [3,56].
The quarry is situated in the southern taiga zone, with an annual precipitation of ~600-700 mm and a mean annual temperature of ~3.5 • C. The quarry, ~10 km long and up to 2 km wide, is a former esker composed of glaciofluvial sands.Mining began in the quarry in 1930.Intensive mining activity was aimed at exploring the ridge (fluvioglacial esker) with an average height ranging from 7 to 10 m.The mining activity resulted in the formation of lakes and ponds, which now occupy up to 70% of the area.The soil parent material comprises homogeneous (in both chemical composition and particle size distribution) glaciofluvial sands.This parent material has promoted rapid pedogenesis, resulting in the formation of a well-developed and complete chronosequence.The vegetation of zonal mature plots Soil Syst.2024, 8, 9 are presented by Scotch pine (Pinus sylvestris) forest stands with shrub (Vaccinium myrtilus, Vaccinium vitis-iadae) and mosses soil cover.Plants of zero stage of succession are presented by sporadic reedgrass plants (Calamagrostis epigejos).At the three-year stage, sandy soils begin to be colonized by pine, alder (Alnus incana) and willow (Salix caprea), while the projective cover of vegetation is 20-30%.Then, the projective coverage increases, alder and willow are replaced by pine, sporadic birch trees (Betula pendula) appear, and a shrub layer is formed.
formation of lakes and ponds, which now occupy up to 70% of the area.The soil parent material comprises homogeneous (in both chemical composition and particle size distribution) glaciofluvial sands.This parent material has promoted rapid pedogenesis, resulting in the formation of a well-developed and complete chronosequence.The vegetation of zonal mature plots are presented by Scotch pine (Pinus sylvestris) forest stands with shrub (Vaccinium myrtilus, Vaccinium vitis-iadae) and mosses soil cover.Plants of zero stage of succession are presented by sporadic reedgrass plants (Calamagrostis epigejos).At the three-year stage, sandy soils begin to be colonized by pine, alder (Alnus incana) and willow (Salix caprea), while the projective cover of vegetation is 20-30%.Then, the projective coverage increases, alder and willow are replaced by pine, sporadic birch trees (Betula pendula) appear, and a shrub layer is formed.Sampling was performed in May 2018, at five morphologically dissimilar sites, representing sequential stages in the soil development: M1-the initial (zero) stage, lacking any evidence of vegetation; M2-sandy 3-year-old primary sod gray humus soil with developed humus layer and first features of initial parent material alteration; M3 and M4-30-and 70-year-old Embryopodzols, respectively; and, finally, M5-mature (benchmark) Podzol beneath mature pine forest.The Embryopodzols were deposited on the flat surfaces of the former spoil heaps that had been naturally revegetated by Scots pine (Pinus sylvestris L.) after the quarrying was ended.Samples were collected from the middle of each morphologically defined genetic horizon.Oi-forest litter with initial features of organic material transformation, Oe-forest litter with medium stage of organic matter transformation, Ay-horizon of humus accumulation, E-eluvial podzol horizon, BFhorizon of iron oxide eluviation, C-initial parent material.
The soil profile morphology (Figure 2) showed that podzolization started at the 15year-old site with the formation of an illuviated iron horizon.In the Russian classification system, this soil can be characterized as an Entic Podzol due to the absence of well-separated eluvial horizons in the vertical profile.The 30-year-old site included a poorly developed E horizon (brackets and a lowercase letter usually mark this horizon).The more mature soils were characterized by a thickening of the profile, particularly in the E and B Sampling was performed in May 2018, at five morphologically dissimilar sites, representing sequential stages in the soil development: M1-the initial (zero) stage, lacking any evidence of vegetation; M2-sandy 3-year-old primary sod gray humus soil with developed humus layer and first features of initial parent material alteration; M3 and M4-30-and 70-year-old Embryopodzols, respectively; and, finally, M5-mature (benchmark) Podzol beneath mature pine forest.The Embryopodzols were deposited on the flat surfaces of the former spoil heaps that had been naturally revegetated by Scots pine (Pinus sylvestris L.) after the quarrying was ended.Samples were collected from the middle of each morphologically defined genetic horizon.Oi-forest litter with initial features of organic material transformation, Oe-forest litter with medium stage of organic matter transformation, Ay-horizon of humus accumulation, E-eluvial podzol horizon, BF-horizon of iron oxide eluviation, C-initial parent material.
The soil profile morphology (Figure 2) showed that podzolization started at the 15-year-old site with the formation of an illuviated iron horizon.In the Russian classification system, this soil can be characterized as an Entic Podzol due to the absence of well-separated eluvial horizons in the vertical profile.The 30-year-old site included a poorly developed E horizon (brackets and a lowercase letter usually mark this horizon).The more mature soils were characterized by a thickening of the profile, particularly in the E and B horizons, and differentiation of the upper soil horizon (O) into accumulation (Oi) and fermentation (Oe) portions.Podzol formation is one of the fastest processes of soil formation, especially on mineralogically poor sands of fluviglacial origin.The differentiation of soil profile by Al-Fe-humus type occurs rapidly, and preclimax soil can be formed in the first hundred years after the beginning of soil formation [3,4].
Each sample was taken in three randomly chosen locations.Subsamples were pooled to create composite samples.Composite samples were taken both for agro-chemical analyses (3 replicates, each containing 200 g of soil, which were immediately transported and analyzed in the laboratory) and soil description and DNA analysis (5 replicates per sample, which were transported to the laboratory in 10 mL sterile tubes and stored at −70 °C).The soils were classified as Podzol Spodosol, which was based on the 30-year-old stage (World Reference Base 2015).

Soil Analyses
Organic and mineral soil materials were ground and sieved through 2 mm meshes, and large root debris was removed manually.A CHN analyzer (model Euro EA3028-НТ) For the further analysis (physicochemical analysis and DNA extraction), the soil was sampled at each horizon of the five sites, according to the following scheme: M1-the initial (zero) stage, lacking any evidence of vegetation (samples of C1-C2-C3-C4); M2-sandy 3-year-old primitive sod gray humus soil (samples of Oi-Ay-BF-C); M3 (samples of Oi-E-BF-BC-C) and M4 (Oi-Oe-E-BHF-BF-C)-30-and 70-year-old Embryopodzols, respectively; and, finally, M5-mature (benchmark) (Oi-E-BF-C).
Each sample was taken in three randomly chosen locations.Subsamples were pooled to create composite samples.Composite samples were taken both for agro-chemical analyses (3 replicates, each containing 200 g of soil, which were immediately transported and analyzed in the laboratory) and soil description and DNA analysis (5 replicates per sample, which were transported to the laboratory in 10 mL sterile tubes and stored at −70 • C).The soils were classified as Podzol Spodosol, which was based on the 30-year-old stage (World Reference Base 2015).

Soil Analyses
Organic and mineral soil materials were ground and sieved through 2 mm meshes, and large root debris was removed manually.A CHN analyzer (model Euro EA3028-HT) was used to determine the C and N contents with the data being converted to the air-dried form.The pH values were determined using a routine method in water and calcium chloride suspensions.The particle size distribution and texture class were determined using a Soil Syst.2024, 8, 9 6 of 25 sedimentometric method with pyrophosphate peptization of the microaggregates [57].The C from the microbial biomass (C mic ) was quantitatively determined according to the formula C mic (g C/g soil) = (emitted CO 2 /g soil per hour) × 40.04 + 0.37.The microbial metabolic rate (the specific respiration of the microbial biomass, qCO2) was determined as the ratio of the basal respiration (BR) to the microbial biomass carbon index: qCO2 (mg CO 2 /C mic mg/h) = BR/C mic [58].ANOVA analysis was performed to evaluate the significant differences among samples in the main physicochemical and microbiological parameters.

Soil DNA Extraction
Soils for DNA extraction were not ground before analysis.The samples were frozen in the field and then transported to the laboratory.DNA was extracted from 0.2 g of soil using a PowerSoil DNA Isolation Kit (Mobio Laboratories, Solana Beach, CA, USA), which included a bead-beating step, based on the manufacturer's specifications.Homogenization of the samples was performed using a Precellys 24 (Bertin Corp, Rockville, MD, USA) at 6.5 m/s twice for 30 s.The purity and quantity of the DNA were tested by electrophoresis in 0.5× TAE buffer on 1% agarose.The quality of DNA extraction was assessed on the electrophoregram by the presence of a clear bright band (the length of the DNA fragments was approximately 20,000 bp compared to the lambda-HindIII marker DNA).The DNA concentrations were measured at 260 nm using a SPECTROStar Nano (BMG LABTECH, Ortenberg, Germany).The average DNA concentrations ranged between 30 and 50 ng/µL.
The relative abundances of bacterial and archaeal small subunit rRNA gene copies were analyzed by quantitative PCR (qPCR).The pairs of primers Eub338/Eub518 [59] and 915f/1059r [60] was used for bacteria and archaea correspondingly [61].ANOVA analysis was performed to evaluate the significance of differences between qPCR samples.

Next-Generation Sequencing (NGS) of Amplicon Libraries of the 16S rRNA Gene
The purified DNA templates were amplified with the universal multiplex primers F515 5 ′ -GTGCCAGCMGCCGCGGTAA-3 ′ and R806 5 ′ -GGACTACVSGGGTATCTAAT-3 ′ [62], which targeted the variable region V4 of the bacterial and archaeal 16S rRNA genes.Each multiplex primer contained the adapter, 4 bp key (The Centre for Applied Genomics, TCAG), 10 bp barcode, and primer sequences.The expected length of the amplification product was 300 bp.The sequencing of the amplicon libraries was conducted using an Illumina MiSeq sequencer in the Genomic Technologies, Proteomics, and Cell Biology Department (All-Russia Research Institute for Agricultural Microbiology, ARRIAM).

Bioinformatics and Statistical Analysis
The raw sequences were processed using the QIIME_1.9.1 environment [63].The raw reads were preliminarily processed using TRIMMOMATIC software [64].The multiplexed reads were first filtered for quality and then grouped according to barcode sequences to reduce sequencing errors.Sequences were omitted from the analysis if they were less than 200 bp, had a quality score of less than 25, contained uncorrectable barcodes, primers, ambiguous characters, or had a homopolymer length equal to or greater than 8 bp.All non-bacterial ribosomal sequences and chimeras were also removed from the libraries.The chimeras were removed using the chimera_slayer.pyscript, which was incorporated in QIIME.For analysis of the operational taxonomic unit (OTU) abundance fluctuations between samples, the decontam package [65] in the R environment (https: //www.R-project.org/, the version R 4.1.2,accessed on 1 November 2021) was used.
In total, 1,023,728 sequences were obtained with an average of 33,023 sequences per library.The dataset was subjected to a rarefaction procedure, reducing the sample outcome to the minimal quantity of sequences in the dataset (16,365 sequences per sample) (Figure S1).OTU picking was performed using a close-reference picking algorithm based on 97% sequence similarity with the SILVA database, version 132, 2018 (https://www.arbsilva.de/download/archive/qiime,accessed on 28 November 2019).Aligned sequences were used to build a distance matrix with a distance threshold of 0.1 and the phylogenetic tree necessary for downstream analysis.
Alpha-and beta-diversity analyses were performed to compare prokaryotic communities.The alpha diversity was determined by calculating the indices for richness (observed species, ChaoI) and evenness (Faith's index, Shannon evenness).Rarefaction analysis was applied to estimate the efficiency of the conducted sequencing effort.An ANOVA was performed to verify the observed differences.For the beta diversity, the weighted UniFrac metric [66] was used to calculate the amount of dissimilarity (distance) between the prokaryotic communities being compared.PERMANOVA was applied to calculate the significant differences between samples in beta-diversity analysis.The results were presented in a principal components analysis (PCoA) using the Emperor program implemented in QIIME.All estimates were measured for the normalized data (normalization was carried out to the smallest number of sequences present in the sample).
The core microbial community (generally defined as the specific taxonomic set presenting in each sample of the current dataset) of the entire profiles (obtained by combining all of the horizons together) in each stage of the soil chronosequence was calculated in the QIIME environment.The OTU abundances were compared among samples by calculating the median relative change values for all groups of replicates.A positive median indicated an increase in abundance, whereas a negative median could be considered evidence for a decrease in abundance.A basic permutation test was used to infer significance, while a jackknife-like resampling approach was applied to test the stability of the median estimates.An assessment of the reliability of differences in the representation of individual taxa in the analyzed samples was conducted using a custom script that performed several pairwise comparisons of the OTU frequency contingency tables in different sample sets, including the replications.The algorithm dynamically selected either the G-test or Fisher's exact test and applied the Bonferroni correction to the p-value.The dynamic behavior of the algorithm allowed it to be applied to all types of datasets, including both the normally and nonnormally distributed data types-that is, the soil microbial community NGS data [67].The 50 most abundant OTUs were chosen to illustrate the transitions in abundance of specific OTUs over time.Each OTU was heat-mapped within the graphically painted soil profiles using the phyloseq [68] and tidyverse (https://CRAN.R-project.org/package=tidyverseversion 2.0.0., the version R 4.1.2,accessed on 1 November 2021) packages in RStudio (http://www.rstudio.com/,the version 2021.09.0, accessed on 29 September 2021).
The same set of dominant OTUs was used to investigate the main physicochemical drivers behind the microbial compositions of the evolving Podzol profiles using canonical correspondence analysis (CCA).The analysis was carried out using the vegan package (https://CRAN.R-project.org/package=vegan,version 2.6-4, accessed on 11 October 2022) in RStudio (http://www.rstudio.com/,the version 2021.09.0, accessed on 29 September 2021).The correlation between the factors was evaluated using the VIFs function.An analysis of variance (ANOVA), including 999 permutations, was used to estimate the significance of certain factors (comprising the model) in shaping the microbial diversity.Visualization was performed using Past3 software (version 1.0.0.0).
All sequences were deposited in the Sequence Read Archive of the National Center for Biotechnology Information in the dataset SUB4606526.

General Soil Characteristics
The initial substrate was slightly acidic (Table S1).As the soil profile developed (from the initial M1 to the mature M5 stages), the acidity of the soil (generally, topsoil-O and E horizons) increased (Tables S1 and S2).This tendency is directly associated with organic matter accumulation and redistribution, and particularly in the observed peak in organic-matter content in the 30-year-old site, leading to a substantial decrease in the corresponding pH values.The accumulation of TOC (as well as soil N) in the O horizons of podzolic soils (stages M3-5) is supported by ANOVA with Fisher's LSD test (Table S2, Figure S2,).The increase in accumulation of soil organic matter in the mineral solum was fixed in the 70-year-old soil (Tables S1 and S2, Figures 3 and S2), indicating the intensification of soil organo-mineral interaction.The development of the corresponding physical and chemical transformations led to a more pronounced differentiation of the vertical soil profile and its division into two spodic sub-horizons-BHF and BF-which had disappeared from the soil profile by the climax stage.Here and afterwards, the climax stage of soil development or corresponding microbial succession will be defined as a stable biome (in terms of energy, diversity estimates, etc.) demonstrating typical characteristics of the certain climatic zone [69].The soils were classified as sandy textured due to the very low clay content.The basal respiration as well as metabolic coefficient demonstrated the tendency of decrease in the mature soil (stage M5) (Table S2).This indicates the stabilization of the soil formation process.

Real-Time PCR
An increase in the number of bacteria and archaea was detected in the mature soil (Figure 4).The bacteria were associated with the upper horizons, whereas the number of archaea increased in the BC mineral horizon.This distribution of prokaryotic communities depended on the organic matter accumulation.The ANOVA analysis results can be seen in Tables S3 and S4.All the differences observed are statistically significant.

Alpha Diversity of Soil Microbiomes
All the measured diversity indices tended to decrease from the initial to the mat stage of soil succession.In Stages M1 and M2, the alpha diversity was not significan different in all soil horizons.In M3, a clear separation of the more diverse upper hori (O) was observed.In M4 and M5, the vertical pattern of microbial diversity was flatte or even biased toward the lower soil horizons (Table 1, Figure S3).

Alpha Diversity of Soil Microbiomes
All the measured diversity indices tended to decrease from the initial to the mature stage of soil succession.In Stages M1 and M2, the alpha diversity was not significantly different in all soil horizons.In M3, a clear separation of the more diverse upper horizon (O) was observed.In M4 and M5, the vertical pattern of microbial diversity was flattened or even biased toward the lower soil horizons (Table 1, Figure S3).C-initial parent material, not altered by pedogeneis, Oi-initial litter, not altered by organic matter transformation process, Oe-litters slightly transformed by enzymatic processes, AY-humus accumulative superficial soil horizon, E-podzolic soil horizon, BF-illuvial iron-oxide horizon, BHF-illuvial iron-oxide-humic complexes horizon, BC-transitional horizon between illuvial horizon and parent material.* The E (eluvial podzol) horizon of M3 stage significantly differed in thickness and physical and chemical properties from the same horizons of more mature soils (M4 and M5 sites).However, its morphological difference from the underlying illuvial (BF) horizon allows it to be identified as a separate horizon and called "enthic podzol horizon" (e).

Beta Diversity of Soil Prokaryotic Communities
The beta diversity demonstrates the correlation of the microbial community structure with soil genesis (Figure 5).
In the initial stages of succession, the prokaryotic communities of the soil layers (the nascent horizons) had a similar structure (Stage M1).With the accumulation of organic matter, the process of vertical diversification and the formation of the predecessors of the soil horizons of the mature soil began (Stage M2).Although the diversification of the prokaryotic communities of the soil horizons was still not significant, according to the PERMANOVA test (p = 0.097, indicating a difference at the trend level), particularly in the second stage, AY and BF horizons formed and took the intermediate position between the mineral and organic soil horizons, which completely separated in the M3 stage of soil formation.In the M3 stage, a clear diversification of the organogenic part of the profile-O and E horizons-was observed, which showed the highest levels of dissimilarity compared to the mineral horizons (p < 0.05), which tended to have quite similar microbial composition.
well as a sharp differentiation in the O horizon, resulting in the complete separation organic and eluvial (E) horizons in the preclimax and climax soil profiles.Except fo profile points (e.g., replicates of the experiment in the case of the humus-accum horizons BF and BHF) that occupied the intermediate position between the mine organic soil horizons (Figure 5), the illuvial humus-accumulative BF horizons we well defined in the mature stage of soil genesis (Figure 6).Stage M4 was characterized by the further diversification of the mineral horizons as well as a sharp differentiation in the O horizon, resulting in the complete separation of the organic and eluvial (E) horizons in the preclimax and climax soil profiles.Except for some profile points (e.g., replicates of the experiment in the case of the humus-accumulative horizons BF and BHF) that occupied the intermediate position between the mineral and organic soil horizons (Figure 5), the illuvial humus-accumulative BF horizons were also well defined in the mature stage of soil genesis (Figure 6).

Taxonomic Composition of the Core Prokaryotic Communities in the Different Stages of Soil Succession
To assess changes in the soil prokaryotic communities' structure during evolution of the soil profile, the core phylotypes were identified for each succession stage (M1-M5) (Figure 7).The composition of the cores was determined from pooled replicates of the samples of each horizon of the site.This procedure was applied to identify all possible trends in the compositions of each site.The core prokaryotic communities contained the following number of identified phylotypes: M1-41; M2-71; M3-59; M4-97; and M5-

Taxonomic Composition of the Core Prokaryotic Communities in the Different Stages of Soil Succession
To assess changes in the soil prokaryotic communities' structure during evolution of the soil profile, the core phylotypes were identified for each succession stage (M1-M5) (Figure 7).The composition of the cores was determined from pooled replicates of the samples of each horizon of the site.This procedure was applied to identify all possible trends in the compositions of each site.The core prokaryotic communities contained the following number of identified phylotypes: M1-41; M2-71; M3-59; M4-97; and M5-94.No phylotypes were present in all five stages simultaneously.Most of the identified phylotypes demonstrated continuous transitions through two or even three successive stages of soil development, such as from M1 to M2 (11 OTUs), M1 to M2 and M3 (11 OTUs), M3 to M4 Soil Syst.2024, 8, 9 13 of 25 (3 OTUs), M3 to M4 and M5 (3 OTUs), and M4 to M5 (the largest group with 55 OTUs).However, a small number of phylotypes had random patterns of distribution, such as one shared OTU between M2 and M4 (ACB, Sg1, Acidobacteria), one between M2 and M5 (Cyanobacteria), and two between M1 and M3.Thus, the overall core diversity could be divided into two large groups-those phylotypes that dominated the initial stages of soil formation (from M1 to M3), and those that dominated in more or less developed soil profiles in Stages M4 and M5.Stages M1 and M2 were dominated by archaea from the family Nitrososphaeraceae (1 OTU) and bacteria from the families Microscillaceae (1 OTU), Xanthomonadaceae (1 OTU), and Nitrosomonadaceae (1 OTU).This transition was also characterized by the presence of various acidobacterial subdivisions, particularly the Pyrinomonadales and Holophaga groups.The transition to Stage M3 contained archaea from the family Nitrososphaeraceae (1 OTU) and bacteria belonging to the Microscillaceae (Bacteroidetes, 1 OTU), Chitinophagaceae (Bacteroidetes, 3 OTUs), Micrococcales (Bacteroidetes, 2 OTUs), Burkholderiaceae (Gammaproteobacteria, 3 OTUs), Gaiellales (Actinobacteria, 1 OTU), Rhizobiales (Alphaproteobacteria, 1 OTU), Nitrosomonadaceae (Gammaproteobacteria, 3 OTUs), and Pyrinomonadales, particularly genus RB41 (Acidobacteria, 3 OTUs) taxonomic subdivisions.Very often, the bacteria that dominated Stages M2 and M3 had been eliminated from the microbial community of Stages M4 and M5.In particular, only Stages M2 and M3 contained bacteria from the phyla Nitrospira and unidentified representatives of the Chloroflexi.The least number of transitions were observed between Stages M3 and M4-M5, including unidentified representatives of the Gemmatimonadetes and bacteria from the family Corynebacteriales, which were absent from M1 and M2.Other bacteria, such as the Rhizobiaceae family and acidobacteria from the genus RB41, remained in the core prokaryotic communities of the mature soils, but it changed in their species and genus compositions.
94.No phylotypes were present in all five stages simultaneously.Most of the identified phylotypes demonstrated continuous transitions through two or even three successive stages of soil development, such as from M1 to M2 (11 OTUs), M1 to M2 and M3 (11 OTUs), M3 to M4 (3 OTUs), M3 to M4 and M5 (3 OTUs), and M4 to M5 (the largest group with 55 OTUs).However, a small number of phylotypes had random patterns of distribution, such as one shared OTU between M2 and M4 (ACB, Sg1, Acidobacteria), one between M2 and M5 (Cyanobacteria), and two between M1 and M3.Thus, the overall core diversity could be divided into two large groups-those phylotypes that dominated the initial stages of soil formation (from M1 to M3), and those that dominated in more or less developed soil profiles in Stages M4 and M5.Stages M1 and M2 were dominated by archaea from the family Nitrososphaeraceae (1 OTU) and bacteria from the families Microscillaceae (1 OTU), Xanthomonadaceae (1 OTU), and Nitrosomonadaceae (1 OTU).This transition was also characterized by the presence of various acidobacterial subdivisions, particularly the Pyrinomonadales and Holophaga groups.The transition to Stage M3 contained archaea from the family Nitrososphaeraceae (1 OTU) and bacteria belonging to the Microscillaceae (Bacteroidetes, 1 OTU), Chitinophagaceae (Bacteroidetes, 3 OTUs), Micrococcales (Bacteroidetes, 2 OTUs), Burkholderiaceae (Gammaproteobacteria, 3 OTUs), Gaiellales (Actinobacteria, 1 OTU), Rhizobiales (Alphaproteobacteria, 1 OTU), Nitrosomonadaceae (Gammaproteobacteria, 3 OTUs), and Pyrinomonadales, particularly genus RB41 (Acidobacteria, 3 OTUs) taxonomic subdivisions.Very often, the bacteria that dominated Stages M2 and M3 had been eliminated from the microbial community of Stages M4 and M5.In particular, only Stages M2 and M3 contained bacteria from the phyla Nitrospira and unidentified representatives of the Chloroflexi.The least number of transitions were observed between Stages M3 and M4-M5, including unidentified representatives of the Gemmatimonadetes and bacteria from the family Corynebacteriales, which were absent from M1 and M2.Other bacteria, such as the Rhizobiaceae family and acidobacteria from the genus RB41, remained in the core prokaryotic communities of the mature soils, but it changed in their species and genus compositions.Some bacterial families, or even bacterial phyla, appeared only in the final stages of soil profile development, among them being entire phyla: WPS-2, Armatimonadetes, Planctomycetes, and almost all (except one OTU in Stage M2) representatives of the Verrucomicrobia; the orders and classes Corynebacteriales, Solirubrobacterales (Actinobacteria), AD3 (Chloroflexi), Bacilli (Firmicutes), Acidobacteriales, Acidobacteria from Subgroups 1 and 2; and the genus Acidothermus (Actinobacteria).Along with the appearance of the new groups of bacteria in the later successional stages, the taxonomic compositions of the previous prokaryotic communities changed dramatically, transforming the components of the Rhizobiales and Chitinophagaceae subgroups (Figure 7).

Dynamic Changes in the Compositions of the Prokaryotic Communities Inhabiting Different Soil Horizons through Time
To study the dynamic changes in the prokaryotic communities during soil development, a set of dominant taxa was selected for each stage (M1-M5), including 47 phylotypes, the abundances of which exceeded 1% in at least one soil sample (Table S5).
Several groups of taxa were found in specific soil horizons.The first group was found exclusively in the BC-C mineral horizons of the mature soil, the second accumulated in the upper and (regularly) in the lower part of the soil profile (in zones of active root growth), and the third occurred in the BF horizon of the mature soil (Figure 8A-C).The first group (Figure 8A) included Myxococcales, MND1, Xanthobacteriaceae (Proteobacteria), and KD 4-26 (Chloroflexi).The second group (Figure 8B) included representatives of Xanthobacteriaceae, Bradyrhizobium, WD260, Burkholderia multivorans (Proteobacteria), Solirubrobacteraceae, Subtercola sp.(Actinobacteria), Chitinophagaceae (Bacteroidetes), and Acidobacteriales (Acidobacteria).The third group (Figure 8C) comprised Subgroup 2, RB41 (Acidobacteria), Massilia sp., Beijerinckiaceae, IS44 (Proteobacteria), Candidatus Udaeobacter (Verrucomicrobia), TK10, AD3 (Chloroflexi), Gemmatimonadaceae (Gemmatimonadetes), and unidentified archaea from the phylum Thaumarchaeota.The analysis included two stages that highlighted the development of the primary Podzol that are known to be the main soil-forming processes under coniferous forests.The first was development of an Embryopodzol (Stage M3), with only a nascent podzolic (e) horizon, and the second was the preclimax stage (M4), where the podzolic process manifested as a clearly pronounced illuvial-eluvial differentiation of the soil profile, along with a thick organogenic horizon, including forest litter (O) (Figure 9).The analysis included two stages that highlighted the development of the primary Podzol that are known to be the main soil-forming processes under coniferous forests.The first was development of an Embryopodzol (Stage M3), with only a nascent podzolic (e) horizon, and the second was the preclimax stage (M4), where the podzolic process manifested as a clearly pronounced illuvial-eluvial differentiation of the soil profile, along with a thick organogenic horizon, including forest litter (O) (Figure 9).The model included the factors total C org , depth, clay content, and pH, and these totally explained 41.9% of the entire variation in the values of the current dataset.The CCA plot allowed us to distinguish several zones.The first zone was characterized by positive values of both the CCA1 and CCA2 axes, explaining 19% of the total variance, which was correlated with the accumulation of C org in the soil and corresponded with diversification in the topsoil of Stage M4.This is the most probable reason for the increase in proteobacteria from the Gammaproteobacteria, Xanthobacteraceae (Rhizobiales), and Bradyrhizobium (Rhizobiales) groups as well as bacteria from the family Chitinophagaceae (Bacteroidetes) and genera Chthoniobacter (Verrucomicrobia) and Mycobacterium (Actinobacteria) in Stage M4.
The negative values of both the CCA1 and CCA2 axes corresponded to the mineral (BF3, BC3, and C3) and illuvial (BF4) soil horizons.The microbial shifts here were likely associated with the more neutral pH values of the soil in solution as well as with the illuviation process and the corresponding increase in the proportion of clay particles (clay factor).Phylotypes of Acidobacteria of Subgroups 6 and 4 (Blastocatellia, Pyrinomonadaceae), Proteobacteria of the family Nitrosomonadaceae (Gammaproteobacteria), Verrucomicrobia phylotypes (Candidatus Udaeobacter), and the Archaea (Thaumarchaeota) were associated with these samples.Negative values of the CCA2 axis determined the overall abundance of taxa associated with the mineral horizons.Of these, Acidobacteria of Subgroup 2, the order Acidobacteriales, and Actinobacteria of the family Solirubrobacteraceae were associated with the podzolic horizon of the preclimax stage (E4, positive values of the CCA1 and negative values of the CCA2 axes).In the mineral horizons of the preclimax stage (BF4 and, to a greater extent, C4), there was an increase in the abundances of Acidobacteria of Subgroup 2, phylotypes of Gemmatimonadetes and Chloroflexi (AD3), and Archaea (Thaumarchaeota) (Figure 9).

Discussion
The investigated soils represented classic Podzol soil evolution through ecogenetic succession.In general, the soil profile development rate was comparable with those studied previously [3][4][5]7].The main reason for this is the origin of the parent materials of the mentioned studies, which are extremely acidic.Due to the specificity of the soil-forming process-the sandy texture of the parent material [70], the acidic reaction of coniferous fallout [71], and the flushing water regime with high precipitation-the soil development in the chronosequence was rapid, and the soil profile became well developed in 70 years, expressing all the genetic horizons that are typical of Podzols of the southern taiga [72].
As the plant cover developed, the biomass of the associated soil microbial community also increased due to fresh organic matter deposition in the organogenic horizons (topsoils) in the form of plant litter and sloughing of root cells during growth, which in turn may explain the increase in bacterial biomass (estimated as amounts of ribosomal operons) in the upper horizons of the mature zonal soils of the preclimax and climax (M4 and M5) stages (Figure 4).When analyzing soil prokaryotic community in the context of soil ecogenesis, the main trend of soil microbiological succession is observed at the level of the quantitative indicators of the microbial diversity, primarily related to the redistribution of organic matter and, in particular, to the shifts in its composition.So, we might suppose the occurrence of bacteria that utilize relatively simple components of the organic matter at the early stages, whereas later stages are characterized by a more diverse community.
Moreover, under intensive leaching of the soil biomass, rates of mineralization of plant residues and humification are suppressed, and the majority of the products of organicmaterial transformation are washed out of the soil profile, accumulating in the illuvial horizon [73].This organic matter concentration can contribute to the maintenance of an oligotrophic microbiota in the lower mineral soil layers, explaining the relative increase in archaeal ribosomal operons in the mineral horizons of the mature soils.Furthermore, the lower portion of the BF horizon is usually poor in organic matter and nutrients, contributing to a relative increase in the abundance of the Archaea, which are more resistant to various environmental stresses [74].For example, an increase in the archaeal-to-bacterial ratio has been reported from the late stages of pedogenesis (more than 120,000 years) in the lower mineral horizons, indicating the adaptation of the Archaea to the nutrient and energy-poor environment typical of the subsoil [42].
The alpha diversity data clearly demonstrated a strong dependence of soil differentiation processes on the organic matter supply and podzolization with the overall level of prokaryotic diversity, sharply decreasing in the latest stages.In general, this trend is common in succession processes with the early stages of biome development characterized by dynamic changes in the architecture and organization of ecological niches that are occupied at the beginning by various copiotrophic microorganisms, which are then replaced by well-specialized oligotrophs, marking the climax stage of soil development [69].During successive differentiation of the vertical soil profile, the diversity indices gradually decreased, marking two hot-spots of relatively high microbial diversity in the organogenic (O) and humus-accumulative (BF), or even mineral BC, horizons (e.g., the M4 stage).Starting from the M3 stage, the tendency toward separation of the podzolic horizon (E) became obvious due to intensification of the leaching processes (i.e., podzolization).Among the named tendencies, the accumulation of microbial diversity in the lower soil layers merits special discussion.Well-established soil mineral horizons are enriched with organic matter from the upper soil layers.In addition, they contain large amounts of clay.A high clay content, in turn, can cause the accumulation of extracellular DNA, which is known to form stable complexes with clay minerals [75].Both factors can stimulate the growth of microbial diversity in these horizons.
The beta-diversity analysis revealed regularities between vertical soil differentiation and the changes that occurred in the microbial community structure.The initial stage demonstrated an absence of differences in the microbial diversity values between samples.Beginning from Stage M2, the process of differentiation gained strength, resulting in the formation of a full set of the predecessors of the mineral and organogenic horizons in the M3 stage.Thus, slightly more than 30 years was enough time for the soil to completely restore its zonal characteristics.The subsequent processes simply reinforced this differentiation without adding anything fundamentally new to the soil structure.Detailed analysis of the transition from the M3 stage of embryopodzol to the mature podzol profile at the M4 stage demonstrated the continuation of the process of acidification of the upper part of the soil profile (Figure 9), which is associated with the accumulation of organic matter in the O horizon of the M4 stage.Also, in relation to mineral horizons, their preferential enrichment with clay particles was observed compared to the upper part of the profile, which was also more characteristic of the more mature preclimax soil (M4).
The continuous development of the soil profile was fully reflected in the taxonomic changes in the microbial community, ending with the formation of the prokaryotic communities, which are typical of podzolic soils in the temperate zone in Russia [76,77].The distribution of specific taxa in the soil profiles showed the dominance of Proteobacteria in the zero stage of soil formation, Chloroflexi in the intermediate stage of soil development, affected by illuviation processes, and Acidobacteria in the preclimax and climax stages.These shifts were quite expected and can be explained by ecological changes.At the beginning of ecosystem development, the dominant positions were occupied by fast-growing copiotrophs (concentrated primarily within the Proteobacteria), which were succeeded by a well-specialized oligotrophic microbiota (e.g., Acidobacteria, Verrucomicrobia, Gemmatimonadetes, AD3, WSP-2, Planctomycetes, and Armatimonadetes) [78][79][80][81].
In addition to these ecological explanations, prokaryotic communities' changes can occur via, for example, intensive acidification of the soil profile, which promotes the growth of Acidobacteria in the latest stages of soil formation; these are known to dominate in soil with pH values of 3 to 5 [82].A change in the plant community, accompanied by a greater amount of fresh organic matter and increased thickness of the organogenic horizon (forest litter), may explain an increase in the abundance of Verrucomicrobia in the microbial community of mature soils.Bacteria belonging to this phylum are known to gravitate to the forest floor [83].
To understand the reasons for, and specific traits of, the adaptive patterns of prokaryotic communities in the course of soil genesis, it is important to consider the changes occurring in the taxonomic structure at lower taxonomic levels.In this respect, the circular representation of changes in the basic (core) community, presented in Figure 7, as well as the gradual change in abundance of the taxa in the soil profiles (Figure 8) is relevant.The adaptation strategies of the microorganisms evolve during the evolution of the soil ecosystem alongside diversification and changes in the dominant components of the microbial community, which are, in turn, responsible for certain soil biological processes.The intensity of these processes also changes during the development of genetic soil horizons.
The first of these processes is root growth and the corresponding migration of the rhizosphere microorganisms to the lower soil layers.For example, Figure 8B presents various bacterial groups that were evenly distributed in the soil profile in the initial stages of succession, while in the latter stages, these bacteria were concentrated in the upper soil layers and mineral BC horizon.The presence of bacteria belonging to the Bradyrhizobium genus or Solirubrobacteraceae and Xanthobacteraceae families points to a rhizospheric origin for this part of the soil prokaryotic community [84,85].
Unfortunately, this study is lacking comprehensive information on the distance between the trees and the points of soil sampling as well as the complete list of the plant species inhabiting this biome.Obviously, it might be the main predictor of the microbial distribution among the soil profile.
The second of these processes involves the N cycle.Analysis of the investigated prokaryotic communities revealed a succession of prokaryotes associated with the transformation of N compounds during soil genesis.In the initial stages, this group was represented by the Archaea (family Nitrososphaeraceae) and a few Nitrosomonadales, the abundance of which increased in subsequent stages, particularly that of the Embryopodzol (M3).These stages were also enriched with cyanobacteria (particularly M2), which were presumably involved in the primary structuring of the soil and its enrichment with N and organic matter [86].The presence of cyanobacteria in the first stages of the regeneration of dump soils has been noted previously [87].Probably, the detection of cyanobacteria in the mineral horizons of the preclimax soils occurred due to the entry of DNA from these microorganisms from the upper horizons during pedogenesis.In the latest stages of succession (M4, M5), Nitrospirae lineages and N-fixators from the Rhizobiales order were included in the N transformation [88].
The final process is bacterial degradation of the main soil biopolymers.Along with the accumulation of fresh organic matter, the proportion and diversity of cellulose, hemicellulose, pectin, and lignin degraders [89,90] belonging to the Bacteroidetes (families Chitinophagaceae, Sphingobacteriaceae), as well as gamma-proteobacteria (family Burkholderiaceae), gradually increased.This group can also include Actinobacteria in the final stages of soil development.
Because the earlier stages (especially M2 and M3) were dominated by the coccoid copiotrophic representatives of this phylum, such as Micrococcales, the latter stages included bacteria with pseudo-(Corynebacteriales) and true (Actinomycetales) mycelium, generally referred to the autochthonous soil microbiota, which utilize humus substances and soil biopolymers, and these are known to integrate the soil community chemically by synthesizing various antibiotics and signal molecules and mechanically aggregating microbial cells and soil particles [90].
In this study, their co-occurrence can be seen in Figure 8A-C, especially Figure 8B, where the bacteria forming the rhizosphere community are represented.The plant-associated nature of this part of the microbial community is marked by indicator phylotypes, such as Bradyrhizobium and bacteria belonging to the families Xanthomonadaceae and Solirubrobacteraceae [84,91,92].So, we might expect the same ecological features or some beneficial function for plant nutrition/development, etc., from the other members of this group, such as from the WD260 phylotype (Proteobacteria).The existing information on the distribution and functional opportunities of these bacteria is quite poor.However, they are likely to have the widest ecological capabilities in both aquatic and terrestrial environments, including disturbed soils [93,94].

Conclusions
The studied model of Podzol soil genesis showed the continuous evolution of the prokaryotic communities' structure, which entirely corresponded with the process of soil organic matter redistribution and formation of the vertical soil profile (podsolization, illuviation of iron oxides).
The results allowed us to trace the relationship between the structure of the soil prokaryotic community and the peculiarities of the evolution of the podzolic soil profile as well as to identify the microbial indicators and drivers of primary pedogenesis.Thus, the direct task was solved-the influence of soil age and the degree of differentiation of its profile on the microbiome was studied.In the future, it is necessary to solve a more complex-inverse problem-how to control soil formation process with the help of functional composition of the microbiome.The solution of this problem will be an important step for soil cover restoration and land reclamation.
Beta-diversity data analysis, as well as analysis of the core prokaryotic communities in each stage, revealed a comparatively large distance and sharp transition, in terms of microbial diversity, between the Embryopodzol (Stage M3) and the later stages of soil development.Thus, in further research, it will be reasonable to complete the chronosequence with at least one intermediate stage of soil formation, ranging from 40 to 45 years old.
The methods applied to studying this multiplex system are generalized and sometimes difficult.In particular, the representation of changes between subsequent stages of soil genesis, in terms of the core prokaryotic communities combined from the horizons and replicates, was adequate but too generalized.Unfortunately, we currently lack an elaborate mathematical model that could highlight the hidden abundance patterns of certain taxa.It would have also been useful to have included different pools of nucleic acids in this study, such as DNA, RNA, and extracellular DNA, and to have been able to divide the soil into its natural microenvironments--that is, soil aggregates.Currently, studying the microbial succession in a developing soil can contribute to various biological and ecological questions, including the discovery of new symbiotic relationships and the matching of specific functional profiles to yet-uncultured prokaryotes.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/soilsystems8010009/s1, Figure S1: Rarefaction curves for alphadiversity analysis.title; Figure S2: Box-plots for the physic-chemical parameters; Figure S3: Box-plots for the alpha-diversity measurements.Table S1: Summary of ANOVA results of the abundance of bacteria and archaea estimated by quantitative PCR; Table S2: Tukey multiple comparisons of means (95% family-wise confidence level) based on ANOVA results of the abundance of bacteria and archaea estimated by quantitative PCR; Table S3: The list of dominant taxa (abundances exceeds 1% of the soil sample community), which are present in at least one of the soil horizons and the corresponding soil samples (the core community).

Figure 1 .
Figure 1.Geographical location of the sample collection sites.Northwest of Russia, Leningrad region, Kirovskiy district.M1-M5-indicating the points of sampling.

Figure 1 .
Figure 1.Geographical location of the sample collection sites.Northwest of Russia, Leningrad region, Kirovskiy district.M1-M5-indicating the points of sampling.

28 Figure 3 .
Figure 3.The values of the soil's main physicochemical characteristics and microbiological parameters of the soil horizons in M1-M5 samples representing their change the course of the soil ecogenesis.Abbreviations: (TOC-total organic carbon, N-total nitrogen, pHH20-water PH measurement, basal respiration, C mic biomass-microbial carbon biomass, metabolic coefficient) of the soil horizons in M1-M5 samples.C-initial parent material, not altered by pedogeneis, Oi-initial litter, not altered by organic matter transformation process, Oe-liters, slightly transformed by enzymatic processes, AY-humus accumulative superficial soil horizon, BF-illuvial iron-oxide horizon, (e)entic podzol horizon, BHF-illuvial iron-oxide-humic complexes horizon.

Figure 3 .
Figure 3.The values of the soil's main physicochemical characteristics and microbiological parameters of the soil horizons in M1-M5 samples representing their change the course of the soil ecogenesis.Abbreviations: (TOC-total organic carbon, N-total nitrogen, pHH 2 O-water PH measurement, basal respiration, C mic biomass-microbial carbon biomass, metabolic coefficient) of the soil horizons in M1-M5 samples.C-initial parent material, not altered by pedogeneis, Oi-initial litter, not altered

Figure 4 .
Figure 4. Number of bacteria and archaea expressed in amounts of 16S rRNA or internal transcri tracer sequences per gram of soil (y-axis), estimated by quantitative PCR.Error bars represent SDs of the values.M1-M5 represent stages in the soil succession.C-initial parent material,altered by pedogeneis, Oi-initial litter, not altered by organic matter transformation process, O litter, slightly transformed by enzymatic processes, AY-humus accumulative superficial soil h zon, BF-illuvial iron-oxide horizon, (e)-entic podzol horizon, BHF-illuvial iron-oxide-hu complexes horizon.

Figure 4 .
Figure 4. Number of bacteria and archaea expressed in amounts of 16S rRNA or internal transcribed tracer sequences per gram of soil (y-axis), estimated by quantitative PCR.Error bars represent the SDs of the values.M1-M5 represent stages in the soil succession.C-initial parent material, not altered by pedogeneis, Oi-initial litter, not altered by organic matter transformation process, Oe-litter, slightly transformed by enzymatic processes, AY-humus accumulative superficial soil horizon, BF-illuvial iron-oxide horizon, (e)-entic podzol horizon, BHF-illuvial iron-oxide-humic complexes horizon.BC-transitional horizon between illuvial horizon and parent material.

Figure 5 .
Figure 5. Beta diversity in the soil horizons of different genesis.Distances between the samp calculated using the weighted UniFrac method.PCoA1, 2, and 3 represent the axes of the PC with the percentage of explained variance.O, E, AY, BHF, BF, BC, and C are the soil horizo M5 represent sequential stages in the soil succession.Gray dots are marking all the samples the range of microbial dynamic shifts in each stage.

Figure 5 .
Figure 5. Beta diversity in the soil horizons of different genesis.Distances between the samples were calculated using the weighted UniFrac method.PCoA1, 2, and 3 represent the axes of the PCoA plot with the percentage of explained variance.O, E, AY, BHF, BF, BC, and C are the soil horizons.M1-M5 represent sequential stages in the soil succession.Gray dots are marking all the samples to show the range of microbial dynamic shifts in each stage.

28 Figure 6 .
Figure 6.Bray-Curtis dissimilarities of the soil samples.O, E, AY, BHF, BF, BC, and C are soil horizons.BF' refers to the lower part of BF horizon in M3 (these parts differed in density, however this difference was within density values of the same horizon).M1-M5 represent stages in the soil succession.

Figure 6 .
Figure 6.Bray-Curtis dissimilarities of the soil samples.O, E, AY, BHF, BF, BC, and C are soil horizons.BF' refers to the lower part of BF horizon in M3 (these parts differed in density, however this difference was within density values of the same horizon).M1-M5 represent stages in the soil succession.

Figure 7 .
Figure 7. Core community of soil microbiomes.Diamonds represent the presence of taxa in the core community at certain stages in the soil succession: Stage M1-yellow; M2-light green; M3-dark green; M4-blue; and M5-violet.

Figure 8 .
Figure 8. Distribution (heat map) of specific taxa in the soil profiles (the pictures were drawn to show the exemplary spatial architecture of the soil profile, using the depth measurement of the soil horizons).Only taxa with significant differences in their abundances (expressed as the percent of the total microbial number of the amplicon library, marked with color range from red-the most abundant to white-the less abundant) are shown.A-C are the phylotype groups, separated according to their tendency to concentrate in particular parts of the soil profile.O, E, AY, BHF, BF, BC, and C are soil horizons.

Figure 8 . 9 3. 8 .
Figure 8. Distribution (heat map) of specific taxa in the soil profiles (the pictures were drawn to show the exemplary spatial architecture of the soil profile, using the depth measurement of the soil horizons).Only taxa with significant differences in their abundances (expressed as the percent of the total microbial number of the amplicon library, marked with color range from red-the most abundant to white-the less abundant) are shown.A-C are the phylotype groups, separated according to their tendency to concentrate in particular parts of the soil profile.O, E, AY, BHF, BF, BC, and C are soil horizons.

Table 1 .
Alpha-diversity indices.The letters denote the significance difference according to ANOVA (p-value < 0.05).Different letters (a-e) mark significant differences in the values of alpha-diversity indices of the correspondence samples within each column.