谷子新品種(系)主要農(nóng)藝性狀及莖桿特性的遺傳多樣性分析

屈 洋，宋 慧，劉 洋，王可珍，高小麗，高金鋒，王鵬科，楊 璞，馮佰利

(1. 旱區(qū)作物逆境生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室/西北農(nóng)林科技大學(xué)，陜西 楊凌 712100；2. 寶雞市農(nóng)業(yè)科學(xué)研究院，陜西 岐山 722499；3.安陽(yáng)市農(nóng)業(yè)科學(xué)院，河南 安陽(yáng) 455000)

作物表型性狀是由遺傳物質(zhì)和環(huán)境綜合作用的結(jié)果[1]，是環(huán)境影響、自然選擇和人工干預(yù)產(chǎn)生的穩(wěn)定遺傳變異[2]。作物表型性狀的多樣性對(duì)環(huán)境適應(yīng)性和作物進(jìn)化具有重要作用。谷子(SetariaitalicBeauv.)耐旱、耐瘠薄，主要分布在北方農(nóng)牧交錯(cuò)地區(qū)，是我國(guó)干旱半干旱地區(qū)的主要作物[3]。了解谷子新品種(系)表型性狀的多樣性，研究農(nóng)藝及莖桿特性與產(chǎn)量性狀的關(guān)系，對(duì)提高谷子新品種(系)資源有效利用和谷子育種水平具有重要意義。表型性狀的描述和多樣性分析是作物種質(zhì)資源有效利用的重要方法和途徑[4]，表型性狀的多樣性分析已應(yīng)用于水稻、小麥、玉米、棉花、大豆種質(zhì)資源或品種的分析[5-9]。此外，甜高粱[10]，燕麥[11]，蕎麥[12]，糜子[13]等也進(jìn)行了表型性狀的多樣性分析，以提高作物種質(zhì)資源的利用效率。近年來(lái)，谷子表型性狀的相關(guān)性分析主要集中在產(chǎn)量性狀上[14-15]，例如出谷率、穗粒重、單穗重、千粒重、穗長(zhǎng)、穗粗等。楊成元[16]利用灰色關(guān)聯(lián)度法分析了農(nóng)藝性狀對(duì)產(chǎn)量的影響；黃英杰等[17]對(duì)谷子主要產(chǎn)量構(gòu)成因素進(jìn)行了穩(wěn)定性分析；趙禹凱等[18]對(duì)谷子主要農(nóng)藝性狀進(jìn)行了相關(guān)和通徑分析；王海崗[19]對(duì)谷子核心種質(zhì)進(jìn)行了表型性狀的多樣性分析并對(duì)資源進(jìn)行了綜合評(píng)價(jià)。目前，谷子表型性狀的多樣性分析主要集中在產(chǎn)量性狀上[20-21]，莖桿特性方面的研究較少。本研究以31份谷子新品種(系)為材料，利用表型性狀和莖桿性狀的鑒定方法對(duì)農(nóng)藝和莖桿性狀的遺傳多樣性進(jìn)行分析，揭示農(nóng)藝和莖桿特性與產(chǎn)量性狀之間的關(guān)系，為谷子品種資源的利用提供參考。

1 材料與方法

1.1 試驗(yàn)地點(diǎn)

試驗(yàn)于2016年在寶雞市農(nóng)科院劉家塬試驗(yàn)站(107.5°E， 34.5°N)進(jìn)行。試驗(yàn)地所在地區(qū)為暖溫帶大陸性季風(fēng)型半濕潤(rùn)氣候，年平均氣溫12℃，極端最高溫度36℃左右，極端最低溫度-4℃左右，年平均日照時(shí)數(shù)2064.8 h，平均無(wú)霜期214 d，年平均降水量623.8 mm。

1.2 參試材料與試驗(yàn)設(shè)計(jì)

