竹亚科牡竹属的巨龙竹(Dendrocalamus sinicus Chia et J. L. Sun)为我国云南特有大型丛生竹种，其胸径可达30 cm，秆高可达30 m以上，是目前已知世界上最高大粗壮的竹种[1]。巨龙竹秆形在自然分布区中存在两种遗传稳定的变型，即“通直型”和“弯曲型”[2]。弯曲型竹秆基部常畸变，呈短缩、鼓胀、歪扭等形状，不利于竹板材加工；通直型巨龙竹秆形通直，平均单株质量100150 kg，单位面积年均产材量比我国主要经济竹种毛竹(Phyllostachys edulis(Carr.) H. de Lehaie.)高58倍[3]。通直型巨龙竹是培育特大型工业用材的优异竹种之一，也是我国南部热区竹产业发展和生态环境建设中最具发展潜力的优良经济竹种之一[3-4]。秆形作为巨龙竹资源培育与利用的关键性状之一，其发育机制尚不清楚，这给巨龙竹优良种源的选择和培育带来了障碍。同时，“通直型”和“弯曲型”巨龙竹作为一对自然的变异类型，也是木本竹种秆形发育研究中理想的对照系。目前，作者正在开展巨龙竹秆形发育的转录组分析，试图解析不同秆形发育过程中关键基因的时空表达模式。

实时荧光定量PCR(quantitative real-time PCR，qRT-PCR)技术因其特异性强、灵敏度高、快速准确、重复性好等[5-6]优点被广泛用于基因表达分析。在植物学研究中，常选用稳定表达的内参基因转录延伸因子-1α(EF-1α)、甘油醛-3-磷酸-脱氢酶(GAPDH)、肌动蛋白基因(Actin)、微管蛋白基因(Tubilin)为内参基因[7-8]。随着基因芯片和基因组学的发展，又挖掘了新的内参基因蛋白磷酸酶2A(PP2A)和液泡膜内嵌蛋白基因(TIP-41)等[8]。各内参基因在不同组织、不同发育时期、不同处理中的稳定性也存在差异[8]。近期有关竹亚科植物内参基因的筛选已有一些相关报道，主要集中在毛竹等温带散生竹中，筛选出PP2A、TIP41等作为毛竹不同组织基因表达定量分析的内参基因[9-11]，但这些内参基因在不同发育时期的表达稳定性是未知的。针对巨龙竹不同秆形发育特点，本研究用qRT-PCR评价EF-1α、GAPDH、Actin、Tubulin、TIP-41、PP2A 6个内参基因在巨龙竹茎秆发育过程中的表达稳定性，筛选在巨龙竹不同生长阶段、不同秆形变异种群中最佳内参基因，为下一步解析参与巨龙竹秆形发育关键基因的表达模式提供技术基础。

1 材料与方法

1.1 材料

1.1.1 供试材料 以云南省西盟县中课种群(通直型)和孟连县景信种群(弯曲型)2个种群内01个月的笋(或幼竹)作为研究对象。巨龙竹发笋期(5月下旬至8月下旬)的定点观察结果显示，秆形变化主要出现在竹秆基部13 m内的竹节，1个月内秆形基本定型。因此，从每个类型种群中选择3个竹丛，每个竹丛选取出笋后5天、15天、30天的竹笋，分别代表秆形发育的初期、中期和后期。样品取自竹笋基部竹节，用干净的刀片切成小块立即放入液氮中，液氮处理后于-80℃贮存备用。

 

表1 采样具体信息Table1 Sampling information

  

地点Site茎秆类型Culmtype经度Longitude纬度Latitude海拔Altitude/m采样时间Time西盟县中课镇窝笼村通直型99°32'06″22°43'48″12542016-07-15孟连县景信乡景信村弯曲型99°36'31″22°26'36″10072016-07-16

1.1.2 主要试剂 RNA提取试剂盒(TransZol Up Plus RNA Kit)，反转录试剂盒(TransScriptⅡ One-Step gDNA Removal and cDNA Synthesis SuperMix)，qRT-PCR荧光染料(TransStart Top Green qPCR SuperMix)均购自北京全式金生物技术有限公司。

1.2 方法

1.2.1 总RNA的提取 每个组织取3个不同的生物学重复。总RNA的提取参照TransZol Up Plus RNA Kit试剂盒说明书进行，提取的总RNA经2%琼脂糖凝胶电泳检测RNA完整性，Nanodrop 8000紫外分光检测仪检测OD值后直接用于反转录或保存于-80℃备用。

