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

    浙江省糯玉米品種穩(wěn)定性、適應(yīng)性和試點(diǎn)綜合評(píng)價(jià)

    2017-07-31 21:09:59李潞潞徐海明
    關(guān)鍵詞:糯玉米區(qū)域試驗(yàn)性狀

    李 燕,林 峰,李潞潞,徐海明*

    (1.浙江省種子管理總站,杭州310020;2.浙江大學(xué)農(nóng)業(yè)與生物技術(shù)學(xué)院作物科學(xué)研究所,杭州310058)

    浙江省糯玉米品種穩(wěn)定性、適應(yīng)性和試點(diǎn)綜合評(píng)價(jià)

    李 燕1,林 峰2,李潞潞2,徐海明2*1

    (1.浙江省種子管理總站,杭州310020;2.浙江大學(xué)農(nóng)業(yè)與生物技術(shù)學(xué)院作物科學(xué)研究所,杭州310058)

    為明確品種、地點(diǎn)、年份因子對(duì)浙江省糯玉米產(chǎn)量性狀和構(gòu)成產(chǎn)量相關(guān)性狀的貢獻(xiàn)率,本研究利用基于混合線性模型的統(tǒng)計(jì)分析方法,對(duì)2013—2015年浙江省糯玉米區(qū)試數(shù)據(jù)進(jìn)行分析,從試點(diǎn)的生產(chǎn)水平、精度、分辨力、代表性和穩(wěn)定性等5方面對(duì)浙江省糯玉米8個(gè)區(qū)試點(diǎn)進(jìn)行初步的綜合評(píng)價(jià)研究,同時(shí)對(duì)29個(gè)糯玉米參試品種進(jìn)行豐產(chǎn)性、穩(wěn)定性和適應(yīng)性評(píng)價(jià)。結(jié)果表明:8個(gè)試點(diǎn)的精確度較高,試驗(yàn)均有效;8個(gè)試點(diǎn)的生產(chǎn)水平有高有低,代表著不同的地點(diǎn)生產(chǎn)水平;嘉善、杭州、淳安等3個(gè)試點(diǎn)對(duì)糯玉米產(chǎn)量性狀的分辨力較高;東陽(yáng)、江山、淳安、杭州等4個(gè)試點(diǎn)對(duì)糯玉米產(chǎn)量性狀具有較高的普遍代表性,宜鑒別廣泛適應(yīng)性的品種;東陽(yáng)、嘉善2個(gè)試點(diǎn)的適宜性較好,穩(wěn)定性較高;品種、地點(diǎn)、品種與地點(diǎn)互作效應(yīng)、年份與地點(diǎn)互作效應(yīng)對(duì)糯玉米產(chǎn)量性狀有顯著影響,品種和地點(diǎn)互作效應(yīng)對(duì)糯玉米產(chǎn)量構(gòu)成的農(nóng)藝性狀沒(méi)有顯著影響。品種產(chǎn)量性狀穩(wěn)定性和豐產(chǎn)性分析結(jié)果表明,糯玉米品種“科甜糯2號(hào)”的產(chǎn)量和穩(wěn)定性都較好。

    糯玉米;農(nóng)藝性狀;區(qū)域試驗(yàn);試點(diǎn)評(píng)價(jià);品種評(píng)價(jià);混合線性模型

    Summary Waxy maize is an importanteconomic crop in Zhejiang Province.Agronomic traits,especially the yield and yieldrelated traits ofwaxy maize are the focus of attention paid by both waxy maize breeders and producers.Most agronomic traits of waxy maize are quantitative traits affected by environment,variety-by-location interaction,variety and location-by-year interaction.It is very important to understand the contribution rates of each factor to agronomic traits of waxy maize,so as to screen appropriate varieties for Zhejiang Province.

    In this study,a regionalvariety testwith random complete block design was implemented in multiple locations from 2013 to 2015.Totally,29 waxy maize varieties and eight test locations were arranged with three blocks in each location.The yield and seven yield-related traits,ear length,ear diameter,net ear rate,rows per ear,row grains,fresh 1 000-grain mass,and seed producing rate were investigated.Based on the location effects,variety effects,year effects,location-by-year interaction effects, variety-by-location interaction effects,predicted by a mixed linear model,the stability and adaptability of the 29 varieties were analyzed,and the eight locations were comprehensively evaluated using productive level,precision,distinguishing ability, representation and stability.SAS and QTModelsoftwares were applied forthe analysis.

    The results of variance analysis showed that the location effects,variety effects,year effects,year-by-location interaction effects and variety-by-location interaction effects were significantforyield trait,and the variety-by-location interaction effects for seven yield-related traits were not significant.Wherein,67.6%of variance for rows per ear was from variety effects,36.3%of variance for yield from location effects.Additionally,the location related effects explained 6.9%-62.0%ofthe totalvariance for each agronomic trait.In accordance with the contribution rate of location?related effects,the traits from high to low contributions were the yield,seed producing rate,row grains,fresh 1 000-grain mass,net ear rate,ear diameter,ear length, rows per ear,respectively.The results for comprehensive evaluation on location showed that productive levels in the eight locations were diverse;the trial precisions in the eight locations were high,showing that the trials were all valid.In these test locations,Jiashan,Hangzhou and Chun’an possessed high distinguishability;Dongyang,Jiangshan,Chun’an and Hangzhou possessed high representativeness;Dongyang and Jiashan possessed high stability and productive level.

    On evaluation of variety stability and yield-ability,according to the results of multiple comparisons by Tukey method for yield trait,the varieties marked by“a”were regarded as appropriate ones,and among them,those with yield mean value equalto or larger than the average ofallappropriate ones were determined as the mostappropriate varieties.Based on the variety effects and variety-by-location interaction effects,the stability of the variety was estimated.As a result,“Ketiannuo 2”was the most appropriate variety for Zhejiang Province.

    作物品種區(qū)域試驗(yàn)是通過(guò)多環(huán)境(多地點(diǎn)或多年多點(diǎn))的試驗(yàn)來(lái)鑒定參試品種的豐產(chǎn)性、穩(wěn)定性和適應(yīng)性,以及其他綜合特性,以確定參試品種推廣應(yīng)用價(jià)值和適宜種植區(qū)域,為品種審定和推廣提供重要依據(jù)。在區(qū)域試驗(yàn)布局中,試點(diǎn)分辨力的大小是一個(gè)重要的考慮因素,選擇具有代表性的試驗(yàn)點(diǎn)準(zhǔn)確評(píng)價(jià)參試品種的特性,不僅可以減少不必要的布點(diǎn)數(shù)量,同時(shí)也節(jié)省了人力和物力[1]。糯玉米是浙江省重要的鮮食旱糧作物,在浙江省糧食生產(chǎn)中占有重要地位。玉米農(nóng)藝性狀主要包括穗行數(shù)、行粒數(shù)、凈穗率、鮮千粒質(zhì)量、出籽率等農(nóng)藝指標(biāo),從不同側(cè)面詮釋了糯玉米的產(chǎn)量性狀。由于這些性狀受到環(huán)境、基因型(品種)以及基因與環(huán)境互作等多因子的共同影響,明確各因子對(duì)糯玉米農(nóng)藝性狀的貢獻(xiàn),對(duì)更好地完成區(qū)試工作,篩選優(yōu)良品種,提高浙江省糯玉米種植的經(jīng)濟(jì)效益具有重要意義。

    品種評(píng)價(jià)和試點(diǎn)指標(biāo)評(píng)價(jià)的準(zhǔn)確性直接影響著優(yōu)良品種的篩選和鑒定評(píng)價(jià)。品種的穩(wěn)定性和適應(yīng)性是評(píng)價(jià)品種優(yōu)劣的重要方面,同時(shí),試點(diǎn)的分辨力、代表性、穩(wěn)定性是衡量一個(gè)地點(diǎn)能否準(zhǔn)確評(píng)價(jià)參試品種特性的重要標(biāo)準(zhǔn),但是由于多環(huán)境(年份和地點(diǎn))下品種和環(huán)境的互作效應(yīng)、年份與地點(diǎn)的互作效應(yīng)形成機(jī)制的復(fù)雜性,關(guān)于區(qū)域試驗(yàn)數(shù)據(jù)的分析目前還沒(méi)有普遍認(rèn)可的最佳模型和方法[2]。長(zhǎng)期以來(lái),作物區(qū)域試驗(yàn)分析中較為常用的分析模型有Eberhart-Russell模型[3]、主效可加互作可乘模型(additive main effects and multiplicative interaction, AMMI)[4]、Finlay-Wilkinson模型[5]等。李本貴等[6]利用AMMI模型分析了小麥區(qū)域試驗(yàn),羅學(xué)芳等[7]采用Eberhart-Russell模型對(duì)四川省一年多點(diǎn)玉米區(qū)試進(jìn)行了聯(lián)合方差分析,何代元等[8]對(duì)黃淮海夏玉米區(qū)試產(chǎn)量采用AMMI模型分析,都表明基因與環(huán)境互作效應(yīng)達(dá)到極顯著水平,篩選出了豐產(chǎn)性高、穩(wěn)定性好的品種。雖然Eberhart-Russell模型和AMMI模型都能解釋品種與環(huán)境的互作效應(yīng),但張群遠(yuǎn)等[9]指出,Eberhart-Russell模型雖然應(yīng)用廣泛,但往往會(huì)碰到模型不適合的問(wèn)題;AMMI模型理論上精確度高,但參數(shù)個(gè)數(shù)不固定,且必須借助圖形分析,不能有效構(gòu)建簡(jiǎn)潔的指標(biāo)。目前,地點(diǎn)分辨力的分析方法主要有變異系數(shù)法、回歸系數(shù)法和AMMI模型Dj值法。樊龍江等[1]、紀(jì)榮昌等[10]采用AMMI模型Dj值法、變異系數(shù)法、回歸系數(shù)法對(duì)區(qū)域試驗(yàn)點(diǎn)進(jìn)行試點(diǎn)分辨力研究,結(jié)果表明:用AMMI模型Dj值法估算試點(diǎn)分辨力更科學(xué),可以得到環(huán)境與基因型獨(dú)立的分量估計(jì);變異系數(shù)法以表型為分析對(duì)象,而表型是基因與環(huán)境互作以及基因型綜合作用的結(jié)果,其參數(shù)估計(jì)不能獨(dú)立評(píng)估環(huán)境分量;回歸系數(shù)法假設(shè)的前提是互作效應(yīng)為線性函數(shù),應(yīng)用效果不太理想。上述模型和方法都不能具體地估算出基因與環(huán)境的互作效應(yīng)值、試點(diǎn)與年份的互作效應(yīng)值、地點(diǎn)效應(yīng)值等,特別是對(duì)多年多點(diǎn)的多環(huán)境區(qū)域試驗(yàn),無(wú)法對(duì)環(huán)境與基因型進(jìn)行獨(dú)立準(zhǔn)確的估計(jì),從而影響對(duì)試點(diǎn)分辨力、試點(diǎn)穩(wěn)定性、品種穩(wěn)定性的準(zhǔn)確性評(píng)價(jià)。

    本研究對(duì)浙江省糯玉米多年多點(diǎn)的田間區(qū)域試驗(yàn)結(jié)果進(jìn)行分析,采用基于混合線性模型的統(tǒng)計(jì)分析方法,解決了上述方法無(wú)法獨(dú)立估算基因與環(huán)境的互作效應(yīng)值、年份與地點(diǎn)互作效應(yīng)值等其他因子效應(yīng)值的問(wèn)題。利用該統(tǒng)計(jì)方法剖析試點(diǎn)、品種、年份與試點(diǎn)互作以及品種與試點(diǎn)互作等因子對(duì)糯玉米產(chǎn)量性狀以及產(chǎn)量相關(guān)性狀的貢獻(xiàn)率,分析各品種的豐產(chǎn)性、穩(wěn)定性、適應(yīng)性,各試點(diǎn)的穩(wěn)定性、代表性、分辨力,為區(qū)試數(shù)據(jù)統(tǒng)計(jì)分析提供一套科學(xué)的統(tǒng)計(jì)方法,并提出理想?yún)^(qū)試點(diǎn)的建議。

    1 材料與方法

    1.1 試驗(yàn)設(shè)計(jì)

    田間試驗(yàn)于2013、2014、2015年實(shí)施,共8個(gè)試點(diǎn):淳安(L1)、東陽(yáng)(L2)、嘉善(L3)、江山(L4)、寧海(L5)、杭州(L6)、嵊州(L7)、仙居(L8)。29個(gè)參試糯玉米品種分別為:麗都1號(hào)(G1)、綠玉糯2號(hào)(G2)、美玉8號(hào)(G3)、浦甜糯1號(hào)(G4)、錢江糯3號(hào)(G5)、天貴糯161(G6)、天糯828 (G7)、新甜糯88(G8)、浙糯1202(G9)、浙糯1302(G10)、花甜糯K10-1(G11)、花甜糯8013(G12)、豇玉糯6號(hào)(G13)、京科糯569(G14)、京科糯929(G15)、科甜糯2號(hào)(G16)、南粵花糯(G17)、天香3號(hào)(G18)、浙糯1201(G19)、彩花糯11-3(G20)、黑糯181(G21)、花甜糯001(G22)、花甜糯2號(hào)(G23)、科糯082(G24)、美玉15號(hào)(G25)、美玉16號(hào)(G26)、錢江糯2號(hào)(G27)、甜彩糯1號(hào)(G28)、浙甜糯615(G29)。其中,2015年8個(gè)試點(diǎn)種植品種為G1~G10共10個(gè)品種,2014年8個(gè)試點(diǎn)(產(chǎn)量性狀不包含L3)種植品種為G11~G19、G3、G5、G8、G9共13個(gè)品種,2013年8個(gè)試點(diǎn)種植品種為G20~G29、G3、G11、G14、G16共14個(gè)品種。

    田間種植采用隨機(jī)區(qū)組設(shè)計(jì),3次重復(fù),小區(qū)面積20 m2,密度為4 000株/667 m2,播種、移栽等田間管理按當(dāng)?shù)亓?xí)慣進(jìn)行。

    1.2 分析方法

    田間試驗(yàn)農(nóng)藝性狀選取糯玉米產(chǎn)量性狀、穗粗、穗長(zhǎng)、穗行數(shù)、行粒數(shù)、凈穗率、鮮千粒質(zhì)量、出籽率等指標(biāo)。產(chǎn)量性狀(其他農(nóng)藝性狀不含區(qū)組效應(yīng))采取以下混合線性模型進(jìn)行方差分析和效應(yīng)值估算[11-12]:

    其中:Yijkh為第h個(gè)年份第j個(gè)試點(diǎn)第k個(gè)區(qū)組內(nèi)第i個(gè)品種的性狀表型值;u為總體均值;Gi為第i個(gè)品種的品種效應(yīng)(隨機(jī)效應(yīng)),Gi~(0,σG2);Lj為第j個(gè)試點(diǎn)的效應(yīng)(隨機(jī)效應(yīng)),Lj~(0,σL2);Yh為第h個(gè)年份的效應(yīng)(隨機(jī)效應(yīng)),Yh~(0,σh2);GLij為第i個(gè)品種與第j個(gè)地點(diǎn)的互作效應(yīng)(隨機(jī)效應(yīng)),GLij~(0,σGL2);YLjh為第h個(gè)年份與第j個(gè)地點(diǎn)的互作效應(yīng)(隨機(jī)效應(yīng)),YLjh~(0,σYL2);Bk(jh)為第h個(gè)年份第j個(gè)地點(diǎn)內(nèi)第k個(gè)區(qū)組效應(yīng)(隨機(jī)效應(yīng)),Bk(jh)~(0,σB2);eijkh為殘差效應(yīng)(隨機(jī)效應(yīng)),eijkh~(0,σ2e)。

    采取最小二次無(wú)偏估計(jì)(MINQUE)法估算模型(1)中各效應(yīng)的方差分量,基于方差分量采用調(diào)整的無(wú)偏預(yù)測(cè)法(AUP)預(yù)測(cè)隨機(jī)效應(yīng)值?;陬A(yù)測(cè)出的效應(yīng)值進(jìn)行試點(diǎn)分辨力、試點(diǎn)穩(wěn)定性、品種穩(wěn)定性、品種豐產(chǎn)性等指標(biāo)分析。對(duì)預(yù)測(cè)出的第j個(gè)試點(diǎn)與品種的互作效應(yīng)值(GLij)、第i個(gè)品種效應(yīng)值(Gi)與該品種和各個(gè)地點(diǎn)互作效應(yīng)值(GLij)的合成值分別進(jìn)行方差分析,估計(jì)σGL2和σi

    2,并作為衡量試點(diǎn)分辨力和品種穩(wěn)定性的指標(biāo)。對(duì)于品種豐產(chǎn)性評(píng)價(jià),對(duì)模型預(yù)測(cè)出的品種效應(yīng)值用Tukey法進(jìn)行多重比較,適宜品種的選擇標(biāo)準(zhǔn)為多重比較結(jié)果中標(biāo)有“a”的品種。在適宜的品種中,選取效應(yīng)值高于所有適宜品種效應(yīng)值的品種作為最適宜品種?;谀P?1),將試點(diǎn)效應(yīng)(Lj)作為固定效應(yīng),品種效應(yīng)(Gi)作為隨機(jī)效應(yīng),對(duì)預(yù)測(cè)出的Lj進(jìn)行多重比較,可以篩選出最適宜的地點(diǎn);利用第j個(gè)試點(diǎn)效應(yīng)值(Lj)和該試點(diǎn)與不同年份的互作效應(yīng)值(YLjh)的合成值進(jìn)行方差分析,得到的方差σj2作為衡量試點(diǎn)穩(wěn)定性的標(biāo)準(zhǔn)。

    分別對(duì)每一試點(diǎn)的性狀觀察值采取以下混合線性模型進(jìn)行試驗(yàn)精確度分析[11-12]:

    其中:Yikh為第h個(gè)年份第k個(gè)區(qū)組內(nèi)第i個(gè)品種的性狀值;GYih為品種與年份的互作效應(yīng)(隨機(jī)效應(yīng)),GYih~(0,σGY2);Bk(h)為第h個(gè)年份內(nèi)第k個(gè)區(qū)組效應(yīng)(隨機(jī)效應(yīng)),Bk(h)~(0,σB2);其余效應(yīng)值與模型(1)中的效應(yīng)值一致。對(duì)模型(2)方差分析得出的誤差項(xiàng)均方(mean square error,MSE)按照陳洪梅等[13]的方法計(jì)算試驗(yàn)誤差變異系數(shù)(coefficient of error variation,CEV),公式如下:

    式中:CEV為衡量試驗(yàn)精確度的指標(biāo);Yˉ為單個(gè)試點(diǎn)所有品種的性狀均值。

    統(tǒng)計(jì)分析采用QTModel軟件[14](http://ibi.zju. edu.cn/software)和SAS 9.0軟件的MIXED模塊。

    2 結(jié)果與分析

    2.1 農(nóng)藝性狀方差分析

    糯玉米產(chǎn)量、穗粗、穗長(zhǎng)、穗行數(shù)、行粒數(shù)、凈穗率、鮮千粒質(zhì)量、出籽率等農(nóng)藝性狀指標(biāo)方差分析見(jiàn)表1。結(jié)果表明:品種效應(yīng)除凈穗率外,其他性狀均達(dá)到極顯著水平,貢獻(xiàn)率在7.4%~67.6%之間,說(shuō)明品種效應(yīng)對(duì)各農(nóng)藝性狀影響差異較大,對(duì)各農(nóng)藝性狀的貢獻(xiàn)程度大小為穗行數(shù)(67.6%)>穗長(zhǎng)(39.9%)>行粒數(shù)(29.4%)>穗粗(23.1%)>產(chǎn)量(14.7%)>鮮千粒質(zhì)量(13.8%)>出籽率(13.6%)>凈穗率(7.4%),可以看出穗行數(shù)和穗長(zhǎng)主要受品種影響,品種間差異較大,凈穗率基本不受品種影響;年份效應(yīng)對(duì)產(chǎn)量、穗長(zhǎng)、鮮千粒質(zhì)量達(dá)到顯著水平,其余農(nóng)藝性狀未達(dá)到顯著水平,對(duì)鮮千粒質(zhì)量貢獻(xiàn)率達(dá)到17.2%,對(duì)其余農(nóng)藝性狀貢獻(xiàn)率很小(基本都在4%以下),說(shuō)明年份效應(yīng)對(duì)鮮千粒質(zhì)量之外的農(nóng)藝性狀影響很?。坏攸c(diǎn)效應(yīng)除穗行數(shù)未達(dá)到顯著水平外,其他性狀均達(dá)到極顯著水平,貢獻(xiàn)率在1.9%~36.3%之間,說(shuō)明地點(diǎn)效應(yīng)對(duì)各農(nóng)藝性狀影響差異比較大,對(duì)各農(nóng)藝性狀的貢獻(xiàn)程度大小為產(chǎn)量(36.3%)>行粒數(shù)(29.3%)>鮮千粒質(zhì)量(27.5%)>出籽率(19.1%)>穗粗(15.4%)>凈穗率(13.0%)>穗長(zhǎng)(8.3%)>穗行數(shù)(1.9%),可以看出產(chǎn)量受地點(diǎn)影響較大,穗行數(shù)和穗長(zhǎng)基本不受地點(diǎn)影響;品種與地點(diǎn)互作效應(yīng)只有產(chǎn)量性狀達(dá)到極顯著水平,貢獻(xiàn)率為10.3%,其他性狀未達(dá)到顯著水平,貢獻(xiàn)率基本在5%以下,說(shuō)明品種與地點(diǎn)互作效應(yīng)對(duì)產(chǎn)量性狀影響較大,對(duì)其他性狀影響較小;年份和地點(diǎn)互作效應(yīng)除穗粗和穗行數(shù)達(dá)到顯著水平外,其余性狀均達(dá)到極顯著水平,貢獻(xiàn)率基本都在15%以上,大小依次為出籽率(37.1%)>產(chǎn)量(25.7%)>凈穗率(20.7%)>行粒數(shù)(16.5%)>穗粗(15.4%)>穗長(zhǎng)(15.1%)>鮮千粒質(zhì)量(10.8%)>穗行數(shù)(5.0%),說(shuō)明年份與地點(diǎn)互作效應(yīng)對(duì)穗行數(shù)影響最小,對(duì)其余性狀都有較大影響,特別是出籽率、產(chǎn)量和凈穗率。地點(diǎn)相關(guān)效應(yīng)(年份與地點(diǎn)互作效應(yīng)+地點(diǎn)效應(yīng))貢獻(xiàn)率在6.9%~62.0%之間,不同農(nóng)藝性狀之間地點(diǎn)相關(guān)效應(yīng)貢獻(xiàn)率差異較大,其中對(duì)產(chǎn)量貢獻(xiàn)率最高,為62.0%,對(duì)構(gòu)成產(chǎn)量相關(guān)農(nóng)藝性狀的貢獻(xiàn)率依次為出籽率(56.2%)>行粒數(shù)(45.8%)>鮮千粒質(zhì)量(38.3%)>凈穗率(33.7%)>穗粗(30.8%)>穗長(zhǎng)(23.4%)>穗行數(shù)(6.9%),說(shuō)明產(chǎn)量、出籽率、行粒數(shù)主要受地點(diǎn)相關(guān)效應(yīng)影響。

    表1 糯玉米農(nóng)藝性狀效應(yīng)方差估計(jì)值及表型方差百分率Table1 Estimated variances and theirpercentages ofphenotypic variance forwaxy maize agronomic traits

    2.2 品種評(píng)價(jià)

    品種評(píng)價(jià)主要是對(duì)品種的穩(wěn)定性、生產(chǎn)水平和適應(yīng)性等指標(biāo)進(jìn)行評(píng)價(jià)。優(yōu)良品種不僅要在有利的條件下保持高產(chǎn),在不利的條件下也應(yīng)具有相對(duì)高產(chǎn)的能力。所以,一個(gè)應(yīng)用范圍廣的品種應(yīng)具有高產(chǎn)穩(wěn)產(chǎn)、適應(yīng)性強(qiáng)的基本特征。根據(jù)分析結(jié)果可以看出,基于產(chǎn)量性狀分析的糯玉米品種數(shù)量多達(dá)29個(gè),所以在實(shí)際生產(chǎn)中建議選取最適宜或者部分適宜的品種,并且參考其他相關(guān)指標(biāo)(穩(wěn)定性)來(lái)決定是否推廣種植。

    2.2.1 品種豐產(chǎn)性

    品種豐產(chǎn)性是指品種的生產(chǎn)水平,即一個(gè)品種在整個(gè)區(qū)域中的平均生產(chǎn)能力。按照模型(1)對(duì)糯玉米產(chǎn)量數(shù)據(jù)進(jìn)行分析,對(duì)品種的產(chǎn)量性狀均值(小區(qū)產(chǎn)量)用Tukey法進(jìn)行多重比較,結(jié)果見(jiàn)表2。根據(jù)多重比較結(jié)果,選取標(biāo)有“a”的品種作為產(chǎn)量性狀較優(yōu)者,適宜在浙江省推廣和種植,其中標(biāo)有“*”的品種為最適宜品種,宜重點(diǎn)推廣和示范。從表2中可以看出,產(chǎn)量性狀適宜性好的品種有天糯828 (G7)、科甜糯2號(hào)(G16)、京科糯569(G14)、新甜糯88(G8)等4個(gè)品種。

    2.2.2 品種穩(wěn)定性

    品種穩(wěn)定性是指品種在不同環(huán)境下表現(xiàn)出的變異特性。用模型(1)預(yù)測(cè)出的各品種效應(yīng)值(Gi)、品種與試點(diǎn)的互作效應(yīng)值(GLij)計(jì)算,評(píng)價(jià)各品種的地點(diǎn)穩(wěn)定性。從表2可以看出,各品種的值有所差別,其中:浙糯1201(G19)、南粵花糯(G17)、科甜糯2號(hào)(G16)、美玉8號(hào)(G3)等4個(gè)品種的值較小,穩(wěn)定性較高;花甜糯8013(G12)、美玉15號(hào)(G25)、天貴糯161(G6)、浙糯1202(G9)等4個(gè)品種的值較大,穩(wěn)定性較差。

    2.3 試點(diǎn)評(píng)價(jià)

    2.3.1 試點(diǎn)總體精度

    試點(diǎn)精度是指試驗(yàn)中相同處理的重復(fù)觀測(cè)值間相互接近的程度,直接反映了區(qū)域試驗(yàn)的質(zhì)量;試驗(yàn)精確度越高,品種間差異的鑒別越準(zhǔn)確。以試點(diǎn)產(chǎn)量為依據(jù)進(jìn)行統(tǒng)計(jì)分析,用模型(2)計(jì)算出誤差項(xiàng)均方(MSE),用式(3)計(jì)算誤差變異系數(shù)(CEV),結(jié)果見(jiàn)表3。誤差變異系數(shù)(CEV)越小,試驗(yàn)精確度越高。一般而言,田間試驗(yàn)的CEV若小于10%,說(shuō)明試驗(yàn)誤差控制得較好,通常在區(qū)域試驗(yàn)中可容忍的CEV為12%。從表3中可以看出,該輪區(qū)域試驗(yàn)(8個(gè)試點(diǎn))的 CEV均低于10%,試驗(yàn)的總體CEV為4.63%,小于10%:說(shuō)明浙江省糯玉米品種區(qū)域試驗(yàn)的整體水平較高,試驗(yàn)誤差控制得較好,數(shù)據(jù)可信。

    表2 各品種產(chǎn)量均值多重比較和穩(wěn)定性分析Table2 Multiple comparison of the mean value and analysis of the Table 2 stability for yield ofvariety

    2.3.2 試點(diǎn)生產(chǎn)水平

    試點(diǎn)生產(chǎn)水平是指一個(gè)試點(diǎn)產(chǎn)量水平的高低,反映了該試點(diǎn)的環(huán)境條件(包括試驗(yàn)條件)是否適宜參試品種的生長(zhǎng),一般用一個(gè)試點(diǎn)所有參試品種的產(chǎn)量均值來(lái)表示。2013、2014和2015年區(qū)試各試點(diǎn)所有參試品種的產(chǎn)量均值見(jiàn)表3。從中可以看出:8個(gè)試點(diǎn)產(chǎn)量均值存在較大差異,說(shuō)明不同試點(diǎn)的生產(chǎn)水平有高有低,平均產(chǎn)量越高,生產(chǎn)水平越高;嵊州的生產(chǎn)水平最高,杭州的生產(chǎn)水平最低,其余試點(diǎn)的生產(chǎn)水平中等。

    2.3.3 試點(diǎn)分辨力

    試點(diǎn)分辨力是指一個(gè)試點(diǎn)上參試品種所表現(xiàn)的變異大小。分辨力高的試點(diǎn),品種間的變異大,表明該試點(diǎn)的環(huán)境條件有利于充分展現(xiàn)參試品種間的差異,有利于不同品種的比較。地點(diǎn)分辨力研究主要是考察基因型(品種)與地點(diǎn)互作效應(yīng)的大小,對(duì)產(chǎn)量性狀單個(gè)試點(diǎn)各個(gè)品種與地點(diǎn)互作效應(yīng)進(jìn)行方差分析求出σGL2,其值越大,該試點(diǎn)對(duì)品種的分辨能力越強(qiáng)。從表3中可以看出:各個(gè)試點(diǎn)對(duì)品種的分辨能力有所差別,嘉善、杭州、淳安等3個(gè)試點(diǎn)的σGL

    2值均較大,具有較強(qiáng)的品種分辨力,其中嘉善的分辨力最強(qiáng);寧海、仙居2個(gè)試點(diǎn)的σGL2值較小,對(duì)參試品種的分辨力較弱。

    2.3.4 試點(diǎn)代表性

    試點(diǎn)代表性分為2種,分別為普遍代表性和特殊代表性。普遍代表性是指一個(gè)試點(diǎn)上的試驗(yàn)結(jié)果與整個(gè)區(qū)試結(jié)果的一致程度,反映了試點(diǎn)對(duì)整個(gè)區(qū)試環(huán)境的代表程度。相反,特殊代表性是指一個(gè)試點(diǎn)上的試驗(yàn)結(jié)果與整個(gè)區(qū)試結(jié)果的不一致程度,反映了試點(diǎn)對(duì)某一類特殊環(huán)境的代表程度。試點(diǎn)代表性用每一個(gè)地點(diǎn)的各品種表型值與各品種總體均值的相關(guān)系數(shù)來(lái)表示,相關(guān)系數(shù)越大,試點(diǎn)普遍代表性越高,反之,試點(diǎn)特殊代表性越高。用模型(1)預(yù)測(cè)的產(chǎn)量性狀品種效應(yīng)值(Gi)、地點(diǎn)效應(yīng)值(Lj)、品種與地點(diǎn)互作效應(yīng)值(GLij)合成值之間的相關(guān)系數(shù)(表3)作為評(píng)價(jià)試點(diǎn)代表性的標(biāo)準(zhǔn)。從中可以看出:東陽(yáng)、江山、淳安、杭州4個(gè)試點(diǎn)的相關(guān)系數(shù)較高,普遍代表性較好,宜鑒別廣泛適應(yīng)性的品種;嵊州、嘉善2個(gè)試點(diǎn)的相關(guān)系數(shù)較低,特殊代表性較高,適于鑒別一些特殊性狀的品種。

    2.3.5 試點(diǎn)穩(wěn)定性和效應(yīng)值綜合評(píng)價(jià)

    試點(diǎn)穩(wěn)定性是指一個(gè)試點(diǎn)在不同年份試驗(yàn)中表現(xiàn)出的變異程度,反映了試點(diǎn)的試驗(yàn)環(huán)境條件的穩(wěn)定性,包括試點(diǎn)環(huán)境對(duì)不同氣候條件的緩沖能力[9]。選取糯玉米產(chǎn)量性狀,用模型(1)預(yù)測(cè)出的地點(diǎn)效應(yīng)均值(Lj)和地點(diǎn)與年份互作效應(yīng)(GLij)計(jì)算,評(píng)價(jià)各試點(diǎn)的穩(wěn)定性,以及對(duì)地點(diǎn)效應(yīng)均值(Lj)用Tukey法進(jìn)行多重比較。值越大,試點(diǎn)穩(wěn)定性越差,越小,試點(diǎn)穩(wěn)定性越高。從表3中可以看出:嵊州產(chǎn)量均值最高,但穩(wěn)定性最差;東陽(yáng)穩(wěn)定性最好,產(chǎn)量均值中等偏高;嘉善穩(wěn)定性和產(chǎn)量均值處于中等水平。

    表3 試點(diǎn)精度、平均產(chǎn)量、分辨力、代表性、穩(wěn)定性分析及試點(diǎn)效應(yīng)的多重比較Table3 Analysis ofprecision,mean yield,distinguishability,representativeness,stability oflocations and multiple comparison oflocation effect

    3 討論

    糯玉米的農(nóng)藝性狀受到基因型、環(huán)境(年份和地點(diǎn))以及基因型與環(huán)境互作等多因子共同影響。前人的研究多注重平均值的比較,不涉及多產(chǎn)區(qū),品種、地點(diǎn)效應(yīng)作為固定效應(yīng)。區(qū)域試驗(yàn)中采用不同的分析模型對(duì)結(jié)果有較大影響,甚至達(dá)到嚴(yán)重差異。通常,區(qū)試的目標(biāo)在于篩選出性狀綜合表現(xiàn)較好、在生態(tài)區(qū)可推廣應(yīng)用的新品種,以及適合區(qū)試的理想試點(diǎn)。因此,一般將品種作為固定效應(yīng),而年份、地點(diǎn)及其與品種的互作效應(yīng)都作為隨機(jī)效應(yīng),當(dāng)參試品種數(shù)目較多時(shí),品種也可作隨機(jī)效應(yīng)處理。無(wú)論是多年多點(diǎn)、多年單點(diǎn)或單年多點(diǎn)的作物區(qū)域試驗(yàn),采用基于混合線性模型的分析策略較為合適。

    目前,多環(huán)境區(qū)試資料分析常采用AMMI模型[7-8],能分析品種與環(huán)境的互作效應(yīng),但AMMI模型是一個(gè)固定模型,只適用于分析單年多點(diǎn)的區(qū)試,不適合分析多年多點(diǎn)的試點(diǎn)資料。胡希遠(yuǎn)等[2,15]指出,Eberhart-Russell模型、AMMI模型并不適合所有試驗(yàn)資料,而混合線性模型可借助多種方差-協(xié)方差結(jié)構(gòu)來(lái)反映各類試驗(yàn)數(shù)據(jù),能很好適應(yīng)方差同質(zhì)和異質(zhì)的情況。本研究是一個(gè)多年多點(diǎn)的區(qū)域試驗(yàn),采用基于混合線性模型的方法,運(yùn)用QTModel軟件和SAS軟件分析試驗(yàn)數(shù)據(jù)。基于預(yù)測(cè)出的各效應(yīng)分量值分析試點(diǎn)穩(wěn)定性、代表性、分辨力、均值差異以及品種穩(wěn)定性、均值差異。對(duì)于各試點(diǎn)的穩(wěn)定性,一般用不同年份間同一組參試品種性狀的相關(guān)系數(shù)表示;試點(diǎn)分辨力,一般用試點(diǎn)對(duì)品種指數(shù)的回歸系數(shù)或各試點(diǎn)參試品種間的遺傳變異系數(shù)來(lái)表示;品種的穩(wěn)定性,采用各試點(diǎn)觀測(cè)值對(duì)環(huán)境指數(shù)的簡(jiǎn)單回歸進(jìn)行分析[3],但是這些分析方法都沒(méi)有將具體的效應(yīng)值分離進(jìn)行計(jì)算,無(wú)法獨(dú)立分析試點(diǎn)穩(wěn)定性、分辨力以及品種穩(wěn)定性。我們通過(guò)估算品種效應(yīng)值、試點(diǎn)效應(yīng)值,以及品種與試點(diǎn)互作效應(yīng)值、試點(diǎn)與年份的效應(yīng)值,對(duì)各具體指標(biāo)進(jìn)行了獨(dú)立分析。例如,基于各試點(diǎn)效應(yīng)值、試點(diǎn)與各年份的互作效應(yīng)值進(jìn)行試點(diǎn)的穩(wěn)定性評(píng)價(jià),基于各品種效應(yīng)值、品種與不同地點(diǎn)的互作效應(yīng)值進(jìn)行品種的穩(wěn)定性評(píng)價(jià),使得分析更為簡(jiǎn)潔直觀。

    根據(jù)浙江省2013—2015年糯玉米區(qū)試結(jié)果,我們對(duì)浙江省8個(gè)區(qū)試點(diǎn)29個(gè)參試品種進(jìn)行了綜合評(píng)價(jià):8個(gè)試點(diǎn)的生產(chǎn)水平有較大差異,其中生產(chǎn)水平較高的試點(diǎn)為嵊州和江山;8個(gè)區(qū)試點(diǎn)的試驗(yàn)誤差變異系數(shù)均較小,試驗(yàn)均有效;嘉善、淳安等2個(gè)試點(diǎn)具有較強(qiáng)的品種鑒別力,寧海、仙居等2個(gè)試點(diǎn)對(duì)品種鑒別力較弱;東陽(yáng)、江山、淳安、杭州等4個(gè)試點(diǎn)普遍代表性高,適于鑒別具有廣泛適應(yīng)性的品種,是理想的糯玉米區(qū)試點(diǎn)和育種點(diǎn);嵊州、嘉善等2個(gè)試點(diǎn)特殊代表性高,適于鑒別一些特殊性狀的品種;東陽(yáng)、江山等2個(gè)試點(diǎn)在產(chǎn)量性狀上穩(wěn)定性好、生產(chǎn)水平高。通過(guò)產(chǎn)量性狀篩選出適宜浙江省種植的品種為天糯828、科甜糯2號(hào)、京科糯569、新甜糯88等4個(gè)品種,穩(wěn)定性較好的品種為浙糯1201、南粵花糯、科甜糯2號(hào)、美玉8號(hào)等4個(gè)品種,綜合考慮穩(wěn)定性及產(chǎn)量因素,認(rèn)為科甜糯2號(hào)為適宜浙江省種植的糯玉米品種。

    [1]樊龍江,胡秉民.作物區(qū)域試驗(yàn)點(diǎn)區(qū)辨力估算方法的比較研究.生物數(shù)學(xué)學(xué)報(bào),2000,15(2):175-179.

    FAN L J,HU B M.Comparison of different methods for estimating location resolution in crop regional trial.Journal of Biomathematics,2000,15(2):175-179.(in Chinese)

    [2]胡希遠(yuǎn),尤海磊,任長(zhǎng)虹,等.基于協(xié)方差陣結(jié)構(gòu)優(yōu)選的作物品種區(qū)域試驗(yàn)分析.作物學(xué)報(bào),2009,35(11):1981-1989.

    HU X Y,YOU H L,REN C H,et al.Analysis of crop variety regional trials based on selection of covariance structures.Acta Agronomica Sinica,2009,35(11):1981-1989.(in Chinese with English abstract)

    [3]EBERHART S A,RUSSELL W A.Stability parameters for comparing varieties.Crop Science,1966,6(1):36-40.

    [4]張群遠(yuǎn),孔繁玲.作物品種區(qū)域試驗(yàn)統(tǒng)計(jì)分析模型的比較.中國(guó)農(nóng)業(yè)科學(xué),2002,35(4):365-371.

    ZHANG Q Y,KONG F L.Comparison of statistical models for regional crop trial analysis.Scientia Agricultura Sinica,2002,35 (4):365-371.(in Chinese with English abstract)

    [5]OMAN S D.Multiplicative effects in mixed model analysis of variance.Biometrika,1991,78(4):729-739.

    [6]李本貴,閻俊,何中虎,等.用AMMI模型分析作物區(qū)域試驗(yàn)中的地點(diǎn)鑒別力.作物學(xué)報(bào),2004,30(6):593-596.

    LI B G,YAN J,HE Z H,et al.Analyzing site discrimination in crop regionalyield trialby AMMImodel.Acta Agronomica Sinica, 2004,30(6):593-596.(in Chinese with English abstract)

    [7]羅學(xué)芳,茍才明.Eberhart和Russell模型在玉米區(qū)域試驗(yàn)中的應(yīng)用.西南農(nóng)業(yè)學(xué)報(bào),2010,23(2):602-604.

    LUO X F,GOU C M.Application of Eberhart and Russellmodel in regional experiment of maize.Southwest China Journal ofAgricultural Sciences,2010,23(2):602-604.(in Chinese with English abstract)

    [8]何代元,胡寧,馬兆錦,等.AMMI模型在玉米區(qū)域試驗(yàn)中的應(yīng)用.玉米科學(xué),2009,17(4):144-147.

    HE D Y,HU N,MA Z J,et al.Application of AMMI model in regional experiment of maize.Journal of Maize Science, 2009,17(4):144-147.(in Chinese with English abstract)

    [9]張群遠(yuǎn),孔繁玲,廖琴,等.作物品種區(qū)域試驗(yàn)的評(píng)價(jià)體系及評(píng)價(jià)方法.農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2000,16(2):81-86.

    ZHANG Q Y,KONG F L,LIAO Q,et al.Evaluation of system and methods for regional crop trials.System Sciences and Comprehensive Studies in Agriculture,2000,16(2):81-86.(in Chinese with English abstract)

    [10]紀(jì)榮昌,盧和頂,陳山虎,等.甜玉米區(qū)域試驗(yàn)點(diǎn)對(duì)品種鑒別能力估計(jì)方法的研究.中國(guó)農(nóng)學(xué)通報(bào),2008,24(10):153-156.

    JI R C,LU H D,CHEN S H,et al.Studies on methods of discriminating ability of location for sweet maize regional tests. Chinese Agricultural Science Bulletin,2008,24(10):153-156.(in Chinese with English abstract)

    [11]朱軍,許馥華,賴?guó)Q岡.作物品種區(qū)域試驗(yàn)非平衡資料的分析方法:?jiǎn)我恍誀畹姆治?浙江農(nóng)業(yè)大學(xué)學(xué)報(bào),1993,19(1):7-13.

    ZHU J,XU F H,LAI M G.Analysis method for unbalanced data from regional trials of crop variety:Analysis for single trait. Journal of Zhejiang Agricultural University,1993,19(1):7-13.(in Chinese with English abstract)

    [12]朱軍,賴?guó)Q岡,許馥華.作物品種區(qū)域試驗(yàn)非平衡資料的分析方法:綜合性狀的分析.浙江農(nóng)業(yè)大學(xué)學(xué)報(bào),1993,19(3):241-247.

    ZHU J,LAIM G,XU F H.Analysis methods for unbalanced data from regional trial of crop variety:Analysis for multiple traits. Journal of Zhejiang Agricultural University,1993,19(3):241-247. (in Chinese with English abstract)

    [13]陳洪梅,楊峻蕓,譚靜.玉米區(qū)域試驗(yàn)精確度分析及品種的灰色綜合評(píng)判.西南農(nóng)業(yè)學(xué)報(bào),2004,17(增刊1):228-232.

    CHEN H M,YANG J Y,TAN J.Analysis of experimental precision and grey integrative evaluation of varieties in maize regional trials.Southwest China Journal of Agricultural Sciences, 2004,17(Suppl.1):228-232.(in Chinese with English abstract)

    [14]YANG J,ZOU Y,ZHU J.Identifying differentially expressed genes in human acute leukemia and mouse brain microarray datasets utilizing QTModel.Functional and Integrative Genomics, 2009,9(1):59-66.

    [15]胡希遠(yuǎn),尤海磊,宋喜芳,等.作物品種穩(wěn)定性分析不同模型的比較.麥類作物學(xué)報(bào),2009,29(1):110-117.

    HU X Y,YOU H L,SONG X F,et al.Comparison of different models for crop stability analysis.Journal of Triticeae Crops, 2009,29(1):110-117.(in Chinese with English abstract)

    Comprehensive evaluation on the stability,adaptability and location for waxy maize varieties in Zhejiang Province.Journal ofZhejiang University(Agric.&Life Sci.),2017,43(3):281-288

    LI Yan1,LIN Feng2,LI Lulu2,XU Haiming2*
    (1.Seed Management Station of Zhejiang Province,Hangzhou 310020,China; 2.Institute of Crop Science,College of Agriculture and Biotechnology,Zhejiang University,Hangzhou 310058,China)

    waxy maize;agronomic trait;regionaltrial;location evaluation;variety evaluation;mixed linearmodel

    ):李燕(http://orcid.org/0000-0002-6078-6944),E-mail:liyan00@163.com

    2016-11-21;接受日期(Accepted):2016-12-08

    S 513;S 338

    A

    10.3785/j.issn.1008-9209.2016.11.211

    浙江省“三農(nóng)六方”科技協(xié)作項(xiàng)目[CTZB-F150922AWZ-SNY1(2)]。

    *通信作者(Corresponding author):徐海明(http://orcid.org/0000-0002-0675-4087),E-mail:hmxu@zju.edu.cn

    猜你喜歡
    糯玉米區(qū)域試驗(yàn)性狀
    寶鐸草的性狀及顯微鑒定研究
    淺論糯玉米苗期田間管理
    9種常用中藥材的性狀真?zhèn)舞b別
    2014-2015國(guó)家甘薯品種南方區(qū)域試驗(yàn)廣西試點(diǎn)報(bào)告
    黑糯玉米果凍的加工工藝研究
    廣西引種臺(tái)灣榿木區(qū)域試驗(yàn)
    黑糯玉米新品種‘華耘黑糯501’的選育
    對(duì)“性狀分離比模擬”實(shí)驗(yàn)的改進(jìn)
    庫(kù)塔墾區(qū)早中熟陸地雜交棉品種區(qū)域試驗(yàn)
    陸地棉數(shù)量性狀的多元統(tǒng)計(jì)分析
    亚洲精品中文字幕一二三四区| 久久婷婷成人综合色麻豆| 麻豆久久精品国产亚洲av | 午夜精品久久久久久毛片777| 女人精品久久久久毛片| 男男h啪啪无遮挡| 女同久久另类99精品国产91| 男女下面进入的视频免费午夜 | 久久久久久人人人人人| 午夜免费观看网址| 中文欧美无线码| svipshipincom国产片| 成人永久免费在线观看视频| av天堂在线播放| 国产黄a三级三级三级人| 午夜福利在线观看吧| 老司机福利观看| av电影中文网址| 久久久久亚洲av毛片大全| 国产一区二区在线av高清观看| 国产色视频综合| 正在播放国产对白刺激| 99久久99久久久精品蜜桃| 日日爽夜夜爽网站| 亚洲五月婷婷丁香| 别揉我奶头~嗯~啊~动态视频| 夜夜爽天天搞| 99精国产麻豆久久婷婷| 91麻豆精品激情在线观看国产 | 麻豆av在线久日| 国产黄色免费在线视频| 亚洲五月婷婷丁香| 99热国产这里只有精品6| 一级作爱视频免费观看| 99热国产这里只有精品6| 免费女性裸体啪啪无遮挡网站| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 免费日韩欧美在线观看| 国产又爽黄色视频| 亚洲aⅴ乱码一区二区在线播放 | 黄色视频不卡| 老司机午夜十八禁免费视频| 超碰成人久久| 免费av毛片视频| 一级,二级,三级黄色视频| 午夜成年电影在线免费观看| 午夜激情av网站| 精品一品国产午夜福利视频| 色精品久久人妻99蜜桃| 淫妇啪啪啪对白视频| 亚洲精品在线美女| 亚洲中文av在线| 免费在线观看黄色视频的| 可以在线观看毛片的网站| 国产成人欧美| 成在线人永久免费视频| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品综合一区在线观看 | 夫妻午夜视频| 欧美成人午夜精品| a级片在线免费高清观看视频| 99riav亚洲国产免费| 国产精品一区二区在线不卡| 精品一区二区三区视频在线观看免费 | 国产野战对白在线观看| 免费人成视频x8x8入口观看| 欧美大码av| 日本黄色视频三级网站网址| 久久久久久久午夜电影 | 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 婷婷丁香在线五月| 亚洲专区国产一区二区| 亚洲精品美女久久av网站| 男人操女人黄网站| 久久精品国产综合久久久| 国产男靠女视频免费网站| 亚洲九九香蕉| 多毛熟女@视频| 亚洲精品成人av观看孕妇| 免费在线观看完整版高清| 欧美成狂野欧美在线观看| 精品国产乱子伦一区二区三区| 亚洲国产欧美日韩在线播放| 天天添夜夜摸| av网站在线播放免费| 老司机在亚洲福利影院| 亚洲久久久国产精品| 久久香蕉激情| 在线国产一区二区在线| 国产成人系列免费观看| 高清毛片免费观看视频网站 | 欧美黑人精品巨大| 两性夫妻黄色片| 一级毛片女人18水好多| 色哟哟哟哟哟哟| 久久人人97超碰香蕉20202| 校园春色视频在线观看| 免费观看精品视频网站| 老司机深夜福利视频在线观看| 水蜜桃什么品种好| 日韩欧美国产一区二区入口| av在线播放免费不卡| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久久毛片微露脸| 日韩高清综合在线| 免费在线观看亚洲国产| 在线播放国产精品三级| 韩国av一区二区三区四区| 一级a爱视频在线免费观看| 欧美日韩亚洲高清精品| 在线视频色国产色| 他把我摸到了高潮在线观看| 久久久国产精品麻豆| 麻豆av在线久日| 在线av久久热| 夜夜躁狠狠躁天天躁| 久久久国产成人精品二区 | 日韩精品中文字幕看吧| 色在线成人网| 黑丝袜美女国产一区| 麻豆国产av国片精品| 久久久久久久久免费视频了| 女性被躁到高潮视频| 神马国产精品三级电影在线观看 | 亚洲一区中文字幕在线| 久久人人精品亚洲av| 日韩免费高清中文字幕av| 亚洲avbb在线观看| 99久久99久久久精品蜜桃| 欧美日本亚洲视频在线播放| 久久人妻福利社区极品人妻图片| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 国产不卡一卡二| 99国产综合亚洲精品| 欧美日韩一级在线毛片| 我的亚洲天堂| 又紧又爽又黄一区二区| 亚洲精品久久成人aⅴ小说| 国产成人精品在线电影| 免费搜索国产男女视频| 精品熟女少妇八av免费久了| 黄色a级毛片大全视频| 亚洲欧美日韩另类电影网站| 国产99久久九九免费精品| 97超级碰碰碰精品色视频在线观看| 黑人欧美特级aaaaaa片| 亚洲午夜理论影院| 午夜日韩欧美国产| 欧美激情极品国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 久久精品91无色码中文字幕| 欧美中文日本在线观看视频| 国产精品一区二区免费欧美| 黄网站色视频无遮挡免费观看| 精品人妻在线不人妻| 后天国语完整版免费观看| 看黄色毛片网站| 丁香欧美五月| 一级作爱视频免费观看| 久久久久国产精品人妻aⅴ院| 亚洲,欧美精品.| 美女 人体艺术 gogo| netflix在线观看网站| 国产又色又爽无遮挡免费看| bbb黄色大片| 老司机亚洲免费影院| 日韩 欧美 亚洲 中文字幕| 国产色视频综合| 在线观看舔阴道视频| avwww免费| 午夜精品在线福利| 欧美人与性动交α欧美软件| 亚洲熟妇熟女久久| 欧美另类亚洲清纯唯美| 国产精品影院久久| 国产色视频综合| 咕卡用的链子| 久久精品国产亚洲av香蕉五月| 后天国语完整版免费观看| av电影中文网址| 国产一卡二卡三卡精品| 狂野欧美激情性xxxx| 亚洲熟女毛片儿| 免费高清在线观看日韩| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 亚洲精品美女久久久久99蜜臀| 首页视频小说图片口味搜索| 黄色 视频免费看| 国产亚洲精品久久久久5区| 国产精品一区二区在线不卡| 国产真人三级小视频在线观看| 欧美日韩视频精品一区| 亚洲国产精品合色在线| 欧美激情极品国产一区二区三区| av天堂久久9| 婷婷精品国产亚洲av在线| 视频在线观看一区二区三区| 国产精品乱码一区二三区的特点 | 国产精品自产拍在线观看55亚洲| 高潮久久久久久久久久久不卡| 久久国产乱子伦精品免费另类| 亚洲人成电影观看| 亚洲国产精品合色在线| 亚洲性夜色夜夜综合| 色尼玛亚洲综合影院| 韩国精品一区二区三区| 中文欧美无线码| 亚洲自偷自拍图片 自拍| 免费在线观看完整版高清| 在线av久久热| 女人被躁到高潮嗷嗷叫费观| 午夜免费激情av| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 国产精品永久免费网站| 男女午夜视频在线观看| 女警被强在线播放| 在线观看www视频免费| 日本免费a在线| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 99香蕉大伊视频| 老司机深夜福利视频在线观看| 国产又爽黄色视频| 中文字幕人妻丝袜制服| 很黄的视频免费| 美女福利国产在线| 欧美av亚洲av综合av国产av| 亚洲精品中文字幕一二三四区| 久久久国产一区二区| 老司机亚洲免费影院| 国产又爽黄色视频| 中文字幕另类日韩欧美亚洲嫩草| 制服诱惑二区| 最新美女视频免费是黄的| 黄网站色视频无遮挡免费观看| 另类亚洲欧美激情| 一个人观看的视频www高清免费观看 | 精品无人区乱码1区二区| 变态另类成人亚洲欧美熟女 | e午夜精品久久久久久久| 国产亚洲av高清不卡| 91av网站免费观看| 在线国产一区二区在线| 十八禁网站免费在线| 高清av免费在线| 国产精品影院久久| 国产欧美日韩一区二区精品| 激情视频va一区二区三区| а√天堂www在线а√下载| 日日爽夜夜爽网站| 一二三四社区在线视频社区8| 两人在一起打扑克的视频| 在线国产一区二区在线| 日本vs欧美在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| av视频免费观看在线观看| 在线观看免费日韩欧美大片| 黄色毛片三级朝国网站| 日韩免费av在线播放| 天天躁夜夜躁狠狠躁躁| 精品免费久久久久久久清纯| 国产日韩一区二区三区精品不卡| 亚洲伊人色综图| 人成视频在线观看免费观看| av国产精品久久久久影院| 色在线成人网| 麻豆av在线久日| 黄色丝袜av网址大全| 成人av一区二区三区在线看| 久久中文字幕人妻熟女| 欧美黄色淫秽网站| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久 | 一区二区三区国产精品乱码| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区 | 身体一侧抽搐| 在线免费观看的www视频| 一区在线观看完整版| 亚洲精品久久午夜乱码| 亚洲精品一区av在线观看| 在线看a的网站| 精品午夜福利视频在线观看一区| 精品一区二区三卡| 亚洲精品中文字幕一二三四区| 婷婷丁香在线五月| 天堂俺去俺来也www色官网| 亚洲人成77777在线视频| 精品国产乱子伦一区二区三区| 黄色女人牲交| 国产精品 欧美亚洲| 无遮挡黄片免费观看| 午夜老司机福利片| 欧美日韩亚洲高清精品| 别揉我奶头~嗯~啊~动态视频| av免费在线观看网站| 免费在线观看完整版高清| av片东京热男人的天堂| 自线自在国产av| 一夜夜www| 最新美女视频免费是黄的| 在线观看免费视频网站a站| 国产成人一区二区三区免费视频网站| 久久精品91蜜桃| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品va在线观看不卡| 亚洲精品美女久久av网站| 成人18禁高潮啪啪吃奶动态图| 国产主播在线观看一区二区| 国产精品久久电影中文字幕| 欧美av亚洲av综合av国产av| 中文字幕色久视频| 成人三级黄色视频| 久久草成人影院| 91成人精品电影| 一进一出抽搐动态| 日韩精品青青久久久久久| 后天国语完整版免费观看| 亚洲人成网站在线播放欧美日韩| 超色免费av| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站 | 久99久视频精品免费| 久久久久久免费高清国产稀缺| 视频在线观看一区二区三区| 国产亚洲精品久久久久久毛片| 天堂中文最新版在线下载| 香蕉久久夜色| 午夜亚洲福利在线播放| 国产成+人综合+亚洲专区| 天天添夜夜摸| 久久精品亚洲精品国产色婷小说| 一夜夜www| 又黄又粗又硬又大视频| 51午夜福利影视在线观看| 日韩免费高清中文字幕av| 久久香蕉国产精品| 亚洲 欧美一区二区三区| a级片在线免费高清观看视频| 亚洲精品国产一区二区精华液| 精品久久久久久,| 日本一区二区免费在线视频| 亚洲国产欧美日韩在线播放| 真人一进一出gif抽搐免费| 亚洲精品一卡2卡三卡4卡5卡| 国产伦人伦偷精品视频| 午夜91福利影院| 亚洲精品美女久久av网站| 久久人人爽av亚洲精品天堂| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区 | 国产成人欧美在线观看| 亚洲精品一区av在线观看| 久99久视频精品免费| 激情视频va一区二区三区| 一级片免费观看大全| 在线十欧美十亚洲十日本专区| 一个人免费在线观看的高清视频| 男男h啪啪无遮挡| 一二三四社区在线视频社区8| 中文字幕最新亚洲高清| 亚洲视频免费观看视频| 级片在线观看| 久久精品国产清高在天天线| 国产人伦9x9x在线观看| a级片在线免费高清观看视频| 多毛熟女@视频| 久久亚洲精品不卡| 美女扒开内裤让男人捅视频| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 日本 av在线| 国产成人啪精品午夜网站| 欧美在线黄色| 精品日产1卡2卡| 在线视频色国产色| 可以免费在线观看a视频的电影网站| av天堂在线播放| aaaaa片日本免费| 成人国产一区最新在线观看| 三级毛片av免费| 国产精品综合久久久久久久免费 | 亚洲熟妇熟女久久| 免费久久久久久久精品成人欧美视频| 91国产中文字幕| 看免费av毛片| 亚洲少妇的诱惑av| 亚洲狠狠婷婷综合久久图片| 最新在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 久热这里只有精品99| 欧美成狂野欧美在线观看| 国产1区2区3区精品| 久久狼人影院| 久9热在线精品视频| 午夜福利影视在线免费观看| 丰满的人妻完整版| 97超级碰碰碰精品色视频在线观看| 天天躁夜夜躁狠狠躁躁| 美女 人体艺术 gogo| 91字幕亚洲| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 真人一进一出gif抽搐免费| 热re99久久精品国产66热6| 久久午夜亚洲精品久久| 国产成人欧美| 国产精品秋霞免费鲁丝片| 久久久国产欧美日韩av| 高清av免费在线| 一进一出好大好爽视频| 亚洲 欧美一区二区三区| 90打野战视频偷拍视频| 欧美精品啪啪一区二区三区| 18美女黄网站色大片免费观看| 天堂俺去俺来也www色官网| 久久久国产成人免费| 自线自在国产av| 亚洲一区高清亚洲精品| 日韩欧美一区视频在线观看| 黄片播放在线免费| 亚洲av电影在线进入| 无限看片的www在线观看| 欧美在线一区亚洲| 欧美日韩黄片免| 99久久国产精品久久久| 日本欧美视频一区| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久av网站| 亚洲午夜理论影院| 亚洲aⅴ乱码一区二区在线播放 | 18美女黄网站色大片免费观看| 久久精品91蜜桃| 国产极品粉嫩免费观看在线| 丁香欧美五月| 国产日韩一区二区三区精品不卡| 琪琪午夜伦伦电影理论片6080| 午夜91福利影院| cao死你这个sao货| 国产无遮挡羞羞视频在线观看| 99热国产这里只有精品6| 男女午夜视频在线观看| 国产精品 欧美亚洲| 1024香蕉在线观看| 丝袜美足系列| 久久久水蜜桃国产精品网| 久久伊人香网站| 丁香欧美五月| 在线av久久热| 国产精品国产av在线观看| 欧美最黄视频在线播放免费 | 国产精品一区二区免费欧美| 一个人免费在线观看的高清视频| 一边摸一边抽搐一进一出视频| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 99久久精品国产亚洲精品| 国产精品香港三级国产av潘金莲| 我的亚洲天堂| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩另类电影网站| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 黄色怎么调成土黄色| 久久中文字幕人妻熟女| 久热爱精品视频在线9| 一区在线观看完整版| 国产亚洲精品第一综合不卡| 国产激情久久老熟女| 深夜精品福利| 国产亚洲精品一区二区www| 亚洲黑人精品在线| 日韩中文字幕欧美一区二区| 免费在线观看亚洲国产| 男女做爰动态图高潮gif福利片 | 9热在线视频观看99| 大香蕉久久成人网| 涩涩av久久男人的天堂| 亚洲精品国产色婷婷电影| 国产一区二区激情短视频| 欧美 亚洲 国产 日韩一| 高清毛片免费观看视频网站 | 后天国语完整版免费观看| 男女午夜视频在线观看| 中文字幕最新亚洲高清| 嫩草影院精品99| av免费在线观看网站| 欧美日韩一级在线毛片| av在线播放免费不卡| 黄色怎么调成土黄色| 日韩欧美免费精品| 亚洲中文字幕日韩| 亚洲免费av在线视频| 午夜免费观看网址| 欧美日韩国产mv在线观看视频| 免费搜索国产男女视频| 三上悠亚av全集在线观看| 亚洲第一欧美日韩一区二区三区| 欧美激情 高清一区二区三区| 精品国产国语对白av| 母亲3免费完整高清在线观看| 午夜视频精品福利| 91麻豆av在线| 成年人免费黄色播放视频| 免费日韩欧美在线观看| 久久久国产成人免费| 巨乳人妻的诱惑在线观看| 91老司机精品| 51午夜福利影视在线观看| 欧美日韩黄片免| 亚洲成国产人片在线观看| 日韩欧美三级三区| 中文字幕精品免费在线观看视频| 两性夫妻黄色片| 波多野结衣av一区二区av| 两个人看的免费小视频| 亚洲九九香蕉| 超碰97精品在线观看| 欧美中文日本在线观看视频| 亚洲久久久国产精品| 超碰成人久久| 69精品国产乱码久久久| avwww免费| 国产精品野战在线观看 | 午夜视频精品福利| 午夜久久久在线观看| 久久久久久久午夜电影 | 女警被强在线播放| 两个人免费观看高清视频| 窝窝影院91人妻| 免费人成视频x8x8入口观看| 男人舔女人的私密视频| www.熟女人妻精品国产| 天天添夜夜摸| 精品国产超薄肉色丝袜足j| 欧美丝袜亚洲另类 | 国产成人影院久久av| 人人妻,人人澡人人爽秒播| 天堂影院成人在线观看| 大陆偷拍与自拍| 精品少妇一区二区三区视频日本电影| 国产国语露脸激情在线看| 国产蜜桃级精品一区二区三区| 亚洲中文日韩欧美视频| 国产精品久久视频播放| 女人高潮潮喷娇喘18禁视频| 精品第一国产精品| а√天堂www在线а√下载| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 日韩视频一区二区在线观看| 黄网站色视频无遮挡免费观看| 美女 人体艺术 gogo| 日韩 欧美 亚洲 中文字幕| 天堂影院成人在线观看| 一二三四社区在线视频社区8| 丝袜美足系列| 亚洲一码二码三码区别大吗| ponron亚洲| 国产av一区在线观看免费| 俄罗斯特黄特色一大片| 热re99久久精品国产66热6| 色综合欧美亚洲国产小说| 五月开心婷婷网| 久久国产精品影院| 首页视频小说图片口味搜索| 欧美日本中文国产一区发布| 九色亚洲精品在线播放| av片东京热男人的天堂| 天堂动漫精品| 亚洲av五月六月丁香网| 日韩一卡2卡3卡4卡2021年| 岛国视频午夜一区免费看| 亚洲国产精品合色在线| 三级毛片av免费| 欧美激情 高清一区二区三区| 99国产精品一区二区三区| 久久香蕉激情| 国产精品国产av在线观看| 一本大道久久a久久精品| 99精品久久久久人妻精品| 亚洲精品久久成人aⅴ小说| 黑人操中国人逼视频| 两性夫妻黄色片| 在线观看免费视频网站a站| 精品国产国语对白av| 久久午夜亚洲精品久久| 久久亚洲精品不卡| 搡老熟女国产l中国老女人| 久9热在线精品视频| 国产精品久久久人人做人人爽| 91在线观看av| 国产欧美日韩一区二区精品| 两个人免费观看高清视频| 成人永久免费在线观看视频| 在线观看免费午夜福利视频| 国产精品国产av在线观看| 超碰成人久久| 国产精品秋霞免费鲁丝片| 日本精品一区二区三区蜜桃| 老熟妇乱子伦视频在线观看| 亚洲成a人片在线一区二区| 日本黄色日本黄色录像| 麻豆成人av在线观看| 18禁黄网站禁片午夜丰满| 免费少妇av软件|