- Open Access
Investigating the microstructure of plant leaves in 3D with lab-based X-ray computed tomography
- Andrew W. Mathers†1,
- Christopher Hepworth†2,
- Alice L. Baillie†2,
- Jen Sloan†2,
- Hannah Jones2,
- Marjorie Lundgren3,
- Andrew J. Fleming2,
- Sacha J. Mooney1 and
- Craig J. Sturrock1Email author
© The Author(s) 2018
- Received: 31 July 2018
- Accepted: 3 November 2018
- Published: 12 November 2018
Leaf cellular architecture plays an important role in setting limits for carbon assimilation and, thus, photosynthetic performance. However, the low density, fine structure, and sensitivity to desiccation of plant tissue has presented challenges to its quantification. Classical methods of tissue fixation and embedding prior to 2D microscopy of sections is both laborious and susceptible to artefacts that can skew the values obtained. Here we report an image analysis pipeline that provides quantitative descriptors of plant leaf intercellular airspace using lab-based X-ray computed tomography (microCT). We demonstrate successful visualisation and quantification of differences in leaf intercellular airspace in 3D for a range of species (including both dicots and monocots) and provide a comparison with a standard 2D analysis of leaf sections.
We used the microCT image pipeline to obtain estimates of leaf porosity and mesophyll exposed surface area (Smes) for three dicot species (Arabidopsis, tomato and pea) and three monocot grasses (barley, oat and rice). The imaging pipeline consisted of (1) a masking operation to remove the background airspace surrounding the leaf, (2) segmentation by an automated threshold in ImageJ and then (3) quantification of the extracted pores using the ImageJ ‘Analyze Particles’ tool. Arabidopsis had the highest porosity and lowest Smes for the dicot species whereas barley had the highest porosity and the highest Smes for the grass species. Comparison of porosity and Smes estimates from 3D microCT analysis and 2D analysis of sections indicates that both methods provide a comparable estimate of porosity but the 2D method may underestimate Smes by almost 50%. A deeper study of porosity revealed similarities and differences in the asymmetric distribution of airspace between the species analysed.
Our results demonstrate the utility of high resolution imaging of leaf intercellular airspace networks by lab-based microCT and provide quantitative data on descriptors of leaf cellular architecture. They indicate there is a range of porosity and Smes values in different species and that there is not a simple relationship between these parameters, suggesting the importance of cell size, shape and packing in the determination of cellular parameters proposed to influence leaf photosynthetic performance.
- X-ray computed tomography (microCT)
- Leaf structure
- Air channels
- Gas exchange
It is estimated that a doubling in agricultural productivity will be required over the next three decades to meet the increasing food demand of a rapidly growing global population . Photosynthesis is an important driver of food production but has thus far been little improved by crop breeding or engineering . Although significant advances have recently begin to be reported via engineering photosynthetic biochemistry [3–5], less progress has been made in the optimisation of internal leaf architecture (the arrangement of cells and airspaces within the leaf) which is also thought to limit photosynthetic carbon assimilation .
For example, the surface area of mesophyll cells exposed to intercellular airspaces (Smes) has been shown to be positively correlated with photosynthetic performance [7, 8] presumably by facilitating increased diffusional flux of CO2. A clearer understanding of how leaf architectural traits, such as Smes and porosity (the proportion of the leaf volume occupied by airspace), influence photosynthetic potential is vital if we are to successfully optimise leaf cellular architecture to maximise carbon assimilation.
However, imaging the microstructure of plant leaves can be challenging due to their low density, fine structure, and sensitivity to desiccation. A number of established stereological approaches are commonly used to quantify leaf structural parameters, such as Smes, from two-dimensional (2D) tissue cross-sections of chemically fixed, resin embedded tissue. Measurements of lengths or areas from the cross-sectional images are transformed using correction factors to generate estimations of three-dimensional (3D) geometry [9–11]. However these techniques are destructive, labour intensive, and in some cases the process of tissue preparation can alter the parameters being measured (e.g. poorly sectioned samples or sectioned at oblique angles), potentially leading to underestimation of values for Smes by as much as 30% .
More recently, 3D imaging techniques have been applied to simplify and improve quantification of plant structures. Tomographic techniques generate non-destructive serial section images through the sample of interest. A range of tomography-based techniques is now available for imaging of low density materials, the majority of which were developed in medical physics as non-invasive diagnostic tools. For example, nuclear magnetic resonance imaging (MRI) allows visualisation of materials based on their water content. It has been successfully applied to studies of seedling germination, plant root growth and architecture in soils [13–18], but its relatively coarse spatial resolution (> 50 μm) makes it unsuitable for imaging the fine microstructure of aerial plant tissues. Positron emission tomography (PET) uses short-lived radioisotopes (typically 11C and 14C) to determine the assimilation of compounds in living organisms with exceptionally high sensitivity (picomolar order of magnitude). However, the spatial resolution of PET is even more coarse than MRI (1–5 mm), so structural information must be gathered independently . Optical projection tomography (OPT) uses visible light (photons) to discriminate between materials. Whilst OPT is capable of capturing high resolution images (ca. 5 μm), like other optical techniques it is limited by sample thickness and requires chemical fixation and staining of tissues  which can often be time consuming and place limitations on throughput. This technique is, however, useful for studying the spatial distribution of marker gene expression in stained plant tissues, as demonstrated by Lee et al. . Finally, X-ray computed microtomography (microCT) combines the advantages of high resolution and excellent depth penetration by using X-rays to visualise structure. The X-ray attenuation coefficient of a material is dictated by its density and atomic number , so the technique is capable of imaging plant tissue structures by discriminating low density intercellular airspaces from denser cellular material .
MicroCT can be conducted in synchrotron facilities or using more compact, lab-based equipment. Synchrotron-based microCT (SRXCT) has the advantage of using a high flux, coherent, monochromatic photon beam permitting collection of both absorption and phase contrast radiographic images at high resolution (e.g. 0.35–5 μm image pixel size. TOMCAT Beamline, Swiss Light Source ). This technique has been successfully used to discriminate individual plant cells and to investigate airspace connectivity in fresh fruit tissue  and leaves . However, the expense and scarcity of such facilities limits the use of synchrotron-based microCT. Although lab-based microCT systems can now achieve similar resolution range to SRXCT, the greater accessibility of lab based equipment has allowed it to be used to study many plant structural features such as trichome distribution on Arabidopsis leaves , leaf venation , panicle development and seed density in rice , floral shape variation in orchids  and volume and surface area measurements of inflorescences of tulips and proteaceae . As benchtop microCT systems have a lower X-ray flux, sample damage is generally considered to be lower compared to Synchrotron based systems. It is therefore possible to perform non-destructive imaging of live plants, allowing repeat measurements on the same individuals over time or before and after a treatment. However, Dhondt et al.  reported inhibition of seedling growth after multiple rounds of scanning, suggesting that there is a limit to the intensity and/or frequency of scanning of live tissue that is possible without affecting development. In microCT systems, the low level of X-ray absorption by plant tissue presents challenges to differentiate cellular level structures such as individual cell types due to insufficient image contrast. This can be overcome to some extent by the use of low energies. The application of contrast agent solutions (e.g. iodine, gadolinium, barium) have also provided promising results to overcome this issue [23, 31], but the use of contrast agents do not appear to increase image quality in all systems , and can lead to longer preparation times compared to scanning fresh tissue.
Several studies have used the image data generated by microCT to calculate quantitative descriptors of plant tissue structure. For example Schneider et al.  used microCT images to calculate vein density in leaf tissue, and Herremans et al.  conducted a very detailed analysis of fruit tissue structure. We have recently used microCT to quantify leaf cellular architecture of Arabidopsis mutants, uncovering relationships between structural parameters and photosynthetic performance [34, 35]. Here, we demonstrate that lab-based X-ray microCT can be used to visualise and quantify differences in leaf intercellular airspace in 3D in a range of species including both dicots and monocots. Our method yields high resolution images (ca. 2.5–2.75 μm) and does not require laborious chemical fixation or staining techniques to prepare the samples. An image analysis pipeline has been developed to provide quantitative descriptors of plant leaf intercellular airspace. We have focused on leaf porosity and Smes, but methods for further structural analyses are included in the Additional file 1. These 3D measurements provide insight into the available pathways for gas flow within the leaf, which in turn influences the potential photosynthetic productivity of the plant.
Arabidopsis seeds (Col-0 ecotype) were sown directly into 60 × 60 × 80 mm pots of damp, lightly compressed soil (3:1 Levington M3 compost:perlite) and stratified at 4 °C for 7 days before transfer into a controlled environment chamber (Conviron, Canada) under short day conditions (12 h light 22 °C/12 h dark 15 °C, 200 µmol m−2 s−1 PAR at rosette level, 60% humidity). Leaf discs were excised from the largest leaves for scanning 30 days after germination. Tomato (Solanum lycopersicum var. Ailsa Craig), pea (Pisum sativum var. Arvense), barley (Hordeum vulgare var. Tipple Fulbourn) and oat (Avena sativa), were sown in 20 × 20 × 30 cm pots of M3 compost and grown under long day conditions (16 h light 22 °C/8 h dark 15 °C, 200 µmol m−2 s−1 PAR, 60% humidity). Leaf discs were excised from the largest, mature leaves for scanning. Rice seeds (Oryza latifolia) were germinated on wet filter paper in 90 mm diameter, 20 mm deep petri plates, and transplanted into water-saturated soil (70% v/v Kettering Loam (Boughton, UK), 23% v/v Vitax John Innes No. 3 (Leicester, UK), 5% v/v silica sand and 2% v/v Osmocote Extract Standard 5–6 month slow release fertiliser (Ipswich, UK)) in 105 × 105 × 185 mm pots, 8 days after germination. Rice plants had a constant water supply from the pot base and were grown in a controlled environment chamber (Conviron, Canada) with 12 h, 30 °C days and 12 h 24 °C nights, 700 µmol m−2 s−1 PAR at canopy level and 60% relative humidity.
For rice and Arabidopsis, n = 5. For the other species, and for the 2D analysis of rice, leaf sections n = 4. To allow comparison of leaves of the same species (or mutants) we selected leaf 5 for analysis in our experimental studies so they are at same developmental growth stage.
Sample preparation for microCT
X-ray microCT scanning
Single microCT scans of leaf discs were performed using a GE Phoenix Nanotom S 180NF (GE Sensing and Inspection Technologies GmbH, Wunstorf, Germany) fitted with a tungsten transmission target and a 5 MP (2304 × 2304 pixel) CMOS digital detector (Hamamatsu Photonics KK, Shizuoka, Japan). A three-point detector calibration was performed, collecting an average of 100 images, with 10 skip images per gain point. Scans were obtained at a spatial resolution of 2.75 μm (2304 × 1400 pixel field of view), with an electron acceleration energy of 85 kV and a current of 100 μA (higher spatial resolutions are possible if a smaller diameter sample can be used). Detector exposure time was 500 ms, collecting 3600 projections in ‘fast scan’ mode (sample rotates continuously), with no averaging or skip images and no pixel binning (1 × 1), resulting in a scan duration of 30 min per sample.
Radiograph reconstruction was carried out using Phoenix Datos|x rec 2 reconstruction software (version 2.3.3; GE Sensing and Inspection Technologies GmbH, Wunstorf, Germany). Radiographs were assessed for sample movement using the autoscan|optimiser module, by comparing the difference between the first and last projection image (0° and 360° rotation) and applying an automatic directional and/or scale correction if movement and/or shrinkage were apparent. Any sample that required more than 3 pixel shift in x or y axis were either rescanned or disregarded as the image quality in these images was low. Beam hardening artefacts were mitigated using the multiple materials function in the BHC + module. A beam hardening correction of 7 was determined to be the most appropriate for plant leaves. Finally, radiographs were manually cropped i.e. resized to remove the scanned area beyond the leaf sample before being reconstructed into 3D volumes using a filtered back-projection algorithm.
An illustration of the image analysis workflow is provided in Fig. 1e–i.
Alignment and cropping
Grayscale volumes were aligned in 3D (adaxial leaf surface facing up), cropped to remove any damaged leaf material at the disc periphery, and converted to stacks of TIFF images in the Z dimension using VG StudioMAX (version 2.2.0; Volume Graphics GmbH, Heidelberg, Germany).
Leaf discs were segmented from the surrounding sample holder by creating material masks in Avizo Fire software (version 6.0.0 Fire; Thermo Fisher Scientific, USA), using the ‘Label Field’ function and then binarising the selection.
Individual grayscale TIFF stacks were thresholded using the ‘Threshold’ function in the open source software package ImageJ (version 1.48; ) and saved as binary TIFF stacks, differentiating solid material from airspace. The automated thresholding algorithm was selected based on comparison between the binarised and the greyscale images, to account for small differences between scans in sample/background contrast, leaf water content and polystyrene elements. Previous research by our group has shown that the IJ Iso-data algorithm proved effective for thresholding Arabidopsis . However, it should be highlighted that a range of automated thresholding algorithms are available within ImageJ and will result in different outputs depending on the grayscale distributions of the image. This unfortunately, results in some level of manual selection of the most appropriate threshold algorithm. We would strongly recommend that the same threshold algorithm is used for all samples within the same study. For the rice and cereal leaves, the Li algorithm was used as they presented a finer pore structure. Material masks were thresholded using the automatic thresholding method ‘MaxEntropy’. All thresholded images were saved as binary TIFF stacks.
Intercellular airspace extraction
Binary material masks were combined with thresholded image stacks using the ‘Image Calculator’ function in ImageJ to create a composite image stack, isolating the extracellular airspace within each leaf disc.
Scans were de-noised using the ‘Remove Outliers’ function in ImageJ. Foreground and background particles < 3 × the spatial resolution were removed.
Region of interest selection
The inclusion of the mid-rib and/or major veins in images subjected to 3D analysis can artificially increase porosity measurements. In monocots, where vasculature is arranged in parallel cell files, regions of interest were selected between major veins. In rice in particular, which has dense vasculature, three 200 × 200 voxel regions were selected for analysis, and all 3D measurements were averaged across these technical replicates to provide representative data for the leaf disc as a whole. In all other species a region of interest (ROI) of ≥ 400 × 400 voxels was used. Due to the non-uniform structure and irregular vasculature of dicot leaves, it was not possible to entirely exclude vasculature, but the largest veins were avoided.
Representative 3D renderings of plant material, with air channel diameters illustrated by heat map, were constructed in VG StudioMAX (version 2.2.0; Volume Graphics GmbH, Heidelberg, Germany) using the isosurface and Phong rendering tools. The heat map data was an output of the ‘Thickness’ function in the ImageJ plugin BoneJ (version 1.3.14; ) which also provides the mean and maximum channel diameter for each stack.
Sample preparation for 2D analysis of fixed tissue sections
Leaf discs of O. latifolia were fixed in 4% v/v formaldehyde in PEM buffer (1.5% w/v Pipes, 0.19% w/v EGTA, 0.124% w/v MgSO4, pH 7) immediately after CT scanning. After no more than 72 h, samples were rinsed in PEM buffer three times for 10 min each. Samples were dehydrated in an ascending ethanol series (10%, 30%, 50%, 70%, 90%, 100% v/v ethanol, 1 h each) then infiltrated with an ascending series of LR white resin (London Resin Company) in ethanol (10%, 20%, 30%, 50%, 70%, 90% v/v 1 h each then 3 × 8 + hours in 100% resin). Samples were kept at 4 °C throughout dehydration and infiltration. Finally samples were stood vertically in gelatine capsules filled with resin and left to polymerise for 5 days at 37 °C. 2 µm sections were cut with a Reichert-Jung Ultracut E ultramicrotome and dried onto vectabond-coated multi-well slides. 4–5 sections were imaged per biological replicate, each of which was at least a cell’s length apart. Sections were stained for 5 min in a 0.1 mg mL−1 solution of propidium iodide in water and rinsed in water before imaging. Samples were imaged using a Leica DM6 microscope and camera equipped with a CoolLED fluorescence system, and images were captured using LASX software. Samples were illuminated with the 535 nm LED line, and visualised through the Y3 filter.
The workflow for stereological analysis is illustrated in Additional file 2: Fig. S1. Masks representing total leaf area (Additional file 2: Fig. S1B) and individual airspaces (Additional file 2: Fig. S1C) were generated using ImageJ (FIJI v1.51u;  with the connection thresholding and edge detection plugins). Masks were smoothed using the Median filter, with a radius of 3 pixels. Airspace area was expressed as a percentage of total leaf area to give an estimate of porosity (the fraction of leaf volume occupied by intercellular airspace).
All statistical analyses were conducted in Graphpad Prism software (version 7.03).
3D analysis of leaves from common dicot and monocot reference species
The quantification of the surface area of mesophyll cells exposed to intercellular airspace (Smes) allows testing of established ideas about the importance of this factor in CO2 uptake. We calculated Smes for each of the six species (Fig. 3b). The dicots with the highest mean values of Smes were pea (18.1 ± 1.2 mm2 mm−2) and tomato (18.7 ± 1.0 mm2 mm−2), significantly higher than Arabidopsis, which had the lowest value of all six species (15.0 ± 0.2 mm2 mm−2). Barley had the greatest Smes value (21.3 ± 0.3 mm2 mm−2), significantly higher than the other two monocots (oat 16.5 ± 1.1 and rice 15.5 ± 0.7 mm2 mm−2). Both porosity and Smes measurements were highly reproducible between biological replicates, as demonstrated by the low standard error values across the range of species.
Comparison of 2D and 3D quantification of rice leaf cellular architecture
Structural variation within leaves
Leaf cellular architecture is known to play an important role in photosynthesis. With the development of more advanced equipment, software and protocols, such as those described here, it is now possible to visualise leaf internal air channels at sub-micron resolution, and to quantify biologically relevant aspects of the air channel network. This method allows rapid imaging of live tissue samples at high resolution. Our previous work has successfully employed this technique for the characterisation of Arabidopsis mutants [34, 35] and here we demonstrate that it can be successfully applied to a wider range of plant species.
The relatively rapid scan time is a key advantage of this method, as it allows the use of live tissue. Crucially, the leaf structure does not change through desiccation during the course of such a rapid scan providing it is adequately supported by the radio opaque polystyrene foam. Fast scanning has the additional advantages of allowing a higher throughput rate than other microCT protocols, and mitigating problems of X-ray induced damage that could occur with prolonged or repeated scanning . The use of live tissue minimises preparation time and removes the risk of artefacts that could be introduced by fixation and staining. Established stereological methods generally use embedded tissue, which risks structural changes during fixation or dehydration stages of the embedding process. Furthermore, the much lower tissue coverage in 2D approaches compared to tomography tends to lead to underestimation of Smes, as demonstrated by Théroux-Rancourt et al.  who sampled 2D slices from their 3D image stacks for a robust comparison. Our comparison of 2D and 3D analysis of identical rice leaf discs confirmed previous reports that 2D approaches lead to lower Smes estimates. It should be stressed that during the image analysis procedure careful testing of the most appropriate automated threshold algorithm must be investigated. An algorithm suitable for one plant species may not be appropriate for another and lead to under or over estimations of porosity and Smes. The very low standard errors for each group of biological replicates suggest that our protocol for microCT analysis is robust and reliable within plant species where the same threshold algorithm was applied.
While the tissue is live at the time of scanning, our method does require destructive sampling. Repeated scans of the same sample after a treatment or over the course of development are therefore not possible. MicroCT scanners are available to scan larger samples, but this comes with a trade-off in resolution that would prevent accurate extraction of such small structures as the leaf airspace network . Furthermore, holding the sample sufficiently stationary to obtain a clear image without detachment from the plant presents a challenge. Excision of leaf discs does result in an area of collapsed tissue around the edge of the sample, but this can be excluded when selecting the ROI for analysis. The disc must be handled with great care during mounting to ensure that no further damage occurs. We make the assumption that any wound response in the tissue does not result in a change in the structure of the center of the disc during the timescale of the scan.
Current limitations on the achievable contrast and resolution of live tissue with lab-based microCT equipment (X-ray absorption based) prevent individual cells from being distinguished in the images. Finding a way to resolve cellular detail would be a challenging but useful target for the future, offering insights into the developmental processes that lead to the formation of the airspace network. Scanning at even greater resolution can produce stacks in which individual cells can be seen in live tissue, but this has only been demonstrated with synchrotron-based microCT to date [25, 41, 42]. Alternatively, increasing the contrast between cells and tissues, or boosting the contrast of the cell outlines (cell wall/cell membrane) using phase contrast techniques have demonstrated improved edge detection of cellular features in plant roots . Dhondt et al.  used iodine as a contrast agent to obtain detail at the individual cell level, but this required more extensive tissue preparation and a much slower scan time. Even if appropriate contrast agents were available for use with live samples, infiltrating them through the full tissue depth would be challenging. Until such a method is available, classical histological techniques will remain useful to complement the microCT data. Recently, combined microCT and histological approaches have provided valuable insights in biomedical studies [44, 45]. The adoption of similar approaches to plants may enable sub-cellular structures to be revealed, such as plastid size and position, which are highly relevant for understanding photosynthesis but cannot be obtained by tomography alone.
Selection of the ROI for computational analysis is critical for obtaining realistic and comparable numerical data. Firstly, damaged tissue areas must be avoided. It is also desirable to avoid veins as far as possible, although this brings a trade-off with ROI size. In these analyses we sought to use the largest ROI possible without inclusion of major veins. In the rice samples, the veins were sufficiently close together that we took multiple ROIs from each scan for the analysis to sample a sufficiently large vein-free area. Taking multiple ROIs per sample is a more labour-intensive approach as each region must be image processed separately. However, a further advantage of smaller ROI is that there is a greater possibility of selecting a flat region of leaf, which in turn makes the separation of distinct leaf layers as sets of z-slices more feasible, such as palisade and spongy mesophyll in dicots. Smaller ROIs also make it possible to avoid other structures, such as large trichomes, which might skew airspace quantification. While structures such as veins and trichomes are currently a complication that we have tried to avoid in our analysis, they are part of the true leaf structure. As our models of leaf development and of gas exchange networks advance, the inclusion of these features in ROIs may become useful and informative.
The image processing workflow presented here allows for the largely automated calculation of many morphological descriptors of the extent and spatial patterning of the leaf airspace network. However, some stages of the analysis still require manual verification by the operator, which are slower and more subjective than the automated steps. Generating the mask to define the tissue volume (as distinct from background and packing elements) is a semi-automated process, but in some species required extensive manual input. In Arabidopsis, for example, masking areas with trichomes requires some manual input if these leaf hairs are to be excluded so as not to affect the quantitative data. Furthermore, the density of some areas of the polystyrene packing discs is similar to that of the Arabidopsis cells (especially if the polystyrene has been compressed), requiring that these regions adjoining the tissue be manually removed from the masks. Defining the automated threshold value for image binarisation is also a manual step, and therefore somewhat subjective. Unfortunately, the availability of suitable ground truthing techniques to support the decision-making process is sadly lacking and therefore the informed ‘expert’ assessment of a trained user is required on a species by species basis. After this initial decision has been reviewed an automated analysis routine can be implemented.
The quantification and spatial mapping of leaf airspace allows us to probe the relationship between structure and function in the leaf by measuring gas exchange in plant lines that vary in their cellular architecture. Dorca-Fornell et al.  reported that, in Arabidopsis plants with altered expression of a cell cycle regulator, an increase in leaf porosity led to a significant increase in stomatal conductance but, perhaps surprisingly, no related increase in carbon assimilation. In contrast, demonstrated a positive correlation between mesophyll cell density and photosynthetic capacity among Arabidopsis mutants with cell cycle gene expression manipulated in targeted tissue layers. These results suggest a complex relationship between leaf structure and photosynthesis, which may involve effects of both gas exchange and light attenuation. Network analysis of the air channels within the leaf may contribute to modelling of airspace arrangements that might increase mesophyll conductance, allowing more effective gas exchange.
In addition to investigating the effects of leaf structure on gas exchange, these morphological data can be used to investigate relationships between leaf structural parameters. We might expect, for example, that porosity and Smes would be inter-dependent, but our data suggest this relationship is not so simple. This becomes a question of cell packing which will be influenced by a range of factors, such as the size and shape of mesophyll cells and the local control of cell separation. For example, increasing the extent of lobing in rice mesophyll cells could elevate Smes without greatly changing porosity. Understanding how to manipulate the development of such elements of leaf structure (and having a robust means of quantifying the output structural parameters) is essential if we are to use such information to manipulate leaf structure with a view to optimising photosynthetic performance.
The method presented here allows for the high resolution imaging of leaf intercellular airspace networks by lab-based microCT, and the largely-automated, quantitative description of those networks. These data can be used to investigate both developmental phenomena, such as the inter-relation of structural parameters, and physiological questions about the effect of leaf structure on gas exchange and photosynthesis. Technical advances in microCT imaging (e.g. improved contrast X-ray detectors with higher sensitivity) may, in the future, offer possibilities for gathering even more detailed information on leaf structure from live tissue, such as resolving individual cells. However, data at the level of detail that can currently be achieved already offer much unexplored potential for testing established ideas and developing new hypotheses to establish which leaf structural features are the most important for photosynthesis. This understanding, combined with knowledge of leaf developmental genetics, could facilitate the re-engineering of the leaf to enhance plant productivity.
AWM, ALB, JS and HJ conducted the experiments and analysed the data. AWM, ALB, JS, CH and CJS wrote the manuscript. All authors reviewed and edited the manuscript. All authors read and approved the final manuscript.
Oat, barley, tomato and pea plants were grown by Alice Mitchell and scanning was performed by Dr. Radek Pajor. The image J macros and R scripts used for analysis of 2D sections were based on those developed by Dr. Rachel Thorley. The authors would like to thank Dr. Michael Pound (Centre for Plant Integrative Biology, University of Nottingham) for his expertise and advice on measurement methods within ImageJ for microCT images.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
This work was supported by the Biotechnology and Biological Sciences Research Council (Grant Nos. BB/J004065/1, BB/N013719/1) the first grant for ‘Optimising photosynthetic efficiency via leaf structure’ and the second for ‘A 3D model of photosynthesis to inform breeding for improved rice performance in a changing climate’.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ray DK, Mueller ND, West PC, Foley JA. Yield trends are insufficient to double global crop production by 2050. PLoS ONE. 2013;8(6):e66428.View ArticleGoogle Scholar
- Long SP, Marshall-Colon A, Zhu X-G. Meeting the global food demand of the future by engineering crop photosynthesis and yield potential. Cell. 2015;161(1):56–66.View ArticleGoogle Scholar
- Glowacka K, Kromdijk J, Kucera K, Xie JY, Cavanagh AP, Leonelli L, et al. Photosystem II Subunit S overexpression increases the efficiency of water use in a field-grown crop. Nat Commun. 2018;9:868.View ArticleGoogle Scholar
- Kromdijk J, Glowacka K, Leonelli L, Gabilly ST, Iwai M, Niyogi KK, et al. Improving photosynthesis and crop productivity by accelerating recovery from photoprotection. Science. 2016;354(6314):857–61.View ArticleGoogle Scholar
- Long BM, Hee WY, Sharwood RE, Rae BD, Kaines S, Lim YL, et al. Carboxysome encapsulation of the CO2-fixing enzyme Rubisco in tobacco chloroplasts. Nat Commun. 2018;9:3570.View ArticleGoogle Scholar
- Mathan J, Bhattacharya J, Ranjan A. Enhancing crop yield by optimizing plant developmental features. Development (Cambridge, England). 2016;143(18):3283–94.View ArticleGoogle Scholar
- Evans JR, Kaldenhoff R, Genty B, Terashima I. Resistances along the CO2 diffusion pathway inside leaves. J Exp Bot. 2009;60(8):2235–48.View ArticleGoogle Scholar
- Park S. Internal leaf area and cellular CO2 resistance: photosynthetic implications of variations with growth conditions and plant species. Physiol Plant. 1977;40:137–44.View ArticleGoogle Scholar
- Turrell FM. The area of the internal exposed surface of dicotyledon leaves. Am J Bot. 1936;23(4):255–64.View ArticleGoogle Scholar
- Thain JF. Curvature correction factors in the measurement of cell surface areas in plant tissues. J Exp Bot. 1983;34(138):87–94.View ArticleGoogle Scholar
- James SA, Smith WK, Vogelmann TC. Ontogenetic differences in mesophyll structure and chlorophyll distribution in Eucalyptus globulus ssp. globulus. Am J Bot. 1999;86(2):198–207.View ArticleGoogle Scholar
- Theroux-Rancourt G, Earles JM, Gilbert ME, Zwieniecki MA, Boyce CK, McElrone AJ, et al. The bias of a two-dimensional view: comparing two-dimensional and three-dimensional mesophyll surface area estimates using noninvasive imaging. New Phytol. 2017;215(4):1609–22.View ArticleGoogle Scholar
- Van As H, Scheenen T, Vergeldt FJ. MRI of intact plants. Photosynth Res. 2009;102(2–3):213–22.PubMedPubMed CentralGoogle Scholar
- Metzner R, Eggert A, van Dusschoten D, Pflugfelder D, Gerth S, Schurr U, et al. Direct comparison of MRI and X-ray CT technologies for 3D imaging of root systems in soil: potential and challenges for root trait quantification. Plant Methods. 2015;11:17.View ArticleGoogle Scholar
- Schmittgen S, Metzner R, Van Dusschoten D, Jansen M, Fiorani F, Jahnke S, et al. Magnetic resonance imaging of sugar beet taproots in soil reveals growth reduction and morphological changes during foliar Cercospora beticola infestation. J Exp Bot. 2015;66(18):5543–53.View ArticleGoogle Scholar
- Metzner R, van Dusschoten D, Bühler J, Schurr U, Jahnke S. Belowground plant development measured with magnetic resonance imaging (MRI): exploiting the potential for non-invasive trait quantification using sugar beet as a proxy. Front Plant Sci. 2014;5:469.View ArticleGoogle Scholar
- Li K, Song W, Zhu L. Observation and measurement of plant root architecture in situ: a review. Shengtaixue Zazhi. 2011;30(9):2066–71.Google Scholar
- Eberhard M, Hardy R, Steffen O-J, Johannes F, André G, Thomas N, et al. A functional imaging study of germinating oilseed rape seed. New Phytol. 2017;216(4):1181–90.View ArticleGoogle Scholar
- Garbout A, Munkholm LJ, Hansen SB, Petersen BM, Munk OL, Pajor R. The use of PET/CT scanning technique for 3D visualization and quantification of real-time soil/plant interactions. Plant Soil. 2012;352(1–2):113–27.View ArticleGoogle Scholar
- Sharpe J. Optical projection tomography. Annu Rev Biomed Eng. 2004;6:209–28.View ArticleGoogle Scholar
- Lee K, Avondo J, Morrison H, Blot L, Stark M, Sharpe J, et al. Visualizing plant development and gene expression in three dimensions using optical projection tomography. Plant Cell. 2006;18(9):2145–56.View ArticleGoogle Scholar
- Flannery BP, Deckman HW, Roberge WG, D’Amico KL. Three-dimensional X-ray microtomography. Science (New York, NY). 1987;237(4821):1439–44.View ArticleGoogle Scholar
- Dhondt S, Vanhaeren H, Van Loo D, Cnudde V, Inzé D. Plant structure visualization by high-resolution X-ray computed tomography. Trends Plant Sci. 2010;15(8):419–22.View ArticleGoogle Scholar
- Marone F, Mokso R, Modregger P, Fife J, Pinzer B, Thuring T, et al. Present and future X-ray tomographic microscopy at TOMCAT. In: McNulty I, Eyberger C, Lai B, editors. 10th international conference on X-ray microscopy. AIP conference proceedings, vol. 1365; 2011. p. 116–9.Google Scholar
- Verboven P, Kerckhofs G, Mebatsion HK, Ho QT, Temst K, Wevers M, et al. Three-dimensional gas exchange pathways in pome fruit characterized by synchrotron X-ray computed tomography. Plant Physiol. 2008;147(2):518–27.View ArticleGoogle Scholar
- Kaminuma E, Yoshizumi T, Wada T, Matsui M, Toyoda T. Quantitative analysis of heterogeneous spatial distribution of Arabidopsis leaf trichomes using micro X-ray computed tomography. Plant J Cell Mol Biol. 2008;56(3):470–82.View ArticleGoogle Scholar
- Schneider JV, Rabenstein R, Wesenberg J, Wesche K, Zizka G, Habersetzer J. Improved non-destructive 2D and 3D X-ray imaging of leaf venation. Plant Methods. 2018;14:7.View ArticleGoogle Scholar
- Jhala VM, Thaker VS. X-ray computed tomography to study rice (Oryza sativa L.) panicle development. J Exp Bot. 2015;66(21):6819–25.View ArticleGoogle Scholar
- Van der Niet T, Zollikofer CP, León MS, Johnson SD, Linder HP. Three-dimensional geometric morphometrics for studying floral shape variation. Trends Plant Sci. 2010;15(8):423–6.View ArticleGoogle Scholar
- Stuppy WH, Maisano JA, Colbert MW, Rudall PJ, Rowe TB. Three-dimensional analysis of plant structure using high-resolution X-ray computed tomography. Trends Plant Sci. 2003;8(1):2–6.View ArticleGoogle Scholar
- Staedler YM, Masson D, Schoenenberger J. Plant tissues in 3D via X-ray tomography: simple contrasting methods allow high resolution imaging. PLoS ONE. 2013;8(9):e75295.View ArticleGoogle Scholar
- Tracy SR, Gomez JF, Sturrock CJ, Wilson ZA, Ferguson AC. Non-destructive determination of floral staging in cereals using X-ray micro computed tomography (microCT). Plant Methods. 2017;13:9.View ArticleGoogle Scholar
- Herremans E, Verboven P, Verlinden BE, Cantre D, Abera M, Wevers M, et al. Automatic analysis of the 3-D microstructure of fruit parenchyma tissue using X-ray micro-CT explains differences in aeration. BMC Plant Biol. 2015;15:264.View ArticleGoogle Scholar
- Dorca-Fornell C, Pajor R, Lehmeier C, Pérez-Bueno M, Bauch M, Sloan J, et al. Increased leaf mesophyll porosity following transient retinoblastoma-related protein silencing is revealed by microcomputed tomography imaging and leads to a system-level physiological response to the altered cell division pattern. Plant J Cell Mol Biol. 2013;76(6):914–29.View ArticleGoogle Scholar
- Lehmeier C, Pajor R, Lundgren MR, Mathers A, Sloan J, Bauch M, et al. Cell density and airspace patterning in the leaf can be manipulated to increase leaf photosynthetic capacity. Plant J Cell Mol Biol. 2017;92(6):981–94.View ArticleGoogle Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.View ArticleGoogle Scholar
- Doube M, Klosowski MM, Arganda-Carreras I, Cordelieres FP, Dougherty RP, Jackson JS, et al. BoneJ Free and extensible bone image analysis in ImageJ. Bone. 2010;47(6):1076–9.View ArticleGoogle Scholar
- Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82.View ArticleGoogle Scholar
- Giuliani R, Koteyeva N, Voznesenskaya E, Evans MA, Cousins AB, Edwards GE. Coordination of leaf photosynthesis, transpiration, and structural traits in rice and wild relatives (genus Oryza). Plant Physiol. 2013;162(3):1632–51.View ArticleGoogle Scholar
- Legland D, Devaux M-F, Guillon F. Quantitative imaging of plants: multi-scale data for better plant anatomy. J Exp Bot. 2018;69(3):343–7.View ArticleGoogle Scholar
- Karunakaran C, Lahlali R, Zhu N, Webb AM, Schmidt M, Fransishyn K, et al. Factors influencing real time internal structural visualization and dynamic process monitoring in plants using synchrotron-based phase contrast X-ray imaging. Sci Rep. 2015;5:12119.View ArticleGoogle Scholar
- Keyes SD, Daly KR, Gostling NJ, Jones DL, Talboys P, Pinzer BR, et al. High resolution synchrotron imaging of wheat root hairs growing in soil and image based modelling of phosphate uptake. New Phytol. 2013;198(4):1023–9.View ArticleGoogle Scholar
- Koebernick N, Daly KR, Keyes SD, George TS, Brown LK, Raffan A, et al. High-resolution synchrotron imaging shows that root hairs influence rhizosphere soil structure formation. New Phytol. 2017;216(1):124–35.View ArticleGoogle Scholar
- Hopkins TM, Heilman AM, Liggett JA, LaSance K, Little KJ, Hom DB, et al. Combining micro-computed tomography with histology to analyze biomedical implants for peripheral nerve repair. J Neurosci Methods. 2015;255:122–30.View ArticleGoogle Scholar
- Girard R, Zeineddine HA, Orsbon C, Tan H, Moore T, Hobson N, et al. Micro-computed tomography in murine models of cerebral cavernous malformations as a paradigm for brain disease. J Neurosci Methods. 2016;271:14–24.View ArticleGoogle Scholar