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

    基于田間性狀和花部性狀的茶用菊品系高產(chǎn)及適應性評價

    2022-01-27 11:04:56李嘉偉李曉宇史亞東陳素梅房偉民張飛趙爽陳發(fā)棣管志勇
    南京農(nóng)業(yè)大學學報 2022年1期
    關鍵詞:品系適應性單株

    李嘉偉,李曉宇,史亞東,陳素梅,房偉民,張飛,趙爽,陳發(fā)棣,管志勇*

    (1.南京農(nóng)業(yè)大學作物遺傳與種質(zhì)創(chuàng)新國家重點實驗室/農(nóng)業(yè)農(nóng)村部景觀設計重點實驗室/園藝學院,江蘇 南京 210095;2.安徽菊泰滁菊草本科技有限公司,安徽 滁州 239000)

    菊花(Chrysanthemummorifolium)是中國的傳統(tǒng)名花,其在中國已有大約三千年的栽培歷史,在人工栽培之前其食用、藥用價值已經(jīng)就被明確提到[1]。目前茶用和藥用兼用的功能性菊花產(chǎn)業(yè)規(guī)模不低于觀賞菊,僅江蘇射陽、湖北麻城、安徽黃山周邊地區(qū)、浙江桐鄉(xiāng)等地就接近20萬hm2。菊花茶的色味、香味和形態(tài)俱佳,具有清心解暑、除煩明目、潤喉生津的功能[2],長期飲用可以增加人體鈣質(zhì)、調(diào)節(jié)心肌功能、降低膽固醇[3]。當前我國茶用菊品種主要是‘杭白菊’‘滁菊’‘貢菊’‘福白菊’等各地方品種及‘金絲黃菊’等。與切花菊的設施生產(chǎn)方式不同,茶用菊的生產(chǎn)為露天栽培,植株在定植后的整個生育期都處于露天環(huán)境,期間不利的環(huán)境因素如高溫、寒冷、干旱、雨澇以及多種病蟲害的逆境侵襲無法避免,這些因素都會導致茶用菊生長發(fā)育異常,甚至減產(chǎn)。目前的主栽品種均不同程度存在對田間多種環(huán)境脅迫適應性差的缺點,各產(chǎn)區(qū)的減產(chǎn)、絕收現(xiàn)象時有發(fā)生,這限制了菊農(nóng)的生產(chǎn)積極性,導致產(chǎn)業(yè)萎縮。選擇對田間自然環(huán)境適應性強的品種進行生產(chǎn)是茶用菊穩(wěn)產(chǎn)高產(chǎn)的基礎。因此迫切需要對現(xiàn)有品種開展生長特性以及田間適應性等方面進行評價,篩選綜合性狀優(yōu)良的茶用菊品種,推動品種創(chuàng)新。因此,構(gòu)建科學、系統(tǒng)的茶用菊綜合評價體系,開展茶用菊品種評價及篩選,可為茶用菊生產(chǎn)、品種選育提供實踐支撐。

    近年來,對觀賞植物品種評價逐步由定性、單一指標向定量、多指標的綜合評價過渡[4]。目前常用于栽培植物品種評價的方法有百分制記分法[5]、層次分析法[6]、灰色關聯(lián)法[7]、模糊數(shù)學法[8]和主成分分析法[9]。其中,層次分析法可以滿足植物品種綜合評價的需要,其突出特點是將定性與定量結(jié)合起來,將主觀判斷和定性分析用數(shù)量分析表述、轉(zhuǎn)換和處理[10],在栽培植物品種評價方面應用較多。目前已在園林小菊[11]、多頭切花菊[12]、單頭切花菊[13]、切花小菊[14]、盆栽多頭小菊[15]以及中國傳統(tǒng)菊[16]品種資源評價上使用,篩選出一系列優(yōu)良品種,而茶用菊品種評價方面的工作則側(cè)重于干花加工后的外觀品質(zhì)、茶湯品質(zhì)和功能性成分的評價[17],但茶用菊的豐產(chǎn)性及適應性則是產(chǎn)業(yè)鏈種植環(huán)節(jié)考慮的首要因素。近年來,低產(chǎn)的品種、適栽性差的品種通常被產(chǎn)區(qū)主動淘汰,啟示茶用菊育種者需給予茶用菊的豐產(chǎn)性及適應性足夠重視。因此本研究基于40個茶用菊品系的田間和花部性狀,采用專家咨詢法、層次分析法以及K-means 聚類分析法構(gòu)建茶用菊的綜合評價體系。

    1 材料與方法

    1.1 試驗材料和試驗地點

    供試材料為40個茶用菊品系,包括‘滁菊’‘杭白菊’‘異種大白花’3個品種,其余37個株系為近年來從茶用菊品種開放授粉的有性后代篩選出的優(yōu)良品系。37個品系的開放授粉育種親本為‘滁菊’‘貢菊’‘杭白菊’。

    試驗在安徽菊泰滁菊草本科技有限公司的生產(chǎn)區(qū)進行。于2018年4月下旬至5月中旬挑選生長強健的單株,取長6~8 cm、長勢整齊、充實健康的穗條,先用多菌靈消毒,再將穗條在生根粉水溶液中蘸取30 s左右后扦插到穴盤中?;|(zhì)采用草炭和石灰?guī)r(體積比3∶1)混合配制。扦插12 d生根完畢后以種植間距為50 cm定植于露地,每個品種或株系定植150株。定植過后20 d進行第1次摘心,再過15 d第2次摘心,以促進側(cè)枝的生長。40個茶用菊品種(系)分2個時間批次即4月20號和5月1號扦插、定植,以獲得2個試驗重復。在各品種進入盛花期時對其性狀觀察測定。

    安徽滁州地區(qū)為北亞熱帶濕潤季風氣候,年平均氣溫15.4 ℃,年平均降水量1 035.5 mm,年日照總時數(shù)2 073.4 h,年無霜期210 d。

    1.2 性狀測定與方法

    1.2.1 性狀篩選基于菊花栽培鑒定試驗或室內(nèi)分析測試(DUS測試)指南[18]完成了對40個茶用菊品系24個性狀的調(diào)查。最終參考茶用菊產(chǎn)業(yè)專家、生產(chǎn)企業(yè)和銷售商的意見,篩選出10個影響茶用菊產(chǎn)業(yè)化生產(chǎn)的關鍵指標,分別為花朵在枝條上的分布(flower distribution,FD)、植株緊湊性(plant compactedness,PC)、株高(plant height,PH)、花序高(inflorescence height,IH)、花徑(inflorescence diameter,ID)、舌狀花數(shù)(number of ligulate flowers,NLF)、花朵重瓣性(polyphyll of flower,PF)、田間適應性(field adaptability,FA)、收獲期的直立性(upright of harvest time,UH)、單株花數(shù)(flower number,FN)。

    1.2.2 性狀測定方法花朵在枝條上的分布:目測確定各品種的花朵在枝條上的位置,花朵集中分布在上部為1級,多數(shù)分布在上部為2級,在枝條的上下都有分布為3級?;ǘ湓谥l上的分布級別見圖1。

    圖1 花朵在枝條上的分布級別圖Fig.1 The grade diagram of the distribution of flowers on the branches A、B和C分別表示在枝頭集中分布、多數(shù)分布于上部以及在枝條上下都有分布,依次賦值3、2和1。A,B and C represent concentrated distribution in the branches,mostly in the upper part,above and below the branches,with a score of 3,2 and 1,respectively.

    植株緊湊性:目測確定各品種的植株緊湊程度,以花枝間緊湊程度確定標準。植株緊湊性級別見圖2。

    圖2 植株緊湊性級別圖Fig.2 The grade diagram of plant compactednessA、B和C分別表示緊湊、中等和松散,依次賦值3、2和1。A,B and C represent compact,medium and loose,with a score of 3,2 and 1,respectively.

    田間適應性:盛花期目測各品種群體缺株數(shù)占總定植數(shù)的比例,以及群體長勢。田間適應性級別見圖3。

    圖3 田間適應性級別圖Fig.3 The grade diagram of field adaptabilityA、B和C分別表示強、中等和弱適應性,依次賦值3、2和1。A,B and C represent strong,medium and weak adaptability,with a score of 3,2 and 1,respectively.

    收獲期花枝的直立性:目測各品種群體發(fā)生花枝倒伏的單株數(shù)量占總定植數(shù)的比例。

    株高、花序高、花徑、舌狀花數(shù)、花重瓣性、單株花數(shù)的測試方法見菊花DUS測試指南[18]?;ㄖ匕晷约墑e見圖4。

    圖4 花重瓣性級別圖Fig.4 The grade diagram of polyphyll of the flower A、B和C分別表示緊湊、中等和松散且重瓣性好、中等和差,依次賦值3、2和1。A,B and C represent compact,medium and loose,and polyphyll of the flower is good,medium and poor,with a score of 3,2 and 1,respectively.

    株高:菊花植株頂部到基部間的距離。

    花序(菊花俗稱的一朵花實質(zhì)上是一個頭狀花序)高:花序頂部到基部間的距離。

    單株花數(shù):破壞性計數(shù),即從植株上摘取花朵并計數(shù)。

    舌狀花數(shù):破壞性計數(shù),即從花序上拔取舌狀花,并計數(shù)。

    數(shù)量性狀通過從株系群體中選取表型性狀較為平均的單株和花朵樣本進行測定,測量10個樣本后取平均值。

    1.3 統(tǒng)計與評價

    用Excel 2019對各性狀的實測數(shù)據(jù)進行整理和統(tǒng)計;用yaahp軟件進行層次分析;用SPSS 20軟件進行K-means 聚類分析[19]。

    1.3.1 評價體系構(gòu)建參考前人的研究方法[11,20],基于茶用菊品種特點,根據(jù)評價指標的建立原則[21]以及體系的框架模型,結(jié)合專家意見征詢確定評價指標以及指標間的隸屬關系,采用1~9標度法(表1)比較指標間的相對重要性,量化各項指標并構(gòu)造判斷矩陣,構(gòu)建出茶用菊層次結(jié)構(gòu)模型(表2),包含目標層(A)、準則層(B)和方案層(C)。

    表1 1~9比率尺度值含義Table 1 The meanings of ratios of 1-9

    表2 各評價體系的層次結(jié)構(gòu)模型Table 2 The hierarchy model of evaluation system

    1.3.2 判斷矩陣及一致性檢驗在層次分析法(analytic hierarchy process,AHP)綜合評價體系中,各評價指標的重要性是相對的,對各個評價指標進行兩兩比較可以衡量指標的重要程度。根據(jù)上面的1~9標度法以及得到的分層結(jié)構(gòu)模型,結(jié)合茶用菊的栽培和生產(chǎn)特性以及相關專家的意見,構(gòu)建出A-B、B-C共 4個判斷矩陣,并進行一致性檢驗(表3)。

    表3 各判斷矩陣及其一致性檢驗Table 3 The results of matrix and its consistency test

    運用yaahp 軟件計算判斷矩陣的最大特征根、相應的特征向量以及權(quán)重值,判斷矩陣計算方法采用“冪法”。一致性指標(consistency index,CI)為度量判斷矩陣偏離一致性的指標,隨機一致性指標(random index,RI)為判斷矩陣的平均隨機一致性指標。以判斷矩陣的隨機一致性比率(consistency ratio,CR)為標準,即判斷矩陣的一般CI與RI的比值。CR<0.10,則說明判斷矩陣一致性滿意,各指標的權(quán)重賦值均符合要求。

    2 結(jié)果與分析

    2.1 茶用菊評分標準的制定

    以挑選可作產(chǎn)業(yè)規(guī)?;a(chǎn)的高產(chǎn)優(yōu)質(zhì)茶用菊為基準,根據(jù)40個茶用菊品種在花期的實測數(shù)據(jù)結(jié)合生產(chǎn)特點、專家意見以及市場需求,將影響茶用菊產(chǎn)量和品質(zhì)的10個指標劃分為3個級別,賦值3、2、1,擬定了茶用菊品種評選的3分制評分標準(表4)。

    表4 各因子評分標準Table 4 Scoring standard of each factor

    2.2 各評價因子的權(quán)重值

    根據(jù)表3的判斷矩陣及一致性檢驗,計算得出方案層C對準則層B以及準則層B對總目標層A的權(quán)重值,再進行加權(quán)綜合,可以得出方案層C對總目標層A的權(quán)重值。由表5可見:單株花數(shù)對茶用菊的影響最大,其權(quán)重值為0.301,其次是田間適應性和舌狀花數(shù),其權(quán)重值分別為0.172和0.119;而花重瓣性(0.113)、花朵在枝條上的分布(0.088)、收獲期花枝的直立性(0.066)、植株緊湊性(0.048)的權(quán)重值對茶用菊的影響位于中等水平;花徑(0.042)、株高(0.027)、花序高(0.024)的權(quán)重值對茶用菊的影響相對較小,即通過計算方案層的權(quán)重值可以得出各個評價指標的相對重要性,從而能夠比較準確、量化并且全面地對茶用菊進行綜合評價。

    表5 各評價因子的權(quán)重值Table 5 The weight value of each factor

    2.3 茶用菊品系的得分及等級

    將茶用菊各指標的評分標準所得的分值與各影響因子C相對于A層的權(quán)重值相乘,計算出每一品種的綜合得分。應用K-means聚類分析法將40個茶用菊品系分成優(yōu)、良、中、差4個級別,依次用數(shù)字1、2、3、4來表示,其聚類中心分別為2.72、2.33、2.01、1.72。40個茶用菊品種的綜合得分以及等級劃分如表6所示。

    表6 茶用菊品系綜合得分及等級Table 6 The score and grade of chrysanthemum for tea

    2.4 評價結(jié)果分析

    茶用菊各性狀的權(quán)重值顯示(表5),適應性的權(quán)重值最高,以單株花數(shù)最為重要,其次是田間適應性和收獲期花枝的直立性;生長特性指標中花朵在枝條上的分布最為重要,其次是植株緊湊性和株高;花部特征指標中舌狀花數(shù)最為重要,其次是花重瓣性和花徑,而花序高重要性程度最低。

    從方案層C對總目標層A的權(quán)重值結(jié)果來看,10個評價指標的權(quán)重值中單株花數(shù)和田間適應性的權(quán)重最高,這2個指標直接反映茶用菊的產(chǎn)量潛力和田間穩(wěn)產(chǎn)性能,故所占權(quán)重值高?;ㄐ蚴遣栌镁丈a(chǎn)的收獲產(chǎn)品,舌狀花數(shù)和花重瓣性與茶用菊的產(chǎn)量及外觀品質(zhì)直接相關,也在品種評價中具有重要位置。花朵在枝條上的分布、收獲期花枝的直立性以及植株緊湊性主要是影響茶用菊的易采收性和花朵的清潔。株高、花徑和花序高3個性狀指標則與產(chǎn)量及品質(zhì)的直接關聯(lián)較小,所占權(quán)重較小。由此可見,對茶用菊產(chǎn)量和品質(zhì)有較大影響的指標權(quán)重值大,體現(xiàn)了這一評價體系的合理性,符合本次評價的目的。

    由綜合評價結(jié)果(表6)可知,40個茶用菊品系按優(yōu)、良、中、差4個等級劃分的比例為12.5%、35.0%、32.5%、20.0%。其中,X1-56、X1-33、X1-25、X1-51和X1-29共5個株系處于1級優(yōu)秀;X1-19、X1-05、‘滁菊’、X1-48等14個品種(系)處于2級良好;X1-01、X1-30、X1-40、X1-57等13個株系處于3級一般;X1-32、X1-03、X1-36、X1-06等8個株系的綜合品質(zhì)較差,處于4級差等株系。其中1級和2級品系在單株花數(shù)、田間適應性和舌狀花數(shù)等權(quán)重較大的指標上表現(xiàn)突出,綜合評價表現(xiàn)優(yōu)良,符合篩選需求,占全部品系的47.5%,對這些品種可以加大推廣和種植力度;3級品系可利用1、2等級的優(yōu)良品種進行品系優(yōu)化;而4級品系綜合評價表現(xiàn)較差,不適合在生產(chǎn)中推廣。

    3 討論與結(jié)論

    層次分析法利用構(gòu)造兩兩對比矩陣的方式來確定各評價指標對綜合性狀的權(quán)重,將主觀判斷的定性問題定量化,使決策過程有很強的客觀性以及科學性,將定性標度與定量標度相結(jié)合有利于客觀全面的評價,而定性信息標度的可信度決定了判斷矩陣質(zhì)量[22]。目前已有不少研究運用層次分析法對不同類型的菊花品種資源進行了評價,篩選出一系列綜合性狀優(yōu)良的品種。朱德寧等[11]運用聚類分析結(jié)合層次分析的方法對多花型園林小菊建立了評價體系,篩選出‘鐘山金陽’等12個優(yōu)等多花型園林小菊品種。張亞瓊等[16]通過測定株高、莖粗、花莖等相關性狀,從觀賞性狀、抗逆性、栽培流程 3 個方面確定了 19 個評價因子的權(quán)重值,使用層次分析法篩選出適合產(chǎn)業(yè)化生產(chǎn)的中國傳統(tǒng)盆栽菊,構(gòu)建了傳統(tǒng)菊花品種的綜合評價體系。對于茶用菊品種而言,單株花數(shù)、田間適應性等與產(chǎn)量相關的指標及舌狀花數(shù)、花重瓣性等與外觀品質(zhì)相關指標,還有涉及生產(chǎn)采收環(huán)節(jié)的花朵在枝條上的分布等指標都是相關的評價指標。因此,在茶用菊品種的綜合評價中,綜合各個指標的信息并科學地確定各個指標的權(quán)重是茶用菊綜合評價的重要內(nèi)容。本研究中,田間適應性、花朵在枝條上的分布、植株緊湊性、收獲期花枝的直立性這4個定性指標與茶用菊品質(zhì)優(yōu)劣判斷的關系明確,定性標度質(zhì)量高;6個定量指標與茶用菊優(yōu)劣關系更加明確,采用層次分析法評價優(yōu)勢明顯。孫明等[23]研究 17 個地被菊品系的株形、開花、適應性等14個性狀,通過性狀分級賦值,利用層次分析法篩選出 5 個觀賞價值較高的地被菊優(yōu)良品系。馬婉茹等[20]采用層次分析法對300個多頭切花菊品種的莖、枝特性建立評價體系,篩選出176個優(yōu)良品種并在生產(chǎn)上推廣應用。與注重個體外觀品質(zhì)的地被菊和切花菊不同的是,茶用菊在主產(chǎn)區(qū)是規(guī)模化露天生產(chǎn),各種非生物逆境及病害、蟲害等生物脅迫,都會影響茶用菊的正常生長,田間甚至發(fā)生菊花對干旱、濕澇、病害綜合因素的不適應而發(fā)生的大范圍缺株、全田絕收的情形,嚴重影響產(chǎn)業(yè)的發(fā)展。因此,使用對田間自然環(huán)境適應性強的品種進行生產(chǎn)是茶用菊穩(wěn)產(chǎn)高產(chǎn)的基礎,通過品種改良從根本上來解決適應性差的問題。群體的田間適應性與長勢、單株花數(shù)是產(chǎn)量的決定性因素,此外,花朵在枝條上的位置、植株緊湊性、收獲期直立性影響了采收的效率和花朵的損耗多少,將這些指標進行量化可以直觀地反映茶用菊品種的田間生產(chǎn)狀況。

    翟麗麗等[24]分別采用盆栽淹水法、自然干旱處理比較國慶小菊的耐澇、旱性;張常青等[25]采用自然失水脅迫和PEG脅迫處理方法,研究地被菊花幼苗的耐寒性;尹冬梅等[26]采用土培模擬澇害的方法比較菊花及其近緣種屬植物耐澇性;孔令接等[27]在苗期和花期進行高溫脅迫處理,比較不同夏菊品種的耐熱性;徐婷婷等[28]利用低溫半致死溫度評價切花菊的耐寒性;李媛媛等[29]利用苗期人工接種法評價了茶用菊的抗病性。上述研究均只對單個生態(tài)因子進行評價,無法獲得涵蓋多種生態(tài)因子在內(nèi)的田間適應性的結(jié)論。在茶用菊的田間生長期間,上述逆境因子都不同程度存在,在這些逆境因子的綜合影響下,茶用菊不同品系的田間成活狀況、田間長勢均會有所差異。充分比較不同品系在大田條件下生長差異就能獲得其對各種生態(tài)因子的適應性差異的結(jié)論,這種方法在花卉[30-31]等園藝作物的種質(zhì)資源篩選上普遍采用。本研究的田間篩選地點位于華東地區(qū)滁州市的茶用菊傳統(tǒng)產(chǎn)區(qū),在露天栽培條件下,篩選群體能充分接受梅雨、濕熱、高溫、干旱、連作病原等因子的影響,具有大田篩選的理想自然條件。以該條件下生長的茶用菊各品系群體最終的成活株數(shù)、生長勢作為生態(tài)適應性共同指標,更具有實踐指導意義。

    史峰霖[32]運用專家咨詢法、層次分析法和K-means聚類分析法構(gòu)建切花菊F1代性狀綜合評價體系,結(jié)果發(fā)現(xiàn)花序直徑、花色和花型3個評價因子對切花菊觀賞性的影響最大,為培育彩色標準型切花菊提供了材料。韓勇等[12]基于AHP對多頭切花菊的品質(zhì)性狀進行了綜合評價,發(fā)現(xiàn)對多頭切花菊品質(zhì)影響最大的性狀是株型,其次是花色和分枝長度,權(quán)重值分別為0.228、0.150、0.132。本研究對茶用菊的10個品質(zhì)性狀指標量化分析的結(jié)果表明,單株花數(shù)、田間適應性、舌狀花數(shù)、花重瓣性、花朵在枝條上的分布和收獲期的直立性等5個指標對于茶用菊綜合水平影響較大,其綜合權(quán)重值依次為0.301、0.172、0.119、0.113、0.088、0.066。這些指標在今后的茶用菊選種過程中應重點關注,單株花數(shù)和田間適應性與茶用菊的產(chǎn)量直接相關,而舌狀花數(shù)和花重瓣性與茶用菊的外觀品質(zhì)直接相關,體現(xiàn)了這一評價體系注重茶用菊產(chǎn)量和品質(zhì)的合理性,符合預期結(jié)果。本研究篩選出的品種與田間目測長勢結(jié)果相一致,該評價體系可以為茶用菊品種生產(chǎn)和選種提供指導。

    本研究從生長特性、花部特征和適應性3個方面出發(fā),結(jié)合田間試驗所得到的各個指標,利用層次分析法構(gòu)建出茶用菊綜合評價體系。應用層次分析法計算得到茶用菊10個評價指標的權(quán)重值,其中單株花數(shù)、田間適應性、舌狀花數(shù)的權(quán)重值最大,其次是花重瓣性、花朵在枝條上的分布和收獲期的直立性,這些指標的權(quán)重值遠遠大于其他指標的權(quán)重值。確定各品種的綜合評分后通過K-means聚類分析法將40個品種劃分為4個等級,其中,1級優(yōu)秀的株系共有5個(12.5%);2級良好的品系共有14個(35.0%);3級一般的株系共有13個(32.5%);4級差等株系共有8個(20.0%)。篩選出的優(yōu)良茶用菊品種田間適應性強,幾乎無病害發(fā)生,花朵繁密,產(chǎn)量高,在枝頭分布集中、易采收,株形直立緊湊、不易倒伏,具有產(chǎn)業(yè)化推廣的潛力。

    猜你喜歡
    品系適應性單株
    貴州黑山羊新品系選育
    10個團豆新品系在綏陽縣的田間性狀及產(chǎn)量表現(xiàn)
    谷子引種適應性鑒定與篩選初報
    無為市太平山楓香樹不同單株葉片性狀多樣性分析
    4個地被菊新品系對濕熱脅迫的耐受性研究
    園林科技(2020年2期)2020-01-18 03:28:18
    健全現(xiàn)代金融體系的適應性之“點論”
    中國外匯(2019年23期)2019-05-25 07:06:20
    種植密度與行距對秋閑田飼用甜高粱單株生產(chǎn)力的影響
    湖南速生、中生、慢生闊葉樹組單株生長模型構(gòu)建
    大型飛機A380-800在既有跑道起降的適應性研究
    固有免疫和適應性免疫與慢性丙肝的研究進展
    日日干狠狠操夜夜爽| 97超级碰碰碰精品色视频在线观看| 国产麻豆69| 在线观看免费日韩欧美大片| 欧美黄色淫秽网站| 日本在线视频免费播放| 国产免费男女视频| 欧美+亚洲+日韩+国产| 国产视频一区二区在线看| 精品电影一区二区在线| 亚洲国产日韩欧美精品在线观看 | 精品国产超薄肉色丝袜足j| 麻豆久久精品国产亚洲av| 亚洲欧美激情在线| 亚洲精品在线观看二区| 一级黄色大片毛片| 欧美日本亚洲视频在线播放| 国产精品二区激情视频| 在线观看午夜福利视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩无卡精品| 亚洲欧美日韩无卡精品| 欧美激情高清一区二区三区| av福利片在线| 高潮久久久久久久久久久不卡| 精品国内亚洲2022精品成人| 9191精品国产免费久久| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 久久亚洲真实| 午夜成年电影在线免费观看| 亚洲午夜精品一区,二区,三区| e午夜精品久久久久久久| 香蕉国产在线看| 性色av乱码一区二区三区2| 99热只有精品国产| 欧美午夜高清在线| 亚洲国产精品sss在线观看| 欧美日韩精品网址| 亚洲 欧美一区二区三区| 精品久久蜜臀av无| 波多野结衣一区麻豆| 操出白浆在线播放| 欧美色欧美亚洲另类二区 | 99久久99久久久精品蜜桃| 黄色成人免费大全| 久久精品91蜜桃| 久久精品国产综合久久久| 午夜亚洲福利在线播放| 日本精品一区二区三区蜜桃| 日本精品一区二区三区蜜桃| 老汉色av国产亚洲站长工具| 中文字幕人成人乱码亚洲影| 亚洲国产日韩欧美精品在线观看 | 少妇的丰满在线观看| 十分钟在线观看高清视频www| 中文字幕另类日韩欧美亚洲嫩草| 两个人免费观看高清视频| 日本欧美视频一区| 夜夜躁狠狠躁天天躁| 国产熟女xx| 成人av一区二区三区在线看| 午夜视频精品福利| 久久精品国产亚洲av香蕉五月| 日本免费a在线| 美女 人体艺术 gogo| 午夜免费鲁丝| 国产欧美日韩一区二区三区在线| 无遮挡黄片免费观看| 亚洲人成77777在线视频| 91成年电影在线观看| 亚洲人成77777在线视频| 叶爱在线成人免费视频播放| 如日韩欧美国产精品一区二区三区| 国产精品亚洲av一区麻豆| 99国产精品99久久久久| 中文字幕久久专区| 国产99白浆流出| 国产一卡二卡三卡精品| 日韩大码丰满熟妇| 国产亚洲欧美精品永久| 国产精品98久久久久久宅男小说| 在线观看一区二区三区| 这个男人来自地球电影免费观看| 亚洲国产中文字幕在线视频| 在线观看66精品国产| 午夜福利欧美成人| 亚洲五月天丁香| 亚洲精品国产一区二区精华液| 久久久久国产精品人妻aⅴ院| 在线视频色国产色| 亚洲精品中文字幕一二三四区| a在线观看视频网站| 成人亚洲精品av一区二区| 亚洲专区字幕在线| 国产真人三级小视频在线观看| 国产伦人伦偷精品视频| 国产成人欧美在线观看| 日韩大尺度精品在线看网址 | 精品免费久久久久久久清纯| 国产激情久久老熟女| 美女国产高潮福利片在线看| 男女做爰动态图高潮gif福利片 | 久久人妻熟女aⅴ| 免费在线观看亚洲国产| 91在线观看av| 精品人妻1区二区| 亚洲人成77777在线视频| 亚洲人成伊人成综合网2020| 91在线观看av| 波多野结衣av一区二区av| 国产真人三级小视频在线观看| 天天添夜夜摸| 69精品国产乱码久久久| 老司机深夜福利视频在线观看| 久久久国产欧美日韩av| 满18在线观看网站| xxx96com| 国产精品影院久久| 久久精品国产亚洲av高清一级| 十分钟在线观看高清视频www| 国产亚洲精品一区二区www| 97人妻精品一区二区三区麻豆 | 亚洲av成人av| 最新在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲国产精品999在线| 国产亚洲精品久久久久5区| 女人精品久久久久毛片| 欧美激情久久久久久爽电影 | 日韩欧美一区二区三区在线观看| 久久精品91无色码中文字幕| 国产精品电影一区二区三区| 91成人精品电影| 亚洲av美国av| 亚洲av第一区精品v没综合| 精品久久蜜臀av无| 亚洲人成网站在线播放欧美日韩| 深夜精品福利| 日本vs欧美在线观看视频| 精品欧美一区二区三区在线| 免费一级毛片在线播放高清视频 | 91九色精品人成在线观看| www.精华液| 成熟少妇高潮喷水视频| 男人操女人黄网站| 黄片大片在线免费观看| 久久久久国内视频| 99在线视频只有这里精品首页| 极品人妻少妇av视频| 精品欧美一区二区三区在线| 91大片在线观看| 亚洲精品av麻豆狂野| 久久精品亚洲熟妇少妇任你| 又紧又爽又黄一区二区| 国产一区二区三区视频了| 国产av在哪里看| 97人妻天天添夜夜摸| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜视频精品福利| 午夜福利,免费看| 天堂√8在线中文| 老司机深夜福利视频在线观看| 日韩精品中文字幕看吧| 一进一出抽搐gif免费好疼| 国产又爽黄色视频| АⅤ资源中文在线天堂| 国产亚洲精品久久久久5区| 久久天堂一区二区三区四区| 国产成人系列免费观看| 精品国产一区二区三区四区第35| 日韩三级视频一区二区三区| 亚洲成人久久性| 又紧又爽又黄一区二区| 大型黄色视频在线免费观看| 欧美+亚洲+日韩+国产| 国产主播在线观看一区二区| 久久狼人影院| 亚洲欧美日韩另类电影网站| 国产欧美日韩一区二区三区在线| 午夜精品国产一区二区电影| 国产三级在线视频| 亚洲成av片中文字幕在线观看| 亚洲av五月六月丁香网| 亚洲欧美日韩无卡精品| 色播亚洲综合网| 伦理电影免费视频| 国产成人欧美| 最好的美女福利视频网| 亚洲av成人不卡在线观看播放网| 亚洲国产日韩欧美精品在线观看 | 欧美激情高清一区二区三区| 久久国产亚洲av麻豆专区| xxx96com| 国产伦一二天堂av在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产熟女xx| 国产精品99久久99久久久不卡| 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 国产精品电影一区二区三区| 中文字幕久久专区| 亚洲自偷自拍图片 自拍| 欧美性长视频在线观看| 国产99白浆流出| 免费看a级黄色片| 亚洲美女黄片视频| 99re在线观看精品视频| 久热这里只有精品99| 天堂√8在线中文| ponron亚洲| 男人舔女人的私密视频| 99国产精品一区二区三区| 国产精品99久久99久久久不卡| 成人国语在线视频| 一进一出抽搐动态| 99精品久久久久人妻精品| 91国产中文字幕| 国产97色在线日韩免费| 啦啦啦韩国在线观看视频| 国产高清有码在线观看视频 | 国产高清videossex| 欧美黄色淫秽网站| e午夜精品久久久久久久| 大型av网站在线播放| 国产极品粉嫩免费观看在线| 国产色视频综合| 午夜福利免费观看在线| 国产精品日韩av在线免费观看 | 国产黄a三级三级三级人| 欧美成人一区二区免费高清观看 | av中文乱码字幕在线| 亚洲va日本ⅴa欧美va伊人久久| 99久久久亚洲精品蜜臀av| 亚洲欧美激情在线| 亚洲激情在线av| 可以在线观看的亚洲视频| 99国产精品一区二区三区| 91麻豆精品激情在线观看国产| 午夜精品久久久久久毛片777| 亚洲国产欧美网| 久久狼人影院| 无人区码免费观看不卡| 国产精品免费视频内射| 午夜福利欧美成人| 国产伦人伦偷精品视频| 人人妻人人澡人人看| aaaaa片日本免费| 麻豆国产av国片精品| 亚洲男人天堂网一区| 欧美午夜高清在线| 久久国产精品男人的天堂亚洲| 天堂√8在线中文| 色播亚洲综合网| 亚洲人成电影观看| 中国美女看黄片| 国产亚洲欧美精品永久| 免费在线观看日本一区| 久久香蕉精品热| 在线免费观看的www视频| 中文字幕久久专区| 亚洲五月婷婷丁香| 免费看美女性在线毛片视频| 性欧美人与动物交配| 精品久久久精品久久久| 12—13女人毛片做爰片一| 国产一卡二卡三卡精品| 色播亚洲综合网| 老司机午夜十八禁免费视频| 国内精品久久久久精免费| 黑人巨大精品欧美一区二区mp4| 亚洲第一av免费看| 女性生殖器流出的白浆| 亚洲五月婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产中文字幕在线视频| 曰老女人黄片| 韩国精品一区二区三区| 久久亚洲精品不卡| 久久人妻熟女aⅴ| 欧美日韩福利视频一区二区| 女性被躁到高潮视频| netflix在线观看网站| 亚洲 欧美 日韩 在线 免费| 国产xxxxx性猛交| 欧美日韩乱码在线| 国产亚洲欧美98| 可以在线观看毛片的网站| 在线观看66精品国产| 99re在线观看精品视频| 99国产综合亚洲精品| 亚洲一区高清亚洲精品| 窝窝影院91人妻| 99国产精品一区二区三区| 此物有八面人人有两片| 88av欧美| 一级片免费观看大全| 日本撒尿小便嘘嘘汇集6| 精品国内亚洲2022精品成人| 中文亚洲av片在线观看爽| 久热这里只有精品99| 欧美成人一区二区免费高清观看 | 成人av一区二区三区在线看| 国产欧美日韩精品亚洲av| 久久久久国内视频| 成人亚洲精品av一区二区| 韩国av一区二区三区四区| 一进一出抽搐gif免费好疼| 国产精品,欧美在线| 亚洲一码二码三码区别大吗| 日本撒尿小便嘘嘘汇集6| 国产99久久九九免费精品| 亚洲av五月六月丁香网| av天堂在线播放| 久久久国产精品麻豆| 真人一进一出gif抽搐免费| 亚洲国产精品久久男人天堂| 女人爽到高潮嗷嗷叫在线视频| 欧美日本亚洲视频在线播放| 91精品三级在线观看| 18禁美女被吸乳视频| 久9热在线精品视频| 午夜久久久久精精品| 色老头精品视频在线观看| av中文乱码字幕在线| 免费av毛片视频| 黑丝袜美女国产一区| 可以在线观看的亚洲视频| 国产av一区在线观看免费| 日本免费a在线| 日韩一卡2卡3卡4卡2021年| 国产成人av激情在线播放| 中文亚洲av片在线观看爽| 欧美激情极品国产一区二区三区| 欧美一级毛片孕妇| 悠悠久久av| 亚洲狠狠婷婷综合久久图片| 国产精品秋霞免费鲁丝片| 手机成人av网站| 天堂动漫精品| aaaaa片日本免费| 亚洲一区中文字幕在线| avwww免费| 欧美成人一区二区免费高清观看 | 欧美色视频一区免费| 亚洲av成人一区二区三| 丁香欧美五月| 久久久久久久精品吃奶| 日本三级黄在线观看| 欧美亚洲日本最大视频资源| 国产一区二区在线av高清观看| 很黄的视频免费| 欧美国产精品va在线观看不卡| 人人妻人人爽人人添夜夜欢视频| 精品欧美国产一区二区三| 我的亚洲天堂| 波多野结衣高清无吗| 一个人观看的视频www高清免费观看 | 精品久久久精品久久久| 此物有八面人人有两片| 国产av精品麻豆| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区中文字幕在线| 男女下面进入的视频免费午夜 | 熟妇人妻久久中文字幕3abv| 国产不卡一卡二| 一个人观看的视频www高清免费观看 | 亚洲最大成人中文| 日本五十路高清| 午夜亚洲福利在线播放| 黄色丝袜av网址大全| 欧美精品亚洲一区二区| 日韩欧美三级三区| 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 大型黄色视频在线免费观看| 69精品国产乱码久久久| 精品久久久久久成人av| 少妇粗大呻吟视频| 18禁美女被吸乳视频| 国产成人av激情在线播放| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 91老司机精品| 乱人伦中国视频| 亚洲欧美一区二区三区黑人| 国产片内射在线| 纯流量卡能插随身wifi吗| 午夜亚洲福利在线播放| 国产av一区在线观看免费| 午夜福利高清视频| 国产精品av久久久久免费| 校园春色视频在线观看| 妹子高潮喷水视频| 欧美激情高清一区二区三区| 黄色丝袜av网址大全| 女人被躁到高潮嗷嗷叫费观| videosex国产| 亚洲av电影不卡..在线观看| 手机成人av网站| 不卡一级毛片| 国产真人三级小视频在线观看| 久久久水蜜桃国产精品网| 日韩av在线大香蕉| 操出白浆在线播放| 欧美乱码精品一区二区三区| 老司机在亚洲福利影院| 无遮挡黄片免费观看| 日韩大码丰满熟妇| 国产精品亚洲美女久久久| 美女午夜性视频免费| 亚洲精品中文字幕一二三四区| 国产精品永久免费网站| 亚洲精品国产区一区二| 天堂√8在线中文| 国产成人一区二区三区免费视频网站| 在线观看舔阴道视频| 国产成人系列免费观看| 久久午夜亚洲精品久久| 亚洲国产日韩欧美精品在线观看 | 日韩成人在线观看一区二区三区| 1024视频免费在线观看| 精品第一国产精品| 成人国产一区最新在线观看| av在线天堂中文字幕| 大码成人一级视频| 国产精品亚洲一级av第二区| 久久精品影院6| 女人精品久久久久毛片| 不卡一级毛片| 亚洲精品中文字幕一二三四区| 久久久久久久久免费视频了| 亚洲欧美日韩无卡精品| 桃红色精品国产亚洲av| 欧美乱色亚洲激情| 久久久久久国产a免费观看| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 91成年电影在线观看| 美女扒开内裤让男人捅视频| 精品国产超薄肉色丝袜足j| 亚洲av成人av| 一卡2卡三卡四卡精品乱码亚洲| 麻豆成人av在线观看| 国产亚洲欧美在线一区二区| 亚洲avbb在线观看| 精品无人区乱码1区二区| 91国产中文字幕| 亚洲色图综合在线观看| 最好的美女福利视频网| 欧美成人一区二区免费高清观看 | 性色av乱码一区二区三区2| 国产单亲对白刺激| 女同久久另类99精品国产91| 国产成人精品在线电影| 国产单亲对白刺激| 久久国产精品影院| 天堂√8在线中文| av片东京热男人的天堂| 国产熟女xx| 一级a爱片免费观看的视频| 午夜免费激情av| 久久久久国产精品人妻aⅴ院| 老汉色av国产亚洲站长工具| 最新在线观看一区二区三区| 变态另类成人亚洲欧美熟女 | 国内精品久久久久久久电影| 国产亚洲精品综合一区在线观看 | 欧美成狂野欧美在线观看| 免费在线观看日本一区| 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 99精品在免费线老司机午夜| av天堂在线播放| 亚洲黑人精品在线| 高清在线国产一区| 午夜精品国产一区二区电影| 亚洲国产欧美一区二区综合| 黄色a级毛片大全视频| 国产野战对白在线观看| 亚洲国产欧美一区二区综合| 亚洲色图综合在线观看| 国产一区在线观看成人免费| 99精品久久久久人妻精品| 超碰成人久久| 亚洲第一青青草原| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 国产精品爽爽va在线观看网站 | 久久精品91蜜桃| 精品久久久久久久毛片微露脸| 午夜久久久久精精品| 色在线成人网| 亚洲精华国产精华精| 人妻丰满熟妇av一区二区三区| 精品人妻1区二区| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 1024香蕉在线观看| 欧美黄色淫秽网站| 成人特级黄色片久久久久久久| 欧美日韩乱码在线| 真人做人爱边吃奶动态| 日韩 欧美 亚洲 中文字幕| a在线观看视频网站| 麻豆一二三区av精品| 18禁观看日本| 亚洲狠狠婷婷综合久久图片| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 美女 人体艺术 gogo| 不卡av一区二区三区| 免费高清视频大片| 国产极品粉嫩免费观看在线| 亚洲av五月六月丁香网| 国产午夜精品久久久久久| 亚洲国产中文字幕在线视频| 成人国产一区最新在线观看| 午夜免费激情av| 高潮久久久久久久久久久不卡| 97碰自拍视频| 亚洲精品久久国产高清桃花| 欧美日本视频| 91成年电影在线观看| 黄色丝袜av网址大全| 女人高潮潮喷娇喘18禁视频| 禁无遮挡网站| 国产精品99久久99久久久不卡| 国产午夜精品久久久久久| svipshipincom国产片| 自线自在国产av| 精品国产乱码久久久久久男人| 成人欧美大片| 在线永久观看黄色视频| 性欧美人与动物交配| 女人被躁到高潮嗷嗷叫费观| 最近最新免费中文字幕在线| 国产亚洲精品av在线| 国产精品二区激情视频| 性欧美人与动物交配| 99久久久亚洲精品蜜臀av| 国产成人精品久久二区二区91| 高清在线国产一区| 99re在线观看精品视频| 极品人妻少妇av视频| 国产一级毛片七仙女欲春2 | 亚洲男人天堂网一区| 日本在线视频免费播放| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美98| 免费看a级黄色片| x7x7x7水蜜桃| 97碰自拍视频| 国产麻豆成人av免费视频| 精品一区二区三区视频在线观看免费| 中亚洲国语对白在线视频| 岛国在线观看网站| 国产高清有码在线观看视频 | 99国产精品99久久久久| 免费高清视频大片| 国产97色在线日韩免费| 国产精品综合久久久久久久免费 | 18禁裸乳无遮挡免费网站照片 | 岛国视频午夜一区免费看| 女人被狂操c到高潮| av中文乱码字幕在线| 欧美黄色淫秽网站| 性色av乱码一区二区三区2| 国产亚洲精品综合一区在线观看 | 精品第一国产精品| 国产一区在线观看成人免费| 美女 人体艺术 gogo| 午夜福利高清视频| 欧美日韩福利视频一区二区| 午夜两性在线视频| 亚洲国产欧美一区二区综合| 亚洲一区中文字幕在线| 国产av精品麻豆| 成人欧美大片| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 搡老岳熟女国产| 熟女少妇亚洲综合色aaa.| 免费看a级黄色片| 国产欧美日韩综合在线一区二区| 不卡一级毛片| 日本一区二区免费在线视频| 黄网站色视频无遮挡免费观看| 电影成人av| 国产亚洲精品一区二区www| 亚洲激情在线av| bbb黄色大片| 又紧又爽又黄一区二区| 色精品久久人妻99蜜桃| 搡老熟女国产l中国老女人| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 日韩三级视频一区二区三区| 制服人妻中文乱码| 此物有八面人人有两片| 久久这里只有精品19| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| 午夜精品国产一区二区电影| 国产精品秋霞免费鲁丝片| 老汉色∧v一级毛片| 又大又爽又粗| 侵犯人妻中文字幕一二三四区| 久久久久国内视频| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频| 亚洲精品av麻豆狂野| 一个人免费在线观看的高清视频| 成人精品一区二区免费|