羊草不同組織實時定量PCR內(nèi)參基因的篩選

胡寧寧，郭慧琴，李西良，孔令琪，武自念，張繼澤，常 春，萬東莉，臧 輝，任衛(wèi)波(.中國農(nóng)業(yè)科學院草原研究所,內(nèi)蒙古 呼和浩特 0000； 2.中國農(nóng)業(yè)科學院研究生院,北京 0008；.內(nèi)蒙古農(nóng)業(yè)大學生命科學學院,內(nèi)蒙古 呼和浩特 0008)

?

植物生產(chǎn)層

羊草不同組織實時定量PCR內(nèi)參基因的篩選

胡寧寧1,2，郭慧琴3，李西良1，孔令琪1，武自念1，張繼澤1，常 春1，萬東莉1，臧 輝1,2，任衛(wèi)波1
(1.中國農(nóng)業(yè)科學院草原研究所,內(nèi)蒙古 呼和浩特 010010； 2.中國農(nóng)業(yè)科學院研究生院,北京 100081；3.內(nèi)蒙古農(nóng)業(yè)大學生命科學學院,內(nèi)蒙古 呼和浩特 010018)

為篩選羊草(Leymuschinensis)不同組織實時定量PCR(qRT-PCR)試驗體系中的最佳內(nèi)參基因，以羊草葉、莖、根、穗為研究材料，利用qRT-PCR技術分析TUA、TUB、18SrRNA、EF-1α、APRT、CYP、Actin和CBP20共8個常用候選管家基因的表達情況，并使用geNorm和NormFinder軟件進行綜合評價，以期篩選出羊草不同組織中最穩(wěn)定的內(nèi)參基因。結(jié)果表明，8個候選內(nèi)參基因在羊草不同組織表達穩(wěn)定性不同，其中，葉片中穩(wěn)定性較高的內(nèi)參基因是Actin；莖中穩(wěn)定性較高的內(nèi)參基因是EF-1a；根中穩(wěn)定性較高的內(nèi)參基因是APRT和18SrRNA；穗中穩(wěn)定性較高的內(nèi)參基因是TUB。本研究結(jié)果將為開展羊草功能基因的表達分析提供重要參考。

羊草；qRT-PCR；內(nèi)參基因；組織

實時定量PCR(real-time quantitative PCR)是一種PCR技術，常用于基因表達檢測，具有特異性強、靈敏度高、重復性好以及高通量等優(yōu)點[1]，已成為當今植物分子生物學領域檢測基因表達與轉(zhuǎn)錄組分析等研究的重要方法之一[2]。qRT-PCR包括絕對定量和相對定量兩種方式，常用到的是相對定量。在相對定量檢測中，通常會選擇內(nèi)參基因來校正數(shù)據(jù)，用來消除試驗過程中如RNA提取質(zhì)量、反轉(zhuǎn)錄cDNA效率及PCR反應條件可能帶來的差異，使目的基因在試驗樣品中表達數(shù)據(jù)更準確。

在qRT-PCR中，使用穩(wěn)定的內(nèi)參基因是精準試驗數(shù)據(jù)的前提[3]。理想的內(nèi)參基因，應該是在所有的生理狀態(tài)下都能穩(wěn)定地表達。然而，有大量研究證實，在同一物種的不同生長時期、不同組織部位以及不同培養(yǎng)條件下，始終穩(wěn)定表達的內(nèi)參基因是不存在的[4]。因此，實時熒光定量PCR中內(nèi)參基因的篩選和穩(wěn)定性評價至關重要。目前，已有不少關于植物內(nèi)參基因篩選評價的研究報道。如小麥(Triticumaestivum)[5]、水稻(Oryzasativa)[6]、玉米(Zeamays)[7]、大豆(Glycinemax)[8]、茶樹(Camelliasinensis)[9-10]、刺槐(Robiniapseudoacacia)[11]、蘋果(Malusdomestica)[12]、松葉豬毛菜(Salsolalaricifolia)[13]、大針茅(Stipagrandis)[14-15]等。

