方法文章

基于转录组学的低成本药物筛选

2.1K 次观看

DOI:

10.3791/65930

2024年2月23日

本文内容

摘要

本方案描述了从ex vivoin vitro细胞培养到转录组数据预处理的工作流程,用于低成本的基于转录组的药物筛选。

摘要

转录组学能够 全面深入地揭示细胞程序及其对扰动的响应。尽管在过去十年中,文库构建和测序的成本显著降低,但将这些技术应用于药物筛选所需的规模时,其费用仍然过高,限制了这些方法巨大潜力的发挥。本研究提出了一种基于转录组的低成本药物筛选系统,结合了微型化扰动培养与微型批量转录组学(mini-bulk transcriptomics)技术。优化后的微型批量转录组学方案可在经济的测序深度下获得具有生物学意义的信号,从而实现对已知药物和新分子的大规模筛选。根据所选处理方式和孵育时间,该方案可在约2天内完成测序文库的制备。由于本方案中设置了多个可暂停的步骤,文库构建和测序可在不同时间点独立进行。该方法支持高通量样本的并行处理,已成功测试同时处理多达384个样本,且未出现数据质量下降。尽管不同药物的最佳孵育时间可能存在差异,但该方法对处理条件和/或药物数量尚无已知限制。

引言

新药研发是一个复杂且耗时的过程,涉及潜在药物及其靶点的识别、候选药物的优化与合成,以及在临床前和临床试验中对其有效性和安全性的测试1。传统的药物筛选方法,即对候选化合物库进行系统性评估以寻找治疗用途的方法,通常采用动物模型或基于细胞的检测手段来测试化合物对特定靶点或通路的作用。尽管这些方法在发现候选药物方面已取得成功,但往往无法充分揭示药物疗效背后的复杂分子机制,同时也难以阐明潜在的毒性及副作用机制。

评估全基因组范围的转录状态为克服当前药物筛选的局限性提供了一种强有力的方法,因为它能够全面评估药物处理后基因表达的变化2。通过以全基因组方式测量特定时间点表达的RNA转录本,转录组学旨在提供药物诱导的转录变化的整体视图,包括基因表达模式的改变、可变剪接以及非编码RNA的表达变化3。这些信息可用于确定药物靶点、预测药物疗效和毒性,并优化药物剂量和治疗方案。

将转录组学与无偏倚药物筛选相结合的主要优势之一,是有可能发现此前未被考虑的新型药物靶点。传统的药物筛选方法通常集中于已确立的靶分子或信号通路,这限制了新靶点的发现,并可能导致药物产生意想不到的副作用以及疗效受限。转录组学可通过揭示药物处理后发生的分子变化,克服这些局限性,从而发现可能此前未被关注的潜在靶点或信号通路2

除了鉴定新的药物靶点外,转录组学还可用于预测药物的有效性和毒性。通过分析与药物反应相关的基因表达模式,可以开发出生物标志物,用于预测患者对特定药物或治疗方案的反应。这也有助于优化药物剂量,降低发生不良副作用的风险4

尽管转录组学具有潜在优势,但其成本仍然是限制其在药物筛选中广泛应用的重要障碍。转录组学分析需要专用设备、专业技术以及数据分析能力,这使得资金有限的小型研究团队或机构难以在药物筛选中应用转录组学。然而,转录组学的成本正在持续下降,使其对科研群体的可及性不断提高。此外,技术进步和数据分析方法的改进使转录组学变得更加高效且经济,进一步提升了其可及性2

在本方案中,我们描述了一种基于转录组的高维度、探索性药物筛选系统,该系统结合了微型化扰动培养与微量批量转录组学分析。5,6通过本方案,可将每个样品的成本降低至原来的1/6th 相较于当前商业化的全长mRNA测序解决方案,本方案成本较低。该流程仅需常规实验室设备,唯一例外是需使用短读长测序技术,若机构内部无测序仪器,可将测序步骤外包。优化后的微量批量测序方案可在经济有效的测序深度下提供富含生物学信息的信号,从而支持对已知药物及新分子进行广泛筛选。

