RGB och spektrala Root Imaging för växt fenotypning och fysiologisk forskning: experiment och Imaging protokoll
Summary
Ett experimentellt protokoll presenteras för bedömning av marken odlas växter rotsystem med RGB och Hyperspektrala imaging. Kombination av RGB-bild, tid serien med chemometric information från Hyperspektrala File optimerar insikter i växten rot dynamics.
Abstract
Bättre förståelse av växten rot dynamics är nödvändigt att förbättra resurseffektiviteten användning av jordbrukssystem och öka motståndet av gröda sorter mot miljöpåfrestningar. Ett experimentellt protokoll presenteras för RGB och Hyperspektrala avbildning av rotsystem. Metoden används rhizoboxes där växter växer i naturlig jord över en längre tid att observera fullt utvecklade rotsystem. Experimentella inställningarna exemplifieras för att bedöma rhizobox växter under vattenstress och studerar roll för rötter. En RGB-imaging installation beskrivs för billig och snabb kvantifiering av roten utveckling över tid. Hyperspektrala imaging förbättrar rot segmentering från jord bakgrunden jämfört med RGB-färg baserat tröskelvärde. Hyperspektrala imaging viss styrka är förvärvet av chemometric information om rot-smutsa systemet för funktionell förståelse. Detta visas med hög upplösning vatten innehåll kartläggning. Spectral imaging är emellertid mer komplex bild förvärvandet, bearbetning och analys jämfört med metoden RGB. En kombination av båda metoderna kan optimera en omfattande bedömning av rotsystemet. Applikationsexempel integrera rot och aboveground egenskaper ges för ramen för fenotypning och växters fysiologiska växtforskning. Ytterligare förbättring av roten imaging kan erhållas genom att optimera RGB bildkvalitet med bättre belysning med olika ljuskällor och i förlängningen av bild analysmetoder att sluta på roten zonegenskaper från spektrala data.
Introduction
Rötterna ger flera viktiga funktioner för växter såsom lagring av assimilerar, ankringen av landlevande växter i jord, och upptag och transport av vatten och näringsämnen1. Ur en evolutionär synvinkel anses bildandet av roten axlar vara en grundläggande förutsättning för beskärningen av mark växter2. Trots denna viktiga roll har historiskt rötter ockuperat endast en marginell position inom biologisk forskning. I nyare tider, men ökar det vetenskapliga intresset växt rotsystem som framgår i figur 1.
Figur 1: Betydelsen av roten studier i botanik.
Antal root relaterade studier som en procentandel av alla publicerade växt studier i SCI tidskrifter under de senaste årtiondena. Sökresultat från Scopus med sökord ”anläggning” och ”planterar och rot”. Klicka här för att se en större version av denna siffra.
Två anledningar kan vara hypotesen att ligger bakom de senaste framstegen inom rot forskning. Första, terrestrial vegetationen utsätts oftare yttre påfrestningar till följd av globala förändringar3. I samband med vegetabiliska jordbruksprodukter produktion uppskattas det att globalt cirka 30% av jordbruksarealen begränsas av vatten och fosfor4,5. Stressreducering av skördarna är en huvudorsak för betydande avkastning luckor som är globalt uppskattas lägre 50 procent av potentiell produktivitet för rainfed jordbruksekosystemen6. Förutom låg resurstillgänglighet, är detta också relaterat till dålig resurseffektivitet användning, dvs otillräcklig kapacitet av en växt att utnyttja tillgängliga resurser7. Detta resulterar i förluster av mobila resurser såsom nitrat som kan inverka negativt på andra ekosystem. Den nuvarande globala kväve effektiviteten beräknas till exempel till 47%8. Bättre resursanvändning effektivitet via metoder för förbättrad hantering och sorter är därför av stor betydelse för både uthållig tillväxt av jordbruks utgångar samt när det gäller miljömässig hållbarhet. I detta sammanhang växten anses rötter vara ett huvudmål för förbättrade grödor och beskärning system9,10.
En andra viktig bakgrund för det senaste intresset för växternas rötter är tekniska framsteg i mätmetoder. Root metoder har länge begränsats av två viktiga utmaningar: för mätning av rötter från växter som växer i jord som de hade att vara isolerad för kvantifiering, mestadels genom att tvätta11, därmed störa det arkitektoniska arrangemanget av roten axlar. In situ – roten observation med utgrävning metoder, därmed bevara den naturliga platsen för rötter i jord, har använts för botaniska Beskrivning12. Fortfarande de är mycket tidskrävande och därmed att kraven inte genomströmning av jämförande strukturella-funktionella rotsystem analys. Å andra sidan var hög genomströmning metoder för roten arkitekturen mätning mestadels gjort på konstgjord medier och för plantor växter13 där en extrapolering av växter naturliga tillväxt miljön är tvivelaktiga14.
Den senaste högkonjunkturen rot forskning är tätt knuten till förskottet i imaging metoder15. Imaging metoder i roten studier kan grovt grupperas i tre typer. Först finns det högupplösta 3D metoder såsom CT och Mr16. Dessa metoder är mest lämpade att studera interaktion processer av växternas rötter med jord, såsom torka inducerad xylem embolism17. Vanligtvis används de på jämförelsevis små prover där de tillåter detaljerade observationer. En jämförelse av CT och MRT för olika stora krukor och fina rot avbildning finns i18. För det andra finns det hög genomströmning bildgivande metoder19,20. Dessa metoder är mestadels baserade på gemensamma 2D RGB avbildning av rötter växa på konstgjorda media (gel, grobarhet papper) där hög kontrast gör jämförelsevis enkel dissektion mellan rötter och bakgrund. De är lämpliga för hög genomströmning jämförelse bland fröplanta rot drag av olika grödor genotyper enligt standardiserade artificiella växande villkor13. Mellan dessa två metoder är metoder som rhizobox: de använder 2D avbildning av rötter som växer i jord under längre tidsperiod och har medelstora genomströmning21,22. En senaste utmaning till (2D) roten imaging är att fånga även indikatorer på roten funktionalitet förutom Beskrivning av struktur23.
I detta dokument presenterar vi de experimentella protokoll för imaging rhizobox växt rot-system som använder en billig och enkel skräddarsydda RGB imaging installationsprogrammet och (ii) en mer komplex NIR imaging. Exempel resultat från dessa två inställningar är visas och diskuteras i samband med anläggningen fenotypning och fysiologiska växtforskning.
Protocol
1. Rhizoboxes for plant growth NOTE: The experimental system uses rhizoboxes to grow plants for root imaging. First the design of the boxes and the substrate used are described, and then details on the filling procedure are given. Design of rhizoboxes Create rhizoboxes (Figure 2) with a back plate and side frames made of grey PVC with 15 mm strength. Use a box size of 300 mm x 1000 mm. For the front window, use 6 mm mineral glass which is attached to the PVC frame by metal rails being screwed into the side walls. Create three holes on the bottom frame to allow drainage of excess water. These holes can be optionally closed by plastic screws. Before filling, adapt the inner diameter of the rhizobox (between 10 mm and 30 mm) by inserting PC multiwall sheets. An inner space of 10 mm is recommended for most crops to reduce the weight (rhizobox weight without soil is 13.2 kg) of the entire system. Figure 2: Rhizobox experimental system and its components. Left figure shows the dimensions of a rhizobox and right figure its single components, being a grey PVC back plate with side frame, a front mineral glass, multiwall sheets for variable inner diameter and metal angles to fix the front glass to the back compartment. Please click here to view a larger version of this figure. Substrate Fill the rhizoboxes with field soil (for this experiment: silt loam top soil from a calcareous chernozem) sieved to 2 mm particle size. Open the rhizoboxes for filling the substrate into the inner compartment (back plate with side frames) in horizontal position using pre-wetted substrate. Fill horizontally to avoid layering and segregation between fine and coarse particles which occurs when filling in vertical position by pouring the substrate through the upper opening. Pre-wet the substrate before filling. Depending on the type of substrate (particularly on its silt and clay content), do not exceed a water content of 0.12-0.18 cm3*cm-3 to avoid smearing and structure degradation. Add the difference between pre-mixing and target water content after filling the substrate into the rhizoboxes. NOTE: Filling with (oven-)dry substrate and subsequent addition of the entire water is not recommended as is can result in strong settlement of the substrate and formation of large cracks. Step by step filling example Define the target water content. Here it is initially set at 80% plant available water (PAW) where plants do not suffer from any water shortage. Determine field capacity (FC) and permanent wilting point (PWP) of the substrate. Here obtain FC using a PVC tube of equivalent height (100 cm) as the rhizoboxes. Close the tube at the bottom with a stubble having small drainage holes, add 1 cm of gravel to avoid closure of the holes from finer substrate and fill the substrate to the same bulk density as used for the rhizoboxes (1.3 g cm-3). Saturate the tube with water until drainage occurs and leave it for two days to equilibrate (the resulting water content is per definition equivalent to field capacity), while covering the upper opening with cling-film to avoid evaporation. The field capacity value achieved for the soil in this experiment was 0.357 cm3 cm-3. NOTE: Water content at PWP has to be known in advance from standard soil physical methods (e.g. pressure plate measurements24) or from texture based pedotranfer functions25. Here it is equal 0.12 cm3 cm-3 for the soil used. When measuring FC by pressure plate extraction too, take the water content at a matrix potential of h=-100 hPa and not h=-330 hPa in order to correspond to the rhizobox geometry (height of 100 cm = 100 hPa). Calculate the water content (WC) at 80% PAW: WC (cm3 cm-3) = 0.80 (FC-PWP) + PWP. For the hydraulic limits of the soil used here this gives a volumetric water content at 80% PAW of 0.31 cm3 cm-3. Calculate the water amount for a rhizobox volume of 2850 cm3 (30 cm width, 1 cm inner space, 95 cm in height, with top 5 cm keeping free of substrate for watering). This gives a water volume of 883.5 cm3 equal 883.5 g for the density of water being 1.0 g cm-3 at 20 °C. Define the bulk density (db) to fill the rhizoboxes. Here, set at 1.3 g cm-3 corresponding to values typically found in agricultural field soils. The amount of dry substrate required to fill a rhizobox volume of 2850 cm3 at this db equals 3705 g of dry soil. Pre-wet the dry soil to a gravimetric water content of 0.108 g g-1 (equal to a volumetric water content of 0.14 cm3 cm-3) by adding 400 g of water for 3705 g of dry soil and mix it gently to obtain a homogeneous water distribution. Manually disrupt larger aggregates to keep the particle size ≤ 2 mm. Fill the pre-wetted soil into the opened rhizoboxes and compact it gently using a polystyrene sheet (30 x 10 x 1.5 cm) to cover the inner volume of the box, thereby resulting in a homogeneous db of 1.3 g cm-3. Add the remaining amount of water (483.2 g) to achieve the target water content of 0.31 cm3 cm-3 by spraying onto the surface with a spray bottle. Ensure small drop size to avoid surface structure degradation and homogeneous wetting. Keep the box on a balance during spraying to monitor the amount of water actually added to substrate. Let the water redistribute for 10 minutes and then press the glass onto the surface and fix it with the side metal rails. The average final weight of rhizoboxes with wetted substrate was 17818 ± 68 g (13230 g rhizobox weight + 3705 g dry soil + 883 g water). NOTE: The homogeneous water content filled in the horizontal boxes will redistribute when the boxes are set in their final position according to the resulting potential gradient. This is a physical process in all plant growth pots according to their geometry (height) and experimenters should be conscious on their pot hydraulics26. 2. Climate room setup Equip the climate room (Figure 3) with 8 LED lamps providing homogeneous illumination of 450 μmol m-2 s-1 with spectral peaks at 440 (blue) and 660 (red) nm for optimum plant growth. Figure 3: Climate chamber with rhizoboxes for stress experiment. (A) Left view of entire chamber with LED illumination, weather station and PC (here for logging leaf hygrometers); right view with a close-up of the metal frame holding rhizoboxes in 45 ° inclination and wooden plates used to shied rhizobox glass window against light. (B) Stress experiment with sugar beet combining four stages with stress due to different atmospheric demand (high/low) and soil water availability (high/low). Green bars of mean stomata conductance give an indication of plant stress response. Please click here to view a larger version of this figure. Set the ambient parameters according to plant/experimental needs. Here, use 14 hours light and 10 hours dark for illumination. During plant establishment and before stress treatments started, set the temperature to 20° C during day and 15° C during night and keep relative humidity at 50 ± 8%. Put the rhizoboxes at an inclination of 45° using an adequate metal framework. This maximizes root growth towards the glass surface due to gravitropism. Cover the glass window by a wooden plate to keep the root zone in dark and avoid algae growth due to light penetrating through the glass surface. 3. Sugar beet example setup and treatments Pre-germinate sugar beet seeds on wet filter paper for three days at 20 °C in an incubator until the radicle emerges. This ensures seeding with viable plants. NOTE: Pre-germination is not required for plants with high germination vigor, thereby avoiding the risk of damaging the radicle at seeding. For sugar beet seeds with a thick pericarp, however the risk of non-viable seeds is high and pre-germination significantly accelerates the emergence of the radicle compared to direct seeding into soil. Drill a small hole to about 1.5 cm soil depth in the middle of the rhizobox with a screwdriver, position one seed in it using tweezers with the radicle oriented downwards and next to the glass window (this improves the initial visibility) and gently cover it with soil. Add a 0.5 cm layer of fine gravel (2-4 mm) on top of soil to protect the soil aggregates from slaking during irrigation and reduce evaporation losses. To facilitate emergence, keep the soil surface free of gravel where the seed has been positioned. Add 10 g of water to enhance establishment. During establishment and early growth until experimental stress treatments start, irrigate the rhizoboxes every 2-4 days to keep the initial moisture content of 80% PAW. Determine the amount of irrigation water required by weighing the rhizoboxes and adding water until achieving the initial weight of each individual box. For manual irrigation use a pipette to avoid surface structure degradation during watering. Arrange the rhizoboxes according to the established design in the climate room. The protocol reported here is based on an experiment with six sugar beet cultivars in five replicates in a completely randomized design (CRD). Reposition rhizoboxes each time when taking them out of the metal holder for weighing and watering. This avoids any effects of residual (light) inhomogeneity inside the climate room. Define the time for onset of stress treatments and the type of stress. The following settings (Figure 4) are used here. Start measurements at BBCH 15 (five leaves unfolded) with roots covering around 75% of the rhizobox depth and the canopy being sufficiently developed for measurements at the leaves. Keep each stage for at least three days to ensure adaptation to the new settings and make measurement. For the non-stress observations keep the initial settings with optimum soil moisture (80% PAW) and ambient conditions (20° C/15° C temperature, 50 ± 8% rH) giving a daytime vapor pressure deficit (VPD) of 1.28 ± 0.1 kPa (cf. Figure 3 C). Raise the atmospheric demand to a daytime VPD of 2.45 ± 0.4 kPa by increasing the temperature to 27°C/20°C and decreasing rH to 35%, while keeping soil moisture at 80% PAW. Subsequently dry down the rhizoboxes to 40% PAW equivalent to a water content of 0.215 cm3 cm-3 by withholding irrigation. Reset the initial ambient conditions with low atmospheric demand (VPD of 1.28 kPa). Combine stresses increasing the atmospheric demand to a VPD of 2.45 kPa and keep soil moisture at 40% PAW. 4. Root imaging methods Combine imaging methods to make use of their respective advantages and depending on the information targeted. Apply RGB imaging in the VIS range to track root growth, architecture and morphology over time which implicitly requires frequent measurement. Advantages of RGB imaging are (i) low costs, (ii) rapid image acquisition, (iii) low requirements of hard-disk space (image size: 48 MB) and (iv) high resolution (3648 x 5472 pixels). Use hyperspectral imaging (HSI) in the NIR range when chemometric features of root and soil are required. Advantages are (i) spectral features for segmentation between roots and soil background and (ii) access to physico-chemical system properties (e.g. soil water content, root age). Disadvantages are (i) higher scanning time (about 16 minutes per rhizobox), (ii) large size of datasets (13.7 GB per rhizobox image) and (iii) lower resolution of the NIR camera (320 x 256 Pixel), and (iv) higher complexity of data analysis. RGB root imaging Use an imaging box (Figure 4) that shields from ambient light and fixes the camera position consisting of a metal frame with a width of 1 m and a height of 1 m with side walls lined with pressboards. At the front end, fix the camera at two positions at a distance of 80 cm from the rhizobox. Attach a tapeline on the rhizobox frame with transparent adhesive tape and position the rhizobox in the holder of the imaging box. Illuminate the rhizobox using four 24 W fluorescent light tubes attached at a distance of 80 cm from the rhizobox. Also, mount four 15 W UV tubes at 20 cm from the rhizobox as alternative illumination making use of root auto-fluorescence in case of low contrast between root and (bright colored) substrate background. Take two images (top and bottom position) to cover the upper and the lower half of a rhizobox with an overlap of about 3 cm. Acquire RGB images with a digital single-lens reflex camera which is fixed by quick release plates on the respective positions of the imaging box. Apply the following settings when using the fluorescent light tubes. Adapt these example settings for any changes in dimensions and illumination as well as the camera model. Turn off autofocus and stabilizer on the camera objective. Set camera to manual mode. Set ISO speed to 500. Set shutter speed to 13. Set Aperture to 5.6. Turn off the mirror lock. Set the white balance to Auto White Balance. Use the following settings for illumination with UV tubes: Turn off autofocus and stabilizer on the camera objective. Set camera in manual mode. Set ISO speed to 1000. Set shutter speed to 13. Set Aperture to 5.6. Turn off the mirror lock. Set the white balance to Fluorescence. Merge the RGB images from top and bottom of the rhizobox into a single image a photo editor (e.g., Adobe Photoshop). Use the tapeline at the side frame of the rhizoboxes and control overlapping objects (root and soil features). Copy the two separate images, each with pixels size of 3648 x 5472, into a new file of size 3648 x 10944 pixels and white background. Reduce the layer opacity for one image to 60% and align overlapping parts of the images (tapeline, objects). Thereafter restore layer opacity to 100%. Based on the ruler on the image add two red lines of exactly 1 cm to the top of the image where no roots are present. These lines are later used at image analysis to scale the image to correct length dimensions. Merge all layers and remove parts of the image outside the soil filled window with the cropping tool. Save the image as tiff-file for further analysis. In case of images with UV illumination, reduce colors to greyscale before saving using the toolbar Adjustments-Black & White, and selecting the predefined High Contrast Blue filter. For segmentation between roots and soil background as well as for quantification of root traits of interest (e.g. length, surface, diameter, branching) use any root analysis software (see www.plant-image-analysis.org for available tools). Here WinRhizo is used. Open the rhizobox tiff-image the software. Calibrate the length scale of the image using the scaling bars added to the image. Select Based on Color in the Analysis menu where selecting Root & Background Distinction. Define a calibration file with color classes corresponding to roots and soil (background). For the example images used here (cf. Figure 8) three root color classes (old laterals, young laterals, tap root) and three soil color classes are defined. For the greyscale images (UV-illumination) select (i) Based on grey levels, and (ii) Pale Root on Black Background in the Root & Background Distinction menu and use a local Lagarde intensity threshold for segmentation. Open a data file where results are saved. Run the analysis and subsequently control whether there are regions (e.g. at the edges) which are mismatched. In this case define an exclusion region and restart the analysis. For roots not classified, add additional color classes and restart the analysis. For elements wrongly classified as roots, activate/increase the Debris & Rough Edges filtering options. Figure 4: Imaging box to acquire RGB rhizoboxes pictures. Left view of front side where rhizoboxes are attached for imaging with light sources inside; right view of backside where the camera is mounted. Please click here to view a larger version of this figure. Hyperspectral root imaging Hardware setup Use a hyperspectral root imaging system (Figure 5) consisting of (i) a thermo-electrically cooled 14-bit monochrome NIR camera with a spectral range from 900 nm to 1700 nm, 320 by 256 pixels and a frame rate of 100Hz and (ii) an imaging spectrograph with a spectral range 900 nm to 2500 nm and a spectral resolution of 3.6 nm. Arrange a halogen line illumination source (four 50 W halogen spots) in a 45°/-45° geometry. Mount the imaging sensor on a two-axis positioning system. The scan window has a size of 240 x 1000 mm, i.e. 30 mm at each edge of the rhizobox are not covered by the image. Control the system by a Matlab script for (i) white and dark standard acquisition, (ii) setting of camera integration time, (iii) selecting spatial (pixel size 0.1 mm; pixel size 1.0 mm) and spectral resolution (all 222 spectral bands with a resolution of 3.6 nm; smoothed spectrum with 54 bands and a resolution of 14.8 nm), and (iv) defining the scan region on the rhizobox. Save images as SIF-files. To avoid problems during saving of large files, subdivide each scanned image stride (9 strides per rhizobox) into four segments (three of 300 mm length, one of 100 mm length) and save separately with a unique file name consisting of stride number (1 to 9) and part (1 to 4) as well as date and time (YYYY.MM.DD HH:MM:SS). An entire rhizobox scan requires a hard-disk space of 13.7 GB. Image acquisition and analysis. NOTE: Figure 6 shows the steps of image acquisition, segmentation and analysis. Image acquisition comprises selection of camera setting for optimum image quality and definition of scan parameters. Determine the camera integration times for the rhizobox scan and the white standard in the camera software. Open the imaging GUI and move the camera to a position of the rhizobox where roots are present. Adjust the integration time of the camera targeting a light object (i.e. root) in a way that approximately 85% of the full dynamic range of the camera is used on the histogram displayed by the software. Repeat for the white standard by moving the camera positioning system to target the white standard. Then close the camera software. Open the Matlab Imaging GUI and make all settings for the current rhizobox scan. For the data reported here, use the following settings: Integration time white standard: 1000 Integration time rhizobox: 4000 Spectral resolution: Full resolution (i.e. 222 narrow-range spectral bands) Full spatial resolution (pixel size of 0.1 mm) Acquire the dark and white standards before each imaging run, e.g. once a day. The dark standard represents the camera noise, while the white standard gives the maximum reflectivity. These data are required for image normalization during pre-processing. Define whether the entire rhizobox or only part of it is scanned. For the present case entire rhizoboxes are imaged. Then start the scan. Process the image with a Matlab script. Operations performed by the script are described. NOTE: Scripts are currently in an undocumented version and can be obtained from the corresponding author. After proper documentation, they will be available for download from the website of the corresponding author's institution (www.dnw.boku.ac.at/pb/). Compose an entire image stride from the rhizobox center (containing roots) merging the four parts of the stride. NOTE: At this stage it is neither necessary nor recommended to use a spectral image of an entire rhizobox (i.e. all 9 strides) as the file size will make each calculation step in Matlab very time consuming and the information contained in one central stride is sufficient for the first steps of image analysis. Normalize the image using the acquired dark and white standards and taking into account the different integration times of white standard and rhizobox scans which are saved automatically during scanning in a file. Optionally apply a smoothing filter to remove noise from the image. The script currently offers 3×3 kernel median filtering and multiple scatter correction. For the image evaluation presented here, no filters are applied. Display the image at all recorded spectral bands to obtain a first insight and decide a wavelength to be displaying for selecting regions of interest (cf. image segmentation). Perform segmentation between roots and soil background in a separate script with the following steps. Select regions of interest (ROI) for root and soil to find spectral features for segmentation. Use the freehand selection tool to mark a ROI on the image displayed at a wavelength previously identified with good contrast between roots and soil. Here, use three ROIs on the root (old and young laterals, tap root) and two ROIs in the soil (dry, wet region). Display a rectangle with the selected ROIs and the remaining part as a black mask and visualize the selected ROI image at all wavelengths. Remove all lines in the image matrix containing pixels of spectral intensity = 0 (black pixels). Fuse the root and soil ROIs into one foreground (root) and one background (soil) matrix for segmentation. Search spectral bands (intensity of single spectra or spectral ratios) providing the best separation between root foreground and soil background27. Quantify the distinction between the resulting pixel histograms for root and soil using Bhattacharyya distance28. Select a threshold intensity value separating the histograms. Create a binary image by applying the selected threshold to the original image. This sets all pixels having smaller intensity than the threshold to zero and those with higher intensity to one (done automatically be the script). Save the binary image as tiff-file. Open the binary image and analyze root traits. Select (i) Based on grey levels, and (ii) Pale Root on Black Background in the Root & Background Distinction menu and use global intensity threshold (the image is already binarized). For mapping the water content from the hyperspectral data acquire a calibration dataset and apply a calibration equation to a rhizobox image. Subdivide a rhizobox into 5 cm compartments using polystyrene sheets to fill them with soil (same substrate and db as used for the experiment) at different water contents (Figure 8). Calculate the respective amounts of water to be mixed with soil and fill the compartments (same procedure as described in 1.3. for the entire rhizobox). Scan the calibration rhizobox with the same settings as used for the planted rhizoboxes. Perform the following steps using a script. Merge the four parts of a stride from the water calibration box to one stride and normalize it with dark and white standards. Select rectangular boxes at each compartment with different water content and save them in a structure array. Determine the spectral feature that best separates the water content compartments. This is done with a search algorithm for the global maximum of the intensity differences between mean spectra of adjacent water contents. Calculate the mean spectral intensity value for this feature for each water content compartment of the calibration rhizobox. Fit a regression equation through the directly measured water content and the respective spectral quantifier. Apply the regression equation to each pixel not classified as root on a rhizobox image and for the spectral feature (band) determined above which best relates to water content. Figure 5: Hyperspectral root scanner. The main components of the scanner are indicated. The small picture shows the camera during imaging of a rhizobox. Please click here to view a larger version of this figure. Figure 6: Steps in hyperspectral root imaging. Hyperspectral root imaging consists of three mains steps being (i) image acquisition, (ii) image segmentation and (iii) analysis of the spectral data. Please click here to view a larger version of this figure. Figure 7: Rhizobox for water calibration. The rhizobox contains compartments with substrate at different water content which are subdivided by polystyrene sheets. Germination paper at the dry compartments ensures that soil particles do not rinse into neighboring compartments. Please click here to view a larger version of this figure. 5. Application examples NOTE: Quantitative root information is applied in the context of plant phenotyping (cultivar comparison) and for plant physiological research. The following aboveground data are reported to exemplify these cases. Leaf area: Measure leaf area non-destructively at selected stages during the experiment via the length and width of leaves as a proxy. Alternatively canopy images can be used29. At the end of the experiment measure leaf length and width together with area of clipped leaves using a leaf area meter. Calibrate the non-destructive method applying a regression equation to the data pairs. Dry matter: At the end of the experiment, measure aboveground dry matter by clipping the plants with scissors and dry them for 24 hours at 105 °C in an oven. Stomata conductance: Measure stomata conductance with a leaf porometer. Before measurement, keep the device for at least one hour in the climate chamber to allow sensors equilibrate with ambient conditions and calibrated the device each time when ambient settings in the climate chamber are changed. Take measurements from at least three leaves per plant.
Representative Results
Example results are presented for root segmentation based on RGB and HS imaging. For spectral imaging an example of high resolution water mapping is provided. Finally results are shown that demonstrate the scientific context where image based root data are applied. RGB based root measurement Figure 8 shows an RGB root image time series of sugar beet cultivar Ferrara. The images reveal some artefacts from inhomogeneous illumination of the rhizoboxes, with brighter areas along the left side and different brightness at the overlapping area between top and bottom images. Figure 8: Root growth time series from the RGB imaging. Pictures show the sugar beet cultivar Ferrara at different days after sowing (DAS). The images show some artefacts due to non-homogeneous illumination at the left side of the image and between top and bottom images. Scale bars, 2 cm. Please click here to view a larger version of this figure. Figure 9 provides details on root segmentation based on color thresholding for cultivar Ferrara at day after sowing (DAS) 35. As a reference (Figure 9A), a binary image is used where all roots were manually tracked with a Graphic Tablet. The time required for manual tracking of the entire, fully developed, dense sugar beet root system was around four hours. Figure 9B gives a detailed view on a selected area at the top of the image where old lateral roots are present. Here several root axes are not classified by the color threshold. At the bottom (Figure 9C) on the contrary, where white young roots are predominant, the color based segmentation properly classifies all root axes. The binarized root system (Figure 9D) shows a black area at the left side from the illumination artefact which was defined as exclusion region before running quantitative analysis. Figure 9E shows the corresponding pixel histograms of selected features (roots vs. soil) for the red channel of the RGB image from Ferrara at DAS 35. The root pixels (blue color) clearly show three peaks corresponding to bright young laterals, dark old laterals and tap root. The overlap between the old laterals and the soil background is very strong, leading to unclassified root axes (cf. Figure 9B). Figure 9: Root segmentation using a color threshold. (A) Manually segmented root system using a Graphic Tablet, (B) area with poorly segmented old root axes in the top and (C) properly segmented young axes in the bottom of the image, and (D) binary image obtained from color based thresholding. (E) Pixel histograms for selected features of the RGB image. Roots are represented by the blue bars with different root types indicated; soil is represented by the red bars. Scale bars in A and D, 2 cm; scale bars in B and C, 1 cm. Please click here to view a larger version of this figure. The resulting total visible root length quantified for the manually segmented reference image is 1534.1 cm, while the automatized, color based segmentation gives a total root length of 1427.6 cm. Greyscale images from UV-illumination do not provide an advantage in the case shown here and performed worse compared to color thresholding (root length: 1679.7 cm). Old roots could not be segmented, and there was more noise in the image, probably due to lower light intensity of the UV lamps. However, in case of young roots with high auto-fluorescence and a bright background substrate, UV-illumination can still be an option as shown by an image obtained from another experiment where sand was used as background substrate (Figure 10). Figure 10: UV illumination to visualize roots on bright background. Example from a durum wheat root system growing in a rhizobox filled with quartz sand. The rhizobox is imaged with illumination for (A) UV light and (B) fluorescent (day) light. Scale bars, 2 cm. Please click here to view a larger version of this figure. HSI based root measurements Figure 11 provides the mean spectra for three root ROIs (old and young lateral, tap root) and two soil ROIs (top and bottom of rhizobox). Figure 11: Mean spectra of root and soil. Spectra from regions of interest (ROIs) on the root (three root types) and in the soil (top and bottom of the rhizobox). The ROIs are selected to determine an optimum segmentation criterion between root and soil. Please click here to view a larger version of this figure. It is evident that the tap and young lateral roots differ substantially from the background in intensity of most spectral bands. For the old laterals the intensity differences are much lower. A feature that can be inferred visually is the different slope of the spectrum around water absorption region (1450 nm). Here the slope of root spectra is higher compared to soil spectra. Furthermore a change of tap and young lateral spectra in the region around 1100 nm can be identified that does not occur in the old laterals. Figure 12A shows the result from the search algorithm identifying a spectral ratio with strongest foreground-background contrast. The ratio of spectra at 1476 nm to 1076 nm provides the best separation between roots and soil. The resulting histogram of root foreground and soil background pixels is shown in Figure 12B. Although there is some overlap, most pixels are clearly separated from the soil background. Fitting a bimodal Gaussian curve through the histogram and using Bhattacharyya distance for quantification, a value of 7.80 is obtained. A value higher 3.0 indicates strong image contrast allowing reliable separation28. Figure 12: Difference in reflectance between root foreground and soil background for different spectral band ratios and pixel histogram at spectral ratio used for segmentation. (A) Bright colors (yellow) show high contrast between foreground and background, dark colors (blue) show low contrast. The first 15 bands have been removed because of noise. The red lines indicate the band ratio with highest contrast. (B) Pixel histogram of roots (blue) and soil (red) at segmentation spectral ratio. Blue bars represent the root and red bars the soil. The intensity value corresponds to the ratio of spectral band 160 to spectral band 49. Please click here to view a larger version of this figure. The binary image (Figure 13) is created by applying a global intensity threshold of the identified spectral ratio at a value of 1.008 calculated from the histogram distance27. Analysis of root length of this image gives a total length of 1557.3 cm which represents an error of only 1.5% compared to the manually tracked reference image. Figure 13: Binary image of the root system of sugar beet cultivar Ferrara. The image is obtained by applying a global spectral threshold. Scale bar (bottom left corner), 2 cm. Please click here to view a larger version of this figure. Although root segmentation has improved using spectral information compared to color based information, the main intention of HS imaging is analysis of chemometric image properties. This is exemplified via mapping the water content of a rhizobox image. Figure 14A shows the mean spectra of the compartments in the calibration rhizobox (cf. Figure 7) filled with soil of different water content. The shape of spectra is similar between the compartments, i.e. here a spectral ratio does not necessarily provide a more stable classification criterion. Thus intensity at a single spectral band (1680 nm), where the average difference between adjacent water contents is maximized, is identified as best separating criterion. The resulting pixel histograms for this spectral wavelength are shown in Figure 14B. Figure 14: Spectral features for water content calibration. (A) Mean spectra of nine water compartments from the calibration rhizobox with different water contents; (B) Pixel histograms for the water compartments at band 216 where average distance between neighboring compartment is maximum. Please click here to view a larger version of this figure. The relation of the average pixel intensity at 1680 nm and the measured water content is shown in Figure 15. Figure 15: Relation of spectral reflectance and volumetric water content. The figure shows data pairs of measured water content and spectral reflectance with empirical curves (linear and exponential) fit to the data excluding the highest water contents (red triangles). Please click here to view a larger version of this figure. Differentiation of higher water contents from spectral intensity becomes difficult. A significant regression (either linear or exponential) with high R2 can be fit to water contents up to around 0.30 cm3 cm-3. Wetter soil conditions cannot be reliably predicted by the intensity value. Similar behavior of an exponential relation between reflectivity and water content with a decreasing response to water contents higher 0.30 cm3 cm-3 was also found in other studies30. A rhizobox image with fine mapping of water content is shown in Figure 16. Four aspects have to be remarked. First, a region of lower water content can be seen in the rooted parts of the rhizobox. Second, strongest depletion is concentrated in the vicinity to single root axes. Third, depletion zones also occur where no root axes are visible on the surface, indicating regions where roots are hidden in soil. Fourth, water mapping without further image-processing results in a patchy appearance due to the aggregated soil. This can indicate inhomogeneous water content distribution at the aggregate scale, but also surface morphology effect on image quality. Chemometric image-processing techniques are an option to overcome such morphological effects in spectral images31, but are not implemented so far in the Matlab scripts used here. Figure 16: Water content mapping on a rhizobox. The dark blue colours represent regions of high water content, green to red areas show regions with low water content. The plant root is overlaid on the image in black. Scale bar (bottom left corner), 2 cm. Please click here to view a larger version of this figure. Application examples Figure 17 relates quantitative root traits from image analysis with aboveground measurements. Figure 17: Typical application examples for root data. (A) and (B) show root information used for aboveground-belowground plant characterization in a phenotyping context. (A) represents root growth from the sugar beet cultivar Ferrara, (B) compares six rhizobox grown sugar beet cultivars using leaf-to-root area ratio (data from one replicate). (C) and (D) are functional relations between traits as found in plant physiological research. (C) shows the influence of leaf-to-root area ratio on dry matter production and (D) the relation of root surface area to stomata conductance. Please click here to view a larger version of this figure. Figure 17A and 17B are relevant for phenotyping focusing on comprehensive aboveground and belowground plant characterization. Figure 17A shows root growth of sugar beet cultivar Ferrara (cf. Figure 8 for images). Expansion of the root system indicates the capacity of a cultivar to explore the soil volume in a given time span of the vegetation period. Figure 17B shows leaf-to-root surface area ratio of six sugar beet cultivars, providing a descriptor for the balance between plant supply (root) and demand (leaf). Figures 17C and 17D give examples for functional relations of interest in physiological research. In Figure 17C leaf-to-root surface area ratio is related to dry matter formed during the experiment, indicating the predominant role of the assimilating surface as a limiting factor for dry matter accumulation. The lack of significance in spite of a comparatively high R2 is related to the low number of paired data (n=6) used here. Figure 17D reveals that cultivars with higher root surface area (improved uptake) have an average higher stomata conductance over the course of the experiment. The higher root area apparently sustains water extraction, thereby prolonging stomata opening.
Discussion
Protokollen ger två kompletterande metoder för jord odlas rotsystem imaging. Ett kritiskt steg för tillförlitlig experimentella resultat är att fylla i den rhizoboxes som har att säkerställa en jämn och homogen substrat lager på frontglaset att ge snäva rot-smutsa kontakten på fönstret observation och undvika luftspalter. Detta är den främsta anledningen att använda jämförelsevis fina sållen jord av < 2 mm: större aggregat resultera i högre ytan morfologi på fönstret observation med håligheter mellan aggregat. Förutom en högre risk för roten tip uttorkning kräver detta också mer komplex bild bearbetning tekniker för vatten mappning31.
Ändringar av protokollet därför fokusera på förbättrad och snabb fyllning av rhizoboxes. För närvarande påfyllningstid är cirka 30 minuter per låda. Vidare användning av rhizoboxes med två glas windows för avbildning från båda sidor och ändringar att optimera belysning homogenitet för bättre RGB-bilder är testade. Ytterligare maskinvara förlängning också överväga integrering av planar optodes32 samt kapacitans imaging33 in i rhizobox-systemet. Detta är dock utöver nuvarande uppgradering verksamhet.
Programvaruändringar fokusera på automatisk bildregistrering till säkring toppen och botten RBG bilder34. För Hyperspektrala metoder imaging avancerad oövervakade funktion utvinning28 samt känsligare övervakade mål identifieringsmetoder såsom SVMs35 testas. De Hyperspektrala data tillåter därmed potentiellt för bedömning av flera jord, odlingsbädden och rot boenden36. Dessutom är avsikten att utveckla en (halv) automatiserade programvara för rhizobox rot bilder baserat på en modifierad version av rotsystemet Analyzer37 att kvantifiera morfologiska (längd, diameter, yta) samt arkitektoniskt drag (förgrening frekvens, förgrening vinklar).
Den största begränsningen av protokollet jämfört 3D imaging tillvägagångssätt är begränsningen till ytan synliga rot och odlingsbädden boenden. Det har dock visats att de synliga rot drag är en pålitlig proxy för den hela rotsystem21. Den rhizobox tekniken är enkelt kombineras med traditionella destruktiva provtagning (tvättning) i slutet av dynamisk tillväxt imaging för att validera förhållandet av synliga vs. totala rotsystem drag. Som denna relation kan variera bland arter21, rekommenderas destruktiva provtagning att garantera tillförlitlig inferens från synliga egenskaper för de nya fenotypning med en annan gröda arter.
Den viktigaste fördelen med av protokollet som presenteras här är kombinationen av realistiska odlingsförhållanden (jord), relativt höga potentiella genomflödet för temporally löst RGB-imaging och inferens på roten funktionalitet (t.ex. upptag av vatten) via chemometric rot och odlingsbädden data från Hyperspektrala imaging. Därmed metoderna som övervinner inferens begränsningar i hög genomströmning plantan och icke-smutsa rot imaging metoder14, medan det kan delvis djupa fenotypning insikter i funktionella processer med mindre experimentella komplexitet och högre dataflöde jämfört med avancerade 3D metoder15.
I kommande experiment protokollet används för att studera effekten av mykorrhiza på rotsystemet utveckling och funktionalitet av baljväxter samt när det gäller fenotypning rot kännetecken för täcka grödan arter i förhållande till markens struktur, kväve och kol cykling.
Declarações
The authors have nothing to disclose.
Acknowledgements
Författarna erkänner finansiering från österrikiska Science fond FWF via den projekt nummer P 25190-B16 (rötterna torka motstånd). Etablering av Hyperspektrala tänkbar infrastruktur stöddes ekonomiskt av den federala regeringen Niederösterreich (Land Niederösterreich) via projektet K3-F-282/001-2012. Ytterligare finansiering för sockerbetor experimentet mottogs från AGRANA forskning & Innovation Center GmbH (Ulf). Författarna tackar Craig Jackson för teknisk support under experimentet och engelska korrigering av manuskriptet. Vi erkänner också Markus Freudhofmaier som bidragit till etableringen inställningen för RGB-imaging och Josef Schodl för byggandet av rhizobox montering.
Materials
| Rhizobox | Technisches Büro für Bodenkultur | Experimental builder | |
| LED Lamps ATUM HORTI 600 | Klutronic | AHI10600F | |
| Fluorescent light tube HiLite T5 Day | Juwel Aquarium | 86324 | |
| UV light tube Eurolite 45cm slim 15 W | Conrad | 593384 – 62 | |
| Canon EOS 6D | Canon Austria GmbH | 8035B024 | |
| Adobe Photoshop CS5 Extended Version 12.0 x 32 | Adobe Systems Software Ireland Ltd. | ||
| WinRhizo Pro v. 2013 | Regent Instruments Inc. | ||
| Xeva-1.7-320 SWIR camera | Xenics | XEN-000105 | |
| Spectrograph Imspector N25E | Specim | ||
| Hyperspectral imaging scanner | Carinthian Tech Research AG | Experimental builder | Design and assemblage of Hyperspectral Imaging Scanner and software |
| Matlab R2106a | Mathworks | Including Toolboxes for Image Processing, Signal Processing and Statistics and Machine Learning | |
| AP4 Poromoeter | Delta-T-Devices | ||
| LI-3100C Area Meter | LI-COR | ||
| BASF Styradur polystyrene sheets | Obi Baumarkt | 9706318 | Different types of polystyrene sheets or other material separating differently moistured soil can be used |
Referências
- Kutschera, L. . Wurzelatlas mitteleuropäischer Ackerunkräuter und Kulturpflanzen. , (1960).
- Kenrick, P., Strullu-Derrien, C. The origin and early evolution of roots. Plant Physiol. 166, 570-580 (2014).
- Franklin, J., Serra-Diaz, J. M., Syphard, A. D., Regan, H. M. Global change and terrestrial plant community dynamics. PNAS. 113, 3725-3734 (2016).
- MacDonald, G. K., Bennett, E. M., Potter, P. A., Ramankutty, N. Agronomic phosphorus imbalances across the world’s croplands. PNAS. 108, 3086-3091 (2011).
- Vörösmarty, C. J., Green, P., Salisbury, J., Lammers, R. B. Global water resources: vulnerability from climate change and population growth. Science. 289, 284 (2000).
- Lobell, D. B., Cassman, K. G., Field, C. B. Crop yield gaps: their importance, magnitudes, and causes. Annu. Rev. Environ. Resour. 34, 179-204 (2009).
- Angus, J. F., Van Herwaarden, A. F. Increasing water use and water use efficiency in dryland wheat. Agron. J. 93, 290-298 (2001).
- Lassaletta, L., Billen, G., Grizzetti, B., Anglade, J., Garnier, J. 50 year trends in nitrogen use efficiency of world cropping systems: the relationship between yield and nitrogen input to cropland. Environ. Res. Lett. 9 (10), 105011 (2014).
- Lynch, J. P. Roots of the second green revolution. Aust. J Bot. 55, 493-512 (2007).
- Comas, L. H., Becker, S. R., Von Mark, V. C., Byrne, P. F., Dierig, D. A. Root traits contributing to plant productivity under drought. Front. Plant Sci. 4, 442 (2013).
- Metcalfe, D. B., et al. A method for extracting plant roots from soil which facilitates rapid sample processing without compromising measurement accuracy. New Phytol. 174, 697-703 (2007).
- Kutschera, L., Lichtenegger, E., Sobotik, M. . Wurzelatlas der Kulturpflanzen gemäßigter Gebiete mit Arten des Feldgemüsebaues. , (2009).
- Gioia, T., et al. GrowScreen-PaGe, a non-invasive, high-throughput phenotyping system based on germination paper to quantify crop phenotypic diversity and plasticity of root traits under varying nutrient supply. Funct. Plant Biol. 44, 76-93 (2017).
- Watt, M., et al. A rapid, controlled-environment seedling root screen for wheat correlates well with rooting depths at vegetative, but not reproductive, stages at two field sites. Ann. Bot. 112, 447-455 (2013).
- Fiorani, F., Schurr, U. Future scenarios for plant phenotyping. Annu. Rev Plant Biol. 64, 267-291 (2013).
- Metzner, R., et al. Direct comparison of MRI and X-ray CT technologies for 3D imaging of root systems in soil: potential and challenges for root trait quantification. Plant Methods. 11, 1 (2015).
- Choat, B., Badel, E., Burlett, R., Delzon, S., Cochard, H., Jansen, S. Non-invasive measurement of vulnerability to drought induced embolism by X-ray microtomography. Plant Physiol. 170, 273-282 (2015).
- Metzner, R., et al. Direct comparison of MRI and X-ray CT technologies for 3D imaging of root systems in soil: potential and challenges for root trait quantification. Plant Methods. 11, 1 (2015).
- 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. 10, 1 (2014).
- Bengough, A. G., et al. Gel observation chamber for rapid screening of root traits in cereal seedlings. Plant Soil. 262, 63-70 (2004).
- Nagel, K. A., 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. 39, 891-904 (2012).
- Price, A. H., et al. Upland rice grown in soil-filled chambers and exposed to contrasting water-deficit regimes: I. Root distribution, water use and plant water status. Field Crops Res. 76, 11-24 (2002).
- Vadez, V. Root hydraulics: the forgotten side of roots in drought adaptation. Field Crops Res. 165, 15-24 (2014).
- Dane, J. H., Hopmans, J. W., Dane, J. H., Topp, G. C. Pressure plate extractor. Methods of soil analysis. Part 4. Physical methods. , 688-690 (2002).
- Wösten, J. H. M., Pachepsky, Y. A., Rawls, W. J. Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics. J. Hydrol. 251, 123-150 (2001).
- Passioura, J. B. The perils of pot experiments. Funct. Plant Biol. 33, 1075-1079 (2006).
- Dorrepaal, R., Malegori, C., Gowen, A. Tutorial: Time series hyperspectral image analysis. J. Near Infrared Spec. 24, 89-107 (2016).
- Kim, D. M., et al. Highly sensitive image-derived indices of water-stressed plants using hyperspectral imaging in SWIR and histogram analysis. Scie rep. 5, (2015).
- Humplík, J. F., Lazár, D., Husičková, A., Spíchal, L. Automated phenotyping of plant shoots using imaging methods for analysis of plant stress responses-a review. Plant Methods. 11, 29 (2015).
- Lobell, D. B., Asner, G. P. Moisture effects on soil reflectance. Soil Sci. Soc. Am. J. 66, 722-727 (2002).
- Esquerre, C., Gowen, A. A., Burger, J., Downey, G., O’Donnell, C. P. Suppressing sample morphology effects in near infrared spectral imaging using chemometric data pre-treatments. Chemometr. Intell. Lab. 117, 129-137 (2012).
- Williams, P. N., et al. Localized flux maxima of arsenic, lead, and iron around root apices in flooded lowland rice. Environ Sci. Technol. 48, 8498-8506 (2014).
- Cseresnyés, I., Takács, T., Végh, K. R., Anton, A., Rajkai, K. Electrical impedance and capacitance method: a new approach for detection of functional aspects of arbuscular mycorrhizal colonization in maize. Eur. J Soil Biol. 54, 25-31 (2013).
- Brown, M., Lowe, D. G. Automatic panoramic image stitching using invariant features. Int. J. Comput. Vision. 74, 59-73 (2007).
- Jiang, Y., Li, C., Takeda, F. Nondestructive detection and quantification of blueberry bruising using near-infrared (NIR) hyperspectral reflectance imaging. Scientific Reports. 6, 35679 (2016).
- Chang, C. -. W., Laird, D. A., Mausbach, M. J., Hurburgh, C. R. Near-infrared reflectance spectroscopy-principal components regression analyses of soil properties. Soil Sci. Soc. Am. J. 65, 480-490 (2001).
- Leitner, D., Felderer, B., Vontobel, P., Schnepf, A. Recovering root system traits using image analysis exemplified by two-dimensional neutron radiography images of lupine. Plant Physiol. 164 (1), 24-35 (2014).
Citar este artigo
Bodner, G., Alsalem, M., Nakhforoosh, A., Arnold, T., Leitner, D. RGB and Spectral Root Imaging for Plant Phenotyping and Physiological Research: Experimental Setup and Imaging Protocols. J. Vis. Exp. (126), e56251, doi:10.3791/56251 (2017).