删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

春小麦籽粒主要品质性状的全基因组关联分析

本站小编 Free考研考试/2021-12-26

严勇亮,1,2, 时晓磊2, 张金波2, 耿洪伟3, 肖菁2, 路子峰2, 倪中福1, 丛花,21中国农业大学农学院,北京 100193
2新疆农业科学院农作物品种资源研究所,乌鲁木齐 830091
3新疆农业大学农学院/农业生物技术重点实验室,乌鲁木齐 830052

Genome-Wide Association Study of Grain Quality Related Characteristics of Spring Wheat

YAN YongLiang,1,2, SHI XiaoLei2, ZHANG JinBo2, GENG HongWei3, XIAO Jing2, LU ZiFeng2, NI ZhongFu1, CONG Hua,21College of Agronomy and Biotechnology, China Agricultural University, Beijing 100193
2Institute of Crop Germplasm Resources, Xinjiang Academy of Agricultural Sciences, Urumqi 830091
3College of Agriculture, Xinjiang Agricultural University/Key Laboratory of Agricultural Biological Technology, Urumqi 830052

通讯作者: 丛花,E-mail: huacong0924@126.com

责任编辑: 李莉
收稿日期:2021-02-1接受日期:2021-03-17
基金资助:国家自然科学基金(31660389)
国家自然科学基金(U1403185)
乌鲁木齐市科学技术计划(Z161210002)


Received:2021-02-1Accepted:2021-03-17
作者简介 About authors
严勇亮,E-mail: yanliang198279@163.com










摘要
【目的】挖掘小麦籽粒品质性状显著相关的SNP位点及候选基因,并揭示其遗传机理,为相关基因克隆和分子标记辅助选择提供理论依据。【方法】通过检测298份国内外春小麦品种(系)5个环境下蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度、出粉率和容重等7个籽粒品质性状的表型,并结合小麦55K SNP芯片,采用Q+K关联混合模型进行全基因组关联分析(genome-wide association study,GWAS)。【结果】外引品种(系)、地方品种(系)和育成品种(系)的7个品质性状在不同环境下的变异系数分别为1.3%—13.4%、1.1%—18.6%和1.0%—13.9%。其中,外引品种(系)的蛋白质含量、湿面筋含量和沉降值的变异系数均为最高;新疆自育品种的淀粉含量、籽粒硬度和出粉率的变异系数最大,而新疆地方品种的蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度和出粉率6个品质性状的变异系数均介于外引品种(系)和新疆自育品种(系)之间。群体结构分析表明,298份小麦品种(系)可分为3个亚群。其中,亚群1包含128份(43.0%)试验材料,主要是来自新疆的地方品种(系);亚群2包含24份(8.1%)试验材料,主要包括外引品种(系)和新疆地方品种;亚群3包含146份(48.9%)试验材料,主要是外引品种(系)。连锁不平衡分析表明A、B和D基因组及全基因组的LD衰减距离分别为10、10、6和8 Mb,依据全基因组的LD衰减距离,将在物理图谱上前后8 Mb区间内的位点认定为一个候选位点。通过GWAS共检测到85个与7个小麦籽粒品质性状显著关联的稳定位点(P<0.001)贡献率为3.7%—10.9%。在1B、1D、2D、3A、3D、4A、4B、5A、6A、6D、7A和7D染色体上均检测到稳定且同时与多个性状关联的位点。其中,7A染色体上的AX-109452823—AX-110545157同时与蛋白质含量、淀粉含量、湿面筋含量、沉降值、出粉率和籽粒硬度相关,且同时在4个环境中均被检测到。对稳定的位点进行候选基因发掘,筛选到10个可能与小麦籽粒品质相关的候选基因。其中TraesCS4A01G299800(阳离子氨基酸转运蛋白)、TraesCS7A01G059500(色氨酸脱羧酶)、TraesCS7A01G331200TraesCS7D01G418700(木葡聚糖内转葡糖基酶/水解酶)对调控小麦籽粒氨基酸含量有重要作用。【结论】检测到85个稳定的且与小麦籽粒品质性状关联的位点,并筛选出10个与小麦籽粒品质性状相关的候选基因。
关键词: 小麦;籽粒品质性状;全基因组关联分析;SNP;候选基因

Abstract
【Objective】Identify SNPs and candidate genes that are significantly related to wheat grain quality traits, and reveal their genetic mechanism. 【Method】In this study, 298 introduced varieties (lines), Xinjiang landrace (lines) and Xinjiang bred varieties (lines) were used for an association population. Seven grain quality traits, including protein content (PRC), wet gluten content (WGC), sedimentation value (SV), starch content (STC), grain hardness (GH), flour yield (FY) and test weight (TW), were measured under five environments. Based on phenotypes of seven quality traits and genotypes of 55K SNP markers in this population, the Q+K association mixed model was used for genome-wide association analysis to obtain significantly associated SNP loci.【Result】The coefficients of variation of the seven grain quality traits of introduced varieties (lines), landraces and bred varieties (lines) under different environments were 1.3%-13.4%, 1.1%-18.6% and 1.0%-13.9%, respectively. Among them, the protein content, wet gluten content and sedimentation value of introduced varieties (lines) have the highest coefficient of variation (CV); Xinjiang bred varieties (lines) have the largest CV of starch content, grain hardness and flour yield. Whereas, for other six grain quality traits, including protein content, wet gluten content, sedimentation value, starch content, grain hardness, and the coefficients of variation of Xinjiang landraces are all between those of the introduced varieties (lines) and Xinjiang bred varieties (lines). Population structure analysis showed that 298 wheat varieties (lines) can be divided into 3 subgroups. Subgroup 1 contains 128 (43.0%) the materials mainly from landrace (lines); Subgroup 2 has 24 (8.1%) materials, mainly including introduced varieties (lines) and landraces; Subgroup 3 contains 146 (48.9%) materials, mainly introduced varieties (lines). The linkage disequilibrium analysis showed that the LD attenuation distances of the A, B and D genomes and the whole genome respectively were 10, 10, 6 and 8 Mb, according to the LD attenuation distance of the whole genome, the loci in the 8 Mb interval after the physical map were identified as a candidate loci. A total of 85 loci were simultaneously detected in two or more environments, that were significantly associated with 7 wheat grain quality traits detected by GWAS, with a contribution rate of 3.7%-10.9%. Stable SNPs associated with multiple traits were detected on chromosomes 1B, 1D, 2D, 3A, 3D, 4A, 4B, 5A, 6A, 6D, 7A and 7D. Among them, AX-109452823-AX-110545157 on chromosome 7A is related to protein content, starch content, wet gluten content, sedimentation value, flour yield and grain hardness, and was detected across four environments. Candidate genes at stable loci associated with multiple traits were searched, and 10 candidate genes that might be related to wheat grain quality were screened. Among them, TraesCS4A01G299800 (cationic amino acid transporter), TraesCS7A01G059500 (tryptophan decarboxylase), TraesCS7A01G331200, TraesCS7D01G418700 (xyloglucan endoglucosylase/hydrolase) play important roles in regulating the amino acid content in wheat grains.【Conclusion】The 85 loci were simultaneously detected in two or more environments, and 10 candidate genes related to wheat grain quality traits were predicted.
Keywords:wheat;quality traits;genome-wide association analysis;SNP;candidate genes


