-
文冠果（Xanthoceras sorbifolium Bunge），又名黄角（Yellow horn），落叶灌木或小乔木，是无患子科文冠果属的唯一植物，广泛种植于我国三北地区，抗性强，尤其抗旱，是我国特有的观赏、药用及油用的高收入木本经济树种[1]。文冠果种子富含油脂，种仁含油率达60%以上，不饱和脂肪酸组成高达92%，富含油酸（约30%）、亚油酸（约45%）以及在植物中极为罕见的神经酸（约3%）[2]，是生产高品质生物柴油和健康食用油的原料，具有优秀的经济价值。由此可见，研究文冠果种子不饱和脂肪酸的合成机制，对培育优质高产的文冠果品种具有现实意义。截止至今，已有文冠果种子发育转录组的文献分析出一些重要的不饱和脂肪酸合成相关候选关键基因[3-4]，探究这些基因在文冠果各组织，尤其不同发育时期种子中的表达模式，尤为关键。
实时荧光定量PCR（qRT-PCR）是大多数科学家常用的检测基因表达丰度的可靠技术，具有高灵敏性、特异性、精确性[5]，而一个可靠的内参基因（RGs）是前提。内参基因一般选择在植物中稳定表达的管家基因，常见的管家基因有actin（ACT）、糖醛酸-3-磷酸脱氢酶基因（GADPH）、延长因子1α（EF-1α）、泛素酶基因（UBQ）、泛素结合酶基因（UCE）以及18S 核糖体RNA（18S），然而这些基因的表达在不同物种的不同组织和其不同发育时期中的稳定性有较大差异，需要深入研究[6-7]。
大量关于文冠果基因表达的研究均采用actin作为内参基因[4,8-13]，王旭等通过对7个内参基因稳定性研究，发现UBQ和eIF-4α为适用于性别分别相关基因的最优内参基因[14]。然而，近年来神经酸的营养价值和药用价值被逐渐挖掘，其种子油脂合成与代谢相关基因的研究亟待深入，需要进一步研究文冠果种子不同发育时期的最适内参基因需要被进一步研究，以便提高其油脂合成与代谢相关基因表达检测的准确性。基于此，本研究首先以SapBase数据库中文冠果各组织转录组数据（http://www.sapindaceae.com/）中常见内参基因的基因表达量数据作为依据[15]，挑选12个在各组织表达量较稳定的候选内参基因；采用qRT-PCR法测试这12个候选内参基因在文冠果种子不同发育时期和不同组织中的表达水平Cq值；利用4种常见的评估内参基因的表达稳定性的数学算法评价这12个候选内参基因在文冠果不同组织和种子发育不同时期中的稳定性，随即利用在线工具RefFinder对这12个候选内参基因进行综合排名选出最佳内参基因。最后，以超长链脂肪酸合成关键酶3-酮酯酰-CoA合酶（KCS）基因验证选择的最佳内参基因的稳定性，本研究旨在选择稳定、可靠的文冠果内参基因。
-
文冠果（Xanthoceras sorbifolium Bunge）成熟实生苗（五年生以上）生长于内蒙古通辽市开鲁县（43°32′ N, 120°30′ E），林间正常管理。文冠果根、茎、叶、花采集于5月10日，种皮采集于6月7日，不同发育时期种子胚乳分别采集于6月18日（谢花后第30 d，30 DAP），6月26日（38 DAP），7月3日（45 DAP），7月11日（53 DAP），7月21日（63 DAP），采集年份均为2022年，5个不同发育时期种子胚乳形态特征见图1。每种文冠果组织采集生物学重复样品3份。
Figure 1. Seed endosperm morphology at different developmental stages (30 DAP～63 DAP)
-
文冠果各植物组织的提取方法参照多糖多酚植物总RNA提取试剂盒（天根，中国）说明书，经过降解酶RQ1 RNase-Free DNase（Promega，美国）降解残留的基因组DNA后，以1 μg RNA为模板，利用GoScriptTM Reverse Transcription System（Promega，美国）试剂盒反转录得到cDNA，稀释10倍后冻于−20 ℃待用
-
参照SapBase数据库（
http://www.sapindaceae.com/ ）中文冠果各组织转录组数据（SRA编号：SRP304628），选择在各组织表达量较稳定的12个内参基因作为候选，同时以超长链脂肪酸合成关键酶3-酮酯酰-CoA合酶（KCS）基因作为验证基因，基因编号均来源于SapBase数据库[15]。荧光定量PCR（qRT-PCR）引物由软件Primer Premier 5.0设计，引物信息列于表1。扩增序列均测序检验其正确性。基因/基因编号
Gene/Gene number基因名称
Gene name引物序列（5′–3′）正向/反向
Primer sequences (5′-3′)
forward/reverse扩增长度
Ampliction
length/bp扩增效率
Amplification
efficiency/%线性相关系数R2
Correlation
coefficient
R2平均值
Cq Average
Cq valueACT/EVM0004930 肌动蛋白 TACCGAGGCACCATTAAATCCC/
AAGGTCCAAACGAAGAATAGCA226 99.23 0.997 21.51 UBQ/EVM0019311 泛素蛋白 TTCCGTCTTGCTCTCTACACTTCC/
GTTTCTGCTTATTGACATCTTTATCC176 98.75 0.995 24.59 UCE/EVM0006084 泛素结合酶 AATGGGGCTTCAGGAACTCTCTAT/
TATCTTCAGGGCGTTGCTTCACAG256 101.14 0.993 24.33 18r/EVM0012402 18S核糖体RNA CAGGCACTACAGTCCTCACAAAA/
GTCACTTCCACCATCCATATCCA117 99.46 0.998 26.68 GADPH/EVM0014701 糖醛酸-3-磷酸脱氢酶 GCAACCTATGAGGACATCAAAGC/
CCAACGAAATCGGTAGAAACAAC107 98.23 0.992 21.59 EF-1α/EVM0013711 延伸因子 GTGCTATGTCACCCCGATTTCCC/
TTTTACCATTGTCGCAGCCTCCT142 98.88 0.995 24.38 α-TUB/EVM0000325 α微管蛋白 CAATCTCAACCGCCTTGTCTCTC/
AGCTGCTCATGGTATGCCTTCTC183 100.72 0.996 22.94 β-TUB/EVM0013594 β微管蛋白 GGGGCAAAATGAGCACTAAAGAA/
CCAGTCGGAGGGATATCACAGAC127 101.35 0.994 23.67 PP2A/EVM0018196 蛋白质磷酸酯酶 TCTTTGACTTCACTGATCCGACA/
TGGACACCATCTTTACCATTTCTT144 99.46 0.998 27.39 Tip41/EVM0004487 Tip41蛋白 AAGATGAATTGGCTGATAACGGA/
AGCTCTCACGAAGAATGACTGGG183 98.58 0.994 25.53 SAM/EVM0014230 甲硫氨酸腺苷转移酶 TGAGGTGAGAAAGAACAAGACGT/
TGCTGGGTTGAAATGAGAACAGTG132 99.73 0.998 25.34 CYP/EVM0006460 亲环蛋白 GAAAGGTGGAGAACAATGATAGAT/
TTGAAAGCACAACACAGGAGAGT152 101.12 0.993 27.15 KCS/EVM0000575 3-酮酯酰-CoA合酶 TCCAAAGGGAGGATGACAACA/
CTCGGACATGGGAAGGACTAG125 99.32 0.997 25.96 Table 1. 12 candidate reference gene primers and amplification efficiency
-
qRT-PCR由ABI StepOne7000（美国）实时荧光定量PCR仪完成，反应体系为30.0 μL，包括SYBR qPCR SuperMix 15.0 μL，上游引物（10.0 μmol·L−1）1.0 μL，下游引物（10.0 μmol·L−1）1.0 μL，ROX 0.5 μL，ddH2O 12.5 μL。反应采用三步法，反应程序为：95 °C. 2 min，1个循环；95 °C 10 s，60 °C，20 s，72 °C，20 s，40个循环；72 °C，20 s，1个循环。熔解曲线在65～95 °C间获得以验证引物的特异性。引物的扩增效率和线性相关系数由cDNA按5个稀释倍数（1、10−1、10−2、10−3、10−4）建立标准曲线计算得到。
4种常见的评估内参基因的表达稳定性的数学算法ΔCT[16]、Best-Keeper[17]、geNorm[18]及NormFinder[19]用于评价12个候选内参基因在文冠果种子不同发育时期和不同组织的稳定性。最后利用在线工具RefFinder（
http://150.216.56.64/referencegene ）综合排名。 -
所有样品RNA均利用Agilent 2100 Bioanalyzer生物分析仪测试RNA浓度和完整性，样本RNA浓度分布于128.5～433.9 ng·μL−1，RIN值均大于7.0，琼脂糖电泳检测均显示有明显的28S和18S条带，RNA质量达到后续实验要求。12个候选内参基因的基因名称、SapBase编号、引物序列、扩增长度、扩增效率、相关系数及Cq平均值均列于表1。扩增片段长度在107 bp（GADPH）～256 bp（UCE），扩增效率在98.23%（GADPH）～101.35（β-TUB）。12个候选内参基因熔解曲线均为单一峰（图2），表明无引物二聚体和非特异性扩增发生。
Figure 2. Melting curves with a single peak for the 12 candidate reference genes
-
qRT-PCR法得到12个候选内参基因在文冠果不同组织和不同种子发育时期中的表达水平Cq值，Cq值越小，基因表达量越高。这些内参基因Cq值在文冠果各个组织有较大差异。这所有检测的组织中，PP2A平均Cq值最高27.39，分布于25.47～28.69；ACT平均Cq值最低21.51，分布18.69～23.65（表1，图3）。这12个候选内参基因的变异系数分别为7.18%（ACT）、3.01%（UBQ）、4.54%（UCE）、3.24%（18r）、6.64%（GADPH）、4.56%（EF-1α）、8.84%（β-TUB）、10.68%（α-TUB）、4.13%（PP2A）、3.87%（Tip41）、9.60%（SAM）、3.11%（CYP）。
Figure 3. The quantification cycle (Cq) values of the 12 candidate reference genes
-
常见的内参基因稳定性算法有ΔCT、BestKeeper、NormFinder以及GeNorm。这4种算法中，ΔCT法是根据Cq值的标准差数值来评价，标准差数值越低的内参基因，稳定性越强；BestKeeper算法是根据Cq值的变异系数和标准差数值联合评价，变异系数和标准差数值越低，内参基因的稳定性越强；NormFinder和GeNorm算法均根据平均表达稳定指数计算，表达稳定指数与内参基因的稳定性呈负相关。
上述4种内参基因稳定性算法用于评价12个候选内参基因在文冠果种子不同发育时期和不同组织中的稳定性。如表2所示，根据ΔCT和BestKeeper算法分析，PP2A和Tip41在文冠果种子不同发育时期和不同组织中稳定性均为最优，NormFinder和GeNorm算法分析的最优结果与ΔCT和BestKeeper算法不同，但PP2A和Tip41均在排名前四的行列。在文冠果种子不同发育时期中，所有算法均显示ACT和GADPH为最不稳定的内参基因；在不同组织中，所有算法均显示ACT和β-TUB为最不稳定的内参基因。
方法
Method稳定性排名（高-低） Ranking order（high-low） 1 2 3 4 5 6 7 8 9 10 11 12 不同发育时期种子 Seeds at different stages ΔCT Tip41 PP2A UCE EF-1α UBQ 18r β-TUB CYP α-TUB SAM ACT GADPH BestKeeper Tip41 PP2A UBQ EF-1α 18r UCE β-TUB CYP SAM α-TUB ACT GADPH NormFinder β-TUB UCE Tip41 PP2A α-TUB EF-1α UBQ 18r CYP SAM ACT GADPH GeNorm EF-1α/ PP2A Tip41 18r UBQ UCE β-TUB CYP SAM α-TUB ACT GADPH RefFinder PP2A Tip41 EF-1α 18r UCE β-TUB UBQ α-TUB CYP SAM ACT GADPH 不同组织 All tissues ΔCt Tip41 PP2A 18r CYP UBQ EF-1α SAM UCE GADPH α-TUB ACT β-TUB BestKeeper Tip41 PP2A 18r Tip41 UBQ EF-1α UCE SAM GADPH α-TUB ACT β-TUB NormFinder Tip41 SAM PP2A 18r CYP UBQ EF-1α GADPH UCE α-TUB ACT β-TUB GeNorm PP2A/CYP Tip41 18r UBQ EF-1α UCE SAM GADPH α-TUB ACT β-TUB RefFinder PP2A Tip41 CYP 18r UBQ SAM EF-1α UCE GADPH α-TUB ACT β-TUB Table 2. Stability ranking of candidate reference genes in seeds at different developmental stages and across diverse tissue samples
在线工具RefFinder根据4种数学算法的排名综合评价这12个候选内参的稳定性。结果如表2所示，在文冠果不同发育时期种子内参基因综合评分排名为PP2A＞Tip41> EF-1α＞18r＞UCE＞β-TUB＞UBQ＞α-TUB＞CYP＞SAM＞ACT＞GADPH；在文冠果不同组织内参基因综合评分排名为PP2A＞Tip41>CYP＞18r＞UBQ＞SAM＞EF-1α＞UCE＞GADPH＞α-TUB＞ACT＞β-TUB。PP2A和Tip41在文冠果不同发育时期种子和不同组织中的排名均为最高，表明为最合适的内参基因。
-
为了评估上述选择的最优内参基因的稳定性，本研究以选择的最优内参基因PP2A和Tip41及文献中最常用的ACT作为后续验证实验的内参基因，利用qRT-PCR法计算种子高表达目标基因KCS在文冠果种子不同发育时期和不同组织中的表达水平。结果如图4所示，以PP2A和Tip41作为内参基因达到大致相同的KCS表达模式（图4A、B）。KCS主要在种子发育期38～53 DAP高表达。在除了种子外的其他组织中，在种皮中几乎不表达，不参与计算，在花中高表达，然而，当采用ACT作为内参基因时，其表达模式有所不同（图4C）。在以PP2A和Tip41作为内参基因时，KCS基因在种子发育30 DAP的表达最低，作为标准值“1”，在茎的表达量高于叶片；而以ACT作为内参基因在根中的表达量最低，且在叶片中的表达量高于茎。在不同发育时期种子中，KCS基因以PP2A和Tip41作为内参基因计算时，种子发育63 DAP的表达量高于30 DAP，而以ACT作为内参基因时相反。上述结果表明，以不同稳定性的内参基因计算得到的组织特异性表达分析结果有所不同。由此可见，选择可靠稳定的内参基因对qRT-PCR实验的准确度非常重要。
Figure 4. Relative expression of the KCS gene in diverse tissue samples using selected reference genes
-
实时荧光定量（qRT-PCR）技术因其具有高灵敏性、特异性、精确性等特点，是分析基因表达水平和核酸含量的应用最为广泛的技术之一。不同物种之间没有固定的内参基因，例如常用的内参基因ACT、GADPH、UBQ等在不同植物物种的不同组织的不同发育时期可表现出表达差异[20]。由此可见，为特定植物的特定组织发育的qRT-PCR实验选择可靠的内参基因极为必要。
文冠果油是在植物种油中极为罕见的含有神经酸的木本油料来源，大量文献均以内参基因ACT作为文冠果基因表达分析的唯一内参基因[4,8-13]，其稳定性是否适宜文冠果种子发育中相关基因的表达分析未有确定。本研究选择了12个常见且表达较为稳定的候选内参基因，利用4种数学算法（ΔCT、BestKeeper、NormFinder、geNorm）评估它们在文冠果种子发育过程中和不同组织中的稳定性，同时采用网络在线工具RefFinder综合这4种算法的结果并排名。
PP2A和Tip41在文冠果不同发育时期种子和不同组织中的排名均为最高，为本研究的最适内参基因，而UBQ和Tip41是油用牡丹（Paeonia ostii T. Hong et J. X. Zhang）种子发育过程中的最适内参基因，PP2A排名中游[6]；UBCE2是草果（Amomum tsaoko Crevost et Lemarie）种子发育的最适内参基因[21]；UCE、RPS13和RPII是美藤果（Plukenetia volubilis L.）种子发育最适内参基因[22]。以上结果表明不同物种不同发育组织最适内参基因大部分均不相同，值得研究。最为关键的是本研究结论中的最适内参基因与前期大部分研究报道所使用的并不一致，前期大量研究均采用ACT基因作为文冠果基因表达分析的内参基因[4,8-13]，而本研究结果表明ACT基因在文冠果种子不同发育时期和不同组织中内参稳定性排名较差，并不适合用于基因表达的研究，同时王旭等研究也表明ACT基因不是研究文冠果性别分化和不同组织的最适内参基因，且排名同样靠后[14]。据此，本研究最重要的意义是更正用于文冠果基因表达分析的内参基因，以期获得更准确的基因表达数据。
-
本研究通过4种数学算法（ΔCT、BestKeeper、NormFinder、geNorm）对文冠果12个候选内参基因稳定性和有效性的评估，发现PP2A和Tip41基因在文冠果种子不同发育时期和不同组织中稳定表达，是文冠果基因表达分析最稳定的内参基因。本研究为文冠果qRT-PCR数据标准化提供了有效依据。
Selection of Reference Genes for Gene Expression Analysis in Seed Development of Xanthoceras sorbifolium Bunge
- Received Date: 2023-07-24
- Accepted Date: 2024-02-21
- Available Online: 2024-04-01
Abstract: :
DownLoad: