农药施用减量增效要求在减少农药施用量的同时提高农药的使用效率，降低农药残留，保证食品安全，这也是当前农药研究者共同努力的方向。徐汉虹等^[1]提出的“导向农药”理念完美契合农药减量增效的目的，对保障食品安全具有非常重要的意义。导向农药能够在植物内源转运蛋白的作用下定向积累至病虫为害部位，它的作用方式可以分为2种：一是当植物受到危害时，导向农药分子能够被植物应激反应产生的酶消解并释放出农药母体分子作用于靶标；二是具有农药活性的导向分子直接作用于靶标从而达到防治病虫害的目的。植物内源糖从源叶合成，在糖转运蛋白介导下源源不断地通过韧皮部转运至植物库组织^[2]；基于导向农药的理念，Yang等^[3]合成了葡萄糖基氟虫腈耦合物，并发现其有可能是通过蓖麻糖转运蛋白介导进入蓖麻韧皮部^[4]。充分利用植物对内外源化合物的吸收与转运机制对于提高农作物产量，探索病虫害防治新策略具有非常重要的意义。

利用植物对内外源糖物质的吸收转运规律，对关键基因进行修饰从而提高作物产量；合成能够被植物转运载体介导的新型农药分子，实现农药分子在植物体内的主动运输和定向积累，对于病虫害防治能起到事半功倍的效果。研究发现草甘膦、百草枯、葡萄糖基氟虫腈、部分农药分子与氨基酸的耦合物L-lysine-2, 4-D等虽然能被磷酸转运蛋白、聚胺转运蛋白、糖转运蛋白、氨基酸转运蛋白所介导^[4-7]，但却无法明确具体的转运子，无法明确化合物与转运载体之间的构效关系，大大降低了这类农药的开发效率。

蓖麻Ricinus communis 是用来研究内外源化合物韧皮部输导性的模式植物，随着其基因组数据库的释放，蓖麻中的各类转运蛋白家族将会被鉴定，利用RT-qPCR技术分析并克隆能被载体介导的化合物处理前后差异表达的基因，是鉴定参与化合物韧皮部装载的相关载体蛋白的重要手段。RT-qPCR自1995年由美国PE公司成功开发以来，便以其操作简便、检测灵敏度高、重复性好、检测速度快等优点，被广泛应用于生命科学、分子诊断学、农学与医学、食品微生物学等领域^[8]。在基因表达分析过程中，RNA的质量和产量、不同组织或细胞间反转录效率的差异等因素均有可能影响到目的基因表达分析结果的准确性^[9]。为了获取准确的结果，除了提高RNA质量，选择合适的内参基因来校正和标准化mRNA是最首要的方式^[10]。理想的内参基因的表达不受试验因素的影响而恒定表达在各种类型的组织或细胞中^[11]。常用的内参基因是构成细胞器骨架的基本组分或是参与生物体的基本生化代谢过程的基因^[12-14]。同一个内参基因在不同的试验条件下表达的稳定性是不一样的，如玉米根萤叶甲Diabrotica virgifera virgifera Leconte在不同的试验条件下β-Actin在其整个生命周期内表达稳定，RPS9则在dsRNA处理后的玉米根萤叶甲中表达稳定性最好，而BT处理后，EF-1在不同的组织中表达最为稳定^[15]；在荔枝果实不同发育阶段用外源生长调节剂处理后β-Actin比GAPDH、18S rRNA、UBQ、eEF和25S rRNA的稳定性更好^[16-17]，由此可见筛选出化合物处理前后都能稳定表达的内参基因是RT-qPCR技术分析并克隆能被载体介导的化合物处理后差异表达的基因的基础。

本研究以蓖麻幼苗子叶为材料，选取常用的内参基因Actin、ARC、SamDC、ef1a、TUA6作为候选基因，用RT-qPCR方法检测上述常用内参基因在不同浓度及不同时间的GTF及其溶剂DMSO处理下的表达量，分别采用常用的内参基因稳定性评价软件geNorm^[10]、NormFinder^[18]、BestKeeper^[19]、Delta CT^[20]进行数据分析，最后用RefFinder^[21]在线工具对上述4种软件的计算结果进行综合评价，为研究GTF相关的转运载体基因的表达分析研究筛选出最为稳定的内参基因。

1 材料与方法 1.1 材料

蓖麻种子由淄博市农业科学院提供。成熟蓖麻种子在湿脱脂棉上于27 ℃下催芽48 h后，选出发芽的种子，播入蛭石中，置于人工气候箱内进行培养。培养条件为：温度为27 ℃，相对湿度为75%~85%，光照度260~350 lx，每天光照14 h, 连续培养6 d的蓖麻幼苗作为供试材料。

1.2 方法 1.2.1 蓖麻幼苗的药物处理

