不同光周期和温度处理下桂花内参基因的筛选

王千千，蒋琦妮，付建新，董 彬，赵宏波

（浙江农林大学 风景园林与建筑学院，浙江 杭州，311300）

桂花Osmanthus fragrans是中国十大传统名花之一，深受人们喜爱。不同学者已经从桂花的花芽分化［1-3］、 开花机制［4-5］、 花色呈现［6-7］、 花香释放［8-9］等不同层面展开了研究。 随着分子水平研究的不断深入，利用实时荧光定量聚合酶链式反应（qRT-PCR）分析基因的表达已成为研究这些过程机理的重要手段，而筛选得到适合不同研究的内参基因格外重要。研究认为［10-11］，内参基因并不存在通用性，不同处理、不同组织中内参基因的表达水平都会发生改变［12］；选用一种内参基因无法准确反映不同植物组织、同一组织不同处理下的基因表达水平［13-14］。因此，根据具体的实验材料和处理条件筛选合适的内参基因十分必要［15］。光周期和温度处理能对桂花的基因表达水平、花芽分化进程等产生影响。研究表明［4，16］：长日照和相对低温都能够加快桂花花芽分化进程，而短日照条件下桂花的花芽分化相对较慢，低温或高温都会抑制桂花的花芽分化。本研究在不同光周期和温度条件下对桂花常用的7个候选内参基因进行比较分析，筛选最佳内参基因，为后期研究光周期和温度影响桂花开花的机理奠定基础。

1 材料与方法

1.1 材料

实验材料为多年生 ‘佛顶珠’O.fragrans‘Foding Zhu’盆栽植株，种植于浙江农林大学桂花资源圃。设置4个处理：长日照高温（光照/黑暗=14 h/10 h，26℃），长日照低温（光照/黑暗=14 h/10 h，19℃），短日照高温（光照/黑暗=10 h/14 h，26℃）和短日照低温（光照/黑暗=10 h/14 h，19℃），光照强度为100%，相对湿度为60%～80%。实验材料在上述环境中培养40 d。处理结束后从8:00起，每隔4 h取当年生枝条的第2片或第3片嫩叶，液氮速冻，放于-80℃保存。

1.2 方法

1.2.1 总RNA的提取和cDNA的合成 每样品取50～100 mg，充分研磨后按照RNAprep pure Plant Kit（天根，北京）说明书步骤提取各样品RNA。通过紫外分光光度计和质量浓度为1%的琼脂糖凝胶电泳检测所提RNA浓度和品质。以所提样品RNA为模板，按照Reverse Transcriptase M-MLV（Takara，大连）说明书合成cDNA第1条链，合成后于-20℃储存备用。

1.2.2 候选内参基因的选择和特异性引物设计 从桂花转录组数据库中选择肌动蛋白基因（OfACT）、延伸因子1α蛋白基因（OfEF1α）、NADP-异柠檬酸脱氢酶基因（OfIDH）、GTP结合蛋白RNA1基因（OfRNA1）、β 微管蛋白基因（OfTUB）、 泛素结合酶 E2 基因（OfUBC2）和 18S 核糖体 RNA 基因（Of18S）［12，17］等7个使用较多的内参基因作为候选内参基因。采用SYBR greenⅡ荧光染料嵌合法进行内参基因的筛选，设计候选内参基因的引物，通过Bioxman和Primer Premier 5软件对引物进行校验，并利用Primer-BLAST进行检测以确认引物序列的特异性。选用的引物序列由上海捷瑞生物工程有限公司合成（表1）。

1.2.3 荧光定量PCR分析 qRT-PCR实验方法参考SYBRPremix Ex TaqTMⅡ（Tli RNaseH Plus）（Takara，大连）试剂盒。采用20.0 μL的反应体系：其中双蒸水 6.4 μL，上下游引物各0.8 μL，cDNA模板2.0 μL，SYBRⅡ荧光染料10.0 μL。反应在ABI 7300实时PCR仪中进行，反应程序为：95℃下预变性30 s；95℃下反应 5 s，60℃反应 31 s，共40个循环；95℃下反应15 s，60℃下反应 1 min，95℃下反应30 s， 60 ℃下反应 15 s［17］。 3次生物学重复。

