組學技術在黑果枸杞資源研究中的應用

李明珠，劉增根

·綜 述·

組學技術在黑果枸杞資源研究中的應用

李明珠1, 3，劉增根1, 2*

1. 中國科學院西北高原生物研究所青海省藏藥研究重點實驗室，青海 西寧 810001 2. 湖北中醫藥大學藥學院，湖北 武漢 430065 3. 煙臺大學生命科學學院，山東 煙臺 264005

黑果枸杞為我國西北特殊生境地區食藥兼用的生態經濟植物，其漿果富含花色苷、花青素、多酚、多糖等物質，是醫藥保健品、食品添加劑（天然色素）的主要成分和原料。隨著測序技術和各組學學科的不斷發展，組學技術已被廣泛應用于黑果枸杞資源的各項研究中，尤其從單一的基因組、轉錄組、代謝組、蛋白質組到多組學聯合分析，對黑果枸杞品質評價、耐鹽堿、果實發育機制以及花青素合成與積累等研究更加深入。綜述了近年來組學技術在黑果枸杞資源綜合評價、抗逆、次生代謝物生物合成代謝調控中的應用，探討了資源品質差異成因和生態環境適應機制，為黑果枸杞資源的保護、種質創新及合理開發利用提供科學依據。

黑果枸杞；次生代謝物；花色苷；花青素；多組學；藥用資源

黑果枸杞Murray為茄科（Solanaceae）枸杞屬L.多棘刺灌木（圖1），其果實味甜多汁，民間常生食或榨汁做飲料，具有滋補強壯、明目及降壓作用，其根可以治咳嗽、哮喘、感冒和發燒，葉民間當茶泡飲[1]。作為我國西北鹽堿、荒漠地區具有重要生態和經濟價值的藥食兼用資源植物，黑果枸杞含有花色苷、花青素、多酚、黃酮、多糖、生物堿、氨基酸、礦質元素、維生素等多種成分，具有補腎益精、養肝明目、生津止渴、潤肺止咳、補血安神等功效[2-3]，尤以其主要次生代謝物花青素類化合物已被廣泛應用于醫藥保健品、食品添加劑、化妝品等領域，并在我國傳統民族醫藥（藏藥、蒙藥、維藥）中占重要地位[4]。黑果枸杞所處的獨特生態地理環境，如生長環境緯度高、鹽堿化程度高、紫外線強、晝夜溫差大、日照時間長等，都有助于次生代謝產物（如花青素、多酚、黃酮、多糖等）的轉化與積累[2,5-7]。近年來，資源過度開發、質量不穩定、產地來源混亂和有效成分地理變異顯著等問題已經明顯制約了黑果枸杞產業的發展[8-9]，因此亟待尋找黑果枸杞品質差異形成原因以及構建資源（道地性）生態環境評價體系。

圖1 黑果枸杞的花和果實

隨著高通量測序技術的出現、色譜/質譜技術的快速發展以及分子生物學新方法的開發[10-11]，促進了以基因組學研究為基礎，涵蓋轉錄組學、蛋白組學、代謝組學、相互作用組學和表型組學等系統生物學研究模式和多組學研究方法的集成（圖2），將為現代藥用植物資源綜合評價、次生代謝物合成調控、種質創新等研究領域帶來新的機遇[12-14]。因此，本文綜述了近年來組學技術在黑果枸杞資源評價、耐鹽堿、花青素合成與代謝、果實發育機制等研究中的應用，為探索黑果枸杞品質形成機制、解析資源生態環境適應機制、科學發展人工種植奠定基礎，為國家特色民族資源保護與利用提供參考。

1 黑果枸杞資源及其綜合評價研究

黑果枸杞是一種多年生灌木，具棘刺，多分枝，漿果球形，成熟后為紫黑色，有很高的光合效率和較強抗逆性，是我國西北荒漠、鹽堿地區一種特有的藥食兼用的民族資源植物；在我國主要分布于內蒙古西部、陜西北部、寧夏、甘肅、青海和新疆等地，中亞、高加索和歐洲等地區亦有分布[15-16]。黑果枸杞蒙名為“喬諾英-哈爾馬格”、藏藥名“旁瑪”，藏醫用它治療心熱病、心臟病、月經不調、停經等病癥。《維吾爾藥志》記載，黑果枸杞果實及根皮常用于治療尿道結石、癬疥、齒齦出血等病癥。黑果枸杞果實味甜多汁且富含花色苷、花青素、生物堿、多酚、黃酮、維生素、有機酸及糖類，民間常生食或榨汁做飲料，可起到滋補強壯、明目、降血壓及預防心腦血管疾病的功效。現代藥理學研究發現，黑果枸杞果實及葉提取物具有增強機體免疫、抗氧化、抗疲勞、調血脂、降尿酸、保肝護腎、抗動脈粥樣硬化等功效[1-2]。2017年2月6日，黑果枸杞被國家衛生健康委員會正式批準納入普通食品管理，從而進入更加廣闊的消費市場。近年來，黑果枸杞種植業發展迅猛，尤以青海、新疆等地的種植面積不斷擴大[17]。據青海省林草局統計，2018年青海省黑果枸杞種植面積達1.20萬hm2，干果產量約1.85萬t，產值達11.67億元；該產業帶動就業人數達3.25萬人，當地農牧民每月人均增收3100多元，成為當地群眾脫貧致富、“鄉村振興”的富民產業[17-19]。

圖2 多組學聯合分析在藥用植物資源研究中的應用

Fig.2 Application of multiomics analysis in medicinal plant resources

