不同组织及干旱胁迫下黄薇内参基因的筛选与验证

赵 雨，林 琳，王 群，张国哲，王 杰，尚林雪，洪思丹，马清清，顾翠花

（1.浙江农林大学 风景园林与建筑学院，浙江 杭州 311300；2.浙江农林大学 浙江省园林植物种质创新与利用重点实验室，浙江 杭州 311300；3.浙江农林大学 南方园林植物种质创新与利用国家林业和草原局重点实验室，浙江 杭州 311300）

黄薇Heimiamyrtifolia是千屈菜科Lythraceae 黄薇属Heimia落叶丛生灌木，原产于巴西，南亚、东亚地区有零星分布，在中国的上海、浙江、广西等地也有引种栽培[1]。黄薇喜温暖湿润且阳光充沛的环境，不耐干旱，夏季开花，花色金黄，花量丰富，可作为园林绿篱、花坛花境背景应用，点缀和美化植物景观，还可以净化城市水质，具有较高的园林价值和应用前景[2]。目前，国内外主要在黄薇的分类和应用[3]、叶内酚类化合物[4]、类黄酮组分[5]、非生物胁迫的抵御机制[6-8]、叶绿体基因组分析[9]等方面进行了研究，但关于黄薇的分子生物学方面的研究则十分匮乏。

实时荧光定量PCR (RT-qPCR)具有灵敏度高、重复性好、特异性强及高通量等优点，近年来已成为研究基因表达分析的常用方法[10]，在分子生物学领域，RT-qPCR 对深入研究植物基因表达调控机制、发现新基因、细胞因子表达分析、预测基因功能、医学诊断以及食品药品检测等至关重要[11]。在实际的应用过程中，为了确保RT-qPCR 结果的准确性，需要引用一个表达量稳定的内参基因作为参照，而内参基因的稳定性也会影响RT-qPCR 结果的准确性[12-13]。因此，筛选到适宜的内参基因是进行RT-qPCR 表达分析时的重要前提[14]。目前暂未见黄薇内参基因筛选的相关研究，这严重制约了对黄薇目标基因的定量分析。因此，本研究基于黄薇转录组数据，对不同组织及不同干旱胁迫下黄薇的内参基因进行研究，为提高不同组织及不同干旱胁迫下基因RT-qPCR 分析的准确性提供依据。

1 材料与方法

1.1 材料

材料取自浙江农林大学风景园林与建筑学院人工气候室内生长健壮、长势一致的黄薇1年生扦插苗。采样时选取黄薇盛花期的根、茎、叶、花，以及模拟自然条件下的对照组 (土壤相对含水量为65%～75%，ck) 和4 个干旱胁迫处理组 (T1、T2、T3、T4土壤相对含水量分别为45%～60%、30%～45%、15%～30%、10%～15%) 的成熟叶片。采样后立即将样品置于液氮中速冻，并保存在-80 ℃冰箱备用。

1.2 总RNA 提取及cDNA 合成

样品总RNA 的提取按诺唯赞公司的RNA isolater Total RNA Extraction Reagent (RC401-01) 试剂盒使用方法进行。利用质量分数为1.2%的琼脂糖凝胶电泳检测RNA 的完整性，使用紫外分光光度计测定浓度。总RNA 纯度D(260)/D(280) 为1.8～2.4，D(260)/D(230) 为1.5～2.4。使用Reverse Transcriptase MMLV 试剂(Takara)合成cDNA 的第1 链，反转录程序为：37 ℃，15 min (反转录过程)；85 ℃，5 s (反转录酶的失活过程)，4 ℃冰箱保存。最后将cDNA 产物于-20 ℃储存备用。

1.3 候选基因的引物设计

基于黄薇转录组数据(序列号：PRJNA804698)，选取9 个基因作为候选内参基因：延伸因子1α 蛋白基因(elongationfactors-1α,EF-1α)、α-微管蛋白基因(α-tubulin,TUA)、亲环蛋白基因(cyclophilin,CYP)、甘油醛-3-磷酸脱氢酶(glyceraldehyde-3-phosphat dehydrogenase,GAPDH)、18S 核糖体RNA(18S ribosomal RNA,18SRNA)、泛素化酶基因(ubiquitin-conjugating enzyme,UBC)、β-维管蛋白基因(tubulin,TUB)、肌动蛋白基因(actin,ACT)和伴侣蛋白基因(chaperone protein,DNAJ)。基于转录组数据，选择与纤维素合成和干旱胁迫相关的差异表达基因CSLD、SOD基因，验证内参基因对基因表达量的影响。利用Primer 5 软件，设计候选内参基因的荧光定量PCR 引物，并由杭州有康生物技术有限公司合成。

1.4 实时荧光定量PCR

