Harimenja Razafintsalama
email: harimenjarazafintsalama@gmail.com

Je Trap
email: jean.trap@ird.fr

Bodovololona Rabary
email: bodo.rabary@gmail.com

Adamso Tsoushim

Ernest Razakatiana

Heriniaina Ramanankierana
email: naina.ramanankierana@yahoo.fr

Lilia Rabeharisoa
email: rabeharisoa.lilia@gmail.com

Thierry Becquer
email: thierry.becquer@ird.fr

Effect of Rhizobium Inoculation on Growth of Common Bean in Low-Fertility Tropical Soil Amended with Phosphorus and Lime

Keywords: common bean, rhizobium, phosphorus, lime, mycorrhizal fungi, biological nitrogen fixation, ferralsol, Madagascar

The cultivation of grain legumes (e.g., common bean) in sub-Saharan Africa contributes to the provision of food for a growing population and delivers environmental benefits such as inputs of nitrogen (N) to crops and soil via symbiotic nitrogen fixation (SNF). However, the success of SNF is constrained by several factors such as the poor efficiency of native rhizobial strains to fix N, the low availability of phosphorus (P) and the acidity of soils. Two trials have been conducted in low-fertility tropical soils at the smallholder farm scale in Madagascar to assess the effects of Rhizobium inoculation together with inputs of P and lime on the growth of the common bean. We showed that inoculation with native strains of Rhizobium had significant effects on bean root nodulation, which was increased by up to 15-fold on plant growth, which increased by 78% and on bean yield, which increased by 126%. Moreover, we observed positive and significant relationships between inoculation with Rhizobium and P fertilization on nodulation, plant growth and yield. However, the addition of dolomite lime did not show any effect in our study. The addition of P decreased the mycorrhization rate of roots. Additional research is still needed to improve our understanding of soil fertility conditions (mainly on nutrient availability, including micronutrients) allowing better efficiency of legume symbionts (rhizobium and mycorrhiza) in such low-fertility soils.

Introduction

The cultivation of grain legumes such as common bean, faba bean, cowpea, etc., in sub-Saharan Africa (SSA) contributes both to providing food for a growing population and to delivering environmental benefits to preserve natural resources. In most SSA countries, cropping systems and staple foods are based on cereals. Grain legumes are recognized as rich in proteins and in some essential amino acids and minerals lacking in cereals [START_REF] Mudryj | Nutritional and health benefits of pulses[END_REF]. Therefore, legumes are an excellent complement and have high potential for the improvement of the nutritional quality of foods [START_REF] Vaz Patto | Achievements and challenges in improving the nutritional quality of food legumes[END_REF]. The common bean also offers ecosystem benefits such as inputs of nitrogen (N) to crops and soil via symbiotic nitrogen fixation (SNF) and assists with a decrease in climate change as N fertilization is responsible for half of all agricultural greenhouse gas emissions [START_REF] Shcherbak | Global metaanalysis of the nonlinear response of soil nitrous oxide (N 2 O) emissions to fertilizer nitrogen[END_REF][START_REF] Jensen | Faba bean in cropping systems[END_REF][START_REF] Köpke | Ecological services of faba bean[END_REF]. Despite these advantages, legume cultivation has declined worldwide in the last 50 years relative to cereals, with the notable exception of soybean production [START_REF] Rubiales | Introduction: Legumes in sustainable agriculture[END_REF]. The low productivity of pulses is a general problem for most farming systems in SSA, leading to a decline in soil fertility and a reduction in N fixation related to biological and environmental factors [START_REF] Chianu | Mineral fertilizers in the farming systems of sub-Saharan Africa. A review[END_REF].

The ability of grain legumes to use N 2 through symbiotic nitrogen fixation (SNF) is constrained by several environmental factors. Low phosphorus (P) availability is considered to be the principal limiting factor for legume growth as it strongly reduces N fixation [START_REF] Augusto | Global assessment of limitation to symbiotic nitrogen fixation by phosphorus availability in terrestrial ecosystems using a meta-analysis approach[END_REF]. P plays a key role in the energy metabolism of N fixation [START_REF] Olivera | Growth, nitrogen fixation and ammonium assimilation in common bean (Phaseolus vulgaris): Effect of phosphorus[END_REF], since this process uses high amounts of adenosine triphosphate. When P availability is limited, as it is often the case in tropical soils, symbioses between crop plants and arbuscular mycorrhizal fungi (AMF) can contribute to plant P nutrition [START_REF] Larimer | The interactive effects of plant microbial symbionts: A review and meta-analysis[END_REF][START_REF] Zhang | Soil phosphorus availability determines the preference for direct or mycorrhizal phosphorus uptake pathway in maize[END_REF]. SNF is also strongly limited by soil acidity and aluminum (Al) toxicity. Highly weathered tropical and subtropical soils are generally P deficient and subject to aluminum toxicity [START_REF] Graham | Nitrogen fixation in perspective: An overview of research and extension needs[END_REF][START_REF] Von Uexküll | Global extent, development and economic impact of acid soils[END_REF]. From an agronomic perspective, the former problem can be alleviated by P fertilization and the latter by liming [START_REF] Ndakidemi | Micronutrient uptake in common bean (Phaseolus vulgaris L.) as affected by Rhizobium inoculation, and the supply of molybdenum and lime[END_REF].

Beside P deficiency and Al toxicity, the lack of efficient native strains of Rhizobium for nitrogen fixation could also limit SNF, since less efficient strains often out-compete the more efficient ones that only occupy a small proportion of nodules [START_REF] Martinez-Romero | Diversity of Rhizobium-Phaseolus vulgaris symbiosis: Overview and perspectives[END_REF]. The beneficial effect of inoculation with efficient strains of Rhizobia on legume production has already been mentioned [START_REF] Goedert | Management of the Cerrado soils of Brazil: A review[END_REF][START_REF] Vanlauwe | The role of legumes in the sustainable intensification of African smallholder agriculture: Lessons learnt and challenges for the future[END_REF]. However, the success of efficient strain inoculation on SNF remains driven by the ability of scientists to select inoculants highly efficient for a specific soil within particular environmental conditions and adapted to bean genotypes. Thus, the inoculation of bean seeds with effective and competitive native Rhizobium strains and the supply of appropriate amounts of P and lime could markedly increase SNF and legume growth.

