Agricultural Practices Modulate the Beneﬁcial Activity of Bacterial-Feeding Nematodes for Plant Growth and Nutrition: Evidence from an Original Intact Soil Core Technique

: Free-living nematodes have beneﬁcial effects on plant growth and nutrition. Exploring how agricultural practices modulate these beneﬁcial effects is still challenging. A study was conducted in Ferralsols from Madagascar from one unmanaged grassland and 16 upland rainfed rice ﬁelds, representative of different agricultural practices: rotation, agroforestry and monoculture. Intact soil cores in plastic cylinders were sampled in the ﬁeld to assess the effects of agricultural practices on changes in plant growth and nutrition induced by the presence of bacterial-feeding nematodes. The soil cores were fumigated to kill the nematodes and moistened with a ﬁltered fresh soil suspension containing only microbial cells. A rice seed was introduced in the core, which was then incubated under natural climatic conditions for 40 days with or without inoculation of the bacterial-feeding nematode Acrobeloides sp. The inoculation of the nematodes induced lower, similar or higher plant biomass and nutrient content in comparison to the control according to the agricultural practices. Positive effects of Acrobeloides sp. on plant functions were frequent in soil cores sampled from ﬁelds with high plant diversity, especially from agroforestry systems. The intact soil core technique appears to be a robust means of mimicking ﬁeld conditions and constitutes a promising tool to assess effects on soil processes of the ecological intensiﬁcation of agricultural practices.


