Next Article in Journal
Advanced Biomethane Production from Biologically Pretreated Giant Reed under Different Harvest Times
Previous Article in Journal
Pseudomonas spp. Producing Antimicrobial Compounds Regulate Fungal Communities Inhabiting Wheat Crown in Southern Chile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Destructive Estimation of the Leaf Weight and Leaf Area in Common Bean

1
Programa de Ingeniería Agroecológica, Facultad de Ingeniería, Universidad de la Amazonia, Florencia 18001, Colombia
2
Centro de Investigaciones Amazónicas CIMAZ Macagual Cesar Augusto Estrada González, Grupo de Investigaciones Agroecosistemas y Conservación en Bosques Amazónicos—GAIA, Florencia 18001, Colombia
3
CATIE—Centro Agronómico Tropical de Investigación y Enseñanza, Turrialba 30501, Costa Rica
4
Facultad de Ciencias Agropecuarias, Universidad Nacional de Córdoba, Cordoba 5016, Argentina
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(3), 711; https://doi.org/10.3390/agronomy12030711
Submission received: 6 December 2021 / Revised: 27 February 2022 / Accepted: 6 March 2022 / Published: 16 March 2022

Abstract

:
Regression models to predict leaf area and leaf weight in common bean (Phaseolus vulgaris) were fitted using the three leaflets of the leaves. A total of 1504 leaves from 40 genotypes were collected, covering a large range of leaf sizes. Width, length, area, and weight were measured for each leaflet. The total leaf area and weight was obtained by the sum of left, central, and right leaflets. The dataset was randomly divided into training and validation sets. The training set was used for model fitting and selection, and the validation dataset was used to obtain statistics for model prediction ability. The leaf area and leaf weight were modeled using different linear regression models based on the length and width of the leaflet. Polynomial regressions involving both length and width of the leaflet provided very good models to estimate the expected area (R2 = 0.978) and weight (R2 = 0.820) of leaves.

1. Introduction

One of the main crops in the diet of the tropics of Latin America and East Africa is the common bean (Phaseolus vulgaris L.). The common bean is produced by smallholders in areas with high climatic vulnerability. Due to their importance, there are currently many issues in relation to the adaptation of improved varieties to different environmental conditions, particularly in the face of imminent challenges posed by global climate change. Leaf area is crucial to evaluating water stress [1], as a larger leaf area leads to greater hydric conductance [2].
To determine the adaptation of bean materials to different conditions, growth should be evaluated and simulated using mathematical models that describe the relationship between plant growth, dry matter production, and expansion of total leaf area. One of the ways to evaluate the interactions between environmental conditions and genotypes is the area and weight of the leaf, which are considered indicators of growth and yield in crops.
Determining leaf area and weight requires the destructive sampling of many leaves, a process that can be time-consuming and labor-intensive, as well as have an impact on plant growth. Regression models allow the prediction of leaf area and weight using variables that can be measured without destructive sampling techniques. In addition, this process allows the repeated measurement of leaves over time, while avoiding the biological alteration characteristic of destructive methods. Phaseolus vulgaris leaves are trifoliate. Therefore, the length and width of each leaflet can be used to estimate the leaf total area and weight.
Several models have been reported to predict leaf area in legumes such as Alfalfa (Medicago sativa) [3], chestnut (Castanea sativa) [4], chickpea (Cicer arietinum L.) [5], cowpea (Vigna unguiculata L.) [6], faba bean (Vicia faba L.) [7,8], peanut (Arachis hypogaea L.) [9,10], pigeon pea (Cajanus cajan L.) [11], and soybean (Glycine max) [12,13,14], where equations were developed from the length or width of some of the leaflets; however, no models are reported to predict leaf weight and area from leaflet measurements for a large number of genotypes of Phaseolus vulgaris.
Specifically for the common bean, a few models have been generated to predict leaf area. The first linear model was reported by Lakitan [15], whose objective was to predict the area of the central and lateral leaflets for three bean genotypes from the length and width of a leaflet. Muñoz et al. [16] proposed to measure leaf area by multiplying the length by the width of the leaflet corrected by a correction factor; then, Iamauti [17] and De Jesus et al. [18] mentioned an empirical model from leaflet length, a situation that was also reported by Figueiredo et al. [19] for the cultivar Perola type Carioca, Toebe et al. [20] for cultivar “Macarrão” and Queiroga et al. [21] for cultivar UEL-2, but these last two specify only the use of the central leaflet to predict leaf area. Bhatt and Chanda [22] reported different models for bean cultivation under different stress conditions. Likewise, Ramírez-Builes et al. [23] reported models for different bean genotypes (BAT 477, Morales, SER 16 and SER 21), some of them resistant to water-stress conditions, based on length, width, and its product, as well as Hara et al. [24]. Finally, Lakitan et al. [25] proposed a simplified procedure for the non-destructive estimation of the trifoliate leaf area using the length and width of the leaflets. Therefore, the objective of this study was to determine an empirical model to estimate the leaf area and weight of different lines of the common bean from the linear dimensions of the leaflets, with an emphasis on 40 advanced modern lines developed with the capacity to adapt to drought, high temperatures, low soil fertility, aluminum toxicity, resistance to pests, and high nutrient content. These lines were obtained by crossing P. vulgaris with other species of genus Phaseolus sp.

2. Materials and Methods

2.1. Plant Material and Experimental Design

