材料和试剂

12 mm圆形盖玻片（Fisher Scientific，目录号：11846933）
6孔板（Corning，目录号：353046）
16% 甲醛（Thermo Fisher Scientific，目录号：28908）
NaCl（VWR International，目录号：10241）
KCl（VWR International，目录号：10198）
KH 2 PO 4 （VWR International，目录号：10203）
Na 2 HPO 4 （VWR International，目录号：10249）
HCl（VWR International，目录号：45002）
Triton-X-100（Sigma-Aldrich，目录号：T8787）
牛血清白蛋白（BSA）（Sigma-Aldrich，目录号：A7906）
4',6-二脒-2'-苯基吲哚二盐酸盐（DAPI）（Thermo Fisher Scientific，目录号：D1306）
甲醇（VWR International，目录号：20847.318）
延长金镶样介质（Thermo Fisher Scientific 目录号：P10144）
Superfrost粘附载玻片（Thermo Fisher Scientific，目录号：10149870）
聚-L-赖氨酸（Sigma-Aldrich，目录号：P6282）
PTEM 缓冲液（见配方）
1×磷酸盐缓冲盐水（1×PBS）（见配方）
1× PBST（见食谱）
1× PBSBT（见食谱）

抗体

一抗：小鼠抗乙酰化微管蛋白（1:2,000 稀释；Sigma-Aldrich，目录号：T7451）
小鼠抗γ-微管蛋白（1:500 稀释；Santa Cruz，目录号：sc-51715）
兔抗 IFT81 抗体（1:200 稀释；Proteintech，目录号：17744-1-AP）
兔抗 Arl13B（1:500；Proteintech，目录号：17711-1-AP）
二抗：Alexa Fluor 488（Thermo Fisher Scientific，目录号：A-21141）、594（Thermo Fisher Scientific，目录号：A11037）或 680（Thermo Fisher Scientific，目录号：A31562）（1:500 稀释）

细胞系

如所示使用A204和RPE-1细胞系。细胞分别在含有 10% FBS 的 Dulbecco's Modified Eagle 培养基 (DMEM) 和 Dulbecco's Modified Eagle 培养基：营养混合物 F-12 (DMEM:F12) 培养基中生长。

设备

轨道振动筛
落射荧光显微镜：正置显微镜（Axio Imager M2；Zeiss），配备 100 倍油镜、1.4 数值孔径 (NA)、相机（ORCA R2；Hamamatsu Photonics）和带有图像处理软件的计算机（Zen）
蔡司 710 共聚焦显微镜与 Zen 软件
计算机规格：运行 MATLAB – Intel 或 AMD x86-64 处理器，具有四个逻辑内核和 AVX2、8 GB RAM、Windows 7 或更新的操作系统

软件

Mathworks MATLAB（R2020b 版）和图像处理工具箱
Excel（微软）
Zen（蔡司图像采集软件）

程序

用于纤毛发生测定的种子细胞和固定

在 140 × g的定轨振荡器中的烧瓶中，在 55°C 下用 1 N HCl洗涤13 mm 盖玻片。
大量洗涤后，在室温下用 0.1 毫克/毫升聚-L-赖氨酸覆盖盖玻片 1 小时。
用水清洗盖玻片并晾干过夜。
电镀前，将盖玻片置于紫外光下 15 分钟。
在 6 孔板（直径 = 3.5 厘米）中，以每孔4 × 10 5的密度将种子细胞接种到盖玻片上。这大约是 60% 的密度，因此第二天细胞将在 80% 左右。这可能因细胞系而异，我们建议进行优化。
允许细胞附着 24 小时。
用无血清培养基 (SFM) 清洗细胞 3 次，并在 SFM 中保持 48 小时。
用 1× PBS（或者用 PTEM 缓冲液）洗涤 2 次，然后进行固定（根据需要，在 -20°C 下用 1× PBS 或甲醇稀释 4% 的甲醛）。

免疫染色方案

注意：所有步骤均在室温下进行。

将盖玻片放在潮湿的室内，用 1× PBST 清洗 2 次以渗透细胞。留在第二次洗涤 5 分钟。
使用 3% (w/v) 牛血清白蛋白在 1× PBST (PBSBT) 中封闭细胞 5 分钟。
添加在 PBSBT 中稀释的主要抗体并孵育 1 小时。
用 1× PBST 清洗盖玻片 3 次，每次 5 分钟。
添加在 PBSBT 中稀释 1:500 的二级抗体并孵育 1 小时。
用 1× PBST 清洗盖玻片 3 次，每次 5 分钟。
添加在 PBSBT 中稀释的 DAPI 并孵育 2 分钟。
清洗 2 次，将 ProLong Gold 滴到载玻片上，然后安装。允许在黑暗中干燥过夜。
24 小时后，用指甲油密封。

