小麦耐低钾性状的全基因组关联分析

摘要

[目的]筛选与小麦耐低钾性状相关的标记,为揭示小麦耐低钾性状的遗传机理奠定基础.[方法]本研究以198份中国黄淮南片麦区的小麦品种(系)作为供试群体.小麦种子发芽后,幼苗在正常钾营养液中生长4天,然后在低钾(0.01 mmol/L)和正常钾(2 mmol/L)营养液中生长25天.调查生物量,分析地上部和根部钾含量,计算小麦7个性状的相对值,利用小麦35K SNP芯片结合Q+K混合线性模型(mixed linear model,MLM)对7个耐低钾性状进行全基因组关联(genome-wide association study,GWAS)分析,筛选位于显著关联位点上的候选基因并进行功能注释,对显著关联位点进行等位变异分析,发掘优异等位变异.[结果]低钾胁迫条件下,地上部、根部、全株钾累积量和钾累积量根冠比显著降低,地上部、根部、全株钾利用指数显著升高.群体结构分析和PCA分析均将供试群体分为2个亚群.通过GWAS分析,共检测到75个显著关联单核苷酸多态性(single nucleotide polymorphism,SNP)标记(P＜0.001),分布在除1A、3B、3D、4D和6A染色体外的16条染色体上.通过候选基因搜寻及注释,共筛选出13个可能与小麦耐低钾胁迫相关的候选基因.等位变异分析共挖掘了14个优异等位变异,其中6个优异等位变异在供试群体中出现的频率较大.[结论]7个耐低钾性状均具有明显的数量性状特征,变异系数范围为20.66％～30.44％.40％的SNP标记分布在2B、5B和6B染色体上,21个SNP位点与多个耐低钾性状显著关联,单个SNP标记解释的表型贡献率的变异范围为5.78％～11.22％.TraesCS4A02G335400、TraesCS2B02G306000、TraesCS5B02G260000可能参与物质转运过程,TraesCS1D02G350600和TraesCSU02G105300可能参与逆境胁迫响应等生理过程,TraesCS2A02G000200可能参与逆境胁迫下的信号转导过程.