Spatial distribution characteristics of stomata at the areole level in Michelia cavaleriei var. platypetala (Magnoliaceae)

A Nonlinear Fitting Method Provides Strong Support for Geometric Series of Stomatal Area in 12 Magnoliaceae Species

Stomatal pore area and density determine the capacity for gas exchange between the leaf interior and the atmosphere. Stomatal area is given by the profile formed by two guard cells, and the cumulative stomatal area characterizes the area of leaf surface occupied by stomata. The areas of all stomata captured in a micrograph are sorted in ascending order to form a sequence, which is referred to as a sequence of stomatal area here. In total, 360 leaves of 12 Magnoliaceae species with 30 leaves for each species were sampled. For each leaf, two 662 μm × 444 μm fields of view (micrographs) of stomata were captured on the right leaf width axis. In each micrograph, the length and width of each stoma were measured, and the area of the stoma was determined using the product of stomatal length and width multiplied by a proportionality coefficient. Stomatal area sequences of Magnoliaceae in the constant field of view were found to follow a geometric series (GS). Prior studies estimated the common ratio of the GS as the mean of the quotients of any two adjacent terms, and estimated the first term as the mean of the first terms (i.e., the smallest stomatal area) represented by the quotient of each term and the estimated common ratio to a power of the order of the term minus 1, which is referred to as Method-1. However, it produced large prediction errors for some stomatal area sequences. In the present study, the nonlinear regression was used to fit the stomatal area sequences using the common ratio and the first term as two model parameters (Method-2). We compared the two methods using the mean absolute percent error (MAPE, ≤5% considered as a good fit) values of the 720 stomatal micrographs from the 12 Magnoliaceae species. The goodness of fit of Method-2 was better than that of Method-1 (52.4% MAPE values were ≤5% for Method-1 and 99.6% for Method-2). There were significant variations in the estimated common ratios, as well as the estimated first terms and the MAPE values across the 12 Magnoliaceae species, but overall, the interspecific differences in the MAPE values were small. We conclude that the GS hypothesis for the stomatal area sequences of the 12 Magnoliaceae species was further strengthened by the new method. This method further provides a valuable approach for the calculation of total stomatal area per unit leaf area.

Geometric series exists in nature: Evidence from sorted area sequences of floral parts and leaves

The concept of a geometric series (GS) plays an important role in mathematics. However, it has been neglected in describing biological size series. Herein, we show that a GS describes the nonreproductive (perianth) parts of the flowers of four Magnoliaceae species and two Rosaceae species and the leaves of 60 Alangium chinense and 60 Shibataea chinensis shoots. The sorted areas of floral parts and leaves formed a sequence that was fitted by a GS with the mean of the quotients of two adjacent members in the sequence as the common ratio of a GS. The mean absolute percent error (MAPE) was used to measure the goodness of fit of each GS. Over 99.7% of the MAPE values (371 out of the 372 tested flowers) were less than 10%, and over 97.8% of the MAPE values were less than 5%. Likewise, over 77.5% of the MAPE values (93 out of the 120 tested shoots) were less than 10%, and over 35% of the MAPE values were less than 5%. These analyses provide empirical evidence that the GS exists in nature, and confirm the usefulness of a classical algebraic formula for the study of plant developmental biology.

Using the Montgomery-Koyama-Smith equation to calculate the stomatal area per unit lamina area for 12 Magnoliaceae species

BACKGROUND AND AIMS

The Montgomery-Koyama-Smith (MKS) equation predicts that total leaf area per shoot is proportional to the product of the sum of individual leaf widths and maximum individual leaf length, which has been validated for some herbaceous and woody plants. The equation is also predicted to be valid in describing the relationship between the total stomatal area per micrograph (AT) and the product of the sum of individual stomatal widths (denoted as LKS) and maximum individual stomatal length (denoted by WKS) in any particular micrograph.

METHODS

To test the validity of the MKS equation, 69 931 stomata (from 720 stomatal micrographs from 12 Magnoliaceae species) were examined. The area of each stoma was calculated using empirical measurements of stomatal length and width multiplied by a constant. Six equations describing the relationships among AT, LKS and WKS were compared. The root mean square (RMSE) and the Akaike information criterion (AIC) were used to measure the goodness of fit and the trade-off between the goodness of fit and the structural complexity of each model, respectively.

KEY RESULTS

Analyses supported the validity of the MKS equation and the power-law equation AT ∝ (LKSWKS)α, where α is a scaling exponent. The estimated values of α at the species level and for the pooled data were all statistically smaller than unity, which did not support the hypothesis that AT ∝ LKSWKS. The power-law equation had smaller RMSE and AIC values than the MKS equation for the data from the 12 individual species and the pooled data.

CONCLUSIONS

These results indicate that AT tends to scale allometrically with LKSWKS, and that increases in AT do not keep pace with increases in LKSWKS. In addition, using LKSWKS is better than using only one of the two variables to calculate AT.

