Skip to main content

RADIX: rhizoslide platform allowing high throughput digital image analysis of root system expansion

Abstract

Background

Phenotyping of genotype-by-environment interactions in the root-zone is of major importance for crop improvement as the spatial distribution of a plant’s root system is crucial for a plant to access water and nutrient resources of the soil. However, so far it is unclear to what extent genetic variations in root system responses to spatially varying soil resources can be utilized for breeding applications. Among others, one limiting factor is the absence of phenotyping platforms allowing the analysis of such interactions.

Results

We developed a system that is able to (a) monitor root and shoot growth synchronously, (b) investigate their dynamic responses and (c) analyse the effect of heterogeneous N distribution to parts of the root system in a split-nutrient setup with a throughput (200 individual maize plants at once) sufficient for mapping of quantitative trait loci or for screens of multiple environmental factors. In a test trial, 24 maize genotypes were grown under split nitrogen conditions and the response of shoot and root growth was investigated. An almost double elongation rate of crown and lateral roots was observed under high N for all genotypes. The intensity of genotype-specific responses varied strongly. For example, elongation of crown roots differed almost two times between the fastest and slowest growing genotype. A stronger selective root placement in the high-N compartment was related to an increased shoot development indicating that early vigour might be related to a more intense foraging behaviour.

Conclusion

To our knowledge, RADIX is the only system currently existing which allows studying the differential response of crown roots to split-nutrient application to quantify foraging behaviour in genome mapping or selection experiments. In doing so, changes in root and shoot development and the connection to plant performance can be investigated.

Background

The spatial distribution of a plant’s root system is crucial for the success of the plant in a given environment because it determines how easily a plant can access water and nutrient resources of the soil. Root system architecture (RSA) is governed by the environment, but also by the genotype. Today it is unclear to what extent genetic variations in root system responses to spatially varying soil resources can be utilized. Closing this knowledge gap is certainly one step forward towards improving resource-use efficient crops. First studies provided evidence for an advantage of certain root architectural characteristics, such as deep rooting under drought [1, 2], shallow rooting under low phosphorus availability [3] or steep rooting angles under low nitrogen availability [47].

One major problem studying root growth is the fact that roots are often hidden in their growth matrix like soil unless the growth medium is transparent such as a gel. In the field, even the most advanced procedures offer only a partial glimpse of RSA. To overcome this limitation, research was done on developing monitoring systems of root growth in greenhouse approaches that provide visibility of the roots. Various phenotyping platforms were developed to monitor root growth non-invasively in soil based systems such as rhizotrons/-boxes containing sand/soil substrates [8], via X-ray micro-tomography [9, 10] or magnetic resonance imaging [11, 12] measurements. However, to date the throughput of these techniques is still not high enough for the assessment of large sets of genotypes as needed for mapping studies of quantitative trait loci (QTL).

For fast and simple screening, germination paper has become state-of-the art in many labs focusing on roots [1316]. The systems using filter paper as substrate are sometimes called “growth pouches” [17]. Rhizoslides are a large-scale refinement of growth pouches, which allow to study root development for a prolonged period of time [18]. The growth pouch system was developed to map QTLs in collections of about 200 genotypes in various populations of maize during the very early seedling development [1922] whereas the rhizoslides allow studying root growth until four-leaf stage.

A limitation of the growth pouches was that only the dynamics of the embryonic roots could be studied, ignoring the most prominent root type of maize, the shoot born crown roots. Albeit embryonic roots and crown roots are under different genetic control [23, 24] and may respond distinctly different to environmental stimuli. For example, Yu et al. [25], reviewed that embryonic roots of maize seedlings respond to nitrogen-rich patches by increasing the length of lateral roots, while crown roots of adult plants increased length and density of lateral roots [26, 27].

It is widely accepted that plants can adjust to localized soil enrichment by physiological changes [28]. Yet, to date there is still a lack of understanding the utility of this response for plant improvement because of a lack of appropriate methods to quantify responses of different genotypes. For breeding research, it would be highly beneficial to select genotypes not only by their final phenotype but also by their dynamic response to temporally changing or spatially varying nutrient availability. For example, Lynch et al. [7] proposed that an “optimal” root system should ignore local resource availability in the top soil and should grow into deep soil layers to be prepared against drought and nitrogen (N) stress later in the season. Thus, a reduced response to spatial variability would be desirable. Others propose the opposite, i.e. that a better response to local variation of nutrients would be beneficial to take advantage of side banded fertilizer placements [27]. A proof-of-concept of the split-nutrient application in rhizoslides was presented by in t’ Zandt et al. [29] observing a strong differential response of crown roots to N. It is conceivable that there is an optimum response to N, which may vary depending on the species, crop management and target environment. With phenotyping systems, such as rhizoslides, it will be possible to test corresponding research hypotheses in this context. One opportunity to evaluate plant performance is to study the shoot development. The measurement of canopy development is an important component of sound root research as the shoot acts as source of carbohydrates and sink for nutrient and water. There is a close link between root and shoot growth in many crops (see [30] for a discussion). Hence, a better understanding of the relationship between root and shoot development is an important need in order to identify heritable root traits for which the additional costs of direct selection are justified as pointed out by Wissuwa et al. [31]. For example, by studying the development of rooting depth in relation to shoot development in a diverse set of maize, Grieder et al. [32] found a strong linear relationship between the two traits under well watered conditions, leaving little genetic variation to alter rooting depth without changing canopy size. The non-destructive measurement of canopy size is well established [33, 34] and the adaptation of linear models is an accepted method [35] whereas for root traits dynamic studies were not often performed and less is known about models fitting this need. Pioneer work was done by Adu et al. [13] modelling the root growth dynamics of Brassica rapa.

With the rhizoslides, a second generation of the growth pouch system was developed [18] and tested for its suitability for split-nutrient application (in t’ Zandt et al. [29]). Here we present the RADIX platform facilitating an improved irrigation, handling and imaging of rhizoslides. The current version contains 100 slides and enables the cultivation of 200 maize plants at the same time, one on each side of the slide. The aim of the study was to demonstrate the value of the platform for phenotyping of natural variation of root traits in maize by investigating (a) root and shoot growth beyond the early seedling stage, and (b) the effect of split-nutrient application on the crown root development. This setup should provide information on genotypic variability of N uptake under spatially heterogeneous N availability.

Methods

Cultivation in the RADIX platform

Plants were grown in so called “rhizoslides” that consist of two PVC bars (600 × 60 × 10 mm) and a plexiglass sheet (650 × 530 × 4 mm) fixed with two screws between the bars (Fig. 1e) [18]. The dimensions of the assembled rhizoslide were 600 × 60 × 38 mm. The plexiglass sheet was covered with wet germination paper (49 × 61 cm) (Anchor Steel Blue Seed Germination Blotter, Anchor Papers Co, USA) on both sides, serving as substrate. These were in turn covered by a transparent, oriented polypropylene (OPP) foil with micro holes of 70 ± 10 µm with a distance of 105 mm (in x and y orientation) to allow for gas exchange (Maag, GmbH, Iserlohn, Germany).

Fig. 1
figure 1

a, b Schematic figures of the RADIX platform with the imaging station moving on rails along a rack holding the rhizoslides. After positioning the station, the rhizoslides can be slid onto the station. Two drip irrigation lines are mounted at the backside of the rack enabling the supply of two different nutrient concentrations. c The imaging station holds two webcams to image the shoot from the side and from the top; a DSLR camera is mounted at 1 m distance from the slide surface to monitor the roots. Roots are illuminated by LED bars (three on each side); illumination is controlled by a microcontroller built and programmed in Arduino 1.0. Rhizoslides on the slide mount can be turned by 180° to image both sides of the rhizoslide. d The dosatron pump system consists of one dosatron supplying base nutrient solution that is split up to two dosatrons allowing for differential nutrient supply. e Scheme of a rhizoslide and the irrigation for the split root setup. High and low N concentrations were dripped into reservoirs from where tubes transported the nutrient solution to the split-nutrient germination paper separated by the wax barrier

