Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 29 October 2025

Sec. Terrestrial Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1654005

This article is part of the Research TopicUnveiling Microbiome Interactions and Functions in Soil HotspotsView all 17 articles

Beyond one-size-fits-all approach: How do various harvesting strategies shape soil microbial diversity in Larix gmelinii (Daxinganling larch) forests of the Greater Khingan Mountains?


Yufeng Wang,Yufeng Wang1,2Na Ta,Na Ta1,2Hao Zhang,Hao Zhang1,2Min Li,Min Li1,2Shengwei Liu,Shengwei Liu1,2Jiaxing Gong,Jiaxing Gong1,2Yake SongYake Song1Rula Sa*Rula Sa1*
  • 1College of Forestry, Inner Mongolia Agricultural University, Hohhot, China
  • 2National Field Scientific Observation and Research Station of Greater Khingan Forest Ecosystem, Genhe, China

This study analyzed the effects of six harvesting strategies, including primary forests (PF), shelterwood cutting (SC), clear cutting (CC), optional cutting with low intensity (OCL), optional cutting with moderate intensity (OCM), and optional cutting with high intensity (OCH), on soil microbial diversity in Daxinganling larch forests of the Greater Khingan Range. The results showed that the dominant bacterial and fungal phyla were similar across the different harvesting strategies. As the intensity of OC increased, the relative abundances of Proteobacteria and Basidiomycota increased, whereas that of Ascomycota decreased. The highest bacterial alpha diversity was observed in the PF sample plots, whereas the highest fungal alpha diversity was observed in the SC sample plots. OCH significantly reduced bacterial alpha diversity (p < 0.05), and a negative correlation was observed between OC intensity and bacterial alpha diversity. Harvesting strategies had no significant effect on bacterial or fungal beta diversity. In the harvesting strategy sample plots, 14 biological markers were enriched, including the bacterial family SC_l_84 and fungal genus Coniochaeta. Soil nitrogen, phosphorus, and other physical and chemical properties were significantly correlated with different microbial markers. Soil bacterial and fungal communities have diverse genetic and ecological functions. The bacterial and fungal community networks in the PF sample plots were the most complex and stable. OC reduced the complexity and stability of soil bacterial community networks but had the opposite effect on fungal communities. This study preliminarily analyzed the effects of different harvesting strategies on soil microbial diversity in the Greater Khingan Range Larix gmelinii forest, which has practical significance for the functional recovery and protection of the Greater Khingan Range forest ecosystem.

1 Introduction

As the largest terrestrial ecosystem on Earth, forests are irreplaceable for maintaining biodiversity, water conservation, and regulating the global carbon balance (Nandal et al., 2023; Hu et al., 2024). Soil microbial communities, dominated by bacteria and fungi, play crucial roles in maintaining forest ecosystems. They profoundly influence ecosystem productivity and stability by driving organic matter decomposition, regulating nutrient cycling, and forming complex symbiotic networks with plant root systems (Delgado-Baquerizo et al., 2020; Li et al., 2025). Harvesting strategies, one of the main forms of anthropogenic disturbance in forests, mainly include SC, OC, CC, and thinning, which aim to obtain timber resources or facilitate forest regeneration (Bont et al., 2022; Repo et al., 2024). However, the impact of harvesting strategies on microbial diversity in forest soils cannot be ignored.

Harvesting activities affect soil microbial community structure and function by directly altering soil physicochemical properties [pH (Melo et al., 2021)], organic matter content (Qiang et al., 2023), temperature, and humidity (Liu et al., 2018a), nutrient substances (Morrison et al., 2016; Liu et al., 2018b), understory microenvironments [light intensity (Zhu et al., 2021)], and plant-microbe interactions (Jog et al., 2014; Cordovez et al., 2021; Wang et al., 2023; Qu et al., 2024).

Different harvesting strategies influence the microbial diversity in forest soils. Reasonable harvesting strategies (moderate thinning and SC) have the potential to maintain or promote soil microbial ecosystems while enhancing forest economic benefits. Thinning contributes to the construction of healthy and stable forest ecosystems by optimizing stand density, wood species composition, and age structure (Zhou et al., 2016). Thinning experiments of different intensities in Chinese fir plantations in the Lishui District, Nanjing Province, China, and oil pine plantations on the Loess Plateau showed that thinning at an appropriate intensity can increase the diversity of plants under the forest and soil microorganisms (Chen et al., 2015; Dang et al., 2018). In a study of secondary forests of Betula platyphylla Suk in the Greater Khingan Mountains, SC sample plots had the highest soil bacterial and fungal beta diversity indices, whereas the most intensively damaged fire sample plots had the lowest fungal beta diversity (Zhai et al., 2023). In contrast, CC generally have a negative impact on soil microbial diversity. A meta-analysis showed that CC resulted in 14–33% and 20–40% declines in soil microbial and fungal biomass, respectively (Holden and Treseder, 2013). A study in British Columbia, Canada, showed that soil bacterial and fungal community composition was significantly lower after CC than before, even after 10 or 15 years of plantation-forest restoration. Whole-plant harvesting (with no residues) results in significantly lower soil microbial biomass and diversity than stem harvesting because of the drastically reduced inputs of litter (Hartmann et al., 2009, 2012; Mushinski et al., 2018).

Daxinganling larch (Larix gmelinii) is a typical forest species in the cold-temperate region of China and a major community-building species in the forest areas of the Greater Khingan Mountains (Yang et al., 2020; Cheng et al., 2023). It is currently planted in an area of 3 × 104km2, accounting for more than 30% of China’s planted forest area (Yang et al., 2010, 2018). However, since the last century, irrational strategies, such as CC and OCH, have proliferated in the region to meet the demand for timber, causing serious damage to forest soil microbial diversity and ecosystem function. In recent years, researchers have conducted preliminary explorations of the impact of harvesting on soil microorganisms in Daxinganling larch. One study showed that CC resulted in a significant decrease in soil microbial carbon and N content, relative abundance of fungi, and Simpson’s diversity index, whereas SC did not (Yang et al., 2020). Another study demonstrated that, compared to PF, Thinning significantly increased the carbon and N content of soil microorganisms in artificial forests, as well as the Shannon and Simpson indices of fungi, but had little effect on bacterial community alpha diversity (p > 0.05) (Li et al., 2021). However, the above studies did not systematically analyze the effects of multiple harvesting strategies (SC, OC, and CC) on the diversity of microbial communities in the PF of Daxinganling larch, and research on different OC intensities remains limited.

This study used Daxinganling larch as the research object and preliminarily analyzed the adaptation and response mechanisms of soil microbial diversity to different harvesting strategies. We assumed that (1) different harvesting strategies would affect the soil microbial community, including changes in community composition and diversity. (2) Key biomarkers that respond to different harvesting strategies are related to soil physicochemical properties. Therefore, this study aims (1) to systematically clarify the effects of different harvesting strategies on the alpha and beta diversity of soil bacterial and fungal communities. (2) To explore the effects of different harvesting strategies on the molecular interaction network, genes, and ecological functions of soil microorganisms. (3) To analyze the correlation between the physical and chemical properties of the soil and key biomarkers in sample plots with different harvesting strategies.

2 Materials and methods

2.1 Location of the study areas

The study was conducted in Inner Mongolia (50° 20’–52° 30’ N, 120° 12’–122° 55’ E), specifically within the Chaocha Forestry. The topography is predominantly characterized by slopes, followed by gentle slopes, with a mean elevation of approximately 950 m (maximum elevation 1210 m). The region experiences a cold-temperate humid monsoon climate, with an average annual temperature of −5.3 °C and an average annual precipitation of 500 mm. The dominant soil type is brown coniferous forest soil, which is sandy, gravelly, and acidic. The forest exhibited high species richness, with Daxinganling larch as the dominant species and broad-leaved conifers, such as Daxinganling white birch (Betula platyphylla Sukaczev) and Daxinganling Scotch pine (Pinus sylvestris var. mongolica litv). The understory vegetation is well-developed, and the main shrubs include marsh labrador tea (R. tomentosum), dahurian rhododendron (R. mucronulatum), and Mongolian meadowsweet (S. pubesscens), which mainly include wintergreen (Pyrola calliantha), Chinese pennisetum (Pennisetum alopecuroides), and great burnet (Sanguisorba officinalis L.).

2.2 Study design

