Skip to main content

Root ethylene mediates rhizosphere microbial community reconstruction when chemically detecting cyanide produced by neighbouring plants



Stress-induced hormones are essential for plants to modulate their microbiota and dynamically adjust to the environment. Despite the emphasis of the role of the phytohormone ethylene in the plant physiological response to heterospecific neighbour detection, less is known about how this activated signal mediates focal plant rhizosphere microbiota to enhance plant fitness. Here, using 3 years of peanut (Arachis hypogaea L.), a legume, and cyanide-containing cassava (Manihot esculenta Crantz) intercropping and peanut monocropping field, pot and hydroponic experiments in addition to exogenous ethylene application and soil incubation experiments, we found that ethylene, a cyanide-derived signal, is associated with the chemical identification of neighbouring cassava and the microbial re-assemblage in the peanut rhizosphere.


Ethylene production in peanut roots can be triggered by cyanide production of neighbouring cassava plants. This gaseous signal alters the microbial composition and re-assembles the microbial co-occurrence network of peanut by shifting the abundance of an actinobacterial species, Catenulispora sp., which becomes a keystone in the intercropped peanut rhizosphere. The re-assembled rhizosphere microbiota provide more available nutrients to peanut roots and support seed production.


Our findings suggest that root ethylene acts as a signal with a dual role. It plays a role in perceiving biochemical cues from interspecific neighbours, and also has a regulatory function in mediating the rhizosphere microbial assembly, thereby enhancing focal plant fitness by improving seed production. This discovery provides a promising direction to develop novel intercropping strategies for targeted manipulations of the rhizosphere microbiome through phytohormone signals.

Video abstract.


Rhizosphere microorganisms are a reservoir of additional functions that extend the plant's ability to adapt to various environmental conditions and changes [1, 2]. Plants determine the rhizosphere microbiome depending on the composition of root-secreted metabolites [3, 4]. In natural environments, plants constantly adjust the composition and concentration of root metabolites in response to various kinds of biotic stressors, such as neighbour competition, pathogen infection and herbivore attack [5,6,7]. This may result in allelochemical responses to neighbouring plants, reduced plant susceptibility to pathogen attack or the prevention of grazing by herbivores [8, 9]. However, relatively few studies have sought to understand the mechanisms and effects of stress-induced root-secreted metabolites on rhizosphere microbiota [10, 11].

Among root-secreted metabolites, stress-derived hormones are a class of small bioactive molecules. In addition to regulating plant physiological and morphological responses [12, 13], a growing number of studies have demonstrated that these phytohormones also sculpt the root microbiome [14, 15]. For instance, the secretion of root salicylic acid (SA) is involved in plant neighbour detection and shapes the root microbiome by modulating taxonomic groups of bacteria [9, 10]. Similar functions have been reported for jasmonic acid (JA) and its derivatives [16]. In addition to JA and SA, ethylene (ET), a volatile signal, can easily diffuse through gas- and water-filled pores in the soil [17]. It has been demonstrated that ET can influence nodulation in legume-rhizobia symbioses and arbuscular mycorrhizal root colonization [18,19,20,21,22]. These findings suggest that ET is a chemical cue used to monitor and interact with soil-specific microbial species near growing roots. Gaseous compounds have a wide effective range in soil for the long-distance attraction of bacteria to roots [23]. However, which rhizobacteria are affected by ET and how ET controls the rhizosphere microbial community assembly remains unknown.

When two species co-exist, heterospecific metabolites are important cues for neighbour detection and subsequently trigger complex plant response strategies [5, 9]. Cyanide commonly occurs in over 3000 plant species, including important crop plants, such as maize, wheat and cassava [24, 25]. Exposure to cyanide can shorten plant embryo dormancy and induce ET production in seedlings [26, 27]. In modern intensive agro-ecosystems, cyanide-containing crops are intercropped with legumes (intercropping: a farming practice involving two or more crop species or genotypes growing together and coexisting for a period of time [28, 29]), but very few studies have focused on the chemical linkage between cyanide-containing plants and legumes [28,29,30]. Whether cyanide from cyanide-rich plants act as a chemical cue to influence ET signalling in non-cyanide-containing plants remains to be uncovered.

To explore whether the belowground chemical interactions between interspecific plants triggers legume rhizosphere microbial re-assembly, we grew cyanide-containing cassava plants together with the legume peanut, which are commonly co-cultivated in subtropical areas [31]. We hypothesized that (1) ET emission from peanut roots can be triggered by the production of cyanide by neighbouring cassava; (2) similar to SA and JA, ET, a volatile hormone, can attract specific microbial species and reshape the microbiota of the plant rhizosphere; (3) the reconstructed microbial community may improve peanut fitness. To evaluate these hypotheses, a 3-year field experiment of peanut and cassava intercropping and peanut monocropping systems and corresponding pot and hydroponic experiments were conducted. We find that neighbour cyanide can be a chemical cue to induce peanut root ET emissions. Using 16S rRNA high-throughput sequencing analysis from the field and ethylene application cultures, we demonstrate that ET modulates the abundance of an actinobacterial Catenulispora species, which functioned as a keystone of the re-assembled microbial network in the intercropped peanut rhizosphere. The reshaped microbial community increased the accumulation of available nutrients in the peanut rhizosphere. Our results reveal a novel chemical dialogue between a focal plant’s roots and its microbiota, entering a partnership that improved the focal plant’s fitness when grown with heterospecific plant neighbours.

Results and discussion

Peanut traits and rhizosphere nutrient characteristics when co-cultured with cassava in the field

To cope with the biotic stresses from heterospecific species, focal plants display a myriad of plastic responses to optimize fitness [32, 33]. Here, the yield (seed production) of peanut in the peanut-cassava intercropping system (Fig. 1a) was similar to that in the monocropping system through an increase in the number of aboveground branches (P = 0.004), the number of pods per plant and the full fruit rate (P < 0.001) (Table 1). However, to gain these advantages, peanut reduced its biomass (P < 0.001), specifically the aboveground biomass (plant height, P < 0.001; and the ratio of aboveground biomass to belowground biomass, P < 0.001) (Table 1). In intercropped peanut individuals, seed production determines the reproductive capacity of the offspring [33, 34]. When peanut was co-cultured with neighbouring cassava, it invested more in seeds than in aboveground tissues. Such physiological alteration in the focal plant community parallels the co-existence strategy for plants with a short stature in Evolutionary Game Theory [35, 36].

Fig. 1

Diagram of intercropping and monocropping sampling sites in the field and pot experiments. a The peanut and cassava intercropping system (i), and the peanut monocropping system (ii) in the field. The intercropping combination included a 2.0-m peanut strip (five rows of peanut, with a 0.4 m interrow distance) and a 0.4-m cassava strip (one row of cassava). The interplant distance within the same column was 0.2 m. In the monocropping system, the interrow distance was 0.4 m, and the interplant distance within the same column was 0.2 m. CRi, BSi and PRi represent soils from the cassava rhizosphere, the bulk soil and the peanut rhizosphere in the intercropping system; PRm and BSm represent soils from the peanut rhizosphere and the corresponding bulk soil in the monocropping system. b Schematic of the intercropping and monocropping pot experiments. Pairs of plants were grown together in a pot separated into compartments with a wire mesh screen filter (red dotted line). Treatment I, cassava monocropping; treatment II, peanut and cassava intercropping, treatment III, peanut monocropping; treatment IV, peanuts were cultured with the addition of cyanide dilution; and control, peanuts were cultured with the addition of water. The rhizosphere and bulk soil collected from the corresponding coloured sites are indicated by the site name plus lowercase “s”; the peanut plants collected from the sites are indicated by the site name plus lowercase “p

Table 1 Peanut growth and production yield in the two cropping systems

In addition to plant physiological changes, significant changes of nutrients in the rhizosphere also attracted our attention (Table S1 in Additional file 1). The rhizosphere is the main micro-domain from which plants acquire resources [37]. The overyielding of peanut individuals did not cause a decline of nutrient supply in the rhizosphere. In contrast, intercropped peanut rhizospheres showed improved soil physiochemical properties, including soil organic carbon (SOC), total nitrogen (TN), ammonia nitrogen (NH4+-N), nitrate nitrogen (NO3-N), total phosphorus (TP) and available phosphorus (AP), compared with monocropped peanut rhizospheres (P < 0.05) (Table S1 in Additional file 1). Since soil microbes play central roles in rhizosphere nutrient release and conversion [37], the differences in nutrient supply between inter- and monocropping peanut rhizospheres may imply an alteration of the microbial communities in the two systems.

Neighbouring cyanide triggers peanut root ethylene production