1.2.4 数据分析 用geNorm，NormFinder和BestKeeper软件分别对qRT-PCR中所得到样品的反应循环数（Ct值）进行分析，比较7个候选内参基因表达的稳定性，并选择最佳内参基因。geNorm软件以各候选内参基因在样品中的表达水平为依据，计算各候选基因的表达稳定性（M），M值越小越稳定；M＜1.5时认为是理想内参［18-19］。geNorm软件还可以根据候选内参基因标准化因子的配对差异性分析（Vn/Vn+1）得出最佳内参基因的个数；其中V表示配对差异值，n表示内参基因个数。软件默认的（Vn/Vn+1）值为0.15，当Vn/Vn＋1＜0.15时，表明引入的n个内参基因已经稳定，不需要引入（n＋1）内参基因，当Vn/Vn＋1＞0.15时，则需要引入（n＋1）个内参基因［17，20］。因geNorm软件最终挑选出2个或2个以上的内参基因组合来校正数据，所以n≥2。为了了解引用1个内参基因的稳定性，需要通过NormFinder和BestKeeper软件做进一步分析。NormFinder软件对内参基因的筛选是基于方差分析的结果，利用△Ct法计算候选内参基因表达的稳定值（S），从而评估候选内参基因的表达稳定性；S越小，基因越稳定［21］。BestKeeper软件通过计算标准偏差（SD）对候选内参基因的表达稳定性进行评价和排序，SD越小，则该基因的稳定性越好；若SD＞1，则认为该基因不稳定［20-21］。最后通过昼夜节律基因GI的相对表达水平，对选择的最佳内参基因进行验证。

表1 qRT-PCR相关内参基因引物信息Table 1 List of primers for qRT-PCR

2 结果和分析

2.1 RNA品质及引物扩增效率检测

提取不同处理桂花样品的总RNA，经紫外分光光度计检测，所有样品总RNA的吸光度值［D（260）/D（280）和D（260）/D（230）］均为1.8～2.1。将所有样品的cDNA等量混合后作为模板，进行普通PCR和定量PCR检测。为检测引物的扩增效率，设置6个cDNA浓度梯度（原液，5-1，5-2，5-3，5-4，5-5原液），根据qRT-PCR结果，利用公式E=（10-1/s-1）×100%得出基因的扩增效率（表2），其中：E（%）表示扩增效率，s表示曲线斜率。经计算，扩增效率为95%～110%，相关系数为0.993～0.998，符合qRT-PCR的基本要求。

2.2 候选内参基因的转录水平分析

根据扩增效率，将各样品的原始cDNA稀释50倍后作为模板进行qRT-PCR。通过比较循环数（Ct）来评估各候选内参基因在所有样品中的转录水平（表2），候选内参基因平均Ct为8.99～28.44，其中Of18S的Ct最小，说明平均表达水平最高；OfTUB最大，说明平均表达水平最低；其他候选内参基因平均表达水平差异则相对较小。

2.3 内参基因表达稳定性分析

geNorm软件分析发现（表3），所有候选内参基因的稳定值均小于1.5，符合内参基因筛选的标准；其中OfRAN1和OfIDH的稳定值相等且最小（M=0.347），在不同样品中表达最为稳定，其次是OfACT基因（M=0.393），Of18S稳定值最大（M=0.601），说明Of18S表达稳定性最差。对候选内参基因的配对差异值测算可知（图1）：V2/3=0.126＜0.15，表明引入2个内参基因已经稳定，无需引入更多的内参基因。geNorm软件分析结果表明，最佳候选内参基因为OfRAN1和OfIDH。

表2 7个候选内参基因扩增参数Table 2 Amplification parameters of seven candidate reference genes

