key: cord-0054293-avr1gqf8 authors: Hugalde, Inés P.; Agüero, Cecilia B.; Barrios-Masias, Felipe H.; Romero, Nina; Viet Nguyen, Andy; Riaz, Summaira; Piccoli, Patricia; McElrone, Andrew J.; Walker, M. Andrew; Vila, Hernán F. title: Modeling vegetative vigour in grapevine: unraveling underlying mechanisms date: 2020-12-22 journal: Heliyon DOI: 10.1016/j.heliyon.2020.e05708 sha: 82759243e281143e2bb6ee0a61bc700a2381fe41 doc_id: 54293 cord_uid: avr1gqf8 Mechanistic modeling constitutes a powerful tool to unravel complex biological phenomena. This study describes the construction of a mechanistic, dynamic model for grapevine plant growth and canopy biomass (vigor). To parametrize and validate the model, the progeny from a cross of Ramsey (Vitis champinii) × Riparia Gloire (V. riparia) was evaluated. Plants with different vigor were grown in a greenhouse during the summer of 2014 and 2015. One set of plants was grafted with Cabernet Sauvignon. Shoot growth rate (b), leaf area (LA), dry biomass, whole plant and root specific hydraulic conductance (k(H) and L(pr)), stomatal conductance (g(s)), and water potential (Ψ) were measured. Partitioning indices and specific leaf area (SLA) were calculated. The model includes an empirical fit of a purported seasonal pattern of bioactive GAs based on published seasonal evolutionary levels and reference values. The model provided a good fit of the experimental data, with R = 0.85. Simulation of single trait variations defined the individual effect of each variable on vigor determination. The model predicts, with acceptable accuracy, the vigor of a young plant through the measurement of L(pr) and SLA. The model also permits further understanding of the functional traits that govern vigor, and, ultimately, could be considered useful for growers, breeders and those studying climate change. Vigor is the propensity to assimilate, store, and/or use non-structural carbohydrates for producing large canopies, and it is associated with intense metabolism and fast shoot growth (Ollat et al., 2003; Rebolledo et al., 2015) . Understanding vigor in grapes is necessary to optimize vineyard management and grape breeding strategies, since a vine needs enough canopy and growth to ripen the grapes. For our purposes, vigor is represented by canopy biomass. Photosynthesis (A) is the ultimate source of carbohydrates that drives plant growth. For A to occur, CO2 must diffuse into the leaf mesophyll, through the stomata, which in turn results in H2O molecules escaping from the leaf to the atmosphere. This inevitable water loss through the stomata (and the depreciable diffusing through cuticle), constitutes transpiration (E). This means that A and stomatal conductance (gs) are tightly correlated (Wong et al., 1979) , and stomata are directly responsible for optimizing E vs. A (Rogiers et al., 2012) . At the cellular level, growth involves expansion and division (Polymenis and Schmidt, 1997) . Cell expansion occurs because cell walls are extensible, meaning they deform under the action of tensile forces, generally caused by turgor (Cosgrove, 2005) . Canopy hydration and turgor maintenance depend on the root system capacity for water uptake (hydraulic conductance, kH) to meet water demands (E). This means that root kH affects the hydration and turgor of the canopy (De Herralde et al., 2006; Lovisolo et al., 2008) , resulting in different levels of cellular extension and plant growth (Di Filippo and Vila, 2011) . Lower kH has been observed in shoots grafted onto dwarfing rootstocks (Syvertsen and Graham, 1985; Atkinson et al., 2003) and in non-grafted, low vigor, cherry plants (Zori c et al., 2012) . This can be explained by the equation of water potential components Ψp ¼ ΨL-Ψπ; and the Van den Honert equation, ΨL ¼ Ψsoil-(E/kH); that state that turgor water potential (Ψp) depends on leaf water potential (ΨL) and osmotic water potential (Ψπ), and that ΨL inversely depends on kH. Thus, Ψp can be directly impacted by changes in kH. In this way, ΨL increases with higher kH, and results in higher leaf turgor (Ψp) and cell expansion (Fan et al., 2012) . However, if kH decreases, as result of decreasing soil moisture and cavitation (Chaves et al., 2003) , or decrease in aquaporin activity (Ehlert et al., 2009; Vandeleur et al., 2009; Gambetta et al., 2012) ; a series of signals may reduce stomatal conductance in the leaves and decrease transpiration rates. According to the mechanistic model of stomatal closure of Buckley et al. (2003) , there is a hydraulic feedback signal between gs and kH through the regulation of the turgor potential of the guard cell, in response to variations in ΨL. This happens because the guard cells directly respond to the epidermal turgor potential (Ψp) (Chaves et al., 2003) , which in turn depends on plant kH, so reductions in kH can be translated into reductions in gs (Lovisolo et al., 2010; Zufferey and Smart, 2012; Pou et al., 2013) . Many hydraulic and metabolic processes are controlled or influenced by hormones, with gibberellins (GAs) being the primary hormones that control growth and development in plants (Huerta et al., 2008) . Gibberellins are tetracyclic diterpenoid acids involved in a number of developmental and physiological processes in plants (Crozier et al., 2000) . Among many effects, GAs activate cell division and control cell elongation (Huerta et al., 2008; Ubeda-Tom as et al., 2009 ), thus controlling shoot growth (amount and direction; Clúa et al., 1996) . They also play an important role in regulating net photosynthesis (Huerta et al., 2008) , chlorophyll degradation (Li et al., 2010) , and carbohydrate production and partitioning (Lambers et al., 1995 Moreno et al., 2011; Murcia et al., 2016) . However, from the þ130 of GA structures characterized up to now (http://www.plant-hormones.info/gainfo.asp), only C19 3-ß and 13-∞ hydroxylated chemical forms, GA1, GA3, GA4, GA7 are active per se (Kobayashi et al., 1994) . High purported concentrations of bioactive GAs have been proposed as responsible for heterosis in plants (Bate et al., 1988; Rood et al., 1988) and therefore to benefit canopy growth, while low concentrations of bioactive GAs benefit root biomass production. The concentration of bioactive GAs in each tissue depends on the organ analyzed, the plant phenology and on external and environmental stimuli. In grape leaves, the amount of GA19, a major GA not active per se but through its conversion to GA1 (Kobayashi et al., 1994) , is about 20 ng g-1 (Acheampong et al., 2015) . They may be site-synthesized by 3-beta oxidation from precursors of the inactive GA pool (Yamaguchi, 2008) , de-conjugated from glycosylated forms (Piccoli et al., 1997; Cassan et al., 2001) or may travel from the roots to the canopy (Hooijdonk et al., 2010) . In addition, GAs affect the biosynthesis and response of ethylene, and the synthesis and transport of auxins. A large number of GA mutants showing growth and development suppressions or dwarf phenotypes have demonstrated that the GA signaling complex is highly conserved in many species, like maize, rice, barley and grapevine (Ikeda et al., 2001; Boss and Thomas, 2002; Chandler et al., 2002; Muangprom and Osborn, 2004) . The surrounding environment also influences biomass production. Soil water status directly affects root hydraulics, while atmospheric humidity, water vapor pressure deficit, light and temperature will affect leaf water potential and stomatal conductance (Buckley et al., 2003; Buckley, 2005) . These responses to external factors will affect the water uptake capacity of the plant, leaf gas exchange (gs and A), hormonal regulation, tissue turgor and in consequence, plant growth. The complexity that characterizes biomass production makes it necessary to find a method that simplifies their study. In this sense, mechanistic models constitute a valuable tool for researchers and biophysicists to study diverse biological phenomena. In contrast with empirical approaches, mechanistic models are based on biophysical/ biochemical mechanisms. These models formalize a hypothesis, given that they are a possible description of the system and its performance. Each aspect of the performance of the dynamic model can be studied in any time frame. Once a model is validated, simulations are quick and inexpensive, whereas empirical models need to rely on experimental approaches (Ingalls, 2013) . To date, several biophysical models have been constructed, interpreting for instance photosynthesis and stomatal conductance (Farquhar et al., 1980; Farquhar and Sharkey, 1982; Pury and Farquhar, 1997; Yin and Struik, 2009; Bagley et al., 2015; Wu et al., 2016; Pleban et al., 2020) . Our paper describes the construction, parametrization and validation of a mechanistic, dynamic model for vigor using grapevine as the biological model. Examples of modeling in grapevine include studies of source sink relationships (Pallas et al., 2010) , stomatal function and cavitation (Hugalde et al., 2018) , rootstock control in scion transpiration , the dynamics of water transport , berry growth (Zhu et al., 2019) and the significance of changing temperatures for grapevine architecture (Schmidt et al., 2019) , among others. Our work considers that growth occurs when biomass is accumulated. This happens when CO2 is assimilated by photosynthesis. Later this biomass is partitioned into canopy and roots. The canopy biomass, by turgor, is then transformed in leaf area, while the root biomass will manage water economy and transport to the whole plant. Our final objective is to elucidate the functional traits that govern vigor, understand the involved mechanisms and, ultimately, hypothesize about the practical use of this knowledge for growers, breeders and ecologists specialized in climate change. This study is organized in four main sections. The first section, model construction, explains the theoretical development of the model based on the integration of a number of already existing models that explain several aspects of plant growth. The second section, which describes model parametrization and model validation, uses data obtained with a population from the cross Ramsey  Riparia to characterize the functions in the model and test its outputs. Ramsey (V. champinii) and Riparia Gloire de Montpellier (V. riparia) are known to confer to the scion high and low vigor, respectively and a population derived from their cross was developed to genetically characterize nematode resistance, lime and salt tolerance and ability to induce vigor (Lowe and Walker, 2006) . Principal Component Analysis (PCA), led the parameterization strategy used to define the input and output variables. PCA allows the identification of the variables that most closely represent the trait, by searching the total variance in the data to reduce the original variables into a smaller set of linear combinations. This technique is mostly used as a tool in exploratory data analysis and for making predictive models. It is often used to visualize genetic distance, relatedness between populations, and for finding patterns in data of high dimension (Jolliffe and Cadima, 2016; Li and Ralph, 2018) . The third section, model simulation and sensitivity analysis, explores the power of the model by predicting scenarios under different conditions. Finally, we discuss our results in the context of related works. This model conceives vigor as the result of the interaction of environment (soil and atmospheric conditions) and plant physiology ( Figure 1 ). It integrates a number of already existing models that explain several aspects of plant growth. Stomatal conductance (gs) is interpreted by the Buckley et al. (2003) model; leaf temperature (TL) is expressed by the Campbell and Norman (2012) equations; osmotic water potential of the guard cell (Ψπg) is empirically expressed as in Taiz et al. (2015) and Tardieu and Simonneau (1998) . Photosynthesis (A) is interpreted by the I.P. Hugalde et al. Heliyon 6 (2020) e05708 Farquhar model (Farquhar et al., 1980) , and biomass production is obtained by the stoichiometric equation that transforms CO2 in carbohydrates and deducts respiration (Salisbury and Ross, 2000) . Soil water status directly affects root hydraulics and plant water status. Root hydraulic conductance (Lp, Equation (1)) depends on specific root hydraulic conductance (Lpr,) and root biomass. Lp constitutes a proxy to plant hydraulic conductance under no cavitation events (kH, Equation (2)). Specific hydraulic conductance (kL, Equation (3)) is estimated from dividing kH by leaf area (LA). Atmospheric variables are formalized by the Campbell and Norman (2012) equations of atmospheric deficit (q 0 , Equation (4)) and guard cell osmotic potential (Ψπg, Equation (5)), which depend on temperature and solar radiation, respectively. Since the purpose of the model is to simulate growth day by day, Ψπg is considered constant and equal to -3 MPa. This value corresponds to the average daily value of the linear fit of Ψπg and osmolites (Kþ and sucrose) that control stomatal aperture during a complete day Zeiger, 1996, 1998) . Equation (4) expresses the difference in density of water vapor between air (ea) and leaf (es(TL)), which is the driving force for transpiration. These variables are then included in the Buckley (2005) equation for gs (eq. 6). An exhaustive list of the variables, their symbols and units, included in the model, and in each equation listed in Table 1 . This gs, along with internal carbon dioxide (Ci) and atmospheric carbon dioxide (Ca), determine A (Equation (7)), in accordance to the classic model of Farquhar (Farquhar and Sharkey, 1982; Farquhar et al., 1980) . A affected by LA, results in whole plant photosynthesis (Aplant, Equation (8)), that is then stoichiometrically transformed in carbohydrates and plant biomass (Equation (9)), assuming 0.648 g of assimilated carbohydrates day -1 m -2 LA, and days with an average of 15 h of sunlight. This was calculated by considering that assimilation of 1 μmol of CO2 s -1 m -2 is transformed into glucose by following 6CO2 þ6H2O ¼ glucose þ O2, and with a glucose molecular weight of 180.156 g mol-1. Respiration (Equation (10)) is deducted from Aplant prior to converting CO2 into Figure 1 . This model conceives vigor, defined as canopy dry weight (DWL þ DWS), as the interaction of environment and plant physiology. Root hydraulic conductance (Lp, 1) is calculated from conductance per unit of root dry weight (Lpr) and root biomass (12). The Lp will then constitute a proxy of the plant hydraulic conductance (kH; under no cavitation events, 2). This kH, when divided by leaf area (LA, 15), will constitute the leaf specific hydraulic conductance (kL, 3). The atmospheric mechanics are formalized by the Campbell and Norman (2012) equations of atmospheric deficit (q 0 , 4) and guard cell osmotic water potential (Ψπg, 5), that are then included in the Buckley et al. (2003) equation for stomatal conductance (gs, 6). Solar radiation, q 0 and temperature (TL) will affect this gs. Later, gs, along with internal carbon (Ci) and atmospheric carbon (Ca), will determine net photosynthesis (A, 7). A affected by LA, will result in whole plant A (Aplant, 8), that will stoichiometrically be transformed in carbohydrates and plant biomass (9), after respiration (10) is deducted from Aplant. Specific leaf area, (SLA, 16) theoretically depends on leaf turgor (ΨpL, 17), but for our purpose, SLA is an input value dependent on the genotype that affects LA. Finally, root biomass and canopy biomass (12 and 13) will be obtained affecting the plant biomass by a factor j dependent on GAs concentration in the tissues (14). GAs also affect growth by modifying biomass with a growth factor i (11). The model has two feedback loops. At time t -1 , LA will be defined by canopy biomass and SLA. Later, at time t, LA will define Aplant that will again define biomass and canopy biomass, restarting the cycle. The same happens with Lp, that at time t -1 , defines kH, that later defines gs, A and biomass. This biomass will re-define Lp at time t. carbohydrates. Biomass also depends on the concentration of bioactive GAs, as this type of hormones promotes growth by controlling cell division and expansion, as mentioned previously. We included a parametrized cell growth factor i, in the biomass equation (Equation (9)), that affects GAs and depends on Lpr (Equation (11)). Root biomass and canopy biomass (Equation (12) and Equation (13)) are obtained by affecting biomass with a partitioning factor, that results from a constant index j, also parameterized by us, that scales with the concentration of bioactive GAs in the tissues (Equation (14)). Finally, LA (Equation (15)) is obtained by multiplying canopy biomass with specific leaf area (SLA, Equation (16)), which, in theory, depends on leaf turgor (ΨPL, Equation (17)) and transpiration (E, Equation (18)) but constitutes an input in our model, and relies on the genotype. Canopy biomass ¼ Biomass Á partitioning factor (13) LA ¼ canopy biomass Á SLA This dynamic model has two feedback loops. At time t -1 , LA (Equation (15)) is defined by canopy biomass (Equation (13)) and SLA (Equation (16)). Later, at time t, LA will define Aplant (Equation (8)) that, in turn, will re-define biomass (Equation (9)). The same happens with Lp, that initially defines kH, which later will define gs, A and biomass. This biomass will re-define Lp in time t. Lpr and SLA constitute the inputs of the model, while canopy dry biomass (or vigor) is the output. Other variables, parameters and constants used in the model and not described in this text, are listed in Table 1 . The model was constructed using Vensim software PLE for Windows 5.4d (Ventana systems Inc, Harvard, MA 01451). 3. Model parametrization and validation One hundred and twenty-five genotypes from the population from a Ramsey  Riparia cross were multiplied through green cuttings (Hugalde et al., 2017) in spring 2014. Each plant was grown in a 4L capacity pot, trained to a single shoot and watered daily. We studied three replicates per genotype in a greenhouse with an average temperature of 25 C. Shoot growth rate was measured for 60 days, while, leaf area, leaf-shoot-root dry biomass, plant hydraulic conductance, stomatal conductance, and water potential were measured at the end of the experiment. Specific leaf area and root specific hydraulic conductivity were calculated at the end of the experiment. Raw data was grouped in clusters, as recommended for the analysis of large, quantitative, continuous and highly scattered data (Goshtasby, 2000; Berkhin, 2006) , following the methodology of Armas (1988, Table 2 ). In 2015, the same experiment was repeated and the variables were measured to parametrize the model. In addition, specific root hydraulic conductance was measured in a subset of 50 genotypes that differed in growth and biomass. Each variable was assessed as described in Hugalde et al. (2017) . For model validation, woody-cuttings from 30 contrasting genotypes were selected from the set used in 2015 for parametrization. They were bench-grafted with Cabernet Sauvignon and grown in the greenhouse as described before. After losing two, 28 genotypes were effectively used. They were divided in 9 groups of increasing vigor (Table 3) . For each group, the average value of the input parameters was used to run the simulations. Constants, parameters and units adopted for model fit and functioning are listed in Table 4 . To obtain the predicted values, the model was run with a timeframe setting of 90 days with initial conditions defined as 0.5 g of initial biomass for all genotypes, no water stress (Ψsoil ¼ -0.01 MPa, q' ¼ 0.03), and a temperature of 25 C. -Growth rate (b, m day-1): shoot length was measured every 3 days, during the first 6 weeks. Parameters (intercept and slopes) were obtained from piece-wise lineal equations. -Stomatal conductance (gs, mmol H2O m-2 s-1): was measured in selected genotypes (40 genotypes in 2014 and 50 in 2015) using a Decagon SC-1 leaf porometer (Decagon Devices Inc., USA). -Leaf and predawn water potential (ΨL and ΨPD, MPa): ΨL and ΨPD were measured at solar midday and predawn with a pressure chamber (PMS Instrument Company, OR, USA), according to Hsiao, 1990) . Predawn water potential (ΨPD) was considered as a proxy for soil water potential (Ψsoil). -Plant hydraulic conductance (kH, mmol m-2 s-1 MPa-1): was calculated through the Van den Honert equation, with Ψsoil, ΨL, and E. -Root hydraulic conductance (Lpr, m3 g DW-1 s-1 MPa-1): was measured after Barrios-Masias et al. (2015) , modified, with custom-made pressure chambers (PMS Instrument Company, OR, USA). Shoots were severed and the entire root systems were washed. Immediately, the root system was placed in plastic containers with deionized water inside the chamber and 1-3 cm stem was left protruding from the chamber lid. Once the system was sealed and secured, pressure was raised to 0.07 MPa and let to stabilize for 10 min. After stabilization, water exudation was measured gravimetrically by weighting a piece of dry gauze placed on top of the protruding stem for another 10 min (first cycle). Pressure flow assessments were conducted at increasing pressures of 0.1, 0.15, 0.2 and 0.25 MPa. At each step the system was stabilized for 7, 5, 5 and 3 min, respectively. Once the 5 cycles were completed, roots were removed from the chamber, pad dried and weighted. Lpr was calculated as the slope of the linear regression (R2!0.98) obtained from the pressure-flow relationship (units: mmol s-1MPa-1) and the dry weight of the entire root system (g DW). -Biomass (DW, g): after 60 days, leaves, shoots and roots, were dried separately at 60 C for 72 h, to obtain their dry weights (DWL, DWS, DWR, respectively). -Leaf area (LA, m2): leaves were scanned with a LI-3100C Area Meter (LiCOR Inc.). -Specific Leaf Area (SLA, m 2 g -1 ): was calculated by dividing LA with DWL. In addition, the pruning weights of two year-old plants growing in the UC Davis vineyard were measured in 108 genotypes of the Ramsey  Riparia population. After studying the variables one by one (Figures 2 and 3) , a multivariate approach was followed in order to analyze the variables and their relations. Principal Component Analyses (PCA) were conducted in 2014 and 2015 using Stat Graphics Plus 4.0 (Statistical Graphics Corp., StatSoft Inc., 2003). The PCA was able to explain 80% of the variability by its two first components. Component 1, accounting for 52.5% of total variability, showed strong positive effects of LA, growth rate (b), DWL and DWR, while strong and negative effect of Lpr ( Figure 4 ). This negative effect shows that in this progeny, more vigor corresponds to lower Lpr, meaning that smaller plants and smaller root systems are more effective in water absorption per biomass weight than vigorous plants. However, this is compensated by DWR, making vigorous root systems, with higher dry biomass, hydraulically superior for water delivery to the canopy. For component 2, accounting for 27.5% of total variability, positive effects were explained by SLA, and the partitioning index between LA and total biomass. The identification of the phenotyped traits that significantly correlated to vigor led the parametrization strategy through Lpr and SLA, the two inputs of the model, both clearly identified in the PCA. Using these input variables and adjusting the purported concentration of bioactive GAs, the model predicted the resulting vigor for each case, allowing the calculation of growth and partitioning indices i and j. The additional measured variables, e.g. E, DWR, DWL, used in the statistical analysis and PCA construction, served as model simulation controls. For parametrization of cell growth factor i and index j, the average value of input variables Lpr and SLA and output variable vigor, were calculated for each group of genotypes (Table 2 ). Significant correlation with L pr was found after analyzing the association between the obtained factor i and the variables in the model. In this way, Equation (19) (R Pearson ¼ 0.73; p-value ¼ 0.016), explained the relation between factor i and L pr . Consequently, cell growth factor i depends on L pr and is required for the calculation of biomass. cell growth factor i ¼ À 2:04 Á L pr þ 0:926 On the other hand, index j was inferred after the analysis of the relation between SLA and the residues calculated from the observed vs. predicted canopy biomass showed no fit (data not shown). It resulted in a constant equal to 0.05, which scaled partitioning with bioactive GA concentration for each time t (days of growing season). An empirical fit of the purported concentrations of bioactive GAs was calculated, based on the seasonal dynamics of some species (Kopcewicz, Table 3 . 1968; Y€ orük et al., 2004) and the measured values of leaves and buds of grapevines (Acheampong et al., 2015) . The equation fit is: The ability of the model to accurately estimate vigor was tested through the correlation between observed and predicted data, and by calculation of the mean-squared prediction error (MSPE), which accounts for the sum of three components, namely mean bias, line bias and random variation (Bibby and Toutenburg, 1977; Dhanoa et al., 1999) . The MSPE is the expected value of the squared difference between the fitted values implied by the predictive function and the values of the (unobservable) function. The root mean square error (RMSE), represents an analogue to standard deviation, and is calculated as square root of MSPE. Statistically significant correlations (R Pearson ¼ 0.85) were found between observed vigor and predicted values ( Figure 5 ). In our experiment, MSPE was 5 g, and the RMSE 2.2 g, being 0.65 g the lowest value measured in the progeny and 13.47 g the highest ( Table 2) . The mean bias, or difference between the observed mean and predicted mean, expressed as proportion of the MSPE (mean bias/MSPE), was 0.89. This means that almost 90% of the MSPE had mean bias. A mean bias different from zero indicates that predicted values can be consistently higher or lower than the observed values. Our estimated values were consistently lower than the observed ones, indicating that the model could strongly predict which plant was bigger or smaller, by sub-estimating the final biomass. The line bias, which is the deviation of the slope of the regression line from the line of unity, scaled to MSPE (line bias/MSPE), was 0.003. Finally, the random variation or the variability that is contained within a process that cannot be determined, scaled to MSPE (random variation/MSPE), was 0.10. The validation of the model allowed simulations of real vs. single variable modification scenarios that helped determine the individual effect of certain variables with respect to others. Since our model has two main plant inputs that control vigor, namely SLA and Lpr, plus the effect of hormone dynamics represented by GAs, we run it keeping two input variables constant while varying one. This made possible simulations of hypothetical situations that resulted in the quantification of the individual effect of the third variable as if "it was the only variable defining vigor differences". Simulations focusing on the individual effect of Lpr, GAs or SLA, using values matching a vigorous plant obtained in the parametrization step, showed that a big plant (high SLA, high GAs and 0.005 mmol H2O m-1 s-1 MPa-1 g-1) would be 42% bigger if it had the Lpr measured in a small plant (0.6 mmol H2O m-1 s-1 MPa-1 g-1) (Figure 6a ). On the other side, the individual effects of purported bioactive GAs (Figure 6b ) and measured SLA (Figure 6c) , using values corresponding to a big plant or a small plant, while keeping the other two variables constant (high SLA and 0.005 mmol H2O m-1 s-1 MPa-1 g-1 or high GAs and 0.005 mmol H2O m-1 s-1 MPa-1 g-1) were much stronger than that of Lpr. In other words, SLA and GAs play in favor of vigor in big plants, masking their less efficient water conductances. The 'compensation effect' is clearly exposed using real numbers of the three variables in the simulation (Figure 6d ). In this research, we studied vigor, growth and the mechanisms involved in biomass accumulation and partitioning. We considered growth a prime topic in relation to agricultural production, the inspiration for this study. The result was the generation of a mechanistic model that integrated both chemical and biophysical mechanisms in relation to growth. The input variables Lpr, SLA and purported amount of bioactive GAs were among the most important variables in defining vegetative vigor. This was proved mechanistically and is consistent with several studies demonstrating that hydraulics are related to growth and vigor (Clearwater et al., 2004; Marguerit et al., 2012; Sivasakthi et al., 2017) , and that active GAs are responsible for vegetative growth (Lambers et al., 1995 Moreno et al., 2011; Murcia et al., 2016) . Our model could also test SLA, and explain its key role in vigor, something previously observed and measured empirically. SLA integrates both most significant vigor measurements: leaf area (LA) and leaf biomass (DWL). This variable has large physiological significance: the higher its value, the bigger the LA that will capture light for photosynthesis per unit of biomass previously synthesized. It has been shown by others that vigorous plants have higher SLA values than small plants (Reich et al., 1998; Bultynck and Lambers, 2004; Nouvellon et al., 2010; Karavin, 2014) , as we confirmed with PCA analysis and our model. The validation of our model makes it a good foundation for future theoretical and applied research. When a theoretical model is validated, it becomes not only an accepted hypothesis, but also an important prediction tool, providing a quicker, confident and low-cost experimental foundation for studying specific phenomena (Fourcard et al., 2008; Ingalls, 2013 ). Our validated model can help study plant growth and development in an easier and less expensive way than the traditional empirical experiment-based methodology. The fact that the model includes key environmental variables that modulate plant physiology, permits to imagine virtual experiments assessing, for instance, plants response to vineyard management practices (i.e. irrigation, shoot thinning, leaf removal, etc.), or even climate change effects on plant growth (i.e., temperature, humidity, water availability, radiation, etc.). The model correctly estimates relative plant size, but the fact that it consistently sub-estimates vigor, should be studied and improved by a fine scaling parametrization. Another characteristic to be observed is that the model does not consider the possible effects of interactions within the canopy because it was developed in plants trained to just one stem. Also, the general nature of some of the equations does not consider aspects like plant nutrition in photosynthesis, or embolism after water scarcity in hydraulics. A future coupling with a previous cavitation model (Hugalde et al., 2018) could cover this issue. The use of Lp as a proxy for kH is another limitation that could be addressed to refine the model, as well as more localized measurements like fine rooting. The lack of hormone measurements was partially overcome by the inclusion of an empirical equation fitting seasonal patterns of GAs previously studied in woody species and grapevines (Kopcewicz, 1968; Y€ orük et al., 2004; Acheampong et al., 2015) . This methodology is widely accepted in modelling, since many parameters that are difficult to assess in situ, are often adopted from bibliography, especially in biological/molecular subjects (Magni and Saparcino, 2014; Cobelli and Carson, 2019) . Peccoux et al. (2018) , for example, adopted parameters from literature when including the Penman-Monteith equation to a model for rootstock control of scion transpiration in grapevine. Measurement of GAs is challenging, costly and difficult, especially because accurate measurement of specific bioactive GAs requires sophisticated equipment and isotopically labeled internal standards ( Simura et al., 2018) . An alternative to validate GAs measurements is using specific inhibitors of the synthesis of bioactive GAs, type Ca-Prohexadione. However, application of Ca-Prohexadione is phytotoxic on grapevines (R. Bottini, personal communication). Nevertheless, the patterns of seasonal changes of GAs in different species including grapevines, is quite similar (Hedden and Sponsel, 2015) , providing confidence in the equation used in this model. Finally, the inclusion of empirical measurements of the model's parameters might need adjustments for other species, such as very small or very vigorous plants, trees or herbs. Model simulations provided clarity on the individual effect of each variable that accounted for large portions of the observed variability. It also allowed us to understand the effect that some variables play over others by comparing the model situations with reality. Simulations with the model using variables from the progeny, found that SLA and GA concentration play vital roles in the development of bigger canopies and root systems with higher whole root hydraulics. In our progeny, the measured Lpr was lower in large plants. This may be understood as an adaptive strategy to improve their root hydraulic conductivity, e.g. plants with small root systems develop higher specific conductivities and plants with low specific conductivity grow larger root systems. Responses of vigor to hydraulic conductivity and its relationship with carbohydrate partitioning have been reported in other systems. For example, genetic transformation of Arabidopsis thaliana with an anti-sense construct targeted to the PIP1b aquaporin gene, resulted in reduced cellular water permeability that was accompanied with an increase of root mass to ensure sufficient water supply to the plant (Kaldenhoff et al., 1998) . Overexpression of the same gene in tobacco improved plant vigor, transpiration rate, stomatal density, and photosynthetic efficiency under favorable conditions (Aharon et al., 2003) . Also, studies in grafted kiwi plants showed that whole plant hydraulic conductance was lower but leaf-area-specific conductance (kL) and gs were both higher in low-vigour rootstocks (Clearwater et al., 2004) . These authors suggest that changes in biomass partitioning to the roots and crown structure are involved in the observed rootstock effects. Another study with two chickpea genotypes contrasting for vigor showed that early vigor genotypes had lower root hydraulic conductivity and transpiration rates than late vigor lines (Sivasakthi et al., 2017) . However, this would be compensated by their larger root system under low vapor pressure deficit conditions. The fact that ungrafted rootstocks grow vigorously does not imply that they will confer more vigor to the scion. For example, relatively vigorous motherlines like 420A, 5C, and 101-14 perform as low vigor rootstocks, while St George, 110R, 140Ru and Ramsey are relatively weak mother vines that confer high vigor to the scions (M. Andrew Walker personal observation). The interaction between scion and rootstock has been widely studied, since the relation between both genotypes can produce varied phenotypes; and both canopy and roots must be involved in whole plant regulation (Tandonnet et al., 2010; Cookson and Ollat, 2013; Albacete et al., 2015; Miele and Rizzon, 2017) . In our study, Ramsey displayed less vigor than Riparia Gloire (Hugalde et al., 2017) . However, it is well documented that they are considered invigorating and devigorating rootstocks respectively. According to our model, this would imply that grafting modifies the relation between variables like Lpr, root biomass, etc. so that we could speculate that the higher Lpr displayed by Ramsey would become a main factor influencing the vigor of grafted plants. This is in accordance with a study in grapevines that showed that fine roots of invigorating and drought resistant 110R had higher Lpr than 101-14, a less invigorating rootstock (Barrios-Masias et al., 2015) . The same was observed before by Gambetta et al. (2012) , who found that invigorating rootstocks tended to have higher L pr than rootstocks that conferred low vigor to the scion. In addition, Peccoux et al. (2018) found that Cabernet Sauvignon grafted onto 110R exhibited higher gs and E than Cabernet Sauvignon grafted onto 101-14 under well-watered and moderate water deficit, probably due to the higher root length area developed by 110R. . But, shoot physiological responses of Cabernet Sauvignon grafted onto Ramsey and Riparia Gloire were similar under drought stress (Barrios-Masias et al., 2018) . It is noteworthy that in our research, model validation was achieved using 3 months old grafted plants, meaning that the model was able to predict vigor in a quite accurate manner in such a different set of plants. In the future, it might be interesting to include a parameter that incorporates grafting effects, in order to better comprehend the relationship between 2 different genotypes and their partitioning. This was not evident in our study since the plants were only 2-3 month old but could be manifested in older plants in the field. The model allowed the assembling of the relevant functional traits that underlay vigor and the comprehension of their interactions and emerging features. We believe that our model could be the initial step in developing a practical diagnosis tool for breeders and growers. It is known that balanced vigor is crucial for grapevine production. Modeling approaches could become an instrument for analysis of crop growth, assisting vineyard management. In addition, models like the one we propose may help ecological and environmental studies address the climate change scenario we are currently facing. The validation of the theoretical model developed for grapevine helped to identified the most important mechanisms underlying vegetative growth. It simulated growth in an acceptable manner, empowering its potential use as a prediction tool that could be applied in the fields of viticulture and ecology. Root hydraulics, theoretical hormone dynamics and specific leaf area were the key variables in our model that controlled vegetative growth. It remains necessary to measure hormones to further refine the model. However, we propose that by measuring Lpr and SLA in young plants, and estimating hormone impacts, our model can predict, with acceptable accuracy, the expected vigor of different genotypes. In es P. Hugalde This work was supported by the Instituto Nacional de Tecnología Agropecuaria (INTA) and the California Grape Rootstock Improvement Commission. Data included in article/supplementary material/referenced in article. The authors declare no conflict of interest. I.P. Hugalde et al. Heliyon 6 (2020) e05708 No additional information is available for this paper. Functional characterization and developmental expression profiling of gibberellin signaling components in Vitis vinifera Overexpression of a plasma membrane aquaporin in transgenic tobacco improves plant vigor under favorable growth conditions but not under drought or salt stress Unravelling rootstock  scion interactions to improve food security Estadística Sencilla: Descriptiva Root and stem hydraulic conductivity as determinants of growth potential in grafted trees of apple (Malus pumila Mill.) The influence of photosynthetic acclimation to rising CO2 and warmer temperatures on leaf and canopy photosynthesis models Differential responses of grapevine rootstocks to water stress are associated with adjustments in fine root hydraulic physiology and suberization Differences in hydraulic traits of grapevine rootstocks are not conferred to a common Vitis vinifera scion Gibberellins and heterosis in poplar A survey of clustering data mining techniques Prediction and Improved Estimation in Linear Models Association of dwarfism and floral induction with a grape 'green revolution' mutation A hydromechanical and biochemical model of stomatal conductance The control of stomata by water balance Effects of applied gibberellic acid and paclobutrazol on leaf expansion and biomass allocation in two Aegilops species with contrasting leaf elongation rates An Introduction to Environmental Biophysics Azospirillum brasilense and Azospirillum lipoferum hydrolyze conjugates of GA20 and metabolize the resultant aglycones to GA1 in seedlings of rice dwarf mutants Mutants at the Slender1 locus of barley cv Himalaya. Molecular and physiological characterization Understanding plant responses to drought -from genes to the whole plant Hydraulic conductance and rootstock effects in grafted vines of kiwifruit Growth habit of Lotus tenuis shoots and the influence of photosynthetic photon flux density, sucrose and endogenous levels of gibberellins A1 and A3 Introduction to Modelling in Physiology and Medicine Grafting with rootstocks induces extensive transcriptional re-programming in the shoot apical meristem of grapevine Growth of the plant cell wall Biosynthesis of hormones and elicitor molecules Effects of rootstock and irrigation regime on hydraulic architecture of Vitis vinifera L. cv Use of mean square prediction error analysis and reproducibility measures to study near infrared calibration equation performance Influence of different rootstocks on the vegetative and reproductive performance of Vitis vinifera L. Malbec under irrigated conditions Aquaporin-mediated reduction in maize root hydraulic conductivity impacts cell turgor and leaf elongation even without changing transpiration Hydraulic conductivity traits predict growth rates and adult stature of 40 Asian tropical tree species better than wood density Stomatal conductance and photosynthesis A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species Plant growth modelling and applications: the increasing importance of plant architecture in growth models The relationship between root hydraulics and scion vigour across Vitis rootstocks: what role do root aquaporins play? Grouping and parameterizing irregularly spaced points for curve fitting A century of gibberellin research Initial alteration of scion architecture by dwarfing apple rootstocks may involve shoot-root-shoot signaling by auxin, gibberellin, and cytokinin Measurements of plant water status Gene expression analysis in citrus reveals the role of gibberellins on photosynthesis and stress The phenomenon of cavitation in grapevine… Unravelling implicated mechanisms Physiological and genetic control of vigour in a 'Ramsey'  'Riparia Gloire de Montpellier' population Slender rice, a constitutive gibberellin response mutant, is caused by a null mutation of the SLR1 gene, an ortholog of the heightregulating gene GAI/RGA/RHT/D8 Mathematical Modelling in Systems Biology: an Introduction Principal component analysis: a review and recent developments Significance of plasmalemma aquaporins for water-transport in Arabidopsis thaliana Effects of leaf and plant age on specific leaf area in deciduous tree species Quercus cerris L. var The metabolism of gibberellin A20 to gibberellin A1 by tall and dwarf mutants of Oryza sativa and Arabidopsis thaliana Seasonal changes of gibberellin-like substances and growth inhibitors in the apical meristems of Pine (Pinus sylvestris L.) Growth and allocation The control of biomass partitioning in plants from ''favourable'' and ''stressful'' environments: a role for gibberellins and cytokinins Local PCA shows how the effect of population structure differs along the genome Gibberellin retards chlorophyll degradation during senescence of Paris polyphylla Drought-induced changes in development and function of grapevine (Vitis spp.) organs and in their hydraulic and non-hydraulic interactions at the whole-plant level: a physiological and molecular update Mercurial inhibition of root hydraulic conductance in Vitis spp. rootstocks under water stress Genetic linkage map of the interspecific grape rootstock cross Ramsey (Vitis champinii)  Riparia Gloire (Vitis riparia) Parameter estimation Rootstock control of scion transpiration and its acclimation to water deficit are controlled by different genes Rootstock-scion interaction: 2. Effect on the composition of Cabernet Sauvignon grape must Gibberellins and abscisic acid promote carbon allocation in roots and berries of grapevines Characterization of a dwarf gene in Brassica rapa, including the identification of a candidate gene ABA and GA3 increase carbon allocation in different organs of grapevine plants by inducing accumulation of non-structural carbohydrates in leaves, enhancement of phloem area and expression of sugar transporters Within-stand and seasonal variations of specific leaf area in a clonal Eucalyptus plantation in the Republic of Vigour conferred by rootstock: hypotheses and direction for research Comparison of three approaches to model grapevine organogenesis in conditions of fluctuating temperature, solar radiation and soil water content Dissecting the rootstock control of scion transpiration using model-assisted analyses in grapevine Hydrolysis of [17,17-2H2] gibberellin A20-glucoside and [17,17-2H2]gibberellin A20 -glucosyl ester by Azospirillum lipoferum cultured in a nitrogen free biotin-based chemically-defined medium Rapid chlorophyll a fluorescence light response curves mechanistically inform photosynthesis modeling Coupling of cell division to cell growth by translational control of the G1 cyclin CLN3 in yeast A putative role for TIP and PIP aquaporins in dynamics of leaf hydraulic and stomatal conductances in grapevine under water stress and re-watering Simple scaling of photosynthesis from leaves to canopies without the errors of big-leaf models Phenotypic and genetic dissection of component traits for early vigour in rice using plant growth modelling, sugar content analyses and association mapping Leaf structure (specific leaf area) modulates photosynthesis-nitrogen relations: evidence from within and across species and functional groups Stomatal response of an anisohydric grapevine cultivar to evaporative demand, available soil moisture and abscisic acid Gibberellins: a phytohormonal basis for heterosis in maize Fisiología de las plantas 1. C elulas: agua, soluciones y superficies Modelling approach for predicting the impact of changing temperature conditions on grapevine canopy architectures Chickpea genotypes contrasting for vigor and canopy conductance also differ in their dependence on different water transport pathways Hydraulic conductivity of roots, mineral nutrition, and leaf gas exchange of citrus rootstocks Plant Physiology and Development Central roles for potassium and sucrose in guard-cell osmoregulation The role of sucrose in guard cell osmoregulation Scion genotype controls biomass allocation and root development in grafted grapevine Variability among species of stomatal control under fluctuating soil water status and evaporative demand: modelling isohydric and anisohydric behaviours Gibberellin signaling in the endodermis controls Arabidopsis root meristem size The role of plasma membrane intrinsic protein aquaporins in water transport through roots: diurnal and drought stress responses reveal different strategies between isohydric and anisohydric cultivars of grapevine Stomatal conductance correlates with photosynthetic capacity Connecting biochemical photosynthesis models with crop models to support crop improvement Gibberellin metabolism and its regulation C 3 and C 4 photosynthesis models: an overview from the perspective of crop modelling Seasonal changes of endogenous plant hormones in Rosa canina A 3-D functional-structural grapevine model that couples the dynamics of water transport with leaf gas exchange Modelling grape growth in relation to whole-plant carbon and water fluxes Anatomical characteristics of cherry rootstocks as possible preselecting tools for prediction of tree vigour Stomatal behaviour of irrigated Vitis vinifera cv. Syrah following partial root removal We thank Karla Huerta, Cassandra Bent and Jake Uretsky for their valuable help in measurements. Special acknowledgements are for Dr. Hebe Cremades, Dr. Rub en Bottini, Dr. Luis Mayorga and Dr. Kent Bradford, who helped discussing key parts of this research paper.This publication documents part of In es P. Hugalde's PhD program at the Universidad Nacional de Cuyo, Mendoza, Argentina and the University of California, Davis.We respectfully remember Dr. Hern an F. Vila, who passed away after fighting COVID 19 in October, 2020, and dedicate this article to his memory.