In the intercropping system, peanut individuals allocated more resources to seed production. This physiological response is always coordinated by phytohormonal signalling [38, 39]. By monitoring phytohormone levels in peanut belowground and aboveground tissues (Fig. 1b), we found that peanut if intercropped with cassava (Pip) reduced the belowground zeatin and indole-3-acetic acid (IAA) concentrations and aboveground jasmonic acid (JA) concentrations while stimulating aboveground gibberellin (GA) and root 1-aminocyclopropane-1-carboxylic acid (ACC) production (P < 0.05) (Fig. S1 in Additional file 1, Fig. 2a). Similar responses could be mimicked by exogenous cyanide addition (CNp in treatment IV) (Fig. S1 in Additional file 1, Fig. 2a). Unlike other food crops, the levels of cyanide production in cassava is depended on plant cultivars and age, with an average of 577 mg kg−1 f.wt in the outer peel (2–5 mm) of the roots [40]. Studies have demonstrated that cyanide can rapidly cross membranes and acts as a regulator to induce ET by modifying ACC synthase and oxidase [26, 41]. However, it has never been reported to participate in natural plant-plant chemical recognition. Here, we found that the soil cyanide concentration formed a steep gradient from the cassava rhizosphere (Cis) to the peanut rhizosphere (Pis) when the two species were co-cultured, and the cyanide content in the peanut rhizosphere (Pis) was still four times higher than that in the peanut rhizosphere of the monocultured system (Pms and PBms, P < 0.05) (Fig. 2b). However, cyanide concentration in peanut tissue was not affected (Fig. 2c). Along with soil cyanide concentration, root ACC in intercropped peanut (Pip) was two times that of the control (Cp) and monocropped peanut (Pmp). Exogenous cyanide addition proved to increase ACC concentration of peanut root (CNp) to a level similar to that of intercropped peanut (Pip) (Fig 2a). Cyanide from neighbouring cassava positively affected the belowground ACC production by peanut. 1-Aminocyclopropane-1-carboxylic acid (ACC) is the direct precursor of ethylene (ET) in plants [42, 43]. Higher ACC in peanut roots corresponded to higher ET production (Fig. 2d).

Fig. 2

Effects of soil cyanide on root ethylene synthesis and release. a The 1-aminocyclopropane-1-carboxylic acid (ACC, ethylene precursor) concentration in aboveground and belowground peanut tissues in different treatments. b Cyanide concentration in cassava and peanut rhizosphere soils and corresponding bulk soils. c Cyanide concentration in peanut roots. d Ethylene (ET) production in peanut roots in different treatments. e qRT-PCR detection of ACS gene (AhACS1 and AhACS2) expression in peanut roots in different treatments. f qRT-PCR detection of ACO gene (AhACO1 and AhACO2) expression in peanut root in different treatments. Treatment I (cassava monocropping), Cms and CBms represent cassava rhizosphere and bulk soil, respectively. Treatment II (cassava and peanut intercropping), Cis, PBis and Pis represent soil from the cassava rhizosphere, the bulk soil and peanut rhizosphere; Pip represents peanut plants from treatment II. Treatment III (peanut monocropping), Pms and PBms represent soil from the peanut rhizosphere and bulk; Pmp represents peanuts from treatment III. Treatment IV (peanut monocropping with exogenous CN addition), CNp represents peanuts from treatment IV. Control (peanut with water addition), Cp represents peanuts from control. Data of soil and plant cyanide are mean values + SD for triplicate; data of plant ACC concentration, ET production and gene expression are values + SD for 9 replicates. Error bars with different letters indicate a significant difference according to one-way analysis of variance (ANOVA) with Tukey’s HSD test (P < 0.05)

Acera et al. [44] reported that some soil microorganisms have the potential to transform cyanide to amino acids (such as β-cyanoalanine), which may be further metabolized to ET [13]. To eliminate the interference of soil microorganisms on ET production, quantitative real-time PCR (qRT-PCR) was used to track gene expression related to 1-aminocyclopropane-1-carboxylic acid synthase (ACS) and 1-aminocyclopropane-1-carboxlic acid oxidase (ACO). The expression of the AhACS1, AhACO1 and AhACO2 transcripts, which probably encode peanut ACS and ACO proteins, was significantly higher in the intercropped and cyanide-treated peanut roots (Pip and CNp) (Fig. 2e, f). Similar differences were found between the cyanide addition treatments and the control in the hydroponic experiments (Fig. S2 in Additional file 1). When plant species co-exist, heterospecific compounds become chemical cues for focal plants to recognize competitors [45]. Here, cyanide from neighbouring cassava was shown to activate the peanut ET production in response to the neighbour belowground.

As the central hormone regulator, ethylene (ET) triggers a range of physiological adaptations [46,47,48]. Higher ET causes plants to invest more in reproduction under harsh conditions [13]. The presence of a heterospecific neighbour maximized peanut seed production at the expense of plant biomass (Table 1). This optimized performance is the objective in agriculture but suppressed naturally [35, 36].

Ethylene chemically regulates rhizosphere microbial diversity and composition

Conventional agricultural wisdom might suggest that the increased peanut yield was associated with high nutrient consumption from the soil [49, 50]; however, in contrast to our expectations, available nutrients, including ammonia nitrogen (NH4+-N) and available phosphorus (AP), were higher in the intercropped peanut rhizosphere than in the monocropped peanut rhizosphere (Table S1 in Additional file 1). Microorganisms are the decomposers to mineralize organic matter in soil, and their activity, diversity and composition determine the amount of available nutrients to plants [37]. Using 16S rRNA high-throughput sequencing, we found that the lowest microbial Shannon and Chao 1 indices were observed in the cassava rhizosphere (CRi) (P < 0.05) (Fig. 3a, b), which may be attributed to the large amounts of cyanide accumulating around cassava roots (Fig. 2b). A high concentration of cyanide has been reported to be toxic to microbes [25, 40]. In contrast, the peanut rhizosphere in the intercropping system (PRi) showed the highest microbial diversity (Shannon and Chao 1 indices, P < 0.05) (Fig. 3a, b). This finding differs from conventional views, because rhizosphere microbial communities are generally less diverse than those in bulk soil [51, 52]. In the intercropping systems, however, root microbiota often show a higher diversity [53, 54]. Depending on microbial substrate preferences, secondary metabolites of plants shape the rhizosphere microbial community [11]. Here, neighbouring cassava increased ethylene production in the belowground portion of peanut. In order to find out if increasing ethylene (ET) concentration directly altered rhizobacterial α-diversity, we applied exogenous ethylene to rhizosphere soil of monocropped peanut in an incubation experiment. We found that compared with an untreated control, ethylene treatment, at the 0.1–0.2 mM concentration level, significantly increased rhizosphere bacterial α-diversity (Shannon and Chao 1 indices, Fig. 3c, d).

Fig. 3

α-Diversities of soil microbiota in the field and ethylene (ET) addition incubation system. a, b Shannon and Chao 1 indices of the microbiota of rhizosphere soil from cassava and peanut and the corresponding bulk soils in the inter- (blue columns) and monocropping (pink columns) field systems. c, d Shannon and Chao 1 indices of the microbiota in different ET concentration treatments. e, f Phylum-level distribution of microbial composition in the field and in the ET addition treatments. Error bars with different letters indicate a significant difference according to one-way analysis of variance (ANOVA) with Tukey’s HSD test (P < 0.05)

Overall, the phyla Acidobacteria, Chloroflexi, Actinobacteria, Alphaproteobacteria, Planctomycetes, Betaproteobacteria and Gammaproteobacteria were present at high relative abundances (average relative abundance > 5%) in the field (Fig. 3e). Interestingly, only the abundance of Actinobacteria and Acidobacteria differed significantly between peanut rhizospheres in intercropping (PRi) and monocropping (PRm) systems (P < 0.001). The same effect could be mimicked by exogenous ethylene addition: as actinobacterial abundance increased, acidobacterial abundance decreased with ethylene addition compared to the control (P < 0.05) (Fig. 3f). This result is in line with Lebeis [10], who showed that the phytohormone salicylic acid (SA) affects the composition of the root microbiome of Arabidopsis thaliana, leading to reduced representation of Acidobacteria and enrichment of Actinobacteria. Here, we paid special attention to the influence of ethylene on the composition of the rhizosphere microbial communities. Similar to SA, ethylene has the potential to alter microbial community assembly of plant roots, which subsequently may feedback on plant phenotypic traits [13]. Such a mechanism could be exploited to improve sustainable plant production systems.

