Rhizocompartmental microbiomes of arrow bamboo (Fargesia nitida) and their relation to soil properties in Subalpine Coniferous Forests

Arrow bamboo (Fargesia nitida) is a pioneer plant in secondary forest succession in the Sichuan Province mountains. To comprehensively investigate the microbial communities and their functional variations in different rhizocompartments (root endosphere, rhizosphere, and root zone) of arrow bamboo (Fargesia nitida), a high-throughput metagenomic study was conducted in the present study. The results showed that the abundances of the dominant bacterial phyla Proteobacteria and Actinobacteria in the bamboo root endosphere were significantly lower than those in the rhizosphere and root zones. In contrast, the dominant fungal phyla, Ascomycota and Basidiomycota, showed the opposite tendency. Lower microbial diversity, different taxonomic composition and functional profiles, and a greater abundance of genes involved in nitrogen fixation (nifB), cellulose degradation (beta-glucosidase), and cellobiose transport (cellulose 1, 4-beta-cellobiosidase) were found in the bamboo root endosphere than in the other rhizocompartments. Greater soil total carbon, total nitrogen, NH4+-N, microbial biomass carbon, and greater activities of invertase and urease were found in the bamboo root zone than in the adjacent soil (spruce root zone). In contrast, the soil microbial community and functional profiles were similar. At the phylum level, invertase was significantly related to 31 microbial taxa, and the effect of NH4+-N on the microbial community composition was greater than that of NO3−-N. The soil physicochemical properties and enzyme activities were significantly correlated with microbial function. These results indicate that the root endosphere microbiomes of arrow bamboo were strongly selected by the host plant, which caused changes in the soil nutrient properties in the subalpine coniferous forest.


INTRODUCTION
Plant-associated microbiomes have an obvious influence on the nutrition, growth, and health of plants, the structure of plant communities, and the resistance and adaptation of plants to environmental changes (Ke, Miki & Ding, 2015).Plants can also regulate microbial activity and community composition in the rhizosphere by secreting root exudates or bioactive molecules (Dhungana, Kantar & Nguyen, 2023;Hu et al., 2018).Another type of root-associated microbiome is endophytes, which colonize the root endosphere to obtain a stable habitat and nutrient supply and exert more lasting effects than rhizospheric microbes (Hong et al., 2015).Endophytes have multiple functions in plant physiology, such as antibiotic activity, phytohormone production, and immunity enhancement, which can promote plant growth and mitigate biotic and abiotic stresses in plants (Dini-Andreote, 2020;Jasim et al., 2014).
Previously, a hierarchical decrease in microbial diversity in rhizocompartments toward the roots has been described, suggesting that plants filter and recruit their microbiome subsets (Bulgarelli et al., 2012).A two-step model for microbiome selection by plants has been suggested (Barajas et al., 2020;Bulgarelli et al., 2013), which includes (1) chemotaxis of microbes from bulk soil to the rhizosphere caused by root exudates, which is the first step in selecting rhizosphere microbes from bulk soil; and (2) selective entry of microbes to the root endosphere regulated by plant genetic factors, such as host plant defense strategies, root structure, and root exudation (Bulgarelli et al., 2012), which further select endosphere microbes from the rhizosphere community.According to this model, the composition of microbial communities in the root endosphere, rhizosphere, and bulk soil could be substantially different from each other, and the filtration effects followed the order of endosphere >rhizosphere >root zone (Xiao et al., 2017).
Arrow bamboo in the genus Fargesia, a semi-woody plant, is a rhizomatous perennial species that often dominates the plant community in the understory of some montane species in East Asia and South America (Liu et al., 2014;Saitoh, Seiwa & Nishiwaki, 2002).The dense leaves of arrow bamboo form a barrier for penetrating sunlight to the lower layers, which could inhibit the survival and regeneration of tree seedlings, saplings, and mature trees in many forests (Takahashi et al., 2003;Wang et al., 2007).Furthermore, Fargesia bamboo is the main food source of the endangered giant panda (Wang, Tao & Zhong, 2009).In southwest China, arrow bamboos are widely distributed understory plants in the subalpine spruce forest, a primary representative subalpine coniferous forest after the clearing of old-age dark coniferous forests since the 1950s (Ma et al., 2007b;Zhang, Liu & Gu, 2011).Arrow bamboo is a pioneer plant in various secondary forest communities in natural succession after local coniferous forests are destroyed (Ma et al., 2007a;Xu et al., 2016).In addition to its negative effects on the diversity of understory species and forest regeneration caused by its high density (Wang, Shi & Tao, 2012), arrow bamboo might also change soil physicochemical traits, as reported for other bamboo species (Kaushal et al., 2020;Umemura & Takenaka, 2015), which might affect the soil microbial communities and subsequent plants.
Reportedly, the bacterial microbiomes of moso bamboo (Phyllostachys edulis) are structurally variable in both the rhizosphere and the endosphere in Guangxi Province, China (Yuan et al., 2021).However, the filtration (selection) of arrow bamboo rhizocompartments and the interactions of the bamboo microbiome with soil characteristics in subalpine coniferous forests have not been described.Therefore, we conducted the present study to (i) investigate the diversity and function of microbiomes in the root zone, rhizosphere, and endosphere of arrow bamboo in the subalpine spruce forest to reveal the filtration function of rhizocompartments and (ii) analyze the linkage between soil properties and microbial communities.

