我们已开发出全面且无偏倚的全基因组筛选方法,以研究基因-药物和基因-环境之间的相互作用。本文介绍了对这些突变体文库进行筛选的方法。
我们已开发出全面且无偏倚的全基因组筛选方法,以研究基因-药物和基因-环境之间的相互作用。本文介绍了对这些突变体文库进行筛选的方法。
得益于新一代测序技术的进步,我们几乎每天都能获得新的基因组序列。这一进展的速度正在加快,预示着更深入、更广泛的应用前景。面对如此迅猛的发展,快速、并行地定义基因功能的方法变得愈发重要。酵母和大肠杆菌(E. coli)的全基因组缺失突变体文库已被广泛用作基因功能表征的重要工具,但该方法难以扩展:目前的基因删除策略要求对基因组中成千上万个基因逐一进行敲除并验证,只有完成这一繁重工作后,才能开展高通量的表型分析。在过去十年中,本实验室优化了一系列可并行操作的、小型化、高通量的全基因组筛选方法。这种并行化得以实现的关键在于将DNA“标签”或“条形码”引入每个突变体中,条形码可作为突变的替代标记,通过检测条形码的丰度即可评估突变体的适应性。在本研究中,我们旨在填补DNA序列信息与条形码标记的突变体文库之间的空白。为此,我们提出了一种转座子破坏与条形码标记相结合的新方法,使新测序但尚未充分表征的微生物也能适用于并行条形码筛选。为展示该方法的应用,我们构建了一个新的白色念珠菌(Candida albicans)条形码标记的破坏突变文库,并描述了如何利用基于微阵列和基于新一代测序的平台,在单次实验中检测10,000至1,000,000个基因-基因及药物-基因相互作用。
1. 背景信息
有多种方法可以构建携带条形码标签的突变体。目前的金标准是由多个实验室组成的联合团队创建的酵母基因敲除(YKO)文库,该文库于2002年完成。 1自最初的酵母基因敲除文库问世以来,研究人员已构建了其他酵母文库;这些文库基于不同的菌株背景,采用过表达载体,并扩展至其他微生物中,例如 E. coli 2与此同时,构建条形码shRNA文库的工作正在快速推进,事实上,这些哺乳动物文库的许多设计原则借鉴自酵母。为了展示条形码转座子如何成为一种快速且广泛适用的系统性突变体文库构建策略,我们重点介绍近期在人类真菌病原体中构建的一个文库, 白色念珠菌. 我们对念珠菌的研究是基于条形码筛选技术的成功应用 S. cerevisiae,此处用作示例生物。该样本方案经少量修改后,可用于筛选任何能够在悬浮培养中生长的生物。由于极少生物具备高效转化和有效有丝分裂重组所必需的高突变率,以产生完全缺失突变体,因此我们开发了一种利用转座子诱变的方案 体外 对基因组DNA文库进行诱变,然后将这些带有条形码的基因组片段转化入 白色念珠菌 3, 4. 受原始YKO文库在基因网络本质基础研究中取得成功及其重要作用的启发 5-8,全基因组范围的单倍剂量不足 9药物靶点与作用机制 10,11以及基因组中所有基因的必需性 12 我们预计,将这一方法拓展至其他微生物将取得极为丰硕的成果。
以下方案假设所需突变体文库已经构建完成(例如YKO或Candida albicans基因敲除文库),并以单独保存的菌株形式提供。有关菌株构建的详细描述,请参见1,13,14。
2. 将单个突变体合并为一个混合池
3. 实验池生长
该步骤如图1所示。
注意:在新构建的菌株混合池中,始终需收集起始细胞样本(即“T0 时间点”),以评估初始菌株的代表性。具体操作为:从冻存的分装样品中直接取 1-2 OD600 的菌液加入 1.5 ml 离心管中,并按以下所述方法进行处理。
4. 基因组DNA提取、PCR及微阵列杂交或测序
5. 芯片分析(使用Candida albicans基因破坏文库获得的结果示例见图2)
6. 通过测序评估条形码酵母菌株的适合度
注:考虑使用Bar-seq作为芯片杂交的替代方法。随着高通量测序成本的降低,利用高通量测序直接检测标签丰度正变得可行,并且在许多情况下更具成本效益18。通过该方法,扩增的PCR产物被直接以“读数”形式定量,而非通过杂交至芯片后检测信号强度。这消除了因标签交叉污染、信号饱和或极高或极低信号强度所导致的假阳性与假阴性结果。此外,通过添加一段4–8个碱基的DNA索引序列,可在测序前将多个实验样本合并处理19。由于酵母条形码长度为20 bp,仅需一次26–28个碱基的双步测序读长,即可同时获取多重索引和唯一性条形码,从而实现超过100重的极高通量多重检测。在撰写本文时,Bar-seq相比条形码微阵列具有成本优势,且Bar-seq本身具备高度灵活性:随着每轮运行的读长数量增加,多重检测水平也可相应提高,从而进一步降低成本。来自所有主要测序平台厂商的若干“中通量”测序仪将进一步推动Bar-seq的普及,测序有望成为首选的检测手段。
该方案也已在Illumina HiSeq2000上验证。
最近,Gresham 等20 的一项研究出色地展示了如何利用Bar-seq来解决酿酒酵母(Saccharomyces cerevisiae)生长调控中的一个基本生物学问题,并概述了若干重要的实验设计与解释原则。
7. 池化筛选数据的验证
任何功能基因组学筛选的结果都应通过在独立培养中使用单个菌株进行验证。由于每个实验在敏感菌株数量上可能不同,因此选择需要验证的候选菌株数量具有一定的任意性。作为参考,可根据 log2 比值或 z 分数对最敏感的菌株进行排序,并对排名前 25–50% 的候选菌株进行验证(通常相当于池中所有菌株均值的 2–3 个标准差),这在成本与效益之间达到了良好的平衡。验证实验可在任意容器中进行,但我们通常在 96 孔板中进行,使用 100 μl 培养基,初始接种量为 0.06 OD600,在摇床分光光度计中培养 5 代,每 15 分钟测量一次(见图 4)。
8. 代表性结果
基因组范围筛选完成后,需对芯片数据进行标准化处理,并将每种菌株在处理组与对照组之间的表型进行比较(例如通过比较微阵列信号强度或测序读数/菌株)。随后,可将数据导入 Excel 文件进行操作,其中基因按对照组/实验组的 log2 比值进行排序。在此类排序中,log2 比值越负,表明该菌株对测试条件越敏感。这些 Excel 文件可使用多种绘图软件进行可视化。我们发现最简便的方式是以 log2 比值作为 Y 轴,以基因或开放阅读框(ORF)名称作为 X 轴进行作图。如图 2a 所示,展示了经克霉唑(一种已知的抗真菌药物)处理后的结果图。所有在 log2 比值为 -2 时显著敏感的菌株以红色标出;通常我们会通过个体生长实验进一步验证这些菌株,检测各突变体在相同药物浓度下的生长情况。在此示例中,共标出 4 个菌株:NCP1、ERG2 以及两个独立的 ERG11 等位基因,后者是克霉唑已知的作用靶蛋白。这 4 个基因均直接参与麦角固醇生物合成途径,该途径在酵母中相当于胆固醇合成。例如,NCP1 编码一种 NADP-细胞色素 P450 还原酶,参与麦角固醇的生物合成,并与 Erg11 相互关联且协同调控。此实例说明,通过这种无偏筛选方法不仅能够识别出已知的药物靶点(Erg11),还能发现靶通路中的其他关键组分。最后,图中以红色标出的若干基因可能参与麦角固醇的生物合成,或涉及其他不同的生物学过程。如前所述,任何在混合筛选中被检测为敏感的菌株,都应通过个体生长实验进一步验证其敏感性。在图 2b 所示的示例中,基于其相较于野生型亲本菌株 BWP17 的生长抑制情况,确认了四个菌株对克霉唑敏感。这些个体生长曲线揭示了此类混合基因-药物筛选的一个重要特征:即某一菌株的绝对排序并不一定准确反映其真实的敏感程度。此外,图 2b 还体现了为每个基因设置多个等位基因的价值;在此例中,两个 erg11 敲除突变体表现出略有差异的敏感性。将这些突变的性质与敏感程度相关联,可为进一步理解药物的作用机制提供额外线索。

