孫吉康 王平 賈浩 周韜 吳艷

摘 要： 為了解蜆殼花椒種子萌發(fā)的分子機制，需要篩選蜆殼花椒種子萌發(fā)時期表達穩(wěn)定的內(nèi)參基因。該研究通過赤霉素處理種子促進萌發(fā)，以不同萌發(fā)階段的蜆殼花椒種子為材料，采用實時熒光定量PCR技術分析了6個候選內(nèi)參基因GAPDH、ACT、18SrRNA、UBQ5、TUA和CYP在蜆殼花椒種子萌發(fā)時期的表達穩(wěn)定性。結果表明：（1）α-淀粉酶基因、DELLA基因和異檸檬酸裂解酶基因分別反應了種子萌發(fā)階段糖、激素和脂肪的代謝活動，因此選擇蜆殼花椒異檸檬酸裂解酶基因（Unigene0032088）、α-淀粉酶基因（Unigene0033597）和DELLA基因（Unigene0058868）作為驗證基因進行相對表達量驗證。（2）綜合geNorm、NormFinder和BestKeeper的分析結果顯示在蜆殼花椒種子萌發(fā)過程中ACT表達穩(wěn)定性最好，UBQ5次之。（3）以ACT、UBQ5基因為內(nèi)參基因的結果顯示驗證基因的表達量與種子萌發(fā)生理狀態(tài)一致，初步揭示了GA處理的種子易于萌發(fā)而清水處理的種子在萌發(fā)第3天容易腐敗這一現(xiàn)象出現(xiàn)的可能原因。綜上所述，ACT是蜆殼花椒種子萌發(fā)時期最合適的內(nèi)參基因，其次是UBQ5。

關鍵詞： 蜆殼花椒， 種子萌發(fā)， 實時熒光定量 PCR， 內(nèi)參基因， 基因表達

中圖分類號： Q943.2 文獻標識碼： A 文章編號： 1000-3142（2018）09-1136-10

Abstract： In order to understand the molecular mechanism of the seed germination， we need to screen the reference genes of stable expression during the seed germination of Zenthoxylum dissitum. The seeds of Z. dissitum were treated with the gibberellins （GA） to promote germination. The samples were collected in 0， 1， 2， 3 and 5 d of germination. The quantitative real-time PCR （qRT-PCR） was performed to analyze the expression stabilities of six candidates such as Glyceraldehyde-3-phosphate （GAPDH）， Actin （ACT）， 18S ribosomal RNA （18SrRNA）， ubiquitin-40S ribosomal （UBQ5）， α-tubulin （TUA） and cyclophilin （CYP） among the seeds of the different germination stages. The α-amylase gene， DELLA gene and isocitrate lyase gene respectively reflected the metabolic activities of the sugars， hormones and lipids during the germination of Z. dissitum seeds. Therefore， isocitrate lyase （Unigene0032088）， α-amylase （Unigene0033597） and DELLA （Unigene0058868） of Z. dissitum were selected as the validation genes for the relative expression verification during the seed germination. The comprehensive analysis of three softwares including geNorm， NormFinder and BestKeeper indicated that the expression stablility of ACT was the best and the second was UBQ5. The expressions of the validation genes were consistent with the states of seeds germination when ACT and UBQ5 were selected as the reference genes. Besides， when the ACT and UBQ5 genes were selected as the reference genes， the relative expression levels of isocitrate lyase， α-amylase and DELLA of zanthoxylum revealed the possible reasons that the seeds treated with GA were able to germinate successfully but the seeds treated with water easily spoiled in the 3rd day of germination. In conclusion， ACT is the preferred reference gene during the seed germination of Z. dissitum， followed by UBQ5.

Key words： Zanthoxylum dissitum， seed germination， quantitative real-time PCR （qRT-PCR）， reference gene， gene expression