Study site
The experiment was conducted at the Maoxian Mountain Ecosystem Research Station of the Chinese Academy of Sciences ( 103• 54 E, 31 • 42 N) on the eastern edge of the Tibetan Plateau of China.Located in a montane temperate climate region, the sampling site presented a mean annual temperature of 9.3 • C, annual precipitation of 850 mm, mean annual evaporation of 796 mm, and Calcic Luvisol soil (Food and Agriculture Organization of the United Nations (FAO), 2006).Local spruce (Picea asperata Mast.) plantations were established in the 1980s as a typical cultivated tree species in the eastern Tibetan Plateau subalpine region.Arrow bamboo, Fargesia nitida (Mitford) Keng f. ex T.P. Yi, is the dominant species scattered throughout the shrub layer of the spruce alpine forest.The other epidemic understory plants in this region are Asparagus filicinus D. Don, Phlomis umbrosa Turcz., Thladiantha davidii Franch., and Sinosenecio oldhamianus (Maxim.)B. Nord., Plantago major Linn., Carpesium divaricatum Sieb.& Zucc., Thalictrum uncinulatum Franch., Paraprenanthes melanantha (Franch) Ze H. Wang, and Rosa sericea Lindl.(Hu et al., 2016).The samples were collected from the Maoxian Mountain Ecosystem Research Station of the Chinese Academy of Sciences with field permit number 2019016.In this study, the treatments (samples) were (1) arrow bamboo root endosphere (BE), (2) arrow bamboo rhizosphere (BR), (3) arrow bamboo root zone (BZ), and (4) spruce root zone (SZ), which also serves as a control of the soil.Soils and roots were sampled at the end of the growing season on October 11, 2019.Nine replicate plots (400 m 2 each) were set-up in the alpine forest (mixed-coniferous forest).
From each plot, the root systems of five randomly selected F. nitida plants were removed.Because most of the roots were distributed in the 0-20-cm soil layer, the root zone soils were also sampled from this layer.After removing the soil particles loosely attached to the roots by vigorous shaking, the soil that tightly adhered to the root surface was brushed off as rhizosphere soil.At the same time, a composite soil sample for each plot was taken as the root zone soil (<20 cm from the collar of the sampled bamboo).Correspondingly, nine spruce root zone soils were collected from the same plots (approximately 20 m from the sampled bamboo to weaken the rhizosphere effect of the bamboo plant).Surface soils (0-20 cm in depth) were sampled from a 20 cm × 20 cm area with a distance of 0.5 m from the spruce tree.All root and soil samples were transported to the laboratory, and parts were stored at 4 • C and −20 • C after passing through a 2-mm sieve for further analysis.Another portion of each soil sample was air-dried for physicochemical analysis.

