Name: AMEBaS:偏振单细胞比例荧光延时的自动中线提取和背景减法 (视频) | JoVE Journal
Uploaded: 2023-06-23T12:30:00.000+00:00
Duration: 6 min 3 s
Description: 学习使用AMEBaS自动提取极化单细胞中线并分析荧光延时数据。该方法通过自动化消除主观偏差，实现高通量细胞动态定量，适用于花粉管、根毛等多种细胞类型的时空动态研究。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

作者感谢圣保罗研究基金会（FAPESP）提供的资助，项目编号为 2015/22308-2、2019/23343-7、2019/26129-6、2020/06744-5、2021/05363-0，以及巴西国家科学技术发展委员会（CNPq）、美国国立卫生研究院（NIH）R01 基金 GM131043 和美国国家科学基金会（NSF）基金 MCB1714993、MCB1930165 提供的经费支持。根毛数据在 Andrea Bassi 教授和 Alex Costa 教授的指导和基础设施支持下完成。

1. 交互式笔记本实验方案
可直接在网页上使用 Google Colab 运行 Jupyter 笔记本，地址为 https://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb，下文所述操作说明即基于此平台。或者，也可从 https://github.com/badain/amebas 获取该 Jupyter 笔记本，下载后在本地 Jupyter 环境中运行（Anaconda 可提供简便且跨平台的安装方式）。完整的测试数据可在 Zenodo 上获取（https://doi.org/10.5281/zenodo.7975350），其中包含表达 pH 或 Ca2+ 报告分子的拟南芥花粉管的单通道和双通道数据14。分析流程已被划分为多个部分，每个步骤均可在设置用户特定选项后，通过点击播放按钮执行。本研究所需文件包含在 AMEBaS- main 压缩文件夹中（补充代码文件 1）。
<ol>
	<li>打开 Jupyter 笔记本并读取延时成像文件。
	<ol>
		<li>访问上述 Google Colab 中交互式笔记本的主页，或下载并打开该笔记本 AMEBaS_Local.ipynb 从 GitHub 获取笔记本。</li>
		<li>为输入和输出数据准备目录结构：
		<ol>
			<li>如果使用本地版本，将荧光延时图像以 TIFF 文件或 DV 文件格式放入名为 数据 必须位于程序的根目录中。一个名为 出 必须创建以接收生成的数据。然后运行 设置 代码块。</li>
			<li>如果在 Google Colab 上使用笔记本，请运行以下代码 设置 自动生成代码块 数据 和 出 文件夹。</li>
		</ol>
		</li>
		<li>运行 文件输入 单击以读取延时数据的代码块 播放 按钮。如果使用的是笔记本的 Google Colab 版本，请点击 选择文件 按钮，以直接上传延时文件至 数据 文件夹 
		注意：通道数量将根据图像的维度自动检测。</li>
		<li>通过将“verbose”参数设置为，选择是否生成每个步骤的额外输出 是 或 错误.</li>
	</ol>
	</li>
	<li>从背景中检测出主要细胞并进行分割（图2).
	<ol>
		<li>运行 单细胞分割 通过点击自动将目标细胞与背景分离的代码块 播放 按钮 
		注意：在通过IsoData阈值法将前景与背景分离并选取最大区域以去除不需要的伪影之前，将中值滤波和高斯滤波作为预处理步骤，用于消除不必要的噪声。
		<ol>
			<li>通过调整“sigma”变量中高斯函数使用的sigma值，可微调分割掩膜的平滑程度。默认值为2.0。</li>
			<li>设置变量 估计 至 错误 将从 Isodata 估计的阈值直接存储或 是 使用局部多项式回归（LOESS）在相邻帧之间进行平滑处理。通过调整其函数参数以优化性能 n_points 变量。默认值为 40。</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>沿细胞突起描记中线图3).
	<ol>
		<li>运行 细胞 中线 示踪 单击代码块 播放 使用 Lee 法自动骨架化细胞的按钮18 并通过线性外推法延长最后一个骨架的末端。
		<ol>
			<li>选择仅在最后一帧或每帧中追踪中线，通过调整 完全骨骼化 论点 
			注意：当所有图像帧均已完成骨架化时，将跳过外推步骤。</li>
			<li>通过调整外推过程中需插值的骨架点比例 插值分数 变量。默认值为 0.25。</li>
			<li>通过更改变量来选择中线外推的长度 外推长度默认值为 -1，表示将骨架延伸至最近的边缘。</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>为每个通道生成明场图像 图4).
	<ol>
		<li>运行第一次 数据可视化 单击代码块 播放 单击按钮可自动为两个通道生成明场显微图像。
		<ol>
			<li>通过调整变量来选择用于平滑处理的高斯核大小 kymograph_kernel. 
			注意：该参数对应于用于平均像素强度的邻域大小（以像素为单位）。默认值为 3 像素 × 3 像素。</li>
			<li>非延展骨架生成的截断Kymograph必须使用自定义色阶来正确显示其强度。选择指定给背景色（黑色）的强度百分比分数，并进行调整 偏移分数 变量。默认值为 0.7。</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>计算通道间的比值（图5).
	<ol>
		<li>运行第二次 数据可视化 点击代码块 播放 按钮，可自动生成比率式显影图（ratiometric kymograph）和比率式时间序列图像（ratiometric timelapse）图6). 
		注意：此步骤仅在使用双通道时间序列时可用。步骤 1.2.1.2 中存储的背景强度阈值将从每个通道中减去。
		<ol>
			<li>调节 转换比率 用于切换在比值计算过程中作为分子和分母的通道顺序的变量。默认值为 错误.</li>
			<li>选择是否通过调整中值滤波器的参数对比率时间序列进行进一步平滑处理 平滑比 变量。默认值为 错误.</li>
			<li>选择是否通过调整分母通道低信号产生的异常值来去除 outliers 剔除离群值 变量。默认为 是 将离群值定义为高于第三四分位数（即75%的数值所在位置）1.5倍四分位距的数值。</li>
			<li>选择是否通过调整变量来导出比例输出中的背景 背景比值。默认值为 错误，用零替换它。</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
