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

TTC染色指數配合Logistic方程鑒定山葡萄種質抗寒性

2018-06-21 09:29:52楊義明王振興劉迎雪劉海雙
農業工程學報 2018年11期

趙 瀅,艾 軍,楊義明,王振興,劉迎雪,何 偉,劉海雙

(中國農業科學院特產研究所,長春 130112)

0 引 言

山葡萄(Vitis amurensis Rupr.)屬多年生木質藤本,是葡萄屬中最抗寒的一個種[1-2]。休眠期枝條的抗寒性具有較高的遺傳力,是葡萄抗寒育種的重要親本材料[3-4]。但田間觀察發現不同山葡萄種質抗寒能力存在差異,而目前適于山葡萄休眠期枝條抗寒性的鑒定方法尚不夠完善,對不同山葡萄種質資源抗寒差異的了解仍不清晰[5],影響了抗寒種質的進一步發掘與利用。

低溫半致死溫度(LT50)可用來反映植物的抗寒性,實驗室常用相對電導率法擬合計算低溫半致死溫度來鑒定葡萄屬的抗寒性[6-7]。然而,在實際測定中發現由于溫度、浸泡時間以及空氣中 CO2的溶解等因素都會對實際測定的電導值產生影響,造成測定結果重復性較差[8-10]。而采用隸屬函數法對作物抗寒性進行綜合評價的結果較為全面、準確[11-14],但對于大量種質抗寒性的篩選則需要尋找更為簡便且準確的鑒定方法。

2,3,5-氯化三苯基四氮唑(TTC)染色法也可用于作物抗寒性鑒定[15-18],具有直觀且操作簡便等優點,但目前關于TTC染色法測定枝條組織活力的研究報道中,均是通過觀察枝條橫截面剪口著色深淺程度進行判定[19]。在實際操作中發現該方法存在一些不足,主要是染色生成紅色的三苯基甲臢(TTF)經常在失去組織活性的枝條橫截面剪口外側堆積,造成觀測誤差;另外,枝條橫切面著色深淺程度目測也難以準確判斷,需量化。

本研究將以40份山葡萄種質1 a生休眠期枝條為試材,對浸染枝段進行縱切,以枝段縱切面染色面積為判斷標準,可避開原有TTC染色法橫截面上出現TTF堆積問題;同時探索TTC染色圖像的可視化評估,以及量化獲得的染色指數能否配合 Logistic方程在山葡萄抗寒性鑒定上進行應用,并與隸屬函數法獲得的綜合抗寒性評價結果及田間鑒定結果進行比較,以此來驗證這種以枝段縱切面染色面積為判定標準的抗寒性鑒定新方法的可行性和準確性,以期建立起適于大量山葡萄種質休眠期抗寒性鑒定的簡便、可靠方法,從而判斷不同山葡萄種質抗寒性強弱,最終為山葡萄抗寒評價研究及抗寒育種材料選擇提供有效鑒定方法和參考依據。

1 材料與方法

1.1 試驗材料

供試材料采自國家果樹種質左家山葡萄圃,種質圃地處吉林省吉林市左家鎮,位于北緯 44°04′,東經 126°05′,海拔190 m。年均氣溫3.6 ℃,絕對低溫-39.8 ℃。圃地為南向緩坡地,土壤為暗棕森林土。圃中每份資源為扦插繁殖的自根苗,籬架生長,株行距1.0 m × 2.5 m,柱間距6 m。樹體肥水管理采用常規管理模式[20]。40份試材包括純種山葡萄種質(V. amurensis)38份,山歐雜種((V.amurensis×V.vinifera)×V.amurensis)2 份。38 份山葡萄種質的休眠期枝條于2017年1月剪取,當年1月最低氣溫為-28.2 ℃;山歐雜種‘北冰紅’和‘左優紅’為對照品種,其休眠枝條于深秋10月底未埋土防寒前剪取,然后放入冷庫中待處理。所用枝條均為長勢和粗度相對一致(約0.8 cm粗)且健康成熟的1年生枝條,架面南側1.0~1.8 m范圍內取樣。每份種質取50根枝條,每根60 cm長。