以31份谷子品種(系)為參試材料(表1)，采用完全隨機(jī)區(qū)組設(shè)計(jì)，3次重復(fù)，小區(qū)面積3.5m×4.6m=16 m2。4月29日播種，密度33萬(wàn)株·hm-2，行距50 cm，株距6 cm。播前施用磷酸二銨600 kg·hm-2。于拔節(jié)期蓋防鳥網(wǎng)，防鳥食。各品種(系)進(jìn)入蠟熟期開始取樣收獲，9月10日完成收獲。

1.3 數(shù)據(jù)采集

調(diào)查每個(gè)小區(qū)的出苗期、拔節(jié)期、抽穗期、成熟期、株高、莖粗、節(jié)長(zhǎng)、苗色、苗葉資、株型、米色、穗型等性狀，參照谷子種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)進(jìn)行[22]，并對(duì)每個(gè)品種(系)的節(jié)長(zhǎng)和莖粗進(jìn)行全節(jié)位測(cè)量。每個(gè)小區(qū)內(nèi)隨機(jī)取樣10株測(cè)量谷子品種(系)各性狀值，各性狀取平均值待分析。

1.4 數(shù)據(jù)處理

分級(jí)性狀主要分析性狀的頻率分布和多樣性指數(shù)，具體分級(jí)標(biāo)準(zhǔn)見表2。利用Excel 2007和SPSS 17.0進(jìn)行描述性和頻次分布統(tǒng)計(jì)，完成對(duì)農(nóng)藝性狀、莖桿特性的最大值、最小值、平均值、標(biāo)準(zhǔn)差、變異系數(shù)和多樣性指數(shù)的計(jì)算。利用Shannon-Weanver遺傳多樣性指數(shù)衡量群體遺傳多樣性大小，計(jì)算公式為：H′=-∑Pi×lnPi[23]。其中H′為遺傳多樣性指數(shù)，Pi為某一性狀第i級(jí)別內(nèi)材料份數(shù)占總份數(shù)的百分比，ln為自然對(duì)數(shù)。利用SPSS 17.0進(jìn)行主成分分析和差異性檢驗(yàn)(Tukey法，P<0.05)，并利用DPS7.05計(jì)算Euclidean(歐式)距離，并用Ward(離差平方和)法進(jìn)行聚類。

2 結(jié)果與分析

2.1 谷子農(nóng)藝性狀的遺傳多樣性分析

2.1.1 分級(jí)性狀的遺傳多樣性 8個(gè)分級(jí)性狀間的多樣性指數(shù)見表3。各分級(jí)性狀間多樣性指數(shù)差異較大，株型、穗型和穗密度頻率在1、3、5、7級(jí)均有分布，離散性較高，且多樣性指數(shù)分別為1.259，1.022和1.237；苗色、苗葉資和花藥色頻率在1、4、7級(jí)有分布，離散程度適中，多樣性指數(shù)分別為0.512、0.689和0.579；谷粒顏色、米色頻率在1和7級(jí)有分布，離散性較差，多樣性指數(shù)最低分別為0.241和0.142。

2.1.2 谷子產(chǎn)量性狀的遺傳多樣性 產(chǎn)量相關(guān)性狀遺傳多樣性指數(shù)較高，變異系數(shù)范圍為3.460～24.852，多樣性指數(shù)范圍為1.539～2.170。單穗重、穗粒重、株高的變異系數(shù)分別為23.794%、24.852%和11.855%，且遺傳多樣性指數(shù)分別為2.170、1.165和2.121，多樣性豐富；穗粗、千粒重的變異系數(shù)為10.980%和20.862%，且遺傳多樣性指數(shù)為2.086和2.025，多樣性較豐富；生育期、出谷率變異系數(shù)為5.135%和3.460%，且多樣性指數(shù)為1.648和1.539，多樣性較差(表4)。

表1 谷子新品種(系)名稱及來(lái)源

表2 谷子8個(gè)分級(jí)性狀及調(diào)查標(biāo)準(zhǔn)

表3 谷子8個(gè)分級(jí)性狀的頻率分布和多樣性指數(shù)

表4 谷子農(nóng)藝性狀的遺傳多樣性

2.2 谷子莖桿特性的遺傳多樣性

