CHAPTER 03
3. 对RNA-seq筛选得到的候选分子进行独立验证和再次收缩
上一阶段利用RNA-seq完成了从全转录组到候选机制的第一次收缩:首先确认X干预能够引起稳定的整体转录响应,随后通过差异表达分析、GO/KEGG富集和GSEA等方法,将大量表达变化重新组织为与核心表型P相关的候选生物学过程和信号通路,并进一步从这些通路中筛选出Y1、Y2、Y3等少量优先候选分子。至此,研究已经逐步聚焦到“哪些分子可能参与X引起P的过程”。

不过,这些Y仍然主要来源于同一批RNA-seq样本及其生物信息学分析。它们之所以被选中,是因为在转录组数据中表现出稳定差异、位于与P相关的候选通路,并且从已有生物学知识来看具有解释P的可能性。这些条件提高了候选分子的研究优先级,但RNA-seq观察到的变化能否在新的实验中再次出现?以及这种转录水平的变化是否真正延伸到了能够发挥生物学作用的蛋白或功能状态?
因此,研究在这一阶段需要从“组学发现”转入“独立实验验证”。其目的是利用新的生物学样本逐层淘汰不稳定或缺乏功能对应关系的候选,将Y1、Y2、Y3进一步收缩为能够进入功能干预的核心候选Y。这个过程通常依次经过mRNA重复验证、蛋白或功能状态验证以及与核心表型P的再次比较,形成“RNA-seq候选分子→独立样本验证→功能层验证→候选重新排序”的第二次收缩。
3.1 用独立生物学样本确认候选分子的mRNA变化
候选分子的第一层验证通常仍然从转录水平开始。RNA-seq已经提供了Y1、Y2、Y3的表达变化信息,但这一结果首先需要回答的是“这些差异是否能够离开最初的测序样本而稳定重复”。因此,更有价值的验证是重新建立X干预组和匹配对照组,在与RNA-seq实验相同或严格对应的条件下获得新一批独立生物学样本,再检测候选分子的表达。
这里需要区分两种不同意义的“验证”。如果使用原RNA-seq样本进行qRT-PCR,主要确认的是两种测量方法对同一样本是否得到一致结果;如果使用新的独立样本,则进一步考察RNA-seq发现的生物学差异能否跨实验批次重复。RNA-seq与qPCR本身并不是天然相互矛盾的检测体系。利用全转录组RT-qPCR数据对多种RNA-seq分析流程进行比较时发现,不同RNA-seq流程获得的基因表达量与qPCR总体具有较高一致性,约85%的基因在组间Fold Change方向上能够得到一致结果。 因此,在机制筛选研究中,qRT-PCR更重要的价值是利用独立实验判断某个被选中的候选是否值得继续进入后续功能研究。有关RNA-seq结果是否需要独立验证的讨论也指出,当后续研究重点集中到少量关键候选基因时,在额外生物学样本中进行正交方法验证具有更明确的意义。
在独立验证过程中,首先需要比较变化方向。若RNA-seq提示Y1在X干预后表达升高,那么在新的实验中也应能够重复观察到同方向变化;如果独立实验中反复出现相反方向,或者多数实验中没有稳定变化,那么即使Y1在原RNA-seq数据中的Fold Change和统计显著性较高,其作为核心候选的优先级也需要下降。相反,某些候选在RNA-seq中的变化幅度并不是最大,却能在多批独立实验中持续呈现相同方向,其稳定性反而使其更适合进入后续机制验证。
qRT-PCR与RNA-seq之间也不需要机械地追求完全相同的Fold Change。两种方法在定量原理、归一化方法和动态范围等方面存在差异,因此同一基因在两种平台上得到不同的具体倍数并不构成矛盾。这里更需要关注的是方向是否一致、差异是否能够在独立样本中重复,以及组内变异是否处于可接受范围。与此同时,qPCR实验本身也需要具备可靠的样本处理、引物特异性、参考基因和数据分析体系。
经过这一轮验证,一部分只在RNA-seq发现阶段出现、但无法在独立实验中重复的候选会被淘汰,剩余候选则获得了一层新的证据:它们成为了“在X干预以后能够稳定重复发生mRNA变化的分子”。由此可以把前一阶段的候选关系进一步收缩为:X干预 → Y的mRNA稳定改变。
3.2 从mRNA变化进一步确认Y的蛋白表达或功能状态
对于蛋白编码基因而言,RNA-seq直接测量的是转录本,而多数细胞功能最终由蛋白的数量、修饰状态、亚细胞定位、相互作用或酶活性执行。因此,当候选Y通过独立mRNA验证以后,下一步需要根据Y本身的分子性质判断:X干预以后,Y是否发生了足以改变其实际生物学功能的变化。
mRNA与蛋白之间并不存在固定的一一对应关系。蛋白丰度除了受到mRNA水平影响,还受到翻译效率、蛋白降解以及其他转录后调控过程的共同控制,因此mRNA丰度只能部分解释蛋白水平的变化。 transcriptomics和proteomics提供的是相互关联但并不冗余的信息,两者不一致本身有时也能够揭示转录后调控或缓冲机制。
因此,对于主要依赖总蛋白丰度发挥作用的候选,可以通过Western blot、流式细胞术或免疫荧光等方法观察蛋白水平是否随X发生相应变化;如果Y属于分泌因子,则检测培养上清或相关生物样本中的分泌量可能比细胞内总蛋白更接近其实际功能状态。对于激酶或其他主要依赖翻译后修饰调节活性的分子,仅检测总蛋白则可能无法回答真正的问题,此时需要进一步观察磷酸化、裂解或其他活化形式。类似地,转录因子的功能可能更多体现为核转位或DNA结合活性,酶类分子的变化可能需要通过酶活性或催化产物判断,而受空间定位调节的蛋白则可能需要比较其胞核、胞质或特定细胞器中的分布。
由此,这一阶段的检测依据是“哪一个指标最能够代表Y真正发挥功能的状态”。这也解释了为什么不同候选不能采用完全相同的淘汰标准。
例如,当Y的mRNA和总蛋白均稳定向同一方向变化时,可以说明转录变化已经延伸到蛋白表达层;如果mRNA明显改变而总蛋白始终稳定,则提示转录变化可能受到翻译效率或蛋白稳定性的进一步调节,此时需要重新评估其作为核心功能节点的优先级;而如果总蛋白变化并不明显,但磷酸化、核转位或其他活化状态发生清晰改变,则真正参与P的可能是Y的功能状态而不是蛋白总量。
经过这一层验证,候选证据进一步推进到“具有潜在功能意义的分子变化”。此前建立的关系可以继续延伸为:X干预 → Y的mRNA发生稳定变化 → Y的蛋白表达或功能状态发生相应变化。
3.3 将分子变化重新放回核心表型P,确定核心候选Y
经过独立mRNA以及蛋白或功能状态验证以后,候选数量通常已经明显减少。此时筛选的重点是比较多个候选中,哪些分子的变化最稳定、最接近前期RNA-seq提示的候选通路,并且与核心表型P表现出较高的一致性。
重复性首先构成候选排序的重要基础。如果某一Y能够在RNA-seq以及多批独立实验中始终保持相同方向的变化,那么这一结果说明它与X干预之间的关系具有较好的稳定性。对于后续机制研究而言,这种跨批次稳定性通常比单次实验中极大的Fold Change更重要,因为后续的loss-of-function、gain-of-function和rescue实验都依赖一个能够可靠重复的分子节点。
在此基础上,还可以进一步比较Y与P之间的动态关系。例如,在不同X干预强度或不同时间点下同步检测X、Y和P,可以观察三者是否形成相对一致的变化趋势。当X变化较弱时Y和P均变化较弱,而随着X干预增强,Y和P也逐渐出现更明显的响应,这种剂量或时间维度上的一致性能够进一步提高Y作为核心候选的优先级。
但这种同步变化属于关联证据。即使X、Y和P在多个条件下高度同步,也存在另一种解释:X可能通过两条相互独立的路径分别影响Y和P,Y只是与P共同响应X,而并非P形成所必需的中间节点。因此,这一阶段可以利用动态一致性帮助“排序候选”。
候选排序还需要重新连接上一阶段获得的通路信息。一个Y如果不仅能够在独立实验中稳定重复,而且位于与P相关候选通路的关键调控位置,其优先级通常高于位于外围、与P缺乏明确功能联系的伴随差异基因。与此同时,还需要考虑下一阶段能否真正对它实施干预,因为机制验证最终必须改变Y本身,再观察P是否受到影响。因此,可进行可靠敲低、敲除、过表达或药理学调节的候选,更适合优先进入首轮因果实验;缺乏有效干预工具的候选可以保留研究价值,但可能并不适合作为最先推进的机制节点。
由此,候选Y的优先级实际上来自几个相互连接的证据层级:RNA-seq提供其最初的表达差异和通路背景,独立qRT-PCR确认这种变化能够跨实验重复,蛋白或功能状态检测判断转录差异是否形成了实际分子效应,而不同时间、剂量或实验批次中的变化则进一步观察Y与核心表型P是否保持稳定的一致关系。最后再结合Y在候选通路中的位置以及后续功能干预的可行性,将Y1、Y2、Y3等逐渐收缩为1—3个优先级较高的核心候选。若多个候选仍然具有相近证据,也不需要在这一阶段依据表达数据人为确定唯一分子,而可以让它们平行进入下一轮初步功能筛选。
因此,本阶段完成的是从“组学候选”到“实验候选”的第二次收缩:RNA-seq筛选Y1、Y2、Y3 → 独立样本mRNA验证 → 蛋白或功能状态验证 → 比较跨实验重复性 → 比较与P的动态一致性 → 综合通路位置及实验可行性 → 确定核心候选Y。
到这里,研究能够支持的结论是:X干预能够稳定伴随Y的表达或功能状态变化,而且这种变化能够在独立实验中重复,并与核心表型P保持较高的生物学一致性,因此Y具有进入后续机制研究的充分优先级。
CHAPTER 04
4. 验证候选分子Y是否具有功能性贡献
上一阶段已经完成了从“组学候选”到“实验候选”的第二次收缩。经过独立样本的qRT-PCR验证以及蛋白表达、活化状态或其他功能指标检测,Y不再只是RNA-seq中一个具有统计差异的基因,而成为能够在不同实验中稳定响应X、并且其变化与核心表型P保持一定一致性的核心候选分子。
但是,到这一位置为止,所有实验仍然遵循同一个观察方向:先改变X,再观察Y和P是否随之发生变化。即使X、Y和P能够在不同批次、不同时间点或不同干预强度下保持稳定的同步关系,也仍然存在多种可能的解释。Y可能位于X影响P的功能路径中,也可能只是X干预后形成的下游伴随变化;还可能存在另一条与Y无关的通路真正驱动P,而Y只是同时受到X调控。因此,前面建立的是一个稳定的候选关系,却还没有直接检验Y本身对P是否具有功能贡献。
研究逻辑在这一阶段因此需要发生一次明显转换:不再继续观察“Y是否跟着X变化”,而是主动改变Y本身,再观察P如何响应。对于基因功能研究而言,loss-of-function和gain-of-function实验的基本价值就在于通过人为扰动某个分子,观察系统表型随之发生的变化。基因功能的推断长期依赖于对基因活性的直接扰动,但不同敲低、敲除或抑制策略产生的效应会受到干预强度、持续时间以及实验体系的影响,功能结论需要建立在确认干预有效并具有足够特异性的基础上。
因此,这一阶段可以沿着三个连续的问题推进。
首先,单独改变Y时,P是否发生变化;随后判断这种变化方向是否与前面提出的X—Y—P模型相符;最后再将X和Y放在同一个实验体系中,通过联合干预和rescue判断X产生的P变化是否依赖Y。由此,前面的观察链:X干预 → Y变化的同时P变化,开始转变为真正的功能扰动链:直接改变Y → P随之变化,并进一步进入:X干预 + 反向调节Y → X引起的P变化被削弱或恢复。
4.1 单独干预Y,判断Y本身能否影响核心表型P
进入功能验证以后,第一步是把X暂时从实验设计中拿开,只针对核心候选Y进行干预。这样可以把问题简化为一个直接的输入—输出关系:人为改变Y以后,前期已经定义好的核心表型P是否发生相应改变。
具体干预方式取决于Y本身的分子性质。对于基因或非编码RNA,可以采用敲低、敲除、CRISPRi、过表达或其他遗传学方法;对于具有成熟药理学工具的受体、酶或信号蛋白,也可以采用抑制剂、激动剂或其他能够改变其功能状态的处理方式。不同方法虽然都能够形成loss-of-function或gain-of-function,但并不具有完全相同的生物学含义。例如,RNA干扰通常降低mRNA水平,基因敲除可能长期消除蛋白功能,而药理学抑制往往只影响某一活性状态。不同干预方式可能产生不同强度和持续时间的效应,因此实验解释首先需要明确“实际改变的是Y的哪一个层面”。
在观察P之前,需要先确认Y本身确实按照预定方向发生了改变。如果采用RNAi,需要确认mRNA和必要的蛋白水平确实下降;如果进行过表达,需要确认新增的Y能够形成预期蛋白或功能状态;如果研究的是激酶、转录因子或酶,则总表达量并不一定代表真正的干预效果,还需要检测磷酸化、核转位、DNA结合能力或酶活等更直接的功能指标。只有先证明Y的功能状态确实被改变,后续出现或未出现的P变化才具有明确的解释基础。
随后需要重新检测第1部分已经建立的核心表型P。这里最好尽可能沿用前面定义P时使用的核心功能实验,因为此时真正需要比较的是两个实验阶段中的同一个生物学结果:最初观察到“改变X能够改变P”,现在则检验“改变Y是否也能够改变这个P”。如果前面P由增殖、迁移、侵袭、凋亡、免疫功能、炎症状态或代谢功能定义,那么这里仍应以相应的功能实验作为主要终点。
如果Y被有效干预以后,P在多次独立实验中都没有发生相应变化,那么需要降低Y作为核心功能节点的优先级,这说明Y的表达变化可能并不是当前P形成所需要的主要功能事件。Y仍可能参与其他表型、反馈调节或与当前实验条件不同的生物过程,但它对这里所研究的P尚未表现出清晰的功能贡献。相反,如果单独改变Y能够稳定改变P,那么Y成为“人为改变以后能够改变P”的功能性候选。
在条件允许时,还可以从相反的方向操纵Y。例如,在降低Y以后观察P是否减弱,同时提高Y观察P是否增强;或者根据Y与P之间实际的负调控关系获得方向相反但逻辑对应的结果。这种双向扰动能够减少某一种干预方式自身造成的解释偏差,并使Y与P之间的功能关系更加完整。
这里常用“必要性”和“充分性”来描述两种实验逻辑。若提高Y能够在当前实验背景下诱导P发生预期变化,可以认为结果支持Y具有驱动该表型的能力,接近对“充分性”的检验;若降低Y能够削弱原本存在的P,则支持Y对该表型的形成或维持具有功能要求,接近对“必要性”的检验。它们并不等同于对绝对充分性或绝对必要性的证明,因为基因冗余、不完全敲低、剂量效应和细胞补偿都可能影响实验结果。例如,loss-of-function实验中的阴性结果有时来自残余功能或遗传补偿,而不一定意味着Y与P完全无关。
经过这一部分,可以把证据推进到:直接改变Y → 核心表型P发生相应改变。
这一步建立了Y与P之间的功能联系,但没有回答Y是否处于X的下游。Y完全可能通过一条独立于X的路径影响同一个P,因此,下一步需要把Y干预产生的表型方向重新放回前面已经提出的X—Y—P关系中进行检验。
4.2 检验Y的功能方向是否与X—Y—P模型一致
当Y已经被证明能够影响P以后,研究开始检验这个功能是否能够解释前面观察到的X干预表型,这里需要先根据已有实验结果明确一个能够被实验检验的方向模型。
例如,前期实验如果观察到:X降低 → Y升高,同时P降低
那么一个可能的模型是:X降低以后引起Y升高,而Y升高参与了P降低。
从这个模型可以事先推导出一个具体预测:如果Y确实参与这一过程,那么在不改变X的情况下单独提高Y,也应该使P向“降低”的方向移动。反过来,如果降低Y,则应使P向相反方向变化,或者至少削弱由Y升高所对应的表型。
这种“先预测方向,再进行干预”的设计与单纯观察相关性不同。相关性只能说明几个变量共同变化,而功能扰动要求模型能够预测变量被人为改变后的结果。如果模型预测“提高Y应该使P下降”,实验结果却稳定显示提高Y反而使P升高,那么前面的X—Y—P解释就出现了冲突。此时需重新考虑Y在系统中的角色:它可能属于X干预后的代偿性反馈,也可能与P形成平行响应,或者前面对Y所在通路位置的判断需要重新调整。
如果Y干预以后不仅核心功能表型P按照预测方向变化,而且与P相关的一组辅助分子指标也发生方向一致的改变,那么可以进一步提高这一功能关系的稳定性。不过,辅助分子指标仍然用于支持P,而不能替代原本的功能表型。例如,如果P定义为迁移能力改变,那么相关迁移蛋白的变化可以作为支持,但核心证据仍然来自迁移功能本身。
在条件允许时,还可以设置不同程度的Y干预,观察Y变化幅度与P之间是否形成一定的剂量—效应关系。若随着Y逐渐改变,P也逐渐增强、减弱或最终进入平台期,可以进一步说明表型响应与Y功能状态之间存在连续关系。
经过前两步以后,现在已经获得了两类彼此连接的结果:X干预后,Y与P按照特定方向共同改变;并且单独改变Y,可以按照这一模型预测的方向改变P。
Y因此已经具有较明确的功能性贡献。但此时仍存在一个尚未排除的模型:X → P,以及另一条独立的Y → P。
也就是说,X和Y都能够影响P,却不一定意味着X对P的作用经过Y传递。要区分“两个都能影响P”和“Y参与了X对P的效应传递”,需要进一步将X与Y同时纳入一个实验体系。
4.3 通过X与Y联合干预判断X产生的表型是否依赖Y
联合干预把前面的两个单变量实验连接到了一起。此时的问题是:
当X已经产生P变化以后,如果人为阻断或恢复Y,X所产生的表型还能否维持?
这一步通常被称为rescue实验。在更广义的遗传学框架中,这类双扰动实验与epistasis分析具有相似的推理结构:通过比较单独扰动与联合扰动的表型,判断两个功能节点之间是否存在层级或依赖关系。在细胞实验中利用RNAi、过表达或药物进行的联合干预能够检验Y是否参与X效应的功能传递,并为上下游关系提供功能证据。
实验分组需要根据前期已经建立的具体方向关系设计。一般至少需要Control、X单独干预、Y单独干预以及X+Y联合干预几个条件。联合组中的Y操作应当与X诱导的Y变化方向相反,其目的不是单纯再改变一次Y,而是把X已经造成的Y变化向基础状态拉回。
例如,如果前面已经建立:X降低 → Y升高 → P降低,
那么联合实验可以在X降低的基础上抑制Y,使Y从升高状态重新向基础水平下降。
如果结果形成:
Control:Y、P处于基础状态
X降低:Y升高,P降低
X降低 + Y抑制:Y被重新拉低,同时P部分或明显恢复
那么与单纯观察X、Y、P共同变化相比,这个结果提供了新的信息:X产生P降低的过程中需要Y的升高参与;当这一变化被人为阻断时,X的表型效应随之减弱。
这一点也说明为什么Y单独干预与rescue不能互相替代。Y单独干预回答的是:Y能否改变P?
而X+Y联合干预进一步回答:X引起的P变化在多大程度上依赖Y的改变?
在RNAi、CRISPR等基因功能研究中,rescue还有另一层方法学价值,即帮助判断干预产生的表型是否真正来自目标基因。RNAi研究方法学综述将使用不受siRNA作用的目标基因重新恢复表型视为验证on-target效应的重要策略;CRISPR功能研究同样常利用重新引入目标基因来提高基因型与表型因果对应关系的可信度。 但rescue实验本身也需要控制表达水平和操作特异性。过量表达Y可能产生正常生理条件下不存在的表型,甚至通过非生理性相互作用形成“表面上的恢复”,因此恢复到接近生理水平通常比极高水平过表达更容易解释。
联合干预实验中还需要同时确认X、Y和P三个层级。首先应验证X仍然保持预定的干预状态,因为如果恢复Y的操作同时意外改变了X,那么就无法判断P的恢复究竟来源于哪个节点。其次,需要直接证明Y确实已经按照设计被恢复或阻断。只有这两点均得到确认,P的变化才能用于回答X是否依赖Y。如果只观察到“P没有被救回来”,却没有证实Y是否真正被恢复,那么rescue失败既可能意味着Y并不是关键中介,也可能仅仅意味着Y的反向干预没有成功。
rescue结果也不需要机械地要求P完全恢复到Control。多数复杂细胞表型受到多个调控分支共同影响,X可能同时通过Y和其他并行机制作用于P。因此,如果恢复Y只能部分削弱X的表型,仍然可以支持Y具有功能贡献。
例如,X干预使P相对对照下降50%,而恢复Y后P只下降20%。这样的结果更适合解释为:Y介导了X对P的部分效应。
剩余的表型提示X可能还存在其他平行分支。只有当恢复Y以后P几乎完全返回基础水平时,才可以认为Y可能承担X对P效应中的主要功能部分;即便如此,也仍然不能单凭这一实验宣称Y是唯一机制。
如果实验体系允许,还可以将“模拟”和“恢复”两个方向结合。首先单独改变Y,看其能否模拟X干预所产生的P;随后在X干预背景下反向调节Y,看能否削弱这一表型。于是可以形成两类互补证据:Y manipulation能够模拟X的部分表型;Y rescue能够削弱X产生的表型。这种组合比“X改变以后Y和P同时变化”提供了更强的功能支持,也比只做一次Y单独干预更能够回答Y在X作用背景中的位置。
经过这一阶段,前面几部分逐步积累的证据可以连接起来:
RNA-seq发现Y发生变化 → 独立样本证明Y能够稳定重复 → 蛋白或功能状态验证Y确实发生生物学变化 → 单独操纵Y能够改变P → Y产生的表型方向符合X—Y—P模型 → X与Y联合干预时,反向调节Y能够部分或明显恢复P。
到这里,已经可以较为谨慎地提出:Y是X影响核心表型P过程中具有功能贡献的中介节点之一。