Analysis of soil physicochemical properties and enzyme activity
The soil samples were passed through sieves with a width of 0.15 mm.Total carbon (TC) and total nitrogen (TN) were determined using an elemental analyzer (Vario MACRO cube CN, Elementar Analysensysteme, Germany) (Wang et al., 2016).Soil NH 4 + -N and NO 3 − -N were extracted using 1 M KCl and measured by SEAL AA3 continuous flow analysis (SealAnalytical, Norderstedt, Germany).Available phosphorus (AP) was evaluated using the Olsen method with inductively coupled plasma combined with optical emission spectroscopy (Olsen & Cole, 1954).The standard method of Mc Lean & Watson (1985) was used to determine available potassium (AK) content.Soil electrical conductivity (EC) and pH were measured in 1:5 and 1:2.5 (w/w) soil: water suspensions, respectively.Soil microbial biomass carbon (MBC) and nitrogen (MBN) were measured using the chloroform fumigation extraction method (Vance, Brookes & Jenkinson, 1987), using Kec and Ken factors of 0.45 and 0.54 for extractions of carbon and nitrogen, respectively.
Invertase activity was measured according to the method of Xiao et al. (2017), based on the absorbance measurement at 508 nm of glucose released, which was reacted with 3,5-dini-trosalicylic acid.To assay the urease activity, 5 g of air-dried soil was incubated with 10 mL of 10% (w/v) urea solution, one mL of toluene, and 20 mL of citrate buffer (pH 6.7) at 37 • C for 24 h.Urease activity was expressed as the ammonium released in mg g −1 soil at 24 h (Xiao et al., 2017).Protease activity was measured using the method described by Ladd and Butler (1972).Protease activity was expressed as the amount of ammonium released (µg g −1 soil h −1 ).Cellulase activity was expressed as milligrams of glucose equivalents per gram of soil per hour (Guan, Zhang & Zhang, 1986).

Sample preparation and metagenomic DNA extraction
Surface sterilization and decontamination of surface DNA were performed on root samples according to a previously described protocol (Woźniak et al., 2019;Yang et al., 2001).The endophytic microbiome DNA (0.5 g of each root sample) was extracted using an EZNA™ plant DNA kit (Omega Bio-Tek, Norcross, GA, USA).Metagenomic DNA (0.5 g) from each rhizosphere and root zone soil sample was extracted using the EZNA™ Omega Mag-bind soil DNA kit in triplicate (Omega Bio-Tek, Norcross, GA, USA), according to the manufacturer's instructions.The quantity and quality of the DNA extracts were evaluated using a NanoDrop spectrophotometer (ND-1000; Thermo Fisher Scientific, Waltham, MA, USA) and 1% (w/v) agarose gel electrophoresis, respectively.

Sequencing and bioinformatic analyses
Initially, 36 metagenomic DNA libraries (four treatments, nine duplicates) for shotgun sequencing were prepared using Illumina TruSeq Nano library preparation kits with insert sizes of 400 bp.Each library was sequenced on an Illumina NovaSeq platform (Illumina, Foster City, CA, USA) using the PE150 strategy at Personal Biotechnology Co., Ltd.(Shanghai, China).Raw sequencing data were processed to obtain quality-filtered reads as follows: adapter sequences were trimmed using Cutadapt (v1.2.1) (Martin, 2011); (2) sequences with low quality (quality value <20 or contain N bases) and short reads (<50 bp) were removed; (3) Chloroplast DNA and mitochondrial DNA sequences were removed using Burrows-Wheeler Aligner (BWA) (Li & Durbin, 2009).The quality-filtered reads remained for de novo assembly using IDBA-UD (Peng et al., 2012) with the following parameter settings: -mink 40, -maxk 120, -min_contig 200.Coding genes were predicted using MetaGeneMark on scaffolds of >300 bp (Zhu, Lomsadze & Borodovsky, 2010).A non-redundant gene catalog was constructed using CD-HIT, with a sequence identity cutoff of 0.95 and a minimum coverage cutoff of 0.9 for the shorter sequences.Gene abundance in each sample was calculated using SOAP alignment (identity = 0.95).Species rarefaction curves were generated to assess the sequencing depth using QIIME 1.9.1.Alpha diversity indices (Shannon, ACE, and Chao1) were compared using QIIME 1.9.1.Microbial profiles of the taxonomy of the non-redundant genes were classified by aligning the obtained sequences against the NCBI NT database using BLASTN (e-value <0.001).DIAMOND (v0.9) was used to obtain functional profiles of the non-redundant genes by aligning them with the KEGG database and the Carbohydrate-Active EnZymes database (CAZy) (Buchfink, Xie & Huson, 2015).All acquired raw sequences in this study have been deposited in the Sequence Read Archive of the Chinese National Microbiology Data Center with the accession number NMDC40026007-NMDC40026042.

