王海岗, 温琪汾, 穆志新^,, 乔治军^,山西省农业科学院农作物品种资源研究所/农业部黄土高原作物基因资源与种质创制重点实验室/杂粮种质资源发掘与遗传改良山西省重点实验室,太原 030031

Population Structure and Association Analysis of Main Agronomic Traits of Shanxi Core Collection in Foxtail Millet

WANG HaiGang, WEN QiFen, MU ZhiXin^,, QIAO ZhiJun^,Institute of Crop Germplasm Resources of Shanxi Academy of Agricultural Sciences/Key Laboratory of Crop Gene Resources and Germplasm Enhancement on Loess Plateau, Ministry of Agriculture/Shanxi Key Laboratory of Genetic Resources and Genetic Improvement of Minor Crops, Taiyuan 030031

通讯作者: 穆志新,E-mail：muzx2008@sina.com乔治军,E-mail：nkypzs@126.com

责任编辑: 李莉
收稿日期:2019-04-22接受日期:2019-06-20网络出版日期:2019-11-16

基金资助:

山西省平台基地和人才专项优秀人才科技创新项目.201705D211026 山西省农业科学院特色农业技术攻关.YGG17057 国家谷子高粱产业技术体系.CARS-06-13.5-A16

Received:2019-04-22Accepted:2019-06-20Online:2019-11-16
作者简介 About authors
王海岗,E-mail：nkywhg@126.com






摘要
【目的】分析山西谷子地方品种遗传多样性和群体遗传结构,筛选与谷子农艺性状相关联的分子标记,为谷子杂交组合亲本选配及分子标记辅助育种提供依据。【方法】 利用96对SSR标记对595份山西谷子核心资源进行全基因组扫描,采用PowerMarker 3.25软件分析群体遗传多样性,利用STRUCTURE 2.3.4软件分析群体遗传结构,使用TASSEL 2.1软件中GLM（general linear model,Q）和MLM（mixed linear model,Q+K）2种方法,进行表型和标记关联分析。【结果】 96对SSR引物共扩增出828个等位变异,平均每对引物扩增到8.6个,变化范围为2—26;基因多样性指数变化范围为0.005—0.941,平均为0.610;多态信息量变化范围为0.005—0.938,平均为0.577;各位点杂合度变化范围为0—0.050,平均位点杂合度仅为0.016。群体结构分析将595份核心资源分为3个亚群。4 560个SSR位点成对组合中,共线性组合和非共线性组合之间都存在一定的连锁不平衡。D′统计概率（P<0.01）支持的LD成对位点1 955个,占全部位点组合的42.9%,D′平均值为0.23。通过GLM方法共检测到12个极显著性位点（P<0.01）,表型变异解释率为2.34%—13.94%,平均为6.33%,贡献率较高的等位变异位点是CAAS2050（R²=13.94%）和B153（R²=11.36%）;通过MLM方法共检测到9个极显著性位点（P<0.01）,表型变异解释率为2.80%—9.22%,平均为5.16%,贡献率较高的等位变异位点是P89（R²=9.22%）和P3*（R²=8.28%）;2种方法共同检测到的极显著性位点有7个。 【结论】 利用SSR标记分析了595份山西谷子核心资源的遗传多样性和群体遗传结构。2种关联分析模型中,GLM方法关联到12个标记与节数、株高、颈长、茎粗、穗长、穗粗、码数、码粒数、蛋白质含量9个性状相关;MLM方法关联到9个标记与节数、颈长、叶宽、茎粗、穗粗、码数、码粒数、千粒重8个性状相关。
关键词： 谷子;山西;地方品种;SSR;遗传多样性;连锁不平衡;关联分析

Abstract
【Objective】The objective of this study is to detect the SSR markers associated with agronomic trait and analyze genetic diversity and genetic structure of foxtail millet landrace in Shanxi province. The results will be helpful for hybridization combination of parent materials and molecular marker assisted breeding.【Method】96 SSR markers on 9 chromosomes were genome-wide screened for polymorphism in core collection of 595 accessions. PowerMarker 3.25 software was used to estimate the polymorphism information of population. Population structure was analyzed using STRUCTURE 2.3.4 software. Then the data were associated with 96 SSR markers by GLM (general linear model, Q) and MLM (mixed linear model, Q+K).【Result】Totally 828 alleles were found with 96 SSR markers and 8.6 alleles were revealed with each marker in average ranged from 2-26. The gene diversity was from 0.005 to 0.941, averagely 0.610. The polymorphism information content (PIC) value ranged from 0.005 to 0.938 with the mean of 0.577. Heterozygosity per locus on average was 0.016, ranging from 0 to 0.050. All the 595 accessions were divided into three subgroups by analysis of population genetic structure. There was linkage disequilibrium (LD) among linked loci and unlinked loci pairs, and 1 955 out of 4 560 loci pairs (42.9%) had significant LD (P < 0.01) with average D′ value of 0.23. A total of 12 locus found by GLM method significantly at the level of P<0.01 which explained 2.34%-13.94% of the phenotypic variance and the mean value was 6.33%. CAAS2050 (R²=13.94%) and B153(R²=11.36%) kept the max value. Meanwhile, 9 loci were found by MLM method significantly at the level of P<0.01 which explained 2.80%-9.22% of the phenotypic variance and the mean value was 5.16%. P89(R²=9.22%) and P3*(R²=8.28%) kept the max value. A total of 7 loci were detected in common by GLM and MLM. 【Conclusion】Genetic diversity and population structure of 595 accessions were analyzed through SSR markers. In the two association analysis models, 12 markers were associated with nine traits including stem node number, plant height, peduncle length, diameter of main stem, panicle length, panicle diameter, primary branch number per panicle, spikelet number per primary branch, protein content by GLM. Nine markers were associated with eight traits including stem node number, peduncle length, leaf width, diameter of main stem, panicle diameter, primary branch number per panicle, spikelet number per primary branch, 1000-grain weight by MLM.
Keywords：foxtail millet;Shanxi;landrace;SSR;genetic diversity;linkage disequilibrium;association analysis


