《Tree Genetics & Genomes》：Genome-wide association analysis of ecologically essential traits in Picea crassifolia
編輯推薦：
為解析青海云杉(Picea crassifolia)關鍵性狀遺傳基礎、克服其傳統育種周期過長瓶頸，本文開展了木材、針葉及生長性狀的全基因組關聯分析(GWAS)。基于祁連山106個無性系的研究發現了與9個性狀顯著關聯的84個SNP標記，鑒定出10個候選基因，部分基因涉及植物形態調控，為未來分子育種提供了寶貴的遺傳工具。
在廣袤的祁連山脈，生長著一種名為青海云杉(Picea crassifolia)的獨特常綠樹種。作為中國的特有物種，它不僅是一道靚麗的風景線，更在調節區域生態平衡中扮演著至關重要的角色。然而，面對氣候變化和人類活動的影響，如何提升這種樹種的經濟與生態價值，成為了科研人員亟待解決的難題。傳統的育種方法，如優選單株雜交，對于青海云杉這類世代周期長、遺傳雜合度高的樹種來說，效率低下，進展緩慢。核心問題在于，我們對其關鍵性狀（如木材密度、針葉形態、生長速度）的遺傳基礎知之甚少。盡管已有研究關注其種群結構和環境適應性，但將木材、針葉和生長等多類性狀整合在一起的系統遺傳研究仍屬空白。這種認知上的局限，嚴重阻礙了青海云杉種子園的建設和分子育種的發展。因此，深入探索青海云杉關鍵性狀的遺傳變異及其背后的分子機制，成為解鎖其育種潛力、應對環境挑戰的關鍵。
為了回答這些問題，一篇發表于《Tree Genetics》的研究《Genome-wide association analysis of ecologically essential traits in Picea crassifolia》應運而生。該研究團隊從祁連山區的核心自然分布區收集了106個青海云杉無性系作為研究材料，并建立了一個無性系種子園。研究旨在通過全基因組關聯分析(GWAS)，系統評估木材、針葉及生長這三類重要性狀，以期闡明其遺傳基礎。研究團隊測量了17項木材性狀（如管胞長度、管胞腔徑、晚材率）、7項針葉性狀（如針葉長度、寬度、干重）和4項生長性狀（樹高、地徑、冠幅、年高生長量），利用基于特異性位點擴增片段測序(SLAF-seq)技術開發了全基因組水平的單核苷酸多態性(SNP)標記，并運用一般線性模型(GLM)和混合線性模型(MLM)進行關聯分析，最終篩選出與性狀顯著關聯的SNP標記并挖掘了鄰近的候選基因。
木材和針葉性狀在青海云杉中表現出高表型變異
研究發現，青海云杉無性系在木材和針葉性狀上存在顯著的表型變異，所有性狀的平均變異系數(CV)為20.79%。其中，木材性狀的平均CV（22.97%）高于針葉性狀（15.49%），表明木材性狀可能受環境影響更大。例如，早材管胞壁腔比的變異系數高達89.47%，而針葉含水量的變異系數最低，為5.52%。對性狀頻率分布的分析表明，大多數性狀符合單峰正態分布，表明它們是典型的數量性狀，受多基因控制，適合進行GWAS分析。相關性分析顯示，大多數木材性狀之間以及大多數針葉性狀之間均存在顯著（P< 0.05）或極顯著（P< 0.01）的相關性，而大多數木材與針葉性狀之間的相關性通常不顯著，暗示這兩類性狀可能由不同的遺傳機制控制。例如，針葉厚度與鮮重、干重、長寬比及含水量呈顯著正相關。
性狀與SNP的關聯分析
通過GWAS，研究團隊使用GLM方法識別出84個SNP分子標記與9個表型性狀在P< 0.05水平上顯著關聯，而MLM方法未檢測到顯著關聯。這些SNP位點對表型變異的解釋率（即R2值）介于24.34%至50.30%之間。其中，13個SNP標記與兩個或更多性狀相關聯，這可能是性狀間表型相關性的遺傳基礎，反映了基因的多效性。特別值得一提的是，有18個SNP標記在GLM和MLM兩種方法中均被一致檢測到，這增加了這些關聯的可靠性。
候選基因的鑒定
為了深入理解這些遺傳標記的功能，研究在顯著關聯的SNP位點上下游100千堿基對(kb)區域內進行了候選基因挖掘，共鑒定出10個候選基因。例如，與生長輪寬度顯著關聯的基因MA_35620g0010被預測編碼生長激素調控的TBC蛋白1；與晚材寬度關聯的基因MA_10430313g0010編碼翻譯起始因子；與針葉鮮重、干重和厚度關聯的基因MA_92129g0010則編碼一個參與植物細胞器基因表達調控的五肽重復序列(PPR)蛋白。這些基因的功能預測為理解木材形成和針葉發育的分子調控網絡提供了線索。
生長性狀的關聯分析
同樣，對生長性狀的關聯分析也取得了豐碩成果。應用GLM方法，共發現109、23和56個SNP標記分別與15年和20年生的樹高(TH)、地徑(GD)和冠幅(CW)的克隆平均值顯著相關。MLM方法則識別出較少但更為保守的關聯標記（TH：18個，GD：8個，CW：9個）。有16個SNP標記與多個生長性狀相關，且在不同年份間被穩定檢測到。最終，10個相關的候選基因被鑒定出來，其中包括預測編碼BolA-like蛋白的MA_129402g0010，該蛋白參與細胞應激反應，可能對生長調控有所貢獻。
綜上所述，本研究揭示了青海云杉在木材、針葉和生長性狀上存在豐富的表型變異，并利用GWAS成功定位了大量與其顯著關聯的SNP標記和候選基因。這些發現為青海云杉的遺傳改良提供了寶貴的分子工具和理論基礎。盡管研究地點（海拔1700米的半干旱氣候）可能篩選出與干旱適應性相關的性狀標記，且研究僅評估了15年和20年生的性狀，未能涵蓋成熟后期（如心材形成）的轉變，但本研究的成果無疑為未來的分子標記輔助選擇育種鋪平了道路。識別出的SNP標記和候選基因，特別是那些涉及植物形態調控的基因，可作為未來育種計劃中有價值的工具，有望顯著縮短育種周期，提高青海云杉的遺傳改良效率，對于培育適應性強、生態與經濟價值更高的青海云杉品種具有重要意義。未來的研究需要在不同生態區（如高海拔與低地）的后代試驗中驗證這些SNP的穩定性，并評估30齡以后性狀與SNP關聯的持久性，以推動研究成果的實際應用。
