《Scientia Horticulturae》:Integrative genome-wide association study and transcriptomic analyses unveil key candidate genes regulating fruit shape diversity in mango
編輯推薦:
為解析芒果果實形狀這一關鍵園藝性狀的遺傳機制,研究人員通過整合基因組組裝、全轉錄組RNA測序(RNA-Seq)、三維(3D)掃描技術以及全基因組關聯研究(GWAS),系統地研究了75份不同果形芒果種質。研究鑒定了與果形指數(FSI)等性狀顯著相關的99個遺傳位點和467個果形相關基因,并最終篩選出8個(如LOC123222869(TCP4-like)等)調控果形變異的關鍵候選基因。該研究為芒果果形改良育種提供了重要的基因資源與理論依據。
芒果作為廣受歡迎的熱帶水果,其果實的形狀千變萬化,從圓形、橢圓形、S形到腎形,不一而足。這份多樣化的“外貌”,不僅是吸引消費者的顏值擔當,更是影響其市場定價、包裝效率乃至運輸損耗的關鍵因素。然而,盡管果形如此重要,其背后的遺傳與分子調控機制卻像一個尚未打開的“黑匣子”。相比之下,西紅柿、黃瓜等作物果形研究的“基因地圖”已相對清晰,但芒果領域的研究相對匱乏。消費者偏愛外形美觀、形狀對稱的水果,市場也對果實的尺寸和均勻度有嚴格標準。因此,揭開芒果果形形成之謎,培育出滿足市場需求的優良品種,對于整個芒果產業的可持續發展和競爭力提升至關重要。
本研究正是為了破解這一難題而展開。研究者們集成了基因組學、轉錄組學、表型組學等多種技術手段,對芒果果形進行了系統性的“基因解碼”,相關成果發表在園藝學領域期刊《Scientia Horticulturae》上。
為系統探究芒果果形的遺傳基礎,研究人員主要采用了以下幾個關鍵技術方法:首先,他們利用PacBio長讀長與Illumina短讀長測序相結合的策略,組裝了圓形果品種‘金蘋芒’ (‘JPM’) 的高質量染色體級別基因組。其次,針對長果形品種‘紅象牙’ (‘HXY’) 與圓形果‘JPM’,在授粉后多個時間點進行采樣,通過高通量RNA測序(RNA-seq)獲取了全轉錄組表達數據,并利用加權基因共表達網絡分析(WGCNA)篩選與果形相關的基因模塊。第三,利用三維(3D)激光掃描技術對75份來自云南麗江芒果種質資源圃、具有不同果形的代表性種質資源進行了9個果形相關性狀(如果重、縱徑、橫徑、側徑、體積、果形指數等)的精確測量。最后,對這75份種質進行了全基因組重測序,獲得了海量的單核苷酸多態性(SNP)數據,并以此為遺傳標記,利用貝葉斯信息和連鎖不平衡迭代嵌套關鍵途徑(BLINK)模型,對關鍵果形性狀進行了全基因組關聯研究(GWAS),以定位相關的遺傳位點。
研究結果
3.1. ‘JPM’品種的基因組組裝與注釋
研究團隊成功組裝了圓形果芒果品種‘JPM’的高質量基因組。最終的染色體級別基因組大小為350.66 Mb,包含34,035個蛋白質編碼基因,重復序列占比39.28%。BUSCO評估顯示基因組完整度達97.70%。通過基因組比較與共線性分析發現,‘JPM’與另一圓形果品種‘Alphonso’基因組間存在強共線性,而與長果形品種‘HXY’的共線性較弱,存在較多倒位事件。這張圖展示了‘JPM’基因組的染色體特征及其與‘Alphonso’、‘HXY’基因組的共線性比較,直觀呈現了不同果形品種間的基因組差異。
3.2. 基于RNA-Seq表達模式的兩種芒果品系比較分析
通過對‘HXY’和‘JPM’果實發育過程(從授粉后1天到110天,共14個時間點)的形態觀察和轉錄組測序,研究人員發現授粉后第3天(DAP3)、第5天(DAP5)和第10天(DAP10)是果形指數(FSI,縱徑/橫徑)趨于穩定的關鍵發育階段。在這三個階段,‘HXY’與‘JPM’之間存在大量差異表達基因。結合WGCNA分析,鑒定出與果形指數最相關的藍色和黃色基因模塊。最終,通過整合差異表達分析與WGCNA,共篩選出467個最可能影響芒果果形發育的關鍵基因。GO富集分析顯示,這些基因主要富集在激素響應調節、細胞骨架依賴性胞質分裂等通路。該圖展示了兩種果形芒果在不同發育階段的果實形態、尺寸變化趨勢、差異基因韋恩圖以及WGCNA模塊與性狀的關聯熱圖,清晰地揭示了果形發育的關鍵時期和相關基因網絡。
3.3. 種質資源重測序揭示了海量SNP
對75份不同果形的芒果種質進行重測序,獲得了超過2600萬個SNP。基于高質量SNP構建的系統發育樹和主成分分析(PCA)將75份種質劃分為4個主要類群,且每個類群在果形上表現出明顯的偏好性,例如中國本土選育品種多呈象牙形,而引進品種更偏向圓形,這表明芒果種群進化與果形存在潛在關聯。這張圖呈現了SNP在染色體上的分布、類型構成以及基于SNP構建的種質系統發育樹,并標注了各分支上的代表性果形,直觀展示了遺傳多樣性與果形表型的對應關系。
3.4. 芒果果形多樣性分析
利用3D掃描技術獲取的9個果形相關性狀數據呈正態分布,且性狀間存在中等到強的相關性。基于這些性狀的PCA分析同樣將75份種質分為4個類群,分別代表了大而圓、小而圓以及兩種不同類型的長形果,證實了研究材料在果形上具有廣泛的代表性。該圖通過密度分布圖和PCA散點圖,展示了所有種質在果形性狀上的變異范圍和自然聚類情況。
3.5. 6個果形相關性狀的GWAS分析
對果重(FW)、果實體積(V)、縱徑(LD)、橫徑(CD)、側徑(SD)和果形指數(FSI)這6個性狀進行GWAS分析。以P值小于10-6為閾值,共鑒定出99個與性狀顯著關聯的SNP位點,其中與FSI相關的有75個,與V相關的有24個。在關聯區域(SNP上下游各100 kb)內共注釋到59個蛋白編碼基因。結合這些基因在‘HXY’和‘JPM’果實發育過程中的表達模式,最終篩選出44個關鍵候選基因。其中,包括LOC123222869(TCP4-like)、LOC123213947(bZIP61-like)、LOC123210819(ABC-G10)、LOC123200572(ARF3) 和LOC123228419(TTL1)在內的8個基因被提議為調控果形變異的關鍵候選者。該圖包含GWAS的曼哈頓圖、QQ圖以及關鍵候選基因的表達熱圖,是定位果形遺傳位點和篩選核心基因的直接證據。
結論與討論
本研究通過整合多組學分析,系統解析了芒果果實形狀的遺傳基礎。主要結論如下:第一,成功組裝了圓形果芒果代表品種‘JPM’的高質量基因組,填補了該果形類型基因組資源的空白。第二,明確了DAP3、DAP5和DAP10是芒果果形發育的關鍵時期,并通過轉錄組分析篩選出467個果形相關基因。第三,基于75份種質的重測序和3D表型數據,通過GWAS定位到99個與果形性狀顯著關聯的遺傳位點。第四,通過整合GWAS與轉錄組數據,最終鑒定出44個關鍵候選基因,并重點提出了8個(如TCP4-like、bZIP61-like等)作為調控果形變異的核心候選基因。
討論部分深入闡述了這些發現的意義。研究指出,芒果果形的多樣性與其在印度和東南亞的兩個獨立馴化中心有關,中國早期引進的品種多源自東南亞,其象牙形果狀與此吻合。本研究所用的75份中國栽培品種在遺傳和表型上均呈現出豐富的多樣性,并可劃分為四個明顯類群,這為針對不同市場需求進行定向育種提供了材料基礎。對于篩選出的關鍵候選基因,研究者結合其在其他物種中的同源基因功能進行了深入探討。例如,TCP4-like基因在擬南芥和番茄中已被證實是調控器官大小和形狀的關鍵因子,其在‘JPM’果實發育過程中持續高表達,可能通過抑制細胞增殖和伸長來促進圓形果形成。bZIP61-like作為一個轉錄因子,可能通過整合生長素和乙烯等激素信號,調控下游細胞壁重塑相關基因的表達來影響果形。ABC-G10屬于細胞分裂素轉運蛋白家族,可能通過調控激素的長距離運輸影響發育。TTL1則可能通過調節微管動力學來協調細胞擴張與增殖,從而影響最終果實大小和形狀。
這項研究的重大意義在于,它首次通過如此系統的多組學整合策略,大規模揭示了芒果果形性狀的遺傳架構和關鍵調控基因。這不僅深化了我們對芒果這一重要熱帶果樹發育生物學的理解,更重要的是,為未來開展基于基因組的精準育種(genome-based breeding)提供了直接的靶點和豐富的遺傳資源。育種者可以利用這些鑒定的遺傳位點和候選基因,通過分子標記輔助選擇(marker-assisted selection)或基因編輯等技術,更高效地培育出果形優美、均勻一致、更符合市場與加工需求的芒果新品種,從而提升產業效益和競爭力。盡管本研究采用的72份種質群體規模在一定程度上限制了對微效基因的檢測,但其所鑒定出的主效基因為后續功能驗證和育種應用奠定了堅實的基礎。未來,擴大種質群體、結合多環境數據和高密度標記將進一步細化果形的遺傳結構,并驗證這些關鍵基因位點的效應。