基于SNP标记技术的糯玉米种质遗传多样性分析
2021-12-22
来源:乌哈旅游
华北农学报・20t 5,30(3):77・82 基于SNP标记技术的糯玉米种质遗传多样性分析 史亚兴,卢柏山,宋 伟,徐 丽,赵久然 (北京市农林科学院玉米研究中心,北京100097) 摘要:为揭示糯玉米种质的遗传多样性,利用新一代分子标记技术一SNPs(Single nucleotide polymorphisms)对39份 不同基因型的糯玉米自交系进行了基因型分析。结果表明,1 059个SNP标记在39份自交系中的多态性信息含量 (PIC)为0.05~0.38,平均含量为0.31;最小等位基因频率(MAF)为0.03—0.50,平均值为0.29;期望杂合度变化范 围为0.05—0.50,平均期望杂合度为0.39。39份自交系之间的亲缘关系(Kinship)系数为0.00~0.91,京6与京糯6 之间亲缘关系最近。通过Neighbor-joining(NJ)聚类分析,将39份自交系划分为5大类群,其中系谱来源不清晰的糯 玉米自交系均被划分至不同类群中,明确了其亲缘关系。因此,SNP标记适用于糯玉米自交系的遗传多样性分析及 亲缘关系研究,可为糯玉米种质资源利用及品种选育提供参考。 关键词:糯玉米;SNP标记;遗传多样性;京糯6;白糯6 中图分类号:S5l3.03 文献标识码:A 文章编号:1000—7091(2015)O3—0077一O6 dol:10.7668/hbnxb.2015.03.O15 Genetic Diversity Analysis of Waxy Corn Inbred Lines by Single Nucleotide Polymorphism(SNP)Markers SHI Ya—xing,LU Bai—shan,SONG Wei,XU Li,ZHAO Jiu—ran (Maize Research Center,Beijing Academy of Agriculture and Forestry Sciences,Beijing 1 00097,China) Abstract:In order to explore genetic diversity of waxy corn varieties,a novel molecular marker technology—Sin- gle nucleotide polymorphism(SNP)markers were adopted in genotypic analysis of 39 elite waxy corn inbred lines. Polymorphic information content(PIC)of 1 059 SNP markers on these inbred lines were 0.05—0.3 8,with average value being 0.31;Minor allele frequency(MAF)were 0.03—0.50,with average value being 0.29;Expected het- erozygosity were 0.05—0.50,with an average of 0.39.Kinship coeficient between 39 inbred lines ranged from f0.00 to 0.91.Jing 6 and Jingnuo 6 had the strongest similarity.Five groups were divided from the 39 lines by Neighbor-joining(NJ)clustering method.Inbred lines with unclear pedigree were divided into different groups,SO their relatedness with other lines were defined.Therefore,SNP markers were applicable to the analysis of genetic diversity and genetic relationships of waxy corn inbred lines,could provide data for utilization and developing new varieties in breeding program. Key words:Waxy corn;SNP marker;Genetic diversity;Jingnuo 6;Bainuo 6 糯玉米(Zea mays L.certaina Kulesh)是普通玉 米发生基因隐性突变(wx突变为WX)而形成的一 种特殊玉米类型。糯玉米胚乳中淀粉主要为支链淀 的类型,因此在我国遗传资源十分丰富 I4 。近几 十年来,育种家们通过多种育种方法和技术,创新培 育了大量优良糯玉米种质,拓宽了我国糯玉米种质 粉,分子量小,食用消化率高,其所含蛋白质比普通 玉米高3%一6%,赖氨酸、色氨酸含量均较高,是一 种高营养价值和高经济价值的特用玉米¨I2 。糯玉 米起源于中国,也是玉米各类型中唯一起源于我国 遗传基础。因此,分析现有糯玉米种质资源的遗传 多样性,研究种质之间的亲缘关系,对选育优良糯玉 米品种有重要作用。 近年来,分子标记技术在玉米种质遗传多样性 收稿日期:2014一l2—25 作者简介:史亚兴(1977一),男,北京人,副研究员,硕士,主要从事鲜食玉米遗传育种研究。 通讯作者:赵久然(1962一),男,北京人,研究员,博士,主要从事玉米栽培与遗传育种研究。 -C T^一 HRICULTURAE 78 BOilEllLi-SiHICl1.....—— 华北农学报 30卷 分析中已得到广泛应用。RFLP标记是最早发展并 被应用的一种标记技术 ,随后RAPD和AFLP标 记也相继发展和应用 。20世纪80年代末,SSR 记对154份普通玉米自交系进行了遗传多样性分 析;wu等 利用1 015个在基因组中平均分布的 SNP标记,对367份普通玉米自交系进行了群体遗 传多样性的分析并进行聚类;吴金凤等 利用 标记迅速发展起来,是目前玉米育种研究中应用较 多的分子标记技术 ,但SSR标记在操作时需要 利用凝胶电泳对基因位点分型,较费时费力,难以实 现大规模、高通量的分析… 。 1996年Lander¨ 首次提出SNP(Single nucleo. 1 041个SNP标记对51份玉米自交系进行了遗传多 样性分析及类群划分,划分结果与系谱来源一致。 但目前利用SNP分子标记技术分析特用玉米种质 遗传多样性的研究还未见报道。因此,本研究选用 39份我国常用及自选糯玉米自交系,利用覆盖全基 tide polymorphism,单核苷酸多态性)为最新一代分 子标记。SNP是广泛存在于基因组中的一类DNA 序列变异,其变异频率为1%或更高。它是由单个 碱基的转换或颠倒引起的点突变,遗传稳定性更高、 可靠性更强,且摆脱了凝胶电泳的瓶颈,通过序列比 因组的SNP标记进行遗传多样性分析和亲缘关系 研究,旨在为糯玉米优良品种选育、杂优模式开发与 利用等提供理论参考。 对来分析遗传物质差异,易于实现高技术含量、大规 模的基因型检测 。 随着现代分子生物学及基因芯片的发展,SNP 分子标记因其优点受到国内外育种家越来越多的重 视,并在黄瓜 、苹果… 、大豆 、小麦 、水 稻¨ 、玉米 等作物的品种鉴定、遗传作图、基因 1 材料和方法 1.1试验材料 供试材料分为两部分;37份常用及自选糯玉米 自交系,编号为N1~N37;2份普通玉米自交系 B73,Mo17,作为标准测验种,分别代表Reid群和 Lancaster群,编号为N38~N39。各自交系名称及 系谱来源见表1。 挖掘与QTL定位等方面均有应用。在玉米种质遗 传多样性研究方面,Yang等¨ 利用884个SNP标 表1 39份自交系名称及系谱来源 Tab.1 List of 39 inbred lines and their pedigrees 编号 No. 名称 Name 系谱来源 Pedigree 编号 No. 名称 Name 系谱来源 Pedigree N1 N2 N3 N4 N5 N6 N7 N8 京糯6 白糯6 京6 BN2 N601 N39 紫糯264 N8131 选自中糯1号 选自紫糯3号 选自中糯1号 选自紫糯3号 京糯968母本 选自中糯309 不祥 不祥 N21 N22 N23 N24 N25 N26 N27 N28 紫糯419 紫糯680 黄糯705 黑糯732 紫糯961 黄糯850 紫糯736 京糯5 不详 不详 不详 不详 不详 不详 不详 选自黄糯种质 N9 N10 N11 N12 N13 京糯72 紫糯5 紫糯3 紫玉2 Z09 昌7—2转糯选系 不详 选自紫糯3号 选自紫香糯1号 引自美国DH系 N29 N30 N31 N32 N33 黄糯6 京糯2 中玉04 糯1 糯2 选自紫糯3号 选自白早糯种质 中糯1号母本 垦粘1号母本 垦粘1号父本 N14 N15 N16 N17 N18 N19 N20 糯70 京糯255 CB1 紫糯258 紫糯263 紫糯877 紫糯196 昌7-2转糯选系 不详 选自美国双隐性种质 不详 不详 不详 不详 N34 N35 N36 N37 N38 N39 京甜糯2 D6644 京甜糯69 D6644—6 B73 Mol7 选自SH一251×白糯6 选自SH一251×白糯6 选自SH一251×白糯6 选自SH一251×白糯6 BSSSC5 l87—2×1O3 1.2试验方法 养发芽,玉米幼苗长至4~5片叶时,取鲜嫩玉米叶 将供试种子在恒温培养箱内培 片,参照Rogers和Bendich 的方法提取叶片基因 1.2.1 DNA提取l C T ● 3期 史亚兴等:基于SNP标记技术的糯玉米种质遗传多样性分析 iORIOULTURA| 79 l明征-tl—SIHICll 赫 童 ∞ ∞∞∞组DNA。所提DNA要求条带单一,完整无降解;紫 1.2.4 数据统计分析 本研究数据利用Power— 外分光光度计检测A: 。 。为1.8~2.0(DNA样品没 Marker V3.25软件 计算其最小等位基因频率 有蛋白、RNA污染),A: 。 。为1.8~2.0(DNA样品 (MAF)、多态性信息含量(PIC)、期望杂合度(Ex. 盐离子浓度低),稀释样品DNA浓度至100 ng/ ̄L。 pected heterozygosity,He)、观测杂合度(Observed 1.2.2 SNP位点选取利用北京市农林科学院玉 heterozygosity,Ho)等。利用SPAGeDi (Spatial 米研究中心开发的Maize SNP 3072芯片对供试自交 pattern analysis of genetic diversity)进行自交系间亲 系进行SNP基因分型 。根据以下指标选用1 059 缘关系(Kinship)分析;采用Nei S 1972算法计算自 个SNP位点:数据质量评估值≥0.7;AA、AB、BB这 交系间遗传距离,并按照邻接法(Neighbor—joining) 3种基因型具备明显界限;位点特异性探针设计评 构建NJ聚类图,进行杂种优势群的划分。 估值≥0.7;在全基因组上均匀分布,用于本研究。 1.2.3基因芯片(Golden gate技术)实验操作流程 2 结果与分析 将提取的玉米基因组DNA与激活的生物素充分 2.1 SNP多态性分析 结合,通过离心沉淀使DNA纯化;将探针与纯化的 利用在全基因组中平均分布、数据质量评估值 目的DNA链杂交,漂洗杂交产物后进行延伸连接反 高的1 059个SNP标记,并通过基因芯片的方法检 应;延伸的产物经漂洗变性处理后作为模板进行 测了39份糯玉米自交系的基因型,并用PowerMark- PCR扩增,PCR产物与磁珠结合后过滤纯化;将纯 er V3.25对基因分型数据进行分析。结果表明, 化的PCR产物与芯片杂交,杂交结束后进行芯片清 1 059个SNP标记在39份自交系中共检测到2 118 洗、真空抽干;风干后的芯片立刻在iScan仪上进行 个等位基因,缺失率为0~21%,平均为4%。平均 扫描,最后利用Genome Studio软件进行数据分析 期望杂合度为0.39,范围为0.05~0.50;平均观测 (详细实验流程参照Illumina公司提供的实验操作 杂合度为0.04,范围为0.00~0.29;平均多态性信 指南Golden Gate ̄Genetyping Assay,Manual Part}} 息含量为0.31,变化范围为0.05~0.38(表2)。 15004070 Rev.B)。 表2 1 059个SNP标记在39份自交系中的多态性 Tab.2 Polymorphisms of 1 059 SNPs in 39 inbred lines 糯玉米自交系上的平均PIC值为0.31,说明所用 SNP标记为中度偏高多态性。期望杂合度常被用来 衡量群体的遗传多样性,平均期望杂合度越高,表明 群体遗传一致性越低,遗传多样性越高 。由表2 可知,本研究中平均期望杂合度为0.39,表明39份 自交系有较好的遗传多态性。 最小等位基因频率(MAF)也叫作次要等位基 o.1O o.14 o.18 0.22 0.26 0.30 0.34 o.38 0.42 o.46 o.50 因频率。本研究中,1 059个SNP标记的MAF值为 最小等位基因频率 IgAF 0.03~0.50,平均值为0.29,其中≥0.20的位点有 835个,占79%。由图1可以看出,MAF在0.10~ 图1 1 059个SNP标记在39份 0.50之间均匀分布。 自交系中的最小等位基因频率 2.2不同供试基因型间的亲缘关系(Kinship)分析 Fig.1 MAF of1 059 SNPs in 39 inbred lines 39份自交系的Kinship分析表明,98.8%的自 PIC值是衡量DNA座位变异程度高低的重要 交系间亲缘关系系数(Kinship系数)分布在0.00~ 指标。某DNA座位的PIC值>0.50时,表明其为 0.50之间。其中60%的Kinship系数为0,表明 高度多态性座位;O.25<PIC值<0.50时,表明其为 60%的成对自交系间无亲缘关系。23.8%的Kin. 中度多态性座位;PIC值<0.25时,表明为有低度多 ship系数分布在0.00~0.10之间,说明这些成对自 态性座位 。本研究中,1 059个SNP标记在39份 交系之间的亲缘关系较弱。15%的Kinship系数分 华《 丑 北农学报 30卷 加 ∞ ∞ ∞ ∞ ∞ 加 0 布在0.10~0.50之间,表明这些自交系之间有不同 2.3不同供试基因型的聚类分析 程度的亲缘关系。另有1.2%的Kinship系数分布 根据39份自交系在1 059个SNP位点的基因 在0.50—0.91之间,说明这些成对自交系间有很大 的亲缘相关性 。39份自交系之间,京6与京糯6 之间Kinship系数最大为0.91,白糯6与BN2之间 Kinship系数次之,为0.82,说明这2对自交系之间 分型数据,运用PowerMarker V3.25软件,并应用 Nei S 1972算法计算自交系之间的遗传距离,构建 NJ聚类图,进行聚类分析(图3)。结果表明,39份 自交系间的遗传距离变化为0.00—0.72,自交系间 平均遗传距离为0.47,其中,自交系白糯6与BN2 亲缘关系最近。此结果与系谱来源一致,也表明我 国糯玉米种质资源有较广泛的遗传基础。 之间遗传距离最小,B73与Mo17之间遗传距离最 远,与Kinship分析结果一致。 聚类分析结果表明,39份自交系可被划分为5 大类群(表3)。第一类群以普通玉米自交系B73为 代表,包括B73以及13个糯玉米自交系,称为Reid 群;第二类群包括9个糯玉米自交系,称为京糯6 群;第三类群包括2个糯玉米自交系,称为糯1群; 第四类群包括7个糯玉米自交系,称为糯2群;第五 0 0.O5 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 类群是以Mo17为代表的Lancaster群,包括Mo17 Kinship系数 Kinship coefficient 以及6个糯玉米自交系(表3)。通过聚类分析和 Kinship分析,将系谱不清晰的糯玉米自交系划分至 不同的类群中,明确了与其余自交系间的亲缘关系, 图2 39份自交系之间的Kinship系数分布 Fig.2 Distribution of Kinship coeficientf of 39 lines 可提高其利用效率。 表3 39份自交系类群划分明细 Tab.3 List of 39 lines of diferent heterotic groups 3 结论与讨论 3.1 SNP遗传多态性 PIC值变化范围为0.05~0.38,平均值为0.31,平 均MAF为0.29。与前人研究结果相比较,本研究 所用SNP标记有较好多态性,属于中度偏高多态性 SNPs是指基因组DNA中由于单个核苷酸的替 换而引起的多态性。以这样的DNA序列多态性特 征做标记即可区分2个个体之间遗传物质的差异。 SNP标记广泛存在于基因组中,多态性数量丰富,与 座位,同时表明糯玉米自交系有较好的遗传多样性。 Yang等¨ 利用884个SNP标记对154份普通玉米 自交系进行研究,所得平均期望杂合度为0.44。本 研究中期望杂合度变化为0.05~0.50,平均值为 0.39,低于Yang等的结果。这可能是由于本研究所 选糯玉米自交系数量较少的原因。 3.2糯玉米种质的聚类分析 新品系的鉴定和筛选是作物育种中的重要环 SSR标记相比,SNP标记稳定性更好,与基因芯片技 术结合,可快速、高通量的检测个体遗传多样性 。 Yang等 利用在全基因组分布的926个SNP 标记分析了527份普通玉米自交系的遗传多样性以 及亲缘关系,结果表明,926个SNP标记的平均PIC 节,利用聚类分析将不同种质资源进行分类,可为品 系(种)间的遗传差异比较、杂交育种亲本组配提供 参考依据。Wu等 叫利用1 015个SNP标记将367 份普通玉米自交系分为2个大群和5个亚群,结果 值为0.31,平均MAF为0.3;Lu等 利用1 034个 SNP标记对770份国内外玉米自交系进行遗传多样 性分析,多态性信息含量平均值为0.26。本研究 中,1 059个SNP标记在39份糯玉米自交系中的 与系谱一致。在糯玉米种质类群划分研究上,前人 3期 史亚兴等:基于SNP标记技术的糯玉米种质遗传多样性分析 白黄紫京 j晕i晕j晕 舭眺树黼瑚啪吼黼躯 蕊锨鲫 紫京京紫中紫京京糯京黑黄黄鼢京糯糯紫 紫紫紫 糯糯糯糯 糯神0●糯糯糯糯∞糯 糯"糯糯幡 3 徽 陀 lNB131 紫 I264 U.趵0. 0.1a u.1U U.∞ UUU 圈3 基于SNP标记的39份自交系NJ聚类图 Fig.3 Neighbor-joining(NJ)trees for 39 inbred lines based on SNP markers 利用SSR标记做了较多研究。王慧等 利用65对 SSR引物将165份糯玉米材料划分为3大类,划分 结果基本符合材料来源;程宇坤等 利用49对 SSR引物将58份糯玉米自交系划分为5个类群。 本研究中,利用覆盖玉米全基因组的1 059个SNP 标记,通过Neighbor—joining(NJ)聚类分析,将39份 自交系划分为5个类群,划分结果与系谱来源基本 一致,说明利用SNP分子标记可用于糯玉米种质遗 传多样性分析及类群划分,并能取得良好效果。 本研究类群划分结果表明,我国优良糯玉米品 种垦粘1号的父母本糯2和糯1被划分至不同类群 中,因此,糯1群×糯2群的杂种优势模式可供借 鉴;国审糯玉米品种京科糯2000(京糯6×BN2)、京 紫糯218(紫糯5×紫糯3)的双亲均分别来自京糯6 群和Reid群,说明京糯6群×Reid群也是我国糯玉 米品种选育的主要杂优模式,在亲本选配中可进行 利用。本研究中系谱来源不清晰的糯玉米自交系被 划分至不同类群中,依据此结果可指导杂交种组配, 避免盲目测配,提高了其利用效率。 参考文献: [1] 梁龙金,董晓君,董长江.糯玉米的特点及其栽培技术 [J].现代化农业,2002(8):16. [2] 罗高玲,吴子恺.甜玉米和糯玉米遗传基因的研究进 展[J].广西农业科学,2003(3):24—26. [3] 谢孝颐,蔡志飞,印志同,等.糯玉米育种概述[J].玉 米科学,2003(专刊):58—67. [4] 黄玉碧,荣廷昭.我国糯玉米种质资源的遗传多样性 和起源进化[J].作物杂志,1998(s1):77—80. [5] 胡丹东,赵久然.DNA分子标记技术及其在玉米育种 中的应用[J].甘肃农业大学学报,2007,12(6):92— 98. [6] 刘新芝,彭泽斌,傅骏骅,等.RAPD在玉米类群划分研 5 m究中的应用[J].中国农业科学,1997, 量詈季}啪 30(3):44—51. [7] 吴敏生,王守才,戴景瑞.AFLP分子标记在玉米优良 自交系优势群划分中的应用[J].作物学报,2001,26 (1):9—13. [8] Smith J S C,Chin E C L,Shu H,et a1.An evolution ofthe utility of SSR loci as molecular markers in maize(Zea mays L.):Comparison with data from RFLPs and pedi— gree[J 3.Theor Appl Genet,1997,95:163—173. [9] 王鹏文,杨 勇.利用SSR标记划分糯玉米的杂种优 势群[J].华北农学报,2007,22(3):35—38. [10] Yan J B,Yang X H,Hector S,et a1.High—throughput SNP genotyping with the golden gate assay in maize[J]. Mol Breed。2010。25:441—451. [11] 唐立群,肖层林,王伟平.SNP分子标记的研究及其应 用进展[J].中国农学通报,2012,28(12):154—158. [12] Lander E S.The new genomics:global views of biology [J].Science,1996,274:536. [13] 李欧静,张桂华,兰青阔,等.基于SNP标记的种子纯 度高效检测分析模型的建立[J].湖南农业科学, 2012(19):9一l1. [14] Han Y P,David C.BAC—end sequence—based SNPs and Bin mapping for rapid integration of physical and genetic maps in apple[J].Genomics,2009(93):282—288. [15] Lam H M,Xu X,Liu X,et a1.Resequencing of 31 wild and cultivated soybean genomes identifies patterns of ge— netie diversity and selection[J].Nature Genetics,2010 (42):1053. [16] Eduard Akhunov,Charles Nicolet,Jan Dvorak,et a1.Sin. gle nucleotide polymorphism genotyping in polyploid wheat with the lllumina GoldenGate assay『J].Theor Appl Genet,2009,119(3):507—517. [17] Yu H H,Xie W B,Wang J,et a1.Gains in QTL detection using an ultra-high density SNP map based on population sequencing relative to traditional RFLP/SSR markers [J].PLoS ONE,2011,6(3):e17595. [18] Kristen L Kump.Peter J Bradbury.Genome.wide associ. ation study of quantitative resistance to southern leaf blight in the maize nested association mapping popula— tion[J].Nature Genetic,2011,43(2):163—168. [19] Xiaohong Yang,Yunbi Xu,Trushar Shah,et a1.Compari. 詈啪I C T ● llORICULTURItE B0睢lltl-SIHI ̄ll 82 华北农学报 30卷 son of SSR and SNPs in assessment of genetic related. length polymorphisms[J].Am J Hum Genet,1980(32) 314—331. ness in maize[J].Genetica,2011,139:1045—1054. [20] Xun Wu,Yongxiang Li,Yunsu Shi,et a1.Fine genetic characterization of elite maize germplasm using hight- [26] Frankham R,Ballou J D,Briscoe D A.Introduction tO conservation genetics[M].Cambridge:Cambridge Uni— versity Press,2002:29—62. hroughput SNP genotyping[J].Theor Appl Genet, 2014,127:624—631. [27] Xiaohong Yang,Shibin Gao,Shutu Xu,et a1.Character— [21] 吴金凤,宋 伟,王 蕊,等.利用SNP标记对5l份 ization of a global germplasm collection and its potential [22] [23] [24] [25] 玉米自交系进行类群划分[J].玉米科学,2014,22 (5):29—34. Rogers S O,Bendich A J.Extraction of DNA from milli— gram amounts of flesh,herbarium and mummified plant tissues[J].Plant Mol Biol,1985(5):69—76. Liu K J,Muse S V.PowerMarker:an integrated analysis environment for genetic marker analysis[J].Bioinformat— ics,2005,21(9):2128—2129. Hardy,Vekemans.SPAGeDi:a versatile computer pro— gram to analyse spatial genetic structure at the individu— al or population levels[J].Molecular Ecology Notes, 2002(2):618—620. Botstein D,White R L,Skolnick M,et a1.Construction of a genetic linkagemap in man using restriction fragment utilization for analysis of complex quantitative traits in maize[J].Mol Breeding,2010,28:511—526. [28] 邹喻苹,葛 颂.新一代分子标记一SNPs及其应用 [J].生物多样性,2003,11(5):370—382. [29] Lu Y,Yan J,Guimaraes C T,et a1.Molecular character ization of global maize breeding germplasm based on ge— nome—wide single nucleotide polymorphisms[J].Theor Appl Genet,2009,120(1):93—115. [30] 王慧,卢有林,孙大鹏,等.糯玉米种质品质性状鉴 定和SSR标记遗传多样性分析[J].植物遗传资源学 报,2013,14(5):800—808. [31] 程宇坤,白建荣,张效梅,等.山西省糯玉米自交系的 遗传多样性分析及类群划分[J].山西农业科学, 2012,40(5):433—438.