本实验的目的是在不同的生物学背景下筛选药物对PBMCs的活性。该方案可应用于任何需要以转录组学读数来测试多种药物的生物学问题,从而全面揭示治疗处理对细胞的转录组水平影响。

方案

本方案遵循波恩大学当地伦理委员会的指导原则。

1. 缓冲液、溶液和仪器的准备

  1. 准备 材料表 中所述的溶液并收集所需材料。
  2. 将水浴加热至 37 °C,并预热完全培养基(RPMI-1640 + 10% 胎牛血清 (FCS) + 1% 青霉素/链霉素)。
  3. 细胞收获时使用冰冻磷酸盐缓冲液 (PBS)。
    ​注:操作细胞时应保持环境洁净,以维持无菌状态。

2. 细胞操作

注意:人血中分离的外周血单个核细胞(PBMC)的详细冻存方案可参见文献7

  1. 细胞解冻与计数
    1. 从液氮中取出冻存管,在水浴中解冻 37 °C 轻柔倒置 2 - 3 分钟。
    2. 将解冻的细胞转移至50 mL离心管中。
    3. 用1 mL温热的完全培养基冲洗冻存管,并将此溶液逐滴加入试管中的细胞(1st 稀释步骤
    4. 以1:1比例重复稀释,直至总体积达到32 mL(共进行5次稀释,分别加入1、2、4、8和16 mL预热的完全培养基)。添加培养基时应逐滴加入,同时轻轻振荡锥形管,以减少对细胞的扰动。
    5. 将细胞悬液离心5分钟(300 x g, 20 °C)并用一个流畅的动作轻轻倾斜锥形管以去除上清液。
    6. 将细胞重悬于3 mL预热的完全培养基中,随后进行细胞计数。
    7. 用于计数时,混合 10 µL 用台盼蓝对细胞悬液进行染色(根据细胞悬液密度按1:2至1:10比例稀释),以在计数时区分活细胞与死细胞。死细胞或受损细胞因摄取染料而呈蓝色,活细胞则不被染色。
      注意:台盼蓝具有轻微细胞毒性,染色后的细胞存放时间不应超过5分钟。
    8. 使用自动细胞计数器或计数板对细胞进行计数 10 µL 染色细胞溶液
      1. 使用改进型血细胞计数板(Neubauer-improved cell chamber)时,计数四个角的大方格内的所有细胞,并利用以下公式计算细胞悬液的浓度。
        静力平衡方程,细胞计数公式;稀释与浓度分析的计算。
    9. 将细胞稀释至终浓度为 1 x 106 细胞/mL,使用预热的完全培养基。仅当所需体积低于3 mL时进行离心,并将细胞沉淀重悬至所需体积。
  2. 细胞接种与处理
    1. 种子 100 µL 每孔加入96孔细胞培养板中的细胞悬液(1 x 10⁵ 个细胞)5 细胞/孔
      注意:接种细胞数量已针对 PBMC 培养进行优化。细胞浓度过低将无法获得足量 RNA 用于测序,而细胞数量过多则会增加细胞裂解不充分及逆转录反应抑制的风险。
    2. 按治疗浓度的两倍(2x)配制药物稀释液。根据化合物的溶解度选择溶剂,优先选用完全培养基。
      注意:如果无法达到足够的溶解度,可使用磷酸盐缓冲液(PBS)或二甲基亚砜(DMSO)。最终孵育体系中 DMSO 的最大浓度不应超过 0.5 %,以防止溶剂引起的细胞毒性。
    3. 添加 100 µL 2倍药物稀释液(孔内终浓度为1X),于 37 °C 指定的时间。
      注意:应进行补充性研究,以确定最佳药物孵育时间,并排除潜在的细胞毒性机制。
  3. 细胞收获与裂解
    1. 将细胞培养板离心 10 分钟(300 x g, 20 °C). 用真空泵小心吸去上清液,同时将培养板保持倾斜角度(30°-45°将吸头对准孔板的底部角落,尽量减少细胞的移除。
    2. 加入洗涤液洗涤细胞 200 µL 加入预冷的 PBS 每孔,并以 300 × g 离心 5 分钟 g, 4 °C). 用真空泵吸去上清液,再次将培养板保持在较陡的角度,以尽可能吸除更多的 PBS。
    3. 按照所述方法制备裂解缓冲液 表1 根据所需反应(rxn)数量,额外增加10%的余量。
    4. 添加 15 µL 加入裂解缓冲液至每孔,并用粘性封板膜密封培养板,以防止样品污染。离心前先涡旋震荡培养板1分钟(1000 x g, 4 °C)并在冰上孵育5分钟。
    5. 收集 6 µL 在PCR板中加入裂解后的细胞溶液,并将细胞裂解液短暂冷冻于-80 °C 以确保最佳的细胞裂解效果。务必避免对培养板进行多次冻融循环。
      注意:暂停点:如果后续不立即进行 cDNA 的合成,细胞裂解液可保存在 -80 °C。80 °C 长达数月而不会显著降低RNA质量。

3. 测序文库制备

  1. 逆转录(RT)反应
    注意:处理 RNA 时,应将所有样品置于冰上,并确保使用无核酸酶的设备(无菌、一次性塑料耗材)和水。
    1. 按照所述方法制备逆转录反应混合液 表2 根据所需反应数量配制,额外增加10%余量。轻轻涡旋混匀,短暂离心。将混合液置于冰上保存,直至使用。
    2. 将细胞裂解液在室温(RT)下解冻,并短暂离心,使所有细胞裂解液聚集于孔板底部。
    3. 在PCR仪中按照以下步骤进行mRNA变性处理 表3.
    4. 从PCR仪中取出反应板,短暂离心以收集可能的冷凝液。加入 6 µL 将RT反应混合液加入各孔,使终体积达到 12 µL. 密封培养板以保护样品并防止蒸发、涡旋和离心。
    5. 将培养板置于热循环仪上,并按照逆转录反应程序开始运行 表3.
  2. 预扩增
    1. 解冻高保真DNA聚合酶以及 原位 PCR(ISPCR)引物在室温下。
    2. 按照说明制备预扩增混合液 表4 根据所需反应数,额外增加10%余量。短暂涡旋混匀后离心。
      注意:该酶混合液在室温下数小时内保持稳定。
    3. 离心PCR板,加入 15 µL 将预扩增混合液加入每孔,并密封板子。
    4. 将其置于热循环仪上,并按照预扩增反应程序开始运行 表5.
      1. 根据 RNA 含量调整循环次数。起始时采用较低的循环数,若 cDNA 产量不足,则逐步增加。
        注意:根据我们的经验,每种细胞类型、实验条件和处理方式的最佳循环数均应单独确定,处理因素不会影响最佳循环数。因此,建议在将本方案用于新细胞类型时,通过实验确定最佳循环数。通常情况下,原代细胞所需的循环次数较细胞系更多。暂停点:预扩增产物可保存于 - 20 °C.
  3. cDNA 纯化与质量控制(QC)
    注意:由于样品在此阶段较为稳定,可依次对每块板进行样品纯化。增加样品数量将延长实验流程的总时长,但不会限制可同时处理的样品数量。
    1. 开始前,将磁性纯化 beads 置于室温,并以高速涡旋振荡 1 min,充分重悬 beads。
    2. 加入0.8倍体积(v/v)的20 µL将磁性纯化珠加入各孔中,室温孵育5分钟。
    3. 将PCR板置于磁力架上静置5分钟,直至磁珠完全分离。小心移除上清液。
    4. 将培养板置于磁力架上,轻轻加入 100 µL 加入新配制的80%乙醇洗涤磁珠,并孵育30秒。用小体积移液器吸除上清液。重复洗涤一次。务必尽可能彻底去除乙醇。
    5. 将磁珠置于磁力架上,室温下晾干,最长不超过5分钟,直至乙醇完全挥发且磁珠不再呈现光泽。
    6. 从磁力架上取下板,重悬珠子 20 µL 加入无核酸酶水并孵育2分钟。
    7. 将培养板再次置于磁力架上,直至磁珠完全分离。
    8. 将洗脱液回收至新的PCR板中,用于cDNA质量检测和片段化。
    9. 使用 TapeStation 或 FragmentAnalyser 检测评估 cDNA 文库的片段大小分布和浓度(推荐:TapeStation D5000 检测)。具体操作详见试剂盒说明书。预期典型产量为 20 ng。
  4. 使用试剂盒进行标签化片段化
    注意:如果已有 established 的其他 Tagmentation 方案,可使用。
    1. 将 cDNA 稀释至终浓度为 150 - 300 pg/μLµL 使用无核酸酶水。
    2. 按照以下设置预先编程热循环仪 表6 以确保在将 cDNA 加入酶混合液后立即启动标签化反应。
    3. 按照说明制备标签化混合液 表7 根据所需反应数,额外增加10%的余量。
    4. 分配 3 µL 每反应将转座混合物加入新的PCR板中,并添加 1 µL 每个反应中加入 cDNA。
    5. 在加入 cDNA 后立即启动标签化反应。
    6. 通过添加使反应失活 1 µL 中和片段化缓冲液(NT;也可使用0.2%十二烷基硫酸钠(SDS))。
  5. 富集PCR
    1. 配制富集PCR反应液 表7 根据反应数量,额外增加10%余量。(推荐:使用Nextera UDI试剂盒,以防止在有序流动池(如Illumina NovaSeq 6000)上发生索引跳跃)
    2. 添加 9 µL 将富集PCR混合液加入各反应体系,使总体积达到 14 µL.
    3. 按照说明书运行富集PCR程序 表8.
      注意:对于扩增PCR,标准循环数为16个。如果cDNA质量较差,可增加额外循环数。
  6. 纯化与质量控制
    1. 实验开始前,将磁性纯化 beads 置于室温,并以高速涡旋震荡 1 分钟,充分重悬 beads。
    2. 加入1.0倍体积(v/v)14 µL)磁性纯化珠,室温孵育5分钟。
    3. 将培养板置于磁力架上,静置5 min,直至磁珠完全分离。小心移除上清液。
    4. 将培养板保持在磁力架上,轻柔地加入 100 µL 加入新配制的80%乙醇洗涤珠子,孵育30秒。使用小体积移液器吸除上清液。再重复洗涤一次。务必尽可能彻底去除乙醇。
    5. 室温下风干珠子3分钟,或直至其表面不再有光泽。
    6. 从磁力架上取下反应板,重悬磁珠 20 µL 加入无核酸酶水并孵育2分钟。
    7. 将平板再次置于磁力架上,直至磁珠分离,收集洗脱液至新的PCR板中,用于文库质控。
    8. 使用 TapeStation 或 Fragment Analyzer 检测以评估 cDNA 文库的片段大小分布和浓度(推荐:TapeStation D1000 高灵敏度检测)。具体操作详见仪器制造商说明书。通常情况下,预期产量约为 10 ng。
      注意:暂停点:PCR 产物可在 - 20 °C.

