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

    基于UPLC-MS/MS代謝組學技術分析不同栽培模式下鐵皮石斛類黃酮化合物差異性

    2022-02-23 12:56:12金傳高任仙櫻潘曉軍吳志剛姜程曦
    中草藥 2022年4期

    李 帆,金傳高,任仙櫻,費 璇,趙 祺,潘曉軍,吳志剛*,姜程曦

    基于UPLC-MS/MS代謝組學技術分析不同栽培模式下鐵皮石斛類黃酮化合物差異性

    李 帆1, 2,金傳高3#,任仙櫻2, 4,費 璇2,趙 祺2, 4,潘曉軍2,吳志剛2*,姜程曦1, 4*

    1. 溫州大學生命科學研究院,浙江 溫州 325035 2. 溫州醫(yī)科大學藥學院,浙江 溫州 325035 3. 浙江高鼻子生物科技有限公司,浙江 樂清 325615 4. 安徽省九華山佛教醫(yī)藥研究所,安徽 池州 247100

    研究不同栽培模式下鐵皮石斛中類黃酮化合物積累的差異,為鐵皮石斛種植方式創(chuàng)新提供依據(jù)。運用超高效液相色譜串聯(lián)質譜(UPLC-MS/MS)廣泛靶向技術對3種栽培模式(大棚栽培、林下活樹附生和巖壁附生栽培)下的鐵皮石斛類黃酮物質進行檢測,使用主成分分析(principal component analysis,PCA)和正交偏最小二乘判別分析法(orthogonal partial least squares-discriminant analysis,OPLS-DA)鑒定差異代謝物。共鑒定202個類黃酮化合物,主要以黃酮(76個)和黃酮醇(60個)化合物為主,其中花青素物質矢車菊素-3--葡萄糖苷-7,3′-di--(6′′--芥子酰葡萄糖苷)為首次在鐵皮石斛中被檢測到;OPLS-DA篩選43個差異化合物,其中大棚組與林下組、大棚組與巖壁組間均具有23個差異物,林下組與巖壁組則具有27個差異物,差異物質對表征不同栽培模式的樣品貢獻率明顯。林下活樹附生和巖壁附生栽培明顯提升黃酮、類黃酮化合物的生物合成水平,保障了藥材質量,研究結果為新型栽培模式的推廣提供了理論基礎。

    鐵皮石斛;UPLC-MS/MS;類黃酮;代謝組學;栽培模式

    鐵皮石斛Kimura et Migo.為石斛屬多年生草本植物,入選“影響世界的中國植物”,被譽為中華“九大仙草”之首,為我國名貴中藥[1-2]。鐵皮石斛地上莖可加工成鐵皮楓斗,《神農(nóng)本草經(jīng)》列為上品,味甘,主傷中、下氣,補五臟虛勞、羸弱,強陰、厚腸胃、輕身延年之功效[3]。藥理研究表明,鐵皮石斛在提高人體免疫能力、抗衰老、抗腫瘤、改善糖尿病癥狀和萎縮性胃炎等方面具有顯著的療效[4-6]。

    浙江是我國鐵皮石斛的道地產(chǎn)區(qū)之一,鐵皮石斛已被列入“新浙八味”。自2000年以來產(chǎn)區(qū)在品種選育、組織培養(yǎng)、栽培措施等關鍵技術方面取得突破性的進展[7-9],促進了石斛產(chǎn)業(yè)快速的發(fā)展。特別是,近年來以林下活樹附生、巖壁附生等近野生栽培技術悄然興起,成為鐵皮石斛的新型的栽培模式,從而逐步使鐵皮石斛種植由大棚向森林回歸。然后,不同的栽培方式的光、溫、濕度等環(huán)境因素千差萬別,一定程度也會對鐵皮石斛生長、有效成分積累造成影響。已有研究表明,光質和低溫顯著鐵皮石斛種苗的發(fā)育、次生代謝產(chǎn)物的調控,但研究主要集中在氨基酸、多糖、生物堿等成分[10-11]。類黃酮是鐵皮石斛一類重要的活性成分,具有抗氧化、降血糖、降低膽固醇、保護心血管系統(tǒng)等藥理活性[12]。光、溫等環(huán)境脅迫是誘導植物類黃酮化合物生長合成的重要環(huán)境因素之一[13]。因而,推測不同栽培模式下因其環(huán)境差異性可能會對類黃酮的生物合成造成一定影響,從而影響鐵皮石斛的品質,而有關這方面未見報道。

    鑒于類黃酮物質重要的藥用價值,本實驗采用超高效液相色譜串聯(lián)質譜(UPLC-MS/MS)結合植物廣泛靶向代謝組學技術整體表征3種栽培模式下鐵皮石斛類黃酮化合物差異性,利用主成分分析(principal component analysis,PCA)和正交偏最小二乘判別分析法(orthogonal partial least squares-discriminant analysis,OPLS-DA)等多元統(tǒng)計手段鑒定主導這種差異性的標志性的代謝物,發(fā)掘潛在的規(guī)律,以為改進鐵皮石斛種植模式及提供高質量產(chǎn)品奠定基礎。

    1 儀器與材料

    超高效液相色譜(UFLC SHIMADZU CBM30,日本島津公司),4500 QTRAP型高分辨串聯(lián)質譜(美國應用生物系統(tǒng)公司)。主要試劑包括甲醇(LC-MS級,美國Merk公司)、乙腈(LC-MS級,美國Merk公司),水為超純水(Milli-Q SP Regent Water system,美國),MM 400型研磨儀(德國萊馳公司)。

    本課題組于2017年4月使用同一批種苗,在雁蕩山石斛谷基地(E121.06°、N28.39°)布置了大棚、林下活樹附生和巖壁栽培實驗。實驗材料于2019年11月采集該實驗條件的樣品,由溫州醫(yī)科大學吳志剛副研究員鑒定為鐵皮石斛Kimura et Migo.鮮條。每個樣品取自5株不同植株的鮮條,液氮下研碎,均勻混合。不同栽培模式均4次生物學重復,共12個樣品。此外,3種栽培模式樣品等量混合構成質控樣本,共3個重復,用于檢測儀器的重復性和精確度。所用樣品放置?80 ℃超低溫冰箱備用。

    2 方法

    2.1 樣品制備

    鐵皮石斛樣品經(jīng)真空冷凍干燥后,研磨儀(30 Hz,1.5 min)研磨至粉末狀;稱取100 mg的粉末,溶解于1.0 mL的70%甲醇水溶液中;溶解后的樣品4 ℃冰箱過夜,期間渦旋3次,以提高提取率;離心(轉速10 000×,10 min)后,上清液用微孔濾膜(0.22 μm)濾過后進樣。

    2.2 色譜和質譜采集條件

    色譜柱:ACQUITY UPLC HSS T3(100 mm×2.1 mm,1.8 μm),十八烷基鍵合硅膠為填充劑。流動相為0.04%乙酸超純水(A)-乙腈(含0.04%的乙酸(B),洗脫梯度:0~11.0 min,5%~95% B;11.0~12.0 min,95% B;12.1~15.0 min,5% B。體積流量0.4 mL/min,柱溫40 ℃,進樣量5 μL。

    采用正負離子模式檢測模式,霧化氣體為高純度氮氣(N2),電噴霧離子源(ESI)溫度550 ℃,離子化電壓5500 V,簾氣(curtain gas,CUR)207 kPa,輔助電壓60 kPa。按參考文獻[15]方法優(yōu)化去簇電壓和碰撞能。質量掃描范圍/100~1200,掃描時間0.2 s。

    2.3 數(shù)據(jù)處理

    MultiaQuant軟件對色譜峰的進行積分和校正,每個色譜峰的峰面積代表對應物質的相對含量,導出數(shù)據(jù)矩陣。利用三重四級桿質譜的多反應監(jiān)測模式(MRM)將獲得的二級質譜數(shù)據(jù)與數(shù)據(jù)庫進行比對,對化合物進行定性、定量分析。使用數(shù)據(jù)庫包括邁維公司自建MWDB、MassBank(http://www. massbank.jp/)、HMDB(http://www.hmdb.ca/)和METLIN(http://metlin.scripps.edu/index.php)等數(shù)據(jù)庫。R軟件ropls包對數(shù)據(jù)矩陣進行PCA和OPLS-DA分析,鑒定差異性代謝產(chǎn)物。使用R軟件pheatmap包進行聚類分析(hierarchical cluster analysis,HCA),以表征不同樣本代謝物的積累模式。

    3 結果與分析

    3.1 鐵皮石斛類黃酮化合物的鑒定

    UPLC-MS/MS首先以正、負離子模式對混樣質控QC樣品進行檢測,質控樣本TIC曲線重疊性高、質譜峰保留時間和強度一致,表明所建方法可靠、儀器穩(wěn)定性較好(圖1-A)。所獲得的代謝物高分辨二級質譜數(shù)據(jù)經(jīng)與公共數(shù)據(jù)庫比對,共鑒定202個類黃酮化合物,其中黃酮化合物76個、黃酮醇60個、黃酮碳糖苷21個、二氫黃酮16個、黃烷醇類10個,此外還包含少量異黃酮、二氫黃酮醇、原花青素和花青素等物質。其中,花青素物質矢車菊素-3--葡萄糖苷- 7,3′-di--(6′′--芥子酰葡萄糖苷)首次在鐵皮石斛中被報道。大棚、林下及巖壁附生栽培的樣品的UPLC-MS/MS基峰離子流圖如圖1-B~D所示,直觀顯示不同栽培模式下類黃酮化合物的存在較為明顯的差異。

    A-混樣質控樣品 B-大棚種植樣品 C-林下活樹附生種植樣品 D-巖壁附生種植樣品,下同

    3.2 PCA和代謝物積累模式分析

    代謝組數(shù)據(jù)經(jīng)標準化處理后進行PCA分析,前2個主成分特征值>1,模型累積解釋率R為0.64。一般認為R大于0.5,表示模型可靠,值越大越能反映樣本間的差異。PCA得分見圖2,表明各組樣本分布于不同區(qū)域,所建立方法可良好表征不同栽培模式下代謝物的差異性。變量載荷圖顯示幾種物質位于坐標邊緣(圖2-B),表明這些物質在組間樣本中差異明顯(有較大/較小的極端值等),它們對排序空間和樣本差異的貢獻較大,可能是重要的標志性代謝物。例如,矢車菊素-3--葡萄糖苷- 7,3′-di--(6′′--芥子酰葡萄糖苷)(Lmtp003544)、原花青素B1(mws0836)、金圣草黃素--葡萄糖醛酸--己糖苷(pmb0587)、金圣草黃素-C-己糖苷(pmb0689)、番石榴苷(mws4183)、多花紫藤苷(Lmdp003994)等。此外,聚類分析顯示(圖2-C),不同組間代謝物積累模式存在明顯的差異,同一組的樣本代謝物積累趨勢較為一致且高度聚合一支。

    A-樣品PCA得分圖 B-代謝物質載荷圖 C-代謝物聚類模式熱圖 S-巖壁種植樣品 T-林下活樹附生種植樣品 D-大棚種植樣品,表1同

    3.3 不同栽培模式下差異代謝物的鑒定

    OPLS-DA為分類監(jiān)督模式識別方法,可降低樣本組內(nèi)差異,能更加準確表征樣本組間特征。為進一步挖掘不同栽培模式下鐵皮石斛差異性,鑒定潛在的差異類黃酮物質,分別對3種栽培模式組間代謝數(shù)據(jù)進行OPLS-DA分析。

    根據(jù)OPLS-DA模型的變量重要性投影(variable importance in project,VIP)鑒定差異代謝物,標準為VIP>1且單維檢驗<0.05(差異倍數(shù)值fold change>2),共鑒定43個差異物質(表1)。表1分析發(fā)現(xiàn),大棚組與林下組共存在23個差異物質,其中7個物質在林下組上調高積累;大棚組與巖壁組也具23個差異物,其中21個在巖壁組上調積累;林下組與巖壁組則存在27個差異物,且均在巖壁組上調積累。值得注意是,3組間差異代謝物共同認定1個標志性物質—矢車菊素-3--葡萄糖苷- 7,3′-di--(6′′--芥子酰葡萄糖苷),暗示該物質在表征3種栽培模式下的代謝物中的貢獻率最大,與PCA結果相互佐證。進一步以VIP>3篩選發(fā)現(xiàn):大棚與林下組間,蘆丁、芹菜素-6,8--二葡萄糖苷及根皮素差異變化大;大棚與巖壁組間,柚皮素、紫鉚素及6-羥基山柰酚-7,6--二葡萄糖苷差異變化大;而在林下與巖壁組中,蘆丁、異金絲桃苷及6-羥基山柰酚-7--葡萄糖苷變化最大。綜合評價,上述物質對于表征不同栽培模式下鐵皮石斛代謝物的差異性貢獻較大。

    OPLS-DA模型優(yōu)化表明大棚組與林下組的2、2、2值,分別為0.665、0.955、0.802;大棚組與崖壁組2、2、2值,分別為0.644、0.979、0.902;林下組與巖壁組2、2、2,分別為0.789、0.973、0.932。OPLS-DA得分圖如圖3所示,表明所建立的模型對組間樣品具有良好的區(qū)分度,相同組內(nèi)樣品高度聚合一類。OPLS-DA排列驗證(=200)實驗表明,3組模型的原始RQ均大于Y置換后的相應的值,暗示3組模型未過度擬合,可用于后續(xù)代謝物的鑒定。

    3.4 差異代謝物代謝通路分析

    為深入解析差異代謝物參與的代謝路徑,對43個差異代謝物進行KEGG代謝途徑注釋(圖4),發(fā)現(xiàn)大部分化合物富集到ko00944(黃酮和黃酮醇生物合成)、ko00941(類黃酮生物物合成)2個代謝途徑中,1個化合物被注釋到ko00942(花青素生物合成)途徑。進一步分析發(fā)現(xiàn),與大棚種植模式相比,林下和巖壁組中上述代謝通路被明顯上調,次生代謝產(chǎn)物合成被顯著提升,尤其以巖壁種植為高。例如柚皮素、圣草酚、紫山葉素等黃酮類前體物質在林下組和巖壁組中含量增加明顯,分別較之大棚組提升了3.6倍和4.9倍、10.6倍和11.1倍、9.3倍和18.0倍?;ㄇ嗨卮x路徑中的矢車菊素-3--葡萄糖苷- 7,3′-di--(6′′--芥子酰葡萄糖苷)在林下組和巖壁組中的含量較之大棚組分別增加了6.8、26.6倍,暗示此物質很可能是近野生栽培過程中鐵皮石斛莖皮紫色表型形成的主要物質基礎。

    表1 不同栽培模式下鐵皮石斛類黃酮差異代謝物

    Table 1 Information of different flavonoes and flavonoids detected in D. officinale from different planting modes

    中文名英文名D vs TD vs ST vs S log2_FCVIPLog2_FCVIPLog2_FCVIP 4,2′,4′,6′-四羥基查耳酮4,2′,4′,6′-tetrahydroxychalcone1.57 2.95 1.52 2.35 ?? 柚皮素naringenin1.85 2.59 2.30 3.08 ?? 圣草酚eriodictyol3.41 2.61 3.48 2.21 ?? 槲皮素-O-蕓香苷-己糖quercetin-O-rutinoside-hexose??2.51 1.41 4.52 1.30 槲皮素-3-O-葡萄糖苷-7-O-鼠李糖苷quercetin-3-O-glucoside-7-O-rhamnoside?2.01 1.55 ??2.77 1.89 槲皮素-O-阿魏酰戊糖苷quercetin-O-feruloyl-pentoside?2.10 1.48 ??2.78 1.64 槲皮素-3-O-鼠李糖基半乳糖苷quercetin-3-O-rhanosylgalactoside?2.03 1.57 ??2.83 1.91 槲皮素-3-O-葡萄糖苷quercetin-3-O-glucoside??2.81 1.65 2.72 1.41 槲皮素-3-O-洋槐糖苷quercetin-3-O-robinobioside?1.73 1.64 ??2.59 2.17 槲皮素-O-葡萄糖苷quercetin-O-glucoside??2.49 1.55 3.16 1.42 槲皮素-3-O-β-D-葡萄糖苷quercetin 3-O-β-D-glucoside?2.70 2.61 ??3.41 2.91 甲基槲皮素葡萄糖-鼠李糖methylquercetin glu-rha??1.57 1.68 3.17 1.74 異金絲桃苷isohyperoside?2.28 2.93 ??2.88 3.23 金絲桃苷hyperin??2.54 2.01 3.19 1.87 橙皮苷hesperidin?1.73 2.83 ??2.70 3.73 橙皮素-5-O-葡萄糖苷hesperetin5-O-glucoside??2.85 1.82 2.78 1.56 橙皮素-O-丙二?;禾擒説esperetin O-malonylhexoside?1.44 1.22 ?2.25 1.24 ?? 繡線菊苷spiraeoside?? 2.64 1.60 2.61 1.39 異鼠李素-3-O-蕓香糖苷isorhamnetin-3-O-rutinoside ?? 1.64 2.13 2.90 2.14 芹菜素-8-C-α-L-阿拉伯糖苷apigenin 8-C-α-L-arabinopyranoside?? 3.76 1.01 ?? 芹菜素-6,8-C-二葡萄糖苷apigenin 6,8-C-diglucoside?1.18 4.24 ???? 矢車菊素-3-O-葡萄糖-7,3′-di-O-(6′′-O-芥子酰葡萄糖苷)cyanidin-3-O-glucoside-7,3′-di-O-[6′′-O-si-pylglucoside]4.74 2.79 2.76 1.24 1.98 1.98 紫杉葉素taxifolin??4.16 1.16 ?? 蘆丁rutin?1.62 4.68 ??2.42 6.08 短葉松素pinobanksin 1.96 2.61 2.31 2.97 ?? 香葉木素-7-O-半乳糖苷diosmetin-7-O-galactoside?1.39 1.15 ?1.89 1.33 ?? 根皮素phloretin3.03 3.19 2.79 2.66 ?? 金圣草黃素-6,8-二-C-葡萄糖苷chrysoeriol-6,8-di-C-glucoside?2.74 2.46 ??1.91 1.18 紫鉚素butin 1.89 2.76 2.37 3.32 ?? 木犀草素 8-C-己糖苷-O-己糖苷luteolin 8-C-hexosyl-O-hexoside?2.11 1.70 ??2.80 1.89 木犀草素-7-O-β-D-蕓香糖苷luteolin-7-O-β-D-rutinoside?1.27 2.29 ???? 木犀草素-6-C-2-葡萄糖醛酸葡萄糖苷luteolin-6-C-2-glucuronylglucoside?2.44 1.21 ???? 木犀草素-7-O-蕓香糖苷luteolin-7-O-rutinoside????4.24 1.29 6-羥基木犀草素 5-葡萄糖苷6-hydroxyluteolin 5-glucoside??2.91 1.65 ?? 牡荊素-2-O-D-吡喃葡萄糖苷vitexin-2-O-D-glucopyranoside?1.34 1.84 ???? 山柰酚-3-O-新橙皮糖苷kaempferol-3-O-neohesperidoside????4.02 2.11 山柰酚-3-O-洋槐糖苷kaempferol-3-O-robinobioside (biorobin)????4.06 1.63 山柰酚-3-O-蕓香糖苷kaempferol-3-O-rutinoside????3.87 1.60 山柰酚-3-O-β-D-(6′′-O-(E)-p-香豆酰)kaempferol-3-O-β-D-(6′′-O-(E)-p-coumaroyl) ????4.09 2.14 6-羥基山柰酚-7-O-葡萄糖苷6-hydroxykaempferol-7-O-glucoside????2.38 3.05 6-羥基山柰酚-3,6-O-二葡萄糖苷6-hydroxykaempferol-3,6-O-diglucoside??8.88 1.10 ?? 6-羥基山柰酚-7,6-O-二葡萄糖苷6-hydroxykaempferol-7,6-O-diglucoside??4.41 3.45 3.85 2.88 6-羥基山柰酚-3,7,6-O-三葡萄糖苷6-hydroxykaempferol-3,7,6-O-triglycoside??4.93 1.85 2.33 1.45

    A-林下附生種植相比大棚種植 B-巖壁附生種植相比大棚種植 C-巖壁附生種植林下附生種植

    紅顏色字體物質代表鑒定的差異物質;代謝物熱圖表示在不同栽培模式中的積累模式,藍色表示低水平,黃色表示高水平

    4 討論

    植物代謝組學具有“高維、海量”的特點,結合多元統(tǒng)計分析方法可多維度揭示植物體在響應不同環(huán)境下的代謝產(chǎn)物變化模式,從而解釋其生命本質問題[15]。UPLC-MS/MS廣泛靶向技術結合了非靶向和靶向代謝組學的優(yōu)點,基于代謝組自建、公共數(shù)據(jù)庫和多反應監(jiān)測模型,可定性、定量且高覆蓋地檢測樣品中代謝物的含量[16]。本實驗利用廣泛靶向組學技術,對來自3種栽培模式的12個鐵皮石斛樣品共檢測到202個黃酮類物質,其中大部分以黃酮(76個)、黃酮醇(60個)為主。不同栽培措施明顯改變了鐵皮石斛體內(nèi)黃酮類物質積累模式(圖3-C),林下和巖壁附生種植較大棚種植明顯提升了黃酮或類黃酮代謝通路(ko00944、ko00941)次生代謝產(chǎn)物合成水平(圖4)。比較發(fā)現(xiàn),柚皮素、圣草酚、紫山葉素等類黃酮合成通路中關鍵前體物質在林下和巖壁栽培組中均明顯上調積累,尤以巖壁附生栽培最高。多項研究表明,紫外線是促使植物合成黃酮類化合物的重要的脅迫環(huán)境之一[13]。例如,UV-A誘導大豆、及等基因的特異表達,且增強了大豆異黃酮、柚皮素等物質在大豆中的積累[17]。因而,分析認為林下活樹附生和巖壁附生等近野生栽培方式因無遮蔭措施,相比大棚栽培,其高UV強度很可能是造成鐵皮石斛高積累黃酮類物質重要因素。深入的原因,仍需通過關聯(lián)基因表達組數(shù)據(jù)予以驗證。

    此外,UV、低溫等逆境脅迫也是誘導植物花青素合成重要因子[18-19]。生產(chǎn)上,林下活樹附生和巖壁附生等近野生栽培方式往往促使鐵皮石斛莖皮形成紫色表型,且莖水溶液呈淡至深紫色,而大棚栽培鐵皮石斛水溶液往往呈深綠色。矢車菊素糖苷配體化合物是花青素中的主要成員,是植物細胞內(nèi)呈現(xiàn)紫色的主要色素之一。本研究中,廣泛靶向技術檢測到唯一個花青苷物質—矢車菊素-3--葡萄糖苷-7,3′- di--(6′′--芥子酰葡萄糖苷)(Lmtp003544),該物質在表征三種栽培模式下的代謝物差異性的貢獻率最大(圖2-A,表1)。比較發(fā)現(xiàn),相比大棚栽培的鐵皮石斛,近野生栽培方式Lmtp003544相對含量至少提升6.8倍,巖壁附生栽培甚至提升了26.6倍,暗示該物質極有可能是紫色表型形成的標志性成分。因而,未來以該指標作為不同來源鐵皮石斛的質量控制方法具有很好的可行性,將為藥材等級的劃分提供重要技術支撐。當前,近野生栽培的鐵皮石斛市場價格是大棚栽培的10~15倍。本研究首次在理論上揭示了近野生栽培能大大提升了類黃酮物質的積累。鑒于該類化合物的重要的醫(yī)藥保健價值,建議生產(chǎn)上大力推廣林下附生與巖壁附生栽培模式,以提高鐵皮石斛的品質和產(chǎn)業(yè)價值。

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

    [1] 中國藥典 [S]. 一部. 2020: 295.

    [2] 李以格, 楊杭, 姜琪夢, 等. 珍稀藥用植物鐵皮石斛的組學及功能基因研究進展 [J]. 生命科學, 2019, 31(9): 959-967.

    [3] 顏美秋, 陳素紅, 呂圭源. 石斛“厚腸胃”相關功效藥理學研究及應用進展 [J]. 中草藥, 2016, 47(21): 3918-3924.

    [4] 孫樂, 陳曉梅, 吳崇明, 等. 鐵皮石斛多糖藥理活性研究進展 [J]. 藥學學報, 2020, 55(10): 2322-2329.

    [5] 奚航獻, 劉晨, 劉京晶, 等. 鐵皮石斛化學成分、藥理作用及其質量標志物 (Q-marker) 的預測分析 [J]. 中草藥, 2020, 51(11): 3097-3109.

    [6] 李毓揚, 呂承豪, 吳廣, 等. 鐵皮石斛干預代謝綜合征的作用分子機制研究進展 [J]. 中國中藥雜志, 2019, 44(23): 5102-5108.

    [7] 陳曉梅, 田麗霞, 單婷婷, 等. 鐵皮石斛種質資源和遺傳育種研究進展 [J]. 藥學學報, 2018, 53(9): 1493-1503.

    [8] 斯金平, 陳梓云, 劉京晶, 等. 鐵皮石斛懸崖附生栽培技術研究 [J]. 中國中藥雜志, 2015, 40(12): 2289-2292.

    [9] 斯金平, 王琦, 劉仲健, 等. 鐵皮石斛產(chǎn)業(yè)化關鍵科學與技術的突破 [J]. 中國中藥雜志, 2017, 42(12): 2223-2227.

    [10] 王珂珂. 不同栽培模式下鐵皮石斛藥用價值的比較研究 [D]. 福州: 福建農(nóng)林大學, 2016.

    [11] Yuan Y D, Tang X G, Jia Z H,. The effects of ecological factors on the main medicinal components ofunder different cultivation modes [J]., 2020, 11(1): 94.

    [12] 呂朝耕, 李佳興, 蔣靖怡, 等. 浙產(chǎn)鐵皮石斛不同生長期石斛酚、有機酸和黃酮類成分積累變化 [J]. 中華中醫(yī)藥雜志, 2019, 34(5): 1989-1992.

    [13] 潘俊倩, 佟曦然, 郭寶林. 光對植物黃酮類化合物的影響研究進展 [J]. 中國中藥雜志, 2016, 41(21): 3897-3903.

    [14] Chen W, Gong L, Guo Z L,. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: Application in the study of rice metabolomics [J]., 2013, 6(6): 1769-1780.

    [15] Villate A, San Nicolas M, Gallastegi M,. Review: Metabolomics as a prediction tool for plants performance under environmental stress [J]., 2021, 303: 110789.

    [16] Wei R, Li G D, Seymour A B. High-throughput and multiplexed LC/MS/MRM method for targeted metabolomics [J]., 2010, 82(13): 5527-5533.

    [17] Lim Y J, Lyu J I, Kwon S J,. Effects of UV-A radiation on organ-specific accumulation and gene expression of isoflavones and flavonols in soybean sprout [J]., 2021, 339: 128080.

    [18] Landi M, Tattini M, Gould K S. Multiple functional roles of anthocyanins in plant-environment interactions [J]., 2015, 119: 4-17.

    [19] 王峰, 王秀杰, 趙勝男, 等. 光對園藝植物花青素生物合成的調控作用 [J]. 中國農(nóng)業(yè)科學, 2020, 53(23): 4904-4917.

    UPLC-MS/MS metabonomics technology revealing difference of flavonoids inunder different planting modes

    LI Fan1, 2, JIN Chuan-gao3, REN Xian-ying2, 4, FEI Xuan2, ZHAO Qi2, 4, PAN Xiao-jun2, WU Zhi-gang2, JIANG Cheng-xi1, 4

    1. Institute of Life Sciences, Wenzhou University, Wenzhou 325035, China 2. School of Pharmacy, Wenzhou Medical University, Wenzhou 325035, China 3. Zhejiang Gaobizi Biological Technology Co., Ltd., Leqing 325615, China 4. Jiuhuashan Buddhist Medicine Research Institute, Chizhou 247100, China

    To reveal the difference of flavonoid components inunder different planting modes, in order to provide a basis for the innovation of planting mode.The extensive target metabonomics tool via UPLC-MS/MS workflow was used to detect the flavonoids instem samples produced from three planting modes (including greenhouse, living tree epiphysis planting, and cliff epiphytic planting). Principal component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA) were performed to identify differential metabolites.A total of 202 compounds were identified, mainly including flavonoids (76) and flavonols (60). A novel anthocyanin, cyanidin-3-- glucoside-7,3′-di--(6′′--Si-pylglucoside) was first detected instem samples. OPLS-DA identified 43 differential compounds, among which 23 flavonoids were found between greenhouse planting and living tree epiphysis planting mode, as well as greenhouse planting and cliff epiphytic planting mode. Additionally, there were 27 differential flavonoids between living tree epiphysis planting and cliff epiphytic planting modes. Overall, the difference of samples from different plant modes was well explained by these differential flavonoids.The living tree epiphysis and cliff epiphytic planting modes significantly unregulated the flavonoids biosynthetic pathway relative to greenhouse planting, thus improving the quality of medicinal materials. This result provides a scientific basis for the promotion of novel planting modes withinproduction.

    Kimura et Migo.;UPLC-MS/MS; flavonoids; metabonomics; planting mode

    R286.2

    A

    0253 - 2670(2022)04 - 1156 - 07

    10.7501/j.issn.0253-2670.2022.04.024

    2021-09-06

    溫州醫(yī)科大學人才項目(89220027);浙江省重點實驗室項目(2021E10013);安徽濟人藥業(yè)有限公司課題(KJHX2009);合肥市未來藥物開發(fā)有限公司(KJHX2008);云南大理藥業(yè)股份有限公司課題(KJHX1603)

    李 帆,碩士研究生,主要從事中藥生物技術研究。E-mail: 1471750585@qq.com

    姜程曦,教授,主要從事中藥生物技術研究。E-mail: jiangchengxi@126.com

    吳志剛,副研究員,主要從事藥用植物次生代謝產(chǎn)物研究。E-mail: wuzhigang177@126.com

    #共同第一作者:金傳高,工程師,從事鐵皮石斛技術開發(fā)研究。E-mail: 850254737@qq.com

    [責任編輯 時圣明]

    日韩三级伦理在线观看| 国产深夜福利视频在线观看| 秋霞伦理黄片| 亚洲av日韩在线播放| 国产成人91sexporn| 黄色怎么调成土黄色| 少妇熟女欧美另类| 久热这里只有精品99| 丝袜美腿诱惑在线| 一级毛片 在线播放| 国产免费视频播放在线视频| 美女视频免费永久观看网站| 老鸭窝网址在线观看| 日产精品乱码卡一卡2卡三| 日本猛色少妇xxxxx猛交久久| 春色校园在线视频观看| 免费少妇av软件| 成人亚洲欧美一区二区av| 婷婷色综合大香蕉| 搡老乐熟女国产| 欧美精品亚洲一区二区| 少妇被粗大的猛进出69影院| 国产黄色免费在线视频| 久久ye,这里只有精品| 欧美精品亚洲一区二区| 久久这里有精品视频免费| 免费观看av网站的网址| 国产成人精品婷婷| 国产精品熟女久久久久浪| 波多野结衣av一区二区av| 成人亚洲欧美一区二区av| 国产精品秋霞免费鲁丝片| 久久综合国产亚洲精品| 啦啦啦啦在线视频资源| www.av在线官网国产| 香蕉国产在线看| 亚洲国产精品一区三区| 亚洲精品在线美女| 中文字幕人妻丝袜一区二区 | 日产精品乱码卡一卡2卡三| 国产亚洲欧美精品永久| 狂野欧美激情性bbbbbb| 国产又色又爽无遮挡免| 久久精品亚洲av国产电影网| 国产av国产精品国产| 丝袜在线中文字幕| 久久久国产精品麻豆| 日本午夜av视频| 人人澡人人妻人| 一区二区三区激情视频| 欧美日本中文国产一区发布| 精品国产一区二区久久| 超碰成人久久| 午夜免费观看性视频| 精品福利永久在线观看| av天堂久久9| 少妇熟女欧美另类| 日本黄色日本黄色录像| 69精品国产乱码久久久| 性色avwww在线观看| 伊人久久大香线蕉亚洲五| 80岁老熟妇乱子伦牲交| 欧美最新免费一区二区三区| 91国产中文字幕| 国产日韩一区二区三区精品不卡| 国产无遮挡羞羞视频在线观看| 伦精品一区二区三区| 国产精品无大码| 黄色怎么调成土黄色| 飞空精品影院首页| 久久久久精品人妻al黑| 日本vs欧美在线观看视频| 好男人视频免费观看在线| 女人高潮潮喷娇喘18禁视频| 亚洲精品aⅴ在线观看| 国产男女内射视频| 妹子高潮喷水视频| 夫妻性生交免费视频一级片| 18禁裸乳无遮挡动漫免费视频| 99久久精品国产国产毛片| 天天躁狠狠躁夜夜躁狠狠躁| 久久久亚洲精品成人影院| 男女国产视频网站| 国产欧美日韩一区二区三区在线| av.在线天堂| 日本色播在线视频| 最新中文字幕久久久久| 丝袜喷水一区| tube8黄色片| 男女无遮挡免费网站观看| 成人毛片a级毛片在线播放| 国产男女内射视频| 建设人人有责人人尽责人人享有的| 亚洲精品久久成人aⅴ小说| 美女大奶头黄色视频| 一级毛片电影观看| 欧美成人精品欧美一级黄| 五月开心婷婷网| 亚洲伊人久久精品综合| 国产精品久久久久久精品电影小说| av视频免费观看在线观看| 日韩av不卡免费在线播放| 欧美国产精品一级二级三级| 欧美成人午夜精品| 欧美av亚洲av综合av国产av | 亚洲人成网站在线观看播放| 91aial.com中文字幕在线观看| 亚洲 欧美一区二区三区| 自线自在国产av| 国产精品.久久久| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕精品免费在线观看视频| 久久久a久久爽久久v久久| 久久久久精品人妻al黑| 日韩人妻精品一区2区三区| 2018国产大陆天天弄谢| 日韩免费高清中文字幕av| 欧美成人精品欧美一级黄| 考比视频在线观看| 免费高清在线观看日韩| 男女免费视频国产| 成人国产麻豆网| 伦理电影免费视频| 一本大道久久a久久精品| 久热久热在线精品观看| 精品一品国产午夜福利视频| www.av在线官网国产| 亚洲精品视频女| 日韩三级伦理在线观看| 国产1区2区3区精品| 激情视频va一区二区三区| 免费黄频网站在线观看国产| 久久精品人人爽人人爽视色| 欧美日韩视频精品一区| 国产精品免费大片| 亚洲综合精品二区| 国产不卡av网站在线观看| 免费观看a级毛片全部| 精品少妇内射三级| 精品人妻在线不人妻| 人人妻人人澡人人爽人人夜夜| 日韩一卡2卡3卡4卡2021年| 中文字幕精品免费在线观看视频| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 亚洲av免费高清在线观看| 亚洲欧美色中文字幕在线| a级毛片在线看网站| 久久ye,这里只有精品| 女人被躁到高潮嗷嗷叫费观| 国产精品av久久久久免费| 大片电影免费在线观看免费| 亚洲成人一二三区av| 亚洲av电影在线进入| 自拍欧美九色日韩亚洲蝌蚪91| 国产乱来视频区| www.av在线官网国产| 各种免费的搞黄视频| 久久精品夜色国产| 欧美最新免费一区二区三区| 国产日韩一区二区三区精品不卡| 久久精品国产鲁丝片午夜精品| 国产福利在线免费观看视频| 一二三四中文在线观看免费高清| 精品人妻在线不人妻| 69精品国产乱码久久久| 又黄又粗又硬又大视频| 捣出白浆h1v1| 精品久久蜜臀av无| 久久久久久久大尺度免费视频| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 韩国精品一区二区三区| 伊人亚洲综合成人网| 黄色怎么调成土黄色| 久久久a久久爽久久v久久| 丝袜人妻中文字幕| 美国免费a级毛片| 国产精品国产av在线观看| 日韩伦理黄色片| 国产精品熟女久久久久浪| 人妻一区二区av| 日日撸夜夜添| 亚洲成人一二三区av| 国产av码专区亚洲av| 黑人猛操日本美女一级片| 飞空精品影院首页| 一级,二级,三级黄色视频| 亚洲色图综合在线观看| 少妇的逼水好多| 精品国产乱码久久久久久小说| 不卡av一区二区三区| 免费在线观看完整版高清| 美女国产视频在线观看| 999精品在线视频| 中文字幕人妻丝袜制服| 最近中文字幕高清免费大全6| 丝袜喷水一区| 国产极品粉嫩免费观看在线| 麻豆乱淫一区二区| 国产免费福利视频在线观看| 最近中文字幕高清免费大全6| 一级爰片在线观看| 亚洲美女黄色视频免费看| 高清欧美精品videossex| 亚洲第一青青草原| 午夜激情久久久久久久| 久久精品国产鲁丝片午夜精品| 亚洲精品自拍成人| 老司机亚洲免费影院| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站 | 老女人水多毛片| 国产在线一区二区三区精| 不卡视频在线观看欧美| 视频在线观看一区二区三区| 精品一区在线观看国产| 成人二区视频| 99国产精品免费福利视频| 中文欧美无线码| 色视频在线一区二区三区| 中文字幕亚洲精品专区| www.精华液| 狂野欧美激情性bbbbbb| 岛国毛片在线播放| 久久久久久伊人网av| 国产黄色视频一区二区在线观看| 欧美日本中文国产一区发布| 成人国产麻豆网| 久久韩国三级中文字幕| 精品亚洲乱码少妇综合久久| 免费高清在线观看日韩| 十分钟在线观看高清视频www| 欧美国产精品va在线观看不卡| kizo精华| 日韩中字成人| 我要看黄色一级片免费的| 精品国产一区二区三区四区第35| 男女啪啪激烈高潮av片| 国产av国产精品国产| 欧美97在线视频| 午夜av观看不卡| 中文字幕精品免费在线观看视频| 蜜桃在线观看..| 久久久精品国产亚洲av高清涩受| 春色校园在线视频观看| 午夜福利乱码中文字幕| 日产精品乱码卡一卡2卡三| 免费观看性生交大片5| 18禁观看日本| 黄频高清免费视频| 久久影院123| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 国产一区二区在线观看av| 国产精品免费视频内射| 日韩一区二区视频免费看| 日韩人妻精品一区2区三区| 色哟哟·www| 国产精品国产三级国产专区5o| 国产亚洲精品第一综合不卡| 久久久久久久久久久久大奶| 久久青草综合色| 巨乳人妻的诱惑在线观看| 午夜日本视频在线| 亚洲欧洲精品一区二区精品久久久 | 黄频高清免费视频| √禁漫天堂资源中文www| 久久人人97超碰香蕉20202| 亚洲国产精品一区二区三区在线| 欧美国产精品一级二级三级| 久久久久久久国产电影| 午夜福利网站1000一区二区三区| 男女边吃奶边做爰视频| 亚洲精品第二区| 精品酒店卫生间| 国产精品免费大片| 亚洲五月色婷婷综合| 老汉色av国产亚洲站长工具| 国产成人精品婷婷| 激情视频va一区二区三区| 看免费成人av毛片| 日本猛色少妇xxxxx猛交久久| 韩国高清视频一区二区三区| 在线天堂中文资源库| 免费播放大片免费观看视频在线观看| 黄色视频在线播放观看不卡| 久久久久久久久久久免费av| 精品一区在线观看国产| 国产野战对白在线观看| 男女边摸边吃奶| 亚洲国产精品一区三区| 乱人伦中国视频| 国产成人精品无人区| 日韩av免费高清视频| 亚洲精品日本国产第一区| 精品人妻在线不人妻| 国产精品一国产av| 色播在线永久视频| 国产成人精品一,二区| 精品人妻偷拍中文字幕| 人人澡人人妻人| 飞空精品影院首页| 国产av精品麻豆| www.精华液| 国产成人欧美| 久久精品人人爽人人爽视色| 精品午夜福利在线看| 日韩中文字幕欧美一区二区 | 久久久久网色| 赤兔流量卡办理| 午夜日本视频在线| 亚洲欧美一区二区三区久久| 久久久国产精品麻豆| 免费观看性生交大片5| 欧美黄色片欧美黄色片| 欧美国产精品一级二级三级| 考比视频在线观看| 日韩大片免费观看网站| 欧美成人精品欧美一级黄| 久久精品国产亚洲av高清一级| 少妇的丰满在线观看| 如日韩欧美国产精品一区二区三区| 亚洲精品久久午夜乱码| av在线app专区| 黄片无遮挡物在线观看| 欧美人与性动交α欧美软件| 在线天堂中文资源库| 黑人猛操日本美女一级片| 极品少妇高潮喷水抽搐| 免费人妻精品一区二区三区视频| 精品卡一卡二卡四卡免费| 欧美变态另类bdsm刘玥| 国产成人91sexporn| 欧美日韩精品网址| 少妇的逼水好多| 欧美日韩精品成人综合77777| 高清欧美精品videossex| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 又大又黄又爽视频免费| 亚洲国产最新在线播放| 成人午夜精彩视频在线观看| 哪个播放器可以免费观看大片| 亚洲综合精品二区| 国产毛片在线视频| 国产97色在线日韩免费| 国产精品国产三级国产专区5o| 夫妻性生交免费视频一级片| 一级毛片电影观看| av线在线观看网站| 亚洲av日韩在线播放| 啦啦啦中文免费视频观看日本| 天天影视国产精品| 黄片无遮挡物在线观看| 老司机亚洲免费影院| 十八禁网站网址无遮挡| 国产成人精品婷婷| 女的被弄到高潮叫床怎么办| 国产成人精品久久二区二区91 | 人妻人人澡人人爽人人| 欧美最新免费一区二区三区| 婷婷色综合www| 国产白丝娇喘喷水9色精品| 国产成人免费无遮挡视频| 国产在线免费精品| 午夜91福利影院| 90打野战视频偷拍视频| 精品卡一卡二卡四卡免费| 成人免费观看视频高清| 国产精品久久久久久av不卡| 国产精品.久久久| www.熟女人妻精品国产| 熟女av电影| 久久热在线av| 久久鲁丝午夜福利片| 中文字幕av电影在线播放| 伦理电影大哥的女人| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久久免费av| 亚洲精品国产av成人精品| 你懂的网址亚洲精品在线观看| 我要看黄色一级片免费的| 欧美xxⅹ黑人| 一二三四在线观看免费中文在| 亚洲精品自拍成人| 韩国高清视频一区二区三区| 制服诱惑二区| 久久精品国产亚洲av天美| 又大又黄又爽视频免费| 精品亚洲乱码少妇综合久久| 777米奇影视久久| 久久婷婷青草| 成人影院久久| 老汉色av国产亚洲站长工具| 国产1区2区3区精品| 女人久久www免费人成看片| 高清不卡的av网站| 成人亚洲欧美一区二区av| 婷婷色综合www| 寂寞人妻少妇视频99o| 男人爽女人下面视频在线观看| 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 日韩av在线免费看完整版不卡| 少妇人妻精品综合一区二区| 日本av免费视频播放| 久久久久久人人人人人| 一级毛片 在线播放| 欧美av亚洲av综合av国产av | 黑人欧美特级aaaaaa片| 女人久久www免费人成看片| 精品国产国语对白av| 在线天堂最新版资源| 热99国产精品久久久久久7| 中文天堂在线官网| 精品国产超薄肉色丝袜足j| 丝袜美足系列| 欧美日韩综合久久久久久| 免费在线观看视频国产中文字幕亚洲 | 黄色怎么调成土黄色| 五月开心婷婷网| 性色avwww在线观看| 久久精品国产亚洲av高清一级| 国产 精品1| 2022亚洲国产成人精品| 精品99又大又爽又粗少妇毛片| 久久久久久伊人网av| 久久国产精品男人的天堂亚洲| 午夜激情久久久久久久| 亚洲国产精品999| 亚洲色图综合在线观看| 人体艺术视频欧美日本| 777久久人妻少妇嫩草av网站| 国产极品天堂在线| 飞空精品影院首页| 国产97色在线日韩免费| 国产日韩欧美在线精品| 国产亚洲最大av| 99国产精品免费福利视频| 久久 成人 亚洲| 人体艺术视频欧美日本| 亚洲欧美一区二区三区国产| 一区在线观看完整版| 十八禁高潮呻吟视频| 777米奇影视久久| 热99久久久久精品小说推荐| 亚洲精品av麻豆狂野| 99国产综合亚洲精品| 在线 av 中文字幕| 成人手机av| 三上悠亚av全集在线观看| 99久久综合免费| 十八禁网站网址无遮挡| 午夜91福利影院| 亚洲精品aⅴ在线观看| 建设人人有责人人尽责人人享有的| 观看av在线不卡| 涩涩av久久男人的天堂| www.精华液| 亚洲久久久国产精品| 亚洲国产欧美网| 满18在线观看网站| 三级国产精品片| 欧美日韩一级在线毛片| 免费少妇av软件| 制服人妻中文乱码| 国产精品.久久久| 热99国产精品久久久久久7| 视频在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 一级片'在线观看视频| 久久av网站| 午夜日韩欧美国产| 中文天堂在线官网| 国产精品蜜桃在线观看| 一本—道久久a久久精品蜜桃钙片| 久热这里只有精品99| 久久青草综合色| 丝袜美腿诱惑在线| 99热全是精品| 亚洲欧洲国产日韩| videos熟女内射| 桃花免费在线播放| 黑人巨大精品欧美一区二区蜜桃| 国产熟女午夜一区二区三区| 亚洲av综合色区一区| 久久这里只有精品19| 国产精品香港三级国产av潘金莲 | 日本猛色少妇xxxxx猛交久久| 高清视频免费观看一区二区| 丰满少妇做爰视频| 婷婷色麻豆天堂久久| 男人舔女人的私密视频| 有码 亚洲区| 欧美日本中文国产一区发布| 在现免费观看毛片| 男人操女人黄网站| tube8黄色片| 丰满迷人的少妇在线观看| 国产男女超爽视频在线观看| 777久久人妻少妇嫩草av网站| 欧美日韩精品成人综合77777| 午夜影院在线不卡| 综合色丁香网| 欧美xxⅹ黑人| 成人二区视频| 777米奇影视久久| 亚洲久久久国产精品| 亚洲欧美清纯卡通| 久久99精品国语久久久| 涩涩av久久男人的天堂| 亚洲av电影在线观看一区二区三区| 午夜福利在线免费观看网站| 亚洲精品第二区| 日本wwww免费看| 成年动漫av网址| 高清不卡的av网站| 母亲3免费完整高清在线观看 | 亚洲精品久久成人aⅴ小说| 天堂8中文在线网| 国产熟女欧美一区二区| 天堂俺去俺来也www色官网| 日日撸夜夜添| 久久av网站| 亚洲欧美清纯卡通| 国产精品无大码| 国产亚洲av片在线观看秒播厂| 大码成人一级视频| 久久婷婷青草| 久久久久人妻精品一区果冻| 99re6热这里在线精品视频| 午夜福利影视在线免费观看| 91国产中文字幕| 乱人伦中国视频| 中文字幕最新亚洲高清| 色网站视频免费| 97人妻天天添夜夜摸| 国产国语露脸激情在线看| 精品一品国产午夜福利视频| 久久韩国三级中文字幕| 麻豆av在线久日| 九九爱精品视频在线观看| 国产男人的电影天堂91| 国产成人一区二区在线| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| 日日啪夜夜爽| 成人18禁高潮啪啪吃奶动态图| 国产熟女欧美一区二区| 蜜桃国产av成人99| 欧美精品国产亚洲| 亚洲久久久国产精品| 不卡av一区二区三区| 婷婷色综合大香蕉| 日韩欧美一区视频在线观看| 毛片一级片免费看久久久久| 汤姆久久久久久久影院中文字幕| av在线老鸭窝| 久久99热这里只频精品6学生| 制服丝袜香蕉在线| 午夜精品国产一区二区电影| 男女边吃奶边做爰视频| 亚洲成人av在线免费| 热re99久久精品国产66热6| 亚洲国产欧美在线一区| 伦理电影大哥的女人| 侵犯人妻中文字幕一二三四区| 99热网站在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人一区二区在线| 在现免费观看毛片| 日韩精品有码人妻一区| 日韩在线高清观看一区二区三区| 人妻人人澡人人爽人人| 亚洲天堂av无毛| 日韩一区二区视频免费看| 9191精品国产免费久久| 一区二区三区精品91| 免费少妇av软件| 国产1区2区3区精品| a级毛片在线看网站| 99久久人妻综合| 久久久久久久久免费视频了| 午夜福利网站1000一区二区三区| 国产深夜福利视频在线观看| 黑人欧美特级aaaaaa片| 亚洲av日韩在线播放| 精品99又大又爽又粗少妇毛片| 亚洲欧洲精品一区二区精品久久久 | 99re6热这里在线精品视频| av女优亚洲男人天堂| 亚洲国产av影院在线观看| 久久亚洲国产成人精品v| www日本在线高清视频| 26uuu在线亚洲综合色| 日韩在线高清观看一区二区三区| 国产精品熟女久久久久浪| 人人妻人人澡人人看| 国产伦理片在线播放av一区| 亚洲一区中文字幕在线| 日韩一卡2卡3卡4卡2021年| 在线精品无人区一区二区三| 777久久人妻少妇嫩草av网站| 色婷婷av一区二区三区视频| 十八禁网站网址无遮挡| 一级a爱视频在线免费观看| 最近最新中文字幕免费大全7| 免费大片黄手机在线观看| 韩国av在线不卡| 精品久久蜜臀av无| 国产免费现黄频在线看| 婷婷色av中文字幕| 欧美少妇被猛烈插入视频|