作物杂志, 2021, 37(4): 51-58 doi: 10.16035/j.issn.1001-7283.2021.04.008

遗传育种·种质资源·生物技术

120份大豆种质资源遗传多样性和亲缘关系分析

李琼,1, 常世豪1, 武婷婷2, 耿臻,1, 杨青春1, 舒文涛1, 李金花1, 张东辉1, 张保亮1

1周口市农业科学院,466001,河南周口

2中国农业科学院作物科学研究所,100081,北京

Analysis of Genetic Diversity and Genetic Relationship for 120 Soybean Germplasms

Li Qiong,1, Chang Shihao1, Wu Tingting2, Geng Zhen,1, Yang Qingchun1, Shu Wentao1, Li Jinhua1, Zhang Donghui1, Zhang Baoliang1

1Zhoukou Academy of Agricultural Sciences, Zhoukou 466001, Henan, China

2Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China

通讯作者: 耿臻,研究方向为大豆遗传育种与栽培,E-mail:gengzhen0616@163.com

收稿日期: 2021-03-14   修回日期: 2021-05-24   网络出版日期: 2021-07-05

基金资助: 河南省甘薯杂粮体系大豆岗位专家(S2020-14-G01)

Received: 2021-03-14   Revised: 2021-05-24   Online: 2021-07-05

作者简介 About authors

李琼,研究方向为大豆遗传育种与栽培,E-mail:15290067998@163.com

摘要

采用覆盖20条染色体的68对多态性引物在120个大豆品种(系)[包括67个中国大豆品种(系)和53个引自美国的品种(系)]进行遗传多样性及亲缘关系分析。结果表明,120份大豆材料具有丰富的遗传多样性(Na=1.9852,Ne=1.4343,H=0.2776,I=0.4361);中国大豆比美国引进大豆遗传多样性水平高;7个群体遗传多样性由高到低为,黄淮海地区组>引种资源组>热带亚热带地区组>长江流域地区组>北方春大豆组>西南山区组>鲜食大豆组;7个群体间总遗传多样度(Ht)为0.2807,群体内遗传多样度(Hs)为0.2361,群体间遗传分化系数(Gst)为0.1589,基因流(Nm)为2.6468,群体间存在中低度遗传分化,遗传变异主要存在于群体内部,群体间Nm较丰富;以群体和单个品种(系)为单位进行UPGMA聚类的结果基本一致,部分材料相互交错。大豆种质遗传多样性与地理来源具有一定相关性的同时,不同地区间存在丰富的基因交流;黄淮海地区、西南山区、热带亚热带地区和长江流域地区的大豆材料遗传距离较近;引进大豆、北方春大豆和鲜食大豆与其他群体遗传距离较远,可作为拓宽中国栽培大豆遗传背景的物质遗传基础。

关键词: 大豆; 种质资源; 遗传多样性; SSR

Abstract

Using 68 SSR markers covering 20 chromosomes, 120 soybean varieties (lines) including 67 Chinese soybean varieties (lines) and 53 introduced soybean varieties (lines) from America were used for genetic diversity and genetic relationship analysis. The results showed that the average values of the number of alleles (Na), effective alleles (Ne), Nei's gene diversity index (H), and Shannon's information index (I) between the varieties (lines) were 1.9852, 1.4343, 0.2776, and 0.4361, respectively, indicating the genetic diversity of 120 soybean germplasms was high. Genetic diversity in Chinese soybeans was higher than introduced soybeans from America. The genetic diversity index order of the seven populations was Huang-Huai-Hai area group > introduction resource group > tropical and subtropical area group > Yangtze River area group > northern spring soybean group > southwest mountainous group > fresh soybean group. The total genetic diversity (Ht) among all seven populations was 0.2807 while the genetic diversity (Hs) within populations was 0.2361, the coefficient of genetic differentiation (Gst) among populations was 0.1589 and the gene flow (Nm) was 2.6468. Low-medium genetic differentiation was identified among populations while genetic variation mainly existed within each group and gene flow in seven populations was abundant. The results of UPGMA clustering analyses using different populations and individual varieties (lines) as the unit were basically the same and some materials were interlaced in different groups. While soybean germplasms had a certain correlation with geographic origin, there were abundant gene exchanges among different regions. The genetic distances were close between varieties (lines) from the Huang-Huai-Hai area, southwest mountainous areas, tropical and subtropical area, and the Yangtze River area. However, the genetic distance between the varieties (lines) from America, northern spring varieties (lines), and fresh soybean varieties (lines) were distant. In consequence, these materials could provide favorable germplasms for broadening the genetic background.

