引用本文
刘新春, 赖运平, 余毅, 袁金娥, 巴桑玉珍, 冯宗云. 青稞籽粒饲草性状遗传特征. 草业科学, 2018,35(6): 1425-1434
Liu Xin-chun, Lai Yun-ping, Yu Yi, Yuan Jin-e, Basangyuzhen, Feng Zong-yun. Genetic dissection of grain forage traits in hulless barley. Pratacultural Science, 2018,35(6): 1425-1434.  
Doi: 10.11829/j.issn.1001-0629.2017-0409
Permissions
Copyright©2018, 《草业科学》编辑部
《草业科学》编辑部
青稞籽粒饲草性状遗传特征
刘新春1, 赖运平1,2, 余毅3, 袁金娥1,2, 巴桑玉珍1,4, 冯宗云1
1.四川农业大学农学院作物遗传育种学系大麦青稞研究中心,四川 成都,611130
2.成都农业科技职业学院,四川 成都,611130
3.四川省农业科学院作物研究所 农业部植物新品种测试(成都)分中心,四川 成都,610066
4.西藏自治区农牧科学院农业研究所,西藏 拉萨,850032
通讯作者:冯宗云(1963-),男,四川广安人,博士,教授,主要从事大麦遗传育种与分子生物学研究。E-mail:[email protected]

第一作者:刘新春(1985-),男,四川乐山人,在读博士生,主要从事青稞品质性状分子生物学研究。E-mail:[email protected]

收稿日期: 2017-08-13 接受日期: 2017-10-11
资助项目: 国家现代农业产业技术体系(大麦青稞体系)建设专项基金(CARS-05)
摘要

为了解青稞( Hordeum vulgare var. nudum Hook. f.)的饲草性状的遗传特征并选育具有优质饲草品质和高产的青稞新品种,对60份来源于中国青藏高原的青稞籽粒饲草性状进行研究,利用TASSEL软件的混合线性模型对其进行非条件和条件关联分析,同时结合分离群体群分离法对影响籽粒饲草性状的分子标记进行检测。结果表明,青稞品种间存在丰富的籽粒饲草性状变异;群体结构分析显示该群体可以分成3个亚群,群体的分类结果与品种的来源地有密切的相关性;在非条件关联分析中共发现23个SRAP分子标记和籽粒饲草品质性状(包括粗蛋白含量、淀粉含量、灰分含量)显著关联,而在条件关联分析中20和24 个SRAP分子标记分别和籽粒粗蛋白含量、淀粉含量显著关联。同时利用分离群体群分离法,检测到37个SRAP标记与籽粒的饲草性状相关,其中2个影响籽粒粗蛋白含量,4个影响淀粉含量SRAP标记在非条件关联分析中检测到。该研究结果为青稞的籽粒饲草品质性状的改良和高产优质饲草青稞育种提供依据。

关键词: 青稞; 籽粒饲草性状; 非条件与条件关联分析; 分离群体群分离法
中图分类号:S816.15;S540.3 文献标志码:A 文章编号:1001-0629(2018)06-1425-10
Genetic dissection of grain forage traits in hulless barley
Liu Xin-chun1, Lai Yun-ping1,2, Yu Yi3, Yuan Jin-e1,2, Basangyuzhen1,4, Feng Zong-yun1
1.Barley and Hulless Barley Research Centre, College of Agronomy, Sichuan Agricultural University, Chengdu 611130, Sichuan, China
2.Chengdu Agricultural College, Chengdu 611130, Sichuan, China
3.Crop Research Institute, Sichuan Academy of Agricultural Sciences/Chengdu Station for Testing Center of New Varieties of Plants, Ministry of Agriculture, Chengdu 610066, Sichuan, China
4.Hulless Barley Research and Development Centre, Tibet Academy of Agriculture and Animal Husbandry Sciences, Lhasa 850032, Tibet, China
Corresponding author: Feng Zong-yun E-mail:[email protected]
Abstract

To gain genetic information on forage traits for breeding new hulless barley varieties with high yield and good forage quality, we characterized the grain forage traits of 60 accessions of hulless barleys originating from the Qing-Tibet plateau. Unconditional and conditional association mapping analyses, conducted using a mixed line model as implemented in TASSEL software, and bulk sample analysis (BSA) methodology were used to identify markers significantly associated with grain forage traits. High genetic variation in grain forage traits was detected in the association mapping panel. These barley accessions were classified into three subgroups based on population genetic structure analysis, which indicated significant relationships between accession subgroups and their region of origin. Twenty-three sequence-related amplified polymorphism (SRAP) markers for grain forage traits, including crude protein content, starch content, and ash content, were determined using unconditional association analysis, whereas 20 and 24 SRAP markers significantly associated with crude protein content and starch content, respectively, were detected in the conditional association study. Additionally, a total of 37 SRAP markers associated with variation in grain forage traits were detected using the BSA methodology, among which, 2 and 4 SRAP markers associated with crude protein content and starch content, respectively, were also identified in the unconditional association study. These results provide potentially valuable information for improving grain forage traits and the cultivation of good quality, high-yielding accessions of hulless barley.