<img alt="显示单细胞数据分析中滤波和信号分离步骤的信号处理示意图。" class="xfigimg" src="/files/ftp_upload/64857/64857fig02.jpg" data-altupdatedat="1788304762000" data-alt="Figure 2"> 
图 2：单细胞分割步骤。 使用滤波、阈值分割和区域标记等图像处理技术来分离目标信号（步骤 1.2）。该组数据的最低强度值为 2556，中位值为 3441，最高强度值为 32125。 <a href="https://www.jove.com/files/ftp_upload/64857/64857fig02large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="流体动力学模拟，流线型物体周围气流可视化的示意图。" class="xfigimg" src="/files/ftp_upload/64857/64857fig03.jpg" data-altupdatedat="1788304762000" data-alt="图3"> 
图3：中轴线追踪概述——单个细胞的中轴线通过计算其骨架（白色）获得。细胞尖端（品红色）由骨架末端的最后几个点线性外推得出（步骤1.3）。在此图中，中轴线及其尖端均叠加在原始细胞图像之上。<a href="https://www.jove.com/files/ftp_upload/64857/64857fig03large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="通道随时间变化的像素强度波动图，用于比较，表明信号分析。" class="xfigimg" src="/files/ftp_upload/64857/64857fig04.jpg" data-altupdatedat="1788304762000" data-alt="图 4"> 
图 4：时间序列的波动图—在“complete_skeletonization”功能关闭时生成的各通道波动图比较（步骤 1.4）。纵轴表示时间进程，横轴表示单个细胞沿外推中线路径的平均强度。对于本组数据，颜色映射表示以下数值：通道 1 最低强度：2886，中位数：3167，最高强度：21021；通道 2 最低强度：3030，中位数：3400，最高强度：29688。<a href="https://www.jove.com/files/ftp_upload/64857/64857fig04large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<img alt="局部回归LOESS图；像素强度与帧数；数据拟合分析；多项式趋势。" class="xfigimg" src="/files/ftp_upload/64857/64857fig05.jpg" data-altupdatedat="1788304762000" data-alt="图5"> 
图5：背景阈值平滑处理— 背景分割阈值通过等量算法（isodata算法）估算，随后通过局部多项式回归进行平滑处理（步骤1.5）。<a href="https://www.jove.com/files/ftp_upload/64857/64857fig05large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="细胞时间序列的比率分析；分割成像、明场图、信号强度曲线。" class="xfigimg" src="/files/ftp_upload/64857/64857fig06.jpg" data-altupdatedat="1788304762000" data-alt="Figure 6"> 
图6：比率分析结果—（A）比率时间序列的最后一帧与原始第一通道分割图像的比较。（B）由比率时间序列生成的明场图（步骤1.5）。<a href="https://www.jove.com/files/ftp_upload/64857/64857fig06large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
2. 批量模式操作流程
<ol>
	<li>从 AMEBaS GitHub 仓库下载文件 pipeline.py，并将其放置在与数据相同的目录中。</li>
	<li>在命令行中程序文件后输入文件路径。</li>
	<li>如需显示流程的内部步骤，可包含位置参数 --v。</li>
	<li>包含 --s 以选择高斯滤波预处理步骤中用于细胞分割的 sigma 值。默认值为 2。</li>
	<li>包含 --a 以对时间序列的每一帧追踪中线。默认情况下，流程仅使用最后一帧。</li>
	<li>包含 --f 以选择在插值过程中使用的骨架比例 [0,1]。默认值为 0.25。</li>
	<li>包含 -e 以选择外推骨架的长度（像素）。默认值为 -1，表示将骨架延伸至最近的边缘。</li>
	<li>包含 --sf 以选择在非外推的显迹图中被移至背景的颜色范围比例。默认值为 0.7。</li>
	<li>包含 --k 以确定显迹图高斯滤波中所用核的大小。默认值为 3。</li>
	<li>包含 --eb 以通过 LOESS 多项式回归对帧特异性背景阈值强度进行全局背景阈值强度估计。
	<ol>
		<li>通过修改参数 --n 自定义用于背景阈值平滑处理的 LOESS 点数。默认值为 40。</li>
	</ol>
	</li>
	<li>如果时间序列包含两个通道，可通过包含 -r 或 --switch_ratio 来切换比值计算中分子和分母所使用的通道。默认情况下，第二个通道为分子，第一个通道为分母。</li>
	<li>通过 --sm 参数选择是否在比值时间序列生成过程中进一步使用中值滤波进行平滑处理。默认值为 False。</li>
	<li>包含 -o 以在生成比值时间序列时排除具有异常强度的像素。</li>
	<li>使用参数 --b 选择是否导出比值输出中的背景。默认值为 False，表示将其替换为零。</li>
	<li>按 enter 运行。输出文件将生成在程序文件所在的同一目录中。</li>
