We provide a novel method to improve the X-ray absorption contrast of maize tissue suitable for ordinary microcomputed tomography scanning. Based on CT images, we introduce a set of image-processing workflows for different maize materials to effectively extract microscopic phenotypes of vascular bundles of maize.
It is necessary to accurately quantify the anatomical structures of maize materials based on high-throughput image analysis techniques. Here, we provide a ‘sample preparation protocol’ for maize materials (i.e., stem, leaf, and root) suitable for ordinary microcomputed tomography (micro-CT) scanning. Based on high-resolution CT images of maize stem, leaf, and root, we describe two protocols for the phenotypic analysis of vascular bundles: (1) based on the CT image of maize stem and leaf, we developed a specific image analysis pipeline to automatically extract 31 and 33 phenotypic traits of vascular bundles; (2) based on the CT image series of maize root, we set up an image processing scheme for the three-dimensional (3-D) segmentation of metaxylem vessels, and extracted two-dimensional (2-D) and 3-D phenotypic traits, such as volume, surface area of metaxylem vessels, etc. Compared with traditional manual measurement of vascular bundles of maize materials, the proposed protocols significantly improve the efficiency and accuracy of micron-scale phenotypic quantification.
The maize vascular system runs through the entire plant, from the root and stem to the leaves, which forms the key transportation paths for delivering water, mineral nutrients, and organic substances1. Another important function of the vascular system is to provide mechanical support for the maize plant. For example, the morphology, number, and distribution of vascular bundles in roots and stems are closely related to the lodging resistance of maize plants2,3. At present, studies on the anatomical structure of vascular bundles mainly utilize microscopic and ultramicroscopic techniques to display the anatomical structures of a certain part of the stem, leaf, or root, and then measure and count these structures of interest by manual investigation. Undoubtedly, manual measurement of various microscopic structures in large-scale microimages is a very tedious and inefficient work and severely limits the precision of microphenotypic traits, owing to its subjectivity and inconsistency4,5.
Maize has no secondary growth, and the cell content essentially consists of water in the primary meristem. Without any pretreatment, fresh samples of maize tissues can be directly scanned using a micro-CT device; however, the scanning results are probably poor and rough. The main reasons are summarized as follows: (1) low attenuation densities of plant tissues, resulting in a low contrast of atomic number and high noise in images; (2) fresh plant materials are prone to dehydrate and shrink during the normal scanning environment, as reported by Du6. The abovementioned problems have become the main constraints for the development and application of microphenotyping technology for maize, wheat, rice, and other monocotyledons. Here, we introduce the 'sample preparation protocol' to pretreat the samples of maize stem, leaf, and root. This protocol avoids the dehydration and deformation of plant materials during the CT scanning; thus, it is beneficial to increase the preservation time of plant samples with nondeformation. Moreover, the dyeing step based on solid iodine also enhances the contrast of plant materials; thus, it makes significant improvements in the imaging quality of micro-CT. Furthermore, we developed image processing software, named VesselParser, to process the CT images of maize stems and leaves. This software integrates a set of image-processing pipelines to perform high-throughput and automatic phenotyping analysis for 2-D CT images of different plant tissues. Vascular bundles in the entire cross-section of the maize stem and leaf are detected, extracted, and identified using an automatic image-processing method. As a result, we obtain 31 microscopic phenotypes of the maize stem and 33 microscopic phenotypes of the maize leaf. For the CT image series of the maize root, we developed an image-processing scheme to acquire 3-D phenotypic traits of metaxylem vessels. This scheme is superior in efficiency of image acquisition and reconstruction compared with traditional methods.
These results indicate that the image processing pipelines considering the imaging characteristics of ordinary X-ray micro-CT provide an effective method for the microscopic phenotyping of vascular bundles; this extremely widens the applications of CT techniques in plant science and improves the automatic phenotyping of plant materials at cellular resolution6,7.
1. Sample Preparation Protocol
2. Micro-CT Scanning Protocol
3. Image Analysis Protocol for a Single CT Image of a Maize Stem or Leaf
NOTE: Use automatic imaging software for vascular bundles to conduct the phenotyping analysis of vascular bundles within the CT slice images of the maize stem and leaf (Figure 2). The software usage steps are described as follows.
4. Image Analysis Protocol for CT Image Series of a Maize Root
NOTE: The CT image series of maize roots are utilized to extract the 3-D structures of metaxylem vessels using image-processing software. The main steps are as follows.
The sample preparation protocol suitable for ordinary micro-CT scanning not only prevents the deformation of plant tissues but also enhances the X-ray absorption contrast. Pretreated plant materials are scanned using a micro-CT system into high-quality slice images, and the highest resolution can reach 2 µm/pixel. Figure 4 shows the scanned micro-CT images of stem, leaf, and root, and the image contrast has a significant improvement compared with the results scanned from the fresh plant materials. In these CT slice images, significant differences in gray-level values are observable among vascular bundles, parenchyma cells, xylem vessels, epidermis, and other tissues.
Based on this imaging software for vascular bundles, the structure and distribution features of vascular bundles within the cross-section of the entire stem or leaf can be automatically analyzed and quantified. Take the image of a maize stem as an example, and the phenotyping algorithm of software consists of five steps: segment the slice image, segment the vascular bundles, analyze the spatial distribution of the vascular bundles, identify and improve the vascular bundles, and compute and output the phenotypic traits of the vascular bundles. The entire epidermis region in the image is, first, segmented based on specific or adaptive threshold values, and contour analysis techniques are performed to extract the epidermis region (Figure 3A and 3B). Next, segmentation of the vascular bundles is conducted. Vascular bundles are surrounded in the epidermis; therefore, the epidermis region with the predefined thickness is removed from the slice image. The rest of the image only consists of scattered vascular bundles with high pixel intensities and parenchyma cells with weak pixel intensities. Then, the image segmentation based on the fixed threshold value is performed to extract all candidate regions of the vascular bundles, and these regions are further determined as valid vascular bundles according to the limitations of area and shape features (Figure 3C and 3D).
After the segmentation of vascular bundles, spatial distribution characteristics of the vascular bundles are extracted. Take the geometric centers of scattered vascular bundles as nodes to generate triangle meshes for all vascular bundles in the slice image, and these meshes are clustered into five types according to their areas. The areas of the vascular bundles manifest a significant decreasing trend from the center to the edge of the maize stem. The triangular and Voronoi meshes describe the spatial distribution and topological connections of vascular bundles, and each mesh is drawn with a specific color according to the clustered results of vascular bundles (Figure 3E – 3H). Vascular bundles that meet the constraints of spatial distribution (mesh area and shape are the important indexes to determine the availability of vascular bundle) are reserved and used to generate the final segmentation results (Figure 3I).
At last, phenotypic traits of vascular bundles, such as geometric, shape, and distribution information, can be computed according to the above analysis, which results in the output of a TXT or CSV file (Figure 3J). Based on the imaging software for vascular bundles, 31 phenotypic traits of the stem can be automatically analyzed; the average computation time for each CT image is ~30 s. The relative phenotypic parameters of the stem are shown in Table 1. Similarly, 33 phenotypic traits of the leaf can be extracted, the average computation time is ~50 s, and these parameter classifications are shown in Table 2. For a list of slice images, the above image analysis pipelines are integrated into a batch processing for automatic execution. This workflow is efficient to analyze the phenotypic traits of all the vascular bundles within an entire slice image of a maize stem and leaf. Notably, most phenotypic traits of vascular bundles, such as the total area, average area, and area ratio of vascular bundles, are significantly difficult to be measured by manual measurements.
Because the metaxylem vessels of maize roots show obvious morphological changes along the direction of root growth, it is more valuable to extract the 3-D structures of metaxylem vessels for phenotypic analysis. Based on the CT image series of maize root, 3-D segmentation, surface reconstruction, and volume visualization are performed. Based on the segmented results, the 3-D structural parameters of metaxylem vessels can be automatically calculated, including the volume, surface area, cross-sectional area (basal), and cross-sectional area (distal) of the total metaxylem vessels and each single metaxylem vessel. This workflow significantly improves the efficiency of 3-D phenotypic trait analysis. The results of segmentation, reconstruction, and volume visualization can directly manifest the spatial structures of metaxylem vessels of maize root, as shown in Figure 5.
Figure 1: Different types of sample baskets printed with a 3-D printer. (A and B) The spiral sample basket for the leaf is composed with (A1) a central oval groove and (A2) the surrounding spiral groove. The width of the spiral groove is set at about 4 mm to accommodate the leaf with the main vein. The sidewall of the spiral groove is created with (A5) square drainage holes and the bottom of the basket with (B6) circular drain holes. (C and D) These two panels show a four-well sample basket suitable for the stem, with (C1) sample holes with a diameter of 25 mm, and the bottom of the basket with (D2) seven drain holes. (E and F) These two panels show a multiwell sample basket suitable for the root, with (E1) four circular holes with a diameter of 10 mm in the central area, and (E2) 13 circular holes with a diameter of 8 mm arranged near the edge of the basket. The bottom of the sample basket has drain holes with a diameter of 1 mm to ensure that the tiny root tissue does not leak out. Please click here to view a larger version of this figure.
Figure 2: Screenshots of automatic imaging software for vascular bundles. (A) This panel shows the data management to import CT images from any file directory and select CT slice images for the subsequent processing. (B) This panel shows the method parameters to identify the organ type of CT slice images and configure the corresponding method parameters. (C) This panel shows the phenotyping computation to perform the batch computation of CT slice images and show the execution progress. (D) This panel shows the statistical analysis to verify the computational results and generate phenotypic traits for all the CT images. (E) This panel shows the computational results to output the analysis results as a TXT or CVS file. Please click here to view a larger version of this figure.
Figure 3: Workflow of image-processing and phenotyping computation based on CT slice images. (A) Import the CT slice image of the maize stem. (B) Segment the slice image with a fixed threshold value. (C) Extract the entire region of the maize stem. (D) Delete the epidermis region of the maize stem. (E) Triangular meshes of vascular bundles. (F) This panel shows a clustering analysis according to the areas of the vascular bundles. (G) This panel shows a clustering analysis according to the areas of the triangular meshes. (H) This panel shows a clustering analysis according to the areas of Voronoi meshes. (I) This panel shows the final segmented results of the vascular bundles. (J) This panel shows the spatial distribution of the vascular bundles. Please click here to view a larger version of this figure.
Figure 4: Multiplanar reconstruction (MPR) images of the maize stem, leaf, and root. The left panel shows an MPR image of a stem cross-section. The middle panel shows an MPR image of a leaf cross-section. The right panel shows an MPR image of a root cross-section. Please click here to view a larger version of this figure.
Figure 5: 3-D visualization of root metaxylem vessels. The bars = 0.2 mm. (A and B) These panels show a 3-D visualization of a root. (C – F) These panels show a 3-D visualization of metaxylem vessels. Please click here to view a larger version of this figure.
Table 1: Microscopic phenotypic traits of the maize stem.
Table 2: Microscopic phenotypic traits of the maize leaf.
With the successful application of CT technology in the fields of biomedicine and materials science, this technology has been gradually introduced into the fields of botany and agriculture, promoting researches in plant life sciences as a promising technical tool. In the late 1990s, CT technology was first used to study the morphological structures and development of plant root systems. In the past decade, synchrotron HRCT has become a powerful, nondestructive tool for plant biologists, and has been successfully used to identify the tissue structures of grape vascular system8, tissue structure of Arabidopsis leaf9,10, and seed structure of rape11. Through the synchrotron HRCT, great progress has been made in the study of structure and function of vascular bundles in woody plants12,13,14. However, little research of HRCT technology for maize, wheat, rice, and other crops has been done15. Maize has no secondary growth, and the cell essentially consists of water in the primary meristem. Although fresh samples can be scanned by micro-CT without any pretreatment, the scanning results are very poor. The main reasons are as follows: (1) low attenuation densities of plant tissues, resulting in a low contrast in atomic number and high noise in images; (2) fresh material tends to dehydrate and shrink during the scanning period, as reported by Du6. The above-mentioned reasons have become the main factors limiting the application of this technology in maize, wheat, rice, and other monocotyledons.
Here, we introduce a simple and practical sample preparation protocol that not only prevents the deformation of plant tissues but also enhances the X-ray absorption contrast. High-quality and high-resolution CT images of root, stem, and leaf were obtained based on the sample preparation protocol and micro-CT imaging system, and the highest resolution was up to 2 µm/pixel. Thus, the sample preparation protocol is suitable for ordinary micro-CT scanning and provides a great opportunity for broader applications in monocot and other plant sciences. This protocol can be easily modified to accommodate other plant materials such as dehydration or drying procedure, and its parameter setting can also be adjusted according to specific plant materials for the best results. Notably, this approach is limited by the size and volume of the plant sample. A too thick sample segment may lead to incomplete drying or deformation of the sample. Therefore, this sample preparation protocol is applicable for small plant materials with a thickness of less than 3 cm and not for much larger maize materials such as maize ear or cob.
The microscopic phenotyping technology of plant materials is one of the hot topics of plant phenotypic studies in recent years, and it is gradually becoming one of the basic support technologies for genetic breeding and plant physiology. Traditional microscopic phenotypic analysis of plants requires a large number of complex sample preparations and tedious manual operations. It is very labor-intensive and time-consuming to count and measure the microscopic traits; the results are also prone to subjective errors. For example, to quantify the phenotypic characteristics of vascular bundles in a maize stem, the stem sample needs to be embedded in paraffin and then sliced, stained, and imaged. For a stained slice image, it is difficult to perform an automated image processing due to the ambiguous definition of cell boundaries; thus, manual identification and segmentation are indispensable16. To satisfy the requirements of large-scale measurements for anatomical traits of the maize stem, Legland and Heckwolf introduced different image-processing methods; however, the structures of vascular bundles in the rind are still a challenge17,18. Hence, high-throughput image analysis and accurate quantification of anatomical traits of maize tissues are necessary. Here, we provide automatic imaging software for vascular bundles, that can automatically extract 31 phenotypic traits of vascular bundles in ~30 s for each CT image of the maize stem and 33 phenotypic traits in ~50 s for each CT image of the maize leaf. The structure and distribution features of vascular bundles within the cross-section of entire stems or leaves can be automatically analyzed and quantified. This software has the following advantages: (1) it automatically processes the CT slice images of the maize stem and leaf and extracts the phenotypic traits of vascular bundles; (2) it has a higher recognition rate of vascular bundles in the CT image, especially for small vascular bundles on the edge; (3) a novel graphic analysis method is used to reveal the distribution characteristics of vascular bundles.
Further, the X-ray micro-CT scanning technology has obvious advantages in the efficiencies of image acquisition and reconstruction compared with traditional reconstruction techniques based on paraffin section images19,20,21. Based on the CT image series of maize root, an image-processing scheme is developed to extract the spatial structures of metaxylem vessels and successfully used for the 3-D measurement of microscopic traits. The main limitation of this scheme is that the 3-D segmented results depend on slight manual interaction. In the future, we aim to develop a set of automated 3-D image analysis software for the CT dataset of maize roots to improve the efficiency of 3-D segmentation and reconstruction.
In conclusion, an ordinary micro-CT scanning based on a practical sample preparation protocol for the maize stem, leaf, and root is constructed to produce high-resolution CT images. The sample preparation protocol provided here not only prevents the deformation of plant tissues but also enhances the X-ray absorption contrast. This protocol is also suitable for other CT scanning applications in wheat, rice, and other monocotyledons. So far, we have developed automatic imaging software for vascular bundles, that is capable of automatically and rapidly extracting the phenotypic traits of vascular bundles from a single CT image of the maize stem and leaf. Based on the CT image series of the maize root, an image-processing scheme is successfully set up to extract the 3-D phenotypic traits of metaxylem vessels. Microscopic phenotyping techniques of plant materials based on X-ray micro-CT provide a new prospect for the accurate and rapid quantification and identification of maize vascular bundles.
The authors have nothing to disclose.
This research was supported by the National Nature Science Foundation of China (No.31671577), the Science and Technology Innovation Special Construction Funded Program of Beijing Academy of Agriculture and Forestry Sciences(KJCX20180423), the Research Development Program of China (2016YFD0300605-01), the Beijing Natural Science Foundation (5174033), the Beijing Postdoctoral Research Foundation (2016 ZZ-66), and the Beijing Academy of Agricultural and Forestry Sciences Grant (KJCX20170404), (JNKYT201604).
Skyscan 1172 X-ray computed tomography system | Bruker Corporation, Belgium | NA | For CT scanning |
CO2 critical point drying system (Leica CPD300) | Leica Corporation, Germany | NA | For sample drying |
Ethanol | Any | NA | For FAA fixation |
Formaldehyde | Any | NA | For FAA fixation |
Acetic acid | Any | NA | For FAA fixation |
Surgical blade | Any | NA | For cutting the sample sgements |
3D printer | Makerbot replicator 2, MakerBot Industries, USA | NA | For printing the sample baskets of maize root, stem, and leaf |
Centrifuge tube | Corning, USA | NA | Place the root, stem, or leaf materials |
Solid iodine | Any | NA | For sample dyeing |
SkyScan Nrecon software | SkyScan NRecon, Version: 1.6.9.4, Bruker Corporation, Belgium | NA | For image reconstruction |
VesselParser software | VesselParser, Version: 3.0, National Engineering Research Center for Information Technology in Agriculture (NERCITA), Beijing, China | NA | Image analysis protocol for single CT image of maize stem or leaf |
ScanIP | ScanIP, Version: 7.0; Simpleware, Exeter, UK | NA | 3D image processing software |
Latex gloves | Any | NA | |
Tweezers | Any | NA |