需要JoVE订阅才能观看此内容。 请登录或开始免费试用

方法文章

用于图像分析的时钟扫描方案:ImageJ 插件

17.5K 次观看

⸱

DOI:

10.3791/55819

⸱

2017年6月19日

 ,  ,  , 

本文内容

摘要

本文介绍了两种用于“时钟扫描”图像分析的新型 ImageJ 插件。这些插件扩展了原始 Visual Basic 6 程序的功能,更重要的是,通过将其集成到免费的图像分析软件包 ImageJ 中,使该程序能够被广大的科研群体所使用。

摘要

用于图像分析的时钟扫描协议是一种高效工具,可量化封闭或分割的凸形感兴趣区域内部、边界处以及外部(背景)的平均像素强度,从而生成平均积分径向像素强度分布曲线。该协议最初于2006年以Visual Basic 6脚本形式开发,但由于其形式所限,传播范围较小。为解决此问题,并响应其他研究者近年来的类似努力,我们将原始的时钟扫描协议代码转换为两个基于Java的插件,可兼容由NIH支持且免费提供的图像分析软件,如ImageJ或Fiji ImageJ。此外,这些插件还新增了多项功能,进一步扩展了原始协议的功能范围,例如可对多个感兴趣区域和图像堆栈进行分析。程序的后一特性在需要确定与时间和空间相关变化的应用中尤为有用。因此,对生物图像堆栈进行时钟扫描分析,可能适用于研究单个细胞内Na+或Ca++的扩散过程,也可用于分析突触连接或间隙连接细胞群体中的扩散性活动(例如,Ca++波)。本文中,我们描述了这些新的时钟扫描插件,并展示了其在图像分析中的一些应用实例。

引言

本研究的目的是提供一种无平台依赖且可自由获取的钟形扫描(Clock Scan)方案,供所有对此类图像分析感兴趣的科研人员使用。钟形扫描方案最初于2006年开发1,旨在改进在凸形感兴趣区域(ROI)内进行像素强度定量的现有方法,该方法具有更强的整合能力以及更高的空间分辨率。在采集过程中,该方案依次收集多条从ROI中心至其边界,或至ROI外预设距离处的径向像素强度分布,以测量位于感兴趣区域外的"背景"像素强度。该方案根据扫描方向上测得的细胞半径对这些分布进行缩放,因此每条径向扫描从中心到ROI边界的距离在X轴上始终为100%。最后,程序将这些单独的分布平均为一条综合的径向像素强度分布曲线。由于进行了缩放处理,由"Clock Scan"方案生成的平均像素强度分布曲线既不受ROI尺寸影响,也在合理范围内不受ROI形状的影响。该方法允许直接比较不同ROI的分布曲线,或在需要时进行平均或相减操作。此外,该方案还允许通过简单减去位于对象外部像素的平均强度,对任意对象的综合像素强度分布进行背景噪声校正。尽管目前仅在生物样本中进行了测试,但本方案为研究围绕某一起始点组织的物理或化学过程图像(例如物质从点源扩散)所使用的其他现有图像分析工具提供了有价值的补充1。

然而,原始图像分析方法的主要局限性在于该方案是基于 Visual Basic 6 (VB6) 编写的(代码),因此具有平台依赖性,且分发困难(需要 VB6 环境)。为解决这一问题,并响应其他研究者近期开展的类似工作2,我们将 VB6 版本的 Clock Scan 程序代码转换为两个基于 Java 的插件,可兼容由美国国立卫生研究院(NIH)支持的、免费开源且跨平台的图像分析软件 ImageJ3 和 Fiji ImageJ4。此外,这些插件还新增了多项功能,扩展了原始方案处理多个感兴趣区域(ROI)和图像堆栈的能力。许多图像分析应用在对多个对象进行统计分析时并不友好,因此通常仅展示代表性数据。借助多对象 Clock Scan ImageJ 插件,可实现对多个对象的同步分析。现在利用该插件扩展,即可对显微图像数据进行可靠的统计评估,特别是针对单个细胞/对象内信号强度分布的分析。本文将介绍 Clock Scan 插件,并展示其在图像分析中的应用实例。