1.2.2 反转录合成cDNA 按照TransScript Ⅱ One-Step gDNA Removal and cDNA Synthesis SuperMix试剂盒说明书，各取样品的总RNA 7μL反转录成cDNA第一链。所得cDNA产物直接用于PCR和qRT-PCR或保存于-20℃备用。

1.2.3 内参基因的引物设计 使用Primer Premier 5.0软件设计候选内参基因的qRT-PCR引物，引物序列见表2。引物由昆明硕擎生物技术有限公司合成，采用PAGE纯化。

 

表2 qRT-PCR检测中所用的6个候选内参基因的引物序列Table2 Primer sequences of six candidatereference genes for qRT-PCR analysis

  

基因Genename引物序列(5’→3’)Primersequence(5’→3’)解链温度Tm/℃扩增产物长度Lengthofamplifiedproduct/bpPP2AF:CAGTCATCCCATGCTCTTGCT60.6158R:GAAGCTGGCTTACGTGAGAT-GAA59.8EF-1αF:GCAACAAAACCACGCTTGAGA59.6186R:GGGTCGTGTTGAGACTGGTG59.5GAPDHF:AGTCCACTGGTGTCTTCACTG59.6178R:CAATTTGTGGTGCAGCTAG-CATT60.0ActinF:TGCCTACATTGCCCTGGACTA59.6162R:TCCATTCCAATGAAGGATG-GCT59.8TubulinF:GTCGATGCGCGAGAATACCT59.4186R:GCAGTTGCAACCATCAAGACC59.6TIP-41F:CAGCGAAATGGAAGT-TCAGAAG59.8169R:TGGTCATCAGTATCCTC-CCAAC59.7

注：正向引物(F)forwards primer，反向引物(R)reverse primer，解链温度(Tm)melting temperature。

1.2.4 候选内参基因的普通PCR 利用表2中的引物对6个候选内参基因进行普通PCR，体系为2×Taq PCR MasterMix 10 μL，混合cDNA 2.0 μL，正向和反向引物(浓度10 μmol·L-1)各1.0 μL，ddH2O补充体积至20 μL。反应温度程序为：95℃预变性2 min；94℃变性30 s，60℃退火30 s，72℃延伸30 s，45个循环；72℃延伸10 min。反应完成后用2%琼脂糖凝胶电泳检测各PCR产物。

2.4.1 GeNorm软件评价结果 GeNorm软件程序采用几何平均算法，通过计算出每个内参基因稳定性的M值来筛选出稳定性较好的内参基因，M值越小内参基因的稳定性就越好，反之，稳定性越差。同时，通过归一化因子法计算出最适合的内参基因搭配组合[12,14]，来确定所需内参基因的最适数目，当Vn/(n+1)<1.5时，说明有n个基因可以作为qRT-PCR定量分析的内参基因，当Vn/(n+1)>1.5时，说明有n+1个基因可以作为qRT-PCR定量分析的内参基因。通过GeNorm分析，6个内参基因的稳定性由高到低的顺序为：EF-1α=GAPDH> Actin>TIP-41>PP2A>Tubulin(图3)。根据标准差数据分析(图4)，Vn/(n+1)均大于1.5，其中V5/6的M值最小，可以选择6个内参基因做组合。

1.2.6 数据分析及内参基因稳定性评价 qRT-PCR扩增完成后，采用Ct值的方法[12]计算6个内参基因的表达量。利用内参基因评价软件GeNorm和NormFinder分别对不同基因进行稳定性评价，以确定理想的内参基因。

1.2.7 内参基因有效性验证 以筛选出的基因为内参，对巨龙竹木质素合成中PTAL(Phenylalanine/tyrosine ammonia-lyase)[13]基因的表达量进行定量分析，验证内参基因的有效性。体系和反应程序同1.2.5，实验设3个生物学重复，每个生物学重复设3个技术重复。根据项目组前期采用相同植物材料进行转录组测序获得的结果(待发表)和实时荧光定量分析之间基因表达量的相关性验证内参基因的有效性。

1.3.1 比较联合组与对照组患者治疗效果。显效:患者在治疗后舒张压降低≥10mm Hg并降至正常范围,或者降低≥20mm Hg;有效:患者在治疗后舒张压降低<10mm Hg以内但血压处于正常值,或者降低10至20mm Hg;无效:患者在治疗后血压降幅未达以上标准。治疗总有效率为显效率、有效率之和。

2 结果与分析

2.1 总RNA样品的质量分析

对于不同样品或组织中基因表达的精确定量分析，选择合适的内参基因是关键，GeNorm软件和NormFinder软件是专业用来筛选内参基因的软件。