Key words: hulless barley; grain forage trait; unconditional and conditional association study; bulk sample analysis (BSA)

青稞(Hordeum vulgare var. nudum Hook. f.)是我国青藏高原藏民族对裸大麦的称呼, 是青藏高原海拔4 500 m以上的地区唯一能正常成熟的粮食作物, 是我国藏族人民的口粮和畜牧的饲料的主要来源[1], 随着近年来藏区畜牧业的发展, 传统的主粮型青稞已不适应藏区畜牧业的发展。因此, 培育具有高产特征、优异主粮和饲草品质的青稞品种对于维护藏区社会稳定和促进该地区经济和社会发展有着重要实践意义。目前有关籽粒饲草品质遗传特征的研究主要集中在大麦(H. vulgare)上, 青稞的相关研究较少, 许多研究者对粗蛋白、酸性洗涤纤维、中性洗涤纤维、酸性洗涤木质素等饲草性状进行了QTL定位研究[2, 3, 4, 5], 但这些QTL都是从连锁作图群体中挖掘出来的, 而且这些QTL通常只代表着少数优质饲草品种的遗传特征, 存在检测时间花费多、育种实践效果差等缺点[6]。

为克服连锁作图方法的缺点, 一种基于群体连锁不平衡的关联分析方法已被大量运用于作物的新基因的发掘与相关优质等位基因变异的鉴定[7]。Cai等[8]利用1 319个DArT标记对158份中国西藏大麦资源材料籽粒总蛋白含量性状进行关联作图, 在6条染色体(除了4H染色体)上检测到了与大麦籽粒总蛋白相关的遗传位点, 这些位点所在区域多数是前人报道的籽粒蛋白QTL的热点区域。Shu等[9]利用大麦的9K iselect SNP芯片对欧洲过去100年的大麦商业品种中的抗性淀粉进行了全基因组关联分析, 在10个与籽粒抗性淀粉相关的SNP的LD衰减区域内, 发现了与籽粒淀粉合成相关的结构基因。随后, 运用相同的方法, 在254份欧洲春性大麦中发现了多个与籽粒淀粉相关性状相关联的SNP标记, 并结合大麦基因组的相关信息, 将这些位点定位于淀粉合成与植物细胞壁合成的相关基因的附近[10]。但这些研究只是对影响单个遗传性状变异的基因进行了遗传解析, 而对这些基因在植物发育过程中和受其他因素影响情况的研究较少, 为更好地解析这些基因的影响类型, Zhu[11]提出的净遗传效应的QTL能够较好地解析这些问题, 并同时提高传统QTL的精确性和完整度。该种方法已在小麦(Triticum aestivum)[12]、玉米(Zea mays)[13]、水稻(Oryza sativa)[14]等多种作物上有大量的运用。最近, Jiang等[15]为了能加快水稻的生育期与产量性状进行分子标记聚合育种, 结合条件QTL和关联分析方法, 发现多个同时影响水稻单株有效穗数和生育期的SSR标记, 而这些SSR标记也在前人的文献报道中出现[16]。而运用条件与非条件关联分析的方法发掘调控青稞籽粒饲草相关性状的QTL(基因)未见相关报道, 因此, 本研究运用该方法解析影响青稞品种籽粒饲草相关性状的遗传基础, 以期为优质饲草青稞品种的选育提供一定的理论基础与实践指导。

1 材料与方法
1.1 材料

供试材料来自中国青藏高原地区60份青稞品种。其中, 25份来自四川, 17份来自西藏, 18份来自青海。其材料编号、名称及来源地参见文献[17], 将试验材料于2013年10月种植在四川崇州国家大麦产业体系基地, 每个品种种植 2 行, 行长 2 m, 株距0.4 m, 行距0.3 m。小区按完全随机区组排布, 共2个重复。肥水管理按正常田间水肥管理进行。试验材料于2014年5月按小区混收。

1.2 青稞品种籽粒饲草品质相关性状考察

将每个小区收获的植株混合脱粒, 随机获取200~300粒, 置于博通谷物近红外品质分析仪, 进行青稞籽粒饲草品质相关性状(包括籽粒粗蛋白质含量、总淀粉含量、灰分含量)的非损伤性考察, 每个小区重复3次。

1.3 青稞品种的基因型与群体结构分析