1.2 材料處理

材料處理參照 Zhang等的方法[21],略有修改。將采集的枝條帶回實驗室,剪成20 cm長的小段,裝于自封袋中,-15 ℃冰箱中貯藏。使用高低溫交變試驗箱(型號:GDJW-225,北京雅士林試驗設備有限公司)對枝條進行低溫處理,處理溫度為-20、-25、-30、-35、-40、-45、-50 ℃,按4.0 ℃/h的速度降溫,至設定溫度后保持12 h,再按4.0 ℃/h的速度升至20 ℃后,取出靜置2 h。去除表皮,避開芽眼,以備用于抗寒性評價及指標測定等研究。

1.3 抗寒性的實驗室鑒定

1.3.1 TTC染色圖像可視化評估鑒定枝條抗寒性

不同低溫處理后,從各枝條上剪取0.5 cm小段,每份種質取10段,放入10 mL 0.5 % TTC染液中,于30 ℃恒溫培養箱中避光染色4 d后取出,然后將枝段縱切,蒸餾水沖洗 3次,濾紙吸干水分,放入連有 WinRHIZOTM圖像分析軟件的 LA2400掃描儀中對枝段縱切面進行掃描;利用圖像分析軟件計算縱切面染色面積并確定染色等級;同時根據染色等級計算染色指數,以染色指數反映各試材的受凍程度。

染色分級標準如下(圖1):4級,縱切面90 %以上的面積呈紅色;3級,縱切面70%~90%的面積呈紅色;2級,縱切面40%~70%的面積呈紅色;1級,縱切面5%~40%的面積呈紅色;0級,縱切面5%以下的面積呈紅色。染色指數= Σ(染色級別×枝段數)/(4×枝段總數)。

圖1 染色分級標準Fig.1 Standard of staining identification

根據不同山葡萄種質各低溫處理下枝段染色指數的變化情況,配合 Logistic方程建立回歸方程,擬合溫度拐點來確定其低溫半致死溫度(LT50),具體參照解越等方法[22]。

1.3.2 隸屬函數法綜合評價枝條抗寒性

采用隸屬函數法對山葡萄種質抗寒性進行綜合評價。根據最適冷凍溫度下枝條 7項抗寒指標(包括相對電導率、丙二醛、脯氨酸、可溶性糖、可溶性蛋白、萌枝率、凍害指數)的平均隸屬度值將山葡萄種質抗寒性分為5個級別,分級標準參照張劍俠等的方法[5],Ⅴ級,平均隸屬度值> 0.70~1.00;Ⅳ級> 0.60~0.70;Ⅲ級>0.40~0.60;Ⅱ級> 0.30~0.40;Ⅰ級 0~0.30。

1) 生理指標測定

相對電導率測定參照 Verslues等[23]的方法。丙二醛含量測定采用硫代巴比妥酸法[24];脯氨酸含量測定采用茚三酮比色法[25];可溶性糖含量測定采用蒽酮比色法[26];可溶性蛋白含量測定采用考馬斯亮藍G-250染色法[27]。

2) 水培萌枝率調查

每份試材選取帶芽眼的枝段共30根,10根為1次重復,共 3次重復。經-35 ℃最適冷凍溫度處理(處理方法同1.2)后插入水中,1個月后調查水培萌枝率,參照Nanami等[28]方法。培養期間每 5 d換一次水。萌枝率=萌芽枝數/總枝數×100 %。

3) 凍害指數調查

枝段經最適冷凍溫度處理后,觀察橫切面組織在低溫處理后的褐變情況來確定凍害級別,凍害級別及冷害指數計算參照呂中偉等的方法[29]。

1.4 抗寒性的田間鑒定

通過調查每份試材的田間萌芽率來對其田間抗寒性進行鑒定。田間萌芽率調查參照劉崇懷等[30]方法,本試驗于2017年5月初調查每份試材在資源圃中經自然越冬后的萌芽情況,統計萌芽率。萌芽率=單株芽眼萌芽數/單株總芽眼數×100%。