2.2.1 谷子基部和首部3節(jié)特性的遺傳多樣性 選取谷子品種(系)的基部和首部的3節(jié)，計(jì)算節(jié)長(zhǎng)和節(jié)粗的遺傳多樣性指數(shù)(表5)。節(jié)長(zhǎng)中，變異系數(shù)變化范圍為23.294%～26.811%，多樣性指數(shù)變化范圍為1.860%～2.184%，首部第3節(jié)變異系數(shù)較大，多樣性指數(shù)最高；基部第2節(jié)變異系數(shù)較低，多樣性指數(shù)較低，多樣性差。節(jié)粗中，變異系數(shù)范圍為16.455%～25.610%，多樣性指數(shù)變化范圍為1.473%～2.023%，首部第1節(jié)變異系數(shù)中等，多樣性指數(shù)較大，多樣性豐富；基部第2節(jié)變異系數(shù)較高，多樣性指數(shù)較低，多樣性差。

2.2.2 谷子基部和首部3節(jié)主成分分析 對(duì)谷子莖桿特性的6個(gè)表型性狀進(jìn)行主成分(PV)分析，前3個(gè)主成分( PV(1)、PV(2)和PV(3))累計(jì)貢獻(xiàn)率達(dá)99.577%(表6)。PV(1)主要包括首部節(jié)長(zhǎng)第1、2、3節(jié)，首部和基部的節(jié)粗第1、2、3節(jié)；PV(2)主要包括首部第2節(jié)節(jié)長(zhǎng)；PV(3)主要包括基部第3節(jié)節(jié)長(zhǎng)。進(jìn)一步以PC1為橫坐標(biāo)，PC2為縱坐標(biāo)，繪制散點(diǎn)圖。不同節(jié)位節(jié)粗具有一定的群集性，不同節(jié)位節(jié)長(zhǎng)位點(diǎn)分布較為分散，但是首部和基部的第2、3節(jié)長(zhǎng)具有一定的群集性(圖1)。

表5 谷子莖桿特性的遺傳多樣性

注：基1、2和3：為基部第一節(jié)、第二節(jié)和第三節(jié)；首1、2和3：為首部第一節(jié)、第二節(jié)和第三節(jié)。

Note：BN(1，2，3)：The first(second，third) node from bottom；TN(1，2，3)：The first(second，third) node from top.

表6 基部和首部莖桿特性的前3個(gè)主成分的特征向量、主成分值，貢獻(xiàn)率和累積貢獻(xiàn)率

注：基1、2和3：為基部第一節(jié)、第二節(jié)和第三節(jié)；首1、2和3：為首部第一節(jié)、第二節(jié)和第三節(jié)。

Note：BN(1，2，3)：the first(second，third) node from bottom；TN(1，2，3)：the first(second，third) node from top.

注：BNL(1、2和3)：基部第一節(jié)、第二節(jié)和第三節(jié)長(zhǎng)；TNL(1、2和3)：首部第一節(jié)、第二節(jié)和第三節(jié)長(zhǎng)；BNW(1、2和3)：基部第一節(jié)、第二節(jié)和第三節(jié)粗；TNW(1、2和3)：首部第一節(jié)、第二節(jié)和第三節(jié)粗。Note：BNL(1，2，3)：the first(second，third) node length from bottom；TNL(1，2，3)：the first(second，third) node length from top；BNW(1，2，3)：the first(second，third) node wide from bottom；TNW(1，2，3)：the first(second，third) node wide from top.圖1 谷子6個(gè)莖桿特性的主成分分析Fig.1 Principal component analysis of 6 node traits

2.3 基于莖桿特性的谷子品種多樣性分析

