Email updates

Keep up to date with the latest news and content from BMC Plant Biology and BioMed Central.

Open Access Highly Accessed Software

GiA Roots: software for the high throughput analysis of plant root system architecture

Taras Galkovskyi1, Yuriy Mileyko1, Alexander Bucksch23, Brad Moore4, Olga Symonova5, Charles A Price6, Christopher N Topp47, Anjali S Iyer-Pascuzzi47, Paul R Zurek47, Suqin Fang478, John Harer17, Philip N Benfey47 and Joshua S Weitz29*

Author Affiliations

1 Department of Mathematics, Duke University, Durham, NC, USA

2 School of Biology, Georgia Institute of Technology, Atlanta, GA, USA

3 School of Interactive Computing, Georgia Institute of Technology, Atlanta, GA, USA

4 Department of Biology, Duke University, Durham, NC, USA

5 , Institute of Science and Technology, Vienna, Austria

6 Department of Plant Biology, University of Western Australia, Perth, Australia

7 Duke Center for Systems Biology, Duke University, Durham, NC, USA

8 Root Biology Center, South China Agricultural University, Guangzhou, China

9 School of Physics, Georgia Institute of Technology, Atlanta, GA, USA

For all author emails, please log on.

BMC Plant Biology 2012, 12:116  doi:10.1186/1471-2229-12-116

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2229/12/116


Received:13 January 2012
Accepted:28 June 2012
Published:26 July 2012

© 2012 Galkovskyi et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Characterizing root system architecture (RSA) is essential to understanding the development and function of vascular plants. Identifying RSA-associated genes also represents an underexplored opportunity for crop improvement. Software tools are needed to accelerate the pace at which quantitative traits of RSA are estimated from images of root networks.

Results

We have developed GiA Roots (General Image Analysis of Roots), a semi-automated software tool designed specifically for the high-throughput analysis of root system images. GiA Roots includes user-assisted algorithms to distinguish root from background and a fully automated pipeline that extracts dozens of root system phenotypes. Quantitative information on each phenotype, along with intermediate steps for full reproducibility, is returned to the end-user for downstream analysis. GiA Roots has a GUI front end and a command-line interface for interweaving the software into large-scale workflows. GiA Roots can also be extended to estimate novel phenotypes specified by the end-user.

Conclusions

We demonstrate the use of GiA Roots on a set of 2393 images of rice roots representing 12 genotypes from the species Oryza sativa. We validate trait measurements against prior analyses of this image set that demonstrated that RSA traits are likely heritable and associated with genotypic differences. Moreover, we demonstrate that GiA Roots is extensible and an end-user can add functionality so that GiA Roots can estimate novel RSA traits. In summary, we show that the software can function as an efficient tool as part of a workflow to move from large numbers of root images to downstream analysis.

Background

Plant roots are essential to the structure and function of plants. They provide access to belowground water and nutrients [1-3], facilitate anchorage of plants in soils [4,5], mediate chemical defense belowground [6,7], and serve as sites of important symbioses with microbiota [8,9]. However, comparatively little is known about the architecture of plant root systems and their relationship to overall plant function given the relative inaccessibility of belowground tissue [10]. Recently, standardized field procedures have been developed for quantifying RSA using excavated and washed roots [11]. In addition, a number of non-destructive approaches to imaging and quantifying RSA have been developed. The approaches include X-ray computed tomography [12-15], nuclear magnetic resonance (NMR) microscopy [16], magnetic resonance imaging [17,18], laser scanning [19], as well as a growing number of imaging methods using plants grown in transparent media [20-25]. These methods are part of a rapidly growing field of “plant phenomics” whose objective is to link plant genotypes to plant phenotypes [26-29], particularly in the service of plant biotechnology and crop improvement [30,31].

Analyzing large numbers of plant root system images requires the use of software tools specifically designed for the high-throughput estimation of RSA traits. A number of different tools are available for the analysis of RSA. For example, some programs are specialized for the analysis of images from a specific apparatus, e.g., images from minirhizotrons [32-35]. Next, general tools are available for in-depth analysis of individual monocot root systems regardless of apparatus. These tools rely on significant user input for processing although they can be used in a batch mode [22,36,37]. Similarly, general purpose image processing programs such as ImageJ may be flexible enough to perform many specialized tasks [38]. In practice, end-users often utilize point-and-click approaches that are not scaleable to large numbers of edges within a root system (for an exception, see the recently released software “SmartRoot” which is designed for semi-automated analysis of the hierarchical structure of root systems [39]). Finally, software is also available that can accurately characterize the simpler dicot root system of Arabidopsis[23,24], in a high-throughput fashion [40], but that is not yet suitable for studying the intricate monocot root systems of rice and maize.

GiA Roots is a software tool that estimates RSA traits from a large number of root system images. The main distinguishing characteristics of GiA Roots are that it is both specifically designed for high-throughput analysis and it is extensible. A typical user begins by interactively setting parameters to enable GiA Roots to identify roots from the background, i.e., segmenting the image. Next, the user selects traits of interest for measurement. Finally, the user instructs GiA Roots to automatically estimate traits for images. The GiA Roots pipeline of steps can be executed from the GUI or command line tool as part of a fully automated workflow. A technical user can add new traits to the GiA Roots pipeline using the application programming interface of GiA Roots (GiA-API).

