大白菜-结球甘蓝易位系实时荧光定量 PCR内参基因的筛选

通过实时荧光定量PCR(Quantitative Real-time PCR，qRT-PCR)技术分析基因表达及调控，对于解析植物各种复杂的生物代谢途径和遗传性状及逆境胁迫生理响应具有重要的意义[1]。在利用qRT-PCR技术中，为了消除或降低不同组织间因RNA的提取质量、初始模板量或者反转录效率等所造成的偏差，必须选择适宜的内参基因进行校正和标准化[2]。选择稳定的内参基因也是保证其分析结果准确可靠的重要前提条件[3-4]。目前，Actin、GAPDH、18S rRNA、Ubiquitin等内参基因已经被广泛应用于qRT-PCR分析。然而，近年来大量的研究结果表明，任何一个内参基因不能够适用于所有条件下的试验研究[5-8]。所谓的持续稳定表达都只是限定在一定类型的细胞或试验因素作用下“有范围”的恒定[9-12]。据 Gutierrez等[13]报道，使用未经验证的内参基因，目标基因表达水平的定量会出现100倍的偏差。因此，在利用qRT-PCR进行定量分析之前，在不同试验材料或不同试验条件中，有必要验证、筛出稳定的内参基因，作为每个试验研究所适用的最佳内参基因稳定的定量因子。

本研究利用geNorm[14]和NormFinder[15]统计软件，比较分析了21个候选内参基因在大白菜-结球甘蓝易位系营养生长阶段的幼苗期、莲座期、包心期、结球期，及包心期IAA和TIBA处理的叶片和生殖生长阶段6个大小不同等级的花蕾的内参基因表达稳定性，旨在为精确分析大白菜-结球甘蓝易位系的基因表达提供保障，同时也为芸薹属植物不同发育时期及激素处理对内参基因的选择提供参考。

1 材料和方法

1.1 试验材料

河北省蔬菜种质创新与利用重点实验室创制的添加结球甘蓝(Brassica oleracea var. capitata)4号染色体片段、基因型纯合的大白菜(Brassica rapa ssp. pekinensis)易位系B3-29，于2015年8月8日播种于温室，进行育苗，9月12日进行田间定植。分别在幼苗期、莲座期、包心期和结球期的幼嫩叶片进行取材。并在包心期分别喷施生长素IAA(0.2，0.6，1.2 mmol/L)和生长素抑制剂TIBA(0.5，1.0，1.5 mmol/L)，处理后12 h进行幼嫩叶片取材。2016年春季，待其抽薹开花取花蕾，分为6个阶段，花发育阶段1(花蕾横径介于1～2 mm)、花发育阶段2(花蕾横径介于2 ～3 mm)、花发育阶段3(花蕾横径大于3 mm，长度大于4 mm)、花发育阶段4(花蕾横径大于4 mm，长度大于5 mm)、花发育阶段5(花冠露出尚未伸展)、花发育阶段6(开花当天)。每个样品3个重复。所有样品液氮速冻后保存于-80 ℃冰箱备用。

1.2 实时荧光定量 PCR

