# Trait Diversity of Pulse Species Predicts Agroecosystem Properties Trade-Offs

Julie Guiguitant, Denis Vile, Hélène Marrou

► To cite this version:

Julie Guiguitant, Denis Vile, Hélène Marrou. Trait Diversity of Pulse Species Predicts Agroecosystem Properties Trade-Offs. Frontiers in Plant Science, 2021, 12, pp.636915. 10.3389/fpls.2021.636915. hal-03212645

HAL Id: hal-03212645

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

Submitted on 29 Apr 2021

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

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

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

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

Check for
updates

# Trait Diversity of Pulse Species Predicts Agroecosystem Properties Trade-Offs

Julie Guiguitant¹,², Denis Vile²* and Hélène Marrou³

¹ SYSTEM, Montpellier SupAgro, INRAE, CIRAD, Institut Agronomique Méditerranéen de Montpellier, CIHEAM, University of Montpellier, Montpellier, France, ² LEPSE, Univ Montpellier, INRAE, Institut Agro, Montpellier, France, ³ African Integrated Plant and Soil Research Group (AiPlaS), AgroBioSciences, Mohammed VI Polytechnic University, Ben Guerir, Morocco

Crop diversity management in agriculture is a fundamental principle of agroecology and a powerful way to promote resilient and sustainable production systems. Pulses are especially relevant for diversification issues. Yet, the specific diversity of legumes is poorly represented in most cropping systems. We used the trait-based approach to quantify the functional diversity of 30 pulses varieties, belonging to 10 species, grown under common field conditions. Our aim was to test relationships between traits, yield, and supporting agroecosystem properties. Our experimental results highlighted trade-offs between agroecosystem properties supported by different combinations of traits. Also, results demonstrated the relevance of leaf nitrogen content (LNC), leaf area ratio (LAR), and reproductive phenology to predict most of the trade-offs observed between agroecosystem properties. A comparison with a previous analysis based on literature data collected in diverse agronomic situations suggested that some traits are more plastic than others and therefore contribute differently to frame legumes diversity depending on the conditions of observation. Present results suggested that the implementation of such trait-based approach would rapidly benefit the selection of species/varieties for specific targeted agroecosystem services provisioning under specific (environmental or management) conditions.

**Keywords:** agroecosystem services, functional diversity, leaf nitrogen content, legumes, pulse crops, trait-based agroecology, traits-services trade-offs

## INTRODUCTION

Crop diversity management in agriculture is a fundamental principle of agroecology and a powerful way to promote resilient and sustainable production systems (Gaba et al., 2015). Agroecosystem management offers multiple alternatives to monoculture, such as diversification with pulses, through intercropping or rotation. Their well-known ability to fix atmospheric nitrogen reduces fossil energy consumption and makes them particularly suitable for low-input systems. They are also particularly important in human nutrition as a source of proteins to complement cereal-based diets (Singh et al., 2003; Tharanathan and Mahadevamma, 2003). Yet, the specific diversity of

**OPEN ACCESS**

**Edited by:**
Eric Von Wettberg,
University of Vermont, United States

**Reviewed by:**
Edward Marques,
University of Vermont, United States
Agnieszka Klimek-Kopyra,
University of Agriculture in Krakow,
Poland
Anoob Prakash,
University of Vermont, United States

*Correspondence:
Denis Vile
denis.vile@inrae.fr

**Specialty section:**
This article was submitted to Crop and Product Physiology, a section of the journal Frontiers in Plant Science

**Received:** 02 December 2020
**Accepted:** 08 March 2021
**Published:** 31 March 2021

**Citation:** Guiguitant J, Vile D and Marrou H (2021) *Trait Diversity of Pulse Species Predicts Agroecosystem Properties Trade-Offs*. *Front. Plant Sci.* 12:636915.
doi: 10.3389/fpls.2021.636915
legumes is poorly represented in most cropping systems. While more than 80 pulses species are known to contribute to human diet, the FAO database includes only 11 of them (Tiwari et al., 2011). This lack of attention creates a risk of continued erosion of knowledge and genetic resources that could further hamper the integration of these crops into more diversified cropping systems (Dansi et al., 2012).

To increase diversity within cropping systems, choosing appropriate species is an important prerequisite (Brooker et al., 2015). The trait-based approach is well suited for designing diversified cropping systems capable of providing specific sets of services. In ecology, the approach has been increasingly employed to link biodiversity to ecosystem functioning and services (Lavorel and Garnier, 2002; Violle et al., 2007; Garnier and Navas, 2012); a methodology that could be applied to cropping systems (Litrico and Violle, 2015; Martin and Isaac, 2015; Wood et al., 2015; Damour et al., 2018). However, trait-based approach requires extensive data collection to develop trait databases for crop species or varieties and document relationships between traits and services (Barot et al., 2017). As a neglected group of species, pulses are particularly difficult to document (Barot et al., 2017), for implementing a trait-based approach.

In a previous study, Guiguitant et al. (2020) set up a database suitable for functional trait approach on pulses. The database gathered disparate data from multiple published and/or public sources and contained data on plant functional traits for 43 pulse species, and agroecosystem properties recorded at canopy level for a subset of these species. The analysis of the functional trait space highlighted three ecological strategies among cultivated pulses that were similar to those already described in spontaneous species in the leaf-seed-height (LSH) scheme (Westoby, 1998). However, reduced data on agroecosystem properties and data disparity did not allow to identify clear link between traits and agroecosystem properties. In an applied perspective of supporting the integration of pulses in cropping systems, a consolidated analysis of trait–properties relationships is still needed to identify which traits would make legume species more productive, competitive, or adapted to low nitrogen or water conditions and favor production services. Trade-offs between plant traits and how they determine the different strategies also require more attention.

The aim of the present study was to identify pulses traits and/or trait combinations that can be used as predictors of yield and supporting agroecosystem services, further accounting for the relationships between traits and between agroecosystem properties. Original data were collected on 30 pulse varieties from 10 different species monitored in similar environmental conditions. Agroecosystem properties were measured at the canopy level, whereas functional traits, i.e., morphological, physiological, or phenological traits supposed to have direct or indirect effects on agroecosystem properties, were measured at the individual or organ level. We combined these original data collected within uniform environmental conditions with literature data collected in diverse agronomic situations (see Guiguitant et al., 2020) to assess the effect of trait measurement conditions on trait values and whether it could distort the relationships between traits and agroecosystem properties.

# MATERIALS AND METHODS

## Plant Material and Growing Conditions