</ol>

<ol><li>Drubin, D. G., Nelson, W. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Origins+of+cell+polarity%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Origins of cell polarity.</a> Cell. 84 (3), 335-344 (1996).</li>
<li>Wodarz, A., Näthke, I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cell+polarity+in+development+and+cancer%5BTitle%5D)+AND+%22Nature+Cell+Biology%22%5BJournal%5D)">Cell polarity in development and cancer.</a> Nature Cell Biology. 9 (9), 1016-1024 (2007).</li>
<li>Palanivelu, R., Preuss, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pollen+tube+targeting+and+axon+guidance%3A+parallels+in+tip+growth+mechanisms%5BTitle%5D)+AND+%22Trends+in+Cell+Biology%22%5BJournal%5D)">Pollen tube targeting and axon guidance: parallels in tip growth mechanisms.</a> Trends in Cell Biology. 10 (12), 517-524 (2000).</li>
<li>Portes, M. T., et al. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPortes%2C+MT+%22The+Pollen+Tube+Oscillator%3A+Integrating+Biophysics+and+Biochemistry+into+Cellular+Growth+and+Morphogenesis.+Rhythms+in+Plants%3A+Dynamic+Responses+in+a+Dynamic+Environment%22">The Pollen Tube Oscillator: Integrating Biophysics and Biochemistry into Cellular Growth and Morphogenesis. Rhythms in Plants: Dynamic Responses in a Dynamic Environment</a>. , Springer. Cham. (2015).</li>
<li>Wudick, M. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((CORNICHON+sorting+and+regulation+of+GLR+channels+underlie+pollen+tube+Ca2%2B+homeostasis%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">CORNICHON sorting and regulation of GLR channels underlie pollen tube Ca2+ homeostasis.</a> Science. 360 (6388), 533-536 (2018).</li>
<li>Hoffmann, R. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Plasma+membrane+H%2B-ATPases+sustain+pollen+tube+growth+and+fertilization%5BTitle%5D)+AND+%22Nature+Communications%22%5BJournal%5D)">Plasma membrane H+-ATPases sustain pollen tube growth and fertilization.</a> Nature Communications. 11 (1), 1-15 (2020).</li>
<li>Michard, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Glutamate+receptor-like+genes+form+Ca2%2B+channels+in+pollen+tubes+and+are+regulated+by+pistil+D-serine%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Glutamate receptor-like genes form Ca2+ channels in pollen tubes and are regulated by pistil D-serine.</a> Science. 332 (6028), 434-437 (2011).</li>
<li>Damineli, D. S., Portes, M. T., Feijó, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Oscillatory+signatures+underlie+growth+regimes+in+Arabidopsis+pollen+tubes%3A+computational+methods+to+estimate+tip+location%2C+periodicity%2C+and+synchronization+in+growing+cells%5BTitle%5D)+AND+%22Journal+of+Experimental+Botany%22%5BJournal%5D)">Oscillatory signatures underlie growth regimes in Arabidopsis pollen tubes: computational methods to estimate tip location, periodicity, and synchronization in growing cells.</a> Journal of Experimental Botany. 68 (12), 3267-3281 (2017).</li>
<li>Li, K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+optimized+genetically+encoded+dual+reporter+for+simultaneous+ratio+imaging+of+Ca2%2B+and+H%2B+reveals+new+insights+into+ion+signaling+in+plants%5BTitle%5D)+AND+%22New+Phytologist%22%5BJournal%5D)">An optimized genetically encoded dual reporter for simultaneous ratio imaging of Ca2+ and H+ reveals new insights into ion signaling in plants.</a> New Phytologist. 230 (6), 2292-2310 (2021).</li>
<li>Sadoine, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Designs%2C+applications%2C+and+limitations+of+genetically+encoded+fluorescent+sensors+to+explore+plant+biology%5BTitle%5D)+AND+%22Plant+Physiology%22%5BJournal%5D)">Designs, applications, and limitations of genetically encoded fluorescent sensors to explore plant biology.</a> Plant Physiology. 187 (2), 485-503 (2021).</li>
<li>Grenzi, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Illuminating+the+hidden+world+of+calcium+ions+in+plants+with+a+universe+of+indicators%5BTitle%5D)+AND+%22Plant+Physiology%22%5BJournal%5D)">Illuminating the hidden world of calcium ions in plants with a universe of indicators.</a> Plant Physiology. 187 (2), 550-571 (2021).</li>
<li>Munglani, G., Vogler, H., Grossniklaus, U. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fast+and+flexible+processing+of+large+FRET+image+stacks+using+the+FRET-IBRA+toolkit%5BTitle%5D)+AND+%22PLoS+Computational+Biology%22%5BJournal%5D)">Fast and flexible processing of large FRET image stacks using the FRET-IBRA toolkit.</a> PLoS Computational Biology. 18 (4), 1009242(2022).</li>
<li>Ridler, T., Calvard, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Picture+thresholding+using+an+iterative+selection+method%5BTitle%5D)+AND+%22IEEE+Transactions+on+Systems%2C+Man%2C+and+Cybernetics%22%5BJournal%5D)">Picture thresholding using an iterative selection method.</a> IEEE Transactions on Systems, Man, and Cybernetics. 8 (8), 630-632 (1978).</li>
<li>Portes, M. T., Feijó, J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPortes%2C+MT+%22Growing+Arabidopsis+pollen+tubes+expressing+genetically+encoded+reporters+for+calcium+and+pH%22">Growing Arabidopsis pollen tubes expressing genetically encoded reporters for calcium and pH</a>. , (2023).</li>
<li>Candeo, A., Doccula, F. G., Valentini, G., Bassi, A., Costa, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Light+sheet+fluorescence+microscopy+quantifies+calcium+oscillations+in+root+hairs+of+Arabidopsis+thaliana%5BTitle%5D)+AND+%22Plant+%26+Cell+Physiology%22%5BJournal%5D)">Light sheet fluorescence microscopy quantifies calcium oscillations in root hairs of Arabidopsis thaliana.</a> Plant & Cell Physiology. 58 (7), 1161-1172 (2017).</li>
<li>Romano Armada, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+vivo+light+sheet+fluorescence+microscopy+of+calcium+oscillations+in+Arabidopsis+thaliana%5BTitle%5D)+AND+%22Methods+in+Molecular+Biology%22%5BJournal%5D)">In vivo light sheet fluorescence microscopy of calcium oscillations in Arabidopsis thaliana.</a> Methods in Molecular Biology. 1925, 87-101 (2019).</li>
<li>Krebs, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((FRET-based+genetically+encoded+sensors+allow+high-resolution+live+cell+imaging+of+Ca2%2B+dynamics%5BTitle%5D)+AND+%22The+Plant+Journal%22%5BJournal%5D)">FRET-based genetically encoded sensors allow high-resolution live cell imaging of Ca2+ dynamics.</a> The Plant Journal. 69 (1), 181-192 (2012).</li>
<li>Lee, T. -C., Kashyap, R. L., Chu, C. -N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Building+skeleton+models+via+3-D+medial+surface+axis+thinning+algorithms%5BTitle%5D)+AND+%22CVGIP%3A+Graphical+Models+and+Image+Processing%22%5BJournal%5D)">Building skeleton models via 3-D medial surface axis thinning algorithms.</a> CVGIP: Graphical Models and Image Processing. 56 (6), 462-478 (1994).</li>
<li>Nunez-Iglesias, J., Blanch, A. J., Looker, O., Dixon, M. W., Tilley, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+new+Python+library+to+analyse+skeleton+images+confirms+malaria+parasite+remodelling+of+the+red+blood+cell+membrane+skeleton%5BTitle%5D)+AND+%22PeerJ%22%5BJournal%5D)">A new Python library to analyse skeleton images confirms malaria parasite remodelling of the red blood cell membrane skeleton.</a> PeerJ. 6, 4312(2018).</li>
</ol>

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

