芸香叶片和花瓣释放挥发性有机物成分及其变化规律

2018-06-22 06:54蒋冬月李永红
浙江农林大学学报 2018年3期
关键词:芸香萜烯酮类

蒋冬月,李永红,沈 鑫

(1.浙江省林业科学研究院,浙江 杭州 310023;2.深圳职业技术学院,广东 深圳518055)

芸香Ruta graveolens又名臭草、香草、百应草、小叶香,是芸香科Rutaceae芸香属Ruta多年生草本植物。植株高达1 m,叶羽状复叶,灰绿或带蓝绿色;花金黄色,花期3-6月;原产地中海沿岸,在中国均有栽培[1]。芸香植株各部分均有浓烈的特殊气味,具有杀虫抑菌效果[2-4];全株均可入药,有清热解毒、凉血散瘀、利尿、消肿、止咳平喘等功效[5-6]。目前,对芸香科的研究主要集中在生理特性[7-8]、分子生物学[9-10]、 精油化学成分[3-4,11-13]及抑菌效果[14-15]等方面, 并未见有关芸香花瓣及叶片挥发性有机物成分随时间变化的相关报道。现有的研究探究芸香的有机物成分主要采用精油作为试验材料[16-18],精油的提取常常会受提取温度和方法的干扰,造成某些挥发性物质受到破坏,不能完全反映芸香释放的挥发性有机物的成分。因此,本研究采用顶空固相微萃取技术吸附采集不同时期芸香叶片释放的挥发性有机成分,结合气质联用仪检测并分析不同季节及1 d中不同时间释放的挥发性有机物成分及相对含量,揭示其叶片挥发性有机物的季节变化规律和日变化规律;并用相同的方法探究叶片和花瓣释放的挥发性有机物成分的差异,为科学利用芸香营造生态型、保健型、芳香园林景观提供理论依据。

1 材料和方法

1.1 试验材料

芸香苗龄为5~6年生,株高70~75 cm,长势繁茂,无病虫害,栽植于深圳职业技术学院西校区芳香植物园内。

1.2 试验方法

1.2.1 植物挥发性有机物的采集和测定 分别于3月、6月、9月、12月的16日上午9:00-10:00在同一位置的植物上摘取向阳侧上、中、下部位的健康无损伤的叶片10片,分别剪碎混匀,称取0.5 g放置5 mL萃取瓶中密封,静置30 min,环境温度为(22.0±3.0)℃。将型号为DVB-CAR-PDMS 100 μm的SPME纤维头(美国Supelco公司)通过聚四氟乙烯瓶垫插入到萃取瓶中,置于样品正上方0.5 cm左右,顶空萃取40 min,然后将纤维头插入6890N/5975气相色谱-质谱联用仪(美国Agilent公司)的气相色谱进样口,解吸10 min。每次收集设置3次重复,吸附空萃取瓶中的气体作为空白对照。花瓣挥发性有机物的测定于3月16日进行,9:00-10:00采摘盛开的花5朵;挥发性有机物日动态变化的测定于6月16日进行,8:00-20:00隔2 h采样1次。色谱条件:HP-5MS弹性石英毛细管柱色谱柱,长为30 m,内径为0.25 mm,液膜厚为0.25 μm,载气为高纯氦气,不分流进样,恒流流速为1.0 mL·min-1,进样口温度为230℃,接口温度为280℃。季节变化和花瓣挥发物测定初始温度为50℃,保持4 min,以6℃·min-1升至150℃,保持2 min,然后以7℃·min-1升至250℃,保持8 min。日动态变化规律测定初始温度为50℃,保持4 min,以4℃·min-1升至150℃,保持2 min,然后以8℃·min-1升至250℃,保持6 min。质谱条件:电子轰击(EI)离子源,电子能量为70 eV,离子阱温度为230℃,四级杆温度为150℃,原子质量扫描范围为30~500。

1.2.2 数据处理 挥发性有机物成分经气相色谱分离,各色谱峰、质谱峰分别利用气相色谱-质谱仪中MSD Productivity ChemStation和NIST Mass Spectral Database 2008谱库进行分析鉴定。各成分在样品气体中的浓度采用峰面积归一法进行计算。相对含量=(该物质峰面积/样品所有气体峰面积之和)×100%。

2 结果