该群体的分子基因型数据为64对显性SRAP标记数据, 具体参见文献[18], 为了降低后续关联分析研究中的假阳性结果出现, 将该群体分子数据中的稀有等位变异频率小于10%的位点剔除。共形成191个SRAP位点。利用Structure 2.3.4软件分析群体结构[19]。软件的参数设定参见文献[20], 并得到相关的Q矩阵。

1.4 青稞籽粒饲草品质相关性状统计分析

利用SPSS21.0软件对群体籽粒饲草相关性状表型数据的描述性统计及相关性分析, 用平均值和标准误表示测定结果。采用基于混合线性模型的数量性状分析软件QGAStation[21]对表型数据进行条件化转化。将某个籽粒饲草品质性状数据表现值转化为基于其他饲草品质性状为同一水平的表现值; 以籽粒粗蛋白含量为例, 将籽粒粗蛋白质含量(Crude protein content)分别转化为基于灰分含量(Ash content)、淀粉含量(Starch content)、纤维含量(Fiber content)为同一水平的条件表型值Crude protein|Ash、Crude protein|Starch、Crude protein|Fiber。籽粒淀粉含量以此处理, 并将各个性状的非转化值命名为该性状的非条件化值。

1.5 青稞籽粒饲草相关性状的条件和非条件值的关联分析

运用TASSEL2.0[22]软件中混合线性模型(MLM)[23]下进行青稞籽粒饲草品质相关性状的条件与非条件关联分析, 以群体分析的Q矩阵和TASSEL软件生成的品种间的亲缘关系矩阵(K)为协变量, 以P=0.05(-lgP=1.3)水平来判定SRAP标记位点与相关性状的条件与非条件值间是否存在关联关系。

1.6 BSA法检测青稞籽粒饲草相关品质性状的关联位点

分别选取青稞籽粒饲草性状具有极端非条件表型值的20个品种, 构建两个品种极端性状表现集团池, 按Takagi等[24]提出的QTL-SEQ理论, 对每个SRAP位点在两个群体池的频率的差值进行统计, 以差值为0.3来判断青稞籽粒饲草品质相关性状与该SRAP位点是否存在着相关性[25]。

2 结果与分析
2.1 青稞品种籽粒饲草品质性状表现分布

利用博通近红外谷物品质分析仪中的大麦分析模块对青稞品种的籽粒饲草品质性状进行相关扫描分析(表1), 结果表明除了粗蛋白在青海地区呈显著差异 (P< 0.05) 外, 所有的饲草品质性状在各个地区内均呈极显著差异(P< 0.01)。青海地区青稞品种的平均粗蛋白含量最高, 为12.78%。西藏地区品种的平均粗蛋白质含量最低, 为12.20%。总群体的粗蛋白平均含量为12.41%, 而最高和最低粗蛋白含量的品种载体分别出现在地区平均值极差的区域, 即最高含量(18.60%)品种出现在青海的门农1号材料中, 而来源于西藏的藏青21号品种拥有最低的粗蛋白含量(8.00%)。淀粉含量以四川地区的青稞品种最大, 为48.55%, 青海地区次之, 为48.24%, 而西藏地区的材料最小, 为47.44%。总群体的平均淀粉含量为48.14%。而西藏地区材料淀粉含量的极差(4.8)远小于其他两个地区, 来源于青海的北青5号表现出最高的淀粉含量(51.70%), 来源于四川的东升黑材料具有最低的淀粉含量(44.50%)。纤维含量总群体的平均值为6.38%, 西藏地区青稞品种的平均纤维含量最高, 四川、青海地区材料的平均纤维含量低于总体平均纤维含量, 表现为西藏> 四川> 青海, 来源于西藏的藏青21号、QB23表现出最高的纤维含量(7.80%), 而来源于青海的东青2号材料具有最低的纤维含量(5.10%)。对于和纤维具有相似性质的灰分, 来源于四川地区的品种平均灰分含量大于总平均值, 西藏和青海地区品种的平均值较为接近, 且小于总体平均值, 3个地区内的灰分极差均为0.5, 60份青稞品种的灰分在1.45%~2.05%, 以四川地区的俄母1号材料最高, 青海地区的北青8号材料最低。总体来看, 总群体粗蛋白含量的变异系数较大(20.31%), 这说明中国青藏高原青稞品种拥有较丰富的控制蛋白质编码位点, 而淀粉和灰分性状的总体变异系数较小, 可能与这两个性状代表的碳水化合物受到一定的人工选择有关。利用SPSS软件中对籽粒的品质性状进行表型相关性分析(表2), 结果发现, 除粗蛋白含量与纤维含量, 其余的品质性状间均呈现显著(P< 0.05)或者极显著(P< 0.01)相关。

表1 青稞品种籽粒饲草品质性状的表型分布 Table 1 Phenotypic distribution of grain forage traits in hulless barley cultivars
表2 青稞品种籽粒饲草品质相关性状的相关性分析 Table 2 Correlation analysis of grain forage traits in hulless barley cultivars
2.2 基于群体结构和条件表型的条件的标记关联作图