Evaluation of leaflet morphometric measurements in bean genotypes was carried out in the Centro de Investigaciones Amazonicas CIMAZ Macagual, Universidad de la Amazonia, Colombia (1°37′ N and 75°36′ W) located in Florencia, Caquetá (Colombia) in a tropical rainforest ecosystem. The location has an average annual precipitation of 3800 mm with 1700 h of sunshine per year, an average temperature of 25.5 °C, and an average relative humidity of 84%. A total of 40 improved lines of common bean from the Mesoamerican and Andean gene pools were used developed by CIAT’s Bean Breeding Program: 18 advanced lines of Phaseolus vulgaris (ALB 349, AMADEUS, BFS 10, BFS 29, ICA QUIMBAYA, ICTA PETEN, SCR 23, SCR 40, SCR 61, SEN 48, SEN 70, SER 323, SER 324, SMC 234, SMN 98, SMR 181, SMR 182, SMR 185), 3 interspecific lines from P. vulgaris × P. coccineus which resulted in a group of Mesoamerican gene pool lines (ALB 350, ALB 352, ALB 60) and Andean (ALB 267, RRA 57), 2 interspecific lines from P. vulgaris × P. acutifolius (SER 213, SIN 351-1), 5 interspecific lines from P. vulgaris × P. dumosus (SMC 101, SMC 33, SMC 41, SMN 24, SMR 84), 4 interspecific lines from P. vulgaris × P. acutifolius × P. coccineus (SEF 1, SEF 10, SEF 49, SEF 70, SEF 73), 1 interspecific lines from P. vulgaris × P. acutifolius × P. dumosus (SMR 195), 4 interspecific lines from P. vulgaris × P. acutifolius × P. dumosus × P. coccineus (SMC 199, SMC 209, SMC 216, SMR 140).
These selected genotypes were resistant to different abiotic and biotic stress factors, high micronutrient content, and other desirable attributes that enhance adaptation and grain yield. The ALB (small red kidney, black kidney) lines have enhanced adaptation to drought and Al toxicity. AMADEUS is a high-temperature tolerant improved line, the BFS lines (small red) have enhanced adaptation to low soil fertility, RRA (various colored) were developed to generate resistance to Andean root rot, genus Pythium sp. and Sclerotium sp. The SEF (red), SEN (black) and SER (small red) lines were improved to adapt to drought and heat. The SIN (black) lines were interspecific with drought tolerance. SCR (red), SMC (colored), SMR (red) and SMN (black) lines are tolerant to drought with a high mineral content in seeds (Fe). ICA QUIMBAYA (red) is a commercial variety resistant to bean golden mosaic virus and Al toxicity, and ICTA PETEN (black) is a commercial variety, it has tolerance to rust and golden mosaic virus with high Fe content.
Approximately twenty plants per plot were sown in a greenhouse for each genotype, and the plots were arranged in a nonreplicated complete randomized design. The same substrate was used in each plot. The plot size was 2 m × 2 m, with three furrows at 0.6 m between them. For each plot, three plants were selected from the middle part of the central furrow to control border effects for a total of 120 plants (3 plants × 40 genotypes). For each plant, between 4 and 19 leaves (depending on the plant size) were systematically selected (one out of three) from the 7th leaf developed from the base to the top of the plant in the growth stage of pod filling (R8), which corresponds to a period of 70 to 80 days after sowing. To predict leaf area and weight, we measured the length and width of each leaflet to see if the morphological relationship between them remains constant for leaves of different ages.
We obtained the measurements of morphometric variables for a total of 4512 leaflets of 1504 leaves. Each leaflet (left, right, and center) on each leaf was weighed independently using an Ohaus Scout electronic scale (100 ± 0.001 g), after which it was scanned using the HP ScanJet Pro 2500 scanner. The length, width, and area of each leaflet was measured from its scanned image using the ImageJ program V1 [26]. From this information, we calculated descriptive statistics for each of the morphometric variables evaluated, as well as calculated the Pearson’s correlation coefficients between morphometric variables.

2.2. Data Analysis