To further elucidate differences in community turnover between the soil and rhizosphere microbiota in the field, we analysed β-diversity using Bray-Curtis distances in principal coordinate analysis (PCoA) based on weighted UniFrac metrics (Fig. 4). In the peanut monocropping system, the rhizosphere (PRm) and bulk soil (BSm) were separated (PANOSIM < 0.01) (Fig. 4a), indicating a clear plant rhizosphere effect on the bacterial community assembly [4, 28, 37]. However, in the presence of neighbouring cassava, the cyanide-influenced bulk soil-associated and peanut rhizosphere bacteria (BSi and PRi) appeared to converge, but were still well separated from those of the monocultures (PRm and BSm) and cassava rhizosphere (CRi) (PANOSIM < 0.01) (Fig. 4a). Similar patterns of clustering by planting system were confirmed by hierarchical clustering (Fig. S3 in Additional file 1). Differences in community composition between PRi and PRm were mainly due to the reduction of Acidobacteria (Fig. 4b). This was consistent with the results of the ethylene addition experiment, where concentrations of 0.1 mM and 0.2 mM ET led to a replacement of Acidobacteria by Actinobacteria clustered away from the control along the first coordinate axis (PANOSIM < 0.01) (Fig. 4c, d). Fu [55] confirmed that many species of Actinobacteria have the potential to utilize ethylene (ET) as a source for their carbon metabolism. This may explain the increasing relative abundance of Actinobacteria in intercropped and ethylene-treated monocropped peanut rhizospheres.

Fig. 4

Principal coordinate analysis (PCoA) of the weighted UniFrac values among the field and ethylene addition incubation samples based on Bray-Curtis distances. a Bacterial community composition in the field by PCoA. The microbiota is clearly clustered into four groups that are marked with four coloured dotted ellipses. b The dominant OTU (relative abundance > 0.1%) scores of the bacterial community in the field according to the PCoA. c Bacterial community composition in the ET addition treatments by PCoA. d The dominant OTU (relative abundance > 0.1%) scores of the bacterial community in the ET addition treatments according to the PCoA. The arrows indicate the centroid of the constrained factor. Circle sizes represent the relative abundances of bacterial OTUs, and colours were assigned to different phyla. CRi, cassava rhizosphere soil in the interspecific experiment; BSi, bulk soil 20 cm away from the peanut and cassava plants; PRi, peanut rhizosphere soil in the interspecific experiment; PRm, peanut rhizosphere soil in the intraspecific experiment; BSm, bulk soil 20 cm away from peanut plants

Ethylene regulates the rhizosphere network through effects on specific microbial taxa

Microorganisms do not exist in isolation but interact with others in nature by positive, negative or neutral ecological interactions [56]. Microbial co-occurrence and network analysis provide a promising approach to investigate these various types of interactions in microbial communities [56, 57]. In this study, we constructed co-occurrence networks using random matrix theory (RMT) to determine the differences in bacterial assemblages in rhizospheres of the different plant species and bulk soil [58, 59]. All values of the calculated modularity index were larger than 0.4 (Table 2), suggesting typical module structures [60]. Overall, crop planting showed a marked effect on the soil microbial network: average path distance (GD), the average clustering coefficient (avgCC) and modularity of the empirical networks were higher than those of corresponding, identically sized random networks (Table 2). Plant rhizosphere network connectivity and complexity were first well described as properties of rhizosphere bacterial assemblages by Shi et al. [58]. Here, we found that plant rhizosphere assemblages (in CRi, PRi and PRm) formed more connected and more complex networks with fewer nodes but more connections (edges) between nodes compared with the bulk soil. Accordingly, the densities of the connections (graph density) in the rhizosphere increased in the plant rhizosphere compared with the corresponding bulk soil. In the networks, the ratio of positive to negative connections (edges) was 1.2 in PRi and 2 in PRm (Table 2), indicating that more competitive or inhibitive connections in the intercropped peanut rhizosphere.

Table 2 Topological properties of the empirical molecular ecological networks of microbial communities in treatments

Microbial communities can harbour keystone taxa, whose removal can cause a dramatic shift in microbiome structure and functioning [57, 60]. Keystone taxa in network analysis can be computationally identified as hubs with a high within-module degree Zi (Zi ≥ 6.2 indicates that the nodes are “well-connected” to other nodes in the module) [59, 61,62,63]. In this study, bacterial module hubs affiliated with Dokdonella (OTU1, H1 belonging to Gammaproteobacteria), Catenulispora (OTU235, H2 belonging to Actinobacteria) and Pseudolabrys (OTU73, H3, belonging to Alphaproteobacteria) were respectively identified in the cassava, intercropped peanut and monocropped peanut rhizosphere samples (CRi, PRi and PRm) (Table 3, Fig. S4 in Additional file 1). No hub was found in the bulk soil (BSi and BSm) (Fig. 5). In the peanut rhizosphere, Catenulispora sp. replaced Pseudolabrys sp. as a keystone taxon in intercropping and showed more negative correlations to other taxa in the network (Table 3, Fig. 5). Catenulispora belongs to Actinobacteria, which is a phylum that serves as a source of highly diverse antibiotics [64, 65]. Catenulispora is known to produce the antibiotics cacibiocin A and B, which are inhibitors of bacterial type II topoisomerases, such as DNA gyrase and DNA topoisomerase IV in Acidobacteria [66]. To a certain extent, this can explain the decrease in the acidobacterial abundance in the rhizosphere of the peanut intercropping system (Fig. 3e, f, Fig. 4b, d).

Table 3 Nodes identified as hubs in the networks in the inter- and monocropping systems
Fig. 5

Plant rhizosphere and the corresponding bulk soil networks in the inter- and monocropping systems. Networks represent random matrix theory co-occurrence models derived from 8 biological replicates at each site, where nodes represent OTUs, and the edges between the nodes indicate significant correlations. A green edge indicates negative covariation between two individual nodes, while a red edge indicates positive covariation. The colours of the nodes indicate major phyla (subphyla for Proteobacteria). In each panel, the size of each node is proportional to the number of connections (i.e. node degree). Nodes marked “H1, H2 and H3” represent the identified module hubs as follows: H1, OTU1 Dokdonella; H2, OTU235 Catenulispora; and H3, OTU73 Pseudolabrys

We further used indicator species analysis to identify operational taxonomic units (OTUs) specifically associated with ethylene addition and found that actinobacterial indicators appeared in the 0.1 mM and 0.2 mM ethylene addition treatments (Table S2 in Additional file 1). Indicator species analysis is a method determined by Dufrene-Legendre [67] to identify species preference for environmental changes. In the control (no ethylene addition) and in high ethylene addition (0.5 mM) treatments, no actinobacterial indicator was found (Table S2 in Additional file 1). Interestingly, among actinobacterial indicators, OTU235 (Catenulispora sp.), which was the keystone in the network of the intercropped peanut rhizosphere, again turned out as indicator taxon in the 0.2 mM ethylene (ET) treatment (Table S2 in Additional file 1). Its relative abundance was significantly dependent on ethylene concentration (Fig. 6a, b). To date, there is no evidence of direct interactions between Catenulispora sp. and plant phytohormones, but species belonging to Catenulisporaceae have been reported to act as plant-associated bacteria strengthening root colonization by producing the enzyme ACC deaminase, which controls the endogenous ethylene (ET) levels [68, 69]. In addition to the increase in Catenulispora sp. relative abundance, the soil ammonium nitrogen (NH4+-N) and available phosphorus (AP) concentrations increased (r = 0.816 and r = 0.733, P < 0.01) (Fig. 6c, d). This may be attributed to the positive effect of the re-assembled microbial community on soil phosphorus (index of acid phosphatase activity, P < 0.05) and organic nitrogen (indices of urease activity and L-glutamate activity, P < 0.05) mineralization (Fig. S5 in Additional file 1).

Fig. 6

The variance of the relative abundance of the keystone and the correlation between its abundance and soil available nutrient content in ethylene addition treatments. a Layout of the exogenous ethylene addition test. Thirty grams of soil from the PRi treatment were placed in the bottom of a sterilized bottle. Then, 200 mL of 0.1 mM, 0.2 mM and 0.5 mM standard ethylene gas (ethylene and air mixed) were injected into the glass bottle (V = 100 mL) to fully replace the air. For the control (0 mM), the injection gas was replaced with air. After injection, the tubes were immediately sealed and incubated at 26 °C for 7 days for microbial community detection. b Dynamics of relative abundance of OTU235 (module hub H2 in RPi network) in different ethylene level treatments. Columns with different colours indicate significant differences (P < 0.05) according to one-way analysis of variance (ANOVA) followed by Tukey's HSD test (P < 0.05). c, d Significant relationship between OTU235 abundance and available nutrients (NH4+-N and AP) in the ethylene addition treatments (P < 0.01)

