王芳,高秋,王杰,馬金星,孫娟*(.青島農(nóng)業(yè)大學經(jīng)濟草本植物應用研究所,山東 青島 6609;.農(nóng)業(yè)部全國草產(chǎn)品質(zhì)量監(jiān)督檢驗測試中心,北京 005)
?
利用SSR標記分析高粱屬種質(zhì)資源的遺傳多樣性
王芳1**,高秋2**,王杰1,馬金星2,孫娟1*
(1.青島農(nóng)業(yè)大學經(jīng)濟草本植物應用研究所,山東 青島 266109;2.農(nóng)業(yè)部全國草產(chǎn)品質(zhì)量監(jiān)督檢驗測試中心,北京 100125)
本研究利用15對SSR引物對161份高粱屬種質(zhì)資源進行了遺傳多樣性分析。結果顯示,15對引物共擴增出118條譜帶,其中88條具有多態(tài)性,多態(tài)性比率(P)為75.21%;多態(tài)性信息含量(PIC值)變幅為0.612~0.806,平均為0.699。161份高粱屬材料遺傳相似性系數(shù)為0.636~0.977,平均基因多樣性指數(shù)(H)為0.696,Shannon信息指數(shù)(I)為1.393。U P G M A聚類分析顯示,161份高粱屬材料在遺傳相似系數(shù)(G S)為0.682處可分為3組;129份高粱和蘇丹草材料在相似系數(shù)為0.671處可分為2組。表明SSR標記可以作為高粱屬種間以及種內(nèi)遺傳分化分析的有力工具,被分析的高粱屬種質(zhì)材料具有較高的遺傳多樣性。
高粱屬;種質(zhì)資源;SSR;遺傳多樣性
http://cyxb.lzu.edu.cn
王芳,高秋,王杰,馬金星,孫娟.利用SSR標記分析高粱屬種質(zhì)資源的遺傳多樣性.草業(yè)學報,2016,25(5):125-133.W A N G Fang,G A O Qiu,W A N G Jie,M A Jin-Xing,S U N Juan.A nalysis of genetic diversity in Sorghum germ plasm collections using SSR m arkers.Acta Prataculturae Sinica,2016,25(5):125-133.
高粱屬(Sorghum)屬于禾本科的須芒草族(A ndropogoneae poaceae)[1-2],全世界大約有29種,主要分布于熱帶到溫帶地區(qū)[3-4]。我國有5種,主要分布在東北、西南以及華南地區(qū)[3]。高粱屬植物多為抗旱作物,在世界各地被廣泛種植,作為糧食、飼料、工業(yè)原料等,具有很高的經(jīng)濟價值[5]。在我國廣泛種植的飼料用高粱屬作物有高粱(Sorghum bicolor)、蘇丹草(Sorghum sudanense)以及高粱和蘇丹草的雜交種高丹草(Sorghum bicolor×Sorghum sudanense),是國家草種質(zhì)資源庫重點收集的種質(zhì)資源之一。國家草種質(zhì)資源庫是草種質(zhì)資源的保存和利用中心,分析種質(zhì)資源庫所保存高粱屬種質(zhì)資源材料的遺傳多樣性分布以及遺傳背景,可以豐富入庫種質(zhì)資源的遺傳信息,有助于發(fā)掘潛在的優(yōu)質(zhì)牧草基因,為育種提供支持,還可了解不同來源的高粱屬資源的遺傳差異,進而提高資源利用效率。
種質(zhì)資源遺傳多樣性分析方法經(jīng)歷了形態(tài)、細胞、生化、分子等發(fā)展階段,目前已有多種分子標記方法被應用到遺傳多樣性評價中[6-7],基于具備多態(tài)性高、重復性好、穩(wěn)定可靠等特點,SS R(sim ple sequence repeat,簡單重復序列)分子標記技術成為其中比較通用的方法[8]。高粱已經(jīng)完成了全基因組測序[9],為開發(fā)利用高粱屬通用SS R引物提供了便利條件,近年來在高粱屬植物中已篩選得到不少擴增效果較好的多態(tài)性引物,現(xiàn)有的研究主要集中在高丹草和蘇丹草一些品系的遺傳分析[10-13],對高粱屬引進以及全國范圍內(nèi)采集的野生、栽培的品種和品系材料的遺傳多樣性分析尚未見報道。
本研究利用SS R標記對國家草種質(zhì)資源庫保存的161份核心高粱屬種質(zhì)材料進行了遺傳多樣性分析,旨在對資源庫中高粱屬種質(zhì)資源的多樣性分布及遺傳背景進行驗證和評價,以期為高粱屬種質(zhì)資源的收集、鑒定、保存及創(chuàng)新利用提供基礎信息。
1.1 試驗材料
供試材料為國家草種質(zhì)資源庫中保存的所有高粱屬材料,共計161份。其中,野生材料12份,品種材料29份,品系材料120份,材料編號及來源地信息見表1。
表1 供試材料Table 1 Plant materials for this study
1.2 樣品制備
2014年9月將試驗材料播種于育苗盤,每份材料保證出苗數(shù)目在100株以上,置于溫室中培養(yǎng)。80%種子的子葉出土2 c m時取樣,每株取0.5 c m左右子葉,每份樣品取100個單株。
1.3 DNA提取
本試驗采用改良C T A B法[14]提取D N A。將材料在液氮中迅速研磨至粉末狀。轉移至已加入預熱(65℃)好的700μL C T A B緩沖液(含10μLβ-疏基乙醇)的離心管中,充分搖勻;置于65℃水浴鍋中,溫育1~2 h,期間每隔15 min搖動一次。加入等體積CI(氯仿∶異戊醇=24∶1),搖動混勻,4℃、12000 r/min離心10 min,取上清液于離心管,重復此步;加入2/3倍體積異丙醇或2倍體積乙醇,-20℃存放1 h;4℃,12000 r/min離心10 min,去上清液;用70%的乙醇清洗沉淀,4℃,12000 r/min離心10 min,去上清液;風干,100μL T E溶解。提取的D N A經(jīng)Nanodrop 2000檢測濃度及質(zhì)量,置于冰箱-20℃長期保存?zhèn)溆?。吸取部分溶液稀釋?0 ng/μL作為工作液。
1.4 PCR擴增及其產(chǎn)物檢測
體系:20μL,其中10×Buffer(含M g2 +)2μL,上游引物(10μm ol/L)與下游引物(10μm ol/L)各加0.5μL,d N T Ps(2.5 m m ol/L)0.4μL,Taq酶(5 U/μL)0.1μL,模板D N A(50 ng/μL)4μL,dd H2O 12.5μL。引物信息見表2,F(xiàn)為上游引物,R為下游引物,由生工生物工程(上海)有限公司合成。所用試劑均為天根科技生化公司產(chǎn)品。
程序:94℃預變性5 min;94℃變性40 s;退火(退火溫度見表2)35 s,72℃延伸45 s,35個循環(huán);72℃延伸5 min,4℃保存。
變性聚丙烯酰胺凝膠電泳:P C R擴增產(chǎn)物中加入6μL Loading buffer(98%甲酰氨,p H 8.0、10 m m ol/L E D T A,0.25%溴酚藍,0.25%二甲苯青),95℃變性5 min,冰上冷卻,在6%的變性聚丙烯酰胺凝膠上電泳50 min左右(80 W恒功率,不含預電泳30 min)。銀染檢測擴增產(chǎn)物。
表2 引物信息及序列Table 2 Primer and sequence information
1.5 數(shù)據(jù)處理
根據(jù)P C R擴增產(chǎn)物在變性聚丙烯酰胺上凝膠電泳的結果,對相同遷移位置或長度一致的條帶進行統(tǒng)計,有帶記為1,無帶記為0,構成標準的0/1矩陣。利用PIC-C A L C 0.6軟件計算位點多態(tài)性含量(poly m orphis m inform ation content,PIC)值。采用N T S Y S 2.10e軟件計算SS R的遺傳相似系數(shù)(G S),并利用U P G M法進行聚類分析,繪制親緣關系樹形圖。
利用Data Trans1.0將構成的0/1矩陣數(shù)據(jù)轉換成基因型數(shù)據(jù),再利用P O P G E N E 1.32軟件計算Shannon值(I,Shannon’s inform ation index)、Nei’s值(H,Nei’s gene diversity)、等位基因數(shù)(Na,observed nu m ber of alleles)與有效等位基因數(shù)(Ne,effective nu m ber of alleles)。
2.1 SS R標記多態(tài)性及性能分析
使用已發(fā)表文獻[11-13,15]中的108對有多態(tài)性的SS R引物進行進一步篩選,得到15對擴增效果好、多態(tài)性高的引物對161份材料進行P C R擴增,分析結果見表3,部分材料的擴增圖譜見圖1。譜帶統(tǒng)計結果顯示,15 對SS R引物共擴增出118條譜帶,其中多態(tài)性譜帶有88條,引物的平均多態(tài)性比率(PPB)為75.21%,各引物擴增譜帶數(shù)目的變幅為6~11,平均每對引物擴增出7.87條譜帶,說明所選引物具有豐富的基因多態(tài)性。不同引物所揭示的供試材料的多態(tài)性信息量(PIC)的變幅為0.612~0.806,均值為0.699,161份高粱屬材料的基因多樣性指數(shù)變幅為0.531~0.827,平均為0.696,有效等位基因數(shù)為2.134~5.777,平均為3.572,Shannon信息指數(shù)為0.824~1.920,平均為1.393,上述統(tǒng)計分析結果表明供試材料遺傳分化豐富。
表3 引物多態(tài)性信息Table 3 The polymorphism of SSR primers
2.2 遺傳相似性分析
根據(jù)擴增結果統(tǒng)計譜帶,計算相似性系數(shù)和遺傳距離。結果顯示161份供試材料的遺傳相似性系數(shù)為0.341,變幅為0.636~0.977。其中,遺傳相似系數(shù)最大的是37和38,達到0.977,其次是26與27、19與21、58 與60、106與107、126與134,達到0.955,表明37與38、26與27、19與21、58與60、106與107、126與134親緣關系很近。登記信息顯示,37來源于北京,38來源于美國,均為高丹草品種樂食;26來源于甘肅,27來源于新疆,均為蘇丹草品種新蘇2號。登記信息與遺傳相似性分析的結果高度一致,不同來源地的同一品種材料聚在了一起,證明所選用引物能夠很好地揭示種質(zhì)材料之間的親緣關系。其他聚在一起的材料由于登記信息不完善,需要開展進一步的田間試驗進行驗證。
圖1 引物C U P57在編號為1~55的高粱屬材料中的擴增結果Fig.1 SSR amplification results of Sorghu m materials coded from 1 to 55 with primer C U P57
2.3 聚類分析
基于遺傳相似性系數(shù),利用U P G M A法構建了161份高粱屬材料的聚類分析圖(圖2)以及129份高粱和蘇丹草的聚類分析圖(圖3)。
在相似性系數(shù)為0.682的水平上,161份材料被分為3組,第1組主要為高丹草以及進口高粱,第2組主要為國內(nèi)的高粱,第3組主要為蘇丹草。其中,編號為123的高粱聚在了蘇丹草一類;編號為124的蘇丹草聚在了高丹草以及進口高粱一類。經(jīng)核查,123和124這兩份材料為同一單位同一批次上交的種子,通過對兩份材料的種子形態(tài)觀察以及與材料上交單位核對,證實材料高粱123應為蘇丹草、蘇丹草124應為高丹草,可能是由于登記、復檢或干燥過程中的失誤導致的。此外,編號為40的高丹草和編號為3的高粱聚在了蘇丹草一支中;編號為22的蘇丹草聚在國內(nèi)高粱一支中,以上3份材料登記信息不完善,還需通過田間試驗進行驗證。高丹草是以高粱為母本、蘇丹草為父本育成的材料,從聚類圖上看其親緣關系更接近高粱。
在129份高粱和蘇丹草的聚類分析中,在相似性系數(shù)為0.671的水平上,供試材料聚成兩大組,第1組所有材料均為高粱(除了存疑材料22),第2組所有材料均為蘇丹草(除了存疑材料3),說明我們選用的SS R引物可較好的將高粱與蘇丹草區(qū)分開。根據(jù)種子入庫登記的信息,種質(zhì)資源庫中的高粱種質(zhì)資源來源廣泛,國內(nèi)來源有東北(吉林、遼寧)、華北(北京、河北、內(nèi)蒙古)、華東(山東、江蘇)、華中(河南、江西)、西北(寧夏、新疆、甘肅)、西南(重慶、四川),同時還收集有不少國外高粱材料。從聚類圖上也可看出高粱種質(zhì)資源地域性強,在高粱這一組中,在相似性系數(shù)為0.693的水平上,55份材料可分為4個亞組(I、II、III、IV),大部分材料(52份)聚在I組和II組中,其中II組主要為來源于東北和國外的高粱材料,來源于國內(nèi)其他5個分布區(qū)的材料主要聚在I組。在I組中,大部分來源于華東、華中和西北的材料也聚成了不同的小支。聚類結果顯示所分析的高粱種質(zhì)資源來源較為豐富,具有較強的地域性。蘇丹草種質(zhì)資源來源比較單一,主要來自甘肅、寧夏、內(nèi)蒙古,在聚類圖中并沒有表現(xiàn)出很強的地域性。
SS R分子標記技術自20世紀90年代出現(xiàn)以來,被各國研究者選用進行多個領域的研究和應用[16-17]。在21世紀初以前由于其高昂的開發(fā)成本和實驗技術的不完善,其優(yōu)越性沒有得到完全發(fā)揮,近十年來D N A測序技術的飛速發(fā)展,在很大程度上解決了SS R分子標記技術遇到的瓶頸問題,因此SS R技術得到了空前的應用和發(fā)展[18-20]。SS R位點具有豐富的多態(tài)性,在品種間甚至于個體間都有變異,適用于進行種內(nèi)遺傳多樣性分析。在種間水平上,由于變異位點兩端保守序列的變異太大,導致一般無法設計通用引物進行分析,限制了其在種間遺傳多樣性分析的應用[21]。本研究利用根據(jù)高粱基因組序列設計的SS R引物,對高粱及其近緣種高丹草和蘇丹草進行了分析,結果表明高粱的SS R引物在高丹草和蘇丹草中具有很好的通用性,可以用來進行這3種材料的遺傳多樣性分析,說明SS R引物不限于種內(nèi)材料的分析,但在多大范圍內(nèi)可以通用還需要對高粱屬其他材料以及高粱屬近緣屬種材料進行分析。詹秋文等[15]使用R A P D標記對高粱屬的兩個物種進行了分析,無法區(qū)分高丹草和蘇丹草,我們選用的SS R引物能夠區(qū)分開高粱和蘇丹草兩個物種,在一定程度上說明SS R分子標記在種間水平的應用可能還有很大的空間。
圖2 161份高粱屬材料的聚類分析Fig.2 Clustering analysis of 161 Sorghum materials
圖3 129份高粱和蘇丹草材料的聚類分析Fig.3 Clustering analysis of 129 Sorghum sudanense and Sorghum bicolor materials
草種質(zhì)資源是進行草品種改良和種質(zhì)創(chuàng)新的基礎,也是環(huán)境保護和生態(tài)建設的重要載體[22],草種質(zhì)資源的遺傳多樣性分析是開發(fā)利用優(yōu)良基因資源、保護基因多樣性的有效手段。全國畜牧總站國家草種質(zhì)資源庫目前已保存草種質(zhì)材料3.2萬份,本研究分析的161份高粱屬材料均來源于該庫。利用SS R分子標記技術對國家草種質(zhì)資源庫的種質(zhì)材料進行遺傳多樣性分析,一是有助于保證種質(zhì)資源信息的準確性,提高入庫種質(zhì)資源的權威性。例如我們通過分析發(fā)現(xiàn)登記為高粱的123號材料和蘇丹草材料聚在一起,登記為蘇丹草的124號材料和高粱材料聚在一起,經(jīng)驗證是登記失誤造成的。還有一些材料沒有聚在合適的位置上,很有可能也是登記錯誤造成的,需要進一步進行驗證,以確保種質(zhì)資源信息的準確性。二是利用遺傳多樣性分析可對種質(zhì)材料的親緣關系進行梳理,為培育新品系提供基礎信息,提高種質(zhì)材料的利用率。目前雜交育種是主流育種方式,育種材料的親緣關系是雜交育種工作的重要指導依據(jù),明晰材料的遺傳背景能夠大大提高育種的成功率。例如遺傳多樣性分析驗證了編號為26和27以及編號為37和38的材料分別為來源地不同的同一品種材料,表明通過SS R標記分析能夠準確識別高粱屬種質(zhì)資源的遺傳距離,提高種質(zhì)資源創(chuàng)新利用效率。此外,遺傳多樣性分析能夠為重點收集提供科學的指導信息。目前草種質(zhì)資源收集的重點已經(jīng)從廣泛收集轉為重點收集,通過對高粱屬種質(zhì)材料的多樣性進行分析,明晰了高丹草和蘇丹草來源比較單一,多樣性不夠高的現(xiàn)狀,可為下一步高粱屬種質(zhì)資源的重點搜集、有效保護工作提供科學指導,有利于進一步明確收集保護的方向。
References:
[1]Clayton W D,Renvoize S A.Genera gra minu m:Grasses of the world:VI.Subfa mily Panicoideae.Kew Bulletin A dditional Series,1986,13(1):256-375.
[2]Liu Q,Peterson P M,Ge X J.Phylogenetic signals in the realized clim ate niches of Chinese grasses(Poaceae).Plant Ecology,2012,212(10):1733-1746.
[3]Chen S L,Phillips S M.Flora of China:Sorghu m M oench[M].Beijing:Science Press & St.Louis:Beijing and Missouri Botanical Garden Press,2006:602-605.
[4]W aston L,Dallwitz M J.The Grass Genera of the W orld[M].O xfordshire:C.B.A.International,1992.
[5]Zhang J Y,Y uan Q H,Zhang W S.The advances in study on forage germ plas m resources and its genetic diversity in China.Grassland of China,2003,25(3):59-65.
[6]W u H,Jian L,Zhu L Q.Current status and advances in Chinese orchids m olecular m arkers research.Biotechnology Bulletin,2010,(9):56-60.
[7]Liang X Y,Ji Y,Bai S Q,etal.Construction of a genetic m ap for chicory using sequence-related a m plified poly m orphis m m arkers.Acta Prataculturae Sinica,2015,24(5):153-158.
[8]Liu X.Scientific evolution of Sorghum taxono m y that was fro m m orphological diversity to D N A diversity.Liaoning A gricultural Science,1999,(5):28-33.
[9]Paterson A H,Bowers J E,Brugg m ann R,etal.The Sorghum bicolor geno m e and the diversification of grasses.Nature,2009,457:551-556.
[10]Li J Q,W ang L H,Zhan Q W,etal.A study on genetic variation of two species of sorghu m by R A P D m arkers.Acta Prataculturae Sinica,2007,16(5):140-144.
[11]Lu J B.The Electronic PC R of SSR Prim ers and Genetic Poly m orphis m in Sorghum bicolor×S.sudanense[D].H efei:A nhui A gricultural U niversity,2010.
[12]Ji M M.Analysis on M ain Agrono mic Characters and SSR of New Strains of Sorghu m-Sudangrass of Super-low H C N Content[D].H ohhot:Inner M ongolia Agricultural University,2011.
[13]Zhu Y Q,Peng J H,Pang L Y,etal.Elite breeding germ plas ms selection and SSR analysis of Sorghum bicolor×S.sudanense induced by space flight.Acta A grestia Sinica,2012,20(6):1150-1155.
[14]Doyle J J,Doyle J L.Isolation of plant D N A fro m fresh tissue.Focus,1990,12:13-15.
[15]Zhan Q W,Li J Q,W ang B H,etal.Establish m ent of D N A fingerprinting for 42 sorghu m and sudangrass accessions and 2 sorghu m-sudangrass hybrids.Acta Prataculturae Sinica,2008,17(6):85-92.
[16]Litt M,Lutyj A.A hypervariable microsatellite revealed by in vitro a m plification of dinucleotide repeats within the cardiac m uscle actin gene.A m erican Journal of H u m an Genetics,1989,44:397-401.
[17]M ccouch S R,Chen X,Panaudo,etal.Microsatellites m arker develop m ent,m apping and applications in rice genetics and breeding.Plant M olecular Biology,1997,35:89-99.
[18]Li L,W ang H G,Zhang X L,etal.A pplication of SSR m olecular m arkers in crop genetic breeding.Journal of Shanxi A gricultural Sciences,2008,36(3):15-18.
[19]Chen F,W ei Z W,Li W M,etal.A nalysis of genetic diversity and population structure in Medicago germ plas m by SSR m arkers.Acta A grestia Sinica,2013,21(4):759-768.
[20]Gao Q,W ang J,M a J X,etal.SSR m olecular m arker technology’s application in grass cultivars for discussion.China Seed Industry,2014,(12):24-27.
[21]Rafalski J A,Vogel J M,M organte M,etal.Generating and using D N A m arkers in plants.In:Birren B,Lai E.Non m a mm alian Geno mic A nalysis:A Practical Guide[M].San Diego:Acade mic Press,1996:75-134.
[22]Chen Z H,Li X F,Y un X J,etal.Diversity and conservation of forage germ plas m resources in China.Pratacultural Science,2009,26(5):1-6.
[5]張吉宇,袁慶華,張文淑.我國牧草種質(zhì)資源及其遺傳多樣性的研究進展.中國草地,2003,25(3):59-65.
[6]武海,蹇黎,朱利泉.中國蘭花資源分子標記鑒定研究進展.生物技術通報,2010,(9):56-60.
[7]梁小玉,季楊,白史且,等.基于SR A P分子標記構建的菊苣遺傳連鎖圖譜.草業(yè)學報,2015,24(5):153-158.
[8]劉欣.高粱屬分類學的科學演變—從形態(tài)多樣性到D N A多樣性.遼寧農(nóng)業(yè)科學,1999,(5):28-33.
[10]李杰勤,王麗華,詹秋文,等.應用R A P D標記對高粱屬兩物種之間遺傳差異的研究.草業(yè)學報,2007,16(5):140-144.
[11]路景標.高丹草SSR引物電子PC R及遺傳多態(tài)性分析[D].合肥:安徽農(nóng)業(yè)大學,2010.
[12]紀明妹.超低氫氰酸高丹草新品系主要農(nóng)藝特性及SSR分析[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學,2011.
[13]朱永群,彭建華,龐良玉,等.高丹草衛(wèi)星搭載材料優(yōu)異種質(zhì)篩選及SSR分析.草地學報,2012,20(6):1150-1155.
[15]詹秋文,李杰勤,汪保華,等.42份高粱與蘇丹草及其2個雜交種D N A指紋圖譜的構建.草業(yè)學報,2008,17(6):85-92.
[18]李麗,王海崗,張曉麗,等.SSR分子標記在作物遺傳育種中的應用.山西農(nóng)業(yè)科學,2008,36(3):15-18.
[19]陳斐,魏臻武,李偉民,等.基于SSR標記的苜蓿種質(zhì)資源遺傳多樣性與群體結構分析.草地學報,2013,21(4):759-768.
[20]高秋,王杰,馬金星,等.SSR分子標記技術在草品種鑒定中的應用探討.中國種業(yè),2014,(12):24-27.
[22]陳志宏,李曉芳,贠旭疆,等.我國草種質(zhì)資源的多樣性及其保護.草業(yè)科學,2009,26(5):1-6.
Analysis of genetic diversity in Sorghum germplasm collections using SSR markers
W A N G Fang1**,G A O Qiu2**,W A N G Jie1,M A Jin-Xing2,S U N Juan1*
1.Applied Economic Herbaceous PlantsInstituteof Qingdao Agricultural University,Qingdao 266109,China;2.NationalQuality Control & Inspection Centre For Grassland Industry Products,M.O.A,Beijing 100125,China
Fifteen SS R m arkers were used to quantify the genetic diversity of the core germ plas m collections of sorghu m.A total of 161 accessions were assayed.118 bands were generated,a m ong w hich 88(75.21%)were poly m orphic.T he average value of poly m orphicinform ation content(PIC)was 0.699(0.612-0.806);the average Nei’s gene diversity index(H)was 0.696;and the average Shannon inform ation index(I)was 1.393.T he genetic similarities a m ong the 161 accessions ranged fro m 0.636 to 0.977.T he U P G M A dendrogra m indicated that the 161 Sorghum germ plas m collections were clustered into three groups with a threshold of 0.682,and that the 129 Sorghum sudanense and Sorghum bicolorx×Sorghum sudanense germ plas m resources were divided into tw o clusters at a threshold of 0.671.T his study de m onstrated that there is considerable genetic diversity conserved in Sorghum germ plas m collections and SS R is an ideal m olecular m arker for assessing this diversity.
Sorghum;germ plas m resource;sim ple sequence repeat(SS R);genetic diversity
.E-m ail:sunjuan@qau.edu.cn
10.11686/cyxb2015336
2015-07-07;改回日期:2015-10-21
國家牧草產(chǎn)業(yè)技術體系項目資助。
王芳(1990-),女,山東濟寧人,在讀碩士。E-m ail:wangfang213097056@163.co m。高秋(1980-),女,云南保山人,畜牧師。E-m ail:gaoqiu1980@g m ail.co m。**共同
T hese authors contributed equally to this work.