图像采集

使用 100 倍油镜 (1.4 NA) 获取图像。在不同条件下为所有图像设置相同的曝光时间。使用蔡司分析软件（Zen）获取图像。

MATLAB 脚本说明和分析

CILIA 1（四通道原始强度）：cilia_linesc_raw_int_v9

该脚本沿手动绘制的线扫描测量荧光强度（通道 2），并沿线扫描输出 10 个相等段的平均强度。这允许用户比较不同长度的纤毛。纤毛使用来自 MATLAB 的 improfile 函数以用户交互方式进行分割，该函数沿用户定义的多线路径检索像素的强度值。该脚本还测量分段纤毛的长度。数据导出为 Excel 电子表格并使用与 TIFF 文件相同的文件名保存。纤毛标记物（即乙酰化微管蛋白或 Arl13B）可用于定义纤毛路径。然后从感兴趣的通道中检索强度分布。为了减少噪声，沿纤毛的每个位置测量三个像素（路径上方、上方和下方）的平均荧光强度（取多个像素的平均值可减少单个噪声像素的影响）。下面详细介绍了如何使用此脚本。

四 (4) 个荧光通道。对于所示示例：通道 1 (C1) 是中心粒标记 (γ-微管蛋白)，通道 2 (C2) 是感兴趣的蛋白质 (IFT81)，通道 3 (C3) 是纤毛标记（乙酰化微管蛋白），通道4 (C4) 是核染色 (DAPI)。纤毛是手动定义的，然后脚本会将纤毛分成 10 个相等的部分（bins）。通道 3 的原始像素强度将沿着定义的纤毛进行测量，并显示为从基部到纤毛尖端的 10 个纤毛部分。在下面的示例中，此脚本用于量化 IFT81（图 1-8）。

图 1. 用于 CILIA 1 脚本 MATLAB 分析的 A204 IFT81 图像。

A204 细胞血清饥饿 48 小时以诱导纤毛发生，然后固定并用 IFT81（红色）、乙酰化微管蛋白（绿色）和γ-微管蛋白（白色）的抗体染色。核染色是 DAPI（蓝色）。比例尺 = 10 μm。

打开MATLAB，导入数据，准备脚本

图像导出和保存

将图像导出为 TIFF 文件格式；每个荧光通道都需要单独的图像。在所示示例中，用于分析的通道已按以下顺序导出：中心粒标记、通道 1 (C1)、要分析的信号、通道 2 (C2)、乙酰化微管蛋白或所选纤毛标记、通道 3 (C3) 和DAPI 或替代核染色，通道 4 (C4)。所有通道 C1 + C2 + C3 + C4 的叠加图像包含在导出的文件中。为了使脚本正常工作，通道 2 (C2) 必须是要分析的通道，通道 3 (C3) 必须是纤毛标记。如果需要，可以更改用于通道 1 和 4 的标记。

打开 MATLAB。
通过拖动文件将四个脚本文件添加到 MATLAB 目录（左侧栏）。
将要分析的 TIFF 图像文件添加到左侧边栏（图 2）。

图 2.将数据导入 MATLAB。

将需要分析的脚本文件和图片拖入MATLAB操作界面左侧边栏。图像必须作为单独的荧光通道作为 TIFF 文件导入。

在脚本的第 12 行，将文件名更改为与要分析的文件相同。从脚本中的文件名中排除 TIFF 文件扩展名（图 3）。

图 3. 更改脚本文件名。

将脚本中的文件名（第 12 行，以蓝色突出显示）更改为与要分析的文件相同的名称。文件名中不应包含 TIFF 扩展名。

单击脚本中第 28 行旁边的破折号 (-)；这将在脚本中添加一个断点，允许脚本在此时暂停（图 4）。

图 4. 选择要分析的文件并在脚本中插入断点。

通过单击第 28 行旁边的破折号在脚本中插入一个断点；这会将破折号转换为红色圆圈（用箭头突出显示）。通过单击“运行”（红色框内的绿色图标）运行脚本。

定义纤毛并分析 IFT81 像素强度

