# Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile

Andrey Aksenov, Liudmila Shirokova, Oksana Ya. Kisil, Sofia Kolesova, Artem Lim, Darya Kuzmina, Sophie Pouillé, Marie Alexis, Maryse Castrec-Rouelle, Sergey Loiko, et al.

► To cite this version:

Andrey Aksenov, Liudmila Shirokova, Oksana Ya. Kisil, Sofia Kolesova, Artem Lim, et al. Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile. Diversity, 2021, 13 (7), pp.328. 10.3390/d13070328 . hal-03974645

HAL Id: hal-03974645

https://cnrs.hal.science/hal-03974645

Submitted on 18 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile

Andrey Aksenov, Liudmila Shirokova, Oksana Ya. Kisil, Sofia Kolesova, Artem Lim, Darya Kuzmina, Sophie Pouillé, Marie Alexis, Maryse Castrec-Rouelle, Sergey Loiko, et al.

► To cite this version:

Andrey Aksenov, Liudmila Shirokova, Oksana Ya. Kisil, Sofia Kolesova, Artem Lim, et al. Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile. Diversity, 2021, 13 (7), pp.328. 10.3390/d13070328 hal-03974645

HAL Id: hal-03974645

https://cnrs.hal.science/hal-03974645v1

Submitted on 18 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile

Andrey S. Aksenov <fnref n="1" />,* ID, Liudmila S. Shirokova <fnref n="2" />,<fnref n="3" />, Oksana Ya. Kisil <fnref n="1" />, Sofiia N. Kolesova <fnref n="1" />, Artem G. Lim <fnref n="4" /> ID,
Darya Kuzmina <fnref n="4" />, Sophie Pouillé <fnref n="5" />, Marie A. Alexis <fnref n="5" />, Maryse Castrec-Rouelle <fnref n="5" /> ID, Sergey V. Loiko <fnref n="4" /> ID
and Oleg S. Pokrovsky <fnref n="3" /> ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Aksenov, A.S.; Shirokova, L.S.; Kisil, O.Y.; Kolesova, S.N.; Lim, A.G.; Kuzmina, D.; Pouillé, S.; Alexis, M.A.; Castrec-Rouelle, M.; Loiko, S.V.; et al. Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile. *Diversity* **2021**, *13*, 328. https://doi.org/10.3390/d13070328

Academic Editors: Michael Wink, Boris A. Levin, Yulia V. Bespalaya and Maria Teresa Ceccherini Guicciardini

Received: 4 June 2021
Accepted: 14 July 2021
Published: 17 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Permafrost peatlands, containing a sizable amount of soil organic carbon (OC), play a pivotal role in soil (peat) OC transformation into soluble and volatile forms and greatly contribute to overall natural CO₂ and CH₄ emissions to the atmosphere under ongoing permafrost thaw and soil OC degradation. Peat microorganisms are largely responsible for the processing of this OC, yet coupled studies of chemical and bacterial parameters in permafrost peatlands are rather limited and geographically biased. Towards testing the possible impact of peat and peat pore water chemical composition on microbial population and diversity, here we present results of a preliminary study of the western Siberia permafrost peatland discontinuous permafrost zone. The quantitative evaluation of microorganisms and determination of microbial diversity along a 100 cm thick peat soil column, which included thawed and frozen peat and bottom mineral horizon, was performed by RT-PCR and 16S rRNA gene-based metagenomic analysis, respectively. Bacteria (mainly *Proteobacteria*, *Acidobacteria*, *Actinobacteria*) strongly dominated the microbial diversity (99% sequences), with a negligible proportion of archaea (0.3–0.5%). There was a systematic evolution of main taxa according to depth, with a maximum of 65% (*Acidobacteria*) encountered in the active layer, or permafrost boundary (50–60 cm). We also measured C, N, nutrients and ~50 major and trace elements in peat (19 samples) as well as its pore water and dispersed ice (10 samples), sampled over the same core, and we analyzed organic matter quality in six organic and one mineral horizon of this core. Using multiparametric statistics (PCA), we tested the links between the total microbial number and 16S rRNA diversity and chemical composition of both the solid and fluid phase harboring the microorganisms. Under climate warming and permafrost thaw, one can expect a downward movement of the layer of maximal genetic diversity following the active layer thickening. Given a one to two orders of magnitude higher microbial number in the upper (thawed) layers compared to bottom (frozen) layers, an additional 50 cm of peat thawing in western Siberia may sizably increase the total microbial population and biodiversity of active cells.

**Keywords:** permafrost; peat; pore water; bacteria; number; metabolism; active layer

<footnote n="1">¹ Arctic Biomonitoring Laboratory, Northern (Arctic) Federal University, 17 Naberezhnaya Sev. Dviny, 163002 Arkhangelsk, Russia; kisil.oksana.jroslavovna@gmail.com (O.Y.K.); sofianvanwyn@gmail.com (S.N.K.)</footnote>

<footnote n="2">² N. Laverov Federal Center for Integrated Arctic Research, Russian Academy of Sciences, 23 Naberezhnaya Sev. Dviny, 163000 Arkhangelsk, Russia; liudmila.shirokova@get.omp.eu</footnote>

<footnote n="3">³ Geosciences and Environment Toulouse, UMR 5563 CNRS, 14 Avenue Edouard Belin, 31400 Toulouse, France; oleg.pokrovsky@get.omp.eu</footnote>

<footnote n="4">⁴ BIO-GEO-CLIM Laboratory, Tomsk State University, 36 Lenina, Build. 1, 634050 Tomsk, Russia; lim_artyom@mail.ru (A.G.L.); kuzmina.d.m.95@gmail.com (D.K.); s.loyko@yandex.ru (S.V.L.)</footnote>

<footnote n="5">5 UMR 7619 Métis, Sorbonne Université, F-75252 Paris, France; 24sophiepouille@gmail.com (S.P.); marie.alexis@sorbonne-universite.fr (M.A.A.); maryse.rouelle@sorbonne-universite.fr (M.C.-R.)</footnote>

<footnote>* Correspondence: a.s.aksenov@narfu.ru</footnote>
Article

# Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile

Andrey S. Aksenov <fnref n="1" />,* , Liudmila S. Shirokova <fnref n="2" />,<fnref n="3" />, Oksana Ya. Kisil <fnref n="1" />, Sofiia N. Kolesova <fnref n="1" />, Artem G. Lim <fnref n="4" />, Darya Kuzmina <fnref n="4" />, Sophie Pouillé <fnref n="5" />, Marie A. Alexis <fnref n="5" />, Maryse Castrec-Rouelle <fnref n="5" />, Sergey V. Loiko <fnref n="4" />, and Oleg S. Pokrovsky <fnref n="3" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Aksenov, A.S.; Shirokova, L.S.; Kisil, O.Y.; Kolesova, S.N.; Lim, A.G.; Kuzmina, D.; Pouillé, S.; Alexis, M.A.; Castrec-Rouelle, M.; Loiko, S.V.; et al. Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile. *Diversity* **2021**, *13*, 328. https://doi.org/10.3390/d13070328

Academic Editors: Michael Wink, Boris A. Levin, Yulia V. Bespalaya and Maria Teresa Ceccherini Guicciardini

Received: 4 June 2021

Accepted: 14 July 2021

Published: 17 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Permafrost peatlands, containing a sizable amount of soil organic carbon (OC), play a pivotal role in soil (peat) OC transformation into soluble and volatile forms and greatly contribute to overall natural CO₂ and CH₄ emissions to the atmosphere under ongoing permafrost thaw and soil OC degradation. Peat microorganisms are largely responsible for the processing of this OC, yet coupled studies of chemical and bacterial parameters in permafrost peatlands are rather limited and geographically biased. Towards testing the possible impact of peat and peat pore water chemical composition on microbial population and diversity, here we present results of a preliminary study of the western Siberia permafrost peatland discontinuous permafrost zone. The quantitative evaluation of microorganisms and determination of microbial diversity along a 100 cm thick peat soil column, which included thawed and frozen peat and bottom mineral horizon, was performed by RT-PCR and 16S rRNA gene-based metagenomic analysis, respectively. Bacteria (mainly *Proteobacteria*, *Acidobacteria*, *Actinobacteria*) strongly dominated the microbial diversity (99% sequences), with a negligible proportion of archaea (0.3–0.5%). There was a systematic evolution of main taxa according to depth, with a maximum of 65% (*Acidobacteria*) encountered in the active layer, or permafrost boundary (50–60 cm). We also measured C, N, nutrients and ~50 major and trace elements in peat (19 samples) as well as its pore water and dispersed ice (10 samples), sampled over the same core, and we analyzed organic matter quality in six organic and one mineral horizon of this core. Using multiparametric statistics (PCA), we tested the links between the total microbial number and 16S rRNA diversity and chemical composition of both the solid and fluid phase harboring the microorganisms. Under climate warming and permafrost thaw, one can expect a downward movement of the layer of maximal genetic diversity following the active layer thickening. Given a one to two orders of magnitude higher microbial number in the upper (thawed) layers compared to bottom (frozen) layers, an additional 50 cm of peat thawing in western Siberia may sizably increase the total microbial population and biodiversity of active cells.

**Keywords:** permafrost; peat; pore water; bacteria; number; metabolism; active layer

<footnote n="1">¹ Arctic Biomonitoring Laboratory, Northern (Arctic) Federal University, 17 Naberezhnaya Sev. Dviny, 163002 Arkhangelsk, Russia; kisil.oksana.jroslavovna@gmail.com (O.Y.K.); sofianvanwyn@gmail.com (S.N.K.)</footnote>

<footnote n="2">² N. Laverov Federal Center for Integrated Arctic Research, Russian Academy of Sciences, 23 Naberezhnaya Sev. Dviny, 163000 Arkhangelsk, Russia; liudmila.shirokova@get.omp.eu</footnote>

<footnote n="3">³ Geosciences and Environment Toulouse, UMR 5563 CNRS, 14 Avenue Edouard Belin, 31400 Toulouse, France; oleg.pokrovsky@get.omp.eu</footnote>

<footnote n="4">⁴ BIO-GEO-CLIM Laboratory, Tomsk State University, 36 Lenina, Build. 1, 634050 Tomsk, Russia; lim_artyom@mail.ru (A.G.L.); kuzmina.d.m.95@gmail.com (D.K.); s.loyko@yandex.ru (S.V.L.)</footnote>

<footnote n="5">5 UMR 7619 Métis, Sorbonne Université, F-75252 Paris, France; 24sophiepouille@gmail.com (S.P.); marie.alexis@sorbonne-universite.fr (M.A.A.); maryse.rouelle@sorbonne-universite.fr (M.C.-R.)</footnote>

<footnote>* Correspondence: a.s.aksenov@narfu.ru</footnote>
# 1. Introduction

In contrast to the extensive research centered on bacterial abundance, diversity and metabolic activity in permanently and seasonally frozen mineral soils from high-latitude Canadian High Arctic [1–3], Spitsbergen [4], Siberia [5–7], Alaskan Arctic [8,9], Greenland [10] and Antarctica [11–13], frozen peat (organic) environments remain poorly characterized in terms of the physiological diversity and the metabolic potential of bacteria. Permafrost peatlands represent the largest reservoir of currently frozen soil organic carbon, susceptible to microbial degradation under ongoing permafrost thaw and the increase of active layer thickness (ALT). For example, the western Siberia Lowland (WSL) exhibits a typical peat thickness of 1 to 3 m [14] and a maximal active layer depth of 30 to 80 cm [15–19]. Here, a significant reservoir of organic carbon and nitrogen is currently frozen but may become available for microbial processing in the next century. However, the microbial metabolic activity and even the number of bacteria in the thawed versus frozen parts of the peat soil cores are only beginning to be explored ([20] and references therein).

Over the past decade, several studies (i.e., [1,3,21–29]) have addressed microbial communities along a soil column that includes both seasonally thawed and permanently frozen horizons; see [30] for a review. However, these studies dealt with samples that contained a relatively small amount of organic carbon (<5–10% of $C_{org}$) and therefore characterized the bacterial activity and diversity in the mineral rather than the organic frozen layer. Much less is known of the microbial diversity in peatland soils subjected to permafrost impact ([31–33]). Moreover, some studies have documented significant differences in the microbial community structure between different soil horizons, including in the organic topsoil of western Siberia [34] and the subarctic soil of Alaska [23,35] and the Canadian High Arctic [1]. A few studies have used a metagenomic approach to address the diversity of microbial communities in arctic peat [31,33,36,37] and mineral [38–41] soils linked to the methane cycle [42,43], and several authors have characterized microbial biogeography in surface waters of permafrost regions [44–47].

However, despite this extensive work on microbial diversity in permafrost environments, the main external (physico-chemical) parameters governing microbial populations remain poorly understood. In this work, we aimed at testing the functional link between microbial diversity along various layers of a peat core and the organic and inorganic chemical composition of peat along the same core. The novelty of the present study consists of relating the genetic diversity of microorganisms to the organic matter quality of the peat core and elementary composition of both solid and liquid compartments of the core sampled during the period of maximal active layer thickness; this includes both the thawed layer and frozen organic horizons, down to the mineral horizon.

In agreement with previous studies [3,4,48–52], we hypothesized a decrease in both the number and genetic diversity of bacteria with the increase in soil depth, with an abrupt drop below the permafrost boundary. We further hypothesized, based on a previous assessment of metabolic diversity of microorganisms along the peat core of the WSL [20], that there will be a functional control, due to the organic chemical composition and trace metal concentration, of the relative abundance of dominant taxa over the surface and deep horizons, with a possible maximum of specific groups at the boundary of the active-permafrost layer. To test these hypotheses, the specific goals of this work are as follows: (i) to quantify the abundance and genetic diversity of bacteria and archaea over the full depth of the peat core, including the underlying mineral horizon and (ii) to identify potential organic substrates and trace elements that may be controlling microbial number and diversity in both the frozen and unfrozen parts of the peat core. Achieving these goals should allow us to predict the possible evolution of the soil microbial community due to downward migration of the active layer boundary, including its reaching of underlying mineral horizons due to ongoing permafrost thaw in western Siberia. This should eventually allow upscaling of these results to other permafrost peatlands of northern Eurasia.
# 1. Introduction

In contrast to the extensive research centered on bacterial abundance, diversity and metabolic activity in permanently and seasonally frozen mineral soils from high-latitude Canadian High Arctic [1–3], Spitsbergen [4], Siberia [5–7], Alaskan Arctic [8,9], Greenland [10] and Antarctica [11–13], frozen peat (organic) environments remain poorly characterized in terms of the physiological diversity and the metabolic potential of bacteria. Permafrost peatlands represent the largest reservoir of currently frozen soil organic carbon, susceptible to microbial degradation under ongoing permafrost thaw and the increase of active layer thickness (ALT). For example, the western Siberia Lowland (WSL) exhibits a typical peat thickness of 1 to 3 m [14] and a maximal active layer depth of 30 to 80 cm [15–19]. Here, a significant reservoir of organic carbon and nitrogen is currently frozen but may become available for microbial processing in the next century. However, the microbial metabolic activity and even the number of bacteria in the thawed versus frozen parts of the peat soil cores are only beginning to be explored ([20] and references therein).

Over the past decade, several studies (i.e., [1,3,21–29]) have addressed microbial communities along a soil column that includes both seasonally thawed and permanently frozen horizons; see [30] for a review. However, these studies dealt with samples that contained a relatively small amount of organic carbon (<5–10% of $C_{org}$) and therefore characterized the bacterial activity and diversity in the mineral rather than the organic frozen layer. Much less is known of the microbial diversity in peatland soils subjected to permafrost impact ([31–33]). Moreover, some studies have documented significant differences in the microbial community structure between different soil horizons, including in the organic topsoil of western Siberia [34] and the subarctic soil of Alaska [23,35] and the Canadian High Arctic [1]. A few studies have used a metagenomic approach to address the diversity of microbial communities in arctic peat [31,33,36,37] and mineral [38–41] soils linked to the methane cycle [42,43], and several authors have characterized microbial biogeography in surface waters of permafrost regions [44–47].

However, despite this extensive work on microbial diversity in permafrost environments, the main external (physico-chemical) parameters governing microbial populations remain poorly understood. In this work, we aimed at testing the functional link between microbial diversity along various layers of a peat core and the organic and inorganic chemical composition of peat along the same core. The novelty of the present study consists of relating the genetic diversity of microorganisms to the organic matter quality of the peat core and elementary composition of both solid and liquid compartments of the core sampled during the period of maximal active layer thickness; this includes both the thawed layer and frozen organic horizons, down to the mineral horizon.

In agreement with previous studies [3,4,48–52], we hypothesized a decrease in both the number and genetic diversity of bacteria with the increase in soil depth, with an abrupt drop below the permafrost boundary. We further hypothesized, based on a previous assessment of metabolic diversity of microorganisms along the peat core of the WSL [20], that there will be a functional control, due to the organic chemical composition and trace metal concentration, of the relative abundance of dominant taxa over the surface and deep horizons, with a possible maximum of specific groups at the boundary of the active-permafrost layer. To test these hypotheses, the specific goals of this work are as follows: (i) to quantify the abundance and genetic diversity of bacteria and archaea over the full depth of the peat core, including the underlying mineral horizon and (ii) to identify potential organic substrates and trace elements that may be controlling microbial number and diversity in both the frozen and unfrozen parts of the peat core. Achieving these goals should allow us to predict the possible evolution of the soil microbial community due to downward migration of the active layer boundary, including its reaching of underlying mineral horizons due to ongoing permafrost thaw in western Siberia. This should eventually allow upscaling of these results to other permafrost peatlands of northern Eurasia.
## 2. Study Site and Methods

### 2.1. Peat Core Collection and Temperature Pattern of the WSL Peat Profile

Peat core samples were collected on 5 August 2017 in the palsa peat bog at the watershed divide of the discontinuous permafrost zone near the Khanymey test site of the northern part of the West Siberian Lowland (Figure 1). A detailed description of the study site is given elsewhere [20]. The dominant soil catenae are flat-mound bogs along the watershed divides, which comprise the Hemic Cryic Histosols on mounds. There, the peat (1.0–1.5 m thick) has developed on sands with minor amounts of clay and silt. The mounds (50 ± 20% of total palsa area) and depressions are less than 1 m in diameter, and the depth of the depressions is 40 to 50 cm. The average thickness of the active layer is 41 ± 5 cm on the mounds.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Ground view of the study site (63°47′50″ N, 10″, 75°33′21,90″ E), a photo of the soil pit and environmental context of the frozen palsa peat bog and a schematic profile of the peat column with relevant horizons. Photo made by Loiko S.V.</figcaption></figure>

Peat core samples from mounds were extracted using a motorized Russian peat corer (UKB-12/25 I, Saint Petersburg, Russia) with a 4 cm diameter corer, sterilized with 40% ethanol prior to each extraction. Core subsections were taken each 5–10 cm over the length of the peat core. The samples were kept at –20 °C during transport and storage. Sterility of the part of the core collected for microbiological analysis was achieved using peat exclusively from the interior of the core (>1 cm from the core liner), following conventional procedures [20,53]. At every 5–10 cm of depth, soil samples were taken with a sharp,
**2. Study Site and Methods**

**2.1. Peat Core Collection and Temperature Pattern of the WSL Peat Profile**

Peat core samples were collected on 5 August 2017 in the palsa peat bog at the watershed divide of the discontinuous permafrost zone near the Khanymey test site of the northern part of the West Siberian Lowland (Figure 1). A detailed description of the study site is given elsewhere [20]. The dominant soil catenae are flat-mound bogs along the watershed divides, which comprise the Hemic Cryic Histosols on mounds. There, the peat (1.0–1.5 m thick) has developed on sands with minor amounts of clay and silt. The mounds (50 ± 20% of total palsa area) and depressions are less than 1 m in diameter, and the depth of the depressions is 40 to 50 cm. The average thickness of the active layer is 41 ± 5 cm on the mounds.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Ground view of the study site (63°47′50″ N,10″, 75°33′21,90″ E), a photo of the soil pit and environmental context of the frozen palsa peat bog and a schematic profile of the peat column with relevant horizons. Photo made by Loiko S.V.</figcaption></figure>

Peat core samples from mounds were extracted using a motorized Russian peat corer (UKB-12/25 I, Saint Petersburg, Russia) with a 4 cm diameter corer, sterilized with 40% ethanol prior to each extraction. Core subsections were taken each 5–10 cm over the length of the peat core. The samples were kept at –20 °C during transport and storage. Sterility of the part of the core collected for microbiological analysis was achieved using peat exclusively from the interior of the core (>1 cm from the core liner), following conventional procedures [20,53]. At every 5–10 cm of depth, soil samples were taken with a sharp,
sterile, single-use knife and placed in sterile PVC doubled-zipped bags. All subsequent
subsampling was performed under sterile, protected conditions in the laminar hood box
(A100).

To collect peat ice, the frozen peat core was cut into 10 cm layers using a ceramic blade and thawed at room temperature. Immediately after full thawing of the frozen core, peat water was extracted through 100 µm pre-cleaned nylon net by applying a local confining pressure (1–5 kg cm⁻², depending on sample humidity). The obtained slurry was centrifuged for 15 min at 3500 rpm, and the supernatant was filtered through a 0.45 µm acetate cellulose filter, with the first 10 mL of filtrate discarded. The water from the thawed part of the peat core, above the permafrost table, was extracted using the same method but in the field and immediately after peat sampling.

The degree of peat decomposition was determined as described previously [20]. The peat is essentially composed of *Sphagnum rubellum*, lichens (*Cladonia stellaris*), dwarf shrubs, green mosses, pinewood, and pine and birch bark, as described in detail elsewhere [20]. The peat color was measured in color coordinates CIE lab using a Rite VS450 (X-Rite, Inc., Grand Rapids, MI, USA).

The radiocarbon age of peat was determined by the liquid scintillation method using a Quantulus 1220 radiometer (Wallac, Turku, Finland) in the laboratory of the Tomsk Collective Use Center (Russian Academy of Sciences). Calibration of the radiocarbon age in the calendar was conducted using the OxCal online program and the IntCal 20 calibration curve [54].

The annual pattern of temperature at different depths was recorded each 10 cm over
a 0–1 m depth via in situ temperature sensors CAM-2 (IMCES SB RAS, Tomsk, Russia),
installed 2 years before sampling of the peat core.

## 2.2. *Total Chemical Analysis*

The peat core was analyzed for concentrations of C, N, P, and K as well as major and trace elements. Nitrogen and carbon concentrations in dry peat samples were measured by Cu-O catalyzed dry combustion at 900 °C with ≤0.5% precision for standard substances (Thermo Flash 2000 CN Analyzer). Sulfur was measured by EMIA-Expert Carbon/Sulfur Analyzer (Horiba, Kyoto, Japan). For the measurement of the element concentrations, the samples were first processed in a clean room (class A 10,000). Approximately 100 mg of peat, previously ground with an agate mortar and pestle to a <1 mm fraction, was placed in Teflon (Savilex®, Eden Prairie, MN, USA) reactors with 6 mL bi-distilled HNO₃, 0.2 mL ultrapure HF and 1 mL ultrapure H₂O₂. Twelve reactors were loaded into a Mars 5 microwave digestion system (CEM, Paris, France) and exposed for 20 min at 150 °C. Each series of reactors was composed of 10 samples of peat, 1 certified lichen standard CRM 482 sample (from BCR, Brussels, Belgium) or other NIST soil standards (NIST 271a Montana and NIST 2709a San Joaquin) and 1 blank sample. The major and trace element concentrations were measured by inductively coupled plasma mass spectrometry (ICP-MS; Agilent 7500 ce) using a three-point calibration against a standard solution of known concentration. Indium and rhenium were used as the internal standards to correct for instrumental drift and eventual matrix effects. Further details of elementary analyses of peat are available in Stepanova et al. [55]. The elementary composition of peat was analyzed with a 5 cm step; the results from individual layers were averaged to match the 10 cm spatial resolution of microbiological analysis. In addition to total chemical analysis, labile P-PO₄ in aqueous leachates was measured using the colorimetric technique [56].

Filtered peat pore water from the active layer and peat ice were analyzed as described elsewhere [57,58]. In brief, filtered solutions were divided into two parts: (1) acidified (pH ~2) with ultrapure double-distilled HNO₃ for trace element (TE) analyses and (2) non-acidified for dissolved organic carbon (DOC), dissolved inorganic carbon (DIC) and anion analyses. DOC and DIC were analyzed using a Carbon Total Analyzer (Shimadzu TOC VSCN) with an uncertainty better than 3%. Major anion (Cl, SO₄²⁻) concentrations were measured by ion chromatography (HPLC, Dionex ICS 2000), with an uncertainty of 2%.
sterile, single-use knife and placed in sterile PVC doubled-zipped bags. All subsequent
subsampling was performed under sterile, protected conditions in the laminar hood box
(A100).

