# Nonlinear Effects Induced by Interactions among Functional Groups of Bacteria and Fungi Regulate the Priming Effect in Malagasy Soils

Benoît Jaillard, Kanto Razanamalala, Cyrille Violle, Laetitia Bernard

► To cite this version:

Benoît Jaillard, Kanto Razanamalala, Cyrille Violle, Laetitia Bernard. Nonlinear Effects Induced by Interactions among Functional Groups of Bacteria and Fungi Regulate the Priming Effect in Malagasy Soils. Microorganisms, 2023, 11, pp.1106. 10.3390/microorganisms11051106. hal-04130630

HAL Id: hal-04130630

https://hal.inrae.fr/hal-04130630

Submitted on 16 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Nonlinear Effects Induced by Interactions among Functional Groups of Bacteria and Fungi Regulate the Priming Effect in Malagasy Soils

Benoît Jaillard<fnref n="1" />, Kanto Razanamalala<fnref n="2" />, Cyrille Violle<fnref n="3" /> and Laetitia Bernard <fnref n="1" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Jaillard, B.; Razanamalala, K.; Violle, C.; Bernard, L. Nonlinear Effects Induced by Interactions among Functional Groups of Bacteria and Fungi Regulate the Priming Effect in Malagasy Soils.
*Microorganisms* **2023**, *11*, 1106.
https://doi.org/10.3390/microorganisms11051106

Academic Editor: Volker Brozel

Received: 21 March 2023

Revised: 11 April 2023

Accepted: 12 April 2023

Published: 23 April 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The priming effect (PE) occurs when fresh organic matter (FOM) supplied to soil alters the rate of decomposition of older soil organic matter (SOM). The PE can be generated by different mechanisms driven by interactions between microorganisms with different live strategies and decomposition abilities. Among those, stoichiometric decomposition results from FOM decomposition, which induces the decomposition of SOM by the release of exoenzymes by FOM-decomposers. Nutrient mining results from the co-metabolism of energy-rich FOM with nutrient-rich SOM by SOM-decomposers. While existing statistical approaches enable measurement of the effect of community composition (linear effect) on the PE, the effect of interactions among co-occurring populations (non-linear effect) is more difficult to grasp. We compare a non-linear, clustering approach with a strictly linear approach to separately and comprehensively capture all linear and non-linear effects induced by soil microbial populations on the PE and to identify the species involved. We used an already published data set, acquired from two climatic transects of Madagascar Highlands, in which the high-throughput sequencing of soil samples was applied parallel to the analysis of the potential capacity of microbial communities to generate PE following a ¹³C-labeled wheat straw input. The linear and clustering approaches highlight two different aspects of the effects of microbial biodiversity on SOM decomposition. The comparison of the results enabled identification of bacterial and fungal families, and combinations of families, inducing either a linear, a non-linear, or no effect on PE after incubation. Bacterial families mainly favoured a PE proportional to their relative abundances in soil (linear effect). Inversely, fungal families induced strong non-linear effects resulting from interactions among them and with bacteria. Our findings suggest that bacteria support stoichiometric decomposition in the first days of incubation, while fungi support mainly the nutrient mining of soil's organic matter several weeks after the beginning of incubation. Used together, the clustering and linear approaches therefore enable the estimation of the relative importance of linear effects related to microbial relative abundances, and non-linear effects related to interactions among microbial populations on soil properties. Both approaches also enable the identification of key microbial families that mainly regulate soil properties.

**Keywords:** assembly motif; bacteria; biodiversity; combinatorial analysis; ecosystem functioning; fungi; organic matter; priming effect; soil

## 1. Introduction

The positive priming effect (PE), i.e., the acceleration of the soil's organic matter (SOM) decomposition by the input of a fresh, easily decomposable, organic matter (FOM), is a key ecosystem process [1]. PE has long been explained by two major processes: (i) stoichiometric decomposition, which results from FOM decomposition, which induces the decomposition of SOM by the release of exoenzymes by FOM-decomposers; (ii) nutrient

<footnote n="1">¹ Eco & Sols, University Montpellier, IRD, INRAE, CIRAD, Institut Agro, 2 place Viala, 34060 Montpellier, France</footnote>

<footnote n="2">² Laboratoire des Radio-Isotopes, Antananarivo 101, Madagascar</footnote>

<footnote n="3">³ CEFE, University Montpellier, CNRS, EPHE, IRD, 1919 Route de Mende, 34293 Montpellier, France</footnote>

<footnote>* Correspondence: laetitia.bernard@ird.fr</footnote>
mining results from the co-metabolism of energy-rich FOM with nutrient-rich SOM by SOM-decomposers [2,3]. Energy-rich FOM can also come from organic matter desorbed by carboxylic acids released by FOM-decomposers [4]. The microbial actors and the biotic and abiotic determinants of the priming effect are therefore multiple [5]. All the proposed processes involve interactions among microbial communities, including competitive interactions among microbial populations for access to FOM and complementarity between the enzymatic activities of different microbial functional groups [6]. In order to improve the understanding of the priming effect, therefore, the identification of the microbial populations that favored each of the different processes is needed, either because of their own activities or because of their interactions with other soil microbial populations. DNA-SIP has been described as a direct method to identify microbial populations feeding on ¹³C-labeled organic matter [7]. However, as the PE generated by nutrient mining results from the co-metabolism between different sources of organic matter (labeled and unlabeled), such a method cannot discriminate the true FOM-decomposers implied in stoichiometric decomposition from the SOM-decomposers benefiting from labeled FOM catabolites and implied in "nutrient mining" [6]. Therefore, an appropriate statistical approach could help to better disentangle the involved processes.

Microbial biodiversity generates diverse effects on soil functions. The effects are most often decomposed into interaction effect, sensu stricto, and composition effect [8]. The composition effect reflects the specific activities of certain key species in the realization of a soil function and in the absence of interaction with other co-occurring species. It is therefore assumed to be cumulative and responds linearly to the presence of these species and their relative abundance within the community. The interaction effect gathers all positive and negative diversity effects other than the composition effect: it is therefore, by definition, non-linear. Several authors [8,9] separate the interaction effect from the composition effect on the basis of the linearity of the response of the composition effect to the relative abundance of key species [10]. The observed function of a community is compared with the sum of properties, characterizing all populations composing the community as they were grown in monoculture, weighted by their relative abundance. Using molecular methods, it is possible to discriminate between different populations of an assemblage and evaluate their relative abundance. However, the properties of each population grown in monoculture cannot be known because of the impossibility to isolate and cultivate many soil-native microorganisms. In this context, it is not possible to separate the two biodiversity effects on the basis of the linearity of the response of the composition effect to community biodiversity and species' relative abundance using classical approaches.

Most statistical methods are based on the linearity of the response of a variable to be explained by a set of explanatory variables: variance analysis, linear regression, and principal component analysis. In the absence of interactions among species, a linear analysis of the property of a community explained by its biodiversity, especially species' relative abundances, would therefore approximate the effect of community composition only. Conversely, a non-linear analysis of the community property would approximate the overall biodiversity effect of a community, i.e., the product of the composition effect and the interaction effect. Jaillard et al. [11,12] recently proposed an approach to explain ecosystem property by the co-occurrence of species in an ecosystem. This approach searches for species combinations that best account for observed variations in the ecosystem property. The approach enables the cluster of species into functional groups whose co-occurrence determines the ecosystem property. This clustering approach is non-linear. It can therefore be used to analyze the overall biodiversity effect of a microbial community. Compared with a linear approach, this non-linear approach allows for the inference of the effect of interactions among species on community properties.

Here, we aim to identify microbial taxa that predominantly contribute to the generation of PE in Malagasy Ferralsols, either by their own activities or through their interactions with other microbial taxa. Our approach is based on the comparison between linear and non-linear analyses of the same dataset. Three soil properties associated with SOM
decomposition activity were analyzed: basal soil respiration, mineralization of ¹³C labeled wheat straw supplied to soil, and the PE generated by the straw supplied [13]. PE was calculated by removing the total carbon mineralized in the part resulting from wheat-straw and that which was registered in the non-amended subsample (i.e., basal respiration). These three soil properties were measured after seven and 42 days of incubation, because those times likely corresponded to different processes of PE generation. The composition of bacterial and fungal communities was assessed in pristine soils before incubation, by high throughput molecular analysis of 16S and 18S ribosomal genes, respectively. The gene sequences were assigned at the family level.

