自动化ROI检测实现细胞培养中数万个突触的突触活动快速量化

《Frontiers in Synaptic Neuroscience》：Automated ROI detection allows rapid quantification of synaptic activity across tens of thousands of synapses in cell culture

钙离子成像（Calcium imaging）在无Mg2+且含TTX（tetrodotoxin，河豚毒素）的培养神经元中可显示NMDAR（N?methyl?D?aspartate receptor，N?甲基?D?天冬氨酸受体）依赖的自发突触钙瞬变（spontaneous synaptic calcium transients），其反映突触前及突触后功能。本研究介绍一种高通量自动化分析流程，结合Suite2p ROI（region of interest，感兴趣区域）检测与Python脚本，对数万个突触进行分析并量化突触前囊泡融合速率（频率）、突触后功能（幅度）及功能性突触数目。研究人员利用该流程测试已知NMDAR激动剂（glycine，甘氨酸）与拮抗剂（ketamine，氯胺酮；memantine，美金刚；APV，DL?2?amino?5?phosphonopentanoic acid）、调节突触前功能的化合物（PDBu，phorbol 12,13?dibutyrate）及抗NMDAR自身抗体（源自脑炎患者血清），证明该流程在单突触水平检测功能障碍方面较其他方法更敏感。利用此流程检测、追踪和量化数万个突触及数百万个突触钙瞬变的活动与变化，将助力于筛选保护突触功能的化合物之药物发现。

论文解读：《Automated ROI detection allows rapid quantification of synaptic activity across tens of thousands of synapses in cell culture》

发表于Frontiers in Synaptic Neuroscience

一、研究背景与目的

突触是神经元间信息传递的基本单元，突触功能障碍是认知疾病与精神障碍的常见诱因。传统电生理记录通量低、操作复杂；已有的自发突触钙瞬变（spontaneous synaptic calcium transients）分析方法依赖人工或半手动框选ROI（region of interest，感兴趣区域），费时费力且易受主观偏差影响，单实验条件通常仅能分析1000–2000个突触，难以满足药物筛选与大规模突触功能研究的需求。因此，研究人员旨在建立一套基于开源工具的自动化、高通量突触钙成像分析流程，实现对培养神经元中数万个突触的自发NMDAR（N?methyl?D?aspartate receptor，N?甲基?D?天冬氨酸受体）依赖钙瞬变的自动检测、追踪与量化（频率、幅度及活性突触数），并验证其对药理学干预及病理源性抗体的敏感性。

二、主要关键技术方法

研究人员取生后第0–2天（P0–P2）Wistar大鼠仔鼠原代皮质?海马混合神经元培养（primary corticohippocampal cultures），经AAV2/1?hSyn?GCaMP6f转导表达基因编码钙指示剂GCaMP6f。在含1 μM TTX（tetrodotoxin，河豚毒素）及0 mM Mg²⁺（去除NMDAR镁阻断）的ACSF（artificial cerebrospinal fluid，人工脑脊液）中，以宽场倒置荧光显微镜20 fps采集DIV19–22培养物影像。采用Suite2p进行自发荧光变化驱动的ROI检测并优化参数（sparse_mode=1，threshold_scaling=1.5，spatial_scale=0判定约6像素≈2 μm，关闭运动校正与去卷积，保留所有检出ROI），后续用Python脚本按荧光偏度（skew＜1）与紧凑度（compactness≤1.4）过滤剔除噪声及树突事件。荧光信号先以70%神经丛（neuropil）荧光减除背景，再用airPLS（adaptive iteratively reweighted penalized least squares）算法基线校正并归一化为ΔF/F₀，以MAD（median absolute deviation，中位数绝对偏差）估算SD（standard deviation，标准差，SD=MAD/0.6745），阈值设为F₀+4.5×SD，用SciPy中find_peaks检测峰（最小间距5帧、最小半宽3帧、 prominence≥F₀+2×SD）。过夜处理组用CellProfiler骨架化神经突起（neurite）长度归一化活性突触数。统计学分析采用R语言中glmmTMB包拟合GLMM（generalized linear mixed?effects model，广义线性混合效应模型），频率假设Gamma分布对数连接，幅度假设log?normal分布，成像区域为随机效应，并以10000次聚类Bootstrap重抽样验证。

三、研究结果

Suite2p与Python可实现突触钙成像分析的自动化

最大?最小投影图像中可辨自发点状钙瞬变，Suite2p可大量自动检出ROI，经偏度与紧凑度过滤区分突触与树突事件（树突事件compactness＞1.4或形态延长予以排除）。airPLS基线校正后ΔF/F₀轨迹可明确标记峰。假阴性率3.7±1.9%，假阳性率0.3±0.7%。每视野检出1700–3300个突触，数据集可达15–34万个个体突触、逾千万瞬变事件。