目前，黑果枸杞因獨特的生物活性和巨大的生態經濟價值使其成為國內外學者競相研究與開發的對象[1-2]，在保健品、食品、化妝品等行業需求量較大，資源的過度利用正嚴重威脅其野生種質資源。此外，隨著整個生態環境日趨惡化和人為過度采摘砍伐，野生資源正在大面積的減少，部分地區甚至出現枯竭[8]。因此，開展黑果枸杞居群遺傳多樣性和分子評價的研究有利于其野生種質資源的保護和可持續利用，并為黑果枸杞的科學育種及種質創新提供依據。劉增根[20]利用序列相關擴增多態性（sequence-related amplified polymorphism，SRAP）分子標記對我國青海、甘肅、寧夏和新疆的14個野生黑果枸杞居群的遺傳多樣性和遺傳結構進行了研究與分析。應用篩選出的31對SRAP引物組合對黑果枸杞基因組DNA進行聚合酶鏈式反應（polymerase chain reaction，PCR）擴增，共檢測到468個有效位點，其中多態性位點398個，多態位點百分率為85.04%。Nei基因多樣性和分子方差分析（molecular variation analysis of variance，AMOVA）顯示，總的遺傳變異中15.55%存在于居群間，84.45%存在于居群內，表明遺傳分化主要存在于種群內。居群間具有中等的遺傳分化（st＝0.215 5）；基因流（m）為1.819 9，表明各居群間有較高水平的基因流動，能有效防止基因空間異質性。阿力同·其米克等[21]采用分子標記技術（inter-simple sequence repeat，ISSR）對新疆南部黑果枸杞6個自然居群及甘肅2個自然居群共115個樣品進行了DNA多態性分析。從60個隨機引物中篩選出7個有效引物，共產生64條DNA片段，其中50條為多態性條帶，多態位點百分率為78.1%。黑果枸杞個體間非加權組平均法聚類結果表明，同一居群的個體不能完全聚在一起，來自新疆和甘肅2區域的黑果枸杞材料也不能完全分開。探討了可能造成上述居群遺傳結構模式的主要因素，同時提出了今后對新疆南部黑果枸杞保護工作中需進一步解決的問題。

本課題組基于前期調查與研究發現，同一產區種植和野生黑果枸杞的主要成分（次生代謝物和營養成分）種類一致，且含量差異不顯著[3,9]，提示發展人工種植黑果枸杞可能是一種更好并值得推廣的方式，在保障與野生黑果枸杞品質一致的情況下實現該資源的健康持續利用。然而，不同產區的黑果枸杞部分表型性狀及有效成分種類和含量差異明顯[3,9,22]（尤以花色苷、多糖、總多酚、總黃酮等次生代謝產物差異顯著），由于缺乏對其次生代謝物差異機制、質量標準與控制、生境適宜性評價等多方面的綜合研究，導致出現種質創新能力差、質量不穩定、市場秩序混亂、種源混雜等諸多問題[2-3,9,18-20]。因此，資源品質差異明顯和適宜性評價體系缺乏嚴重制約了黑果枸杞產業發展，如何保護黑果枸杞資源并在保護的基礎上合理開發利用這筆寶貴的財富，同時探索資源品質形成過程和生境適宜性評價與區劃研究，探究組學技術在資源研究中的應用，以及科學發展人工種植緩解野生資源危機，已是從事黑果枸杞研究和開發人員所關注的重大課題。

2 基因組學在黑果枸杞研究中的應用

基因組學是一門研究生物基因組以及如何利用基因的學科[23]，主要通過對個體與群體的基因進行定性定量分析，挖掘基因型與表型之間的關系，揭示生物機體可能發生什么。基因組學研究方法有基因組從頭測序、重測序和簡化基因組測序，該研究主要分功能基因組學和結構基因組學2個方面，后者的目標是通過基因組測序與拼接獲得全基因組序列。結構基因組學研究包括遺傳作圖、基因組de novo測序、葉綠體和線粒體基因組、比較基因組學[24]。目前有關黑果枸杞結構基因組學研究主要集中在葉綠體基因組學分析。

葉綠體是植物特有的細胞器，葉綠體基因組的特點是拷貝數高、進化速率適中、單親遺傳和不存在基因重組等。隨著現代測序技術和高通量、高分辨率技術平臺的不斷發展，葉綠體基因組分析在系統發育樹構建和物種親緣關系研究中得到廣泛應用[25-26]。劉晶星[27]基于Illumina HiSeq 2 000測序平臺對黑果枸杞葉片基因組進行分析，獲得葉綠體基因組總長154 965 bp，GC含量為37.91%，基因區所占的比例接近60%；在130個注釋基因中，蛋白基因86個，rRNA基因8個，tRNA基因36個。另外，基因編碼乙酰輔酶A羧化酶的CT亞基，是乙酰輔酶A羧化酶表達水平的限制因子；而黑果枸杞的基因有3個位置的突變，造成其編碼的蛋白只有480 aa。而Yisilam等[28]學者測得整個葉綠體基因組的長度為154 996 bp，由2個反向重復區域（25 395 bp）、1個85 993 bp的大單拷貝區和1個18 213 bp的小單拷貝區組成。基因組的總體GC含量為37.90%。葉綠體基因組由111個獨特基因組成，包括30個tRNA基因、4個rRNA基因和77個蛋白質編碼基因。黑果枸杞的序列與其他茄科葉綠體基因組一致，系統發育樹由MEGA6構建，以旋花科（Convolvulaceae）物種為外群，系統發育分析表明，黑果枸杞與茄科植物顛茄L.、馬尿泡Maxim.和(Dunn) Nakai聚在一起，具有密切的親緣關系。

Wang等[29]對寧夏植物園采集的黑果枸杞鮮葉進行完整葉綠體基因組測序，分析得到葉綠體基因組的大小為154 979 bp，包括2個反向重復區域（25 395 bp）、1個85 984 bp的大單拷貝區和1個18 205 bp的小單拷貝區組成。共預測了132個基因，包括37個tRNA、8個rRNA和86個蛋白質編碼基因。此外，9個PcG基因具有單個內含子，74個PcG基因沒有內含子，另外3個基因具有2個內含子，6個tRNA基因是單個內含子；葉綠體基因組的總GC含量為37.90%。對黑果枸杞和其他23個完整的茄科物種葉綠體基因組進行了系統發育分析，包括茄亞科和煙亞科以及旋花科的2個物種牽牛L、大星牽牛(Kunth) G. Don和另外2個物種(Lour.) Kuntze和Nakai為外群，用RAxML構建的最大似然樹，系統發育分析進一步驗證了黑果枸杞歸入茄科植物。

