《Theoretical and Applied Genetics》：From QTL to candidate genes: a data-driven approach to unravel the genetic architecture of yellow rust resistance in central European wheat
编辑推荐：
关键信息摘要：一种利用高分辨率SNP（Single-Nucleotide Polymorphism，单核苷酸多态性）信息、序列变异注释及来自三个群体的QTL（Quantitative Trait Loci，数量性状位点）信息的数据挖掘策略，促进了YR（Yell
关键信息摘要：一种利用高分辨率SNP（Single-Nucleotide Polymorphism，单核苷酸多态性）信息、序列变异注释及来自三个群体的QTL（Quantitative Trait Loci，数量性状位点）信息的数据挖掘策略，促进了YR（Yellow Rust，条锈病）抗性候选基因的高效提名，并得到了Yr27及功能注释的支持。全基因组关联研究（GWAS，Genome-wide Association Studies）在许多作物中已成为常规手段，但候选基因的优先级排序仍具挑战性。研究人员开发了一种新方法，利用GWAS识别环境特异性QTL，分析了分布于三个实验群体中的5,840份小麦基因型，包括5,243份单交杂交种和597份精英品系。对亲本基因型进行重测序（Resequencing）后，经过滤鉴定出超过640,000个SNP。杂交组合面板在19个环境中检测了对Puccinia striiformis f. sp. tritici（Pst，小麦条锈病菌）的感病性，该病原菌尤其自2012年以来引发了严重的YR流行病。不同环境中的QTL模式存在显著差异，反映了时空动态性。结合通过新型GWAS方法获得的高分辨率SNP信息、序列变异注释及来自不同群体的QTL信息，使得能够高效提名与中欧小麦特别相关的YR抗性位点的候选基因。所开发策略从不同实验群体中挖掘数据的效力得到了证实，已验证的抗性基因Yr27被鉴定为一个QTL区域的唯一候选基因。
研究背景与意义
条锈病由Puccinia striiformis f. sp. tritici（Pst）引起，过去十年在中欧引发多次流行，植物抗病基因与病原菌毒力进化的军备竞赛是农业核心挑战。全基因组关联研究（GWAS，Genome-wide Association Studies）虽已常规化，但候选基因优先级排序困难，且已有数百个QTL（Quantitative Trait Loci，数量性状位点）报道存在冗余。目前仅克隆了少数YR抗性（R）基因，且多数NLR（Nucleotide-binding site-Leucine-rich repeat，核苷酸结合位点-富含亮氨酸重复序列）类基因符合基因对基因模型易被克服，而成株期抗性（APR，Adult Plant Resistance）基因更具持久性。现有抗性基因图谱多基于平均表现或靶向RenSeq（Resistance gene enrichment sequencing，抗性基因富集测序），忽略了环境特异性及非NLR类基因，且基因型与环境（G×E）互作显著。为此，研究人员开展本研究，旨在利用重测序数据与新型GWAS框架解析中欧小麦YR抗性的时空遗传架构并高效提名候选基因，该论文发表于《Theoretical and Applied Genetics》。
主要关键技术方法
研究人员利用三个实验系列（Exp I、II、III）共计5,243份单交杂交种及597份精英冬小麦品系，在19个环境中进行YR严重度表型评价；对亲本进行全基因组重测序（WGS，Whole-Genome Sequencing）获超640,000个高质量SNP（Single-Nucleotide Polymorphism，单核苷酸多态性）；开发三种线性混合模型（LMM，Linear Mixed Model）分别检测加性、显性及标记-环境互作效应，采用Bonferroni校正；基于PLINK按LD（Linkage Disequilibrium，连锁不平衡）r2>0.3且间距<50 Mbp划定QTL，区分单环境（SEQ）与多环境（MEQ）QTL；整合中国春参考基因组v2.1注释、NLR-Annotator及ANNOVAR变异注释，通过跨群体共有显著关联与错义突变（Nonsynonymous SNPs）双重标准筛选候选基因，并利用近等基因系九年田间数据验证。
研究结果
Strong genotype-by-environment interactions
研究人员在19个环境评估5,840份基因型成株期YR严重度，基因组重复性多高于0.3，广义遗传力超0.7。环境间表型Pearson相关聚类分四组，同地点不同年份相关性低（0.05–0.65），表明基因型×年份及基因型×地点×年份互作是G×E重要组分。
QTL consistency was low across experimental series
跨环境检测QTL主效应发现Exp I显著MTA（Marker-Trait Associations，标记-性状关联）多于其他系列，联合检验加性与显性增加MTA数。仅Q2A_0.47_39.79在三个系列共有，加性与显性QTL重叠极少，Exp I显性QTL（32个）远多于加性（2个），低重叠源于G×E及群体遗传组成差异。
Substantial QTL-by-environment interactions
检测标记×环境互作发现Exp II、III显著交互MTA多于Exp I，加性×环境互作多于显性×环境。共界定246个QTL，208个SEQ与38个MEQ（多环境QTL），68.42% MEQ与既往研究共定位。MEQ解释表型方差（PVE，Phenotypic Variance Explained）显著高于SEQ，联合效应检出83.7% QTL，多数QTL显性度（d/|a|）为负值显示抗性完全显性或超显性。
QTL-by-year interactions indicate potential breakdown of resistance genes
Hadmersleben地点五年分析显示SEQ仅现于2016年，18个MEQ现于2012与2016年。仅Q2A_0.47_39.79在2012与2016显著，Q1A_3.92_10.73等仅2016检出。九年Avocet与21份载不同Yr基因基因型田间验证显示Yr5、Yr27等七年稳定高抗，Yr17、Yr18具年份敏感性，证实年份效应致抗性基因失效。
Identification of candidate genes
38个MEQ涵盖5,100个高置信基因，171个NLR基因。跨实验系列共有显著关联与非同义SNP筛选后，Q1A_3.92_10.73提名NLR基因TraesCS1A03G0032200；Q2B_112.61_207.39唯一候选为Yr27（TraesCS2B03G0442800）；Q6A_607.21_617.66提名两LRR域基因；Q2A_0.47_39.79含33候选含6个NLR，对应Yr17渗入区段。
讨论与结论总结
讨论指出需表征田间Pst种群以解析QTL×E，本研究LMM框架兼顾加性、显性及互作优于既往仅加性模型；严格候选基因筛选利用跨群体LD与错义突变降假阳性，成功锁定Yr27验证策略有效，但参考基因组Chinese Spring可能影响外源渗入区（如Aegilops ventricosa来源Yr17）变异检出；中欧育种池中Q2A（Yr17渗入）、Q2B（Yr27/Lr13）、Q6A（潜在APR位点）、Q1A（成株抗性）为稳定位点；杂交育种中部分显性抗性仅需固定一方亲本池即可实现基因累加（Stacking）。
结论部分译为：研究人员构建了含超5,000份小麦杂交种及597份品系、跨越19环境超30,000条YR评分及WGS数据集，满足解析G×E需求；开发了关联作图方法识别环境特异性QTL并解析中欧小麦YR抗性G×E互作，拓展复杂性状遗传解析方法学；QTL模式反映时空动态，九年21份Yr基因系佐证G×E；QTL显性效应较加性跨环境更稳定，杂交育种或可提高抗性稳定性；结合高分辨率SNP、变异注释及三群体QTL的数据挖掘策略高效提名候选基因，与克隆基因Yr27及功能注释重叠证实策略效力。
