中国农业科学 ›› 2023, Vol. 56 ›› Issue (8): 1429-1443.doi: 10.3864/j.issn.0578-1752.2023.08.001

• 作物遗传育种·种质资源·分子遗传学 • 上一篇    下一篇

结合基因关联和转录组分析鉴定小麦成株期抗条锈病位点YrZ501-2BL的候选基因

张旭1(), 韩金妤1(), 李晨晨1, 张丹丹1, 吴启蒙1, 刘胜杰1, 焦韩轩1, 黄硕1, 李春莲1, 王长发1, 曾庆东2, 康振生2, 韩德俊1(), 吴建辉1()   

  1. 1 西北农林科技大学农学院/旱区作物逆境生物学国家重点实验室,陕西杨凌 712100
    2 西北农林科技大学植物保护学院/旱区作物逆境生物学国家重点实验室,陕西杨凌 712100
  • 收稿日期:2022-12-08 接受日期:2023-01-19 出版日期:2023-04-16 发布日期:2023-04-23
  • 联系方式: 张旭,E-mail:zhxu@nwafu.edu.cn。韩金妤,E-mail:1021979347@qq.com。 张旭和韩金妤为同等贡献作者。
  • 基金资助:
    国家重点研发计划(2021YFD1200600); 青海省重点研发与转化计划(2022-NK-125)

Identification of Adult Plant Stripe Rust Resistance Candidate Genes of YrZ501-2BL by Gene Association and Transciptome Analysis in Wheat (Triticum aestivum L.)

ZHANG Xu1(), HAN JinYu1(), LI ChenChen1, ZHANG DanDan1, WU QiMeng1, LIU ShengJie1, JIAO HanXuan1, HUANG Shuo1, LI ChunLian1, WANG ChangFa1, ZENG QingDong2, KANG ZhenSheng2, HAN DeJun1(), WU JianHui1()   

  1. 1 College of Agronomy, Northwest A&F University/State Key Laboratory of Crop Stress Biology for Arid Areas, Yangling 712100, Shaanxi
    2 College of Plant Protection, Northwest A&F University/State Key Laboratory of Crop Stress Biology for Arid Areas, Yangling 712100, Shaanxi
  • Received:2022-12-08 Accepted:2023-01-19 Published:2023-04-16 Online:2023-04-23

摘要:

【目的】小麦条锈病是小麦的主要病害之一,每年都会对小麦产量安全造成严重危害,挖掘小麦抗条锈病基因,为小麦抗条锈病种质创新和揭示小麦抗条锈病遗传机制奠定基础。【方法】利用多组学手段结合全基因关联分析(GWAS)开展对小麦成株期抗条锈病性状的解析。首先对411份来自CIMMYT和ICARDA的春小麦进行全基因组关联分析,在小麦2BL染色体上定位到一个主效的成株期抗条锈病位点,并利用含有该位点的抗病材料Z501及感病亲本晋麦79的双亲群体进行连锁作图,成功验证了该位点抗性的稳定性,暂命名为YrZ501-2BL。在此基础上,通过基因注释、比较基因组分析、转录组分析和候选基因的关联分析对目标区间筛选候选基因。【结果】综合GWAS和连锁作图结果,将YrZ501-2BL锁定在小麦2B染色体0.26 Mb(575.706—576.587 Mb)范围内,根据中国春参考基因组注释信息分析,该区间含有12个基因,其中,高可信基因6个;利用在线网站,将目标区间所在的中国春参考基因组与其他已公布的不同倍性小麦基因组进行比较,发现该区间的6个高可信小麦基因基本都能在其他小麦材料中找到同源基因,且基因排列顺序相同,说明该区间可能不存在较大片段的插入、缺失、倒位等现象,可以利用参考基因组信息进行候选基因预测;结合抗病亲本Z501和感病亲本晋麦79的成株期接种条锈菌后的转录组数据,发现只有TraesCS2B02G406400TraesCS2B02G406500TraesCS2B02G406600受诱导表达,且抗感病亲本存在表达差异。根据中国春基因组注释,三者分别编码GATA转录因子、SH3P2蛋白和锌指蛋白。通过进一步的候选基因关联分析,发现只有SH3P2蛋白中存在与条锈病表型显著性差异的SNP位点(G/A),虽然该位点(G1369A)在2个可变剪切的转录本中均未引起氨基酸编码变化(TCG和TCA均编码丝氨酸),但可能与可变剪切有关,同时该位点(G1369A)的不同单倍型表型之间也存在极显著差异,进一步分析发现,在G1369A位点下游还有2个引起氨基酸改变的变异位点G1377A和G1431A,分别引起缬氨酸(GTT)到异亮氨酸(ATT)和缬氨酸(GTG)到甲硫氨酸(ATG)的改变,但这两个位点在455份重测序材料中所占比例只有0.87%,属于稀有变异,因此,未进行显著性检验。综上,推测TraesCS2B02G406500YrZ501的重要抗病候选基因;此外,针对YrZ501候选区间的差异SNP开发了相应的AQP标记,可用于辅助选择,为下一步小麦抗锈病分子育种应用提供了标记资源。【结论】利用多组学整合关联分析的方法,成功在小麦2B染色体上挖掘到1个抗条锈病候选基因TraesCS2B02G406500

