本手稿使用基于 Fiji 的开源软件包 VirusMapper,对超分辨率显微镜图像进行单颗粒分析,以生成纳米尺度结构的精确模型。
本手稿使用基于 Fiji 的开源软件包 VirusMapper,对超分辨率显微镜图像进行单颗粒分析,以生成纳米尺度结构的精确模型。
超分辨率荧光显微镜目前正彻底改变细胞生物学研究。其突破约300 nm分辨率极限的能力,使得对纳米尺度生物复合物和生物过程的常规成像成为可能。分辨率的提升也意味着电子显微镜中常用的方法(如单颗粒分析)可直接应用于超分辨率荧光显微镜。通过将这种分析方法与超分辨率光学成像相结合,可以利用荧光显微镜对分子的特异性标记能力,生成亚稳态结构内分子组分的结构图谱。为此,我们开发了一种新型算法——VirusMapper,并将其封装为易于使用、高性能、高通量的ImageJ插件。本文对该软件进行了深入介绍,展示了其揭示生物分子复合物中新结构特征的能力。文中详细说明了如何准备兼容的数据,并提供了使用该算法对超分辨率图像进行单颗粒分析的逐步操作流程。
超分辨率(SR)显微技术通过实现对关键分子过程的成像,并结合理解这些过程所必需的分子特异性标记,对细胞生物学产生了重大影响。如今,SR使得光学显微镜能够达到以往仅电子显微镜(EM)可实现的分辨率水平(20–150 nm),同时保留了光学显微镜的主要优势,例如对活细胞进行成像的潜力1,2。此外,纳米尺度下存在的结构保守性使得单颗粒分析(SPA)可应用于SR数据,这一概念在电子显微镜中已被广泛应用3。利用SPA,可以对同一结构的多个高度保守拷贝进行成像并加以平均,从而提高所观察对象的分辨率、精确度或信噪比。当与SR结合使用时,SPA已被证明是高精度绘制核孔复合体4,5、中心体6以及病毒(如HIV7和HSV-18)组分分布的有力工具。
然而,由于缺乏可用的软件,SR 与 SPA 的常规联合应用一直面临挑战。为此,我们开发了 VirusMapper,这是对广受欢迎的图像处理软件 ImageJ/Fiji9 的一个插件。这是首个可免费获取的、适用于荧光图像的通用 SPA 软件包10,旨在为采用 SR 显微镜成像的结构提供快速、用户友好、多通道的朴素平均分析。尽管该软件最初为病毒研究设计,但也可应用于任何可通过成像识别并定位不同分子种类的大分子复合物。
VirusMapper 可用于生成任何已知结构的高精度分子模型,从而计算平均尺寸及其他参数。该算法设计特别适用于分离不同结构群体,能够确定不同的取向或不同的形态状态。此外,在基础结构已知的情况下,可采用多通道成像并利用参考通道,从而实现基于参考的结构发现。软件的下载与安装说明请访问 https://bitbucket.org/rhenriqueslab/nanoj-virusmapper。示例数据也可在该网站获取,建议用户先在示例数据上练习使用该软件,再尝试应用于自身数据。
本文描述了使用该插件从原始数据生成单粒子分析(SPA)模型的步骤。该软件以包含单标记或多标记结构的原始图像作为输入,根据运行过程中调整的多个参数,输出显示成像结构内标记组分平均分布情况的SPA模型。
本实验方案的目标是根据图1所示流程,生成精确的单粒子平均(SPA)模型,以反映成像结构中各组分的平均定位信息。如图1所示,软件工作流程可有效划分为三个阶段。第一阶段是对大尺寸图像进行分割,从而为每个通道生成粒子堆栈。这些粒子是用于创建模型以及生成模型种子的基本单位。第二阶段是生成种子图像,这些种子图像将在最终阶段用于对整套粒子进行配准。该过程通过选择一个参考通道,并在该通道中手动选取将用于生成种子的粒子来完成。种子在参考通道中选定,但可为所有通道生成。粒子首先通过在该通道中拟合二维高斯函数进行初始重对齐。所有被选中并完成重对齐的粒子将被平均,以生成一个种子。对于数据中需要建模的每一种常见结构,应选择能够清晰且准确代表该结构的粒子作为种子。此阶段的界面也有助于在数据中扫描和识别此类结构。
最后阶段是通过模板匹配生成模型。该过程通过将最初提取的颗粒与上一节生成的种子图像进行互相关配准来实现。选取一部分已配准的颗粒进行平均,若需要,可进一步迭代此过程以降低模型的均方误差。该子集的确定依据是设定与种子图像之间的最小相似性阈值,只有满足该阈值的颗粒才被纳入。当在多个通道中同时创建模型时,使用联合相似性,即各通道相似性的平均值。所得模型及其贡献颗粒可进一步用于后续分析。
注意:本方案及视频是对原始论文10的补充,原始论文更详细地描述了该软件包。建议读者仔细查阅该论文,以获得有关软件使用的更多指导。整个流程主要包括三个阶段:颗粒提取,即将大图像分割为单个颗粒;种子选择,在此阶段识别数据中的共同结构并进行比对以生成种子,这些种子将用于最后阶段;以及模型生成,在此阶段基于这些种子通过模板匹配对提取的颗粒进行比对,并对部分颗粒进行平均以生成单颗粒分析(SPA)模型。
1. 运行软件包前的准备工作
2. 提取颗粒
3. 选择种子
4. 生成模型
在此,我们以模式痘病毒——痘苗病毒为例来演示该软件的应用。痘苗病毒是结构最复杂的哺乳动物病毒之一,其呈砖形的病毒粒子大小约为 350 × 270 × 250 nm3,内部包裹约 80 种不同的蛋白质13,14。通过电子显微镜可分辨出三种亚结构:包含双链 DNA 基因组的中央核心;位于核心两侧的两个蛋白质结构,称为侧体;以及包裹在外层的单层蛋白脂质双分子层包膜15。由于其体积较大、结构复杂,且易于进行重组荧光蛋白标记,痘苗病毒是展示 VirusMapper 工作流程的理想体系。
使用此处所述软件,可对多种蛋白质在痘苗病毒粒子上的分布进行建模。对目标蛋白进行标记和成像,可能同时结合另一种已知分布的蛋白作为参考,然后按照所述方法使用该软件生成该蛋白在病毒颗粒上定位的平均模型。在本示例中,建模了两种蛋白:内层核心蛋白L4和主要侧体组分F17。
使用一种F17标记GFP、L4标记mCherry的重组痘苗病毒16。将纯化的病毒稀释于1 mM Tris(pH 9)中,并通过在室温下孵育30分钟将其结合到经洗涤的高性能盖玻片上。随后用含4%多聚甲醛的PBS固定样品20分钟。盖玻片立即用抗荧光淬灭封片剂封固于载玻片上。采用商用SIM显微镜进行SIM成像。选择包含数百个病毒的视野,使用561 nm(32 µm光栅周期)和488 nm(32 µm光栅周期)激光,结合5个相位移动和3个光栅旋转采集图像。使用sCMOS相机采集图像,并通过显微镜配套软件进行处理。根据使用相同成像参数采集的多色荧光微珠样品滑片对各通道进行对齐。SIM重构及通道对齐后,将图像在Fiji中打开并合并为单个图像序列。
使用 L4 通道作为参考,从图像中提取病毒颗粒,未应用任何高斯模糊,因为这些颗粒具有中心最大值。本实验中约提取了 15,000 个颗粒。
由于痘苗病毒的几何结构,其侧体根据病毒的取向呈现出明显不同的外观。我们观察到两种取向,分别可区分出一个或两个侧体,将其称为正面取向和矢状取向。
在“生成种子”阶段,通过浏览粒子列表,分别选择前额向和矢状向的种子(图4和图5);选择那些明显属于某一取向的粒子。以L4通道作为参考通道,对种子进行相互对齐。同样,无需进行高斯模糊处理。每种取向各选取5个粒子,并对其进行平均以生成种子。
基于这些种子为每种取向生成了模型。未使用参考通道或平方强度值。最初将最大迭代次数设置为1,最小相似性设定为在每种情况下包含约1,000个颗粒,从而使得每种取向的模型外观保持一致。随后增加最大迭代次数,以允许模型收敛。因此,针对两个通道中的两种取向生成了模型(图7).

图1: VirusMapper 工作流程。 该插件分为三个主要阶段。首先从大尺寸图像中提取病毒颗粒,然后从数据中半手动选择模板图像或种子,最后参照这些种子从数据中生成最终的单颗粒分析(SPA)模型。 请点击此处查看该图的放大版本。

图 2: “提取病毒结构”对话框。 选择“提取病毒结构”时,将弹出此对话框。应使用初始估计值填写参数,以实现最佳分割效果。随后可选择“显示预览”,以便预览感兴趣区域(ROI),并对参数进行精细调整。 请点击此处查看该图的放大版本。

图3: 设置提取参数。 在预览将要提取的感兴趣区域(ROI)后,调整ROI半径、ROI数量和最大ROI重叠度,以达到如图所示的效果。ROI略大于颗粒,所有颗粒均被包含在某个ROI内,且ROI之间可充分重叠,以便分离成簇的颗粒。请点击此处查看该图的放大版本。

图4: 生成模板匹配种子。 “生成种子”对话框(1)用于设置需指定的参数。参考粒子序列(2)允许用户浏览参考通道中的粒子。当在参考粒子序列中查看某个粒子时,可在重新对齐粒子预览窗口(3)中查看所有通道的重新对齐后的粒子。请点击此处查看该图的放大版本。

图5: 添加种子图像。 当向“要使用的帧”框中添加种子时,系统将显示所有种子的平均图像(4)以及所涉及的帧(5)。对话框中会推荐与当前平均种子相似的颗粒(6)。请点击此处查看此图的放大版本。

