張尚雄，荊 照，吳菲菲，欒瑞濤，王建麗，張玉霞，王顯國

(1.中國農業(yè)大學動物科技學院，北京 100193； 2.通遼星圣農業(yè)有限公司，內蒙古 霍林郭勒 029200；3.黑龍江農業(yè)科學院，黑龍江 哈爾濱 624400； 4.內蒙古民族大學農學院，內蒙古 通遼 028000)

飼草高粱作為優(yōu)質牧草，具有產量高、營養(yǎng)豐富、適口性好等優(yōu)點，可作為泌乳奶牛日糧中的玉米(Zeamays)飼草的替代飼料，且不會降低奶牛生產性能[1]，但高粱(Sorghumbicolor)木質化程度高，消化率低是其不足之處[2]。褐色中脈(brown mid-rib,BMR)飼草高粱能有效降低木質素含量，并且顯著提高消化率[3-4]，其難以消化的木質素含量比傳統(tǒng)高粱低40%～60%[5]。目前，褐色中脈飼草高粱改良傳統(tǒng)飼草高粱品質已成為國內外學者的研究熱點之一。有研究報道，BMR飼草高粱比傳統(tǒng)高粱明顯開花早和成熟早，且能提高秸稈的消化率，但籽粒產量顯著低于傳統(tǒng)高粱[6]。BMR飼草高粱年刈割兩次比刈割一次能增加葉莖比，降低莖粗和株高，有效改善BMR品種的飼草品質[7]。BMR-12為BMR飼草高粱的品種之一，具有更低的中性洗滌纖維、酸性洗滌纖維、酸性洗滌木質素含量，相對飼用價值較高，育種潛力較大[8]。BMR飼草高粱不僅可以作為飼草高粱，而且可以作為能源型高粱，有研究表明，BMR飼草高粱低木質素含量與生物乙醇產量無關，但其高葡萄糖含量與生物乙醇正相關[9]。因此，無論作為飼草高粱，還是能源型高粱，褐色中脈極具潛力。

序列相關擴增多態(tài)性(sequence related amplified polymorphism,SRAP)是一種簡單的分子標記系統(tǒng)，具有多態(tài)性高、重復性好，穩(wěn)定可靠等特點[10]。國內外SRAP標記研究植物遺傳多樣性應用廣泛，研究發(fā)現，SRAP能高效地分析垂穗披堿草(Elymusnutans)、老芒麥(E.sibiricus)遺傳變異，并區(qū)分遺傳關系很近的品種。SRAP標記技術的高頻率共顯性以及在基因組中均勻分布的特性優(yōu)于AFLP(amplified fragment length polymorphism)標記，并且相對于SSR(simple sequence repeats)，SRAP特異條帶比例較高，并且更加快捷，且成本低[10-12]。研究種質資源遺傳多樣性有利于解析材料的遺傳背景和親緣關系，對正確評價和有效利用種質資源、服務育種實踐等具有重要意義[13]。為此，本研究對國內首次引進的21份BMR褐色中脈不育系材料進行農藝性狀與遺傳多樣性分析，明確不同種質材料之間的遺傳關系，以期為引進BMR高粱不育系材料的高效利用及種質創(chuàng)新提供依據。

1 材料與方法

21份試驗材料引自美國[14]，均為BMR飼草高粱不育系(表1)。

1.1 供試材料農藝性狀測定方法

試驗地位于內蒙古民族大學，122°27′E，43°62′ N，年均氣溫6～7 ℃，年降水量350～420 mm，土壤類型為沙壤土。試驗小區(qū)隨機排列，各材料播種5行，行距50 cm，行長6 m，株距15 cm，穴播2～3粒，小區(qū)面積為2.5 m×6 m，保苗1.3×105株·hm-2。各材料之間設1 m隔離帶，小區(qū)四周均設保護行。適時進行田間雜草防除，乳熟期測定株高、莖粗、葉寬、莖葉比，采用李源等[8]的方法。

1.2 供試材料遺傳多樣性研究方法