Statistical analyses
All statistical analyses were performed using R version 4.0.2(R Core Team, 2020).Differences in the composition and diversity of microbiomes among the rhizocompartments of bamboo were statistically compared by ANOVA using Duncan's multiple range test.The differences between the root zone soils of bamboo and spruce were evaluated using the function ''t.test'' in the stats package.Taxonomic and functional structure was visualized using principal coordinate analysis (PCoA) on Bray-Curtis distances.PERMANOVA was conducted to test the differences in beta diversity among rhizocompartments using the ''adonis'' function in ''vegan.''The correlation between the microbial community (at the phylum level) and soil physicochemical properties was tested using the Pearson coefficient.The effects of soil physicochemical properties on the soil microbial community were tested via redundancy analysis (RDA) using the vegan package in R. The correlations among the soil properties, genus-level taxonomic composition of microbiomes, and gene family distances (CAZy and KEGG) in the metagenome were determined using the Mantel test in the R vegan package (Dixon, 2003).

Comparison of soil physicochemical properties
The results in Fig. 1 show that TN, TC, NH 4 + -N, MBC content, and the activities of invertase and urease were significantly higher in the bamboo root zone than in the spruce root zone soils.No obvious differences in the MBN and NO3-N contents, pH and EC values, or protease and cellulase activities were observed between these two soil samples.

Microbial community composition
A total of 3,068,211,656 bp of DNA sequences were obtained after quality control, and a summary of the sequencing data is presented in Table S1.Rarefaction curves (Fig. S1) indicated that sufficient sequencing depth was achieved in this study.There were significant differences in the taxon numbers at the domain level of the microbiomes among the root endosphere and soil samples (Table S2).Bacterial taxa accounted for 74.45% of the microbes in root endosphere of bamboo, which was significantly lower than that in the rhizosphere and root zone, whereas eukarya showed the opposite pattern.Archaea accounted for 0.07% of the microbes in root endosphere, being significant lower in the root endosphere than that in the rhizosphere and root zone of bamboo.No significant differences were detected in archaea, bacteria, or eukarya between the root zone soils of bamboo and spruce.The root endosphere and soil microbiomes were dominated by Proteobacteria (59.81%-71.84%),Actinobacteria (6.45%-13.01%),Acidobacteria (1.90%-5.03%),Ascomycota (0.16%-3.12%), and Planctomycetes (1.20%-1.74%),although the relative abundances varied among the samples (Fig. 2 and Table S3).The abundance of dominant bacterial phyla (such as Proteobacteria and Actinobacteria) in the root endosphere was significantly lower than that in the rhizosphere and root zones.Abundances of the nine most dominant fungal phyla, Ascomycota (comprising 56.65%-59.71% of fungal sequences) and Basidiomycota (21.98%-32.46% of fungal sequences) were significantly greater in the endosphere of bamboo than in the other rhizocompartments (Fig. 2).Representatives of archaea, such as Euryarchaeota, Thaumarchaeota, and Crenarchaeota, were significantly less abundant in the root endosphere than in the other rhizocompartments.At the phylum level, the abundance of bacteria in Aquificae, Thermotogae, Thermodesulfobacteria, Deferribacteres, and fungi in Mucoromycota significantly differed between the root zone soils of bamboo and spruce.As the dominant genera among the samples, Pseudomonas and Burkholderia exhibited significantly higher abundances in the root endosphere than in the other rhizocompartments.Similar results were observed in Klebsiella, Massilia, and Collimonas, whereas Bradyrhizobium, Rhizobium, and Streptomyces values were reversed (Table 1).Higher abundances of Pseudomonas and Massilia were observed in the bamboo root zone than in the spruce root zone.As estimated by the Shannon, ACE, and Chao1 indices, microbial diversity in the root endosphere was significantly lower than that in the other bamboo rhizocompartments (Fig. 3).ACE richness was significantly higher in the bamboo root zone than in the spruce root zone.The composition of the microbial communities in the bamboo root endosphere differed from that in the other rhizocompartments.In contrast, the root zone soil compositions of bamboo and spruce were similar, as shown in the PCoA plot and PERMANOVA analysis (Fig. 4A and Table S4).Compared with the rhizosphere and root zone of bamboo, 243 and 242 significantly enriched genera, and 900 and 909 significantly depleted genera were detected in the bamboo root endosphere (Fig. S2).Compared with the microbiomes in the spruce root zone, 247, five, and three significantly enriched genera and 947, 55, and 59 significantly depleted genera were detected in the bamboo endosphere, rhizosphere, and root zone, respectively (Fig. S3).Furthermore, 77 genera were specific to spruce soils, and 116 genera were specific to bamboo soil, but all of them were minor groups with relative abundances ≤1 ppm (J.Xiang, 2021, pers.obs.).