PDF (2476KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
王海岗, 温琪汾, 穆志新, 乔治军. 山西谷子核心资源群体结构及主要农艺性状关联分析[J]. 中国农业科学, 2019, 52(22): 4088-4099 doi:10.3864/j.issn.0578-1752.2019.22.013
WANG HaiGang, WEN QiFen, MU ZhiXin, QIAO ZhiJun. Population Structure and Association Analysis of Main Agronomic Traits of Shanxi Core Collection in Foxtail Millet[J]. Scientia Acricultura Sinica, 2019, 52(22): 4088-4099 doi:10.3864/j.issn.0578-1752.2019.22.013

0 引言

【研究意义】谷子（Setaria italica（L.）Beauv.）是起源于中国黄河流域的粮饲兼用作物,又称粟,去壳后为小米,在中国的栽培历史可追溯到11 500年前^[1],是世界范围内最古老的禾谷类作物之一。山西谷子种质资源丰富,占全国谷子种质资源的1/5^[2],山西传统的名优谷子品种较多,如东方亮、沁州黄、晋谷21号、晋谷10号,这些品种蕴藏着丰富的优异基因资源。优异基因发掘是作物种质资源研究的重要部分,对作物分子育种具有重要的实践意义。关联分析是基因发掘的有效途径之一,以连锁不平衡为基础,通过对自然群体表型数据和基因型数据的统计检测,寻找性状可遗传变异和标记多态性之间的关联^[3,4,5]。【前人研究进展】目前,关联分析已在水稻库源相关性状^[6]、纹枯病抗性^[7]等;小麦株高^[8]、旗叶叶绿素^[9]、穗发芽抗性^[10]等性状;大豆农艺性状^[9]、加工和品质性状^[12,13,14]、幼苗期耐淹性^[15]及棉花农艺和纤维品质性状^[16]、叶绿素含量^[17]、适宜机采相关性状^[18]等已进行了大量研究。且在大麦^[19,20,21]、谷子^{[22,23,24,25]}、小豆^[28]、豇豆^[29]等小宗作物也进行了一定研究。山西谷子年播种面积在22.46×10⁴ hm²左右,占全国的1/4^[30],但是作为山西传统杂粮作物的谷子分子辅助育种相对缓慢,充分挖掘山西本土谷子资源中的优良等位变异,利用分子标记辅助选择,创制谷子新种质是谷子资源研究急需开展的工作。【本研究切入点】本研究中所用材料主要来自20世纪50—80年代征集的山西谷子地方品种^[31,32],这些材料未经过育种的选择压力,多数资源是通过自然选择或突变后筛选而保留下来的种质,是环境和历史的载体材料。传统育种主要根据作物的表型和产量性状进行杂交亲本组合配制,对遗传组成的差异考虑较少,导致亲本来源单一,育成品种间遗传基础狭窄,很难培育出具有突破性的优良品种。【拟解决的关键问题】本研究对595份山西谷子地方品种初选核心资源进行遗传多样性分析,明确山西谷子地方品种的遗传结构,寻找与农艺性状显著相关联的分子标记,为谷子种质创新中的优异等位变异发掘、亲本组配和标记辅助选择等提供参考。

1 材料与方法

1.1 试验材料及表型鉴定

所用595份山西谷子地方品种来源于5 627份山西谷子资源构建的初选核心资源^[30]及近年新收集的谷子资源。2017—2018年在山西省农业科学院东阳试验基地种植,对抽穗期、节数、株高、颈长、叶长、叶宽、穗长、穗粗、茎粗、单穗重、穗粒重、码数、码粒数、千粒重、蛋白质含量等15个性状进行鉴定。

1.2 基因组DNA提取

谷子抽穗后取倒二叶,液氮速冻,-80℃保存备用。采用植物基因组DNA提取试剂盒（DP305-03）提取DNA,利用全功能酶标仪（美国伯腾SynergyH1）测定DNA浓度,-20℃保存备用。

1.3 SSR分析

选取分布于9条染色体上的96对引物进行多态性扫描（表1）^{[22,23,24,25,26,27]}。SSR反应体系包括模板DNA 1.0 μL、2×Tap PCR Mastermix 5.0 μL、10 mmol·L^-1 Forward primer 0.4 μL、10 mmol·L^-1 Reverse primer 0.4 μL和ddH₂O 3.2 μL。PCR反应程序为95℃ 5 min;95℃ 30 s,58℃（不同引物退火温度根据温度梯度试验确定）30 s,72℃ 60 s,32个循环;72℃ 10 min;4℃保存。扩增产物用8%聚丙烯酰胺凝胶电泳分离,银染显色。

Table 1
表1
表196对SSR引物遗传多样性
Table 1Genetic diversity at 96 SSR markers

标记 Marker	染色体 Chromosome	基因型数 Genotype No.	等位基因数 Allele No.	基因多样性指数 Gene diversity	多态性信息含量 PIC	杂合度 Heterozygosity
In1-1	Chr.1	3	2	0.4955	0.3727	0.0183
In1-3	Chr.1	5	4	0.1810	0.1663	0.0200
In1-7	Chr.1	7	4	0.3955	0.3521	0.0267
In1-9	Chr.1	3	2	0.1113	0.1051	0.0083
CAAS1054	Chr.1	34	19	0.9051	0.8976	0.0283
P3*	Chr.1	27	15	0.8264	0.8072	0.0233
P88	Chr.1	27	14	0.8852	0.8756	0.0233
B153	Chr.1	13	10	0.7789	0.7477	0.0050
P58	Chr.1	14	8	0.7565	0.7184	0.0333
B218	Chr.1	18	9	0.8303	0.8084	0.0250
CAAS2006	Chr.2	11	5	0.7162	0.6675	0.0167
CAAS2010	Chr.2	13	6	0.6777	0.6269	0.0250
In2138	Chr.2	3	2	0.4867	0.3682	0.0267
In2184	Chr.2	4	4	0.1995	0.1815	0.0000
b242	Chr.2	10	9	0.7997	0.7730	0.0017
CAAS2050	Chr.2	11	6	0.6971	0.6435	0.0300
CAAS2030	Chr.2	6	4	0.5026	0.4165	0.0283
In2-11	Chr.2	3	2	0.0535	0.0521	0.0017
B169	Chr.2	23	18	0.9235	0.9182	0.0083
B249	Chr.2	24	22	0.9286	0.9241	0.0033
MPGD13	Chr.2	8	5	0.5607	0.4671	0.0067
b101	Chr.3	22	11	0.8607	0.8449	0.0233
B163	Chr.3	10	7	0.7193	0.6760	0.0067
B186	Chr.3	37	17	0.9182	0.9123	0.0383
P61*	Chr.3	6	4	0.6897	0.6299	0.0100
B225	Chr.3	26	16	0.8963	0.8869	0.0233
P98	Chr.3	9	6	0.5330	0.4808	0.0267
MPGD32	Chr.3	9	6	0.4217	0.3865	0.0150
MPGD44	Chr.3	8	6	0.4247	0.3914	0.0100
B224	Chr.3	47	25	0.9411	0.9379	0.0383
P78	Chr.3	14	9	0.7941	0.7657	0.0150
In4-3	Chr.4	5	4	0.4948	0.3796	0.0317
CAAS4019	Chr.4	7	5	0.4388	0.4096	0.0083
CAAS4033	Chr.4	14	8	0.8181	0.7936	0.0200
IN4-6	Chr.4	3	2	0.2608	0.2268	0.0117
IN4-7	Chr.4	6	4	0.2004	0.1909	0.0117
B109	Chr.4	38	26	0.8654	0.8533	0.0267
P2	Chr.4	28	18	0.9045	0.8970	0.0183
B247*	Chr.4	22	14	0.8866	0.8762	0.0167
P100**	Chr.4	8	7	0.5368	0.4883	0.0017
P89	Chr.4	24	17	0.8428	0.8246	0.0267
CAAS5048	Chr5	9	6	0.7318	0.6894	0.0050
In5-7	Chr.5	2	2	0.0099	0.0099	0.0000
In5-8	Chr.5	3	2	0.4382	0.3422	0.0283
CAAS5032	Chr.5	26	12	0.8911	0.8810	0.0300
b111	Chr.5	17	12	0.8842	0.8729	0.0100
P17X*	Chr.5	12	6	0.5827	0.5436	0.0200
B223	Chr.5	32	19	0.9141	0.9078	0.0317
B237	Chr.5	19	8	0.8512	0.8329	0.0267
B117	Chr.5	11	7	0.5355	0.4950	0.0233
MPGA51	Chr.5	5	4	0.2576	0.2367	0.0017
MPGD4	Chr.5	4	3	0.3705	0.3400	0.0117
In6-2	Chr.6	3	2	0.2733	0.2360	0.0133
CAAS6023	Chr.6	7	5	0.2109	0.1972	0.0067
In6-9	Chr.6	3	2	0.4986	0.3743	0.0400
In6-11	Chr.6	4	3	0.3183	0.2697	0.0017
CAAS6018	Chr.6	10	8	0.7108	0.6608	0.0117
CAAS6007	Chr.6	17	10	0.7759	0.7517	0.0183
B159	Chr.6	23	11	0.8445	0.8263	0.0233
P32	Chr.6	8	6	0.3231	0.2957	0.0050
P12X	Chr.6	18	12	0.8195	0.8004	0.0183
B250	Chr.6	35	20	0.8997	0.8917	0.0267
MPGA50	Chr.6	8	5	0.6601	0.6232	0.0150
In7-1	Chr.7	3	2	0.0050	0.0050	0.0017
CAAS7021	Chr.7	7	5	0.4033	0.3782	0.0050
SIMS14450	Chr.7	6	5	0.2846	0.2612	0.0117
CAAS7036	Chr.7	3	2	0.1569	0.1446	0.0150
In7-13	Chr.7	3	2	0.4911	0.3705	0.0500
P45	Chr.7	5	5	0.5788	0.4933	0.0000
B142	Chr.7	16	12	0.8185	0.7957	0.0100
B123	Chr.7	14	11	0.6491	0.6170	0.0083
P29	Chr.7	21	13	0.8830	0.8713	0.0150
B200	Chr.7	9	5	0.5637	0.5071	0.0083
B180	Chr.7	10	7	0.7220	0.6758	0.0067
SI 227	Chr.7	12	10	0.8304	0.8084	0.0050
In8-1	Chr.8	7	4	0.6511	0.5845	0.0183
P6	Chr.8	5	5	0.3672	0.3143	0.0000
CAAS8015	Chr.8	10	7	0.7853	0.7523	0.0100
CAAS8011	Chr.8	19	8	0.8140	0.7888	0.0250
B185	Chr.8	32	22	0.9234	0.9181	0.0183
B258	Chr.8	39	21	0.9211	0.9157	0.0300
J2（GTDZ7-1）	Chr.8	2	2	0.0099	0.0099	0.0000
P14	Chr.8	25	W	0.7974	0.7791	0.0150
B128	Chr.8	25	22	0.9309	0.9267	0.0050
b246	Chr.9	23	12	0.8797	0.8676	0.0333
CAAS9036	Chr.9	24	15	0.8811	0.8703	0.0183
b166	Chr.9	20	9	0.8320	0.8107	0.0267
CAAS9081	Chr.9	17	8	0.8058	0.7790	0.0250
CAAS9082	Chr.9	12	8	0.7926	0.7633	0.0083
CAAS9058	Chr.9	7	6	0.1620	0.1555	0.0100
P20*	Chr.9	15	11	0.7462	0.7074	0.0100
P41	Chr.9	11	9	0.7927	0.7623	0.0033
SIMS 298	Chr.9	4	3	0.4286	0.3423	0.0150
B145	Chr.9	9	6	0.5104	0.4788	0.0133
B212	Chr.9	3	3	0.0394	0.0390	0.0000
SI 110	Chr.9	5	5	0.1118	0.1093	0.0000
平均值 Mean		13	8.6	0.6097	0.5773	0.0159

新窗口打开|下载CSV

1.4 数据分析

根据SSR扩增产物的电泳结果,参照DNA Ladder记录扩增产物条带的分子量大小。使用Powermarker 3.25软件计算每个位点等位基因数、基因型数、基因多样性指数、多态性信息含量指数和杂合度,并计算Nei遗传距离采用MEGA 5.0软件构建NJ聚类图。

采用Structure2.3.4软件分析群体遗传结构,K取值为1—10,3次重复;将MCMC（markov chain monto carlo）开始时的不作数迭代设为100 000次,再将不作数迭代后的MCMC设为100 000次,其余参数采用软件默认的设置。根据lnP(D)计算ΔK,并依据ΔK值选择一个合适的K值。使用Tassel 2.1软件计算SSR位点组合间的连锁不平衡（linkage disequilibrium,LD）水平及支持概率,绘制LD配对检测矩阵图,使用一般线性模型（general linear model,GLM）和混合线性模型（mixed liner model,MLM）结合基因型数据、表型值和Q值,对SSR标记与表型性状进行关联分析,计算P<0.01时标记位点对表型变异的贡献率（R²）。采用STRUCTURE软件计算Q值,用软件SPAGeDi1.5软件计算Kinship值。

2 结果

2.1 SSR多样性

96对SSR引物在595份山西谷子地方品种核心资源中共检测到828个等位变异,变化范围2—26个,平均为8.6个,其中,标记B109检测到的等位变异最多为26个（表1）。基因型数变化范围在2—47,平均每个位点为13。基因多样性指数变化范围为0.0050—0.9411,平均为0.6097;多态信息量变化范围为0.0050—0.9379,平均为0.5773;各位点杂合度变化范围为0—0.050,平均位点杂合度仅为0.0159,说明谷子异交率很小接近于纯系。每条染色体上标记平均数有10对,平均等位变异为92个;第3和第8染色体检测到等位变异最多,均为107个,第7染色体最少为79个。第3染色体基因多样性指数和多态性信息含量最高,分别为0.720和0.6912,第7染色体最低,分为0.5322和0.4940。第7染色体杂合度最小,为0.0114,第1染色体杂合度最大,为0.0212。

2.2 群体结构分析

采用Structure 2.3.4软件对595份山西谷子核心资源进行群体遗传结构分析,计算每份材料的Q值（每份种质被分配到对应亚群中的概率）。随着K值的增大,LnP(D)递增,未出现拐点,因此无法准确判断K值（图1）。参照EVANNO等^[33]方法,当K=3时,ΔK出现最大值,因此该群体材料被划分为3个亚群（图2）。按照亚群中个体的概率≥0.6时,个体被分配到相应聚类的亚群;亚群中个体的概率≤0.6时,个体被分配到一个混合类群的分配原则,分析每份核心资源对应的Q值。结果表明,595份材料归属为3个亚群,分别包含77、222和187份材料,3个亚群各占比例为12.9%、37.3%和31.4%;另有109份材料在任意一个亚群内Q值均<0.6,单独划为一个混合群,占比例为18.3%。

图1

新窗口打开|下载原图ZIP|生成PPT
图1595份材料的基于群体结构分析的K值与lnP(D)值和ΔK值变化图

A：ln P(D)值与K值的变化图;B：ΔK值随K值的变化图
Fig. 1lnP(D) and ΔK based on population structure analysis for 595 foxtail millet landrace

A: Line chart of K and ln P(D); B: Magnitude of ΔK as a function of K

图2

新窗口打开|下载原图ZIP|生成PPT
图2595份山西谷子地方品种的群体遗传结构

Fig. 2Population structure of 595 foxtail millet landrace



基于SSR遗传距离的11个地市聚类分析可以看出（图3）,运城、晋城、临汾、长治、阳泉聚为第一类;大同、朔州、忻州、太原聚为第二类;晋中、吕梁聚为第3类。划分的3个类群按照纬度高低自北向南依次分开,这与山西复杂多样的地理生态类型有一定关系,说明山西谷子种质资源存在地理生态型的分化。

图3

新窗口打开|下载原图ZIP|生成PPT
图3基于SSR标记数据的595份核心资源不同地市间遗传距离聚类图

Fig. 3Dendrogram of inter-city group of the 595 core collection based on SSR marker

2.3 连锁不平衡分析

关联分析以不同位点等位基因间的连锁不平衡为基础,D′描述群体的重组史。对分布于谷子全基因组96对SSR标记进行连锁不平衡分析,4 560个SSR位点成对组合中,共线性组合（同一染色体）和非共线性组合（不同染色体）之间都存在一定的连锁不平衡（图4）。D′统计概率（P<0.01）支持的LD成对位点1 955个,占全部位点组合的42.9%,D′平均值为0.23。LD成对位点的D′值范围主要集中在0—0.2和0.2—0.4（表2）,在第2和第4染色体上未发现D′在0.5以上的LD位点,说明这两条染色体较其他7条染色体受到的选择压力小;分析共线性位点组合的LD,发现第4、第6和第9染色体上的LD成对位点分别有31、35和36个,高于其他染色体,且第6染色体上D′平均值最高为0.28,说明群体材料中,该染色体上发生重组选择要大于其他染色体。

图4

新窗口打开|下载原图ZIP|生成PPT
图49条染色体SSR位点的连锁不平衡分布

Fig. 4Distribution of LD among SSR locus on 9 chromosomes



Table 2
表2
表2不同D'值范围LD成对位点数
Table 2Number of LD pairs in different D' value ranges

D' 值范围 D' value ranges	LD成对点数 Number of LD pairs
0—0.2	946
0.2—0.4	870
0.4—0.6	103
0.6—0.8	29
0.8—1	7

新窗口打开|下载CSV

2.4 SSR标记与性状关联分析

采用Tassel 2.1软件分析了595份山西谷子核心资源在一点2年2种模型的15个性状表型和基因型,结果显示,在P<0.01下,共检测到14个SSR标记与11个性状相关联（表3）。GLM方法检测到12个标记与节数、株高、颈长、茎粗、穗长、穗粗、码数、码粒数和蛋白质9个性状相关联。表型变异解释率为2.34%—13.94%,平均值为6.33%。表型变异解释率较大的标记是CAAS2050和B153,解释率分别为13.94%和11.36%;MLM方法检测到9个标记与节数、颈长、叶宽、茎粗、穗粗、码数、码粒数和千粒重8个性状相关联。表型变异解释率为2.80%—9.22%,平均为5.16%。贡献率较高的位点是P89和P3*,解释率分别为9.22%和8.28%。比较GLM和MLM结果,CAAS6023、P3*、In6-2、CAAS2050和B117 5个标记在2种模型条件下同时被检测到。B153、P58、B117、In6-2、CAAS6023和B142标记同时与2个或多个性状相关联,这可能与性状相关有一定关系,说明控制该性状的等位基因可能连锁或一因多效。

Table3
表3
表3GLM和MLM关联分析
Table3Association analysis by the GLM and MLM (P<0.01)

性状 Trait	位点 SSR loci	染色体 Chromosome	GLM						MLM
			2017			2018			2017			2018
			R2	P		R2	P		R2	P		R2	P
节数SN	P58	Chr.1	0.0572	6.42E-08		0.0597	1.73E-07
	B142	Chr.7										0.052	6.39E-05
株高PH	P58	Chr.1	0.0541	1.06E-07		0.0547	2.19E-07
颈长PL	CAAS6023	Chr.6	0.0558	3E-09		0.0517	1.77E-08					0.0298	8.8E-05
叶宽LW	In6-2	Chr.6										0.0444	4.34E-05
茎粗DMS	P3*	Chr.1	0.0774	1.04E-05		0.0977	4.95E-08					0.0828	3.73E-06
	In6-2	Chr.6	0.0458	2.91E-05		0.0498	8.95E-06					0.0445	5.1E-05
	P58	Chr.1							0.0522	3.73E-05
	P89	Chr.4							0.0922	1.3E-06
	B142	Chr.7										0.0628	6.19E-05
穗长PL	B153	Chr.1	0.0605	7.03E-07		0.0517	7.02E-06
	CAAS7036	Chr.7	0.0234	3.08E-05		0.0306	1.16E-06
穗粗PD	B142	Chr.7	0.068	1.37E-05		0.0646	5.05E-05
	B117	Chr.5							0.0457	6.4E-05
码数PBNP	CAAS2050	Chr.2	0.0514	3.48E-06		0.1394	3.59E-18					0.053	1.61E-06
	In2-11	Chr.2										0.028	1.82E-05
	In6-2	Chr.6	0.0258	5.04E-05		0.0277	3.55E-05
	B142	Chr.7	0.0583	9.36E-05		0.0772	1.14E-06
码粒数SNPB	B117	Chr.5	0.0559	4.73E-06		0.0439	9.4E-05		0.046	6.78E-05
千粒重GW	CAAS6023	Chr.6							0.0375	5.05E-05
蛋白含量PC	B153	Chr.1	0.1136	1.12E-12		0.1081	4.84E-11
	In1-7	Chr.1	0.0512	3.08E-07		0.0425	1.09E-05
	P14	Chr.8	0.0901	3.77E-07		0.1108	7.03E-09

SN: Stem node number; PH: Plant height; PL: Peduncle length; LW: Leaf width; DMS: Diameter of main stem; PL: Panicle length; PD: Panicle diameter; PBNP: Primary branch number per panicle; SNPB: Spikelet number per primary branch; GW: 1000-grain weight; PC: Protein content; GLM: General linear model; MLM: Mixed linear model
SN：节数;PH：株高;PL：颈长;LW：叶宽;DMS：茎粗;PL：穗长;PD：穗粗;PBNP：码数;SNPB：码粒数;GW：千粒重;PC：蛋白含量;GLM：一般线性模型;MLM：混合线性模型

新窗口打开|下载CSV

与茎粗相关的标记有5个,表型变异解释率为4.45%—9.77%,平均值为6.72%,贡献率最高的位点为P89;与码数相关的标记有4个,表型变异解释率为2.58%—13.94%,平均值为5.76%,贡献率最高的位点为CAAS2050;与蛋白质相关的位点有3个,表型变异解释率为4.25%—11.36%,平均值为8.61%,贡献率最高的位点为B153;与节数相关的位点有2个,表型变异解释率为5.20%—5.97%,平均值为5.63%,贡献率最高的位点为P58;与穗长相关的位点有2个,表型变异解释率为2.34%—6.05%,平均值为4.16%,贡献率最高的位点为B153;与穗粗相关的位点有2个,表型变异解释率为4.57%—6.46%,平均值为5.94%,贡献率最高的位点为B142。

3 讨论

3.1 山西谷子资源的遗传多样性

优异基因发掘是种质资源高效利用和作物遗传改良的基础,种质资源的遗传多样性评价和群体遗传结构的分析是优异基因发掘的前提条件。关联分析将蕴藏在不同种质中的等位变异与表型性状建立了联系,通过寻找显著关联位点,可以发掘种质资源中优异等位基因。山西谷子资源丰富,栽培历史悠久,名优品种众多,但是传统谷子育种只是简单的品种杂交,通过应用分子标记辅助选择聚合有利基因,改良谷子地方品种,是山西谷子育种和种质创制的新途径。本研究采用SSR标记对山西谷子地方品种核心资源进行遗传多样性评价,分析群体遗传结构,挖掘与农艺性状关联的分子标记,对指导山西谷子资源的利用和育种亲本选择具有指导意义。

通过对595份山西谷子核心资源进行分子标记扫描,发现平均等位基因8.6个,平均基因多样性指数0.610,平均PIC为0.577,低于WANG等^[22]和JIA等^[23]对全国谷子地方品种和育成品种研究结果,同时也低于VETRIVENTHAN等^[24]的研究,说明山西谷子资源的遗传多样性要低于全国和世界的谷子资源。但与GUPTA等^[25]和JIA等^[34]研究结果相近,这与等位变异统计方法和试验选择引物有一定关系。山西行政区域从北纬34°34′到40°44′横跨6个纬度,海拔从180 m到3 058 m,地貌类型复杂多样,全省各县均有谷子资源分布,形成了多类型谷子种植区,按照纬度和海拔分为春播早熟区、春播中熟区、春播晚熟区和夏播区^[35]。山西谷子核心资源按照11个地市来源划分成3个类群,基本按照纬度高低由北到南划分开来^[36],说明山西谷子种质资源存在地理生态型分化,也从分子水平证明了山西省是春、夏谷2种类型同时存在的省份。同时,明确了各地市谷子资源的遗传距离,为山西谷子生态区域划分和育种亲本选配提供了分子证据。

3.2 SSR位点连锁不平衡及关联分析

已有研究表明^[22,24],谷子LD衰减距离在15—40 cM,说明谷子的LD水平较高,有利于进行全基因组关联分析,进而发掘种质资源中的有利基因。本研究利用覆盖全基因组的96对SSR引物进行标记扫描,4 560个位点成对组合中,不论是共线性组合还是非共线性组合,都有LD存在。D' 统计概率（P<0.01）支持的LD成对位点占42.9%,较高水平的LD成对位点（D' >0.5）主要分布在第1、第2、第6和第9染色体上及与其组合的位点,说明这4条染色体面临的选择压力大于其他染色体,这些材料多经历了较多农艺性状的选择。本研究发现核心资源的整体连锁不平衡水平较低（D′ 平均值为0.23）,说明山西谷子资源在经历人为重组和自然突变的历程较少,且谷子是自花授粉作物,自然异交率低,地方品种中携带的大量优异等位变异仍单独存在,所以山西谷子育种仍有进一步提升空间。关联分析在极显著（P<0.01）条件下共检测到14个与农艺性状和蛋白质含量相关的位点,对表型变异解释率平均值为5.98%,在谷子育种中可利用这些标记进行材料的等位变异标记辅助选择研究;另外,有6个标记位点同时与2个以上性状相关联,其中B142位点与4个性状相关联,这可能是性状相关和一因多效的遗传基础,说明这些标记可以同时进行改良多个目标性状,可用于多个性状的聚合育种研究。

4 结论

595份谷子核心资源遗传结构存在3个类群,11个地市按照纬度由高到低划分成3类,山西谷子资源存在地理生态型分化。共线性或非共线性位点组合都有LD存在,核心资源群体有利于关联分析和优异等位变异挖掘。2种关联模型下分别找到12个与节数、株高、颈长、茎粗、穗长、穗粗、码数、码粒数、蛋白质含量相关联的标记和9个与节数、颈长、叶宽、茎粗、穗粗、码数、码粒数、千粒重相关联的标记。

参考文献 View Option 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1] 贾冠清, 刁现民 . 谷子(Setaria italica (L.) P. Beauv.)作为功能基因组研究模式植物的发展现状及趋势
生命科学, 2017,29(3):292-301.
[本文引用: 1]