按图 4 中突出显示的“run”顶部的绿色图标运行脚本。乙酰化微管蛋白的图像将出现在新窗口中。单击“工具>放大”并在要分析的纤毛周围绘制一个正方形。
按脚本上的“继续”（“继续”与运行图标位于同一位置）。
使用光标，沿纤毛中心画一条线，从根部到纤毛尖端。如果纤毛弯曲，请停止、单击并继续。当线到达纤毛尖端时，用鼠标左键双击（图 5）。

图 5. 确定要分析的纤毛。

沿着放大的纤毛的中心画一条线，确保线从根部到纤毛的尖端。当线位于纤毛尖端时双击。本示例使用 A204 单元格。

沿着纤毛绘制的线将自动叠加到 IFT81 染色图像上，并沿着这条线测量像素强度（测量在整个纤毛的尖端和横向延伸几个像素，以确保纤毛的像素强度是测量）（图 6）。

图 6. 覆盖在 IFT81 图像上的纤毛线。

沿着纤毛画的线会自动转移到IFT81图像上，测量这部分图像的像素强度。

显示了 IFT81 染色的直方图（图 7）。

图 7. IFT81 像素强度直方图。

使用脚本生成的直方图示例。显示了纤毛（1-10，从近端到远端）的每个部分的平均像素强度。y 轴是平均像素强度，x 轴显示纤毛分成 10 个相等的部分。第 1 部分代表纤毛基部，第 10 部分代表尖端。使用中心体标记物γ-微管蛋白确定纤毛的基部。

当前文件夹侧栏中显示了纤毛 IFT81 染色的 Excel 电子表格。该电子表格将纤毛的平均像素强度分为 10 个部分（1 个近端-10 个远端），还包括纤毛长度。如果要分析每个图像的多个纤毛，更改文件的名称很重要，以便在分析后续纤毛时不会覆盖文件（图 8）。

图 8. 包含数据的 Excel 电子表格。

IFT81 像素强度的Excel 电子表格将出现在当前文件夹的侧栏中（以蓝色突出显示）。如果每个图像要分析多个纤毛，请更改文件名以避免覆盖此文件。

数据分析

在典型的实验中，每个实验条件使用 150 根纤毛分析像素强度。然后将数据呈现为沿纤毛近端到远端部分的像素强度。使用Student's t-检验计算统计显着性。在这个例子中，我们的分析表明，达沙替尼抗性 A204 细胞沿着纤毛的近端和远端具有较低的 IFT81 蛋白水平（Jenks等，2018）。该方法可用于分析由可使用免疫荧光成像的蛋白质或生物分子给出的任何信号的分布。这突出了该分析工具的多功能性，它能够准确量化纤毛定位蛋白质的分布。

这个详细的协议将允许用户：

分析可以使用免疫荧光观察到的任何蛋白质的纤毛分布。
比较不同细胞系组（例如，耐药细胞系与药物敏感细胞系）中纤毛蛋白的分布。
定量分析蛋白质定位到纤毛的水平。
测量纤毛长度。
提供发布数据。

其他 MATLAB 脚本：CILIA 2、CILIA 3、CILIA 4 和 CILIA 5

四个额外的 MATLAB 脚本被设计用于分析纤毛蛋白定位。每个脚本都使用与所描述示例类似的工作流程（图 1-8）。CILIA 2 可以自动识别和量化纤毛长度和强度。CILIA 2 和 CILIA 3 要求组合图像具有 4 个荧光通道。为了使 CILIA 2 (cilia_linesc_aut_v6) 从图像中自动量化纤毛，导出的 TIFF 荧光通道的顺序必须如下：通道 1-中心标记、通道 2-要分析的信号、通道 3-Arl13B 或纤毛标记选择，和通道 4-DAPI 或替代核标记。必须包含所有通道的叠加层，以便此脚本自动识别和定向纤毛。CILIA 3 还要求通道 2 是要分析的信号，通道 3 是纤毛标记；然而，其余通道可以更改，因为纤毛是手动定义的。CILIA 2 (cilia_linesc_aut_v6)（图 9-15 中详述的协议）和 CILIA 3（cilia_linesc_int_v9）（图 9、16-22 中详述的协议）将染色强度表示为 10 个纤毛部分中每一个的平均像素强度除以平均像素整个纤毛的强度。其余两个额外的脚本，CILIA 4 (cilia_linesc_raw_int_3c_v9) 和 CILIA 5 (cilia_linesc_int_3c_v9)，只需要三个荧光通道，并且需要按如下方式组织导出的 TIFF 荧光通道：通道 1-Arl13B 或选择的纤毛标记和通道 2-信号到进行分析。对于所示示例，通道 3 中使用 DAPI 来定义细胞核（如果需要，可以更改此染色）。CILIA 4 (cilia_linesc_raw_int_3c_v9) 和 CILIA 5 (cilia_linesc_int_3c_v9) 仅在染色强度的呈现方式上有所不同。CILIA 4 (cilia_linesc_raw_int_3c_v9)（图 23-30 中详述的协议）产生的结果是平均像素强度，而 CILIA 5 (cilia_linesc_int_3c_v9)（图 31-37 中详述的协议）产生的结果是被划分的 10 个纤毛部分中的每一个的平均像素强度通过整个纤毛的平均像素强度。这些额外的脚本在下面进一步描述，使用乙酰化微管蛋白和 Arl13B 染色的纤毛复合图像（图 9）。