2017年1月16日將試驗材料播種于育苗盤，置于溫室中培養(yǎng)，保證每份材料出苗數10株以上，于2017年2月22日取樣，每份材料取10個單株幼嫩葉片，剪取等量混合，液氮迅速冷凍保存，放于-80 ℃冷凍保存?zhèn)溆谩９┰嚥牧螪NA提取、PCR反應條件參考文獻[15-16]，用1%的瓊脂糖凝膠檢測DNA的質量和濃度。選取試驗材料農藝特性差異較大的5份材料進行引物篩選，且PCR擴增產物用2%瓊脂糖凝膠電泳分離，并用自動凝膠成像系統(tǒng)照相保存。引物信息如表2所列。

1.3 數據統(tǒng)計與分析

將SRAP擴增產物相同遷移位置或長度一致的條帶視為1個位點，統(tǒng)計位點總數和多態(tài)位點數，有帶記為1，無帶記為0。試驗采用R 3.4.1進行供試材料農藝性狀分析，遺傳結構聚類(unweighted pair-group method with arithmetic means,UPGMA)和主成分分析，利用R數據包有ggplot2、cluster、vegan、gclus、scatterplot3d、ape等。

擴增帶總數(total number of bands,TNB)、多態(tài)性條帶數(number of polymorphic bands,NPB)、多態(tài)性比率(percentage of polymorphic bands,PP)、多態(tài)性信息含量(polymorphic information content,PIC)和標記指數(marker index,MI)的計算方法采用顧曉燕等[11]。

2 結果與分析

2.1 21份不育系材料農藝性狀

試驗表明，BMP448(材料3)株高最高，為138.75 cm，BMP450(材料5)最低，為88.00 cm；BMP17(材料16)莖粗最粗、葉寬最寬，分別為29.97 mm和9.88 cm，說明BMP17的莖粗、葉寬特異性突出。BMP20(材料18)莖葉比顯著高于其他供試材料(P<0.05)，為8.69(表3)。將農藝性狀進行聚類分析，表明遺傳距離為0.129時，21份材料可分為3類，第Ⅰ類為材料3，材料1、2、5、6、8、9、12、13、15、17為第Ⅱ類，剩余材料為第Ⅲ類(圖1)。

表1 供試材料Table 1 Experimental materials of 21 sorghum germplasm and their pedigrees

表2 SRAP引物序列Table 2 Primer sequences of SRAP

2.2 SRAP擴增產物的多態(tài)性

130對SRAP引物中，篩選出21對擴增效果好，多態(tài)性較高的引物，共擴增出74個條帶，有多態(tài)性的64條，每個引物平均擴增出3.52個條帶，其中引物的平均多態(tài)性比率為86.5%。21個引物組合的平均PIC值為0.49，其中me13+em4最高，為0.70，而me8+em9的PIC最小，為0.10。引物組合me6+em9標記指數最大，為3.61，而me10+em9的標記指數最小，為0.19(表4)。

2.3 高粱各材料的遺傳相似系數

通過R語言計算供試材料Jaccard遺傳相似系數(genetic similarity coefficient，GS)，結果表明，GS的變化范圍為0.086～0.977，平均GS為0.653。材料BMP447(材料2)和BMP450(材料5)的遺傳相似系數最近，為0.977，親緣關系最近，材料BMP455(材料10)和BMP21(材料19)的遺傳相似系數最小，為0.086，親緣關系最遠。供試材料GS頻率表明，GS值在0.6～0.9的共占46.67%，說明近1/2的供試材料遺傳基礎較窄(圖2)。

同列不同小寫字母表示不同材料間差異顯著(P<0.05)。

Different lowercase letters indicate significant difference among different experimental materials at the 0.05 level.

表4 供試材料SRAP標記引物擴增結果Table 4 Amplification results of SRAP primers

圖2 供試材料遺傳相似系數頻率Fig. 2 Frequency of sample genetic similarity coefficient

2.4 供試材料聚類分析

聚類分析能較好地反映種質材料之間的親緣關系，根據供試材料遺傳距離進行聚類，遺傳距離為0.488時，可將供試材料分為3類，第Ⅰ類為材料19，第Ⅱ類為材料16和20，剩余材料為第Ⅲ類(圖3)。