3）果园悬挂糖醋液，诱杀苹小食心虫，可减轻其危害。糖醋液配制比例：糖0.5 kg、醋1 kg、水10 kg。

  

图1 巨龙竹不同秆形各发育时期总RNA电泳图Fig.1 Electrophoresis results of total RNA extracted from different types of D.sinicus at three developmental stages

2.2 内参基因的普通PCR分析

将苦杏仁苷标准溶液与 20%甲醇溶液混合，在190～400 nm波段下扫描。该溶液在210 nm波长附近有最大吸收，确定210 nm为苦杏仁苷的检测波长。

我有一个朋友，年富力强，精力充沛。他除了在我们当地有一份工作，在外地还有一份跨界领域的工作，年头到年尾，两地跑，两头忙。由于他年年出色完成“两下里”的工作任务，便理所当然地享受双份报酬，买车便买车，买房便买房，没有这“奴”那“奴”的困厄和压力。我深为他过人的才气和不俗的业绩所折服，可又“眼前有景道不得”，一直不知道该怎样来概括他的这种进击状态，每每用“你活成了两个人”由衷地称赞他。

一是全面实施绩效管理。财政部门制定全面实施绩效管理意见，在预算编制、执行和监督环节全面引入绩效管理理念及要求。扩大绩效评价范围，健全绩效评价指标库，完善财政预算管理工作绩效评价机制，将部门预算管理评价结果纳入省级机关绩效考核范围，实现财政管理评价结果与相关支持政策挂钩。

  

图2 6个候选内参基因常规PCR的电泳结果Fig.2 Electrophoresis results of six candidateerence genes

2.3 标准曲线与熔解曲线的分析

2.4.2 NormFinder软件评价结果 NormFinder软件采用方差估计算法，先获得内参基因表达稳定值，再根据稳定值(M值)大小来筛选最适合的内参基因，M值最小的内参基因为最合适的内参基因[14-15]。通过NormFinder分析，6个内参基因的稳定性由高到低的顺序为：Actin> GAPDH>EF-1α> PP2A>Tubulin> TIP-41。

2.4 内参基因稳定性评价

提取的总RNA经2%琼脂糖凝胶电泳结果显示如图1，所提总RNA有28S、18S两条清晰的条带，条带无拖尾，RNA完整。经Nanodrop 8000紫外分光检测仪检测，A260/280值在1.72.0范围内，A260/230值大于2.0，RNA的纯度高，表明所制备的总RNA质量较好，满足试验需求。

1.2.5 候选内参基因的qRT-PCR分析 qRT-PCR采用20 μL的反应体系：2×TransStart Top Green qPCR SuperMix 10 μL，各样品的cDNA 1.0 μL(阴性对照加ddH2O 1.0 μL)，正向和反向引物(浓度10 μmol·L-1)各0.4 μL，用ddH2O补足体积至20 μL。反应程序为：94℃预变性30 s；94℃变性5 s，60℃退火30 s，40个循环；6595℃溶解曲线分析。实验设计3个技术重复。加样完成后迅速放入BioRad CFX 96 real-time PCR仪上进行反应。数据读取由仪器自动完成。

常规PCR结果显示，各内参基因均可见预期大小的目标条带(如图2)，且条带单一，无非特异性扩增，表明引物特异性良好，可用于后续qRT-PCR实验。

  

图3 6个候选内参基因的表达稳定性GeNorm分析Fig.3 Expression stabilities of six candidate reference genes analyzed by GeNorm

  

图4 6个候选内参基因标准化因子的配对差异分析Fig.4 Pairwise variance analyses of six candidate reference genes

设置5个浓度梯度，3个技术重复，标准曲线相关系数R2均大于0.99，说明了线性关系好。扩增效率E在90%105%。熔解曲线有单一峰值，峰线窄而尖，说明引物特异性好，扩增产物单一。

《国家中长期教育改革和发展规划纲要（2010-2020年）》指出：“重点发展农村学前教育。努力提高农村学前教育普及程度。”2018年5月15日颁布的《山东省政府办公厅关于加快学前教育改革发展的意见》指出：“实施公费师范生乡村幼儿教师培养计划，按照公费师范生培养政策，2018～2020年每年招收1500名公费师范生，为乡镇中心幼儿园培养高层次幼儿教师。”由此看来，农村学前教育日益受到国家的重视。

  

图5 6个候选内参基因的表达稳定性NormFinder分析Fig.5 Expression stabilities of six candidate reference genes analyzed by NormFinder

综上所述：通过使用GeNorm与NormFinder 2种软件对6个候选内参基因的稳定性评估结果得知GAPDH、Actin、EF-1α的稳定性最佳，可作为研究巨龙竹秆形发育过程基因表达分析使用的内参基因。

