Summary
健身定量分析（QFA）是一系列估计微生物培养适应度的实验和计算方法的互补。 QFA估计的基因突变，药物或其他适用于微生物的生长治疗的效果。可以设计实验缩放从单一文化的集中分析数以千计的并行文化。
Abstract
定量健身分析（QFA）是为适应度比较在1,2,3,4平行生长的微生物培养的实验和计算的工作流程。 QFA可以应用于单一文化的集中的意见，而且是最有用的全基因组的遗传相互作用或药物屏幕调查到数以千计的独立文化。中央实验方法是独立的，稀的液体微生物固体琼脂平板上培养和定期拍下文化的接种。从每个时间点的照片进行了分析，生产定量细胞密度估计，这是用来兴建的生长曲线，从而得到定量的健身措施。文化适应度可以比较量化和排序遗传相互作用的优势或药物的敏感性。任何治疗的健身文化的影响补充到基板琼脂（如小分子物质，抗生素或营养）或适用于板的外部LY（如紫外线照射，温度），可量化的QFA。
QFA工作流程产生类似于光度法测量平行的液体培养在96孔或200井板读者获得的增长率估计。重要的是，QFA与这些方法相比有显着更高的吞吐量。 QFA文化上生长了坚实的琼脂表面，因此在增长，而不需要进行搅拌或摇晃充气。
QFA吞吐量不高，一些合成基因阵列（SGA）的筛选方法5,6。然而，由于严重之前被接种到琼脂稀释QFA文化，QFA可以捕捉更完整的生长曲线，包括指数和饱和阶段3。例如，生长曲线观察允许直接与高精度培养倍增时间估计，作为讨论以前1。
在这里，我们提出了具体QFA协议适用于成千上万的S.酵母是由机器人自动处理文化在孕育，孵化和成像。这些自动化的任何步骤可以由具有同等作用，手动程序取代，吞吐量减少，我们也提出一个较低的吞吐量手册协议。 QFA相同的软件工具可以应用在任何工作流程中捕获的图像。
我们有丰富的经验，申请QFA芽殖酵母培养学酵母 ，但我们预计，，QFA将证明裂殖酵母为研究文化同样有用S。酵母和细菌培养。
Protocol
手动，QFA协议（ 酵母菌株 ）
1。培养酵母菌株
- 96个独立的酵母菌株在200μL丰富的液体介质，在96培养皿培养。株生长饱和控制孵化温度。
- 一个96针，1/8“直径手册副本柏列达（西格玛的R-2508）第一副本柏列达在100％的乙醇浸泡，燃烧，然后离开降温消毒。
- 200μl无菌水的排列在96孔培养板。饱和的文化振荡（Eppendorf公司MixMate，30秒，750转），浸入饱和株三次消毒的针工具板含有水一次，然后稀释。
- 针工具重新消毒为1.2。
- 用于灭菌的针工具发现上出售的琼脂平板稀释文化浸入板中含有稀释的文化三倍列印转移到一个长方形的固体琼脂平板上的针工具。
- 长方形的固体琼脂平板培养在适当的温度前手动使用标准普尔机器人spImager成像。可替代的图像采集方法也可以使用，如富士LAS4000，或更便宜的数码单反相机，如果已采取审慎措施确保即使照明和成像相机相对一致的车牌定位。
- 该方法是适应平行的文化，如48，可装上一轮培养皿，如果已采取审慎措施确保在每张照片成像相机相同的角度相对较小的数字。
全自动QFA协议（ 酵母菌株 ）
2。培养酵母菌株
- 启动一个独立的酵母菌株固体琼脂上生长的矩形阵列。在这个例子中，出发菌株过温度（YDC）由合成基因阵列（SGA）的酵母删除收集敏感的查询突变。最后SGA的板块是在1536年的格式。这是16行24列代表四个殖民地/文化，384个独立的酵母菌株的复制。板布局见图1。
- 使用了96针，针直径1毫米无菌pintool成4个96孔培养板含有选择性培养基200μL（ 图1）S＆P的机器人公司BM3-SC机器人，机器人被转移了1536年的384株。 pintool洁净和无菌实现顺序在无菌水清洗两次（一次一个旋转的刷子，清除杂物），70％的乙醇（超声）和最后100％的乙醇。
- 文化发展到饱和（3天在20°C在这个例子中，参见图1）。
3。发现酵母培养
- 使用1贝克曼Biomek FX，饱和文化重新悬浮涡旋上的Variomag Teleshake（1000RPM逆时针20秒）和稀释成200μL无菌水寄希望于用无菌的96针机器人针工具（V＆P科学磁性安装pintool，2毫米针直径）约1:70。针工具的清洁和无菌，无菌水冲洗，在70％乙醇洗涤（用刷子清除杂物），最后经干燥70％的乙醇浸泡。
- 察觉到384格式使用相同的机器人针工具（ 见图1）清洗和消毒后的固体琼脂板稀释的文化。
- 板固定后，被转移到一个Cytomat温度控制，加湿的孵化器与自动传送带和出入舱口。
4。图像捕捉
- 标准普尔机器人自动化成像反复删除板块从Cytomat孵化器，消除了板盖，根据一项附上他们正是D，均匀照射下佳能EOS反叛钛135数码单反的空间，在5184 x 3456像素分辨率图像捕获，取代了板盖，并返回他们的孵化器传送带。紧随最初转移到孵化时间点零增长的测量（背景）允许一个图像被捕获。机器人处理和摄影每盘需要2分钟。
- 随着一辆满载的传送带，达到最大的板图像采集频率是一张照片，每六个小时。板培养6天（直到菌落生长饱和。 图1显示了三种典型的图像从一个时间过程，在20°C
- 用于跟踪目的，板被命名为根据孵化器的名称，温度，独特的序批数和在孵化器的位置。图片名称是串联板的名称和时间戳和自动图像采集后兴建（如K000011_030_010_2011-09-22_09-47-54.jpg，文件类型FT1， 表1）。
5。实验元数据捕获
本节中描述的元数据文件是制表符分隔的文本文件手动产生（例如，使用电子表格软件）。每次实验，实验的描述文件（FT2， 表1）是独一无二的，但图书馆描述氨酸（FT3， 表1）可以重新广泛使用一次构建。
- 实验中介绍了实验的描述文件，包含条码列（或自动生成的板块名称）（FT2， 表1），实验开始时间戳，板治疗，固体琼脂培养基，屏幕上的名字，图书馆，车牌号码（内容与图书馆多板）和重复象限数量从1536格式缩减到384格式（ 见图1）。
- 酵母菌株库库描述文件（FT3，
- 可以提供一个可选的标准基因名称文件（FT4水平， 表1）描述每个系统的ORF Y-号码识别正在筛选菌株（如YDR217C）相关标准的基因名称（如RAD9）。此文件包含两列：ORF的基因名称。
6。数据分析
QFA计算的工作流程需要访问一个相当强大的多核计算机工作站（如至强四核2.67 GHz处理器和12GB内存的戴尔Precision T3500）上安装的Colonyzer图像分析软件工具，3,4和QFA R包7 ，这两者都是网上记录，都可以免费上运行的操作系统范围。
- 运行捕获板作为输入图像的每一个Colonyzer，为每个捕获的图像中产生一个Colonyzer输出文件（FT5 表1）。 colonyzer输出文件指定的养殖密度估计，文化领域，为每个384位置上的成像板的形状和颜色。输出文件名 会自动复制源图像文件名 （如板照片K000011_030_010_2011-09-22_09-47-54.jpg（FT1加热， 表1）对应的Colonyzer输出文件K000011_030_010_2011-09-22_09-47-54.dat（ FT5， 表1））。
- 使用QFA R包，加载实验的元数据（FT2， 表1）和Colonyzer输出文件（FT5 表1）。这些数据合并到R的数据帧（QFA原始数据的文本文件“（FT6， 表1）可作为出口）作进一步的分析。
- QFA R包中包含的功能，以收集每一种文化细胞密度timecourses的，适合一概而论劳的GISTIC人口的意见和绘图模型（见图2和图3的例子）。拟合参数值被写入QFA物流参数文件（FT7 表1）。见QFA R为进一步的细节包的文档。
- QFA R包还包含质量控制功能：边缘殖民地板边缘和困难，由于提供更多的营养物质被丢弃在附近的钢板墙，文化未能SGA和基因型显示屏幕的特定标记基因的连锁与图像分析从分析中被剥离。
- 来自Logistic人口模型参数为每个文化几个定量的健身措施。其中包括最大的人口倍增率（MDR）和部门的数量从接种到饱和度（MDP）。为每个复制的文化集（例如独特的基因型），并为每个健身的定义，健身估计数汇总统计，排版uted：平均数和中位数的健身，健美的标准偏差和数量的重复观测。汇总统计数据输出QFA健身概要文件（FT8， 表1），作进一步的分析，如计算遗传相互作用分数（FT9 表1）。
7。代表结果
图2显示了一个典型的增长384 QFA增长在20°C固体琼脂由Colonyzer量化，随着细胞密度的增加大约2天之后首次成为检测文化的曲线。这种延迟可能是文化的滞后期后接种到固体培养基和细胞密度较低的检测限制的综合效应。 图3显示了从一个单一的板捕获的308类似的增长曲线图4显示了两个全基因组的比较QFA屏幕推断遗传互动优势（从1改编）。 QFA R包7还包括功能产生图3和输出研究基因型和遗传相互作用的强度（GIS）的估计排名名单，一个Q值（错误发现率（FDR）的校正p值）观测到的地理信息系统的意义。
图1。在384点的格式为机器人的酵母菌株接种计划。此过程始于1536独立文化（左）每盘。在这个典型的例子，殖民地位1,1; 1,2; 2,1和2,2（红色）四个重复相同的基因型HIS3 ::黄色KANMX文化，生长板边缘。 ，有一个增长的优势，由于缺乏竞争，因此不QFA检查。这些复制（如1,1）在96孔板接种到液体培养基使用96针的工具，接种1536菌落每次96。接种，以接种四个不同的“象限”（显示为红色，蓝色，绿色和紫色）的384个基因的缺失，每一个复制到四个不同的96孔含生长介质板。文化增长饱和（如3天在20°C）后，用水稀释，然后从一个重复的四个象限察觉到了坚实的琼脂板（右）在SGA的原板相同的模式，在384格式（颜色表示）。这个过程可以重复测试其他的复制：1,2，2,1和2,2。例如时间推移右侧的图像被抓获接种后0.5年，2年和3.5天。调整补充图1 2。 点击这里查看大图 。
018/4018fig2.jpg“/>
图2。一个单一的机器人捕获的QFA生长曲线A）QFA增长20琼脂中含半乳糖的酵母菌株的生长曲线°C。图像被抓获机器人大约每2个小时。指数相在此温度下，观察到约1.5天，与养殖密度的增加，首次成为检测接种后约2天。广义Logistic模型是适合的观测数据（灰色曲线）。 QFA R包自动装配模型的参数是：K（承载能力（非盟）），R（增长率（四-1）），G（接种密度（非盟）），V（增长对称）。 b）一个面板，日志规模与 细胞密度的绘制。 点击这里查看大图 。
图3。自动基因额定生长曲线捕获从一个单一的384格式板平行。每个面板显示细胞密度估计（红叉）和独立的文化标记基因名称或框ŸQFA过程中生长模型拟合（黑色曲线）。这个例子板成像和机器人在自动孵化成长，在20°C实验的元数据将自动列入图中的标题：名盘，板治疗，琼脂培养基和图书馆车牌号码。拟合模型的参数值也自动打印在每个面板（见图2）。边缘文化被剥夺QFA分析，留在这个小区的308文化（见QFA议定书5.4）。由于边缘文化不QFA分析，这些文化的地点通常是充满了中性的控制株（在这种情况下pGAL-HIS3）。 QFA R包的输出，这个数字的原始版本，是PDF格式，因此是无限可缩放和可检索的文本查询。:/ / www.jove.com/files/ftp_upload/4018/4018fig3.pdf“的目标=”_blank“>点击这里查看大图。
图4。两个QFA屏幕比较适应度来推断基因的相互作用。此图显示了芽殖酵母结合中性ura3Δ突变或对温度敏感的cdc13-1突变，生长在半放任温度缺失（yfgΔ）的健身比较cdc13 1（27℃）。健身计算作为最大的速度提高一倍的产品和最大倍增潜力1。大多数缺失株生长良好，与的中性ura3Δ突变结合，如预期，并具有较高的健身，但结合cdc13-1的缺失，有广泛的适应度。淡蓝色的线交叉在日ê的的中性his3Δ突变的位置，实线是线性回归模型的遗传独立性（cdc13-1预计健身yfgΔ突变体健身的ura3ΔyfgΔ突变）和虚线是平等健身行。缺失色绿，显着提高健身cdc13-1株的缺陷。删除红色显着抑制相同的缺陷。任何突变cdc13-1yfgΔ实线的垂直距离显示板条件下cdc13-1和yfgΔ之间的遗传相互作用的强度。注cdc13-1最强的抑制一些先前发现的右上角附近的检查点基因，RAD17，DDC1和RAD24和核酸EXO1。基因的缺失降低健身包括附近YKU80（编码DNA修复和端粒封盖蛋白Yku80）右下角。还注意到，一些团体共同发挥作用的基因，往往会被位于该地块上合作。强调在蓝色SPE1，SPE2和SPE3的三个例子集群：亚精胺合成的基因，RAD17，DDC1和RAD24：检查点滑动钳和检查站钳装载机，MRE11，RAD50和XRS2编码组件的MRX复杂，在DDR和EST1的参与EST3：基因参与端粒长度调节。改编自补充图1B 1，从这个情节产生的原始数据，可以从这里下载： http://research.ncl.ac.uk/colonyzer/AddinallQFA/ 。 点击这里查看大图 。
表1。电子QFA文件。列出的所有文件（除FT1）是制表符分隔的文本文件，可以使用任何文本编辑器或电子表格软件建造或阅读。关于构建QFA元数据文件和QFA输出的精确解释的进一步详细信息，请参阅的QFA R包文件7和Colonyzer软件文档3。
表2。评估文化健身方法。总结的特点和要求ments的评估微生物培养健身方法。审查布隆贝格8包含其中一些更详细的比较。 点击这里查看大图 。
Discussion
QFA是在许多感官的直接后裔三立台规模的微生物技术：文化系列稀释现场测试9，生长曲线的测定充气液体培养和副本电镀13。这三种方法进行了总结和与QFA和其他高通量技术在表2比较。而稀释系列现场测试定义为应变能力，形成一系列种植范围从接种稀释的文化殖民地的健身，QFA措施，由一个单一的文化反复观察应变健身构建一个增长曲线。量化健身从单一的文化，让更多株进行检查，同时在相同条件下。文化复制阵列允许重复的应变集合，最有效的测试时，测试在不同的环境或遗传背景的文化健身差异。 QFA协议介绍这里使用昂贵的机器设备，复制，接种，孵育和图像琼脂平板上，适当的观察数以千计的独立文化（机器人QFA， 表2）的增长。然而，由于QFA是既定的，传统的实验室技术为基础，它也可以进行更便宜，取代手工步骤（手动QFA， 表2）机器人援助。手动QFA涉及复制和接种到琼脂平板上培养使用手册针板和摄影的孵化器的工具和手动传输。可应用于相同的计算分析工作流程，无论从实验设计产生的生长曲线。
QFA健身的估计似乎是比较精确。在图4中，突出四个功能相关的基因缺失，例如集群的平均适应度。每个功能相关的基因类成员的接近两independ耳鼻喉科遗传背景（ura3Δ和cdc13-1）表示QFA健身估计的可重复性。例如，只有3个基因缺失（走出了一条可能的4300）分开保守的MRX复杂的三名成员。
QFA为测试微生物菌种的生长特性，高吞吐量的替代品，包括SGA的5,6，竞争屏幕11,12条码库和屏幕捕捉光密度动力学分光板读者10总结在表2和更全面的描述别处8 。 QFA和SGA板，文化固体琼脂表面（琼脂为基础的方法， 表2）成长由机器人可以快速，轻松地处理。固体琼脂表面培养良好的透气性，整个生长和细胞可以生长在固定的文化，让交际的微生物生长在一个社区14。原封不动，微生物群落ffect自己的微环境，如乙醇分泌毒素，并可能在细胞之间的信号。然而，连续的混合液体培养，要达到足够的通风，人为破坏微生物群落和它们的微环境，这可能会影响他们的成长模式。交叉污染是在固体琼脂方法的关注，因为有较少的机会进行不同文化之间的细胞的液体飞溅。如果外国空气传播，在固体琼脂检测微生物污染时，它往往可以通过目测固体琼脂平板进行检测，并占或删除。
QFA吞吐量明显比两种方式并行的液体方法。首先，QFA（SGA）的板，接种培养挤在一起更密集，给每盘更独立的文化。在QFA通常有308文化，不计非实验性的边缘文化（ 图1）COM相比，每盘96或100文化与平行的液体培养。其次，QFA实验还可以扩展到更大的板材数量。 QFA在一个单一的实验分析板的数量，而只由孵化器或温暖的房间里，允许的最小图像采集频率和达到的最大捕获频率的可用空间的限制，在液体中的生长实验板的数量是强烈有限酶标仪（通常是一个或两个板块），或连接到酶标仪（通常为25-50板）的堆的容量。最近，我们进行完全自动化的QFA上123板，每个实验308实验的文化，同时给予37884文化。我们估计，在液体中越来越多的独立文化的最大实现数是96（文化/板）×50（板/堆叠器）= 4800，这是QFA吞吐量比我们少十倍左右。液体培养屏幕的另一种方法是利用许多自动化GROWTH设备并行（ 表从审查布隆贝格8 1），每有一个或两个板块的能力。在每个设备的温度控制是独立的，这样的条件是不相同的，但假设它们，这个工作流程将需要至少190个此类设备QFA吞吐量（假设200％液体培养设备）相匹配。
总之，QFA是一个高质量的工作流程，可以有效地聚集在数量增长的表型集中两个小规模试验和高通量筛选。它是足够灵活，有效地应用于机器人设备的要求不一。 QFA工作流程的计算组件的基础上的自由，开放的源代码。
Disclosures
没有利益冲突的声明。
Acknowledgments
我们非常感谢我们的实验室和老龄化，支持和有益的讨论营养（CISBAN）集成系统生物学中心的所有成员。这项研究是由生物技术和生物科学研究理事会（BBSRC）（BB/C008200/1）和威康信托基金会（075294，093088）。
Materials
| Name | Company | Catalog Number | Comments |
| Replica Plater | Sigma | R-2508 | 96 pin manual pintool with 1/8 " diameter pins |
| Mix Mate | Eppendorf | 5353 000.014 | |
| Biomek FX | Beckman | A31842 | |
| Teleshake | Thermo Scientific | 50095890 | Installed on Biomek FX |
| BM3-SC | S&P Robotics Inc | BM3-SC | 192 shelf rotating carousel |
| spImager | S&P Robotics Inc | spImager | High resolution manual imaging |
| spImager with Cytomat | S&P Robotics Inc | Custom | Temperature controlled automated high resolution imaging |
| Cytomat 6001 with heat exchanger | Thermo Scientific | 51022222 | Attached to spImager |
| Ecoline RE207 | Lauda | RE207 | Attached to Cytomat |
| spImager with carousel | S&P Robotics Inc | Custom | Automated high resolution imaging |
| Robot pin tool fixture for FP12pins | V&P Scientific | AFIX96FP12 | N/A |
| 96 x FP12 pins | V&P Scientific | FP12 | 50.4 mm long, 17 mm exposed pin length |
| Docking Station for Pin Tool | V&P Scientific | VP425 | Docking station base removed to allow fan drying of pins |
| Pin Cleaning Brush | V&P Scientific | VP425 | N/A |
| Replica Plater | Sigma | R-2508 | Manual pin tool |
| Nunc OmniTray with Lid | Nunc | 734-0490 | N/A |
| 96 well sterile polystyrene plates | Greiner BioOne | 655161 | N/A |
| 96 well sterile polystyrene lids | Greiner BioOne | 656171 | N/A |
Table 3. Specific reagents and equipment. |
|||
References
- Addinall, S. G., Holstein, E., Lawless, C., Yu, M., Chapman, K., et al. Quantitative Fitness Analysis shows that NMD proteins and many other protein complexes suppress or enhance distinct telomere cap defects. PLoS Genetics. 7 (4), e1001362 (2011).
- Addinall, S. G., Downey, M., Yu, M., Zubko, M. K., Dewar, J., et al. A genomewide suppressor and enhancer analysis of cdc13-1 reveals varied cellular processes influencing telomere capping in Saccharomyces cerevisiae. Genetics. 180, 2251-2266 (2008).
- Lawless, C., Wilkinson, D., Young, A., Addinall, S., Lydall, D. Colonyzer: automated quantification of micro-organism growth characteristics on solid agar. BMC Bioinformatics. 11, 287-28 (2010).
- Colonyzer - Image analysis software [Internet]. , Newcastle University. Newcastle, UK. Available from: http://research.ncl.ac.uk/colonyzer/ (2012).
- Tong, A. H., Evangelista, M., Parsons, A. B., Xu, H., Bader, G. D., et al. Systematic genetic analysis with ordered arrays of yeast deletion mutants. Science. 294, 2364-2368 (2001).
- Costanzo, M., Baryshnikova, A., Bellay, J., Kim, Y., Spear, E. D., et al. The genetic landscape of a cell. Science. 327, 425-431 (2010).
- qfa - An R Package for Quantiative Fitness Analysis [Internet]. , University of Vienna. Vienna, Austria. Available from: http://qfa.r-forge.r-project.org/ (2012).
- Blomberg, A. Measuring growth rate in high-throughput growth phenotyping. Curr. Opin. Biotechnol. 22 (1), 94-102 (2011).
- Maringele, L., Lydall, D. EXO1-dependent single-stranded DNA at telomeres activates subsets of DNA damage and spindle checkpoint pathways in budding yeast yku70Δmutants. Genes and Dev. 16 (15), 1919-1933 (2002).
- Warringer, J., Blomberg, A. Automated screening in environmental arrays allows analysis of quantitative phenotypic profiles in Saccharomyces cerevisiae. Yeast. 20, 53-67 (2003).
- Hillenmeyer, M., Fung, E., Wildenhain, J., Pierce, S., Hoon, S., Lee, W., Proctor, M., St. Onge, R., Tyers, M., Koller, D., Altman, R., Davis, R., Nislow, C., Giaever, G. The chemical genomic portrait of yeast: uncovering a phenotype for all genes. Science. 320 (5874), 362-365 (2008).
- Smith, A. M., Durbic, T., Oh, J., Urbanus, M., Proctor, M., Heisler, L. E., Giaever, G., Nislow, C. Competitive Genomic Screens of Barcoded Yeast Libraries. J. Vis. Exp. (54), e2864 (2011).
- Lederberg, J., Lederberg, E. M. Replica plating and indirect selection of bacterial mutants. J. Bacteriol. 63, 399-406 (1952).
- Queller, C. Behavioural ecology: The social side of wild yeast. Nature. 456, 589-590 (2008).