Influence of soil properties, plant traits and the rhizosphere microbial community on peanut production

To investigate the potentially important predictors of peanut seed production, we conducted random forest modelling with the soil properties (including TN, SOC, AP and pH), plant traits (including ethylene emission and plant biomass), and bacterial community (including β-diversity, relative abundance of keystone taxa and average network connectivity) [70]. The model indicated that the most important predictor of seed production was average network connectivity, followed by bacterial β-diversity, ethylene (ET) emission, relative abundance of keystone taxa, SOC, TN, plant biomass and AP (P < 0.05) (Fig. 7a). Soil pH had no influence on plant seed production (P = 0.33) (Fig. 7a). We then used structural equation modelling (SEM) to identify the potential direct and indirect effects of soil properties, plant traits and the microbial community on seed production (Fig. 7b). Structural equation modelling (SEM) is a powerful statistical modelling technique that is widely used in the behavioural sciences and ecological research [71]. It can be viewed as a combination of factor analysis and regression or path analysis. In this study, soil properties, ethylene (ET) production and peanut biomass were positively related to peanut seed production (P < 0.01). Although ethylene slightly reduced plant biomass (P < 0.05), it modulated the microbial community structure by affecting keystone species abundance (P < 0.01). The resulting microbial re-assembly showed a direct positive correlation with the soil properties (P < 0.001), suggesting that microbial community reconstruction did benefit nutrient mineralization (Fig. 7b).

Fig. 7

Direct and indirect effects of ethylene on peanut seed production. a Mean predictor importance of factors related to plant fitness based on random forest analysis. b Direct and indirect effects of ethylene emission on seed production using structural equation modelling (SEM). The significance levels of each predictor in the random forest analysis are as follows: *P < 0.05 and **P < 0.01. The blue and red arrows in the structural equation model indicate positive and negative relationships, respectively, and dotted arrows represent nonsignificant paths (P > 0.05). Numbers adjacent to arrows are standardized path coefficients, and the path width indicates the strength of significant standardized path coefficients. The first principal coordinate (PCoA1, which explained 59.3% of the variation) is used to represent the composition of the bacterial β-diversity. ET, ethylene emission; Bio, plant biomass; KRA, keystone relative abundance; avgK, average network connectivity; β-div, bacterial β-diversity; SOC, soil organic carbon; TN, soil total nitrogen; AP, soil available phosphorus

Keystone species serve as gatekeepers in the ecological functions of the bacterial community [59, 60]. When detecting cyanide from a neighbour, root ethylene from the peanut directly affects the abundance of the gatekeeper, which leads to a dramatic shift in the community composition of the peanut rhizosphere microbiota, and likely further increases seed production by enhancing the accumulation of available nutrients in the rhizosphere (Fig. 8). During the long-term co-evolution of plants and microorganisms, plants have developed mechanisms to process signals, thus triggering an optimized response to maximize resistance at minimal costs [13]. Signal integration between plants and microbes is a result of killing two birds with one stone, stimulating plant physiological responses and causing plants and microbes to form a coherent unit to shift plant responses towards or away from adaptation to specific situations [15]. In the past, these views were mostly speculative [13, 15]. Here, using a series of experiments, we demonstrated that ethylene not only participates in the recognition of heterospecific plant neighbours, but also helps plants to regulate rhizosphere microbial communities to improve plant fitness in the process of plant interspecific interactions.

Fig. 8

The overview of the mechanism by which peanut mediates rhizosphere microbiota to increase available nutrients in the cassava-peanut intercropping system. Cyanide exposure belowground induces ethylene production in peanut roots, and ethylene mediates actinobacterial species to reshape rhizosphere microbiota, which promote peanut seed production by increasing available nutrient content


Cyanide produced by neighbouring plants is a chemical cue that increases root ethylene emission from the legume peanut. This phytohormone signal led to an increased relative abundance of specific actinobacterial species, thereby reshaping the whole rhizosphere microbial network. This led to an improved availability of essential nutrients and a shift in plant allocation towards increased seed production at the expense of plant biomass. Our results reveal a novel function of a stress-derived plant signal to re-shape microbial community assembly in the rhizosphere and establishing strategic plant-microbe partnership enhancing plant fitness in intercropping agro-ecosystems.

Materials and methods

Field experiment site

The field site is located at the Red Soil Ecological Experimental Station of the Chinese Academy of Sciences in Yingtan, Jiangxi Province, China (28° 15′ N, 116° 5′ E). The site is located at an altitude of 79 m a.s.l., with a mean annual temperature of 17.6 °C and a mean annual precipitation of 1795 mm. The total solar radiation is 6514.2 MJ m−2 year−1, and the potential evaporation is 1318 mm. The frost-free period is 262 days. The soil type is classified as acid loamy clay derived from Quaternary red clay (Udic Ferralsols in the Chinese Soil Taxonomy and Ferric Acrisols in the FAO classification system) with a pH (water) of 4.9. In the top 25 cm of the soil profile, the soil organic carbon (SOC) content was 10.23 g kg−1, the total nitrogen (TN) content was 0.90 g kg−1, the available phosphorus (AP) content was 34.15 mg kg−1, and the available potassium (AK) content was 235.11 mg kg−1.

Inter- and monocropping field experimental design

The experiment was conducted from 2013 to 2016 and included the following two treatments: (i) PC, a peanut (Arachis hypogaea, a legume) and cyanide-containing cassava (Manihot esculenta Crantz) intercropping system, and (ii) PP, a peanut monocropping system (Fig. 1a). The design of the PC treatment (i) included a 2.0-m peanut strip (five rows of peanut, with a 0.4 m interrow distance) and a 0.4 m cassava strip (one row of cassava). The interplant distance within the same row was 0.2 m for peanut and 0.8 m for cassava. In the PP treatment (ii), the interrow and interplant distances were 0.4 m and 0.2 m, respectively, which made the peanut density identical to that in a comparable area in the PC treatment (Fig. 1a). Each treatment was replicated four times. All plots received 120 kg ha−1 nitrogen fertilizer (urea containing 46% N), 90 kg ha−1 phosphorus (calcium superphosphate containing 12.5% P2O5) and 135 kg ha−1 K2O (potassium chloride containing 60.0% K2O). Each individual plot was 10 × 8 m (length × width), and a ridge (with a width of 0.4 m and a height of 0.3 m) separated adjacent plots. Cassava stems (20 cm in length) were transplanted on 15–18 April and harvested on 10–15 November, while the peanuts were sown on 15–20 April and harvested on 15–20 August. All plots were irrigated and weeded during the growing period. The yields of peanut and cassava were determined at harvest in 2016.

Soil and plant sampling in the field

Field soil samples were collected as follows: (1) PC intercropping: peanut rhizosphere soil (PRi), cassava rhizosphere soil (CRi), and bulk soil (BSi) 20 cm away from the peanut and cassava roots; (2) PP monocropping: peanut rhizosphere soil (PRm) and bulk soil (BSm) 20 cm away from peanut roots (Fig. 1a). Soil that tightly adhered to the plant roots was defined as rhizosphere soil [63]. Due to the small amount of rhizosphere soil per plant, we randomly selected six plant roots that were pooled into a single rhizosphere soil sample. In total, 40 soil samples (5 sites × 2 replicates for each plot × 4 plot replicates) were immediately sieved (4 mm) in the laboratory. Any visible living plant material was manually removed. Five grams of each soil sample was stored at − 80 °C for microbial molecular analysis, and the rest was stored at 4 °C for chemical analyses. Standard methods were used to characterize soil chemical properties (see the details in Methods S1 of Additional file 1). An additional sample of 500 kg of soil was collected at a 5–20-cm depth within a 1.5 × 2 m area in the PP treatment after the peanuts were harvested and stored at 20 °C for the soil pot experiment.

Sixteen peanut plants in each treatment (4 biological replicates × 4 plot replicates) were randomly selected to determine the chlorophyll content in situ (from a single leaf at the top of the plant) using SPAD 502 plus (Konica, Tokyo, Japan) and to measure plant height and aboveground branching. Plant biomass was determined after drying at 65 °C to a constant weight. All the soil and plant samples were collected on 29 May 2016 at the peanut flowering stage.

Pot experiment to detect peanut ethylene production