21份材料聚類后用雙排序矩陣熱圖進行驗證(圖4)，從方塊區(qū)域和聚類樹可以看出兩種分類結果基本一致，不同顏色代表遺傳距離，根據顏色的深淺度可以看出，遺傳距離較小的顏色為深紅色。聚類樹與聚類熱圖組內連接方式不同，導致其分類圖有差異，但對分類結果無影響。

2.5 供試份材料的主成分分析

根據主成分特征值大于1表明，21份材料在遺傳組成上分為3個成分，第Ⅲ主成分為19，特征值方差百分比為4.85%，第Ⅱ主成分為16和20，特征值方差百分比為6.25%，其他材料為第Ⅰ主成分，特征值方差百分比為69.22%。3個主成分特征值方差百分比貢獻率(80.32%)大于80%，因此，本研究結果可用主成分進行分析，并且主成分分析與聚類分析的結果基本一致(圖5)。

2.6 供試份材料農藝性狀與遺傳特性的相關性

21份材料農藝性狀各指標，以及綜合性狀與分子標記的Mantel分析表明，農藝性狀及各指標與分子標記沒有顯著的相關性，其中莖葉比的相關性較高，為0.219(表5)。

圖3 21份高粱種質的SRAP標記聚類結果 Fig. 3 SRAP markers clusters of the 21 sorghum germplasms

圖4 21份供試材料聚類結果重排距離矩陣熱圖Fig. 4 Heat map of distance matrix reorder according to the dendrogram of the 21 sorghum germplasms

1，2，3，…21分別代表材料編號，如表1所列，下圖同。

1, 2, 3, …21, material code, at listed in Table 1, similarly for Fig. 5.

3 討論與結論

飼草高粱農藝性狀受遺傳背景、生長環(huán)境、刈割制度的影響，所以遺傳背景不同的BMR飼草高粱及傳統(tǒng)高粱的農藝性狀會表現出差異，進而影響飼草品質性狀，綜合評價農藝性狀與品質性狀，才能篩選出適合特定生長環(huán)境的優(yōu)質飼草[2,17]。物種的長期進化和適應造成遺傳物質的改變，產生遺傳變異，而遺傳多樣性是遺傳變異高低最直接的表現形式[13]。本研究引進的21份BMR不育系高粱種質譜系表明，親本包含OK11、bmr12、bmr6時可將材料分為兩類，材料19、20為一類，其他材料為另一類，這與分子聚類結果基本相似，由此可見親本差異導致了材料19、20與其他材料的遺傳距離較遠。而材料15與16親本相同，但遺傳分類較遠，系譜來源顯示其培育過程不同，可能是導致其遺傳距離較遠因子之一。因本研究21份BMR不育系高粱種質均由人工長期栽培而成，所以不同親本以及復雜的培育過程使其產生遺傳變異，從而表現出遺傳多樣性。農藝性狀各指標顯著性方差分析中，材料16莖粗、葉寬最大，材料19與材料20莖粗顯著差異(P<0.05)，而分子標記的聚類分析結果中，材料16、19、20與其他材料不同類，因此推斷莖粗、葉寬可能是影響21份不育系高粱種質材料遺傳距離的關鍵因子。飼草高粱不育系生物產量、莖粗、株高、分蘗等農藝性狀受遺傳力加性和非加性效應的雙重影響，但遺傳物質的作用相對較弱，栽培條件及環(huán)境因子起主導作用[2]。由此可見，本研究21份引進材料農藝性狀與分子標記聚類結果幾乎不同，并且Mantel分析也證明二者的相關性較差，這可能與本研究栽培條件和環(huán)境因子有關。本研究遺傳相似系數表明，21份引進材料近一半遺傳基礎較窄，但材料16、19、20的親緣關系較遠，并且這3份材料在莖粗、葉寬等農藝性狀方面有突出表現，因此可以根據分子標記和農藝性狀聚類對現有資源充分、合理利用，在育種實踐中培育出新的種質，拓寬BMR飼草高粱的遺傳基礎，增強BMR種質的利用效率和效果。

遺傳多樣性的探究不僅受物種本身的影響，還有研究方法造成的差異。本研究中21對引物擴增出74條帶，多態(tài)性帶的64條，多態(tài)性比率為86.5%。而采用相同標記方法，高建明等[18]得出甜高粱的多態(tài)性比率為45.5%，李金旺等[19]得出為97.9%，可見不同材料對SRAP標記的特異性表現不同。瓊脂糖電泳和變性聚丙稀酰胺電泳也會影響條帶的顯影，從而影響多態(tài)性分析。

