Search tips
Search criteria 


Logo of eceLink to Publisher's site
Ecol Evol. 2017 September; 7(18): 7503–7514.
Published online 2017 August 11. doi:  10.1002/ece3.3307
PMCID: PMC5606866

Rhododendron aureum Georgi formed a special soil microbial community and competed with above‐ground plants on the tundra of the Changbai Mountain, China


Rhododendron aureum Georgi is a perennial evergreen dwarf shrub that grows at all elevations within the alpine tundra of northern China. Previous research has investigated the plant communities of R. aureum; however, little information is available regarding interspecific competition and underground soil microbial community composition. The objective of our study was to determine whether the presence of R. aureum creates a unique soil microbiome and to investigate the relationship between R. aureum and other plant species. Our study site ranged from 1,800 to 2,600 m above sea level on the northern slope of the Changbai Mountain. The results show that the soil from sites with an R. aureum community had a higher abundance of nitrogen‐fixing bacteria and a higher resistance to pathogens than soils from sites without R. aureum. We emphasize that R. aureum promotes a unique soil microbial community structure that is distinct from those associated with other plants. Elevation and microbial biomass were the main influencing factors for plant community structure. Analysis of interspecific relationships reveals that R. aureum is negatively associated with most other dominant shrubs and herbs, suggesting interspecific competition. It is necessary to focus on other dominant species if protection and restoration of the R. aureum competition is to occur. In the future, more is needed to prove whether R. aureum decreases species diversity in the tundra ecosystems of Changbai Mountain.

Keywords: function microbial, interspecific competition, microbial biomass, microbial community, Rhododendron aureum, soil microbial


Soil microbes play a key role in ecosystems and mediate many ecological processes that are critical to decomposition, nutrient cycling, and plant community dynamics(Fierer et al., 2012; Van der Putten, Klironomos, & Wardle, 2007). In addition, they can act as a source and sink for plant nutrients (Nannipieri et al., 2012). The interaction between plants and soil microbes is at the forefront of international ecological research (Mendes et al., 2011; Park et al., 2006). Studies have shown that different plant species host‐specific microbial communities, suggesting that plant species play an important role in shaping the soil microbiome (Berendsen, Pieterse, & Bakker, 2012; Oh et al., 2012). Microbial biomass and activity are typically thought to be constrained by soil nutrients (Liu et al., 2015). Studies reveal that plant production and diversity increases the soil microbial biomass by altering the environmental factors, for instance by increasing the concentration of nitrogen and organic matter (Eisenhauer et al., 2013; Lange et al., 2015). Furthermore, the microbial community composition of soil has a large impact on plant–plant interactions (Bever, Platt, & Morton, 2012) and, consequently, on plant diversity and composition (Heijden et al., 2006; Schlatter, Bakker, Bradeen, & Kinkel, 2015).

So far, there were many researches have reported about the shifting in microbial community structure and functional genes along a large‐scale alpine climosequence (Collins, Carey, Aronson, Kopp, & Diez, 2016; Cui et al., 2016; Zhang, Liang, He, & Zhang, 2013). At relatively small scales, individual plants can also alter the composition of the soil community (Bever et al., 2010; Putten et al., 2013); however, there is currently a relatively limited understanding of how different plants affect the microbial community structure and function in alpine tundra that was not formed with a vertical distribution of vegetation.

In the alpine environment, interspecies relationships within the community are important for plants. Many studies have confirmed widespread positive interactions (facilitation) within alpine plant communities (Baumeister & Callaway, 2006; Schob et al., 2014). Previous research has also demonstrated the existence of competition or “parasitic” interactions between benefactor and beneficiary species in alpine ecosystems (Kikvidze et al., 2005; Maestre, Callaway, Valladares, & Lortie, 2009; Schob et al., 2014). Moreover, soil microbes have a substantial impact on plant productivity, community composition, and diversity (van der Heijden, Bardgett, & van Straalen, 2008). Because plants cannot fix atmospheric nitrogen, which is a main limiting factor for their growth, nitrogen‐fixing bacteria are important regulators of plant productivity (Iii, 2003). Soil microbes can decompose soluble and insoluble organic matter and convert it into organic and plant‐effective forms (Raab, Lipson, & Monson, 1999). Accelerated carbon and nitrogen cycling, as well as decomposition in the soil, stimulate microbial activity, thereby enhancing shrub development (Bret‐Harte et al., 2001; Hobbie, Nadelhoffer, & Högberg, 2002; Naito & Cairns, 2011; Sturm et al., 2005; Syndonia, Shaver, & Stuart, 2002).

Alpine tundra of the Changbai Mountain is one of two rare alpine tundra ecosystems in China (Zhao, Qi, Lyu, Yu, & Chen, 2016). Rhododendron aureum Georgi is a perennial evergreen shrub found in Northeast Asia that is able to endure the cold alpine climate (Kudo, 1993; Zhao et al., 2016). It is the only woody plant distributed at all elevations throughout the tundra and it plays an important role in maintaining the ecological balance, while preventing and controlling soil erosion. These unique plant communities can produce specific physical and chemical properties within the soil through their litter and root exudates, thus influencing the soil microorganisms (Cheng et al., 2015). Past studies have focused on the plant communities (Jin et al., 2015; Liu et al., 2012; Zhao et al., 2016), with little research into interspecific competition and underground soil community composition.

Rhododendron aureum is highly adaptive and can thrive at various altitudes. It is crucial to understand the complexity of this adaptive processes and the main forces driving population structure (Joost et al., 2007; Mosca, Gonzálezmartínez, & Neale, 2014). In this study, we (1) hypothesize that on a small‐scale, R. aureum influences the soil microbial community, distinct from other species; (2) explore the relationship between R. aureum and other species in alpine tundra; (3) discuss the interaction between plant community structure and soil properties. To test our hypothesis, we set an altitudinal gradient ranging from 1,800 m to 2,600 m above sea level using a systematic sampling design. We established pairwise treatments at the same elevation, consisting of sites with R. aureum and sites with other species but without R. aureum. We then analyzed the soil from both treatment types with regards to soil properties, N‐cycling, microbial function, and microbial community composition. Finally, we evaluated the diversity variation within the native above‐ground communities and predicted the biotic interactions.


2.1. Study area

Sampling sites were established on the northern slope of the Changbai Mountain which is acknowledged as the highest mountain in northeastern China and eastern Eurasia. The Changbai Mountain has clearly defined vertical vegetation zones: deciduous broad‐leaved forest, mixed deciduous broad‐leaved/conifer forest, dark coniferous forest, Betula ermine forest, and tundra (Zheng, Wallin, & Hao, 1997). The study area was located on the northern slope at elevations from 1,800 m to 2,600 m above sea level (Figure 1). The climate is characterized by low temperatures, heavy precipitation, and a short‐growing season. The average annual temperature is 3–7°C, with an annual precipitation is over 600 mm (Yang & Wu, 1998). In accordance with the alpine climate, common plants are of low stature or prostrate shrubs, with R. aureum being the most dominant evergreen shrub in the study area. Other plants grow here include Dryas octopetala, Vacciniumu liginosum L., Ligularia jamesii (Hemsl.) Kom, Trollius chinensis, Rhodiola cretinii, and Sanguisorba tenuifolia var. (Yang & Wu, 1998; Zong et al., 2014).

Figure 1

(a) Distribution of eight sites sampled along an elevation gradient on the tundra of Changbai Mountain. We extracted the elevational information from DIVAGIS Free Spatial Data ( and created the map by the program ...

2.2. Experimental design and sampling

Fieldwork was conducted from 2014 to 2015 along an altitudinal gradient ranging from 1,800 to 2,600 m, using a systematic sampling design. A total of 16 plots from eight sample sites were established in different altitudinal sites covering all the dominant species associated with R. aureum in the alpine ecosystem of Changbai Mountain. In each site, we established plots consisting of at least 90% R. aureum and <10% non‐R. aureum species (with R. aureum; N1‐N8) and plots with 100% non‐R. aureum species (without R. aureum; X1–X8). Each plot measured 20 m × 20 m and contained 12 shrub subplots, measuring 5 m × 5 m, and 12 herb subplots, measuring 1 m × 1 m. The foliage coverage of vascular plant species in each plot was recorded as percentages, along with their height, by visual estimation of each species.

Soil samples were collected in mid‐July 2015 at the 16 established plots. At each plot, 12 soil cores at depths of 0–10 cm were collected, using a core sampler with a 10 cm inner diameter at a distance of 8 m around the plot center. The 12 cores were pooled to obtain one composite sample for each plot and immediately stored on ice in insulated containers. After returning to the laboratory, each sample was processed by passing it through 2‐mm mesh and dividing it into two groups: one group was stored at 4°C for chemical analyses and the other one at −80°C for microbial community analysis.

2.3. Analysis of soil characteristics

Soil moisture was determined after drying at 105°C; pH was determined by mixing soil and deionized water at a 1:5 mass to volume ratio. Soil total nitrogen (TN) and alkaline hydrolysis nitrogen (AN) were measured by the semimicro‐Kjeldahl (KDY‐9820) digestion and alkali diffusion method, respectively (Bao, 2000). Available N (ammonium and nitrate) was extracted with 2 mol/L KCl and then measured by indophenol‐blue and phenoldisulfonic acid colorimetry, respectively (Bao, 2000). Soil total phosphorus (TP) and total potassium (TK) were determined by the molybdate method (Murphy & Riley, 1962) and extracted by incubation with sodium hydroxide, respectively. Available phosphorus (AP) was determined by the Olsen method (Olsen et al., 1954), and available potassium (AK) was determined using the ammonium acetate extraction procedure (Bao, 2000). Microbial biomass, C and N (MBC and MBN), was extracted with 0.5 mol/L K2SO4 and analyzed using the chloroform fumigation extraction method (Brookes, Landman, Pruden, & Jenkinson, 1985). Phosphatase activity was determined with the modified Schinner and von Mersi method (Schinner & Mersi, 1990). Catalase activity and invertase activity were determined using the 0.1 N KMnO4 titration method (Johnson & Temple, 1964) and the 3,5‐dinitrosalicylic acid technique (Vaughan & Ord, 1980), respectively. Urease activity was determined according to Klose and Tabatabai (Klose & Tabatabai, 2000).

2.4. Soil DNA extraction

DNA was extracted from 0.3 g of freeze‐dried soil using the Power Soil DNA isolation kit (MoBio). Extracted DNA was purified using the GV‐High‐Efficiency Agarose Gel DNA Purification Kit. DNA concentrations were checked with the Qubit quantification platform. DNA was diluted to 10 ng/μl using TE buffer and then stored at −80°C prior to molecular analysis.

The V3‐V4 hypervariable region of the bacterial 16S rRNA gene was amplified using the primers U341F (5′‐ACT CCT ACG GGA GGC AGC AG‐3′) and U806R (5′‐GGA CTA CHV GGG TWT CTA AT‐3′). An equal amount of PCR product for each sample, as measured by Qubit quantification platform, was pooled in a single tube to be sequenced using the Illumina MiSeq platform (Realbio Genomics Institute, Shanghai, China). The raw data were then subjected to quality control using UPARSE (Edgar, 2013). The sequences were quality‐filtered and clustered into operational taxonomic units (OTUs) at a 97% similarity level using USEARCH (Edgar, 2010). A representative sequence for each OTU was assigned to a taxonomic level using the RDP (Ribosomal Database Project) database by the RDP classifier (Maidak et al., 1996).

The relative abundance of genes encoding key enzymes for biological N‐cycling (i.e., nifH, archaeal amoA, bacterial amoA, and nosZ) was assessed by qPCR (ABI 7500). The gene fragments were amplified using the primers PolF/PolR for nifH, Arch‐amoAF/Arch‐amoAR for archaeal amoA, amoA‐1F/amoA‐1R for bacterial amoA, and nosZ‐F/nosZ‐1622R for nosZ(Levy‐Booth, Prescott, & Grayston, 2014; Mao, Yannarell, Davis, & Mackie, 2013). Primers and qPCR conditions are described in Table S1. Standards for the qPCR assays were generated with PCR amplicons from the pooled DNA sample (Chen, Yu, Jr, Wittum, & Morrison, 2007). The copy number of the genes in each standard was calculated by dividing the DNA concentration (ng/μl) by the average molecular weight of the amplified gene fragment (Mao et al., 2013). qPCR amplification efficiencies ranged from ~90% to 110%, with R 2 values between 0.992 and 0.999.

2.5. Statistical analysis

The coverage (%) of plant communities was transformed to an ordinal scale and assigned to one of four coverage classes according to the 1–4 Lagerberg–Raunkiaer scale (Ludwig & Reynolds, 1988) and used for the interspecies association analysis (Guozhen Du & Wang, 1989). The importance value of a species, which is considered the most realistic and important element in assessing vegetation inventory, is useful in comparing the ecological significance of species in a given area (Curtis & McIntosh, 1951; Lamprecht, 1989). This value was calculated from the sum of the relative density, relative frequency, and relative dominance, as recommended by Kent and Coker (1992). The species diversity was determined using the Shannon–Wiener diversity index (H’) and Pielou index (Pielou, 1969).

General linear model analysis (mean ± SD, n = 5) of variance was performed to test differences between sites with and without R. aureum. Tukey's LSD post hoc test was used to identify significant differences. Statistical analyses were performed using SPSS (v. 19.0 for Windows, SPSS Inc., Chicago).

The number of OTUs in a community was used to estimate its richness and structure. Alpha diversity was measured by the Chao1 and Shannon indices. Beta diversity of the microbial communities was investigated using multivariate analysis and nonmetric multidimensional scaling (NMDS). A Bray–Curtis distance matrix was calculated using the soil microbial OTU data using the vegan package in R statistical software (v.3.1.1; Oksanen, Kindt, Legendre, & O'Hara, 2010). Treatment differences among these metrics were assessed using analysis of similarities (ANOSIM) in R. Additionally, redundancy analysis (RDA) was performed, which is an extension of principal components analysis, where the main components are constrained to linear combinations of environmental variables. To determine the environmental factors that significantly correlate with community composition (abundance of OTUs or relative abundance of species), the envfit function (999 permutations) in the vegan package of R was used.


3.1. Soil nutrient dynamics

The thirteen soil parameters varied across the two treatment sites (with and without R. aureum; Table 1). Sites with R. aureum were generally more nutrient rich in terms of nutrient pools (TN, TOC, TP, and TK), moisture, and soil pH; however, these values were not significantly different when compared to the sites without R. aureum (> .05). Only microbial biomass carbon and nitrogen (MBC and MBC) were significantly higher (< .05) than sites without R. aureum. To investigate the cause of these changes, we performed correlation analysis between soil properties and elevation (Table S2). Out of nine parameters were significantly correlated with elevation (Table S2). Among them, TK and pH positively correlated with elevation, while other parameters negatively correlated. Furthermore, four soil enzymatic activities were measured (Figure 2). When comparing the two treatment types, sites with R. aureum contained a higher amount of enzymatic activity. Specifically, sucrase and urease were significantly higher in sites with R. aureum (< .05, Figure 2).

Figure 2

Soil enzyme activities difference between the two treatments sites (with and without R. aureum). (a): Urease, p = .03; (b): Sucrase, p = .015; (c): Acid phosphatase, p = .637; (d) Catalase, p = .165. ...

Table 1

Soil variables (mean ± SD) were from two treatments

The abundance of N‐cycling genes was not significantly different (> .05) between treatment sites (Figure 3), with the exception of the bacterial N‐fixing gene nifH. The abundance of nifH was higher in sites with R. aureum than without (< .05; Figure 3). Correlation analysis showed that the abundance of N‐cycling genes was significantly correlated (< .05) with some soil parameters (Table 2). Abundance of the nifH gene was positively correlated with total nutrients (TOC and TN), available nutrients (AN and AP), microbial biomass (MBC and MBN), and soil moisture (Table 2). The abundance of AOA (archaeal amoA) was negatively correlated with soil nitrate nitrogen (NO3 ); however, the abundance of N‐cycling genes was not significantly correlated with ammonium nitrogen (NH4 +) and AK. The abundance of AOB (bacterial amoA) was negatively correlated with AK, but positively correlated with soil pH and elevation; however, abundance of other N‐cycling genes was not significantly correlated with elevation. The abundance of nosZ was positively correlated with TP and AP (Table 2).

Figure 3

Real‐time PCR quantification of the N‐cycling function genes (nifH, archaeal amoA, bacterial amoA, and nosZ) in the bulk soil under two treatments. The copy number of genes in each gram of dry soil was estimated based on the results of ...

Table 2

Relationships between the abundance of N‐cycling function genes and soil properties

3.2. Soil microbial communities

Sequencing of the 16S rRNA gene resulted in 431,560 sequences from 16 samples following the removal of low‐quality sequences, chimeras, and denoising. The number of sequences per sample ranged from 19,009 to 31,653, with an average sequence count of 26,973. These sequences represented 2,022 OTUs at 97% similarity. The bacterial community composition was relatively consistent between soil samples from sites with and without R. aureum. At the phylum level, the majority of sequences corresponded to Acidobacteria (30.6%), Proteobacteria (29.6%), and Actinobacteria (18.7%; Figure 4a). Species annotation revealed that R. aureum‐associated soil contained 245 bacterial genera and soil from sites without R. aureum contained 236 bacterial genera, with 190 shared genera between them (Figure 4b).

Figure 4

(a) Relative abundances of the dominant bacterial phyla in soils separated according to elevation and species categories. Relative abundances are based on the proportional frequencies of those DNA sequences that could be classified at the phylum level. ...

Alpha diversity analysis shows that the microbial diversity of soils with R. aureum was lower than that of soil samples without R. aureum, but not significantly (Table S3); however, the relative abundance of genera associated with nitrogen fixation and pathogenesis were significantly higher within R. aureum‐associated soil (Figure 5).

Figure 5

Relative abundances of the dominant function bacterial at gene level. (a) Pathogenic properties bacterial, < .001; (b) nitrogen fixation bacteria, p = .04

The heatmap shows pairwise comparisons of bacterial community structures between all samples (Figure 6a). Results indicate that the sites with R. aureum cluster into one clade, while the sites without R. aureum cluster into another. NMDS ordination also reveals that microbial communities were significantly different between soils with and without R. aureum (= .01; Figure 6b). In order to better understand the influence of soil microbial communities, we investigated the interdependence between 16S rDNA transcript abundance of individual phylotypes and single environmental factors using RDA (Figure 7a). Of all the environmental variables tested, TN (r 2 = 0.43, p = .02) and pH (r 2 = 0.54, p = .03) demonstrated the highest correlation with community composition (Table S4). Other factors, such as elevation (r 2 = 0.32, p = .09) and AP (r 2 = 0.32, p = .09) also showed a high correlation with bacterial community composition (Table S4). The abundance of Gemmatimonadetes (r = .686, p = .003), Bacteroidetes (r = .706, p = .002), and Nitrospirae (r = .694, p = .003) was positively correlated with soil pH. TN was negatively correlated with Gemmatimonadetes (r = −.541, p = .03) and positively correlated with Fusobacteria (r = .776, p = .001; Fig. S1).

Figure 6

(a) Pairwise comparison of β‐diversity with all samples and annotations. Clustering and heatmap were computed using R package. (b) Nonmetric multidimensional scaling (NMDS) plot of community composition based on pyrosequencing. N indicates ...

Figure 7

(a) Redundancy analysis (RDA) triplots of 16S rDNA fingerprint patterns, showing the contribution of 14 environmental parameters to variability. (b) Redundancy analysis (RDA) triplots of plant species richness showing the contribution of 14 environmental ...

3.3. Interspecific relationships

Within the 16 sampling plots, a total of 81 plant species were encountered, belonging to 31 families and 58 genera. Based on the overall importance values, a relative measure accounting for the sum and frequency of coverage, the most important species were R. aureum, Vaccinium uliginosum, Saussurea tomentosa, and Streptopus koreanus Ohwi. The vast majority (77%) of the species had importance values below 1%, indicating that they were locally rare. Variance analysis shows that the overall interspecific correlation of the nine dominant shrubs was significantly positive, (VR = 1.55, W = 148.73, χ2 (0.95)96 = 74.4), and that the distribution of the species was not independent; however, when using individuals to assess the abundance index, there is a significant negative association between R. aureum and other accompanying species at the 99% confidence level (Table 3).

Table 3

Interspecies association between Rhododendron aureum and its coexisting species

RDA was used to investigate the impact of plant community structure (Figure 7b). The results show that elevation, MBC, and MBN significantly correlate with community composition (Table S5). Many of the dominant species positively correlate with each other (Figure 7b), thereby contributing to the overall interspecific association observed; however, R. aureum negatively correlates with half of the dominant species, including both shrub and herbaceous species (Figure 7b). This suggests that R. aureum either competes for resources or interferes with other species.


A species can survive and become widely distributed in an alpine environment not only based on its ability to adapt to abiotic stresses, but also its ability to cope with biotic influences. Both competition and facilitation must be considered in order to obtain a realistic view of the interactions within alpine plant communities (Kikvidze, 2002; Olofsson, 2004). The overall interspecific association among dominant shrubs on the Changbai Mountain exhibits a significant positive association. Many studies propose the stress gradient hypothesis, in which plants exhibit positive relationships with increasing stress in alpine ecosystems, especially in high altitude areas (Callaway et al., 2002; Choler, Michalet, & Callaway, 2001; Reid, Lamarque, & Lortie, 2010). Changbai Mountain is a typical alpine environment with low temperatures, heavy precipitation, short‐growing seasons, strong winds, high solar radiation, and low resource availability. Facilitation appears to be driven by architecture‐mediated protection from winds and low temperatures, as well as avoiding resource competition (Baumeister & Callaway, 2006). Nevertheless, R. aureum shows significant negative associations with most dominant shrubs and herbs, suggesting interspecific competition. A negative relationship between species can occur when resource abundance drives abiotic stress (Kikvidze et al., 2005; Maestre et al., 2009), especially in the infertile alpine soil.

In this study, we first examined whether the R. aureum plant community can form its own unique soil environment, which is distinct from other plant communities. A deeper understanding of the effects of vegetation types on the soil microbiome in alpine ecosystems can help to develop better strategies to protect the harsh environment of the alpine tundra (Klimeš, 2006). We started with the hypothesis that R. aureum supports a unique soil system allowing for its expansion and succession within the subalpine ecosystem; however, contrary to this, no significant difference was seen for most soil nutrients between sites with and without R. aureum. Many of the soil properties significantly correlated with elevation, which is consistent with previous studies (Wu et al., 2015). This suggests that elevation strongly influences soil properties, such as total nitrogen, organic matter, and some enzymatic activities, which decrease as elevation increases. Some soil properties were significantly different between the vegetation types in this study. For instance, soil microbial biomass (MBN and MBC), urease, and sucrase were higher at sites with R. aureum. A large amount of organic carbon (C) conversion, storage, and respiration by microorganisms in the soil (Spohn, 2015), and in soil microbial biomass, is an important determinant for soil carbon dynamics (Thakur et al., 2015). Thakur et al. (2015) suggested that higher plant diversity has a strong effect on soil microbial communities and could increase microbial biomass. Zak, Holmes, White, Peacock, and Tilman (2003) also revealed that, within a given region, microbial biomass increases with plant productivity. Furthermore, recent studies have highlighted that higher soil microbial biomass is driven by an input of organic matter, nitrogen, and phosphorus, as well as regulation of soil moisture (Lange et al., 2015; Liu et al., 2015; Thakur et al., 2015). Greater total nutrients were observed in soils with R. aureum than in soils without R. aureum, indicating that plant community establishment plays a crucial role in microbial biomass accumulation.

A greater abundance of nitrogen‐fixing bacteria, extracellular enzymes, and nitrogen‐cycling functional groups in soils with R. aureum could increase nutrient availability. Nitrogen is one of the most important limiting factors for organism productivity; thus, nitrogen‐fixing organisms play a key role in microbial and plant biomass production, as well as organic matter accumulation in soils (Jhp et al., 1997). Nitrogen‐fixing bacteria are important functional microorganisms in soil quality (Ahmad et al., 2016) because they promote the establishment and growth of plants by increasing the nutrient supply (Jasper, Abbott, & Robson, 1989; Liang, Pan, He, Chen, & Su, 2016). Our results show that the abundance of N‐fixing genes has a significant positive correlation with soil nutrients, especially total nitrogen and available nitrogen. This is consistent with our hypothesis that the R. aureum‐associated soil harbors microbes that provide nutrition to promote R. aureum growth.

Soil pathogens have distinct spatial and temporal patterns within the natural plant community (Van der Putten, 2003). Plant succession can also be driven by soil pathogens when the accumulation of pathogens leads to the extinction of the dominant plant and replacement with other plant species that are immune or insensitive to the pathogens (Putten, Dijk, & Peters, 1993). R. aureum is relatively pathogen tolerant, as demonstrated by the higher abundance of pathogens and lower abundance of antibiotic‐producing bacteria in its associated soil. The observed pathogens would likely arrest the growth of other species that are sensitive to them. Due to the competitive relationship between R. aureum and other species, the higher relative abundance of pathogens may threaten the survival of other species, so as to reduce species diversity. In the future, it is necessary to pay attention to the influence of R. aureum growth on the plant diversity in alpine tundra.

In‐depth analysis showed that the microbial community in soils with R. aureum had lower diversity than soils without R. aureum; however, the difference was not significant, even though there was a difference between the microbial community structures in the different vegetative soils. Moreover, NMDS ordination revealed that the soil microbial structure was distinct between sites with and without R. aureum. The heatmap showing pairwise comparisons of the bacterial community structures between all samples agrees with the NMDS results. Furthermore, based on the results of the RDA (Figure 6a), the bacterial community composition was significantly affected by soil pH and total nitrogen. Correlation analysis indicates a significant positive relationship between soil pH and the three dominant microbial phyla (Gemmatimonadetes, Bacteroidetes, and Nitrospirae). Furthermore, total nitrogen significantly correlated with two of the dominant microbial phyla (Gemmatimonadetes and Fusobacteria). Soil bacteria are more sensitive to factors such as pH, moisture, nutrient availability, and C/N ratios (Liu et al., 2014; Uroz, Tech, Sawaya, Frey‐Klett, & Leveau, 2014). Shen et al. (2013) suggested that the bacterial community structure in Changbai Mountain soils at different altitudes is mainly affected by pH. Different plant types have different effects on the soil (Wu et al., 2015; Zhang, Liu, Xue, & Xiao, 2013), and soil microbial community structure over altitudinal gradients is mainly governed by changes in the vegetation cover type (Singh et al., 2014).

Plants can change the physical, chemical, and biological characteristics of the surrounding soil environment, thereby affecting the future growth of plants in that area. In order to study the interaction between plant community structure and soil properties, RDA was used to demonstrate that microbial biomass (MBN and MBC) and elevation significantly affected plant community structure (Figure 6b). Many past studies have shown that microbial biomass contributes to plant productivity; however, our research provides evidence that it also affects the structure of the plant community. In this analysis, not all soil properties significantly affect the plant community, which may be related to plant growth characteristics. The vertical distribution of vegetation in Changbai Mountain is known to mirror the horizontal vegetation zonation from temperate to frigid zones on the Eurasian continent (Xu, He, Chen, & Liu, 2004). This study also shows that elevation drives the spatial distribution of plant communities on the Changbai Mountain.

In summary, we demonstrate that the nitrogen‐fixing and pathogenic bacteria are more abundant in soil with an R. aureum community than in soils without this community. Despite the lack of significant variation in microbial diversity, the microbial communities in the two groups could be clearly divided into two clusters. Therefore, we conclude that R. aureum facilitates a unique soil microbial community structure, which is distinct from other plant species. The plant community structure significantly correlates with elevation, as well as microbial biomass carbon and nitrogen. Elevation may drive the spatial distribution of plant communities on Changbai Mountain. In addition, R. aureum exhibits a negative association with most other dominant shrubs and herbs, suggesting interspecific competition.


None declared.


W.X.L. and L.L. conceived the study. L.L., W.X.L., Z.W., and Z.J.X. performed the experiments. W.X.L. and C.X. interpreted the results and W.X.L. wrote the manuscript.

Supporting information


We are grateful to the staff of Changbai Mountain Academy of Sciences for study area access and vital logistic support during fieldwork. We thank Dr. Ming Xing for ordering the reagents. We thank two anonymous reviewers provided suggestions that greatly improve an earlier version of this manuscript.


Wang X, Li L, Zhao W, Zhao J, Chen X. Rhododendron aureum Georgi formed a special soil microbial community and competed with above‐ground plants on the tundra of the Changbai Mountain, China. Ecol Evol. 2017;7:7503–7514.


  • Ahmad P., Mishra R. K., Mishra V., Pandey H., Pandey A. C., Sharma S., & Dikshit A. (2016). Mycorrhizal symbiosis. Chichester, UK: John Wiley & Sons Ltd.
  • Bao S. D. (2000). Analysis of soil and agricultural chemistry. Beijing: Agriculture Press.
  • Baumeister D., & Callaway R. M. (2006). Facilitation by Pinus flexilis during succession: A hierarchy of mechanisms benefits other plant species. Ecology, 87, 1816–1830. [PubMed]
  • Berendsen R. L., Pieterse C. M. J., & Bakker P. A. H. M. (2012). The rhizosphere microbiome and plant health. Trends in Plant Science, 17, 478–486. [PubMed]
  • Bever J. D., Dickie I. A., Facelli E., Facelli J. M., Klironomos J., Moora M., … Zobel M. (2010). Rooting theories of plant community ecology in microbial interactions. Trends in Ecology & Evolution, 25, 468. [PubMed]
  • Bever J. D., Platt T. G., & Morton E. R. (2012). Microbial population and community dynamics on plant roots and their feedbacks on plant communities. Annual Review of Microbiology, 66, 265–283. [PMC free article] [PubMed]
  • Bret‐Harte M. S., Shaver G. R., Zoerner J. P., Johnstone J. F., Wagner J. L., Chavez A. S., … Laundre J. A. (2001). Developmental plasticity allows Betula nana to dominate tundra subjected to an altered environment. Ecology, 82, 18–32.
  • Brookes P. C., Landman A., Pruden G., & Jenkinson D. S. (1985). Chloroform fumigation and the release of soil nitrogen: A rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biology & Biochemistry, 17, 837–842.
  • Callaway R. M., Brooker R. W., Choler P., Kikvidze Z., Lortie C. J., Michalet R., … Cook B. J. (2002). Positive interactions among alpine plants increase with stress. Nature, 417, 844–848. [PubMed]
  • Chen J., Yu Z., Jr M. F., Wittum T., & Morrison M. (2007). Development and application of real‐time PCR assays for quantification of erm genes conferring resistance to macrolides‐lincosamides‐streptogramin B in livestock manure and manure management systems. Applied & Environmental Microbiology, 73, 4407–4416. [PubMed]
  • Cheng F., Wei X., Lin H., Shang Z., Peng X., Zhao P., … Zhang S. (2015). Soil fungal communities of montane natural secondary forest types in China. Journal of Microbiology, 53, 379–389. [PubMed]
  • Choler P., Michalet R., & Callaway R. M. (2001). Facilitation and competition on gradients in alpine plant communities. Ecology, 82, 3295–3308.
  • Collins C. G., Carey C. J., Aronson E. L., Kopp C. W., & Diez J. M. (2016). Direct and indirect effects of native range expansion on soil microbial community structure and function. Journal of Ecology, 104, 1271–1283.
  • Cui H.‐J., Wang G.‐X., Yang Y., Yang Y., Chang R.‐Y., & Ran F. (2016). Soil microbial community composition and its driving factors in alpine grasslands along a mountain elevational gradient. Journal of Mountain Science, 13, 1013–1023.
  • Curtis J. T., & McIntosh R. P. (1951). An upland forest continuum in the prairie‐forest border region of Wisconsin. Ecology, 32, 476–496.
  • Edgar R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics, 26, 2460–2461. [PubMed]
  • Edgar R. C. (2013). UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nature Methods, 10, 996–998. [PubMed]
  • Eisenhauer N., Dobies T., Cesarz S., Hobbie S. E., Meyer R. J., Worm K., & Reich P. B. (2013). Plant diversity effects on soil food webs are stronger than those of elevated CO2 and N deposition in a long‐term grassland experiment. Proceedings of the National Academy of Sciences of the United States of America, 110, 6889. [PubMed]
  • Fierer N., Leff J. W., Adams B. J., Nielsen U. N., Bates S. T., Lauber C. L., … Caporaso J. G. (2012). Cross‐biome metagenomic analyses of soil microbial communities and their functional attributes. Proceedings of the National Academy of Sciences, 109, 21390–21395. [PubMed]
  • Guozhen Du D. Z., & Wang G. (1989). Technology of interspecies association of herbaceous plant community. Journal of Ecology, 8, 59–61.
  • Heijden M. G., Bakker R., Verwaal J., Scheublin T. R., Rutten M., Van Logtestijn R., & Staehelin C. (2006). Symbiotic bacteria as a determinant of plant community structure and plant productivity in dune grassland. FEMS microbiology ecology, 56, 178–187. [PubMed]
  • Hobbie S. E., Nadelhoffer K. J., & Högberg P. (2002). A synthesis: The role of nutrients as constraints on carbon balances in boreal and arctic regions. Plant and Soil, 242, 163–170.
  • Iii F. S. C. (2003). The mineral nutrition of wild plants. Annual Review of Ecology & Systematics, 11, 233–260.
  • Jasper D. A., Abbott L. K., & Robson A. D. (1989). Soil disturbance reduces the infectivity of external hyphae of vesicular‐arbuscular mycorrhizal fungi. New Phytologist, 112, 93–99.
  • Jhp K., Woomer P., George T., Fmdes M., Karanja N. K., & Giller K. E. (1997). Agricultural intensification, soil biodiversity and ecosystem function in the tropics: The role of nitrogen‐fixing bacteria. Applied Soil Ecology, 6, 55–76.
  • Jin H., Zhao Y., Zhao W., Yin H., Qin L., Cui K., … Huang X. (2015). Elevational changes in species diversity and similarity analysis of a Rhododendron chrysanthum community on Changbai Mountain. Acta Ecologica Sinica, 35, 125–133.
  • Johnson J. L., & Temple K. L. (1964). Some variables affecting the measurement of “catalase activity” in soil. Soil Science Society of America Journal, 28, 207–209.
  • Joost S., Bonin A., Bruford M., Despres L., Conord C., Erhardt G., & Taberlet P. (2007). A spatial analysis method (SAM) to detect candidate loci for selection: Towards a landscape genomics approach to adaptation. Molecular Ecology, 16, 3955–3969. [PubMed]
  • Kent M. C. P., Coker P. (1992). Vegetation description and analysis: A practical approach. New York, NY, USA: John Wiley and Sons.
  • Kikvidze Z., Pugnaire F. I., Brooker R. W., Choler P., Lortie C. J., Michalet R., & Callaway R. M. (2005). Linking patterns and processes in alpine plant communities: A global study. Ecology, 86, 1395–1400.
  • Klimeš L. (2006). Alpine plant life. Functional plant ecology of high mountain ecosystems by C. Körner. Folia Geobotanica, 41, 454–455.
  • Klose S., & Tabatabai M. A. (2000). Urease activity of microbial biomass in soils as affected by cropping systems. Biology & Fertility of Soils, 31, 191–199.
  • Kudo G. (1993). Relationships between flowering time and fruit set of the entomophilous alpine shrub, Rhododendron aureum (Ericaceae), inhabiting snow patches. American Journal of Botany, 80, 1300–1304.
  • Lamprecht H. (1989). Silviculture in the tropics. Tropical forest ecosystems and their tree species—possibilities and methods for their long‐term utilization. Eschborn: Gtz.
  • Lange M., Eisenhauer N., Sierra C. A., Bessler H., Engels C., Griffiths R. I., … Scheu S. (2015). Plant diversity increases soil microbial activity and soil carbon storage. Nature Communications, 6, 6707. [PubMed]
  • Levy‐Booth D. J., Prescott C. E., & Grayston S. J. (2014). Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems. Soil Biology & Biochemistry, 75, 11–25.
  • Liang Y., Pan F., He X., Chen X., & Su Y. (2016). Effect of vegetation types on soil arbuscular mycorrhizal fungi and nitrogen‐fixing bacterial communities in a karst region. Environmental Science & Pollution Research, 23, 18482. [PubMed]
  • Liu L., Gundersen P., Zhang W., Zhang T., Chen H., & Mo J. (2015). Effects of nitrogen and phosphorus additions on soil microbial biomass and community structure in two reforested tropical forests. Scientific Reports, 5, 14378. [PubMed]
  • Liu J., Sui Y., Yu Z., Shi Y., Chu H., Jin J., … Wang G. (2014). High throughput sequencing analysis of biogeographical distribution of bacterial communities in the black soils of northeast China. Soil Biology & Biochemistry, 70, 113–122.
  • Liu Y. F., Xing M., Zhao W., Fan R. J., Luo S., & Chen X. (2012). Genetic diversity analysis of Rhododendron aureum Georgi (Ericaceae) located on Changbai Mountain using ISSR and RAPD markers. Plant Systematics and Evolution, 298, 921–930.
  • Ludwig J. A., & Reynolds J. F. (1988). Statistical ecology‐a primer on methods and computing. New York, NY, USA: Wiley.
  • Maestre F. T., Callaway R. M., Valladares F., & Lortie C. J. (2009). Refining the stress‐gradient hypothesis for competition and facilitation in plant communities. Journal of Ecology, 97, 199–205.
  • Maidak B. L., Olsen G. J., Larsen N., Overbeek R., Mccaughey M. J., & Woese C. R. (1996). The ribosomal database project (RDP). Nucleic Acids Research, 24, 82. [PubMed]
  • Mao Y., Yannarell A. C., Davis S. C., & Mackie R. I. (2013). Impact of different bioenergy crops on N‐cycling bacterial and archaeal communities in soil. Environmental Microbiology, 15, 928–942. [PubMed]
  • Mendes R., Kruijt M., De Bruijn I., Dekkers E., Van der Voort M., Schneider J. H., … Bakker P. A. (2011). Deciphering the rhizosphere microbiome for disease‐suppressive bacteria. Science, 332, 1097–1100. [PubMed]
  • Kikvidze Z. (2002). Facilitation and competition in alpine plant communities. Global Environmental Research, 6, 53‐58.
  • Mosca E., Gonzálezmartínez S. C., & Neale D. B. (2014). Environmental versus geographical determinants of genetic structure in two subalpine conifers. New Phytologist, 201, 180–192. [PubMed]
  • Murphy J., & Riley J. P. (1962). A modified single solution method for the determination of phosphate in natural waters. Analytica Chimica Acta, 27, 31–36.
  • Naito A. T., & Cairns D. M. (2011). Patterns and processes of global shrub expansion. Progress in Physical Geography, 35, 423–442.
  • Nannipieri P., Giagnoni L., Renella G., Puglisi E., Ceccanti B., Masciandaro G., … Marinari S. (2012). Soil enzymology: Classical and molecular approaches. Biology and Fertility of Soils, 48, 743–762.
  • Oh Y. M., Kim M., Lee‐Cruz L., Lai‐Hoe A., Go R., Ainuddin N., … Adams J. M. (2012). Distinctive bacterial communities in the rhizoplane of four tropical tree species. Microbial Ecology, 64, 1018–1027. [PubMed]
  • Oksanen B. J., Kindt R., Legendre P., & O'Hara B. (2010). vegan: Community Ecology Package. R package version 1.8‐6, URL
  • Olofsson J. (2004). Positive and negative plant–plant interactions in two contrasting arctic‐alpine plant communities arctic. Antarctic, and Alpine Research, 36, 464–467.
  • Olsen S. R., Cole C. V., Watanabe F. S., Dean L. A., Olsen S. R., Cole C. V., … Dean L. A. (1954). Estimation of available phosphorus in soils by extraction with sodium bicarbonate. Washington DC, USA: United States Department Of Agriculture.
  • Park S., Ku Y. K., Mi J. S., Kim D. Y., Ji E. Y., Lee K. M., … Kim H. M. (2006). Principal component analysis and discriminant analysis (PCA–DA) for discriminating profiles of terminal restriction fragment length polymorphism (T‐RFLP) in soil bacterial communities. Soil Biology & Biochemistry, 38, 2344–2349.
  • Pielou E. C. (1969). An introduction to mathematical ecology. New York, NY, USA: Wiley Interscience. John Wiley & Sons.
  • Putten W. H., Bardgett R. D., Bever J. D., Bezemer T. M., Casper B. B., Fukami T., … Schweitzer J. A. (2013). Plant–soil feedbacks: The past, the present and future challenges. Journal of Ecology, 101, 265–276.
  • Putten W. H. V. D., Dijk C. V., & Peters B. A. M. (1993). Plant‐specific soil‐borne diseases contribute to succession in foredune vegetation. Nature, 362, 53–56.
  • Raab T. K., Lipson D. A., & Monson R. K. (1999). Soil amino acid utilization among species of the Cyperaceae: Plant and soil processes. Ecology, 80, 2408–2419.
  • Reid A.M., Lamarque Laurent J., & Lortie C. J. (2010). A systematic review of the recent ecological literature on cushion plants: Champions of plant facilitation. Web Ecology, 10, 44–49.
  • Schinner F., & Mersi W. V. (1990). Xylanase‐, CM‐cellulase‐ and invertase activity in soil: An improved method. Soil Biology & Biochemistry, 22, 511–515.
  • Schlatter D. C., Bakker M. G., Bradeen J. M., & Kinkel L. L. (2015). Plant community richness and microbial interactions structure bacterial communities in soil. Ecology, 96, 134–142. [PubMed]
  • Schob C., Michalet R., Cavieres L. A., Pugnaire F. I., Brooker R. W., Butterfield B. J., … Callaway R. M. (2014). A global analysis of bidirectional interactions in alpine plant communities shows facilitators experiencing strong reciprocal fitness costs. New Phytologist, 202, 95–105. [PubMed]
  • Shen C., Xiong J., Zhang H., Feng Y., Lin X., Li X., … Chu H. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biology & Biochemistry, 57, 204–211.
  • Singh D., Lee‐Cruz L., Kim W. S., Kerfahi D., Chun J. H., & Adams J. M. (2014). Strong elevational trends in soil bacterial community composition on Mt. Halla, South Korea. Soil Biology & Biochemistry, 68, 140–149.
  • Spohn M. (2015). Microbial respiration per unit microbial biomass depends on soil litter carbon‐to‐nitrogen ratio. Biogeosciences, 12, 817–823.
  • Sturm M., Schimel J., Michaelson G., Welker J. M., Oberbauer S. F., Liston G. E., … Romanovsky V. E. (2005). Winter biological processes could help convert arctic tundra to shrubland. BioScience, 55, 17–26.
  • Syndonia B. H. M., Shaver G. R., & Stuart C. F. (2002). Primary and secondary stem growth in arctic shrubs: Implications for community response to environmental change. Journal of Ecology, 90, 251–267.
  • Thakur M. P., Milcu A., Manning P., Niklaus P. A., Roscher C., Power S., … Ai F. (2015). Plant diversity drives soil microbial biomass carbon in grasslands irrespective of global environmental change factors. Global Change Biology, 21, 4076. [PubMed]
  • Uroz S., Tech J. J., Sawaya N. A., Frey‐Klett P., & Leveau J. H. J. (2014). Structure and function of bacterial communities in ageing soils: insights from the Mendocino ecological staircase. Soil Biology & Biochemistry, 69, 265–274.
  • van der Heijden M. G., Bardgett R. D., & van Straalen N. M. (2008). The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecology Letters, 11, 296–310. [PubMed]
  • Van der Putten W. H. (2003). Plant defense belowground and spatiotemporal processes in natural vegetation. Ecology, 84, 2269–2280.
  • Van der Putten W. H., Klironomos J. N., & Wardle D. A. (2007). Microbial ecology of biological invasions. The ISME Journal, 1, 28–37. [PubMed]
  • Vaughan D., & Ord B. G. (1980). An effect of soil organic matter on invertase activity in soil. Soil Biology & Biochemistry, 12, 449–450.
  • Wu Z., Lin W., Li B., Wu L., Fang C., & Zhang Z. (2015). Terminal restriction fragment length polymorphism analysis of soil bacterial communities under different vegetation types in subtropical area. PLoS ONE, 10, e0129397. [PubMed]
  • Xu W., He X., Chen W., & Liu C. (2004). Characteristics and succession rules of vegetation types in Changbai Mountain. Chinese Journal of Ecology, 23, 162–174.
  • Yang X., & Wu G. (1998). The strategy for conservation and sustainable utilization of biodiversity in Changbaishan Biosphere Reserve. Journal of Forestry Research, 9, 217–222.
  • Zak D. R., Holmes W. E., White D. C., Peacock A. D., & Tilman D. (2003). Plant diversity, soil microbial communities, and ecosystem function: Are there any links? Ecology, 84, 2042–2050.
  • Zhang B., Liang C., He H., & Zhang X. (2013). Variations in soil microbial communities and residues along an altitude gradient on the northern slope of Changbai mountain, china. PLoS ONE, 8, e66184. [PubMed]
  • Zhang C., Liu G. B., Xue S., & Xiao L. (2013). Effect of different vegetation types on the rhizosphere soil microbial community structure in the loess plateau of china. Journal of Integrative Agriculture, 12, 2103–2113.
  • Zhao W., Qi X., Lyu J., Yu Z., & Chen X. (2016). Characterization of microbial community structure in rhizosphere soils of Cowskin Azalea (Rhododendron aureum Georgi) on northern slope of Changbai Mountains, China. Chinese Geographical Science, 26, 78–89.
  • Zheng D., Wallin D. O., & Hao Z. (1997). Rates and patterns of landscape change between 1972 and 1988 in the Changbai Mountain area of China and North Korea. Landscape Ecology, 12, 241–254.
  • Zong S., Jiawei X. U., Zhengfang W. U., Qiao L., Wang D., Meng X., … Wang P. (2014). Analysis of the process and impacts of Deyeuxia angustifolia invasion on the Alpine Tundra, Changbai Mountain. Acta Ecologica Sinica, 34, 87–104.

Articles from Ecology and Evolution are provided here courtesy of Wiley-Blackwell