北京大学未来技术学院席鹏团队开发保真优先的荧光显微镜反卷积方法
- 2024-09-06 08:46:253524
来源:北京大学
【仪器网 科技成果】荧光显微镜是活细胞成像中的重要工具。超分辨率显微镜的提出进一步打破了荧光显微镜的衍射极限,使得观察更为精细的细胞结构成为可能。然而,成像过程中的随机噪声以及光学模糊影响了荧光成像质量。在活细胞成像中,高光子剂量对细胞状态造成影响,并且会导致光漂白,增加了获得高信噪比图像的难度。在超快成像中,曝光时间不足也制约着高信噪比图像的获得。综合而言,荧光成像中的成像速度、持续时间和分辨率等关键参数受到了极大限制。
反卷积算法(deconvolution)在荧光成像领域得到了很大应用,因为它们具有通过以低光子预算从丰富的噪声中恢复荧光信号来提高成像持续时间和时空分辨率的潜力。尽管算法补偿在增强这些荧光成像的关键方面显示出了巨大的潜力,但其保真度,尤其是在分辨率提升方面仍然受到质疑。此外开发能够工作在更低信噪比的反卷积算法,可以进一步地提高成像时程以及速度,为发现新生物现象提供更多可能。
为了解决上述问题,北京大学未来技术学院席鹏教授团队重塑了荧光成像反卷积算法的噪声控制模型以及分辨率提升机制,提出了多分辨率(Multi-resolution analysis, MRA)反卷积算法。该方法有效提高了反卷积算法可以恢复的最低信噪比极限,并解决了传统反卷积算法在分辨率提高时容易产生伪影的缺陷。该文章以“Multi-resolution analysis enables fidelity-ensured deconvolution for fluorescence microscopy”为题发表在eLight。
研究者首先分析了传统的基于空间变化的全变分(total variation)以及海森(Hessian)正则的缺陷,认为虽然其提取的连续性特征部分符合荧光成像的特点,但其也会混淆荧光信号与背景边缘,并且难以从极高噪声中提取信号。根据荧光成像特异性标记以及生物结构的特点,研究者首次提出两个基本属性来区分荧光信号和噪声:(1)跨边缘的高对比度,以及(2)沿边缘的高连续性。在数学实现上,作者利用多尺度framelet以及curvelet变换来提取荧光信号中的这两个关键特征,称为MRA正则化方法。
在分辨率提升机制方面,研究者首先通过模拟实验证明了Richardson-Lucy迭代由于其统计估计特性会自发产生伪影。基于物理模型的迭代逆滤波虽然可以在无噪声情况下完美复原模糊图像,但其高频细节恢复需要大量迭代次数,容易被噪声淹没且与基于空间变化的噪声控制正则相容较差。而研究者所提出的MRA正则化方法由于提供了更精确的荧光噪声模型,可以有效控制逆滤波迭代中的噪声并保留其推断的高频细节。研究者进一步提出了一种加速FISTA迭代的策略,可以在短计算时间内完成上千次约束逆滤波迭代来恢复高频细节。
通过多种成像条件的高低信噪比对图像恢复实验,研究者验证了MRA正则化方法相较基于空间变化的正则方法可以将信噪比进一步提高2—10dB。
图1. MRA反卷积实现更优的噪声控制和保真分辨率提升。a. 利用SIM高低信噪比对MRA 噪声控制的验证; b. MRA与Hessian正则在连续帧图像处理中的对比; c. U2OS细胞线粒体的宽敞图像、Huygens、Sparse以及MRA解卷积结果和作为真值参考的SR-SIM结果; d. 从开源BioSR数据集中获得的线性非线性SIM图像对,非线性SIM作为Huygens、Sparse以及MRA对线性SIM解卷积结果的比对真值;e. 沿d中两个白色箭头指示的线的强度分布
研究者利用宽场-SIM、线性SIM-非线性SIM、Confocal-SoRa以及Confocal-STED多种高低分辨率数据对,验证了MRA反卷积相较以Richardson-Lucy为分辨率提升机制的反卷积算法具有显著优越的保真性,解析的高频细节结构与超分辨显微镜几乎一致。
研究者进一步提出了层切多分辨率分析反卷积算法(SecMRA),引入了全新设计的偏置稀疏正则,该方法可以有效地提高在强背景图像处理中的表现。SecMRA反卷积后的宽场图像可以达到类似共聚焦显微镜的层切效果,可以在提供高质量图像的同时进一步降低活细胞成像的光毒性。作者在宽场、Lightsheet、Confocal、SoRa、SIM、STED和Expansion显微镜等多种成像模态中验证了该方法的有效性。
研究者展示了利用MRA反卷积算法来支持低光照低毒性荧光成像以研究细胞活动的实例。利用MRA反卷积,研究者实现了对微管和溶酶体共标的长时程超分辨成像,观察到了溶酶体运动与微管结构的紧密联系。作者还利用宽场显微镜实现了细胞核周半小时以上的线粒体和内质网共标成像,观察到了线粒体迅速延伸后收缩最后绑定在内质网上的现象。
为了推广该技术,研究者还提出了解决正则反卷积算法中超参数调整的方案,通过用曲线波系数稀疏度来估计噪声强度,实现了参数自动确定。研究者编写了用户友好的GUI软件,并且开源了所有代码以及荧光图像原始数据,以便于用户和开发者使用。研究者希望该工作可以为便于生物学家和显微学家可以一手使用的反卷积工具。
图2. SecMRA反卷积助力长时程活细胞成像。a、b. 微管溶酶体双色成像,SecMRA反卷积有效去除图像中的噪声、离焦背景并提高分辨率;c—g. SecMRA反卷积实现内质网与线粒体互作成像
北京大学未来技术学院博士生侯宜伟为论文的第一作者。席鹏、北大李美琪博士为该论文的共同通信作者。论文主要贡献者还包括北京大学博士生付允哲以及北京艾锐精仪公司工程师王文熠、盖希川。该研究工作得到了国家自然科学基金区域联合重点项目、面上项目等的大力支持。