To collect peat ice, the frozen peat core was cut into 10 cm layers using a ceramic blade and thawed at room temperature. Immediately after full thawing of the frozen core, peat water was extracted through 100 µm pre-cleaned nylon net by applying a local confining pressure (1–5 kg cm⁻², depending on sample humidity). The obtained slurry was centrifuged for 15 min at 3500 rpm, and the supernatant was filtered through a 0.45 µm acetate cellulose filter, with the first 10 mL of filtrate discarded. The water from the thawed part of the peat core, above the permafrost table, was extracted using the same method but in the field and immediately after peat sampling.

The degree of peat decomposition was determined as described previously [20]. The peat is essentially composed of *Sphagnum rubellum*, lichens (*Cladonia stellaris*), dwarf shrubs, green mosses, pinewood, and pine and birch bark, as described in detail elsewhere [20]. The peat color was measured in color coordinates CIE lab using a Rite VS450 (X-Rite, Inc., Grand Rapids, MI, USA).

The radiocarbon age of peat was determined by the liquid scintillation method using a Quantulus 1220 radiometer (Wallac, Turku, Finland) in the laboratory of the Tomsk Collective Use Center (Russian Academy of Sciences). Calibration of the radiocarbon age in the calendar was conducted using the OxCal online program and the IntCal 20 calibration curve [54].

The annual pattern of temperature at different depths was recorded each 10 cm over
a 0-1 m depth via in situ temperature sensors CAM-2 (IMCES SB RAS, Tomsk, Russia),
installed 2 years before sampling of the peat core.

## 2.2. *Total Chemical Analysis*

The peat core was analyzed for concentrations of C, N, P, and K as well as major and trace elements. Nitrogen and carbon concentrations in dry peat samples were measured by Cu-O catalyzed dry combustion at 900 °C with ≤0.5% precision for standard substances (Thermo Flash 2000 CN Analyzer). Sulfur was measured by EMIA-Expert Carbon/Sulfur Analyzer (Horiba, Kyoto, Japan). For the measurement of the element concentrations, the samples were first processed in a clean room (class A 10,000). Approximately 100 mg of peat, previously ground with an agate mortar and pestle to a <1 mm fraction, was placed in Teflon (Savilex®, Eden Prairie, MN, USA) reactors with 6 mL bi-distilled HNO₃, 0.2 mL ultrapure HF and 1 mL ultrapure H₂O₂. Twelve reactors were loaded into a Mars 5 microwave digestion system (CEM, Paris, France) and exposed for 20 min at 150 °C. Each series of reactors was composed of 10 samples of peat, 1 certified lichen standard CRM 482 sample (from BCR, Brussels, Belgium) or other NIST soil standards (NIST 271a Montana and NIST 2709a San Joaquin) and 1 blank sample. The major and trace element concentrations were measured by inductively coupled plasma mass spectrometry (ICP-MS; Agilent 7500 ce) using a three-point calibration against a standard solution of known concentration. Indium and rhenium were used as the internal standards to correct for instrumental drift and eventual matrix effects. Further details of elementary analyses of peat are available in Stepanova et al. [55]. The elementary composition of peat was analyzed with a 5 cm step; the results from individual layers were averaged to match the 10 cm spatial resolution of microbiological analysis. In addition to total chemical analysis, labile P-PO₄ in aqueous leachates was measured using the colorimetric technique [56].

Filtered peat pore water from the active layer and peat ice were analyzed as described elsewhere [57,58]. In brief, filtered solutions were divided into two parts: (1) acidified (pH ~2) with ultrapure double-distilled HNO₃ for trace element (TE) analyses and (2) non-acidified for dissolved organic carbon (DOC), dissolved inorganic carbon (DIC) and anion analyses. DOC and DIC were analyzed using a Carbon Total Analyzer (Shimadzu TOC VSCN) with an uncertainty better than 3%. Major anion (Cl, SO₄²⁻) concentrations were measured by ion chromatography (HPLC, Dionex ICS 2000), with an uncertainty of 2%.
Major cations (Ca, Mg, Na, K), Si and ~40 trace elements were determined with an Agilent
ce 7500 ICP-MS, with In and Re as internal standards.

2.3. *Organic Matter Isotopic Composition and Quality*

In freeze-dried samples of different layers from the peat core, the δ¹³C (expressed
in ‰ of international PDB standard) and the δ¹⁵N (expressed in ‰ of atmospheric N₂)
were measured using an isotopic ratio mass spectrometer IRMS (Micromass Iso Prime) in
the iEES, Paris, France. An analytical standard (tyrosine) was measured for each of the
4 samples. The analytical precision was 0.2 and 0.1‰ for the δ¹³C and δ¹⁵N, respectively.

The organic matter of samples was characterized by Nuclear Magnetic Resonance (NMR) of solid-state ¹³C using a Bruker Avance 500 NMR (Bruker Analytik GmbH, Rheinstetten, Germany) in IMPC (Sorbonne Université, Paris, France), with a resonance frequency of 125 MHz for ¹³C. A CP-MAS (cross-polarization magic angle spinning) was used, with a rotation speed of 14 KHz. The contact and relaxation time were 1 ms and 1 s, respectively. The chemical shift was quantified with respect to trimethylsilane (TMS) and adjusted with adamantane (1.85 ppm for ¹H and 38.52 ppm for ¹³C) as external references. Spectra deconvolution was performed with the Dmfit 2015 program [59]. The spectra were divided into 4 regions of chemical shift, as defined by [60]: alkyl C (0–45 ppm), O-alkyl C (45–110 ppm), aromatic C (110–160 ppm) and carboxylic C (160–210 ppm). The proportion of these functional groups, after deconvolution and integration, was obtained via adding the areas of peaks corresponding to the defined regions.

2.4. *Analysis of Bacterial Number and Genetic Diversity*

The extraction of genomic DNA from frozen peat samples was conducted at the Research Institute for Agricultural Microbiology (Pushkin, St.-Petersburg, Russia) using the Macherey-Nagel nucleic spin kit (Macherey-Nagel, Düren, Germany), following the manufacturer's instructions. Taxonomic analysis of the bacterial community was carried out using universal primers F515 (GTGCCAGCMGCCGCGGTAA) and R806 (GGAC-TACVSGGGTATCTAAT) on the variable region V3–V4 of the 16S rRNA gene, which is specific for a wide range of microorganisms, including bacteria and archaea. A rough estimate of the number of microorganisms (counting clones) in the samples was carried out using real-time PCR using the following primers: Eub338/Eub518 for bacteria and arc915f/arc1059r for archaea. As a control for bacteria, cloned fragments of the ribosomal operon *Escherichia coli* (Sigma) were used, and for the archaea, strain FG-07 *Halobacterium salinarum* (the strain was provided by G. Jurgens, University of Helsinki, Helsinki, Finland). The determination for each sample was carried out in triplicate. Quantitative estimates are given as to the number of 16S rRNA gene copies per gram of soil.

2.5. *Statistical Treatment*

Statistical treatment of the data included best-fit functions based on the Pearson
correlation and one-way ANOVA with the STATISTICA version 8 software (StatSoft Inc.,
Tulsa, OK, USA). Linear regressions were used to examine the relationships among the
total and culturable cell numbers, the total average well color development (AWCD) and
the consumption of an individual group of substrates, and the chemical composition of
the peat layer (C, N, major and trace elements). The criterion for a significant correlation
between chemical and biological parameters was that the Pearson coefficients were higher
than 0.5 (R ≥ 0.5 and *p* < 0.05). The ANOVA method was used to test the differences in the
cell number, elementary composition, organic matter quality parameters and shares of the
16S rRNA gene fragments of individual phylogenetic groups. Further statistical treatment
of a complete set of a share of phylogenetic groups, inorganic (ICP-MS) and organic (NMR)
chemical and isotopic composition and physical parameters of peat layer included principal
components analysis (PCA). The PCA allowed testing of the effect of various parameters,
the differences between peat layers in particular, on the chemical and physical parameters
of peat and spatial distribution of microorganisms. Note that the selected peat core is highly
Major cations (Ca, Mg, Na, K), Si and ~40 trace elements were determined with an Agilent
ce 7500 ICP-MS, with In and Re as internal standards.

2.3. *Organic Matter Isotopic Composition and Quality*

In freeze-dried samples of different layers from the peat core, the δ¹³C (expressed in ‰ of international PDB standard) and the δ¹⁵N (expressed in ‰ of atmospheric N₂) were measured using an isotopic ratio mass spectrometer IRMS (Micromass Iso Prime) in the iEES, Paris, France. An analytical standard (tyrosine) was measured for each of the 4 samples. The analytical precision was 0.2 and 0.1‰ for the δ¹³C and δ¹⁵N, respectively.

The organic matter of samples was characterized by Nuclear Magnetic Resonance (NMR) of solid-state ¹³C using a Bruker Avance 500 NMR (Bruker Analytik GmbH, Rheinstetten, Germany) in IMPC (Sorbonne Université, Paris, France), with a resonance frequency of 125 MHz for ¹³C. A CP-MAS (cross-polarization magic angle spinning) was used, with a rotation speed of 14 KHz. The contact and relaxation time were 1 ms and 1 s, respectively. The chemical shift was quantified with respect to trimethylsilane (TMS) and adjusted with adamantane (1.85 ppm for ¹H and 38.52 ppm for ¹³C) as external references. Spectra deconvolution was performed with the Dmfit 2015 program [59]. The spectra were divided into 4 regions of chemical shift, as defined by [60]: alkyl C (0–45 ppm), O-alkyl C (45–110 ppm), aromatic C (110–160 ppm) and carboxylic C (160–210 ppm). The proportion of these functional groups, after deconvolution and integration, was obtained via adding the areas of peaks corresponding to the defined regions.

2.4. *Analysis of Bacterial Number and Genetic Diversity*

The extraction of genomic DNA from frozen peat samples was conducted at the Research Institute for Agricultural Microbiology (Pushkin, St.-Petersburg, Russia) using the Macherey-Nagel nucleic spin kit (Macherey-Nagel, Düren, Germany), following the manufacturer's instructions. Taxonomic analysis of the bacterial community was carried out using universal primers F515 (GTGCCAGCMGCCGCGGTAA) and R806 (GGAC-TACVSGGGTATCTAAT) on the variable region V3–V4 of the 16S rRNA gene, which is specific for a wide range of microorganisms, including bacteria and archaea. A rough estimate of the number of microorganisms (counting clones) in the samples was carried out using real-time PCR using the following primers: Eub338/Eub518 for bacteria and arc915f/arc1059r for archaea. As a control for bacteria, cloned fragments of the ribosomal operon *Escherichia coli* (Sigma) were used, and for the archaea, strain FG-07 *Halobacterium salinarum* (the strain was provided by G. Jurgens, University of Helsinki, Helsinki, Finland). The determination for each sample was carried out in triplicate. Quantitative estimates are given as to the number of 16S rRNA gene copies per gram of soil.

2.5. *Statistical Treatment*

Statistical treatment of the data included best-fit functions based on the Pearson
correlation and one-way ANOVA with the STATISTICA version 8 software (StatSoft Inc.,
Tulsa, OK, USA). Linear regressions were used to examine the relationships among the
total and culturable cell numbers, the total average well color development (AWCD) and
the consumption of an individual group of substrates, and the chemical composition of
the peat layer (C, N, major and trace elements). The criterion for a significant correlation
between chemical and biological parameters was that the Pearson coefficients were higher
than 0.5 (R ≥ 0.5 and *p* < 0.05). The ANOVA method was used to test the differences in the
cell number, elementary composition, organic matter quality parameters and shares of the
16S rRNA gene fragments of individual phylogenetic groups. Further statistical treatment
of a complete set of a share of phylogenetic groups, inorganic (ICP-MS) and organic (NMR)
chemical and isotopic composition and physical parameters of peat layer included principal
components analysis (PCA). The PCA allowed testing of the effect of various parameters,
the differences between peat layers in particular, on the chemical and physical parameters
of peat and spatial distribution of microorganisms. Note that the selected peat core is highly
representative of the region of WSL permafrost peatlands. Multiple cores sampled within the same test site of Khanymey (discontinuous permafrost zone), from different mounds, demonstrated a high similarity of the chemical composition of both pore water and peat ice [58]. Due to the high homogeneity of the peat and underlaying mineral substrate of the region, we used a single soil core for representing the chemical diversity of solid and liquid phases and the metabolic activity of microorganisms ([20]).

## 3. Results

### 3.1. Peat Formation History, Thermal Regime, Physical and Elementary Characteristics of the Peat Core and Chemical Composition of the Pore Water

Accumulation of peat started 8878 ± 74 BP (depth of 100 cm). The age of peat at the ALT boundary (41 cm) is 1886 ± 84 BP, whereas at a depth of 19 cm, the peat is 374 ± 89 years old. This corresponds to the rate of peat accumulation of 0.084 mm y⁻¹ for the depth of 41 to 100 cm and 0.15 mm y⁻¹ for the 19–41 cm layer. The upper 19 cm of peat was growing with a rate of 0.51 mm y⁻¹. This trend in accumulation rate increasing with time reflects the dominance in the upper part of the core of sphagnum mosses, which exhibit better preservation due to their slower degradation rate.

The thermal regime of the peat mound allowed for the determination of the average thickness of the active layer as 45–50 cm (Figure 1). The degree of peat decomposition increased from the top to the bottom of the sampled column, from 5% in the upper 10–20 cm layer to 30% below 70 cm depth, whereas the bulk peat density increased from 20 kg m⁻³ in the upper 0–2 cm layer to 200–250 kg m⁻³ below 60 cm depth.

The physical parameters of peat as well as the major and trace element compositions of the peat and peat pore water are listed in Table 1 and Table S1 of the Supplementary Materials. The C and N concentrations in the peat core are illustrated in Figure 2. Non-systematic variations in the C and N contents were observed over the peat core length, with an abrupt decrease in C and N in the mineral layer. The N<sub>tot</sub> concentration exhibited a local maximum at 75 cm depth. The C:N ratio varied between 50 and 100 in the active layer, decreased to 40–45 in the 40–85 cm layer, then increased to 63–64 in the bottom layers.

<table><caption>**Table 1.** Physical and chemical parameters of the studied core. The temperature (T) of the core was measured via in situ sensors over the whole year (2018–2019). The isotopic and RMN data corresponded to slightly different depths, as listed in Table 2, and they were averaged to match the depth horizons of genomic analysis.</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td>10–20</td><td>peat</td><td>-15</td><td>-0.047</td><td>-8.4</td><td>11.3</td><td>0.16</td><td>5</td></tr><tr><td>30–40</td><td>peat</td><td>-35</td><td>-0.63</td><td>-6.1</td><td>4.4</td><td>0.29</td><td>10</td></tr><tr><td>50–60</td><td>peat</td><td>-55</td><td>-0.8</td><td>-4.1</td><td>0.2</td><td>0.29</td><td>17</td></tr><tr><td>70–80</td><td>peat</td><td>-75</td><td>-0.74</td><td>-2.8</td><td>-0.1</td><td>0.23</td><td>30</td></tr><tr><td>90–100</td><td>peat</td><td>-95</td><td>-0.65</td><td>-2.3</td><td>-0.1</td><td>0.23</td><td>32</td></tr><tr><td>100–110</td><td>mineral</td><td>-105</td><td>-0.58</td><td>-1.3</td><td>-0.1</td><td>1.30</td><td></td></tr><tr><td>Depth, cm</td><td>substrate</td><td>N, %</td><td>C, %</td><td>C:N</td><td>L* by CIE lab</td><td>a* by CIE lab</td><td>b* by CIE lab</td></tr><tr><td>10–20</td><td>peat</td><td>0.79</td><td>44.6</td><td>56.1</td><td>43.7</td><td>9.9</td><td>16.5</td></tr><tr><td>30–40</td><td>peat</td><td>0.61</td><td>46.0</td><td>75.4</td><td>39.2</td><td>9.1</td><td>17.5</td></tr><tr><td>50–60</td><td>peat</td><td>0.83</td><td>40.6</td><td>48.9</td><td>33.2</td><td>8.7</td><td>17.5</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>N, %</td><td>C, %</td><td>C:N</td><td>L* by CIE lab</td><td>a* by CIE lab</td><td>b* by CIE lab</td></tr><tr><td>70–80</td><td>peat</td><td>1.36</td><td>56.1</td><td>41.3</td><td>25.9</td><td>7.8</td><td>14.0</td></tr><tr><td>90–100</td><td>peat</td><td>0.84</td><td>54.41</td><td>64.8</td><td>25.3</td><td>6.9</td><td>12.3</td></tr><tr><td>100–110</td><td>mineral</td><td>0.03</td><td>1.905</td><td>63.5</td><td>66.2</td><td>3.9</td><td>10.6</td></tr></tbody></table>
representative of the region of WSL permafrost peatlands. Multiple cores sampled within the same test site of Khanymey (discontinuous permafrost zone), from different mounds, demonstrated a high similarity of the chemical composition of both pore water and peat ice [58]. Due to the high homogeneity of the peat and underlaying mineral substrate of the region, we used a single soil core for representing the chemical diversity of solid and liquid phases and the metabolic activity of microorganisms ([20]).

## 3. Results

### 3.1. Peat Formation History, Thermal Regime, Physical and Elementary Characteristics of the Peat Core and Chemical Composition of the Pore Water

Accumulation of peat started 8878 ± 74 BP (depth of 100 cm). The age of peat at the ALT boundary (41 cm) is 1886 ± 84 BP, whereas at a depth of 19 cm, the peat is 374 ± 89 years old. This corresponds to the rate of peat accumulation of 0.084 mm y⁻¹ for the depth of 41 to 100 cm and 0.15 mm y⁻¹ for the 19–41 cm layer. The upper 19 cm of peat was growing with a rate of 0.51 mm y⁻¹. This trend in accumulation rate increasing with time reflects the dominance in the upper part of the core of sphagnum mosses, which exhibit better preservation due to their slower degradation rate.

The thermal regime of the peat mound allowed for the determination of the average thickness of the active layer as 45–50 cm (Figure 1). The degree of peat decomposition increased from the top to the bottom of the sampled column, from 5% in the upper 10–20 cm layer to 30% below 70 cm depth, whereas the bulk peat density increased from 20 kg m⁻³ in the upper 0–2 cm layer to 200–250 kg m⁻³ below 60 cm depth.

The physical parameters of peat as well as the major and trace element compositions of the peat and peat pore water are listed in Table 1 and Table S1 of the Supplementary Materials. The C and N concentrations in the peat core are illustrated in Figure 2. Non-systematic variations in the C and N contents were observed over the peat core length, with an abrupt decrease in C and N in the mineral layer. The N<sub>tot</sub> concentration exhibited a local maximum at 75 cm depth. The C:N ratio varied between 50 and 100 in the active layer, decreased to 40–45 in the 40–85 cm layer, then increased to 63–64 in the bottom layers.