Cui等[30]也獲得了黑果枸杞葉綠體基因長度（15 486 9 bp），總GC含量為37.90%。葉綠體基因組共編碼113個差異基因，包括蛋白編碼基因85個、tRNA 37個、rRNA 8個和假基因3個，基因在黑果枸杞中發生突變，且包含18個具有內含子的基因。通過構建最大似然法系統進化樹發現，寧夏枸杞L和中華枸杞Mill.有較近的親緣關系，與黑果枸杞聚為姐妹分支。上述葉綠體基因組的測序工作，豐富了黑果枸杞的基因組學數據，為黑果枸杞資源表型多元化、生態適宜性和種質評價研究提供了科學依據。

3 轉錄組學在黑果枸杞研究中的應用

轉錄組學是從核酸水平研究特定細胞和組織在某一發育階段或功能狀態下轉錄出來的所有RNA的總和[31]，揭示生命機體正在發生什么。轉錄組的定義中包含時間和空間限定，通過測序分析可檢測所有正在表達基因的變化動態，并幾乎可獲得所有轉錄本，并闡明細胞或組織的生命活動狀態。高通量測序主要包括Roche 454測序技術、Illumina Sanger測序技術和ABI SOLID測序技術，目前應用最廣泛的測序技術是由Illumina公司開發的MiSeq和HiSeq平臺，該技術測序具有通量高、準確率高、速度快和成本低等特點[32-33]。

3.1 鹽脅迫轉錄組學分析

目前，全世界超過1/3的農業面積受到鹽堿化的影響。土壤鹽分是一個日益嚴重的全球性問題，它會影響植物生長和作物產量，給現代農業帶來嚴重問題[34]。滲透脅迫伴隨著鹽分脅迫會導致氣孔迅速關閉，從而降低植物吸收二氧化碳的能力[35]。鹽脅迫已經造成對枸杞屬物種光合作用和蒸騰作用的負面影響[36]。

Qin等[37]對黑果枸杞轉錄組學分析，在高鹽度處理下，有141個基因差異表達，其中80個基因上調，61個基因下調。黑果枸杞鹽敏感基因主要富集于內質網的碳代謝和蛋白質加工有關的信號轉導途徑。Chen等[38]利用對照和混合鹽堿處理黑果枸杞幼苗的mRNA構建混合cDNA文庫，在優化組裝后，共獲得68 063個獨特的轉錄序列，平均長度為877 bp；在這些序列中，4096個單基因在鹽堿混合處理后上調，4381個單基因在鹽堿混合處理后下調。

Mo等[39]開展了鹽脅迫處理黑果枸杞幼苗（21 d生，200 mol/L）和成年樹（3年生，300 mol/L）的生理表型分析，采用加權基因共表達網絡分析轉錄組數據挖掘與花青素合成、鹽脅迫響應相關的候選LrMYB112類轉錄因子，并對其調控花青素、鹽脅迫相關基因表達的功能進行了初步探究。研究結果發現，鹽脅迫促進黑果枸杞花青素生物合成，啟動耐鹽脅迫的生理反應機制，并證實了LrMYB78a、LrMYB78b和LrMYB112轉錄因子參與黑果枸杞花青素生物合成與響應鹽脅迫。

3.2 果實發育過程中轉錄組學分析

Ma等[40]通過研究與花青素積累相關的綠色成熟階段、表皮顏色變化階段和完全成熟階段的黑果枸杞轉錄組，結果發現，在43 573個組裝的單基因中，有12 734個在黑果枸杞的果實成熟過程中差異表達；鑒定出25個顯著差異表達的結構基因（、、、、、、、、、、）可能與花青素生物合成相關。

Li等[41]基于Illumina HiSeq測序平臺對黑果枸杞5個不同發育階段果實顏色進行轉錄組分析，結果顯示，在轉錄組中共鑒定出31 302個差異基因，72個關鍵候選基因對參與類黃酮和花青素生物合成的14種酶產生直接影響。朱雪冰等[42]在黑果枸杞中測得1.167×108和9.695×107reads，篩選出了13個與花青素生物合成有關的基因，克隆黑果枸杞中基因的轉錄因子Sl AN2，該基因全長774 bp，編碼257個氨基酸，相對分子質量為29 775.84，等電點為7.79。嚴莉等[43]基于轉錄組數據分析基因在黑果枸杞果實不同發育期的差異表達模式，結果表明，基因可能參與了果實不同發育時期花青素變化的調控，利用熒光定量PCR的差異表達數據進一步驗證了部分MYB轉錄因子在果實不同發育時期的花青素合成中可能起到調控作用；同時在黑果枸杞基因家族注釋得到83個MYB類轉錄因子基因，為進一步研究MYB家族的基因結構和生物學功能奠定了基礎。

彭勇等[44]對黑果枸杞3個不同發育時期的果實進行轉錄組測序分析表明，京都基因與基因組百科全書（Kyoto encyclopedia of genes and genomes，KEGG）信號通路富集顯示亮氨酸生物合成途徑在變色期對比青果期有7個差異表達基因，黑熟期對比變色期有35個差異表達基因，且全部為上調表達；在糖酵解途徑中3個不可逆的反應步驟限速酶基因在黑熟期均為上調表達。孟小偉等[45]對黑果枸杞果實不同發育過程進行轉錄組測序，共得到43 573條Unigene，有23 723條Unigene在基因本體（gene ontology，GO）、KEGG、真核同源群簇（euKaryotic orthologous groups，KOG）、非冗余蛋白（nonredundant protein，NR）和SwissProt 5個數據庫中得以注釋，篩選出與類黃酮生物合成有關的Unigene 45個，這為與黑果枸杞果實品質有關基因篩選、克隆和功能分析等提供研究基礎和理論依據。

3.3 花青素生物合成轉錄組學分析

包雪梅等[46]基于轉錄組測序和生物信息學分析，從黑果枸杞和白果枸杞（黑果變異）中共檢測到25 279個基因差異表達，其中12 381個基因在黑果枸杞表達上調，12 898個基因表達下調；KEGG富集結果顯示，共有571個差異表達基因參與植物次生代謝產物合成，尤其是在黃酮類和花青素生物合成代謝通路中差異最為明顯。黑果枸杞中共有8個參與花青素合成代謝的結構基因的表達量均高于白果枸杞，表明花青素的生物合成可能與黑果枸杞黑果性狀相關。調控花青素生物合成的bHLH轉錄因子CL8159.Contig5_All的表達水平比白果枸杞增強了2 139.57倍，推測這個基因很有可能是控制黑果枸杞黑果性狀的候選基因。