图 9. CILIA 2 和 CILIA 3 脚本分析中使用的 RPE-1 纤毛图像。

将 RPE-1 细胞血清饥饿 48 小时以诱导纤毛发生，然后固定并用 Arl13B（红色）、乙酰化微管蛋白（绿色）和γ-微管蛋白（白色）的抗体染色。核染色是 DAPI（蓝色）。比例尺 = 10 μm。

纤毛2 —— 自动化四通道自动化脚本（cilia_linesc_aut_v6）：

四 (4) 个荧光通道。对于所示示例，C1 是中心粒标记（γ-微管蛋白），C2 是感兴趣的蛋白质（乙酰化微管蛋白），C3 是纤毛标记（Arl13B），C4 是核染色（DAPI）。纤毛由脚本自动定义并分成 10 个相等的部分。数据表示为每个纤毛部分的通道 2 像素强度除以沿整个纤毛的平均通道 2 像素强度。该比率将显示为 10 个纤毛部分（垃圾箱），从底部到纤毛尖端。在图 10-15 所示的示例中，它用于量化乙酰化微管蛋白。要检测所有纤毛，可能需要在脚本的第 10 行（级别）上调整通道 2 像素强度阈值。

图 10. 将数据导入 MATLAB。

将要分析的脚本文件和图像拖入操作界面最左侧的MATLAB窗口中。图像必须作为单独的荧光通道作为 TIFF 文件导入。

图 11. 更改脚本文件名。

将脚本中的文件名（第 13 行，以蓝色突出显示）更改为与要分析的文件相同的名称。文件名中不应包含 TIFF 扩展名。

图 12.运行脚本。

按“运行”（顶部突出显示的绿色图标）启动脚本。

图 13. 覆盖在乙酰化微管蛋白图像上的纤毛线。

纤毛将使用 Arl13B 染色和γ- tub自动定义以定位纤毛方向。这些纤毛将被转移到乙酰化微管蛋白图像上，并且将分析图像这些部分的像素密度。

图 14. 乙酰化微管蛋白像素强度直方图。

从包含六个纤毛的图像生成的直方图示例。y 轴将像素强度表示为平均像素强度除以整个纤毛的平均像素强度的比率。第 1 节代表纤毛基部，第 10 节代表纤毛尖端。使用中心体标记物γ-微管蛋白确定纤毛的基部。

图 15. 像素比 Excel 电子表格示例。

乙酰化微管蛋白像素比的 Excel 电子表格将出现在当前文件夹的侧栏中（以蓝色突出显示）。

CILIA 3 – 四通道像素比 (cilia_linesc_int_v9)：

四 (4) 个荧光通道。对于所示示例，C1 是中心粒标记 ( γ- tub)，C2 是感兴趣的蛋白质（乙酰化微管蛋白），C3 是纤毛标记 (Arl13B)，C4 是核染色 (DAPI)。手动定义纤毛后，脚本会将纤毛分成 10 个相等的部分。每个部分的通道 2 像素强度将除以整个纤毛的平均通道 2 像素强度。该比率将显示为从基部到纤毛尖端的 10 个纤毛部分。在图 16-22 所示的示例中，它用于量化乙酰化微管蛋白。

图 16. 将数据导入 MATLAB。

将需要分析的脚本文件和图片拖入操作界面左侧的MATLAB目录下。图像必须作为单独的荧光通道作为 TIFF 文件导入。

图 17. 更改脚本文件名。

将脚本中的文件名（第 12 行，以蓝色突出显示）更改为与要分析的文件相同的名称。文件名中不应包含 TIFF 扩展名。

图 18. 选择要分析的文件并在脚本中插入断点。