羊草(Leymuschinensis)又名堿草，為禾本科多年生C3草本植物，是生長在歐亞大陸草原的一種優(yōu)質(zhì)牧草，具有抗寒、抗旱、耐鹽堿等特性。目前，有關羊草的研究報道有很多，但主要集中在放牧利用[16]、功能性狀[17]以及生理生化處理[18]等方面，還沒有任何關于羊草實時熒光定量PCR內(nèi)參基因篩選的報道。本研究選取8個候選內(nèi)參基因，包括α-微管蛋白基因TUA、β-微管蛋白基因TUB、核糖體18SrRNA、轉(zhuǎn)錄延伸因子EF-1α、腺嘌呤磷酸核糖轉(zhuǎn)化酶基因APRT、親環(huán)蛋白基因CYP、肌動蛋白基因Actin、類胡蘿卜素結(jié)合蛋白基因CBP20；采用qRT-PCR技術，結(jié)合geNorm 和NormFinder軟件對候選的8個內(nèi)參基因在羊草不同組織間的表達穩(wěn)定性進行評價，以期篩選出在羊草各組織相對穩(wěn)定表達的內(nèi)參基因，為羊草研究提供參考。

1 材料與方法

1.1 材料處理

本試驗選用“吉生4號”羊草種子，先用5%次氯酸鈉對種子進行消毒處理，后用無菌水沖洗干凈，種植在20 cm (高)×16 cm(內(nèi)徑)花盆中，在恒溫培養(yǎng)箱(Percival、美國)中進行25 ℃，16 h光照、8 h黑暗培養(yǎng)，正常澆水管理。待生長至抽穗期，選取生長良好的羊草，取不同生長器官(根、莖、葉、穗)，用蒸餾水清洗并擦干，錫紙包裹，液氮速凍后置于-80 ℃超低溫冰箱中保存?zhèn)溆?。每組樣品設3個生物學重復。

1.2 試驗方法

1.2.1 羊草總RNA的提取及cDNA的合成 取100 mg羊草不同組織樣品，置于研缽中，液氮條件下充分研磨，使用Trizol法提取總RNA，1%瓊脂糖凝膠電泳檢測核酸質(zhì)量。取1 μg總RNA，使用Prime ScriptTMRT reagent kit with gDNA Eraser (Perfect Real Time)試劑盒(TaKaRa)反轉(zhuǎn)錄羊草cDNA第一鏈，并將合成的cDNA置于-20 ℃冰箱中保存?zhèn)溆谩?/p>

1.2.2 內(nèi)參基因的定量PCR引物設計 本研究選取了8個常見內(nèi)參(TUA、TUB、18SrRNA、EF-1α、APRT、CYP、Actin和CBP20)作為羊草的候選內(nèi)參基因，其中，Actin序列來源于GenBank(HM623326.1)，其它7個基因均為農(nóng)業(yè)部牧草資源與利用重點實驗室自行擴增,由華大生物技術有限公司測序所得，尚未登錄信息。利用PerlPrimer v.1.1.21軟件，結(jié)合qRT-PCR引物設計原則設計各內(nèi)參基因的qPCR引物。所有引物由華大生物技術有限公司合成(表1)。

1.2.3 內(nèi)參基因的普通PCR擴增 PCR擴增反應體系為20 μL，包括10 μL 2×Taq Master Mix(康為世紀)，10 μmol·L-1的上、下游引物各1 μL，0.5 μL cDNA和7.5 μL RNase-Free Water。反應條件：94 ℃預變性2 min；94 ℃變性30 s，60 ℃退火30 s，72 ℃延伸30 s，35個循環(huán)；72 ℃終延伸2 min。最后用2%瓊脂糖凝膠進行電泳檢測擴增產(chǎn)物。

1.2.4 內(nèi)參基因的熒光定量PCR分析 qRT-PCR反應體系為20.0 μL：SYBR?Select Master Mix(2×) 10.0 μL，上、下游引物(10 mol·L-1)各0.4 μL， cDNA模板1 μL，8.2μL ddH2O，每個樣品技術重復3次。反應在Applied Biosystems StepOnePlusTM實時定量PCR儀(USA)上進行，qRT-PCR反應程序參照試劑盒說明書，95 ℃預變性30 s；95 ℃變性5 s，60 ℃退火延伸30 s，40個循環(huán)；擴增結(jié)束后，在95 ℃獲取熔解曲線，通過分析熔解曲線來確定擴增產(chǎn)物特異性。