## 2. Materials and Methods

### 2.1. Dataset

The dataset used here was already partly published by [13]. Ferralsols from the Highlands of Madagascar were collected (May 2014, 57 sites) under herbaceous savannah dominated by the Poaceae *Aristida* sp., along gradients of temperature and rainfall. At each site, five soil cores were sampled from 0–10 cm depth and mixed to form a composite soil sample. The 57 composite soil samples were sieved to 2 mm, stored at 4 °C, then sent to the Genosol platform in Dijon, France. A 50 g aliquot of each composite soil sample was subsampled, then freeze-dried prior to genetic analysis. Soils have also been characterized as fully described in [13] and characteristics can be consulted in Table S1.

Bacterial and fungal community compositions were assessed on 57 soil samples. Banks of 16S and 18S ribosomal sequences were prepared prior to pyrosequencing analyses as fully described in [14]. Pyrosequencing was carried out on a GS FLX Titanium (Roche 454 Sequencing System). Biocomputing analysis of the sequences was performed using the GnS-PIPE [15] as described in detail in Tardy et al. [16]. In order to compare the datasets efficiently and avoid biased community comparisons, the sample reads were reduced by random selection closed to the lowest datasets (3000 reads for 16S- and 18S-rRNA gene sequences, respectively, for each soil sample). The retained high-quality reads were used for taxonomy-based analysis using similarity approaches and dedicated reference databases from SILVA [17]. Assignments published in Razanamalala et al. [13] were made at the highest taxonomic level, i.e., the phylum level. However, in the present study, we preferred to perform assignments at a more functionally relevant level, i.e., the family level. Actually, we tested the class, order, and family phylogenetic levels, and family appeared as the best compromise to improve functional inference without drastically reducing the dataset due to the lack of known species and even genera in the SILVA database. In order to compare this work with that of Razanamalala et al. [13], the links between families and phyla are given in Table S2. The relative abundance matrices had a dimension of *n* = 57 soils for *p* = 60 bacterial or fungal families, respectively. This number of families was chosen in order to be at the limit of the recommendations of the linear approaches in terms of the number of variables, which must be lower compared with the number of samples (see Section 2.2.3). Relative abundances per soil aliquot ranged from 1 to 1486, and from 5 to 1048, out of 3000 sequences, with a median of 4.0 and 11.5, for bacteria and fungi, respectively (Figure S1). The observation frequency of each family ranges from 32 to 100%, and 35 to 100% of soils, for bacteria and fungi, respectively (Figure S2).

In parallel, the soils were incubated in the absence or presence of ¹³C-enriched wheat straw (for details, see [13]). CO₂ release was measured after seven and 42 days of incubation. Those two times were chosen because they likely corresponded to different processes of PE generation [13]. Three measurements were made: the CO₂ released by soil in the absence of straw, the CO₂ released by soil with straw supplied, and the ¹³C-labeled CO₂ release due to the mineralization of the straw supplied, which was previously enriched with ¹³C. The measurements enabled the determination of three soil properties concerning SOM decomposition activity carried out by microbial communities: (i) basal soil respiration, equal to the release of CO₂ by soil in absence of straw, (ii) straw mineralization, equal to the release of ¹³C-labeled CO₂ due to the mineralization of the straw supplied only,
and (iii) the priming effect (PE) generated by the straw supply, equal to the CO₂ release by soil with straw supplied, minus the CO₂ resulting from soil basal respiration and straw mineralization.

## 2.2. Linear Analyses

### 2.2.1. Modeling SOM Decomposition by the Relative Abundance of Microbial Families

A multi-linear analysis without interactions assumes that each microbial family contributes to SOM decomposition in proportion to its relative abundance. The relative abundance of each family is then assumed to be independent of the relative abundance of other families. A multi-linear analysis therefore models SOM decomposition as a linear combination of the relative abundances of families that co-occur in soil. The model consisted of weighing the relative abundance of each family to optimize the fit to SOM decomposition.

### 2.2.2. Evaluation of the Model Quality

The explanatory power of the model was measured by its coefficient of determination, R², that is the ratio of the variance explained by multi-linear analysis on the total observed variance. The model quality, i.e., the compromise between the model precision and the number of explanatory variables, hereafter, the number of families (plus one for the regression intercept) was measured by its Akaike Information Criterion AICc, corrected for small sample sets [18].

### 2.2.3. Model Reduction by Variable Selection

A multi-linear analysis is improved by reducing the number of explanatory variables. We used a backward stepwise variable selection, based on the progressive removal of microbial families considered as irrelevant explanatory variables because of non-significant contributions to SOM decomposition. Each family was successively removed from the set of explanatory variables, and the effect of each of these removals on model quality was evaluated. The family whose removal maximized the model quality was considered irrelevant and was definitively removed from the set of explanatory variables. We repeated the procedure while the removal of a family increased the model quality. We stopped when removing any remaining families decreased model quality. This means that all remaining families were then identified as relevant explanatory variables, because their contribution to SOM decomposition was significant and because their inclusion in the set of explanatory variables improved model quality. Throughout the variable selection procedure, the quality of the model was measured by its AICc.

A linear model consists of a system of *n* equations with *p* unknowns, hereafter, *n* observed soils for *p* microbial families: it requires a non-zero number of degrees of freedom, i.e., a number *n* of observations greater than the number *p* of families, plus one (*n* > *p* + 1). The bacteria and fungus relative abundance matrices comprised *n* = 57 soils for *p* = 60 families: this mathematical constraint imposed to keep only *p* = 55 families for *n* = 57 soils. When analyzing only the bacterial families, or only the fungal families, we therefore removed the five least abundant families from the bacteria and fungi relative abundance matrices, respectively.

## 2.3. Clustering Analysis

### 2.3.1. Modeling SOM Decomposition by the Co-Occurring of Relative Abundance Classes of Microbial Families

A clustering analysis is based on the occurrence of families in microbial communities, i.e., qualitative data. The presence of each microbial family in soil, however, was measured by their relative abundance, i.e., quantitative data. The relative abundance of each family was therefore converted into classes of relative abundance, then each quantitative matrix of relative abundance was converted into a qualitative one-shot encoding table. Clustering methods are all the more effective when each family to be clustered is equitably observed. The least frequent bacterial and fungal families were observed in 32% and 35% of sampled
soils, respectively. The distribution of relative abundance of each family was therefore segmented into classes of relative abundance, centered on an observation frequency of 30% (Figure S2). The ubiquitous families, i.e., the families observed in 100% of soils, were segmented into three classes of relative abundance, noted 1/3, 2/3, and 3/3 from the least to the most abundant classes. The least frequent families were not segmented and were noted 1/1 as a reminder (Figure S2a,c). Families in between were segmented into two classes and noted as 1/2 and 2/2. The segmentation thresholds were determined in such a way that the numbers of observations in each class were equal to, within one unit. The resulting one-hot encoding tables contained 153 and 175 classes of family relative abundance, with an average observation frequency of 31.7% and 32.4% (median 32.2% and 33.3%) for bacteria and fungi, respectively (Figure S2b,d). The clustering analysis then aimed to group classes of family relative abundance whose co-occurrence best accounted for variations of SOM decomposition. The method then looked to maximize the inter-group variance of SOM decomposition associated with different combinations of functional groups, i.e., the different assembly motifs [11]. As in the linear modelling, the explanatory power of the model was measured by the coefficient of determination, R², that is the ratio of the variance explained by the clustering analysis on the total observed variance.

### 2.3.2. Model Reduction by Variable Selection