本研究共检测了21个候选内参基因的表达稳定性：腺嘌呤转磷酸核糖基酶基因(Adeninephosphoribosyl transferase，Apr)、液泡膜内在蛋白基因(Tonop last intrinsic proteins，BcTIP41)、未知功能基因(Unknown protein，U34559)、延伸因子1基因(Elongationfactor1α，EF1α)、微管蛋白基因(Tubulin beta，TUB4)、亲环蛋白基因(Cyclophilin，CYP)、锌指蛋白基因(DnaJ-Like，DNAJ)、组蛋白基因(Histone h3，HIS)、α-微管蛋白基因(Alpha tubulin 5，TUA5)、假设蛋白基因(Hypothetical，UKN1)、互作蛋白基因(Ask-interacting protein 16，SKIP16)、网格衔接蛋白复合体基因(Clathrin adaptor complex，CAC)、肌动蛋白基因3个同源基因(ACTIN、ACTIN-1、ACTIN-2)、甘油醛-3-磷酸-脱氢酶基因(Glyc-eraldehyde-3-phosphate dehydrogenase，GAPDH)、泛素结合酶基因(Ubiquitin-conjugating enzyme，UBC30)、泛素连接酶基因(Poly ubiquitinenzyme，UBQ)、三角状五肽重复结构域基因(Pentatricopeptide repeat，PPR)、蛋白磷酸酶2A基因(Protein phosphatase 2A gene，PP2A)、苹果酸脱氢酶基因(Malate dehydrogenase，MDH)所有引物序列见表1。将样品的cDNA进行5个浓度梯度稀释，优化反应条件，确定检测样品的稀释浓度。qRT-PCR反应体系(20 μL)：cDNA模板 2 μL，正反向引物(50 ng/μL)各0.5 μL，ddH2O 7 μL，THUN-DERBIRD SYBR qPCR Mix(QPS-201(-)TOYOBO公司)10 μL。反应程序：98 ℃预变性3 min；98 ℃变性30 s，58 ℃退火10 s，68 ℃延伸10 s，45个循环；之后65～95 ℃熔解曲线分析。每个样品均3次重复。试验使用罗氏LightCycler96实时荧光定量PCR仪定量分析系统进行数据采集。本试验所用试剂IAA和TIBA购于SIGMA公司，纯度大于98%。用RNA提取试剂盒ESAYspin Plus(北京艾德莱公司)分别提取样品总RNA。反转录根据ReverTra Ace qPCR RT Master Mix 试剂盒(东洋纺生物科技有限公司)说明进行cDNA合成。

1.3 数据处理和分析

采用常用的内参分析软件geNorm(ver. 3.5)和NormFinder(ver. 0.953)进行数据处理和分析。根据qRT-PCR结果中Ct值，利用公式Q=E-ΔCt，计算出各基因的相对表达量Q[20]。利用geNorm软件对各候选内参基因的相对表达量Q进行表达稳定性统计学分析。其中，M值越小，稳定性越高。利用标准化因子的配对差异分析(V)，确定合适的内参基因数目。该程序默认V值为0.15，若Vn/n+1>0.15，则引入第n+1个内参基因[14]。NormFinder软件通过方差分析得出基因的表达稳定值M，直接评价基因的稳定性[15]。

2 结果与分析

2.1 候选内参基因表达水平

对于21个候选内参基因，通过计算Ct值来比较基因的表达丰度，结果表明，试验涉及的21个内参基因在大白菜-结球甘蓝易位系不同类型材料中的平均Ct值为17.28～31.98(表2，3)，并且大部分Ct值都集中在20.00～25.00，可用于后续试验。其中，GAPDH和MDH表达量较高，Ct值分别在19.00和20.00左右，PPR表达量较低，Ct值在30左右。从所有材料中的21个内参基因平均Ct值的箱式图中可以看出(图1)，GAPDH和SKIP16基因的平均Ct值的最大值与最小值浮动范围较大。

2.2 候选内参基因表达稳定性分析

2.2.1 geNorm 分析通过geNorm软件计算在大白菜-结球甘蓝易位系不同样品中各候选内参基因表达的M值，进行稳定性分析。结果显示(表4)，在大白菜-结球甘蓝易位系营养生长阶段和生殖生长的不同时期或是在IAA和TIBA激素处理后，各候选内参基因表达稳定性有明显差异。在大白菜-结球甘蓝易位系营养生长阶段幼苗期、莲座期、包心期和结球期UKN1和TUB4基因表达较为稳定。在大白菜-结球甘蓝易位系包心期用不同浓度生长素IAA处理后，均表现出最稳定的内参基因为BcTIP41和ACTIN。用不同浓度生长素抑制剂TIBA处理后，最稳定的内参基因均为UKN1和BcTIP41。生殖生长时期花发育阶段的6个等级大小花蕾中最稳定内参基因为DNAJ和ACTIN，其次为PP2A。在总样品分析中认为，内参基因为BcTIP41和U34559较为稳定，其次为ACTIN，而内参基因SKIP16、GAPDH和MDH表现不稳定。

