《Plant Stress》：Identification of genetic loci and candidate genes for salt tolerance in barley (
Hordeum vulgare L.) at the germination stage using a genome-wide association study and multi-omics analysis
編輯推薦：
為解決土壤鹽漬化嚴重制約作物生產的問題，本研究通過整合全基因組關聯(lián)研究(GWAS)、轉錄組和代謝組學，深入解析了大麥(Hordeum vulgare L.)萌芽期耐鹽性的復雜調控網(wǎng)絡。研究在250份大麥種質中鑒定出306個顯著SNPs，揭示了鹽脅迫響應中氨基酸與苯丙烷類生物合成等關鍵通路，并篩選出24個與抗氧化酶、膜蛋白轉運、轉錄因子等相關的候選基因，為培育耐鹽大麥品種提供了寶貴的遺傳資源和理論依據(jù)。
在全球范圍內，有近十億公頃的土地受到鹽漬化的威脅，不合理的灌溉和劇烈的氣候變化更使鹽漬地面積持續(xù)擴張。土壤鹽分不僅直接抑制植物生長，更對全球糧食安全構成了嚴峻挑戰(zhàn)。大麥，作為全球第四大谷類作物，因其相對較強的耐鹽性而被視為開發(fā)利用鹽堿地的“先鋒作物”。然而，植物在種子萌發(fā)階段對鹽脅迫最為敏感，而作物的耐鹽性又是一個由眾多基因相互作用的復雜數(shù)量性狀，其背后的遺傳基礎和分子機制仍有待深入解析。為了挖掘大麥萌芽期耐鹽性的關鍵基因資源，并為分子育種提供靶點，這項研究應運而生。
為了系統(tǒng)解答上述問題，研究人員開展了一項綜合性研究，成果發(fā)表在期刊《Plant Stress》上。該研究通過對一個包含250份不同來源大麥種質的自然群體，在連續(xù)兩年內施加200 mM NaCl鹽脅迫，評估了相對發(fā)芽率、相對根/芽長、相對根/芽鮮重等五個關鍵形態(tài)性狀。隨后，整合了全基因組關聯(lián)研究(GWAS)、非靶向代謝組學(基于LC-MS/MS)和轉錄組學(RNA-seq)分析，并結合實時定量PCR(RT-qPCR)驗證，系統(tǒng)鑒定了與耐鹽性相關的遺傳位點、差異代謝物、差異表達基因及核心候選基因，最終提出了大麥萌芽期響應鹽脅迫的潛在分子機制模型。
本研究主要運用了以下幾項關鍵技術方法：首先，基于250份大麥自然種質（部分材料來自中國江蘇沿海地區(qū)農業(yè)科學研究所）的表型數(shù)據(jù)，利用混合線性模型(MLM)進行了全基因組關聯(lián)研究(GWAS)。其次，以耐鹽品種“Liu Leng Zi Da Mai”(LLZDM)的根組織為材料，通過液相色譜-串聯(lián)質譜(LC-MS/MS)進行非靶向代謝組學分析。同時，對鹽處理和對照條件下的LLZDM根組織進行了轉錄組測序(RNA-seq)。最后，通過整合GWAS和轉錄組數(shù)據(jù)篩選候選基因，并使用RT-qPCR在多個品種中對部分基因的表達進行了驗證。
研究結果如下：
3.1. 表型變化分析與性狀間相關性
鹽處理顯著抑制了所有大麥基因型的芽和根生長，各性狀的相對值均顯著降低。250份大麥種質在耐鹽性上表現(xiàn)出顯著的遺傳變異。其中，地方品種如“Liu Leng Zi Da Mai”(LLZDM)比栽培品種如“Nasi Nijo”(NN)表現(xiàn)出更強的耐鹽性。在正常條件下兩者無顯著差異，但在鹽脅迫下，NN的性狀下降幅度(50.00%至78.03%)遠大于LLZDM(21.16%至39.93%)，清晰證明了LLZDM的耐鹽優(yōu)勢。
3.2. 萌芽期耐鹽性的全基因組關聯(lián)分析
利用106,131個高質量SNP進行GWAS分析，以-log10P≥ 4.0為閾值，在兩年間共檢測到306個與五個性狀顯著相關的SNP。其中，有23個SNP被不同性狀同時檢測到，8個SNP在兩年中均被穩(wěn)定檢測到，這些穩(wěn)定SNP主要定位于2號和6號染色體。單個SNP所能解釋的表型變異率(PVE)范圍為0.09%至3.91%。
3.3. 大麥根響應鹽脅迫的代謝組學分析
對耐鹽品種LLZDM根組織的非靶向代謝組學分析表明，鹽脅迫誘導了324個差異積累代謝物(DAMs)，其中167個上調，157個下調。這些代謝物主要屬于脂質和類脂分子、有機酸及其衍生物等類別。京都基因與基因組百科全書(KEGG)通路富集分析顯示，DAMs顯著富集于生物素代謝、氨基糖和核苷酸糖代謝等通路。
3.4. 鹽脅迫條件下大麥根的轉錄組學分析概覽
RNA-seq分析在LLZDM根組織中鑒定出2351個差異表達基因(DEGs)，其中1544個上調，807個下調。基因本體(GO)富集分析顯示，DEGs主要富集于離子轉運、對激素響應、脅迫響應等功能類別。KEGG通路分析進一步表明，DEGs顯著富集于苯丙烷生物合成、丙氨酸、天冬氨酸和谷氨酸代謝、谷胱甘肽代謝等通路。
3.5. 代謝組與轉錄組的整合分析
整合分析發(fā)現(xiàn)，DEGs和DAMs共同顯著富集于“氨基酸的生物合成”和“苯丙烷生物合成”等關鍵通路。在氨基酸生物合成通路中，鹽脅迫導致谷氨酸、脯氨酸等含量顯著上升，而賴氨酸、組氨酸等含量下降，同時有12個相關基因上調。在苯丙烷生物合成通路中，鹽脅迫誘導了咖啡酰輔酶A、圣草酚等代謝物的積累，并上調了包括5個肉桂酰輔酶A還原酶(COMT)基因在內的多個關鍵基因的表達。
3.6. 候選基因的鑒定與驗證
通過整合GWAS和RNA-seq數(shù)據(jù)，最終鑒定出24個與耐鹽性相關的候選基因。這些基因在染色體上分布不均，其中17個上調，7個下調。根據(jù)功能可分為幾類：與抗氧化酶相關，如過氧化氫酶(CAT)、谷胱甘肽S-轉移酶(GST)；與膜蛋白轉運相關，如跨膜蛋白184C(TP184C)、鉀轉運蛋白(KT)；轉錄因子，如一個WRKY轉錄因子顯著上調；以及受體蛋白，如受體樣蛋白激酶(RLKs)。RT-qPCR對其中12個基因的表達驗證結果與轉錄組數(shù)據(jù)基本一致，證實了篩選結果的可靠性。
歸納研究結論與討論：
本研究通過多組學整合策略，系統(tǒng)揭示了大麥萌芽期耐鹽性的復雜調控網(wǎng)絡。結論表明，鹽脅迫下大麥通過協(xié)調氨基酸代謝和苯丙烷類代謝來應對氧化脅迫和離子失衡。其中，脯氨酸、谷氨酸等滲透調節(jié)物質的積累，以及苯丙烷類途徑中COMT等基因的上調和類黃酮等抗氧化物質的積累，是重要的耐受機制。
更重要的是，研究篩選出的24個核心候選基因，包括GST、TP184C、KT、WRKY和RLKs等，分別從抗氧化防御、離子穩(wěn)態(tài)、脅迫信號感知與轉導等不同層面發(fā)揮作用。基于這些發(fā)現(xiàn)，論文提出了一個耐鹽機制模型：鹽脅迫信號首先被細胞膜上的RLKs等受體感知，進而激活WRKY等轉錄因子，后者調控下游靶基因（如抗氧化酶基因、離子轉運蛋白基因等）的表達，最終通過增強ROS清除能力、維持離子平衡和促進保護性代謝物合成，來改善大麥在鹽脅迫下的發(fā)芽和幼苗生長。
這項研究的意義重大。它不僅為理解大麥萌芽期耐鹽性的分子機制提供了新的見解，更重要的是，所鑒定的306個顯著SNP和24個候選基因，特別是那些在兩年中穩(wěn)定出現(xiàn)的位點，為后續(xù)通過基因組輔助選擇培育耐鹽大麥品種提供了寶貴的遺傳資源和分子靶點。盡管基于有限樣本的轉錄組分析可能遺漏其他種質中的有利等位基因，且單個標記的解釋率有限，但將多個位點的優(yōu)良等位基因聚合，有望顯著提升品種的耐鹽性。未來，通過增加測序深度、擴大群體規(guī)模以及利用基因編輯等技術對候選基因進行功能驗證，將進一步推動大麥耐鹽分子育種的發(fā)展，對于利用廣大鹽堿地、保障糧食安全具有重要的實踐價值。