In 2005, six harvesting strategies were set up in the primary forest of Daxinganling larch in Chaocha Forestry: PF, SC, OCL, OCM, OCH, and CC. Each harvesting strategy corresponded to an independent 50 × 50 m plot, totaling six sample plots in the study area. All sample plots had the same ground conditions. Eighteen years after the resumption of harvesting (2023), soil samples were collected, and their physical and chemical properties and microbial diversity were analyzed. Litter was collected and weighed. The characteristics of the sample plots are listed in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of sample plots information for different harvesting strategies.

2.3 Soil sample collection

Six soil sampling points were evenly distributed diagonally across each plot. First, three sampling points were randomly selected, and during sampling, withered material and debris on the soil surface were removed. Soil samples were collected at a depth of 0–20 cm using sterile equipment. After removing visible plant roots, stones, and gravel, the soil samples were air-dried under ventilated and light-free conditions. The air-dried samples were ground and mixed in a sterile mortar, passed through a 2 mm aperture nylon sieve, and then transferred to sterile self-sealing bags and numbered for storage. Part of the samples was used to determine total N, total phosphorus, total potassium, soil organic carbon (SOC), and pH of the soil, while the other part of the samples was immediately placed in an ultra-low-temperature refrigerator at -80°C for storage. Finally, high-throughput sequencing was performed. The remaining three sampling points were used to determine the mass of the fallen leaves in the forest. The area of each sample was 1 m × 1 m in size. All litter in the sample plots was harvested and placed in self-sealing bag.

2.4 High-throughput sequencing analysis of soil microorganisms

Total DNA was extracted from the soil samples using the cetyltrimethylammonium bromide method. The purity and concentration were assessed using a NanoDrop 2000 Ultra-Micro Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and qualified samples were stored at -20°C for subsequent use. The V3 highly variable region of the bacterial genomic 16S rRNA was amplified using PCR (Polymerase Chain Reaction) with the primers 515F (5′-GTGCCAGCMGCCGCGGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). The PCR mix was added to 15 μL Phusion High-Fidelity PCR Master Mix, 0.2 μM primer, and 10 ng genomic DNA template. The PCR conditions were as follows: pre-denaturation at 98 °C for 1 min; 30 cycles at 98 °C for 10 s, 50 °C for 30 s, 72 °C for 30 s, and finally, at 72 °C for 5 min. Primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′) were used for fungal PCR amplification of the ITS1 region of the genome. The reaction conditions were consistent with those used for the bacteria. The PCR amplicons were electrophoresed on a 2% agarose gel at 100–120 V for 20 min, and the target fragments were recovered using the AxyPrep PCR Purification Kit (E.Z.N.A.® Soil DNA Kit). The purified PCR products were quantified with the Quant-iT PicoGreen dsDNA Assay Kit on a Qubit fluorometric quantification system. Subsequently, library construction was performed. PCR products were quantified using the QuantiFluor™-ST Fluorometer (Promega Biotech, Beijing, China), and the samples were adjusted as needed for sequencing. Sequencing was conducted by Shanghai Majorbio Bio-pharm Technology (Shanghai, China), using an Illumina MiSeq platform (San Diego, CA, USA).

Raw sequencing data were analyzed using the QIIME2 (v2021.11) platform, and each sample data was split from the downstream data according to the barcode and primer sequences using FLASH (version 1.2.11) and Cutadapt software. The data were quality-controlled using fastp software (version 0.23.1), and tag data were obtained from the quality-controlled data. Data quality control was performed using Fastp software (version 0.23.1) to obtain high-quality tagged data. Chimeric sequences were removed by comparison with the Silva and Unite species annotation databases, and effective tags were obtained (Edgar et al., 2011). The DADA2 module (Wang et al., 2021) was used to reduce the noise of the effective tags to obtain Amplicon Sequence Variants (ASVs) and characterization tables. Species annotation was performed using QIIME2, with Silva 138.1 for the bacterial species annotation database and Unite v9.0 for the fungal species annotation database. Rapid multiple sequence comparisons were performed using QIIME2 to determine the phylogenetic relationships of all ASV sequences. Finally, the least amount of data in each sample was used as the criterion for homogenization, and homogenized data were obtained.

2.5 Physical and chemical properties of soil and litter quality analysis

Kay type nitrogen determination method, acid-soluble molybdenum antimony colorimetry, and flame photometry were used, respectively, to measure N, P, K content of soil. The electrode potential method was used to determine the soil pH and potassium dichromate oxidation-Spectrophotometry was used to determine the total SOC. The collected litter was brought to the laboratory, dried in an oven at 85 °C until a constant mass was achieved, and then weighed. All the above-mentioned indicators were measured using methods established in previous literature (Liu et al., 2021).

2.6 Data processing and statistical analysis

Species abundance histograms of relative abundance in Perl were plotted using the SVG function; species abundance clustering heatmap was plotted using the pheatmap function in R language; Venn diagrams and flower diagrams were generated by the VennDiagram function in R language and SVG function in Perl language, respectively; and the chao1, observed_features, Pielou _e, Shannon, and Simpson diversity indices were generated using QIIME2 software. One-way analysis of variance was performed using SPSS software (version 26.0) to detect the differences between diversity indices (p < 0.05); Non-metric multidimensional (NMDS) were performed using R software with the ade4 and ggplot2 packages. Based on the selected data set, Adonis test between different groups was calculated using the psych package (Revelle, 2017). Effect size measurements (LEfSe) software was used to calculate and draw the linear discriminant analysis (LDA) value distribution bar chart and the evolution branch chart (Segata et al., 2011). The functional potential of bacterial communities across different samples was predicted using PICRUSt2 (V2.3.0) based on ASV abundance and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, whereas the ecological functions of fungal communities were predicted using the FunGuild database. Spearman correlation indices were calculated for the target samples, and filtering conditions were set to (1) remove connections with correlation coefficients <0.6, (2) filter out node self-connections, (3) remove connections with node abundance less than 0.005%, and (4) remove connections with p-values of correlation coefficients >0.05. Finally, graphviz-2.38.0 software was used to visualize the microbial community network diagram.

The soil physical and chemical properties and litter mass data were processed using Microsoft Excel 2016. Duncan’s multiple range test was used to analyze the significance of differences between groups, with a significance level of p < 0.05. The Mantel test and Pearson’s correlation analysis were performed using the dplyr, ggcor, and ggplot2 packages in R.

3 Results

3.1 Analysis of soil physical and chemical properties in Daxinganling larch under different harvesting strategies

Soil physical and chemical properties and litter quality are important indicators for evaluating forest soil health. As shown in Figure 1, the total N and SOC content in SC sample plots were significantly higher than those in other sample plots (p < 0.05); the total P content in OCL sample plots was significantly higher than that in other sample plots (p < 0.05), and the total P content decreased with increasing OC intensity; the total K content in OCH sample plots was significantly higher than in other sample plots (p < 0.05), with total K content increasing with increasing OC intensity; the pH of the PF sample plots was significantly higher than SC, OCL, OCH, and CC sample plots (p < 0.05). An analysis of litterfall mass in different harvesting sample plots showed that, as depicted in Figure 1, litterfall mass in PF sample plots was significantly higher than that in other sample plots (p < 0.05), and litterfall mass increased with increasing OC intensity. The litterfall mass in the CC sample plots was significantly lower than that in the other sample plots (p < 0.05).

FIGURE 1
“Bar graphs showing soil properties across different land types: - Nitrogen (N) is highest in SC, lowest in OCH. - Phosphorus (P) peaks in OCL, least in CC. - Potassium (K) is most in OCH, lowest in CC. - Soil Organic Carbon (SOC) is highest in SC, lowest in CC. - pH is similar in OCL and OCH, lowest in CC. - Litter is most in PF, least in CC. Categories include PF, SC, OCL, OCM, OCH, CC.”

Figure 1. Analysis of soil physical and chemical properties and litter mass in plots with different harvesting strategies. PF, primary forest plots (no harvesting); SC, shelterwood cutting plot; OCL, optional cutting with low intensity plots; OCM, optional cutting with moderate intensity plots; OCH, optional cutting with high intensity plots; CC, clear cutting plots.

3.2 Soil microbial ASV counts in Daxinganling larch under different harvesting strategies

Soil bacterial and fungal ASV counts were analyzed in the sample plots of the different harvesting strategies (Figure 2). At the bacterial level, the number of ASVs was 1717 in the PF sample plots, 1375 in the SC sample plots, 1131 in the OCL sample plots, 885 in the OCM sample plots, 595 in the OCH sample plots, and 939 in the CC sample plots; 554 ASVs were noted in sample plots with different harvesting strategies; and 888 ASVs were detected in the OCL, OCM and OCH sample plots. At the fungal level, the number of ASVs was 822 in the PF sample plots, 849 in the SC sample plots, 281 in the OCL sample plots, 428 in the OCM sample plots, 602 in the OCH sample plots, and 436 in the CC sample plots. In the sample plots with different harvesting strategies, 56 ASVs were recorded, while the total number of ASVs in the OCL, OCM, and OCH sample plots was 140. Across all sample plots, the number of bacterial ASVs was higher than that of fungal ASVs, indicating a greater bacterial abundance than that of fungi.