图1. 混合培养实验与条形码检测的工作流程。将解冻的混合细胞等分试样接种至培养物中(步骤1),随后在机器人操作(选项A)或手动操作(选项B)条件下培养至所需增殖代数(步骤2)。通过离心收集细胞(步骤3),从收集的细胞中提取基因组DNA(步骤4),分别扩增上游标签(uptags)和下游标签(downtags)(步骤5),然后杂交至芯片(步骤6a)或直接进行测序(步骤6b)。

图 2. 在实验方案的某些时间点采集的示例数据。(A)筛选结果的示例数据(改编自13)。在含有克霉唑或DMSO(对照)的条件下,将带标签的突变体文库培养20代。计算Log2比值(对照强度/处理强度)并以基因为横坐标作图。高度敏感的菌株(红色)包括克霉唑的已知靶点ERG11p。需要注意的是,该检测通常除了化合物的实际靶点外,还会发现其他敏感突变体。这些突变体通常为与靶点存在合成致死效应的突变体、参与普遍应激/处理反应通路的突变体,或为验证失败的假阳性结果。(B)验证实验的示例数据(改编自13)。可通过单独培养菌株并将其生长情况与野生型(黑色)进行比较,以验证混合培养实验的结果。

图3. 从条形码混合检测中扩增出的扩增子结构,用于微阵列杂交或条形码测序。文库中每个突变体所产生的扩增子包含用于基因组整合的同源序列(蓝色区域,标记为ATG和TAA)、独特的条形码(标记为AG,以黑色短横线表示)。对于微阵列杂交,使用蓝色的通用引物扩增出一段60 bp的探针用于微阵列杂交。对于条形码测序,在PCR反应中使用延伸引物,上游引物由编码Illumina接头的序列(红色条带)、6个碱基的索引(交叉条纹)和蓝色通用引物组成,下游引物则使用相同的复合引物(不含6碱基索引)。