2.3.1 基于莖桿特性的主成分分析 根據(jù)測(cè)量節(jié)位的節(jié)長(zhǎng)和節(jié)粗對(duì)不同品種進(jìn)行主成分分析，見表7。前4個(gè)主成分的特征值均大于1，累計(jì)貢獻(xiàn)率達(dá)到79.84%。PV(1)的貢獻(xiàn)率最大(28.20%)，其次是PV(2)，最后是PV(5)。品種1對(duì)PV(1)貢獻(xiàn)較大，品種25對(duì)PV(2)貢獻(xiàn)較大，品種18對(duì)PV(3)的貢獻(xiàn)率較大，品種2對(duì)PV(4)的貢獻(xiàn)率較大。進(jìn)一步對(duì)31個(gè)對(duì)谷子品種(系)計(jì)算歐式距離，利用離差平方和法進(jìn)行聚分析，如圖2可將谷子品種分為兩類。第一類包括品種1、16、18、5、29、17、6、10、2、9、30、26、27、12；第二類包括品種3、19、14、15、23、4、22、7、13、25、31、8、20、21、24、28，分類結(jié)果說明31個(gè)谷子新品種(系)在莖桿特性上存在差異，且不同品種之間莖桿特性存在一定的相似性，可作為31份谷子新品種(系)莖桿特性及抗倒伏能力差異的依據(jù)。

表7 不同品種莖桿特性的主成分分析

圖2 谷子品種莖桿特性的聚類分析Fig.2 Cluster analysis of node traits

2.3.2 不同品種莖桿多樣性與產(chǎn)量性狀的分析 對(duì)參試材料的莖桿不同節(jié)位節(jié)粗和節(jié)長(zhǎng)計(jì)算變異系數(shù)和多樣性指數(shù)，并與單穗粒重進(jìn)行相關(guān)分析，見表8。單穗粒重與莖桿粗H′和莖桿長(zhǎng)CV呈正相關(guān)(P<0.05)，且與莖桿粗CV和莖桿長(zhǎng)H′呈負(fù)相關(guān)(P<0.05)；莖桿粗H′與莖桿粗CV呈負(fù)相關(guān)(P<0.05)，莖桿長(zhǎng)CV與莖桿長(zhǎng)H′呈負(fù)相關(guān)(P<0.05)。利用單穗粒重與莖桿特性進(jìn)行回歸分析(圖3)，莖桿長(zhǎng)H′越高，單穗粒重越低(圖3A)，莖桿長(zhǎng)CV越大，單穗粒重越大(圖3B)，莖桿粗H′越高，單穗粒重越高(圖3C)，莖桿粗CV越高，單穗粒重越低(圖3D)。

3 討 論

3.1 谷子農(nóng)藝性狀的多樣性

谷子農(nóng)藝性狀的多樣性研究主要集中在產(chǎn)量性狀[24]，例如單穗粒重、穗長(zhǎng)、穗粗等，這些性狀共同構(gòu)成產(chǎn)量的主要因素，其遺傳變異的信息可以作為親本材料選擇的依據(jù)[25]。本研究中31個(gè)谷子新品種(系)的8個(gè)分級(jí)性狀和8個(gè)農(nóng)藝性狀表現(xiàn)出豐富的變異，各性狀的變化幅度、變異系數(shù)和多樣性指數(shù)均較高，其中單穗重和穗粒重與產(chǎn)量密切相關(guān)，變異系數(shù)為23.794%和24.852%，多樣性指數(shù)分別為2.170和2.165，穗型和穗密度的多樣性指數(shù)為1.022和1.237，說明31個(gè)谷子新品種(系)農(nóng)藝性狀變異豐富，多樣性較高，可以進(jìn)一步在育種中進(jìn)行應(yīng)用，提高谷子新品種(系)的效率。然而，31個(gè)新品種(系)的群體小，遺傳基礎(chǔ)狹窄，而且谷子產(chǎn)量等相關(guān)性狀的變異受環(huán)境影響較大[26]，不同的產(chǎn)區(qū)應(yīng)根據(jù)農(nóng)藝性狀的實(shí)際水平進(jìn)行種質(zhì)改良和創(chuàng)制，有條件的產(chǎn)區(qū)可引入外來(lái)種質(zhì)資源，拓寬遺傳基礎(chǔ)。

表8 莖桿多樣性與產(chǎn)量性狀的相關(guān)性分析

