植物代谢组学研究论文

摘要：从技术步骤、分析方法以及实际应用三个方面对当前药用植物代谢组学研究领 域的一些理论问题和实践中面临的挑战进行综述。 关键词：药用植物；代谢组学；功能基因组学 代谢组学是对生物体内代谢物进行大规模分析的一项技术［1］，它是系统生物学的 重要组成部分（如图 1 所示），药用植物代谢组学主要研究外界因素变化对植物所造成的 影响，如气候变化、营养胁迫、生物胁迫，以及基因的突变和重组等引起的微小变化，是 物种表型分析最强有力的工具之一。在现代中药研究中，代谢组学在药物有效性和安全性、 中药资源和质量控制研究等方面具有重要理论意义和应用价值。另外，在对模式植物突变 体文库或转基因文库进行分析之前，代谢组学往往是首先考虑采用的研究方法之一。目前， 国外已有成功利用代谢组学技术对拟南芥突变株进行大规模基因筛选的例子，这为与重要 性状相关基因功能的阐明和选育可供商业化利用的转基因作物奠定了基础。 图 1 系统生物学研究的四个层次略 目前，还有许多经济作物的全基因组测序计划尚未完成，由于代谢组学研究并不要求 对基因组信息的了解，所以在与这些作物有关的研究领域具有更大的利用价值，这也是其 与转录组学和蛋白组学研究相比的优势之一。代谢组学研究涉及与生物技术、分析化学、 有机化学、化学计量学和信息学相关的大量知识，Fiehn［2］对代谢组学有关的研究方向 进行了分类（见表 1）。 1 代谢组学研究的技术步骤 代谢组学研究涉及的技术步骤主要包括植物栽培、样本制备、衍生化、分离纯化和数 据分析 5 个方面（见图 2）。 1.1 植物栽培 对研究对象进行培育的目的是为了对样本的稳定性进行控制，相对于微生物和动物而 言，植物的人工栽培需要考 表 1 代谢组学的分类及定义略 虑更多的问题，如中药材在不同年龄、不同发育阶段、不同部位以及光照、水肥、耕 作等环境因素的微小差异都可引起生理状态的变化，而这些非可控及可控双重因素的影响 很难进行精确的控制，从而影响药用植物代谢组研究的重复性。为了解决以上问题，推荐 使用大容量的培养箱［3］，定时更换培养箱中栽培对象的位置，以及使用无土栽培技术 等，FukusakiE［4］利用无土栽培系统将水和养分直接引入植物根部，并且对供给量进 行精确地控制，大大提高了实验的重复性。 1.2 样本制备 为了获得稳定的实验结果，样本制备需要考虑样本的生长、取样的时间和地点、取样 量以及样本的处理方法等问题，并根据分析对象的分子结构、溶解性、极性等理化性质及 其相对含量大小对提取和分离的方法进行选择，逐一优化试验方案。MaharjanRP 等 ［5］用 6 种方法分别对大肠杆菌中代谢产物进行提取，发现用－40℃甲醇进行提取的效 果最好。现阶段代谢组学的分析对象主要集中在亲水性小分子，尤其是初级代谢产物，气 相色谱质谱联用质谱联用质谱质谱联用联用(GCMS)GCMS)质谱联用MS))和毛细管电泳质谱（质谱联用质谱质谱联用（CMS)E质谱联用MS)）联用都是分析亲水小分子的重要 技术。FiehnO 等［6］使用 GCMS)质谱联用MS) 对拟南芥叶片中的亲水小分子进行了分析，发现酒石 酸半缩醛、柠苹酸、别苏氨酸、羟基乙酸等 15 种植物代谢物。 1.3 衍生化处理 对目标代谢产物的衍生化处理取决于所使用的分析设备，GCMS)质谱联用MS) 系统只适合对挥发性 成分进行分析，高效液相色谱质谱联用法（HPLCMS)）一般则使用紫外或荧光标记的方法对样本进行 衍生处理，BlauK［7］对酯化、酰化、烷基化、硅烷化、硼烷化、环化和离子化等衍生方 法进行了详细的说明。然而离子化抑制常使得质谱质谱联用分析过程中目标代谢产物的离子化效率 降低，这主要是由于分离过程中污染物与目标代谢物难以完全分离开所引起的，优化色谱质谱联用 分离时间可有效缓解离子化抑制，然而在实际操作中不可能对上百种代谢产物的分离时间 进行优化，利用非放射性同位素稀释法进行相对定量可以很好的解决该问题。HanDK 等 ［8］应用同位素编码的亲和标记（ICMS)AT），根据经诱导分化的微粒蛋白及其同位素标记 物的峰面积比，对该蛋白的相对含量进行分析。ZhangR 等［9］发现同位素标记技术也 可用于代谢组学的研究，但是却存在许多困难。活体的同位素标记方法对于同位素的洗脱 是一种非常有潜力的技术，目前关于使用 34s 的研究已有报道［10］。 图 2 代谢组学研究技术步骤略 1.4 分离和定量 分离是代谢组学研究中的重要步骤，与质谱质谱联用联用的色谱质谱联用和电泳质谱（分析技术都是使用紫外 或电化学检测的方法进行定量，其对代谢组数据的分辨率与定量能力都有一定的影响。 TomitaM 等［11］总结了各种色谱质谱联用分离法中经常遇到的技术问题，认为毛细管电泳质谱（和气 相色谱质谱联用法由于具有较高的分辨率，已成为代谢组学研究的常规技术手段之一，液相色谱质谱联用因 其适用范围广，应用也相当广泛。 TanakaN 等［12］用高效液相色谱质谱联用对样品进行分离，认为使用硅胶基质填充毛细管 整体柱的高效液相色谱质谱联用系统具有用量少、灵敏性高、低压降高速分离等优势；同时， TolstikovV 等［13］也使用硅胶填充的毛细管液相色谱质谱联用方法对聚戊烯醇类异构体进行了 有效分离，获得了很好的分辨率。TanakaN 等［14］发现二维毛细管液相色谱质谱联用法的分辨 率比传统的高效液相法高 10 倍。相对于其他色谱质谱联用方法而言，超临界流体色谱质谱联用（S)FCMS)）是 分离疏水代谢物最具潜力的技术之一，特别适用于分离那些传统 HPLCMS) 难以分析的疏水聚 合物，BambaT 等［15］通过 S)FCMS) 对聚戊烯醇进行分析，证明其具有较好的分离能力。 针对质谱质谱联用中存在的共洗脱现象，HalketJM 等［16］发明了一种适用于 GCMS)质谱联用MS) 的反褶积 系统，对共洗脱的代谢产物进行分离与识别。AharoniA 等［17］使用傅立叶变换离子回 旋共振质谱质谱联用（FT质谱联用ICMS)R质谱联用MS)）对非目标代谢物进行分析，快速扫描植物突变样品，获得了一 定量的代谢成分。 与分离一样，定量能力也是代谢组学研究中的重要因素，其取决于各分析系统的线性 范围。傅立叶转换核磁共振（FT质谱联用NMR)、傅立叶红外光谱质谱联用（FT质谱联用IR）以及近场红外光谱质谱联用法 （NIR）等技术由于敏感性低，重复性受共洗脱现象影响较小也被用于检测中。近年来， FT质谱联用NMR 技术常被用于植物代谢组的指纹图谱质谱联用研究［18］，但由于 NMR 分析需要样品量 较大，分析结果易受污染，GriffinJL［19］发现将统计模式识别与 FT质谱联用NMR 相结合可以对 代谢物进行全面分析。除 FT质谱联用NMR 之外，FT质谱联用IR 通过对有机成分的结构进行常规光谱质谱联用测定， 也可适用于代谢组学的研究，特别是应用于构建代谢组学的指纹图谱质谱联用。尽管它不能对代谢 物进行全面分析，但对具有特定功能的组分却有很好的定量效果，对从工业及食品原材料 中分离的代谢混合物也可以进行全面分析，目前，已有学者将其成功地应用于拟南芥 ［20］和番茄［21］代谢产物指纹图谱质谱联用的研究中。 1.5 数据转换 为阐明代谢物复杂的线性或非线性关系，需要进行多变量分析，将原始的色谱质谱联用图数据 转换为数字化的矩阵数据，通过对色谱质谱联用峰鉴定和整合从而进行多变量分析。由于环境等因 素的干扰，光谱质谱联用数据需要通过适当的数据加工方法进行校正，包括：①降低噪声；②校正 基线；③提高分辨率；④数据标准化。JonssonP 等［22］报道了一种关于 GCMS)质谱联用MS) 色谱质谱联用 图数据处理的方法，可以对大量代谢产物样品进行有效的识别。 2 代谢组学中的数据分析方法 2.1 主成分分析法（PCMS)A） 主成分分析法，将实测的多个指标用少数几个潜在的相互独立的主成分指标线性组合 来表示，反映原始测量指标的主要信息。使得分析与评价指标变量时能够找出主导因素， 切断其他相关因素的干扰，作出更为准确的估量与评价。PCMS)A 数据矩阵通常来自于 GCMS)质谱联用MS)，LCMS)质谱联用MS) 或 CMS)E质谱联用MS)，因此将目标代谢产物作为自变量，而相应的代谢产物含量作 为因变量，定义与最大特征值方向一致的特征向量为第一主成分，依此类推，PCMS)A 便能通 过对几个主要成分的分析，从代谢组中识别出有效信息。主成分分析有助于简化分析和多 维数据的可视化，但是该方法可能导致一部分有用信息的丢失。 2.2 层次聚类分析法（HCMS)A） 层次聚类分析法也常用于代谢组学的研究中，它是将 n 个样品分类，计算两两之间的 距离，构成距离矩阵，合并距离最近的两类为一新类，计算新类与当前各类的距离。再合 并、计算，直至只有一类为止。进行层次聚类前首先要计算相似度（similarity），然后 使用最短距离法（NearestNeighbor）、最长距离法（FurthestNeighbor）、类间平均 链锁法（Between质谱联用groupsLinkage）或类内平均链锁法（Within质谱联用groupsLinkage）四 种方法计算类与类之间的距离。该方法虽然精确，但计算机数据密集，对大量数据点进行 分析时，更适合选用 K质谱联用均值聚类法(GCMS)KMCMS))或批次自组织映射图法(GCMS)BL质谱联用S)OM)，而 HCMS)A 适合 将数据转换为主成分后使用。2.3 自组织映射图法(GCMS)S)OM) 神经网络中邻近的各个神经元通过侧向交互作用相互竞争，发展成检测不同信号的特 殊检测器，这就是自组织特征映射的含义。其基本原理是将多维数据输入为几何学节点， 相似的数据模式聚成节点，相隔较近的节点组成相邻的类，从而使多维的数据模式聚成二 维节点的自组织映射图。除 PCMS)A 和 HCMS)A 外，S)OM 同样也可应用于包括基因组和转录组等 组学研究中［23］。最初 S)OM 计算时间长，依靠数据输入顺序决定聚类结果，近年来 S)OM 逐渐发展成为不受数据录入顺序影响的批次自组织映射图法(GCMS)BL质谱联用S)OM)。由于 BL质谱联用S)OM 可以对类进行调整，且有明确的分类标准，优化次序优于其他聚类法，已在基因 组学和转录组学数据分析中得到广泛的应用。 2.4 其他数据采矿方法 除 PCMS)A、HCMS)A 和 S)OM 外，很多变量分析方法都可用于植物代谢组学的分析。软独立 建模分类法(GCMS)S)IMCMS)A)是利用主成分模型对未知样品进行分类和预测，适合对大量样本进行 分析；近邻分类法(GCMS)KNN)和 K 平均值聚类分析法(GCMS)KMN)也可用于样品分类；主成分回归法 (GCMS)PCMS)R)或偏最小二乘回归法(GCMS)PLS))在某些情况下也可使用。然而到目前为止由于还没有建立 一个标准的数据分析方法，代谢组学仍然是一门有待完善的学科。 3 代谢组学在药用植物中的实践 植物药材来源于药用植物体，而药用植物体的形态建成是其体内一系列生理、生化代 谢活动的结果。植物代谢活动分为初生代谢和次生代谢，初生代谢在植物生命过程中始终 都在发生，其通过光合作用、柠檬酸循环等途径，为次生代谢的发生提供能量和一些小分 子化合物原料。次生代谢往往发生在植物生命过程中的某一阶段，其主要生物合成途径有 莽草酸途径、多酮途径和甲瓦龙酸途径等。植物药材含有的生物碱、胺类、萜类、黄酮类、 醌类、皂苷、强心苷等活性物质的绝大多数属于次生代谢产物，因此探讨次生代谢产物在 药用植物体内的合成积累机制及其影响因素，对于提高活性物质含量、保证药材质量、稳 定临床疗效等具有重要意义。孙视等［24］通过对银杏叶中黄酮类成分积累规律的研究， 提出了选择具有一定环境压力的次适宜生态环境解决药用植物栽培中生长和次生产物积累 的矛盾。王昆等［25］以人参叶组织为材料,总结了构建人参叶 cDNA 文库过程中存在的 一些关键问题和应采取的对策,为今后关于人参有效成分如人参皂苷的生物合成途径及其调 控的基础研究提供技术参考和理论指导。最近，美国加利福尼亚大学伯克利分校的 Keasling 等［26］采用一系列的转基因调控方法，通过基因工程酵母合成了青蒿素的前 体物质——青蒿酸，其产量超过 100mg/L，为有效降低抗疟药物的成本提供了机遇。经 过长期的研究积累，人们对代谢途径的主干部分（为次生代谢提供底物的初生代谢途径） 已经基本了解，例如酚类的莽草酸途径，萜类的异戊二烯二磷酸(GCMS)IPP)途径等。被子植物中 一些相对保守的次生代谢途径也得到了很好的研究，如黄酮类、木质素的生物合成与调控。 然而，对次生代谢最丰富最神奇的部分——特定产物合成与积累的过程，还所知甚少 ［27］。 4 展望 近年来，代谢组学正日益成为研究的热点，越来越多的人已加入到代谢组学的研究中。 随着代谢组学积累的数据和信息量的增大，其在药用植物学各个领域的应用价值也与日俱 增。它将不仅能对单个代谢物进行全方面的分析，更能寻找其代谢过程中的关键基因、通 过代谢指纹分析对药用植物进行快速分类、进一步研究药用植物有效成分代谢途径以及环 境因子对植物代谢和品质的影响与调控机制。 然而依据传统中医药学和系统生物学的指导思想，目前急待解决的是中药种质资源的 代谢组学研究和中药体内作用的代谢组学研究。同时，代谢组学在分析平台技术、方法学 手段和应用策略等方面相对于其他组学技术还需要进一步发展和完善，还需要其他学科的 配合和介入。相信随着更有力的成分分析设备的使用及代谢组数据库的建立，药用植物代 谢组学将对中医药学产生深远的影响。