NormFinder软件分析发现，OfRAN1基因的表达稳定性最高，其次是OfIDH，Of18S基因表达稳定最差，与geNorm软件分析的结果一致（表3）。BestKeeper软件的分析结果显示（表3），在7个候选内参基因中，OfIDH基因的表达稳定性最高，其次是OfRAN1，Of18S基因的表达稳定性最差。分析结果尽管与geNorm和NormFinder软件的结果有些差异，但最佳候选内参基因结果一致。

表3 不同软件分析7个候选内参基因表达稳定性Table 3 Gene expression stability of seven reference genes by different software

综合3个软件的分析结果可以得出，不同处理条件下，OfRAN1和OfIDH的表达稳定性均为最高，可作为最佳内参基因，而Of18S是表达稳定性最差的基因。

图1 geNorm软件分析候选内参基因的配对差异值Figure 1 Pairwise variationof candidate reference genes as calculated by geNorm

2.4 桂花GI基因对最佳内参基因稳定性的验证

GIGANTEA基因（GI）是光周期途径的重要调控基因，其表达呈现明显的昼夜节律。选用桂花昼夜节律相关基因GI对候选内参基因进行验证，可提高内参基因筛选的准确性，保证最佳内参基因对光周期途径基因表达水平的校准作用。选择OfRAN1/OfIDH，OfRAN1，OfIDH以及其他候选内参基因对GI基因的相对表达水平进行校准（图2），结果显示：相比其他内参基因，不同处理下OfRAN1/OfIDH校准的GI基因的表达模式表现为光照条件下积累，黑暗条件下降低的周期性表达，符合GI基因的昼夜表达规律。说明3个内参基因筛选软件对候选内参基因筛选的结果是正确的，即OfRAN1和OfIDH基因可作为桂花不同光周期和温度处理的最佳内参基因。

3 讨论

荧光定量PCR因其众多优点被广泛运用于基因表达的研究中，成为量化基因转录表达，揭示基因表达规律的有效方法［22-23］；但引用适当的内参基因对qRT-PCR的结果进行校准是非常必要的［20］。理想内参基因在不同组织、不同发育阶段、不同生长条件下均能稳定表达［17，24］，但事实上并没有哪一个基因能够满足所有条件［25］。因此，根据不同实验材料和处理条件，筛选出合适的、特定的内参基因是确保实验成功的必要条件。

图2 不同内参基因或组合校正桂花不同条件下GI基因的相对表达Figure 2 Relative expression levels of GI in different conditions using different reference genes or reference gene combination for normalization

研究发现［26］：选用2个或更多内参基因组合可以提高内参基因校准的精确度，弱化单一内参基因带来的不准确性。RAN（ras-related nuclear protein）是小G蛋白的一类，在细胞生理进化中起到重要的作用，它的很多同源基因是不同光周期、温度条件下常用的内参基因［12，17］。IDH的同源基因被证实为不同温度条件下桉树Eucalyptus robusta的最佳内参基因［27］；ACT和18S同源基因分别可作为不同光周期［28］、温度条件［29］下大豆Glycine max的内参基因；EF1a和TUB基因可作为不同光周期、温度条件下玉米Zea mays的内参基因［30］；UBC2的同源基因可作为不同温度条件下欧芹Petroselinum crispum的内参基因［31］。本研究选用了桂花中常用的7个候选内参基因，对其在不同光周期和温度处理下嫩叶中的表达稳定性进行研究；利用geNorm，NormFinder和BestKeeper等3个常用于内参基因筛选的软件进行评估，发现这7个候选内参基因在所有样品中表达稳定性的排序基本相同，均显示OfRAN1和OfIDH基因是最稳定的2个内参基因，将两者进行组合可以准确的校准不同处理下桂花叶片中目的基因的表达水平；而Of18S表达最不稳定，在所有实验材料中都存在较高的表达丰度，因而不适合作为本研究的内参基因。

光周期和温度途径在调控植物成花过程中发挥重要作用。筛选得到的桂花不同光周期和温度条件下的内参基因，能够为研究桂花在不同温度和光周期处理条件下的成花机制提供参考依据，为了解相关成花基因表达模式提供保障。