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

    高低氮處理下小麥旗葉性狀的遺傳分析

    2022-02-22 07:31:56武亞瑞劉錫建楊國敏劉紅偉孔文超吳永振孫晗秦冉崔法趙春華
    中國農業(yè)科學 2022年1期
    關鍵詞:加性旗葉葉面積

    武亞瑞,劉錫建,楊國敏,劉紅偉,孔文超,吳永振,孫晗,秦冉,崔法,趙春華

    高低氮處理下小麥旗葉性狀的遺傳分析

    武亞瑞1,劉錫建1,楊國敏2,劉紅偉1,孔文超1,吳永振1,孫晗1,秦冉1,崔法1,趙春華1

    1魯東大學農學院,山東煙臺 264025;2菏澤市土壤肥料工作站,山東菏澤 274000

    【】旗葉是小麥光合碳固定的重要場所,對小麥產量起十分重要的作用。研究小麥旗葉在高、低氮環(huán)境下的遺傳特性,分析其遺傳機制,為優(yōu)異株型育種、高產育種提供參考依據。以科農9204和京411為親本所構建的188個RIL群體為材料,分別種植在6個不同的高、低氮環(huán)境下,通過對群體旗葉性狀調查并進行遺傳分析,從而確定控制各性狀的基因數目,估計遺傳效應值及遺傳率,并對小麥旗葉性狀與產量之間的關系進行分析。在遺傳估測中,低氮環(huán)境下:旗葉長在E3環(huán)境的最適遺傳模型為2MG-CE,即2對互補作用主基因遺傳模型,其加性×加性上位性互作效應值為1.098,主基因遺傳率為31.35%,在其他低氮環(huán)境下均表現為多基因遺傳;旗葉寬均表現為多基因遺傳;旗葉面積(除E5)的最適遺傳模型均為2MG-CE,加性×加性上位性互作效應值為1.884,主基因遺傳率為36.7%,在E5為多基因遺傳。高氮環(huán)境下:旗葉長(除E4)的最適遺傳模型為2MG-CE,加性×加性上位性互作效應值為1.133,主基因遺傳率為32.6%,在E4環(huán)境的最適遺傳模型為2MG-ER,即2對隱性上位主基因遺傳模型,其第一對主基因的加性效應值為1.431,第二對主基因的加性效應值為1.108,主基因遺傳率為51.77%;旗葉寬(除E2)的最適遺傳模型為2MG-CE,加性×加性上位性互作效應值為0.119,主基因遺傳率為37.29%,在E2表現為多基因遺傳;旗葉面積的最適遺傳模型為2MG-CE,加性×加性上位性互作效應值為3.067,主基因遺傳率為44.42%。旗葉性狀在不同環(huán)境的遺傳模型不同,在高氮環(huán)境下遺傳較為穩(wěn)定,在低氮下受環(huán)境影響較大。在旗葉與產量性狀的相關性分析中,旗葉性狀與穗粒數、穗粒重、單株產量之間呈顯著正相關,且在不同環(huán)境下的影響程度不同。旗葉性狀易受外界環(huán)境影響,在高、低氮環(huán)境下的表現不同。旗葉在低氮環(huán)境下表現為不同的主基因遺傳和多基因遺傳;在高氮環(huán)境下大多表現為主基因遺傳,由2對基因控制,并且存在基因之間的相互作用,且可能存在效應較大的主效QTL。改善旗葉性狀可以提高小麥的單株產量、穗粒重等產量性狀。

    小麥;旗葉性狀;遺傳模型;產量;相關分析

    0 引言

    【研究意義】植物葉片的光合作用是轉換能量生產營養(yǎng)物質的基礎[1],可作為驗證植物生長狀況的指標[2-4]。旗葉是小麥進行光合作用的主要場所,對產量起非常重要的作用[5-7]。氮素作為作物生長需要的大量營養(yǎng)元素之一,施用氮肥較少或過量均可能對作物產量和品質產生一定的影響[8-10]。因此,利用主基因+多基因數量性狀混合遺傳模型[11-14],分析小麥旗葉在不同氮素環(huán)境下的遺傳特點,可為合理施用氮肥,提高小麥產量提供參考。【前人研究進展】植物的大多性狀是數量性狀。隨著科技的發(fā)展,將數量性狀與植物遺傳育種結合,利用數學統(tǒng)計和遺傳模型分析群體數量性狀遺傳規(guī)律,可為現代植物育種提供科學依據[15-18]。目前,對小麥旗葉性狀QTL定位有較多的研究[19-23],但對其遺傳模型的研究較少。主基因+多基因數量性狀混合遺傳模型[11-14]已經在大量作物中廣泛應用[24-28]。如江建華等[29]對水稻柱頭性狀進行研究,測得該性狀由2對主效基因和微效基因共同控制;黃冰艷等[30]利用F2群體分析花生油酸和亞油酸含量的遺傳方式,結果表明,這兩個性狀均為2對主基因加性-顯性-上位性遺傳模型,并且存在多基因效應和基因互作;解松峰等[31-32]利用小麥重組自交系(RIL群體)對產量相關性狀進行研究,結果表明,穗部性狀以多基因遺傳為主,其他產量性狀(單株產量、千粒重、株高、穗下節(jié)間長、旗葉上節(jié)間長和分蘗數)在不同環(huán)境下受到主基因+多基因混合遺傳控制?!颈狙芯壳腥朦c】迄今為止,對于高、低氮環(huán)境下小麥旗葉性狀的數量性狀遺傳分析仍鮮見報道,中國科學院遺傳與發(fā)育生物學研究所農業(yè)資源研究中心李俊明研究員課題組以科農9204和京411為親本構建了重組自交系群體?!緮M解決的關鍵問題】本研究通過對構建的RIL群體(F7:8和F8:9)旗葉性狀的表型數據進行分析,運用主基因+多基因混合遺傳模型的方法確定旗葉性狀在高、低氮環(huán)境下的遺傳模型和基因作用方式,為小麥旗葉性狀的遺傳規(guī)律、育種研究提供理論參考。

    1 材料與方法

    1.1 試驗材料

    供試材料為科農9204和京411作為親本組成的重組自交系群體(KJ-RIL),由中國科學院遺傳與發(fā)育生物學研究所農業(yè)資源研究中心李俊明研究員課題組提供。田間種植包括6個環(huán)境,將親本和RIL群體分別在2011—2012年石家莊低氮環(huán)境(E1)、2011—2012年石家莊高氮環(huán)境(E2)、2012—2013年石家莊低氮環(huán)境(E3)、2012—2013年石家莊高氮環(huán)境(E4)、2012—2013年北京低氮環(huán)境(E5)、2012—2013年北京高氮環(huán)境(E6)進行種植。每個環(huán)境中0—20 cm土層內的土壤硝態(tài)氮及總氮含量見文獻[33]。每個高氮小區(qū)播種前施磷酸二胺300 kg·hm-2、尿素150 kg·hm-2,每年伸長期施尿素150 kg·hm-2。在低氮小區(qū),小麥生長期不施氮肥(缺氮)。每個材料種植2行,每行40粒種子,2 m行長,行距0.25 m。各環(huán)境隨機區(qū)組設計,設置2個重復。

    1.2 試驗方法

    每個環(huán)境在小麥開花后15 d進行旗葉長和旗葉寬的測定,每個家系選擇5株長勢一致的單株進行主莖穗旗葉長和旗葉寬的測量,分別量取旗葉最長和最寬部位。旗葉面積為旗葉長×旗葉寬×0.83。其他農藝性狀的測量參考Fan等[34]方法。

    1.3 數據分析

    應用EXCEL2010和SPSS 24.0軟件對旗葉性狀表型數據分析極值、變異系數、偏度和峰度等,并分析不同環(huán)境下旗葉與產量之間的相關性。根據植物數量性狀主基因+多基因混合遺傳模型分析方法,利用章元明教授課題組開發(fā)的Windows軟件包中的SEA對小麥旗葉性狀進行遺傳分析[11-14]。通過極大似然法和ECM算法等對混合分布中的有關成分分布參數做出估計,然后通過AIC值的判別選出AIC值最小的模型或AIC值相近的幾個模型再根據均勻性檢驗(12、22、32)、Smirnov檢驗(2)和Kolmogorov檢驗()進行適合性檢驗選擇出最適遺傳模型。最后根據最小二乘法估計各性狀的主基因和多基因效應值等遺傳參數。

    2 結果

    2.1 旗葉性狀的表型數據分析

    對不同環(huán)境下2個親本進行方差分析(表1)。結果表明,環(huán)境對旗葉長、旗葉寬和旗葉面積的影響都極顯著;親本的遺傳對旗葉寬和旗葉面積的影響極顯著,對旗葉長的影響不顯著。對小麥RIL群體旗葉性狀的表型進行測定,并分析了該群體旗葉長、寬、面積在6個環(huán)境下的表現(表2)。低氮環(huán)境下,親本科農9204和京411的旗葉長分別為13.08和13.81 cm,旗葉寬分別為1.49和1.34 cm,旗葉面積分別為16.10和15.67 cm2。群體的旗葉長介于9.33—19.45 cm,均值為14.47 cm,超過兩親本均值;群體的旗葉寬介于1—1.73 cm,均值為1.33 cm,偏向于低值親本;群體的旗葉面積介于9.26—25.88 cm2,均值為16.04 cm2,略低于高值親本。高氮環(huán)境下,親本科農9204和京411的旗葉長分別為18.26和17.87 cm,旗葉寬分別為1.64和1.42 cm,旗葉面積分別為25.04和21.85 cm2。群體的旗葉長介于12.96—24.62 cm,均值為18.24 cm,略低于高值親本;群體的旗葉寬介于1.07—2.05 cm,均值為1.49 cm,偏向于低值親本;群體的旗葉面積介于13.88—36.44 cm2,均值為22.91 cm2,略高于低值親本。旗葉性狀在高氮環(huán)境下表型值整體高于低氮環(huán)境。

    表1 不同環(huán)境下兩親本旗葉性狀的方差分析

    FLL:旗葉長;FLW:旗葉寬;FLA:旗葉面積。下同

    FLL: Flag leaf length; FLW: Flag leaf width; FLA: Flag leaf area. The same as below

    表2 不同生境下旗葉性狀的統(tǒng)計值

    LN:低氮;HN:高氮。下同 LN: Low nitrogen; HN: High nitrogen. The same as below

    6個環(huán)境下,各性狀極值變化范圍很大,變異幅度大,適合進行遺傳分析。旗葉面積的變異系數(13%—19%)略高于旗葉長(9%—12%)和旗葉寬(9%—11%),說明旗葉面積在形態(tài)變異上比旗葉長和旗葉寬豐富。整體旗葉性狀的變異系數為9%—19%,表明群體表型變異穩(wěn)定。整體旗葉長、寬和面積的偏度和峰度系數絕對值均小于1,表現為數量性狀遺傳特征。

    由圖1可知,旗葉性狀在不同環(huán)境下均表現為正態(tài)分布,符合遺傳分析的要求。旗葉性狀在高、低氮環(huán)境下均表現為單峰分布,符合主基因+多基因表型分布特征。

    LN:低氮;HN:高氮;FLL:旗葉長;FLW:旗葉寬;FLA:旗葉面積。E1:2011—2012年石家莊低氮;E2:2011—2012年石家莊高氮;E3:2012—2013年石家莊低氮;E4:2012—2013年石家莊高氮;E5:2012—2013年北京低氮;E6:2012—2013年北京高氮

    2.2 旗葉性狀高低氮遺傳模型分析

    利用植物數量性狀主基因+多基因混合遺傳模型分析小麥旗葉性狀,選取AIC值最小或與AIC最小值比較接近的幾個模型作為備選模型,對備選模型進行適合性檢驗從而得到旗葉性狀的最適遺傳模型(表3)。結果表明,低氮環(huán)境下,旗葉長在E3的最適遺傳模型為2MG-CE,該模型具有最小的AIC值,為1 464.11;旗葉面積在E1、E3的最適遺傳模型2MG-CE,該模型具有最小的AIC值,分別為1 714.21和1 890.73;旗葉長在E1、E5與旗葉寬在低氮環(huán)境下以及旗葉面積在E5的最適遺傳模型均為0MG,具有最小的AIC值。高氮環(huán)境下,旗葉長在E2、E6的最適遺傳模型為2MG- CE,具有最小的AIC值,分別為1 470.30和1 493.63,在E4的最適遺傳模型2MG-ER,該模型具有最小的AIC值,為1 684.05;旗葉寬在E2的最適遺傳模型為0MG,該模型有最小的AIC值,為-82.91,旗葉寬在E4、E6環(huán)境下的最適遺傳模型均為2MG-CE,具有最小的AIC值;旗葉面積的最適遺傳模型為2MG-CE。

    表3 旗葉性狀在重組自交系群體中的最適遺傳模型分析

    B模型表示2對主基因,無多基因遺傳模型;—:空缺;MG:主基因模型;ER:隱性上位;CE:互補作用。下同

    B: Two major genes without polygene; —: the value is absent; MG: Major gene model; ER: Epistasis recessively; CE: Complementary effect. The same as below

    2.3 旗葉性狀高低氮遺傳參數估計

    通過R語言的SEA軟件包對旗葉性狀最適遺傳模型估計其一階及二階遺傳參數(表4)(部分一階遺傳參數、多基因遺傳方差及遺傳率為空值)。低氮環(huán)境下,旗葉長在E3的主基因遺傳率為31.35%,加性×加性上位性互作效應值為1.098;旗葉面積在E1的主基因遺傳率為38.35%,加性×加性上位性互作效應值為1.726;旗葉面積在E3的主基因遺傳率為34.97%,加性×加性上位性互作效應值為2.042。高氮環(huán)境下,旗葉長在E2、E6的主基因遺傳率分別為42.41%和22.79%,加性×加性上位性互作效應值分別為1.302和0.963,因此,旗葉長表現為2對互補作用主基因遺傳時,主基因遺傳率相差較大,多基因遺傳率未能檢出;旗葉長在E4的主基因遺傳率為51.77%,第一對主基因的加性效應值為1.431,第二對主基因的加性效應值為1.108,旗葉長在該環(huán)境表現為2對隱性上位主基因遺傳,主基因遺傳率較高,遺傳較為穩(wěn)定;旗葉寬在E4、E6的主基因遺傳率分別為38.58%和35.99%,加性×加性上位性互作效應值分別為0.122和0.116;旗葉面積在E2、E4、E6的主基因遺傳率分別為42.77%、44.54%和45.96%,加性×加性上位性互作效應值分別為2.39、3.716和3.095。

    旗葉性狀在其余環(huán)境下都表現為多基因遺傳模型,但是檢測不到基因遺傳的數據特征,沒有辦法顯示多基因效應,其是否存在多基因修飾需要進行多世代的聯(lián)合分離分析,從而確定數量性狀的遺傳特征。

    2.4 小麥RIL群體旗葉性狀與產量性狀的相關性分析

    對旗葉性狀與產量性狀進行相關性分析(表5)。旗葉長在4個環(huán)境(除E5、E6)與穗粒數呈顯著正相關;所有環(huán)境下與穗粒重顯著正相關;4個環(huán)境(除E2、E6)與單株產量顯著正相關;在E4和E5與有效穗數顯著正相關;在E5環(huán)境下與千粒重顯著正相關。旗葉寬在4個環(huán)境下(除E1、E2)與穗粒數正相關性顯著;在5個環(huán)境下(除E2)與穗粒重顯著正相關;在2個環(huán)境下(E3、E4)與單株產量顯著正相關;在E2環(huán)境與有效穗數顯著正相關;在E6環(huán)境下與有效穗數、千粒重顯著負相關。旗葉面積在所有環(huán)境下與穗粒數、穗粒重顯著正相關;在5個環(huán)境下(除E6)與單株產量正相關性顯著;在E4環(huán)境與有效穗數顯著正相關;在E6環(huán)境與有效穗數顯著負相關。

    表4 小麥RIL群體旗葉性狀的各階遺傳參數

    m:群體均方;d:主基因效應值;da:第一對主基因的加性效應值;db:第二對主基因的加性效應值;iab(i*):第一對主基因的加性效應×第二對主基因的加性上位性互作效應;2:群體方差;2:主基因遺傳方差;2:誤差方差(環(huán)境方差);2(%):主基因遺傳率

    m: The mean value of P 1 generation; d: Main gene effect value; da: Addictive effect of the first pair major gene; db: Addictive effect of the second pair major gene; iab(i*): Addictive effect plus addictive effect of the 1st pair major genes × the 2ndpair major gene;2: Phenotypic variance;2: Major gene variance;2: error variance;2(%): Heritability of major gene

    3 討論

    前人對旗葉的研究多側重于QTL定位分析[19-23]、旗葉光合特性及其生理特性[35-37]等,而對其遺傳模型的分析較少。Khan等[38]對2個小麥雜交組合的P1、P2、F1、F2、B1和B26個基本代的旗葉面積進行了混合遺傳分析,結果表明,該性狀主要受1個主基因與多基因組合的控制,主基因遺傳力高于多基因遺傳力,遺傳變異受環(huán)境波動的影響。本研究旗葉面積在大多數環(huán)境下表現為2對互補作用的主基因遺傳,且由環(huán)境引起的變異占總變異的幅度較大。Fan等[34]利用該RIL群體的旗葉性狀進行QTL定位分析,在8個不同環(huán)境中共檢測到38個與旗葉性狀相關的加性QTL,其中,旗葉寬、與旗葉面積的QTL在5個環(huán)境以上均表現穩(wěn)定,為主效穩(wěn)定QTL,旗葉長沒有檢測到主效穩(wěn)定的QTL。

    本研究表明旗葉性狀在高氮環(huán)境下受2對主基因及基因間的作用共同控制,旗葉長的最適遺傳模型存在差異。在低氮環(huán)境下,旗葉性狀有許多為0MG的最適遺傳模型,表現為多基因遺傳,環(huán)境對2對主基因表達的影響較大。旗葉性狀在各環(huán)境下的主基因遺傳率為31.35%—51.77%,且由環(huán)境引起的變異較大,表明在旗葉性狀的遺傳過程中環(huán)境與基因之間存在互作,因此,在實踐過程中應充分考慮基因和環(huán)境效應的相互影響。旗葉長主要表現為主基因遺傳,且在E4環(huán)境下的主基因遺傳率較大,受環(huán)境影響較小,可以作為改善旗葉性狀早期選擇的直接指標,提高育種效率。旗葉寬在低氮環(huán)境下表現為多基因遺傳,受環(huán)境影響較大。旗葉面積在不同環(huán)境下的最佳遺傳模型相同,可能與一因多效或基因連鎖有關。在小麥育種過程中,要加強理論與實踐的結合,并且考慮不同生態(tài)環(huán)境的影響。

    表5 小麥RIL群體的旗葉性狀與產量性狀間的相關性

    *:<0.05水平相關性顯著;**:<0.01水平相關性顯著

    *: significant correlation at<0.05; **: significant correlation at<0.01

    鈕力亞等[39]研究表明旗葉面積與穗粒數、產量之間存在顯著或極顯著的正相關;王敏等[40]發(fā)現旗葉性狀與穗粒重、穗粒數和小區(qū)產量呈現極顯著正相關;黃杰等[41]研究表明旗葉寬、旗葉面積與穗粒數顯著正相關;姜莉莉等[42]研究表明旗葉寬與千粒重呈負相關。本研究結果與其一致。因此,在小麥品種選育中,要充分考慮旗葉對產量的影響,選擇合適的旗葉長和旗葉寬以提高小麥產量。

    4 結論

    明確了小麥旗葉性狀是數量性狀并且可以進行遺傳分析,群體的正態(tài)分布圖可以較好地反映旗葉表型分布特征。旗葉主要由2對主基因控制,并且存在基因間的相互作用,主基因遺傳率為31%—51%。由于環(huán)境和遺傳的共同作用,使小麥旗葉的遺傳模型也有所差異。明確了小麥旗葉與產量之間存在明顯聯(lián)系,與穗粒數、穗粒重顯著相關。

    [1] 夏江寶, 張光燦, 許景偉, 李傳榮. 干旱脅迫下常春藤凈光合速率日變化及其影響因子分析. 水土保持通報, 2010, 30(3): 78-82.

    XIA J B, ZHANG G C, XU J W, LI C R. Analysis on diurnal variation of net photosynthetic rate and its influencing factors under drought stress. Bulletin of Soil and Water Conservation, 2010, 30(3): 78-82. (in Chinese)

    [2] 馬富舉, 李丹丹, 蔡劍, 姜東, 曹衛(wèi)星, 戴廷波. 干旱脅迫對小麥幼苗根系生長和葉片光合作用的影響. 應用生態(tài)學報, 2012, 23(3): 724-730.

    MA F J, LI D D, CAI J, JIANG D, CAO W X, DAI T B. Effects of drought stress on root growth and leaf photosynthesis of wheat seedlings. Chinese Journal of Applied Ecology, 2012, 23(3): 724-730. (in Chinese)

    [3] 賀安娜, 姚奕. 虎耳草冬季凈光合速率、蒸騰速率日變化及其影響因子分析. 西南農業(yè)學報, 2011, 24(4): 1298-1302.

    HE A N, YAO Y. Analysis of diurnal variation of net photosynthetic rate, transpiration rate and its influencing factors in winter. SouthwestChina Journal of Agricultural Sciences, 2011, 24(4): 1298-1302. (in Chinese)

    [4] 張治安, 楊福, 陳展宇, 徐克章. 菰葉片凈光合速率日變化及其與環(huán)境因子的相互關系. 中國農業(yè)科學, 2006, 39(3): 502-509.

    ZHANG Z A, YANG F, CHEN Z Y, XU K Z. Diurnal variation of net photosynthetic rate in leaves of wild rice and its relationship with environmental factors. Scientia Agricultura Sinica, 2006, 39(3): 502-509. (in Chinese)

    [5] 唐曉培, 楊麗, 馮冬雪, 高壯壯, 張文杰, 劉海軍. 非充分灌溉下8個小麥品種旗葉光合與產量及水分利用效率的關系. 干旱地區(qū)農業(yè)研究, 2020, 38(4): 245-252+265.

    TANG X P, YANG L, FENG D X, GAO Z Z, ZHANG W J, LIU H J. The relationship between flag leaf photosynthesis and yield and water use efficiency of eight wheat cultivars under inadequate irrigation. Agricultural Research in the Arid Areas, 2020, 38(4): 245-252+265. (in Chinese)

    [6] 徐恒永, 趙君實. 高產冬小麥的冠層光合能力及不同器官的貢獻. 作物學報, 1995, 21(2): 204-209.

    XU H Y, ZHAO J S. Canopy photosynthetic capacity and contribution of different organs in high-yielding winter wheat. Acta Agronomica Sinica, 1995, 21(2): 204-209. (in Chinese)

    [7] LIU L P, SUN G L, REN X F, LI C D, SUN D F. Identification of QTL underlying physiological and morphological traits of flag leaf in barley. BMC Genetics, 2015, 16(1):29.

    [8] 熊淑萍, 吳克遠, 王小純, 張捷, 杜盼, 吳懿鑫, 馬新明. 不同氮效率基因型小麥根系吸收特性與氮素利用差異的分析. 中國農業(yè)科學, 2016, 49(12): 2267-2279.

    XIONG S P, WU K Y, WANG X C, ZHANG J, DU P, WU Y X, MA X M. Analysis of root absorption characteristics and nitrogen utilization of wheat genotypes with different n efficiency.Scientia Agricultura Sinica,2016, 49(12): 2267-2279. (in Chinese)

    [9] LIDIYA M, OREST F, VALERIY P, et al. Nitrogen balance of crop production in Ukraine. Journal of Environmental Management, 2019, 246: 860-867.

    [10] LASSALETTA L, BILLEN G, GARNIER J, BOUWMAN L, VELAZQUEZ E, MUELLER N D, GERBER J S. Nitrogen use in the global food system: Past trends and future trajectories of agronomic performance, pollution, trade, and dietary demand. Environmental Research Letters, 2016, 11(9):095007.

    [11] 蓋鈞鎰, 章元明, 王建康. QTL混合遺傳模型擴展至2對主基因+多基因時的多世代聯(lián)合分析. 作物學報, 2000, 26(4): 385-391.

    GAI J Y, ZHANG Y M, WANG J K. A joint analysis of multiple generations for QTL models extended to mixed two major genes plus polygene.Acta Agronomica Sinica, 2000, 26(4): 385-391. (in Chinese)

    [12] WANG J K, PODLICH D W, COOPER M, DELACY I H. Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits. Theoretical and Applied Genetics, 2001, 103: 804-816.

    [13] WANG J K, GAI J Y. Mixed inheritance model for resistance to agromyzid beanfly (Zehntner) in soybean. Euphytica, 2001, 122(1): 9-18.

    [14] GAI J Y. Segregation analysis on genetic system of quantitative traits in plants. Frontiers of Biology, 2006, 1(1): 85-92.

    [15] KHAN M I, KHATTAK G S S, KHAN A J, SUBHAN F, MOHAMMAD T, ALI A. Genetic control of flag leaf area in wheat () crosses. African Journal of Agricultural Research, 2012, 7(27): 3978-3990.

    [16] KUMAR S. Quantitative genetics, molecular markers, and plant improvement.Scholarly Journal of Agricultural Science, 2014, 4(10): 502-511.

    [17] BECHE E, BENIN G, SILVA C L, MUNARO L B, MARCHESE J A. Genetic gain in yield and changes associated with physiological traits in Brazilian wheat during the 20th century. European Journal of Agronomy, 2014, 61: 49-59.

    [18] 蓋鈞鎰, 章元明, 王建康. 植物數量性狀遺傳體系. 北京: 科學出版社, 2003: 96-102.

    GAI J Y, ZHANG Y M, WANG J K. Genetic system of quantitative traits in plants. Beijing: Science Press, 2003: 96-102. (in Chinese)

    [19] MA J, TU Y, ZHU J, LUO W, LIU H, LI C, LI S, LIU J J, DING P Y, AHSAN H. Flag leaf size and posture of bread wheat: genetic dissection, QTL validation and their relationships with yield-related traits. Theoretical and Applied Genetics, 2020, 133(1): 297-315.

    [20] FAROKHZADEH S, FAKHERI B A, NEZHAD N M, TAHMASEBI S, MIRSOLEIMANI A. Mapping QTLs of flag leaf morphological and physiological traits related to aluminum tolerance in wheat (L.). Physiology and Molecular Biology of Plants, 2019, 25(4): 975-990.

    [21] LIU K Y, XU H, LIU G, GUAN P F, ZHOU X Y, PENG H R, YAO Y Y, NI Z F, SUN Q X, DU J K. QTL mapping of flag leaf-related traits in wheat (L.).Theoretical and Applied Genetics, 2018, 131(4): 839-849.

    [22] FRANCKI M G, WAIKER E, LI D A, FORREST K. High-density SNP mapping reveals closely linked QTL for resistance toblotch (SNB) in flag leaf and glume of hexaploid wheat. Canadian Journal of Genetics and Cytology, 2018, 61(2): 145-149.

    [23] HUSSAIN W, BAENZIGER P S, BELAMKAR V, GUTTIERI M J, VENEGAS J P, EASTERLY A, SALLAM A, POLAND J. Genotyping-by-sequencing derived high-density linkage map and its application to QTL mapping of flag leaf traits in bread wheat. Scientific Reports, 2017, 7(1): 16394.

    [24] 趙倩茹, 鐘興華, 張飛, 房偉民, 陳發(fā)棣, 滕年軍. 切花小菊綠心性狀雜種優(yōu)勢與混合遺傳分析. 中國農業(yè)科學, 2018, 51(5): 964-976.

    ZHAO Q R, ZHONG X H, ZHANG F, FANG W M, CHEN F D, TENG N J. Heterosis and mixed genetic analysis of green-center trait of spray cut chrysanthemum. Scientia Agricultura Sinica, 2018, 51(5): 964-976. (in Chinese)

    [25] 葉紅霞, 呂律, 海睿, 胡雨晴, 汪炳良. 甜瓜果實糖含量的主基因+多基因遺傳分析. 浙江大學學報(農業(yè)與生命科學版), 2019, 45(4): 391-400.

    YE H X, Lü L, HAI R, HU Y Q, WANG B L. Genetic analysis of main gene + polygene in sugar content of melon fruit. Journal of Zhejiang University(Agriculture and Life Sciences), 2019, 45(4): 391-400. (in Chinese)

    [26] 曹齊衛(wèi), 張允楠, 王永強, 楊桂蘭, 孫小鐳, 李利斌. 黃瓜節(jié)間長的主基因+多基因混合遺傳模型分析. 農業(yè)生物技術學報, 2018, 26(2): 205-212.

    CAO Q W, ZHANG Y N, WANG Y Q, YANG G L, SUN X L, LI L B. Genetic analysis of internode length using mixed major-gene plus polygene inheritance model in. Journal of Agricultural Biotechnology, 2018, 26(2): 205-212. (in Chinese)

    [27] 汪文祥, 胡瓊, 梅德圣, 李云昌, 周日金, 王會成, 洪濤, 付麗, 劉佳. 甘藍型油菜分枝角度主基因+多基因混合遺傳模型及遺傳效應. 作物學報, 2016, 42(8): 1103-1111.

    WANG W X, HU Q, MEI D S, LI Y C, ZHOU R J, WANG H C, HONG T, FU L, LIU J. Genetic effects of branch angle using mixture model of major gene plus polygene in. Acta Agronomica Sinica, 2016, 42(8): 1103-1111. (in Chinese)

    [28] YE Y J, WU J Y, FENG L, JU Y Q, CAI M, CHENG T R, PAN H T, ZHANG Q X. Heritability and gene effects for plant architecture traits of crape myrtle using major gene plus polygene inheritance analysis.Scientia Horticulturae, 2017, 225: 335-342.

    [29] 江建華, 張武漢, 黨小景, 榮慧, 葉琴, 胡長敏, 張瑛, 何強, 王德正. 水稻核不育系柱頭性狀的主基因+多基因遺傳分析.作物學報, 2021, 47(7): 1215-1227.

    JIANG J H, ZHANG W H, DANG X J, RONG H, YE Q, HU C M, ZHANG Y, HE Q, WANG D Z. Genetic analysis of the stigma trait in rice nuclear male sterile lines. Acta Agronomica Sinica, 2021, 47(7): 1215-1227. (in Chinese)

    [30] 黃冰艷, 張新友, 苗利娟, 劉華, 秦利, 徐靜, 張忠信, 湯豐收, 董文召, 韓鎖義, 劉志勇. 花生油酸和亞油酸含量的遺傳模式分析. 中國農業(yè)科學, 2012, 45(4): 617-624.

    HUANG B Y, ZHANG X Y, MIAO L J, LIU H, QIN L, XU J, ZHANG Z X, TANG F S, DONG W Z, HAN S Y, LIU Z Y. Inheritance analysis of oleic acid and linoleic acid content of. Scientia Agricultura Sinica, 2012, 45(4): 617-624. (in Chinese)

    [31] 解松峰, 吉萬全, 王長有, 胡衛(wèi)國, 李俊, 張耀元, 師曉曦, 張俊杰, 張宏, 陳春環(huán). 小麥穗部性狀的主基因+多基因混合遺傳模型分析. 中國農業(yè)科學, 2019, 52(24): 4437-4452.

    XIE S F, JI W Q, WANG C Y, HU W G, LI J, ZHANG Y Y, SHI X X, ZHANG J J, ZHANG H, CHEN C H. Analysis of mixed genetic model of major gene and polygene for wheat panicle traits. Scientia Agricultura Sinica,2019, 52(24): 4437-4452. (in Chinese)

    [32] 解松峰, 吉萬全, 張耀元, 張俊杰, 胡衛(wèi)國, 李俊, 王長有, 張宏, 陳春環(huán). 小麥重要產量性狀的主基因+多基因混合遺傳分析. 作物學報, 2020, 46(3): 365-384.

    XIE S F, JI W Q, ZHANG Y Y, ZHANG J J, HU W G, LI J, WANG C Y, ZHANG H, CHEN C H. Genetic analysis of major genes and polygenes in important yield traits of wheat.Acta Agronomica Sinica, 2020, 46(3): 365-384. (in Chinese)

    [33] CUI F, FAN X L, ZHAO C H, ZHANG W, CHEN M, JI J, LI J M. A novel genetic map of wheat: utility for mapping QTL for yield under different nitrogen treatments. BMC Genetics, 2014, 15: 57-74.

    [34] FAN X L, CUI F, ZHAO C H, ZHANG W, YANG L J, ZHAO X Q, HAN J, SU Q N, JI J, ZHAO Z W, TONG Y P, LI J M. QTLs for flag leaf size and their influence on yield-related traits in wheat (L.). Molecular Breeding, 2015, 35(1): 1-16.

    [35] 胡燕美, 蘇慧, 朱玉磊, 李金鵬, 李金才, 宋有洪. 花后早期增溫對小麥旗葉光合和抗氧化特性及籽粒發(fā)育的影響. 麥類作物學報, 2020, 40(10): 1247-1256.

    HU Y M, SU H, ZHU Y L, LI J P, LI J C, SONG Y H. Effects of early post-anthesis warming on photosynthesis and antioxidant characteristics of flag leaves and grain development of wheat. Journal of Triticeae Genetics, 2020, 40(10): 1247-1256. (in Chinese)

    [36] 高娣, 吳宏亮, 劉根紅, 康建宏, 堅天才, 李鑫. 不同追氮時期對花后高溫脅迫下春小麥旗葉蛋白質和核酸代謝的影響. 農業(yè)科學研究, 2020, 41(3): 21-25.

    GAO D, WU H L, LIU G H, KANG J H, JIAN T C, LI X. Effects of nitrogen topdressing on protein and nucleic acid metabolism in flag leaves of spring wheat under post-anthesis heat stress. Research in Agricultural Sciences, 2020, 41(3): 21-25. (in Chinese)

    [37] 尹嘉德, 侯慧芝, 張緒成, 王紅麗, 于顯楓, 方彥杰, 馬一凡, 張國平, 雷康寧. 全膜覆土下施有機肥對春小麥旗葉碳氮比、光合特性和產量的影響. 應用生態(tài)學報, 2020, 31(11): 3749-3757.

    YIN J D, HOU H Z, ZHANG X C, WANG H L, YU X F, FANG Y J, MA Y F, ZHANG G P, LEI K N. Effects of organic fertilizer application on carbon nitrogen ratio, photosynthetic characteristics and yield of flag leaf of spring wheat. Chinese Journal of Applied Ecology, 2020, 31(11): 3749-3757. (in Chinese)

    [38] Khan M I, Khattak G S S, KHAN A J, KHAN A J, SUBHAN F, MOHAMMAD T, ALI A. Genetic control of flag leaf area in wheat () crosses.African Journal of Agricultural Research, 2012, 7(27): 3978-3990.

    [39] 鈕力亞, 王偉偉, 王偉, 王奉芝, 趙松山, 于亮. 小麥功能葉對產量及其構成因素的影響. 作物研究, 2018, 32(4): 295-298.

    NIU L Y, WANG W W, WANG W, WANG F Z, ZHAO S S, YU L. Effects of functional leaves on yield and its components in wheat. Crop Research, 2018, 32(4): 295-298. (in Chinese)

    [40] 王敏, 張從宇. 小麥旗葉性狀與產量因素的相關與回歸分析. 種子, 2004, 23(3): 17-21.

    WANG M, ZHANG C Y. Correlation and regression analysis of flag leaf traits and yield factors in wheat. Seed, 2004, 23(3): 17-21. (in Chinese)

    [41] 黃杰, 喬冀良, 苗運武, 張振永, 葛昌斌, 廖平安. 小麥產量與旗葉性狀的相關性分析. 中國種業(yè), 2018(1): 63-64.

    HUANG J, QIAO J L, MIAO Y W, ZHANG Z Y, GE C B, LIAO P A. Correlation analysis of wheat yield and flag leaf traits. China Seed Industry, 2018(1): 63-64.(in Chinese)

    [42] 姜莉莉, 王東. 不同年代小麥產量性狀與農藝性狀的相關性比較. 華北農學報, 2017, 32(S1): 130-134.

    JIANG L L, WANG D. Comparison of yield and agronomic traits of wheat in different ages. Acta Agriculturae Boreali-Sinica, 2017, 32(S1): 130-134. (in Chinese)

    Genetic Analysis of Flag Leaf Traits in Wheat Under High and Low Nitrogen

    WU YaRui1, LIU XiJian1, YANG GuoMin2, LIU HongWei1, KONG WenChao1, WU YongZhen1, SUN Han1, QIN Ran1, CUI Fa1, ZHAO ChunHua1

    1College of Agriculture, Ludong University, Yantai 264025, Shandong;2Soil and Fertilizer Workstation in Heze City, Heze 274000, Shandong

    【】Flag leaf is an important place for wheat photosynthetic carbon fixation,which plays an important role in wheat yield. The genetic characteristics and the genetic mechanism were analyzed under high and low nitrogen for flag leaf traits of wheat, which will provide a reference for excellent plant-type breeding and high-yield breeding. 【】188 recombinant inbred line (RIL) populations derived from a cross between Kenong9204 and Jing411 was used in this study, which were planted in low nitrogen (LN) and high nitrogen (HN), respectively. The flag leaf traits of 188 RILs were investigated in 6 different environments, then the genetic analysis was conducted to determine the number of genes controlling each trait, and to estimate the genetic effect value and the heritability. In addition, the relationship between flag leaf characters and yield related traits of wheat was also studied.【】Under LN environment: the optimal genetic model of flag leaf length was 2MG-CE (two pairs of interaction major genes) in E3. The additive × additive epistatic interaction value was 1.098, and the heritability of major genes was 31.35%. The flag leaf length was polygenic in another LN environment. The width of flag leaf was polygenic in all the LN environment. The optimal genetic model for flag leaf area (except E5) was 2MG-CE. The additive × additive epistatic interaction value was 1.884 and the heritability of major genes was 36.7%, while it was polygenic inheritance in E5. Under HN environment: The optimal genetic model for flag leaf length (except E4) was 2MG-CE, the additive × additive epistatic interaction value was 1.133, and the heritability of major genes was 32.6%. The optimal genetic model was 2MG-ER (two pairs of recessive epistatic major genes) in E4, which the additive effect value was 1.431 and 1.108 for the first and the second major genes respectively, and the heritability of the major gene was 51.77%. The optimal genetic model for flag leaf width (except E2) was 2MG-CE, the additive × additive epistatic interaction value was 0.119, and the heritability of major genes was 37.29%, while it showed polygenic inheritance in E2. The optimal genetic model for flag leaf area was 2MG-CE, which the additive × additive epistatic interaction value was 3.067 and the heritability of the main gene was 44.42%. The genetic models of flag leaf traits were different in different environments, which the genetic model was more stable under HN than that in LN. The correlation analysis of flag leaf and yield traits showed that flag leaf traits were significantly positively correlated with grain number per spike, grain weight per spike and yield per plant, and the influence degree was different in the 6 environments. 【】Flag leaf traits are easily affected by environment, and the performance of flag leaf traits is different in HN and LN. Flag leaf traits exhibited different major gene inheritance and polygene inheritance in LN, while they showed major gene inheritance which controlled by two pairs of interactions genes in most of HN environment, which might be major QTLs. Yield per plant and grain weight per spike could be increased by improving flag leaf traits.

    wheat; flag leaf traits; genetic model; production; correlation analysis

    10.3864/j.issn.0578-1752.2022.01.001

    2021-06-07;

    2021-08-16

    山東省農業(yè)良種工程(2019LZGC016)、山東省自然科學基金重大基礎研究(ZR2019ZD16)、山東省高等學校青創(chuàng)科技計劃(2019KJF002)、煙臺市新舊動能研究院(2019XJDN007)

    武亞瑞,E-mail:1684813375@qq.com。通信作者崔法,E-mail:sdaucf@126.com。通信作者趙春華,E-mail:sdauzch@126.com

    (責任編輯 李莉)

    猜你喜歡
    加性旗葉葉面積
    不同品種小麥灌漿期旗葉光合特性及光合基因表達對臭氧濃度升高的響應
    作物學報(2022年9期)2022-07-18 00:58:42
    ?2?4[u]-加性循環(huán)碼
    水分虧缺對小麥芒和旗葉光合特性及蔗糖、淀粉合成的影響
    旗葉衰老產量差異 可作小麥優(yōu)選依據
    作物葉面積測量的研究進展
    馬奶子葡萄葉面積評估模型的建立
    抽穗后不同時期去除旗葉對不同穗型小麥產量的影響
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關系:超可加性、次可加性或不可加性
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關系:超可加性、次可加性或不可加性
    “光合作用與細胞呼吸”重點分析和突破
    两性夫妻黄色片| 搡老乐熟女国产| 91国产中文字幕| 久久中文字幕人妻熟女| 久久国产精品男人的天堂亚洲| 亚洲av片天天在线观看| 欧美黄色淫秽网站| 亚洲av电影在线进入| 久久狼人影院| 亚洲成人免费电影在线观看| 香蕉久久夜色| 50天的宝宝边吃奶边哭怎么回事| 热99re8久久精品国产| 性色av乱码一区二区三区2| 亚洲成人国产一区在线观看| av在线播放免费不卡| 日韩中文字幕欧美一区二区| 亚洲午夜理论影院| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩一级在线毛片| 精品无人区乱码1区二区| 桃红色精品国产亚洲av| 中国美女看黄片| 久久精品熟女亚洲av麻豆精品| 欧美一级毛片孕妇| 国产午夜精品久久久久久| 亚洲 欧美一区二区三区| 久久精品国产亚洲av高清一级| 麻豆乱淫一区二区| 久久精品亚洲av国产电影网| 好男人电影高清在线观看| 国产单亲对白刺激| 精品久久久精品久久久| 亚洲七黄色美女视频| 午夜福利乱码中文字幕| 自线自在国产av| 国产区一区二久久| 国产主播在线观看一区二区| 老熟女久久久| 黄网站色视频无遮挡免费观看| 人人妻人人添人人爽欧美一区卜| 久久这里只有精品19| 国产精品久久久久久人妻精品电影| 午夜福利一区二区在线看| 大型黄色视频在线免费观看| 久久精品亚洲精品国产色婷小说| 国产aⅴ精品一区二区三区波| 国产国语露脸激情在线看| 免费人成视频x8x8入口观看| 夜夜爽天天搞| 国产精华一区二区三区| 日韩视频一区二区在线观看| 黄色女人牲交| 老熟女久久久| 久热爱精品视频在线9| 亚洲午夜理论影院| 亚洲情色 制服丝袜| 免费在线观看黄色视频的| 最近最新中文字幕大全免费视频| 久久中文看片网| 在线视频色国产色| 不卡一级毛片| 亚洲国产精品合色在线| 脱女人内裤的视频| 夜夜躁狠狠躁天天躁| 国产在线精品亚洲第一网站| 欧美成人免费av一区二区三区 | 在线视频色国产色| 婷婷成人精品国产| 国产视频一区二区在线看| 国产精品免费视频内射| 久久久精品区二区三区| 9热在线视频观看99| 人妻丰满熟妇av一区二区三区 | 婷婷丁香在线五月| 久久九九热精品免费| 亚洲av熟女| 成人免费观看视频高清| 久久天堂一区二区三区四区| 国产在线观看jvid| 日韩中文字幕欧美一区二区| 丁香六月欧美| 91老司机精品| 老司机午夜福利在线观看视频| 香蕉国产在线看| 老鸭窝网址在线观看| 国产成人精品无人区| 黄频高清免费视频| 欧美激情高清一区二区三区| 久久99一区二区三区| 最近最新中文字幕大全免费视频| 91av网站免费观看| 色精品久久人妻99蜜桃| 亚洲欧美日韩高清在线视频| 91字幕亚洲| 香蕉丝袜av| 国产麻豆69| 一本综合久久免费| 国产真人三级小视频在线观看| а√天堂www在线а√下载 | 一a级毛片在线观看| 国产精品 国内视频| 国产精品永久免费网站| ponron亚洲| av有码第一页| 亚洲综合色网址| 手机成人av网站| av视频免费观看在线观看| 18在线观看网站| 亚洲欧美日韩高清在线视频| 无遮挡黄片免费观看| 51午夜福利影视在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜免费鲁丝| 欧美日韩黄片免| 精品人妻在线不人妻| 久久国产精品大桥未久av| 女人被狂操c到高潮| 在线十欧美十亚洲十日本专区| 国产又爽黄色视频| 亚洲五月婷婷丁香| 在线观看舔阴道视频| 亚洲欧美色中文字幕在线| 人成视频在线观看免费观看| 久久午夜综合久久蜜桃| 少妇粗大呻吟视频| 亚洲欧美一区二区三区黑人| 成人影院久久| 人人妻,人人澡人人爽秒播| 欧美日韩乱码在线| 精品一区二区三区av网在线观看| 999久久久国产精品视频| 天天添夜夜摸| 成人三级做爰电影| 热re99久久精品国产66热6| 欧美激情极品国产一区二区三区| 国产麻豆69| 免费一级毛片在线播放高清视频 | 首页视频小说图片口味搜索| 亚洲一区中文字幕在线| 亚洲熟妇中文字幕五十中出 | 黄色怎么调成土黄色| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 黄片播放在线免费| 日韩欧美一区视频在线观看| 亚洲熟女精品中文字幕| 老司机亚洲免费影院| 国产免费av片在线观看野外av| 日韩精品免费视频一区二区三区| 久9热在线精品视频| 亚洲午夜理论影院| 亚洲,欧美精品.| 亚洲人成电影免费在线| 激情在线观看视频在线高清 | 国产aⅴ精品一区二区三区波| 啦啦啦在线免费观看视频4| 韩国av一区二区三区四区| 亚洲国产精品合色在线| 悠悠久久av| 黄色女人牲交| а√天堂www在线а√下载 | 国产一区二区三区视频了| 国产又色又爽无遮挡免费看| 亚洲精品美女久久av网站| 国产男女超爽视频在线观看| 12—13女人毛片做爰片一| 亚洲国产精品合色在线| 黄色女人牲交| 国产精品久久视频播放| 国产aⅴ精品一区二区三区波| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 1024视频免费在线观看| 天天躁日日躁夜夜躁夜夜| 精品一品国产午夜福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 咕卡用的链子| 午夜福利一区二区在线看| 久久久国产成人免费| 精品久久久精品久久久| 国产av一区二区精品久久| 极品教师在线免费播放| 另类亚洲欧美激情| 国产精品久久久久成人av| 亚洲免费av在线视频| 亚洲在线自拍视频| 激情在线观看视频在线高清 | 日本黄色日本黄色录像| 免费看十八禁软件| 人妻久久中文字幕网| 好看av亚洲va欧美ⅴa在| ponron亚洲| 中文字幕人妻丝袜一区二区| 视频在线观看一区二区三区| 久久草成人影院| 亚洲中文日韩欧美视频| 亚洲精品国产区一区二| 99精品欧美一区二区三区四区| 中国美女看黄片| 久久久久精品国产欧美久久久| 国产精品国产高清国产av | 91精品国产国语对白视频| 亚洲成人免费av在线播放| 一级片免费观看大全| 一级毛片精品| 国产又色又爽无遮挡免费看| 很黄的视频免费| 老司机在亚洲福利影院| 天天添夜夜摸| 国产精品一区二区在线不卡| 国产日韩欧美亚洲二区| 别揉我奶头~嗯~啊~动态视频| 午夜福利影视在线免费观看| 在线免费观看的www视频| 日本vs欧美在线观看视频| 欧美不卡视频在线免费观看 | 啪啪无遮挡十八禁网站| 啦啦啦视频在线资源免费观看| 波多野结衣av一区二区av| 亚洲精品自拍成人| 在线观看舔阴道视频| 久久天堂一区二区三区四区| 99久久综合精品五月天人人| 日韩制服丝袜自拍偷拍| 成人手机av| 国产精品久久久久成人av| 亚洲av美国av| 多毛熟女@视频| 韩国精品一区二区三区| 制服诱惑二区| 欧美不卡视频在线免费观看 | 欧美色视频一区免费| 热99国产精品久久久久久7| 99re6热这里在线精品视频| 国产麻豆69| 欧美成人午夜精品| 欧美日韩亚洲综合一区二区三区_| 色尼玛亚洲综合影院| 韩国精品一区二区三区| 亚洲熟女精品中文字幕| 亚洲精品成人av观看孕妇| 老司机影院毛片| av天堂久久9| 国产亚洲精品第一综合不卡| 两性夫妻黄色片| 在线十欧美十亚洲十日本专区| 亚洲久久久国产精品| 色综合欧美亚洲国产小说| 在线观看www视频免费| 51午夜福利影视在线观看| 免费观看精品视频网站| 国产不卡av网站在线观看| 老司机午夜十八禁免费视频| 午夜视频精品福利| 久久性视频一级片| 一级黄色大片毛片| 精品福利永久在线观看| 超碰成人久久| 一个人免费在线观看的高清视频| 看免费av毛片| 日韩欧美三级三区| 精品午夜福利视频在线观看一区| 日本五十路高清| 欧美乱妇无乱码| 亚洲 欧美一区二区三区| 悠悠久久av| 亚洲第一青青草原| 欧美成狂野欧美在线观看| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲| 一级,二级,三级黄色视频| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| 90打野战视频偷拍视频| 亚洲第一欧美日韩一区二区三区| 精品福利观看| 国产欧美日韩精品亚洲av| 国产亚洲欧美98| 99国产精品一区二区三区| 久久天堂一区二区三区四区| 亚洲精华国产精华精| 欧美最黄视频在线播放免费 | 国产精品乱码一区二三区的特点 | 中文字幕高清在线视频| 亚洲成人免费av在线播放| 亚洲七黄色美女视频| 高潮久久久久久久久久久不卡| 亚洲专区中文字幕在线| 亚洲国产精品一区二区三区在线| 自拍欧美九色日韩亚洲蝌蚪91| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 久久99一区二区三区| 国产精品欧美亚洲77777| 亚洲专区字幕在线| 新久久久久国产一级毛片| 国产精品自产拍在线观看55亚洲 | 日本黄色视频三级网站网址 | 黑人猛操日本美女一级片| 曰老女人黄片| 他把我摸到了高潮在线观看| 久久精品亚洲av国产电影网| 精品一区二区三区av网在线观看| 不卡一级毛片| av有码第一页| 久久久精品国产亚洲av高清涩受| 精品高清国产在线一区| 麻豆乱淫一区二区| 亚洲,欧美精品.| 亚洲国产精品一区二区三区在线| 精品久久蜜臀av无| 手机成人av网站| 国产不卡一卡二| 视频在线观看一区二区三区| 久久久久久久久免费视频了| 超碰97精品在线观看| 国产91精品成人一区二区三区| 亚洲av成人不卡在线观看播放网| 久久精品人人爽人人爽视色| 性少妇av在线| 在线观看舔阴道视频| 久久狼人影院| 在线观看午夜福利视频| 黄色丝袜av网址大全| 免费在线观看影片大全网站| 日韩欧美三级三区| 午夜精品在线福利| 少妇猛男粗大的猛烈进出视频| 悠悠久久av| 久久人妻福利社区极品人妻图片| 大码成人一级视频| av免费在线观看网站| 免费在线观看完整版高清| 欧美乱妇无乱码| 美女福利国产在线| 涩涩av久久男人的天堂| 亚洲七黄色美女视频| 国产精品免费大片| 国产精品一区二区精品视频观看| 国产不卡av网站在线观看| 别揉我奶头~嗯~啊~动态视频| netflix在线观看网站| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 国产精品一区二区在线不卡| 老司机靠b影院| 欧美日韩亚洲综合一区二区三区_| 亚洲成a人片在线一区二区| 日日夜夜操网爽| av中文乱码字幕在线| 久久久久久久精品吃奶| 黄色丝袜av网址大全| 女人被躁到高潮嗷嗷叫费观| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 国产精品1区2区在线观看. | 麻豆av在线久日| 亚洲 欧美一区二区三区| 无遮挡黄片免费观看| 久久精品91无色码中文字幕| 黄片播放在线免费| av线在线观看网站| 一进一出抽搐动态| 亚洲伊人色综图| 久久这里只有精品19| 国产精品综合久久久久久久免费 | 另类亚洲欧美激情| 国产高清国产精品国产三级| 国产精品免费一区二区三区在线 | 国产男靠女视频免费网站| 视频区图区小说| 午夜视频精品福利| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 一级作爱视频免费观看| 嫩草影视91久久| 在线观看免费午夜福利视频| 亚洲av成人不卡在线观看播放网| 亚洲欧美激情综合另类| 亚洲三区欧美一区| 少妇裸体淫交视频免费看高清 | 天天影视国产精品| 久久精品成人免费网站| xxx96com| 女人高潮潮喷娇喘18禁视频| a级毛片在线看网站| 一级毛片女人18水好多| 国产成人系列免费观看| 久久性视频一级片| 无限看片的www在线观看| 亚洲第一青青草原| 首页视频小说图片口味搜索| 精品一区二区三区视频在线观看免费 | 91国产中文字幕| 91在线观看av| 国产91精品成人一区二区三区| 极品人妻少妇av视频| 国产精品1区2区在线观看. | 99久久精品国产亚洲精品| 午夜两性在线视频| 国产有黄有色有爽视频| 中出人妻视频一区二区| 亚洲精品在线美女| 老司机影院毛片| 国产片内射在线| 久久久久久久国产电影| 性色av乱码一区二区三区2| 欧洲精品卡2卡3卡4卡5卡区| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 午夜成年电影在线免费观看| 18禁观看日本| 亚洲av第一区精品v没综合| 精品国产乱子伦一区二区三区| 波多野结衣av一区二区av| 嫁个100分男人电影在线观看| 91大片在线观看| 国产亚洲欧美精品永久| 亚洲精品在线美女| 叶爱在线成人免费视频播放| 久久热在线av| 女人高潮潮喷娇喘18禁视频| 国产精品久久久久久人妻精品电影| 精品久久久久久久久久免费视频 | 性少妇av在线| 国产亚洲精品久久久久久毛片 | 制服人妻中文乱码| 国产精品影院久久| 久久久久国内视频| 香蕉久久夜色| 国产精品二区激情视频| 少妇的丰满在线观看| 欧美日韩亚洲高清精品| 1024视频免费在线观看| 咕卡用的链子| 老汉色av国产亚洲站长工具| 校园春色视频在线观看| 久久精品国产a三级三级三级| 国产精品久久视频播放| 日韩免费高清中文字幕av| 99在线人妻在线中文字幕 | 一区二区日韩欧美中文字幕| 很黄的视频免费| 中文字幕制服av| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 香蕉丝袜av| 国产1区2区3区精品| 成年动漫av网址| 国产国语露脸激情在线看| 女性被躁到高潮视频| 亚洲熟妇中文字幕五十中出 | 99久久精品国产亚洲精品| 久久香蕉精品热| 在线看a的网站| 国产成人精品久久二区二区免费| 他把我摸到了高潮在线观看| 国产97色在线日韩免费| 精品福利永久在线观看| av一本久久久久| 亚洲国产毛片av蜜桃av| 男人操女人黄网站| 国产精品99久久99久久久不卡| 成年人免费黄色播放视频| 两人在一起打扑克的视频| 亚洲专区国产一区二区| а√天堂www在线а√下载 | 国产精品二区激情视频| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 亚洲欧洲精品一区二区精品久久久| 99精品欧美一区二区三区四区| 免费一级毛片在线播放高清视频 | 亚洲av成人av| 精品国产国语对白av| 亚洲一码二码三码区别大吗| 在线观看www视频免费| 无人区码免费观看不卡| 久久狼人影院| 亚洲精品国产一区二区精华液| 女人高潮潮喷娇喘18禁视频| 别揉我奶头~嗯~啊~动态视频| 在线观看66精品国产| 欧美激情高清一区二区三区| 人妻丰满熟妇av一区二区三区 | 亚洲男人天堂网一区| e午夜精品久久久久久久| 国产在视频线精品| 在线天堂中文资源库| 久久精品国产99精品国产亚洲性色 | 精品高清国产在线一区| 涩涩av久久男人的天堂| 99精品久久久久人妻精品| 不卡av一区二区三区| 日韩欧美三级三区| 久久久国产一区二区| 久久精品国产a三级三级三级| 别揉我奶头~嗯~啊~动态视频| 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 黄片小视频在线播放| 免费一级毛片在线播放高清视频 | 久久中文看片网| 国产aⅴ精品一区二区三区波| 日本精品一区二区三区蜜桃| 日韩免费高清中文字幕av| 精品免费久久久久久久清纯 | 国产亚洲精品久久久久5区| 午夜激情av网站| 欧美人与性动交α欧美精品济南到| 久久久久久久国产电影| 久久久国产一区二区| 免费在线观看完整版高清| 久久影院123| 国产主播在线观看一区二区| xxxhd国产人妻xxx| 51午夜福利影视在线观看| 高清欧美精品videossex| 国内毛片毛片毛片毛片毛片| 欧美日韩av久久| tocl精华| 天天影视国产精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲九九香蕉| 日韩中文字幕欧美一区二区| 日本vs欧美在线观看视频| 国产蜜桃级精品一区二区三区 | 91麻豆精品激情在线观看国产 | 校园春色视频在线观看| 99久久人妻综合| 亚洲精品久久成人aⅴ小说| 亚洲精品乱久久久久久| 国产在线观看jvid| 亚洲人成电影免费在线| 男人舔女人的私密视频| 精品国产一区二区久久| 村上凉子中文字幕在线| 亚洲免费av在线视频| 日本欧美视频一区| 精品久久蜜臀av无| 精品人妻在线不人妻| 久久影院123| 亚洲成av片中文字幕在线观看| 久久ye,这里只有精品| 国产91精品成人一区二区三区| 黑人欧美特级aaaaaa片| 亚洲熟妇中文字幕五十中出 | 成人18禁高潮啪啪吃奶动态图| 欧美午夜高清在线| 伦理电影免费视频| 丰满的人妻完整版| 脱女人内裤的视频| 精品久久久久久久毛片微露脸| 欧美人与性动交α欧美软件| a级毛片在线看网站| 亚洲精品国产区一区二| 免费在线观看黄色视频的| 十八禁高潮呻吟视频| 狂野欧美激情性xxxx| 精品人妻在线不人妻| 高潮久久久久久久久久久不卡| 久久久国产精品麻豆| 国产亚洲精品一区二区www | 99国产精品一区二区三区| 最近最新免费中文字幕在线| 热re99久久精品国产66热6| 国产极品粉嫩免费观看在线| 精品亚洲成国产av| 国产成人啪精品午夜网站| 久久中文看片网| 午夜日韩欧美国产| 免费人成视频x8x8入口观看| 国产亚洲精品久久久久久毛片 | 飞空精品影院首页| 国产区一区二久久| 色尼玛亚洲综合影院| 亚洲精品在线观看二区| 亚洲综合色网址| 80岁老熟妇乱子伦牲交| 精品一区二区三区av网在线观看| 国产欧美亚洲国产| 午夜福利影视在线免费观看| 亚洲黑人精品在线| 国产一区二区三区视频了| 极品人妻少妇av视频| 精品国产一区二区三区四区第35| 正在播放国产对白刺激| 女人久久www免费人成看片| 好看av亚洲va欧美ⅴa在| 久久久久久久久久久久大奶| 久久国产精品影院| 久9热在线精品视频| 亚洲成人免费av在线播放| 国产精品免费大片| 日韩欧美一区二区三区在线观看 | 久久久久久人人人人人| e午夜精品久久久久久久| 国产又色又爽无遮挡免费看| 亚洲欧洲精品一区二区精品久久久| 一本一本久久a久久精品综合妖精| 丰满人妻熟妇乱又伦精品不卡| 欧美黄色片欧美黄色片| 亚洲中文日韩欧美视频| 国产麻豆69| 久久久久久人人人人人| 免费看十八禁软件| 欧美日韩精品网址| 一夜夜www| 久久精品亚洲熟妇少妇任你| 国产亚洲精品久久久久5区| av电影中文网址| 99精国产麻豆久久婷婷| 亚洲国产欧美网| 天堂√8在线中文|