FIGURE 2
Two diagrams labeled A and B illustrate overlapping circles. A contains two Venn diagrams: one with overlapping circles labeled OCL, OCM, and OCH, showing intersections with numbers; another flower-shaped diagram includes SC, OCL, OCM, OCH, PF, and CC, with a core labeled 554. B mirrors A’s layout, with different intersection numbers and a core labeled 56. A color key identifies circle labels.

Figure 2. Analysis of soil bacterial and fungal Amplicon Sequence Variants (ASV) counts in sample plots with different harvesting strategies. (A) ASV counts of soil bacteria in different harvesting practices. (B) ASV counts of soil fungi in different harvesting practices. The CORE numbers in the middle of the petal diagram represent the number of feature sequences common to all sample plots, and the numbers on the petals represent the number of feature sequences specific to the sample plots. The numbers in the overlapping part of the circles in the Wayne diagram represent the number of feature sequences shared among the sample plots, and the numbers without overlapping part represent the number of feature sequences specific to the sample plots.

3.3 Soil microbial relative abundances in Daxinganling larch under different harvesting strategies

Based on the results of species annotation at the bacterial and fungal phylum levels, the top 10 phyla in terms of relative abundance in each plot were selected, and the remaining phyla were grouped as “Others” to create a bar chart of the relative abundance of species.

At the phylum level (Figure 3), the dominant phyla in the bacterial communities were similar among the sample plots with different logging methods; however, their relative abundances differed. Proteobacteria, Verrucomicrobiota, and Acidobacteria were highly abundant as co-dominant phyla in all the sample plots. Proteobacteria had the highest relative abundance in all sample plots, Acidobacteria had the highest relative abundance in OCM sample plots, and Verrucomicrobiota had the second highest relative abundance in all sample plots after Proteobacteria and Aidobacteria. Desulfobacterota and Methylomirabilota had lower relative abundances in all sample plots; Desulfobacterota had the highest relative abundance in the PF sample plots and the lowest in CC sample plots, and Methylomirabilota had the highest relative abundance in PF sample plots and the lowest in CC sample plots. The relative abundance of Proteobacteria tended to increase with increasing intensity of OC.

FIGURE 3
Stacked bar chart displaying the relative abundance of various bacterial phyla across six categories: PF, SC, OCL, OCM, OCH, and CC. Each bar is segmented by color-coded phyla including Proteobacteria, Verrucomicrobiota, Acidobacteriota, and others, as indicated in the legend on the right. The vertical axis represents relative abundance from 0 to 1.

Figure 3. Relative abundance of soil bacteria at the gate level in sample plots with different harvesting practices.

At the fungal phylum level (Figure 4), the dominant phyla of the fungal communities in the different harvesting method sample plots were the same, Ascomycota and Basidiomycota, but their relative abundance varied among the sample plots. The relative abundance of Ascomycota was highest in the CC, OCL, OCM, and OCH sample plots, and that of Basidiomycota was highest in the PF and SC. The sum of the relative abundances of Ascomycota and Basidiomycota was greater than 0.93 in the PF, SC, OCL, OCM, and OCH sample plots, and was 0.73 in the CC sample plots. Mortierellomycota was second only to Ascomycota and Basidiomycota in terms of relative abundance across sample plots, with the highest relative abundance in the CC sample plots and the lowest in the OLM. The relative abundances of fungal phyla, such as Rozellomycota, Mucoromycota, and Chytridiomycota, were low in all sample areas. The relative abundance of Ascomycota decreased with increasing harvesting intensity in the OC, whereas that of Basidiomycota increased with increasing harvesting intensity.

FIGURE 4
Bar chart showing the relative abundance of various fungal phyla across six different sample groups labeled PF, SC, OCL, OCM, OCH, and CC. Ascomycota (red) and Basidiomycota (blue) are predominant in all groups. Other phyla such as Glomeromycota, Olpidiomycota, and Aphelidiomycota are represented in smaller proportions, indicated by various colors in the legend on the right.

Figure 4. Relative abundance of soil fungi at the phylum level in sample plots with different harvesting practices.

To further investigate the similarities and differences in the composition of bacterial and fungal communities in the sample plots under the different harvesting strategies, the top 35 genera ranked by relative abundance were selected for species abundance clustering analysis based on species annotation and abundance information at the genus level. As shown in Figure 5, the relative abundance of bacteria at the genus level in each sample plot was analyzed using a clustering analysis. The relative abundances of Rhodoplanes and Terrimonas in the PF sample plots were higher than those in the other sample plots. In the CC sample plots, the relative abundance of Candidatus_Udaeobacter was higher than that in the other samples, whereas the relative abundance of unidentified_TK10 was lower than that in the other samples. In the SC sample plots, the relative abundances of Janthinobacterium, unidentified_Mitochondria, and Pseudomonas were higher than those in the other samples, whereas the relative abundance of Candidatus_Solibacter was lower than that in the other samples. The relative abundances of Haliangium and Mycobacterium were higher than those in the other samples in the OCL sample plots. The relative abundance of unidentified_Subgroup_2 was higher than that in the other samples, whereas the relative abundances of Gemmatimonas and Chthoniobacter were lower than those in the other samples in the OCM sample plots. The relative abundances of Massilia, Ellin6067, and IS-44 were higher than those of the other samples, and the relative abundances of Haliangium and Candidatus_Koribacter were lower than those of the other samples in the OCH sample plots.

FIGURE 5
Heatmap showing bacterial taxa abundance across different categories labeled PF, SC, OCL, OCM, OCH, and CC. Colors represent abundance from low (blue) to high (red). A side legend identifies phyla with colors: Acidobacteriota, Actinobacteriota, Bacteroidota, Chloroflexi, Desulfobacterota, Gemmatimonadota, Myxococcota, Proteobacteria, and Verrucomicrobiota. Taxa are listed on the right.

Figure 5. Clustering diagram of species abundance at the bacterial genus level in the sample plots of different harvesting strategies. The horizontal coordinates in the figure are the sample plots names and the vertical coordinates are the species annotation information; the clustering tree on the left side of the figure is the species clustering tree; the values corresponding to the heatmap are the Z-values obtained by standardizing the relative abundance of species in each row.

The relative abundance of fungi at the genus level in each sample plot was analyzed using clustering analysis (Figure 6). The relative abundance of Microbotryales_gen_Incertae_sedis, Saitozyma and Myxozyma in the PF sample plots was higher than that of the other sample plots, while the relative abundance of Serendipita was lower than that of the other sample plots. The relative abundances of Mortierella, Podila, and Inocybe were higher in the CC sample plots than in the other samples. In the SC sample plots, the relative abundances of Penicillium, Leucosporidium, Tricholoma, Coniochaeta, Neurospora, Paracladophialophora, Hymenogaster, Sebacina and Paratritirachium were higher than those in the other samples. The relative abundances of GS_ord_Incertae_sedis_gen_Iencertae_sedis and Tomentella were higher than those in the other OCL sample plots, whereas the relative abundance of Piloderma was lower than that of the other samples. The relative abundances of Suillus, Archaeorhizomyces, and Archaeorhizomyces were higher in the OCM sample plots, whereas the relative abundances of Ceminibasidum and Lachnellula were lower than those in the other samples. The relative abundances of Hebeloma, Gymnopus, and Clavulinopsis were higher than those in OCH sample plots.

FIGURE 6
Heatmap depicting the abundance of different fungal genera across various samples labeled PF, SC, OCL, OCM, OCH, and CC. Color intensity ranges from blue (lower abundance) to red (higher abundance), indicating the abundance level. The right color key and phylum labels show Ascomycota (red), Basidiomycota (green), and Mortierellomycota (blue). Fungal genera listed on the right include Serendipita, Thelephora, and others, with hierarchical clustering dendrograms on the left.

Figure 6. Clustering of species abundance at the fungal genus level in the sample plots of different harvesting strategies.

3.4 Alpha diversity of soil microbial communities in Daxinganling larch under different harvesting strategies

