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

    南方春大豆品種秋播條件下的綜合評(píng)價(jià)及聚類分析

    2014-05-30 16:50李清華等
    安徽農(nóng)學(xué)通報(bào) 2014年8期
    關(guān)鍵詞:秋播聚類分析

    李清華等

    摘 要:通過測(cè)定17個(gè)南方春大豆品種在莆田地區(qū)秋播下的生物學(xué)性狀,變異分析結(jié)果表明,17個(gè)春大豆品種底莢高度變異系數(shù)最大;相關(guān)性分析表明,春大豆秋播條件下,單株粒重與單株粒數(shù)呈顯著正相關(guān);主成分分析結(jié)果表明,前4個(gè)主成分的總方差遺傳累積貢獻(xiàn)率達(dá)84.43%,基本上反映了所有性狀的變異信息;系統(tǒng)聚類分析表明,17份品種主要可分為有籽粒大、中、小,且產(chǎn)量高、中、低交叉等7大類型。

    關(guān)鍵詞:春大豆;秋播;主成份分析;聚類分析

    中圖分類號(hào) S565.1 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731(2014)08-52-04

    Abstract:Tested the main biological characters of 17 south spring soybean varieties sowed in autumn in Putian region. The results showed that the variance coefficient of height of pod was the largest. The correlation analysis indicated that seed weight per plant positive significant correlation with number of seeds per plant when south spring soybean sowed in autumn;Principal components analysis showed that the first four principal components variance accumulation contribution rate amounted to 84.43%,which basically reflected the variance information as listed of all characteristics above. Cluster analysis showed that 17 south spring soybean could be divided into 7 types which included big、medium and little seed,and high、medium and low yield,and so on.

    Key words:Spring soybean;Autumn sowing;Principal component analysis;Cluster analysis

    大豆是糧油兼用作物,是人們獲取植物蛋白的主要來源,也是發(fā)展畜牧業(yè)的重要精飼料和許多輕工業(yè)的重要原料。福建省屬亞熱帶濕潤(rùn)季風(fēng)區(qū),大豆栽培歷史悠久,適宜一年兩播,但是由于不同播種季節(jié)所處的氣候環(huán)境存在一定差異,造成不同播種季節(jié)下種植產(chǎn)量、效益的差距。春播大豆生育期一般處于3月上旬至7月上旬,而在5、6月份結(jié)莢鼓粒期至成熟期福建地區(qū)多臺(tái)風(fēng)暴雨,易造成大豆落花落莢、植株倒伏、爛莢爛粒等現(xiàn)象。秋播大豆生育期一般處于7月下旬至10月中下旬,生育前期多高溫,植株生長(zhǎng)過快,導(dǎo)致后期植株變矮、粒重下降。以上不利因素導(dǎo)致了大豆產(chǎn)量降低、效益減少,進(jìn)而影響了農(nóng)民群眾種植春大豆積極性,導(dǎo)致福建省春大豆種植面積的不斷減少。

    近年來,數(shù)量分類學(xué)在農(nóng)作物品種研究上得到廣泛應(yīng)用,科學(xué)準(zhǔn)確地鑒定篩選新品種,是大豆育種、推廣程序中的一個(gè)重要環(huán)節(jié)。特別是主成分分析和聚類分析在花生[1-2]、玉米[3]、水稻[4]、甘薯[5]等作物種質(zhì)資源分類方面的應(yīng)用,對(duì)于科學(xué)評(píng)價(jià)品系好壞真實(shí)性、品種適應(yīng)性起到了較好的效果,同時(shí)在大豆育種中品系鑒定和品種推廣應(yīng)用也有較多的報(bào)道[6-8]。本試驗(yàn)通過借鑒以往研究方法的基礎(chǔ)上,采用主成分分析和聚類分析方法對(duì)17份南方春大豆品種在莆田地區(qū)秋種的適應(yīng)性進(jìn)行評(píng)價(jià),以期為福建省秋季大豆種植品種更新和篩選及春大豆育種中就地秋播加代對(duì)性狀有效選擇提供理論依據(jù),為大豆品種的區(qū)域適應(yīng)性評(píng)價(jià)提供方法參考。

    1 材料與方法

    1.1 試驗(yàn)材料 本試驗(yàn)的材料包括有:福建省莆田市農(nóng)業(yè)科學(xué)研究所育成的莆豆8008、莆豆10號(hào)和莆豆5號(hào);福建省農(nóng)科院作物所育成的福豆310和福豆234;福建省泉州市農(nóng)業(yè)科學(xué)研究所育成的泉豆6號(hào)和泉豆7號(hào)。福建省華南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院育成的華春1號(hào)、華春2號(hào)、華春3號(hào)、華春5號(hào)和華春6號(hào);廣西壯族自治區(qū)玉米研究所育成的桂春1號(hào)、桂春2號(hào)、桂春5號(hào)、桂春6號(hào)和桂春8號(hào)。

    1.2 試驗(yàn)設(shè)計(jì) 試驗(yàn)在福建省莆田市農(nóng)業(yè)科學(xué)研究所的試驗(yàn)基地進(jìn)行,該地塊土層深厚,地勢(shì)平坦,肥力均勻中上,半粘壤土,排灌方便,前茬水稻,多為水稻-蔬菜輪作。2013年8月6日播種,小區(qū)隨機(jī)區(qū)組排列,3次重復(fù),小區(qū)面積6.67m2(6.67m×1m)。每小區(qū)雙行穴播,行距0.40m,株距0.238m,種植密度16 800株/667m2。

    1.3 數(shù)據(jù)收集與處理 觀察與考種記載項(xiàng)目包括全生育期(T)、主莖高(X1)、主莖節(jié)數(shù)(X2)、有效分枝(X3)、底莢高度(X4)、單數(shù)莢數(shù)(X5)、單株粒數(shù)(X6)、百粒重(X7)、單株粒重(X8)。采用DPS軟件進(jìn)行數(shù)據(jù)處理[9]。數(shù)據(jù)處理包括變異分析、相關(guān)性分析、主成分分析和聚類分析。聚類分析先按性狀組合進(jìn)行聚類,然后進(jìn)行綜合性狀聚類;聚類時(shí)將原始數(shù)據(jù)先進(jìn)行標(biāo)準(zhǔn)化變換,然后在歐氏距離水平上采用離差平方和法進(jìn)行系統(tǒng)聚類分析[10]。

    2 結(jié)果與分析

    2.1 生物學(xué)性狀的變異分析 試驗(yàn)材料的9個(gè)主要觀測(cè)值及變異情況見表1。由表1可以看出,有效分枝、底莢高度和單株粒數(shù)變異系數(shù)較大,分別為21.76%、25.21%、14.16%,變異明顯;主莖高、主莖節(jié)數(shù)、單株莢數(shù)、百粒重和單株粒重變異系數(shù)分別為7.21%、6.02%、8.52%、9.92%、6.63%,變異較明顯;生育期變異系數(shù)最小,僅為4.06%。endprint

    2.2 生物學(xué)性狀的相關(guān)性分析 對(duì)17分春大豆品種主要農(nóng)藝性狀的相關(guān)分析表明(表2):?jiǎn)沃炅?shù)與單株莢數(shù)呈極顯著正相關(guān),百粒重與單株莢數(shù)、單株粒數(shù)呈極顯著負(fù)相關(guān);單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)、單株莢數(shù)及單株粒數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān);主莖高與單株粒數(shù)呈負(fù)相關(guān),與百粒重呈正相關(guān);主莖節(jié)數(shù)與底莢高度呈負(fù)相關(guān)、與單株粒數(shù)呈正相關(guān);有效分枝與底莢高度呈正相關(guān)。

    2.3 主成分分析 對(duì)供試的17份春大豆品種的9個(gè)主要觀測(cè)值采用DPS數(shù)據(jù)處理軟件進(jìn)行主成分分析,計(jì)算得特征值Yi、特征值累計(jì)百分率和主成份特征向量。各特征值大小代表各綜合指標(biāo)遺傳方差大小,各特征值累計(jì)百分率代表各有關(guān)綜合指標(biāo)對(duì)總遺傳方差貢獻(xiàn)的百分率。為提高分析精度,只選用對(duì)總遺傳方差貢獻(xiàn)率達(dá)85.43%的前4個(gè)特征值(表3)。

    從表3可以看出,第一主成分的特征向量中,單株莢數(shù)、單株粒數(shù)和單株粒重的權(quán)重系數(shù)為正值,百粒重的權(quán)重系數(shù)為負(fù)值,均較大,構(gòu)成了變異主要來源,說明單株莢數(shù)、單株粒數(shù)與單株粒重呈同步增長(zhǎng)態(tài)勢(shì),而單株莢數(shù)、單株粒數(shù)同時(shí)影響了籽粒的大小,故可稱為產(chǎn)量關(guān)鍵因子;第二主成分的特征向量中,有效分枝、底莢高度權(quán)重系數(shù)為正值,百粒重、單株粒重權(quán)重系數(shù)為負(fù)值,且較大,構(gòu)成了變異主要來源,說明有效分枝過多、底莢高度過高將影響大豆籽粒大小,進(jìn)而影響了大豆的單株產(chǎn)量,故可稱為產(chǎn)量影響因子;第三主成分的特征向量中,主莖高、主莖節(jié)數(shù)和有效分枝權(quán)重系數(shù)均為正數(shù),且較大,這些主要反映了植株生長(zhǎng)勢(shì),故稱為生長(zhǎng)勢(shì)因子;第四主成份的特征向量中,生育期與主莖高、底莢高度權(quán)重系數(shù)較大,可稱之為生育期與植株高度綜合因子。

    2.4 聚類分析

    2.4.1 不同性狀組合的聚類分析 對(duì)參試品種分別按植株性狀與經(jīng)濟(jì)性狀2種性狀組合分別進(jìn)行少類性狀聚類。其中植株性狀包括:主莖高、主莖節(jié)數(shù)、有效分枝與底莢高度;經(jīng)濟(jì)性狀包括:?jiǎn)沃昵v數(shù)、單株粒數(shù)、百粒重與單株粒重。聚類結(jié)果表明:參與聚類的性狀組合不同,其聚類結(jié)果各異。按植株性狀方面聚類(圖1),17個(gè)品種在聚類水平D2=1.07時(shí)可分成12類,類型較為多樣復(fù)雜;按經(jīng)濟(jì)性狀方面聚類(圖2),17個(gè)品種在聚類水平D2=4.06時(shí)只分成了四大類,且從圖中可以看出福建省育成的多數(shù)春大豆品種在經(jīng)濟(jì)性狀方面較為相似。

    2.4.2 綜合性狀聚類分析 對(duì)參試品種9個(gè)觀測(cè)值進(jìn)行標(biāo)準(zhǔn)化變換后在歐氏距離水平上采用離差平方和法進(jìn)行系統(tǒng)聚類分析(圖3),根據(jù)統(tǒng)一類群內(nèi)大豆品種類群間距離接近且綜合性狀差異值較小的原則,17份大豆品種在聚類水平D2=4.24時(shí)可以劃分為7大類群。

    第1類群有1、2和7號(hào)3個(gè)品種,該類群品種有效分枝、底莢高度兩性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于單株莢數(shù)、粒數(shù)偏少,籽粒較大的低產(chǎn)型品種;第2類群有4、5、6和17號(hào)4個(gè)品種,該類群品種在植株性狀和單株莢數(shù)存在明顯差異,其他性狀均表現(xiàn)相近,屬于生育期較長(zhǎng),粒少、粒大的中產(chǎn)型品種;第3類群有3和9號(hào)2個(gè)品種,兩品種在產(chǎn)量性狀表現(xiàn)上存在較小差異,均屬于生育期適中,植株旺盛,結(jié)莢數(shù)適宜,粒少、粒大的中高產(chǎn)型品種;第4類群有8、16、11和15號(hào)4個(gè)品種,該類群生育期、主莖節(jié)數(shù)和有效分枝性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于底莢高度較低,莢多粒多,籽粒較小的高產(chǎn)型品種;第5類群有10號(hào)1個(gè)品種,該品種屬于株高適中,分枝旺盛,底莢高度較高,粒多、粒小的早熟低產(chǎn)型品種;第6類群12和14號(hào)3個(gè)品種,該類群品種屬于生育期較長(zhǎng),底莢高度偏高,籽粒大小中等,其他性狀均表現(xiàn)優(yōu)異的中產(chǎn)型品種。第7類群有13號(hào)1個(gè)品種,該品種屬于生育期偏長(zhǎng),主莖節(jié)數(shù)少,底莢高度高,莢數(shù)粒數(shù)多,籽粒小的低產(chǎn)型品種。

    從不同角度進(jìn)行的聚類分析結(jié)果可以看出,綜合性狀的聚類與兩種性狀組合聚類的結(jié)果存在明顯不同。少數(shù)性狀組合的聚類只考慮到少數(shù)性狀表現(xiàn)的相似性,各類群內(nèi)只有被聚類的性狀表現(xiàn)一致,而大多性狀表現(xiàn)明顯差異。綜合性狀的聚類使類群內(nèi)大多性狀表現(xiàn)一致,僅個(gè)別性狀有差異,結(jié)果更具整體性和合理性[3]。

    3 結(jié)論和討論

    本研究結(jié)果表明,供試的17份南方春大豆品種秋播下的生物學(xué)性狀變異系數(shù)以底莢高度最大、生育期最小。單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)及單株莢數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān)。主成份分析結(jié)果表明,秋播條件下不同大豆品種的生物學(xué)性狀可分為產(chǎn)量關(guān)鍵因子、產(chǎn)量影響因子、生長(zhǎng)勢(shì)因子和生育期與株高綜合因子4個(gè)主成分,它們基本上代表了所有性狀的變異信息。系統(tǒng)聚類分析結(jié)果表明,17份品種類群較為復(fù)雜,其中有莢少粒少、粒大的低產(chǎn)型品種,粒少、粒大的中高產(chǎn)型品種,莢多粒多,籽小的高產(chǎn)型品種等;同時(shí)聚類分析結(jié)果還表明了來自同一單位育成的品種不一定類型相似。總之,主成分分析結(jié)果反映了17份春大豆秋播下生物學(xué)性狀的主導(dǎo)影響因子,聚類結(jié)果反映了產(chǎn)量的差異,這為正確評(píng)價(jià)大豆品種的適應(yīng)性及科學(xué)地選擇品種提供了依據(jù)[6]。

    參考文獻(xiàn)

    [1]殷冬梅,李銓柱,崔黨群.花生主要農(nóng)藝性狀的相關(guān)性及聚類分析[J].中國(guó)油料作物學(xué)報(bào),2010,32(2):212-216.

    [2]李清華,黃金堂,陳海玲,等.27份花生種質(zhì)資源的主成份分析與遺傳距離測(cè)定[J].植物遺傳資源學(xué)報(bào),2011,12(4):519-524.

    [3]李淑華,荊紹凌,楊丹萍,等.聚類分析法在玉米種質(zhì)資源中的應(yīng)用[J].農(nóng)業(yè)技術(shù),2008,28(5):43-45.

    [4]馮永祥,劉沐江,汪秀志,等.寒地水稻產(chǎn)量性狀相關(guān)及聚類分析[J].安徽農(nóng)業(yè)科學(xué),2008,36(23):9 925-9 928.

    [5]蘇春華,李育明,黃迎冬,等.甘薯主要親本材料的主成分分析及聚類分析[J].雜糧作物,2007,27(6):405-409.

    [6]張玉革,胡緒彬.基于主成分和聚類分析的大豆品種生物學(xué)性狀的比較研究[J].大豆科學(xué),2004,23(3):178-183[J].

    [7]李向華,常釹鎮(zhèn).中國(guó)春大豆品種聚類分析及主成分分析[J].作物學(xué)報(bào),1998,24(3):325-332.

    [8]劉燦洪,林榮輝,劉德金.福建省大豆地方品種多元統(tǒng)計(jì)分析[J].福建農(nóng)業(yè)學(xué)報(bào),1999,14(4):11-15.

    [9] 唐啟義,馮光明. 實(shí)用統(tǒng)計(jì)分析及DPS其數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [10]王林輝.基于主成分分析的棉花品種綜合評(píng)價(jià)及聚類分析[J].廣東農(nóng)業(yè)科學(xué),2009,(1):29-32. (責(zé)編:張宏民)endprint

    2.2 生物學(xué)性狀的相關(guān)性分析 對(duì)17分春大豆品種主要農(nóng)藝性狀的相關(guān)分析表明(表2):?jiǎn)沃炅?shù)與單株莢數(shù)呈極顯著正相關(guān),百粒重與單株莢數(shù)、單株粒數(shù)呈極顯著負(fù)相關(guān);單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)、單株莢數(shù)及單株粒數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān);主莖高與單株粒數(shù)呈負(fù)相關(guān),與百粒重呈正相關(guān);主莖節(jié)數(shù)與底莢高度呈負(fù)相關(guān)、與單株粒數(shù)呈正相關(guān);有效分枝與底莢高度呈正相關(guān)。

    2.3 主成分分析 對(duì)供試的17份春大豆品種的9個(gè)主要觀測(cè)值采用DPS數(shù)據(jù)處理軟件進(jìn)行主成分分析,計(jì)算得特征值Yi、特征值累計(jì)百分率和主成份特征向量。各特征值大小代表各綜合指標(biāo)遺傳方差大小,各特征值累計(jì)百分率代表各有關(guān)綜合指標(biāo)對(duì)總遺傳方差貢獻(xiàn)的百分率。為提高分析精度,只選用對(duì)總遺傳方差貢獻(xiàn)率達(dá)85.43%的前4個(gè)特征值(表3)。

    從表3可以看出,第一主成分的特征向量中,單株莢數(shù)、單株粒數(shù)和單株粒重的權(quán)重系數(shù)為正值,百粒重的權(quán)重系數(shù)為負(fù)值,均較大,構(gòu)成了變異主要來源,說明單株莢數(shù)、單株粒數(shù)與單株粒重呈同步增長(zhǎng)態(tài)勢(shì),而單株莢數(shù)、單株粒數(shù)同時(shí)影響了籽粒的大小,故可稱為產(chǎn)量關(guān)鍵因子;第二主成分的特征向量中,有效分枝、底莢高度權(quán)重系數(shù)為正值,百粒重、單株粒重權(quán)重系數(shù)為負(fù)值,且較大,構(gòu)成了變異主要來源,說明有效分枝過多、底莢高度過高將影響大豆籽粒大小,進(jìn)而影響了大豆的單株產(chǎn)量,故可稱為產(chǎn)量影響因子;第三主成分的特征向量中,主莖高、主莖節(jié)數(shù)和有效分枝權(quán)重系數(shù)均為正數(shù),且較大,這些主要反映了植株生長(zhǎng)勢(shì),故稱為生長(zhǎng)勢(shì)因子;第四主成份的特征向量中,生育期與主莖高、底莢高度權(quán)重系數(shù)較大,可稱之為生育期與植株高度綜合因子。

    2.4 聚類分析

    2.4.1 不同性狀組合的聚類分析 對(duì)參試品種分別按植株性狀與經(jīng)濟(jì)性狀2種性狀組合分別進(jìn)行少類性狀聚類。其中植株性狀包括:主莖高、主莖節(jié)數(shù)、有效分枝與底莢高度;經(jīng)濟(jì)性狀包括:?jiǎn)沃昵v數(shù)、單株粒數(shù)、百粒重與單株粒重。聚類結(jié)果表明:參與聚類的性狀組合不同,其聚類結(jié)果各異。按植株性狀方面聚類(圖1),17個(gè)品種在聚類水平D2=1.07時(shí)可分成12類,類型較為多樣復(fù)雜;按經(jīng)濟(jì)性狀方面聚類(圖2),17個(gè)品種在聚類水平D2=4.06時(shí)只分成了四大類,且從圖中可以看出福建省育成的多數(shù)春大豆品種在經(jīng)濟(jì)性狀方面較為相似。

    2.4.2 綜合性狀聚類分析 對(duì)參試品種9個(gè)觀測(cè)值進(jìn)行標(biāo)準(zhǔn)化變換后在歐氏距離水平上采用離差平方和法進(jìn)行系統(tǒng)聚類分析(圖3),根據(jù)統(tǒng)一類群內(nèi)大豆品種類群間距離接近且綜合性狀差異值較小的原則,17份大豆品種在聚類水平D2=4.24時(shí)可以劃分為7大類群。

    第1類群有1、2和7號(hào)3個(gè)品種,該類群品種有效分枝、底莢高度兩性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于單株莢數(shù)、粒數(shù)偏少,籽粒較大的低產(chǎn)型品種;第2類群有4、5、6和17號(hào)4個(gè)品種,該類群品種在植株性狀和單株莢數(shù)存在明顯差異,其他性狀均表現(xiàn)相近,屬于生育期較長(zhǎng),粒少、粒大的中產(chǎn)型品種;第3類群有3和9號(hào)2個(gè)品種,兩品種在產(chǎn)量性狀表現(xiàn)上存在較小差異,均屬于生育期適中,植株旺盛,結(jié)莢數(shù)適宜,粒少、粒大的中高產(chǎn)型品種;第4類群有8、16、11和15號(hào)4個(gè)品種,該類群生育期、主莖節(jié)數(shù)和有效分枝性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于底莢高度較低,莢多粒多,籽粒較小的高產(chǎn)型品種;第5類群有10號(hào)1個(gè)品種,該品種屬于株高適中,分枝旺盛,底莢高度較高,粒多、粒小的早熟低產(chǎn)型品種;第6類群12和14號(hào)3個(gè)品種,該類群品種屬于生育期較長(zhǎng),底莢高度偏高,籽粒大小中等,其他性狀均表現(xiàn)優(yōu)異的中產(chǎn)型品種。第7類群有13號(hào)1個(gè)品種,該品種屬于生育期偏長(zhǎng),主莖節(jié)數(shù)少,底莢高度高,莢數(shù)粒數(shù)多,籽粒小的低產(chǎn)型品種。

    從不同角度進(jìn)行的聚類分析結(jié)果可以看出,綜合性狀的聚類與兩種性狀組合聚類的結(jié)果存在明顯不同。少數(shù)性狀組合的聚類只考慮到少數(shù)性狀表現(xiàn)的相似性,各類群內(nèi)只有被聚類的性狀表現(xiàn)一致,而大多性狀表現(xiàn)明顯差異。綜合性狀的聚類使類群內(nèi)大多性狀表現(xiàn)一致,僅個(gè)別性狀有差異,結(jié)果更具整體性和合理性[3]。

    3 結(jié)論和討論

    本研究結(jié)果表明,供試的17份南方春大豆品種秋播下的生物學(xué)性狀變異系數(shù)以底莢高度最大、生育期最小。單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)及單株莢數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān)。主成份分析結(jié)果表明,秋播條件下不同大豆品種的生物學(xué)性狀可分為產(chǎn)量關(guān)鍵因子、產(chǎn)量影響因子、生長(zhǎng)勢(shì)因子和生育期與株高綜合因子4個(gè)主成分,它們基本上代表了所有性狀的變異信息。系統(tǒng)聚類分析結(jié)果表明,17份品種類群較為復(fù)雜,其中有莢少粒少、粒大的低產(chǎn)型品種,粒少、粒大的中高產(chǎn)型品種,莢多粒多,籽小的高產(chǎn)型品種等;同時(shí)聚類分析結(jié)果還表明了來自同一單位育成的品種不一定類型相似??傊?,主成分分析結(jié)果反映了17份春大豆秋播下生物學(xué)性狀的主導(dǎo)影響因子,聚類結(jié)果反映了產(chǎn)量的差異,這為正確評(píng)價(jià)大豆品種的適應(yīng)性及科學(xué)地選擇品種提供了依據(jù)[6]。

    參考文獻(xiàn)

    [1]殷冬梅,李銓柱,崔黨群.花生主要農(nóng)藝性狀的相關(guān)性及聚類分析[J].中國(guó)油料作物學(xué)報(bào),2010,32(2):212-216.

    [2]李清華,黃金堂,陳海玲,等.27份花生種質(zhì)資源的主成份分析與遺傳距離測(cè)定[J].植物遺傳資源學(xué)報(bào),2011,12(4):519-524.

    [3]李淑華,荊紹凌,楊丹萍,等.聚類分析法在玉米種質(zhì)資源中的應(yīng)用[J].農(nóng)業(yè)技術(shù),2008,28(5):43-45.

    [4]馮永祥,劉沐江,汪秀志,等.寒地水稻產(chǎn)量性狀相關(guān)及聚類分析[J].安徽農(nóng)業(yè)科學(xué),2008,36(23):9 925-9 928.

    [5]蘇春華,李育明,黃迎冬,等.甘薯主要親本材料的主成分分析及聚類分析[J].雜糧作物,2007,27(6):405-409.

    [6]張玉革,胡緒彬.基于主成分和聚類分析的大豆品種生物學(xué)性狀的比較研究[J].大豆科學(xué),2004,23(3):178-183[J].

    [7]李向華,常釹鎮(zhèn).中國(guó)春大豆品種聚類分析及主成分分析[J].作物學(xué)報(bào),1998,24(3):325-332.

    [8]劉燦洪,林榮輝,劉德金.福建省大豆地方品種多元統(tǒng)計(jì)分析[J].福建農(nóng)業(yè)學(xué)報(bào),1999,14(4):11-15.

    [9] 唐啟義,馮光明. 實(shí)用統(tǒng)計(jì)分析及DPS其數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [10]王林輝.基于主成分分析的棉花品種綜合評(píng)價(jià)及聚類分析[J].廣東農(nóng)業(yè)科學(xué),2009,(1):29-32. (責(zé)編:張宏民)endprint

    2.2 生物學(xué)性狀的相關(guān)性分析 對(duì)17分春大豆品種主要農(nóng)藝性狀的相關(guān)分析表明(表2):?jiǎn)沃炅?shù)與單株莢數(shù)呈極顯著正相關(guān),百粒重與單株莢數(shù)、單株粒數(shù)呈極顯著負(fù)相關(guān);單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)、單株莢數(shù)及單株粒數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān);主莖高與單株粒數(shù)呈負(fù)相關(guān),與百粒重呈正相關(guān);主莖節(jié)數(shù)與底莢高度呈負(fù)相關(guān)、與單株粒數(shù)呈正相關(guān);有效分枝與底莢高度呈正相關(guān)。

    2.3 主成分分析 對(duì)供試的17份春大豆品種的9個(gè)主要觀測(cè)值采用DPS數(shù)據(jù)處理軟件進(jìn)行主成分分析,計(jì)算得特征值Yi、特征值累計(jì)百分率和主成份特征向量。各特征值大小代表各綜合指標(biāo)遺傳方差大小,各特征值累計(jì)百分率代表各有關(guān)綜合指標(biāo)對(duì)總遺傳方差貢獻(xiàn)的百分率。為提高分析精度,只選用對(duì)總遺傳方差貢獻(xiàn)率達(dá)85.43%的前4個(gè)特征值(表3)。

    從表3可以看出,第一主成分的特征向量中,單株莢數(shù)、單株粒數(shù)和單株粒重的權(quán)重系數(shù)為正值,百粒重的權(quán)重系數(shù)為負(fù)值,均較大,構(gòu)成了變異主要來源,說明單株莢數(shù)、單株粒數(shù)與單株粒重呈同步增長(zhǎng)態(tài)勢(shì),而單株莢數(shù)、單株粒數(shù)同時(shí)影響了籽粒的大小,故可稱為產(chǎn)量關(guān)鍵因子;第二主成分的特征向量中,有效分枝、底莢高度權(quán)重系數(shù)為正值,百粒重、單株粒重權(quán)重系數(shù)為負(fù)值,且較大,構(gòu)成了變異主要來源,說明有效分枝過多、底莢高度過高將影響大豆籽粒大小,進(jìn)而影響了大豆的單株產(chǎn)量,故可稱為產(chǎn)量影響因子;第三主成分的特征向量中,主莖高、主莖節(jié)數(shù)和有效分枝權(quán)重系數(shù)均為正數(shù),且較大,這些主要反映了植株生長(zhǎng)勢(shì),故稱為生長(zhǎng)勢(shì)因子;第四主成份的特征向量中,生育期與主莖高、底莢高度權(quán)重系數(shù)較大,可稱之為生育期與植株高度綜合因子。

    2.4 聚類分析

    2.4.1 不同性狀組合的聚類分析 對(duì)參試品種分別按植株性狀與經(jīng)濟(jì)性狀2種性狀組合分別進(jìn)行少類性狀聚類。其中植株性狀包括:主莖高、主莖節(jié)數(shù)、有效分枝與底莢高度;經(jīng)濟(jì)性狀包括:?jiǎn)沃昵v數(shù)、單株粒數(shù)、百粒重與單株粒重。聚類結(jié)果表明:參與聚類的性狀組合不同,其聚類結(jié)果各異。按植株性狀方面聚類(圖1),17個(gè)品種在聚類水平D2=1.07時(shí)可分成12類,類型較為多樣復(fù)雜;按經(jīng)濟(jì)性狀方面聚類(圖2),17個(gè)品種在聚類水平D2=4.06時(shí)只分成了四大類,且從圖中可以看出福建省育成的多數(shù)春大豆品種在經(jīng)濟(jì)性狀方面較為相似。

    2.4.2 綜合性狀聚類分析 對(duì)參試品種9個(gè)觀測(cè)值進(jìn)行標(biāo)準(zhǔn)化變換后在歐氏距離水平上采用離差平方和法進(jìn)行系統(tǒng)聚類分析(圖3),根據(jù)統(tǒng)一類群內(nèi)大豆品種類群間距離接近且綜合性狀差異值較小的原則,17份大豆品種在聚類水平D2=4.24時(shí)可以劃分為7大類群。

    第1類群有1、2和7號(hào)3個(gè)品種,該類群品種有效分枝、底莢高度兩性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于單株莢數(shù)、粒數(shù)偏少,籽粒較大的低產(chǎn)型品種;第2類群有4、5、6和17號(hào)4個(gè)品種,該類群品種在植株性狀和單株莢數(shù)存在明顯差異,其他性狀均表現(xiàn)相近,屬于生育期較長(zhǎng),粒少、粒大的中產(chǎn)型品種;第3類群有3和9號(hào)2個(gè)品種,兩品種在產(chǎn)量性狀表現(xiàn)上存在較小差異,均屬于生育期適中,植株旺盛,結(jié)莢數(shù)適宜,粒少、粒大的中高產(chǎn)型品種;第4類群有8、16、11和15號(hào)4個(gè)品種,該類群生育期、主莖節(jié)數(shù)和有效分枝性狀表現(xiàn)上存在較大差異外,其他性狀均表現(xiàn)相近,屬于底莢高度較低,莢多粒多,籽粒較小的高產(chǎn)型品種;第5類群有10號(hào)1個(gè)品種,該品種屬于株高適中,分枝旺盛,底莢高度較高,粒多、粒小的早熟低產(chǎn)型品種;第6類群12和14號(hào)3個(gè)品種,該類群品種屬于生育期較長(zhǎng),底莢高度偏高,籽粒大小中等,其他性狀均表現(xiàn)優(yōu)異的中產(chǎn)型品種。第7類群有13號(hào)1個(gè)品種,該品種屬于生育期偏長(zhǎng),主莖節(jié)數(shù)少,底莢高度高,莢數(shù)粒數(shù)多,籽粒小的低產(chǎn)型品種。

    從不同角度進(jìn)行的聚類分析結(jié)果可以看出,綜合性狀的聚類與兩種性狀組合聚類的結(jié)果存在明顯不同。少數(shù)性狀組合的聚類只考慮到少數(shù)性狀表現(xiàn)的相似性,各類群內(nèi)只有被聚類的性狀表現(xiàn)一致,而大多性狀表現(xiàn)明顯差異。綜合性狀的聚類使類群內(nèi)大多性狀表現(xiàn)一致,僅個(gè)別性狀有差異,結(jié)果更具整體性和合理性[3]。

    3 結(jié)論和討論

    本研究結(jié)果表明,供試的17份南方春大豆品種秋播下的生物學(xué)性狀變異系數(shù)以底莢高度最大、生育期最小。單株粒重與單株粒數(shù)呈顯著正相關(guān),與生育期、主莖節(jié)數(shù)及單株莢數(shù)呈正相關(guān),與主莖高、有效分枝、底莢高度、百粒重呈負(fù)相關(guān)。主成份分析結(jié)果表明,秋播條件下不同大豆品種的生物學(xué)性狀可分為產(chǎn)量關(guān)鍵因子、產(chǎn)量影響因子、生長(zhǎng)勢(shì)因子和生育期與株高綜合因子4個(gè)主成分,它們基本上代表了所有性狀的變異信息。系統(tǒng)聚類分析結(jié)果表明,17份品種類群較為復(fù)雜,其中有莢少粒少、粒大的低產(chǎn)型品種,粒少、粒大的中高產(chǎn)型品種,莢多粒多,籽小的高產(chǎn)型品種等;同時(shí)聚類分析結(jié)果還表明了來自同一單位育成的品種不一定類型相似。總之,主成分分析結(jié)果反映了17份春大豆秋播下生物學(xué)性狀的主導(dǎo)影響因子,聚類結(jié)果反映了產(chǎn)量的差異,這為正確評(píng)價(jià)大豆品種的適應(yīng)性及科學(xué)地選擇品種提供了依據(jù)[6]。

    參考文獻(xiàn)

    [1]殷冬梅,李銓柱,崔黨群.花生主要農(nóng)藝性狀的相關(guān)性及聚類分析[J].中國(guó)油料作物學(xué)報(bào),2010,32(2):212-216.

    [2]李清華,黃金堂,陳海玲,等.27份花生種質(zhì)資源的主成份分析與遺傳距離測(cè)定[J].植物遺傳資源學(xué)報(bào),2011,12(4):519-524.

    [3]李淑華,荊紹凌,楊丹萍,等.聚類分析法在玉米種質(zhì)資源中的應(yīng)用[J].農(nóng)業(yè)技術(shù),2008,28(5):43-45.

    [4]馮永祥,劉沐江,汪秀志,等.寒地水稻產(chǎn)量性狀相關(guān)及聚類分析[J].安徽農(nóng)業(yè)科學(xué),2008,36(23):9 925-9 928.

    [5]蘇春華,李育明,黃迎冬,等.甘薯主要親本材料的主成分分析及聚類分析[J].雜糧作物,2007,27(6):405-409.

    [6]張玉革,胡緒彬.基于主成分和聚類分析的大豆品種生物學(xué)性狀的比較研究[J].大豆科學(xué),2004,23(3):178-183[J].

    [7]李向華,常釹鎮(zhèn).中國(guó)春大豆品種聚類分析及主成分分析[J].作物學(xué)報(bào),1998,24(3):325-332.

    [8]劉燦洪,林榮輝,劉德金.福建省大豆地方品種多元統(tǒng)計(jì)分析[J].福建農(nóng)業(yè)學(xué)報(bào),1999,14(4):11-15.

    [9] 唐啟義,馮光明. 實(shí)用統(tǒng)計(jì)分析及DPS其數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2002.

    [10]王林輝.基于主成分分析的棉花品種綜合評(píng)價(jià)及聚類分析[J].廣東農(nóng)業(yè)科學(xué),2009,(1):29-32. (責(zé)編:張宏民)endprint

    猜你喜歡
    秋播聚類分析
    “冷季型草+秋播”生草模式總結(jié)
    北欒秋播育苗技術(shù)
    秋播大蒜優(yōu)質(zhì)高產(chǎn)栽培技術(shù)
    新媒體用戶行為模式分析
    農(nóng)村居民家庭人均生活消費(fèi)支出分析
    棗陽 黨組織服務(wù)秋播生產(chǎn)
    基于省會(huì)城市經(jīng)濟(jì)發(fā)展程度的實(shí)證分析
    基于聚類分析的互聯(lián)網(wǎng)廣告投放研究
    “縣級(jí)供電企業(yè)生產(chǎn)經(jīng)營(yíng)統(tǒng)計(jì)一套”表輔助決策模式研究
    拔節(jié)期追氮對(duì)春播和秋播糯玉米淀粉膠凝和回生特性的影響
    欧美+亚洲+日韩+国产| 丰满饥渴人妻一区二区三| 午夜两性在线视频| 精品视频人人做人人爽| 午夜福利在线免费观看网站| 国产熟女午夜一区二区三区| 我的亚洲天堂| 国产一区二区三区视频了| 麻豆乱淫一区二区| 热re99久久国产66热| 久久av网站| 热re99久久精品国产66热6| 亚洲午夜精品一区,二区,三区| 色综合婷婷激情| av片东京热男人的天堂| 中亚洲国语对白在线视频| 亚洲精品中文字幕一二三四区 | 十八禁人妻一区二区| 一区二区三区精品91| 美女主播在线视频| 啦啦啦视频在线资源免费观看| 亚洲精品一卡2卡三卡4卡5卡| 丝袜喷水一区| 国产精品av久久久久免费| 亚洲熟女精品中文字幕| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| 亚洲精品一二三| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 成人永久免费在线观看视频 | 狂野欧美激情性xxxx| 亚洲人成伊人成综合网2020| 性少妇av在线| 伦理电影免费视频| 日本av手机在线免费观看| 黄色怎么调成土黄色| 高潮久久久久久久久久久不卡| 亚洲欧美日韩高清在线视频 | 十八禁网站免费在线| 国产精品久久久人人做人人爽| 精品人妻在线不人妻| 国产主播在线观看一区二区| 美女高潮喷水抽搐中文字幕| 老汉色av国产亚洲站长工具| 极品教师在线免费播放| 亚洲精华国产精华精| 啪啪无遮挡十八禁网站| 欧美日韩成人在线一区二区| 久久久国产成人免费| 国产熟女午夜一区二区三区| 可以免费在线观看a视频的电影网站| 伦理电影免费视频| 亚洲精品中文字幕在线视频| 亚洲男人天堂网一区| 夜夜爽天天搞| 伦理电影免费视频| 超碰97精品在线观看| 一进一出好大好爽视频| 国产精品九九99| 巨乳人妻的诱惑在线观看| 十八禁人妻一区二区| 国产福利在线免费观看视频| 国产欧美亚洲国产| kizo精华| 午夜福利一区二区在线看| 国产免费av片在线观看野外av| 日本精品一区二区三区蜜桃| 亚洲精华国产精华精| 久久久国产一区二区| 又黄又粗又硬又大视频| 99re在线观看精品视频| 欧美精品人与动牲交sv欧美| 亚洲精品国产一区二区精华液| 久久精品亚洲精品国产色婷小说| 建设人人有责人人尽责人人享有的| 99久久99久久久精品蜜桃| 人妻 亚洲 视频| 好男人电影高清在线观看| 黄色成人免费大全| 黄色毛片三级朝国网站| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区视频了| 日本黄色视频三级网站网址 | 黄色毛片三级朝国网站| 亚洲欧美一区二区三区久久| 精品亚洲成国产av| 精品久久久久久久毛片微露脸| bbb黄色大片| 99国产精品99久久久久| 午夜福利视频在线观看免费| 久久亚洲真实| 国产在线精品亚洲第一网站| 国产精品久久久久久精品古装| 日韩欧美三级三区| 别揉我奶头~嗯~啊~动态视频| 精品国产乱子伦一区二区三区| 久久久久久免费高清国产稀缺| 亚洲成人国产一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看舔阴道视频| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 一边摸一边做爽爽视频免费| 热99久久久久精品小说推荐| 久久性视频一级片| 窝窝影院91人妻| 动漫黄色视频在线观看| 性高湖久久久久久久久免费观看| 久久久久久亚洲精品国产蜜桃av| 久久久久网色| 欧美午夜高清在线| 高清毛片免费观看视频网站 | 精品免费久久久久久久清纯 | 又黄又粗又硬又大视频| 久久精品aⅴ一区二区三区四区| 欧美精品一区二区大全| 久久午夜综合久久蜜桃| 激情视频va一区二区三区| 久久中文看片网| 亚洲人成电影观看| 午夜免费鲁丝| 久久精品国产a三级三级三级| 午夜福利在线免费观看网站| 欧美精品一区二区免费开放| 一级毛片电影观看| 成人特级黄色片久久久久久久 | 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 久久狼人影院| 欧美乱妇无乱码| 丁香六月欧美| 午夜福利欧美成人| 91国产中文字幕| 日本av手机在线免费观看| 好男人电影高清在线观看| 99精国产麻豆久久婷婷| 久久久国产成人免费| 高清av免费在线| 久久久久久久久免费视频了| 精品久久久精品久久久| 一区二区三区激情视频| 精品亚洲乱码少妇综合久久| 日韩欧美一区视频在线观看| videosex国产| 男女边摸边吃奶| 国产极品粉嫩免费观看在线| 久久精品亚洲熟妇少妇任你| 日韩欧美三级三区| 夜夜夜夜夜久久久久| 757午夜福利合集在线观看| 考比视频在线观看| 午夜老司机福利片| 国产精品99久久99久久久不卡| 久久毛片免费看一区二区三区| 成人手机av| 久久 成人 亚洲| 女人高潮潮喷娇喘18禁视频| 久久久久国内视频| 欧美另类亚洲清纯唯美| 午夜福利影视在线免费观看| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 久久久国产精品麻豆| 国产精品美女特级片免费视频播放器 | 丝瓜视频免费看黄片| 免费高清在线观看日韩| 久久ye,这里只有精品| 亚洲国产成人一精品久久久| 69av精品久久久久久 | 午夜福利免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久精品古装| 露出奶头的视频| 777米奇影视久久| 午夜两性在线视频| 香蕉丝袜av| 亚洲国产中文字幕在线视频| 久久人人97超碰香蕉20202| 国产有黄有色有爽视频| 亚洲国产中文字幕在线视频| 久久久久久久国产电影| 性色av乱码一区二区三区2| 多毛熟女@视频| 超碰97精品在线观看| 久久精品人人爽人人爽视色| 国产精品欧美亚洲77777| 另类精品久久| 久久99一区二区三区| 91大片在线观看| 大型黄色视频在线免费观看| 一级片免费观看大全| 精品国产乱子伦一区二区三区| 日本五十路高清| 久久亚洲真实| 在线观看舔阴道视频| 两个人免费观看高清视频| 成人黄色视频免费在线看| 国产亚洲精品一区二区www | 成人手机av| 亚洲午夜理论影院| 亚洲视频免费观看视频| 久久九九热精品免费| svipshipincom国产片| 久久中文看片网| 久久久久精品人妻al黑| 亚洲av日韩在线播放| 免费观看a级毛片全部| 精品午夜福利视频在线观看一区 | 午夜激情av网站| 巨乳人妻的诱惑在线观看| 久久国产精品大桥未久av| 一本色道久久久久久精品综合| 一级黄色大片毛片| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 美女视频免费永久观看网站| 电影成人av| 国产1区2区3区精品| 少妇精品久久久久久久| 超色免费av| 午夜日韩欧美国产| 99re6热这里在线精品视频| 操出白浆在线播放| 侵犯人妻中文字幕一二三四区| 成人亚洲精品一区在线观看| 国产精品一区二区在线不卡| 美女午夜性视频免费| 久久午夜亚洲精品久久| 水蜜桃什么品种好| 免费观看av网站的网址| 午夜激情av网站| 亚洲,欧美精品.| 日本撒尿小便嘘嘘汇集6| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| √禁漫天堂资源中文www| 嫩草影视91久久| 啪啪无遮挡十八禁网站| 咕卡用的链子| 亚洲男人天堂网一区| 国产精品亚洲av一区麻豆| 午夜两性在线视频| 久久人妻熟女aⅴ| 高潮久久久久久久久久久不卡| 国产精品久久久av美女十八| 国产极品粉嫩免费观看在线| 在线观看www视频免费| 久久人妻熟女aⅴ| 精品少妇久久久久久888优播| tube8黄色片| 成人手机av| 欧美性长视频在线观看| 人人妻人人澡人人爽人人夜夜| 欧美国产精品一级二级三级| 国产精品麻豆人妻色哟哟久久| 成年人黄色毛片网站| 免费日韩欧美在线观看| 亚洲av电影在线进入| 最新美女视频免费是黄的| 91成人精品电影| 欧美精品一区二区免费开放| 侵犯人妻中文字幕一二三四区| 黄色视频,在线免费观看| 久久久久久久久免费视频了| 国产不卡一卡二| 69精品国产乱码久久久| 变态另类成人亚洲欧美熟女 | 国产精品熟女久久久久浪| 国产欧美亚洲国产| 伦理电影免费视频| 久久精品熟女亚洲av麻豆精品| 国产成人av教育| 黄色成人免费大全| 大陆偷拍与自拍| 淫妇啪啪啪对白视频| 国产日韩欧美在线精品| 免费不卡黄色视频| 国产高清视频在线播放一区| 亚洲精华国产精华精| 国产亚洲精品一区二区www | 精品国产一区二区久久| 婷婷成人精品国产| 欧美日韩成人在线一区二区| 一级a爱视频在线免费观看| 国产精品自产拍在线观看55亚洲 | 精品一区二区三区视频在线观看免费 | 五月天丁香电影| 岛国毛片在线播放| 人人澡人人妻人| 国产不卡一卡二| 亚洲人成伊人成综合网2020| 18禁观看日本| 十八禁人妻一区二区| 在线十欧美十亚洲十日本专区| 国产精品九九99| 最新美女视频免费是黄的| 国产不卡一卡二| 性少妇av在线| 日韩欧美国产一区二区入口| 国产欧美日韩综合在线一区二区| 丰满饥渴人妻一区二区三| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影 | 超碰成人久久| 精品国产一区二区三区久久久樱花| 淫妇啪啪啪对白视频| av视频免费观看在线观看| 18禁裸乳无遮挡动漫免费视频| 飞空精品影院首页| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 国产免费现黄频在线看| 日日夜夜操网爽| 国产欧美日韩一区二区三| 一本—道久久a久久精品蜜桃钙片| 12—13女人毛片做爰片一| 最新的欧美精品一区二区| 成年人午夜在线观看视频| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 欧美日韩中文字幕国产精品一区二区三区 | av不卡在线播放| 日韩熟女老妇一区二区性免费视频| 国产一区二区 视频在线| 精品一品国产午夜福利视频| 人妻久久中文字幕网| 国产在线一区二区三区精| 制服诱惑二区| 亚洲精品中文字幕在线视频| 制服诱惑二区| 久久婷婷成人综合色麻豆| 又大又爽又粗| 国产精品亚洲一级av第二区| 18在线观看网站| 黄色视频在线播放观看不卡| 电影成人av| 色婷婷久久久亚洲欧美| 亚洲精品中文字幕一二三四区 | 中亚洲国语对白在线视频| 免费在线观看视频国产中文字幕亚洲| 肉色欧美久久久久久久蜜桃| 老司机深夜福利视频在线观看| 久久久国产精品麻豆| 国产99久久九九免费精品| 久久中文看片网| 高清在线国产一区| 蜜桃在线观看..| 国产精品欧美亚洲77777| 日韩视频在线欧美| av网站在线播放免费| 丁香欧美五月| 婷婷丁香在线五月| 亚洲第一av免费看| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦视频在线资源免费观看| 亚洲少妇的诱惑av| 露出奶头的视频| 女人精品久久久久毛片| a级毛片黄视频| 18禁国产床啪视频网站| 精品熟女少妇八av免费久了| 天天躁日日躁夜夜躁夜夜| 精品一品国产午夜福利视频| 嫁个100分男人电影在线观看| 成人精品一区二区免费| 亚洲美女黄片视频| av超薄肉色丝袜交足视频| 看免费av毛片| 在线观看人妻少妇| 汤姆久久久久久久影院中文字幕| 精品国产乱子伦一区二区三区| 一个人免费看片子| 免费人妻精品一区二区三区视频| 午夜日韩欧美国产| 少妇被粗大的猛进出69影院| av福利片在线| 老熟妇仑乱视频hdxx| 欧美人与性动交α欧美精品济南到| 一区二区三区激情视频| 制服诱惑二区| 国产欧美日韩一区二区三区在线| 国产成人av教育| 亚洲 国产 在线| bbb黄色大片| 好男人电影高清在线观看| 精品国产乱码久久久久久小说| 久久久久久人人人人人| 人人妻人人澡人人看| 18禁观看日本| 国产视频一区二区在线看| a级毛片在线看网站| 国产精品98久久久久久宅男小说| 18在线观看网站| 又紧又爽又黄一区二区| 中文欧美无线码| 亚洲精品国产区一区二| 久9热在线精品视频| 丝袜在线中文字幕| 最黄视频免费看| 午夜精品国产一区二区电影| 国产麻豆69| 国产在线一区二区三区精| 99热网站在线观看| 巨乳人妻的诱惑在线观看| 人妻久久中文字幕网| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产一区二区精华液| 在线观看免费视频网站a站| 欧美国产精品一级二级三级| 久久久欧美国产精品| 国产aⅴ精品一区二区三区波| 欧美激情 高清一区二区三区| 久久久久精品国产欧美久久久| 日韩制服丝袜自拍偷拍| 最新的欧美精品一区二区| 黄网站色视频无遮挡免费观看| 亚洲欧美一区二区三区黑人| 91精品国产国语对白视频| 久久久水蜜桃国产精品网| 免费一级毛片在线播放高清视频 | 久久精品亚洲精品国产色婷小说| 99国产精品一区二区蜜桃av | 久久精品国产亚洲av香蕉五月 | 麻豆乱淫一区二区| 久久久国产成人免费| 亚洲全国av大片| 国产一区二区在线观看av| videos熟女内射| 大码成人一级视频| 欧美黑人精品巨大| 久久人人爽av亚洲精品天堂| av有码第一页| 又黄又粗又硬又大视频| 电影成人av| 国产日韩欧美视频二区| 免费观看a级毛片全部| 久久人人97超碰香蕉20202| av片东京热男人的天堂| 成年动漫av网址| 久久久精品94久久精品| 日韩欧美三级三区| 99国产精品99久久久久| 天堂中文最新版在线下载| 黄色毛片三级朝国网站| 男女床上黄色一级片免费看| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清 | 老司机午夜福利在线观看视频 | 国产成人免费观看mmmm| 精品视频人人做人人爽| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 一级毛片精品| 免费观看av网站的网址| 又紧又爽又黄一区二区| 美女主播在线视频| 在线观看人妻少妇| 丰满人妻熟妇乱又伦精品不卡| 午夜福利,免费看| 久久 成人 亚洲| 99国产综合亚洲精品| 欧美国产精品va在线观看不卡| 一级片免费观看大全| 热re99久久精品国产66热6| 女警被强在线播放| 搡老岳熟女国产| 免费久久久久久久精品成人欧美视频| 午夜福利,免费看| 91精品国产国语对白视频| 少妇裸体淫交视频免费看高清 | 婷婷丁香在线五月| a级片在线免费高清观看视频| 亚洲精品久久午夜乱码| 丝袜美腿诱惑在线| 久热爱精品视频在线9| 亚洲熟女精品中文字幕| 新久久久久国产一级毛片| 美女扒开内裤让男人捅视频| 在线永久观看黄色视频| 9色porny在线观看| 黑丝袜美女国产一区| 一边摸一边抽搐一进一出视频| 一夜夜www| 国产亚洲精品一区二区www | 亚洲情色 制服丝袜| 精品亚洲成国产av| 久久九九热精品免费| 成人18禁在线播放| 国产xxxxx性猛交| 日韩一卡2卡3卡4卡2021年| 亚洲精品中文字幕在线视频| 麻豆av在线久日| 老司机亚洲免费影院| 欧美日韩亚洲综合一区二区三区_| 免费一级毛片在线播放高清视频 | 啪啪无遮挡十八禁网站| 欧美性长视频在线观看| 国内毛片毛片毛片毛片毛片| 色播在线永久视频| 五月开心婷婷网| 免费一级毛片在线播放高清视频 | 1024视频免费在线观看| 老司机在亚洲福利影院| 亚洲七黄色美女视频| 免费在线观看视频国产中文字幕亚洲| 黑丝袜美女国产一区| 色播在线永久视频| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 亚洲黑人精品在线| 久久久久久人人人人人| 午夜激情久久久久久久| 99精品欧美一区二区三区四区| 91精品国产国语对白视频| 精品亚洲成国产av| 国产区一区二久久| 12—13女人毛片做爰片一| 黑人巨大精品欧美一区二区蜜桃| 高清黄色对白视频在线免费看| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区久久久樱花| 精品高清国产在线一区| 国产精品九九99| 国产精品电影一区二区三区 | 悠悠久久av| 男人操女人黄网站| 久久国产精品影院| 欧美在线一区亚洲| 亚洲色图av天堂| 久久久久久久久久久久大奶| 在线 av 中文字幕| 51午夜福利影视在线观看| 欧美人与性动交α欧美软件| 一夜夜www| 女性被躁到高潮视频| 69av精品久久久久久 | 大型av网站在线播放| 美女国产高潮福利片在线看| av网站免费在线观看视频| 欧美精品亚洲一区二区| 国产精品久久久久久人妻精品电影 | 午夜免费成人在线视频| 国产麻豆69| 亚洲美女黄片视频| 亚洲成国产人片在线观看| 色在线成人网| 日韩中文字幕视频在线看片| 国产精品久久久人人做人人爽| 少妇猛男粗大的猛烈进出视频| 人成视频在线观看免费观看| 天天躁夜夜躁狠狠躁躁| 色尼玛亚洲综合影院| 色综合欧美亚洲国产小说| 国产欧美日韩精品亚洲av| 欧美精品啪啪一区二区三区| 欧美精品人与动牲交sv欧美| 麻豆乱淫一区二区| 成人黄色视频免费在线看| 在线av久久热| 黑人巨大精品欧美一区二区mp4| 天堂中文最新版在线下载| 欧美激情久久久久久爽电影 | 久久久水蜜桃国产精品网| 精品卡一卡二卡四卡免费| 叶爱在线成人免费视频播放| 欧美性长视频在线观看| 亚洲国产欧美一区二区综合| 成人特级黄色片久久久久久久 | 777米奇影视久久| 国产日韩欧美在线精品| 国产成人欧美在线观看 | 成年动漫av网址| 亚洲精品国产色婷婷电影| 一本一本久久a久久精品综合妖精| 精品国内亚洲2022精品成人 | 国产成人系列免费观看| 自拍欧美九色日韩亚洲蝌蚪91| av欧美777| 国产人伦9x9x在线观看| 99国产精品免费福利视频| 国产成人免费无遮挡视频| 成在线人永久免费视频| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 国产精品免费视频内射| 青草久久国产| 亚洲 国产 在线| 18禁裸乳无遮挡动漫免费视频| 亚洲第一av免费看| 欧美精品av麻豆av| 成年人免费黄色播放视频| 大片免费播放器 马上看| 大型av网站在线播放| 好男人电影高清在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲九九香蕉| 国产成人欧美| 亚洲欧洲日产国产| 露出奶头的视频| 成人特级黄色片久久久久久久 | av线在线观看网站| 9色porny在线观看| 亚洲一码二码三码区别大吗| 大型黄色视频在线免费观看| 国产高清国产精品国产三级| 国产有黄有色有爽视频| 老司机靠b影院| 成人18禁在线播放| 日韩欧美免费精品| 欧美亚洲日本最大视频资源| 高清av免费在线| 免费高清在线观看日韩| 热re99久久精品国产66热6| 国产男女内射视频|