Diverse Physiological and Physical Responses among Wild, Landrace and Elite Barley Varieties Point to Novel Breeding Opportunities

: Climate change from elevated [CO 2 ] may reduce water availability to crops through changes in precipitation and higher temperatures. However, agriculture already accounts for 70% of human consumption of water. Stomata, pores in the leaf surface, mediate exchange of water and CO 2 for the plant. In crops including barley, the speed of stomatal response to changing environmental conditions is as important as maximal responses and can thus affect water use efﬁciency. Wild barleys and landraces which predate modern elite lines offer the breeder the potential to ﬁnd unexploited genetic diversity. This study aimed to characterize natural variation in stomatal anatomy and leaf physiology and to link these variations to yield. Wild, landrace and elite barleys were grown in a polytunnel and a controlled environment chamber. Physiological responses to changing environments were measured, along with stomatal anatomy and yield. The elite barley lines did not have the fastest or largest physiological responses to light nor always the highest yields. There was variation in stomatal anatomy, but no link between stomatal size and density. The evidence suggests that high photosynthetic capacity does not translate into yield, and that landraces and wild barleys have unexploited physiological responses that should interest breeders.


Introduction
Rising concentrations of greenhouse gases including CO 2 ([CO 2 ]) are expected to raise global mean surface temperatures by 3 • C above the baseline by the end of this century [1,2]. Extreme weather is also becoming more common, with heat waves, droughts and precipitation events expected more often, that are more severe and for longer periods than in the past [1]. Meanwhile, the global population is expected to rise from 7.2 to 9.6 billion by 2050 [3] putting increasing pressure on food production. It has been predicted that crop yields have to double by 2050 as a result [4]. Widespread use of biofuels and changing dietary preferences along with rising urbanisation [5,6] are also adding pressure to raise yields in the key crops including rice, wheat and maize that are responsible for the vast majority of all calories consumed globally [7].
As a result of climate-induced surface temperature increases, photosynthesis could be reduced in some temperature-sensitive plants (dependent on species and variety), decreasing productivity [8]. Direct thermal damage to plants will rise, while soil water availability will fall [8]. Greater evapotranspiration (the process by which liquid water enters the gas phase either from the soil (evaporation) or plants (transpiration)) may further reduce soil moisture content via higher stomatal conductance (g s ) (and thus transpiration rates) in some C3 plants which demand greater evaporative cooling, although elevated [CO 2 ] tends to reduce stomatal aperture [9][10][11]. High vapour pressure deficit (VPD) acts in addition to direct temperature effects on evapotranspiration by reducing g s [10,12].
Meanwhile changed patterns of precipitation will additionally reduce replenishment rates of water stores [1, 3,13,14]. As agriculture already accounts for up to 70% of freshwater withdrawals worldwide, the combined threat of climate change and yield pressure on water resources is clear, with inevitable consequences expected on crop yields [1, 3,13,15,16].
One of the key determinants of the rate of g s and transpiration is the number, size and position of stomatal pores on the leaf external surfaces [17]. Stomata must simultaneously manage mesophyll demands for CO 2 diffusion into the leaf for assimilation (A), while minimizing water loss through transpiration [18]. The ratio of water loss through stomatal conductance to assimilation rate is known as intrinsic water use efficiency (WUE i = A/g s ) [19]. Reducing stomatal density and thereby g s has been shown to enhance crop biomass with lower water use in rice and barley [20,21], while greater WUE i can also be achieved by increasing assimilation rates without using more water, through genetic manipulation of photosynthetic pathways [22][23][24][25]. Steady-state gas exchange measurements have shown a close proportional relationship between A and g s [26] and several studies have demonstrated that g s is also correlated with final yield in the field [27,28]. However under ambient conditions, such as those observed in the field, 'steady-state' is rarely achieved; light intensity can change rapidly (seconds to minutes) over a large range of intensities [29]. The assimilation rate adjusts within seconds in response to changing light intensity or other environmental cues, while g s is often an order of magnitude slower to respond [30][31][32][33]. Therefore, a disconnect between changes in A and g s can lead to periods of reduced CO 2 assimilation relating to diffusional constraints on A from low g s , as well as periods of unnecessary water loss when the slow speed of stomatal closure does not match a rapid decline in A [30,31]. The kinetic responses of stomata including the rapidity and magnitude of change in g s provide a potentially unexploited target for increasing crop yield [30,34]. Small stomata are reported to be faster to respond to changing environmental cues owing to greater surface area to volume ratio [17], enhanced further in dumbbell-shaped guard cells, which are surrounded by subsidiary cells, the latter facilitating speedy stomatal movements by providing a reservoir of solutes that can be rapidly transferred to guard cells where they act as an osmotic to alter turgidity [29,35,36]. Therefore, before we can exploit these potential novel breeding targets we must understand the underlying mechanisms responsible for g s dynamics, and the processes and pathways that coordinate A and g s [37].
Historically, breeding has been directed toward yield improvement with a consequent loss of genetic diversity [4]. Heterogenous stomatal kinetic responses do however exist within and between C3 crop species and cultivars [38][39][40][41][42]. Exploiting such variation in landraces, wild relatives and elite varieties of major crops, including barley, could provide an avenue to recapture genetic diversity and identify new traits for exploitation [31,43,44].
The aim of this study was to characterise natural variation in stomatal anatomy and leaf physiology of wild barleys and landraces relative to elite varieties grown in both field and controlled environments to understand the extent of variation in g s , and A in non-elite cultivars compared to commercial varieties and assess the impact of stomatal responses on these. By linking physiological to anatomical data we hope to provide a greater understanding of underlying mechanisms that coordinate between the two, and how these may be exploited for improved yield potential. This study showed that there was no link between stomatal size and density and that the non-elite barley lines possessed physiological traits that would be attractive to breeders.

Plant Material and Growth Conditions
A wide combination of varieties of barley Hordeum vulgare L and H. vulgare ssp. Spontaneum (Table 1) were sourced from KWS (KWS UK Ltd., 56 Church St, Thriplow, UK) and John Innes Centre Germplasm Resources Unit (JIC GRU, Norwich Research Park, Norwich, UK) encompassing two and six-row ears, spring or winter habits and malt/feed end-uses. Table 1. List of elite, wild and landrace barleys, their key selection criteria, the experiments in which they were used and the source of the seed material. Selection criteria were the varietal type (wild, landrace or elite cultivars), the number of rows on the ear (2 or 6), the habit or vernalization requirement (winter or spring) and the typical end use of the variety (feed or malt). A Polytunnel and growth chamber were used to permit a degree of control over the environment in which the experiments were conducted. Seeds were sown on 24 March 2016 into modular pots in Levington F2 + S at the KWS UK Ltd. field site at Thriplow, Cambs. Following an initial 5 days stratification in the dark at 5 • C, the seeds were transferred to a glasshouse. After three weeks they were transplanted into a polytunnel for the duration of the experiment in order to assure rapid establishment and growth, in a randomised 6-block design with one biological replicate per block. All six varieties (Table 1) were grown in each block, with double guard rows surrounding them, between March and July 2016 and were watered weekly. No additional fertiliser was used. Three groups of varieties were selected: 2 wild barleys originating in Central Asia (B3733 and B3745); 2 Central/Southern European landraces (Golden Archer and Hatif de Grignon) and 2 elite cultivars (KWS Orwell and KWS Sassy).

Variety
In addition, a selection of landrace barleys (Dea, Eire 6-Row, Golden Archer and Hatif de Grignon) as well as the elite variety KWS Irina were germinated on damp paper towels in two separate experiments in 2016 in the glasshouse at the University of Essex. On germination, they were transplanted into Levington F2 + S in 1 L pots and transferred to a controlled growth chamber (Conviron, Isleham, Cambs, UK) at 23 • C day/15 • C night at ambient humidity and CO 2 with 188 +/− 6 µmol m −2 s −1 PPFD.