Functional gene profiles of microbial communities
Using the KEGG and CAZy databases, 12,465 KO genes (Fig. 5) among the metagenomic reads and 332 CAZy (Fig. 6) functional genes related to carbon and nitrogen cycles were screened, respectively, in all the samples.Functional changes between root and soil microbiomes were investigated via PCoA and PERMANOVA analysis (Figs. 4B and 4C, Table S4), which showed a significant difference between the endosphere and the other bamboo rhizocompartments but was similar between the bamboo root zone and spruce root zone (KEGG and CAZy databases).The relative abundances of genes involved in the carbon and nitrogen cycles in the root endosphere were different from those in the other bamboo rhizocompartments (Fig. 5).For example, the abundance of genes encoding the nitrogen fixation protein NifB, beta-glucosidase, and cellulose 1, 4-beta-cellobiosidase in the root endosphere was 10-, 200-and 39-fold higher, respectively, than those in the other rhizocompartments (Table S5).In addition, the abundance of some genes involved in multiple sugar transporters in the root endosphere was significantly different (greater or less depending on the proteins) from those in the rhizosphere and root zone of bamboo.The abundance of nitrate reductase changed significantly between the root zone soils of bamboo and spruce.By annotating metagenomic sequences to the CAZy database, the most abundant enzyme class in the soil samples was the glycosyl transferases (GTs) (relative abundance, 37.15%-48.55%),followed by glycoside hydrolases (GHs) (relative abundance, 31.31-33.80%;Table S6).The relative abundances of GH, GT, and polysaccharide lyase (PL) classes were greater in the root endosphere than in the other rhizocompartments of bamboo.In contrast, those of auxiliary activities (AA), carbohydrate-binding modules (CBM), and carbohydrate esterases (CE) classes were reversed (P < 0.05).The 50 dominant CAZy families in the root endosphere were significantly different from those in the rhizosphere and root zone of bamboo (Fig. 6).In contrast, most of these families did not vary between the root zone soils of bamboo and spruce (Table S7).

Relationship between soil properties and microbiomes
Mantel tests between soil properties and microbial community distance matrices demonstrated that the soil properties were not correlated with shifts in the microbial community but were significantly correlated with the functional composition of microbiomes in the bamboo and spruce root zone soils (Table 2).The Spearman correlation coefficient revealed significant correlations between the soil properties and significantly varied microbial taxa (Fig. 7).In this study, 651 sequences belonging to 42 genera (including Pseudacidobacterium and Flavihumibacter as the dominant genera) were correlated with invertase (K01193), and 5,050 sequences belonging to more than 100 genera (including Pseudomonas, Pseudonocardia, Streptomyces, Rhodoplanes, and Microvirga) were correlated with urease (K01427, K01428, K01429, K01430, and K14048).Cellulases and proteases also correspond to many KOs belonging to over 100 genera.Invertase was significantly related to 31 microbial taxa (at the phylum level), and the effect of NH 4 + -N on microbial taxa was greater than that of NO 3 − -N.Soil physicochemical properties and enzyme activities were negatively associated with Crenarchaeota (archaea) and Elusimicrobia (bacteria).In the CAZy database, pH, TN, TC, NH 4 + -N, MBC, urease, and invertase significantly influenced CAZyme families (Table S8).The RDA results showed that soil pH and invertase were the main contributors to microbial taxa and KEGG compositions, and invertase and TN were the main contributors to CAZyme families (Fig. S4).