将蓖麻幼苗小心剥去胚乳，用手轻轻擦洗子叶表面，切勿折叠或弄破子叶。用子叶培养缓冲液(20 mmol·L^-1 MES, 0.25 mmol·L^-1 MgCl₂, 0.5 mmol·L^-1 CaCl₂, 2 mmol·L^-1 EDTA, pH 5.0)配制相应浓度的GTF或者DMSO后，将蓖麻子叶浸泡于相应含药缓冲液中，同时将蓖麻根部置于根部培养缓冲液(0.5 mmol·L^-1 CaCl₂, 1 mmol·L^-1 MES)中。浓度梯度试验设定7个GTF和DMSO浓度梯度，分别为10、25、50、100、200、500、1 000 μmol·L^-1，于人工气候箱内培养6 h后收集样品；时间梯度试验设置为以100 μmol·L^-1的GTF或DMSO处理蓖麻幼苗，分别于0、3、6、9、12 h后收集子叶。每个处理设3个重复，以不含GTF及DMSO的缓冲液为对照组，完成相应的处理后在蓖麻幼苗胚柄处剪断，将子叶装入锡纸袋，用液氮速冻后置于-80℃冰箱保存备用。

1.2.2 蓖麻子叶总RNA的提取及第一链cDNA的合成

使用Omega植物RNA提取试剂盒提取经过处理的蓖麻子叶总RNA，并用核酸蛋白定量仪(ND1000，Eppendorf)检测RNA样品的D_{260 nm}/D_{280 nm}和D_{260 nm}/D_{230 nm}值，确定RNA的浓度和纯度，同时用含甲醛(w=12.5%)的琼脂糖凝胶电泳检测RNA的完整性。随后RNA用RQ1 RNase-free DNase处理，并用反转录试剂盒(Bio-Rad)进行第一链cDNA的反转录。转录完成后以合成的cDNA为模板克隆内参基因片段以进行RT-qPCR试验。

1.2.3 候选内参基因的处理

根据水稻内参基因的氨基酸序列^[20], 在NCBI上获得对应的蓖麻内参基因登录号和基因序列，并通过Primer Premier 5.0软件设计特异性引物，由北京六合华大基因科技有限公司合成。选取的内参基因引物序列见表 1。

以含候选基因片段的质粒为模板构建其对应的标准曲线。利用cDNA为模板获得候选内参基因片段。通过20 g·L^-1凝胶电泳检测并用DNA凝胶纯化试剂盒回收目的片段，将回收产物连接pMD20-T载体上，蓝白斑筛选阳性克隆并测序。经验证正确的阳性克隆扩大培养后，用质粒小提试剂盒(Omega)提取质粒。取500 ng质粒用Dilution Buffer(Takara)稀释, 每个梯度稀释10倍，依次稀释5~8个梯度，则浓度分别为50、5、1/2、1/20、1/200和1/2 000、1/20 000、1/200 000 ng·μL^-1。以稀释后的质粒为模板进行RT-qPCR反应，PCR仪自带的计算软件CFX Manager TM Software自动计算获得熔解曲线，并根据所得C_t值绘制标准曲线，得到斜率(k)和线性相关系数(R)。根据公式E=5^-1/k-1，计算得出扩增效率(E)。候选内参基因引物特异性则通过RT-qPCR的熔解曲线的峰形图来评价。

1.2.4 RT-qPCR分析

RT-qPCR反应用Bio-Rad CFX-96 Real-Time PCR System完成。反应体系按照IQ^TM SYBR^® Green Supermix说明书进行配制，反应总体积为20 μL；反应条件如下：95℃ 30 s，95 ℃ 5 s，62 ℃ 40 s，40个循环后进行65~95 ℃熔解曲线分析。每个样品设定3个重复。

1.3 数据分析

根据Bio-Rad CFX-96 Real-Time PCR System获得的各样品的C_t值，采用geNorm，NormFinder，BestKeeper, Delta CT 4种算法对各候选内参基因的稳定性进行比较，同时参考RefFinder的计算结果，筛选出在GTF及DMSO处理下稳定性最好的内参基因。

2 结果与分析 2.1 候选内参基因引物特异性及扩增效率

以反转获得的cDNA为模板扩增候选内参基因片段，电泳检测结果显示所得目的条带大小与理论值一致, 且没有发现其他条带，扩增特异性好(图 1)。

RT-PCR的结果表明，各候选内参基因的相关系数(R)≥0.996, 引物的扩增效率为95.0%~107.7%，均符合RT-qPCR对扩增效率的要求(表 2)；进一步分析熔解曲线，发现各候选内参基因都只产生单一的熔解峰。上述结果说明所设计的引物特异性良好，RT-qPCR反应的专一性高，准确性好。

2.2 候选内参基因的表达谱分析