Infra-Red Gas Exchange Measurements of Photosynthesis.
Assimilation rate (A) was measured as a function of light intensity to assess photosynthetic performance. The mid-portion of the youngest fully expanded leaf (rather than a particular leaf succession, to account for developmental differences in the timing of leaf emergence) of each replicate was placed in a LiCOR 6400 (Li-COR Inc, Lincoln, NE, USA) chamber at a light intensity of 1500 µmol m −2 s −1 PPFD and left until A was stable (15-30 min). PPFD was then decreased in a stepwise manner and A recorded when a new steady state had been achieved (ca. 2 min). The light intensities used were 1500, 1250, 1100, 900, 700, 600, 500, 400, 200, 100, 50 and 0 µmol m −2 s −1 PPFD. VPD was maintained at 1 KPa, temperature at 22 • C [45] and [CO 2 ] at 400 µmol mol −1 . Data were modelled as Michaelis-Menten [46] enzyme kinetic parameters with PPFD as the substrate with concentration [S].
The rate of assimilation A was thus a function of R d , the respiration rate in the dark, A sat the maximal assimilation rate at saturating light intensity and K m the Michaelis-Menten constant where the reaction rate was half maximum. The initial slope (i.e., quantum Assimilation rate was also measured as a function of [CO 2 ] to assess the impact of internal CO 2 concentration on photosynthetic capacity. The mid-portion of the youngest fully expanded leaf of each replicate was placed in a LiCOR 6400 (Li-COR Inc, Lincoln, NE, USA) chamber at a light intensity of 1500 µmol m −2 s −1 PPFD and left until A was stable (15-30 min). Partial pressure of CO 2 was then decreased in a stepwise manner and A recorded when a new steady state had been achieved (ca. 2 min). The CO 2 concentrations used were 400, 250, 150, 100, 50, 400, 550, 700, 900, 1100 and 1300 µmol CO 2 mol −1 air. VPD was maintained at 1 KPa and temperature at 22 • C. Data were modelled using the methods outlined by Duursma et al. [47] by curve-fitting the Farquhar-Berry-von Caemmerer model of stomatal conductance and expressing the results as three parameters of a series of nonlinear regressions; V c max , the maximum velocity of Rubisco for carboxylation, J max , the maximum rate of electron transport demand for Ribulose 1,5 bisphosphate and R d , the rate of respiration in the dark.

Stomatal Kinetic Responses
Kinetic responses of A and g s were measured by placing the youngest fully expanded leaf in the cuvette of an infra-red gas analyser (LiCOR 6400, Li-COR Inc., Lincoln, NE, USA). Leaves were first equilibrated at a PPFD of 0 µmol m −2 s −1 until both A and g s were stable (ca. 20-30 min). PPFD was then increased to 1000 µmol m −2 s −1 for 40 min. The leaf cuvette was maintained at 400 µmol mol −1 CO 2 concentration (C a ), a block temperature of 25 • C and a VPD of 1. Assimilations and stomatal conductance were recorded every 1 min. Intrinsic water use efficiency (WUE i ) was calculated as WUE i = A/g s . Data were modelled according to the method in Vialet-Chabrand et al.  [39]. The model described the temporal response of g s at time t using a time constant (τ, min), an initial time lag (λ, min) and a steady-state g s (G s1000 , mmol m −2 s −1 ) reached at given PPFD: while the change in A was modelled as a simple exponential without the lag, λ: Time t = 0 was the point at which PPFD was increased from 0 to 1000 µmol m −2 s −1 ; g 0 (mmol m −2 s −1 ), was the initial value of stomatal conductance before the change in PPFD, A 0 (µmol m −2 s −1 ) the initial value of A before a change in PPFD and A 1000 (µmol m −2 s −1 ) the value of A at 1000 µmol m −2 s −1 PPFD.
A second parameter combining rapidity and amplitude of the response, the maximum slope (Sl max ), was used to describe the maximal slope of the g s and A responses to the step-change in PPFD:

Stomatal Anatomical Characteristics
Stomatal impressions of the abaxial leaf surfaces were taken at the same locations measured during gas exchange. A negative impression was made using dental polymer (Xantoprene, Heraesus Kulzer Ltd., Hanau, Germany) following the methods of Weyers and Johansen [49]. A positive impression was made from the dried polymer by painting with nail varnish and placing the dried film on a microscope slide. Three pseudo-replicates from each impression (with a field of view of~1250 µm 2 ) were averaged to give stomatal density and stomatal pore length determined using a calibrated 5 MP eyepiece camera (MicroCAM 5 MP, Bresser Optics, Rhede, Germany).

Growth Measurements
A number of measurements important in partitioning growth and development were taken at harvest 162 days after sowing when the plants were visibly mature and dry: tiller number, fertile ears tiller −1 , harvest Index (grain mass/above-ground biomass, HI [50], thousand grain weight (TGW), grain number plant −1 and grain mass plant −1 [51]. Thousand grain weight is a measure of yield quality that is known to have a genetic component [52].

Statistical Analyses
Statistical analyses were conducted using R and RStudio [53]. A two-way analysis of variance (ANOVA) was used for gas exchange data when two factors (genotype × block) were present (i.e., for the variables A, g s , τ, A 1000 , Gs 1000 and Sl max of A and g s for fieldwork) with outliers greater than 3 sd from the mean removed. Single factor analyses were carried out using one-way ANOVA (i.e., for A, g s , A sat , g s sat , V c max , J max in the controlled environment chamber). Curve fitting in R used the plantecophys package [47]. Shapiro-Wilk's and Breusch-Pagan's tests were used to test data for normality and homogeneity of variance, respectively. The strength of trait associations (between stomatal density and size, and with kinetic and maximal responses) were measured using Pearson's correlation coefficient.

Analysis of Steady-State and Kinetics of g s and A in Field-Grown Varieties
To evaluate the responsiveness of A and g s to changing light intensity in the field, we applied a step change in light intensity from low (0 µmol m −2 s −1 PPFD) to high (1000 µmol m −2 s −1 PPFD) and measured the responses using infra-red gas exchange analysis ( Figure 1). The data were modelled as per the exponential Equations (3)-(5) described in the methods section ( Figure 2). Assimilation was negative at the start of the measurements at 0 µmol m −2 s −1 PPFD ( Figure 1A), and initially rose rapidly after PPFD was increased to 1000 µmol m −2 s −1 before reaching a plateau toward the end of the measurement period. There was a lag in the response of g s to increasing light intensity, but it too rose, albeit more slowly, toward a plateau toward the end of the period. Maximal rates (taken as the last five minutes of the measurement period) for A and gs at 1000 μmol m −2 s −1 PPFD (Figure 2A,B) were significantly different across genotypes (F(5, 70) = 23.2 and F(5, 70) = 9.8, respectively, p < 0.05). The wild barley B3733 had the lowest Maximal rates (taken as the last five minutes of the measurement period) for A and gs at 1000 μmol m −2 s −1 PPFD (Figure 2A,B) were significantly different across genotypes (F(5, 70) = 23.2 and F(5, 70) = 9.8, respectively, p < 0.05). The wild barley B3733 had the lowest maximal assimilation (11.0 μmol m −2 s −1 ) while KWS Sassy had the highest (21.7 μmol m −2 s −1 ); similarly there was a 2-fold difference in maximal gs between the landrace Golden Archer and the elite cultivar KWS Sassy.  Maximal rates (taken as the last five minutes of the measurement period) for A and g s at 1000 µmol m −2 s −1 PPFD (Figure 2A,B) were significantly different across genotypes (F (5, 70) = 23.2 and F (5, 70) = 9.8, respectively, p < 0.05). The wild barley B3733 had the lowest maximal assimilation (11.0 µmol m −2 s −1 ) while KWS Sassy had the highest (21.7 µmol m −2 s −1 ); similarly there was a 2-fold difference in maximal g s between the landrace Golden Archer and the elite cultivar KWS Sassy.
We further characterised the speed of responses to changing light intensity from 0 µmol m −2 s −1 to 1000 µmol m −2 s −1 , separately from maximal values ( Figure 2C-F). There was no genotypic variation in the time constant of assimilation responses (p = 0.51, Figure 2C), although the overall model was significant for the maximum rate of change of the assimilation response (F (5, 70) = 3.45, p < 0.05, Figure 2E) and while A for KWS Orwell increased by 55.1 µmol m −2 s −1 min −1 compared to 19.1 µmol m −2 s −1 min −1 for B3733, Tukey's test was unable to distinguish between the varieties. In addition, there was a clear effect of variety on stomatal responses to changing light intensity. The time constant of g s response to increasing light was affected by variety (F (5, 71) = 2.41, p < 0.05, Figure 2B), although not the rate of response of g s (p = 0.48, Figure 2D). The time constants for g s ( Figure 2B) were longest for the landraces Golden Archer (6.3 min) and Hatif de Grignon (5.6 min), and shortest for the wild barley, B3733 (3.3 min). Meanwhile, there was trend for the wild barleys (B3733 = 0.39 µmol m −2 s −1 min −1 and B3745 = 0.5937 µmol m −2 s −1 min −1 ) to have the slowest rate of change in response of g s to increasing light intensity ( Figure 2F), while landraces and elite varieties were the fastest to react by 43% or more.

Field Trial Stomatal Anatomy
To explore whether the responses to changing light intensity resulted from stomatal anatomical variability, stomatal impressions were taken of the leaves subjected to the step change measurements (Figure 3A-C). There was a significant genotypic variation in stomatal density ( Figure 3A, F (5, 22) = 7.0, p < 0.05) but not size ( Figure 3B, p = 0.056) and no correlation between the two ( Figure 3C, p = 0.20), even when accounting for the effect of variety on that relationship (p = 0.52). The wild barley variety B3745 had the lowest stomatal density (41.9 mm −2 ) and along with the landrace Hatif de Grignon (42.6 mm −2 ) had density that was significantly lower than B3733 (61.8 mm −2 ), Golden Archer (66.2 mm −2 ) and KWS Sassy (59.4 mm −2 ). Elsewhere, KWS Sassy had a pore length that was only 13.5% longer than that of the variety with the shortest pores, B3733 ( Figure 3C). creased by 55.1 μmol m s min compared to 19.1 μmol m s min for B3733, Tukey's test was unable to distinguish between the varieties. In addition, there was a clear effect of variety on stomatal responses to changing light intensity. The time constant of gs response to increasing light was affected by variety (F(5, 71) = 2.41, p < 0.05, Figure 2B), although not the rate of response of gs (p = 0.48, Figure 2D). The time constants for gs ( Figure  2B) were longest for the landraces Golden Archer (6.3 min) and Hatif de Grignon (5.6 min), and shortest for the wild barley, B3733 (3.3 min). Meanwhile, there was trend for the wild barleys (B3733 = 0.39 μmol m −2 s −1 min −1 and B3745 = 0.5937 μmol m −2 s −1 min −1 ) to have the slowest rate of change in response of gs to increasing light intensity ( Figure 2F), while landraces and elite varieties were the fastest to react by 43% or more.

Field Trial Stomatal Anatomy
To explore whether the responses to changing light intensity resulted from stomatal anatomical variability, stomatal impressions were taken of the leaves subjected to the step change measurements (Figure 3A-C). There was a significant genotypic variation in stomatal density ( Figure 3A, F(5, 22) = 7.0, p < 0.05) but not size ( Figure 3B, p = 0.056) and no correlation between the two ( Figure 3C, p = 0.20), even when accounting for the effect of variety on that relationship (p = 0.52). The wild barley variety B3745 had the lowest stomatal density (41.9 mm −2 ) and along with the landrace Hatif de Grignon (42.6 mm −2 ) had density that was significantly lower than B3733 (61.8 mm −2 ), Golden Archer (66.2 mm −2 ) and KWS Sassy (59.4 mm −2 ). Elsewhere, KWS Sassy had a pore length that was only 13.5% longer than that of the variety with the shortest pores, B3733 ( Figure 3C). To characterize the relationship between stomatal anatomy and g s kinetic responses further, the relationships between stomatal density and A or g s at high light intensity were explored. There was no correlation between stomatal density and A 1000 (p = 0.64), or with gs 1000 (p = 0.57) even accounting for variety. Similarly, there was no evidence that stomatal size was linked to the kinetic responses to changing light intensity including any of τ A (p = 0.84), τ gs (p = 0.61), Sl max A (p = 0.86) or Sl max g s (p = 0.92) (Supplementary Materials).

Field Trial Harvest
At the completion of the field trial, a range of data relating to harvest were recorded ( Figure 4) to understand the variability in yield and yield components among the different varieties, and to attempt to link physiological and anatomical traits to yield. The elite variety KWS Sassy outperformed other varieties across all metrics compared to landrace and wild varieties ( Figure 4A-H). The landrace Golden Archer delivered yield and yield component outcomes that were comparable with those for KWS Orwell on all variables measured and contrasted strongly with the performance of the landrace Hatif de Grignon ( Figure 4A-H). gs1000 (p = 0.57) even accounting for variety. Similarly, there was no evidence that stomatal size was linked to the kinetic responses to changing light intensity including any of τA (p = 0.84), τgs (p = 0.61), Slmax A (p = 0.86) or Slmax gs (p = 0.92) (Supplementary Materials).

Field Trial Harvest
At the completion of the field trial, a range of data relating to harvest were recorded ( Figure 4) to understand the variability in yield and yield components among the different varieties, and to attempt to link physiological and anatomical traits to yield. The elite variety KWS Sassy outperformed other varieties across all metrics compared to landrace and wild varieties ( Figure 4A-H). The landrace Golden Archer delivered yield and yield component outcomes that were comparable with those for KWS Orwell on all variables measured and contrasted strongly with the performance of the landrace Hatif de Grignon (Figure 4A-H).   Figure 4A). This disparity in biomass was reflected in the number of ears, where KWS Sassy managed 27.5 ears plant −1 compared to just 11.0 in the case of Hatif de Grignon ( Figure 4C). KWS Sassy also stood in contrast to Hatif de Grignon and B3745 in grain number plant −1 delivering 2.4 and 3.3 times more grains, respectively, ( Figure 4D). Likewise, KWS Sassy was superior to the other varieties in TGW and grain mass plant −1 (Figure 4E,F). One reason for the relatively poor harvest outcomes for Hatif de Grignon could have been the number of infertile tillers ( Figure 4H), where over half the latter's tillers did not bear ears compared to under 2% in the case of KWS Sassy.

Photosynthetic Capacity of Elite and Landrace Barley Lines
To examine the possibility that heterogeneity in kinetic responses to changing light intensity was driven by biochemical processes, the photosynthetic capacities of a range of landraces and the elite variety KWS Irina were characterized ( Figure 5). By measuring assimilation at a range of light intensities from 0 to 1300 µmol m −2 s −1 PPFD, it was possible to assess the quantum efficiency (i.e., the initial slope of assimilation response to increasing light above 0 µmol m −2 s −1 PPFD) as well as the level of assimilation when light was saturating ( Figure 5A-C). The Central/Southern European landraces Dea and Hatif de Grignon had greater quantum efficiency (0.061 and 0.055 µmol µmol −1 , respectively) compared to the northern European cultivars Golden Archer and Eire 6-Row (0.029 and 0.033 µmol µmol −1 ; Figure 5C). Furthermore, Dea and Hatif de Grignon had higher A sat (31.5 and 28.9 µmol m −2 s −1 , respectively) compared to Eire 6-Row and Golden Archer (14.5 and 15.9 µmol m −2 s −1 , respectively). The response of A to changing Ci can be used to generate insights into biochemical activity in the chloroplast; notably the maximum rate of carboxylation of Rubisco (Vc max) and the maximum rate of electron transfer (Jmax). As Ci increases, so does assimilation which asymptotically approaches a maximum level ( Figure 6A). The response of A to changing C i can be used to generate insights into biochemical activity in the chloroplast; notably the maximum rate of carboxylation of Rubisco (V c max ) and the maximum rate of electron transfer (J max ). As C i increases, so does assimilation which asymptotically approaches a maximum level ( Figure 6A). Here, the both the Vc max ( Figure 6B) and J max ( Figure 6C) responses followed a similar pattern to those observed in the light response curves. Variety had a significant effect on Vc max (F (4, 37) = 2.63, p < 0.05), and the landraces Dea and Hatif de Grignon had the highest Vc max (92.2 and 78.5 µmol m −2 s −1 mol −1 ) compared to Golden Archer and Eire 6-Row (53.3 and 42.6 µmol m −2 s −1 mol −1 ). It was a similar pattern for J max , where there was also a significant effect of variety (F (4, 37) = 9.62, p < 0.05), and once again, Dea and Hatif de Grignon had higher rates of electron transport (142.9 and 116.2 µmol m −2 s −1 , respectively) than either Golden Archer or Eire 6-Row (75.9 and 71.6 µmol m −2 s −1 , respectively). There was also a close correlation between Vc max and A sat overall (σ = 0.41, p < 0.05) and between J max and A sat (σ = 0.93, p < 0.05).

Discussion
The aim of this study was to assess the extent of natural variation in anatomical and physiological characteristics among a diverse range of barleys. Despite many decades of breeding for yield [54,55], we did not fully confirm our expectation that the elite varieties (KWS Orwell, KWS Sassy and KWS Irina) would have higher assimilation rates and g s and faster g s kinetic responses to changes in light intensity leading to greater growth relative to landraces and wild barleys (Figures 1, 2 and 4) as shown for the elite wheat [38].
Previous research examining variation in stomatal anatomy has suggested a tradeoff between stomatal size and stomatal density [17,31,56,57], with species displaying higher stomatal density usually having smaller stomata, while lower stomatal density coincides with larger stomata. This relationship has been shown to hold in closely related species [17,56] and different environments [57], although this has not always been demonstrated across diverse species [39,58]. Our findings here do not support a relationship between stomatal density and size, and furthermore, we observed no relationship between anatomical characters and physiological responses (Figure 3, Supplemental Materials).
Previous studies have associated the speed of stomatal responses with the size of stomata, with smaller stomata facilitating rapid changes in pore area in a range of different species although the relationship may not be true of all species [17,34,39,43,44,56,57,59,60]. The majority of the above studies have focused on stomatal opening in response to environmental change, and the size-speed relationship has been queried by Elliott-Kingston et al. [60], who reported that the speed of closure was not related the stomatal length in a range of species including ferns, gymnosperms, monocots and eudicots. Crops such as barley have dumbbell-shaped stomata and it is well-established that such morphology promotes rapid stomatal responses compared with the more-common kidney-shaped anatomy [39,57,[61][62][63].
Rapidity in stomatal behavior has gained considerable attention recently due to the potential for greater co-ordination between A and g s which would optimize gaseous exchange under dynamic environmental conditions [15,30,31,39,44]. Slow stomatal opening to increasing light intensity for example has been reported to limit A by ca. 10% [39] while slow closure during a shade fleck could decrease water use efficiency by 20% [31]. Even though the mechanisms behind the rapidity of stomatal movements are still to be fully elucidated [31,34,56,[63][64][65], studies that have successfully manipulated the speed of stomatal response have resulted in improved CO 2 uptake and/or water use efficiency [63], and therefore this represents an unexploited opportunity to finding novel genetic traits for future breeding programs [55]. Recently Salter et al. [66] identified QTLs for dynamic photosynthesis in barley, with novel QTLs for stomatal conductance described, illustrating the potential of such an approach.
Here we demonstrate significant variation in the time constants for increases in g s (τ gs ) in response to a step change in light intensity, as well as overall steady state g s values, although no relationship between the two was found (Figure 2). Both the kinetic responses and light response curves demonstrate a clear correlation between A and g s (Figures 1 and 2, Supplemental Materials), with high g s values correlated with A, supporting the close relationship between A and g s [54]. Although dynamic responses are clearly important and have implications for the rates of g s and A realized in the field, photosynthetic capacity (which provides an indication of maximum potential A) can also demonstrate significant variation between crop cultivar [67,68] and has been correlated with yield [68]. Here we measured photosynthesis as a function of light intensity to assess variation in light use efficiency between cultivars, while photosynthesis as a function of C i permits the biochemistry of Rubisco and the Calvin cycle to be evaluated. As expected both curves showed a hyperbolic response (Figures 5A and 6A). However it is clear that Dea and Hatif de Grignon, both landraces, showed considerably higher rates of photosynthesis compared to the three other cultivars, including an elite line, which grouped together at lower rates. Between cultivars we also found significant variation in light-saturated rate of photosynthesis (A sat ), quantum efficiency (QE), maximum rate of carboxylation (V c max ) and maximum electron transport (J max ) rates determined from light and CO 2 response curves ( Figures 5 and 6). Interestingly, the greatest values in the above parameters were observed in the landraces Dea and Hatif de Grignon, while the elite cultivars grouped together at lower values.
McAusland et al. [69] report significant differences in photosynthetic rates and potential in wild relatives of modern wheat and suggested that this represents an underutilized source of genetic and phenotypic diversity, while the same is also true for landraces and older elite varieties [38,67] that have yet to be evaluated. The greater QE observed in Dea and Hatif ( Figure 5) suggests a greater efficiency in use of absorbed light in CO 2 fixation, as well as enhanced carboxylation capacity of Rubisco. This superior physiological performance could be due to underlying biochemistry [70], as well as anatomical differences or modifications in diffusional constraints. Irrespective of the mechanism they represent an unexploited pool of potential genetic targets for increasing photosynthetic capacity and efficiency even if the evidence for yield enhancement is as yet unproven.
A couple of caveats accompany the findings presented here. The ambient temperature in the polytunnel was somewhat variable, with a mean of 26.6 • C +/− 1.1 • C and a peak of 44 • C during the time gas exchange measurements were being taken (15-17 June 2016) while RH was also variable through the day (44.2% +/− 2.6%). Such stress is known to affect yield negatively in cereals [71] with adverse impacts of high temperatures not just during stem elongation, but also during booting and heading [72]. All cultivars would be affected by extreme temperatures, but we would expect that landraces bred for cool, moist climates such as Eire 6-Row would be affected more than other cultivars [73].
There was also the lack of additional N application in this study, which in modern varieties appears deleterious [74,75]. However, not providing additional N provided us with the opportunity to explore physiological variation and interesting traits in wild barleys and landraces that could be masked if modern N application rates were used. We have shown in the absence of a substantial push from N fertilizers, landraces can deliver very creditable yields that are in line with those seen in some elite cultivars (Figure 3).
Therefore the physiological and biochemical traits identified above offer the breeder interesting avenues to explore in the search for new germplasm [76,77]. Tiller numbers are known to be relatively invariant under selection, so finding a similarity across varieties may have been expected [52,78]. However, the strong yield performance of the landrace Golden Archer in relation to the elite cv KWS Orwell should be of real interest.
Furthermore, the contrast in light and C i responses and in stomatal density between the two landraces Golden Archer and Hatif de Grignon is a trade-off in anatomy and physiology that warrants further investigation, particularly since similar contrasts were consistently observed in the light and C i responses in the other landraces, Dea and Eire 6-Row.
Varieties which tend not to conserve water also deliver higher yields through improved temperature control, particularly around the mid-day canopy temperature peak [20,54,79]. Such non-conservative varieties of barley should have slower kinetic responses delivering larger yields which was broadly confirmed here [80]. Hatif de Grignon and Golden Archer had relatively higher time constants for opening stomata suggesting Agronomy 2021, 11, 921 13 of 16 a non-conservative phenotype compared to the elite cultivars, while yields in the former tended to be lower than KWS Sassy, but at a level similar to KWS Orwell.
We might have expected the Central European landrace Dea and the Southern European Hatif de Grignon [80][81][82] to be optimized for higher WUE i based on a need to conserve water compared to Golden Archer and Eire 6-Row which are typically found in the wetter British Isles [80]. As climactic conditions faced by UK and other farmers force changes in varietal selection [83], the results presented above point to candidates for breeding programs under a wide variety of possible environmental conditions. In our study, high photosynthetic capacity did not translate directly into yield, with the elite varieties outperforming the landraces and wild relatives in terms of tillers, biomass and grain mass [67]. However, the increase in performance was not as high as might have been expected. Yields for elite crop cultivars have clearly improved recent decades [4]. However, in the absence of additional N, as in this study, yield performance by the elite varieties varied and therefore the higher photosynthetic capacity from some of the landraces could potentially be beneficial in the future or in certain regions when inputs may become limited. While high photosynthetic capacity does not always translate into increased yield [67], it does provide an opportunity to breed for specific traits that could be extremely beneficial in future crop breeding programme [68], that will require the identification of new traits to produce crops with greater resilience or the ability to grow in certain regions or under different climate scenarios.

Conclusions
The results presented here offer opportunities for breeders to identify and include additional genetic material to improve diversity and reduce bottlenecks to yield improvement but that opportunity surprisingly may not occur through directly trading off stomatal size and density. The evidence in this study suggests that even relatively modern landraces have many features of interest including rapid stomatal kinetics and increased photosynthetic capacity. These are features that could be exploited for the improvement of elite lines through enhanced physiology, as agronomists increasingly tune varietal performance to regional environmental requirements. Data Availability Statement: The datasets used or analyzed during the current study, as well as the materials used, are available from the corresponding author on reasonable request.