Summary
我们已经开发了一种用于研究异源细胞群体动力学以响应扰动的新方案。该手稿描述了基于成像的平台，其产生定量数据集，以稳健的方式同时表征异源细胞群体的多种细胞表型。
Abstract
细胞过程是复杂的，并且是由多种细胞类型与其环境之间的相互作用产生的。现有的细胞生物学技术往往不能准确地解释这种相互作用。使用基于定量成像的方法，我们提出了一个高含量的方案来表征异质细胞群体对环境刺激变化的动态表型反应（ 即形态学变化，增殖，凋亡）。我们强调我们根据应用来区分基于荧光强度或固有形态特征的细胞类型的能力。该平台允许更全面地表征亚群对扰动的反应，同时利用更短的时间，更少量的试剂，并且具有比传统细胞生物学测定更低的误差可能性。然而，在一些情况下，细胞群体可能难以基于复合纤维素进行鉴定和定量lar功能，需要额外的故障排除;我们强调协议中的某些情况。我们在癌症模型中使用对药物的反应来证明这一应用;然而，它可以更容易地应用于其他生理过程。该协议允许人们识别共培养系统内的亚群体，并表征每个人对外部刺激的特定反应。
Introduction
基于细胞的测定在基础研究和药物开发环境中已经成为一个主力。然而，这些标准测定的局限性随着体外和临床数据的不一致以及大多数药物未能获得FDA批准而变得越来越明显。在这里，我们提出了一种利用定量成像同时分析异源细胞表型以响应相关的共同环境刺激的新方法。
用于测量细胞活力的传统的基于细胞的测定法包括：台盼蓝排除测定，MTT / MTS和Annexin V-FITC流式细胞术染色。台盼蓝排除测定虽然简单而便宜，但需要大量的细胞，这是耗时的，并且经常受到用户偏倚的影响1 。 MTT和MTS测定通过测量线粒体代谢率间接测量细胞活力。然而，代谢活性细胞的数量可能受到不同培养条件（如培养基或氧浓度）的影响，导致不准确的结果，并阻止细胞类型和条件的标准化2,3 。这些技术的另一个主要缺点是它们无法区分多种细胞类型 - 大多数生物系统是异源细胞的。虽然流式细胞术方法具有区分多细胞群体的能力，但需要细胞标记，动态取样具有挑战性，而当使用贴壁细胞时，该应用程序变得耗时且容易出错。
其他重要的细胞表型，包括形态学变化，发生在响应于环境刺激，但不被传统的基于细胞的测定法捕获。通过形态表征和映射样本之间的相似性，分析细胞状态是一个强大的，无偏见的工具能够提供基础和翻译研究的许多方面的新见解，包括基础细胞生物学和药物发现4 。此外，已经显示肿瘤细胞形态与肿瘤亚型5和侵袭性6相关 。因此，研究这些细胞特征以及它们如何与特定的环境扰动有关是非常有意义的。此外，人们可以使用形态特征的差异来区分共培养系统中的亚群体。荧光标记细胞具有垮台（ 即改变固有的细胞特性，耗时），因此分类细胞类型的其它方法是有利的。
基于显微镜的成像是以多重，定量和鲁棒的方式分析细胞表型的替代方法。在这份手稿中，我们运用定量成像管道来突出演变肿瘤内异质细胞群的动态变化。我们专注于非小细胞肺癌（NSCLC）细胞和肿瘤相关成纤维细胞（CAF）之间的相互作用，这是肿瘤中最常见的基质细胞类型。 CAFs参与了肿瘤的起始，进展和治疗反应;因此，在不存在CAF的情况下对肿瘤细胞进行表型测定可能会误导7,8,9 。具体来说，我们评估了CAF对肿瘤细胞对厄洛替尼的影响，厄洛替尼是靶向经常用于NSCLC临床治疗的表皮生长因子受体（EGFR）的小分子。我们利用高分辨率的筛选平台及其附带的图像分析软件进行评估;然而，为了使其他研究人员可以使用这种方法，我们还使用开源软件开发了可比较的下游协议：CellProfiler 10和CellProfiler Analyst 11 。大多数基于图像的高含量筛选测定用特定于给定仪器模型的商业化软件进行分析。结果很难在其他实验室用不同的软件复制，因为底层算法通常是专有的。使用这种基于图像的管道，测量了使用基于荧光和形态的分类来响应药物治疗的异源细胞培养物的每个亚群的细胞增殖，死亡和形态。以下协议提供了一种鲁棒的方法来探测复杂的蜂窝过程。
Protocol
细胞培养
注意：在与肺成纤维细胞（CCD-19Lu GFP）共同培养的情况下，研究NSCLC细胞（H3255）对EGFR靶向药物埃罗替尼的表型反应。对于相同的数据集，进行了两个细胞群体（H3255和CCD-19Lu GFP）的基于荧光和基于形态的分类，以证明它们的一致性。然而，只需要使用一种分类方法，并应根据应用选择。
- 准备500毫升补充有10％热灭活FBS和1％青霉素/链霉素的RPMI-1640。
注意：生长培养基和补充剂可以根据其他细胞类型的要求进行更换。 - 培养细胞在10cm 3平板中作为纯种群在5％CO 2中于37℃，并以适当的比例每3-4天通过。
2.细胞的制备
- 要制备细胞悬液，取出c并用5mL 1X磷酸盐缓冲盐水（PBS）洗涤细胞。
- 从PBS中吸出，加入1 mL加入到37°C的0.05％胰蛋白酶，并在37℃下孵育5分钟（或直到细胞不再粘附于平板）。
- 为了中和胰蛋白酶，向H3255和CCD-19Lu GFP板中加入4mL RPMI，并将细胞溶液移至15mL锥形管中。
- 在室温下以150×g离心细胞5分钟。
- 将上清液吸出并将细胞沉淀重悬于10mL RPMI培养基的锥形管中以准备播种。
电池电镀
- 计数细胞以使种子标准化。
- 反转管将溶液中的细胞混合。将10μL每个细胞悬浮液移至微量离心管中，并加入10μL台盼蓝。
- 将10μL此溶液移至细胞计数载玻片并插入自动细胞计数器。
注意：其他细胞计数方法（ 即。血细胞计数器）。 - 重复细胞计数至少三次，或直到获得的计数一致。
- 将种子细胞置于多孔板上
注意：种子的数量取决于将被成像的时间点数。或者，如果细胞用组蛋白-2B-GFP（或提供足够的核分割的其他稳定的慢病毒核染色体）转导，则可以随时间重新成像一个平板。- 在100μLRPMI中种子总共1500个细胞/孔。准备三种细胞悬浮液来平板三种细胞群体：H3255，CCD-19Lu GFP和50％H3255 + 50％CCD-19Lu GFP。对于每个时间点，使用相同的板设置执行一式三份。
注意：应针对每个细胞系和板尺寸优化初始播种密度。 - 使用多通道移液器在3×96孔板中的种子细胞。在37℃，5％CO 2下孵育细胞O / N。
- 在100μLRPMI中种子总共1500个细胞/孔。准备三种细胞悬浮液来平板三种细胞群体：H3255，CCD-19Lu GFP和50％H3255 + 50％CCD-19Lu GFP。对于每个时间点，使用相同的板设置执行一式三份。
药物剂量
- 将DMSO加入厄洛替尼粉末中以制备10mM浓度的最终埃罗替尼储备溶液。
注意：药物可以根据需要进行取代。 - 用细胞培养基稀释药物至最终浓度的2倍，最终浓度为10μM，以1:10的比例连续稀释4次，共5种药物浓度和无药物对照。
注意：如果最高药物剂量的DMSO浓度等于或超过0.1％v / v，则无药物对照应含有等量的DMSO，以确保观察到的观察效果不是由于DMSO毒性所致。 - 将100μL药物溶液移至相应的孔中，以最终药物浓度1×稀释于培养基中。
图像采集
注意：在第0,2和3天获得图像。根据正在研究的细胞类型和细胞过程，o可以研究期望的时间点。
- 染色细胞准备成像。
- 以下列终浓度补充染料溶液。
- 对于基于荧光的分类，在PBS中制备5μg/ mL核染色和5μM死细胞染色（ 参见表材料 ）。
- 对于基于形态学的分类，在PBS中制备5μg/ mL核染色剂，5μM死细胞染色剂和5μM细胞染色剂（ 参见表材料 ）。
- 向每个孔中加入20μL染料溶液。在37°C孵育30分钟防光。
- 以下列终浓度补充染料溶液。
- 优化和获取图像。
- 从培养箱中取出孔板，用70％EtOH擦拭板的底部，并将板放置在成像室中。
- 在“设置”选项卡中，单击“通道选择”下的“+”按钮添加相应的通道图像（ 即明场，核染色，死细胞染色，GFP和RFP信号）。
- 单击“布局选择”，并从0μm开始，以2μm开始，以2μm为单位进行z层叠的测试图像，并以20μm结束，以识别对焦平面。为“高度”下的每个通道输入此距离。
- 点击每个通道下的“快照”来评估强度和优化曝光时间。根据需要在“时间”下输入更高或更低的值。
- 在屏幕的右侧，在板原理图上突出显示适当的孔（待成像）。在下面的井示意图中，以类似的方式突出显示要成像的二十五个场。
- 在“运行实验”下，在左侧选项卡上的“板名称”中标识板，然后单击“开始”按钮（下面），以10X物镜运行图像采集协议，以在上述通道中生成灰度TIFF图像。
注意：分步骤成像协议的拍卖在仪器之间可能不同，但参数值仍应进行优化。 - 在第二天和第三天重复成像。
图像分析
注意：所有图像分析都是使用专有软件进行的。然而，由于这不是公开的，因此还在CellProfiler 2.2和CellProfiler Analyst 2.0上设计了可比较的分析，下面列出了一个简短的协议，并提供了详细的协议作为补充材料（测试图像已经被加载到管道中以测试工作流程）。来自该实验的图像在每个井基础上分组在一起，使得每个孔可以单独加载。下面的CellProfiler管道包含一个正则表达式，从每个图像文件名分析元数据信息，并允许图像按渠道进一步分组。
- 下载并安装开源软件：'CellProfiler9;和“CellProfiler Analyst”。在安装期间接受所有默认选项和加载项。
- 将异源细胞培养物分类为亚群体。
- 打开“CellProfiler”软件，单击“文件”| “进口管道”| '从文件'，并选择相应的文件（“荧光分类cppipe”进行基于荧光的分类，或“Morphology_Classification.cppipe”用于基于形态的分类）。
注意：这些文件包含基本图像处理和核/细胞分割的协议。由于这些细胞上存在不均匀的Hoechst染色体，细胞核被分割和扩大，然后用于掩蔽图像以允许以更低的强度分割细胞核。- 选择基于荧光的分类流程，根据eGFP强度将细胞分类为H3255或CCD-19Lu，并计算形态特征。
注意：CCD-19Lu细胞用GFP转导。 - 选择基于形态的分类管线来细分核和细胞，并计算形态特征。
注意：需要进一步下游分析“CellProfiler Analyst”进行分类。死细胞通过基于荧光的分类进行分类。
- 选择基于荧光的分类流程，根据eGFP强度将细胞分类为H3255或CCD-19Lu，并计算形态特征。
- 点击屏幕左下角的“查看输出设置”。选择要分析的图像文件的位置（“默认输入文件夹”）和提取数据的目标位置（“默认输出文件夹”）。
- 点击“分析图像”开始分析。分析完成后，点击“确定”，进入“默认输出文件夹”，查看计算的数据。
注： 补充文件1-CellProfilerProtocol.pdf中提供了示例参数值和图像。 - 对于基于形态的分类（仅），请执行以下操作。
- 打开“CellProfiler Analyst”软件，选择CellProfiler中生成的数据库属性文件，然后单击“打开”。
- 点击屏幕左上角的“分类器”选项卡。要从实验中调用随机单元格图像，请单击“提取”;图像将显示在未分类的窗口中。
- 通过选择并将细胞拖动到适当的储存箱（见补充文件2-Pipeline.pdf ： 图B ），将细胞手动分类为阳性（H3255）或阴性（CCD-19Lu）。每个子群体分类至少50个细胞后，单击“培训分类器”，然后单击“检查进度”。
注意：细胞通过用户监督的随机森林机器学习算法12进行分类 。如果准确度不高于批准阈值（> 90％），可能需要返回并优化“CellProfiler”上的细胞分割，否则细胞可能不适合分类通过形态学。 - 点击“全部评分”为每个子群体生成一个具有细胞计数的表格。
- 打开“CellProfiler”软件，单击“文件”| “进口管道”| '从文件'，并选择相应的文件（“荧光分类cppipe”进行基于荧光的分类，或“Morphology_Classification.cppipe”用于基于形态的分类）。
Representative Results
我们生成了一个由25个字段/孔，54个孔/板（3个细胞群体×6个药物浓度×3个重复）组成的图像组，跨越三个平板，共有4,050个单独的图像。在实验过程中生成的图像集使用专有软件（参见材料表）进行分析，以提取细胞的各种定量性质（ 即形态学，荧光），然后可用于分类细胞亚群。然而，由于所使用的商业软件的访问有限，因此创建了CellProfiler和CellProfiler Analyst中可比较的下游管道。
异种细胞分类为亚群
基于DNA染色（这里是Hoechst）鉴定和分割细胞核，并且基于荧光或莫尔将细胞群分类形态学（ 图1 ）。对于基于荧光的分类，成纤维细胞（CCD-19Lu）预先用GFP-慢病毒转导。测量每个核的GFP强度水平，并且计算高于允许阈值（基于背景信号）的那些被分类为CCD-19Lu，而下面被鉴定为肿瘤细胞（H3255）。对于基于形态学的分类，细胞先前用无毒细胞染色剂染色（参见材料表），并将其用于鉴定和分割细胞质。训练每个人口约50-100个细胞的机器学习算法。鉴定出种群间的形态特征差异显着，然后将其用于设计线性分类器以区分CCD-19Lu和H3255细胞。荧光和形态分类方案在区分两个细胞群体时一致为97.4％（n = 1403）在药物治疗条件（1μM埃罗替尼）中，92.5％（n = 916）与未治疗条件相同（ 图2 ）。
子群体表型分析
除了区分细胞类型之外，我们旨在表征每个亚群的表型性质。多重测试可节省时间和试剂，增加一致性，并提供有关正在研究的系统的其他信息。有许多潜在的表型产出，应根据感兴趣的问题选择它们。在这里，研究了响应于厄洛替尼治疗的细胞形态和活力状态的变化。经过三天的药物治疗后，观察到H3255细胞的核面积减少和细胞面积的增加（ 图3A ）。核区域之间的平均差异通过双侧型-2（相等方差） t检验（ p = 7.92×10 -16 ），发现“无药物”和“药物”治疗的人群具有统计学显着性。我们假设这种观察是对药物治疗施加的压力的细胞反应。
研究药物是否具有细胞毒性（ 即 ，随着时间的推移死亡细胞数量的增加）或细胞抑制（ 即随着时间的推移而减少细胞出生数）对细胞的影响，因为这具有深刻的临床影响。例如，细胞抑制药物作用诱导生长停滞，但不能消除肿瘤中的细胞，因此一旦药物被除去，癌细胞就有可能重新启动细胞增殖。药物作用常常是上下文，浓度和细胞类型依赖性。我们以前观察到厄洛替尼在一种细胞类型中引发细胞毒性反应，同时显示ac在另一个13的 ytostatic反应。
传统的活力测定法输出相对细胞数，因此，不区分生长停滞和细胞死亡。在此，基于碘化丙啶染色鉴定死细胞（ 图3B ）。观察到厄洛替尼对H3255细胞的细胞毒性和细胞生长抑制作用，药物治疗后死亡人数增多，出生人数减少（ 图3C ）。值得注意的是，由于细胞清除，死亡细胞数量在第1天后下降。 CCD-19Lu细胞不受药物的影响。该平台的另外一个优点就是生成定量数据。例如，在我们的共培养实验中，在没有或与厄洛替丁组合的三天后，发现1,118（75.8％）H3255细胞的初始亚群为2,817（87.9％）或396（57.2％）ib处理（ 图4 ）。因为我们可以产生实际细胞计数而不是相对百分比（如流式细胞术方法），我们得出结论，药物治疗期间组成的变化是由于H3255细胞减少而不是CCD-19Lu的增加。由于细胞清除，死亡率可能被低估是值得的，这是难以通过实验评估，并且可能在细胞类型之间不同。
图1：图像分析协议概述。两个潜在的下游图像分析管线可以使用基于形态的分类或基于荧光的分类对异源细胞群进行分类。刻度棒=100μm。 请点击她e查看这个数字的较大版本。
图2：形态与基于荧光的分类之间的一致性。 （ A ）显示两个分类协议重叠的一致图。使用基于形态和荧光的分类将相同的细胞分类为H3255。两种方案是一致的，分类为未处理细胞的97.4％（n = 1403）和用厄洛替尼处理的细胞的92.5％（n = 916）（注意：白色区域太小而不可视化）。 （ B ）10X图像描绘了基于荧光和基于形态的分类之间良好和差的一致性的示例。白色箭头指出在平台之间不一致的细胞。输入图像：蓝 - 核（Hoechst）;绿色 - CCD19Lu（GFP）。 CLASSI化妆图像：红色 - H3255;绿色 - CCD-19Lu。刻度棒=100μm。 请点击此处查看此图的较大版本。
图3：单个实验设置的复用表型测量。 （ A ）在存在和不存在药物时，在单细胞水平上计算细胞核和细胞面积等形态特征。注意：测量小于100μm2的细胞区域被认为是碎片，并排除在分析之外。箱体图描绘了具有第一和第三四分位数范围和95％置信区间误差条的中值。 （ B ）共同培养H3255（蓝色）和CCD-19Lu（绿色）细胞，基于丙酸的强度鉴定死细胞碘化镉染色（红色），并使用10X物镜成像。刻度棒= 1 mm（顶面）; 100μm（底部图像）。 （ C ）用或不用药物治疗三天计算活细胞死亡细胞数，活细胞数量明显减少，加上厄洛替尼死亡细胞数量增加。误差条代表基于三次重复的平均值的标准误差。 请点击此处查看此图的较大版本。
图4：随时间推移的亚群动态。在含有和不含药物的第0天或第3天，含有H3255（蓝色）和CCD-19Lu（绿色）的孔的代表性10X图像。计数属于每个亚群的细胞，比例饼图显示a样本中人口组成的变化趋势。刻度棒= 1毫米（中间面板，闪烁;最后面板），1毫米（顶部图像），100μm（底部图像）。 请点击此处查看此图的较大版本。
Disclosures
作者没有什么可以披露的。
Acknowledgments
这项工作由美国国家癌症研究所（NCI）拨款U54CA143798和U54CA143907资助，分别在达纳 - 法伯癌症研究所和南加州大学建立物理科学 - 肿瘤学中心（PS-OC）。 SM Mumenthaler获得了PS-OC跨网络奖，支持了这项工作。
衷心感谢我们的慈善支持者，特别是Stephenson家族，Emmet，Toni和Tessa，他们捐赠了Operetta HCS平台。我们还要感谢J. Foo指导，应用分子医学中心团队成员：D. Agus临床指导和指导，K. Patsch与实验设计进行有意义的讨论，R. Rawat为图像分析协议提供帮助，J. Katz与Operetta的技术协助，以及P. Macklin和D. Ruderman进行了有益的讨论和反馈。
Materials
| Name | Company | Catalog Number | Comments |
| RPMI-1640 | Corning | 10-040-CV | cell culture medium |
| Fetal Bovine Serum (FBS) | Gemini | 100-106 | medium supplement |
| Penicillin/streptomycin | Gibco | 15-140-122 | medium supplement |
| TC20 | Biorad | 1450102 | cell counter |
| TC20 slides with trypan blue | Biorad | 1450003 | cell counter slides |
| 96-well plates | Corning | 3904 | clear bottom black plates |
| Erlotinib | LC Laboratories | E-4007 | |
| DMSO | VWR | 317275-100ML | solution to resuspend drug |
| Hoechst | Invitrogen | H21491 | nuclear dye |
| Propidium Iodide (PI) | Invitrogen | P1304MP | dead cell stain |
| Cell Tracker Orange CMRA | Life Technologies | C34551 | whole cell stain |
| Operetta high content imaging system | Perkin Elmer | ||
| CellProfiler | Broad Institue | version 2.2.0 (rev 9969f42) | http://cellprofiler.org/releases/ |
| Cell culture incubator | Any cell culture incubator will be suitable - cells were cultured under 37 ºC at 5% CO2. | ||
| 15 mL Falcon conical tubes | Falcon | 14-959-53A | |
| 10 cm2 cell culture plates | TPP | 93040 |
References
- Kanemura, Y., et al. Evaluation of in vitro proliferative activity of human fetal neural stem/progenitor cells using indirect measurements of viable cells based on cellular metabolic activity. J Neurosci Res. 69 (6), 869-879 (2002).
- Plumb, J. A., Milroy, R., Kaye, S. B. Effects of the pH dependence of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide-formazan absorption on chemosensitivity determined by a novel tetrazolium-based assay. Cancer Res. 49 (16), 4435-4440 (1989).
- Hsu, S., et al. Green tea polyphenols induce differentiation and proliferation in epidermal keratinocytes. J Pharmacol Exp Ther. 306 (1), 29-34 (2003).
- Caicedo, J. C., Singh, S., Carpenter, A. E.
- Sero, J. E., et al. Cell shape and the microenvironment regulate nuclear translocation of NF-kappaB in breast epithelial and tumor cells. Mol Syst Biol. 11 (3), 790 (2015).
- Chen, J. F., et al. Subclassification of prostate cancer circulating tumor cells by nuclear size reveals very small nuclear circulating tumor cells in patients with visceral metastases. Cancer. 121 (18), 3240-3251 (2015).
- Ohlund, D., Elyada, E., Tuveson, D.
- Paraiso, K. H., Smalley, K. S.
- Kalluri, R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 16 (9), 582-598 (2016).
- Carpenter, A. E., et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 7 (10), R100 (2006).
- Jones, T. R., et al. CellProfiler Analyst: data exploration and analysis software for complex image-based screens. BMC Bioinformatics. 9, 482 (2008).
- Dao, D., et al. CellProfiler Analyst: interactive data exploration, analysis and classification of large biological image sets. Bioinformatics. 32 (20), 3210-3212 (2016).
- Garvey, C. M., et al. A high-content image-based method for quantitatively studying context-dependent cell population dynamics. Sci Rep. 6, 29752 (2016).
- Juarez, E. F., et al. Quantifying differences in cell line population dynamics using CellPD. BMC Syst Biol. 10 (1), 92 (2016).
- Mumenthaler, S. M., et al. The Impact of Microenvironmental Heterogeneity on the Evolution of Drug Resistance in Cancer Cells. Cancer Inform. 14, Suppl 4. 19-31 (2015).