1.5 數據統計分析

采用SAS9.0軟件對數據進行方程及相關性分析,利用 Duncan多重比較進行檢驗;作圖采用 Sigmaplot10.0軟件。

2 結果與分析

2.1 TTC染色圖像可視化評估鑒定山葡萄種質抗寒性

山葡萄經不同低溫處理后染色指數與冷凍溫度呈典型的倒‘S’型曲線,說明可以利用非線性回歸分析結合Logistic方程擬合溫度拐點,確定低溫半致死溫度(LT50);低溫半致死溫度越低,抗寒性越強[31],‘雙優’和‘北冰紅’低溫半致死溫度分別為-38.63和-29.75 ℃(表1),可見 2個品種抗寒性差異明顯;圖像可視化評估所求得的染色指數也能將‘雙優’和‘北冰紅’的抗寒性有效分開(圖2a);同時不同低溫處理下枝條的水培萌芽情況與2個品種抗寒能力相符(圖2b;表2)。說明TTC染色圖像可視化評估配合 Logistic方程鑒定山葡萄種質抗寒性具有一定的可行性,40份山葡萄種質抗寒性鑒定結果詳見表1。

表1 40份山葡萄種質抗寒性鑒定結果Table 1 Identification of cold resistance in 40 V. amurensis germplasms

表2 低溫處理下2個山葡萄品種枝條水培萌枝率Table 2 Shoots sprout rate of in two varieties of V. amurensis in hydroponic culture under cold temperatures %

通過對表1中同一溫度下40份山葡萄種質的染色指數與其低溫半致死溫度的相關性分析表明,-20、-25、-30、-35和-40 ℃處理下的枝段染色指數與其低溫半致死溫度之間均呈極顯著(P<0.000 1)負相關,其中-20 ℃下,相關系數r=-0.632 55,P<0.000 1;-25 ℃下,r=-0.791 16,P<0.000 1;-30 ℃下,r=-0.832 20,P<0.000 1;-35 ℃下,r=-0.861 38,P<0.000 1;-40 ℃下,r=-0.726 03,P<0.000 1。可見兩者在-35 ℃下的相關性最強,因此本研究將-35 ℃作為最適冷凍溫度,進一步利用隸屬函數法對最適冷凍溫度下40山葡萄種質枝條的抗寒性進行綜合評價。

圖2 染色指數及水培萌枝率在冷凍溫度下的變化Fig.2 Change of staining index and hydroponic culture sprout rate under cold temperatures

2.2 隸屬函數法綜合評價山葡萄種質的抗寒性

枝條經-35 ℃最適冷凍溫度處理后,對40份山葡萄種質抗寒指標進行相關性分析結果表明,指標變量間存在一定的相關性(表 3),說明多個指標提供的抗寒信息有交叉重疊現象,進一步利用主成分分析對所測指標進行綜合評價。

表3 最適冷凍溫度下山葡萄種質抗寒性指標的相關性分析Table 3 Correlation analysis of cold resistance indexes in V.amurensis germplasms under optimum freezing temperature

以累計貢獻率大于85%為原則選擇主成分,其中第1主成分中特征向量絕對值較大的是相對電導率、丙二醛及可溶性蛋白含量,貢獻率達41.48%;第2主成分則是萌枝率和凍害指數,貢獻率達25.70%;第3主成分中可溶性糖含量特征向量絕對值較大,貢獻率為15.53%;第4主成分的游離脯氨酸含量特征向量絕對值較大,貢獻率為10.90%。4個獨立主成分的累計貢獻率為93.61%(表4),利用隸屬函數法對以上指標進行綜合分析,以評價山葡萄種質的綜合抗寒性。

表4 抗寒指標主成分分析Table 4 Principal components analysis for cold resistance indexes