蜆殼花椒（Zanthoxylum dissitum）屬于蕓香科花椒屬木質(zhì)藤本植物，主要生長于我國南方海拔500～1 200 m的山地林下、林緣以及灌木叢中，是南方地區(qū)重要的林藥資源。蜆殼花椒種子自然萌發(fā)率極低，目前野生植株采伐數(shù)量超過自然生長率，導致野生蜆殼花椒資源受到嚴重破壞，藥材原料匱乏（馬英姿和王平，2008；程鵬等，2013）。因此開展蜆殼花椒種子萌發(fā)分子機理的研究有助于促進對蜆殼花椒資源的保護。

定量實時聚合酶鏈反應（qRT-PCR）技術常用于基因表達水平分析，具有定量準確、高靈敏度、低成本以及高通量特性。內(nèi)參基因用于標準化轉(zhuǎn)錄本的豐度，是決定qRT-PCR結果準確性的主要因素（Huggett et al，2005）。理想的內(nèi)參基因不受任何實驗處理的影響并且在不同組織以及不同發(fā)育階段表達穩(wěn)定（Czechowski et al，2005）。實際上，影響基因表達的因素很多，如實驗材料、樣品量、細胞活性、RNA完整性、cDNA質(zhì)量以及逆轉(zhuǎn)錄效率等（Bustin，2002）。相關研究也證實使用表達不穩(wěn)定的內(nèi)參基因?qū)@著錯估基因表達量并導致對數(shù)據(jù)的錯誤分析（Ferguson et al，2010；Mafra et al，2012）。基因表達具有顯著的組織特異性，同時基因表達也會因為實驗組織所處生理狀態(tài)的不同或?qū)嶒炋幚矸绞降牟煌兓?。目前，還沒有發(fā)現(xiàn)有基因能作為通用的內(nèi)參基因（Nicot et al，2005；Gutierrez et al，2008；Gimeno et al，2014）。RNA測序（RNA-seq）廣泛應用于各種物種的轉(zhuǎn)錄組分析（Huggett et al，2005；Czechowski et al，2005；Ferguson et al，2010），同時也用于搜索內(nèi)參基因。Zhuang et al（2015）結合黃花棘豆轉(zhuǎn)錄組數(shù)據(jù)庫從18SrRNA、actin、tubulin、GAPDH 和 histone H3（HIS）等12個候選內(nèi)參基因中篩選出HIS 和 actin是黃花棘豆在非生物脅迫下表達最穩(wěn)定的基因。Li et al（2016）通過轉(zhuǎn)錄組數(shù)據(jù)庫在牡丹花的不同發(fā)育時期從actin、tubulin、ubiquitin、GAPDH、EF-1α 和 cyclophilin（CYP）等10個候選參考基因中篩選出GAPDH和ubiquitin在花的不同發(fā)育時期表達穩(wěn)定性最好。然而關于種子萌發(fā)過程中表達穩(wěn)定的內(nèi)參基因研究很少，目前僅發(fā)現(xiàn)Li et al（2012）以大豆發(fā)芽種子為材料，篩選出在大豆種子萌發(fā)過程最適合的內(nèi)參基因是Glyma05g37470和Glyma08g28550。迄今為止，沒有發(fā)現(xiàn)蕓香科花椒屬內(nèi)參基因研究的相關報道。

基于我們的蜆殼花椒種子萌發(fā)過程中的基因轉(zhuǎn)錄組數(shù)據(jù)并參考相關文獻，我們選擇6個常用內(nèi)參基因GAPDH、ACT、18SrRNA、UBQ5、TUA和CYP作為候選內(nèi)參基因。這些候選內(nèi)參基因通過qRT-PCR進行表達穩(wěn)定性分析，采用NormFinder，BestKeeper和geNorm（Andersen et al，2004；Vandesompele et al，2002；Pfaffl，2001）軟件對多個候選內(nèi)參基因的相對表達穩(wěn)定性并進行評價、排序，并通過種子萌發(fā)過程中具有重要功能的異檸檬酸裂解酶基因、α-淀粉酶基因和DELLA基因為驗證基因進行相對表達量的驗證分析，以確定蜆殼花椒種子萌發(fā)時期表達最穩(wěn)定的內(nèi)參基因，為更準確分析蜆殼花椒種子萌發(fā)過程中的基因表達數(shù)據(jù)提供理論依據(jù)。

1 材料與方法

1.1 材料

本實驗所用蜆殼花椒種子于2016年12月取自張家界地區(qū)采集種子，晾干后置于4 ℃保存。開始萌發(fā)前，將去殼種子浸泡在赤霉素溶液（1% w/v）3.5 h。以無菌水浸泡的種子作為對照。在9 cm有機玻璃培養(yǎng)皿中，將種子置于無菌雙蒸去離子水浸透的無菌濾紙上，然后在光照培養(yǎng)箱中（白天溫度18 ℃ 16 h，夜間溫度10 ℃ 8 h，光強103 W·m-2）萌發(fā)。清水浸泡的種子在萌發(fā)第3天通常腐壞，而GA浸泡種子可以繼續(xù)萌發(fā)并且一周后胚根突出種皮0.2～0.5 cm。在GA處理種子萌發(fā)的第0天、第1天、第2天、第3天、第5天取樣作為內(nèi)參基因篩選的樣本（第3天重復取樣），標記為0D、1D、2D、3-1D、3-2D、5D。驗證基因的樣本則在萌發(fā)開始的第1天和第2天取樣，得到5個驗證樣本，即C0-萌發(fā)0 d，W1-清水處理后萌發(fā)1 d，W2-清水處理后萌發(fā)2 d，GA1-GA處理后萌發(fā)1 d，GA2-GA處理后萌發(fā)2 d。

1.2 試劑與儀器

1.2.1 試劑 RNA提取試劑盒購自天根生化科技有限公司。反轉(zhuǎn)錄試劑盒TaKaRa PrimeScriptTM II 1st Strand cDNA Synthesis Kit 、DNA Marker及實時定量PCR試劑盒購自大連寶生物工程公司。引物由北京全式金生物技術有限公司合成。

1.2.2 儀器 智能光照培養(yǎng)箱-GZH-328A，ABI 7500實時熒光定量PCR儀。

1.3 方法

1.3.1 蜆殼花椒種子RNA的提取、純化及cDNA合成 將去除種皮的蜆殼花椒種子迅速置于液氮中研磨，使用RNA植物試劑盒（Tiangen Biotech，Beijing，China）提取總RNA?？俁NA的濃度和純度用NanoDrop 2000分光光度計檢測，總RNA的完整性用1.2%的瓊脂糖凝膠電泳檢驗。為進行qRT-PCR，使用PrimeScriptTM RT試劑盒和gDNA Eraser（TaKaRa Bio Inc.，Dalian，China）按照使用說明書配制20 μL的反應體系（總RNA 1.5 μg，1 μL dNTP mix（10 mmol·L-1），1 μL隨機引物（100 μmol·L-1））用于cDNA合成。最后用無核酸酶的水將cDNA稀釋40倍用于qRT-PCR。

1.3.2 內(nèi)參基因的選擇、引物設計及PCR反應 6個候選內(nèi)參基因分別是GAPDH，ACT，UBQ5，CYP，TUA和18SrRNA基因（表1）。引物通過網(wǎng)站IDTdna的primerquest工具設計，參數(shù)進行如下設置： GC含量45%～55%，擴增序列長度100～200個堿基對，引物長度為17～28個堿基，前導鏈與后隨鏈引物的溶解溫度值相差不超過5度。熒光定量分析采用大連寶生物工程有限公司的SYBR Premix ExTaqTM II試劑盒。使用ABI 7500實時定量PCR儀和96孔板。每個反應為20 μL體系： 4 μL模板cDNA，前導鏈與后隨鏈引物均0.8 μL（10 μmol·L-1），10 μL SYBR Premix Ex TaqTM II，4.4 μL ddH2O。PCR反應條件設定為95 ℃預變性30 s；95 ℃變性5 s，58 ℃退火20 s，72 ℃延伸15 s，40個循環(huán)反應。每個樣品3次重復。

1.3.3 內(nèi)參基因引物擴增效率（E）計算及特異性鑒定 將種子萌發(fā)0、1、2、3、5 d 5個樣品的cDNA模板按照相同分量進行混合，然后以5倍比例稀釋成4個梯度，則得到混合的cDNA模板濃度分別是1、1/5、1/25、1/125倍。每個反應重復3次。溶解溫度通過軟件7500 Software v2.3得到。在EXCEL中根據(jù)Ct值繪制標準曲線，回歸分析得到斜率（K）和線性相關系數(shù)（R2）；利用公式E=5-1/K-1，計算引物的擴增效率（E）。引物的擴增特異性鑒定則將qRT-PCR反應產(chǎn)物進行2%瓊脂糖凝膠電泳以觀察非特異性擴增條帶是否存在。

1.3.4 內(nèi)參基因篩選的數(shù)據(jù)處理與分析 geNorm、

NormFinder和BestKeeper 三個軟件被用來評價候選內(nèi)參基因在蜆殼花椒種子萌發(fā)時期的表達穩(wěn)定性。通過他們的綜合比較分析篩選出表達最穩(wěn)定的內(nèi)參基因。在geNorm和NormFinder軟件處理數(shù)據(jù)之前，需要先將Ct值轉(zhuǎn)化成Q值后進行分析。轉(zhuǎn)化公式：Q=ECtmin-Ctsample，E為引物的擴增效率，當引物擴增效率接近100%時，E值一般為2。Ctsample指的是基因的Ct值，Ctmin指的是基因在所有組織中的最小Ct值。BestKeeper軟件直接輸入Ct值進行計算分析。

2 結果與分析

2.1 蜆殼花椒候選內(nèi)參基因的選擇

以擬南芥中表達最穩(wěn)定的基因和在qRT-PCR研究中常用的內(nèi)參基因為參照，與蜆殼花椒轉(zhuǎn)錄組數(shù)據(jù)進行BLASTX后得到蜆殼花椒同源的候選內(nèi)參基因。表1描述了6個候選內(nèi)參基因的qRT-PCR引物序列和擴增序列特征。候選內(nèi)參基因按照擬南芥同源基因以及蜆殼花椒轉(zhuǎn)錄組Nr注釋命名。

2.2 擴增效率和擴增特異性

根據(jù)已知的候選內(nèi)參基因序列分析，設計溶解溫度值在58～62 ℃之間的特異性引物。以5倍濃度梯度稀釋的cDNA為模板，進行qRT-PCR擴增，根據(jù)獲得的數(shù)據(jù)作6個候選內(nèi)參基因的標準曲線。結果顯示各候選內(nèi)參基因的決定系數(shù)R2≥0.969，引物的擴增效率在91.99%～104.0%之間，均符合qRT-PCR對擴增效率的要求（表2）。分析溶解曲線發(fā)現(xiàn)各內(nèi)參基因都只產(chǎn)生單一的溶解峰，各重復樣品間的曲線重復性好（圖1）。上述結果說明所設計引物的特異性良好，qRT-PCR反應專一性高，結果準確可靠。

2.3 內(nèi)參基因Ct值分析