为了获得更加准确可靠的RT-qPCR分析结果，通常会使用2个或更多内参基因进行校正。利用geNorm软件进行标准化因子的配对差异分析(V)，确定各组合样品中最佳内参基因数目，由图2结果可知，对总样品分析中，V3/4的值是0.12，小于界限值0.15。因此，在总样品分析中，一组最优的稳定的内参基因至少需要3个，分别为：BcTIP41、U34559和ACTIN。同时花发育阶段组合V3/4的值是0.14，小于界限值0.15，表明在花发育阶段也需选用3个最稳定的内参基因进行标准化，分别为：DNAJ、ACTIN和PP2A。营养生长阶段和激素处理后材料的V2/3的值均小于界限值0.15，故只需选用2个最稳定的内参基因进行标准化。在营养生长阶段(从幼苗到结球)最稳定的内参基因为UKN1和TUB4。在包心期用不同浓度生长素IAA处理后，内参基因BcTIP41和ACTIN最为稳定，用不同浓度生长素抑制剂TIBA处理后，内参基因UKN1和BcTIP41最为稳定。

2.2.2 NormFinder软件分析为了进一步验证geNorm分析的结果，笔者也使用了NormFinder软件进行稳定性分析。2个软件推荐的最稳定和最不稳定的基因比较一致(表4)，如：在2个软件的分析显示，在总样品、营养生长阶段的幼苗期到结球期、生长素抑制剂TIBA处理后推荐的最为稳定的和最不稳定内参基因相一致。而在生长素IAA处理后最稳定内参基因BcTIP41的排名在2个软件存在一些差别，在geNorm软件中与ACTIN基因的M值相同排在第1位，而NormFinder软件分析却排在第4位，但是BcTIP41的M值与排在第2位基因仅相差0.06。在花发育阶段也出现了类似的现象，但是M值之间差距较小。

3 讨论与结论

目前，在多数植物中都鉴定出了各自在不同的试验条件下表达最为稳定的内参基因，包括大豆(Glycine max)、小拟南芥(Arabidopsis pumila)、野生条斑紫菜(Pyropia yezoensis)、牡丹(Paeonia suffruticosa)、菊花(Chrysanthemum lavandulifolium)、鼠尾草(Salvia hispanica)、茄子(Solanum melongena L.)[21-27]。研究结果表明，内参基因GAPDH在葡萄(Vitis amurensis)和荔枝(Litchi chinensis Sonn.)中表达均较稳定，而在黑麦草(Lolium perenne)、芥菜(Brassica juncea)、柳枝稷(Panicum virgatum L.)中的表达稳定性则较差 [28-32]。在菊花中能够稳定表达的SKIP16，在大白菜花的发育中则稳定性较低[18]。李晗等[33]报道，在羽衣甘蓝的不同组织中GAPDH为最不稳定内参基因。本研究中，GAPDH和SKIP16基因在大白菜-结球甘蓝易位系营养生长时期叶球发育过程中及IAA和TIBA处理后都不能稳定表达。表明在不同物种中，同一内参基因的表达稳定性各不相同。

本研究大白菜-结球甘蓝易位系不同发育阶段、不同组织部位以及不同激素处理条件下，21个候选内参基因的表达稳定性也各不相同。Zhang等[34]对麻风树(Jatropha curcas)内参基因分析，认为营养生长阶段ACT和TUB8内参基因最稳定，在生殖生长阶段内参基因GAPDH和EF-1α最为稳定。吝月爱[35]在玉米的不同组织中认为，内参基因TUB、UBQ9和EF1lα是茎中最稳定；EF1α、UBQ9和GAPDH是叶片中最稳定；在激素处理条件下，EF1α、TUB和GAPDH在脱落酸处理条件下最稳定；ACT2和TUB是赤霉素处理条件下最稳定。刘洪峰等[24]在牡丹内参基因筛选时，认为不同组织(茎、叶片、花萼、花瓣)中内参基因AMPDS和PUF1639最为稳定，在不同发育时期的种子中认为RPS9和PUF1639基因最为稳定。进一步证明了同一植物中的不同发育阶段、不同组织部位及不同处理条件下，内参基因的稳定性不尽相同。

Cheng等[18]对13个内参基因在不结球大白菜花的不同发育阶段及不同花器官组织部位中进行内参基因筛选，其中，DNAJ、UKN1和PP2A被确定为最稳定的内参基因。本研究对大白菜-结球甘蓝易位系花发育阶段的6个等级大小的花蕾进行21个内参基因筛选发现，最稳定的内参基因为DNAJ、ACTIN和PP2A，与Cheng等[18]研究结果相似。本研究选择的内参基因在包含了Cheng等[18]研究的12个内参基因基础上，又增加了9个内参基因进行筛选，共21个内参基因，其中，DNAJ和PP2A在不结球大白菜、大白菜-结球甘蓝易位系的花器官中都表现了高稳定性，表明该研究不仅为大白菜-结球甘蓝易位系基因表达量的精准分析提供了必要保障，同时也为芸薹属其他植物内参基因的选择提供了参考。

在大白菜-结球甘蓝易位系营养生长阶段(从幼苗到结球)内参基因UKN1和TUB4表达最为稳定。生长素IAA处理后，最稳定的内参基因为BcTIP41和ACTIN，生长素抑制剂TIBA处理后，最稳定的内参基因为UKN1和BcTIP41。在花发育阶段6个大小等级花蕾中DNAJ、ACTIN和PP2A基因表达稳定。内参基因GAPDH和MDH在大白菜-结球甘蓝易位系的营养生长阶段、激素IAA和TIBA处理及花发育阶段表达均不稳定。

[1] Mascia T，Santovito E，Gallitelli D，et al. Evaluation of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in infected tomato plants[J]. Molecular Plant Pathology，2010，11(6)：805-816.

[2] Derveaux S，Vandesompele J，Hellemans J. How to do successful gene expression analysis using Real-time PCR[J]. Methods，2010，50(4)：227-230.

[3] Kozera B，Rapacz M. Reference genes in Real-time PCR[J]. Journal of Applied Genetics，2013，54(4)：391-406.

[4] Vanguilder H D，Vrana K E，Freeman W M. Twenty-five years of quantitative PCR for gene expression analysis[J]. BioTechniques，2008，44(5)：619-626.

[5] Dong M，Zhang X，Chi X，et al. The validity of a reference gene is highly dependent on the experimental conditions in green alga Ulva linza[J]. Current Genetics，2012，58(1)：13-20.

[6] 林榕燕，钟淮钦，黄敏玲，等.杂交兰实时荧光定量PCR内参基因的筛选[J].中国细胞生物学学报，2018，40(3)：381-389.

[7] Zhou X H，Liu J，Zhuang Y. Selection of appropriate reference genes in eggplant for quantitative gene expression studies under Different experimental conditions[J]. Scientia Horticulturae，2014，176(2)：200-207.

[8] Chen H，Yang Z Q，Hu Y，et al. Reference genes selection for quantitative gene expression studies in Pinus massoniana L.[J]. Trees，2016，30(3)：685-696.

[9] Long X Y，Wang J R，Ouellet T，et al. Genome-wide identification and evaluation of novel internal control genes for Q-PCR based transcript normalization in wheat[J]. Plant Molecular Biology，2010，74(3)：307-311.

[10] Hua W，Zhu J，Shang Y，et al. Identification of suitable reference genes for barley gene expression under abiotic stresses and hormonal treatments[J]. Plant Molecular Biology Reporter，2015，33(4)：1002-1012.

[11] Xu Y，Hui L，Li X，et al. Systematic selection and validation of appropriate reference genes for gene expression studies by quantitative real-time PCR in pear[J]. Acta Physiologiae Plantarum，2015，37(2)：40.

[12] Ma S，Niu H，Liu C，et al. Expression stabilities of candidate reference genes for RT-qPCR under Different stress conditions in soybean[J]. PloS One，2013，8(10)：e75271.

[13] Gutierrez L，Mauriat M，Gué nin S，et al. The lack of a systematic validation of reference genes：a serious pitfall undervalued in reverse transcription-polymerase chain reaction(RT-PCR)analysis in plants[J]. Plant Biotechnology Journal，2008，6(6)：609-618.

[14] Vandesompele J，Preter K D，Pattyn F，et al. Accurate normalization of Real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes[J]. Genome Biology，2002，3(7)：research0034.1.

