黑果枸杞 Lycium ruthenicum Murray为茄科（Solanaceae）枸杞属 Lycium L.多棘刺灌木（图1），其果实味甜多汁，民间常生食或榨汁做饮料，具有滋补强壮、明目及降压作用，其根可以治咳嗽、哮喘、感冒和发烧，叶民间当茶泡饮[1]。作为我国西北盐碱、荒漠地区具有重要生态和经济价值的药食兼用资源植物，黑果枸杞含有花色苷、花青素、多酚、黄酮、多糖、生物碱、氨基酸、矿质元素、维生素等多种成分，具有补肾益精、养肝明目、生津止渴、润肺止咳、补血安神等功效[2-3]，尤以其主要次生代谢物花青素类化合物已被广泛应用于医药保健品、食品添加剂、化妆品等领域，并在我国传统民族医药（藏药、蒙药、维药）中占重要地位[4]。
黑果枸杞所处的独特生态地理环境，如生长环境纬度高、盐碱化程度高、紫外线强、昼夜温差大、日照时间长等，都有助于次生代谢产物（如花青素、多酚、黄酮、多糖等）的转化与积累[2,5-7]。近年来，资源过度开发、质量不稳定、产地来源混乱和有效成分地理变异显著等问题已经明显制约了黑果枸杞产业的发展[8-9]，因此亟待寻找黑果枸杞品质差异形成原因以及构建资源（道地性）生态环境评价体系。
随着高通量测序技术的出现、色谱/质谱技术的快速发展以及分子生物学新方法的开发[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]。
目前，黑果枸杞因独特的生物活性和巨大的生态经济价值使其成为国内外学者竞相研究与开发的对象[1-2]，在保健品、食品、化妆品等行业需求量较大，资源的过度利用正严重威胁其野生种质资源。此外，随着整个生态环境日趋恶化和人为过度采摘砍伐，野生资源正在大面积的减少，部分地区甚至出现枯竭[8]。
因此，开展黑果枸杞居群遗传多样性和分子评价的研究有利于其野生种质资源的保护和可持续利用，并为黑果枸杞的科学育种及种质创新提供依据。本课题组基于前期调查与研究发现，同一产区种植和野生黑果枸杞的主要成分（次生代谢物和营养成分）种类一致，且含量差异不显著[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个。
另外，基因 accD 编码乙酰辅酶A羧化酶的CT亚基，是乙酰辅酶A羧化酶表达水平的限制因子；而黑果枸杞的 accD 基因有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）物种为外群，系统发育分析表明，黑果枸杞与茄科植物颠茄 Atropa belladonna L.、马尿泡 Przewalskia tangutica Maxim.和 Scopolia parviflora (Dunn) Nakai聚在一起，具有密切的亲缘关系。
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个显著差异表达的结构基因（ PAL 、 C4H 、 4CL 、 CHS 、 CHI 、 F3H 、 F3’H 、 F3’5’H 、 DFR 、 ANS 、 UFGT ）可能与花青素生物合成相关。
Li等[41]基于Illumina HiSeq测序平台对黑果枸杞5个不同发育阶段果实颜色进行转录组分析，结果显示，在转录组*共中**鉴定出31 302个差异基因，72个关键候选基因对参与类黄酮和花青素生物合成的14种酶产生直接影响。朱雪冰等[42]在黑果枸杞中测得1.167×10^8和9.695×10^7 reads，筛选出了13个与花青素生物合成有关的基因，克隆黑果枸杞中 MYB 基因的转录因子Sl AN2，该基因全长774 bp，编码257个氨基酸，相对分子质量为29 775.84，等电点为7.79。
严莉等[43]基于转录组数据分析 MYB 基因在黑果枸杞果实不同发育期的差异表达模式，结果表明， MYB 基因可能参与了果实不同发育时期花青素变化的调控，利用荧光定量PCR的差异表达数据进一步验证了部分MYB转录因子在果实不同发育时期的花青素合成中可能起到调控作用；同时在黑果枸杞 MYB 基因家族注释得到83个MYB类转录因子基因，为进一步研究MYB家族的基因结构和生物学功能奠定了基础。
彭勇等[44]对黑果枸杞3个不同发育时期的果实进行转录组测序分析表明，京都基因与基因组百科全书（KEGG）信号通路富集显示亮氨酸生物合成途径在变色期对比青果期有7个差异表达基因，黑熟期对比变色期有35个差异表达基因，且全部为上调表达；在糖酵解途径中3个不可逆的反应步骤限速酶基因在黑熟期均为上调表达。
孟小伟等[45]对黑果枸杞果实不同发育过程进行转录组测序，共得到43 573条Unigene，有23 723条Unigene在基因本体（GO）、KEGG、真核同源群簇（KOG）、非冗余蛋白（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个不同发育阶段花青素生物合成相关的基因序列，分离出参与黑果枸杞花青素生物合成的结构基因（ C4H 、 CHS 、 CHI 、 F3’H 、 F3’5’H 、 DFR 和 ANS ），从转录组数据库中检索到编码 4CL 、 F3’H 、 UF3GT 、 GST 、 AT 、 AN2 、 JAF13 、 AN1B 和 AN11 的序列，并应用Mega 4.0对相关基因进行系统发育分析。上述基因在果实发育过程中的表达谱表明，在果实中转录丰度与花青素积累之间存在显著的正相关； AN2 和 JAF13 或 AN1b 可能共同调控 DFR 、 ANS 、 UF3GT 、 GST 和 AT 的转录，从而调控黑果果实中花色苷的合成。此外，调控基因和结构基因的表达模式以及分支节结构基因 F3’5’H / F3’H 的转录比率可能决定了黑果枸杞和宁夏枸杞果实之间花青素生物合成的表型差异。
Zong等[48]研究 AN2 等位基因存在与果实颜色的关系，根据转录组分析结果从黑果枸杞和宁夏枸杞中分离得到编码MYB转录因子的 LrAN2 和 LbAN2 。在系统发育树中， LrAN2 和 LbAN2 与调节*草烟**花青素合成的 NTAN2 关系密切，而 LrAN2 的过度表达诱导*草烟**所有组织中花青素的生物合成，且明显强于 LbAN2 。此外， AN2 转录因子仅在黑果果实中检测到，并在果实发育过程中增加，伴随着花色苷的积累；而 LrAN2 的功能多样性和高表达水平可能是黑果枸杞果实花青素含量高的原因。上述研究结果为黑果枸杞资源改良与利用提供了重要基因资源，为花青素的生物合成与代谢提供了坚实的理论基础。
4 蛋白质组学在黑果枸杞研究中的应用
蛋白组学是以蛋白质组为研究对象，研究细胞、组织或生物体蛋白质组成及其变化规律的科学[49]，它主要利用高效液相色谱-串联质谱法技术获取某一时刻生物机体内所有表达蛋白，揭示生物体已经发生什么。蛋白组学是基因表达调控机制研究的一个重要工具。当前，随着高精确度、高分辨率、高灵敏度、高通量的生物质谱仪器的快速发展，以及各物种基因数据库、蛋白质数据库的不断完善，配合基于复杂算法的生物信息学工具，蛋白质组学的研究内容、研究方法和研究领域都发生了翻天覆地的变化[50]。
目前的蛋白质组学研究方法有双向凝胶电泳技术、质谱技术、同位素的相对和绝对定量标记技术、蛋白质芯片技术、酵母双杂交系统技术和免疫共沉淀等技术；其中，同位素标记相对和绝对定量（iTRAQ）、串联质谱标签（TMT）是最常用的基于同位素标记质谱分析的定量蛋白质组学分析技术，可对多达16（8）个不同样本进行蛋白定量比较，标记效率和质谱检测灵敏度高，适用范围广泛[51]。
同一生物体在不同生长时间和生长环境下，基因表达情况并不完全相同，进行多器官、多时期、多维度、多水平的药用植物蛋白质组学研究，建立药用植物功能基因库，可为药用植物种质资源评价、中药材质量控制和标准建立、新品种选育、有效成分的合成生物学等领域提供分子基础，推动药用植物产业的健康可持续发展[52]。
目前，蛋白质组学技术在植物的研究领域使用广泛，然而关于黑果枸杞蛋白质组学的研究较少。汪亚娟等[53]采用iTRAQ技术鉴定出862个黑果枸杞蛋白质，经100 mmol/L盐处理7 d后黑果枸杞检测出165个表现差异的蛋白质，其中有72个蛋白质上调，93个蛋白质下调。并了解到盐胁迫主要影响了黑果枸杞的能量代谢与离子运输，其中参与离子运输的差异蛋白质共4个，全都是上调的膜孔蛋白；参与能量代谢的差异蛋白质共有13个，包括9个上调蛋白和4个下调蛋白。
此外，黑果枸杞受到盐胁迫时有利于黄酮类化合物的大量合成，并且以增加对香豆酸进入木质素生物合成途径的方式来促进黄酮类化合物的积累，从而缓解盐胁迫带来的氧化损伤。因此，体内黄酮类化合物的积累可能是盐生植物黑果枸杞对抗高盐环境的机制之一，而且该研究鉴定出的一些有价值的蛋白质，为进一步研究盐生植物耐盐机制和为作物耐盐驯化提供了新的依据。
5 代谢组学在黑果枸杞研究中的应用
代谢组学是对内源性及外源性小分子代谢产物进行定性和定量检测，并寻找与机体生命活动、病理变化相对应的代谢物质的研究方法[54]。代谢组学揭示了生命体活动事件确实发生了，它根据研究目的不同分为靶向代谢组学和非靶向代谢组学，两者各有各的优势和特点。靶向代谢组定性定量准确，但对物质的覆盖率有限；而非靶向代谢组对物质的覆盖率广泛，但缺乏绝对定性定量的数据[55]。
代谢组学是对基因组学和蛋白质组学的补充，常见的分析手段包括液质谱联用（LC-MS）、气质谱联用（GC-MS）、超高效液质联用（UPLC-MS）、毛细管电泳质谱联用（CE-MS）、核磁共振（NMR）等，四极杆线性离子阱（Q-Trap）、四极杆飞行时间（Q-TOF）及离子阱飞行时间（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种花青素转运机制，即谷胱甘肽 S -转移酶介导的转运、多耐药相关蛋白和多药和有毒化合物外排转运蛋白介导的转运以及膜囊泡介导的转运。
李丽等[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转录因子 LrAN2 ，其表达量随果实的发育而增加，类黄酮物质含量也随之增加，而该基因在宁夏枸杞中几乎不表达[9,48]。自然群体中 LrAN2 的等位变异与黑色果实性状紧密关联，在*草烟**里过量表达 LrAN2 ，除根部外，其他组织部位均表现出紫黑色。通过多组*联学**合分析，初步发现 LrAN2 转基因*草烟**中185个基因（146上调、39个下调）产生差异表达，其中类黄酮生物合成代谢途径差异最显著。
代谢组学可以批量获得基因表达的最终产物，反映了基因和蛋白表达的微小变化以及机体系统的生理和病理状态；蛋白质组学能够动态地描述基因调节，可获取某一时刻机体内所有表达蛋白。因此，将蛋白质组学和代谢组学进行关联分析，对基因表达的蛋白质和代谢物进行定量测定和差异比较，有助于进一步揭示基因表达调控网络和次生代谢物的生物合成代谢机制。在黑果枸杞耐盐机制研究中，学者利用蛋白质组学和代谢组*联学**合分析方法，分析了黑果枸杞愈伤组织在不同浓度盐胁迫下的蛋白质和化学成分，阐释了差异性蛋白质基因表达主要富集于离子运输和能量代谢途径，并且盐胁迫有助于植株黄酮类化合物的合成与积累，缓解逆境带来的氧化损伤[53]。
尽管目前组学技术已经在黑果枸杞中得以广泛应用，尤其集中在抗逆（耐盐）和果实发育过程中的花青素合成代谢（表1），但是研究学者们对其进行的多组*联学**合研究还为之甚少。显然，单凭一种组学技术已经无法满足黑果枸杞适应与进化、花色苷差异成因解析等相关深层次的研究，因此，基于多组*联学**合分析阐明黑果枸杞次生代谢物（花青素、黄酮、生物碱等）合成调控机制、揭示黑果枸杞耐旱和抗盐碱机制将是今后枸杞属植物资源的研究重点。
7 结语及展望
随着多组学和相应分子生物学技术的发展与渗透，药用植物资源品质形成的探索不论在生态环境等外在因素，还是基因、蛋白质与遗传分化的分子水平都有了一定的推进[12,60-61]。基于多组*联学**合分析是药用植物次生代谢产物合成途径和资源品质差异成因研究的重要技术[62-63]。在进入大数据时代的背景下，基于系统生物学研究模式和多组学研究思路，能够为药用植物次生代谢物生物合成代谢调控、资源品质综合评价以及种质创新等相关科学问题的解决打下坚实的理论基础。
从已有的研究成果可以看出，各组学应用在黑果枸杞资源研究中已经取得了一定进展，为资源的品质差异、次生代谢产物合成途径、有效成分确定以及生境适宜性等研究提供了重要参考。然而，有关黑果枸杞资源评价、组学研究、品质差异及其与生态环境关系方面的研究仍然存在一些不足。主要表现在：（1）很少有研究关注不同产区黑果枸杞次生代谢产物差异机制及生态环境对有效成分合成积累的影响，资源生境适宜性系统评价和区划研究基本处于空白。（2）前期有关黑果枸杞组学研究主要集中在叶绿体基因组学和转录组学，缺乏蛋白质组和代谢组及其组学间的整合分析，组学间数据交集的挖掘深度不够，且未实现多组学的真正融合。
综上所述，随着后续黑果枸杞全基因组数据的公布，并基于化学生物学分析和多组*联学**合应用，鉴定黑果枸杞次生代谢物生物合成过程中的关键酶和调控因子，深入解析主要有效成分花色苷合成调控机制（图3），结合生态环境因子的相关性分析，将对探究黑果枸杞花色苷差异成因、资源品质形成过程和生境适宜性区划评价具有重要的理论指导作用，从而促进黑果枸杞资源健康可持续开发与利用。
利益冲突 所有作者均声明不存在利益冲突
参考文献（略）
来 源：李明珠，刘增根．组学技术在黑果枸杞资源研究中的应用 [J]. 中草药, 2023, 54(1):272-282.