Significant correlation between leaf vein length per unit area and stomatal density: evidence from Red Tip and Chinese photinias

The vascular veins in photosynthetic leaves play an important role in transporting water and sugars throughout the plant body, and their venation pattern and vein density determine the hydraulic efficiency of the leaf. Likewise, stomatal density (SD) can influence photosynthetic gas exchange. However, the correlation between leaf vein density and SD is seldom reported. Herein, we examined 16 leaves from the hybrid Photinia × fraseri and 16 leaves from one of its parents, P. serratifolia, to explore the correlation between leaf vein density and SD. For each leaf, equidistant lamina quadrats were excised along two longitudinal transects (one along the midrib and another along the leaf margin). For each quadrat, micrographs of 1.2 mm × 0.9 mm stomatal imprints, and 2.51 mm × 1.88 mm micrographs of leaf veins were used to measure total vein area per leaf unit area (VAA) and total vein length per unit area (VLA), as indicators of leaf vein density, to determine the correlation between SD and leaf vein density. For each taxon, there was no significant correlation between SD and VAA, but there was a significant correlation between SD and VLA. The data indicate that SD is not positively correlated with VAA but positively correlated with VLA for both the hybrid and the parent species. This study indicates that future work should focus on the relationships between SD and total vein length per unit area rather than on total leaf vein area per unit area within and across species.

Vertical variation in leaf functional traits of Parashorea chinensis with different canopy layers

Introduction

Canopy species need to shift their ecological adaptation to improve light and water resources utilization, and the study of intraspecific variations in plant leaf functional traits based at individual scale is of great significance for evaluating plant adaptability to climate change.

Methods

In this study, we evaluate how leaf functional traits of giant trees relate to spatial niche specialization along a vertical gradient. We sampled the tropical flagship species of Parashorea chinensis around 60 meters tall and divided their crowns into three vertical layers. Fourteen key leaf functional traits including leaf morphology, photosynthetic, hydraulic and chemical physiology were measured at each canopy layer to investigate the intraspecific variation of leaf traits and the interrelationships between different functional traits. Additionally, due to the potential impact of different measurement methods (in-situ and ex-situ branch) on photosynthetic physiological parameters, we also compared the effects of these two gas exchange measurements.

Results and discussion

In-situ measurements revealed that most leaf functional traits of individual-to-individual P. chinensis varied significantly at different canopy heights. Leaf hydraulic traits such as midday leaf water potential (MWP) and leaf osmotic potential (OP) were insignificantly correlated with leaf photosynthetic physiological traits such as maximal net assimilation rate per mass (A mass). In addition, great discrepancies were found between in-situ and ex-situ measurements of photosynthetic parameters. The ex-situ measurements caused a decrease by 53.63%, 27.86%, and 38.05% in A mass, and a decrease of 50.00%, 19.21%, and 27.90% in light saturation point compared to the in-situ measurements. These findings provided insights into our understanding of the response mechanisms of P. chinensis to micro-habitat in Xishuangbanna tropical seasonal rainforests and the fine scale adaption of different resultant of decoupled traits, which have implications for understanding ecological adaption strategies of P. chinensis under environmental changes.

An Inverse Scaling Relationship between Stomatal Density and Mean Nearest Neighbor Distance: Evidence from a Photinia Hybrid and One of Its Parents

Stomata are involved in transpiration and CO2 uptake by mediating gas exchange between internal plant tissues and the atmosphere. The capacity for gas exchange depends on stomatal density (SD), stomatal size, and pore dimensions. Most published work on stomatal quantification has assumed that stomatal distribution and stomatal density are spatially homogeneous across the leaf, but this assumption has been seldom tested. We selected 32 leaves from a Photinia hybrid, Photinia × fraseri 'Red Robin', and one of its parents, P. serratifolia. For each leaf, the leaf surface was divided into three or four equidistant layers along the apical-basal axis, and, in each layer, two positions, one closer to the midrib and the other closer to the leaf margin, were further selected. We calculated SD and mean nearest neighbor distance (MNND) for each lamina section and tested the scaling relationship between SD and MNND of the sampled stomatal centers using reduced major axis protocols. In addition, we calculated the stomatal aggregation index (SAI) for each lamina section to examine the spatial arrangement of stomata at the given size of field of view of 1.2 mm × 0.9 mm. We observed that SD decreased from the lamina apex towards the base for central lamina areas but varied little at leaf margins. An inverse scaling relationship between SD and MNND was observed for both species. This relationship could be used for SD estimation using the rapidly estimated trait, MNND. SAI did not vary significantly throughout leaf lamina, and the numerical values of SAI for all fields of view were greater than one, which indicates significant spatial repulsion between stomata. The study suggests that SD varies across leaf lamina to fine-tune plant water use and maximize carbon gain. However, spatial structures of stomata from different lamina sections exhibit similar patterns (i.e., spatial inhibition between stomata at small scales), probably due to hierarchical leaf vein patterns.

Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.)

Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower-order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Limiting resource and leaf functional traits jointly determine distribution patterns of leaf intrinsic water use efficiency along aridity gradients

