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

    我國(guó)北部冬麥區(qū)小麥區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化設(shè)計(jì)

    2020-07-21 08:29:10許乃銀郭利磊楊子光張笑晴楊曉妮
    作物學(xué)報(bào) 2020年8期
    關(guān)鍵詞:冬麥區(qū)遺傳力單點(diǎn)

    張 毅 許乃銀 郭利磊 楊子光 張笑晴 楊曉妮

    我國(guó)北部冬麥區(qū)小麥區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化設(shè)計(jì)

    張 毅1許乃銀2郭利磊1楊子光3張笑晴1楊曉妮2

    1全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心, 北京 100125;2江蘇省農(nóng)業(yè)科學(xué)院, 江蘇南京 210014;3洛陽(yáng)農(nóng)林科學(xué)院, 河南洛陽(yáng) 471022

    農(nóng)作物品種區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的合理配置有利于提高試驗(yàn)的成本效率和品種選擇效率。本研究分析了2010—2019年期間北部冬麥區(qū)小麥品種區(qū)域試驗(yàn)的重復(fù)次數(shù)和試點(diǎn)數(shù)量設(shè)置的合理性, 依據(jù)小麥品種試驗(yàn)的信噪比和遺傳力水平隨重復(fù)次數(shù)和試點(diǎn)數(shù)量的變化規(guī)律, 提出了重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化設(shè)計(jì)方案。結(jié)果表明: (1) 北部冬麥區(qū)小麥單點(diǎn)試驗(yàn)的遺傳力平均達(dá)到0.87, 需要的重復(fù)次數(shù)平均值僅為1.4, 說(shuō)明3次重復(fù)可以充分保證試驗(yàn)精確度的需求。(2)北部冬麥區(qū)水地組和旱地組小麥區(qū)域試驗(yàn)達(dá)到0.75的遺傳力水平時(shí), 需要的試點(diǎn)數(shù)量分別為11個(gè)和13個(gè), 目前有效試點(diǎn)數(shù)量分別約為11個(gè)和8個(gè), 分別達(dá)到0.75和0.60的遺傳力水平。(3)小麥品種區(qū)域試驗(yàn)結(jié)果對(duì)品種的審定和應(yīng)用十分重要, 而每年都可能有少數(shù)試驗(yàn)點(diǎn)因?yàn)楦鞣N異常情況而報(bào)廢, 為保證試驗(yàn)結(jié)果的可靠性, 可按= 0.75的水平需求安排試驗(yàn)點(diǎn)數(shù)量和重復(fù)次數(shù), 即重復(fù)次數(shù)可保持當(dāng)前的3次; 水地組的試點(diǎn)數(shù)量可保持在11個(gè)左右; 旱地組可將試點(diǎn)增加到13個(gè); 如要將遺傳力提高到0.80的水平, 則需約16個(gè)試點(diǎn)。

    小麥(L.); 區(qū)域試驗(yàn); 遺傳力; 信噪比; 重復(fù)次數(shù); 試點(diǎn)數(shù)量

    小麥(L.)是我國(guó)第三大糧食作物, 種植面積僅次于玉米和水稻。北方冬麥區(qū)主要包括河北大部、河南、山東、山西南部、江蘇和安徽北部、陜西和甘肅部分地區(qū), 是我國(guó)小麥主要產(chǎn)區(qū), 種植面積約占全國(guó)小麥總面積的70%, 總產(chǎn)占75%左右[1]。北部冬麥區(qū)主要包括遼寧南端的營(yíng)口和大連市; 河北省境內(nèi)長(zhǎng)城以南的廊坊、保定、滄州、唐山和秦皇島市全部; 京、津兩市全部; 山西省朔州以南陽(yáng)泉、太原、晉中、長(zhǎng)治和呂梁等市全部和臨汾市北部地區(qū); 陜西省延安市全部、榆林長(zhǎng)城以南大部, 咸陽(yáng)、寶雞和銅川市部分縣; 甘肅省隴東慶陽(yáng)市全部和平?jīng)鍪械牟糠挚h[2]。北部冬麥區(qū)作為北方冬麥區(qū)的重要產(chǎn)區(qū), 提供了其中22%以上的產(chǎn)量[3]。加強(qiáng)適應(yīng)于該區(qū)域生態(tài)條件的小麥新品種選育和推廣應(yīng)用對(duì)保障小麥生產(chǎn)水平和我國(guó)糧食安全具有重要意義。小麥品種區(qū)域試驗(yàn)在多環(huán)境條件下對(duì)新品種的豐產(chǎn)性、穩(wěn)產(chǎn)性、適應(yīng)性、抗逆性、品質(zhì)等性狀進(jìn)行測(cè)試與綜合評(píng)價(jià), 對(duì)促進(jìn)我國(guó)小麥品種更新和保障小麥生產(chǎn)水平起著重要的作用[4]。區(qū)域試驗(yàn)的可靠性主要決定于試驗(yàn)效率, 而試驗(yàn)效率又決定于試驗(yàn)遺傳力[5-6]。遺傳力是遺傳方差占表型方差的比率[7], 區(qū)域試驗(yàn)的遺傳力越高, 對(duì)品種的選擇效率就越高[8]。遺傳力隨著試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的增加速率呈先快后慢的“邊際收益遞減”規(guī)律[6]。因此, 如何合理設(shè)置試驗(yàn)試點(diǎn)數(shù)量和重復(fù)次數(shù), 既可保持可靠的遺傳力水平和試驗(yàn)效率, 又可合理利用有限的人力、物力和財(cái)力資源, 是優(yōu)化區(qū)域試驗(yàn)設(shè)置方案中亟待解決的問(wèn)題。

    嚴(yán)威凱等[8]利用區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量與遺傳力的函數(shù)關(guān)系, 估算出燕麥品種試驗(yàn)在特定誤差水平下不同遺傳力水平需要的試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量, 并指出通過(guò)增加試點(diǎn)數(shù)和重復(fù)次數(shù)可以快速提高遺傳力的上限為75%, 其后遺傳力提高的速率急劇降低, 故認(rèn)為試驗(yàn)遺傳力達(dá)到0.75水平時(shí)需要的重復(fù)次數(shù)和試點(diǎn)數(shù)量是較適宜的數(shù)量。許乃銀等[5]對(duì)國(guó)家棉花區(qū)試的遺傳力水平及其隨著試點(diǎn)數(shù)量增加的規(guī)律進(jìn)行了定量化分析, 并提出針對(duì)各主產(chǎn)棉區(qū)的區(qū)域試驗(yàn)設(shè)置優(yōu)化方案。關(guān)于小麥品種區(qū)域試驗(yàn)適宜重復(fù)次數(shù)和試點(diǎn)數(shù)量的相關(guān)研究尚未見(jiàn)報(bào)道。本研究的目的是綜合利用重復(fù)次數(shù)和試點(diǎn)數(shù)量的估計(jì)方法, 全面分析和評(píng)價(jià)2010—2019年期間我國(guó)北部冬麥區(qū)國(guó)家小麥品種試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量設(shè)置的科學(xué)性和優(yōu)化方案, 為國(guó)家小麥區(qū)域試驗(yàn)的優(yōu)化資源配置提供依據(jù), 并為其他作物品種區(qū)域試驗(yàn)的優(yōu)化設(shè)置提供參考。

    1 材料與方法

    1.1 數(shù)據(jù)來(lái)源

    2010—2019年期間北部冬麥區(qū)水地組和旱地組小麥品種試驗(yàn)設(shè)置的試點(diǎn)數(shù)量和參試品種數(shù)詳見(jiàn)表1。所有試驗(yàn)均采用隨機(jī)區(qū)組排列, 重復(fù)3次, 小區(qū)面積12~14 m2。2009—2010年度(簡(jiǎn)稱為“2010年”, 下同)水地組區(qū)域試驗(yàn)以“京冬8號(hào)”為對(duì)照, 以“中麥175”輔助對(duì)照品種。其余年份水地組和旱地組小麥品種試驗(yàn)分別以中麥175和“長(zhǎng)6878”為對(duì)照品種。水地組常年在北京市的昌平和順義; 天津市的武清和寶坻; 河北省的保定、遵化、灤南和固安; 山西省的太原、屯留、介休等地共設(shè)置約11個(gè)試點(diǎn)(2015年前在新疆設(shè)有阿拉爾試點(diǎn))。旱地組常年在山西省的高平、長(zhǎng)子、長(zhǎng)治和平?jīng)? 甘肅省的莊浪、定西、鎮(zhèn)原和慶陽(yáng); 寧夏的固原、隆德、彭陽(yáng)等地共設(shè)置約11個(gè)試點(diǎn)。田間栽培管理按當(dāng)?shù)亓?xí)慣和區(qū)域試驗(yàn)方案要求進(jìn)行。本研究分別采用2010—2019年期間北部冬麥區(qū)小麥品種試驗(yàn)的196個(gè)單年單點(diǎn)試驗(yàn)(其中, 水地組112個(gè), 旱地組84個(gè))和20組一年多點(diǎn)試驗(yàn)的產(chǎn)量數(shù)據(jù)分析試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量與遺傳力的關(guān)系, 估算在特定遺傳力水平下對(duì)重復(fù)次數(shù)和試點(diǎn)數(shù)量的需求。

    1.2 統(tǒng)計(jì)分析

    采用數(shù)據(jù)分析和管理系統(tǒng)GGEbiplot Pattern Explorer[9]對(duì)2010—2019年期間北部冬麥區(qū)小麥品種試驗(yàn)中196個(gè)單點(diǎn)試驗(yàn)分別進(jìn)行方差分析, 得出遺傳方差和誤差方差, 再依據(jù)Yan等[6,8]和許乃銀等[5]介紹的方法分別計(jì)算單點(diǎn)試驗(yàn)的信噪比、遺傳力及遺傳力為0.75時(shí)需要的試驗(yàn)重復(fù)次數(shù), 并根據(jù)196次單點(diǎn)試驗(yàn)在不同重復(fù)次數(shù)需求范圍內(nèi)的累積頻數(shù)(cumulative frequency, C)分布規(guī)律, 初步分析小麥品種區(qū)域試驗(yàn)的適宜重復(fù)次數(shù)。其次, 對(duì)20組一年多點(diǎn)試驗(yàn)分別進(jìn)行方差分析, 得出遺傳方差、品種與試點(diǎn)互作方差和誤差方差, 計(jì)算一年多點(diǎn)試驗(yàn)的信噪比、遺傳力及遺傳力為0.75時(shí)需要的試點(diǎn)數(shù)量,分別估計(jì)水地組和旱地組試驗(yàn)對(duì)試點(diǎn)數(shù)量的需求。然后, 依據(jù)水地組和旱地組小麥品種試驗(yàn)的多次單點(diǎn)試驗(yàn)平均信噪比(r)和一年多點(diǎn)試驗(yàn)平均信噪比(e)與重復(fù)次數(shù)、試點(diǎn)數(shù)及遺傳力的關(guān)系式, 對(duì)不同重復(fù)次數(shù)和試點(diǎn)數(shù)量條件下的遺傳力水平進(jìn)行定量分析。最后, 依據(jù)遺傳力與信噪比、重復(fù)次數(shù)和試點(diǎn)數(shù)的關(guān)系式, 在設(shè)定區(qū)間內(nèi)擬合遺傳力增量與重復(fù)次數(shù)和試點(diǎn)數(shù)量增量的響應(yīng)曲線, 從而提出小麥品種試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量設(shè)置的優(yōu)化方案。單年單點(diǎn)和一年多點(diǎn)區(qū)域試驗(yàn)的遺傳力及其增量、信噪比、重復(fù)次數(shù)和試點(diǎn)數(shù)量的計(jì)算公式如下[5-6]:

    單年單點(diǎn)試驗(yàn)的遺傳力=g/(g+e/r) (1)

    單年單點(diǎn)試驗(yàn)的信噪比r=e/g(2)

    單年單點(diǎn)試驗(yàn)所需的重復(fù)數(shù)r= (e/g)×/ (1?) =r×/(1?) (3)

    單年單點(diǎn)試驗(yàn)遺傳力與重復(fù)次數(shù)及信噪比的關(guān)系式=r/(r+r) (4)

    單年單點(diǎn)試驗(yàn)遺傳力增量D= [r/(r+r)][(r?Dr)/((r?Dr)+r)] (5)

    一年多點(diǎn)試驗(yàn)的遺傳力=g/(g+ge/e+e/er) (6)

    一年多點(diǎn)試驗(yàn)的信噪比Q= (ge+e/e)/g(7)

    一年多點(diǎn)試驗(yàn)所需的試點(diǎn)數(shù)e= [(ge+e/e)/g] ×/(1?) =e×/(1?) (8)

    一年多點(diǎn)試驗(yàn)遺傳力與試點(diǎn)數(shù)及信噪比的關(guān)系式=e/(e+e) (9)

    一年多點(diǎn)試驗(yàn)遺傳力增量D= [e/(e+e)][(e?De)/((e?De)+e)] (10)

    式中,為遺傳力,g為遺傳方差,e為誤差方差,r為重復(fù)數(shù),r為單點(diǎn)試驗(yàn)信噪比,r為單點(diǎn)試驗(yàn)需要的重復(fù)數(shù),ge為品種與試點(diǎn)互作方差,e為試點(diǎn)數(shù)量,e為一年多點(diǎn)試驗(yàn)的信噪比,e為一年多點(diǎn)試驗(yàn)需要的試點(diǎn)數(shù)量,D為試驗(yàn)遺傳力增量,Dr和De分別為單位重復(fù)數(shù)和試點(diǎn)數(shù),(r?Dr)和(e?De)分別為在r和e基礎(chǔ)上減少一個(gè)單位的重復(fù)數(shù)和試點(diǎn)數(shù)。

    表1 2010–2019年我國(guó)北部冬麥區(qū)國(guó)家小麥品種區(qū)域試驗(yàn)組數(shù)、品種數(shù)和產(chǎn)量統(tǒng)計(jì)表

    2 結(jié)果與分析

    2.1 單年單點(diǎn)小麥區(qū)域試驗(yàn)的重復(fù)次數(shù)分析

    表2表明: (1)在北部冬麥區(qū)共196次單點(diǎn)試驗(yàn)中, 在保證遺傳力達(dá)到0.75水平下需要1次、2次和3次重復(fù)以內(nèi)的試驗(yàn)累積頻率分別約為69%、83%和88%, 需要3~4次重復(fù)、4~5次重復(fù)和5次以上重復(fù)的試驗(yàn)比例分別約為4%、2%和7%。其中, 水地組和旱地組需要3次重復(fù)以內(nèi)的試驗(yàn)累積頻率分別約為87%和89%, 差異不大。可見(jiàn), 采用3次重復(fù)可以保證近90%的單點(diǎn)試驗(yàn)遺傳力達(dá)到0.75的水平, 僅有少數(shù)試點(diǎn)因試驗(yàn)誤差偏大等原因?qū)е逻z傳力偏低; (2)北部冬麥區(qū)共196次單點(diǎn)試驗(yàn)需要的重復(fù)次數(shù)總均值僅為1.4次, 其中, 水地組和旱地組需要的重復(fù)次數(shù)平均都在2次以下, 說(shuō)明北部冬麥區(qū)小麥品種試驗(yàn)采用3次重復(fù)可以充分保證試驗(yàn)遺傳力和試驗(yàn)效率符合試驗(yàn)質(zhì)量要求。

    表2 2010–2019年北部冬麥區(qū)單年單點(diǎn)小麥品種試驗(yàn)在遺傳力為0.75時(shí)所需要重復(fù)數(shù)的次數(shù)分布

    重復(fù)數(shù)范圍是指在=0.75時(shí)所需要的試驗(yàn)重復(fù)次數(shù)的分布范圍。北部冬麥區(qū)包括水地組和旱地組。N、(%)、C(%)、r和r分別為次數(shù)、頻數(shù)、累積頻數(shù)、需要重復(fù)數(shù)和信噪比。

    Range ofrmeans the replicates needed to achieve 0.75 of heritability in half-open intervals. The Northern winter wheat region covers Irrigated region and Rainfed region. N,(%), C(%),r,r, and(%) stand for number of times, frequency, cumulative frequency, replicate needed and noise-signal ratio in each interval, respectively.

    2.2 小麥品種試驗(yàn)適宜重復(fù)數(shù)和試點(diǎn)數(shù)的定量分析

    依據(jù)2010—2019年水地組、旱地組及北部冬麥區(qū)小麥品種試驗(yàn)的單點(diǎn)試驗(yàn)平均信噪比(r)和一年多點(diǎn)試驗(yàn)平均信噪比(e)與重復(fù)次數(shù)或試點(diǎn)數(shù)及遺傳力的關(guān)系式, 在不同重復(fù)次數(shù)和試點(diǎn)數(shù)量條件下對(duì)遺傳力水平的定量分析(圖1和表3)表明: (1)試驗(yàn)遺傳力隨著重復(fù)次數(shù)和試點(diǎn)數(shù)的增加, 表現(xiàn)為先快后慢的曲線關(guān)系。隨著重復(fù)次數(shù)和試點(diǎn)數(shù)的增加, 遺傳力在0.75以下時(shí)提升速度非??? 為遺傳力的“快速增長(zhǎng)期”; 在0.75和0.90期間提升速率明顯減緩, 但仍有顯著的增量, 為遺傳力的“緩慢增長(zhǎng)期”; 達(dá)0.90之后, 再增加重復(fù)次數(shù)或試點(diǎn)數(shù)量對(duì)遺傳力的提高作用微小, 是遺傳力的“增長(zhǎng)停滯期”。(2)在遺傳力為0.75水平下, 水地組、旱地組及北部冬麥區(qū)小麥品種試驗(yàn)需要的重復(fù)次數(shù)都在2次以下; 重復(fù)次數(shù)為3次的情況下, 水地組、旱地組及北部冬麥區(qū)小麥品種試驗(yàn)的遺傳力分別達(dá)到了約0.85、0.88和0.87的水平, 已經(jīng)充分保證了試驗(yàn)的精度要求; 如需將遺傳力提高在0.90水平下, 則需要約4~5次重復(fù), 其后繼續(xù)增加重復(fù)次數(shù)對(duì)提高遺傳力沒(méi)有意義。(3)在遺傳力為0.75水平下, 水地組、旱地組及北部冬麥區(qū)小麥品種試驗(yàn)需要的試點(diǎn)數(shù)量分別約為11個(gè)、13個(gè)和12個(gè); 水地組、旱地組及北部冬麥區(qū)小麥品種試驗(yàn)平均匯總的試點(diǎn)數(shù)分別約為11個(gè)、9個(gè)和10個(gè), 其試驗(yàn)遺傳力分別達(dá)到了約0.75、0.68和0.71的水平, 水地組已經(jīng)符合小麥品種試驗(yàn)的精確度要求, 而旱地組需要增加4個(gè)試點(diǎn)才可以達(dá)到試驗(yàn)的精確度要求; 在條件許可的情況下, 將旱地組試點(diǎn)數(shù)增加到13個(gè), 可將整體試驗(yàn)的遺傳力水平提升到0.75; 如需將遺傳力提高到0.85水平, 則水地組和旱地組分別需要約21個(gè)和24個(gè)試點(diǎn); 如需將遺傳力提高到0.90水平, 則水地組和旱地組分別需要約33個(gè)和38個(gè)試點(diǎn), 其后繼續(xù)增加試點(diǎn)數(shù)沒(méi)有意義。

    圖1 北部冬麥區(qū)小麥品種試驗(yàn)的遺傳力(H)及其所需重復(fù)(Nr)或試點(diǎn)數(shù)量(Ne)的關(guān)系

    Irrigated和Rainfed分別代表北部冬麥區(qū)的水地組和旱地組小麥品種試驗(yàn);r和e分別代表重復(fù)數(shù)信噪比和試點(diǎn)數(shù)信噪比。

    Irrigated and Rainfed stand for the irrigated group and rainfed group in winter wheat variety trials of Northern winter wheat region.randestand for the noise-signal ratio of replicate and test location number, respectively.

    表3 北部冬麥區(qū)小麥品種試驗(yàn)在不同遺傳力水平下需要的重復(fù)數(shù)和試點(diǎn)數(shù)量定量分析

    重復(fù)的信噪比(r)為單年單點(diǎn)試驗(yàn)的平均值; 試點(diǎn)的信噪比(e)一年多點(diǎn)試驗(yàn)校正平均值; 合格數(shù)量為2010–2019年期間北部冬麥區(qū)水地組和旱地組實(shí)際匯總的重復(fù)或試點(diǎn)數(shù)量。

    The noise-signal ratio of replicate (r) is the mean of all single trials in each group, while locationeis the mean of all one-year multi-locational trials in the group.Nis the mean of replicates or test locations summarized in the annual wheat cultivar trial report of the irrigated group and rainfed group trials in Northern winter wheat region (NWWR) from 2010 to 2019.

    2.3 小麥品種試驗(yàn)中重復(fù)數(shù)和試點(diǎn)數(shù)對(duì)遺傳力增量的影響

    依據(jù)2010—2019年北部冬麥區(qū)水地組和旱地組小麥品種試驗(yàn)的單點(diǎn)試驗(yàn)平均信噪比(r)和一年多點(diǎn)試驗(yàn)平均信噪比(e)與重復(fù)次數(shù)或試點(diǎn)數(shù)及遺傳力的關(guān)系式, 計(jì)算在設(shè)定區(qū)間內(nèi)與單位增量的重復(fù)次數(shù)和試點(diǎn)數(shù)量相對(duì)應(yīng)的遺傳力增量。圖2表明, (1)水地組和旱地組小麥品種試驗(yàn)重復(fù)數(shù)和試點(diǎn)數(shù)增量對(duì)遺傳力增量的影響趨勢(shì)基本一致, 可合并分析。(2)重復(fù)次數(shù)以0.2為單位, 在區(qū)間[0,1]中每增加一個(gè)單位, 遺傳力依次提升0.30、0.16、0.10、0.07和0.05, 當(dāng)重復(fù)次數(shù)為1時(shí)遺傳力達(dá)到0.68; 在區(qū)間[1,5]中每增加一次重復(fù), 遺傳力依次提升0.13、0.05、0.03和0.02; 其后繼續(xù)增加重復(fù)數(shù), 遺傳力基本不再提升。重復(fù)數(shù)從1次增加到2次、3次和4次, 遺傳力依次從0.68提升到0.81、0.86、0.89和0.91, 可見(jiàn), 第2次重復(fù)的遺傳力增量較大, 第3次重復(fù)的遺傳力增加已經(jīng)很小, 后續(xù)重復(fù)數(shù)增加對(duì)遺傳力的提升效應(yīng)十分微小。(3)試點(diǎn)數(shù)以1為單位, 在區(qū)間[1,5]中每增加1個(gè)試點(diǎn), 遺傳力增量依次為0.20、0.13、0.10、0.07和0.06, 從0.20提高到0.60; 其后繼續(xù)增加試點(diǎn)數(shù), 遺傳力的提升速度減緩, 在試點(diǎn)數(shù)區(qū)間[6,10]中每增加1個(gè)試點(diǎn)約提高遺傳力0.02, 在試點(diǎn)數(shù)區(qū)間[11,15]中每增加1個(gè)試點(diǎn)約提高遺傳力0.01, 其后試點(diǎn)數(shù)增加對(duì)遺傳力的提升效應(yīng)較小。試點(diǎn)數(shù)從5個(gè)增加到10個(gè)、15個(gè)和20個(gè), 遺傳力依次從0.56提升到0.72、0.79和0.84, 后續(xù)試點(diǎn)數(shù)增加對(duì)遺傳力的提升效應(yīng)基本上沒(méi)有意義。

    圖2 北部冬麥區(qū)小麥品種試驗(yàn)的遺傳力增量(ΔH)與單位重復(fù)數(shù)或試點(diǎn)數(shù)的關(guān)系

    3 討論

    3.1 小麥品種區(qū)域試驗(yàn)效率的評(píng)價(jià)指標(biāo)探討

    農(nóng)作物多環(huán)境品種試驗(yàn)(品種區(qū)域試驗(yàn))是客觀評(píng)價(jià)農(nóng)作物品種的豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性及生產(chǎn)應(yīng)用價(jià)值必要環(huán)節(jié)[10]。作物區(qū)域試驗(yàn)的科學(xué)性、可靠性和試驗(yàn)效率主要決定于其對(duì)品種間差異的鑒別能力, 而鑒別力又決定于試驗(yàn)的精確度[11]。試驗(yàn)精確度在區(qū)域試驗(yàn)質(zhì)量評(píng)價(jià)中一直受到廣泛關(guān)注[12-16]。試驗(yàn)的精確度主要與試驗(yàn)誤差有關(guān), 試驗(yàn)誤差越大, 則精確度越低。區(qū)域試驗(yàn)中通常用試驗(yàn)誤差變異系數(shù)(CV, coefficient of error variation)來(lái)表示試驗(yàn)的精確度, 并將CV大于15%的單年單點(diǎn)試驗(yàn)數(shù)據(jù)在匯總報(bào)告中剔除。作為更直接的試驗(yàn)對(duì)品種鑒別能力的指標(biāo), 孔繁玲等[11]將品種間多重比較時(shí)的最小顯著差數(shù)(least significant difference, LSD)與總體均值的比率稱為相對(duì)最小顯著差數(shù)(relative least significant difference, RLSD)作為試驗(yàn)精確度的指標(biāo)[17]。然而, 一年多點(diǎn)品種區(qū)域試驗(yàn)的精確度和有效性不僅受到遺傳效應(yīng)和誤差大小的影響, 還受到品種與環(huán)境互作效應(yīng)等因素的影響[18-20]。一組由多個(gè)試驗(yàn)誤差變異都很小的單點(diǎn)試驗(yàn)組成一年多點(diǎn)品種區(qū)域試驗(yàn)并不一定是有效的試驗(yàn)。例如, 2017年北部冬麥區(qū)旱地組小麥區(qū)域試驗(yàn)中全部9個(gè)單點(diǎn)試驗(yàn)中的CV值都在8%以下, 多點(diǎn)聯(lián)合方差分析的CV值僅為3.8%, 但其試驗(yàn)遺傳力卻只有0.30。從CV值來(lái)看這是一個(gè)精確度很好的試驗(yàn), 而從遺傳力來(lái)看又是一個(gè)“壞試驗(yàn)”, 因?yàn)樗鼘?duì)品種的鑒別能力太差, 無(wú)法分辨品種間的真實(shí)遺傳差異。品種與環(huán)境互作方差大而遺傳方差相對(duì)很小是產(chǎn)生這種現(xiàn)象的主要原因。當(dāng)遺傳力接近于0時(shí), 無(wú)論品種間表型差異有多大, 基于表型變異的品種選擇都是完全無(wú)效的[8]。遺傳力表達(dá)了試驗(yàn)對(duì)品種間遺傳差異的鑒別能力和對(duì)品種的選擇效率, 也就是體現(xiàn)了品種區(qū)域試驗(yàn)的效率。在農(nóng)作物品種多環(huán)境試驗(yàn)質(zhì)量評(píng)價(jià)體系中除了用CV表達(dá)試驗(yàn)的誤差控制水平外, 還應(yīng)當(dāng)將試驗(yàn)遺傳力作為試驗(yàn)有效性的重要評(píng)價(jià)指標(biāo)[5]。Yan等[6]研究發(fā)現(xiàn)= 0.75是比較適當(dāng)?shù)脑囼?yàn)評(píng)價(jià)標(biāo)準(zhǔn)。許乃銀等[5]研究證實(shí)了遺傳力保持0.75的水平是我國(guó)棉花品種區(qū)域試驗(yàn)質(zhì)量評(píng)價(jià)比較適合的標(biāo)準(zhǔn), 但也需要依據(jù)具體品種試驗(yàn)發(fā)展水平確定其適宜的目標(biāo)遺傳力水平。本研究發(fā)現(xiàn)我國(guó)北部冬麥區(qū)小麥品種試驗(yàn)水地組和旱地組在正常年份的試驗(yàn)遺傳力平均分別為0.75和0.68, 同時(shí)也考慮到小麥品種區(qū)域試驗(yàn)對(duì)品種推薦審定和應(yīng)用十分重要, 每年都可能有部分試驗(yàn)點(diǎn)因?yàn)楦鞣N異常原因而報(bào)廢, 為充分保證試驗(yàn)匯總結(jié)果的可靠性, 可按遺傳力為0.85的水平需求安排試驗(yàn)點(diǎn)數(shù)量和重復(fù)次數(shù)。

    3.2 小麥品種區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化策略

    小麥品種區(qū)域試驗(yàn)從試驗(yàn)設(shè)計(jì)、田間管理和試驗(yàn)統(tǒng)計(jì)分析等各項(xiàng)工作都是為了提高小麥品種試驗(yàn)的遺傳力, 從而提高試驗(yàn)的品種選擇效率[5]。品種多點(diǎn)試驗(yàn)的遺傳力受到試驗(yàn)誤差、品種與試點(diǎn)互作效應(yīng)、試點(diǎn)數(shù)量和試驗(yàn)重復(fù)次數(shù)等因素的綜合影響。因此, 加強(qiáng)試驗(yàn)田間管理以減少試驗(yàn)誤差、開(kāi)展適當(dāng)?shù)钠贩N生態(tài)區(qū)劃分以減少品種與環(huán)境互作效應(yīng)、適當(dāng)增加試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量區(qū)域試驗(yàn)等措施有利于提高試驗(yàn)遺傳力和試驗(yàn)效率。試驗(yàn)遺傳力隨著重復(fù)次數(shù)和試點(diǎn)數(shù)量增加而提高的速率符合先快后慢的“邊際效應(yīng)遞減”規(guī)律, 在達(dá)到適當(dāng)?shù)臄?shù)量后再增加重復(fù)次數(shù)和試點(diǎn)數(shù)量對(duì)遺傳力的增量貢獻(xiàn)很小, 試驗(yàn)效率提高不大。同時(shí), 重復(fù)次數(shù)和試點(diǎn)數(shù)量的增加也意味著試驗(yàn)投入的人力、物力和財(cái)力成本的增加, 成本投入對(duì)試驗(yàn)效率的提高符合報(bào)酬遞減規(guī)律。因此, 品種區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化就是依據(jù)遺傳力隨著重復(fù)次數(shù)和試點(diǎn)數(shù)量增加的遞增規(guī)律, 設(shè)計(jì)適當(dāng)?shù)闹貜?fù)次數(shù)和試點(diǎn)數(shù)量在充分保證試驗(yàn)的可靠性和科學(xué)性的前提下又能節(jié)省試驗(yàn)成本, 減少不必要的資源浪費(fèi)。嚴(yán)威凱等[6]研究指出遺傳力達(dá)到0.75的水平是試驗(yàn)成本和試驗(yàn)效率的產(chǎn)投比“拐點(diǎn)”, 遺傳力大于0.75時(shí)的試點(diǎn)數(shù)量增加對(duì)遺傳力增量效應(yīng)將快速衰退。本研究也證實(shí)了試驗(yàn)遺傳力隨著重復(fù)次數(shù)和試點(diǎn)數(shù)的增加表現(xiàn)為先快后慢的曲線關(guān)系, 遺傳力在0.75以下時(shí)表現(xiàn)為“快速增長(zhǎng)期”, 在0.75和0.90期間表現(xiàn)為“緩慢增長(zhǎng)期”, 在0.90之后表現(xiàn)為“增長(zhǎng)停滯期”。本研究發(fā)現(xiàn)我國(guó)北部冬麥區(qū)小麥品種試驗(yàn)第2次重復(fù)的遺傳力增量較大, 第3次重復(fù)的遺傳力增加已經(jīng)很小, 后續(xù)重復(fù)數(shù)增加對(duì)遺傳力的提升效應(yīng)十分微小; 試點(diǎn)數(shù)從5個(gè)增加到10個(gè)、15個(gè)和20個(gè), 遺傳力依次從0.56提升到0.72、0.79和0.84, 后續(xù)試點(diǎn)數(shù)增加對(duì)遺傳力的提升效應(yīng)基本上沒(méi)有意義。因此, 北部冬麥區(qū)小麥品種試驗(yàn)可保持當(dāng)前的3次重復(fù); 水地組的試點(diǎn)數(shù)量可保持在11個(gè)左右; 旱地組可將試點(diǎn)數(shù)增加到13個(gè), 使遺傳力達(dá)到0.75; 在條件許可的前提下, 旱地組和水地組的試點(diǎn)數(shù)可分別增加到15個(gè)和17個(gè), 使遺傳力達(dá)到0.80的水平。

    4 結(jié)論

    我國(guó)北部冬麥區(qū)小麥品種試驗(yàn)采用3次重復(fù)可以充分保證試驗(yàn)遺傳力和試驗(yàn)效率符合試驗(yàn)質(zhì)量要求, 目前設(shè)置的試點(diǎn)數(shù)量也可以充分滿足區(qū)域試驗(yàn)的精確度要求。針對(duì)小麥品種試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的優(yōu)化方案包括,重復(fù)次數(shù)仍可設(shè)計(jì)為3次; 水地組可保持當(dāng)前11個(gè)左右的試點(diǎn)數(shù)量, 旱地組和在國(guó)家財(cái)力、物力和人力等條件許可的情況下, 可以將試點(diǎn)數(shù)量增加到13個(gè)左右, 從而使北部冬麥區(qū)小麥品種試驗(yàn)遺傳力達(dá)到到0.75的水平; 如需將遺傳力提高在0.80水平, 則需要約16個(gè)試點(diǎn); 增加的試點(diǎn)數(shù)量可與增加重復(fù)次數(shù)相比較, 優(yōu)先選擇增加小區(qū)數(shù)量較少的方案, 有利于節(jié)約試驗(yàn)成本。

    [1] 胡學(xué)旭, 孫麗娟, 周桂英, 吳麗娜, 陸偉, 李為喜, 王爽, 楊秀蘭, 宋敬可, 王步軍. 2000–2015年北部、黃淮冬麥區(qū)國(guó)家區(qū)試品種的品質(zhì)特征. 作物學(xué)報(bào), 2017, 43: 501–509. Hu X X, Sun L J, Zhou G Y, Wu L N, Lu W, Li W X, Wang S, Yang X L, Song J K, Wang B JQuality characteristics of winter wheat varieties tested in national trials in northern region and Yellow-Huai River valley winter wheat zone from 2000 to 2015., 2017, 43: 501–509 (in Chinese with English abstract).

    [2] 趙廣才. 小麥優(yōu)質(zhì)高產(chǎn)栽培理論與技術(shù). 北京: 中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社, 2018. Zhao G C. Good Quality and High Yield Cultivation of Wheat: Theory and Technology. Beijing: China Agricultural Science and Technology Press, 2018.

    [3] 孟繁圓, 馮利平, 張豐瑤, 張祎, 伍露, 王春雷, 閆錦濤, 彭明喜, 莫志鴻, 余衛(wèi)東. 北部冬麥區(qū)冬小麥越冬凍害時(shí)空變化特征. 作物學(xué)報(bào), 2019, 45: 1576–1585. Meng F Y, Feng L P, Zhang F Y, Zhang W, Wu L, Wang C L, Yan J T, Peng M X, Mo Z H, Yu W D. Temporal and spatial variations of winter wheat freezing injury in northern winter wheat region., 2019, 45: 1576–1585 (in Chinese with English abstract).

    [4] 胡學(xué)旭, 孫麗娟, 周桂英, 吳麗娜, 陸偉, 李為喜, 王爽, 楊秀蘭, 宋敬可, 王步軍. 2000–2015年國(guó)家黃淮和北部冬麥區(qū)域試驗(yàn)品種品質(zhì)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2016, 49: 4677–4686. Hu X X, Sun L J, Zhou G Y, Wu L N, Lu W, Li W X, Wang S, Yang X L, Song J K, Wang B JQuality variation of national tested varieties in northern winter wheat region and Yellow-Huai River valley winter wheat region from 2000 to 2015., 2016, 49: 4677–4686 (in Chinese with English abstract).

    [5] 許乃銀, 金石橋, 李健. 我國(guó)棉花品種區(qū)域試驗(yàn)重復(fù)次數(shù)和試點(diǎn)數(shù)量的設(shè)計(jì). 作物學(xué)報(bào), 2016, 42: 43–50.Xu N Y, Jin S Q, Li J. Design of test location number and replicate frequency in the regional cotton variety trials in China., 2016, 42: 43–50 (in Chinese with English abstract)

    [6] Yan W K, Frégeau-Reid J, Martin R, Pageau D, Mitchell-Fetch J. How many test locations and replications are needed in crop variety trials for a target region?, 2015, 202: 361–372.

    [7] 翟虎渠. 應(yīng)用數(shù)量遺傳. 北京: 中國(guó)農(nóng)業(yè)出版社, 2001. Zhai H Q. Applied Quantitative Genetics. Beijing: Chinese Agriculture Press, 2001 (in Chinese).

    [8] Yan W K. Crop Variety Trials: Data Management and Analysis. Oxford, UK: Wiley-Blackwell, 2014.

    [9] Yan W K. GGEbiplot-A windows application for graphical analysis of multienvironment trial data and other types of two-way data., 2001, 93: 1111–1118.

    [10] 許乃銀, 李健. 我國(guó)棉花品種區(qū)域試驗(yàn)精確度的演變分析. 棉花學(xué)報(bào), 2016, 28: 34–41 Xu N Y, Li J. Evolution analysis of the precision of regional cotton variety trials in China., 2016, 28: 34–41 (in Chinese with English abstract).

    [11] 孔繁玲, 張群遠(yuǎn), 楊付新, 郭恒敏. 棉花品種區(qū)域試驗(yàn)的精確度探討. 作物學(xué)報(bào), 1998, 24: 601–607. Kong F L, Zhang Q Y, Yang F X, Guo H M. Studies on the precision of regional cotton variety trial., 1998, 24: 601–607 (in Chinese with English abstract).

    [12] 許乃銀, 金石橋. 棉花品種區(qū)域試驗(yàn)適宜試驗(yàn)點(diǎn)數(shù)量的抽樣估計(jì). 棉花學(xué)報(bào), 2013, 25: 57–62. Xu N Y, Jin S Q. Sampling estimation of suitable quantity of test sites in cotton variety regional trials., 2013, 25: 57–62 (in Chinese with English abstract).

    [13] Gouy M, Rousselle Y, Bastianelli D, Lecomte P, Bonnal L, Roques D, Efile J C, Rocher S, Daugrois J, Toubi L, Nabeneza S, Hervouet C, Telismart H, Denis M, Thong-Chane A, Glaszmann J C, Hoarau J Y, Nibouche S, Costet L. Experimental assessment of the accuracy of genomic selection in sugarcane., 2013, 126: 2575–2586.

    [14] Iwata H, Jannink J. Accuracy of genomic selection prediction in barley breeding programs: a simulation study based on the real single nucleotide polymorphism data of barley breeding lines., 2011, 51: 1887–1902.

    [15] Kelly A M, Smith A B, Eccleston J A, Cullis B R. The accuracy of varietal selection using factor analytic models for multi-environment plant breeding trials., 2007, 47: 1063–1070.

    [16] Gauch H G, Zobel R W. Accuracy and selection success in yield trial analyses., 1989, 77: 473–481.

    [17] 王潔, 廖琴, 胡小軍, 萬(wàn)建民. 北方稻區(qū)國(guó)家水稻品種區(qū)域試驗(yàn)精確度分析. 作物學(xué)報(bào), 2010, 36: 1870–1876. Wang J, Liao Q, Hu X J, Wan J M. Precision evaluation of rice variety regional trials in Northern China., 2010, 36: 1870–1876 (in Chinese with English abstract).

    [18] Brancourt-Hulmel M, Lecomte C. Effect of environmental variates on genotype × environment interaction of winter wheat: a comparison of biadditive factorial regression to AMMI., 2003, 43: 608–617.

    [19] Brennan P S, Byth D E, Drake D W, DeLacy I H, Butler D G. Determination of the location and number of test environments for a wheat cultivar evaluation program., 1981, 32: 189–201.

    [20] Imrie B C, Drake D W, Delacy I H, Byth D E. Analysis of genotypic and environmental variation in international mungbean trials., 1981, 30: 301–311.

    Optimization of test location number and replicate frequency in regional winter wheat variety trials in northern winter wheat region in China

    ZHANG Yi1, XU Nai-Yin2, GUO Li-Lei1, YANG Zi-Guang3, ZHANG Xiao-Qing1, and YANG Xiao-Ni2

    1National Extension and Service Center of Agricultural Technology, Beijing 100125, China;2Jiangsu Academy of Agricultural Sciences, Nanjing 210014, Jiangsu, China;3Luoyang Academy of Agriculture and Forestry Sciences, Luoyang 471022, Henan, China

    The rational allocation of the replicate frequency and test location number in regional crop trials are highly beneficial to the enhancement of both trial cost efficiency and new cultivar selection efficiency. The rationality of test location number and replicate frequency for national wheat regional trials in the Northern winter wheat region (NWWR) in China was evaluated by using historical experimental datasets during the last 10 years according to the dynamics of trial heritability and noise-signal quotient () with the increase of test locations and replicates within trials in 2010–2019, proposing an optimal design scheme of the replicate frequency and test location number for the wheat planting region. The result indicated three replicates at current achieve 0.87 of the averaged within-trial heritability of single-trials in NWWR, and only 1.4 times of replicates on average was needed to achieve 0.75 of within-trial heritability, so three replicates were obviously sufficient to maintain enough test accuracy. The needed test location number to achieve 0.75 of cross-trial heritability was estimated as 11 for the irrigated group and 13 for the rainfed group in the one-year multi-locational trials in NWWR, while in the current wheat trials that was about 11 and 8 effective test locations, achieving a heritability level of 0.75 and 0.60, respectively. The irrigated group was exactly meet the trial accuracy requirement (= 0.75), while the rainfed group was somewhat insufficient. Considering the importance of the regional wheat trials in recommending new varieties for registration and production, and the possible trial cancellation due to various abnormalities, the practical management should set a heritability level around 0.75, and maintain three replicates, with around 11 test locations for the irrigated group, and 13 locations for the rain-fed group. In order to increase heritability to the level of 0.80, about 16 test locations should be required.

    wheat (L.); regional trial; heritability; noise-signal quotient; number of replicate; number of test location

    10.3724/SP.J.1006.2020.91069

    本研究由國(guó)家農(nóng)作物品種區(qū)域試驗(yàn)專項(xiàng)資助。

    This study was supported by the Special Project of National Crop Regional Trials.

    E-mail: zhangy@oilcrops.cn, Tel: 010-64194512

    2019-11-29;

    2020-03-24;

    2020-05-08.

    URL: http://kns.cnki.net/kcms/detail/11.1809.S.20200507.1620.010.html

    猜你喜歡
    冬麥區(qū)遺傳力單點(diǎn)
    某國(guó)家生豬核心育種場(chǎng)加系大白母豬繁殖性狀遺傳參數(shù)估計(jì)
    美系大白種豬生長(zhǎng)性狀和繁殖性狀遺傳力估計(jì)
    養(yǎng)豬(2022年4期)2022-08-17 07:07:02
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測(cè)速算法
    黃淮和長(zhǎng)江中下游冬麥區(qū)小麥品種(系)農(nóng)藝性狀 及其聚類(lèi)分析*
    超薄異型坯連鑄機(jī)非平衡單點(diǎn)澆鑄實(shí)踐與分析
    山東冶金(2019年5期)2019-11-16 09:09:10
    國(guó)內(nèi)外301份小麥品種(系)種子萌發(fā)期抗旱性鑒定及評(píng)價(jià)
    小麥栽培的特點(diǎn)及不同冬麥區(qū)存在的問(wèn)題
    數(shù)字電視地面?zhèn)鬏斢脝晤l網(wǎng)與單點(diǎn)發(fā)射的效果比較
    16噸單點(diǎn)懸掛平衡軸的優(yōu)化設(shè)計(jì)
    通過(guò)胎盤(pán)效率改良母豬繁殖性能的研究進(jìn)展
    精品人妻在线不人妻| 国产无遮挡羞羞视频在线观看| 亚洲伊人色综图| 精品少妇黑人巨大在线播放| 男人操女人黄网站| svipshipincom国产片| 亚洲婷婷狠狠爱综合网| 成年av动漫网址| 欧美人与善性xxx| 黑人猛操日本美女一级片| 日韩欧美精品免费久久| 男女午夜视频在线观看| 90打野战视频偷拍视频| 亚洲天堂av无毛| 中国国产av一级| 下体分泌物呈黄色| 日韩欧美一区视频在线观看| 制服丝袜香蕉在线| 国产精品免费大片| 国产精品久久久久久精品电影小说| 嫩草影视91久久| 免费少妇av软件| 色网站视频免费| 欧美日韩福利视频一区二区| 91精品伊人久久大香线蕉| 亚洲精品成人av观看孕妇| 一区福利在线观看| 丝袜人妻中文字幕| 观看美女的网站| 麻豆精品久久久久久蜜桃| 一区二区日韩欧美中文字幕| 婷婷色麻豆天堂久久| 国产亚洲欧美精品永久| 老司机深夜福利视频在线观看 | 少妇的丰满在线观看| 少妇猛男粗大的猛烈进出视频| 美女视频免费永久观看网站| videosex国产| 女人精品久久久久毛片| 老司机亚洲免费影院| 涩涩av久久男人的天堂| 男女下面插进去视频免费观看| av又黄又爽大尺度在线免费看| 无限看片的www在线观看| 亚洲少妇的诱惑av| 丝袜美腿诱惑在线| 亚洲一级一片aⅴ在线观看| 丰满少妇做爰视频| 男人添女人高潮全过程视频| 亚洲精品自拍成人| 免费av中文字幕在线| a 毛片基地| 久久婷婷青草| 一区二区三区精品91| 成人漫画全彩无遮挡| 亚洲欧美激情在线| 国产在线视频一区二区| 性少妇av在线| 欧美日韩av久久| 国产精品香港三级国产av潘金莲 | 精品久久蜜臀av无| 久热爱精品视频在线9| 日韩制服丝袜自拍偷拍| 精品人妻在线不人妻| 一级爰片在线观看| 国产精品久久久av美女十八| 成年动漫av网址| 捣出白浆h1v1| 制服丝袜香蕉在线| 国产欧美日韩综合在线一区二区| 可以免费在线观看a视频的电影网站 | e午夜精品久久久久久久| 国产精品欧美亚洲77777| 校园人妻丝袜中文字幕| 国产精品嫩草影院av在线观看| 青春草国产在线视频| 国产爽快片一区二区三区| avwww免费| 久久久久精品人妻al黑| 两个人免费观看高清视频| 成人毛片60女人毛片免费| 最新的欧美精品一区二区| 国产精品欧美亚洲77777| 亚洲精品第二区| 国产成人午夜福利电影在线观看| 大香蕉久久成人网| 99久久综合免费| 极品少妇高潮喷水抽搐| 老鸭窝网址在线观看| 亚洲欧洲日产国产| 天堂俺去俺来也www色官网| 免费在线观看完整版高清| 亚洲精品自拍成人| 国产1区2区3区精品| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 欧美少妇被猛烈插入视频| 精品少妇一区二区三区视频日本电影 | 亚洲精品aⅴ在线观看| 亚洲精品国产一区二区精华液| 777久久人妻少妇嫩草av网站| 80岁老熟妇乱子伦牲交| 久久99精品国语久久久| 久久韩国三级中文字幕| 日韩视频在线欧美| 欧美日韩成人在线一区二区| av在线播放精品| 免费黄网站久久成人精品| 夫妻性生交免费视频一级片| 中文字幕色久视频| 免费黄色在线免费观看| 五月开心婷婷网| 桃花免费在线播放| 久久久久国产精品人妻一区二区| 18禁国产床啪视频网站| 极品人妻少妇av视频| 国产在视频线精品| 亚洲国产精品999| 免费观看a级毛片全部| 久久 成人 亚洲| 无遮挡黄片免费观看| 欧美在线一区亚洲| 国产99久久九九免费精品| 日韩制服丝袜自拍偷拍| 日韩人妻精品一区2区三区| 久久影院123| 欧美亚洲日本最大视频资源| 成年美女黄网站色视频大全免费| 久久久久久久久免费视频了| 在线观看人妻少妇| 汤姆久久久久久久影院中文字幕| 亚洲色图综合在线观看| 高清在线视频一区二区三区| 交换朋友夫妻互换小说| 黑人欧美特级aaaaaa片| 18禁动态无遮挡网站| 多毛熟女@视频| 欧美人与性动交α欧美精品济南到| 免费av中文字幕在线| av在线播放精品| 久久综合国产亚洲精品| 亚洲欧美精品综合一区二区三区| 成年动漫av网址| 欧美黑人欧美精品刺激| 电影成人av| 久久久久视频综合| 黄色毛片三级朝国网站| 老司机深夜福利视频在线观看 | 精品久久蜜臀av无| 免费日韩欧美在线观看| 亚洲国产毛片av蜜桃av| 亚洲七黄色美女视频| 777米奇影视久久| 国产1区2区3区精品| 国产爽快片一区二区三区| 亚洲久久久国产精品| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜爱| 少妇 在线观看| 国产日韩欧美在线精品| 久久人妻熟女aⅴ| 少妇人妻久久综合中文| 丁香六月欧美| 亚洲精品国产区一区二| 80岁老熟妇乱子伦牲交| 亚洲人成电影观看| 成人漫画全彩无遮挡| 免费黄色在线免费观看| 亚洲精品一二三| 亚洲精品视频女| 激情五月婷婷亚洲| 久久久久精品性色| 欧美97在线视频| 国产黄色免费在线视频| 亚洲精品国产色婷婷电影| 精品少妇久久久久久888优播| 一二三四在线观看免费中文在| a级片在线免费高清观看视频| 只有这里有精品99| 人人妻人人添人人爽欧美一区卜| 在线 av 中文字幕| 老熟女久久久| 国产伦理片在线播放av一区| 久久久久久久久久久久大奶| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频精品| 中文字幕制服av| 精品人妻熟女毛片av久久网站| 18在线观看网站| 男女国产视频网站| 女性被躁到高潮视频| 亚洲欧美色中文字幕在线| 精品国产露脸久久av麻豆| 精品国产一区二区久久| 久久青草综合色| 免费看av在线观看网站| 久久av网站| 免费高清在线观看视频在线观看| 亚洲精品国产色婷婷电影| 曰老女人黄片| 美女国产高潮福利片在线看| 久久青草综合色| 国产精品99久久99久久久不卡 | 天天添夜夜摸| 最黄视频免费看| 看免费av毛片| 18禁动态无遮挡网站| 欧美日本中文国产一区发布| 亚洲人成电影观看| 各种免费的搞黄视频| 亚洲国产中文字幕在线视频| 亚洲自偷自拍图片 自拍| 老司机在亚洲福利影院| 久久狼人影院| 男女下面插进去视频免费观看| 2018国产大陆天天弄谢| 亚洲一区中文字幕在线| 欧美在线黄色| 国产成人91sexporn| a 毛片基地| 成人亚洲精品一区在线观看| 久久综合国产亚洲精品| 国产又爽黄色视频| 操美女的视频在线观看| 久久久久久人人人人人| 汤姆久久久久久久影院中文字幕| 国产老妇伦熟女老妇高清| 欧美日本中文国产一区发布| 国产精品一区二区在线不卡| 国产淫语在线视频| 99热网站在线观看| 国产欧美日韩一区二区三区在线| 欧美亚洲 丝袜 人妻 在线| 丰满少妇做爰视频| 性高湖久久久久久久久免费观看| 我的亚洲天堂| 亚洲人成网站在线观看播放| 久久久久视频综合| 精品亚洲乱码少妇综合久久| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 久久国产精品大桥未久av| 欧美另类一区| 午夜影院在线不卡| 少妇人妻久久综合中文| 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 好男人视频免费观看在线| 久热这里只有精品99| 99热网站在线观看| 熟妇人妻不卡中文字幕| 男人操女人黄网站| 黄色视频在线播放观看不卡| 满18在线观看网站| 亚洲自偷自拍图片 自拍| 黄色视频不卡| 日本av手机在线免费观看| 秋霞伦理黄片| 99热国产这里只有精品6| 侵犯人妻中文字幕一二三四区| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 精品人妻在线不人妻| 桃花免费在线播放| 亚洲伊人色综图| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 亚洲第一区二区三区不卡| 啦啦啦中文免费视频观看日本| 欧美变态另类bdsm刘玥| 一边亲一边摸免费视频| 国产精品99久久99久久久不卡 | 国产在视频线精品| 女性生殖器流出的白浆| 亚洲精品av麻豆狂野| 热99久久久久精品小说推荐| 国产精品亚洲av一区麻豆 | 少妇的丰满在线观看| 999久久久国产精品视频| 视频在线观看一区二区三区| 黄色视频在线播放观看不卡| 精品福利永久在线观看| 亚洲四区av| 亚洲av电影在线进入| 精品一区二区三区av网在线观看 | 国产国语露脸激情在线看| 亚洲五月色婷婷综合| 中文字幕另类日韩欧美亚洲嫩草| 咕卡用的链子| 亚洲av电影在线进入| 日韩成人av中文字幕在线观看| 一区二区日韩欧美中文字幕| 国产成人啪精品午夜网站| 极品少妇高潮喷水抽搐| 99久久精品国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 一级毛片我不卡| 五月开心婷婷网| 久久久久久久久免费视频了| 国产毛片在线视频| 制服人妻中文乱码| 国产精品久久久人人做人人爽| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 黄片无遮挡物在线观看| 国产乱来视频区| 成人毛片60女人毛片免费| 美国免费a级毛片| 久久久国产精品麻豆| 久久久久人妻精品一区果冻| 午夜福利乱码中文字幕| www.精华液| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 老熟女久久久| 精品一区二区免费观看| 韩国精品一区二区三区| 亚洲国产av新网站| 精品亚洲乱码少妇综合久久| 日日撸夜夜添| av在线app专区| 国产国语露脸激情在线看| 国产精品久久久久久久久免| 别揉我奶头~嗯~啊~动态视频 | 又黄又粗又硬又大视频| 国产男女超爽视频在线观看| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 中文字幕av电影在线播放| 亚洲成av片中文字幕在线观看| 国产亚洲最大av| 久热这里只有精品99| 高清av免费在线| 国产xxxxx性猛交| 99热网站在线观看| 免费黄频网站在线观看国产| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 日韩不卡一区二区三区视频在线| 99久久人妻综合| 视频区图区小说| 女人被躁到高潮嗷嗷叫费观| 日本一区二区免费在线视频| 999久久久国产精品视频| 美女主播在线视频| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 国产精品久久久久久久久免| 日韩伦理黄色片| 亚洲国产成人一精品久久久| 国产极品粉嫩免费观看在线| 晚上一个人看的免费电影| 亚洲色图 男人天堂 中文字幕| 在线观看www视频免费| 大香蕉久久成人网| 婷婷成人精品国产| 国产欧美日韩综合在线一区二区| 久久精品国产a三级三级三级| 少妇人妻 视频| 亚洲美女搞黄在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久人人人人人| 国产麻豆69| 乱人伦中国视频| 国产精品秋霞免费鲁丝片| 老熟女久久久| 各种免费的搞黄视频| 少妇精品久久久久久久| 十分钟在线观看高清视频www| 日日啪夜夜爽| 国产成人免费观看mmmm| 久久国产精品男人的天堂亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲第一av免费看| 高清av免费在线| 欧美97在线视频| 国产男女内射视频| 在线观看免费视频网站a站| 亚洲国产中文字幕在线视频| 国产不卡av网站在线观看| av天堂久久9| 亚洲成人国产一区在线观看 | 热re99久久精品国产66热6| 青草久久国产| 中文字幕亚洲精品专区| 亚洲国产av新网站| 在线观看一区二区三区激情| 国产精品免费大片| a级片在线免费高清观看视频| 黄片播放在线免费| 国产精品.久久久| 亚洲国产欧美网| 亚洲男人天堂网一区| 久久久精品94久久精品| 1024香蕉在线观看| av国产精品久久久久影院| 欧美日本中文国产一区发布| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 色婷婷久久久亚洲欧美| 一级片免费观看大全| 国产精品一区二区精品视频观看| 亚洲国产中文字幕在线视频| 精品一品国产午夜福利视频| 精品第一国产精品| 99国产精品免费福利视频| 久久综合国产亚洲精品| 成人免费观看视频高清| 中文字幕精品免费在线观看视频| 纵有疾风起免费观看全集完整版| 日韩一区二区三区影片| 亚洲熟女精品中文字幕| 99久久人妻综合| 成人国语在线视频| 一区福利在线观看| 久久人人爽av亚洲精品天堂| 国产成人精品福利久久| 亚洲国产最新在线播放| 亚洲精品aⅴ在线观看| 在线免费观看不下载黄p国产| 久久精品亚洲av国产电影网| 女人被躁到高潮嗷嗷叫费观| 日日撸夜夜添| 亚洲精品久久成人aⅴ小说| www.熟女人妻精品国产| 99国产精品免费福利视频| 卡戴珊不雅视频在线播放| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 国产精品 国内视频| 亚洲欧美日韩另类电影网站| 别揉我奶头~嗯~啊~动态视频 | 飞空精品影院首页| 国产成人免费观看mmmm| 秋霞伦理黄片| 久久久久久久精品精品| 国产av精品麻豆| 超碰成人久久| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| av免费观看日本| 自线自在国产av| 桃花免费在线播放| 久久精品熟女亚洲av麻豆精品| 久久久久人妻精品一区果冻| 成人免费观看视频高清| 老司机影院成人| 国产黄色免费在线视频| 三上悠亚av全集在线观看| 国产片特级美女逼逼视频| 欧美精品人与动牲交sv欧美| 啦啦啦啦在线视频资源| 久久久久国产一级毛片高清牌| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 国产色婷婷99| 久久综合国产亚洲精品| 国产男女超爽视频在线观看| 丝袜脚勾引网站| 街头女战士在线观看网站| 韩国精品一区二区三区| 日韩av免费高清视频| 1024香蕉在线观看| 91精品伊人久久大香线蕉| 午夜福利在线免费观看网站| 久久久久精品性色| 大香蕉久久成人网| 伊人久久国产一区二区| 国产精品免费大片| 精品人妻熟女毛片av久久网站| 久久精品国产综合久久久| 欧美精品一区二区大全| 妹子高潮喷水视频| 亚洲国产av新网站| 精品人妻在线不人妻| 校园人妻丝袜中文字幕| 色网站视频免费| 欧美精品亚洲一区二区| 日本色播在线视频| tube8黄色片| 国产男人的电影天堂91| 免费看不卡的av| 一二三四中文在线观看免费高清| 国产免费福利视频在线观看| 精品国产超薄肉色丝袜足j| 国产一区二区在线观看av| 日韩视频在线欧美| 校园人妻丝袜中文字幕| 免费观看av网站的网址| 性色av一级| 亚洲国产av影院在线观看| 青春草视频在线免费观看| 久久热在线av| 天天添夜夜摸| 亚洲欧美日韩另类电影网站| 一区二区av电影网| 两个人免费观看高清视频| 综合色丁香网| 国产伦人伦偷精品视频| 性色av一级| 国产成人免费观看mmmm| 国产一卡二卡三卡精品 | 黄片小视频在线播放| 精品卡一卡二卡四卡免费| 老鸭窝网址在线观看| 欧美日韩亚洲高清精品| 又粗又硬又长又爽又黄的视频| 一本一本久久a久久精品综合妖精| 欧美日韩国产mv在线观看视频| 亚洲av欧美aⅴ国产| 美女视频免费永久观看网站| 久久久久国产一级毛片高清牌| 午夜福利,免费看| 精品久久久精品久久久| 18在线观看网站| 乱人伦中国视频| 欧美少妇被猛烈插入视频| 国产精品.久久久| 欧美日韩精品网址| 日韩精品有码人妻一区| 午夜福利乱码中文字幕| 美女福利国产在线| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线| 国产男女内射视频| av有码第一页| 777久久人妻少妇嫩草av网站| 别揉我奶头~嗯~啊~动态视频 | 精品亚洲成国产av| 国产人伦9x9x在线观看| 精品国产一区二区三区久久久樱花| 一二三四在线观看免费中文在| 少妇人妻 视频| 欧美变态另类bdsm刘玥| av国产久精品久网站免费入址| 男女边摸边吃奶| 久久久久久人妻| 大码成人一级视频| 不卡视频在线观看欧美| 国产老妇伦熟女老妇高清| av福利片在线| 人妻 亚洲 视频| 嫩草影院入口| 汤姆久久久久久久影院中文字幕| 涩涩av久久男人的天堂| 天堂8中文在线网| 国产成人精品久久久久久| 蜜桃国产av成人99| 一级片'在线观看视频| 黄色一级大片看看| 曰老女人黄片| 国产麻豆69| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 如日韩欧美国产精品一区二区三区| 久久久国产欧美日韩av| 精品第一国产精品| 制服丝袜香蕉在线| 18禁动态无遮挡网站| 男女边吃奶边做爰视频| 人人妻人人澡人人爽人人夜夜| 国产亚洲av高清不卡| 日本vs欧美在线观看视频| 免费在线观看完整版高清| 无限看片的www在线观看| 波多野结衣一区麻豆| 精品一区二区三卡| 美女福利国产在线| 日韩大码丰满熟妇| 在线观看人妻少妇| 亚洲 欧美一区二区三区| 久久久久久久久久久久大奶| 伊人久久国产一区二区| 久久久久久久久免费视频了| 亚洲av日韩在线播放| 男人舔女人的私密视频| av一本久久久久| 你懂的网址亚洲精品在线观看| 黄色 视频免费看| 亚洲精品国产一区二区精华液| 亚洲成色77777| 丰满饥渴人妻一区二区三| 国产男人的电影天堂91| 久久久久久久久久久免费av| 免费黄色在线免费观看| 国产成人欧美在线观看 | 日本wwww免费看| 久久av网站| 成年人午夜在线观看视频| 精品国产一区二区久久| 久久av网站| 一边摸一边抽搐一进一出视频| 亚洲欧美成人综合另类久久久| 国产不卡av网站在线观看| 99久国产av精品国产电影| 精品国产乱码久久久久久小说| 国产成人午夜福利电影在线观看| 欧美人与性动交α欧美软件| 最黄视频免费看| 亚洲精华国产精华液的使用体验| 最近手机中文字幕大全| 最黄视频免费看| 亚洲欧美成人精品一区二区| 老司机影院成人| 少妇人妻久久综合中文| 捣出白浆h1v1| 亚洲五月色婷婷综合| 最新在线观看一区二区三区 | 男女之事视频高清在线观看 | 午夜福利一区二区在线看| 夜夜骑夜夜射夜夜干| 日本一区二区免费在线视频| 久久av网站| 国产精品99久久99久久久不卡 |