2.5 内参基因应用的验证

为了验证筛选的内参基因在巨龙竹茎秆发育过程中应用的有效性，选取木质素合成过程中的酶基因PTAL，分别以Actin、EF-1α和GAPDH为内参，对两种秆型巨龙竹不同发育时期的材料中的PTAL基因进行定量分析。结果如图6所示，以Actin定量的PTAL基因表达量最高，GAPDH定量的PTAL基因表达量最低。根据前期转录组测序的结果，在弯曲型巨龙竹中，PTAL基因的FPKM值在发育5天、15天、30天的材料中分别为24.15、73.52和1 029.38；在通直型巨龙竹中分别为245.31、53.82和363.4。以不同内参基因获得的表达量与转录组测序获得的表达量进行相关性分析，相关性从大到小为Actin>EF-1α>GAPDH，相关系数分别为0.99，0.95和0.92，说明筛选出的3个内参基因有效性均较高。

  

注：以筛选出的3个内参基因Actin、EF-1α和GAPDH，对两种秆形巨龙竹3个不同发育时期笋中PTAL基因定量的结果，FPKM表示转录组测序结果中PATL基因在不同样品中的表达量。Note: PTAL expression of two types of D.sinicus at three develo pment stages were quantitated by three reference genes Actin、EF-1α and GAPDH. FPKM represents PATL expression level in different samples from the transcriptome sequencing result.图6 不同内参基因对PTAL基因的定量分析Fig.6 Quantitative analysis of different reference genes to PTAL

3 讨论

近年来，随着测序技术的快速突破，“基因组学”的研究越来越便捷，极大地促进了转录组等技术手段的迅猛发展，使之成为探索基因调控网络研究的得力工具。目前，荧光定量PCR(qRT-PCR)技术由于其更高的敏感性和特异性，已经成为基因表达深入研究的基本方法[6]。而合适稳定表达的内参基因的筛选将决定表达量数据的可靠性。本研究对不同变异类型、不同发育时期的巨龙竹竹笋样品进行内参基因筛选的结果对解析“通直型”和“弯曲型”巨龙竹差异发育的分子机理有重要意义。

在qRT-PCR反应中，合适内参基因的选择是决定目的基因精确定量分析的先决条件。理想的内参基因应在生物体不同的生长发育阶段、组织及各种外界因素条件下均恒定表达。然而，内参基因的适用性并不是普遍的，需要从多个内参基因中筛选出表达相对稳定的一个或一个以上作为研究对象的内参基因[15]。Zeng等[16]在对枸杞属两个种Lycium barbarum L.和L. ruthenicum Murray内参基因的筛选中，确定在不同的组织细胞中，表达最稳定的内参基因并不统一。对牡丹不同发育阶段种子和花瓣进行qRT-PCR内参基因筛选的结果也显示，在不同发育阶段及不同组织，最稳定的内参基因也不同[17]。而相同内参基因在禾本科的大麦、小麦和燕麦间的通用性研究结果表明，GAPDH、TUBB、18S rRNA在三者中都表达稳定，而EIF4A和TUBA分别在小麦和燕麦中表达稳定，说明并不是所有内参基因在亲缘关系相近的物种中都通用[18]。在竹亚科植物中，相关内参基因筛选的研究十分有限，而且并非所有的竹亚科植物都能通用一种内参基因。对毛竹成花基因PheTFL1表达分析的研究中，选择了PP2A作为毛竹PheTFL1基因定量分析的内参基因[9]。Fan[10]等从14个候选内参基因中选择TIP41和核苷酸通道结合蛋白基因(NTB)作为毛竹开花前后表达量分析的内参基因。

本研究中，选择EF-1α、GAPDH、Actin、Tubulin、TIP-41、PP2A作为巨龙竹的候选内参基因，其中EF-1α、GAPDH、Actin、Tubilin均为植物学研究中经典的内参基因。EF-1α主要参与快速增殖细胞的生长与增殖[19]，GAPDH、Actin和Tubilin均与细胞骨架结构的形成有关[20-22]。PP2A和TIP-41是新挖掘的内参基因，PP2A主要参与细胞功能的调节[23]，TIP-41主要参与细胞水分的转运[24]。这些内参基因已在多种植物研究中证明了其有效性[9-10, 18]。

4 结论

通过候选内参基因的稳定性和有效性的验证分析与评估，在不同变异类型的巨龙竹笋期及其不同发育时期，表达稳定性最佳的内参基因是Actin、EF-1α、GAPDH。下一步将结合巨龙竹秆形发育转录组学的分析，利用筛选出的内参基因对后续筛选的目的基因进行定量分析，研究巨龙竹秆型发育过程中相关基因的表达模式。