JIA G Q, DIAO X M . Current status and perspectives of researches on foxtail millet (Setaria italica(L.) P. Beauv.): A potential model of plant functional genomics studies.
Chinese Bulletin of Life Sciences, 2017,29(3):292-301. (in Chinese)
[本文引用: 1]

[2] 王海岗, 秦慧彬, 田翔, 吕建珍, 陈凌, 王君杰, 曹晓宁, 刘思辰, 王纶, 温淇汾, 穆志新, 乔治军 . 山西谷子地方品种表型多样性分析
中国农学通报, 2018,34(32):19-24.
[本文引用: 1]

WANG H G, QIN H B, TIAN X, Lü J Z, CHEN L, WANG J J, CAO X N, LIU S C, WANG L, WEN Q F, MU Z X, QIAO Z J . Phenotypic diversity of local foxtail millet varieties in Shanxi
Chinese Agricultural Science Bulletin, 2018,34(32):19-24. (in Chinese)
[本文引用: 1]

[3] 王荣焕, 王天宇, 黎裕 . 关联分析在作物种质资源分子评价中的应用
植物遗传资源学报, 2007,8(3):366-372.
[本文引用: 1]

WANG R H, WANG T Y, LI Y . Application of association analysis in molecular evaluation of crop germplasm resources
Journal of Plant Genetic Resources, 2007,8(3):366-372. (in Chinese)
[本文引用: 1]

[4] 杨小红, 严建兵, 郑艳萍, 余建明, 李建生 . 植物数量性状关联分析研究进展
作物学报, 2007,33(4):523-530.
[本文引用: 1]

YANG X H, YAN J B, ZHENG Y P, YU J M, LI J S . Reviews of association analysis for quantitative traits in plants
Acta Agronomica Sinica, 2007,33(4):523-530. (in Chinese)
[本文引用: 1]

[5] 冯建英, 温阳俊, 张瑾, 章元明 . 植物关联分析方法的研究进展
作物学报, 2016,42(7):945-956.
[本文引用: 1]

FENG J Y, WEN Y J, ZHANG J, ZHANG Y M . Advances on methodologies for genome-wide association studies in plants
Acta Agronomica Sinica, 2016,42(7):945-956. (in Chinese)
[本文引用: 1]

[6] 赵宏亮, 陈凯, 张强, 徐正进, 徐建龙, 黎志康 . 基于连锁不平衡水稻源库相关性状的关联分析
核农学报, 2015,29(4) : 674-684.
[本文引用: 1]

ZHAO H L, CHEN K, ZHANG Q, XU Z J, XU J L, LI Z K . Association analysis on source sink-related traits in rice based on linkage disequilibrium
Journal of Nuclear Agricultural Sciences, 2015,29(4):674-684. (in Chinese)
[本文引用: 1]

