在分析功能基因的表达水平时，为了消除不同样品在RNA的提取质量和得率及在RNA反转录成cDNA的效率上存在的差异，通常选用适当的内参基因进行校正，以确定目的基因的相对表达量。理想的内参基因要求表达量适中，表达水平稳定，不受内源性(如生长周期或不同器官)或外界条件(如实验措施或环境因素)的影响。

植物肌动蛋白(Actin)是一类古老的、由单一多肽构成的球状蛋白质^[1]，构成细胞的微丝系统，参与细胞分裂和形态控制、信号传导、细胞器运动等细胞生理活动^[2-3]；3-磷酸甘油醛脱氢酶(Glyceraldehyde-3-phosphate dehyrogenase, GAPDH)是糖酵解过程中的关键酶，氧化3-磷酸甘油醛形成1, 3-二磷酸甘油酸^[4]，与生命活动紧密相关；微管蛋白(tubulin)不仅是细胞骨架的组成部分，而且与胞内物质运输、信号识别、细胞运动等关系密切，主要由α、β两种微管蛋白以二聚体的形式聚合而成^[5]。Actin、GAPDH和tubulin表达相对稳定，常作为内参基因用于基因表达分析^[6-9]。但是，已有研究发现，持家基因的表达水平也存在不稳定性^[10-11]，应根据不同试验目的筛选合适的内参基因。

金线莲Anoectochilusroxburghii (Wall.) Lindl是兰科Orchidaceae开唇兰属Anoectochilus植物，又称金线兰、金草，具有保肝、降血糖、抗肿瘤、提高免疫力等作用^[12-13]，主要分布在我国的亚热带地区，福建和台湾是主产地^[14]，因野生资源匮乏，主要靠种植组培苗满足市场的需求。金线莲富含黄酮、生物碱、皂苷、甾体化合物、三萜类等活性成分^[15-17]，而组培苗的种植时间会影响金线莲活性成分含量^[18]。姜福星等^[19]对台湾金线莲进行转录组测序，分析发现与黄酮类、生物碱和萜类物质代谢相关的Unigenes分别有130、114和441条，而有关金线莲活性成分合成功能基因的研究还处于起步阶段。研究金线莲药用成分合成调控关键酶基因的功能，分析关键酶基因表达模式，必须选择理想的内参基因，而目前有关金线莲内参基因的研究还鲜有报道。本研究从金线莲转录组测序数据中选择GAPDH、Actin和α-tubulin等3个Unigenes，根据序列信息设计引物，应用半定量RT-PCR技术分析不同生长期和不同温度处理条件下3个基因在金线莲茎、叶的表达水平，以期为金线莲药用活性成分合成关键酶基因的研究选择合适的内参基因。

1.   材料与方法

1.1   试验材料

金线莲(采自福建省龙岩市梅花山)生根瓶苗为本课题组培养。pMD19-T连接试剂盒、Ex Taq DNA Polymerase、PrimeScript^TM Reverse Transcriptase、RNAiso Plus购自宝生物公司。UNIQ-10柱式DNA胶回收试剂盒购自上海生工生物工程有限公司，引物合成、克隆测序委托上海生工。

1.2   金线莲温度处理

将光照培养箱的温度设置为5℃(T5)、15℃(T15)、25℃(T25)、35℃(T35)，光照强度4 000 lx，待温度稳定后放入可以移栽的组培生根瓶苗，培养10 d，随机剪取15株金线莲的茎、叶，用铝箔分开包装，液氮速冻后置超低温冰箱保存。

1.3   金线莲不同生长阶段处理

2017年11月3日种植，将金线莲组培生根苗从瓶子取出，洗净根部的培养基后种植于育苗盘，每盘30株，光照强度在3 000~5 000 lx，栽培基质为泥炭土，栽培过程根据基质的干湿情况及时补水。2017年12月至2018年5月，每月取样1次(分别用P1~P6表示)，随机剪取15株金线莲的茎、叶，用铝箔分开包装，液氮速冻后置超低温冰箱保存。

1.4   总RNA的提取和cDNA合成

利用RNAiso Plus提取金线莲茎、叶总RNA，经电泳检测后，用超微量紫外可见光分光光度计(ND-1000)测定RNA纯度并定量为500 ng·μL^-1。