Zeng等[47]采用Illumina Hiseq 2000平臺測定了黑果枸杞果實在5個不同發育階段花青素生物合成相關的基因序列，分離出參與黑果枸杞花青素生物合成的結構基因（、、、、、和），從轉錄組數據庫中檢索到編碼、、、、、、、和的序列，并應用Mega 4.0對相關基因進行系統發育分析。上述基因在果實發育過程中的表達譜表明，在果實中轉錄豐度與花青素積累之間存在顯著的正相關；和或可能共同調控、、、和的轉錄，從而調控黑果果實中花色苷的合成。此外，調控基因和結構基因的表達模式以及分支節結構基因/的轉錄比率可能決定了黑果枸杞和寧夏枸杞果實之間花青素生物合成的表型差異。

Zong等[48]研究等位基因存在與果實顏色的關系，根據轉錄組分析結果從黑果枸杞和寧夏枸杞中分離得到編碼MYB轉錄因子的和。在系統發育樹中，和與調節煙草花青素合成的關系密切，而的過度表達誘導煙草所有組織中花青素的生物合成，且明顯強于。此外，轉錄因子僅在黑果果實中檢測到，并在果實發育過程中增加，伴隨著花色苷的積累；而的功能多樣性和高表達水平可能是黑果枸杞果實花青素含量高的原因。上述研究結果為黑果枸杞資源改良與利用提供了重要基因資源，為花青素的生物合成與代謝提供了堅實的理論基礎。

4 蛋白質組學在黑果枸杞研究中的應用

蛋白組學是以蛋白質組為研究對象，研究細胞、組織或生物體蛋白質組成及其變化規律的科學[49]，它主要利用高效液相色譜-串聯質譜法技術獲取某一時刻生物機體內所有表達蛋白，揭示生物體已經發生什么。蛋白組學是基因表達調控機制研究的一個重要工具。當前，隨著高精確度、高分辨率、高靈敏度、高通量的生物質譜儀器的快速發展，以及各物種基因數據庫、蛋白質數據庫的不斷完善，配合基于復雜算法的生物信息學工具，蛋白質組學的研究內容、研究方法和研究領域都發生了翻天覆地的變化[50]。目前的蛋白質組學研究方法有雙向凝膠電泳技術、質譜技術、同位素的相對和絕對定量標記技術、蛋白質芯片技術、酵母雙雜交系統技術和免疫共沉淀等技術；其中，同位素標記相對和絕對定量（isobaric tags for relative and absolute quantification，iTRAQ）、串聯質譜標簽（tandem mass tags，TMT）是最常用的基于同位素標記質譜分析的定量蛋白質組學分析技術，可對多達16（8）個不同樣本進行蛋白定量比較，標記效率和質譜檢測靈敏度高，適用范圍廣泛[51]。

同一生物體在不同生長時間和生長環境下，基因表達情況并不完全相同，進行多器官、多時期、多維度、多水平的藥用植物蛋白質組學研究，建立藥用植物功能基因庫，可為藥用植物種質資源評價、中藥材質量控制和標準建立、新品種選育、有效成分的合成生物學等領域提供分子基礎，推動藥用植物產業的健康可持續發展[52]。

目前，蛋白質組學技術在植物的研究領域使用廣泛，然而關于黑果枸杞蛋白質組學的研究較少。汪亞娟等[53]采用iTRAQ技術鑒定出862個黑果枸杞蛋白質，經100 mmol/L鹽處理7 d后黑果枸杞檢測出165個表現差異的蛋白質，其中有72個蛋白質上調，93個蛋白質下調。并了解到鹽脅迫主要影響了黑果枸杞的能量代謝與離子運輸，其中參與離子運輸的差異蛋白質共4個，全都是上調的膜孔蛋白；參與能量代謝的差異蛋白質共有13個，包括9個上調蛋白和4個下調蛋白。此外，黑果枸杞受到鹽脅迫時有利于黃酮類化合物的大量合成，并且以增加對香豆酸進入木質素生物合成途徑的方式來促進黃酮類化合物的積累，從而緩解鹽脅迫帶來的氧化損傷。因此，體內黃酮類化合物的積累可能是鹽生植物黑果枸杞對抗高鹽環境的機制之一，而且該研究鑒定出的一些有價值的蛋白質，為進一步研究鹽生植物耐鹽機制和為作物耐鹽馴化提供了新的依據。

5 代謝組學在黑果枸杞研究中的應用

代謝組學是對內源性及外源性小分子代謝產物進行定性和定量檢測，并尋找與機體生命活動、病理變化相對應的代謝物質的研究方法[54]。代謝組學揭示了生命體活動事件確實發生了，它根據研究目的不同分為靶向代謝組學和非靶向代謝組學，兩者各有各的優勢和特點。靶向代謝組定性定量準確，但對物質的覆蓋率有限；而非靶向代謝組對物質的覆蓋率廣泛，但缺乏絕對定性定量的數據[55]。

代謝組學是對基因組學和蛋白質組學的補充，常見的分析手段包括液質譜聯用（liquid chromatograph mass spectrometer，LC-MS）、氣質譜聯用（gas chromatograph mass spectrometer，GC-MS）、超高效液質聯用（ultra performance liquid chromatography mass spectrometer，UPLC-MS）、毛細管電泳質譜聯用（capillary electrophoresis mass spectrometer，CE-MS）、核磁共振（nuclear magnetic resonance，NMR）等，四極桿線性離子阱（quadrupole-trap，Q-Trap）、四極桿飛行時間（quadrupole-time of flight，Q-TOF）及離子阱飛行時間（ion trap-time of flight，IT-TOF）等串聯質譜也被廣泛應用。植物代謝物種類繁多（已超過20萬種），很多醫藥資源、食品及工業原料都是植物代謝產物，因此，代謝組學在植物研究領域的應用尤為重要。通過植物代謝組學的定性和定量分析，可闡明不同生態地理環境下代謝物的變異規律、研究同一植物不同時期或不同部位或不同產區代謝物種類及含量的差異、推測相應的合成與代謝途徑以及代謝調控網絡等，為藥用植物資源的品質差異成因解析、物種鑒別以及活性成分合成分子機理研究提供理論依據[56]。目前，有關黑果枸杞代謝組學的研究主要集中在鹽脅迫和果實發育過程方面。