<table><caption>**Table 1.** Physical and chemical parameters of the studied core. The temperature (T) of the core was measured via in situ sensors over the whole year (2018–2019). The isotopic and RMN data corresponded to slightly different depths, as listed in Table 2, and they were averaged to match the depth horizons of genomic analysis.</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td>10–20</td><td>peat</td><td>−15</td><td>−0.047</td><td>−8.4</td><td>11.3</td><td>0.16</td><td>5</td></tr><tr><td>30–40</td><td>peat</td><td>−35</td><td>−0.63</td><td>−6.1</td><td>4.4</td><td>0.29</td><td>10</td></tr><tr><td>50–60</td><td>peat</td><td>−55</td><td>−0.8</td><td>−4.1</td><td>0.2</td><td>0.29</td><td>17</td></tr><tr><td>70–80</td><td>peat</td><td>−75</td><td>−0.74</td><td>−2.8</td><td>−0.1</td><td>0.23</td><td>30</td></tr><tr><td>90–100</td><td>peat</td><td>−95</td><td>−0.65</td><td>−2.3</td><td>−0.1</td><td>0.23</td><td>32</td></tr><tr><td>100–110</td><td>mineral</td><td>−105</td><td>−0.58</td><td>−1.3</td><td>−0.1</td><td>1.30</td><td></td></tr><tr><td>Depth, cm</td><td>substrate</td><td>N, %</td><td>C, %</td><td>C:N</td><td>L* by CIE lab</td><td>a* by CIE lab</td><td>b* by CIE lab</td></tr><tr><td>10–20</td><td>peat</td><td>0.79</td><td>44.6</td><td>56.1</td><td>43.7</td><td>9.9</td><td>16.5</td></tr><tr><td>30–40</td><td>peat</td><td>0.61</td><td>46.0</td><td>75.4</td><td>39.2</td><td>9.1</td><td>17.5</td></tr><tr><td>50–60</td><td>peat</td><td>0.83</td><td>40.6</td><td>48.9</td><td>33.2</td><td>8.7</td><td>17.5</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>N, %</td><td>C, %</td><td>C:N</td><td>L* by CIE lab</td><td>a* by CIE lab</td><td>b* by CIE lab</td></tr><tr><td>70–80</td><td>peat</td><td>1.36</td><td>56.1</td><td>41.3</td><td>25.9</td><td>7.8</td><td>14.0</td></tr><tr><td>90–100</td><td>peat</td><td>0.84</td><td>54.41</td><td>64.8</td><td>25.3</td><td>6.9</td><td>12.3</td></tr><tr><td>100–110</td><td>mineral</td><td>0.03</td><td>1.905</td><td>63.5</td><td>66.2</td><td>3.9</td><td>10.6</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td>Depth, cm</td><td>substrate</td><td>S, %</td><td>P-lab, mg/kg</td><td>δ<sup>13</sup>C</td><td>δ<sup>15</sup>N</td><td>Alkyl C</td><td>O-alkyl C</td></tr><tr><td>10-20</td><td>peat</td><td>0.38</td><td>0.76</td><td>-24.7</td><td>4.4</td><td>24.05</td><td>71.9</td></tr><tr><td>30-40</td><td>peat</td><td>0.38</td><td>0.34</td><td>-25.9</td><td>4.8</td><td>15.3</td><td>77.3</td></tr><tr><td>50-60</td><td>peat</td><td>0.29</td><td>0.06</td><td>-26.4</td><td>4.25</td><td>26.2</td><td>63.9</td></tr><tr><td>70-80</td><td>peat</td><td>0.39</td><td>0.04</td><td>-26.9</td><td>3.7</td><td>37.1</td><td>50.5</td></tr><tr><td>90-100</td><td>peat</td><td>0.22</td><td>0.03</td><td>-26.2</td><td>4.1</td><td>47.1</td><td>39.3</td></tr><tr><td>100-110</td><td>mineral</td><td>0.05</td><td>0.00</td><td>-26.9</td><td></td><td></td><td></td></tr><tr><td>Depth, cm</td><td>substrate</td><td>aromatic C</td><td>carboxylic C</td><td>Li peat</td><td>B peat</td><td>Na peat</td><td>Mg peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10-20</td><td>peat</td><td>1.6</td><td>2.3</td><td>577</td><td>4518</td><td>372,422</td><td>296,970</td></tr><tr><td>30-40</td><td>peat</td><td>5.2</td><td>2.2</td><td>294</td><td>4544</td><td>103,287</td><td>327,233</td></tr><tr><td>50-60</td><td>peat</td><td>7.8</td><td>2.1</td><td>498</td><td>4434</td><td>72,538</td><td>194,873</td></tr><tr><td>70-80</td><td>peat</td><td>10.4</td><td>2.0</td><td>294</td><td>4892</td><td>29,929</td><td>169,643</td></tr><tr><td>90-100</td><td>peat</td><td>10.1</td><td>3.5</td><td>424</td><td>11,608</td><td>83,755</td><td>145,879</td></tr><tr><td>100-110</td><td>mineral</td><td></td><td></td><td>2441</td><td>3359</td><td>235,580</td><td>74,206</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Al peat</td><td>P peat</td><td>K peat</td><td>Ca peat</td><td>Ti peat</td><td>V peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10-20</td><td>peat</td><td>2,302,515</td><td>448,947</td><td>938,435</td><td>1,049,413</td><td>183,028</td><td>3059</td></tr><tr><td>30-40</td><td>peat</td><td>1,392,605</td><td>295,471</td><td>232,556</td><td>1,525,343</td><td>52,626</td><td>1199</td></tr><tr><td>50-60</td><td>peat</td><td>1,813,247</td><td>209,406</td><td>284,074</td><td>911,927</td><td>86,250</td><td>1987</td></tr><tr><td>70-80</td><td>peat</td><td>1,925,618</td><td>200,513</td><td>192,876</td><td>921,297</td><td>59,999</td><td>1936</td></tr><tr><td>90-100</td><td>peat</td><td>2,797,844</td><td>172,804</td><td>327,006</td><td>832,506</td><td>116,491</td><td>2942</td></tr><tr><td>100-110</td><td>mineral</td><td>1,748,772</td><td>22,357</td><td>991,020</td><td>279,474</td><td>218,869</td><td>1677</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Cr peat</td><td>Mn peat</td><td>Fe peat</td><td>Co peat</td><td>Ni peat</td><td>Cu peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10-20</td><td>peat</td><td>2504</td><td>10,234</td><td>1,311,464</td><td>281</td><td>1169</td><td>2138</td></tr><tr><td>30-40</td><td>peat</td><td>1966</td><td>5228</td><td>1,433,817</td><td>426</td><td>1083</td><td>594</td></tr><tr><td>50-60</td><td>peat</td><td>2745</td><td>5699</td><td>472,019</td><td>218</td><td>942</td><td>1152</td></tr><tr><td>70-80</td><td>peat</td><td>2793</td><td>2139</td><td>451,913</td><td>284</td><td>1841</td><td>2718</td></tr><tr><td>90-100</td><td>peat</td><td>3871</td><td>5225</td><td>519,426</td><td>522</td><td>1573</td><td>5064</td></tr><tr><td>100-110</td><td>mineral</td><td>1895</td><td>8516</td><td>290,073</td><td>177</td><td>263</td><td>914</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Zn peat</td><td>As peat</td><td>Rb peat</td><td>Sr peat</td><td>Y peat</td><td>Zr peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10-20</td><td>peat</td><td>15,946</td><td>1354</td><td>3086</td><td>9573</td><td>618</td><td>3614</td></tr><tr><td>30-40</td><td>peat</td><td>15,957</td><td>315</td><td>748</td><td>13,479</td><td>345</td><td>1356</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Zn peat</td><td>As peat</td><td>Rb peat</td><td>Sr peat</td><td>Y peat</td><td>Zr peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>50-60</td><td>peat</td><td>7367</td><td>261</td><td>1185</td><td>8224</td><td>463</td><td>2777</td></tr><tr><td>70-80</td><td>peat</td><td>9325</td><td>391</td><td>694</td><td>8406</td><td>332</td><td>1843</td></tr><tr><td>90-100</td><td>peat</td><td>17,238</td><td>380</td><td>1084</td><td>8394</td><td>414</td><td>5941</td></tr><tr><td>100-110</td><td>mineral</td><td>2355</td><td>299</td><td>2854</td><td>7518</td><td>913</td><td>14,104</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Nb peat</td><td>Mo peat</td><td>Cd peat</td><td>Sb peat</td><td>Cs peat</td><td>Ba peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10-20</td><td>peat</td><td>476</td><td>155</td><td>130</td><td>260</td><td>189</td><td>27,243</td></tr><tr><td>30-40</td><td>peat</td><td>462</td><td>87</td><td>132</td><td>39</td><td>46</td><td>20,066</td></tr><tr><td>50-60</td><td>peat</td><td>863</td><td>88</td><td>71</td><td>29</td><td>72</td><td>22,801</td></tr><tr><td>70-80</td><td>peat</td><td>579</td><td>211</td><td>86</td><td>42</td><td>41</td><td>29,938</td></tr><tr><td>90-100</td><td>peat</td><td>887</td><td>108</td><td>72</td><td>37</td><td>65</td><td>39,452</td></tr><tr><td>100-110</td><td>mineral</td><td>924</td><td>48</td><td>14</td><td>139</td><td>83</td><td>69,278</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td>Depth, cm</td><td>substrate</td><td>S, %</td><td>P-lab, mg/kg</td><td>δ<sup>13</sup>C</td><td>δ<sup>15</sup>N</td><td>Alkyl C</td><td>O-alkyl C</td></tr><tr><td>10–20</td><td>peat</td><td>0.38</td><td>0.76</td><td>-24.7</td><td>4.4</td><td>24.05</td><td>71.9</td></tr><tr><td>30–40</td><td>peat</td><td>0.38</td><td>0.34</td><td>-25.9</td><td>4.8</td><td>15.3</td><td>77.3</td></tr><tr><td>50–60</td><td>peat</td><td>0.29</td><td>0.06</td><td>-26.4</td><td>4.25</td><td>26.2</td><td>63.9</td></tr><tr><td>70–80</td><td>peat</td><td>0.39</td><td>0.04</td><td>-26.9</td><td>3.7</td><td>37.1</td><td>50.5</td></tr><tr><td>90–100</td><td>peat</td><td>0.22</td><td>0.03</td><td>-26.2</td><td>4.1</td><td>47.1</td><td>39.3</td></tr><tr><td>100–110</td><td>mineral</td><td>0.05</td><td>0.00</td><td>-26.9</td><td></td><td></td><td></td></tr><tr><td>Depth, cm</td><td>substrate</td><td>aromatic C</td><td>carboxylic C</td><td>Li peat</td><td>B peat</td><td>Na peat</td><td>Mg peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10–20</td><td>peat</td><td>1.6</td><td>2.3</td><td>577</td><td>4518</td><td>372,422</td><td>296,970</td></tr><tr><td>30–40</td><td>peat</td><td>5.2</td><td>2.2</td><td>294</td><td>4544</td><td>103,287</td><td>327,233</td></tr><tr><td>50–60</td><td>peat</td><td>7.8</td><td>2.1</td><td>498</td><td>4434</td><td>72,538</td><td>194,873</td></tr><tr><td>70–80</td><td>peat</td><td>10.4</td><td>2.0</td><td>294</td><td>4892</td><td>29,929</td><td>169,643</td></tr><tr><td>90–100</td><td>peat</td><td>10.1</td><td>3.5</td><td>424</td><td>11,608</td><td>83,755</td><td>145,879</td></tr><tr><td>100–110</td><td>mineral</td><td></td><td></td><td>2441</td><td>3359</td><td>235,580</td><td>74,206</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Al peat</td><td>P peat</td><td>K peat</td><td>Ca peat</td><td>Ti peat</td><td>V peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10–20</td><td>peat</td><td>2,302,515</td><td>448,947</td><td>938,435</td><td>1,049,413</td><td>183,028</td><td>3059</td></tr><tr><td>30–40</td><td>peat</td><td>1,392,605</td><td>295,471</td><td>232,556</td><td>1,525,343</td><td>52,626</td><td>1199</td></tr><tr><td>50–60</td><td>peat</td><td>1,813,247</td><td>209,406</td><td>284,074</td><td>911,927</td><td>86,250</td><td>1987</td></tr><tr><td>70–80</td><td>peat</td><td>1,925,618</td><td>200,513</td><td>192,876</td><td>921,297</td><td>59,999</td><td>1936</td></tr><tr><td>90–100</td><td>peat</td><td>2,797,844</td><td>172,804</td><td>327,006</td><td>832,506</td><td>116,491</td><td>2942</td></tr><tr><td>100–110</td><td>mineral</td><td>1,748,772</td><td>22,357</td><td>991,020</td><td>279,474</td><td>218,869</td><td>1677</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Cr peat</td><td>Mn peat</td><td>Fe peat</td><td>Co peat</td><td>Ni peat</td><td>Cu peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10–20</td><td>peat</td><td>2504</td><td>10,234</td><td>1,311,464</td><td>281</td><td>1169</td><td>2138</td></tr><tr><td>30–40</td><td>peat</td><td>1966</td><td>5228</td><td>1,433,817</td><td>426</td><td>1083</td><td>594</td></tr><tr><td>50–60</td><td>peat</td><td>2745</td><td>5699</td><td>472,019</td><td>218</td><td>942</td><td>1152</td></tr><tr><td>70–80</td><td>peat</td><td>2793</td><td>2139</td><td>451,913</td><td>284</td><td>1841</td><td>2718</td></tr><tr><td>90–100</td><td>peat</td><td>3871</td><td>5225</td><td>519,426</td><td>522</td><td>1573</td><td>5064</td></tr><tr><td>100–110</td><td>mineral</td><td>1895</td><td>8516</td><td>290,073</td><td>177</td><td>263</td><td>914</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Zn peat</td><td>As peat</td><td>Rb peat</td><td>Sr peat</td><td>Y peat</td><td>Zr peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10–20</td><td>peat</td><td>15,946</td><td>1354</td><td>3086</td><td>9573</td><td>618</td><td>3614</td></tr><tr><td>30–40</td><td>peat</td><td>15,957</td><td>315</td><td>748</td><td>13,479</td><td>345</td><td>1356</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Zn peat</td><td>As peat</td><td>Rb peat</td><td>Sr peat</td><td>Y peat</td><td>Zr peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>50–60</td><td>peat</td><td>7367</td><td>261</td><td>1185</td><td>8224</td><td>463</td><td>2777</td></tr><tr><td>70–80</td><td>peat</td><td>9325</td><td>391</td><td>694</td><td>8406</td><td>332</td><td>1843</td></tr><tr><td>90–100</td><td>peat</td><td>17,238</td><td>380</td><td>1084</td><td>8394</td><td>414</td><td>5941</td></tr><tr><td>100–110</td><td>mineral</td><td>2355</td><td>299</td><td>2854</td><td>7518</td><td>913</td><td>14,104</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Nb peat</td><td>Mo peat</td><td>Cd peat</td><td>Sb peat</td><td>Cs peat</td><td>Ba peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td></tr><tr><td>10–20</td><td>peat</td><td>476</td><td>155</td><td>130</td><td>260</td><td>189</td><td>27,243</td></tr><tr><td>30–40</td><td>peat</td><td>462</td><td>87</td><td>132</td><td>39</td><td>46</td><td>20,066</td></tr><tr><td>50–60</td><td>peat</td><td>863</td><td>88</td><td>71</td><td>29</td><td>72</td><td>22,801</td></tr><tr><td>70–80</td><td>peat</td><td>579</td><td>211</td><td>86</td><td>42</td><td>41</td><td>29,938</td></tr><tr><td>90–100</td><td>peat</td><td>887</td><td>108</td><td>72</td><td>37</td><td>65</td><td>39,452</td></tr><tr><td>100–110</td><td>mineral</td><td>924</td><td>48</td><td>14</td><td>139</td><td>83</td><td>69,278</td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td></td><td></td><td>La peat</td><td>Ce peat</td><td>Dy peat</td><td>Yb peat</td><td>Pb peat</td><td>U peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td colspan="3">ppb</td></tr><tr><td>10–20</td><td>peat</td><td>1184</td><td>2365</td><td>120</td><td>69</td><td>6485</td><td>80</td></tr><tr><td>30–40</td><td>peat</td><td>537</td><td>1167</td><td>72</td><td>32</td><td>1009</td><td>50</td></tr><tr><td>50–60</td><td>peat</td><td>775</td><td>1559</td><td>90</td><td>54</td><td>647</td><td>92</td></tr><tr><td>70–80</td><td>peat</td><td>534</td><td>1090</td><td>65</td><td>35</td><td>532</td><td>88</td></tr><tr><td>90–100</td><td>peat</td><td>727</td><td>1413</td><td>83</td><td>47</td><td>837</td><td>111</td></tr><tr><td>100–110</td><td>mineral</td><td>850</td><td>1498</td><td>108</td><td>86</td><td>1036</td><td>162</td></tr><tr><td>Depth, cm</td><td>substrate</td><td colspan="6">Humidity, %</td></tr><tr><td>10–20</td><td>peat</td><td colspan="6">317</td></tr><tr><td>30–40</td><td>peat</td><td colspan="6">337</td></tr><tr><td>50–60</td><td>peat</td><td colspan="6">700</td></tr><tr><td>70–80</td><td>peat</td><td colspan="6">458</td></tr><tr><td>90–100</td><td>peat</td><td colspan="6">650</td></tr><tr><td>100–110</td><td>mineral</td><td colspan="6">82</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Total C<sub>org</sub> (a) and N (b) in the peat core as a function of depth. The actual depth of the active layer on the day of sampling is 43 cm, shown by dashed line, and the mean maximal depth at the end of the active season is ~54 cm, shown by solid line. The concentrations of C and N in the peat and mineral layer are represented by triangles (upper X axis) and circles (bottom X axis), respectively.</figcaption></figure>

The upper 20 cm layer of peat was enriched in macro- (P, K, Mg) and micro-nutrients (Mn, Zn, Rb) compared to the rest of the peat core, without systematic variations in the depth below 20 cm (Table 1). Some trace elements of atmospheric origin (As, Cd, Sb, Cs, Tl, Pb) also enriched the upper 10–20 cm of the peat core. The other major and trace elements
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Depth, cm</th><th>Substrate</th><th>Average Depth, cm</th><th>Average T, °C</th><th>Minimal T, °C</th><th>Maximal T, °C</th><th>Peat, Bulk Density, g/cm³</th><th>Degree of Decomposition, %</th></tr></thead><tbody><tr><td></td><td></td><td>La peat</td><td>Ce peat</td><td>Dy peat</td><td>Yb peat</td><td>Pb peat</td><td>U peat</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td>ppb</td><td></td><td></td></tr><tr><td>10–20</td><td>peat</td><td>1184</td><td>2365</td><td>120</td><td>69</td><td>6485</td><td>80</td></tr><tr><td>30–40</td><td>peat</td><td>537</td><td>1167</td><td>72</td><td>32</td><td>1009</td><td>50</td></tr><tr><td>50–60</td><td>peat</td><td>775</td><td>1559</td><td>90</td><td>54</td><td>647</td><td>92</td></tr><tr><td>70–80</td><td>peat</td><td>534</td><td>1090</td><td>65</td><td>35</td><td>532</td><td>88</td></tr><tr><td>90–100</td><td>peat</td><td>727</td><td>1413</td><td>83</td><td>47</td><td>837</td><td>111</td></tr><tr><td>100–110</td><td>mineral</td><td>850</td><td>1498</td><td>108</td><td>86</td><td>1036</td><td>162</td></tr><tr><td>Depth, cm</td><td>substrate</td><td>Humidity, %</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>10–20</td><td>peat</td><td>317</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>30–40</td><td>peat</td><td>337</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>50–60</td><td>peat</td><td>700</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>70–80</td><td>peat</td><td>458</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>90–100</td><td>peat</td><td>650</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>100–110</td><td>mineral</td><td>82</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Total C<sub>org</sub> (a) and N (b) in the peat core as a function of depth. The actual depth of the active layer on the day of sampling is 43 cm, shown by dashed line, and the mean maximal depth at the end of the active season is ~54 cm, shown by solid line. The concentrations of C and N in the peat and mineral layer are represented by triangles (upper X axis) and circles (bottom X axis), respectively.</figcaption></figure>

The upper 20 cm layer of peat was enriched in macro- (P, K, Mg) and micro-nutrients (Mn, Zn, Rb) compared to the rest of the peat core, without systematic variations in the depth below 20 cm (Table 1). Some trace elements of atmospheric origin (As, Cd, Sb, Cs, Tl, Pb) also enriched the upper 10–20 cm of the peat core. The other major and trace elements
varied non-systematically across the sampled depth and did not demonstrate a significant (at *p* < 0.05) difference between the active and the frozen layer.

The peat pore water from the active layer and dispersed ice from the permafrost hori-
zons exhibited a vertical pattern of elementary composition that was different from that of
the solid phase (Figure 3). The phosphorus demonstrated a maximum at 70–80 cm, whereas
the majority of other elements exhibited a downward increase in concentration, with some
of them, such as micro-nutrients (B, K, Mn, Ni, Zn, Mo, Ba) and atmospherically-deposited
toxicants (As, Sb, Pb), also showing enrichment in the shallowest surface (0–10 cm) layer
(Table S2).

<figure><img src="image_6.png" /></figure>

**Figure 3.** Concentrations of DOC (a), specific UV absorbance at 254 nm (SUVA<sub>254</sub>) (b), P (c), Mn (d), Zn (e), Mo (f), As (g), Cd (h), and Pb (i) in the fluid phase of the peat core (pore water of the active layer and the permafrost ice). The actual depth of the active layer on the day of sampling is 43 cm, shown by dashed line, and the mean maximal depth at the end of the active season is ~54 cm, shown by solid line.
varied non-systematically across the sampled depth and did not demonstrate a significant (at *p* < 0.05) difference between the active and the frozen layer.

The peat pore water from the active layer and dispersed ice from the permafrost hori-
zons exhibited a vertical pattern of elementary composition that was different from that of
the solid phase (Figure 3). The phosphorus demonstrated a maximum at 70–80 cm, whereas
the majority of other elements exhibited a downward increase in concentration, with some
of them, such as micro-nutrients (B, K, Mn, Ni, Zn, Mo, Ba) and atmospherically-deposited
toxicants (As, Sb, Pb), also showing enrichment in the shallowest surface (0–10 cm) layer
(Table S2).

<figure><img src="image_6.png" /></figure>

**Figure 3.** Concentrations of DOC (a), specific UV absorbance at 254 nm (SUVA<sub>254</sub>) (b), P (c), Mn (d), Zn (e), Mo (f), As (g), Cd (h), and Pb (i) in the fluid phase of the peat core (pore water of the active layer and the permafrost ice). The actual depth of the active layer on the day of sampling is 43 cm, shown by dashed line, and the mean maximal depth at the end of the active season is ~54 cm, shown by solid line.
## 3.2. Organic Matter Isotopic Composition and Functional Group Analysis by NMR

The isotopic and organic matter (OM) quality characteristics are listed in Table 2. The $\delta^{13}C$ decreased significantly with depth ($p = 0.03$), showing a maximum at the 0–5 cm layer (–23.6‰) and exhibiting non-systematic variations between –26 and –27‰ below 60 cm depth. The $\delta^{15}N$ varied between +3 and +6‰, without any tendency with respect to depth. Based on Cross Polarization Magic-Angle Spinning (CP-MAS) solid-state $^{13}C$ NMR spectra, the ratio alkyl C:O-alkyl C and the proportion of aromatic C significantly increased with sample depth ($p = 0.05$ and 0.049, respectively). The proportion of carboxylic C (band at 170–175 ppm) over the sampled horizons is highly variable (between 1.2 and 3.5%) and did not show any trend with depth.

<table><caption>**Table 2.** Isotopic and NMR results of the CKH18-1T peat core. n.m. * stands for not measured.</caption><thead><tr><th rowspan="2">Depth</th><th>δ<sup>13</sup>C</th><th>δ<sup>15</sup>N</th><th>Alkyl C</th><th>O-Alkyl C</th><th>Aromatic C</th><th>Carbonyl C</th><th>Alkyl C: O-Alkyl C</th></tr><tr><th>‰</th><th>‰</th><th></th><th></th><th>%</th><th></th><th></th></tr></thead><tbody><tr><td>(0–5)</td><td>–23.6</td><td>2.8</td><td>16.7</td><td>79.8</td><td>0.0</td><td>3.5</td><td>0.21</td></tr><tr><td>(19–21)</td><td>–25.8</td><td>6.0</td><td>31.4</td><td>64.1</td><td>3.3</td><td>1.2</td><td>0.49</td></tr><tr><td>(30–35)</td><td>–25.9</td><td>4.8</td><td>15.3</td><td>77.3</td><td>5.2</td><td>2.2</td><td>0.20</td></tr><tr><td>(65–75)</td><td>–26.7</td><td>3.3</td><td>30.5</td><td>57.5</td><td>10.4</td><td>1.7</td><td>0.53</td></tr><tr><td>(75–85)</td><td>–27.1</td><td>4.1</td><td>43.7</td><td>43.5</td><td>10.4</td><td>2.4</td><td>1.00</td></tr><tr><td>(85–95)</td><td>–26.2</td><td>4.1</td><td>47.1</td><td>39.3</td><td>10.1</td><td>3.5</td><td>1.20</td></tr><tr><td>(105–115)</td><td>–26.9</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td></tr></tbody></table>

## 3.3. Cell Number and Microbial Diversity as a Function of Depth

The total number of microorganisms according to the real-time PCR, integrated over 0–110 cm depth, amounted to $2.26 \times 10^{10}$ (Table S2). Most of the peat microorganisms were present in the active layer, where environmental conditions are relatively milder than in the permafrost. The number of microorganisms decreased almost 1000-fold with depth, from $1.35 \times 10^{10}$ in the upper 0–20 cm layer to $2.27 \times 10^{7}$ in the mineral 100–110 cm layer.

Analysis of the composition of microbial communities by 16S rRNA sequencing demonstrated that the members of domain bacteria predominated in the soil profile (Figure 4). Based on the 16S rRNA gene fragments, the share of individual phylogenetic groups varied with soil depth (Table 3). The members of *Proteobacteria* (38.5%) and *Acidobacteria* (30.6%) predominated the sequence. The share of *Firmicutes* (3.25%) and *Actinobacteria* (14.1%) members was relatively high. Overall, the diversity of bacterial phyla gradually decreased with depth (Figure 4a). However, at the depth of 90–100 cm, the diversity sharply increased, but then it decreased again in the mineral layer (100–110 cm). The percentage of *Acidobacteria* and *Actinobacteria* increased until the depth of 50–60 cm, then decreased, and increased again at the depth of 90–100 cm (*Actinobacteria*) or in the mineral layer (*Acidobacteria*). The share of *Proteobacteria* decreased with depth, from 26.4% (10–20 cm) to 10.8% (50–60 cm). From the depth of 70–80 cm, the proportion of *Proteobacteria* began to rise again, up to 89.2% at the depth of 100–110 cm. The phylum *Verrucomicrobia* was present in significant quantities at two non-adjacent depths: 10.3% at 10–20 cm and 5.3% at 90–100 cm. The share of *Firmicutes* was insignificant up to a depth of 90 cm. In the deeper peat layers, the percentage of this phylum increased sharply to 15.9%, whereas in the underlaying mineral layer (100–110 cm), *Firmicutes* accounted for 2.5% of total phyla of microorganisms.
## 3.2. Organic Matter Isotopic Composition and Functional Group Analysis by NMR

The isotopic and organic matter (OM) quality characteristics are listed in Table 2. The $\delta^{13}C$ decreased significantly with depth ($p = 0.03$), showing a maximum at the 0–5 cm layer (–23.6‰) and exhibiting non-systematic variations between –26 and –27‰ below 60 cm depth. The $\delta^{15}N$ varied between +3 and +6‰, without any tendency with respect to depth. Based on Cross Polarization Magic-Angle Spinning (CP-MAS) solid-state $^{13}C$ NMR spectra, the ratio alkyl C:O-alkyl C and the proportion of aromatic C significantly increased with sample depth ($p = 0.05$ and 0.049, respectively). The proportion of carboxylic C (band at 170–175 ppm) over the sampled horizons is highly variable (between 1.2 and 3.5%) and did not show any trend with depth.

<table><caption>**Table 2.** Isotopic and NMR results of the CKH18-1T peat core. n.m. * stands for not measured.</caption><thead><tr><th rowspan="2">Depth</th><th>δ<sup>13</sup>C</th><th>δ<sup>15</sup>N</th><th>Alkyl C</th><th>O-Alkyl C</th><th>Aromatic C</th><th>Carbonyl C</th><th>Alkyl C: O-Alkyl C</th></tr><tr><th>‰</th><th>‰</th><th></th><th></th><th>%</th><th></th><th></th></tr></thead><tbody><tr><td>(0–5)</td><td>–23.6</td><td>2.8</td><td>16.7</td><td>79.8</td><td>0.0</td><td>3.5</td><td>0.21</td></tr><tr><td>(19–21)</td><td>–25.8</td><td>6.0</td><td>31.4</td><td>64.1</td><td>3.3</td><td>1.2</td><td>0.49</td></tr><tr><td>(30–35)</td><td>–25.9</td><td>4.8</td><td>15.3</td><td>77.3</td><td>5.2</td><td>2.2</td><td>0.20</td></tr><tr><td>(65–75)</td><td>–26.7</td><td>3.3</td><td>30.5</td><td>57.5</td><td>10.4</td><td>1.7</td><td>0.53</td></tr><tr><td>(75–85)</td><td>–27.1</td><td>4.1</td><td>43.7</td><td>43.5</td><td>10.4</td><td>2.4</td><td>1.00</td></tr><tr><td>(85–95)</td><td>–26.2</td><td>4.1</td><td>47.1</td><td>39.3</td><td>10.1</td><td>3.5</td><td>1.20</td></tr><tr><td>(105–115)</td><td>–26.9</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td><td>n.m *</td></tr></tbody></table>

## 3.3. Cell Number and Microbial Diversity as a Function of Depth

The total number of microorganisms according to the real-time PCR, integrated over 0–110 cm depth, amounted to $2.26 \times 10^{10}$ (Table S2). Most of the peat microorganisms were present in the active layer, where environmental conditions are relatively milder than in the permafrost. The number of microorganisms decreased almost 1000-fold with depth, from $1.35 \times 10^{10}$ in the upper 0–20 cm layer to $2.27 \times 10^{7}$ in the mineral 100–110 cm layer.