cDNA合成：在200 μL PCR管中顺序加入DEPC水3 μL，总RNA 2 μL，随机引物(50 μmol·L^-1)1 μL，dNTP Mixture(各2.5 mmol·L^-1)4 μL，65℃保温5 min后，迅速在冰上急冷5 min，离心数秒后加入预先混合好的5×PrimeScript^TM Buffer 4 μL，RNase inhibitor (40 U·μL^-1)0.5 μL，PrimeScript^TM Reverse 0.5 μL，DEPC水5 μL，混匀后放入PCR仪，PCR仪反应程序为：30℃ 10 min，42℃ 60 min，70℃ 15 min，4℃结束反应。反应结束后利用超微量紫外可见光分光光度计(ND-1000)测定吸光值并定量为200 ng·μL^-1。

1.5   引物设计及扩增特异性分析

根据金线莲转录组测序数据中的Actin、GAPDH和α-tubulin等3个Unigenes序列信息，利用primer 5.0软件设计扩增引物(表 1)。以茎、叶混合的cDNA为模板PCR扩增，验证引物特异性。扩增体系20 μL：模板1 μL，引物(10 μmol·L^-1)各1 μL，10×Buffer 2 μL，Ex Taq DNA聚合酶(5 U·μL^-1) 0.1 μL，dNTP Mixture(各2.5 mmol·L^-1) 1.6 μL，ddH₂O 13.3 μL。反应程序为：94℃预变性5 min，94℃变性30 s，退火(Actin 50℃、GAPDH 55℃、α-tubulin 55℃)30 s，72℃延伸20 s，30个循环，最后72℃延伸10 min。胶回收目的片段，与克隆载体pMD19-T连接后转化大肠杆菌DH5α感受态细胞，经验证后送上海生工测序。

表  1  PCR引物及其序列

Table  1.  PCR primers and sequences

靶基因	引物名称	引物序列(5′-3′)	产物长度/bp
Actin	ArACT-F	AGATGAGGCACAGTCCAAGA	228
	ArACT-R	GCTGGAACATTGAAGGTCTC
GAPDH	ArGAPDH-F	CCCTTTGGACCTGTGCTACC	162
	ArGAPDH-R	GATTTGAACTGTGCCCGTCTC
α-tubulin	ArATUB-F	GCTGTTGGAGGTGGGACTG	162
	ArATUB-R	AGAGTGTGTGGACAGGACGC

下载: 导出CSV 

| 显示表格

1.6   循环参数优化与表达分析

以不同处理样品总RNA反转录获得的cDNA为模板进行PCR，优化循环参数。扩增体系40 μL：模板2 μL，引物(10 μmol·L^-1)各2 μL，10×Buffer 4 μL，Ex Taq DNA聚合酶(5 U·μL^-1) 0.2 μL，dNTP Mixture(各2.5 mmol·L^-1) 3.2 μL，ddH₂O 26.6 μL。反应程序为：94℃预变性5 min，94℃变性30 s，退火30 s，72℃延伸20 s，20个循环开始取样5 μL，3个循环取样1次，取样5次，确定循环数。以优化的循环再进行PCR，重复3次，分析3个基因的表达。

2.   结果与分析

2.1   引物扩增特异性分析

为了验证设计引物扩增的特异性，以合成的引物ArACT-F和ArACT-R、ArGAPDH-F和ArGAPDH-R、ArATUB-F和ArATUB-R进行PCR扩增，从图 1可以看出，3对引物均能扩增特异性条带，无非特异性扩增，目的条带的测序结果经DNAMAN V6.0比对后确认均是目的片段，可作进一步分析。

图  1  PCR电泳结果

注：M为Marker；1为GAPDH；2为α-tubulin；3为Actin。

Figure  1.  PCR products by electrophoresis

下载: 全尺寸图片 幻灯片

2.2   Actin在不同处理温度和不同生长期的表达分析

以不同处理温度和不同生长期的金线莲茎、叶提取的总RNA反转录的cDNA为模板，优化PCR扩增循环参数，确定Actin扩增循环数为26。以优化的循环数进行PCR扩增，分析Actin在金线莲茎、叶中的表达情况。从图 2可以看出，不同温度处理和不同生长期，Actin在茎、叶中的表达量基本一致，表达稳定，可作为温度处理和生长期金线莲功能基因表达分析的内参基因。

图  2  不同温度处理和不同生长期Actin的表达分析