表1 8個內(nèi)參基因的定量PCR引物Table 1 Primers used for 8 reference genes in qRT-PCR

1.3 數(shù)據(jù)分析

使用geNorm和NormFinder軟件，對8個常用內(nèi)參基因在羊草不同組織器官的表達穩(wěn)定性進行統(tǒng)計學分析。geNorm軟件是通過分析內(nèi)參基因在不同試驗樣品中的表達穩(wěn)定性值(expression stability value，M)來確定內(nèi)參基因穩(wěn)定性的，其中，M值越小,表示基因表達越穩(wěn)定；反之，M值越大,則基因表達越不穩(wěn)定[19]。NormFinder軟件通過計算不同試驗樣品各內(nèi)參基因的表達穩(wěn)定值(stability value，SV)來確定最優(yōu)的內(nèi)參基因，SV值越小表示該基因表達越穩(wěn)定[19]。

2 結(jié)果與分析

2.1 樣品總RNA提取與質(zhì)量檢測

羊草各組織樣品總RNA的提取采用Trizol法，將提取得各組織總RNA進行1%的瓊脂糖凝膠電泳檢測，所提取的樣品總RNA譜帶清晰，較為完整(圖1)。用NanoDrop2000c檢測總RNA，其A260/A280均在1.8～2.1，A230/A260在2.0左右。以上檢測結(jié)果證實所提取羊草樣品總RNA完整性好，純度高，可進行反轉(zhuǎn)錄合成cDNA。

2.2 內(nèi)參基因引物特異性分析

使用羊草葉片總RNA反轉(zhuǎn)錄cDNA為模板進行PCR擴增，得到的擴增產(chǎn)物進行2%的瓊脂糖凝膠電泳檢測，得到的擴增基因片段特異性良好，片段大小也和預測的結(jié)果一致(圖2)。再分別以羊草葉、莖、根和穗4種組織樣品的cDNA為模板進行qRT-PCR檢測，所有基因均得到明顯的單峰熔解曲線，擴增曲線重復性較好，說明各內(nèi)參基因產(chǎn)物特異性好，不存在引物二聚體，因此證實qRT-PCR的檢測結(jié)果準確。

圖1 羊草不同組織樣品總RNA電泳結(jié)果Fig. 1 Electrophoresis of total RNA from different tissues of Leymus chinensis

圖2 8個內(nèi)參基因的PCR產(chǎn)物Fig. 2 PCR produces of 8 reference genes

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

在qRT-PCR中，基因Ct值(Cycle threshold，循環(huán)閾值)的大小反映其表達量，即Ct值越大，基因表達量越低；Ct值越小，基因表達量越高。本研究中8個內(nèi)參基因Ct值都處于10(18SrRNA)至28(TUB)之間，其中TUB的Ct值最大，在24～30；18SrRNA的Ct值最小，在10～13；TUA、CYP和Actin的Ct值均在20～25；APRT、EF-1α和CBP20的Ct值均處于23～27(圖3)。該結(jié)果表明，8個內(nèi)參基因的表達量存在差異。

圖3 8個內(nèi)參基因Ct值分布圖Fig. 3 Ct values of 8 reference genes

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

2.4.1 geNorm軟件分析 經(jīng)geNorm軟件分析，比較各基因在不同樣品所得表達穩(wěn)定性(M)值，結(jié)果發(fā)現(xiàn)以羊草各組織為試驗材料，其內(nèi)參基因表達穩(wěn)定性存在明顯差異。其中，8個內(nèi)參基因表達在羊草葉、莖、根和穗中穩(wěn)定性由高到底排序分別為(圖4)CYP=TUA>Actin>18SrRNA>APRT>EF-1α>TUB>CBP20；EF-1α=APRT>18SrRNA>TUB>TUA>Actin>CBP20>CYP；CYP=APRT>18SrRNA>TUB>TUA>EF-1α>Actin>CBP20；APRT=Actin>TUB>CYP>CBP20>EF-1α>18SrRNA>TUA。即在羊草葉片表達最穩(wěn)定的內(nèi)參基因是CYP和TUA、最不穩(wěn)定的內(nèi)參基因是CBP20；在莖中表達最穩(wěn)定的內(nèi)參基因是EF-1a和APRT，最不穩(wěn)定的內(nèi)參基因是CYP；在根部表達最穩(wěn)定的內(nèi)參基因是CYP和APRT，最不穩(wěn)定的內(nèi)參基因是CBP20；在穗中表達最穩(wěn)定的內(nèi)參基因是APRT和Actin，最不穩(wěn)定的內(nèi)參基因是TUA。

2.4.2 NormFinder軟件分析 利用NormFinder軟件將所有處理的樣品進行內(nèi)參基因表達分析，通過比較各內(nèi)參基因的表達穩(wěn)定值，得到的結(jié)果如表2所示：在羊草葉、莖、根和穗4種組織中，8個內(nèi)參基因表達穩(wěn)定性由高到底排序分別為APRT=18SrRNA>Actin>EF-1α>CYP>TUA>TUB>CBP20；EF-1α=APRT>18SrRNA>TUB>TUA>CBP20>Actin>CYP；TUB>18SrRNA>APRT>CYP=EF-1α>TUA>Actin>CBP20；CYP>TUB>CBP20>Actin>APRT>EF-1α>18SrRNA>TUA。即在羊草葉片表達最穩(wěn)定的內(nèi)參基因是APRT和18SrRNA；在羊草莖中表達最穩(wěn)定的內(nèi)參基因是EF-1a和APRT；在羊草根部表達最穩(wěn)定的內(nèi)參基因是TUB，其次是18SrRNA和APRT；在羊草穗中表達最穩(wěn)定的內(nèi)參基因是CYP，其次是TUB、CBP20。

圖4 geNorm軟件評價內(nèi)參基因表達穩(wěn)定值Fig. 4 Expression stability values of the reference genes evaluated by geNorm

表2 羊草不同組織器官中8個內(nèi)參基因穩(wěn)定性的NormFinder分析Table 2 Analysis of the stabilities of 8 reference genes in different tissues value of Leymus chinensis by NormFinder

2.5 羊草各組織中最佳內(nèi)參基因的選擇

綜合比較geNorm和NormFinder兩種軟件對羊草葉片、莖、根和穗中8個候選內(nèi)參基因表達穩(wěn)定性評估結(jié)果(表3)，發(fā)現(xiàn)geNorm和NormFinder的評價結(jié)果存在差異，導致內(nèi)參基因的穩(wěn)定性排序不一致。但表達穩(wěn)定性排名前三的基因存在部分重疊，故選擇重疊基因作為每個組織的最佳管家基因。在葉片中選擇Actin作為內(nèi)參基因；在莖中選擇EF-1a作為內(nèi)參基因；在根中選擇APRT和18SrRNA作為內(nèi)參基因；在穗中選擇TUB作為內(nèi)參基因。

3 討論與結(jié)論

近年來大量研究發(fā)現(xiàn)，內(nèi)參基因穩(wěn)定性根據(jù)環(huán)境和品種不同而變化，穩(wěn)定內(nèi)參基因的選擇是基因表達定量分析的關鍵因素。目前，國內(nèi)外關于內(nèi)參基因篩選的報道有很多[20-24]。最初，Actin被認為是在植物各個生長時期及各組織部位均穩(wěn)定表達的內(nèi)參基因[25]。在茶樹[9]和刺槐[11]的不同組織定量PCR穩(wěn)定內(nèi)參篩選中，Actin基因表達最穩(wěn)定。在楊樹(Populus)葉片與根系中Actin也是表達最穩(wěn)定的內(nèi)參基因[26]。Actin基因在鹽生草(Halogetonglomeratus)材料中高度保守且表達量恒定，可作為內(nèi)參[27]。

表3 8個內(nèi)參基因表達穩(wěn)定性geNorm和NormFinder比較Table 3 Comparison of expression stability of 8 reference genes generated by geNorm and NormFinder

