亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于UPLC-MS/MS的林下三七不同部位生物堿類成分分析

        2023-07-30 14:23:25張澳嚴星茹熊冰杰黃佑國何舒周發(fā)金王澍施蕊
        南方農(nóng)業(yè)·上旬 2023年5期
        關鍵詞:代謝組學

        張澳 嚴星茹 熊冰杰 黃佑國 何舒 周發(fā)金 王澍 施蕊

        摘 要 為了挖掘林下三七地上部分的利用價值,給提高三七資源的綜合利用率提供新思路,采用超高效液相色譜和串聯(lián)質(zhì)譜法分析了林下三七的果實、葉片、根三個部位的生物堿類成分。結果表明,林下三七果實、葉片和根共檢測到53個生物堿類化合物,但各部位總含量不相同,果實>葉片>根;多元統(tǒng)計分析發(fā)現(xiàn)果實與根、葉片與根成分差異顯著,分別有26個和24個差異化合物,果實中3-吲哚丙酸、甲氧基吲哚乙酸和油酰單乙醇胺的含量較高,葉片中3-羥基吡啶、N-甲基煙酰胺和5-羥基吲哚-3-乙醇的含量較高;篩選得到鑒別果實、葉片、根的生物堿類差異標志物,分別是油酰單乙醇胺、3-羥基吡啶和乙酰色胺。

        關鍵詞 林下三七;生物堿類成分;代謝組學;差異化合物;植株部位

        中圖分類號:S567.23+6 文獻標志碼:A DOI:10.19415/j.cnki.1673-890x.2023.09.004

        五加科人參屬三七[Panax notoginseng (Burk) F. H. Chen]作為云藥的代表聞名杏林,《本草綱目》譽之“金不換”“補血第一”,即具備散瘀止血、消腫止痛的功效[1]。近年來云南三七產(chǎn)量全國占比高達98%以上,正逐漸形成種植加工銷售一體化的龐大產(chǎn)業(yè)體系[2]。國家中醫(yī)藥發(fā)展戰(zhàn)略提出需要全面提升中藥產(chǎn)業(yè)發(fā)展水平即實現(xiàn)產(chǎn)業(yè)轉型升級[3],而三七產(chǎn)業(yè)轉型升級過程發(fā)現(xiàn)了一個嚴重的門檻,三七產(chǎn)品開發(fā)深度和廣度不足導致三七產(chǎn)值低、產(chǎn)業(yè)鏈短[4-6]。因此,為增加三七產(chǎn)值和延長三七產(chǎn)業(yè)鏈,提高三七資源的綜合利用率是十分必要的。

        現(xiàn)階段,三七的研究領域主要集中于三七的傳統(tǒng)入藥部位根部及其皂苷類成分[7-8],有關三七不同部位和非皂苷類成分的研究較少。三七的果實和葉片作為其生長發(fā)育過程中的伴生產(chǎn)物,常常順應自然凋落腐爛,實際是一種資源浪費。生物堿類成分是植物體內(nèi)的一類次生代謝產(chǎn)物,根據(jù)化學結構可分為吡咯、吡啶、吲哚等十余種[9],具備豐富的生物活性,王騰飛等研究發(fā)現(xiàn)其具有抗癌、抗病毒、抗炎和抗菌等多種藥理作用,在植物的適應性、抗逆性等方面也發(fā)揮著重要作用[10]。因此,為增加三七資源綜合利用的多樣性,可以從三七不同部位的生物堿類成分入手探尋新思路。

        植物代謝組學,是一門采用液相-質(zhì)譜聯(lián)用、氣相-質(zhì)譜聯(lián)用等多種手段全面、高通量對植物化合物進行定量、定性分析的新興學科[11]。代謝組學方法正逐漸成為應用在藥用植物的中藥基原鑒別[12]、野生與栽培品的鑒別[13]、藥用部位鑒別[14-15]、不同發(fā)育期的成分鑒定及不同產(chǎn)地藥材的品質(zhì)差異[16]等方面的有效手段。

        本實驗采用廣泛靶向代謝組學方法對三七的果實、葉片、根的生物堿類成分進行分析,

        通過主成分分析(PCA)結合正交偏最小二乘判別分析(OPLS-DA)的方法,篩選三七不同部位的差異化合物,探尋鑒別三七部位的生物堿類標志物,挖掘三七果實、葉片和根部的利用價值,以期為提高三七資源的綜合利用率提供新思路。

        1 ?材料與方法

        1.1 ?實驗材料

        本實驗樣品林下三七采集于瀾滄竹塘(N22.74°,E99.82°,H1 457.39 m),按不同部位分為果實、葉片、根,即三組樣品,每組3次生物學重復,樣品信息如表1所示。

        儀器:超高效液相色譜(Ultra Performance Liquid Chromatography,UPLC)和串聯(lián)質(zhì)譜儀(Tandem mass spectrometry,MS/MS);凍干機(Scientz-100F);研磨儀(MM 400, Retsch);離心機;渦旋儀,精密天平。

        試劑:甲醇(色譜純,Merck);乙腈(色譜純,Merck);標準品(BioBioPha/Sigma-Aldrich);純水。

        1.2 ?實驗方法

        1.2.1 ?樣品制備和提取

        將林下三七完整植株鮮品真空冷凍干燥成干品,然后分揀成果實、葉片、根三部分,分部位放置研磨儀中研磨成粉末(30 Hz,1.5 min),分裝編號成F(果實)、L(葉片)、R(根)備用。精密稱定100 mg粉末浸泡于1.2 mL的70%甲醇提取液中,使用借助渦旋儀每30 min渦旋1次,每次30 s,共渦旋6次,然后充分混勻后,放置于4 ℃冰箱冷凍過夜?;旌弦河秒x心機離心(12 000 r·min-1)10 min,吸取上清液,用微孔濾膜(0.22 μm)過濾得到濾液,裝入進樣瓶進樣分析。

        1.2.2 ?超高效液相色譜和串聯(lián)質(zhì)譜法(UPLC-MS/MS)檢測

        液相條件:采用Agilent SB-C18 1.8?m(2.1 mm×100 mm)的色譜柱進行色譜分離;設置柱溫40 ℃,流速0.35 mL·min-1,進樣量4 μL;流動相A為加入0.1%甲酸的超純水,B相為加入0.1%甲酸的乙腈,梯度洗脫條件見表2。

        質(zhì)譜條件:采用三重四極桿線性離子阱質(zhì)譜儀(Q TRAP)進行質(zhì)譜分析,配備ESI Turbo離子噴霧接口;渦輪噴霧離子源操作參數(shù)設置源溫度550 °C,離子噴霧電壓(IS)5 500 V(正離子模式)/-4 500 V(負離子模式),離子源氣體I(GSI)、氣體II(GSII)和簾氣(CUR)分別設置為50、60和25.0 psi,碰撞誘導電離參數(shù)高。

        1.2.3 ?數(shù)據(jù)處理

        采用UPLC-MS/MS方法獲得質(zhì)譜峰圖,利用Analyst 1.6.3軟件處理三七樣品質(zhì)譜數(shù)據(jù),利用MultiaQuant軟件對質(zhì)譜峰進行積分校正,完成定性定量處理。對通過混合樣本檢測得到的總離子流圖(TIC圖)進行重疊展示分析,判斷數(shù)據(jù)的重復性、可靠性和儀器的穩(wěn)定性。處理后的數(shù)據(jù)歸一化后采用多元統(tǒng)計分析中的主成分分析方法,對龐大復雜的數(shù)據(jù)進行“簡化和降維”,建立PCA得分圖,根據(jù)各樣本組間、組內(nèi)的分離趨勢判斷其差異程度。采用正交偏最小二乘法判別分析(OPLS-DA)彌補PCA方法的不足,在原始數(shù)據(jù)進行l(wèi)og2轉換后進行中心化處理,利用R軟件中的MetaboAnalystR包OPLSR.Anal函數(shù)進行分析,建立OPLS-DA正、負模型,得到OPLS-DA得分圖、OPLS-DA S-plo圖、VIP(投影中變量重要性)圖和聚類熱圖,使組間區(qū)分最大化,篩選出差異化合物。結合PCA分析和OPLS-DA分析,采用VIP值和FC值(差異倍數(shù)值)組合的方法判斷三七不同部位之間的差異情況,得到核心差異化合物和明確三七不同部位的開發(fā)利用價值[17]。

        2 ?結果與分析

        2.1 ?生物堿類化合物的定性定量分析

        本實驗采用三七的三個部位即果實、葉片、根,共選取9個樣品進行代謝組學分析。

        通過UPLC-MS/MS方法檢測所得MRM化合物檢測多峰圖(見圖1),峰積分校正處理后經(jīng)標準品對照和數(shù)據(jù)庫自查,結果表明一共檢測到53個生物堿類化合物,經(jīng)含量對比發(fā)現(xiàn)果實、葉片和根的生物堿總含量具有差異性,果實>葉片>根。后續(xù)通過多元統(tǒng)計分析方法對53個生物堿類化合物進行具體分析,找出果實、葉片和根的差異化合物。

        2.2 ?主成分分析

        對檢測所得的化合物進行主成分分析,建立PCA得分圖,每個點表示一個樣品,同樣的圖形表示同一組樣品,點與點之間的距離遠近表示分離效果,距離越遠分離效果越好,表示樣品間差異性越大。

        如圖2所示,第一主成分的解釋率為42.51%,第二主成分的解釋率為24.51%,表明兩個主成分能夠基本反映檢測樣的主要特征信息;mix樣本即質(zhì)控樣本,位于三個部位樣本的中間且組內(nèi)趨于同一點,說明數(shù)據(jù)重復性良好、可靠;F位于PC1的負半軸和PC2的正半軸,L位于PC1的正半軸和PC2的負半軸,R位于PC1的正半軸和PC2的正半軸,F(xiàn)、L、R三者間明顯分離,說明果實、葉片和根之間的生物堿類成分差異明顯。

        2.3 ?OPLS-DA分析

        正交偏最小二乘法判別分析方法可以與PCA方法互補,更全面地計算成分間的相關性,去除不相關的差異信息,使組間區(qū)分最大化篩選差異化合物。對三個部位作兩兩對比即FvsR(果實和根對比)、LvsR(葉片和根對比)、FvsL(果實和葉片對比),進行OPLS-DA分析區(qū)分不同部位的組間差異,結果再一次論證果實、葉片和根之間的生物堿類成分差異明顯,與PCA分析所得結果一致。其中FvsL差異性比FvsR和LvsR較低,加之與傳統(tǒng)藥用部位根作對比更有意義,所以繼而將重點分析FvsR和LvsR的差異化合物。

        以VIP ≥ 1,F(xiàn)C(差異倍數(shù)值) ≥2和FC ≤ 0.5作為篩選標準明確不同部位間差異顯著的化合物,F(xiàn)vsR篩選得到26個差異化合物,下調(diào)22個,上調(diào)4個;LvsR篩選得到24個差異化合物,下調(diào)15個,上調(diào)9個。如圖3所示,正值(或空心圓)代表上調(diào)差異表達化合物,負值(或實心圓)代表下調(diào)差異表達化合物,上調(diào)表示該化合物在根中的相對含量較高,下調(diào)表示該化合物在根中的相對含量較低。果實和根對比中發(fā)現(xiàn),上調(diào)最明顯的物質(zhì)是乙酰色胺(pmb0131),下調(diào)最明顯的物質(zhì)是油酰單乙醇胺(mws0983);葉片和根對比中發(fā)現(xiàn),上調(diào)最明顯的物質(zhì)是乙酰色胺(pmb0131),下調(diào)最明顯的物質(zhì)是N-甲基煙酰胺(pma6298)。

        2.4 ?差異化合物分析

        多元統(tǒng)計分析篩選得到FvsR和LvsR組間差異明顯,F(xiàn)vsR的下調(diào)物占比為84.6%,LvsR的下調(diào)物占比為62.5%,表明果實、葉片的生物堿類成分比根中更多,同時可以作為鑒別部位的一類標志物。下調(diào)表示該化合物在果實或葉片中的相對含量比根高,因此重點關注下調(diào)化合物,繼而可以發(fā)現(xiàn)不同部位開發(fā)利用價值的多樣性。

        對比分析組間差異化合物中下調(diào)明顯的前三個物質(zhì)的相對含量,結果如表3和圖4所示,相對于根而言,3-吲哚丙酸、甲氧基吲哚乙酸和油酰單乙醇胺在果實中相對含量更高,油酰單乙醇胺下調(diào)最明顯,因為其在根中未檢測出;3-羥基吡啶、N-甲基煙酰胺和5-羥基吲哚-3-乙醇在葉片中相對含量更高,3-羥基吡啶和N-甲基煙酰胺在根中未檢測出;同時發(fā)現(xiàn)果實、葉片兩者和根相比,兩組的上調(diào)差異化合物中最明顯的都是乙酰色胺,因為其只在根中檢測出。

        有關報道表明,三七果實和葉片中的3-吲哚丙酸、油酰單乙醇胺、3-羥基吡啶和N-甲基煙酰胺物質(zhì)具有進一步研究的價值,可以應用于保健品、藥品等產(chǎn)品開發(fā)。3-吲哚丙酸是一種強大的抗氧化劑,對治療阿爾茨海默病有一定作用[18],也是一種腸道微生物群產(chǎn)生的代謝產(chǎn)物,是2型糖尿?。═2D)發(fā)展的潛在生物標志物,可通過保持β細胞功能來介導其保護作用[19]。油酰單乙醇胺可通過過氧化物酶體增殖物激活受體-α(PPAR-α)介導的機制調(diào)節(jié)脂質(zhì)代謝,可以在肥胖癥治療中發(fā)揮作用[20]。3-羥基吡啶的衍生物是維生素B6的結構類似物,具有抗氧化作用等廣泛的藥理特性[21]。N-甲基煙酰胺是煙酰胺的化合物之一,可通過調(diào)節(jié)ADMA-DDAH軸改善內(nèi)皮功能障礙并減輕動脈粥樣硬化,也影響血管血栓形成[22]。而對甲氧基吲哚乙酸和5-羥基吲哚-3-乙醇有關藥理作用的研究甚少,有待發(fā)掘空間大。甲氧基吲哚乙酸可以根據(jù)其化學式C11H11NO3,參考其同分異構體的化合物——5-甲氧基吲哚-3-乙酸,是一種褪黑素的代謝產(chǎn)物[23]。5-羥基吲哚-3-乙醇曾在哺乳動物中被發(fā)現(xiàn),是5-羥色胺的化合物,可以作為一類急性酒精攝入的標志,可以借為參考[24]。根的特有物質(zhì)乙酰色胺也具備利用價值,色胺是一種單胺生物堿,在大腦中的色胺作為神經(jīng)調(diào)節(jié)劑或神經(jīng)遞質(zhì)發(fā)揮作用,在植物中的色胺是影響植物生長和微生物群落代謝途徑的原料,它曾被發(fā)現(xiàn)可能是植物生長調(diào)節(jié)劑吲哚-3-乙酸生物合成途徑中的一種中間產(chǎn)物[25]。

        3 ?討論與結論

        本研究在林下三七的果實、葉片和根中共檢測出53個生物堿類化合物,從生物堿總含量多少看,

        果實>葉片>根;從成分差異性看,果實和根對比有26個差異化合物,果實中3-吲哚丙酸、甲氧基吲哚乙酸和油酰單乙醇胺的含量較高,葉片和根對比有24個差異化合物,葉片中3-羥基吡啶、N-甲基煙酰胺和5-羥基吲哚-3-乙醇的含量較高;篩選得到果實特有物油酰單乙醇胺,葉片特有物3-羥基吡啶,根特有物乙酰色胺,可以作為鑒別這三個部位的生物堿類差異標志物。

        與常用的高效液相色譜、超高效液相色譜、氣相色譜-質(zhì)譜聯(lián)用等成分分析方法相比,本研究采用的超高效液相色譜和串聯(lián)質(zhì)譜法是一種更為快速、準確、可靠、直觀的方法,為中藥材成分檢測、質(zhì)量評價、部位鑒別等提供新思路[26]。三七的傳統(tǒng)藥用部位主要是根,本研究發(fā)現(xiàn)三七果實和葉片的生物堿總量高于根,與崔秀明等分析總結發(fā)現(xiàn)三七不同藥用部位的藥用成分含量具有差異性的結論一致[27],這說明三七不同部位的功效也可能不相同,揭示三七的果實和葉片這類地上部分也可能有利用價值。曲正義研究三七果實成分,發(fā)現(xiàn)三七果實成分含有抗腫瘤活性[28],孫璠和黃積武等分析三七葉片成分,發(fā)現(xiàn)三七葉片成分含有抗炎活性[29]和抗氧化活性[30],但發(fā)揮功效的具體物質(zhì)及其作用機理尚未明確。本研究發(fā)現(xiàn)果實中含量較高的3-吲哚丙酸和油酰單乙醇胺有抗氧化、調(diào)節(jié)脂質(zhì)代謝的作用,葉片中含量較高的3-羥基吡啶和N-甲基煙酰胺也有抗氧化、保護心血管的作用,結合前人的研究成果說明三七果實和葉片有巨大的開發(fā)利用潛力。

        本研究從林下三七不同部位的生物堿類成分分析入手,旨在發(fā)現(xiàn)提高三七資源綜合利用率的更多可能性,結果為三七不同部位鑒別提供了新的差異標志物,為三七果實和葉片的科學開發(fā)利用提供參考價值。針對性了解三七的生物堿類成分,能夠為開發(fā)創(chuàng)新產(chǎn)品提供科學依據(jù),促進三七資源開發(fā)利用多樣性發(fā)展。

        參考文獻:

        [1] 中華人民共和國藥典2015年版.一部[S].2015:50.

        [2] 劉立紅,劉英,王芬,等.云南三七產(chǎn)業(yè)發(fā)展現(xiàn)狀及發(fā)展建議[J].中國現(xiàn)代中藥,2017,19(9):

        1331-1335.

        [3] 姜藝佼,張思文,姜慶丹,等.“健康中國”視域下2006—2020年中醫(yī)藥健康產(chǎn)業(yè)發(fā)展的研究熱點與趨勢分析[J].中國醫(yī)藥導報,2021,18(31):28-33.

        [4] 詹鵬,韋美膛,徐波,等.基于產(chǎn)業(yè)鏈的我國三七產(chǎn)業(yè)發(fā)展現(xiàn)狀與對策[J].北方園藝,2021(14):158-164.

        [5] 李耿,陳潔,李振坤,等.三七專利數(shù)據(jù)分析與對產(chǎn)業(yè)發(fā)展啟示[J].中國中藥雜志,2020,45(5):

        1011-1022.

        [6] 周家明,張文斌,趙愛,等.文山三七產(chǎn)品加工現(xiàn)狀和發(fā)展趨勢分析[J].中國藥業(yè),2010,19(24):72-73.

        [7] 夏鵬國,張順倉,梁宗鎖,等.三七化學成分的研究歷程和概況[J].中草藥,2014,45(17):2564-2570.

        [8] 唐卿雁,陳庚,盧迎春,等.HPLC同時測定三七不同部位4種皂苷含量[J].熱帶農(nóng)業(yè)科學,2018,38(9):91-95.

        [9] 黃玉香,譚何新,于劍,等.藥用植物生物堿次生代謝工程研究進展[J].中草藥,2016,47(23):

        4271-4281.

        [10] 王騰飛,李信瑜,任佳俐,等.擬康定烏頭中二萜生物堿成分的研究[J].華西藥學雜志,2021,36(2):122-128.

        [11] 馮敏,張文智,張雯霞,等.植物代謝組學技術在中藥資源研究中的應用[J].藥物評價研究,2020,43(4):785-789.

        [12] DELFIN J C, WATANABE M, TOHGE T. Understanding the function and regulation of plant secondary metabolism through metabolomics approaches[J]. Theoretical and Experimental Plant Physiology, 2019, 31: 127-138.

        [13] LIN H Q, ZHU H L, TAN J, et al. Non-targeted metabolomic analysis of methanolic extracts of wild-simulated and field-grown American ginseng[J]. Molecules, 2019, 24(6): 1053.

        [14] WANG C, ZHANG N, WANG Z, et al. Nontargeted metabolomic analysis of four different parts of Platycodon grandiflorum grown in northeast China[J]. Molecules, 2017, 22(8): 1280.

        [15] DAN M, SU M, GAO X, et al. Metabolite profiling of Panax notoginseng using UPLC-ESI-MS[J]. Phytochemistry, 2008, 69(11): 2237-2244.

        [16] ZHU H, LIN H, TAN J, et al. UPLC-QTOF/MS-bsed anontargeted metabolomic analysis of mountain-and garden-cultivated Ginseng of different ages in Northeast China[J]. Molecules, 2019, 24(1): 33.

        [17] 何靜娟,范燕萍.觀賞植物花色相關的類胡蘿卜素組成及代謝調(diào)控研究進展[J].園藝學報,2022,49(5):1162-1172.

        [18] WIKOFF W R, ANFORA A T, LIU J, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites[J]. Proceedings of the national academy of sciences, 2009, 106(10): 3698-3703.

        [19] DE M V D, PAANANEN J, LINDSTR?M J, et al. Indolepropionic acid and novel lipid metabolites are associated with a lower risk of type 2 diabetes in the Finnish Diabetes Prevention Study[J]. Scientific reports, 2017, 7(1): 46337.

        [20] CHEN L, LI L, CHEN J, et al. Oleoylethanolamide, an endogenous PPAR-α ligand, attenuates liver fibrosis targeting hepatic stellate cells[J]. Oncotarget, 2015, 6(40): 42530.

        [21] LUKYANOVA L D, GERMANOVA E L, TSYBINA T A, et al. Energotropic effect of succinate-containing derivatives of 3-hydroxypyridine[J]. Bulletin of experimental biology and medicine, 2009, 148: 587-591.

        [22] MOGIELNICKI A, KRAMKOWSKI K, PIETRZAK L, et al. N-methylnicotinamide inhibits arterial[J]. Journal of Physiology and Pharmacology, 2007, 58(3): 515-527.

        [23] KVEDER S, MCISAAC W M. The metabolism of melatonin (N-acetyl-5-methoxytryptamine) and

        5-methoxytryptamine[J]. Journal of Biological

        Chemistry, 1961, 236(12): 3214-3220.

        [24] BECK O, HELANDER A. 5-Hydroxytryptophol as a marker for recent alcohol intake[J]. Addiction, 2003, 98: 63-72.

        [25] JONES R S G. Tryptamine: a neuromodulator or neurotransmitter in mammalian brain?[J]. Progress in neurobiology, 1982, 19(1-2): 117-139.

        [26] 鐘貴,王元忠,張霽,等.人參屬植物代謝組學研究進展[J].時珍國醫(yī)國藥,2013,24(9):2234-2237.

        [27] 崔秀明,徐洛珊,王強,等.三七中三七素的含量測定[J].中國藥學雜志,2005,40(13):1017-101.

        [28] 曲正義.三七果化學成分及抗腫瘤活性研究[D].北京:中國農(nóng)業(yè)科學院,2021.

        [29] 孫璠.三七葉化學成分及其體外抗炎活性的研究[D].天津:天津中醫(yī)藥大學,2020.

        [30] 黃積武,李創(chuàng)軍,楊敬芝,等.三七葉的化學成分及其抗氧化活性研究[J].中草藥,2017,48(21):

        4381-4386.

        (責任編輯:丁志祥)

        猜你喜歡
        代謝組學
        何首烏致肝損傷大鼠的動態(tài)血清代謝組學研究
        基于UPLC—Q—TOF—MS技術的牛血清白蛋白誘導過敏反應的代謝組學研究
        基于UPLC—Q—TOF—MS技術的牛血清白蛋白誘導過敏反應的代謝組學研究
        藥用植物代謝組學研究
        藏藥鑒定及質(zhì)量控制研究現(xiàn)狀
        新疆維吾爾族和漢族大學生尿液代謝輪廓差異分析
        分析化學(2015年10期)2015-11-03 07:42:21
        枯草芽孢桿菌代謝組樣品前處理方法的比較研究
        分析化學(2015年8期)2015-08-13 07:23:40
        鎘超富集植物東南景天根系分泌物的代謝組學研究
        分析化學(2015年1期)2015-01-20 06:52:54
        国产精品久久久久久福利| 男人的天堂av一二三区| 亚洲高清av一区二区| 五月婷婷开心六月激情| 国产人妖乱国产精品人妖| 精品国产一二三产品区别在哪| 亚洲欧美日本| 五月天无码| 久久精品国产亚洲av专区| 日韩精品在线视频一二三| 精品一区二区三区免费视频| 理论片午午伦夜理片影院| 白丝美女被狂躁免费视频网站| 亚洲一区二区三区一区| 国产av一级黄一区二区三区| 吃奶呻吟打开双腿做受视频| 又黄又爽又色又刺激的视频| 国产精品乱子伦一区二区三区| 中文乱码字幕人妻熟女人妻| 久久九九精品国产av| 人人妻人人爽人人澡欧美一区 | 牛鞭伸入女人下身的真视频| 91青青草久久| 香蕉蜜桃av一区二区三区| 国产精品国产高清国产专区| 一本色道久久88精品综合| 国产午夜无码视频免费网站| 亚洲av熟女天堂系列| 亚洲激情一区二区三区不卡| 国产亚洲精品a片久久久| 无码a∨高潮抽搐流白浆| 91极品尤物国产在线播放| 快射视频网站在线观看| 亚洲av色香蕉一区二区三区| 少妇做爰免费视频网站| 亚洲九九九| 亚洲精彩av大片在线观看| 射精区-区区三区| 五月婷婷俺也去开心| 国产精品欧美成人片| 日韩精品极品系列在线免费视频|