本文介绍的新方法是一种强大的工具，可用于简化和自动化极化细胞荧光显微镜图像堆栈的分析。目前文献中描述的方法（如ImageJ的Kymograph插件）需要手动追踪目标极化细胞的中线，这一过程不仅耗时，而且容易引入人为误差。本流程中，中线的定义由一种数值方法支持18,19，该方法执行骨架化处理，从而消除了主观判断，为该过程引入了定量标准。这对于拥有大量数据的研究人员尤为有用，且该流程可根据不同需求进行定制，包括为每一帧提取中线。此外，背景值的扣除通常是任意的，往往由人工为某一通道的所有帧选择单一背景值。本文采用等值区域（isodata）分割法，客观地确定每一帧的背景阈值，并可选地通过局部多项式回归对其进行平滑处理，以捕捉由荧光衰减（光漂白）引起的长期变化。尽管在背景选择中通过像素面积选取较大对象的方式忽略了可能的伪影和次要成分，但也可采用其他方法（如FRET-IBRA12）来消除阴影等效应。空间伪影（如阴影）可能影响通过阈值分割得到的细胞形态，这可能是AMEBaS标志所示比率成...

细胞极性是一种基本的生物学过程，其中一系列空间聚集的分子和结构协同作用，最终形成特化的形态学子细胞区域1。细胞分裂、生长和迁移均依赖于此类极性位点，而极性的丧失则与上皮组织相关疾病中的癌症发生有关2。
顶端生长的细胞是极性现象的一个显著例子，其细胞顶端的极性位点通常会对外部胞外信号重新定向3。这类细胞包括正在发育的神经突、真菌菌丝、根毛和花粉管，其中多种细胞过程从细胞顶端向主干方向表现出明显差异。特别是在花粉管中，肌动蛋白聚合、囊泡运输以及离子浓度均呈现显著的极性分布，形成以顶端为中心的梯度分布模式4。花粉管是开花植物的雄配子体，负责通过细胞顶端进行快速生长，将精子细胞递送至胚珠，其生长速率在单细胞中属于已知最快之列。钙离子5（Ca2+）和质子6（H+）等离子的顶端聚焦梯度在维持花粉管生长中起着关键作用，这对于实现其主要生物学功能——最终完成双受精过程——至关重要5,6。因此，定量分析顶端生长细胞中轴线上时空动态变化的方法，对于研究极性生长背后的细胞与分子机制至关重要7,8,9。研究人员常使用时空图（kymograph），即一种表示细胞中轴像素强度（如列为位置）随时间变化（如行为时间）的矩阵，从而通过斜线模式可视化细胞的生长与迁移过程(图1)。尽管该方法具有实用性，但时空图通常依赖手动追踪中轴线获得，容易引入主观偏差和人为误差，且操作繁琐耗时。因此，亟需一种自动化的中轴提取方法。本文介绍的 AMEBaS 流程正是以此为首要功能：Automatic Midline Extraction and Background Subtraction of ratiometric fluorescence time lapses of polarized single cells（极化单细胞比率荧光延时成像的自动中轴提取与背景扣除）。
在实验方法方面，利用基因编码的荧光探针可实现对单个细胞中目标离子/分子/物种的定量成像10。在日益增多的探针选择中，比率型探针因其在结合与未结合目标分子时可发射不同荧光波长而成为最准确的方法之一11。该特性允许通过计算两个通道的荧光强度比值（每个通道均扣除其特异性背景信号）来校正探针在细胞内分布的空间异质性。然而，由于阴影效应（图像边缘区域相对于中心区域存在亮度差异）导致的空间变化以及荧光团随时间衰减（光漂白）引起的时间变化，为每个通道和时间点估计背景阈值可能是一项复杂任务12。尽管存在多种可能的解决方法，本文提出采用Isodata算法获得的分割阈值自动确定背景强度13，并进一步通过多项式回归对帧间阈值进行平滑处理作为标准方法。此前研究中已去除与目标细胞无关的荧光异质性所导致的空间成分12，但本方法未考虑此类因素。自动阈值分割可通过多种方法实现，但经验表明Isodata算法效果最佳。因此，自动背景值扣除与比率计算构成了AMEBaS的第二个主要功能（图1），该功能整体接收双通道荧光显微图像序列作为输入，估算细胞中线位置及各通道特异性背景，经背景扣除、平滑处理和异常值剔除后，输出两个通道及其比值的明场图（主要输出#1），同时生成一组比率图像序列（主要输出#2）。
AMEBaS 在显微镜下对生长中的拟南芥花粉管进行荧光时间序列测试，花粉管中分别表达了花粉特异性启动子 LAT52 驱动的 Ca2+（CaMeleon）8 或 pH（pHluorin）6 比值型传感器。使用倒置显微镜，每 4 秒采集一次各通道图像，配备前照式相机（2560 像素 × 2160 像素，像素大小 6.45 μm）、荧光照明装置以及 63 倍水浸物镜（1.2 NA）。CaMeleon 使用的滤光片设置为：激发光 426–450 nm（CFP）和 505–515 nm（YFP），发射光 458–487 nm（CFP）和 520–550 nm（YFP）；pHluorin 的滤光片设置为：激发光 318–390 nm（DAPI）和 428–475 nm（FITC），发射光 435–448 nm（DAPI）和 523–536 nm（FITC）。完整数据集已上传至 Zenodo 供测试使用（DOI: 10.5281/zenodo.7975350）14。
此外，该流程还使用根毛数据进行了测试，成像采用光片显微镜（SPIM）进行，实验方法如前所述15,16，其中拟南芥根毛在UBQ10启动子调控下表达遗传编码的Ca2+探针NES-YC3.6。17 自制的LabView软件控制相机采集、样品平移以及光片显微镜的快门，可同时观察cpVenus和CFP两个荧光通道，并实现实时比值成像。每个延时比值图像均为对样品每隔3 μm采集的15个层面图像中cpVenus与CFP荧光通道图像进行最大强度投影（MIP）所得。延时采集的MIP比值图像序列被保存后直接用于AMEBaS分析。
尽管该流程可用于多种生长和迁移的细胞类型，但它专门设计用于分析仅在顶端生长的细胞，例如花粉管、根毛和真菌菌丝，这些细胞在连续图像帧之间具有非生长胞质区域的对应关系。当这种对应关系不存在时，用户应在步骤 1.3.1.1 中选择 complete_skeletonization 选项（详见讨论部分）。
<img alt="展示过滤、骨架化、中线追踪和显微图分析的细胞分割示意图" class="xfigimg" src="/files/ftp_upload/64857/64857fig01.jpg" data-altupdatedat="1788304762000" data-alt="图 1"> 
图 1：流程工作概述。 AMEBaS 流程对显微镜下的延时图像进行分析和处理，主要包括三个步骤：单细胞分割、中线追踪和显微图生成。<a href="https://www.jove.com/files/ftp_upload/64857/64857fig01large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>

