蔡 萍，何 丹，陳 林，梁雪娟，唐純玉，張水寒，黃惠勇

（1．湖南省中醫(yī)藥研究院中藥研究所，湖南長沙410013；2．湖南省中醫(yī)藥大學，湖南長沙410208；3．湖南時代陽光藥業(yè)股份有限公司，湖南永州425000；4．湖南省藥食同源功能性食品工程技術研究中心，湖南長沙410208）

2015版《中國藥典》收載的藥用百合為百合科植物卷丹（LiliumlancifoliumThunb．）、百合（Lilium browniiF．E．Brown var．viridulumBaker）或細葉百合（LiliumpumilumDC．）的干燥肉質鱗葉［1］。百合最早收錄于《神農(nóng)本草經(jīng)》，是藥食兩用之品。生百合及蜜百合為臨床常用飲片，煉蜜有“補脾氣，潤肺燥”之功，故蜜炙后百合潤肺止咳的作用增強，多用于肺虛久咳、肺癆咳血等癥，是輔料與藥物的協(xié)同作用所致。

高效液相（HPLC）指紋圖譜強調對樣品組成的系統(tǒng)表征及穩(wěn)定性組分的鑒定及評估，被廣泛用于中藥質量控制研究［2－3］。但單一的指紋圖譜無法明確樣品間差異來源，在精確篩選質量控制的指標性成分方面存在缺陷?；瘜W計量學方法，如相似性分析（SA）、主成分分析（Principal Component Analysis，PCA）和偏最小二乘法－判別分析（PLS－DA）等已經(jīng)被世界衛(wèi)生組織（WHO）［4］及其他部門［5］引入，并接受用于復雜多分量系統(tǒng)的定性分析和質量控制［6－9］。本研究收集23批不同來源的生百合及蜜炙品，建立百合蜜炙前后的HPLC指紋圖譜；利用中藥相似性分析軟件，找到兩類百合的共有峰；應用隨機森林算法對指紋進行分類識別，為百合的質量控制和指標成分的選擇提供依據(jù)，為全面控制百合質量提供參考。

1 實驗部分

1．1 儀器及試劑

Agilent1200高效液相色譜儀，包括二極管陣列檢測器、自動脫氣機、四元泵、自動進樣器，柱溫箱；KQ5200DE型數(shù)控超聲波清洗器（昆山市超聲儀器有限公司）；AL204電子天平（梅特勒－托利多儀器上海有限公司）；布氏漏斗。

乙腈、甲醇（色譜純，德國Merck公司）；H3PO4（分析純，湖南匯虹試劑有限公司）；無水乙醇（分析純，湖南匯虹試劑有限公司）；純凈水（華潤怡寶飲料有限公司）。

1．2 實驗樣本

本研究通過購買和到主產(chǎn)區(qū)采集方式，獲得了23份蜜炙百合與生百合樣本，并經(jīng)湖南省中醫(yī)藥研究所劉浩助理研究員鑒定為百合科植物卷丹（Lilium lancifoliumThunb．）的干燥肉質鱗葉，樣品包括13批蜜炙百合：S1（豪州市貢藥飲片廠）、S2（安徽大別山）、S3（湖南衡陽）、S4（湖南龍山石牌鎮(zhèn)）、S5（湖南龍山冼洛鄉(xiāng)）、S6（安徽毫藥千草國藥股份有限公司）、S7（湖南隆回縣北山鄉(xiāng)）、S8（亳州華云飲片有限公司）、S9（河北祁新中藥顆粒飲片有限公司）、S10（保定商宇藥業(yè)有限公司）、S11（湖南龍山縣茅坪鄉(xiāng)干壩村）、S12（保定商宇藥業(yè)有限公司）、S13（湖南白羊鄉(xiāng)）；以及10批生品百合：S14（豪州市貢藥飲片廠）、S15（湖南衡陽）、S16（湖南邵東縣界嶺）、S17（湖南龍山石牌鎮(zhèn)）、S18（云南灣塘鄉(xiāng)）、S19（湖南龍山冼洛鄉(xiāng)）、S20（亳州華云飲片有限公司）、S21（湖南龍山縣石牌鎮(zhèn)紅橋村）、S22（湖南龍山縣三元鄉(xiāng)青崗村）、S23（湖南省龍山縣桶車鄉(xiāng)）。

1．3 實驗條件