[7] 孙晓棠, 卢冬冬, 欧阳林娟, 胡丽芳, 边建民, 彭小松, 陈小荣, 傅军如, 贺晓鹏, 贺浩华, 朱昌兰 . 水稻纹枯病抗性关联分析及抗性等位变异发掘
作物学报, 2014,40(5):779-787.
[本文引用: 1]

SUN X T, LU D D, OUYANG L J, HU L F, BIAN J M, PENG X S, CHEN X R, FU J R, HE X P, HE H H, ZHU C L . Association mapping and resistant alleles analysis for sheath blight resistance in rice
Acta Agronomica Sinica, 2014,40(5):779-787. (in Chinese)
[本文引用: 1]

[8] 魏添梅, 昌小平, 闵东红, 景蕊莲 . 小麦抗旱品种的遗传多样性分析及株高优异等位变异挖掘
作物学报, 2010,36(6):895-904.
[本文引用: 1]

WEI T M, CHANG X P, MIN D H, JING R L . Analysis of genetic diversity and tapping elite alleles for plant height in drought-tolerant wheat varieties
Acta Agronomica Sinica, 2010,36(6):895-904. (in Chinese).
[本文引用: 1]

[9] 张冬玲, 张洪娜, 郝晨阳, 王兰芬, 李甜, 张学勇 . 花后冠层温度对小麦产量的影响及几个关联SSR位点的效应分析
作物学报, 2015,41(4):548-556.
[本文引用: 2]