2.1 芸香叶片释放挥发性有机物成分和其相对含量的季节动态变化

芸香叶片1 a中不同季节释放的挥发性有机物中,共鉴定出29种化合物(图1,表1),其中酯类11种,酮类7种,烷类4种,萜烯类,酸类和醇类各2种,以及少量的酚类和其他化合物;主要成分为乙酸壬酯、4-(3,4-亚甲基二氧基苯基)-2-丁酮、2-壬酮和2-十一(烷)酮等。1 a中检测到叶片释放的挥发性有机物种类数量呈先增多后减少的趋势;6月最多,为18种,9月最少,仅为10种(图2)。3月挥发性有机物成分检测到酮类、酯类和烷类3类化合物,以酮类最多,其次为酯类;其中2-壬酮的相对含量最高,为48.98%。6月检测到酯类、酮类、萜烯类、酸类、烷类、酚类和醇类7类化合物,以酯类最多,7种。9月检测到的挥发性有机物种类仅为6月的55.56%,酮类和酯类占50.00%,其中2-十一(烷)酮的相对含量最高,为34.13%。12月检测到的挥发物中以酮类为主,占当月总种类的54.55%。1 a中每个季节都可以检测到乙酸壬酯、2-壬酮和2-十一(烷)酮。

图1 芸香叶片不同季节释放挥发性有机物的总离子流图Figure 1 TIC of VOCs released from leaves of Ruta graveolens in different seasons

1 a中不同季节,芸香叶片释放的酯类物质的相对含量呈先上升后下降的趋势,在6月最高,为44.78%,12月最低(图2)。烷类和酮类化合物均呈先下降后上升的趋势,6月其相对含量最低,12月达最高。萜烯类、酸类和醇类化合物只在6月和9月检测到,且6月检测到的萜烯类和酸类化合物相对含量均高于9月;酚类物质仅在6月检测到。叶片释放的主要成分季节变化趋势如下:2-壬酮和2-十一(烷)酮的相对含量随季节呈高-低-高的趋势,而乙酸壬酯则呈现低—高—低的变化趋势。

2.2 芸香叶片释放挥发性有机物的日动态变化

芸香叶片1 d中不同时段释放的挥发性有机物共鉴定出41种 (表2)。其中酯类化合物10种,酮类7种,酸类和醛类各6种,烷类和醇类各4种,以及少量的萜烯类、炔类等化合物;主要成分为乙酸仲辛酯,2-壬酮,2-十一(烷)酮,十六烷酸,油酸等。8:00-20:00叶片释放的挥发性有机物种类呈先增多后减少的趋势(图3),12:00最多,为22种,20:00最少,为14种。8:00检测到酮类、酯类、酸类和醛类等4类化合物,以酮类最多;但挥发物的相对含量以酸类最高,其中十六烷酸和油酸,分别为39.90%和19.31%。10:00挥发物种类较8:00增加了6种,出现了烷类、醇类和炔类;酯类数量最多,其次为酮类,2类物质的相对含量较8:00大幅上升,酮类占此时挥发物含量的58.11%,其中2-壬酮占35.10%;而酸类化合物的相对含量急剧下降,仅为3.31%。12:00-14:00挥发物种类较10:00变化不大,14:00出现了萜烯类,但相对含量较低。16:00-20:00种类逐渐减少,主要成分和相对含量均以酮类和酯类化合物为主;18:00酸类相对含量在1 d中最低,而此时萜烯类和醛类含量达最高。1 d中,每个时段均可检测到(Z)-乙酸-3-己烯-1-醇酯,乙酸壬酯,2-壬酮,2-癸酮,2-十一(烷)酮,2-十三烷酮和十六烷酸7种物质。

图2 1 a中不同季节芸香叶片释放的各类挥发物数量和相对含量的变化Figure 2 Component numbers and relative contents of classified VOCs from leaves in different seasons of a year

表1 芸香叶片释放的挥发性有机物成分及其相对含量的季节变化Table 1 Seasonal variations of components and relative contents of VOCs released from leaves of Ruta graveolens

表2 芸香叶片1 d中不同时间释放的挥发性有机物成分及其相对含量Table 2 Components and relative contents of VOCs released from leaves of Ruta graveolens in one day

