Structural and functional leaf diversity lead to variability in photosynthetic capacity across a range of Juglans regia genotypes

Abstract Similar to other cropping systems, few walnut cultivars are used as scion in commercial production. Germplasm collections can be used to diversify cultivar options and hold potential for improving crop productivity, disease resistance and stress tolerance. In this study, we explored the anatomical and biochemical bases of photosynthetic capacity and response to water stress in 11 Juglans regia accessions in the U.S. department of agriculture, agricultural research service (USDA‐ARS) National Clonal Germplasm. Net assimilation rate (A n) differed significantly among accessions and was greater in lower latitudes coincident with higher stomatal and mesophyll conductances, leaf thickness, mesophyll porosity, gas‐phase diffusion, leaf nitrogen and lower leaf mass and stomatal density. High CO2‐saturated assimilation rates led to increases in A n under diffusional and biochemical limitations. Greater A n was found in lower‐latitude accessions native to climates with more frost‐free days, greater precipitation seasonality and lower temperature seasonality. As expected, water stress consistently impaired photosynthesis with the highest % reductions in lower‐latitude accessions (A3, A5 and A9), which had the highest A n under well‐watered conditions. However, A n for A3 and A5 remained among the highest under dehydration. J. regia accessions, which have leaf structural traits and biochemistry that enhance photosynthesis, could be used as commercial scions or breeding parents to enhance productivity.