Intrinsic water use efficiency (iWUE) is a critical eco-physiological function allowing plants to adapt to water- and nutrient-limited habitats in arid and semi-arid regions. However, the distribution of iWUE in coexisting species along aridity gradients and its controlling factors are unknown. We established two transects along an aridity gradient in the grasslands of Losses Plateau (LP) and Inner Mongolia Plateau (MP) to elucidate the patterns and underlying mechanisms of iWUE distribution in coexisting species along aridity gradient. We determined leaf carbon (δ¹³C) and oxygen (δ¹⁸O) stable isotopes, functional traits related to carbon fixation, and limiting resources. Bulk leaf δ¹³C and δ¹⁸O were used as proxies for time-integrated iWUE and stomatal conductance (gs) during the growing season. Our results showed that variability in iWUE within transect was primarily controlled by species, sampling sites and an interactive effect between species and sampling sites. Mean values of iWUE (iWUE_Mean) increased and coefficient of variation (CV) in iWUE (iWUE_CV) decreased with an increase in aridity, demonstrating that increases in aridity lead to conservative and convergent water use strategies. Patterns of iWUE_Mean and iWUE_CV were controlled primarily by the ratio of soil organic carbon to total nitrogen in LP and soil moisture in MP. This revealed that the most limited resource drove the distribution patterns of iWUE along aridity gradients. Interspecific variation in iWUE within transect was positively correlated with Δ¹⁸O, indicating that interspecific variation in iWUE was primarily regulated by gs. Furthermore, relationship between iWUE and multi-dimensional functional trait spectrum indicated that species evolved species-specific strategies to adapt to a harsh habitat by partitioning limiting resources. Overall, these findings highlighted the interactive effects of limiting resources and leaf functional traits on plant adaptation strategies for iWUE, and emphasized the importance of considering biological processes in dissecting the underlying mechanisms of plant adaptation strategies at large regional scales.

Scaling relationships of leaf vein and areole traits versus leaf size for nine Magnoliaceae species differing in venation density

PREMISE

Across species, main leaf vein density scales inversely with leaf area (A). Yet, minor vein density manifests no clear relationship with respect to A, despite having the potential to provide important insights into the trade-off among the investments in leaf mechanical support, hydraulics, and light interception.

METHODS

To examine this phenomenon, the leaves of nine Magnoliaceae leaves were sampled, and the scaling relationships among A and midrib length (ML), total vein length (TVL), total vein area (TVA), total areole area (TAA), and mean areole area (MAA) were determined. The scaling relationships between MAA and areole density (the number of areoles per unit leaf area) and between MAA and A were also analyzed.

RESULTS

For five of the nine species, A was proportional to ML² . For eight of the nine species, TVL and TVA were both proportional to A. The numerical values of the scaling exponents for TAA vs. A were between 1.0 and 1.07 for eight species; i.e., as expected, TAA was isometrically proportional to A. There was no correlation between MAA and A, but MAA scaled inversely with respect to areole density for each species.

CONCLUSIONS

The correlation between midrib "density" (i.e., ML/A) and A, and the lack of correlation between total leaf vein density and A result from the A ∝$\propto $ ML² scaling relationship and the proportional relationship between TVL and A, respectively. Leaves with the same size can have widely varying MAA. Thus, leaf size itself does not directly constrain leaf hydraulic efficiency and redundancy.

Scaling relationships of leaf vein and areole traits versus leaf size for nine Magnoliaceae species differing in venation density

PREMISE

Across species, main leaf vein density scales inversely with leaf area (A). Yet, minor vein density manifests no clear relationship with respect to A, despite having the potential to provide important insights into the trade-off among the investments in leaf mechanical support, hydraulics, and light interception.

METHODS

To examine this phenomenon, the leaves of nine Magnoliaceae leaves were sampled, and the scaling relationships among A and midrib length (ML), total vein length (TVL), total vein area (TVA), total areole area (TAA), and mean areole area (MAA) were determined. The scaling relationships between MAA and areole density (the number of areoles per unit leaf area) and between MAA and A were also analyzed.

RESULTS

For five of the nine species, A was proportional to ML² . For eight of the nine species, TVL and TVA were both proportional to A. The numerical values of the scaling exponents for TAA vs. A were between 1.0 and 1.07 for eight species; i.e., as expected, TAA was isometrically proportional to A. There was no correlation between MAA and A, but MAA scaled inversely with respect to areole density for each species.

CONCLUSIONS

The correlation between midrib "density" (i.e., ML/A) and A, and the lack of correlation between total leaf vein density and A result from the A ∝$\propto $ ML² scaling relationship and the proportional relationship between TVL and A, respectively. Leaves with the same size can have widely varying MAA. Thus, leaf size itself does not directly constrain leaf hydraulic efficiency and redundancy.