茶，作為世界上最古老的飲料之一，不僅有著令人愉悅的香味，更是承載著千年文化的象征。茶樹的種植面積和產量在全球范圍內十分龐大，經濟價值相當可觀，是茶葉出口國和多個發(fā)展中國家消除貧困和糧食安全的重要保障。茶樹的育種工作也一直備受關注，然而傳統(tǒng)茶樹一直依賴于對自然變異和有價值的基因型重組的長期選擇。如何更加精確、高效地育種一直是茶樹育種家關注的焦點。
為此，科學家們提出了利用泛基因組（Pan-genome）輔助育種技術，有望實現(xiàn)更快速地植物性狀改良。泛基因組包含更加全面的遺傳信息，可以有效降低參考基因組偏差對遺傳變異檢測的影響。通過該技術，人們能夠更全面地了解植物的遺傳特性。泛基因組育種方法已經在水稻、番茄等一些常見模式作物上取得了成功，并且正逐漸被應用于小麥、高粱、土豆等非模式作物的遺傳改良。
圖1|茶的多樣性
為了更好地理解茶樹的遺傳特性和表型多樣性，基因組所張興坦課題組主導，聯(lián)合國內多家單位成功構建茶樹首個泛基因組圖譜，完成了22個茶樹品種的測序并成功構建茶樹泛基因組，并在此基礎上揭示了茶樹基因組中近期基因大規(guī)模擴張的遺傳基礎，并鑒定了多個參與茶樹葉色差異，芽期，香氣等多個重要農藝性狀相關的遺傳變異，為茶樹重要基因挖掘和分子設計育種提供了重要參考。研究成果于2023年11月28日以“Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis”為題發(fā)表在國際期刊《自然植物（Nature Plants）》上。
研究選取了22個能最大化代表茶樹遺傳多樣性和表型多樣性的茶樹品種，包括大葉茶（C.sinensis var.assamica,CSA）,小葉茶（C.sinensis var.sinensis，CSS）和白毛茶（C.sinensis var.pubilimba(CSP)）進行測序組裝。分析發(fā)現(xiàn)，茶樹基因組中有高達5萬至6萬個蛋白編碼基因，相較于已發(fā)表的版本(1–6)，編碼基因數(shù)目增加了近20%。茶樹基因組包含大量的重復序列，其中長末端重復轉座子（LTRs）占據了基因組序列的50%以上。茶樹基因組一次近期的LTR爆發(fā)事件（LTRburst，約30-50萬年前）導致蛋白編碼基因的迅速擴張（圖2）。
圖2.22個泛基因組測序的代表性茶樹品種。a.736份茶樹重測序主成分分析(PCA)，揭示其遺傳多樣性。彩色點圖是用于泛基組分析的茶樹品種。b.22個代表性茶樹品種單拷貝基因（single-copygenes）的系統(tǒng)發(fā)育分析。C.sinensis var.assamica(CSA)樣品為紅色背景，C.sinensis var.pubilimba(CSP)樣品為黃色背景，C.sinensis var.sinensis(CSS)樣品為淺藍色背景。右圖顯示了CSA和CSS的形態(tài)學差異。
通過茶樹基因組的結構變異分析，發(fā)現(xiàn)茶樹基因組中的結構變異序列（SVs）大多源自于轉座元件（TEs），TE的高度活躍性產生了基因組中大量SV。此外，本研究揭示了SV如何影響茶葉的品質特征，如葉色的多樣性。通過對茶樹的轉錄組和代謝數(shù)據的分析，發(fā)現(xiàn)不同品種之間的葉色差異主要與葉綠素、類黃酮、花青素等相關代謝物的含量有關。通過基因的SV分析，發(fā)現(xiàn)了結構變異可能影響相應代謝物的生物合成，從而造成了葉色的差異。比如，MYB114是調控花青素合成的一個重要轉錄因子，其啟動子區(qū)域的一段LTR插入是形成茶樹紫葉的關鍵因素（圖3）。
圖3.影響茶樹葉色相關的功能性結構變異。f.CsMYB114的示意圖。藍色矩形表示基因外顯子，紅色三角形指啟動子的一段LTR插入。g.不同葉色品種CsMYB114啟動子區(qū)域LTR插入的Read深度圖。紫葉品種（‘ZJ’,‘JMZ’，‘JGY’）的代表品種以及綠葉茶品種（‘FDDB’,‘BHZ’，‘GH3H’）。
芽期是（TBF）影響茶樹的一個關鍵經濟性狀，培育早生優(yōu)質茶樹品種是一個重要的育種目標。本研究通過圖形泛基因組關聯(lián)分析（pan-GWAS）深入挖掘了與茶樹芽期相關的遺傳變異，鑒定出了與早芽萌發(fā)相關的QTL(qSPI4)(7)。同時，新發(fā)現(xiàn)一個潛在參與早芽性狀（EBF）基因CYP72A9（圖4）。該研究揭示了pan-GWAS在識別茶葉性狀相關遺產變異的有效性。
圖4.基于圖形泛基因組的GWAS在挖掘茶樹TBF性狀相關關鍵基因中的應用。a.曼哈頓圖顯示了TBF性狀的GWAS結果，包括SNP-GWAS和SV-GWAS。b.不同品種間DREB2A-like等位基因的序列比對。具有LBF性狀的個體表現(xiàn)為C基因型，而具有EBF性狀的樣本表現(xiàn)出A基因型。c.遺傳變異與染色體4上的EBF性狀之間的關聯(lián)分析。藍色矩形表示在染色體4上顯著標記-性狀關聯(lián)的候選基因區(qū)域（227.490–227.520Mb）。d.EBF和LBF品種的DREB2A-like等位基因的系統(tǒng)發(fā)育分析。e.一個61kb的結構變異（SV_2_83963），位于CsCYP72A9（TGY012913）的上游。f.群體分析支持CsCYP72A9上游的SV_2_83963與TBF性狀相關。
此外，結合圖形泛基因組、群體基因組學以及代謝組數(shù)據的綜合分析，本研究還鑒定到多個影響茶葉風味性狀的等位變異，包括兒茶素合成，茶香氣合成的核心基因的等位變異，為茶樹分子育種改良提供了重要靶標。
茶樹泛基因組圖譜的構建更新了大眾對茶樹基因組的認識，包括蛋白編碼基因數(shù)以及遺傳變異對茶樹表型的影響。該研究不僅為茶樹基因組和重要性狀的演化提供新見解，而且為茶樹重要基因挖掘和遺傳育種提供了重要平臺資源和信息，為加速茶樹育種進入全基因組設計時代奠定了重要基礎。同時也為木本植物的生物學研究提供重要參考。
基因組所博士研究生陳帥、西北農林科技大學王鵬杰教授、基因組所孔維龍副研究員以及客座碩士研究生柴琨為論文共同第一作者，基因組所張興坦研究員、福建農林大學葉乃興教授和廣東省農業(yè)科學院茶葉研究所吳華玲教授為論文共同通訊作者。此外，博士研究生張晟鋮、王毅斌、余嘉鑫、高瑩瑩、蔣夢薇、雷文龍、王文玲、陳曉、屈伸洋、王芳、王英豪、谷夢雅，以及基因組所博士后張清、廣東省農業(yè)科學院茶葉研究所方開星博士、中國農科院茶葉研究所馬春雷副研究員、福建農林大學孫威江教授也參與了此項工作。該研究得到了廣東省重點領域研發(fā)計劃子項目、深圳市優(yōu)秀青年科學基金、國家優(yōu)秀青年科學基金的資助。
來源：中國農科院基因組所，信息貴在分享，如涉及版權問題請聯(lián)系刪除