Keywords: Soybean; Germplasm resources; Genetic diversity; SSR

PDF (1130KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李琼, 常世豪, 武婷婷, 耿臻, 杨青春, 舒文涛, 李金花, 张东辉, 张保亮. 120份大豆种质资源遗传多样性和亲缘关系分析. 作物杂志, 2021, 37(4): 51-58 doi:10.16035/j.issn.1001-7283.2021.04.008

Li Qiong, Chang Shihao, Wu Tingting, Geng Zhen, Yang Qingchun, Shu Wentao, Li Jinhua, Zhang Donghui, Zhang Baoliang. Analysis of Genetic Diversity and Genetic Relationship for 120 Soybean Germplasms. Crops, 2021, 37(4): 51-58 doi:10.16035/j.issn.1001-7283.2021.04.008

大豆起源于中国,已有5000年栽培历史。大约于秦代自我国华北引入朝鲜,又传入日本,18世纪后逐渐传播到美国和欧洲的一些国家,成为世界上主要的油料和饲料作物之一[1,2,3]。大豆种质资源是新品种选育的物质基础,挖掘丰富的遗传变异资源,探索国内外大豆种质资源多样性,对拓宽大豆遗传背景和创制大豆优良新种质具有重要意义。大豆育成品种在我国大豆生产中发挥着重要作用,是重要的大豆种质资源。栽培大豆是自然选择与人工选择共同作用的结果[4]。人工选育栽培品种过程中对优良品种的推广和应用,造成地方品种和农家种的逐步流失。对关键性亲本种质的重复利用导致栽培大豆遗传基础越来越窄,造成大豆突破性品种难以育成。

基于SSR分子标记技术对各类大豆种质资源进行发掘、研究,采用野生大豆和国外引进品种(系)等亲缘关系较远的大豆材料来拓宽国内大豆品种(系)的遗传基础,打破我国大豆资源丰富但遗传基础狭隘的局面[5,6]。钟文娟等[7]利用53对SSR标记对190份我国四川大豆和引进大豆资源遗传多样性和群体结构进行分析,证明我国四川大豆可以利用国外资源与直立型野生大豆资源来丰富本地的基因遗传多样性。Wang等[8,9]利用SSR标记分别对中国大豆资源七大生态区域类型的遗传多样性和群体结构进行分析,得出中国大豆种质资源遗传多样性在某种程度上与地理距离相关。Iquira等[10]利用39对SSR引物对100份加拿大本地品种和200份外来种质进行遗传多样性分析,得知外来种质比本地材料具有更高的遗传多样性。本研究利用68对SSR引物对120份大豆种质材料(中国和美国大豆种质)遗传多样性和群体亲缘关系进行探索,为大豆新品种选育提供相应的理论参考。

1 材料与方法

1.1 试验材料

供试材料共120份,其中67份选自2014-2015年国家大豆区域6个参试试验组,53份育成品种(系)来自美国,均由中国农业科学院作物科学研究所提供。67份国内参试材料按照参试区域分为6个群体,群体1(POP-1)为北方春大豆组,群体2(POP-2)为黄淮海地区组,群体3(POP-3)为长江流域地区组,群体4(POP-4)为鲜食大豆组,群体5(POP-5)为西南山区组,群体6(POP-6)为热带亚热带地区组。将53份引进材料设为群体7(POP-7)。详情见表1

表1   供试大豆种质详情

Table 1  Details of the tested soybean germplasms

群体Population品种(系)Cultivar (line)来源Source
北方大豆组合交2001-336-7、绥交08-5262、牡602中国黑龙江
Northern spring soybean group (POP-1)公交04-489-4中国吉林
齐黄35中国山东
黄淮海地区组K117-3、中品12585、中作J10153、中黄13中国北京
Huang-huai-hai area group (POP-2)晋科1号、汾豆90中国山西
宁豆5号中国宁夏
沧豆11、石豆6号、邯豆5号中国河北
郑豆0689、秋乐1205、濮豆1802、周豆22、许豆0406、驻豆03-16、安豆5156中国河南
蒙0804、皖宿1128、益科豆112-1中国安徽
潍豆80031、鲁97013-1、菏01-7、济087129、齐黄34、圣豆LU05、南圣43中国山东
徐9418-2中国江苏
长江流域地区组苏HT012、泗豆209、苏夏13-13、南农56-10、通豆7号中国江苏
Yangtze River basin area group (POP-3)兴豆2号、中豆8号、油07-73、天隆一号中国湖北
秋乐1302中国河南
蒙11-1中国安徽
鲜食大豆组Fresh soybean group (POP-4)杭豆3-2、浙9815、浙鲜豆5号、浙09006、杭鲜豆2号中国浙江
K丰77-2、开科源特早、辽04M05-4中国辽宁
西南山区组安07109中国贵州
Southwest mountain group (POP-5)滇豆7号中国云南
黔豆2012-1中国贵州
贡2026R中国四川
粤春2011-2中国广东
热带亚热带地区组贡豆723-2中国四川
Tropical and subtropical area group (POP-6)华春6号、粤春2010-1、华春2号、粤夏2012-1中国广东
桂早1号、桂春豆103、桂147中国广西
泉豆7号、泉豆5号中国福建
引进材料Imported materials (POP-7)AD01~AD53美国

新窗口打开| 下载CSV


1.2 试验方法

1.2.1 大豆DNA的提取 每份材料选取2粒饱满的种子于育苗盘中种植,放置在培养箱内(25℃,光照8h/d)培育。采集鲜嫩叶片,按照CTAB法[11]提取DNA。采用紫外分光光度计检测DNA浓度,并稀释至50ng/μL,放置于-80℃冰箱备用。

1.2.2 引物 参考http://soybase.org/公布的引物标记,在标记中选取均匀分布于20条染色体上的68对SSR引物。SSR标记引物由上海生物工程有限公司合成。

1.2.3 PCR扩增及电泳 采用10μL PCR反应体系:2×TranTaq-PCRSuper-Mix 5μL,上、下游引物各0.5μL,DNA模板1μL,ddH2O 3μL。PCR反应程序:94℃预变性5min;94℃变性40s,55℃退火50s,72℃延伸50s,循环34次;72℃延伸10min,4℃保存。PCR产物用6%变性聚丙烯酰胺凝胶电泳检测,电压200V,时长2.5h,银染法染色,即时拍照。

1.3 数据分析

对扩增条带进行人工读带并记录带型,同一位置有条带赋值“1”,无条带则赋值“0”。利用POPGEN 32软件对所选群体的等位基因数(allele number,Na)、有效等位基因数(effective number of allele,Ne)、Shannon’s信息指数[12](Shannon’s information index,I)、遗传相似度(genetic identity,GI)和遗传距离(genetic distance,GD)、多态条带百分比(percentage of polymorphic bands,PPB)、Nei’s基因多样性指数(H)、总遗传多样度(Ht)、群体内遗传多样度(Hs)、群体间遗传分化系数(Gst)和基因流(Nm)进行分析。计算多态性信息含量(polymorphism information content,PIC[13]。利用NTSYS-PC软件计算遗传相似系数(genetic similarity,Gs),按照非加权配对法UPGMA和SHAN程序进行聚类分析。

2 结果与分析

2.1 SSR引物标记多态性分析

利用均匀分布于20条染色体的68对多态性SSR引物(表2)对120份供试材料进行扩增,共得到203条多态性条带,扩增片段长度在117~340bp之间,平均每个引物扩增条带3.0个,多态性条带为200个。其中,satt343和satt373多态性位点最多。多态性信息含量(PIC)变幅为0.290~0.991,平均0.803。朱振东等[14]提出PIC大于0.5时,SSR引物为高度信息引物。由表2可知,有64对(占比94.12%)引物位点为高度多态性位点。

表2   供试材料SSR引物多态性信息

Table 2  Polymorphism information of SSR primers in the tested materials

编号
No.
引物
Primer
连锁群
Linkage group
多态位点数
Polymorphic locus
PIC
PIC value
编号
No.
引物
Primer
连锁群
Linkage group
多态位点数
Polymorphic locus
PIC
PIC value
1Sat_143G40.91435satt272B220.470
2Sat_292D240.87136satt279H30.865
3Sat_219I40.91137satt286C230.865
4Sat_304N40.89938satt300A130.822
5Sat_282O20.81039satt302H20.760
6Sat_271A120.82640satt307C230.869
7Satt130G20.99141satt309G30.874
8Satt132J20.40042satt334F20.823
9Satt142H20.29043satt335F20.786
10Satt146F30.88344satt338C140.904
11Sat_020K30.90145satt339N30.891
12Sat_076C240.77446satt343F50.862
13sat_207C130.51847satt346M20.725
14Sat_342B240.91548satt373L50.931
15Sat_292D230.87149satt386D220.560
16Satt001J40.84750satt423F20.805
17Satt002D230.82251satt429A220.687
18Satt005D1b30.82252satt431J30.776
19Satt022N20.80053satt429A230.875
20Satt168B230.57754satt434H40.879
21Satt182L20.88055satt453B130.893
22Satt184D1a30.77556satt463M40.899
23Satt187A220.73057satt530N30.873
24Satt191G30.74858satt553D230.869
25Satt194C120.46759satt556B230.642
26Satt216D1b40.90060satt565C140.893
27Satt199G30.85561satt571I30.867
28Satt197B130.84062satt573D220.768
29Satt228A230.74463satt577B230.848
30Satt239I20.75264satt588K30.821
31Satt243O30.82365satt590M40.843
32satt254D1a30.76466satt599A130.895
33satt257N20.78967satt631N40.919
34satt268D240.92068satt708C230.842

新窗口打开| 下载CSV


2.2 不同大豆种质群体间的遗传变异分析

表3可知,供试的120份大豆种质等位基因数(Na)为1.9852,有效等位基因数(Ne)为1.4343,Nei’s基因多样性指数(H)为0.2776,Shannon’s信息指数(I)为0.4361,多态条带百分比(PPB)为98.52%。遗传多样性分析结果表明,120份大豆材料总体遗传多样性水平较高。其中,国内材料与美国引进大豆材料遗传多样性水平差异较为明显,国内大豆种质资源遗传多样性水平更高(Ne=1.4196,H=0.2710,I=0.4282),美国引进资源的遗传多样性略低(Ne=1.3985,H=0.2501,I=0.3916)。

表3   7个群体间大豆种质材料的遗传多样性指标

Table 3  Genetic diversity indexes of soybean germplasm materials between seven populations

群体
Population
NaNeHI多态位点数
Number of polymorphic sites
PPB
(%)
HtHsGst基因流Nm
POP-1~POP-61.98031.41960.27100.428219998.03
POP-71.92611.39850.25010.391618892.61
总体Total1.98521.43430.27760.436120098.520.28070.23610.15892.6468

新窗口打开| 下载CSV


通过Nei’s分析法将总的遗传多样性分为群体内遗传多样性和群体间遗传多样性,遗传分化系数是指群体间变异占总变异量的比重[15]。7个群体间总遗传多样度(Ht)为0.2807,群体内遗传多样度(Hs)为0.2361。群体间遗传分化系数(Gst)为0.1589,基因流(Nm)为2.6468。由此可见,供试大豆资源具有丰富的多样性。大豆材料7个群体遗传变异中,84.11%的遗传变异存在于7个群体内部个体间,15.89%的遗传变异存在7个群体之间。7个群体间存在中低度遗传分化,遗传变异主要存在于群体内部,群体间基因流较为丰富。

对7个群体遗传多样性进一步分析(表4),不同群体各遗传多样指标的排列顺序为POP-2(Ne=1.4208,H=0.2678,I=0.4209)>POP-7(Ne=1.3985,H=0.2501,I=0.3916)>POP-6(Ne=1.3949,H=0.2424,I=0.3753)>POP-3(Ne=1.3850,H=0.2403,I=0.3741)>POP-1(Ne=1.3734,H=0.2254,I=0.3430)>POP-5(Ne=1.3611,H=0.2193,I=0.3329)>POP-4(Ne=1.3429,H=0.2075,I=0.3194)。

表4   7个群体大豆种质材料的遗传多样性指标

Table 4  Genetic diversity indexes of seven populations of soybean germplasm materials

群体PopulationNaNeHI
POP-11.67491.37340.22540.3430
POP-21.96061.42080.26780.4209
POP-31.82271.38500.24030.3741
POP-41.67981.34290.20750.3194
POP-51.64531.36110.21930.3329
POP-61.81771.39490.24240.3753
POP-71.92611.39850.25010.3916

新窗口打开| 下载CSV


2.3 不同大豆种质群体间的遗传相似度和UPGMA聚类分析

表5可知,整体来看,群体间遗传相似度(GI)变幅为0.8950~0.9749,群体间遗传距离(GD)变幅为0.0254~0.1109,7个群体间相似度高,遗传距离小。从群体之间遗传关系来看,POP-1与POP-5 GD最大,POP-2与POP-3最小。

表5   不同大豆群体间的遗传相似度(GI)和遗传距离(GD

Table 5  Genetic similarity (GI) and genetic distance (GD) between different soybean populations

群体
Population
POP-1POP-2POP-3POP-4POP-5POP-6POP-7
POP-1*0.93990.92020.89920.89500.92360.9485
POP-20.0620*0.97490.92570.95330.96220.9633
POP-30.08320.0254*0.90680.95720.95760.9345
POP-40.10620.07720.0978*0.89720.89880.9205
POP-50.11090.04780.04380.1084*0.93800.9160
POP-60.07950.03850.04340.10670.0640*0.9475
POP-70.05290.03740.06780.08290.08780.0539*

右上为遗传相似度(GI),左下为遗传距离(GD

Nei’s genetic identity (above diagonal) and genetic distance (below diagonal)

新窗口打开| 下载CSV


根据7个群体间的Nei’s遗传相似度进行UPGMA聚类分析结果(图1)可知:7个群体分为两大类,不同群体按照1、2、3、4、5和6的步骤演示出群体之间的亲缘关系。其中,类群I包括POP-1、POP-2、POP-3、POP-5、POP-6和POP-7;类群I可分为2个亚类,亚类I-i包括POP-1和POP-7,亚类I-ii包括POP-2、POP-3、POP-6和POP-5。类群II仅包括POP-4 1个群体。总体而言,黄淮海地区组(POP-2)、长江流域地区组(POP-3)、西南山区组(POP-5)和热带亚热带地区组(POP-6)亲缘关系较近。北方春大豆(POP-1)、引进大豆(POP-7)和鲜食大豆(POP-4)与这4个群体的遗传距离较远。

图1

图1   7个群体的大豆种质材料的聚类图

Fig.1   UPGMA tree plot of soybean germplasm materials of seven populations


2.4 120份大豆材料聚类分析

利用NTSYS-PC软件分析68对SSR引物检测出的200个多态性位点数据,计算并建立遗传相似度矩阵,进行UPGMA聚类。由图2可知,120份黄淮海大豆种质的遗传相似系数(Gs)变化范围为0.62~0.98。遗传相似系数(Gs)在0.62处,120份材料被分为2个类群。类群I包括71份材料,分为2个亚类。其中,亚类I-i包括美国引进材料51份、北方春大豆全部材料5份、黄淮海地区材料5份(鲁97013-1、沧豆11、淮豆80031、中品12585和中作J10153)、西南山区组材料1份(贡2026R);亚类I-ii包括鲜食大豆组全部材料8份和美国引进材料1份(AD31)。类群II包括49份材料,分为2个亚类。其中,亚类II-i包括黄淮海地区组材料23份、长江流域地区组材料10份、西南山区材料4份和热带亚热带地区组材料4份;亚类II-ii包括热带亚热带地区组材料6份、长江流域地区材料1份(南农56-10)和引进材料1份(AD48)。以大豆品种(系)为单位的聚类分析进一步印证了2.3部分的结论。

图2

图2   120份大豆种质资源SSR标记UPGMA聚类分析

Fig.2   UPGMA tree of SSR marker cluster analysis of 120 soybean germplasm resources


3 讨论

3.1 供试大豆种质资源遗传多样性和亲缘关系

本试验选取覆盖大豆全基因组的68对SSR标记,对120份大豆种质材料进行遗传多样性分析,以探索中国和美国大豆资源中各群体的遗传背景和群体间的亲缘关系。本研究所选取的68对引物具有较高的多态性。7个群体大豆资源(Na=1.9852,Ne=1.4343,H=0.2776,I=0.4361)总体具有丰富的多样性,多态条带百分比为98.52%。国内大豆种质资源(Ne=1.4196,H=0.2710,I=0.4282)比引进资源(Ne=1.3985,H=0.2501,I=0.3916)遗传多样性水平高。我国是大豆的起源地[16],大豆种质资源较美国更丰富。

3.2 供试大豆种质资源不同群体遗传多样性比较

文自翔[17]研究表明群体遗传多样性是长期形成的产物,遗传多样性越高,遗传变异越丰富,对环境变化的适应能力就越强。国内的6个群体中黄淮海地区组大豆材料的遗传多样性最高,热带亚热带地区组次之,长江流域地区组、北方春大豆组、西南山区组和鲜食大豆组遗传多样性较低。这与以往研究结果一致。朱申龙等[18]认为我国黄淮海地区的大豆较其他地区具有较高水平的遗传多样性。周新安等[19]对国家种质库保存的9个形态性状的遗传多样性分析,判定黄淮夏大豆遗传多样性较为丰富。王彩洁等[20]利用与大豆性状相关的125对SSR标记对东北和黄淮海地区大面积推广的89个大豆品种进行遗传多样性分析,得出黄淮海地区品种的遗传多样性最高,黄淮海流域位于我国南、北大豆产区的交汇处,无霜期较长,大豆品种多样,是不同来源大豆品种相互交融的理想地带。由7个群体遗传多样性比较可知,美国引进大豆资源遗传多样性指数紧随我国黄淮海地区之后,位居第二。美国引进大豆种植区域主要分布在我国西部和北部,育种所用亲本主要与我国少数几个原始亲本有关,比如Mandarin(黑龙江绥化四粒黄)、Manchu(黑龙江宁安黄豆)和Mukedn(沈阳小金黄)等,推测美国大豆资源与中国大豆东北材料具有密切的亲缘关系。栽培大豆起源于中国,且大豆黄河中下游起源学说得到广泛的支持。因此,美国引进资源的遗传多样性相对于黄淮海地区组较低也就有迹可循。美国引进资源的国外的祖先亲本量占祖先总亲本量的50%,而我国各地区大豆育成品种的亲本大部分来源于本地区[21]。推测这就是美国亲本杂交繁育出品种的遗传多样性比除黄淮海地区以外中国其他各地区组均丰富的原因。但各地区参试品种(系)数量存在差异,造成7个群体样本量有所不同,对各群体遗传多样性分析结果存在一定程度上的影响,因此仍需搜集更多的种质资源进一步分析验证。

群体间总遗传多样度(Ht=0.2807)大于群体内遗传多样度(Hs=0.2361)。群体间存在中低度遗传分化(Gst=0.1589),遗传变异主要存在于群体内部,各群体间有频繁基因交流(Nm=2.6468)。我国育成的新品种群体遗传基础越来越窄,导致育成品种群体的遗传脆弱性和单一性[22]。Burnham等[12]和Ude等[23]研究也发现中国和美国大豆亲缘关系较近。Nichols等[24]对中国及北美大豆基因型的变异性进行研究,认为中、美栽培大豆品种的遗传基础比较狭窄。

3.3 供试大豆种质资源两种聚类结果比较

本研究分别以群体和单个品种(系)为单位进行UPGMA聚类比较。结果相似点有:整体来看,中国北方春大豆全部材料(100%)和美国引进资源大部分材料(96.23%)均聚在一起;中国黄淮海地区大部分材料(82.14%)、西南山区大部分材料(80%)、长江流域地区和热带亚热带地区全部材料(100%)聚在一起。2种聚类结果的不同点有:(1)鲜食大豆组与其他6个群体亲缘关系均远,单独分为一类。按品种(系)聚类可知,中国鲜食大豆组材料与美国引进材料、中国北方春大豆材料和中国部分黄淮海地区、中国西南山区材料聚为一类。2次聚类结果不一致,可能是由于鲜食大豆群体取材数较少,品种单一,遗传变异小,多态性低。(2)品种(系)个体聚类可知,2份美国引进材料(AD31和AD48)分布在I-ii亚类(以鲜食大豆为主)中和II-ii亚类(以热带亚热带地区大豆材料为主)中;中国5份黄淮海地区材料(鲁97013-1、沧豆11、淮豆80031、中品12585和中作J10153)和1份西南山区大豆材料(贡2026R)分布在亚类I-i(以美国进口大豆材料为主)中;1份长江流域地区大豆材料(南农56-10)分布在亚类II-ii(以热带亚热带地区大豆材料为主)中。(3)品种(系)个体聚类,黄淮海地区、西南山区、热带亚热带地区和长江流域地区4个群体中大豆材料,只有黄淮海地区部分材料和热带亚热带地区部分材料分别聚到一起,其余材料相互交错,这4个群体大豆品种(系)的亲缘关系与地理来源分布较为分散,表明这4个群体间大豆育成品种种质的基因交流较为频繁。

2次聚类分析结果基本一致,体现出大豆种质资源与地理分布在一定程度上具有显著的相关性。但是部分材料也存在相互融合,从而展现出不同大豆群体之间广泛的基因交流。

3.4 供试大豆种质资源不同群体遗传相似度和遗传距离比较

对120份大豆材料遗传相似度和遗传距离进行研究,得出7个大豆群体的平均遗传距离仅为0.0700。也许与所选引物个数少和品种(系)材料多样性有限有关。材料组间的遗传距离值差异虽小,但仍可明显地发现各群体之间具有一定的差异。大豆种质不同群体遗传相似度、遗传距离的分析结果表明,黄淮海地区材料和长江流域地区材料之间遗传相似度最高(0.9749),遗传距离最小(0.0254)。北方春大豆材料和西南山区材料遗传相似度最低(0.8950),遗传距离最大(0.1109)。推测黄淮海和长江流域地区因地理位置相邻,大豆栽培品种(系)频繁交流造成该片区遗传相似度高,遗传距离缩小的局面。北方地区与西南山区地理位置相距较远,北方地区地域辽阔且西南山区地形复杂均能构成天然屏障,减少了遗传背景相似的可能。除北方春大豆材料以外,其他6个群体均与黄淮海地区材料间遗传距离较近。Wang等[8]鉴定出黄河流域下游地区的种质等位基因丰富度最高,成对遗传多样性估计值最低,并且散布在其他5个群体中,这表明黄河流域可能是中国栽培大豆的多样性中心。然而,盖钧镒等[25]则认为从遗传多样性反推大豆起源中心是有一定风险的。我国北方春大豆材料与美国材料遗传距离相距最近。田清震等[26]曾利用AFLP数据将栽培大豆在地理分布上明显区分南北类型,北方地区和黄淮海区大豆可明确区分,表明品种间遗传关系与地理分布密切相关。盖钧镒[21]曾表示北方地区特别是东北地区大豆种质的南移对拓宽国内大豆遗传背景具有重要意义。盖钧镒等[25,27]研究表明美国及南美洲的大豆在遗传和亲缘关系上与中国北方春大豆最接近。鲜食大豆材料与其他群体材料平均遗传距离较远,但与我国黄淮海地区和美国材料遗传距离较近。我国南方(长江中下游、台湾和东南沿海地区)和日本是世界菜用大豆的主要生产和消费地区。中国传统菜用大豆生产主要采用地方品种(系)。目前,我国东南沿海地区主要种植中国台湾和日本品种(系)。20世纪90年代后期,美国逐渐引进菜用大豆,大部分来自中国,少部分来自泰国和印度尼西亚[28]。推测鲜食大豆材料与中国黄淮海地区和美国个别材料遗传背景相似。

4 结论

用SSR技术对来自中国和美国的120份大豆品种(系)及由其组成的7个群体的遗传背景进行分析,得出中国大豆总体遗传多样性较美国大豆高,美国引进大豆资源遗传多样性稍低于中国黄淮海地区。我国黄淮海地区、长江流域地区、西南山区和热带亚热带地区组这4个群体的大豆种质资源交流频繁,遗传背景极为相似。我国北方春大豆、美国引进大豆和我国鲜食大豆与这4个群体遗传距离最远,可作为拓宽我国栽培大豆遗传背景的遗传基础。

参考文献

赵团结, 盖钧镒.

栽培大豆起源与演化研究进展

中国农业科学, 2004,37(7):954-962.

[本文引用: 1]

Hymowitz T, Newell C A.

Taxonomy of the genus Glycine,domestication and uses of soybeans

Economic Botany, 1981,35(3):272-288.

[本文引用: 1]

夏正俊.

大豆基因组解析与重要农艺性状基因克隆研究进展

植物学报, 2017,52(2):148-158.

[本文引用: 1]

郭文韬.

略论中国栽培大豆的起源

南京农业大学学报(社会科学版), 2004,4(1):60-69.

[本文引用: 1]

傅旭军.

大豆种质资源遗传多样性分析及其利用研究

杭州:浙江大学, 2009.

[本文引用: 1]

蒲艳艳, 宫永超, 李娜娜.

中国大豆种质资源遗传多样性研究进展

大豆科学, 2018,37(2):315-321.

[本文引用: 1]

钟文娟, 袁灿, 周永航, .

基于SSR标记的四川大豆与引进大豆资源遗传多样性和群体结构分析

大豆科学, 2017,36(5):657-668.

[本文引用: 1]

Wang L X, Guan R X, Liu Z X, et al.

Genetic diversity of Chinese cultivated soybean revealed by SSR markers

Crop Science, 2006,46(3):1032-1038.

[本文引用: 2]

Wang L X, Guan R X, Liu Z X, et al.

Genetic diversity of Chinese spring soybean germplasms revealed by SSR markers

Plant Breeding, 2008,127(1):56-61.

[本文引用: 1]

Iquira E, Gagnon E, Belzile F.

Comparison of genetic diversity between Canadian adapted genotypes and exotic germplasm of soybean

Genome, 2010,53(5):337-345.

[本文引用: 1]

Doyle J J.

Isolation of plant DNA from fresh tissue

Focus, 1990,12(1):13-15.

[本文引用: 1]

Burnham K D, Francis D M, Dorrance A E, et al.

Genetic diversity patterns among phytophthora resistant soybean plant introductions based on SSR markers

Crop Science, 2002,42(2):338-343.

[本文引用: 2]

Zhang J, Hu K, Li K J, et al.

Simulating the effects of long-term discontinuous and continuous fertilization with straw return on crop yields and soil organic carbon dynamics using the DNDC model

Soil and Tillage Research, 2017,165:302-314.

[本文引用: 1]

朱振东, 王化波, 王晓鸣, .

黑龙江主要栽培大豆品种(系)对大豆疫霉根腐病的多抗性评价

植物遗传资源学报, 2004,5(1):22-25.

[本文引用: 1]

王衍莉, 杨义明, 范书田, .

基于SSR分子标记的73份山葡萄及杂交后代的遗传多样性分析

生物技术通报, 2021,37(1):189-197.

[本文引用: 1]

石慧, 王思明.

大豆在中国的历史变迁及其动因探究

农业考古, 2019(3):32-39.

[本文引用: 1]

文自翔.

中国栽培和野生大豆的遗传多样性、群体分化和演化及其育种性状QTL的关联分析

南京:南京农业大学, 2008.

[本文引用: 1]

朱申龙, Morti M L, Rao R.

应用AFLP方法研究中国大豆的遗传多样性

浙江农业学报, 1998,10(6):302-309.

[本文引用: 1]

周新安, 彭玉华, 王国勋, .

中国栽培大豆遗传多样性和起源中心初探

中国农业科学, 1998,31(3):37-43.

[本文引用: 1]

王彩洁, 孙石, 金素娟, .

中国大豆主产区不同年代大面积种植品种的遗传多样性分析

作物学报, 2013,39(11):1917-1926.

[本文引用: 1]

盖钧镒. 作物遗传育种学各论:第2版. 北京: 中国农业出版社, 1997.

[本文引用: 2]

宋喜娥, 李英慧, 常汝镇, .

中国栽培大豆[Glycine max (L.) Merr.]微核心种质的群体结构与遗传多样性

中国农业科学, 2010,43(11):2209-2219.

[本文引用: 1]

Ude G N, William J, Kenworthy , et al.

Genetic diversity of soybean cultivars from China,Japan,North America,and North American ancestral lines determined by amplified fragment length polymorphism

Crop Science, 2003,43(5):1858-1867.

[本文引用: 1]

Nichols D M, Wang L Z, Yan L P, et al.

Variability among Chinese Glycine soja and Chinese and North American soybean genotypes

Crop Science, 2007,47(3):1289-1298.

[本文引用: 1]

盖钧镒, 赵团结.

中国大豆育种的核心祖先亲本分析

南京农业大学学报, 2001,24(2):20-23.

[本文引用: 2]

田清震, 盖钧镒, 喻德跃, .

我国野生大豆与栽培大豆AFLP指纹图谱研究

中国农业科学, 2001,34(5):480-485.

[本文引用: 1]

盖钧镒.

植物种质群体遗传结构改变的测度

植物遗传资源学报, 2005,6(1):1-8,14.

[本文引用: 1]

韩天富, 盖钧镒.

世界菜用大豆生产、贸易和研究的进展

大豆科学, 2002,21(4):278-284.

[本文引用: 1]

/