As a linear model, a clustering model is improved by reducing the number of explanatory variables. We used a backward stepwise selection similar to the one applied to linear analysis. The first difference was that each family was segmented into several classes of relative abundance: removing a family here therefore meant removing all the classes of relative abundance from the family. A second difference is that a clustering model enables the clustering of families into functional groups and observations into assembly motifs, i.e., combinations of functional groups: we imposed that for a family to be considered as an irrelevant variable, its removal must change neither the functional groups nor the assembly motifs. A third difference is that all the least significant classes of family relative abundance are clustered in the same functional group by default, with the residual group labeled by the letter of highest rank: only the removal of families belonging to this group are therefore tested.

### 2.3.3. Hierarchy of Functional Groups

A clustering analysis ultimately provided a hierarchical tree of the functional groups of classes of family relative abundance that best accounted for SOM decomposition. The functional groups are sorted from left to right, from those that best to least explain SOM decomposition. Within each functional group, the classes of family relative abundance were also sorted from left to right, based on their ability to explain SOM decomposition. The residual group that clusters the least significant classes of family relative abundance is therefore on the right, labeled by the letter of highest rank, and the least significant classes of family relative abundance belonging to this group is on the right of this functional group.

### 2.3.4. Evaluation of the Model Quality

As for a multi-linear model, the explanatory capacity of a clustering model was measured by its coefficient of determination, R², that is the ratio of variance explained by the clustering analysis on the total observed variance. The quality of the model was measured by its AICc, i.e., the trade-off between model precision and the number of explanatory variables, hereafter, the number of families were considered relevant. As explained above, the variable selection procedure was based on constant numbers of functional groups and observed assembly motifs. However, the number of families on which the clustering analysis was based decreased. The number of variables taken into account in the AICc calculation were the family numbers, instead of the numbers of observed assemblage motifs that stayed constant [11].
### 2.3.5. Linear Effects Associated with Each Functional Group

The clustering analysis looks to maximize the inter-group variance of SOM decomposition associated with the different assembly motifs, which are the combinations of functional groups. It is also possible to assess the linear effects explained by clustering analysis using a variance analysis of SOM decomposition explained by each functional group of classes of family relative abundance and their statistical interactions.

### 2.4. Statistical Computations

All computations were performed using R-software [19]. Linear and multi-linear regressions were performed using the lm function, and variance analyses using the aov function of the R-package stats. Clustering analyses were performed using the fclust function of the R-package functClust. Backward stepwise selections were performed using the lm_step and fclust_step functions. Comparison of medians was performed using a non-parametric Kruskal-Wallis test (function kruskal.test of the R-package stats). All statistical tests were performed at the probability threshold of *p* < 0.001. All scripts are available on the INRAE dataverse at https://doi.org/10.15454/LRY3M2 (accessed on 20 March 2023).

## 3. Results

### 3.1. Linear Analyses of the Relationship between SOM Decomposition Activity and Environmental Conditions of Sampling Sites

Basal soil respirations after seven and 42 days of incubation were positively correlated ($R^2$ = 0.773; *p* < 10<sup>-10</sup>). Straw mineralization after 42 days did not depend on that observed after seven days of incubation ($R^2$ = 0.007; *p* = 0.238), and the priming effects after seven and 42 days of incubation were negatively correlated ($R^2$ = 0.203; *p* = 0.256 × 10<sup>-3</sup>). This means that soils with a high PE after seven days had a low PE after 42 days of incubation.

In general, SOM decomposition functions did not vary with precipitation at sampling sites. In contrast, basal soil respirations after seven and 42 days were negatively correlated ($R^2$ = 0.496; *p* = 0.601 × 10<sup>-9</sup> and $R^2$ = 0.391; *p* = 0.118 × 10<sup>-7</sup>, respectively) with temperatures, and PE after seven and 42 days were negatively and positively correlated, respectively ($R^2$ = 0.043; *p* = 2.064 × 10<sup>-5</sup>; and $R^2$ = 0.180; *p* = 0.588 × 10<sup>-3</sup>), with temperatures at the sample sites.

### 3.2. Multi-Linear and Clustering Analyses of the Relationship between Initial Composition of Bacterial and Fungal Communities and SOM Decomposition Activity

Both multi-linear and clustering approaches to the relative abundance of bacterial families accurately explained SOM decomposition functions, i.e., basal soil respiration, straw mineralization, and the priming effect, after seven and 42 days of incubation (Table 1). The six multi-linear models selected, in median, 26 key families out of 55, with model degrees of freedom of 28 and $R^2$ of 0.959. The Fisher ratios were around 25.0, all associated with a probability lower than 10<sup>-10</sup>. The median AICc of the six models was -316.6. The clustering analyses selected, in median, only 14 key families out of 60, with model degrees of freedom of 46 and $R^2$ of 0.965. The Fisher ratio equaled 49.3 and the AICc -356.2. Regarding fungal communities (Table 1), the variable selection procedure applied in the multi-linear and clustering models was almost as efficient as for bacteria communities, with 18 and 16 selected key families, respectively. While the median $R^2$ of the multi-linear models were around 0.632, those of clustering models reached 0.933, indicating the same accuracy for fungal communities as bacterial communities.
<table><caption>**Table 1.** Statistics of multi-linear and clustering analyses of the relationships between relative abundance of bacterial and fungal families and basal soil respiration, straw mineralization, and the priming effect. Analyses of relative abundance of bacteria, fungi, and bacteria and fungi get together.</caption><thead><tr><th></th><th>Mineralization</th><th>Duration of Incubation</th><th>Statistical Model</th><th>Initial Number of Families</th><th>Number of Key Families</th><th>Degrees of Freedom</th><th>R<sup>2</sup></th><th>F-Ratio</th><th>AICc</th></tr></thead><tbody><tr><td rowspan="10">bacteria</td><td rowspan="2">soil</td><td rowspan="2">7 days</td><td>lm</td><td>55</td><td>32</td><td>22</td><td>0.993</td><td>112.1</td><td>-319.1</td></tr><tr><td>fclust</td><td>60</td><td>17</td><td>43</td><td>0.966</td><td>39.3</td><td>-328.2</td></tr><tr><td rowspan="2">42 days</td><td>lm</td><td>55</td><td>28</td><td>26</td><td>0.981</td><td>50.8</td><td>-304.9</td></tr><tr><td>fclust</td><td>60</td><td>12</td><td>48</td><td>0.966</td><td>84.6</td><td>-354.8</td></tr><tr><td rowspan="2">straw</td><td>7 days</td><td>55</td><td>38</td><td>16</td><td>0.985</td><td>30.4</td><td>-91.3</td></tr><tr><td>fclust</td><td>60</td><td>13</td><td>47</td><td>0.963</td><td>59.2</td><td>-246.9</td></tr><tr><td rowspan="2">42 days</td><td>lm</td><td>55</td><td>19</td><td>35</td><td>0.883</td><td>14.6</td><td>-314.0</td></tr><tr><td>fclust</td><td>60</td><td>15</td><td>45</td><td>0.927</td><td>24.9</td><td>-357.5</td></tr><tr><td rowspan="2">priming effect</td><td>7 days</td><td>55</td><td>22</td><td>32</td><td>0.879</td><td>11.2</td><td>-351.6</td></tr><tr><td>fclust</td><td>60</td><td>20</td><td>40</td><td>0.985</td><td>99.0</td><td>-480.5</td></tr><tr><td rowspan="10">fungi</td><td rowspan="2">soil</td><td rowspan="2">7 days</td><td>lm</td><td>55</td><td>11</td><td>43</td><td>0.598</td><td>6.1</td><td>-209.3</td></tr><tr><td>fclust</td><td>60</td><td>26</td><td>34</td><td>0.977</td><td>30.7</td><td>-302.4</td></tr><tr><td rowspan="2">42 days</td><td>lm</td><td>55</td><td>23</td><td>31</td><td>0.666</td><td>4.3</td><td>-193.4</td></tr><tr><td>fclust</td><td>60</td><td>14</td><td>46</td><td>0.929</td><td>18.2</td><td>-307.1</td></tr><tr><td rowspan="2">straw</td><td>7 days</td><td>55</td><td>13</td><td>41</td><td>0.585</td><td>4.7</td><td>-109.6</td></tr><tr><td>fclust</td><td>60</td><td>24</td><td>36</td><td>0.982</td><td>40.4</td><td>-236.9</td></tr><tr><td rowspan="2">42 days</td><td>lm</td><td>55</td><td>7</td><td>47</td><td>0.557</td><td>8.8</td><td>-280.6</td></tr><tr><td>fclust</td><td>60</td><td>10</td><td>50</td><td>0.540</td><td>4.2</td><td>-269.9</td></tr><tr><td rowspan="2">priming effect</td><td>7 days</td><td>55</td><td>24</td><td>30</td><td>0.724</td><td>3.5</td><td>-293.1</td></tr><tr><td>fclust</td><td>60</td><td>17</td><td>43</td><td>0.937</td><td>23.3</td><td>-412.7</td></tr><tr><td rowspan="10">bacteria and fungi</td><td rowspan="2">soil</td><td rowspan="2">7 days</td><td>lm</td><td>55</td><td>23</td><td>31</td><td>0.848</td><td>8.0</td><td>-304.2</td></tr><tr><td>fclust</td><td>60</td><td>8</td><td>52</td><td>0.838</td><td>20.5</td><td>-361.0</td></tr><tr><td rowspan="2">42 days</td><td>lm</td><td>55</td><td>120</td><td>13</td><td>107</td><td>0.949</td><td>181.4</td><td>-320.0</td></tr><tr><td>fclust</td><td>60</td><td>120</td><td>14</td><td>106</td><td>0.971</td><td>59.3</td><td>-358.2</td></tr><tr><td rowspan="2">straw</td><td>7 days</td><td>55</td><td>120</td><td>14</td><td>106</td><td>0.938</td><td>37.1</td><td>-214.7</td></tr><tr><td>fclust</td><td>60</td><td>120</td><td>20</td><td>100</td><td>0.979</td><td>29.2</td><td>-407.4</td></tr><tr><td rowspan="2">priming effect</td><td>7 days</td><td>55</td><td>120</td><td>19</td><td>101</td><td>0.959</td><td>49.3</td><td>-428.1</td></tr><tr><td>fclust</td><td>60</td><td>120</td><td>13</td><td>107</td><td>0.933</td><td>145.1</td><td>-395.7</td></tr></tbody></table>