Introduction
Nematodes are cosmopolite and abundant animals [1]. In the soil, they contribute to key functions [2], such as organic matter decomposition [3], nutrient cycling [4], stimulation of microbial activity [5,6] and plant nutrition and growth [7]. Free-living nematodes are known to greatly improve nitrogen (N) and phosphorus (P) availability for plants through the soil microbial loop [8][9][10]. Because N and P are the main macronutrients limiting crop productivity in many terrestrial ecosystems [11,12], improving nutrient availability for crops by promoting the beneficial activities of free-living nematodes is a promising technique for sustainable agriculture.
Identifying agricultural practices that promote the beneficial effects of free-living nematodes, beyond measurements of field nematode abundance and composition, is however still challenging. To date, various methodological tools have been used to study the effects of bacterial-feeding organisms on nutrient availability [9,[13][14][15][16], from the simplest ones, appropriate for disentangling mechanisms, to more complex devices that attempt to mimic field conditions. The simplest experimental models were liquid cultures that could be a soil suspension or a mineral saline solution [14,[17][18][19]. For instance, Griffiths [18] showed, using liquid culture, that N mineralization was greater when bacterial-feeding nematodes (Rhabditis sp.) and protozoan ciliates (Colpoda steinii) were combined than with C. steinii alone. Some studies used sieved soil without plants [17,[20][21][22][23]. More recent work introduced plants in their assays, either with autoclaved sand or a soil-sand mixture [24], or using a humus-sand mixture [25], or sometimes with agarose [16,26]. For instance, Clarholm [10] showed that the presence of protists (from a soil suspension) increased wheat biomass and N content in plant tissues by 78% after six weeks of incubation. Similar patterns have been observed for plant P content [27]. These designs allowed a good control of experimental systems in which the form and the amounts of each nutrient are well known [7,10,28,29].
Whatever the experimental design, quantifying the beneficial activity of soil bacterivores for soil and plant functions is a technical challenge [30]. Most of the studies inoculate soil organisms in disturbed soils, usually 2 mm-sieved soils. Sieving the soil is indeed required to homogenize the substrate, to kill bacterial-feeding organisms and to avoid confounding effects in controlled experiments. It is still however questionable whether these methods reflect actual field conditions and whether they alter the effect of bacterialfeeding organisms on plant growth and nutrition. More recently, Gebremikael et al. [15] used soil cores in plastic tubes with sieved soil compacted to the bulk density found in the field to assess the effects of nematodes on plant functions in conditions close to the field. Fauna were killed by applying gamma irradiation with a dose proven to eradicate nematodes and other fauna without significantly altering the native microbial community. This approach is technically relevant, even if the soil is sieved, but remains less effective for the monitoring and screening of the beneficial activity of bacterivores in the field. In addition, this method is difficult to apply in less developed countries where use of gamma irradiation is not feasible.
It appears important however to be able to measure to what extent a given agricultural practice influences the beneficial activity of soil nematodes for the growth and nutrition of crop plants. For example, by providing fertilizers to the soil, the farmer strongly influences the abundance and composition of nematode communities [31] but also the nitrogen mineralizing activity of bacterial-feeding nematodes [8]. Practices that can alter soil porosity can also influence nematode dispersal and nematode-root interactions [32]. The beneficial activities of nematodes on plants can also be influenced by the presence of other soil organisms with which they interact in trophic and non-trophic manners [33][34][35][36]. Thus, the effects of free-living soil nematodes on plant growth are likely to vary depending on how agricultural practices alter the chemical, physical and biological properties of the soil nematode habitat.
In the highlands of Madagascar, the soils of hills cropped with rainfed rice are acidic nutrient-depleted Ferralsols [37]. Poor smallholder farmers do not have access to expensive mineral fertilizers. In these conditions, an effective intensification of the soil functions involved in improved soil nutrient availability is required. In this study, we aimed at testing a new method based on intact, undisturbed soil cores from various rainfed rice fields to assess the potential beneficial activity of the bacterial-feeding nematode Acrobeloides sp. for rice growth and nutrition. This method was developed to limit disturbance. We first hypothesized that the presence of Acrobeloides sp. in the rhizosphere of rice has a positive effect on plant growth and nutrition as reported in the literature [38]. However, the size of the beneficial effect of the nematode on plant nutrition and growth is likely to vary depending on the management of the field from which the soil cores are sampled.

Selection of Upland Rice Fields
The study was carried out in upland rainfed rice (Oryza sativa) systems, in farmers' fields. A total of 17 fields were selected in Madagascar in the region of Itasy ( Figure S1 in Supplementary Materials). We selected only upland rainfed rice fields and a natural grassland as a reference non-cultivated area. Rice was grown as a monoculture, as an intercrop or in rotation (Table 1). For each field, we collected agronomic data (e.g., crop history, crop management, soil management). The 17 fields were gathered in five agronomic groups depending on crop diversity: (i) agroforestry (F1, F2, F4, F5, F6 and F7), (ii) rotation (F3, F8, F10, F12 and F13), (iii) rice monoculture (F11, F14, F15 and F16), (iv) fallow (F9) and (v) grassland (F17). The F3 field, which presents intercropping, was classified as in rotation. These fields were thus representative of the cropping practices developed by farmers in the highlands of Madagascar, i.e., association with legumes or trees, organic inputs, fertilization and soil tillage ( Table 1). The soils from the fields were all ferrallitic soils. Nine fields were sampled in 2016, and eight were sampled in 2017.

Soil Sampling
For each campaign, seven days before the beginning of the assay, three intact soil cores (top 10 cm depth) were sampled in the middle of each field, near a plant, with PVC cylinders (Figure 1). The PVC cylinders (4.5 cm diameter and 12 cm height) were open from the top and from the bottom and drilled on the sides (12 holes~5 mm in diameter throughout the height of the tube) to facilitate air exchange between the soil inside the cylinder and the atmosphere. These three soil cores were labeled "T -7 -cores". A composite sample made of 5 soil samples, also in the 0-10 cm depth, was collected for the analysis of soil abiotic parameters. This composite soil sample was collected in the middle of each field. The soil cores and the composite sample were carefully transported to the laboratory. One week later, 10 soil cores, designated "T 0 -cores", were sampled per field in the same area using identical PVC cylinders. They were used to conduct the assay. Soil bulk density was determined in the field using a soil core sampler, which was pushed into the soil at 10 cm depth. The soil was sieved at 2 mm and oven dried at 105 • C for 24 h in the laboratory.  (Figure 1). The PVC cylinders (4.5 cm diameter and 12 cm height) were open from the top and from the bottom and drilled on the sides (12 holes ~5 mm in diameter throughout the height of the tube) to facilitate air exchange between the soil inside the cylinder and the atmosphere. These three soil cores were labeled "T-7-cores". A composite sample made of 5 soil samples, also in the 0-10 cm depth, was collected for the analysis of soil abiotic parameters. This composite soil sample was collected in the middle of each field. The soil cores and the composite sample were carefully transported to the laboratory. One week later, 10 soil cores, designated "T0-cores", were sampled per field in the same area using identical PVC cylinders. They were used to conduct the assay. Soil bulk density was determined in the field using a soil core sampler, which was pushed into the soil at 10 cm depth. The soil was sieved at 2 mm and oven dried at 105 °C for 24 h in the laboratory.

Soil Analyses
The three "T-7-cores" were used to assess the soil water-holding capacity (WHC) within the cylinder (Figure 1). The soil cores were saturated with water in a basin for 2 h and then placed on a sand tray for 1 h. Soil moisture was determined after oven drying an aliquot (20 g) at 105°C for 24 h. The composite soil sample was air-dried until constant weight and stored at room temperature for further analyses. Three sub-samples were randomly taken to perform the physicochemical measurements. Soil pH was measured in water (pH H2O) (1:2.5 mass ratio) after 30 min shaking on a back-and-forth agitator. Particle size distribution was assessed for five fractions: clay (0-2 µm), fine silt (2-20 µm), coarse silt (20-50 µm), fine sand (50-200 µm) and coarse sand (200-2000 µm). Cation exchange capacity (CEC) was measured after saturation with a cobaltihexamine chloride solution. CEC was estimated by measuring the Cobalt remaining in solution by flame atomic absorption spectrometry (iCE 3000 Series, Thermo Scientific, Japan). We also measured three exchangeable cations: Ca 2+ , Mg 2+ and K + . Soil organic carbon (C) and N contents were assessed using Walkley and Black [39] and Kjeldahl [40] methods, respectively. Total phosphorus (P) content was determined after perchloric acid digestion [41]. Labile inorganic P in water was extracted with an ion-exchange resin filled with bicarbonate (0.5M) [42]. P concentration was analyzed by colorimetry with malachite green [43]. Nematodes were extracted using a modified Seinhorst [44] elutriation method from 250 g moist soil.
The number of free-living nematodes was counted in 5 mL of the final volume of water collected.

Soil Analyses
The three "T -7 -cores" were used to assess the soil water-holding capacity (WHC) within the cylinder (Figure 1). The soil cores were saturated with water in a basin for 2 h and then placed on a sand tray for 1 h. Soil moisture was determined after oven drying an aliquot (20 g) at 105 • C for 24 h. The composite soil sample was air-dried until constant weight and stored at room temperature for further analyses. Three sub-samples were randomly taken to perform the physicochemical measurements. Soil pH was measured in water (pH H 2 O) (1:2.5 mass ratio) after 30 min shaking on a back-and-forth agitator. Particle size distribution was assessed for five fractions: clay (0-2 µm), fine silt (2-20 µm), coarse silt (20-50 µm), fine sand (50-200 µm) and coarse sand (200-2000 µm). Cation exchange capacity (CEC) was measured after saturation with a cobaltihexamine chloride solution. CEC was estimated by measuring the Cobalt remaining in solution by flame atomic absorption spectrometry (iCE 3000 Series, Thermo Scientific, Japan). We also measured three exchangeable cations: Ca 2+ , Mg 2+ and K + . Soil organic carbon (C) and N contents were assessed using Walkley and Black [39] and Kjeldahl [40] methods, respectively. Total phosphorus (P) content was determined after perchloric acid digestion [41]. Labile inorganic P in water was extracted with an ion-exchange resin filled with bicarbonate (0.5M) [42]. P concentration was analyzed by colorimetry with malachite green [43]. Nematodes were extracted using a modified Seinhorst [44] elutriation method from 250 g moist soil. The number of free-living nematodes was counted in 5 mL of the final volume of water collected.

Biological Materials
We used the upland rice cultivar B22, which was introduced in Madagascar in the early 1980s. This cultivar is widely used by farmers in Itasy [45,46]. Its optimal altitude of cropping is between 700 to 1200 m. Monoxenic populations of nematodes (Acrobeloides sp.), isolated from Lazaina (18 • 46 55 59 • S, 47 • 32 46 3 • N, 1274 m.a.s.l. Madagascar), were grown on E. coli following the procedure described by Irshad et al. [26]. Acrobeloides is a ubiquitous genus that dominates bacterivorous nematode communities in upland rice systems in the Malagasy highlands [47]. The nematodes were collected by rinsing a Petri dish with sterile water. The nematode suspension was centrifuged at 4000 rpm for 10 min and washed with sterile water. The number of individuals in the washed nematode suspension was counted under a binocular stereoscope. Half of this suspension was sampled, and the volume was adjusted to provide an expected number of nematodes per mL (solution A). The remaining half was filtered at 8 µm (solution B) and used for the control. The absence of nematodes (including eggs) in solution B was confirmed.

Soil Core Preparation for the Assay
The first step aimed to kill the nematodes. The T 0 -cores were moistened with distilled water at 50% of their water-holding capacity and fumigated with chloroform vapor for 2 h ( Figure 1). We repeated the 2 h chloroform fumigation after 24 h. Because fumigation with chloroform vapor induces a lysis of microbial cell membranes [48], a microbial suspension was prepared for each field by vigorously mixing 100 g of non-fumigated soil in 400 mL of distilled sterile water for 10 min. The soil from the suspension was decanted by sedimentation for 1 h, and the suspension was filtered at 8 µm in order to eliminate freeliving nematodes. Two days after the second sampling, the soil cores were moistened at 100% of their WHC with a fresh microbial suspension, i.e., each soil core was re-inoculated with its own microbial bacterial community. The soil cores were incubated for seven days in the dark at 24 • C.

Experimental Design and Setup
The soil cores were split into two groups (5 cores per group): (i) without the nematode Acrobeloides sp. corresponding to the control (-N) and (ii) with the nematodes (+N). For each field, five replicates were thus set up for each nematode treatment. Seven days after the addition of the microbial suspension, one rice seed was sown in each soil core ( Figure 1). For the +N treatment, 10 individuals of Acrobeloides per gram of dry soil were added by adding 10 mL of solution A. For the -N treatment, 10 mL of solution B was added to the core. The soil cores were then incubated in the experimental station of the Laboratoire des Radio Isotopes (University of Antananarivo, Madagascar) on benches located outside. Soil water content was adjusted twice per day for optimal nematode and plant growth. For the first campaign (February to April 2016), the temperature was 17-25 • C, the mean air humidity was 91% and the cloud coverage was 51%. For the second campaign (April and May 2017), the temperature was 19-23 • C, the air humidity was 74% and the cloud coverage was 25% (Ampandrianomby Weather Center, Antananarivo, Madagascar).

Soil and Plant Analyses at the End of the Assay
After 40 days of growth, shoots and roots of each plant were collected separately. The roots were gently washed with distilled water to remove the adhering soil. Dry root and shoot biomasses were determined after oven drying at 60 • C for 48 h. Total C and N contents in plant tissues were determined using a CHNS/O elemental analyzer. Total P content in plant tissues was determined after crushing and mineralization of plant tissues with nitric acid 65% (US EPA Method 3051). P concentration was analyzed in the mineralized solution by colorimetry with the vanado-molybdate method. A soil aliquot of 150 g from each core was used to extract the nematodes using the same method described earlier. The total nematode density was determined using a stereoscope.

Statistical Analyses
To explore the variability in soil properties among the fields, we conducted a principal component analysis using soil abiotic data [49]. For each field and nematode treatment, the mean and standard deviation (n = 5) of each variable were computed. We then calculated the effect size of the nematode treatment on plant growth and nutrition with Cohen's d (1988) and the following formula (1): where M1 is the mean of each variable for the +N treatment, M2 is the mean of each variable for the -N treatment and pooled SD, the square root of the common variance of the populations. Cohen's d effect size can be described as being "very small, small, medium, large, very large or huge" according to the magnitude of d. The d value can thus be used as a proxy of the beneficial activity of Acrobeloides sp. To determine if the value of the mean effect size was significantly different from zero, we performed a Student test (statistical significance set at 0.05). One-Way ANOVA and post hoc Tukey HSD tests were used to test for significant differences among the fields in final soil nematode abundance. The normality of model residuals was verified using the Shapiro test. Finally, partial least squares regression (PLS-R) models were used to evaluate the relationships between the agronomic and soil abiotic variables and the effect of Acrobeloides sp. on rice growth and nutrition [50,51]. Cohen's d was used as the response variable. Model strength was assessed by the proportion of variance in the dependent variable that can be predicted by the model (Q 2 ). At the 5% level, the test is considered significant when Q 2 exceeds a critical value of 0.097 [50]. The variable importance of projection (VIP) reflects the relative influence of each explanatory variable in the construction of the model. Explanatory variables with VIP > 1 are the most relevant and significant for explaining the dependent variable. We used three PLS-R models to evaluate the relationships between the agronomic and soil abiotic variables as predictors of the effect of Acrobeloides sp. on plant total biomass (first model), plant total nitrogen (second model) and plant total phosphorus (third model). For each model, we first included all the explanatory variables. Then, a second model was achieved with only the explanatory variables with VIP > 1 in order to improve the quality of the models. All tests were done with R, except those with the PLS-R models, which were performed using TANAGRA 1.4.5 [52].

Soil Properties of the Selected Fields
All fields exhibited soil properties typical of Ferralsols, with low organic matter contents and acidic pH H 2 O (Table S1 in Supplementary Materials). On average, the soils contained 26 g.kg −1 total C, 2.1 g.kg −1 total N, 0.72 cmol + .kg −1 Ca and 0.26 cmol+.kg −1 Mg. The mean pH H 2 O was 5.04. The sampled fields were also characterized by a low C/N ratio (12) and high C/P ratio (31). However, although all the soils were Ferralsols, we observed a large variation for some soil parameters among the fields with a coefficient of variation (CV) higher than 15%. Olsen P (CV = 83%) and K contents (CV = 71%) were the soil parameters with the highest variations. The principal component analysis performed on soil parameters reveals the large dispersion on the first two axes, explaining 44.6% of the total inertia, of soil properties ( Figure 2). The first axis (relative inertia = 22.7%) opposes fields with high (coarse and fine) sand content, bulk density and C/P ratio (negative scores) to fields with high P, N, K and clay contents (positive scores) ( Figure 2B). The second axis (relative inertia = 21.9%) opposes fields with high Ca, Mg, Pi and CEC contents (positive scores) to fields with high total C and N, fine silt contents, Ca/Mg ratio and pH KCl .

Final Soil Nematode Abundance
Forty days after inoculation, we did not observe nematodes in the controls. However, the final abundance of Acrobeloides sp. in the inoculated soil cores varied between fields ( Figure 3). Nematode abundance was the highest in F5 (around 110 individuals per g dry soil) and the lowest in grassland (~3 individuals per g dry soil).

Effects of Bacterial-Feeding Nematodes on Plant Growth
The inoculation of Acrobeloides sp. in the soil cores induced negative, neutral or positive effects on plant growth and nutrition (Figure 4). Nematode inoculation induced significantly higher shoot biomass than the control in 10 of the 17 fields with an average positive effect of 35% more than the control ( Figure 4A).

Final Soil Nematode Abundance
Forty days after inoculation, we did not observe nematodes in the controls. However, the final abundance of Acrobeloides sp. in the inoculated soil cores varied between fields ( Figure 3). Nematode abundance was the highest in F5 (around 110 individuals per g dry soil) and the lowest in grassland (~3 individuals per g dry soil).

Final Soil Nematode Abundance
Forty days after inoculation, we did not observe nematodes in the controls. However, the final abundance of Acrobeloides sp. in the inoculated soil cores varied between fields ( Figure 3). Nematode abundance was the highest in F5 (around 110 individuals per g dry soil) and the lowest in grassland (~3 individuals per g dry soil).

Effects of Bacterial-Feeding Nematodes on Plant Growth
The inoculation of Acrobeloides sp. in the soil cores induced negative, neutral or positive effects on plant growth and nutrition (Figure 4). Nematode inoculation induced significantly higher shoot biomass than the control in 10 of the 17 fields with an average positive effect of 35% more than the control ( Figure 4A).  Figure 3. Final nematode abundance in the inoculated soil cores after 40 days of plant growth. Whiskers correspond to standard deviation. Different letters (a, b and c) indicate significant differences between the fields according to Tukey HSD test (p value < 0.05).

Effects of Bacterial-Feeding Nematodes on Plant Growth
The inoculation of Acrobeloides sp. in the soil cores induced negative, neutral or positive effects on plant growth and nutrition (Figure 4). Nematode inoculation induced significantly higher shoot biomass than the control in 10 of the 17 fields with an average positive effect of 35% more than the control ( Figure 4A). More precisely, we recorded the most positive effect sizes of Acrobeloides sp. on shoot biomass in soil cores sampled from fields characterized by the presence of trees. The effect size for F1 was very large (d= 1.9), whereas the effect size for F5 was small (d= 0.4). On the other hand, the inoculation of nematodes had a negative effect on shoot biomass for three of the 17 fields, with a mean negative effect of around 19%. These fields were characterized by the presence of legumes and mostly corresponded to rotated crops. The patterns of variation for root biomass were close to those found for shoot biomass ( Figure 4B). Nematode inoculation induced significantly higher and lower root biomasses than the control in seven and five fields, respectively. Again, the highest positive effects were found for the soil cores sampled in agroforestry fields, while fields F10 (sweet potato) and F8 (maize with jugo bean) exhibited the highest negative effect sizes.

Effects of Bacterial-Feeding Nematodes on Plant Nutrition
The bacterial-feeding nematode increased rice total N uptake in seven fields with an average positive effect of 39% compared to the control ( Figure 4C). The effect size was very large for F2 (d = 1.94) and medium for F12 (d = 0.52). The effect of the bacterivorous nematodes was significantly negative for five of the 17 fields with a mean negative effect of around 19%. More precisely, nematodes significantly decreased rice total N uptake, with medium to large effect sizes, for the microcosms from fields under a rice monoculture system (F14, F15 and F16) and in the rotated crops (F8 and F10).
Rice total P uptake increased following nematode inoculation for nine of the 17 fields with an average positive effect of 62% more than the control ( Figure 4D). The effect size was very large for F3 (d = 1.5) and medium for F2 (d = 0.6). The nematodes had a significantly negative effect for six of the 17 fields resulting in an average of 30% less total P uptake than the control. The bacterivorous nematodes significantly decreased rice total P uptake for the microcosms from fields under a rice monoculture system (F14, F15 and F16) and in the rotated crops (F8 and F10).

Agronomic and Soil Drivers of the Beneficial Activity of Bacterial-Feeding Nematodes for Plant Growth and Nutrition
According to the PLS regression models, both agronomic and soil abiotic properties were significantly related to the effects of bacterial-feeding nematodes on plant growth  More precisely, we recorded the most positive effect sizes of Acrobeloides sp. on shoot biomass in soil cores sampled from fields characterized by the presence of trees. The effect size for F1 was very large (d = 1.9), whereas the effect size for F5 was small (d = 0.4). On the other hand, the inoculation of nematodes had a negative effect on shoot biomass for three of the 17 fields, with a mean negative effect of around 19%. These fields were characterized by the presence of legumes and mostly corresponded to rotated crops. The patterns of variation for root biomass were close to those found for shoot biomass ( Figure  4B). Nematode inoculation induced significantly higher and lower root biomasses than the control in seven and five fields, respectively. Again, the highest positive effects were found for the soil cores sampled in agroforestry fields, while fields F10 (sweet potato) and F8 (maize with jugo bean) exhibited the highest negative effect sizes.

Effects of Bacterial-Feeding Nematodes on Plant Nutrition
The bacterial-feeding nematode increased rice total N uptake in seven fields with an average positive effect of 39% compared to the control ( Figure 4C). The effect size was very large for F2 (d = 1.94) and medium for F12 (d = 0.52). The effect of the bacterivorous nematodes was significantly negative for five of the 17 fields with a mean negative effect of around 19%. More precisely, nematodes significantly decreased rice total N uptake, with medium to large effect sizes, for the microcosms from fields under a rice monoculture system (F14, F15 and F16) and in the rotated crops (F8 and F10).
Rice total P uptake increased following nematode inoculation for nine of the 17 fields with an average positive effect of 62% more than the control ( Figure 4D). The effect size was very large for F3 (d = 1.5) and medium for F2 (d = 0.6). The nematodes had a significantly negative effect for six of the 17 fields resulting in an average of 30% less total P uptake than the control. The bacterivorous nematodes significantly decreased rice total P uptake for the microcosms from fields under a rice monoculture system (F14, F15 and F16) and in the rotated crops (F8 and F10).

Agronomic and Soil Drivers of the Beneficial Activity of Bacterial-Feeding Nematodes for Plant Growth and Nutrition
According to the PLS regression models, both agronomic and soil abiotic properties were significantly related to the effects of bacterial-feeding nematodes on plant growth and nutrition ( Figure 5). The effect size of nematodes on plant growth was well explained by agronomic and soil data, with an R 2 of 0.66 (quality of fit) and a Q 2 of 0.34 (predictive quality). The slope (positive or negative) indicates the sign of the relationship between the predictors and the response variable. The presence of trees and rice in the field the year we sampled the soil and soil pH (positive slope) were the most important predictors. We also observed lower, but still high, variable importance in projection (VIP) values for soil Mg content, soil K/Mg ratio, coarse sand content, Ca/Mg and soil CEC ( Figure 5A). Plant N nutrition after the inoculation of Acrobeloides sp. was the best model with an R 2 of 0.94 and Q 2 of 0.79 ( Figure 5B). The predictors were mostly the presence of trees, the soil Mg and Ca contents and the soil pH H 2 O. Lower VIP values (~0.7-0.8) were observed for the soil C and N contents, the soil K/Mg ratio, the presence of rice or rotated crops at the sampling date and coarse sand content ( Figure 5B). Plant P nutrition was well explained (R 2 = 0.66, Q 2 = 0.49) by soil total C and N contents (negative slopes), the rotated crops at sampling date and soil Ca content (positive slopes, Figure 5C). We also found high VIP values for the soil K/Mg ratio, the presence of rice, soil Mg content and coarse sand content, pH and the presence of trees ( Figure 5C). Sustainability 2021, 13, x FOR PEER REVIEW 9 of 14 and nutrition ( Figure 5). The effect size of nematodes on plant growth was well explained by agronomic and soil data, with an R 2 of 0.66 (quality of fit) and a Q 2 of 0.34 (predictive quality). The slope (positive or negative) indicates the sign of the relationship between the predictors and the response variable. The presence of trees and rice in the field the year we sampled the soil and soil pH (positive slope) were the most important predictors. We also observed lower, but still high, variable importance in projection (VIP) values for soil Mg content, soil K/Mg ratio, coarse sand content, Ca/Mg and soil CEC ( Figure 5A). Plant N nutrition after the inoculation of Acrobeloides sp. was the best model with an R² of 0.94 and Q² of 0.79 ( Figure 5B). The predictors were mostly the presence of trees, the soil Mg and Ca contents and the soil pH H2O. Lower VIP values (~0.7-0.8) were observed for the soil C and N contents, the soil K/Mg ratio, the presence of rice or rotated crops at the sampling date and coarse sand content ( Figure 5B). Plant P nutrition was well explained (R 2 = 0.66, Q 2 = 0.49) by soil total C and N contents (negative slopes), the rotated crops at sampling date and soil Ca content (positive slopes, Figure 5C). We also found high VIP values for the soil K/Mg ratio, the presence of rice, soil Mg content and coarse sand content, pH and the presence of trees ( Figure 5C). "Rice" = rice as the crop at soil sampling date; "Rot" = rotated crops at the soil sampling date; "pHw" = soil water pH; "CSa" = the percentage in coarse sand; "C" = total soil C content; "N" = total soil N content; "Ca" = soil Ca content; "Mg" = soil Mg content; "Ca/Mg" = soil Ca/Mg ratio; "K/Mg" = soil K/Mg ratio.

Diversity of Practices in Upland Rainfed Rice Systems
We aimed at sampling fields with contrasting agricultural practices within the same cropping system, i.e., upland rainfed rice in high altitude tropical conditions, to test the effect of nematode inoculation in various agronomic conditions but within the same pedoclimatic context. We obtained a large panel of land and soil management practices Black circles are fields. Grey area corresponds to 95% confidence interval. "Trees" = presence of trees; "Rice" = rice as the crop at soil sampling date; "Rot" = rotated crops at the soil sampling date; "pHw" = soil water pH; "CSa" = the percentage in coarse sand; "C" = total soil C content; "N" = total soil N content; "Ca" = soil Ca content; "Mg" = soil Mg content; "Ca/Mg" = soil Ca/Mg ratio; "K/Mg" = soil K/Mg ratio.

Diversity of Practices in Upland Rainfed Rice Systems
We aimed at sampling fields with contrasting agricultural practices within the same cropping system, i.e., upland rainfed rice in high altitude tropical conditions, to test the effect of nematode inoculation in various agronomic conditions but within the same pedoclimatic context. We obtained a large panel of land and soil management practices within the same area and soil type. Principal component analysis is a relevant global visualization tool to summarize data set variance in low-dimensional space. The PCA allowed us to determine which soil properties best explained the total variance. The PCA performed on the soil parameters showed a wide dispersion in some soil abiotic parameters, even if all the fields were on Ferralsols. We observed a significant textural gradient not associated with cropping systems but related to soil total P content. Unsurprisingly, high P content in soil was associated with high clay content, in line with the high P sorption ability of clay [53]. The presence of trees was related to higher Ca, Mg and K contents and soil pH as already observed in the literature [54]. We can thus conclude that our final field sampling fits well the criteria we defined initially.

Variable Effects of Bacterial-Feeding Nematodes on Plant Growth and Nutrition
Using an original intact soil core technique, we observed variable effects of the inoculation of the nematode Acrobeloides sp. on rice growth and nutrition. We observed significant positive and negative effects of the nematode on rice biomass in 50% and 23% of all the agronomic situations, respectively. The mean overall effects of Acrobeloides sp. on rice shoot and root biomass were +19% and +5%, respectively. In a meta-analysis, Trap et al. [38] reported larger effects of bacterivores (nematodes and protists) on plant growth. They showed that bacterivores significantly increased shoot and root biomass with an average positive effect of 27% and 21%, respectively, in comparison to the control. They mostly reported positive or neutral effects of bacterivores on plant growth. Similarly, in our study Acrobeloides sp. had either significant positive (47% of the fields) or negative (32% of the fields) effects on rice nutrition. The mean overall effects of Acrobeloides sp. on total plant N amount and P amount were around +13% and +21%, respectively. Again, the meta-analysis conducted by Trap et al. (2016) showed greater overall bacterivore effects on the total plant N and P amounts in shoots and roots. The authors found that bacterivores increased shoot total N amount, root total N amount, shoot total P amount and root total P amount by 59%, 28%, 38% and 55% compared to the control, respectively. They reported very few negative effects on plant total N amount (only 3% of their observations).
We supposed that the use of intact soil cores should better allow the conservation of soil properties and differences in soil functioning between fields. Conversely, soil sampling, sieving, homogenization before filling the microcosms and other technical procedures must certainly converge towards similar soil situations limiting the differences between plots. Effect variance is indeed larger in our study, with more negative and null effect sizes, while the use of sieved soil, agar or liquid culture results in more homogeneous experimental conditions, likely to promote fewer variable patterns. The negative effects of the nematodes on plant growth and nutrition are more difficult to explain. However, the B22 rice cultivar is known to attract nematodes to a lesser extent than other cultivars [55]. It is possible that in specific soil conditions, the presence of the nematodes in high density close to the plant roots caused high rhizodeposition rates as observed by Sundin et al. [56] in comparison with plants without nematode inoculation. This excess of carbon loss may not be compensated by higher nutrient availability following nematode inoculation in comparison to the control [10], resulting in a net carbon loss and thus a lower final plant biomass. This hypothesis could be tested by measuring root exudation rate in the presence and absence of the nematode in various soil conditions. Nevertheless, in this case, the causes behind an absence in soil nutrient availability increase remain to be elucidated. To conclude, our method revealed that (1) the effect of bacterial-feeding nematodes on plant growth and nutrition may be more variable than expected from the literature according to environmental conditions, and (2) agricultural practices can modify the effect size of the beneficial activity of nematodes for plant growth and nutrition.

Soil Drivers of the Beneficial Activity of Bacterial-Feeding Nematodes for Plant Growth and Nutrition
The PLS models showed that soil pH was one of the best predictors of the effect of Acrobeloides sp. on plant growth and nutrition, especially N ( Figure 5). More precisely, the effects of Acrobeloides sp. on plant growth and nutrition increased with increasing pH H2O. Soil pH is a well-known regulatory factor of soil nutrient availability [57]. Nutrients released by bacterivores can be absorbed by plants, re-immobilized within the microbial biomass or adsorbed to soil constituents according to the pH [8,58]. Ferralsols are characterized by high contents of Al and Fe oxides able to precipitate inorganic P as minerals such as variscite (AlPO 4 ·2H 2 O) and strengite (FePO 4 ·2H 2 O), whose solubility increases with pH [59]. On the other hand, free Al concentration, and thus Al toxicity, decrease dramatically at pH 5.0 [60]. We thus expected greater effects of bacterial-feeding nematodes on nutrient availability for soil pH values above 5.5-6, i.e., improving nutrient availability, especially inorganic P, and reducing Al toxicity for rice.
We also found that the effect of Acrobeloides sp. on plant total N and P amount increased with increasing values of soil Mg and Ca content. A nutrient-omission trial conducted on Ferralsols collected from the same region showed high Mg and Ca deficiencies for rice and revealed a N-P-Ca-Mg co-limitation [37]. Furthermore, Mg is known to be a crucial element for the uptake of N, especially nitrate [61]. It is thus possible that the increase in N and P availability caused by nematodes was more pronounced when Mg and Ca limitations were reduced. In contrast to the two previous models, the effect of nematodes on plant P nutrition was mostly explained by soil C and N contents, i.e., the effects of Acrobeloides sp. on plant total P amount increased with decreasing values of soil total N and C. Because soil C is a limiting factor of microbial growth and activity, it is possible that a large number of P ions released by the bacterial-feeding nematodes were rapidly re-immobilized within the microbial biomass.
An interesting observation is that, despite high variations in final abundance of soil nematodes among the fields, this variable was not selected by the PLS model to explain the effect sizes. We thus concluded that the beneficial activity of nematodes for plant growth and nutrition was not related to nematode abundance. We noted that the final abundance of nematodes was higher in soil from the second campaign. An important soil parameter explaining nematode survival rates is the soil texture [62]. The soil particle size distribution of the fields sampled during the second campaign was coarser in comparison to the first one: clay −21%, fine silt −17%, coarse silt +6%, fine sand +14% and coarse sand +83%. As shown by Jiang et al. [63], the average number of total nematodes increases with increasing aggregate size. It is thus possible that the final abundance of Acrobeloides sp. was related to changes in soil texture during the two campaigns. We exclude the effect of variable climate between the two campaigns, even if the soil cores were not incubated in a greenhouse. In both campaigns, the temperature reached the temperature niche breadth for Acrobeloides sp. ranging from 15 to 30 • C [64], and soil moisture was monitored and corrected every day. To conclude, our PLS models revealed that agricultural practices may affect the beneficial activity of nematodes for plant nutrition and growth mostly by changing soil pH and cation (Ca, Mg and K) contents but also total C and N contents.

Agronomic Drivers of the Beneficial Activity of Bacterial-Feeding Nematodes for Plant Growth and Nutrition
The PLS models revealed that Acrobeloides sp. mostly increased plant growth and total N amount in the microcosms from fields under agroforestry. We do not have precise elements to explain this pattern, and we can only speculate on certain mechanisms. Because soil bulk density was not a significant predictor in the PLS models, we can exclude the hypothesis that tree roots might improve nematode dispersal rate within the soil cores by increasing soil porosity and creating macropores. In contrast, as discussed earlier, soil Ca, K and Mg contents were good predictors; Pardon et al. [54] showed that the presence of trees increases soil nutrient content, resulting from the input of leaf litter or from deeper soil layers [65]. We also observed a high positive association between agroforestry and Ca, K and Mg contents on the PCA (Figure 2). It is thus possible that the soil solution and microbial biomass nutrient contents near the tree roots are improved. It would be informative (i) to investigate spatial variation in soil carbon, nutrient, microbial biomass and nematodes within agroforestry fields according to tree distance and (ii) to conduct similar nematode assay with soil core sampling at various distances from the trees.
We also found that rotation (the rotated crop) was a good predictor for the effect of Acrobeloides on plant total P amount. The precise quantity and type of fertilizers used by the farmers depend on the type of crop. The addition of crop residues stimulated bacterial growth by increasing the availability of labile carbon derived during initial phases of litter decomposition [30]. Thus, higher quantities of nutrients are stored in the microbial biomass, potentially unlocked by the activity of the bacterivores. To conclude, plant diversity in the field, and especially the presence of trees, was a key practice that affected the beneficial activity of nematodes for plant nutrition and growth.

Conclusions: Limits and Perspectives
The new method based on intact soil cores proposed in this study appears relevant and robust to test the effect of the bacterial-feeding nematode Acrobeloides sp. on rice growth and nutrition. Using this method, we were able to assess the effect of agricultural practices on the soil microbial loop driven by Acrobeloides sp. Indeed, a main finding of our study is the variable effect of the bacterial-feeding nematode Acrobeloides sp. on plant growth and function according to agricultural practices. We thus believe that driving the beneficial activity of free-living nematodes by appropriate agricultural practices is feasible. Agricultural practices based on improved plant diversity were highly correlated with beneficial activity of the nematodes. Furthermore, higher values in soil pH, and cation contents were associated with positive effect sizes of the nematode on plant growth and function. The nematode assay was time-consuming and required specific materials and skills. We quickly reached the maximum number of fields to investigate per year due to technical constraints. This is the reason why only 17 fields were investigated (170 soil cores) and why we split the samples in two campaigns. In this study, we did not measure the bacterial community structure or biomass, although it may explain the effects of nematodes on plants. It is likely that practices influence this microbial community with possible consequences on the effects of nematodes on rice. It would therefore be interesting to measure this aspect with a view to improving our understanding of the determination of the microbial loop. Finally, the identification of indicators of the beneficial activity of free-living nematodes based on rapid and cheap soil parameters will be a promising way to increase our ability to monitor the ecological intensification of agricultural practices.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/su13137181/s1. Figure S1. Localization of the fields in Madagascar. Table S1. Raw data of soil properties for each field.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.