This is a protocol preprint. Preprints are preliminary reports that have not undergone peer review. They should not be considered conclusive, used to inform clinical practice, or referenced by the media as validated information.
Method Article
Probabilistic mapping and image clustering for quantitative assessment of fluorescent protein localizations in Arabidopsis guard cells
Takumi Higaki
Graduate School of Frontier Sciences, The University of Tokyo
Corresponding Author
License:
This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License
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.
Keywords
Image processing, Microscopy, Protein localization, Organelle distribution, Plant stomata
Figure 1
Figure 2
Figure 3
Figure 4
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) analyses1,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 metrics4. 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 clustering5. 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":http://hasezawa.ib.k.u-tokyo.ac.jp/lips/ 5. We prefer self-organizing mapping (SOM)6 as an image clustering method because of its superiority for image inspection.
Computer
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 (http://rsbweb.nih.gov/ij/download.html). For general operations, please refer the official user guide7 "http://rsbweb.nih.gov/ij/docs/guide/userguide.html":http://rsbweb.nih.gov/ij/docs/guide/userguide.html.
Download Kbi plug-ins and Hig macros
- Download our in-house KBI JAR file ('kbi_ij_plugins...') and scala JAR file ('scala-library') in the plugin folder of ImageJ "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/ImageJKbiPlugins":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/ImageJKbiPlugins.
- Download our in-house HIG IJM files in the macro folder of ImageJ "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering.
- Restart the ImageJ.
Image processing for probabilistic mapping
- 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":http://hasezawa.ib.k.u-tokyo.ac.jp/lips/ (Figure 1).
- Run the ImageJ software "http://rsbweb.nih.gov/ij/":http://rsbweb.nih.gov/ij/.
- Open the target image using the ImageJ menu “File-Open”. "http://rsbweb.nih.gov/ij/docs/guide/userguide-23.html#toc-Subsection-23.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-23.html#toc-Subsection-23.2.
- Draw the major stomatal pore line on the bright field image using the ImageJ tool bar “Straight Line Selection Tool” "http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.2 (Figure 2a).
- Copy the segmented straight line from the bright field image to the fluorescence images using the ImageJ menu “Edit-Selection-Restore Selection” "http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.12":http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.12 (Figure 2b).
- Rotate the fluorescence image using the ImageJ menu “Plugins-kbi-Kbi_registration (mode: horizoner)” (Figure 2c).
- Get the maximum intensity projection image using the ImageJ menu “Image-Stacks-Z project…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.6":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.6 (Figure 2d).
- Put a bounding box around the cell using the ImageJ tool bar menu “Rectangular Selection Tool” "http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.1":http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.1 (Figure 2e).
- Crop the boxed cell region using the ImageJ menu “Image-Crop” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.8":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.8 (Figure 2f).
- 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…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.12":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.12.
- Resize the images to the same size using the ImageJ menu “Image-Adjust-Size…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.2. Here, we resized all images to 304×119 pixels, which was the mean size of the guard-cell bounding boxes.
- Make a stack image of the grouped images using the ImageJ menu “Plugins-kbi-Kbi_StkJoin (mode: all axis: z)”. (Figure 2g).
- Run the macro hig_Pdmapping.ijm "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering. You will see the probability map (Figure 2h).
Save the probability map image using the ImageJ menu “File-Save”.
Measurement of raster-scanned intensity profiles
- Run the ‘ImageJ’ software.
- Open the stack image of the probability maps using the ImageJ menu “File-Open”.
- 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.
- Run the macro hig_Lineprofile.ijm "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering.
- You will see the Results table of the raster-scanned intensity profiles (Figure 3).
- Uncheck the checkbox ‘Save column headers’ and ‘Save row numbers’ using the Results table menu “Results-Options…”.
Save the raster-scanned intensity profiles as a csv file using the Results table menu “File-Save As…”.
SOM clustering with ImageJ
- Run the ‘ImageJ’ software.
- Open the stack image of the probability maps using the ImageJ menu “File-Open”.
- Run the plug-in using the ImageJ menu “Plugins-kbi-Kbi_Clustering” and select ‘som’ in the modeClustering drop-down list.
- Select ‘sampleZ_featureCsv_patchImg’ in the somMode drop-down list.
- Load the csv file of the raster-scanned intensity profiles.
- 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” "http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.13":http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.13
Display error message ‘Stack required’:
The macros hig_Pdmapping.ijm and hig_Lineprofile.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 hig_Pdmapping.ijm and hig_Lineprofile.ijm do not support RGB images. Please prepare gray-scale images.
- 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).
- 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).
- Jones, D.K. et al. Spatial normalization and averaging of diffusion tensor MRI data sets. Neuroimage, 17, 592-617 (2002).
- 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).
- Higaki, T. et al. Statistical organelle dissection of Arabidopsis guard cells using image database LIPS. Sci. Rep. in press.
- Kohonen, T. Self-Organizing Maps. Springer-Verlag, ISBN 3-540-67921-9, New York, Berlin, Heidelberg (2001).
- ImageJ User Guide website. Available: "http://rsbweb.nih.gov/ij/docs/guide/userguide.html":http://rsbweb.nih.gov/ij/docs/guide/userguide.html. 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.
Version History
CURRENT STATUS: Posted
Version 1
Posted 21 May, 2012
No community comments so far
Metrics
Comments: 0
HTML Views: ...
Subject Areas
Cell biology
Computational biology and bioinformatics
Biological techniques
Plant sciences
Associated Publications
Statistical organelle dissection of Arabidopsis guard cells using image database LIPS
by Takumi Higaki, Natsumaro Kutsuna, Yoichiroh Hosokawa, +5
Scientific Reports (21 May, 2012)
Method Article
Probabilistic mapping and image clustering for quantitative assessment of fluorescent protein localizations in Arabidopsis guard cells
Takumi Higaki
Graduate School of Frontier Sciences, The University of Tokyo
Corresponding Author
Version History
CURRENT STATUS: Posted
Version 1
Posted 21 May, 2012
No community comments so far
Metrics
Comments: 0
HTML Views: ...
License:
This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License
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.
Figure 1
Figure 2
Figure 3
Figure 4
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) analyses1,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 metrics4. 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 clustering5. 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":http://hasezawa.ib.k.u-tokyo.ac.jp/lips/ 5. We prefer self-organizing mapping (SOM)6 as an image clustering method because of its superiority for image inspection.
Computer
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 (http://rsbweb.nih.gov/ij/download.html). For general operations, please refer the official user guide7 "http://rsbweb.nih.gov/ij/docs/guide/userguide.html":http://rsbweb.nih.gov/ij/docs/guide/userguide.html.
Download Kbi plug-ins and Hig macros
- Download our in-house KBI JAR file ('kbi_ij_plugins...') and scala JAR file ('scala-library') in the plugin folder of ImageJ "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/ImageJKbiPlugins":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/ImageJKbiPlugins.
- Download our in-house HIG IJM files in the macro folder of ImageJ "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering.
- Restart the ImageJ.
Image processing for probabilistic mapping
- 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":http://hasezawa.ib.k.u-tokyo.ac.jp/lips/ (Figure 1).
- Run the ImageJ software "http://rsbweb.nih.gov/ij/":http://rsbweb.nih.gov/ij/.
- Open the target image using the ImageJ menu “File-Open”. "http://rsbweb.nih.gov/ij/docs/guide/userguide-23.html#toc-Subsection-23.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-23.html#toc-Subsection-23.2.
- Draw the major stomatal pore line on the bright field image using the ImageJ tool bar “Straight Line Selection Tool” "http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.2 (Figure 2a).
- Copy the segmented straight line from the bright field image to the fluorescence images using the ImageJ menu “Edit-Selection-Restore Selection” "http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.12":http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.12 (Figure 2b).
- Rotate the fluorescence image using the ImageJ menu “Plugins-kbi-Kbi_registration (mode: horizoner)” (Figure 2c).
- Get the maximum intensity projection image using the ImageJ menu “Image-Stacks-Z project…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.6":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.6 (Figure 2d).
- Put a bounding box around the cell using the ImageJ tool bar menu “Rectangular Selection Tool” "http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.1":http://rsbweb.nih.gov/ij/docs/guide/userguide-18.html#toc-Subsection-18.1 (Figure 2e).
- Crop the boxed cell region using the ImageJ menu “Image-Crop” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.8":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.8 (Figure 2f).
- 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…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.12":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.12.
- Resize the images to the same size using the ImageJ menu “Image-Adjust-Size…” "http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.2":http://rsbweb.nih.gov/ij/docs/guide/userguide-25.html#toc-Subsection-25.2. Here, we resized all images to 304×119 pixels, which was the mean size of the guard-cell bounding boxes.
- Make a stack image of the grouped images using the ImageJ menu “Plugins-kbi-Kbi_StkJoin (mode: all axis: z)”. (Figure 2g).
- Run the macro hig_Pdmapping.ijm "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering. You will see the probability map (Figure 2h).
Save the probability map image using the ImageJ menu “File-Save”.
Measurement of raster-scanned intensity profiles
- Run the ‘ImageJ’ software.
- Open the stack image of the probability maps using the ImageJ menu “File-Open”.
- 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.
- Run the macro hig_Lineprofile.ijm "http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering":http://hasezawa.ib.k.u-tokyo.ac.jp/zp/Kbi/HigPDMClustering.
- You will see the Results table of the raster-scanned intensity profiles (Figure 3).
- Uncheck the checkbox ‘Save column headers’ and ‘Save row numbers’ using the Results table menu “Results-Options…”.
Save the raster-scanned intensity profiles as a csv file using the Results table menu “File-Save As…”.
SOM clustering with ImageJ
- Run the ‘ImageJ’ software.
- Open the stack image of the probability maps using the ImageJ menu “File-Open”.
- Run the plug-in using the ImageJ menu “Plugins-kbi-Kbi_Clustering” and select ‘som’ in the modeClustering drop-down list.
- Select ‘sampleZ_featureCsv_patchImg’ in the somMode drop-down list.
- Load the csv file of the raster-scanned intensity profiles.
- 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” "http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.13":http://rsbweb.nih.gov/ij/docs/guide/userguide-24.html#toc-Subsection-24.13
Display error message ‘Stack required’:
The macros hig_Pdmapping.ijm and hig_Lineprofile.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 hig_Pdmapping.ijm and hig_Lineprofile.ijm do not support RGB images. Please prepare gray-scale images.
- 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).
- 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).
- Jones, D.K. et al. Spatial normalization and averaging of diffusion tensor MRI data sets. Neuroimage, 17, 592-617 (2002).
- 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).
- Higaki, T. et al. Statistical organelle dissection of Arabidopsis guard cells using image database LIPS. Sci. Rep. in press.
- Kohonen, T. Self-Organizing Maps. Springer-Verlag, ISBN 3-540-67921-9, New York, Berlin, Heidelberg (2001).
- ImageJ User Guide website. Available: "http://rsbweb.nih.gov/ij/docs/guide/userguide.html":http://rsbweb.nih.gov/ij/docs/guide/userguide.html. 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.
Associated Publications
Statistical organelle dissection of Arabidopsis guard cells using image database LIPS
by Takumi Higaki, Natsumaro Kutsuna, Yoichiroh Hosokawa, +5
Scientific Reports (21 May, 2012)