A field experiment was conducted at UCOSEM's experimental farm in Lectoure (lat./long 43.911158/0.666182) on a deep clay-limestone soil. Three varieties for each one of the 10 pulses species (Table 1) were planted in two sequential batches/sub-experiments on contiguous fields of approximately 1 ha each. According to their temperature requirements: Cool season legumes were sown after the last negative temperatures, on March 14, 2018, and warm season legumes were sown on June 7, 2018, when daily temperature was over 15°C. Each batch/sub-experiment was laid out in split plot design with four replicates, varieties being subplots and species the main plot. Plots were 6 m long and counted four rows. Preceding crop was sunflower, and fields were fertilized with 0/22/8 NPK (350 kg ha⁻¹) fertilizer before sowing. Each species was manually inoculated before sowing with an adequate strain of rhizobia (provider: UMR LSTM, Montpellier, France). Sowing density was chosen according to agronomical recommendations for each species and ranged from 10 to 150 plants m⁻². Manual weeding and drip irrigation were performed on a regular basis from sowing. Soil moisture was recorded weekly in each block and for each species at three depths (0–20, 20–40, and 40–60 cm) using TDR probes installed vertically at sowing (mini TRASETM system with buriable probes, Soilmoisture Equipment Corp., CA) and confirmed that the plants did not experiment any water deficit in none of the two sub-experiments.

## Measurements of Individual Plant Traits

Occurrence of phenological stages (date of 50% flowering, date of physiological maturity, i.e., when more than 50% of the pods were yellow-brown) was recorded in calendar days and then converted into thermal units using a base temperature of 0°C for *Cicer arietinum* and *Vicia faba*, 2°C for *Lens culinaris*, 3.5°C for *Lathyrus sativus*, 4.1°C for *Trigonella foenum-graecum* (Covell et al., 1986; Lamichhane et al., 2020), 4.7°C for *Lotus tetragonolobus* (Moot et al., 2000), 6°C for *Phaseolus vulgaris* (López et al., 2003), and 8.5°C for *Vigna* sp. (Craufurd et al., 1996). This allowed to compute degree days from sowing to flowering (FLO) and degree days from sowing to maturity (MAT).

Biomass was harvested at flowering stage and maturity. At each sampling date, plants were collected on 0.5 m² quadrats, in the two central rows to avoid border effect. Plant number was recorded, and then leaves, stems, flowers, and pods were detached and then dried separately at 60°C until constant weight. Leaf and stem dry mass ratio over total aboveground plant mass (LMR and SMR, respectively) was calculated at flowering, and harvest index (HI) was calculated at maturity. Leaf area ratio (LAR; calculated as total plant leaf area (LA) over total aboveground mass), total LA (mm² plant⁻¹), specific LA (SLA, cm² g⁻¹), and leaflet length (LL, mm) as well as leaflet width (LW, mm) were measured at flowering (mean of the measurement of every leaf of a plant sub-sample), whereas thousand seed weight (TSW, g) was determined at harvest.
<table><caption>**TABLE 1** | List of the 10 species included in the experiment with a description of the three varieties associated to each species, the abbreviation (Abb.) used in the study, sources of the seeds, and the β value used for %Ndfa computation.</caption><thead><tr><th>Species</th><th>Varieties</th><th>Abb.</th><th>Source</th><th>β (‰)</th></tr></thead><tbody><tr><td>Vigna mungo</td><td>ISRA 58-77 (V1), ISRA 66-41 (V2), ISRA 67-30 (V3)</td><td>Vmu</td><td>ISRA, Sénégal</td><td>-1,75</td></tr><tr><td>Vigna unguiculata</td><td>YACINE (V1), SAM (V2), IT98K-1092-1-1 (V3)</td><td>Vun</td><td>ISRA, Sénégal</td><td>-1,61</td></tr><tr><td>Vigna aconitifolia</td><td>SH VM (V1), SH VA (V2), Tvu 45 35 (V2)</td><td>Vac</td><td>ISRA, Sénégal</td><td>-0,91</td></tr><tr><td>Lotus tetragonolobus</td><td>Landrace</td><td>Lte</td><td>Palmbeach medicinal herb, rareexoticseed, hobbyseed</td><td>-0,12</td></tr><tr><td>Phaseolus vulgaris</td><td>Red kidney (V1), zorro (V2), michelet (V3)</td><td>Pvu</td><td>Epi de Gascogne, France</td><td>-2,16</td></tr><tr><td>Cicer arietinum</td><td>Elmo (V1), billy bean (V2), orion (V3)</td><td>Car</td><td>Epi de Gascogne, France; USDA Pullman, WS</td><td>-1,75</td></tr><tr><td>Lens culinaris</td><td>Belezana (V1), Richela (V2), Anicia (V3)</td><td>Lcu</td><td>Epi de Gascogne, USDA Pullman, WS</td><td>-0,56</td></tr><tr><td>Trigonella foenum-graecum</td><td>FENUSOL, landrace</td><td>Tfo</td><td>Epi de Gascogne, France; ICARDA, Morroco</td><td>-0,42</td></tr><tr><td>Lathyrus sativus</td><td>Wassie (V1), Almenaza (V2), Chicharo (V3)</td><td>Lsa</td><td>Epi de Gascogne, France</td><td>-0,38</td></tr><tr><td>Vicia faba</td><td>Aguadulce (V1), alfia (V2), Iobaba (V3)</td><td>Vfa</td><td>ICARDA, Morroco</td><td>-0,5</td></tr></tbody></table>

Total nitrogen content was determined (Vario PyroCube auto-analyzer and Precision mass spectrometer) on leaves, stems, and grain samples collected at maturity and at flowering. Total plant nitrogen content was computed both at flowering and maturity. Leaf nitrogen content (LNC, mg g⁻¹), and seed nitrogen content (SNC, mg g⁻¹) were determined at flowering and at maturity, respectively.

Phyllochron was calculated from observations of three tagged plants on each plot for which emerged leaves were counted weekly (cool season legumes) or twice a week (warm season legumes).

The rate of soil coverage was estimated by measuring soil cover weekly, using nadiral photos taken at the same location within each plot, with a large angle camera at chest height, and then analyzed with ImageJ software (Schneider et al., 2012). Soil cover was modeled as a logistic function of thermal time. The maximal rate of soil coverage was defined as the slope at the inflexion point of the logistic curve.

## Measurement of Agroecosystem Properties

Agroecosystem properties were measured on the same experimental plots (Table 2). Grain yield (GY, t ha⁻¹), biomass yield (BY, t ha⁻¹), and biomass measured at flowering (t ha⁻¹) were determined from plant sampling on 0.5 m² quadrats. Water use efficiency (WUE) was computed as grain mass produced per mm of water input (rainfall and irrigation). Maximum soil cover (%), LA index (LAI), and the duration (in thermal units) of the phase when soil coverage is above 80% of the maximum (soil cover duration) were recorded. Finally, the capacity to fix nitrogen from the atmosphere was estimated through the percent of nitrogen derived from biological fixation at flowering and at maturity, measured with the natural isotopic dilution method, using barley and maize grown in the experiment outer borders as control plants for cool season and warm season legumes, respectively. From the ¹⁵N abundance (d¹⁵N) of the samples, the proportion of plant N derived from atmosphere (%Ndfa) was calculated from Eq. 1 (Shearer and Kohl, 1986):