Phorbol esters增加，APV减少自发突触钙瞬变频率

急性给予PDBu（1 μM，促进Munc13/PKC介导的突触前囊泡释放）使频率由0.046 Hz升至0.11 Hz（2.34倍，p＜0.0001），活性突触数无显著变化；继加APV（50 μM NMDAR竞争性拮抗剂）使频率降至0.031 Hz（66%对照），活性突触数下降约7倍（p＜0.0001），幅度在PDBu处理后略降（0.38 vs 0.48 ΔF/F₀），APV不再进一步降低幅度——证实该钙瞬变为突触前谷氨酸释放触发、NMDAR介导的突触后Ca²⁺内流。

Glycine增加自发突触钙瞬变频率与活性突触数目

100 μM glycine（NMDAR共激动剂）使活性突触数由~2503增至~2832（p＜0.0001），频率由0.051 Hz升至0.079 Hz（1.55倍），平均幅度略降（0.43 vs 0.50 ΔF/F₀），表明NMDAR共激动可提升功能性突触检出率及事件频率。

AMPAR对自发突触钙瞬变贡献有限

10 μM NBQX（AMPA receptor拮抗剂）不改变频率（0.98倍，CI跨1），活性突触数有增多趋势（p=0.088），幅度轻微升高（0.52 vs 0.48 ΔF/F₀，p＜0.05），提示AMPAR不直接介导该Ca²⁺信号，阻断后可能因减弱NMDAR?AMPAR交互抑制或增大Ca²⁺驱动力使NMDAR来源瞬变幅度增大。

急性NMDAR靶向化合物ketamine与memantine降低突触活动

1 μM ketamine或1.5 μM memantine急性给药降低活性突触数（ketamine：~2179→~1769；memantine：~2376→~1427，均p＜0.0001）、频率（ketamine：0.70倍；memantine：0.63倍）及幅度（memantine显著降低至0.89倍，ketamine呈下降趋势）。随时间推移抑制效应渐增（与通道打开后阻滞特性相符），而APV与NBQX作用迅速且稳定。

过夜ketamine处理降低突触钙瞬变频率与幅度

过夜1 μM或10 μM ketamine不影响单位神经突起长度活性突触数（10 μM有减少趋势p=0.023 vs 1 μM），但10 μM显著降低频率（0.044 vs 0.058 Hz，0.76倍），1 μM与10 μM同等降低幅度（0.50 vs 0.65 ΔF/F₀，0.77倍）。

过夜memantine处理降低突触钙瞬变幅度

过夜1.5 μM与10 μM memantine不影响活性突触数及频率，但10 μM显著降低幅度（0.80倍，p＜0.05），与含1 mM Mg²⁺培养基中memantine亲和力较低致较弱作用相符。

突触钙成像检测患者源性NMDAR自身抗体的差异效应

过夜15 μg/mL抗NMDAR自身抗体003?102与008?218（vs对照抗体MGO?53）降低活性突触数（008?218显著减少，003?102边缘显著）、频率（分别降21%与23%）及幅度（均约0.93倍）；007?124（低亲和力）无显著影响——说明该流程可在功能水平区分不同自身抗体对突触的作用差异（内化vs构象封闭）。

突触钙成像作为高通量筛选的潜力

各层级变异中突触内变异最大，生物学重复间变异最小。以median/MAD计算Z'（Z?prime），频率/重复水平Z'=0.874，活性突触数/重复Z'=0.592，达细胞学高通量筛选可接受标准（0–0.5可用，＞0.5优秀），具药物筛选应用前景。

四、讨论与结论翻译

研究人员建立的自动化钙成像ROI检测与分析流程可报告包括功能性突触数目、突触前囊泡融合概率（频率）及突触后NMDAR功能（幅度）在内的多维度突触功能参数，在无人工选点下采集逾119万活性突触及逾1020万独立钙瞬变。该流程Z'值在重复与孔水平达高通量筛选可用范围。NMDAR激动/拮抗剂同时改变频率与活性突触数，而突触前调节剂（PDBu）不改变活性突触数，提示需联合活性突触计数以区分突触前/后效应。APV消除～92%事件验证NMDAR依赖性；NBQX不抑制频率但小幅升高幅度，符合AMPAR?NMDAR交互抑制模型。患者源性抗NMDAR抗体功能筛查揭示不同抗体异质性效应，优于固定组织免疫染色。该流程亦适用于阿尔茨海默病、精神分裂症、自闭症及帕金森病等突触病变模型的药物筛选与突触可塑性研究。

综上，本研究通过将Suite2p适配于突触尺度ROI检测并结合定制Python后处理，实现了培养神经元中数万突触NMDAR依赖自发钙瞬变的高通量、无偏、敏感量化，克服传统手动分析局限，为突触功能基础研究及突触保护性化合物的高通量药物发现提供了可靠平台。