2. 实验设计

做时间序列RNA-seq，最怕什么？

我最怕实验做完了，拿到数据才发现——

• 生物学重复不够，差异结果全靠运气
• 时间点选得太稀，关键拐点完美错过
• 样本混到一起，批次效应把真实信号淹了

这些坑，我全踩过。今天咱们就聊聊，怎么把实验设计这张图画好。

2.1 生物学重复：三头猪还是十头牛？

生物学重复，不是技术重复。技术重复是把同一管RNA上机两次，看的是机器噪音。生物学重复是取不同个体的样本，反映的是群体差异。

时间序列RNA-seq的重复设计，我一般遵循以下原则：

样本来源	最低重复	建议重复	原因
细胞系（同一批次培养）	3	3~4	遗传背景一致，变异小
小鼠/大鼠（近交系）	3	4~6	个体间仍有微小差异
人/其他动物（个体差异大）	5	6~10	遗传背景混杂，噪音大
临床组织（活检/尸检）	6	8~12	批次、病程、用药史混杂

时间序列还有一个特殊问题：每个时间点的个体必须是独立采样，不能在同一只动物身上反复穿刺。为什么呢？

举个例子——

你想看小鼠肝损伤后的基因表达变化，每6小时从同一只小鼠取血。结果发现某个炎症基因持续升高。但这是真正的响应，还是手术应激 + 反复抓取的累积效应？你分不清。

2.2 时间点选择：密一点还是松一点？

时间点选少了，你会错过调控峰。选多了，经费撑不住。这里有几个经验法则：

• 稳态 vs 动态过程：如果是昼夜节律这种周期性变化，4~6小时间隔足够。如果是药物响应这种快速爆发，可能每30分钟取一次点。
• 先查文献：你研究的通路，别人用过什么时间窗口？直接拿来用，再在关键区间加密1~2个点。
• 做预实验：我一般会提前做一次低深度测序（1M reads/样），跑3~5个粗时间点。看看哪个时间窗口表达波动最大，再细化。
• 考虑一下”零时刻”：很多人只设处理后的时间点。但你知道处理前是什么基线吗？没有t0，后面的差异都是相对的，你敢信？

举个例子，看一个典型的昼夜节律时间序列设计：

# 时间点设置示例（昼夜节律：小鼠）
# 开灯时间 ZT0，关灯时间 ZT12
# 每4小时采样，共6个时间点
time_points = [ZT0, ZT4, ZT8, ZT12, ZT16, ZT20]
# 每个时间点 4 只独立小鼠
replicates = 4
# 总样本数 = 6 × 4 = 24

2.3 采样策略：稳定、快速、一致

采样质量直接决定RNA降解程度。我自己总结了一套”三固定”原则：

1. 固定流程：从取组织到液氮速冻，时间控制在2分钟以内。超过5分钟的样本，我觉得可以直接扔了。
2. 固定时间：所有时间点的采样操作，最好在一天的同一时间段完成（比如都在上午9~11点）。避免昼夜节律干扰。
3. 固定人员：不同人的操作习惯不同，这可能是隐含的批次效应。我习惯同一个技术人员负责同一批实验的采样。

还有一个容易被忽略的细节：样本保存液。如果无法立即提取RNA，使用RNAlater这类稳定液。但注意——样本在RNAlater中放太久也不好，一般4°C保存不超过1个月，-80°C不超过6个月。

2.4 批次效应：你防得住吗？

批次效应，说白了就是”非生物学变异”。

按时间序列RNA-seq的特点，批次效应主要来自三个方面：

• 不同天数的建库：今天建一批，下周建一批，试剂批次、操作手差异都跑不掉。
• 不同测序lane/flowcell：同一个样本不同lane测，也会出现系统性偏差。
• 不同时间点采集时的环境变化：室温、湿度、操作疲劳程度等。

怎么防？我的策略是：

1. 随机化样本分配：不要按时间顺序上机！把所有样本拆散，随机分配到不同的建库批次和lane中。保证每个批次里各个时间点的样本比例相似。
2. 加入“桥梁样本”：用一个稳定的参考RNA（比如混合所有样本的pool），在每个批次中都加一管。后期可以通过这个样本校正批次效应。
3. 记录所有批次变量：建库日期、试剂盒货号、操作人员、测序仪器、flowcell编号等。记录下来，分析时作为协变量。
4. 尽量单批次完成：如果样本量 < 96，我强烈建议一次性完成建库和测序。多批次的坑，能避就避。

2.5 知识结构总览

2.6 检查清单与避坑指南

最后，我给你一个我自己用的检查清单。每次设计时间序列RNA-seq实验前，逐条确认：

• ☐ 每个时间点是否设了独立的生物学重复？（不得纵向采样同一只/同一孔）
• ☐ 重复数是否符合最低要求？（细胞≥3，动物≥4，临床≥6）
• ☐ 时间点是否覆盖了关键生理窗口？是否包含t0基线？
• ☐ 预实验是否已经跑过粗时间窗？
• ☐ 采样流程是否标准化？RNA降解风险是否可控？
• ☐ 所有时间点的样本是否已随机分配到建库批次？
• ☐ 是否准备了桥梁样本用于批次校正？
• ☐ 批次变量（日期、操作员、试剂盒、仪器）是否已记录在元数据表格中？