$$ \%Ndfa = 100 \frac{\delta15N_{ref} - \delta15N_{legume}}{\delta15N_{ref}^{-\beta}} \quad (1) $$

<table><caption>**TABLE 2** | List of agroecosystem properties and traits measured in the experiment with their respective abbreviation and unit.</caption><thead><tr><th>Variable</th><th>Abbreviation</th><th>Unit</th></tr></thead><tbody><tr><td colspan="3"><strong>Agroecosystem property</strong></td></tr><tr><td>Grain yield</td><td>GY</td><td>t ha<sup>-1</sup></td></tr><tr><td>Biomass yield</td><td>BY</td><td>t ha<sup>-1</sup></td></tr><tr><td>Leaf area index</td><td>LAI</td><td></td></tr><tr><td>% Nitrogen derived from atmosphere at maturity</td><td>Ndfa_MAT</td><td>%</td></tr><tr><td>% Nitrogen derived from atmosphere at flowering</td><td>Ndfa_FLO</td><td>%</td></tr><tr><td>Water use efficiency</td><td>WUE</td><td>kg ha<sup>-1</sup> mm<sup>-1</sup></td></tr><tr><td>Biomass at flowering</td><td>B_FLO</td><td>t ha<sup>-1</sup></td></tr><tr><td>Maximal soil cover</td><td>Maximum soil cover</td><td>%</td></tr><tr><td>Soil cover duration</td><td>Soil cover duration</td><td>Thermal unit</td></tr><tr><td colspan="3"><strong>Traits</strong></td></tr><tr><td>Leaf area</td><td>LA</td><td>mm<sup>2</sup></td></tr><tr><td>Specific leaf area</td><td>SLA</td><td>cm<sup>2</sup> g<sup>-1</sup></td></tr><tr><td>Leaflet length</td><td>LL</td><td>mm</td></tr><tr><td>Leaflet width</td><td>LW</td><td>mm</td></tr><tr><td>Leaf nitrogen content</td><td>LNC</td><td>mg g<sup>-1</sup></td></tr><tr><td>Plant N content at flowering</td><td>Nplant_FLO</td><td>g plant<sup>-1</sup></td></tr><tr><td>Plant N content at maturity</td><td>Nplant_MAT</td><td>g plant<sup>-1</sup></td></tr><tr><td>Phyllochron</td><td>Phyllochron</td><td>Thermal unit</td></tr><tr><td>Soil cover rate</td><td>Cover rate</td><td></td></tr><tr><td>Degree days to maturity</td><td>MAT</td><td>Thermal unit</td></tr><tr><td>Degree days to flowering</td><td>FLO</td><td>Thermal unit</td></tr><tr><td>Thousand seed weighed</td><td>TSW</td><td>G</td></tr><tr><td>Leaf area ratio</td><td>LAR</td><td>–</td></tr><tr><td>Leaf mass ratio</td><td>LMR</td><td>–</td></tr><tr><td>Stem mass ratio</td><td>SMR</td><td>–</td></tr><tr><td>Harvest index</td><td>HI</td><td>–</td></tr><tr><td>Seed nitrogen content</td><td>SNC</td><td>mg g<sup>-1</sup></td></tr></tbody></table>

Where β is the ¹⁵N natural abundance of the N derived from biological nitrogen fixation (BNF) in the aerial tissue of the legume, δ15N<sub>ref</sub> is the ¹⁵N abundance in the control plant (maize or barley) tissues, and δ15N<sub>legume</sub> is the ¹⁵N abundance in legume tissues. The β values used are detailed in **Table 1** (Unkovich et al., 2008).
## Data Analysis

Data were first analyzed through the prism of the framework presented in Guiguitant et al. (2020) in order to compare the variability observed with the global functional diversity of pulses species. Data of the present field measurements were projected as supplementary individuals in a probabilistic principal component analysis (PPCA) (Tipping and Bishop, 1999) computed with nine common functional traits (DM, DF, HI, LL, LW, LA, SLA, LNC, and TSW) values taken from Guiguitant et al. (2020) database. We used Euclidian distance to quantify the distance between the scores of barycenters of each species as calculated from field data and previously published in Guiguitant et al. (2020). The consistency of species ranking in trait values across literature and experimental data was tested with Spearman's rank correlation analysis.

Second, a principal component analysis (PCA) was performed on experimental data including all traits measured. Agroecosystem properties were added as supplementary variables, in order to investigate how they fit into the phenotypic space defined by the traits. The appropriate number of PCs was determined with Kaiser' criterion which selects components that correspond to eigenvalues larger than 1. A between and within-class correspondence analysis was carried out to evaluate the percentage of variance captured by the intra and inter-specific diversity (Doledec, 1987).

Finally, trade-offs among agroecosystem properties were analyzed using a hierarchical clustering analysis (HCA). A classification and regression tree (CART) analysis was done on raw data (n = 120) to identify how functional traits could allow to predict the clusters of varieties obtained from the HCA. Regression trees are prediction models obtained through machine-learning algorithms that recursively partition the data space in order to fit the simplest prediction model within each partition. The resulting partitioning can be represented graphically as a decision tree. The random forest aggregation method, which consists in adding a random component to the choice of variables in the model, was also used. This method does not allow to generate a visual decision tree, but it provides indices proportional to the importance of each variable in the aggregated model and thus its participation to the discrimination between groups of individuals (here, varieties). Mean decrease in Gini impurity (MDG) can be defined as the total decrease in node impurity (weighted by the proportion of samples reaching a given node) averaged across all of the trees that make up the forest. The most important variables to the model will have the largest MDG values; conversely, the least important variable will have the smallest MDG values.

All statistical analyses were performed in the computing environment R 3.5.2 (R Core Team, 2018) using *pcaMethods* package (Stacklies et al., 2007), *rpart* package (Therneau et al., 2015), *ade4* package (Dray and Dufour, 2007), and *randomForest* package (Breiman, 2001).

## RESULTS

### Consistency of Pulses Functional Trait Space

A PPCA was performed on the 43-species-database from Guiguitant et al. (2020) using only traits that were also measured in the present field experiment. The first three PC axes of this PPCA explained 70% of total variance observed in the literature (Figure 1). PC1 (34%) was associated with morphological traits, especially leaf traits (LL, LW, LA, and to a lesser extent TSW). PC2 (24%) was positively associated with late flowering (DF) and maturity (DM) and to a lesser extent HI. PC3 (12%) was strongly associated to SLA and LNC. Experimental data variability on PC axes was in the same range as the literature variability. Experimental variance was close to 100 and 60% of literature variance on PC1 and PC2, respectively. On the contrary, the variance of supplementary individuals on the third axis (driven by SLA and LNC) was almost three times higher than the variance of active data from the literature (Appendix 1). For all other traits, field measurements range was always narrower than the one observed in Guiguitant et al. (2020), for a larger set of species (Appendix 1).

