【Nat.Biotechnol.】突破单光子极限!清华戴琼海院士团队推出EFLIM技术,将荧光寿命成像灵敏度提升超2个数量级
文章标题:High-fidelity fast fluorescence lifetime imaging by event-based denoising
通讯作者:Jiamin Wu (吴嘉敏), Qionghai Dai (戴琼海)
文章概要
清华大学自动化系戴琼海院士在荧光成像领域取得重大突破,针对传统荧光寿命成像 microscopy (FLIM) 长期受制于高光子需求、难以实现活体深层组织快速成像的科学瓶颈,提出了一种名为基于事件的第一光子荧光寿命成像 (EFLIM) 的自监督去噪方法。该研究创新性地放弃了传统的每个像素累积光子计数直方图的固有范式,而是将每一次脉冲激发过程精确建模为离散的二元事件。在极弱光照条件下,EFLIM成功实现了在平均每像素光子数低于1个的极端场景下对荧光寿命的精准反演,将FLIM的成像灵敏度直接提升了两个数量级以上。这一通用成像框架在攻克光毒性限制、抗干扰活体深层组织脉冲记录等方面展现出巨大潜力,为细胞生物学、免疫学以及病理诊断提供了强有力的全新数字化观测工具。
引言
荧光寿命作为荧光分子在辐射跃迁前处于激发态的平均停留时间,是一项对微环境变化高度敏感的物理量。不同于传统的荧光强度成像,荧光寿命成像 microscopy (FLIM) 能够摆脱探针浓度波动、激发光不均匀以及光漂白等物理干扰,从而对细胞内的微环境变化、分子间相互作用以及蛋白质构象演变进行定量测量。在神经科学、结构生物学及细胞代谢监测等前沿领域,FLIM已经逐渐成为不可或缺的分析手段。
然而,受限于荧光发射本身固有的量子随机特性,单分子在每次激发后的光子发射时间存在显著涨落。为了在离散的直方图中构建出符合指数衰减规律的统计曲线,传统分析算法通常要求每个像素必须收集数百甚至数千个有效的发射光子。这一严苛的光子预算直接限制了成像的采样速度与探测效率,在面对极其脆弱的活体标本或处于深度组织的微弱荧光信号时,高强度的激发光极易引发严重的光毒性和光漂白,进而导致常规算法彻底失效。尽管近年来涌现出许多基于监督学习的深度学习去噪网络,但在平均每像素光子数小于1的极稀疏光子预算下,由于无法建立鲁棒的直方图统计特征,且受困于监督学习模型的泛化性不足,这使得极弱光环境下的高保真荧光寿命观测依然是一个巨大的科学挑战。

Fig. 1: Principle of EFLIM.
a, In conventional photon-counting FLIM, each pixel is repeatedly excited to generate abundant photons with distinct arrival times. These photons are grouped into discrete time bins to form a histogram, from which fluorescence lifetime is estimated through methods such as LS fitting. b, Under sparse excitation with only a few PPP, histograms cannot be reliably constructed, leading to failure in conventional fitting methods. c, EFLIM models each excitation event as a binary process, detecting either no photon or the first-arrival photon. d, EFLIM leverages information from temporally adjacent frames to enable self-supervised network training without the need for ground truth. e, During inference, EFLIM uses a sliding 3D window to estimate lifetime of the central frame. f, The overall workflow of EFLIM. Intensity and lifetime are processed from the raw photon data, followed by an integration for visualization.
主要实验及结论
为了从根本上解决极弱光下的寿命估计难题,研究团队独辟蹊径地提出了EFLIM框架。如图1所示,团队不再使用累积直方图,而是将飞秒激光的每一次脉冲激发电极化为一个离散的二元过程,即要么完全没有光子被检测到,要么仅记录下第一个到达光子的精准时间戳。在自监督网络训练的物理逻辑上,研究人员充分利用了生物样品在空间和时间维度上的高相关性先验,同时确保光子噪声在时空上的独立性,构建了类似于Noise2Noise和Noise2Void的自监督学习模型。该模型能够直接把三维的时空稀疏光子序列映射为高保真的寿命图像。在技术细节上,为了防范空像素的无效信息对网络参数造成系统性干扰,EFLIM的损失函数被精细配置为仅针对检测到有效光子的像素起作用,实验证实 omission 这一掩膜操作会导致寿命估算发生严重的系统性偏低。