To identify the metabolites involved in interspecific interactions, inter- and monocropping pot treatments (treatments I–III) were performed as shown in Fig. 1b. Each pot was 80 × 40 × 15 cm (length × width × depth). A wire mesh screen (0.05 mm) was placed in the centre of a pot to create two compartments (40 × 40 × 15 cm in length × width × depth). In treatment I, two cassava individuals were planted in the two compartments in a single pot; in treatment II, cassava and peanut individuals were planted in the two compartments; and in treatment III, two peanut individuals were planted in the two compartments (Fig. 1b). Based on the variation in the soil cyanide concentration and peanut phytohormones in the inter- and monocropping pot treatments, an additional exogenous cyanide (CN) treatment (treatment IV) was performed to identify the hormonal response of peanut to soil cyanide concentrations. For treatment IV, 50 mL of 4 mg L−1 CN in water (diluted cyanide standard) was sprayed on the soil surface before 24 h of destructive plant sampling. For the control (C), the solution was replaced with distilled water. Each treatment comprised nine biological replicates.

Plants were germinated separately and transplanted to the pots when their aboveground height reached 10 cm. They were grown at a temperature of 25–26 °C, 65–70% relative humidity, and a 12/12 h light/dark photoperiod with a light intensity of 10,000 lx. Each pot contained 16 kg (dry soil weight) of soil from the peanut monocropping system in the field (PP). After 20 days of transplantation, plants were destructively sampled.

For soil cyanide measurement, 10 g of each soil sample from seven sites (Fig. 1b) in treatments I–III (treatment I: Cms, cassava rhizosphere soil, and CBms, the bulk soil 20 cm away from the cassava roots; treatment II: Cis, cassava rhizosphere soil; Pis, peanut rhizosphere soil; and PBis, the bulk soil 20 cm away from the cassava and peanut roots; and treatment III: Pms, peanut rhizosphere soil, and PBms, the bulk soil 20 cm away from the peanut roots) were collected and determined based on the methods of the U.S. EPA [72]. To detect plant cyanide, 5 g of fresh plant tissue was ground and measured based on the methods of the U.S. EPA [72]. Among nine replicates of each treatment, three replicates were randomly chosen (7 sampling sites × 3 replicates) to collect soil from specific sites.

For peanut phytohormone detection (see the details in Methods S1 of Additional file 1), a total of 36 plant individuals collected from four treatments (treatment II: Pip, intercropped peanut; treatment III: Pmp, monocropped peanut; treatment IV: CNp, 24 h of cyanide stressed peanut; and control: Cp, peanut in the control) were immediately divided into below- and aboveground tissues and then frozen in liquid nitrogen for phytohormone determination, including zeatin, gibberellin (GA3), abscisic acid (ABA), auxin (IAA), salicylic acid (SA), jasmonic acid (JA) and ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC). Each treatment comprised nine biological replicates, including 36 peanut plants (from 4 treatments × 9 replicates) and 21 soil samples (7 sampling sites × 3 replicates).

Measurement of ethylene production from peanut roots

Due to strong variations of root 1-aminocyclopropane-1-carboxylic acid (ACC), we examined ethylene (ET) production and its synthesis by gene expression in peanut roots. We repeated treatments II, III, and IV and the control (C) mentioned in the “Pot experiment to detect peanut ethylene production” section. After culturing for 20 days, a total of 36 peanut roots (4 treatments × 9 replicates) were cut and placed in 15 mL glass vials containing 1 mL 0.6% water agar and rapidly sealed with a gas-proof septum according to Wu [73]. After 4 h of incubation in darkness at 30 °C, 1 mL of the gas was withdrawn from the airspace of each vial using a gas-tight syringe (Focus GC, Thermo, Massachusetts, USA) and injected into a gas chromatograph (Focus GC, Thermo) equipped with a capillary column (CP-CarboPLOT P7, California, USA) and flame-ionization detector for ET determination. The production of ethylene (ET) was calculated on the basis of the fresh weight (f.wt) of the root samples [74]. Moreover, 36 other peanut roots (4 treatments × 9 replicates) were collected and frozen immediately in liquid nitrogen for RNA extraction and ACS/ACO gene expression analysis.

Expression analysis of ACS and ACO transcripts in peanut roots