The aim of this study was to assess the effects of inoculation of bean seeds with suitable Rhizobium spp. together with inputs of P and dolomite lime on root nodulation and production of the common bean (Phaseolus vulgaris L.), growing on acid and P-depleted tropical soils. We hypothesized that (i) inoculation with efficient Rhizobium strains, well adapted to common bean, and (ii) improved soil conditions, more specifically a better soil P availability and a lower soil acidity, were needed to improve common bean growth and productivity. Experiments were conducted in two trials in field conditions of resource-poor, smallholder farms in the highlands of Madagascar.

Materials and Methods

Site Location and Description

The study was carried out in Ankazomiriotra, Vakinankaratra region, in the midwest of the Malagasy Highlands. The region has a humid subtropical climate, with a rainy season from October to April. Two sites belonging to different farmers were selected: Site 1 (S: 19 C. The soils of the two sites were classified as Ferralic Cambisol [START_REF]International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports 106[END_REF].

Experimental Design

The field experiments were conducted during the second half of the rainy season, from February to May. The experiment design (Figure 1) included three treatments: (i) input of mineral P as triple super phosphate (TSP) with three application rates (0, 50 and 200 kg P 2 O 5 ha -1 , respectively, designated as P0, P1 and P2); (ii) input of dolomite lime with two application rates (0 or 500 kg ha -1 , respectively, designated as L0 and L1) and (iii) Rhizobium inoculation (absent (I0) or present (I1)). Thus, for each site, we had a 3 × 2 × 2 factorial scheme replicated in 4 blocks, corresponding to 48 plots. Each plot measured 4 m × 2 m and was interspaced by 50 cm.

sowing. The seeds were inoculated with an inoculum provided by the Centre National d Recherches sur l'Environnement (CNRE), Laboratoire de Microbiologie de l'Environne ment (LME), Antananarivo, Madagascar. It corresponded to a cocktail of ten Rhizobium spp. strains selected from bioassays carried out on a collection of infective isolates for thei higher symbiotic effectiveness [START_REF] Razakatiana | Benefits of dual inoculation with arbuscular mycorrhizal fungi and rhizobia on Phaseolus vulgaris planted in a low-fertility tropical soil[END_REF]. The strains are native strains from the highlands o Madagascar, collected in a rice/common bean intercropping trial conducted on a Ferralso located in Lazaina (18° 46′ S; 47° 32′ E), near Antananarivo.

Plant and Soil Sampling

Before the start of the trials, bulk soil sampling was conducted on the 0-10 cm so layer to characterize some physico-chemical properties and determine the mineralogy o the soils. The soil samples were air dried, sieved at 2 mm for the determination of soil contents and pH and ground at 0.2 mm for the determination of total carbon and nitroge contents and soil mineralogical content. Common bean (Phaseolus vulgaris L.) was planted in the fields. We used a local cultivar, called "RI-5-2", resulting from the crossing between two varieties, Ranjonomby and Ikinimba, which has a complete growth cycle of 80 days. Seeds were planted singly at a density of approximately 100,000 plants per ha. The space between rows was 50 cm and the interplant spacing was 20 cm.

Phosphorus and lime were applied after ploughing along the sowing line just before sowing. The seeds were inoculated with an inoculum provided by the Centre National de Recherches sur l'Environnement (CNRE), Laboratoire de Microbiologie de l'Environnement (LME), Antananarivo, Madagascar. It corresponded to a cocktail of ten Rhizobium spp. strains selected from bioassays carried out on a collection of infective isolates for their higher symbiotic effectiveness [START_REF] Razakatiana | Benefits of dual inoculation with arbuscular mycorrhizal fungi and rhizobia on Phaseolus vulgaris planted in a low-fertility tropical soil[END_REF]. The strains are native strains from the highlands of Madagascar, collected in a rice/common bean intercropping trial conducted on a Ferralsol located in Lazaina (18 • 46 S; 47 • 32 E), near Antananarivo.

Plant and Soil Sampling

Before the start of the trials, bulk soil sampling was conducted on the 0-10 cm soil layer to characterize some physico-chemical properties and determine the mineralogy of the soils. The soil samples were air dried, sieved at 2 mm for the determination of soil P contents and pH and ground at 0.2 mm for the determination of total carbon and nitrogen contents and soil mineralogical content.

During field trials, samplings were conducted at two sampling dates for each experiment: at the flowering stage of the common bean (45 days after sowing), in which plants and soils were collected, and at the maturity of the plant (80 days after sowing), for the evaluation of the bean yield.

During the first sampling (at flowering), sixteen plants were randomly selected from a subplot in the middle of each plot. The plants were carefully dug out with their entire root system. The rhizospheric soil, defined as the soil attached to roots after gentle shaking by hand, was carefully separated from the plants and mixed together to form a composite sample. A subsample was frozen for the determination of mineral N content and another subsample was air dried and sieved to 2 mm for other soil analysis. From eight of the sixteen plants, the nodules were separated from the roots and counted to determine the nodule number. Shoots and roots were separated and washed to remove the remaining soil. Plant subsamples (shoot, root and nodule) were oven dried at 60 • C for 48 h and weighed to determine the shoot and root biomasses and the nodule weight. Another four plants with their root systems were washed and oven dried as previously and ground into a fine powder to assess N and P plant contents and SNF. Fresh roots from the remaining four plants were kept in the fridge for mycorrhization-rate determination.

Two other plant species (Urena lobata and Richardia scabra), common on the 2 sites, were also collected as reference plant for the assessment of the SNF. The samples' treatment was the same as for the common bean samples.

At maturity, the bean pods were harvested by hand from ten plants and sun dried for three days before threshing and weighing.

Plant and Soil Analysis

Soil mineralogical contents (i.e., kaolinite, gibbsite and iron oxides (Fe 2 O 3 )) of the bulk soil samples were estimated by near infrared reflectance spectroscopy (NIRS), using the models developed by Ramaroson et al. [START_REF] Ramaroson | Mineralogical analysis of ferralitic soils in Madagascar using NIR spectroscopy[END_REF]. The total carbon and nitrogen contents were determined by dry combustion in a Flash 2000 CHN Analyzer (Thermo Fisher Scientific GmbH., Dreieich, Germany). The remaining phosphorus (Prem), defined as the P concentration that remains in solution after shaking a soil sample for a certain period with a solution with known initial P concentration, was extracted according to the methodology of Alvarez et al. [START_REF] Alvarez | Determinação e uso do fósforo remanescente[END_REF]. In brief, P content was determined at equilibrium after soil shaking in a solution of 0.01 mol L -1 CaCl 2 containing 60 mg L -1 P, using a 1:10 soil-to-solution ratio for 1 h. Available P (Resin P) was extracted with anion exchange membrane (BDH no. 55164, 6 cm × 2 cm) in bicarbonate form [START_REF] Kuo | In Methods of Soil Analysis, Part 3. Chemical Methods[END_REF]. Inorganic P was determined spectrophotometrically for Prem and Resin P, using the molybdenum blue-ascorbic acid method [START_REF] Murphy | A modified single solution method for the determination of phosphate in natural waters[END_REF]. The soil pH was measured in water and in 1 M KCl suspensions using 1:5 (w:v) soil-to-solution ratio. The main mineralogical and physico-chemical soil properties of the two sites are shown in Table 1. The mineral N content of the rhizospheric soil was determined by shaking 5 g of fresh soil in 50 mL 1 M KCl for one hour. After centrifugation and filtration, the solution was analyzed by colorimetry using a segmented flow analyzer (SAN++ SA3000/5000 Skalar Analytical, Breda, The Netherlands). The Resin P content and the soil pH water were determined as described previously. Total microbial activity in soil was estimated by measuring fluorescein diacetate (FDA) hydrolysis according to the method of Schnürer and Rosswall [START_REF] Schnürer | Fluorescein diacetate hydrolysis as a measure of total microbial activity in soil and litter[END_REF]. The acid phosphatase activity was measured by the method of Tabatabai [START_REF] Tabatabai | Methods of Soil Analyses Part 2, Microbiological and Biochemical Properties[END_REF] using hydrolysis of p-nitrophenylphosphate tetrahydrate (p-NPP) buffered at 6.0. The total N content of a plant was determined by dry combustion in a Flash 2000 CHN Analyzer (Thermo Fisher Scientific GmbH, Dreieich, Germany). The P content was determined after calcination at 550 • C. After cooling, the ash was dissolved in warm 2% HCl. Total P was determined colorimetrically using the molybdenum blue-ascorbic acid method. The SNF was assessed by the natural abundance method. The 15 N natural abundance (δ 15 N) of plant samples was measured with a mass spectrometer (Isoprime Precision, Elementar, Langenselbold, Germany) coupled with a CN elemental analyzer (Vario Pyro Cube, Elementar). The percentage of nitrogen derived by the atmosphere (%Ndfa) was calculated as follows [START_REF] Samago | Grain yield of common bean (Phaseolus vulgaris L.) varieties is markedly increased by rhizobial inoculation and phosphorus application in Ethiopia[END_REF]: %Ndfa = ((δ 15 Nrefδ 15 Nleg)/(δ 15 Nref -B)) × 100 where δ 15 Nleg corresponds to the legume δ 15 N, δ 15 Nref is the nonfixing reference plant's δ 15 N and the B value is the δ 15 N of the legume grown in a substrate without N. The δ 15 Nleg ref was the mean of two reference plants δ 15 N, Richardia scabra and Urena lobata. The δ 15 N for the reference plants were 6.62 and 6.16‰, respectively, for Site 1 and Site 2. The B-value was -0.482‰ [START_REF] Samago | Grain yield of common bean (Phaseolus vulgaris L.) varieties is markedly increased by rhizobial inoculation and phosphorus application in Ethiopia[END_REF]. The amount of N fixed by the crop was deduced by multiplying % Ndfa by the shoot N content and shoot biomass.

To determine the mycorrhization rates from the roots, fresh roots were stained with Trypan Blue [START_REF] Phillips | Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection[END_REF] and the percentage of root length colonized by the mycorrhizal fungus was quantified using 15 × 1 cm root segments by the method of Giovannetti and Mosse [START_REF] Giovannetti | An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots[END_REF].

Statistical Analysis

Statistical significance was set at p-value. All statistical tests were performed with the R software. Mean and standard deviation values were computed per treatment for all variables. A principal component analysis (PCA) was conducted on plant and soil variables using the package "FactoMineR". Then, we performed a three-way ANOVA on the first two PCA axes as response variables (explaining 46% of total inertia) and with rhizobial inoculation, liming and P supply as factors. The normality and homogeneity of variance of the residuals were tested using the Shapiro's test and Levene's test (packages "stats" and "cars"). Because it was not significant, we did not include the effect of site in the ANOVA. We calculated the effect size for ANOVA using the eta-squared from the "lsr" package.

Results

Legume Symbioses

The inoculation with rhizobia and the supply of P had significant effects on the number of nodules. In the absence of inoculation, P and lime bean nodulation was nil (I0P0L0). After inoculation, the mean number of nodules increased from 6 nodules plant -1 for I0 to 91 nodules plant -1 for I1 (Figure 2a). However, even with inoculation, the number of nodules remained low without P (39 nodules plant -1 ), and it increased sharply with P to reach 116 to 118 nodules plant -1 for P1 and P2, respectively, with no significant difference between P1 and P2 (Figure 2b). A positive interaction was also observed between P and dolomite lime inputs, leading to a slight increase in the number of nodules per plant with dolomite lime (L1) (100 nodules plant -1 for I1L1) compared with L0 (82 nodules plant -1 for I1L0) but without significant differences between treatments with or without lime (see ANOVA and comparisons of means for all treatments on the data repository).

The nodule weight varied similarly to the number of nodules (Figure 2c,d). Nevertheless, the positive effect of the application of P was reinforced due to the increase in the average weight of the nodules with additional P, with an average weight of 0.90 mg per nodule for I0P0 and 1.2 mg per nodule for I0P1 and I0P2. Lime had no effect on nodule weight. I1L0) but without significant differences between treatments with or without lime (see ANOVA and comparisons of means for all treatments on the data repository).

The nodule weight varied similarly to the number of nodules (Figure 2c,d). Nevertheless, the positive effect of the application of P was reinforced due to the increase in the average weight of the nodules with additional P, with an average weight of 0.90 mg per nodule for I0P0 and 1.2 mg per nodule for I0P1 and I0P2. Lime had no effect on nodule weight. The rate of mycorrhization was not impacted by inoculation with rhizobia (see ANOVA and comparisons of means for all treatments on the data repository). However, an interaction was observed between P and dolomite inputs (Figure 3). In the absence of dolomite, the addition of P seemed to reduce mycorrhization, which decreased from 37% to 21%. With dolomite, the intermediate dose of P (P1) seemed to promote mycorrhization, with a mycorrhization rate of 43%.

Legume Growth and Biological N2 Fixation

The sole significant effect was the interaction of inoculation and P supply, which had positive effects on biomass production at flowering and on bean yield (Figure 4a,b). The The rate of mycorrhization was not impacted by inoculation with rhizobia (see ANOVA and comparisons of means for all treatments on the data repository). However, an interaction was observed between P and dolomite inputs (Figure 3). In the absence of dolomite, the addition of P seemed to reduce mycorrhization, which decreased from 37% to 21%. With dolomite, the intermediate dose of P (P1) seemed to promote mycorrhization, with a mycorrhization rate of 43%. The rate of mycorrhization was not impacted by inoculation w ANOVA and comparisons of means for all treatments on the data repo an interaction was observed between P and dolomite inputs (Figure 3). dolomite, the addition of P seemed to reduce mycorrhization, which de to 21%. With dolomite, the intermediate dose of P (P1) seemed to promot with a mycorrhization rate of 43%.