在RT-qPCR试验中，C_t值与基因的表达量成反比。本研究结果表明同一个扩增反应的3个重复的标准差很小，C_t值非常接近，说明扩增反应重复性好，能准确反映基因的转录水平；5个候选内参基因的C_t值均位于17~24，同时都处于相应候选内参基因标准曲线的C_t跨距之内；同一个基因，浓度处理对其表达量影响较大；其中对ef1a 的表达量最大(17~21)，SamDC及ARC的表达量相对较小(18~24)；在所有的样品中，ef1a的C_t值变化范围最小，ΔC_t为2.77，其次是Actin，SamDC的C_t值变化最大，ΔC_t为4.23。同一个基因在不同的样品中表达量的差异说明有必要对这些常用内参基因进行筛选，以选择合适的内参基因或组合来进行接下来的定量试验。

2.3 不同处理候选内参基因稳定性分析 2.3.1 4种软件分析内参基因稳定性

将RT-qPCR试验所获得的C_t值处理后导入geNorm(V3.5)计算出各候选内参基因的稳定值(M值)，M值越大，表示稳定性越差。计算结果(图 2a)表明，5个候选内参基因均相对较稳定(M < 1.5)，Actin和ef1a的M值最低，稳定性最好；ARC的M值最高，稳定性最差；这5个候选内参基因的稳定性从高到低依次为：Actin=ef1a > SamDC > ARC > TUA6，即在geNorm的算法中，Actin与ef1a的稳定性最好。

按NormFinder软件要求，将处理好的数据导入软件进行计算，得出各候选内参基因的稳定值M，M值越大，说明基因的稳定性越差；反之则说明基因的表达稳定性越好。经计算获得5个候选内参基因的M值，按稳定性由高到低依次为：SamDC > ARC  > Actin > ef1a > TUA6(图 2b)，即在NormFinder算法中，SamDC的稳定性最好。

将RT-qPCR反应获得的C_t值导入BestKeeper软件，根据各候选基因标准偏差值(SD)来确定候选基因的相对稳定性，SD值越小，稳定性越好。在本次试验中，所有基因的SD值均小于1，其SD大小顺序依次为：Actin > ef1a > SamDC > ARC > TUA6(图 2c)，即在BestKeeper软件的算法中Actin的稳定性最好。

利用Delta CT软件计算5个候选内参基因的C_t值平均标准偏差(SD), SD越小，表明基因稳定性越高。将计算结果按稳定性由高到低进行排序，依次为：SamDC > Actin > ARC > ef1a > TUA6(图 2d)，即在Delta CT软件的算法中SamDC的稳定性最好。

2.3.2 RefFinder综合分析

由以上分析可知，geNorm，NormFinder，Delta CT和BestKeeper 4种软件由于不同的算法分析出来的结果也存在差异。此时可以利用RefFinder在线分析工具来进行综合分析从而选出表达相对最为稳定的内参基因。RefFinder可以根据各种算法所得的候选内参基因稳定表型排名的几何平均数作为综合排名，从而筛选出最稳定的内参基因。基因综合稳定参数与基因表达稳定性成反比。根据计算结果可知GTF及DMSO处理下，候选内参基因表达稳定性由高到低依次为：Actin > SamDC > ef1a > ARC > TUA6(图 4)。

2.3.3 DMSO处理下内参基因的稳定性分析

DMSO在蓖麻体系中的韧皮部输导性研究中常被用作溶剂，有必要对其处理的样品单独进行内参基因稳定性分析，以作参考。RefFinder分析结果(图 4)表明，在10~1 000 μmol·L^-1的浓度梯度及0~12 h的时间梯度的DMSO处理下ef1a是最稳定的内参基因。

3 结论

本研究比较了Actin、ARC、SamDC、ef1a、TUA6基因在不同浓度及不同时间的GTF及其溶剂DMSO处理前后的表达量，并分别采用geNorm、NormFinder、BestKeeper、Delta CT 4个软件及RefFinder在线分析工具对各候选内参基因的稳定性进行计算分析，以期获得相对稳定的用于校正和标准化RT-qPCR试验的一个或是一组内参基因。geNorm软件的计算结果中，Actin与ef1a为稳定性最好的组合；在NormFinder与Delta CT的算法中，SamDC是稳定性最好、最适合的内参基因；而通过BestKeeper获得的结果表明，Actin是最稳定的内参基因。由于各软件算法的不同导致计算结果不尽相同，难以取舍。而RefFinder在线工具可以根据以上4种算法所得的候选内参基因稳定性的几何平均数综合排名，筛选出表达最为稳定的基因，结果表明GTF、DMSO共同处理综合分析时Actin的稳定性最好。单独对DMSO处理后的内参基因稳定性分析时发现ef1a是最稳定的内参基因。本文分析了5个常用内参基因在GTF和DMSO共同处理及DMSO单独处理下内参基因的稳定性，可为相关研究提供参考。