Only 28 genotypes from a total of 40 had enough observations (n ≥ 14) to fit their own complete polynomial model (Model 1, see below). Nevertheless, the observations of these 28 genotypes (n = 1398) represented 93.2% of the complete dataset (n = 1500). To compare the performance of genotype-specific models against a unique one, we used both strategies to obtain predictions. Then, we calculated correlation concordance coefficients for the predictions under the two strategies. For total leaf area, the correlation concordance coefficient was 0.99, whereas for total leaf weight it was 0.93. According to these preliminary results, we decided to fit unique models for all genotypes which are more suitable for practical applications.
We used multiple linear regression models to predict leaf area and weight from leaflets length and width. Due to the curvilinear trends between the dependent (leaf weight and area) and predictor variables (leaflet length and width), multiple linear regression models including second-order polynomials on length, width or both variables, and the interaction (length × width) were used to estimate the expected value of leaflet weight and area. The best model was selected using AIC and BIC criteria. The smaller the better, for both criteria. Because none of them were uniformly the best, we decided to choose a model in which both criteria agree.
We evaluated five different models according to the morphometric measures used (c = central, r = right, and l = left leaflets; W = width, WW = width2, L = length, LL = length2 and LW = length × width):
Model 1: to estimate leaf area and weight using all the leaflet data:
Yi = β0 + β1 cLi + β2 cWi + β3 cLLi + β4 cWWi + β5 lLi + β6 lWi + β7 lLLi + β8 lWWi + β9 rLi + β10 rWi + β11 rLLi + β12 rWWi + ei
Model 2: to estimate leaf area and weight using the central leaflet data:
Yi = β0 + β1 cLi + β2 cWi + β3 cLLi + β4 cWWi + ei
Model 3: to estimate leaf area and weight using the left leaflet data:
Yi = β0 + β1 lLi + β2 lWi + β3 lLLi + β4 lWWi + ei
Model 4: to estimate leaf area and weight using the right leaflet data:
Yi = β0 + β1 rLi + β2 rWi + β3 rLLi + β4 rWWi + ei
Model 5: to estimate leaf area and weight using only the leaflet width data:
Yi = β0 + β1 cWi + β2 cWWi + β3 lWi + β4 lWWi + β5 rWi + β6 rWWi + ei
For all models, ei (i = 1, …, 1200) is the error term assumed to be normal and independently distributed, with zero mean and a common variance.
The models using only length information from the three leaflets were not reported because they were much less predictive than those based only on width.
From the total of the 1504 leaves recorded, we obtained 1500 valid data values. We randomly selected 1200 leaves to conform the training set used to estimate the models. The remaining 300 leaves were used as a validation set to measure the predictive ability of a fitted model, calculating, from this set, the coefficient of determination (R2) and prediction mean square error (PMSE). We used the adjusted R2 in multiple linear regression models and R2 in the case of simple regression models. Adjusted R2 (or R2) and PMSE are summary statistics commonly used for the comparison of models. However, they do not say anything about patterns of departure of observed values against predicted. A simple and effective way to visualize any trouble regarding this issue is to draw a scatter plot of observed vs. predicted values and overlap a reference line (y = x) as well as the regression line of observed vs. predicted values. Departures of the regression line from the reference line can suggest problems of over or underestimation [27]. All the statistical analyses were done in InfoStat version 2021 [28] and R version 3.6.1 [29].

3. Results

3.1. Summary Measures of Morphometric Variables Evaluated in the Different Bean Lines

The summary statistics are shown for each of the morphometric variables evaluated (Table 1). The maximum and minimum statistics of the regressors are useful to approximate the range of prediction space in which the prediction model is reliable. Pairwise Pearson’s correlation coefficients were also calculated among the set of morphometric variables (Figure 1). The correlation between leaf total area and the leaflet measurements (Figure 1) was higher than the correlation between leaf total weight and leaflet measurements (Figure 2) (ranged between 0.886–0.931 and 0.78–0.86, respectively). For leaf total area, the width measurements of the leaflets were more correlated than those of leaflet length (ranged between 0.886–0.9 and 0.91–0.93, respectively). The leaf total weight and the width measurements of the leaflets were more correlated than with those of leaflet length (ranged between 0.78–0.826 and 0.804–0.813, respectively).

3.2. Regression Models

In all the models, evaluating the inclusion of interaction (length × width) was not necessary (p > 0.10), but the second-order polynomial term on leaflet length and leaflet width was included. In addition, the inclusion of an interaction term in any model increases both AIC and BIC values. Regression coefficients for the five models described above were estimated for both dependent variables: weight and leaf area (Table 2). In total area estimation, Model 1, which is a second-order polynomial on all leaflet length and leaflet width, without their interaction, has the lowest PMSE and the greatest adjusted R2 (90.95 and 0.9781, respectively) and has the least values for both AIC and BIC. The models to estimate the total area from the individual leaflets (Models 2 to 4) have greater PMSE and lesser adjusted R2 (ranged between 345.65 and 368.13, and ranged between 0.9108 and 0.9154, respectively). The model to estimate leaf area using only the width information (Model 5) has a PMSE = 177.72 and adjusted R2 = 0.9571 (Table 2).
In total weight estimation, Model 1 has the lowest PMSE and the greater adjusted R2 (0.2033 and 0.8205, respectively) and has the least values for both AIC and BIC. The models to estimate the total weight from the individual leaflets (Models 2 to 4) have the greatest PMSE and lower adjusted R2 (PMSE ranged between 0.2160 and 0.3412, and R2 ranged between 0.6988 and 0.8092, respectively). The model to estimate leaf weight using only the width information (Model 5) has a PMSE = 0.2142 and adjusted R2 = 0.8109 (Table 2).
In order to evaluate the prediction quality, we used scatter plots of observed against predicted values, adding the reference line y = x. A departure of the regression line of observed vs. predicted from the reference line is evidence of bias. For leaf area estimation (Figure 3), the model with the least bias is Model 2. However, because it only considers the central leaflet, its R2 is lower than the R2 in Model 1. Model 1 shows a little departure from the reference line for large values of area, showing a little underestimation. Model 4 is similar in bias to Model 1 but shows a smaller R2. Finally, Model 5 exhibits a similar behavior to Model 1, also underestimating for large-area leaves but showing a high R2 in both cases. When modeling the area and considering the PMSE (smaller is better), these two models have lower PMSE (90.95 and 177.72, respectively) and higher R2 than the other models. For this reason, we recommend using Model 1 or Model 5, but it is important to consider that Model 5 only needs width measures of leaflets (more parsimonious model).
For leaf weight estimation (Figure 4), the model with the least bias is Model 4. However, because it only considers the right leaflet, its R2 is lower than the Model 1 R2 (considering the three leaflets). In the same way as for area estimation, Model 1 shows a little departure from the reference line for large values of weight, showing a little underestimation. Model 2 and Model 3 show a little bias for large weight values but underestimate and overestimate, respectively. Finally, Model 5 shows a similar behavior to Model 1, also underestimating for large leaf weights but showing a high R2 in both cases (0.81 and 0.82, respectively). Considering the PMSE, these two models have lower PMSE (0.20 for Model 1 and 0.21 for Model 5) and higher R2 than the other models. In weight estimation, Model 5 has almost the same PMSE and R2 as Model 1, but only uses width measures of leaflets (more parsimonious model).

