《Theoretical and Applied Genetics》：Genome-wide association study uncovers the genetic basis of stem rust resistance in the southern US oat (Avena sativa L.) germplasm at different growth stages
編輯推薦：
本研究關注由Puccinia graminis f. sp. avenae引起的燕麥莖銹病(SR)對全球燕麥生產的嚴重威脅。為支持美國東南部的抗性育種，研究人員利用包含440個品系的美國南部燕麥關聯分析面板(SOAP)，進行了多環境、多年份的全基因組關聯研究(GWAS)，在成株期和苗期共檢測到17個數量性狀位點(QTL)，其中染色體1A的Pg13區間信號最為顯著。效應呈加性，堆疊5個穩定的成株期QTL可使病害嚴重度(SV)降低近一半。該研究不僅定位了6個已知的抗性區域，還新報告了11個位點，為燕麥莖銹病的分子標記輔助選擇(MAS)和抗性基因的持久性利用提供了關鍵的遺傳資源和理論依據。
燕麥是一種重要的全球性谷物，不僅為人類和動物提供食物，其富含的蛋白質、β-葡聚糖和植物化學物質還具有促進心血管和代謝健康的價值。然而，在過去幾十年里，美國的燕麥種植面積急劇萎縮，這為滿足國內市場對高品質燕麥的需求帶來了挑戰。開發結合優良品質和對關鍵生產限制因素具有抗性的新品種，是扭轉這一趨勢的關鍵。莖銹病(SR)正是制約燕麥生產的主要病害之一，其病原菌Puccinia graminis f. sp. avenae能夠侵染燕麥的莖稈、葉鞘和穗，導致光合作用減弱、籽粒灌漿受損，在嚴重流行年份可造成高達85%的產量損失。歷史上，燕麥對莖銹病的抗性主要由稱為Pg基因的主效基因提供，這類基因通常表現為小種專化性抗性，但易受病原菌小種組成變化的影響而失效。例如，Pg2、Pg9、Pg13等基因曾在北美燕麥品種中被廣泛應用，但新的病原菌小種（如TJJ、TGN）的出現，已能夠克服這些抗性基因，導致抗性喪失。相較于小種�；钥剐裕芍昶诓糠挚剐酝ǔ８鼮槌志�，但選擇起來也更困難。通過聚合多個有效抗性位點來延長抗性的持久性，是當前抗病育種的重要策略，但這一進程一直受到缺乏緊密連鎖標記以及能夠區分目標基因的病原菌分離物的限制。為了系統解析美國南部燕麥種質資源中莖銹病抗性的遺傳基礎，并支持該地區的抗性育種項目，本研究應運而生。
研究人員開展這項研究運用了幾個關鍵技術方法。首先，利用名為“美國南部燕麥關聯分析面板(SOAP)”的包含440個來自多個美國南部公立育種項目的精英品系、品種和歷史重要材料的群體，在2022-2024年的多個地點（如佛羅里達州的Citra、德克薩斯州的Castroville、路易斯安那州的Baton Rouge和Winnsboro）進行了成株期莖銹病的多年多點田間表型鑒定，測量了病害嚴重度(SV)和侵染反應(IR)。其次，在苗期，利用四種病原菌小種(DBD, SGD, TGN, TJS)的分離物，對SOAP群體進行了苗期抗性評估，以感染類型(IT)作為表型指標。在基因型分析方面，研究采用簡化基因組測序(GBS)技術對群體進行了基因分型，獲得了12914個高質量的單核苷酸多態性(SNP)標記用于后續分析。最后，通過全基因組關聯分析(GWAS)，結合混合線性模型(MLM)校正群體結構和親緣關系，系統檢測了與苗期和成株期莖銹病抗性相關的基因組區域。
Phenotypic analysis and heritability (表型分析與遺傳力)
研究在成株期觀察到了莖銹病抗性性狀的顯著變異。在大多數環境中，基因型效應高度顯著，成株期病害嚴重度(SV)和侵染反應(IR)的廣義遺傳力(H2)為中等至高，表明這些性狀受較強的遺傳控制。結合七個環境的數據分析顯示，SV和IR的H2分別為0.87和0.81，而狹義遺傳力(h2)分別為0.79和0.70，這表明成株期莖銹病抗性的遺傳方差有很大比例可歸因于加性效應。在苗期，SOAP群體對四個病原菌小種表現出廣泛的感染類型反應。對于DBD和SGD小種，IT值覆蓋了全范圍，存在抗、感和中間型反應。而對于TGN和TJS小種，平均IT值很高，反應以感病為主，表明多數品系對這些當前流行的小種高度感病。
Population structure analysis and marker coverage (群體結構分析與標記覆蓋度)
主成分分析(PCA)揭示了燕麥群體中存在微弱但可檢測的群體分層，可分為三個亞群。前兩個主成分分別解釋了總遺傳變異的6.3%和4.9%。群體結構與育種地理來源相關，例如來自德克薩斯的品系形成了一個獨特的簇，而來自北卡羅來納州的品系則單獨聚類。結構分析也支持將群體劃分為三個組。對全基因組12914個SNP標記的分析顯示，標記在21條染色體上分布不均勻，其中C亞基因組的標記密度最高，其次是A和D亞基因組。
Linkage disequilibrium (連鎖不平衡)
不同亞基因組間的連鎖不平衡(LD)衰減模式差異顯著。A和D亞基因組的LD衰減最快，分別在5.37 Mb和4.59 Mb處降至r2= 0.2的閾值。相比之下，C亞基因組的LD衰減非常緩慢，持續至51.42 Mb。全基因組水平的LD衰減距離為13.35 Mb，這凸顯了C亞基因組相對于A和D亞基因組具有更長的LD。
Genome-wide association study (GWAS) (全基因組關聯研究)
通過GWAS，研究共檢測到475個超過顯著性閾值的SNP，通過獨立性檢驗，這些信號可歸納為17個不同的QTL。其中，10個在成株期特異性表達，5個在苗期特異性表達，2個在兩個階段均有表達。在染色體1A上鑒定出一個由四個峰值組成的簇，它們都位于已知的Pg13基因區間(370.8–464.4 Mb)內。條件分析表明，這些峰值不能用單一效應解釋，表明該區間內存在多個獨立的決定因子。其中兩個位點同時影響苗期對TGN小種的IT反應和成株期田間反應，這與至少存在一個全生育期小種�；越M分相一致，而另外兩個位點則僅限于成株期。此外，在染色體1D和6D上還檢測到對成株期田間抗性有加性貢獻的較小效應位點。相比之下，在染色體2A、2D和3C上檢測到了小種專化性的苗期抗性位點。
Post-hoc evaluation of additivity (加性效應的事后評估)
對五個在至少兩個環境中穩定檢測到的成株期QTL進行的等位基因堆疊分析顯示，抗性效應具有累積性。當攜帶三個或更多有利等位基因時，品系的病害嚴重度(SV)相較于不攜帶任何有利等位基因的品系顯著降低（降低25-32%），而攜帶五個有利等位基因可使SV降低43.6%。對于侵染反應(IR)，也觀察到了類似的顯著降低模式。這些結果表明，有意義的成株期莖銹病抗性增益主要是通過聚合多個有利等位基因實現的。
Candidate gene identification (候選基因鑒定)
在顯著SNP兩側±1 Mb的保守區間內，研究鑒定出多個與植物免疫相關的候選基因。在成株期，候選基因富集了經典的抗性相關家族，例如染色體1A上的F-box/LRR重復蛋白和酪蛋白激酶樣蛋白，染色體1D上的疾病抗性蛋白簇(RGA5, Pik-2, At1g63350)和受體樣胞質激酶，以及染色體6D上的富含半胱氨酸的受體樣蛋白激酶和MAPK8。在苗期，在特定小種相關的位點附近也鑒定出強候選基因，如染色體2A上的多個RPS2樣抗性蛋白和可能的LRR受體樣絲氨酸/蘇氨酸蛋白激酶，染色體2D上的FLS2樣激酶，以及染色體3C上的凝集素受體激酶。這些基因在病原菌識別和下游防御信號傳導中扮演關鍵角色。
綜上所述，本研究通過對美國南部燕麥種質資源進行系統的GWAS分析，全面揭示了苗期和成株期莖銹病抗性的復雜遺傳結構。研究成功鑒定出17個QTL，其中染色體1A的Pg13區間是主導信號區域，且該區域內存在多個獨立的抗性因子。研究證實了抗性效應的加性特性，并證明通過聚合多個有利等位基因可顯著增強成株期抗性水平。此外，研究還定位了多個新的抗性位點，并鑒定了富含NBS-LRR和受體樣激酶等典型防御相關基因家族的候選區間。這些發現不僅深化了對燕麥-莖銹菌互作遺傳機制的理解，更重要的是，為美國東南部燕麥育種項目提供了可直接用于標記輔助選擇(MAS)的分子標記和候選基因資源。通過利用這些新發現的遺傳位點，育種家可以更高效地聚合抗性基因，培育出具有更廣譜、更持久莖銹病抗性的燕麥新品種，這對于應對不斷演變的病原菌威脅、穩定和提高燕麥產量、保障全球燕麥生產的可持續性具有重要的理論和實踐意義。
