谷子（Setaria italica Beauv.）作为中国农耕文化的代表作物，具有耐贫瘠、抗旱及广适等特性，其籽粒富含膳食纤维、维生素、酚类、活性肽和类胡萝卜素等功能活性成分，具有抗氧化、降血糖、降血脂和降血压的保健功能[1-2]，不仅是旱作生态农业的理想栽培作物之一[3]，而且是独特的“食药两用”营养优势作物。我国谷子的优势产区主要在内蒙古东部、河北及山西大部、陕西北部、甘肃东部和宁夏中部等地，种植面积常年稳定在133万hm2左右[4]，约占世界的90%以上；此外，谷子的秸秆和籽粒具有其他谷类不具备的饲用价值优势[5]，常作为优质饲草在农牧交接地区广泛种植。近年来，谷子作为低耗水肥的环境友好型作物，迅速发展成为抗旱、耐逆及C4光合途径研究的模式作物[6-7]。
新石器时代谷子逐渐成为我国主栽作物[8]，陕西、山西及甘肃东部一带是其栽培品种的起源中心[9]，经过长期的自然选择、人工驯化和外扩传播，形成了遍及世界各地的丰富遗传资源。谷子种质资源作为种质创新和遗传改良的物质基础，摸清其遗传本底是谷子育种突破的关键[10]。目前，利用形态性状[11]、品质性状[9,12]、细胞学标记[13]和分子标记[14]等手段，开展了谷子种质的遗传多样性、抗旱抗病性鉴定评价与筛选优异种质研究[15-16]。随着分子生物学技术的发展，也相继开展了谷子的抗锈病基因鉴定、分蘖相关QTL定位、耐旱性miRNAs分析及QTL发掘等工作[17⇓-19]。这些工作的开展对谷子种质资源科学认识和优异种质挖掘具有实际意义。
1 材料与方法
1.1 试验材料
供试材料为中国农业科学院作物科学研究所提供的谷子核心种质资源120份，编号G1~G120（表1）。
1.2 试验地概况
试验于2020-2021年在甘肃省定西市农业科学研究院科技创新试验地（104°42′ E，35°32′ N）进行，海拔1920m。日照时数2437.7h，年均气温6.4℃，无霜期140d，年均降水量400mm。试验采用顺序排列，不设重复，每份材料种植2行，行长3m，理论留苗密度为37.5万株/hm2，2020-2021年谷子生育期间平均气温和降水量未发生异常现象（图1）。
图1
图1
谷子种质生育期间降水量及气温状况
Fig.1
Precipitation and temperature during the growth and development of foxtail millet
1.3 试验方法
表2 谷子种质资源8个质量性状的描述标准
Table 2
| 性状 Trait | 描述标准 Description standard |
|---|---|
| 幼苗叶色LC | 1.绿色，2.黄绿，3.紫绿 |
| 叶鞘色LSC | 1.绿色，2.红色，3.紫色 |
| 幼苗叶姿SL | 1.上举，2.半上举，3.平展，4.下披 |
| 开花期叶姿SLF | 1.上举，2.半上举，3.平展，4.下披 |
| 穗颈形状SPN | 1.直立，2.中弯，3.弯曲，4.勾形 |
| 穗形PT | 1.鸡嘴形，2.纺锤形，3.圆筒形，4.棍棒形，5.鸭嘴形，6.猫爪形，7佛手形 |
| 穗松紧度CS | 1.松，2.中，3.紧 |
| 谷粒色HC | 1.白色，2.黄色，3.橙色，4.红色，5.青色，6.褐色，7.黑色 |
LC：叶色，LSC：叶鞘色，SL：幼苗叶姿，SLF：开花期叶姿，SPN：穗颈形状，PT：穗形，CS：穗松紧度，HC：谷粒色，下同
LC: leaf color, LSC: leaf sheath color, SL: shape of leaf, SLF: shape of leaf in flowering, SPN: shape of panicle neck, PT: panicle type, CS: compactness of spike, HC: hull color, the same below
1.4 数据处理
利用Microsoft Excel 2007对试验数据进行整理，计算各个性状的变异系数和遗传多样性指数（H′），采用Shannonʼs遗传多样性指数进行描述和评价遗传多样性指数，其计算公式为：
2 结果与分析
2.1 谷子种质的变异分析及表型多样性分析
120份谷子种质8个质量性状的H′范围为0.32~1.26，变异系数范围为23.99%~58.20%（表3）；幼苗叶色以绿色为主（频数为0.908），叶鞘色以绿色为主（频数为0.725），幼苗叶姿以半上举和平展为主（频数分别为0.416和0.425），开花期叶姿以半上举和下披为主（频数分别为0.391和0.383），穗颈形状以弯曲和勾形为主（频数分别为0.384和0.475），穗形以纺锤形为主（频数为0.317），穗松紧度以中和紧为主（频数分别为0.464和0.383），谷粒色以黄色为主（频数为0.716），表明供试材料的质量性状变异丰富。
表3 参试材料8个质量性状的遗传变异及遗传多样性
Table 3
| 性状 Trait | 遗传多样性指数 H′ | 变异系数 CV (%) | 频率分布Frequency distribution | ||||||
|---|---|---|---|---|---|---|---|---|---|
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | |||
| 幼苗叶色LC | 0.32 | 34.46 | 0.908 | 0.051 | 0.041 | － | － | － | － |
| 叶鞘色LSC | 0.39 | 53.68 | 0.725 | 0.141 | 0.134 | － | － | － | － |
| 幼苗叶姿SL | 0.96 | 25.28 | 0.075 | 0.416 | 0.425 | 0.084 | － | － | － |
| 开花期叶姿SLF | 1.03 | 42.61 | 0.134 | 0.391 | 0.167 | 0.308 | － | － | － |
| 穗颈形状SPN | 1.08 | 23.99 | 0.033 | 0.108 | 0.384 | 0.475 | － | － | － |
| 穗形PT | 1.26 | 58.20 | 0.267 | 0.317 | 0.116 | 0.250 | 0.017 | 0.025 | 0.008 |
| 穗松紧度CS | 1.01 | 30.99 | 0.150 | 0.464 | 0.383 | － | － | － | － |
| 谷粒色HC | 1.11 | 55.85 | 0.084 | 0.716 | 0.034 | 0.041 | 0.050 | 0.041 | 0.034 |
120份谷子的13个数量性状的H′范围为0.63~2.02，变异系数范围为7.64%~ 64.12%（表4）。其中，单株干草重、分蘖粒重和分蘖穗重变异系数均高于60%，分别为64.12%、62.83%和62.51%，说明植株分蘖相关性状在当地生态环境条件下表现出更高的变异性，进行饲草型谷子遗传改良的潜力较大；主穗长、主穗重、主穗粒重和主穗粒数变异系数分别为47.03%、32.75%、36.96%和34.10%，说明与籽粒产量相关性状表现出较高的变异性；株高、主茎长、主穗粗、千粒重和出谷率的变异系数在11.55%~18.74%，生育期的变异系数低于10%。以上结果表明，120份谷子种质资源性状差异较大，具有丰富的遗传多样性。
表4 参试材料13个数量性状的变异分析及遗传多样性分析
Table 4
| 性状 Trait | 生育期 GP (d) | 株高 PH (cm) | 主茎长 MSL (cm) | 主穗长 MPL (cm) | 主穗直径 MPD (cm) | 主穗重 MPW (g) | 分蘖穗重 TPW (g) | 主穗粒重 GWMP (g) | 分蘖粒重 GWTP (g) | 千粒重 TGW (g) | 单株干草重 HWP (g) | 主穗 粒数 NSMP | 出谷率 MMP (%) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 最小值Minimum | 109.00 | 90.67 | 70.17 | 12.33 | 1.40 | 8.33 | 0.00 | 5.00 | 0.00 | 1.60 | 13.33 | 1960.78 | 50.00 |
| 最大值Maximum | 165.00 | 179.67 | 154.33 | 28.33 | 4.50 | 46.67 | 70.00 | 40.00 | 53.33 | 4.75 | 213.33 | 1 3725.49 | 91.00 |
| 极差Range | 56.00 | 89.00 | 84.17 | 16.00 | 3.10 | 38.33 | 70.00 | 35.00 | 53.33 | 3.15 | 200.00 | 1 1764.71 | 41.00 |
| 平均值Mean | 139.74 | 138.15 | 115.75 | 23.27 | 2.96 | 27.03 | 13.66 | 20.89 | 9.82 | 3.11 | 49.72 | 6754.67 | 76.00 |
| 标准差SD | 10.68 | 15.95 | 15.29 | 10.94 | 0.62 | 8.85 | 8.54 | 7.72 | 6.17 | 0.58 | 31.88 | 2303.66 | 9.00 |
| 变异系数CV (%) | 7.64 | 11.55 | 13.21 | 47.03 | 20.89 | 32.75 | 62.51 | 36.96 | 62.83 | 18.74 | 64.12 | 34.10 | 12.19 |
| 遗传多样性指数H′ | 1.78 | 1.74 | 1.83 | 0.63 | 1.79 | 2.02 | 1.41 | 2.02 | 1.43 | 1.80 | 1.45 | 1.87 | 1.97 |
GP：生育期，PH：株高，MSL：主茎长，MPL：主穗长，MPD：主穗直径，MPW：主穗重，TPW：分蘖穗重，GWMP：主穗粒重，GWTP：分蘖粒重，TGW：千粒重，HWP：单株干草重，NSMP：主穗粒数，MMP：出谷率，下同
GP: growth period, PH: plant height, MSL: main stem length, MPL: main panicle length, MPD: main panicle diameter, MPW: main panicle weight, TPW: tiller panicle weight, GWMP: grain weight of main panicle, GWTP: grain weight of tiller panicle, TGW: 1000-grain weight, HWP: hay weight per plant, NSMP: number of seeds in main panicle, MMP: milled millet percentage, the same below
2.2 谷子种质资源21个表型性状的相关性分析
120份种质资源21个表型性状相关性分析结果（表5）表明，大部分性状之间存在极显著（P<0.01）或显著（P<0.05）相关性，其中生育期与籽粒产量相关的因素均呈负相关，与单株干草重呈极显著正相关，这与谷子生长期间当地气候条件有关；主穗重与主穗粒重、主穗粒数及出谷率呈极显著正相关，与分蘖穗重、分蘖粒重及单株干草重呈极显著负相关；单株干草重与主穗重、主穗粒重及千粒重呈极显著负相关，与生育期、分蘖穗重及分蘖粒重呈极显著正相关。分析结果反映各表型性状间既相互独立又相互影响、相互制约的关系。
表5 谷子种质资源各性状的相关性
Table 5
| 性状 Trait | 生育期 GP | 幼苗 叶色 LC | 叶鞘色 LSC | 幼苗 叶姿 SL | 开花期 叶姿 SLF | 穗姿 SPN | 穗型 PT | 穗松 紧度 CS | 谷粒 色 HC | 株高 PH | 主茎 长 MSL | 主穗 长 MPL | 主穗 粗 MPD | 主穗 重 MPW | 分蘖 穗重 TPW | 主穗 粒重 GWMP | 分蘖 粒重 GWTP | 千粒 重 TGW | 单株 干草重 HWP | 主穗 粒数 NSMP | 出谷 率 MMP |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 生育期 GP | 1.00 | ||||||||||||||||||||
| 幼苗叶色 LC | 0.03 | 1.00 | |||||||||||||||||||
| 叶鞘色 LSC | -0.29** | 0.38** | 1.00 | ||||||||||||||||||
| 幼苗叶姿 SL | -0.01 | -0.11 | -0.01 | 1.00 | |||||||||||||||||
| 开花期 叶姿 SLF | -0.22* | 0.00 | -0.07 | 0.07 | 1.00 | ||||||||||||||||
| 穗姿SPN | -0.36** | 0.00 | 0.16 | -0.06 | 0.32** | 1.00 | |||||||||||||||
| 穗型PT | 0.05 | -0.04 | -0.15 | 0.15 | 0.09 | 0.14 | 1.00 | ||||||||||||||
| 穗松紧度 CS | -0.11 | -0.03 | -0.06 | -0.12 | 0.07 | -0.04 | -0.01 | 1.00 | |||||||||||||
| 谷粒色HC | 0.06 | -0.02 | 0.13 | 0.14 | 0.01 | 0.02 | -0.04 | -0.04 | 1.00 | ||||||||||||
| 株高PH | -0.03 | -0.03 | 0.04 | 0.17 | 0.05 | 0.20* | 0.08 | -0.11 | 0.02 | 1.00 | |||||||||||
| 主茎长 MSL | -0.03 | -0.04 | 0.05 | 0.16 | 0.03 | 0.12 | 0.09 | -0.04 | 0.06 | 0.96* | 1.00 | ||||||||||
| 主穗长 MPL | 0.06 | -0.02 | 0.03 | 0.05 | 0.04 | 0.11 | 0.14 | -0.08 | 0.04 | 0.22 | 0.07 | 1.00 | |||||||||
| 主穗粗 MPD | -0.03 | 0.05 | 0.05 | 0.08 | 0.08 | 0.19* | 0.24** | -0.10** | -0.02 | 0.03 | 0.01 | 0.10 | 1.00 | ||||||||
| 主穗重 MPW | -0.22* | -0.07 | -0.01 | -0.01 | 0.14 | 0.30** | 0.22* | 0.07* | -0.12 | 0.03 | -0.06 | 0.13 | 0.58** | 1.00 | |||||||
| 分蘖穗重 TPW | -0.06 | -0.05 | -0.04 | -0.03 | -0.12 | -0.02 | -0.10 | -0.03 | -0.04 | 0.02 | 0.09 | -0.07 | -0.35** | -0.46** | 1.00 | ||||||
| 主穗粒重 GWMP | -0.32** | -0.10 | -0.01 | -0.01 | 0.14 | 0.31** | 0.22** | 0.14* | -0.15 | 0.10 | 0.03 | 0.08 | 0.54** | 0.95** | -0.41** | 1.00 | |||||
| 分蘖粒重 GWTP | -0.05 | -0.07 | -0.01 | -0.04 | -0.12 | -0.04 | -0.13 | 0.02 | -0.07 | 0.03 | 0.11 | -0.08 | -0.39** | -0.50** | 0.98** | -0.43** | 1.00 | ||||
| 千粒重 TGW | -0.23* | -0.13 | -0.14 | -0.11 | 0.13 | 0.21* | 0.09* | 0.16 | -0.11 | 0.01 | -0.05 | 0.04 | 0.11 | 0.38** | -0.21* | 0.41** | -0.21* | 1.00 | |||
| 单株 干草重 HWP | 0.50** | -0.05 | -0.16 | 0.01 | -0.10 | -0.19* | -0.04* | 0.02 | 0.00 | 0.07 | 0.09 | 0.12 | -0.20* | -0.30** | 0.51** | -0.35** | 0.50** | -0.23* | 1.00 | ||
| 主穗粒数 NSMP | -0.23* | -0.03 | 0.09 | 0.06 | 0.08 | 0.22* | 0.17* | 0.05 | -0.09 | 0.11 | 0.07 | 0.06 | 0.54** | 0.82** | -0.33** | 0.86** | -0.35** | -0.09 | -0.26** | 1.00 | |
| 出谷率 MMP | -0.39** | -0.21* | -0.04 | -0.04 | 0.08 | 0.16 | 0.09 | 0.26** | -0.16 | 0.19* | 0.21* | -0.09 | 0.11 | 0.26** | 0.06 | 0.53** | 0.09 | 0.23* | -0.25** | 0.47** | 1.00 |
“*”代表在0.05水平上显著相关，“**”代表在0.01水平上显著相关
“*”represents significant correlation at P < 0.05,“**”represents significant correlation at P < 0.01
2.3 谷子种质资源农艺性状的主成分分析
对120份谷子种质资源的21个表型性状进行主成分分析（表6），遵循主成分累计贡献率≥80%的原则，前10个主成分累计贡献率达81.777%，包含了质量性状和数量性状的绝大部分信息，可作为谷子表型性状评价的综合指标。
表6 谷子种质资源主要农艺性状的主成分分析
Table 6
| 性状 Trait | 主成分Principal component | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
| 生育期GP | -0.210 | -0.235 | 0.332 | -0.233 | 0.000 | 0.000 | -0.134 | 0.000 | 0.000 | -0.161 |
| 幼苗叶色LC | 0.000 | -0.102 | 0.000 | 0.392 | -0.313 | -0.150 | -0.345 | -0.175 | -0.351 | 0.135 |
| 叶鞘色LSC | 0.000 | 0.000 | 0.000 | 0.527 | -0.337 | 0.000 | 0.000 | -0.104 | 0.120 | 0.233 |
| 幼苗叶姿SL | 0.000 | 0.000 | 0.278 | 0.000 | 0.150 | 0.125 | 0.561 | -0.132 | -0.181 | 0.287 |
| 开花期叶姿SLF | 0.120 | 0.000 | 0.000 | 0.000 | 0.314 | -0.395 | 0.103 | -0.379 | -0.330 | -0.464 |
| 穗姿SPN | 0.196 | 0.218 | 0.000 | 0.143 | 0.000 | -0.493 | 0.000 | 0.000 | 0.000 | -0.169 |
| 穗形PT | 0.122 | 0.000 | 0.225 | -0.201 | 0.000 | -0.260 | 0.131 | 0.000 | -0.403 | 0.542 |
| 穗松紧度CS | 0.000 | 0.000 | -0.240 | -0.217 | 0.000 | 0.145 | -0.271 | -0.729 | 0.000 | 0.199 |
| 谷粒色HC | 0.000 | 0.000 | 0.143 | 0.256 | 0.157 | 0.000 | 0.332 | -0.370 | 0.538 | 0.000 |
| 株高PH | 0.000 | 0.387 | 0.440 | 0.000 | 0.101 | 0.146 | -0.283 | 0.000 | 0.000 | 0.000 |
| 主茎长MSL | 0.000 | 0.413 | 0.402 | 0.000 | 0.000 | 0.231 | -0.243 | 0.000 | 0.000 | 0.000 |
| 主穗长MPL | 0.000 | 0.000 | 0.300 | 0.000 | 0.000 | -0.404 | -0.161 | 0.000 | 0.431 | 0.313 |
| 主穗粗MPD | 0.277 | -0.136 | 0.195 | 0.000 | -0.224 | -0.113 | 0.113 | 0.000 | 0.000 | 0.000 |
| 主穗重MPW | 0.403 | -0.103 | 0.000 | -0.141 | -0.150 | 0.000 | 0.000 | 0.000 | 0.130 | 0.000 |
| 分蘖穗重TPW | -0.271 | 0.332 | -0.118 | -0.196 | -0.271 | -0.211 | 0.150 | 0.000 | 0.000 | 0.000 |
| 主茎粒重GWMP | 0.423 | 0.000 | 0.000 | -0.183 | -0.158 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
| 分蘖粒重GWTP | -0.278 | 0.35 | -0.139 | -0.188 | -0.254 | -0.171 | 0.117 | 0.000 | 0.000 | 0.000 |
| 千粒重TGW | 0.189 | 0.000 | -0.206 | -0.146 | 0.403 | -0.151 | -0.282 | 0.190 | 0.160 | 0.244 |
| 单株干草重HWP | -0.262 | 0.000 | 0.212 | -0.306 | -0.196 | -0.220 | 0.000 | -0.180 | 0.125 | -0.144 |
| 主穗粒数NSMP | 0.360 | 0.000 | 0.000 | -0.101 | -0.400 | 0.105 | 0.150 | -0.110 | 0.000 | -0.160 |
| 出谷率MMP | 0.220 | 0.287 | -0.186 | -0.237 | -0.104 | 0.236 | 0.000 | -0.105 | 0.000 | 0.110 |
| 特征值Eigen value | 2.180 | 1.661 | 1.425 | 1.344 | 1.188 | 1.120 | 1.080 | 1.001 | 0.975 | 0.926 |
| 贡献率Contribution rate (%) | 21.600 | 12.544 | 9.227 | 8.217 | 6.411 | 5.701 | 5.306 | 4.552 | 4.322 | 3.897 |
| 累计贡献率Cumulative contribution rate (%) | 21.600 | 34.144 | 43.371 | 51.588 | 57.999 | 63.700 | 69.006 | 73.557 | 77.879 | 81.777 |
由表6可知，第1主成分贡献率为21.600%，主穗粒重、主穗重和主穗粒数特征向量绝对值较高，因此第1主成分为主穗产量因子；第2主成分贡献率为12.544%，分蘖粒重和分蘖穗重特征向量绝对值较高，因此称第2主成分为分蘖产量因子；第3主成分贡献率为9.227%，株高和主茎长特征向量绝对值较高，因此称第3主成分为株高构成因子；第4主成分贡献率为8.217%，叶鞘色、幼苗叶色和谷粒色特征向量绝对值较高，因此称第4主成分为叶鞘色因子；第5主成分贡献率为6.411%，千粒重特征向量绝对值较高，因此称第5主成分为千粒重因子；第6主成分贡献率为5.701%，穗姿、主穗长和开花期叶姿特征向量绝对值较高，因此称第6主成分为穗姿因子；第7主成分贡献率为5.306%，幼苗叶姿特征向量绝对值较高，因此称第7主成分为幼苗叶姿因子；第8主成分贡献率为4.552%，穗松紧度特征向量绝对值最高，因此称第8主成分为穗松紧度因子；第9主成分贡献率为4.322%，谷粒色特征向量绝对值较高，因此称第9主成分为谷粒色因子；第10主成分贡献率为3.897%，穗形特征向量绝对值较高，因此称第10主成分为穗形因子。
2.4 谷子种质资源21个表型性状的聚类及类群特征
对21个表型性状进行聚类分析（图2），将120份谷子种质资源划分为6大类，同时对6个类群谷子种质资源表型性状进行统计（表7）。类群I包括17份材料，其主要特征为株高、主茎长、主茎粒重和主穗粗最大，出谷率最高，为高秆、高产谷子优异种质，可作为选育高产品种的亲本材料。类群Ⅱ包括22份材料，其主要特征为生育期最短，分蘖穗重、分蘖粒重和单株干草重最小，属于早熟、少分蘖特异种质材料。类群Ⅲ包括21份材料，其主要特征为主穗长、主穗粗、主茎穗重和主茎粒重最小，分蘖穗重、分蘖粒重和单株干草重最大，为多穗、多分蘖特异种质材料，可作为饲草谷子选育亲本材料。类群Ⅳ包括16份材料，其主要特征为株高较矮，主茎长较小，主穗长较大，属于矮秆、长穗特异种质材料，可作为选育抗倒伏品种亲本。类群Ⅴ包括17份材料，其主要特征为株高较高、主穗重、千粒重和出谷率较高，属于具有增产潜力种质材料。类群Ⅵ包括27份材料，其主要特征为生育期最长，千粒重最大，属于晚熟、大粒特异种质材料。
图2
表7 6个谷子类群表型性状的平均值及特征
Table 7
| 性状 Trait | 种质群Germplasm group | |||||
|---|---|---|---|---|---|---|
| I | Ⅱ | Ⅲ | Ⅳ | Ⅴ | Ⅵ | |
| 幼苗叶色LC | 绿 | 绿，紫绿，黄绿 | 绿，黄绿 | 绿，紫绿，黄绿 | 绿 | 绿 |
| 叶鞘色LSC | 绿 | 绿，紫 | 绿，紫 | 绿，紫 | 绿，紫 | 绿，紫 |
| 幼苗叶姿SL | 平展，半上举 | 平展，半上举 | 半上举，平展 | 半上举，平展 | 半上举，平展 | 半上举，平展 |
| 开花期叶姿SLF | 半上举，上举 | 半上举，下披 | 下披，半上举 | 半上举，下披 | 半上举，下披 | 半上举，下披 |
| 穗姿SPN | 弯曲，中弯，勾形 | 弯曲，勾形 | 勾形，弯曲，中弯 | 勾形，弯曲，中弯 | 弯曲，勾形，直立 | 勾形，弯曲，中弯 |
| 穗形PT | 棍棒形，鸡嘴形 | 鸡嘴形，纺锤形 | 鸡嘴形，棍棒形 | 棍棒形，纺锤形 | 棍棒形，鸡嘴形 | 棍棒形，鸡嘴形 |
| 穗松紧度CS | 紧，中 | 紧，中 | 紧，中 | 中，紧 | 中，紧 | 中，松 |
| 谷粒色HC | 黄，青，红 | 黄，青，红，黑 | 黄，白，红，青 | 黄，橙，褐，黑 | 黄，青，褐 | 黄，黑，褐 |
| 生育期GP (d) | 137.53 | 135.33 | 141.64 | 141.81 | 136.65 | 144.24 |
| 株高PH (cm) | 145.39 | 140.86 | 136.65 | 133.16 | 136.31 | 139.45 |
| 主茎长MSL (cm) | 122.33 | 117.25 | 116.39 | 110.59 | 114.63 | 115.90 |
| 主穗长MPL (cm) | 22.78 | 24.04 | 20.16 | 27.34 | 21.18 | 22.48 |
| 主穗粗MPD (cm) | 3.42 | 3.21 | 2.43 | 2.90 | 3.13 | 2.82 |
| 主穗重MPW (g) | 36.67 | 33.71 | 14.91 | 27.10 | 30.10 | 22.06 |
| 分蘖穗重TPW (g) | 7.55 | 0.92 | 18.47 | 9.69 | 3.73 | 8.63 |
| 主茎粒重GWMP (g) | 29.86 | 27.37 | 10.30 | 19.69 | 24.12 | 16.69 |
| 分蘖粒重GWTP (g) | 4.90 | 0.83 | 14.00 | 6.66 | 2.55 | 6.41 |
| 千粒重TGW (g) | 2.94 | 3.12 | 3.07 | 3.06 | 3.19 | 3.34 |
| 单株干草重HWP (g) | 51.96 | 33.08 | 69.55 | 51.85 | 36.08 | 51.67 |
| 主穗粒数NSMP | 10 199.17 | 8790.41 | 3368.84 | 6440.34 | 7549.91 | 5000.83 |
| 出谷率MMP (%) | 82.00 | 82.00 | 69.00 | 74.00 | 80.00 | 76.00 |
2.5 优异谷子种质的筛选
综合考虑120份谷子种质资源各农艺性状表现，并从早熟、高产和单株干草重均高的3个方面考虑，筛选出11份优异谷子种质资源（表8）。
表8 优异谷子种质筛选
Table 8
| 特性 Trait | 编号 Number | 生育期 GP (d) | 株高 PH (cm) | 主茎长 MSL (cm) | 主穗长 MPL (cm) | 主穗粗 MPD (cm) | 主穗重 MPW (g) | 分蘖穗重 TPW (g) | 主穗粒重 GWMP (g) | 分蘖粒重 GWTP (g) | 千粒重 TGW (g) | 单株干草重 HWP (g) | 出谷率 MMP (%) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 早熟 Early ripen | G27 | 114 | 139.00 | 105.67 | 32.33 | 2.50 | 33.33 | 0.00 | 30.00 | 0.00 | 3.30 | 16.67 | 90 |
| G34 | 123 | 104.33 | 89.00 | 17.83 | 3.00 | 23.33 | 5.00 | 20.00 | 6.67 | 2.35 | 20.00 | 86 | |
| G67 | 124 | 137.33 | 117.00 | 21.67 | 3.00 | 25.00 | 15.00 | 18.33 | 11.67 | 3.75 | 33.33 | 73 | |
| G68 | 119 | 149.67 | 125.00 | 23.67 | 2.50 | 30.00 | 15.00 | 26.67 | 10.00 | 2.95 | 93.33 | 89 | |
| 高产 High yield | G4 | 129 | 146.67 | 125.00 | 20.67 | 3.67 | 16.67 | 50.00 | 13.33 | 26.67 | 3.20 | 73.33 | 80 |
| G23 | 155 | 141.00 | 126.00 | 17.67 | 1.67 | 13.33 | 70.00 | 10.00 | 53.33 | 2.70 | 213.33 | 75 | |
| G66 | 126 | 130.67 | 110.00 | 19.00 | 2.00 | 13.33 | 50.00 | 10.00 | 36.67 | 3.35 | 60.00 | 75 | |
| 干草重 HWP | G1 | 140 | 154.33 | 127.67 | 24.00 | 2.50 | 33.33 | 30.00 | 23.33 | 20.00 | 2.40 | 160.00 | 70 |
| G31 | 155 | 154.67 | 125.67 | 28.33 | 3.50 | 28.33 | 8.33 | 20.00 | 6.67 | 2.95 | 118.33 | 71 | |
| G106 | 147 | 135.67 | 121.00 | 21.33 | 2.33 | 20.00 | 66.67 | 16.67 | 43.33 | 2.80 | 100.00 | 83 | |
| G113 | 152 | 161.00 | 128.33 | 128.33 | 3.17 | 30.00 | 8.33 | 20.00 | 5.00 | 3.05 | 115.00 | 67 |
早熟材料4份，分别为78-0625（G27）、早谷（G34）、昆昌谷（G67）和黄沙突出谷（G68），生育期均小于125d（平均生育期120d）。其中，G68属于类群Ⅲ，株高中等（149.67cm），茎长中等（125.00cm），主穗长中等（23.67cm），主穗粗中等（2.50cm），单株穗重中等（45.00g），单株粒重中等（36.67g），千粒重中等（2.95g），单株草重较高（93.33g），出谷率较高（81%），可作为选育早熟材料的亲本。
高产材料3份，分别为338A（G4）、中牟一保谷（G23）和金德（G66），单株穗重均大于60.00g（平均穗重71.00g）。3份材料均属于类群Ⅱ。其中，G23株高中等（141.00cm），茎长中等（126.00cm），主穗长偏小（17.67cm），主穗粗偏小（1.67cm），单株穗重较高（83.33g），单株粒重较高（63.33g），千粒重偏小（2.70g），单株草重较高（213.33g），出谷率中等（75%），可作为选育高产材料的亲本。
单株干草重较高材料4份，分别为公谷68号（G1）、延谷13号（G31）、大洋谷（G106）和红酒谷（G113），单株干草重均大于100.00g（平均107.50g）。4份材料均属于类群Ⅵ。其中，G1株高中等（154.33cm），茎长中等（127.67cm），主穗长中等（24.00cm），主穗粗中等（2.50cm），单株穗重中等（60.33g），单株粒重中等（43.33g），千粒重偏小（2.40g），单株干草重较高（160.00g），出谷率中等（70%），可作为选育饲草材料的亲本。
3 讨论
种质资源是品种遗传改良的基础，形态学标记是鉴定种质资源遗传变异最简单直观的方法，在遗传多样性研究及育种材料的选择中有着重要作用[27]。王海岗等[11]对878份谷子核心种质的15个农艺性状进行形态标记分析，筛选出8个性状作为谷子种质评价指标，其中叶鞘色、株高、穗长、穗粒重和粒色在本研究中作为鉴定种质资源的主要指标。对21个性状进行调查，将主茎与分蘖分别标记，这与王晓娟等[28]以单株分析有所不同；遗传多样性结果表明，除穗形和主穗长多样性指数外，质量性状H′低于数量性状，这与前人在谷子[29]和苦荞[30]中的研究结果相一致；13个数量性状的遗传变异范围7.64%~64.12%，其中以单株干草重、分蘖穗重和分蘖粒重的变异系数最高，依次为64.12%、62.83%和62.51%；H′范围为0.63~2.02，其中以主穗重和主穗粒重H′最高，为2.02，其他性状也表现出较高的变异系数和H′，表明谷子种质遗传基础较广，存在丰富的遗传变异[31]；120份材料在同一生态种植条件下，不同性状在不同材料间表现出不同程度的多样性，充分说明种质资源间的遗传多样性[30⇓-32]，可以为谷子品种改良提供丰富的亲本材料。
在21个表型性状相关性分析中，数量性状与质量性状间相关系数较小，部分显著或极显著，而与数量性状间的相关系数均较高，且多数显著或极显著，说明各性状间紧密关联，遗传复杂性有待于在育种实践中验证；分蘖穗重与分蘖粒重、株高与主茎长之间的相关系数分别高达0.98和0.96，且达到极显著与显著水平，说明可以通过分蘖粒重和主茎长有效选择分蘖穗重和株高[33]，这将大幅度提高谷子选育效率；生育期与籽粒产量构成因素均呈负相关，与单株干草重呈极显著正相关。7-8月正值谷子抽穗灌浆期，降雨多以雷雨形式出现，水分利用率低，致使其灌浆期缩短[34]；9月阴雨天气较多，有效积温低，导致谷子贪青晚熟，延迟其刈割和成熟[35]，说明在甘肃干旱高寒区谷子的生育期过长是造成谷子籽粒产量减少的主要原因。
相关性分析表明，分蘖的穗重及粒重与主穗长、主穗粗及主穗重、主穗粒重呈显著负相关，这与赵会杰等[36]和佟汉文等[37]对小麦分蘖发育、分蘖对产量贡献研究结果一致，说明主茎对分蘖存在较强的竞争优势，主茎主要贡献于籽粒产量；而分蘖的穗重及粒重与单株干草重呈显著正相关关系，表明分蘖与植株的光吸收能力、营养物质的转运、分配密切相关，影响茎蘖间营养物质的运输与干物质积累[38-39]，分蘖主要贡献于植株生物产量。康晨等[40]对饲用高粱分蘖相关性状的分析研究表明，其生物产量与籽粒产量间存在负相关，分蘖增多可提高生物产量而降低籽粒产量，与本研究结果相似；而洛桑江白[41]对15个谷子品种的籽粒产量与单株草产量的相关分析与本研究有所不同。因此，关于谷子茎蘖间同化物的转运与分配、干物质积累未见明确报道，还需进一步探索。
本研究对120份谷子种质资源的21个表型性状进行主成分分析，10个主成分包含了性状的大部分信息，累计贡献率达到81.777%，且性状信息具有一定的相关性，其中主穗产量因子中主穗粒重和主穗粒数的特征向量大于其他农艺性状，吕丹等[42]在苦荞高产种质的筛选中发现株粒重和株粒数是苦荞籽粒产量的主要因素，与本研究结果一致。为了克服以个别性状进行分类的弊端，采用聚类分析将性状相近的种质分为6大类群，明确了谷子种质资源的不同类型，各类群在生育期、株高、分蘖性、籽粒产量和干草产量等性状间差异显著，有明显的形态学特征，根据育种目标选择亲本配制组合，为亲本的合理选配及选育不同用途谷子类型提供参考依据。
4 结论
本研究通过对120份谷子种质资源农艺性状的鉴定分析，明确饲草型谷子是甘肃干旱高寒地区谷子育种的主要方向。同时，基于当地气候条件，根据熟性、籽粒产量和干草产量等表型筛选出78-0625、早谷、昆昌谷和黄沙突出谷等早熟材料4份，338A、中牟一保谷和金德等高产材料3份，公谷68号、延谷13号、大洋谷和红酒谷等干草产量较高材料4份，从而为谷子品种改良和亲本选择提供理论依据。
参考文献
谷子(Setaria italica (L.) P. Beauv.)作为功能基因组研究模式植物的发展现状及趋势
谷子分蘖相关QTL定位及紧密连锁标记开发
DOI:10.7668/hbnxb.2018.06.005
[本文引用: 1]
分蘖是与农作物产量相关的重要农艺性状之一，分蘖对于谷子产量的提高具有重要意义，但是目前谷子分蘖遗传分子机制还不清楚。利用简化基因组测序技术（RAD-seq），以2个F<sub>2</sub>群体（命名为Cross AJ和Cross HC）为研究对象，分别对2个F<sub>2</sub>群体的543个和131个单株及其亲本进行RAD-seq，利用MSTmap和WinQTLCart 2.5软件进行遗传图谱构建和QTL分析，结果共鉴定出8个控制分蘖相关的QTL。其中，利用Cross AJ共鉴定了6个QTL，分别为qAJTN1、qAJTN5、qAJTN7-1、qAJTN7-2、qAJTN7-3和qAJTN9，可解释表型变异0.7%~9.8%；利用Cross HC群体共鉴定了2个QTL，分别为qHCTN5和qHCTN7，可解释表型变异1.4%~8.3%。在这8个QTL中，qAJTN1、qAJTN5、qAJTN7-2、qAJTN9和qHCTN5为该研究新鉴定的位点，其余3个QTL（qAJTN7-1、qAJTN7-3和qHCTN7）与前人研究结果相一致。另外，通过将亲本的测序结果与参考基因组进行比对，搜索插入缺失位点InDel（Insertion-deletion），新开发了谷子分蘖相关QTL（qAJTN7-3和qHCTN7）紧密相关的InDel标记。结果为谷子分蘖机制的研究奠定一定基础，所开发的分子标记对于分蘖型谷子的分子标记辅助育种具有重要意义。
西北半干旱区引种燕麦品种产量与品质的关联分析及评价
DOI:10.11733/j.issn.1007-0435.2018.01.015
[本文引用: 1]
为了明确不同引种燕麦品种资源的生产特性和利用潜力，应用相关和灰色关联度对22个燕麦品种的主要农艺性状和营养品质进行了关联分析和综合评价。结果表明：不同燕麦品种对同一生态环境的适应性差异较大，营养品质与农艺性状间存在显著的（P<0.05）的相关性。生育期与籽粒产量的关联度最大（0.8350），所占权重最高（0.1146），分蘖数与籽粒产量的关联度最低（0.6042），所占权重最小（0.0829），11个性状的权重大小顺序为：生育期 > 株高 > 籽粒产量 > 粗蛋白 > 淀粉 > 小穗数 > 千粒重 > 粗脂肪 > 穗粒重 > 穗粒数 > 分蘖数。通过加权关联度最优评价得出综合性状优良的皮燕麦为‘定燕2号’、‘冀张燕4号’、‘坝燕4号’，‘白燕14号’与‘冀张燕5号’加权关联度最小，裸燕麦中‘远杂2号’、‘白燕15号’和‘燕科2号’的加权关联度较大，表现突出。结果可为甘肃中部乃至西北黄土高原半干旱地区燕麦种质创新和引种推广提供科学依据。
不同生态区谷子种质资源表型比较分析
DOI:10.13304/j.nykjdb.2020.0175
[本文引用: 1]
为明确不同生态区谷子种质资源的表型，从春播特早熟区、春播中熟区、春播晚熟区引进435个谷子品种，在赤峰地区进行表型鉴定及比较分析。结果表明，3个生态区谷子品种的生育期、农艺性状、产量性状差异显著。春播特早熟区谷子品种生育期平均值为（98.35±13.00）d，分蘖多、植株矮、主茎粗、旗叶短、主穗小、穗粒重、出谷率高、产量中等，表型性状遗传多样性整体较低；春播中熟区谷子品种生育期平均值为（11307±8.63）d，分蘖少、植株高、穗下节长、主茎细、旗叶短、主穗大、穗粒重、出谷率高、产量较低，产量性状遗传多样性整体较高；春播晚熟区谷子品种生育期平均值为（125.80±6.76） d；穗下节短、主茎粗、旗叶长、主穗大、穗粒轻、出谷率低、产量较高，农艺性状遗传多样性整体较高。利用标准化处理后的16个表型进行聚类分析，可将参试品种划分为4大类群，类群Ⅰ品种中熟、生物量大、出谷率和产量较高，类群Ⅱ品种晚熟、生物量小、出谷率和产量较低，类群Ⅲ品种早熟、生物量小、产量较高，类群Ⅳ品种中晚熟、生物量大、产量较高。
The mathematical theory of communication
A method for the effi-cient management and utilization of large germplasm collections
DOI:10.2135/cropsci1999.0011183X003900040046x URL [本文引用: 1]
谷子分蘖相关QTL定位及紧密连锁标记开发
DOI:10.7668/hbnxb.2018.06.005
[本文引用: 1]
分蘖是与农作物产量相关的重要农艺性状之一，分蘖对于谷子产量的提高具有重要意义，但是目前谷子分蘖遗传分子机制还不清楚。利用简化基因组测序技术（RAD-seq），以2个F<sub>2</sub>群体（命名为Cross AJ和Cross HC）为研究对象，分别对2个F<sub>2</sub>群体的543个和131个单株及其亲本进行RAD-seq，利用MSTmap和WinQTLCart 2.5软件进行遗传图谱构建和QTL分析，结果共鉴定出8个控制分蘖相关的QTL。其中，利用Cross AJ共鉴定了6个QTL，分别为qAJTN1、qAJTN5、qAJTN7-1、qAJTN7-2、qAJTN7-3和qAJTN9，可解释表型变异0.7%~9.8%；利用Cross HC群体共鉴定了2个QTL，分别为qHCTN5和qHCTN7，可解释表型变异1.4%~8.3%。在这8个QTL中，qAJTN1、qAJTN5、qAJTN7-2、qAJTN9和qHCTN5为该研究新鉴定的位点，其余3个QTL（qAJTN7-1、qAJTN7-3和qHCTN7）与前人研究结果相一致。另外，通过将亲本的测序结果与参考基因组进行比对，搜索插入缺失位点InDel（Insertion-deletion），新开发了谷子分蘖相关QTL（qAJTN7-3和qHCTN7）紧密相关的InDel标记。结果为谷子分蘖机制的研究奠定一定基础，所开发的分子标记对于分蘖型谷子的分子标记辅助育种具有重要意义。
饲用高粱分蘖相关性状的遗传分析