Ct值與基因的表達量成反比關系，即Ct值越小，表達量越大。qRT-PCR的結果表明，6個候選內(nèi)參基因的Ct值在10～30之間，其中18SrRNA的Ct值最低，在10～14之間；TUA的Ct值較高，在26～30之間；ACT、GAPDH、UBQ5和CYP的Ct值則在22～28之間（表3）。表3結果表明，18SrRNA在蜆殼花椒種子萌發(fā)早期的不同發(fā)育時期表達豐度最高，其次是CYP、UBQ5、ACT，而TUA的表達豐度最低。同時，GAPDH、TUA在不同發(fā)育時期表達豐度變異較大?？梢?，這些內(nèi)參基因在不同發(fā)育時期表達豐度的變化規(guī)律較為復雜，有必要采用內(nèi)參基因分析軟件對內(nèi)參基因表達的穩(wěn)定性進行評估。

2.4 候選內(nèi)參基因穩(wěn)定性分析

BestKeeper、geNorm和NormFinder軟件常被用于對候選內(nèi)參基因進行評價。其中BestKeeper程序直接輸入Ct值進行表達的穩(wěn)定性分析； 而geNorm和NormFinder軟件在處理數(shù)據(jù)前需要先將Ct值轉(zhuǎn)換成相對表達量再進行穩(wěn)定性分析。我們用這3個軟件對GAPDH、ACT、UBQ5、CYP、TUA和18SrRNA候選內(nèi)參基因在蜆殼花椒種子萌發(fā)早期的表達穩(wěn)定性進行統(tǒng)計學分析。

2.4.1 geNorm分析 geNorm V3.5軟件根據(jù)平均表達穩(wěn)定指數(shù)M值確定在不同條件下表達最穩(wěn)定的內(nèi)參基因。M值越大，基因的表達穩(wěn)定性越低，反之越高。M值小于1.5的候選內(nèi)參基因才被認為表達相對穩(wěn)定。本研究分析結果表明，6個候選內(nèi)參基因中3個的M值均小于1.5，分別是ACT、UBQ5、18SrRNA，說明這3個候選內(nèi)參基因都較為穩(wěn)定。在3個候選內(nèi)參基因中，ACT的M值最低，說明它的穩(wěn)定性最好。6個候選內(nèi)參基因的穩(wěn)定性從高到低依次為ACT >UBQ5>18SrRNA>CYP>GAPDH>TUA（圖2）。另外，M值都偏大，估計與種子萌發(fā)時期表達豐度的變異程度較大有關。

2.4.2 BestKeeper分析 BestKeeper分析的原理是直接將候選內(nèi)參基因的Ct值用EXCEL軟件計算標準差（s）。一般s值越小，基因的表達穩(wěn)定性越好。當s值>1時表明該基因表達不穩(wěn)定。本研究分析結果（表4）顯示，GAPDH和TUA的s值大于1，其余都小于1；其中ACT的s最小，表明ACT的表達穩(wěn)定性最好。按照穩(wěn)定性從高到低排序依次為ACT>UBQ5> CYP>18SrRNA>TUA > GAPDH。

2.4.3 NormFinder分析 NormFinder程序依據(jù)相對表達量用EXCEL軟件計算基因的穩(wěn)定值M。M值越低，則該基因越穩(wěn)定。6個候選內(nèi)參基因的M值如表5所示。表5結果表明，UBQ5的M值最低，說明它是所選內(nèi)參基因中穩(wěn)定性最好的；TUA的M值最高，說明它的穩(wěn)定性最差。6個候選內(nèi)參基因按照穩(wěn)定性從高到低排序依次為UBQ5>ACT>18SrRNA>CYP>GAPDH>TUA。

2.5 內(nèi)參基因的綜合穩(wěn)定性分析

使用BestKeeper、geNorm和NormFinder 3個軟件對蜆殼花椒中6個候選內(nèi)參基因的穩(wěn)定性進行統(tǒng)計學綜合分析（表6）。geNorm軟件分析結果顯示，ACT和UBQ5的M值最低，說明它們的穩(wěn)定性較好，6個候選內(nèi)參基因按照穩(wěn)定性從高到低排序依次為ACT>UBQ5>18SrRNA>CYP>GAPDH>TUA。BestKeeper軟件分析結果顯示，ACT的s值最低，穩(wěn)定性較好；GAPDH和TUA的s值大于1，穩(wěn)定性最差，6個候選內(nèi)參基因的穩(wěn)定性排序為ACT>UBQ5>CYP>18SrRNA>TUA>GAPDH。NormFinder軟件分析結果顯示，UBQ5的M值最低，說明其穩(wěn)定性較好，穩(wěn)定性排序為UBQ5>ACT>18SrRNA>CYP>GAPDH>TUA。3個軟件的分析結果一致性較好，結果可靠度高。從綜合分析結果來看，ACT、UBQ5和18SrRNA適合作為內(nèi)參基因，但ACT和UBQ5表達更穩(wěn)定且表達豐度也適中。因此，優(yōu)先選擇ACT和UBQ5作為內(nèi)參基因。

2.6 內(nèi)參基因穩(wěn)定性驗證

以ACT、UBQ5、18SrRNA和GAPDH作為內(nèi)參基因，分析異檸檬酸裂解酶基因Unigene0032088、α-淀粉酶基因Unigene0033597、DELLA基因Unigene0058868在種子萌發(fā)初始階段的相對表達量，以研究不同內(nèi)參基因?qū)嶒灲Y果的影響。3個驗證基因分別是與種子萌發(fā)時期糖代謝、脂代謝、激素密切相關的重要基因，能夠較好地反應種子萌發(fā)時期的生理狀態(tài)。

圖3結果表明，以ACT、UBQ5、18SrRNA為內(nèi)參基因時3個驗證基因在各樣本中的表達趨勢基本一致，而且在GA處理種子萌發(fā)的第1天、第2天都活躍表達，這與種子萌發(fā)的生理過程相一致。以GAPDH為內(nèi)參基因時驗證基因的表達趨勢則明顯不同。進一步分析發(fā)現(xiàn)以GAPDH為內(nèi)參基因時驗證基因表達量在GA種子萌發(fā)的第1天、第2天普遍偏小，反而在第0天高表達，這顯然與GA種子萌發(fā)趨勢不一致，而且表3也顯示GAPDH的表達趨勢是在萌發(fā)的第0天、第1天、第2天表達量明顯增大，因此當以GAPDH為內(nèi)參基因時驗證基因在第1天、第2天的相對表達量會明顯偏小。GAPDH（甘油醛-3-磷酸脫氫酶）是參與糖酵解的一種關鍵酶。在種子萌發(fā)早期， 外源GA促進α-淀粉酶表達以加強對糖的分解利用，可見GAPDH表達量并不穩(wěn)定，因此GAPDH不適合為內(nèi)參基因。綜上所述，在用qRT-PCR分析基因表達時，選擇不恰當?shù)膬?nèi)參基因?qū)⒖赡軐е洛e估目標基因的相對表達量，由此可見內(nèi)參基因選擇的重要性。

3 討論與結論

種子萌發(fā)是一個復雜的植物生理發(fā)育過程，其中許多植物激素在種子萌發(fā)過程中具有重要作用：如赤霉素GA、油菜素內(nèi)脂BR促進種子萌發(fā)，而脫落酸ABA、茉莉酸JA抑制萌發(fā)，但可以增強植物抗氧化和抗微生物感染能力。GA是促進種子萌發(fā)的激素，在谷物種子萌發(fā)早期GA誘導糊粉層中的α-淀粉酶基因轉(zhuǎn)錄表達（Ye et al，2011）。DELLA蛋白抑制GA信號傳導通路（徐恒恒等，2014），并負調(diào)節(jié)BR信號途徑（Iliev et al，2002），同時促進JA信號通路應答基因表達（Aleman et al，2016）。由此可見，DELLA 蛋白是多種激素間相互作用的關鍵調(diào)控因子。蜆殼花椒種子富含油脂，油脂在種子萌發(fā)過程中可以提供碳源與能量。種子油脂一般先分解為脂肪酸，并經(jīng)過β-氧化途徑轉(zhuǎn)變?yōu)橐阴］o酶A，然后通過乙醛酸循環(huán)生成糖加以利用。異檸檬酸裂解酶是乙醛酸循環(huán)途徑關鍵酶，通過對異檸檬酸裂解酶的表達情況可以了解種子萌發(fā)過程中油脂動員情況。因此，α-淀粉酶基因、DELLA基因和異檸檬酸裂解酶基因可以較好地反映種子萌發(fā)階段糖、激素和脂肪的相關代謝活動。