Distance between the barycenters of each species as tested in the present experiment and from the literature database was short in the first plane of the PPCA (Figure 1). Species ranking in literature versus present data remained conserved on PC1 ($ρ = 0.85, p = 0.003$) and to a lesser extent on PC2 ($ρ = 0.53, p = 0.10$). Species rankings were rather inconsistent on PC3 ($ρ = 0.24, p = 0.12$). Cool season species had barycenters closely positioned along the three dimensions of the PPCA (Euclidian distance lower than 2; Table 3). Among them, *C. arietinum* and *T. foenum-graecum* had the highest Euclidian distance mostly because of the divergence with literature on the third axis, while *V. faba* diverged mostly on the second axis. For warm season species, distances between barycenters from the literature and from the present experiment were larger, on PC3 but also on PC2 for *Vigna* sp.

### Characterization of Functional Trait Diversity

The first four axes of the PCA performed on 30 individuals and 16 traits from the present experiment (Table 2) captured 80% of the total variance (Figure 2). Interspecies variance explained 81% of the total variance, while intra-species explained 19%.

Leaf attributes were strongly associated with PC1 (40%). On the first principal component (PC1), LA, LL, and plant N content at flowering were negatively correlated with SLA, LNC, SNC, and LAR. This axis discriminated cool season species which were characterized by small leaflets and small LA with high leaf and seed nitrogen content and a large LAR, from warm season species characterized by large leaflet, high LA, and high plant nitrogen content. PC2 (20%) was positively correlated to LMR and negatively correlated to SMR and cover rate. The cover rate was significantly negatively correlated to MAT ($ρ = -0.72, p < 0.001$). SMR was also significantly negatively correlated to the phyllochron ($ρ = -0.60, p < 0.001$). *L. tetragonolobus* and
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Principal component analysis (PCA) performed on nine functional traits collected on 43 pulses species. (A,C) Individual species projection on the first three axes of the PCA; light gray points represent the 43 species from literature (Guiguitant et al., 2020) and colored points are from the current experimental study projected as supplementary individuals. Transparent points are the data of the four plots for the three varieties of the 10 species, the larger opaque points represent the barycenter of each species, and the horizontal and vertical bars represent the standard deviation along axes. Visualization of the variables and correlation circles on (B) PC1–PC2 and (D) PC2–PC3 planes. See **Table 2** for abbreviations of traits.</figcaption></figure>

P. vulgaris were apart from all other species on PC2, and both expressed a high LMR at flowering. PC3 (11%) and PC4 (9%) were associated to FLO and TSW, respectively. PC3 separated V. faba from L. culinaris and L. tetragonolobus, whereas PC4 discriminated C. arietinum and P. vulgaris.

All agroecosystem properties, except LAI, biomass at flowering and Ndfa at maturity were not well represented in the factor planes formed by the fourth PCs constructed with traits (Figure 2). LAI and biomass at flowering were relatively well represented in the first plane and were significantly positively correlated to LL ($\rho$ = 0.46, $p$ < 0.001; $\rho$ = 0.63, $p$ = 0.01, respectively) and LA ($\rho$ = 0.67, $p$ < 0.001; $\rho$ = 0.83, $p$ < 0.001, respectively) as well as cover rate ($\rho$ = 0.59, $p$ < 0.001; $\rho$ = 0.60, $p$ < 0.001, respectively). Although Ndfa at flowering was poorly represented on the first PC, this property was significantly

positively correlated to plant N content at flowering ($\rho$ = 0.41, $p$ = 0.02). Ndfa at maturity was not correlated to Ndfa at flowering ($\rho$ = 0.04, $p$ = 0.83) but was significantly positively correlated to LNC ($r$ = 0.49, $p$ = 0.005) and SLA ($\rho$ = 0.57, $p$ = 0.001) as well as with SMR ($\rho$ = 0.43, $p$ = 0.02). The soil cover duration was inversely correlated to LAI ($\rho$ = -0.21, $p$ = 0.0428) and biomass at flowering ($\rho$ = -0.22, $p$ = 0.033) and significantly positively correlated with FLO ($\rho$ = 0.47, $p$ = 0.009) and to a lesser extent with MAT ($\rho$ = 0.33, $p$ = 0.08). Maximum soil cover was not correlated to soil cover duration ($\rho$ = -0.03, $p$ = 0.86) and was slightly associated with SMR ($\rho$ = 0.35, $p$ = 0.05). Finally, GY and WUE were positively correlated to PC3 which was mostly determined by FLO, whereas BY, and to a lesser extent GY and WUE, were positively associated to PC4 alongside TSW.
<table><caption>**TABLE 3** | Euclidian distance between literature and current study coordinates of the barycenter of each species on the firsts three PCA axes.</caption><thead><tr><th>Species</th><th>PC1</th><th>PC2</th><th>PC3</th><th>Three PCs</th></tr></thead><tbody><tr><td>Cicer arietinum</td><td>1.1</td><td>0.39</td><td>1.25</td><td>1.68</td></tr><tr><td>Lathyrus sativus</td><td>1.0</td><td>0.18</td><td>0.11</td><td>1.03</td></tr><tr><td>Lens culinaris</td><td>0.4</td><td>0.25</td><td>1.29</td><td>1.39</td></tr><tr><td>Lotus tetragonolobus</td><td>0.9</td><td>1.41</td><td>1.74</td><td>2.40</td></tr><tr><td>Phaseolus vulgaris</td><td>0.4</td><td>0.92</td><td>3.74</td><td>3.87</td></tr><tr><td>Trigonella foenum-graecum</td><td>0.8</td><td>0.44</td><td>1.73</td><td>1.97</td></tr><tr><td>Vicia faba</td><td>1.1</td><td>1.20</td><td>0.22</td><td>1.64</td></tr><tr><td>Vigna aconitifolia</td><td>1.0</td><td>2.61</td><td>2.96</td><td>4.08</td></tr><tr><td>Vigna mungo</td><td>0.6</td><td>1.21</td><td>1.76</td><td>2.23</td></tr><tr><td>Vigna unguiculata</td><td>0.7</td><td>1.23</td><td>2.39</td><td>2.77</td></tr></tbody></table>

## Trade-Offs Between Agroecosystem Properties

The classification analysis revealed seven clusters of species and varieties based on their agroecosystem properties (Figure 3). All varieties of a species generally grouped into the same cluster. *V. aconitifolia* varieties were all in cluster 3 together with one variety of *Vigna unguiculata*. This cluster is apart from a group of close clusters composed essentially by cool season species; *L. sativus* varieties were grouped with *L. tetragonolobus* (cluster 4), *L. culinaris* was alone in cluster 5, cluster 6 was dominantly occupied by *T. foenum-graecum*, and cluster 7 grouped both *V. faba* and *Vigna mungo*. These groups were distant from cluster 1 essentially represented by *C. arietinum* and cluster 2 that contained varieties of *P. vulgaris*, and of *V. unguiculata*.

