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

    基于GC-MS研究不同產(chǎn)地三色藜麥種子的代謝物差異

    2019-06-11 08:27時(shí)羽杰李興龍唐媛董雨鄔曉勇
    關(guān)鍵詞:抗逆性

    時(shí)羽杰 李興龍 唐媛 董雨 鄔曉勇

    摘要:為了研究不同生境對(duì)藜麥發(fā)育的影響,以不同產(chǎn)地的藜麥種子為研究對(duì)象,篩選顯著差異代謝物質(zhì)和相關(guān)代謝通路,利用GC-MS平臺(tái)進(jìn)行非靶向代謝組學(xué)分析,分析了青海都蘭縣和四川鹽源縣兩地產(chǎn)出的黑藜、紅藜、白藜種子的代謝輪廓,將原始數(shù)據(jù)經(jīng)XCMS online平臺(tái)預(yù)處理后進(jìn)行后續(xù)多元統(tǒng)計(jì)分析。結(jié)果表明,青海都蘭縣和四川鹽源縣兩地的黑、紅、白藜麥種子的代謝特征均不相同,其中紅藜麥種子與白藜麥種子差異較大,而黑藜麥種子差異最小;紅藜存在顯著差異,共有特征峰570個(gè),上調(diào)物少于下調(diào)物;黑藜差異最小,上、下調(diào)趨勢(shì)與紅藜一致;白藜差異較大,上、下調(diào)趨勢(shì)與紅藜和黑藜相反,共有特征峰486個(gè),上調(diào)物質(zhì)多于下調(diào)物質(zhì);不同生境對(duì)藜麥的影響較大,黑藜和紅藜更適宜在青海都蘭縣種植,而白藜更適宜在四川鹽源縣種植。研究結(jié)果可以指導(dǎo)青海都蘭縣和四川鹽源縣兩地藜麥的栽培育種,也為深入研究藜麥的抗逆性機(jī)理提供了理論依據(jù)。

    關(guān)鍵詞:作物栽培學(xué);藜麥;GC-MS;XCMS;代謝差異;抗逆性

    中圖分類號(hào):TS210.1文獻(xiàn)標(biāo)志碼:A

    Abstract:In order to study the effects of different habitats on the development of quinoa, taking quinoa seeds from different habitats as research objects, significant differences in metabolites and related metabolic pathways are screened. The metabolic profiles of the seeds of black quinoa, red quinoa and white quinoa from Dulan County, Qinghai Province and Yanyuan County, Sichuan Province are analyzed through non-targeted metabonomics analysis with GC-MS platform. The original data are pretreated by XCMS online platform for subsequent multivariate statistical analysis. The results show that the metabolic characteristics of black, red and white seeds in Dulan county and Yanyuan County are different. Among them, the difference between red and white seeds from the two counties is larger, while the difference of the black quinoa seeds is the smallest. There are 570 characteristic peaks in red seeds which have significant differences, and the up-regulated substances are less than the down-regulated. The difference between black quinoa from two habitats is the smallest, and the up- and down-regulated trends are like that of red quinoa. There is secondly significant difference between the white seeds from the two habitats, and the up- and down-regulated trends are the opposite of that of the red and black quinoas, and 486 characteristic peaks are detected, with more up-regulated substances than down-regulated. Different habitats have great impact on quinoa. Black and red quinoas are more suitable to be planted in Dulan County, Qinghai Province, while white quinoa is more suitable to be planted in Yanyuan County, Sichuan Province. The conclusions could guide the cultivation and breeding of quinoa in Dulan county and Yanyuan County, and provide theoretical basis for further study of the stress resistance mechanism of quinoa.

    Keywords:crop cultivation;quinoa;GC-MS;XCMS;metabolic differences;stress resistance

    代謝組學(xué)(metabonomics/metabolomics)是繼基因組學(xué)、蛋白組學(xué)之后又一門(mén)新興的學(xué)科,是系統(tǒng)生物學(xué)的重要組成部分,旨在研究生物體或組織,甚至單個(gè)細(xì)胞的全部小分子代謝物(MV<1 000)成分及其動(dòng)態(tài)變化[1],根據(jù)分析目的的不同又可以分為靶向代謝組學(xué)和非靶向代謝組學(xué)。靶向代謝組學(xué)主要針對(duì)幾種目標(biāo)化合物或某條通路上涉及的部分代謝物,利用標(biāo)準(zhǔn)品,構(gòu)建特異性強(qiáng)、靈敏度高、重復(fù)性好的檢測(cè)方法,對(duì)目標(biāo)物質(zhì)進(jìn)行分析。非靶向代謝輪廓分析方法面向所有能檢測(cè)到的化合物,可以較全面地反映整體代謝狀態(tài)的變化,并通過(guò)多元統(tǒng)計(jì)分析篩選差異代謝物,對(duì)差異代謝物進(jìn)行通路分析,揭示其變化的生理機(jī)制,在轉(zhuǎn)基因植物非預(yù)期效應(yīng)研究[2]、非生物與生物脅迫研究中顯示出了極大的潛力。

    藜麥,古代印加人稱之為“糧食之母”[3]。藜麥含有蛋白質(zhì)、礦物質(zhì)、氨基酸、纖維素、維生素等營(yíng)養(yǎng)成分,含量明顯高于大多數(shù)常見(jiàn)谷物[4]。聯(lián)合國(guó)糧農(nóng)組織于2013年正式將藜麥推薦為最適宜人類的“全營(yíng)養(yǎng)食品”[5]。藜麥還含有豐富的多酚類、黃酮類、皂苷等活性成分[6],具有耐寒、耐旱、耐瘠薄、耐鹽堿等生理特點(diǎn)。藜麥種子主要有白、黑、紅、黃等幾種顏色系,營(yíng)養(yǎng)品質(zhì)差異較大,藜麥中含有生物堿、多酚等物質(zhì)。黑藜中總多酚的含量最高,明顯高于紅藜和白藜,而紅藜和白藜中總多酚含量差異不大[7]。趙亞?wèn)|[8]研究發(fā)現(xiàn)青海的不同色系藜麥中,紅藜含皂苷最多,抗氧化性極強(qiáng),白藜抗氧化性最弱。這些研究均表明,藜麥的生境對(duì)籽粒發(fā)育影響顯著。鑒于目前國(guó)內(nèi)從代謝組學(xué)方面對(duì)不同產(chǎn)地藜麥的研究甚少,本研究主要以非靶向代謝組學(xué)研究青海都蘭縣和四川鹽源縣兩地的黑、白、紅藜麥,以獲取較為全面的代謝物質(zhì)信息,全面、準(zhǔn)確地分析其差異物質(zhì)和代謝通路,分析3種藜麥種子代謝物在不同生境下的異同,為揭示藜麥在不同栽培地區(qū)造成的差異等后續(xù)研究提供理論依據(jù),指導(dǎo)藜麥的耕作與栽培。

    1主要材料與儀器

    1.1材料

    分別選2018年采收自青海省海西蒙古族藏族自治州都蘭縣(平均海拔3 000 m)和四川省涼山州的鹽源縣(平均海拔2 000 m)的黑、紅、白色藜麥種子,經(jīng)發(fā)芽試驗(yàn)篩選,取發(fā)芽率和發(fā)芽勢(shì)相同的一批進(jìn)行實(shí)驗(yàn),每一個(gè)樣品均重復(fù)6次。文中所提到的青海都蘭縣黑藜、紅藜、白藜均以A,C,E表示,四川鹽源縣黑藜、紅藜、白藜均以B,D,F(xiàn)表示。

    1.2試劑

    核糖醇,甲氧氨基鹽酸鹽,N-甲基-N-三甲基硅基三氟乙酰胺(MSTFA),購(gòu)自Sigma Aldrich公司;吡啶,甲醇,氯仿,三甲基氯硅烷均為色譜純,成都科隆化學(xué)品有限公司提供。

    1.3儀器

    Agilent GC-MS 7890B 5977氣質(zhì)聯(lián)用儀,產(chǎn)地美國(guó);DB-5MS毛細(xì)管氣相色譜柱(30 m×0.25 mm,0.25 μm),購(gòu)自安捷倫公司(美國(guó));THZ-82水浴恒溫振蕩器,常州榮華儀器制造有限公司提供;通風(fēng)櫥、H1850R臺(tái)式高速冷凍離心機(jī),產(chǎn)地湖南;ZLS-2真空離心濃縮儀,湖南赫西儀器有限公司提供。

    2實(shí)驗(yàn)方法

    2.1樣品制備

    分別取來(lái)自青海都蘭縣和四川鹽源縣兩地的黑、紅、白藜麥種子,參考LISEC等[9]的研究方法對(duì)植物樣本進(jìn)行提取,條件稍作優(yōu)化。對(duì)樣品進(jìn)行6次生物學(xué)重復(fù),并做質(zhì)控樣本QC,QC樣本是將所有提取物的均等混合樣隨機(jī)插入檢測(cè)的樣本序列中,主要用于實(shí)驗(yàn)方法和儀器穩(wěn)定的校正。將所有樣品置于EP管,于真空離心濃縮儀中進(jìn)行揮干備用。

    2.2樣品衍生化

    在上述干燥后的樣品管中加入60 μL甲氧氨基鹽酸鹽(20 mg/mL吡啶溶液),于37 ℃振蕩反應(yīng)2 h。反應(yīng)完畢后再加入100 μL的N-甲基-N-三甲基硅基三氟乙酰胺和1%(體積分?jǐn)?shù),下同)三甲基氯硅烷(MSTFA+1%TMCS),于37 ℃振蕩反應(yīng)30 min。將衍生化完成后的樣品用0.22 μm濾膜過(guò)濾[10]后轉(zhuǎn)移至適合GC-MS分析的內(nèi)襯管中。

    2.3GC-MS分析條件

    2.3.1色譜條件

    采用DB-5MS(30 m×0.25 mm,0.25 μm)毛細(xì)管氣相色譜柱;升溫程序:起始溫度保持在80 ℃,2 min后以10 ℃/min升至190 ℃,再以3 ℃/min升至220 ℃,最后以20 ℃/min升至280 ℃,保持10 min;載氣:氦氣(體積分?jǐn)?shù)為99.999%);恒流,流速為1 mL/min;進(jìn)樣模式:不分流進(jìn)樣,進(jìn)樣量為1 μL。

    2.3.2質(zhì)譜條件

    電子碰撞內(nèi)能:70 eV;離子源溫度:220 ℃;傳輸線溫度:270 ℃;注射劑溫度:280 ℃;掃描速度:每20 s光譜全波段掃描法捕獲;溶劑延遲:420 s;質(zhì)量掃描范圍m/z:70~500。

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

    將GC-MS原始數(shù)據(jù)[11]上傳至XCMS online程序,進(jìn)行行峰識(shí)別(peaks identification)、峰過(guò)濾(peaks filtration)、峰對(duì)齊(peaks alignment)等[12],同時(shí)設(shè)置相應(yīng)參數(shù)[13],最后得到包括質(zhì)核比(m/z)和保留時(shí)間(rt)及峰面積等信息的數(shù)據(jù)矩陣。將數(shù)據(jù)導(dǎo)出至Excel進(jìn)行后續(xù)整理,用SIMCA-P(v14.0) 軟件進(jìn)行主成分分析(PCA)和正交-偏最小二乘判別分析(OPLS-DA),結(jié)合XCMS分析篩選差異代謝物,進(jìn)行后續(xù)的差異代謝物統(tǒng)計(jì)分析。

    3結(jié)果與分析

    3.1主成分分析

    原始數(shù)據(jù)經(jīng)XCMS進(jìn)行峰識(shí)別、過(guò)濾、對(duì)齊等處理后,結(jié)果見(jiàn)圖1。將處理后的數(shù)據(jù)整理后導(dǎo)入SIMCA-P軟件中進(jìn)行主成分分析(PCA),用UV處理數(shù)據(jù)進(jìn)行擬合,擬合后方程的R2X=0.855>0.5,說(shuō)明擬合性好,且Q2=0.648,與R2X相差較小,說(shuō)明擬合方程的穩(wěn)定性好。圖2為得到的PCA得分,由圖2可知,QC的重復(fù)樣本間較為聚集,說(shuō)明實(shí)驗(yàn)方法和檢測(cè)儀器較為穩(wěn)定,準(zhǔn)確性較好。此外,6組樣品中除了B6,D3和E3外的所有點(diǎn)基本都處在95%的置信區(qū)間內(nèi),且每一組內(nèi)的重復(fù)樣品較為聚集,而各組間的樣品又相對(duì)較為分散,組與組之間有明顯的界限,說(shuō)明樣品組內(nèi)的重復(fù)性好,數(shù)據(jù)較為準(zhǔn)確;而樣品組間的明顯界限說(shuō)明,各組間存在著明顯的差異。還可以看出,在A組、C組、E組間,A組與C組、E組在第一主成分上被明顯分開(kāi),其中A組被分在了第一主成分的左邊,而C組和E組被分在了第一主成分的右邊,說(shuō)明A組與C組、E組間差異較大,C組和E組間差異較小。B組、D組、F組間也存在著一樣的差異趨勢(shì),說(shuō)明同一產(chǎn)地的不同色系藜麥間也存在著明顯差異。用DModX和置信區(qū)間進(jìn)行初步篩選,剔除B6,D3,E2這3個(gè)異常點(diǎn)[14],以保證后續(xù)數(shù)據(jù)分析的準(zhǔn)確性。

    PCA只能反映樣品的整體趨勢(shì),不利于發(fā)現(xiàn)組間差異和差異化合物,所以需要通過(guò)正交-偏最小二乘法判別(OPLS-DA)來(lái)更為準(zhǔn)確地分析數(shù)據(jù),尋找組間的差異[15]。首先根據(jù)種子顏色進(jìn)行分組:A-B組、C-D組和E-F組,用Par算法進(jìn)行數(shù)據(jù)轉(zhuǎn)換,用OPLS-DA方程進(jìn)行擬合,得到擬合方程的模型參數(shù)見(jiàn)表1。其中所有的R2X,R2Y和Q2都大于0.5且各參數(shù)間差異較小,說(shuō)明擬合方程較為可靠。由方程擬合得到OPLS-DA得分圖,見(jiàn)圖3中的a),b)和c)。由圖3可知,A-B組、C-D組和E-F組的分離程度較大,說(shuō)明各組間有明顯差異,且組內(nèi)點(diǎn)相對(duì)較為聚集,表明樣品間重復(fù)性好且較為穩(wěn)定。再經(jīng)過(guò)OPLS-DA置換檢驗(yàn)來(lái)判斷方程是否過(guò)擬合,得到檢驗(yàn)圖(見(jiàn)圖3中的d),e)和f))。圖3中所有的Q2點(diǎn)從左到右均低于最右邊的原始Q2點(diǎn),且Q2<0,R2和Q2的回歸線與橫坐標(biāo)交叉或者小于0,說(shuō)明評(píng)估模型可靠有效[16]。

    3.3篩選顯著性差異代謝物

    采用XCMS的Pairwise中的對(duì)比分析,上傳A-B組、C-D組和E-F組數(shù)據(jù)進(jìn)行分組分析處理,得到Metabolomic Cloud Plot圖,設(shè)置參數(shù)p-value<0.01,fold change>1.5,進(jìn)行顯著性差異代謝物的篩選,得到代謝物的云圖,見(jiàn)圖4。圖中M線上部為上調(diào),M線下部為下調(diào),橫坐標(biāo)為保留時(shí)間,縱坐標(biāo)是m/z。圖4顯示,A-B組有140個(gè)差異代謝物,上調(diào)的有90個(gè),下調(diào)的有50個(gè);C-D組有570個(gè)差異代謝物,上調(diào)的有209個(gè),下調(diào)的有361個(gè);E-F組有486個(gè)代謝物,上調(diào)的有379個(gè),下調(diào)的有107個(gè)。其中A-B組中上調(diào)物質(zhì)占總物質(zhì)的64.28%,下調(diào)物質(zhì)占35.72%;E-F組中上調(diào)物質(zhì)占總物質(zhì)的77.98%,下調(diào)物質(zhì)占總物質(zhì)的22.02%,兩組趨勢(shì)和比例相差不大。C-D組中,上調(diào)物質(zhì)占總物質(zhì)的36.67%,下調(diào)物質(zhì)占總物質(zhì)的6333%,與A-B組和E-F組完全相反。其中C-D組間差異物質(zhì)最多,差異最大,這一結(jié)果也與胡一晨等[7]的報(bào)道一致。

    3.4差異代謝物通路分析

    將上述篩選鑒定的差異代謝物上傳至Metaboanalyst,利用MetPA數(shù)據(jù)庫(kù)對(duì)其進(jìn)行超幾何檢驗(yàn)和拓?fù)浣Y(jié)構(gòu)分析,根據(jù)代謝通路的濃縮,識(shí)別出可能的受生物擾動(dòng)的代謝通路[17]。對(duì)差異代謝物和通路進(jìn)行相關(guān)性分析,得到多條相關(guān)通路。根據(jù)-log p值和Impact值綜合分析,每一組篩選出最為顯著的2條代謝通路,如表2所示。

    由表2可知,黑藜中主要是吲哚生物堿合成途徑,合成的生物堿可以降血糖和血脂,有效治療糖尿病,與房娜[18]的研究結(jié)果一致;泛酸與CoA生物合成途徑與植物的碳水化合物、脂肪酸、蛋白質(zhì)和能量代謝密切相關(guān)。紅藜中主要是硫代謝途徑,它與植物抗逆性密切相關(guān),其途徑中SO2-4經(jīng)過(guò)植物一系列還原與同化反應(yīng)后進(jìn)入有機(jī)骨架,生成半胱氨酸(Cys),植物再以Cys為前體,合成眾多具有重要生物學(xué)功能的代謝產(chǎn)物,直接關(guān)系到植物耐逆境脅迫和農(nóng)作物的產(chǎn)量與品質(zhì),苯丙氨酸代謝途徑與植物的防御系統(tǒng)密切相關(guān)。白藜中主要是谷胱甘肽代謝途徑,它與植物直接和間接的抗逆功能活動(dòng)密切相關(guān)[19];丙氨酸、天冬氨酸和谷氨酸代謝途徑與糖代謝密切相關(guān)。

    4討論與結(jié)論

    青海省海西蒙古族藏族自治州都蘭縣是典型的高原大陸性氣候,日照時(shí)間長(zhǎng),太陽(yáng)輻射強(qiáng),晝夜溫差大,常年干旱、多風(fēng)、少雨,年平均氣溫4 ℃左右;而四川省涼山州鹽源縣的日溫差較大,年均氣溫12.1 ℃,最高溫度30.7 ℃,最低溫度-11.3 ℃;雖然兩地都較適宜藜麥的種植,但不同藜麥品種間的適應(yīng)性不同。

    本研究對(duì)青海都蘭縣和四川鹽源縣兩地的黑、紅、白藜麥種子做了代謝組學(xué)分析,發(fā)現(xiàn)兩地種子之間都存在著明顯差異,其中紅藜麥種子與白藜麥種子差異較大,而黑藜麥種子差異最小。差異物質(zhì)主要是脂肪酸(棕櫚酸、肉桂酸、油酸)、糖類(乳糖、葡萄糖、脫氧糖苷)和氨基酸(脯氨酸、絲氨酸、組氨酸)等。通過(guò)比較發(fā)現(xiàn),四川鹽源縣與青海都蘭縣的黑藜相比,差異代謝物纈氨酸表達(dá)下調(diào),而纈氨酸是合成α-酮異戊酸的前體,α-酮異戊酸在酶作用下生成酮泛解酸,再與β-Ala合成泛酸。纈氨酸的下調(diào)導(dǎo)致泛酸合成受阻,使得植物體的碳水化合物、脂肪酸、蛋白質(zhì)和能量代謝受阻。四川鹽源縣和青海都蘭縣的紅藜相比,差異代謝物O-Acetylserine表現(xiàn)為顯著下調(diào),而這與半胱氨酸(Cys)密切相關(guān)。O-Acetylserine是Cys合成的前體,Cys合成受阻會(huì)導(dǎo)致硫代謝受阻。這些結(jié)果均表明,青海都蘭縣的黑藜與紅藜均比四川鹽源縣有較強(qiáng)的耐受性,均能通過(guò)調(diào)節(jié)體內(nèi)能量代謝通路的變化更好地去適應(yīng)青海相對(duì)低溫的生境。而四川鹽源縣和青海都蘭縣的白藜相比,差異代謝物谷胱甘肽(還原型)表達(dá)上調(diào)。因?yàn)楣入赘孰挠泻軓?qiáng)的抗氧化性,其上調(diào)增加直接影響到植物的抗逆性,使得四川鹽源縣的白藜能更快地通過(guò)代謝通路的變化去消除體內(nèi)過(guò)氧化物以及自由基等,使藜麥進(jìn)行正常的生長(zhǎng),因此四川鹽源縣的白藜比青海都蘭縣的白藜有更好的適應(yīng)性,而這一研究結(jié)果也與青海都蘭縣白藜麥栽培的結(jié)果一致——溫度的變化會(huì)導(dǎo)致白藜生長(zhǎng)緩慢[21]。

    本實(shí)驗(yàn)通過(guò)代謝組學(xué)分析四川鹽源縣與青海都蘭縣的黑、紅、白藜,發(fā)現(xiàn)兩地的同種藜麥間代謝物質(zhì)有很大差異。在此基礎(chǔ)上通過(guò)分析差異代謝物的代謝通路,發(fā)現(xiàn)青海都蘭縣的黑藜和紅藜比四川鹽源縣具有更好的抗逆性,而青海都蘭縣的白藜較四川鹽源縣的白藜抗逆性弱??寡趸耘c抗逆性密切相關(guān),這與趙亞?wèn)|等[22]的研究結(jié)果一致。所以紅藜和黑藜更適宜在青海都蘭縣種植,而白藜更適宜在四川鹽源縣種植。研究結(jié)果可為后續(xù)在代謝途徑層面上深入研究藜麥的抗逆性機(jī)理提供理論依據(jù)。

    本研究所選取的樣本種類過(guò)少,在對(duì)不同生境下藜麥種子進(jìn)行分析時(shí)只進(jìn)行了代謝層面的分析,缺少對(duì)通路中關(guān)鍵蛋白質(zhì)和關(guān)鍵酶的分析,對(duì)整個(gè)機(jī)理的解釋存在著一定局限性。后續(xù)將結(jié)合多組學(xué)技術(shù)包括基因組學(xué)、轉(zhuǎn)錄組學(xué)和蛋白質(zhì)組學(xué)等,對(duì)不同生境下黑、紅、白藜麥的差異進(jìn)行更加全面深入的分析,以期選育出優(yōu)質(zhì)、高產(chǎn)、抗逆性強(qiáng)的藜麥新品種,為藜麥的栽培種植提供理論依據(jù)。

    參考文獻(xiàn)/References:

    [1]FIEHN O. Metabolomics the Link Between Genotypes and Phenotypes[M]. Berlin:Springer Netherlands, 2002:155-171.

    [2]DUBOUZET J G, ISHIHARA A, MATSUDA F, et al. Integrated metabolomic and transcriptomic analyses of high-tryptophan rice expressing a mutant anthranilate synthase alpha subunit[J]. Journal of Experimental Botany, 2007, 58(12):3309-3321.

    [3]FERREIRA D S, PALLONE J A L, POPPI R J. Direct analysis of the main chemical constituents in Chenopodium quinoa, grain using Fourier transform near-infrared spectroscopy[J]. Food Control, 2015, 48(48):91-95.

    [4]VEGAGLVEZ A, MIRANDA M, VERGARA J, et al. Nutrition facts and functional potential of quinoa (Chenopodium quinoa Willd.), an ancient Andean grain: A review[J]. Journal of the Science of Food & Agriculture, 2010, 90(15):2541-2547.

    [5]JAMES L E A. Chapter 1 quinoa (Chenopodium quinoa, Willd.):Composition, chemistry, nutritional, and functional properties[J]. Advances in Food & Nutrition Research, 2009, 58(9):1-10.

    [6]肖正春, 張廣倫. 藜麥及其資源開(kāi)發(fā)利用[J]. 中國(guó)野生植物資源, 2014, 33(2):62-66.XIAO Zhengchun, ZHANG Guanglun. Development and utilization of Chenopodium quinoa Willd.[J].Chinese Wild Plant Resources, 2014, 33(2):62-66.

    [7]胡一晨,趙鋼,秦培友,等. 藜麥活性成分研究進(jìn)展[J].作物學(xué)報(bào),2018,44(11):1579-1591.

    HU Yichen, ZHAO Gang, QIN Peiyou, et al. Research progress on bioactive components of quinoa (Chenopodium quinoa Willd.)[J].Acta Agronomica Sinica,2018,44(11):1579-1591.

    [8]趙亞?wèn)|. 青海藜麥資源營(yíng)養(yǎng)品質(zhì)評(píng)價(jià)及功能成分與抗氧化活性研究[D].西寧:青海大學(xué),2018.

    ZHAO Yadong. Nutritional Quality Evaluation and Study of the Functional Compositions and Antioxidant Activities in Quinoa from Qinghai[D]. Xining:Qinghai University, 2018.

    [9]LISEC J, SCHAUER N, KOPKA J, et al. Gas chromatography mass spectrometry-based metabolite profiling in plants[J]. Nature Protocols, 2006, 1(1):387-396.

    [10]哈婧, 王嬌嬌,趙帆帆,等. 石韋配方顆粒HPLC-DAD特征圖譜研究[J]. 河北科技大學(xué)學(xué)報(bào), 2017, 38(1):39-45.

    HA Jing,WANG Jiaojiao,ZHAO Fanfan,et al.Study of HPLC-DAD characteristic chromatogram of Pyrrosia leaf formula granules[J]. Journal of Hebei University of Science and Technology, 2017, 38(1):39-45.

    [11]COLIN A S, ELIZABETH J W, GRACE O, et al. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification[J]. Analytical Chemistry, 2006, 78(3):779-787.

    [12]FORSBERG E M, TAO H, RINEHART D, et al. Data processing, multiomic pathway mapping, and metabolite activity analysis using XCMS online[J]. Nature Protocols, 2018, 13(4):633-651.

    [13]POGACIC T, MAILLARD M B, LECLERC A, et al. A methodological approach to screen diverse cheese-related bacteria for their ability to produce aroma compounds[J]. Food Microbiology, 2015, 46:145-153.

    [14]阿基業(yè). 代謝組學(xué)數(shù)據(jù)處理方法——主成分分析[J]. 中國(guó)臨床藥理學(xué)與治療學(xué), 2010, 15(5):481-489.A Jiye. Analysis of metabolomic data:Principal component analysis[J]. Chinese Journal of Clinical Pharmacology and Therapeutics, 2010, 15(5):481-489.

    [15]JOHAN A W, EWOUD J J, HUUB C J H, et al. Multivariate paired data analysis: Multilevel PLSDA versus OPLSDA[J]. Metabolomics, 2010, 6(1):119-128.

    [16]BOULESTEIX A L, STRIMMER K. Partial least squares: A versatile tool for the analysis of high-dimensional genomic data[J]. Briefings in Bioinformatics, 2007, 8(1):32-37.

    [17]OGATA H, GOTO S, SATO K, et al. KEGG: Kyoto encyclopedia of genes and genomes[J]. Nucleic Acids Research, 2000, 27(1):29-34.

    [18]房娜. 觀察藜麥飲食對(duì)2型糖尿病患者餐后血糖影響[J]. 中國(guó)藥物與臨床, 2017, 17(8):1180-1182.

    [19]張潔琳, 付暢. 半胱氨酸蛋白酶及半胱氨酸蛋白酶抑制劑對(duì)逆境脅迫的應(yīng)答[J]. 分子植物育種, 2018, 16(13):4453-4459.ZHANG Jielin, FU Chang. Response of cysteine protease and cysteine protease inhibitor to adverse stress[J].Molecular Plant Breeding, 2018, 16(13):4453-4459.

    [20]孫國(guó)祥,張旭,張成國(guó),等.谷胱甘肽在動(dòng)植物領(lǐng)域的應(yīng)用研究進(jìn)展[J].安徽農(nóng)業(yè)科學(xué),2018,46(8):42-45.SUN Guoxiang, ZHANG Xu, ZHANG Chengguo,et al. Application of glutathione in animal and plant fields[J]. Journal of Anhui Agricultural Sciences, 2018,46(8): 42-45.

    [21]劉洋, 閆殿海, 張玉清, 等. 青海藜麥的發(fā)展現(xiàn)狀與發(fā)展中存在的問(wèn)題[J]. 青海農(nóng)林科技, 2015(4):31-34.LIU Yang, YAN Dianhai, ZHANG Yuqing, et al. The current situation and problems in development of Chenopodium quinoa in Qinghai[J].Science and Technology of Qinghai Agriculture and Forestry, 2015(4):31-34.

    [22]趙亞?wèn)|, 黨斌, 楊希娟, 等. 青海藜麥皂苷超聲提取工藝及抗氧化活性[J]. 食品工業(yè)科技, 2017,38(19):45-52.

    猜你喜歡
    抗逆性
    球孢白僵菌植物內(nèi)生性及其應(yīng)用特性研究進(jìn)展
    2018年貴州省玉米區(qū)域試驗(yàn)L組六枝點(diǎn)總結(jié)
    體育教育對(duì)高職院校學(xué)生心理資本干預(yù)研究
    植物生長(zhǎng)調(diào)節(jié)劑誘導(dǎo)植物抗逆性研究進(jìn)展
    亚洲精品日韩在线中文字幕| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区 | √禁漫天堂资源中文www| 在线观看免费高清a一片| 亚洲精品成人av观看孕妇| 美女中出高潮动态图| 成人亚洲精品一区在线观看| 成人国产麻豆网| 青春草亚洲视频在线观看| 久久国内精品自在自线图片| 国产高清三级在线| 99精国产麻豆久久婷婷| 日韩三级伦理在线观看| 午夜福利影视在线免费观看| 在线观看人妻少妇| 中文字幕人妻丝袜制服| 久久这里只有精品19| 最近中文字幕高清免费大全6| 国产精品一区www在线观看| 国产亚洲av片在线观看秒播厂| 亚洲在久久综合| 丰满少妇做爰视频| 国产色爽女视频免费观看| 亚洲av电影在线进入| 中文字幕人妻熟女乱码| 观看美女的网站| 麻豆乱淫一区二区| 欧美激情国产日韩精品一区| 纵有疾风起免费观看全集完整版| 最近最新中文字幕免费大全7| 黄色怎么调成土黄色| 免费黄色在线免费观看| 欧美老熟妇乱子伦牲交| 街头女战士在线观看网站| 91久久精品国产一区二区三区| 婷婷色综合www| 97在线人人人人妻| 亚洲经典国产精华液单| 色婷婷av一区二区三区视频| 一区二区三区乱码不卡18| 九九爱精品视频在线观看| 成年人免费黄色播放视频| 成人手机av| 国产视频首页在线观看| 久久99热6这里只有精品| 亚洲精品一二三| 一级,二级,三级黄色视频| 精品久久国产蜜桃| 国产日韩欧美视频二区| 亚洲精华国产精华液的使用体验| 蜜臀久久99精品久久宅男| 97在线视频观看| 国产亚洲精品久久久com| 91成人精品电影| 大陆偷拍与自拍| 国产一区二区在线观看av| 99久久综合免费| 婷婷成人精品国产| 久久久国产精品麻豆| 99热国产这里只有精品6| 一区二区三区四区激情视频| 中文字幕制服av| 伦理电影大哥的女人| 国产一区二区在线观看日韩| 一级毛片我不卡| 亚洲伊人色综图| 色94色欧美一区二区| 又大又黄又爽视频免费| 美女内射精品一级片tv| 天天影视国产精品| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| 国产乱来视频区| 日韩在线高清观看一区二区三区| 国产精品偷伦视频观看了| 成人18禁高潮啪啪吃奶动态图| 午夜激情av网站| 另类亚洲欧美激情| 蜜臀久久99精品久久宅男| 国产高清国产精品国产三级| 晚上一个人看的免费电影| 免费大片黄手机在线观看| 69精品国产乱码久久久| 99九九在线精品视频| 亚洲av在线观看美女高潮| 免费av不卡在线播放| 七月丁香在线播放| 国产成人aa在线观看| 国产av一区二区精品久久| 亚洲av国产av综合av卡| 一二三四中文在线观看免费高清| 成年av动漫网址| 看十八女毛片水多多多| 草草在线视频免费看| 毛片一级片免费看久久久久| 十八禁高潮呻吟视频| 九九在线视频观看精品| 国产精品欧美亚洲77777| 午夜激情av网站| 国产精品不卡视频一区二区| 99热国产这里只有精品6| 久久精品国产综合久久久 | 国产精品三级大全| 蜜桃国产av成人99| 99国产综合亚洲精品| 久久久久视频综合| 少妇熟女欧美另类| 人人妻人人澡人人爽人人夜夜| 99热国产这里只有精品6| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久成人aⅴ小说| 这个男人来自地球电影免费观看 | 大码成人一级视频| 麻豆精品久久久久久蜜桃| 亚洲性久久影院| 免费在线观看黄色视频的| 18禁在线无遮挡免费观看视频| 高清毛片免费看| 亚洲人成77777在线视频| 一本色道久久久久久精品综合| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久| 九九在线视频观看精品| 青春草亚洲视频在线观看| 欧美日韩亚洲高清精品| 中文字幕另类日韩欧美亚洲嫩草| 精品久久国产蜜桃| 天天躁夜夜躁狠狠躁躁| 精品熟女少妇av免费看| 久久久精品94久久精品| 另类精品久久| 日本-黄色视频高清免费观看| 五月玫瑰六月丁香| 亚洲av综合色区一区| 黄色一级大片看看| 久久人人爽人人爽人人片va| 久久久久久人人人人人| 成人亚洲欧美一区二区av| 日韩中文字幕视频在线看片| 在线观看www视频免费| 黄色配什么色好看| 国产亚洲欧美精品永久| 一级黄片播放器| 不卡视频在线观看欧美| 欧美精品国产亚洲| 久热久热在线精品观看| 国产有黄有色有爽视频| 在线观看国产h片| 18禁动态无遮挡网站| 免费高清在线观看日韩| 中文字幕最新亚洲高清| 老女人水多毛片| 边亲边吃奶的免费视频| 亚洲精品久久午夜乱码| 精品国产一区二区三区四区第35| 亚洲精品,欧美精品| 91精品国产国语对白视频| 精品亚洲成国产av| 亚洲精品乱码久久久久久按摩| 国产视频首页在线观看| 精品人妻一区二区三区麻豆| 这个男人来自地球电影免费观看 | 色哟哟·www| 激情视频va一区二区三区| 伦理电影大哥的女人| 777米奇影视久久| 国产精品国产三级国产专区5o| 国产探花极品一区二区| 十分钟在线观看高清视频www| 中文乱码字字幕精品一区二区三区| h视频一区二区三区| av播播在线观看一区| 9191精品国产免费久久| 午夜福利,免费看| 亚洲av成人精品一二三区| 女性被躁到高潮视频| 日韩成人伦理影院| 桃花免费在线播放| 婷婷成人精品国产| 天天躁夜夜躁狠狠久久av| 最近中文字幕高清免费大全6| 伊人久久国产一区二区| 欧美精品人与动牲交sv欧美| 好男人视频免费观看在线| 久久久国产欧美日韩av| 欧美精品一区二区免费开放| 日韩成人av中文字幕在线观看| 九色亚洲精品在线播放| 黄色配什么色好看| 久久精品国产自在天天线| 18禁在线无遮挡免费观看视频| 一级爰片在线观看| 观看av在线不卡| 亚洲欧美清纯卡通| 日韩,欧美,国产一区二区三区| 精品一区二区三区四区五区乱码 | 久久国产精品大桥未久av| 国产精品免费大片| 国产日韩欧美亚洲二区| 欧美精品一区二区免费开放| 99久久精品国产国产毛片| 中国美白少妇内射xxxbb| 国产一区二区在线观看av| 欧美精品亚洲一区二区| 成人亚洲欧美一区二区av| 久久ye,这里只有精品| 久久精品国产鲁丝片午夜精品| 一级毛片黄色毛片免费观看视频| 国产高清国产精品国产三级| 国产色婷婷99| 成人国产麻豆网| 最近2019中文字幕mv第一页| 在线天堂最新版资源| 久久久久人妻精品一区果冻| 日日啪夜夜爽| 中文字幕人妻熟女乱码| 人妻系列 视频| 夫妻性生交免费视频一级片| 亚洲精品久久午夜乱码| 免费看不卡的av| 日韩不卡一区二区三区视频在线| 欧美xxxx性猛交bbbb| videosex国产| 久久久精品94久久精品| 国产又爽黄色视频| 国产精品秋霞免费鲁丝片| 97人妻天天添夜夜摸| 亚洲中文av在线| 久久久精品94久久精品| 夜夜爽夜夜爽视频| 精品久久久久久电影网| 国产一级毛片在线| 精品福利永久在线观看| 大香蕉97超碰在线| 免费高清在线观看日韩| 国产日韩一区二区三区精品不卡| 一个人免费看片子| 少妇人妻 视频| 一级黄片播放器| 国产精品久久久久久久电影| 黄色毛片三级朝国网站| 99国产精品免费福利视频| 国产一区亚洲一区在线观看| 亚洲精品一区蜜桃| 观看av在线不卡| 日韩 亚洲 欧美在线| videosex国产| 国产av国产精品国产| 一区二区日韩欧美中文字幕 | 亚洲精品第二区| 全区人妻精品视频| 久久韩国三级中文字幕| 亚洲精华国产精华液的使用体验| 我的女老师完整版在线观看| 欧美日韩亚洲高清精品| 黄片无遮挡物在线观看| 最近的中文字幕免费完整| 美女脱内裤让男人舔精品视频| 搡老乐熟女国产| 久久久久久久亚洲中文字幕| 看非洲黑人一级黄片| 久久人人爽人人片av| av.在线天堂| 久久久久视频综合| 欧美亚洲日本最大视频资源| 免费高清在线观看视频在线观看| 9191精品国产免费久久| 免费大片黄手机在线观看| 男女高潮啪啪啪动态图| 秋霞伦理黄片| 2018国产大陆天天弄谢| av在线老鸭窝| 视频区图区小说| 韩国高清视频一区二区三区| 最新的欧美精品一区二区| 中文字幕精品免费在线观看视频 | 久久精品国产综合久久久 | 好男人视频免费观看在线| 深夜精品福利| 黑人欧美特级aaaaaa片| 久久久欧美国产精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲美女搞黄在线观看| 欧美日韩精品成人综合77777| 日本欧美国产在线视频| 国产av精品麻豆| av国产精品久久久久影院| 欧美变态另类bdsm刘玥| 免费av中文字幕在线| av在线老鸭窝| 一本色道久久久久久精品综合| 蜜桃国产av成人99| 如日韩欧美国产精品一区二区三区| 我要看黄色一级片免费的| 日韩熟女老妇一区二区性免费视频| √禁漫天堂资源中文www| 国产高清三级在线| 亚洲精品乱久久久久久| 亚洲图色成人| 亚洲av成人精品一二三区| 婷婷成人精品国产| 亚洲熟女精品中文字幕| 欧美丝袜亚洲另类| 91aial.com中文字幕在线观看| 久久久久网色| 亚洲久久久国产精品| 国产精品人妻久久久影院| 男人舔女人的私密视频| 日韩中字成人| 女人被躁到高潮嗷嗷叫费观| 日本色播在线视频| 亚洲成人一二三区av| 18在线观看网站| 国产 精品1| av.在线天堂| 黑人猛操日本美女一级片| tube8黄色片| 亚洲伊人久久精品综合| 亚洲图色成人| 亚洲国产av新网站| av黄色大香蕉| 色吧在线观看| 国产精品久久久av美女十八| 各种免费的搞黄视频| 国产深夜福利视频在线观看| 亚洲国产av影院在线观看| 女的被弄到高潮叫床怎么办| 观看av在线不卡| 在线观看www视频免费| 亚洲精品日韩在线中文字幕| 久久人人爽人人片av| 天天躁夜夜躁狠狠久久av| 国产老妇伦熟女老妇高清| 99视频精品全部免费 在线| 国产免费视频播放在线视频| 五月开心婷婷网| 久久97久久精品| 人人妻人人爽人人添夜夜欢视频| www.熟女人妻精品国产 | 久久ye,这里只有精品| 人妻 亚洲 视频| 欧美精品国产亚洲| 国产淫语在线视频| 免费av中文字幕在线| 秋霞伦理黄片| 国产亚洲最大av| h视频一区二区三区| 一级片'在线观看视频| 欧美精品高潮呻吟av久久| 国产精品国产三级国产av玫瑰| 18禁动态无遮挡网站| 男人操女人黄网站| 欧美成人精品欧美一级黄| 亚洲色图 男人天堂 中文字幕 | 成年美女黄网站色视频大全免费| 丰满少妇做爰视频| 精品熟女少妇av免费看| 五月天丁香电影| 满18在线观看网站| 九色亚洲精品在线播放| 最新的欧美精品一区二区| 在线免费观看不下载黄p国产| 亚洲一码二码三码区别大吗| 9热在线视频观看99| 天天操日日干夜夜撸| 老司机亚洲免费影院| 久久人妻熟女aⅴ| 亚洲四区av| 免费观看av网站的网址| 国产精品熟女久久久久浪| 亚洲av日韩在线播放| 日本欧美国产在线视频| 欧美精品人与动牲交sv欧美| 一级黄片播放器| 欧美亚洲 丝袜 人妻 在线| 黄网站色视频无遮挡免费观看| 男女啪啪激烈高潮av片| 麻豆精品久久久久久蜜桃| 一区二区三区精品91| 久久久精品区二区三区| 最黄视频免费看| 国产爽快片一区二区三区| 久久精品国产自在天天线| 天天躁夜夜躁狠狠躁躁| 我的女老师完整版在线观看| 夜夜爽夜夜爽视频| 91久久精品国产一区二区三区| 秋霞伦理黄片| 亚洲国产精品一区三区| 在线天堂中文资源库| 91在线精品国自产拍蜜月| 亚洲精品美女久久久久99蜜臀 | 久久久久久久久久人人人人人人| av卡一久久| 国产男女超爽视频在线观看| 又黄又粗又硬又大视频| 天天操日日干夜夜撸| 狂野欧美激情性bbbbbb| 街头女战士在线观看网站| 国产成人欧美| 国产精品秋霞免费鲁丝片| av在线app专区| 免费不卡的大黄色大毛片视频在线观看| 亚洲图色成人| 日韩三级伦理在线观看| 日韩电影二区| 欧美性感艳星| 女人久久www免费人成看片| 亚洲高清免费不卡视频| 狠狠婷婷综合久久久久久88av| 一区二区三区乱码不卡18| 天堂中文最新版在线下载| 久久久久网色| 欧美日韩成人在线一区二区| 欧美日韩亚洲高清精品| 精品国产一区二区久久| 国产高清三级在线| 久久久精品94久久精品| 亚洲久久久国产精品| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频 | 男人舔女人的私密视频| 国产精品国产三级国产av玫瑰| 91成人精品电影| 啦啦啦在线观看免费高清www| 亚洲国产毛片av蜜桃av| av一本久久久久| 欧美+日韩+精品| 亚洲精品av麻豆狂野| 亚洲av电影在线进入| 国产亚洲精品久久久com| 成人毛片a级毛片在线播放| 秋霞在线观看毛片| 日韩伦理黄色片| 亚洲av在线观看美女高潮| 国产欧美另类精品又又久久亚洲欧美| 巨乳人妻的诱惑在线观看| av一本久久久久| 亚洲久久久国产精品| av卡一久久| 精品熟女少妇av免费看| 亚洲精品成人av观看孕妇| 国产 一区精品| 久久人人97超碰香蕉20202| 七月丁香在线播放| 久久久精品免费免费高清| 在线观看美女被高潮喷水网站| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久电影网| 丁香六月天网| 男女国产视频网站| 亚洲性久久影院| 久久国产亚洲av麻豆专区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品成人久久小说| 免费高清在线观看视频在线观看| 黄色 视频免费看| 日本vs欧美在线观看视频| 男人添女人高潮全过程视频| 久久人人爽av亚洲精品天堂| 久久国产精品大桥未久av| 久久久久网色| 日本av免费视频播放| 老熟女久久久| 精品久久国产蜜桃| 精品人妻熟女毛片av久久网站| 高清不卡的av网站| 天堂中文最新版在线下载| 免费黄色在线免费观看| 国产欧美日韩综合在线一区二区| 国产深夜福利视频在线观看| 国产精品免费大片| 久久鲁丝午夜福利片| 国产免费一区二区三区四区乱码| 全区人妻精品视频| 欧美成人午夜精品| 精品一区二区三区四区五区乱码 | 国产精品一国产av| 亚洲精品一区蜜桃| 欧美另类一区| 精品酒店卫生间| 亚洲精品国产av成人精品| 久久久亚洲精品成人影院| 国产精品久久久久久久久免| 一本色道久久久久久精品综合| 成年女人在线观看亚洲视频| 久久青草综合色| 你懂的网址亚洲精品在线观看| 免费日韩欧美在线观看| 免费人妻精品一区二区三区视频| 三上悠亚av全集在线观看| 五月玫瑰六月丁香| 亚洲久久久国产精品| 人妻人人澡人人爽人人| 日韩av在线免费看完整版不卡| 欧美丝袜亚洲另类| av免费在线看不卡| 另类亚洲欧美激情| 午夜免费鲁丝| 欧美人与善性xxx| 亚洲色图综合在线观看| 日韩,欧美,国产一区二区三区| 国产精品女同一区二区软件| 国产又色又爽无遮挡免| 男女高潮啪啪啪动态图| 午夜福利乱码中文字幕| 精品久久国产蜜桃| h视频一区二区三区| 国产精品免费大片| 日韩一区二区视频免费看| 亚洲欧洲国产日韩| 日韩制服丝袜自拍偷拍| 少妇被粗大猛烈的视频| 亚洲人与动物交配视频| 韩国av在线不卡| 久久久久久久久久成人| 伦理电影免费视频| 另类亚洲欧美激情| 人人妻人人添人人爽欧美一区卜| 国产欧美亚洲国产| 久久久久精品性色| 在线观看免费日韩欧美大片| 午夜福利视频在线观看免费| 青青草视频在线视频观看| 一级毛片黄色毛片免费观看视频| 国产 精品1| 成人国产麻豆网| 一区二区日韩欧美中文字幕 | 少妇人妻精品综合一区二区| 搡女人真爽免费视频火全软件| 久久久a久久爽久久v久久| 国产日韩一区二区三区精品不卡| 90打野战视频偷拍视频| 99热全是精品| 我的女老师完整版在线观看| 深夜精品福利| 欧美日韩国产mv在线观看视频| 久久久久久久久久久久大奶| 18禁观看日本| 国产亚洲午夜精品一区二区久久| 久久久久精品人妻al黑| 男人爽女人下面视频在线观看| 欧美日韩成人在线一区二区| tube8黄色片| 精品少妇久久久久久888优播| 亚洲九九香蕉| 动漫黄色视频在线观看| 99久久综合精品五月天人人| 亚洲熟妇熟女久久| 美女国产高潮福利片在线看| 视频区图区小说| 黄色女人牲交| 午夜91福利影院| 久久精品成人免费网站| 婷婷精品国产亚洲av在线 | 在线天堂中文资源库| 欧美 亚洲 国产 日韩一| 亚洲熟女精品中文字幕| 啦啦啦免费观看视频1| 精品国产乱子伦一区二区三区| 免费观看人在逋| 身体一侧抽搐| 一级片'在线观看视频| 国产精品九九99| 久久国产精品影院| 亚洲一区高清亚洲精品| 亚洲av成人av| 国产精品98久久久久久宅男小说| 一进一出抽搐动态| 国产一卡二卡三卡精品| 欧美 亚洲 国产 日韩一| 国产熟女午夜一区二区三区| 精品一品国产午夜福利视频| 成人精品一区二区免费| 桃红色精品国产亚洲av| 搡老岳熟女国产| 一级毛片精品| 在线观看午夜福利视频| tube8黄色片| 国精品久久久久久国模美| 国产精品久久视频播放| 无遮挡黄片免费观看| 国产激情欧美一区二区| 在线观看午夜福利视频| 国产精品久久久av美女十八| 免费一级毛片在线播放高清视频 | 男女之事视频高清在线观看| 国产欧美日韩一区二区三区在线| 在线av久久热| 在线天堂中文资源库| 麻豆成人av在线观看| 亚洲精品乱久久久久久| 在线天堂中文资源库| 亚洲精品成人av观看孕妇| 18禁裸乳无遮挡动漫免费视频| 深夜精品福利| 777久久人妻少妇嫩草av网站| 99久久国产精品久久久| 午夜福利,免费看| 成年人午夜在线观看视频| 精品一区二区三区四区五区乱码| 久久精品国产清高在天天线| av不卡在线播放| 一进一出抽搐动态| 99热只有精品国产| 成年版毛片免费区| 99热网站在线观看| 在线天堂中文资源库| 国产亚洲av高清不卡| 成人三级做爰电影| 国产精品自产拍在线观看55亚洲 | e午夜精品久久久久久久| 亚洲av成人av| 久热爱精品视频在线9|