本研究結果表明，蜆殼花椒種子萌發(fā)時期的合適內(nèi)參基因是ACT、UBQ5和18SrRNA。當以ACT、UBQ5和18SrRNA為內(nèi)參基因時，3個驗證基因具有相似的相對表達量，且在清水處理種子萌發(fā)第2天表達量都極低，表明此時種子糖、脂、激素相關代謝活動幾乎處于停滯狀態(tài)，顯示清水處理種子的萌發(fā)過程在萌發(fā)第2天被中斷。這一現(xiàn)象可能與清水處理種子在第3天容易腐壞有關。考慮到18SrRNA表達豐度遠高于大多數(shù)基因，因此ACT是蜆殼花椒種子萌發(fā)時期最合適的內(nèi)參基因，其次是UBQ5。

本研究中，當以ACT、UBQ5為內(nèi)參基因時，異檸檬酸裂解酶基因的表達量在GA處理的種子萌發(fā)第2天具有極大提高，表明在GA處理的種子萌發(fā)第2天乙醛酸循環(huán)旺盛，此時種子油脂的動員也非?；钴S。α-淀粉酶基因在GA處理種子萌發(fā)第1天、第2天均有高表達，同時清水種子在萌發(fā)第1天也有一定程度的表達，表明無論是GA處理種子或是清水種子，在萌發(fā)早期糖代謝都較為活躍。DELLA基因在GA處理種子萌發(fā)的第1天高表達，但是在第2天顯著下調(diào)，表明DELLA對于GA信號通路的抑制作用減弱，有利于誘導谷物種子的糊粉層中α-淀粉酶的轉(zhuǎn)錄和種子萌發(fā)?？梢姡訟CT、UBQ5為內(nèi)參基因的相對表達量分析初步揭示了GA處理的種子易于萌發(fā)而清水處理的種子在萌發(fā)第3天容易腐敗這一現(xiàn)象出現(xiàn)的可能原因。

總之，本研究以不同萌發(fā)階段的蜆殼花椒種子為材料，分析了6個候選內(nèi)參基因GAPDH、ACT、18SrRNA、UBQ5、TUA和CYP在蜆殼花椒種子萌發(fā)時期的表達穩(wěn)定性。綜合geNorm、NormFinder和BestKeeper的分析結果并結合基因的相對表達量分析驗證，結果表明蜆殼花椒種子萌發(fā)時期最合適的內(nèi)參基因是ACT，其次是UBQ5。這個結果為蜆殼花椒種子萌發(fā)過程中的熒光定量PCR分析提供了詳盡的內(nèi)參基因參考。

參考文獻：

ALEMAN F， YAZAKI J， LEE M， et al， 2016. An ABA-increased interaction of the PYL6 ABA receptor with MYC2 transcription factor： A putative link of ABA and JA signaling [J]. Sci Rep， 30（6）：28941-28958.

ANDERSEN CL， JENSEN JL， ORNTOFT TF， 2004. 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 Res， 64（15）：5245-50.

BUSTIN S， 2002. Quantification of mRNA using real-time reverse transcription PCR （RT-PCR）： trends and problems [J]. J Mol Endocrinol， 29（1）：23-39.