Analysis of the composition of microbial communities by 16S rRNA sequencing demonstrated that the members of domain bacteria predominated in the soil profile (Figure 4). Based on the 16S rRNA gene fragments, the share of individual phylogenetic groups varied with soil depth (Table 3). The members of *Proteobacteria* (38.5%) and *Acidobacteria* (30.6%) predominated the sequence. The share of *Firmicutes* (3.25%) and *Actinobacteria* (14.1%) members was relatively high. Overall, the diversity of bacterial phyla gradually decreased with depth (Figure 4a). However, at the depth of 90–100 cm, the diversity sharply increased, but then it decreased again in the mineral layer (100–110 cm). The percentage of *Acidobacteria* and *Actinobacteria* increased until the depth of 50–60 cm, then decreased, and increased again at the depth of 90–100 cm (*Actinobacteria*) or in the mineral layer (*Acidobacteria*). The share of *Proteobacteria* decreased with depth, from 26.4% (10–20 cm) to 10.8% (50–60 cm). From the depth of 70–80 cm, the proportion of *Proteobacteria* began to rise again, up to 89.2% at the depth of 100–110 cm. The phylum *Verrucomicrobia* was present in significant quantities at two non-adjacent depths: 10.3% at 10–20 cm and 5.3% at 90–100 cm. The share of *Firmicutes* was insignificant up to a depth of 90 cm. In the deeper peat layers, the percentage of this phylum increased sharply to 15.9%, whereas in the underlaying mineral layer (100–110 cm), *Firmicutes* accounted for 2.5% of total phyla of microorganisms.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Phylum-level (a) and class-level (b) bacterial and archaeal community structure in the peat core based on 16S rRNA.</figcaption></figure>

The phyla *Armatimonadetes* (0.2%), *Bacteroidetes* (1.6%), *Patescibacteria* (0.3%) and *Verru-
comicrobia* (2.8%) were less diverse and were represented by one class each. The members
of the class *Alphaproteobacteria* predominated and varied depending on the depth. The
largest share of this class (50.2%) was related to the depth of 70–80 cm, and then it de-
creased (Figure 4b). Classes *Acidobacteriia* (30.3%) and *Gammaproteobacteria* (20.4%) were
also represented by a large percentage. *Acidobacteriia* increased in number with depth
from 20.5% at the surface to 64.7% in the 50–60 cm layer, and then decreased. However,
in the mineral horizon (100–110 cm), *Acidobacteriia* were sizably higher (4.9%) than in the
bottom peat layer (90–100 cm, 1.0%). The 16S rRNA gene sequences belonging to the
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Phylum-level (a) and class-level (b) bacterial and archaeal community structure in the peat core based on 16S rRNA.</figcaption></figure>

The phyla *Armatimonadetes* (0.2%), *Bacteroidetes* (1.6%), *Patescibacteria* (0.3%) and *Verru-
comicrobia* (2.8%) were less diverse and were represented by one class each. The members
of the class *Alphaproteobacteria* predominated and varied depending on the depth. The
largest share of this class (50.2%) was related to the depth of 70–80 cm, and then it de-
creased (Figure 4b). Classes *Acidobacteriia* (30.3%) and *Gammaproteobacteria* (20.4%) were
also represented by a large percentage. *Acidobacteriia* increased in number with depth
from 20.5% at the surface to 64.7% in the 50–60 cm layer, and then decreased. However,
in the mineral horizon (100–110 cm), *Acidobacteriia* were sizably higher (4.9%) than in the
bottom peat layer (90–100 cm, 1.0%). The 16S rRNA gene sequences belonging to the
class *Gammaproteobacteria* were concentrated at the depth of 90–110 cm (17.7–87.5%). The class *Actinobacteria* exhibited sizable and non-systematic variations over depth, reaching the largest proportion (18.4%) at the depth of 90 to 100 cm. Classes *Bacteroidia* (6.1%), *Bacilli* (10.7%) and *Clostridia* (4.9%) were mostly present also at the same depth. Classes *Chthonomonadetes*, *Phycisphaerae*, *Planctomycetacia* and other non-identified bacteria predominated in the active layers (between 10 and 40 cm depth). Classes *Deltaproteobacteria* and *Verrucomicrobiae* were found in the greatest number at the surface and bottom horizon of the peat column: 2.9% and 10.3%, respectively, in the 10–20 cm layer; 1.6% and 5.3%, respectively, in the 90–100 cm layer. Class *Thermoleophilia* was almost uniformly distributed over the entire peat profile (Figure 4b). The classes *Gammaproteobacteria*, *Actinobacteria* and *Alphaproteobacteria* were most diverse and included six, five, and four orders, respectively. Classes *Acidobacteriia*, *Bacilli*, *Planctomycetacia* and *Verrucomicrobia* included two orders each. Other classes were represented only by one order.

<table><caption>**Table 3.** Relative proportion of class units of peat microorganisms over the peat column.</caption><thead><tr><th rowspan="3">Abbreviation</th><th rowspan="3">Taxon Unit (Phylum: Class)</th><th colspan="6">Peat Horizon Depth, cm</th></tr><tr><th colspan="6">Share of 16S rRNA Gene Sequences in % of the Number of Obtained Sequences</th></tr><tr><th>10–20</th><th>30–40</th><th>50–60</th><th>70–80</th><th>90–100</th><th>100–110</th></tr></thead><tbody><tr><td>EuMeth</td><td>Euryarchaeota; Methanobacteria</td><td>0.0</td><td>0.1</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td></tr><tr><td>ThNit</td><td>Thaumarchaeota; Nitrososphaeria</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.5</td><td>0.1</td></tr><tr><td>Ac</td><td>Acidobacteria; Acidobacteriia</td><td>20.5</td><td>46.0</td><td>64.7</td><td>44.8</td><td>1.0</td><td>4.9</td></tr><tr><td>Ac6</td><td>Acidobacteria; Subgroup 6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.4</td><td>0.1</td></tr><tr><td>AcAci</td><td>Actinobacteria; Acidimicrobiia</td><td>2.1</td><td>2.8</td><td>1.3</td><td>0.1</td><td>0.6</td><td>0.1</td></tr><tr><td>AcAct</td><td>Actinobacteria; Actinobacteria</td><td>10.6</td><td>1.8</td><td>12.1</td><td>0.4</td><td>18.4</td><td>1.0</td></tr><tr><td>AcOt</td><td>Actinobacteria; Others</td><td>0.6</td><td>20.7</td><td>4.6</td><td>0.3</td><td>0.2</td><td>0.2</td></tr><tr><td>AcTh</td><td>Actinobacteria; Thermoleophilia</td><td>2.1</td><td>0.4</td><td>1.8</td><td>1.1</td><td>1.1</td><td>0.3</td></tr><tr><td>Ar</td><td>Armatimonadetes; Armatimonadia</td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ArC</td><td>Armatimonadetes; Chthonomonadetes</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ArFi</td><td>Armatimonadetes; Fimbriimonadia</td><td>0.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.1</td><td>0.0</td></tr><tr><td>Bac</td><td>Bacteroidetes; Bacteroidia</td><td>2.6</td><td>0.0</td><td>0.2</td><td>0.0</td><td>6.1</td><td>0.4</td></tr><tr><td>Chl</td><td>Chloroflexi; Chloroflexia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.7</td><td>0.0</td></tr><tr><td>CyOt</td><td>Cyanobacteria; Oxyphotobacteria</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td></tr><tr><td>FirB</td><td>Firmicutes; Bacilli</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td><td>10.7</td><td>0.7</td></tr><tr><td>FirCl</td><td>Firmicutes; Clostridia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.7</td><td>4.9</td><td>1.8</td></tr><tr><td>Gem</td><td>Gemmatimonadetes; Gemmatimonadetes</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.6</td><td>0.0</td></tr><tr><td>GLon</td><td>Gemmatimonadetes; Longimicrobia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td></tr><tr><td>Ot</td><td>Others; Others</td><td>13.3</td><td>17.5</td><td>2.1</td><td>0.6</td><td>1.9</td><td>0.4</td></tr><tr><td>PatS</td><td>Patescibacteria; Saccharimonadia</td><td>0.0</td><td>0.6</td><td>1.3</td><td>0.0</td><td>0.1</td><td>0.0</td></tr><tr><td>PIPh</td><td>Planctomycetes; Phycisphaerae</td><td>3.9</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td></tr><tr><td>PIPl</td><td>Planctomycetes; Planctomycetacia</td><td>5.7</td><td>0.0</td><td>0.7</td><td>0.0</td><td>0.7</td><td>0.0</td></tr><tr><td>ProA</td><td>Proteobacteria; Alphaproteobacteria</td><td>9.4</td><td>6.3</td><td>9.7</td><td>50.2</td><td>24.6</td><td>1.7</td></tr><tr><td>ProD</td><td>Proteobacteria; Deltaproteobacteria</td><td>2.9</td><td>0.6</td><td>0.0</td><td>0.0</td><td>1.6</td><td>0.1</td></tr><tr><td>ProG</td><td>Proteobacteria; Gammaproteobacteria</td><td>12.7</td><td>1.6</td><td>1.1</td><td>1.8</td><td>17.7</td><td>87.5</td></tr><tr><td>ProOt</td><td>Proteobacteria; Others</td><td>1.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.4</td><td>0.0</td></tr><tr><td>Ver</td><td>Verrucomicrobia; Verrucomicrobiae</td><td>10.3</td><td>1.0</td><td>0.1</td><td>0.0</td><td>5.3</td><td>0.3</td></tr></tbody></table>

The descriptions of microbial variety at the order, family and genus level are provided in Figure 5 and Tables S3–S5 and briefly summarized below. However, these levels were not used for correlations with peat chemical parameters. The order *Acidobacteriales* predominated among other orders and constituted 29.6% of all varieties. The orders *Clostridiales* and *Rhizobiales* were the most diverse among other orders and included three families; however, they were present in minor numbers. The order *Pseudomonadales* included two families. The other orders were represented by one family. *Acidobacteriaceae* (Subgroup 1) predominated and increased in number with depth, up to 50–60 cm, where the share reached 64.7% and then started to decrease. The next in the number of families was *Burkholderiaceae* (12.6%),
class *Gammaproteobacteria* were concentrated at the depth of 90–110 cm (17.7–87.5%). The class *Actinobacteria* exhibited sizable and non-systematic variations over depth, reaching the largest proportion (18.4%) at the depth of 90 to 100 cm. Classes *Bacteroidia* (6.1%), *Bacilli* (10.7%) and *Clostridia* (4.9%) were mostly present also at the same depth. Classes *Chthonomonadetes*, *Phycisphaerae*, *Planctomycetacia* and other non-identified bacteria predominated in the active layers (between 10 and 40 cm depth). Classes *Deltaproteobacteria* and *Verrucomicrobiae* were found in the greatest number at the surface and bottom horizon of the peat column: 2.9% and 10.3%, respectively, in the 10–20 cm layer; 1.6% and 5.3%, respectively, in the 90–100 cm layer. Class *Thermoleophilia* was almost uniformly distributed over the entire peat profile (Figure 4b). The classes *Gammaproteobacteria*, *Actinobacteria* and *Alphaproteobacteria* were most diverse and included six, five, and four orders, respectively. Classes *Acidobacteriia*, *Bacilli*, *Planctomycetacia* and *Verrucomicrobia* included two orders each. Other classes were represented only by one order.

<table><caption>**Table 3.** Relative proportion of class units of peat microorganisms over the peat column.</caption><thead><tr><th rowspan="3">Abbreviation</th><th rowspan="3">Taxon Unit (Phylum: Class)</th><th colspan="6">Peat Horizon Depth, cm</th></tr><tr><th colspan="6">Share of 16S rRNA Gene Sequences in % of the Number of Obtained Sequences</th></tr><tr><th>10–20</th><th>30–40</th><th>50–60</th><th>70–80</th><th>90–100</th><th>100–110</th></tr></thead><tbody><tr><td>EuMeth</td><td>Euryarchaeota; Methanobacteria</td><td>0.0</td><td>0.1</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td></tr><tr><td>ThNit</td><td>Thaumarchaeota; Nitrososphaeria</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.5</td><td>0.1</td></tr><tr><td>Ac</td><td>Acidobacteria; Acidobacteriia</td><td>20.5</td><td>46.0</td><td>64.7</td><td>44.8</td><td>1.0</td><td>4.9</td></tr><tr><td>Ac6</td><td>Acidobacteria; Subgroup 6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.4</td><td>0.1</td></tr><tr><td>AcAci</td><td>Actinobacteria; Acidimicrobiia</td><td>2.1</td><td>2.8</td><td>1.3</td><td>0.1</td><td>0.6</td><td>0.1</td></tr><tr><td>AcAct</td><td>Actinobacteria; Actinobacteria</td><td>10.6</td><td>1.8</td><td>12.1</td><td>0.4</td><td>18.4</td><td>1.0</td></tr><tr><td>AcOt</td><td>Actinobacteria; Others</td><td>0.6</td><td>20.7</td><td>4.6</td><td>0.3</td><td>0.2</td><td>0.2</td></tr><tr><td>AcTh</td><td>Actinobacteria; Thermoleophilia</td><td>2.1</td><td>0.4</td><td>1.8</td><td>1.1</td><td>1.1</td><td>0.3</td></tr><tr><td>Ar</td><td>Armatimonadetes; Armatimonadia</td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ArC</td><td>Armatimonadetes; Chthonomonadetes</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ArFi</td><td>Armatimonadetes; Fimbriimonadia</td><td>0.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.1</td><td>0.0</td></tr><tr><td>Bac</td><td>Bacteroidetes; Bacteroidia</td><td>2.6</td><td>0.0</td><td>0.2</td><td>0.0</td><td>6.1</td><td>0.4</td></tr><tr><td>Chl</td><td>Chloroflexi; Chloroflexia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.7</td><td>0.0</td></tr><tr><td>CyOt</td><td>Cyanobacteria; Oxyphotobacteria</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td></tr><tr><td>FirB</td><td>Firmicutes; Bacilli</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td><td>10.7</td><td>0.7</td></tr><tr><td>FirCl</td><td>Firmicutes; Clostridia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.7</td><td>4.9</td><td>1.8</td></tr><tr><td>Gem</td><td>Gemmatimonadetes; Gemmatimonadetes</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.6</td><td>0.0</td></tr><tr><td>GLon</td><td>Gemmatimonadetes; Longimicrobia</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td></tr><tr><td>Ot</td><td>Others; Others</td><td>13.3</td><td>17.5</td><td>2.1</td><td>0.6</td><td>1.9</td><td>0.4</td></tr><tr><td>PatS</td><td>Patescibacteria; Saccharimonadia</td><td>0.0</td><td>0.6</td><td>1.3</td><td>0.0</td><td>0.1</td><td>0.0</td></tr><tr><td>PIPh</td><td>Planctomycetes; Phycisphaerae</td><td>3.9</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td></tr><tr><td>PIPl</td><td>Planctomycetes; Planctomycetacia</td><td>5.7</td><td>0.0</td><td>0.7</td><td>0.0</td><td>0.7</td><td>0.0</td></tr><tr><td>ProA</td><td>Proteobacteria; Alphaproteobacteria</td><td>9.4</td><td>6.3</td><td>9.7</td><td>50.2</td><td>24.6</td><td>1.7</td></tr><tr><td>ProD</td><td>Proteobacteria; Deltaproteobacteria</td><td>2.9</td><td>0.6</td><td>0.0</td><td>0.0</td><td>1.6</td><td>0.1</td></tr><tr><td>ProG</td><td>Proteobacteria; Gammaproteobacteria</td><td>12.7</td><td>1.6</td><td>1.1</td><td>1.8</td><td>17.7</td><td>87.5</td></tr><tr><td>ProOt</td><td>Proteobacteria; Others</td><td>1.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.4</td><td>0.0</td></tr><tr><td>Ver</td><td>Verrucomicrobia; Verrucomicrobiae</td><td>10.3</td><td>1.0</td><td>0.1</td><td>0.0</td><td>5.3</td><td>0.3</td></tr></tbody></table>

The descriptions of microbial variety at the order, family and genus level are provided in Figure 5 and Tables S3–S5 and briefly summarized below. However, these levels were not used for correlations with peat chemical parameters. The order *Acidobacteriales* predominated among other orders and constituted 29.6% of all varieties. The orders *Clostridiales* and *Rhizobiales* were the most diverse among other orders and included three families; however, they were present in minor numbers. The order *Pseudomonadales* included two families. The other orders were represented by one family. *Acidobacteriaceae* (Subgroup 1) predominated and increased in number with depth, up to 50–60 cm, where the share reached 64.7% and then started to decrease. The next in the number of families was *Burkholderiaceae* (12.6%),
which exhibited a maximum at the peat bottom and mineral horizons (100–110 cm), reaching 87.5% abundance, whereas at the depth of 30–80 it was represented in minor quantities. Families *Sphingomonadaceae*, *Bacillaceae* and *Propionibacteriaceae* were concentrated at the bottom of the frozen peat, in the 90–100 cm layer. *Acetobacteraceae* was evenly distributed over first 60 cm of depth. The families *Chitinophagoceae*, *Acidothermaceae* and *Xanthobacteraceae* were mostly concentrated within the active layer, at the 10–20 cm depth level. The family *Clostridiaceae* 1 was mainly present in the mineral layer, at a depth of 100–110 cm. Among all orders, *Clostridiales* and *Rhizobiales* were most diverse and included three families; however, they were present in minor numbers. The order *Pseudomonadales* included two families. The other orders were represented by one family. Microorganisms of the family *Acidobacteriaceae* (Subgroup 1) from the order *Acidobacteriales* constituted the major part of the phylum *Acidobacteria* (17%). The abundance of family members increased in number from the surface to the 50–60 cm layer, where their share reached 58%, and then started to decrease with depth. Among *Proteobacteria*, the maximum abundance was encountered for the family *Burkholderiaceae* (12.6%) from the order *Betaproteobacteriales*. The abundance of microorganisms of this family increased sharply at a depth of 90–110 cm (up to 83.8%), whereas at a depth of 30–80 cm it was present only in minor quantities. The bacteria of families *Sphingomonadaceae*, *Bacillaceae* and *Propionibacteriaceae* were concentrated at a depth of 90 to 100 cm. *Acetobacteraceae* was evenly distributed at a depth of 10–60 cm (first three layers). The remaining families were represented in minor quantities (less 1%) and mainly occurred in the upper layers of the peat core.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Order-level bacterial and archaeal community structure in the peat core based on 16S rRNA.</figcaption></figure>
which exhibited a maximum at the peat bottom and mineral horizons (100–110 cm), reaching 87.5% abundance, whereas at the depth of 30–80 it was represented in minor quantities. Families *Sphingomonadaceae*, *Bacillaceae* and *Propionibacteriaceae* were concentrated at the bottom of the frozen peat, in the 90–100 cm layer. *Acetobacteraceae* was evenly distributed over first 60 cm of depth. The families *Chitinophagoceae*, *Acidothermaceae* and *Xanthobacteraceae* were mostly concentrated within the active layer, at the 10–20 cm depth level. The family *Clostridiaceae* 1 was mainly present in the mineral layer, at a depth of 100–110 cm. Among all orders, *Clostridiales* and *Rhizobiales* were most diverse and included three families; however, they were present in minor numbers. The order *Pseudomonadales* included two families. The other orders were represented by one family. Microorganisms of the family *Acidobacteriaceae* (Subgroup 1) from the order *Acidobacteriales* constituted the major part of the phylum *Acidobacteria* (17%). The abundance of family members increased in number from the surface to the 50–60 cm layer, where their share reached 58%, and then started to decrease with depth. Among *Proteobacteria*, the maximum abundance was encountered for the family *Burkholderiaceae* (12.6%) from the order *Betaproteobacteriales*. The abundance of microorganisms of this family increased sharply at a depth of 90–110 cm (up to 83.8%), whereas at a depth of 30–80 cm it was present only in minor quantities. The bacteria of families *Sphingomonadaceae*, *Bacillaceae* and *Propionibacteriaceae* were concentrated at a depth of 90 to 100 cm. *Acetobacteraceae* was evenly distributed at a depth of 10–60 cm (first three layers). The remaining families were represented in minor quantities (less 1%) and mainly occurred in the upper layers of the peat core.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Order-level bacterial and archaeal community structure in the peat core based on 16S rRNA.</figcaption></figure>
Generally, the richness of life forms gradually decreased with depth. However, at the bottom of the frozen peat core (90–100 cm layer), it sharply increased. Most genera were represented in minor quantities. Dominant protists belonged to the genera *Granulicella* (13.3%) and *Paucibacter* (10.1%). The genus *Granulicella* was mainly distributed within the 30 to 80 cm depth. The genus *Paucibacter*, on the contrary, was concentrated in the mineral layer, at 100–110 cm depth (87.8%). More than 19% of protists belonged to unknown genera of phyla *Actinobacteria*, *Bacteroidetes*, *Planctomycetes*, *Proteobacteria*, and *Verrucomicrobia*.

### 3.4. Pair Correlation and Multiparametric Statistics of the Full Data Set