访问受限。请登录或开始试用以查看此内容。

方案

1. 软件安装

  1. 根据各自网站上的建议,安装最新版本的捆绑 Java 以及 ImageJ 或 Fiji ImageJ(参见材料表中的相应网站链接)。在下文中,这两个程序均称为 "ImageJ"。
  2. 使用材料表中提供的链接复制 "Clock_Scan-1.0.1. jar" 和 "Multi_Clock_Scan-1.0.1.jar" 插件文件,并将其粘贴到 ImageJ 插件目录中。或者,在计算机硬盘上保存这些文件后,使用 "Plugins | Install plugin" 菜单选项安装这些文件。

2. 时钟扫描分析

  1. 标准时钟扫描插件 (图1):
    1. 使用 ImageJ "文件 | 打开" 菜单命令以打开感兴趣的图像。
    2. 单击“多边形”工具或“折线选择”工具,然后在图像上绘制,以勾勒出整个感兴趣区域(ROI)或该区域的一部分。参见 图1A 以多边形选择为例(内侧虚线轮廓)。
      注意:软件中提供的其他选择工具(矩形、椭圆和手绘线条选择)也可使用。
    3. 选择 "插件 | 时钟扫描" 从菜单中选择以打开标准时钟扫描协议的弹出选项窗口。请注意,此命令还将同时打开ROI管理器窗口,并自动添加轮廓。
    4. 使用插件选项窗口执行以下操作。
      1. 通过滚动条或修改相应输入框中的数值,审阅并调整 ROI 中心的 X 和 Y 坐标(系统自动计算为物理质心坐标)。参见 图1B.
      2. 根据扫描时需要覆盖物体外部背景区域的范围,通过使用以下方法调整扫描范围 "扫描范围" 滚动条。参见 图1A.
        注意:扫描限值(scan limit)是指在任意给定方向上,扫描超出物体边界的分数值;默认值为1.20,表示扫描长度将在扫描方向上比物体半径长20%;参见 图1A,外侧虚线)。
      3. 使用插件修改输出 "实际半径", "减去背景", "极坐标变换" 和/或 "标准差图" 复选框
      4. 点击 "好的" 以运行该插件。参见 图1C-H.
        注意:本方案输出结果的示例包括 "标准差图" 和 "极坐标变换" 或 "实际半径" 和 "极坐标变换" 已选选项显示在 图1C 和 1D 和 图1E 和 1F,分别。请注意,计算所得的标准差(SD)值表示对象各次径向像素强度扫描之间的变异情况。同时请注意 "ROI 选择长度" 插件窗口中的线条,显示以像素为单位测量的感兴趣区域(ROI)轮廓长度信息。
    5. 在生成的 "时钟扫描轮廓图",使用 "列表" 用于绘制灰度图像数据中两列(X 和 Y)数值的命令,以及绘制 RGB 图像数据中一列 X 和四列 Y 的数值的命令,其中 Y0、Y1、Y2 和 Y3 列将分别填入整体积分及各个颜色通道(红、绿、蓝)的像素强度值。
  2. 多区域兴趣时钟扫描插件 - 用于处理多个区域兴趣(ROI)图2):
    1. 打开包含多个感兴趣区域(ROI)的图像。
    2. 通过点击打开 ROI 管理器 "分析 | 工具 | 感兴趣区域管理器".
    3. 依次勾画(参见步骤 2.1.2)每个 ROI,并通过点击将其添加到 ROI 管理器中 "添加" 在 ROI Manager 窗口中执行此操作;对图像中的所有 ROI 均进行相同处理。使用 "分析|测量" 如果感兴趣的是感兴趣区域(ROI)的指标,请使用该命令。
      1. 参见 图2A 多段线选择示例 图2E 以查看多个多边形选择的示例。
    4. 选择 "多时钟扫描" 在 "插件" 打开协议选项弹出窗口的菜单。
    5. 使用协议选项窗口执行以下操作。
      1. 如有需要,按照步骤 2.1.4.2 重置扫描限值;默认值为 1.20。
      2. 如有需要,可通过勾选选项来绘制带有标准差(SD)误差棒的平均时钟扫描曲线 "标准差图" 盒。参见 图2C 和 D.
        注意:所计算的标准差(SD)值将表示不同对象的整体时钟扫描轮廓之间的变异。同时,注意插件窗口中显示信息的那条线 "选定的 ROI 数量".
      3. 点击 "好的" 执行该实验方案。
    6. 在生成的 "时钟扫描轮廓图",使用 "列表" 绘制所显示数值的命令 "绘图数值" 窗口。参见 "多时钟扫描轮廓图" 按颜色通道标注柱状图的图例窗口。
    7. 请注意,ROI 已编号,且针对任意给定颜色通道,其时钟扫描轮廓按照 ROI 被勾画并添加到图像中的顺序依次绘制。 "ROI 管理器".
  3. 多区域时钟扫描插件——处理图像堆栈(图3):
    1. 打开感兴趣的图像堆栈。
    2. 通过点击打开 ROI 管理器 "分析 | 工具 | 感兴趣区域管理器".
    3. 在图像堆栈中勾画出感兴趣区域(ROI),并按照步骤 2.1.2 和 2.2.3 中所述将其添加到 ROI 管理器中。使用 "分析 | 测量" 如果对 ROI 指标感兴趣,请使用该命令。
    4. 选择 "多时钟扫描" 在 "插件" 打开协议选项弹出窗口的菜单。
    5. 使用协议选项窗口执行以下操作。
      1. 按照步骤 2.1.4.2 中所述重置扫描范围;默认值为 1.20。
      2. 选择“绘制标准差图”选项,以绘制带有标准差(SD)误差棒的平均时钟扫描曲线" 盒子
        注意:所计算的标准差(SD)值将表示图像堆栈中选定对象不同实例之间的变异情况。同时,请注意插件窗口中显示相关信息的行 "图像堆栈中的图像数量".
      3. 点击 "好的" 执行该实验方案。
    6. 在 "时钟扫描轮廓图" 窗口,单击 "列表" 绘制所显示的数值 "绘图数值" 窗口中,Y 列编号表示图像在堆栈中的位置减 1。

