韓彥卿，鄭潔，武彩娟，元香梅，韓淵懷*

(1.山西農(nóng)業(yè)大學 農(nóng)學院，山西 太谷 030801；2.山西農(nóng)業(yè)大學 農(nóng)業(yè)生物工程研究所，山西 太谷 030801)

谷子(Setariaitalica(L.) P. Beauv.)，去殼后稱為小米，是起源于我國黃河流域的古老作物，由青狗尾草馴化而來，現(xiàn)廣泛栽培于歐亞大陸的溫帶和熱帶地區(qū)[1]。2015年我國谷子栽培面積1 200萬畝，占世界80%，總產(chǎn)量占90%，位居雜糧作物之首。然而，近幾年來，谷子各種病害發(fā)生和流行，嚴重降低了谷子的產(chǎn)量和品質(zhì)。例如，谷子白發(fā)病在大發(fā)生年份可造成產(chǎn)量減產(chǎn)70%以上[2]。由此可見，選育抗病高產(chǎn)優(yōu)良谷子品種是解決谷子免受病害侵襲的首選[3]?；诒狙芯渴姨镩g多年調(diào)查統(tǒng)計發(fā)現(xiàn)，晉谷21雖然品質(zhì)優(yōu)良，但田間表現(xiàn)極易感白發(fā)病。然而其親本晉汾52在田間表現(xiàn)為抗病。這為我們研究2個品種間的抗感分子機制提供了非常好的研究材料。

隨著測序技術(shù)的突飛猛進和發(fā)展，基因組重測序技術(shù)已經(jīng)成為發(fā)掘重要抗病相關(guān)基因的重要手段,已經(jīng)廣泛應(yīng)用到各個研究領(lǐng)域[4～6]。Jia等2013年對916份來自世界各地的谷子核心資源進行了基因組重測序，分析獲得了200多萬個SNP標記和50萬個插入缺失(InDel)的變異圖譜，為谷子的遺傳研究和抗病改良奠定了堅實基礎(chǔ)[7]。施陽選用高抗白粉病片段代換系和高感的受體親本2個黃瓜品種進行全基因組重測序，鑒定出了具有LRR結(jié)構(gòu)域的抗病基因，初步明確了黃瓜抗白粉病機制[6]。本研究基于本實驗室對晉谷52 和晉汾21全基因組重測序的基礎(chǔ)上，依據(jù)已鑒定NBS家族基因的SNP、InDel、CNV和SV位點，對NBS(核苷酸結(jié)合位點)類抗性基因進行比較分析，篩選出6個突變位點較多，接著克隆差異較大的基因并測序分析，擬挖掘鑒定出導致晉谷21及其親本晉汾52抗病性差異的關(guān)鍵基因，旨在為谷子抗白發(fā)病分子輔助育種的培育提供理論指導。

1 材料與方法

1.1 試驗材料

本試驗選用感病谷子品種晉谷21(JG21)，抗病谷子品種晉汾52(JF52)，對照品種豫谷1號(YG1)為試材。3個品種在溫室盆栽種植，谷子生長條件模擬田間條件進行，待谷子長到3～5葉期時，取其幼嫩葉片，立即放入液氮速凍備用。

1.2 試驗方法

1.2.1 谷子基因組DNA提取及PCR擴增

2個品種基因組DNA的提取，采用北京康為世紀 CW0553植物基因組提取試劑盒。所提取DNA的純度和質(zhì)量分別用Nanodrop2000檢測核酸濃度和純度，1.0%瓊脂糖凝膠電泳檢測DNA質(zhì)量。用Primer Premier5.0設(shè)計PCR擴增引物。引物合成由生工生物工程(上海)股份有限公司北京DNA合成部合成。PCR擴增體系：采用50 μL的總反應(yīng)體系，其中2×GC PCR Buffer 25 μL，dNTP Mixture 8 μL，LA Taq酶0.5 μL，正反向引物各2 μL，模板DNA 2 μL ddH2O 10.5 μL。PCR程序設(shè)置為：94 ℃預(yù)變性5 min，94 ℃變性30 sec，退火溫度57～59 ℃不等1 min，72 ℃延伸1～4 min，72 ℃復延伸10 min，循環(huán)38次，4 ℃下保存(表1)。

1.2.2 NBS類型基因比對分析