依據平均隸屬度值將40份山葡萄種質抗寒性分為5個級別,從Ⅴ級到Ⅰ級抗寒性依次降低。其中Ⅴ級抗寒種質1份;Ⅳ級抗寒種質10份;Ⅲ級抗寒種質24份;Ⅱ級抗寒種質3份;Ⅰ級抗寒2份為山歐雜種(表5)。結果表明,40份山葡萄種質抗寒性存在明顯差異。整體來看,純種山葡萄抗寒性 > 山歐雜種(對照)。

表5 隸屬函數法綜合評價山葡萄種質抗寒性Table 5 Comprehensive analysis of cold resistances in V. amurensis germplasms

2.3 田間萌芽率

本試驗連續3 a對40份山葡萄種質自然越冬后枝條的萌芽情況進行了調查,發現同一資源各年份間的田間萌芽率無明顯差異。但由于田間抗寒性鑒定數據的年份要與實驗室抗寒性鑒定數據的年份相對應,所以本試驗田間萌芽率為2017年的調查結果(表6)。圖3為2017年1—5月的氣象數據,最低氣溫出現在1月下旬,最低為-28.2 ℃。經自然越冬后,40份山葡萄種質因抗寒性存在差異,其田間萌芽率也表現出不同(表6)。

2.4 相關性分析

利用SAS相關性分析結果表明,TTC染色圖像可視化評估量化獲得的染色指數配合Logistic方程求得的40份山葡萄種質低溫半致死溫度與其抗寒平均隸屬度、田間萌芽率均呈極顯著負相關(P<0.000 1),相關系數分別為-0.883 03(圖4a)和-0.876 70(圖4b),說明該方法與隸屬函數法對山葡萄種質抗寒性鑒定結果較為一致;且染色指數擬合計算出的低溫半致死溫度與實際田間萌芽率的線性關系為 y=-5.576 5x–127.36,為驗證模型準確性,另從源圃中隨機選取了15份山葡萄種質,利用該線性關系方程計算出了 15份種質田間萌芽率的理論值(表7),與田間萌芽率的實際值相關性為-0.985 84,說明TTC染色圖像可視化評估配合Logistic方程求得的山葡萄種質低溫半致死溫度能夠反映其在田間的實際抗寒性,且建模模型準確。

表6 山葡萄種質田間萌芽率調查結果Table 6 Investigation on germination rate of V. amurensis germplasms in field

圖3 2017年氣象數據Fig.3 Meteorological data in 2017

圖4 抗寒性鑒定結果的相關性分析Fig.4 Correlation analysis of results in cold resistances identification

表7 抗寒性數據模型驗證Table 7 Proof of cold resistance data model

3 結 論

1)TTC染色枝段圖像的可視化評估實現了染色面積的準確量化,量化得出的染色指數與冷凍溫度呈典型的倒‘S’型曲線,說明可以利用非線性回歸分析結合Logistic 方程擬合溫度拐點,確定低溫半致死溫度。

2)該方法量化得出的染色指數配合Logistic方程確定的40份山葡萄種質低溫半致死溫度與其抗寒性綜合評價的平均隸屬度值及田間萌芽率值均呈極顯著負相關(P<0.000 1),相關系數分別為-0.883 03和-0.876 70。

3)染色指數擬合計算出的低溫半致死溫度與實際田間萌芽率的線性關系為y= - 5.576 5 x-127.36。利用該線性關系方程計算了建模數據以外隨機選取的山葡萄種質田間萌芽率的理論值,與田間萌芽率的實際值相關性為-0.985 84,說明TTC染色圖像可視化評估配合Logistic方程求得的山葡萄種質低溫半致死溫度能夠反映其在田間的實際抗寒性。

4)40份山葡萄種質低溫半致死溫度在-39.41~-29.75 ℃范圍內,種質抗寒性存在明顯差異;總體而言,該研究的純種山葡萄低溫半致死溫度的平均值為-37.25 ℃,山歐雜種山葡萄低溫半致死溫度的平均值為-29.81 ℃,可見純種山葡萄的抗寒性強于山歐雜種山葡萄。