访问受限。请登录或开始试用以查看此内容。

结果

本文用于说明目的的图像来自我们先前的细胞和组织生物学研究中创建的数据库5,6,7,以及来自艾伦小鼠脑图谱数据库8。这两个插件已在 ImageJ 1.50i/Java 1.8.0_77、ImageJ 2.0.0-rc-44/1.50e/Java 1.8.9_66 和 Fiji ImageJ 2.0.0-rc54/1.51g/Java 1.8.0_66 程序环境中成功测试。

图1展示了使用标准 Clock Scan 插件进行图像分析的代表性结果。对于这两个插件,时钟扫描程序的基本代码和主要步骤与原始方案中描述的基本相同1

访问受限。请登录或开始试用以查看此内容。

讨论

时钟扫描方案: 时钟扫描法是一种快速且简单的图像分析工具。与现有的常见图像分析方法(如线性像素强度扫描或感兴趣区域(ROI)的平均像素强度计算)相比,该方法的优势已在先前的出版物中得到详细描述。1,9简而言之,该方案通过量化位于感兴趣区域(ROI)中心不同距离处的像素强度(例如物体边界,或物体外预设位置(背景))来生成完整的径向像素强度分布曲线。由于可对每个ROI的“时钟扫描”(clock scan)曲线进行即时背景校正,因此在生物学应用中,该方法对样本内或样本间标记/染色的不均匀性,以及显微镜光源强度或荧光曝光时间的不稳定性具有较低的依赖性。时钟扫描曲线对物体大小和形状的独立性进一步拓展了该方案的应用范围,使其能够用于不同物体之间的比较,并可通过逐点减去背景曲线的方式进行校正。 "阳性" 和 "阴性" 对照物体

时钟扫描插件: 分发和共享原始方案的主要限制在于其代码...

访问受限。请登录或开始试用以查看此内容。

披露

作者声明不存在竞争性财务利益或其他利益冲突。

致谢

我们感谢德国柏林分子药理学莱布尼茨研究所的Tanja Maritzen博士和Fabian Feutlinske博士与我们分享他们版本的Fuji ImageJ Clock Scan插件,并激励我们开发本版本程序。我们还感谢马克斯·普朗克感染生物学研究所淋巴细胞发育系的Fritz Melchers博士慷慨允许我们使用其部门数据库中的图像,用于测试和改进该插件。资助:转化神经科学中心;美国国立卫生研究院基金:P30-GM110702-03。

访问受限。请登录或开始试用以查看此内容。

材料

本文使用的材料清单
姓名公司目录编号评论
计算机任意与下文所列软件兼容
ImageJ 或 Fiji ImageJNIHhttps://imagej.nih.gov/ij/ 或 https://fiji.sc/需捆绑 Java 1.8 或更高版本
Clock-scan 插件免费软件https://sourceforge.net/projects/clockscan/Clock_Scan-1.0.1 jar 和 Multi_Clock_Scan-1.0.1/jar
Origin 9.0OriginLabNorthampton, MA, USA本程序用于生成原始 Clock Scan 数据的部分图表。其他任意图形软件也可用于实现此功能

参考文献

  1. Dobretsov, M., Romanovsky, D. "Clock-scan" protocol for image analysis. Am J Physiol Cell Physiol. 291, 869-879 (2006).
  2. Feutlinske, F., Browarski, M., Ku, M. C., et al. Stonin1 mediates endocytosis of the proteoglycan NG2 and regulates focal adhesion dynamics and cell motility. Nat Commun. 6, 8535(2015).
  3. Schneider, C. A., Rasband, W. S., Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 9, 671-675 (2012).
  4. Schindelin, J., Arganda-Carreras, I., Frise, E., et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 9, 676-682 (2012).
  5. Dobretsov, M., Hastings, S. L., Stimers, J. R. Non-uniform expression of alpha subunit isoforms of the Na+/K+ pump in rat dorsal root ganglia neurons. Brain Res. 821, 212-217 (1999).
  6. Hayar, A., Gu, C., Al-Chaer, E. D. An improved method for patch clamp recording and calcium imaging of neurons in the intact dorsal root ganglion in rats. J Neurosci Methods. 173, 74-82 (2008).
  7. Dobretsov, M., Pierce, D., Light, K. E., Kockara, N. T., Kozhemyakin, M., Wight, P. A. Transgenic mouse model to selectively identify alpha3 Na,K-ATPase expressing cells in the nervous system. Society for Neuroscience. , Online Program No. 123.01/B54 1(2015).
  8. Lein, E. S., Hawrylycz, M. J., Ao, N., et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature. 445, 168-176 (2007).
  9. Romanovsky, D., Mrak, R. E., Dobretsov, M. Age-dependent decline in density of human nerve and spinal ganglia neurons expressing the alpha3 isoform of Na/K-ATPase. Neuroscience. 310, 342-353 (2015).
  10. Campbell, J., Singh, D., Hollett, G., et al. Spatially selective photoconductive stimulation of live neurons. Front Cell Neurosci. 8, 142(2014).
  11. Yuryev, M., Pellegrino, C., Jokinen, V., et al. In vivo Calcium Imaging of Evoked Calcium Waves in the Embryonic Cortex. Front Cell Neurosci. 9, 500(2015).
  12. Qiao, M., Sanes, J. R. Genetic Method for Labeling Electrically Coupled Cells: Application to Retina. Front Mol Neurosci. 8, 81(2015).

访问受限。请登录或开始试用以查看此内容。

重印与许可

探索更多文章

像素强度分析感兴趣区域图像堆栈分析径向像素强度多时钟扫描背景扣除标准差图感兴趣区域管理器