The Influence of Leaf Anatomical Traits on Photosynthesis in Catimor type Arabica coffee

Leaf photosynthesis is largely determined by anatomical features. This study aimed to reveal the quantitative effects of the anatomical structure of Coffea arabica leaves on photosynthesis. Pearson’s correlation and path analysis were used to explore the correlation between leaf structure and photosynthesis. To calculate the comprehensive evaluation value of the correlation between leaf anatomical traits and photosynthetic parameters, the Criteria Importance Through Intercriteria Correlation (CRITIC) method was used to obtain the objective weight of each photosynthetic parameter. The study revealed that leaf anatomical traits were highly significant ( p <0.01) and correlated with photosynthetic parameters, suggesting that anatomical traits greatly influenced photosynthesis in C. arabica leaf. Similarly, path coefficient analysis strongly showed direct and indirect correlation of photosynthetic capacity of stomatal conductance of the leaves of C. arabica . The results of the comprehensive evaluation also indicated that leaf thickness (LT) and stomatal density ( d ) were the anatomical characteristics most closely related to photosynthesis. In these results, understanding the effects of the anatomical structure of coffee leaves on photosynthesis, may provide useful information for coffee breeding programs and the management of coffee plantations to optimize Photosynthetic capacity.


Introduction
Coffee (Coffea arabica) considered an important tropical crop, is one of the most traded agricultural products in the world.Under suboptimal growing conditions, shade trees, which could reduce the stress on coffee by ameliorating adverse climatic conditions, have been recommended [1,2] .Yunnan province, located in Southwest China, is the most suitable region for coffee production in China [3] .Although the mountainous landscape and mild climate in this region are well-suited to coffee production, many areas are suboptimal for coffee cultivation [3,4] .Thus, agroforestry systems have been recommended for these areas, especially at altitudes below 1000 m [4] .
Since leaves are specialized organs that are responsible for interception of light and uptake of CO 2 through the stomata for photosynthesis [5] , light and CO 2 availability in the mesophyll of plant leaves largely determin the net carbon assimilation rate (A) of leaves [6] .To analyze the response of photosynthesis to irradiance, the light response curves (LRCs) were used to describe the relationship between photosynthesis and light intensity [7] .It is an important method to elucidate the response mechanism of photosynthesis and evaluate photosynthetic efficiency by fitting of light-response model [5] .Often, LRCs are mathematically described using empirical models.The light response curves (LRCs) can be used to derive key photosynthetic parameters including the maximum net photosynthetic capacity (P max ), initial and apparent quantum yields (IQY and AQY), light compensation point (LCP) and light saturation point (LSP), as well as the dark respiration rate (R D ).In addition, mech-anistic information about the key biochemical parameters of carbon assimilation can be estimated by fitting A/Ci curves with the mechanistic Farquhar-von Caemmerer-Berry (FvCB) biochemical model described by Farquhar et al. (1980)  [8−10] .This model has been widely used to simulate CO 2 assimilation and applied for decades [9,11] .It clearly explains the physiological properties of photosynthesis of intact leaves [5] , and the parameters reflecting the biochemical properties of leaves as being easier to estimate from this model with gas exchange measurements [8] .
The structural organization of plant leaves can either facilitate or constrain the photochemical reactions underlying photosynthesis [12,13] .Therefore, the leaf structural changes may affect photosynthetic efficiency [7,14] .However, although a variety of ecological studies have revealed a correlation between leaf structural parameters and photosynthetic performance [15−17] , there were contradictory differences between foliage structural and functional variables in the previous work [18] .Therefore, in order to understand the effect of leaf structure on each of the photosynthetic performance parameters, it is necessary to analyze the relationship between leaf structure and leaf photosynthesis and to comprehensively evaluate the importance of each component in leaf photosynthesis.Furthermore, evaluation of the relationship between leaf structure and leaf photosynthesis could provide ecophysiological evidence for understanding the structural properties of leaves and how they affect photosynthetic performance.Even though a number of studies analyze the effects of leaf structure on photosynthetic performance [12,15] , there has been limited synthetic analysis of the correlations between leaf structural