First, we checked for linear correlations between specific classes of microorganisms and chemical parameters of peat and peat fluids (water and ice). It can be seen from Table S6 that numerous groups of microorganisms (EuMeth (*Methanobacteria*), ThNit (*Nitrosphaeria*), Chl (*Chloroflexia*), CytOt (*Oxyphotobacteria*), FirB (*Bacilli*), FirCl (*Clostridia*), Gem (*Gemmatimonadetes*)) strongly correlated with macro- (K, Fe, Mg) and micro- (Rb, Mn, Mo, Zn, Ni, Co...) nutrients as well as indifferent geochemical tracers (trivalent and tetravalent hydrolysates) in the peat water and ice. Only a few classes exhibited strong correlations with solid peat chemical composition, such as Fe, Mg, δ¹⁵N (AcAci), Ca, Sr (AcOt), macro- and micronutrients, and δ¹³C (Ar (*Armatimonadia*), ArC (*Chthonomonadetes*), PIPh (*Phycisphaerae*), PPII (*Planctomycetacia*), ProOt (*Proteobacteria*; other), and Ver (*Verrucomicrobiae*).

In order to achieve further insights into the multiple correlations between physical, chemical and biological parameters, we performed multiparametric statistical treatment of the full dataset, considering all available temperature data as well as peat physical and chemical properties. For this treatment, we considered quantitative data on microorganism numbers and the percentage of different classes in each of the five organic (peat) layers and the one mineral layer. The temperature in various peat layers, the inorganic chemical composition of peat, pore water and ice, the organic (C, N, δ¹³C, δ¹⁵N) composition, and the organic matter quality in each measured layer were also considered. Note that for this treatment, the depth of each layer was considered as an independent variable.

Two PCA treatments of the peat core were performed: with and without the mineral layer. However, due to insufficient data on the mineral layer (lack of ice chemistry and OM quality), the explicatory capacity of PCA for both the organic and mineral part of the core was much lower than that for the peat part solely. As such, we present only the PCA results for the organic part of the peat column.

There were three main factors potentially responsible for the observed variability in the physical, chemical and microbial parameters of the peat core (Table 4, Figure 6). The first factor comprised bottom layers of the column with minimal temperature variations, high peat degradation degree, high contribution of alkyl C, aromatic C and carboxylic C and elevated concentrations of B, Cr and Cu in the peat and DOC, Ca, Mg, K, Si, Fe and various micronutrients (Mn, Ba, Zn, Co, Rb, Mo, B) in the peat water and ice. This factor provided the highest fraction of 16S r-RNA of *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Nitrososphaeria*, *Acidobacteria* (Subgroup 6) and *Chloroflexia*. These classes of microorganisms were primarily controlled by peat water and ice chemistry. The second factor acted on the chemical composition of solid peat and marked the minimal proportion of genes of *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteriia*. These classes avoided the medium part of the peat core, where there was an elevated concentration of lithogenic elements (Ti, Y, REE) and a minimal concentration of some nutrients (K, Rb, V, Mn) and toxicants (As, Pb, Sb). Presumably, these elements marked the presence of silicate dust, and thus the classes of microorganisms of the second factor preferred the peat layers where the proportion of mineral dust was minimal. Finally, the third factor acted on *Actinobacteria*, which were linked to low N and Mo concentration and the highest δ¹⁵N and C:N ratio in the peat and low P concentration in water and ice.
Generally, the richness of life forms gradually decreased with depth. However, at the bottom of the frozen peat core (90–100 cm layer), it sharply increased. Most genera were represented in minor quantities. Dominant protists belonged to the genera *Granulicella* (13.3%) and *Paucibacter* (10.1%). The genus *Granulicella* was mainly distributed within the 30 to 80 cm depth. The genus *Paucibacter*, on the contrary, was concentrated in the mineral layer, at 100–110 cm depth (87.8%). More than 19% of protists belonged to unknown genera of phyla *Actinobacteria*, *Bacteroidetes*, *Planctomycetes*, *Proteobacteria*, and *Verrucomicrobia*.

### 3.4. Pair Correlation and Multiparametric Statistics of the Full Data Set

First, we checked for linear correlations between specific classes of microorganisms and chemical parameters of peat and peat fluids (water and ice). It can be seen from Table S6 that numerous groups of microorganisms (EuMeth (*Methanobacteria*), ThNit (*Nitrosphaeria*), Chl (*Chloroflexia*), CytOt (*Oxyphotobacteria*), FirB (*Bacilli*), FirCl (*Clostridia*), Gem (*Gemmatimonadetes*)) strongly correlated with macro- (K, Fe, Mg) and micro- (Rb, Mn, Mo, Zn, Ni, Co...) nutrients as well as indifferent geochemical tracers (trivalent and tetravalent hydrolysates) in the peat water and ice. Only a few classes exhibited strong correlations with solid peat chemical composition, such as Fe, Mg, δ¹⁵N (AcAci), Ca, Sr (AcOt), macro- and micronutrients, and δ¹³C (Ar (*Armatimonadia*), ArC (*Chthonomonadetes*), PIPh (*Phycisphaerae*), PPII (*Planctomycetacia*), ProOt (*Proteobacteria*; other), and Ver (*Verrucomicrobiae*).

In order to achieve further insights into the multiple correlations between physical, chemical and biological parameters, we performed multiparametric statistical treatment of the full dataset, considering all available temperature data as well as peat physical and chemical properties. For this treatment, we considered quantitative data on microorganism numbers and the percentage of different classes in each of the five organic (peat) layers and the one mineral layer. The temperature in various peat layers, the inorganic chemical composition of peat, pore water and ice, the organic (C, N, δ¹³C, δ¹⁵N) composition, and the organic matter quality in each measured layer were also considered. Note that for this treatment, the depth of each layer was considered as an independent variable.

Two PCA treatments of the peat core were performed: with and without the mineral layer. However, due to insufficient data on the mineral layer (lack of ice chemistry and OM quality), the explicatory capacity of PCA for both the organic and mineral part of the core was much lower than that for the peat part solely. As such, we present only the PCA results for the organic part of the peat column.

There were three main factors potentially responsible for the observed variability in the physical, chemical and microbial parameters of the peat core (Table 4, Figure 6). The first factor comprised bottom layers of the column with minimal temperature variations, high peat degradation degree, high contribution of alkyl C, aromatic C and carboxylic C and elevated concentrations of B, Cr and Cu in the peat and DOC, Ca, Mg, K, Si, Fe and various micronutrients (Mn, Ba, Zn, Co, Rb, Mo, B) in the peat water and ice. This factor provided the highest fraction of 16S r-RNA of *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Nitrososphaeria*, *Acidobacteria* (Subgroup 6) and *Chloroflexia*. These classes of microorganisms were primarily controlled by peat water and ice chemistry. The second factor acted on the chemical composition of solid peat and marked the minimal proportion of genes of *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteriia*. These classes avoided the medium part of the peat core, where there was an elevated concentration of lithogenic elements (Ti, Y, REE) and a minimal concentration of some nutrients (K, Rb, V, Mn) and toxicants (As, Pb, Sb). Presumably, these elements marked the presence of silicate dust, and thus the classes of microorganisms of the second factor preferred the peat layers where the proportion of mineral dust was minimal. Finally, the third factor acted on *Actinobacteria*, which were linked to low N and Mo concentration and the highest δ¹⁵N and C:N ratio in the peat and low P concentration in water and ice.
<table><caption>**Table 4.** Correlation of factors with peat parameters and share of bacterial taxa (class level).</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>Coordinate a* according to CIE lab (red color)</td><td>a*</td><td>-0.96</td><td>Verrucomicrobia; Verrucomicrobiae</td><td>* Ver</td><td>-0.99</td><td>N peat</td><td>N</td><td>-0.86</td></tr><tr><td>Depth peat (negative values)</td><td>Depth</td><td>-0.94</td><td>Proteobacteria; Deltaproteobacteria</td><td>* ProD</td><td>-0.98</td><td>Mo peat</td><td>MoP</td><td>-0.86</td></tr><tr><td>O-alkyl C</td><td>COalk</td><td>-0.93</td><td>Ti peat</td><td>TiP</td><td>-0.95</td><td>P ice water</td><td>PW</td><td>-0.74</td></tr><tr><td>SUVA peat ice/water</td><td>IW SUVA</td><td>-0.88</td><td>Proteobacteria; Others</td><td>* ProOt</td><td>-0.95</td><td>Sr peat</td><td>SrP</td><td>0.74</td></tr><tr><td>Coordinate b* according to CIE lab (yellow color)</td><td>b*</td><td>-0.86</td><td>Armatimonadetes; Fimbriimonadia</td><td>* ArFi</td><td>-0.89</td><td>Actinobacteria; Others</td><td>* AcOt</td><td>0.81</td></tr><tr><td>L* CIE lab</td><td>L*</td><td>-0.85</td><td>K peat</td><td>KP</td><td>-0.88</td><td>δ<sup>15</sup>N</td><td>15N</td><td>0.84</td></tr><tr><td>Mg peat</td><td>MgP</td><td>-0.82</td><td>La peat</td><td>LaP</td><td>-0.86</td><td>C/N</td><td>C/N</td><td>0.93</td></tr><tr><td>S peat</td><td>S</td><td>-0.82</td><td>Proteobacteria; Gammaproteobacteria</td><td>* ProG</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>P peat</td><td>PP</td><td>-0.77</td><td>Rb peat</td><td>RbP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Cd peat</td><td>CdP</td><td>-0.75</td><td>V peat</td><td>VP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Plab peat</td><td>Plab</td><td>-0.71</td><td>Cs peat</td><td>CsP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Maximum temperature in the horizon</td><td>MaxAn</td><td>-0.70</td><td>Planctomycetes; Phycisphaerae</td><td>* PlPh</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Zr peat</td><td>ZrP</td><td>0.71</td><td>Mean temperature in the horizon</td><td>MeAn</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Mo ice water</td><td>MoW</td><td>0.72</td><td>Ce peat</td><td>CeP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Nb peat</td><td>NbP</td><td>0.76</td><td>Mn peat</td><td>MnP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Be ice water</td><td>BeW</td><td>0.76</td><td>Na peat</td><td>NaP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>U peat</td><td>UP</td><td>0.76</td><td>Dy peat</td><td>DyP</td><td>-0.84</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 4.** Correlation of factors with peat parameters and share of bacterial taxa (class level).</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>Coordinate a* according to CIE lab (red color)</td><td>a*</td><td>-0.96</td><td>Verrucomicrobia; Verrucomicrobiae</td><td>* Ver</td><td>-0.99</td><td>N peat</td><td>N</td><td>-0.86</td></tr><tr><td>Depth peat (negative values)</td><td>Depth</td><td>-0.94</td><td>Proteobacteria; Deltaproteobacteria</td><td>* ProD</td><td>-0.98</td><td>Mo peat</td><td>MoP</td><td>-0.86</td></tr><tr><td>O-alkyl C</td><td>COalk</td><td>-0.93</td><td>Ti peat</td><td>TiP</td><td>-0.95</td><td>P ice water</td><td>PW</td><td>-0.74</td></tr><tr><td>SUVA peat ice/water</td><td>IW SUVA</td><td>-0.88</td><td>Proteobacteria; Others</td><td>* ProOt</td><td>-0.95</td><td>Sr peat</td><td>SrP</td><td>0.74</td></tr><tr><td>Coordinate b* according to CIE lab (yellow color)</td><td>b*</td><td>-0.86</td><td>Armatimonadetes; Fimbriimonadia</td><td>* ArFi</td><td>-0.89</td><td>Actinobacteria; Others</td><td>* AcOt</td><td>0.81</td></tr><tr><td>L* CIE lab</td><td>L*</td><td>-0.85</td><td>K peat</td><td>KP</td><td>-0.88</td><td>δ<sup>15</sup>N</td><td>15N</td><td>0.84</td></tr><tr><td>Mg peat</td><td>MgP</td><td>-0.82</td><td>La peat</td><td>LaP</td><td>-0.86</td><td>C/N</td><td>C/N</td><td>0.93</td></tr><tr><td>S peat</td><td>S</td><td>-0.82</td><td>Proteobacteria; Gammaproteobacteria</td><td>* ProG</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>P peat</td><td>PP</td><td>-0.77</td><td>Rb peat</td><td>RbP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Cd peat</td><td>CdP</td><td>-0.75</td><td>V peat</td><td>VP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Plab peat</td><td>Plab</td><td>-0.71</td><td>Cs peat</td><td>CsP</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Maximum temperature in the horizon</td><td>MaxAn</td><td>-0.70</td><td>Planctomycetes; Phycisphaerae</td><td>* PlPh</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Zr peat</td><td>ZrP</td><td>0.71</td><td>Mean temperature in the horizon</td><td>MeAn</td><td>-0.85</td><td></td><td></td><td></td></tr><tr><td>Mo ice water</td><td>MoW</td><td>0.72</td><td>Ce peat</td><td>CeP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Nb peat</td><td>NbP</td><td>0.76</td><td>Mn peat</td><td>MnP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Be ice water</td><td>BeW</td><td>0.76</td><td>Na peat</td><td>NaP</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>U peat</td><td>UP</td><td>0.76</td><td>Dy peat</td><td>DyP</td><td>-0.84</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>aromatic C</td><td>Carom</td><td>0.79</td><td>δ<sup>13</sup>C peat</td><td>13C</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>B ice water</td><td>BW</td><td>0.80</td><td>pH ice water</td><td>IW pH</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Euryarchaeota; Methanobacteria</td><td>* EuMeth</td><td>0.81</td><td>Y peat</td><td>YP</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>As ice water</td><td>AsW</td><td>0.81</td><td>Planctomycetes; Planctomycetacia</td><td>* PlPh</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>carboxylic C</td><td>Ccarb</td><td>0.82</td><td>Pb peat</td><td>PbP</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>Dy ice water</td><td>DyW</td><td>0.82</td><td>As peat</td><td>AsP</td><td>-0.81</td><td></td><td></td><td></td></tr><tr><td>Ce ice water</td><td>CeW</td><td>0.83</td><td>Armatimonadetes; Chthonomonadetes</td><td>* ArC</td><td>-0.80</td><td></td><td></td><td></td></tr><tr><td colspan="9"><b>Table 3, continued.</b></td></tr><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr><tr><td>Y ice water</td><td>YW</td><td>0.83</td><td>Armatimonadetes; Armatimonadia</td><td>* Ar</td><td>-0.80</td><td></td><td></td><td></td></tr><tr><td>Minimum temperature in the horizon</td><td>MinAn</td><td>0.84</td><td>Sb peat</td><td>SbP</td><td>-0.79</td><td></td><td></td><td></td></tr><tr><td>Sb ice water</td><td>SbW</td><td>0.84</td><td>Yb peat</td><td>YbP</td><td>-0.77</td><td></td><td></td><td></td></tr><tr><td>Ba peat</td><td>BaP</td><td>0.85</td><td>Li peat</td><td>LiP</td><td>-0.74</td><td></td><td></td><td></td></tr><tr><td>Cu ice water</td><td>CuW</td><td>0.85</td><td>Bacteroidetes; Bacteroidia</td><td>* Bac</td><td>-0.71</td><td></td><td></td><td></td></tr><tr><td>La ice water</td><td>LaW</td><td>0.85</td><td>Acidobacteria; Acidobacteriia</td><td>* Ac</td><td>0.76</td><td></td><td></td><td></td></tr><tr><td>Cu peat</td><td>CuP</td><td>0.86</td><td>Bulk density peat</td><td>BD</td><td>0.79</td><td></td><td></td><td></td></tr><tr><td>U ice water</td><td>UW</td><td>0.88</td><td>DIC ice water</td><td>IW DIC</td><td>0.83</td><td></td><td></td><td></td></tr><tr><td>Decomposition of peat</td><td>DP</td><td>0.89</td><td>UV<sub>254</sub></td><td>UV<sub>254</sub></td><td>0.84</td><td></td><td></td><td></td></tr><tr><td>Alkyl C</td><td>Calkyl</td><td>0.90</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>aromatic C</td><td>Carom</td><td>0.79</td><td>δ<sup>13</sup>C peat</td><td>13C</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>B ice water</td><td>BW</td><td>0.80</td><td>pH ice water</td><td>IW pH</td><td>-0.84</td><td></td><td></td><td></td></tr><tr><td>Euryarchaeota; Methanobacteria</td><td>* EuMeth</td><td>0.81</td><td>Y peat</td><td>YP</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>As ice water</td><td>AsW</td><td>0.81</td><td>Planctomycetes; Planctomycetacia</td><td>* PlPh</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>carboxylic C</td><td>Ccarb</td><td>0.82</td><td>Pb peat</td><td>PbP</td><td>-0.82</td><td></td><td></td><td></td></tr><tr><td>Dy ice water</td><td>DyW</td><td>0.82</td><td>As peat</td><td>AsP</td><td>-0.81</td><td></td><td></td><td></td></tr><tr><td>Ce ice water</td><td>CeW</td><td>0.83</td><td>Armatimonadetes; Chthonomonadetes</td><td>* ArC</td><td>-0.80</td><td></td><td></td><td></td></tr><tr><td colspan="9"><b>Table 3, continued.</b></td></tr><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr><tr><td>Y ice water</td><td>YW</td><td>0.83</td><td>Armatimonadetes; Armatimonadia</td><td>* Ar</td><td>-0.80</td><td></td><td></td><td></td></tr><tr><td>Minimum temperature in the horizon</td><td>MinAn</td><td>0.84</td><td>Sb peat</td><td>SbP</td><td>-0.79</td><td></td><td></td><td></td></tr><tr><td>Sb ice water</td><td>SbW</td><td>0.84</td><td>Yb peat</td><td>YbP</td><td>-0.77</td><td></td><td></td><td></td></tr><tr><td>Ba peat</td><td>BaP</td><td>0.85</td><td>Li peat</td><td>LiP</td><td>-0.74</td><td></td><td></td><td></td></tr><tr><td>Cu ice water</td><td>CuW</td><td>0.85</td><td>Bacteroidetes; Bacteroidia</td><td>* Bac</td><td>-0.71</td><td></td><td></td><td></td></tr><tr><td>La ice water</td><td>LaW</td><td>0.85</td><td>Acidobacteria; Acidobacteriia</td><td>* Ac</td><td>0.76</td><td></td><td></td><td></td></tr><tr><td>Cu peat</td><td>CuP</td><td>0.86</td><td>Bulk density peat</td><td>BD</td><td>0.79</td><td></td><td></td><td></td></tr><tr><td>U ice water</td><td>UW</td><td>0.88</td><td>DIC ice water</td><td>IW DIC</td><td>0.83</td><td></td><td></td><td></td></tr><tr><td>Decomposition of peat</td><td>DP</td><td>0.89</td><td>UV<sub>254</sub></td><td>UV<sub>254</sub></td><td>0.84</td><td></td><td></td><td></td></tr><tr><td>Alkyl C</td><td>Calkyl</td><td>0.90</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>Nb ice water</td><td>NbW</td><td>0.90</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cyanobacteria; Oxyphotobacteria</td><td>* CyOt</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gemmatimonadetes; Gemmatimonadetes</td><td>* Gem</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Thaumarchaeota; Nitrososphaeria</td><td>* ThNit</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Acidobacteria; Subgroup 6</td><td>* Ac6</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chloroflexi; Chloroflexia</td><td>* Chl</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Zr ice peat</td><td>ZrW</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Yb ice peat</td><td>YbW</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Firmicutes; Bacilli</td><td>* FirB</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cr ice peat</td><td>CrW</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>B peat</td><td>BP</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cr peat</td><td>CrP</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ti ice peat</td><td>TiW</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Specific conductivity of peat ice waters</td><td>IW EC</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mn ice peat</td><td>MnW</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ba ice peat</td><td>BaW</td><td>0.94</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Firmicutes; Clostridia</td><td>* FirCl</td><td>0.94</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.1 (52.0%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.2 (27.3%)</th><th>Parameter</th><th>Abbreviated in Figure 6</th><th>F.3 (13.2%)</th></tr></thead><tbody><tr><td>Nb ice water</td><td>NbW</td><td>0.90</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cyanobacteria; Oxyphotobacteria</td><td>* CyOt</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gemmatimonadetes; Gemmatimonadetes</td><td>* Gem</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Thaumarchaeota; Nitrososphaeria</td><td>* ThNit</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Acidobacteria; Subgroup 6</td><td>* Ac6</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chloroflexi; Chloroflexia</td><td>* Chl</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Zr ice peat</td><td>ZrW</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Yb ice peat</td><td>YbW</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Firmicutes; Bacilli</td><td>* FirB</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cr ice peat</td><td>CrW</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>B peat</td><td>BP</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cr peat</td><td>CrP</td><td>0.92</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ti ice peat</td><td>TiW</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Specific conductivity of peat ice waters</td><td>IW EC</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mn ice peat</td><td>MnW</td><td>0.93</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ba ice peat</td><td>BaW</td><td>0.94</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Firmicutes; Clostridia</td><td>* FirCl</td><td>0.94</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><td>Parameter</td><td>Abbreviated in Figure 6</td><td>F.1 (52.0%)</td></tr></thead><tbody><tr><td>V ice peat</td><td>VW</td><td>0.95</td></tr><tr><td>Cs ice peat</td><td>CsW</td><td>0.95</td></tr><tr><td>Na ice peat</td><td>NaW</td><td>0.96</td></tr><tr><td>Ca ice peat</td><td>CaW</td><td>0.96</td></tr><tr><td>Li ice peat</td><td>LiW</td><td>0.96</td></tr><tr><td>Fe ice peat</td><td>FeW</td><td>0.96</td></tr><tr><td>Mg ice peat</td><td>MgW</td><td>0.97</td></tr><tr><td>Cd ice peat</td><td>CdW</td><td>0.97</td></tr><tr><td>K ice peat</td><td>KW</td><td>0.97</td></tr><tr><td>Rb ice peat</td><td>RbW</td><td>0.98</td></tr><tr><td>Zn ice peat</td><td>ZnW</td><td>0.98</td></tr><tr><td>Co ice peat</td><td>CoW</td><td>0.99</td></tr><tr><td>Ni ice peat</td><td>NiW</td><td>0.99</td></tr><tr><td>Sr ice peat</td><td>SrW</td><td>0.99</td></tr><tr><td>DOC ice peat</td><td>IW DOC</td><td>0.99</td></tr><tr><td>Si ice peat</td><td>SiW</td><td>1.00</td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><td>Parameter</td><td>Abbreviated in Figure 6</td><td>F.1 (52.0%)</td></tr></thead><tbody><tr><td>V ice peat</td><td>VW</td><td>0.95</td></tr><tr><td>Cs ice peat</td><td>CsW</td><td>0.95</td></tr><tr><td>Na ice peat</td><td>NaW</td><td>0.96</td></tr><tr><td>Ca ice peat</td><td>CaW</td><td>0.96</td></tr><tr><td>Li ice peat</td><td>LiW</td><td>0.96</td></tr><tr><td>Fe ice peat</td><td>FeW</td><td>0.96</td></tr><tr><td>Mg ice peat</td><td>MgW</td><td>0.97</td></tr><tr><td>Cd ice peat</td><td>CdW</td><td>0.97</td></tr><tr><td>K ice peat</td><td>KW</td><td>0.97</td></tr><tr><td>Rb ice peat</td><td>RbW</td><td>0.98</td></tr><tr><td>Zn ice peat</td><td>ZnW</td><td>0.98</td></tr><tr><td>Co ice peat</td><td>CoW</td><td>0.99</td></tr><tr><td>Ni ice peat</td><td>NiW</td><td>0.99</td></tr><tr><td>Sr ice peat</td><td>SrW</td><td>0.99</td></tr><tr><td>DOC ice peat</td><td>IW DOC</td><td>0.99</td></tr><tr><td>Si ice peat</td><td>SiW</td><td>1.00</td></tr></tbody></table>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Results of PCA on peat (and peat water/ice), without mineral layer. A plot of F1 × F2 and F2 × F3 factors (upper and low panel, respectively). The abbreviations used in the Figure are detailed in Table 3.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Results of PCA on peat (and peat water/ice), without mineral layer. A plot of F1 × F2 and F2 × F3 factors (upper and low panel, respectively). The abbreviations used in the Figure are detailed in Table 3.</figcaption></figure>
# 4. Discussion

## 4.1. Elementary Peat and Peat Fluid Composition and Organic Matter Quality

As also noted in the previous study of the peat core from a neighboring mound [20], the inter-layer variations of most major and trace elements in the frozen and thawed peat of the peat bog mound in this study are within the average elementary peat composition of western Siberia across the more than 1500 km latitudinal transect [55]. In particular, an enrichment of the upper 0–20 cm horizon in macro- (P, K, Mg) and micro-nutrients (Mn, Zn, Rb) as well as Cs and potential toxicants of atmospheric origin (As, Cd, Sb, Tl, Pb) reflects the role of ground vegetation in providing these elements via vertical downward water migration. As such, the sampled peat core may be considered representative of the permafrost-affected part of western Siberia. For most elements, the bottom (frozen) part of the peat core exhibited a composition that was quite similar to the upper (thawed) part of the core. Using a Mann-Whitney test, significant ($p < 0.05$) differences between frozen and thawed layers were observed only for C:N ratio, Na, Mg, K, Al, Ca, Zn, Ge, Hf, W, Nb, Pb, U, Ba, Cd, Sb and Th. The differences ($p < 0.05$) in element concentration between peat ice and peat pore water were observed for S.C., SUVA, Li, Al, Si, P, Fe, Co, Rb, Cs, DOC, Ca, Sr, Mg, Ni, Cd and Tl.

Concerning peat pore water and ice chemical composition, a local maximum of DOC, UV absorbency, specific conductivity and some trace element concentrations at 40–50 cm depth (Figure 3) suggested some solute and organic colloid concentrations via freezing front migration, as demonstrated in the peat ice study across the WSL [58]. However, this pattern of concentration was different for other nutrients, such as phosphorus, which exhibited a peak at 80 cm depth. Finally, the downward increase in the concentration of micro-nutrients (B, K, Mn, Ni, Zn, Mo, Ba) suggested their possible use for biota in case of a massive ice peat thaw. In contrast, enrichment of peat pore waters from the surface layer (0–10 cm) in toxicants (As, Sb, Pb) suggested deposition of these elements from the atmosphere.

The increase of alkyl C:O alkyl C ratio and of aromatic C proportion with depth along the peat profile may reflect an advanced degree of OM degradation [61]. From 0 to 35 cm depth, lichen and sphagnum are the main contributors, whereas, from 65 to 95 cm depth, wood pieces and shrub dominate the plant remains in the peat [20]. As a result, an increase in aromatic C (in relation to phenols from lignin) and a decrease in O-alkyl C proportions may also result from the variable contribution of different plant species to peat formation [62,63]. Except for the 19–21 cm depth, which presents slightly higher value, the $\delta^{15}$N along the peat core is in the range given by Amundson et al. [64] for the soils of this climatic region, and the absence of isotopic enrichment with depth is also consistent with data from water-saturated soils [65].

## 4.2. Bacterial Number and Genetic Diversity of Soil Microorganisms

Numerous data available for the bacterial counts in Siberian soils demonstrated that the total number of cells generally decreases with the depth ([66] and references therein). In accordance with recent studies of enzyme patterns and microbial community compositions in western Siberian soils [34], as well as metabolic activity with respect to various organic substrates assessed by the Biolog Ecoplate technique [20], we demonstrated that each peat soil horizon harbored different microbial consortia.

The total number of cells of microorganisms quantified in this study (from $1.35 \times 10^{10}$ at the surface to $2.27 \times 10^7$ cell $g_{\text{peat}}^{-1}$ in the mineral layer) is comparable to the total viable bacterial count of a similar peat core sampled three years before, as measured by the DAPI technique ($73 - 1.6) \times 10^7$ cell $g_{\text{peat}}^{-1}$ [20]. The latter values are within the range measured in the Antarctic and High Canadian Arctic permafrost ($10^5-10^7$ cells/g$_{\text{soil}}$ [48,67] and ($2-4) \times 10^7$ cell/g$_{\text{soil}}$ [1]) and Siberian permafrost ($10^7-10^8$ cells/g$_{\text{soil}}$ [68,69]). Total cell count ranged from $9.32 \times 10^6$ cells per g dry weight to $1.37 \times 10^7$ cell/g$_{\text{soil}}$ in the permafrost samples from Alaska [70]. Note, however, the quite high number of microbial
# 4. Discussion

## 4.1. Elementary Peat and Peat Fluid Composition and Organic Matter Quality

As also noted in the previous study of the peat core from a neighboring mound [20], the inter-layer variations of most major and trace elements in the frozen and thawed peat of the peat bog mound in this study are within the average elementary peat composition of western Siberia across the more than 1500 km latitudinal transect [55]. In particular, an enrichment of the upper 0–20 cm horizon in macro- (P, K, Mg) and micro-nutrients (Mn, Zn, Rb) as well as Cs and potential toxicants of atmospheric origin (As, Cd, Sb, Tl, Pb) reflects the role of ground vegetation in providing these elements via vertical downward water migration. As such, the sampled peat core may be considered representative of the permafrost-affected part of western Siberia. For most elements, the bottom (frozen) part of the peat core exhibited a composition that was quite similar to the upper (thawed) part of the core. Using a Mann-Whitney test, significant ($p < 0.05$) differences between frozen and thawed layers were observed only for C:N ratio, Na, Mg, K, Al, Ca, Zn, Ge, Hf, W, Nb, Pb, U, Ba, Cd, Sb and Th. The differences ($p < 0.05$) in element concentration between peat ice and peat pore water were observed for S.C., SUVA, Li, Al, Si, P, Fe, Co, Rb, Cs, DOC, Ca, Sr, Mg, Ni, Cd and Tl.

Concerning peat pore water and ice chemical composition, a local maximum of DOC, UV absorbency, specific conductivity and some trace element concentrations at 40–50 cm depth (Figure 3) suggested some solute and organic colloid concentrations via freezing front migration, as demonstrated in the peat ice study across the WSL [58]. However, this pattern of concentration was different for other nutrients, such as phosphorus, which exhibited a peak at 80 cm depth. Finally, the downward increase in the concentration of micro-nutrients (B, K, Mn, Ni, Zn, Mo, Ba) suggested their possible use for biota in case of a massive ice peat thaw. In contrast, enrichment of peat pore waters from the surface layer (0–10 cm) in toxicants (As, Sb, Pb) suggested deposition of these elements from the atmosphere.

The increase of alkyl C:O alkyl C ratio and of aromatic C proportion with depth along the peat profile may reflect an advanced degree of OM degradation [61]. From 0 to 35 cm depth, lichen and sphagnum are the main contributors, whereas, from 65 to 95 cm depth, wood pieces and shrub dominate the plant remains in the peat [20]. As a result, an increase in aromatic C (in relation to phenols from lignin) and a decrease in O-alkyl C proportions may also result from the variable contribution of different plant species to peat formation [62,63]. Except for the 19–21 cm depth, which presents slightly higher value, the $\delta^{15}$N along the peat core is in the range given by Amundson et al. [64] for the soils of this climatic region, and the absence of isotopic enrichment with depth is also consistent with data from water-saturated soils [65].

## 4.2. Bacterial Number and Genetic Diversity of Soil Microorganisms

Numerous data available for the bacterial counts in Siberian soils demonstrated that the total number of cells generally decreases with the depth ([66] and references therein). In accordance with recent studies of enzyme patterns and microbial community compositions in western Siberian soils [34], as well as metabolic activity with respect to various organic substrates assessed by the Biolog Ecoplate technique [20], we demonstrated that each peat soil horizon harbored different microbial consortia.

The total number of cells of microorganisms quantified in this study (from $1.35 \times 10^{10}$ at the surface to $2.27 \times 10^7$ cell $g_{\text{peat}}^{-1}$ in the mineral layer) is comparable to the total viable bacterial count of a similar peat core sampled three years before, as measured by the DAPI technique ($73 - 1.6) \times 10^7$ cell $g_{\text{peat}}^{-1}$ [20]. The latter values are within the range measured in the Antarctic and High Canadian Arctic permafrost ($10^5-10^7$ cells/g$_{\text{soil}}$ [48,67] and ($2-4) \times 10^7$ cell/g$_{\text{soil}}$ [1]) and Siberian permafrost ($10^7-10^8$ cells/g$_{\text{soil}}$ [68,69]). Total cell count ranged from $9.32 \times 10^6$ cells per g dry weight to $1.37 \times 10^7$ cell/g$_{\text{soil}}$ in the permafrost samples from Alaska [70]. Note, however, the quite high number of microbial
cells in the surface layer of peat core, which exceeds the typical values of mineral soils by a factor of 100 to 1000.

In agreement with Jansson and Taş [69], the diversity of bacteria in the permafrost peat assessed in this study is much higher than that of archaea (99% of all sequences in this study). The dominant phylum *Proteobacteria*, *Firmicutes*, *Acidobacteria* and *Actinobacteria* identified in the present study are also reported in other permafrost environments [36,69]. The microbial population of the peat soil of Central Siberia was dominated by bacteria (96.6–98.7% of all sequences) and included *Proteobacteria* (34–25%), *Actinobacteria* (20–18%), *Acidobacteria* (8–14%) and *Verrucomicrobia* (13–14%), as reported by Grodnitskaya et al. [33]. Similarly, polygonal peat bogs of the Canadian Arctic contained *Proteobacteria* (40–50%), *Bacteroidetes* (20–40%) and *Acidobacteria* (10–15%) [71].

In studied peat core, *Proteobacteria* were abundant in all layers and accounted for 47% of all sequences, on average. This was represented by *Alphaproteobacteria* (5–50%) and *Gammaproteobacteria*, with maximal abundance at the mineral layer (84%). The possibility of *Proteobacteria* synthesizing DNA at a temperature from 0 to -20 °C [72] may explain their abundance at the 60–110 cm depth, where the temperature ranged from -0.1 to -0.5 °C. It is noteworthy that *Alphaproteobacteria* are capable of growing after nitrogen addition [35], which can explain the dominance of this class at 70–80 cm depth, where the N concentration achieves its maximum. The main representative of *Alphaproteobacteria* is clade *Sphingomonades*, which is abundant at the depth of 90 to 100 cm (20.3%). This clade is known to degrade lignocellulose debris [73], which is consistent with the botanical composition of studied peat: the wood debris is highly abundant below 80 cm depth. The representatives of clade *Rhizobiales* are highly abundant at 0–20 cm depth. This order includes gram-negative soil bacteria, capable of fixing nitrogen in symbiosis with plants [74–79]. Overall, the near-surface layer of peat is the major location for root exudates, most likely in the form of labile C substances, which are known to favor the growth of copiotrophic microorganisms [23] such as *Gamma*- and *Alphaproteobacteria*, as also evidenced from a study of catabolic activity of peat microorganisms [20].

Class *Acidobacteria* is highly abundant in soils and may constitute up to 52% of the microbial community [80]. In studied peat core, oligotrophic *Acidobacteria* were highly abundant and mostly present in the thawed layer or permafrost boundary (50–60 cm depth) and strongly decreased in the deep-frozen and mineral horizon. This may reflect the need of these bacteria for nutrients, which are concentrating at this boundary. Moreover, the peat water from 50–60 cm depth exhibited the lowest pH of all sampled horizons, consistent with the preference of *Acidobacteria* for acidic environments.

The *Actinobacteria* were present over the full depth of the peat core, which might be due to their ability to assimilate a broad spectrum of organic substrates and their adaptation to low temperature, for example, due to elevated content of guanidine and cytosine [80,81]. The *Bacteroidetes* are the main representatives of gram-negative anaerobic bacteria. Their localization at 90–100 cm depth reflects both the anaerobic conditions at the bottom of the peat core and the fact that this phylum includes several psychrophilic bacteria. Note that *Bacteroidetes* is one of the dominant phyla in polygonal tundra soils of the Lena Delta [71]. Finally, the deep 90–100 layer of frozen peat also harbored other anaerobic spore-forming bacteria, such as *Clostridia* (4.9%) and *Bacilli* (10.7%).

### 4.3. Multiple Factors Governing the Microbial Diversity in the Peat Core

The chemical composition of peat porewaters and suprapermafrost waters sampled in different sites of the Khanymey region (i.e., mounds, depressions) was shown to be rather similar, and overall in western Siberia, the impact of latitude and permafrost was much stronger than the difference between micro-landscapes (Raudina et al., 2018). Because the fluid phase of the peat layer represents maximal possible variations in chemical composition among soil micro-environments, we believe that the difference in the chemical and microbial composition between adjacent cores is quite small.
cells in the surface layer of peat core, which exceeds the typical values of mineral soils by a factor of 100 to 1000.

In agreement with Jansson and Taş [69], the diversity of bacteria in the permafrost peat assessed in this study is much higher than that of archaea (99% of all sequences in this study). The dominant phylum *Proteobacteria*, *Firmicutes*, *Acidobacteria* and *Actinobacteria* identified in the present study are also reported in other permafrost environments [36,69]. The microbial population of the peat soil of Central Siberia was dominated by bacteria (96.6–98.7% of all sequences) and included *Proteobacteria* (34–25%), *Actinobacteria* (20–18%), *Acidobacteria* (8–14%) and *Verrucomicrobia* (13–14%), as reported by Grodnitskaya et al. [33]. Similarly, polygonal peat bogs of the Canadian Arctic contained *Proteobacteria* (40–50%), *Bacteroidetes* (20–40%) and *Acidobacteria* (10–15%) [71].

In studied peat core, *Proteobacteria* were abundant in all layers and accounted for 47% of all sequences, on average. This was represented by *Alphaproteobacteria* (5–50%) and *Gammaproteobacteria*, with maximal abundance at the mineral layer (84%). The possibility of *Proteobacteria* synthesizing DNA at a temperature from 0 to -20 °C [72] may explain their abundance at the 60–110 cm depth, where the temperature ranged from -0.1 to -0.5 °C. It is noteworthy that *Alphaproteobacteria* are capable of growing after nitrogen addition [35], which can explain the dominance of this class at 70–80 cm depth, where the N concentration achieves its maximum. The main representative of *Alphaproteobacteria* is clade *Sphingomonades*, which is abundant at the depth of 90 to 100 cm (20.3%). This clade is known to degrade lignocellulose debris [73], which is consistent with the botanical composition of studied peat: the wood debris is highly abundant below 80 cm depth. The representatives of clade *Rhizobiales* are highly abundant at 0–20 cm depth. This order includes gram-negative soil bacteria, capable of fixing nitrogen in symbiosis with plants [74–79]. Overall, the near-surface layer of peat is the major location for root exudates, most likely in the form of labile C substances, which are known to favor the growth of copiotrophic microorganisms [23] such as *Gamma*- and *Alphaproteobacteria*, as also evidenced from a study of catabolic activity of peat microorganisms [20].

Class *Acidobacteria* is highly abundant in soils and may constitute up to 52% of the microbial community [80]. In studied peat core, oligotrophic *Acidobacteria* were highly abundant and mostly present in the thawed layer or permafrost boundary (50–60 cm depth) and strongly decreased in the deep-frozen and mineral horizon. This may reflect the need of these bacteria for nutrients, which are concentrating at this boundary. Moreover, the peat water from 50–60 cm depth exhibited the lowest pH of all sampled horizons, consistent with the preference of *Acidobacteria* for acidic environments.

The *Actinobacteria* were present over the full depth of the peat core, which might be due to their ability to assimilate a broad spectrum of organic substrates and their adaptation to low temperature, for example, due to elevated content of guanidine and cytosine [80,81]. The *Bacteroidetes* are the main representatives of gram-negative anaerobic bacteria. Their localization at 90–100 cm depth reflects both the anaerobic conditions at the bottom of the peat core and the fact that this phylum includes several psychrophilic bacteria. Note that *Bacteroidetes* is one of the dominant phyla in polygonal tundra soils of the Lena Delta [71]. Finally, the deep 90–100 layer of frozen peat also harbored other anaerobic spore-forming bacteria, such as *Clostridia* (4.9%) and *Bacilli* (10.7%).

### 4.3. Multiple Factors Governing the Microbial Diversity in the Peat Core

The chemical composition of peat porewaters and suprapermafrost waters sampled in different sites of the Khanymey region (i.e., mounds, depressions) was shown to be rather similar, and overall in western Siberia, the impact of latitude and permafrost was much stronger than the difference between micro-landscapes (Raudina et al., 2018). Because the fluid phase of the peat layer represents maximal possible variations in chemical composition among soil micro-environments, we believe that the difference in the chemical and microbial composition between adjacent cores is quite small.
Overall, in our approach to microbial diversity, we followed the methodology of a previous study of the peat core from the Khanymey site of the WSL (Morgalev et al., 2017). The peat on the watershed divide of the Khanymey Research Station is highly homogeneous.

However, we admit that a special study of microbial diversity in peat cores sampled in different micro-landscapes and across the permafrost gradient is needed to extend the obtained results to a much larger territory of western Siberia.

The decrease in the bacterial number (this study) and AWCD [20] with depth, primarily within the active layer, probably occurred due to a decrease in the number of consumed substrates. In the 0–20 cm upper layer of peat, which contained lichen debris, the most metabolically active bacteria were those that consumed amines such as *Proteobacteria* [82]. In a previous study of metabolic activity of peat microorganisms, we demonstrated a peak of activity at the border of the thawed layer and permafrost [20]. It is not excluded that the maximal variety of *Acidobacteria*, also occurring at this layer (50–60 cm) in this study, corresponds to the maximal capacity of these oligotrophic bacteria to process organic substrates. It could not be proven, however, that the peak of metabolic activity and the highest diversity of *Acidobacteria* at the permafrost–thaw layer boundary are linked to the specificity of organic carbon quality in this layer or the chemical composition of peat pore water, since no direct correlation between these parameters has been established.

An increase in Ni concentration in peat ice below 70–80 cm depth might be consistent with its involvement in the metabolism of anaerobic methane-producing bacteria, which become abundant at the bottom of the peat core. The deep layers of strongly degraded peat contain elevated B, Cr, Cu and Ba concentrations, aromatic C, carboxylic C and alkyl C and are also enriched in DOC, Ca, Mg, K, Fe, Si, Zn, Co, Ni, Mn, V and Rb in the liquid (interstitial water and ice) phase. The latter favored the elevated proportions of bacterial classes *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Chloroflexia* and *Archaea Thau-marchaeota*, and *Nitrososphaeria*. Several classes of bacteria demonstrated a tendency to avoid the medium part of the peat profile, located just below the active layer boundary, where there were minimal concentrations of lithogenic elements linked to mineral silicate dust. These were *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteriia*. Finally, a single group of *Actinobacteria* was associated with minimal N concentration and the highest δ¹⁵N. Presumably, it actively consumed N-rich compounds for guanine and cytosine synthesis, which is consistent with a metabolic diversity of the studied peat column [20]. It is also possible that, after cell lysis, the bacterial N is removed and translocated from the relevant peat layer.

Placing the obtained results in the context of permafrost thaw and active layer deepening, we hypothesize that freezing front propagation downwards leads to mobilization of essential macro- and micro-nutrients at the active layer boundary; as a result, the layer of maximal genetic diversity of the microbial population will move down the soil profile. Moreover, the total microbial number strongly decreased (×10) from the surface to the active layer and further to the mineral horizon (×10). We, therefore, estimate that, under a climate warming and permafrost thaw scenario, in the case of full thawing of the frozen peat and underlying mineral horizons, the total number of bacteria over the 100 cm thick layer may significantly increase.

**5. Conclusions**

Unlike the total bacterial number, which decreased almost two orders of magnitude from the surface down to mineral frozen layer of the peat core, the genetic diversity of microorganisms increased with depth. Moreover, consistent with previously reported peaks of the intensity of substrate utilization right at the boundary of the active layer and frozen peat, we evidenced a dominance of *Acidobacteria* in the genetic diversity of microorganisms of the peat core at the thawed peat–permafrost boundary.
Overall, in our approach to microbial diversity, we followed the methodology of a previous study of the peat core from the Khanymey site of the WSL (Morgalev et al., 2017). The peat on the watershed divide of the Khanymey Research Station is highly homogeneous.

However, we admit that a special study of microbial diversity in peat cores sampled in different micro-landscapes and across the permafrost gradient is needed to extend the obtained results to a much larger territory of western Siberia.

The decrease in the bacterial number (this study) and AWCD [20] with depth, primarily within the active layer, probably occurred due to a decrease in the number of consumed substrates. In the 0–20 cm upper layer of peat, which contained lichen debris, the most metabolically active bacteria were those that consumed amines such as *Proteobacteria* [82]. In a previous study of metabolic activity of peat microorganisms, we demonstrated a peak of activity at the border of the thawed layer and permafrost [20]. It is not excluded that the maximal variety of *Acidobacteria*, also occurring at this layer (50–60 cm) in this study, corresponds to the maximal capacity of these oligotrophic bacteria to process organic substrates. It could not be proven, however, that the peak of metabolic activity and the highest diversity of *Acidobacteria* at the permafrost–thaw layer boundary are linked to the specificity of organic carbon quality in this layer or the chemical composition of peat pore water, since no direct correlation between these parameters has been established.

An increase in Ni concentration in peat ice below 70–80 cm depth might be consistent with its involvement in the metabolism of anaerobic methane-producing bacteria, which become abundant at the bottom of the peat core. The deep layers of strongly degraded peat contain elevated B, Cr, Cu and Ba concentrations, aromatic C, carboxylic C and alkyl C and are also enriched in DOC, Ca, Mg, K, Fe, Si, Zn, Co, Ni, Mn, V and Rb in the liquid (interstitial water and ice) phase. The latter favored the elevated proportions of bacterial classes *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Chloroflexia* and *Archaea Thau-marchaeota*, and *Nitrososphaera*. Several classes of bacteria demonstrated a tendency to avoid the medium part of the peat profile, located just below the active layer boundary, where there were minimal concentrations of lithogenic elements linked to mineral silicate dust. These were *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteria*. Finally, a single group of *Actinobacteria* was associated with minimal N concentration and the highest $\delta^{15}$N. Presumably, it actively consumed N-rich compounds for guanine and cytosine synthesis, which is consistent with a metabolic diversity of the studied peat column [20]. It is also possible that, after cell lysis, the bacterial N is removed and translocated from the relevant peat layer.

Placing the obtained results in the context of permafrost thaw and active layer deepening, we hypothesize that freezing front propagation downwards leads to mobilization of essential macro- and micro-nutrients at the active layer boundary; as a result, the layer of maximal genetic diversity of the microbial population will move down the soil profile. Moreover, the total microbial number strongly decreased (×10) from the surface to the active layer and further to the mineral horizon (×10). We, therefore, estimate that, under a climate warming and permafrost thaw scenario, in the case of full thawing of the frozen peat and underlying mineral horizons, the total number of bacteria over the 100 cm thick layer may significantly increase.

**5. Conclusions**

Unlike the total bacterial number, which decreased almost two orders of magnitude from the surface down to mineral frozen layer of the peat core, the genetic diversity of microorganisms increased with depth. Moreover, consistent with previously reported peaks of the intensity of substrate utilization right at the boundary of the active layer and frozen peat, we evidenced a dominance of *Acidobacteria* in the genetic diversity of microorganisms of the peat core at the thawed peat–permafrost boundary.
A multi-parametric statistical analysis demonstrated that three factors were responsible for the main variability of chemical and microbial parameter distribution. There was a strongly pronounced preference of *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Thaumarchaeota*, *Nitrososphaeria* and *Chloroflexia* for deeper peat horizons, which were enriched in alkyl C, carboxyl C, aromatic C, B, Cr, and Cu in the solid phase and DOC and various macro- and micro-nutrients in the fluid (peat porewater and ice) phase. In contrast, classes *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteriia* preferentially avoided the intermediate part of the peat core, where there was a minimal concentration of mineral dust. *Actinobacteria* demonstrated a strong link with the highest C:N ratio and $\delta^{15}$N value.

Given the largest diversity of microorganisms in the bottom organic (80–100 cm) and mineral (100–110 cm) layers, the frozen peat and underlying mineral substrates of western Siberia present a high potential for microbial diversity; in addition to its sizeable nutrient (mineral and organic) resources, this area may represent an important source of active microorganisms under continuous permafrost thaw.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/d13070328/s1, Table S1: Chemical composition of peat interstitial water and ice; Table S2: Number of microorganisms as a function of peat depth; Table S3: Relative proportion of order units of peat microorganisms over the peat column; Table S4: Family units of peat microorganisms; Table S5: Genus units of peat microorganisms; Table S6: Most pronounced pairwise correlations between class of microorganisms and peat and peat water/ice chemical components.

**Author Contributions:** Conceptualization, O.S.P., A.S.A. and L.S.S.; writing—original draft preparation, O.S.P.; sample collection, A.G.L., D.K. and S.V.L.; organic matter quality and isotopic composition, S.P., M.A.A., M.C.-R. and S.V.L.; molecular genetic analysis, O.Y.K.; writing—review and editing (microbiology part), S.N.K.; visualization and supervision (genetic part), A.S.A.; chemical analyses, A.G.L., S.V.L., D.K. and L.S.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research was funded by the Government of the Russian Federation in compliance with the Resolution of 09 April 2010 No. 220 (the contract from 14.03.2017 No. 14.Y26.31.0009). The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** The research was done using equipment of the Core Centrum 'Genomic Technologies, Proteomics and Cell Biology' in ARRIAM.

**Conflicts of Interest:** The authors declare that there is no conflict of interest related to the publication of this article, which should be reported.

## References

1. Wilhelm, R.C.; Niederberger, T.D.; Greer, C.; Whyte, L.G. Microbial diversity of active layer and permafrost in an acidic wetland from the Canadian High Arctic. *Can. J. Microbiol.* **2011**, *57*, 303–315. [CrossRef]

2. Steven, B.; Briggs, G.; McKay, C.P.; Pollard, W.H.; Greer, C.W.; Whyte, L.G. Characterization of the microbial diversity in a permafrost sample from the Canadian high Arctic using culture-dependent and culture-independent methods. *FEMS Microbiol. Ecol.* **2007**, *59*, 513–523. [CrossRef]

3. Steven, B.; Pollard, W.H.; Greer, C.W.; Whyte, L.G. Microbial diversity and activity through a permafrost/ground ice core profile from the Canadian high Arctic. *Environ. Microbiol.* **2008**, *10*, 3388–3403. [CrossRef] [PubMed]

4. Hansen, A.A.; Herbert, R.A.; Mikkelsen, K.; Jensen, L.L.; Kristoffersen, T.; Tiedje, J.M.; Lomstein, B.A.; Finster, K. Viability, diversity and composition of the bacterial community in a high Arctic permafrost soil from Spitsbergen, Northern Norway. *Environ. Microbiol.* **2007**, *9*, 2870–2884. [CrossRef] [PubMed]
A multi-parametric statistical analysis demonstrated that three factors were responsible for the main variability of chemical and microbial parameter distribution. There was a strongly pronounced preference of *Clostridia*, *Bacilli*, *Oxyphotobacteria*, *Gemmatimonadetes*, *Thaumarchaeota*, *Nitrososphaeria* and *Chloroflexia* for deeper peat horizons, which were enriched in alkyl C, carboxyl C, aromatic C, B, Cr, and Cu in the solid phase and DOC and various macro- and micro-nutrients in the fluid (peat porewater and ice) phase. In contrast, classes *Verrucomicrobiae*, *Deltaproteobacteria*, *Fimbriimonadia*, *Gammaproteobacteria*, *Phycisphaerae*, *Planctomycetacia*, *Chthonomonadetes*, *Armatimonadia*, *Bacteroidia* and *Acidobacteriia* preferentially avoided the intermediate part of the peat core, where there was a minimal concentration of mineral dust. *Actinobacteria* demonstrated a strong link with the highest C:N ratio and δ¹⁵N value.

Given the largest diversity of microorganisms in the bottom organic (80–100 cm) and mineral (100–110 cm) layers, the frozen peat and underlying mineral substrates of western Siberia present a high potential for microbial diversity; in addition to its sizeable nutrient (mineral and organic) resources, this area may represent an important source of active microorganisms under continuous permafrost thaw.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/d13070328/s1, Table S1: Chemical composition of peat interstitial water and ice; Table S2: Number of microorganisms as a function of peat depth; Table S3: Relative proportion of order units of peat microorganisms over the peat column; Table S4: Family units of peat microorganisms; Table S5: Genus units of peat microorganisms; Table S6: Most pronounced pairwise correlations between class of microorganisms and peat and peat water/ice chemical components.

**Author Contributions:** Conceptualization, O.S.P., A.S.A. and L.S.S.; writing—original draft preparation, O.S.P.; sample collection, A.G.L., D.K. and S.V.L.; organic matter quality and isotopic composition, S.P., M.A.A., M.C.-R. and S.V.L.; molecular genetic analysis, O.Y.K.; writing—review and editing (microbiology part), S.N.K.; visualization and supervision (genetic part), A.S.A.; chemical analyses, A.G.L., S.V.L., D.K. and L.S.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research was funded by the Government of the Russian Federation in compliance with the Resolution of 09 April 2010 No. 220 (the contract from 14.03.2017 No. 14.Y26.31.0009). The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** The research was done using equipment of the Core Centrum 'Genomic Technologies, Proteomics and Cell Biology' in ARRIAM.

**Conflicts of Interest:** The authors declare that there is no conflict of interest related to the publication of this article, which should be reported.

## References

1. Wilhelm, R.C.; Niederberger, T.D.; Greer, C.; Whyte, L.G. Microbial diversity of active layer and permafrost in an acidic wetland from the Canadian High Arctic. *Can. J. Microbiol.* **2011**, *57*, 303–315. [CrossRef]

2. Steven, B.; Briggs, G.; McKay, C.P.; Pollard, W.H.; Greer, C.W.; Whyte, L.G. Characterization of the microbial diversity in a permafrost sample from the Canadian high Arctic using culture-dependent and culture-independent methods. *FEMS Microbiol. Ecol.* **2007**, *59*, 513–523. [CrossRef]

3. Steven, B.; Pollard, W.H.; Greer, C.W.; Whyte, L.G. Microbial diversity and activity through a permafrost/ground ice core profile from the Canadian high Arctic. *Environ. Microbiol.* **2008**, *10*, 3388–3403. [CrossRef] [PubMed]

4. Hansen, A.A.; Herbert, R.A.; Mikkelsen, K.; Jensen, L.L.; Kristoffersen, T.; Tiedje, J.M.; Lomstein, B.A.; Finster, K. Viability, diversity and composition of the bacterial community in a high Arctic permafrost soil from Spitsbergen, Northern Norway. *Environ. Microbiol.* **2007**, *9*, 2870–2884. [CrossRef] [PubMed]
5. Vorobyova, E.; Soina, V.; Gorlenko, M.; Minkovskaya, N.; Zalinova, N.; Mamukelashvili, A.; Gilichinsky, D.; Rivkina, E.; Vishnivetskaya, T. The deep cold biosphere: Facts and hypothesis. *FEMS Microbiol. Rev.* **1997**, *20*, 277–290. [CrossRef]

6. Schnecker, J.; Wild, B.; Hofhansl, F.; Alves, R.J.E.; Bárta, J.; Capek, P.; Fuchslueger, L.; Gentsch, N.; Gittel, A.; Guggenberger, G.; et al. Effects of Soil Organic Matter Properties and Microbial Community Composition on Enzyme Activities in Cryoturbated Arctic Soils. *PLoS ONE* **2014**, *9*, e94076. [CrossRef] [PubMed]

7. Wild, B.; Schnecker, J.; Alves, R.J.E.; Barsukov, P.; Bárta, J.; Capek, P.; Gentsch, N.; Gittel, A.; Guggenberger, G.; Lashchinskiy, N.; et al. Input of easily available organic C and N stimulates microbial decomposition of soil organic matter in arctic permafrost soil. *Soil Biol. Biochem.* **2014**, *75*, 143–151. [CrossRef]

8. Wallenstein, M.D.; McMahon, S.; Schimel, J. Bacterial and fungal community structure in Arctic tundra tussock and shrub soils. *FEMS Microbiol. Ecol.* **2007**, *59*, 428–435. [CrossRef] [PubMed]

9. Mackelprang, R.; Waldrop, M.; DeAngelis, K.; David, M.M.; Chavarria, K.L.; Blazewicz, S.J.; Rubin, E.M.; Jansson, J. Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw. *Nature* **2011**, *480*, 368–371. [CrossRef]

10. Gittel, A.; Bárta, J.; Kohoutová, I.; Mikutta, R.; Owens, S.; Gilbert, J.; Schnecker, J.; Wild, B.; Hannisdal, B.; Maerz, J.; et al. Distinct microbial communities associated with buried soils in the Siberian tundra. *ISME J.* **2013**, *8*, 841–853. [CrossRef] [PubMed]

11. Yergeau, E.; Bokhorst, S.; Huiskes, A.H.L.; Boschker, H.T.S.; Aerts, R.; Kowalchuk, G.A. Size and structure of bacterial, fungal and nematode communities along an Antarctic environmental gradient. *FEMS Microbiol. Ecol.* **2006**, *59*, 436–451. [CrossRef]

12. Yergeau, E.; Kang, S.; He, Z.; Zhou, J.; Kowalchuk, G.A. Functional microarray analysis of nitrogen and carbon cycling genes across an Antarctic latitudinal transect. *ISME J.* **2007**, *1*, 163–179. [CrossRef]

13. Yergeau, E.; Schoondermark-Stolk, S.A.; Brodie, E.L.; Déjean, S.; DeSantis, T.Z.; Gonçalves, O.; Piceno, Y.M.; Andersen, G.L.; Kowalchuk, G.A. Environmental microarray analyses of Antarctic soil microbial communities. *ISME J.* **2009**, *3*, 340–351. [CrossRef]

14. Kremenetski, K.; Velichko, A.; Borisova, O.; MacDonald, G.; Smith, L.; Frey, K.; Orlova, L. Peatlands of the Western Siberian lowlands: Current knowledge on zonation, carbon content and Late Quaternary history. *Quat. Sci. Rev.* **2003**, *22*, 703–723. [CrossRef]

15. Tyrtikov, A.P. Thawing of soils in tundra of western Siberia. In *Natural Environment of Western Siberia*; Popov, A.I., Ed.; MG Izd-vo: Moscow, Russia, 1973; pp. 160–169.

16. Tyrtikov, A.P. *Dynamics of Vegetation Coverage and Permafrost Development in Western Siberia*; Nauka: Moscow, Russia, 1979; 114p.

17. Vasil'evskaya, V.D.; Ivanov, V.V.; Bogatyrev, L.G. *Soils of North of Western Siberia*; Moscow University Publ. House: Moscow, Russia, 1986; 228p.

18. Novikov, S.M.; Moskvin, Y.P.; Trofimov, S.A.; Usova, L.I.; Batuev, V.I.; Tumanovskaya, S.M.; Smirnova, V.P.; Markov, M.L.; Korotkevich, A.E.; Potapova, T.M. *Hydrology of Bog Territories of the Permafrost Zone of Western Siberia*; BBM publ. House: St. Petersburg, Russia, 2009; 535p.

19. Khrenov, V.Y. *Soils of Cryolithozone of Western Siberia: Morphology, Physico-Chemical Properties and Geochemistry*; Nauka: Moscow, Russia, 2011; 214p.

20. Morgalev, Y.N.; Lushchaeva, I.V.; Morgaleva, T.G.; Kolesnichenko, L.G.; Loiko, S.; Krickov, I.V.; Morgalev, S.Y.; Lim, A.; Raudina, T.V.; Volkova, I.I.; et al. Bacteria primarily metabolize at the active layer/permafrost border in the peat core from a permafrost region in western Siberia. *Polar Biol.* **2017**, *40*, 1645–1659. [CrossRef]

21. Yergeau, E.; Hogues, H.; Whyte, L.G.; Greer, C.W. The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. *ISME J.* **2010**, *4*, 1206–1214. [CrossRef]

22. Gittel, A.; Barta, J.; Lacmanova, I.; Schnecker, J.; Wild, B.; Capek, P.; Kaiser, C.; Torsvik, V.; Richter, A.; Schleper, C.; et al. Site- and horizon-specific patterns of microbial community structure and enzyme activities in permafrost-affected soils of Greenland. *Front. Microbiol.* **2014**, *5*, 541. [CrossRef]

23. Deng, J.; Gu, Y.; Zhang, J.; Xue, K.; Qin, Y.; Yuan, M.; Yin, H.; He, Z.; Wu, L.; Schuur, E.A.G.; et al. Shifts of tundra bacterial and archaeal communities along a permafrost thaw gradient in A laska. *Mol. Ecol.* **2015**, *24*, 222–234. [CrossRef] [PubMed]

24. Johnston, E.; Rodriguez-R, L.M.; Eluo, C.; Yuan, M.M.; Ewu, L.; He, Z.; Konstantinidis, K.T.; Luo, Y.; Tiedje, J.M.; Zhou, J.; et al. Metagenomics Reveals Pervasive Bacterial Populations and Reduced Community Diversity across the Alaska Tundra Ecosystem. *Front. Microbiol.* **2016**, *7*, 579. [CrossRef] [PubMed]

25. Inglese, C.N.; Christiansen, C.; Lamhonwah, D.; Moniz, K.; Montross, S.N.; Lamoureux, S.; Lafrenière, M.; Grogan, P.; Walker, V.K. Examination of Soil Microbial Communities After Permafrost Thaw Subsequent to an Active Layer Detachment in the High Arctic. *Arct. Antarct. Alp. Res.* **2017**, *49*, 455–472. [CrossRef]

26. Malard, L.A.; Pearce, D.A. Microbial diversity and biogeography in Arctic soils. *Environ. Microbiol. Rep.* **2018**, *10*, 611–625. [CrossRef]

27. Monteux, S.; Weedon, J.T.; Blume-Werry, G.; Gavazov, K.; Jassey, V.E.J.; Johansson, M.; Keuper, F.; Olid, C.; Dorrepaal, E. Long-term in situ permafrost thaw effects on bacterial communities and potential aerobic respiration. *ISME J.* **2018**, *12*, 2129–2141. [CrossRef] [PubMed]

28. Malard, L.A.; Anwar, M.Z.; Jacobsen, C.S.; Pearce, D.A. Biogeographical patterns in soil bacterial communities across the Arctic region. *FEMS Microbiol. Ecol.* **2019**, *95*, fiz128. [CrossRef] [PubMed]
5. Vorobyova, E.; Soina, V.; Gorlenko, M.; Minkovskaya, N.; Zalinova, N.; Mamukelashvili, A.; Gilichinsky, D.; Rivkina, E.; Vishnivetskaya, T. The deep cold biosphere: Facts and hypothesis. *FEMS Microbiol. Rev.* **1997**, *20*, 277–290. [CrossRef]

6. Schnecker, J.; Wild, B.; Hofhansl, F.; Alves, R.J.E.; Bárta, J.; Capek, P.; Fuchslueger, L.; Gentsch, N.; Gittel, A.; Guggenberger, G.; et al. Effects of Soil Organic Matter Properties and Microbial Community Composition on Enzyme Activities in Cryoturbated Arctic Soils. *PLoS ONE* **2014**, *9*, e94076. [CrossRef] [PubMed]

7. Wild, B.; Schnecker, J.; Alves, R.J.E.; Barsukov, P.; Bárta, J.; Capek, P.; Gentsch, N.; Gittel, A.; Guggenberger, G.; Lashchinskiy, N.; et al. Input of easily available organic C and N stimulates microbial decomposition of soil organic matter in arctic permafrost soil. *Soil Biol. Biochem.* **2014**, *75*, 143–151. [CrossRef]

8. Wallenstein, M.D.; McMahon, S.; Schimel, J. Bacterial and fungal community structure in Arctic tundra tussock and shrub soils. *FEMS Microbiol. Ecol.* **2007**, *59*, 428–435. [CrossRef] [PubMed]

9. Mackelprang, R.; Waldrop, M.; DeAngelis, K.; David, M.M.; Chavarria, K.L.; Blazewicz, S.J.; Rubin, E.M.; Jansson, J. Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw. *Nature* **2011**, *480*, 368–371. [CrossRef]

10. Gittel, A.; Bárta, J.; Kohoutová, I.; Mikutta, R.; Owens, S.; Gilbert, J.; Schnecker, J.; Wild, B.; Hannisdal, B.; Maerz, J.; et al. Distinct microbial communities associated with buried soils in the Siberian tundra. *ISME J.* **2013**, *8*, 841–853. [CrossRef] [PubMed]

11. Yergeau, E.; Bokhorst, S.; Huiskes, A.H.L.; Boschker, H.T.S.; Aerts, R.; Kowalchuk, G.A. Size and structure of bacterial, fungal and nematode communities along an Antarctic environmental gradient. *FEMS Microbiol. Ecol.* **2006**, *59*, 436–451. [CrossRef]

12. Yergeau, E.; Kang, S.; He, Z.; Zhou, J.; Kowalchuk, G.A. Functional microarray analysis of nitrogen and carbon cycling genes across an Antarctic latitudinal transect. *ISME J.* **2007**, *1*, 163–179. [CrossRef]

13. Yergeau, E.; Schoondermark-Stolk, S.A.; Brodie, E.L.; Déjean, S.; DeSantis, T.Z.; Gonçalves, O.; Piceno, Y.M.; Andersen, G.L.; Kowalchuk, G.A. Environmental microarray analyses of Antarctic soil microbial communities. *ISME J.* **2009**, *3*, 340–351. [CrossRef]

14. Kremenetski, K.; Velichko, A.; Borisova, O.; MacDonald, G.; Smith, L.; Frey, K.; Orlova, L. Peatlands of the Western Siberian lowlands: Current knowledge on zonation, carbon content and Late Quaternary history. *Quat. Sci. Rev.* **2003**, *22*, 703–723. [CrossRef]

15. Tyrtikov, A.P. Thawing of soils in tundra of western Siberia. In *Natural Environment of Western Siberia*; Popov, A.I., Ed.; MG Izd-vo: Moscow, Russia, 1973; pp. 160–169.

16. Tyrtikov, A.P. *Dynamics of Vegetation Coverage and Permafrost Development in Western Siberia*; Nauka: Moscow, Russia, 1979; 114p.

17. Vasil'evskaya, V.D.; Ivanov, V.V.; Bogatyrev, L.G. *Soils of North of Western Siberia*; Moscow University Publ. House: Moscow, Russia, 1986; 228p.

18. Novikov, S.M.; Moskvin, Y.P.; Trofimov, S.A.; Usova, L.I.; Batuev, V.I.; Tumanovskaya, S.M.; Smirnova, V.P.; Markov, M.L.; Korotkevich, A.E.; Potapova, T.M. *Hydrology of Bog Territories of the Permafrost Zone of Western Siberia*; BBM publ. House: St. Petersburg, Russia, 2009; 535p.

19. Khrenov, V.Y. *Soils of Cryolithozone of Western Siberia: Morphology, Physico-Chemical Properties and Geochemistry*; Nauka: Moscow, Russia, 2011; 214p.

20. Morgalev, Y.N.; Lushchaeva, I.V.; Morgaleva, T.G.; Kolesnichenko, L.G.; Loiko, S.; Krickov, I.V.; Morgalev, S.Y.; Lim, A.; Raudina, T.V.; Volkova, I.I.; et al. Bacteria primarily metabolize at the active layer/permafrost border in the peat core from a permafrost region in western Siberia. *Polar Biol.* **2017**, *40*, 1645–1659. [CrossRef]

21. Yergeau, E.; Hogues, H.; Whyte, L.G.; Greer, C.W. The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. *ISME J.* **2010**, *4*, 1206–1214. [CrossRef]

22. Gittel, A.; Barta, J.; Lacmanova, I.; Schnecker, J.; Wild, B.; Capek, P.; Kaiser, C.; Torsvik, V.; Richter, A.; Schleper, C.; et al. Site- and horizon-specific patterns of microbial community structure and enzyme activities in permafrost-affected soils of Greenland. *Front. Microbiol.* **2014**, *5*, 541. [CrossRef]

23. Deng, J.; Gu, Y.; Zhang, J.; Xue, K.; Qin, Y.; Yuan, M.; Yin, H.; He, Z.; Wu, L.; Schuur, E.A.G.; et al. Shifts of tundra bacterial and archaeal communities along a permafrost thaw gradient in A laska. *Mol. Ecol.* **2015**, *24*, 222–234. [CrossRef] [PubMed]

24. Johnston, E.; Rodriguez-R, L.M.; Eluo, C.; Yuan, M.M.; Ewu, L.; He, Z.; Konstantinidis, K.T.; Luo, Y.; Tiedje, J.M.; Zhou, J.; et al. Metagenomics Reveals Pervasive Bacterial Populations and Reduced Community Diversity across the Alaska Tundra Ecosystem. *Front. Microbiol.* **2016**, *7*, 579. [CrossRef] [PubMed]

25. Inglese, C.N.; Christiansen, C.; Lamhonwah, D.; Moniz, K.; Montross, S.N.; Lamoureux, S.; Lafrenière, M.; Grogan, P.; Walker, V.K. Examination of Soil Microbial Communities After Permafrost Thaw Subsequent to an Active Layer Detachment in the High Arctic. *Arct. Antarct. Alp. Res.* **2017**, *49*, 455–472. [CrossRef]

26. Malard, L.A.; Pearce, D.A. Microbial diversity and biogeography in Arctic soils. *Environ. Microbiol. Rep.* **2018**, *10*, 611–625. [CrossRef]

27. Monteux, S.; Weedon, J.T.; Blume-Werry, G.; Gavazov, K.; Jassey, V.E.J.; Johansson, M.; Keuper, F.; Olid, C.; Dorrepaal, E. Long-term in situ permafrost thaw effects on bacterial communities and potential aerobic respiration. *ISME J.* **2018**, *12*, 2129–2141. [CrossRef] [PubMed]

28. Malard, L.A.; Anwar, M.Z.; Jacobsen, C.S.; Pearce, D.A. Biogeographical patterns in soil bacterial communities across the Arctic region. *FEMS Microbiol. Ecol.* **2019**, *95*, fiz128. [CrossRef] [PubMed]
29. Xue, Y.; Jonassen, I.; Øvreås, L.; Taş, N. Metagenome-assembled genome distribution and key functionality highlight importance of aerobic metabolism in Svalbard permafrost. *FEMS Microbiol. Ecol.* **2020**, *96*, fiaa057. [CrossRef]

30. Bland, A.; Lerch, T.Z.; Phoenix, G.K.; Osborn, A.M. Arctic soil microbial diversity in a changing world. *Res. Microbiol.* **2015**, *166*, 796–813. [CrossRef]

31. Lipson, D.A.; Haggerty, J.M.; Srinivas, A.; Raab, T.K.; Sathe, S.; Dinsdale, E.A. Metagenomic Insights into Anaerobic Metabolism along an Arctic Peat Soil Profile. *PLoS ONE* **2013**, *8*, e64659. [CrossRef]

32. Ivanova, A.A.; Kulichevskaya, I.S.; Merkel, A.Y.; Toshchakov, S.; Dedysh, S.N. High Diversity of Planctomycetes in Soils of Two Lichen-Dominated Sub-Arctic Ecosystems of Northwestern Siberia. *Front. Microbiol.* **2016**, *7*, 2065. [CrossRef] [PubMed]

33. Grodnitskaya, I.D.; Trusova, M.Y.; Syrtsov, S.N.; Koroban, N.V. Structure of microbial communities of peat soils in two bogs in Siberian tundra and forest zones. *Microbiology* **2018**, *87*, 89–102. [CrossRef]

34. Schnecker, J.; Wild, B.; Takriti, M.; Alves, R.J.E.; Gentsch, N.; Gittel, A.; Mikutta, R.; Richter, A. Microbial community composition shapes enzyme patterns in topsoil and subsoil horizons along a latitudinal transect in Western Siberia. *Soil Biol. Biochem.* **2015**, *83*, 106–115. [CrossRef]

35. Kim, H.M.; Jung, J.Y.; Yergeau, E.; Hwang, C.; Hinzman, L.; Nam, S.; Hong, S.G.; Kim, O.-S.; Chun, J.; Lee, Y.K. Bacterial community structure and soil properties of a subarctic tundra soil in Council, Alaska. *FEMS Microbiol. Ecol.* **2014**, *89*, 465–475. [CrossRef]

36. Tveit, A.; Schwacke, R.; Svenning, M.M.; Urich, T. Organic carbon transformations in high-Arctic peat soils: Key functions and microorganisms. *ISME J.* **2013**, *7*, 299–311. [CrossRef] [PubMed]

37. Tveit, A.T.; Urich, T.; Frenzel, P.; Svenning, M.M. Metabolic and trophic interactions modulate methane production by Arctic peat microbiota in response to warming. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, E2507–E2516. [CrossRef] [PubMed]

38. Martineau, C.; Pan, Y.; Bodrossy, L.; Yergeau, E.; Whyte, L.G.; Greer, C.W. Atmospheric methane oxidizers are present and active in Canadian high Arctic soils. *FEMS Microbiol. Ecol.* **2014**, *89*, 257–269. [CrossRef] [PubMed]

39. Lee, S.-H.; Jang, I.; Chae, N.; Choi, T.; Kang, H. Organic Layer Serves as a Hotspot of Microbial Activity and Abundance in Arctic Tundra Soils. *Microb. Ecol.* **2012**, *65*, 405–414. [CrossRef] [PubMed]

40. Taş, N.; Prestat, E.; Wang, S.; Wu, Y.; Ulrich, C.; Kneafsey, T.; Tringe, S.G.; Torn, M.; Hubbard, S.S.; Jansson, J.K. Landscape topography structures the soil microbiome in arctic polygonal tundra. *Nat. Commun.* **2018**, *9*, 777. [CrossRef] [PubMed]

41. Vishnivetskaya, T.A.; Buongiorno, J.; Bird, J.; Krivushin, K.; Spirina, E.V.; Oshurkova, V.; Shcherbakova, V.A.; Wilson, G.; Lloyd, K.G.; Rivkina, E.M. Methanogens in the Antarctic Dry Valley permafrost. *FEMS Microbiol. Ecol.* **2018**, *94*, 109. [CrossRef] [PubMed]

42. Trotsenko, Y.A.; Khmelenina, V.N. Aerobic methanotrophic bacteria of cold ecosystems. *FEMS Microbiol. Ecol.* **2005**, *53*, 15–26. [CrossRef] [PubMed]

43. Shcherbakova, V.; Yoshimura, Y.; Ryzhmanova, Y.; Taguchi, Y.; Segawa, T.; Oshurkova, V.; Rivkina, E. Archaeal communities of Arctic methane-containing permafrost. *FEMS Microbiol. Ecol.* **2016**, *92*, fiw135. [CrossRef]

44. Comte, J.; Monier, A.; Crevecoeur, S.; Lovejoy, C.; Vincent, W.F. Microbial biogeography of permafrost thaw ponds across the changing northern landscape. *Ecography* **2015**, *39*, 609–618. [CrossRef]

45. Roiha, T.; Peura, S.; Cusson, M.; Rautio, M. Allochthonous carbon is a major regulator to bacterial growth and community composition in subarctic freshwaters. *Sci. Rep.* **2016**, *6*, 34456. [CrossRef]

46. Peura, S.; Wauthy, M.; Simone, D.; Eiler, A.; Einarsdóttir, K.; Rautio, M.; Bertilsson, S. Ontogenic succession of thermokarst thaw ponds is linked to dissolved organic matter quality and microbial degradation potential. *Linnol. Oceanogr.* **2019**, *65*, S248–S263. [CrossRef]

47. Przytulska, A.; Comte, J.; Crevecoeur, S.; Lovejoy, C.; Laurion, I.; Vincent, W. Phototrophic pigment diversity and picophyto-plankton in permafrost thaw lakes. *Biogeosciences* **2016**, *13*, 13–26. [CrossRef]

48. Steven, B.; Léveillé, R.; Pollard, W.H.; Whyte, L.G. Microbial ecology and biodiversity in permafrost. *Extremophiles* **2006**, *10*, 259–267. [CrossRef]

49. Liebner, S.; Harder, J.; Wagner, D. Bacterial diversity and community structure in polygonal tundra soils from Samoylov Island, Lena Delta, Siberia. *Int. J. Microbiol.* **2008**, *1*, 195–202.

50. Barbier, B.A.; Dziduch, I.; Liebner, S.; Ganzert, L.; Lantuit, H.; Pollard, W.; Wagner, D. Methane-cycling communities in a permafrost-affected soil on Herschel Island, Western Canadian Arctic: Active layer profiling ofmcrAandpmoAgenes. *FEMS Microbiol. Ecol.* **2012**, *82*, 287–302. [CrossRef]

51. Ganzert, L.; Bajerski, F.; Wagner, D. Bacterial community composition and diversity of five different permafrost-affected soils of Northeast Greenland. *FEMS Microbiol. Ecol.* **2014**, *89*, 426–441. [CrossRef]

52. Zhou, L.; Zhou, Y.; Yao, X.; Cai, J.; Liu, X.; Tang, X.; Zhang, Y.; Jang, K.S.; Jeppesen, E. Decreasing diversity of rare bacterial subcommunities relates to dissolved organic matter along permafrost thawing gradients. *Environ. Int.* **2020**, *134*, 105330. [CrossRef] [PubMed]

53. Coolen, M.J.L.; Van De Giessen, J.; Zhu, E.Y.; Wuchter, C. Bioavailability of soil organic matter and microbial community dynamics upon permafrost thaw. *Environ. Microbiol.* **2011**, *13*, 2299–2314. [CrossRef] [PubMed]
29. Xue, Y.; Jonassen, I.; Øvreås, L.; Taş, N. Metagenome-assembled genome distribution and key functionality highlight importance of aerobic metabolism in Svalbard permafrost. *FEMS Microbiol. Ecol.* **2020**, *96*, fiaa057. [CrossRef]

30. Bland, A.; Lerch, T.Z.; Phoenix, G.K.; Osborn, A.M. Arctic soil microbial diversity in a changing world. *Res. Microbiol.* **2015**, *166*, 796–813. [CrossRef]

31. Lipson, D.A.; Haggerty, J.M.; Srinivas, A.; Raab, T.K.; Sathe, S.; Dinsdale, E.A. Metagenomic Insights into Anaerobic Metabolism along an Arctic Peat Soil Profile. *PLoS ONE* **2013**, *8*, e64659. [CrossRef]

32. Ivanova, A.A.; Kulichevskaya, I.S.; Merkel, A.Y.; Toshchakov, S.; Dedysh, S.N. High Diversity of Planctomycetes in Soils of Two Lichen-Dominated Sub-Arctic Ecosystems of Northwestern Siberia. *Front. Microbiol.* **2016**, *7*, 2065. [CrossRef] [PubMed]

33. Grodnitskaya, I.D.; Trusova, M.Y.; Syrtsov, S.N.; Koroban, N.V. Structure of microbial communities of peat soils in two bogs in Siberian tundra and forest zones. *Microbiology* **2018**, *87*, 89–102. [CrossRef]

34. Schnecker, J.; Wild, B.; Takriti, M.; Alves, R.J.E.; Gentsch, N.; Gittel, A.; Mikutta, R.; Richter, A. Microbial community composition shapes enzyme patterns in topsoil and subsoil horizons along a latitudinal transect in Western Siberia. *Soil Biol. Biochem.* **2015**, *83*, 106–115. [CrossRef]

35. Kim, H.M.; Jung, J.Y.; Yergeau, E.; Hwang, C.; Hinzman, L.; Nam, S.; Hong, S.G.; Kim, O.-S.; Chun, J.; Lee, Y.K. Bacterial community structure and soil properties of a subarctic tundra soil in Council, Alaska. *FEMS Microbiol. Ecol.* **2014**, *89*, 465–475. [CrossRef]

36. Tveit, A.; Schwacke, R.; Svenning, M.M.; Urich, T. Organic carbon transformations in high-Arctic peat soils: Key functions and microorganisms. *ISME J.* **2013**, *7*, 299–311. [CrossRef] [PubMed]

37. Tveit, A.T.; Urich, T.; Frenzel, P.; Svenning, M.M. Metabolic and trophic interactions modulate methane production by Arctic peat microbiota in response to warming. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, E2507–E2516. [CrossRef] [PubMed]

38. Martineau, C.; Pan, Y.; Bodrossy, L.; Yergeau, E.; Whyte, L.G.; Greer, C.W. Atmospheric methane oxidizers are present and active in Canadian high Arctic soils. *FEMS Microbiol. Ecol.* **2014**, *89*, 257–269. [CrossRef] [PubMed]

39. Lee, S.-H.; Jang, I.; Chae, N.; Choi, T.; Kang, H. Organic Layer Serves as a Hotspot of Microbial Activity and Abundance in Arctic Tundra Soils. *Microb. Ecol.* **2012**, *65*, 405–414. [CrossRef] [PubMed]

40. Taş, N.; Prestat, E.; Wang, S.; Wu, Y.; Ulrich, C.; Kneafsey, T.; Tringe, S.G.; Torn, M.; Hubbard, S.S.; Jansson, J.K. Landscape topography structures the soil microbiome in arctic polygonal tundra. *Nat. Commun.* **2018**, *9*, 777. [CrossRef] [PubMed]

41. Vishnivetskaya, T.A.; Buongiorno, J.; Bird, J.; Krivushin, K.; Spirina, E.V.; Oshurkova, V.; Shcherbakova, V.A.; Wilson, G.; Lloyd, K.G.; Rivkina, E.M. Methanogens in the Antarctic Dry Valley permafrost. *FEMS Microbiol. Ecol.* **2018**, *94*, 109. [CrossRef] [PubMed]

42. Trotsenko, Y.A.; Khmelenina, V.N. Aerobic methanotrophic bacteria of cold ecosystems. *FEMS Microbiol. Ecol.* **2005**, *53*, 15–26. [CrossRef] [PubMed]

43. Shcherbakova, V.; Yoshimura, Y.; Ryzhmanova, Y.; Taguchi, Y.; Segawa, T.; Oshurkova, V.; Rivkina, E. Archaeal communities of Arctic methane-containing permafrost. *FEMS Microbiol. Ecol.* **2016**, *92*, fiw135. [CrossRef]

44. Comte, J.; Monier, A.; Crevecoeur, S.; Lovejoy, C.; Vincent, W.F. Microbial biogeography of permafrost thaw ponds across the changing northern landscape. *Ecography* **2015**, *39*, 609–618. [CrossRef]

45. Roiha, T.; Peura, S.; Cusson, M.; Rautio, M. Allochthonous carbon is a major regulator to bacterial growth and community composition in subarctic freshwaters. *Sci. Rep.* **2016**, *6*, 34456. [CrossRef]

46. Peura, S.; Wauthy, M.; Simone, D.; Eiler, A.; Einarsdóttir, K.; Rautio, M.; Bertilsson, S. Ontogenic succession of thermokarst thaw ponds is linked to dissolved organic matter quality and microbial degradation potential. *Linnol. Oceanogr.* **2019**, *65*, S248–S263. [CrossRef]

47. Przytulska, A.; Comte, J.; Crevecoeur, S.; Lovejoy, C.; Laurion, I.; Vincent, W. Phototrophic pigment diversity and picophyto-plankton in permafrost thaw lakes. *Biogeosciences* **2016**, *13*, 13–26. [CrossRef]

48. Steven, B.; Léveillé, R.; Pollard, W.H.; Whyte, L.G. Microbial ecology and biodiversity in permafrost. *Extremophiles* **2006**, *10*, 259–267. [CrossRef]

49. Liebner, S.; Harder, J.; Wagner, D. Bacterial diversity and community structure in polygonal tundra soils from Samoylov Island, Lena Delta, Siberia. *Int. J. Microbiol.* **2008**, *1*, 195–202.

50. Barbier, B.A.; Dziduch, I.; Liebner, S.; Ganzert, L.; Lantuit, H.; Pollard, W.; Wagner, D. Methane-cycling communities in a permafrost-affected soil on Herschel Island, Western Canadian Arctic: Active layer profiling ofmcrAandpmoAgenes. *FEMS Microbiol. Ecol.* **2012**, *82*, 287–302. [CrossRef]

51. Ganzert, L.; Bajerski, F.; Wagner, D. Bacterial community composition and diversity of five different permafrost-affected soils of Northeast Greenland. *FEMS Microbiol. Ecol.* **2014**, *89*, 426–441. [CrossRef]

52. Zhou, L.; Zhou, Y.; Yao, X.; Cai, J.; Liu, X.; Tang, X.; Zhang, Y.; Jang, K.S.; Jeppesen, E. Decreasing diversity of rare bacterial subcommunities relates to dissolved organic matter along permafrost thawing gradients. *Environ. Int.* **2020**, *134*, 105330. [CrossRef] [PubMed]

53. Coolen, M.J.L.; Van De Giessen, J.; Zhu, E.Y.; Wuchter, C. Bioavailability of soil organic matter and microbial community dynamics upon permafrost thaw. *Environ. Microbiol.* **2011**, *13*, 2299–2314. [CrossRef] [PubMed]
54. Reimer, P.J.; Austin, W.E.N.; Bard, E.; Bayliss, A.; Blackwell, P.G.; Ramsey, C.B.; Butzin, M.; Cheng, H.; Edwards, R.L.; Friedrich, M.; et al. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). *Radiocarbon* **2020**, *62*, 1–33. [CrossRef]

55. Stepanova, V.M.; Pokrovsky, O.S.; Viers, J.; Mironycheva-Tokareva, N.P.; Kosykh, N.P.; Vishnyakova, E.K. Major and trace elements in peat profiles in Western Siberia: Impact of the landscape context, latitude and permafrost coverage. *Appl. Geochem.* **2015**, *53*, 53–70. [CrossRef]

56. Koroleff, F. Determination of phosphorus. In *Methods for Seawater Analysis*; Grasshoff, K., Ehrhardt, M., Kremling, K., Eds.; Wiley-VCH Verlag GmbH: Weinheim, Germany, 1983; pp. 125–136.

57. Raudina, T.V.; Loiko, S.V.; Lim, A.; Manasypov, R.M.; Shirokova, L.S.; Istigechev, G.I.; Kuzmina, D.M.; Kulizhsky, S.P.; Vorobeyev, S.N.; Pokrovsky, O.S. Permafrost thaw and climate warming may decrease the CO₂, carbon, and metal concentration in peat soil waters of the Western Siberia Lowland. *Sci. Total. Environ.* **2018**, *634*, 1004–1023. [CrossRef]

58. Lim, A.G.; Loiko, S.V.; Kuzmina, D.M.; Krickov, I.V.; Shirokova, L.S.; Kulizhsky, S.P.; Vorobeyev, S.N.; Pokrovsky, O.S. Dispersed ground ice of permafrost peatlands: Potential unaccounted carbon, nutrient and metal sources. *Chemosphere* **2020**, *266*, 128953. [CrossRef]

59. Massiot, D.; Fayon, F.; Capron, M.; King, I.; Le Calvé, S.; Alonso, B.; Durand, J.O.; Bujoli, B.; Gan, Z.; Hoatson, G. Modelling one- and two-dimensional solid-state NMR spectra. *Magn. Reson. Chem.* **2002**, *40*, 70–76. [CrossRef]

60. Kögel-Knabner, I. 13C and 15N NMR spectroscopy as a tool in soil organic matter studies. *Geoderma* **1997**, *80*, 243–270. [CrossRef]

61. Dymov, A.; Zhangurov, E.; Hagedorn, F. Soil organic matter composition along altitudinal gradients in permafrost affected soils of the Subpolar Ural Mountains. *Catena* **2015**, *131*, 140–148. [CrossRef]

62. Gundelwein, A.; Müller-Lupp, T.; Sommerkorn, M.; Haupt, E.T.; Pfeiffer, E.M.; Wiechmann, H. Carbon in tundra soils in the Lake Labaz region of arctic Siberia. *Eur. J. Soil Sci.* **2007**, *58*, 1164–1174. [CrossRef]

63. Maksimova, V.; Klavina, L.; Bikovens, O.; Zicmanis, A.; Purmalis, O. Structural Characterization and Chemical Classification of Some Bryophytes Found in Latvia. *Chem. Biodivers.* **2013**, *10*, 1284–1294. [CrossRef] [PubMed]

64. Amundson, R.; Austin, A.T.; Schuur, E.A.G.; Yoo, K.; Matzek, V.; Kendall, C.; Uebersax, A.; Brenner, D.; Baisden, W.T. Global patterns of the isotopic composition of soil and plant nitrogen. *Glob. Biogeochem. Cycles* **2003**, *17*. [CrossRef]

65. Craine, J.M.; Brookshire, E.N.J.; Cramer, M.D.; Hasselquist, N.J.; Koba, K.; Marin-Spiotta, E.; Wang, L. Ecological interpretations of nitrogen isotope ratios of terrestrial plants and soils. *Plant Soil* **2015**, *396*, 1–26. [CrossRef]

66. Parinkina, O.M. *Microflora of Tundra Soils: Ecologico-Geographical Features and Productivity*; Nauka: St. Petersburg Leningrad, Russia, 1989; 159p.

67. Cowan, D.A.; Russell, N.J.; Mamais, A.; Sheppard, D.M. Antarctic Dry Valley mineral soils contain unexpectedly high levels of microbial biomass. *Extremophiles* **2002**, *6*, 431–436. [CrossRef]

68. Altshuler, I.; Goordial, J.; Whyte, L.G. Microbial life in permafrost. In *Psychrophiles: From Biodiversity to Biotechnology*; Margesin, R., Schinner, F., Marx, J.C., Gerday, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 153–179.

69. Jansson, J.; Taş, N. The microbial ecology of permafrost. *Nat. Rev. Microbiol.* **2014**, *12*, 414–425. [CrossRef]

70. Mackelprang, R.; Burkert, A.; Haw, M.; Mahendrarajah, T.; Conaway, C.H.; Douglas, T.A.; Waldrop, M.P. Microbial survival strategies in ancient permafrost: Insights from metagenomics. *ISME J.* **2017**, *11*, 2305–2318. [CrossRef]

71. Frank-Fahle, B.A.; Yergeau, É.; Greer, C.W.; Lantuit, H.; Wagner, D. Microbial Functional Potential and Community Composition in Permafrost-Affected Soils of the NW Canadian Arctic. *PLoS ONE* **2014**, *9*, e84761. [CrossRef]

72. Tuorto, S.J.; Darias, P.; McGuinness, L.R.; Panikov, N.; Zhang, T.; Haggblom, M.; Kerkhof, L.J. Bacterial genome replication at subzero temperatures in permafrost. *ISME J.* **2013**, *8*, 139–149. [CrossRef]

73. Huhe, Chen, X.; Hou, F.; Wu, Y.; Cheng, Y. Bacterial and Fungal Community Structures in Loess Plateau Grasslands with Different Grazing Intensities. *Front. Microbiol.* **2017**, *8*, 606. [CrossRef] [PubMed]

74. Carvalho, F.M.; Souza, R.C.; Barcellos, F.G.; Hungria, M.; Vasconcelos, A.T.R. Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales. *BMC Microbiol.* **2010**, *10*, 37. [CrossRef] [PubMed]

75. Erlacher, A.; Cernava, T.; Cardinale, M.; Soh, J.; Sensen, C.W.; Grube, M.; Berg, G. Rhizobiales as functional and endosymbiontial members in the lichen symbiosis of *Lobaria pulmonaria*. *Front. Microbiol.* **2015**, *6*, 53. [CrossRef] [PubMed]

76. Garrido-Oter, R.; Nakano, R.T.; Dombrowski, N.; Ma, K.-W.; Tea, T.A.; McHardy, A.C.; Schulze-Lefert, P. Modular Traits of the Rhizobiales Root Microbiota and Their Evolutionary Relationship with Symbiotic Rhizobia. *Cell Host Microbe* **2018**, *24*, 155–167.e5. [CrossRef]

77. Cole, J.K.; Morton, B.R.; Cardamone, H.C.; Lake, H.R.R.; Dohnalkova, A.C.; Kim, Y.-M.; Romine, M.F.; Kyle, J.E.; Maezato, Y.; Dana, K.L.; et al. *Salinivirga fredricksonii* gen. nov., sp. nov., a heterotrophic halophile isolated from a photosynthetic mat, a member of a novel lineage (*Salinarimonadaceae* fam. nov.) within the order *Rhizobiales*, and reclassification of the genus *Salinarimonas* Liu et al. 2010 into *Salinarimonadaceae*. *Int. J. Syst. Evol. Microbiol.* **2018**, *68*, 1591–1598. [CrossRef]

78. Taylor, B.N.; Simms, E.L.; Komatsu, K.J. More Than a Functional Group: Diversity within the Legume-Rhizobia Mutualism and Its Relationship with Ecosystem Function. *Divers.* **2020**, *12*, 50. [CrossRef]

79. Wang, S.; Meade, A.; Lam, H.-M.; Luo, H. Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales. *mSystems* **2020**, *5*, e00438-20. [CrossRef] [PubMed]

80. Kielak, A.M.; Barreto, C.C.; Kowalchuk, G.A.; Van Veen, J.A.; Kuramae, E.E. The Ecology of Acidobacteria: Moving beyond Genes and Genomes. *Front. Microbiol.* **2016**, *7*, 744. [CrossRef]
54. Reimer, P.J.; Austin, W.E.N.; Bard, E.; Bayliss, A.; Blackwell, P.G.; Ramsey, C.B.; Butzin, M.; Cheng, H.; Edwards, R.L.; Friedrich, M.; et al. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). *Radiocarbon* **2020**, *62*, 1–33. [CrossRef]

55. Stepanova, V.M.; Pokrovsky, O.S.; Viers, J.; Mironycheva-Tokareva, N.P.; Kosykh, N.P.; Vishnyakova, E.K. Major and trace elements in peat profiles in Western Siberia: Impact of the landscape context, latitude and permafrost coverage. *Appl. Geochem.* **2015**, *53*, 53–70. [CrossRef]

56. Koroleff, F. Determination of phosphorus. In *Methods for Seawater Analysis*; Grasshoff, K., Ehrhardt, M., Kremling, K., Eds.; Wiley-VCH Verlag GmbH: Weinheim, Germany, 1983; pp. 125–136.

57. Raudina, T.V.; Loiko, S.V.; Lim, A.; Manasypov, R.M.; Shirokova, L.S.; Istigechev, G.I.; Kuzmina, D.M.; Kulizhsky, S.P.; Vorobeyev, S.N.; Pokrovsky, O.S. Permafrost thaw and climate warming may decrease the CO₂, carbon, and metal concentration in peat soil waters of the Western Siberia Lowland. *Sci. Total. Environ.* **2018**, *634*, 1004–1023. [CrossRef]

58. Lim, A.G.; Loiko, S.V.; Kuzmina, D.M.; Krickov, I.V.; Shirokova, L.S.; Kulizhsky, S.P.; Vorobeyev, S.N.; Pokrovsky, O.S. Dispersed ground ice of permafrost peatlands: Potential unaccounted carbon, nutrient and metal sources. *Chemosphere* **2020**, *266*, 128953. [CrossRef]

59. Massiot, D.; Fayon, F.; Capron, M.; King, I.; Le Calvé, S.; Alonso, B.; Durand, J.O.; Bujoli, B.; Gan, Z.; Hoatson, G. Modelling one- and two-dimensional solid-state NMR spectra. *Magn. Reson. Chem.* **2002**, *40*, 70–76. [CrossRef]

60. Kögel-Knabner, I. 13C and 15N NMR spectroscopy as a tool in soil organic matter studies. *Geoderma* **1997**, *80*, 243–270. [CrossRef]

61. Dymov, A.; Zhangurov, E.; Hagedorn, F. Soil organic matter composition along altitudinal gradients in permafrost affected soils of the Subpolar Ural Mountains. *Catena* **2015**, *131*, 140–148. [CrossRef]

62. Gundelwein, A.; Müller-Lupp, T.; Sommerkorn, M.; Haupt, E.T.; Pfeiffer, E.M.; Wiechmann, H. Carbon in tundra soils in the Lake Labaz region of arctic Siberia. *Eur. J. Soil Sci.* **2007**, *58*, 1164–1174. [CrossRef]

63. Maksimova, V.; Klavina, L.; Bikovens, O.; Zicmanis, A.; Purmalis, O. Structural Characterization and Chemical Classification of Some Bryophytes Found in Latvia. *Chem. Biodivers.* **2013**, *10*, 1284–1294. [CrossRef] [PubMed]

64. Amundson, R.; Austin, A.T.; Schuur, E.A.G.; Yoo, K.; Matzek, V.; Kendall, C.; Uebersax, A.; Brenner, D.; Baisden, W.T. Global patterns of the isotopic composition of soil and plant nitrogen. *Glob. Biogeochem. Cycles* **2003**, *17*. [CrossRef]

65. Craine, J.M.; Brookshire, E.N.J.; Cramer, M.D.; Hasselquist, N.J.; Koba, K.; Marin-Spiotta, E.; Wang, L. Ecological interpretations of nitrogen isotope ratios of terrestrial plants and soils. *Plant Soil* **2015**, *396*, 1–26. [CrossRef]

66. Parinkina, O.M. *Microflora of Tundra Soils: Ecologico-Geographical Features and Productivity*; Nauka: St. Petersburg Leningrad, Russia, 1989; 159p.

67. Cowan, D.A.; Russell, N.J.; Mamais, A.; Sheppard, D.M. Antarctic Dry Valley mineral soils contain unexpectedly high levels of microbial biomass. *Extremophiles* **2002**, *6*, 431–436. [CrossRef]

68. Altshuler, I.; Goordial, J.; Whyte, L.G. Microbial life in permafrost. In *Psychrophiles: From Biodiversity to Biotechnology*; Margesin, R., Schinner, F., Marx, J.C., Gerday, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 153–179.

69. Jansson, J.; Taş, N. The microbial ecology of permafrost. *Nat. Rev. Microbiol.* **2014**, *12*, 414–425. [CrossRef]

70. Mackelprang, R.; Burkert, A.; Haw, M.; Mahendrarajah, T.; Conaway, C.H.; Douglas, T.A.; Waldrop, M.P. Microbial survival strategies in ancient permafrost: Insights from metagenomics. *ISME J.* **2017**, *11*, 2305–2318. [CrossRef]

71. Frank-Fahle, B.A.; Yergeau, É.; Greer, C.W.; Lantuit, H.; Wagner, D. Microbial Functional Potential and Community Composition in Permafrost-Affected Soils of the NW Canadian Arctic. *PLoS ONE* **2014**, *9*, e84761. [CrossRef]

72. Tuorto, S.J.; Darias, P.; McGuinness, L.R.; Panikov, N.; Zhang, T.; Haggblom, M.; Kerkhof, L.J. Bacterial genome replication at subzero temperatures in permafrost. *ISME J.* **2013**, *8*, 139–149. [CrossRef]

73. Huhe, Chen, X.; Hou, F.; Wu, Y.; Cheng, Y. Bacterial and Fungal Community Structures in Loess Plateau Grasslands with Different Grazing Intensities. *Front. Microbiol.* **2017**, *8*, 606. [CrossRef] [PubMed]

74. Carvalho, F.M.; Souza, R.C.; Barcellos, F.G.; Hungria, M.; Vasconcelos, A.T.R. Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales. *BMC Microbiol.* **2010**, *10*, 37. [CrossRef] [PubMed]

75. Erlacher, A.; Cernava, T.; Cardinale, M.; Soh, J.; Sensen, C.W.; Grube, M.; Berg, G. Rhizobiales as functional and endosymbiontial members in the lichen symbiosis of *Lobaria pulmonaria*. *Front. Microbiol.* **2015**, *6*, 53. [CrossRef] [PubMed]

76. Garrido-Oter, R.; Nakano, R.T.; Dombrowski, N.; Ma, K.-W.; Tea, T.A.; McHardy, A.C.; Schulze-Lefert, P. Modular Traits of the Rhizobiales Root Microbiota and Their Evolutionary Relationship with Symbiotic Rhizobia. *Cell Host Microbe* **2018**, *24*, 155–167.e5. [CrossRef]

77. Cole, J.K.; Morton, B.R.; Cardamone, H.C.; Lake, H.R.R.; Dohnalkova, A.C.; Kim, Y.-M.; Romine, M.F.; Kyle, J.E.; Maezato, Y.; Dana, K.L.; et al. *Salinivirga fredricksonii* gen. nov., sp. nov., a heterotrophic halophile isolated from a photosynthetic mat, a member of a novel lineage (*Salinarimonadaceae* fam. nov.) within the order *Rhizobiales*, and reclassification of the genus *Salinarimonas* Liu et al. 2010 into *Salinarimonadaceae*. *Int. J. Syst. Evol. Microbiol.* **2018**, *68*, 1591–1598. [CrossRef]

78. Taylor, B.N.; Simms, E.L.; Komatsu, K.J. More Than a Functional Group: Diversity within the Legume-Rhizobia Mutualism and Its Relationship with Ecosystem Function. *Divers.* **2020**, *12*, 50. [CrossRef]

79. Wang, S.; Meade, A.; Lam, H.-M.; Luo, H. Evolutionary Timeline and Genomic Plasticity Underlying the Lifestyle Diversity in Rhizobiales. *mSystems* **2020**, *5*, e00438-20. [CrossRef] [PubMed]

80. Kielak, A.M.; Barreto, C.C.; Kowalchuk, G.A.; Van Veen, J.A.; Kuramae, E.E. The Ecology of Acidobacteria: Moving beyond Genes and Genomes. *Front. Microbiol.* **2016**, *7*, 744. [CrossRef]
81. Ventura, M.; Canchaya, C.; Tauch, A.; Chandra, G.; Fitzgerald, G.F.; Chater, K.F.; van Sinderen, D. Genomics of Actinobacteria: Tracing the Evolutionary History of an Ancient Phylum. *Microbiol. Mol. Biol. Rev.* **2007**, *71*, 495–548. [CrossRef] [PubMed]

82. Lladó, S.; Žifčáková, L.; Větrovský, T.; Eichlerová, I.; Baldrian, P. Functional screening of abundant bacteria from acidic forest soil indicates the metabolic potential of Acidobacteria subdivision 1 for polysaccharide decomposition. *Biol. Fertil. Soils* **2016**, *52*, 251–260. [CrossRef]
81. Ventura, M.; Canchaya, C.; Tauch, A.; Chandra, G.; Fitzgerald, G.F.; Chater, K.F.; van Sinderen, D. Genomics of Actinobacteria: Tracing the Evolutionary History of an Ancient Phylum. *Microbiol. Mol. Biol. Rev.* **2007**, *71*, 495–548. [CrossRef] [PubMed]

82. Lladó, S.; Žifčáková, L.; Větrovský, T.; Eichlerová, I.; Baldrian, P. Functional screening of abundant bacteria from acidic forest soil indicates the metabolic potential of Acidobacteria subdivision 1 for polysaccharide decomposition. *Biol. Fertil. Soils* **2016**, *52*, 251–260. [CrossRef]