图6: “生成模型”对话框。 当选择“基于种子生成模型”时,将出现此对话框。应填入初始估计值以实现最优模型生成,并选择在计算过程中显示的模型生成步骤。随后可选择“显示预览”,以运行模型生成过程并对参数进行精细调整。请点击此处查看该图的放大版本。

图7: 使用VirusMapper生成的模型。 对L4核心蛋白标记mCherry、F17侧体蛋白标记EGFP的痘苗病毒颗粒,采用SIM成像,并按照本方案所述方法利用软件生成模型。两个方向(正面和矢状面)可通过侧体的形态加以区分。比例尺 = 100 nm。 请点击此处查看该图的放大版本。
通过该方法,研究人员能够结合单粒子分析(SPA)和超分辨率显微技术的力量,以生成病毒及其他大分子复合物蛋白质结构的高精度、多通道二维模型。然而,一些重要的因素需要加以考虑。
应选择能够代表数据中一致存在的结构的种子。因此,在选择种子之前,必须仔细检查原始数据。这对于防止产生有偏倚的模型至关重要。可以通过检查为包含特定数量粒子所需达到的最低相似性阈值来验证种子的选择。显然,对于某一特定种子,若要包含给定数量的粒子所需设定的阈值越高,则说明该结构在数据中越明显。
当数据存在异质性时,模板匹配的概念尤为有用。所有可见的不同结构都应被识别,并为每种情况创建不同的模型。通过在一个通道中分离异质结构,同时在另一个通道中创建模型,可能揭示出原本并不明显存在的模式。
使用该算法时还需注意,迭代过程会最大化随机不对称性。例如,在对具有两个对称极大值的结构进行建模时,迭代过程中这两个极大值之间的所有微小不对称性将被相互对齐,最终模型因而会呈现最大程度的不对称性。如果这与所建模结构中已知的对称性不符,则应予以考虑。目前,避免这种最大化效应的唯一方法是将迭代次数限制为1次,但未来可能的改进方向是使VirusMapper在模型生成过程中纳入对称轴。VirusMapper的任何新版本都将在引用的网站上提供(参见材料表)。用户还可在此找到常见问题解答,以回应常见疑问。
所描述的软件适用于任何能够以足够分辨率成像的结构,以便可视化用户希望建模的特征。尽管单颗粒分析(SPA)可以提高分辨率,但它显然无法增强原本不可见特征的可见性。因此,本方案并非用于提高数据质量的方法。与任何技术一样,仔细的样品制备和成像策略的优化将提供最清晰的数据和最佳的最终模型。
超分辨率成像方式的选择同样重要,通常取决于具体的样本。VirusMapper 已经被验证可与 SIM 和 STED10 良好兼容,也可用于高质量的定位显微镜数据,但在这种情况下需格外注意,稀疏标记可能导致类似于不对称性最大化所产生的问题。
目前,VirusMapper 是唯一可免费获取的用于荧光图像单颗粒分析的算法,也是唯一通用的二维单颗粒分析平均软件。其他采用相同原理的研究4,6,8均使用了为特定研究定制的专用软件。尽管已有通用算法用于三维数据重构的发表5,18,但均未提供相应的软件。
按照本文所述方法使用时,VirusMapper 可用于生成病毒及其他复合物大分子蛋白质结构的精确、准确且稳健的模型。利用这些模型,研究人员能够对所研究结构的平均尺寸进行精确测量,从而有可能得出原本无法获得的生物学结论。
此外,利用该技术的多通道能力,可以对复合物内无限数量的蛋白质和组分进行定位,并发现新的蛋白质组织模式。在不同的生物学相关条件下(例如病毒生命周期的不同阶段)研究纳米尺度结构的变化,有望为生物学提供有价值的见解。
作者无任何利益冲突需要披露。
我们感谢 Corina Beerli、Jerzy Samolej、Pedro Matos Pereira、Christopher Bleck 和 Kathrin Scherer 对 VirusMapper 最初开发与验证工作的贡献。我们还要感谢 Artur Yakimovich 对手稿的审阅。本研究由生物技术和生物科学研究中心基金(BB/M022374/1)(R.H.);伦敦大学学院医学研究理事会分子细胞生物学实验室的核心资助(J.M.);欧洲研究理事会基金(649101-UbiProPox)(J.M.);以及医学研究理事会基金(MR/K015826/1)(R.H. 和 J.M.)资助。R.G. 由工程与物理科学研究理事会(EP/M506448/1)资助。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| Fiji | 开源图像分析软件 | ||
| NanoJ-VirusMapper | 由 Henriques 实验室开发 | 开源 Fiji 插件(https://bitbucket.org/rhenriqueslab/nanoj-virusmapper) | |
| VectaShield 抗荧光淬灭封片剂 | Vector Labs | H-100 | |
| Elyra PS1 | Zeiss | ||
| ZEN BLACK | Zeiss | SIM 图像处理软件 | |
| 高性能盖玻片 | Zeiss | 474030-9000-000 | |
| TetraSpeck 微球 | ThermoFisher | T7279 |