Structural organization of the spongy mesophyll

Three-dimensional reconstruction and multiomics analysis reveal a unique pattern of embryogenesis in Ginkgo biloba

Ginkgo (Ginkgo biloba L.) is one of the earliest extant species in seed plant phylogeny. Embryo development patterns can provide fundamental evidence for the origin, evolution, and adaptation of seeds. However, the architectural and morphological dynamics during embryogenesis in G. biloba remain elusive. Herein, we obtained over 2,200 visual slices from 3 stages of embryo development using micro-computed tomography imaging with improved staining methods. Based on 3-dimensional (3D) spatiotemporal pattern analysis, we found that a shoot apical meristem with 7 highly differentiated leaf primordia, including apical and axillary leaf buds, is present in mature Ginkgo embryos. 3D rendering from the front, top, and side views showed 2 separate transport systems of tracheids located in the hypocotyl and cotyledon, representing a unique pattern of embryogenesis. Furthermore, the morphological dynamic analysis of secretory cavities indicated their strong association with cotyledons during development. In addition, we identified genes GbLBD25a (lateral organ boundaries domain 25a), GbCESA2a (cellulose synthase 2a), GbMYB74c (myeloblastosis 74c), GbPIN2 (PIN-FORMED 2) associated with vascular development regulation, and GbWRKY1 (WRKYGOK 1), GbbHLH12a (basic helix-loop-helix 12a), and GbJAZ4 (jasmonate zim-domain 4) potentially involved in the formation of secretory cavities. Moreover, we found that flavonoid accumulation in mature embryos could enhance postgerminative growth and seedling establishment in harsh environments. Our 3D spatial reconstruction technique combined with multiomics analysis opens avenues for investigating developmental architecture and molecular mechanisms during embryogenesis and lays the foundation for evolutionary studies of embryo development and maturation.

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

Beauty is more than epidermis deep: How cell division and expansion sculpt the leaf spongy mesophyll

As the main location of photosynthesis, leaf mesophyll cells are one of the most abundant and essential cell types on earth. Forming the bulk of the internal tissues of the leaf, their size, shape, and patterns of interconnectivity define the internal structure and surface area of the leaf, which in turn determines the efficiency of light capture and carbon fixation. Understanding how these cellular traits are controlled and translated into tissue- and organ-scale traits, and how they influence photosynthetic performance will be key to our ability to improve crop plants in the face of a changing climate. In contrast to the extensive literature on the anatomical and physiological aspects of mesophyll function, our understanding of the cell-level morphogenetic processes underpinning mesophyll cell growth and differentiation is scant. In this review, we focus on how cell division, expansion, and separation are coordinated to create the intricate architecture of the spongy mesophyll.

Leaf morpho-anatomical adjustments in a Quercus pubescens forest after 10 years of partial rain exclusion in the field

In the Mediterranean region, a reduction of annual precipitation and a longer and drier summer season are expected with climate change by the end of the century, eventually endangering forest survival. To cope with such rapid changes, trees may modulate their morpho-anatomical and physiological traits. In the present study, we focused on the variation in leaf gas exchange and different leaf morpho-anatomical functional traits of Quercus pubescens Willd. in summer using a long-term drought experiment in natura consisting of a dynamic rainfall exclusion system where trees have been submitted to amplified drought (AD) (~-30% of annual precipitation) since April 2012 and compared them with trees under natural drought (ND) in a Mediterranean forest. During the study, we analyzed net CO2 assimilation (An), stomatal conductance (gs), transpiration (E), water-use efficiency (WUE), stomatal size and density, density of glandular trichomes and non-glandular trichomes, thickness of the different leaf tissues, specific leaf area and leaf surface. Under AD, tree functioning was slightly impacted, since only An exhibited a 49% drop, while gs, E and WUE remained stable. The decrease in An under AD was regulated by concomitant lower stomatal density and reduced leaf thickness. Trees under AD also featured leaves with a higher non-glandular trichome density and a lower glandular trichome density compared with ND, which simultaneously limits transpiration and production costs. This study points out that Q. pubescens exhibits adjustments of leaf morpho-anatomical traits which can help trees to acclimate to AD scenarios as those expected in the future in the Mediterranean region.

Plant microphenotype: from innovative imaging to computational analysis

The microphenotype plays a key role in bridging the gap between the genotype and the complex macro phenotype. In this article, we review the advances in data acquisition and the intelligent analysis of plant microphenotyping and present applications of microphenotyping in plant science over the past two decades. We then point out several challenges in this field and suggest that cross-scale image acquisition strategies, powerful artificial intelligence algorithms, advanced genetic analysis, and computational phenotyping need to be established and performed to better understand interactions among genotype, environment, and management. Microphenotyping has entered the era of Microphenotyping 3.0 and will largely advance functional genomics and plant science.

Integrated Ca, Mg, Cu, and Zn supply upregulates leaf anatomy and metabolic adjustments in Eucalyptus seedlings

Adaptive responses to abiotic stresses such as soil acidity in Eucalyptus-the most widely planted broad-leaf forest genus globally-are poorly understood. This is particularly evident in physiological and anatomical disorders that inhibit plant development and wood quality. We aimed to explore how the supply of Ca and Mg through liming (lime), combined with Cu and Zn fertilization (CZF), influences physiological and anatomical responses during Eucalyptus grandis seedlings growth in tropical acid soil. Therefore, related parameters of leaf area and leaf anatomy, stomatal size, leaf gas exchange, antioxidant system, nutrient partitioning, and biomass allocation responses were monitored. Liming alone in Eucalyptus increased specific leaf area, stomatal density on the abaxial leaf surface, and Ca and Mg content. Also, Eucalyptus exposed only to CZF increased Cu and Zn content. Lime and CZF increased leaf blade and adaxial epidermal thickness, and improved the structural organization of the spongy mesophyll, promoting increased net CO2 assimilation, and stomatal conductance. Fertilization with Ca, Mg, Cu, and Zn positively affects plant nutrition, light utilization, photosynthetic rate, and antioxidant performance, improving growth. Our results indicate that lime and CZF induce adaptive responses in the physiological and anatomical adjustments of Eucalyptus plantation, thereby promoting biomass accumulation.

A rapid freezing method to determine tissue layer thickness in drought-stressed leaves

Plants have been affected by water stress ever since they settled on dry land. In severe and persisting drought, plant leaves are wilting. However, a documentation at the anatomical level of the minute changes that occur before wilting is challenging. On the other hand, understanding the anatomical alteration in plant leaves with respect to water stress provides a stronger basis to study molecular and submolecular processes through which plants enhance drought tolerance. In this work, we applied an affordable method to visualise mesophyll layers of Arabidopsis thaliana cell lines without preparation steps that would alter the volume of the cells. We rapidly plunge-froze the leaves in liquid nitrogen, cut them while in the N2 bath, and immediately imaged the mesophyll cross sections in a scanning electron microscope. We applied a reduction of watering from 60 to 40 to 20 mL per day and investigated two time points, 7 and 12 days, respectively. Interestingly, the overall thickness of leaves increased in water stress conditions. Our results showed that the palisade and spongy layers behaved differently under varying watering regimes. Moreover, the results showed that this method can be used to image leaf sections after drought stress without the risk of artefacts or swelling caused by contact to liquids as during chemical fixation.

Intraspecific variations of leaf hydraulic, economic, and anatomical traits in Cinnamomum camphora along an urban-rural gradient

Urbanization brings numerous benefits to residents, but it also introduces complex, variable, and heterogeneous habitat conditions to urban plants, resulting in an arid and hot urban environment that decreases tree growth and the ecological service capacity of trees. In this study, we evaluated leaf hydraulic, economic, and anatomical traits and their covariations of Cinnamomum camphora along an urban-rural gradient in Hefei, Eastern China. We found that Cinnamomum camphora in urban adopted a conservative hydraulic strategy with low leaf turgor loss point (Tlp), leaf hydraulic conductance (Kleaf), and leaf water potential resulting in 50 % loss of hydraulic conductance (P50), as well as a quick investment-return economic strategy with low unit leaf dry matter content (LMA) and high leaf nitrogen content (Leaf N). P50, Kleaf and LMA were significantly positively correlated with the urban-rural gradient (PC1urban-rural gradient), while Leaf N exhibited a negative correlation with it. The results showed a trade-off between intraspecific safety and efficiency in leaf hydraulic traits along the urban-rural gradient and an intraspecific coordinated variation in leaf hydraulic and economic traits. In addition, based on the analysis of a trait coordination network, it was revealed that leaf mesophyll and stomata were key structures for trait adjustment and coordination. Furthermore, our findings offer a significant theoretical underpinning for the effective management of landscape trees and the strategic planning of urban tree species.

Conserved cellular patterning in the mesophyll of rice leaves