2.2.1 群体结构分析 运用Structure软件对60份中国青藏高原青稞品种的群体分化进行分析(图1), 可以看出, 随着K值的不断增加(2~11), 最大似然值的对数值在不断地增加, 在K=3明显出现一个拐点, 初步认为该群体由3个小亚群构成, 结合Evanno等[26]在2005年提出的Δ K的算法对该群体进行了Δ K发掘(图2)。在K=3时, Δ K峰值出现, 认为该群体应由3个亚群构成。结合其来源地信息, 以Q=0.6为阈值判断该品种所属的遗传类群, 除了阿青4号、阿青5号和拉萨勾芒3个品种属于混合类群外, 其余的57份品种按其来源地, 适当地分布于各自的遗传类群里, 遗传类群与来源地的分布相关性高达100%, 这说明青藏高原的青稞品种的选育与种植地适应性有紧密的联系(图3)。

图1 K值与lnP(D)值的关系图Fig. 1 Relationship between K and change in lnP(D) values

图2 Structure软件分析推定的种群分类数K值所对应的Δ K散点图Fig. 2 Δ K for different numbers of populations assumed (K)in the Structure analys

图3 利用Structure软件获得的60份青稞品种的群体遗传结构Fig. 3 Population structure of hulless barley cultivars by Structure software

2.2.2 籽粒饲草品质相关性状的条件与非条件关联分析 运用Tassel软件中MLM模型对青稞品种的饲草品质性状进行非条件关联分析(表3-5)。以P=0.05为阈值, 基于一般混合线性模型共发现了5个控制籽粒粗蛋白含量变异的非条件SRAP标记位点(表3), 5个位点的平均解释表型变异率为9.41%, 其中me11 em14-1位点的R2值最大(13.62%), 其余4个位点的R2相差较小, 在7.34%~8.73%, 以me7 em17-4位点的R2最小。9个影响籽粒总淀粉含量的非条件SRAP标记位点被发掘出来, 9个位点的R2平均值为8.11%(表4), 最高的是me7 em17-4, 为15.13%, 最低的是me10 em15-1, 为5.59%, 除了me7 em17-4位点外, me11 em14-1位点的R2值也大于10%。与淀粉具有同源性的粗灰分含量性状, 也同样发现了9个影响灰分含量的非条件SRAP位点, 9个位点的平均解释表型变异率为7.89%(表5), 范围在5.72%~10.04%, me11 em19-2为最大R2的载体位点, me10 em15-1为最小R2的载体位点。其中me10 em15-1、me6 em13-3和me7 em17-4这3个位点也与籽粒淀粉含量相关联。在P=0.05时, 没有检测到与籽粒纤维含量显著相关的SRAP位点。

表3 青稞品种中与籽粒粗蛋白显著相关的SRAP位点及其对表型变异的解释率 Table 3 Sequence-related amplified polymorphism markers significantly associated with grain crude protein and corresponding explained-phenotypic variation in hulless barley cultivars

运用Zhu[11] 提出的条件QTL的观点, 将相关表型数据进行条件化, 运算得到青稞品种的籽粒饲草品质性状的条件关联标记。对于蛋白含量性状来说(表3), 分别将各种影响籽粒粗蛋白含量饲草品质性状调至同一水平时, 发现了多个与籽粒粗蛋白显著关联的条件SRAP位点, 在将籽粒淀粉含量的影响调整至同一水平时, 分别发现了8个影响籽粒粗蛋白含量的条件性SRAP位点, 这两个性状是所有品质性状中检测到最多的条件性关联位点的性状。纤维则是在粗蛋白条件性关联分析中发掘的最少位点个数(5个)性状。在将淀粉含量调至同一水平时, 发现的8个影响粗蛋白含量的条件性位点, 仅有1个位点在非条件下被检测出来, 而将纤维含量调到同一水平时, 5个条件关联性位点有4个能在非条件下检测出来。对于籽粒淀粉性状(表4), 当消除籽粒灰分含量变异时, 能够检测到最多数目(12个)的影响籽粒淀粉含量的条件关联性位点。而当将籽粒纤维调至同一水平时, 仅有4个条件性位点被检测出来, 2个位点在非条件情况下也能检测, 说明在这2个影响籽粒淀粉的SRAP位点不是通过改变籽粒纤维含量的变异来改变籽粒淀粉含量, 同时也发现在12个去掉籽粒灰分变异而影响籽粒淀粉含量的条件性SRAP位点中, 有8个位点在非条件环境下没有被检测出来, 说明这些调控籽粒淀粉含量是被灰分含量的变异所抑制的, 当灰分变异拉至同一水平, 即灰分因素的影响去掉后, 这些被抑制的位点才能消除抑制影响, 表现出影响籽粒淀粉的表现型出来。