4. Discussion

The estimation of the leaf total area and weight by non-destructive methods is widely used in studies of physiology (photosynthetic capacity) and management intensity (fertilization levels, water availability), among others [30,31,32,33]. The use of destructive methods has a disadvantage when the aim is to evaluate growth, because the use of different leaves through time increases the experimental error [21,34].
Some research proposes statistical models to estimate leaf area and weight in legumes [3,4,5,6,7,8,9,10,11,12,13,14]. In all cases, the linear regressions used involve measurements of length and width (L and W) of one or more leaflets of the leaf. Furthermore, due to the non-linear relationship between these measurements and leaf area or weight, their squares (LL and WW) or their interactions (LW) are incorporated (e.g., Lakitan et al. [25]).
In the case of Phaseolus vulgaris, the leaf is trifoliate, therefore, the length and width of each leaflet can be obtained to estimate the leaf’s total area and weight. A few models have been proposed to estimate the leaf area in beans using linear regressions [20,21,22,25]. These models were made using one, two or three cultivars. The major novelty of our research lies in the dataset covering a large variety of genotypes and leaf sizes and ages, which makes it a general-use model. In addition, we propose models to estimate leaf weight based on L and W measurements of the leaflets. In our models, estimates were made using a training dataset and a validation dataset, which allowed us to obtain independent measurements of the PMSE and R2 statistics [35]. Among the models found for beans, only Lakitan et al. [25] and Queiroga et al. [21] used a validation dataset.
The number of leaves measured in this research (1500) far exceeds the sample size used by other authors. Due to the number of intervening genotypes, compared to other research, the range of variation of the length and width measurements of each leaflet was much greater, increasing the range of validity of the models that we propose in this research. The minimum and maximum W (1.02 cm; 13.34 cm) and the minimum and maximum L (2.40 cm; 17.88 cm) in this study have both a larger range than reported by Lakitan et al. [25] (W: 0.9 cm; 9 cm, L: 1.5 cm; 15.5 cm). The range of regressor values should always be reported [36], however, sometimes the ranges are not provided (e.g., Queiroga et al. [21]).
Several studies have determined that some species present more regularity in the area of the central leaf compared with the lateral ones, finding that there is more correlation between length and width of the central leaf with respect to the other two [21]. In our research, the correlation between L and W for central leaflet was 0.87, while the correlation between L and W for right and left leaflets was 0.90 in both cases, showing a similar behavior for the shape of the three leaflets. This different result may be due to the fact that in our study, there were 40 genotypes instead of one, as with Queiroga et al. [21]. Marshall [37] and Chanda et al. [38] point out that this correlation can change according to the plant age and environmental conditions.
The best models in our research do not include the LW interaction, but they include the second-order polynomial term. Lakitan et al. [25] developed models using L, W, and LW. However, they use zero intercept in their models, because for them “it makes sense to force intercept to zero, since if L, W or LW is zero, then leaf area should also be zero”. This is a common error in regression models. When the zero-intercept regression model is fitted, it should not be based only on conceptual issues but mainly on the data. If there are not near-zero-regressor values available, when assuming zero intercept, we are also assuming the linearity of the model in a range of values were the regressors were not observed. Forcing the model to have zero-intercept because of a non-sense estimated intercept (for example negative area) leads, at the same time, to miss-fitting the model for observed values.
Pompelli et al. [39] and Lakitan et al. [25] found that the power model predicted the leaf area with good accuracy but increasing heteroscedastic residual dispersion compared with the use of LW. We could not verify such a behavior.
We obtained the best model to estimate leaf area and weight using all the leaflet information (L and W) and its second-order polynomial (LL and WW, Model 1). These models showed only a little underestimation for large leaf area and weight values. In addition, estimates of leaf total area and weight are provided based on the length and width measurements of leaflets, achieving that with only two measurements, there are models with R2 greater than 0.9 in the case of total area. In addition, for total area estimation, the models using the central, left, or right leaflet measurements are very similar in terms of their PMSE and R2.
For weight, there is a marked difference between the estimates obtained from the individual leaflets. Those from the length and width of the right leaflet shows the highest R2 and lowest PMSE than the central and left leaflets. The correlation between leaf area and weight with W measurements was higher than the correlation with L measures for all leaflets. Queiroga et al. [21] estimated the leaf area by measuring only the width of the central leaflet. Using a logarithmic regression with zero intercept they obtained a good percentage of explained variability (R2 = 0.98), but they do not evaluate the presence of under or overestimation of their models. The selection of the best model must be based not only on R2 and PMSE but also use some measure of accuracy [27].
As in Queiroga et al. [21], in our research, we found a good estimation of leaf total area and weight only based on the measurement of the width of the three leaflets. Other authors found more correlation between width and leaf area than length (for instance, Zaffaroni [40] in Vigna unguiculata L. and Bianco [41] in Amarantus retroflexus L.). The width of the leaflets on the leaf is easier to measure than the length, due, among other things, to the boundary between the lamina and the petiole of the leaflet. The results using only the width of the leaflets gave estimates with R2 values very similar to those corresponding to the use of the two dimensions for both total area and weight, but unlike the area estimation, the weight showed a very similar PMSE between Models 1 and 5, meaning that only measuring the width of the leaflets is better in the case of leaf weight estimation. We also evaluated a model using only central W as in Queiroga et al. [21], but the estimated model (Area = 1.19 + 5.86 W + 1.66 WW) showed an R2 = 0.86. This result is a consequence of the number of genotypes in our study, and the corresponding differences between W in the central, leaf, and right leaflets are because of the genotype effect.