ZHANG D L, ZHANG H N, HAO C Y, WANG L F, LI T, ZHANG X Y . Influence of canopy temperature (ct) during grain-filling period on yield and effects of several ct-associated SSR loci
Acta Agronomica Sinica, 2015,41(4):548-556. (in Chinese)
[本文引用: 2]

[10] 朱玉磊, 王升星, 赵良侠, 张德新, 胡建帮, 曹雪连, 杨亚杰, 常成, 马传喜, 张海萍 . 以关联分析发掘小麦整穗发芽抗性基因分子标记
作物学报, 2014,40(10):1725-1732.
[本文引用: 1]

ZHU Y L, WANG S X, ZHAO L X, ZHANG D X, HU J B, CAO X L, YANG Y J, CHANG C, MA C X, ZHANG H P . Exploring molecular markers of preharvest sprouting resistance gene using wheat intact spikes by association analysis
Acta Agronomica Sinica, 2014,40(10):1725-1732. (in Chinese)
[本文引用: 1]

[11] 杨胜先, 牛远, 李梦, 魏世平, 刘晓芬, 吕海燕, 章元明 . 栽培大豆农艺性状的关联分析及优异等位变异挖掘
中国农业科学, 2014,47(20):3941-3952.


YANG S X, NIU Y, LI M, WEI S P, LIU X F, Lü H Y, ZHANG Y M . Association mapping of agronomic traits in soybean ( Glycine max L. Merr.) and mining of novel alleles
Scientia Agricultura Sinica, 2014,47(20):3941-3952. (in Chinese)