图4. 个体生长实验,用于:1)对野生型酵母进行化合物初筛,以确定全基因组筛选的合适剂量;2)验证全基因组筛选的结果。(A)在96孔平底板中加入100 μl、OD值为0.062的细胞悬液。每个孔可包含相同菌株(用于剂量确定)或不同菌株与药物组合(用于验证实验)。加入2 μl化合物(通常溶于DMSO),在30°C下持续振荡培养16–20小时。DMSO的终浓度不应超过2%。在此示例中,每孔的生长曲线以黑色表示,对照生长曲线以红色表示。(B)多个初筛结果的高分辨率图像,使用一种示例药物进行测试,并相互叠加。在此滴定系列中,紫色剂量可获得IC10-15,适用于缺失型谱型分析(HIP和HOP)。由于在较高光密度下存在非线性效应,Tecan(或任何类似酶标仪)测得的OD值必须使用“传统”1 mm光程比色皿获得的OD值进行校准。
本文概述了一种实验方案,经过适度修改后,该方案可轻松适用于多种不同微生物的条形码突变体文库,从而构建带标签的突变体文库。需要强调的是,尽管我们报道的是针对致病性酵母菌C. albicans的标签化转座子诱变方案,但非常相似的方案亦可适用于多种单细胞真菌。经修改后,该方案在细菌中同样有效13,目前多个其他真菌和细菌基因组的文库正在构建中。目前,该检测方法是唯一可用于全基因组范围、无偏倚筛选基因与小分子相互作用的综合性方法。该方法尤为突出的一个特点是无需事先了解目标基因或小分子的任何信息。尽管此类检测方法具有广泛的适用性和强大的功能,但其在其他实验室的推广应用在一定程度上受到初始资金投入以及结果分析所需生物信息学工具的限制。随着下一代测序读出技术的采用以及强大分析工具的配套使用,我们预期该方法的应用将日益广泛。
未声明任何利益冲突。
感谢多伦多大学 HIP HOP 实验室的 Ron Davis、Adam Deutschbauer 及全体成员提供的讨论与建议。C.N. 的研究得到了美国国家人类基因组研究中心(资助号 HG000205)、RO1 HG003317、加拿大卫生研究院(CIHR MOP-84305)以及加拿大癌症学会(#020380)的资助支持。J.O. 得到了美国国家人类基因组研究中心资助的斯坦福基因组培训项目(资助号 T32 HG00044)以及美国国立卫生研究院(资助号 P01 GH000205)的支持。GG 的研究由 NHGRI RO1 HG003317 和加拿大癌症学会(资助号 #020380)提供支持。TD 和 Donnelly 测序中心部分得到了加拿大创新基金会向 Brenda Andrews 博士和 Jack Greenblatt 博士提供的资助支持。A.M.S. 得到了多伦多大学开放奖学金的支持。