注：*代表在P<0.05水平差異顯著，下同。

Note:* A significant level atP<0.05， the same as below.

圖3 莖桿長(zhǎng)和粗與單穗粒重的相關(guān)性分析(A：莖桿長(zhǎng)H′，B：莖桿長(zhǎng)CV， C： 莖桿粗H′，D:莖桿粗CV)Fig.3 Seed weight and genetic diversity of node from different varieties(A:H′ of node length， B: CV of node length， C: H′ of node wide， D: CV of node wide)

3.2 谷子莖桿性狀的多樣性

作物表型性狀的多樣性分析包括物候期、產(chǎn)量構(gòu)成因素和植株性狀[27-28]等，而對(duì)莖桿相關(guān)方面的研究重要集中在新種質(zhì)創(chuàng)新上[29]。本研究中精確測(cè)量谷子新品種(系)不同節(jié)位節(jié)長(zhǎng)和節(jié)粗，節(jié)長(zhǎng)中平均變異系數(shù)25.586%，變化幅度較大，平均多樣性指數(shù)1.972，多樣性豐富；節(jié)粗中，平均變異系數(shù)20.550%，平均多樣性指數(shù)1.798，多樣性豐富；不同節(jié)位節(jié)粗群集性較強(qiáng)，節(jié)長(zhǎng)群集性較差；不同品種的莖桿性狀對(duì)31個(gè)谷子品種(系)進(jìn)行聚類，可將谷子群體分為兩類，說明谷子莖桿性狀具有豐富的遺傳多樣性，且不同谷子新品種(系)在莖桿特性上具有一定的相似性。

3.3 谷子莖桿特性與產(chǎn)量性狀

作物產(chǎn)量是由品種、環(huán)境、栽培管理和病蟲害防控等因素綜合作用的結(jié)果，而品種是影響產(chǎn)量的重要因素[30-31]。在品種農(nóng)藝性狀中，地上部分的植株性狀和地下部分的根系作為重要的產(chǎn)量影響因素常被用作新品種選育的重要參考指標(biāo)[32-33]，其中株高、穗長(zhǎng)、穗粗、穗密度、出谷率等是谷子育種中進(jìn)行單株選擇的重要依據(jù)[34]。本研究對(duì)產(chǎn)量主要性狀單穗粒重與莖桿特性進(jìn)行相關(guān)分析，單穗粒重與莖桿特性的H′和CV存在顯著的相關(guān)關(guān)系(P<0.05)，單穗粒重與莖桿粗H′和莖桿長(zhǎng)CV呈正相關(guān)，單穗粒重與莖桿粗CV和莖桿長(zhǎng)H′呈負(fù)相關(guān)。這些結(jié)果說明莖桿特性也是影響谷子產(chǎn)量構(gòu)成的重要因素之一，區(qū)別于以往農(nóng)藝性狀與產(chǎn)量關(guān)系的研究[35-36]，谷子莖桿特性在種質(zhì)創(chuàng)制和育種實(shí)踐中應(yīng)給于足夠的重視。值得注意的是31個(gè)谷子新品種(系)群體莖桿粗H′和CV呈負(fù)相關(guān)(r=-0.40，P<0.05)，莖桿長(zhǎng)H′和CV呈負(fù)相關(guān)(r=-0.41，P<0.05)，可能與莖桿的遺傳基礎(chǔ)存在一定的關(guān)聯(lián)，其內(nèi)在的機(jī)制還需進(jìn)一步的研究。

4 結(jié) 論

31個(gè)谷子新品種(系)8個(gè)分級(jí)性狀(苗色、幼苗葉資、株型、花藥顏色、穗型、穗碼密度、籽粒顏色和米色)和8個(gè)農(nóng)藝性狀(生育期、株高、穗長(zhǎng)、穗粗、單穗重、穗粒重、出谷率和千粒重)變異幅度較大，群體多樣性豐富，依據(jù)莖桿特性可將谷子新品種(系)劃分為兩大類群，且莖桿長(zhǎng)和粗具有較高的遺傳變異水平，莖稈粗H′和莖稈長(zhǎng)CV與單穗粒重呈正相關(guān)，莖稈粗CV和莖稈長(zhǎng)H′與單穗粒重呈負(fù)相關(guān)，可作為未來(lái)谷子抗倒伏和高產(chǎn)育種重要農(nóng)藝性狀的補(bǔ)充。

參考文獻(xiàn)：

[1] 李斌， 顧萬(wàn)春， 盧寶明. 白皮松天然種群種實(shí)性狀表型多樣性研究[J]. 生物多樣性， 2002，10(2):181-188.

[2] 尹明宇，姜仲茂，朱緒春，等.內(nèi)蒙古山杏種群表型變異[J].植物生態(tài)學(xué)報(bào)，2016，40(10):1090-1099.

[3] Jia G Q，Huang X，Zhi H，Diao X，et al. A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet (Setariaitalica) [J].Nature Genetic, 2013, 45: 957-961.

[4] 李潤(rùn)枝，陳晨，張培培．我國(guó)燕麥種質(zhì)資源與遺傳育種研究進(jìn)展[J]．現(xiàn)代農(nóng)業(yè)科技，2009，17：44-45.

[5] 代攀虹，孫君靈，何守樸，等. 陸地棉核心種質(zhì)表型性狀遺傳多樣性分析及綜合評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)科學(xué)，2016，49(19):3694-3708.

[6] 陳雪燕，工亞娟，雒景吾，等. 陜西省小麥地方品種主要性狀的遺傳多樣性研究[J].麥類作物學(xué)報(bào)，2007，27(3):456-460.

[7] 王果，胡正，張保缺，等.山西省野生大豆資源遺傳多樣性分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2008，41(7):2182-2190.

[8] 胡標(biāo)林，萬(wàn)勇，李霞，等.水稻核心種質(zhì)表型性狀遺傳多樣性分析及綜合評(píng)價(jià)[J].作物學(xué)報(bào)，2012，38(5):829-839.

[9] 蔡一林，劉志齋，王天宇，等.國(guó)內(nèi)部分玉米地方品種的品質(zhì)與農(nóng)藝性狀的表型多樣性分析[J].植物遺傳資源學(xué)報(bào)，2011，12(1)：31-36.

[10] 王繼師，劉祖昕，樊帆，等.24個(gè)甜高粱品種主要農(nóng)藝性狀與品質(zhì)性狀遺傳多樣性分析[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2012，17(6)：83-91.

[11] 南銘，馬寧，劉彥明，等.燕麥種質(zhì)資源農(nóng)藝性狀的遺傳多樣性分析[J].干旱地區(qū)農(nóng)業(yè)研究，2015，33(1)：262-267.

[12] 高金鋒，張慧成，高小麗.西藏苦蕎種質(zhì)資源主要農(nóng)藝性狀分析[J].河北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2008，31(2):1-5.

[13] 董俊麗，王海崗，陳凌，等．糜子骨干種質(zhì)遺傳多樣性和遺傳結(jié)構(gòu)分析[J]．中國(guó)農(nóng)業(yè)科學(xué)，2015，48(16):3121-3131.

[14] 欒素榮，王占廷，李青松.谷子產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].河北農(nóng)業(yè)科學(xué)，2010，14(11):115-116，118.

[15] 趙禹凱，王顯瑞，張立媛，等.谷子產(chǎn)量與主要農(nóng)藝性狀的灰色關(guān)聯(lián)度分析[J].吉林農(nóng)業(yè)科學(xué)，2014，39(2)：9-12.

[16] 楊成元.春谷產(chǎn)量與其相關(guān)性狀的關(guān)聯(lián)度分析[J]. 雜糧作物，2002，22(5):259-261.

[17] 黃英杰，張巖.谷子品種產(chǎn)量及主要產(chǎn)量構(gòu)成因素穩(wěn)定性的分析[J].作物雜志，2002，5:43-44.

[18] 趙禹凱，王顯瑞，陳高勛，等.谷子主要農(nóng)藝性狀的相關(guān)和通徑分析[J].內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)，2014，35(2)：35-38.

[19] 王海崗，賈冠清，智慧，等.谷子核心種質(zhì)表型遺傳多樣性分析及綜合評(píng)價(jià)[J].作物學(xué)報(bào)，2016，42(1):19-30.

[20] 刁現(xiàn)民. 基礎(chǔ)研究提升傳統(tǒng)作物谷子和黍稷的科研創(chuàng)新水平[J].中國(guó)農(nóng)業(yè)科學(xué)，2016，49(17):3261-3263.

[21] 孟慶立，關(guān)周博，馮佰利，等.谷子抗旱相關(guān)性狀的主成分與模糊聚類分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2009，42(8):2667-2675.

[22] 陸平.谷子種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M].北京：中國(guó)農(nóng)業(yè)出版社，2006.

[23] 屈洋，周瑜，王釗，等.苦蕎產(chǎn)區(qū)種質(zhì)資源遺傳多樣性和遺傳結(jié)構(gòu)分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2016，49(11)：2049-2062.

[24] 田伯紅.谷子地方品種和育成品種的遺傳多樣性研究[J].植物遺傳資源學(xué)報(bào)，2010，11(2):224-228.

[25] 楊天育，黃相國(guó)，何繼紅，等.谷子遺傳資源多樣性研究進(jìn)展[J].西北農(nóng)業(yè)學(xué)報(bào)，2003，12(1):43-47.

[26] Qu Y， Zhou Y， Chao G， et al. Water use efficiency of foxtail millet (PanicumitalicumL.) under climate change conditions in Northwest regions of China[J]. Agrociencia，2016, 50: 665-676.

[27] 刁現(xiàn)民，程汝宏.十五年區(qū)試數(shù)據(jù)分析展示谷子糜子育種現(xiàn)狀[J].中國(guó)農(nóng)業(yè)科學(xué)，2017，50(23)：4469-4474.

[28] 張艾英，郭二虎，刁現(xiàn)民，等.2005-2015年西北春谷中晚熟區(qū)谷子育成品種評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)科學(xué)，2017，50(23)：4486-4495.

[29] 劉正理，程汝宏，黃文勝，等.甜稈多穗型超早熟谷子新種質(zhì)的創(chuàng)新[J].中國(guó)農(nóng)業(yè)科學(xué)，2005，38(1):17-21.

[30] 劉正理，程汝宏，李香月，等.不同類型谷子品種葉片對(duì)產(chǎn)量形成效應(yīng)的研究[J].中國(guó)農(nóng)學(xué)通報(bào)，1995，11(3)：12-16.

[31] 夏雪巖，師志剛，張婷，等.不同粒色谷子品系的鳥害程度研究[J].河北農(nóng)業(yè)科學(xué)，2014，18(2)：4-6.

[32] 張文英，智慧，柳斌輝，等.豫谷1號(hào)和青狗尾草RIL群體根系變異和垂直分布[J].作物學(xué)報(bào)，2014，40(10):1717-1724.

[33] 劉正理，程汝宏.谷子目標(biāo)性狀基因庫(kù)育種技術(shù)體系的構(gòu)建及其應(yīng)用[J].中國(guó)農(nóng)業(yè)科學(xué)，2005，38(7):1306-1311.

[34] 程汝宏，籍貴蘇，李香月.夏谷高產(chǎn)育種單株選擇模式[J].華北農(nóng)學(xué)報(bào)，1996，11(1)：123-127.

[35] 陳素省，趙國(guó)順，郝艷敏.應(yīng)用灰色關(guān)聯(lián)度法評(píng)價(jià)谷子新品種[J].河北農(nóng)業(yè)科學(xué)，2006，10(2)：62-65.

[36] 劉正理，程汝宏，張鳳蓮，等.華北夏谷區(qū)主要谷子品種遺傳基礎(chǔ)分析及改良途徑的探討[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2006，21(增刊)：84-90，99.