A “Do-It-Yourself” phenotyping system: measuring growth and morphology throughout the diel cycle in rosette shaped plants
Plant Methods volume 13, Article number: 95 (2017)
Improvements in high-throughput phenotyping technologies are rapidly expanding the scope and capacity of plant biology studies to measure growth traits. Nevertheless, the costs of commercial phenotyping equipment and infrastructure remain prohibitively expensive for wide-scale uptake, while academic solutions can require significant local expertise. Here we present a low-cost methodology for plant biologists to build their own phenotyping system for quantifying growth rates and phenotypic characteristics of Arabidopsis thaliana rosettes throughout the diel cycle.
We constructed an image capture system consisting of a near infra-red (NIR, 940 nm) LED panel with a mounted Raspberry Pi NoIR camera and developed a MatLab-based software module (iDIEL Plant) to characterise rosette expansion. Our software was able to accurately segment and characterise multiple rosettes within an image, regardless of plant arrangement or genotype, and batch process image sets. To further validate our system, wild-type Arabidopsis plants (Col-0) and two mutant lines with reduced Rubisco contents, pale leaves and slow growth phenotypes (1a3b and 1a2b) were grown on a single plant tray. Plants were imaged from 9 to 24 days after germination every 20 min throughout the 24 h light–dark growth cycle (i.e. the diel cycle). The resulting dataset provided a dynamic and uninterrupted characterisation of differences in rosette growth and expansion rates over time for the three lines tested.
Our methodology offers a straightforward solution for setting up automated, scalable and low-cost phenotyping facilities in a wide range of lab environments that could greatly increase the processing power and scalability of Arabidopsis soil growth experiments.
There is an urgent need for novel technologies that monitor and predict the impact of abiotic (e.g. light, temperature) and biotic (e.g. pests, disease) stresses on plant growth and productivity at a large scale. Understanding the dynamic interactions between genotype, phenotype and the growth environment are key to predicting plant performance, resource use efficiency, stress tolerance and yield (for review see [1, 2]). Although molecular profiling technologies now enable the generation of large amounts of data with decreasing costs, the significantly slower process of phenotypic trait characterisation has become a bottleneck in advancing plant and agronomic science. To address this challenge, several automated non-destructive phenotyping platforms have been developed . Government investment in expensive, commercial high-throughput plant phenotyping tools (e.g. ) has led to the development of several state-of-the-art research hubs (e.g. the National Plant Phenomics Centre (www.plant-phenomics.ac.uk) in the UK). However, shared access is required and often at a premium, while hub locations can make direct research interaction challenging. Several academic and commercial lab-based tools also have been published for characterising plants, particularly the model species Arabidopsis thaliana (hereafter Arabidopsis). However, the wider uptake of such tools has been hampered by the local expertise required for hardware and software development. The lack of accessible hardware and appropriate algorithms for trait extraction is now widely considered the major bottleneck in the plant phenotyping field . Several online tools are available for measuring plant traits (e.g. www.plant-image-analysis.org ) to supplement either commercial systems or support affordable solutions (e.g. www.phenotiki.com ). However, many image analysis systems lack robust validation and are not well supported .
Developing capacity to quantify plant growth phenotypes throughout the plant life cycle is important for understanding dynamic growth traits. For example, plants grow differently during the day light–dark cycle (i.e. the diel cycle), and show significant changes in phenotypic traits such as carbon partitioning to biomass, expansion rates, and leaf movement [7,8,9,10,11,12,13,14]. Regulation of such traits is co-ordinated by internal metabolic status (e.g. carbon availability, source-sink dynamics), developmental status and organ age, and the response of signalling mechanisms (e.g. the circadian clock, plant hormones) to the external environment (e.g. light quality, temperature, the length of the photoperiod). Significant research efforts are focused on investigating and modelling these complex relationships to understand how plants optimise performance in response to a dynamic growth environment [15,16,17,18,19]. To date, published reports that monitor dynamic growth throughout the diel cycle have been restricted to expensive and/or custom built-hardware systems [20,21,22,23]. The most commonly used method of assessing Arabidopsis rosette growth rates is overhead (top-view) 2D image acquisition (Table 1) (although 3D systems are now becoming available [13, 17]). Following image acquisition, image processing enables the extraction of plant features from the image background, either through self-calibrated threshold analysis of colour and brightness or more complex computer-vision based methods (e.g. [5, 24]).
Here, we have developed an affordable image capture system (ICS) based on the Phenotiki approach  for the high-throughput phenotyping of soil-grown Arabidopsis plants throughout the diel cycle. The hardware is based on commonly available parts that are straightforward to order, and assembly requires some basic soldering and tools available at most fabrication facilities. The associated software module (iDIEL Plant) retains a high degree of versatility by allowing users to choose suitable parameters for experimental analysis rather than requiring adherence to a specific growth condition. To test the capacity of the ICS to capture growth data, and iDIEL Plant to identify and accurately measure growth for different plant phenotypes, we first performed a robust validation exercise and then tracked the growth of wild-type (WT) plants and two Rubisco-deficient mutants (1a2b and 1a3b) that differed in growth rate and leaf colour (i.e. chlorophyll content) . Our results demonstrated that this system offers a robust, low-cost solution for the uninterrupted capture of plant growth data throughout the diel cycle.
Construction of the image capture system (ICS)
The near infra-red (NIR) LED frame was designed in qCAD [www.qcad.org]) and fabricated using a 4 mm thick clear acrylic sheet [400 × 280 mm (L × W)] (Fig. 1a, b). A total of 173 NIR 1.3 V through-hole 5 mm LEDs with a peak emittance of 940 nm [Kingbright L-7113F3C, RS (www.uk.rs-online.com)] were mounted on the acrylic sheet through paired circular holes (1.4 mm in diameter, 2.6 mm apart) using an Epilog 60 W CO2 Laser Cutter (www.epiloglaser.co.uk). The NIR LEDs were positioned 30 mm apart along the length and 20 mm apart along the width of the acrylic sheet and were soldered together using solder wire (0.5 mm, 60/40 tin/lead) and a standard soldering iron with a fine point tip. A total of 28 axial fixed resistors were soldered in line with each parallel NIR LED circuit (see Additional file 1 for further details). The soldered wires were attached to the acrylic sheet using acrylic glue. The NIR LED circuit was powered via a 12 V DC 2A power supply adaptor. We designed four supporting legs [two perpendicular acrylic strips, 420 × 20 × 4 mm (L × W × B), glued together] and attached an acrylic foot base [40 × 40 × 4 mm (L × W × B)] to each leg (Epilog 60 W CO2 Laser Cutter). The legs were attached to the NIR LED frame using 3D printed corner braces [Wanhao Duplicator 4 (www.wanhao3dprinter.eu)]. The uniformity of NIR illumination provided by the rig and light transmittance through the acrylic frame were measured using a UV–Vis–NIR spectrometer [USB2000 + UV–Vis (www.oceanoptics.com)].
Images were acquired using a Raspberry Pi (RPi) camera module (RPi NoIR with a 5 MP OmniVision OV5647 sensor and IR filter removed) connected to a model 2B RPi computer (www.raspberrypi.org). The RPi NoIR was mounted centrally in the acrylic frame of the rig above the plant tray (Fig. 1a) and was operated and configured remotely using the RPi Cam-Web-Interface (www.elinux.org/RPi-Cam-Web-Interface) via a web browser on a Windows computer. The images were transferred remotely using Filezilla (www.filezilla-project.org/). The total weight of the ICS was approximately 400 g.
Plant materials and growth conditions
Initial validation of the ICS was made using Arabidopsis [Arabidopsis thaliana (L.) Heynh. Col-0] WT plants 15 days after germination (DAG) grown at 21 °C, ambient CO2, 70% relative humidity and 150 μmol photons m−2 s−1 in a 12: 12 h light: dark cycle. To further test the ICS, the growth phenotypes of Arabidopsis WT plants and two Arabidopsis Rubisco small subunit mutant lines, 1a2b and 1a3b (Col-0 background) , were compared. All plants in the latter study were grown from seeds of the same age and storage history, and were harvested from plants grown in the same environmental conditions to have a robust comparison of the three different genotypes. The seeds were sown in pots containing an organic compost-soil mix [F2 + S, Levington® (www.icl-sf.com)], stratified for 3 days at 4 °C, transferred to a side-lit growth cabinet [Snijders Scientific model ECD01E, (www.snijderslabs.com)] and grown in the same conditions described above. Seedlings (8 DAG) were transplanted to individual pots in a plant tray (each pot 50 × 50 mm) containing F2 + S. The transplanted plants were maintained in the same growth cabinet conditions.
During growth experiments, images were captured every 20 min throughout the diel cycle for 16 days (9–24 DAG). The NIR LED array was synchronised with the growth cabinet to turn on at the beginning of the dark period and turn off and the beginning of the next light period. The ICS and the plant tray were not moved throughout the experiment to preserve background consistency.
Segmentation of VIS and NIR images
For image analysis we developed a software module called iDIEL Plant. iDIEL Plant can segment Arabidopsis rosettes from both VIS (light) and NIR (dark) images and can extract quantitative estimates of projected rosette area (PRA). The software automatically numbers plants from left to right as they appear in the image, or manually by clicking on the centre of each plant. Thus, plants could be analysed regardless of arrangement in the image. ImageJ [v1.5n (imagej.nih.gov/ij)] was used to validate VIS and NIR image analyses. We further validated iDIEL Plant by examining segmentation accuracy from a published dataset of VIS and NIR Arabidopsis images . The iDIEL Plant module will be made available for free download.
For growth experiments, each light and dark period contained 35 images, which were batch analysed. Plants remained separated (i.e. rosette leaves of different plants did not touch or overlap) during the growth period. Segmentation of rosettes from VIS images was done by thresholding each rosette from the background (Fig. 1c, d). iDIEL Plant offers the flexibility of choosing between three different colour spaces (RGB, HSV and Lab) and grayscale Otsu  (i.e., only luminance) for thresholding, which increases capacity for finding a suitable colour space for any given image dataset. Each colour space is composed of three channels: RGB is composed of red, green and blue; HSV is composed of hue, saturation (S) and value (V); and Lab is composed of lightness, a (representing green/red) and b (representing blue/yellow) channels [37, 38]. Thresholding can be done by modifying the range of each channel in a colour space until the plants remain in the foreground and the background is segmented away. The effects of the modifications can be seen in real time to give a better indication of the effectiveness of segmentation. The VIS images in this study were segmented in the HSV colour space as it provided the best segmentation results for our setup.
NIR images of the plant tray contain information saturated on just the grayscale spectrum, preventing the use of HSV colour segmentation. Thus, to batch process the dark image sets, each NIR image was segmented using a Chan-Vese active contour method to separate foreground (plant) from the background . The Chan–Vese active contour algorithm is a region-based intensity driven model for delineating objects in images. Two regions are defined in the image by one or more curves that wrap around the chosen objects, dividing the image into foreground (inside the curve) and background (outside the curve). Each region is represented by a constant energy term, which is given by the average pixel value intensity within the region. The algorithm allows the curves to expand (or contract) to minimise the energy difference between the regions, resulting in segmentation of the chosen objects (i.e. rosettes). Since the Chan–Vese algorithm requires initial curves, the segmented rosettes from the last VIS image of each light period was taken as the initialization region of interest (ROI) mask for the active contour algorithm (Fig. 1c–e). The segmented rosettes from the following NIR image mask then was taken as the initialization mask for the subsequent NIR image for the remaining images in each set. As images were captured every 20 min, differences in plant appearance between sequential images were relatively small. The reverse approach can also be applied by using the opposite energy constraint (i.e. to contract) on the Chan-Vese formulation and using the first image of the light period as the starting mask and processing the NIR images in reverse starting from the end of the dark period. In this study both forward (day1-night-day2) and reverse (day2-night-day1) approaches were used and the data averaged.
LED placement and growth chamber characteristics can affect the quality of the NIR images.
iDIEL Plant is able to estimate and correct for uneven illumination (Fig. 1e). Firstly, grayscale pixel values are brought closer to the mean intensity values of the whole image. Secondly, a Gaussian blur is applied to the image to obtain a pixel intensity map of the lighter and darker regions. The pixel intensity values in the original image are then increased or decreased, respectively, based on the difference between the pixel intensity values of the intensity map and the mean intensity of the image map. The uneven illumination correction was used when analysing single channel grayscale images (i.e. the dark period and thresholding in the grayscale channel).
When segmentation was completed, the segmented image was correlated with the given plant indices and the area of each individual plant estimated. The data was exported to a Microsoft Excel file. Relative expansion rate (RER) was calculated in terms of PRA expansion according to Eq. 1 [13, 40]:
Daily RER was estimated using the PRAs obtained at the beginning and end of each light period. RER was also estimated over diel cycle using a sliding median window of 10 frames (i.e. approximately 3 h) .
Results and discussion
Performance of the image capture system (ICS)
The ICS consisted of a low-cost custom-built NIR LED frame and RPi-based image acquisition system. As the material costs of the rig were < £100, the ICS proved to be a relatively affordable but powerful tool for Arabidopsis phenotyping . The relatively small size and weight of the ICS simplified handling and transport. Furthermore, the qCAD design is malleable—the LED frame could be modified or arrayed to the end-users requirements.
The use of a clear acrylic frame minimised reflection and shading effects when images were taken during the light cycle. No measurable decrease in photosynthetically active radiation (PAR) was observed when the rig was set up in the side-lit growth cabinet used for growth experiments in this study. When light levels were measured under the rig in a vertically lit growth chamber (Snijders Scientific model MC1000), an average decrease of 14% in PAR was observed (180–154 μmol photons m−2 s−1). The reduction in PAR could be compensated for by adjusting the growth cabinet light output. Thus, we recommend testing light intensity under the rig for top-light or glass house experiments.
To test the uniformity of NIR illumination we measured light quality and quantity in the area under the ICS (Fig. 2). We observed a clear peak at 940 nm with a range of 850–1000 nm, and an average intensity of 30 μmol photons m−2 s−1. A small decrease in light intensity was observed in the central area around the camera, which had a reduced density of NIR LEDS. As expected, the light intensity also showed a drop off towards the edges of the ICS frame. However, the plant tray was located well within this zone. Despite the minor variations in NIR light intensity, image capture quality was sufficient to accurately segment all plant images throughout the growth experiment.
Raspberry Pi cameras have previously been shown to produce good quality VIS images for plant trait extraction and are a low-cost and easy-to-setup alternative to expensive cameras used in other plant imaging systems . Here, we found that the RPi NoIR camera was capable of capturing VIS and NIR images for image analysis. The ICS was able to acquire images at a maximum rate of ten images per minute. Furthermore, the camera could be controlled and image data transmitted remotely, giving significant flexibility to the user.
Data processing and software validation
VIS images can be batch processed by other available software tools for plant segmentation, such as Phenotiki and Rosette Tracker [5, 21]. However, we found that automatic segmentation of NIR images for soil-grown plants was challenging using available software tools. For example, software from Zhang et al.  suffered from low signal-to-noise ratios and removed night images during analysis. Rosette Tracker  relies on corresponding chlorophyll fluorescence or VIS images for manual annotation of dark images taken with a thermal IR camera, and Dhondt et al.  is optimised for imaging plants grown in agar-containing Petri dishes under relatively low PAR levels (60 μmol photons m−2 s−1). Therefore, we developed a MatLab-based software module called iDIEL Plant for segmentation and improved batch extraction of Arabidopsis rosette images from both VIS and NIR images.
Processing NIR images can be difficult due to increases in background noise (e.g. soil reflectance) when using thresholding compared to VIS images. Furthermore, NIR images are grayscale and cannot use the VIS image colour spaces and channels for thresholding (Fig. 1c, e), which makes it more challenging to define contours for segmentation. To overcome this problem, we relied on an active contour algorithm that uses a VIS image obtained shortly before/after the NIR image as an initialization mask . Previously, De Vylder et al.  developed a similar approach for the Rosette Tracker software, which uses a correspondent VIS image as a projected mask to analyse thermal infrared images (TIR). Rosette Tracker utilises a warping method that requires the user to manually identify and click on several regions in each VIS and IR corresponding image. For iDIEL Plant, image analysis was batched, such that the segmentation mask of each image was used as the initialization mask for the subsequent image. We found that this approach reduced background noise and provided a better initialisation point for the active contour algorithm.
The iDIEL Plant module currently extracts PRA estimates from VIS and NIR images, and could export segmentation masks and write raw and processed data in the database format required by Phenotiki . This permits the use of the Phenotiki software to extract additional growth traits (e.g. compactness, stockiness) and utilisation of other modules within Phenotiki (e.g. leaf counting and semi-automated leaf segmentation) .
Validation using ground truth measurements and/or published datasets is now considered a key requirement for developing plant imaging systems . The accuracy of iDIEL Plant was initially validated by comparing PRA estimates with two different manual measurement methods (Fig. 3a). Firstly, 2D image data (VIS and NIR) were manually contoured in ImageJ to calculate rosette areas. Secondly, we used a destructive (ground truth) approach, where whole rosettes were dissected and individual leaves of each plant were scanned. The destructive method excluded the presence of features that might cause discrepancies in the estimated area obtained by 2D imaging methods, such as leaf curvature, hyponasty, hidden petioles or overlapped leaves. iDIEL Plant was able to measure PRA with high accuracy for VIS and NIR images, showing estimates close to those obtained by manual segmentation and ground truth measurements. On average, the estimates from ground truth measurements were closer to those obtained from iDIEL Plant than those from ImageJ. The manual PRA estimates from ImageJ were not significantly different from values given by iDIEL Plant. To further validate the capacity of iDIEL Plant to accurately extract PRA estimates from available datasets, we used VIS and NIR plant images provided by Cruz et al.  (Fig. 3b, c). The dataset contained images of Arabidopsis Col-0 WT plants grown in a 16: 8 h light: dark cycle from 15 to 24 DAG on soil (i.e. a black background). Cruz et al.  provided raw images and rosette segmentation masks at a temporal resolution of four images per day. The masks were used as a representative proxy for PRA and were compared with the PRA estimates obtained from iDIEL Plant following analysis of the raw images. The raw VIS images and NIR images differed in resolutions and were not co-registered, thus we analysed them independently. VIS image datasets were batch processed in iDIEL Plant using thresholding in the HSV colour space. For processing NIR images, we manually calibrated the first image of each plant in the dataset (15 DAG) using grayscale Otsu thresholding , and then batch processed further images using the Chan-Vese active contour approach (see “Methods” section). Comparison between the provided masks and PRA estimates produced by iDIEL Plant showed no significant difference for both VIS and IR image datasets. These validation results also demonstrated that iDIEL Plant can accurately predict PRA from images obtained using different parameters from those used in our study.
Furthermore, we tested the segmentation performance of iDIEL Plant with different soil backgrounds (Additional file 2). Comparison of the PRAs of plants placed on four common soil mixture components (F2 + S, F2 + S mixed with perlite, coco peat, and vermiculite) showed that the software could accurately segment rosettes in VIS and NIR images with a variety of different backgrounds. However, we advise against using components with reflective particles, such as perlite of vermiculite, as these can sometimes be included in segmentation if touching the rosette, potentially leading to an artificial increase in PRA. If these are required for an experiment, we recommend using a black felt fabric cover or a layer of a homogeneous dark top soil to maximize segmentation accuracy.
Diel analyses of plant growth traits
Despite the reduced growth and pale leaf phenotypes of the two Rubisco mutants, iDIEL Plant successfully segmented and analysed 1a2b and 1a3b plants during the growth period tested. Following segmentation and trait extraction, we initially compared the end of day PRA for each genotype (Fig. 4a). WT plants showed a typical exponential growth pattern, with a PRA of 470 mm2 at 24 DAG. Both 1a2b and 1a3b showed a significantly decreased PRA compared to the WT plants throughout the growth experiment. At 24 DAG, 1a2b had a PRA of 150 mm2, while 1a3b was eightfold smaller than WT plants, reaching only 62 mm2. The reduced growth rates observed for 1a2b and 1a3b were in accordance with previous reports [35, 42].
Tracking PRA estimates at an increased temporal resolution (i.e. 20 min intervals) revealed additional growth details (Fig. 4b). We captured evidence of rhythmic leaf movements, characterised by the elevation (hyponasty) and lowering (epinasty) of the leaves during the diel cycle . Leaf movement was evident for WT plants by depressions in the growth curve at the beginning of the light period (i.e. the PRA of the 2D rosette images was reduced). Video analysis revealed that WT leaves moved downwards before moving upwards again during this period (Additional file 3). Our observations correlate well with the temporal leaf movement patterns seen in other studies with WT Arabidopsis plants [11, 13]. In contrast, similar leaf movement patterns were much reduced for 1a2b and not detected for 1a3b at the onset of the light period. This observation may indicate that leaf movement is inhibited for both Rubisco mutants. Alternatively, leaf movements may have been less pronounced compared to WT plants and were not detected using the ICS. Previous work has shown that leaves of the starchless pgm mutant have reduced and delayed hyponastic responses, suggesting that changes in carbon metabolism may impact on the rhythmicity of leaf movements . Thus, the Rubisco mutants used in the present study may be similarly affected by reduced carbon availability.
Next, we examined the 3-day average RER of each genotype over the growth period to investigate if the observed growth impairment in the Rubisco mutants was a result of a sustained reduction in RER or a decrease in RER during a particular developmental phase (Fig. 4a). WT plants had an increased RER compared to both mutants from 9 to 20 DAG. The RER of 1a3b was significantly lower than WT plants throughout this period. In contrast, 1a2b was on average lower than WT plants, but the observed reduction was not consistently significant. These dissimilarities in RER could be explained by the relative differences in Rubisco content for 1a2b and 1a3b, which are reduced by 50 and 30%, respectively [35, 42]. Rubisco content in C3 leaves is usually found in excess to facilitate dynamic adjustment to changes in the environment, such as increased irradiance or temperature. Even when Rubisco activity is reduced by up to 50%, plants may not show a significant change in growth phenotype under non-limiting growth conditions (i.e. the light levels used in this study were non-saturating) . Nevertheless, both Rubisco mutants showed a progressive improvement in RER as they matured and were more similar to WT values at 21 DAG (Figs. 4a, 5), indicating that the observed slow growth phenotypes of 1a2b and 1a3b was primarily due to a reduced RER in the early principal growth stages of rosette development .
The RER of each genotype was further investigated over the diel cycle at a higher temporal resolution during three different growth phases (Fig. 5). WT plants showed a relatively consistent RER pattern at all three growth phases, which included a rise in RER during the light period followed by a decline, and then a brief peak during the beginning of the dark period followed by a steady decrease until the following light period. RER fluctuations were much less pronounced in older WT plants. These growth rate patterns were similar to those observed for Arabidopsis rosettes in previous 2D studies and high precision 3D analyses of RER in diel cycles [13, 17, 20]. Both Rubisco mutants also showed an increase in RER during the light period, but the peak was much more discrete than WT plants, particularly for 1a3b. This observation is in accordance with the lower daily RER observed for these plants from 9 to 11 DAG (Fig. 4a). Surprisingly, younger Rubisco mutants showed an increase in RER during the dark period, typically with a higher peak than in the light period (Fig. 5a, b). In general, both Rubisco mutants showed a similar RER pattern, the difference being 1a3b had a comparatively reduced RER. This pattern persisted in the second growth phase (15–17 DAG), but gradually switched to a WT-like pattern during the third growth phase (21–23 DAG) (Fig. 5b, c). Previously, transgenic tobacco plants with reduced Rubisco levels were characterised by a decreased accumulation of sugars and starch during the light period compared to WT plants, but typically maintained WT-like rates of starch utilisation in the dark period, such that starch reserves in the transgenic plants were more depleted by dawn [9, 43]. The latter phenomenon was attributed to an increase in sink demand for carbon compared to WT plants  and may account for the RER patterns observed for the younger 1a3b and 1a2b mutants.
Our goal was to develop a simple, robust and affordable phenotyping system capable of tracking the growth of multiple plants during development throughout the diel cycle. The use of low cost components, easy-to-source hardware, and simplified construction requirements should significantly lower the entry barrier for plant researchers and facilitate scaling to more high-throughput phenotyping analyses . Our ICS design enabled the acquisition of dynamic and continuous image data, while the iDIEL Plant module could batch process image data for Arabidopsis genotypes that differed in size and hue. As a caveat, 2D imaging approaches should be considered only a proxy for 3D RER, as 2D data can convolute leaf growth, leaf movement and, in older rosettes, leaf overlap [5, 20, 21]. Furthermore, leaf tissue characteristics can vary between genotypes or in response to the growth environment, thus it is important to match 2D imaging of novel or uncharacterised genotypes with destructive analyses (e.g. fresh and dry rosette weight) to determine parameters such as specific leaf area, leaf thickness and density [16, 35]. Nevertheless, 2D RER appears to be a sensitive method for detecting relative changes in rosette shape, and could be used as a proxy to indicate differences between genotypes and/or environmental factors.
Furbank RT, Tester M. Phenomics–technologies to relieve the phenotyping bottleneck. Trends Plant Sci. 2011;16:635–44.
Fahlgren N, Gehan MA, Baxter I. Lights, camera, action: high-throughput plant phenotyping is ready for a close-up. Curr Opin Plant Biol. 2015;24:93–9.
Humplík JF, Lazár D, Husičková A, Spíchal L. Automated phenotyping of plant shoots using imaging methods for analysis of plant stress responses—a review. Plant Methods. 2015;11:29.
Arvidsson S, Pérez-Rodríguez P, Mueller-Roeber B. A growth phenotyping pipeline for Arabidopsis thaliana integrating image analysis and rosette area modeling for robust quantification of genotype effects. N Phytol. 2011;191:895–907.
Minervini M, Giuffrida MV, Perata P, Tsaftaris SA. Phenotiki: an open software and hardware platform for affordable and easy image-based phenotyping of rosette-shaped plants. Plant J. 2017;90:204–16.
Lobet G. Image analysis in plant sciences: publish then perish. Trends Plant Sci. 2017;22:559–66.
Wiese A, Christ M, Virnich O, Schurr U, Walter A. Spatio-temporal leaf growth patterns of Arabidopsis thaliana and evidence for sugar control of the diel leaf growth cycle. N Phytol. 2007;174:752–61.
Jansen M, Gilmer F, Biskup B, Nagel KA, Rascher U, Fischbach A, et al. Simultaneous phenotyping of leaf growth and chlorophyll fluorescence via Growscreen Fluoro allows detection of stress tolerance in Arabidopsis thaliana and other rosette plants. Funct Plant Biol. 2009;36:902–14.
Graf A, Schlereth A, Stitt M, Smith AM. Circadian control of carbohydrate availability for growth in Arabidopsis plants at night. Proc Natl Acad Sci. 2010;107:9458–63.
Bours R, van Zanten M, Pierik R, Bouwmeester H, van der Krol A. Antiphase light and temperature cycles affect PHYTOCHROME B-controlled ethylene sensitivity and biosynthesis, limiting leaf movement and growth of Arabidopsis. Plant Physiol. 2013;163:882–95.
Dornbusch T, Michaud O, Xenarios I, Fankhauser C. Differentially phased leaf growth and movements in Arabidopsis depend on coordinated circadian and light regulation. Plant Cell. 2014;26:3911–21.
Sulpice R, Flis A, Ivakov AA, Apelt F, Krohn N, Encke B, et al. Arabidopsis coordinates the diurnal regulation of carbon allocation and growth across a wide range of photoperiods. Mol Plant. 2014;7:137–55.
Apelt F, Breuer D, Nikoloski Z, Stitt M, Kragler F. Phytotyping4D: a light-field imaging system for non-invasive and accurate monitoring of spatio-temporal plant growth. Plant J. 2015;82:693–706.
Kölling K, Thalmann M, Müller A, Jenny C, Zeeman SC. Carbon partitioning in Arabidopsis thaliana is a dynamic process controlled by the plants metabolic status and its circadian clock. Plant Cell Environ. 2015;38:1965–79.
Chew YH, Wenden B, Flis A, Mengin V, Taylor J, Davey CL, et al. Multiscale digital Arabidopsis predicts individual organ and whole-organism growth. Proc Natl Acad Sci. 2015;112:E2556.
Weraduwage SM, Chen J, Anozie FC, Morales A, Weise SE, Sharkey TD. The relationship between leaf area growth and biomass accumulation in Arabidopsis thaliana. Front Plant Sci. 2015;6:1–21.
Apelt F, Breuer D, Olas JJ, Annunziata MG, Flis A, Nikoloski Z, et al. Circadian, carbon, and light control of expansion growth and leaf movement. Plant Physiol. 2017. doi:10.1104/pp.17.00503.
Chew YH, Seaton DD, Millar AJ. Multi-scale modelling to synergise plant systems biology and crop science. Field Crops Res. 2017;202:77–83.
Zhu XG, Lynch JP, Lebauer DS, Millar AJ, Stitt M, Long SP. Plants in silico: Why, why now and what? An integrative platform for plant systems biology research. Plant Cell Environ. 2016;39:1049–57.
Dhondt S, Gonzalez N, Blomme J, De Milde L, Van Daele T, Van Akoleyen D, et al. High-resolution time-resolved imaging of in vitro Arabidopsis rosette growth. Plant J. 2014;80:172–84.
De Vylder J, Vandenbussche F, Hu Y, Philips W, Van Der Straeten D. Rosette tracker: an open source image analysis tool for automatic quantification of genotype effects. Plant Physiol. 2012;160:1149–59.
Arend D, Lange M, Pape J-M, Weigelt-Fischer K, Arana-Ceballos F, Mücke I, et al. Quantitative monitoring of Arabidopsis thaliana growth and development using high-throughput plant phenotyping. Sci Data. 2016;3:160055.
Flood PJ, Kruijer W, Schnabel SK, van der Schoor R, Jalink H, Snel JFH, et al. Phenomics for photosynthesis, growth and reflectance in Arabidopsis thaliana reveals circadian and long-term fluctuations in heritability. Plant Methods. 2016;12:14.
Tomé F, Jansseune K, Saey B, Grundy J, Vandenbroucke K, Hannah MA, et al. rosettR: protocol and software for seedling area and growth analysis. Plant Methods. 2017;13:13.
Granier C, Aguirrezabal L, Chenu K, Cookson SJ, Dauzat M, Hamard P, et al. PHENOPSIS, an automated platform for reproducible phenotyping of plant responses to soil water deficit in Arabidopsis thaliana permitted the identification of an accession with low sensitivity to soil water deficit. N Phytol. 2006;169:623–35.
Skirycz A, Vandenbroucke K, Clauw P, Maleux K, De Meyer B, Dhondt S, et al. Survival and growth of Arabidopsis plants given limited water are not equal. Nat Biotechnol. 2011;29:212–4.
Bours R, Muthuraman M, Bouwmeester H, van der Krol A. OSCILLATOR: A system for analysis of diurnal leaf growth using infrared photography combined with wavelet transformation. Plant Methods. 2012;8:29.
Green JM, Appel H, Rehrig EM, Harnsomburana J, Chang J-F, Balint-Kurti P, et al. PhenoPhyte: a flexible affordable method to quantify 2D phenotypes from imagery. Plant Methods. 2012;8:45.
Zhang X, Hause RJ, Borevitz JO. Natural genetic variation for growth and development revealed by high-throughput phenotyping in Arabidopsis thaliana. G3 (Bethesda). 2012;2:29–34.
Ispiryan R, Grigoriev I, Zu Castell W, Schäffner AR. A segmentation procedure using colour features applied to images of Arabidopsis thaliana. Funct Plant Biol. 2013;40:1065–75.
Tessmer OL, Jiao Y, Cruz JA, Kramer DM, Chen J. Functional approach to high-throughput plant growth analysis. BMC Syst Biol. 2013;7(Suppl6):S17.
Easlon HM, Bloom AJ. Easy Leaf Area: automated digital image analysis for rapid and accurate measurement of leaf area. Appl Plant Sci. 2014;2:1400033.
Awlia M, Nigro A, Fajkus J, Schmoeckel SM, Negrão S, Santelia D, et al. High-throughput non-destructive phenotyping of traits that contribute to salinity tolerance in Arabidopsis thaliana. Front Plant Sci. 2016;7:1–15.
Cruz JA, Yin X, Liu X, Imran SM, Morris DD, Kramer DM, et al. Multi-modality imagery database for plant phenotyping. Mach Vis Appl. 2016;27:735–49.
Atkinson N, Leitão N, Orr DJ, Meyer MT, Carmo-Silva E, Griffiths H, et al. Rubisco small subunits from the unicellular green alga Chlamydomonas complement Rubisco-deficient mutants of Arabidopsis. N Phytol. 2017;214:655–67.
Otsu N. A threshold selection method from gray-level histograms. IEEE Trans Syst Man Cybern. 1979;9:62–6.
Philipp I, Rath T. Improving plant discrimination in image processing by use of different colour space transformations. Comput Electron Agric. 2002;35:1–15.
Tkalčič M, Tasič JF. Colour spaces—perceptual, historical and applicational background. In: Proceedings of IEEE Region 8 EUROCON 2003 Computer as a Tool 2003;A:304–8.
Chan TF, Vese LA. Active contours without edges. IEEE Trans Image Process. 2001;10:266–77.
Hunt R. Growth analysis, individual plants. In: Thomas B, Murphy DJ, Murray BG, editors. Encyclopedia of applied plant sciences. London: Academic Press; 2003.
Giuffrida MV, Minervini M, Tsaftaris S. Learning to count leaves in rosette plants. In: Tsaftaris SA, Scharr H, Pridmore T, editors. The computer vision problems in plant phenotyping. Guildford: BMVA Press; 2015. p. 1.1–1.13.
Izumi M, Tsunoda H, Suzuki Y, Makino A, Ishida H. RBCS1A and RBCS3B, two major members within the Arabidopsis RBCS multigene family, function to yield sufficient Rubisco content for leaf photosynthetic capacity. J Exp Bot. 2012;63:2159–70.
Quick WP, Schurr U, Fichtner K, Schulze ED, Rodermel SR, Bogorad L, Stitt M. The impact of decreased Rubisco on photosynthesis, growth, allocation and storage in tobacco plants which have been transformed with antisense rbcS. Plant J. 1991;1:51–8.
Boyes DC. Growth stage-based phenotypic analysis of Arabidopsis: a model for high-throughput functional genomics in plants. Plant Cell. 2001;13:1499–510.
AD, AJM designed and constructed the phenotyping rig and software. AD conducted the growth experiment. AD, LCTS conducted the validation experiments. LCTS, AD, AJM analysed and interpreted the data. LCTS, AJM, AD, SAT wrote the manuscript. AJM and SAT supervised the study. All authors read and approved the final manuscript.
We thank the UoE FabLab+ (http://www.fablab.eng.ed.ac.uk/) for assistance in hardware design and construction.
The authors declare that they have no competing interests.
Availability of data and materials
The iDiel Plant software module, source code and CAD files for the ICS design are hosted on the Edinburgh DataShare (http://datashare.is.ed.ac.uk/). Links can be found at http://mccormick.bio.ed.ac.uk/software/ and on the Edinburgh Predictive Plant website (http://predictiveplant.uk/).
Consent for publication
Ethics approval and consent to participate
LCTS was funded by the UK Biotechnology and Biological Sciences Research Council (Grants BB/M006468/1 and BB/N02334X/1 to AJM). AD was partly funded by the University of Edinburgh (UoE) Innovation Initiative Grant Scheme. AD is currently supported by an EPSRC DTP Ph.D. fellowship (EP/N509644/1). ST is partly supported by the UK Biotechnology and Biological Sciences Research Council (BB/N02334X/1 and BB/P023487/1).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Detailed schematics of the image capturing system. a) Dimensions of the acrylic sheet used to construct the NIR LED array (assembled in qCAD, [www.qcad.org]). Paired holes (1.4 mm in diameter) were positioned 2.6 mm apart to allow for mounting of standard through-hole 5 mm NIR LEDs (940 nm, Kingbright L-7113F3C, RS (www.uk.rs-online.com)]. The acrylic sheet was cut using an Epilog 60 W CO2 Laser Cutter (www.epiloglaser.co.uk). b) Schematic of the NIR LED circuit. The array is composed of 28 parallel rows of NIR LEDs powered by a 12 V/2A power supply. The double arrows next to the LED symbols (triangle and line) represent the direction of current flow from the cathode (positive terminal) to the anode (negative terminal). All LEDs must be soldered in the right orientation, such that the cathode of one LED (represented by a grey 1) is connected to the anode (represented by a grey 2) of the next LED. The resistor values are given in ohms (Ω).
Additional file 2. Projected rosette area estimations of Arabidopsis rosettes on different soil backgrounds. Three mature Arabidopsis rosettes were carefully excised and imaged on four common soil mixture components under VIS and NIR illumination. a) Example of an Arabidopsis rosette under the conditions specified. The NIR images show the outline of the active contour mask based on the segmentation of the corresponding VIS image. b) Projected rosette areas (PRAs) of rosettes on four different soil backgrounds (bars represent the mean ± SE; n = 3). No significant differences in PRA were observed (one-way ANOVA; p ≤ 0.05). However, note that for soils containing reflective particles (e.g. F2 + S + Perlite and Vermiculite) PRA can be slightly higher than more homogeneous soils, as the segmentation can sometimes include the particles that are directly under a leaf.
Additional file 3. Video showing growth and leaf movement for three Arabidopsis genotypes. Six plants for WT (middle) and Rubisco mutants 1a3b (top) and 1a2b (bottom) are shown from 18 to 19 DAG at 5 frames per second.
About this article
Cite this article
Dobrescu, A., Scorza, L.C.T., Tsaftaris, S.A. et al. A “Do-It-Yourself” phenotyping system: measuring growth and morphology throughout the diel cycle in rosette shaped plants. Plant Methods 13, 95 (2017). https://doi.org/10.1186/s13007-017-0247-6
- Arabidopsis thaliana
- Near-infrared LED array
- Image analysis
- Low cost