[1] 賀普超, 晁無疾. 我國葡萄屬野生種質資源的抗寒性分析[J]. 園藝學報,1982,9(3):17-21.He Puchao,Chao Wuji. Studies on the cold hardiness of eight species of vitis L.growing wild in China[J]. Acta Horticulturae Sinica, 1982, 9(3): 17-21. (in Chinese with English abstract)

[2] 沈育杰,史貴文,徐浩,等. 我國山葡萄種質資源研究與利用[J]. 特產研究,1992(3):29-32.Shen Yujie, Shi Guiwen, Xu Hao, et al. Research and utilization of Vitis amurensis germplasm resources in China[J]. Special Wild Economic Animal and Plant Research,1992(3): 29-32. (in Chinese with English abstract)

[3] Welling A, Moritz T, Palvae T, et al. Independent activation of cold acclimation by low temperature and short photoperiod in hybrid aspen[J]. Plant Physiology, 2002, 129(4): 1633-1641.

[4] Junttila O, Kaurin L. Environmental control of cold acclimation in Salix pentandra[J]. Scandinavian Journal of Forest Research, 2008, 5(1): 195-204.

[5] 張劍俠,吳行昶,楊亞州. 葡萄種間雜種的抗寒性評價[J].北方園藝,2014(13):1-6.Zhang Jianxia, Wu Xingchang, Yang Yazhou. Cold-resistance evaluation in grape interspecific hybrids[J]. Northern Horticulture, 2014(13): 1-6. (in Chinese with English abstract)

[6] Deans J D, Billington F J, Harvey F J. Assessment of frost damage to leafless stem tissues of Quercus petraea: A reappraisal of the method of relative conductivity[J]. Forestry,1995, 68: 25-34.

[7] 張倩,劉崇懷,郭大龍,等. 5個葡萄種群的低溫半致死溫度與其抗寒適應性的關系[J]. 西北農林科技大學學報(自然科學版),2013,41(5):149-154.Zhang Qian, Liu Chonghuai, Guo Dalong, et al. Relationship between LT50 and cold adaptability of five grape varieties[J].Journal of Northwest A & F University (Nat. Sci. Ed.), 2013,41(5): 149-154. (in Chinese with English abstract)

[8] Zhang M I N, Willison J H M. An improved conductivity method for the measurement of frost hardiness[J]. Canadian Journal of Botany, 1987, 65: 710-715.

[9] Stergios B G, Howell G S. Evaluation of viability tests for cold stressed plants[J]. Journal of the American Society for Horticultural Science, 1973, 98(4): 325-330.

[10] 安徽農業大學. 一種茶樹抗寒性早期鑒定方法:201110045014.1[P]. 2012-08-29.

[11] 張文娥,王飛,潘學軍. 應用隸屬函數法綜合評價葡萄種間抗寒性[J]. 果樹學報,2007,24(6):849-853.Zhang Wene, Wang Fei, Pan Yiingjun. Comprehensive evaluation on cold hardiness of Vitis species by subordinate function[J]. Journal of Fruit Science, 2007, 24(6): 849-853.(in Chinese with English abstract)

[12] 劉杜玲,張博勇,孫紅梅,等. 早實核桃不同品種抗寒性綜合評價[J]. 園藝學報,2015,42(3):545-553.Liu Duling, Zhang Boyong, Sun Hongmei, et al. Comprehensive evaluation on drought resistance of early fruiting walnut cultivars[J]. Acta Horticulturae Sinica, 2015, 42(3): 545-553. (in Chinese with English abstract)

[13] 令凡,焦健,李朝周,等. 不同油橄欖品種對低溫脅迫的生理響應及抗寒性綜合評價[J]. 西北植物學報,2015,35(3):508-515.Ling Fan, Jiao Jian, Li Chaozhou, et al. Physiological response and comprehensive evaluation of cold resistance under cold stress for different varieties of Olea europaea[J].Acta. Bot. Boreal-Occident Sin, 2015, 35(3): 508-515. (in Chinese with English abstract)