Legume Growth and Biological N2 Fixation

The sole significant effect was the interaction of inoculation and P s positive effects on biomass production at flowering and on bean yield (

Legume Growth and Biological N 2 Fixation

The sole significant effect was the interaction of inoculation and P supply, which had positive effects on biomass production at flowering and on bean yield (Figure 4a,b). The biomass was multiplied by a factor of 2.8 from I0P0 to I1P2, from 0.92 g plant -1 to 2.61 g plant -1 . The bean yield increased by a factor of 3.5, from 82 to 287 kg ha -1 .

biomass was multiplied by a factor of 2.8 from I0P0 to I1P2, from 0.92 g plant -1 . The bean yield increased by a factor of 3.5, from 82 to 287 kg ha -1 . The inoculation significantly increased the N content of the biomass 30 to 35 g kg -1 (Figure 5a). However, the lime decreased the N content of th had a significant effect on plant P content, which almost doubled between 2.7 to 5 g kg -1 (Figure 5b). Inoculation and interactions between inoculation, P and lime inputs effects on the %Ndfa and the amount of N fixed by the plant (Figure 6). W The inoculation significantly increased the N content of the biomass (by 17%), from 30 to 35 g kg -1 (Figure 5a). However, the lime decreased the N content of the plant. P input had a significant effect on plant P content, which almost doubled between P0 and P2, from 2.7 to 5 g kg -1 (Figure 5b). The inoculation significantly increased the N content of the biomass (by 17%), from 30 to 35 g kg -1 (Figure 5a). However, the lime decreased the N content of the plant. P input had a significant effect on plant P content, which almost doubled between P0 and P2, from 2.7 to 5 g kg -1 (Figure 5b). Inoculation and interactions between inoculation, P and lime inputs had significant effects on the %Ndfa and the amount of N fixed by the plant (Figure 6). Without inoculation, the %Ndfa was low (%Ndfa < 30% for I0). Inoculation without P and lime inputs had no effect (%Ndfa = 35% for I1P0L0). After inoculation, the inputs of lime and P had positive effects, the %Ndfa reaching 68% for I1P0L1, 80% for I1P1L0 and I1P2L0 and up to 89% for I1P1L1 and I1P2L1. The amount of fixed N, which was negligible in the absence of inoculation (<1 kg ha -1 ), reached 9.5 kg ha -1 for I1P2L1. Inoculation and interactions between inoculation, P and lime inputs had significant effects on the %Ndfa and the amount of N fixed by the plant (Figure 6). Without inoculation, the %Ndfa was low (%Ndfa < 30% for I0). Inoculation without P and lime inputs had no effect (%Ndfa = 35% for I1P0L0). After inoculation, the inputs of lime and P had positive effects, the %Ndfa reaching 68% for I1P0L1, 80% for I1P1L0 and I1P2L0 and up to 89% for I1P1L1 and I1P2L1. The amount of fixed N, which was negligible in the absence of inoculation (<1 kg ha -1 ), reached 9.5 kg ha -1 for I1P2L1.

Rhizospheric Soil

Inoculation had no effect on soil properties. The main effects of the treatments were the effects of P and lime on Resin P and pH. The Resin P was very low without P addition (3.7 mg P kg -1 , on average) and increased rapidly after the addition of P1 (96 mg P kg -1 ) and P2 (308 mg P kg -1 ). However, an interaction between P and lime was observed, resulting in a decrease in Resin P content at P0, after lime input (see ANOVA and comparison of means for all treatments on the data repository). Lime also had a noticeable effect on soil pH, which increased from 5 to 6 (Figure 7b). The treatments did not have marked effects either on the soil mineral nitrogen contents (N-NO3 -and N-NH4 + ) or on the soil biological activities (total microbial activity and phosphatase acid activity). All treatments combined, the average nitrate and ammonium contents of the rhizospheric soil were, respectively, 15.6 mg N-NO3 -kg -1 and 7.88 N-NH4 + kg -1 soil. The total microbial activity of the rhizospheric soil was estimated at 161 μg fluorescein g -1 h -1 and the acid phosphatase activity at 191 μg p-NP g -1 h -1 .

Effects of Inoculation of Bean Seed and Inputs of P and Dolomite on Common Bean and Soil

The factorial map of the PCA (Figure 8a) showed a clear distinction between the noninoculated treatments, positioned on the positive scores of the first axis, and the inoculated treatments with the negative scores. However, only the inoculated treatments that received P had negative scores, whereas the treatments inoculated without supply of P

Rhizospheric Soil

Inoculation had no effect on soil properties. The main effects of the treatments were the effects of P and lime on Resin P and pH. The Resin P was very low without P addition (3.7 mg P kg -1 , on average) and increased rapidly after the addition of P1 (96 mg P kg -1 ) and P2 (308 mg P kg -1 ). However, an interaction between P and lime was observed, resulting in a decrease in Resin P content at P0, after lime input (see ANOVA and comparison of means for all treatments on the data repository). Lime also had a noticeable effect on soil pH, which increased from 5 to 6 (Figure 7b). The treatments did not have marked effects either on the soil mineral nitrogen contents (N-NO 3 -and N-NH 4 + ) or on the soil biological activities (total microbial activity and phosphatase acid activity). All treatments combined, the average nitrate and ammonium contents of the rhizospheric soil were, respectively, 15.6 mg N-NO 3 -kg -1 and 7.88 N-NH 4 + kg -1 soil. The total microbial activity of the rhizospheric soil was estimated at 161 µg fluorescein g -1 h -1 and the acid phosphatase activity at 191 µg p-NP g -1 h -1 .

Sustainability 2022, 14, x FOR PEER REVIEW 8 of 14

Rhizospheric Soil

Inoculation had no effect on soil properties. The main effects of the treatments were the effects of P and lime on Resin P and pH. The Resin P was very low without P addition (3.7 mg P kg -1 , on average) and increased rapidly after the addition of P1 (96 mg P kg -1 ) and P2 (308 mg P kg -1 ). However, an interaction between P and lime was observed, resulting in a decrease in Resin P content at P0, after lime input (see ANOVA and comparison of means for all treatments on the data repository). Lime also had a noticeable effect on soil pH, which increased from 5 to 6 (Figure 7b). The treatments did not have marked effects either on the soil mineral nitrogen contents (N-NO3 -and N-NH4 + ) or on the soil biological activities (total microbial activity and phosphatase acid activity). All treatments combined, the average nitrate and ammonium contents of the rhizospheric soil were, respectively, 15.6 mg N-NO3 -kg -1 and 7.88 N-NH4 + kg -1 soil. The total microbial activity of the rhizospheric soil was estimated at 161 μg fluorescein g -1 h -1 and the acid phosphatase activity at 191 μg p-NP g -1 h -1 .

Effects of Inoculation of Bean Seed and Inputs of P and Dolomite on Common Bean and Soil

The factorial map of the PCA (Figure 8a) showed a clear distinction between the noninoculated treatments, positioned on the positive scores of the first axis, and the inoculated treatments with the negative scores. However, only the inoculated treatments that received P had negative scores, whereas the treatments inoculated without supply of P

Effects of Inoculation of Bean Seed and Inputs of P and Dolomite on Common Bean and Soil

The factorial map of the PCA (Figure 8a) showed a clear distinction between the noninoculated treatments, positioned on the positive scores of the first axis, and the inoculated treatments with the negative scores. However, only the inoculated treatments that received P had negative scores, whereas the treatments inoculated without supply of P had positive scores similar to the noninoculated treatments. The second axis of the PCA opposed the treatments following the level of P (the positive scores corresponding to the absence of P supply) and, to a lesser extent, according to the supply of lime (the negative scores corresponding to the absence of lime supply). The individual effects of the treatments were highlighted in Figure 8c-e and by the statistical results of the ANOVA carried out on the first two axes (Figure 8f). had positive scores similar to the noninoculated treatments. The second axis of the PCA opposed the treatments following the level of P (the positive scores corresponding to the absence of P supply) and, to a lesser extent, according to the supply of lime (the negative scores corresponding to the absence of lime supply). The individual effects of the treatments were highlighted in Figure 8c-e and by the statistical results of the ANOVA carried out on the first two axes (Figure 8f).

Treatments inoculated and fertilized with P, associated with the negative scores of the first axis, corresponded to an increase in nodulation in %Ndfa as well as an increase in biomass at flowering and in bean yield (Figure 8b). On the second axis, negative values corresponded to an increase in soil P availability and an increase in plant P content. Positive values corresponded to an increase in soil pH. Treatments inoculated and fertilized with P, associated with the negative scores of the first axis, corresponded to an increase in nodulation in %Ndfa as well as an increase in biomass at flowering and in bean yield (Figure 8b). On the second axis, negative values corresponded to an increase in soil P availability and an increase in plant P content. Positive values corresponded to an increase in soil pH.

Discussion

Our results indicate that inoculation with native strains of Rhizobium enhanced bean nodule production, biomass production and bean yield (Figures 2, 3 and8). Many studies have highlighted the positive effects of rhizobial symbioses on legumes [START_REF] Díaz-Ambrona | Cereal-legume rotations in a Mediterranean environment: Biomass and yield production[END_REF][START_REF] Lupwayi | Grain legumes in northern Great Plains: Impacts on selected biological soil processes[END_REF][START_REF] Zhong | Nitrous oxide emissions associated with nitrogen fixation by grain legumes[END_REF]. A metaanalysis showed an increase in bean yield of 59% in pot experiments and 16% in the field due to rhizobial inoculation [START_REF] Kaschuk | Three decades of soil microbial biomass studies in Brazilian ecosystems: Lessons learned about soil quality and indications for improving sustainability[END_REF]. Here, we showed an increase of 126% for the bean yield with inoculation in the field. The use of native rhizobia to increase SNF has been shown to be beneficial to the common bean [START_REF] Ouma | Elucidating the potential of native rhizobial isolates to improve biological nitrogen fixation and growth of common bean and soybean in smallholder farming systems of Kenya[END_REF][START_REF] Koskey | Potential of native rhizobia in enhancing nitrogen fixation and yields of climbing beans (Phaseolus vulgaris L.) in contrasting environments of Eastern Kenya[END_REF], soybean [START_REF] Ouma | Elucidating the potential of native rhizobial isolates to improve biological nitrogen fixation and growth of common bean and soybean in smallholder farming systems of Kenya[END_REF], or cowpea [START_REF] Mathu | Potential of indigenous bradyrhizobia versus commercial inoculants to improve cowpea (Vigna unguiculata L. walp.) and green gram (Vigna radiata L. wilczek.) yields in Kenya[END_REF]. According to Ouma et al. [START_REF] Ouma | Elucidating the potential of native rhizobial isolates to improve biological nitrogen fixation and growth of common bean and soybean in smallholder farming systems of Kenya[END_REF], it is preferable to use native rhizobia because they are better suited to local pedoclimatic conditions. The common bean is a relatively permissive host, nodulating with many rhizobial species [START_REF] Michiels | Phaseolus vulgaris is a non-selective host for nodulation[END_REF]. The rhizobial strains we used were isolated in a companion study aiming to describe the diversity of the community of rhizobia associated with common beans cultivated in the Madagascar Highlands [START_REF] Razakatiana | Rôle des Symbioses Mycorhizienne et Bactérienne: Biodisponibilité du Phosphore et Développement des Plantes de Haricot Sur le Sol Ferralitique de Madagascar[END_REF]. The selected strains were able to infect common bean root (infectivity) efficiently and had good effectiveness, i.e., a capacity to produce a larger amount of bean biomass compared with the noninoculated control [START_REF] Razakatiana | Benefits of dual inoculation with arbuscular mycorrhizal fungi and rhizobia on Phaseolus vulgaris planted in a low-fertility tropical soil[END_REF]. At this time, the efficient rhizobia strains we used are not fully characterized, but they were able to multiply bean biomass by a factor of 2.6 to 3.1 [START_REF] Razakatiana | Rôle des Symbioses Mycorhizienne et Bactérienne: Biodisponibilité du Phosphore et Développement des Plantes de Haricot Sur le Sol Ferralitique de Madagascar[END_REF]. Indeed, as most of Rhizobium native strains of soils are poorly infective or not efficient [START_REF] Zhong | Nitrous oxide emissions associated with nitrogen fixation by grain legumes[END_REF][START_REF] Catroux | Trends in rhizobial inoculant production and use[END_REF][START_REF] Hynes | Rhizobium population dynamics in the pea rhizosphere of rhizobial inoculant strain applied in different formulations[END_REF], there is a huge need to select efficient strains which are well adapted to each different soil condition.

The effectiveness of inoculation was very limited in the absence of P input. However, inoculation with rhizobia and fertilization with P had a highly significant effect on common bean nodulation and productivity (Figures 2 and3). While mean nodulation was only 8.7 ± 2.1 mg plant -1 in the absence of inoculation and 17.0 ± 3.1 mg plant -1 in the absence of P, it increased to 84.1 ± 13.8 mg plant -1 after inoculation and fertilization (Figure 2c,d). The increase in P availability (Resin P) (Figure 7a) and shoot P content (Figure 5b) suggests that SNF was strongly limited by the low P availability. It is widely accepted that the low availability of P in terrestrial ecosystems limits SNF [START_REF] Augusto | Global assessment of limitation to symbiotic nitrogen fixation by phosphorus availability in terrestrial ecosystems using a meta-analysis approach[END_REF]. The low content of available P in soils in Madagascar is one of the main constraints on crop production [START_REF] Raminoarison | Multiple-nutrient limitation of upland rainfed rice in ferralsols: A greenhouse nutrient-omission trial[END_REF]. A survey of farm fields showed that 100% of rice flag leaf P content was below the critical value (2.4 g P kg -1 leaf biomass) in rainfed rice in highly weathered soils of Madagascar [START_REF] Rabeharisoa | Larger bioavailability of soil phosphorus for irrigated rice compared with rainfed rice in Madagascar: Results from a soil and plant survey[END_REF]. In these soils, without P input, the %Ndfa was very low ( ∼ =10%) and the amount of fixed N < 1 kg ha -1 , whereas the %Ndfa was up to 95% and the amount of fixed N increased up to 10 kg ha -1 for I1P2L1. Indeed, legumes depending on symbiotic N 2 fixation have a higher internal P requirement than plants dependent on combined N [START_REF] Ribet | The phosphorus requirement of N 2 -fixing and urea-fed Acacia mangium[END_REF]. The conversion of N 2 into ammonium by nitrogenase requires a large amount of energy in the form of ATP [START_REF] Amâncio | Nitrogen Acquisition and Assimilation in Higher Plants; Plant ecophysiology[END_REF][START_REF] Neila | Phaseolus vulgaris-Rhizobia symbiosis increases the phosphorus uptake and symbiotic N 2 fixation under insoluble phosphorus[END_REF][START_REF] Zoundji | Inoculation and phosphorus application effects on soybean [Glycine max (L.) Merrill] productivity grown in farmers' fields of Benin[END_REF]. Then, P deficiency has a negative effect on the energy metabolism of legume nodules [START_REF] Sa | Energy status and functioning of phosphorus-deficient soybean nodules[END_REF]. However, P deficiency is not the only nutrient deficiency that can impact nodulation effectiveness. Other nutrients, such as Ca, are very deficient in Malagasy soils [START_REF] Raminoarison | Multiple-nutrient limitation of upland rainfed rice in ferralsols: A greenhouse nutrient-omission trial[END_REF]. Micronutrients are also often deficient in tropical soils. The supply of highquality organic matter could therefore have a beneficial effect on nodulation through an increase in ecological processes involved in the supply of bioavailable P to plants and soil organisms [START_REF] Trap | Effects of the earthworm Pontoscolex corethrurus on rice P nutrition and plant-available soil P in a tropical Ferralsol[END_REF][START_REF] Trap | Agricultural practices modulate the beneficial activity of bacterial-feeding nematodes on plant growth and nutrition: Evidence from an original intact soil core technique[END_REF].

Liming had no marked effect on bean nodulation or biomass. In fact, liming can have contrasting effects on the soil or living organisms (plants, microorganisms). On the one hand, the positive effects of liming on the fertility of tropical soils are well known [START_REF] Fageria | Ameliorating soil acidity of tropical Oxisols by liming for sustainable crop production[END_REF]. Benefits are also provided to soil microorganisms; liming of acid soils provides nutrients (Ca and Mg) and improves soil conditions for bacterial growth by increasing soil pH and availability of P and molybdenum [START_REF] Andrade | The diversity of Phaseolus-nodulating rhizobial populations is altered by liming of acid soils planted with Phaseolus vulgaris L. in Brazil[END_REF]. At pH 4.9 to 5.1 (Figure 7b), due to the close relationship between pH and exchangeable aluminum content, toxic effects of aluminum to biota are expected. Unfortunately, the exchangeable aluminum was not analyzed in our study. Nevertheless, at pH 4.7 in the bulk soil (Table 1) or pH 5.0 in the rhizospheric soil (Figure 7b), the conditions are highly favorable to the presence of exchangeable Al in tropical soils [START_REF] Vendrame | Acidity control in Latosols under long-term pastures in the Cerrado region, Brazil[END_REF] and therefore to the solubilization of Al in toxic forms in the soil solution [START_REF] Boudot | Aluminium toxicity in declining forests: A general overview with a seasonal assessment in the Vosges mountains (France)[END_REF]. The rhizobial population can decline rapidly in soils with high aluminum saturation [START_REF] Andrade | The diversity of Phaseolus-nodulating rhizobial populations is altered by liming of acid soils planted with Phaseolus vulgaris L. in Brazil[END_REF]. On the other hand, the increase of soil pH after liming influences soil Pi sorption, leading to the increase in soil P availability [START_REF] Barrow | The effects of pH on phosphate uptake from the soil[END_REF]. However, in our study, P availability (Resin P) was not increased by liming (Figure 7a). Indeed, the key role often attributed to pH on P availability can be questioned. In two recent studies [START_REF] Achat | Soil properties controlling inorganic phosphorus availability: General results from a national forest network and a global compilation of the literature[END_REF][START_REF] Khaledian | Developing global pedotransfer functions to estimate available soil phosphorus[END_REF] focusing on soil properties that control P availability, pH did not appear to be a major explanatory property. Moreover, phosphate uptake by roots increased when pH decreased [START_REF] Barrow | The effects of pH on phosphate uptake from the soil[END_REF], leading to a higher plant P content in the most acidic conditions (Figure 8b). Thus, the alleviation of the toxic effects of aluminum at pH 6 (Figure 7b) could be counteracted by a decrease in plant P uptake.

Interactions between rhizobia and arbuscular mycorrhizal fungi (AMF) can also be expected. AMF are symbionts that provide various compounds to plants, mainly P, in exchange for plant photosynthates [START_REF] Smith | Mycorrhizal Symbiosis[END_REF]. Arbuscular mycorrhiza can enhance P uptake due to colonization by hyphae of soil volumes poorly accessible to roots, i.e., through a better mycorrhization rate and through mineralization of organic P by a higher phosphatase activity [START_REF] Joner | Phosphatase activity of extra-radical arbuscular mycorrhizal hyphae: A review[END_REF][START_REF] Hinsinger | Bioavailability of soil inorganic P in the rhizosphere as affected byroot-induced chemical changes: A review[END_REF]. Our results did not show any clear effect of inoculation with Rhizobium on the mycorrhization rate (Figure 3). However, high soil available P was associated to a lower mycorrhization rate. Indeed, high P levels are known to affect AMF symbiosis negatively [START_REF] Balzergue | The regulation of arbuscular mycorrhizal symbiosis by phosphate in pea involves early and systemic signalling events[END_REF]. Anyway, to optimize the P plant uptake by the plant, it could be necessary to inoculate with both AMF and Rhizobium. The benefits of co-inoculation are reported by several authors [START_REF] Larimer | The interactive effects of plant microbial symbionts: A review and meta-analysis[END_REF][START_REF] Aryal | Rhizobia and AM fungal inoculation improve growth and nutrient uptake of bean plants under organic fertilization[END_REF]. Nevertheless, to optimize the beneficial effect of the co-inoculation, and increase the supply of P to rhizobia by mycorrhizae, it is necessary to manage the phosphate fertilization in an adequate way. Indeed, too-low or too-high availability of P both limit mycorrhization and thus promote plant nutrition by the direct route (roots) instead of the indirect route (mycorrhiza) [START_REF] Zhang | Soil phosphorus availability determines the preference for direct or mycorrhizal phosphorus uptake pathway in maize[END_REF].

Conclusions

The low crop productivity in grain legumes in sub-Saharan Africa limits the increase in their production in most farming systems. Low productivity is due to biological and environmental factors, i.e., the poor efficiency of native rhizobial strains, the low availability of P and the acidity of soils. Our results show that rhizobial inoculation with native strains enhanced bean root nodulation, and as a consequence, plant productivity through a better FSN. In addition, we showed a positive significant relationship between rhizobial inoculation and P input. However, we showed an antagonistic effect of P input on the rate of mycorrhization. We did not show a positive effect of the supply of dolomite lime in our study, the positive effects of Ca and Mg inputs being able to be offset by negative effects on the absorption of P by the plants. According to our results, we suggest that relatively low rates of P (50 kg P 2 O 5 ha -1 ) and lime (200-300 kg dolomite ha -1 , i.e., below the tested level L1 that led to a high soil pH increase), applied in the fertilizer banding line or in the sowing hole, are sufficient to increase P availability and reduce Al toxicity for the legume and rhizobium, with weak detrimental effects. Nevertheless, it remains for us to investigate further how the co-inoculation in the field, with both rhizobial and mychorrhiza, can contribute to the success of legumes, especially in harsh environments and low-fertility soils.

Figure 1 .

 1 Figure 1. Localization of the field trials in Ankazomiriotra, Vakinankaratra region, Madagascar (a and schematic representation of the trial (b).

Figure 1 .

 1 Figure 1. Localization of the field trials in Ankazomiriotra, Vakinankaratra region, Madagascar (a) and schematic representation of the trial (b).

Figure 2 .

 2 Figure 2. Rhizobial inoculation (I) (a,c) and interaction effect between inoculation (I) and phosphorus (P) (b,d) on nodule number and nodule dry weight of common bean. Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 3 .

 3 Figure 3. Interaction between phosphorus (P) and lime (L) inputs on the rate of mycorrhization. Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 2 .

 2 Figure 2. Rhizobial inoculation (I) (a,c) and interaction effect between inoculation (I) and phosphorus (P) (b,d) on nodule number and nodule dry weight of common bean. Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 2 .

 2 Figure 2. Rhizobial inoculation (I) (a,c) and interaction effect between inoculatio rus (P) (b,d) on nodule number and nodule dry weight of common bean. Diffe significant difference among the treatments according to post hoc Tukey HSD te

Figure 3 .

 3 Figure 3. Interaction between phosphorus (P) and lime (L) inputs on the rate of Different letters indicate significant difference among the treatments according HSD test at p-value < 0.05.

Figure 3 .

 3 Figure 3. Interaction between phosphorus (P) and lime (L) inputs on the rate of mycorrhization. Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 4 .

 4 Figure 4. Interaction between inoculation (I) and phosphorus (P) on the biomass (a common bean (b). Different letters indicate significant difference among the trea to post hoc Tukey HSD test at p-value < 0.05.

Figure 5 .

 5 Figure 5. Interaction between inoculation (I) and lime (L) on the total plant N (a) (P) effect on the total plant P of common bean (b). Different letters indicate sign among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 4 .

 4 Figure 4. Interaction between inoculation (I) and phosphorus (P) on the biomass (a) and the yield of common bean (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Sustainability 2022 ,

 2022 14, x FOR PEER REVIEW 7 of 14 biomass was multiplied by a factor of 2.8 from I0P0 to I1P2, from 0.92 g plant -1 to 2.61 g plant -1 . The bean yield increased by a factor of 3.5, from 82 to 287 kg ha -1 .

Figure 4 .

 4 Figure 4. Interaction between inoculation (I) and phosphorus (P) on the biomass (a) and the yield of common bean (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 5 .

 5 Figure 5. Interaction between inoculation (I) and lime (L) on the total plant N (a) and phosphorus (P) effect on the total plant P of common bean (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 5 .

 5 Figure 5. Interaction between inoculation (I) and lime (L) on the total plant N (a) and phosphorus (P) effect on the total plant P of common bean (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 6 .

 6 Figure 6. Inoculation (I), phosphorus (P) and lime (L) interaction on the %Ndfa (a) and the amount of fixed N (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 7 .

 7 Figure 7. Interaction between phosphorus (P) and lime (L) on the soil Resin P content (a) and soil pH (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 6 .

 6 Figure 6. Inoculation (I), phosphorus (P) and lime (L) interaction on the %Ndfa (a) and the amount of fixed N (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 6 .

 6 Figure 6. Inoculation (I), phosphorus (P) and lime (L) interaction on the %Ndfa (a) and the amount of fixed N (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 7 .

 7 Figure 7. Interaction between phosphorus (P) and lime (L) on the soil Resin P content (a) and soil pH (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 7 .

 7 Figure 7. Interaction between phosphorus (P) and lime (L) on the soil Resin P content (a) and soil pH (b). Different letters indicate significant difference among the treatments according to post hoc Tukey HSD test at p-value < 0.05.

Figure 8 .Figure 8 .

 88 Figure 8. Principal component analysis performed on plant and soil variables. (a) Individual factorial map with clustering based on rhizobial inoculation, liming and P treatments. (b) Circle of correlations. (c) Individual clustering based on rhizobial inoculation treatments. (d) Individual clustering based on P treatments. (e) Individual clustering based on liming treatments. (f) ANOVA performed on the first two PCA axes as response variables and with rhizobial inoculation, liming and P treatment as factors. Abbreviations in (b): plt.mass: plant biomass at flowering; yield: bean yield; N.plt: plant N content; P.plt: plant P content; Nod.number: nodule number; Nod.mass: nodule weight; Myc: rate of mycorrhization; pH: pH of the rhizospheric soil; nitrate: nitrate content of the Figure 8. Principal component analysis performed on plant and soil variables. (a) Individual factorial map with clustering based on rhizobial inoculation, liming and P treatments. (b) Circle of correlations. (c) Individual clustering based on rhizobial inoculation treatments. (d) Individual clustering based on P treatments. (e) Individual clustering based on liming treatments. (f) ANOVA performed on the first two PCA axes as response variables and with rhizobial inoculation, liming and P treatment as factors. Abbreviations in (b): plt.mass: plant biomass at flowering; yield: bean yield; N.plt: plant N content; P.plt: plant P content; Nod.number: nodule number; Nod.mass: nodule weight; Myc: rate of mycorrhization; pH: pH of the rhizospheric soil; nitrate: nitrate content of the rhizospheric soil; ammonium: ammonium content of the rhizospheric soil; Res.P: resin P content of the rhizospheric soil; FDA: total microbial activity of the rhizospheric soil; Aci.P: acid phosphatase activity of the rhizospheric soil. Principal component analysis performed on plant and soil variables. (a) Individual factorial map with clustering based on rhizobial inoculation, liming and P treatment. (b) Circle of correlations. (c) Individual clustering based on rhizobial inoculation treatments. (d) Individual clustering based on P treatments. (e) Individual clustering based on liming treatments. (f) ANOVA performed on the first two PCA axes as response variables and with rhizobial inoculation, liming and P treatment as factors.

  • 40 36.1 ; E: 46 • 33 32.4 ; 1171 m asl) and Site 2 (S: 19 • 39 40.6 ; E: 046 • 33 24.2 ; 1181 m asl). The trials were set up in 2015. The temperature and rainfall were measured daily with an automatic weather station (Cimel Electronique S.A.S, Paris, France) installed at Ivory, also in the midwest of the Malagasy Highlands (19 • 33 S, 46 • 24 E, 900 m asl) 20 km from the studied sites. The average annual rainfall and temperature (2007-2015) were, respectively, 1326 mm year -1 and 23

Table 1 .

 1 Physico-chemical properties and mineralogy of the soils of the two sites.

	Site 1	Site 2

Acknowledgments:

The authors are grateful to the two farmers from the Vakinankaratra region which allowed us to set up our trials in their plots. We thank Jean-Jacques Drevon for his leadership of the Fabatropimed project (Ecological services of legumes for nitrogen and phosphorus biogeochemical cycles and C sequestration in cereal cropping systems in Africa and the Mediterranean Basin), which funded this work. We are also thankful to M.P. Razafimanantsoa and A.H.D. Razafimahafaly (LRI, Antananarivo, Madagascar) for their technical assistance in the laboratory. We would also like to thank the DP SPAD (Dispositif en partenariat "Systèmes de Production d'Altitude et Durabilité à Madagascar), an interinstitutional scientific platform on highland agricultural systems and sustainability in Madagascar, for the meteorological data collected through the automatic weather station installed in the site of Ivory.

Funding: This research was supported by the Agropolis Fondation (France) under the reference ID 321 1001-009 through the "Investissements d'Avenir" program (Labex Agro: ANR-10-LABX-322 0001-01).

Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.

Data Availability Statement: Data were deposited on a figshare repository at the following DOI: 10.6084/m9.figshare.19433318 (Licence CC BY 4.0).

Conflicts of Interest:

The authors declare no conflict of interest.