5.1 鹽脅迫代謝組學分析

Qin等[37]發現黑果枸杞在鹽脅迫下有69種代謝產物的表達水平發生了變化，其中34種表達上調，35種表達下調。研究還發現正常條件下黑果枸杞中的脫落酸含量顯著低于枸杞，而在鹽脅迫處理下黑果枸杞中的脫落酸含量急劇上調，但是枸杞中并未發生顯著變化。鹽脅迫在枸杞葉片中誘導的代謝物變化主要通過以下途徑富集：嘌呤代謝、半胱氨酸和蛋氨酸代謝以及蛋白質的消化和吸收；而在鹽度壓力下，黑果枸杞葉代謝產物差異主要集中于不同環境中的微生物代謝、嘌呤代謝、色氨酸代謝、苯丙烷類生物合成以及黃酮和黃酮醇生物合成代謝等途徑。此外，在正常條件下，黑果枸杞比枸杞含有更高的黃酮類化合物含量，提示黑果枸杞比枸杞更耐鹽。綜上，該研究表明脫落酸和黃酮類化合物在枸杞屬植物耐鹽中發揮重要作用，該研究結果對未來作物抗鹽育種和資源馴化有重要參考價值。

5.2 成熟漿果的代謝組學分析

Li等[41]采用LC-ESI-Q-Trap-MS/MS技術對黑果枸杞和白果枸杞黃酮代謝產物進行差異分析，研究鑒定出15種花青素，包括天竺葵素、矮牽牛素、錦葵素和飛燕草素等，在黑果枸杞果實中發現23種黃酮代謝物；2種漿果的代謝物在早期苯丙烷途徑上沒有差異，從黃酮類化合物合成代謝途徑開始差異顯著。黑果枸杞成熟果實主要積累飛燕草素衍生物，白果枸杞成熟果實中的黃酮醇含量是黑色果實的2～3倍；另外，黑果枸杞成熟果實具有更多的上游底物柚皮素查耳酮和柚皮素來產生下游黃酮化合物，且有更多的二氫楊梅素和二氫槲皮素來合成下游花青素。同時還發現了黑果枸杞漿果成熟過程中存在3種花青素轉運機制，即谷胱甘肽-轉移酶介導的轉運、多耐藥相關蛋白和多藥和有毒化合物外排轉運蛋白介導的轉運以及膜囊泡介導的轉運。

李麗等[57]采用UPLC-Orbitrap MS技術對植物乳桿菌GH-6發酵前后黑果枸杞漿的成分進行分析，共鑒定出61個代謝物。隨后，在發酵液、原漿以及添加蔗糖及乳清蛋白后的發酵液中篩選出22個差異性代謝物，其中2種糖類、3種有機酸類、3種嘌呤類、3種氨基酸類、4種花青素類和7種小肽類，這也說明乳桿菌GH-6發酵可顯著影響黑果枸杞漿中的代謝變化。Liu等[3,9]基于黑果枸杞植物化學成分（如飛燕草素、錦葵素、矮牽牛素及其花色苷類以及多糖、多酚、黃酮、生物堿等）的色譜和質譜研究以及代謝組學分析，建立了有效成分測定與分析技術平臺；同時分析了不同產區及野生和種植的黑果枸杞花青素、花色苷、黃酮、酚酸、多糖等成分含量差異以及與生態環境因子的相關性，構建了黑果枸杞化學指紋圖譜與有效成分定量分析相結合的化學多樣性評價體系。

彭玉嬌等[58]采用LC-MS非靶向代謝組技術對黑果枸杞和寧夏枸杞漿果進行代謝物差異分析，共檢測出6982種代謝物并獲得501種差異代謝物，黑果枸杞中表達上調的差異代謝物有240種，表達下調的代謝物有261種，前50的差異代謝物中有6種是氨基酸和多肽類物質、7種是黃酮類化合物。對總黃酮的含量分析發現，黑果枸杞總黃酮含量較高，有5種差異代謝物參與到黃酮合成代謝途徑中，另有15種差異代謝物參與4個氨基酸代謝相關通路，提示這些代謝物的差異表達可能導致寧夏枸杞與黑果枸杞氨基酸的差異。該研究為深入了解兩枸杞物種中的氨基酸代謝以及對它們的差異利用提供科學基礎。

6 多組學聯合分析在黑果枸杞研究中的應用

多組學聯合分析方法將基因、mRNA、調控因子、蛋白、代謝等不同層面之間信息進行整合，構建基因調控網絡，深層次理解各個分子之間的調控及因果關系，從而更深入探索和認知機體生長發育進程和抵御逆境過程中復雜性狀的分子機制和遺傳基礎[12-13,59]。因此，通過轉錄組學和代謝組學聯合分析，有助于解析黑果枸杞中代謝物合成的表達調控網絡及代謝網絡，為挖掘到與類黃酮（花青素）生物合成和抗逆（耐鹽堿）相關的關鍵酶基因提供支撐。通過轉錄組學和代謝組學的聯合分析顯示，正常條件下黑果枸杞中黃酮和類黃酮含量較高，相比于黑果枸杞鹽脅迫顯著提高了枸杞中的黃酮和類黃酮含量；另外，結果還表明枸杞中的鹽響應基因主要富集在絲裂原活化蛋白激酶信號、氨基糖和核苷酸糖代謝、碳代謝和植物激素信號轉導通路，而黑果枸杞中的鹽響應基因主要與內質網中的碳代謝和蛋白質加工有關[37]。基于轉錄組和代謝組數據聯合分析，證實花青素類化合物在黑果枸杞成熟漿果中的含量要明顯高于其他枸杞物種中的含量，而花青素生物合成相關基因在黑果枸杞中的表達量也要顯著性高于其他枸杞物種[41,47]。