Based on the annotated ASV data, the alpha diversity of the soil bacterial communities in the sample plots of the different harvesting strategies was analyzed (Table 2). The highest index values (chao1, observed_features, Pielou, Shannon, and Simpson) were noted in the PF sample plots, which were higher than those in the other sample plots (p < 0.05) and significantly higher than those in the OCH sample plots. The index values showed a decreasing trend with increasing OC intensity. The PF sample plots had the most abundant and uniformly distributed bacterial species and the highest alpha diversity, whereas OCH significantly reduced the bacterial alpha diversity (p < 0.05); its alpha diversity was the lowest, and the intensity of OC was negatively correlated with the alpha diversity of the soil bacterial community.

TABLE 2
www.frontiersin.org

Table 2. Alpha diversity index of soil bacterial communities in sample plots of different harvesting strategies.

Based on the annotated ASV data, the alpha diversity of the soil fungal communities in the sample plots under different harvesting strategies was analyzed (Table 3). The SC sample plots had the highest values for all indices (chao1, observed_features, Pielou, Shannon, and Simpson), which were higher than those of the other sample plots. The OCL sample plots had significantly lower values of chao1 and observed_features than the SC sample plots, and the OCM sample plots had the lowest Pielou, Shannon, and Simpson indices. As OC intensity increased, the values of chao1 and observed_features tended to increase, and the values of Pielou, Shannon, and Simpson tended to decrease and then increase again. This indicates that the fungal species in the SC sample plots were the most abundant and evenly distributed, and the alpha diversity was the highest.

TABLE 3
www.frontiersin.org

Table 3. Alpha diversity index of soil fungal communities in sample plots of different harvesting strategies.

3.5 Beta diversity analysis of soil microbial communities in Daxinganling larch under different harvesting strategies

To study the similarity of soil bacterial and fungal community structures in the sample plots of the different harvesting strategies, NMDS and Adonis analyses were performed. As shown in Figures 7A, B, the stress was less than 0.2, indicating that NMDS could accurately reflect the degree of difference in the bacterial and fungal community structures under different harvesting strategies. As shown in Figure 7A, the different OC intensity, CC, and PF sample plots were far apart, and R2 was greater than 0.4 (Table 4), whereas the SC and PF sample plots were closer to each other, and R2 was 0.24, indicating that the bacterial community structure of the different intensity of OC and CC sample plots was more different from that of PF sample plots, whereas the SC and PF sample plots were similar in terms of their bacterial community structure. As shown in Table 4, the R2 of each harvested sample plot was in the range of 0.3–0.4 compared with the PF sample plots, and the CC and OCM were farther away from the PF sample plots, whereas the other sample plots were closer to the PF sample plots (Figure 7B). This indicates that the fungal community structures of the CC and OCM sample plots differed greatly from those of the PF sample plots, whereas the fungal community structures of the other sample plots were similar to those of the PF sample plots. As shown in Table 4, the Pr values of the differences between the bacterial and fungal groups were greater than 0.05, indicating that the harvesting strategies had no significant effect on the beta diversity of either bacteria or fungi.

FIGURE 7
Two NMDS plots labeled A and B display different colored and shaped data points representing groups: SC (red squares), QCL (purple circles), PF (blue circles), CC (orange diamonds), OCM (aqua triangles), and OCH (green X’s). Plot A has a stress value of 0.12, and Plot B has a stress value of 0.14. The axes are labeled NMDS1 and NMDS2.

Figure 7. NMDS analysis of bacterial and fungal community structure in sample plots of different harvesting strategies. (A) NMDS analysis of bacterial community structure in the sample plots with different harvesting strategies. (B) NMDS analysis of fungal community structure in the sample plots of different harvesting strategies. Each point in the graph indicates a sample, the distance between points indicates the degree of difference, and sample plots of the same group are indicated by the same color. When stress is less than 0.2, NMDS can accurately reflect the degree of difference between sample plots.

TABLE 4
www.frontiersin.org

Table 4. Bacterial and fungal community structure of sample plots with different harvesting strategies Adonis intergroup difference analysis.

3.6 Analysis of soil bacterial and fungal biomarkers in Daxinganling larch under different harvesting strategies

LEfSe multilevel species difference discriminant analysis was used to identify biomarkers associated with differences in the structure of bacterial and fungal communities across sample plots. At the bacterial level, LEfSe analyses of the bacterial communities in sample plots with different OC intensities (OCL, OCM, and OCH) did not reveal biomarkers with significantly different abundances. LEfSe analysis of bacterial communities in all sample plots (Figure 8) showed that A21b was significantly enriched in the OCH sample plots, Methylomirabilia, Rokubacteriales, SC_l_84, and Methylomirabilota were significantly enriched in the PF sample plots, whereas the remaining sample plots with other harvesting strategies were not significantly enriched for the biomarkers (Figure 8A). The absolute LDA values of the five biomarkers were all greater than four, indicating highly significant differences (p < 0.01), which may be the key biomarkers leading to differences in the bacterial community structure among the PF, CC, and SC, and different intensities of OC sample plots. Cluster analysis of the evolutionary relationships of the above five biomarkers (Figure 8B) showed that the biomarkers in the PF sample plots were categorized into two evolutionary branches (Methylomirabilota and SC_l_84). Methylomirabilia and Rokubacteriales were the class and order belonging to Methylomirabilota, respectively, whereas SC_l_84 formed an independent evolutionary branch. A21b formed an independent evolutionary branch in the OCH sample plots.

FIGURE 8
Bar chart and cladogram illustrate microbial communities. In part A, the bar chart shows LDA scores for microbial taxa in two groups: OCH (red) and PF (green). Taxa such as f_A21b (red) and others (green) are differentiated. In part B, the cladogram graphically represents the phylogenetic relationships of these taxa, with colored nodes corresponding to the same groups: OCH and PF.

Figure 8. LEfSe analysis of soil bacterial communities in sample plots with different harvesting strategies. (A) Histogram of the distribution of bacterial LDA values. (B) Branching diagram of bacterial evolution. Species with LDA score greater than a set value (default setting of 4) are shown in the histogram of LDA value distribution.

At the fungal level, LEfSe analyses of the fungal communities were performed in sample plots with different OC intensities (OCL, OCM, and OCH) (Figure 9). The results showed that Hyaloscyphaceae and Helotiales were significantly enriched in the OCL sample plots, Omphalotaceae, Gymnopus, and Gymnopus_junquilleus were significantly enriched in the OCH sample plots, and no biomarkers were significantly enriched in the OCM sample plots (Figure 9A). The absolute LDA values of the five biomarkers were all greater than four, indicating highly significant differences (p < 0.01), which may be the key biomarkers leading to differences in the fungal community structure among the OCL, OCM, and OCH sample plots. Cluster analysis of the evolutionary relationships of the five biomarkers showed that two in the OCL sample plots were within the evolutionary branch of Helotiales and three in the OCH sample plots were within the evolutionary branch of Omphalotaceae (Figure 9B).

FIGURE 9
Graph A shows an LDA score chart comparing OCL (red) and OCH (green) groups with labeled taxa. Graph B displays a cladogram highlighting specific taxa with OCL in red and OCH in green.

Figure 9. LEfSe analysis of soil fungal communities in sample plots with different selection intensities. (A) Histogram of the distribution of fungal LDA values. (B) Branching diagram of fungal evolution.

LEfSe analysis of the soil fungal communities in the sample plots for each harvesting strategy (Figure 10) showed that Sordariomycetes and Coniochaeta were significantly enriched in the SC sample plots. Archaeorhizomycetales, Archaeorhizomycetaceae, Archaeorhizomyces_sp, Archaeor- hizomyces and Archaeorhizomycetes were significantly enriched in OCM sample plots. Russula_cascadensis and Saitozyma_podzolica were significantly enriched in the PF sample plots, whereas the other harvesting strategy sample plots were not significantly enriched in biomarkers (Figure 10A). The absolute LDA values of the nine biomarkers were all greater than four, indicating highly significant differences (p < 0.01), which may be the key biomarkers leading to differences in the fungal community structure among the PF, CC, and OC sample plots. Cluster analysis of the evolutionary relationships of the nine biomarkers showed that five of them in the OCM sample plots were within the evolutionary branch of Archaeorhizomycetes, two biomarkers in the SC sample plots were within the evolutionary branch of Sordariomycetes, and two biomarkers in the PF sample plots belonged to two separate evolutionary branches (Figure 10B).

FIGURE 10
Panel A shows a horizontal bar chart representing LDA scores (log 10) for different taxonomic groups with OCM (red), PF (green), and SC (blue) categories. Panel B is a cladogram illustrating the phylogenetic relationships among the same groups, with highlighted sections in red, green, and blue corresponding to the categories from the bar chart.

