Imaging Techniques Plants Computational Biology

scientificprotocols authored about 8 years ago

Authors: Takumi Higaki, Natsumaro Kutsuna & Seiichiro Hasezawa


The protocol reported here describes a method to quantitatively evaluate fluorescently-tagged protein localizations from fluorescent microscopic images with a combination of probabilistic mapping and image clustering. We demonstrate the use of this protocol using kidney-shaped guard cells of plants.


Microscopic assessment of protein localizations with fluorescent protein tagging is an important assay in cell biological studies. We believe that protein localizations should be statistically assessed based on many observations, because the interpretation of a single microscopic image often leads to biased information. As proved in previous magnetic resonance imaging (MRI) analyses (1,2,3), probabilistic mapping is a simple and powerful method to determine average distributions of visualized objects. Clustering is a kind of unsupervised learning based on distances between each data point. Similarities between microscopic images can be evaluated using clustering with appropriate image metrics (4). Using the fluorescent image of kidney-shaped guard cells, we recently developed the image processing framework to quantitatively evaluate protein localizations with probabilistic mapping and image clustering (5). This protocol explains how to evaluate the localization similarity of marker proteins from the 3-D images of Arabidopsis guard cells that are freely available in our image database LIPS (5). We prefer self-organizing mapping (SOM) (6) as an image clustering method because of its superiority for image inspection.



A laptop computer (CF-J9, Panasonic) with 64-bit Windows operating system, 6 GiB of random access memory (RAM), and the Intel Core i5 CPU was used here.

ImageJ software

