本方案描述了从ex vivo或in vitro细胞培养到转录组数据预处理的工作流程,用于低成本的基于转录组的药物筛选。
本方案描述了从ex vivo或in vitro细胞培养到转录组数据预处理的工作流程,用于低成本的基于转录组的药物筛选。
转录组学能够 全面深入地揭示细胞程序及其对扰动的响应。尽管在过去十年中,文库构建和测序的成本显著降低,但将这些技术应用于药物筛选所需的规模时,其费用仍然过高,限制了这些方法巨大潜力的发挥。本研究提出了一种基于转录组的低成本药物筛选系统,结合了微型化扰动培养与微型批量转录组学(mini-bulk transcriptomics)技术。优化后的微型批量转录组学方案可在经济的测序深度下获得具有生物学意义的信号,从而实现对已知药物和新分子的大规模筛选。根据所选处理方式和孵育时间,该方案可在约2天内完成测序文库的制备。由于本方案中设置了多个可暂停的步骤,文库构建和测序可在不同时间点独立进行。该方法支持高通量样本的并行处理,已成功测试同时处理多达384个样本,且未出现数据质量下降。尽管不同药物的最佳孵育时间可能存在差异,但该方法对处理条件和/或药物数量尚无已知限制。
新药研发是一个复杂且耗时的过程,涉及潜在药物及其靶点的识别、候选药物的优化与合成,以及在临床前和临床试验中对其有效性和安全性的测试1。传统的药物筛选方法,即对候选化合物库进行系统性评估以寻找治疗用途的方法,通常采用动物模型或基于细胞的检测手段来测试化合物对特定靶点或通路的作用。尽管这些方法在发现候选药物方面已取得成功,但往往无法充分揭示药物疗效背后的复杂分子机制,同时也难以阐明潜在的毒性及副作用机制。
评估全基因组范围的转录状态为克服当前药物筛选的局限性提供了一种强有力的方法,因为它能够全面评估药物处理后基因表达的变化2。通过以全基因组方式测量特定时间点表达的RNA转录本,转录组学旨在提供药物诱导的转录变化的整体视图,包括基因表达模式的改变、可变剪接以及非编码RNA的表达变化3。这些信息可用于确定药物靶点、预测药物疗效和毒性,并优化药物剂量和治疗方案。
将转录组学与无偏倚药物筛选相结合的主要优势之一,是有可能发现此前未被考虑的新型药物靶点。传统的药物筛选方法通常集中于已确立的靶分子或信号通路,这限制了新靶点的发现,并可能导致药物产生意想不到的副作用以及疗效受限。转录组学可通过揭示药物处理后发生的分子变化,克服这些局限性,从而发现可能此前未被关注的潜在靶点或信号通路2。
除了鉴定新的药物靶点外,转录组学还可用于预测药物的有效性和毒性。通过分析与药物反应相关的基因表达模式,可以开发出生物标志物,用于预测患者对特定药物或治疗方案的反应。这也有助于优化药物剂量,降低发生不良副作用的风险4。
尽管转录组学具有潜在优势,但其成本仍然是限制其在药物筛选中广泛应用的重要障碍。转录组学分析需要专用设备、专业技术以及数据分析能力,这使得资金有限的小型研究团队或机构难以在药物筛选中应用转录组学。然而,转录组学的成本正在持续下降,使其对科研群体的可及性不断提高。此外,技术进步和数据分析方法的改进使转录组学变得更加高效且经济,进一步提升了其可及性2。
在本方案中,我们描述了一种基于转录组的高维度、探索性药物筛选系统,该系统结合了微型化扰动培养与微量批量转录组学分析。5,6通过本方案,可将每个样品的成本降低至原来的1/6th 相较于当前商业化的全长mRNA测序解决方案,本方案成本较低。该流程仅需常规实验室设备,唯一例外是需使用短读长测序技术,若机构内部无测序仪器,可将测序步骤外包。优化后的微量批量测序方案可在经济有效的测序深度下提供富含生物学信息的信号,从而支持对已知药物及新分子进行广泛筛选。
本实验的目的是在不同的生物学背景下筛选药物对PBMCs的活性。该方案可应用于任何需要以转录组学读数来测试多种药物的生物学问题,从而全面揭示治疗处理对细胞的转录组水平影响。
本方案遵循波恩大学当地伦理委员会的指导原则。
1. 缓冲液、溶液和仪器的准备
2. 细胞操作
注意:人血中分离的外周血单个核细胞(PBMC)的详细冻存方案可参见文献7。

3. 测序文库制备
4. 测序与数据预处理

按照已报道的方案,将人外周血单个核细胞(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)。

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

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

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

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

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

