《The Plant Genome》：Leaf transcriptome differences between diploid and tetraploid bahiagrass
編輯推薦：
本研究通過對亞熱帶牧草巴哈雀稗（Paspalum notatum）的二倍體與四倍體性系植株進行葉片轉錄組比較，系統揭示了全基因組復制（WGD）引起的廣泛轉錄組重塑。研究發現，多倍體化影響了發育、氧化還原穩態、光合作用等關鍵通路，鑒定出大量差異表達轉錄本（DET），包括轉錄因子（TF）和激素相關基因等關鍵育種靶點。共表達網絡分析進一步鑒定出532個受基因組加倍影響的主調控基因（switch genes），并發現其在基因組特定區域（如5號、10號染色體）存在非隨機聚集的“熱點”。該研究為理解多倍體化的分子后果提供了新視角，所建立的轉錄本目錄及鑒定的候選基因，為旨在提高脅迫耐受性、活力等農藝性狀的分子育種項目提供了寶貴資源。
引言：多倍體與巴哈雀稗研究模型
多倍體，即全基因組復制（WGD）導致的多套染色體組現象，是植物進化中的常見事件，在增強遺傳多樣性和環境適應性方面扮演關鍵角色。亞熱帶牧草巴哈雀稗是研究多倍體、生殖模式與脅迫關聯的理想模型。自然界中，該物種的二倍體（2n = 2x = 20）基因型通過自交不親和的性生殖繁殖，而主要的四倍體（2n = 4x = 40）基因型則通過假配子無孢子生殖（apomixis）進行無性繁殖。盡管此前有研究關注花器官中與生殖相關的基因表達變化，但葉片對多倍體化的轉錄組響應尚未被專門探究。本研究旨在通過比較性系二倍體與四倍體巴哈雀稗的葉片轉錄組，排除無孢子生殖控制位點（ACL）的干擾，明確鑒定出與倍性水平（如脅迫耐受性、巨型效應）特異性相關的候選基因，為基于分子的育種策略提供支持。
材料與方法：轉錄組構建與分析流程
研究使用了來自NCBI數據庫的轉錄組測序原始數據，涉及兩個性系二倍體（2xsex）和一個性系四倍體（4xsex）基因型，每個基因型包含三個生物學重復。首先，利用已有的巴哈雀稗R1參考基因組進行參考轉錄本組裝，構建了二倍體和四倍體葉片轉錄組。通過主成分分析（PCA）和聚類分析驗證了樣本按倍性水平清晰分離。差異表達（DE）分析基于先前構建的參考葉片轉錄本組，使用Kallisto進行表達量定量，并采用DESeq2和limma-voom流程（后者結合duplicateCorrelation功能以考慮基因型內相關性）進行差異表達鑒定，設定錯誤發現率（FDR）< 0.05且|log2FC| > 2為標準。對鑒定出的差異表達轉錄本（DET）進行了基因本體（GO）富集分析、KEGG通路映射、轉錄因子（TF）和激素相關基因分類。此外，利用SWItch Miner（SWIM）工具構建基因共表達網絡，識別受基因組加倍影響的主調控基因（switch genes），并通過泊松富集檢驗分析其在基因組上的分布。
結果
1. 樣本聚類與轉錄組組裝
PCA和聚類分析顯示，所有生物學重復緊密聚集，并清晰地按倍性（二倍體 vs. 四倍體）分組。組裝得到的二倍體葉片轉錄本組包含76,681個轉錄本，四倍體轉錄本組包含69,996個轉錄本。BUSCO評估顯示轉錄本組完整性良好。
2. 二倍體與四倍體葉片中的差異表達轉錄本
初步DESeq2分析在二倍體與四倍體樣本間鑒定出40,844個具有顯著差異表達的轉錄本（DET），其中18,785個在二倍體中過表達，22,059個在四倍體中過表達。為驗證倍性效應并減少假陽性，后續采用包含更多基因型的分析以及limma-voom流程進行驗證，最終獲得28,054個重疊的“真實DET”用于后續分析。
3. 受倍性增加影響的分子通路
GO富集分析顯示，在細胞組分層面，最富集的類別包括質體膜、外膜、葉綠體類囊體等，表明多倍體化對光合裝置有顯著影響。在生物過程層面，受影響最大的通路包括蛋白質運輸、非編碼RNA加工、單羧酸代謝過程和葉片發育等。分子功能類別則富集了蛋白質復合物結合、ATP水解酶活性等。KEGG通路富集分析表明，卡爾文循環（M00165）、糖異生（M00003）、單木質醇生物合成（M00039）等中心碳水化合物代謝和木質素合成通路在多倍體植物中被顯著重組。MapMan分析進一步證實，全基因組復制（WGD）影響了蛋白質穩態、RNA生物合成、染色質組織、脂質代謝和光合作用等主要類別。與脅迫響應和信號傳導相關的類別，如植物激素作用、氧化還原穩態和DNA損傷反應也顯著富集。特別是在氧化還原穩態類別中，多個活性氧清除酶（如超氧化物歧化酶、過氧化氫酶、谷胱甘肽過氧化物酶等）的表達發生失調。
4. 二倍體與多倍體植物中差異表達的轉錄因子和激素相關轉錄本
轉錄因子分析顯示，受WGD影響最顯著的家族包括bHLH、bZIP、C2H2、MYB-related、NAC和WRKY等。其中，一些TF家族（如B3、Dof、TCP等）在四倍體中持續上調，而另一些（如DBB、YABBY）則在二倍體中過表達。激素相關轉錄本分析表明，脫落酸（ABA）、生長素（Auxin）和茉莉酸（JA）是受影響最顯著的類別。熱圖清晰展示了與赤霉素、茉莉酸、脫落酸、乙烯、細胞分裂素和生長素相關的多個轉錄本群的差異表達模式。
5. 共表達網絡分析
SWIM相關性網絡分析從2835個節點中鑒定出532個主調控基因（switch genes）。網絡穩健性通過留一法（LOO）和80%子抽樣分析得到證實。在這些主調控基因中，有17個編碼轉錄因子（如bHLH、AP2、NAC等家族成員），14個編碼蛋白激酶，10個與光合裝置特異性相關（如PsbP），35個被鑒定為脅迫響應基因（如谷胱甘肽S-轉移酶GST）。這支持了多倍體在脅迫耐受和穩態維持中的作用。
6. 巴哈雀稗基因組中的主調控基因熱點
將532個主調控基因轉錄本映射到二倍體參考基因組上，泊松概率分布分析揭示了它們在特定基因組區域存在統計上顯著的聚集（熱點）。主要熱點位于5號染色體（10-15 Mb區域）和10號染色體（35-40 Mb的著絲粒周圍區域）。其他非隨機熱點分布在1-4、7和9號染色體上。這表明倍性變化可能影響包含主要調控基因的特定基因組區域，從而大規模改變表達網絡。
討論與意義
巴哈雀稗是南美天然牧場的主要草種，其天然的無孢子生殖特性已催生了新穎、快速的育種計劃。本研究通過比較性系二倍體與四倍體基因型的葉片轉錄組，提供了一個全面的葉片表達基因目錄，并鑒定了在基因組復制后上調或下調的基因。本體論分析揭示了與多倍體表型典型相關的類別富集，如光合作用、細胞周期和脅迫響應。轉錄因子和激素相關轉錄本的鑒定為靶向育種提供了關鍵候選目標。值得注意的是，主調控基因在基因組中并非隨機分布，而是聚集在特定熱點區域，這些區域可能是協調從2倍體到4倍體轉變過程中復雜生物途徑的基因組樞紐。
綜上所述，本研究的結果深化了對巴哈雀稗基因組復制轉錄組后果的理解，并突出了一系列可用于未來分子育種或基因組編輯計劃的候選基因和調控熱點，旨在改善該物種及其他重要農藝物種的脅迫耐受性、生產力和恢復力等理想性狀。結合已有的參考轉錄本組、參考基因組和成熟的轉化平臺，這些資源將助力設計靶向策略，推動分子輔助育種和遺傳工程的發展。