Download and install the ImageJ software ( For general operations, please refer the official user guide7

Download Kbi plug-ins and Hig macros

  1. Download our in-house KBI JAR file (‘kbiijplugins…’) and scala JAR file (‘scala-library’) in the plugin folder of ImageJ
  2. Download our in-house HIG IJM files in the macro folder of ImageJ
  3. Restart the ImageJ.


Image processing for probabilistic mapping

  1. Prepare the target stack image. Here, we prepare the 16-bit TIFF stack images (bright field and fluorescence images) that are freely available in the LIPS database (Figure 1).
  2. Run the ImageJ software
  3. Open the target image using the ImageJ menu “File-Open”.
  4. Draw the major stomatal pore line on the bright field image using the ImageJ tool bar “Straight Line Selection Tool” (Figure 2a).
  5. Copy the segmented straight line from the bright field image to the fluorescence images using the ImageJ menu “Edit-Selection-Restore Selection” (Figure 2b).
  6. Rotate the fluorescence image using the ImageJ menu “Plugins-kbi-Kbi_registration (mode: horizoner)” (Figure 2c).
  7. Get the maximum intensity projection image using the ImageJ menu “Image-Stacks-Z project…” (Figure 2d).
  8. Put a bounding box around the cell using the ImageJ tool bar menu “Rectangular Selection Tool” (Figure 2e).
  9. Crop the boxed cell region using the ImageJ menu “Image-Crop” (Figure 2f).
  10. If the is pore located in the upper half of the cropped image, rotate the image 180 degrees using the ImageJ menu “Image-Transform-Rotate…”
  11. Resize the images to the same size using the ImageJ menu “Image-Adjust-Size…” Here, we resized all images to 304×119 pixels, which was the mean size of the guard-cell bounding boxes.
  12. Make a stack image of the grouped images using the ImageJ menu “Plugins-kbi-Kbi_StkJoin (mode: all axis: z)”. (Figure 2g).
  13. Run the macro hig_Pdmapping.ijm You will see the probability map (Figure 2h).
  14. Save the probability map image using the ImageJ menu “File-Save”.

Measurement of raster-scanned intensity profiles

  1. Run the ‘ImageJ’ software.
  2. Open the stack image of the probability maps using the ImageJ menu “File-Open”.
  3. To reduce the data size, shrink the image size using the ImageJ menu “Image-Adjust-Size…”. Here, we reduced the size to 70×27 pixels.
  4. Run the macro hig_Lineprofile.ijm
  5. You will see the Results table of the raster-scanned intensity profiles (Figure 3).
  6. Uncheck the checkbox ‘Save column headers’ and ‘Save row numbers’ using the Results table menu “Results-Options…”.
  7. Save the raster-scanned intensity profiles as a csv file using the Results table menu “File-Save As…”.

SOM clustering with ImageJ

  1. Run the ‘ImageJ’ software.
  2. Open the stack image of the probability maps using the ImageJ menu “File-Open”.
  3. Run the plug-in using the ImageJ menu “Plugins-kbi-Kbi_Clustering” and select ‘som’ in the modeClustering drop-down list.
  4. Select ‘sampleZfeatureCsvpatchImg’ in the somMode drop-down list.
  5. Load the csv file of the raster-scanned intensity profiles.
  6. You will see the SOM image (Figure 4). In the case that multiple images were allocated in same node in SOM, the images will be allocated in next slice of the stack image (SOM output) (Figure 4).


Display error message ‘Out of memory’:

Allocate more memory using the ImageJ menu “Edit-Options-Memory & Threads…-Maximum memory”

Display error message ‘Stack required’:

The macros higPdmapping.ijm and higLineprofile.ijm do not support single-slice images. Please prepare a stack image.

Display error message ‘This macro does not work with RGB images’:

The macros higPdmapping.ijm and higLineprofile.ijm do not support RGB images. Please prepare gray-scale images.


  1. Penhune, V.B., Zatorre, R.J., MacDonald, J.D. & Evans, A.C. Interhemispheric anatomical differences in human primary auditory cortex: probabilistic mapping and volume measurement from magnetic resonance scans. Cereb Cortex, 6, 661-672 (1996).
  2. McCarthy, G., Luby, M., Gore, J. & Goldman-Rakic, P. Infrequent events transiently activate human prefrontal and parietal cortex as measured by functional MRI. J Neurophysiol. 77, 1630-1634 (1997).
  3. Jones, D.K. et al. Spatial normalization and averaging of diffusion tensor MRI data sets. Neuroimage, 17, 592-617 (2002).
  4. Higaki, T., Kutsuna, N., Sano, T., Kondo, N., Hasezawa, S. Quantification and cluster analysis of actin cytoskeletal structures in plant cells: role of actin bundling in stomatal movement during diurnal cycles in Arabidopsis guard cells. Plant J. 61, 156-165 (2010).
  5. Higaki, T. et al. Statistical organelle dissection of Arabidopsis guard cells using image database LIPS. Sci. Rep. in press.
  6. Kohonen, T. Self-Organizing Maps. Springer-Verlag, ISBN 3-540-67921-9, New York, Berlin, Heidelberg (2001).
  7. ImageJ User Guide website. Available: Accessed 2012 May 1.


This study was supported by a Grant-in-Aid for Research Activity start-up to T.H. (No. 22870004); by a Grant-in-Aid for Young Scientists (B) to N.K. (No. 24770038) from the Japanese Ministry of Education, Culture, Sports, Science and Technology; and Grants-in-Aid from Development of Systems and Technology for Advanced Measurement and Analysis, Japan Science and Technology Agency to S.H.


Figure 1: Examples of target images.

Download Figure 1

The 3-D images of Arabidopsis guard cells are freely available from LIP database.

Figure 2: Workflow of image processing for probability mapping.

Download Figure 2

(a) An original bright-field image with a stomatal pore line. (b) An original green-fluorescence protein (GFP) stack image with a stomatal pore line. (c) A rotated GFP stack image. (d) A rotated maximum intensity projection (MIP) image. (e) Selection of the cell region with a bounding box. (f) A cropped guard cell MIP image. (g) Multiple cropped MIP images. Note that all images are the same size. (h) A probabilistic map.

Figure 3: Screenshot of output of the macro hig_Pdmapping.ijm

Download Figure 3

Figure 4: An example of SOM output.

Download Figure 4

The images allocated in the same node are displayed in the same positions of different slices.

Associated Publications

Statistical organelle dissection of Arabidopsis guard cells using image database LIPS. Takumi Higaki, Natsumaro Kutsuna, Yoichiroh Hosokawa, Kae Akita, Kazuo Ebine, Takashi Ueda, Noriaki Kondo, and Seiichiro Hasezawa. Scientific Reports 2 () 11/05/2012 doi:10.1038/srep00405

Author information

Takumi Higaki, Natsumaro Kutsuna & Seiichiro Hasezawa, Graduate School of Frontier Sciences, The University of Tokyo

Correspondence to: Takumi Higaki ([email protected])

Source: Protocol Exchange (2012) doi:10.1038/protex.2012.013. Originally published online 21 May 2012.

Average rating 0 ratings