[14] 劉海卿,孫萬倉,劉自剛,等. 北方寒旱區白菜型冬油菜抗寒性與抗旱性評價及其關系[J]. 中國農業科學,2015,48(18):3743-3756.Liu Haiqing, Sun Wancang, Liu Zigang, et al. Evaluation of drought resistance and cold resistance and research of their relationship at seedling stage of winter rapeseed (Brassica campestris L.) in cold and arid regions in North China[J].Scientia Agricultura Sinica, 2015, 48(18): 3743-3756. (in Chinese with English abstract)

[15] 張海晶,張金梅,辛霞,等. 馬鈴薯莖尖超低溫保存流程TTC活力響應[J]. 植物遺傳資源學報,2015,16(3):555-560.Zhang Haijing, Zhang Jinmei, Xin Xia, et al. Study on potato shoot tips vitality responses during cryopreservation by TTC staining method[J]. Journal of Plant Genetic Resources, 2015,16 (3) :555-560. (in Chinese with English abstract)

[16] 張啟翔. 梅花及其雜交種根系抗寒性研究初報[J]. 北京林業大學學報,1992,14(增刊4):83-86.Zhang Qixiang. A preliminary study on cold resistance of root system of Mei and its hybrid cultivars[J]. Journal of Beijing Forestry University, 1992, 14 (Supp.4): 83-86. (in Chinese with English abstract)

[17] Thomas F M, Ahlers U. Effects of excess nitrogen on frost hardiness and freezing injury of above-ground tissue in you6ng oaks (Quercus petraea and robur)[J]. New phytologist,2010, 144(1): 73-83

[18] 董麗,黃亦工,賈麥娥,等. 北京園林主要常綠闊葉植物抗凍性及其測定方法[J]. 北京林業大學學報,2002,24(3):70-73.Dong Li, Huang Yigong, Jia Maie, et al. Freezing tolerance and methods of evaluation of introduced evergreen broad-leaf plants in open gardens in Beijing of China[J]. Journal of Beijing Forestry University, 2002, 24 (3): 70-73. (in Chinese with English abstract)

[19] 田景花,王紅霞,高儀,等. 核桃屬植物休眠期的抗寒性鑒定[J]. 園藝學報,2013,40(6):1051-1060.Tian Jinghua, Wang Hongxia, Gao Yi, et al. Assessment of freezing tolerance of Juglans Germplasms by using annual dormant branches[J]. Acta Horticulturae Sinica, 2013, 40 (6):1051-1060. (in Chinese with English abstract)

[20] 艾軍,沈育杰. 山葡萄規范化栽培與釀酒技術[M]. 北京:中國農業出版社,2017.

[21] Zhang J,Wu X,Niu R,et al. Cold-resistance evaluation in 25 wild grape species[J]. Vitis, 2012, 51(4):153-160.

[22] 解越,張敏,梁飛俠,等. 電導法配合Logistic方程鑒定西葫蘆和黃瓜果實抗寒性[J]. 農業工程學報,2015,31(7):292-295.Xie Yue,Zhang Min,Liang Feixia, et al. Identification of cold resistance of squash and cucumber fruits by electrical conductivity and Logistic equation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(7): 292-295. (in Chinese with English abstract)

[23] Verslues P E, Agarwal M, Katiyar-Agarwal S, et al. Methods and concepts in quantifying resistance to drought, salt and freezing, abiotic stress that affect plant water status[J]. Plant Journal, 2006, 45(4): 523-539.

[24] Yang S H, Wang J, Li S H, et al. The effects of UV-B radiation on photosynthesis in relation to Photosystem II photochemistry, thermal dissipation and antioxidant defenses in winter wheat (Triticum aestivum L.) seedlings at different growth temperatures[J]. Functional Plant Biology, 2007,34(10): 907-917.

[25] 李合生. 植物生理生化實驗原理和技術[M]. 北京:高等教育出版社,2000.

[26] 高俊鳳. 植物生理學實驗指導[M]. 北京:高等教育出版社,2006.

[27] Bradford M M. A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding[J]. Analytical Biochemistry,1976, 72, 248-254.

[28] Nanami S, Kawaguchi H, Tateno R, et al. Sprouting traits and population structure of co-occurring Castanopsis species in an evergreen broad- leaved forest in southern China[J].Ecological Research, 2004, 19(3): 341-348.

[29] 呂中偉,王鵬,許領軍,等. 河南省葡萄凍害調查及應對措施[J]. 山西果樹,2011(1):37-38.Lu Zhongwei, Wang Peng, Xu Lingjun, et al. Investigation and countermeasures of grape freeze injury in Henan[J].Shanxi Fruits, 2011(1): 37-38. (in Chinese with English abstract)

[30] 劉崇懷,沈育杰,陳俊. 葡萄種質資源描述規范和數據標準[M]. 北京: 中國農業出版社,2006.

[31] 于立洋,李政,韓佩堯,等. 8個新疆野蘋果優良無性系抗寒性比較[J]. 核農學報,2017,31(9):1827-1835.Yu Liyang, Li Zheng, Han Peiyao, et al. Comparison of cold hardiness of 8 Malus sieversii Clones[J]. Journal of Nuclear Agricultural Sciences, 2017, 31(9):1827-1835. (in Chinese with English abstract)

主站蜘蛛池模板: 久久香蕉国产线看观| 精品视频一区二区三区在线播| 黄色网站不卡无码| 亚洲精品无码不卡在线播放| 2022国产91精品久久久久久| 91小视频在线观看免费版高清| 亚洲Av综合日韩精品久久久| 欧美精品1区| 国产香蕉在线| 麻豆精品在线播放| 国产精品亚欧美一区二区| 爱色欧美亚洲综合图区| 欧美区一区二区三| 97影院午夜在线观看视频| 亚洲免费成人网| 亚洲aⅴ天堂| 美女无遮挡拍拍拍免费视频| 沈阳少妇高潮在线| 日韩AV无码免费一二三区| 亚洲天堂免费在线视频| 成年片色大黄全免费网站久久| 久久综合久久鬼| 免费看一级毛片波多结衣| 精品一区二区久久久久网站| 国产小视频免费| 日韩欧美中文字幕在线韩免费 | 久久久久国产精品熟女影院| 亚洲欧美日韩高清综合678| 国产成人综合在线视频| 色吊丝av中文字幕| 亚洲欧美人成人让影院| 亚洲欧美另类久久久精品播放的| 亚洲美女高潮久久久久久久| 精品国产成人高清在线| 国产精品熟女亚洲AV麻豆| 日韩大片免费观看视频播放| 国产成人艳妇AA视频在线| 成人午夜网址| 欧美午夜视频在线| 成年女人a毛片免费视频| 久久久噜噜噜久久中文字幕色伊伊 | 国产一级无码不卡视频| 国产va在线观看| 国产无码高清视频不卡| 亚洲成av人无码综合在线观看| 亚洲日本一本dvd高清| 国产中文一区a级毛片视频 | 久久午夜影院| 国产精品视频久| 最新日本中文字幕| 亚洲无码免费黄色网址| 日韩经典精品无码一区二区| 91原创视频在线| 日本免费一级视频| а∨天堂一区中文字幕| 国产免费羞羞视频| 久久久久久尹人网香蕉| 伊人久久大香线蕉aⅴ色| 精品无码专区亚洲| 国产在线91在线电影| 无码精油按摩潮喷在线播放 | 五月丁香伊人啪啪手机免费观看| 亚洲国产91人成在线| 亚洲美女一区二区三区| 在线观看无码a∨| 久久久91人妻无码精品蜜桃HD| 亚洲精品视频免费看| 日本a级免费| 欧美精品不卡| 国产va视频| 国产真实乱人视频| 久久中文电影| 国产情侣一区二区三区| 国产日韩精品欧美一区灰| 国产高清在线观看| 青青青视频蜜桃一区二区| 国产精品林美惠子在线播放| 强奷白丝美女在线观看| 99免费视频观看| 在线国产91| 国产av无码日韩av无码网站| 韩国v欧美v亚洲v日本v|