1 d不同时段,芸香叶片释放的酯类化合物的相对含量呈先上升后下降的趋势,10:00达最大;酮类呈先上升后略微下降再上升的趋势,20:00达最大,为65.49%;酸类呈先下降后上升再下降的趋势,8:00最大,为59.80%。烷类和醛类相对含量变化不明显;萜烯类、炔类和醇类的相对含量较少,仅在少数几个时间段出现(表3)。其中,主要成分日变化趋势如下:乙酸仲辛酯和2-十一(烷)酮的相对含量呈现低-高的变化趋势,十六烷酸则呈现高—低的变化趋势,2-壬酮的相对含量呈现低—高—低—高的变化趋势,油酸呈现出高—低—高—低的变化趋势。

2.3 芸香花瓣释放的挥发性有机成分及其相对含量

通过对芸香花瓣的SPME/GC/MS总离子流图分析(图4),扣除本底空气中的杂质,3月花期时,花瓣释放的挥发性有机物中检测出8种挥发性成分,其中酮类化合物最多,有5种,其相对含量占花瓣总挥发物的82.24%;其次为萜烯类化合物,其相对含量为11.20%(表4)。花瓣中检测到的相对含量最多的挥发性成分是2-十一(烷)酮,为66.10%,其次为2-壬酮和壬烯。

图3 1 d中不同时段芸香叶片释放的各类挥发物数量的变化Figure 3 Component numbers of classified VOCs from leaves in different times of a day

3 讨论

目前,学者对芸香或芸香科其他植物挥发性有机物的研究多采用水蒸气蒸馏法和溶剂提取法[3-5,11]。蒸馏法简单易行,适于难溶或不溶于水的成分的提取;但属于非活体提取,长时间高温提取会使某些挥发性物质受到破坏;溶剂提取法是根据相似相溶原理,选择适合的溶剂非常关键,但溶剂中的微量杂质沉淀在提取物中,影响提取物的纯度[19]。通过对比各种提取方法,本研究采用了顶空-固相微萃取采集法,此方法所需样品量少,无需溶剂,操作简化,集采样、萃取、浓集、进样于一体,非常适合芳香植物挥发性成分的研究,已在刺山柑Capparis spinosa,茶Camellia sinensis,黄兰 Michelia chmpaka,麦秆菊 Helichrysum bracteatum等多种植物挥发物的研究中得到应用[20-23]。

表3 芸香叶片1 d中不同时段挥发性有机物成分及其相对含量分类统计Table 3 Classified statistic about components and relative contents of VOCs released from leaves in different time of a day

表4 芸香花瓣释放的挥发性有机物成分及其相对含量Table 4 Components and relative contents of VOCs released from flowers of Ruta graveolens

植物释放挥发性有机物的成分和相对含量在同种之间存在差异,与植物的品种、自身结构特点、生理状态和株龄等有密切关系[24-25]。 如 HADDOUCHI 等[4]研究得出阿尔及利亚生长的芸香含有10种挥发性有机物,唐祖年等[15]从芸香挥发油中鉴定出21种成分,徐汉虹等[2]从芸香精油中鉴定出39种成分。本研究在芸香叶片中共鉴定出41种,花瓣中检测到8种,虽然与前人鉴定的结果不尽相同,但通过对比发现,芸香挥发性有机物的主要成分均为2-十一(烷)酮和2-壬酮。主要原因是植物的自身特性决定了不同植物所释放的有机物成分、含量和释放速率不同。

图4 芸香花瓣挥发性有机物总离子图Figure 4 TIC of VOCs released from flowers of Ruta graveolens

前人[26-31]对侧柏 Platycladus orientalis, 薄荷 Mentha×piperita, 红松 Pinus koraiensis, 桂花 Osmanthus fragrans等植物挥发性有机物释放规律的研究表明,挥发物释放的种类和相对含量与采样时温度、湿度、气孔导度、光照以及蒸腾速率等关系密切。一般而言,随着温度、湿度、光照等的提高,植物释放的有机物也会增多。本研究表明,不同季节和1 d中不同时段芸香释放的挥发性有机物的种类均呈先上升后下降的趋势;其中,6月最多,为18种;1 d中12:00最多,为22种。在测试的4个月中,6月光照时间较长,温度较高,达30℃左右,植物进入生长旺盛阶段,植物自身旺盛的光合作用为有机物合成提供充足的碳源和还原力,体内有机物合成、代谢加快,强烈的蒸腾作用也为有机物的扩散创造了条件,使挥发性有机物的种类增多[32-34]。说明随着温度、光照等的提高,芸香叶片释放的挥发性有机物的种类也随之增多。另外,研究发现芸香花瓣释放的挥发性有机物种类少于叶片,叶片和花瓣挥发性有机物的主要成分均为酮类化合物,叶片释放的物质中以2-壬酮为主,花瓣中以2-十一(烷)酮为主。芸香释放的挥发性有机物种类多少、相对含量的高低与外界环境以及植物自身条件的密切相关,两者之间关系有待进一步深入研究。

植物释放的挥发性有机物具有一定的药用保健作用。前人[35-36]研究已证实,针叶树释放的萜烯类化合物使人情绪趋于放松状态,使人感觉清新、舒爽和愉悦;柠檬油、雪松醇、龙脑和1,8-桉叶油能够降低心率和血压;萜烯类和生物碱类物质具有一定的抗菌和消炎作用[3-6,16]。芸香叶片和花瓣释放的挥发性有机物的主要成分均具有一定的药用价值,可以起到杀虫抑菌效果,在某种程度上对调节人体情绪具有一定的作用。如:2-壬酮具有果香、甜香、皂香及椰子、奶油的气味,2-十一(烷)酮是芸香特有的香气,低浓度时有类似桃子的香气;这些有机物的气味或许可以使人趋于放松状态,从而间接起到调节人体情绪的作用。因此,在植物景观设计中,可以适当栽植芸香,营造保健型的生态植物景观。

4 参考文献

[1]曾建飞.中国植物志:第43卷第2分册[M].北京:科学出版社,1997.

[2]徐汉虹,赵善欢,周俊,等.芸香精油的化学成分和杀虫活性初探[J].天然产物研究与开发,1994,6(4):56-61.XU Hanhong,ZHAO Shanhuan,ZHOU Jun,et al.Preliminary studies on insecticidal activity of the rue oil and analysis of its components [J].Nat Product Res Dev,1994,6(4):56 - 61.

[3]REDDY D N,ALRAJAB A J.Chemical composition,antibacterial and antifungal activities of Ruta graveolens L.volatile oils [J].Cog Chem,2016,2(1):1 - 11.

[4]HADDOUCHI F,CHAOUCHE T M,ZAOUALI Y,et al.Chemical composition and antimicrobial activity of the essential oils from four Ruta species growing in Algeria [J].Food Chem,2013,141(1):253 - 258.

[5]顾关云,蒋昱.芸香的化学成分与药理活性[J].国外医药:植物药分册,2003,18(2):47-50.GU Guanyun,JIANG Yu.Rue chemical constituents and pharmacological activities [J].World Phytom,2003,18(2):47-50.

[6]袁海梅,邱露,谢贞建,等.花椒属植物生物碱类成分及其药理活性研究进展[J].中国中药杂志,2015,40(23): 4573 - 4584.YUAN Haimei,QIU Lu,XIE Zhenjian,et al.Research progress on alkaloids constituents from Zanthoxylum and their pharmacological activities [J].China J Chin Mater Med,2015,40(23):4573 - 4584.

[7]刘文哲,胡正海.中国芸香科植物叶分泌囊比较解剖学研究[J].植物分类学报,1998,36(2):119-127.LIU Wenzhe,HU Zhenghai.Comparative anatomy of secretory cavities in leaves of the Rutaceae in China [J].Acta Phytotaxon Sin,1998,36(2):119 - 127.

[8]刘文哲,胡正海.芸香科植物茎分泌囊的研究[J].西北植物学报,1999,19(3):456-460.LIU Wenzhe,HU Zhenghai.Studies on the secretory cavity of stems in Rutaceae [J].Acta Bot Boreal-Occident Sin,1999,19(3):456 - 460.

[9]罗焜,陈士林,陈科力,等.基于芸香科的植物通用DNA条形码研究[J].中国科学C辑:生命科学,2010, 40(4):342 - 358.LUO Kun,CHEN Shilin,CHEN Keli,et al.Assessment of candidate plant DNA barcodes using the Rutaceae family[J].Sci China Life Sci,2010,40(4):342 - 358.

[10]柴素芬,陈兆贵,洪彩霞.芸香科药用植物ISSR-PCR反应体系的建立及优化[J].安徽农业科学,2008,36(33): 14433 - 14435.CHAI Sufen,CHEN Zhaogui,HONG Caixia.Construction and optimization of ISSR reaction system of medicinal Rutaceae plants [J].J Anhui Agric Sci,2008,36(33):14433 - 14435.

[11]王兆玉,郑家欢,黎恩立,等.3种芸香科植物果皮挥发油成分GC-MS分析与比较[J].中药材,2016,39(5):1071 - 1074.WANG Zhaoyu,ZHENG Jiahuan,LI Enli,et al.GC-MS analysis and comparison of VOCs of 3 Rutaceae plants peel[J].J Chin Med Mater,2016,39(5):1071 - 1074.

[12]刘向前,魏圣淑,白完淑,等.中国产和韩国产芸香科植物的成分比较研究[J].中草药,2003,34(7):586-589.LIU Xiangqian,WEI Shengshu,BAI Wanshu,et al.Studies on constituents of Rutaceae plants growing in China comparied with those in Korea [J].Chin Trad Herb Drug,2003,34(7):586 - 589.

[13]杨序成,侯娜.顶空固相微萃取-气质联用法测定日本花椒挥发油成分分析[J].贵州科学,2017,35(1):94-96.YANG Xucheng,HOU Na.Composition analysis of volatile oil of Zanthoxylum japonica by SPMEGC-MS [J].Guizhou Sci,2017,35(1):94 - 96.

[14]吴娇,苗辉,邝芷琪,等.芸香科植物杀虫及抑菌活性研究进展[J].中国植保导刊,2015,35(10):18-26.WU Jiao,MIAO Hui,KUANG Zhiqi,et al.Research progress on insecticidal and antifungal activity of Rutaceae plants [J].China Plant Prot,2015,35(10):18 - 26.

[15]唐祖年,杨月,杨扬,等.芸香挥发油GC-MS分析及其生物活性研究[J].中国现代应用药学,2011,28(9): 834 - 838.TANG Zunian,YANG Yue,YANG Yang,et al.Chemical composition and biological activity of the essential oil of Ruta graveolens [J].Chin J Mod Appl Pharm,2011,28(9):834 - 838.

[16]IVANOVA A,MIKHOVA B,NAJDENSKI H,et al.Antimicrobial and cytotoxic activity of Ruta graveolens [J].Fitoterapia,2005,76(3/4):344 - 347.

[17]de FEO V,de SIMONE F,SENATORE F.Potential allelochemicals from the essential oil of Ruta graveolens [J].Phytochemistry,2002,61(5):573-578.

[18]DZHURMANSKI A,ZHEKOVA G,ANGELOVA D.Accumulation dynamic of Ruta graveolens L.essential oil [J].Agric Sci Technol,2011,3(4):343 - 345.

[19]蒋冬月,李永红.植物挥发性有机物的研究进展[J].黑龙江农业科学,2011(11):143-149.JIANG Dongyue,LI Yonghong.Review on volatile organic compounds of plant [J].Heilongjiang Agric Sci,2011(11):143 - 149.

[20]ROMEO V,ZIINO M,GIUFFRIDA D,et al.Flavour profile of capers (Capparis spinosa L.) from the Eolian Archipelago by HS-SPME/GC-MS [J].Food Chem,2007,101(3):1272 - 1278.

[21]WU Yuanshuang,LÜ Shidong,LIAN Ming,et al.Study of characteristic aroma components of baked Wujiatai green tea by HS-SPME/GC-MS combined with principal component analysis [J].CyTA J Food,2016,14(3):1 - 10.

[22]蒋冬月,李永红,夏兵,等.黄兰叶片和花瓣挥发性成分及其抑菌效果[J].东北林业大学学报,2012,40(5): 71 - 74.JIANG Dongyue,LI Yonghong,XIA Bing,et al.Components and antibacterial functions of volatile organic compounds from leaves and flowers of Michelia chmpaka [J].J Northeast For Univ,2012,40(5):71 - 74.

[23]GIULIANI C,LAZZARO L,CALAMASSI R,et al.A volatolomic approach for studying plant variability:the case of selected Helichrysum species (Asteraceae) [J].Phytochemistry,2016,130:128.

[24]BENJAMIN M T,WINER A M.Estimating the ozone-forming potential of urban trees and shrubs [J].Atmos Environ,1998,32(1):53 - 68.

[25]OWEN S M,HARLEY P,GUENTHER A,et al.Light dependency of VOC emissions from selected Mediterranean plant species [J].Atmos Environ,2002,36(19):3147 - 3159.

[26]杨伟伟,王成,郄光发,等.北京西山春季侧柏游憩林内挥发物成分及其日变化规律[J].林业科学研究,2010, 23(3): 462 - 466.YANG Weiwei,WANG Cheng,QIE Guangfa,et al.Compounds and diurnal variation of VOCs of Platycladus orientalis recreation forest in Beijing western hills in spring [J].For Res,2010,23(3):462 - 466.

[27]李娟,王成,彭镇华,等.侧柏春季挥发物浓度日变化规律及其影响因子研究[J].林业科学研究,2011,24(1): 82 - 90.LI Juan,WANG Cheng,PENG Zhenhua,et al.The diurnal variation and influence factors of VOC of Platycladus orientalis in spring [J].For Res,2011,24(1):82 - 90.

[28]GERSHENZON J,McCONKEY M E,CROTEAU R B.Regulation of monoterpene accumulation in leaves of peppermint[J].Plant Physiol,2000,122(1):205 - 213.

[29]SON Y S,KIM K J,JUNG I H,et al.Seasonal variations and emission fluxes of monoterpene emitted from coniferous trees in East Asia:focused on Pinus rigida and Pinus koraiensis [J].J Atmos Chem,2015,72(1):27 - 41.

[30]SAATHOFF H,NAUMANN K H,HLER O M,et al.Temperature dependence of yields of secondary organic aerosols from the ozonolysis of α-pinene and limonene [J].Atmos Chem Physics,2009,9(5):1551 - 1577.

[31]蔡宙霏,陈雅奇,许馨露,等.4个桂花品种开花进程释放VOCs动态变化分析[J].浙江农林大学学报,2017,34(4): 608 - 619.CAI Zhoufei,CHEN Yaqi,XU Xinlu,et al.Changes of volatile organic compounds released during flowering in four Osmanthus fragrans cultivar groups [J].J Zhejiang A&F Univ,2017,34(4):608 - 619.

[32]STAUDT M,LHOUTELLIER L.Monoterpene and sesquiterpene emissions from Quercus coccifera exhibit interacting responses to light and temperature [J].Biogeosci Discuss,2011,8(9):5691 - 5728.

[33]SHARKEY T D,LORETO F.Water stress,temperature,and light effects on the capacity for isoprene emission and photosynthesis of Kudzu leaves [J].Oecologia,1993,95(3):328 - 333.

[34]HARLEY P,GUENTHER A,ZIMMERMAN P.Effects of light,temperature and canopy position on net photosynthesis and isoprene emission from sweetgum (Liquidambar styraciflua) leaves [J].Tree Physiol,1996,16(1/2):25 -32.

[35]高岩.北京市绿化树木挥发性有机物释放动态及其对人体健康的影响[D].北京:北京林业大学,2005.GAO Yan.Releasing Variation and Effects on Human Health of Volatile Organic Compounds from Landscape Trees in Beijing [D].Beijing:Beijing Forestry University,2005.

[36]王艳英,王成,蒋继宏,等.侧柏、香樟枝叶挥发物对人体生理的影响[J].城市环境与城市生态,2010,23(3): 30 - 32, 37.WANG Yanying,WANG Cheng,JIANG Jihong,et al.Effect of VOCs from branch and leaf of Piatycladus orientalis and Cinnamomum camphora on human physiology [J].Urban Environ Urban Ecol,2010,23(3):30 - 32,37.

猜你喜欢
芸香萜烯酮类
广州市环境空气中醛酮类化合物的污染特征研究
割手密萜烯合成酶(TPS)基因家族分析及其在生物胁迫下的表达分析
大气中醛酮类污染现状研究
基于网络药理学的2-(2-苯乙基)色酮类成分潜在作用靶点及作用机制预测
漫步在森林当中为何让人感觉心情舒畅?
书香究竟是什么香
8个杧果品种胶乳挥发性成分及胶乳的灼伤机制初探
书香门第中的“书香”是一种什么香
高脂肪低碳水化合物有助对抗流感
破壁灵芝孢子粉的灵芝清香及其成分来源探讨