Figure 10. LEfSe analysis of soil fungal communities in sample plots with different harvesting strategies. (A) Histogram of the distribution of fungal LDA values. (B) Branching diagram of fungal evolution.

3.7 Correlation between key microbial markers and physical and chemical properties of the soil

To further clarify the relationship between key microbial markers and soil physicochemical indicators in Xing’an larch forests under different logging methods, the Mantel Test was used for analysis. From a bacterial perspective (Figure 11A), total soil N was significantly correlated with the bacterial family SC_l_84 (p < 0.05). SOC was significantly correlated with the bacterial family SC_l_84 (p < 0.01). From a fungal perspective (Figure 11B), total soil N was significantly correlated with the fungal class Sordariomycetes and the fungal genus Coniochaeta (p < 0.01). Soil total P was significantly correlated with the fungal class Archaeorhizomycetes and the fungal order Archaeorhizomycetales (p < 0.05), and was significantly correlated with the fungal family Archaeorhizomycetaceae, fungal species Archaeorhizomyces_sp, and fungal genus Archaeorhizomyces (p < 0.01). SOC was highly significantly correlated with the fungal class Archaeorhizomycetes and the fungal genus Coniochaeta (p < 0.01).

FIGURE 11
Side-by-side correlograms labeled A and B. Both display relationships among variables like N, P, K, SOC, and pH, with colored tiles indicating correlation and Mantel’s p-value levels. Diagram A focuses on classes like Methylomirabilota and Rokubacteria, while diagram B includes groups such as Archaeorhizomycetes and Sordariomycetes. Lines connect variables, indicating interactions.

Figure 11. Correlation analysis between key microbial markers in soil and soil physical and chemical properties. (A) Bacteria. (B) Fungi.

3.8 Predictive analysis of soil microbial functions in Daxinganling larch under different harvesting strategies

Based on the prediction of bacterial gene function using PICRUSt2, the top 35 functions in terms of abundance and abundance information in each sample plots were selected to draw heat maps for comparison with the KEGG database. As shown in Figure 12, the relative abundance of gene functions of ABC-2 type transport system permease protein and elongation factor G was higher in the PF sample plots than in the other sample plots; the relative abundance of gene functions of the putative ABC transport system ATP-binding protein in the SC sample plots were higher than that of other sample plots; the relative abundance of acetyl-CoA C-acetyltransferase, peptide/nickel transport system substrate-binding protein, and long-chain- C-acetyltransferase, peptide/nickel transport system substrate-binding protein, and long-chain- fatty-acid-CoA ligase in the OCL sample plots; in OCM sample plots, the relative abundance of gene function of periplasmic protein TonB, basic amino acid/polyamine antiporter, APA family, serine/threonine protein kinase, bacterial, putative ABC transport system permease protein, and thioredoxin reductase were higher than in the other sample plots; and in the OCH sample plots, the relative abundance of branched- chain amino acid transport system substrate-binding protein, glutathione S-transferase, iron complex outer membrane receptor protein, methyl-accepting chemotaxis protein were higher. The relative abundance of the gene functions of the biopolymer transport protein ExbD and the uncharacterized protein was higher in the CC sample plots than in the other sample plots. Overall, compared to PF, SC, OC, and CC induced the enrichment of different functional bacteria to a certain extent to adapt to the effects of harvesting on soil microbiology in the understory.

FIGURE 12
Heatmap displaying the expression levels of various proteins across different conditions, labeled PF, SC, OCL, OCM, OCH, and CC. Protein names are listed on the right, and color gradients from blue to red indicate expression levels, with blue for low and red for high. Color bar on the right shows values ranging from negative one point five to positive one point five. Dendrograms indicate clustering of similar expression patterns.

Figure 12. Predictive analysis of soil bacterial functions in sample plots with different harvesting strategies. The horizontal coordinate is the sample plots name, the vertical coordinate is the functional annotation information, and the clustering tree on the left side of the figure is the functional clustering tree; the value corresponding to the heat map is the Z-value obtained by normalizing the relative abundance of each row of function.

Fungal ecological functions were predicted based on FunGuild, and the top 35 functions in terms of abundance and their abundance in each sample plots were selected to draw heat maps. As shown in Figure 13, the fungal ecological functions enriched in the sample plots of different harvesting strategies were different, and there were multifunctional combinations of fungi, such as Animal_Parasite-Fungal_ Paraplot and Leaf_Saprotroph, indicating that the fungi possessed multiple ecological functions in the soil to cope with changes in the soil microenvironment. The relative abundance of eight functions, such as Animal_Parasite-Fungal_Paraplotand Leaf_Saprotroph in the PF sample plots was higher than that in the other sample plots; the relative abundance of seven functions, such as Dung_Saprotroph-Ectomycorrhizal and Lichenized in the SC sample plots was higher than other sample plots; the relative abundance of Bryophyte_Parasite-Undefined_Saprotroph in the OCL sample plots; the relative abundance of four functions, including Soil_Saprotroph and Lichenized_Paraplotin the OCM sample plots; and the relative abundance of the Undefined Saprotroph-Undefined_Saprotroph function and six other functions in the OCH sample plots were higher than other sample plots; the relative abundance of Undefined_Saprotroph-Undefined_Symbiotroph and Ericoid_Mycorrhizal in CC sample plots was higher than other sample plots. Different harvesting strategies induced different ecological functions in the enriched fungal community.

FIGURE 13
Heatmap illustrating fungal diversity across different sample types, labeled PF, SC, OCL, OCM, OCH, and CC. Color gradient from blue to red indicates data values ranging from -1.5 to 2. Various fungal categories are listed on the right, showing their relative abundance or presence in each sample type.

Figure 13. Predictive analysis of soil fungal function in sample plots with different harvesting practices.

3.9 Analysis of soil microbial network interactions in Daxinganling larch under different harvesting strategies

To investigate the covariance characteristics of soil microorganisms, the microbial community covariance network of the sample plots with different harvesting strategies was determined using a one-way covariance network analysis.

Analysis of the bacterial networks in the sample plots subjected to various harvesting strategies (Figure 14 and Table 5) revealed that the PF sample plots exhibited higher node and edge numbers, clustering coefficients, network path densities, and average degrees, along with shorter average path lengths. These findings suggest that the PF bacterial community demonstrated the greatest complexity, largest network size, highest connectivity, and most robust network stability. By analyzing sample plots with different OC intensities, the number of nodes, edges, clustering coefficient, network path density, and average degree of the bacterial network showed a decreasing trend, whereas the average path length showed an increasing trend. This indicates that OC can reduce the complexity and stability of soil bacterial community network. The number of nodes, edges, network path density, and average degree in the CC sample plots were higher than those in the OC sample plots, indicating that the bacterial network in the CC sample plots had more species associations, higher network complexity, and a more stable community.

FIGURE 14
Six network graphs display different 16S networks labeled as PF, SC, OCL, OCM, OCH, and CC. Each graph shows nodes and connecting lines in purple, blue, or red shades, representing relationships among microbial communities. Nodes are labeled by various phyla, with a legend indicating color-coded classifications such as Firmicutes and Proteobacteria. The visual complexity varies, with densely connected networks showing intricate patterns of connectivity.

Figure 14. Soil bacterial network analysis of sample plots with different harvesting strategies. Different nodes represent different genera, node size represents the average relative abundance of the genus, nodes of the same phylum are of the same color, the thickness of the connecting lines between the nodes is positively correlated with the absolute value of the correlation coefficient of the species interactions, and the color of the connecting lines corresponds to the positive and negative of the correlation (red is positively correlated, blue is negatively correlated).

TABLE 5
www.frontiersin.org

Table 5. Topological parameters of soil bacterial and fungal co-occurrence network in sample plots of different harvesting strategies.

Fungal networks in sample plots with different harvesting strategies were analyzed (Figure 15 and Table 5). Within the PF sample plots, the fungal network exhibited a greater number of nodes and edges, along with higher clustering coefficients, network path densities, and average degrees than those in the other sample plots. Conversely, it demonstrated a shorter average path length. These characteristics indicate that the fungal community within the PF sample plots possessed the highest complexity, largest network dimensions, optimal connectivity, and most stable network. By analyzing the sample plots with different OC intensities, the number of nodes, edges, clustering coefficient, network path density, and average degree of the fungal network showed an increasing trend with increasing intensity of OC. This indicates that OC can improve the complexity and stability of the soil fungal community network.

FIGURE 15
Six interconnected network graphs titled “ITS Network: RF,” “ITS Network: SC,” “ITS Network: DCL,” “ITS Network: DCM,” “ITS Network: OCH,” and “ITS Network: CC” showing nodes and colorful lines indicating relationships. Nodes are color-coded by phylum: Ascomycota in red, Basidiomycota in green, Chytridiomycota in blue, and Olpidiomycota as Fungi_SPP_unclass_zotab. Lines between nodes vary in color intensity, suggesting differing connection strengths. Each graph represents complex relationships, with distinct clusters interconnected through various pathways.

Figure 15. Soil fungal network analysis of sample plots with different harvesting strategies.

4 Discussion

4.1 Composition of soil microbial communities in sample plots with different harvesting strategies varies

In this study, six harvesting strategies were used to treat Daxinganling larch, and the soil microbial diversity in the sample plots of each harvesting method was analyzed after an 18-year recovery period. This study showed that in all sample plots, the number of bacterial ASVs was higher than that of fungi (Figure 2). This indicates that bacterial abundance was higher than that of fungi in the inter-root soil of Daxinganling larch, which may be because the root system of Daxinganling larch has a stronger symbiotic interaction with bacteria than with fungi and a higher demand for bacterial populations (Santhanam et al., 2015).

The dominant bacterial phyla in the soil of Daxinganling larch under different harvesting strategies were Proteobacteria, Verrucomicrobiota, and Acidobacteriota. The results of the thinning showed that Proteobacteria, Actinobacteria, and Acidobacteria were the dominant phyla in both the thinning and control sample plots, followed by Verrucomicrobiota (Li et al., 2021), which is consistent with the results of our study. Several studies on soil bacterial diversity in Daxinganling larch have shown that Proteobacteria, Verrucomicrobiota, Acidobacteria, and Verrucomicrobiota are the dominant phyla (Ping et al., 2019; Jia et al., 2020). Our research findings confirmed this result. The relative abundance of Proteobacteria showed an increasing trend with the increase of OC intensity, which might be due to the fact that with the increase of OC, the content of litter left over from felling increases (Figure 1), and the eutrophic phylum Proteobacteria reproduces faster under the high nutrient environmental conditions (Fierer et al., 2007). Ascomycota and Basidiomycota were the dominant fungal phyla in the soils of Daxinganling larches under different harvesting strategies. Ascomycota and Basidiomycota were the dominant fungal phyla in the thinning sample plots of the Daxinganling larch forest (Li et al., 2021), which were the dominant fungal phyla in different forest types of Daxinganling larch (Yang et al., 2017). This indicates that Ascomycota and Basidiomycota, as eutrophic and oligotrophic fungal phyla, respectively, were the primary fungal phyla in Daxinganling larch, and their dominant ecological niches may not have been affected by harvesting strategies.

4.2 Diversity of soil microbial communities in Daxinganling larch is closely related to the harvesting strategies

Comparative analysis of the alpha diversity of bacterial communities in the sample plots with different harvesting strategies showed that the values of all alpha diversity indices of bacterial communities in the PF sample plots were significantly higher than those in the OCH sample plots (p < 0.05) and higher than those in the other sample plots. This indicates that different harvesting strategies can reduce the number of soil bacterial species. A previous study showed that the Shannon and Simpson values of soil bacterial communities in Daxinganling larch decreased after reasonable thinning, but no significant difference was observed compared with those before the intervention (Li et al., 2021). This suggests that moderate harvesting had little effect on the alpha diversity of soil bacteria in Daxinganling larch, whereas OCH significantly reduced bacterial alpha diversity (p < 0.05). In this study, the alpha diversity index values of the CC sample plots were higher than those of the OCH sample plots, which may be due to the large number of Daxinganling larch were planted after harvesting (Table 1), and the number of bacterial species was higher than that of the OCH sample plots after a recovery period of 18 years. A comparative analysis of the alpha diversity of the fungal communities in the different harvested sample plots showed that the alpha diversity indices of the soil fungal communities in the SC sample plots were higher than those in the PF sample plots. This suggests that SC can effectively increase the number of fungal species in the soil. This may be because SC is less damaging to the soil environment of Daxinganling larch than other harvesting strategies, and it can improve the diversity of understory plants by improving understory temperature, humidity, and light (Rottstock et al., 2014), which in turn promotes an increase in the number of soil fungal species and numbers (Qu et al., 2024).

Analysis of the similarity of the bacterial and fungal community structures in the sample plots of different harvesting strategies showed that the bacterial community structures in the SC sample plots did not differ significantly from those in the PF sample plots, consistent with a previous study (Yang et al., 2020). The bacterial community structure in the other sample plots differed from that in the PF sample plots, but the degree of difference decreased with an increase in OC intensity (Table 4), the reasons for which need further investigation. Our study results showed that the fungal community structure of all harvested sample plots differed from that of the PF sample plots, and the degree of difference decreased with an increase in OC intensity, indicating that the response mechanism of the fungal community structure to harvesting was more complex and sensitive. However, the specific reasons for this need to be elucidated in further studies (Zhou et al., 2020).

4.3 Biomarkers associated with differences in bacterial and fungal community structure between sample plots with different harvesting strategies

After analyzing the biomarkers related to the differences in the bacterial community structure of different harvesting strategies, it was observed that the bacterial genus SC_l_84, which has a denitrification effect, was significantly enriched in the PF sample plots. After long-term natural succession, the humus layer in the PF sample plots is relatively deep, and it is easy to form an anoxic environment suitable for heterotrophic anaerobic bacterial SC_l_84, and the continuous large amount of apoplastic input provides a stable carbon source (Holden and Treseder, 2013). However, when PF was harvested, the anaerobic environment of the humus layer was destroyed, the litter was reduced (Hartmann et al., 2009), and the amount of SC_l_84 decreased significantly. We analyzed the biomarkers associated with differences in fungal community structure in the sample plots under different harvesting practices. The results showed that the five OCM biomarkers belong to Archaeorhizomycetes, which are often distributed in plant roots and have saprophytic functions but do not form mycorrhizal structures with roots (Rosling et al., 2011) and may form a mutualistic network with tufted mycorrhizal Glomerales to affect plant growth (Choma et al., 2016). The biomarker Russula_cascadensis, an ectomycorrhizal fungus at the PF sample plots, can promote plant root growth and nutrient uptake rate (Yang et al., 2025), and the mechanism of its effect on the fungal community structure in the unharvested sample plots requires further exploration.

4.4 Effects of different harvesting strategies on soil microbial networks in Daxinganling larch

Microbial co-occurrence network analysis can reveal the interactions between microorganisms, which, in turn, reflect the complexity and stability of the community (De Vries et al., 2018; Ramirez et al., 2018). In our study, the bacterial and fungal community networks were the most complex and stable in the PF sample plots, and the complexity and stability of the bacterial community in the CC sample plots were higher than those in the OC. Several studies on Hainan forests in China have shown that the bacterial and fungal community networks in OC sample plots are more structured than those in PF sample plots, whereas the network structure of CC sample plots is poorer (Chen et al., 2019; Yu et al., 2023). This finding differs from that of our study. First, this is because of the different study areas and forest types. Additionally, all harvested sample plots in this study were artificially planted with a large amount of Daxinganling larch after CC, which recovered the soil bacterial community faster than the OC. In this study, OC reduced the network complexity and stability of soil bacterial communities; however, the opposite was observed for fungal communities. This may be because OCH accelerates the formation of forest windows, promotes apoplastic inputs and nutrient cycling (Figure 1), provides ecological niches for rare fungal taxa, drives the shift from competition to mutualistic symbiosis, and enhances the stability of fungal community networks (Qiang et al., 2023).

5 Conclusion

This study provides new insights into the molecular response mechanisms of different harvesting strategies that affect microbial diversity. The highest bacterial alpha diversity was observed in the PF sample plots, whereas the highest fungal alpha diversity was observed in SC sample plots. OCH significantly reduced bacterial alpha diversity (p < 0.05), and a negative correlation was observed between OC intensity and bacterial alpha diversity. Harvesting strategies had no significant effect on bacterial or fungal beta diversity. Soil nitrogen, phosphorus, and other physical and chemical properties were significantly correlated with different microbial markers. The bacterial and fungal community networks in PF sample plots exhibited the highest complexity and stability. OC reduced the complexity and stability of soil bacterial networks, while it had the opposite effect on fungal communities, enhancing their network complexity and stability. This study provides a preliminary analysis of how different harvesting strategies influence soil microbial diversity in L. gmelinii forests of the Greater Khingan Range, offering practical insights for the functional restoration and conservation of this forest ecosystem.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: http://www.ncbi.nlm.nih.gov/bioproject/PRJNA1347439/.

Author contributions

YW: Conceptualization, Writing – original draft, Writing – review & editing. NT: Data curation, Investigation, Project administration, Writing – review & editing. HZ: Data curation, Investigation, Writing – review & editing. ML: Investigation, Writing – review & editing. SL: Investigation, Writing – review & editing. JG: Investigation, Writing – review & editing. YS: Investigation, Writing – review & editing. RS: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Natural Science Foundation of China (Grant number:32460388) and the Natural Science Foundation of Inner Mongolia Autonomous Region (Grant number:2023MS03051).

Acknowledgments

We are grateful to the reviewers for their valuable comments that helped improve the manuscript’s quality. Thanks to the Natural Science Foundation of China (Grant number: 32460388) and the Natural Science Foundation of Inner Mongolia Autonomous Region (Grant number: 2023MS03051) for funding this study.

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.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Abbreviations

SC, shelterwood cutting; OC, optional cutting; CC, clear cutting; OCH, optional cutting with high intensity; PF, primary forests; OCL, optional cutting of low intensity; OCM, optional cutting of moderate intensity; SOC, soil organic carbon; ASVs, Amplicon Sequence Variants.

References

Bont, L. G., Fraefel, M., Frutig, F., Holm, S., Ginzler, C., and Fischer, C. (2022). Improving forest management by implementing best suitable timber harvesting strategies. J. Environ. Manag. 302:114099. doi: 10.1016/j.jenvman.2021.114099

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Xu, H., He, D., Li, Y. D., Luo, T. S., Yang, H. G., et al. (2019). Historical harvesting alters soil fungal community composition and network in a tropical rainforest. Forest Ecol. Manag. 433, 228–239. doi: 10.1016/j.foreco.2018.11.005

Crossref Full Text | Google Scholar

Chen, X. L., Wang, D., Chen, X., Wang, J., Diao, J. J., Zhang, J. Y., et al. (2015). Soil microbial functional diversity and biomass as affected by different thinning intensities in a Chinese fir plantation. Appl. Soil Ecol. 92, 35–44. doi: 10.1016/j.apsoil.2015.01.018

Crossref Full Text | Google Scholar

Cheng, Z. C., Wu, S., Du, J., Liu, Y. Z., Sui, X., and Yang, L. B. (2023). Reduced arbuscular mycorrhizal fungi (AMF) diversity in light and moderate fire plots in Taiga forests, Northeast China. Microorganisms 11:1836. doi: 10.3390/microorganisms11071836

PubMed Abstract | Crossref Full Text | Google Scholar

Choma, M., Bárta, J., Šantrüčková, H., and Urich, T. (2016). Low abundance of Archaeorhizomycetes among fungi in soil metatranscriptomes. Sci. Rep. 6:38455. doi: 10.1038/srep38455

PubMed Abstract | Crossref Full Text | Google Scholar

Cordovez, V., Rotoni, C., Dini-Andreote, F., Oyserman, B., Carrión, V. J., and Raaijmakers, J. M. (2021). Successive plant growth amplifies genotype-specific assembly of the tomato rhizosphere microbiome. Sci. Total Environ. 772:144825. doi: 10.1016/j.sc.10.1016/j.scitotenv.2020.144825

Crossref Full Text | Google Scholar

Dang, P., Gao, Y., Liu, J. L., Yu, S. C., and Zhao, Z. (2018). Effects of thinning intensity on understory vegetation and soil microbial communities of a mature Chinese pine plantation in the Loess Plateau. Sci. Total Environ. 630, 171–180. doi: 10.1016/j.scitotenv.2018.02.197

PubMed Abstract | Crossref Full Text | Google Scholar

De Vries, F. T., Griffiths, R. I., Bailey, M., Craig, H., Girlanda, M., Gweon, H. S., et al. (2018). Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 9:3033. doi: 10.1038/s41467-018-05516-7

PubMed Abstract | Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., Abades, S., Alfaro, F. D., et al. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat. Ecol. Evol. 4, 210–220. doi: 10.1038/s41559-019-1084-y

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | Crossref Full Text | Google Scholar

Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839

PubMed Abstract | Crossref Full Text | Google Scholar

Hartmann, M., Howes, C. G., VanInsberghe, D., Yu, H., Bachar, D., Christen, R., et al. (2012). Significant and persistent impact of timber harvesting on soil microbial communities in Northern coniferous forests. ISME J. 6, 2199–2218. doi: 10.1038/ismej.2012.84

PubMed Abstract | Crossref Full Text | Google Scholar

Hartmann, M., Lee, S., Hallam, S. J., and Mohn, W. W. (2009). Bacterial, archaeal and eukaryal community structures throughout soil horizons of harvested and naturally disturbed forest stands. Environ. Microbiol. 11, 3045–3062. doi: 10.1111/j.1462-2920.2009.02008.x

PubMed Abstract | Crossref Full Text | Google Scholar

Holden, S. R., and Treseder, K. K. (2013). A meta-analysis of soil microbial biomass responses to forest disturbances. Front. Microbiol. 4:163. doi: 10.3389/fmicb.2013.00163

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, Y., Zhang, X., Chen, H., Jiang, Y., and Zhang, J. (2024). Effects of forest age and season on soil microbial communities in Chinese fir plantations. Microbiol. Spectr. 12, e4075–e4023. doi: 10.1128/spectrum.04075-23

PubMed Abstract | Crossref Full Text | Google Scholar

Jia, D., Dong, X. B., and Li, Y. H. (2020). Effect of strip harvesting on bacterial diversity of forest soils in the Daxing’an Mountains. Soil Sci. Soc. Am. J. 84, 512–521. doi: 10.1002/saj2.20046

Crossref Full Text | Google Scholar

Jog, R., Pandya, M., Nareshkumar, G., and Rajkumar, S. (2014). Mechanism of phosphate solubilization and antifungal activity of Streptomyces spp. isolated from wheat roots and rhizosphere and their application in improving plant growth. Microbiology 160, 778–788. doi: 10.1099/mic.0.074146-0

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Yu, J. B., Zhang, X. Z., Liu, Y., and Yang, L. X. (2021). Differential responses of soil bacterial and fungal community to short-term crop tree management in a Larix gmelinii plantation. Forests 12:1411. doi: 10.3390/f12101411

Crossref Full Text | Google Scholar

Li, Y., Sun, L., Zhang, J., Zhao, H., Su, T., Li, W., et al. (2025). Tea plant/Ophiopogon japonicus intercropping drives the reshaping of soil microbial communities in terraced tea plantation’s micro-topographical units. Agriculture 15:1150. doi: 10.3390/agriculture15111150

Crossref Full Text | Google Scholar

Liu, J. L., Dang, P., Gao, Y., Zhu, H. L., Zhu, H. N., Zhao, F., et al. (2018a). Effects of tree species and soil properties on the composition and diversity of the soil bacterial community following afforestation. Forest Ecol. 427, 342–349. doi: 10.1016/j.foreco.2018.06.017

Crossref Full Text | Google Scholar

Liu, J. L., Ha, V. N., Shen, Z., Zhu, H. L., Zhao, F., and Zhao, Z. (2018b). Characteristics of bulk and rhizosphere soil microbial community in an ancient Platycladus orientalis forest. Appl. Soil Ecol. 132, 91–98. doi: 10.1016/j.apsoil.2018.08.014

Crossref Full Text | Google Scholar

Liu, S., Yin, H., Li, X., Li, X., Fan, C., Chen, G., et al. (2021). Short-term thinning influences the rhizosphere fungal community assembly of Pinus massoniana by altering the understory vegetation diversity. Front. Microbiol. 12:620309. doi: 10.3389/fmicb.2021.620309

PubMed Abstract | Crossref Full Text | Google Scholar

Melo, V. F., Barros, L. S., Silva, M. C. S., Veloso, T. G. R., Senwo, Z. N., Matos, K. S., et al. (2021). Soil bacterial diversities and response to deforestation, land use and burning in North Amazon, Brazil. Appl. Soil Ecol. 158:103775. doi: 10.1016/j.apsoil.2020.103775

Crossref Full Text | Google Scholar

Morrison, E. W., Frey, S. D., Sadowsky, J. J., van Diepen, L. T., Thomas, W. K., and Pringle, A. (2016). Chronic nitrogen additions fundamentally restructure the soil fungal community in a temperate forest. Fungal Ecol. 23, 48–57. doi: 10.1016/j.funeco.2016.05.011

Crossref Full Text | Google Scholar