Each cluster was characterized by a combination of values of agroecosystem properties. Cluster 1 grouped species with high production (GY and BY) as well as a high maximum soil cover and soil cover duration, but a low Ndfa and a low biomass and LAI at flowering. Cluster 2 was characterized by a low biomass at flowering as well as a low LAI at flowering and maximum soil cover but a high Ndfa at flowering and a long soil cover duration. Cluster 3 varieties expressed the lowest yield as well as the lowest WUE but high Ndfa at maturity, a high LAI at flowering, and a long cover duration. Cluster 4 varieties expressed the highest WUE concomitantly with a high GY and a low Ndfa at both stages. Varieties of cluster 5 expressed high Ndfa both at flowering and maturity and high maximum soil cover. Yet, they showed a low maximum soil cover and soil cover duration. Cluster 6 was characterized by a low BY and GY, low Ndfa both at flowering and maturity, and poor soil covering capacity, yet LAI and soil cover at flowering were high. Finally, cluster 7 was composed of varieties producing the highest biomass at flowering, along with the highest maximum soil cover, LAI, and Ndfa at flowering.

## Relevant Traits for Prediction of Agroecosystem Properties and Structuration of Functional Diversity

Classification and regression tree analysis showed very diverse trait combinations to predict the cluster affiliation (Figure 4A). LNC was identified as the most important trait to discriminate the clusters as it had the highest MDG value, which means that it

allowed to gain purity in the process of classification for most of the nodes in the CART it was associated with (Figure 4B). LNC was assigned to the first division node and allowed to sort cluster 1 varieties apart from the remaining ones. Cluster 1 varieties have an LNC between 2.7 and 3 mg N g⁻¹. Varieties of cluster 4 were separated. Half of them fell on the same side of the tree as cluster 1 but exhibited lower LNC value (below 2.7 mg N g⁻¹), while the other half had an average LNC above 3 mg N g⁻¹. Varieties in this second cluster 4 half also had an LMR above 43% concomitantly with low plant N content at flowering and at maturity and a longer cycle. LMR, DDR, and plant N content at maturity had strong MDG unlike the plant N content at flowering. Varieties with an LMR under 0.43 and having high MAT were assigned either to cluster 3 (characterized by a high plant N content but only at flowering), or cluster 2 or 5 (both characterized by high plant N content at flowering and maturity). Varieties in cluster 2 flowered later than those in cluster 5. Short cycle varieties were dominantly associated to cluster 6 but were also found in cluster 5 if they also expressed a high LA. Finally, cluster 7 was characterized by an LMR under 43%.

Relying on MDG, plant N content at maturity was the second most important trait to discriminate the clusters. Phyllochron, SNC, SLA, and LL were not selected in the tree and had very low mean decrease in Gini value. LAR was not used in the trees, but had a higher mean decrease in Gini value than LA and LMR. Yet, LA and LMR captured more variability in the PCA than LAR. Finally, cover rate did not appear on the tree despite its high mean decrease in Gini value, probably due to its negative correlation with DDR.

## DISCUSSION

### Traits Related to Resource Acquisition Are More Plastic Across Growing Conditions and Varieties Than Architectural Traits

Through functional traits, researchers usually seek to make generalizable predictions across organizational and spatial scales (Adler et al., 2013) and therefore usually use data collected in diverse environment and management contexts. On the contrary, field measurements sharing similar environment as in the present study deliver a local picture of functional diversity which variations can be solely attributed to genotypic determinants. Several studies have demonstrated that despite the significant plasticity of most plant traits across environments, the classification of species for a given trait remains fairly consistent in space and time (Garnier et al., 2001; Kazakou et al., 2014; Tribouillois et al., 2015). Yet, other studies have reported a significant effect of plasticity on species ranking across environments, particularly with respect to traits associated with light acquisition and nitrogen nutrition (Lipowsky et al., 2015; Siebenkäs et al., 2015; Keenan and Niinemets, 2017). These inconsistencies may impair the use of simple aggregation functions (means, medians, etc.) for generalized trait-based approaches at the species level.
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Principal component analysis (PCA) performed on traits measured on three varieties of 10 species. Visualization of the variables and correlation circles on (A,C) PC1–PC2 and (B,D) PC3–PC4 planes. Gray scale represents the contribution of each variable to axes construction. Orange variables are agroecosystem properties, used as supplementary variables. See **Table 2** for abbreviations of traits and properties. Soil cover stands for maximal soil cover.</figcaption></figure>

In the present study, we compared common field measurements to the meta-analysis by Guiguitant et al. (2020) to analyze how environmental and management conditions could alter trait values and distort their relationships. When using the experimental dataset, where all data were acquired in a common environment, species rankings on the PC axis summarizing the variation of architectural traits were fairly consistent with literature ranking. However, rankings on the axis summarizing SLA and LNC variation diverged between the two datasets. This contrast mainly concerned warm-season species, particularly *Vigna* sp. for which divergence between the two datasets was also noted in terms of phenology. Number of days required to reach maturity was generally lower in literature than in experimental data. This comparison between local measurements and literature averages suggests that some

traits are more plastic than others and therefore contribute differently to frame legumes diversity depending on the conditions of observation.

## Intraspecific Variation Contributed Only Marginally to Global Variability

Results from the PCA showed that the total trait variance was mainly explained by the variance between species. Variation in intraspecific traits is generally considered of marginal importance compared to differences in interspecific traits (McGill et al., 2006; Hulshof and Swenson, 2010). However, several recent studies have shown that intraspecific variation in traits accounts for a substantial proportion of overall plant community variability (Jung et al., 2010; Albert et al., 2011; Roscher et al., 2015;
<figure><img src="image_5.png" /></figure>

Siefert et al., 2015). The small number of varieties considered in this study may have not been able to highlight such diversity.

## Leaf Morphological Traits Confirmed as Variability Drivers Among Pulses

Guiguitant et al. (2020) found that seed and leaf traits exhibited the greatest variability among 43 pulse species. In the present study, leaf traits (LA and LL) remained an important variability driver that was correlated to the first PCA axis, whereas seed size, associated to the first axis in the literature analysis, poorly contributed to the total variance across varieties (associated only with the fourth PCA axis) and was decoupled from leaf size, which is consistent with LHS scheme (Westoby, 1998).