有研究證明，常用的傳統(tǒng)內(nèi)參基因也會出現(xiàn)表達不穩(wěn)定的情況。如在矮牽牛(Petuniahybrida)葉與花組織研究中，Actin基因表達不夠穩(wěn)定，EF-1a和CYP被確定為最佳內(nèi)參基因組合[24]。在白菜(Brassicarapassp.pekinensis)[28]、黃瓜(Cucumissativus)[29]內(nèi)參篩選中，EF-1a基因也被證實表達最穩(wěn)定。該基因在小麥不同組織[3]以及在茶樹不同成熟葉片[9]中均表現(xiàn)出較高的表達穩(wěn)定性。這說明EF-1a基因在某些組織研究中比Actin基因更適合選作內(nèi)參基因。由此可見，在實際研究中各內(nèi)參基因均在特定的生長時期、組織或條件下相對表達穩(wěn)定。從本研究對羊草不同組織中8個常用內(nèi)參基因的表達分析可以看出，并沒有一個基因在所有組織中都穩(wěn)定表達。EF-1a在羊草莖組織表達最穩(wěn)定；而Actin基因在羊草葉片組織中表達較穩(wěn)定，在其它組織中表達穩(wěn)定性不是最佳。同樣，18SrRNA基因在羊草根、莖組織中表達相對穩(wěn)定，但在羊草穗中表達穩(wěn)定性較差。有研究認為，在qRT-PCR檢測中，18SrRNA的目標mRNA轉(zhuǎn)錄豐度不精確，故18SrRNA基因不可作為內(nèi)參基因[13]。

本研究中，選用兩個不同的內(nèi)參穩(wěn)定性評價軟件分析得到的結(jié)果存在差異，導致內(nèi)參基因的穩(wěn)定性排序不一致。該現(xiàn)象在其它研究中也有存在，如在大豆[30]，柑橘(Citrusreticulata)[31]，二穗短柄草(Brachypodiumdistachyon)[32]等內(nèi)參基因的研究中，均出現(xiàn)類似結(jié)果，推測可能是分析軟件的統(tǒng)計學原理不同造成了這種差異[30]。綜合數(shù)據(jù)分析，雖然評價結(jié)果存在差異，但表達穩(wěn)定性較高的基因存在部分重疊，故選擇重疊基因作為每個組織的最佳管家基因。即在羊草葉片組織中，Actin基因表達相對穩(wěn)定；莖組織中，EF-1a基因表達最穩(wěn)定；根系組織中，APRT和18SrRNA基因表達相對穩(wěn)定；而在穗組織中，基因TUB表達較穩(wěn)定。本研究結(jié)果將為開展羊草功能基因的表達分析提供參考。

References:

[1] Huggett J,Dheda K,Bustin S,Zumla A.Real-time RT-PCR normalisation;strategies and considerations.Genes Immun,2005,6(4):279-284.

[2] Gachon C,Mingam A,Charrier B.Real-time PCR:What relevance to plant studies.Journal of Experimental Botany,2004,55:1445-1454.

[3] 袁偉,萬紅建,楊悅儉.植物實時熒光定量PCR內(nèi)參基因的特點及選擇.植物學報,2012,47(4):427-436. Yuan W,Wan H J,Yang Y J.Characterization and selection of reference genes for real-time quantitative RT-PCR of plants.Chinese Bulletin of Botany,2012,47(4):427-436.(in Chinese)

[4] Xu H,Bao J,Dai J,Li Y,Zhu Y.Genome-wide identification of new reference genes for qRT-PCR normalization under high temperature stress in rice endosperm.PLoS One,2015,10(11):e0142015.

[5] Paolacci A R,Tanzarella O A,Porceddu E,Ciaffi M.Identification and validation of reference genes for quantitative RT-PCR normalization in wheat.BMC Molecular Biology,2009,10(11):1471-2199.

[6] 魏毅東,陳玉,郭海萍,謝華安,張建福,王宗華.水稻缺素脅迫下實時熒光定量RT-PCR的內(nèi)參基因的選擇.農(nóng)業(yè)生物技術學報,2013,21(11):1302-1312. Wei Y D,Chen Y,Guo H P,Xie H A,Zhang J F,Wang Z H.Selection of reference genes for real-time quantitative RT-PCR in rice(OryzasativaL. ssp.japonica) under nutrient deficiency.Journal of Agricultural Biotechnology,2013,21(11):1302-1312.(in Chinese)