In this manuscript, we describe the implementation of the software, including the main computational steps. We discuss the use of GiA Roots, including validating its measurements via estimation of 19 RSA traits on a large-scale data set of 2393 images from rice root systems. We demonstrate the segmentation capabilities of GiA Roots on root system images. We also demonstrate how GiA Roots can be extended to include a novel trait presented in a recent RSA analysis of 3D reconstructions of root networks [22] but not contained in our earlier study of 2D images of root networks [21]. Finally, we discuss the limitations of the software and plans for further development.

Implementation

GiA Roots is implemented in C++ and is comprised of two applications: a GUI and a command-line tool. The core component of GiA Roots is a trait estimation pipeline. The trait estimation pipeline includes: an optional user-assisted processing of images, scale calibration, trait selection, image segmentation (i.e. the separation of root pixels from background), trait measurement, and output. The GUI and command-line tool can be used together. The GUI workflow provides an interactive means to utilize the trait estimation pipeline of the command line tool (see Figure 1). Hence, the GiA Roots GUI helps users fine-tune the performance of GiA Roots’s image processing steps. The command-line tool facilitates the high throughput processing of large data sets.

thumbnailFigure 1. Annotated snapshot of the GiA Roots GUI. The GiA Roots GUI is a standalone application that provides a user interface to manage input images, define a processing pipeline and manage output. It features a main window with several task windows that can be accessed through the left sidebar. Sequential access to each window page accomplishes major tasks: managing data; selecting traits to measure; tweaking parameters; performing processing; reviewing the results. This linear design helps users to keep track of progress and proceed intuitively with processing of the data. Main computational steps are described in greater detail in Implementation.

Trait estimation workflow

Importing imaging data

A user starts GiA Roots by specifying a set of root images to load. Two different import modes allow users to import either (i) sets of images from one folder; or (ii) the whole hierarchy of folders that contain images. Images are automatically organized into stacks grouped by the folder they were imported from. It is possible to browse all imported images, or filter only a subset of images (i.e., “image-stack”), based on certain search criteria.

Image pre-processing and scale calibration

GiA Roots can apply optional pre-processing transformations to images. The transformations are cropping, rotation and scale calibration. Transformations can be applied on a per-image, per-image stack and an all-images basis. Cropping and rotation are standard image transformations. Scale calibration allows users to set the scale in terms of pixels/cm. This functionality is accessible from the detailed view of any image in the project. Users can either draw a line that corresponds to a specified distance of the scale bar (in cm) or enter the scale value. The scale value can be applied to the current image, current image stack, or all images in the project. GiA Roots retains dimensional information when processing, so that trait estimates will be reported in the proper units, such as cm, cm2 or cm3. When a scale is not specified for an image, its traits will not be scaled and will be presented in units of pixels.

Trait selection

Users can select traits to be estimated from a list of 19 established RSA traits. Illustrated descriptions are provided in the GUI for each trait to detail the relationship between the trait and the structure of the plant. The traits included in this distribution of GiA Roots are based upon two sources [21,22], and include the following: aspect ratio, average width of roots, bushiness, convex area, network depth, network length distribution, major ellipse axis, maximum number of roots, maximum width of root system, median number of roots, minimum ellipse axis, network area, network length, network solidity, network surface area, network volume, network width to depth ratio, perimeter, and specific root length. Specifications for these traits can be found in Table 1 and more information can be found in [21,22]. Additionally, each intermediate processing step can also be selected for output.

Table 1. Traits description

Image processing via segmentation – extracting roots from background

Trait estimation in GiA Roots is done on “segmented” images, in which root pixels are separated from background pixels. Segmentation can be applied on a per-image, per-image stack and an all-images basis. GiA Roots offers three thresholding methods to segment the image into foreground (root) and background. The methods are: global thresholding, adaptive thresholding, and double adaptive thresholding. Global thresholding identifies all roots above a critical pixel intensity as part of the root and all other pixels as background. Adaptive thresholding divides the image into small square regions, calculates the mean pixel intensity in each region, and then decides whether a pixel is part of the root based on whether its intensity exceeds that of the mean intensity in each region by a user-defined offset. Finally, double adaptive thresholding evaluates whether or not a given pixel is part of the root by growing a square centered around the pixel and checking if mean intensity in the square decreases rapidly before the square reaches a critical size. This decrease is common in cases where a group of root pixels are brighter than the background. In all cases, we also provide an additional step in which small connected groups of pixels are removed from the segmented image.