PDF (702KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
严勇亮, 时晓磊, 张金波, 耿洪伟, 肖菁, 路子峰, 倪中福, 丛花. 春小麦籽粒主要品质性状的全基因组关联分析. 中国农业科学, 2021, 54(19): 4033-4047 doi:10.3864/j.issn.0578-1752.2021.19.001
YAN YongLiang, SHI XiaoLei, ZHANG JinBo, GENG HongWei, XIAO Jing, LU ZiFeng, NI ZhongFu, CONG Hua. Genome-Wide Association Study of Grain Quality Related Characteristics of Spring Wheat. Scientia Acricultura Sinica, 2021, 54(19): 4033-4047 doi:10.3864/j.issn.0578-1752.2021.19.001


开放科学(资源服务)标识码(OSID):

0 引言

【研究意义】小麦是全球种植面积最广、总产量和营养价值最高、加工食品种类最多的粮食作物[1,2,3,4]。随着中国经济的发展和人们生活水平的不断提高,消费者对面制品、保健食品和营养食品的需求不断增加,并且对小麦相关产品的品质要求也越来越高。在供给侧改革需求和产业内循环的新时代要求下,在加工品质研究的基础上,营养和健康品质已成为小麦品质改良的重要方向和重要育种目标[5]。蛋白质含量、湿面筋含量、淀粉含量和沉降值等都是重要的小麦品质性状。其中,蛋白质含量和湿面筋含量是小麦品质分级的重要指标;淀粉含量对面粉品质和加工品质具有重要意义[6,7,8],研究此类性状关键基因的挖掘及遗传机理对于提高小麦品质的遗传改良具有重要的意义。【前人研究进展】小麦籽粒品质性状是由多基因控制的复杂的数量性状,关二旗等[9]通过研究小麦籽粒品质与基因型及环境条件的关系,发现小麦籽粒品质主要受基因型控制,但也受环境条件的影响。在环境因素中,气温和降水量的变化对小麦籽粒品质影响尤为明显。信志红等[10]研究了小麦籽粒品质性状对气象因子的响应,结果表明,小麦籽粒品质主要受基因型控制,其中,淀粉含量和蛋白质含量受环境影响最小。全基因组关联分析(genome-wide association study,GWAS)和数量性状位点(quantitative trait locus,QTL)定位技术的不断发展为小麦籽粒品质性状的分子检测和相关数量性状的研究提供了有力支撑[6, 11]。较QTL定位来说,用于GWAS的群体构建耗时短,群体遗传背景较为丰富,更容易发现重要的数量性状基因位点,挖掘出更多的功能基因,在小麦数量性状研究方面的应用越来越广泛[12]。LI等[13]、ECHEVENY-SOLARTE等[14]、吴云鹏等[15]、郭利建等[6]、黄梦豪等[16]利用RIL群体结合SSR标记、DArT标记、SNP芯片等,在1A、1B、1D、2A、2B、2D、3A、3B、3D、4B、5A、5B、6A、6B、7B和7D染色体上共定位到31个与小麦籽粒蛋白质含量相关的QTL,表型贡献率为0.7%—16.9%。LI等[13]、郭利建等[6]、黄梦豪等[16]利用RIL群体分别结合SSR标记、SNP芯片等,在1A、1B、1D、2A、4A和6B染色体上共定位到20个与湿面筋含量相关的QTL,表型贡献率为0.6%—13.0%。ECHEVENY-SOLARTE等[14]以RIL群体结合包含939个DArT标记的图谱,在1A、1B、2B、3D、4A和6A染色体上共定位到6个与小麦籽粒出粉率相关的QTL,表型贡献率为4.9%—19.0%。吴云鹏等[15]、郭利建等[6]、黄梦豪等[16]利用RIL群体结合SSR标记、SNP芯片等,在1B、1D、2A、3A、3B、5A和7D染色体上共定位到11个与小麦籽粒沉降值相关的QTL,表型贡献率为2.7%—20.2%。郭利建等[6]、黄梦豪等[16]利用RIL群体结合SSR标记、SNP芯片等,在1A、2A、2D、4B、5A和6A染色体上共定位到10个与小麦籽粒淀粉含量相关的QTL,表型贡献率为2.9%—11.2%。【本研究切入点】SNP标记具有分布广、数量多、定位精度高等特点,是目前最具潜力的分子标记[17]。应用GWAS挖掘新疆春小麦品质性状相关基因的研究仍鲜见报道。【拟解决的关键问题】本研究以298份春小麦品种(系)构成的自然群体为材料,利用小麦55K SNP基因芯片分型数据结合蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度,出粉率和容重7个小麦品质相关性状表型鉴定结果,进行全基因组关联分析,发掘小麦品质相关性状的遗传位点,筛选相关候选基因,为相关基因克隆和分子标记辅助选择提供理论依据。

1 材料与方法

1.1 试验材料与表型鉴定

以298份春小麦品种(系)为材料,其中包含国外引进品种(系)110份、新疆地方种141份、新疆育成种47份,该批材料由国家作物种质资源库——新疆分库提供。该群体一定程度上反映了不同时期适合新疆地区种植的小麦种质的特性。2018—2020年连续3个年度将参试材料种植于新疆乌鲁木齐市新疆农业科学院安宁渠基地,分别简称为E1、E2和E4;2019—2020年连续2个年度将参试材料种植于新疆喀什市新疆农业科学院泽普小麦育种家基地,分别简称为E3和E5。采用随机区组设计,每份材料种植3行,行长2 m,行距20 cm,3次重复,人工播种。在试验田四周设置保护行。试验田土壤肥力、灌溉等田间管理基本一致,播种前施基肥25 kg磷酸二铵和5 kg尿素,拔节期和抽穗期均使用滴灌追肥尿素10 kg和磷酸一铵5 kg。由于地方品种白粉病等发病率较高,在抽穗期人工喷施己唑醇进行防治,喷施2次,间隔10 d。抽穗期小麦植株会出现轻微倒伏现象,使用竹竿进行支撑固定以防止出现大面积倒伏。

为保证测量数据的准确性,将小麦材料收获脱粒后,在适宜的环境保存2个月,采用瑞典波通仪器公司生产的DA-7250近红外分析仪对蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度,出粉率和容重7个品质相关性状进行测定,3次重复,取平均值;5个环境下各重复试验每份材料各品质性状结果取平均值。

1.2 表型数据的统计与分析

采用Excel 2016、SPSS 21.0软件对蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度,出粉率和容重7个品质相关性状进行统计分析,并计算各性状的平均值、标准差和变异系数。采用SAS软件进行描述性统计分析和方差分析[18]。广义遗传力用公式H2g2/(σg2ge2Ɛ2/n)计算,其中,σg2为遗传方差,σge2为基因型与环境互作方差,σƐ2为环境方差,n为环境数。

1.3 DNA提取及SNP分型

小麦生长至三叶期时取10株幼苗,将叶片混合后利用改良的CTAB法提取DNA[19]。用1.0%琼脂糖凝胶电泳检验DNA质量和浓度。利用小麦基因组55K SNP芯片对供试材料进行基因分型。小麦55K SNP芯片分型分析及样品质检工作由中玉金标记(北京)生物技术股份有限公司完成,共检测到覆盖小麦全基因组的53 063个SNP标记。采用Tassel 5.0软件对标记信息去冗余,剔除最小等位基因频率(minor allele frequency,MAF)小于5%和缺失率(missing rate)大于25%的SNP,最终得到23 632个多态性高、稳定性好的优质SNP标记用于后续关联分析。

1.4 群体结构分析

应用Power Maker软件计算多态性信息量(polymorphic information content,PIC,PIC=1-ΣP2ij,P2ij表示第i个位点的第j个等位变异出现的频率)[20]。从筛选过的标记中,选取2 000个最小等位基因频率大于10%且在染色体上均匀分布的SNP标记,利用Structure 2.3.4软件进行群体结构分析。参数设置:Length of Burn-Period=10 000,MCMC Reps after Burn- in=100 000,选择Admixture Ancestry模型和Dependent Allele Frequencies模式,令K=2—12,每个K值重复运行5次。将不同K值下重复运行5次的结果上传至Structure Harvester(http://taylor0.biology.ucla.edu/struct_harvest/[21],使用Tassel 5.0软件绘制群体结构分析图。

1.5 连锁不平衡分析

以位点间的相关系数平方(r2)作为衡量多态性位点两两之间的连锁不平衡(linkage disequilibrium,LD)参数。采用TASSEL5.0软件计算r2,以第95百分位的r2值作为阈值估测LD衰减距离。全基因组关联分析中,超过LD衰减距离的2个位点则认为是2个不同的位点,在LD衰减距离内的位点均视为同一位点。

1.6 全基因组关联分析与候选基因筛选

使用TASSEL5.0软件中的Q+K混合线性模型对298份试验材料在5个不同环境下的8个籽粒品质性状的表型值进行性状与标记之间的GWAS,以P=1.0×10-3为阈值,判定SNP标记与目标性状关联的显著性[22],以LD衰减距离为依据将获得的SNP转化为与目标性状显著关联的遗传位点,将在2个及2个以上的环境中发现的位点视为稳定的位点。

将多个环境中稳定出现的SNP标记的延伸序列在小麦TGACv1.0数据库(http://plants.ensembl.org/hmmer/index.html)和NCBI网站(https://www.ncbi.nlm.nih.gov/)进行BLAST比对,对候选位点进行功能注释[23]

2 结果

2.1 小麦籽粒品质性状表型数据的统计分析

通过分析298份供试材料在5个不同环境下的7个小麦籽粒品质性状的表型,发现国外引进品种(系)、新疆地方品种(系)和新疆育成品种(系)的品质性状有一定的差异(表1),3种类型小麦在5个环境下的变异系数分别为1.3%—13.4%、1.1%—18.6%和1.0%—13.9%。其中,容重的变异幅度最小;沉降值的变异幅度最大。不同环境中,国外引进品种(系)的蛋白质含量、湿面筋含量和沉降值平均值最高,分别为17.9%、40.4%和41.0%;新疆育成品种(系)的淀粉含量、籽粒硬度和出粉率最高,分别为67.4%、63.2%和71.3%;新疆地方品种(系)的上述6个品质性状在5个环境下均介于国外引进品种(系)和新疆育成品种(系)之间;容重在各环境中的表现不具有一致性,这个性状可能受环境的影响较大。供试材料7个品质性状的基因型×年份、环境×年份和基因型×环境×年份互作效应间差异均达极显著水平(表2),说明基因型×年份、环境×年份和基因型×环境×年份对于小麦品质性状均有极显著影响。所测品质性状的遗传力为0.61—0.95,其中,容重的遗传力最小;籽粒硬度的遗传力最大,说明籽容重受环境因素的影响较大。供试材料7个品质性状5个环境下的相关性分析显示(电子附表1),容重在5个环境两两之间相关系数平均值较小,分别为0.48;蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度和出粉率5个环境两两之间相关系数平均值为0.64—0.81,且各性状5个环境两两之间相关系数平均值与遗传力呈正相关。

Table 1
表1
表1春小麦品质性状表型变异
Table 1Phenotypic variation of quality traits of wheat in spring wheat
性状
Trait
环境
Environment
最大值 Max最小值 Min平均值 Mean标准差 SD变异系数 CV (%)
外引种
Introduced varieties
地方种
Xinjiang landrace
育成种
Xinjiang bred varieties
外引种
Introduced varieties
地方种
Xinjiang landrace
育成种
Xinjiang bred varieties
外引种
Introduced varieties
地方种
Xinjiang landrace
育成种
Xinjiang bred varieties
外引种
Introduced varieties
地方种
Xinjiang landrace
育成种
Xinjiang bred varieties
外引种
Introduced varieties
地方种
Xinjiang landrace
育成种
Xinjiang bred varieties
蛋白质
含量
PRC
(%)
E122.6022.0319.7913.9613.4014.2218.6716.5316.811.701.711.519.1310.368.98
E221.4023.2419.9714.7213.1812.9317.9516.1615.801.421.871.157.8911.607.27
E322.0222.7918.8413.9213.2413.4817.6116.1715.451.381.481.147.849.127.40
E420.7921.8319.7113.9912.3713.1816.6016.1815.851.421.581.158.549.747.28
E523.1322.7218.5115.0315.1614.5918.5917.6816.451.351.390.917.267.875.54
湿面筋
含量
WGC
(%)
E151.2750.0545.1531.9029.1831.4342.0336.7637.933.884.193.299.2311.408.68
E248.5252.9045.4633.1928.4427.9740.3935.8835.803.334.572.728.2312.737.61
E350.0552.4443.0931.7628.9729.8839.9036.0935.353.213.582.668.059.937.52
E446.7449.7444.8431.8426.9229.1337.6136.1736.143.183.762.678.4610.397.39
E552.2651.9840.8034.3133.3332.4041.8439.5437.673.113.412.077.428.615.51
沉降值
SV
(mL)
E162.2763.1251.6931.0626.6430.2346.6939.7839.176.006.404.8212.8616.0812.31
E258.0360.4745.3830.5723.3025.3842.6135.9334.054.746.683.6311.1318.5810.67
E356.5259.7442.3827.6022.3723.7538.6234.8630.385.195.904.2213.4416.9313.88
E450.5657.7843.5625.2118.6625.0938.5435.7232.194.836.403.9812.5217.9312.37
E556.4864.2147.1226.4229.0828.0638.3239.6635.164.905.363.6312.7813.5110.32
淀粉
含量
STC
(%)
E172.6272.2472.8761.0657.7665.0166.3568.7169.512.462.351.513.713.422.17
E272.3572.9273.9361.9061.8165.9868.2270.0770.781.891.741.322.772.491.86
E371.7371.5072.6660.3362.5466.0567.4568.2770.241.771.631.282.622.391.82
E473.2173.2472.8964.2660.9465.6369.2769.1469.981.841.851.332.652.681.91
E570.0770.2570.5459.8061.6366.9765.9066.7568.691.751.750.932.652.621.36
籽粒
硬度
GH
(%)
E175.0269.3775.1751.9949.5155.7261.6658.7665.974.994.364.098.107.426.20
E276.1672.2973.2053.6051.1353.7566.0259.5167.245.145.364.097.799.006.07
E369.4467.0770.0753.3549.5454.1461.5656.2764.423.923.973.346.377.055.18
E477.1070.2668.8953.0349.7052.7868.6757.6363.445.384.573.577.837.925.62
E570.3073.1270.6545.9848.6453.2957.9755.5563.515.564.283.669.597.705.76
出粉率
FY (%)
E174.3673.4475.2262.3957.4766.4668.4269.5772.062.712.191.823.963.142.52
E280.2477.0676.7368.3768.3771.2873.0972.4674.782.001.551.272.732.141.70
E377.4174.3474.1667.1066.8269.3471.1269.8972.511.591.411.222.242.021.68
E476.5575.4076.9368.1867.3369.9372.4370.8473.181.771.551.342.442.191.82
E575.3874.6573.7267.3864.4567.5371.2068.5570.511.681.561.342.362.281.90
容重
TW (g·L-1)
E1823.25821.55828.75751.62732.95779.65790.40794.85804.7215.0713.2011.081.911.661.38
E2840.12825.69824.17770.63777.00778.85805.11801.52799.7513.329.8110.821.651.221.35
E3818.22817.23825.58756.64771.87788.60796.16795.10804.4010.679.088.391.341.141.04
E4839.55817.75818.78763.24758.88758.43812.36793.52783.9314.499.7911.881.781.231.52
E5824.60821.85822.48774.59764.13771.23796.36796.72797.7711.279.849.691.421.241.21
PRC:蛋白质含量;WGC:湿面筋含量;SV:沉降值;STC:淀粉含量;GH:籽粒硬度;FY:出粉率;TW:容重;SD:标准差;CV:变异系数;E1:2018年泽普环境点;E2:2019年安宁渠环境点;E3:2019年泽普环境点;E4:2020年安宁渠环境点;E5:2020年泽普环境点。下同
PRC: Protein content; WGC: Wet gluten content; SV: Sedimentation value; STC: Starch content; GH: Grain hardness; FY: Flour yield; TW: Test weight; SD: Standard deviation; CV: Variable coefficient; E1: 2018 Zepu environmental point; E2: 2019 Anningqu environmental point; E3: 2019 Zepu environmental point; E4: 2020 Anningqu environmental point; E5: 2020 Zepu environmental point. The same as below

新窗口打开|下载CSV

Table 2
表2
表2春小麦品质性状方差分析
Table 2Analysis of variance of quality traits of spring wheat
性状
Trait
基因型
Genotype
环境
Environment
基因型×环境
G×E
重复
Replicate
误差
Error
遗传力
H2
蛋白质含量PRC30.50***323.25***2.71***28.92***0.94***0.92
湿面筋含量WGC173.71***1575.75***14.49***138.56***5.07***0.92
沉降值SV396.29***5503.76***38.53***482.08***12.76***0.91
淀粉含量STC42.33***1067.05***4.58***50.85***1.38***0.90
籽粒硬度GH406.19***4375.56***20.43***20.52***3.83***0.95
出粉率FY41.49***1793.73***4.76***5.17***0.91***0.89
容重TW1308.20***7599.86***267.46***247.03***54.04***0.82
***:在P<0.001水平差异显著
***: Significant at P<0.001

新窗口打开|下载CSV

2.2 SNP多态性及分布

从小麦55K SNP芯片中筛选出具有多态性的SNP标记23 632个,分布在A、B和D基因组染色体上的SNP位点数目分别为11 559(48.9%)、10 366(43.9%)和1 707(7.2%)个,可以看出,B基因组的多态性最高;D基因组的多态性明显低于A和B 2个基因组的(电子附表2)。在小麦21条染色体中,2A染色体上的标记数目最多(2 033个),4D染色体上的标记数目最少(117个)。3个染色体组的PIC表现为B(0.3119)>A(0.3036)>D(0.2319)。全基因组的PIC变异范围在0.0068—0.3750,平均值为0.3108。其中2A染色体的PIC最大(0.3423),4D染色体的最小(0.2204)。说明选用的标记在各染色体组上分布较均匀,且标记的等位变异间具有高度的多态性。

2.3 群体结构分析和连锁不平衡分析

利用Structure2.3.4软件对群体遗传结构进行分析,对每个可能的K值模拟运算,用K值与ΔK值做图(图1-A),在K=3处ΔK值最大,曲线变化程度最大。由图1-B可以看出,供试材料被分为3个亚群,其中亚群1包含128(43.0%)份,主要来自新疆的地方品种(系);亚群2包含24(8.1%)份,主要包括外引品种(系)和新疆地方种(系);亚群3包含146(48.9%)份,主要是外引品种(系),还包括新疆育成品种(系)。因此将K=3时生成的Q矩阵,用于性状与标记的关联分析。

图1

新窗口打开|下载原图ZIP|生成PPT
图1298份小麦品种(系)群体结构分析

A:群体的∆K值;B:群体结构示意图
Fig. 1Population structure analysis of 298 wheat varieties (lines)

A: Estimation of ∆K value in population; B: Group structure diagram


考虑连锁群(linkage group)间连锁不平衡背景的影响,对连锁群间的r2值进行平方根转换,以大于此分布95%的参数值为阈值,用来截取同一连锁群内LD的衰减距离[24],经计算得到A、B和D基因组及全基因组的LD衰减距离分别为10、10、6和8 Mb,依据全基因组的LD衰减距离,将在物理图谱上前后8 Mb区间内的位点认定为一个候选位点(图2)。

图2

新窗口打开|下载原图ZIP|生成PPT
图2连锁不平衡衰减图

Fig. 2Linkage disequilibrium plot



2.4 全基因组关联分析

对298份供试材料在5个环境下的籽粒品质性状表型结合55K基因芯片分型的23 632个SNP标记进行GWAS,并采用MLM(mixed linear model,MLM)Q+K关联混合模型,共有85个位点同时在2个及2个以上环境中检测到,分布在除4D和6B染色体外的19条染色体上,是比较稳定的位点(表3)。在1B、1D、2D、3A、3D、4A、4B、5A、6A、6D、7A和7D染色体上检测到稳定的、且同时与多个性状关联的位点(表4)。其中,7A染色体上的8.91—16.64 Mb和26.47—32.70 Mb区段、7D染色体上的534.44—540.87 Mb和566.20—566.20 Mb区段检测到的稳定且同时与多个性状关联的位点可作为后续研究的重点。

Table 3
表3
表3春小麦品质性状显著关联的稳定的位点信息
Table 3Stable locus information with significant correlation of spring wheat quality traits
性状
Trait
位点
Marker
染色体
Chr.
物理位置
Position (Mb)
P
P value
表型贡献率
R2 (%)
环境
Environment
沉降值
SV (mL)
AX-1099237771B335.511.25E-05—4.48E-056.96—7.83E4, E5
AX-1093893041B445.676.44E-05—4.57E-045.49—6.79E4, E5
AX-949351571D470.73—478.862.57E-05—3.85E-045.42—7.54E2, E5
AX-1094866953A725.151.15E-04—5.59E-045.11—6.23E2, E5
AX-1089861463B547.592.89E-05—3.18E-045.56—7.22E4, E5
AX-1115835894A121.515.34E-05—7.95E-044.95—6.94E4, E5
AX-1104173164A130.421.05E-05—6.22E-056.78—7.98E4, E5
AX-1112648924A241.731.40E-04—2.24E-045.79—6.17E4, E5
AX-1096087304A413.762.96E-05—3.05E-045.67—7.25E4, E5
AX-1092851364A594.90—598.799.75E-05—8.96E-043.82—6.39E2, E3, E4
AX-1098340404A609.822.7E-04—6.83E-044.98—5.65E4, E5
AX-1100027194B43.562.69E-04—9.30E-043.79—4.69E2, E4
AX-945135184B556.062.13E-07—5.41E-045.24—10.89E2, E3
AX-1100232365A30.667.47E-04—9.28E-044.80—4.95E2, E4
AX-1093741285A600.177.45E-05—4.00E-045.43—6.61E4, E5
AX-1101498095B554.322.12E-04—3.56E-045.45—5.92E3, E5
AX-1099839325B604.615.60E-04—5.87E-045.15—5.18E2, E3
AX-1094025635D28.883.63E-05—1.86E-046.05—7.11E4, E5
AX-1100829185D529.225.79E-05—3.09E-045.57—6.72E4, E5
AX-1095123426A25.80—28.333.91E-05—4.89E-045.25—7.11E2, E3, E4, E5
AX-945643076D21.984.64E-05—7.11E-044.98—6.93E2, E4
AX-1116148756D70.623.43E-05—1.81E-045.96—7.10E4, E5
AX-1094528237A26.47—32.705.93E-05—9.46E-043.72—6.34E1, E2, E4, E5
AX-1094775467B155.48—159.711.99E-05—7.71E-045.03—10.36E2, E3, E4, E5
AX-944207397D540.879.04E-04—9.43E-044.83—5.08E2, E4
出粉率
FY (%)
AX-1094280411A22.484.14E-04—7.46E-045.02—5.24E3, E5
AX-948022451D73.73—77.361.14E-04—1.00E-034.88—6.34E3, E4
AX-1110995552B18.94—24.135.37E-05—3.64E-044.42—5.72E3, E4
AX-1110400454A632.76-659.661.51E-04—9.74E-043.83—5.96E3, E5
AX-1106714786A3.80—6.573.32E-06—9.73E-043.78—7.77E3, E4
AX-1115348756A5.51—12.281.44E-05—3.98E-045.61—7.91E1, E3, E4
AX-1117455186A17.54—24.411.51E-04—4.85E-046.14—6.15E1, E4
AX-1116118856D2.49—5.821.57E-04—7.67E-045.00—6.14E3, E4
AX-1110764187A8.91—12.441.70E-04—8.12E-043.75—5.14E4, E5
蛋白质含量
PRC (%)
AX-949351571D470.73—478.188.52E-04—9.52E-043.85—4.77E2, E4
AX-946948982B620.642.57E-04—8.16E-044.87—5.71E2, E4
AX-1112182822D523.363.71E-04—8.94E-044.82—5.47E2, E4
AX-945893282D624.106.02E-04—8.85E-045.10—5.12E2, E4
AX-1104423133A510.31—510.602.91E-04—8.05E-044.88—5.62E2, E4
AX-1100027194B43.561.20E-04—2.83E-044.59—5.13E2, E4, E5
AX-1115608745B10.501.38E-04—5.85E-045.11—6.15E2, E4
AX-1093007336A12.28—17.545.03E-05—8.33E-044.87—6.98E2, E4
AX-951945866D77.523.12E-04—9.12E-044.97—5.68E2, E4
AX-1116547437A32.00—32.708.79E-05—8.48E-043.83—5.30E1, E2, E3, E5
AX-1087978507B453.772.79E-04—7.67E-044.93—5.66E2, E4
AX-1087531317D534.44—540.873.72E-04—6.37E-045.17—5.70E2, E4
AX-1089743577D566.20—566.234.25E-04—8.84E-044.99—5.36E2, E4
AX-1109459867D609.24—610.232.75E-04—8.28E-044.97—5.68E2, E3
淀粉含量
STC (%)
AX-951801001D359.672.12E-05—2.57E-045.75—7.55E4, E5
AX-945893282D624.10—629.218.13E-05—8.67E-044.86—6.63E1,E2
AX-1099347633D559.16—561.013.24E-05—9.76E-044.78—7.56E1, E2
AX-1104173164A130.42—136.974.48E-04—8.12E-044.93—5.59E1, E2
AX-1093329134A701.532.32E-04—7.15E-045.26—5.79E1, E2
AX-1101694145A575.262.31E-04—5.72E-044.07—4.90E1, E2
AX-1100465176A8.09—12.282.58E-05—8.16E-045.14—7.43E1, E3
AX-1117455186A17.545.73E-04—8.12E-044.90—5.42E1, E2
AX-1088884436D16.23—21.986.92E-04—9.04E-044.82—5.26E1, E2
AX-1105451577A32.702.86E-04—7.03E-044.14—4.54E1, E2
AX-1089743577D566.20—566.232.60E-04—7.07E-045.00—5.80E1, E2
容重
TW (g·L-1)
AX-1093016332A84.154.89E-04—9.53E-044.50—5.44E1, E3
AX-1101735802A697.04—698.012.69E-04—4.71E-044.23—6.47E1, E2
AX-861845112B158.15—164.432.03E-04—9.14E-044.53—4.79E1, E4
AX-1109466104A608.03—616.897.28E-04—7.84E-043.93—4.71E1, E2
AX-1099474315B672.48—678.512.29E-04—6.26E-045.26—5.84E3, E5
AX-1104028767D626.923.38E-04—4.10E-045.43—5.56E2, E5
湿面筋含量
WCG (%)
AX-1100027194B43.561.27E-04—3.97E-044.32—5.06E2, E4, E5
AX-1115608745B10.502.83E-04—9.94E-044.71—5.61E2, E4
AX-1117455186A17.54—25.803.01E-04—6.74E-044.97—5.75E2, E4
AX-1116547437A32.00—32.707.64E-05—7.13E-044.72—5.35E1, E2
AX-944207397D540.876.47E-04—6.52E-045.04—-5.22E2, E4
籽粒硬度
GH (%)
AX-1111076791A476.08—483.941.32E-04—6.54E-044.92—5.01E4, E5
AX-1109182371A487.48—489.185.71E-05—8.41E-044.76—6.59E4, E5
AX-1094994081A496.31—500.661.77E-04—6.78E-043.94—5.91E3, E5
AX-1103690381B627.15—628.513.72E-04—9.45E-043.81—4.74E1, E2, E3
AX-945346661D457.15—457.953.16E-06—5.85E-044.01—7.42E2, E3, E5
AX-1116364402A505.33—507.361.45E-04—9.12E-043.92—5.34E1, E3
AX-944474032D439.445.38E-04—7.59E-044.84—5.39E1, E3
AX-1116404323A728.42—731.081.70E-04—5.00E-044.12—4.76E2, E3, E5
AX-1109517563B7.087.16E-06—8.58E-044.73—8.11E2, E3
AX-1110182093D553.59—560.753.97E-04—8.39E-044.98—5.28E1, E5
AX-1110314575A571.71—571.949.82E-05—1.82E-044.72—5.12E3, E5
AX-1089459486A9.14—13.061.00E-04—9.09E-045.00—6.02E1, E5
AX-1088028717A10.09—16.642.96E-04—8.09E-044.32—5.85E1, E5
AX-1116644637A123.20—127.558.90E-05—6.93E-043.86—6.29E4, E5
AX-1104954967A484.39—489.952.30E-04—9.39E-043.72—5.65E3, E4, E5

新窗口打开|下载CSV

Table 4
表4
表4同时与多个品质性状显著关联的稳定位点信息
Table 4Stable loci significantly associated with multiple quality traits
性状
Trait
位点
Marker
染色体
Chr.
物理位置
Position (Mb)
P
P value
表型贡献率
R2 (%)
环境
Environment
沉降值,蛋白质含量
SV, PRC
AX-949351571D470.73—478.862.57E-05—9.52E-043.85—7.54E2, E4, E5
淀粉含量,蛋白质含量
STC, PRC
AX-945893282D624.10—629.218.13E-05—8.85E-044.86—6.63E1, E2, E4
籽粒硬度,沉降值
GH, PRC
AX-1116404323A725.15—731.081.15E-04—5.59E-044.12—6.23E2, E3, E5
籽粒硬度,淀粉含量
GH, STC
AX-1110182093D553.59—561.013.24E-05—9.76E-044.78—7.56E1, E2, E5
淀粉含量,沉降值
STC, SV
AX-1104173164A130.42—136.971.05E-05—8.12E-044.93—7.98E1, E2, E4, E5
容重,沉降值
TW, SV
AX-1109466104A608.03—616.892.70E-04—7.84E-043.93—5.65E1, E2, E4, E5
沉降值,蛋白质含量,湿面筋含量
SV, PRC, WGC
AX-1100027194B43.561.20E-04—9.30E-043.79—5.13E2, E4, E5
籽粒硬度,淀粉含量
GH, STC
AX-1110314575A571.71—575.269.82E-05—5.72E-044.07—5.12E1, E2, E3, E5
蛋白质含量,湿面筋含量
PRC, WGC
AX-1115608745B10.501.38E-04—9.94E-044.71—6.15E2, E4
出粉率,淀粉含量,籽粒硬度
FY, STC, GH
AX-1115348756A5.51—12.281.44E-05—9.09E-045.00—7.91E1, E3, E4, E5
淀粉含量,蛋白质含量
STC, PRC
AX-1117455186A12.28—17.545.03E-05—8.33E-044.87—6.98E1, E2, E4
出粉率,湿面筋含量
FY, WGC
AX-1117455186A17.54—25.801.51E-04—6.74E-044.97—6.15E1, E2, E4
淀粉含量,沉降值
STC, SV
AX-1088884436D16.23—21.984.64E-05—9.04E-044.82—5.69E1, E2, E4
沉降值,蛋白质含量
SV, PRC
AX-1116148756D70.62—77.523.43E-05—9.12E-044.97—7.10E2, E4, E5
籽粒硬度,出粉率
GH, FY
AX-1088028717A8.91—16.641.70E-04—8.12E-043.75—5.85E1, E4, E5
沉降值,淀粉含量,蛋白质含量,湿面筋含量 SV, STC, PRC, WGCAX-1094528237A26.47—32.705.93E-05—9.46E-043.72—6.34E1, E2, E3, E4, E5
沉降值,湿面筋含量,蛋白质含量
SV, WGC, PRC
AX-944207397D534.44—540.873.72E-04—9.43E-044.83—5.70E2, E4
蛋白质含量,淀粉含量
PRC, STC
AX-1089743577D566.20—566.232.60E-04—8.84E-044.99—5.80E1, E2, E4
加粗的位点为发现的疑似新的位点
The thickened loci is a suspected new loci

新窗口打开|下载CSV

位于1B染色体445.67 Mb检测到与沉降值相关的位点;3A染色体510.31—510.60 Mb与蛋白质含量相关的位点;6A染色体5.51—12.28 Mb同时与出粉率、淀粉含量和籽粒硬度相关的位点;6A染色体12.28—17.54 Mb同时与淀粉含量和蛋白质含量相关的位点;6A染色体17.54—25.80 Mb同时与出粉率和湿面筋含量相关的位点;7A染色体26.47—32.70 Mb同时与沉降值、淀粉含量、蛋白质含量和湿面筋含量相关的位点;以上6个位点可能为新的位点(表4)。

2.5 候选基因筛选

将与小麦籽粒品质性状显著关联的稳定位点SNP标记序列在普通小麦中国春基因组数据库搜索,获取基因序列,在NCBI数据库中进行BLASTx[25],共得到10个与小麦籽粒品质性状相关的基因(表5)。其中,TraesCS4A01G299800编码阳离子氨基酸转运蛋白(cationic amino acid transporter,putative);TraesCS6A01G054400编码F-box家族蛋白(F-box family protein);TraesCS7A01G059500编码色氨酸脱羧酶(tryptophan decarboxylase);TraesCS6A01G012800编码富含亮氨酸的重复受体样蛋白激酶家族蛋白(leucine-rich repeat receptor-like protein kinase family protein);TraesCS4B01G054200编码组蛋白H2B(Histone H2B);TraesCS1D01G382000编码锌指状蛋白(zinc finger protein-like);TraesCS3A01G507800编码含螺旋结构域的蛋白(coiled- coil domain-containing protein);TraesCS7A01G331200TraesCS7D01G418700编码木葡聚糖内转葡糖基酶/水解酶(xyloglucan endotransglucosylase/hydrolase);TraesCS7D01G445400对应结构为ATP依赖性Clp蛋白酶ATP结合亚基(ATP-dependent Clp protease ATP-binding subunit)。

Table 5
表5
表5稳定位点筛选获得候选基因信息
Table 5Selection of stable sites and information of candidate genes
性状
Trait
位点
Marker
染色体
Chr.
物理位置
Position (Mb)
基因
Gene
基因注释或编码蛋白
Gene annotation or coding protein
沉降值
SV
AX-1092851364A594.90—598.79TraesCS4A01G299800阳离子氨基酸转运蛋白,推定
Cationic amino acid transporter, putative
AX-1095123426A25.80—28.33TraesCS6A01G054400box家族蛋白 F-box family protein
蛋白质含量,沉降值
PRC, SV
AX-1094528237A26.47—32.70TraesCS7A01G059500色氨酸脱羧酶
Tryptophan decarboxylase
出粉率
FY
AX-1115348756A5.51—12.28TraesCS6A01G012800富含亮氨酸的重复受体样蛋白激酶家族蛋白
Leucine-rich repeat receptor-like protein kinase family protein
湿面筋含量,
蛋白质含量 WGC, PRC
AX-1100027194B43.56TraesCS4B01G054200组蛋白H2B
Histone H2B
籽粒硬度
GH
AX-945346661D457.15—457.95TraesCS1D01G382000锌指状蛋白 Zinc finger protein-like
AX-1116404323A728.42—731.08TraesCS3A01G507800含螺旋结构域的蛋白Coiled-coil domain-containing protein
AX-1104954967A484.39—489.95TraesCS7A01G331200木葡聚糖内转葡糖基酶/水解酶
Xyloglucan endotransglucosylase/hydrolase
蛋白质含量
PRC
AX-1089743577D566.20—566.23TraesCS7D01G445400ATP依赖性Clp蛋白酶ATP结合亚基
ATP-dependent Clp protease ATP-binding subunit
AX-1087531317D534.44—540.87TraesCS7D01G418700木葡聚糖内转葡糖基酶/水解酶
Xyloglucan endotransglucosylase/hydrolase

新窗口打开|下载CSV

3 讨论

3.1 春小麦籽粒品质性状表型变异

小麦是全球主要粮食作物之一。自20世纪80年代以来,小麦品质遗传改良一直是育种家比较关注的领域[6]。本研究以298份春小麦为材料,通过近红外谷物分析仪测定5个环境下小麦籽粒品质性状。结果显示,不同环境下国外引进品种(系)的蛋白质含量、湿面筋含量和沉降值均最高;新疆育成品种(系)的淀粉含量、籽粒硬度和出粉率最高,这可能是由于新疆地区对于小麦品质的关注较晚一些,之前新疆小麦育种家在小麦品种选育时关注小麦的产量、出粉率等比较直观的性状,而忽略了肉眼不可见的品质性状。近些年,小麦品质逐渐成为新疆小麦育种的主要目标之一。本研究和张金波等[26]研究均表明,近十几年来,育种家以新疆地方种和国外引进品种(系)为亲本选育了一批适合在新疆种植的品质较优的春小麦品种,这些品种的蛋白质含量和湿面筋含量较之前主要种植的地方品种都有显著提高,这也说明了小麦育种家不再单一注重小麦的产量和出粉率,也逐渐关注小麦品质遗传改良的结果。本研究与金欣欣等[27]研究均表明,小麦籽粒品质性状主要由遗传因素控制,因此,可以通过挖掘与小麦品质性状相关的遗传位点和候选基因对春小麦进行遗传改良,为优质小麦新品种选育提供参考。

3.2 小麦品质的生态适应性

本研究表明,7个小麦籽粒品质性状的基因型×年份、环境×年份和基因型×环境×年份互作效应间差异均达极显著水平(表2),说明基因型×年份、环境×年份和基因型×环境×年份对于小麦品质性状均有极显著影响,基因型和环境是决定小麦品质性状的重要因素;但是,7个品质性状在5个环境的相关性分析显示,容重在5个环境两两之间相关系数平均值较小(0.48),同时其遗传力也相对较小(0.82);蛋白质含量、湿面筋含量、沉降值、淀粉含量、籽粒硬度和出粉率在5个环境两两之间相关系数平均值为0.64—0.81,遗传力为0.89—0.95;且遗传力大的性状其5个环境两两之间相关系数平均值也较大,说明基因型是决定小麦品质性状的主要因素,这与燕丽等[3]的研究结果一致。与北疆地区种植的(E1、E2和E4)小麦品质性状相比,南疆地区种植的小麦籽粒蛋白质含量和湿面筋含量较高;而沉降值、淀粉含量、籽粒硬度、出粉率和容重相对较低,这是由于不同的品质性状对环境变化的响应有所差异所导致的。光照和温湿度对面筋含量和蛋白质含量影响显著,这与张影全等[28]的研究结果一致。信志红等[10]和吴新元等[29]的研究也表明,光照和温度与小麦籽粒蛋白质含量呈正相关。同一生育期南疆地区的日照和温度均高于北疆地区,这可能是导致南疆地区小麦籽粒蛋白质含量和湿面筋含量较高的主要原因。可以通过挖掘调控小麦品质的遗传位点和功能基因,选育出适宜南、北疆地区种植的春小麦品种,扩大优质小麦的种植面积,以期解决现阶段人们高质量的饮食需求。

3.3 春小麦品质性状全基因组关联分析

全基因组关联分析在小麦遗传改良方面的应用越来越广泛[8]。LI等[13]在1B、1D、6B和7D染色体上检测到与湿面筋含量相关的位点,与本研究检测到的位点一致。本研究在7D染色体上的540.87 Mb处还检测到与湿面筋含量显著关联的环境稳定位点。黄梦豪等[16]、杨林等[30]也在7D染色体上检测到与小麦品质性状相关的位点,进一步验证了7D染色体对于调控小麦品质性状的重要作用。本研究与吴云鹏等[15]均在3A和3B染色体上检测到与蛋白质含量相关的位点,在1B、1D和3B染色体上检测到与沉降值相关的位点,表明这些染色体上可能存在调控小麦籽粒蛋白质含量和沉降值的基因。本研究和LOU等[31]在2B染色体上18.94—24.13 Mb这一区段内同时检测到与小麦籽粒品质性状相关的位点,说明2B染色体的这一区段在调控小麦籽粒品质性状方面起重要作用,可以作为下一步研究的重点。本研究与杨林等[30]、沈玮囡等[32]均在1B和3B染色体上检测到与小麦籽粒品质性状相关的位点。本研究结合前人的研究结果筛选出多个与小麦籽粒品质性状相关的相同位点和重要染色体区段,这些位点具有深入研究的价值。

LOU等[31]在3A染色体上的484.64 Mb处检测到与蛋白质含量相关的稳定位点。本研究在3A染色体的510.31—510.60 Mb区段检测到与蛋白质含量相关的稳定位点。这两个位点在染色体上的距离为25.67—25.96 Mb,可以确定不是同一位点,本研究发现的这一位点可能是新的控制小麦籽粒品质性状的遗传位点。本研究在3B染色体上的547.59 Mb处检测到的位点与LOU等[31]在3B染色体554.03 Mb处检测到的位点在染色体上的距离为6.44 Mb,根据全基因组LD衰减距离(8 Mb)可以确定这两个位点为同一位点;但位于1B染色体上的445.67 Mb区段与沉降值相关的位点与LOU等[31]在1B染色体上的471.52 Mb处检测到的位点在染色体上的距离为25.85 Mb,也可以确定不是同一位点,同样也未见相关报道,可能是新的位点。本研究与郭利建等[6]、杨林等[30]在5A和6A染色体上检测到的与蛋白质含量、湿面筋含量等相关的位点。其中本研究在6A染色体上的3.80—25.80 Mb区间内发现多个稳定的位点,且未见有关于这一区段内控制品质性状位点报道。LOU等[31]在7A染色体的6.85 Mb处检测到与淀粉含量相关的位点与本研究在7A染色体上的8.91—12.44 Mb处检测到的与出粉率和籽粒硬度相关的位点在染色体上的距离为2.06—5.59 Mb,根据全基因组LD衰减距离(8 Mb)可以确定这两个位点为同一位点。本研究在7A染色体上的26.47—32.70 Mb区段检测到的与沉降值、淀粉含量、蛋白质含量和湿面筋含量相关的位点未见报道,这一位点可能是新的控制小麦品质性状的遗传位点。

3.4 候选基因预测

将全基因组关联分析检测到的与小麦籽粒品质性状显著关联、且稳定的位点及重要区段进行候选基因筛选,在7A染色体上的484.39—489.95 Mb区段和7D染色体上的534.44—540.87 Mb区段处筛选出与籽粒硬度相关的候选基因TraesCS7A01G331200和与蛋白质含量相关的候选基因TraesCS7D01G418700编码木葡聚糖内转葡糖基酶/水解酶(xyloglucan endotransglucosylase/hydrolase,XTH)。前人研究表明,XTH基因在植株体内表现出非常复杂的功能[33]。它参与了种子萌发[34]、根系建成[35]、花形成开放以及果实发育成熟[36]等重要生物学过程。关于XTH基因在小麦籽粒品质方面的作用还需进一步研究。本研究还筛选出编码F-box家族蛋白(TraesCS6A01G054400)、富含亮氨酸的重复受体样蛋白激酶家族蛋白(TraesCS6A01G012800)、锌指状蛋白(TraesCS1D01G382000)等结构的候选基因,这些基因在调控小麦抗逆性方面有重要作用[37,38,39]。关于这些基因在品质方面的作用,还需要进一步研究。小麦中含有多种人体必需氨基酸,同时,小麦籽粒中赖氨酸等氨基酸的含量对小麦的营养品质有重要作用[40],本研究在4A染色体上的594.90—598.79 Mb区段筛选到编码阳离子氨基酸转运蛋白(TraesCS4A01G299800)的候选基因;在7A染色体上的26.47—32.70 Mb区段内筛选到编码色氨酸脱羧酶(TraesCS7A01G059500)的候选基因;在6A染色体上的5.51—12.28 Mb区段内筛选到编码富含亮氨酸的重复受体样蛋白激酶家族蛋白(TraesCS6A01G012800)的候选基因;这些候选基因可能与小麦籽粒中氨基酸含量有关,需要进一步对其进行深入研究,为小麦品质改良和新品种选育奠定基础。

4 结论

通过对298份春小麦品种(系)3年5点的7个籽粒品质性状结合55K SNP基因芯片进行全基因组关联分析,发现85个稳定的显著关联位点,包括18个与多个性状关联的多效性位点;在1B(445.67 Mb)、3A(510.31—510.60 Mb)、6A(5.51—12.28 Mb、12.28—17.54 Mb和17.54—25.80 Mb)和7A(26.47—32.70 Mb)染色体上的位点可能是新的品质性状位点。

参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

王立祥, 廖允成. 中国粮食问题. 北京: 阳光出版社, 2012.
[本文引用: 1]

WANG L X, LIAO Y C. Food Issues in China. Beijing: Sunshine Press, 2012. (in Chinese)
[本文引用: 1]

赵广才, 常旭虹, 王德梅, 陶志强, 王艳杰, 杨玉双, 朱英杰. 小麦生产概况及其发展
作物杂志, 2018(4):1-7.

[本文引用: 1]

ZHAO G C, CHANG X H, WANG D M, TAO Z Q, WANG Y J, YANG Y S, ZHU Y J. General situation and development of wheat production
Crops, 2018(4):1-7. (in Chinese)

[本文引用: 1]

燕丽, 王志忠, 郑文寅, 张文明, 郭文善, 姚大年. 基因型和环境对安徽小麦品质性状的影响
麦类作物学报, 2016, 36(11):1497-1501.

[本文引用: 2]

YAN L, WANG Z Z, ZHENG W Y, ZHANG W M, GUO W S, YAO D N. Effects of genotype and environment on wheat quality in Anhui province
Journal of Triticeae Crops, 2016, 36(11):1497-1501. (in Chinese)

[本文引用: 2]

关二旗, 魏益民, 张波, 郭进考, 张国权, 刘彦军, 罗勤贵, 班进福. 黄淮冬麦区部分区域小麦品种构成及品质性状分析
中国农业科学, 2012, 45(6):1159-1168.

[本文引用: 1]

GUAN E Q, WEI Y M, ZHANG B, GUO J K, ZHANG G Q, LIU Y J, LUO Q G, BAN J F. Analysis of the variety composition and quality properties of wheat in a part of the Yellow-Huai river zone
Scientia Agricultura Sinica, 2012, 45(6):1159-1168. (in Chinese)

[本文引用: 1]

张爱民, 李欣, 刘冬成, 孙家柱, 阳文龙. 品质支撑农作物产业与未来发展
中国农业科学, 2016, 49(22):4265-4266.

[本文引用: 1]

ZHANG A M, LI X, LIU D C, SUN J Z, YANG W L. Quality-the future of crop production
Scientia Agricultura Sinica, 2016, 49(22):4265-4266. (in Chinese)

[本文引用: 1]

郭利建, 王竹林, 汪世娟, 刘振华, 刘香利, 胡胜武, 赵惠贤. 基于SRAP和SSR标记的小麦品质相关性状的QTL定位
麦类作物学报, 2016, 36(10):1275-1282.

[本文引用: 8]

GUO L J, WANG Z L, WANG S J, LIU Z H, LIU X L, HU S W, ZHAO H X. QTL mapping of wheat grain quality traits based on SRAP and SSR marker
Journal of Triticeae Crops, 2016, 36(10):1275-1282. (in Chinese)

[本文引用: 8]

刘建军, 何中虎, 赵振东, 宋建民, 刘爱峰. 小麦面条加工品质研究进展
麦类作物学报, 2001(2):81-84.

[本文引用: 1]

LIU J J, HE Z H, ZHAO Z D, SONG J M, LIU A F. Review of noodle industrial quality of wheat
Journal of Triticeae Crops, 2001(2):81-84. (in Chinese)

[本文引用: 1]

陈锋, 何中虎, 崔党群, 赵武善, 张艳, 王德森. 利用近红外透射光谱技术测定小麦品质性状的研究
麦类作物学报, 2003(3):1-4.

[本文引用: 2]

CHEN F, HE Z H, CUI D Q, ZHAO W S, ZHANG Y, WANG D S. Measurement of wheat quality traits by near infrared transmittance spectroscopy
Journal of Triticeae Crops, 2003(3):1-4. (in Chinese)

[本文引用: 2]

关二旗, 魏益民, 张波. 小麦籽粒品质与基因型及环境条件的关系
麦类作物学报, 2010, 30(5):963-969.

[本文引用: 1]

GUAN E Q, WEI Y M, ZHANG B. Relationships between wheat kernel quality and genotype as well as environmental conditions
Journal of Triticeae Crops, 2010, 30(5):963-969. (in Chinese)

[本文引用: 1]

信志红, 郭建平, 谭凯炎, 刘凯文, 杨荣光, 张利华, 孙义. 冬小麦籽粒品质评价及其对气象因子的响应研究
中国生态农业学报, 2019, 27(8):1205-1217.

[本文引用: 2]

XIN Z H, GUO J P, TAN K Y, LIU K W, YANG R G, ZHANG L H, SUN Y. Evaluation of grain quality of winter wheat and its response to meteorological factors
Chinese Journal of Eco-Agriculture, 2019, 27(8):1205-1217. (in Chinese)

[本文引用: 2]

NELSON J C, ANDREESCU C, BRESEGHELLO F, FINNEY P L, GUALBERTO D G, BERGMAN C, PENA R J, PERRETANT M R, LEROY P, QUALSET C O, SORRELLS M E. Quantitative trait locus analysis of wheat quality traits
Euphytica, 2006, 149:145-159.

DOI:10.1007/s10681-005-9062-7URL [本文引用: 1]

TURUSPEKOV Y, BAIBULATOVA A, YENMEKBAYEY K, TOKHETOVA L, CHUDINOV V, SEREDA G, GANAL M, GRIFFITHS S, ABUGALIEVA S. GWAS for plant growth stages and yield components in spring wheat (Triticum aestivum L.) harvested in three regions of Kazakhstan
BMC Plant Biology, 2017, 17(S1):51-61.

DOI:10.1186/s12870-017-0999-1URL [本文引用: 1]

LI H M, TANG Z X, ZHANG H Q, YAN B J, REN Z L. Major quality trait analysis and QTL detection in hexaploid wheat in humid rain-fed agriculture
Genetics and Molecular Research, 2013, 12(2):1740-1751.

DOI:10.4238/2013.May.21.5URL [本文引用: 3]

ECHEVENY-SOLARTE M, KUMAR A, KIANIAN S, SIMSEK S, ALAMRI M S, MANTOVANI E E, MCCLEAN P E, DECKARD E L, ELIAS E, SCHATZ B, XU S S, MERGOUM M. New QTL alleles for quality-related traits in spring wheat revealed by RIL population derived from supernumerary×non-supernumerary spikelet genotypes
Theoretical and Applied Genetics, 2015, 128(5):893-912.

DOI:10.1007/s00122-015-2478-0URL [本文引用: 2]

吴云鹏, 张业伦, 肖永贵, 阎俊, 张勇, 张晓科, 张利民, 夏先春, 何中虎. 小麦重要品质性状的QTL定位
中国农业科学, 2008, 41(2):331-339.

[本文引用: 3]

WU Y P, ZHANG Y L, XIAO Y G, YAN J, ZHANG Y, ZHANG X K, ZHANG L M, XIA X C, HE Z H. QTL mapping for important quality traits in common wheat
Scientia Agricultura Sinica, 2008, 41(2):331-339. (in Chinese)

[本文引用: 3]

黄梦豪, 刘天相, 强琴琴, 李春莲, 王中华. 基于SNP和SSR标记的小麦品质性状的QTL定位
分子植物育种, 2019, 17(12):3966-3973.

[本文引用: 5]

HUANG M H, LIU T X, QIANG Q Q, LI C L, WANG Z H. QTL mapping of wheat grain quality traits based on SNP and SSR marker
Molecular Plant Breeding, 2019, 17(12):3966-3973. (in Chinese)

[本文引用: 5]

WURSCHUM T, LANGER S M, LONGIN C F, KORZUN V, AKHUNOV E, EBMEYER E, SCHACHSCHNEIDER R, SCHACHT J, KAZMAN E, REIF J C. Population structure, genetic diversity and linkage disequilibrium in elite winter wheat assessed with SNP and SSR markers
Theoretical and Applied Genetics, 2013, 126(6):1477-1486.

DOI:10.1007/s00122-013-2065-1URL [本文引用: 1]

MENG L, LI H, ZHANG L, WANG J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations
The Crop Journal, 2015, 3(3):269-283.

DOI:10.1016/j.cj.2015.01.001URL [本文引用: 1]

陈昆松, 李方, 徐昌杰, 张上隆, 傅承新. 改良CTAB法用于多年生植物组织基因组DNA的大量提取
遗传, 2004, 26(4):529-531.

[本文引用: 1]

CHEN K S, LI F, XU C J, ZHANG S L, FU C X. Modified CTAB method for extraction of genome DNA from perennial plant tissues
Hereditas, 2004, 26(4):529-531. (in Chinese)

[本文引用: 1]

WANG S X, ZHU Y L, ZHANG D X, SHAO H, LIU P, HU J B, ZHANG H, ZHANG H P, CHANG C, LU J, XIA X C, SUN Z L, MA C X. Genome-wide association study for grain yield and related traits in elite wheat varieties and advanced lines using SNP markers
PLoS ONE, 2017, 12(11):1-14.

[本文引用: 1]

ZHU C S, GORE M, BUCKLER E S, BUCKLER, YU J M. Status and prospects of association mapping in plants
The Plant Genome, 2008, 1(1):5-20.

[本文引用: 1]

YU J, PRESSOIR G, BRIGGS W H, BI I V, YAMASAKI M, DOEBIEY J F, MCMULLEN 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, 2006, 38(2):203-208.

DOI:10.1038/ng1702URL [本文引用: 1]

翟俊鹏, 李海霞, 毕惠惠, 周思远, 罗肖艳, 陈树林, 程西永, 许海霞. 普通小麦主要农艺性状的全基因组关联分析
作物学报, 2019, 45(10):1488-1502.

[本文引用: 1]

ZHAI J P, LI H X, BI H H, ZHOU S Y, LUO X Y, CHEN S L, CHENG X Y, XU H X. Genome-wide association study for main agronomic traits in common wheat
Acta Agronomica Sinica, 2019, 45(10):1488-1502. (in Chinese)

[本文引用: 1]

BRESEGHELLO F, SORRELLS M E. Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars
Genetics, 2006, 172(2):1165-1177.

DOI:10.1534/genetics.105.044586URL [本文引用: 1]

周思远, 毕惠惠, 程西永, 张旭睿, 闰永行, 王航辉, 毛培钧, 李海霞, 许海霞. 小麦耐低磷相关性状的全基因组关联分析
植物遗传资源学报, 2020, 21(2):431-445.

[本文引用: 1]

ZHOU S Y, BI H H, CHENG X Y, ZHANG X R, RUN Y H, WANG H H, MAO P J, LI H X, XU H X. Genome-wide association study of low-phosphorus tolerance related traits in wheat
Journal of Plant Genetic Resources, 2020, 21(2):431-445. (in Chinese)

[本文引用: 1]

张金波, 严勇亮, 王小波, 路子峰, 肖菁, 彭惠茹, 丛花. 新疆春小麦育成品种遗传演变分析
新疆农业科学, 2020, 57(3):418-426.

[本文引用: 1]

ZHANG J B, YAN Y L, WANG X B, LU Z F, XIAO J, PENG H R, CONG H. Analysis of the genetic evolution of cultivated spring wheat varieties in Xinjiang
Xinjiang Agricultural Sciences, 2020, 57(3):418-426. (in Chinese)

[本文引用: 1]

金欣欣, 姚艳荣, 贾秀领, 姚海坡, 申海平, 崔永增, 李谦. 基因型和环境对小麦产量、品质和氮素效率的影响
作物学报, 2019, 45(4):635-644.

[本文引用: 1]

JIN X X, YAO Y R, JIA X R, YAO H B, SHEN H P, CUI Y Z, LI Q. Effects of genotype and environment on wheat yield, quality, and nitrogen use efficiency
Acta Agronomica Sinica, 2019, 45(4):635-644. (in Chinese)

[本文引用: 1]

张影全, 唐娜, 张波, 景东林, 马永安, 魏益民. 冀南地区小麦籽粒品质现状及利用潜力分析
麦类作物学报, 2018, 38(2):157-163.

[本文引用: 1]

ZHANG Y Q, TANG N, ZHANG B, JING D L, MA Y A, WEI Y M. Study on kernel quality property and processing potency of wheat in the south of Hebei province
Journal of Triticeae Crops, 2018, 38(2):157-163. (in Chinese)

[本文引用: 1]

吴新元, 芦静, 张新忠, 黄天荣, 李建疆, 周安定, 梁晓东, 曹俊梅, 高永红, 曾潮武. 新疆小麦品质生态区划研究
新疆农业科学, 2017, 54(8):1373-1383.

[本文引用: 1]

WU X Y, LU J, ZHANG X Z, HUANG T R, LI J J, ZHOU A D, LIANG X D, CAO J M, GAO Y H, ZENG C W. Study of ecological division for wheat quality in Xinjiang
Xinjiang Agricultural Sciences, 2017, 54(8):1373-1383. (in Chinese)

[本文引用: 1]

杨林, 吴青霞, 邵慧, 冉从福, 余静, 李立群, 李学军. 小麦籽粒品质性状的QTL分析
西北植物学报, 2013, 33(8):1574-1583.

[本文引用: 3]

YANG L, WU Q X, SHAO H, RAN C F, YU J, LI L Q, LI X J. QTL mapping for grain quality traits in wheat
Acta Botanica Boreali-Occidentalia Sinica, 2013, 33(8):1574-1583. (in Chinese)

[本文引用: 3]

LOU H Y, ZHANG R Q, LIU Y T, GUO D D, ZHAI S S, CHEN A Y, ZHANG Y F, XIE C J, YOU M S, PENG H R, LIANG R Q, NI Z F, SUN Q X, LI B Y. Genome-wide association study of six quality-related traits in common wheat (Triticum aestivum L.) under two sowing condition
Theoretical and Applied Genetics, 2021, 134:399-418

DOI:10.1007/s00122-020-03704-yURL [本文引用: 5]

沈玮囡, 王竹林, 杨睿, 李美霞, 梁子英, 奚亚军, 孙风丽, 刘曙东. 波兰小麦品系XN555×普通小麦品系中13衍生重组自交系(RILs)群体中籽粒品质相关性状QTL定位
农业生物技术学报, 2014, 22(5):561-571.

[本文引用: 1]

SHEN W N, WANG Z L, YANG R, LI M X, LIANG Z Y, XI Y J, SUN F L, LIU S D. QTL analysis of grain quality related traits using recombinant inbred lines (RILs) derived from the cross of Triticum polonicum L. line XN555×T. aestivum L. line Zhong 13
Journal of Agricultural Biotechnology, 2014, 22(5):561-571. (in Chinese)

[本文引用: 1]

贾鑫磊, 何贝轩, 郭丹丹, 郭美丽. 膨胀素和木葡聚糖内转葡糖基酶/水解酶基因的功能研究进展
植物生理学报, 2018, 54(11):1659-1668.

[本文引用: 1]

JIA X L, HE B X, GUO D D, GUO M L. Research progress in the function of expansins and xyloglucan endotransglucosylase/hydrolase
Plant Physiology Journal, 2018, 54(11):1659-1668. (in Chinese)

[本文引用: 1]

MIEDES E, ZARRA I, HOSON T, HERBERS K, SONNEWALD U, LORENCES E P. Xyloglucan endotransglucosylase and cell wall extensibility
Journal of Plant Physiology, 2011, 168(3):196-203.

DOI:10.1016/j.jplph.2010.06.029URL [本文引用: 1]

XU P, CAI X T, WANG Y, XING L, CHEN Q, XIANG C B. HDG11 upregulates cell-wall-loosening protein genes to promote root elongation in Arabidopsis
Journal of Experimental Botany, 2014, 65(15):4285-4295.

DOI:10.1093/jxb/eru202URL [本文引用: 1]

HAN Y, HAN S, BAN Q, HE Y H, JIN M J, RAO J P. Overexpression of persimmon DkXTH1 enhanced tolerance to abiotic stress and delayed fruit softening in transgenic plants
Plant Cell Reports, 2017, 36:583-596.

DOI:10.1007/s00299-017-2105-4URL [本文引用: 1]

茹京娜, 于太飞, 陈隽, 陈明, 周永斌, 马有志, 徐兆师, 闵东红. 小麦锌指转录因子TaDi19A对低温的响应及其互作蛋白的筛选
中国农业科学, 2017, 50(13):2411-2422.

[本文引用: 1]

RU J N, YU T F, CHEN Y, CHEN M, ZHOU Y B, MA Y Z, XU Z S, MIN D H. TaDi19A of zinc finger transcription factors in wheat in response to low temperature
Scientia Agricultura Sinica, 2017, 50(13):2411-2422. (in Chinese)

[本文引用: 1]

魏春茹, 孟钰玉, 范润侨, 赵梦伊, 于秀梅, 赵伟全, 康振生, 刘大群. 小麦F-box/Kelch类基因TaFKOR23的抗逆相关表达模式及分子互作蛋白鉴定
植物遗传资源学报, 2020, 21(3):695-705.

[本文引用: 1]

WEI C R, MENG Y Y, FAN R Q, ZHAO M Y, YU X M, ZHAO W Q, KANG Z S, LIU D Q. Stress-related expression profile of F-box/Kelch gene TaFKOR23 in wheat and molecular characterization of the interacting target protein
Journal of Plant Genetic Resources, 2020, 21(3):695-705. (in Chinese)

[本文引用: 1]

邱丽丽, 赵琪, 张玉红, 戴绍军. 植物质膜蛋白质组的逆境应答研究进展
植物学报, 2017, 52(2):128-147.

DOI:10.11983/CBB16001 [本文引用: 1]
质膜作为原生质体与外界环境的屏障, 除了维持正常的细胞内稳态和营养状况, 还参与感知和应答各种环境胁迫。近年来, 植物质膜蛋白质组学研究为深入分析植物应答不同生物和非生物胁迫的分子机制提供了重要信息, 已经报道了模式植物拟南芥(Arabidopsis thaliana)和水稻(Oryza sativa)等10种植物质膜应对生物胁迫(白叶枯病菌(Xanthomonas oryzae pv. oryzae)感染)与非生物胁迫(冷、盐、水淹、渗透、高pH值、Fe缺乏及过量、氮素、脱落酸、壳聚糖和壳寡糖)过程的蛋白质丰度模式变化。通过整合分析植物质膜响应逆境的蛋白质组学研究结果, 揭示了质膜在植物应答逆境胁迫过程中的重要作用。植物通过调节转运蛋白、通道蛋白及膜泡运输相关蛋白的丰度变化促进细胞内外的信号传递、物质交换与运输; 同时利用膜相关的G蛋白、Ca<sup>2+</sup>信号、磷酸肌醇信号途径及BR信号途径等多种信号通路, 通过蛋白质可逆磷酸化作用感知和传递胁迫信号, 调节植物抵御胁迫。研究结果为从蛋白质水平认识质膜逆境应答分子调控机制提供了新线索。
QIU L L, ZHAO Q, ZHANG Y H, DAI S J. Advances in adversity response in proteome of plant plasma membrane
Chinese Bulletin of Botany, 2017, 52(2):128-147. (in Chinese)

DOI:10.11983/CBB16001 [本文引用: 1]
质膜作为原生质体与外界环境的屏障, 除了维持正常的细胞内稳态和营养状况, 还参与感知和应答各种环境胁迫。近年来, 植物质膜蛋白质组学研究为深入分析植物应答不同生物和非生物胁迫的分子机制提供了重要信息, 已经报道了模式植物拟南芥(Arabidopsis thaliana)和水稻(Oryza sativa)等10种植物质膜应对生物胁迫(白叶枯病菌(Xanthomonas oryzae pv. oryzae)感染)与非生物胁迫(冷、盐、水淹、渗透、高pH值、Fe缺乏及过量、氮素、脱落酸、壳聚糖和壳寡糖)过程的蛋白质丰度模式变化。通过整合分析植物质膜响应逆境的蛋白质组学研究结果, 揭示了质膜在植物应答逆境胁迫过程中的重要作用。植物通过调节转运蛋白、通道蛋白及膜泡运输相关蛋白的丰度变化促进细胞内外的信号传递、物质交换与运输; 同时利用膜相关的G蛋白、Ca<sup>2+</sup>信号、磷酸肌醇信号途径及BR信号途径等多种信号通路, 通过蛋白质可逆磷酸化作用感知和传递胁迫信号, 调节植物抵御胁迫。研究结果为从蛋白质水平认识质膜逆境应答分子调控机制提供了新线索。

陈雪燕, 王灿国, 程敦公, 李豪圣, 宋健民, 刘爱峰, 王利彬, 董爽爽, 赵振东, 刘建军, 曹新有. 小麦加工品质相关贮藏蛋白、基因及其遗传改良研究进展
植物遗传资源学报, 2018, 19(1):1-9.

[本文引用: 1]

CHEN X Y, WANG C G, CHENG D G, LI H S, SONG J M, LIU A F, WANG L B, DONG S S, ZHAO Z D, LIU J J, CAO X Y. Research progress on wheat processing quality related storage proteins, genes and genetic improvement
Journal of Plant Genetic Resources, 2018, 19(1):1-9. (in Chinese)

[本文引用: 1]

相关话题/环境 基因 新疆 遗传 植物