For sterilization, the germination paper was heated in three cycles from room temperature to 80 °C and kept at this temperature for at least 120 min. Between the heating periods the paper was kept for 20–22 h in an oven at 37 °C and 50 % relative humidity [36]. To prevent fungal growth during the experiment, a method described by Bohn et al. [37] was used. The germination paper was moistened with deionized water containing 2.5 gL−1 Malvin (Syngenta Agro AG, Dielsdorf, Switzerland) containing the active component Captan. Once a week, each rhizoslide was watered with 15 mL of the Malvin solution to balance for leached fungicide.

The plant material consisted of the core EURoot maize panel (www.euroot.eu), with 25 dent inbred lines B73 (1), EC169 (4), EZ11A (13), F98902 (5), FV353 (6), LH38 (14), MS153 (7), Oh33 (16), EZ37 (10), MS71 (22), Os420 RootABA− (11), Os420 RootABA+ (12), Mo17 (3), FC1890 (18), LAN496 (23), W64A (20), F7028 (8), EZ47 (9), B100 (24), N6 (25), N25 (26), A347 (32), PH207 (33), F1808 (34), B97 (36),) crossed to UH007 as common flint tester supplied by Delley seeds and plants Ltd. (DSP Ltd., Switzerland). The number in brackets indicated the EU_ID.

Seeds were surface sterilized with 0.3 % PrevicurN (Bayer CropScience GmbH, Monheim, Germany) for 30 min and became touch dry without rinsing the seed. Subsequently, seeds were incubated in Petri dishes lined on paper soaked with a soil bacteria mixture (0.0001 % RhizoPlus 42, Andermatt Biocontrol AG, Grossdietwil, Switzerland and 0.01 % FZB24, Bayer AG CropScience, Zollikofen, Switzerland) to promote the development of a healthy rhizosphere and to prevent seed-borne infections. Seeds were kept for 48 h at 26 °C in the dark for germination and were then transferred into the rhizoslides.

For this study, 100 Rhizoslides were prepared with one plant in each rhizoslide. The seed was placed on top of the germination paper enabling simultaneous study of seminal and crown roots (for details, see [24]). A split-nutrient system was created by ironing a vertical water impermeable wax barrier into the germination paper and by a drip irrigation system (for details, see [29]). The roots could freely develop on both germination papers and on both sides of the wax barrier. Each paper was a split nutrient setup; irrigated with a high and low N solution, respectively. All 100 rhizoslides were supplied with a solution with low N (10 % of full strength ~ high N solution) for 14 days (2–3 leaf stage). Subsequently, one side was irrigated with a thousand times lower concentration of N compared to the other side of the wax barrier (Fig. 1e). The split root treatment remained for 14 days until the end of the experiment (3–4 leaf stage). For studying root growth, images were taken every second day.

Composition of the nutrient solutions

To prevent precipitation, two stock solutions were prepared. The solution A was identical for both treatments containing: KH2PO4: 0.5 mmol L−1; MgSO4 × 7 H2O: 2 mmol L−1; MnCl2 × 4 H2O: 9.15 µmol L−1; CuSO4 × 5 H2O: 0.2 µmol L−1; H3BO3: 46.25 µmol L−1; Na2MoO4 × 2 H2O: 0.58 µmol L−1; ZnSO4 × 7 H2O: 0.77 µmol L−1; FeEDTA (C10H12FeN2NaO8 × H2O): 0.04 mmol L−1, (pH ~ 6). The solution B contained: KNO3, Ca(NO3)2 × 4H2O and NH4NO3, (pH ~ 6). The final concentrations in the high N treatment were: KNO3: 4.5 mmol L−1; Ca(NO3)2 × 4H2O: 4 mmol L−1; NH4NO3: 1 mmol L−1. In the low N treatment, the ammonium and nitrate concentrations were reduced. The solution contained: KNO3: 5 µmol L−1; Ca(NO3)2 × 4 H2O: 5 µmol L−1; NH4NO3: 1 µmol L−1and to balance the ion activities additionally KCl: 5 µmol L−1 and CaCl2 × 2H2O: 5 µmol L−1. The solutions were calculated using GEOCHEM-EZ [38]. pH was adjusted by adding KOH.

Construction of the RADIX platform

The rhizoslides were hanging in a rack (500 × 65 × 80 cm) constructed with aluminium profiles (KANYA AG, Rüti, Schweiz) with 100 rhizoslide compartments. Each compartment was 60 cm long and 4 cm wide (Fig. 1a, Additional file 1). The distance between two neighbouring plants was 5 cm. A declining drip pan was built under the rack to collect the dripping nutrient solution (Fig. 1b) and a pump with automatic water level sensors (Gardena Comfort Tauchpumpe 9000 Aquasenor, Gardena Ulm, Germany) was used to pump down the solution. The sides of the rack were covered with a black polyethylene foil to prevent the roots against incidence of light. Each rhizoslide was covered on top with a cover plate containing holes for the shoot.

The irrigation system of every rhizoslide consisted of two PVC tubes (5/3 mm) (GVZ-Gossart AG, Otelfingen, Switzerland) and 25 mL PE tubes (Semadeni AG, Ostermundigen, Switzerland) acting as water/nutrient solution reservoirs. The reservoirs were filled via an irrigation system consisting of three Dosatron pumps (Diaphragm range 2.5 m3 h−1) (Dosatron, International S.A.S, Tresses, France) to mix the stock solutions (100 times higher than final concentration) with deionized water. The stock solutions were homogenized by circulation starting ½ h before the mixing of the solutions. One Dosatron was placed in series with two parallel switched Dosatrons (Fig. 1d).

The first Dosatron in the line mixed the identical stock solution A. In the following, the line was separated and the parallel switched Dosatrons were adding solution B either with a high/or low N concentration. The two lines were supplying two different pipes, each containing 100 micro drips (Netafilm CNL-Tropfer Junior, GVZ-Rossat, Otelfingen, Schweiz); dripping into one of the two nutrient solution reservoirs with a flow of 25 ml min−1 (Fig. 1a, b). Irrigation was controlled by an irrigation computer (Gardena, C1060 Profi, Ulm, Germany). The irrigation cycle was 4 h with duration of one minute. The amounts to 25 mL per micro drip summing up to a daily rate of 300 mL per slide.

The RADIX platform was placed in a cultivation room (720 × 330 × 250 cm). Environmental settings were a day period of 14 h light, a temperature of 28/24 °C (day/night) at seed level, approximately 50 % air humidity and a light intensity of 600 µmol m−2 s−1 photosynthetic active radiation at plant canopy level supplied by GreenPower LED production modules (Koninklijke Philips Electronics N.V., Amsterdam, Netherlands). Sixteen LED modules were arranged in four subunits with four modules each. The subunits were hanging from the ceiling and were equipped with a motor to adjust the distance of LED lamps and plant canopy level. Each subunit was controlled separately. The dimension of the lighting system was 500 × 70 cm and each subunit was 125 × 70 cm.

Imaging