Coefficient of determination (R²), Fisher ratio (F-ratio), and Akaike Information Criterion corrected for small samples (AICc) of the models.

The variable selection procedure retained, in median, only 14 key families out of 120, for the six models (Table 1). The number of degrees of freedom was thus 106, the coefficient of determination R² 0.954, and the F-ratio 54.3. All F-ratio were associated with a probability lower than 10⁻¹⁰. The resulting AICc was -377.0 in median, i.e., the best mean AICc among all models. Applied to soil bacteria and fungus communities combined together, the clustering analyses provided parsimonious and accurate models.

### 3.3. Identifying Bacterial and Fungal Families That Regulate SOM Decomposition Activity

Both multi-linear and clustering approaches enable the identification of bacterial or fungal families as the key actors in basal soil respiration (Figures S3–S6), straw mineralization (Figures S7–S10), and the priming effect (Figures S11–S14) after seven and 42 days of incubation, by separately analyzing each of the two communities. The results are summarized in Table 2.
For each SOM decomposition function and time of incubation, families identified as key actors in initial soil communities by either the multi-linear or the clustering approach accounted for 47% and 42% of the total number of bacterial and fungal families, respectively. Only 13% and 7% of bacterial and fungal families were commonly identified as key for using both approaches, while 37% and 51% were commonly rejected as non-significant (Table 2). This indicates that multi-linear and clustering analyses modeled different aspects of the relationship between microbial families of initial soil communities and SOM decomposition, and this difference was stronger for fungi than for bacteria.

The clustering analyses of bacterial and fungal families gathered together in the same dataset confirmed the results (Figures S15–S20, summarized in Table 3). The variable selection procedure retained less than 14% of total families for a given SOM decomposition function and a given incubation duration, and also rejected 44% of total families showing no relationship for any SOM decomposition function at any incubation duration. This indicates that a small number of microbial families of initial soil communities explained a large part of SOM decomposition after seven and 42 days of incubation.

Although each family was observed in about 30% of soils, the extent of the relative abundances of bacterial and fungal families were very broad, ranging from 1 to 1486, and from 5 to 1048 out of 3000 sequences, respectively (Tables 2 and 3, Figure S1). However, a comparison of medians confirmed that key families were not related to their relative abundance, whatever the SOM decomposition function, the time of duration, or the modelling approach (Table S3).

### 3.4. *Identity of Potentially Key Bacterial and Fungal Controlling the Priming Effect*

The main results of the clustering analyses of bacterial and fungal families gathered together in the same dataset are reported in Figures 1–3. Figure 1 shows the clustering of relative abundance classes of the microbial families of initial soil communities identified as key for basal soil respiration (Figure 1a,b, see Figures S15a and S16a) and straw mineralization (Figure 1c,d, see Figures S17a and S18a) after seven and 42 days of incubation. Figures 2 and 3 report the clustering analyses of key bacterial and fungal families for the priming effect after seven and 42 days of incubation, respectively (see Figures S19 and S20).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Clustering analysis of basal soil respiration (a,b) and straw mineralization (c,d) after seven (a,c) and 42 days of incubation, explained by the relative abundance classes of bacterial and fungal families in soil. The functional groups and the classes of family relative abundance inside each functional group are sorted from left to right by their decreasing effects on the soil property. Each functional group has a color attributed: red for A, dark blue for B, yellow for C, light blue for D, pink for E, and green for F.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Clustering analysis of the priming effect after seven days of incubation, explained by the relative abundance classes of bacterial and fungal families in soil. (a) Clustering tree of classes of family relative abundance. The functional groups and the classes of family relative abundance inside each functional group are sorted from left to right by their decreasing effects on the priming effect. Each functional group has a color attributed: red for A, dark blue for B, yellow for C, light blue for D, pink for E. The above values indicate the significant effects (at *p* < 0.001) of the different functional groups and their interactions. (b) Mean observed priming effect of soils containing a given class of family relative abundance. The linked points correspond to the different classes of increasing relative abundance (from left to right) for each family. The used symbols are the same as those used in (a). (c) Boxplots of the observed priming effect sorted by assembly motifs, i.e., combinations of functional groups. The used symbols are the same as those used in (b). (d) Modeled versus observed priming effect. Different symbols correspond to different assembly motifs, i.e., combinations of functional groups. The red line is the bisector. (b–d) The dotted line is the mean observed priming effect.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Clustering analysis of the priming effect after 42 days of incubation, explained by the relative abundance classes of bacterial and fungal families in soil. (a) Clustering tree of classes of family relative abundance. The functional groups and the classes of family relative abundance inside each functional group are sorted from left to right by their decreasing effects on the priming effect. The above values indicate the significant effects (at *p* < 0.001) of the different functional groups and their interactions. (b) Mean observed priming effect of soils containing a given class of family relative abundance. The linked points correspond to the different classes of increasing relative abundance (from left to right) for each family. The used symbols are the same as those used in (a). (c) Boxplots of the observed priming effect sorted by assembly motifs, i.e., combinations of functional groups. The used symbols are the same as those used in (b). (d) Modeled versus observed priming effect. Different symbols correspond to different assembly motifs, i.e., combinations of functional groups. The red line is the bisector. (b–d) The dotted line is the mean observed priming effect.</figcaption></figure>
Ten bacterial and nine fungal families of initial soil communities were identified as key for the priming effect after seven days of incubation. The 19 microbial families were clustered into five functional groups (Figure 2). Functional groups were labeled from A to E by their decreasing contribution to PE (Figure 2a). Within each group, the relative abundance classes of microbial families are sorted from left to right, by their decreasing contribution to PE (Figure 2a). The specific contribution of the relative abundance classes of family ranged from −15% to +20% of the median PE (Figure 2b). Some families showed a linear contribution, such as the bacteria *Caulobacteraceae*, while others displayed a non-linear contribution, such as the fungus *Sarcosomataceae* and the bacteria *Bacillaceae* (Figure 2b). The first main functional group, noted A, included three bacterial families, *Caulobacteraceae* in its highest relative abundance (class 3/3), *Methylocystaceae* in its lowest relative abundance (class 1/2), and *Bacillaceae* in its medium relative abundance (class 2/3), and three fungal families, *Sarcosomataceae* (class 2/3), *Pluteaceae* (class 2/3), and *Teratosphaeriaceae* (class 3/3) (Figure 2a). The functional group B contained two bacterial families, *Phyllobacteriaceae* (class 1/1) and *Acidobacteria*-Gp5 (class 2/3), and the fungal family *Ambisporaceae* (class 2/3) (Figure 2a). A variance analysis indicates that the functional groups A and B are associated with a PE that is 39% and 17% higher on average, respectively, than the median priming effect (P(A) < $10^{-10}$; P(B) = 0.63 × $10^{-9}$). Moreover, their statistical interaction is associated with an even higher PE of 23% (P(A:B) = 0.51 × $10^{-6}$). As a corollary, all assembly motifs containing the functional group A were associated with high PE, and the assembly motifs containing both the functional groups A and B were associated with the highest PE (Figure 2c). At the opposite, the functional groups C and D are associated with a PE of an average 14% and 12%, respectively (P(C) = 0.37 × $10^{-6}$; P(D) = 0.53 × $10^{-6}$), lower than the median priming effect. However, the functional group C interacts positively with the A and B groups and is associated with a PE of 20% and 3%, respectively (P(A:C) = 0.44 × $10^{-7}$; P(B:C) = 0.62 × $10^{-3}$), higher than the median. The functional group C clustered three bacterial families, *Planococcaceae* (class 1/2), *Acidobacteria*-Gp3 (class 2/2), and *Xanthobacteraceae* (class 2/2), and two fungal families, *Boletaceae* (class 1/3) and *Tubeufiaceae* (class 2/3). The functional group D included two bacterial families, *Acidobacteria*-Gp2 (class 1/3) and *Solirubobacteraceae* (class 2/3), and the fungal family *Cordycipitaceae* (class 1/3): it did not interact with other functional groups. The last functional group E gathered all other bacterial and fungal families identified as key: it is not associated with any significant variation in the priming effect. Finally, the clustering of 19 families of initial soil communities identified as key accurately explained the priming effect after seven days of incubation (R² = 0.959; AICc = −428.1; *p* < $10^{-10}$) (Figure 2d). The clustering was related neither to precipitation (*p* = 0.503) nor to temperature (*p* = 0.027) at sampling sites.

