利用SNP标记对100份玉米种质资源的遗传多样性分析
Genetic Diversity Analysis of 100 Maize Germplasm Resources by SNP Markers
通讯作者:
收稿日期: 2024-01-16 修回日期: 2024-03-22
基金资助: |
|
Received: 2024-01-16 Revised: 2024-03-22
作者简介 About authors
孙艳杰,主要从事玉米遗传育种研究,E-mail:
杂种优势是玉米育种的重要理论基础,杂种优势类群和杂种优势模式的深入研究为玉米品种选育和杂交种的组配方式提供了理论指导。为了提高育种效率,使用60K SNP芯片对100份玉米种质资源进行基因型分析。通过对遗传距离进行聚类分析,将100份玉米自交系划分为3个类群,分别为Reid、Non-Reid和Dom群。遗传相似度的主成分分析结果与遗传距离的聚类分析结果一致。
关键词:
Heterosis is a key theoretical basis of maize breeding, and through research on heterosis groups and heterosis patterns provide theoretical guidance for the selection of maize varieties and combination methods of hybrids. In order to improve breeding efficiency, the 100 maize germplasm resources were genotyped using a 60K SNP chip. By cluster analysis of genetic distance, the 100 maize inbred lines were divided into three groups, including Reid, Non-Reid and Dom. The results of principal component analysis of genetic similarity were consistent with those of cluster analysis of genetic distance.
Keywords:
本文引用格式
孙艳杰, 魏国才, 吴雨恒, 石运强, 邵勇, 刘英蕊, 南元涛, 张维耀.
Sun Yanjie, Wei Guocai, Wu Yuheng, Shi Yunqiang, Shao Yong, Liu Yingrui, Nan Yuantao, Zhang Weiyao.
玉米杂种优势类群的划分方法有系谱分析、表型聚类分析、配合力分析、同工酶法以及分子标记法。现在主要用分子标记法进行种质资源类群的划分,分子标记法包括SSR标记、SNP标记及芯片技术[4]。作为第三代分子标记的SNP标记,其优点为遗传稳定性好、分布密度高和可自动化分析等,因此SNP标记常被应用于玉米杂种优势类群划分与遗传多样性分析[5]。卢媛等[6]用34 257个SNP标记把44份糯玉米自交系分为5个类群,分别为P群、塘四平头群、旅大红骨群、Reid群和Lancaster群。高嵩等[7]利用6973个SNP位点将205份玉米自交系划分为7个类群,分别为Lancaster群、塘四平头群、旅大红骨群、Reid、PA、PB群和热带类群。姜思奇等[8]利用46 899个高质量的SNP标记将36份玉米自交系分为4个杂种优势群,分别为Lancaster群、Reid群、塘四平头与旅大红骨混合群和类PH4CV群。师亚琴等[9]利用4550个SNP位点将80份玉米自交系划分为7个类群,包括Reid、黄改群、P群、Lancaster、Iodent、迪卡选系和先锋改良群。Li等[10]和Wang等[11]对不同玉米杂种优势类群的基因组和表型分化的遗传基础进行了研究,从分子层面对玉米优势类群的机理进行解析,这对利用分子标记划分类群提供了理论支撑。
杂种优势类群划分结果的种类数量差异是由类群划分的标准不同造成的,类群划分结果均可在育种工作中作为参考。杂种优势类群划分的种类过多会导致杂种模式过于复杂,对育种工作会有所干扰而降低育种效率。目前在我国北方的生产应用中,育种家主要利用Reid、Non-Reid(以Lancaster为主)和Dom群(塘四平头群和旅大红骨类群)这三大优势类群的材料,对具有国外血缘的自交系可划分为Reid和Non-Reid两大杂种优势群,Dom群为国内种质,包括塘四平头群和旅大红骨类群[12]。本试验利用玉米60K SNP芯片鉴定100份玉米种质资源并进行遗传多样性分析,为品种选育和提高育种效率提供一定的借鉴。
1 材料与方法
1.1 试验材料
供试材料为黑龙江省农业科学院绥化分院提供的100份玉米资源,编号、名称及来源详见表1。其中以6个已知来源的骨干自交系作为测验种,其中SH060、SH075和SH083属于Reid群,SH085属于Non-Reid群,SH084和SH097属于Dom群。
表1 100份玉米种质资源
Table 1
编号Number | 名称Name | 来源Source | ||
---|---|---|---|---|
SH001 | KWS10/KWS73/ 欧早 | 黑龙江省农业科学院绥化分院 | ||
SH002 | SX618/欧早 | 黑龙江省农业科学院绥化分院 | ||
SH003 | S1029 | 黑龙江省农业科学院绥化分院 | ||
SH004 | SI15412 | 黑龙江省农业科学院绥化分院 | ||
SH005 | S1035 | 黑龙江省农业科学院绥化分院 | ||
SH006 | S1037 | 黑龙江省农业科学院绥化分院 | ||
SH007 | S1040 | 黑龙江省农业科学院绥化分院 | ||
SH008 | SI15111 | 黑龙江省农业科学院绥化分院 | ||
SH009 | SI15215 | 黑龙江省农业科学院绥化分院 | ||
SH010 | LK702 | 黑龙江省龙科种业集团有限公司 | ||
SH011 | LK708 | 黑龙江省龙科种业集团有限公司 | ||
SH012 | LK725 | 黑龙江省龙科种业集团有限公司 | ||
SH013 | LK727 | 黑龙江省龙科种业集团有限公司 | ||
SH014 | LK729 | 黑龙江省龙科种业集团有限公司 | ||
SH015 | LK738 | 黑龙江省龙科种业集团有限公司 | ||
SH016 | SD3♀/欧早 | 黑龙江省农业科学院绥化分院 | ||
SH017 | N3527 | 黑龙江省农业科学院绥化分院 | ||
SH018 | N3528 | 黑龙江省农业科学院绥化分院 | ||
SH019 | N3529-2 | 黑龙江省农业科学院绥化分院 | ||
SH020 | N3532-1 | 黑龙江省农业科学院绥化分院 | ||
SH021 | N3533-4 | 黑龙江省农业科学院绥化分院 | ||
SH022 | N3535-1 | 黑龙江省农业科学院绥化分院 | ||
SH023 | N3538-3 | 黑龙江省农业科学院绥化分院 | ||
SH024 | N3540-2 | 黑龙江省农业科学院绥化分院 | ||
SH025 | N3542-3 | 黑龙江省农业科学院绥化分院 | ||
SH026 | N3543-2 | 黑龙江省农业科学院绥化分院 | ||
SH027 | N3545-4 | 黑龙江省农业科学院绥化分院 | ||
SH028 | N3546-3 | 黑龙江省农业科学院绥化分院 | ||
SH029 | N3547-2 | 黑龙江省农业科学院绥化分院 | ||
SH030 | N3547-4 | 黑龙江省农业科学院绥化分院 | ||
SH031 | N3528-2 | 黑龙江省农业科学院绥化分院 | ||
SH032 | N3540-3 | 黑龙江省农业科学院绥化分院 | ||
SH033 | N3542-2 | 黑龙江省农业科学院绥化分院 | ||
SH034 | N3543-2 | 黑龙江省农业科学院绥化分院 | ||
SH035 | SYN12 | 黑龙江省农业科学院绥化分院 | ||
SH036 | K3/郑58 | 黑龙江省农业科学院绥化分院 | ||
SH037 | B73/郑58 | 黑龙江省农业科学院绥化分院 | ||
SH038 | 30715-3 | 黑龙江省农业科学院绥化分院 | ||
SH039 | 30716-4 | 黑龙江省农业科学院绥化分院 | ||
SH040 | K3/341 | 黑龙江省农业科学院绥化分院 | ||
SH041 | 8941/340 | 黑龙江省农业科学院绥化分院 | ||
SH042 | 344/9229 | 黑龙江省农业科学院绥化分院 | ||
SH043 | SDM7 | 黑龙江省农业科学院绥化分院 | ||
SH044 | 加群 | 中国农业科学院作物科学研究所 | ||
SH045 | 98-4/郑58 | 黑龙江省农业科学院绥化分院 | ||
SH046 | 4207-5/C8605 | 黑龙江省农业科学院绥化分院 | ||
SH047 | 9176/郑58 | 黑龙江省农业科学院绥化分院 | ||
SH048 | 扎461/(330/ M67) | 黑龙江省农业科学院绥化分院 | ||
SH049 | 30716 | 黑龙江省农业科学院绥化分院 | ||
SH050 | 15D86 | 吉林省宏宇种业有限公司 | ||
SH051 | 15D90 | 吉林省宏宇种业有限公司 | ||
SH052 | W08 | 吉林省宏宇种业有限公司 | ||
SH053 | 93622 | 黑龙江省农业科学院绥化分院 | ||
SH054 | 15D98 | 吉林省宏宇种业有限公司 | ||
SH055 | 15D109 | 吉林省宏宇种业有限公司 | ||
SH056 | 6171/SX711 | 黑龙江省农业科学院绥化分院 | ||
SH057 | 98-4/扎461 | 黑龙江省农业科学院绥化分院 | ||
SH058 | 344/SX711 | 黑龙江省农业科学院绥化分院 | ||
SH059 | SX711/SX709/ 81162 | 黑龙江省农业科学院绥化分院 | ||
SH060(Reid) | SX718 | 黑龙江省农业科学院绥化分院 | ||
SH061 | SX718/8941/ C7-2 | 黑龙江省农业科学院绥化分院 | ||
SH062 | RA1 | 中国农业科学院作物科学研究所 | ||
SH063 | RA3 | 中国农业科学院作物科学研究所 | ||
SH064 | RA7 | 中国农业科学院作物科学研究所 | ||
SH065 | RA8 | 中国农业科学院作物科学研究所 | ||
SH066 | RA12 | 中国农业科学院作物科学研究所 | ||
SH067 | RA10 | 中国农业科学院作物科学研究所 | ||
SH068 | RA13 | 中国农业科学院作物科学研究所 | ||
SH069 | S1025 | 黑龙江省农业科学院绥化分院 | ||
SH070 | S007♀ | 黑龙江省农业科学院绥化分院 | ||
SH071 | SY203♂ | 黑龙江省农业科学院绥化分院 | ||
SH072 | SD3♀ | 黑龙江省农业科学院绥化分院 | ||
SH073 | SWM76 | 黑龙江省农业科学院绥化分院 | ||
SH074 | CA193 | 中国农业科学院作物科学研究所 | ||
SH075(Reid) | KWS10/KWS73 | 黑龙江省农业科学院绥化分院 | ||
SH076 | L201 | 国内骨干 | ||
SH077 | SD3♂ | 黑龙江省农业科学院绥化分院 | ||
SH078 | L237 | 国内骨干 | ||
SH079 | SL202 | 黑龙江省农业科学院绥化分院 | ||
SH080 | S65♀ | 黑龙江省农业科学院绥化分院 | ||
SH081 | S65♂ | 黑龙江省农业科学院绥化分院 | ||
SH082 | WC007 | 吉林省中玉农业有限公司 | ||
SH083(Reid) | 郑58 | PA群对照 | ||
SH084(Dom) | 444 | 唐四平头群对照 | ||
SH085(Non-Reid) | 合344 | 兰卡斯特对照 | ||
SH086 | 9176/C8605 | 黑龙江省农业科学院佳木斯分院 | ||
SH087 | SX717 | 黑龙江省农业科学院绥化分院 | ||
SH088 | SX711 | 黑龙江省农业科学院绥化分院 | ||
SH089 | SX721 | 黑龙江省农业科学院绥化分院 | ||
SH090 | SX718/30747 | 黑龙江省农业科学院绥化分院 | ||
SH091 | 4207/C8605/ SX718 | 黑龙江省农业科学院绥化分院 | ||
SH092 | SX619 | 黑龙江省农业科学院绥化分院 | ||
SH093 | 20837 | 黑龙江省农业科学院绥化分院 | ||
SH094 | SD302♂ | 黑龙江省农业科学院绥化分院 | ||
SH095 | SD302♀ | 黑龙江省农业科学院绥化分院 | ||
SH096 | SX706 | 黑龙江省农业科学院绥化分院 | ||
SH097(Dom) | C7-2 | 唐四平头群对照 | ||
SH098 | SD1♀ | 黑龙江省农业科学院绥化分院 | ||
SH099 | 531 | 东北农业大学 | ||
SH100 | WY1M1 | 吉林省中玉农业有限公司 |
1.2 DNA提取
大田中幼苗3叶1心时取嫩叶,采用磁珠法提取100份样品DNA。分别通过琼脂糖电泳和紫外分光光度计检测DNA纯度和浓度等质量状况。选择条带单一、明亮、完整无降解且无拖尾现象的DNA检测琼脂糖电泳,紫外分光光度计检测OD260/OD280值介于1.8~2.0,DNA浓度介于30~120 ng/μL。
1.3 基因型鉴定
采用北京市农林科学院玉米研究中心研制的Maize6H-60K芯片进行基因型分析。Maize6H-60K芯片中有61 214个SNP位点在染色体上的分布,染色体1号到10号上SNP位点数分别是9552、7334、7264、6591、6423、5068、4991、4886、4873、4232个。Affymetrix芯片使用DQC和QC call-rate进行试验质控。DQC是使用已知的、非多态性的marker计算AT通道和GC通道的区分度,DQC是对芯片通道信号值的质控。QC call-rate是选取已测试过的分型效果好的约20K探针评估call rate,QC call-rate是对样品的质控。试验要求QC>90,DQC>0.82。使用Axiom Analysis Suite 3.1.51软件分别依据PolyHighResolution、NoMinorHom、MonoHighResolution、OTVCallRateBelowThreshold、Other分类法对CEL文件进行分型[11]。
1.4 数据处理
利用TASSEL软件完成基因型数据描述性分析和玉米种质资源间遗传距离计算,采用UPGMA方法进行聚类分析,利用MEGA 11软件绘制聚类图,利用R软件计算玉米种质资源间遗传相似度。用GCTA 1.94.0进行主成分分析。
2 结果与分析
2.1 基因型数据的描述性分析
采用Maize6H-60K芯片进行扫描,获得100个样本的61 224个SNP位点的原始数据,将这些数据利用Axiom Analysis进行基因型分析,结果(表2)显示,61 224个SNP位点有6种类型,其中分型较好的有PHR类型,为39 381个,MHR为178个,NMH为127个,这些数据可进行各种遗传应用分析;其中分型不好的有OTV类型,为11 216个,CRBT类型1373个,其他类型8939个,这些SNP位点未能获得准确可区分的基因型数据。
表2 61 224个SNP位点的6种类型
Table 2
转换型Conversion type | 数量Number | 百分比Percentage (%) |
---|---|---|
PHR | 39 381 | 64.333 |
OTV | 11 216 | 18.323 |
其他Other | 8939 | 14.603 |
CRBT | 1373 | 2.243 |
MHR | 178 | 0.291 |
NMH | 127 | 0.207 |
PHR类型的39 381个SNP在100份玉米资源的基因型最小等位变异频率变化范围为0~0.5,平均值为0.3146;缺失数据比例变化范围0~0.1200,平均值为0.01095;杂合度比例变化范围为0~ 0.70213,平均值为0.03846。
2.2 供试材料聚类分析
聚类分析结果(图1)表明,以6份测验种为参考,100份玉米资源被划分为3大类,其中类群Ⅰ为Reid群,占77%,包含KWS10/KWS73/欧早、SX618/欧早、S1029、SI15412、S1035等77份玉米资源;类群Ⅱ为Non-Reid群,占19%,包含344/9229、98-4/郑58、6171/SX711等19份玉米资源;类群Ⅲ为Dom群,共4份资源,占4%,包含W08、L237、444和C7-2。
图1
图1
基于遗传距离的100份玉米自交系聚类分析图
Fig.1
Cluster analysis diagram of the 100 maize inbred lines based on genetic distance
2.3 供试材料遗传相似度分析
图2
2.4 供试材料主成分分析
根据遗传相似矩阵进行主成分分析(图3),PC1和PC2的贡献率分别为30.43%和20.51%,累计贡献率50.94%。根据100份自交系的集中程度,可清晰地鉴定出3大类群,这与群体遗传距离聚类分析结果基本一致。
图3
图3
根据遗传相似矩阵对100份玉米自交系进行主成分分析
Fig.3
Principal component analysis of the 100 maize inbred lines based on genetic similarity matrix
3 讨论
研究者[16]常根据材料的表型、系谱来源和分子标记等方面对玉米种质资源进行优势类群划分,目前利用分子标记技术划分玉米杂种优势类群更方便可靠[17]。在分子标记技术中,SNP标记技术由于具有高通量检测和高密度分布等优势,常被应用于划分玉米杂种优势类群[18⇓⇓⇓⇓⇓-24]。目前在我国北方的生产应用中,育种家主要利用Reid、Non-Reid(以Lancaster为主)和Dom群(塘四平头群和旅大红骨类群)这三大优势类群的材料,本研究利用Maize6H-60K SNP芯片对100份玉米自交系进行基因型分析,计算出100份玉米自交系之间的遗传距离,根据材料间遗传距离进行聚类分析,将这100份种质资源划分为3个类群(Reid、Non-Reid和Dom)。其中划分为Reid类群的自交系最多,占77%,而Dom类群只有4%,这是因为近些年黑龙江地区的玉米育种目标主要为高产、抗病、抗倒伏、脱水快和适于机械化收获。Reid类群的特点多为植株叶片较上冲、果穗大、籽粒马齿型、抗病、抗倒伏和配合力高,因此该类群材料在育种工作中使用较多,血缘偏向于Reid类群的种质资源占比较大。Dom类群为国内玉米种质资源,包括塘四平头群和旅大红骨类群,在黑龙江地区的玉米育种中的使用相对较少,但是该类群有很大的利用潜力,比如塘四平头群在耐密方面有优势,旅大红骨群的果穗大、配合力强。
遗传相似度分析可以更进一步研究种质资源间的遗传关系。相同杂种优势类群的玉米自交系的遗传相似度可达到60%以上,遗传相似度为50%左右或以下代表血缘关系不明显[25]。当种质间遗传相似度小于0.9时,这2个种质材料可当不同的自交系应用于育种工作。当遗传相似度大于0.9或接近1.0,代表2份材料遗传背景相似度过高,育种上应尽量避免使用[13]。本研究遗传相似度结果表明,100份玉米资源遗传相似度分布范围为0.3380~0.9928。有4对玉米材料之间的遗传相似度超过了0.9,其中SH071和SH100之间的遗传相似度大于0.99。说明SH71和SH100可能为同一来源。通过遗传相似分析,对这100份玉米自交系的血缘关系的了解更加清晰,为以后的育种工作避免了一些无效组合。为了得到更加直观的遗传相似度分析结果,对这100份自交系的遗传相似度进行主成分分析,根据100份自交系的集中程度,可清晰地鉴定出3个优势类群,这与群体遗传距离的聚类分析结果基本一致。通过聚类分析、遗传相似度分析和主成分分析,对这100份种质资源遗传关系的认知变得更加清晰,为利用和改良这些自交系提供参考。
4 结论
通过SNP芯片标记进行100份玉米种质资源遗传多样性分析,根据遗传距离进行聚类分析,划分为3个类群(Reid、Non-Reid和Dom)。利用R软件计算100份玉米资源间遗传相似度,分布范围为0.3380~0.9928,发现有4对玉米材料之间的遗传相似度超过了0.9,这些材料在育种工作中使用时需要注意。对玉米材料间遗传相似度进行主成分分析,与遗传距离的聚类分析结果基本一致。
参考文献
SSR标记在玉米杂种优势类群划分中的应用
Genomic insights into historical improvement of heterotic groups during modern hybrid maize breeding
De novo genome assembly and analyses of 12 founder inbred lines provide insights into maize heterosis
The wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat
Development of a maize 55K SNP array with improved genome coverage for molecular breeding
Potential to improve oilseed rape and canola breeding in the genomics era.
Crop breeding chips and genotyping platforms: progress, challenges, and perspectives
/
〈 |
|
〉 |