[12] 文自翔, 赵团结, 郑永战, 刘顺湖, 王春娥, 王芳, 盖钧镒 . 中国栽培和野生大豆农艺品质性状与SSR标记的关联分析: I.群体结构及关联标记
作物学报, 2008,34(7):1169-1178.
[本文引用: 1]

WEN Z X, ZHAO T J, ZHENG Y Z, LIU S H, WANG C E, WANG F, GAI J Y . Association analysis of agronomic and quality traits with SSR markers in glycine max and glycine soja in China: I. Population structure and associated markers
Acta Agronomica Sinica, 2008,34(7):1169-1178. (in Chinese)
[本文引用: 1]

[13] 文自翔, 赵团结, 郑永战, 刘顺湖, 王春娥, 王芳, 盖钧镒 . 中国栽培和野生大豆农艺及品质性状与SSR标记的关联分析: Ⅱ.优异等位变异的发掘
作物学报, 2008,34(8):1339-1349.
[本文引用: 1]

WEN Z X, ZHAO T J, ZHENG Y Z, LIU S H, WANG C E, WANG F, GAI J Y . Association analysis of agronomic and quality traits with SSR markers in glycine max and glycine soja in china: ii. exploration of elite alleles
Acta Agronomica Sinica, 2008,34(8):1339-1349. (in Chinese)
[本文引用: 1]

[14] 范虎, 文自翔, 王春娥, 王芳, 邢光南, 赵团结, 盖钧镒 . 中国野生大豆群体农艺加工性状与SSR关联分析和特异材料的遗传构成
作物学报, 2013,39(5):775-788.
[本文引用: 1]

FAN H, WEN Z X, WANG C E, WANG F, XING G N, ZHAO T J, GAI J Y . Association analysis between agronomic-processing traits and SSR markers and genetic dissection of specific accessions in Chinese wild soybean population
Acta Agronomica Sinica, 2013,39(5):775-788. (in Chinese)
[本文引用: 1]

[15] 孙慧敏, 张军, 赵团结, 盖钧镒 . 亚洲地区中、外大豆品种幼苗期耐淹性与SSR标记的关联分析
作物学报, 2010,36(10):1615-1623.
[本文引用: 1]

SUN H M, ZHANG J, ZHAO T J, GAI J Y . Association analysis between submergence tolerance and SSR markers in domestic and foreign soybean cultivars in Asia
Acta Agronomica Sinica, 2010,36(10):1615-1623. (in Chinese)
[本文引用: 1]

[16] 聂新辉, 尤春源, 鲍健, 李晓方, 惠慧, 刘洪亮, 秦江鸿, 林忠旭 . 基于关联分析的新陆早棉花品种农艺和纤维品质性状优异等位基因挖掘
中国农业科学, 2015,48(15):2891-2910.
[本文引用: 1]

NIE X H, YOU C Y, BAO J, LI X F, HUI H, LIU H L, QIN J H, LIN Z X . Exploration of elite alleles of agronomic and fiber quality traits in Xinluzao cotton varieties by association analysis
Scientia Agricultura Sinica, 2015,48(15):2891-2910. (in Chinese)
[本文引用: 1]

[17] 刘其宝, 李黎贝, 张驰, 宿俊吉, 魏恒玲, 王寒涛, 喻树迅 . 陆地棉叶片叶绿素含量与SSR标记的关联分析及优异等位变异的挖掘
中国农业科学, 2017,50(18):3439-3449.
[本文引用: 1]

LIU Q B, LI L B, ZHANG C, SU J J, WEI H L, WANG H T, YU S X . Association analysis of leaf chlorophyll content with SSR markers and exploration of superior alleles in upland cotton
Scientia Agricultura Sinica, 2017,50(18):3439-3449. (in Chinese)
[本文引用: 1]

[18] 王娟, 董承光, 刘丽, 孔宪辉, 王旭文, 余渝 . 棉花适宜机采相关性状的SSR标记关联分析及优异等位基因挖掘
作物学报, 2010,36(10):1615-1623.
[本文引用: 1]

WANG J, DONG C G, LIU L, KONG X H, WANG X W, YU Y . Association analysis and exploration of elite alleles of mechanical harvest-related traits with SSR markers in upland cotton cultivars ( Gossypium hirsutum L.)
Acta Agronomica Sinica, 2010,36(10):1615-1623. (in Chinese)
[本文引用: 1]

[19] 赖勇, 王鹏喜, 范贵强, 司二静, 王晋, 杨轲, 孟亚雄, 李葆春, 马小乐, 尚勋武, 王化俊 . 大麦SSR标记遗传多样性及其与农艺性状关联分析
中国农业科学, 2013,46(2):233-242.
[本文引用: 1]

LAI Y, WANG P X, FAN G Q, SI E J, WANG J, YANG K, MENG Y X, LI B C, MA X L, SHANG X W, WANG H J . Genetic diversity and association analysis using SSR markers in barley
Scientia Agricultura Sinica, 2013,46(2):233-242. (in Chinese)
[本文引用: 1]

[20] 姜晓东, 郭刚刚, 张京 . 中国大麦地方品种的遗传多样性及α-淀粉酶活性的全基因组关联分析
中国农业科学, 2013,46(4):668-677.
[本文引用: 1]

JIANG X D, GUO G G, ZHANG J . Analysis on the genetic diversity in Chinese barley landrace and genomic wide association of α-amylase activity genes
Scientia Agricultura Sinica, 2013,46(4):668-677. (in Chinese)
[本文引用: 1]

