小麦籽粒砷含量全基因组关联分析及候选基因预测

doi:10.15933/j.cnki.1004-3268.2025.10.003

摘要/Abstract

摘要： 以包含207个小麦品种（系）的自然群体为材料，利用小麦660K SNP芯片进行基因分型，并结合不同环境下表型数据和最佳线性无偏预测（BLUP）值对小麦籽粒砷（As）含量进行全基因组关联分析（GWAS），鉴定其遗传位点，预测候选基因，并利用单倍型分析、线性回归分析方法对小麦籽粒As含量稳定遗传SNP位点的优势、劣势单倍型及其聚合效应进行分析，为小麦籽粒As生物强化育种提供理论依据。结果表明，小麦自然群体籽粒As含量在多个环境下均存在明显变异，各环境下籽粒As含量均呈正态分布。采用混合线性模型（MLM）在全基因组水平上共鉴定到67个SNP与小麦籽粒As含量显著关联，其中，位于1A、2A、2D、3B、6B染色体上的6个SNP（AX-110393422、AX-110370420、AX-109339465、AX-109455432、AX-109359598、AX-111655266）在多个环境下表现稳定，为稳定遗传数量性状核苷酸位点（QTN）；同时发现在6B染色体末端存在1个SNP热点区域，峰值SNP AX-111655266在所有环境下均被检测到，且与小麦籽粒As含量关联程度最高，可解释表型变异的15.14%~21.53%，为主效QTN。基于小麦表达谱公共数据和自然群体籽粒灌浆20 d的候选基因差异表达分析结果，结合物理位置、基因功能注释等信息进行候选基因预测，发现与SNP AX-111655266 紧密连锁的候选基因TraesCS6B02G418600编码类磷酸转运蛋白PHO1，可能参与As的跨膜运输和信号转导，为控制小麦籽粒As含量的主效基因；其余5个SNP位点的候选基因分别为TraesCS1A02G346700、TraesCS2A02G106800、TraesCS2D02G171500、TraesCS3B02G131500、TraesCS3B02G252500，分别编码F-box蛋白、硫胺素焦磷酸激酶1、类GEM蛋白1、类泛素结合酶E2蛋白、蛋白磷酸酶2C家族蛋白，为控制小麦籽粒As含量的候选基因。6个稳定遗传SNP位点在所有环境下，优势单倍型群体籽粒As含量均极显著低于劣势单倍型群体。此外，随着优势单倍型数量的增加，小麦籽粒As含量逐渐降低，具有明显的加性效应。

关键词: 小麦, 籽粒As含量, 全基因组关联分析, 优势和劣势单倍型, 聚合育种, 候选基因预测

Abstract: In order to provide a theoretical basis for the arsenic（As） biofortification breeding in wheat grains，an association population comprising 207 wheat varieties（lines）was used as material，genotyping was done by wheat 660K SNP array，combining with phenotypic data under different environments and best linear unbiased prediction（BLUP）value，a genome‐wide association study（GWAS）was conducted on As content in wheat grains，genetic loci were identified，and candidate genes associated with As content in wheat grains were predicted. The haplotype analysis and linear regression analysis methods were used to analyze the superior and inferior haplotypes of stable SNP loci associated with As content in wheat grains and their aggregated effect. The results showed that there were obvious variation in grain As content under multiple environments，and the grain As content exhibited a normal distribution.By using a mixed linear model（MLM），a total of 67 SNPs were identified to be significantly associated with grain As content at the genome‐wide level. Among these，six SNPs located on chromosomes 1A，2A，2D，3B，and 6B（AX‐110393422，AX‐110370420，AX‐109339465，AX‐109455432，AX‐109359598，AX‐111655266）were consistently detected under multiple environments，and were considered as stable quantitative trait nucleotides（QTNs）.Additionally，a SNP hotspot was identified at the end of chromosome 6B.The peak SNP AX‐111655266 was detected under all environments，which exhibited the strongest association with grain As content，explained the phenotypic variation of 15.14%—21.53% ，and was a major‐effect QTN.Based on public wheat expression profile data and differential expression analysis of candidate genes at 20 days post‐anthesis in the association population，combined with physical location and functional annotation information，candidate genes were predicted.The candidate gene TraesCS6B02G418600，tightly linked to SNP AX‐111655266，encoding a phosphate transporter PHO1‐like protein，which may be involved in transmembrane transport and signal transduction of As，was a major‐effect gene controlling grain As content in wheat. The candidate genes for the other five SNP loci were TraesCS1A02G346700，TraesCS2A02G106800，TraesCS2D02G171500，TraesCS3B02G131500，and TraesCS3B02G252500，encoding F‐box protein，thiamine pyrophosphokinase 1，GEM‐like protein 1，ubiquitin‐onjugating enzyme E2‐like protein，and protein phosphatase 2C family protein，respectively. They were candidate genes involved in regulating As accumulation in wheat grains.The grain As content in the superior haplotype group was significantly lower than that in the inferior haplotype group under all environments.Furthermore，as the number of superior haplotypes increased，the grain As content gradually decreased，demonstrating a clear additive effect.

Key words: Wheat, As content in grain, Genome‐wide association study, Superior and inferior haplotypes, Pyramid breeding, Candidate gene prediction

中图分类号:

S512

时夏, 孙东玲, 常阳, 李文旭, 代资举, 周正富, 吴政卿, 雷振生. 小麦籽粒砷含量全基因组关联分析及候选基因预测[J]. 河南农业科学, 2025, 54(10): 21-32.

SHI Xia, SUN Dongling, CHANG Yang, LI Wenxu, DAI Ziju, ZHOU Zhengfu, WU Zhengqing, LEI Zhensheng. Genome‐Wide Association Study and Candidate Gene Prediction of Arsenic Content in Wheat Grains[J]. Journal of Henan Agricultural Sciences, 2025, 54(10): 21-32.

图/表 10

参考文献

［1］MONDAL P， MAJUMDER C B， MOHANTY B.Laboratory based approaches for arsenic remediation from contaminated water：Recent developments［J］.Journal of Hazardous Materials，2006，137（1）：464‐479.

［2］SMITH A H，LOPIPERO P A，BATES M N，et al. Public health. Arsenic epidemiology and drinking water standards［J］. Science，2002，296（5576）：2145‐2146.

［3］WU G R，HONG H L，YAN C L. Arsenic accumulation and translocation in mangrove（Aegiceras corniculatum L.）grown in arsenic contaminated soils［J］.International Journal of Environmental Research and Public Health，2015，12（7）：7244‐7253.
［4］SU Q T，DU Z X，HUANG X Y，et al. Managing arsenic pollution from soil‐plant systems：Insights into the role of biocha［r J］. Plants，2025，14（10）：1553.

［5］FENDORF S，MICHAEL H A，VAN GEEN A. Spatial and temporal variations of groundwater arsenic in South and Southeast Asia［J］. Science，2010，328（5982）：1123‐1127.

［6］HUG S J，WINKEL L H E，VOEGELIN A A，et al.Arsenic and other geogenic contaminants in groundwater：A global challenge［J］. Chimia，2020，74（7/8）：524‐537.

［7］YU L X，LIU X C，BOGE W，et al. Genome‐wide association study identifies loci for salt tolerance during germination in autotetraploid alfalfa（Medicago sativa L.）using genotyping‐by‐sequencing［J］. Frontiers in Plant Science，2016，7：956.

［8］WAN X M，LEI M，CHEN T B. Review on remediation technologies for arsenic‐contaminated soil［J］. Frontiers of Environmental Science & Engineering，2019，14（2）：24.

［9］DAI Y C，LÜ J L，LIU K，et al.Major controlling factors and prediction models for arsenic uptake from soil to wheat plants［J］.Ecotoxicology and Environmental Safety，2016，130：256‐262.

［10］MUNIR M，AHMAD K，KHAN Z I，et al. Trophic transfer of arsenic in the food chain from selected crops grown in long‐term wastewater contaminated soil［J］. Agricultural Water Management，2025，309：109345.

［11］BASU A，SAHA D，SAHA R，et al. A review on sources，toxicity and remediation technologies for removing arsenic from drinking water［J］.Research on Chemical Intermediates，2014，40（2）：447‐485.

［12］YAMAJI N，SAKURAI G，MITANI‐UENO N，et al.Orchestration of three transporters and distinct vascular structures in node for intervascular transfer of silicon in rice［J］. Proceedings of the National Academy of Sciences of the United States of America，2015，112（36）：11401‐11406.

［13］CHEN Y，SUN S K，TANG Z，et al. The Nodulin 26‐like intrinsic membrane protein OsNIP3；2 is involved in arsenite uptake by lateral roots in rice［J］.Journal of Experimental Botany，2017，68（11）：3007‐3016.

［14］ROSAS‐CASTOR J M，GUZMÁN‐MAR J L，HERNÁNDEZ‐RAMÍREZ A，et al. Arsenic accumulation in maize crop（Zea mays）：A review［J］.Science of the Total Environment，2014，488：176‐187.

［15］陈国梁，冯涛，陈章，等. 砷在农作物中的累积及其耐受机制研究综述［J］.生态环境学报，2017，26（11）：1997‐2002.

CHEN G L，FENG T，CHEN Z，et al. Research review on the mechanisms of arsenic uptake and its resistance in crops［J］. Ecology and Environmental Sciences，2017，26（11）：1997‐2002.

［16］ZAHEER M S，AIJAZ N，HAMEED A，et al. Cultivating resilience in wheat：Mitigating arsenic toxicity with seaweed extract and Azospirillum brasilense［J］.Frontiers in Microbiology，2024，15：1441719.

［17］ANGON P B，ISLAM M S，SHREEJANA K，et al.Sources，effects and present perspectives of heavy metals contamination：Soil，plants and human food chain［J］. Heliyon，2024，10（7）：e28357.

［18］SHI X，GENG S H，HOU J N，et al. Unveiling novel genetic loci and superior alleles for nickel accumulation in wheat via genome‐wide association study［J］. Plants，2025，14（8）：1262.

［19］BHATTA M，STEPHEN BAENZIGER P，WATERS B M，et al. Genome‐wide association study reveals novel genomic regions associated with 10 grain minerals in synthetic hexaploid wheat［J］. International Journal of Molecular Sciences，2018，19（10）：3237.

［20］ZENG Z Y，OU D D，YIN L，et al. Genetic identification and characterization of two novel loci for grain arsenic concentration in tetraploid wheat at various field environments［J/OL］. Plant and Soil，2025［2025‐05‐20］.https：//doi.org/10.1007/s11104‐025‐07303‐y.

［21］WANG W，GUO H，WU C N，et al. Identification of novel genomic regions associated with nine mineral elements in Chinese winter wheat grain［J］. BMC Plant Biology，2021，21（1）：311.

［22］CUI F，ZHANG N，FAN X L，et al. Utilization of a wheat 660K SNP array‐derived high‐density genetic map for high‐resolution mapping of a major QTL for kernel number［J］. Scientific Reports，2017，7（1）：3788.

［23］PANAGIOTOU O A，EVANGELOU E，IOANNIDIS J P A. Genome‐wide significant associations for variants with minor allele frequency of 5% or less：An overview：A HuGE review［J］. American Journal of Epidemiology，2010，172（8）：869‐889.

［24］LIPKA A E，TIAN F，WANG Q S，et al. GAPIT：Genome association and prediction integrated tool［J］.Bioinformatics，2012，28（18）：2397‐2399.

［25］BRADBURY P J，ZHANG Z W，KROON D E，et al.TASSEL：Software for association mapping of complex traits in diverse samples［J］. Bioinformatics，2007，23 （19）：2633‐2635.

［26］RAJ A， STEPHENS M， PRITCHARD J K.fastSTRUCTURE：Variational inference of population structure in large SNP data sets［J］. Genetics，2014，197（2）：573‐589.

［27］VANRADEN P M. Efficient methods to compute genomic predictions［J］.Journal of Dairy Science，2008，91（11）：4414‐4423.

［28］YANG X，ZHONG S B，ZHANG Q J，et al. A loss‐of‐function of the dirigent gene TaDIR‑B1 improves resistance to Fusarium crown rot in wheat［J］. Plant Biotechnology Journal，2021，19（5）：866‐868.

［29］SUN C W，ZHANG F Y，YAN X F，et al. Genome‐wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China［J］. Plant Biotechnology Journal，2017，15（8）：953‐969.

［30］WANG H，JIA J，CAI Z D，et al. Identification of quantitative trait loci（QTLs）and candidate genes of seed iron and zinc content in soybean［Glycine max（L.）Merr. ］［J］. BMC Genomics，2022，23（1）：146.

［31］LIU H X，ZHANG X，SHANG Y P，et al. Genome‐wide association study reveals genetic loci for ten trace elements in foxtail millet（Setaria italica）［J］.Theoretical and Applied Genetics，2024，137（8）：186.

［32］MOUNTJOY E，SCHMIDT E M，CARMONA M，et al.An open approach to systematically prioritize causal variants and genes at all published human GWAS trait‐associated loci［J］.Nature Genetics，2021，53 （11）：1527‐1533.

［33］ZHAO Z，ZHANG H S，FU Z J，et al. Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method［J］. Plant Biotechnology Journal，2018，16（5）：1085‐1093.

［34］ARSHAD K，WANG Y，HAN S Q，et al. GWAS on HTP‐enabled dynamic traits unravels novel genetic architecture of salt tolerance in soybean［J］. The Plant Journal，2025，122（4）：e70177.

［35］LIU J，HUANG L，WANG C Q，et al. Genome‐wide association study reveals novel genomic regions associated with high grain protein content in wheat lines derived from wild emmer wheat［J］. Frontiers in Plant Science，2019，10：464.

［36］LIU X L，HUANG M，FAN B，et al. Iterative usage of fixed and random effect models for powerful and efficient genome‐wide association studies［J］. PLoS Genetics，2016，12（2）：e1005767.

［37］ZHOU Z F，SHI X，ZHAO G Q，et al. Identification of novel genomic regions and superior alleles associated with Zn accumulation in wheat using a genome‐wide association analysis method［J］. International Journal of Molecular Sciences，2020，21（6）：1928.

［38］WANG J S，SHI X，ZHOU Z F，et al. Genetic dissection of grain iron concentration in hexaploid wheat（Triticum aestivum L.）using a genome‐wide association analysis method［J］. PeerJ，2022，10：e13625.
［39］SANDHU G，KHAN A，TRIVEDI P K. Transport channels enabling uptake， translocation and detoxification of arsenic in plants［J］.Plant Physiology and Biochemistry，2025，225：109994.
［40］AMNA S，QAMAR S，ABU TURAB NAQVI A，et al.Role of sulfur in combating arsenic stress through upregulation of important proteins，and in‐silico analysis to study the interaction between phosphate transporter（PHO1），arsenic and phosphate in spinach［J］. Plant Physiology and Biochemistry，2020，157：348‐358.

［41］MAITY J P，CHEN C Y，BHATTACHARYA P，et al.Advanced application of nano‐technological and biological processes as well as mitigation options for arsenic removal［J］. Journal of Hazardous Materials，2021，405：123885.
［42］CLEMENS S，MA J F. Toxic heavy metal and metalloid accumulation in crop plants and foods［J］. Annual Review of Plant Biology，2016，67：489‐512.
［43］GIORDANO D，BLANDINO M. Arsenic，lead and cadmium distribution in the pearled fractions of different winter wheat cultivars（Triticum aestivum L.）［J］. Journal of Cereal Science，2018，80：94‐101.
［44］SHI G L，LI D J，WANG Y F，et al. Accumulation and distribution of arsenic and cadmium in winter wheat（Triticum aestivum L.） at different developmental stages［J］. Science of the Total Environment，2019，667：532‐539.

[1]	张素瑜, 刘红杰, 任德超, 杨明达, 葛君, 赵敬领, 朱培培, 郑东方. 外源调节剂对小麦强、弱势粒灌浆及产量的影响[J]. 河南农业科学, 2025, 54(9): 23-33.
[2]	王皓东, 逯一方, 杨珊, 刘凤英, 王刚. 绿针假单胞菌gacS 基因对小麦全蚀病生防能力的影响[J]. 河南农业科学, 2025, 54(9): 112-118.
[3]	刘星麟, 刘园, 杨凡, 刘布春, 韩锐. 基于机器学习算法的河南省县级冬小麦产量模拟研究[J]. 河南农业科学, 2025, 54(8): 167-180.
[4]	马尚英, 夏婷婷, 韩鹏彬, 张梦姣, 毛应杰, 王志强, 辛泽毓, 林同保, 连延浩, 任永哲. 施氮量和施氮方式对宽窄行种植冬小麦氮素吸收利用、产量及土壤氮残留的影响[J]. 河南农业科学, 2025, 54(8): 38-50.
[5]	黄世鹏, 何虎强, 李光照, 陈传鑫, 钟拓, 谢剑波, 周湘晚, 王丹, 肖应辉. 水稻苗期耐旱性全基因组关联分析及候选基因预测[J]. 河南农业科学, 2025, 54(7): 12-20.
[6]	臧贺藏, 周萌, 王亚辉, 彭一龙, 赵晴, 张杰, 李国强. 基于改进YOLOv8n 的麦穗自动检测方法研究[J]. 河南农业科学, 2025, 54(7): 162-169.
[7]	王海洋, 靳海洋, 宋航, 潘秀燕, 闫雅倩, 杨习文, 曾昭海, 臧华栋, 郑念, 李向东, 贺德先. 不同茬口对冬小麦干物质积累和转运、氮素吸收利用及产量的影响[J]. 河南农业科学, 2025, 54(6): 1-10.
[8]	伍革民, 韩雪, 吴松成, 冯宇隆. 三穗鸭体尺与屠宰性状的全基因组关联分析[J]. 河南农业科学, 2025, 54(6): 135-143.
[9]	李世行, 高巍, 李洋洋, 孙世龙, 班佳浩, 聂兆君, 耿月华, 汪敏, 赵鹏. 高去除性耐镉真菌的筛选及其对小麦幼苗生长和镉积累的影响[J]. 河南农业科学, 2025, 54(5): 59-67.
[10]	李晓德, 柏心彤, 吴恩召, 吕春蕾, 陈树林, 殷贵鸿, 王道文, 张坤普. 小麦骨干亲本周8425B 衍生品种（系）春化和光周期基因的组成分析[J]. 河南农业科学, 2025, 54(5): 23-30.
[11]	赵立尚, 杨天聪, 张小美, 宋全昊, 马知遥, 冯伟, 段剑钊. 耕作与秸秆还田方式对冬小麦产量及N₂O排放的影响[J]. 河南农业科学, 2025, 54(4): 27-36.
[12]	张博文, 靳海洋, 许海霞, 郑飞, 何宁, 苏亚中, 赵国建, 程红建, 程西永, 李向东. 花后喷施氨基酸螯合锌对小麦籽粒淀粉和蛋白质合成积累的影响[J]. 河南农业科学, 2025, 54(3): 8-19.
[13]	高翠民, 田媛, 张瑞卿, 何方, 韩伟锋, 张运红, 潘晓莹, 杨永辉. 灌溉方式对小麦-花生轮作体系作物产量及水氮利用率的影响[J]. 河南农业科学, 2025, 54(3): 31-39.
[14]	张福彦, 孙丛苇, 陈晓杰, 王嘉欢, 程仲杰, 赵婉, 马旭辉, 张建伟, 范家霖. 小麦矮秆小粒突变体smk56表型分析与基因初步定位[J]. 河南农业科学, 2025, 54(2): 21-29.
[15]	韩锐锋, 任荣魁, 董向阳, 胡泽彬, 郭方君, 方甜甜, 马文豪, 张书红. 不同类型肥料施用深度和施用量对冬小麦生长发育的影响[J]. 河南农业科学, 2025, 54(2): 30-39.