1．3．1 供試品溶液的制備 取百合樣品粉末（過40目標準檢驗篩）約4．0g，精密稱定，加入80%甲醇溶液50mL，靜置20min，超聲30min，過濾，濾液經(jīng)0．22μm微孔濾膜濾過，制得供試品溶液。

1．3．2 色譜條件 色譜柱：Eclipse Plus C18柱（250×4．6mm，5μm）；流動相：乙腈（A）－0．3%H3PO4溶液（B），梯度洗脫：0～5min，5%～8%A；5～7min，8%A；7～18min，8%～16%A；18～20min，16%A；20～25min，16%～18%A；25～28min，18%～20%A；28～38min，20%A；38～43min，20%～24%A；43～50min，24%A。流速：1．0mL／min；柱溫：25℃；進樣體積：20μL。檢測波長：250nm。

1．3．3 質譜條件 電噴霧電離（ESI）；檢測方式：負離子模式；干燥氣溫度325℃；干燥氣體積流量11．0 L／min；鞘氣溫度350℃；毛細管電壓4．0V。負離子模式獲得電噴霧電離質譜圖，全掃描質譜范圍m／z100～1 000。

1．4 方法學考察

1．4．1 精密度實驗 取湖南龍山冼洛鄉(xiāng)蜜百合供試品溶液（樣品5），按優(yōu)化色譜條件連續(xù)進樣6次，記錄指紋圖譜。以6號峰為參照峰計算各峰相對保留時間和相對峰面積的精密度，其相對標準偏差（RSD）值均小于3，符合指紋圖譜技術要求，表明儀器等整個檢測系統(tǒng)精密度高。

1．4．2 重復性實驗 平行制備6份供試品溶液，按“1．3．2”色譜條件進行分析。計算各峰相對保留時間和相對峰面積考察方法重復性，其RSD均小于3%，說明方法重現(xiàn)性良好。

1．4．3 穩(wěn)定性實驗 取同一供試品溶液（樣品5），按優(yōu)化色譜條件，分別在0、6、12、24、48h檢測指紋圖譜。計算各峰相對保留時間和相對峰面積，其RSD值均小于3%，表明樣品在48h內穩(wěn)定性良好。

1．5 統(tǒng)計分析

經(jīng)《中藥色譜指紋圖譜相似度評價（2004藥典版）》軟件分析百合HPLC指紋圖譜數(shù)據(jù)，確定共有峰，生成對照圖譜并進行相似度評價。以共有峰面積為變量，將各批次樣品數(shù)據(jù)導入Matlab軟件，采用隨機森林算法來分析研究不同來源蜜炙百合和生百合的差異，篩選導致樣品間差異的化學成分。

2 結果與討論

2．1 提取條件選擇

考察了80%乙醇、50%甲醇、80%甲醇、100%甲醇等不同提取溶劑對提取效率的影響。結果顯示，以80%甲醇作為提取溶劑獲得的成分最豐富，含量相對較高。

2．2 色譜條件優(yōu)化

使用二極管陣列檢測器在200～400nm波長范圍對樣品進行掃描，分析比較各波長下色譜圖。結果顯示，250nm下檢測到各成分相對較好，分離度良好，峰數(shù)目最多，基線平穩(wěn)，特征峰明顯且峰形較好，能最大程度的反映樣品信息。其次對水相（0．3%甲酸水溶液和0．3%H3PO4水溶液）進行考察，結果表明流動相pH值對百合酚酸糖苷類類化合物色譜峰的分離有一定影響，0．3%H3PO4水溶液作為水相可獲得更好的分離效果。

2．3 百合HPLC指紋圖譜

將23批生百合和蜜百合按優(yōu)化色譜條件分別進樣，得到各樣品的HPLC色譜圖。采用《中藥色譜指紋圖譜相似度評價（2004藥典版）》軟件分析，選擇分離度好且出峰數(shù)多的18號樣品作為參照譜圖，時間寬度0．5，中位數(shù)為相關系數(shù)對樣品進行自動匹配，得到生百合和蜜百合不同來源指紋圖譜分析，見圖1。

圖1 23批百合的HPLC指紋圖譜Fig．1 HPLC fingerprint of 23batches of Lilium

2．4 共有峰確認

23批蜜百合和生百合共提取到13個共有峰，生成對照指紋圖譜（圖2）。按照優(yōu)化質譜條件，采用負離子模式分析，對13個共有峰進行分析鑒定，根據(jù)其質譜信息和裂解規(guī)律，參考文獻報道共解析得到6個酚酸糖苷類成分，分別為王百合苷C（3號峰）、王百合苷A（5號峰）王百合苷D（6號峰）、王百合苷F（7號峰）王百合苷B（10號峰）和4－乙酰王百合苷D（11號峰）。質譜結果見圖3。

圖2 百合對照指紋圖譜及共有峰Fig．2 Lilium control fingerprint and common peaks3．Regaloside C；5．Regaloside A；6．Regaloside D；7．Regaloside F；10．Regaloside B，11．4－Acetyl derivative of regalosid e D．

化合物3，tR＝16．278min，m／z：415．1273［M－H］－，二級 掃 描 裂 解 產(chǎn) 生 m／z：253［M－H－C6H10O5］－，179［M－H－C9H16O7］－，161［M－H－C9H16O7－H2O］－，135［M－H－C9H16O7－CO2］－（圖3），分子式為 C18H24O11，相對分子質量為415．1273，結合文獻報道［10］鑒別化合物3為王百合苷C。

化合物5和化合物6的相對保留時間分別為：21．153min與24．198min，它們具有相同m／z：399［M－H］－，二級掃描裂解產(chǎn)生相同裂解碎片 m／z：237［M－H－C6H10O5］－、163［M－H－C9H16O7］－、145［M－H－C9H16O7－H2O］－、119［M－H－C9H16O7－CO2］－（圖3），分子式為 C18H24O10，相對分子質量分別為：399．1317與399．1316，結合出峰先后及文獻報道［11］，鑒別化合物5和化合物6分別為王百合苷A、王百合苷D。

化合物7，tR＝26．280min，m／z：429．1424［M－H］－，二級掃描裂解產(chǎn)生 m／z：429［M－H］－、267［M－H－C6H10O5］－、193［M－H－C9H16O7］－、175［M－H－C9H16O7－H2O］－、149［M－H－C9H16O7－CO2］－（圖3），分子式為C19H26O11，相對分子質量為429．1424，結合文獻報道［12］，鑒別化合物7為王百合苷F。

化合物10和化合物11的相對保留時間分別為：36．929min與37．614min，它們具有相同的m／z：441［M－H］－，二 級 掃 描 裂 解 產(chǎn) 生 m／z：381［M－H－C2H2O－H2O］－、163［M－H－C11H18O8］－、145［M－H－C11H18O8－H2O］－、119［M－H－C11H18O8－CO2］－（圖3），分子式為C20H26O11，結合出峰先后及文獻報道［13］，鑒別化合物10和化合物11分別為王百合苷B、4－乙酰王百合苷D。

2．5 無監(jiān)督模式識別的分析結果

圖3 負離子模式下百合中5個成分的二級質譜圖Fig．3 Mass spectra of 5species in Lilium in negative ion mode

利用23批樣品中的13個共有峰及其峰面積，構建出23×13的數(shù)據(jù)矩陣，將得到的數(shù)據(jù)矩陣輸入隨機森林算法進行分類識別，用于判別生品和蜜炙百合間的差異。在構建分類識別模型的過程中，可得到樣品間的相似性矩陣，并進一步利用多維特征尺度變換（MSD），將其映射到特征空間從中找到樣品間的相互關系（相似或相近程度）（圖4）?？梢娚泛兔壑税俸祥g存在較大的差異，并且在構建分類模型的過程中，可得到各共有峰對分類的貢獻度，以貢獻度的高低來判別各共有峰的重要性，見圖5。從圖5中可見，峰6（王百合苷D）和峰10（王百合苷B）具有最高的貢獻度，為百合品類的鑒別提供了更加確切的依據(jù)。

圖4 生百合與蜜百合分類圖Fig．4 Classification map of candied Lilium and raw Lilium

圖5 各變量的分類重要度Fig．5 Classification importance of each variable

3 結論

本文利用高效液相色譜－飛行時間－串聯(lián)質譜聯(lián)用技術和模式識別方法，對不同來源的生百合與蜜百合樣本進行分析，能有效區(qū)分生品與蜜炙的樣品。本研究是對現(xiàn)有方法的有效補充，隨機森林算法更加適合多類數(shù)據(jù)的分類分析，能為百合的質量評估提供新的方法和依據(jù)。