表4 青稞品种中与籽粒淀粉显著相关的SRAP位点及其对表型变异的解释率 Table 4 Sequence-related amplified polymorphism markers significantly associated with grain starch and corresponding explained-phenotypic variation in hulless barley cultivars
表5 青稞品种中与灰分性状显著相关的SRAP位点及其对表型变异的解释率 Table 5 Sequence-related amplified polymorphism markers significantly associated with ash content and corresponding explained-phenotypic variation in hulless barley cultivars
2.3 基于BSA思路发掘籽粒饲草品质相关的SRAP位点

利用QTL-SEQ思路结合群分离基因的观点发掘控制籽粒饲草品质相关的SRAP位点, 以阈值0.3为标准(表6), 分别发现了3个、18个、4个和12个影响籽粒粗蛋白含量、淀粉含量、纤维含量、灰分含量的SRAP位点, me11 em14-1和me11 em17-3两个影响蛋白含量的位点也在非条件环境检测到, me12 em17-1等4个调控籽粒淀粉含量的SRAP位点也在非条件关联分析中检测到, 而影响灰分的SRAP位点均未在非条件关联分析中找到。

表6 青稞品种中基于类Δ SNP index 值发掘与籽粒饲草品质相关性状显著相关的SRAP标记 Table 6 Sequence-related amplified polymorphism markers significantly associated with grain forage traits based on Δ SNP index-like value in hulless barley cultivars
3 讨论
3.1 青稞群体的基于贝叶斯模型的群体结构分析

了解青稞品种的遗传组成与遗传背景是改良现有低产品种和选育新的优质青稞品种的前提。随着现代社会的发展, 作物育种家多使用自己熟悉的原始品种和能具有广泛适应性的品种[27], 这客观上造成了当前作物品种的遗传基础狭窄, 以及改良旧有品种的困难。Yang等[18]利用64对SRAP引物对中国青藏高原青稞品种的遗传多样性进行考察, 得到了青稞品种的遗传多样性低, 遗传基础狭窄的结论。这些不同来源地的品种在群体遗传学上呈现什么样的特征, 对于后期的青稞分子标记育种有重要的影响。本研究针对前期的分子数据, 利用基于贝叶斯模型的Structure 软件[19]对青稞品种的遗传构成进行划分, 发现来自于3个地区的青稞品种由3个亚群构成, 除了3个品种外, 其余57个材料分别分成3个亚群。分群的个体与其来源地存在100%相关性。与原来基于遗传距离的聚类分群结果相比, 按群体遗传学分类更能反映出青稞品的遗传基础, 这与魏世平等[28]利用3种分群法探索中国栽培大豆(Glycine max)群体结构时, 得到的基于贝叶斯模型的Structure软件最适宜反映中国栽培大豆的遗传结构的结论一致, 同时也支持了Kline等[29]认为的群体结构分析更能清楚地反映单个个体趋向群体的趋势和各个亚群体之间的基因交流情况。这些结果也反映出中国当代青稞品种的趋势, 即在当地的适应性良好的品种之间杂交选育品种, 而外来资源利用较少, 不同于中国小麦品种在选育过程大量使用少数几个如南大2419等骨干品种[30], 导致品种与来源地存在没有一定的相关性。但同时从Q值矩阵来看(因数据太多, 结果未列出), 发现在归于某些亚群的个体在归类的概率比较接近, 说明在亚群内部的遗传差异较小, 暗示群体内部的遗传多样性较小, 也验证利用SRAP位点发现该群体遗传狭窄的观点。孟亚雄等[1]、赖勇等[31]、巴桑玉珍等[20]分别利用SSR分子标记解析了中国青稞种质资源材料的遗传多样性和其遗传群体结构, 均发现我国青稞种质资源材料的遗传基础差异小, 群体结构与品种的来源地没有联系(即没有清晰来源地亚群的分类), 但是本研究结果(亚群的分类与地理来源密切相关)却与上述观点相反, 可能与本研究利用的材料为品种, 与前面利用的资源材料有所不同。另外, 本研究利用的SRAP标记是一种基于外显子和内含子边界保守性来设计的目的引物[32], 与SSR主要扩增的基因组中的重复区域有所不同, 不同的基因区域一般反映的基因组的信息不同, 功能区域的变异更能准确地反映植物在进化上的适应情况, 因此本研究的分子数据能更良好地反映中国青藏高原青稞品种的差异。

3.2 青稞籽粒饲草品质性状条件关联分析