圖5 21份供試材料主成分分析三維排序圖Fig. 5 Three-dimensinal diagram of the 21 samples after PC

分子標記特性GeneticcharactersrP株高Plantheight0.0460.326莖粗Diameterofstem0.0330.380葉寬Widthofleaf0.0140.405莖葉比Stem-leafratio0.2190.107農藝性狀Agronomiccharacters0.0850.240

本研究結果顯示，供試材料中BMP17莖粗、葉寬最大，分別為29.97 mm和9.88 cm。BMP20莖葉比顯著高于其他供試材料(P<0.05)。供試21份褐色中脈不育系高粱種質材料共擴增出74條譜帶，其中64條有多態(tài)性，多態(tài)性比率為86.5%，多態(tài)性信息含量平均值為0.49。遺傳相似系數的變化范圍為0.086～0.977。UPGMA聚類分析顯示，遺傳距離為0.488時，可將供試材料分為3類。這表明，所引進的BMR飼草高粱不育系材料遺傳基礎較窄。

參考文獻References:

[1] Theurer C B,Huber J T,Delgado-Elorduy A,Wanderley R.Invited review:Summary of steam-flaking corn or sorghum grain for lactating dairy cows.Journal of Dairy Science,1999,82:1950-1959.

[2] 呂鑫,平俊愛,張福耀,杜志宏,李慧明,楊婷婷,牛皓,姚琳.新選飼草高粱恢復系農藝性狀配合力效應分析.草業(yè)科學,2016,33(7):1361-1366.

Lyu X,Ping J A,Zhang F Y,Du Z H,Li H M,Yang T T,Niu H,Yao L.Effect analysis on the combining ability of main agronomic traits for new breeding restorer lines derived from forage sorghum.Pratacultural Science,2016,33(7):1361-1366.(in Chinese)

[3] Bucholtz D L,Cantrell J R,Axtell P D,Victor L.Lignin biochemistry of normal and brown midrib mutant sorghum.Journal of Agriculture.1980,28(6):1239-1241.

[4] Pedersen J F,Vogel K P,Funnell D L.Impact of reduced lignin on plant fitness.Crop Science,2005,45(3):812-819.

[5] Casler M D,Pedersen J F,Undersander D J.Forage yield and economic losses associated with the brown-midrib trait in Sudangrass.Crop Science,2003,43(3):782-789.

[6] Umakanth A V,Bhat B V,Blummel M,Aruna C,Seetharama N,Patil J V.Yield and stover quality of brown mid-rib mutations in different genetic backgrounds of sorghum.Indian Journal of Animal Sciences,2014,84(2):181-185.

[7] 許能祥,顧洪如,董臣飛,程云輝,張文潔,丁成龍.刈割對不同品種褐色中脈飼用高粱飼草品質及農藝性狀的影響.草地學報,2014,22(3):623-630.

Xu N X,Gu H R,Dong C F,Cheng Y H,Zhang W J,Ding C L.Effects of cutting on the feed quality and agronomic traits of brown midrib(bmr) forage sorghum varieties.Acta Agrestia Sinica,2014,22(3):623-630.(in Chinese)

[8] 李源,游永亮,趙海明,謝楠,翟蘭菊,劉貴波.褐色中脈飼草高粱品種產量品質與抗逆性分析.草地學報,2014,22(4):889-896.

Li Y,You Y L,Zhao H M,Xie N,Zhai L J,Liu G B.Analysis on the yield,quality and stress resistance of brown midrib forage sorghum varieties.Acta Agrestia Sinica,2014,22(4):889-896.(in Chinese)

[9] Guragaina Y N,Ganeshc K M,Bansala S,Sathishc R S,Rao N,Vadlani P V.Low-lignin mutant biomass resources:Effect of compositional changes on ethanol yield.Industrial Crops and Products,2014,61:1-8.