The experimental data confirmed that LL is an important determinant of legume diversity. This trait was strongly related to LNC and to a lesser extent SLA. This finding can be related to the well known leaf economic spectrum (LES; Wright et al., 2004). The LES reflects a trade-off between a leaf's lifespan and its maximum net photosynthetic rate. Yet, relations between LNC and photosynthesis in grain legumes have been questioned (Adams et al., 2016). Thus, the fact that this axis reflects a gradient in leaf photosynthetic capacity is uncertain. However, a high LNC could have other biological implications such as remobilization to the reproductive organs (Davies et al., 2000; Schiltz et al., 2005). Indeed, SNC was correlated with LNC in our experiment. Small leaflet associated with a high LNC could support the growth of N-rich legume seeds through remobilization.
<figure><img src="image_6.png" /></figure>

**FIGURE 4 | (A)** Classification and regression tree (CART) regression for prediction of membership to one of the seven clusters based a hierarchical clustering analysis of the 30 pulse varieties based on agroecosystem properties (see **Figure 3**). Root nodes represent single input variables (traits) and related split point used to make the prediction; the best discriminant is chosen at each step of the classification procedure. Leaf nodes contain mean, number, and percentage of observations of the predicted variable (cluster membership). (B) Gini decrease value for each trait obtained through a random forest procedure. A high Gini decrease value means the variable considered allowed to gain purity in the process of classification for most of the nodes in the CART it was associated with during the random forest procedure.
Yield Is Inversely Related With Ndfa and Soil Cover

Nitrogen fixation by legumes is generally expected to contribute to improving soil fertility through the mineralization of crop residues such as senescent leaves or stems, or roots. In the current study, Ndfa at flowering and Ndfa at harvest in grain tissue varied almost independently. Clustering reveled three types of "high fixing species": cluster 7, cluster 3, and cluster 5. On the other hand, cluster 1 or 6 expressed a low N fixation capacity. "High yielding" species were clustered in cluster 1 and, to a lesser extent, cluster 4 both concerned by low fixation. Cluster 5 fixed a lot of N during flowering and produced a decent yield. However, yield in cluster 5 was lower than that of clusters 1 or 4. Conversely, the very low yielding species (cluster 3) expressed a significant increase in relative Ndfa between flowering and maturity. Other multiple species comparisons confirmed the negative relation between BNF and seed yield (Liu et al., 2019). Some studies reported BNF over a certain level results in a reduction of yield (Mesfin et al., 2020).

The multi-faceted benefits of cover crops are well known (Blanco-Canqui et al., 2015), but land cover is often not considered for cash crops and the interaction between cover capacity and yield is rarely discussed. Maximum cover and cover duration appeared to be related to biomass at flowering and LAI as well as to yield properties. Cluster 7 expressed an important cover at flowering as a result of an important biomass and a large LA. However, cluster 1 and 3 had a low biomass at flowering and a relatively high maximum cover that lasts long, probably because of a delayed senescence. Except for cluster 1, most of those species had low yield. Delayed senescence, by prolonging photosynthesis, is usually expected to contribute to increasing the yield of high-carbon crops such as cereals (Gregersen et al., 2013). However, studies suggest that the stay-green trait may have limited value for high N species (Ismail et al., 2000; Kumudini, 2002) as its value lies in an improvement of remobilization capacity. This corroborates the absence of relation we found between yield and duration of ground cover.

Trait-Based Predictions of Agroecosystem Properties

Understanding the links between functional traits and agroecosystem functioning is fundamental for optimal management of agroecosystem services (Zhang et al., 2007; Litrico and Violle, 2015; Wood et al., 2015). As noted earlier, leaf characteristics and especially LNC were crucial in describing pulse diversity. Other traits such as LAR, total nitrogen, seed size, and cycle length are also important, which is consistent with the study of Guiguitant et al. (2020). The decision tree did not highlight all these traits but discriminated specific functional characteristics for each cluster. Clusters 1 and 4 (high yielding species) were characterized by low LNC values which is consistent with their low N fixation. Cluster 7 (high biomass at flowering) was characterized by a large LA but a low biomass investment in leaves. Clusters 6 and 5 (medium maximum soil cover and low value for every other properties) were characterized by a short cycle, with cluster 5 (high nitrogen

fixation) having a higher LA. Clusters 2 (medium value for every properties) and 3 (high N fixation at maturity and long cover duration) were characterized by a long cycle and high nitrogen content at maturity or flowering, respectively.

To conclude on the highlighted linkages, a high LNC value does not translate into higher yields but is generally associated with a high N fixation. Investment in LA is important for N₂ fixation at flowering but is associated with rapid senescence and low yield. Biomass investment in leaves and a short cycle coincide with low values for almost every property except maximum soil cover. Finally, long cycle and high N content at maturity are related to delayed senescence and maintenance of N₂ fixation along pod filling.

CONCLUSION

Integrating pulses in resilient and sustainable production systems requires a better understanding of their functional diversity and how this functional diversity scales up to provide specific sets of agroecosystem services. Previous studies demonstrated the usefulness of trait value aggregation at the species-level to implement trait-based approaches for a better knowledge of the plant-agroecosystem functioning relationships and increasing our predictive capacity (Lavorel and Garnier, 2002; Violle et al., 2007; Garnier and Navas, 2012; Kazakou et al., 2014; Martin and Isaac, 2015; Wood et al., 2015; Guiguitant et al., 2020). Our experimental results highlighted trade-offs and synergies between agroecosystem properties supported by different pulses species. We were able to demonstrate the relevance of LNC, LAR, and phenology to predict most of the trade-offs observed for agroecosystem properties in those conditions.

However, the present study advocates further for a systematic characterization of the functional diversity (and its plasticity), including intra-specific variability. We argue that the trait-based approach may be easier to implement in crop species where varieties origin and genetics are often better described and accessible, compared with the huge genetic variation encountered in populations of spontaneous species. Application of such trait-based approach would rapidly benefit the selection of species/varieties for specific targeted agroecosystem services provisioning under specific (environmental or management) conditions.

DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS

All authors conceived the ideas and approach together. JG conducted the experiments and data analyses and led the writing
of the manuscript. DV and HM contributed substantially to the manuscript through extensive revision. HM designed the experiment. DV participated to data analyses.

## FUNDING

This work was funded by the European Programme ARIMNet/ERA-NET, FP7-KBBE, session 2017–2020, SEMIARID, grant agreement no. 618127.

## REFERENCES

Adams, M. A., Turnbull, T. L., Sprent, J. I., and Buchmann, N. (2016). Legumes are different: leaf nitrogen, photosynthesis, and water use efficiency. *Proc. Natl. Acad. Sci. U.S.A.* 113, 4098–4103. doi: 10.1073/pnas.1523936113

Adler, P. B., Fajardo, A., Kleinhesselink, A. R., and Kraft, N. J. B. (2013). Trait-based tests of coexistence mechanisms. *Ecol. Lett.* 16, 1294–1306. doi: 10.1111/ele.12157

