Beat Frey
Gilda Varliero
Weihong Qi3
Beat Stierli
Lorenz Walthert
Ivano Brunner- 1Forest Soils and Biogeochemistry, Swiss Federal Institute for Forest, Snow and Landscape Research (WSL), Birmensdorf, Switzerland
- 2Centre for Microbial Ecology and Genomics, Genetics and Microbiology, University of Pretoria, Pretoria, South Africa
- 3Functional Genomics Center Zurich (FGCZ), ETH Zürich/University of Zurich, Zurich, Switzerland
Soil microorganisms such as Bacteria and Archaea play important roles in the biogeochemical cycling of soil nutrients, because they act as decomposers or are mutualistic or antagonistic symbionts, thereby influencing plant growth and health. In the present study, we investigated the vertical distribution of soil metagenomes to a depth of 1.5 m in Swiss forests of European beech and oak species on calcareous bedrock. We explored the functional genetic potential of soil microorganisms with the aim to disentangle the effects of tree genus and soil depth on the genetic repertoire, and to gain insight into the microbial C and N cycling. The relative abundance of reads assigned to taxa at the domain level indicated a 5–10 times greater abundance of Archaea in the deep soil, while Bacteria showed no change with soil depth. In the deep soil there was an overrepresentation of genes for carbohydrate-active enzymes, which are involved in the catalyzation of the transfer of oligosaccharides, as well as in the binding of carbohydrates such as chitin or cellulose. In addition, N-cycling genes (NCyc) involved in the degradation and synthesis of N compounds, in nitrification and denitrification, and in nitrate reduction were overrepresented in the deep soil. Consequently, our results indicate that N-transformation in the deep soil is affected by soil depth and that N is used not only for assimilation but also for energy conservation, thus indicating conditions of low oxygen in the deep soil. Using shotgun metagenomics, our study provides initial findings on soil microorganisms and their functional genetic potential, and how this may change depending on soil properties, which shift with increasing soil depth. Thus, our data provide novel, deeper insight into the “dark matter” of the soil.
Introduction
In agricultural and forested areas, studies on soil microbial communities rarely include soil depths exceeding 30 cm (Yost and Hartemink, 2020). Most of the organic C in soils from various biomes is located in the topsoil at 0–30 cm soil depth (Crowther et al., 2019). Jobbágy and Jackson (2000) estimated in a worldwide study that 50% of the organic C in forests is in the top 20 cm. Xu et al. (2013) similarly calculated that about 70% of the soil microbial C and N in forested ecosystems is located in the upper 30 cm of the soils. In a few recent studies on soils, however, measurement went beyond the topsoil and included the deep soil down to a depth of 1 m in investigations of extracellular enzyme activities (Dove et al., 2020) or metagenomic attributes of novel bacterial taxa (Brewer et al., 2019). However, studies of the soil microbiome including soils at depths greater than 1 m are still particularly rare, leaving this soil compartment largely unexplored as a terra incognita (Andújar et al., 2017).
Deep soil layers harbor poorly known bacterial and fungal communities (Brewer et al., 2019; Robin et al., 2019; Liu et al., 2020; Frey et al., 2021). It is assumed that living conditions in deeper soils for soil organisms and plant roots are harsher in deeper soils, as a result of higher soil density, lower oxygen concentrations, and lower carbon (C) and nutrient availability (Lennon, 2020). Because soil processes, soil properties, and microbial communities are depth-dependent, soils should be studied at greater depths for a more comprehensive understanding of their relationship and interaction (Yost and Hartemink, 2020). In addition, it is important to understand deep soil microbial communities because deep soils are poorly accounted for in current models of biogeochemical processes.
Microorganisms are the key drivers of both C- and nutrient-cycling processes in the soil (Falkowski et al., 2008). Most soil microorganisms in forest ecosystems gain their energy by decomposing different types of C from soil organic C (SOC), which vary in degradability and require different mechanisms for decomposition (Lladó et al., 2019). Carbon and nutrients typically are more available in the topsoil, mainly because of larger input of leaf and fine root litter and root exudates, as well as higher biotic activity (Eldridge et al., 2016; Liu et al., 2020). Consequently, microbial biomass per unit soil mass is usually one to two orders of magnitude lower in the deep soil than in the topsoil (Eilers et al., 2012; Spohn et al., 2016). Nevertheless, the total biomass of Bacteria and fungi inhabiting deeper soil horizons can be as large as in the topsoil, and these microorganisms play similarly important roles in biogeochemical processes (Brewer et al., 2019), although we know very little about the microbial genetic potential including soil C and nutrient cycling, soil formation, iron redox reactions, and pollutant degradation.
Because C usually occurs in lower amounts in deeper soil layers compared to the topsoil, chemolithotrophic processes could become prominent in deep soil layers. This means that microorganisms could gain their energy from the oxidation of inorganic compounds instead of organic compounds (Kelly and Wood, 2013). The most common substrates that are reduced are nitrogen (N), sulfur and iron compounds as well as hydrogen (Thiel, 2011). Besides of inorganic C, oxygen concentrations also tend to be lower in deeper soil layers. Thus, anaerobic processes can also occur, where energy is generated through the reduction of nitrate, nitrite or sulfate (Hooper and Dispirito, 2013). Some microorganisms are also able to use CO2 as a C source. In this case, the resulting products are methane and acetic acid (Hooper and Dispirito, 2013). Ammonia-oxidizing Bacteria (AOB) and ammonia-oxidizing Archaea (AOA) are central players in the global N cycle. These microbial groups can perform the first part of the nitrification process, ammonia oxidation. This is an important rate-limiting step in the various N-cycling processes, including N fixation, mineralization, nitrification, and denitrification (Kowalchuk and Stephen, 2001; Zhang et al., 2013). Since AOB gain only little energy by oxidizing ammonia, a lower soil NH4+-N content in the deep soil leads to a decrease in AOB abundance, whereas AOA appear to be less constrained by the availability of N substrates and prosper in oligotrophic environments in the deep soil (Martens-Habbena et al., 2009; Trivedi et al., 2019).
In deeper soil layers, weathering processes can additionally contribute to the presence of essential nutrient elements such as Mg, P and Fe, through the mineralization of the bedrock (Uroz et al., 2009, 2015; Xi et al., 2018). Bacteria and fungi play a role in rock weathering by producing organic acids, hydrogen cyanide or siderophores (Frey et al., 2010; Brunner et al., 2011; Wongfun et al., 2014; Xi et al., 2018; Wang et al., 2019). In addition, specific environmental enzymes involved in S and Fe metabolisms could be relevant (Jiang et al., 2019; Li et al., 2020). Genes related to such weathering processes could be those involved in oxalate biosynthesis, in cyanide synthesis, and in siderophore synthesis and transport (Varliero et al., 2021).
In a first survey of the deep soil horizons in Swiss beech and oak forests, we observed the occurrence of poorly known bacterial taxa belonging to Nitrospirae, Chloroflexi, Rokubacteria, Gemmatimonadetes, and Firmicutes (Frey et al., 2021). Furthermore, archaeal phyla such as Thaumarchaeota and Euryarchaeota were more abundant in the deep soil than in the topsoil. Previous research has indicated that members of Chloroflexi might be adapted to nutrient-poor environmental conditions (Rime et al., 2016; Adamczyk et al., 2019) whereas Gemmatimonadetes might be specifically adapted to low-moisture availability (DeBruyn et al., 2011; Rime et al., 2015). It is noteworthy, that some members of Euryarchaeota produce methane (Bodelier and Steenbergh, 2014), but also oxidize methane, fix N, and reduce nitrates (Cabello et al., 2004).
The aim of the present study was to improve our understanding of the metabolic capabilities of the microbial communities in deep soil layers. While we focused on changes in taxonomic composition with increasing soil depth in our first survey (Frey et al., 2021), here we used shotgun metagenomics to assess changes in microbial C- and N-cycling potential and link these data to the properties in deep soil layers. In particular, we hypothesized that in deep soil layers: (1) the genus of the dominant tree species in a site does not play a crucial role in shaping microbial metabolic capabilities in soil, (2) the biogeochemical and chemolithotrophic capabilities of the microbial communities are enhanced relative to shallower soil, and (3) soil N cycling is particularly affected by the metabolic capabilities of the microbial communities. We specifically focused on genes in the Carbohydrate-Active enZymes (CAZy) database, which are responsible for the catabolism of various C-complexes with varying decomposability, ranging from labile C (e.g., monosaccharides and polysaccharides) to recalcitrant C (e.g., lignin). Our study is one of the first to characterize the functional genetic potential and metabolic capabilities of microorganisms in deep soil layers and contributes to an integrated perspective on the soil microbial communities of beech and oak forests.
Materials and Methods
Forest Sites
The six investigated forest sites are located in the relatively warm and dry Rhone Valley (Chamoson and Saillon) and in the Jura Mountains (Neunkirch) of Switzerland with mean annual precipitation sums of 815–955 mm and mean annual temperatures of 8.5–9.4°C (Supplementary Table 1; Remund et al., 2014). The sites are covered mainly by European beech (Fagus sylvatica L.) or by oak species (Quercus spp.), with oak sites having warmer and drier soils than European beech sites (see also Frey et al., 2021). All of the forest stands have a near natural tree species composition originating from natural regeneration and have been unmanaged for several decades (Frey et al., 2021). The forest sites are part of the research platforms WSL Forest Soil Database. The soil classification (Supplementary Table 1) used in this study followed the system of the Food and Agriculture Organization of the United Nations (IUSS Working Group WRB, 2007).
Soil Sampling and Soil Chemical and Physical Analyses
A soil profile down to 2 m depth was excavated in the center of a group of dominant trees at each study site in spring 2014. Three soil samples per soil depth were taken with 1 L steel cylinders (10.7 cm height and 11.1 cm diameter) at depths of 0–10 cm (L1), 15–25 cm (L2), 45–55 cm (L3), 75–85 cm (L4), 110–125 cm (L5), 140–155 cm (L6), and 180–200 cm (L7) (for details see Frey et al., 2021). For soil chemical and physical analyses, approximately 3–4 kg of soil was sampled separately from the various depths (L1–L7). All soil samples were then transported in plastic bags to the laboratories and stored at 4°C until further processing. The soil samples for chemical and particle size analyses were dried a 60°C and sieved through a 2 mm mesh sieve. Soil pH was determined in 0.01 M CaCl2. For C and N analyses, an aliquot of the dried and sieved samples was ground for 3 min using a vibrating ball mill (MM2000, Retsch, Haan, Germany) with zircon-grinding tools. Organic C (Corg) and total N (Ntot) were determined with a CN elemental analyzer (NC2500; CE Instruments Ltd., Hindley Green, United Kingdom) whereas inorganic carbonates were removed by HCl fumigation of the samples prior to combustion (Walthert et al., 2010). Exchangeable cations (Na, K, Mg, Ca, Mn, Al, and Fe) were extracted (in triplicate) in an unbuffered solution of 1 M NH4Cl. The effective cation-exchange capacity (CEC) was obtained by summing the charge equivalents of exchangeable Na, K, Mg, Ca, Mn, Al, Fe, and H. Base saturation (BS) was calculated as the ratio of the sum of exchangeable Na, K, Mg, and Ca to the CEC. Concentrations of the elements in the extracts were determined by ICP-OES (Optima 3000, PerkinElmer Inc., Waltham, MA, United States). Soil particles were fractionated into sand, silt and clay using the pipette method according to Gee and Bauder (1986). Plant-available water storage capacity (AWC), i.e., the storable amount of water between field capacity and permanent wilting point, was calculated with a pedotransfer function according to Teepe et al. (2003) for each sampled soil layer for a thickness of 10 cm considering fine-earth density and soil texture, i.e., particle size distribution. Further details of the soil chemical and soil physical analysis are described in Frey et al. (2021).
In order to estimate fine-root biomass (<2 mm in diameter) and fine-earth densities, soil samples from the 1 L steel cylinders were sieved with a 4 mm sieve and the fine roots were sorted out (see also Frey et al., 2021). For soil microbial DNA sampling, fresh soil aliquots of about 5 g were put into a Ziploc® bag and stored frozen at −80°C until DNA analyses. Finally, the remaining soil was dried at 105°C, and volumetric stone content and fine-earth density were determined.
To distinguish between topsoil and deep soil, data from soil layers L1 and L2 (topsoil; 0 – 25 cm) and from L5 and L6 (deepsoil; 110–155 cm) were summed and averaged. Soil measurements represent the mean for the topsoil and the deep soil, respectively. Data from L3, L4, and L7 were not included in this study.
DNA Extraction and Shotgun Sequencing
DNA was extracted using the PowerSoil DNA Isolation Kit (Qiagen, Hilden, Germany) and was quantified using the high-sensitivity Qubit assay (Thermo Fisher Scientific, Reinach, Switzerland). Library preparation using the TruSeq DNA Library Prep Kit (Illumina Inc., San Diego, CA, United States) and shotgun sequencing of the eluted DNA samples was performed at Microsynth AG (Balgach, Switzerland), using the Illumina NovaSeq system (2 × 150 bp; Illumina Inc., San Diego, CA, United States). The twelve metagenomes were from three different forest sites (Chamoson, Neunkirch, and Saillon), each with two different forest types (beech forest and oak forest) and two different soil depths (topsoil and deepsoil). Extracted DNA was pooled from soil layers L1 and L2 for topsoil samples, and from layers L5 and L6 for deepsoil samples. Raw sequences were uploaded in the NCBI Sequence Read Archive under the accession number PRJNA783873.
Assembly and Functional Annotation of Assembled Contigs
Pre-processing of metagenomic reads, assembly of reads into contigs, contig binning, and functional and phylogenetic annotation of contigs and bins were achieved using a customized pipeline (Donhauser et al., 2021; Perez-Mon et al., 2021). Briefly, raw reads were quality checked using FastQC v.0.11.91. They were quality filtered and trimmed (i.e., pre-processed reads) using Trimmomatic v0.36 (Q = 20, minimum read length = 40; Bolger et al., 2014). Pre-processed read pairs and singleton reads from all samples were assembled into contigs (>200 bp) by iteratively building de Bruijn graphs using k-mers of increasing size with the de novo assembler MEGAHIT v1.1.3 (–k-min 27 –k-step 10; Li et al., 2015).
Genes were annotated to the coassembly (i.e., contigs obtained from the assembly of all the samples together) as this approach considerably increases the quality of the coding sequences (CDS) and gene annotation (e.g., Tamames et al., 2019). CDSs were predicted with MetaGeneMark v3.38 (Zhu et al., 2010). To uncover the potential metabolic capabilities of the soil metagenomes, CDS and genes were assigned to functions (i.e., functional genes). About 45% of the predicted genes (22,941,682) were assigned to general metabolic and cellular functions through EggNOG v4.5 (Evolutionary genealogy of genes: Non-supervised Orthologous Groups), which classifies the genes into clusters of orthologous groups (COGs) of proteins and organizes the COGs into general functional categories (Jensen et al., 2008; Huerta-Cepas et al., 2016). Annotation to EggNOG v4.5 was performed using the EggNOG-mapper v1.0.3 with the DIAMOND v.2.0.7 search mode against all protein sequences (Huerta-Cepas et al., 2017). About 1% of the protein-coding genes (507,956) were annotated to carbohydrate-active enzymes using the CAZy database (Carbohydrate-Active enZymes: release of July 2017 version; Cantarel et al., 2009). About 0.2% of the protein-coding genes (82,022) were annotated to N-cycling families using the N-cycling genes (NCyc) database (syn. NCycDB: curated integrative database for fast and accurate metagenomic profiling of NCyc clustered at 100% sequence identity; Tu et al., 2019). Annotations against the CAZy and NCyc databases were performed using SWORD v1.0.3 (Vaser et al., 2016) (−v 10–5; Anwar et al., 2019).
Kaiju v1.7.4, a program for sensitive taxonomic classification of high-throughput sequencing reads from metagenomic whole genome sequencing (Menzel et al., 2016), was used for the taxonomic classification of the assembled contigs using the prebuilt “nr + euk” database containing the bacterial, archaeal, viral, fungal, and microbial eukaryotic protein sequences from the NCBI BLAST non-redundant protein database (version 2021-02-24).
Abundance Quantification of Protein-Coding Genes and Taxonomic Profiles
Pre-processed read pairs from each of the samples were mapped back to the assembled contigs, using the BWA aligner v0.7.15 (bwa-mem; Li, 2013). The function featureCounts from the package Subread v1.5.1 (-minOverlap 10, Q = 10, -primary, Liao et al., 2014) was used to count the reads mapped within the predicted protein coding. Taxonomy was also assigned by counting reads that mapped back to the taxonomic assigned contigs.
Quantitative PCR of Functional Genes
Abundances of C- and NCyc were quantified using primers and thermocycling conditions as reported in Frey et al. (2020). Functional marker genes encoding for enzymes catalyzing major processes during methane oxidation (pmoA), methanogenesis (mcrA), denitrification (nirS, nosZ), nitrification (bacterial amoA, archaeal amoA, nxrB), and N fixation (nifH), are targeted. The specificity of the amplification products was confirmed by melting-curve analysis, and the expected sizes of the amplified fragments were checked in a 1.5% agarose gel stained with ethidium bromide. Three standard curves per target region (correlations ≥ 0.997) were obtained using tenfold serial dilutions (10–1 to 10–9 copies) of plasmids generated from cloned targets (Frey et al., 2011). Data was converted to represent average copy number of targets per μg DNA.
Data Analyses
Statistical analyses were completed using the open-source software R v4.1.0 (R Core Team, 2019), and graphical representations of results were created with the R package ggplot2 v3.3.5 (Wickham, 2016), unless specified otherwise. Protein-coding genes for which the sum of the reads over all soil samples was <10 were excluded from the analyses to lower the false discovery rate caused by stochasticity between samples.
To evaluate the dissimilarities in metabolic capabilities between the different soil habitats, we identified the functional genes - annotated against EggNOG, CAZy and NCyc databases – that were differentially abundant between the soils using the R package DESeq2 v 1.34 (Love et al., 2014). Data were normalized using DESeq2 and corrected for sequencing depth and annotated gene number variation, following the standard protocol used by the DESeq2 pipeline. All analyses were performed on the normalized dataset. To calculate the differentially abundant genes for each contrast, effect size shrinkage was applied using the lfcShrink() function with the ‘apeglm’ method (Zhu et al., 2019). Genes were reported as significantly over- or underrepresented only for pairwise comparisons with an adjusted P-value lower than 0.01; P-values were adjusted for multiple testing using the Benjamini–Hochberg method. Bray–Curtis dissimilarities between samples were visualized with non-metric multidimensional scaling (NMDS) (vegan R package; Oksanen et al., 2019). The Bray–Curtis dissimilarity matrix was calculated on the vst transformed gene dataset [varianceStabilizingTransformation() function from the DESeq2 package]. The envfit() function was then used to calculate and plot the environmental vectors in the NMDS orientation plot.
The statistical significance of observed differences was assessed with permutational analyses of variance (PERMANOVA, 105 permutations, Monte Carlo approximated p-value). Multivariate homogeneity of group dispersions was checked prior to the PERMANOVAs to ensure that detected significant differences were associated with the tested factors and not with differences in the within-group variabilities. To calculate alpha diversity indices, the raw gene count dataset was resampled three times (n = 5,000,000) with the ‘rrarefy’ function from the R package ‘vegan.’ We calculated richness with the function ‘specnumber,’ and the Shannon index with the diversity() function.
Results
Forest Site Property Change With Depth
Analysis of soil chemical and physical parameters at the six forest sites showed that they were not dependent on the dominant tree genus at the site, as there were no significant differences in any of these variables. In contrast, some of the soil properties showed significant changes with soil depth (Table 1). Soil pH increased significantly with soil depth, whereas soil base saturation overall was high, around 100%. Both, Corg and Ntot decreased significantly with soil depth, Corg about tenfold and Ntot about fourfold. In contrast, soil density increased significantly with soil depth by about a factor of 1.5 (Table 1).
Table 1. Mean values of soil biological, chemical, and physical properties (n = 3 forest sites).
Analysis of fine-root variables indicated, that oak sites had no significantly greater fine-root biomass compared with beech sites (Table 1). However, the fine-root biomass decreased significantly with increasing soil depth in both tree species by about tenfold, similarly to the microbial biomass, which also decreased significantly by about tenfold (Table 1).
Overall Metagenome Sequencing Results
After quality filtering of the metagenomes of twelve soil samples we obtained in total 1.6 × 109 high-quality reads (ranging from 1.2 to 1.7 × 108 reads per soil sample). The assembly consisted of 3.7 × 107 MEGAHIT contigs of 568 bp on average, ranging from 202 bp to 4.1 × 105 bp, with a N50 value of 581 bp and a GC-content of 64% (Supplementary Table 2). We found a significantly higher percentage of reads mapped to contigs and to protein coding genes in the deep soil (contigs: 83.3%; CDS: 66.2%) than in the topsoil (contigs: 71.7%; CDS: 53.6%) (Table 2), whereas the percentages did not differ significantly between the tree genera. In total, 5.0 × 107 predicted genes were found among the contigs, 2.3 × 107 of which we could annotate with the EggNOG database, 5.1 × 105 with the C-cycling gene database CAZy, and 8.2 × 104 with the N-cycling gene database NCyc (Supplementary Table 2). The fraction of the predicted genes annotated by each database were 45.1% (for EggNOG), 1.0% (for CAZy), and 0.2% (for NCyc).
Table 2. Mean number of sequences and percentage of protein-coding genes (CDS genes), and the relative abundance of contigs assigned to taxa at the domain level for topsoil or deepsoil (n = 3 forest sites).
Taxonomic Composition of the Metagenomes
Contigs were classified to different taxa, whereas reads mapped to classified contigs were used to quantify the taxa abundance. Bacteria were always the most common group in all samples, with relative abundances of >90% in each sample, whereby soil depth had no significant influence on abundance. Archaea were 5–9 times more abundant in the deep soil than in the topsoil, whereas both Bacteria and Archaea were unaffected by the tree genus (Table 2). In contrast to Archaea, Eukarya were more abundant in the topsoil than in the deep soil (Table 2). Viruses were always the least abundant group and did not differ between soil depths or tree genera (Table 2).
The most abundant archaeal phyla, such as Crenarchaeota, Euryarchaeota, and Thaumarchaeota, tended to be significantly more abundant in the deep soil than in the topsoil (Table 3). In particular, the classes of the anaerobic Methanomicrobia and the chemoautotrophic nitrifying Nitrososphaeria were significantly more abundant in the deep soil. Among the Bacteria, the phyla Actinobacteria, Chloroflexi, and Nitrospirae were significantly more abundant in the deep soil, whereas Proteobacteria were more abundant in the topsoil (Table 3). Among them the bacterial classes such as Actinomycetia, Anaerolineae, Chloroflexia, and Nitrospira were significantly more abundant in the deep soil, whereas the classes Alphaproteobacteria and Gammaproteobacteria were significantly more abundant in the topsoil.
Table 3. Mean number of reads assigned to taxonomic contigs of selected phyla and classes of the domains Archaea and Bacteria (n = 3 forest sites).
Diversity of the Protein-Coding Genes Changes With Soil Depth
There was an overall significant effect of the soil depth on the alpha-diversity of functional genes calculated over the entire dataset (Pseudo-F = 17.7, R2 = 0.34, P < 0.001 for richness, Pseudo-F = 17.2, R2 = 0.34, P < 0.001 for Shannon index), as well as for the CAZy dataset (Pseudo-F = 6.33, R2 = 0.16, P = 0.016 for richness, Pseudo-F = 13.9, R2 = 0.29, P < 0.001 for Shannon index) and the NCyc dataset (Pseudo-F = 12.7, R2 = 0.27, P < 0.001 for Shannon index; Supplementary Table 3). The tree genus, in contrast, did not show any significant effect on the alpha-diversity for any of the three datasets. In the analysis of the entire dataset, however, site had a weak significant effect on Shannon index (Pseudo-F = 3.30, R2 = 0.17, P = 0.048) but not on richness (Supplementary Table 3). This significance was mainly attributed to the fact that the Shannon indices at the site Neunkirch site varied only marginally, whereas at Chamoson and Saillon they varied rather strongly (Supplementary Figure 1).
Soil Depth but Not Tree Genus Change the Functional Gene Structure
For the functional genes annotated with C- and NCyc, topsoil was clearly separated from deep soil. The NMDS ordination for the CAZy and NCyc gene datasets in topsoil and deep soil of beech and oak sites indicated a strong separation of the functional genes (stress = 0.028 and 0.032, respectively; Figure 1). We further identified the individual soil environmental variables explaining the functional gene structures. Correlations of the genes with soil environmental variables revealed that fine-root biomass, clay content, and plant-available water storage capacity (AWC) were related in the topsoil, whereas in the deep soil the genes correlated closely with pH, base saturation (BS), and fine-earth density.
Figure 1. Non-metric multidimensional scaling (NMDS) ordination based a Euclidean distance matrix calculated on the variance stabilizing transformation (VST) transformed (A) CAZy gene and (B) NCyc gene dataset. Abiotic environmental variables were projected as arrow vectors onto the NMDS ordinations. Plant-available water storage capacity (AWC) per 10 cm soil thickness; BS, base saturation (see also Frey et al., 2021).
The functional gene structure of all predicted genes showed a clear separation between the metagenomes of the topsoil and the deepsoil. This was confirmed by the PERMANOVA, where soil depth explained most of the variance and was a main factor driving sample clustering (Pseudo-F = 5.64, R2 = 0.36, P = 0.003; Table 4). In contrast, site and tree genus were not significant explanatory factors. Similarly, for the CAZy and NCyc datasets, the functional gene structure of the metagenomes in the deep soil was significantly different from that in the topsoil (CAZy genes: Pseudo-F = 4.89, R2 = 0.33, P = 0.002, NCyc genes: Pseudo-F = 5.38, R2 = 0.35, P = 0.002; Table 4).
Table 4. Results of a PERMANOVA testing the effects of various factors for all predicted genes, for the CAZy genes, and for the NCyc genes.
Differentially Abundant Genes in Deep Soils
In order to investigate changes in the abundance of functional genes with site, tree genus and soil depth, we calculated log2-fold changes in all predicted genes, and for the genes annotated with the CAZy and NCyc databases (Table 5). The comparison of all genes between deepsoil and topsoil agglomerated over the three forest sites resulted in 5.4 × 105 over- and 3.1 × 105 underrepresented genes (P < 0.01; Table 5 and Figure 2). Pairwise comparisons of the three forest sites indicated that there were about 10-times fewer over- or underrepresented genes than found in the comparison between topsoil and deepsoil (overrepresented: 9–25 × 103; underrepresented: 9–75 × 103; Table 5). The genes annotated to the CAZy and the NCyc databases showed the most over- and underrepresented genes in the comparison of the topsoil with the deep soil. Thus, most of the differentially abundant genes explained most of the differences between deepsoil and topsoil samples, with the majority of genes were overrepresented in deep soil samples. In contrast, comparisons of the two tree genera indicated <10 over- or underrepresented genes when considering all genes, as well as for the genes annotated with the CAZy and NCyc databases (Table 5).
Table 5. Pairwise comparison and number of significantly (P < 0.01) over- and underrepresented differentially abundant genes out of all predicted genes, of the CAZy genes, and of the NCyc genes.
Figure 2. Volcano plots performed on pairwise comparisons for the factors (A) “tree” beech vs. oak, and (B) “soil depth” deepsoil vs. topsoil. Blue data points represent the genes that are differentially abundant (adjusted P-value < 0.01).
Change in the Abundance of EggNOG, CAZy and NCyc Superfamilies With Depth
Functional genes annotated against EggNOG database were used to assess overall changes in functional genes in response to soil depth and tree genus. The most abundant categories were amino acid transport (E), energy production and conversion (C), signal transduction (T), and replication, recombination, and repair (L). The first three categories had significantly more abundant genes in the deep soil than in the topsoil, whereas soil depth had no effect on the abundance of genes in the L category (Supplementary Table 4). Aggregating over all three levels of functional categories, the categories that increased most significantly in the deep soil were on the levels of “cellular processes and signaling genes” (7 of 10 categories) and “metabolism genes” (5 of 8 categories; Supplementary Table 4). The only category with genes that significantly less abundant in the deep soil was cell motility (N). In contrast, tree genus did not have any significant influence on the number of genes of any category.
Overall, the most abundant genes annotated to the CAZy database were found in the category of glycoside transferases (GT) and glycoside hydrolases (GH), followed by carbohydrate binding modules (CBM). In particular, GT and CBM were significantly more abundant in the deep soil than in the topsoil (Table 6). This indicates, that genes involved in the degradation of cellulose and chitin (CBM) were overrepresented in the deep soil, as well as genes for anabolic processes to build new tissues (GT), which play essential roles in biosynthesis pathways of oligo- and polysaccharides, as well as protein glycosylation and the formation of valuable natural products.
Table 6. Number of CAZy and NCyc genes (n = 3 forest sites).
The most abundant of N-cycling families were those involved in the organic degradation and synthesis (Ods), followed by assimilatory nitrate reduction (Anr), denitrification (Den), and dissimilatory nitrate reduction (Dnr). The family Ods was significantly more abundant in the deep soil, whereas the Anr family only tended to increase in the deep soil (Table 6). Although genes associated to nitrogen fixation (Nif) were significantly abundant in the deep soil, the number of genes was very low overall in the forest soils. No N-cycling families were affected by the tree genus (Table 6).
Functional Categories of Differentially Abundant Genes Annotated to EggNOG, CAZy, and NCyc
In order to investigate changes in the abundance of functional genes with soil depth, we calculated log2-fold changes (LFC) for the genes annotated with EggNOG. Significantly differently abundant microbial functional genes (>1 or <1 LFC) with greater numbers of genes in the deep soil were observed in all functional categories, with the highest values for nuclear structure, translation, and cell cycle control (Supplementary Figure 2). Furthermore, we were particularly interested in the changes in the abundance of functional genes involved in carbohydrate transport and metabolism (Supplementary Figure 3) and amino acid transport and metabolism (Supplementary Figure 4). Although not significant as a whole functional category, a majority of the genes associated with carbohydrate transport and metabolisms were highly significantly more abundant in the deep soil, with the highest positive LFC values corresponding to the genes for glucanotransferase (COG1640), sugar and pyruvate kinases (COG2379, COG0469, COG1080), hydrolase of glycosyl (COG1543), dehydrogenase of phosphogluconate (COG1023), citrate lyase (COG2301), and mannosidase (COG0383; Supplementary Figure 3). The large majority of the genes of the category amino acid transport and metabolism were significantly enhanced in the deep soil, with highest positive LFC of the genes for ABC-type amino acid transporters (COG4177), aminopeptidases (COG2234), amidohydrolases (COG1228), dehydrogenases of alanine and aspartate (COG0686, COG0136), and synthetases (COG0128, COG0337; Supplementary Figure 4).
The total number of differentially abundant genes annotated to CAZy 35,493 (6.99% of the total number of genes), and the abundance of most of these genes differed between soil depths. The comparison between topsoil and deepsoil indicated 10,135 under- and 17,723 overrepresented differentially abundant genes (Table 5). In the family of glycosyl hydrolases (GH family), the genes with the highest positive LFC in the deep soil were involved in the activity of alpha-amylase (GH13 and GH57; involved in starch degradation), alpha-L-rhamnosidase (GH28; pectin), beta-xylosidase (GH43 and GH116; oligosaccharides), 2-mannosidase (GH99; oligosaccharides), alpha-galactosidase (GH27; hemicellulose), alpha-N-acetylgalactosaminidase (GH109; hemicellulose), chitinase (GH18, chitin), and endo-1-4-glucanase (GH5, cellulose) (Figure 3). In the family of the auxiliary activities (AA family), the main activities favored in the deep soil were laccase (AA1; lignin), manganese peroxidase (AA5_2; lignin), and oxidase (AA7; lignin). Genes of the family of carbohydrate binding modules family (CBM family) were also favored in the deep soil, and these genes were mainly involved in the binding of starch (CBM20, CBM26, and CBM41), pectin (CBM32), chitin (CBM12 and CMB14), and cellulose (CBM2, CBM3, CBM6, CBM10, CBM16, and CBM56). In the polysaccharide lyase family (PL family), the genes with the highest positive LFC in the deep soil were pectate lyase (PL10_3; pectin) and heparin lyase (PL12; pectin). Genes of the glycosyl transferase family (GT family) were significantly more abundant in the deep soil, and these genes were mainly involved in the catalyzation of polypeptide transferase (GT27), sialyltransferase (GT29), or arabinofuranosyltransferase (GT95), or in the N-glycosylation of proteins (GT66). In contrast, some genes had a significantly higher abundance in the topsoil, and these genes were mainly involved in the degradation of pectin (PL1: pectin lyase), oligosaccharides (GH43_11; xylan 1,4-β-xylosidase), murein (GH103, peptidoglycan lytic transglycosylase), hemicellulose (GH115; xylan α-1,2-glucuronidase), cellulose (GH5_13; GH17, GH128; cellulase) and cutin (CE5; cutinase).
Figure 3. Under- and overrepresented genes annotated to the CAZy database for the pairwise comparison of topsoil vs. deepsoil. Only significantly (P < 0.01) differentially abundant genes between the two soil depths whose log2-fold change was lower than –2 or higher than +2 are displayed.
For NCyc, the total number of differentially abundant genes was 4,701 (the 5.87% of the total number of genes), and most of these genes varied with soil depth. The comparison between topsoil and deepsoil indicated 1,171 under- and 2,655 overrepresented differentially abundant genes (Table 5). In the comparisons between deepsoil and topsoil, the family of organic degradation and synthesis genes (Ods: ammonification, assimilation) was overrepresented (positive LFC) in the deep soil, with the LFC highest in the deep soil corresponding to the genes ureB, ureA, gs_K00284, glnA and gdh_K00261 (Figure 4). In the family of nitrification genes (Nit), the genes nxrB, amoA_A, amoB_A, amoC_A, amoB_B, and amoC_B had higher LFC values in the deep soil. In the denitrification (Den) and dissimilatory nitrate reduction (Dnr) gene families, norC, nirS, nirK, nirD, narJ, narH, and napC were overrepresented in the deep soil. Furthermore, nirA and narC in the family of assimilatory nitrate reductase (Anr) and the gene pmoA (ammonia/methane oxidation) were more abundant in the deep soil (Figure 4). In contrast, only a few genes were overrepresented in the topsoil, i.e., nrfB, norB, narl, and napB, as well as hao (Figure 4).
Figure 4. Under- and overrepresented genes annotated to the NCyc database for the pairwise comparison of topsoil vs. deepsoil. Only significantly (P < 0.01) differentially abundant genes between the two soil depths whose log2-fold change was lower than –1.5 or higher than +1.5 are displayed.
A summary of the differentially abundant genes annotated to the NCyc database for topsoil and deep soil is shown in Figure 5. It is apparent that mainly the pathways for nitrification (Nit), assimilatory and dissimilatory nitrate reduction (Anr and Dnr), and organic degradation and synthesis (Ods) are increased in the deep soil (Figure 5).
Figure 5. Summary of microbial transformations of N compounds in (red) topsoil and (blue) deepsoil according to the significantly differentially abundant genes between topsoil and deepsoil annotated to the NCyc database following the scheme of Kuypers et al. (2018). Dashed arrows show the hypothetical direction without any differentially abundant genes, thin arrows have one or two increased differentially abundant genes, bold arrows have more than two increased differentially abundant genes between the two soil depths. Nitrogen-transforming processes: Ana, anammox; Anr, assimilatory nitrate reduction; Den, denitrification; Dnr, dissimilatory nitrate reduction; Nit, nitrification; Nif, nitrogen fixation; Ods, organic degradation and synthesis (13)ammonification, 14)assimilation). Enzymes performing N transformations according to Tu et al. (2019): (1) napB, narI; (2) nrfB; (3) norB; (4) hao; (5) narC; (6) napA, napC, narG, narH, narJ; (7) nirA; (8) nirD, nrfC; (9) nirK, nirS; (10) norC; (11) amoA_A, amoB_A, amoC_A, amoB_B, amoC_B; (12) nxrB; (13) gdh_K00260; gdh_K00261, gdh_K00262, gdh_K15371, nmo, ureA, ureB, ureC; (14) asnB, glnA, gs_K00265, gs_K00266, gs_K00284.
Quantitative PCR of C- and N-Transforming Processes
Gene abundances of C- and N-transforming processes were mainly driven by soil depth and not the tree genus. Genes involved in nitrification (archaeal amoA: P = 0.039 and nxrB: P < 0.001) were considerably more abundant in the deep soil compared with the topsoil (Table 7). In contrast, genes involved in denitrification (nosZ: P = 0.002), methane oxidation (pmoA: P = 0.039) and methanogenesis (mrcA: P = 0.002) were significantly less abundant in the deep soil than in the topsoil (Table 7).
Table 7. Mean values of gene abundances of C- and N-transforming processes, and assessed by quantitative PCR (n = 3 forest sites).
Discussion
Altered Microbial Communities in Deep Soil Layers
In deep soil layers, the relative abundance of contigs assigned to Archaea was 5–9 times higher than in the topsoil, while the relative abundance of contigs assigned to Bacteria was not influenced by soil depth. Similarly, both Bacteria and Archaea were unaffected by the tree genus. In deep soil layers, the classes of the chemoautotrophic nitrifying Nitrososphaeria of the archaeal phylum Thaumarchaeota and the anaerobic Methanomicrobia of the archaeal phylum Euryarchaeota were particularly abundant. Archaea have previously been regarded as rare microorganisms that play a limited role in biogeochemical cycles in non-extreme environments (Moissl-Eichinger et al., 2018). Recently, however, Archaea have been found to be distributed widely and comprise a broad diversity of gene repertoires and lifestyles (Martínez-Núñez and Pérez-Rueda, 2016). These authors suggested that Archaea contain more promiscuous enzymes, i.e., the ability of an enzyme to catalyze a random side reaction in addition to its main reaction, providing them with an enzymatic repertoire that enables them to face multiple ecological challenges in harsher environments such as the deep soil.
It is not unexpected to find Nitrososphaera, a class within the phylum Thaumarchaeota, in deep soils. Thaumarchaeota are generally negatively correlated with clay and Corg. Our deep soil layers are characterized by a high sand and low Corg. content, indicating well-drained and oligotrophic conditions. However, ecological niche preferences for ammonia-oxidizing Archaea (AOA) are only found in the literature for topsoil (e.g., Saghaï et al., 2021) and data from deep soils are not available as they are currently very poorly explored. Among the Bacteria, the phyla Nitrospirae, Chloroflexi, and Actinobacteria were significantly more abundant in deep soil. In general, our shotgun metagenomics data are consistent with the results of Frey et al. (2021), who used the same forest sites as in the present study using amplicon sequencing of 16S rRNA genes.
The increased occurrence of some groups of less-known prokaryotes, e.g., Archaea, Chloroflexi, or Nitrospirae at greater soil depths indicates that living conditions for microorganisms change to an oligotrophic and a more oxygen limited environment in deep soil. Such organisms gain their energy from the oxidation of inorganic atoms or molecules as a growth-supporting reductant such as ammonia or nitrite (Hooper and Dispirito, 2013; Kelly and Wood, 2013). Taxa of Thaumarchaeota couple ammonia oxidation with C fixation (Könneke et al., 2014), and taxa of Nitrospirae are able to oxidize nitrite to nitrate or even facilitate the complete cycle from ammonia to nitrite to nitrate (commamox; Daims et al., 2015). Moreover, AOA are favored by oligotrophic soil conditions (Ollivier et al., 2011). Taxa of Methanomicrobia are able to produce methane by using CO2 (Enzmann et al., 2018), and taxa of Actinobacteria form one of the few bacterial groups able to degrade lignin (Abdel-Hamid et al., 2013). However, it also indicates that deep soil microorganisms have to invest more effort into the mobilization of C and N than microorganisms in the topsoil, because they gain their energy from the oxidation of inorganic compounds (lithotrophs), which has a higher cost than required to gain energy from the oxidation of organic compounds (organotrophs) (Hooper and Dispirito, 2013).
Functional Gene Structure and Differentially Abundant Genes
Using shotgun metagenomics, we primarily investigated the potential of functional profiles. The functional gene structure of all predicted genes showed a clear separation between the metagenomes of the topsoil and the deepsoil, and similarly, for the EggNOG, CAZy, and NCyc datasets. Moreover, most of the differentially abundant genes contributed to the difference in functional gene structure observed between the deep soil and the topsoil, with the majority of genes being overrepresented in the deep soil samples (>5 × 105 genes). In contrast, the comparison of the two tree genera beech and oak showed only a few over- or underrepresented genes when all genes were considered, as well as for the genes annotated with the above-mentioned databases. The three forest sites also differed in functional gene structure when compared to each other, but the differentially abundant genes were more than ten times less abundant than when comparing the two soil depths. Therefore, we concluded that the comparison between the topsoil and the deepsoil is the most relevant.
Functional genes annotated against the EggNOG database, which maps predicted genes of general metabolic and cellular functions, showed that the majority of the functional gene categories were significantly overrepresented in the deep soil. In particular, genes of the functional categories such as cell cycle control, nucleotide transport, or amino acid transport and metabolism were highly significantly more abundant. This overall could indicate that in the deep soil altered cellular and metabolic processes occur in the deep soil, which affect the growth of cells and the mobilization of nutrients. Among the most overrepresented genes of the carbohydrate and amino acid transport and metabolism categories were ABC-type amino acid transporters, peptidases, synthetases, transferases, kinases, hydrolases, dehydrogenases, lyases, and mannosidases. Such a dominance of overrepresented functional gene categories is surprising but not unusual, as Perez-Mon et al. (2021) found a similar pattern when comparing permafrost soils with active-layer soils in a high mountain soil in the Swiss Alps. Likewise, in a comparison of cultivated and uncultivated chernozems in Russia, the majority of the functional gene categories were overrepresented in the virgin soil (Gorbacheva et al., 2018). Similarly, Chen et al. (2021) reported, by using the KEGG (Kyoto Encyclopedia of Genes and Genomes) database, that microorganisms in bare soils, compared with those in vegetated soils, have greater relative abundances of genes associated with metabolic functions. They concluded that bare soils harbor a higher proportion of genes that are not available in public reference databases. This is in accordance to our findings that microbial communities in deepsoils are less studied compared to those of topsoils. We assume that oligotrophic microbial communities in deep soils harbor more unannotated genes than copiotrophic microbial communities. These results can be attributed to the fact that oligotrophs are less readily culturable than copiotrophs, and thus their taxonomic and functional information are underrepresented in current reference databases.
Differentially Abundant Genes of C and N-Transforming Processes
Specific information on metabolic properties can be gained from the differentially abundant genes of the CAZy and NCyc databases. The composition of C-degradation genes can reflect the ability of microorganisms to utilize various C-complexes, thereby affecting the accumulation of SOC pools. Gene relative abundance is an adequate predictor of the associated enzyme activity (Trivedi et al., 2016).
With the CAZy database we identified a significant increase in the deep soil in genes for glycosyl transferases (GT) and the carbohydrate binding modules (CBM). The GT genes play essential roles in the biosynthesis pathways of oligo- and polysaccharides, as well as in protein glycosylation and the formation of valuable natural products, and the CBM genes are involved in the binding of carbohydrates, mainly starch, pectin, chitin, and cellulose. It has been considered that CBM genes enhance enzymatic hydrolysis, especially for insoluble substrates (Sidar et al., 2020). In addition, there was a distinct dominance of overrepresented genes in the CAZy database in the deep soil, indicating high metabolic activity of anabolic and catabolic processes. The degradation of oligosaccharides, chitin, pectin, starch, cellulose, and hemicellulose were very common and dominant in the deep soil, whereas those for degradation of lignin did not differ from the topsoil, with the exception of overrepresented genes for multicopper oxidases (AA1, e.g., laccase). It seems that lignin degradation mainly occurs in the topsoil, where most of the fungi live and are active (e.g., Žifčáková et al., 2017; Frey et al., 2021), if we consider fungi as the main lignin decomposers. Thus, we conclude that Bacteria and fungi in deep soil horizons mainly access easily degradable carbohydrates for their survival and growth, as oxygen and nutrients may be limited and at the edge of suitable concentrations (e.g., Lennon, 2020). Microorganisms cannot only promote the release of C into the atmosphere through their catabolic activities but also synthesize labile C into a stable form through their anabolic functions (Liang et al., 2017). We observed that genes for microbial anabolic activities (GTs) were more abundant in the deep soil, however, further studies are required to investigate the contribution of microbial-derived C to SOC sequestration in the deep soil.
With the NCyc database we found a significant increase in the deep soil in genes involved in organic degradation and synthesis. The significant increase in the abundance of N-fixing genes, however, can be neglected because the number of DNA reads was extremely low and, therefore, not relevant. Overrepresented genes in the deep soil were mainly involved in the degradation and synthesis of N compounds such as urea or glutamate, in nitrification and denitrification, and in nitrate reduction. Interestingly, genes involved in nitrification via ammonia monooxygenases (amo) were dominant, with all three subunit types (A, B, C) of Archaea and two subunit types (B, C) of Bacteria overrepresented in the deep soil. It has previously been reported that ammonia-oxidizing Archaea rather prefer areas with low ammonium and have a higher survival rate under conditions of low oxygen than ammonia-oxidizing Bacteria (He et al., 2018). This finding was also confirmed by our quantitative PCR analysis, which indicated that the deep soil harbors significantly higher gene abundances of archaeal amoA. However, such strong seasonal shifts in ammonia and pH are not expected in the investigated calcareous deep soil layers.
In the deep soil, genes for both assimilatory and dissimilatory nitrate reduction were overrepresented. Thus, nitrate is used not only for the assimilation but also for the conservation of energy, indicating that it is used as an electron acceptor for respiration in the (near) absence of oxygen (Kamp et al., 2015). Moreover, dissimilatory nitrate reduction and nitrate storage are physiological life traits that provide the microorganisms with flexibility and resource independence when environments are temporarily anoxic and/or nitrate free (Kamp et al., 2015). Clearly, the overrepresentation of the dissimilatory pathway is an indication of environments that are limited in oxygen even that the investigated soils did not show redoximorphic features. Consequently, microorganisms living in deeper soil horizons have to put more efforts into the mobilization, transformation, and formation of N compounds, which results in greater abundance of genes which can use nitrate not only as a N source but also in energy conservation.
Overall, there are several indications that soils in deep layers may be limited in oxygen and therefore represent a challenging environment for oxic-adapted microbial life (MacDonald et al., 1993). Among the functional genes involved in methane cycling, no methyl-coenzyme M reductase genes (the primary gene involved in methanogenesis, mcrA) were detected in the metagenomes. By using quantitative PCR, however, we found an overall low abundance of mcrA genes at both soil depths, with a higher abundance of mcrA genes in the topsoil. Methanogens use CO2 as an electron acceptor during anaerobic respiration and produce methane. Therefore, these low abundances of mcrA genes could indicate oxygen-limited but not fully anaerobic conditions (Frey et al., 2011). This corresponds to the findings that the abundance of the gene for methane monooxygenase (pmoA), an enzyme that metabolizes methane for C-metabolism and energy, was also significantly less abundant in the deep soil. Our findings suggest that methane as a substrate for pmoA may not be produced at higher concentrations in the deep soils studied, although methane could be oxidized anaerobically using alternative electron acceptors, e.g., nitrite, nitrate, and ferric iron (Cai et al., 2018; Shen et al., 2019).
As no measurements of oxygen concentrations in topsoils and deepsoils were carried out, one can only speculate about the interpretation of the complex interactions between methane-cycling organisms, AOA and oxygen concentrations. However, two points need to be considered. Firstly, due to the sandy and well-drained soils in deep soils, it can be assumed that deepsoils are not strongly oxygen-limited and therefore contain a lower abundance of mcrA genes but a higher abundance of AOA due to the more oligotrophic conditions. Such a condition is not unusual. Secondly, due to the higher microbial activity in topsoils resulting in higher oxygen consumption, and in combination with a high clay content, there are more anaerobic microniches in the topsoil. Thus, a higher abundance of mcrA genes can be expected in topsoils. As for methanogens, which are strictly anaerobic organisms, the presence of mcrA genes has previously been shown in oxic soils in anaerobic microniches (Frey et al., 2011; Li et al., 2021). In addition, anaerobic methanotrophs have been reported that can oxidize methane produced by methanogens (Luesken et al., 2012). Anaerobic methanotrophs might be present in the deeper soil layer, most likely due to their sensitivity to oxygen. This also explains the relative high abundance of pmoA genes in both top and deep soils.
Conclusion
The metagenomic analysis of deep soil layers of beech and oak forests in Switzerland revealed that tree genus and forest site factors have only a minor effect on the composition of the microbial gene repertoire, while soil depth had strong effects. Soil depth decisively changes chemical conditions and resource availability leading to more favorable conditions for Archaea than for Bacteria. The strong increase in the archaeal domains in the deep soil layers also indicates that biogeochemical processes and cycling are likely to be affected by soil depth. In deep soils, genes for carbohydrate-active enzymes involved in catalyzing the transfer of saccharide moieties are overrepresented, as are enzymes involved in binding carbohydrates such as chitin or cellulose. The greater abundance of these genes in the deep soil demonstrates that greater efforts have been made for enzymatic hydrolysis, especially for insoluble substrates. Furthermore, overrepresented genes of the N cycle in the deep soils are involved in the degradation and synthesis of N compounds, in nitrification, in denitrification, and in nitrate reduction. Consequently, the entire N-cycle is affected and N is not only used for assimilation but also for energy conservation, indicating conditions of low oxygen, even if the soils were not hydromorphic. Overall, metagenomics proved to be an appropriate approach to gain insights into biogeochemical processes that may change with altered soil properties, such as those associated with greater soil depths. Our study is one of the first to offer insight into the functional diversity of heterogeneous microbial assemblages in deep forest soils. Our data on the functional genetic potential of the soil microbiomes in deep soils provide information about their metabolic capacities, enabling modeling of depth-specific biogeochemical processes that may change as a result of altered soil conditions.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/783873, bioproject/PRJNA783873.
Author Contributions
BF and IB designed the microbial study and wrote the main parts of the manuscript. LW established and maintained the forest sites. LW and IB provided soil and fine-root data. BS and WQ performed genetic analyses in the lab. BF, GV, and IB performed statistical analyses. All authors contributed to the final version of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge Roger Köchli and Rafael Boduryan for processing soils and roots, the WSL Central Laboratory for conducting soil chemical analyses, and Stephan Zimmermann for providing data from the WSL soil database. We also acknowledge Microsynth AG, Balgach, Switzerland for completing shotgun sequencing. We thank Melissa Dawes for help editing the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.828977/full#supplementary-material
Footnotes
References
Abdel-Hamid, A. M., Solbiati, J. O., and Cann, I. K. O. (2013). Insights into lignin degradation and its potential industrial applications. Adv. Appl. Microbiol. 82, 1–28. doi: 10.1016/B978-0-12-407679-2.00001-6
Adamczyk, M., Hagedorn, F., Wipf, S., Donhauser, J., Vittoz, P., Rixen, C., et al. (2019). The soil microbiome of Gloria mountain summits in the Swiss Alps. Front. Microbiol. 10:1080. doi: 10.3389/fmicb.2019.01080
Andújar, C., Arribas, P., and Vogler, A. P. (2017). Terra incognita of soil biodiversity: unseen invasions under our feet. Mol. Ecol. 26, 3087–3089. doi: 10.1111/mec.14112
Anwar, M. Z., Lanzen, A., Bang-Andreasen, T., and Jacobsen, C. S. (2019). To assemble or not to resemble–a validated Comparative Metatranscriptomics Workflow (CoMW). Gigascience 8:giz096. doi: 10.1093/gigascience/giz096
Bodelier, P. L. E., and Steenbergh, A. K. (2014). Interactions between methane and the nitrogen cycle in light of climate change. Curr. Opin. Environ. Sustain. 9-10, 26–36. doi: 10.1016/j.cosust.2014.07.004
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Brewer, T. E., Aronson, E. L., Arogyaswamy, K., Billings, S. A., Botthoff, J. K., Campbell, A. N., et al. (2019). Ecological and genomic attributes of novel bacterial taxa that thrive in subsurface soil horizons. mBio 10:e01318-19. doi: 10.1128/mBio.01318-19
Brunner, I., Plötze, M., Rieder, S., Zumsteg, A., Furrer, G., and Frey, B. (2011). Pioneering fungi from the Damma glacier forefield in the Swiss Alps can promote granite weathering. Geobiology 9, 266–279. doi: 10.1111/j.1472-4669.2011.00274.x
Cabello, P., Roldan, M. D., and Moreno-Vivian, C. (2004). Nitrate reduction and the nitrogen cycle in archaea. Microbiology 150, 3527–3546. doi: 10.1099/mic.0.27303-0
Cai, C., Leu, A. O., Xie, G.-J., Guo, J., Feng, Y., Zhao, J.-X., et al. (2018). A methanotrophic archaeon couples anaerobic oxidation of methane to Fe(III) reduction. ISME J. 12, 1929–1939. doi: 10.1038/s41396-018-0109-x
Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The carbohydrate-active enzymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663
Chen, Y., Neilson, J. W., Kushwaha, P., Maier, R. M., and Barberán, A. (2021). Life-history strategies of soil microbial communities in an arid ecosystem. ISME J. 15, 649–657. doi: 10.1038/s41396-020-00803-y
Crowther, T. W., van den Hoogen, J., Wan, J., Mayes, M. A., Keiser, A. D., Mo, L., et al. (2019). The global soil community and its influence on biogeochemistry. Science 365:eaav0550. doi: 10.1126/science.aav0550
Daims, H., Lebedeva, E. V., Pjevac, P., Han, P., Herbold, C., Albertsen, M., et al. (2015). Complete nitrification by Nitrospira bacteria. Nature 528, 504–509. doi: 10.1038/nature16461
DeBruyn, J. M., Nixon, L. T., Fawaz, M. N., Johnson, A. M., and Radosevich, M. (2011). Global biogeography and quantitative seasonal dynamics of Gemmatimonadetes in soil. Appl. Environ. Microbiol. 77, 6295–6300. doi: 10.1128/AEM.05005-11
Donhauser, J., Qi, W., Bergk-Pinto, B., and Frey, B. (2021). High temperature enhance the microbial genetic potential to recycle C and N from necromass in high-mountain soils. Glob. Change Biol. 27, 1365–1386. doi: 10.1111/gcb.15492
Dove, N. C., Arogyaswamy, K., Billings, S. A., Botthoff, J. K., Carey, C. J., Cisco, C., et al. (2020). Continental-scale patterns of extracellular enzyme activity in the subsoil: an overlooked reservoir of microbial activity. Environ. Res. Lett. 15:1040a1. doi: 10.1088/1748-9326/ABB0B3
Eilers, K. G., Debenport, S., Anderson, S., and Fierer, N. (2012). Digging deeper to find unique microbial communities: the strong effect of depth on the structure of bacterial and archaeal communities in soil. Soil Biol. Biochem. 50, 58–65. doi: 10.1016/j.soilbio.2012.03.011
Eldridge, D. J., Delgado-Baquerizo, M., Woodhouse, J. N., and Neilan, B. A. (2016). Mammalian engineers drive soil microbial communities and ecosystem functions across a disturbance gradient. J. Anim. Ecol. 85, 1636–1646. doi: 10.1111/1365-2656.12574
Enzmann, F., Mayer, F., Rother, M., and Holtmann, D. (2018). Methanogens: biochemical background and biotechnological applications. AMB Express 8:1. doi: 10.1186/s13568-017-0531-x
Falkowski, P. G., Fenchel, T., and Delong, E. F. (2008). The microbial engines that drive Earth’s biogeochemical cycles. Science 320, 1034–1039. doi: 10.1126/science.1153213
Frey, B., Carnol, M., Dharmarajah, A., Brunner, I., and Schleppi, P. (2020). Only minor changes in the soil microbiome of a sub-alpine forest after 20 years of moderately increased nitrogen loads. Front. For. Glob. Change 3:77. doi: 10.3389/ffgc.2020.00077
Frey, B., Niklaus, P. A., Kremer, J., Lüscher, P., and Zimmermann, S. (2011). Heavy-machinery traffic impacts methane emissions as well as methanogen abundance and community structure in oxic forest soils. Appl. Environ. Microbiol. 77, 6060–6068. doi: 10.1128/AEM.05206-11
Frey, B., Rieder, S. R., Brunner, I., Plötze, M., Koetzsch, S., Lapanje, A., et al. (2010). Weathering-associated bacteria from the Damma glacier forefield: physiological capabilities and impact on granite dissolution. Appl. Environ. Microbiol. 76, 4788–4796. doi: 10.1128/AEM.00657-10
Frey, B., Walthert, L., Perez-Mon, C., Stierli, B., Köchli, R., Dharmarajah, A., et al. (2021). Deep soil layers of drought-exposed forests harbor poorly known bacterial and fungal communities. Front. Microbiol. 12:674160. doi: 10.3389/fmicb.2021.674160
Gee, G. W., and Bauder, J. W. (1986). “Particle-size analysis,” in Methods of Soil Analysis, Part 1, Physical and Mineralogical Methods, Vol. 9, ed. A. Klute (Madison WI: American Society of Agronomy), 383–411. doi: 10.2136/sssabookser5.1.2ed.c15
Gorbacheva, M. A., Melnikova, N. V., Chechetkin, V. R., Kravatsky, Y. V., and Tchurikov, N. A. (2018). DNA sequencing and metagenomics of cultivated and uncultivated chernozems in Russia. Geoderma Reg. 14:e00180. doi: 10.1016/j.geodrs.2018.e00180
He, H., Zhen, Y., Mi, T., Fu, L., and Yu, Z. (2018). Ammonia-oxidizing archaea and bacteria differentially contribute to ammonia oxidation in sediments from adjacent waters of Rushan Bay, China. Front. Microbiol. 9:116. doi: 10.3389/fmicb.2018.00116
Hooper, A. B., and Dispirito, A. A. (2013). Chemolithotrophy. Encycl. Biol. Chem. 1, 486–492. doi: 10.1016/B978-0-12-378630-2.00219-X
Huerta-Cepas, J., Forslund, K., Coelho, L. P., Szklarczyk, D., Jensen, L. J., von Mering, C., et al. (2017). Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol. Biol. Evol. 34, 2115–2122. doi: 10.1093/molbev/msx148
Huerta-Cepas, J., Szklarczyk, D., Forslund, K., Cook, H., Heller, D., Walter, M. C., et al. (2016). eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, D286–D293. doi: 10.1093/nar/gkv1248
IUSS Working Group WRB (2007). World Reference Base for Soil Resources 2006, First Update 2007. World Soil Resources Reports No. 103. Rome: FAO.
Jensen, L. J., Julien, P., Kuhn, M., Mering, C., von Muller, J., Doerks, T., et al. (2008). eggNOG: automated construction and annotation of orthologous groups of genes. Nucleic Acids Res. 36, D250–D254. doi: 10.1093/nar/gkm796
Jiang, Y., Lei, Y., Qin, W., Korpelainen, H., and Li, C. (2019). Revealing microbial processes and nutrient limitation in soil through ecoenzymatic stoichiometry and glomalin-related soil proteins in a retreating glacier forefield. Geoderma 338, 313–324. doi: 10.1016/j.geoderma.2018.12.023
Jobbágy, E. G., and Jackson, R. B. (2000). The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 10, 423–436. doi: 10.1890/1051-0761(2000)010[0423:tvdoso]2.0.co;2
Kamp, A., Høgslund, S., Risgaard-Petersen, N., and Stief, P. (2015). Nitrate storage and dissimilatory nitrate reduction by eukaryotic microbes. Front. Microbiol. 6:1492. doi: 10.3389/fmicb.2015.01492
Kelly, D. P., and Wood, A. P. (2013). “The chemolithotrophic prokaryotes,” in The Prokaryotes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer). doi: 10.1007/978-3-642-30123-0_63
Könneke, M., Schubert, D. M., Brown, P. C., Hügler, M., Standfest, S., Schwander, T., et al. (2014). Ammonia-oxidizing archaea use the most energy-efficient aerobic pathway for CO2 fixation. Proc. Natl. Acad. Sci. U.S.A. 111, 8239–8244. doi: 10.1073/pnas.1402028111
Kowalchuk, G. A., and Stephen, J. R. (2001). Ammonia-oxidizing bacteria: a model for molecular microbial ecology. Annu. Rev. Microbiol. 55, 485–529. doi: 10.1146/annurev.micro.55.1.485
Kuypers, M. M. M., Marchant, H. K., and Kartal, B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol. 16:263. doi: 10.1038/nrmicro.2018.9
Li, D., Liu, C.-M., Luo, R., Sadakane, K., and Lam, T.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033
Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 3:13033997.
Li, Q., Liu, Y., Gu, Y., Guo, L., Huang, Y., Zhang, J., et al. (2020). Ecoenzymatic stoichiometry and microbial nutrient limitations in rhizosphere soil along the Hailuogou Glacier forefield chronosequence. Sci. Total Environ. 704:135413. doi: 10.1016/j.scitotenv.2019.135413
Li, Q., Peng, C., Zhang, J., Yongfu, L., and Song, X. (2021). Nitrogen addition decreases methane uptake caused by methanotroph and methanogen imbalances in a Moso bamboo forest. Sci. Rep. 11:5578. doi: 10.1038/s41598-021-84422-3
Liang, C., Schimel, J. P., and Jastrow, J. D. (2017). The importance of anabolism in microbial control over soil carbon storage. Nat. Microbiol. 2:17105 doi: 10.1038/nmicrobiol.2017.105
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Liu, G., Chen, L., Deng, Q., Shi, X., Lock, T. R., Kallenbach, R. L., et al. (2020). Vertical changes in bacterial community composition down to a depth of 20 m on the degraded Loess Plateau in China. Land Degrad. Dev. 31, 1300–1313. doi: 10.1002/ldr.3542
Lladó, S. F., Vìtrovskı, T., and Baldrian, P. (2019). Tracking of the activity of individual bacteria in temperate forest soils shows guild-specific responses to seasonality. Soil Biol. Biochem. 135, 275–282. doi: 10.1016/j.soilbio.2019.05.010
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Luesken, F. A., Wu, M. L., den Camp, H. J. M. O., Keltjens, J. T., Stunnenberg, H., Francoijs, K.-J., et al. (2012). Effect of oxygen on the anaerobic methanotroph ‘Candidatus Methylomirabilis oxyfera’: kinetic and transcriptional analysis. Environ. Microbiol. 14, 1024–1034. doi: 10.1111/j.1462-2920.2011.02682.x
MacDonald, J. D., Costello, L. R., and Berger, T. (1993). An evaluation of soil aeration status around healthy and declining oaks in an urban environment in California. J. Arboricul. 19, 209–219.
Martens-Habbena, W., Berube, P. M., Urakawa, H., de la Torre, J. R., and Stahl, D. A. (2009). Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461, 976–979. doi: 10.1038/nature08465
Martínez-Núñez, M. A., and Pérez-Rueda, E. (2016). Do lifestyles influence the presence of promiscuous enzymes in bacteria and Archaea metabolism? Sustain. Chem. Process. 4:3. doi: 10.1186/s40508-016-0047-8
Menzel, P., Ng, K. L., and Krogh, A. (2016). Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat. Commun. 7:11257. doi: 10.1038/ncomms11257
Moissl-Eichinger, M., Pausan, M., Taffner, J., Berg, G., Bang, C., and Schmitz, R. A. (2018). Archaea are interactive components of complex microbiomes. Trends Microbiol. 26, 70–85. doi: 10.1016/j.tim.2017.07.004
Oksanen, J., Blanchet, G. F., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). vegan: Community Ecology Package. R package version 2.5-6.
Ollivier, J., Töwe, S., Bannert, A., Hai, B., Kastl, E.-M., Meyer, A., et al. (2011). Nitrogen turnover in soil and global change. FEMS Microbiol. Ecol. 78, 3–16. doi: 10.1111/j.1574-6941.2011.01165.x
Perez-Mon, C., Qi, W., Vikram, S., Frossard, A., Makhalanyane, T., Cowan, D., et al. (2021). Shotgun metagenomics reveals distinct functional diversity and metabolic capabilities between 12?000-year-old permafrost and active layers on Muot da Barba Peider (Swiss Alps). Microb. Genom. 7:000558. doi: 10.1099/mgen.0.000558
R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Remund, J., Rihm, B., and Huguenin-Landl, B. (2014). Klimadaten für die Waldmodellierung für das 20. und 21. Jahrhundert. Birmensdorf: Eidg. Forschungsanstalt für Wald, Schnee und Landschaft WSL. doi: 10.3929/ethz-a-010693673
Rime, T., Hartmann, M., Brunner, I., Widmer, F., Zeyer, J., and Frey, B. (2015). Vertical distribution of the soil microbiota along a successional gradient in a glacier forefield. Mol. Ecol. 24, 1091–1108. doi: 10.1111/mec.13051
Rime, T., Hartmann, M., and Frey, B. (2016). Potential sources of microbial colonizers in an initial soil ecosystem after retreat of an alpine glacier. ISME J. 10, 1625–1641. doi: 10.1038/ismej.2015.238
Robin, A., Pradier, C., Sanguin, H., Mahé, F., Lambais, G. R., de Araujo Pereira, A. P., et al. (2019). How deep can ectomycorrhizas go? A case study on Pisolithus down to 4 meters in a Brazilian eucalypt plantation. Mycorrhiza 29, 637–648. doi: 10.1007/s00572-019-00917-y
Saghaï, A., Banjeree, S., Degrune, F., Edlinger, A., García-Palacios, P., Garland, G., et al. (2021). Diversity of archaea and niche preferences among putative ammonia-oxidizing Nitrososphaeria dominating across European arable soils. Environ. Microbiol. doi: 10.1111/1462-2920.15830
Shen, L.-D., Ouyang, L., Zhu, Y., and Trimmer, M. (2019). Active pathways of anaerobic methane oxidation across contrasting riverbeds. ISME J. 13, 752–766. doi: 10.1038/s41396-018-0302-y
Sidar, A., Albuquerque, E. D., Voshol, G. P., Ram, A. F. J., Vijgenboom, E., and Punt, P. J. (2020). Carbohydrate binding modules: diversity of domain architecture in amylases and cellulases from filamentous microorganisms. Front. Bioeng. Biotechnol. 8:871. doi: 10.3389/fbioe.2020.00871
Spohn, M., Klaus, K., Wanek, W., and Richter, A. (2016). Microbial carbon use efficiency and biomass turnover times depending on soil depth – implications for carbon cycling. Soil Biol. Biochem. 96, 74–81. doi: 10.1016/j.soilbio.2016.01.016
Tamames, J., Cobo-Simón, M., and Puente-Sánchez, F. (2019). Assessing the performance of different approaches for functional and taxonomic annotation of metagenomes. BMC Genomics 20:960. doi: 10.1186/s12864-019-6289-6
Teepe, R., Dilling, H., and Beese, F. (2003). Estimating water retention curves of forest soils from soil texture and bulk density. J. Soil Sci. Plant Nutr. 166, 111–119. doi: 10.1002/jpln.200390001
Thiel, V. (2011). “Chemolithotrophy,” in Encyclopedia of Geobiology, eds J. Reitner and V. Thiel (Dordrecht: Springer). doi: 10.1007/978-1-4020-9212-1_53
Trivedi, C., Reich, P. B., Maestre, F. T., Hu, H. W., Singh, B. K., and Delgado-Baquerizo, M. (2019). Plant-driven niche differentiation of ammonia-oxidizing bacteria and archaea in global drylands. ISME J. 13, 2727–2736. doi: 10.1038/s41396-019-0465-1
Trivedi, P., Delgado-Baquerizo, M., Trivedi, C., Hu, H. W., Anderson, I. C., Jeffries, T. C., et al. (2016). Microbial regulation of the soil carbon cycle: evidence from gene-enzyme relationships. ISME J. 10, 2593–2604. doi: 10.1038/ismej.2016.65
Tu, Q., Lin, L., Cheng, L., Deng, Y., and He, Z. (2019). NCycDB: a curated integrative database for fast and accurate metagenomic profiling of nitrogen cycling genes. Bioinformatics 35, 1040–1048. doi: 10.1093/bioinformatics/bty741
Uroz, S., Calvaruso, C., Turpault, M. P., and Frey-Klett, P. (2009). Mineral weathering by bacteria: ecology, actors and mechanisms. Trends Microbiol. 17, 378–387. doi: 10.1016/j.tim.2009.05.004
Uroz, S., Kelly, L. C., Turpault, M. P., Lepleux, C., and Frey-Klett, P. (2015). The mineralosphere concept: mineralogical control of the distribution and function of mineral-associated bacterial communities. Trends Microbiol. 23, 751–762. doi: 10.1016/j.tim.2015.10.004
Varliero, G., Anesio, A. M., and Barker, G. L. A. (2021). A taxon-wise insight into rock weathering and nitrogen fixation functional profiles of proglacial systems. Front. Microbiol. 12:627437. doi: 10.3389/fmicb.2021.627437
Vaser, R., Dario Pavlović, D., and Šikić, M. (2016). SWORD-a highly efficient protein database search. Bioinformatics 32, i680–i684. doi: 10.1093/bioinformatics/btw445
Walthert, L., Graf, U., Kammer, A., Luster, J., Pezzotta, D., Zimmermann, S., et al. (2010). Determination of organic and inorganic carbon, d13C, and nitrogen in soils containing carbonates after acid fumigation with HCl. J. Soil Sci. Plant Nutr. 173, 207–216. doi: 10.1002/jpln.200900158
Wang, Y. L., Wang, Q., Yuan, R., Sheng, X. F., and He, L. Y. (2019). Isolation and characterization of mineral-dissolving bacteria from different levels of altered mica schist surfaces and the adjacent soil. World J. Microbiol. Biotechnol. 35:2. doi: 10.1007/s11274-018-2573-x
Wongfun, N., Plötze, M., Furrer, G., and Brandl, H. (2014). Weathering of granite from the Damma Glacier area: the contribution of cyanogenic bacteria. Geomicrobiol. J. 31, 93–100. doi: 10.1080/01490451.2013.802396
Xi, J., Wei, M., and Tang, B. (2018). Differences in weathering pattern, stress resistance and community structure of culturable rock-weathering bacteria between altered rocks and soils. RSC Adv. 8, 14201–14211. doi: 10.1039/C8RA01268G
Xu, X., Thornton, P. E., and Post, W. M. (2013). A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Glob. Ecol. Biogeogr. 22, 737–749. doi: 10.1111/geb.12029
Yost, J. L., and Hartemink, A. E. (2020). How deep is the soil studied – an analysis of four soil science journals. Plant Soil 452, 5–18. doi: 10.1007/s11104-020-04550-z
Zhang, X., Liu, W., Schloter, M., Zhang, G., Chen, Q., Huang, J., et al. (2013). Response of the abundance of key soil microbial nitrogen-cycling genes to multi-factorial global changes. PLoS One 8:e76500. doi: 10.1371/journal.pone.0076500
Zhu, A., Ibrahim, J. G., and Love, M. I. (2019). Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 5, 2084–2092. doi: 10.1093/bioinformatics/bty895
Zhu, W., Lomsadze, A., and Borodovsky, M. (2010). Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 38:e132. doi: 10.1093/nar/gkq275
Keywords: forest, C and N cycles, CAZy, subsoil, metagenomics
Citation: Frey B, Varliero G, Qi W, Stierli B, Walthert L and Brunner I (2022) Shotgun Metagenomics of Deep Forest Soil Layers Show Evidence of Altered Microbial Genetic Potential for Biogeochemical Cycling. Front. Microbiol. 13:828977. doi: 10.3389/fmicb.2022.828977
Received: 04 December 2021; Accepted: 11 February 2022;
Published: 01 March 2022.
Edited by:
Marcus A. Horn, Leibniz University Hannover, GermanyReviewed by:
Boyke Bunk, Leibniz Institute DSMZ, GermanyAchim Quaiser, University of Rennes 1, France
Copyright © 2022 Frey, Varliero, Qi, Stierli, Walthert and Brunner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Beat Frey, beat.frey@wsl.ch