CHENG P， WANG P， SUN JK， et al， 2013. Research progresses on regulation mechanisms of plant seeds dormancy and germination [J]. J Centr S Univ For Technol， 33（5）：52-58. [程鵬， 王平， 孫吉康， 等， 2013. 植物種子休眠與萌發(fā)調(diào)控機制研究進展 [J]. 中南林業(yè)科技大學學報， 33（5）：52-58.]

CZECHOWSKI T， STITT M， ALTMANN T， et al， 2005. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis [J]. Plant Physiol， 139（1）：5-17.

FERGUSON BS， NAM H， HOPKINS RG， et al， 2010. Impact of reference gene selection for target gene normalization on experimental outcome using real-time qRT-PCR in Adipocytes [J]. PLoS ONE， 5（12）：e15208-19.

GIMENO J， EATTOCK N， DEYNZE AV， et al， 2014. Selection and validation of reference genes for gene expression analysis in switchgrass （Panicum virgatum） using quantitative real-time RT-PCR [J]. PLoS ONE， 9（3）：e91474-80.

GUTIERREZ L， MAURIAT M， GUENIN S， et al， 2008. 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 Biotechnol J， 6（6）：609-618.

HUGGETT J， DHEDA K， BUSTIN S， et al， 2005. Real-time RT-PCR normalisation， strategies and considerations [J]. Genes Immun， 6（4）：279-284.

ILIEV EA， XU W， POLISENSKY DH， et al， 2002. Transcriptional and posttranscriptional regulation of Arabidopsis TCH4 expression by diverse stimuli. Roles of cis regions and brassinosteroids [J]. Plant Physiol， 130（2）：770-83.

LI J， HAN J， HU Y， et al， 2016. Selection of reference genes for quantitative real-time pcr during flower development in tree peony （Paeonia suffruticosa Andr.） [J]. Front Plant Sci， 21（7）：516-527.

LI Q， FAN CM， ZHANG XM， et al， 2012. Validation of reference genes for real-time quantitative PCR normalization in soybean developmental and germinating seeds [J]. Plant Cell Rep， 31（10）：1789-1798.

MA YZ， WANG P， 2008. Evaluation of threatened degree and sustainable utilization countermeasures of Zanthoxylum Dissitum resource [J]. Nonw For Res， 26 （1）： 91-94. [馬英姿， 王平， 2008. 蜆殼花椒受威脅程度評價及其可持續(xù)利用對策 [J]. 經(jīng)濟林研究， 26（1）： 91-94.]

MAFRA V， KUBO KS， ALCWS-FERREIRA M， et al， 2012. Reference genes for accurate transcript normalization in citrus genotypes under different experimental conditions [J]. PLoS ONE， 7（2）：e31263-70.

NICOT N， HAUSMAN JF， HOFFMANN L， et al， 2005. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress [J]. J Exp Bot， 56 （421）：2907-2914.

PFAFFL MW， 2001. A new mathematical model for relative quantification inrealtime RT-PCR [J]. Nucl Acids Res， 29（9）：e45-53.

VANDESOMPLETE J， PRETER KD， PATTYN F， et al， 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes [J]. Genome Biol， 3（7）：1-12.

XU HH， LI N， LIU SJ， et al， 2014. Research progress in seed germination and its control [J]. Acta Agron Sin， 40 （7） ： 1141-1156. [徐恒恒， 黎妮， 劉樹君， 等， 2014. 種子萌發(fā)及其調(diào)控的研究進展 [J]. 作物學報， 40 （7） ： 1141-1156.]

YE N， ZHU G， LIU Y， et al， 2011. ABA controls H2O2 accumulation through the induction of OsCATB in rice leaves under water stress [J]. Plant Cell Physiol， 52（4）：689-98.

ZHUANG H， FU Y， HE W， et al， 2015. Selection of appropriate reference genes for quantitative real-time PCR in Oxytropis ochrocephala Bunge using transcriptome datasets under abiotic stress treatments [J]. Front Plant Sci， 30（6）：475-487.