[15] Andersen C L，Jensen J L，Orntoft 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] Qi J，Yu S，Zhang F，et al. Reference gene selection for Real-time quantitative polymerase chain reaction of mRNA transcript levels in Chinese cabbage(Brassica rapa L. ssp. pekinensis)[J]. Plant Molecular Biology Reporter，2010，28(4)：597-604.

[17] Xu X，Yang Z，Sun X，et al. Selection of reference genes for quantitative Real-time PCR during flower bud development in CMS7311 of heading Chinese cabbage(Brassica rapa L. ssp. pekinensis)[J]. Acta Physiologiae Plantarum，2014，36(3)：809-814.

[18] Cheng W，Cui H M，Huang T H，et al. Identification and validation of reference genes for RT-qPCR analysis in non-heading Chinese cabbage flowers[J]. Frontiers in Plant Science，2016，7(651)：811.

[19] Xiao D，Zhang N W，Zhao J J，et al. Validation of reference genes for Real-time quantitative PCR normalisation in non-heading Chinese cabbage[J]. Functional Plant Biology，2012，39(4)：342-350.

[20] 李钱峰，蒋美艳，于恒秀，等.水稻胚乳RNA定量RT-PCR分析中参照基因选择[J].扬州大学学报：农业与生命科学版，2008，29(2)：61-66.

[21] 侯高杰.铝毒胁迫下大豆内参基因的筛选及相关基因表达分析[D].南宁：广西大学，2015.

[22] 郑丽洁，林军，黄先忠.小拟南芥实时荧光定量PCR(qRT-PCR)内参基因的筛选[J].基因组学与应用生物学，2017(2)：774-783.

[23] Kong F N，Cao M，Sun P P，et al. Selection of reference genes for gene expression normalization in Pyropia yezoensis using quantitative Real-time PCR[J]. Journal of Applied Phycology，2015，27(2)：1003-1010.

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

[25] Shuai Q，Yang L，Wen X，et al. Reference gene selection for RT-qPCR analysis of flower development in Chrysanthemum morifolium and Chrysanthemum lavandulifolium [J]. Frontiers in Plant Science，2016，7(287).

[26] Gopalam R，Rupwate S D，Tumaney A W. Selection and validation of appropriate reference genes for quantitative Real-time PCR analysis in Salvia hispanica[J]. PloS One，2017，12(11)：e0186978.

[27] 庞强强，李植良，罗少波，等.高温胁迫下茄子qRT-PCR内参基因筛选及稳定性分析[J].园艺学报，2017，44(3)：475-486.

[28] Peng F，James S，Niclas O，et al. An optimized grapevine RNA isolation procedure and statistical determination of reference genes for Real-time RT-PCR during berry development[J]. Bmc Plant Biology，2006，6(1)：27.

[29] Zhong H Y，Chen J W，Li C Q，et al. Selection of reliable reference genes for expression studies by reverse transcription quantitative Real-time PCR in litchi under Different experimental conditions[J]. Plant Cell Reports，2011，30(4)：641-653.

[30] Ruby C，Rehna A，Bisht N C. Evaluation of candidate reference genes for gene expression normalization in Brassica juncea using real time quantitative RT-PCR[J]. PloS One，2012，7(5)：e36918.

[31] 严海东，蒋晓梅，张新全，等.非生物胁迫下多年生黑麦草qRT-PCR分析中内参基因的选择[J].农业生物技术学报，2014，22(12)：1494-1501.

[32] 蒋晓梅，张新全，严海东，等.柳枝稷根组织实时定量PCR分析中内参基因的选择[J].农业生物技术学报，2014，22(1)：55-63.

[33] 李晗，李治龙，李晓屿，等.羽衣甘蓝不同组织及柱头发育实时荧光定量PCR内参基因的筛选[J].植物研究，2016，36(4)：565-572.

[34] Zhang L, He L L, Fu Q T, et al. Selection of reliable reference genes for gene expression studies in the biofuel plant Jatropha curcas using Real-time quantitative PCR[J].International Journal of Molecular Sciences, 2013, 14(12):24338.

[35] 吝月爱.玉米在非生物胁迫和激素处理条件下实时荧光定量PCR内参基因的选择[D].雅安：四川农业大学，2012.