1-Aminocyclopropane-1-carboxlic acid synthase (ACS) and (1-aminocyclopropane-1-carboxlic acid oxidase (ACO) are key enzymes controlling plant ethylene synthesis. These genes are organ-specific and are differentially regulated by different environmental signals [75]. To date, these genes in peanut have not been reported. Fortunately, based on the ACS/ACO protein sequences in Arabidopsis thaliana (TSA contigs GDKN01000001 - GDKN01102303 in NCBI, we mined three hypothetical and expressed AhACS/AhACO transcripts in allotetraploid peanut transcript data [75]. Total root RNA was extracted with E.Z.N.A. Total RNA Kit I (OMEGA, GA, USA) and reverse transcribed into cDNA with cDNA Synthesis Kit (Thermo Fisher Scientific, MA, USA) following the manufacturer’s instructions. Quantitative real-time PCR (qRT-PCR) was performed in a 20 μL reaction volume using a CFX connect Real-Time System (Thermo Fisher Scientific, MA, USA) and TB GreenTM Premix Ex TaqTM (Takara, Kusatsu, Japan). The peanut Actin gene (Aradu. W2Y55) was used as the internal control, and the relative quantitation of target gene expression among the different experimental conditions was calculated using the comparative Ct method. All the qRT-PCR primers are listed in Table S3 of Additional file 1. All the mRNA data were expressed as percent of the corresponding Actin transcript levels.

Effect of exogenous cyanide on peanut ethylene production in a hydroponic experiment

To eliminate the interference of soil microbiota with rhizosphere ethylene production, a hydroponic culture was set up (Fig. S2a of Additional file 1). Seeds were surface-sterilized with 0.1% HgCl2 for 5 min and washed twice with sterilized-distilled water and then germinated on sterilized water-moistened gauze. Seedlings were cultivated with 1/4 Hoagland’s nutrient solution until the aboveground height reached 10 cm. Plant growth conditions were consistent with those in the pot experiment “Pot experiment to detect peanut ethylene production”. Then, the seedlings were transferred to 4 mg L−1 cyanide-containing solutions (cyanide) and 4 mg L−1 cyanide + 4 mg L−1 CoCl2 (cobalt chloride, an inhibitor of ethylene biosynthesis) solutions (cyanide + CoCl2). For the control (H2O), solutions were replaced with distilled water. Twenty-four and 48 h after incubation, all roots were cut for the measurement of ethylene production and ACS/ACO transcript expression. Each culture was grown in eight replicates for each time point sampling. After incubation for 24 h and 48 h, four replicates were used to detect root ethylene production, and the other four replicates were used for gene expression.

Soil incubation with the addition of exogenous ethylene

To further determine whether rhizosphere ethylene had the potential to mediate soil-specific microbiota, we incubated peanut rhizosphere soil with different concentrations of ethylene. Thirty grams of soil from the field monocropping system (PP) was placed in a sterilized glass bottle (V = 100 mL). Three concentrations, (A) 0.1 mM, (B) 0.2 mM and (C) 0.5 mM ethylene standard, were injected into the bottle to replace the air (Fig. 6a). A control (0 mM) with air injection was established in an identical manner. After injection, the bottles were immediately sealed and incubated in triplicate at 26 °C for 7 days (Fig. 6a). After 7 days, 1 g soil samples were collected for 16S high-throughput sequencing analysis. The rest of the soil was used for the measurements of chemical properties and enzymatic activities (see the details in Methods S1 of Additional file 1). Each treatment was incubated in triplicates.

Soil DNA extraction and 16S rRNA high-throughput sequencing

Soil samples (1.0 g) from field and incubation experiments were extracted using the FastDNA SPIN Kit (MP Biomedical, California, USA) according to the manufacturer’s instructions. The quantity and purity of DNA were examined with a Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Delaware, USA). The V4-V5 region of the bacterial 16S rRNA gene was amplified using the primers 515F and 907R [76]. Each sample was amplified in a 20-μl reaction system, which contained 0.5 μM forward and reverse primers, 1 × Premix Taq DNA polymerase (Takara, Kusatsu, Japan) and 20 ng DNA templates. After an initial denaturation at 95 °C for 3 min, the targeted region was amplified by 20 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, followed by a final extension at 72 °C for 1 min in a thermal cycler (GeneAmp PCR system 2700; Applied Biosystems, New York, USA). Amplicon sequencing libraries were constructed using the MiSeq Reagent Kit v3 according to the manufacturer’s instructions. High-throughput paired-end sequencing was performed on the Illumina MiSeq PE250 platform.

The raw data were screened and trimmed by the QIIME pipeline (version 1.9.0) [77]. To minimize the effects of random sequencing errors, low-quality and ambiguous reads (Phred quality score Q < 25 or sequence shorter than 150 bp) were eliminated. Chimeras were filtered with the UCHIME algorithm in the USEARCH package [78, 79]. High-quality sequences were clustered into operational taxonomic units (OTUs) using UCLUST with a similarity threshold of 97% [78]. The sequences were then phylogenetically assigned to taxonomic classifications using the RDP (Ribosomal Data Project database) classifier and were allocated to different levels [80]. Singletons were removed, and all samples were rarefied to 20,000 sequences per sample for further analysis. The sequencing data were deposited in the European Nucleotide Archive of EMBL under the accession number PRJEB22658.

Soil microbial network construction and keystone identification

Microbial networks were constructed for rhizosphere and bulk soil communities based on OTU relative abundance in the field experiment. Covariations were determined across eight biological replicates to create each network. Random matrix theory (RMT) was used to construct co-occurrence networks by calculating all pairwise Spearman’s rank correlations (P < 0.01). Random networks were generated based on the Maslov and Sneppen model [81]. Global network properties were characterized according to Deng et al. [59]. All network analyses were performed using the Molecular Ecological Network Analyses pipeline ( written in Java and Perl scripts [59, 82, 83]. To reduce network complexity, we only considered bacterial OTUs with an average abundance > 0.1%. The OTUs detected in more than 75% of the samples were retained in the network. Various indices, including the average clustering coefficient (avgCC), average geodesic distance (GD), and size and modularity of the network, were calculated to describe network topologies. Average connectivity (avgK) was calculated to measure the complexity of the network structure [59]. The topological role of each node was determined based on two properties: the within-module connectivity (Zi) and the among-module connectivity (Pi) [61]. All species were sorted into four subcategories on the basis of these simple criteria: peripherals (nodes in the modules with few outside connections, Zi < 2.5 and Pi < 0.62), connectors (nodes that connect modules, Pi > 0.62), module hubs (highly connected nodes within modules, Zi > 2.5) and network hubs (highly connected nodes within the entire network, Zi > 2.5 and Pi > 0.62) [59, 63, 82, 84]. Correlation networks were visualized using Gephi software [85].

Structural equation modelling

Random forest modelling was conducted to quantitatively assess the important predictors of plant productivity [70, 86], including soil properties (total nitrogen, TN; soil organic carbon, SOC; available phosphorus, AP and pH), plant physiology (ethylene production, ET; peanut biomass, Bio), and the bacterial community (bacterial β-diversity, β-div; keystone relative abundance, KRA; average connectivity, avgK). The importance of each factor was evaluated by the decrease in prediction accuracy (that is, an increase in the mean square error (MSE) between observations and predictions) when the data were randomly permuted [86]. This accuracy importance measure was computed for each tree and averaged over the forest (500 trees). These analyses were conducted using the randomForest package [86], and the significance of the model and predictor importance were determined using the rfUtilities and rfPermute packages in the R software, respectively [87, 88].

We then used structural equation modelling to evaluate the direct and indirect relationships between root ethylene and soil properties and plant physiology and bacterial community. We use the χ2 test (χ2; the model has a good fit when χ2 ≤ 2) and P values (traditionally ≥ 0.05), goodness-of-fit index and root mean square error of approximation (RMSEA, ≤ 0.05) and P values (traditionally ≥ 0.05) to evaluate the structural equation model fit [89]. The best-fitting and most parsimonious model was obtained after excluding all non-significant parameters. All SEM analyses were conducted using the Amos 17.0 software package (Smallwaters, IL, USA).

Statistical and network analyses

The differences in the peanut traits and the abundance of bacterial phyla in mono- and intercropping systems (two groups) in the field were analysed by one-way ANOVA. Significant differences (P < 0.05) in the soil chemical properties in the field and ethylene addition experiments, soil cyanide and plant hormone concentrations, the relative abundances of soil specific phyla, and ACS/ACO transcript expression between treatments in the pot or hydroponic experiments were evaluated by Tukey’s honest significant difference (HSD) test with SPSS 18.0 (SPSS, Chicago, IL, USA).

In this study, we used Shannon and Chao 1 indices to characterize bacterial ɑ-diversity [90], and the data were subjected to ANOVA using Tukey's HSD test at P < 0.05. To assess the influence of the different experimental factors on β-diversity, we calculated Bray-Curtis distances and then performed a canonical analysis of principal coordinates (CAP) constrained by the factor of interest and conditioned by the remaining variables [63, 91]. We employed the “capscale” and “permutest” permutation-based testing functions for the CAP analysis and the calculation of the significance values, respectively [63]. Permutational multivariate analysis of variance (PERMANOVA) was conducted to separate and quantitatively evaluate the effects of driving factors on the composition of the soil bacterial community using the “Anosims” function [92, 93]. Bootstrapped trees were constructed using the weighted pair group method with arithmetic mean (UPGMA) based on 1000 hierarchical clusters in the QIIME package [77]. Indicator species that were specifically associated with the different sampling sites were determined using the “labdsv” package in R (version 3.2.1) [67, 94].

Availability of data and materials

The sequences of 16S rRNA genes were deposited in the European Nucleotide Archive of EMBL under the accession number PRJEB22658.


  1. 1.

    Berendsen RI, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Lau JA, Lennon JT. Rapid responses of soil microorganisms improves plant fitness in novel environments. P Natl Acad Sci USA. 2012;109:14058–62.

    CAS  Article  Google Scholar 

  3. 3.

    Chagas FO, Pessotti RD, Caraballo-Rodriguez AM, Pupo MT. Chemical signaling involved in plant-microbe interactions. Chem Soc Rev. 2018;47:1652–704.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Huang AC, Jiang T, Liu YX, Bai YC, Reed J, et al. A specialized metabolic network selectively modulates Arabidopsis root microbiota. Science. 2019;6440:546.

    Google Scholar 

  5. 5.

    Bais HP, Weir TL, Perry LG, Gilroy S, Vivanco JM. The role of root exudates in rhizosphere interactions with plants and other organism. Annu Rev Plant Biol. 2006;57:233–66.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Lanoue A, Burlat V, Henkes GJ, Koch I, Schurr U, Rose USR. De novo biosynthesis of defense root exudates in response to Fusarium attack in barley. New Phytol. 2010;185:577–88.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Kostenko O, Mulder PPJ, Courbois M, Bezemer TM. Effects of plant diversity on the concentration of secondary plant metabolites and the density of arthropods on focal plants in the field. J Ecol. 2017;105:647–60.

    Article  Google Scholar 

  8. 8.

    Dudenhoffer JH, Scheu S, Jousset A. Systemic enrichment of antifungal traits in the rhizosphere microbiome after pathogen attack. J Ecol. 2016;104:566–1575.

    Article  CAS  Google Scholar 

  9. 9.

    Kong CH, Zhang SZ, Li YH, Xia ZC, Yang XF, et al. Plant neighbor detection and allelochemical response are driven by root-secreted signaling chemicals. Nat Commun. 2018;9:3867.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. 10.

    Lebeis SL, Paredes SH, Lundberg DS, Breakfield N, Gehring J, et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science. 2015;349:860–4.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Zhalnina K, Louie KB, Hao Z, Mansoori N, da Rocha UN, Shi SJ, et al. Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly. Nat Microbiol. 2018;3:470–80.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Pieterse CMJ, Leon-Reyes A, Van der Ent S, Van Wees SCM. 2009. Networking by small-molecule hormones in plant immunity. Nat Chem Biol. 2009;5:308–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Ravanbakhsh M, Sashmi S, Voesenek LA, Kowalchuk G, Jousset A. Microbial modulation of plant ethylene signaling: ecological and evolutionary consequences. Microbiome. 2018;6:52.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Cahill JF, McNickle GG, Haag JJ, Lamb EG, Nyanumba SM, Clair CCS. Plants integrate information about nutrients and neighbors. Science. 2010;328:1657.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Vandenkoornhuyse P, Quaiser A, Duhamel M, Le Van A, Dufresne A. The importance of the microbiome of the plant holobiont. New Phytol. 2015;206:1196–206.

    PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Carvalhais LC, Schenk PM, Dennis PG. Jasmonic acid signalling and the plant holobiont. Curr Opin Microbiol. 2017;37:42–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Broekgaarden C, Caarls L, Vos IA, Pieterse CMJ, Van Wees SCM. Ethylene: traffic controller on hormonal crossroads to defense. Plant Physiol. 2015;169:2371–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Oldroyd GED, Engstrom EM, Long SR. Ethylene inhibits the Nod factor signal transduction pathway of Medicago truncatula. Plant Cell. 2001;13:1835–49.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Barry CS, Fox EA, Yen HC, Lee S, Ying TJ, Grierson D, et al. Analysis of the ethylene response in the epinastic mutant of tomato. Plant Physiol. 2001;127:58–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Penmetsa RV, Uribe P, Anderson J, Lichtenzveig J, Gish JC, Nam YW, et al. 2008. The Medicago truncatula ortholog of Arabidopsis EIN2, sickle, is a negative regulator of symbiotic and pathogenic microbial associations. Plant J. 2008;55:580–95.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Riedel T, Groten K, Baldwin IT. Symbiosis between Nicotiana attenuata and Glomus intraradices: ethylene plays a role, jasmonic acid does not. Plant Cell Environ. 2008;31:1203–13.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Camehl I, Sherameti I, Venus Y, Bethke G, Varma A, Lee J, et al. Ethylene signalling and ethylene-targeted transcription factors are required to balance beneficial and nonbeneficial traits in the symbiosis between the endophytic fungus Piriformospora indica and Arabidopsis thaliana. New Phytologist. 2010; 185: 1062-73 (2010).

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Schulz-Bohm K, Gerards S, Hundscheid M, Melenhorst J, de Boer W, Garbeva P. Calling from distance: attraction of soil bacteria by plant root volatiles. ISME J. 2018;12:1252–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Conn EE. Cyanogenic glycosides. In: Stumpf PK, Conn EE, editors. The biochemistry of plants: a comprehensive treatise-secondary plant products. New York: Academic; 1981. p. 479–500.

    Google Scholar 

  25. 25.

    Jones DA. Why are so many food plants cyanogenic? Phytochemistry. 1998;47:155–62.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Gniazdowska A, Krasuska U, Bogatek R. Dormancy removal in apple embryos by nitric oxide or cyanide involves modifications in ethylene biosynthetic pathway. Planta. 2010;232:1397–407.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Xu F, Zhang DW, Zhu F, Tang H, Lv X, Cheng J, et al. A novel role for cyanide in the control of cucumber (Cucumis sativus L.) seedlings response to environmental stress. Plant Cell Environ. 2012;35:1983–97.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Li L, Tilman D, Lambers H, Zhang FS. Plant diversity and overyielding: insights from belowground facilitation of intercropping in agriculture. New Phytol. 2014;203:63–9.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  29. 29.

    Brooker RW, Bennett AE, Cong WF, Daniell TJ, George TS, Hallett PD, et al. Improving intercropping: a synthesis of research in agronomy, plant physiology and ecology. New Phytol. 2015;206:107–17.

    PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Li B, Li YY, Wu HM, Zhang FF, Li CJ, Li XX, et al. Root exudates drive interspecific facilitation by enhancing nodulation and N2 fixation. P Natl Acad Sci USA. 2016;113:6496–501.

    CAS  Article  Google Scholar 

  31. 31.

    Mason SC, Leihner DE, Vorst JJ. Cassava-cowpea and cassava-peanut intercropping. I. Yield and land use efficiency. Agron J. 1986;78:43–6.

    Article  Google Scholar 

  32. 32.

    Pierik R, Mommer L, Voesenek LACJ. Molecular mechanisms of plant competition: neighbor detection and response strategies. Funct Ecol. 2013;27:841–53.

    Article  Google Scholar 

  33. 33.

    Abakumova M, Zobel K, Lepik A, Semchenko M. Plasticity in plant functional traits is shaped by variability in neighbourhood species composition. New Phytol. 2016;211:455–63.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Sletvold N. Effect of plant size on reproductive output and offspring performance in the facultative biennial Digitalis purpurea. J Ecol. 2002;6:958–66.

    Article  Google Scholar 

  35. 35.

    Anten NPR, Vermeulen PJ. Tragedies and crops: understanding natural selection to improve cropping systems. Trends Ecol Evol. 2016;31:429–39.

    PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    McNickle GG, Dybzinski R. Game theory and plant ecology. Ecol Lett. 2013;16:545–55.

    PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Philippot L, Raaijmakers JM, Lemanceau P, Putten WH. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol. 2013;11:789–99.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Forlani S, Masiero S, Mizzotti C. Fruit ripening: the role of hormones, cell wall modifications, and their relationship with pathogens. J Exp Bot. 2019;11:2993–3006.

    Article  CAS  Google Scholar 

  39. 39.

    Bennett EJ, Roberts JA, Wagstaff C. The role of the pod in seed development: strategies for manipulating yield. New Phytol. 2011;4:838–53.

    Article  Google Scholar 

  40. 40.

    Gallo M, Sayre R. Removing allergens and reducing toxins from food crops. Curr Opin Biotech. 2009;20:191–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Siegieñ I, Bogatek R. Cyanide action in plants - from toxic to regulatory. ACTA Physiol Plant. 2006;28:483–97.

    Article  Google Scholar 

  42. 42.

    Yang Z, Tian LN, Latoszek-Green M, Brown D, Wu KQ. Arabidopsis ERF4 is a transcriptional repressor capable of modulating ethylene and abscisic acid responses. Plant Mol Biol. 2005;58:585–96.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Yang CY, Hsu FC, Li JP, Wang NN, Shih MC. The AP2/ERF transcription factor AtERF73/HRE1 modulates ethylene responses during hypoxia in Arabidopsis. Plant physiol. 2011;156:202–12.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Acera F, Carmona MI, Castillo F, Quesada A, Blasco R. A cyanide-induced 3- cyanoalanine nitrilase in the cyanide-assimilating bacterium Pseudomonas pseudoalcaligenes strain CECT 5344. App Environ Microb. 2017;83:e00089–17.

    CAS  Article  Google Scholar 

  45. 45.

    Semchenko M, Saar S, Lepik A. Plant root exudates mediate neighbour recognition and trigger complex behavioural changes. New Phytol. 2014;204:631–7.

    PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Clark DG, Gubrium EK, Barrett JE, Nell TA, Klee HJ. Root formation in ethylene-insensitive plants. Plant Physiol. 1999;121:53–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Hahn A, Zimmermann R, Wanke D, Harter K, Edelmann HG. The root cap determines ethylene-dependent growth and development in maize roots. Mol Plant. 2008;1:359–67.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Song QX, Ando A, Xu DQ, Fang L, Zhang TZ, Huq E, et al. Diurnal down-regulation of ethylene biosynthesis mediates biomass heterosis. P Natl Acad Sci USA. 2018;115:5606–11.

    CAS  Article  Google Scholar 

  49. 49.

    Pretty J, Bharucha ZP. Sustainable intensification in agricultural systems. Ann Bot. 2014;8:1571–96.

    Article  Google Scholar 

  50. 50.

    Vavrova K, Knapek J, Weger J. Modeling of biomass potential from agricultural land for energy utilization using high resolution spatial data with regard to food security scenarios. Renew Sust Energ Rev. 2014;35:436–44.

    Article  Google Scholar 

  51. 51.

    Bakker PAHM, Pieterse CMJ, de Jonge R, Berendsen RL. The soil-borne legacy. Cell. 2018;6:1178–80.

    Article  CAS  Google Scholar 

  52. 52.

    Peiffer JA, Spor A, Koren O, Jin Z, Tringe SG, Dangl JL, et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. P Natl Acad Sci USA. 2013;6:1178–80.

    Google Scholar 

  53. 53.

    Bargaz A, Noyce GL, Fulthorpe R, Carlsson G, Furze JR, Jensen ES, et al. Species interactions enhance root allocation, microbial diversity and P acquisition in intercropped wheat and soybean under P deficiency. Appl Soil Ecol. 2017;120:179–88.

    Article  Google Scholar 

  54. 54.

    Wu X, Wu FZ, Zhou XG, Fu XP, Tao Y, Xu WH, et al. Effects of intercropping with potato onion on the growth of tomato and rhizosphere alkaline phosphatase genes diversity. Front Plant Sci. 2016;7:846.

    PubMed  PubMed Central  Google Scholar 

  55. 55.

    Fu YM, Shao LZ, Liu H, Tong L, Liu H. Ethylene removal evaluation and bacterial community analysis of vermicompost as biofilter material. J Hazard Mater. 2011;192:658–66.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  56. 56.

    Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2012;8:538–50.

    Article  CAS  Google Scholar 

  57. 57.

    Banerjee S, Schlaeppi K, van der Heijden MGA. Keystone taxa as drivers of microbiome structure and functioning. Nat Rev Microbiol. 2018;16:567–76.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  58. 58.

    Shi SJ, Nuccio EE, Shi ZJ, He ZL, Zhou JZ, Firestone MK. The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages. Ecol Lett. 2016;19:926–36.

    PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Deng Y, Jiang YH, Yang YF, He ZL, Luo F, Zhou JZ. Molecular ecological network analyses. BMC Bioinformatics. 2012;13:113.

    PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Newman MEJ. Modularity and community structure in networks. P Natl Acad Sci USA. 2006;103:8577–82.

    CAS  Article  Google Scholar 

  61. 61.

    Shetty SA, Hugenholtz F, Lahti L, Smidt H, de Vos WM. Intestinal microbiome landscaping: insight in community assemblage and implications for microbial modulation strategies. FEMS Microbiol Rev. 2017;41:182–99.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Guimera R, Amaral L. Functional cartography of complex metabolic networks. Nature. 2005;459:193–9.

    Google Scholar 

  63. 63.

    Jiang YJ, Liu MQ, Zhang JB, Chen Y, Chen XY, Chen LJ, et al. Nematode grazing promotes bacterial community dynamics in soil at the aggregate level. ISME J. 2017;11:2705–17.

    PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Hazen EL, Brown R. Antifungal agent produced by a soil Actinomycete. Science. 1950;112:423.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Gordon JJ, Miller GA, Kelly BK. Actinonin-an antibiotic substance produced by an Actinomycete. Nature. 1962;195:701.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  66. 66.

    Zettler J, Xia HY, Burkard N, Kulik A, Grond S, Heide L, et al. New aminocoumarins from the rare Actinomycete Catenulispora acidiphila DSM 44928: identification, structure elucidation, and heterologous production. Chembiochem. 2014;15:612–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  67. 67.

    Dufrene M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monog. 1997;67:345–66.

    Google Scholar 

  68. 68.

    Qin S, Zhang YJ, Yuan B, Xu PY, Xing K, Wang J, et al. Isolation of ACC deaminase-producing habitat-adapted symbiotic bacteria associated with halophyte Limonium sinense (Girard) Kuntze and evaluating their plant growth-promoting activity under salt stress. Plant Soil. 2014;374:753–66.

    CAS  Article  Google Scholar 

  69. 69.

    Li ZY, Chang SP, Ye ST, Chen MY, Lin L, Li YY, et al. Differentiation of 1-aminocyclopropane-1-carboxylate (ACC) deaminase from its homologs is the key for identifying bacteria containing ACC deaminase. FEMS Microbiol Ecol. 2015; 91: fiv112.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  70. 70.

    Delgado-Baquerizo M, Maestre FT, Reich PB, Jeffries TC, Gaitan JJ, Encinar D, et al. Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat Comms. 2015;7:10541.

    Article  CAS  Google Scholar 

  71. 71.

    Shao Y, Bao W, Chen D, Eisenhauer N, Zhang W, Pang X, et al. Using structural equation modeling to test established theory and develop novel hypotheses for the structuring forces in soil food webs. Pedobiologia. 2015;58:137–45.

    Article  Google Scholar 

  72. 72.

    United States Environmental Protection Agency (EPA). Titrimetric and manual spectrophotometric determinative methods for cyanide. Method, 1996; p.9014.

  73. 73.

    Wu JJ, Wang CA, Zheng LQ, Wang L, Chen YL, Whelan J, et al. Ethylene is involved in the regulation of iron homeostasis by regulating the expression of iron-acquisition-related genes in Oryza sativa. J Exp Bot. 2011;62:667–74.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  74. 74.

    Zhu XF, Zhu CQ, Zhao XS, Zheng SJ. Shen RF. Ethylene is involved in root phosphorus remobilization in rice (Oryza sativa) by regulating cell-wall pectin and enhancing phosphate translocation to shoots. Ann Bot. 2016; 118: 645-653.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. 75.

    Clevenger J, Chu Y, Scheffler B, Ozias-Akins P. A developmental transcriptome map for allotetraploid Arachis hypogaea. Front Plant Sci. 2016;7:1446.

    PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2013;8:790–803.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  77. 77.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. 78.

    Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  79. 79.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

    CAS  Article  Google Scholar 

  80. 80.

    Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Microbiol Biotech. 2007;73:5261–7.

    CAS  Google Scholar 

  81. 81.

    Maslov S, Sneppen K. Specificity and stability in topology of protein networks. Science. 2002;296:910–3.

    CAS  Article  Google Scholar 

  82. 82.

    Zhou JZ, Deng Y, Luo F, He ZL, Tu QC. Zhi XY. Functional molecular ecological networks. mBio. 2010;1:e00169–10.

    PubMed  PubMed Central  Google Scholar 

  83. 83.

    Zhou JZ, Deng Y, Luo F, He ZL, Yang YF. Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. mBio. 2011; 2: e00122-11.

  84. 84.

    Olesen J, Bascompte J, Dupont Y, Jordano P. The modularity of pollination networks. P Natl Acad Sci USA. 2007;104:19891–6.

    CAS  Article  Google Scholar 

  85. 85.

    Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. San Jose, California: In International AAAI Conference on weblogs and social media; 2009.

    Google Scholar 

  86. 86.

    Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2:18–22.

    Google Scholar 

  87. 87.

    Archer E. rfPermute: estimate permutation p-values for random forest importance metrics. .

  88. 88.

    Evans JS, Murphy MA. rfUtilities: random forests model selection and performance evaluation. Retrieved from https://CRAN. = rfUtilities. .

  89. 89.

    Hooper D, Coughlan J, Mullen M. Structural equation modeling: guidelines for determining model fit. Electron J Bus Res Methods. 2008;6:53–60.

    Google Scholar 

  90. 90.

    Gotelli NJ, Colwell RK. Estimating species richness. In: Magurran AE, McGill BJ, editors. Biological diversity: frontiers in measurement and assessment. UK: Oxford University; 2011. p. 39–54.

    Google Scholar 

  91. 91.

    Anderson MJ, Willis TJ. Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology. 2003;84:511–25.

    Article  Google Scholar 

  92. 92.

    Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.

    Google Scholar 

  93. 93.

    Warton DI, Wright ST, Wang Y. Distance based multivariate analyses confound location dispersion effects. Methods Ecol Evol. 2012;3:89–101.

    Article  Google Scholar 

  94. 94.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references


We thank Renfeng Tu (Anhui Agricultural University), Jianbo Fan and Zhiying Yu (Red Soil Ecological Experimental Station of Institute of Soil Science) for assistance in the management of the field experiment and soil sampling and Ming Zhong and Hong Chen (Soil and Environmental Analysis Center of Institute of Soil Science) for assistance with the GC-MS analysis. We thank Wenzhi Zhou (Shanghai Sanshu Biotechnology Co., LTD) for assistant with the phytohormone detection.


This work was supported by the National Key R&D Program (2016YFD0200309), National Natural Science Foundation of China (NSFC No. 41977098), Distinguished Youth Scholar Program of Jiangsu Province (BK20180049), Modern Agricultural Industry Technology System of the Ministry of Agriculture of China (CARS-22-G-14) and the Natural Science Foundations of Jiangsu Province (BK20191510 and BK20171520).

Author information




Y.C. and B.S. designed the research. Y.C., Y.S, Y.J. and X.W. performed experiments and conducted fieldwork. Y.C., M.B., B.S.G. and B.S. analysed the data. Y.C. M.B. and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bo Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Methods S1.

Soil chemical property and enzymatic activity determination, and peanut phytohormone measurements. Fig. S1. Results of the pot experiments conducted for the detection of peanut phytohormone production. Fig. S2. Hydroponic culturing to detect the regulation of root ethylene production in peanut by exogenous cyanide. Fig. S3. Hierarchical clustering of soil bacterial community in cassava and peanut intercropping (PC) and peanut monocropping (PP) systems based on pairwise Bray-Curtis distances. Fig. S4. Zi-Pi plot showing the distribution of OTUs based on their topological roles in the intercropping and monocropping networks. Fig.S5. Soil enzymatic activity of organic phosphorus and nitrogen mineralization. Table S1. Soil chemical properties in the rhizosphere and bulk soil in the intercropping and monocropping systems. Table S2. Information regarding the indicators in ethylene incubating samples. Table S3. Gene-specific primers used for qRT-PCR.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, Y., Bonkowski, M., Shen, Y. et al. Root ethylene mediates rhizosphere microbial community reconstruction when chemically detecting cyanide produced by neighbouring plants. Microbiome 8, 4 (2020).

Download citation


  • Ethylene signal
  • Neighbour detection
  • Chemical cue
  • Interspecific interaction
  • Cyanide
  • Rhizosphere microbial assemblage
  • Plant fitness