5. Conclusions

Second-order polynomial regression involving both the length and width of each leaflet provided a very good fit to estimate the expected area and weight of leaves in Phaseolus vulgaris. However, it is possible to estimate total leaf area and weight with good accuracy using a second-order polynomial regression involving only the width of each leaflet. The loss of accuracy due to the use of this simplified model is less when estimating weight than when estimating area. Because the purpose of the predictive models we are presenting is not specific to a given genotype, using any model depending on only one leaflet is less accurate than a model using all of them.

Author Contributions

Conceptualization, J.C.S. and F.C.; methodology, J.C.S., F.C. and J.D.R.; software, J.C.S., F.C. and J.D.R.; validation, J.C.S., F.C. and J.D.R.; formal analysis, J.C.S., F.C. and J.D.R.; investigation, J.C.S.; resources, J.C.S.; data curation, J.C.S. and F.C.; writing—original draft preparation, J.C.S. and F.C.; writing—review and editing, J.C.S. and F.C.; visualization, J.C.S., F.C. and J.D.R.; supervision, J.C.S.; project administration, J.C.S.; funding acquisition, J.C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the authors upon request.

Acknowledgments

We acknowledge the financial support from the CGIAR Research Program on Grain Legumes and Dryland Cereals for the development of breeding lines of common bean. We would also like to thank all donors who supported this work through their contributions to the CGIAR Fund.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, C.; Zhang, J.; Zhang, H.; Zhao, J.; Wu, Q.; Zhao, Z.; Cai, T. Mechanisms for the relationships between water-use efficiency and carbon isotope composition and specific leaf area of maize (Zea mays L.) under water stress. Plant Growth Regul. 2015, 77, 233–243. Available online: https://link.springer.com/article/10.1007/s10725-015-0056-8 (accessed on 12 November 2021). [CrossRef]
  2. Sellin, A.; Õunapuu-Pikas, E.; Kaurilind, E.; Alber, M. Size-dependent variability of leaf and shoot hydraulic conductance in silver birch. Trees 2011, 26, 821–831. Available online: https://link.springer.com/content/pdf/10.1007/s00468-011-0656-5.pdf (accessed on 12 November 2021). [CrossRef]
  3. Bianco, S.; De Carvalho, L.B.; Bianco, M.S. Determination of alfalfa leaf area by non-destructive method. Commun. Plant Sci. 2011, 1, 17–20. Available online: https://www.researchgate.net/publication/224884475 (accessed on 12 November 2021). [CrossRef]
  4. Serdar, Ü.; Demirsoy, H. Non-destructive leaf area estimation in chestnut. Sci. Hortic. 2006, 108, 227–230. [Google Scholar] [CrossRef]
  5. Nahas, S.B.; Arce, O.E.A.; Ricci, M.; Romero, E.R. Leaf area estimation of individual leaf and whole plant of chickpea (Cicer arietinum L.) by means of regression methods. Rev. Agron. Noroeste Argent. 2019, 39, 99–106. Available online: https://www.researchgate.net/publication/338867825 (accessed on 12 November 2021).
  6. de Oliveira, R.L.L.; Moreira, A.R.; da Costa, A.V.A.; de Souza, L.C.; Lima, L.G.S.; da Silva, R.T.L. Modelos de determinação não destrutiva de área foliar de feijão caupi Vigna unguiculata (L.). Glob. Sci. Technol. 2015, 8, 17–27. [Google Scholar] [CrossRef]
  7. Erdoğan, C. A Leaf Area Estimation Model for Faba Bean (Vicia faba L.) Grown in the Mediterranean Type of Climate. SDU J. Fac. Agric. Ziraat Fakültesi Derg. 2012, 7, 58–63. [Google Scholar]
  8. Peksen, E. Non-destructive leaf area estimation model for faba bean (Vicia faba L.). Sci. Hortic. 2007, 113, 322–328. [Google Scholar] [CrossRef]
  9. Ma, L.; Gardner, F.P.; Selamat, A. Estimation of Leaf Area from Leaf and Total Mass Measurements in Peanut. Crop Sci. 1992, 32, 467–471. Available online: https://acsess.onlinelibrary.wiley.com/doi/full/10.2135/cropsci1992.0011183X003200020036x (accessed on 14 November 2021). [CrossRef]
  10. Kathirvelan, P.; Kalaiselvan, P. Groundnut (Arachis hypogaea L.) Leaf Area Estimation Using Allometric Model. Res. J. Agric. Biol. Sci. 2007, 3, 59–61. [Google Scholar]
  11. Pezzini, R.V.; Filho, A.C.; Alves, B.M.; Follmann, D.N.; Kleinpaul, J.A.; Wartha, C.A.; Silveira, D.L. Models for leaf area estimation in dwarf pigeon pea by leaf dimensions. Bragantia 2018, 77, 221–229. Available online: http://dx.doi.org/10.1590/1678-4499.2017106 (accessed on 19 November 2021). [CrossRef] [Green Version]
  12. Nehbandani, A.; Soltani, A.; Zeinali, E.; Raeisi, S.; Najafi, A. Allometric relationships between leaf area and vegetative characteristics in soybean. Int. J. Agric. Crop. Sci. 2013, 6, 1127–1136. [Google Scholar]
  13. Bakhshandeh, E.; Kamkar, B.; Tsialtas, J.T. Application of linear models for estimation of leaf area in soybean [Glycine max (L.) Merr]. Photosynth 2011, 49, 405. [Google Scholar] [CrossRef]
  14. Bakhshandeh, E.; Ghadiryan, R.; Kamkar, B. A rapid and non-destructive method to determine the leaflet, trifoliate and total leaf area of soybean. Aust. J. Plant Sci. Biotechnol. 2010, 4, 19–23. [Google Scholar]
  15. Lakitan, B. Statistical Design and Analysis for Intercropping Experiments. Annu. Rep. Bean. Improv. Coop. 1989, 32, 19–21. [Google Scholar]
  16. Muñoz, G.; Giraldo, G.; Fernández de Soto, J. Descriptores Varietales: Arroz, Fríjol, Maíz, Sorgo; Publication No. 177; Ciat: Cali, Colombia, 1993. [Google Scholar]
  17. Iamauti, M.T. Avaliação de Danos Causados por Uromyces appendiculatus No Feijoeiro. Ph.D. Thesis, Universidade de São Paulo, Piracicaba, Brazil, 1995. [Google Scholar]
  18. de Jesus, W.C.; Vale, F.X.R.D.; Coelho, R.R.; Costa, L.C. Comparison of two methods for estimating leaf area index on common bean. Agron. J. 2001, 93, 989–991. Available online: https://acsess.onlinelibrary.wiley.com/doi/full/10.2134/agronj2001.935989x (accessed on 1 November 2021). [CrossRef]
  19. da Silva Figueiredo, E.; dos Santos, M.E.; Garcia, A. Modelos de determinação não destrutivo da área foliar do feijoeiro comum (Phaseolus vulgaris L.). Nucleus 2012, 9, 79–84. Available online: https://dialnet.unirioja.es/servlet/articulo?codigo=3988133&info=resumen&idioma=POR (accessed on 19 November 2021). [CrossRef] [Green Version]
  20. Toebe, M.; Filho, A.C.; Loose, L.H.; Heldwein, A.B.; Zanon, A. Área foliar de feijão-vagem (Phaseolus vulgaris L.) em função de dimensões foliares. Semin. Agrar. 2012, 33 (Suppl. S1), 2491–2500. Available online: https://www.researchgate.net/publication/278121722 (accessed on 19 November 2021). [CrossRef] [Green Version]
  21. Queiroga, J.L.; Romano, E.D.U.; Souza, J.R.P.; Miglioranza, É. Estimativa da área foliar do feijão-vagem (Phaseolus vulgaris L.) por meio da largura máxima do folíolo central. Hortic. Bras. 2003, 21, 64–68. [Google Scholar] [CrossRef] [Green Version]
  22. Bhatt, M.; Chanda, S.V. Prediction of leaf area in Phaseolus vulgaris by non-destructive method. Bulg. J. Plant Physiol. 2003, 29, 96–100. [Google Scholar]
  23. Ramírez-Builes, V.H.; Porch, T.G.; Harmsen, E.W. Development of linear models for estimation of leaflet area in common bean (Phaseolus vulgahs L.). J. Agric. Univ. Puerto. Rico. 2008, 92, 171–182. Available online: https://pubag.nal.usda.gov/catalog/31098 (accessed on 22 November 2021).
  24. Hara, A.T.; Gonçalves, A.C.A.; Maller, A.; Hashiguti, H.T.; de Oliveira, J.M. Ajuste de modelo de predição de área foliar do feijoeiro em função de medidas lineares. Rev. Eng. Agric. 2019, 27, 179–186. Available online: https://www.researchgate.net/publication/332774019 (accessed on 4 November 2021). [CrossRef]
  25. Lakitan, B.; Widuri, L.I.; Meihana, M. Simplifying procedure for a non-destructive, inexpensive, yet accurate trifoliate leaf area estimation in snap bean (Phaseolus vulgaris). J. Appl. Hortic. 2017, 19, 15–21. Available online: http://horticultureresearch.net/jah/2017_19_1_15_21.PDF (accessed on 4 November 2021). [CrossRef]
  26. Ferreira, T.; Ferreira, T.; Rasb, W. ImageJ User Guide: IJ 1.46r. 2012. Available online: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.731.215 (accessed on 19 November 2021).
  27. Suárez, J.C.; Melgarejo, L.M.; Durán Bautista, E.H.; Di Rienzo, J.A.; Casanoves, F. Non-destructive estimation of the leaf weight and leaf area in cacao (Theobroma cacao L.). Sci. Hortic. 2018, 229, 19–24. [Google Scholar] [CrossRef]
  28. Di Rienzo, J.; Casanoves, F.; Balzarini, M.; Gonzalez, L.; Tablada, M.; Robledo, C.W. InfoStat versión 2021. Grupo InfoStat, FCA, Universidad Nacional de Córdoba, Argentina. Available online: http://www.infostat.com.ar.2021 (accessed on 10 November 2021).
  29. R Development Core Team. R: The R Project for Statistical Computing. Foundation for Statistical Computing, V., Austria. 2021. Available online: https://www.r-project.org/ (accessed on 9 November 2021).
  30. Blanco, F.F.; Folegatti, M.V. Leaf area for greenhouse cucumber under salinity and grafting. Sci. Agric. 2005, 62, 305–309. [Google Scholar] [CrossRef] [Green Version]
  31. De Swart, E.A.M.; Groenwold, R.; Kanne, H.J.; Stam, P.; Marcelis, L.F.M.; Voorrips, R.E. Non-destructive estimation of leaf area for different plant ages and accessions of Capsicum annuum L. J. Hortic. Sci. Biotechnol. 2004, 79, 764–770. Available online: https://www.tandfonline.com/doi/abs/10.1080/14620316.2004.11511840 (accessed on 19 November 2021). [CrossRef]
  32. Rouphael, Y.; Rivera, C.M.; Cardarelli, M.; Fanasca, S.; Colla, G. Leaf area estimation from linear measurements in zucchini plants of different ages. J. Hortic. Sci. Biotechnol. 2006, 81, 238–241. Available online: https://www.tandfonline.com/doi/abs/10.1080/14620316.2006.11512056 (accessed on 19 November 2021). [CrossRef]
  33. Lizaso, J.I.; Batchelor, W.D.; Westgate, M.E. A leaf area model to simulate cultivar-specific expansion and senescence of maize leaves. Field Crops Res. 2003, 80, 1–17. [Google Scholar] [CrossRef]
  34. Silva NFDa Ferreira, F.A.; Fontes, P.C.R.; Cardoso, A.A. Modelos para estimas a área foliar de abóbora por meio de medidas lineares. Revista Ceres. 1998, 45, 287–291. Available online: http://www.ceres.ufv.br/ojs/index.php/ceres/article/download/2488/482 (accessed on 19 November 2021).
  35. Draper, N.R.; Smith, H. Applied Regression Analysis Applied Regression Analysis; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2014; pp. 1–716. Available online: https://books.google.es/books?hl=es&lr=&id=d6NsDwAAQBAJ&oi=fnd&pg=PR13&dq=DRAPER,+N.+R.%3B+SMITH,+H.+A.+Applied+regression+analysis.+3.+ed.+New+York:+Wiley,+1998.+736+p&ots=Byo8lbj2PR&sig=z3y9N4qbsQqWP5OLAFfLXjtQ1Gg (accessed on 19 November 2021).
  36. Balzarini, M.; Di Rienzo, J.; Tablada, M.; Gonzalez, L.; Bruno, C.; Córdoba, M.; Casanoves, F. Estadística y Biometrías. Ilustraciones del Uso de Infostat en Problemas de Agronomía; Segunda Edición; Editorial Brujas: Córdoba, Argentina, 2012; 380p. [Google Scholar]
  37. Marshall, J.K. Methods for leaf area maeasurement of large and small leaf samples. Photosynthetica 1968, 2, 41–47. [Google Scholar]
  38. Chanda, S.; Joshi, A.; Vaishnav, P.; Singh, Y. Leaf area determination in pearl millet using linear measurements-area and matter-area relationships. Photosynthetica 1985, 19, 424–427. [Google Scholar]
  39. Pompelli, M.F.; Antunes, W.C.; Ferreira, D.T.R.G.; Cavalcante, P.G.S.; Wanderley-Filho, H.C.L.; Endres, L. Allometric models for non-destructive leaf area estimation of Jatropha curcas. Biomass Bioenergy 2012, 36, 77–85. [Google Scholar] [CrossRef]
  40. Zaffaroni, E. Estimativa da área foliar em feijão macassar (Vigna unguiculata (L.) Walp) a partir de medidas dos folíolos. Agropecuária Técnica 1981, 2, 80–85. [Google Scholar]
  41. Bianco, S.; Pitelli, R.; Pavani, M.; Silva, R. Estimativa da área foliar de plantas daninhas. XIII-Amaranthus retroflexus L. Ecossistema 1995, 20, 5–9. [Google Scholar]
Figure 1. Scatter plots matrix between leaf total area versus length and width of the central, right, and left leaflet of leaves in common bean. Validation observations in red, training observations in black. Corr: Pearson correlation coefficients.
Figure 1. Scatter plots matrix between leaf total area versus length and width of the central, right, and left leaflet of leaves in common bean. Validation observations in red, training observations in black. Corr: Pearson correlation coefficients.
Agronomy 12 00711 g001
Figure 2. Scatter plots matrix between leaf total weight versus length and width of the central, right, and left leaflet of leaves in common bean. Validation observations in red, training observations in black. Corr: Pearson correlation coefficients.
Figure 2. Scatter plots matrix between leaf total weight versus length and width of the central, right, and left leaflet of leaves in common bean. Validation observations in red, training observations in black. Corr: Pearson correlation coefficients.
Agronomy 12 00711 g002
Figure 3. Relationship between the observed leaf area and the predicted leaf total area of improved lines of common bean for the five evaluated models. Regression line of observed versus predicted values: red line, reference line (x = y): black line.
Figure 3. Relationship between the observed leaf area and the predicted leaf total area of improved lines of common bean for the five evaluated models. Regression line of observed versus predicted values: red line, reference line (x = y): black line.
Agronomy 12 00711 g003
Figure 4. Relationship between the observed weight and the predicted leaf total area of improved lines of common bean for the five evaluated models. Regression line of observed versus predicted values: red line, reference line (x = y): black line.
Figure 4. Relationship between the observed weight and the predicted leaf total area of improved lines of common bean for the five evaluated models. Regression line of observed versus predicted values: red line, reference line (x = y): black line.
Agronomy 12 00711 g004
Table 1. Summary statistics for morphometric variables of leaflets of improved lines of common bean. The sample was taken to cover a large range of leaf sizes.
Table 1. Summary statistics for morphometric variables of leaflets of improved lines of common bean. The sample was taken to cover a large range of leaf sizes.
VariableLeafletMeanStd. DevMinimumMaximum
Area (cm2)Right38.2420.621.70132.78
Central39.6520.092.33132.76
Left39.5020.582.48146.99
Total117.4059.348.94412.53
Weight (g)Right0.570.350.022.17
Central0.600.350.032.13
Left0.580.340.032.28
Total1.751.010.096.58
Length (cm)Right8.772.302.5116.84
Central9.132.312.5317.88
Left8.912.342.4017.14
Mean8.932.212.4816.94
Width (cm)Right6.271.891.0212.86
Central6.591.951.1812.65
Left6.401.931.0313.34
Mean6.421.841.2012.95
Table 2. Regression coefficients of the five models evaluated to predict leaf area and weight from leaflet measures of improved lines of common bean. Prediction mean square error (PMSE) and determination coefficient (adjusted R2) are calculated on the validation dataset.
Table 2. Regression coefficients of the five models evaluated to predict leaf area and weight from leaflet measures of improved lines of common bean. Prediction mean square error (PMSE) and determination coefficient (adjusted R2) are calculated on the validation dataset.
Leaf AreaLeaf Weight
ParameterCoefficienttp-ValueCoefficienttp-ValueModel
β01.3603850.390.69580.11510.710.47521. All the leaflets’ measures
β1 cL−1.217674−0.780.4375−0.0025−0.030.9730
β2 cW−1.851119−1.140.25400.03100.410.6800
β3 cLL0.1929132.480.0134−0.0004−0.110.9114
β4 cWW0.7199126.670.00000.00440.870.3828
β5 lL2.2645341.330.1848−0.1135−1.440.1512
β6 lW1.8585300.990.32090.17041.970.0496
β7 lLL0.0577700.660.51130.00812.000.0460
β8 lWW0.3510002.680.0075−0.0050−0.820.4113
β9 rL1.7325120.970.3303−0.0600−0.730.4670
β10 rW−3.661254−1.920.05530.01440.160.8706
β11 rLL0.0956021.030.30390.00771.790.0736
β12 rWW0.8456746.170.00000.01282.020.0439
R20.9781 0.8205
PMSE90.9593 0.2033
AIC8447.72 1492.26
BIC8518.09 1562.63
β013.6148362.160.03130.21871.200.23172. Only central leaflet
β1 cL−6.117789−2.630.0087−0.1971−2.930.0035
β2 cW6.7016262.590.00970.24653.290.0010
β3 cLL0.8668857.510.00000.01755.230.0000
β4 cWW0.8373294.710.00000.00360.700.4834
R20.9164 0.6988
PMSE345.6574 0.3412
AIC9930.04 1944.81
BIC9960.20 1974.97
β0−3.350078−0.520.60070.09720.570.56693. Only left leaflet
β1 lL2.2107640.820.4139−0.1692−2.360.0185
β2 lW2.0945040.670.50020.25753.120.0018
β3 lLL0.3196712.280.02300.01494.010.0001
β4 lWW1.3640226.090.00000.00560.950.3434
R20.9153 0.7545
PMSE350.0354 0.2784
AIC10,045.55 1858.97
BIC10,075.51 1889.13
β0−8.177508−1.420.1562−0.0085−0.050.95704. Only right leaflet
β1 rL5.8096522.390.0170−0.0901−1.360.1734
β2 rW−1.468532−0.540.58790.16792.280.0229
β3 rLL0.1528911.190.23480.01143.250.0012
β4 rWW1.6765788.440.00000.01332.460.0140
R20.9108 0.8092
PMSE368.1351 0.2160
AIC9755.55 1616.08
BIC9785.71 1646.24
β07.928129−7.160.0000−0.0813−4.670.00005. All the leaflets width
β1 cW1.8952573.540.00040.04582.470.0136
β2 rW−3.8454770.570.5674−0.03390.540.5915
β3 lW2.6343872.140.03260.05070.710.4752
β4 cWW0.54365110.160.00000.0022−0.100.9212
β5 rWW1.21242111.620.00000.02306.210.0000
β6 lWW0.6089819.140.00000.00762.920.0035
R20.9571 0.8109
PMSE177.72 0.2142
AIC9001.98 1556.72
BIC9042.19 1596.93
β0 = intercept, β1, …, β12 = correlation coefficients, c = central, r = right, and l = left leaflets; W = width, WW = width2, L = length, LL = length2 and LW = length × width, AIC = Akaike information criteria, BIC = Bayesian information criteria.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suárez, J.C.; Casanoves, F.; Di Rienzo, J. Non-Destructive Estimation of the Leaf Weight and Leaf Area in Common Bean. Agronomy 2022, 12, 711. https://doi.org/10.3390/agronomy12030711

AMA Style

Suárez JC, Casanoves F, Di Rienzo J. Non-Destructive Estimation of the Leaf Weight and Leaf Area in Common Bean. Agronomy. 2022; 12(3):711. https://doi.org/10.3390/agronomy12030711

Chicago/Turabian Style

Suárez, Juan Carlos, Fernando Casanoves, and Julio Di Rienzo. 2022. "Non-Destructive Estimation of the Leaf Weight and Leaf Area in Common Bean" Agronomy 12, no. 3: 711. https://doi.org/10.3390/agronomy12030711

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop