基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析.pdf

天然产物研究与开发NatProdResDev2023,35:1305-1313,1392基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析吴姣姣,刘东,张莉,夏佳,杨成都中医药大学药学院，成都6 11137摘要：建立白术挥发油指纹图谱,并结合网络药理学对白术挥发油质量标志物（Q-Marker)进行初步预测分析。采用气相色谱-质谱联用技术（GC-MS)对白术挥发油进行成分分析,对11批白术饮片挥发性成分进行相似度评价、主成分分析（PCA)和聚类分析（HCA）,筛选挥发油中的Q-Marker候选成分；对筛选的Q-Marker候选成分进行靶点搜集和网络药理学分析，构建“成分-靶点-通路”网络，并预测白术挥发油的Q-Marker。结果显示11批白术饮片挥发油的指纹图谱相似度均大于0.9 9 0,提示不同产地和批次的白术饮片挥发油成分整体组成上具有一致性;PCA分析显示浙江产区白术分布更集中,HCA分析显示当分类距离为15时，浙江产区的4批白术被聚为一类，这表明道地产区的白术挥发油质量更稳定；根据GC-MS共有成分得到白术挥发油中10 个Q-Marker候选成分，网络药理学分析提示这些成分涉及19 个关键靶点和6 6 条通路。本研究初步预测棕榈酸、-桉叶醇、-榄香烯、榄香醇和苍术酮是白术挥发油的Q-Marker,为白术的质量控制提供了理论依据，也为其他类似中药材的质量标志物研究提供了理论参考。关键词：白术；挥发油；指纹图谱；网络药理学；质量标志物中图分类号：R932D01:10.16333/j.1001-6880.2023.8.003Predictive analysis of quality markers of Atractylodis Macrocephalae Rhizomavolatile oil based on GC-MS fingerprint and network pharmacology雨，唐飞，陈东璐，敖慧*文献标识码：A文章编号：10 0 1-6 8 8 0(2 0 2 3)8-130 5-10WU Jiao-jiao,LIU Dong,ZHANG Li,XIA Jia,YANG Yu,TANG Fei,CHEN Lu,AO HuiCollege of Pharmacy,Chengdu University of Traditional Chinese Medicine,Chengdu 611137,ChinaAbstract:To analyze and predict the potential quality markers of Atractylodis Macrocephalae Rhizoma volatile oil(AMO)based on fingerprint and network pharmacology,11 batches of AMO was analyzed by gas chromatography-mass spectrometry(GC-MS),and similarity evaluation,principal component analysis(PCA)and cluster analysis were carried out,and the Q-Marker candidate components in the volatile oil were screened.Target collection and network pharmacological analysis werecarried out to build a component-target-pathway network and predict the Q-Marker of AMO.The fingerprints of AMO wereestablished with similarities greater than 0.990,indicating that consistency in the composition of AMO from different originsand batches.PCA analyses showed a more concentrated distribution of AMO from the Zhejiang,and HCA analyses showed thatwhen the classification distance was 15,four batches of AMO from the Zhejiang region were clustered into one group,indica-ting that the volatile oil quality of AMO from the Daoji region was more stable.Ten Q-Marker candidates were obtained afterGC-MS analysis,19 key targets and 66 pathways were obtained from the network pharmacological analysis,and palmitic acid,-eudesmol,y-elemene,elemene and atractylone were the Q-Markers of AMO according to the preliminary prediction,provi-ding a theoretical reference for quality control of Atractylodis Macrocephalae Rhizoma and other Traditional Chinese herbs.Key words:Atractylodis Macrocephalae Rhizoma;volatile oi;fingerprint;network pharmacology;quality markers中药材普遍存在多基原、多产地等复杂情况,质收稿日期:2 0 2 3-0 1-18接受日期：2 0 2 3-0 4-2 5基金项目：国家青年科学基金（8 150 32 7 2）；四川省教育厅项目(23NSFSC2057)*通信作者E-mail:量差异较大，函待建立科学合理的中药质量保障体系。目前多数中药材质量控制指标选择多为其含有的化学成分,但其与药物功效关联不足 1。基于此,刘昌孝院士提出用“中药质量标志物”的概念 2 。中药质量标志物是与中药的功能属性密切1306相关的化学物质，反映了中药质量控制的发展方向 1-3。故研究中药质量标志物是提升我国中药质量控制水平,为临床提供安全有效中药产品的重要手段。白术是常用大宗中药材,具有健脾燥湿、利水止汗等功效,广泛应用于中医临床。然而,目前白术并没有业内公认的质量标志物 4，这将影响白术的质量提升和进一步开发利用。现有研究表明挥发油和多糖是白术的主要有效成分 5.6 。多糖由于结构特殊，目前很难作为质量标志物或质量控制成分 7 。研究表明白术挥发油具有调节胃肠道、修复胃黏膜、抑菌、保肝等药理作用，与白术药理作用基本类似 8.9 ,并且可用GC-MS等技术实现定性定量检测 10 。故从白术挥发油中筛选白术质量标志物是可行有效的思路。但白术挥发油受产地、土壤、气候、采收等因素影响,化学成分种类差异较大 10.1,因此有必要建立白术挥发油GC-MS指纹图谱寻找其共有化学特征，初步确定候选质量标志物。网络药理学具有系统性和整体性的特点,是一种建立在系统生物学和网络生物学基础上的新兴研究方法，已被广泛应用于单味药及复方药治疗疾病的物质基础和多靶点作用机制预测 12 。将指纹图谱与网络药理学相结合寻找中药质量标志物的方法已成功应用于江积壳、连翘、白花蛇舌草、山银花等中药的质量标志物预测 12-17 ,但尚未见GC-MS指纹图谱和网络药理学相结合的方法应用于白术挥发油质量标志物预测分析。因此,本研究采用GC-MS对白术挥发油成分进行分析和指纹图谱建立，并结合利用网络药理学方法对白术挥发油成分、靶点、通路进行系统性分析,并预测白术挥发油的潜在质量标志物，旨在为全面控制和评价白术片质量提供科学依据和理论参考。1材料与方法1.1试验材料1.1.1药材本研究11批白术药材（编号S1S11)均购于成都国际商贸城荷花池中药材市场,经成都中医药大学陈璐副教授鉴定为菊科植物白术Atractylodesmacrocephala Koidz.的干燥根茎，样品信息详见表1。1.1.2试剂正已烷（分析纯，成都市科隆化学品有限公司，批号：2 0 17 112 40 1）；无水硫酸钠（分析纯，成都市科隆化学品有限公司，批号：2 0 2 0 0 3310 1）；水为蒸馏水。天然产物研究与开发编号No.S1S2S3S4S5S61.1.33仪器轻油提取器（四川蜀牛玻璃仪器有限公司）；Agilent气相色谱-质谱联用仪、Agilent GC-MSD工作站（7 8 9 0 A-5975C，美国安捷伦公司）；Sartorius电子天平(BP211D AG,德国赛多利斯公司）;高速中药粉碎机（HX-200型，浙江省永康市溪岸五金药具厂。1.2实验方法1.2.1白术挥发油提取采用水蒸气蒸馏法提取白术挥发油,称取白术粗粉8 0 g，装人10 0 0 mL圆底烧瓶，加人8 倍量水浸泡3h，电热套加热至微沸，回流提取6 h经无水硫酸钠脱水后,收集挥发油于已称重的EP管中,再次称重,密封后置于4冰箱保存。按公式(1)计算挥发油含量。挥发油提取率=挥发油质量（g)/药材质量（g）10 0%1.2.2供试品溶液制备取挥发油30 L溶于1mL正已烷中,微孔滤膜过滤，滤液装人进样瓶,为GC-MS测定样品。1.2.3GC-MS 条件色谱条件：色谱柱为Agilent19091S-433色谱柱(30 m250 m0.25 m）;载气为氮气（He）,载气速度1mL/min;采取程序升温，初始柱温6 0，保持2 min，以5/min升至16 0,保持5min,再以4/min升至2 0 0,保持8 min,然后以6/min升至2 50，保持10 min；分流进样，分流比为20:1;进样量:1 L;溶剂延迟:3 min。质谱条件：电子轰击（EI）离子源，电子能量为70eV,离子源温度：2 30，四级杆温度:150，数Vol.35表1样品信息Table 1 Sample information采集地Place of collection安徽安徽安徽安徽浙江浙江编号No.Place of collectionS7浙江S8浙江S9四川S10云南S11江西(1)采集地Vol.35据采集扫描模式为全扫描。1.2.4定性定量方法各色谱峰对应的质谱图经联用仪的计算机谱库检索进行定性,质谱数据为NIST14.0谱库。各组分的相对含量是根据总离子流图由计算机采用峰面积归一化法计算获得。2白术挥发油的GC-MS 分析结果2.1白术挥发油得率样品S1S11挥发油得量及提取率见表2。2.2GC-MS 指纹图谱建立取11批白术挥发油，按“1.2.2”项下条件制备供试品溶液,按1.2.3”项下GC-MS条件依次进样，并记录色谱峰面积,11批白术挥发油指纹图谱见图1。吴姣姣等：基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析挥发油质量No.Volatile oil weight(g)S10.193S20.145S30.193S40.290S50.145S60.386S70.193S80.289S90.195S100.241S110.2861307表2 11批白术药材挥发油提取率Table 2Extraction rate of 11 batches of AMO编号0.2410.1810.2410.3620.1810.4830.2410.3610.2440.3010.358提取率Extraction rate(%)456783.0913.09图111批白术挥发油指纹图谱和共有峰标记Fig.1 The fingerprints and common peak markers of 11 batches of AMO2.3指纹图谱共有峰标记及相似度评价将11批白术挥发油总离子流图导人中药指纹图谱相似度评价系统,以 S8为参照图谱,采用平均数法,时间宽度设为0.1min,确定了10 个共有峰，峰序及所对应的化合物见图1和表2。以共有成分的色谱峰为参照进行多点校正,对11批白术挥发油序号tRNo.(min)116.47219.57295910-36-4328.4776989-21-51023.0933.09时间Time（mi n)整体相似度进行评价,结果显示11批白术样品相似度均大于0.9 9 2,表明11批白术样品挥发油相似度较高,质量较为稳定，但各样品挥发油中同一色谱峰的峰面积不同,说明所含成分的含量存在一定差异，如浙江产白术的挥发油中苍术酮整体含量高于其他产区,榄香醇含量整体低于其他产区,详见表3、表4。表311批白术挥发油共有成分信息及相对含量Table 3Common components and relative content of 11 batches of AMO分子式化合物CASMolecularformula29873-99-2Cis H24Cis H24IsophoreneCis H26 0苍术酮 AtractyloneS11S10S9S8S7S6S5S4S3S2S143.0953.09相对含量 Relative content(%)CompoundS1y-榄香烯-Elemene3.4523.713.057 1.807 5.982 8.223 5.644 3.028 2.813(-)-异喇叭烯(-)-2.416 6.429 4.482 7.928 3.803 4.715 4.273 7.152 2.166 8.16 2.03760.378 58.05 58.818 61.416 72.168 64.212 64.044 65.457 53.09 60.449 62.82963.09S2S373.09S4S5S6S7S8S94.042.493S10S111308续表 3(Continued Tab.3)序号tRNo.(min)431.493532.088632.31618594-05-3733.475201358-24-9835.4465101628-22-2937.562101527-74-61051.868天然产物研究与开发分子式化合物CASMolecularformulaCis H2%0473-15-4639-99-657-10-3Vol.35相对含量 Relative content(%）CompoundS1-桉叶醇-Eu-2.7564.6852.2290.3892.3663.5552.3251.1632.7511.4322.694desmolCis H26 0榄香醇 Elemene4-环已基苯乙酮 4-Cy-Ci4H1s 0clohexylacetophenone24-王二烯-3,12-二烯C29 H4624-Noroleana-3,dieneCis H24 0缬草醇Valerenol反式戊烯醇醋酸酯Ci7 H26 02trans-Valerenylace-tateC16 H32 02棕榈酸 Palmitic acidS21.56 0.824 1.663 1.07 0.658 0.422 0.908 0.808 1.448 1.066 1.7331.4431.7521.841 2.115 0.9760.611.482 1.1771.51.2581.59112-8.4415.6264.765 5.725 3.6224.7274.1473.635.8243.7395.031.518 0.832 1.291 0.548 0.926 0.828 1.141 0.4922.566 0.5971.7785.186 4.5926.784.103 3.548 3.341 5.249 3.365 5.551 3.998 6.0243.275 3.696 5.284 2.854 1.6971.391.881 2.983 8.285 3.192 6.76S3S4S5S6S7S8S9S10S11表411批白术挥发油GC-MS指纹图谱相似度评价Table 4Similarity of GC-MS fingerprints of 11 batches of AMO编号No.S1S2S3S4S5S62.4主成分分析采用帕诺米克分析平台（http：/v 2.b i o d e e p.相似度编号SimilarityNo.0.997S70.998S80.998S90.998S100.998S110.997相似度Similarity0.9990.9990.9920.9990.996cn/)对11批白术挥发油GC-MS共有成分含量信息16,17 进行PCA主成分分析,相应得分结果见图2。由PCA分析可见，批次S2、S4、S10 分布于第一象限,S1、S3、S9 分布于第二象限,S11 在第三象限,S5S8 浙江产白术挥发油主要分布在第四象限。3Group回AC?N0535图2 11批白术挥发油PCA得分图Fig.2PCA score of the 11 batches of AMO注：A：安徽;C：四川;J：江西;Y：云南;Z:浙江，下同。Note:A：A n h u i;C:Si c h u a n;JJi a n g x i;Y：Yu n n a n;Z:Z h e j i a n g,t h e s a me b e l o w.2.5聚类分析采用帕诺米克分析平台（http:/ 16.17 5PC1(39.8%)进行系统聚类分析,聚类结果见图3。由聚类分析可见当分类距离为15时,样品被分为2 类，S4、S5S8 及 S10 为第一类,S1 S3、S9、S11 为第二类。Vol.35吴姣姣等：基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析1309NO:S10385888858S1120151050图311批白术挥发油聚类分析结果Fig.3Cluster analysis results of the 11 batches of AMO3Q-Marker成分-靶点-通路网络药理学分析3.1Q-Marker候选成分选择根据GC-MS指纹图谱确定的10 个共有成分，采用峰面积归一化法计算得到共有成分分别占S1S11批次白术样品挥发油总量的9 0.42%、90.2%、9 0.2 1%、8 7.9 6%、9 5.7 5%、9 2.0 2%、91.09%、8 9.2 6%、8 5.9 9%、8 7.9 3%、9 2.9 7%。表明这10 个共有成分在白术挥发油中的代表性较强 15-18 ,可作为白术挥发油Q-Marker的候选成分。3.2成分靶点预测利用 Pub Chem（h t t p s:/p u b c h e m.n c b i.n l m.nih.gov/）、T CM SP（h t t p s:/o l d.t c ms p-e.c o m/t c m-sp.php)、Sw i s s T a r g e t Pr e d i c t i o n(h t t p:/w w w.s w i s-stargetprediction.ch/）和 PharmMapper(http:/www.lilab- Target Prediction数据库中取Probability大于0 的靶点。除2 4-noroleana-3，12-diene 无符合条件靶点外,共搜集了9 个白术挥发油成分靶点,除去重复靶点后使用uniprot转化为基因名，得到17 3个相关靶点。RHPATG5XRTNESRARARLOXH图4靶点PPI网络图Fig.4 Target PPI network diagram13103.3PPI网络分析将获得的17 3个相关靶点导人STRING11.5（h t t p s:/c n.s t r i n g-d b.o r g/）数据库中,选择“Homosapiens物种,去除单一节点,导出结果为.tsv格式，并将建立的蛋白相互作用关系导入Cytoscape3.9.1靶点简称Abbreviation of targetTNFPPARGPPARAESR1PTCS2RXRAPTENMAPK1IL10NCoA2NOS3ABCB1CYP2C19SLC6A4CYP19A1FABP1GRM5NROB2PGR天然产物研究与开发软件进行可视化分析，得到蛋白互作（protein-proteininteraction,PPI)网络，结果见图4（图中degree值越大，节点越大，颜色越深)；选取网络中degree值大于2 倍中位数的靶点作为关键靶点 16.17 ，结果共获得19个关键靶点，结果见表5。表5关键靶点筛选结果Table 5KKey target screening results靶点全称Full name of target肿瘤坏死因子Tumor necrosis factor过氧化物酶体增殖物激活受体Peroxisome proliferator-activated receptor过氧化物酶体增殖物激活受体Peroxisomeproliferator-activatedreceptor前列腺素G/H合成酶2 Prostaglandin G/H synthase 2视黄醇 X受体 Retinoid X receptor 磷脂酰肌醇-3,4,5-三磷酸3-磷酸酶和双特异性蛋白磷酸酶PTEN Phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase and dual-specificity protein phosphatase PTENMAP激酶 ERK2 MAPkinase ERK2白细胞介素-10 Interleukin-10核受体辅激活因子2 Nuclear receptor coactivator2内皮一氧化氮合酶Nitric-oxide synthase，e n d o t h e l i a lP-糖蛋白 1 P-glycoprotein 1细胞色素P4502C19CytochromeP4502C19血清素转运体 Serotonin transporter细胞色素P45019A1CytochromeP45019A1肝脏脂肪酸结合蛋白Liver fatty acid-binding protein代谢型谷氨酸受体5Metabotropic glutamate receptor5核受体亚家族0 B组成员2 Nuclear receptor subfamilyO group Bmember 2孕激素受体Progesterone receptorVol.35度值Degree544639雌激素受体 Estrogen receptor 373429272725252422222222212019193.4CGO与KEGG通路富集分析利用DAVID数据库（https:/david.ncifcrf.gov/)对获得的17 3 个靶点进行基因本体（gene on-tology,GO)富集分析。选择“Homo sapiens物种下的“Gene_Ontology进行富集分析,根据P0.05对富集结果进行筛选，得到397 个生物过程（biologicalprocess,BP）条目,55个细胞组成（cellular compo-nent,CC）条目,和116 个分子功能（molecular func-tion,MF)条目,按P值从小到大选取排名前2 0 的条目,得到图5。如图5所示,BP主要与药物反应、G蛋白偶联受体、突触传递、激素信号传导、脂质代谢等有关；CC主要涉及细胞质、细胞膜、内质网膜组成、染色质等;MF主要与RNA转录因子活性、转录因子结合、锌离子结合、DNA结合、G蛋白偶联受体活性、氧化还原酶活性等相关。利用DAVID数据库进行KEGG_PATHWAY富集分析,并根据P0.05,筛选得到KEGG信号通路66条,按P值从小到大选取排名前2 0 的通路绘制得到图6。图6 表明白术挥发油主要通过神经活性配体-受体相互作用、代谢途径、化学致癌-受体激活、血清素能突触、氮代谢、PPAR信号通路、花生四烯酸代谢、类固醇激素生物合成、钙信号通路、胆汁分泌、胆碱能突触等发挥作用。Vol.35吴姣姣等：基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析1311BPCCMF70一POPRuodturooonduosuenayuseidopusnbaL3XplounaprxoLLOoouanbastoosuen5Biological-processCellular-component图5GO功能富集分析结果Fig.5GO function enrichment analysis resultsMolecular-functionMetabolic pathwaysNeuroactiveligand-receptorinteractionChemical carcinogenesis-receptoractivationCalcium signaling pathwayPPARsignalingpathwaySerotonergic synapse-logio(pvalue)Estrogen signaling pathway10.0Cholinergic synapse7.5Arachidonicacidmetabolism5.0Thl7celldifferentiation2.5Steroid hormonebiosynthesiscountRetrogradeendocannabinoid signaling10Bile secretion2030Nitrogenmetabolism40MorphineadlictionInflammatorymediator regulation of TRPchannelsLeishmaniasis-Thyroid cancerNicotineaddictionAntifolateresistanceFig.6Enrichment results of KEGG pathway10图6KEGG通路富集结果2013123.5“成分-靶点-通路”网络构建和Q-Marker预测分析将9个搜集到靶点的白术Q-Marker候选成分、19个degree值大于2 倍中位数 16.17 的关键靶点和20条通路导入Cytoscape3.9.1软件中构建“白术挥发油Q-Marker候选“成分-靶点-通路”网络，详见图7（图中三角形节点代表活性成分,菱形节点代表关键靶点,方形节点代表通路，节点大小代degree值大小）。以化合物、靶点蛋白、信号通路之间的连接度degree值为参考（图中degree值越大，节点越大,颜色越深）,发现棕榈酸、-桉叶醇、-榄香烯、榄香醇天然产物研究与开发和苍术酮等成分的连接度相对较高，提示这些成分可能是白术挥发油中发挥药效的主要活性物质;靶点MAPK1、PT G S2、R X R A、CY P19A 1、PPA R A、NO S3、TNF、CY P2 C19 的连接度较高,可能是白术挥发油发挥作用的关键靶点；代谢途径（Metabolicpathways）、化学致癌-受体激活（chemical carcinogenesis-receptoractivation）、PPA R 信号通路（PPAR signaling path-way）、血清素能突触（serotonergic synapse）、雌激素信号通路（estrogen signaling pathway）的连接度较大,提示这些通路可能是白术挥发油发挥作用的关键信号通路。Vol.35MotabolicNouroetiveplthwaysnd-recoptorntoractonPPARsignalRerogradescahnattsianaingNROB2ishmomaPGREauogenwgnalingDthwayChoinerpkBZB25NOS3PPARAPPARGsyapseNCOAMAPKTIL10rotonegSteroild hormoneDICosymthesTh17celdiferentiasoinoguhetrecnpiorBZ8CalciuBZ9pathwayHonaingAadhidonicThyroidcanceridAntfolatemetabolianrusistanonPTENPTOS2RXRASLCBAGRMSFABPYESR1CYP2C19TNFCYP1RASABCB1图7“成分-靶点-通路”网络图Fig.7 Composition-target-pathway network diagram注：BZ1：缬草醇；BZ2:（-）-异喇叭烯；BZ3:-榄香烯；BZA:-桉叶醇；BZ5：棕榈酸；BZ6:反式戊烯醇醋酸酯；BZ7：榄香醇；BZ8：苍术酮；BZ9：4-环己基苯乙酮。Note:BZ1:Valerenol;BZ2:(-)-Isotrumpetene;BZ3:-Elemene;BZ4:-Eucalyptol;BZ5:Palmitic acid;BZ6:trans-Valerenylacetate;BZ7:Elemental alcohol;BZ8:Atractylone;BZ9:4-Cyclohexylacetophenone.4讨论与结论本研究通过GC-MS的方法建立了11批次不同产地白术挥发油类成分的指纹图谱,并指认了10 个共有峰,共有成分的总含量均占挥发油总量的8 5%以上。共有成分PCA分析结果显示浙江产白术挥发油与安徽产白术挥发油有明显区分,浙江白术分布更集中；聚类分析显示当分类距离为15时，4批浙江白术可被聚为一类,这表明与安徽白术相比,浙江白术挥发油类化学成分含量更为相似稳定，而四川、云南、江西产区白术由于批次单一，暂无法确定其整体分布以及与浙江、安徽白术的差异。此外,本研究以11批白术饮片挥发油共有成分为出发点,结合网络药理学的分析方法,构建白术挥发油“成分-靶点-通路”作用网络，预测白术挥发油的主要作用通路及靶点。网络药理学结果显示棕榈酸-桉叶醇、-榄香烯、榄香醇和苍术酮可能是白术挥发油的潜在质量标志物。白术挥发油具有调节胃肠道、修复胃黏膜、抑菌、保肝等作用18.9。现代研究表明棕榈酸、-桉叶醇、-榄香烯、榄香醇和苍术酮等成分具有抗菌 17.18 、镇静镇痛 19、调节胃肠运动 2 0.2 1、降压 2 、保肝抗炎 2 3、抗病毒 2 4、抗感染 2 5 等药理作用,与白术挥发油具有相似之处,故这些成分可能是白术挥发油发挥药理作用的主要物质基础,与中药的功能属性密切相关。“成分-靶点-Vol.35通路”网络构建结果显示，MAPK1、PT G S2、R X R A、CYP19A1、PPA R A、NO S3、T NF、CY P2 C19 和代谢途径、化学致癌-受体激活、PPAR信号通路、血清素能突触和雌激素信号通路则可能是白术挥发油发挥作用的关键靶点和信号通路。综上所述,我们的研究发现不同产地的白术饮片挥发油类成分的组成以苍术酮、（-)-异喇叭烯、-榄香烯、-桉叶醇、棕榈酸、trans-valerenyl acetate、榄香醇、4-cyclohexylacetophenone等为主,但含量有所不同。GC-MS 指纹图谱结合网络药理学确定棕榈酸、-桉叶醇、-榄香烯、榄香醇和苍术酮可能是白术挥发油的潜在质量标志物。本研究进一步完善了白术的质量控制体系，为白术的质量控制提供理论依据，也为其他类似中药材的质量标志物研究提供了理论参考。参考文献1 Liu CX,Chen SL,Xiao XH,et al.A new concept on qualitymarker of Chinese materia medica:quality control for Chinesemedicinal productsJ.Chin Tradit Herb Drugs（中草药）,2016,47:1443-1457.2Zhang TJ,Bai G,Liu CX.The concept,core theory and re-search methods of Chinese medicine quality markers J.Ac-ta Pharm Sin(药学学报）,2 0 19,54:18 7-196.3Xu Q,Yang YJ.Quality marker(Q-marker)proposal:a newperspective on Chinese medicine C/Proceedings of the5th National Chinese Medicine Commodity Conference of theChineseCommodity Society（中国商品学会第五届全国中药商品学术大会论文集），2 0 17：94-96.4 Xu S,Xu WF,Jiang WQ,et al.Research progress on qualityevaluation of Atractylodis Macrocephalae Rhizoma J.NorthwestPharm J（西北药学杂志）,2 0 2 2,37:152-155.5Wei HT,Fang QH,Zhang BS,et al.Research progress on theeffective ingredients and pharmacological effects of Atrac-tylodes macrocephala J.Guangdong J Anim Vet Sci(广东畜牧兽医科技）,2 0 2 2,47：40-44.6Yang DY,Yu H,Wu XY,et al.Research progress on chemi-cal constituents and their biological activities of AtractylodisMacrocephalae Rhizoma J.Chin Arch Tradit Chin Med(中华中医药学刊)：1-18 2 0 2 3-0 4-11.http:/ Y,Wei MX,Wu YY,et al.Research progress on ex-traction separation,chemical constitution and pharmacologi-cal activities of polysaccharide extracted from AtractylodesmacrocephalaJ.Chin Tradit Herb Drugs（中草药）,2 0 2 1,吴姣姣等：基于GC-MS指纹图谱和网络药理学的白术挥发性成分质量标志物预测分析究与开发）,2 0 2 3,35:2 0 8-2 2 0.16 Wang H,Liang YQ,Yu XT,et al.Predictive analysis of qual-ity markers of Fructus Auranti Immaturus based on finger-print and network pharmacology J.Chin J Exp Tradit MedForm（中国实验方剂学杂志），2 0 2 2,2 8:198-2 0 8.17 Jing FT,Feng S,Wang J,et al.Predictive