<table><tbody><tr><td>Github</td><td>Github</td><td></td><td>https://github.com/badain/amebas</td></tr><tr><td>Google Colab</td><td>谷歌</td><td></td><td>https://colab.research.google.com/github/badain/amebas/blob/main/AMEBAS_Colab.ipynb</td></tr></tbody></table>

amebas：极化单细胞比率荧光延时成像的自动中线提取与背景扣除

AMEBaS: Automatic Midline Extraction and Background Subtraction of Ratiometric Fluorescence Time-Lapses of Polarized Single Cells

细胞极性是一种由空间上聚集的分子和结构共同建立的宏观现象，最终在亚细胞水平上形成特化的功能区域。该现象与产生不对称形态结构相关，这些结构是细胞分裂、生长和迁移等关键生物学功能的基础。此外，细胞极性的破坏已被证实与多种组织相关疾病有关，例如癌症和胃部发育异常。
目前评估单个极化细胞中荧光报告分子时空动态的方法通常需要手动沿细胞长轴描绘中线，这一过程耗时且容易引入显著偏差。此外，尽管比率分析可通过两个荧光通道校正报告分子分布不均的问题，但背景扣除技术往往具有任意性，且缺乏统计学支持。
本文介绍了一种新颖的计算流程，用于通过细胞极性模型——花粉管/根毛生长及胞质离子动态——自动分析并量化单个细胞的时空行为。该流程开发了一种三步算法，用于处理比率成像图像，并提取细胞内动态与生长的定量表征。第一步通过像素强度空间中的阈值技术，将细胞与背景分离，生成二值化掩膜。第二步通过骨架化操作，沿细胞中线追踪一条路径。最后，第三步将处理后的数据以比率时间序列形式呈现，并生成比率显影图（即随时间变化的一维空间剖面图）。本方法利用来自生长中花粉管的基因编码荧光报告分子所获取的比率成像数据进行了验证。该流程可更快速、更少偏差且更准确地表征极化细胞中线上的时空动态，从而推进了研究细胞极性可用的定量工具集。AMEBaS Python 源代码位于：&#160;https://github.com/badain/amebas.git