Study site characterization and plant material
The experiment was conducted in a shade coffee plantation of the Dehong Tropical Agriculture Research Institute of Yunnan (DTARI), Ruili, southwest China.
The study was conducted from 2019 to 2021.The experimental site was situated at Latitude 24.025 south and longitude 97.855 west at an altitude of 890 meters above sea level with an average temperature of 22.2 °C and precipitation of 1260 mm.The soil type in the plantation is acidic lateritic red soil.
The cultivars used for this study was 'Catimor' which is an interspecies hybrid derived from the cross between Caturra (C.a arabica L.) and Timor (C.arabica × C.canephora Pierre ex Froehner) The coffee trees were approximately five years old and spaced 1 × 2 m.Shade trees (Dimocarpus longan Lour.) were planted in East to West row orientation at 12 × 4 m dimension.
Fourteen plants with consistent growth were selected in the coffee plantation.Three leaves were taken from fourteen cohort plants in different directions for measurement of leaf gas exchange parameters.Measurements were taken from the third or fourth completely expanded leaf pair at the apex of the plagiotropic (lateral) branches, located in the upper third of the plant.

Gas exchange measurements
The leaf gas exchange parameters ( net CO 2 assimilation rate (A), stomatal conductance to water vapor (g s ), transpiration rate (E), internal CO 2 concentration (C i ), vapor pressure deficit (VPD), and leaf water use efficiency (WUE=A•E −1 ) and environmental indicators (leaf temperature (TL), relative humidity (RH) and ambient temperature (TA)) were determined simultaneously by using a portable photosynthesis system CIRAS 3 (PP Systems, USA).Leaf water use efficiency (WUE), representing the units assimilated for CO 2 per unit of water lost through transpiration, was determined via the user manual version 1.06.All measurements were taken on sunny days from 09:00 to 11:30 and 15:00 to 18:00.The referenced CO 2 concentration was set at 390 ppm using a buffer bottle, the temperature at 27 °C, flow rate at 300 µmol•s −1 , relative humidity at 50% ± 10%, and irradiance at 800 µmol m −2 •s −1 .The leaf photosyntheticLRCs) and internal CO 2 response (A/C i ) curves were constructed by using leaf gas exchange measurement.To collect data for building LRCs, net photosynthesis was measured at 0, 50, 100, 150, 200, 250, 300, 350, 400, 500, 600, 700, 800, 1,000, and 1,200 µmol•m −2 •s −1 (PAR) with a cuvette environment that maintained leaf temperature at 27 °C, reference CO 2 concentration at 400 ppm.The gas exchange rate usually stabilizes within approximately 3-4 min after the leaf tissue is clamped in the leaf chamber.Data for the construction of the A/Ci curves were collected immediately following the LRC measurements from the same sampled leaf.P n was measured at CO 2 concentrations of 400, 300, 200, 100, 150, 250, 350, 450, 550, 650, 750, 850, 950, 1,050, 1,250, 1,450 and 1,650 µmol•m −2 •s −1 at photosynthesis saturating irradiance 800 µmol•m −2 •s −1 ，which light intensity is sufficiently high to saturate the photosynthetic machinery without causing photoinhibition.A total of 41 LRC and 42 A/C i curves were generated.

Measurement of leaf anatomical traits
Following the gas exchange measurement, five samples from leaves of each plant were harvested, sealed in polyethylene bags, and immediately brought back to the laboratory.Two segments (3 mm × 5 mm) were cut from each leaf halfway between the leaf apex and base.The segments were dehydrated in sucrose solutions, followed by the inclusion with Gum Arabic aqueous solution.Transverse sections (7-10 µm) thickness was cut using a rotary microtome (Leica CM3050S, Germany) and stained with toluidine blue.Surface impressions were made manually and stained with0.1% safranine aqueous solution.For each leaf, five transversal and surface sections were selected systematically to measure the leaf thickness (LT), upper cuticle thickness (UCT), upper epidermal thickness (UET), palisade tissue thickness (PT), spongy tissue thickness (ST), lower epidermal thickness (LET), lower cuticle thickness (LCT), stomatal density (d), and guard cell length (l) and width (w) using a binocular microscope binocular with 10 × and 40 × objectives (Leica DM5000 B).Subsequently, five micrographs per section were captured with a photomicroscope equipped with a Leica DFC500 camera (Leica Microsystems, Rijswijk, The Netherlands).Images were analyzed using Leica Leica LAS X software.Total stomach pore area (TSP) was calculated as d × 1/4π × l × w (7), assuming pore was an ellipse [16,19,20] .Here d is the stomatal frequency calculated by dividing the stomatal count by t field area, l is the guard cell length, and w is the guard cell width.