通过单击第 28 行旁边的破折号在脚本中插入一个断点。这会将破折号转换为红色圆圈（用箭头突出显示）。按“运行”（顶部突出显示的绿色图标）启动脚本。

图 19. 识别要分析的纤毛。

沿着放大纤毛的中心画一条线，确保这条线从根部到纤毛的尖端。当线位于纤毛尖端时双击。

图 20. 纤毛线叠加在乙酰化微管蛋白图像上。

沿着纤毛绘制的线将转移到乙酰化微管蛋白图像。然后，将分析这部分图像的像素强度。

图 21. 乙酰化微管蛋白像素强度直方图。

为纤毛每个部分的平均像素强度生成的直方图示例（1 个近端到 10 个远端）。y 轴将像素强度显示为每个部分的平均像素强度除以整个纤毛的平均像素强度的比率。第 1 节代表纤毛基部，第 10 节代表纤毛尖端。使用中心体标记物γ-微管蛋白确定纤毛的基部。

图 22. 平均像素强度 excel 电子表格示例。

乙酰化微管蛋白像素强度的 Excel 电子表格将出现在当前文件夹的侧栏中（以蓝色突出显示）。如果每个图像要分析多个纤毛，请更改文件名以避免覆盖此文件。

最后，对于 CILIA 4 和 CILIA 5 脚本的示例，我们使用了图 9 中所示的图像，没有γ-微管蛋白染色（图 23）。

图 23. CILIA 4 和 CILIA 5 脚本分析中使用的 RPE-1 纤毛图像。

将 RPE-1 细胞血清饥饿 48 小时以诱导纤毛发生，然后固定并用 Arl13B（红色）和乙酰化微管蛋白（绿色）的抗体染色。核染色是 DAPI（蓝色）。比例尺 = 10 μm。

CILIA 4 – 三通道原始强度脚本 (cilia_linesc_raw_int_3c_v9)：

该脚本使用三 (3) 个荧光通道。对于所示示例，C1 是纤毛标记 (Arl13B)，C2 是感兴趣的蛋白质（乙酰化微管蛋白），C3 是核染色剂 (DAPI)。手动定义纤毛后，脚本会将纤毛分成 10 个相等的部分，通道 2 的原始像素强度将沿着定义的纤毛进行测量，并显示为从底部到纤毛尖端的 10 个纤毛部分。在图 24-30 所示的示例中，这用于量化乙酰化微管蛋白。

图 24. 将数据导入 MATLAB。

将要分析的脚本文件和图像拖入操作界面最左侧的MATLAB窗口中。图像必须作为单个荧光通道的 TIFF 文件导入。

图 25. 更改脚本文件名。

将脚本中的文件名（第 12 行，以蓝色突出显示）更改为与要分析的文件相同。文件名中不应包含 TIFF 扩展名。

图 26. 选择要分析的文件并在脚本中插入断点。

通过单击第 28 行旁边的破折号在脚本中插入一个断点，这会将破折号转换为红色圆圈（使用箭头突出显示）。按“运行”启动脚本（顶部的绿色图标“突出显示”）。

图 27. 识别要分析的纤毛。

沿着放大纤毛的中心画一条线，从根部到纤毛的尖端。当线位于纤毛尖端时双击。

图 28. 覆盖在乙酰化微管蛋白图像上的纤毛线。

沿着纤毛绘制的线将转移到乙酰化微管蛋白图像，并分析图像中该部分的像素强度。

图 29. 乙酰化微管蛋白像素强度直方图。

为纤毛每个部分（从近端 1 到远端 10）中的平均像素强度生成的直方图示例。y 轴显示像素强度，x 轴显示纤毛分成 10 个相等的部分。第 1 节代表纤毛基部，第 10 节代表纤毛尖端。

图 30. 平均像素强度 excel 电子表格示例。

乙酰化微管蛋白像素强度的 Excel 电子表格将出现在当前文件夹的侧栏中（以蓝色突出显示）。如果每个图像要分析多个纤毛，请更改文件名以避免覆盖。

CILIA 5 – 三通道像素比 (cilia_linesc_int_3c_v9)：

该脚本使用三 (3) 个荧光通道。对于所示示例，C1 是纤毛标记 (Arl13B)，C2 是感兴趣的蛋白质（乙酰化微管蛋白），C3 是核染色 (DAPI)。手动定义纤毛后，脚本会将纤毛分成 10 个相等的部分，每个部分的通道 2 像素强度将除以整个纤毛的平均通道 2 像素强度。该比率将显示为从基部到纤毛尖端的 10 个纤毛部分。在图 31-37 所示的示例中，这用于量化乙酰化微管蛋白。