[7] 姜婷,蘇喬,安利佳.多重脅迫下玉米實時定量PCR內(nèi)參基因的篩選與驗證.植物生理學報,2015,51(9):1457-1464. Jiang T,Su Q,An L J.Screening and validation of reference genes of qPCR in maize under multiple stresses.Plant Physiology Journal,2015,51(9):1457-1464.(in Chinese)

[8] 劉偉燦,王騏,周永剛,鄧宇,趙利旦,王興超,靳京,董園園,王南,王法微,陳歡,李曉薇,李海燕.大豆干旱脅迫下miRNA與mRNA熒光定量PCR內(nèi)參基因的篩選.西北農(nóng)林科技大學學報：自然科學版,2016,44(2):61-67. Liu W C,Wang Q,Zhou Y G,Deng Y,Zhao L D,Wang X C,Jin J,Dong Y Y,Wang N,Wang F W,Chen H,Li X W,Li H Y.Selection of reference genes for quantitative polymerase chain reaction of miRNA and mRNA in soybean under drought stress. Journal of Northwest A & F University：Natural Science Edition，2016,44(2):61-67.(in Chinese)

[9] 孫美蓮,王云生,楊冬青,韋朝領,高麗萍,夏濤,單育,駱洋.茶樹實時熒光定量PCR分析中內(nèi)參基因的選擇.植物學報,2010,45(5):579-587. Sun M L,Wang Y S,Yang D Q,Wei C L,Gao L P,Xia T,Shan Y,Luo Y.Reference genes for real-time fluorescence quantitative PCR inCamelliasinensis.Chinese Bulletin of Botany,2010,45(5):579-587.(in Chinese)

[10] 劉圓,王麗鴛,韋康,成浩,張芬,吳立赟,胡娟.不同氮處理茶樹實時定量PCR內(nèi)參基因篩選和驗證.茶葉科學,2016,36(1):92-101. Liu Y,Wang L Y,Wei K,Cheng H,Zhang F,Wu L Y,Hu J.Screening and validation of reference genes for quantitative real-time PCR analysis in tea plant(Camelliasinensis) under different nitrogen nutrition.Journal of Tea Science,2016,36(1):92-101.(in Chinese)

[11] 王金星,張利軍,廖資億,張運根,邱乾棟,孫鵬,孫宇涵,胡瑞陽,盧楠,李云.刺槐實時定量PCR分析中內(nèi)參基因的選擇.林業(yè)科學,2014,50(9):167-172. Wang J X,Zhang L J,Liao Z Y,Zhang Y G,Qiu Q D,Sun P,Sun Y H,Hu R Y,Lu N,Li Y.The selection of reference genes for real-time quantitative PCR normalization in black locust(Robiniapseudoacacia).Scientia Silvae Sinicae,2014,50(9):167-172.(in Chinese)

[12] 樊連梅,王超,劉更森,原永兵.蘋果著色期實時定量PCR內(nèi)參基因的篩選和驗證.植物生理學報,2014,50(12):1903-1911. Fan L M,Wang C,Liu G S,Yuan Y B.Screening and validation of reference genes for real-time fluorescence quantitative PCR during coloring period in apple(Malusdomestica).Plant Physiology Journal,2014,50(12):1903-1911.(in Chinese)

[13] 聞志彬,張明理.松葉豬毛菜干旱脅迫下實時定量PCR內(nèi)參基因的篩選.植物生理學報,2015,51(11):2031-2038. Wen Z B,Zhang M L.Reference gene selection for real-time quantitative PCR inSalsolalaricifoliaunder soil drought stress.Plant Physiology Journal,2015,51(11):2031-2038.(in Chinese)

[14] 焦志軍.大針茅干旱脅迫下qRT-PCR內(nèi)參基因的篩選及PIP2s基因的表達分析.呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學碩士學位論文,2014. Jiao Z J.Selection of reference gene ofStipagrandisin qRT-PCR under drought stress and analysis of expression ofPIP2sgene.Master Thesis.Huhhot:Inner Mongolia Agricultural University,2014.(in Chinese)