Calculations and statistical analysis
−23] .The modified rectangular hyperbola model is described by Equation (1).
where P n is the net photosynthetic rate (µmol•CO 2 •m −2 •s −1 ), I is the photosynthetic photon flux density (µmol•m −2 •s −1 ), is compensation irradiance, β and γ are the coefficients which are independent of I, and α is the absolute value of slope between I = 0 and I = I c .
For I = 0, the quantum yield at this point is defined as intrinsic quantum yield (IQY), which is given by Equation (2).
The apparent quantum yield (AQY) ignored the Kok effect is obtained by Equation (3).
For I = 0, the rate of dark respiration (R d ) is Equation (4).
The saturation irradiance I sat is obtained by Equation (5).
The maximum photosynthetic rate P nmax is given by Equa- I n P r e s s tion (6).
A/Ci curves were fit by uploading the datasets to the website (http://Leafweb.ornl.gov) to estimate the parameters (R d , V cmax , J max , TPU, Γ * , g i ).In this study, the A/C i curves were taken at the saturation light level.Thus, the rate of electron transport (J) was to assume the maximum rate of electron transport (J max ).
The relationship between photosynthetic capacity and leaf anatomical traits was analyzed by calculating the Pearson correlation coefficient.The correlation coefficients were partitioned into direct and indirect effects using path analysis.The path diagram in Fig. 1 was utilized to examine the causal pathways between nine selected leaf anatomical traits and each photosynthetic parameter.In Fig. 1, the direct effects of leaf anatomical traits on photosynthetic parameters are represented by single-headed arrows, while coefficients of intercorrelations between leaf anatomical traits are represented by double-headed arrows.
Path coefficient analysis was calculated as the method described by Dewey and Lu (1959) [24] .
where r ij represents the simple correlation coefficient between the independent character (i) and dependent character (j) as measured by the correlation coefficient, P ij is the component of direct effects of the independent character (i) and dependent character (j) as measured by the path coefficient and, ∑r ik p kj is the summation of components of the indirect effect of a given independent character (i) on the given dependent character (j) via all other independent characters (k).An uncorrelated residue (U) was estimated by using the Equation (9).
where R 2 is the coefficient of determination, R 2 =∑r ij p ij .
The multiple criteria decision analysis (MCDA) was used to evaluate comprehensive evaluation of the relationship between leaf anatomical traits and photosynthetic parameters.The objective weight of each photosynthetic parameter (W j ) is derived by using the Criteria Importance Though Intercrieria Correlation (CRITIC) method, which reflects the relative importance by applying the comparative and conflict information The Influence of Leaf Anatomical Traits on Photosynthesis in Catimor type Arabica coffee

Xiao et al. Beverage Plant Research 2023, in press
Page 3 of 10

I n P r e s s
among the indicators [25] .Two calculation steps are used to evaluate the weights of each criterion using the CRITIC method.
Step 1 is to normalize the decision matrix using Equation (10).
Step 2 is to calculate the weights for each criterion where, σ j was the standard deviation of each criterion j, and R ij is the linear correlation coefficient between the vectors X i and X j.The comprehensive evaluation value of the effects of leaf anatomical traits on photosynthetic parameters were estimated from Equation (13).
where D ij was the correlation (direct effects, indirect effects, and Pearson's correlation, respectively, respectively) between leaf anatomical trait i and photosynthetic parameter j.All statistical analyses were conducted with IBM SPSS Statistics 24.

Results
The photosynthetic light-response curve (Fig. 2) was well fitted by the modified rectangular hyperbola model (R 2 = 0.945 ± 0.006) described by Ye (2007) [23] .The P n increased with I below the I sat .As I rise above the I sat , P n decreases as I increase, implying a photoinhibition phenomenon.The calculated AQY would be 0.028 ± 0.003 if the Kok effect was ignored.This value was lower than α (0.031 ± 0.003) and the IQY (0.033 ± 0.003).The maximum P n values calculated by the modified rectangular hyperbola model at the I sat of 621.762 ± 37.358 µmol (photon) m -2 •s -1 was 5.131 ± 0.478 µmol (CO 2 ) m -2 •s -1 for coffee leaf.The values of I c and R D were 60.238 ± 11.332 µmol (photon) m -2 •s -1 and 0.930 ± 0.098 µmol (CO 2 ) m -2 •s -1 , respectively.
Cross-sections observed under a light microscope revealed that the coffee leaves were bifacial, with a typical well-differentiated single layer of elongated palisade mesophyll cells on the adaxial side and large spongy mesophyll cells on the abaxial side.The large spongy mesophyll cells at the abaxial side comprised almost 60% of the total blade thickness (Table 1).Of all the tissues, the palisade tissue thickness (PT) showed the highest variation in dimensions (CV = 13.243%).The stomatal density (d) and total stomatal pore area (TSP) were 152.028 and 0.439, respectively (Table 1).

Q i
To identify the combined effect of leaf anatomical traits on photosynthetic parameters, a comprehensive evaluation analysis was conducted.Based on the evaluation values ( ), the ranking of the comprehensive scores of relationships between photosynthetic capacity and leaf anatomical traits was shown in Table 3.For the direct effect, the LT had the highest comprehensive score, followed by LET, d, ST, PT, UET, TSP LCT, and UCT.For the comprehensive score of indirect effect of leaf anatomical traits on photosynthetic capacity, the UET held the highest comprehensive score, and UCT was the second, followed by ST, TSP, PT, d, LCT, LT, and LET.Similar results from the Pearson correlation analyses of comprehensive assessments revealed that the first comprehensive assessment score was d and that the second assessment score was also ST.Then, the ranking order after the third was UET, LT, TSP, PT, LET, LCT, and UCT.

Discussion
Several studies involving gas exchange measurements of coffee leaves have been performed.These studies indicate that the photosynthetic capacity of coffee leaves varies with varying environmental conditions [27−30] .For single coffee leaves, the saturating irradiance is between 300 and 700 µmol m −2 •s −1 , with shade leaves showing a lower value than lower sun leaves [28,31] .This study found that the saturating irradiance of sun leaves was about 620 µmol photons m −2 s −1 , which was consistent with the result of the above mention.Quantum yields revealed the relationship between a given light-dependent product and the number of absorbed photons [32] .Quantum yield is 0 when none of the light energy is used in photo-    I n P r e s s synthesis.Quantum yield is 1 when all the absorbed light is used.Based on the results of YE ( 2007), the quantum yield decreased with I increasing [23] .Thus, the IQY was higher than AQY.In the past, there were several pertinent points about the maximal photosynthetic rates of coffee leaves.For example, Kumar and Tieszen (1980) and Cannell (1985) have pointed out that the maximal photosynthetic rates of sun leaves of coffee are around 7 µmol CO 2 m −2 s −1 and 8.8 µmol CO 2 m −2 s −1 , respectively [27,33] .However, Bote et al. (2018) reported that the maximal photosynthetic rates of sun leaves of coffee were lower than 7 µmol CO 2 m −2 s −1 with a larger N supply [34] .In this study, the maximal photosynthetic rates were lower than the observations of Kumar and Tieszen (1980) and Cannell (1985)  and consistent with the results of Bote et al. (2018).For coffee leaves, light saturation would lead to an excess of electron transport capacity [35] .Therefore, a higher J max was observed under light saturation.In addition, the values of V cmax and J max were larger than the Martins et al. (2014) reported result.The correlation analysis is a great utility method for understanding the relationships between variables.In this study, significant correlations between photosynthetic capacity and leaf anatomical traits were observed, indicating that leaf anatomical traits could affect photosynthetic capacity.Despite the correlation coefficients could provide a good way to understand the relationships between variables, the cause and effect relationships between variables could not be estimated [36] .Path analysis, a method that investigates the causal relationship, gives detailed understanding of positive and negative correlations between traits [36] .The result of the path analysis was shown in figure 3, indicating that there were high-level direct and indirect effects of leaf anatomical traits on photosynthetic capacity.This may explain the importance of leaf structure in photosynthesis.However, the estimated residual effect in this study was high (0.661-0.908), indicating that additional characters which affect photosynthetic capacity are not included in the investigation.For leaf photosynthesis, the overall photosynthetic limitation can be partitioned into different components, such as stomatal, mesophyll, and biochemical limitations [37] .The stomatal and mesophyll limitations in coffee leaves accounted for 0.30 and 0.38, respectively [35] .This is consistent with the results, where about 30% of the variability in photosynthetic capacity was contributed by the characters studied in the path analysis in this study.In order to quantitatively evaluate the effect of each structure of the leaf on photosynthesis, the comprehensive evaluation values calculated by using the formula (13) introduced in the methods.The ordering Pearson's correlation of comprehensive score showed that stomatal density (d) and spongy tissue thickness (ST) were the most important leaf anatomical traits for coffee leaf photosynthesis.For Arabica coffee, leaf photosynthesis is greatly limited by CO 2 diffusion [6,38] .The CO 2 diffusion from the air to the leaf mesophyll is mainly modulated by the stomatal conductance (g s ), which is associated with stomatal density (d) [5,39] .According to the comprehensive score of indirect effect, the stomatal density (d) also had the highest comprehensive score.This indicates that stomatal density (d) is essential for coffee leaf photosynthesis.In fact, higher stomatal density will increase g s by increasing the same total pore area [39] and allow more CO 2 diffusion.This may partly explain why the total stomach pore area (TSP) had higher comprehensive evaluation values in the effect of coffee leaf photosynthesis.Therefore, the result of a comprehensive evaluation revealed that Arabica coffee leaf of stomatal characteristics, such as guard cell width and length, stomatal density, and stomatal pore surface, plays a significant role in leaf photosynthetic capacity.This finding was also consistent with previous work, indicating that stomatal characteristics of the leaf are closely associated with photosynthesis by controlling the water loss and CO 2 uptake [40,41] .
The leaf thickness was strongly positively correlated with the fraction of intercellular air space in the leaves [18] .The result obtained in this study showed that leaf thickness (LT) affects photosynthesis mainly through direct and indirect effects, given that the CO 2 diffusion in the gaseous phase is about three orders of magnitude larger than in the liquid [12] .Thus, the diffusive resistance should not necessarily increase considerably with increasing leaf thickness [42] .The result of direct effects of leaf thickness (LT) on photosynthetic capacity confirmed that increasing leaf thickness may increase g i .Moreover, due to the chloroplasts adhering to the inner surface of exposed mesophyll cell walls, thicker leaves provide space for more chloroplasts per unit leaf area [5] and possess more chloroplast surface area, which results in a substantial increase in carbon dioxide absorption because of a larger surface area for diffusion [42] .The direct effects of the path analysis indicate that greater leaf thickness would facilitate CO 2 diffusion in the mesophyll (g i ).Furthermore, thicker leaves may have higher total contents of Rubisco and leaf N per unit area [43] , so it advantageous to realize high photosynthetic rates in highlight environments [44] .In this study, leaf thickness (LT) affects photosynthetic indexes that were IQY, α, AQY, R d , g i , and I sat , and, consequently, photo- I n P r e s s synthesis rate.Therefore, the result of a comprehensive evaluation revealed that leaf thickness (LT) was the most important role in the direct effect on leaf photosynthetic capacity.The structure of the leaf cuticle consists of epicuticular wax, birefringent wax embedded in cutin, cutin matrix, and pectic substance [45] .The epicuticular wax can decrease light by scattering [34] .Similar results have been found in this study where an increase in upper cuticle thickness (UCT) could result in a decrease in IQY, AQY, and α, which were related to light absorption in photosynthesis.The most cutin matrix in the cuticle is heterogeneous, and the areas of embedded waxes exhibit birefringence in polarized light [45] .This heterogeneous structure of the leaf cuticle may affect epidermal focusing, allowing actinic light to penetrate deeper into the mesophyll [34] .These findings most likely explained the significant relationship between cuticle thickness and quantum yield of photosynthesis, such as IQY, AQY, and α.In addition, the cuticle thickness would increase resistance to CO 2 diffusion [46] , and the negative relationships between UCT and LCT and photosynthetic parameters were also observed in this study.
The epidermal cell usually has a large transparent central vacuole and lacks mature chloroplasts, thus the light can directly pass from these cells into the mesophyll [34] .However, these epidermal cells have various shapes, such as planoconvex, spherical, conical, or some other convex shape [47] .The curved outermost epidermal cell wall can affect the focal point within the leaf [48] .However, the result of Brodersen and Vogelmann (2007) showed that there was no relation between the convexity of epidermal cells and the absorptance of diffuse light [49] .Although two contradictory hypotheses have been proposed concerning the impacts of epidermal cells on photosynthesis, the significant relationship between UET and photosynthetic capacity parameters in this study examined epidermal lens effects on many photosynthetic parameters and elucidated epidermal cells property plays an important role in photosynthetic capacity.
The columnar palisade tissue minimizes light scattering when the incident light collimated with the columnar palisade [50] .As a result, palisade tissue allows large amounts of light to penetrate the chloroplasts within the leaf [50] .The present study has revealed a remarkable relation between quantum yields and PT.However, the results of this study indicated that PT increase would directly decrease the IQY, AQY, and J max .This circumstance is very difficult to comprehend.Generally, palisade tissue has much more chloroplasts compared with other leaf tissue [51] .Thus, the photosynthetic activity of the palisade tissue would be higher.Indeed, the shape of cells in the leaves is strongly associated with photosynthetic performance by affecting the movement and distribution of chloroplasts [52] .There was a great variety of dimensions of the palisade tissue thickness (PT) in this study.Coordinated regulation of leaf cell shape and chloroplast motion according to light conditions is essential for efficient leaf photosynthesis [52] .This work could partially explain the result of a negative direct effect of PT on photosynthetic parameters.
The concentration of CO 2 in the mesophyll is affected by air temperature, and each degree rise in temperature above 24°C results in a 20 ppm rise in CO 2 concentration [53] .However, Khairi and Hall's (1976) studies of citrus photosynthesis revealed that mesophyll conductance to CO 2 decreased as temperature increased from 22 °C to around 40 °C (Khairi and   Hall 1976). [54].In the analysis of the effects of leaf anatomical traits on photosynthetic capacity, the spongy tissue thickness (ST) had larger indirect effects and Pearson's correlation comprehensive score, suggesting that spongy tissue thickness (ST) plays an important role in the photosynthetic process.Furthermore, Martins et al. (2014) demonstrated that mesophyll thickness is related to water flux and gas exchange per leaf area [35] , which is consistent with the findings of this study.The spongy tissue thickness (ST) had direct effects on WUE and VPD.On the other hand, the irregularly shaped spongy cells alternating with air spaces increase light absorption due to the increased optical path length through a leaf [12] , and this absorbed light will be advanced for photosynthesis [33,49] .Therefore, the high direct spongy tissue thickness (ST) effects upon I sat were also observed in this study.In fact, the increase in spongy tissue thickness (ST) would increase the surface area of chloroplasts exposed to intercellular airspace, which was associated with CO 2 transfer conductance [55] .Although there was little relationship between spongy tissue thickness (ST) and photosynthetic parameters, high direct spongy tissue thickness (ST) effects upon I sat , WUE, and VPD were observed.
Stomatal density (SD) is a factor associated with the photosynthetic rates of the leaf [56] .The CO 2 diffusion conductance is influenced by the density and arrangement of the stomata [12] .For the coffee leaf, the net carbon assimilation rate (A) is greatly governed by the diffusive conductance of CO 2 , thus the photosynthesis of the coffee leaf was largely limited by a diffusive factor [57] .Stomata have traditionally been thought to play an important role in controlling CO 2 diffusion [57] and were a key parameter in C. arabica photosynthesis [30] .Among the various stomatal characters, stomatal density plays a major role in gas exchange and photosynthesis [58] .In this study, it is shown that stomatal density had the highest comprehensive scores compared with other anatomical traits on indirect effects and Pearson's correlation comprehensive evaluation.In addition, structural modification of the stomata is one of the key factors affecting stomatal conductance [32] .The result of this study also showed that the guard cell length (l) significantly correlated positively with C i .
For light absorption, although the result of Martins et al.  (2014) showed that the guard cell length did not change in response to light treatments [35] , there was a significant correlation with IQY and AQY.The coordination between increased initiation of stomata cells and expansion of epidermal cells was observed in coffee leaves, implying that there was an optimization of the trade-off between transpiration costs and CO 2 assimilation [35] .In addition, the size and density of the stomata largely affect the diffusive conductance of the leaves to CO 2 , and the balance between the amount of CO 2 required for photosynthesis and the level of water availability is also determined by the stomatal traits [39] .Thus, lower epidermal stomatal density had a significant correlation with photosynthetic parameters or direct effects and indirect effects on photosynthetic parameters.Moreover, these traits also had larger comprehensive scores in a comprehensive evaluation, which means that stomatal density was a key parameter in the regulation of photosynthesis of C. arabica.I n P r e s s thetic parameters.This result provided holistic relationships between leaf anatomical traits and photosynthetic parameters in coffee.It is tempting to suggest that these relationships might, to some extent, explain the photosynthetic behavior of coffee (C.arabica L.).The relationships between leaf anatomical traits and photosynthetic parameters emphasize the importance of leaf anatomy in determining photosynthesis at the tissue level.Moreover, these correlations between leaf anatomical traits and photosynthetic parameters were useful to explore the impact of different anatomical features on photosynthetic efficiency and accordingly provide some information to design leaf anatomy for enhanced photosynthetic efficiency.

Conclusion
To a certain extent, path analysis and comprehensive evaluation could indicate an interrelationship between photosynthetic capacity and leaf anatomical traits.However, the causal interrelationship between photosynthetic parameters and leaf anatomical traits is less clear.In addition, This study only analyzed the effect of anatomical structure on photosynthesis for one Catimor type of Arabica coffee.Effects on photosynthesis are typically genotype dependent, with some varieties being more tolerant to partial shading than others.Future studies should compare different Arabica coffee varieties for their anatomical traits and photosynthetic response to suitability in agroforestry-based production systems.Despite this, the result of the study also suggests that thinner leaves and higher stomatal densities might support a higher photosynthetic capacity.This result may be useful information for coffee breeding programs and plantation management.
Photosynthetic capacity parameters Photosynthetic capacity parameters
Many leaf anatomical traits significantly affect photosyn-The Influence of Leaf Anatomical Traits on Photosynthesis in Catimor type Arabica coffee Page 8 of 10 Xiao et al.Beverage Plant Research 2023, in press

Table 1 .
The components of anatomical tissues in Arabica coffee leaf

Table 2 .
The correlation between leaf anatomical traits and photosynthetic parameters LT -leaf thickness, UCT -upper cuticle thickness, UET -upper epidermal thickness, PT -palisade tissue thickness, ST -spongy tissue thickness, LET -lower epidermal thickness, LCT -lower cuticle thickness, w -guard cell width, l -guard cell length, d -stomatal density, TSP -total stomach pore area.Ns = nonsignificant.*, ** represent significant; highly significant differences at 5 % and 1 % probability levels, respectively.The Influence of Leaf Anatomical Traits on Photosynthesis in Catimor type Arabica coffee Xiao et al.Beverage Plant Research 2023, in press Page 5 of 10 I n P r e s s

Table 3 .
Comprehensive ordering of the effects of leaf anatomical traits on photosynthetic capacity with CRITIC method.LT -leaf thickness, UCT -upper cuticle thickness, UET -upper epidermal thickness, PT -palisade tissue thickness, ST -spongy tissue thickness, LET -lower epidermal thickness, LCT -lower cuticle thickness, w -guard cell width, l -guard cell length, d -stomatal density, TSP -total stomach pore area.The Influence of Leaf Anatomical Traits on Photosynthesis in Catimor type Arabica coffee Xiao et al.Beverage Plant Research 2023, in press Page 7 of 10