| INTRODUCTION
Common walnut, Juglans regia L., is an important and widely grown agronomic species with major production areas concentrated in the northern hemisphere. Its natural range encompasses mountains from western China to central Asia (McGranahan & Leslie, 2009) and was extended by humans spreading the species throughout eastern and southwestern Europe from Central Asia (Leslie & McGranahan, 1998). Throughout its natural habitat, J. regia grows under a range of climatic conditions with mean monthly maximum temperatures ranging from −9 to +30°C and annual cumulative precipitation from 175 to 1150 mm during the growing season (Geospatial Data; Duke, 1978;FAO, 2021 [https://aquastat.fao.org/climate-informationtool/]; Figure S1). J. regia has been utilized to develop commercial scion cultivars (e.g., Chandler) and hybrid rootstocks with resistance to abiotic and biotic stresses (Kluepfel et al. 2015;Leslie, et al. 2015).
Global walnut production is highly dependent on the limited genetic diversity of the commonly used scions. For example, only four scion cultivars account for~80% of total yields. This results in orchard susceptibility to abiotic stresses, disease and pathogens. In California, the Chandler cultivar accounts for 50% of productive acreage (USDA-NASS, 2020) but has limited capacity in dealing with high temperatures and water deficits, which can also increase susceptibility to plant pathogens and low-quality kernel production (Grant & Shackel, 1998;Lampinen et al., 2005;Rosati et al., 2006). Wild germplasm collections can serve as a valuable resource to increase genetic variability to improve tolerance to abiotic and biotic stressors and crop productivity. The diverse J. regia collection at the U.S. department of agriculture, agricultural research service (USDA-ARS) National Clonal Germplasm Repository (NCGR) located in Winters, CA, USA, holds such potential. However, to date, this collection has not been exploited to identify genotypes with increased abiotic stress tolerance and physiological traits related to enhanced yield.
Photosynthesis is a key determinant of crop productivity and positively related to biomass accumulation and yield production (Faralli & Lawson, 2020;Fischer et al., 1998;Kruger & Volin, 2006;Long et al. 2006;Simkin et al. 2019). Photosynthetic CO 2 response curves (A n − C i ; net assimilation, A n vs. CO 2 inside the leaf, C i ) can be used to assess the biochemical and diffusive limitations that determine photosynthetic rates (Long & Bernacchi, 2003;Sharkey 2016). The biochemical limitations are determined from the maximum carboxylation rate of Rubisco (V cmax ), the maximum rate of electron transport (J max ) and the maximum rate of triose phosphate use (TPU), all derived from A n − C i curves. When combined with chlorophyll fluorescence measurements, A n − C i curves can also provide information on diffusive limitations associated with mesophyll conductance (g m ; Harley et al., 1992), which is a measure of the ease with which CO 2 diffuses from the substomatal cavity to the site of carboxylation inside chloroplasts. Recently, these biochemical characteristics have been used to evaluate germplasm in crop breeding programmes (De Souza & Long, 2018;De Souza et al., 2020).
Identifying variation in diffusive limitations which are strongly linked to g m and leaf structure can also be used in breeding programs to improve photosynthesis (Tomás et al., 2013;Tosens et al., 2012). g m involves a complex pathway and a series of resistances in both the gas and liquid phases (Flexas et al., 2008;Tosens & Laanisto, 2018), and is impacted by various leaf structural traits including intercellular airspace (IAS) volume (i.e., porosity), mesophyll surface area exposed to the IAS (SA mes /V mes ), mesophyll cell diameter and density, and cell wall thickness (Evans, 2021;Flexas et al., 2008Flexas et al., , 2012. Leaf structure of some species exhibits plasticity in response to the growth environment (Salk, 2012), resulting in functional variation, which can help optimize resource use (Muir et al., 2017;Wright et al., 2004). Any inherent variation in leaf structural and physiological traits, as a function of the habitat environment to which it has been adapted, may play an important role in regulating photosynthetic capacity in J. regia accessions. We also recently found that two Juglans spp. exhibit changes in mesophyll structure under dehydration associated with changes in cell volume, orientation and arrangement that increases porosity (Momayyezi et al., 2022). Desiccation influences CO 2 diffusion and water relations as mesophyll cell turgor changes (Buckley et al., 2017;Scoffoni et al., 2014). How genotypic diversity in mesophyll cell packing and distribution across J. regia accessions may link with photosynthetic performance and susceptibility to drought is yet to be investigated. X-ray microcomputed tomography (microCT) provides an in-depth assessment of leaf mesophyll traits (i.e., porosity and tortuosity) and cells orientation and geometry (Earles et al., 2018(Earles et al., , 2019Lundgren & Fleming, 2020;Théroux-Rancourt et al., 2020).
In this study, we combined gas exchange physiological analysis with microCT imaging of leaves to explore: (1) the photosynthetic capacity of numerous J. regia wild accessions originating from habitats with varied climatic conditions; (2) links between leaf structural diversity and physiological features that enhance photosynthetic capacity; and (3) whether genotypic differences hold up under water stress conditions. Based on our previous observations for two Juglans species (Momayyezi et al. 2022), we hypothesize that greater A n will be associated with thicker leaves and higher mesophyll porosity and gas-phase diffusion (g IAS ; Han et al., 2019;Tomás et al., 2013). We also expect accessions originating from lower latitudes, characterized by warmer and wetter habitats, would exhibit higher inherent A n to support growth through enhanced biochemical activities as predicted for plants under warmer conditions (Moore et al., 2021;Ruiz-Vera et al., 2018). We anticipate that greater A n would be concurrent with higher IAS and conductance as temperature and precipitation are known to strongly influence functional diversity of photosynthesis across species (Harrison et al., 2020;Ordonez & Svenning, 2017).

| MATERIALS AND METHODS
Stems were collected from 11 genetically unique J. regia accessions ( Figure 1a,b) at the USDA-ARS-NCGR in Wolfskill Experimental Orchard, Winters, CA, USA, to use as scion and were grafted by Sierra Gold Nursery onto a commonly used commercial rootstock, RX1 (Juglans microcarpa × J. regia). We used a common rootstock to eliminate any own-root effects and to simulate conditions for a commercial walnut orchard setting, where rootstocks are commonly used. The grafted saplings were repotted and transferred to the Armstrong lathe house facility at the University of California, Davis, in June 2019, and kept under natural light, temperature and relative humidity (~50%). The budburst for accessions occurred between April and May.
Measurements were initiated on nonstressed plants in August 2019 under well-watered conditions for all accessions. The measurements on the same plants were repeated under dehydrated condition through a gradual dry down procedure using the methods as described by Knipfer et al. (2020). Briefly, water loss from pots via transpiration and soil evaporation were quantified during the experiment by weighing pots to calculate the required amount of water needed per pot under each treatment, daily. After completion of measurements in nonstressed conditions, water application was reduced to 75% of full irrigation during the first week and then to 50% of full irrigation in the second week of drying. This watering regime was then maintained for the dehydration treatment until completion of the experiment. Plants were maintained under ambient natural light with a~15 h photoperiod during the experiment, maximum temperature of 35°C during day and minimum of 15°C during night, in 2.65 L pots containing a 40% pine bark, 40% sphagnum peat moss and 20% vermiculite, and a thin layer of the slow-release fertilizer (Osmocote Smart-Release Plus) was added to the top soil at the plantation time. The two irrigation treatments were maintained for 2 weeks before the measurements.

| Illumina sequencing and phylogeny tree construction
Genotyping-by-sequencing libraries were prepared by restricting genomic DNA with HindIII-HF (NEB) with Ampure XP bead (Beckman-Coulter) cleanup. Sequence data were obtained on an Illumina Hiseq 4000 instrument at the UC Davis Genome Center.
Single nucleotide polymorphism (SNP) calling was performed using the TASSEL GBS pipeline (Glaubitz et al., 2014) with alignment to the Chandler genome (Marrano et al., 2020), resulting in 13,320 polymorphic SNPs. The phangorn package in R (R Core Team, 2017; Schliep et al., 2017) was used to construct a neighbour-joining tree from the distance matrix.

| Photosynthetic measurements
Net assimilation rate (A n ), stomatal conductance (g s ) and the IAS CO 2 concentration (C i ) were measured on one leaflet of the leaf from each accession in five replications using a LI-COR 6800 system fitted with 6800-01A fluorometer. To provide a general leaf developmental assessment before the measurements, the leaf plastochrone index (LPI) was calculated as described by Erickson & Michelini (1957). The 4th or 5th leaflet was counted down from the tip leaflet on the leaves at LPI between 5 and 6. All measurements were done under photosynthetic photon flux density (PPFD) of 1500 µmol m −2 s −1 (10% blue vs. 90% red), a saturating PPFD for all accessions ( Figure S2), chamber temperature at 25°C, ambient chamber CO 2 concentration (C a ) at 400 (µmol mol −1 ), flow rate at 500 (µmol air s −1 ) and vapour pressure deficit between 1.5 and 2.0 kPa. To calculate the F I G U R E 1 (a) Geographic distribution map for 11 J. regia accessions. Location data were found in the USDA-ARS National Plant Germplasm System for individual accessions (Germplasm Resources Information Network, GRIN). (b) Unrooted neighbourjoining tree for 11 J. regia accessions. Thirteen thousand three hundred and twenty polymorphic single nucleotide polymorphisms (SNPs) among these 11 trees were discovered by Illumina sequencing for the tip leaflet and used to construct NJ trees in R using the "phangorn" package. *Genotype data from C4 29, different individual but same accession as A9 (C4 28; see Table S1) was used, as data for original A9 was not available.

| Calculation of g m by chlorophyll fluorescence and of CO 2 concentration in the chloroplast (C c )
The "variable J method" was used to estimate g m based on calculation of electron transport rate (J flu ) from measurements of chlorophyll fluorescence (Bongi & Loreto, 1989;Harley et al., 1992): where β (0.5 for C 3 plants) is the fraction of absorbed quanta reaching photosystem II (Bernacchi et al., 2002). The leaf absorbance, α, was previously measured for Juglans species to be 85.3% based on the F I G U R E 2 Leaf cross sections from representative scans for 11 J. regia accessions under well-watered condition obtained using X-ray microcomputed tomography.
We used this approach previously for two Juglans species with contrasting leaf anatomy (Momayyezi et al., 2022) and applied it to the F I G U R E 3 Photosynthetic CO 2 response curves were constructed using Sharkey's fitting calculator version 2.0 (Sharkey 2016), averaged for 5 replications in 11 J. regia accessions under well-watered condition. A − C i curves are shown with coloured circles and error bars measured directly. A − C c curves were used to generate V cmax and J max , and averaged over five replications for each accession (±SE, n = 5). Black dashed horizontal line indicates assimilation rate at saturating CO 2 (A max ) at triose phosphate use (TPU) limitation state, and Rubisco and RuBP regeneration limitations are indicated for each accession by red and black curves, respectively. Colour scheme is consistent with accession numbers presented in Figure 3 and in Table 1 11 J. regia accessions in this study. We assumed any potential small changes in C i * and R d under dehydration will be negligible, and the same values measured for Juglans species under well-watered conditions were used to estimate g m under dehydration. Having obtained g m by the chlorophyll fluorescence method, the CO 2 concentration in the chloroplast (C c ) was found according to Harley et al. (1992): Before running A n − C i curves, the seal of the fluorometer chamber was tested for CO 2 leakage by running a full curve at different C a concentration in empty cuvette and the measurements were corrected using the Licor's equation (Bernacchi et al., 2002;Flexas et al., 2007). The A n and corresponding C i values for each accession in five replications were introduced to the Sharkey's fitting calculator version 2.0 (Sharkey 2016), a Farquhar-von Caemmerer-Berry-based model to estimate V cmax , J max and g m using differences in CO 2 partial pressure between C i and C c and predict assimilation rates (A) at the limiting states of Rubisco, RuBP-regeneration and TPU (Sharkey et al., 2007). We used photosynthesis at the TPU state to set the maximum assimilation rate (A max ) at saturating CO 2 as in Sharkey et al. (2007).
g m obtained from the chlorophyll fluorescence method was verified against g m found using A n − C i method ( Figure S3). g m obtained from the chlorophyll fluorescence method was used to present correlations, prepare figures and tables and make statistical comparisons for each accession under both well-watered and dehydrated conditions.

| Leaflet water potentials
Leaflet water potential (Ψ leaflet ) was measured using a pressure chamber (PMS Instrument Company, Model 1505D) immediately after gas exchange measurements between 10 AM to 12 PM (Williams & Araujo, 2002). A leaflet directly opposite the one used for gas-exchange measurements was used to measure leaflet water potential from each accession at five replications. Each leaflet was cut at petiolule base and bagged for 10-15 min to allow equilibration within the leaflet. Then, the petiolule was placed inside the pressure chamber gasket. Chamber pressure was increased slowly until the balancing pressure was reached.

| Leaf mass and nitrogen content per area
Immediately after gas-exchange measurements, the leaflet adjacent to the one used for gas exchange was punched and collected discs were dried at 72°C oven for 48 h. Using the leaflet discs area and dry weight, leaf mass per area (LMA) was calculated for each accession using five replications. Later, the dried samples were ground and 4.5-5.0 mg of them were encapsulated in tin cups and sent to the Stable Isotope Facility, at University of California Davis for total nitrogen content elemental analysis. Using LMA data, the nitrogen content results per unit mass were converted to leaf nitrogen content per unit area (Leaf N).

| X-ray microCT imaging
Leaves from each accession and treatment in five replications were scanned using X-ray microCT at beamline 8.
as described in Earles et al., (2018). The L geo and L Euc were mapped and quantified for all voxels along the mesophyll surface and τ was calculated for the whole three-dimensional (3D) image array as in Earles et al. (2018). Then, leaf-level tortuosity (τ leaf ) was calculated as the mean of τ values at the edge of mesophyll cells. The lateral path lengthening λ (m m −1 ) was calculated using L Euc and a second distance map as described by Earles et al. (2018) to measure the shortest unobstructed distance in a straight line between the abaxial epidermis and all points along the mesophyll surface L epi (Legland et al., 2016): Similarly, leaf-level lateral path lengthening (λ leaf ) was then calculated as the mean of λ values at the edge of mesophyll cells.
Vein volume relative to the leaf volume ratio (V vein /V leaf , m 3 m −3 ) was calculated as a fraction of total leaf volume as described by Trueba et al. (2022) from the same scans used for θ IAS , τ and λ leaf calculations.

| IAS conductance and stomatal density
The τ leaf , λ leaf and θ IAS were used to calculate leaf-level IAS conductance (g IAS ), where D m is the diffusivity of CO 2 in air (m 2 s −1 ). Diffusion path length in gas phase was equal to half of the mesophyll thickness (L mes ) for hypostomatous leaves (Earles et al., 2018;Niinemets & Reichstein, 2003;Tomás et al., 2013):  Figure S4).
To quantify the stomatal density, the grid reconstructions were used at the paradermal direction. The stomata density and size were measured in 0.04 mm 2 of leaflet area for all accessions (n = 5).

| Climatic data for accessions' native habitats
Coordinates for each accession's native habitats were extracted from the USDA-ARS Germplasm Resources Information Network database (https://npgsweb.ars-grin.gov/gringlobal/search). Temperature and precipitation data were obtained from FAO climate information tools (https://aquastat.fao.org/climate-information-tool/) for each of these native habitat locations ( Figure S1). The temperature and precipitation seasonality were calculated based on the ratio of the SD for monthly mean temperature or precipitation to the monthly mean temperature or precipitation (known as the coefficient of variation) and multiplied by 100. The temperature or precipitation values for each month were averaged and the SD was calculated using the mean of the 12 months (Donnell & Ignizio, 2012).

| Statistics
Linear regression and Pearson correlation coefficients were used to examine relationships between latitude, temperature and precipitation seasonality, frost-free days and A n , A max , g s , g m , g IAS , L leaf , θ IAS , λ leaf , τ, V vein /V leaf , WUE i , Leaf N, LMA, J max , V cmax and Ψ leaflet using GraphPad prism 9 software (GraphPad Software, Inc.). Paired t-test was used to check for systematic differences between the chlorophyll fluorescence and A n − C i curve methods for estimating g m and C c . Mixed linear models were used to compare relative changes in percent for A n , g m , g s , L leaf , g IAS , θ IAS and Ψ leaflet under dehydration for all accessions using SAS 9.4 (SAS Institute Inc., 2013). The p required for significance (0.002) was adjusted by dividing α (0.05) by the number of comparisons per test (25 here).
Logarithm or squared transformations were performed to meet normality and equal variance assumptions where needed.
Similarly, A n was extracted from A n − C i curves at C a of 40.4 (Pa) and was positively correlated with chlorophyll fluorescence g m (p < 0.0001) and g s (p = 0.034; Figure 4) across the accessions.

| Climate-driven photosynthetic capacity
Photosynthetic capacity and associated leaf physiological and anatomical characteristics for the accessions were associated with the climatic conditions in native habitats. Accessions from habitats with lower temperature seasonality in lower latitudes had higher g ias and Leaf N (p ≤ 0.05), which may support higher their A n through greater g m and θ IAS ( Figure 5; p = 0.05; Figure S5). Despite nonsignificant relationships, parallel decreases in g m , g IAS , Leaf N and θ IAS (p > 0.05; Figure S5) may suggest a pattern for decline in mesophyll CO 2 diffusion rate with latitude, while they accumulate more LMA (p = 0.015; Figure 5). Higher variability in precipitation seasonality and more frost-free days were significantly related to increased g IAS concurrent with lower stomatal density ( Figure 5). Decreases in LMA and frost-free days were associated with increases in Leaf N ( Figure 5).

| Responses under dehydration
As expected, dehydration impaired photosynthesis and altered leaf structure with reduced g s , g m and L leaf and increased θ IAS and g IAS in all accessions. The percent reduction in A n was significantly correlated with percent reductions in g s (p = 0.007), g m (p = 0.002) and Ψ leaflet (p = 0.002) under dehydration (Table 1). Accessions A3, A5 and A9, which had the highest A n and A max under well-watered conditions, exhibited the greatest percent reductions amongst accessions in these parameters under drought stress (i.e., >50% reduction for all 3; Table 1). However, the absolute values of A n for accessions A3 and A5 were not significantly lower than other accessions under dehydration, whereas they were among the lowest for A9 under dehydration (Table S2). The reduction in L leaf was linked with decreases in g ias (p = 0.02) and g m (p = 0.08). The concurrent reduction in Ψ leaflet was significantly correlated with percent decline in g s (p = 0.004; Table 1). Under dehydration, absolute A n , g m , θ IAS , g IAS and Leaf N remained negatively correlated with latitude (p < .05, Net assimilation rate (A n , μmol CO 2 m −2 s −1 ) relationship with mesophyll conductance obtained from chlorophyll fluorescence method (g m , mol CO 2 m −2 s −1 ), stomatal conductance (g s , mol m −2 s −1 ), mesophyll porosity (θ IAS , m 3 m −3 ), leaf thickness (L leaf , μm), leaf nitrogen per unit area (Leaf N, mg cm −2 ) and leaf mass per unit area (LMA, mgcm −2 ) in 11 J. regia accessions using mean values (±SE, n = 5). A n was extracted from A n − C i curves at C a of 40.4 Pa. Colour scheme is consistent with accession numbers presented in Figure 3 and in Table 1. [Color figure can be viewed at wileyonlinelibrary.com] significantly higher A n , A max and CO 2 diffusion capacity under the well-watered condition with the highest V cmax and J max . Greater photosynthetic capacity is typically linked to carboxylation capacity via increased Rubisco protein abundance and activity (Hikosaka & Shigeno, 2009;von Caemmerer & Farquhar, 1981;Kattge et al., 2011).
Higher photosynthetic capacity was strongly linked to leaf thickness and mesophyll structure and supported by higher Leaf N. As expected, J. regia accessions with higher A n had thicker leaves and exhibited greater θ IAS and g m (Figure 4). This agrees with previous findings for J. regia cv. Chandler (Momayyezi et al., 2022). J. regia accessions with higher A n had greater A max , V cmax and Leaf N, suggesting higher carboxylation capacity and performance (Figures 3 and S6). Although previous work reported a positive relationship between A n with LMA (i.e. across Quercus ilex provenances; Peguero-Pina et al., 2017), we found J. regia accessions exhibited greater A n with lower LMA. Increasing cell density would reduce mesophyll surface area exposed to IAS as a result of high cell packing and could also be impacted by cell wall thickness (Niinemets et al., 2009a;Tomás et al., 2013;Tosens et al., 2012). A more porous mesophyll and thicker leaves with shorter λ leaf ( Figure S7), contributed to higher A max across J. regia accessions ( Figure S6), highlighting the fact that thickness and cell density may not change in the same direction (Niinemets, 1999;Syvertsen et al., 1995). Increases in leaf porosity are known to reduce diffusive resistance and lateral path lengthening in other species (Earles et al., 2018). Additionally, leaf mesophyll geometry and IAS are known to impact stomatal patterning, and water supply to replace losses due to transpiration Zwieniecki et al., 2007).
+15.5 ± 0.9 (A8,10) +28.5 ± 3.6 (A1,7,10) A 10 −44.9 ± 1.2 (A1,2,4,5,6,7,9,11)  suggesting a potential positive impact of irrigation on leaf performance by improving CO 2 diffusion. The accessions with lower diffusion resistance (e.g., more porous leaves) benefit from a higher CO 2 diffusion capacity, and exhibit inherently higher performance and biochemical activities (i.e., V cmax ) under the common garden conditions of the experiment. Comparing A − C i and A − C c curves provides insights into the limitations imposed by mesophyll g m on assimilation rate associated with the drawdown of CO 2 from IAS (C i ) to the chloroplast (C c ; Figure 3). J. regia accessions from lower latitude exhibit enhanced performance through increasing V cmax at low CO 2 concentrations, found by the net assimilation rate response to C c (A − C c curve), when Rubisco is limiting A (Figure 3). That was in line with greater leaf nitrogen content for those accessions ( Figure 4) and more investment toward Rubisco activity (Niinemets et al., 2009b;Sharkey et al., 2007;Warren & Adams, 2004), enhancing their A n and CO 2 diffusion capacity at lower CO 2 concentrations ( Figure 4).
Greater maximum electron transport rate (J max ) also found using the A − C c curve, where RuBP-regeneration is limiting A, was associated with lower diffusion resistance through mesophyll (e.g., higher θ IAS ) and higher g m , and higher enzymatic activity during CO 2 fixation and carbohydrate formation (i.e., Calvin cycle) for low latitude J. regia accessions. A max was strongly associated to V cmax and J max ( Figure S6), supporting greater biochemical and diffusional capacity for low latitude accessions which exhibit thicker leaves with greater θ IAS and g m , and more nitrogen accumulation within leaves.
On the other hand, the photosynthetic performance seems to be unrelated to the phylogenetic history. Accessions with greater A n and A max , such as A3 and A9, share a close evolutionary background with low-performance accessions, like A4 and A11, but not with each other ( Figure 1B). Therefore, unlike studies that reported strong phylogenetic support for water stress resistant traits like xylem cavitation vulnerability in stem and root (Willson et al., 2008) or vein development, patterning and hydraulic conductance , our results suggest that geographical variability is more strongly linked with differences in photosynthetic rate in Juglans accessions. Therefore, high intraspecific variability in leaf anatomical characteristics and high variation in photosynthetic performance between J. regia accessions are in line with their establishment across a wide range of habitats with contrasting environmental conditions. This is consistent with results of He et al. (2018b), who showed that species with lower variability in traits like specific leaf area have narrower habitat range, contrasting other studies with strong influence of phylogenetic history on traits variation (Homeier et al., 2021).
Leaf nitrogen content per unit area (Leaf N) was negatively related to latitude and supported higher A n with greater biochemical activity. Increases in A n related to changes in Leaf N and chlorophyll content have been reported for Populus balsamifera and P. angustifolia populations as an adaptive response to growing season length (Soolanayakanahally et al., 2009;Kaluthota et al., 2015). Latitudinal variation in photosynthetic variables have been reported more broadly across various species, however, the patterns were opposite in Populus spp. For example, latitudinal increases in A n and for P. trichocarpa genotypes was accompanied by greater g m through higher carbonic anhydrase activity and aquaporins functioning that were attributed to growth under shorter growing season in northern habitats (Gornall & Guy, 2007;McKown et al., 2014aMcKown et al., , 2014bMomayyezi & Guy, 2017.

| Dehydration-induced responses
As expected, dehydration negatively impacted photosynthesis in all accessions, but some exhibited greater reductions in A n with decreases in g m , g s and Ψ leaflet (p = 0.007). Accessions A3, A5 and A9 exhibited similar absolute values to each other (Table S2) and had the highest photosynthetic capacity than other accessions, with the highest % reductions due to drought. A reduction in photosynthesis under dehydration was associated with decreases in PSII efficiency, which could be due to increases in photorespiration associated with increased resistance to CO 2 diffusion through stomata and mesophyll (Busch 2020;Neto et al., 2017;Sharkey 1988). Even under dehydration, accessions A3 and A5 from lower latitudes maintained A n , g m , θ IAS , g IAS and Leaf N higher than other accessions ( Table 2), suggesting that these accessions hold potential for commercial production without increasing susceptibility to stress in absolute terms. Note: Pearson correlation coefficients between the absolute values of the physiological and anatomical variables and climatic data for 11 J. regia accessions under dehydration treatment using mean values (±SE, n = 5). Bold indicates significance at p < 0.05 and * indicates significance after Bonferroni corrections (p < 0.0025).
LEAF STRUCTURE-FUNCTION AND PHOTOSYNTHETIC CAPACITY | 2361 Dehydration reduced leaf thickness and increased θ IAS and g IAS by shrinking mesophyll cell size more than IAS (Table 1 and S2); this is consistent with our earlier observations for J. regia cv. Chandler and J. microcarpa (Momayyezi et al., 2022). In the current study, g IAS increased by 13-35% in different J. regia accessions under dehydration. However, its contribution to g m (calculated as described by Niinemets & Reichstein, 2003) was 7-23% under well-watered but decreased to 3-8% under drought across accessions, which is within the expected limiting range (3-37%) for woody perennial species with hypostomatous leaves (Harwood et al., 2021;Niinemets & Reichstein, 2003;Parkhurst & Mott, 1990;Tomás et al., 2013;Tosens et al., 2012). Additionally, reductions in the g IAS contribution to g m was closely and positively related to A n (R 2 = 0.63; p =;0.003) and A max (R 2 = 0.37; p = 0.044), suggesting reduced gas phase diffusion under stress may decrease photosynthesis further by limiting CO 2 diffusion in liquid phase through chloroplast re-positioning and activity of carbonic anhydrases and aquaporins (Evans et al., 2009;Miyazawa et al., 2008;Momayyezi et al., 2020;Tholen et al., 2008;Tomás et al., 2013).

| Conclusion
We found that photosynthetic capacity in J. regia accessions was associated with leaf anatomical and biochemical components that impact CO 2 diffusion. Leaves with greater porosity and g IAS contribution to g m exhibited the highest photosynthetic capacity at ambient and saturating CO 2 . Improved photosynthesis was supported by increased carboxylation capacity and leaf nitrogen accumulation. Higher photosynthesis across accessions was associated with frost-free days and precipitation and temperature seasonality patterns in low-latitude native habitats. Although J. regia has a limited resilience under dehydration, two of the low-latitude accessions (e.g. A3 and A5) with the highest inherent photosynthetic capacity, sustained performance under stress. These accessions hold promise for high productivity and use in breeding programs for commercial walnut production ( Figure 2).