关键词: 小麦, 条锈病抗性, YrZ501, 多组学整合分析, SH3P2蛋白

Abstract:

【Objective】Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), significantly reduced wheat production worldwide. Identification of stripe rust resistance genes is the foundation of improving wheat resistance breeding and revealing its genetic mechanism.【Method】A multi-omics approach combined with genome-wide association study (GWAS) was used for dissecting adult plant stripe rust resistance for wheat advanced breeding lines collected from International Maize and Wheat Improvement Center (CIMMYT) and International Centre for Agricultural Research in the Dry Areas (ICARDA) bread-wheat breeding programs. In the present study, a diversity panel of 411 wheat lines from CIMMYT and ICARDA was used for genome-wide association study and a major locus on chromosome arm 2BL was identified. In order to verify the stability of the locus, the resistant line Z501 with the resistance allele of the locus was crossed by the susceptible line Jinmai 79, and the locus tentatively named YrZ501 was successfully confirmed using linkage mapping based on F2:3 genetic population of Jinmai 79×Z501. Then we performed candidate gene analysis based on gene annotation, comparative genome, transcriptome and gene-based association analysis. 【Result】Combining GWAS and linkage mapping results, the YrZ501-2BL was located in the physical interval of 0.26 Mb (575.706-576.587 Mb) on chromosome 2B. According to the annotation information of Chinese Spring reference genome IWGSC v1.1, there were six high confidence genes of 12 genes in this region. Using online website, the target interval in the Chinese spring reference genome was compared with other published different ploidy wheat genomes. The six high-confidence genes within this interval can basically be found homologous in other wheat lines, and the genes arranged in the same order, indicating that the interval may not have large fragment insertions, deletions and inversions. The above results showed that we can perform candidate gene prediction analysis based on the reference genome information. After analysis of their transcriptomic data between the resistant parent Z501 and susceptible parent Jimai 79, only three genes, TraesCS2B02G406400, TraesCS2B02G406500 and TraesCS2B02G406600 showed variable expression levels and were induced by stripe rust infection. Further, they encode GATA transcription factor, SH3 domain-containing protein 2 and zinc finger protein, respectively. Gene-based association analysis revealed that there was a significant SNP (G1369A) in TraesCS2B02G406500 that was associated with stripe rust responses. Although this SNP (G1369A) did not cause amino acid coding changes (both TCG and TCA encode serine), it may be associated with alternative splicing. Moreover, it showed significant differences of the stripe rust responses between the different haplotypes (G1369A). Further analysis revealed two other variants G1377A and G1431A, that caused amino acid changes, i. e. valine (GTT) to isoleucine (ATT) and valine (GTG) to methionine (ATG), respectively. However, the two SNPs were rare variants as they accounting for only 0.87% of the 455 re-sequencing wheat accessions and they were not tested for significance. In summary, TraesCS2B02G406500 was preliminarily considered as an important candidate gene of YrZ501-2BL. In addition, the corresponding AQP markers were developed based on the SNPs among the YrZ501 candidate regions, which can be used to marker-assisted selection in molecular breeding application of wheat stripe rust resistance.【Conclusion】A candidate causal gene TraesCS2B02G406500 associated with stripe rust resistance was successfully identified on wheat chromosome 2B using an integrated method of multi-omics and association analysis, which laid a solid foundation for further gene cloning and functional verification.

Key words: Triticum aestivum, stripe rust resistance, YrZ501, integrative analysis of multiple omics, SH3P2