此外，基于枸杞屬植物中的轉錄組和代謝組數據，通過組裝、拼接、比對等分析發現在黑果枸杞中存在一個調控花色苷合成代謝的MYB轉錄因子，其表達量隨果實的發育而增加，類黃酮物質含量也隨之增加，而該基因在寧夏枸杞中幾乎不表達[9,48]。自然群體中的等位變異與黑色果實性狀緊密關聯，在煙草里過量表達，除根部外，其他組織部位均表現出紫黑色。通過多組學聯合分析，初步發現轉基因煙草中185個基因（146上調、39個下調）產生差異表達，其中類黃酮生物合成代謝途徑差異最顯著。

代謝組學可以批量獲得基因表達的最終產物，反映了基因和蛋白表達的微小變化以及機體系統的生理和病理狀態；蛋白質組學能夠動態地描述基因調節，可獲取某一時刻機體內所有表達蛋白。因此，將蛋白質組學和代謝組學進行關聯分析，對基因表達的蛋白質和代謝物進行定量測定和差異比較，有助于進一步揭示基因表達調控網絡和次生代謝物的生物合成代謝機制。在黑果枸杞耐鹽機制研究中，學者利用蛋白質組學和代謝組學聯合分析方法，分析了黑果枸杞愈傷組織在不同濃度鹽脅迫下的蛋白質和化學成分，闡釋了差異性蛋白質基因表達主要富集于離子運輸和能量代謝途徑，并且鹽脅迫有助于植株黃酮類化合物的合成與積累，緩解逆境帶來的氧化損傷[53]。

盡管目前組學技術已經在黑果枸杞中得以廣泛應用，尤其集中在抗逆（耐鹽）和果實發育過程中的花青素合成代謝（表1），但是研究學者們對其進行的多組學聯合研究還為之甚少。顯然，單憑一種組學技術已經無法滿足黑果枸杞適應與進化、花色苷差異成因解析等相關深層次的研究，因此，基于多組學聯合分析闡明黑果枸杞次生代謝物（花青素、黃酮、生物堿等）合成調控機制、揭示黑果枸杞耐旱和抗鹽堿機制將是今后枸杞屬植物資源的研究重點。

表1 組學技術在黑果枸杞資源研究中的應用

a-葉綠體基因組 b-逆境脅迫（鹽脅迫） c-果實發育涉及的花青素生物合成

a-chloroplast genome b-stress (salt stress) c-anthocyanin biosynthesis in fruit development

7 結語及展望

隨著多組學和相應分子生物學技術的發展與滲透，藥用植物資源品質形成的探索不論在生態環境等外在因素，還是基因、蛋白質與遺傳分化的分子水平都有了一定的推進[12,60-61]。基于多組學聯合分析是藥用植物次生代謝產物合成途徑和資源品質差異成因研究的重要技術[62-63]。在進入大數據時代的背景下，基于系統生物學研究模式和多組學研究思路，能夠為藥用植物次生代謝物生物合成代謝調控、資源品質綜合評價以及種質創新等相關科學問題的解決打下堅實的理論基礎。

從已有的研究成果可以看出，各組學應用在黑果枸杞資源研究中已經取得了一定進展，為資源的品質差異、次生代謝產物合成途徑、有效成分確定以及生境適宜性等研究提供了重要參考。然而，有關黑果枸杞資源評價、組學研究、品質差異及其與生態環境關系方面的研究仍然存在一些不足。主要表現在：（1）很少有研究關注不同產區黑果枸杞次生代謝產物差異機制及生態環境對有效成分合成積累的影響，資源生境適宜性系統評價和區劃研究基本處于空白。（2）前期有關黑果枸杞組學研究主要集中在葉綠體基因組學和轉錄組學，缺乏蛋白質組和代謝組及其組學間的整合分析，組學間數據交集的挖掘深度不夠，且未實現多組學的真正融合。

綜上所述，隨著后續黑果枸杞全基因組數據的公布，并基于化學生物學分析和多組學聯合應用，鑒定黑果枸杞次生代謝物生物合成過程中的關鍵酶和調控因子，深入解析主要有效成分花色苷合成調控機制（圖3），結合生態環境因子的相關性分析，將對探究黑果枸杞花色苷差異成因、資源品質形成過程和生境適宜性區劃評價具有重要的理論指導作用，從而促進黑果枸杞資源健康可持續開發與利用。

PAL-苯丙氨酸裂解酶 C4H-肉桂酸-4-羥化酶 4CL-4-香豆酰CoA連接酶 CHS-查耳酮合成酶 CHI-查耳酮異構酶 F3H-黃烷酮-3-羥化酶 F3′H-類黃酮-3′-羥化酶 F3′5′H-類黃酮-3′,5′-羥化酶 DFR-二氫黃酮醇還原酶 ANS-花青素合成酶 3GT-3-糖基轉移酶 UFGT-類黃酮糖基轉移酶 MT-甲基轉移酶 AT-酰基轉移酶

利益沖突 所有作者均聲明不存在利益沖突

[1] 馬麗娟, 霍鵬超, 孫夢茹, 等. 黑果枸杞化學成分和藥理活性的研究進展 [J]. 中草藥, 2020, 51(22): 5884-5893.

[2] Wang H Q, Li J N, Tao W W,.studies: Molecular biology, phytochemistry and pharmacology [J]., 2018, 240: 759-766.

[3] Liu Z G, Liu B L, Wen H X,. Phytochemical profiles, nutritional constituents and antioxidant activity of black wolfberry (Murr.) [J]., 2020, 154: 112692.

[4] 張弓, 陳莎莎, 周武, 等. 黑果枸杞的功用考證及研究進展 [J]. 華西藥學雜志, 2019, 34(6): 638-642.

[5] Miller J S. Phylogenetic relationships and the evolution of gender dimorphism in(Solanaceae) [J]., 2002, 27(2): 416-428.

[6] Guo Y Y, Yu H Y, Kong D S,. Effects of drought stress on growth and chlorophyll fluorescence ofMurr. seedlings [J]., 2016, 54(4): 524-531.