图6:探索性数据分析。 探索性数据分析的代表性结果,显示所选免疫调节药物对PBMCs的影响。(A)按基因类型(X轴)分类的检测到的基因数量(Y轴)柱状图。(B)根据不同药物处理着色的数据集中所有基因的PCA分析。相同颜色的点为生物学重复。 请点击此处查看该图的放大版本。
| 试剂 | 浓度 | 体积 [µL] /反应 |
| 盐酸胍 | 80 mM | 7.50 |
| 脱氧核苷三磷酸混合物 (dNTPs) | 每种 10 mM | 6.52 |
| SMART dT30VN 引物 | 100 µM | 0.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 °C | 2 min |
| 冰上放置 | 2 min | |
| 逆转录 | ||
| 步骤 | 温度 | 时间 |
| 逆转录 | 42 °C | 90 min |
| 酶失活 | 70 °C | 15 min |
| 4 °C | 保持 | |
表3:mRNA变性与逆转录(RT)的PCR仪程序。
| 试剂 | 体积 [µL] /反应 |
| 高保真DNA聚合酶 | 12.50 |
| ISPCR引物(10 µM) | 0.15 |
| 无核酸酶水 | 2.35 |
| 总体积 | 15.00 |
表4:预扩增混合液。
| 步骤 | 温度 | 持续时间 | |
| 初始变性 | 98 °C | 3 min | |
| 变性 | 98 °C | 20 s | 16 – 18 个循环 |
| 退火 | 67 °C | 20 s | |
| 延伸 | 72 °C | 6 min | |
| 4 °C | 保持 |
表5:预扩增PCR仪程序。
| 步骤 | 温度 | 持续时间 |
| 片段化与接头连接(Tagmentation) | 55 °C | 8 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 °C | 5 min | |
| 初始变性 | 98 °C | 30 s | |
| 变性 | 98 °C | 10 s | 16 个循环 |
| 退火 | 60 °C | 30 s | |
| 延伸 | 72 °C | 1 min | |
| 最终延伸 | 72 °C | 5 min | |
| 4 °C | 保持 |
表8:富集PCR扩增仪程序。
| 仪器 | 上样浓度 |
| MiSeq v2 | 10 pM |
| NextSeq 500/550 | 1.4 pM |
| NovaSeq 6000 | 1250 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 scientific | 10203001 | |
| 粘性 PCR 板密封膜 | Thermo Fisher Scientific | AB0558 | |
| 扩增子片段化混合液(ATM) | Illumina | FC-131-1096 | Nextera XT DNA 文库制备试剂盒(96 样本) |
| AMPure XP 磁珠 | Beckman Coulter | A 63881 | |
| Betaine | Sigma-Aldrich | 61962 | |
| 细胞培养级 96 孔板 | Thermo Fisher Scientific | 260860 | |
| 细胞培养真空泵(VACUSAFE) | Integra Bioscience | 158300 | |
| 脱氧核苷三磷酸(dNTPs)混合液(每种 10 mM) | Fermentas | R0192 | |
| DMSO | Sigma-Aldrich | 276855 | |
| DTT(100 mM) | Invitrogen | 18064-014 | |
| EDTA | Sigma-Aldrich | 798681 | 用于贴壁细胞 |
| 乙醇 | Sigma-Aldrich | 51976 | |
| 胎牛血清 | Thermo Fisher Scientific | 26140079 | |
| 滤芯吸头(10 µL) | Gilson | F171203 | |
| 滤芯吸头(100 µL) | Gilson | F171403 | |
| 滤芯吸头(20 µL) | Gilson | F171303 | |
| 滤芯吸头(200 µL) | Gilson | F171503 | |
| 盐酸胍 | Sigma-Aldrich | G3272 | |
| ISPCR 引物(10 µM) | Biomers.net GmbH | SP10006 | 5′-AAGCAGTGGTATCAACGCAGAG T-3′ |
| KAPA HiFi HotStart ReadyMix(2X) | KAPA Biosystems | KK2601 | |
| 氯化镁(MgCl2) | Sigma-Aldrich | M8266 | |
| 96 孔磁力架 | Ambion | AM10027 | |
| 中和片段化(NT)缓冲液 | Illumina | FC-131-1096 | Nextera XT DNA 文库制备试剂盒(96 样本),或替代使用 0.2% SDS |
| 与 Nextera 兼容的索引引物 | Illumina | ||
| 无核酸酶水 | Invitrogen | 10977049 | |
| PBS | Thermo Fisher Scientific | AM9624 | |
| PCR 96 孔板 | Thermo Fisher Scientific | AB0600 | |
| PCR 板封膜仪 | Thermo Fisher Scientific | HSF0031 | |
| 青霉素/链霉素 | Thermo Fisher Scientific | 15070063 | |
| Qubit 4 荧光计 | Invitrogen | 15723679 | |
| 重组 RNase 抑制剂(40 U/µL) | TAKARA | 2313A | |
| RPMI-1640 细胞培养基 | Gibco | 61870036 | 若非使用 PBMCs,请根据细胞类型调整 |
| SMART dT30VN 引物 | Sigma-Aldrich | 5' Bio-AAGCAGTGGTATCAACGCAGAG TACT30VN-3 | |
| 常规实验室设备 | 多种 | 多种 | 例如离心机、制冰机、冰桶、蒸馏水、水浴锅 |
| SuperScript II 逆转录酶(SSRT II) | Thermo Fisher Scientific | 18064-014 | |
| SuperScript II 逆转录酶(SSRT II)缓冲液(5x) | Thermo Fisher Scientific | 18064-014 | |
| 片段化 DNA 缓冲液(TD) | Illumina | FC-131-1096 | Nextera XT DNA 文库制备试剂盒(96 样本) |
| TapeStation 4200 系统 | Agilent | G2991BA | |
| 热循环仪(S1000) | Bio-Rad | 1852148 | |
| TSO-LNA(100 µM) | Eurogentec | 5' Biotin AAGCAGTGGTATCAACGCAGAG TACAT(G)(G){G | |
| Vortex-Genie 2 混匀仪 | Sigma-Aldrich | Z258415 |