# Effects of habitat constraints on soil microbial community function

Naoise Nunan, Julie Leloup, Leo S. Ruamps, Valerie Pouteau, Claire Chenu

► To cite this version:

Naoise Nunan, Julie Leloup, Leo S. Ruamps, Valerie Pouteau, Claire Chenu. Effects of habitat con-straints on soil microbial community function. Scientific Reports, 2017, 7, pp.4280. 10.1038/s41598-017-04485-z . hal-01567140

HAL Id: hal-01567140

https://agroparistech.hal.science/hal-01567140

Submitted on 21 Jul 2017

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>
SCIENTIFIC REPORTS

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

OPEN

# Effects of habitat constraints on soil microbial community function

Naoise Nunan<fnref n="1" />, Julie Leloup<fnref n="1" />, Léo S. Ruamps<fnref n="1" />, Valérie Pouteau² & Claire Chenu²

An underlying assumption of most soil carbon (C) dynamics models is that soil microbial communities are functionally similar; in other words, that microbial activity under given conditions is not dependent on the composition or diversity of the communities. Although a number of studies have indicated that microbial communities are not intrinsically functionally similar, most soil C dynamics models can adequately describe C dynamics without explicitly describing microbial functioning. Here, we provide a mechanistic basis for reconciling this apparent discrepancy. In a reciprocal transplant experiment, we show that the environmental context (soil and pore-network properties) of microbial communities can constrain the activity of functionally different communities to such an extent that their activities are indistinguishable. The data also suggest that when microbial activity is less constrained, the intrinsic functional differences among communities can be expressed. We conclude that soil C dynamics may depend on microbial community structure or diversity in environments where their activity is less constrained, such as the rhizosphere or the litter layer, but not in oligotrophic environments such as the mineral layers of soil.

Received: 9 November 2016

Accepted: 16 May 2017

Published online: 27 June 2017

An underlying assumption of soil C dynamics models is that microbial communities are all similar from a functional viewpoint, regardless of their composition or diversity. Influential soil C dynamics models such as RothC<fnref n="1" /> or Century² for example, do not explicitly describe soil microbial communities. In these models microbial activity is represented by the decay rate constants of the various C pools, but the rate constants are not linked microbial composition, diversity or physiology. The assumption of functional similarity is central to how these models describe the responses of soil C dynamics to perturbations, such as those encountered with global climate change. Indeed, models in which the physiological responses of microbial communities to external perturbations are linked to microbial decomposition of organic C can result in different predictions of C dynamics relative to models that don't and provide a better prediction of the apparent attenuation of the soil respiratory response to warming³, ⁴. Changes in the carbon use efficiency of soil microbial communities can explain this apparent attenuation³ and it has been suggested that the carbon use efficiency of microbial communities is related to the communities' diversity⁵, although this remains to be demonstrated.

However, a substantial body of literature lends credence to the assumption of functional similarity amongst microbial communities: a meta-analysis of the relationship between microbial community structure and ecosystem process rates has suggested that abiotic variables are better predictors of processes such as organic C or N mineralisation than microbial community structure⁶ and little evidence of a relationship between microbial diversity and C dynamics has been uncovered to date in mineral soils, other than when communities have extremely low levels of diversity⁷⁻¹². The lack of relationship between process rates and microbial communities has been explained by a high degree of functional redundancy within microbial communities¹³, due to the very high levels of diversity and the physiological flexibility of microbial communities¹⁴. However, the studies investigating the relationship between microbial diversity and process rates established the different diversity treatments by diversity erosion using either serial dilutions of microbial suspensions¹¹ or by differential fumigation⁷. These approaches reduce diversity by preferentially removing the least abundant species or the species most sensitive to fumigation. Therefore, all the diversity levels may have contained the same active microbial groups, resulting in similar levels of activity. Furthermore, significant relationships between microbial community diversity and process rates have been detected in a number of studies, although these have been in the organic or litter layers¹⁵, ¹⁶, with low levels of diversity¹⁶ or the diversity effect appears to have been confounded with a biomass effect¹⁷. Philippot et al.¹⁸ showed a significant relationship between the diversity of denitrifiers and denitrification rates in