AMEBaS 流程可自动从荧光显微镜图像序列中提取极化单细胞的中线动态信息，从而减少耗时并降低人为误差。该方法通过在生长的单细胞中生成明暗图（kymographs）和比率图像序列（图1）来量化这些时间变化过程。该流程可调整用于迁移中的单细胞分析，但尚需进一步实验验证。AMEBaS 以 Python 编写，提供交互式 Jupyter Notebook 实现形式（详见交互式笔记本协议部分），便于无编程经验的用户使用；同时也提供命令行工具形式（见批量模式协议部分），可用于对多个图像序列使用相同参数集进行分析。尽管可使用一个或两个荧光通道，但建议采用具有两个发射通道的比率探针，因为探针未结合状态的荧光发射可减轻细胞质中蛋白质分布不均所引起的空间异质性。
该流程首先在最强通道上生成最大细胞的二值化掩膜，并导出为名为 Filename_binary_mask.tiff (图2...

目前用于分析极化单细胞内动态的方法通常依赖人工操作，且缺乏标准化。本文介绍了一种新颖的图像分析流程，可在用户友好的在线界面中自动提取单个极化细胞的中线，并对延时成像中的时空行为进行量化。

AMEBaS：极化单细胞比率荧光延时成像的自动中线提取与背景扣除

细胞极性是由一系列空间上富集的分子和结构共同作用所建立的一种宏观现象，最终导致……的出现

该方法通过表达特定目标（如肌动蛋白、离子或囊泡动态）的荧光报告分子或染料，定量分析极化细胞沿其中轴线的时空动态。自动化版本消除了人为偏差，并提高了可扩展性，尤其适用于那些原本需手动进行数字化操作、耗时且易受主观影响的任务。首先，打开 Jupyter Notebook 并读取延时成像文件，可通过访问 Google Colab 交互式笔记本主页，或从 GitHub 下载并打开 AMEBaS_local IP YNB 文件来实现。使用本地版本准备输入和输出数据的目录结构时，将荧光延时成像文件（TIF 或 DV 格式）放入程序根目录下名为 data 的文件夹中。创建一个 out 文件夹用于接收生成的数据，然后运行设置代码块。如果在 Google Colab 中使用笔记本，请执行设置代码块以自动生成 data 和 out 文件夹，然后运行文件输入代码块，通过点击播放按钮读取延时成像数据。如果使用 Google CoLab，请通过点击“选择文件”按钮将延时成像文件上传至数据文件夹。接下来，通过将 verbose 参数设置为 true 或 false，选择是否为每一步生成额外输出。为了检测主要细胞并将其与背景分离，请点击播放按钮运行单细胞分割代码块，以自动将目标细胞从背景中分离出来。调整 sigma 变量中的 sigma 值，以精细调节分割掩膜的平滑程度。默认值为 2.0。现在将变量 estimate 设置为 false，以直接从 ISO 数据中存储估计的阈值；或设置为 true，以通过局部多项式回归在相邻帧之间对阈值进行平滑处理。调整 n_points 变量以微调函数，默认值为 40。要沿细胞突起追踪中线，请通过点击播放按钮运行细胞中线追踪代码块，使用 Lee 方法对细胞进行骨架化，并通过线性外推法延伸最后一个骨架的末端。接下来，通过调整 complete_skeletonization 参数选择中线追踪选项，以在最后一帧上进行追踪，或每帧追踪一次。调整 interpolation_fraction 变量以设置在延伸过程中用于插值的骨架点比例，默认值为 0.25。接下来，修改 extrapolation_length 变量以确定中线延伸的长度。默认值为负一，这会将骨架延伸至最近的边缘。现在通过点击播放按钮运行第一个数据可视化代码块，以自动生成两个通道的显微轨迹图。通过调整 kymograph_kernel 变量来选择用于平滑处理的高斯核大小。为了正确显示非延伸骨架的强度，需要为其加帽的显影图（kymograph）设置自定义的颜色映射。调整 shift_fraction 变量，以选择将被指定为背景色（黑色）的强度所占的百分比。然后点击播放按钮运行第二个数据可视化代码块，自动生成比率度量显影图和时间序列图像。调整 switch_ratio 变量以更改在比值计算过程中作为分子和分母的通道顺序，默认值为 false。通过调整 smooth_ratio 变量，检查比值时间序列是否需要通过中值滤波进行额外平滑处理，默认情况下该值为 false。接下来，通过调整 reject_outliers 变量选择是否移除或保留由于分母通道信号过低而产生的离群值。离群值定义为高于第三四分位数1.5倍四分位距的数值。最后，通过调整 background_ratio 变量选择是否导出背景值和比值指标结果。默认情况下为 false，表示将背景替换为零。AMEBaS 在多个数据集上进行了测试，例如表达钙离子探针 cameleon 的花粉管、表达 pH 指示剂 florin 的花粉管，以及表达钙离子探针 NESYC 3.6 的根毛，尽管这些样本在生长方向、成像技术、荧光报告分子和细胞类型方面存在差异，但仍成功实现了分析。请记得执行该代码块，以安装所需的软件包。此外，必须针对您的图像专门调整sigma值，以便获得更准确的结果。生成的显微轨迹图可用于分析目标信号的时空动态。此外，您还可以通过计算其寻找方法、生长过程和时间序列来分析正确的解决方案。

amebas-automatic-midline-extraction-background-subtraction

Research

科研

JoVE Journal

Developmental Biology

发育生物学

AMEBaS:偏振单细胞比例荧光延时的自动中线提取和背景减法

通过追踪光转化的基质靶向GFP在随时间变化的线粒体网络中的平衡情况，来测量线粒体融合。迄今为止，每次仅能对一个细胞进行长达一小时的动力学分析。我们提出了一种可同时测量多个细胞的方法，从而加快了数据采集过程。...

一种基于mtPA-GFP的快速、高分辨率线粒体融合检测方法：利用共聚焦显微镜并行获取多个细胞的动力学数据

我们描述了一种超分辨率成像方法，用于研究细菌FtsZ环的结构组织，该结构是细胞分裂所必需的重要装置。该方法基于对光激活定位显微镜（PALM）图像的定量分析，可应用于其他细菌细胞骨架蛋白的研究。...

细菌分裂机器的超分辨率成像

质膜蛋白和脂质的空间分布与时间动态是生物学中的热点问题。本文通过时空图像涨落分析来探讨这一问题，该方法在概念上可提供与单粒子追踪相同的物理量，但使用小分子标记物和常规显微镜装置即可实现。...

从快速荧光成像到商用显微镜下活细胞膜上的分子扩散规律研究

生物分子研究 体内 对于理解生物体系中分子功能至关重要。本文描述了一种新方法，可实现荧光生物分子（如 DNA 或蛋白质）向活体微生物内部的导入。对 体内...

通过电穿孔法将荧光生物分子导入活体微生物并进行观察

详细介绍了人脐静脉内皮细胞（HUVEC）培养、微管生长荧光标记物的瞬时转染、活细胞成像以及间期微管生长动态的自动化分析方案。...

活体人脐静脉内皮细胞中微管动态的高分辨率延时成像与自动分析

准确量化囊泡运输事件通常可为特定蛋白质的功能及突变效应提供关键见解。本文介绍了利用对pH敏感的GFP变体——超 ecliptic...

pHluorin在芽殖酵母内吞作用的定量、动力学及高通量分析中的应用

我们介绍了如何利用微操作和光操控技术（如荧光漂白恢复技术FRAP和光激活技术），来测定迁移细胞内蛋白质的运动参数及其时空动态。实验读出结果包括运动调控因子或其下方肌动蛋白细胞骨架的亚细胞动态行为与更新速率。...

允许分析细胞骨架调节因子形态发生动力学与周转的显微操作技术

本方案描述了一种基于荧光涨落光谱的方法，用于研究介导细胞间相互作用的蛋白质之间的相互作用。 即...

一种用于检测细胞间接触处蛋白质-蛋白质相互作用的荧光涨落光谱分析方法

三维单分子定位显微技术被用于研究活体细菌细胞中荧光标记蛋白的空间位置和运动轨迹。本文所述的实验与数据分析方案基于整合的单分子轨迹，确定胞质蛋白的主要扩散行为。...

单分子追踪显微技术——用于确定细胞质分子扩散状态的工具

本研究概述了利用免费提供的图像处理工具，对C. elegans中的突触大小与定位、肌肉形态以及线粒体形态进行定量测量的方法。该方法可使未来在 C.

用于研究细胞结构与细胞器形态的定量方法 秀丽隐杆线虫

本文介绍了一种利用宽场荧光显微镜对基于移动单分子荧光共振能量转移（smFRET）探针进行时空分析的方法。新开发的软件工具包能够确定运动探针的smFRET时间轨迹，包括随时间变化的准确FRET效率和分子位置。...

自动化二维时空分析可移动单分子FRET探针

学习使用AMEBaS自动提取极化单细胞中线并分析荧光延时数据。该方法通过自动化消除主观偏差，实现高通量细胞动态定量，适用于花粉管、根毛等多种细胞类型的时空动态研究。

AMEBaS：极化单细胞比率荧光延时成像的自动中线提取与背景扣除

探索更多视频