Albert, C. H., Grassein, F., Schurr, F. M., Vieilledent, G., and Violle, C. (2011). When and how should intraspecific variability be considered in trait-based plant ecology? *J. Perspect. Plant Ecol. Evol. Syst. Sources* 13, 217–225. doi: 10.1016/j.ppees.2011.04.003

Barot, S., Allard, V., Cantarel, A., Enjalbert, J., Gauffreteau, A., Goldringer, I., et al. (2017). Designing mixtures of varieties for multifunctional agriculture with the help of ecology. A review. *Agron. Sustain. Dev.* 37:13. doi: 10.1007/s13593-017-0418-x

Blanco-Canqui, H., Shaver, T. M., Lindquist, J. L., Shapiro, C. A., Elmore, R. W., Francis, C. A., et al. (2015). Cover crops and ecosystem services: insights from studies in temperate soils. *Agron. J.* 107, 2449–2474. doi: 10.2134/agronj15.0086

Breiman, L. (2001). Random forests. *Mach. Learn.* 45, 5–32. doi: 10.1023/A:1010933404324

Brooker, R. W., Jones, H. G., Paterson, E., Watson, C., Brooker, R. W., Bennett, A. E., et al. (2015). Improving intercropping?: a synthesis of research in agronomy, plant physiology and ecology. *New Phytol.* 206, 107–117. doi: 10.1111/nph.13132

Covell, S., Ellis, R. H., Roberts, E. H., and Summerfield, R. J. (1986). The influence of temperature on seed germination rate in grain legumes. I. A comparison of chickpea, lentil, soyabean and cowpea at constant temperatures. *J. Exp. Bot.* 37, 705–715. doi: 10.1093/jxb/37.5.705

Craufurd, P. Q., Ellis, R. H., Summerfield, R. R. J., and Menin, L. (1996). Development in cowpea (*Vigna unguiculata*). I. The influence of temperature on seed germination and seedling emergence. *Exp. Agric.* 32, 1–12. doi: 10.1017/S001447970025801

Damour, G., Navas, M. L., and Garnier, E. (2018). A revised trait-based framework for agroecosystems including decision rules. *J. Appl. Ecol.* 55, 12–24. doi: 10.1111/1365-2664.12986

Dansi, A., Vodouhè, R., Azokpota, P., Yedomonhan, H., Assogba, P., Adjatjan, A., et al. (2012). Diversity of the neglected and underutilized crop species of importance in benin. *Sci. World J.* 2012:932947. doi: 10.1100/2012/932947

Davies, S. L., Turner, N. C., Palta, J. A., Siddique, K. H. M., and Plummer, J. A. (2000). Remobilisation of carbon and nitrogen supports seed filling in chickpea subjected to water deficit. *Aust. J. Agric. Res.* 51, 855–866. doi: 10.1071/AR00018

Doledec, S. (1987). *Rythmes Saisonniers et Composantes Stationnelles en Milieu Aquatique. 1 - Description d'un Plan D'observations Complet par Projection de Variables*, Vol. 83, 403–426. Available online at: http://www.cbn-alpin-biblio.fr/Record.htm?idlist=1&record=19178667124919968499 (accessed February 17, 2021).

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

Gaba, S., Lescourret, F., Boudsocq, S., Enjalbert, J., Hinsinger, P., Journet, E. P., et al. (2015). Multiple cropping systems as drivers for providing multiple ecosystem services: from concepts to design. *Agron. Sustain. Dev.* 35, 607–623. doi: 10.1007/s13593-014-0272-z

## ACKNOWLEDGMENTS

The authors thank Laurent Palau and UCOSEM team in Lectreou for having the field experiment on their farm and ensuring proper crop husbandry, and for providing genetic resources. The authors are also grateful to the LSTM lab (B. Brunel) and the USDA ARS winter legume breeding program at Pullman (WS) (R. McGuee) for providing, respectively, rhizobial inoculum as well as plant genetic resources. The authors also thank ISRA Senegal for sharing Vigna accessions.

Garnier, E., Laurent, G., Bellmann, A., Debain, S., Berthelier, P., Ducout, B., et al. (2001). Consistency of species ranking based on functional leaf traits. *New Phytol.* 152, 69–83. doi: 10.1046/j.0028-646x.2001.00239.x

Garnier, E., and Navas, M. L. (2012). A trait-based approach to comparative functional plant ecology: concepts, methods and applications for agroecology. *Agron. Sustain. Dev.* 32, 365–399. doi: 10.1007/s13593-011-0036-y

Gregersen, P. L., Culetic, A., Boschian, L., and Krupinska, K. (2013). Plant senescence and crop productivity. *Plant Mol. Biol.* 82, 603–622. doi: 10.1007/s11103-013-0013-8

Guiguitant, J., Vile, D., Ghanem, M. E., Wery, J., and Marrou, H. (2020). Evaluation of pulse crops' functional diversity supporting food production. *Sci. Rep.* 10:3416. doi: 10.1038/s41598-020-60166-4

Hulshof, C. M., and Swenson, N. G. (2010). Variation in leaf functional trait values within and across individuals and species: an example from a Costa Rican dry forest. *Funct. Ecol.* 24, 217–223. doi: 10.1111/j.1365-2435.2009.01614.x

Ismail, A. M., Hall, A. E., and Ehlers, J. D. (2000). Delayed-leaf-senescence and heat-tolerance traits mainly are independently expressed in cowpea. *Crop Sci.* 40, 1049–1055. doi: 10.2135/cropsci2000.4041049x

Jung, V., Violle, C., Mondy, C., Hoffmann, L., and Muller, S. (2010). Intraspecific variability and trait-based community assembly. *J. Ecol.* 98, 1134–1140. doi: 10.1111/j.1365-2745.2010.01687.x

Kazakou, E., Violle, C., Roumet, C., Navas, M.-L., Vile, D., Kattge, J., et al. (2014). Are trait-based species rankings consistent across data sets and spatial scales? *J. Veg. Sci.* 25, 235–247. doi: 10.1111/jvs.12066

Keenan, T. F., and Niinemets, Ü (2017). Global leaf trait estimates biased due to plasticity in the shade. *Nat. Plants* 3:16201. doi: 10.1038/nplants.2016.201

Kumudini, S. (2002). Trials and tribulations: a review of the role of assimilate supply in soybean genetic yield improvement. *Field Crops Res.* 75, 211–222. doi: 10.1016/S0378-4290(02)00027-8

Lamichhane, J. R., Constantin, J., Schoving, C., Maury, P., Debaeke, P., Aubertot, J.-N., et al. (2020). Analysis of soybean germination, emergence, and prediction of a possible northward establishment of the crop under climate change. *Eur. J. Agron.* 113:125972. doi: 10.1016/j.eja.2019.125972

Lavorel, S., and Garnier, E. (2002). Predicting changes in community composition and ecosystem functioning from plant traits: revisiting the Holy Grail. *Funct. Ecol.* 16, 545–556. doi: 10.1046/j.1365-2435.2002.00664.x