An imaging station (~168 × 164 × 110 cm) was constructed allowing for imaging the root system and the shoot at once (Fig. 1c, Additional file 2). The imaging station was moving on a rail parallel to the slide rack and the rhizoslides were manually slid onto the imaging mount (Fig. 1a, b). The mount was rotatable allowing imaging the root system growing on both sides of the slides with only one camera (Fig. 1c, Additional file 2). The camera was fixed at a distance of 1 m from the center of the acrylic sheet (Figs. 1c, 2a, Additional file 2). For shoot imaging, two webcams (Logitech, Business) were fixed at a distance of 65 cm (increased to 80 cm in the meantime in front of the rhizoslide) and 65 cm (increased to 115 cm) on top of the imaging car (Fig. 1c). The backside of the rhizoslide mount was covered with a blue Komatex plate for an optimal contrast between shoot and background (Fig. 1c). The root zone was hanging in a box constructed with black and white Komatex plates (Röhm AG, Schweiz) at the sides to prevent light scattering in the root zone; white plates in the back of the rhizoslide and black plates in front of it. The white plates in the back area were built as a light box, reflecting the light, whereas the black plates in the front inhibited reflections into the camera focus area (Figs. 1c, 2a). Only the front view was uncovered for root imaging. Black cloth was used to cover the area between camera and box (the cloth is not shown in Figs. 1c, 2a, Additional file 2).

Fig. 2
figure 2

a Imaging car subunit which facilitated analysis of the root system. The root system on a split-root rhizoslide was placed approx. in 1 m distance to the 24 megapixel DSLR-camera. LED bars were placed between camera and rhizoslide on the right and left side, respectively. b Root system developed on one side of the rhizoslide and treated with a split-nitrogen application 10 days after solution change. Roots highlighted in yellow are the crown roots which were used in this example image to investigate crown root growth. Waved brackets indicate the two segments used to investigate the number of lateral roots and the medium (MedLat) and longest (MaxLat) lateral root length in the first (MedLat 1st, MaxLat 1st) and in the second (MedLat 2nd, MaxLat 2nd) segment. The first segment started at the position of the the last developed lateral root (lateralfirst1st) at solution change and the second segment continued after the first segment (laterallast1st). The lateral roots used to investigate the medium and longest lateral root length in the two segments are highlighted in green

Images were taken with a 22.3 mega pixel full-frame digital single-lens reflex camera (EOS 5D Mark III, Canon, Tokyo, Japan) equipped with a 50 mm lens (compact macro 50 mm f/2.5, Canon, Tokyo, Japan). The resolution of the images was around 0.13 mm pixel−1 (for details, see [18]). The flash lights were replaced by six LED bars (TRENOVA -Highpower- LED, PUR-LED GmbH & Co. KG, Selzen, Germany) allowing for a more homogenous illumination (Figs. 1a–c, 2a). Shoots were imaged with the two webcams simultaneously with the root system; one image from the front (maximal extent of the rhizoslide parallel to image plane) one image from the backside (180° turn on the horizontal axis) and images from the top view simultaneously with the side view images (in total four images). As the maximal extension of the leaf area differs between individual plants, imaging from different focal axis allows minimizing the error resulting from shoot orientation. In our setup, leaves were developing in parallel to the rhizoslide axis and therefore, imaging from the front and top rendered sufficient results and no additional imaging from the side was necessary.

Image processing

The two root images taken with either illumination from the right or from the left side were combined to one 24 bit RGB image using Matlab (Version 7.12 The Mathworks, Natick, MA, USA) by keeping only the minimum tonal value present in each image (minimum tonal image) as described by Le Marié et al. [18].

In the current setup, the number of pixels of the front, back and top view of the maize shoot was used as proxy for canopy size. All four shoot images were segmented and used for leaf area calculation. Best results were obtained by calculating an average pixel number over the two sides and top views. Shoot images were segmented from the background using background subtraction and thresholding adaption to the image material written in Matlab R2013a (The Mathworks, Natick, MA, USA). Shoot segmentation was done in a two-step method on a feature map f: f = 2 * G-R-B. Segmenting with a high threshold 55 (in the result pixels of f with values ≥ 55 are set to 1 all others are set to 0) yields a coarse segmentation sc and identifies all plant regions but does not nicely determine their borders. A fine segmentation sf with Otsu’s method [39] resulted in more accurate sizes of the shoot regions but also yields additional image regions which were separated from the shoot. Therefore both segmentations were combined by selecting all regions of sc and enlarging them to all therewith connected regions of sf.

Calculation of the root growth rate

Crown roots of the first whorl were traced using SmartRoot [39]. Crown roots of the first whorl were not always placed on the high and low N side. Accordingly also roots of the second whorl were used. For tracing, one crown root was selected in a high N compartment and one crown root in a low N compartment, ideally with comparable developmental stages at the start of the split root treatment. Furthermore, only crown roots were selected that did not cross the wax barrier or reach the side or bottom of the germination paper at the end of the experiment. This step was necessary to ensure that zero growth in subsequent analyses was not due to the fact that a root reached the borders of the paper and to exclude compensatory growth of lateral roots. All selected roots were healthy and undamaged until the end of the experiment. Next, root elongation of crown roots over time was calculated using R (Foundation for Statistical Computing, Vienna, Austria). Single plants did not establish traceable crown roots on both sides of the split-nutrient germination paper as sometimes crown roots on one side reached the sides or bottom of the rhizoslide before the end of the experiment or no crown roots were formed on one of the split root sides. Though, the root system (seminal and crown roots) developed over the whole germination paper and was exposed to inhomogeneous N distribution. Therefore, it was expected that plant response should be comparable to plants showing traceable crown root growth in both compartments.

Modelling the dynamics of root and shoot elongation

For a better understanding of responses towards a heterogeneous environment, a linear model was fitted to describe changes in root and shoot growth dynamics over time.

A simple linear regression was fit to the length development of each individual crown axile root using a linear regression fit lm() in R [40]:

$$L_{Cr} (t) = a + bt$$
(1)

where L Cr (t) is the crown root length as a function of time t, a the estimated crown axile root length at the start of the split-nitrogen application (y-axis) and b is the slope of the regression line. The shoot growth was calculated simultaneously to the crown root length.

Lateral root length and number

To investigate the effect of genotype and treatment on lateral root growth and number, medium and maximal lateral root length were measured and maximal lateral root number was counted at the last day of the experiment in two branching zones (Fig. 2b). The branching zone, developing after solution change, was defined by the youngest, most distal lateral root at solution change (proximal end of the branching zone) and the youngest, most distal lateral root at harvest (distal end of the branching zone) (Fig. 2b) [41]. The length of the branching zone was measured using the ruler tool in Image J (SmartRoot) and divided into two sections of equal length (first and second branching zone). A medium sized lateral root and the longest lateral root were traced with SmartRoot in both segments and the number of formed lateral roots was counted (Fig. 2b).

Chlorophyll meter (SPAD) measurements

Chlorophyll measurements (SPAD-502, Minolta Corporations, Ramsey, NJ, USA) were done over the experimental period to control for plant health and chlorophyll content. Measurements were done on the last fully developed leaf (usually leaf 4) of every plant measuring at the shoot base, middle and tip of the leaf and averaging these values.

Harvest

After 26 days in total, plants were harvested. The shoot was cut off at the base and leaf area and fresh weight were determined. Leaf area was determined using a leaf-area meter (LI-COR 3100, LI-COR, Lincoln, NE, USA). In order to determine shoot dry weight, shoots were dried at 60 °C for 4 days. Crown and seminal roots were separately harvested depending on the side they were growing (high or low N treatment), bagged and dried at 60 °C for 4 days. Only single lateral roots were growing into the germination paper and could not be harvested. To ensure a correct classification of crown roots, all crown and seminal roots were marked at harvest and an image was taken to compare with the last image taken in the RADIX setup.

CN measurement

After drying, shoots were ground using a swing mill (Retsch MM200, Retsch GmbH, Haan, Germany) with a frequency of 25 s−1 for one minute. N content of the shoots was measured with a CN analyser (Flash EA, 1112 series) [42]. As reference plant material Alfalfa was used.

Verification of the split nutrient system