在大量作物QTL定位报道研究中, 发现表型上具有紧密相关性的多个性状的QTL位点多以紧密连锁形式或以一因多效在染色体上分布[33]。Zhao等[34]认为这种一因多效的QTL是一种对于两个关系密切、且又相互独立的性状间的遗传关系的条件QTL的表现形式, Jiang等[15]结合条件QTL定位和关联分析的思路, 发现影响水稻生育期的SSR标记中, 有些是通过影响群体的株高或单株有效穗数的表达来调控水稻生育期的表现, 而有些SSR标记则是直接调控水稻生育期的表现。本研究对青稞籽粒饲草相关性状进行遗传解析, 也发现相似的现象; 即影响籽粒饲草品质性状的条件化关联性状的SRAP位点也大量在非条件化关联分析中找到。因此, 条件性关联分析的结果具有进一步更好地对相关性紧密的多性状的遗传特点进行深层次的解析的能力, 随着高通量测序技术和基因组信息的完善, 这种基于条件化的关联分析结果可以与转录组和基因共表达等数据结合起来, 为动植物复杂数量性状的研究提供一种新的思路。

3.3 BSA法检测青稞籽粒饲草性状相关性位点

BSA(bulked sample analysis)是在传统的分离群体群分离法发掘基因的基础上, 结合新一代测序技术发展而来的检测目的性状的功能性分子标记方法, 该方法具有花费少和精确性高等优良特点; 同时运用这种方法的检测标记也不局限于DNA标记, 可以与相关材料的基因表达数据或蛋白表达相结合, 从转录或翻译水平上, 发现与验证与性状相关的标记所联系的功能候选基因[35]。Jin等[36]利用BSA+RNA-SEQ(BSR)方法发掘出了影响茶叶(Camellia sinensis)儿茶素含量变异的重要候选基因F3’ 5’ H基因, 并在儿茶素代谢基因的表达谱数据, 验证了BSR的结果。Zeng等[25]运用BSA的方法, 发现了多个抗柳枝稷(Panicum virgatum)白粉病的SLAF标记, 并在这些标记中发现4个具有NBS结构的R基因。本研究也利用该方法找到69个与青稞籽粒饲草品质相关的SRAP位点, 但与非条件关联分析结果相比, 仅有8个位点是相同的。说明这些位点很可能与青稞品种的饲草品质性状有一定的相关性。但同时也发现了多个未在关联分析筛选出的, 而在BSA方法中获得与籽粒饲草性状相关的SRAP位点, 这可能是由于BSA方法检测相关受到标记种类、混合群体个体数目等多方面影响。针对这些问题, 下一步将对该群体进行多种标记扫描, 并对该群体的饲草品质数据进行多年多点检测, 以期能够找到影响青稞籽粒饲草品质性状稳定表达的遗传位点, 为进一步改良青稞和培育优质饲草青稞品种提供理论依据。

4 结论

本研究利用遗传群体结构分析了中国青藏高原青稞品种的群体结构, 60份品种分成3个亚群, 群体的类群与品种的来源地密切相关, 非条件关联分析一共发现23个影响籽粒饲草品质性状变异的SRAP标记, 影响淀粉或粗灰分含量的SRAP标记最多, 而在条件关联分析研究中发现44个SRAP标记能影响籽粒饲草品质净遗传变异值。发现6个SRAP标记在BSA分析中和非条件关联分析中均能影响籽粒的饲草品质变异, 这些标记有助于加快青稞籽粒饲草性状的分子育种进程。

The authors have declared that no competing interests exist.