Lipowsky, A., Roscher, C., Schumacher, J., Michalski, S. G., Gubsch, M., Buchmann, N., et al. (2015). Plasticity of functional traits of forb species in response to biodiversity. *Perspect. Plant Ecol. Evol. Syst.* 17, 66–77. doi: 10.1016/j.ppees.2014.11.003

Litrico, I., and Violle, C. (2015). Diversity in plant breeding: a new conceptual framework. *Trends Plant Sci.* 20, 604–613. doi: 10.1016/j.tplants.2015.07.007

Liu, L., Knight, J. D., Lemke, R. L., and Farrell, R. E. (2019). A side-by-side comparison of biological nitrogen fixation and yield of four legume crops. *Plant Soil* 442, 169–182. doi: 10.1007/s11104-019-04167-x

López, J. C., Pñrez-Parra, J., Baille, A., and Bonachela, S. (2003). Effects of heating strategies on earliness and yield of snap beans (*Phaseolus vulgaris* l.) grown under *Acta Hortic.* 614, 439–444. doi: 10.17660/ActaHortic.2003.614.66

Martin, A. R., and Isaac, M. E. (2015). Plant functional traits in agroecosystems: a blueprint for research. *J. Appl. Ecol.* 52, 1425–1435. doi: 10.1111/1365-2664.12526
McGill, B. J., Enquist, B. J., Weiher, E., and Westoby, M. (2006). Rebuilding community ecology from functional traits. *Trends Ecol. Evol.* 21, 178–185. doi: 10.1016/j.tree.2006.02.002

Mesfin, S., Gebresamuel, G., Haile, M., Zenebe, A., and Desta, G. (2020). Mineral fertilizer demand for optimum biological nitrogen fixation and yield potentials of legumes in Northern Ethiopia. *Sustainability* 12:6449. doi: 10.3390/su12166449

Moot, D. J., Scott, W. R., Roy, A. M., and Nicholls, A. C. (2000). Base temperature and thermal time requirements for germination and emergence of temperate pasture species. *N. Z. J. Agric. Res.* 43, 15–25. doi: 10.1080/00288233.2000.9513404

R Core Team (2018). *R: The R Project for Statistical Computing*. Available online at: https://www.r-project.org/ (accessed January 7, 2020).

Roscher, C., Schumacher, J., Schmid, B., and Schulze, E.-D. (2015). Contrasting effects of intraspecific trait variation on trait-based niches and performance of legumes in plant mixtures. *PLoS One* 10:e0119786. doi: 10.1371/journal.pone.0119786

Schiltz, S., Munier-Jolain, N., Jeudy, C., Burstin, J., and Salon, C. (2005). Dynamics of exogenous nitrogen partitioning and nitrogen remobilization from vegetative organs in pea revealed by ¹⁵N in vivo labeling throughout seed filling. *Plant Physiol.* 137, 1463–1473. doi: 10.1104/pp.104.056713

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH image to ImageJ: 25 years of image analysis. *Nat. Methods* 9, 671–675. doi: 10.1038/nmeth.2089

Shearer, G., and Kohl, D. H. (1986). N2-fixation in field settings: estimations based on natural ¹⁵N abundance. *Funct. Plant Biol.* 13, 699–756. doi: 10.1071/pp9860699

Siebenkäs, A., Schumacher, J., and Roscher, C. (2015). Phenotypic plasticity to light and nutrient availability alters functional trait ranking across eight perennial grassland species. *AoB Plants* 7:lv029. doi: 10.1093/aobpla/plv029

Siefert, A., Violle, C., Chalmandrier, L., Albert, C. H., Taudiere, A., Fajardo, A., et al. (2015). A global meta-analysis of the relative extent of intraspecific trait variation in plant communities. *Ecol. Lett.* 18, 1406–1419. doi: 10.1111/ele.12508

Singh, B. B., Ajeigbe, H. A., Tarawali, S. A., Fernandez-Rivera, S., and Abubakar, M. (2003). Improving the production and utilization of cowpea as food and fodder. *Field Crops Res.* 84, 169–177. doi: 10.1016/S0378-4290(03)00148-5

Stacklies, W., Redestig, H., Scholz, M., Walther, D., and Selbig, J. (2007). pcaMethods a bioconductor package providing PCA methods for incomplete data. *Bioinformatics* 23, 1164–1167. doi: 10.1093/bioinformatics/btm069

Tharanathan, R. N., and Mahadevamma, S. (2003). Grain legumes – a boon to human nutrition. *Trends Food Sci. Technol.* 14, 507–518. doi: 10.1016/j.tifs.2003.07.002

Therneau, T., Atkinson, B., Ripley, B., and Ripley, M. B. (2015). *rpart: Recursive Partitioning and Regression Trees*. R Package Version 41-10.

Tipping, M. E., and Bishop, C. M. (1999). Probabilistic principal component analysis. *J. R. Stat. Soc. B* 61, 611–622.

Tiwari, B. K., Gowen, A., and McKenna. B. M. (2011). *Pulse Foods: Processing, Quality and Nutraceutical Applications*. Amsterdam: Academic Press.

Tribouillois, H., Fort, F., Cruz, P., Charles, R., Flores, O., Garnier, E., et al. (2015). A functional characterisation of a wide range of cover crop species: growth and nitrogen acquisition rates, leaf traits and ecological strategies. *PLoS One* 10:e0122156. doi: 10.1371/journal.pone.0122156

Unkovich, M., Herridge, D., Peoples, M., Cadisch, G., Boddey, B., Giller, K., et al. (2008). *Measuring Plant-Associated Nitrogen Fixation in Agricultural Systems*. Canberra, ACT: Australian Centre for International Agricultural Research (ACIAR).

Violle, C., Navas, M. L., Vile, D., Kazakou, E., Fortunel, C., Hummel, I., et al. (2007). Let the concept of trait be functional! *Oikos* 116, 882–892. doi: 10.1111/j.2007.0030-1299.15559.x

Westoby, M. (1998). A leaf-height-seed (LHS) plant ecology strategy scheme. *Plant Soil* 199, 213–227. doi: 10.1023/A:1004327224729

Wood, S. A., Karp, D. S., DeClerck, F., Kremen, C., Naeem, S., and Palm, C. A. (2015). Functional traits in agriculture: agrobiodiversity and ecosystem services. *Trends Ecol. Evol.* 30, 531–539. doi: 10.1016/j.tree.2015.06.013

Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. *Nature* 428, 821–827. doi: 10.1038/nature02403

Zhang, W., Ricketts, T. H., Kremen, C., Carney, K., and Swinton, S. M. (2007). Ecosystem services and dis-services to agriculture. *Ecol. Econ.* 64, 253–260. doi: 10.1016/j.ecolecon.2007.02.024

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Guiguitant, Vile and Marrou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
# APPENDIX

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