参考文献：

[1] 耿伯介，王正平．中国植物志[M]．北京：科学出版社，1996．

[2] 谷志佳，杨汉奇，孙茂盛，等．巨龙竹资源分布特点及其开花结实现象[J]．林业科学研究，2012，25(1)：1-5．

[3] 辉朝茂，张国学，李在留，等．珍稀竹种巨龙竹种群特性及其保护研究[J]．竹子研究会刊，2004，23(4)：4-9．

[4] 杨汉奇，孙茂盛，阮桢媛，等．云南4种典型热带丛生竹种源分化[J]．林业科学研究，2014，27(2)：172-176．

[5] Bustin S A，Beaulieu J F，Huggett J，et al.MIQE précis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments[J]．BMC Molecular Biology，2010，11(1)：1-5．

[6] Chi X，Hu R，Yang Q，et al.Validation of reference genes for gene expression studies in peanut by quantitative real-time RT-PCR[J]．Molecular Genetics and Genomics，2012，287(2)：167-176．

[7] 孙美莲，王云生，杨冬青，等．茶树实时荧光定量PCR分析中内参基因的选择[J]．植物学报，2010，45(5)：579-587．

[8] 袁 伟，万红建，杨悦俭．植物实时荧光定量PCR内参基因的特点及选择[J]．植物学报，2012，47(4)：427-436．

[9] 齐飞燕，胡 陶，彭镇华，等．毛竹实时荧光定量PCR内参基因的筛选及成花基因PheTFL1表达分析[J]．西北植物学报，2013，33(1)：48-52．

[10] Fan C， Ma J，Guo Q，et al.Selection of Reference Genes for Quantitative Real-Time PCR in Bamboo (Phyllostachys edulis)[J]．Plos One，2013，8(2)：e56573．

[11] Cui K，Wang H，Liao S，et al.Transcriptome Sequencing and Analysis for Culm Elongation of the World's Largest Bamboo (Dendrocalamussinicus)[J]．Plos One，2016，11(6)：e0157362．

[12] 魏洁书，杨锦芬．应用荧光定量比较Ct值法测定基因相对表达量[J]．中国科技论文在线精品论文，2013，6(5)：390-395．

[13] Fraser C M，Chapple C．The Phenylpropanoid Pathway in Arabidopsis[J]．American Society of Plant Biologists，2011，9(1)：e0152．

[14] 吴建阳，何 冰，杜玉洁，等．利用geNorm、NormFinder和BestKeeper软件进行内参基因稳定性分析的方法[J]．现代农业科技，2017，(5)：278-281．

[15] Andersen C L，Jensen J L，rntoft T F．Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets[J]．Cancer Research，2004，64(15)：5245-5250．

[16] Zeng S，Liu Y，Wu M，et al.Identification and validation of reference genes for quantitative real-time PCR normalization and its applications in Lycium[J]．Plos One，2014，9(5)：e97039．

[17] 刘洪峰，高乐旋，胡永红．牡丹不同发育阶段种子和花瓣组织实时荧光定量PCR中内参基因的挖掘与筛选[J]．农业生物技术学报，2015，23(12)：1639-1648．

[18] Jaroová J，Kundu J K．Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR[J]．BMC Plant Biology，2010，10(1)：1-9．

[19] 覃迎姿，黄先益，叶兴枝，等．植物延伸因子eEF1A研究进展[J]．南方农业学报，2009，40(5)：472-477．

[20] 李 洁，谢文光，沈倍奋，等．GAPDH功能多样性研究进展[J]．军事医学，2006，30(5)：483-486．

[21] 郭景康，陈青云，戢 茜，等．拟南芥、水稻和杨树ACTIN家族全基因组分析[J]．上海大学学报:自然科学版，2009，18(4)：426-431．

[22] 董 辉，李越中，胡 玮．γ-微管蛋白研究进展[J]．生物化学与生物物理进展，2002，29(5)：686-690．

[23] Chao Z，Yan W，Fu M，et al.Molecular cloning, expression and functional analysis of three subunits of protein phosphatase 2A (PP2A) from black tiger shrimps (Penaeus monodon)[J]．Comparative Biochemistry & Physiology Part B Biochemistry & Molecular Biology，2017，204(1)：77-89．

[24] 邢朝斌，刘 岩，周 秘，等．刺五加液泡膜内在蛋白基因的克隆与表达分析[J]．中草药，2014，45(2)：250-254．