[15] 武志娟,王照蘭,韓冰,焦志軍,趙鴻彬.大針茅干旱脅迫下最佳內(nèi)參基因組合的篩選及驗證.中國草地學報,2016,38(4):8-12. Wu Z J,Wang Z L,Han B,Jiao Z J,Zhao H B.Screening and verification of the best combination of reference genes ofStipagrandisunder drought stress.Chinese Journal of Grassland,2016,38(4):8-12.(in Chinese)

[16] Li X,Wu Z,Liu Z,Hou X,Badgery W,Guo H,Zhao Q,Hu N,Duan J,Ren W.Contrasting effects of long-term grazing and clipping on plant morphological plasticity:Evidence from a rhizomatous grass.PLoS One,2015,10(10):e0141055

[17] 宋彥濤,李強,王平,周道瑋,烏云娜.羊草功能性狀和地上生物量對氮素添加的響應.草業(yè)科學,2016,33(7):1383-1390. Song Y T,Li Q,Wang P,Zhou D W,Wu Y N.Response ofLeymuschinensisfunctional traits and aboveground biomass to nitrogen addition in Songnen grassland in northeast China.Pratacultural Science,2016,33(7):1383-1390.(in Chinese)

[18] 郭慧琴,任衛(wèi)波,李平,武自念,萬東莉.2,4-表油菜素內(nèi)酯和赤霉素互作對羊草種子萌發(fā)及幼苗生長的影響.草業(yè)科學,2014,31(6):1097-1103. Guo H Q,Ren W B,Li P,Wu Z N,Wan D L.Effect of epi-brassinosteroid and gibberellin on seed germination and seedling growth ofLeymuschinensis.Pratacultural Science,2014,31(6):1097-1103.(in Chinese)

[19] Wu Z,Tian C,Jiang Q,Li X,Zhuang J.Selection of suitable reference genes for qRT-PCR normalization during leaf development and hormonal stimuli in tea plant (Camelliasinensis).Scientific Reports,2016,6:19748.

[20] 史彩華,張友軍.內(nèi)參基因在昆蟲實時熒光定量PCR中的研究進展.應用昆蟲學報,2016(2):237-246. Shi C,Zhang Y J.Advances in reference gene for real-time quantitative reverse transcription PCR(qRT-PCR)of insects research.Chinese Journal of Applied Entomology,2016(2):237-246.(in Chinese)

[21] Koramutla M K,Aminedi R,Bhattacharya R.Comprehensive evaluation of candidate reference genes for qRT-PCR studiesof gene expression in mustard aphid,Lipaphiserysimi(Kalt).Scientific Reports,2016,6:25883.

[22] Kong Q,Yuan J,Gao L,Zhao S,Jiang W,Huang Y,Bie Z.Identification of suitable reference genes for gene expression normalization in qRT-PCR analysis in watermelon.PLoS One,2014,9(2):e90612.

[23] Zhao X,Zhang X,Guo X,Li S,Han L,Song Z,Wang Y,Li J,Li M.Identification and validation of reference genes for qRT-PCR studies of gene expression inDioscoreaopposita.BioMed Research International,2016,2016:3089584.

[24] Mallona I,Lischewski S,Weiss J,Hause B,Egea-Cortines M.Validation of reference genes for quantitative real-time PCR during leaf and flowers development inPetuniahybrida.BMC Plant Biology,2010,10(4):1471-2229.

[25] An Y Q,McDowell J M,Huang S R,McKinney E C,Chambliss S,Meagher R B.Strong,constitutive expression of theArabidopsisACT2/ACT8 actin subclass in vegetative tissues.Plant Journal,1996,10(1):107-121.

[26] Xu M,Zhang B,Su X,Zhang S,Huang M.Reference gene selection for quantitative real-time polymerase chain reaction inPopulus.Analytical Biochemistry,2011,408(2):337-339.

[27] 馬艷紅,徐先良,汪軍成,任盼榮,楊柯,孟亞雄,李葆春,馬小樂,王化俊.鹽生草Actin基因片段的克隆及表達.草業(yè)科學,2015,32(9):1432-1437. Ma Y H,Xu X L,Wang J C,Ren P R,Yang K,Meng Y X,Li B C,Ma X L,Wang H J.Cloning and expression analysis ofActingene fragment from halophyteHalogetonglomeratus.Pratacultural Science,2015,32(9):1432-1437.(in Chinese)

[28] Qi J,Yu S,Zhang F,Shen X,Zhao X,Yu Y,Zhang D.Reference gene selection for real-time quantitative polymerase chain reaction of mRNA transcript levels in chinese cabbage (BrassicarapaL. ssppekinensis).Plant Molecular Biology Reporter,2010,28(4):597-604.

[29] Wan H,Zhao Z,Qian C,Sui Y,Malik A A,Chen J.Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber.Analytical Biochemistry,2010,399(2):257-261.

[30] 曾文韜,柴春月,竇道龍.適用于大豆實時熒光定量PCR分析的內(nèi)參基因的篩選和驗證.南京農(nóng)業(yè)大學學報,2015,38(5):787-795. Zeng W T,Chai C Y,Dou D L.Selection and validation of reference genes for quantitative RT-PCR analysis in soybean.Journal of Nanjing Agricultural University,2015,38(5):787-795.(in Chinese)

[31] Mafra V,Kubo K S,Alves-Ferreira M,Ribeiro-Alves M,Stuart R M,Boava L P,Rodrigues C M,Machado M A.Reference genes for accurate trancript normalization in citrus genotypes under different experimental conditions.PLoS One,2012,7(2):e31263.

[32] Hong S,Seo P J,Yang M,Xiang F,Park C.Exploring valid reference genes for gene expression studies inBrachypodiumdistachyonby real-time PCR.BMC Plant Biology,2008,8:112.

(責任編輯 王芳)

Selection of reference genes for quantitative real-time PCR ofLeymuschinensisin different tissues

Hu Ning-ning1,2, Guo Hui-qin3, Li Xi-liang1, Kong Ling-qi1, Wu Zi-nian1, Zhang Ji-ze1, Chang Chun1, Wan Dong-li1, Zang Hui1,2, Ren Wei-bo1
(1. Grassland Research Institute, Chinese Academy of Agricultural Sciences, Hohhot 010010, China; 2.Graduate School of Chinese Academy of Agriculture Sciences, Beijing 100081, China; 3.College of Life Science Inner Mongolia Agricultural University, Hohhot 010018, China)

To explore the optimal reference genes in different tissues ofLeymuschinensisfor a quantitative real-time PCR experiment, eight general housekeeping genes includingTUA,TUB, 18SrRNA,EF-1α,APRT,CYP,ActinandCBP20 were selected as candidate reference genes. Using roots, stems, leaves, and spikes ofL.chinensis, we analyzed the 8 candidate reference gene expressions using qRT-PCR and evaluated their expression stability using geNorm and NormFinder software. The results showed that stability differed significantly among the eight candidate reference genes in different tissues. The most stable reference gene wasActinin leaves,EF-1ain stems,APRTand 18SrRNAin root andTUBin spikes. The present study has provided an important reference for analysis the expression of functional genes inL.chinensis.

Leymuschinensis; qRT-PCR; reference gene; tissue

Ren Wei-bo E-mail:rppcaucau@163.com

2016-10-08 接受日期：2016-12-12

國家重點基礎研究計劃973項目(2014CB138804)；國家自然科學基金(31502008)；中央級公益性科研院所基本科研業(yè)務費

胡寧寧(1991-)，女，陜西榆林人，在讀碩士生，主要從事羊草功能基因研究。E-mail：huningning317@163.com

任衛(wèi)波(1979-)，男，山西晉城人，副研究員，博士，主要從事牧草分子生物學研究。E-mail：rppcaucau@163.com

10.11829/j.issn.1001-0629.2016-0510

S543+.9;Q943.2

A

1001-0629(2017)07-1434-08

胡寧寧,郭慧琴,李西良,孔令琪,武自念,張繼澤,常春,萬東莉,臧輝,任衛(wèi)波.羊草不同組織實時定量PCR內(nèi)參基因的篩選.草業(yè)科學,2017,34(7)：1434-1441.

Hu N N,Guo H Q,Li X L,Kong L Q,Wu Z N,Zhang J Z,Chang C,Wan D L,Zang H,Ren W B.Selection of reference genes for quantitative real-time PCR ofLeymuschinensisin different tissues.Pratacultural Science，2017,34(7)：1434-1441.