Filtration of rhizocompartments for microbiomes
The results obtained in the present study (Figs.2-4) and in previous studies (Woźniak et al., 2019;Xiao et al., 2017;Zhou et al., 2020) demonstrated that root filtration significantly reduced biodiversity, represented by Shannon, ACE, and Chao1 indices (Fig. 3), as well as shifted the composition (Figs. 2 and 4) and functional patterns (Figs. 5 and 6) of microbiomes in the root endosphere.Furthermore, hierarchical selection in the other  rhizocompartments (bamboo rhizosphere, bamboo root zone, and spruce root zone) was not reflected in the biodiversity indices, community composition, and functions of the microbiomes (Figs.3-5).Previously, the bacteriome structure in the endosphere, which is different from that in the rhizosphere, was also reported for microbiomes associated with Moso bamboo (Yuan et al., 2021).It was estimated that the difference was related to the rhizosphere physicochemical properties, the effects of roots on nutrient supply, and the environmental variation from the soil to the endosphere (Gottel et al., 2011).The lower abundance of the predominant bacterial phyla Proteobacteria and Acidobacteria in the root endosphere of arrow bamboo than in the rhizosphere might also result from root filtration.Although the abundance of dominant bacterial phyla decreased in the endosphere, the abundance of the main fungal phyla increased, suggesting the benefit of fungal endophytic colonization (Yan et al., 2019) by arrow bamboo roots from the surrounding soil fungal communities (Bulgarelli et al., 2012).Another possible reason for the increased abundance of fungi in the root endosphere could be related to mycorrhizal symbiosis, as F. nitida forms mycorrhiza (Luo et al., 2022), possibly with enrichment of Ascomycota, Basidiomycota, and Mucoromycota in the root endosphere (Van der Heijden et al., 2015;Bonfante & Venice, 2020).
In the present study, Pseudomonas was the predominant genus in all samples, but its abundance in the arrow bamboo root endosphere was significantly greater than that in the rhizosphere and root zone soils (Table 1).This finding was consistent with our isolation results (Zhang et al., 2023), indicating that plants of arrow bamboo prefer these bacteria as their endophytes.The endophytic Pseudomonas isolated from roots of bamboo (species name not reported) grown in Aceh Province of Indonesia presented good potential as a biological control agent against nematodes of Meloidogyne spp.(Maulidia & Sriwati, 2020).Endophytic Pseudomonas members isolated from bamboo (Bambusa blumeana) in Guangdong Province of China are renowned for their plant growth-promoting and nitrogen-fixing properties (Wei et al., 2010).In addition, many endophytes of arrow bamboo roots contain yolk-degrading bacteria and produce IAA (Zhang et al., 2023).Therefore, endophytic Pseudomonas may directly or indirectly stimulate bamboo growth and pathogen resistance (Guzmán-Guzmán & Santoyo, 2022).
The dominant bacterial genera in arrow bamboo, Pseudomonas, Burkholderia, and Klebsiella, are frequent endophytes of bamboo species (Liu et al., 2017;Singh et al., 2021).Bamboo endophytes may also colonize the guts of arrow bamboo eaters, such as giant pandas and insects (Yao et al., 2021).Being colonized in guts, these endophytes could help animals digest plant fibers or antimetabolites and synthesize metabolites with bioactivities (Martínez-Romero et al., 2021).It has been reported that endophytic bacteria in Pseudomonas, Burkholderias, and Enterobacter can digest some antimetabolites in their host plants (Shanmuganandam et al., 2019).In the gut metagenomes of the panda, an arrow bamboo eater, most genes related to the degradation of plant secondary metabolites are associated with Pseudomonas (Zhu et al., 2018).Identifying genes related to microbe-plant interactions and degradation of cellulose in Klebsiella strains isolated from panda faeces also indicated their endophytic origin and possible contribution to the digestion of cellulose in the panda gut (Lu et al., 2015).These findings enlarged the impact of plant endophytes on herbivore digestion and implied that endophytes might be a determinant of the feeding habits of some animals.This is the first report on the primary understanding of microbial function in the root and root-associated soils of arrow bamboo by analyzing the abundance of gene families and the functional classification schemes of the KEGG and CAZy databases.Endophytes commonly contain enzymes that can degrade plant cell wall composition, such as endoglucanase, cellulase (Sharma et al., 2020), and chitin degradation enzymes (Proença et al., 2018).Overall, the functional gene contents of soil microbial communities were significantly different between the endosphere and the rhizosphere (Figs.5-6), indicating that the filtration effects for microbes by the rhizocompartments also change the metabolism, biosynthesis, and degradation processes in the carbon and nitrogen cycles of arrow bamboo.

Soil and microbial properties of arrow bamboo and spruce root zone
Arrow bamboo is the pioneer and companion species of secondary forest communities in the subalpine spruce plantations of southwest China.Different plant species in aboveground systems can affect soil properties and root-associated microbes, mainly through root activity and litter input (Diao et al., 2020;Prescott & Vesterdal, 2013).The fact that higher values of soil MBC, invertase activity, urease activity, and contents of soil TC, TN, and NH 4 + -N were observed in the arrow bamboo root zone than in the spruce root zone (Fig. 1) demonstrated the effects of plant species on soil physicochemical traits.These results were similar to those of a previous study, in which greater activities of invertase and urease and greater contents of TN and microbial carbon were recorded in soils supplied with manure and mineral fertilizers (Mikanová et al., 2015).Nitrogen, as the most important nutrient, commonly constrains plant growth due to its low supply in forests, and variations in nitrogen content in soil might affect the ecological succession of forests (Högberg et al., 2017).It has been reported that the interaction of nitrogen-fixing microbes with roots could help the plant achieve an efficient NH 4 + -N supply in the rhizosphere of spruce-fir forests (Zhu, Liu & Yin, 2021).In the present study, the higher NH 4 + -N and TN concentrations in the bamboo root zone soil than in the spruce root zone were consistent with previous results reviewed by Fuke et al. (2021).This finding indicates that bamboo roots might stimulate NH 4 + -N production and immobilization and benefit the subsequent colonization of spruce trees.Urease catalyzes the hydrolysis of urea to ammonia and is an important soil enzyme that mediates the conversion of organic nitrogen to inorganic nitrogen.Invertase is a ubiquitous enzyme in soils that catalyzes the hydrolysis of sucrose to glucose (Kandeler et al., 1999).Urease and invertase are essential soil enzymes for transforming organic carbon and nitrogen into simple carbon and nitrogen sources (Li et al., 2018).In this study, the higher activity of urease and invertase in arrow bamboo might lead to greater efficacy of soil carbon and nitrogen transformation, which is also consistent with the higher microbial biomass in the root zone of arrow bamboo (Diao et al., 2020).Arrow bamboo regeneration in spruce forests may improve soil properties through the rhizosphere effect, limiting or decreasing soil degradation in spruce plantations.It is well known that the soil microbial community structure can be altered by plants via direct and indirect pathways through root exudates, litter, and nutrient requirements (Zhou et al., 2017).In the present study, soil microbial community compositions and functional profiles in the bamboo and spruce root zones were very similar, with differences only in a few phyla, genera, KO, and CAZy gene families.Therefore, the highly similar microbial community and function in the root zones of arrow bamboo and spruce are the result of selection by plant and environmental factors during forest succession.

Linkage between soil properties and microbial community
Soil properties are determinants of microbial community establishment, which could explain the presence or absence of determined taxa in the rhizosphere, even at very small distances (Lladó, López-Mondéjar & Baldrian, 2018;Zhong et al., 2018).Soil properties were also tightly associated with microbial functional genes (Table 2, Fig. S4) by studies of secondary successional chronosequences on the Loess Plateau of China (Zhong et al., 2018).Our study showed that NH 4 + -N concentration was a more important environmental factor than NO 3 − -N concentration, significantly influencing microbial functions and microbial taxa (Fig. 7 and Fig. S4 and Table S8).The greater NH 4 + -N in the root zone of arrow bamboo evidenced the effects of bamboo roots on NH 4 + -N production and immobilization, which in turn affected soil microbes (Zhu, Liu & Yin, 2021).As shown by the RDA results (Fig. S4), soil pH, TN, and invertase were the main contributors to the microbial taxa and functions.Soil pH can affect soil nutrient solubility and availability to plants and the soil microbial community (Naz et al., 2022;Zheng et al., 2023).As an essential nutrient, soil nitrogen is generally the limiting factor for plants in terrestrial ecosystems, resulting in strong competition between microorganisms and plants (Geisseler et al., 2010).Nitrogen enrichment shifts the functional genes of microbes related to nitrogen and carbon acquisition (Treseder et al., 2018).The trade-off between soil carbon and nitrogen nutrients may have a greater impact on genes in CAZy families (Cardenas et al., 2018).

CONCLUSIONS
The filtration of arrow bamboo roots significantly decreased the diversity of the endophytic microbiome.It made the endomicrobiome structure different from that in its rhizosphere and root zone, with a reduced abundance of bacterial phyla Proteobacteria, Actinobacteria, and Acidobacteria and an increased abundance of fungal phyla Ascomycota, Basidiomycota, and Mucoromycota.Pseudomonas was the most abundant and enriched bacterium in the root endosphere.Genes involved in nitrogen fixation, cellulose degradation, and multiple sugar transporters in the root endosphere differ from those in the rhizosphere and root zones of bamboo.The soil TC, TN, NH 4 + -N, MBC, and invertase and urease activities in the bamboo root zone were greater than those in the spruce root zone, while the soil microbial community and functional profiles were similar.Soil properties were not related to soil microbial communities but were significantly correlated with microbial functions.The insights into the microbial community and functional structuring gained in this study provide a basis for understanding that the host plant and root zone soil properties select

Figure 2
Figure 2 Relative abundances of the bacterial (A), fungal (B) and archaeal (C) phylotypes in the metagenomes of microbiome in each sample, and their major variations in root endosphere, rhizosphere and root zone soils.The other Eucaryota and viruses occupied the remaining parts of the relative abundances of the sequencing data (TableS2).Relative abundances are calculated as the proportion of individual taxa to the total number of corresponding phylotypes reads.BE, arrow bamboo root endosphere; BR, arrow bamboo rhizosphere; BZ, arrow bamboo root zone; and SZ, spruce root zone.Full-size DOI: 10.7717/peerj.16488/fig-2

Figure 3 Figure 4
Figure 3 The microbial community diversity in the roots and soils.Different lowercase letters represent significant differences among bamboo samples using Duncan's multiple range test.T test was used for comparing root zone soil of bamboo and spruce * p < .05.BE, bamboo root endosphere; BR, bamboo rhizosphere; BZ, bamboo root zone; SZ, spruce root zone.Full-size DOI: 10.7717/peerj.16488/fig-3

Figure 5
Figure 5 Heatmaps indicating functional genes abundances of the C and N cycles in the KEGG database that differ between the roots and soils.BE, arrow bamboo root endosphere; BR, arrow bamboo rhizosphere, BZ, arrow bamboo root zone; and SZ, spruce root zone.Full-size DOI: 10.7717/peerj.16488/fig-5

16488 7/22 Table 1 Comparative analysis of relative abundance of the domain (top 15) genera in the roots and soils.
Zhang et al. (2023), PeerJ, DOI 10.7717/peerj.Notes.Different lowercase letters represent significant differences among bamboo samples using Duncan's multiple range test.T test was used for comparing root zone soil of bamboo and spruce.*p < .05.BE, bamboo root endosphere; BR, bamboo rhizosphere; BZ, bamboo root zone; SZ, spruce root zone.