Huadong Xu · Wenjing Chen · Yanjun Cao ·Lihai Wang

Abstract The decay rate of standing Korean pine ( Pinus koraiensis) in natural forests can be as high as 50% and is likely inf luenced by the soil properties and nutrient and water status of the site. To clarify the relationship between the severity of tree decay and soil properties in order to prevent decay in a natural mixed forest in the Xiaoxing’an Mountains, wood strength of standing trees was nondestructively assessed, and the severity of decay of extracted wood cores was quantif ied based on diff erences in mass between two decayed increment cores extracted at breast height and an intact increment core near the decayed ones. Soil samples from the critical root zone (non-rhizosphere) of each tree were analyzed for chemical properties and microbial composition. The abundance of chemical elements (especially total N and K) and the species richness of soil microbes increased as decay severity increased. Fungal number (FN)and actinomycetes number (AN) were related to decay severity ( R 2 = 0.504). Bacterial number (BN) was higher than FN or AN, but had a minor Effect on tree decay. Path analysis showed BN might indirectly inhibit decay by aff ecting FN.Decay severity was not signif icantly correlated with either soil fungal or bacterial diversity. These results suggest that forest managers need to monitor levels of fungi and total N and total K levels to reduce the decay of Korean pine.

Keywords Pinus koraiensis · Soil chemical properties ·Soil microbial properties · Trunk decay degree

Introduction

Korean pine (Pinus koraiensisSiebold & Zucc.), one of the dominant tree species in natural forests of the Xiaoxing’an Mountains in North China, is irreplaceable for providing structural stability to local forest ecosystems. Its wood is prized for its high quality and strength and easy processability (Lv and Qu 2010). However, wood decay ofPinus koraiensisin natural forests has become a serious problem(Li et al. 2014; Wang et al. 2015). Tree decay is an essential biological process that produces carbon dioxide and water with a release of energy to maintain forest processes (Bednarz et al. 2013), but wood quality is diminished by wood rot fungi during the decay process (Shortle and Dudzik 2012). Such decay aff ects tree health and timber production by degrading wood cell walls and reducing wood mechanical and physical properties. Factors such as forest density,site conditions (e.g., elevation, slope, humidity, temperature,light, and soil), hydrology and human and animal activities(Chen 1959) need to be assessed for their contributions to tree decay and loss so that appropriate management strategies can be designed.

Among the site conditions, soil properties are critical to the availability and absorption of water and nutrients by trees and aff ect the rhizosphere of standing trees and thus inf luence tree decay, as shown by our previous studies on soil physicochemical properties in relation to tree decay (Sun et al 2015). We found that enhanced C/N ratios and pH in the soil of the root zone was associated with more decay in Korean pine trunks. Hietala et al. ( 2016) noted that the incidence of decay in Norway spruce also increased as the pH increased from 3.8 to 6.4. These studies indicate a close relationship between tree decay and soil chemical properties.

These properties in turn are closely related to soil microbial communities (Hu et al. 2006; Martiny et al. 2006; Pietri and Brookes 2008); the abundance and composition of the communities can inf luence soil chemical properties, which can then aff ect chemical composition (Cai and Huang 2016).In addition, soil microbes can infect roots which could eventually cause trunk decay. However, little information is available on direct relationships between tree decay and soil microbial properties and how soil microbial quantity and diversity contribute to the incidence and severity of decay.Such information is critical to control tree decay; when the relationships are clearer, management practices can be designed to adjust soil properties to improve forest soil and reduce decay in Korean pine.

Therefore, the objective of this study was to analyze the relationship between the degree of tree decay and soil microbial properties based on our previous research and the Effects of soil chemical and microbial properties on decay severity in standing trees ofPinus koraiensis. This basic data will lead to the development of forest management practices that reduce decay in this important species.

Materials and methods

Study area