4. 测序与数据预处理

  1. 测序
    注意:以下指南适用于所有用于短读长测序的Illumina仪器。如果无法获得相关仪器,可将测序工作委托给外部测序机构完成。也可采用其他测序方法。为简化说明,本文仅报告目前使用最广泛的测序技术。
    注意:以下与软件操作相关的步骤描述基于Illumina NovaSeq6000测序仪的操作流程。
    1. 根据步骤3.6.8获得的结果,将带有唯一索引的文库按等摩尔比例混合。
    2. 使用高灵敏度检测方法测定最终文库混合物的浓度,以计算样品的摩尔浓度,公式如下:
      DNA定量公式,核酸浓度计算方程,包含ng/µL,示意图。
    3. 根据仪器说明书及实验优化条件加载流动池。常见仪器的加载浓度示例见表9
    4. 触摸屏幕选择测序以启动运行设置。
    5. 按照屏幕提示加载流动池、边合成边测序试剂盒、簇生成试剂盒、缓冲液试剂盒,并确保废液容器为空。
    6. 当仪器识别所有试剂后,点击运行设置。在此处定义运行名称和用于存储数据的输出文件夹。
    7. 设定测序参数为双端测序,两端读长均为51 bp,同时对两个索引序列进行测序,每个索引读长为8 bp。
    8. 点击查看,确认所有测序参数无误后,点击启动运行
      注意:推荐的测序深度为每个样品5 × 106条读段,最低不少于每个样品1 × 106条读段。
  2. 数据预处理
    1. 使用Bcl2Fastq2工具将原始测序数据转换为FASTQ格式,并根据样品索引进行解复用。解复用时采用默认参数设置。有关Bcl2Fastq2的详细操作说明,请参阅其参考手册。
      注意:若测序不在本实验室进行,FASTQ格式转换和解复用通常由测序机构完成。测序机构通常会提供已完成解复用的FASTQ文件供后续分析使用。
    2. 目前有多种可用于测序读段比对和丰度定量的分析方法(推荐:nf-core RNA-seq流程(https://nf-co.re/rnaseq))。该流程提供多种选项,建议采用默认设置,使用STAR8进行序列比对,Salmon9进行转录本丰度定量。
    3. 注意:在后续生物信息学分析方面,存在多种可用方法。本实验方案不涵盖所有分析流程(推荐:DESeq2分析流程10)。我们开发的基于DESeq2工作流的标准分析脚本可在GitHub上获取(https://github.com/jsschrepping/RNA-DESeq2)。

结果

按照已报道的方案,将人外周血单个核细胞(PBMCs)接种后,用不同的免疫调节药物处理,并在不同孵育时间后收集细胞,采用测序方案进行大规模转录组分析(图1)。

应在本方案之前,借助互补的实验策略并根据具体的科学问题,确定待测化合物的理想药物浓度和孵育时间。在大多数情况下,孵育2 - 4小时和24小时可分别代表药物处理后早期和晚期的转录响应情况。

评估本实验方案正确执行的最重要结果是cDNA和文库的质控(QC)图2图3)。cDNA 谱应具有较宽的分布,平均大小为 > 1000 bp (图2),较低的平均分子量或低分子量下分子的累积图4)可能表明RNA投入量较低或RNA降解。

构建高质量的测序文库也是本实验方案中的关键步骤;片段化后的文库片段长度应主要集中分布在约 250 bp 左右(图 3);过长的片段在测序过程中表现较差(图 5)。

测序完成后,原始的 FASTQ 文件将比对至相应的参考基因组(例如人或小鼠),并对每个样本的转录本丰度进行定量(参见 nf-core RNA-seq 分析流程(https://nf-co.re/rnaseq))。接下来将进行探索性数据分析,以评估数据的整体质量。比对后的数据应能检测到大量蛋白编码基因,因为本实验方案仅捕获带有聚腺苷酸尾的 RNA(图 6A)。在人源样本中,预期可检测到 15000 至 20000 个转录本(该数值基于 GENCODE 27 版人类参考基因组注释11)。进一步的探索性数据分析可包括主成分分析(PCA)。在此类分析中,数据的内在结构可通过二维图进行可视化展示。图 6B 展示了一个典型的 PCA 图,其中各点按处理条件着色,显示出三种不同实验条件形成各自的转录组谱聚类。同时可见,生物学重复样本(同色的点)在转录水平上具有高度相似性,表明本实验方案具有良好的稳健性。本图所示结果由基于 DESeq2 分析流程的 GitHub 公开工具生成(https://github.com/jsschrepping/RNA-DESeq2)。

细胞处理至测序示意图,逆转录PCR,cDNA制备,工作流程步骤,实验室设置。
图1:时间估算与工作流程。A)本实验方案在96个样本运行情况下的时间估算。(B)从细胞解冻到数据预处理各步骤的流程示意图。请按照箭头方向从上至下阅读该图。带圆圈的箭头表示步骤的重复。红色圆点表示本实验方案中进一步说明的暂停点。请点击此处查看该图的放大版本。

电泳结果图;DNA片段大小分布;大小范围15至10000 bp。
图2cDNA文库的典型结果。 微型电泳结果展示了一个代表性cDNA文库的片段大小分布。上下信号代表用于样本比对的标记物。蓝色线条指示平均片段大小(蓝色括号)。cDNA谱显示了广泛的分布,平均大小为 > 1000 bp 请点击此处以查看此图的放大版本。

凝胶电泳结果;DNA 片段大小分布;图示;峰分析;分子大小测定。
图 3:测序文库的典型结果。 微型化电泳结果显示成功构建的测序文库具有较窄的片段大小分布,平均长度约为 250 bp。上下信号代表用于样本比对的标记物。蓝色线条表示平均片段大小(蓝色括号)。请点击此处查看该图的放大版本。

电泳结果,DNA片段化图谱,显示片段大小分布,片段分析。
图 4:cDNA文库的典型非理想结果。 亚最优cDNA文库的微型电泳结果,显示其平均大小的片段分布情况 < 200 bp。上下信号代表用于样本比对的标记。蓝线表示平均片段大小(蓝色方括号)。 请点击此处以查看该图的放大版本。

电泳结果,样本强度与片段大小关系,DNA Ladder峰,用于遗传分析的图谱。
图5:测序文库中典型的非理想结果示例。 微型电泳结果显示一个非理想测序文库的片段大小分布,其中包含200 - 1000 bp的较长片段。上方和下方的信号代表用于样本校准的标记物。蓝色线条表示平均片段大小(蓝色括号)。请点击此处查看该图的放大版本。

基因表达数据分析;显示处理差异的柱状图和PCA散点图。
图6:探索性数据分析。 探索性数据分析的代表性结果,显示所选免疫调节药物对PBMCs的影响。(A)按基因类型(X轴)分类的检测到的基因数量(Y轴)柱状图。(B)根据不同药物处理着色的数据集中所有基因的PCA分析。相同颜色的点为生物学重复。 请点击此处查看该图的放大版本。

试剂浓度体积 [µL] /反应
盐酸胍80 mM7.50
脱氧核苷三磷酸混合物 (dNTPs)每种 10 mM6.52
SMART dT30VN 引物100 µM0.33
无核酸酶水0.65
总体积15.00

表1:裂解缓冲液。

试剂体积 [µL] /反应
SSRT II buffer (5x)2.00
DTT (100 mM)0.50
Betaine (5 M)2.00
MgCl2 (1 M)0.14
SSRT II (200 U/µL)0.25
RNA酶抑制剂 (40 U/µL)0.25
TSO-LNA (100 µM)0.20
无核酸酶水0.66
总体积6.00

表2:逆转录(RT)反应体系。

mRNA 变性
步骤温度时间
mRNA 变性95 °C2 min
冰上放置2 min
逆转录
步骤温度时间
逆转录42 °C 90 min
酶失活70 °C15 min
4 °C保持

表3:mRNA变性与逆转录(RT)的PCR仪程序。

试剂体积 [µL] /反应
高保真DNA聚合酶12.50
ISPCR引物(10 µM)0.15
无核酸酶水2.35
总体积15.00

表4:预扩增混合液。

步骤温度持续时间
初始变性98 °C3 min
变性98 °C20 s16 – 18 个循环
退火67 °C20 s
延伸72 °C 6 min
4 °C保持

表5:预扩增PCR仪程序。

步骤温度持续时间
片段化与接头连接(Tagmentation)55 °C8 min
4 °C保持

表6:打断片段化反应热循环程序。

片段化反应混合液
试剂体积 [µL] /反应
Amplicon Tagment Mix (ATM)1.0
Tagment DNA Buffer (TD)2.0
总体积3.0
扩增富集 PCR 混合液
试剂体积 [µL] /反应
高保真 DNA 聚合酶7.0
与 Nextera 兼容的索引引物2.0
总体积9.0

表7:片段化反应混合液与富集PCR混合液。

步骤温度持续时间
热启动72 °C5 min
初始变性98 °C30 s
变性98 °C10 s16 个循环
退火60 °C30 s
延伸72 °C1 min
最终延伸72 °C5 min
4 °C保持

表8:富集PCR扩增仪程序。

仪器上样浓度
MiSeq v210 pM
NextSeq 500/5501.4 pM
NovaSeq 60001250 pM

表9:常见测序仪器上样浓度的示例。

讨论

药物发现与药物开发可从批量转录组学所提供的细胞过程整体视角中显著受益。然而,由于标准批量RNA测序方案的实验成本较高,该方法常受到限制,从而阻碍了其在学术环境中的应用以及在工业规模上的扩展潜力。

该实验方案中最关键的步骤是细胞复苏以及文库构建的初始步骤。复苏后确保细胞具有较高的存活率,对于后续处理和转录组分析的成功至关重要。从细胞收集到cDNA合成的文库构建初期步骤,对于维持RNA的完整性至关重要。在此阶段,必须始终将细胞裂解物置于冰上,并尽快处理样品。若观察到RNA出现严重降解,应使用特定产品清洁实验台面和仪器,以抑制任何RNase活性。

本文所述方案目前针对健康供体外周血单个核细胞(PBMCs)的药物处理进行了优化,最长处理时间不超过24小时。若使用不同细胞类型或延长孵育时间,可能需要相应优化接种细胞数量和培养条件。

通过本方案,我们提供了一种使用常规实验室设备对药物处理后的外周血单个核细胞(PBMCs)进行转录组分析的工作流程,且无需商业试剂盒。采用该方法并省去RNA纯化步骤,我们显著降低了成本,从而能够对大量化合物进行平行分析。

由于该方案基于体外实验,其主要局限性在于无法评估药物可能经历的任何代谢加工过程,而这些过程可能导致不同的生物活性。此外,所捕获的转录本数量和测序深度将低于标准的批量全长转录组学方法,因此这些数据不适用于需要更高信息含量的应用,例如差异剪接分析或单核苷酸多态性的定量分析。

披露

作者声明无竞争利益。

致谢

J.L.S. 受德国研究基金会(DFG)在德国卓越战略(EXC2151-390873048)以及 SCHU 950/8-1;GRK 2168,TP11;CRC SFB 1454 Metaflammation,IRTG GRK 2168,WGGC INST 216/981-1,CCU INST 217/988-1,德国联邦教育和研究部(BMBF)资助的卓越项目 Diet-Body-Brain(DietBB);以及欧盟项目 SYSCID(资助编号 733100)的支持。M.B. 受 DFG(IRTG2168-272482170,SFB1454-432325352)支持。L.B. 受 DFG(ImmuDiet BO 6228/2-1,项目编号 513977171)和德国卓越战略(EXC2151-390873048)支持。图片使用 BioRender.com 创建。

材料

本文使用的材料清单
姓名公司目录编号评论
50 mL 锥形管fisher scientific10203001
粘性 PCR 板密封膜Thermo Fisher ScientificAB0558
扩增子片段化混合液(ATM)IlluminaFC-131-1096Nextera XT DNA 文库制备试剂盒(96 样本)
AMPure XP 磁珠Beckman CoulterA 63881
Betaine Sigma-Aldrich61962
细胞培养级 96 孔板Thermo Fisher Scientific260860
细胞培养真空泵(VACUSAFE)Integra Bioscience158300
脱氧核苷三磷酸(dNTPs)混合液(每种 10 mM)FermentasR0192
DMSOSigma-Aldrich276855
DTT(100 mM)Invitrogen18064-014
EDTASigma-Aldrich798681用于贴壁细胞
乙醇Sigma-Aldrich51976
胎牛血清Thermo Fisher Scientific26140079
滤芯吸头(10 µL)Gilson F171203
滤芯吸头(100 µL)Gilson F171403
滤芯吸头(20 µL)Gilson F171303
滤芯吸头(200 µL)Gilson F171503
盐酸胍Sigma-AldrichG3272
ISPCR 引物(10 µM)Biomers.net GmbHSP100065′-AAGCAGTGGTATCAACGCAGAG
T-3′
KAPA HiFi HotStart ReadyMix(2X)KAPA BiosystemsKK2601
氯化镁(MgCl2) Sigma-AldrichM8266
96 孔磁力架AmbionAM10027
中和片段化(NT)缓冲液 IlluminaFC-131-1096Nextera XT DNA 文库制备试剂盒(96 样本),或替代使用 0.2% SDS
与 Nextera 兼容的索引引物Illumina
无核酸酶水Invitrogen10977049
PBSThermo Fisher ScientificAM9624
PCR 96 孔板Thermo Fisher ScientificAB0600
PCR 板封膜仪Thermo Fisher ScientificHSF0031
青霉素/链霉素 Thermo Fisher Scientific15070063
Qubit 4 荧光计Invitrogen15723679
重组 RNase 抑制剂(40 U/µL)TAKARA2313A
RPMI-1640 细胞培养基 Gibco61870036若非使用 PBMCs,请根据细胞类型调整 
SMART dT30VN 引物Sigma-Aldrich5' Bio-AAGCAGTGGTATCAACGCAGAG
TACT30VN-3
常规实验室设备多种多种例如离心机、制冰机、冰桶、蒸馏水、水浴锅
SuperScript II 逆转录酶(SSRT II)Thermo Fisher Scientific18064-014
SuperScript II 逆转录酶(SSRT II)缓冲液(5x)Thermo Fisher Scientific18064-014
片段化 DNA 缓冲液(TD)IlluminaFC-131-1096Nextera XT DNA 文库制备试剂盒(96 样本)
TapeStation 4200 系统AgilentG2991BA
热循环仪(S1000)Bio-Rad1852148
TSO-LNA(100 µM)Eurogentec5' Biotin AAGCAGTGGTATCAACGCAGAG
TACAT(G)(G){G
Vortex-Genie 2 混匀仪Sigma-AldrichZ258415

参考文献

  1. Hughes, J. P., Rees, S., Kalindjian, S. B., Philpott, K. L. Principles of early drug discovery. Br J Pharmacol. 162 (6), 1239-1249 (2011).
  2. Yang, X., et al. High-throughput transcriptome profiling in drug and biomarker discovery. Front Genet. 11, 19(2020).
  3. Bonaguro, L., et al. A guide to systems-level immunomics. Nat Immunol. 23 (10), 1412-1423 (2022).
  4. Carraro, C., et al. Decoding mechanism of action and sensitivity to drug candidates from integrated transcriptome and chromatin state. ELife. 11, 78012(2022).
  5. Picelli, S., et al. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods. 10 (11), 1096-1098 (2013).
  6. Picelli, S., et al. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc. 9 (1), 171-181 (2014).
  7. De Domenico, E., et al. Optimized workflow for single-cell transcriptomics on infectious diseases including COVID-19. STAR Protoc. 1 (3), 100233(2020).
  8. Dobin, A., et al. ultrafast universal RNA-seq aligner. Bioinformatics. 29 (1), 15-21 (2013).
  9. Patro, R., et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 14 (4), 417-419 (2017).
  10. Love, M. I., Huber, W., Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550(2014).
  11. Frankish, A., et al. GENCODE reference annotation for the human and mouse genomes. Nucleic Acids Res. 47, D766-D773 (2019).

重印与许可

标签

PCR