注：M为Marker，T₅、T₁₅、T₂₅、T₃₅分别表示不同处理温度；P₁~P₆表示不同生长时期，A为茎，B为叶。图 3、4同。

Figure  2.  Expression of Actin gene in response to temperature changes in different growth stages

下载: 全尺寸图片 幻灯片

图  3  不同温度处理和不同生长期GAPDH的表达分析

Figure  3.  Expression of GAPDH gene in response to temperature changes in different growth stages

下载: 全尺寸图片 幻灯片

图  4  不同温度处理和不同生长期α-tubulin的表达分析

Figure  4.  Expression of α-tubulin gene in response to temperature changes in different growth stages

下载: 全尺寸图片 幻灯片

2.3   GAPDH在不同处理温度和不同生长期的表达分析

以各种处理金线莲茎、叶提取的总RNA反转录的cDNA为模板，经过PCR扩增循环参数优化，确定GAPDH扩增循环数为27。以优化的循环数进行PCR扩增，分析GAPDH在金线莲茎、叶中的表达情况。从图 3可以看出，不同温度处理GAPDH表达存在明显差异，随着温度的降低，GAPDH在茎、叶中的表达量越来越高，5℃的表达量最高，说明低温能上调GAPDH的表达。不同生长期GAPDH的表达量也存在差异，随着种植时间的延长，GAPDH在茎、叶中的表达量越来越高，种植4~5个月表达量最高，种植6个月的表达量下降。

2.4   α-tubulin在不同处理温度和不同生长期的表达分析

以金线莲茎、叶各种样品提取的总RNA反转录的cDNA为模板，经优化α-tubulin扩增循环数确定为28。以优化的循环数进行PCR扩增，分析α-tubulin基因在金线莲茎、叶中的表达情况。从图 4可以看出，金线莲茎15℃处理的表达量较高，其他3个温度基本一致，叶中15℃和25℃表达量比5℃和35℃稍高。随着种植时间的延长，α-tubulin在茎、叶中的表达也越来越高，种植5个月的表达量均最高，种植6个月的表达量下降。

3.   讨论

持家基因不仅参与细胞的代谢过程，还具有其他的生理功能，不同样品、不同试验因素的内参基因不同，要根据试验目的不同，分析选择适宜的持家基因做内参，提高研究的可能性和准确性。

梁云等^[20]克隆发现百合lilyActin在花蕾、叶片和鳞片中表达稳定。黄连香等^[21]克隆发现三叶青肌动蛋白基因ThAct1和ThAct2在茎、普通根和块根中都稳定表达，但ThAct1比ThAct2表达量高；叶与茎、根的表达量存在差异，且叶中ThAct2的表达量稍高，而ThAct1表达量较低。亓希武等^[22]克隆了金银花LjActin在幼蕾期、绿蕾期、白蕾期、银花期和金花期等5个发育阶段的花中稳定表达。敖特根白音等^[23]克隆了蒙古冰草MwACT2，在高盐、干旱、低温等逆境胁迫下表达稳定。本研究中金线莲肌动蛋白基因在不同温度处理和不同生长期均能稳定表达，是研究金线莲在温度胁迫或不同生长期基因表达分析的理想内参基因。

已有研究发现植物GAPDH表达并不稳定，可被高温、低温和盐胁迫等诱导。如卜娇迪等^[24]克隆了金丝小枣的ZjGAPDH，定量分析发现该基因在果实不同发育阶段表达量存在差异，全红期表达量最高，半红期最低。张霞等^[25]研究发现盐穗木GAPDH受盐胁迫和ABA上调表达，盐胁迫6 h、ABA处理12 h时达到最高值。本研究也发现低温上调金线莲GAPDH的表达，不同生长期GAPDH在茎、叶中的表达量也不一致，不适宜在分析温度处理和生育期金线莲功能基因表达时作为内参基因。

Jian等^[26]分析了大豆不同发育时期10个内参基因，发现α-tubulin表达稳定性最好。然而，本研究中发现金线莲α-tubulin在不同的温度和生长期的表达并不稳定，不适宜在分析温度处理和生育期金线莲功能基因表达时作为内参基因。

本研究选择金线莲的GAPDH、Actin和α-tubulin等3个持家基因，应用半定量RT-PCR技术分析了不同生长期和不同温度处理条件下3个基因在茎、叶的表达水平，没有进行准确的表达定量分析，其他环境条件下金线莲持家基因的表达稳定性也有待进一步研究。