基于前期基因組重測序及測序拼接的基礎(chǔ)上，本實驗室獲得了完整的NBS類型基因序列。隨后利用序列比對軟件ClustalX對基因的核苷酸和氨基酸序列進行同源比對，尋找差異位點。并對含有突變位點進行多次測序驗證，確認確實是由于輻射誘變導致的突變。蛋白特性分析采用Compute pI/Mw軟件，蛋白的二級結(jié)構(gòu)分析用SOPMA和ScanProsite軟件，三級結(jié)構(gòu)預(yù)測和分析分別軟件SignalP和ProtScale軟件。

表1 引物序列Table 1 Primer sequences

2 結(jié)果與分析

2.1 NBS類基因的擴增片段比較

利用上述設(shè)計的6對PCR引物，對3個不同谷子品種進行PCR擴增，分析擴增條帶發(fā)現(xiàn)，4個基因Si028867、Si033294-1、Si033294-2和Si028382在3個谷子品種中均擴增出單一條帶；對Si027260基因擴增發(fā)現(xiàn)，在YG1號和JG21中擴增出單一清晰的片段，JF52中沒有擴增出條帶；而Si027455中均未擴增出條帶。通過對3個品種的Si027260基因片段擴增結(jié)果我們推測認為，JF52未擴增出條帶可能與其抗病性有關(guān)(圖1)。

2.2 NBS類基因的核苷酸序列比對

基于擴增結(jié)果，我們試圖分析3個品種中均擴增出條帶的2個基因Si028867和Si033294在DNA序列上存在哪些堿基差異。經(jīng)核苷酸序列比對后，在Si028867基因中發(fā)現(xiàn)有5處單核苷酸突變(圖2A)，依次分別是T-C、G-A、C-G、T-C、A-C，其中包含3個位置的轉(zhuǎn)換和兩個位置的顛換，前3個位置的堿基突變發(fā)生在中間部位，位于1 000～1 500 bp之間，剩余的2處堿基突變發(fā)生末端；比對Si033294基因堿基序列發(fā)現(xiàn)，在2個品種中存在1處顛換，為T-G，該突變發(fā)生在基因的起始位置，此外，在基因的260 bp處有一個單堿基A的插入(圖2B)，其余部分序列未發(fā)現(xiàn)突變。

圖1 NBS類基因PCR產(chǎn)物電泳圖譜Fig.1 Electrophoretic map of PCR product of NBS gene

圖2 (2A, 2B)Si028867和Si033294兩個核苷酸序列同源比對Fig.2 (2A and 2B)Homologous comparison of two nucleotide sequence in Si028867 and Si033294A.Si028867；B.Si033294

2.3 NBS類基因的的氨基酸序列分析

為明確JF52和JG21品種的2個NBS類基因Si028867和Si033294差異，對其氨基酸序列同源比對分析。發(fā)現(xiàn)2個品種的Si028867蛋白中有4處氨基酸發(fā)生突變，分別在338、416、704和737位氨基酸處，依次是L-S、L-V、F-L和H-P；發(fā)現(xiàn) 2個品種的Si033294蛋白的氨基酸差異較大，同源性較低，推測可能是由于誘變而來的JG21可能由于該基因中核苷酸片段的插入導致(圖3A和3B)。

2.4 差異基因的蛋白結(jié)構(gòu)分析與預(yù)測

蛋白二級結(jié)構(gòu)特征分析發(fā)現(xiàn)，Si028867蛋白在2個品種中，無論是蛋白特性，還是二級結(jié)構(gòu)特點均沒有明顯差異。而對Si033294的蛋白分析發(fā)現(xiàn)存在較大不同，其在JF52中蛋白為親水性，其二級結(jié)構(gòu)為無規(guī)則卷曲，其次為α-螺旋；在JG21中的蛋白為疏水性，二級結(jié)構(gòu)則與其相反(表2)。對Si028867蛋白的三級結(jié)構(gòu)預(yù)測發(fā)現(xiàn)，JF52和JG21中蛋白的三級結(jié)構(gòu)十分類似。而Si033294蛋白的三級結(jié)構(gòu)分析發(fā)現(xiàn)，2個谷子品種之間的差異較大(圖4)。

2.5 抗性基因的系統(tǒng)進化分析

為了明確抗性基因的遺傳進化，我們將Si028867和Si033294基因編碼的氨基酸序列同水稻、高粱、小麥等禾本科植物蛋白序列比對，發(fā)現(xiàn)這2個蛋白與所選植物蛋白序列的同源性較高，Si028867與比對的幾種植物序列的高度同源區(qū)段相對集中，而Si033294則比較分散，分布在幾個區(qū)段(圖5A和5B)。隨后，我們分別構(gòu)建了Si028867和Si033294進化發(fā)育樹，谷子中Si028867蛋白與玉米的同源蛋白的親緣關(guān)系最近，其次是水稻和高粱；而谷子中Si033294蛋白與谷子和高粱的同源蛋白聚在一起，而其它物種親緣關(guān)系較遠(圖6)。

圖3 (3A，3B)Si028867和Si033 294兩個基因氨基酸序列同源比對Fig.3 (3A and 3B)Homologous comparison of two gene’s amino acid sequence in Si028867和Si033294A.Si028867；B.Si033294

蛋白Protein特性Characteristics二級結(jié)構(gòu)Secondary structure等電點pI分子量Mw(kD)信號肽Signal peptide親水性Hydrophilyα-螺旋/%α-helixβ-轉(zhuǎn)角/%β-turn無規(guī)則卷曲Randomcoil延伸鏈/%ExtendedJF52-Si0288677.0490.95None-0.25526 433.1749 6.1631 539.5716 821.11JG21-Si0288677.0392.24None-0.23127 433.8746 5.6932 339.9316 620.52JF52-Si0332949.44153.48None-0.23340 729.951 047.6558 442.9726 419.43JG21-Si0332949.80155.97None0.02654 440.241 188.7337 527.7431 523.30

圖4 Si028867和Si033294的三級結(jié)構(gòu)預(yù)測Fig.4 Tertiary structure prediction of Si028867and Si033294

3 討論與結(jié)論

植物中的NBS類型基因是抗病基因家族中數(shù)目最龐大的一類家族，抗病基因多成簇分布?？共』虻腖RR(富含亮氨酸的重復序列)位點可與病原菌相互識別。相關(guān)研究表明，LRR區(qū)域的核苷酸突變是導致抗感差異的重要因素[8～10]。Bryan 等研究發(fā)現(xiàn)在水稻NBS-LRR序列中單核苷酸變

異是導致水稻抗/感稻瘟病的原因所在[11]。

我們基于前期全基因組重測序數(shù)據(jù)分析獲得的NBS類型抗病基因，進一步克隆比對篩選，最后發(fā)現(xiàn)2個基因：Si028867和Si033294。經(jīng)對晉谷21、晉汾52的Si028867和Si033294序列比對發(fā)現(xiàn)，Si028867在LRR區(qū)域有5個核苷酸位點發(fā)生突變；對2個品種抗病基因的氨基酸序列比對發(fā)現(xiàn)，氨基酸在4個位置上發(fā)生了非同義突變：T-C、C-G、T-A、A-C。非同義突變很可能造成蛋白結(jié)構(gòu)發(fā)生一定變化，蛋白構(gòu)象變化會直接影響抗病性表現(xiàn)。由此推測氨基酸的非同義突變可能是造成2個親緣關(guān)系較近的品種抗病和感病的原因。這一研究與施陽的研究結(jié)果一致[6]。對另一個基因Si033294分析，其突變也由核苷酸的插入而導致非同義突變，由于該堿基的插入導致2個氨基酸序列發(fā)生了很大差異，這些突變是否導致2個谷子品種對病原菌抗性差異的主要原因有待于進一步研究。

本研究通過對前期重測序結(jié)果分析篩選出的6個NBS類型基因通過鑒定和比較分析，晉谷21和晉汾52兩個品種中發(fā)現(xiàn)Si028867和Si033294兩個基因的突變位點較多，主要表現(xiàn)為非同義突變，二級結(jié)構(gòu)和三級結(jié)構(gòu)中也存在蛋白特征和構(gòu)象變化，推測可能是2個品種間的抗病基因序列發(fā)生非同義突變，是造成抗病性差異的主要原因。本研究結(jié)果為研究谷子抗病分子機理和抗病分子輔助育種提供了理論指導。

圖6 (6A, 6B)Si028867和Si033294的系統(tǒng)發(fā)育進化樹Fig.6 (6A and 6B)Phylogenetic tree of protein in millet and other plant注：左圖為Si028867的系統(tǒng)發(fā)育進化樹；右圖為Si033294的系統(tǒng)發(fā)育進化樹Note：Phylogenetic tree on the left map was Si028867; Phylogenetic tree on the right map was Si033294