[21] 赖勇, 王鹏喜, 范贵强, 司二静, 王晋, 杨轲, 孟亚雄, 李葆春, 马小乐, 尚勋武, 王化俊 . 大麦遗传多样性及连锁不平衡分析
作物学报, 2013,39(12):2154-2161.
[本文引用: 1]

LAI Y, WANG P X, FAN G Q, SI E J, WANG J, YANG K, MENG Y X, LI B C, MA X L, SHANG X W, WANG H J . Genetic diversity and linkage disequilibrium analysis in barley
Acta Agronomica Sinica, 2013,39(12):2154-2161. (in Chinese)
[本文引用: 1]

[22] WANG C F, JIA G Q, ZHI H, NIU Z G, CHAI Y, LI W, WANG Y F, LI H Q, LU P, ZHAO B H, DIAO X M . Genetic diversity and population structure of Chinese foxtail millet [Setaria italica (L.) Beauv.] landraces.
Genes Genomes Genetics, 2012,7(2):769-777.
[本文引用: 4]

[23] JIA G Q, LIU X T, SCHNABLE J C, NIU Z G, WANG C F, LI Y H, WANG S J, WANG S Y, LIU J R, GUO E H, ZHI H, DIAO X M . Microsatellite variations of elite Setaria varieties released during last six decades in China
PLoS ONE, 10(5):e0125688.
[本文引用: 3]

[24] VETRIVENTHAN M, UPADHYAYA HARI D, ANANDAKUMAR C R, SENTHILVEL S, VARSHNEY R K, PARZIES H K . Population structure and linkage disequilibrium of ICRISAT foxtail millet (Setaria italica (L.) P. Beauv.) core collection.
Euphytica, 2014,196:423-435
[本文引用: 4]

[25] GUPTA S, KUMARI K, MUTHAMILARASAN M, PARIDA W K, PRASAD M . Population structure and association mapping of yield contributing agronomic traits in foxtail millet
Plant Cell Report, 2014,33:881-893
[本文引用: 3]

[26] 王春芳 . 利用微卫星标记分析中国谷子地方品种的群体结构与遗传多样性
[D]. 石家庄: 河北师范大学, 2011.
[本文引用: 1]

WANG C F . Population structure and genetic diversity of foxtail millet landraces in China assessed by microsatellites
[D]. Shijiazhuang: Hebei Normal University, 2011. (in Chinese)
[本文引用: 1]

[27] ZHANG S, TANG C J, ZHAO Q, LI J, YANG L F, QIE L F, FAN X K, LI L, ZHANG N, ZHAO M C, LIU X T, CHAI Y, ZHANG X, WANG H L, LI Y T, LI W, ZHI H, JIA G Q, DIAO X M . Development of highly polymorphic simple sequence repeat markers using genome- wide microsatellite variant analysis in foxtail millet [Setaria italica (L.) P. Beauv.].
BMC Genomics, 2014,15:78.
[本文引用: 1]

[28] 白鹏, 程须珍, 王丽侠, 王素华, 陈红霖 . 小豆遗传差异、群体结构和连锁不平衡水平的SSR分析
作物学报, 2014,40(5):788-797.
[本文引用: 1]

BAI P, CHENG X Z, WANG L X, WANG S H, CHEN H L . Genetic diversity, population structure and linkage disequilibrium in adzuki bean by using SSR markers
Acta Agronomica Sinica, 2014,40(5):788-797. (in Chinese)
[本文引用: 1]

[29] 熊海铮 . 豇豆遗传多样性及若干农艺性状关联分析
[D]. 杭州: 浙江大学, 2016.
[本文引用: 1]

XIONG H Z . Genetic diversity and association analysis of agronomic traits in Cowpea
[D]. Hangzhou: Zhejiang University, 2016. (in Chinese)
[本文引用: 1]

[30] 王海岗, 温琪汾, 乔治军, 穆志新 . 山西谷子地方品种初选核心种质构建
农学学报, 2019,9(4):26-31.
[本文引用: 2]

WANG H G, WEN Q F, QIAO Z J, MU Z X . Construction of candidate core collection local varieties of foxtail millet in Shanxi
Journal of Agriculture, 2019,9(4):26-31. (in Chinese)
[本文引用: 2]

[31] 温琪汾, 王纶, 赵卫红, 畅建武 . 山西省谷子种质资源的繁种入库
山西农业科学, 2002 , 30(4):32-34.
[本文引用: 1]

WEN Q F, WANG L, ZHAO W H, CHANG J W . Reproduction and storage of foxtail millet germplasm in Shanxi
Journal of Shanxi Agricultural Sciences, 2002,30(4):32-34. (in Chinese)
[本文引用: 1]

[32] 温琪汾, 王纶, 王星玉 . 山西省谷子种质资源及抗旱种质的筛选利用
山西农业科学, 2005,33(4):32-33.
[本文引用: 1]

WEN Q F, WANG L, WANG X Y . The foxtail millet germplasm resources and screening and utilization of drought resistance germplasm in Shanxi
Agricultural Sciences, 2005,33(4):32-33. (in Chinese)
[本文引用: 1]

[33] EVANNO G, REGNAUTS S, GOUDET J . Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study
Molecular Ecology, 2005,14:2611-2620.
[本文引用: 1]

[34] JIA X P, ZHANG Z B, LIU Y H, ZHANG C W, SHI Y S, SONG Y C, WANG T Y, LI Y . Development and genetic mapping of SSR markers in foxtail millet [Setaria italica (L.) P. Beauv.].
Theoretical and Applied Genetics, 2009,118(4):821-829.
[本文引用: 1]

[35] 古兆明, 古世禄 . 山西谷子起源与发展研究. 北京: 中国农业科学技术出版社, 2007: 166-168.
[本文引用: 1]

GU Z M, GU S L. Origin and Development of Shanxi Foxtail Millet. Beijing: Agricultural Science and Technology Press, 2007: 166-168. (in Chinese)
[本文引用: 1]

[36] 任继海, 温琪汾, 王星玉, 王纶 . 山西省谷子品种类型及其分布
华北农学报, 1999,14(3):25-30.
[本文引用: 1]

REN J H, WEN Q F, WANG X Y, WANG L . Varieties of foxtail millet and their distribution in Shanxi province
Acta Agriculturae Boreali-Sinica, 1999,14(3):25-30. (in Chinese)
[本文引用: 1]