[7] 袁海靜, 安巍, 李立會, 等. 中國枸杞種質資源主要形態學性狀調查與聚類分析 [J]. 植物遺傳資源學報, 2013, 14(4): 627-633.

[8] Liu Z G, Shu Q Y, Wang L,. Genetic diversity of the endangered and medically importantMurr. revealed by sequence-related amplified polymorphism (SRAP) markers [J]., 2012, 45: 86-97.

[9] Liu Z G, Dong B, Liu C,. Variation of anthocyanin content in fruits of wild and cultivated[J]., 2020, 146: 112208.

[10] Mosher J J, Bowman B, Bernberg E L,. Improved performance of the PacBio SMRT technology for 16S rDNA sequencing [J]., 2014, 104: 59-60.

[11] Zhang F, Wen Y, Guo X. CRISPR/Cas9 for genome editing: Progress, implications and challenges [J]., 2014, 23(R1): R40-R46.

[12] 黃璐琦, 高偉, 周潔, 等. 系統生物學方法在藥用植物次生代謝產物研究中的應用 [J]. 中國中藥雜志, 2010, 35(1): 8-12.

[13] 劉義飛, 胡志剛, 黃必勝, 等. 組學技術在中藥種質資源遺傳評價與創新中的應用 [J]. 世界科學技術—中醫藥現代化, 2022, 24(4): 1315-1324.

[14] 豐美靜, 張愷愷, 陳段芬, 等. 多組學技術在紅豆杉研究中的應用 [J]. 分子植物育種, 2022, 20(3): 840-846.

[15] 甘青梅, 左振常, 昌也平, 楊文蓮. 藏藥“旁瑪”的考證及生藥學研究 [J]. 中國民族民間醫藥雜志, 1995, 4(1): 31-33.

[16] 王曉宇, 陳鴻平, 銀玲, 等. 中國枸杞屬植物資源概述 [J]. 中藥與臨床, 2011, 2(5): 1-3.

[17] 劉增根, 康海林, 岳會蘭, 等. 黑果枸杞資源調查及其原花青素含量差異分析 [J]. 時珍國醫國藥, 2018, 29(7): 1713-1716.

[18] 楊麗萍. 提檔升級做強黑枸杞產業 [J]. 民生周刊, 2019(5): 78-79.

[19] 張佳琪, 趙英力, 王浩寧, 等. 黑果枸杞研究進展及產業發展對策 [J]. 北方園藝, 2020(19): 134-140.

[20] 劉增根. 黑果枸杞遺傳多樣性及可持續利用研究 [D]. 北京: 中國科學院大學, 2014.

[21] 阿力同·其米克, 王青鋒, 楊春鋒, 等. 新疆產藥用植物黑果枸杞遺傳多樣性的ISSR分析 [J]. 植物科學學報, 2013, 31(5): 517-524.

[22] Wang Z C, Yan Y Z, Nisar T,. Comparison and multivariate statistical analysis of anthocyanin composition inMurray from different regions to trace geographical origins: The case of China [J]., 2018, 246: 233-241.

[23] Marian A J. Sequencing your genome: What does it mean? [J]., 2014, 10(1): 3-6.

[24] 張聞婷, 焦萌, 王繼華. 藥用植物基因組學研究進展 [J]. 廣東農業科學, 2021, 48(12): 138-150.

[25] Clarke J. Progress of chloroplast genome analysis [J]., 2008, 35(1): 1-8.

[26] 樊守金, 郭秀秀. 植物葉綠體基因組研究及應用進展 [J]. 山東師范大學學報: 自然科學版, 2022, 37(1): 22-31.

[27] 劉晶星. 寧夏枸杞和黑果枸杞葉綠體全基因組測序及基因注釋分析 [D]. 北京: 中國科學院北京基因組研究所, 2013.

[28] Yisilam G, Mamut R, Li J,. Characterization of the complete chloroplast genome of(Solanaceae) [J]., 2018, 3(1): 361-362.

[29] Wang C P, Jia G L, Zhu L Z,. The complete chloroplast genome sequence of(Solanaceae), a traditional medicinal plant in China [J]., 2019, 4(2): 2495-2496.

[30] Cui Y X, Zhou J G, Chen X L,. Complete chloroplast genome and comparative analysis of three(Solanaceae) species with medicinal and edible properties [J]., 2019, 17: 100464.

[31] McGettigan P A. Transcriptomics in the RNA-seq era [J]., 2013, 17(1): 4-11.

[32] 陳惠, 辛麗麗, 龔婕寧. 基于全轉錄組測序技術的轉錄組學在中醫藥領域的應用前景分析 [J]. 環球中醫藥, 2013, 6(10): 759-763.

[33] 劉歡, 黃麗娜, 張奮強, 等. RNA-Seq高通量測序技術在果樹功能基因組學研究的應用進展 [J]. 生物技術進展, 2017, 7(2): 144-148.

[34] Earthscan. The state of the world’s land and water resources for food and agriculture: Managing systems at risk [J]., 2012, 66(4): 418-419.

[35] Mansour M M F, Ali E F. Evaluation of proline functions in saline conditions [J]., 2017, 140: 52-68.

[36] Zhang Z Z, He K N, Zhang T,. Physiological responses of Goji berry (L.) to saline-alkaline soil from Qinghai region, China [J]., 2019, 9(1): 12057.

[37] Qin X Y, Yin Y, Zhao J H,. Metabolomic and transcriptomic analysis ofchinese andunder salinity stress [J]., 2022, 22(1): 8.

[38] Chen J H, Zhang D Z, Zhang C,. Physiological characterization, transcriptomic profiling, and microsatellite marker mining of[J]., 2017, 18(11): 1002-1021.

[39] Mo S Q, A B, Wang Z Q,. Spatio transcriptome uncover novel insight into theseedling tolerant to salt stress [J]., 2022, 177: 114502.

[40] Ma Y J, Duan H R, Zhang F,. Transcriptomic analysis ofMurr. during fruit ripening provides insight into structural and regulatory genes in the anthocyanin biosynthetic pathway [J]., 2018, 13(12): e0208627.