[10] 苗佳敏,張新全,陳智華,鐘金城,陳仕勇,馬嘯,白史且.青藏高原和新疆地區(qū)垂穗披堿草種質的SRAP及RAPD分析.草地學報,2011,19(2):306-316.

Miao J A,Zhang X Q,Chen Z H,Zhong J C,Chen S Y,Ma X,Bai S Q.SRAP and RAPD analysis ofElymusnutansGriseb. germplasm from the Qinghai-Tibetan Plateau and Xinjiang.Acta Agrestia Sinica,2011,19(2):306-316.(in Chinese)

[11] 顧曉燕,郭志慧,張新全,周永紅,白史且,張昌兵,蔣忠榮,劉新,周朝杰,馬嘯.老芒麥種質資源遺傳多樣性的SRAP分析.草業(yè)學報,2014,23(1):205-216.

Gu X Y,Gu Z H,Zhang X Q,Zhou Y H,Bai S Q,Zhang C B,Jiang Z R,Liu X,Zhou C J,Ma X.Genetic diversity ofElymussibiricusgermplasm resources revealed by SRAP markers.Acta Pratacultural Sinica,2014,23(1):205-216.(in Chinese)

[12] Ahmad R,Potter D,Southwick S M.Genotyping of peach and nectarine ultivars with SSR and SRAP molecular marker.Journal of American Society for Horticultural Science,2004,129(2):204-210.

[13] 倪先林,趙甘霖,汪小楷,劉天朋,陳國民,丁國祥.42份糯高粱種質資源的SSR標記遺傳多樣性分析.分子植物育種,2015,13(9):1962-1969.

Ni X L,Zhao G L,Wang X K,Liu T P,Chen G M,Ding G X.The genetic diversity analysis of 42 Glutinous sorghum germplasm resources using SSR makers.Molecular Plant Breeding,2015,13(9):1962-1969.(in Chinese)

[14] Norris R A.Plant Inventory No.207.Plant Materials Introduced January 1 to December 31,1998,Nos.601817 to 606707.Washington,DC：U.S.Department of Agriculture,Agricultural Research Service，1998.

[15] Li G,Quiros C F. Sequence-related amplified polymorphism(SRAP),a new marker system based on a simple PCR reaction:Its application to mapping and gene tagging inBrassica.Theoretical and Applied Genetics,2001,103(2):455-461.

[16] 陳智華,苗佳敏,鐘金城,馬嘯,陳仕勇,張新全.野生垂穗披堿草種質遺傳多樣性的SRAP研究.草業(yè)學報,2009,18(5):192-200.

Chen Z H,Miao J M,Zhong J C,Ma X,Chen S Y,Zhang X Q.Genetic diversity of wildElymusnutansgermplasm detected by SRAP markers.Acta Agrestia Sinica,2009,18(5):192-200.(in Chinese)

[17] 韓云華,王顯國,呼天明,薛建國.褐色中脈(bmr)高粱農藝性狀與飼用價值.草業(yè)科學,2012,29(10):1596-1603.

Han Y H,Wang X G,Hu T M,Xue J G.Agronomic performance and feeding value of brown midrib forage sorghum.Pratacultural Science,2012,29(10):1596-1603.(in Chinese)

[18] 高建明,羅峰,裴忠有,陳秋玲,魏進招,李子芳,桂枝,孫守鈞.甜高粱重要種質材料的SRAP指紋分析.華北農學報,2010,25(2):93-98.

Gao J M,Luo F,Pei Z Y,Chen Q L,Wei J Z,Li Z F,Gui Z,Sun S J.SRAP finger printing of 47 elite sweet sorghum germplasm.Acta Agriculture Boreali-Sinica,2010,25(2):93-98.(in Chinese)

[19] 李金旺,佟德清,李歐靜,陳秋玲,呂建澎,段霞飛,羅峰,高建明,孫守均,裴忠有.甜高粱含糖量與出汁率的SRAP分子標記.分子植物育種,2017,15(9):3553-3559.

Li J W,Tong D Q,Li O J,Chen Q L,Lyu J P,Duan X F,Luo F,Gao J M,Sun S J,Pei Z Y.Reserch on SRAP molecular marker of sugar content and iuice yield in sorghum.Molecular Plant Breeding,2017,15(9):3553-3559.(in Chinese)