999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

水氮互作對陜南烤煙煙堿合成及代謝的影響

2025-08-09 00:00:00郝蘇娟張傲杰王平平張明誠艾綏龍張立新
中國種業 2025年6期
關鍵詞:水氮煙堿氮量

Effects of Water and Nitrogen Interaction on Nicotine Synthesis and Metabolism of Flue-Cured Tobacco in Southern Shaanxi

HAO Sujuan',ZHANG Aojie2, WANG Pingping^,ZHANG Mingcheng2,AI Suilong2,ZHANG Lixin2 ('Lianhu DistrictMeteorologicalServiceCenter,Xi'an7o;CollegeofLifeScience,NorhwestAamp;FUniversity, Xianyang712ioo,Shaanxi;YantaDistrict InstituteofTobaccoScience,Xi'an710000)

煙草是我國的主要葉用經濟作物之一,是卷煙生產的主要原料,還可用于提取尼古丁、木質素等,具有較高的商業價值[1]。此外,煙草在國家稅收中有著舉足輕重的地位,同時在精準扶貧方面發揮著重要作用。烤煙生產對提高農民收人、推動地方經濟發展有重要意義。

水分和肥料是影響煙草生長的最為重要的兩大栽培因素。水分是肥料發揮效用、被作物利用的基礎,肥料是水土系統中影響作物生產效率的關鍵[3]。氮是作物生產的主要驅動因素之一,與土壤水分條件有明顯互作[4-5]。施氮和灌水在提高烤煙水分利用率上存在顯著的互作效應,在土壤水分不足時,施加氮肥能夠促進烤煙對水分的利用,而在施氮量較多時,灌水又能促進烤煙對肥料的吸收利用。水氮之間的互作效應顯著促進了烤煙干物質的積累,影響烤煙株高、莖圍、葉面積等農藝性狀,從而提高煙葉產量。 15N 穩定同位素試驗表明,水氮互作可以顯著提高烤煙的氮肥利用率[。合理的水氮組合能協調烤煙各種化學成分含量,提高煙葉品質。此外,施氮量增加可以提升烤煙的氧化防御系統,減輕水分脅迫下烤煙的氧化損傷[7]。

代謝物是細胞調控終產物,參與多種生化反應及生理調控[,最接近生物的表型或者生理狀態,其種類、含量變化可以看作是生物體對表型或基因變化響應的結果,更能反映生物的基因、環境與表型三者之間的關系。代謝組學應運而生,主要通過高通量檢測對生物體在限定條件下的代謝物進行定性、定量分析,探究其代謝產物差異及代謝途徑的變化規律來揭示內在機制[9-10]。代謝組學作為系統生物學新的分支,目前已經成為繼基因組學、轉錄組學和蛋白質組學之后新的研究熱點。Zhao等[運用代謝組學方法分析不同產地來源的烤煙代謝物差異,發現不同產地的烤煙代謝物呈現明顯的分離,不同生態區域可影響烤煙氮代謝等多種代謝途徑,進一步篩選出多種與不同產地煙葉香味特征有關的差異代謝產物,為煙葉香型的研究提供了理論依據。隨著相關技術的發展,在玉米、水稻等作物上通過代謝組學方法進行機理探索的研究逐步深入。代謝組學技術在煙草研究中的應用也逐步展開[12-15],成為全面、系統研究煙草對環境或外界刺激響應機理的重要手段[12]。

本研究針對陜南煙區烤煙生產中的水氮供應匹配性,以云煙99為材料,探究水氮互作對烤煙煙堿含量、煙堿合成過程中關鍵酶活性的影響,同時采用代謝組學方法探究不同水氮處理影響煙堿合成的內在機制,旨在為烤煙產質量調控提供參考。

1材料與方法

1.1試驗地概況試驗于2021年4-8月在西北農林科技大學北校區曹新莊試驗站( 34°1723′′N 108°412′′E ,海拔 520m )日光大棚內完成。種植土壤為陜西省商洛市洛南縣植煙黃棕壤土,基本理化性質: pH 值7.22,有機質含量 10.25g/kg ,速效氮含量 30.25mg/kg ,速效磷含量 34.25mg/kg ,速效鉀含量 127.48mg/kg 。

1.2試驗設計烤煙品種為云煙99,由陜西省煙草科學研究所提供。本試驗涉及水、氮2個因素,其中土壤水分條件3個水平:W1(土壤相對含水量75%~85% )W2(土壤相對含水量 55%~65% )、W3(土壤相對含水量 35%~45% ),施氮量3個水平:N1( 4g/ 株)N2( 5g/ 株)N3( 6g/ 株),雙因素完全組合設計共計9個處理(表1。煙苗成苗后移入桶中,每桶定植1株,移栽前保持桶中土壤相對含水量60% ,移栽后36d進人旺長期,開始進行水分處理。采用稱重法補水,控制土壤相對含水量,控水處理進行 30d 。按照試驗處理,稱取相應克數的氮肥給每株煙苗進行基施。

表1試驗處理

1.3樣品采集于移栽后36d(旺長初期)51d(旺長中期)66d(旺長末期)85d(成熟期)分別對根部及葉片進行取樣。葉片于 105% 殺青 30min 65°C 烘干至恒重,研磨并過100目篩,用于煙堿含量測定。取根尖3cm迅速放入液氮,后保存于 -80% 冰箱,用于測定煙堿合成途徑中關鍵酶的活性、代謝產物含量。取樣方式為各處理隨機3株樣品混合,每個處理3次重復。

1.4測定項目與方法煙堿含量測定參考鄭璞帆[的方法并稍作改動,采用HCI浸提法。稱取粉碎過篩烘干的煙葉樣品 0.03g ,加人 0.5mol/L 鹽酸溶液 25mL ,活性炭1/3角匙,震蕩 10min 后過濾,定容至 100mL ,用紫外分光光度計測定 236nm 、259nm 和 282nm 波長的吸光度。

煙堿合成途徑中關鍵酶活性測定煙草根樣添加液氮后進行粉碎,取 0.1g 加人 0.05mol/L 磷酸緩沖液(pH值7.4) 1.5mL ,充分混勻后 4°C 下8000r/min 離心 20min ,上清液用于精氨酸脫羧酶(ADC)鳥氨酸脫羧酶(ODC)腐胺 N- 甲基轉移酶(PMT)的測定,其測定步驟按照上海通蔚生物科技有限公司的酶聯免疫試劑盒說明書進行。

代謝產物含量測定樣品前處理:煙草根樣添加液氮后進行粉碎,取根樣 0.2g ,加入 75% 甲醇溶液 400μL ,充分混勻后超聲提取 30min,4°C 下 14000r/min 離心 20min ,上清液用 0.22μm 有機相過濾器過濾后加入 2mL 進樣瓶中,用于高分辨離子淌度液質聯用儀HRLCMS( LC-30A+ TripleTOF5 5600+ )上機進行正、負離子模式分析。

數據庫匹配將提取的峰信息在MassBank、Respect、GNPS這3個庫進行全庫檢索與比對。數據庫匹配參數見表2。

表2數據庫匹配參數

1.5數據統計與分析采用Excel2013、IBMSPSSStatistics26和RStudio1.3.1093.0對數據進行整理、統計分析和做圖。

2 結果與分析

2.1水氮互作對烤煙煙堿含量的影響由表3可以看出,不同水氮處理的烤煙中部葉煙堿含量在不同時間段均逐漸上升。移栽后36d,各處理的中部葉煙堿含量差異不大,此后水分控制開始,各處理的中部葉煙堿含量差異逐漸增大,在移栽后66d各處理中部葉煙堿含量差異最為顯著。整體來看,同一時期相同氮肥處理下,中部葉煙堿含量隨著土壤相對含水量的增加而降低;W1N1和W1N2處理在整個過程中保持相對平穩的增長,在移栽后85d煙堿含量均較低;W3水分條件下,各時期中部葉煙堿含量均較高。同一時期相同水分處理下,中部葉煙堿含量隨著氮肥的增加而增加。

2.2水氮互作對煙堿合成途徑關鍵酶活性的影響

由圖1可知,移栽后36d即控水之前各處理的精氨酸脫羧酶(ADC)活性無顯著差異。移栽后51d,相同氮肥處理下ADC活性呈現隨土壤相對含水量降低而升高的趨勢,且W3處理下ADC活性達到W1條件下的 1.64~1.81 倍;W1水分條件下,W1N2處理的ADC活性略低于其余處理,但無顯著差異;W2、W3水分條件下,ADC活性隨施氮量增加而增加,W2N2、W2N3處理的ADC活性較W2N1處理分別提高 10.7%.21.1% ;W3N2、W3N3處理的ADC活性較W3N1處理分別提高 0.5%.0.8% ;因此,隨土壤相對含水量降低施氮量對ADC活性的提升效果減弱。移栽后66d,W2水分條件下ADC活性最高,且隨施氮量增長ADC活性提高;W3水分條件下,ADC活性較移栽后51d下降,且隨施氮量增加呈現先增加后降低的趨勢。移栽后85d,所有處理的ADC活性均處于較高水平,施氮量對ADC活性影響效果減弱。

由圖2可知,W1、W2水分條件下,各處理的鳥氨酸脫羧酶(ODC)活性隨移栽后時間增加而增加,W3水分條件下,各處理ODC活性呈現隨移栽后時間增加先增加后降低再增加的趨勢。移栽后36d,各處理的ODC活性無顯著差異;移栽后51d,W2、W3條件下ODC活性隨施氮量增加而升高,其中W2條件下增幅較高,W2N2、W2N3較W2N1分別提高了 17.6%.28.6% ;移栽后66d,W3水分條件下,各處理較移栽后51d分別降低了 2.1%.6.9%.2.5% 移栽后85d,W1水分條件下ODC活性仍隨施氮量增加而提高,W2、W3水分條件下則隨施氮量增加呈現先增加后降低的趨勢。

表3中部葉煙堿含量 (%)
同列不同小寫字母表示在0.05水平上存在顯著差異,下同
圖1水氮互作對精氨酸脫羧酶活性的影響

由圖3可知,所有處理的腐胺N-甲基轉移酶(PMT)活性均隨著移栽時間的增加而增加,移栽后36d各處理之間的PMT活性無顯著差異;進行水分控制后,各處理PMT活性出現差異,整體表現為隨土壤相對含水量的降低PMT活性升高。移栽后36\~66d,均以W3N3處理的PMT活性最高;移栽后85d,W1、W3水分條件下不同施氮量處理的PMT活性相近;W2水分條件下,W2N3處理的PMT活性最高,比W2N1提高了 17.6% 。

2.3水氮互作下烤煙代謝組學分析

2.3.1總體代謝物層次聚類分析對煙堿含量差距較大的W3N3和W1N1兩組處理進行代謝組學分析,由PCA主成分分析結果可知進樣過程中儀器運行較為穩定,可以用于后續的分析。對鑒定出的全部代謝物的相對定量值歸一轉化后進行層次聚類得到圖4,列為所有樣本,行為全部代謝物,紅色代表高表達代謝物,藍色代表低表達代謝物。由圖4可以看出,有部分代謝物高、低表達交互在一起,但高、低表達區域較明確。同一組別樣品的高、低表達區域較為相近,不同組別樣品間高、低表達區域差異明確。

圖2水氮互作對鳥氨酸脫羧酶活性的影響圖3水氮互作對腐胺 N- 甲基轉移酶活性的影響
2.3.2 差異代謝物分析對代謝產物以Fold
圖4總體代謝物層次聚類分析圖圖5代謝物火山圖

產物中包括10種有機酸,7種胺類,5種醇類,3種酚類,2種酮類,2種酯類,1種氨基酸類,1種堿類,其他類8種。

2.3.3代謝通路富集分析如圖6所示,差異代謝產物主要富集的代謝通路為:氨酰tRNA生物合成,ABC轉運器,氨基酸生物合成,嘌呤代謝,代謝途徑,咖啡因代謝,纈氨酸、亮氨酸和異亮氨酸生物合成,丙氨酸、天冬氨酸、谷氨酸代謝,谷胱甘肽代謝,單酰胺類生物合成代射通路。氣泡圖中,氣泡顏色越深代表該通路差異顯著性越高,直徑越大代表富集到的代謝物越多。因此,在富集到的通路中,氨酰tRNA生物合成途徑富集的代謝物較少,但差異顯著性較高;代謝途徑富集的代謝物最多。

2.3.4煙堿代謝相關產物將已識別的代謝產物中與煙堿合成相關的代謝產物進行聚類分析,結果如圖7,有部分代謝物高、低表達交互在一起,但有相對明確的高、低表達區域。

圖6代謝通路富集分析圖

a為柱狀圖,b為氣泡圖;Valine,leucineandisoleucinebiosynthesis為氨酸亮氨酸和異亮氨酸生物合成,Purinemetabolism為嘌呤代謝,Monobactam biosynthesis 為單酰胺類生物合成代謝通路,Metabolic pathways為代謝途徑,Glutathione metabolism為谷胱甘肽代謝,

Cafeine metabolism為咖啡因代謝,Biosynthesis ofaminoacids為氨基酸生物合成,Aminoacy1-tRNAbiosynthesis為氨酰tRNA生物合成,Alanine,aspartate and glutamate metabolism為丙氨酸、天冬氨酸和谷氨酸代謝,ABC transporter為ABC轉運器

對煙堿合成相關代謝產物含量進行相關性分析,結果如圖8。有9對代謝產物呈顯著正相關,有9對代謝產物呈顯著負相關。天冬氨酸、精氨酸含量與煙堿含量呈顯著正相關。初步判斷水氮處理對烤煙煙堿含量造成的影響可能是由于天冬氨酸、精氨酸含量的變化。

3討論

常規化學成分是判斷煙葉品質好壞的關鍵因素[17],煙堿是煙葉最重要的化學成分之一,煙堿是煙草特有的生物堿,影響烤煙的勁頭和刺激性,其含量高低是衡量烤煙品質的重要指標[18-22]。土壤水分條件會影響煙堿含量,本研究中中部葉煙堿含量表現為土壤相對含水量 35%~45% 最高, 75%~85%最低,這與汪耀富等[3研究結果一致。氮作為煙堿的重要組成成分,其施用量與煙堿含量密切相關。水分和氮肥在作物生長中存在相互促進、相互制約的關系。馬振奇[23]進行水氮耦合研究發現,不同水氮處理下中部葉的煙堿含量均呈現隨土壤相對含水量降低而升高,隨施氮量增加而升高的規律。本研究中中部葉煙堿含量整體變化規律與之一致。

圖7煙堿合成相關代謝產物熱圖圖8煙堿合成相關代謝產物含量相關性分析

張瑞霞[24]測定了不同滲透脅迫程度下24h內烤煙葉片精氨酸脫羧酶和鳥氨酸脫羧酶的酶活性變化,發現精氨酸脫羧酶活性隨脅迫程度的加劇而顯著升高,這與本研究一致;不同的是,在該研究中隨脅迫程度加劇、時間延長,烤煙鳥氨酸脫羧酶活性無顯著變化,而在本研究中移栽后36d、51d烤煙根部鳥氨酸脫羧酶活性受土壤相對含水量影響較小,在移栽后66d土壤相對含水量 35%~45% 條件下其活性降低,這種差異是否是由取樣部位不同、處理時長不同或者是由植株生長水平不同帶來的尚需進一步研究。本研究中,不同土壤水分條件下施氮量的增加不同程度上提高了精氨酸脫羧酶、鳥氨酸脫羧酶活性,但其作用機制尚不明確。黃化剛等[25]研究表明,同等用氮量下有機氮相比于無機氮可以提高鳥氨酸脫羧酶的表達量,導致烤煙上部葉煙堿含量增加。不同氮素形態也會影響精氨酸脫羧酶、鳥氨酸脫羧酶活性,從而導致煙堿積累量的差異[26]。烤煙根系腐胺N-甲基轉移酶活性與葉片煙堿含量呈顯著正相關[27-28]。

近年來,代謝組學不斷發展,為煙草的提質增香及代謝機理研究奠定了基礎。本研究中,精氨酸含量與煙堿含量呈正相關。對比分析發現,這種差異的形成可能是由于取樣部位不同造成,本研究所用材料為合成煙堿的根系部位,而向海英等[29]研究對象為烤煙葉片部位,根部與葉片部位的這種差異與煙堿合成、轉運之間的關系仍需進一步研究;此外,研究共識別出4086個含有二級碎片的一級質譜峰,鑒別出245個差異代謝產物。本研究中,對不同水氮處理后煙堿含量差異較高的兩個處理進行代謝物提取及組學分析,共鑒別出190種代謝物,這種差異可能是由于取樣部位、提取效率、儀器分辨率等不同導致。識別出的代謝產物數量及精準度是代謝組學精確分析的基礎,完善植物組織中代謝產物提取方法,提高代謝產物提取效率及比對精準度是代謝組學發展的重要內容。

4結論

施氮量對煙堿合成的作用效果因土壤相對含水量不同而有所差異,中部葉煙堿含量在土壤相對含水量 35%~45% 條件下最高,且以施氮量 6g/ 株處理煙堿含量最高,煙堿含量以W1N1處理最低。不同水分條件下適宜的施氮量可以提高煙堿合成關鍵酶活性,代謝物差異分析結果表明,差異代謝物在氨酰tRNA生物合成途徑、ABC轉運器通路、氨基酸生物合成途徑及嘌呤代謝中顯著富集;煙堿合成途徑代謝產物存在顯著差異,其中天冬氨酸、精氨酸含量與煙堿含量呈顯著正相關。

參考文獻

[1]陳雪嬌,李云飛,李剛,宗凱,鄭海松,孫娟娟,姚劍,余曉峰,徐業平.煙草主要病害及其相關檢測方法研究進展.河南農業科學,2013,42(11):1-6

[2]周繼洲.基于GIS的陜西省環秦嶺生態煙區煙草生態適宜性區劃.楊凌:西北農林科技大學,2014

[3]汪耀富,楊天旭,孫德梅,葉紅潮.灌水及氮素形態對烤煙營養元素含量及產量品質的影響.河南農業大學學報,2006,40(5):477-481

[4]ZhaoSM,ZhengH,LuX,LiuY,SuBL,XuGW.Metabonomicsand phospholipid metabolic profilingof abnormal glucose metabolismbasedon high performance liquidchromatography-electrospraymassspectrometry.Chinese Journal of Chromatography,2011,29(4):307-313

[5]LiQH,ZhaoCX,LiY,ChangYW,WuZM,PangT,LuX,WuY,XuGW.Liquid chromatography/mass spectrometry-based metabolicprofiling to elucidate chemical differences of tobacco leaves betweenZimbabwe and China.Journal of Separation Science,2010,34(2):119-126

[6]袁仕豪.水氮互作對烤煙氮素吸收利用及煙葉產量和品質的影響.鄭州:河南農業大學,2008

[7]崔保偉,陸引罡,張振中,任琨.不同生育期水分脅迫對烤煙生理特性及化學品質的影響.中國煙草科學,2009,30(3):19-23

[8]高婭北,李成剛,陳二龍,張明剛,孔德輝,孫占偉,史龍飛,宋朝鵬.烤煙煙葉類黃酮合成關鍵酶活性及其調控基因的表達分析.中國煙草學報,2019,25(1):86-92

[9]謝良文,路曉崇,裴曉東,蔣博文,李生棟,魏碩,李帆.烤煙烘烤過程中氨基酸含量變化對主要含氮化合物的影響.湖南農業大學學報:自然科學版,2016,42(2):152-156

[10]PiornosJA,BalagiannisDP,Methven L,KoussissiE,BrouwerE,ParkerJK.Elucidating the odor-activearoma compoundsinalcohol-freebeerand theircontributiontotheworty flavor.JournalofAgricultural andFoodChemistry,2020,68(37):10088-10096

[11]ZhaoJY,HuCX,ZengJ,ZhaoYN,ZhangJJ,ChangYW,LiLL,ZhaoCX,LuX,XuGW.Studyofpolarmetabolitesintobacco from different geographical originsby usingcapillaryelectrophoresis-mass spectrometry.Metabolomics,2014,10(5):805-815

[12]王小莉,付博,趙銘欽,賀凡,王鵬澤,劉鵬飛.代謝組學技術在煙草研究中的應用進展.中國煙草科學,2016,37(1):89-96

[13]梁珂馨,奚家勤,陳若星,牟文君,劉萍萍,鄭慶霞,王得強,賈中林,李曉闖,張艷玲.基于非靶向代謝組學的片煙醇化進程判斷及代謝物差異分析.煙草科技,2023,56(4):17-25

[14]虞昕磊,何結望,林國平,李金海,王大愛,袁躍斌,劉圣高,李志豪,陶德欣.夏冬兩季發酵雪茄煙葉的代謝組差異分析.生物技術通報,2024,40(6):260-270

[15]陸岸芷,張樂毅,王平平,陳新穎,杜姍姍,李先鋒,王輝,張強,艾綏龍,張立新.半干旱區2種雪茄煙葉發酵后差異代謝物分析.干旱地區農業研究,2024,42(4):22-28

[16]鄭璞帆.陜南地區不同成熟度煙葉生理特性和烤后品質研究.楊凌:西北農林科技大學,2017

[17]康崢,唐衛迪,楊繼利,曾皓,趙曉軍,李俊花,孔云鑫,張廷金,鄧綱,張海平.不同生物菌劑對根結線蟲病害地區煙草生長及品質的影響.中國種業,2024(8):79-86

[18]鄧金華.三種有機物料配施對土壤水分特性及烤煙生理與品質的影響.楊凌:西北農林科技大學,2022

[19]雷亞芳.不同酶制劑與香料處理對煙葉品質及微生物多樣性的影響.楊凌:西北農林科技大學,2021

[20]向羽,周成,黃楚天,何輝,柳洪入.結構優化對不同葉位煙葉品質及揮發性物質含量的影響.安徽農業科學,2024,52(10):15-20

[21]吳鑫穎.煙葉醇化過程的動態解析及功能微生物研究.無錫:江南大學,2022

[22]陳愛國,劉光亮,周道金,謝鳳標,周培祿,蔣代兵,周東新,林天然,沈燦柱,曾文龍.龍巖鮮煙葉特征香氣成分前體物代謝組學的通路分析.安徽農業大學學報,2018,45(3):545-550

[23]馬振奇.水氮耦合對烤煙生長發育以及土壤養分的影響.鄭州:河南農業大學,2016

[24]張瑞霞.烤煙激素與多胺對水分迫和氮肥形態的響應及其調節研究.鄭州:河南農業大學,2004

[25]黃化剛,申燕,吳飛躍,楊振智,劉世碧,喻奇偉,賈宏昉.不同氮素來源對烤煙上部葉碳氮代謝及其基因表達的影響.南方農業學報,2018,49(3):462-468

[26]王學龍.礦質營養對烤煙煙堿的調控及其生理機制研究.武漢:華中農業大學,2012

[27]蒲媛媛.煙堿合成途徑中關鍵酶活性對雜種優勢形成的影響.貴陽:貴州大學,2020

[28]MizusakiS,TanabeY,NoguchiM,TamakiE.Changesintheactivitiesof ornithine decarboxylase,putrescine N-methyltransferaseandN-methylputrescine oxidase in tobaccorootsin relationtonicotinebiosynthesis.PlantandCellPhysiology,1973,14(1):103-110

[29]向海英,劉欣,逢濤,曾婉俐,許力,高茜,鄧樂樂,張建鐸,李雪梅,李勇,孔維松.基于非靶向代謝組學的紅花大金元不同生育期煙葉代謝圖譜差異分析.中國煙草學報,2022,28(4):77-84

(收稿日期:2025-02-19)

猜你喜歡
水氮煙堿氮量
煙草品質構成要素研究進展
機器學習算法在煙葉烘烤過程化學成分估測中的應用
煙葉原料組合加工工藝研究
施氮水平對植絲式混合草坪生長及氮磷養分去向的影響
草地學報(2025年7期)2025-08-03 00:00:00
不同施氮量對大寧曬煙生長發育和品質的影響
不同水氮處理對超晚播冬麥生長、水氮利用及產量的影響
不同施氮量對節水抗旱稻產量及氮肥利用效率的影響
北方水稻(2025年4期)2025-07-30 00:00:00
水氮耦合對水稻生長發育、碳氮代謝和氮吸收利用的調控研究進展
北方水稻(2025年4期)2025-07-30 00:00:00
水氮調控對滴灌香梨生長及產量品質的影響
不同水氮調控對綠洲灌區膜下滴灌西瓜生長及產量的影響
中國瓜菜(2025年6期)2025-07-22 00:00:00
主站蜘蛛池模板: 天天做天天爱夜夜爽毛片毛片| 精品国产成人a在线观看| 国产男女免费完整版视频| 国产精品无码一二三视频| 亚洲免费福利视频| 福利小视频在线播放| 国产一级无码不卡视频| 亚洲日韩国产精品综合在线观看| 四虎永久免费在线| 亚洲精品在线91| 国产迷奸在线看| 秋霞午夜国产精品成人片| 欧美日韩激情在线| 亚洲国产日韩一区| 尤物成AV人片在线观看| 亚洲自偷自拍另类小说| 午夜国产不卡在线观看视频| 亚洲福利网址| 国产成人乱无码视频| 片在线无码观看| 国产美女91视频| 四虎国产永久在线观看| 免费不卡视频| 广东一级毛片| a级毛片免费看| 国产农村妇女精品一二区| 日本在线欧美在线| 中国一级特黄视频| 日韩123欧美字幕| 久久精品无码专区免费| 国产视频你懂得| 99久久精品免费观看国产| 欧美一级在线看| 99久久国产综合精品女同| 国产交换配偶在线视频| 看你懂的巨臀中文字幕一区二区| 国内精品伊人久久久久7777人| 亚洲视频在线青青| 波多野结衣无码视频在线观看| 国产精品一区二区在线播放| 亚洲国产成人无码AV在线影院L| 国产精品大白天新婚身材| 日本在线亚洲| 露脸真实国语乱在线观看| 91年精品国产福利线观看久久 | 中文字幕自拍偷拍| 国产亚洲精品在天天在线麻豆| 欧美不卡视频在线观看| 日韩一级二级三级| 欧美国产日韩在线| 黄色网站不卡无码| 中文字幕在线永久在线视频2020| 91免费国产在线观看尤物| 亚洲日本中文综合在线| 中文字幕亚洲乱码熟女1区2区| 色婷婷亚洲综合五月| 亚洲无码视频一区二区三区| 久久夜色撩人精品国产| 在线无码av一区二区三区| 日本国产精品一区久久久| 四虎亚洲国产成人久久精品| 国产成人免费观看在线视频| 欧美日韩精品综合在线一区| 国产尤物视频网址导航| 91香蕉视频下载网站| 日本不卡在线视频| 亚洲av无码久久无遮挡| 男人天堂亚洲天堂| 亚洲精品色AV无码看| 欧美午夜网站| 欧美日韩亚洲综合在线观看 | 日韩精品无码免费一区二区三区| 1024你懂的国产精品| 午夜毛片福利| 69国产精品视频免费| 亚洲国产精品无码AV| 亚洲色成人www在线观看| 欧美一级特黄aaaaaa在线看片| аⅴ资源中文在线天堂| 精品精品国产高清A毛片| 99久久人妻精品免费二区| 综1合AV在线播放|