参考文献
[1] 孟亚雄, 孟祎林, 汪军成, 司二静, 张海娟, 任盼荣, 马小乐, 李葆春, 杨轲, 王化俊. 青稞遗传多样性及其农艺性状与SSR标记的关联分析. 作物学报, 2016, 42(2): 180-189.
Meng Y X, Meng Y L, Wang J C, Si E J, Zhang H J, Ren P R, Ma X L, Li B C, Yang K, Wang H J. Genetic diversity and association analysis of agronomic characteristics with SSR markers in hulless barley. Acta Agronomica Sinica, 2016, 42(2): 180-189. (in Chinese) [本文引用:2]
[2] Han F, Ullrich S E, Romagosa I, Clancy J A, Froseth J A, Wesenberg D M. Quantitative genetic analysis of acid detergent fiber content in barley grain. Journal of Cereal Science, 2003, 38(2): 167-172. [本文引用:1]
[3] Siahsar B, Peighambari S, Taleii A, Naghavi M R, Nabipour A, Sarrafi A. QTL analysis of forage quality traits in barley ( Hordeum vulgare L. ). Cereal Research Communications, 2009, 37(4): 479-488. [本文引用:1]
[4] Surber L, Abdel-Haleem H, Martin J, Hensleigh P, Cash D, Bowman J, Blake T. Mapping quantitative trait loci controlling variation in forage quality traits in barley. Molecular Breeding, 2011, 28(2): 189-200. [本文引用:1]
[5] Abdel-Haleem H, Bowman J, Surber L, Blake T. Variation in feed quality traits for beef cattle in Steptoe× Morex barley population. Molecular Breeding, 2012, 29(2): 503-514. [本文引用:1]
[6] Gupta P K, Rustgi S, Kulwal P L. Linkage disequilibrium and association studies in higher plants: Present status and future prospects. Plant Molecular Biology, 2005, 57(4): 461-485. [本文引用:1]
[7] Huang X, Han B. Natural variations and genome-wide association studies in crop plants. Annual Review of Plant Biology, 2014, 65(1): 531-551. [本文引用:1]
[8] Cai S, Gang Y, Chen X, Huang Y, Jiang X, Zhang G, Jin X. Grain protein content variation and its association analysis in barley. BMC Plant Biology, 2013, 13(1): 35. [本文引用:1]
[9] Shu X, Backes G, Rasmussen S K. Genome-wide association study of resistant starch (RS) phenotypes in a barley variety collection. Journal of Agricultural and Food Chemistry, 2012, 60(41): 10302-10311. [本文引用:1]
[10] Shu X, Rasmussen S K. Quantification of amylose, amylopectin, and β-glucan in search for genes controlling the three major quality traits in barley by genome-wide association studies. Frontiers in Plant Science, 2014, 5: 197. [本文引用:1]
[11] Zhu J. Analysis of conditional genetic effects and variance components in developmental genetics. Genetics, 1995, 141: 1633-1639. [本文引用:2]
[12] Zhang H, Wang H G. QTL mapping for traits related to P-deficient tolerance using three related RIL populations in wheat. Euphytica, 2015, 203(3): 505-520. [本文引用:1]
[13] Zhang Z, Wu X, Shi C, Wang R, Li S, Wang Z, Liu Z, Xue Y, Tang G, Tang J. Genetic dissection of the maize kernel development process via conditional QTL mapping for three developing kernel-related traits in an immortalized F2 population. Molecular Genetics and Genomics, 2016, 291(1): 437-454. [本文引用:1]
[14] Sun J, Zou D T, Luan F S, Zhao H W, Wang J G, Liu H L, Xie D W, Su D Q, Ma J, Liu Z L. Dynamic QTL analysis of the Na+ content, K+ content, and Na+/K+ ratio in rice roots during the field growth under salt stress. Biologia Plantarum, 2014, 58(4): 689-696. [本文引用:1]
[15] Jiang J H, Chen L, Liu Q, He Y J, Hong D L. Mining elite alleles of growth duration and productive panicle number per plant by association mapping with conditional phenotypic value in japonica rice. Rice Science, 2013, 20(3): 200-206. [本文引用:2]
[16] Jiang J H, Zhao Q B, Liu Q M, Chen L, Chen F L, Qiao B J, Hong D L. Mining applicable elite alleles of growth duration, plant height and panicle number per plant by conditional QTL mapping in japonica rice. Rice Science, 2011, 18(3): 196-203. [本文引用:1]
[17] 刘新春, 苟琳, 杨平, 刘仙俊, 王希文, 何守朴, 李刚, 冯宗云. 青藏高原青稞品种醇溶蛋白的遗传多样性. 植物遗传资源学报, 2008, 9(2): 180-185.
Liu X C, Gou L, Yang P, Liu X J, Wang X W, He S P, Li G, Feng Z Y. Genetic diversity of Hordein on Qingke(hulless barley) varieties from the Qinghai-Tibet plateau of China. Journal of Plant Genetic Resources, 2008, 9(2): 180-185 (in Chinese) [本文引用:1]
[18] Yang P, Liu X J, Liu X C, Yang W Y, Feng Z Y. Diversity analysis of the developed qingke (hulless barley) cultivars representing different growing regions of the Qinghai-Tibet Plateau in China using sequence-related amplified polymorphism (SRAP) markers. African Journal of Biotechnology, 2010, 9(50): 8530-8538. [本文引用:2]
[19] Pritchard J K, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics, 2000, 155(2): 945-959. [本文引用:2]
[20] 巴桑玉珍, 刘新春, 付国勇, 李东梅, 王丹丹, 强小林, 冯宗云. 基于SSR标记揭示青藏高原地区青稞品种种质资源的遗传多样性. 麦类作物学报, 2017, 37(1): 40-47.
Basangyuzhen, Liu X C, Fu G Y, Li D M, Wang D D, Qiang X L, Feng Z Y. Genetic diversity and population structure analysis of hulless barley with cold tolerance from the Qinghai-Tibetan plateau revealed by SSR markers. Journal of Triticeae Crops, 2017, 37(1): 40-47. (in Chinese) [本文引用:2]
[21] Chen G B, Zhu Z X, Zhang F T, Zhu J. Quantitative genetic analysis station for the genetic analysis of complex traits. Chinese Science Bulletin, 2012, 57(21): 2721-2726. [本文引用:1]
[22] Bradbury P J, Zhang Z, Kroon D E, Casstevens T M, Ramdoss Y, Buckler E S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics, 2007, 23(19): 2633-2635. [本文引用:1]
[23] Yu J, Pressoir G, Briggs W H, Vroh B I, Yamasaki M, Doebley J F, cMullen M D, Gaut B S, Nielsen D M, Holland J B, Kresovich S, Buckler E S. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nature Genetics, 2005, 38(2): 203-208. [本文引用:1]
[24] Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, Uemura A, Utsushi H, Tamiru M, Takuno S, Innan H, Cano L M, Kamoun S, Terauchi R. QTL-seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant Journal for Cell and Molecular Biology, 2013, 74(1): 174-183. [本文引用:1]
[25] Zeng B, Yan H D, Liu X C, Zang W, Zhang A, Zhou S, Huang L, Liu J. Genome-wide association study of rust traits in orchardgrass using SLAF-seq technology. Hereditas, 2017, 154(1): 5. [本文引用:2]
[26] Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: A simulation study. Molecular Ecology, 2005, 14(8): 2611-2620. [本文引用:1]
[27] Sharma S, Upadhyaya H D, Varshney R K, Gowda C L. Pre-breeding for diversification of primary gene pool and genetic enhancement of grain legumes. Frontiers in Plant Science, 2013, 4: 309. [本文引用:1]
[28] 魏世平, 刘晓芬, 杨胜先, 吕海燕, 牛远, 章元明. 中国栽培大豆群体结构不同分类方法的比较. 南京农业大学学报, 2011, 34(2): 13-17.
Wei S P, Liu X F, Yang S X, Lyu H Y, Niu Y, Zhang Y M. Comparison of various clustering methods for population structure in Chinese cultivated soybean( Glycine max (L. )Merr. ). Journal of Nanjing Agricultural University, 2011, 34(2): 13-17. (in Chinese) [本文引用:1]
[29] Kline J B, Moore D J, Clevenger C V. Activation and association of the Tec tyrosine kinase with the human prolactin receptor: Mapping of a Tec/Vav1-receptor binding site. Molecular Endocrinology, 2001, 15(5): 832-841. [本文引用:1]
[30] 刘新春, 赖运平, 刘仙俊, 李俊, 袁金娥, 万洪深, 冯宗云, 杨武云. Lr1基因在四个小麦骨干亲本衍生系中的分布及选择效应. 麦类作物学报, 2014, 34(5): 597-602.
Liu X C, Lai Y P, Liu X J, Li J, Yuan J E, Wan H S, Feng Z Y, Yang W Y. Distribution and selective effects of Lr1 gene in the wheat varieties derived from four founder breeding parents in China. Journal of Triticeae Crops, 2014, 34(5): 597-602. (in Chinese) [本文引用:1]
[31] 赖勇, 王晋民, 任龙, 朱惠琴, 马辉, 王金贵, 巨霞, 李宗仁. 大麦SSR标记遗传多样性及连锁不平衡分析. 核农学报, 2016, 30(10): 1889-1897.
Lai Y, Wang J M, Ren L, Zhu H Q, Ma H, Wang J G, Ju X, Li Z R. Genetic diversity and linkage disequilibrium analysis of barley using SSR markers. Journal of Nuclear Agricultural Sciences, 2016, 30(10): 1889-1897. (in Chinese) [本文引用:1]
[32] Li G, Quiros C F. Sequence-related amplified polymorphism (SRAP), a new marker system based on a simple PCR reaction: Its application to mapping and gene tagging in Brassica. Theoretical and Applied Genetics, 2001, 103(2-3): 455-461. [本文引用:1]
[33] 郭蕾蕾. 大麦分子遗传图谱构建及其主要农艺性状和功能成分的QTL定位. 成都: 四川农业大学硕士学位论文, 2012.
Guo L L. Construction of molecular genetic map and QTL mapping for main agronomic characters and functional compounds of barley. Master Thesis. Chengdu: Sichuan Agricultural University, 2012. (in Chinese) [本文引用:1]
[34] Zhao J, Becker H D, Zhang Y, Zhang Y, Ecke W. Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield. Theoretical and Applied Genetics, 2006, 113(1): 33-38. [本文引用:1]
[35] Zou C, Wang P, Xu Y. Bulked sample analysis in genetics, genomics and crop improvement. Plant Biotechnology Journal, 2016, 14(10): 1941-1955. [本文引用:1]
[36] Jin J Q, Ma J Q, Yao M Z, Ma C L, Chen L. Functional natural allelic variants of flavonoid 3', 5'-hydroxylase gene governing catechin traits in tea plant and its relatives. Planta, 2017, 245(3): 523-538. [本文引用:1]