Fig. 2: EFLIM achieves over two orders of magnitude improvement in sensitivity compared with existing methods.
a, Cross-section of pine stem with its intensity and lifetime distributions. A total of 200,000 pulses per pixel were used for excitation. Lifetimes were estimated using CMM. The white dashed box indicates the region enlarged in b. The lifetime profile along the indicated line is shown. Scale bar, 100 μm. b, Comparison of lifetime estimation by CMM, FPFLI and EFLIM. We randomly selected different numbers of excitation pulses to demonstrate varying photon budgets. Scale bar, 50 μm. c, Comparison of lifetime SNR under different photon budgets. Faded dots represent n = 10 independent training runs; solid lines and error bars show the mean ± s.d. d, Phasor plots at 10 and 10,000 pulses per pixel. Conventional phasor clusters were dispersed under low photon counts (top), whereas EFLIM reduced noise and sharpened cluster separation (bottom).
为了量化评估这一算法在实际硬件平台上的技术优势,研究人员使用搭建的二光子荧光寿命成像系统对松茎横切面标本进行了对比测试。如图2所示,在收集了高达20万次脉冲的超高信噪比参考图像中,传统的中心矩算法能够清晰分辨微环境差异带来的寿命空间异质性。随后,研究团队通过对该原始数据集进行随机抽样,模拟了不同光子预算下的成像质量。实验结果令人瞩目,即便在平均每个像素仅分到10个激发脉冲(对应每像素仅0.025到0.15个有效光子)的极端稀疏条件下,EFLIM预测的荧光寿命图依然能够保持高达17分贝的图像信噪比,其图像质量和相量图的聚类清晰度甚至超越了常规算法在1万个脉冲预算下的表现,充分印证了其提升两个数量级灵敏度的绝对技术优势。

Fig. 3: EFLIM enables high-speed deep-tissue lifetime imaging with strong robustness to intensity fluctuations.
We recorded calcium activity in awake mice expressing GCaMP6s, a genetically encoded calcium indicator that exhibits large intensity changes but negligible lifetime variation upon calcium binding. Motion introduced additional intensity fluctuations. a, Lifetime estimated by CMM, FPFLI and EFLIM. The white dashed box highlights a representative cell undergoing a calcium event, which is enlarged below as a time series. Scale bars, 50 μm (main) and 10 μm (enlarged views). b, Raw soma-averaged intensity trace corresponding to the CMM analysis. c, Raw soma-averaged intensity trace corresponding to the FPFLI analysis. d, Soma-averaged intensity trace after denoising, corresponding to the EFLIM result. e–g, Soma-averaged lifetime trace estimated using CMM (e), FPFLI (f) and EFLIM (g). h, PPP of a single pixel, after intensity denoising. i, Lifetime of a single pixel in the soma, estimated by EFLIM. j, Comparison of soma pixel lifetimes estimated by CMM, FPFLI and EFLIM at a single frame. Data are the mean ± s.d. The average number of photons collected in 3,470 frames is plotted as a red dashed line.
在攻克极弱光适应性后,EFLIM展现出极强的抗强度伪影干扰能力。在活体深层组织成像中,散射、组织吸收和动物运动往往会给强度信号带来剧烈抖动。如图3所示,研究人员对清醒小鼠的大脑皮层深层神经元进行了高速钙离子寿命成像。由于实验采用的钙指示剂GCaMP6s在结合钙离子时仅发生强烈的荧光强度跃迁而荧光寿命基本维持恒定,加之小鼠呼吸与微动带来的伪影,传统分析算法在单像素级别输出了剧烈震荡的错误寿命轨迹。相反,EFLIM成功解耦了寿命估计与强度变化,有效抑制了强度伪影,在细胞体及单像素层面上均输出了一条水平稳定的寿命曲线,证明了其在动态活体成像中的极高鲁棒性。

Fig. 4: Quantitative intracellular calcium signaling.
a, Lifetime estimated by CMM, FPFLI and EFLIM. The white dashed box highlights a representative cell undergoing a calcium event, which is enlarged below as a time series. Points 1 (orange) and 2 (blue) mark locations used for temporal analysis shown in b–d. Scale bars, 20 μm (main) and 5 μm (enlarged views). b–d, Lifetime traces at points 1 and 2 over 60 s using CMM (b), FPFLI (c) and EFLIM (d).
这种高时空分辨与定量无创的特性,使得EFLIM在捕捉细胞内瞬态化学信号方面大放异彩。如图4所示,在评估HeLa细胞受凝血酶受体激动剂肽诱导的内源钙离子释放实验中,每像素仅能分配到区区零点几个光子。在此场景下常规算法完全丢失了寿命波动的物理信号,而EFLIM则以高清视觉效果重构出了荧光寿命从2.5纳秒到3.3纳秒的动态跃迁过程。得益于极高的灵敏度,EFLIM甚至首次在胞内不同微小区域之间捕捉到了明显的空间传递延迟,这与局部钙微结构域先于全局激活的生物学理论完美契合。