The priming effect intensity was lower after 42 days than seven days of incubation. Seven bacterial and six fungal families of initial soil communities were identified as key for the PE. The 13 microbial families were clustered into five functional groups (Figure 3). The specific effect of relative abundance classes of family ranged from −43% to +60% of the median PE (Figure 3b). The main first functional group A included the group of unknown bacteria (class 3/3), two bacterial families, *Polyangiaceae* (class 2/3) and *Actinospicaceae* (class 2/3), and three fungal families, *Ajellomycetaceae* (class 3/3), *Phaeosphaeriaceae* (class 2/3), and *Montagnulaceae* (class 1/3): the group A is associated with a PE 140% higher than the median (P(A) < $10^{-10}$) (Figure 3a). This means that the co-occurrence of at least two families belonging to the group A doubled the PE variation associated with only a single family (Figure 3b); all assembly motifs containing the functional group A were therefore associated with high PE (Figure 3c). The functional group B is also associated with a PE 51% higher than the median (P(B) = 0.85 × $10^{-9}$); yet it contained only one fungal family, *Chaetomiaceae*, in its lowest relative abundance (Figure 3a). The functional group C clustered the bacterial families *Burkholderiaceae* (class 3/3), *Phyllobacteriaceae* (class 1/1), *Sinobacteraceae* (class 2/2), and *Rhodospirillaceae* (class 2/3), and the fungal family *Montagnulaceae* (class 3/3) (Figure 3a). This group is associated with a PE 86% lower than
the median (P(C) < 10⁻¹⁰). It strongly interacted with the functional group B, associated with a PE 48% even lower than the median (P(B:C) = 0.42 × 10⁻³). As a consequence, the assembly motifs containing the functional group C, and even more when associated with B, showed the lowest PE (Figure 3c). In addition, functional group D is associated with no specific variation of PE, but it strongly interacted with A, which is associated with a PE 63% higher than the median (P(A:D) = 0.34 × 10⁻³): it included the fungal family *Ajellomycetaceae* (class 2/3) and the bacterial family *Burkholderiaceae* (class 1/3) (Figure 3a). Finally, the clustering of only 13 families identified as key in initial soil communities accurately and parsimoniously modeled the PE after 42 days of incubation (R² = 0.933; AICc = -395.7; *p* < 10⁻¹⁰) (Figure 3d). The clustering was related neither to precipitation (*p* = 0.148) nor to temperature (*p* = 0.003) at sample sites.

Except for the bacterial family *Phyllobacteriaceae* and the fungal family *Cordycipitaceae*, all families identified as key in initial soil communities for the priming effect were different after seven and 42 days of incubation (Figures 2b and 3b).

# 4. Discussion

## 4.1. A Few Microbial Families Were Associated with Significant Variations in SOM Decomposition Functions

A striking result is that whatever the soil function, a small number of microbial families of initial soil communities explain most of the SOM decomposition functions after seven and 42 days of incubation. When analyzing bacterial and fungal families gathered together, the clustering analysis identified less than sixteen families as key out of 120, i.e., less than 14% of observed families. These findings confirm that, for a given ecosystem property, biodiversity is most often a saturating function; only a fraction of the identified species contribute significantly to the emergence of the property studied [20–22]. They highlight the importance of procedures of variable selection for analyzing the relationships between biodiversity and ecosystem functioning; these procedures reduce the communities to only the species that are key, at one time or another, for the property studied [23,24]. Note that our analyses are based on microbial composition of initial soil communities before incubation; the microbial communities obviously changed without our knowledge during incubation, mainly by changing their relative abundances. The initial composition of soil microbial communities is therefore key in anticipating SOM decomposition functions after seven and 42 days of incubation, especially the priming effect.

Another striking result is that when bacterial and fungal families were analyzed together, all functional groups contained both bacteria and fungi. This means that some bacteria and fungi are functionally redundant. This functional redundancy suggests that, despite their morphological and physiological differences, bacteria and fungi may be associated with similar processes. However, it is possible that this functional redundancy is only a façade. Like any statistical analysis, a clustering analysis is not causal; it only relates a variable to be explained with explanatory variables. The clustering of families or classes of families in relative abundance in the same functional group therefore does not prejudice the role played by each of the clustered families. It only indicates that different combinations of microbial families are associated with different variations in SOM decomposition activity which, moreover, take place seven and 42 days after analysis of the initial composition of soil microbial communities. The functional groups generated by clustering analysis can therefore be composed of microbial families that directly and actively affect SOM decomposition activity after seven or 42 days of incubation, and families that contribute to variations in biotic or abiotic conditions are associated with variations in SOM decomposition.

For instance, the bacterial *Methylocystaceae* family is identified as key in the initial soil community for basal soil respiration, straw mineralization, and the priming effect after seven days of incubation. As the *Methylocystaceae* family is a strictly methanotrophic family that releases CO₂ by oxidizing CH₄ [25], they could not have contributed to any heterotrophic respiration. Therefore, our results highlight that soil conditions favorable
to the *Methylocystaceae* family promote the priming effect after seven days of incubation. As another example, the bacterial *Polyangiaceae* family is identified as key in the initial soil community for basal soil respiration and the priming effect after 42 days of incubation; it then belongs to the functional group A, which is associated with a high priming effect. *Polyangiaceae* are bacteria that prey on other bacteria [26]; this family possibly plays a positive role in regulating specific populations having a negative effect on PE.

