Summary
这篇手稿描述了如何区分不同的线虫使用远场衍射签名。我们比较了139种野生类型和 108 "滚子" C. 线虫的运动, 通过平均频率与瞬时霍夫费衍射签名在一个位置使用连续波激光。
Abstract
这篇手稿描述了如何使用世俗远场衍射特征对线虫进行分类。单个C. 线虫悬浮在光学试管内的水柱中。632 nm 连续波氖激光器是通过试管使用前表面反射镜。在光线穿过试管后, 至少有20-30 厘米的距离可以确保一个有用的远场 (弗劳恩霍夫) 衍射图样。当线虫在激光束内游动时, 衍射模式会发生实时变化。光电二极管被放置在衍射图样的偏心位置。光电二极管的电压信号是实时观测的, 并使用数字示波器记录。此过程重复139野生类型和 108 "滚子" C. 线虫。野生型蠕虫在溶液中表现出快速振荡模式。"滚子" 蠕虫有一个突变的角质层的关键组成部分, 干扰平滑的运动。不含饱和度和不活动的时间间隔将被丢弃。将每个平均值除以其最大值以比较相对强度是可行的。每个蠕虫的信号是傅里叶变换, 使每个蠕虫的频率模式出现。每种类型的蠕虫的信号都是平均的。对野生类型和 "滚子"线虫的平均傅里叶谱进行了明显的区别, 揭示了两种不同的蠕虫病毒的动态蠕虫形态可以用傅立叶分析来区分。每个蠕虫应变的傅里叶谱匹配一个近似模型使用两个不同的二进制蠕虫形状对应于运动时刻。实际和模拟蠕虫的平均频率分布的包络确认模型与数据匹配。这种方法可以作为许多微观物种的傅里叶分析的基线, 因为每个微生物都有其独特的傅里叶谱。
Introduction
该方法比较了两种不同运动模式的菌株的实验和模拟的频率谱. 结果表明, 当线虫在水柱中游动时, 频谱依赖于时间变化, 因此不需要清晰的显微图像进行分析。此方法允许定量 real-time 分析, 并提供与传统显微镜获得的图像/视频的补充信息。弗劳恩霍夫衍射, 也称为远场衍射, 为获得活衍射数据提供了依据1,2。在衍射图样的任何一个点上的光强度是从线虫的轮廓中的每一点的叠加光的结果3。因此, 随着时间的推移收集的光强度携带的信息关于线虫的运动。分析了随时间推移的衍射信号可以识别出相应突变体的特征运动, 因为分析了运动中涉及的所有频率, 补充了传统的视频分析。在这种情况下, 通过比较两种不同菌株的频率谱, 证实了 "滚子" 与野生型线虫的运动之间的特征差异。
一些以前的特征已被证实使用频率分析的衍射信号, 如游泳频率2,4。更重要的是, 这种方法可以作为传统显微镜的一种补充方法, 在收集数据的同时, 在计算机屏幕上实时观察运动。通过考虑衍射信号的傅里叶变换信号, 可以对具有不同运动模式的蠕虫的频谱进行量化。
Fourier-based 衍射在这项工作中的多学科性质涉及生物学和物理领域。在生物的晶体结构的研究中, 长期使用的衍射在生物学5和其他领域。然而, 在这个实验中, 采样6,7创建了远场衍射模式, 使生物体以激光束为中心。采样通常用于与重建原始对象的图像的相位检索算法相结合的无透镜成像8 。当散射像线虫那样存在时, 相位检索很难实现。时间衍射特征足以评估蜗杆运动的关键频率。这种方法是较少的计算征税和提供一种光学的方式来量化的运动。这项技术可以很容易地被改编为分析突变或环境条件, 改变行为。
Protocol
1. C. 线虫 生长和维护
- 准备线虫培养皿。
- 用琼脂溶液填充培养皿, 然后让它们凝固, 然后用 OP50 菌株的 大肠杆菌 区域性进行播种 9 , 10 .
- 在每个板上准备一个成年线虫的起始种群, 将几个成虫移动到新鲜琼脂填充的培养皿中, 并带有一个 大肠杆菌 补丁。保持20和 #176 的线虫培养; C 在孵化器中.
注意: 线虫菌株可从 线虫线虫 基因组中心获得。为这项研究, 狂放的类型, N2, 劳损和 OH7547 (otls199 [猫-2:: GFP + rgel-1 (F25B3.3)::d 失调 + 卷-6 (su1006)]) 应变, 陈列滚轮表型, 被使用了. - 为未来的区域性传播蠕虫。
- 将包含 C. 线虫 的 petri 盘和温度控制箱中未使用的生境培养皿移除。把它们放在解剖显微镜的舞台上.
- 点燃本生燃烧器, 并通过将金属放入火焰中以使其呈红色, 来杀菌铂线虫的选择。允许选择冷却到室温。不要设置拾取或让选择接触到污染物.
- 轻轻触碰尖端的拾取到细菌圈的边缘。这种物质是粘性的, 将使采摘个别成年线虫更容易.
- 将多达4的妊娠线虫转移到线虫生长培养基 (NGM) 琼脂填充的培养皿中, 并在20和 #176 孵育;蠕虫将产卵, 将在四天内成熟.
- 在移动四成虫线虫后, 将剩余的线虫返回孵化箱.
2。光学设置 ( 图 1 )
- 在光学工作台的后左角附近保护氦氖激光器, 并将其连接到电源.
注: 激光和 #39 的光束必须满足采样的要求。 C. 线虫 的长度约为1毫米, 因此, 在线虫发生的情况下, 激光束的直径应大于2毫米, 但不大于5毫米, 因此不难找到衍射图案. - 在氦氖激光器和样品之间放置一个中性密度的过滤器, 使激光光束在到达样品之前通过过滤器进行传输.
- 使用两个前表面铝制转向镜, 通过在中性密度过滤器后保护第一个镜像来建立潜望镜。将第二个镜子固定在第一个镜子下面10厘米的地方, 给出空间来引导激光束, 并在镜子之间插入试管。将激光束和试管对准, 使激光束通过试管垂直移动.
注意: 从衍射生物体到光电二极管的距离必须比生物体本身要大得多, 才能达到远场衍射。在这个实验中, 从试管到光电二极管的距离是20厘米. - 将光电二极管直接从第二个镜像中保护, 并将其传感器朝向镜像.
注: 含有线虫的试管将用化学钳放置在两个镜子之间。参见4和5节. - 将一个充水的试管放在支架上。调整展位高度。调整镜1和镜2的高度和角度, 使激光光束通过靠近试管的方向移动, 而不是直接对准光电二极管.
- 使用一个级别来确保展台形成试管的平整曲面。如有必要, 请进一步调整镜像.
- 使用数字示波器提供的 USB 电缆将光电二极管连接到数字示波器。将数字示波器连接到将用于记录和保存数据的计算机.
3。示波器设置
- 使用计算机上的示波器软件, 将采样速率设置为至少 8 Hz, 以充分解决蠕虫的颠簸循环.
注意: 采样率应该比预期的抖动频率的两倍还多, 因此奈奎斯特定理 11 满足.
4。准备数据收集的蠕虫和试管
- 将四条成虫线虫转移到新鲜的 NGM 10 琼脂填充的 Petri 板, 使用薄的扁平铂线拾取 (见第1.3 节).
- 从包装上卸下一次性塑料试管, 小心只触及其脊边上的试管.
- 使用微将蒸馏水注入试管, 直到试管大约80% 的蒸馏水.
注: 在处理线虫时, 只有使用蒸馏水或电离缓冲器, 如 M9 10 或磷酸盐缓冲盐水 (PBS), 因为自来水中含有微生物杀死的化合物. - 将含有 C. 线虫的培养皿放置在解剖范围内使用 .
- 使用铂金镐, 从培养皿中取出一条成熟的 C. 线虫 , 然后将其浸入试管中, 如有必要, 将其移入圆圈以去除线虫.
- 为了防止气泡在试管中形成, 请用水填充试管, 直到它在试管和 #39 上略微凸出。将试管和 #39 的盖子全部装满水, 然后迅速把盖子放到试管上.
- 使用光学清洁布除去可能溢出的水滴和光学清洁纸, 以清除任何小的残留液滴.
5。实时数据采集衍射图样强度变化
- 打开氦氖激光器并调整频率/颜色设置, 使其产生红色光束。打开传感器.
警告: 使用低能量光束在 632 nm, 作为 C. 线虫 避免高频 (蓝色) 光 12 。 - 在试管中定位蠕虫。将试管在其脊边上, 轻轻地倾斜试管直到线虫近似地位于试管部分的中心.
注意: 摇晃或倾斜试管猛烈地导致蠕虫与试管的墙壁碰撞。这会破坏线虫。- 将试管放置在光学系统中的支架上, 将其置于激光和 #39 中, 从镜像1反射到镜像2。一定要消除杂散光.
- 将蠕虫集中在激光束中。
- 将光电二极管置于衍射图样中, 以使光电二极管的位置和衍射图样的中心最大值不重合.
- 调整中性密度过滤器以防止光电二极管饱和。旋转中性密度过滤轮调节光强。
- 旋转中性密度滤镜, 使光电二极管的电压输出增加.
注: 使用数字光电二极管的软件观察电压输出。如果电压不变, 光电二极管就会饱和。在这种情况下, 旋转中性密度的过滤器, 直到电压降低, 而不会在最小读数时变平。确保电压信号在峰值读数不平坦, 指示光电二极管的饱和度。如果观测到饱和度, 则通过旋转中性密度过滤轮来减小光强.
- 旋转中性密度滤镜, 使光电二极管的电压输出增加.
- 一旦移动的衍射模式可见, 通过单击控制示波器的软件上的 "启动" 按钮来收集数据, 同时监控蠕虫和 #39 的移动。继续测量直到蠕虫从激光束和 diffrac 中移出模式消失, 通常需要大约二十年代。
- 通过单击示波器软件上的 "停止" 按钮来停止数据收集过程。以. csv 或 .txt 格式保存每个试用和 #39; s 数据.
- 重复步骤 5.2-5.3, 直到至少为每个表型收集了50数据集。每次试验使用八至十只动物.
- 如果试管被划伤处置它和蠕虫, 并重复步骤4。如果蠕虫在转移中受到损害, 则在重复步骤4使用相同的试管之前, 将其处理并用蒸馏水冲洗试管.
- 重复步骤5使用 OH7547 和 #34; 轧辊和 #34; 应变, 小心标签的数据, 以表明蠕虫的应变.
6。数据的傅立叶频谱
- 将获得的数据导入到能够执行离散傅立叶变换的数据分析程序 3 。
- 使用软件的快速傅立叶变换 (FFT) 选项对每个数据集执行傅立叶变换.
- N2 野生类型蠕虫的每个振幅的平均频率.
- 使用 OH7547 和 #34 中的 FFTs 重复步骤 6.3; 滚轮和 #34; 蠕虫.
7。傅立叶谱的建模
- 对线虫运动的二进制模型进行编程 (参见补充材料中包含的程序)
注意: 该模型是一个粗略的近似, 它首先显示了蜗杆运动的显著特征。模型可以被提炼, 因为结果与实际蠕虫比较。蠕虫形状是使用显微镜图像 13 近似的。- 创建二进制模型的连续帧, 以移动至少两个蠕虫周期 ( 图 2a )。查看补充材料中的视频;C 蠕虫视频 (CWorm) 和 W 蠕虫视频 (WWorm. avi).
- 生成连续的衍射图案。请参阅补充材料中的视频。C 蠕虫衍射视频 (CWormDiff. avi) 和 W 蠕虫视频 (WWormDiff. avi)。
- 傅立叶变换蠕虫图像的每个二进制帧。围绕蠕虫的帧的填充越大, 衍射图像的分辨率就越高.
注: 每个傅里叶变换帧的绝对值与相应的衍射图样成正比 ( 图 2b ). - 通过将衍射图样的强度映射到对数刻度来调整衍射图案的对比度.
注意: 相机和眼睛的功能在一个非线性的规模。对数刻度可以模拟人们眼中的衍射图样.
- 傅立叶变换蠕虫图像的每个二进制帧。围绕蠕虫的帧的填充越大, 衍射图像的分辨率就越高.
- 提取模型化的衍射信号。
- 选取与光电二极管在衍射图样中的位置对应的偏心位置.
- 在光电二极管的位置周围添加相邻的矩阵元素以模拟光电二极管的大小。光电二极管的大小通常是0.1% 的模拟衍射模式.
- 记录和绘制衍射信号的大小序列。检查信号的物理合理性 ( 即 , 在周期性颠簸的情况下, 光电二极管发出的信号也应该是周期性的).
- 傅立叶变换在7.3 获得的衍射信号, 并与实验数据进行比较.
- 重复不同的蠕虫病毒并进行比较.
Representative Results
在图 1中显示的光学实验装置允许对微生物进行研究, 而不被绑在一个焦平面上。在采集数据的同时, 可以在计算机屏幕上实时观察到光电二极管的抖动信号。不寻常的模式将立即可见, 而不必详细分析视频。
模型化顺序蠕虫移动和相应的衍射模式的示例如图 2所示。模型化的衍射模式定性地类似于实验模式1 , 是模拟成功地对线虫进行建模的初步迹象。
本文所研究的两种类型的C. 线虫的时间衍射样本在图 3中显示。从定性上看, 每种线虫的以率和振幅不同。某些差异可以通过曲线拟合进行量化, 如以前的出版物1所做的那样。然而, 离散傅立叶变换揭示了有关嵌入频率的更多细节:
, (1)
其中, Fk 是数字傅立叶变换 (FT) 和Fn 是与离散时间变量n和离散频率变量 n 的原始衍射信号, 其总数为数据点。平均数字傅立叶变换允许线虫被确定由它的衍射频谱的振幅 (图 4)。野生型谱的控制频率比滚筒运动谱低。
接近野生类型的模型与滚子C. 线虫注意, 野生类型倾向于在波 (W 或 S 形状) 运动中进行捶打 (图 2a), 而滚轮倾向于偏向于一个大致类似于振荡 C 形状的一侧 (图 5)。这为不同的光谱提供了一些解释。滚子将主要形成一个 c 到一侧, 而 W 振荡可以被认为是两个相反的 c 运动。因此, W 运动比 C 运动更复杂地揭示出更多的次低频频率。这一结果在计算模型中得到证实。W 形状的频率密度比 C 形状高得多 (图 6)。这在图 4的 FFT 中得到确认, 在这里, 滚子频率更集中, 而不是完全离散。轧辊的统计数据是倾斜的, 因为滚轮可以暂时恢复为野生类型的运动。
滚子类型的平滑功率谱c. 线虫显示了一个宽峰值在 ~ 1.5 赫兹, 而游泳狂放的类型c. 线虫陈列多式联运光谱 (包括峰顶在 ~ 1.0 hz 和 1.75 hz)。光电二极管 (PD) 有一个有限的大小扩散在几个矩阵元素。由于构造性和破坏性干扰不同, 衍射图样上的单个矩阵元素或点的强度各不相同;然而, 强度变化的频率对于所有矩阵元素都是相同的, 如图 7所示。考虑到时间导数 Eq 1, 可以看出频率波动不依赖于相位矩阵, 而只取决于原始对象的波动:
, (2)
当 PD 在多个矩阵元素上传播时, 峰值位置平均到一致的频率剖面。一些变化可以预期, 可以提供线索的方向蠕虫。当蠕虫的步态改变时, 频率分布会发生变化。目前的模型是一个简单的模型, 只允许评估峰值位置, 而不是相对峰值高度。不同的运动模式将平均到不同的峰值位置。
图1。实验装置.低功率激光束通过中性密度滤光片, 通过反射镜 M1 向下穿过试管, 将蜗杆载入镜面 M2, 并向光电二极管方向移动。请单击此处查看此图的较大版本.
图2。顺序蠕虫形状和相应的衍射图案.(a) 一些选定的 W 型线虫的序列二进制图像和 (b) 相应的连续衍射模式。请单击此处查看此图的较大版本.
图3。实验样品衍射特征.为 (a) OH7547 "滚子" 和 (b) N2 的衍射特征码在衍射图样中使用单光电二极管 .请单击此处查看此图的较大版本.
图4。实验平均功率谱的轧辊和野生型霍夫琅衍射系列.该谱显示了在光电二极管记录的时间序列的平均傅里叶变换中存在的频率。标准偏差的高斯滤波器0.075 赫兹, 被截断在3标准偏差, 用于平滑。请注意平滑滚子频谱中的宽光谱峰值 (1.5 赫兹), 与多式平滑的野生类型频谱 (包括在 ~ 1.0 和1处的峰值) 比较75赫兹)。请单击此处查看此图的较大版本.
图5。衍射形成例证.衍射模式可以建模的每一个线段的思维无穷小直线 (左)。叠加这些线 (右) 演示了由 C 形线虫产生的远场衍射图样的构造。请单击此处查看此图的较大版本.
图6。滚子和野型霍夫弗衍射系列的模拟功率谱.(a) C 形状和 (b) W 形状的蠕虫与光电二极管中心在矩阵元素 200 (垂直) 和 175 (水平)。W 形状显示更高的频率密度由于更复杂的运动。请单击此处查看此图的较大版本.
图7。不同光电二极管位置的滚子和野型霍夫弗衍射级数的模拟功率谱.(a) W 形状蠕虫和 (b) C 形状的蠕虫, 用于在不同位置模拟光电二极管的不同位置的单个矩阵元素。不同地点的高峰高度各不相同;但是, 特定形状的峰值位置仍然保持不变。请单击此处查看此图的较大版本.
Discussion
包括不活动的数据的延伸将扭曲结果, 因为人工较低的频率将被平均纳入结果。在原始数据中, 可以通过平坦峰值或 "切断" 峰值来识别光电二极管的饱和。按峰值强度划分每个原始数据将有助于计算激光强度的波动。
峰值频率是总颠簸频率的指示器;然而, 复杂的运动会引起衍射图样中跳动频率的干扰, 需要仔细检查。
该方法可用于研究其他线虫的运动。可以将环境更改为其他媒体。波长也可能改变。工作在可见光范围内的电磁波谱是最简单和最安全的。
更精确的模型将在未来更真实地模拟衍射光谱。未来的模型可能包括可以改变方向的蠕虫, 这不会影响频率位置, 而是相对峰值高度。一个更现实的模型将允许概率分布的颠簸频率, 这将扩大峰值的实验数据。频率的传播会考虑到颠簸频率的变化。
当前的蠕虫形状是粗的, 特别是在头部和尾部区域, 这应该比目前的模型更锥形。对信号的时间序列进行详细的分析可能会很有趣, 因为它可以给出不同突变体的运动复杂性的线索。
将该技术推广到同时对多线虫进行定性是值得考虑的。这种方法应该被理解为使用传统显微镜的现有方法的一种补充方法。这种方法在数据获取过程中不需要显微镜, 因此蠕虫可能会移出焦平面。平均频率谱在蜗杆运动中显示出明显的差异, 并可通过流行的频率峰值进行量化, 这是一种新的定量的蠕虫运动方法。衍射特征的数据分析正在进一步发展, 并有望导致多个突变体和个体的自动识别过程。
Disclosures
作者没有什么可透露的。
Acknowledgments
我们感谢胡安。我们感谢瓦萨大学本科生研究暑期研究所 (熊胆)、露西-梅纳德三文鱼研究基金和 NSF 奖 No. 1058385 的支持。
Materials
| Name | Company | Catalog Number | Comments |
| Tunable Helium-Neon laser | Research Electro-Optics | 30602 | Four wavelengths can be selected between 543 nm and 633 nm. |
| 2 Front Surface Aluminum Mirrors | Thorlabs | PF10-03-F01 | |
| Photodiode: SI Amplified Detector | Thorlabs | PDA 100A | |
| Quartz Cuvette | Starna Cells | 21/G/5 | Plastic cells may be used as well. |
| MatLab (Software) | MathWorks | R2016b (9.1.0.441655) | Use the fft command to simulate diffraction |
| Excel | Microsoft | 14.7.1 | Used for data analysis of Figure 4 |
| Caenorhabditis elegans Roller | University of Minnesota Caenorhabditis elegans Center (CGC) | Strain: OH7547 Genotype: otIs199. |
https://cbs.umn.edu/cgc/home |
| Caenorhabditis elegans Wild Type | University of Minnesota Caenorhabditis elegans Center (CGC) | Strain:N2 Genotype: C. elegans wild isolate | https://cbs.umn.edu/cgc/home |
References
- Magnes, J., et al. Analysis of Freely Swimming C. elegans Using Laser Diffraction. Open J. Biophys. 2, 101-107 (2012).
- Magnes, J., Raley-Susman, K. M., Eells, R. Quantitative Locomotion Study of Freely Swimming Micro-organisms Using Laser Diffraction. J. Vis. Exp. (68), (2012).
- James, J. F. A Student's Guide to Fourier Transforms with Applications in Physics and Engineering. , Cambridge University Press. Cambridge. (1995).
- Korta, J., Clark, D. A., Gabel, C. V., Mahadevan, L., Samuel, A. D. T. Mechanosensation and mechanical load modulate the locomotory gait of swimming C. elegans. J. Exp. Biol. 210, (2007).
- Martin-Garcia, J. M., Conrad, C. E., Coe, J., Roy-Chowdhury, S., Fromme, P. Serial femtosecond crystallography: A revolution in structural biology. Arch. Biochem. Biophys. 602, 32-47 (2016).
- Thibault, P., Rankenburg, I. C. Optical diffraction microscopy in a teaching laboratory. Amer. J. Phys. 75 (9), 827-832 (2007).
- Miao, J., Ishikawa, T., Anderson, E. H., Hodgson, K. O. Phase retrieval of diffraction patterns from non crystalline samples using the oversampling method. Phys. Rev. B. 67, 174104 (2003).
- Zhang, Y. P., Zhang, J. Q., Xu, W. Method for eliminating zero-order diffraction in lensless Fourier transform digital holography. Optik - International Journal for Light and Electron Optics. 124 (21), 4873-4875 (2013).
- Brody, A. H., Chou, E., Gray, J. M., Pokrywka, N. J., Raley-Susman, K. M. Mancozeb-induced behavioral deficits precede structural neural degeneration. NeuroToxicology. 34, 74-81 (2013).
- Stiernagle, T. The C. elegans Research Community. Maintenance of C. elegans. WormBook. , (2006).
- Dasalukunte, D., Öwall, V., Rusek, F., Anderson, J. B. Faster than Nyquist Signaling. Algorithms to Silicon. , Springer International Publishing. Switzerland. (2014).
- Edwards, S. L., et al. A novel molecular solution for ultraviolet light detection in Caenorhabditis elegans. PLoS Biol. 6 (8), e198 (2008).
- Bilbao, A., Wajnryb, E., Vanapalli, S. A., Blawzdziewicz, J. Nematode Locomotion in Confined and Unconfined Fluids. Phys. Fluids. 25, 081902 (2013).