Each thresholding algorithm has default parameter values that facilitate segmentation of foreground and background (see Table 2 and additional information in the GiA Roots manual). Multiple algorithms are provided given the commonly noted inadequacy of any single thresholding algorithm to serve as a universal best choice depending on the application [41,42]. For non-technical users, GiA Roots provides an interactive means to facilitate rapid review and improvement of segmentation. Sliders are available in the GUI in which users can modify thresholding parameters and view the result of segmentation immediately. Hence, users can tweak parameters until they are satisfied with the expected output. Users can save thresholding configurations for later use in the GUI or in the command line tool. GiA Roots saves all configurations along with intermediate images and intermediate numerical results to a chosen folder. GiA Roots therefore enables reproducibility of results with known configuration sets.

Table 2. Segmentation parameters

Trait estimation

After thresholding parameters are specified, the user directs GiA Roots to compute traits from a set of specified images. GiA Roots computes traits directly from the image mask or by computing properties of the skeleton of the image mask (see Table 1). We have chosen the medial axis transformation as a skeletonization method; this transformation is equivalent to a morphological thinning on a binary mask that preserves the topology of the original image [43]. Different skeletonization algorithms can be implemented into GiA Roots via the GiA-API. The GUI workflow follows the processing pipeline of the GiA Roots command line tool. Hence, there is no user involvement in trait computation. In the GUI, partial results of trait estimation can be accessed while computing the RSA traits of all selected images. Preliminary checks on image stacks can be used to test if a given set of configuration parameters yields satisfactory output. In practice, we recommend that users save configuration files in the GUI when analyzing a large number of images. Subsequently, a user may open a saved configuration file which speeds up the workflow prior to executing the computation. Configuration files can also be utilized as part of the command line user interface.

Output

There are two methods to export measurements obtained by GiA Roots. First, numerical measurements can be exported in the form of a comma-separated value (.csv) file suitable for analysis in spreadsheet or other statistical analysis package. Second, image data of all processing steps can be exported to verify results. For example, segmented root images or root image skeletons can be selected for export. All processing intermediates can also be selected for export to verify the performance of GiA Roots.

Automation and integration

GiA Roots can handle small to large data sets, from a single image to many thousands of images. The GUI is particularly well-suited for identifying processing parameters that facilitate estimation of root network traits. Large-scale analysis can also be conducted using the GUI. The command-line tool provides an automated means for processing large numbers of images given configuration parameters (e.g., determined using the GUI). Specifications are contained in an XML file (see Manual). The XML file allows for control of each step of the processing pipeline from algorithm parameters to selected traits to compute. Practically, this enables the GiA Roots pipeline to be integrated with 3rd party scripts or applications. GiA Roots also exposes the intermediate output of all processing steps. This allows users to redirect parts of the GiA Roots pipeline into their workflows. For example, a user could use the output of the segmentation step as input into a 3D reconstruction application. It is important to note that both the GUI and command-line tool output a project XML file that contains enough information to reproduce any output or intermediate outputs.

In practice, the integration of GUI and command-line tool works as follows. Consider a scenario in which users take images of root systems using the same imaging platform and setup. First, the user runs the GiA Roots GUI to fine-tune the image processing steps for a sample of images and determines which RSA traits to measure. Once confirmed, these settings can be exported from the GUI as an XML configuration file. RSA traits of the entire image set can be estimated using the GUI or via the command line tool using the saved XML configuration file. When new data is available, the user executes the command-line tool with the exported configuration file and a job file describing the new incoming data (see Manual for more details). The execution of the command-line tool will be non-interactive (except for displaying a progress log). The command line tool produces the same results as the GUI. Furthermore, since the tool is non-interactive, multiple instances of it can be run on multiple servers. Each instance can also be configured to exploit the parallel architecture of a single machine using multi-threading.

Extending GiA Roots

Technical users may be interested in extending GiA Roots so that it can estimate novel traits. Such extensions are made possible by the general processing scheme in GiA Roots (see Figure 2). In GiA Roots, there are a pre-defined set of data types (e.g., raw image, cropped images, transformed image, etc.) and a set of transformations from one data type into another. Each of the transformations is realized as a software plugin. However, multiple algorithms can be associated with each transformation. For example, a “segmentation” plugin is necessary in order to transform a cropped image into a segmented image. GiA Roots is pre-bundled with three different segmentation plugins: global thresholding, adaptive thresholding, and double adaptive thresholding (see prior text and GiA Roots manual for more details). GiA Roots enables users to specify which plugin to utilize in either the GUI or the command-line tool. Further, users can specify parameters specific to each plug-in to optimize the performance for their image cases.

thumbnailFigure 2. GiA Roots processing chart. Data types are enclosed in ellipses, interactions are enclosed in rectangles. Interactions are realized by plugins, and have several variants. They can also be configured.

Importantly, all transformations are independent of one another. This independence is retained by a software framework design in which data types are the only ways that plugins can communicate. Hence, developers can add a plugin to function as an alternative algorithm to a pre-existing transformation (e.g., segmentation) or trait estimation step (e.g., calculating the total length of the network). Of course, the outputs of algorithms do interact with each other, so modification of a segmentation algorithm that leads to a different segmented image may lead to a different image skeleton. In developing a new feature, developers must ensure that their plug-ins fit into a specific location in a processing pipeline and strictly comply with uniqueness requirements for both the input and output data type. Developers must specify the dependency of the plugin based on the data types it requires and the data types it outputs. Description of the APIs for adding new features can be found in the manual (see project homepage).