Our results also show that key families are neither the most frequent nor the most
abundant families. For instance, when the functional groups A and B co-occurred in the
initial soil community, the highest PE was observed after seven days of incubation. All
bacterial families identified, *Caulobacteraceae*, *Methylocystaceae*, and *Bacillaceae*, in functional
group A, and *Phyllobacteriaceae* and *Acidobacteria*-Gp5 in functional group B, are not highly
abundant. Likewise, the fungal families *Sarcosomataceae* and *Pluteaceae* in functional group
A, and *Ambisporaceae* in functional group B, are not very abundant. Only the ubiquist
fungus *Teratosphaeriaceae* was moderately abundant and co-occurred in its highest relative
abundance class in its initial soil community. All the microbial families were, however,
associated with high PE. When the functional groups A and B co-occurred in the initial
soil community, the highest PE was observed after 42 days of incubation. The unknown
bacteria were very abundant and co-occurred in their most abundant relative class, although
they correspond to a phylogenetically and functionally heterogeneous group. The fungus
families *Phaeosphaeriaceae* and *Ajellomycetaceae*, however, were much less abundant in soils
and co-occurred in the medium relative abundance class. It is possible that the most
abundant families play an active role in SOM mineralization functions after seven or
42 days of incubation, while the less abundant families were rather associated with biotic
or abiotic conditions favoring SOM decomposition.

We should, however, remember that (i) first, next generation sequencing technology
is not a quantitative molecular method [14]; (ii) second, the occurrence of each microbial
family was determined by only a few sequences, although it was observed in about 30%
of soils. This lack of relationship between observed effects and frequency or abundance
of microbial families should nevertheless be a question to us. In ecology, it is common to
eliminate the least frequent or least abundant families because these rare species are a priori
assumed to be inefficient [27,28]. This practice tends to underestimate the contribution
of rare families to community functioning, even though some of them may play a key
role [29–31]. In bacterial communities, the clustering analysis confirmed that most of the
families were eliminated in the multi-linear analysis because of their low frequency or
abundance not being key. In fungal communities, clustering analysis identified more than
half of the families as key that were eliminated in the multi-linear analysis due to their low
frequency or abundance. Our results thus highlight that taxonomic rarity is not a relevant
criterion. Clustering analysis enables us to more deeply explore the relationship between
biodiversity and ecosystem functioning than multi-linear analysis. It should be favored
in looking for species that play a key role in the relationships between biodiversity and
ecosystem functioning.

## 4.2. Multi-Linear and Clustering Analyses Enable the Identification of Microorganisms That Stimulate SOM Decomposition Activities after Seven and 42 Days of Incubation

Our results show that SOM decomposition activities are accurately explained by variations in the relative abundance of bacterial families in pristine soil. It is more accurately and parsimoniously explained by the co-occurrence of bacterial families at different levels of relative abundance. However, the improvement is marginal. Multi-linear and combinatorial approaches are therefore equally relevant. These results support the observations of Razanamalala et al. [13] made at the phylum level on the same samples, who found that bacteria were strongly linearly correlated with PE. The effect induced by bacteria on SOM decomposition activity would therefore be close to a composition effect.

Contrastingly, the relationship between fungi occurrence and any SOM decomposition activity at any time duration are hardly explained by multi-linear analyses, in contrast to the
clustering analysis, which remains accurate. The relationships between fungi occurrence and SOM decomposition activities are thus non-linear; their occurrence moderately explains the PE after seven days of incubation, but much more explains the PE after 42 days. At this time, the occurrence of some functional groups composed of both bacteria and fungi in initial soil communities is related to a double PE, while the occurrence of other functional groups is related to significantly reduced PE. The main bacterial and fungal families selected for their potential involvement in PE recorded after 42 days of incubation are not the same as at the beginning of incubation, and non-linear processes of high intensity have taken over linear processes.

### 4.3. Linear and Clustering Analyses Enlighten the Priming Effect

Several authors have suggested that mechanisms underlying the priming effect can occur simultaneously: one may dominate for a period of time and then be replaced by another depending on changes in environmental conditions, including the availability of mineral nitrogen in soil solution [2,32,33]. Based on the literature, Razanamalala et al. [13] proposed that the PE measured after seven days was mainly generated by stoichiometric decomposition, and that the PE measured after 42 days was mainly the result of nutrient mining. Strikingly, we showed that microbial families favoring PE measured after seven and 42 days of incubation were not the same.

After seven days of incubation, the functional group A associated with the highest priming effect seems functionally convergent as it is mainly composed of bacterial families endowed with specific enzymatic capacities for the decomposition of lignin (Caulobacteraceae, [34]) and cellulose (Bacillaceae, [35]). This group also contains saprotrophic fungi well known for their lignocellulolytic capacities [36]. These results support the stoichiometric decomposition hypothesis after seven days of incubation, yet none of these families appeared to have a preponderant role, either in basal soil respiration or in straw mineralization. Although they might play a key role in PE, they are not among the most abundant nor the most frequent families. These families could therefore benefit from the exoenzymatic activity of other populations directly involved in straw mineralization as suggested in the stoichiometric hypothesis. Our results also show that PE depended on positive and negative interactions between families of group A and other functional groups. For instance, group B, which increased the effect of group A, is composed of two bacterial families' characteristics of rich media, the Phyllobacteriaceae [37], or at least enriched in nutrients, the Acidobacteria-GP5 [38], and whose common feature is their ability to denitrify. Group C has a rather negative effect on PE as well as on group A, and is functionally more diversified than groups A and B.

When the labile part of FOM and the available nitrogen is exhausted, SOM-decomposers can use the energy of FOM catabolites to mine nitrogen in SOM [5]. Fontaine et al. [36] suggested that fungi are the main decomposers of cellulose and are thus the main actors of the nutrient mining PE. This important role in nutrient mining is also linked to the ability of fungi to explore a large volume of soil, having their mycelium on FOM and SOM simultaneously. It is indeed reasonable to assume that interactions between bacteria are limited due to their relative immobility. In contrast, fungi develop an extensive mycelial network that allows them to interact over long distances with other populations of bacteria and fungi. After 42 days of incubation, unlike after seven days of incubation, the functional group A that increased PE by 140% is functionally quite heterogeneous. It is composed of two families of saprotrophic Ascomycetes (Phaeosphaeriaceae and Ajellomycetaceae, [39,40]), a family of acidophilic Actinomycetes (Actinospicaceae, [41]) known for their ability to produce secondary metabolites [42], and unknown bacteria, which are functionally heterogenous but generally considered oligotrophic because they do not grow in the too-rich conditions of laboratory environments. Against all expectations, it is the bacteria (unknown and Actinospicaceae) that seem to engage in co-metabolism between FOM and SOM, not the fungi. However, as the group of unknown bacteria brings together various families under the same name, it can actually be divided into families specialized in either the breakdown
of FOM or SOM, without co-metabolism. Enrichment of the databases should gradually reduce the size of this group in future studies. It has recently been suggested that PE could also be generated by the release of molecules with a high affinity for soil minerals, making it possible to desorb the initially associated labile organic matter rich in N and P [5]. Representatives of *Burkholderiaceae* are able to weather minerals by the release of gluconic acid [43]. However, in our study, *Burkholderiaceae* have a strong inhibitory role on PE after 42 days of incubation, linked to a density effect (their lowest and highest abundance classes are found in functional groups C and D). At the same time, they have an inhibiting role in straw mineralization. *Burkholderiaceae* have other capacities, such as exopolysaccharide production, antifungal release, and organic and inorganic *p* solubilization via siderophores production [44,45]. All those capacities can interfere directly or indirectly in the PE generation process. Therefore, the role of this family in the late priming effect has to be more deeply investigated in the future.

**5. Conclusions**

The clustering approach makes it possible to highlight taxa with non-linear effects on a flow resulting from complex processes, such as the priming effect, but also the effect on this flow of interactions between certain taxa. This approach has the advantage of not being constrained by the ratio between the number of variables and the number of samples, unlike linear statistical approaches. Normally applied to occurrence matrices, it has been adapted here to an abundance matrix by dividing each taxon into one, two or three abundance classes according to their respective densities. The results obtained are surprising and open up new challenges for their interpretation. This approach has proven its interest and deserves to be applied to more recent datasets resulting from sequencing techniques and more advanced bioinformatics analysis, which could only improve the quality and precision of the results. It is obvious that the biofunctioning of the soil is subservient to interactions between various functional entities and that this type of method considerably increases the field of possibilities in terms of understanding its various workings.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11051106/s1, Figure S1: Relative abundance of bacterial and fungal families; Figure S2: Conversion of abundances into classes of relative abundance for clustering analysis; Figure S3: Clustering analysis of basal soil respiration after seven days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S4: Clustering analysis of basal soil respiration after 42 days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S5: Clustering analysis of basal soil respiration after seven days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S6: Clustering analysis of basal soil respiration after 42 days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S7: Clustering analysis of straw mineralization after seven days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S8: Clustering analysis of straw mineralization after 42 days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S9: Clustering analysis of straw mineralization after seven days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S10: Clustering analysis of straw mineralization after 42 days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S11: Clustering analysis of the priming effect after seven days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S12: Clustering analysis of the priming effect after 42 days of incubation, explained by the relative abundance classes of bacterial families in soil; Figure S13: Clustering analysis of the priming effect after seven days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S14: Clustering analysis of the priming effect after 42 days of incubation, explained by the relative abundance classes of fungal families in soil; Figure S15: Clustering analysis of basal soil respiration after seven days of incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Figure S16: Clustering analysis of basal soil respiration after 42 days of incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Figure S17: Clustering analysis of straw mineralization after seven days of
incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Figure S18: Clustering analysis of straw mineralization after 42 days of incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Figure S19: Clustering analysis of the priming effect after seven days of incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Figure S20: Clustering analysis of the priming effect after 42 days of incubation, explained by the relative abundance classes of bacterial and fungal families together in soil; Table S1: Physicochemical description of soil samples; Table S2: Relationship between families, codenames, and phyla of soil bacteria and fungi; Table S3: Comparison of abundance medians between key and non-key families using non-parametric Kruskal-Wallis test.

**Author Contributions:** Data acquisition, L.B. and K.R.; statistical analyses, software, writing original draft preparation, B.J.; Validation, writing—review and editing, B.J., K.R., C.V. and L.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** The dataset treated in the present study was funded by the French Foundation for Research on Biodiversity (FRB-AAP-SCEN-2013 II—CAMMiSolE project).

**Data Availability Statement:** DNA sequences were deposited in the European Nucleotide Archive, under the study accession number PRJEB19651. The package functClust is available on R-CRAN at https://CRAN.R-project.org/package=functClust (accessed on 20 March 2023). The datasets and R-codes used will be available on the INRAE dataverse at https://doi.org/10.15454/LRY3M2 (accessed on 20 March 2023).

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Kuzyakov, Y.; Friedel, J.; Stahr, K. Review of mechanisms and quantification of priming effects. *Soil Biol. Biochem.* **2000**, *32*, 1485–1498. [CrossRef]

2. Fontaine, S.; Mariotti, A.; Abbadie, L. The priming effect of organic matter: A question of microbial competition? *Soil Biol. Biochem.* **2003**, *35*, 837–843. [CrossRef]

3. Chen, R.; Senbayram, M.; Blagodatsky, S.; Myachina, O.; Dittert, K.; Lin, X.; Blagodatskaya, E.; Kuzyakov, Y. Soil C and N availability determine the priming effect: Microbial N mining and stoichiometric decomposition theories. *Glob. Chang. Biol.* **2014**, *20*, 2356–2367. [CrossRef] [PubMed]

4. Keiluweit, M.; Bougoure, J.J.; Nico, P.S.; Pett-Ridge, J.; Weber, P.K.; Kleber, M. Mineral protection of soil carbon counteracted by root exudates. *Nat. Clim. Chang.* **2015**, *5*, 588–595. [CrossRef]

5. Bernard, L.; Basile-Doelsch, I.; Derrien, D.; Fanin, N.; Fontaine, S.; Guenet, B.; Karimi, B.; Marsden, C.; Maron, P.-A. Advancing the mechanistic understanding of the priming effect on soil organic matter mineralisation. *Funct. Ecol.* **2022**, *36*, 1355–1377. [CrossRef]

6. Bernard, L.; Maron, P.-A.; Mougel, C.; Nowak, V.; Lévêque, J.; Marol, C.; Balesdent, J.; Gibiat, F.; Ranjard, L. Contamination of Soil by Copper Affects the Dynamics, Diversity, and Activity of Soil Bacterial Communities Involved in Wheat Decomposition and Carbon Storage. *Appl. Environ. Microbiol.* **2009**, *75*, 7565–7569. [CrossRef]

7. Friedrich, M.W. Stable-isotope probing of DNA: Insights into the function of uncultivated microorganisms from isotopically labeled metagenomes. *Curr. Opin. Biotechnol.* **2006**, *17*, 59–66. [CrossRef]

8. Wilson, J.B. Shoot competition and root competition. *J. Appl. Ecol.* **1988**, *25*, 279–296. [CrossRef]

9. Loreau, M.; Hector, A. Partitioning selection and complementarity in biodiversity experiments. *Nature* **2001**, *412*, 72–76. [CrossRef]

10. Loreau, M.; Hector, A. Not even wrong: Comment by Loreau and Hector. *Ecology* **2019**, *100*, e02794. [CrossRef]

11. Jaillard, B.; Deleporte, P.; Loreau, M.; Violle, C. A combinatorial analysis using observational data identifies species that govern ecosystem functioning. *PLoS ONE* **2018**, *13*, e0201135. [CrossRef] [PubMed]

12. Jaillard, B.; Deleporte, P.; Isbell, F.; Loreau, M.; Violle, C. Consistent functional clusters explain the effects of biodiversity on ecosystem productivity in a long-term experiment. *Ecology* **2021**, *102*, e03441. [CrossRef] [PubMed]

13. Razanamalala, K.; Razafimbelo, T.; Maron, P.-A.; Ranjard, L.; Chemidlin, N.; Lelièvre, M.; Dequiedt, S.; Ramaroson, V.H.; Marsden, C.; Becquer, T.; et al. Soil microbial diversity drives the priming effect along climate gradients: A case study in Madagascar. *ISME J.* **2018**, *12*, 451–462. [CrossRef]

14. Maron, P.-A.; Mougel, C.; Ranjard, L. Soil microbial diversity: Methodological strategy, spatial over-view and functional interest. *Comptes Rendus Biol.* **2011**, *334*, 403–411. [CrossRef] [PubMed]

15. Terrat, S.; Christen, R.; Dequiedt, S.; Lelièvre, M.; Nowak, V.; Regnier, T.; Bachar, D.; Plassart, P.; Wincker, P.; Jolivet, C.; et al. Molecular biomass and MetaTaxogenomic assessment of soil microbial communities as influenced by soil DNA extraction procedure. *Microb. Biotechnol.* **2012**, *5*, 135–141. [CrossRef]
16. Tardy, V.; Spor, A.; Mathieu, O.; Lévèque, J.; Terrat, S.; Plassart, P.; Regnier, T.; Bardgett, R.D.; van Der Putten, W.H.; Roggero, P.P.; et al. Shifts in microbial diversity through land use intensity as drivers of carbon mineralization in soil. *Soil Biol. Biochem.* **2015**, *90*, 204–213. [CrossRef]

17. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. *Nucleic Acids Res.* **2013**, *41*, 590–596. [CrossRef]

18. Burnham, K.P.; Anderson, D.R. Multimodel Inference—Understanding AIC and BIC in Model Selection. *Sociol. Methods Res.* **2004**, *33*, 261–304. [CrossRef]

19. R Development Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2009; ISBN 3-900051-07-0. Available online: https://www.R-project.org (accessed on 20 March 2023).

20. Isbell, F.; Gonzalez, A.; Loreau, M.; Cowles, J.; Díaz, S.; Hector, A.; Mace, G.M.; Wardle, D.A.; O'Connor, M.I.; Duffy, J.E.; et al. Linking the influence and dependence of people on biodiversity across scales. *Nature* **2017**, *546*, 65–72. [CrossRef]

21. Yu, X.; Polz, M.F.; Alm, E.J. Interactions in self-assembled microbial communities saturate with diversity. *ISME J.* **2019**, *13*, 1602–1617. [CrossRef]

22. Clements, C.S.; Hay, M.E. Biodiversity has a positive but saturating effect on imperiled coral reefs. *Sci. Adv.* **2021**, *7*, eabi8592. [CrossRef]

23. Maire, E.; Villéger, S.; Graham, N.A.; Hoey, A.S.; Cinner, J.; Ferse, S.C.; Aliaume, C.; Booth, D.J.; Feary, D.A.; Kulbicki, M.; et al. Community-wide scan identifies fish species associated with coral reef services across the Indo-Pacific. *Proc. R. Soc. B Biol. Sci.* **2018**, *285*, 20181167. [CrossRef] [PubMed]

24. Delalandre, L.; Gaüzère, P.; Thuiller, W.; Cadotte, M.; Mouquet, N.; Mouillot, D.; Munoz, F.; Denelle, P.; Loiseau, N.; Morin, X.; et al. Functionally distinct tree species support long-term productivity in extreme environments. *Proc. R. Soc. B* **2022**, *289*, 20211694. [CrossRef] [PubMed]

25. Chowdhury, T.R.; Dick, R.P. Ecology of aerobic methanotrophs in controlling methane fluxes from wetlands. *Appl. Soil Ecol.* **2013**, *65*, 8–22. [CrossRef]

26. Petters, S.; Groß, V.; Söllinger, A.; Pichler, M.; Reinhard, A.; Bengtsson, M.M.; Urich, T. The soil microbial food web re-visited: Predatory myxobacteria as keystone taxa? *ISME J.* **2021**, *15*, 2665–2675. [CrossRef] [PubMed]

27. Gaston, K.J.; Lawton, J.H. The population ecology of rare species. *J. Fish Biol.* **1990**, *37*, 97–104. [CrossRef]

28. Marchant, R. How important are rare species in aquatic community ecology and bioassessment? A comment on the conclusions of Cao et al. *Limnol. Oceanogr.* **1999**, *44*, 1840–1841. Available online: http://pascal-francis.inist.fr/vibad/index.php?action=getRecordDetail&idt=1209822 (accessed on 20 March 2023).

29. Cao, Y.; Larsen, D.P. Rare species in multivariate analysis for bioassessment: Some considerations. *J. N. Am. Benthol. Soc.* **2001**, *20*, 144–153. [CrossRef]

30. Cao, Y.; Williams, D.D.; Williams, N.E. How important are rare species in aquatic community ecology and bioassessment? *Limnol. Oceanogr.* **1998**, *43*, 1403–1409. [CrossRef]

31. Violle, C.; Thuiller, W.; Mouquet, N.; Munoz, F.; Kraft, N.J.; Cadotte, M.W.; Livingstone, S.W.; Mouillot, D. Functional Rarity: The Ecology of Outliers. *Trends Ecol. Evol.* **2017**, *32*, 356–367. [CrossRef]

32. Derrien, D.; Plain, C.; Courty, P.E.; Gelhaye, L.; Moerdijk-Poortvliet, T.C.; Thomas, F.; Versini, A.; Zeller, B.; Koutika, L.S.; Boschker, H.T.; et al. Does the addition of labile substrate destabilise old soil organic matter? *Soil Biol. Biochem.* **2014**, *76*, 149–160. Available online: https://www.sciencedirect.com/science/article/pii/S0038071714001606 (accessed on 20 March 2023). [CrossRef]

33. Blagodatskaya, E.; Khomyakov, N.; Myachina, O.; Bogomolova, I.; Blagodatsky, S.; Kuzyakov, Y. Microbial interactions affect sources of priming induced by cellulose. *Soil Biol. Biochem.* **2014**, *74*, 39–49. [CrossRef]

34. Wilhelm, R.C.; Singh, R.; Eltis, L.D.; Mohn, W.W. Bacterial contributions to delignification and lignocellulose degradation in forest soils with metagenomic and quantitative stable isotope probing. *ISME J.* **2019**, *13*, 413–429. [CrossRef]

35. Mandic-Mulec, I.; Stefanic, P.; van Elsas, J. Ecology of Bacillaceae. *Microbiol. Spectr.* **2015**, *3*, TBS-0017–2013. [CrossRef]

36. Fontaine, S.; Hénault, C.; Aamor, A.; Bdioui, N.; Bloor, J.M.G.; Maire, V.; Mary, B.; Revaillot, S.; Maron, P.A. Fungi mediate long term sequestration of carbon and nitrogen in soil through their priming effect. *Soil Biol. Biochem.* **2011**, *43*, 86–96. [CrossRef]

37. Xing, X.; Xu, H.; Zhang, W.; Hou, H.; Qin, H.; Liu, Y.; Zhang, L.; Fang, Y.; Wei, W.; Sheng, R. The characteristics of the community structure of typical nitrous oxide-reducing denitrifiers in agricultural soils derived from different parent materials. *Appl. Soil Ecol.* **2019**, *142*, 8–17. [CrossRef]

38. Diamond, S.; Andeer, P.F.; Li, Z.; Crits-Christoph, A.; Burstein, D.; Anantharaman, K.; Lane, K.R.; Thomas, B.C.; Pan, C.; Northen, T.R.; et al. Mediterranean grassland soil C–N compound turnover is dependent on rainfall and depth, and is mediated by genomically divergent microorganisms. *Nat. Microbiol.* **2019**, *4*, 1356–1367. [CrossRef]

39. Karunarathna, A.; Papizadeh, M.; Senanayake, I.C.; Jeewon, R.; Phookamsak, R.; Goonasekara, I.D. Novel fungal species of Phaeosphaeriaceae with an asexual/sexual morph connection. *Mycosphere* **2017**, *8*, 1818–1834. [CrossRef]

40. Marin-Felix, Y.; Stchigel, A.M.; Cano-Lira, J.F.; Sanchis, M.; Mayayo, E.; Guarro, J. Emmonsiellopsis, a new genus related to the thermally dimorphic fungi of the family Ajellomycetaceae. *Mycoses* **2015**, *58*, 451–460. [CrossRef]

41. Cavaletti, L.; Monciardini, P.; Schumann, P.; Rohde, M.; Bamonte, R.; Busti, E.; Sosio, M.; Donadio, S. *Actinospica robiniae* gen. nov., sp. nov. and *Actinospica acidiphila* sp. nov.: Proposal for *Actinospicaceae* fam. nov. and *Catenulisporinae* subord. nov. in the order *Actinomycetales*. *Int. J. Syst. Evol. Microbiol.* **2006**, *56*, 1747–1753. [CrossRef]
42. Guo, X.; Liu, N.; Li, X.; Ding, Y.; Shang, F.; Gao, Y.; Ruan, J.; Huang, Y. Red Soils Harbor Diverse Culturable Actinomycetes That Are Promising Sources of Novel Secondary Metabolites. *Appl. Environ. Microbiol.* **2015**, *81*, 3086–3103. [CrossRef] [PubMed]

43. Uroz, S.; Christophe, C.; Turpault, M.-P.; Frey-Klett, P. Mineral weathering by bacteria: Ecology, actors and mechanisms. *Trends Microbiol.* **2009**, *17*, 378–387. [CrossRef]

44. Emmett, B.D.; Lévesque-Tremblay, V.; Harrison, M.J. Conserved and reproducible bacterial communities associate with extraradical hyphae of arbuscular mycorrhizal fungi. *ISME J.* **2021**, *15*, 2276–2288. [CrossRef] [PubMed]

45. Wang, S.; Walker, R.; Schicklberger, M.; Nico, P.S.; Fox, P.M.; Karaoz, U.; Chakraborty, R.; Brodie, E.L. Microbial Phosphorus mobilization strategies across a natural nutrient limitation gradient and evidence for linkage with Iron solubilization traits. *Front. Microbiol.* **2021**, *12*, 572212. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.