<footnote n="1">¹iEES Paris, UMR 7618 CNRS-UPMC-INRA-IRD-Paris 7-UPEC, 4 place Jussieu, 75005, Paris, France. ²EcoSys, 1402 UMR EGC-ECOSYS, INRA- AgroParisTech-Université Paris Saclay, 78850, Thiverval-Grignon, France. Correspondence and requests for materials should be addressed to N.N. (email: naoise.nunan@upmc.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Water-filled pore volumes in the silty (black) and sandy (grey) soils at saturation and at the two incubation matric potentials (a), and the proportion of water-filled pores in the pore network of the silty (black) and sandy (grey) soils (b). Error bars are standard error of the mean (n = 4).</figcaption></figure>

mineral soils. However, the levels of denitrifier diversity were much lower¹⁸ than the levels of diversity of microbial organisms involved in organic C decomposition.

A common feature of all of these studies and models is that they do not account for the microbial abiotic environment. A number of recent reviews of the literature have concluded that soil C dynamics cannot be fully apprehended without considering the constraints imposed by the abiotic environment on decomposer communities¹⁹, ²⁰. Although it is well established that the mineralisation of soil C is strongly related to the total organic C content of the soil, it is believed that the processes and constraints that regulate this activity occur at fine scales²¹⁻²³. The mineralisation of organic C added to large pores has been found to be greater than in smaller pores²⁴, ²⁵ and Strong et al.²⁶ showed that the decomposition of added plant material was most strongly correlated with pores with a maximum neck diameter of 15–60 µm. However, it is now well established that microbial community structure is related to the micro-environment²⁷, ²⁸ and it has been shown that there is a tight relationship between the structure of microbial communities and the pore network²⁵, ²⁹: the structure and diversity of microbial communities can change as a function of the neck diameter and connectivity of the pores in which they are located. Therefore, the differences in organic C decomposition in different regions of the soil pore network cannot necessarily be attributed to abiotic constraints on microbial activity as these putative constraints are confounded with differences within the microbial communities themselves.

It is necessary to understand the mechanisms that regulate process rates and the relative importance of the different regulatory mechanisms under different conditions. The most powerful experimental approach for distinguishing the effects of microbial community properties from abiotic constraints on activity, as well as their interactions, on process rates is the “reciprocal transplant” approach, where communities from different environments are exposed to each other’s native environments, effectively decoupling the relationship between microbial communities and their environment³⁰. A number of reciprocal transplant studies have been carried out. Strickland et al.¹⁴ found that microbial communities accounted for a significant portion of the variation in litter decomposition and concluded that the implicit assumption of functional equivalence in C dynamics models was incorrect for litter decay. The activity of extracellular enzymes also appeared to be dependent, at least in part, on the type of microbial community³¹. However, in both studies, the properties of the abiotic environment also explained a significant portion of the variation in activity and Griffiths et al.³² found that the resistance and resilience of microbial decomposition to external stresses was more closely related to the soil properties than to the microbial communities.

Here, we present the results of a reciprocal transplant experiment in which sterile samples of soil were inoculated with their native microbial communities or the communities from the other soil at two different matric potentials. The objective was to determine the contribution of intrinsic microbial properties and of abiotic soil properties, at both the pore and soil scales, to the variation of soil organic C mineralisation rates.

## Results

The silty soil retained more water at each matric potential and had a lower total porosity than the sandy soil (Fig. 1). Although the same volume of the pore network was inoculated in each treatment (0.05 cm³ g⁻¹ soil), the proportion of the total pore volume inoculated was lower in the sandy soil (Fig. 1). The microbial inocula
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Mean carbon mineralisation curves of the microbial communities from the sandy soil (left panel) and from the silty soil (right panel) when inoculated into the silty soil (square symbols) and into the sandy soil (round symbols). Solid symbols are soils that received the microbial inoculum in the more-connected micro-environments and open symbols are soils that received the inoculum in the less-connected micro-environments. Lines represent model fit. Bars represent standard-error of the mean, where the SEM is greater than the size of the symbol.</figcaption></figure>

occupied 19 and 16% of the total pore volume in the silty and sandy soils, respectively. The water-filled pore volume represented 46 and 18% of the total pore volume in the silty and sandy soils, respectively, in the less connected treatment. In the more connected treatment, the water-filled pore volume represented 74 and 40% of the total pore volume in the silty and sandy soils, respectively (Fig. 1).

The average cumulative mineralisation curves for each treatment combination are presented in Fig. 2. At the end of the incubation significantly more C had been mineralised in the silty soil ($P < 0.01$; Table S1, Fig. 2 and S2), regardless of the microbial community or of the micro-environment. Microbial communities in the more-connected micro-environments also produced more $CO_2$ than those in less-connected micro-environments ($P < 0.01$; Table S1, Fig. 2 and S2), regardless of soil. No significant difference between the amounts of $CO_2$ produced by the different microbial communities during the incubation was observed (Table S1, Fig. 2 and S2). However, there was a significant microbial inoculum x micro-environment interaction ($P < 0.01$; Table S1). This was because microbial communities from the silty soil in the more-connected micro-environment of the silty soil did not mineralise more C than when in the less-connected micro-environment of the silty soil, contrary to all the other treatment combinations, where more C was mineralised in the more-connected micro-environments (Fig. 2 and S2).

In order to understand the effects of the different treatments in finer detail, the cumulative mineralisation curves were fitted with two-compartment first-order models and the parameters derived from these models were analysed. An ANOVA of the model parameters showed that the rate at which the labile pool of organic C was mineralised (parameter $\alpha$) depended significantly on the micro-environment ($P < 0.01$) and microbial inoculum ($P < 0.01$), but not on the soil (Fig. 3 and Table S2). Although the $\alpha$ values associated with communities from the sandy soil were always higher than those associated with the silty soil communities, there was a small but significant microbial community x micro-environment interaction for the parameter $\alpha$ (Table S2). This was due to the difference between the communities being greater in the less-connected micro-environment, possibly as a result of the high variability in the measurements in one of the more-connected micro-environment treatments (Fig. 3). The size of the labile pool of organic C (parameter a) was significantly affected by soil type ($P < 0.01$), micro-environment ($P < 0.01$) and microbial inoculum ($P < 0.01$) (Fig. 3 and Table S3). The size of the pool of organic C that was mineralised at a slower rate (parameter b) and the rate at which the pool was mineralised (parameter $\beta$) were significantly affected by soil ($P < 0.01$) and micro-environment ($P < 0.01$), but not by microbial inoculum (Fig. 3 and Tables S4 and S5). There were also significant microbial inoculum x micro-environment interactions for parameters $\beta$ and b. These were due to the rate constant $\beta$ of the microbial communities from the sandy soil being higher than that of microbial communities from the silty soil in the less-connected micro-environment of both soils, but not in the more-connected micro-environment. The size of the C pool b associated with the microbial communities from the sandy soil was smaller than that associated with the microbial communities from the silty soil in the less-connected micro-environment of both soils
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Treatment effects on the rate constant (parameter α) of the labile pool of organic C (top left panel), on the rate constant (parameter β) of the slow pool of organic C (top right panel), on the size of the labile pool of organic C (parameter a, bottom left panel) and on the size of the slow pool of organic C (parameter b, bottom right panel). The box-plots on the left of the dashed line are parameters for the silty soil (Silt) and those on the right for the sandy soil (Sand). The box-plots on the white background are the parameters for the less-connected micro-environments (Less) and those on the grey background for the more-connected micro-environments (More). The white box-plots are the parameters for microbial communities from the silty soil (SiltMC) and the dark grey box-plots for the communities from the sandy soil (SandMC).</figcaption></figure>

but not in the more-connected micro-environments. However, the differences observed in the less-connected micro-environment were not significant when analysed on their own (*P* > 0.05).

The time taken for 99% of the labile pool of organic C to be consumed was estimated using the values of the model parameters (Table S6). Less than 1% of the pool *a* remained after 18.5 days in the sandy soil that had been inoculated in the more-connected micro-environment with the sandy microbial communities and after 60 days in the sandy soil that had also received a silty inoculum in the less-connected environment. In all the other treatment combinations the pool *a* of organic matter was 99% consumed after intermediate durations (Table S6).

The structure of the microbial communities in different micro-environment x soil x day of analysis combinations was significantly different (*P* = 0.01; Fig. 4). The difference among the different communities represented 54% of the total variation in the B-ARISA data. The structure of the microbial communities in their respective native soils were the most different, whilst the communities inoculated into the non-native soil tended to converge towards the sandy communities in the sandy soil.

## Discussion

The data presented here suggest that there is a significant environmental regulatory control on C mineralisation in soil and that this control is significantly greater than the role played by the microbial communities, the effectors of the mineralisation, as others have found for enzyme production³¹ and for the resistance and resilience of microbial decomposition³². The relationship between the organic C content and heterotrophic respiration in soil is well established³³, ³⁴. Here, the silty soil had higher mineralisation rates for both pools of C (Figs 2 and 3, Table S5), suggesting that the organic C of the silty soil was more mineralisable and/or that the conditions in the silty soil were more suitable for microbial activity.

At each incubation matric potential, a larger proportion of the pore network was water-filled in the silty than in the sandy soil (Fig. 1), meaning that a larger proportion of the overall resources in the silty soil were potentially available to the microbial communities, through longer solute and enzyme diffusion pathways, to the microbial communities²³, ³⁵. Assuming that the organic C was homogeneously distributed in the soils and that all the organic C in the water-connected pore space was potentially available to the microbial communities, then 6.21 and 2.65 mg C g⁻¹ soil (total organic C x proportion of pore-network that was water-filled) were potentially available to the microbial communities in the more-connected micro-environment in the silty and sandy soils,
<figure><img src="image_6.png" /><figcaption>**Figure 4.** BGA ordination plot of bacterial community ARISA profiles. The large symbols indicate the group centroids and the small symbols the individual samples. The pink symbols are the microbial communities from the silty soil and the blue symbols the microbial communities from the sandy soil. The square centroid symbols indicate the silty soil and the round centroid symbols the sandy soil. The number in the centroid symbols indicate the day of analysis. The circles are envelopes that incorporate 66% of the individual samples in a given group.</figcaption></figure>

respectively. The amounts of C mineralised during the incubation were therefore equivalent to 6.2 and 6.9% of the potentially available organic in the silty and sandy soils, respectively. The same calculation for the less-connected micro-environment treatment shows that 8.8 and 11.2% of the potentially available organic C in the silty and sandy soils, respectively, were mineralised during the incubation. Were the assumptions made here correct, then, despite the large differences in total C mineralisation in the two soils (Fig. 2), the proportion of available C mineralised was similar, suggesting that it is not necessary to invoke microbial community properties or organic matter composition in order to accurately describe C mineralisation; microbial access to organic C may suffice. Of course, the assumption that organic C is homogeneously distributed in soil is unlikely to be true³⁶, ³⁷. Furthermore, all the C in the water-connected pore space is unlikely to be accessible to the microbial communities, particularly in the silty soil. The greater surface area of the silty soil, due to the higher clay content, is likely to have resulted in lower water film thicknesses and therefore longer solute diffusion pathways³⁸, which may have reduced the potential availability of the organic C. The higher clay content may also have resulted in a higher proportion of the organic matter being adsorbed and hence unavailable to microorganisms²¹. Nevertheless, mapping the availability of organic C to microbial communities may be useful for mechanistic models of soil organic C mineralisation. This would require that the spatial distribution of organic C in relation to connected diffusion pathways to microbial communities at different matric potentials be accounted for.
Differences in C mineralisation may also have arisen because of an incomplete colonisation of the targeted micro-environments, especially in the case of the less-connected micro-environment. If the more connected micro-environment was more widely colonised than the less-connected micro-environment, then one would expect higher mineralisation rates because a higher probability of encounter between microorganisms and substrate³⁹⁻⁴¹. The relationship between matric potential and soil respiration is also believed to be influenced by osmotic stress, gas diffusion and the availability of O₂³⁵. However, the matric potentials used in this study are unlikely to have induced any significant osmotic stress on the microbial communities⁴².

The two microbial communities used in this study were not intrinsically functionally similar. At the beginning of the incubation, when the rate at which the labile pool of organic matter was mineralised was very high, the microbial communities showed significant differences both in the rates at which they mineralised the available organic C (differences in the parameter α) and in their capacities to utilise the range of molecular types available as suggested by the differences in the sizes of the labile organic C pool *a*. The rate at which the labile pool of organic C was mineralised was an order of magnitude greater than the rates that are generally found in bulk mineral soils⁴³, ⁴⁴ but similar to rates found for added labile organic matter²². This suggests that the differences between the microbial communities were expressed when microbial activity was higher than is generally found in bulk soils and that the activity was higher because of a high availability of organic substrate, possibly made available during the sterilisation process⁴⁵.

As the incubation progressed however, the differences in functioning between the two communities diminished (i.e. neither of the parameters related to the mineralisation of the organic C pool *b* were significantly related to the microbial communities). Two studies have shown that the composition of complex microbial inocula that are initially different tend to converge when placed in the same soil environment³¹, ³², which might explain the convergence in functioning of the different communities³⁰. This was not the case here however, as the differences among communities did not converge. In fact, the communities from the silty soil that were inoculated into the silty soil tended to diverge during the incubation (Fig. 4). Therefore, even though the functioning of the communities converged after the readily mineralisable pool of organic C had been consumed, the composition of the microbial communities did not.

Despite the fact that the composition of the communities remained distinct, it is plausible that the active component of the two communities became more similar during the incubation, resulting in a similar mineralisation activity. There may also have been what has been called a “portfolio effect”⁴⁶, where the sum of the different individual activities contained in the communities was the same, regardless of the communities’ compositions. Although plausible, neither of these explanations are likely, as they would require that the active component of the two communities converge in two different soils or that the sum of the constituent activities of the microbial communities be similar in the two different soils.

The mineralization rates of the two microbial communities converged at rates that were an order of magnitude below the mineralization rates of the labile pool of organic C (parameter β versus parameter α; Fig. 3; Table S6) and therefore at rates well below what both communities were capable of. It is possible that the low mineralisation rates at which the communities converged were due to the remaining available organic C being of a complex nature, that required the production of a broad spectrum of enzymes to be produced for decomposition to proceed, thus resulting in similar low mineralisation rates for both communities. However, a number of recent studies have suggested that stable C in soil is not chemically recalcitrant but tends to be composed of relatively simple molecular types⁴⁷, ⁴⁸ that are protected by adsorption to minerals or are inaccessible to microbial degraders.

An equally plausible explanation for the convergence of mineralisation rates is that it occurred due to the environment constraining the intrinsic differences that existed between the communities. The dramatic difference in mineralisation rates of the two organic C pools suggests that constraints on activity during the latter part of the incubation were so great that the manifest differences between the microbial communities' functioning could no longer be expressed. Were this scenario to be true, then the functional redundancy that has often been attributed to microbial communities may, in fact, be the result of environmental constraints: some form of habitat or environment-induced functional equivalence in soil microbial communities. The nature of the constraint or constraints cannot be ascertained, however, it is appealing to suggest that low levels of access to organic substrate constrained the mineralisation activity of the microbial communities. Studies that have identified a relationship between microbial communities and C mineralisation tend to have been carried out in the litter or organic layers of soil¹⁴⁻¹⁶, whilst those that have not have been carried out in mineral soil⁸, ⁹, ¹¹, with lower levels of organic C present.

Habitat filtering is known to be a dominant structuring agent in microbial communities⁴⁹ and other studies that have carried out reciprocal transplant experiments in mineral soil have shown that the structure of microbial communities is principally governed by the soil environment in which they are resident³¹, ³². Here, although the structures of the microbial communities inoculated into the sandy soil were similar, this was not the case for the communities inoculated into the silty soil (Fig. 4). It is known that the pH of soil is one the primary drivers of microbial community structure⁵⁰. The pH of the sandy soil was lower than that of the silty soil (5.2 vs 6.8). These data suggest therefore, that microbial communities from the low pH sandy soil were more resilient to change or that the acidic environment of the sandy soil exerted more pressure on the microbial communities.

We further showed that the micro-environment, in the form of the connectivity of the local environment, also affected the structure of the microbial communities, albeit to a much lesser extent. This tends to corroborate the conclusion drawn in Ruamps *et al.*²⁵ that there is a pore scale microbial biogeography. The differences among treatments and dates accounted for 54% of the total variability in the B-ARISA profiles and the two axes presented in Fig. 4 only 16% of the total variability. Other BGA axes also show differences among treatments and dates (data not shown). Furthermore, due to stochastic processes that can occur during the inoculation and the rapid growth stages at the beginning of the incubation, high inter-replicate variability can be expected in this type of experiment⁸, ¹¹.
## Conclusions

The key take-home messages that emerge from this study are: (1) that microbial respiration in soil may be dependent on the spatial organisation and connectivity of diffusion pathways between substrate and decomposers and (2) that the relationship between microbial communities and soil C dynamics may be dependent on trophic conditions. This is one of the first empirical studies to indicate that the relative contribution of species to community functioning changes with environmental context⁵¹. In environments were the activity of microbial communities is not as restricted as in bulk mineral soil, such as in the litter or organic layers or the rhizosphere, C mineralisation may possibly be related to the composition or diversity of the resident microbial communities. However, when microbial activity is restricted to levels well below the rates at which the communities can metabolise, as might be the case in mineral soil layers, then the composition or diversity of the microbial communities may no longer be consequential because all communities mineralise at the rate at which the abiotic constraints allow. This scenario provides an elegant mechanistic basis for soil C dynamics models.

## Materials and Methods

**Soil sampling and biocidal treatment.** Soil was collected from the surface 10 cm at two sites in Ile de France region (France). Samples were taken from plots under wheat at the “Closeaux” field experiment (INRA research centre in Versailles, France) and from a natural grassland at the CEREEP (Centre de Recherche en Écologie Expérimentale et Prédictive) experimental station (Saint-Pierre-lès-Nemours, France). The Closeaux soil is classified as an Eutric Cambisol (17.4% clay, 53% silt and 29.6% sand) with a pH of 6.8, and organic C and total N contents of 13.5 and 1.23 g.kg⁻¹ soil, respectively. The CEREEP soil is classified as a Sandy Cambisol (6.9% clay, 19% silt and 74.1% sand). The pH of the soil was 5.2 and the organic C and total N contents were 14.7 and 1.19 g.kg⁻¹, respectively. The soils were sieved (<5 mm) to remove stones and plants residues. The soils are referred to as silty (Silty Eutric Cambisol) and sandy (Sandy Cambisol) in the following. The soils had different pore size distributions, as demonstrated by the measured moisture release curves (Fig. S1). The choice of two soils with similar C contents but differing in many other respects was deliberate, as there is a well established relationship between total C and the mineralisation rates of organic C in soil. This relationship *per se* was not the object of this study, hence the choice of soils with similar total C contents.

The soils were sterilised using gamma-irradiation as it is an effective biocidal treatment, but does not disrupt the physical structure of soil nor have a dramatic effect on non-biomass organic matter⁵². It has been shown that gamma irradiation affects primarily molecules of biomass origin⁵³. The soil was air-dried to a matric potential of -1000 kPa and then sealed in polypropylene bags for sterilisation by gamma-irradiation (IONISOS, ZA). The soil was irradiated in a relatively dry state to minimise the impact of irradiation on the chemical stability of soil organic matter⁵². The soil (20 kg) was exposed to a 70 kGy dose (with a minimum guaranteed exposure of 45 kGy). After sterilisation, the soil was placed at 4 °C for 15 days to ensure that the free radicals that may have been generated during gamma irradiation had disappeared prior to the start of the incubation.

**Extraction of microbial communities.** Microbial communities from both soils (silty-MC and sandy-MC) were extracted from non-sterile samples by adding 200 mL ¼ strength Ringer's solution to 20 g soil and shaking the suspension with glass beads for 30 minutes. The soil suspension was left to settle for 10 minutes and the supernatant, containing the microbial suspension, was removed and used to inoculate the sterilized soils.

**Moisture retention curves.** The moisture retention curves of the soils were established in triplicate using Richard's pressure plates (Fig. S1). The water volume at saturation was considered to represent the total pore volume in the soils. The volume of water in the soils at the incubation matric potentials (see below) represented the water-filled volume during the incubation of the soils. The percentage of the total pore volume occupied by the microbial inocula was determined by dividing the volume of the microbial suspension used for inoculation (0.045 cm³ g⁻¹ soil) by the total pore volume.

**Microbial inoculation of sterilized soil.** Sterilized samples of the silty and the sandy were inoculated with either their native microbial suspensions or with microbial suspensions from the other soil. The soils were inoculated at two different matric potentials in order to place the microbial communities in distinct regions of the sterilized soils' pore networks, as described previously⁵⁴⁻⁵⁶. Half the samples were equilibrated at a matric potential below -1600kPa and brought to a matric potential of -100kPa with the microbial suspensions whilst the other half were equilibrated at a matric potential slightly above -100kPa and brought to -3.15kPa with the microbial suspensions. The microbial communities in the first half of the samples were more likely to be found in pores or crevices with a neck diameter <3 µm or along surfaces of larger pores where the water-film thickness was not sufficient to completely immerse a bacterial cell²³. The microbial communities in the second half of the samples were more likely to have been placed in pores with a maximum neck diameter ranging from 6 to 97 µm and were more likely to have been immersed in thicker water-films with longer, more connected diffusion pathways. The microbial communities in the first half of the samples were therefore located in a physical micro-environment in which the diffusion pathways, necessary for substrate and enzymes to come into contact and, ultimately, for microbial activity, were less connected²³,⁵⁷. For simplicity, these two treatments are referred to as “less-connected micro-environment” and “more-connected micro-environment” in the following, although any effects of this treatment on microbial respiration were likely due to more than just the connectivity of the diffusion pathways, such as the local chemical environment and O₂ availability. The same volume of microbial suspension (4.5 mL) was added to each 90 g sterile soil aliquot in all treatments. The same volume of inoculant was used to ensure that the same volume of pore space was inoculated in all treatments. There were 8 treatment combinations (2 soils x 2 communities x 2 micro-environments) and four replicates per treatment combination, resulting in 32 microcosms being prepared for incubation (see next section). In order to accommodate the destructive sampling
necessary for the microbial structure analyses, 4 subs-samples (90 g each) were prepared per replicate. Microbial community structure was determined in 1 of the sub-samples per replicate 6 h after inoculation and the other sub-samples were incubated as described below.

**Incubation and respiration measurements.** The inoculated samples (3 × 90 g) were incubated in 1 L jars fitted with a septum for headspace sampling. All the incubation equipment was autoclaved prior to the addition of the soil samples. The samples were incubated at matric potentials of −3.15 kPa (more-connected micro-environment) or −100 kPa (less-connected micro-environment) at 20 °C in the dark for 310 days. The production of CO₂ was measured 10 times (days 1, 3, 8, 14, 29, 36, 42, 85, 92 and 310) during the incubation by gas chromatography (Agilent 3000 A, Massy, France). The microcosms were flushed under sterile conditions if the headspace CO₂ concentration reached 2% by leaving the jars open under a sterile laminar flow hood for 2 hours. The CO₂ concentration was measured again after flushing. The moisture content of the samples was verified gravimetrically and corrected by adding filter-sterilised water when necessary.

**Microbial community structure analysis.** Microbial community structure was analysed on days 0, 42 and 92 using a bacterial automated rRNA intergenic spacer analysis (ARISA) fingerprinting approach based on the length polymorphism of the bacterial InterGenic Spacer (IGS)⁵⁸ (GenoSol platform). It has been shown that ARISA fingerprinting has as much explanatory power of process rates as next generation sequencing technology⁶.

Total nucleic acids were extracted from 0.5 g soil samples (dry weight equivalent) with the FastDNA Spin kit for soil in combination with the FastPrep FP120 bead beating system (MP-Biomedicals, CA, USA), according to the manufacturer's instructions. Bulk total DNA was purified by elution through Geneclean Turbo columns according to the manufacturer's instructions (MP Biomedicals, CA, USA). The concentration and purity of the resulting DNA were determined spectrophotometrically by measuring the absorbance at 260 and 280 nm and calculating the ratio A260/A280 (NanoDrop ND-1000 spectrophotometer).

The bacterial intergenic spacers were amplified with the primer set: S-D-Bact-1522-b-S-20/ L-D-Bact-132-a-A-18, with 50 ng of DNA as template⁵⁸. The S-D-Bact-1522-b-S-20 primer was labeled at the 5'-end with the IRD800 dye fluorochrome (MWG SA Biotech, Ebersberg, Deutschland) to allow for the detection of the PCR fragments by the LiCor DNA sequencer system (ScienceTec, Les Ulis, France). Data obtained from the 1D-Scan software (Sciencetec) were converted into a table summarizing band presence (that is, peaks) and intensity (that is, height or area of peak) using the PrepRISA program²⁸; 100 peaks, 2 bp resolution and Gaussian peak area were used to provide a robust analysis of bacterial communities. The resulting bacterial-ARISA data matrix (bacterial communities as rows and bands as columns) accounted for the presence/absence and relative intensity of bands.

**Statistical analysis and modelling.** The cumulative respiration curves were described by fitting a two-compartment first-order model (Eq. 1) to the experimental data using a modification of the Levenberg-Marquardt nonlinear least-squares algorithm in the R minpack.lm package⁵⁹:

$$CO_2 = a(1 - e^{-\alpha t}) + b(1 - e^{-\beta t}) \quad (1)$$

were a is a pool of organic C that decomposes rapidly, b is an organic C pool that decomposes at a slower rate and α and β are the respective decomposition rates. The two-compartment first-order model provided a good fit to the respiration curves because the process of sterilisation releases readily mineralisable C, primarily of microbial origin, into the soil solution, which results in a rapid respiratory flush as the microbial inocula colonise the dead biomass⁴⁵,⁵⁵. This is the basis for the fumigation-incubation measurement of the soil microbial biomass⁴⁵. The pool b of organic C can be considered to represent the organic matter unaffected by the biocidal treatment.

In order to determine the effects of soil, micro-environment (matric potential) or microbial community origin on microbial community structure, a between-group analysis (BGA) was carried out on the B-ARISA profiles using the ade4 package⁶⁰ in the R environment (Studio Team, 2015; RStudio: Integrated Development for R. RStudio, Inc., Boston, MA. URL http://www.rstudio.com/). BGA is a supervised classification method that is based on carrying out an ordination of groups of samples rather than of samples. Whilst principal component analysis (PCA) maximises the variance amongst individual samples, BGA maximises the differences among the centroids of groups of samples. The individual sample locations are subsequently projected onto the resulting axes. The B-ARISA profiles were grouped according to treatment combination and day of analysis. The ARISA profiles were transformed using the Hellinger distance prior to carrying out the BGA⁶¹.

**Data availability statement.** The datasets generated during the current study are available from the corresponding author on reasonable request.

## References

1. Coleman, K. et al. Simulating trends in soil organic carbon in long-term experiments using RothC-26.3. *Geoderma* **81**, 29–44 (1997).

2. Parton, W. J., Stewart, J. W. B. & Cole, C. V. Dynamics of C, N, P and S in grassland soils: a model. *Biogeochemistry* **5**, 109–131 (1988).

3. Allison, S. D., Wallenstein, M. D. & Bradford, M. A. Soil-carbon response to warming dependent on microbial physiology. *Nat. Geosci.* **3**, 336–340 (2010).

4. Cheng *et al.* Synthesis and modeling perspectives of rhizosphere priming. *New Phytol.* **201**, 31–44 (2014).

5. Louis *et al.* Microbial diversity indexes can explain soil carbon dynamics as a function of carbon source. *PLOS One* **11**, e0161251 (2016).

6. Graham, E. B. *et al.* Microbes as Engines of Ecosystem Function: When Does Community Structure Enhance Predictions of Ecosystem Processes? *Front. Microbiol.* **7** (2016).

7. Griffiths, B. S. *et al.* Ecosystem response of pasture soil communities to fumigation-induced microbial diversity reductions: an examination of the biodiversity-ecosystem function relationship. *Oikos* **90**, 279–294 (2000).
8. Griffiths, B. S. et al. An examination of the biodiversity-ecosystem function relationship in arable soil microbial communities. *Soil Biol. Biochem.* **33**, 1713–1722 (2001).

9. Juarez, S., Nunan, N., Duday, A.-C., Pouteau, V. & Chenu, C. Soil carbon mineralisation responses to alterations of microbial diversity and soil structure. *Biol. Fertil. Soils* **49**, 939–948 (2013).

10. Nielsen, U. N., Ayres, E., Wall, D. H. & Bardgett, R. D. Soil biodiversity and carbon cycling: a review and synthesis of studies examining diversity-function relationships. *Eur. J. Soil Sci.* **62**, 105–116 (2011).

11. Wertz, S. et al. Maintenance of soil functioning following erosion of microbial diversity. *Environ. Microbiol.* **8**, 2162–2169 (2006).

12. Wertz, S. et al. Decline of soil microbial diversity does not influence the resistance and resilience of key soil microbial functional groups following a model disturbance. *Environ. Microbiol.* **9**, 2211–2219 (2007).

13. Nannipieri, P. et al. Microbial diversity and soil functions. *Eur. J. Soil Sci.* **54**, 655–670 (2003).

14. Strickland, M. S., Lauber, C., Fierer, N. & Bradford, M. A. Testing the functional significance of microbial community composition. *Oecologia* **90**, 441–451 (2009).

15. Salonius, P. O. Metabolic capabilities of forest soil microbial populations with reduced species diversity. *Soil Biol. Biochem.* **13**, 1–10 (1981).

16. Setälä, H. & McLean, M. A. Decomposition rate of organic substrates in relation to the species diversity of soil saprophytic fungi. *Oecologia* **139**, 98–107 (2004).

17. Bonkowski, M. & Roy, J. Soil microbial diversity and soil functioning affect competition among grasses in experimental microcosms. *Oecologia* **143**, 232–240 (2005).

18. Philippot *et al*. Loss in microbial diversity affects nitrogen cycling in soil. *ISME J.* **7**, 1609–1619 (2013).

19. Dungait, J. A. J., Hopkins, D. W., Gregory, A. S. & Whitmore, A. P. Soil organic matter turnover is governed by accessibility not recalcitrance. *Glob. Change Biol.* **18**, 1781–1796 (2012).

20. Schmidt, M. W. I. et al. Persistence of soil organic matter as an ecosystem property. *Nature* **478**, 49–56 (2011).

21. Chenu, C. & Stotzky, G. In Interactions Between Soil Particles and Microorganisms - Impact on the Terrestrial Ecosystem 3–40 (John Wiley & Sons, Ltd, 2002).

22. Juarez, S. et al. Effects of different soil structures on the decomposition of native and added organic carbon. *Eur. J. Soil Biol.* **58**, 81–90 (2013).

23. Or, D., Smets, B. F., Wraith, J. M., Dechesne, A. & Friedman, S. P. Physical constraints affecting bacterial habitats and activity in unsaturated porous media - a review. *Adv. Water Resour.* **30**, 1505–1527 (2007).

24. Killham, K., Amato, M. & Ladd, J. N. Effect of substrate location in soil and soil pore-water regime on carbon turnover. *Soil Biol. Biochem.* **25**, 57–62 (1993).

25. Ruamps, L. S., Nunan, N. & Chenu, C. Microbial biogeography at the soil pore scale. *Soil Biol. Biochem.* **43**, 280–286 (2011).

26. Strong, D. T., De Wever, H., Merckx, R. & Recous, S. Spatial location of carbon decomposition in the soil pore system. *Eur. J. Soil Sci.* **55**, 739–750 (2004).

27. Constancias *et al*. Microscale evidence for a high decrease of soil bacterial density and diversity by cropping. *Agron. Sustain. Dev.* **34**, 1–10 (2014).

28. Blaud, A. et al. Dynamics of bacterial communities in relation to soil aggregate formation during the decomposition of 13C-labelled rice straw. *Appl. Soil Ecol.* **53**, 1–9 (2012).

29. Carson, J. K. et al. Low Pore Connectivity Increases Bacterial Diversity in Soil. *Appl. Environ. Microbiol.* **76**, 3936–3942 (2010).

30. Reed, H. E. & Martiny, J. B. H. Testing the functional significance of microbial composition in natural communities: Functional significance of microbial composition. *FEMS Microbiol. Ecol* **62**, 161–170 (2007).

31. Xun, W. et al. Environmental conditions rather than microbial inoculum composition determine the bacterial composition, microbial biomass and enzymatic activity of reconstructed soil microbial communities. *Soil Biol. Biochem.* **90**, 10–18 (2015).

32. Griffiths, B. S. et al. Functional resilience of soil microbial communities depends on both soil structure and microbial community composition. *Biol. Fertil. Soils* **44**, 745–754 (2008).

33. Franzluebbers, A. J., Haney, R. L., Hons, F. M. & Zuberer, D. A. Active fractions of organic matter in soils with different texture. *Soil Biol. Biochem.* **28**, 1367–1372 (1996).

34. Hassink, J. Effects of soil texture and grassland management on soil organic C and N and rates of C and N mineralization. *Soil Biol. Biochem.* **26**, 1221–1231 (1994).

35. Moyano, F. E., Manzoni, S. & Chenu, C. Responses of soil heterotrophic respiration to moisture availability: An exploration of processes and models. *Soil Biol. Biochem.* **59**, 72–85 (2013).

36. Chenu, C. & Plante, A. F. Clay-associated and clay-bound organic matter in a cultivation chronosequence: abundance and morphology. *Eur. J. Soil Sci.* **57**, 596–607 (2006).

37. Peth *et al*. Localization of soil organic matter in soil aggregates using synchrotron-based X-ray microtomography. *Soil Biol. Biochem.* **78**, 189–194 (2014).

38. Moldrup, P. *et al*. Diffusivity, and Permeability in the Soil Liquid and Gaseous Phases. *Soil Sci. Soc. Am. J.* **65**, 613–623 (2001).

39. Dechesne, A. et al. Biodegradation in a Partially Saturated Sand Matrix: Compounding Effects of Water Content, Bacterial Spatial Distribution, and Motility. *Environ. Sci. Technol.* **44**, 2386–2392 (2010).

40. Kim, H. et al. The spatial distribution of exoenzyme activities across the soil micro-landscape, as measured in micro- and macro-aggregates, and ecosystem processes. *Soil Biol. Biochem.* **91**, 258–267 (2015).

41. Pallud, C. et al. Modification of spatial distribution of 2,4-dichloro- phenoxyacetic acid degrader microhabitats during growth in soil columns. *Appl. Environ. Microbiol.* **70**, 2709–2716 (2004).

42. Manzoni, S. et al. A theoretical analysis of microbial eco-physiological and diffusion limitations to carbon cycling in drying soils. *Soil Biol. Biochem.* **73**, 69–83 (2014).

43. Lerch, T. Z., Nunan, N., Dignac, M.-F., Chenu, C. & Mariotti, A. Variations in microbial isotopic fractionation during soil organic matter decomposition. *Biogeochemistry* **106**, 5–21 (2011).

44. Zak, D. R., Holmes, W. E., MacDonald, N. W. & Pregitzer, K. S. Soil temperature, matric potential, and the kinetics of microbial respiration and nitrogen mineralization. *Soil Sci. Soc. Am. J.* **63**, 575–584 (1999).

45. Jenkinson, D. S. The effects of biocidal treatments on metabolism in soil—IV. The decomposition of fumigated organisms in soil. *Soil Biol. Biochem.* **8**, 203–208 (1976).

46. Bradford, M. A. & Fierer, N. in *Soil Ecology and Ecosystem Services* (eds. Wall, D. H. *et al.*) 189–200 (Oxford University Press, 2012).

47. Kleber, M. *et al*. Old and stable soil organic matter is not necessarily chemically recalcitrant: implications for modeling concepts and temperature sensitivity: slow turnover of labile soil organic matter. *Glob. Change Biol.* **17**, 1097–1107 (2011).

48. Nunan *et al*. Metabolising old soil carbon: Simply a matter of simple organic matter? *Soil Biol. Biochem.* **88**, 128–136 (2015).

49. Hughes Martiny, J. B. *et al*. Microbial biogeography: putting microorganisms on the map. *Nat. Rev. Microbiol.* **4**, 102–112 (2006).

50. Fierer, N. & Jackson, R. B. From the Cover: The diversity and biogeography of soil bacterial communities. *Proc. Natl. Acad. Sci.* **103**, 626–631 (2006).

51. Cardinale, B. J., Nelson, K. & Palmer, M. A. Linking species diversity to the functioning of ecosystems: on the importance of environmental context. *Oikos* **91**, 175–183 (2000).

52. McNamara, N. P., Black, H. I. J., Beresford, N. A. & Parekh, N. R. Effects of acute gamma irradiation on chemical, physical and biological properties of soils. *Appl. Soil Ecol.* **24**, 117–132 (2003).
53. Berns, A. E. et al. Effect of gamma-sterilization and autoclaving on soil organic matter structure as studied by solid state NMR, UV and fluorescence spectroscopy. *Eur. J. Soil Sci.* **59**, 540–550 (2008).

54. Postma, J., van Veen, J. A. & Walter, S. Influence of different initial soil moisture contents on the distribution and population dynamics of introduced Rhizobium leguminosarum biovar Trifolii. *Soil Biol. Biochem.* **21**, 437–442 (1989).

55. Ruamps, L. S. et al. Regulation of soil organic C mineralisation at the pore scale. *FEMS Microbiol. Ecol.* **86**, 26–35 (2013).

56. Wright, D. A., Killham, K., Glover, L. A. & Prosser, J. I. Role of Pore-Size Location in Determining Bacterial-Activity During Predation by Protozoa in Soil. *Appl. Environ. Microbiol.* **61**, 3537–3543 (1995).

57. So, H. B. & Nye, P. H. The effect of bulk density, water content and soil type on the diffusion of chloride in soil. *J. Soil Sci.* **40**, 743–749 (1989).

58. Ranjard et al. Characterisation of bacterial and fungal soil communities by automated ribosomal intergenic spacer analysis fingerprints: biological and methodological variability. *Appl. Environ. Microbiol.* **67**, 4479–4487 (2001).

59. Elzhov T.V., Mullen K.M., Spiess A.N., Bolker B. minpack.pm: R Interface to the Levenberg-Marquardt Nonlinear Least-Squares Algorithm Found in MINPACK, Plus Support for Bounds. R package version 1.1-9. http://CRAN.R-project.org/package=minpack.pm (2015)

60. Dray, S. & Dufour, A. B. The ade4 package: implementing the duality diagram for ecologists. *J. Stat. Soft.* **22**, 1–20 (2007).

61. Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. *Oecologia* **129**, 271–280 (2001).

**Acknowledgements**

This work was funded by the Agence Nationale de la Recherche – project numbers: ANR-09-SYSC-007 and ANR-15-CE01-006-03.

**Author Contributions**

N.N., C.C. and L.S.R. designed the experiment, L.S.R., V.P. and J.L. carried out the experiment, N.N. analysed the data, N.N., C.C. and J.L. wrote the manuscript.

**Additional Information**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-04485-z

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017