图 31. 将数据导入 MATLAB。

将需要分析的脚本文件和图像拖入操作界面的MATLAB左侧边栏窗口中。图像必须从各个荧光通道导入为 TIFF 文件。

图 32. 更改脚本文件名。

将脚本中的文件名（第 12 行，以蓝色突出显示）更改为与要分析的文件相同的名称。文件名中不应包含 TIFF 扩展名。

图 33. 选择要分析的文件并在脚本中插入断点。

通过单击第 28 行旁边的破折号在脚本中插入一个断点，将破折号转换为红色圆圈（使用箭头突出显示）。按“运行”启动脚本（顶部的绿色图标“突出显示”）。

图 34. 识别要分析的纤毛。

沿着放大纤毛的中心画一条线，确保这条线从根部到纤毛的尖端。当线位于纤毛尖端时双击。

图 35. 覆盖在乙酰化微管蛋白图像上的纤毛线。

沿着纤毛绘制的线将转移到乙酰化微管蛋白图像，并分析图像这部分的像素强度。

图 36. 乙酰化微管蛋白像素强度直方图。

使用 CILIA 5 脚本计算的直方图示例。y 轴将像素强度显示为每个部分的平均像素强度除以整个纤毛的平均像素强度的比率。第 1 节代表纤毛基部，第 10 节代表纤毛尖端。

图 37. 像素强度 excel 电子表格示例。

乙酰化微管蛋白像素强度的 Excel 电子表格将出现在当前文件夹的侧栏中（以蓝色突出显示）。如果每个图像要分析多个纤毛，请更改文件名以避免覆盖。

食谱

PTEM 缓冲液

20 mM 管道，pH 6.8

0.2% 海卫 X-100

10 mM EGTA

1 毫米氯化镁2

1 ×磷酸盐缓冲盐水 (1 × PBS)

溶解 10 g/L NaCl、0.25 g/L KCl、0.25 g/L KH 2 PO 4和 1.437 g/L Na 2 HPO 4 。

使用 HCl 将 pH 值调节至 7.2。

1×PBST

1×PBS

0.1% v/v Triton-X-100

1×PBSBT

1×PBS

0.1% (v/v) Triton-X-100

3% (w/v) 牛血清白蛋白

0.05% (w/v) NaN 3

致谢

这项研究部分由癌症研究所和英国肉瘤基金（给 BET）资助。特别感谢 Adam Tyson (UCL) 对这份手稿的批判性阅读。

利益争夺

不存在利益冲突或竞争利益。

参考

Anvarian, Z.、Mykytyn, K.、Mukhopadhyay, S.、Pedersen, LB 和 Christensen, ST（2019 年）。初级纤毛在发育、器官功能和疾病中的细胞信号传导。Nat Rev Nephrol 15（4）：199-219。
Christensen, ST, Clement, CA, Satir, P. 和 Pedersen, LB (2012)。初级纤毛和受体酪氨酸激酶 (RTK) 信号的协调。J Pathol 226(2): 172-184。
Goetz, SC 和 Anderson, KV (2010)。初级纤毛：脊椎动物发育过程中的信号中心。Nat Rev Genet 11(5): 331-344。
Hansen, JN、Rassmann, S.、Stuven, B.、Jurisch-Yaksi, N. 和 Wachten, D. (2021)。CiliaQ：一个简单的开源软件，用于自动量化 2D、3D 和 4D 图像中的纤毛形态和荧光。Eur Phys JE Soft Matter 44(2): 18。
Jenks, AD, Vyse, S., Wong, JP, Kostaras, E., Keller, D., Burgoyne, T., Shoemark, A., Tsalikis, A., de la Roche, M., Michaelis, M., Cinatl, J., Jr.、Huang, PH 和 Tanos, BE (2018)。原发性纤毛介导癌症中的多种激酶抑制剂抗性机制。细胞代表23(10)：3042-3055。
Lauring, MC, Zhu, T., Luo, W., Wu, W., Yu, F. 和 Toomre, D. (2019)。用于细胞中自动纤毛检测的新软件 (ACDC)。纤毛8：1。
Schou, KB, Pedersen, LB 和 Christensen, ST (2015)。初级纤毛中 GPCR 信号的来龙去脉。EMBO 代表16(9)：1099-1113。