Results and discussion

GiA Roots can estimate RSA traits for thousands (and more) root images

Here, we demonstrate that the output of trait estimation algorithms in GiA Roots is equivalent to a previously validated set of algorithms for estimating RSA traits [21]. In the previous study [21], RSA traits of a manufactured “root” wire model were estimated using a MATLAB-based processing pipeline. The estimates were shown to be in close agreement with hand measurements (0.5% to 6.7%differences) [21]. GiA Roots leverages the OpenCV libary to estimate many of the same traits (in addition to novel traits) as estimated in the prior MATLAB-based pipeline. Hence, we anticipate that estimates of traits in GiA Roots should coincide with prior estimates. We compared GiA Roots against the output of the prior MATLAB-based processing pipeline given a set of 2393 segmented rice root images from 12 genotypes (see http://www.rootnet.biology.gatech.edu/data/PlantPhysData.zip webcite for the complete set of images). In Figure 3 we show the correspondence between the MATLAB-based output (x-axis) and the GiA Roots-based output (y-axis) for two traits where each point in a panel denotes a given image. The two traits we include are median number of roots and network volume. All comparisons can be found in Additional file 1. The results correlated as expected with an R2≈1.0. Note that small differences in trait estimates can arise due to implementation differences in the medial axis transformation. There are two caveats involving refinement of trait definitions between the study of [21] and the current study. Previously, network length distribution was defined as the ratio of “the total root length in the upper one-third of the root depth to the total root length located in the lower two-thirds of the root depth” [21]. Whereas here, we re-define this trait as the fraction of network pixels found in the lower 2/3 of the network, where the lower 2/3 of the network is defined based on the network depth (see Table 1 for trait descriptors). In addition, we previously reported the thickness of roots as the average root radius [21], whereas now we define it as the average root diameter (hence there is a factor of 2 difference which we correct in evaluating the correspondence). Note that GiA Roots also includes three novel traits not contained in the previous study [21], these are properties of the best-fit ellipse: (i) minor ellipse axis length; (ii) major ellipse axis length; and (iii) ellipse axis ratio. We validated our implementation of these features in GiA Roots against a MATLAB implementation. All ellipse-related traits are highly correlated with R2≈1 as before. The comparisons are included in Additional file 1. Processing of the dataset from [21] took 2 hours and 28 minutes for 2393 images on a MacBook Pro, 2.2 Ghz Intel i7, 4GB DDR3 RAM, with Mac OSX 10.6.8. Runtime performance of GiA Roots is comparable to the MATLAB implementation. However, the ability to interactively adjust parameters is available in GiA Roots and is not available in the MATLAB implementation.

Additional file 1. Comparison of GiA Roots to previous benchmark. Compilation of statistical correlation of the output of GiA Roots against a previous benchmark code written in Matlab [21]. The comparisons are for 16 different traits estimated from 2393 previously thresholded rice root images taken from 12 genotypes. All 2393 images are available at http://www.rootnet.biology.gatech.edu/data/PlantPhysData.zip webcite.

Format: PDF Size: 555KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 3. Comparison of GiA Roots estimation of trait values against a prior benchmark. The two plots represent a comparison of GiA Roots estimation of trait values compared against a previously validated set of algorithms [21]. Each point represents a trait estimate from one of 2393 images. Note that the median number of roots is defined to be an integer so nearly all comparisons in the left panel exactly coincide, but when plotted they appear to give rise to a ‘gridded’ pattern. The R2values confirm the strong correspondence of the two implementations of the same trait.

GiA Roots identifies roots from noisy images

GiA Roots utilizes a thresholding algorithm to segment foreground (roots) from background. The three algorithms available in GiA Roots are global thresholding (GT), adaptive thresholding (AT), and double-adaptive thresholding (DAT). As stated before, the choice of the thresholding algorithm is highly dependent on the provided image quality [41,42]. Here, we visually compare the thresholding algorithms provided in GiA Roots with two root images (see Figure 4). The two images are selected from the publicly available data set in [21] and correspond to individuals on the 14th day after planting from the rice genotypes: Jefferson (top) and Teqing (bottom). We compare the default and manually optimized parameters for each of the two images as summarized in Table 3. The annotations in Figure 4 illustrate differences between the default parameters and those improved via user-assisted optimization.

thumbnailFigure 4. Comparison of the three available thresholding algorithms in GiA Roots. The three algorithms are global thresholding (GT), adaptive thresholding (AT), and double adaptive thresholding (DAT). Each algorithm is demonstrated on two relevant examples with default parameter settings and manually optimized parameter settings. Annotations 1-11 point at significant differences which are detailed in the text.

Table 3. Parameters used for comparison of thresholding

In the first example, annotation (1) shows improved root separation of adaptive thresholding (AT) in comparison to global thresholding (GT). Annotation (2) shows a tendency of double adaptive thresholding (DAT) to omit dense clusters of roots; in the given example this omission occurs in the upper root crown. The optimized parameter set shows additional crown roots (3) achieved with GT and even additional lateral roots (4) with AT. Although lateral roots are less resolved for DAT compared to AT, the root networks exhibit less noisy edges (5), which improves estimates of traits that rely on skeleton extraction.

In the second example, annotation (6) shows a major drawback of GT. The dense part of the root crown is not resolved into single roots. By comparison, AT provides a clear separation of the individual roots (7) and DAT shows the loss of root structure (8) for the default parameter set. The manually optimized parameters again show additional fine roots in the upper root crown and slightly better root separation (9). Annotations (10) and (11) show comparable results of AT and DAT. Both AT and DAT were able to resolve lateral roots. Although the default parameters may be suitable for many use cases, we recommend the modification of thresholding parameters to suit specific needs.

Custom traits can be integrated into the GiA Roots pipeline

GiA Roots can be extended by copying compiled modules and accompanying information into the GiA Roots folder. Hence, GiA Roots does not require recompilation or any change of the software processing pipeline to include new traits or algorithms. In order to demonstrate this feature, we show how GiA Roots can be augmented to analyze 3D reconstructions of root networks. Although these features are still in development, we intend to release them in future versions of the software. As can be seen in Figure 5, a set of new traits was added to the software. The GUI of GiA Roots will automatically detect the new traits and a new checkbox will appear in the trait selection panel.

thumbnailFigure 5. Schematic of GiA Roots extended with new trait estimation algorithms. A set of new traits is highlighted in the “Features & Algorithms” panel of GiA Roots. Installing new plugins is accomplished by copying the compiled plugin and documentation into the GiA Roots folder.

The extension of GiA Roots to include 3D trait estimation was accomplished as follows: (i) A new series of algorithms were written in C++ following the protocol for APIs (see Manual). (ii) The algorithm was compiled as a library. (iii) The compiled library was copied into the library directory of the GiA Roots directory structure. (iv) Visual depiction and textual description of the trait was added into the info directory including both an .html file description of the trait and a .png image file for display in the GUI – this step is not necessary, but it is helpful and encouraged. (v) GiA Roots was run, without recompilation, and the new trait functionality and its description are available to the user.

Caveats and limitations

The objective of GiA Roots is to improve the high-throughput estimation of RSA traits from root network images. Hence, the caveats and limitations of the software reflect the design choices associated with meeting this objective. First, GiA Roots does not include the ability to interactively measure properties of individual networks (e.g., distance between two arbitrary points on the root network). Hence, if image quality varies significantly from image to image, then a more interactive software that allows greater user-interaction with images is likely to be more appropriate. Next, the quality of trait estimation in GiA Roots depends on the quality of segmentation. Our solution is to have the user select amongst thresholding algorithms with pre-defined parameter settings that we have found to be efficient in extracting root networks from background. These algorithms work whether roots are light on a dark background or vice-versa. However, parameters can also be changed by the user which are then store for full reproducibility of output. Care should be taken by users to validate the performance of these segmentation algorithms on their data. Finally, GiA Roots is extensible like some, but not all other software tools (see Table 4). However, the API is currently only available for libraries written in C++. Hence, developers require specific expertise in order to add new functionality.

Table 4. Comparison of root system analysis software

Future plans for development

GiA Roots is a general framework for analyzing the shapes of root systems. We anticipate extending GiA Roots to address problems that are direct extensions of the current software as well as explore possible applications to closely related problem domains.

First, a number of technological innovations now permit the 3D reconstruction of root systems. For example, an imaging platform and software tool for 3D reconstruction of rice root systems was recently introduced that utilizes multiple images of the same root system grown in a transparent gel [22]. As we have already shown, it is possible to extend GiA Roots to include traits that operate directly on 3D voxel data rather than sets of 2D images. We plan to build a robust set of trait estimation tools for 3D voxel data as part of a future version.

Second, root systems are dynamic, they grow heterogeneously in space and in time. Hence, we envision linking images of root systems (whether in 2D or 3D) in order to track the growth of the entire root system or components thereof. A similar approach has been espoused using point and click approaches [22], however we will attempt to fully automate the root tracking problem.

Third, we will explore how GiA Roots can be extended to analyze other plant phenomics problems. For example, a subset of us (CAP, YM, OS and JSW) collaborated on the development of LEAF GUI [44], a user-assisted software for the extraction of leaf venation networks. LEAF GUI enables users to utilize a set of transformations in a non-linear fashion to clean images prior to analysis. Hence, it is ideally suited for analyzing archival image data of varying quality. However, large-scale analysis of hundreds, thousands or more leaf images may benefit from matching the capabilities of interactive software like LEAF GUI with the more restrictive pipeline of GiA Roots. We note that other software is already available to characterize leaf shape [45], and so we emphasize that our efforts will focus on physical plant networks for which many of the tools of GiA Roots will apply.

Conclusions

GiA Roots is a software framework for the high-throughput analysis of root system architecture. GiA Roots can analyze a single image up to many thousands of images, and for each image, extract the root network, estimate its traits, and report quantitative trait estimates and processing intermediates back to the end-user. GiA Roots is designed for end-users with limited technical background and for developers who wish to integrate a root trait estimation framework into their workflows. As such, GiA Roots is comprised of two tools: a GUI and a command line tool. The two tools perform all trait estimation steps using the same image processing pipeline. We anticipate that most users will utilize the GUI exclusively, whereas technical users will also utilize the command line tool. Importantly, GiA Roots is extensible unlike some root trait estimation packages currently available. Hence, developers can add on new algorithms and trait estimation steps using plugins and have these plugins interface with GiA Roots using documented APIs. This extensible step will help shorten the development time between identification of a possibly important root trait and its inclusion in processing pipelines.

The ability to link genotype to phenotype in the service of plant biotechnology will require the contributions of many groups and the utilization of many molecular, physiological and imaging techniques. We envision that GiA Roots can be utilized to help link genotypic with phenotypic information by improving the pipeline for extracting RSA traits. The trait algorithms currently bundled in GiA Roots are designed to quantify the RSA of monocot plants. However, in the future we anticipate adding modules that can characterize RSA for dicots, quantify temporal growth, and estimate traits from 3D reconstructions.

Availability and requirements

Project name: GiA Roots

Project home page:http://www.giaroots.org webcite

Accessibility Academic users are encouraged to includetheir contact information before downloading thesoftware. Commercial users should contact JSW(http://jsweitz@gatech.edu webcite) for more information on obtainingthe software.

Operating system(s): Binary distributions available forMicrosoft Windows 7 and Apple Mac OS X 10.6.x.

Programming language: C++

Frameworks and libraries used:

· Qt framework (http://qt.nokia.com/products webcite) - a cross-platform application and UI framework. Copyright (c) 2008-2010 Nokia Corporation and/or its subsidiaries. Licensed under LGPL v2.1.

· OpenCV (http://opencv.willowgarage.com/wiki/ webcite) - a library of programming functions for real time computer vision. BSD License.

· tinyxml (http://sourceforge.net/projects/tinyxml webcite) - a simple, small, minimal, C++ XML parser. Licensed under zlib.

License: Academic use is governed by the followinglicense. For commercial licensing contacthttp://jsweitz@gatech.edu webcite.

GiA Roots - “Software for the High Throughput Analysis of Plant Root System Images”. Copyright (C) 2010-2012 Georgia Tech Research Corporation and Duke University.By downloading and/or utilizing this Program, you agree to become bound by the terms and conditions of this license. If you do not agree with the terms and conditions setforth below, do not use this Program or any portion thereof in any form or manner.

This Program is licensed, not sold to you (“User”), by Georgia Tech Research Corporation (“GTRC”), owner of all rights and title to the code and accompanying documentation (“Program”) for use under the terms and conditions of this license. As such, GTRC reserves any and all rights not expressly granted to User under this license. With respect to the Program to which GTRC has exclusionary rights, GTRC hereby grants to User a nontransferable, non-exclusive license to use the Program for User’s own educational and non-commercial research purposes only. GTRC shall have the right to terminate this agreement and/or license without cause at any time.

User accepts the Program on an “as is” basis. GTRC makes no representation or warranty that the Program will be accessible or downloadable. GTRC makes no warranty that all errors can be or have been eliminated from the Program. GTRC shall not be responsible for losses of any kind resulting from the use of the Program and its accompanying document, and can in no way provide compensation for any losses sustained, including but not limited to any obligation, liability, right, claim or remedy for tort, or for any actual or alleged infringement of patents, copyrights, trade secrets, or similar rights of third parties, nor any business expense, machine downtime or damages caused User by any deficiency, defect or error in the Program or malfunction thereof, nor any incidental or consequential damages, however caused. GTRC disclaims all warranties, both express and implied respecting the use and operation of the Program and its accompanying documentation, including all implied warranties of merchantability and fitness for particular purpose and any implied warranty arising from course of performance, course of dealing or usage of trade. The User of the Program is expected to make the final evaluation of the Program’s usefulness in User’s own environment.

Abbreviations

GiA Roots, General Image analysis of Roots; API, Application programming interface; RSA, Root system architecture.

Competing interests

PNB is CEO and co-founder of GrassRoots Biotechnology, which uses the gel platform as part of its root phenotyping efforts. JH and CNT consulted for GrassRoots Biotechnology. The remaining authors declare that they have no competing interests.

Authors’ contributions

TG was the primary developer and coder of GiA Roots, co-designed the GiA Roots software framework, and assisted in writing the manuscript. YM co-designed the GiA Roots framework. AB performed analysis of GiA Roots on data and assisted with software implementation. OS designed and developed algorithms. CAP assisted with algorithm design and GUI requirements. BM assisted with user-interface design and usability for both the GUI and command line interface. PZ assisted with the GUI. ASIP contributed root images and contributed to user-interface and algorithm design. CT contributed to user-interface and algorithm design. SF contributed to user-interface and algorithm design. JH contributed to algorithm design. PNB contributed root images and contributed to algorithm design. JSW developed the project idea, designed and developed algorithms, and contributed to software design. JSW wrote the manuscript with contributions from TG, YM, AB, OS, CAP, CT, ASIP, and PNB. All authors read and approved the final manuscript.

Acknowledgements

The authors thank R. Clark, L. Kochian, J. Lynch, and T. Schaprian for feedback on the software, T. Hilley for technical assistance with the release of the software, and three anonymous reviewers for feedback on the manuscript. This work was supported by a grant from the National Science Foundation (DBI-0820624-02 to JH, PNB and JSW) and from the Defense Advanced Research Projects Agency (HR0011-09-1-0055 to JH, PNB and JSW). Joshua S. Weitz, Ph.D., holds a Career Award at the Scientific Interface from the Burroughs Wellcome Fund. SF is supported by National Science Foundation of China under grant 30890131.

References

  1. Lynch J: Root architecture and plant productivity.

    Plant Physiol 1995, 109(1):7-13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Casper BB, Jackson RB: Plant competition underground.

    Ann Rev Ecol Syst 1997, 28:545-570. Publisher Full Text OpenURL

  3. Hodge A, Robinson D, Griffiths B, Fitter A: Why plants bother: root proliferation results in increased nitrogen capture from an organic patch when two grasses compete.

    Plant Cell and Dev 1999, 22(7):811-820. OpenURL

  4. Ennos AR: The scaling of root anchorage.

    J Theor Biol 1993, 161(1):61-75. Publisher Full Text OpenURL

  5. Niklas KJ, Spatz HC: Allometric theory and the mechanical stability of large trees: Proof and conjecture.

    Am J Bot 2006, 93(6):824-828. PubMed Abstract | Publisher Full Text OpenURL

  6. Van Der Putten WH: Plant defense belowground and spatiotemporal processes in natural vegetation.

    Ecology 2003, 84(9):2269-2280. Publisher Full Text OpenURL

  7. Bertin C, Yang XH, Weston LA: The role of root exudates and allelochemicals in the rhizosphere.

    Plant and Soil 2003, 256(1):67-83. OpenURL

  8. Lynch JM, Whipps JM: Substrate flow in the rhizosphere.

    Plant and Soil 1990, 129(1):1-10. Publisher Full Text OpenURL

  9. Freiberg C, Fellay R, Bairoch A, Broughton W, Rosenthal A, Perret X: Molecular basis of symbiosis between Rhizobium and legumes.

    Nature 1997, 387(6631):394-401. PubMed Abstract | Publisher Full Text OpenURL

  10. Waisel Y, Eshel A, Kafkafi, U (Eds): Plant Roots: The Hidden Half. Marcel Dekker Inc, New York, Basel; 2002. OpenURL

  11. Trachsel S, Kaeppler SM, Brown KM, Lynch JP: Shovelomics: high throughput phenotyping of maize (Zea mays L.) root architecture in the field.

    Plant and Soil 2011, 341(1-2):75-87. Publisher Full Text OpenURL

  12. Heeraman DA, Hopmans JW, Clausnitzer V: Three dimensional imaging of plant roots in situ with x-ray computed tomography.

    Plant and Soil 1997, 189(2):167-179. OpenURL

  13. Perret JS, Al-Belushi ME, Deadman M: Non-destructive visualization and quantification of roots using computed tomography.

    Soil Biol & Biochem 2007, 39(2):391-399. PubMed Abstract | Publisher Full Text OpenURL

  14. Gregory PJ, Hutchison DJ, Read DB, Jenneson PM, Gilboy WB, Morton EJ: Non-invasive imaging of roots with high resolution X-ray micro-tomography.

    Plant and Soil 2003, 255(1):351-359. OpenURL

  15. Tracy SR, Roberts JA, Black CR, McNeill A, Davidson R, Mooney SJ: The X-factor: visualizing undisturbed root architecture in soils using X-ray computed tomography.

    J Exp Bot 2010, 61(2):311-313. PubMed Abstract | Publisher Full Text OpenURL

  16. van der Weerd L, Claessens MMAE, Ruttink T, Vergeldt FJ, Schaafsma TJ, Van As H: Quantitative NMR microscopy of osmotic stress responses in maize and pearl millet.

    J Exp Bot 2001, 52(365):2333-2343. PubMed Abstract | Publisher Full Text OpenURL

  17. Van As H: Intact plant MRI for the study of cell water relations, membrane permeability, cell-to-cell and long distance water transport.

    J Exp Bot 2007, 58(4):743-756. PubMed Abstract | Publisher Full Text OpenURL

  18. Van As H, Scheenen T, Vergeldt FJ: MRI of intact plants.

    Photosynth Res 2009, 102(2-3, Sp. Iss. SI):213-222. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Fang S, Yan X, Liao H: 3D reconstruction and dynamic modeling of root architecture in situ and its application to crop phosphorus research.

    Plant J 2009, 60(6):1096-1108. PubMed Abstract | Publisher Full Text OpenURL

  20. Yazdanbakhsh N, Fisahn J: High throughput phenotyping of root growth dynamics, lateral root formation, root architecture and root hair development enabled by PlaRoM.

    Funct Plant Biol 2009, 36(10-11, Sp. Iss. SI):938-946. OpenURL

  21. Iyer-Pascuzzi AS, Symonova O, Mileyko Y, Hao Y, Belcher H, Harer J, Weitz JS, Benfey PN: Imaging and analysis platform for automatic phenotyping and trait ranking of plant root systems.

    Plant Physiol 2010, 152(3):1148-1157. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Clark RT, MacCurdy RB, Jung JK, Shaff JE, McCouch SR, Aneshansley DJ, Kochian LV: Three-dimensional root phenotyping with a novel imaging and software platform.

    Plant Physiol 2011, 156(2):455-465. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. French A, Ubeda-Tomas S, Holman TJ, Bennett MJ, Pridmore T: High-throughput quantification of root growth using a novel image-analysis tool.

    Plant Physiol 2009, 150(4):1784-1795. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Armengaud P, Zambaux K, Hills A, Sulpice R, Pattison RJ, Blatt MR, Amtmann A: EZ-Rhizo: integrated software for the fast and accurate measurement of root system architecture.

    Plant J 2009, 57(5):945-956. PubMed Abstract | Publisher Full Text OpenURL

  25. Hund A, Trachsel S, Stamp P: Growth of axile and lateral roots of maize: I development of a phenotying platform.

    Plant and Soil 2009, 325(1-2, Sp. Iss. SI):335-349. Publisher Full Text OpenURL

  26. Furbank RT: Plant phenomics: from gene to form and function.

    Funct Plant Biol 2009, 36(10-11, Sp. Iss. SI):V-VI. OpenURL

  27. Sozzani R, Benfey PN: High-throughput phenotyping of multicellular organisms: finding the link between genotype and phenotype.

    Genome Biol 2011, 12(3):219. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Finkel E: Imaging with ‘phenomics,’ plant scientists hope to shift breeding into overdrive.

    Science 2009, 325(5939):380-381. PubMed Abstract | Publisher Full Text OpenURL

  29. Zhu J, Ingram PA, Benfey PN, Elich T: From lab to field, new approaches to phenotyping root system architecture.

    Curr Opin Plant Biol 2011, 14:310-317. PubMed Abstract | Publisher Full Text OpenURL

  30. de Dorlodot S, Forster B, Pages L, Price A, Tuberosa R, Draye X: Root system architecture: opportunities and constraints for genetic improvement of crops.

    Trends Plant Sci 2007, 12(10):474-481. PubMed Abstract | Publisher Full Text OpenURL

  31. Tester M, Langridge P: Breeding technologies to increase crop production in a changing world.

    Science 2010, 327(5967):818-822. PubMed Abstract | Publisher Full Text OpenURL

  32. Zeng G, Birchfield ST, Wells CE: Automatic discrimination of fine roots in minirhizotron images.

    New Phytologist 2008, 177(2):549-557. PubMed Abstract OpenURL

  33. Rootfly: Software for Minirhizotron Image Analysis Http://www.ces.clemson.edu/stb/rootfly/ webcite.

  34. RootView Http://www.mv.helsinki.fi/aphalo/RootView.html webcite.

  35. WinRHIZO Tron Http://www.regent.qc.ca/products/rhizo/RHIZOTron.html webcite.

  36. Le Bot J, Serra V, Fabre J, Draye X, Adamowicz S, Pages L: DART: a software to analyse root system architecture and development from captured images.

    Plant and Soil 2010, 326(1-2, Sp. Iss. SI):261-273. Publisher Full Text OpenURL

  37. WinRHIZO Http://www.regent.qc.ca/products/rhizo/WinRHIZO.html webcite.

  38. Rasband WS: ImageJ, US National Institutes of Health, Bethesda MD USA. http://imagej.nih.gov/ij/ webcite.

  39. Lobet G, Pagés L, Draye X: A novel image-analysis toolbox enabling quantitative analysis of root system architecture.

    Plant Physiol 2011, 157:29-39. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Naeem A, French AP, Wells DM, Pridmore TP: High-throughput feature counting and measurement of roots.

    Bioinformatics 2011, 27(9):1337-1338. PubMed Abstract | Publisher Full Text OpenURL

  41. Pal NR, Pal SK: A review on image segmentation techniques.

    Pattern Recognit 1993, 26:1277-1294. Publisher Full Text OpenURL

  42. Zhang YJ: A survey on evaluation methods for image segmentation.

    Pattern Recognit 1996, 29:1335-1346. Publisher Full Text OpenURL

  43. Serra J: Image analysis and mathematical morphology. Academic Press, London; 1982. OpenURL

  44. Price CA, Symonova O, Mileyko Y, Hilley T, Weitz JS: Leaf extraction and analysis framework graphical user interface: segmenting and analyzing the structure of leaf veins and areoles.

    Plant Physiol 2011, 155(1):236-245. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Bylesjoe M, Segura V, Soolanayakanahally RY, Rae AM, Trygg J, Gustafsson P, Jansson S, Street NR: LAMINA: a tool for rapid quantification of leaf size and shape parameters.

    BMC Plant Biol 2008, 8:82. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL