Skip to main content
  • Methodology article
  • Open access
  • Published:

The RootScope: a simple high-throughput screening system for quantitating gene expression dynamics in plant roots

Abstract

Background

High temperature stress responses are vital for plant survival. The mechanisms that plants use to sense high temperatures are only partially understood and involve multiple sensing and signaling pathways. Here we describe the development of the RootScope, an automated microscopy system for quantitating heat shock responses in plant roots.

Results

The promoter of Hsp17.6 was used to build a Hsp17.6p:GFP transcriptional reporter that is induced by heat shock in Arabidopsis. An automated fluorescence microscopy system which enables multiple roots to be imaged in rapid succession was used to quantitate Hsp17.6p:GFP response dynamics. Hsp17.6p:GFP signal increased with temperature increases from 28°C to 37°C. At 40°C the kinetics and localization of the response are markedly different from those at 37°C. This suggests that different mechanisms mediate heat shock responses above and below 37°C. Finally, we demonstrate that Hsp17.6p:GFP expression exhibits wave like dynamics in growing roots.

Conclusions

The RootScope system is a simple and powerful platform for investigating the heat shock response in plants.

Background

A systems level understanding of plant responses to the environment requires quantitative gene expression data with high spatial and temporal resolution from multiple individuals [1, 2]. This type of data cannot be readily obtained using conventional gene expression assays. RT-PCR, microarray, and RNAseq approaches are limited by low spatial resolution unless tissues are either dissected manually or sorted using FACS. While in situ and GUS reporter based approaches can provide high spatial resolution, they require sacrificing the sampled tissue and thus cannot be used to easily monitor gene expression dynamics in individual plants over time. Both Luciferase and fluorescent protein (FP) based reporters expressed in living tissues can provide dynamic expression information, but only FP based reporters allow for both high spatial and temporal resolution.

Several systems for assaying gene expression dynamics in plant roots using FP reporters have previously been developed. All imaging approaches involve tradeoffs between spatial resolution, temporal resolution, and the number of plants that can be imaged at a time. Cost and ease of use are also important considerations. Light sheet microscopy has been used to generate high resolution 3D expression data in an Arabidopsis root at 10 minute intervals, but only one root can be imaged at a time [3]. Two systems have been developed for imaging FPs in multiple roots in rapid succession. The RootChip is a microfluidics based system which allows extended imaging of eight (and potentially more) roots and allows rapid exchanges of the growth media [4, 5]. A second microfluidics device, the RootArray, uses confocal microscopy to image up to 64 roots in three dimensions and at high resolution. However, long acquisition times in this system limit the temporal resolution; imaging even a subset of the roots at high resolution takes several hours [6]. While extremely powerful, both of these systems require plants to be grown in the microfluidic device so that the roots grow close to a cover slip which enables high resolution imaging. This is a labor-intensive procedure and requires specialized and relatively expensive microfluidics devices. We are interested in using gene expression dynamics as a phenotype in a genetic screen for high temperature stress 'thermostat’ genes. Because the screen involves characterizing large numbers of individuals we wanted to develop a system which maximizes the number of plants imaged simultaneously while retaining high spatial and temporal resolution and providing accurate control of plant temperature.

The motivation for this screen is to understand the genetic network(s) involved in heat shock (HS) sensing and signal transduction. Plants are sessile organisms and must endure widely varying environmental conditions, many of which result in stress and reductions in yield in agricultural systems. Heat is an abiotic stress that disrupts protein and membrane homeostasis at a cellular level. Heat stress is often accompanied by drought stress which exacerbates its deleterious effects. Global climate change is predicted to result in more frequent and severe episodes of both of these stresses [7, 8]. Plants have evolved a complex set of mechanisms for dealing with both stressful as well as non-stressful high temperatures [9–12]. At stressful temperatures the heat shock response (HSR) results in global changes in gene regulation. An important component of the HSR is the synthesis of Heat Shock Proteins (HSPs) which protect normal protein and membrane function. HSP function has been extensively characterized biochemically and genetically and is reasonably well understood [11, 13]. Much less is known about the high temperature sensing and signal transduction mechanisms, the 'thermostat’, which regulates the HSR [9].

In Arabidopsis (as in all organisms) HSP expression is strongly induced in response to heat stress. HSP regulation is mediated by the activity of a family of temperature and cAMP regulated Ca2+ channels (CNGC2,4, and 6), a calmodulin (CAM3), a kinase (CBK3) and a phosphatase (PP7) [14–19]. These genes are predicted to function upstream of a family of highly redundant Heat Shock transcription Factors (HSFs) which directly regulate much, but not all, HSP expression [20–23]. HSF activity is also thought to be regulated through direct physical interactions with the HSPs HSP70 and HSP90 and a HS induced HSF binding protein (HSBP), creating feedback loops which finely tune the HSR [24–27]. Notably, none of these thermostat genes were identified based on their mutant phenotypes, perhaps due to redundancy in the HSR network. Instead, these genes were identified based on homology with other systems, expression profiling, or biochemical interactions and they have all been characterized using reverse genetics approaches.

Many details of the HS thermostat’s function are complex and not fully understood. The CNGC and HSF gene families exemplify this complexity. Based on studies in mammalian systems, CNGCs are thought to function as hetero-tetramers and the HSFs as homo- and hetero-trimers [14, 27, 28]. The in vivo subunit composition of these complexes has not been described in plants. Since the function of these genes is likely influenced by their co-complex interaction partners the analysis of their mutant phenotypes is complicated. Both the CNGC and HSF families contain members which, when mutated, exhibit opposite phenotypes. cngc2 and 4 mutants increase Ca2+ influx into cells and result in elevated thermotolerance while cngc6 mutants decrease Ca2+ influx and exhibit decreased thermotolerance [14, 17]. Similarly, HSFA1a,b,d, and e are positive regulators of the HSR and are required for thermotolerance while HSFB1 and HSFB2b attenuate the HSR [20, 21, 23]. Not only are the physical interactions among Arabidopsis family members not well understood, but their interactions with and regulation by other molecules are also poorly described. For instance, the CNGCs are activated by the second messenger cAMP yet only a few genes involved in cAMP production have been identified in Arabidopsis. None have demonstrated biological functions in plants [14, 29, 30]. Given these ambiguities it is likely that multiple components of the HSR network remain unidentified.

Unbiased forward genetics approaches based on thermotolerance defects have identified only a few HSR genes including HSP101 (HOT1), a chitinase (HOT2), and an exportin (HIT2) [31–33]. None of these genes have obvious sensing or signaling functions and, to our knowledge, no thermostat genes have been identified by a high temperature thermotolerance screen. This suggests that screening for thermotolerance defects using conventional assays may be of limited value in identifying HS thermostat genes.

Forward genetic approaches have failed to identify the signaling and sensing components of the HSR. This may be attributed to: 1) redundancy in the network of genes that makes up the thermostat, and 2) phenotyping approaches which assay the end point of the HSR as opposed to its dynamics. Evidence for the masking effects of redundancy and the limitations of conventional phenotyping approaches comes from studies on the HSF transcription factors. To uncover thermotolerance phenotypes due to a loss of HSFA1 function, at least three family members must be knocked out [21, 23]. Similarly, hsfb1 single mutants have no HSP expression or thermotolerance phenotypes while hsfb1/hsfb2b double mutants have both [20]. Figure 1 illustrates why phenotypic assays with detection thresholds, especially assays which take place at a single late time point, are likely to overlook the effects of genes with redundant family members (Figure 1).

Figure 1
figure 1

A model of how overlapping functions in a gene family might obscure thermotolerance phenotypes. Single and double mutants of the HSFA1 transcription factors have no thermotolerance phenotypes [21, 23], presumably because they still cross a response threshold before the phenotype is assayed. Triple and quadruple mutant phenotypes are detected because they do not cross the threshold. A system that allows the kinetics of the response to be monitored could uncover these single and double mutant phenotypes.

To overcome these limitations we have developed the RootScope which generates the quantitative data required to identify new and characterize existing components of the HSR thermostat. The RootScope consists of a robust HSR reporter and an automated imaging system for monitoring the reporter in multiple plants with high temporal and spatial resolution.

Methods

Construction of the HSP17.6p:GFP reporter line

The Hsp17.6 C-CI (AT1G53540) promoter was amplified from genomic DNA using primers AT1G53540-FWD (5′-CACCATTCAGGTAATCAGGTTGTCTGC-3′) and AT1G53540-REV (5′-CGTTTCACTTCCTCTTGTGATTGC-3′). These primers amplify 713 bp starting in the 3′ UTR of the upstream gene AT1G53530 and ending at the ATG start codon of AT1G53540. CACC was added to the 5′ end of AT1G53540-FWD to provide directionality for TOPO cloning into pENTR/D-TOPO (Invitrogen, Grand Island, NY). TOPO cloning of the PCR product into pENTR/D-TOPO was followed by Gateway mediated LR recombination (Invitrogen) into destination vector pFAST-R07 (Shimada et al., 2010) to create a C-terminal promoter:GFP fusion (HSP17.6p:GFP). This construct was transformed into Col-O plants using Agrobacterium tumefaciens strain GV3101 (Clough and Bent, 1998). Primary transformants were identified using the RFP seed marker contained in the pFAST-R07 construct and were allowed to self-fertilize to generate homozygous lines. Homozygous lines were screened for GFP expression at basal temperatures (22°C) and after a 1 hour 37°C heat shock as 3 day old seedlings. Fluorescence imaging was performed using a Leica MZ16F microscope equipped with a 2x Planapo lens, a GFP3 filter set, and an EL6000 light source (Leica Microsystems).

Construction of the automated imaging system

The heated stage consists of an Inheco CPAC Ultraflat HT microplate heater with a temperature range of 4°C to 110°C (Inheco). Microplates were placed on the CPAC and held in place by a custom 3D-printed chamber which incorporates the heated glass plate from a Tokai Hit INU-GSI stage top incubator. For horizontal root growth the heated stage was mounted on an ASI LS-50 linear stage which provided Z (focus) movement. The LS-50 linear stage was then mounted on an ASI MS-2000 XY stage to create a heated XYZ stage. For vertical root growth a heated XY stage was built by directly mounting the heated stage on the MS-2000 XY stage. All ASI stages were controlled by a LX-4000 control unit (ASI).

Images were acquired using a Hamamatsu ORCA II BT-512 camera cooled to -55°C in both horizontal and vertical arrangements. Exposure times ranged from 75–200 msec depending on the plant line being imaged. For horizontal root growth the heated XYZ stage and camera were mounted under a Leica MZ16F microscope equipped with a 2x Planapo lens, a GFP3 filter set, and an EL6000 light source (Leica Microsystems). For vertical root growth the heated XY stage was mounted vertically on a 1.5″ mounting post (Thorlabs). A Leica Z16 APO macroscope with a 2x Planapo lens and a 1x tube lens (Leica) was coupled to a DFM dichroic filter cube with a MDF-GFP filter set (Thorlabs) and the ORCA II camera using custom adapters. Illumination was provided by a M470-L2 high power LED collimated using an ACL2520-A lens (Thorlabs). These components were mounted on a 8″ x 24″ x 1/2″ aluminum breadboard which in turn was mounted on a Parker Daedal CR4900-8 crossed roller bearing stage (Parker Daedal). The LS-50 linear stage was connected to the breadboard to provide focus control by moving the entire microscope relative to the XY stage. All of these components were built up on a 12″ x 36″ optical breadboard (Vere).

Both the LED and EL6000 light sources were shuttered using a USB connected Arduino UNO driving a SainSmart 2-Channel relay (SainSmart). For the EL6000 light source the relay was directly connected to the light source shutter inputs. For the high power LED light source the relay was connected to the DIM input of a Luxdrive 1400 mA Buckblock (LUXdrive) which was used to drive the LED.

μManager microscopy automation software [34] running on a 64 bit Windows 7 (Microsoft) system was used to control the Arduino based shutter, the LX-4000 stage control unit, and the ORCA II camera. Root images were quantitated and, in some cases straightened, using custom ImageJ [35] and Perl scripts. Alignment of straightened roots was performed using the Template Matching and Slice Alignment ImageJ plugin [36].

Plant growth