Fig. 5: Multiplexed imaging of lymphocytes during immune response with a single spectral channel.
a, Experimental setup. T cells were labeled with CFP and B cells were labeled with PA-GFP; photons below 700 nm were collected in a single spectral channel. b, Lifetime-resolved imaging of GC response in inguinal lymph nodes over 2 h, showing densely aggregated GC B cells (orange) surrounded and intermittently infiltrated by follicular helper T cells (green). Scale bar, 100 μm. c, Cell velocities at T = 1 min. Scale bar, 100 μm. d, Lifetime histograms distinguishing three populations: T cells (green), B cells (orange) and a group with low lifetimes (blue). e, Instantaneous speed plotted against lifetime, showing faster migration of T cells relative to B cells. f, Time series of a representative T cell secreting a vesicle and subsequent contacts with other T cells. Scale bar, 10 μm. g, Lifetime decay of secreted vesicles, with spots marking contacts with other T cells. Vesicle lifetimes decrease from ~2.3 ns to ~1.6 ns over 25–30 min. h, Spatial distribution of vesicles, with triangle size indicating residence time. Vesicles were predominantly in follicles and only rarely detected inside the GC. Scale bar, 100 μm. i, Comparison of track-averaged speeds and lifetimes. Each data point represents one track; n = 1,116 tracks in total, classified by track-averaged lifetime as vesicles (<2.1 ns), T cells (2.1–2.4 ns) and B cells (2.4–3.3 ns). Box plots show the medians, 25th and 75th percentiles, whiskers extending to 1.5× the interquartile range and outliers as crosses. Statistical significance was assessed between vesicles and T cells using a two-sided Welch’s t-test (P < 0.0001, Cohen’s d = 2.51).
不仅如此,EFLIM还展现出强大的单通道多组分多路复用能力。如图5所示,在为期两小时的生发中心免疫反应长时程追踪实验中,研究团队仅通过单 spectral 通道收集光子,利用EFLIM算法将T细胞和B细胞基于其特异性的荧光寿命特征进行了精准的分离与伪彩标注。实验不仅清晰刻画了免疫细胞相互缠绕碰撞的动力学轨迹,更意外发现并追踪了一类由T细胞释放的微小囊泡结构。这些囊泡在脱离母体后,其荧光寿命在半小时内表现出从2.3纳秒向1.6纳秒持续缩短的线性特征,空间分布高度集中在卵泡区域,这一发现为研究免疫微环境中的旁分泌介导机制提供了直观的微观物理证据。

Fig. 6: High-speed label-free fluorescence lifetime imaging of human glioma tissue with histological validation.
a,b, Label-free fluorescence lifetime maps of a human glioma section acquired using CMM (a) and EFLIM (b). For CMM, tenfold more photons are collected. Scale bars, 1 mm. c, Corresponding H&E staining of the same tissue, used for morphological reference. Scale bar, 1 mm. d, Lifetime maps estimated using CMM with different numbers of excitation pulses per pixel (700, 2,100, 12,600 and 37,800) at the center FOV. Scale bar, 100 μm. e, Lifetime map reconstructed by EFLIM using 70 pulses per pixel at the center FOV. Scale bar, 100 μm. f, Enlarged CMM lifetime map of one boxed region in a. g, Enlarged EFLIM lifetime map of the corresponding boxed region in b. h, Enlarged H&E image of the corresponding boxed region in c. i, Enlarged CMM lifetime map of another boxed region in a. j, Enlarged EFLIM lifetime map of the corresponding boxed region in b. k, Enlarged H&E image of the corresponding boxed region in c. EFLIM preserved tissue-scale structural features and reveals lifetime variations that correspond to morphological differences observed in H&E. Scale bars, 100 μm. l, Higher-magnification EFLIM lifetime map of the dashed box region in g. m, Corresponding H&E image from h. n, Higher-magnification EFLIM lifetime map of the dashed box region in j. o, Corresponding H&E image from k. p, Higher-magnification EFLIM lifetime map of another dashed box region in j. q, Corresponding H&E image from k. Scale bars, 50 μm.
最后,为推动该技术的临床转化,团队探讨了其在医疗无标记诊断中的应用可能。如图6所示,研究人员针对人类胶质瘤组织切片进行了大面积的内源性自体荧光寿命成像,并与后续的HE染色组织病理学金标准进行了严格的形态学空间对齐。EFLIM在单视场仅需0.33秒的高速扫描下,便精准还原了基质区、肿瘤核心区以及坏死边界的复杂异质性代谢轮廓,其血管壁的长寿命特征和红细胞的极短寿命特征皆被完美保留,证明该方案完全具备在手术切缘评估等时间敏感型临床场景中落地的可行性。
总结及展望
综上所述,自监督EFLIM技术的诞生彻底颠覆了传统的直方图拟合思维,通过将光子到达事件转化为精细的二元概率逆问题求解,首次在单光子级预算下实现了兼具高时间分辨率与空间保真度的荧光寿命成像。这项技术大幅降低了对昂贵高效率探测器的硬件依赖,极大扩展了生物大分子长时程三维活体动态观测的科学边界。
展望未来,尽管EFLIM目前作为一阶矩估计器,对于复杂的多组分混合衰减和极端高背景噪声环境仍需结合相量分析法等数学工具进行校正与分离,但该架构表现出的超高光子利用率,已经为其在眼底视网膜无创超微成像、胚胎发育长时程安全监测、超分辨单分子寿命定位以及光纤内窥诊断等光敏感及空间几何受限的前沿医学场景中应用奠定了坚实的基础,有望成为跨越生物物理与临床医学的下一代核心成像底座。