To verify the separation of the two sides, a white germination paper was soacked in bromocresol green and afterwards watered with a basic (NH4HPO4) and an acidic [Ca(NO3)2] solution. A webcam was installed to image the elution. After 48 h the bromocresol was washed out and the paper was sprayed with it again to visualize the pH gradient (Additional file 3).

Statistics

The rhizoslide experiment was a complete randomized block design with four replications. Each experimental unit consisted of one rhizoslide containing one plant. Linear mixed models were calculated in ASREML-R [43] according to the following full model for root observation in split N application:

$$y_{ijkl} = \mu + G_{i} + N_{j} + GN_{ij} + R_{k} + S_{l} + GR_{ik} + \varepsilon_{ijkl}$$
(2)

where Y ijk is the trait measurement (slope or intercept of crown root length, length lateral root, number of children), of the root of the ith genotype (G = 1, 2, 3 … 24) within the jth N treatment (N = HN, LN) in the rth replication (k = 1, 2, 3, 4) on the lth side of the slide (l = front or back), µ is the general mean, Gi is the effect of genotype, Nj is the effect of the N treatment, GNij is the genotype-treatment interaction, Sl is the effect of the side of the slide the observed root grew on (front or back), R k is the effect of the kth replicate, GRik is the genotype-replication interaction identifying individual slides and ε ijkl is the residual error. The genotype-replication interaction was considered as random whereas all other effects were considered as fixed factors. A different variance was adjusted for each treatment level. For traits representing only one measurement per slide (SPAD, dry weight, leaf area, pixel number, slope or intercept of shoot growth), the model was reduced by only keeping the genotype and replicate as fixed factors. For root traits (Elongation rate crown roots (ERCr), length of crown roots at solution change (intercept; ICCr), length of representative lateral root (MedLat), maximal lateral root length (MaxLat), number of lateral roots (NoLat) and root dry weight (DWR) a square root transformation was performed to achieve normal distribution and no outlier were removed. Root traits that were not normally distributed (length branching zone, branching density in zone 1 or 2, start branching zone) were excluded from further analysis. For shoot traits neither a transformation nor removing of outlier was necessary to achieve normal distribution.

Best linear unbiased estimators (BLUPs) were used to investigate the phenotypic diversity in the set [44]. Tukey’s honestly significant differences (Tukey HSD) were used as post hoc test and calculated as:

$$TukeyHSD = q*\sqrt {\frac{MSE}{n}}$$
(3)

where q is the critical value according to the chosen significance level and degrees of freedom, MSE is the mean square error calculated from the average standard error of the difference (avsed) supplied by the predict function of ASReml-R and n is the number of treatment levels.

Heritability

For heritability estimations across genotypes the model described by Eq. (2) was used, but the genotype and the genotype-treatment interaction were set as random factors. To estimate the heritability depending on the nitrogen treatment applied to the root system, a heterogeneous variance model with regard to the nitrogen-effect and its interactions was fit in ASReml-R following the recommendations of Butler et al. [43].

Mean-based heritability was calculated for each root trait within each nitrogen level and for each shoot trait as follows:

$$H^{2} = \frac{{\sigma_{g}^{2} }}{{\sigma_{g}^{2} + \frac{{\sigma_{e}^{2} }}{r}}}$$
(4)

where \(\sigma_{g}^{2}\) and \(\sigma_{e}^{2}\) are the genotypic and residual error variance respectively obtained from the linear mixed model fitting (Eq. 2) and r is the number of replications.

Results and discussion

RADIX: high throughput analysis of root elongation and dynamics of root topology

The established RADIX platform, which is based on rhizoslides as individual monitoring units, enables a high-throughput screening of RSA and root topology dynamics. One major advantage compared to existing phenotyping platforms such as GROWSCREEN-Rhizo [8] is the visibility of all targeted roots and a sufficient contrast for image processing. No targeted roots are hidden in deeper soil layers like it is the case in rhizotrons and limitations due to contrast between soil and root are usually a minor problem. The advantage of soil-based systems like rhizotrons and CT is that they are closer to natural conditions [8, 10]. Another advantage of rhizoslides is the high achievable throughput. The definition of high throughput usually depends on the analysed traits. Root accessibility and the measurement of growth and topology are certainly the most rate-limiting aspects studying root elongation and dynamics. Here, we define high throughput as the ability to evaluate QTL or association mapping populations of at least 200 genotypes within a minimum of three independent replications within half a year. The imaging process in rhizoslides takes approximately 1 min per slide (front and back side). Image processing in SmartRoot is still the time-limiting step (approximately 80 h to mark axile roots and the segments of lateral roots for a time series of 1400 images) but delivers detailed information of the development of root topology. Adu et al. [13] developed a screening method with medium throughput based on germination paper (n = 24). A scanner unit is fixed in front of each paper with roots growing on it. With this setup, higher temporal resolutions can be realized. However, a mobile imaging unit as realized in the RADIX balances between temporal resolution and throughput. As the response of the root system to stimuli, such as N, happens within days rather than minutes [29], daily or half daily scans are considered sufficient. With the current setup of the RADIX, the number of genotypes that can be evaluated remains the bottleneck. The highest throughput realized to date to study root dynamics in existing platforms is achieved in hydro- or aeroponic setups, yet without the possibility of split-nutrient application. It remains to be elucidated, under which conditions and for which root traits one or the other phenotyping systems is most suitable.

To our knowledge, RADIX is the only system currently existing which allows studying root growth responses under heterogeneous nutrient availabilities with a throughput high enough for QTL and association studies. The importance of root system architecture for N-uptake is supported by several positive QTL collocations between the two traits [45, 46]. However, to date it is still unclear, if increased uptake efficiency, increased N utilization efficiency [47, 48] or an optimal housekeeping with internal N resources [49], is most promising for optimal plant performance [46]. For example, under high N, the genetic variation in nutrient uptake seems to be related to uptake efficiency, whereas at low N, it seems to be related to utilization efficiency [48]. The RADIX setup could shed light onto this topic by studying individual plants under control as well as under stress conditions. In doing so, changes in root and shoot development and the connection to plant performance can be investigated.

Dynamics of shoot development: an imperative for sound root research

The platform was developed to allow assessing both root and shooting growth. In the current setup, the number of pixels of the front and back view of the maize shoot was used as proxy for canopy size. Shoots were successfully segmented from the background using a two-step method of coarse and fine segmentation (Fig. 3c). The number of pixels correlated well with leaf area r2  = 0.91***; Fig. 3a) as well as with shoot biomass (r2 = 0.87**; Fig. 3a). Best correlation could be achieved by calculating the mean value of the two images taken from the sides (rhizoslide parallel to camera focal axis in 0° (front) and 180° (back) and from the top (in total two front views and two top images). In a previous experiment, different imaging methods for the shoot were investigated. Best results were obtained by using the here presented combination of side and top view. This is in contradiction with investigations made in the GROWSCREEN-Rhizo facility where best results were obtained with images from the front side (0°) and a 90° horizontal rotation [8]. An explanation for this phenomenon could be the combination of the alternate leaf orientation in monocot species and a parallel orientation to the rhizoslide influenced by neighbouring plants (distance between two neighbouring plants of 5 cm) [45]. As leaves exceeded the imaging area 10 days after start of the treatment, the shoots dynamics could be assessed only up to this point. For future experiments, this was improved by constructing a larger vertical screen for the imaging station (already displayed in Fig. 1).

Fig. 3
figure 3

a Correlation between the average number of pixel (Pixel) of the top and side images and the leaf area (cm−2) or shoot dry weight (g) measured destructively at the end of experiment. b Correlation of shoot pixel count development (Pixel day−1) expressed as the slope best fitting the linear model and the number of pixels (Pixel), leaf area (cm−2) and shoot dry weight (g). a, b Numbers represent the average of all replicates of a single genotype. For the corresponding genotypes see “Cultivation in the RADIX platform” in the “Methods” section. The standard error of the difference was 3392 pixels for the average number of pixel, 18.01 cm−2 for the leaf area, 0.05 g for the shoot dry weight and 315 pixels for the shoot pixel count development. c Growth series of an exemplary shoot and images of the segmented shoot

Evaluation of shoot traits

Genotypes differed with respect to leaf area expansion (measured as pixel numbers), leaf area at harvest and shoot dry weight at harvest (Table 1a, b). For example, genotype A347 grew almost two times faster than genotype Oh33 (Additional file 4, Additional file 5B). As the traits, measured at final harvest, were strongly correlated with leaf area expansion, it can be concluded that final leaf area at this stage was driven by differences in vigour rather than differences in germination (Fig. 3b). Furthermore, the shoot pixel count at solution change was not correlated with the shoot pixel count development (Fig. 4, Additional file 6). No genotype-specific differences and no heritability were observed for chlorophyll content investigated by means of SPAD measurements (Table 1a, b). However, the genotypes differed in their N content (Table 1a). CN measurements revealed significant differences in N content between genotypes with a maximal N proportion of 5.07 % and a minimal proportion of 3.61 % (Table 1a, b). A high formation of biomass resulted in a reduction of leaf N content (Fig. 4) although this correlation was not significant (Additional file 6).

Table 1 (a) Results of the analysis of variance, (b) maximal, minimal and median values and heritability of traits
Fig. 4
figure 4

Principal component analysis (PCA) of traits based on best linear unbiased estimates (BLUPS). Abbreviations: Number of lateral roots in the first segment (NoLat 1st) or second segment (NoLat 2nd), length of representative lateral root in the first segment (MedLat 1st) or second segment (MedLat 2nd), maximal lateral root length in the first segment (MaxLat 1st) or in the second segment (MaxLat 2nd), total number of lateral roots (NoLatTot), elongation rate crown roots (ERCr), length of crown roots at solution change (intercept; ICCr), embryonic root dry weight (DWER), crown root dry weight (DWCR), leaf greenness (SPAD), leaf area measured (LAm), dry weight shoot (DWS), leaf area pixel based (LAPix), shoot pixel count at solution change intercept (ICS), shoot pixel count development (ERS), N content in the leaf in % of total dry weight (N). The numbers in the PCA plot correspond to the EU_IDs. For the corresponding genotypes see “Cultivation in the RADIX platform” in the “Methods” section

Crown root development

The elongation of crown axile roots was described best by linear models based on the observation of original data and residual plots (Additional files 7, 8, 9). In a previous study, best results were obtained by fitting a logistic model for the elongation rates of the crown roots on the high-N side and a regression model with segmented relationships for roots grown under no N [29]. The need for two different modelling approaches in the previous study did not enable the examination of interaction effects in a common linear mixed model. As already observed by in ’t Zandt et al. [29], crown root elongation tended to decrease in the low-N compartment (Additional file 8) while the elongation rate tended to increase in the high-N compartment (Additional file 7). However, the trend was far less pronounced in the present experiment enabling the fit of a common linear model (Additional file 9). The lower effects were possibly due to the replacement of the no-N treatment [29] by a low-N treatment. Unfortunately, the standard genotype UH007 × B73, used in the study of in ’t Zandt et al. [29], could not be investigated here, as it did not germinate.

First crown roots established on the germination paper after 2 weeks and their response to split-nutrient application could be measured for another 2 weeks. By then, they had reached the bottom or sides of the rhizoslides and the plants had reached four fully developed leaves. The realization of longer cultivation phases allowing the study of crown root dynamics is a major advantage compared to the growth pouch system used in past studies [17]. Genotypes differed in their developmental speed. While some genotypes started to form crown roots already seven days after start of the experiment, no crown roots were visible on the germination paper for others at the time point of solution change. The evaluation of the dynamic development of individual roots is certainly a way to account for such developmental differences.

For the genotypes F1808, EZ37, FC1890 and W64A only one or two replicates could be evaluated under low or high N as only those did develop crown roots until the time point of solution change (Additional files 7, 8). Twelve individual plants did not develop a traceable crown root on the high N side until the time point of solution change and seventeen on the low N side (Additional files 7, 8). We used a linear mixed model approach which handles missing values better than classical linear models [50]. One opportunity to increase the number of traceable crown roots would be to postpone the treatment to later developmental stages and another option would be to increase the number of replications.

Differences in the amount of roots growing at the low or high N side might influence the morphological responses of the roots to local high nitrate. However, in ‘t Zandt et al. [29] found no direct dependence of the dynamic of roots exposed to differential N supply in the rhizoslides concluding that the responses on either side occurred largely independent from each other.

Genotype-by-nitrogen interaction for crown root development

The crown root elongation under low N was on average 26 % reduced compared to the high N side (Table 1b). The minimum observed elongation in the high N compartment was 1.44 cm d−1 and the maximal elongation was 1.77 cm d−1 (Table 1b). There were differences for the intercept, i.e. the crown root length at the time of solution change.

There was a significant genotype-by-treatment interaction for the response of crown root growth to split-nutrient application. More specifically, the crown axile roots at different N-concentrations showed genotype-dependent differential response. For example, crown roots of the genotype PH207 grew 39 % slower under low N conditions than under high N conditions whereas the genotype LAN 496 showed a much lower reduction in growth (15 %) under low N conditions (Table 1b, Additional file 7, 8). Interestingly, the heritability was high (h2 = 0.41) for roots grown under high N but zero for roots grown under low N conditions (Table 1b). This means that it is not possible to select genotypes based on their rooting behaviour on the low-N side of a split-nutrient setup. By contrast, the fast exploration of N rich patches and selective root placement is highly heritable and could be an advantageous trait with respect to inhomogeneity in the field or patch wise application of fertilizers e.g. by row fertilization. According to a model by Dunbabin et al. [51], selective root placement could lead to a more than two fold higher N uptake efficiency throughout the whole growing season, compared to a root system that is only poorly capable to respond to N. Already the study of in ’t Zandt et al. [29] demonstrated a oppositional response in the low N and high N compartment. However the study was limited to one genotype. To our knowledge, this is the first time that such a selective response could be observed for a larger set of maize genotypes.

Length and density of lateral roots

Investigations of lateral root number and length were only done at the end of the experiment as in the study of in ’t Zandt [29] the maximal number of lateral roots in the segment was already reached two days after solution change in the high N compartment and after four days in the low N compartment. Thus, it was assumed that the final number of lateral roots is representative for the branching intensity in these segments. Furthermore, the number of laterals in the zone present before solution change did not significantly increase during the two weeks of treatment, hypothesizing that only newly developed root tissues are able to respond to N availability (Additional file 10).

In general, lateral roots became longer under high N compared to low N conditions (Table 1b, Additional file 5A-C). This finding is in close agreement with previous studies reporting a selective root placement within the high N compartment by stronger lateral root formation and growth [52, 53]. However, some genotypes were characterized by an opposed trend highlighting the opportunity to identify contrasting genotypes. Past studies already observed different responses for different species and even different genotypes within a species to the same environmental stimulus [54], but to date no study tried to use contrasting genotypic responses for mapping of quantitative traits.

The length of the branching zone (Tab 1. LBrZ) did not differ among genotypes but the number of lateral roots (Table 1, NoLat1st and NoLat2nd) in this zone was heritable and correlated with the axile root growth [r = 0.47* (NoLat1st) and r = 0.49* (NoLat2nd)]. Interestingly, the genotypes did not differ for the linear density of lateral roots per unit axile root length. Hence, the formation of lateral roots seemed to be mainly driven by a stronger axile root elongation rather than a higher branching density as no significant treatment effects could be observed for the branching density or the length of the branching zone. Although the branching intensity was the only trait that was correlated with leaf N content (Additional file 6), the missing heritability exclude the integration of this trait into breeding schemes.

Genotypic differences were observed for all root traits except the maximal lateral root length in the first segment and the dry weight of embryonic roots. Accordingly, the heritability of the number and length of lateral roots was consistently moderate to high in both segments except for the maximal length in the first segment [h2 = 0 (low N); h2 = 0.41 (high N)] whereas the heritability was higher in the high N (0.41 ≤ h2 ≥ 0.72) than in the low N compartment (0.03 ≤ h2 ≥ 0.59) in general (Table 1b).

Genotype-by-nitrogen interaction for lateral root characteristics

A significant genotype-by-N placement interaction was observed for the medium and maximum length of lateral roots. The ability of genotypes to form long first-order lateral roots and its genetic control is a very interesting research topic. Eventually such roots are able to replace some of the functions of their parental axile root, thus allowing a more flexible response of a root system.

Relationship between root and shoot development

We found a strong connection between lateral root formation and shoot performance. We used correlations based on genotypic mean values to evaluate the dependencies between roots formed on the high and low N side and between roots and shoots. Root traits were strongly correlated among themselves (Fig. 4) and two major clusters of root traits could be identified by principal component analysis: the lateral root length and the number of lateral roots (Fig. 4). Interestingly, the number of laterals was closely correlated with the elongation of crown roots (Fig. 4) supporting the hypothesis that the number of laterals was mainly driven by axial root elongation rather than the branching density.

The leaf area was positively correlated to the medium length and number of lateral roots (Figs. 4, 5, Additional file 6) whereas the medium length was only under high N conditions correlated to shoot growth (Fig. 5). Lengths and densities of lateral roots measured in the second segment were stronger positively correlated with shoot traits than the first segment (Additional file 6). This inconsistency might be related to a feedback loop. N starvation during the establishment phase may have led to a limited link between shoot and root growth whereas in the second segment, the roots already profit from a higher N uptake resulting in a positive feedback from the shoot. This observation supports the hypothesis of de Kroon et al. [55] suggesting a coupling of nutrient sensing and coordinated growth of the root components and is in line with observations of in ’t Zandt et al. [29] demonstrating a dynamic differential response of crown and lateral root growth under either zero N or high N conditions. Furthermore, it indicates a signalling between roots and shoot transmitting the information about N availability that itself is transformed into differential growth e.g. by the provision of carbohydrates. Indeed, Tabata et al. [56] observed a signalling from the root to the shoot triggering the activity of nitrate transporters within root regions with high nitrate availability in Arabidopsis thaliana. An alternative explanation of this phenomenon could be a better N supply by the roots after the formation of laterals in the first segment resulting in an increasing root surface area and a simultaneously increased uptake of N after the change of the nutrient solution. Furthermore, a very high N stress level can lead to an inhibition of lateral root formation whereas a moderate stress can induce it [57].

Fig. 5
figure 5

Correlation of shoot pixel count development (Pixel day−1) expressed as the slope best fitting the linear model and medium lateral root length (cm) or the number of lateral roots (counts) either under high N (a) or low N conditions (b). Significance levels: ≤0.001***, ≤0.01**, ≤0.05*. Correlation was done based on best linear unbiased predictors (BLUPS). The standard error of the difference was 323 pixels for the shoot and 0.1 cm for the lateral root length under high as well as under low N conditions

Considerations concerning a maize ideotype for improved foraging behaviour

There are many open questions with respect to the optimal root ideotype of maize. Foraging of N in nutrient-rich patches is mainly achieved by changing root system architecture. While, in general, the root length within the nutrient-rich patch increases, root growth outside the patch decreases [52, 53, 58, 59]. Thus, there is a trade-off between intensive foraging in patches and the overall direction of the root system to an available resource which might become important at critical stages of development, e.g. water and N at depth during grain filling. Accordingly, the proposed steep, cheap and deep ideotype to optimize water and N acquisition at depth to capture leached nitrate [7] includes “unresponsiveness of lateral branching to localized resource availability”. Paper-based systems offer a huge potential to establish rapid screens for N responsiveness.

This screening already provides an insight into the diversity of responses of genotypes towards inhomogeneous nutrient distributions and gives a strong indication that the responsiveness is closely linked to shoot development. However, further investigations with contrasting genotypes under field conditions are necessary, before drawing a conclusion on the utility of the presented screening method for selection purposes.

Conclusion

The RADIX platform allowed studying dynamic changes in root-system architecture to split-root application of nitrogen. This is the first study evaluating such a differential response using a larger set of maize genotypes. A stronger selective root placement in the high N-compartment was related to an increased shoot development. This indicates that high early vigour might be related to a more intense foraging behaviour. In ongoing experiments, we aim to (1) verify these results under field conditions and (2) map the genomic regions controlling these responses using an association panel. Apart from studying nutrient-use efficiency, the system may be also used to evaluate responses to stresses like extreme pH, agrochemicals or microbes.

Abbreviations

RSA:

root system architecture

N:

nitrogen

QTL:

quantitative trait loci

References

  1. Hund A, Ruta N, Liedgens M. Rooting depth and water use efficiency of tropical maize inbred lines, differing in drought tolerance. Plant Soil. 2009;318:311–25.

    Article  CAS  Google Scholar 

  2. Watt M, Magee LJ, McCully ME. Types, structure and potential for axial water flow in the deepest roots of field-grown cereals. New Phytol. 2008;178:135–46.

    Article  PubMed  Google Scholar 

  3. Zhu J, Kaeppler SM, Lynch JP. Topsoil foraging and phosphorus acquisition efficiency in maize (Zea mays). Funct Plant Biol. 2005;32:749–62.

    Article  CAS  Google Scholar 

  4. Smith S, De Smet I. Root system architecture: insights from Arabidopsis and cereal crops. Philos Trans R Soc Lond B Biol Sci. 2012;367:1441–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Burton AL, Brown KM, Lynch JP. Phenotypic diversity of root anatomical and architectural traits in species. Crop Sci. 2013;53:1042–55.

    Article  Google Scholar 

  6. Worku M, Bänziger M, Schulte auf’m Erley G, Friesen D, Diallo AO, Horst WJ. Nitrogen efficiency as related to dry matter partitioning and root system size in tropical mid-altitude maize hybrids under different levels of nitrogen stress. Field Crop Res. 2012;130:57–67.

    Article  Google Scholar 

  7. Lynch JP. Steep, cheap and deep: an ideotype to optimize water and N acquisition by maize root systems. Ann Bot Lond. 2013;112:347–57.

    Article  CAS  Google Scholar 

  8. Nagel KA, Putz A, Gilmer F, Heinz K, Fischbach A, Pfeifer J, Faget M, Blossfeld S, Ernst M, Dimaki C, et al. GROWSCREEN-Rhizo is a novel phenotyping robot enabling simultaneous measurements of root and shoot growth for plants grown in soil-filled rhizotrons. Funct Plant Biol. 2012;39:891–904.

    Article  Google Scholar 

  9. Heeraman D, Hopmans J, Clausnitzer V. Three dimensional imaging of plant roots in situ with X-ray computed tomography. Plant Soil. 1997;189:167–79.

    CAS  Google Scholar 

  10. Gregory PJ, Hutchison D, Read D, Jenneson P, Gilboy W, Morton E. Non-invasive imaging of roots with high resolution X-ray micro-tomography. Plant Soil. 2003;255:351–9.

    Article  CAS  Google Scholar 

  11. Blossfeld S, Le Marié C, Van Dusschoten D, Suessmilch S, Kuhn A. Non-invasive investigation of root growth via NMR imaging. Commun Agric Appl Biol Sci. 2011;76:11.

    PubMed  Google Scholar 

  12. Rascher U, Blossfeld S, Fiorani F, Jahnke S, Jansen M, Kuhn AJ, Matsubara S, Märtin LLA, Merchant A, Metzner R. Non-invasive approaches for phenotyping of enhanced performance traits in bean. Funct Plant Biol. 2011;38:968–83.

    Article  CAS  Google Scholar 

  13. Adu MO, Chatot A, Wiesel L, Bennett MJ, Broadley MR, White PJ, Dupuy LX. A scanner system for high-resolution quantification of variation in root growth dynamics of Brassica rapa genotypes. J Exp Bot. 2014;65:2039–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Planchamp C, Balmer D, Hund A, Mauch-Mani B. A soil-free root observation system for the study of root-microorganism interactions in maize. Plant Soil. 2013;367:605–14.

    Article  CAS  Google Scholar 

  15. Zhu J, Kaeppler SM, Lynch JP. Mapping of QTLs for lateral root branching and length in maize (Zea mays L.) under differential phosphorus supply. Theor Appl Genet. 2005;111:688–95.

    Article  CAS  PubMed  Google Scholar 

  16. Liao H, Rubio G, Yan X, Cao A, Brown KM, Lynch JP. Effect of phosphorus availability on basal root shallowness in common bean. Plant Soil. 2001;232:69–79.

    Article  CAS  PubMed  Google Scholar 

  17. Hund A, Trachsel S, Stamp P. Growth of axile and lateral roots of maize: I development of a phenotying platform. Plant Soil. 2009;325:335–49.

    Article  CAS  Google Scholar 

  18. Le Marié C, Kirchgessner N, Marschall D, Walter A, Hund A. Rhizoslides: paper-based growth system for non-destructive, high throughput phenotyping of root development by means of image analysis. Plant Methods. 2014;10:13.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Ruta N, Liedgens M, Fracheboud Y, Stamp P, Hund A. QTLs for the elongation of axile and lateral roots of maize in response to low water potential. Theor Appl Genet. 2009;120:621–31.

    Article  PubMed  Google Scholar 

  20. Ruta N, Stamp P, Liedgens M, Fracheboud Y, Hund A. Collocation of QTLs for seedling traits and yield components of tropical maize under limited water availability. Crop Sci. 2010;50:1385–92.

    Article  CAS  Google Scholar 

  21. Trachsel S, Messmer R, Stamp P, Hund A. Mapping of QTLs for lateral and axile root growth of tropical maize. Theor Appl Genet. 2009;119:1413–24.

    Article  CAS  PubMed  Google Scholar 

  22. Reimer R. Fesponses of maize (Zea mays L.) seedlings to low and high temperature: association mapping of root growth and photosynthesis-related traits. PhD Thesis. Diss. ETH No. 18807, AGRL; 2010.

  23. Hochholdinger F, Park WJ, Sauer M, Woll K. From weeds to crops: genetic analysis of root development in cereals. Trends Plant Sci. 2004;9:42–8.

    Article  CAS  PubMed  Google Scholar 

  24. Hund A, Reimer R, Messmer R. A consensus map of QTLs controlling the root length of maize. Plant Soil. 2011;344:143–58.

    Article  CAS  Google Scholar 

  25. Yu P, White PJ, Hochholdinger F, Li C. Phenotypic plasticity of the maize root system in response to heterogeneous nitrogen availability. Planta. 2014;240:667–78.

    Article  CAS  PubMed  Google Scholar 

  26. Yu P, Li X, Yuan L, Li C. A novel morphological response of maize (Zea mays) adult roots to heterogeneous nitrate supply revealed by a split-root experiment. Physiol Plant. 2014;150:133–44.

    Article  CAS  PubMed  Google Scholar 

  27. Peng Y, Li X, Li C. Temporal and spatial profiling of root growth revealed novel response of maize roots under various nitrogen supplies in the field. PLoS ONE. 2012;7:e37726.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Jackson R, Manwaring J, Caldwell M. Rapid physiological adjustment of roots to localized soil enrichment. Nature. 1990;344:58–60.

    Article  CAS  PubMed  Google Scholar 

  29. in ’t Zandt D, Le Marié C, Kirchgessner N, Visser E, Hund A. High-resolution quantification of root dynamics in split-nutrient rhizoslides reveals rapid and strong proliferation of maize roots in response to local high nitrogen. J Exp Bot. 2015;66:5507–17.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Vadez V. Root hydraulics: the forgotten side of roots in drought adaptation. Field Crop Res. 2014;165:15–24.

    Article  Google Scholar 

  31. Wissuwa M, Mazzola M, Picard C. Novel approaches in plant breeding for rhizosphere-related traits. Plant Soil. 2009;321:409–30.

    Article  CAS  Google Scholar 

  32. Grieder C, Trachsel S, Hund A. Early vertical distribution of roots and its association with drought tolerance in tropical maize. Plant Soil. 2014;377:295–308.

    Article  CAS  Google Scholar 

  33. Hartmann A, Czauderna T, Hoffmann R, Stein N, Schreiber F. HTPheno: an image analysis pipeline for high-throughput plant phenotyping. BMC Bioinform. 2011;12:148.

    Article  Google Scholar 

  34. Hairmansis A, Berger B, Tester M, Roy S. Image-based phenotyping for non-destructive screening of different salinity tolerance traits in rice. Rice. 2014;7:16.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Reymond M, Muller B, Leonardi A, Charcosset A, Tardieu F. Combining quantitative trait loci analysis and an ecophysiological model to analyze the genetic variability of the responses of maize leaf growth to temperature and water deficit. Plant Physiol. 2003;131:664–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Freund H. Sterilisieren der fertigen Ampullen. In: Freund H, editor. Die Ampullenfabrikation. Berlin: Springer; 1916. p 61–72.

  37. Bohn M, Novais J, Fonseca R, Tuberosa R, Grift T. Genetic evaluation of root complexity in maize. Acta Agron Hung. 2006;54:291–303.

    Article  CAS  Google Scholar 

  38. Shaff JE, Schultz BA, Craft EJ, Clark RT, Kochian LV. GEOCHEM-EZ: a chemical speciation program with greater power and flexibility. Plant Soil. 2010;330:207–14.

    Article  CAS  Google Scholar 

  39. Otsu N. A threshold selection method from gray-level histograms. IEEE Trans Syst Man Cybern. 1979;9:62–6.

    Article  Google Scholar 

  40. R Core Team: R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria 2015. http://www.R-project.org/.

  41. NPTEL. http://nptel.ac.in/courses/102103016/module1/lec1/3.html. Accessed 17 June 2016.

  42. Abbeddou S, Diekmann J, Rischkowsky B, Kreuzer M, Oberson A. Unconventional feeds for small ruminants in dry areas have a minor effect on manure nitrogen flow in the soil–plant system. Nutr Cycl Agroecosyst. 2013;95:87–101.

    Article  Google Scholar 

  43. Butler D: asreml: asreml() fits the linear mixed model. R package version 3.0. www.vsni.co.uk. 2009.

  44. Piepho H-P, Möhring J. Computing heritability and selection response from unbalanced plant breeding trials. Genetics. 2007;177:1881–8.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Coque M, Martin A, Veyrieras J, Hirel B, Gallais A. Genetic variation for N-remobilization and postsilking N-uptake in a set of maize recombinant inbred lines. 3. QTL detection and coincidences. Theor Appl Genet. 2008;117:729–47.

    Article  CAS  PubMed  Google Scholar 

  46. Garnett T, Conn V, Kaiser BN. Root based approaches to improving nitrogen use efficiency in plants. Plant Cell Environ. 2009;32:1272–83.

    Article  CAS  PubMed  Google Scholar 

  47. Clark RB. Physiology of cereals for mineral nutrient uptake, use, and efficiency. In: Baligar VC, Duncan RR, editors. Crops as enhancers of nutrient use. San Diego: Academic Press, Inc; 1990. p. 131–209.

  48. Moll R, Kamprath E, Jackson W. Analysis and interpretation of factors which contribute to efficiency of nitrogen utilization. Agron J. 1982;74:562–4.

    Article  Google Scholar 

  49. Lynch J. Root architecture and plant productivity. Plant Physiol. 1995;109:7.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Piepho H, Möhring J, Melchinger A, Büchse A. BLUP for phenotypic selection in plant breeding and variety testing. Euphytica. 2008;161:209–28.

    Article  Google Scholar 

  51. Dunbabin V, Rengel Z, Diggle A. Simulating form and function of root systems: efficiency of nitrate uptake is dependent on root system architecture and the spatial and temporal variability of nitrate supply. Funct Ecol. 2004;18:204–11.

    Article  Google Scholar 

  52. Drew M. Comparison of the effects of a localised supply of phosphate, nitrate, ammonium and potassium on the growth of the seminal root system, and the shoot, in barley. New Phytol. 1975;75:479–90.

    Article  CAS  Google Scholar 

  53. Drew MC, Saker LR, Ashley TW. Nutrient supply and the growth of the seminal root system in barley: I. The effect of nitrate concentration on the growth of axes and laterals. J Exp Bot. 1973;24:1189–202.

    Article  CAS  Google Scholar 

  54. Forde B, Lorenzo H. The nutritional control of root development. Plant Soil. 2001;232:51–68.

    Article  CAS  Google Scholar 

  55. De Kroon H, Visser EJ, Huber H, Mommer L, Hutchings MJ. A modular concept of plant foraging behaviour: the interplay between local responses and systemic control. Plant, Cell Environ. 2009;32:704–12.

    Article  Google Scholar 

  56. Tabata R, Sumida K, Yoshii T, Ohyama K, Shinohara H, Matsubayashi Y. Perception of root-derived peptides by shoot LRR-RKs mediates systemic N-demand signaling. Science. 2014;346:343–6.

    Article  CAS  PubMed  Google Scholar 

  57. Forde BG. Nitrogen signalling pathways shaping root system architecture: an update. Curr Opin Plant Biol. 2014;21:30–6.

    Article  CAS  PubMed  Google Scholar 

  58. Granato T, Raper C. Proliferation of maize (Zea mays L.) roots in response to localized supply of nitrate. J Exp Bot. 1989;40:263–75.

    Article  CAS  PubMed  Google Scholar 

  59. Hackett C. A method of applying nutrients locally to roots under controlled conditions, and some morphological effects of locally applied nitrate on the branching of wheat roots. Aust J Biol Sci. 1972;25:1169–80.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

CLM and AH, PF JP, NK and AW designed the RADIX system, PF constructed the system, NK designed the illumination unit and programmed the imaging software. CLM and JP performed the experiment. CLM and NK did the image analysis. AH and CLM did the statistical analysis and drafted the manuscript supported by AW. All authors read and approved the final manuscript.

Acknowledgements

The authors thank Brigitta Herzog, Hansueli Zellweger and Frank Liebisch from ETH Zürich for their support handling the experiment. Thanks to Delley Seeds and Plants Ltd. for hybrid production. We kindly thank the donors of the genetic material: Department of Agroenvironmental Science and Technologies (DiSTA), University of Bologna, Italy (RootABA lines); Estación Experimental de Aula Dei (CSIC), Spain (EZ47, EZ11A, EZ37); Centro Investigaciones Agrarias de Mabegondo (CIAM), Spain (EC169); University of Hohenheim, Versuchsstation für Pflanzenzüchtung, Germany (UH007, UH250); and INRA CNRS UPS AgroParisTech, France (supply of the remaining, public lines).

Competing interests

The authors declare that they have no competing interests.

Funding

This research received funding from the European Community Seventh Framework Programme FP7-KBBE-2011-5 under Grant Agreement No. 289300.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Hund.

Additional files

Additional file 1. Constructional drawing of the RADIX platform (units:mm).

Additional file 2. Constructional drawing of the imaging station (units:mm).

13007_2016_140_MOESM3_ESM.avi

Additional file 3. The video shows the eluation of the germination paper soacked in bromocresol green and afterwards watered with a basic (NH4HPO4) and an acidic (Ca(NO3)2) solution. A webcam was installed to image the elution. After 48 hours the bromocresol was washed out and the paper was sprayed with it again to visualize the pH gradient.

13007_2016_140_MOESM4_ESM.pdf

Additional file 4. Shoot growth of 24 different genotypes quantified by the increase in number of pixels after the start of the split root treatment.

13007_2016_140_MOESM5_ESM.pdf

Additional file 5. A) Best linear unbiased prediction (BLUPS) of parameters with significant treatment effect, but no genotype effect or genotype: treatment interaction. B) Prediction of mean values of parameters with significant genotype effect, but no genotype: treatment interaction. C) Prediction of mean values of parameters with significant genotype: treatment interaction. Abbreviations: Number of lateral roots in the first segment (NoLat 1st) or second segment (NoLat 2nd), length of representative lateral root in the first segment (MedLat 1st) or second segment (MedLat 2nd), maximal lateral root length in the first segment (MaxLat 1st) or in the second segment (MaxLat 2nd), elongation rate crown roots (ERCr), length of crown roots at solution change (intercept; ICCr), embryonic root dry weight (DWER), crown root dry weight (DWCR), chlorophyll measurements (SPAD, measured leaf area (LAm), shoot dry weight (DWS), number of pixels specifying the leaf area (LAPix), shoot pixel count at solution change (intercept; ICS), shoot pixel count development (ERS), N content in the leaf in % of total dry weight (N in %).

13007_2016_140_MOESM6_ESM.pdf

Additional file 6. Correlation coefficients (Pearson’s r) of traits. Correlation was done based on best linear unbiased estimates (BLUPS). Abbreviations: Number of lateral roots in the first segment (NoLat 1st) or second segment (NoLat 2nd), length of representative lateral root in the first segment (MedLat 1st) or second segment (MedLat 2nd), maximal lateral root length in the first segment (MaxLat 1st) or in the second segment (MaxLat 2nd), branching density in the first segment (BrLat1st) or in the second segment (BrLat2nd), branching density across both segments (BrLatTot), length of the branching zone (LBrZone), total number of lateral roots (NoLatTot), elongation rate crown roots (ERCr), length of crown roots at solution change (intercept; ICCr), embryonic root dry weight (DWER), crown root dry weight (DWCR), leaf greenness (SPAD), leaf area measured (LAm), dry weight shoot (DWS), leaf area pixel based (LAPix), shoot pixel count at solution change intercept (ICS), shoot pixel count development (ERS), N content in the leaf in % of total dry weight (N). Significance level: ≤ 0.001 ***, ≤ 0.01**, ≤ 0.05*.

13007_2016_140_MOESM7_ESM.pdf

Additional file 7. Increase in crown root length of 24 genotypes under high nitrogen after the start of the split root treatment.

13007_2016_140_MOESM8_ESM.pdf

Additional file 8. Increase in crown root length of 24 genotypes under low nitrogen after the start of the split root treatment.

13007_2016_140_MOESM9_ESM.pdf

Additional file 9. Residual vs. fitted crown root length of multiple linear models to determine intercept and slope of crown root development after solution change.

Additional file 10. Increase of lateral roots after solution change above the last formed lateral.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Le Marié, C., Kirchgessner, N., Flütsch, P. et al. RADIX: rhizoslide platform allowing high throughput digital image analysis of root system expansion. Plant Methods 12, 40 (2016). https://doi.org/10.1186/s13007-016-0140-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13007-016-0140-8

Keywords