The study area in northeastern China is located in the Liangshui National Nature Reserve on the southern slope of the Dailing region, the largest branch of the Xiaoxing’an Mountains. It has a cold continental monsoon climate with an annual average temperature of 1.4 °C. Field tests and surveys ofPinus koraiensiswere carried out in a 30-ha sample plot in the 18th block of the Liangshui Forest District. The main soil type is dark brown forest soil; the dominant tree species isPinus koraiensisaccompanied by species such asUlmus pumilaTrautv.,Abies fabri, Betula costata Trautv., andTilia tuanSzyszyl.

Determination of tree decay severity

In the 18th block, the sample plot was divided into 750 quadrats, each 20 m × 20 m, and 100 quadrats containingP.koraensiswere then selected along an S-shaped transect in the plot. For ensuring the selected tree types were distributed over a wider area, no more than three trees were tested in each quadrat. A total of 200 trees were visually observed for decay symptoms and the internal physical structure of the wood tested using 2D stress wave tomography (Arbotom,Rinntech, Heidelberg, Germany) and a Resistograph (model 4453, Rinntech) according to the manufacturer’s instructions. Based on the nondestructive test results, nine trees with decay and one healthyP. koraensistrees were chosen as sample trees, each in a diff erent quadrat. For each sample tree, two decaying increment cores were extracted at breast height, and an intact core extracted near the decaying ones.Each wood core was oven-dried (105 °C) to a constant mass and weighed. The mass of the decayed wood cores (m1) and intact wood cores (m2) were used to calculate the loss in mass of wood cores (S) asS= [(m2−m1)/m2] × 100 (Liu et al.( 2015) to describe the severity of tree decay.

Soil chemical properties measurement

Soil samples of 500 g were collected from four sampling sites in the critical root zone (rhizosphere) of the selected trees. Surface vegetation, branches and leaves were removed before drilling with an auger to obtain samples from the A horizon at 6-20 cm depth. The 44 soil samples were placed in a sterile bag, which was sealed and taken immediately to the laboratory. Half of the samples were used for chemical analyses, and the other half for microbe isolations and stored at − 20 °C until used.

For chemical analyses, soil samples were air-dried,screened with a 2 mm mesh, and f ine roots removed. Total and hydrolyzed N, total and available P, total and available K contents were determined according to forest industry standards in China (LY). Total N was measured using the Kjeldahl method (LY/T 1228-1999). Total K was tested using acid dissolution-f lame photometry (LY/T 1234-1999). Total P was measured using anti-colorimetric determination of acid-soluble molybdenum and antimony (LY/T1232-1999).Hydrolyzed N was measured using alkali hydrolysis diff usion (LY/T 1231-1999). Available K was extracted using ammonium acetate-f lame photometry (LY/T 1236-1999).Available P was tested using the anti-colorimetric method and molybdenum and antimony with sodium hydroxide extraction (LY/T 1233-1999).

Soil microbial isolation and quantif ication

From each fresh soil sample, 1 g was diluted to 10 -2 − 10 -7 with distilled water and shaken 15-20 min. Each dish of medium received 0.05 mL of the 10 -5 to 10 -7 dilutions to plate bacteria on beef extract peptone, 10 -3 to 10 -5 for fungi on Martin’s medium and 10 -2 to 10 -4 for actinomycetes on Gaoshi No. 1 medium (Table 1; Xu et al. 2019). Each concentration was plated out on three dishes.

Table 1 Composition of media to isolate bacteria, fungi and actinomycetes from soil

The plates were inverted and incubated at 37 °C. Colonies were counted using an automatic colony counter (Interscience scan300) after two days for bacteria, three days for fungi and seven days for actinomycetes. The number of colonies per gram of fresh soil was then calculated.

Fungal and bacterial diversity estimation

DNA was extracted from soil microorganisms for PCR amplification of the internal transcribed spacer (ITS)sequence to estimate fungal diversity and the 16S rDNA sequence to estimate bacterial diversity.

Extraction of DNA

The Fast DNA Spin Kit for Soil Kit (MP Biomedicals Biomedical Co., Irvine, CA, USA) was used to extract total DNA from 500 mg in a Lysing Matrix E tube with 978 ml sodium phosphate buff er and 122 ml MT buff er.The tube was placed in the Fastprep instrument and rotated at 6.0 m s −1 speed for 40 s, centrifuged for 10 min at 14,000×g. After centrifugation, the supernatant was transferred to 250 μL [protein precipitation solution (PPS)] in a 2.0 mL catch tube which was centrifuged for 5 min at 14,000×gto precipitate the protein. The supernatant was transferred to a 15 ml clean test tube with 1 ml binding matrix solution and rotated for 2 min, and then precipitated for 3 min to obtain the silica gel matrix. After removing the supernatant and adding SEWS-M and air-drying, DNA was extracted and stored in a refrigerator at − 20 °C.

PCR amplif ication

The extracted DNA was separated electrophoretically in 1% agarose to detect whether the extracted DNA was in accordance with the sequencing requirements. Before amplif ication of the target genes, a random bacterial and fungal DNA sample was selected to test whether the target gene was amplif ied and with a minimum number of cycles.

For bacteria, the 16S V3-V4 region was amplified using primer pair 338F (ACT CCT ACG GGA GGC AGC AG)/806R (GGA CTA CHVGGG TWT CTAAT), Trans-Start Fastpfu DNA polymerase (Beijing TransGen Biological Co., Beijing) with 30 ng DNA sample, 1 μL forward primer (5 μM), 1 μL reverse primer (5 μM), 3 μL BSA (2 ng μL −1 ), 12.5 μL 2 × Taq PCR MasterMix and 7.5 μL ddH2O. Cycling conditions were 95 °C for 5 min;28 cycles of 95 °C for 45 s, 55 °C for 45 s and extended at 72 °C for 45 s; 10 min at 72 °C. PCR amplif ication products (target genes) were preserved at 4 °C and sequenced as soon as possible.

For fungi, the ITS1-ITS2 region was amplif ied using primer pair ITS1-F (CTT GGT CAT TTA GAG GAA GTAA)/ITS2 (TGC GTT CTT CAT CGA TGC) and the same polymerase and PCR system used for bacteria, except 32 cycles were used.

DNA high-throughput sequencing

After electrophoretic detection and quality control, the amplif ied fragments were used to construct a Miseq library using the NEBNext Ultr DNA Library Prep Kit for Illumina and the standard protocol. The amplif ied fragment and the expanded sublibrary were then sequenced using the Illumina MiSeq PE300 platform (Illumina, San Diego, CA, USA).The PE reads obtained were spliced according to the overlap relationship between PE reads, and the sequence quality was controlled and f iltered. According to the barcode and primer information at both ends of the sequence, the Effective sequence was distinguished, the sequence direction was corrected, and the optimized sequence was obtained by using software FLASH and Trimmomatic. These steps were carried out by the Allwegene Technology Company in Beijing.

Statistical analyses

We used Pearson correlation analysis to assess correlations between the severity of decay and each of the soil property indicators, and a one-samplet-test to analyze diff erences between soil fungal diversity in the root zone of decaying and healthy sample trees. Multiple linear regression was used to construct the optimum regression equation relating the severity of tree decay with soil microbial quantity indicators, which were the dependent and independent variables, respectively. Path analysis was employed to compare the Effect of diff erent kinds of soil microorganisms on tree decay. All statistical analyses were done using SPSS 22.0(IBM Co., Armonk, NY, USA).

Results

Tree trunk decay severity

Decay changes wood properties and decreases its quality.Mass loss in the wood core is usually used as a true value of wood decay (Liu et al. 2015). The mass loss data for 10 sample trees is considered as decay severity (Table 2). The mass loss of the wood core C2 was 0 i.e., from a healthy tree.Cores R3 and S2 had the most severe decay with mass loss of 47.1% and 48.3%, respectively.

Soil chemical properties versus trunk decay severity of standing trees

As the degree of decay increased for the sample trees, total N and K in the root zone soil rhizosphere gradually increased(Fig. 1 a). Each was higher than the total P which changed slightly as decay severity increased. Similarly, hydrolytic N and available K were higher than available P in the root zone soil and showed a slightly increasing trend (Fig. 1 b).Thus, these elements were more abundant in the root soil around decaying standing trees than around healthy trees,and might be due to the their release from decaying litter and their incorporation into soil organic matter. In addition, there was more decaying litter under decaying trees. Reduced uptake of nutrients from the soil might also contribute to the increase in these elements.

Fig. 1 Relation between tree decay severity and soil chemical properties

Table 2 Decay severity as percentage loss in mass in 10 wood cores from sample trees

Soil microbial quantity versus trunk decay severity of standing trees

The distribution of soil microorganisms around the roots of sample trees with diff erent decay severities was: bacteria > actinomycetes > fungi, with the number of bacteria(BN) 10 − 10 4 times higher than actinomycetes, and 10 3 -10 5 times higher than fungi (Fig. 2). Thus, bacteria were clearly the key contributors to soil microbial activity. As the severity of decay increased, the number of microbes of all microbial types increased gradually (Fig. 2), indicating decay severity might be closely related to the community of root soil microorganisms.

Multiple linear regression analysis was used to analyze the relationship between decay severity and root soil microbial abundance. The regression result isy= − 18.457 + 3.564 × 10-4 × 2 + 7.683 × 10-5 × 3 (R2 = 0.664,p= 0.000),where y is tree decay severity (Wood core WLR), × 2 is the number of fungi (FN) and × 3 is the number of actinomycetes (AN). The equation shows that decay severity was mainly related to the number of FN and AN in the soil.However, a low determination coeffi cient (R2 = 0.664) and a large residual factor [e = 0.5(1 −R) = 0.336] were also found,which indicates that other factors might exist in addition to FN and AN.

Fig. 2 Relation between tree decay severity and soil microbial quantity

In the path analysis to further study the Effect of soil microorganisms on tree decay, the direct path coeffi cient(DPC,Piy) between any independent variablexiand dependent variableyis the standard coeffi cient of the linear regression equation, and the indirect path coeffi cient (IPC,Pxi-xj)is the product of the Person correlation coeffi cient (r ij) and DPC (P jy) of any independent variablexi toy. The DPCs of BN, FN and AN on decay severity was − 0.305, 0.397 and 0.413, respectively, indicating that the direct Effect of bacteria on decay severity was negative and not signif icant(p= 0.113), but the direct Effects of fungi and actinomycetes on decay severity were positive.Based on the Pearson correlation coeffi cient, the IPCs between soil microbe numbers and decay severity were calculated. The greatest negative IPC, − 0.209, was found between BN and FN, indicating that BN could aff ect FN in the soil and then might indirectly inhibit tree decay. All other IPCs between soil microbes were positive; the greatest positive IPC, 0.211, was between AN and FN, which meant that soil AN could accelerate tree decay by aff ecting FN. The IPC between AN and BN was the smallest (0.117).

Soil microbial diversity versus trunk decay severity of standing trees

Soil microbial species diversity

α-Diversity represents the richness and evenness of diff erent species of microorganisms. PD_whole_tree curve and Rank_abundance curve were applied to characterize the α-diversity of species (Figs. 3, 4).

Fig. 3 Rarefaction measure PD_whole_tree curves for a fungi and b bacteria

PD_whole_tree is the sum of the branching lengths in the phylogenetic trees constructed using all groups in the sample (Goedert et al. 2015). The higher the PD_whole_tree value, the greater the species richness in the sample. For both fungal and bacteria species, F3 was the most abundant (fungi: PD_whole_tree value > 175; bacteria: > 80) and D3 was the lowest in species number and diversity.

For the Rank_abundance curve to analyze species diversity, the shape reflects the homogeneity of species composition. The slower the downward trend of the curve, the higher the uniformity of the species composition (Bates et al. 2013). For fungi and bacteria, the curve decreased rapidly and steeply, indicating that the dominant flora in the samples accounted for a higher proportion and the species distribution was not uniform (Fig. 4).

Statistical analyses for relation between soil microbial diversity and tree decay

α-Diversity represents the richness of organisms in a particular region or ecosystem. The commonly used indices are the Chao1 index, PD_whole_tree index and Shannon index. In the Pearson correlation analysis of the relationship between decay severity ofP. koraiensisand soil microbial α-diversity, there was no signif icant correlation between fungal α-diversity and decay severity. The same results were found for bacterial α-diversity (Table 3).

One-samplet-test was employed to analyze the relation between fungal PD_whole_tree index and bacterial PD_whole_tree index; the results show that there was no signif icant diff erence between decay and healthy tree with a value in fungal diversity (2-tailed,p= 0.062) and in bacterial diversity (2-tailed,p= 0.195), shown in Table 4.

Fig. 4 Rank_abundance curves of fungi and bacteria

Table 3 Pearson correlation analysis of microbial diversity and decay severity of sample trees

Table 4 One-sample t-test of PD_whole_tree for bacteria

Fig. 5 The unifrac_distance graph for species in various fungal phyla, Cercozoa, and others

Identif ication and analysis of dominant microbial species

Species richness, the number of individuals of each species in the samples, is obtained from the operational taxonomic units (OTUs), while the abundance of OTUs can be calculated from the number of sequences in OTU. Through high-throughput sequencing, seven known phyla and 203 known species were found in soil samples from the rhizosphere around the roots. Among all species, 21 had a relative abundance > 1%. The top three fungal phyla were Zygomycota (6.3%), Ascomycota (35.4%) and Basidiomycota(55.4%), accounting for 97.1% of all fungi (Fig. 5). The relative abundance of agaricomycetes (Basidiomycota), the primary group associated with wood decay (Shang 2008),was 46.6%, but further statistical analysis showed that their abundance was not correlated with decay severity of standing trees (R= -0.229,p= 0.525). For other species, some are known soil inhabitants (Cryptococcus_podzolicus; Liu 2016), others ectomycorrhizal fungi (Inocybe_maculata;Larsson et al. 2009), and other pathogenic fungi (Nectria_ramulariae). But their distribution in 10 samples was not correlated with decay severity.

Table 5 Fungal species with ＞ 1% abundance and their relative abundance in 10 core samples (C2 is from a healthy tree)

The relationship between relative abundance of the 203 species and the decay severity ofP. koraiensiswas further analyzed. Of all species, the number with more than 1%of relative abundance was 21, includingAmanita_argentea,Tricholoma_terreum,Russula_viscida,Russula_aff_integra_r_03014. Their relative abundance in 10 samples was shown in Table 5.

In similar analyses for bacteria, species with relative abundance ＞ 1% included Nitrospirae (1.78%), Bacteroidetes(4.32%), Gemmatimonadetes (4.34%), Chlorof lexi (5.16%),Verrucomicrobia (7.36%), Actinobacteria (10.82%), Acidobacteria (30.03%) and Proteobacteria (33.34%), accounting for 97.15% of all bacteria. The relative abundance of diff erent bacterial species was not correlated with various decay severities.

Conclusions

Higher levels of N, P, and K, especially total N and K,were found in the rhizosphere from decaying standing trees than from healthy standing trees.

As tree decay severity increased, the quantity of microorganisms in the soil around the tree roots increased gradually. Statistical analysis showed that decay severity was related to the number of fungi and actinomycetes. Soil bacteria may indirectly inhibit tree decay by aff ecting the number of fungi in the soil.

High-throughput sequencing and statistical analysis showed that decay severity ofPinus koraiensiswas not correlated with soil microbial (fungal and bacterial)diversity.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Publisher’s NoteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affi liations.