The mesophyll cells of grass leaves, such as rice, are traditionally viewed as displaying a relatively uniform pattern, in contrast to the clear distinctions of palisade and spongy layers in typical eudicot leaves. This quantitative analysis of mesophyll cell size and shape in rice leaves reveals that there is an inherent pattern in which cells in the middle layer of the mesophyll are larger and less circular and have a distinct orientation of their long axis compared to mesophyll cells in other layers. Moreover, this pattern was observed in a range of rice cultivars and species. The significance of this pattern with relation to potential photosynthetic function and the implication of the widespread use of middle layer mesophyll cells as typical of the rice leaf have been investigated and discussed.

Creating a virtual leaf

When microscopy meets modelling the exciting concept of a 'virtual leaf' is born. The goal of a 'virtual leaf' is to capture complex physiology in a virtual environment, resulting in the capacity to run experiments computationally. One example of a 'virtual leaf' application is capturing 3D anatomy from volume microscopy data and estimating where water evaporates in the leaf and the proportions of apoplastic, symplastic and gas phase water transport. The same 3D anatomy could then be used to improve established 3D reaction-diffusion models, providing a better understanding of the transport of CO₂ across the stomata, through the airspace and across the mesophyll cell wall. This viewpoint discusses recent progress that has been made in transitioning from a bulk leaf approach to a 3D understanding of leaf physiology, in particular, the movement of CO₂ and H₂O within the leaf.

Hydraulic differences between flowers and leaves are driven primarily by pressure-volume traits and water loss

Flowers are critical for successful reproduction and have been a major axis of diversification among angiosperms. As the frequency and severity of droughts are increasing globally, maintaining water balance of flowers is crucial for food security and other ecosystem services that rely on flowering. Yet remarkably little is known about the hydraulic strategies of flowers. We characterized hydraulic strategies of leaves and flowers of ten species by combining anatomical observations using light and scanning electron microscopy with measurements of hydraulic physiology (minimum diffusive conductance (g _min) and pressure-volume (PV) curves parameters). We predicted that flowers would exhibit higher g _min and higher hydraulic capacitance than leaves, which would be associated with differences in intervessel pit traits because of their different hydraulic strategies. We found that, compared to leaves, flowers exhibited: 1) higher g _min, which was associated with higher hydraulic capacitance (C _T); 2) lower variation in intervessel pit traits and differences in pit membrane area and pit aperture shape; and 3) independent coordination between intervessel pit traits and other anatomical and physiological traits; 4) independent evolution of most traits in flowers and leaves, resulting in 5) large differences in the regions of multivariate trait space occupied by flowers and leaves. Furthermore, across organs intervessel pit trait variation was orthogonal to variation in other anatomical and physiological traits, suggesting that pit traits represent an independent axis of variation that have as yet been unquantified in flowers. These results suggest that flowers, employ a drought-avoidant strategy of maintaining high capacitance that compensates for their higher g _min to prevent excessive declines in water potentials. This drought-avoidant strategy may have relaxed selection on intervessel pit traits and allowed them to vary independently from other anatomical and physiological traits. Furthermore, the independent evolution of floral and foliar anatomical and physiological traits highlights their modular development despite being borne from the same apical meristem.

Leaf anatomy does not explain the large variability of mesophyll conductance across C3 crop species

Increasing mesophyll conductance of CO₂ (g_m ) is a strategy to improve photosynthesis in C₃ crops. However, the relative importance of different anatomical traits in determining g_m in crops is unclear. Mesophyll conductance measurements were performed on 10 crops using the online carbon isotope discrimination method and the 'variable J' method in parallel. The influences of crucial leaf anatomical traits on g_m were evaluated using a one-dimensional anatomical CO₂ diffusion model. The g_m values measured using two independent methods were compatible, although significant differences were observed in their absolute values. Quantitative analysis showed that cell wall thickness and chloroplast stroma thickness are the most important elements along the diffusion pathway. Unexpectedly, the large variability of g_m across crops was not associated with any investigated leaf anatomical traits except chloroplast thickness. The g_m values estimated using the anatomical model differed remarkably from the values measured in vivo in most species. However, when the species-specific effective porosity of the cell wall and the species-specific facilitation effect of CO₂ diffusion across the membrane and chloroplast stoma were taken into account, the model could output g_m values very similar to those measured in vivo. These results indicate that g_m variation across crops is probably also driven by the effective porosity of the cell wall and effects of facilitation of CO₂ transport across the membrane and chloroplast stroma in addition to the thicknesses of the elements.

Analyzing anatomy over three dimensions unpacks the differences in mesophyll diffusive area between sun and shade Vitis vinifera leaves

Leaves grown at different light intensities exhibit considerable differences in physiology, morphology and anatomy. Because plant leaves develop over three dimensions, analyses of the leaf structure should account for differences in lengths, surfaces, as well as volumes. In this manuscript, we set out to disentangle the mesophyll surface area available for diffusion per leaf area (S m,LA) into underlying one-, two- and three-dimensional components. This allowed us to estimate the contribution of each component to S m,LA, a whole-leaf trait known to link structure and function. We introduce the novel concept of a 'stomatal vaporshed,' i.e. the intercellular airspace unit most closely connected to a single stoma, and use it to describe the stomata-to-diffusive-surface pathway. To illustrate our new theoretical framework, we grew two cultivars of Vitis vinifera L. under high and low light, imaged 3D leaf anatomy using microcomputed tomography (microCT) and measured leaf gas exchange. Leaves grown under high light were less porous and thicker. Our analysis showed that these two traits and the lower S m per mesophyll cell volume (S m,Vcl) in sun leaves could almost completely explain the difference in S m,LA. Further, the studied cultivars exhibited different responses in carbon assimilation per photosynthesizing cell volume (A Vcl). While Cabernet Sauvignon maintained A Vcl constant between sun and shade leaves, it was lower in Blaufränkisch sun leaves. This difference may be related to genotype-specific strategies in building the stomata-to-diffusive-surface pathway.

Localized growth and remodelling drives spongy mesophyll morphogenesis

The spongy mesophyll is a complex, porous tissue found in plant leaves that enables carbon capture and provides mechanical stability. Unlike many other biological tissues, which remain confluent throughout development, the spongy mesophyll must develop from an initially confluent tissue into a tortuous network of cells with a large proportion of intercellular airspace. How the airspace in the spongy mesophyll develops while the tissue remains mechanically stable is unknown. Here, we use computer simulations of deformable polygons to develop a purely mechanical model for the development of the spongy mesophyll tissue. By stipulating that cell wall growth and remodelling occurs only near void space, our computational model is able to recapitulate spongy mesophyll development observed in Arabidopsis thaliana leaves. We find that robust generation of pore space in the spongy mesophyll requires a balance of cell growth, adhesion, stiffness and tissue pressure to ensure cell networks become porous yet maintain mechanical stability. The success of this mechanical model of morphogenesis suggests that simple physical principles can coordinate and drive the development of complex plant tissues like the spongy mesophyll.

Monitoring plant water status via static uniaxial compression of the leaf lamina

Turgor pressure is an essential, but difficult to measure indicator of plant water status. Turgor has been quantified by localized compression of cells or tissues, but a simple method to perform these measurements is lacking. We hypothesized that changes in leaf turgidity can be monitored by uniaxially compressing the leaf lamina and measuring the mechanical stress under a constrained thickness (stress relaxation) and that changes in leaf water content can be monitored by measuring the leaf thickness under constant mechanical stress. Using a simple, custom-built leaf squeeze-flow rheometer, we performed different compression tests on leaves from 13 plant species. The mechanical stress measured during stress relaxation was correlated with leaf bulk turgor pressure (R²  > 0.95) and thus with balancing pressure (R²  > 0.94); the leaf thickness measured under constant mechanical stress was correlated with relative water content (R²  > 0.74). The coefficients of these relationships were related to the leaf bulk osmotic pressure at the turgor-loss point. An idealized average-cell model suggests that, under isothermal conditions, the stationary bulk modulus during compression is largely determined by the bulk osmotic pressure. Our study presents an inexpensive, accessible and automatable method to monitor plant water status noninvasively.

Leaf cell-specific and single-cell transcriptional profiling reveals a role for the palisade layer in UV light protection

Like other complex multicellular organisms, plants are composed of different cell types with specialized shapes and functions. For example, most laminar leaves consist of multiple photosynthetic cell types. These cell types include the palisade mesophyll, which typically forms one or more cell layers on the adaxial side of the leaf. Despite their importance for photosynthesis, we know little about how palisade cells differ at the molecular level from other photosynthetic cell types. To this end, we have used a combination of cell-specific profiling using fluorescence-activated cell sorting and single-cell RNA-sequencing methods to generate a transcriptional blueprint of the palisade mesophyll in Arabidopsis thaliana leaves. We find that despite their unique morphology, palisade cells are otherwise transcriptionally similar to other photosynthetic cell types. Nevertheless, we show that some genes in the phenylpropanoid biosynthesis pathway have both palisade-enriched expression and are light-regulated. Phenylpropanoid gene activity in the palisade was required for production of the ultraviolet (UV)-B protectant sinapoylmalate, which may protect the palisade and/or other leaf cells against damaging UV light. These findings improve our understanding of how different photosynthetic cell types in the leaf can function uniquely to optimize leaf performance, despite their transcriptional similarities.

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.