Stomata are portals in plant leaves that control gas exchange for photosynthesis, a process fundamental to life on Earth. Gas fluxes and plant productivity depend on external factors such as light, water and CO2 availability and on the geometrical properties of the stoma pores. The link between stoma geometry and environmental factors has informed a wide range of scientific fields-from agriculture to climate science, where observed variations in stoma size and density are used to infer prehistoric atmospheric CO2 content. However, the physical mechanisms and design principles responsible for major trends in stomatal patterning are not well understood. Here, we use a combination of biomimetic experiments and theory to rationalize the observed changes in stoma geometry. We show that the observed correlations between stoma size and density are consistent with the hypothesis that plants favour efficient use of space and maximum control of dynamic gas conductivity, and that the capacity for gas exchange in plants has remained constant over at least the last 325 Myr. Our analysis provides a new measure to gauge the relative performance of species based on their stomatal characteristics.

The mesophyll surface area exposed to intercellular air space per leaf area (S_m ) is closely associated with CO₂ diffusion and photosynthetic rates. S_m is typically estimated from two-dimensional (2D) leaf sections and corrected for the three-dimensional (3D) geometry of mesophyll cells, leading to potential differences between the estimated and actual cell surface area. Here, we examined how 2D methods used for estimating S_m compare with 3D values obtained from high-resolution X-ray microcomputed tomography (microCT) for 23 plant species, with broad phylogenetic and anatomical coverage. Relative to 3D, uncorrected 2D S_m estimates were, on average, 15-30% lower. Two of the four 2D S_m methods typically fell within 10% of 3D values. For most species, only a few 2D slices were needed to accurately estimate S_m within 10% of the whole leaf sample median. However, leaves with reticulate vein networks required more sections because of a more heterogeneous vein coverage across slices. These results provide the first comparison of the accuracy of 2D methods in estimating the complex 3D geometry of internal leaf surfaces. Because microCT is not readily available, we provide guidance for using standard light microscopy techniques, as well as recommending standardization of reporting S_m values.

Maintaining high rates of photosynthesis in leaves requires efficient movement of CO₂ from the atmosphere to the mesophyll cells inside the leaf where CO₂ is converted into sugar. CO₂ diffusion inside the leaf depends directly on the structure of the mesophyll cells and their surrounding airspace, which have been difficult to characterize because of their inherently three-dimensional organization. Yet faster CO₂ diffusion inside the leaf was probably critical in elevating rates of photosynthesis that occurred among angiosperm lineages. Here we characterize the three-dimensional surface area of the leaf mesophyll across vascular plants. We show that genome size determines the sizes and packing densities of cells in all leaf tissues and that smaller cells enable more mesophyll surface area to be packed into the leaf volume, facilitating higher CO₂ diffusion. Measurements and modelling revealed that the spongy mesophyll layer better facilitates gaseous phase diffusion while the palisade mesophyll layer better facilitates liquid-phase diffusion. Our results demonstrate that genome downsizing among the angiosperms was critical to restructuring the entire pathway of CO₂ diffusion into and through the leaf, maintaining high rates of CO₂ supply to the leaf mesophyll despite declining atmospheric CO₂ levels during the Cretaceous.