• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    結合非靶向代謝組學和指紋圖譜對通體結香技術產沉香與野生沉香化學成分的比較分析

    2022-09-20 09:16:46侯文成劉洋洋魏建和
    中草藥 2022年18期
    關鍵詞:結香通體乙基

    郁 萌,侯文成,馮 劍,陳 蘭,楊 云,劉洋洋*,魏建和*

    結合非靶向代謝組學和指紋圖譜對通體結香技術產沉香與野生沉香化學成分的比較分析

    郁 萌1, 2,侯文成1,馮 劍1,陳 蘭1,楊 云1,劉洋洋1, 2*,魏建和1, 2*

    1. 中國醫(yī)學科學院 北京協和醫(yī)學院 藥用植物研究所海南分所,海南省南藥資源保護與開發(fā)重點實驗室,國家中醫(yī)藥管理局沉香可持續(xù)利用重點研究室,海南 ???570311 2. 中國醫(yī)學科學院 北京協和醫(yī)學院 藥用植物研究所,中草藥物質基礎與資源利用教育部重點實驗室,瀕危藥材繁育國家工程實驗室,北京 100193

    系統全面分析比較通體結香技術產沉香(簡稱通體香)與野生沉香化學成分,為通體結香技術進一步的應用和推廣提供科學可靠的依據。采用超高效液相色譜串聯四級桿飛行時間質譜(UPLC-Q/TOF-MS)和氣相色譜-質譜(GC-MS)結合非靶向代謝組學和指紋圖譜,對通體香與野生沉香的乙醇提取物和揮發(fā)性成分進行分析。16批通體香與16批野生沉香的共有成分相同,但通體香指紋圖譜相似度較高,且4個倍半萜、3個小分子芳香物質和8個2-(2-苯乙基)色酮等共17個共有成分在通體香中含量較高。通體香和野生沉香主要化學成分組成相同,但不同批次間通體香的質量更穩(wěn)定。

    沉香;通體結香技術;2-(2-苯乙基)色酮;非靶向代謝組學;指紋圖譜;2,4-二叔丁基苯酚;6,7-二甲基-1,2,3,5,8,8a-六氫萘;去氫蜂斗菜酮

    沉香為名貴香藥,《中國藥典》2020年版規(guī)定其為瑞香科植物白木香(Lour.) Gilg含有樹脂的木材[1]。自然情況下,白木香樹須受到外界傷害后,在其生長過程中緩慢形成沉香[2]。為了減少對野生沉香資源的依賴,我國從20世紀90年代就開始大力發(fā)展白木香種植,并探索人工結香方法。其中,魏建和等[3]發(fā)明的通體結香技術已成為目前我國乃至全世界范圍內主要沉香產區(qū)應用最廣泛的沉香結香技術。李浩洋等[4]對10批通體結香技術所產沉香(簡稱通體香)的質量進行評價,認為通體香的質量符合《中國藥典》的規(guī)定,可代替野生沉香入藥。相關研究證明通體香對實驗動物的止嘔[5]、保肝[6]、抗心肌缺血[7]、抗胃潰瘍[8]和抗結腸炎[9]等作用與野生沉香相當或優(yōu)于野生沉香。Li等[10]采用UPLC-Q/TOF-MS結合多元變量統計分析初步發(fā)現通體香和野生沉香中6,7-二甲氧基- 2-(2-苯乙基)色酮等14個2-(2-苯乙基)色酮類成分在含量上存在一定差異。而劉洋洋[11]研究發(fā)現通體香乙醇提取物中主要成分為2-(2-苯乙基)色酮類化合物(圖1),其揮發(fā)性成分主要為倍半萜類化合物和小分子芳香物質。已有研究尚未從全化學組成角度分析比較通體香與野生沉香化學成分的差異。

    為了系統全面地分析和評價通體香藥材質量,本研究收集我國沉香2大主要產區(qū)海南和廣東的具有代表性的通體香和野生沉香樣品,采用超高效液相色譜串聯四級桿飛行時間質譜(UPLC-Q/TOF- MS)和氣相色譜質譜(GC-MS)分析沉香化學成分,基于中藥指紋圖譜相似度評價方法和非靶向代謝組學方法對通體香與野生沉香的化學成分進行比較分析,從而評價通體香與野生沉香化學成分的組成差異,以期為通體結香技術的推廣以及通體香在醫(yī)藥和大健康產業(yè)中的應用提供依據。

    圖1 2-(2-苯乙基)色酮類化合物結構

    1 材料、儀器與試劑

    1.1 材料

    32批沉香樣品詳細信息見表1。其中16批為通體香樣品,16批為野生沉香樣品,所有樣品經中國醫(yī)學科學院藥用植物研究所海南分所劉洋洋研究員鑒定為瑞香科植物白木香(Lour.) Gilg含有樹脂的木材。

    表1 樣品信息

    Table 1 Information of samples

    序號來源分類序號來源分類 tt1廣東茂名通體香ys1海南五指山野生沉香 tt2廣東化州通體香ys2海南定安野生沉香 tt3廣東茂名通體香ys3海南定安野生沉香 tt4廣東湛江通體香ys4市售(購買地廣東江門)野生沉香 tt5廣東湛江通體香ys5市售(購買地廣東江門)野生沉香 tt6廣東化州通體香ys6市售(購買地廣東江門)野生沉香 tt7廣東化州通體香ys7廣東化州野生沉香 tt8廣東電白通體香ys8廣東化州野生沉香 tt9廣東化州通體香ys9海南定安野生沉香 tt10海南瓊中通體香ys10海南定安野生沉香 tt11海南??谕w香ys11海南定安野生沉香 tt12海南??谕w香ys12海南定安野生沉香 tt13海南定安通體香ys13市售(購買地廣東江門)野生沉香 tt14海南瓊中通體香ys14海南定安野生沉香 tt15海南儋州通體香ys15海南定安野生沉香 tt16海南儋州通體香ys16市售(購買地廣東江門)野生沉香

    1.2 儀器

    ASE 350型加速溶劑萃取儀(賽默飛世爾科技有限公司);7890A-5975C氣相色譜質譜聯用儀(GC-MS,安捷倫科技有限公司);Xevo G2-XS型超高效液相色譜串聯四極桿飛行時間質譜儀(美國沃特世科技有限公司)。

    1.3 試劑

    硅藻土(賽默飛世爾科技有限公司,粒徑大于三號篩),乙腈(賽默飛世爾科技有限公司,質譜純),屈臣氏蒸餾水,其余試劑均為國產分析純。

    2 方法

    2.1 UPLC-Q/TOF-MS分析方法

    2.1.1 樣品溶液制備 精密稱取沉香樣品粉末0.1 g,加50%乙醇10 mL,稱定質量,超聲(40 kHz,600 W)提取30 min,靜置、冷卻,稱定質量,用50%乙醇補足質量,取上清液作為UPLC分析用樣品溶液,待用。從以上32個樣品溶液中各取200 μL混合,制成質控樣品[12]。

    2.1.2 分析條件 采用UPLC-Q/TOF-MS分析32個樣品溶液,每6個樣品間插入1個質控樣品。

    (1)UPLC色譜條件:Waters ACQUITY UPLC BEH C18色譜柱(100 mm×3 mm,1.7 μm),流動相乙腈(A)-水(B)梯度洗脫(0~20 min,20%~50%A;20~30 min,50%~100%A),體積流量0.25 mL/min,柱溫30 ℃,檢測波長254 nm,進樣量2 μL。

    (2)質譜條件:Q/TOF-MS采用電噴霧電離源(ESI)在正離子模式下以MSE模式采集,掃描范圍/100~1500,采用亮氨酸-腦啡肽(/556.277 1 [M+H]+)實時校正。毛細管電離電壓3.0 kV;錐孔電壓40 V;離子源溫度100 ℃;脫溶劑氣溫度350 ℃;脫溶劑N2體積流量600 L/h;錐孔反吹N2體積流量50 L/h。低能量掃描無碰撞能,高能量MSE模式掃描碰撞能35~45 eV,碰撞氣為氬氣。

    2.2 GC-MS分析方法

    2.2.1 樣品溶液制備 精密稱取0.1 g沉香樣品粉末,與適量的硅藻土研磨混勻,裝入34 mL不銹鋼萃取池中,以硅藻土填滿萃取池,置于加速溶劑萃取儀中。加速溶劑萃取儀以水為提取劑,加熱溫度120 ℃,靜態(tài)時間5 min,沖洗體積為萃取池體積的30%,循環(huán)2次。所得提取液用5 mL正己烷萃取,制得GC-MS分析用樣品溶液,待用。從32個樣品溶液中各取200 μL混合,制成質控樣品[12]。

    2.2.2 分析條件 采用GC-MS以液體進樣方式分析32個樣品溶液,記錄TIC圖。每6個樣品間插入1個質控樣品。

    (1)色譜條件:色譜柱為HP-5MS石英毛細管柱(30 m×250 μm),進樣口溫度240 ℃;柱溫50 ℃保持1 min,以15 ℃/min至140 ℃并保持8 min;以1 ℃/min至155 ℃并保持8 min;以10 ℃/min至175 ℃并保持7.5 min;以5 ℃/min至200 ℃并保持9.5 min;最后以20 ℃/min至260 ℃并保持5 min。進樣量1 μL。

    (2)質譜條件:EI離子源,離子阱溫度230 ℃,四級桿溫度150 ℃,接口溫度260 ℃,質量范圍/50~300。

    2.3 共有峰分析

    將16批通體香和16批野生沉香樣品的GC-MS總離子流和UPLC-Q/TOF-MS基峰圖(BPC,base peak chromatogram)分別導入《中藥色譜指紋圖譜相似度評價系統軟件》(2012版),以質控樣品QC指紋圖譜為參照圖譜,時間窗寬度為0.5 s,經多點校正、自動匹配后,分別生成通體香的GC-MS、UPLC-Q/TOF-MS指紋圖譜和野生沉香GC-MS、UPLC-Q/TOF-MS指紋圖譜,并依次提取共有峰。

    2.4 多元變量統計分析

    采用MS-DIAL軟件[13]對GC-MS原始數據進行峰提取、識別、對齊、歸一化、尺度化和缺失值填充等數據預處理,采用Waters Progenesis QI軟件[14]對UPLC-Q/TOF-MS原始數據進行峰提取、識別、對齊、歸一化、尺度化和缺失值填充等數據預處理,并分別輸出每個樣品標準化后得到的保留時間、精確質荷比和峰面積組成的三維數據矩陣,導入至SIMCA 14.0軟件分別對GC-MS和UPLC-Q/TOF- MS數據進行多元變量統計分析。

    3 結果與分析

    3.1 指紋圖譜方法學考察

    3.1.1 精密度試驗 精密稱取質控樣品,分別按“2.1.1”和“2.2.1”項下方法制得樣品溶液,并分別按“2.1.2”和“2.2.2”項下色譜條件連續(xù)進樣測定6 次。UPLC-Q/TOF-MS以2號峰為參照峰、GC-MS以3號峰為參照峰,計算各共有峰相對保留時間和相對峰面積RSD值分別小于0.59%和3.30%、0.079%和2.49%,表明儀器精密度良好。

    3.1.2 重復性試驗 精密稱取質控樣品6份,分別按“2.1.1”和“2.2.1”項下方法制備樣品溶液,并分別按“2.1.2”和“2.2.2”項下色譜條件進樣,UPLC-Q/TOF-MS以2號峰為參照峰、GC-MS以3號峰為參照峰計算各共有峰相對保留時間和相對峰面積RSD值分別小于0.93%和3.51%、0.084%和2.54%,表明該方法重復性良好。

    3.1.3 穩(wěn)定性試驗 精密稱質控樣品,分別按“2.1.1”和“2.2.1”項下方法制得樣品溶液,并分別按“2.1.1”和“2.2.1”項下色譜條件分別于0、2、4、6、12、24 h進樣,UPLC-Q/TOF-MS以2號峰為參照峰、GC-MS以3號峰為參照峰計算各共有峰相對保留時間和相對峰面積RSD值分別小于0.76%和3.74%、0.05%和2.33%,表明24 h內樣品溶液在室溫下的穩(wěn)定性良好。

    3.2 指紋圖譜的建立及共有峰分析

    采用UPLC-Q/TOF-MS對16批通體香和16批野生沉香50%乙醇提取物的化學成分指紋圖譜進行分析,見圖2。結果表明16批通體香樣品之間相似度為0.598~0.869,而16批野生沉香樣品間的相似度僅為0.092~0.815,可見16批通體香醇提物化學成分的相似度較16批野生沉香醇提物化學成分相似度大。同時,基于MSE質譜采集模式高能量場得到的碎片離子信息,并參考已有文獻數據[10, 15-16]確認了2-(2-苯乙基)色酮類化合物色酮部分(A環(huán))和苯乙基部分(B環(huán))取代基的類型和數目,進而指認了32個2-(2-苯乙基)色酮類化合物,見表2。

    根據指紋圖譜共有峰分析結果可知,通體香的UPLC-Q/TOF-MS指紋圖譜共有峰(RLC-tt,圖3-A)和野生沉香的UPLC-Q/TOF-MS指紋圖譜共有峰(RLC-ys,圖3-B)均為24個??梢?,通體香和野生沉香50%乙醇提取物的主要化學成分一致,均為2-(2-苯乙基)色酮類成分,且根據分子離子峰的相對豐度比較,各成分在通體香中的含量均高于野生沉香。

    采用GC-MS分析16批通體香和16批野生沉香的揮發(fā)性成分,結果發(fā)現16批通體香樣品之間相似度為0.846~0.995,16批野生沉香樣品間的相似度僅為0.393~0.974,可見16批通體香揮發(fā)性成分的相似度較16批野生沉香揮發(fā)性成分相似度大。

    圖2 通體香(I) 和野生沉香(II) 的GC-MS (A) 和UPLC-Q/TOF-MS (B) 指紋圖譜

    表2 通體香和野生沉香UPLC-Q/TOF-MS共有成分及差異化合物信息

    Table 2 Overview of common compounds and potential markers from agarwood induced by Agar-Wit and natural agarwood analyzed in UPLC-Q/TOF-MS total ion chromatograms common models fingerprint

    序號tR/min[M+H]+誤差(×10?6)MSE碎片離子峰取代基分子式類型名稱差異成分A環(huán)B環(huán) 1*5.06340.259 0 未鑒定 2*5.37319.118 20301.108 0, 283.096 8, 255.102 04×OH C17H18O6THPEC2-(2-phenylethyl)-5,6,7,8-tetrahydroxy-5,6,7,8-tetrahydrochromone(isomer) 3*5.55407.790 8 未鑒定 4*6.04319.118 0?0.6301.108 9, 283.096 8, 255.105 44×OH C17H18O6THPEC2-(2-phenylethyl)-5,6,7,8-tetrahydroxy-5,6,7,8-tetrahydrochromone(isomer) 5*6.30319.118 0?0.6301.108 4, 283.096 5, 255.101 94×OH C17H18O6THPEC2-(2-phenylethyl)-5,6,7,8-tetrahydroxy-5,6,7,8-tetrahydrochromone(isomer) 6*8.58329.102 3?0.6301.107 8, 301.107 8, 137.059 82×-O-OH; OCH3C18H16O6DEPECoxidoagarochromone C(isomer) 7*11.81331.117 2?3.0313.107 3, 285.112 0, 121.064 92×OH; -O-OCH3C18H18O6EPEC2-[2-(4'-methoxyphenyl)ethyl]-6,7-epoxy-5,8-dihydroxy-5,6,7,8-tetrahydrochromone(isomer) 8*12.04301.109 04.6283.097 8, 255.103 8, 227.109 12×OH; -O- C17H16O5EPEC2-(2-phenylethyl)-6,7-epoxy-5,8-dihydroxy-5,6,7,8-tetrahydrochromone(isomer) 9*14.33313.107 60283.112 5, 257.120 1, 121.068 52×-O-OCH3C18H16O5DEPECoxidoagarochromone B(isomer)L1 10*14.69283.097 10.3255.102 0, 227.107 1, 192.042 02×-O- C17H14O4DEPECoxidoagarochromone A(isomer)L2 11*17.47297.113 83.7137.060 0 OH; OCH3C18H16O4FTPEC2-[2-(2′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 12*17.47327.123 20191.071 0, 137.059 9OCH3OH; OCH3C19H18O5FTPEC6-methoxy-2-[2-(3′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 13*17.72327.123 40.6191.070 6, 137.060 2OCH3OH; OCH3C19H18O5FTPEC6-methoxy-2-[2-(3′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 14*18.08297.113 83.7137.060 0 OH; OCH3C18H16O4FTPEC2-[2-(2′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 15*19.06297.112 70137.060 6 OH; OCH3C18H16O4FTPEC2-[2-(2′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 16*19.71267.102 72.2137.023 9, 176.047 4OH C17H14O3FTPEC2-[2-(4′-methoxyphenyl)ethyl]chromone(isomer)L3 17*20.94341.139 10.6121.069 12×OCH3OCH3C20H20O5FTPEC6,7-dimethoxy-2-[2-(4′-methoxyphenyl)ethyl]chromone(isomer)L4 18*21.61311.128 50.6121.067 8OCH3OCH3C19H18O4FTPEC6-methoxy-2-[2-(4′-methoxyphenyl)ethyl]chromone(isomer)L5 19*24.16251.107 61.6173.060 1, 121.029 3 C17H14O2FTPEC2-(2-phenylethyl)chromoneL6 20*24.47311.128 10.6121.065 0OCH3OCH3C19H18O4FTPEC6-methoxy-2-[2-(4′-methoxyphenyl)ethyl]chromone(isomer)L7 21*24.88281.117 90.3121.069 2 OCH3C18H16O3FTPEC2-[2-(4'-methoxyphenyl)ethyl]chromone(isomer)L8 22*25.14233.154 5 未鑒定 23*25.40255.154 6 未鑒定 L9 24*26.55549.191 1?0.3531.183 4, 458.136 8, 440.126 0; 283.097 5 C34H28O7PEC dimers 259.26301.107 60283.097 4, 255.102 52×OH; -O- C18H16O5DEPEC2-(2-phenylethyl)-6,7-epoxy-5,8-dihydroxy-5,6,7,8-tetrahydrochromone(isomer) 2610.93367.094 5?0.8349.084 2, 331.118 2, 313.107 63×OH; ClOCH3C18H20ClO6THPEC8-chloro-5,6,7-trihydroxy-5,6,7,8-tetrahydro-2-[2-(4′-methoxyphenyl)ethyl]chromone(isomer)

    續(xù)表2

    序號tR/min[M+H]+誤差(×10?6)MSE碎片離子峰取代基分子式類型名稱差異成分 A環(huán)B環(huán) 279.82361.129 72.8137.060 7-O-; OH; OCH3OH; OCH3C19H21O7EPEC2-[2-(3′-hydroxy-4′-methoxyphenyl)ethyl]-7,8-epoxy-5-methoxy-6-hydroxy-5,6,7,8- tetrahydrochromone(isomer) 2811.24337.084 1?0.6319.073 7, 301.062 6, 283.096 73×OH; Cl C17H18ClO5THPEC8-chloro-5,6,7-trihydroxy-5,6,7,8-tetrahydro-2-(2-phenylethyl)chromone(isomer) 2912.42313.107 70.3285.112 8 OCH3C18H16O5DEPECoxidoagarochromone B(isomer) 3016.75267.102 51.5189.073 3, 176.047 5OH C17H14O3FTPEC2-[2-hydroxy-2-(2-phenylethyl)]chromone(isomer) 3116.05267.102 10173.056 2, 161.059 9 OHC17H14O3FTPEC2-[2-(3′-hydroxyphenyl)ethyl]chromone(isomer) 3216.69297.113 01.0137.060 5 OH; OCH3C18H16O4FTPEC2-[2-(2′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 3317.11297.112 90.7137.059 6 OH; OCH3C18H16O4FTPEC2-[2-(2′-hydroxy-4′-methoxyphenyl)ethyl]chromone(isomer) 3418.34327.125 77.6137.059 3OCH3OH; OCH3C19H18O5FTPEC6-methoxy-2-[2-(3′-hydroxy-4′-methoxyphenyl) ethyl]chromone(isomer) 3523.57281.117 4?1.4151.039 3OCH3 C18H16O3FTPEC6-methoxy-2-(2-phenylethyl)chromone(isomer) 3625.88297.112 90.7206.057 5, 191.034 5, 167.034 6OH; OCH3 C18H16O4FTPEC6-hydroxy-7-methoxy-2-(2-phenylethyl)chromone(isomer)

    DEPEC-雙氧環(huán)四氫色酮 EPEC-單氧環(huán)四氫色酮 THPEC-四羥基四氫2-(2-苯乙基)色酮 FTPEC-簡單2-(2-苯乙基)色酮 PEC dimers-2-(2-苯乙基)色酮二聚體;“*”為指紋圖譜共有成分;L1~L9代表顯著差異性成分

    DEPEC-diepoxy-2-(2-phenylethyl)chromones EPEC-epoxy-2-(2-phenylethyl)chromones THPEC-5,6,7,8-tetrahydro-2-(2-phenylethyl)chromones FTPEC-flindersia type 2-(2-phenylethyl)chromones PEC dimers 2-(2-phenylethyl)chromone dimers;“*”the common components in the fingerprint;L1—L9 represent significantly different components

    圖中編號與表2對應

    根據指紋圖譜共有峰分析結果可知,通體香的GC-MS指紋圖譜共有峰(RGC-tt,圖4-A)和野生沉香的GC-MS指紋圖譜共有峰(RGC-ys,圖4-B)均為32個。通過檢索NIST17譜庫并結合已有文獻對共有峰進行解析,指認了其中20個共有成分,見表3。包括9個倍半萜類化合物、8個小分子芳香化合物、2個2-(2-苯乙基)色酮類化合物以及1個酚酸類化合物。由此可見,通體香和野生沉香含有的主要揮發(fā)性成分一致。

    圖中編號與表3對應

    3.3 差異性分析

    為了進一步探究通體香和野生沉香化學組成的差異,采用有監(jiān)督的正交偏最小二乘判別分析(orthogonal partial least squares discriminant analysis,OPLS-DA)模型分析數據。結合OPLS-DA模型與獨立樣本檢驗,根據VIP值、S-Plot(即PC1方差貢獻率(p[1])和變量相關系數(p(corr)[1])以及值(閾值0.05)篩選差異變量。

    根據32批沉香樣品的UPLC-Q/TOF-MS分析結果,以VIP值>8、分布在S-plot左下角(p[1]<?0.1且p(corr)[1]<?0.5)以及<0.05為篩選條件(圖5-A),得到9個顯著差異性成分L1~L9,其中,L1和L2屬于三元氧環(huán)四氫色酮(DEPECs),L3、L4、L5、L6、L7以及L8屬于2-(2-苯乙基)色酮(FTPECs)。但是由表2可知,該9個成分均為通體香和野生沉香UPLC-Q/TOF-MS指紋圖譜中的共有成分,且從S-plot圖(圖5)可發(fā)現,L1~L9等9個差異成分均位于X軸負半軸,表明其在通體香中的含量均顯著高于野生沉香。根據32批沉香樣品的GC-MS分析結果,以VIP值>1且分布在S-plot左下角(p[1]<?0.07且p(corr)[1]<?0.15)和右上角(p[1]>0.045且p(corr)[1]>0.5)以及<0.05為篩選條件(圖5-B),得到12個顯著差異性成分G1~G12,包括3個小分子芳香化合物(G1、G2、G5)、5個倍半萜類化合物(G4、G7~G10)、2個2-(2-苯乙基)色酮類化合物(G11、G12)、1個烯烴類(G6)以及1個酚酸類化合物(G3),其中顯著差異性成分2-(2-苯乙基)色酮(G15)與UPLC-Q/TOF-MS分析所得顯著差異性成分2-(2-苯乙基)色酮(L6)相同。但是由表3可知,該12個成分均為通體香和野生沉香GC指紋圖譜中的共有成分,且從S-plot圖(圖5)可發(fā)現,2,4-二叔丁基苯酚(G3)、6,7-二甲基- 1,2,3,5,8,8a-六氫萘(G6)、以及去氫蜂斗菜酮(G7)等3個成分在野生沉香中的含量顯著高于通體香,其余9個成分在通體香中的含量均顯著高于野生沉香。

    表3 通體香和野生沉香GC-MS指紋圖譜已鑒定共有成分及差異化合物信息

    Table 3 Overview of identified common compounds and potential markers from agarwood induced by Agar-Wit and natural agarwood analyzed in GC-MS total ion chromatograms common models’ fingerprint

    峰號tR/min化合物名稱碎片離子 (m/z)匹配度分子式差異成分 16.02苯甲醛105.1, 77.1, 51.198.3C7H6OG1 39.75芐基丙酮148.1, 105.1, 91.1, 43.199.3C10H12OG2 510.744-甲氧基苯甲醛135.0, 77.0, 63.096.0C8H8O2 715.152,4-二叔丁基苯酚206.0, 191.0, 57.095.5C14H22OG3 818.51α-檀香醇206.0, 191.0, 57.080.1C15H24OG4 920.624-(4-甲氧基苯基)-2-丁酮178.0, 121.0, 43.093.8C11H14O2G5 1021.75沉香螺旋醇204.2, 161.1, 119.1, 59.193.1C15H26O 1122.57(–)-馬兜鈴烯161.1, 105.191.190.1C15H24 1223.422-((2R,8R,8aS)-8,8a-二甲基-1,2,3,4,6,7,8,8a-八氫萘-2-甲基)異丙基-2-醇204.2, 161.1, 59.190.7C15H26O 1323.96(4aS,8R)-4a,8-二甲氧基-4,4a,5,6,7,8-六氫萘-2(3H)-酮136.1, 121.0, 107.185.5C12H18O 1426.121,3-二(1-甲基乙基)苯162.0, 147.0, 105.0, 59.080.9C12H18 1526.656,7-二甲基-1,2,3,5,8,8a-六氫萘162.2, 147.1, 59.196.4C12H14O4G6 1940.32(4aR-cis)-4,4a,5,6,7,8-六氫-4a,5-二甲基-3-(1-甲基亞乙基)-2(3H)-萘175.0, 161.0, 147.089.0C15H22O 2143.08[3S-(3α,4aα,5α)]-3,4,4a,5,6,7-六氫-4a,5-二甲基-3-(1-甲基乙烯基)-1(2H)-萘177.1, 147.0, 91.179.9C15H22OG7 2653.68(4aR,5S)-1-羥基-4a,5-二甲基-3-亞丙基-4,4a,5,6-四氫萘-2(3H)-酮232.2, 217.0, 91.1, 77.183.1C15H20O2G8 2753.945,8-二羥基-4a-甲基-4,4a,4b,5,6,7,8,8a,9,10-十氫-2(3H)-菲酮185.2, 105.1, 91.182.0C15H22O3G9 2957.88(2aS,3aR,5aS,9bR)-2a,5a,9-三甲基-2a,4,5,5a,6,7,8,9b-八氫-2H-萘并 [1,2-b]環(huán)氧基[2,3-c]呋喃261.0, 105.0, 91.067.9C15H22O2G10 3067.022,2′-亞甲基雙[6-(1,1-二甲基乙基)-4-甲基苯酚340.0, 177.0, 161.097.0C23H32O2 3168.022-(2-苯乙基)色酮250.1, 91.194.7C17H14O2G11 3272.436-甲氧基-2-(2-苯乙基)色酮280.0, 91.083.4C18H16O3G12

    G1~G12代表顯著差異性成分

    G1-G12 represent significantly different components

    圖5 通體香和野生沉香UPLC-Q/TOF-MS和GC-MS的S-Plot圖

    將以上20個差異成分導入到Metaboanalyst軟件中繪制熱圖,距離度量方法為皮爾森積矩相關系數,聚類方法為平均連接法,結果見圖6,通體香和野生沉香各聚為一類??梢?,20個差異成分可以區(qū)分通體香和野生沉香,即通體香和野生沉香組間化學成分的差異遠大于其組內差異。通過熱圖的顏色變化發(fā)現,除了揮發(fā)性成分中G3、G6和G7在野生沉香中的含量顯著高于通體香,其余17個成分在通體香中的含量均顯著高于野生沉香,且具有顯著差異的2-(2-苯乙基)色酮類化合物均在通體香中具有較高的含量。

    4 討論

    4.1 樣品前處理方法的選擇

    本研究在UPLC-Q/TOF-MS的樣品前處理過程中對不同濃度甲醇溶液和乙醇溶液的熱回流提取效果進行了比較。結果發(fā)現,50%乙醇溶液和無水甲醇的提取效果較好,但采用無水甲醇溶解待測樣品時對其極性較大成分有溶劑效應,導致在UPLC上分離度較差。因此,本研究選擇50%乙醇溶液進行熱回流作為UPLC-Q/TOF-MS的樣品前處理方式。此外,本研究GC-MS的樣品前處理采用了加速溶劑萃取法代替水蒸氣蒸餾法,以水作為提取溶劑,通過提高溫度與壓力,在120 ℃提取沉香揮發(fā)性成分,可實現名貴藥材沉香的微量取樣,及其揮發(fā)性成分的高效、標準化提取。

    4.2 通體香與野生沉香的主要化學成分相同

    中藥指紋圖譜是采用一定分析手段得到的反應該中藥主要化學成分的共有峰圖譜,作為一種半定量鑒別手段,常用于中藥質量控制[17]。本研究發(fā)現,通體香和野生沉香UPLC-Q/TOF-MS指紋圖譜共有峰均為24個,主要為2-(2-苯乙基)色酮類化合物;而二者GC-MS指紋圖譜共有峰均為32個,主要為小分子芳香化合物和倍半萜。由于倍半萜和2-(2-苯乙基)色酮類等成分均是因白木香樹受到外界傷害產生防御反應后產生的防御性物質[2]??梢?,通體香和野生沉香在形成過程中產生的主要化學成分相同,這可能也是通體香乙醇提取物抗胃潰瘍[8]、抗結腸炎[9]等作用與野生沉香作用相當的原因。

    中藥指紋圖譜可全面反映中藥所含的物質群,其相似度系數反映了不同批次樣品間化學組成的整體波動程度,即中藥材質量的穩(wěn)定性[17-18]。沉香化學成分與其結香時間、生長環(huán)境和結香原因等因素有關[19]。Gao等[20]基于GC-MS指紋圖譜分析發(fā)現野生沉香和不同人工結香方法所產沉香的主要揮發(fā)性成分一致,但是不同人工結香方法所產沉香的指紋圖譜相似度較低,僅為0.180~0.811。而本研究發(fā)現16批采用通體結香技術誘導白木香樹所產沉香UPLC-Q/TOF-MS指紋圖譜和GC-MS指紋圖譜相似度較高,分別為0.598~0.869和0.846~0.995,且相似度高于野生沉香??梢?,不同批次通體香的化學組成相對于野生沉香更穩(wěn)定,可能與通體結香技術所采用的標準化結香、種植方式有關[21]。而16批野生沉香化學成分的相似度較16批通體香化學成分相似度小,可能與野生沉香生長環(huán)境、結香原因和產地加工等復雜因素有關,且其結香時間通常為數年至數十年不等[22-23]。

    圖6 通體香和野生沉香差異成分熱圖分析

    4.3 通體香中部分主要化學成分的含量高于野生沉香

    中藥材化學成分是其質量優(yōu)劣和藥理活性的物質基礎[24-25]。本研究基于UPLC-Q/TOF-MS和GC-MS并采用多元變量統計分析方法,從通體香和野生沉香中找出20個差異成分,但該20個成分均為二者的共有成分,在通體香和野生沉香中的含量存在差異。其中α-santalol等5個倍半萜和5個小分子芳香物質以及8個2-(2-苯乙基)色酮類成分在通體香中的含量均顯著高于野生沉香。這可能是通體香乙醇提取物對家鴿的止嘔作用[5]、對四氯化碳致小鼠急性肝損傷的保護作用[26]和改善異丙腎上腺素誘導的小鼠心肌缺血作用[7]優(yōu)于野生沉香的主要原因。其中,沉香乙醇提取物可通過降低多巴胺D2和5-羥色胺2A發(fā)揮止嘔作用[5]。α-santalol是多巴胺D2和5-羥色胺2A受體拮抗劑[27],而本研究發(fā)現通體香中的α-santalol比野生沉香高,可進一步解釋通體香止嘔作用優(yōu)于野生沉香。但是通體香中其他含量較高的化學成分及其藥理作用與機制還有待開展深入研究。

    此外,以本研究篩選出的20個差異成分通過聚類分析可將通體香和野生沉香各自聚為一類,可見20個主要成分在通體香和野生沉香組間的差異遠大于其組內差異,因此可作為通體香和野生沉香的鑒定依據。

    利益沖突 所有作者均聲明不存在利益沖突

    [1] 中國藥典 [S]. 一部2020: 192-193.

    [2] 張爭, 楊云, 魏建和, 等. 白木香結香機制研究進展及其防御反應誘導結香假說 [J]. 中草藥, 2010, 41(1): 156-159.

    [3] Liu Y Y, Chen H Q, Yang Y,. Whole-tree agarwood-inducing technique: An efficient novel technique for producing high-quality agarwood in cultivatedtrees [J]., 2013, 18(3): 3086-3106.

    [4] 李浩洋, 張憲臣, 張小花, 等. “通體結香技術”所結人工沉香的質量評價 [J]. 理化檢驗: 化學分冊, 2017, 53(4): 470-473.

    [5] 王燦紅, 彭德乾, 劉洋洋, 等. 沉香醇提物對順鉑誘導家鴿嘔吐的止嘔作用研究 [J]. 中國現代中藥, 2020, 22(10): 1611-1617.

    [6] 王燦紅, 王帥, 彭德乾, 等. 通體沉香醇提物對氟尿嘧啶致小鼠肝損傷的保護作用 [J]. 國際藥學研究雜志, 2018, 45(3): 187-197.

    [7] 王燦紅, 彭德乾, 劉洋洋, 等. 沉香醇提物對大鼠氣滯血瘀心肌缺血的保護作用 [J]. 藥物評價研究, 2019, 42(7): 1279-1285.

    [8] 劉洋洋, 王帥, 周岳, 等. “通體結香技術”所產沉香提取物對胃腸動力及胃潰瘍的作用 [J]. 國際藥學研究雜志, 2016, 43(6): 1076-1081.

    [9] Wang C H, Wang S, Peng D Q,. Agarwood extract mitigates intestinal injury in fluorouracil-induced mice [J]., 2019, 42(7): 1112-1119.

    [10] Li Y B, Sheng N, Wang L L,. Analysis of 2-(2-phenylethyl)chromones by UPLC-ESI-QTOF-MS and multivariate statistical methods in wild and cultivated agarwood [J]., 2016, 17(5): 771.

    [11] 劉洋洋. 通體結香技術產沉香的質量特征研究 [D]. 北京: 北京協和醫(yī)學院, 2018.

    [12] Want E J, Masson P, Michopoulos F,. Global metabolic profiling of animal and human tissues via UPLC-MS [J]., 2013, 8(1): 17-32.

    [13] Tsugawa H, Cajka T, Kind T,. MS-DIAL: Data-independent MS/MS deconvolution for comprehensive metabolome analysis [J]., 2015, 12(6): 523-526.

    [14] Zhang J, Yang W, Li S,. An intelligentized strategy for endogenous small molecules characterization and quality evaluation of earthworm from two geographic origins by ultra-high performance HILIC/QTOF MS(E) and Progenesis QI [J]., 2016, 408(14): 3881-3890.

    [15] 張倩, 霍會霞, 孫慧, 等. LC-IT-TOF-MSn快速鑒別沉香中2-(2-苯乙基)色酮類成分 [J]. 中國現代中藥, 2019, 21(7): 869-877.

    [16] Yang J L, Dong W H, Kong F D,. Characterization and analysis of 2-(2-phenylethyl)-chromone derivatives from agarwood () by artificial holing for different times [J]., 2016, 21(7): E911.

    [17] 任德權. 中藥指紋圖譜質控技術的意義與作用 [J]. 中藥新藥與臨床藥理, 2001, 12(3): 135-140.

    [18] 于燕喬, 史大卓, 曲華, 等. 定性研究方法在中醫(yī)藥領域的應用與思考[J]. 世界中醫(yī)藥, 2021, 16(13): 2058-2060.

    [19] Naziz P S, Das R, Sen S. The scent of stress: Evidence from the unique fragrance of agarwood [J]., 2019, 10: 840.

    [20] Gao X X, Xie M R, Liu S F,. Chromatographic fingerprint analysis of metabolites in natural and artificial agarwood using gas chromatography-mass spectrometry combined with chemometric methods [J]., 2014, 967: 264-273.

    [21] 周亞非, 魏建和, 孫佩文, 等. 中國沉香科技創(chuàng)新與產業(yè)發(fā)展的現狀及思考 [J]. 中國現代中藥, 2017, 19(8): 1051-1056.

    [22] Tan C S, Isa N M, Ismail I,. Agarwood induction: Current developments and future perspectives [J]., 2019, 10: 122.

    [23] 黃俊卿, 魏建和, 張爭, 等. 沉香結香方法的歷史記載、現代研究及通體結香技術 [J]. 中國中藥雜志, 2013, 38(3): 302-306.

    [24] 張鐵軍, 許浚, 申秀萍, 等. 基于中藥質量標志物(Q-Marker)的元胡止痛滴丸的“性-效-物”三元關系和作用機制研究 [J]. 中草藥, 2016, 47(13): 2199-2211.

    [25] 張萌, 封亮, 賈曉斌. 基于生物活性與效應基準的中藥質量評價技術發(fā)展現狀與展望[J]. 世界中醫(yī)藥, 2020, 15(15): 2234-2239.

    [26] 王燦紅, 王帥, 彭德乾, 等. 沉香醇提物對四氯化碳致小鼠急性肝損傷的保護作用 [J]. 中國現代中藥, 2017, 19(8): 1091-1096.

    [27] Okugawa H, Ueda R, Matsumoto K,. Effects of sesquiterpenoids from “Oriental incenses” on acetic acid-induced writhing and D2 and 5-HT2A receptors in rat brain [J]., 2000, 7(5): 417-422.

    Comparative analysis of chemical constituents between the agarwood produced by Whole-tree agarwood-inducing technique and natural agarwood based on the combination of non-targeted metabolomics and fingerprint

    YU Meng1, 2, HOU Wen-cheng1, FENG Jian1, CHEN Lan1, YANG Yun1, LIU Yang-yang1, 2, WEI Jian-he1, 2

    1. Hainan Provincial Key Laboratory of Resources Conservation and Development of Southern Medicine & Key Laboratory of State Administration of Traditional Chinese Medicine for Agarwood Sustainable Utilization, Hainan Branch of the Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences and Peking Union Medical College, Haikou 570311, China 2. Key Laboratory of Bioactive Substances and Resources Utilization of Chinese Herbal Medicine, Ministry of Education & National Engineering Laboratory for Breeding of Endangered Medicinal Materials, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100193, China

    To provide scientific basis for the advanced application and popularization of agarwood produced by Whole-tree agarwood-inducing technique (Agar-Wit), the compositional differences between agarwood induced by Agar-Wit and natural agarwood were analyzed and compared in the study.The ethanol extracts and volatile components of samples were analyzed by ultra-performance liquid chromatography tandem quadrupole time-of-flight mass spectrometry (UPLC-Q/TOF-MS) and gas chromatography-mass spectrometry (GC-MS) combined with non-targeted metabolomics and fingerprint in this study.The common components of agarwood induced by Agar-Wit were consistent with natural agarwood. The similarity coefficients among Agar-Wit agarwood were higher than that in natural agarwood. In addition, the content of three low-molecular-weight aromatic compounds, four kinds of sesquiterpenes and eight 2-(2-phenylethyl)chromones were richer in agarwood induced by Agar-Wit.The main chemical composition of agarwood induced by Agar-Wit, which is more stable, is proved to be in accordance with wild agarwood.

    agarwood; whole-tree agarwood-inducing technique (Agar-Wit); wild agarwood; non-targeted metabolomics; fingerprint

    R284.1

    A

    0253 - 2670(2022)18 - 5625 - 11

    10.7501/j.issn.0253-2670.2022.18.005

    2022-02-20

    國家重點研發(fā)計劃項目(2018YFC1706403);海南省自然科學基金高層次人才項目(2019RC340);現代農業(yè)產業(yè)技術體系建設專項(CARS-21);中國醫(yī)學科學院醫(yī)學與健康科技創(chuàng)新工程項目(2016-I2M-2-003);2021年國家藥品標準提高課題(2021Z14)

    郁 萌,碩士研究生,研究方向為沉香質量分析。E-mail: jokeryml@163.com

    劉洋洋,研究員,碩士生導師,研究方向為中藥質量控制技術及其應用。Tel/Fax: (0898)31589000 E-mail: eadchris@163.com

    魏建和,研究員,博士生導師,研究方向為藥用植物基因資源、分子育種及次生代謝產物調控研究。Tel: (010)57833016 E-mail: wjianh@263.net

    [責任編輯 王文倩]

    猜你喜歡
    結香通體乙基
    痛風
    通體結冰的球
    結香的奧秘
    一種通體大理石瓷磚的制作方法
    佛山陶瓷(2018年6期)2018-09-14 10:54:46
    結香屬植物化學成分及藥理作用的研究進展
    中成藥(2018年1期)2018-02-02 07:20:11
    中國沉香基地及通體結香技術
    結香:枝條常被打結的神奇樹種
    硫酸鋅電解液中二(2-乙基己基)磷酸酯的測定
    2-羧乙基苯基次膦酸的胺化處理及其在尼龍6中的阻燃應用
    雙[2-(5-硝基-2H-四唑基)-2,2-二硝乙基]硝胺的合成與量子化學計算
    火炸藥學報(2014年5期)2014-03-20 13:17:51
    久久久久久人人人人人| 亚洲av成人av| 18禁国产床啪视频网站| 啦啦啦免费观看视频1| av欧美777| 久热这里只有精品99| 99精品欧美一区二区三区四区| 精品乱码久久久久久99久播| 国产精品一区二区三区四区久久 | 久久婷婷人人爽人人干人人爱| 国产又色又爽无遮挡免费看| 亚洲第一青青草原| 香蕉av资源在线| 国产av一区二区精品久久| 久久精品国产清高在天天线| 欧美激情极品国产一区二区三区| 亚洲精品国产精品久久久不卡| 免费女性裸体啪啪无遮挡网站| 午夜福利在线观看吧| av福利片在线| 黄色片一级片一级黄色片| 久久精品国产99精品国产亚洲性色| 性欧美人与动物交配| 成人18禁在线播放| 国内精品久久久久精免费| 欧美大码av| 日日摸夜夜添夜夜添小说| 熟妇人妻久久中文字幕3abv| 黑人巨大精品欧美一区二区mp4| 成人av一区二区三区在线看| 亚洲成人精品中文字幕电影| 午夜福利一区二区在线看| 波多野结衣高清作品| 久久中文字幕一级| 性色av乱码一区二区三区2| 国产精品 国内视频| 久9热在线精品视频| 观看免费一级毛片| 免费看日本二区| 亚洲 欧美一区二区三区| 两个人看的免费小视频| 淫妇啪啪啪对白视频| 人人妻人人看人人澡| 中文字幕高清在线视频| 成年版毛片免费区| 老司机靠b影院| 久久久国产欧美日韩av| 国产不卡一卡二| 在线天堂中文资源库| 18禁美女被吸乳视频| 十八禁人妻一区二区| 日韩一卡2卡3卡4卡2021年| 国产亚洲欧美精品永久| 精品久久久久久久人妻蜜臀av| 欧美绝顶高潮抽搐喷水| 亚洲电影在线观看av| 美女免费视频网站| 国产黄a三级三级三级人| 69av精品久久久久久| 国产欧美日韩一区二区精品| 久久香蕉激情| 国产精品久久电影中文字幕| 啦啦啦观看免费观看视频高清| 日日爽夜夜爽网站| 亚洲一区二区三区不卡视频| 侵犯人妻中文字幕一二三四区| 亚洲精品中文字幕在线视频| 97碰自拍视频| 婷婷六月久久综合丁香| 亚洲无线在线观看| 国产精品久久视频播放| 国内精品久久久久久久电影| 国产精品 欧美亚洲| 久久人妻av系列| 黑人操中国人逼视频| 在线观看66精品国产| 精品人妻1区二区| 色综合亚洲欧美另类图片| 国产精品自产拍在线观看55亚洲| 亚洲国产精品成人综合色| 18禁黄网站禁片午夜丰满| 久久伊人香网站| 国产精品电影一区二区三区| 天天一区二区日本电影三级| 麻豆久久精品国产亚洲av| av中文乱码字幕在线| 日韩欧美一区二区三区在线观看| 国产成年人精品一区二区| 免费看美女性在线毛片视频| 午夜精品久久久久久毛片777| 成年免费大片在线观看| 午夜免费鲁丝| 欧美日韩瑟瑟在线播放| 黄色 视频免费看| 午夜福利18| www.999成人在线观看| 久久亚洲真实| 少妇的丰满在线观看| 精品一区二区三区av网在线观看| 久久青草综合色| 亚洲午夜精品一区,二区,三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧美日韩在线播放| 美女免费视频网站| 国产一区二区在线av高清观看| 搡老妇女老女人老熟妇| 黄色视频,在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美黑人巨大hd| 热99re8久久精品国产| 国产精品永久免费网站| 久久国产亚洲av麻豆专区| 久久国产精品影院| 女人高潮潮喷娇喘18禁视频| 两个人免费观看高清视频| 免费av毛片视频| 国产精品爽爽va在线观看网站 | 18禁裸乳无遮挡免费网站照片 | 香蕉久久夜色| 精品一区二区三区视频在线观看免费| 亚洲专区中文字幕在线| 久久狼人影院| 国产av不卡久久| 男女之事视频高清在线观看| 精品熟女少妇八av免费久了| 精品久久久久久成人av| 精品国内亚洲2022精品成人| 精品福利观看| 久久精品人妻少妇| 婷婷六月久久综合丁香| videosex国产| 大型av网站在线播放| 最近在线观看免费完整版| 亚洲欧美精品综合久久99| 久久久久久久久中文| 女同久久另类99精品国产91| 国产野战对白在线观看| 91麻豆精品激情在线观看国产| 哪里可以看免费的av片| 男女视频在线观看网站免费 | 可以在线观看的亚洲视频| 中文在线观看免费www的网站 | 亚洲国产精品合色在线| 亚洲专区中文字幕在线| 国产精品,欧美在线| 久久久久免费精品人妻一区二区 | 亚洲精品在线美女| 欧美日本亚洲视频在线播放| 亚洲av成人av| 99国产精品一区二区蜜桃av| 日本一区二区免费在线视频| 91大片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 午夜两性在线视频| 国产又色又爽无遮挡免费看| 国产精品精品国产色婷婷| 欧美色视频一区免费| 法律面前人人平等表现在哪些方面| 一级毛片精品| 精品国产超薄肉色丝袜足j| 18禁黄网站禁片免费观看直播| 亚洲精品美女久久av网站| 18禁美女被吸乳视频| 亚洲成人久久爱视频| 午夜久久久久精精品| 亚洲va日本ⅴa欧美va伊人久久| 国产高清视频在线播放一区| 男女做爰动态图高潮gif福利片| 琪琪午夜伦伦电影理论片6080| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 一本综合久久免费| 一区二区三区高清视频在线| 他把我摸到了高潮在线观看| 18禁美女被吸乳视频| 狂野欧美激情性xxxx| 国产伦一二天堂av在线观看| 每晚都被弄得嗷嗷叫到高潮| 午夜福利一区二区在线看| 亚洲av电影在线进入| 国产精华一区二区三区| 国产精品永久免费网站| 成人国产综合亚洲| 国产av在哪里看| 欧美激情极品国产一区二区三区| 久久亚洲真实| 老司机午夜福利在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 女性被躁到高潮视频| 午夜影院日韩av| 淫秽高清视频在线观看| 一级作爱视频免费观看| 黑人欧美特级aaaaaa片| 女性被躁到高潮视频| 日韩免费av在线播放| 免费在线观看完整版高清| 天堂√8在线中文| 禁无遮挡网站| 成人18禁在线播放| 欧美不卡视频在线免费观看 | 一区二区三区激情视频| 村上凉子中文字幕在线| 国产精品九九99| av福利片在线| 人人澡人人妻人| 少妇 在线观看| 男女下面进入的视频免费午夜 | 亚洲成a人片在线一区二区| 亚洲人成网站高清观看| 欧美精品亚洲一区二区| 美女免费视频网站| 日韩欧美免费精品| 国产精品美女特级片免费视频播放器 | 欧美zozozo另类| 最近在线观看免费完整版| 制服人妻中文乱码| 亚洲欧美精品综合一区二区三区| 国产亚洲精品第一综合不卡| 亚洲男人天堂网一区| 大型av网站在线播放| 十分钟在线观看高清视频www| 波多野结衣巨乳人妻| 悠悠久久av| 亚洲第一青青草原| 久久久久久久精品吃奶| 人妻久久中文字幕网| 女性被躁到高潮视频| 国产一区在线观看成人免费| 中文字幕另类日韩欧美亚洲嫩草| 真人一进一出gif抽搐免费| 人人澡人人妻人| 熟女电影av网| 黄色女人牲交| 一个人观看的视频www高清免费观看 | 在线av久久热| 国产成+人综合+亚洲专区| 嫩草影视91久久| 亚洲熟妇熟女久久| 99精品久久久久人妻精品| 制服人妻中文乱码| 亚洲最大成人中文| 午夜福利18| 他把我摸到了高潮在线观看| 丝袜美腿诱惑在线| 久久久久国内视频| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| 亚洲电影在线观看av| 久久久久久久久久黄片| 久久精品91无色码中文字幕| 熟女少妇亚洲综合色aaa.| 真人做人爱边吃奶动态| 好男人电影高清在线观看| 久久久久久久午夜电影| 欧美人与性动交α欧美精品济南到| xxx96com| 不卡一级毛片| 亚洲精品色激情综合| tocl精华| 亚洲五月天丁香| 国产精品久久电影中文字幕| 国产精品免费视频内射| 日本 av在线| 国产成人精品久久二区二区91| 性欧美人与动物交配| 久久久久国产精品人妻aⅴ院| 波多野结衣av一区二区av| 叶爱在线成人免费视频播放| av片东京热男人的天堂| 热re99久久国产66热| 国语自产精品视频在线第100页| 亚洲自拍偷在线| 大型av网站在线播放| 日本三级黄在线观看| 成人国产综合亚洲| 日韩精品免费视频一区二区三区| 成在线人永久免费视频| 亚洲成a人片在线一区二区| 99久久99久久久精品蜜桃| 亚洲精品一区av在线观看| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 亚洲专区字幕在线| 搡老岳熟女国产| 亚洲av熟女| 亚洲精品国产一区二区精华液| 天天一区二区日本电影三级| 精品国产一区二区三区四区第35| 欧美激情高清一区二区三区| 精品久久蜜臀av无| 波多野结衣巨乳人妻| 女人高潮潮喷娇喘18禁视频| 国产精品自产拍在线观看55亚洲| 99久久精品国产亚洲精品| 天天一区二区日本电影三级| 亚洲片人在线观看| 一区二区三区精品91| 久久 成人 亚洲| 亚洲成人精品中文字幕电影| 亚洲黑人精品在线| 国产精品久久电影中文字幕| 久久香蕉精品热| 女生性感内裤真人,穿戴方法视频| 中文字幕人妻熟女乱码| 窝窝影院91人妻| 男人舔奶头视频| 久久精品国产亚洲av高清一级| 欧美黄色淫秽网站| 欧美精品啪啪一区二区三区| 男女下面进入的视频免费午夜 | 午夜福利成人在线免费观看| 国产黄片美女视频| 亚洲国产欧洲综合997久久, | 999精品在线视频| 99在线人妻在线中文字幕| 99热6这里只有精品| 国产伦一二天堂av在线观看| av在线播放免费不卡| 一区二区三区精品91| 丁香六月欧美| 亚洲自拍偷在线| 亚洲国产欧美一区二区综合| 亚洲在线自拍视频| 日韩免费av在线播放| 亚洲欧美一区二区三区黑人| 国产精品永久免费网站| 亚洲国产欧美日韩在线播放| av有码第一页| 美女 人体艺术 gogo| 又黄又爽又免费观看的视频| 宅男免费午夜| 在线观看免费日韩欧美大片| xxx96com| 欧美日本亚洲视频在线播放| 亚洲国产欧美一区二区综合| 久久午夜综合久久蜜桃| 国产一级毛片七仙女欲春2 | 一卡2卡三卡四卡精品乱码亚洲| 精品不卡国产一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲一区高清亚洲精品| 欧美精品啪啪一区二区三区| 一区福利在线观看| a级毛片在线看网站| 久久热在线av| 午夜福利高清视频| а√天堂www在线а√下载| 最新美女视频免费是黄的| 好男人电影高清在线观看| 亚洲aⅴ乱码一区二区在线播放 | 黑人欧美特级aaaaaa片| 欧美一区二区精品小视频在线| 一个人免费在线观看的高清视频| 欧美黄色淫秽网站| 女同久久另类99精品国产91| 成人亚洲精品一区在线观看| 亚洲成a人片在线一区二区| 在线观看免费视频日本深夜| 啦啦啦 在线观看视频| 午夜福利成人在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲人成伊人成综合网2020| 亚洲午夜精品一区,二区,三区| 亚洲成人国产一区在线观看| 久久这里只有精品19| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 欧美成狂野欧美在线观看| 中文字幕精品免费在线观看视频| 国产又色又爽无遮挡免费看| 精品欧美国产一区二区三| 久久精品亚洲精品国产色婷小说| 一级a爱片免费观看的视频| 1024手机看黄色片| 日韩国内少妇激情av| 免费av毛片视频| 国产国语露脸激情在线看| 亚洲专区字幕在线| 久久狼人影院| 亚洲精品久久国产高清桃花| 国产亚洲精品综合一区在线观看 | 成人三级黄色视频| 精品久久久久久久末码| 黄色毛片三级朝国网站| 亚洲精品一区av在线观看| 啪啪无遮挡十八禁网站| 欧美一级a爱片免费观看看 | 亚洲国产欧美日韩在线播放| 久久人妻福利社区极品人妻图片| 欧美另类亚洲清纯唯美| 男人舔女人的私密视频| av视频在线观看入口| 女性生殖器流出的白浆| www国产在线视频色| 在线观看www视频免费| 国产午夜精品久久久久久| 婷婷丁香在线五月| 精品不卡国产一区二区三区| 露出奶头的视频| 一本精品99久久精品77| 久久久久亚洲av毛片大全| 亚洲五月天丁香| 久久国产亚洲av麻豆专区| 午夜a级毛片| 亚洲成人久久爱视频| 一二三四在线观看免费中文在| avwww免费| 三级毛片av免费| 777久久人妻少妇嫩草av网站| 国产亚洲欧美在线一区二区| www国产在线视频色| 999久久久精品免费观看国产| 亚洲狠狠婷婷综合久久图片| 亚洲第一电影网av| 人人妻人人澡欧美一区二区| 91国产中文字幕| 熟妇人妻久久中文字幕3abv| 亚洲激情在线av| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三 | 国产熟女午夜一区二区三区| 麻豆一二三区av精品| 香蕉av资源在线| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| 制服丝袜大香蕉在线| 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 一本久久中文字幕| 免费在线观看视频国产中文字幕亚洲| 免费电影在线观看免费观看| 免费看日本二区| 日韩免费av在线播放| 久久天躁狠狠躁夜夜2o2o| 欧美乱码精品一区二区三区| 免费av毛片视频| 国产97色在线日韩免费| 亚洲色图av天堂| 一个人免费在线观看的高清视频| 精品久久久久久成人av| 天天躁夜夜躁狠狠躁躁| 一级片免费观看大全| 午夜激情av网站| 欧美成人午夜精品| 免费观看人在逋| 最近最新中文字幕大全电影3 | 少妇粗大呻吟视频| 女人高潮潮喷娇喘18禁视频| 此物有八面人人有两片| 91在线观看av| 两性午夜刺激爽爽歪歪视频在线观看 | 国产三级黄色录像| 亚洲成国产人片在线观看| 白带黄色成豆腐渣| 亚洲性夜色夜夜综合| 成人一区二区视频在线观看| 亚洲片人在线观看| 欧美黄色淫秽网站| 91字幕亚洲| 青草久久国产| 日韩国内少妇激情av| 久久久久久久午夜电影| 高清在线国产一区| 一级a爱片免费观看的视频| 99热6这里只有精品| 不卡一级毛片| 在线观看午夜福利视频| 我的亚洲天堂| 波多野结衣高清无吗| 久久午夜亚洲精品久久| 欧美不卡视频在线免费观看 | 欧美黑人欧美精品刺激| 中出人妻视频一区二区| 精品人妻1区二区| 日韩有码中文字幕| 精品国产亚洲在线| 久久精品91无色码中文字幕| 国产精品国产高清国产av| 免费高清在线观看日韩| 国产人伦9x9x在线观看| 精品人妻1区二区| 1024手机看黄色片| 日本成人三级电影网站| 午夜亚洲福利在线播放| 人人妻人人看人人澡| 精品乱码久久久久久99久播| 精品日产1卡2卡| 波多野结衣高清无吗| 18禁黄网站禁片午夜丰满| 亚洲中文字幕一区二区三区有码在线看 | 非洲黑人性xxxx精品又粗又长| 亚洲五月色婷婷综合| 丝袜美腿诱惑在线| 黄色a级毛片大全视频| 丰满的人妻完整版| 亚洲成a人片在线一区二区| 男人舔奶头视频| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 日韩av在线大香蕉| 777久久人妻少妇嫩草av网站| 日韩视频一区二区在线观看| 精品国内亚洲2022精品成人| 天堂√8在线中文| 国产激情欧美一区二区| 女生性感内裤真人,穿戴方法视频| 免费高清在线观看日韩| 亚洲熟妇熟女久久| 91老司机精品| 十八禁网站免费在线| 亚洲在线自拍视频| 亚洲国产看品久久| 久久性视频一级片| 岛国在线观看网站| 日本a在线网址| 免费在线观看视频国产中文字幕亚洲| 听说在线观看完整版免费高清| 亚洲国产中文字幕在线视频| 日韩欧美国产在线观看| av在线播放免费不卡| 亚洲国产中文字幕在线视频| 中文字幕人成人乱码亚洲影| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| cao死你这个sao货| 岛国在线观看网站| 狠狠狠狠99中文字幕| 色综合站精品国产| www.精华液| 久久国产精品男人的天堂亚洲| 99热6这里只有精品| 国产男靠女视频免费网站| 国产97色在线日韩免费| 大型黄色视频在线免费观看| 在线观看日韩欧美| 久久久久久九九精品二区国产 | 亚洲精品美女久久av网站| 一本精品99久久精品77| 国产私拍福利视频在线观看| 久久热在线av| 国产精品,欧美在线| 成年免费大片在线观看| 午夜影院日韩av| 精品欧美一区二区三区在线| 波多野结衣av一区二区av| xxxwww97欧美| 最好的美女福利视频网| 最新在线观看一区二区三区| 国产亚洲精品第一综合不卡| xxxwww97欧美| 在线视频色国产色| 午夜福利成人在线免费观看| 两个人免费观看高清视频| 亚洲avbb在线观看| 亚洲片人在线观看| 久久伊人香网站| 久久婷婷人人爽人人干人人爱| 99热这里只有精品一区 | 国产av又大| 欧美一区二区精品小视频在线| 精品熟女少妇八av免费久了| 亚洲全国av大片| 久久久精品国产亚洲av高清涩受| 伊人久久大香线蕉亚洲五| 国产精品美女特级片免费视频播放器 | 在线免费观看的www视频| 欧美久久黑人一区二区| 长腿黑丝高跟| 真人做人爱边吃奶动态| 欧美性猛交╳xxx乱大交人| 又黄又粗又硬又大视频| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品av在线| 免费女性裸体啪啪无遮挡网站| 国产伦在线观看视频一区| 亚洲七黄色美女视频| 国产黄片美女视频| 久久性视频一级片| 999久久久精品免费观看国产| 午夜精品在线福利| 老熟妇仑乱视频hdxx| 成年免费大片在线观看| 在线观看一区二区三区| 最近最新中文字幕大全电影3 | 国产精华一区二区三区| 欧美一区二区精品小视频在线| 18禁黄网站禁片午夜丰满| 亚洲美女黄片视频| 亚洲国产中文字幕在线视频| 欧美午夜高清在线| 我的亚洲天堂| 国产黄a三级三级三级人| netflix在线观看网站| 欧美成人性av电影在线观看| 欧美精品亚洲一区二区| www.www免费av| 少妇粗大呻吟视频| 日本撒尿小便嘘嘘汇集6| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美激情性xxxx| 少妇被粗大的猛进出69影院| 国产av一区在线观看免费| 深夜精品福利| 精品一区二区三区av网在线观看| 变态另类丝袜制服| 欧美日本视频| 免费搜索国产男女视频| 亚洲中文字幕一区二区三区有码在线看 | 窝窝影院91人妻| 欧美激情高清一区二区三区| 国产精品一区二区精品视频观看| 日本 欧美在线| 免费在线观看影片大全网站| 国产成人一区二区三区免费视频网站|