以所有试样的cDNA 为模板依次稀释5 个浓度梯度(5-1、5-2、5-3、5-4、5-5)，设3 次重复，制作9 个内参基因的标准曲线。反转录产物cDNA 稀释10 倍，RT-qPCR 反应使用诺唯赞公司出品的SYBR®Primix ExTaqTM 试剂盒(Q712)。荧光定量反应在Roche LightCycler 480Ⅱ荧光定量仪上进行，反应体系为20 μL：2×TaqPro Universal SYBR qPCR Master Mix 荧光染料10 μL、cDNA 模板2 μL、上下游引物各0.8 μL、ddH2O 6.4 μL。PCR 扩增反应程序为：预变性95 ℃，10 s；PCR 反应95 ℃，10 s，60 ℃，30 s，共40 个循环。采集溶解曲线荧光信号进行分析，反应体系为：从60 ℃到95 ℃，15 s；60 ℃，60 s；95 ℃，15 s。

1.5 数据分析

使用Excel 2016 统计基因的Ct值，用Origin 9.1 制图。采用geNorm[15]、NormFinder[16]、BestKeeper[17]软件以及在线分析软件RefFinder (https://www.heartcure.com.au/reffinder/)，综合分析不同情况下最稳定的内参基因。

2 结果与分析

2.1 RNA 的质量检测和引物特异性验证

所得RNA 产物D(260)/D(280)为1.8～2.2，D(260)/D(230)为1.8～2.4，电泳条带清晰(图1A)，说明提取的总RNA 质量良好，无明显降解，满足后续实验要求。9 个候选内参基因琼脂糖凝胶电泳检测后出现与预期大小一致的单一条带(图1B)，说明引物不存在二聚体，可进行特异性扩增。根据结果绘制标准曲线并计算斜率，如表1显示：候选基因的线性相关系数R2＞0.990，引物扩增效率为95.6%～109.7%。综上，候选内参基因的扩增效率、产物特异性都符合RT-qPCR 的条件，可用于后续的分析。

表1 黄薇候选内参基因引物信息Table 1 Primer information of candidate reference gene in H. myrtifolia

图1 黄薇总RNA 琼脂糖凝胶电泳分析(A)和RT-qPCR 扩增产物特异性(B)Figure 1 Agarose gel electrophoresis analysis of total RNA (A) and specificity of products amplified by RT-qPCR (B) in H. myrtifolia

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

根据Ct表示不同条件下基因的相对表达量，Ct越小，表达丰度越高[18]。所有样品中，9 个候选内参基因的Ct为21.12～33.74 (图2)，其中平均表达量最高的基因为GAPDH，Ct平均值为25.01；18SRNA的Ct平均为29.90，其表达量最低。EF-1α基因的表达水平变化最小，Ct为23.97～29.06；18SRNA 的Ct为26.64～33.47，表达变化最大。从Ct的大小和变化情况来看，GAPDH基因的表达相对稳定且丰度最高，可初步评价候选内参基因的表达水平和稳定性。

图2 9 个候选内参基因的Ct 值Figure 2 Ct values of the 9 candidate reference genes

2.3 geNorm 分析

geNorm 软件根据平均稳定值M来判定结果，M越小，基因表达越稳定，反之越不稳定[19]。使用2-ΔΔCt法计算内参基因的表达量：在黄薇不同组织中，候选内参基因的表达稳定性从高到低依次为CYP=GAPDH、18SRNA、UBC、DNAJ、ACT、EF-1α、TUB、TUA(图3A)。在全部样品和干旱胁迫中，黄薇GAPDH和UBC基因最稳定，TUB基因最不稳定(图3B 和图3C)。geNorm 软件还可以根据成对变异值，分析最适的内参基因数目，软件默认的配对差异值(Vn/Vn+1，V为配对变异指数，n为可使RT-qPCR 结果准确的最少基因数目)为0.15[20]。当Vn/Vn+1＞0.15，需要加入第n+1 基因；当Vn/Vn+1＜0.15，不需要引入新的基因。黄薇在不同组织和干旱胁迫下(图3D)，除不同组织中V2/V3＜0.15 外，其余情况均大于0.15。该现象可能是因为在不同组织和胁迫处理下，黄薇基因的表达量变化差异大，但软件说明书上也提出Vn/Vn+1最大值设定为0.15 并非总是必须的[21]。

图3 geNorm 分析9 个候选内参基因的表达稳定值和配对变异值Figure 3 Stable expression values and paired variation values of nine candidate reference genes were analyzed by geNorm

2.4 NormFinder 分析

NormFinder 软件按照2-ΔΔCt法计算稳定性，其中表达稳定值越小，候选内参基因越稳定[22]。结果表明(表2)：在不同组织中，GAPDH和CYP基因表达最稳定，TUA基因最不稳定。在不同干旱胁迫处理下，GAPDH基因是最稳定的内参基因，TUB基因最不稳定。在黄薇的全部样品中，最稳定的内参基因是GAPDH基因，这与geNorm 软件得出的结果相同。

表2 NormFinder 分析内参基因的表达稳定值Table 2 Expression stability values of nine candidate reference genes calculated by the NormFinder

2.5 BestKeeper 分析

BestKeeper 软件以内参基因的标准差(SD)和变异系数(CV)为基准，其数值越小越稳定[23]。从表3可见：在不同组织中，GAPDH基因的CV±SD 为2.65 ± 0.69，排名第一，最稳定；在干旱胁迫和全部样品中，EF-1α基因的CV±SD 分别为1.97±0.53 和2.58±0.71，表达最稳定。综合不同条件下的分析结果，EF-1α基因被鉴定为最稳定的内参基因。

2.6 RefFinder 综合分析

上述3 种软件的结果存在一定的差异，使用在线分析软件RefFinder，整合上述分析结果并进行综合分析。由表4显示：不同组织中，最稳定的内参基因是GAPDH和CYP，最不稳定的内参基因是TUA；在干旱胁迫下，最稳定的内参基因是GAPDH和TUA，最不稳定的内参基因是TUB。在全部样本中，最稳定的内参基因为GAPDH，最不稳定的内参基因为18SRNA。

表4 RefFinder 分析内参基因的表达稳定值Table 4 Expression stability of candidate reference genes ranked by RefFinder

2.7 内参基因稳定性的验证

为验证候选内参基因的稳定性，采用RT-qPCR 技术，根据筛选的稳定内参基因，对干旱胁迫下的靶基因CSLD和SOD的表达模式进行验证。CSLD基因是一种纤维素合成酶的结构蛋白，在茎中的表达量一般较高[24]。SOD基因属于抗氧化酶功能基因，在植物遭遇干旱胁迫时上调表达以抵御外界伤害[25]。从图4可见：在不同组织中，GAPDH和CYP基因在单独和组合校正下，CSLD均在茎中表达量最高，在其他组织中表达变化不大。但使用TUA作内参基因时，CSLD在茎中的表达量最低，在叶片中表达却最高。干旱胁迫下，使用GAPDH和基因组合GAPDH+TUA时，SOD基因在T3处理表达量最高，表达趋势大致相同；而利用TUB基因进行验证时，SOD基因在T2处理表达量最高，与其他结果差异较大。以上结果进一步验证了GAPDH基因作为内参基因的稳定性与可靠性。

图4 不同组织和干旱胁迫下CSLD (A)和SOD 基因(B)利用不同内参基因的表达分析Figure 4 Expression analysis for CSLD gene (A) and SOD gene (B) using different reference genes in different tissues and drought stress

3 讨论与结论

RT-qPCR 技术是一种分析基因表达量的研究工具，而使用相对定量法检测的可靠性，极大依赖内参基因的稳定性[26]。理想的内参基因在不同条件下，基因的表达变化较小，常用的内参基因一般都是维持植物生命活动必须的看家基因[27]。本研究基于黄薇转录组数据，筛选了9 个常用内参基因(EF-1α、TUA、CYP、GAPDH、18SRNA、UBC、TUB、ACT、DNAJ)，其中，GAPDH在不同组织和干旱胁迫下都较稳定，最不稳定的是18SRNA。geNorm 和NormFinder 结果中，得出的最适内参基因组合为GAPDH+UBC，BestKeeper 得出最佳基因为EF-1α，但GAPDH的稳定值排名综合也较靠前。3 个筛选软件得出的结果基本一致，但也存在一定的差异，这可能是由于算法的不同所导致的[28]。与本研究结果不同的是，张海洋等[29]探究了菠菜Spinacia oleracea在不同胁迫处理下稳定的内参基因，3 个软件的结果一致，在氯化钠和高温胁迫下，G6PD的表达稳定性最好，聚乙二醇胁迫下ELF1B的表达稳定性最好。王蕊等[30]研究了不同发育时期大豆Glycinemax在不同组织、不同非生物胁迫下稳定表达的最适内参基因，4 个软件在全部组织及全部胁迫中综合分析结果均一致，在全部组织中最适内参基因为ACT，全部胁迫中最适内参基因为EF-1α。

非生物胁迫会改变植物正常的代谢过程，筛选在胁迫处理下稳定表达的内参基因，有利于后续开展目的基因的鉴定，以便进一步了解目的基因的功能[31]。目前已有大量研究筛选不同非生物胁迫下的最佳内参基因。如紫鸭跖草Tradescantiapallida在不同浓度的铜离子胁迫下，最佳内参基因组合为18SRNA +TUB+UBI[32]。马铃薯Solanumtuberosum在干旱和渗透胁迫条件下，EF-1α和sec3 基因表达最稳定[33]。上述结果与本研究的结果存在一定差异，说明不同植物在不同非生物胁迫下最适内参基因也在变化，特定条件下的最适内参基因仍需通过试验筛选。