Arabidopsis seeds were surface sterilized with 30% bleach and 0.1% SDS for 15 minutes, rinsed three times in sterile water, and plated on OmniTray single well microplates (Thermo Scientific product #242811) containing 60 ml of 0.5x MS media. The plated seeds were stratified for three days at 4°C and then transferred to Percival E-30B growth chambers (Percival Scientific) set at 22°C with 24 hours of light for three days.

Results

Generation of a Hsp17.6p:GFP reporter for the heat shock response

Assaying the HSR in multiple plants with high spatial and temporal resolution requires a live cell reporter with a high dynamic range and a strong signal when induced. We chose Hsp17.6 (AT1G53540) because this small HSP is consistently among the most highly induced HSPs in both our own and published HS microarray experiments [21, 37, 38]. We cloned the AT1G53540 promoter upstream of GFP in the pFAST-R07 vector [39] and transformed the resulting Hsp17.6p:GFP construct into Arabidopsis. We used this vector because it contains a RFP seed marker which allows for easy identification of plants which contain the reporter. This feature will be useful for identifying plants homozygous for the reporter in the future. We selected 46 independent transformants and allowed each to self-fertilize to generate plants homozygous for the reporter. Homozygous lines were screened for GFP expression in 3 day old roots at basal temperatures and an hour after a one hour 37°C heat shock. None of the lines exhibited GFP expression at basal temperatures, consistent with a lack of Hsp17.6 expression in the absence of cellular stress. The intensity of the GFP signal after the heat shock varied among lines and we chose three of the brightest lines (8, 11, and 26) for further characterization (data not shown).

A simple high-throughput screening system for quantitating gene expression dynamics in plant roots

To monitor the Hsp17.6p:GFP reporter we built an automated fluorescence microscopy system controlled by μManager software [34]. The microscope has a temperature-controlled chamber consisting of a Peltier heating/cooling device that accepts a petri dish with the same size format as a 96 well plate. Because condensation on the lid of the petri dish would prevent imaging we incorporated a heated glass lid into the chamber (Figure 2A-F). The chamber was mounted on a motorized XY stage to allow multiple roots to be imaged in rapid succession. The chamber can be used horizontally on the stage of a conventional upright fluorescence microscope (Figure 2G). In this orientation roots grow horizontally across the surface of the media. The chamber can also be mounted vertically on a custom fluorescence microscope (Figure 2H-I), which allows the roots to grow vertically on the surface of the media. Vertical growth is advantageous because, when grown horizontally, about half of the roots stop growing horizontally, attempt to grow into the media, and thus cannot be imaged.

Figure 2
figure 2

An overview of RootScope construction. A heating/cooling block is used to provide temperature control (A). A chamber is placed on the heating/cooling block (B) and a microplate containing the plants to be imaged is placed in the chamber (C). A gasket is placed on top of the chamber (D) allowing the heated glass plate to sit on top of and seal the chamber (E). The resulting stack of components (side view, F) can then be mounted either horizontally (G) or vertically (H-I) on an automated microscope.

The advantage of this setup compared to a microfluidic device is that we can easily and quickly plate 100 (or more) seeds on MS media on a single, cheap, and commercially available microplate. Seed placement is not critical as the roots do not have to grow in defined micro-channels. Root positions are identified manually using the μManager slide explorer plugin once the plate has been transferred to the microscope for imaging. The seeds are germinated in an incubator, grown for 3 or 4 days, and then transferred to the microscope where they can be simultaneously heat shocked and imaged. The system allows us to quantitate the HSR at high temporal (every 4 minutes with 100 plants) and spatial (~2.4μm with our current camera and optics) resolution with a 12 bit dynamic range.

Characterization of the induction and attenuation kinetics of the Hsp17.6p:GFP reporter

Seeds homozygous for Hsp17.6p:GFP from lines 8, 11, and 26 were plated on the same plate, germinated at 22°C, and at four days placed in the growth chamber on the microscope and heat shocked continuously at 37°C. The roots were imaged every four minutes for 800 minutes. Induction was detectable at one hour and peaked at approximately eight hours in all lines. The intensity of the GFP signal varied between lines; line 26 was the brightest and line 8 was the dimmest (Figure 3A). To determine the variability of the induction kinetics of these three lines we normalized their maximum fluorescence intensity. All three lines exhibit very similar kinetics (Figure 3B) suggesting that line 26, while brighter than the other lines we screened, exhibits similar kinetics to the other lines and is not brightest because of an abnormal expression pattern. Because it has the highest dynamic range and normal induction kinetics, line 26 was used exclusively from this point on and all references to Hsp17.6p:GFP refer to line 26.

Figure 3
figure 3

Comparison of Hsp17.6 p :GFP expression levels and dynamics in three independent transformants. Four day old seedlings from lines 8 (green), 11 (red), and 26 (blue) homozygous for Hsp17.6p:GFP were heat shocked at 37°C starting at time point 0 and imaged every four minutes. Induction of the reporter over 13:20 is shown without (A) and with normalization (B). Error bars are +/- standard deviation.

To characterize the attenuation of Hsp17.6 expression after a short heat shock we heat shocked Hsp17.6p:GFP for two hours and then reduced the temperature to 22°C, imaging the plants every 4 minutes. This HS regimen results in the induction of Hsp17.6p:GFP at approximately one hour (starting half way through the heat shock) and attenuation at five hours. We performed this analysis with plants both homozygous and heterozygous for Hsp17.6p:GFP to determine if the dosage of the reporter affected the intensity or the kinetics of the reporter. Although the induction kinetics for heterozygotes and homozygotes were similar, plants containing a single copy of the reporter were on average only 62% as bright as plants homozygous for the reporter and attenuation occurred earlier in heterozygotes (Figure 4). This result highlights the importance of ensuring matched copy numbers when using the Hsp17.6p:GFP reporter to characterize the HSR.

Figure 4
figure 4

Induction kinetics of Hsp17.6 p :GFP homozygous and heterozygous plants. Plants homozygous (blue) and heterozygous (red) for Hsp17.6p:GFP were heat shocked for two hours at 37°C and imaged every four minutes over a 12:00 period. Error bars are +/- standard deviation.

Hsp17.6p:GFP expression correlates with the heat shock temperature

To determine if the reporter could be used as a read out of root temperature we exposed plants to heat shock temperatures between 28 and 40°C for one hour. Consistent with existing array data for Hsp17.6, very little response was detected below 30°C, a temperature at which Arabidopsis root growth is largely inhibited. The magnitude of the induction is linear from 31–37 degrees (R2 = 0.9977) and the response over the temperature range from 28–37 degrees can be modeled with a 2° order polynomial (R2 = 0.984). This demonstrates that the reporter can be used to accurately report high temperatures in the root (Figure 5). The maximal response was observed at 37°C, a result similar to those of a previous experiment in which acquired thermotolerance (AT) was maximally induced at 37°C [22]. These data suggest that Hsp17.6 is induced by a HS thermostat that is distinct from the thermostat which function at non-stressful temperatures (<28°C) to regulate growth and development [12, 40].

Figure 5
figure 5

Induction kinetics of Hsp17.6 p :GFP in response to heat shock temperatures from 28°C to 40°C. Plants were heat shocked for one hour at temperatures from 28-40°C and imaged every four minutes for at least 11:30 (A). The average intensity of the GFP signal from time points between 2:00–2:20 (1:00–1:20 after the end of heat shock) are plotted with either a linear regression (31-37°C, red) or a 2° order polynomial regression (28-37°C) (B). Error bars are +/- standard deviation.

Both the kinetics and the localization of the Hsp17.6p:GFP response are markedly different at 40°C compared to lower temperatures. At 40°C induction does not occur until after the heat shock ends, suggesting that at temperatures above 37°C the HSR machinery itself may be inhibited (Figure 5). Again, this result is consistent with the observation that acclimation temperatures above 37°C result in lower levels of AT [22]. In addition to differences in the kinetics of the response, the localization of the response maximum is different below and above 37°C. Below 37°C the response maximum occurs in the differentiation zone at approximately 5 hours. In contrast, at 40°C the maximum signal occurs after 10 hours and is initially located at the root tip and subsequently in the root tip and the vasculature (Figure 6 and Additional file 1). The marked difference in kinetics and localization supports the idea that multiple HS thermostats exist in plants and demonstrates the temporal and spatial resolution of our system.

Figure 6
figure 6

The location of Hsp17.6 p :GFP expression maxima is temperature dependent. Hourly images of representative roots heat shocked for an hour at 28-37°C are shown. Arrowheads indicate the location of the brightest signal for the first two time points where the maximum is clear (either the differentiation zone in roots heat shocked from 28-37°C or the root tip in roots heat shocked at 40°C). Arrows indicate late vascular expression in roots heat shocked at 40°C. The contrast of all images in each temperature series was uniformly enhanced using the ImageJ 'enhance contrast’ tool so that both low and high expression patterns can be seen. As a result, intensity should not be compared between temperature series. Quantification of signal intensity can be found in Figure 5.

Hsp17.6p:GFP expression exhibits wave like dynamics

Several thousand genes expressed in the root (3493) have been shown to oscillate in two phases with a periodicity of ~6 hours in the meristematic and elongation zones. These oscillations were characterized at relatively low spatial and temporal resolution using luciferase based markers and microarray data and are associated with developmental events including root bending and branching. The periodicity of these oscillations was stable over a temperature range of 18-24°C, however the effects of stressful high temperatures (>30°C) on these oscillations were not characterized [41]. Plants analyzed in these studies were grown under normal conditions and thus did not express the HSR genes, which are expressed only in response to cellular stress. To determine if Hsp17.6 expression oscillates we examined the time lapse series from Hsp17.6p:GFP roots which had been heat shocked for one hour at 37°C. After induction, Hsp17.6p:GFP exhibited wave like expression dynamics in the epidermis, sweeping from the differentiation zone through the elongation zone to the meristem and back again. The periodicity of these waves was 3.9 +/- 0.8 hours (n = 15). The waves were associated with developmental events, sweeping in front of cell differentiation as visualized by the appearance of root hairs on epidermal cells (Figure 7 and Additional file 2). Because Hsp17.6 is not expressed at basal temperatures it is unclear whether the differences in periodicity between the previously described genes and Hsp17.6p:GFP are due to multiple oscillating circuits of gene expression in the root, to the effects of the heat shock on the wavelength of the oscillations, or due to differences between the folding and degradation rates of luciferase and GFP.

Figure 7
figure 7

Oscillating waves of Hsp17.6 p :GFP expression during root growth after a 1 hour 37°C heat shock. A kymograph display of the position of maximal Hsp17.6p:GFP expression in the epidermis (A). Images of a straightened and aligned root exhibiting Hsp17.6p:GFP oscillations at 20 minute intervals (B).

Discussion

The RootScope is a simple and flexible system that makes it possible to determine the dynamics of gene expression in large numbers of plant roots with high temporal and spatial resolution. The time and cost required to set up plants for screening are minimized since seeds are plated on standard media in microplates and plating seeds for an experiment takes only minutes. The most time consuming step in a RootScope experiment is the manual identification of root positions once the plate is transferred to the microscope. This step takes less than 20 minutes per plate for an experienced undergraduate student. It should be possible to automate this task using image recognition algorithms; however, since automated localization may need to be confirmed and fine-tuned manually, the advantages of automating this step may not be significant in practice. In this paper we report using the RootScope as a fluorescence microscope for quantifying a GFP based transcriptional reporter. Similar to existing non-fluorescent phenotyping systems such as Phytomorph, the microscope is flexible and could also be used to generate high resolution incident light images to characterize growth dynamics in multiple plants simultaneously [42, 43]. Because the optics we have used are based on a zoom macroscope both larger and smaller fields of view can be selected (with concomitant changes in resolution). Finally, although we have chosen to image roots, there is no technical reason why the aerial portions of plants could not be imaged using this system.

The induction kinetics of Hsp17.6 in this system is slower than previous reports of HSP protein and transcript accumulation in which HSPs accumulate within half an hour of a heat shock [44, 45]. The sensitivity of our system and the maturation kinetics of GFP are both likely to contribute to this difference. Both RT-PCR and enzymatic Western blot detection techniques amplify small signals significantly making it possible to detect small numbers of molecules while our imaging approach does not. This limits the sensitivity of our system at very low Hsp17.6 expression levels. The folding, chromophore maturation, and degradation dynamics of GFP are also likely to influence the kinetics reported here, most likely by introducing delays in both the appearance and attenuation of the GFP signal. Hsp17.6p:GFP was built using pFAST-R07 which contains EGFP [39, 46]. EGFP is a GFP variant which exhibits approximately a four fold increase in the rate of chromophore maturation compared to wild type GFP. GFP variants with even faster maturation rates such as GFPmut2 or GFPm have been demonstrated to mature over twice as fast as EGFP in vitro[47] although no in planta maturation rate data exist in the literature. While the use of a rapidly maturing GFP variant has the potential to decrease the delay between transcription and detection of the GFP signal using the RootScope, there is no reporter that we are aware of which can completely eliminate this delay. Since we are not concerned with the absolute timing of these events and given the relatively rapid expression dynamics we have observed, we do not anticipate that these artifacts will prohibit the use of this system for identifying HSR thermostat mutants. In addition to using a rapidly maturing GFP to decrease delays in the onset of the signal, the use of a translational as opposed to a transcriptional reporter might result in more accurate degradation kinetics.

The ability to monitor the kinetics of the HSR using Hsp17.6p:GFP should allow us to take a forward genetics approach to identifying HSR genes based on changes in Hsp17.6p:GFP induction and decay kinetics. We have EMS mutagenized 3,000 Hsp17.6p:GFP seeds and harvested pools of 4–5 self-pollinated plants. 80 seeds from each pool will be screened on plates with 20 non-mutagenized seeds as controls. This pooled approach should allow us to identify most of the mutants in each family and strikes a balance between throughput (each run takes 12 hours) and coverage. Assuming 200 days of screening in the first year we could screen 800–1000 families which, while it would not reach saturation, should uncover several HSR thermostat mutants [48]. We will inspect the data both manually and computationally (based on extracting the times at which each plant’s HSR signal rises above a baseline and reaches a maximum and on the slope of the induction and attenuation phases) to identify mutants that have altered HSR kinetics, localization, or respond as if they are hotter or colder than they actually are. In addition to identifying HSR thermostat mutants we also anticipate discovering mutations which disrupt the wave like expression dynamics of Hsp17.6p:GFP. These mutants may provide insights into the co-ordination of developmental events and stress responses. This system will also be useful for uncovering the effects of previously identified HSR genes (many of which do not display single mutant phenotypes) on HSR dynamics.

Conclusions

We have demonstrated that the RootScope can be used to generate gene expression data with high temporal and spatial resolution in multiple plants simultaneously. We have used the system to uncover evidence that supports the hypothesis that multiple distinct high temperature responses exist in plants. Additionally, we have documented previously unknown waves of HS induced gene expression. Apart from characterizing the HSR, the RootScope should be of general use for quantitating other types of gene expression dynamics. Using other types of reporters [2, 49], it should also be useful for quantitating metabolite and small molecule dynamics in plant roots.

References

  1. Megason SG, Fraser SE: Imaging in systems biology. Cell. 2007, 130 (5): 784-795. 10.1016/j.cell.2007.08.031.

    Article  PubMed  CAS  Google Scholar 

  2. Wee CW, Dinneny JR: Tools for high-spatial and temporal-resolution analysis of environmental responses in plants. Biotechnology letters. 2010, 32 (10): 1361-1371. 10.1007/s10529-010-0307-8.

    Article  PubMed  CAS  Google Scholar 

  3. Sena G, Frentz Z, Birnbaum KD, Leibler S: Quantitation of cellular dynamics in growing Arabidopsis roots with light sheet microscopy. PLoS ONE. 2011, 6 (6): e21303-10.1371/journal.pone.0021303.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Grossmann G, Guo WJ, Ehrhardt DW, Frommer WB, Sit RV, Quake SR, Meier M: The RootChip: an integrated microfluidic chip for plant science. The Plant cell. 2011, 23 (12): 4234-4240. 10.1105/tpc.111.092577.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Grossmann G, Meier M, Cartwright HN, Sosso D, Quake SR, Ehrhardt DW, Frommer WB: Time-lapse fluorescence imaging of Arabidopsis root growth with rapid manipulation of the root environment using the RootChip. J Vis Exp. 2012, 65: e4290-10.3791/4290

    Google Scholar 

  6. Busch W, Moore BT, Martsberger B, Mace DL, Twigg RW, Jung J, Pruteanu-Malinici I, Kennedy SJ, Fricke GK, Clark RL, et al: A microfluidic device and computational platform for high-throughput live imaging of gene expression. Nature methods. 2012, 9 (11): 1101-1106. 10.1038/nmeth.2185.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Ahuja I, de Vos RCH, Bones AM, Hall RD: Plant molecular stress responses face climate change. Trends Plant Sci. 2010, 15 (12): 664-674. 10.1016/j.tplants.2010.08.002.

    Article  PubMed  CAS  Google Scholar 

  8. Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, Kitzberger T, Rigling A, Breshears DD, Hogg EH, et al: A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For Ecol Manage. 2010, 259 (4): 660-684. 10.1016/j.foreco.2009.09.001.

    Article  Google Scholar 

  9. McClung CR, Davis SJ: Ambient thermometers in plants: from physiological outputs towards mechanisms of thermal sensing. Curr Biol. 2010, 20 (24): R1086-R1092. 10.1016/j.cub.2010.10.035.

    Article  PubMed  CAS  Google Scholar 

  10. Yeh CH, Kaplinsky NJ, Hu C, Charng YY: Some like it hot, some like it warm: phenotyping to explore thermotolerance diversity. Plant Sci : Int J Exp Plant Biol. 2012, 195: 10-23.

    Article  CAS  Google Scholar 

  11. Kotak S, Larkindale J, Lee U, von Koskull-Döring P, Vierling E, Scharf K-D: Complexity of the heat stress response in plants. Curr Opin Plant Biol. 2007, 10 (3): 310-316. 10.1016/j.pbi.2007.04.011.

    Article  PubMed  CAS  Google Scholar 

  12. Kumar SV, Wigge PA: H2A.Z-Containing Nucleosomes Mediate the Thermosensory Response in Arabidopsis. Cell. 2010, 140 (1): 136-147. 10.1016/j.cell.2009.11.006.

    Article  PubMed  CAS  Google Scholar 

  13. Wang W, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9 (5): 244-252. 10.1016/j.tplants.2004.03.006.

    Article  PubMed  CAS  Google Scholar 

  14. Gao F, Han X, Wu J, Zheng S, Shang Z, Sun D, Zhou R, Li B: A heat-activated calcium-permeable channel–Arabidopsis cyclic nucleotide-gated ion channel 6–is involved in heat shock responses. Plant J. 2012, 70 (6): 1056-1069. 10.1111/j.1365-313X.2012.04969.x.

    Article  PubMed  CAS  Google Scholar 

  15. Liu H-T, Sun D-Y, Zhou R-G: Ca2+ and AtCaM3 are involved in the expression of heat shock protein gene in Arabidopsis. Plant Cell Environ. 2005, 28 (10): 1276-1284. 10.1111/j.1365-3040.2005.01365.x.

    Article  CAS  Google Scholar 

  16. Liu H-T, Li G-L, Chang HUI, Sun D-Y, Zhou R-G, Li B: Calmodulin-binding protein phosphatase PP7 is involved in thermotolerance in Arabidopsis. Plant Cell Environ. 2007, 30 (2): 156-164. 10.1111/j.1365-3040.2006.01613.x.

    Article  PubMed  CAS  Google Scholar 

  17. Finka A, Cuendet AF, Maathuis FJ, Saidi Y, Goloubinoff P: Plasma membrane cyclic nucleotide gated calcium channels control land plant thermal sensing and acquired thermotolerance. The Plant Cell. 2012, 24 (8): 3333-3348. 10.1105/tpc.112.095844.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Zhang W, Zhou RG, Gao YJ, Zheng SZ, Xu P, Zhang SQ, Sun DY: Molecular and genetic evidence for the key role of AtCaM3 in heat-shock signal transduction in arabidopsis. Plant Physiol. 2009, 149 (4): 1773-1784. 10.1104/pp.108.133744.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Liu HT, Gao F, Li GL, Han JL, Liu DL, Sun DY, Zhou RG: The calmodulin-binding protein kinase 3 is part of heat-shock signal transduction in arabidopsis thaliana. Plant J. 2008, 55 (5): 760-773. 10.1111/j.1365-313X.2008.03544.x.

    Article  PubMed  CAS  Google Scholar 

  20. Ikeda M, Mitsuda N, Ohme-Takagi M: Arabidopsis HsfB1 and HsfB2b Act as repressors of the expression of heat-inducible hsfs but positively regulate the acquired thermotolerance. Plant Physiol. 2011, 157 (3): 1243-1254. 10.1104/pp.111.179036.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Liu HC, Liao HT, Charng YY: The role of class A1 heat shock factors (HSFA1s) in response to heat and other stresses in arabidopsis. Plant Cell Environ. 2011, 34 (5): 738-751. 10.1111/j.1365-3040.2011.02278.x.

    Article  PubMed  CAS  Google Scholar 

  22. Liu HC, Charng YY: Acquired thermotolerance independent of heat shock factor A1 (HsfA1), the master regulator of the heat stress response. Plant Signal Behav. 2012, 7 (5): 547-550. 10.4161/psb.19803.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Yoshida T, Ohama N, Nakajima J, Kidokoro S, Mizoi J, Nakashima K, Maruyama K, Kim J-M, Seki M, Todaka D, et al: Arabidopsis HsfA1 transcription factors function as the main positive regulators in heat shock-responsive gene expression. Mol Genet Genomics. 2011, 286 (5–6): 321-332.

    Article  PubMed  CAS  Google Scholar 

  24. Baler R, Welch W, Voellmy R: Heat shock gene regulation by nascent polypeptides and denatured proteins: hsp70 as a potential autoregulatory factor. J Cell Biol. 1992, 117 (6): 1151-1159. 10.1083/jcb.117.6.1151.

    Article  PubMed  CAS  Google Scholar 

  25. Abravaya K, Myers MP, Murphy SP, Morimoto RI: The human heat shock protein hsp70 interacts with HSF, the transcription factor that regulates heat shock gene expression. Genes Dev. 1992, 6 (7): 1153-1164. 10.1101/gad.6.7.1153.

    Article  PubMed  CAS  Google Scholar 

  26. Hahn A, Bublak D, Schleiff E, Scharf K-D: Crosstalk between Hsp90 and Hsp70 chaperones and heat stress transcription factors in tomato. The Plant Cell Online. 2011, 23 (2): 741-755. 10.1105/tpc.110.076018.

    Article  CAS  Google Scholar 

  27. Hsu S-F, Lai H-C, Jinn T-L: Cytosol-localized heat shock factor-binding protein, AtHSBP, functions as a negative regulator of heat shock response by translocation to the nucleus and is required for seed development in arabidopsis. Plant Physiol. 2010, 153 (2): 773-784. 10.1104/pp.109.151225.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Scharf KD, Heider H, Hohfeld I, Lyck R, Schmidt E, Nover L: The tomato Hsf system: HsfA2 needs interaction with HsfA1 for efficient nuclear import and may be localized in cytoplasmic heat stress granules. Mol Cell Biol. 1998, 18 (4): 2240-2251.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Wong A, Gehring C: Computational Identification of Candidate Nucleotide Cyclases in Higher Plants. Cyclic Nucleotide Signaling in Plants. Edited by: Gehring C. New York, NY, USA: Humana Press; 2013:1016: 195-205. 10.1007/978-1-62703-441-8_13.

    Chapter  Google Scholar 

  30. Ruzvidzo O, Dikobe B, Kawadza D, Mabadahanye G, Chatukuta P, Kwezi L: Recombinant Expression and Functional Testing of Candidate Adenylate Cyclase Domains. Cyclic Nucleotide Signaling in Plants. Edited by: Gehring C. New York, NY, USA: Humana Press; 2013:1016: 13-25. 10.1007/978-1-62703-441-8_2.

    Chapter  Google Scholar 

  31. Hong S-W, Vierling E: Hsp101 is necessary for heat tolerance but dispensable for development and germination in the absence of stress. Plant J. 2001, 27 (1): 25-35. 10.1046/j.1365-313x.2001.01066.x.

    Article  PubMed  CAS  Google Scholar 

  32. Kwon Y, Kim S-H, Jung M-S, Kim M-S, Oh J-E, Ju H-W, Kim K-i, Vierling E, Lee H, Hong S-W: Arabidopsis hot2 encodes an endochitinase-like protein that is essential for tolerance to heat, salt and drought stresses. Plant J. 2007, 49 (2): 184-193. 10.1111/j.1365-313X.2006.02950.x.

    Article  PubMed  CAS  Google Scholar 

  33. Wu SJ, Wang LC, Yeh CH, Lu CA, Wu SJ: Isolation and characterization of the arabidopsis heat-intolerant 2 (hit2) mutant reveal the essential role of the nuclear export receptor EXPORTIN1A (XPO1A) in plant heat tolerance. New Phytol. 2010, 186 (4): 833-842. 10.1111/j.1469-8137.2010.03225.x.

    Article  PubMed  CAS  Google Scholar 

  34. Edelstein A, Amodaj N, Hoover K, Vale R, Stuurman N: Current Protocols in Molecular Biology. Computer Control of Microscopes Using μManager. Hoboken, NJ, USA: John Wiley & Sons, Inc; 2010.

    Google Scholar 

  35. Schneider CA, Rasband WS, Eliceiri KW: NIH image to ImageJ: 25 years of image analysis. Nature methods. 2012, 9 (7): 671-675. 10.1038/nmeth.2089.

    Article  PubMed  CAS  Google Scholar 

  36. Template Matching and Slice Alignment--- ImageJ Plugins. https://sites.google.com/site/qingzongtseng/template-matching-ij-plugin,

  37. Busch W, Wunderlich M, Schöffl F: Identification of novel heat shock factor-dependent genes and biochemical pathways in arabidopsis thaliana. Plant J. 2005, 41 (1): 1-14.

    Article  PubMed  CAS  Google Scholar 

  38. Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K: The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007, 50 (2): 347-363. 10.1111/j.1365-313X.2007.03052.x.

    Article  PubMed  CAS  Google Scholar 

  39. Shimada TL, Shimada T, Hara-Nishimura I: A rapid and non-destructive screenable marker, FAST, for identifying transformed seeds of arabidopsis thaliana. Plant J. 2010, 61 (3): 519-528. 10.1111/j.1365-313X.2009.04060.x.

    Article  PubMed  CAS  Google Scholar 

  40. Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136( ((1): 2621-2632.

    Article  Google Scholar 

  41. Moreno-Risueno MA, Van Norman JM, Moreno A, Zhang J, Ahnert SE, Benfey PN: Oscillating gene expression determines competence for periodic arabidopsis root branching. Science. 2010, 329 (5997): 1306-1311. 10.1126/science.1191937.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. Miller ND, Parks BM, Spalding EP: Computer-vision analysis of seedling responses to light and gravity. Plant J. 2007, 52 (2): 374-381. 10.1111/j.1365-313X.2007.03237.x.

    Article  PubMed  CAS  Google Scholar 

  43. Subramanian R, Spalding E, Ferrier N: A high throughput robot system for machine vision based plant phenotype studies. Mach Vis Appl. 2013, 24 (3): 619-636. 10.1007/s00138-012-0434-4.

    Article  Google Scholar 

  44. Charng YY, Liu HC, Liu NY, Hsu FC, Ko SS: Arabidopsis Hsa32, a novel heat shock protein, is essential for acquired thermotolerance during long recovery after acclimation. Plant Physiol. 2006, 140 (4): 1297-1305. 10.1104/pp.105.074898.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Swindell WR, Huebner M, Weber AP: Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics. 2007, 8: 125-10.1186/1471-2164-8-125.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Karimi M, Inze D, Depicker A: GATEWAY vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 2002, 7 (5): 193-195. 10.1016/S1360-1385(02)02251-3.

    Article  PubMed  CAS  Google Scholar 

  47. Iizuka R, Yamagishi-Shirasaki M, Funatsu T: Kinetic study of de novo chromophore maturation of fluorescent proteins. Anal Biochem. 2011, 414 (2): 173-178. 10.1016/j.ab.2011.03.036.

    Article  PubMed  CAS  Google Scholar 

  48. Page DR, Grossniklaus U: The art and design of genetic screens: arabidopsis thaliana. Nat Rev Genet. 2002, 3 (2): 124-136. 10.1038/nrg730.

    Article  PubMed  CAS  Google Scholar 

  49. Okumoto S, Jones A, Frommer WB: Quantitative imaging with fluorescent biosensors. Annu Rev Plant Biol. 2012, 63: 663-706. 10.1146/annurev-arplant-042110-103745.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank John Kelly and Steve Palmer for their support in building the RootScope and Bill Gresh for providing plant care. Catherine Crouch and Carl Grossman provided valuable advice about microscope optics. We would like to thank two anonymous reviewers for their comments and suggestions which helped us to improve this paper. The research reported in this publication was supported by the Friedman Fund for Physiology at Swarthmore College and by the National Institute of General Medical Sciences of the National Institutes of Health under award number 1R15GM093933. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health or of Swarthmore College.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nicholas J Kaplinsky.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

NK, EK, DN, and RL designed experiments. NK built the automated microscope. RL and CR built the Hsp17.6p:GFP reporter and transformed it into plants. DN and CR characterized transformed plant lines and generated homozygous stocks. NK, EK, and DN performed and analyzed RootScope experiments. All authors helped to draft and approved the final manuscript.

Erin J Kast, Minh-Duyen T Nguyen, Rosalie E Lawrence contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

12870_2013_1364_MOESM1_ESM.avi

Additional file 1: Induction kinetics of Hsp17.6 p :GFP in response to heat shock temperatures at 37°C and 40°C. Roots heat shocked for one hour at either 37°C (top) or 40°C (bottom) were imaged every four minutes for 11 hours. (AVI 2 MB)

12870_2013_1364_MOESM2_ESM.avi

Additional file 2: Oscillating waves of epidermal Hsp17.6 p :GFP expression during root growth after a 1 hour 37°C heat shock. Roots were heat shocked for one hour and imaged every four minutes for 11 hours. Root images were computationally straightened and aligned to allow analysis of the position of the wave. (AVI 568 KB)

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Kast, E.J., Nguyen, MD.T., Lawrence, R.E. et al. The RootScope: a simple high-throughput screening system for quantitating gene expression dynamics in plant roots. BMC Plant Biol 13, 158 (2013). https://doi.org/10.1186/1471-2229-13-158

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-13-158

Keywords