作者简介:王国荣(1988—),男,云南玉溪人,硕士研究生,主要从事植物分子育种研究。E-mail: 294828485@qq.com
为揭示裸大麦育成品种与种质资源材料的遗传结构,利用反转录转座子标记(反转录转座子-微卫星扩增多态性,REMAP;反向反转录转座子扩增多态性,IRAP)与简单重复序列(SSR)标记,分析63份裸大麦的遗传多样性。IRAP、REMAP、SSR标记分别检测到315、143、38个等位变化,变异范围分别为9~58、7~25、2~4个,平均每对引物检测24.23、14.30、2.38个等位变化。REMAP和IRAP单一标记多态性位点贡献率高于SSR标记。反转录转座子标记揭示材料间遗传相似性系数(GS)为0.452~0.937,平均为0.674,在GS值0.620水平上将63份材料分为2大类,分别包含20份和43份材料。SSR标记结果显示材料间GS为0.351~0.973,平均为0.716,在GS值0.620水平上将63份材料区分为2大类,分别包含2份和61份材料。反转录转座子标记聚类结果在GS值0.740水平上能区分西藏野生资源和栽培资源,但SSR标记不能有效区分这2类材料。反转录转座子标记的主坐标及群体结构分析结果与GS聚类分析结果一致性较高;SSR标记群体结构与主坐标分析、GS聚类分析结果存在明显差异。本研究表明,裸大麦材料遗传背景简单,亚类间缺少基因交流,反转录转座子标记在裸大麦品种亲缘关系鉴定中有一定的优势。反转录转座子标记与SSR标记的协同使用可为裸大麦遗传多样性分析、种质鉴定及亲本选配揭示更多有用信息。
To evaluate the genetic structure of hulless barley, retrotransposons-based markers (retrotransposon-microsatellite amplified polymorphism, REMAP; inter-retrotransposon amplified polymorphism, IRAP) and simple sequence repeat (SSR) markers were used to obtain the genetic data of 63 hulless barley materials. The IRAP, REMAP and SSR markers detected 315, 143 and 38 alleles, ranged from 9-58, 7-25, and 2-4, with an average of 24.23, 14.30, 2.38 alleles per marker, respectively. The contribution rate of single retrotransposons-based marker polymorphic site was higher than that of SSR marker. As revealed by retrotransposons-based markers, the genetic similarity (GS) among these 63 accessions was from 0.452 to 0.937, with the average of 0.674. The accessions were divided into two subgroups at the GS level of 0.620, which contained 20 and 43 materials, respectively. As revealed by SSR markers, GS among these 63 accessions was from 0.351 to 0.973, with the average of 0.716. The accessions were divided into two subgroups at the GS level of 0.620, which contained 2 and 61 materials, respectively. The wild hulless barley accessions was clustered into a small category by the retrotransposons-based markers at the GS level of 0.740. However, it was unable to distinguish the wild hulless barley from the cultured by SSR markers in the present study. The population structure revealed by retrotransposons-based markers was consistent with the principal coordinate analysis and genetic similarity analysis. However, the population structure revealed by SSR markers showed difference. It was shown that most of the accessions were with single genetic components. The retrotransposons-based markers had advantage in the identification of the relationship of hulless barleys. And more genetic information could be explored by combining the retrotransposons-based markers and SSR markers.
裸大麦(青稞)是青藏高原地区最重要的粮食作物, 但其育成品种与种质资源材料的遗传多样性水平存在争议。孟凡磊等[1]、杨平等[2]的研究认为, 川藏地区育成青稞品种的遗传多样性较低。曾兴权等[3]的研究表明, 同一来源地区的青稞育成品种遗传基础较为狭窄, 但不同地区间的遗传差异较大。Feng等[4]、赖勇等[5]和巴桑玉珍等[6]利用简单重复序列(simple sequence repeat, SSR)标记对青藏地区青稞品种与资源材料的群体遗传结构、地理分化及与重要农艺性状的关联分析进行研究, 认为青稞育成品种遗传基础较为狭窄, 但种质资源的遗传多样性较高。
反转录转座子是高等植物基因组中广泛存在的一类可移动的DNA元件[7], 是植物基因组进化的重要动力[8, 9, 10], 参与植物性染色体形成[10]。反转录转座子是大麦基因组的重要组成部分[11], 是大麦基因组响应外界胁迫的重要元件[12, 13], 并能介导基因缺失[13], 在大麦基因组进化中具有重要作用[13, 14]。反转录转座子标记具有灵敏度高、基因组覆盖度广、多态性高等特点, 适用于遗传多样性研究及种质资源鉴定[15, 16, 17]。利用分子标记进行遗传多样性分析、指纹图谱构建、关联分析研究的报道在动植物中都较为常见[18, 19, 20], 然而利用反向反转录转座子扩增多态性(inter-retrotransposon amplified polymorphism, IRAP)、反转录转座子-微卫星扩增多态性(retrotransposon-microsatellite amplified polymorphism, REMAP)等反转录转座子标记进行大麦遗传群体结构研究的报道较少[17], 至今未见利用IRAP、REMAP揭示裸大麦育成品系及资源材料遗传多样性和遗传结构的报道。
鉴于此, 本研究选取63份来自不同地区的裸大麦品系及西藏裸大麦资源材料, 利用10对REMAP标记、13对IRAP标记对这些材料的遗传多样性及群体遗传结构进行研究, 并与SSR标记揭示的遗传多样及群体结构进行比较分析, 从反转录转座子角度丰富裸大麦的遗传多样性研究, 提供新颖的遗传多样性资料, 以期更好地解析裸大麦的遗传基础, 为裸大麦及其种质资源的高效利用提供更多依据。
供试63份裸大麦材料包含36份青海品系、9份西藏野生资源材料、4份西藏品系、5份云南品系、3份江苏品系、2份国外材料, 以及湖北、甘肃、黑龙江、浙江材料各1份(表1)。供试青海、西藏品系由拉萨市农业科学研究所提供。供试西藏野生材料由武汉大学生命科学学院丁毅教授提供。其他材料由浙江省农业科学院作物与核技术利用研究所国家大麦改良中心提供。
PCR体系中所用10× buffer、Mg2+、dNTPs、TaqDNA聚合酶均购自上海博彩生物科技有限公司, 引物由南京金斯瑞生物科技有限公司合成。
取三叶期新鲜叶片, 用CTAB小样法提取基因组DNA, 琼脂糖凝胶电泳检测所提DNA的质量, NanoDrop 2000c超微量分光光度计测定DNA浓度, -80 ℃保存备用。
IRAP与REMAP标记引物序列见表2。用于IRAP标记的共有13对组合, 分别为A1、A2、A3、A4、A10、A1+A11、A2+A11、A6+A13、A12+A13、A1+A2、A1+A5、A12+A14、A13+A14; 用于REMAP标记的共有10对组合, 分别为A2+A15、A5+A16、A5+A15、A6+A16、A9+A15、A11+A16、A12+A15、A13+A15、A14+A15、A2+A16。
PCR扩增体系(20 μ L):10× buffer 2 μ L, 25 mmol· L-1 Mg2+ 1.6 μ L, 10 mmol· L-1 dNTPs 0.4 μ L, 5 U· μ L-1Taq酶0.3 μ L, 10 mmol· L-1引物各1 μ L, 50 ng· μ L-1模板DNA 2 μ L, ddH2O 11.7 μ L。PCR扩增程序:94 ℃预变性5 min; 94 ℃变性30 s, 54~58 ℃退火40 s, 72 ℃延伸2 min, 36个循环; 72 ℃延伸10 min。
REMAP标记PCR产物使用质量分数为2.5%的琼脂糖凝胶电泳检测, 使用凝胶成像系统(Syngene, G:BOX-CHEMI-XRQ)记录结果; IRAP标记PCR产物使用质量分数为5%的非变性聚丙烯酰胺凝胶电泳分离, 硝酸银染色检测带纹, 数码相机拍照记录结果。
SSR标记引物及其染色体位置信息由澳大利亚莫道克大学李承道教授提供(表3)。SSR扩增体系(10 μ L):10× buffer 1.0 μ L, 25 mmol· L-1 Mg2+ 0.8 μ L, 10 mmol· L-1 dNTPs 0.2 μ L, 5 U· μ L-1Taq DNA聚合酶0.15 μ L, 10 μ mol· L-1正反向引物各0.5 μ L, 50 ng· μ L-1模板DNA 2 μ L, ddH2O 4.85 μ L。SSR标记PCR扩增程序:95 ℃预变性4 min; 94 ℃变性30 s, 54~61 ℃退火50 s, 72 ℃延伸50 s, 33个循环; 72 ℃延伸10 min。SSR标记PCR产物使用质量分数为6%的非变性聚丙烯酰胺凝胶电泳分离, 硝酸银染色检测带纹, 数码相机拍照记录结果。
IRAP、REMAP和SSR标记的结果以二进制和基因型记录, 同一位点上具有相同迁移率的条带记为1, 无带记为0, 同时以数字1、2、3 等代表不同等位基因。通过Popgene 32统计分析计算基因频率、Shannon指数等遗传参数。以NTSYS-PC软件计算遗传相似性系数(genetic similarity coefficient, GS), 按照非加权平均法(unweighted pair-group method with arithmetic means, UPGMA)和SHAN 程序进行聚类分析作图。利用EIGEN程序进行主坐标分析。以Structure 2.3.1软件分析群体遗传结构, 估计最佳群体组群数(K), 取值范围为1~11, 将参数iterations设为10 000次, 再将burn-in period设为100 000次, 重复11次, 计算Q参数。当K值持续增大时, 通过计算Δ K确定K值[21]。
REMAP标记片段大小集中于500~3 000 bp之间(图1-A), 10对标记共检测到143个等位变异, 每标记平均14.30个, 变幅7~25个, 其中有131个多态性位点, 占总条带数的91.61%, 平均多样性指数为0.303 3, 每个位点的有效等位基因数为1.514 6。IRAP标记条带大小集中于500~2 000 bp(图1-B), 13对标记共检测到315个等位变异, 每标记平均24.23个, 变幅9~58个, 其中有303个多态性位点, 占总条带数的96.19%, 平均多样性指数为0.268 2, 每个位点的有效等位基因数为1.440 5。16对SSR标记共检测到38个等位变异(图1-C), 每标记平均2.38个, 条带集中于100~400 bp, 其中有16个多态性位点, 多态性位点百分率 100%, 平均多样性指数为0.342 8, 每个位点的有效等位基因数为1.610 2。
标记类型不同, 单标记所得条带数及有效变异数差异明显。SSR标记单标记平均贡献2.38条条带, 平均有效变异数为2.38; REMAP标记单标记平均贡献14.30条条带, 平均有效变异数为13.10; IRAP标记单标记平均贡献24.23条条带, 平均有效变异数为23.31。较之SSR标记, 反转录转座子标记的单一标记多态性位点贡献率较高。
反转录转座子标记聚类结果显示, 63份裸大麦材料GS变幅为0.452~0.937, 平均值0.674(图2-A)。其中:紫光芒与紫光芒166间的GS值最大(0.931), 表明二者亲缘关系最近; 青海21与野生材料80322、80402的GS值最小(均为0.452), 表明其亲缘关系较远(图2-A)。在GS值0.620水平上, 63份祼大麦材料分为2大类。第1大类可进一步分为2个亚类, 第1亚类由鄂507、苏裸2号及Inbyf-08-11组成, 第2亚类中包含8份栽培材料和9份野生材料。在GS值0.740水平上, 第2亚类中的9份野生材料单独聚为一个小类(图2-A), 说明反转录转座子标记可以很好地区分野生裸大麦材料与其他材料。第2大类在GS值0.680水平上以青海1为一个亚类, 剩余42份材料为另一亚类。
反转录转座子标记主坐标分析结果显示(图3-A), 反转录转座子标记将供试材料分为2大类。第1大类共20份材料包含3个亚类:第1个亚类包含鄂507、苏裸2号、Inbyf-08-11; 第2个亚类是9份野生材料; 第3个亚类为其余8份材料。第2大类共43份材料, 由青海1和其余42份材料组成。主坐标分析结果与聚类分析结果一致度较高, 材料所属类群界限明显。
SSR标记聚类结果(图2-B)显示, 63份裸大麦材料GS变幅为0.351~0.973, 平均值0.716。其中:门农1与青海7的GS值最大(1.000), 表明二者亲缘关系较近, 其次为北青3与青海7、北青3与门农1, GS值均达到0.973; 青海6与苏裸2号的GS值最小(0.351), 表明二者亲缘关系较远。聚类分析表明, 在GS值0.620水平上63份裸大麦材料可分为2大类, 其中苏裸1号、苏裸2号为第1类, 其余的61份材料为第2类。第2类在GS值0.660水平上可分为2个亚类, 大部分野生裸大麦(8份)包括在第1亚类的21份材料中, 裸大麦1579和青海4等40份材料为第2亚类(图2-B)。SSR标记聚类结果中, 野生裸大麦材料没有被单独聚为一类。
SSR标记主坐标分析结果(图3-B)显示, SSR标记将供试材料分为2大类:第1大类包括黑麦10号等22份材料; 第2大类包含2个亚类, 即由紫光芒等21份材料组成的第1亚类和由青海8等20份材料组成的第2亚类。主坐标分析中第1大类材料与GS聚类区分的第1大类及第2大类下的第1亚类(紫光芒除外)所属材料一致, 其他GS区分的第2类(紫光芒除外)所属材料同主坐标分析的第2类材料一致。
基于反转录转座子标记的数据群体结构分析表明, 样本的等位变异频率特征类型数K呈持续增大趋势, 当K=2时, Δ K值最大, 据此将63份材料分成2个亚群(图4-A), 分别包含20、43份材料, 其中60份材料Q值大于0.7, 表明这些裸大麦材料来源单一, 遗传背景简单, 亚类间缺少基因交流(图5-A)。同反转录转座子标记GS聚类及主坐标分类2个类群的20、43份材料相比, 分类一致。
基于SSR标记的数据群体结构分析显示, 样本的等位变异频率特征类型数K呈持续增大趋势, 当K=2时, Δ K值最大, 据此将63份材料分成2个亚群(图4-B), 分别包括23、40份材料, 其中58份材料Q值大于0.7。这再次证明了这些裸大麦材料来源单一, 遗传背景简单, 亚类间缺少基因交流(图5-B)。但与主坐标分析、聚类分析所区分的材料大类、亚类所属类群相比, 无明显重合性。
反转录转座子标记是鉴定和区分种质资源材料的有效方法。基于IRAP与REMAP标记可以区分葡萄牙历史上育成的所有小麦品系[15]。Biswas等[22]在柑橘中的研究也表明, 反转录转座子标记是区分柑橘及其近缘种的有效手段。Kim等[23]建立了基于反转录转座子标记的日本梨品质识别方法。杜晓云等[24]利用IRAP 建立柿属植物指纹图谱, 能有效区分芽变品种。肖炳光等[25]用 IRAP标记的方法对烤烟品种扩增, 平均多态性比率达96%, 显著优于其他测试标记分析结果。
在本研究中, IRAP、REMAP、SSR标记分别检测到315、143、38个等位变化; 变异范围分别为9~58、7~25、2~4个, 平均每对引物检测到24.23、14.30、2.38个。赖勇等[5, 26]利用SSR标记对国内外大麦材料进行检测, 检测出的等位基因变异范围为2~7个, 平均每对标记检测2.9个。同这些研究相比, 反转录转座子标记(IRAP、REMAP)检测到的等位基因变异范围更大, 平均每个标记检测到变异数更多, 说明反转录转座子水平上的遗传多样性分析能揭示更为丰富的裸大麦遗传结构信息, 可为后续育种利用提供更多信息。
赖勇等[5]用SSR标记检测青稞材料的GS变幅为0.497~0.970, 平均0.761; 巴桑玉珍等[6]利用SSR标记分析青稞耐寒资源的GS变化范围为0.469~0.924, 平均0.745; 尚毅等[27]的研究表明, 浙江省裸大麦地方品种遗传多样性水平较高。较以上研究, 本研究揭示的裸大麦遗传相似性系数变异幅度更大, GS平均值更小, 亲缘关系更远, 这可能与本研究中试验材料地理来源丰富, 并且使用了反转录转座子及SSR两种标记有关。
基于数学模型的群体结构分析中, 本研究认为这些裸大麦材料来源单一, 遗传背景简单, 亚类间几乎没有基因交流, 在育种利用的过程中, 可以考虑加大不同类群裸大麦材料之间的交流。本研究SSR标记的群体结构分析与GS聚类结果差异较大, 这可能是因为群体结构分析可以降低人为主观分类对关联分析的影响, 能更准确地了解这些材料的遗传背景[28]。群体结构分析是关联作图的基础, 本研究基于2种标记的群体结构分析结果存在较大差异。在后续关联分析中, 应充分考虑标记类型对关联分析结果的影响。
本研究中, 反转录转座子与SSR标记的GS聚类都将63份裸大麦划分为2个大的类群, 但反转录转座子标记的GS聚类及主坐标分析皆能将野生祼大麦单独聚为一个小类, 而SSR标记聚类结果中野生裸大麦材料未能单独聚为一类。这种差异可能是由不同标记揭示了基因组不同类型和DNA多样性所造成的。IRAP和REMAP反映基因组反转录转座子区域的结构变异, SSR标记揭示的是基因组简单重复序列的多样性。反转录转座子在植物非生物胁迫应答及基因组进化中起着重要作用[29]。已有充分的证据表明, 反转录转座子驱动了大麦基因组进化, 是大麦基因组响应外界胁迫的重要元件[12, 13, 14]。西藏大麦材料基因组反转录转座子区域出现了不同于其他材料的进化模式是个值得进一步研究的问题。
本研究表明, 多标记协同检测可以揭示材料间更丰富的遗传背景差异。但本研究未能对不同标记间检测结果差异化的存在做更深入的探究。基于反转录转座子的标记较常规的SSR标记, 提供的信息丰富, 用时少, 费用低[30], 而且具有检测区域不同、灵敏度高、多态性好、基因组覆盖度广等诸多优势。对以上2种标记协同区分裸大麦遗传背景差异化的研究, 将有利于发掘常规标记检测局限区域外有差异化的亲本材料, 对多标记协同在遗传多样性分析、种质鉴定、亲缘关系分析及亲本选配中的应用提供了实例。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|