-
黄薇Heimia myrtifolia是千屈菜科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年生扦插苗。采样时选取黄薇盛花期的根、茎、叶、花，以及模拟自然条件下的对照组 (土壤相对含水量为65%~75%，ck)和4个干旱胁迫处理组 (T1、T2、T3、T4土壤相对含水量分别为45%~60%、30%~45%、15%~30%、10%~15%)的成熟叶片。采样后立即将样品置于液氮中速冻，并保存在−80 ℃冰箱备用。
-
样品总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 M-MLV试剂(Takara)合成cDNA的第1链，反转录程序为：37 ℃，15 min (反转录过程)；85 ℃，5 s (反转录酶的失活过程)，4 ℃冰箱保存。最后将cDNA产物于−20 ℃储存备用。
-
基于黄薇转录组数据(序列号：PRJNA804698)，选取9个基因作为候选内参基因：延伸因子1α蛋白基因(elongationfactors-1α, EF-1α)、α-微管蛋白基因(α-tubulin, TUA)、亲环蛋白基因(cyclophilin, CYP)、甘油醛-3-磷酸脱氢酶(glyceraldehyde-3-phosphat dehydrogenase, GAPDH)、18S核糖体RNA(18S ribosomal RNA, 18S RNA)、泛素化酶基因(ubiquitin-conjugating enzyme, UBC)、β-维管蛋白基因(tubulin, TUB)、肌动蛋白基因(actin, ACT)和伴侣蛋白基因(chaperone protein, DNA J)。基于转录组数据，选择与纤维素合成和干旱胁迫相关的差异表达基因CSLD、SOD基因，验证内参基因对基因表达量的影响。利用Primer 5软件，设计候选内参基因的荧光定量PCR引物，并由杭州有康生物技术有限公司合成。
-
以所有试样的cDNA为模板依次稀释5个浓度梯度(5−1、5−2、5−3、5−4、5−5)，设3次重复，制作9个内参基因的标准曲线。反转录产物cDNA稀释10倍，RT-qPCR反应使用诺唯赞公司出品的SYBR®Primix Ex TaqTM试剂盒(Q712)。荧光定量反应在Roche LightCycler 480Ⅱ荧光定量仪上进行，反应体系为20 μL：2×Taq Pro 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。
-
使用Excel 2016统计基因的Ct值，用Origin 9.1制图。采用geNorm[15]、NormFinder[16]、BestKeeper[17]软件以及在线分析软件RefFinder (
https://www.heartcure.com.au/reffinder/ )，综合分析不同情况下最稳定的内参基因。 -
所得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的条件，可用于后续的分析。
Figure 1. Agarose gel electrophoresis analysis of total RNA (A) and specificity of products amplified by RT-qPCR (B) in H. myrtifolia
基因 正向引物序列(5′→3′) 反向引物序列(5′→3′) 扩增大小/bp 退火温度/℃ 扩增效率/% R2 EF-1α TGGTTTTGAGGCTGGTATCTCC TTTGCTTGACACCAAGGGTGA 80 56.5 95.6 0.998 TUA TCTCTGCCTTGACCGAGTGA ACCACCAACGGCACTGAAAA 82 56.5 96.6 0.998 CYP ACCCCGACTCGTCCTACAAG TCGGTGTTCCGCTCCAAATG 130 58.0 106.8 0.999 GAPDH AGAAGGTCGTCATTTCTGCCC TGGTTGTGCAGCTAGCGTTG 114 57.5 108.4 0.999 18S RNA CAGGGCCTAGGATTTCGTCC GCCTTCAATCTTAGTCGTGGC 113 58.5 100.5 0.992 UBC GACCTGATGACACTCCCTGG TCACAGTTGGTGGTTTGTTCG 87 57.5 99.1 0.999 TUB GGGTGCTGAGCTTATTGATGC TGAGCAATGTCCCCATGCCT 131 57.5 96.8 0.996 ACT AGGGAATGCCTTTTGATTGATCC AAACATAAGCTCCACTGCCCTC 102 57 109.7 0.999 DNA J CGGAGCTATCACCCCGATG CGGCCTCACCATACCTGTCA 127 59.5 100.2 0.997 CSLD TACCTTGTCCCTTTCGGCG TCAGCGTCCTCATCCCGATA 149 57.4 95.1 0.997 SOD GTTGACGCAAGACGAGGGA CCGTTGGTCGTGTCACCAT 108 57.3 96.0 0.997 Table 1. Primer information of candidate reference gene in H. myrtifolia
-
根据Ct表示不同条件下基因的相对表达量，Ct越小，表达丰度越高[18]。所有样品中，9个候选内参基因的Ct为21.12~33.74 (图2)，其中平均表达量最高的基因为GAPDH，Ct平均值为25.01；18S RNA的Ct平均为29.90，其表达量最低。EF-1α基因的表达水平变化最小，Ct为23.97~29.06；18S RNA的Ct为26.64~33.47，表达变化最大。从Ct的大小和变化情况来看，GAPDH基因的表达相对稳定且丰度最高，可初步评价候选内参基因的表达水平和稳定性。
Figure 2. Ct values of the 9 candidate reference genes
-
geNorm软件根据平均稳定值M来判定结果，M越小，基因表达越稳定，反之越不稳定[19]。使用2−ΔΔCt法计算内参基因的表达量：在黄薇不同组织中，候选内参基因的表达稳定性从高到低依次为CYP=GAPDH、18S RNA、UBC、DNA J、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]。
Figure 3. Stable expression values and paired variation values of nine candidate reference genes were analyzed by geNorm
-
NormFinder软件按照2−ΔΔCt法计算稳定性，其中表达稳定值越小，候选内参基因越稳定[22]。结果表明(表2)：在不同组织中，GAPDH和CYP基因表达最稳定，TUA基因最不稳定。在不同干旱胁迫处理下，GAPDH基因是最稳定的内参基因，TUB基因最不稳定。在黄薇的全部样品中，最稳定的内参基因是GAPDH基因，这与geNorm软件得出的结果相同。
基因 不同组织 干旱胁迫 全部样品 表达稳定值 排序 表达稳定值 排序 表达稳定值 排序 GAPDH 0.069 2 0.142 1 0.359 1 UBC 0.593 5 0.452 3 0.578 2 EF-1α 0.904 6 0.800 5 0.862 3 TUA 1.483 9 0.239 2 0.910 4 CYP 0.069 1 0.611 4 0.911 5 DNA J 0.579 4 1.479 6 1.094 6 18S RNA 0.494 3 1.542 8 1.366 7 ACT 1.095 7 1.507 7 1.390 8 TUB 1.294 8 1.816 9 1.505 9 Table 2. Expression stability values of nine candidate reference genes calculated by the NormFinder
-
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α基因被鉴定为最稳定的内参基因。
基因 不同组织 干旱胁迫 全部样品 CV±SD 排序 CV±SD 排序 CV±SD 排序 EF-1α 3.46 ± 0.95 2 1.97 ± 0.53 1 2.58 ± 0.71 1 GAPDH 2.65 ± 0.69 1 3.78 ± 0.94 4 3.59 ± 0.91 2 UBC 3.73 ± 1.04 4 4.32 ± 1.18 5 4.21 ± 1.16 3 TUB 3.47 ± 1.06 3 6.14 ± 1.86 9 4.98 ± 1.51 4 TUA 6.65 ± 1.89 8 3.52 ± 0.95 3 5.26 ± 1.45 5 ACT 5.03 ± 1.46 5 6.11 ± 1.84 8 5.72 ± 1.69 6 CYP 6.36 ± 1.99 6 3.34 ± 0.96 2 6.11 ± 1.83 7 DNA J 7.31 ± 2.18 9 6.11 ± 1.81 7 6.81 ± 2.03 8 18S RNA 6.57 ± 2.05 7 5.68 ± 1.64 6 7.61 ± 2.28 9 Table 3. Expression stability values of nine candidate reference genes calculated by the BestKeeper
-
上述3种软件的结果存在一定的差异，使用在线分析软件RefFinder，整合上述分析结果并进行综合分析。由表4显示：不同组织中，最稳定的内参基因是GAPDH和CYP，最不稳定的内参基因是TUA；在干旱胁迫下，最稳定的内参基因是GAPDH和TUA，最不稳定的内参基因是TUB。在全部样本中，最稳定的内参基因为GAPDH，最不稳定的内参基因为18S RNA。
排序 不同组织 干旱胁迫 全部样品 基因 表达稳
定值基因 表达稳
定值基因 表达稳
定值1 GAPDH 1.41 GAPDH 1.19 GAPDH 1.19 2 CYP 2.11 TUA 1.86 UBC 2.11 3 18S RNA 3.57 EF-1α 3.34 EF-1α 2.45 4 UBC 4.76 UBC 3.41 TUA 4.28 5 TUB 4.76 CYP 4.00 CYP 5.23 6 EF-1α 5.05 ACT 6.70 DNA J 6.45 7 ACT 5.12 DNA J 6.96 ACT 6.73 8 DNA J 5.44 18S RNA 7.20 TUB 6.84 9 TUA 9.00 TUB 9.00 18S RNA 7.45 Table 4. Expression stability of candidate reference genes ranked by RefFinder
-
为验证候选内参基因的稳定性，采用RT-qPCR技术，根据筛选的稳定内参基因，对干旱胁迫下的靶基因CSLD和SOD的表达模式进行验证。CSLD基因是一种纤维素合成酶的结构蛋白，在茎中的表达量一般较高[24]。SOD基因属于抗氧化酶功能基因，在植物遭遇干旱胁迫时上调表达以抵御外界伤害[25]。从图4可见：在不同组织中，GAPDH和CYP基因在单独和组合校正下，CSLD均在茎中表达量最高，在其他组织中表达变化不大。但使用TUA作内参基因时，CSLD在茎中的表达量最低，在叶片中表达却最高。干旱胁迫下，使用GAPDH和基因组合GAPDH+TUA时，SOD基因在T3处理表达量最高，表达趋势大致相同；而利用TUB基因进行验证时，SOD基因在T2处理表达量最高，与其他结果差异较大。以上结果进一步验证了GAPDH基因作为内参基因的稳定性与可靠性。
Figure 4. Expression analysis for CSLD gene (A) and SOD gene (B) using different reference genes in different tissues and drought stress
-
RT-qPCR技术是一种分析基因表达量的研究工具，而使用相对定量法检测的可靠性，极大依赖内参基因的稳定性[26]。理想的内参基因在不同条件下，基因的表达变化较小，常用的内参基因一般都是维持植物生命活动必须的看家基因[27]。本研究基于黄薇转录组数据，筛选了9个常用内参基因(EF-1α、TUA、CYP、GAPDH、18S RNA、UBC、TUB、ACT、DNA J)，其中，GAPDH在不同组织和干旱胁迫下都较稳定，最不稳定的是18S RNA。geNorm和NormFinder结果中，得出的最适内参基因组合为GAPDH + UBC，BestKeeper得出最佳基因为EF-1α，但GAPDH的稳定值排名综合也较靠前。3个筛选软件得出的结果基本一致，但也存在一定的差异，这可能是由于算法的不同所导致的[28]。与本研究结果不同的是，张海洋等[29]探究了菠菜Spinacia oleracea在不同胁迫处理下稳定的内参基因，3个软件的结果一致，在氯化钠和高温胁迫下，G6PD的表达稳定性最好，聚乙二醇胁迫下ELF1B的表达稳定性最好。王蕊等[30]研究了不同发育时期大豆Glycine max在不同组织、不同非生物胁迫下稳定表达的最适内参基因，4个软件在全部组织及全部胁迫中综合分析结果均一致，在全部组织中最适内参基因为ACT，全部胁迫中最适内参基因为EF-1α。
非生物胁迫会改变植物正常的代谢过程，筛选在胁迫处理下稳定表达的内参基因，有利于后续开展目的基因的鉴定，以便进一步了解目的基因的功能[31]。目前已有大量研究筛选不同非生物胁迫下的最佳内参基因。如紫鸭跖草Tradescantia pallida在不同浓度的铜离子胁迫下，最佳内参基因组合为18S RNA + TUB + UBI[32]。马铃薯Solanum tuberosum在干旱和渗透胁迫条件下，EF-1α和sec3基因表达最稳定[33]。上述结果与本研究的结果存在一定差异，说明不同植物在不同非生物胁迫下最适内参基因也在变化，特定条件下的最适内参基因仍需通过试验筛选。
Screening and validation of reference genes in Heimia myrtifolia in different tissues and under drought stress
doi: 10.11833/j.issn.2095-0756.20220342
- Received Date: 2022-05-06
- Accepted Date: 2022-12-07
- Rev Recd Date: 2022-11-30
- Available Online: 2023-05-22
- Publish Date: 2023-05-20
Abstract:
| Citation: | ZHAO Yu, LIN Lin, WANG Qun, ZHANG Guozhe, WANG Jie, SHANG Linxue, HONG Sidan, MA Qingqing, GU Cuihua. Screening and validation of reference genes in Heimia myrtifolia in different tissues and under drought stress[J]. Journal of Zhejiang A&F University, 2023, 40(3): 665-672. doi: 10.11833/j.issn.2095-0756.20220342
|
DownLoad: