Puchang Wang
Leilei Ding
Chao Zou3- 1School of Life Sciences, Guizhou Normal University, Guiyang, China
- 2Guizhou Institute of Prataculture, Guizhou Academy of Agricultural Sciences, Guiyang, China
- 3College of Animal Science, Guizhou University, Guiyang, China
Microbes, especially abundant microbes in bulk soils, form multiple ecosystem functions, which is relatively well studied. However, the role of rhizosphere microbes, especially rhizosphere rare taxa vs. rhizosphere abundant taxa in regulating the element circling, multifunctionality, aboveground net primary productivity (ANPP) and the trade-offs of multiple functions remains largely unknown. Here, we compared the multiple ecosystem functions, the structure and function of rhizosphere soil bacterial and fungal subcommunities (locally rare, locally abundant, regionally rare, regionally abundant, and entire), and the role of subcommunities in the Zea mays and Sophora davidii sole and Z. mays/S. davidii intercropping ecosystems in subtropical China. Results showed that intercropping altered multiple ecosystem functions individually and simultaneously. Intercropped Z. mays significantly decreased the trade-off intensity compared to sole Z. mays, the trade-off intensity under intercropped S. davidii was significantly higher than under intercropped Z. mays. The beta diversities of bacterial and fungal communities, and fungal functions in each subcommunity significantly differed among groups. Network analysis showed intercropping increased the complexity and positive links of rare bacteria in Z. mays rhizosphere, but decreased the complexity and positive links of rare bacteria in S. davidii rhizosphere and the complexity and positive links of fungi in both intercropped plants rhizosphere. Mantel test showed significant changes in species of locally rare bacteria were most strongly related to nitrogen-cycling multifunctionality, ANPP and trade-offs intensity, significant changes in species of locally rare fungus were most strongly related to carbon-cycling multifunctionality, phosphorus-cycling multifunctionality, and average ecosystem multifunctionality. This research highlights the potential and role of rare rhizosphere microorganisms in predicting and regulating system functions, productivity, and trade-offs.
Introduction
Losses in taxonomic and functional diversities are pervasively at both global and local scales (Radchuk et al., 2016; Baldrighi et al., 2017; Huang et al., 2019). This trend is predicted to continue over this century (Huang et al., 2019), and raises increasing concerns on the influence of biodiversity on ecosystem functions (Karolína et al., 2014; Radchuk et al., 2016; Baldrighi et al., 2017). Unlike flora (Niklaus et al., 2017; Wang et al., 2020b) or fauna (Gagic et al., 2015; Tonin et al., 2018) both which have been well studied in the studies of the relationship between biodiversity and ecosystem function, soil microorganisms represent the richest and highest diverse life (Delgado-Baquerizo et al., 2017; Luo et al., 2018; Chen et al., 2020a), nevertheless, the relationships between soil microorganisms and ecosystem functions are not fully understood (Delgado-Baquerizo et al., 2017). Previous studies have found that soil microbes contribute to driving multiple ecosystem functions simultaneously (MF) (Fanin et al., 2017; Jiao et al., 2018; Chen et al., 2020a), such as carbon and nutrient cycling (Delgado-Baquerizo et al., 2017; Jiao et al., 2018) and productivity (George et al., 2019; Zheng et al., 2019). However, most studies have centered on the temperate communities (Guerra et al., 2020) and the bulk soil communities (Luo et al., 2018; Wen et al., 2020). In contrast, subtropical data on soil microbial diversity and ecosystem functions are particularly scant (Guerra et al., 2020). Furthermore, the rhizosphere is a true hot point of plant-microbial-soil interactions, and the rhizosphere microbial communities are clearly distinctive from the surrounding bulk soil’s (Fan et al., 2017; Jiao et al., 2018), because plants filter microbes for special structures and function (Del Galdo Jos et al., 2003), to benefit their growth, nutrition (Mendes et al., 2014), and function (Lu et al., 2018). Notwithstanding, how rhizosphere bacterial and fungal communities participate in the nutrient cycling of rhizosphere soil, and drive the rhizosphere multiple functions, aboveground primary productivity and trade-offs among functions remains less explicitly acknowledged. Although bacteria and fungi are the most frequently studied communities in soil biodiversity and ecosystem function research (Guerra et al., 2020), their diversity effects on ecosystem functions have not yet been fully explored for the rhizosphere soils.
Besides, the majority of studies were based on abundant taxa of microbiome (Chen et al., 2020a; Liang et al., 2020). A large number of low abundance taxa (rare microbial taxa) were deleted before analyzing data of microbiome (Lynch and Neufeld, 2015; Jousset et al., 2017; Chen et al., 2020a). However, there is a theoretical contradiction for this deletion. On the one hand, the influential mass ratio hypothesis states (Grime, 1998) that the influences of species/functions on an ecosystem function/process is in proportion to their biomass/relative abundance (Karolína et al., 2014; Bagousse-Pingueta et al., 2019), in this sense, it may be reasonable to focus only on the role of abundant taxa. On the other hand, the most of taxa in almost all ecosystems are low abundance (rare) (David et al., 2013; Jain et al., 2014; Lynch and Neufeld, 2015), and dominant species account for most of the total abundance (Dee et al., 2019), resulting in high diversity in the rare subcommunity (Jain et al., 2014; Lynch and Neufeld, 2015) and low diversity in the dominant subcommunity. Generally, the selection effect [selection/occurrence of some particular species (or functions)/identity effects] and complementary effect (niche partitioning/different resource utilization/facilitation) are recognized mechanisms to interpret the role of biodiversity in shaping ecosystem functions (Loreau, 2000; Mensah et al., 2020; Ding and Wang, 2021). Since high species diversity should increase the likelihood of the selection effect (Mensah et al., 2020) and/or the complementary effect (Mensah et al., 2018), rare subcommunities or their functions is inferred to contribute in greater proportion to a given ecosystem function than the dominant. This contradicts the mass ratio hypothesis. Besides, the deletion changed substantially the profile of rare taxa (Chen et al., 2020a). Unfortunately, rare species are often more sensitive (Jain et al., 2014; Guo et al., 2020; Zhou et al., 2020b) and vulnerable to vanish firstly (David et al., 2013). Therefore, if the above inference is true, conservation and use of rare species will be more imperative than those of the abundance species to maintain ecosystem functions. However, the ecological role of rare species is poorly known (Säterberg et al., 2019). Although the complementarity and selection effects contribute to MF, their relative contributions remains controversial (Mensah et al., 2020). Moreover, the relative importance of these two effects of rare rhizosphere microbes in explaining rhizosphere functions, plant productivity and trade-offs are understudied.
In this study, our hypothesis was that rare rhizosphere taxa contributed to multiple functions and trade-offs in larger proportion than the abundant taxa did. To examine this hypothesis, we characterized the rhizosphere abundant and rare bacteria and fungi, rhizosphere C/N/P-cycling multifunctionality, ecosystem multifunctionality, aboveground net primary productivity and trade-offs, as well as explored how bacteria and fungi with different abundance differently link to the multiple rhizosphere functions, aboveground net primary productivity and trade-offs in subtropical sole and intercropping systems in China.
Materials and methods
Study area and design
Sampling was conducted in a cropping common garden in the Guizhou Academy of Agricultural Sciences (26°30′N, 106°39′E, 1,100 m a.s.l.). This region undergoes a north subtropical monsoon climate characterized by annual precipitation of 1,130 mm, annual average air temperature of 15.3°C, and Haplic alisol soil. The sole and intercropping systems were selected based on the fact that monoculture systems can reduce biodiversity (Shao et al., 2020), intercropping systems can enhance biodiversity (Martin-Guay et al., 2018; Dingha et al., 2021), and both systems are widely used around the world. The experimental design was a randomized block design (Luo et al., 2018) established using three cropping systems, including Sophora davidii (Franch.) monocropping (MS), Zea mays L. monocropping (MZ), and intercropping of both (IS and IZ), in May 1st 2019. Five repeated plots (2 m × 5 m per repeat) were arranged for each system with intervals of 1 m. MS was planted with row distances and plant distances of 0.6 m, MZ was planted with row distances of 0.6 m and plant distance of 0.3 m. Row and plant distances of S. davidii (Franch.) and Z. mays L. in the intercropping system were identical to in monocropping system and the distance between adjacent rows was 0.3 m in the intercropping system. Since this small-scale design can minimize the variation in either of none-design factors, such as topographies, climate, and soil type (Chen et al., 2020c; Ding and Wang, 2021), this approach is convenient to compare the contributions of diversity vs. other drivers in shaping ecosystem function (van der Plas, 2019), allowing the resolution required to capture relationships between the plant-selected microbes and the rhizosphere element circling, aboveground productivity and trade-offs.
Fresh aboveground S. davidii and Z. mays were harvested using a 1 m2 (1 × 1) quadrat in the monoculture and intercropping plots, respectively, in August 11th 2019. After the fresh was killed at 105°C and oven-dried at 65°C, the ANPP was obtained. Soil that tightly adhering to roots (rhizosphere soils) was collected by trembling the roots (Xiong et al., 2020), and was mixed to obtain a composite sample for each replicate. The soil was sieved using 2-mm meshes. A total of 20 composite samples [(two monocropping + two intercropping) × five replicates] were obtained. The sample was divided into subsamples for downstream analysis.
Individual rhizosphere functions
The rhizosphere physicochemical properties and C, N, and P-circling enzyme activities were assayed using the methods listed in previous studies (Zheng et al., 2019; Ding et al., 2020a,b) and the Supplementary Description of Method. Briefly, rhizosphere soil pH was measured with a suspension (soil: water = 1:2.5 w/v); Water content (WC, %) was obtained by oven-drying at 105°C (Wang et al., 2022). Rhizosphere C-circling functions were characterized as follows (Bowker et al., 2013; Bastida et al., 2016; Luo et al., 2018; Ding et al., 2020a): organic carbon (g kg−1, OC) was assayed using potassium dichromate volumetric method, β-glucosidase (C-circling enzyme, μmol d−1 g−1 dry soil, βG) was detected using an ELISA test kit (Shanghai Enzyme-linked Biotechnology Co., Ltd., China). Microbial biomass carbon (mg kg−1, MBC) was determined by the fumigation–extraction method. Rhizosphere N-circling functions were characterized as follows (Jiao et al., 2018; Luo et al., 2018; Bagousse-Pingueta et al., 2019; Ding et al., 2020a): ammoniacal nitrogen (mg kg−1, NH4_N) was assayed using ultraviolet spectrophotometry, nitrate nitrogen (mg kg−1, NO3_N) was using the colorimetric method based on extraction of potassium chloride, inorganic nitrogen (mg kg−1, IN) was the sum of ammoniacal and nitrate nitrogen, five N-circling enzymes [N-acetylglucosaminidase (NAG); Leucine aminopeptidase (LAP), Nitrogenase, Nitric oxide synthetase (NOS), Glutamine synthetase (GS); μmol d−1 g−1 dry soil] were detected using ELISA kits. Microbial biomass nitrogen (mg kg−1, MBN) was determined by the fumigation–extraction method. Rhizosphere P-circling functions were characterized as follows (Bowker et al., 2013; Bastida et al., 2016; Jiao et al., 2018; Luo et al., 2018; Bagousse-Pingueta et al., 2019; Ding et al., 2020a): total phosphorus (mg kg−1, TP) was determined by NaOH digestion, available phosphorus (mg kg−1, AP) was assayed via the NaHCO3-ultraviolet spectrometer, and acid phosphatase (P-circling enzyme, μmol d−1 g−1 dry soil, ACP) was detected using ELISA test kits. These properties were used because they either measure real functions or are good surrogates of C, N, and P-cycling (Bowker et al., 2013; Delgado-Baquerizo et al., 2017; Bagousse-Pingueta et al., 2019; Chen et al., 2020a,b).
DNA extraction and PCR amplification
Genomic DNA was extracted using HiPure Soil DNA Kits (Guangzhou Meiji Biotechnology Co., Ltd., China) following the manufacturer’s instructions. The quality of extracted DNA was checked with electrophoresis in a 1.2% (wt/vol) agarose gel and a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, United States). The V3–V4 region of the bacterial 16S rRNA gene and the internal transcribed spacer (ITS) regions of the fungal rRNA gene were amplified with primers 341F and 806R and ITS3_KYO2 and ITS4, respectively (detailed PCR conditions are described in Supplementary Description of Method).
Amplicon sequencing and data processing
Amplicons were extracted from 2% agarose gels (Gong et al., 2021) and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) according to the manufacturer’s instructions and quantified using ABI StepOnePlus Real-Time PCR System (Life Technologies, Foster, CA, United States). Purified amplicons were pooled in equimolar and paired-end sequenced (PE250) on an Illumina platform by Guangzhou Genedenovo Biotechnology Co., Ltd. (Guangzhou, China). The raw data were deposited in NCBI (Accession number: PRJNA731978 for 16S, PRJNA731989 for ITS).
Raw data were processed as follows: briefly, reads containing more than 10% of unknown nucleotides and containing less than 50% of bases with quality (Q-value) > 20 were removed using FASTP (V 0.18.01). Paired-end clean reads were merged as raw tags using FLSAH (V 1.2.112) with a minimum overlap of 10 bp and mismatch error rates of 2%. The unique tag sequence was selected using Mothur (v1.39.13). Chimeras were removed using UCHIME algorithm in USEARCH4. Sequences were then assigned to operational taxonomic unit (OTUs) at a 97% similarity threshold (Ding et al., 2020a; Ding and Wang, 2021; Jin et al., 2022b) using UPARSE (USEARCH v9.2.64, see Text Footnote 4). The representative OTU sequences were classified into organisms by a naive Bayesian model using RDP classifier (version 2.2) based on SILVA database (version 1325) for bacteria or UNITE database (version 8.06) for fungi. Bioinformatics analysis was done using an online platform7 by Guangzhou Genedenovo Biotechnology Co., Ltd. (Guangzhou, China).
Definition of abundant and rare microbial taxa
The locally abundant and rare OTUs were defined using abundance thresholds of 0.1 and 0.01% in a sample (Lynch and Neufeld, 2015), and regionally abundant and rare OTUs were defined by abundance thresholds of 0.01 and 0.001% across samples (Mo et al., 2018). The subsequent analyses were done at five levels: whole, locally abundant, locally rare, regionally abundant, and regionally rare OTUs. The richness, Simpson diversity, and Shannon diversity indices were applied to describe the taxonomic diversity.
Predicted functions of bacteria and fungi
Tax4Fun and FAPROTAX 1.2.3 were used to predict bacterial functional profiles (Aßhauer et al., 2015; Louca et al., 2016). FUNGuild (Guilds_v1.1) was used to ppredicting fungal functional profiles (Nguyen et al., 2016). The ‘‘Vegan’’ package was applied to calculate the richness, Shannon diversity, and Simpson diversity of bacterial and fungal functions at the commonly used level (i.e., Tax4Fun KEGG 2 level and FUNGuild Guild level) in R v3.5.38.
Multifunctionality and trade-off intensity
Multifunctionality is a crucial management and ecological index (Delgado-Baquerizo et al., 2020), and is defined as the synthesis of different ecosystem properties at a small scale (Bagousse-Pingueta et al., 2019; Ding and Wang, 2021). Rhizosphere C/N/P-circling multifunctionality (CCMF, NCMF, and PCMF) was calculated based on rhizosphere C/N/P-circling functions (Garibotti et al., 2018) in the Section “Individual rhizosphere functions,” and rhizosphere multifunctionality was calculated based on all rhizosphere functions (Bagousse-Pingueta et al., 2019) listed in Section “Individual rhizosphere functions.” Z-score transformation was performed to standardize the data of microbial community and ecosystem functions (Bastida et al., 2016; Bagousse-Pingueta et al., 2019) using the “standardizeZScore” function (Byrnes et al., 2014) in R before calculation of multifunctionality (Chen et al., 2020a). Average approach (the mean of all standardized functions) is intuitive and easily interpretable and widely used in the multifunctionality studies (Delgado-Baquerizo et al., 2016, 2017), therefore, average approaches-based multifunctionality index for each sample were calculated (Shi et al., 2021). Trade-off intensity was defined as the geometric distance from a point to a straight line of 1:1 to quantify the trade-off intensity between two properties of ecosystem (Zhong et al., 2020).
Data analysis
The statistical analysis followed the work flow and methods in previous studies (Ding et al., 2020a; Ding and Wang, 2021) and in the Supplementary Description of Method. The Shapiro–Wilk normality test and Levene’s test were used to test the normality and homoscedasticity of the data, respectively (Gao et al., 2021). When the data could meet the normality and homoscedasticity criteria simultaneously, ANOVA and t-test were used to test the significance of difference among groups. When the data could not meet, Kruskal–Wallis rank sum test and Wilcoxon test were used. The ηp2 (partial eta-squared) statistic was conducted to test the relative influences of plant identity (S. davidii and Z. mays), system type (sole and intercropping) and interaction on each individual rhizosphere functions in IBM SPSS (version 25, IBM, Armonk, NY, United States). Principal component analysis (PCA) (Merino-Martín et al., 2021) was used to determine which functions are the main ecosystem functions in the rhizosphere. The ANOSIM (Analysis of Similarities) test with 9,999 permutations was used to determine significant differences in rhizosphere ecosystem functions between plant species (S. davidii and Z. mays), system type (sole and intercropping) and among groups (MZ, IZ, MS, and IS) and in species/functions among groups for each subcommunity. Correlation relationships among individual rhizosphere ecosystem functions were visualized based Spearman using “igraph” package (Csardi and Nepusz, 2006). To retain the number of edges, vertices, and positive edges of microbial association networks, robust correlations were built based absolute value of Pearson correlation’s r > 0.8 and false discovery rate-corrected p < 0.001 using ‘‘WGCNA’’ and ‘‘igraph’’ packages in R. Gephi 0.9.29 were applied to visualize the networks.
Non-metric multidimensional scaling (NMDS) (Jiao et al., 2017; Khashi u Rahman et al., 2021) was applied to simplify samples or OUTs in high-dimensional to low-dimensional space for location using the “metaMDS” function based on the Bray–Curtis distance (He et al., 2021). Linear discriminant analysis (LDA) effect size (LEfSe) was run to determine the significantly enriched clades (LDA scores ≥ 2 and p < 0.05).
According to the definition of the selection effect (Mensah et al., 2020), we used Spearman rank correlations (Zheng et al., 2019) to select the species, functions, and diversity of the whole community, and abundant and rare subcommunities that were significantly related to C/N/P-circling multifunctionality (CCMF, NCMF, and PCMF), average ecosystem multifunctionality (AEMF), and aboveground net primary productivity (ANPP), respectively. We used Kruskal–Wallis rank sum test to select the species, functions, and diversity of the whole, abundant, and rare subcommunities that significantly changed among groups. We then calculated the intersection of species/functions/diversities that significantly changed and were significantly related to the above ecosystem functions. To decipher how microbes link multifunctionality (Delgado-Baquerizo et al., 2017), the Mantel test (Ding et al., 2020a) with 9,999 permutations was performed to detect the correlation between the intersection and above ecosystem functions, based on the Spearman method and Euclidean distance matrix. The greater the correlation of the Mantel test, the greater the potential impact intensity of species/functions/diversities on the above ecosystem functions (Jiao et al., 2017; Mo et al., 2018; Zheng et al., 2019; Zhou et al., 2019b; Xiong et al., 2020; Xue et al., 2020; Ding and Wang, 2021). The “ggplot2” (Wickham, 2016) and “circlize” (Gu, 2014) packages were used to visualize the results.
Results
Intercropping altered individual functions, multifunctionality, aboveground net primary productivity, and trade-offs
Compared with S. davidii, regardless of system type, the remarkably higher rhizosphere OC of Z. mays was observed (71% higher under monoculture, t-test p < 0.05; 40% higher under intercropping, t-test p < 0.05, Figure 1A). Compared with monoculture, intercropped Z. mays did not significantly change OC, intercropped S. davidii significantly elevated OC by 33% (t-test p < 0.05). Regardless of system type, statistical difference was not observed in MBC between Z. mays and S. davidii (t-test p > 0.05), (Figure 1B). Compared with under sole S. davidii, MBC under intercropping has increasing trend (by 14%, t-test p > 0.05), but MBC under intercropped Z. mays significantly decreased by 17% (t-test p < 0.05). Compared with monoculture, intercropped Z. mays did not significantly change the βG activity in rhizosphere (t-test p > 0.05, Figure 1C), but intercropped S. davidii increased the βG activity by 20% (t-test p < 0.05), therefore, despite statistical difference was not found in the βG activity between sole Z. mays and S. davidii (t-test p > 0.05), the βG activity under intercropped S. davidii was higher 23% than under intercropped Z. mays (t-test p < 0.05). WC was higher under sole S. davidii than under sole Z. mays by 22% (t-test p < 0.05, Figure 1D). Compared with monoculture, intercropped Z. mays decreased WC by 14% (t-test p < 0.05), intercropped S. davidii increased WC by 9% (t-test p < 0.05), resulting in insignificant difference in WC between intercropped S. davidii and Z. mays (t-test p > 0.05). NH4_N was higher under sole Z. mays than under sole S. davidii by 189% (Wilcoxon test p < 0.05, Figure 1E), NO3_N was higher under sole S. davidii than under sole Z. mays by 53% (t-test p < 0.05, Figure 1F), statistical difference was not detected between two plants, and between sole and intercropping (t-test p > 0.05). MBN was lower under sole S. davidii than under sole Z. mays by 49% (t-test p < 0.05, Figure 1G), but NAG, Nitrogenase, NOS, and GS under sole S. davidii were higher than under sole Z. mays by 36, 26, 34, and 33%, respectively (t-test p < 0.05, Figures 1H,I,K,L), the LAP activity under sole S. davidii were higher than under sole Z. mays by 12% (t-test p > 0.05, Figure 1J). Compared with monoculture, intercropped Z. mays decreased MBN by 47% (t-test p > 0.05), but increased NAG, LAP, and NOS by 40, 18, and 23%, respectively (t-test p < 0.05). Intercropped S. davidii increased MBN by 65% (t-test p > 0.05). MBN, Nitrogenase, NOS, and GS under intercropped S. davidii were higher than under intercropped Z. mays by 59, 21, 16, and 30%, respectively (t-test p < 0.05). TP and AP under sole Z. mays were higher than under sole S. davidii by 50 and 163%, respectively (t-test p < 0.05, Figures 1M,N). Intercropped Z. mays did not significantly change TP, AP, and ACP, in comparison to sole Z. mays (t-test p < 0.05, Figures 1M–O), intercropped S. davidii increased TP and AP in comparison to sole S. davidii (t-test p < 0.05), and statistical difference was not found in TP between intercropped S. davidii and Z. mays (t-test p > 0.05). AP under intercropped Sophora davidii was lower than under intercropped Z. mays (t-test p < 0.05), ACP under intercropped S. davidii was higher than under intercropped Z. mays (t-test p < 0.05). Statistical differences were not observed in IN and pH among group (Supplementary Figure 1, t-test p > 0.05).
Figure 1. Intercropping altered multiple rhizosphere functions. “ns,” p > 0.05; “*,” p < 0.05; “**,” p < 0.01; “***,” p < 0.001; MZ, monoculture Zea mays; IZ, intercropping Zea mays; MS, monoculture Sophora davidii; IZ, intercropping Sophora davidii. (A) Intercropping altered the organic carbon (OC), (B) intercropping altered the microbial biomass carbon content (MBC), (C) intercropping altered the β-glucosidase activity (βG), (D) intercropping altered the water content (WC), (E) intercropping altered the ammoniacal nitrogen content (NH4_N), (F) intercropping altered the nitrate nitrogen content (NO3_N), (G) intercropping altered the microbial biomass nitrogen content (MBN), (H) intercropping altered the N-cetylglucosaminidase activity (NAG), (I) intercropping altered Nitrogenase activity, (J) intercropping altered the Leucine aminopeptidase activity (LAP), (K) intercropping altered the Nitric oxide synthetase activity (NOS), (L) intercropping altered the Glutamine synthetase activity (GS), (M) intercropping altered the total phosphorus content (TP), (N) intercropping altered the available phosphorus content (AP), and (O) intercropping altered the acid phosphatase activity (ACP).
The ηp2 statistic showed that plant species was the main factor that impacted the OC, βG, NH4_N, NO3_N, Nitrogenase, NOS, GS, TP, AP, ACP, WC, and ANPP, and system type was the main factor that impacted the LAP, the interaction between plant species and system type was the main factor that impacted the MBC, MBN, and NAG (Supplementary Table 1). PCA showed that sampling points were well separated between plant species (ANOSIM r = 0.41, p = 0.046) and among groups (ANOSIM r = 0.38, p < 0.001), but sampling points were not separated between system type (ANOSIM r = 0.01, p = 0.352), indicating that plant species was the main factor that impacted the rhizosphere functions (Figure 2A). The main contributors to these separations were WC, MBN, AP, and MBC. The sole and intercropped Z. mays had higher WC, TP, AP, and OC, the sole and intercropped S. davidii had higher NO3_N, Nitrogenase, NOS, and GS. Correlation network analysis showed that intercropping increased the complexity of functional relationships and the frequency of functional synergy, while decreased the frequency of functional trade-off compared to monoculture (Figures 2B–E). The trade-off analysis showed that the intensity of ecosystem function trade-off was highly variable when low number of paired functions was included, but tends to be stable as the number of included paired functions increased, indicated that it is necessary to consider sufficient number of function pairs when evaluating the strength of function trade-offs (Figures 3A–D). Wilcoxon test showed that intercropped Z. mays significantly decreased the trade-off intensity compared to sole Z. mays, trade-off intensity under intercropped S. davidii was significantly higher than under intercropped Z. mays (Figures 3E–H). Multifunctionality analysis showed intercropping significantly increased the CCMF, and PCMF and AEMF (p = 0.032, 0.0079, and 0.0079) of the S. davidii rhizosphere; however, it did not change the NCMF and ANPP of S. davidii (Wilcoxon test, p = 0.056 – 0.22, Figure 4). Additionally, intercropping significantly decreased the ANPP of Z. mays (Wilcoxon test, p = 0.008), but did not alter the CCMF, NCMF, and PCMF and AEMF of Z. mays rhizosphere (Wilcoxon test, p = 0.15 – 0.84, Figure 4).
Figure 2. (A) Principal component analysis of multiple rhizosphere functions and in the different groups, showing overall intuitive distribution of multiple functions and aboveground net primary productivity, and the contribution of each function to this distribution. (B–E) Relationships among individual rhizosphere functions and aboveground net primary productivity (ANPP). MZ, monoculture Zea mays; IZ, intercropping Zea mays; MS, monoculture Sophora davidii; IZ, intercropping Sophora davidii.
Figure 3. (A–D) Boxplots indicating that the trade-off intensity of functions changed with the number of paired functions. (E–H) Difference of the trade-off intensity between groups at the number of paired functions (Wilcoxon test at α = 0.05). MZ, monoculture Zea mays; IZ, intercropping Zea mays; MS, monoculture Sophora davidii; IZ, intercropping Sophora davidii.
Figure 4. Intercropping altered rhizosphere C, N, P circling multifunctionality (CCMF, NCMF, and PCMF), average rhizosphere ecosystem multifunctionality (AEMF) and aboveground net primary productivity (ANPP). “ns,” p > 0.05; “*,” p < 0.05; “**,” p < 0.01; MZ, monoculture Zea mays; IZ, intercropping Zea mays; MS, monoculture Sophora davidii; IZ, intercropping Sophora davidii.
Intercropping shifted the diversity of rhizosphere microbes
The average Good coverage for all samples was >98.7%, suggesting that the sampling was adequate (Jin et al., 2022a). Intercropping made whole bacterial richness, whole bacterial Shannon diversity, locally abundant bacterial Simpson diversity, regionally rare bacterial richness, regionally rare bacterial Shannon diversity, regionally abundant fungal richness evidently distinct between two plants (ANOVA p < 0.05), made no differences between two plants for locally rare fungal richness, and whole fungal richness (ANOVA p > 0.05) compared to monoculture. Intercropping increased the whole bacterial richness and Shannon diversity of the Z. mays rhizosphere but decreased the whole bacterial Shannon diversity, whole fungal Shannon and Simpson diversity, locally abundant fungal richness and Simpson diversity, and regionally abundant fungal Shannon and Simpson diversity of the Z. mays rhizosphere. Intercropping decreased the whole fungal Shannon diversity, locally abundant fungal richness, and regionally abundant fungal richness and Shannon diversity of S. davidii rhizosphere (ANOVA p < 0.05, Supplementary Table 2). Intercropping made no differences between two plants for locally abundant fungal FUNGuild functional richness and regionally rare fungal FUNGuild functional Simpson diversity (ANOVA p > 0.05). Intercropping increased regionally rare fungal FUNGuild functional Shannon and Simpson diversity of Z. mays rhizosphere; however, intercropping decreased whole fungal FUNGuild functional Shannon diversity and whole bacterial Tax4Fun functional Shannon and Simpson diversity of the Z. mays rhizosphere, increased regionally rare fungal FUNGuild functional Simpson diversity of the S. davidii rhizosphere. Intercropping also decreased whole bacterial FAPROTAX functional Simpson diversity (ANOVA p < 0.05, Supplementary Table 2).
Intercropping shifted the taxonomic and functional composition of rhizosphere microbes
The NMDS and ANOSIM test showed bacterial species, fungal species, and fungal functions significantly differed among groups in each subcommunity (p < 0.01, Figure 5), and bacterial functions significantly differed among groups in each subcommunity (p < 0.05), except for the locally abundant species subcommunity. LEfSe analysis showed that, in the Z. mays rhizosphere, four bacterial clades and 36 fungal phylotypes were evidently enriched under monoculture (p < 0.05, Supplementary Figures 3A,B), and 18 bacterial phylotypes and one fungal phylotype were evidently enriched under intercropping. In the S. davidii rhizosphere, 10 bacterial phylotypes and nine fungal phylotypes were evidently enriched under monoculture, and five bacterial phylotypes and four fungal phylotypes were evidently enriched under intercropping (p < 0.05, Supplementary Figures 3C,D). Across groups, one bacterial phylotype and two fungal phylotypes were evidently enriched in the S. davidii rhizosphere under monoculture, six bacterial phylotypes and one fungal phylotype were evidently enriched in the S. davidii rhizosphere under intercropping, two bacterial phylotypes and 13 fungal phylotypes were evidently enriched in the Z. mays rhizosphere under monoculture, and one bacterial phylotype and no fungal phylotype were evidently enriched in the Z. mays rhizosphere under intercropping (p < 0.05). We also examined any significant changes in species or function among groups using the Kruskal–Wallis rank sum test. A total of 1,383 bacterial species, 34 bacterial functions, 314 fungal species, and 16 fungal functions with significant changes were detected (p < 0.05), 0–50.0% of which were locally abundant species, 30.6–50.0% were locally rare species, 4.9–68.8% were regionally abundant species, and 8.3–29.4% were regionally rare species.
Figure 5. Non-metric multidimensional scaling (NMDS) ordinations based on Bray–Curtis distance matrices of taxonomy (A,B) and functions (C,D) of bacterial (A,C) and fungal (B,D) subcommunities for all samples (n = 20). A stress <0.1 indicated that the result of NMDS is reliable. ANPP, aboveground net primary productivity; MZ, monoculture Zea mays; IZ, intercropping Zea mays; MS, monoculture Sophora davidii; IZ, intercropping Sophora davidii. R and p were reported by ANOSIM.
The networks of the rare subcommunity and whole community have higher complexity than that of the abundant (Supplementary Figures 4, 5). For Z. mays, intercropping increased the number of edges, vertices, and positive edges of whole, locally rare, regionally rare bacteria networks, but decreased those of locally abundant bacteria networks, and for S. davidii, intercropping decreased those of whole, locally abundant, locally rare, regionally rare bacteria networks. This indicated that the effect of intercropping on the network of bacterial subcommunity was regulated by plant species. However, consistent decline in the edges, vertices, and positive edges by intercropping were observed in most of fungal subcommunity, regardless of plant species (Table 1). These results indicated that intercropping increased the complexity and positive links of rare bacteria in Z. mays rhizosphere, but decreased the complexity and positive links of rare in S. davidii rhizosphere, and the complexity and positive links of fungi in both intercropped plants rhizosphere.
Table 1. Properties of subcommunity microbial association network in different groups.
Rhizosphere rare taxa were related to rhizosphere soil functions and aboveground net primary productivity, and changed trade-offs
The Spearman correlation analysis showed different species, functions or diversity indices were significantly related to CCMF, NCMF, and PCMF, AEMF, and ANPP (p < 0.05), respectively. Mantel test and Hierarchical cluster analysis showed three clusters of factors that can impact the system functions (Figure 6). A higher r indicates a stronger relationship. Significant changes in species of locally rare bacteria (Mantel r = 0.8185, p = 0.0001) had a stronger relationship with ANPP than others did (Mantel r = 0.2107 – 0.8024, p = 0.0001–0.0228, Figure 6). Significant changes in species of locally rare fungus (Mantel r = 0.5069, p = 0.0001) had a stronger relationship with PCMF than others did (Mantel r = −0.0127 – 0.4767, p = 0.0001 – 1). Significant changes in species of locally rare bacteria (Mantel r = 0.5340, p = 0.0001) had a stronger relationship with NCMF than others did (Mantel r = 0 – 0.5118, p = 0.0001 – 1). Significant changes in species of locally rare fungus (Mantel r = 0.3295, p = 0.009) had a stronger relationship with CCMF than others did (Mantel r = −0.0800 – 0.2890, p = 0.0092 – 1). Significant changes in species of locally rare fungus (Mantel r = 0.3984, p = 0.0017) had a stronger relationship with AEMF than others did (Mantel r = 0.0079 – 0.2604, p = 0.0151 – 1). Significant changes in species of locally rare bacteria (Mantel r = 0.6026, p = 0.0001) had a stronger relationship with trade-offs intensity (Mantel r = 0.1983 – 0.5962, p = 0.0001 – 1). The significant changes in rhizosphere species and functions had a stronger relationship with CCMF (Mantel r = 0.2776 vs. 0, Kruskal–Wallis rank sum test p = 0.0277), NCMF (Mantel r = 0.3478 vs. 0.3285, Kruskal–Wallis rank sum test p = 0.9199), PCMF (Mantel r = 0.4325 vs. 0, Kruskal–Wallis rank sum test p = 0.0119), AEMF (Mantel r = 0.2786 vs. 0, Kruskal–Wallis rank sum test p = 0.0277), ANPP (Mantel r = 0.5074 vs. 0.2904, Kruskal–Wallis rank sum test p = 0.0365), and trade-offs (Mantel r = 0.3931 vs. 0.2084, Kruskal–Wallis rank sum test p = 0.0036) than the significantly diversity indices.
Figure 6. Heatmap showing the intensity of the selected species (or functions) or the diversity indices impacting rhizosphere functions and aboveground net primary productivity (ANPP). The intensity was indicated by r of Mantel test. The selected species (or functions) were the intersection of species or functions that significantly differed among the groups (MZ, IZ, MS, and IS) and were significantly correlated with system function (C/N/P-circling multifunctionality, AEMF, and ANPP, respectively). Only the significant intensity was showed. Hierarchical cluster analysis showed the clusters of factors. CCMF, carbon-circling multifunctionality; NCMF, nitrogen-circling multifunctionality; PCMF, phosphorus-circling multifunctionality; AEMF, average ecosystem multifunctionality.
Discussion
The selection effect was stronger than the complementary effect on multiple functions
In agreement with findings from macroecology (Karolína et al., 2014; Baldrighi et al., 2017; Fanin et al., 2017; Bagousse-Pingueta et al., 2019) and bulk soils (Bastida et al., 2016; Delgado-Baquerizo et al., 2016; Mori et al., 2016), the diversities of the microbes in the rhizosphere played an essential role in driving multiple functions (Supplementary Table 3). This indicated that maintaining functioning needs protection and use of rhizosphere microbial diversity. Furthermore, the positive effects dominated CCMF, but the negative effects dominated AEMF. The positive and negative effects on NCMF and ANPP were detected. The negative effects (negative correlation, 15/17) were more frequent than the positive effects (positive correlation, 2/17) in driving NCMF (Supplementary Table 3), indicating that the negative effects dominated the NCMF. However, the positive effects (positive correlation, 8/9) were more frequent than the negative effects (negative correlation, 1/9) in driving ANPP (Supplementary Table 3), suggested that the positive effects dominated ANPP.
Here, we showed for the first time that (Figure 5) the significant changes in rhizosphere species and functions had a stronger relationship with CCMF, NCMF, PCMF, AEMF, ANPP, and trade-offs than the significantly diversity indices, indicating that the selection effect played a chief role in driving multiple functions. This finding was different from the findings of studies in the macroecology (Tonin et al., 2018; Woodcock et al., 2019; Chun et al., 2020) and bulk soils (Wen et al., 2020), which indicated that the complementarity effect played a chief role in driving multiple functions whereas the selection effect had a limited role (Emily et al., 2018; Li et al., 2020a). This inconsistency indicated that the selection effect in the rhizosphere could not be ignored and might be a novel way to regulate productivity. Our result could be explained by the following potential causes: (1) The rhizosphere community was under strong selective pressure for special microbes (Poole et al., 2018) based on crucial functions associated with the metabolism of N and P, which are associated with plant growth promotion and nutrition (Mendes et al., 2014), therefore, plant productivity relied on the microbiome for the uptake of nutrients (Chen et al., 2020b). (2) The rhizosphere microbes mediated the root exudation (Korenblum et al., 2020), changed the nutrient supply and absorption of plants and rhizosphere NCMF. In turn, rhizosphere NCMF negatively modulated plant productivity (Spearman r = −0.70, p < 0.001, Supplementary Figure 2). Collectively, the strong selectivity of plants to rhizosphere microorganisms (Philippot et al., 2013) determined which microorganisms or functions appear; therefore, the multiple functions of the rhizosphere and ANPP were more dependent on the selection effect than the complementary effect, indicating that community assembly processes might determine the diversity and composition of community, and it was the result of these processes that determined how selection and complementary effects occur (Leibold et al., 2017). Therefore, our study provided new insight into the rhizosphere that differ from bulk soils.
Rare rhizosphere taxa might contribute over proportionately to multiple functions
Consistent with our previous findings from bulk soils (Ding and Wang, 2021), distinct microbes (bacteria or fungi) dominated the given functions. In this study, fungi likely dominated the CCMF and PCMF, bacteria likely dominated the NCMF and ANPP (Figure 6). This might be due to bacteria and fungus had different metabolic niches (Ding and Wang, 2021). For instance, fungus rather than most bacteria could secrete lignin-degrading enzymes (Ding et al., 2020a), enhanced rhizodeposition, and suppressed the organic matter degradation (Zhou et al., 2020a). Fungi dominated P uptake in symbiotic plants (Wang et al., 2020a; Jiang et al., 2021). The denitrification, mineralization, and assimilation of N were mainly driven by bacteria (Starke et al., 2016; Li et al., 2019; Ding and Wang, 2021). N is the major nutrient limiting plant growth (Moreau et al., 2019; Wang et al., 2020a), thus, if bacteria dominated the NCMF, bacteria were expected to dominate ANPP (Figure 6). Furthermore, this was the first study to report on locally rare bacteria and fungi were overlooked keystone taxa shaping ecosystem functions and trade-offs. These findings supported established results (Radchuk et al., 2016; Jousset et al., 2017; Dee et al., 2019; Zhang et al., 2019; Chen et al., 2020a).
Rare taxa could act as keystone species through several mechanisms. We grouped the potential mechanisms accounting for the greater contribution of the rare species to ecosystem functions into three pathways that may operate simultaneously: (1) Rare species occupied the majority of species in ecosystems (Jain et al., 2014; Lynch and Neufeld, 2015; Jousset et al., 2017; Zhang et al., 2019; Guo et al., 2020) and had extremely high diversity (Mo et al., 2018). Rare taxa had more stronger impact on the multifaceted diversity of the community than abundant taxa (Zhou et al., 2019a). In this study, the species richness of locally rare bacteria was 433- to 677-fold that of locally abundant bacteria, the species richness of regionally rare bacteria was six–seven fold that of regionally abundant bacteria, the species richness of locally rare fungus was 20- to 52-fold that of locally abundant fungus, the species richness of regionally rare bacteria was 1.2- to 1.8- fold that of regionally abundant fungus (Supplementary Table 2). On the one hand, rare species provided insurance effects (Jousset et al., 2017; Chen et al., 2020a) as implied by the insurance hypothesis (Jiao et al., 2017), and recruitment from the persistent rare microbial seed bank provided a broad reservoir of ecological function (Lynch and Neufeld, 2015). On the other hand, most distinct traits combinations were supported predominantly by rare species (David et al., 2013). Similar observations were detected in our study. Both locally and regionally, the function richness of rare bacteria and fungus occupied 57–100% of that of whole bacteria and fungus (Supplementary Table 2). According to the hypothesis of complementary and selection hypotheses (Mensah et al., 2020), high diversity elevated the chance of rare taxa to contributing to ecosystem functionality; in this sense, rare species had the potential to change an ecosystem’s multifunctionality via the enhancement of biodiversity (Angelini et al., 2015; Jousset et al., 2017). (2) Many studies had supported that “being different” was crucial for the influence of rare taxa/. In plant communities, rare species had a higher effect on ecosystem functioning because the of the rare taxa individual mass is higher than that of the abundant taxa (Radchuk et al., 2016). Similar effects had been observed in decomposition systems (Guo et al., 2020) and microbial system (Wei et al., 2019). Essentially, microbes changed system functions through their metabolic functions (Ding and Wang, 2021). In this study, there were obvious differences in functions (Figures 5C,D) and functional diversity (Supplementary Table 2) of the rare versus abundant species, suggested that rare species had the potential to act significant roles in ecosystem functioning via providing different functions. (3) Affecting the interactions. The effect of some taxa on ecosystem functions were not independent from their interactions with other taxa (Wagg et al., 2019). The disappearance of rare caused an obvious or bad alteration in the composition or function of the community (Nannipieri et al., 2020). In contrast, the occurrence of rare could have a good effect on the community (Xue et al., 2020). They could create circumstances that supported the co-occurrence of high densities of different functional organisms, thereby, enhancing MF (Angelini et al., 2015). Rare microbes could also heighten the role of abundant microbes (Jousset et al., 2017). Rare dissimilated low content compounds into materials needed by other microbes or synthesized effective bioactive compounds (Harrison et al., 2021). The rare could reshape rhizosphere community, thereby promote crop growth (Li et al., 2020b). Rare taxa could affect species interactions (Xue et al., 2020). In this study, rare contributed more to the microbe’s positive interaction than abundant (Table 1 and Supplementary Figures 4, 5). Higher frequent facilitations than competitions possibly yielded the complementarity effects (Ding and Wang, 2021). Therefore, rare contributed to ecosystem functions through species interactions (Jousset et al., 2017; Dee et al., 2019). Since rare species are most vulnerable to be lost (Dee et al., 2019) and are largely unexplored, we ascribed great importance to the rare species and suggested to optimize their taxa for maintaining a high ecosystem functionality.
Rare rhizosphere taxa might contribute to trade-offs of multiple functions
Monoculture had higher ANPP than intercropping, this implied that monoculture brought higher benefits than intercropping, which was consistent with the recent findings in tropics (Grass et al., 2020). However, intercropping also shifted the ANPP-CCMF relationship from being none to being positive for Z. mays (Figures 2B,C), shifted the ANPP-NCMF relationship from being negative to being positive for S. davidii (Figures 2D,E). Intercropping decreased the trade-offs intensity compared to sole for Z. mays (Figure 3E). MF and yields were not always synergistic (Supplementary Figure 2), confirming recent findings (Garland et al., 2021). Our study also suggested, for the first time, that the locally rare bacteria species were most strongly related to the trade-offs of multiple functions, and indicated that the trade-offs would likely be reduced by optimizing the taxa of locally rare bacteria.
Conclusion
This study investigated the rhizosphere abundant and rare bacteria and fungi, rhizosphere C/N/P-cycling multifunctionality, ecosystem multifunctionality, aboveground net primary productivity and trade-offs in the Z. mays and S. davidii sole and Z. mays/S. davidii intercropping ecosystems. Results demonstrated that intercropping altered multiple ecosystem functions individually and simultaneously. Intercropped Z. mays significantly decreased the trade-off intensity compared to sole Z. mays, the trade-off intensity under intercropped S. davidii was significantly higher than under intercropped Z. mays. Moreover, both rhizosphere abundant and rare could predict and might affect rhizosphere elements circling, multifunctionality, aboveground productivity and trade-offs, whereas, the significant changes in species of locally rare microbes were the best predictor of rhizosphere elements circling, multifunctionality, aboveground productivity and trade-offs. We thus ascribe great importance to the rare species. Indeed, our results may help in driving a high functionality by directing future efforts on collection, conservation and manipulation of rhizosphere rare species. Further research with more ecosystems and operation of rare microbe combinations will facilitate to maintain a higher ecosystem function and a better understanding of cause-and-effect mechanisms.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
PW and LD: conceptualization. PW: methodology and writing—original draft preparation. LD: data curation and visualization. CZ, LD, YZ, PW, and MW: investigation and writing—reviewing and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by China National Nature Science Foundation (31960341), Guizhou Province Science and Technology Foundation of China (Qiankehejichu-ZK[2021]yiban157), Qiannongkeyuanguojihoubuzhu Project [Qiannongkeyuanguojihoubuzhu (2021)03], and National Grass Industry Technology System (CARS-34).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.985574/full#supplementary-material
Footnotes
- ^ https://github.com/OpenGene/fastp
- ^ http://ccb.jhu.edu/software/FLASH/
- ^ https://www.mothur.org/
- ^ http://drive5.com/usearch
- ^ https://www.arb-silva.de/
- ^ http://unite.ut.ee
- ^ http://www.omicsmart.com
- ^ https://cran.r-project.org/bin/windows/base/old/3.5.3/R-3.5.3-win.exe
- ^ https://gephi.org/
References
Angelini, C., van der Heide, T., Griffin, J. N., Morton, J. P., Derksen-Hooijberg, M., Lamers, L. P., et al. (2015). Foundation species’ overlap enhances biodiversity and multifunctionality from the patch to landscape scale in southeastern United States salt marshes. Proc. Biol. Sci. B 282:20150421. doi: 10.1098/rspb.2015.0421
Aßhauer, K. P., Wemheuer, B., Daniel, R., and Meinicke, P. (2015). Tax4fun: predicting functional profiles from metagenomic 16s rrna data. Bioinformatics 31, 2882–2884. doi: 10.1093/bioinformatics/btv287
Bagousse-Pingueta, Y. L., Soliveresc, S., Grossa, N., Toricesa, R., Berdugoa, M., and Maestre, F. T. (2019). Phylogenetic, functional, and taxonomic richness have both positive and negative effects on ecosystem multifunctionality. PNAS 116, 8419–8424. doi: 10.1073/pnas.1815727116
Baldrighi, E., Giovannelli, D., D’Errico, G., Lavaleye, M., and Manini, E. (2017). Exploring the relationship between macrofaunal biodiversity and ecosystem functioning in the deep sea. Front. Mar. Sci. 4:198. doi: 10.3389/fmars.2017.00198
Bastida, F., Torres, I. F., Moreno, J. L., Baldrian, P., Ondoño, S., Ruiz-Navarro, A., et al. (2016). The active microbial diversity drives ecosystem multifunctionality and is physiologically related to carbon availability in Mediterranean semi-arid soils. Mol. Ecol. 25, 4660–4673. doi: 10.1111/mec.13783
Bowker, M. A., Maestre, F. T., and Mau, R. L. (2013). Diversity and patch-size distributions of biological soil crusts regulate dryland ecosystem multifunctionality. Ecosystems 16, 923–933. doi: 10.1007/s10021-013-9644-5
Byrnes, J. E. K., LarsGamfeldt Isbell, F., Lefcheck, J. S., Griffin, J. N., Hector, A., et al. (2014). Investigating the relationship between biodiversity and ecosystemmultifunctionality: challenges and solutions. Methods Ecol. Evol. 5, 114–124. doi: 10.1111/2041-210X.12143
Chen, Q. L., Ding, J., Zhu, D., Hu, H. W., Delgado-Baquerizo, M., Ma, Y. B., et al. (2020a). Rare microbial taxa as the major drivers of ecosystem multifunctionality in long-term fertilized soils. Soil Biol. Biochem. 141:107686. doi: 10.1016/j.soilbio.2019.107686
Chen, Q. L., Ding, J., Zhu, Y. G., He, J. Z., and Hu, H. W. (2020b). Soil bacterial taxonomic diversity is critical to maintaining the plant productivity. Environ. Int. 140:105766. doi: 10.1016/j.envint.2020.105766
Chen, W., Wang, J., Meng, Z., Xu, R., Chen, J., Zhang, Y., et al. (2020c). Fertility-related interplay between fungal guilds underlies plant richness–productivity relationships in natural grasslands. New Phytol. 226, 1129–1143. doi: 10.1111/nph.16390
Chun, J.-H., Ali, A., and Lee, C.-B. (2020). Topography and forest diversity facets regulate overstory and understory aboveground biomass in a temperate forest of South Korea. Sci. Total Environ. 744:140783. doi: 10.1016/j.scitotenv.2020.140783
Csardi, G., and Nepusz, T. (2006). The igraph software package for complex network research. Inter. J. Complex Systems 1695, 1–9.
David, M., David, R. B., Christopher, B., Jerome, C., Rene, G., Mireille, H.-V., et al. (2013). Rare species support vulnerable functions in high-diversity ecosystems. PLoS Biol. 11:e1001569. doi: 10.1371/journal.pbio.1001569
Dee, L. E., Cowles, J., Isbell, F., Pau, S., Gaines, S. D., and Reich, P. B. (2019). When do ecosystem services depend on rare species? Trends Ecol. Evol. 34, 746–758. doi: 10.1016/j.tree.2019.03.010
Del Galdo Jos, I., Peressotti, A., and Cotrufo, M. F. (2003). Assessing the impact of land-use change on soil C sequestration in agricultural soils by means of organic matter fractionation and stable C isotopes. Global Change Biol. 9, 1204–1213. doi: 10.1046/j.1365-2486.2003.00657.x
Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 7:10541. doi: 10.1038/ncomms10541
Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., and Singh, B. K. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat. Ecol. Evol. 4, 210–220. doi: 10.1038/s41559-019-1084-y
Delgado-Baquerizo, M., Trivedi, P., Trivedi, C., Eldridge, D. J., Reich, P. B., Jeffries, T. C., et al. (2017). Microbial richness and composition independently drive soil multifunctionality. Funct. Ecol. 31, 2330–2343. doi: 10.1111/1365-2435.12924
Ding, L., Shang, Y., Zhang, W., Zhang, Y., Li, S., Wei, X., et al. (2020a). Disentangling the effects of driving forces on soil bacterial and fungal communities under shrub encroachment on the Guizhou Plateau of China. Sci. Total Environ. 709:136207. doi: 10.1016/j.scitotenv.2019.136207
Ding, L., Wang, P., Zhang, W., Zhang, Y., Li, S., Wei, X., et al. (2020b). Soil stoichiometry modulates effects of shrub encroachment on soil carbon concentration and stock in a subalpine grassland. iForest Biogeosci. For. 13, 65–72. doi: 10.3832/ifor3091-012
Ding, L., and Wang, P. (2021). Afforestation suppresses soil nitrogen availability and soil multifunctionality on a subtropical grassland. Sci. Total Environ. 761, 143663. doi: 10.1016/j.scitotenv.2020.143663
Dingha, B. N., Jackai, L. E., Amoah, B. A., and Akotsen-Mensah, C. (2021). Pollinators on cowpea Vigna unguiculata: implications for intercropping to enhance biodiversity. Insects 12:54. doi: 10.3390/insects12010054
Emily, G., Chad, R. Z., Tyler, B., and Lars, A. B. (2018). Ecosystem multifunctionality increases with beta diversity in restored prairies. Oecologia 188, 837–848. doi: 10.1007/s00442-018-4248-6
Fan, K., Cardona, C., Li, Y., Shi, Y., Xiang, X., Shen, C., et al. (2017). Rhizosphere-associated bacterial network structure and spatial distribution differ significantly from bulk soil in wheat crop fields. Soil Biol. Biochem. 113, 275–284. doi: 10.1016/j.soilbio.2017.06.020
Fanin, N., Nilsson, M.-C., Gundale, M. J., Farrell, M., Kardol, P., Ciobanu, M., et al. (2017). Consistent effects of biodiversity loss on multifunctionality across contrasting ecosystems. Nat. Ecol. Evol. 2, 269–278. doi: 10.1038/s41559-017-0415-0
Gagic, V., Bartomeus, I., Jonsson, T., Taylor, A., Winqvist, C., Fischer, C., et al. (2015). Functional identity and diversity of animals predict ecosystem functioning better than species-based indices. Proc. Biol. Sci. 282:20142620. doi: 10.1098/rspb.2014.2620
Gao, D., Pan, X., Khashi u Rahman, M., Zhou, X., and Wu, F. (2021). Common mycorrhizal networks benefit to the asymmetric interspecific facilitation via K exchange in an agricultural intercropping system. Biol. Fertility Soils 57, 959–971. doi: 10.1007/s00374-021-01561-5
Garibotti, I. A., Gonzalez Polo, M., Tabeni, S., and Rasmann, S. (2018). Linking biological soil crust attributes to the multifunctionality of vegetated patches and interspaces in a semiarid shrubland. Funct. Ecol. 32, 1065–1078. doi: 10.1111/1365-2435.13044
Garland, G., Edlinger, A., Banerjee, S., Degrune, F., García-Palacios, P., Pescador, D. S., et al. (2021). Crop cover is more important than rotational diversity for soil multifunctionality and cereal yields in European cropping systems. Nat. Food 2, 28–37. doi: 10.1038/s43016-020-00210-8
George, P. B. L., Lallias, D., Creer, S., Seaton, F. M., Kenny, J. G., Eccles, R. M., et al. (2019). Divergent national-scale trends of microbial and animal biodiversity revealed across diverse temperate soil ecosystems. Nat. commun. 10:1107. doi: 10.1038/s41467-019-09031-1
Gong, X., Shi, J., Zhou, X., Yuan, T., Gao, D., and Wu, F. (2021). Crop rotation with cress increases cucumber yields by regulating the composition of the rhizosphere soil microbial community. Front. Microbiol. 12:631882. doi: 10.3389/fmicb.2021.631882
Grass, I., Kubitza, C., Krishna, V. V., Corre, M. D., Musshoff, O., Putz, P., et al. (2020). Trade-offs between multifunctionality and profit in tropical smallholder landscapes. Nat. Commun. 11:1186. doi: 10.1038/s41467-020-15013-5
Grime, J. P. (1998). Benefits of plant diversity to ecosystems: immediate, fillter and founder effects. J. Ecol. 86, 902–910. doi: 10.1046/j.1365-2745.1998.00306.x
Gu, Z. (2014). circlize implements and enhances circular visualization in R. Bioinformatics 30, 2811–2812. doi: 10.1093/bioinformatics/btu393
Guerra, C. A., Heintz-Buschart, A., Sikorski, J., Chatzinotas, A., Guerrero-Ramirez, N., Cesarz, S., et al. (2020). Blind spots in global soil biodiversity and ecosystem function research. Nat. Commun. 11:3870. doi: 10.1038/s41467-020-17688-2
Guo, C., Cornelissen, J. H. C., Tuo, B., Ci, H., Yan, E. R., and Gomez-Aparicio, L. (2020). Non-negligible contribution of subordinates in community-level litter decomposition: deciduous trees in an evergreen world. J. Ecol. 108, 1713–1724. doi: 10.1111/1365-2745.13341
Harrison, J. G., Beltran, L. P., Buerkle, C. A., Cook, D., Gardner, D. R., Parchman, T. L., et al. (2021). A suite of rare microbes interacts with a dominant, heritable, fungal endophyte to influence plant trait expression. ISME J. 15, 2763–2778. doi: 10.1038/s41396-021-00964-4
He, X., Xie, H., Gao, D., Khashi, U. R. M., Zhou, X., and Wu, F. (2021). Biochar and intercropping with potato-onion enhanced the growth and yield advantages of tomato by regulating the soil properties, nutrient uptake, and soil microbial community. Front. Microbiol. 12:695447. doi: 10.3389/fmicb.2021.695447
Huang, X., su, J., Li, S., Liu, W., and Lang, X. (2019). Functional diversity drives ecosystem multifunctionality in a Pinus yunnanensis natural secondary forest. Sci. Rep. 9:6979. doi: 10.1038/s41598-019-43475-1
Jain, M., Flynn, D. F., Prager, C. M., Hart, G. M., Devan, C. M., Ahrestani, F. S., et al. (2014). The importance of rare species: a trait-based assessment of rare species contributions to functional diversity and possible ecosystem function in tall-grass prairies. Ecol. Evol. 4, 104–112. doi: 10.1002/ece3.915
Jiang, F., Zhang, L., Zhou, J., George, T. S., and Feng, G. (2021). Arbuscular mycorrhizal fungi enhance mineralization of organic phosphorus (P) by carrying bacteria along their extraradical hyphae. New Phytol. 230, 304–315. doi: 10.1111/nph.17081
Jiao, S., Chen, W., Wang, J., Du, N., Li, Q., and Wei, G. (2018). Soil microbiomes with distinct assemblies through vertical soil profiles drive the cycling of multiple nutrients in reforested ecosystems. Microbiome 6, 1–13. doi: 10.1186/s40168-018-0526-0
Jiao, S., Luo, Y., Lu, M., Xiao, X., Lin, Y., Chen, W., et al. (2017). Distinct succession patterns of abundant and rare bacteria in temporal microcosms with pollutants. Environ. Pollut. 225, 497–505. doi: 10.1016/j.envpol.2017.03.015
Jin, X., Bai, Y., Khashi u Rahman, M., Kang, X., Pan, K., Wu, F., et al. (2022a). Biochar stimulates tomato roots to recruit a bacterial assemblage contributing to disease resistance against Fusarium wilt. iMeta 1:e37. doi: 10.1002/imt2.37
Jin, X., Wang, Z., Wu, F., Li, X., and Zhou, X. (2022b). Litter mixing alters microbial decomposer community to accelerate tomato root litter decomposition. Microbiol. Spectrum 10:e0018622. doi: 10.1128/spectrum.00186-22
Jousset, A., Bienhold, C., Chatzinotas, A., Gallien, L., Gobet, A., Kurm, V., et al. (2017). Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 11, 853–862. doi: 10.1038/ismej.2016.174
Karolína, B., Moretti, M., Bello, F. D., Dias, A. T., and Berg, M. P. (2014). Disentangling community functional components in a litter-macrodetritivore model system reveals the predominance of the mass ratio hypothesis. Ecol. Evol. 4, 408–416. doi: 10.1002/ece3.941
Khashi u Rahman, M., Wang, X., Gao, D., Zhou, X., and Wu, F. (2021). Root exudates increase phosphorus availability in the tomato/potato onion intercropping system. Plant Soil 464, 45–62. doi: 10.1007/s11104-021-04935-8
Korenblum, E., Dong, Y., Szymanski, J., Panda, S., Jozwiak, A., Massalha, H., et al. (2020). Rhizosphere microbiome mediates systemic root metabolite exudation by root-to-root signaling. PNAS 117, 3874–3883. doi: 10.1073/pnas.1912130117
Leibold, M. A., Chase, J. M., and Ernest, S. K. M. (2017). Community assembly and the functioning of ecosystems: how metacommunity processes alter ecosystems attributes. Ecology 98, 909–919. doi: 10.1002/ecy.1697
Li, C., Hoffland, E., Kuyper, T. W., Yu, Y., Li, H., Zhang, C., et al. (2020a). Yield gain, complementarity and competitive dominance in intercropping in China: a meta-analysis of drivers of yield gain using additive partitioning. Euro. J. Agron. 113:125987. doi: 10.1016/j.eja.2019.125987
Li, F., Zhang, S., Wang, Y., Li, Y., and Han, Y. (2020b). Rare fungus, Mortierella capitata, promotes crop growth by stimulating primary metabolisms related genes and reshaping rhizosphere bacterial community. Soil Biol. Biochem. 151:108017. doi: 10.1016/j.soilbio.2020.108017
Li, Z., Tian, D., Wang, B., Wang, J., Wang, S., Chen, H. Y. H., et al. (2019). Microbes drive global soil nitrogen mineralization and availability. Global Change Biol. 25, 1078–1088. doi: 10.1111/gcb.14557
Liang, Y., Xiao, X., Nuccio, E. E., Yuan, M., Zhang, N., Xue, K., et al. (2020). Differentiation strategies of soil rare and abundant microbial taxa in response to changing climatic regimes. Environ. Microbiol. 22, 1327–1340. doi: 10.1111/1462-2920.14945
Loreau, M. (2000). Biodiversity and ecosystem functioning: recent theoretical advances. Oikos 91, 3–17. doi: 10.1034/j.1600-0706.2000.910101.x
Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507
Lu, T., Ke, M., Lavoie, M., Jin, Y., Fan, X., Zhang, Z., et al. (2018). Rhizosphere microorganisms can influence the timing of plant flowering. Microbiome 6, 1–12. doi: 10.1186/s40168-018-0615-0
Luo, G., Rensing, C., Chen, H., Liu, M., Wang, M., Guo, S., et al. (2018). Deciphering the associations between soil microbial diversity and ecosystem multifunctionality driven by long-term fertilization management. Funct. Ecol. 32, 1103–1116. doi: 10.1111/1365-2435.13039
Lynch, M. D. J., and Neufeld, J. D. (2015). Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol. 13:217. doi: 10.1038/nrmicro3400
Martin-Guay, M.-O., Paquette, A., Dupras, J., and Rivest, D. (2018). The new green revolution: sustainable intensification of agriculture by intercropping. Sci. Total Environ. 615, 767–772. doi: 10.1016/j.scitotenv.2017.10.024
Mendes, L. W., Kuramae, E. E., Navarrete, A. A., van Veen, J. A., and Tsai, S. M. (2014). Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J. 8, 1577–1587. doi: 10.1038/ismej.2014.17
Mensah, S., du Toit, B., and Seifert, T. (2018). Diversity–biomass relationship across forest layers: implications for niche complementarity and selection effects. Oecologia 187, 783–795. doi: 10.1007/s00442-018-4144-0
Mensah, S., Valère Salako, K., Assogbadjo, A., Glèlè Kaka, R., Brice, S., and Seifert, T. (2020). Functional trait diversity is a stronger predictor of multifunctionality than dominance: evidence from an afromontane forest in south africa. Ecol. Indic. 115:106415. doi: 10.1016/j.ecolind.2020.106415
Merino-Martín, L., Stokes, A., Gweon, H. S., Moragues-Saitua, L., Staunton, S., Plassard, C., et al. (2021). Interacting effects of land use type, microbes and plant traits on soil aggregate stability. Soil Biol. Biochem. 154:108072. doi: 10.1016/j.soilbio.2020.108072
Mo, Y., Zhang, W., Yang, J., Lin, Y., Yu, Z., and Lin, S. (2018). Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes. ISME J. 12, 2198–2210. doi: 10.1038/s41396-018-0153-6
Moreau, D., Bardgett, R. D., Finlay, R. D., Jones, D. L., and Philippot, L. (2019). A plant perspective on nitrogen cycling in the rhizosphere. Funct. Ecol. 33, 540–552.
Mori, A. S., Isbell, F., Fujii, S., Makoto, K., Matsuoka, S., and Osono, T. (2016). Low multifunctional redundancy of soil fungal diversity at multiple scales. Ecol. Lett. 19, 249–259. doi: 10.1111/ele.12560
Nannipieri, P., Ascher-Jenull, J., Ceccherini, M. T., Pietramellara, G., Renella, G., and Schloter, M. (2020). Beyond microbial diversity for predicting soil functions: a mini review. Pedosphere 30, 5–17. doi: 10.1016/s1002-0160(19)60824-6
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). Funguild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Niklaus, P. A., Baruffol, M., Jin-Sheng, H., Ma, K., and Schmid, B. (2017). Can niche plasticity promote biodiversity–productivity relationships through increased complementarity? Ecology 98, 1104–1116. doi: 10.1002/ecy.1748
Philippot, L., Raaijmakers, J. M., Lemanceau, P., and van der Putten, W. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Poole, P., Ramachandran, V., and Terpolilli, J. (2018). Rhizobia: from saprophytes to endosymbionts. Nat. Rev. Microbiol. 16, 291–303. doi: 10.1038/nrmicro.2017.171
Radchuk, V., De Laender, F., Van den Brink, P. J., and Grimm, V. (2016). Biodiversity and ecosystem functioning decoupled: invariant ecosystem functioning despite non-random reductions in consumer diversity. Oikos 125, 424–433. doi: 10.1111/oik.02220
Säterberg, T., Jonsson, T., Yearsley, J., Berg, S., and Ebenman, B. (2019). A potential role for rare species in ecosystem dynamics. Sci. Rep. 9:11107. doi: 10.1038/s41598-019-47541-6
Shao, S., Chen, M., Liu, W., Hu, X., Wang, E. T., Yu, S., et al. (2020). Long-term monoculture reduces the symbiotic rhizobial biodiversity of peanut. Syst. Appl. Microbiol. 43:126101. doi: 10.1016/j.syapm.2020.126101
Shi, X., Wang, J., Lucas- Borja, M. E., Wang, Z., Li, X., and Huang, Z. (2021). Microbial diversity regulates ecosystem multifunctionality during natural secondary succession. J. Appl. Ecol. 58, 2833–2842. doi: 10.1111/1365-2664.14015
Starke, R., Kermer, R., Ullmann-Zeunert, L., Baldwin, I. T., Seifert, J., Bastida, F., et al. (2016). Bacteria dominate the short-term assimilation of plant-derived N in soil. Soil Biol. Biochem. 96, 30–38. doi: 10.1016/j.soilbio.2016.01.009
Tonin, A. M., Pozo, J., Monroy, S., Basaguren, A., Perez, J., and Goncalves, J. F. Jr., et al. (2018). Interactions between large and small detritivores influence how biodiversity impacts litter decomposition. J. Anim. Ecol. 87, 1465–1474. doi: 10.1111/1365-2656.12876
van der Plas, F. (2019). Biodiversity and ecosystem functioning in naturally assembled communities. Biol. Rev. 94, 1220–1245. doi: 10.1111/brv.12499
Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E. E., and van der Heijden, M. G. A. (2019). Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat. Commun. 10:4841. doi: 10.1038/s41467-019-12798-y
Wang, S., Chen, A., Xie, K., Yang, X., Luo, Z., Chen, J., et al. (2020a). Functional analysis of the OsNPF4.5 nitrate transporter reveals a conserved mycorrhizal pathway of nitrogen acquisition in plants. PNAS 117, 16649–16659. doi: 10.1073/pnas.2000926117
Wang, X., Li, F. Y., Wang, Y., Liu, X., Cheng, J., Zhang, J., et al. (2020b). High ecosystem multifunctionality under moderate grazing is associated with high plant but low bacterial diversity in a semi-arid steppe grassland. Plant Soil 448, 265–276. doi: 10.1007/s11104-020-04430-6
Wang, J., Shi, X., Lucas-Borja, M. E., Lam, S. K., Wang, Z., and Huang, Z. (2022). Plants, soil properties and microbes directly and positively drive ecosystem multifunctionality in a plantation chronosequence. Land Degrad Dev. 33, 1–9. doi: 10.1002/ldr.4371
Wei, X., Hu, Y., Razavi, B. S., Zhou, J., Shen, J., Nannipieri, P., et al. (2019). Rare taxa of alkaline phosphomonoesterase-harboring microorganisms mediate soil phosphorus mineralization. Soil Biol. Biochem. 131, 62–70. doi: 10.1016/j.soilbio.2018.12.025
Wen, Z., Zheng, H., Zhao, H., Xie, S., Liu, L., and Ouyang, Z. (2020). Land-use intensity indirectly affects soil multifunctionality via a cascade effect of plant diversity on soil bacterial diversity. Glob. Ecol. Conserv. 23:e01061. doi: 10.1016/j.gecco.2020.e01061
Woodcock, B. A., Garratt, M. P. D., Powney, G. D., Shaw, R. F., Osborne, J. L., Soroka, J., et al. (2019). Meta-analysis reveals that pollinator functional diversity and abundance enhance crop pollination and yield. Nat. Commun. 10:1481. doi: 10.1038/s41467-019-09393-6
Xiong, C., He, J. Z., Singh, B. K., Zhu, Y. G., Wang, J. T., Li, P. P., et al. (2020). Rare taxa maintain the stability of crop mycobiomes and ecosystem functions. Environ. Microbiol. 23, 1907–1924. doi: 10.1111/1462-2920.15262
Xue, M., Guo, Z., Gu, X., Gao, H., Weng, S., Zhou, J., et al. (2020). Rare rather than abundant microbial communities drive the effects of long-term greenhouse cultivation on ecosystem functions in subtropical agricultural soils. Sci. Total Environ. 706:136004. doi: 10.1016/j.scitotenv.2019.136004
Zhang, Y., Dong, S., Gao, Q., Ganjurjav, H., Wang, X., and Geng, W. (2019). “Rare biosphere” plays important roles in regulating soil available nitrogen and plant biomass in alpine grassland ecosystems under climate changes. Agr. Ecosyst. Environ. 279, 187–193. doi: 10.1016/j.agee.2018.11.025
Zheng, Q., Hu, Y., Zhang, S., Noll, L., and Wanek, W. (2019). Soil multifunctionality is affected by the soil environment and by microbial community composition and diversity., 136, 107521-. Soil Biol. Biochem. 136, 107521. doi: 10.1016/j.soilbio.2019.107521
Zhong, L., Wang, J., Zhang, X., and Ying, L. (2020). Effects of agricultural land consolidation on ecosystem services: trade-offs and synergies. J. Cleaner Product. 264:121412. doi: 10.1016/j.jclepro.2020.121412
Zhou, J., Zang, H., Loeppmann, S., Gube, M., Kuzyakov, Y., and Pausch, J. (2020a). Arbuscular mycorrhiza enhances rhizodeposition and reduces the rhizosphere priming effect on the decomposition of soil organic matter. Soil Biol. Biochem. 140:107641. doi: 10.1016/j.soilbio.2019.107641
Zhou, Z., Wang, C., and Luo, Y. (2020b). Meta-analysis of the impacts of global change factors on soil microbial diversity and functionality. Nat. Commun. 11:3072. doi: 10.1038/s41467-020-16881-7
Zhou, Q., Zhang, X., He, R., Wang, S., Jiao, C., Huang, R., et al. (2019a). The composition and assembly of bacterial communities across the rhizosphere and phyllosphere compartments of phragmites Australis. Diversity 11:98. doi: 10.3390/d11060098
Keywords: abundant taxa, rare taxa, diversity, trade-off, multifunctionality
Citation: Wang P, Ding L, Zou C, Zhang Y and Wang M (2022) Rhizosphere element circling, multifunctionality, aboveground productivity and trade-offs are better predicted by rhizosphere rare taxa. Front. Plant Sci. 13:985574. doi: 10.3389/fpls.2022.985574
Received: 04 July 2022; Accepted: 25 August 2022;
Published: 08 September 2022.
Edited by:
Xingang Zhou, Northeast Agricultural University, ChinaReviewed by:
Sai Shiva Krishna Prasad Vurukonda, University of Modena and Reggio Emilia, ItalyCopyright © 2022 Wang, Ding, Zou, Zhang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Leilei Ding, peterding2007gy@163.com
†These authors have contributed equally to this work and share first authorship