[41] Li T T, Fan Y F, Qin H,. Transcriptome and flavonoids metabolomic analysis identifies regulatory networks and hub genes in black and white fruits ofMurray [J]., 2020, 11: 1256.

[42] 朱雪冰. 黑枸杞花青素合成代謝相關基因克隆及再生體系建立 [D]. 西寧: 青海師范大學, 2017.

[43] 嚴莉, 王翠平, 陳建偉, 等. 基于轉錄組信息的黑果枸杞MYB轉錄因子家族分析 [J]. 中國農業科學, 2017, 50(20): 3991-4002.

[44] 彭勇, 陳尚武, 馬會勤. 黑果枸杞果實成熟發育過程表達譜差異分析 [J]. 生物技術通報, 2016, 32(11): 144-151.

[45] 孟小偉, 牛赟, 馬彥軍. 黑果枸杞果實發育過程中轉錄組測序分析 [J]. 中南林業科技大學學報, 2020, 40(9): 147-155.

[46] 包雪梅, 胡娜, 宗淵, 等. 轉錄組學挖掘黑果枸杞中花青素生物合成代謝關鍵bHLH調控基因 [J]. 分子植物育種, 2021, 7: 1-20.

[47] Zeng S H, Wu M, Zou C Y,. Comparative analysis of anthocyanin biosynthesis during fruit development in twospecies [J]., 2014, 150(4): 505-516.

[48] Zong Y, Zhu X B, Liu Z G,. Functionaltranscription factor encoding gene AN2 is associated with anthocyanin biosynthesis inMurray [J]., 2019, 19(1): 169.

[49] Golaz O, Wilkins M R, Sanchez J C,. Identification of proteins by their amino acid composition: An evaluation of the method [J]., 1996, 17(3): 573-579.

[50] Larance M, Lamond A I. Multidimensional proteomics for cell biology [J]., 2015, 16(5): 269-280.

[51] 江小羊, 翟曉巧. 植物蛋白質組學技術的研究進展 [J]. 河南林業科技, 2021, 41(2): 17-20.

[52] 付強, 趙杰宏, 劉育辰, 等. 藥用植物蛋白質組學的研究進展 [J]. 貴州農業科學, 2018, 46(10): 16-19.

[53] 汪亞娟. 基于iTRAQ技術的鹽 (NaCl) 脅迫下黑果枸杞蛋白質組學研究 [D]. 蘭州: 蘭州大學, 2013.

[54] 許國旺, 路鑫, 楊勝利. 代謝組學研究進展 [J]. 中國醫學科學院學報, 2007, 29(6): 701-711.

[55] 范仕成, 高悅, 張慧貞, 等. 非靶向和靶向代謝組學在藥物靶點發現中的應用 [J]. 藥學進展, 2017, 41(4): 263-269.

[56] 王斌, 張騰霄, 趙倩, 等. 植物代謝組學在藥用植物中的應用進展 [J]. 中華中醫藥學刊, 2021, 39(5): 28-31.

[57] 李麗, 馮華峰, 周淳, 等. 植物乳桿菌發酵黑果枸杞的代謝組學研究 [J]. 化學試劑, 2022, 44(8): 1088-1096.

[58] 彭玉嬌, 崔學宇, 邵元元, 等. 利用LC-MS技術解析黑果枸杞、紅果枸杞代謝物的差異 [J]. 食品與生物技術學報, 2021, 40(9): 56-63.

[59] 劉鋒. 不同種重樓間的差異蛋白質組學和代謝組學研究 [D]. 北京: 軍事科學院, 2019.

[60] Saito K, Matsuda F. Metabolomics for functional genomics, systems biology, and biotechnology [J]., 2010, 61: 463-489.

[61] 王程成, 趙慧, 嚴穎, 等. 道地藥材品質形成機制的組學研究思路 [J]. 中國中藥雜志, 2018, 43(11): 2407-2412.

[62] Zhu G T, Wang S C, Huang Z J,. Rewiring of the fruit metabolome in tomato breeding [J]., 2018, 172(1/2): 249-261.e12.

[63] Wang Z, Cui Y, Vainstein A,. Regulation of fig (L.) fruit color: Metabolomic and transcriptomic analyses of the flavonoid biosynthetic pathway [J]., 2017, 8: 1990.

Application of omics technology inresources

LI Ming-zhu1,3, LIU Zeng-gen1, 2

1. Key Laboratory of Tibetan Medicine Research of Qinghai Province, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810001, China 2. School of Pharmacy, Hubei University of Chinese Medicine, Wuhan 430065, China 3. School of Life Science, Yantai University, Yantai 264005, China

is an edible and medicinal fruit tree, which widely distributes in the arid desert and saline-alkali region of Northwest China.berries have rich secondary metabolites (such as anthocyanins, anthocyanidins, polyphenols, and polysaccharides) and high economic and ecological values, which were widely used in the field of nutritional foods, biochemical drugs, and food additives. With the continuous development of sequencing technology and omics, omics technology (genomics, transcriptomics, proteomics, and metabolomics) and integrative multi-omics approach have become important means for the research ofresources (especially in the areas of quality evaluation, saline-alkali resistance, fruit development, and synthesis and accumulation of anthocyanins). The application of omics in recent years in comprehensive evaluation of, stress resistance, and biosynthesis and regulation of secondary metabolites was reviewed in this paper. The causes of differences in theresource quality and ecogeographic adaptation mechanism were also investigated, which provide scientific basis for theconservation, germplasm improvement and innovation and rational development and utilization of.

Murray; secondary metabolites; anthocyanin; anthocyanidin; multi-omics; medicinal resources

R282.71

A

0253 - 2670(2023)01 - 0272 - 11

10.7501/j.issn.0253-2670.2023.01.029

2022-10-09

國家自然科學基金面上項目（32170394）；中國科學院西部之光青年學者A類項目（2019）

李明珠，碩士研究生，研究方向為中藥資源化學及品質綜合評價。E-mail: lmz1559725@163.com

通信作者：劉增根，副研究員。Tel: (0971)6143747 E-mail: lzg2005sk@126.com

[責任編輯 崔艷麗]