Mushinski, R. M., Gentry, T. J., and Boutton, T. W. (2018). Organic matter removal associated with forest harvest leads to decade scale alterations in soil fungal communities and functional guilds. Soil Biol. Biochem. 127, 127–136. doi: 10.1016/j.soilbio.2018.09.019

Crossref Full Text | Google Scholar

Nandal, A., Yadav, S. S., Rao, A. S., Meena, R. S., and Lal, R. (2023). Advance methodological approaches for carbon stock estimation in forest ecosystems. Environ. Monit. 195:315. doi: 10.1007/s10661-022-10898-9

PubMed Abstract | Crossref Full Text | Google Scholar

Ping, L., Rongjiu, S., Feng, Z., Hua, Y. J., Yang, C. X., Gui, H. J., et al. (2019). Soil bacterial community structure and predicted functions in the larch forest during succession at the Greater Khingan Mountains of Northeast China. Chin. J. Appl. Ecol. 30, 95–107. doi: 10.13287/j.1001-9332.201901.010

PubMed Abstract | Crossref Full Text | Google Scholar

Qiang, W., Gunina, A., Kuzyakov, Y., Luo, R., Zhang, Y., Liu, B., et al. (2023). Shifts of understory vegetation induced by thinning drive the expansion of soil rare fungi. J. Environ. Manag. 342:118119. doi: 10.1016/j.jenvman.2023.118119

PubMed Abstract | Crossref Full Text | Google Scholar

Qu, X. J., Li, X. G., Bardgett, R. D., Kuzyakov, Y., Revillini, D., Sonne, C., et al. (2024). Deforestation impacts soil biodiversity and ecosystem services worldwide. Proc. Natl. Acad. Sci. U.S.A. 121:e2318475121. doi: 10.1073/pnas.2318475121

PubMed Abstract | Crossref Full Text | Google Scholar

Ramirez, K. S., Geisen, S., Morriën, E., Snoek, B. L., and van der Putten, W. H. (2018). Network analyses can advance above-belowground ecology. Trends Plant Sci. 23, 759–768. doi: 10.1016/j.tplants.2018.06.009

PubMed Abstract | Crossref Full Text | Google Scholar

Repo, A., Albrich, K., Jantunen, A., Aalto, J., Lehtonen, I., and Honkaniemi, J. (2024). Contrasting forest management strategies: Impacts on biodiversity and ecosystem services under changing climate and disturbance regimes. J. Environ. Manag. 371:123124. doi: 10.1016/j.jenvman.2024.123124

PubMed Abstract | Crossref Full Text | Google Scholar

Revelle, W. R. (2017). psych: Procedures for personality and psychological research. Evanston, IL: Northwestern University.

Google Scholar

Rosling, A., Cox, F., Cruz-Martinez, K., Ihrmark, K., Grelet, G.-A., Lindahl, B. D., et al. (2011). Archaeorhizomycetes: Unearthing an ancient class of ubiquitous soil fungi. Science 333, 876–879. doi: 10.1126/science.1206958

PubMed Abstract | Crossref Full Text | Google Scholar

Rottstock, T., Joshi, J., Kummer, V., and Fischer, M. (2014). Higher plant diversity promotes higher diversity of fungal pathogens, while it decreases pathogen infection per plant. Ecology 95, 1907–1917. doi: 10.1890/13-2317.1

PubMed Abstract | Crossref Full Text | Google Scholar

Santhanam, R., Luu, V. T., Weinhold, A., Goldberg, J., Oh, Y., and Baldwin, I. T. (2015). Native root-associated bacteria rescue a plant from a sudden-wilt disease that emerged during continuous cropping. Proc. Natl. Acad. Sci. U.S.A. 112, E5013–E5020. doi: 10.1073/pnas.1505765112

PubMed Abstract | Crossref Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, H. R., Du, X. R., Zhang, Z. Y., Feng, F. J., and Zhang, J. M. (2023). Rhizosphere interface microbiome reassembly by arbuscular mycorrhizal fungi weakens cadmium migration dynamics. Imeta 2:e133. doi: 10.1002/imt2.133

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y. Y., Guo, H., Gao, X. G., and Wang, J. H. (2021). The intratumor microbiota signatures associate with subtype, tumor stage, and survival status of esophageal carcinoma. Front. Oncol. 11:754788. doi: 10.3389/fonc.2021.754788

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, J. Y., Karunarathna, S. C., Patabendige, N., Tarafder, E., Lou, D. J., Zhou, Y. Q., et al. (2025). Unveiling the bioactive compounds and therapeutic potential of Russula: A comprehensive review. J. Fungi 11:341. doi: 10.3390/jof11050341

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, K., Zhu, J. J., Xu, S., and Zheng, X. (2018). Conversion from temperate secondary forests into plantations (Larix spp.): Impact on belowground carbon and nutrient pools in northeastern China. Land Degrad. Dev. 29, 4129–4139. doi: 10.1002/ldr.3169

Crossref Full Text | Google Scholar

Yang, K., Zhu, J. J., Yan, Q. L., and Sun, O. J. (2010). Changes in soil P chemistry as affected by conversion of natural secondary forests to larch plantations. Forest Ecol. Manag. 260, 422–428. doi: 10.1016/j.foreco.2010.04.038

Crossref Full Text | Google Scholar

Yang, L. B., Sui, X., Zhu, D. G., Cui, F. X., Li, J. B., Song, R. Q., et al. (2017). Study on fungal communities characteristics of different Larix gmelini forest typesin cold temperate zone. Cent. South Univ. For. Technol. 37, 76–84. doi: 10.14067/j.cnki.1673-923x.2017.12.013

Crossref Full Text | Google Scholar

Yang, Y., Qiu, Y. M., Wang, Z. B., and Qu, L. Y. (2020). Effect of different harvesting strategies on physicochemical properties and microbial community of Larix gmelinii rhizosphere soil. Acta Ecol. Sin. 40, 7621–7629. doi: 10.5846/stxb202001070053

Crossref Full Text | Google Scholar

Yu, J. J., Cong, W., Ding, Y., Jin, L. X., and Zhang, Y. G. (2023). Natural restoration characteristics and assembly mechanisms of soil microbial community in tropical rainforest under different disturbance types. Chin. J. Ecol. 42:534. doi: 10.13292/j.1000-4890.202303.007

Crossref Full Text | Google Scholar

Zhai, K. T., Hua, Y. C., Liang, J. W., Li, J., Wang, Z., Liu, L., et al. (2023). Soil microbial diversity under different types of interference in birch secondary forest in the Greater Khingan Mountains in China. Front. Microbiol. 14:1267746. doi: 10.3389/fmicb.2023.1267746

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, L., Cai, L., He, Z., Wang, R., Wu, P., and Ma, X. (2016). Thinning increases understory diversity and biomass, and improves soil properties without decreasing growth of Chinese fir in southern China. Environ. Sci. Pollut. Res. 23, 24135–24150. doi: 10.1007/s11356-016-7624-y

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, Z. H., Wang, C. K., Ren, C. J., and Sun, Z. H. (2020). Effects of thinning on soil saprotrophic and ectomycorrhizal fungi in a Korean larch plantation. Forest Ecol. Manag. 461:117920. doi: 10.1016/j.foreco.2020.117920

Crossref Full Text | Google Scholar

Zhu, W. J., Zhang, H., and Wang, Y. X. (2021). Research progress on effects of harvesting on forest soil respiration. J. Zhejiang A F Univers. 38, 1000–1011. doi: 10.11833/j.issn.2095-0756.20210365

Crossref Full Text | Google Scholar

Keywords: harvesting, Larix gmelinii, Greater Khingan Mountains/Daxinganling, microbial diversity, functional prediction

Citation: Wang Y, Ta N, Zhang H, Li M, Liu S, Gong J, Song Y and Sa R (2025) Beyond one-size-fits-all approach: How do various harvesting strategies shape soil microbial diversity in Larix gmelinii (Daxinganling larch) forests of the Greater Khingan Mountains? Front. Microbiol. 16:1654005. doi: 10.3389/fmicb.2025.1654005

Received: 25 June 2025; Accepted: 06 October 2025;
Published: 29 October 2025.

Edited by:

Weronika Goraj, The John Paul II Catholic University of Lublin, Poland

Reviewed by:

Sunil Kumar Behera, Berhampur University, India
Sunanda Mishra, Odisha University of Agriculture & Technology, India
Ananya Kuanar, Siksha ‘O’ Anusandhan University, India

Copyright © 2025 Wang, Ta, Zhang, Li, Liu, Gong, Song and Sa. 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: Rula Sa, c2FydWxhMjEzQDE2My5jb20=

Disclaimer: 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.