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

    水稻核不育系柱頭性狀的主基因+多基因遺傳分析

    2021-04-27 01:03:08江建華張武漢黨小景胡長敏王德正
    作物學(xué)報(bào) 2021年7期
    關(guān)鍵詞:加性柱頭顯性

    江建華 張武漢 黨小景 榮 慧 葉 琴 胡長敏 張 瑛 何 強(qiáng) 王德正,*

    水稻核不育系柱頭性狀的主基因+多基因遺傳分析

    江建華1,**張武漢2,**黨小景1榮 慧3葉 琴3胡長敏1張 瑛1何 強(qiáng)2,*王德正1,*

    1安徽省農(nóng)業(yè)科學(xué)院水稻研究所/ 水稻遺傳育種安徽省重點(diǎn)實(shí)驗(yàn)室, 安徽合肥 230031;2雜交水稻國家重點(diǎn)實(shí)驗(yàn)室 / 湖南雜交水稻研究中心, 湖南長沙 410125;3安徽科技學(xué)院, 安徽鳳陽 233100

    柱頭性狀是影響水稻不育系異交繁殖和雜交水稻制種產(chǎn)量的重要性狀。為創(chuàng)制長柱頭、高外露率的水稻溫敏核不育系提供遺傳信息, 調(diào)查了短柱頭、低外露率的粳型光溫敏核不育系7001S和長柱頭、高外露率的溫敏核不育系紫泰S及其雜交、自交獲得的F1、F2群體(350個(gè)株系)和F2:3群體(320個(gè)株系)的4個(gè)柱頭性狀, 分析了4個(gè)性狀之間的相關(guān)性, 并運(yùn)用主基因+多基因混合遺傳模型, 對(duì)2個(gè)世代4個(gè)性狀進(jìn)行了遺傳分析。結(jié)果表明, 4個(gè)柱頭性狀間均表現(xiàn)出極顯著正相關(guān), 相關(guān)系數(shù)介于0.262和0.895之間。柱頭長度、花柱長度、柱頭和花柱總長度(以下簡稱柱花總長度)均表現(xiàn)出受2對(duì)主效基因和微效基因共同控制, 除F2群體中柱花總長度的2對(duì)主基因表現(xiàn)為等加性效應(yīng)和等顯性效應(yīng)外, 其余均表現(xiàn)為加性-顯性-上位性效應(yīng), 3個(gè)性狀均表現(xiàn)出以主基因間的上位性效應(yīng)為主; F2群體柱頭外露率受2對(duì)加性-顯性-上位性主基因+多基因控制, 而F2:3群體則表現(xiàn)為受1對(duì)加性-顯性主基因+多基因控制, 以主基因間的加性效應(yīng)為主。2個(gè)世代中的4個(gè)柱頭性狀均以主基因遺傳為主。

    水稻; 核不育系; 柱頭性狀; 主基因+多基因模型; 遺傳分析

    我國65%以上人口以稻米為主食, 在人口持續(xù)增加而耕地卻以數(shù)萬公頃逐年減少的情況下, 大幅提高單位面積產(chǎn)量是保障糧食安全的必然選擇[1-2]。始于20世紀(jì)70年代的雜交水稻較同熟期常規(guī)稻增產(chǎn)20%, 其商業(yè)化種植使水稻生產(chǎn)力有了大幅提高[3]。我國雜交水稻年種植面積約1670萬公頃, 年制種面積高達(dá)18萬公頃[4]。雖然雜交水稻的制種產(chǎn)量從1976年的0.27 t hm–2增加到現(xiàn)在的2.7 t hm–2, 部分高產(chǎn)田塊達(dá)4 t hm–2 [5], 但隨著城市化進(jìn)程加快,勞動(dòng)力短缺, 制種成本不斷增加。提高制種產(chǎn)量是降低雜交水稻成本、增加種植效益的決定性因素。研究表明, 阻礙雜交水稻制種產(chǎn)量增加的主要因素是不育系異交結(jié)實(shí)率低[6]。不育系柱頭外露是限制雜交水稻制種異交結(jié)實(shí)率的關(guān)鍵因素[7]。代貴金等[8]認(rèn)為柱頭外露率和制種產(chǎn)量成極顯著正相關(guān); 田大成等[9]發(fā)現(xiàn)水稻閉穎后外露柱頭穎花平均結(jié)實(shí)率為64%, 遠(yuǎn)高于非外露穎花的15.7%; 李慶榮[10]發(fā)現(xiàn)柱頭外露穎花的結(jié)實(shí)率占制種異交結(jié)實(shí)率的70%~80%; 楊保漢[11]發(fā)現(xiàn)自然條件下不育系柱頭外露率每提高1個(gè)百分點(diǎn), 其結(jié)實(shí)率可提高0.74~0.92個(gè)百分點(diǎn), 每公頃制種產(chǎn)量可增加47~68 kg。可見, 水稻異交結(jié)實(shí)率與其柱頭性狀密切相關(guān)。

    水稻柱頭性狀主要包括柱頭長度、花柱長度、柱花總長度以及柱頭外露率。柱花總長度與柱頭外露率呈顯著[12-14]或極顯著[15-16]正相關(guān), 柱頭長度、花柱長度以及柱花總長度間也表現(xiàn)出極顯著正相關(guān)[16-17]。相比較于水稻其他性狀, 作為水稻重要異交性狀的柱頭性狀, 其研究進(jìn)展是嚴(yán)重滯后的。當(dāng)前對(duì)于柱頭性狀的研究主要集中于柱頭外露率, 研究者利用F2群體[12,18-21]、RIL群體[22-27]和自然群體[16,28]等在水稻12條染色體上共檢測(cè)到112個(gè)控制該性狀的QTL, 并將[18]、[29]和[30]分別精細(xì)定位于28.4 kb、322.9 kb和350.7 kb的區(qū)間。而關(guān)于柱頭長度性狀的研究則較少, 研究者利用F2群體[31]、RIL群體[22,31]、BIL群體、DH群體[31]以及自然群體[16-17,28,32]等在水稻的11條染色體(除11號(hào)染色體外)上檢測(cè)到28個(gè)控制STL的QTL, 在水稻1號(hào)、2號(hào)、3號(hào)、4號(hào)、6號(hào)、7號(hào)、10號(hào)和11號(hào)染色體上檢測(cè)到21個(gè)控制SYL的QTL[16-17,22,31-32], 在水稻的9條染色體(除8號(hào)、9號(hào)和10號(hào)染色體外)上檢測(cè)到16個(gè)控制TSSL的QTL[15,17,32-34], 但僅Liu等[34]將精細(xì)定位于水稻3號(hào)染色體短臂19.8 kb區(qū)間。說明當(dāng)前對(duì)于柱頭性狀的研究較少且未深入。這也可能是導(dǎo)致至今未見利用相關(guān)研究結(jié)果輔助選育長柱頭、高外露率不育系等育種材料報(bào)道的原因。

    本研究以課題組自行選育的短柱頭、低外露率的粳型光溫敏核不育系7001S[35]為母本, 長柱頭、高外露率的秈型溫敏核不育系紫泰S為父本, 在陵水進(jìn)行雜交和自交, 以獲得的7001S/紫泰S組合F1、F2和F2:3構(gòu)成的遺傳群體為研究對(duì)象, 利用廣泛應(yīng)用于水稻[36]、油菜[37]、大豆[38]、番茄[39]、小麥[40]等各種農(nóng)作物的主基因+多基因混合模型[41]對(duì)2個(gè)世代的4個(gè)柱頭性狀進(jìn)行遺傳分析, 同時(shí), 對(duì)2個(gè)世代(不同年份) 4個(gè)柱頭性狀間的相關(guān)性進(jìn)行分析, 旨在揭示柱頭性狀的遺傳基礎(chǔ)及其相關(guān)性, 以期為水稻柱頭性狀的QTL定位和全基因組關(guān)聯(lián)分析提供互補(bǔ)驗(yàn)證, 同時(shí)也為長柱頭、高柱頭外露率的水稻溫敏核不育系的創(chuàng)制提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    無色短柱頭、低柱頭外露率的粳型光溫敏核不育系7001S (P1) (圖1-A和圖2-A)、紫色長柱頭、高柱頭外露率的秈型溫敏核不育系紫泰S (P2) (圖1-B和圖2)及其雜交、自交所衍生的7001S/紫泰S組合F1、F2和F2:3構(gòu)成的遺傳群體材料。

    1.2 試驗(yàn)方法

    2016年冬季在海南陵水以7001S (P1)為母本, 紫泰S (P2)為父本, 雜交獲得7001S/紫泰S F1, 2017年正季將F1種植于安徽省農(nóng)業(yè)科學(xué)院水稻研究所廬江基地(31o30¢23¢¢N, 117o18¢9¢¢E), 自交獲得F2種子, 2018年正季將P1、P2、F1和F2種植于廬江基地(E1), 5月12日播種, 6月17日移栽。2019年正季將P1、P2、F1和F2:3(由2018年F2各單株收2粒種子混合而成)種植于廬江基地(E2), 5月13日播種, 6月14日移栽。P1、P2和F1分別種植4行, F2和F2:3群體均種植50行, 每行9株, 株行距16.7 cm × 26.7 cm, 常規(guī)栽培管理。

    柱頭長度(stigma length, STL)是指水稻雌蕊上毛刷狀部分的長度, 花柱長度(style length, SYL)是指毛刷狀底部與子房交界處之間的長度[16], 柱花總長度(the sum of stigma and style length, TSSL)是指柱頭長度和花柱長度之和(圖2-B)。

    柱頭長度相關(guān)性狀的測(cè)量方法: 植株盛花期從主莖穗(最高穗)上取10個(gè)正張開的穎花(每天取樣時(shí)間為10:00—14:00); 從穎花中剝?nèi)?個(gè)完整的雌蕊, 置于帶有攝像頭(MDX4)的體式顯微鏡(MZ11, 廣州市明美光電技術(shù)有限公司)視野中拍照; 利用圖像處理系統(tǒng)(MicroShot v1.2)在圖片上測(cè)量柱頭長度相關(guān)性狀(圖2-A, B); 計(jì)算單株所有16個(gè)柱頭(每個(gè)雌蕊2個(gè)柱頭)長度相關(guān)性狀的平均值, 作為該單株(株系)柱頭長度相關(guān)性狀的表型值。P1、P2和F1各株系調(diào)查8個(gè)雌蕊, 即16個(gè)柱頭長度相關(guān)性狀的數(shù)據(jù)作為該株系表型值; F2和F2:3分別調(diào)查了350個(gè)和320個(gè)單株, 每單株調(diào)查16個(gè)柱頭長度相關(guān)性狀的數(shù)據(jù)作為該單株表型值。

    A: 7001S和紫泰S的柱頭長度比較; B: 水稻雌蕊各部位名稱; C: 水稻柱頭無外露、單外露以及雙外露的表現(xiàn); D: 紫泰S田間柱頭外露表現(xiàn)。

    A: stigma morphology of 7001S and Zitai S; B: names of rice pistil parts defined in this study; C: phenotypes of single, dual, and no stigma exsertion in a spikelet; D: phenotype of exserted stigma in Zitai S. STL: stigma length; SYL: style length; TSSL: the sum of stigma and style length.

    柱頭外露率(percentage of exserted stigma, PES)是在家系開花結(jié)束后10 d調(diào)查, 具體方法參照Miyata等[12]略作改動(dòng), 簡言之, 柱頭外露率等于主莖穗的外露柱頭數(shù)除以總柱頭數(shù)(總穎花數(shù)的2倍)所得的百分?jǐn)?shù)(圖2-C, D)。P1、P2和F1分別調(diào)查10個(gè)單株的主莖穗, 10個(gè)單株主莖穗的均值即為該株系柱頭外露率表型值; F2和F2:3分別調(diào)查了350個(gè)和320個(gè)單株, 各單株分別調(diào)查主莖穗和2個(gè)分蘗穗, 3個(gè)穗子柱頭外露率的均值為該單株柱頭外露率的表型值。

    1.3 數(shù)據(jù)分析

    用SPSS軟件(SPSS, Inc., Chicago, IL, USA)進(jìn)行柱頭性狀統(tǒng)計(jì)分析及性狀間的相關(guān)性分析。遺傳力的計(jì)算方法:2= [F2-1/3(P1+P2+F1)]/F2。

    F2和F2:3群體4個(gè)柱頭相關(guān)性狀的次數(shù)分布按蓋鈞鎰[42]的方法進(jìn)行整理, 采用曹錫文等[41]研制的植物數(shù)量性狀分離分析軟件SEA-G4F2 (P1、P2、F1和F2)和SEA-G4F3 (P1、P2、F1和F2:3)中的極大似然法和IECM算法, 估算2個(gè)世代4個(gè)柱頭相關(guān)性狀獲得1對(duì)主基因(one major gene model, 1 MG)、2對(duì)主基因(two major genes model, 2 MG)、多基因(polygene model, PG)、1對(duì)主基因+多基因(mixed one major gene and polygene model, MX1)和2對(duì)主基因+多基因(mixed two major genes and polygene model, MX2)共5類24種遺傳模型的極大對(duì)數(shù)似然函數(shù)值和AIC (Akaike’s information criterion)值。根據(jù)最佳模型選擇原則, AIC值最小的模型為相對(duì)最適模型, 若出現(xiàn)模型間AIC值差異不大時(shí), 則會(huì)出現(xiàn)幾個(gè)備選模型, 此時(shí)需通過樣本分布與理論分布間的均勻性檢驗(yàn)和Kolmogorov-Smirnov檢驗(yàn)來比較各備選模型適合性概率的大小, 從而確定最適模型, 同時(shí)可得到最適模型相應(yīng)的各成分分布參數(shù), 包括各主基因的加性、顯性、上位性效應(yīng), 多基因的加性、顯性、上位性效應(yīng)以及主基因遺傳率和多基因遺傳率等。其中主基因遺傳率2mg(%) =2mg/2p′100; 多基因遺傳率2pg(%) =2pg/2p′100; 式中2mg、2pg和2p分別為主基因方差、多基因方差和表型方差。

    2 結(jié)果與分析

    2.1 水稻柱頭相關(guān)性狀的遺傳分析

    2.1.1 柱頭相關(guān)性狀表型特征及次數(shù)分布 E1和E2中紫泰S的STL、SYL、TSSL和PES均極顯著高于7001S (表1和圖3)。7001S / 紫泰S F1的4個(gè)柱頭相關(guān)性狀均介于雙親之間, 其中STL、SYL和TSSL的接近于雙親均值, 而PES的偏向于高值親本。4個(gè)性狀在分離世代群體(F2和F2:3)中的變異系數(shù)介于14.29%~54.94%, 說明該組合中的目標(biāo)性狀有較好的遺傳多樣性。2個(gè)環(huán)境中4個(gè)性狀的遺傳力均超過78%, 說明這4個(gè)性狀以遺傳為主。但是從表2的P1、P2和F1聯(lián)合方差分析的結(jié)果來看, 除TSSL外, 環(huán)境效應(yīng)以及基因型與環(huán)境之間的互作效應(yīng)對(duì)柱頭性狀的表型值也存在一定影響。分離群體中, 除F2中的STL外, 其余性狀在2個(gè)分離群體中均出現(xiàn)雙向超親個(gè)體, 表現(xiàn)出連續(xù)正態(tài)分布, 適于進(jìn)行遺傳分析(表1)。

    表1 7001S、紫泰S組合F2和F2:3群體柱頭相關(guān)性狀描述

    E1: 2018年廬江; E2: 2019年廬江。

    E1: Lujiang in 2018; E2: Lujiang in 2019. STL: stigma length; SYL: style length; TSSL: the sum of stigma and style length; PES: percentage of exserted stigma; SD: standard deviation; CV: coefficient of variation; H: heritability.

    表2 7001S、紫泰S和7001S/紫泰S F1 4個(gè)柱頭性狀的聯(lián)合方差分析

    **表示0.01水平顯著。

    STL: stigma length; SYL: style length; TSSL: the sum of stigma and style length; PES: percentage of exserted stigma; DF: degrees of freedom; SS: sum of squares; MS: mean square.**: significant differences at< 0.01.

    縮寫同表1。Abbreviations are the same as those given in Table 1.

    2.1.2 柱頭相關(guān)性狀遺傳模型的建立 根據(jù)5類24種模型以及柱頭性狀的次數(shù)分布, 算出該性狀所有24種遺傳模型的極大似然值和AIC值(表略)。AIC值最小的模型為最佳可能模型, 當(dāng)出現(xiàn)與最小AIC值接近的模型時(shí), 將這些模型與最佳可能模型一起都作為該性狀的備選模型。再對(duì)備選模型的樣本分布和理論分布進(jìn)行適合性檢驗(yàn)(12,22,32,2,D)以確定最佳模型。結(jié)果E1中STL的MX2-ADI-AD模型AIC值最小, 為-127.50, 但MX1-AD-ADI、MX2-ADI-ADI和MX2-AD-AD模型的AIC值與之相近, 故上述4個(gè)模型作為E1中STL的備選模型。應(yīng)用與E1中STL相同的方法進(jìn)行其他柱頭性狀備選模型的選擇。2個(gè)環(huán)境下4個(gè)柱頭性狀備選模型的極大似然值和AIC值見表3。

    表3 4個(gè)柱頭性狀備選模型配合表型分布的極大似然函數(shù)值和AIC值(IECM算法)

    a“適合性檢驗(yàn)”中5個(gè)數(shù)字分別指在1222322和D適合性檢驗(yàn)中達(dá)到顯著水平的統(tǒng)計(jì)量個(gè)數(shù)。MX: 主基因+多基因模型; MG: 主基因模型; PG: 多基因模型; A: 加性效應(yīng); D: 顯性效應(yīng); E: 等效; I: 上位性。E1: 2018年廬江; E2: 2019年廬江。粗體表示最適模型。

    aFive digitals in test of goodness-of-fit column denote numbers of significant statistic parameters among12,22,32,2, andD. MX: mixed major gene and polygene model; MG: major gene model; PG: polygene model; A: additive effect; D: dominance effect; E: equal; I: interaction. E1: 2018 Lujiang; E2: 2019 Lujiang. The model in bold in this table is the optimal model. STL: stigma length; SYL: style length; TSSL: the sum of stigma and style length; PES: percentage of exserted stigma.

    經(jīng)適合性檢驗(yàn), E1中STL的MX2-AD-AD模型20個(gè)適合性檢驗(yàn)統(tǒng)計(jì)量有4個(gè)達(dá)到顯著(α=0.05), 而MX1-AD-ADI、MX2-ADI-ADI和MX2-ADI-AD均只有2個(gè)統(tǒng)計(jì)量達(dá)到顯著, 但MX2-ADI-AD模型理論分布配合實(shí)際分布適合性的概率值比MX1-AD-ADI和MX2-ADI-ADI模型高, 說明MX2-ADI-AD模型較其他3個(gè)備選模型的適合性好。因此, 7001S/紫泰S組合E1中STL的最佳遺傳模型為MX2-ADI-AD。說明該性狀E1中表現(xiàn)出受2對(duì)主效基因和微效基因共同控制, 其中2對(duì)主效基因表現(xiàn)出加性-顯性-上位性效應(yīng), 而微效基因表現(xiàn)出加性-顯性效應(yīng)。E2中STL的PG-ADI和MX2-ADI-ADI模型20個(gè)適合性檢驗(yàn)統(tǒng)計(jì)量均有2個(gè)達(dá)到顯著(α=0.05); 而MX1-AD-ADI和MX2-ADI-AD模型均只有1個(gè)統(tǒng)計(jì)量達(dá)到顯著(α=0.05), 但MX2-ADI-AD模型理論分布配合實(shí)際分布適合性的概率值比MX1-AD-ADI模型高, 說明MX2-ADI-AD模型比MX1-AD-ADI模型適合性好。因此, 7001S/紫泰S組合E2中STL的最佳遺傳模型也為MX2-ADI-AD。說明該性狀E2中也表現(xiàn)出受2對(duì)加性-顯性-上位性主基因+加性-顯性多基因控制。

    應(yīng)用與STL相同方法, 對(duì)其他3個(gè)性狀進(jìn)行最佳遺傳模型的選擇。對(duì)于SYL, 2個(gè)環(huán)境中的最佳遺傳模型均為MX2-ADI-AD, 說明SYL在E1和E2中均表現(xiàn)出受2對(duì)加性-顯性-上位性主基因+加性-顯性多基因控制。對(duì)于TSSL, E1中最佳遺傳模型為MX2-EEAD-AD, 即TSSL在E1中表現(xiàn)出受2對(duì)主效基因和微效基因共同控制, 其中2對(duì)主效基因的加性效應(yīng)值和顯性效應(yīng)值分別相等, 而微效基因表現(xiàn)出加性-顯性效應(yīng); E2中TSSL最佳遺傳模型為MX2-ADI-AD, 說明TSSL在E2中表現(xiàn)為2對(duì)加性-顯性-上位性主基因+加性-顯性多基因控制。對(duì)于PES, E1中最佳遺傳模型為MX2-ADI-ADI, 即PES在E1中表現(xiàn)出受2對(duì)主效基因和微效基因共同控制,2對(duì)主效基因和微效基因均表現(xiàn)出加性-顯性-上位性效應(yīng); E2中PES最佳遺傳模型為MX1-AD-ADI模型, 即表現(xiàn)出受1對(duì)主效基因和微效基因共同控制, 主效基因表現(xiàn)出加性-顯性效應(yīng), 而微效基因則表現(xiàn)出加性-顯性-上位性效應(yīng)。

    E1中TSSL和E2中PES分別有5個(gè)和3個(gè)成分分布, 其余性狀在2個(gè)環(huán)境中均表現(xiàn)出9個(gè)成分分布。4個(gè)柱頭性狀在各自最佳遺傳模型下理論分布擬合曲線與實(shí)際分布曲線如圖4所示。

    2.1.3 柱頭性狀的遺傳參數(shù) 表4列出了4個(gè)柱頭性狀在2個(gè)環(huán)境中各自最佳模型下的遺傳參數(shù)估計(jì)值。從一階參數(shù)看, E1中STL的||+|j|+|j|+||>|a|+|b|>|h|+|h|, 說明主效基因中的上位性效應(yīng)要大于加性效應(yīng)和顯性效應(yīng), 以上位性效應(yīng)為主; E2中STL的||+|j|+|j|+||>|h|+|h|>|a|+|b|, 說明主效基因中的上位性效應(yīng)要大于顯性效應(yīng)和加性效應(yīng), 以上位性效應(yīng)為主; 2個(gè)環(huán)境中的||+|j|+|j|+||+ |h|+|h|+|a|+|b|與|[]|+|[]|值基本一致, 說明STL中微效基因也起重要作用。對(duì)于SYL, E1和E2中||+|j|+|j|+||>|h|+|h|>|a|+|b|, 說明主基因效應(yīng)中是以上位性效應(yīng)為主; 與STL一樣, 主效基因的各效應(yīng)值之和與微效基因的各效應(yīng)值之和基本一致,說明SYL中微效基因效應(yīng)也起重要作用。E1中TSSL的||>|[]|+|[]|, 說明主效基因的加性效應(yīng)大于微效基因的加性效應(yīng)和顯性效應(yīng)之和; E2中||+|j|+|j|+||>|h|+|h|>|a|+|b|, 說明2對(duì)主基因以上位性效應(yīng)為主; ||+|j|+|j|+||+|h|+|h|+|a|+ |b|=2.87, |[]|+|[]|=2.97, 說明TSSL中微效基因也起重要作用。對(duì)于PES, E1中|a|+|b|>||+|j|+|j|+||>|h|+|h|, E2中||>||, 說明PES是以加性效應(yīng)為主。

    從表4二階遺傳參數(shù)看, 2個(gè)環(huán)境中4個(gè)柱頭性狀的主基因遺傳率介于71.13%~97.09%, 多基因遺傳率介于0~22.21%, 說明這4個(gè)柱頭性狀均以主基因遺傳為主。

    2.2 水稻柱頭性狀間的相關(guān)性分析

    由表5可知, 4個(gè)柱頭性狀在7001S/紫泰S F2群體和F2:3群體中均表現(xiàn)出極顯著正相關(guān), 其中F2群體4個(gè)性狀間的相關(guān)系數(shù)介于0.262 (SYL和PES)和0.895 (STL和TSSL)之間, F2:3群體4個(gè)性狀間的相關(guān)系數(shù)介于0.281 (STL和PES)和0.885 (STL和TSSL)之間。柱花總長度和柱頭外露率在2個(gè)世代中的表型相關(guān)系數(shù)分別為0.377 (F2)和0.339 (F2:3), 呈極顯著正相關(guān)(2個(gè)性狀的散點(diǎn)圖見圖5), 決定系數(shù)(2)分別達(dá)到14.24%和11.51%, 說明在柱頭外露率的變異中, 有11%以上是由柱頭長度的變異引起的。這種極顯著的線性正相關(guān)具有明顯的生物學(xué)意義。

    3 討論

    本研究利用由2個(gè)核不育系雜交衍生的未經(jīng)人工選擇的F2(350個(gè)單株)和F2:3(320個(gè)單株)群體對(duì)4個(gè)柱頭性狀間的相關(guān)性進(jìn)行了分析。結(jié)果表明, 4個(gè)性狀兩兩之間均表現(xiàn)出極顯著正相關(guān), 這說明柱頭性狀是一個(gè)復(fù)雜性的綜合性狀。其中STL與SYL和TSSL間的相關(guān)系數(shù)均達(dá)到0.76以上, 呈緊密相關(guān), 這與Zhou等[16]和Dang等[17]利用不同水稻材料構(gòu)成的自然群體進(jìn)行分析所得的結(jié)果一致, 說明STL、SYL和TSSL可能是由相同的基因位點(diǎn)(QTL)控制的。TSSL與PES間的相關(guān)系數(shù)分別為0.377 (F2)和0.339 (F2:3), 達(dá)極顯著正相關(guān), 這與前人的研究結(jié)果一致[12-15]。同時(shí), 本研究還發(fā)現(xiàn)STL與PES間也表現(xiàn)出極顯著正相關(guān)。說明育種家可以通過增加水稻溫敏核不育系的柱花總長度(TSSL)來達(dá)到提高柱頭外露率(PES)的目的。同時(shí), 柱花總長度的增加, 必然會(huì)導(dǎo)致柱頭長度(STL)的增加。柱頭毛刷狀部分是柱頭接收花粉的器官, 柱頭長度的增加、柱頭外露率的提高必然導(dǎo)致雜交水稻異交結(jié)實(shí)率的增加, 從而增加雜交水稻的制種產(chǎn)量[12,43]。

    表4 4個(gè)柱頭性狀遺傳參數(shù)估計(jì)值

    a: 第1對(duì)主基因的加性效應(yīng);b: 第2對(duì)主基因的加性效應(yīng);a: 第1對(duì)主基因的顯性效應(yīng);b: 第2對(duì)主基因的顯性效應(yīng);: 加性×加性互作效應(yīng);ab: 加性×顯性互作效應(yīng);ba: 顯性×加性互作效應(yīng);: 顯性×顯性互作效應(yīng); []: 多基因加性效應(yīng); []: 多基因顯性效應(yīng);2p: 表型方差;2mg: 主基因方差;2pg: 多基因方差;2mg(%): 主基因遺傳率;2pg(%): 多基因遺傳率。

    a: additive effect of the first major gene;b: additive effect of the second major gene;a: dominant effect of the first major gene;b: dominant effect of the second major gene;: epistatic effect value betweenaandb;ab:epistatic effect value betweenaandb;ba: epistatic effect value betweenaandb;: epistatic effect value betweenhandh; []: additive effect of polygene; []: dominant effect of polygene;2p: phenotypic variance;2mg: major gene variance;2pg: polygene variance;2mg(%): heritability of major gene;2pg(%): heritability of polygene. STL: stigma length; SYL: style length; TSSL: the sum of stigma and style length; PES: percentage of exserted stigma.

    表5 4個(gè)柱頭性狀間的相關(guān)性系數(shù)

    左下角和右上角分別表示F2和F2:3的相關(guān)性系數(shù);**表示0.01水平顯著??s寫同表1。

    The lower-left corner and the upper-right corner indicate correlation coefficient in F2and F2:3, respectively;**indicates significant differences at the 0.01 probability level. Abbreviations are the same as those given in Table 1.

    解松峰等[40]認(rèn)為與傳統(tǒng)的遺傳分析方法相比, 主基因+多基因混合遺傳模型進(jìn)行分析所得的結(jié)果更加客觀, 育種家可通過表型值快速判斷目標(biāo)性狀的遺傳組成。近些年主基因+多基因遺傳模型更是被頻繁用于各種農(nóng)作物[36-40,44]。本研究利用植物數(shù)量性狀分離分析軟件中的SEA-G4F2 (P1、P2、F1和F2)和SEA-G4F3 (P1、P2、F1和F2:3)對(duì)4個(gè)柱頭性狀進(jìn)行了主基因+多基因遺傳分析。結(jié)果發(fā)現(xiàn), STL、SYL和TSSL均表現(xiàn)出受2對(duì)主效基因+微效基因控制, 且除E1中TSSL的2對(duì)主基因表現(xiàn)出等效的加性效應(yīng)和等效的顯性效應(yīng)外, 其余的2對(duì)主效基因間均表現(xiàn)出加性-顯性-上位性效應(yīng), 這也進(jìn)一步從遺傳上證實(shí)了STL、SYL和TSSL之間的顯著相關(guān)性, 說明在該群體中這3個(gè)性狀極可能是由相同的主效位點(diǎn)(QTL)控制的。而2個(gè)世代PES的最適模型并不一致, E1中PES表現(xiàn)出受2對(duì)主基因控制, 2對(duì)主基因間表現(xiàn)為加性-顯性-上位性效應(yīng), 而E2中則表現(xiàn)出受1對(duì)加性-顯性主效基因控制, 說明較之于柱頭長度相關(guān)性狀, PES更易受到環(huán)境因素的影響[18]。結(jié)合P1、P2和F1聯(lián)合方差分析結(jié)果發(fā)現(xiàn), 4個(gè)柱頭性狀均是以基因型遺傳為主, 但同時(shí)也受環(huán)境效應(yīng)的影響(表2)。除TSSL外, 環(huán)境效應(yīng)對(duì)STL、SYL和PES的影響均達(dá)到了極顯著水平, 且STL和SYL同時(shí)受到基因型′環(huán)境互作效應(yīng)的極顯著影響。從而導(dǎo)致2個(gè)環(huán)境TSSL和PES檢測(cè)到的主效基因數(shù)量以及主效基因間作用方式的不同, 需要通過其他方法, 如多環(huán)境下的QTL定位等, 進(jìn)行進(jìn)一步的驗(yàn)證。

    當(dāng)前, 對(duì)于柱頭性狀的研究主要集中于PES, 而對(duì)STL、SYL和TSSL的研究較少[15-17,22,28,31-34]。因此, 對(duì)于4個(gè)柱頭性狀遺傳效應(yīng)的研究也主要集中于PES。從一階遺傳參數(shù)來看, PES以加性效應(yīng)為主, 這與Yan等[29]和Li等[45]研究結(jié)果是一致的。而STL、SYL和TSSL則以上位性效應(yīng)遺傳為主, 說明這3個(gè)性狀的遺傳較PES復(fù)雜。同時(shí), 本研究發(fā)現(xiàn)的TSSL以上位性效應(yīng)為主的結(jié)果與劉強(qiáng)明[15]利用CSSL群體所得出的以加性效應(yīng)遺傳為主的結(jié)果是不一致的, 這可能是由于材料類型的不同造成的。從二階遺傳參數(shù)來看, 2個(gè)環(huán)境中4個(gè)性狀主基因遺傳率介于71.13% (E2中PES)和97.09% (E2中TSSL)之間, 說明這4個(gè)性狀均以主效基因遺傳為主, 而除F2:3群體中的PES受1對(duì)主基因+多基因控制外, 其余均表現(xiàn)出受2對(duì)主基因+多基因控制; 且除F2群體中TSSL的2對(duì)主基因表現(xiàn)為等效的加性效應(yīng)和等效的顯性效應(yīng)外, 其余2對(duì)主效基因間均表現(xiàn)出加性-顯性-上位性效應(yīng)。綜上結(jié)果一方面, 說明4個(gè)柱頭性狀為數(shù)量性狀, 研究材料以及環(huán)境等因素可能會(huì)對(duì)檢測(cè)到的主基因數(shù)量以及作用方式產(chǎn)生影響[28-30]; 但另一方面, 本研究同時(shí)利用F2和F2:3群體在2個(gè)環(huán)境下對(duì)4個(gè)柱頭性狀進(jìn)行遺傳分析, 其結(jié)果可以互相驗(yàn)證。從本研究結(jié)果分析, 7001S/紫泰S組合中的STL、SYL和TSSL是由2個(gè)主效位點(diǎn)控制的, 且控制這3個(gè)性狀的主效位點(diǎn)可能為相同位點(diǎn), 且這2個(gè)位點(diǎn)中至少有1個(gè)位點(diǎn)同時(shí)也是控制PES的主效位點(diǎn)。王春娥等[44]曾同時(shí)利用主基因+多基因模型和QTL定位對(duì)大豆豆腐和豆乳得率進(jìn)行了遺傳分析和QTL定位, 認(rèn)為這2種方法檢測(cè)到的主基因數(shù)量具有相對(duì)一致性。說明主基因+多基因混合遺傳模型與QTL定位間存在緊密聯(lián)系, 兩者相輔相成。實(shí)踐中可以利用遺傳分析的結(jié)果對(duì)QTL定位的結(jié)果進(jìn)行校對(duì), 便于育種家根據(jù)表型快速得到目標(biāo)性狀的遺傳組成; 同時(shí), 也可以用QTL定位的結(jié)果來校正遺傳分析的結(jié)果, 這樣可進(jìn)一步更加準(zhǔn)確了解QTL定位所反映目標(biāo)性狀的遺傳組成。本研究利用主基因+多基因模型對(duì)柱頭4個(gè)相關(guān)性狀進(jìn)行遺傳分析, 其中柱頭外露易于觀察, 水稻開花后田間可直接觀察, 但柱頭性狀較小, 田間肉眼很難直接進(jìn)行比較, 需取樣室內(nèi)測(cè)量。育種家在高異交率溫敏核不育選育的過程中, 可將田間柱頭外露率的觀察和室內(nèi)柱頭長度性狀的調(diào)查相結(jié)合, 通過長柱頭(毛刷長度和柱花總長度)、高外露率的選擇進(jìn)而改良水稻的異交率。同時(shí), 育種家在進(jìn)行長柱頭、高外露率溫敏核不育系選育過程中, 要充分考慮主效基因和微效基因的影響, 及時(shí)開展早期有效的單株選擇, 重點(diǎn)關(guān)注高代單株性狀的后期穩(wěn)定表現(xiàn), 排除或降低環(huán)境因素對(duì)目標(biāo)性狀選擇造成的影響。主基因+多基因遺傳分析對(duì)4個(gè)柱頭相關(guān)性狀的遺傳組成、遺傳規(guī)律特征及遺傳參數(shù)進(jìn)行了解析, 這將會(huì)為水稻柱頭性狀的改良提供有意義的遺傳信息, 為下一步通過QTL初定位及精細(xì)定位、闡明性狀間的遺傳調(diào)控網(wǎng)絡(luò)和分子設(shè)計(jì)育種提供重要理論依據(jù)。

    4 結(jié)論

    2個(gè)環(huán)境中4個(gè)柱頭性狀兩兩間均表現(xiàn)出極顯著正相關(guān), 相關(guān)系數(shù)介于0.262 (E1中PES和SYL)和0.895 (E2中STL和TSSL)之間。STL、SYL和TSSL均表現(xiàn)出受2對(duì)主基因+多基因控制, 除E1中TSSL表現(xiàn)為等加性-等顯性主基因+加性-顯性多基因控制外, 其余均表現(xiàn)為加性-顯性-上位性主基因+加性-顯性多基因。E1中PES表現(xiàn)出受2對(duì)加性-顯性-上位性主基因+加性-顯性-上位性多基因控制, 而E2中則表現(xiàn)出受1對(duì)加性-顯性主基因+加性-顯性-上位性多基因。2個(gè)環(huán)境中4個(gè)性狀均以主基因遺傳為主。本研究結(jié)果為創(chuàng)制長柱頭、高柱頭外露率核不育系材料提供了理論支持, 同時(shí)也為柱頭長度相關(guān)性狀主效位點(diǎn)的挖掘以及分子標(biāo)記輔助選擇提供理論依據(jù)。

    [1] 袁隆平. 發(fā)展雜交水稻保障糧食安全. 第一屆中國雜交水稻大會(huì)論文集. 湖南長沙, 2010. pp 1–2. Yuan L P. Developing hybrid rice to ensure food security. In: Collection of the 1st Chinese Hybrid Rice Conference. Changsha, Hunan, China, 2010. pp 1–2 (in Chinese).

    [2] Huang X H, Han B. Rice domestication occurred through single origin and multiple introgressions.,2016, 2: 16207.

    [3] Qian Q, Guo L B, Smith S M, Li J Y. Breeding high-yield superior quality hybrid super rice by rational design., 2016, 3: 283–294.

    [4] Xia Y M, Tang N, Hu Y Y, Li D, Li S C, Bu X L, Yu M L, Qi S W, Yang Y S, Zhu H J, Cao C Y, Li P, Yuan L P, Cao M L. A method for mechanized hybrid rice seed production using female sterile rice., 2019, 12: 39.

    [5] 劉文炳. 三系雜交水稻超高產(chǎn)制種配套組裝技術(shù). 第一屆中國雜交水稻大會(huì)論文集. 湖南長沙, 2010. pp 465–468. Liu W B. Packaging techniques for super-high-yield seed production of hybrid rice. In: Collection of the 1st Chinese Hybrid Rice Conference. Changsha, Hunan, China, 2010. pp 465–468 (in Chinese).

    [6] 舒志芬, 陳勇, 劉釗, 袁露, 熊朝, 張海清. 5個(gè)水稻光溫敏核不育系柱頭特性與異交結(jié)實(shí)率的關(guān)系. 作物研究, 2015, 29: 343–347. Shu Z F, Chen Y, Liu Z, Yuan L, Xiong C, Zhang H Q. Correlation between outcrossing characteristics and stigma feature of 5 PTGMS lines., 2015, 29: 343–347 (in Chinese with English abstract).

    [7] 田大成, 黃三奎, 段永國, 王友紅. 水稻不育系花時(shí)和受粉時(shí)間與異交結(jié)實(shí)率的關(guān)系. 雜交水稻, 2004, 19(3): 50–54. Tian D C, Huang S K, Duan Y G, Wang Y H. The relationship between flowering and pollination time and outcrossing rate of male sterile lines in hybrid rice seed production., 2004, 19(3): 50–54 (in Chinese with English abstract).

    [8] 代貴金, 華澤田, 蔡偉, 王妍, 張忠旭, 王大為. 激素在雜交粳稻制種中的應(yīng)用研究. 雜交水稻, 1999, 14(增刊): 15–20. Dai G J, Hua Z T, Cai W, Wang Y, Zhang Z X, Wang D W. Application of plant growth substances in the seed production ofhybrid rice., 1999, 14(S): 15–20 (in Chinese).

    [9] 田大成, 張素英, 秦春林. 提高柱頭外露率是雜交稻制種穩(wěn)產(chǎn)高產(chǎn)的關(guān)鍵措施. 四川農(nóng)業(yè)科技, 1990, (2): 16–18. Tian D C, Zhang S Y, Qin C L. Increasing stigma exsertion rate is a key measure for hybrid rice seed production., 1990, (2): 16–18 (in Chinese).

    [10] 李慶榮. 柱頭外露率的環(huán)境變異及其在異交結(jié)實(shí)中的作用. 種子科技, 1995, (2): 32–33. Li Q R. Environmental variation of stigma exsertion rate and its role in outcrossing., 1995, (2): 32–33 (in Chinese).

    [11] 楊保漢. 不育系柱頭外露率及其結(jié)實(shí)率研究. 雜交水稻, 1997, 12(1): 13–15. Yang B H. Studies on stigma exsertion rate and outcrossing rate of CMS lines in rice., 1997, 12(1): 13–15 (in Chinese with English abstract).

    [12] Miyata M, Yamamoto T, Komori T, Nitta N. Marker-assisted selection and evaluation of the QTL for stigma exsertion underrice genetic background., 2007, 114: 539–548.

    [13] Virmani S S, Athwal D S. Genetic variability in floral characteristics influencing outcrossing in Oryza sativa L. Crop Sci, 1973, 13: 66–67.

    [14] Takano-Kai N, Doi K, Yoshimura A. GS3participates in stigma exsertion as well as seed length in rice. Breed Sci, 2011, 61: 244–250.

    [15] 劉強(qiáng)明. 水稻異交相關(guān)性狀的QTL遺傳剖析與柱頭長度QTL的精細(xì)定位. 南京農(nóng)業(yè)大學(xué)博士學(xué)位論文, 江蘇南京, 2015. Liu Q M. QTL Dissection of Out-crossing Related Traits and Fine Mapping of a Stigma Length QTLin Rice (L.). PhD Dissertation of Nanjing Agricultural University, Nanjing, Jiangsu, China, 2015 (in Chinese with English abstract).

    [16] Zhou H, Li P, Xie W, Hussain S, Li Y, Xia D, Zhao H, Sun S, Chen J, Ye H, Hou J, Zhao D, Gao G, Zhang Q, Wang G, Lian X, Xiao J, Yu S, Li X, He Y. Genome-wide association analyses reveal the genetic basis of stigma exsertion in rice., 2017, 10: 634–644.

    [17] Dang X, Liu E, Liang Y, Liu Q, Breria C M, Hong D. QTL detection and elite alleles mining for stigma traits inby association mapping., 2016, 7: 1188.

    [18] Liu Y, Zhang A, Wang F, Kong D, Li M, Bi J, Zhang F, Wang J, Luo X, Pan Z, Yu X, Liu G, Luo L. Fine mapping a quantitative trait locus,, that controls stigma exsertion rate in rice (L.)., 2019, 12: 46.

    [19] Xiong L X, Liu K D, Dai X K, Xu C G, Zhang Q. Identification of genetic factors controlling domestication-related traits of rice using an F2 population of a cross between Oryza sativa and O. rufipogon., 1999, 98: 243–251.

    [20] 李威, 圣忠華, 朱子亮, 魏祥進(jìn), 石磊, 鄔亞文, 唐紹清, 王建龍, 胡培松. 粳稻柱頭外露率QTL定位. 中國水稻科學(xué), 2017, 31: 23–30. Li W, Sheng Z H, Zhu Z L, Wei X J, Shi L, Wu Y W, Tang S Q, Wang J L, Hu P S. QTL mapping ofrice stigma exsertion rate., 2017, 31: 23–30 (in Chinese with English abstract).

    [21] Bakti C, Tanaka J. Detection of dominant QTLs for stigma exsertion ratio in rice derived from Oryza rufipogon accession ‘W0120’., 2019, 69: 143–150.

    [22] Uga Y, Fukuta Y, Cai H W, Iwata H, Ohsawa R, Morishima H, Fujimura T. Mapping QTLs influencing rice floral morphology using recombinant inbred lines derived from a cross betweenL. andGriff., 2003, 107: 218–226.

    [23] Yu X Q, Mei H W, Luo L J, Liu G L, Zou G H, Hu S P, Li M S, Wu J H. Dissection of additive, epistatic effect and Q × E interaction of quantitative trait loci influencing stigma exsertion under water stress in rice., 2006, 33: 542–550.

    [24] Hu S P, Zhou Y, Zhang L, Zhu X D, Wang Z G, Li L, Luo L J, Zhou Q M. QTL analysis of floral traits of rice (L.) under well-watered and drought stress conditions., 2009, 31: 173–181.

    [25] Li P B, Feng F C, Zhang Q L, Chao Y, Gao G J, He Y Q. Genetic mapping and validation of quantitative trait loci for stigma exsertion rate in rice., 2014, 34: 2131–2138.

    [26] 尹成, 李平波, 高冠軍, 張慶路, 羅利軍, 何予卿. 水稻柱頭外露率QTL定位. 分子植物育種, 2014, 12: 43–49. Yin C, Li P B, Gao G J, Zhang Q L, Luo L J, He Y Q. QTL analysis of percentage of exserted stigma in rice., 2014, 12: 43–49 (in Chinese with English abstract).

    [27] Rahman M H, Yu P, Zhang Y X, Sun L P, Wu W X, Shen X H, Zhan X D, Chen D B, Cao L Y, Cheng S H. Quantitative trait loci mapping of the stigma exertion rate and spikelet number per panicle in rice (L.)., 2016, 15: gmr15048432.

    [28] Yan W G, Li Y, Agrama H A, Luo D G, Gao F Y, Lu X J, Ren G J. Association mapping of stigma and spikelet characteristics in rice (L.)., 2009, 24: 277–292.

    [29] Zhang K, Zhang Y, Wu W, Zhan X, BakrAnis G, Rahman M H, Hong Y, Riaz A, Zhu A, Cao Y, Sun L, Yang Z, Yang Q, Cao L, Cheng S. qSE7 is a major quantitative trait locus (QTL) influencing stigma exsertion rate in rice (L.)., 2018, 8: 14523.

    [30] Rahman M H, Zhang Y, Zhang K, Rahman M S, Barman H N, Riaz A, Chen Y, Wu W, Zhan X, Cao L, Cheng S. Genetic dissection of the major quantitative trait locus (), and its validation as the major influence on the rate of stigma exsertion in rice (L.)., 2017, 8: 1818.

    [31] Uga Y, Siangliw M, Nagamine T, Ohsawa R, Fujimura T, Fukuta Y. Comparative mapping of QTLs determining glume, pistil and stamen sizes in cultivated rice (Oryza sativa L.). Plant Breed, 2010, 129: 657–669.

    [32] Marathi B, Ramos J, Hechanova S L, Oane R H, Jena K K. SNP genotyping and characterization of pistil traits revealing a distinct phylogenetic relationship among the species of., 2015, 201: 131–148.

    [33] 陳蘭, 張紅, 張啟武, 江建華, 王洋, 陳杰麗, 洪德林. 水稻6個(gè)異交相關(guān)性狀的SSR關(guān)聯(lián)分析. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 35(2): 1–9. Chen L, Zhang H, Zhang Q W, Jiang J H, Wang Y, Chen J L, Hong D L. Association analysis of six outcrossing related traits in rice (L.) with SSR markers., 2012, 35(2): 1–9 (in Chinese with English abstract).

    [34] Liu Q, Qin J, Li T, Liu E, Fan D, Edzesi W M, Liu J, Jiang J, Liu X, Xiao L, Liu L, Hong D. Fine mapping and candidate gene analysis of qSTL3, a stigma length-conditioning locus in rice (L.). PLoS One, 2015, 10: e0127938.

    [35] 李成荃, 許克農(nóng), 王守海, 羅彥長, 王德正. 粳型水稻光敏核不育系7001S的育性與利用研究. 安徽農(nóng)業(yè)科學(xué), 1994, 22(1): 11–15. Li C Q, Xu K L, Wang S H, Luo Y C, Wang D Z. Study on the fertility ofPGMR 7001S and its utilization., 1994, 22(1): 11–15 (in Chinese with English abstract).

    [36] 張啟武, 江建華, 姚瑾, 洪德林. 穞稻與粳稻恢復(fù)系C堡籽粒灌漿速率的特征及遺傳分析. 作物學(xué)報(bào), 2009, 35: 1229–1235. Zhang Q W, Jiang J H, Yao J, Hong D L. Characterization and genetic analysis of grain filling rate of Ludao and restorer line C-bao inrice (L.)., 2009, 35: 1229–1235 (in Chinese with English abstract).

    [37] 汪文祥, 胡瓊, 梅德圣, 李云昌, 周日金, 王會(huì), 成洪濤, 付麗, 劉佳. 甘藍(lán)型油菜分枝角度主基因+多基因混合遺傳模型及遺傳效應(yīng). 作物學(xué)報(bào), 2016, 42: 1103–1111.Wang W X, Hu Q, Mei D S, Li Y C, Zhou R J, Wang H, Cheng H T, Fu L, Liu J. Genetic effects of branch angle using mixture model of major gene plus polygene inL., 2016, 42: 1103–1111 (in Chinese with English abstract).

    [38] 梁慧珍, 余永亮, 許蘭杰, 楊紅旗, 董薇, 譚政偉, 李磊, 裴新涌, 劉新梅. 大豆α-生育酚的遺傳與QTL分析. 中國農(nóng)業(yè)科學(xué), 2019, 52: 11–20. Liang H Z, Yu Y L, Xu L J, Yang H Q, Dong W, Tan Z W, Li L, Pei X Y, Liu X M. Inheritance and QTL mapping for α-tocopherol in soybean., 2019, 52: 11–20 (in Chinese with English abstract).

    [39] 翟英. 番茄果實(shí)色澤和色素含量的遺傳特征. 分子植物育種, 2019, 17: 264–269. Zhai Y. Genetic characteristics of color and pigment content in Tomato fruits., 2019, 17: 264–269 (in Chinese with English abstract).

    [40] 解松峰, 吉萬全, 張耀元, 張俊杰, 胡衛(wèi)國, 李俊, 王長有, 張宏, 陳春環(huán). 小麥重要產(chǎn)量性狀的主基因+多基因混合遺傳分析. 作物學(xué)報(bào), 2020, 46: 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 effects of important yield traits analysed by mixture model of major gene plus polygene in wheat., 2020, 46: 365–384 (in Chinese with English abstract).

    [41] 曹錫文, 劉兵, 章元明. 植物數(shù)量性狀分離分析Windows軟件包SEA的研制. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2013, 36(6): 1–6. Cao X W, Liu B, Zhang Y M. SEA: a software package of segregation analysis of quantitative traits in plants., 2013, 36(6): 1–6 (in Chinese with English abstract).

    [42] 蓋鈞鎰. 試驗(yàn)統(tǒng)計(jì)方法. 北京: 中國農(nóng)業(yè)出版社, 2000. pp 99–126, 157–189. Gai J Y. Methods of Experiments and Statistics. Beijing: China Agriculture Press, 2000. pp 99–126, 157–189 (in Chinese).

    [43] 吳爽, 江建華, 湯修竹, 王德正. 水稻溫敏核不育系廣茉S的異交特性研究. 雜交水稻, 2015, 30(3): 35–38. Wu S, Jiang J H, Tang X J, Wang D Z. Studies on outcrossing characteristics of TGMS line Guangmo S in rice., 2015, 30(3): 35–38 (in Chinese with English abstract).

    [44] 王春娥, 蓋鈞鎰, 傅三雄, 喻德躍, 陳受宜. 大豆豆腐和豆乳得率的遺傳分析與QTL定位. 中國農(nóng)業(yè)科學(xué), 2008, 41: 1274–1282. Wang C E, Gai J Y, Fu S X, Yu D Y, Chen S Y. Inheritance and QTL mapping of Tofu and soymilk output in soybean., 2008, 41: 1274–1282 (in Chinese with English abstract).

    [45] Li T, Chen Y W. Genetics of stigma exsertion in rice., 1985, 2: 84–85.

    Genetic analysis of stigma traits with genic male sterile line by mixture model of major gene plus polygene in rice (L.)

    JIANG Jian-Hua1,**, ZHANG Wu-Han2,**, DANG Xiao-Jing1, RONG Hui3, YE Qin3, HU Chang-Min1, ZHANG Ying1, HE Qiang2,*, and WANG De-Zheng1,*

    1Institute of Rice Research, Anhui Academy of Agricultural Sciences / Rice Genetic Breeding Key Laboratory of Anhui Province, Hefei 230031, Anhui, China;2State Key Laboratory of Hybrid Rice / Hunan Hybrid Rice Research Center, Changsha 410125, Hunan, China;3Anhui Science and Technology University, Fengyang 233100, Anhui, China

    Rice stigma traits are important traits that affects the yield of male sterile lines and hybrid rice seed production. To provide genetic information for the creation of thermo-sensitive genic male sterile lines with long stigmas and high exsertion rates in rice, we investigated the correlation among four stigma traits of 7001S (photoperiod-thermo-sensitive genic male sterile line with short stigma length and low exsertion rate, P1), Zitai S (thermo-sensitive genic male sterile line with long stigma length and high exsertion rate, P2), and F1, F2(350 lines), and F2:3(320 lines) from the cross P1× P2. Coefficients of linear correlation in stigma traits were calculated, and genetic patterns were analyzed by mixed major-gene plus polygene inheritance models. The results showed that there were extremely significant linear positive correlations among four stigma traits. The correlation coefficients were between 0.262 and 0.895. Genetic analysis revealed that stigma length (STL), style length (SYL), the sum of stigma, and style length (TSSL) were controlled by two major genes plus polygenes. The two major genes expressed additive-dominant-epistatic effects, but they expressed equal additive effects and equal dominant effects of TSSL in F2population. The epistatic effect was dominant in STL, SYL, and TSSL in F2and F2:3. Percentage of exserted stigma (PES) was controlled by two major genes plus polygenes with additive-dominant-epistatic effect in F2. However, PES was controlled by one major gene with additive-dominant effect plus polygenes in F2:3. The additive effect was dominant of PES in F2and F2:3. Four stigma traits were mainly governed by major genes.

    rice (L.); genic male sterile line; stigma trait; major gene plus polygene model; genetic analysis

    10.3724/SP.J.1006.2021.02057

    本研究由雜交水稻國家重點(diǎn)實(shí)驗(yàn)室(湖南雜交水稻研究中心)開放課題基金(2018KF01), 水稻遺傳育種安徽省重點(diǎn)實(shí)驗(yàn)室開放基金(SDKF-201904)和國家自然科學(xué)基金項(xiàng)目(31601374)資助。

    This study was supported by the Open Research Fund of National Key Laboratory of Hybrid Rice (Hunan Hybrid Rice Research Center) (2018KF01), the Program of Rice Genetic Breeding Key Laboratory of Anhui Province (SDKF-201904), and the National Natural Science Foundation of China (31601374).

    何強(qiáng), E-mail: heqiang@ hhrrc.ac.cn, Tel: 0731-82872990; 王德正, E-mail: wdzhengqin@163.com, Tel: 0551-62160454

    **同等貢獻(xiàn)(Contributed equally to this work)

    江建華, E-mail: peanutlabjjh@163.com, Tel: 0551-62160454

    2020-08-18;

    2020-12-01;

    2020-12-31.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20201231.1007.004.html

    猜你喜歡
    加性柱頭顯性
    BS型小麥光溫敏雄性不育系柱頭外露規(guī)律研究
    ?2?4[u]-加性循環(huán)碼
    川鄂鹽道白菜柱頭圖形符號(hào)意義分析
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關(guān)系:超可加性、次可加性或不可加性
    企業(yè)家多重政治聯(lián)系與企業(yè)績效關(guān)系:超可加性、次可加性或不可加性
    顯性激勵(lì)與隱性激勵(lì)對(duì)管理績效的影響
    卷柱頭薹草草坪建植技術(shù)探析
    秈型光溫敏核不育系柱頭性狀表現(xiàn)
    中國稻米(2017年6期)2017-12-28 02:13:01
    社會(huì)權(quán)顯性入憲之思考
    基于加性指標(biāo)的網(wǎng)絡(luò)斷層掃描的研究
    国产成人精品在线电影| 无人区码免费观看不卡| 丝袜美足系列| 99精品久久久久人妻精品| 丝瓜视频免费看黄片| 午夜成年电影在线免费观看| 精品福利永久在线观看| 欧美成人午夜精品| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 美国免费a级毛片| 成人18禁高潮啪啪吃奶动态图| 一级,二级,三级黄色视频| 国产欧美日韩一区二区三区在线| 丁香六月欧美| 欧美精品人与动牲交sv欧美| 99热只有精品国产| 老熟妇乱子伦视频在线观看| 我的亚洲天堂| 午夜精品在线福利| 1024香蕉在线观看| 欧美乱色亚洲激情| 桃红色精品国产亚洲av| 涩涩av久久男人的天堂| 黄片播放在线免费| 久久九九热精品免费| 国产高清激情床上av| 成年人午夜在线观看视频| 两个人看的免费小视频| 天天添夜夜摸| 欧美+亚洲+日韩+国产| 看片在线看免费视频| 极品人妻少妇av视频| 五月开心婷婷网| 日韩免费av在线播放| 50天的宝宝边吃奶边哭怎么回事| 99香蕉大伊视频| 捣出白浆h1v1| 欧美乱色亚洲激情| 男女下面插进去视频免费观看| 午夜视频精品福利| av不卡在线播放| 日韩大码丰满熟妇| 亚洲国产毛片av蜜桃av| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区 | 国产免费现黄频在线看| 很黄的视频免费| 欧美精品高潮呻吟av久久| 亚洲午夜理论影院| 免费av中文字幕在线| 亚洲一区二区三区欧美精品| av在线播放免费不卡| 久久久久精品国产欧美久久久| 亚洲aⅴ乱码一区二区在线播放 | 一级毛片高清免费大全| 18在线观看网站| 一a级毛片在线观看| 国产亚洲精品一区二区www | 久久精品亚洲精品国产色婷小说| x7x7x7水蜜桃| 国产又色又爽无遮挡免费看| 女警被强在线播放| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 99国产精品免费福利视频| 操美女的视频在线观看| 丝袜在线中文字幕| 两性夫妻黄色片| 激情视频va一区二区三区| 国产国语露脸激情在线看| 一本一本久久a久久精品综合妖精| 麻豆乱淫一区二区| 多毛熟女@视频| 亚洲精华国产精华精| av欧美777| 日本黄色视频三级网站网址 | 国产精品电影一区二区三区 | 国产精品永久免费网站| 久久久久久久久久久久大奶| 国产欧美日韩一区二区三| 岛国在线观看网站| 国产欧美日韩综合在线一区二区| 亚洲第一av免费看| 成人永久免费在线观看视频| 很黄的视频免费| 国产亚洲精品久久久久5区| 精品久久久久久,| av不卡在线播放| 18禁黄网站禁片午夜丰满| 性色av乱码一区二区三区2| 亚洲国产看品久久| 国产无遮挡羞羞视频在线观看| 啦啦啦视频在线资源免费观看| 久久久久精品国产欧美久久久| a在线观看视频网站| 午夜两性在线视频| 热99re8久久精品国产| 亚洲黑人精品在线| 少妇被粗大的猛进出69影院| 超碰97精品在线观看| 免费在线观看视频国产中文字幕亚洲| 欧美乱妇无乱码| 亚洲av欧美aⅴ国产| 99国产精品一区二区蜜桃av | 欧美日韩乱码在线| 国产在线一区二区三区精| 国产高清视频在线播放一区| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩黄片免| 久久国产精品影院| 男男h啪啪无遮挡| 精品国产乱子伦一区二区三区| 日韩人妻精品一区2区三区| 亚洲性夜色夜夜综合| 热99久久久久精品小说推荐| 女同久久另类99精品国产91| 免费av中文字幕在线| 久久午夜亚洲精品久久| 69av精品久久久久久| 国产精品一区二区在线观看99| 国产日韩欧美亚洲二区| 亚洲专区国产一区二区| 两个人看的免费小视频| 国产男女内射视频| 天天操日日干夜夜撸| 宅男免费午夜| ponron亚洲| 9热在线视频观看99| 国产高清videossex| 欧美日韩成人在线一区二区| 国产亚洲精品久久久久5区| 天堂√8在线中文| 久久精品国产亚洲av香蕉五月 | www.自偷自拍.com| 欧美最黄视频在线播放免费 | 成人免费观看视频高清| 欧美日韩瑟瑟在线播放| 天天添夜夜摸| 在线免费观看的www视频| 久久国产精品男人的天堂亚洲| 一二三四在线观看免费中文在| 国产在线精品亚洲第一网站| 欧美日韩视频精品一区| 免费观看人在逋| 国产av一区二区精品久久| 国产亚洲一区二区精品| 一级毛片高清免费大全| www.精华液| 69精品国产乱码久久久| 国内毛片毛片毛片毛片毛片| 他把我摸到了高潮在线观看| 亚洲男人天堂网一区| 久久国产精品影院| 欧美大码av| 正在播放国产对白刺激| 午夜老司机福利片| 啪啪无遮挡十八禁网站| 91精品三级在线观看| 亚洲少妇的诱惑av| 女同久久另类99精品国产91| 男男h啪啪无遮挡| √禁漫天堂资源中文www| 亚洲一卡2卡3卡4卡5卡精品中文| 男女下面插进去视频免费观看| 国产蜜桃级精品一区二区三区 | 天天添夜夜摸| 午夜免费观看网址| 欧美 日韩 精品 国产| 久久精品国产综合久久久| 欧美国产精品一级二级三级| 69精品国产乱码久久久| 男女免费视频国产| 久久天堂一区二区三区四区| 激情在线观看视频在线高清 | 免费av中文字幕在线| 热99久久久久精品小说推荐| 亚洲精品国产区一区二| 国产在视频线精品| 欧美大码av| 两人在一起打扑克的视频| 国产精品偷伦视频观看了| 99国产综合亚洲精品| 精品国产美女av久久久久小说| 伊人久久大香线蕉亚洲五| 岛国在线观看网站| 最近最新中文字幕大全电影3 | 18禁黄网站禁片午夜丰满| 亚洲av熟女| 91精品三级在线观看| 国产淫语在线视频| 悠悠久久av| 久久久久国产精品人妻aⅴ院 | 免费日韩欧美在线观看| √禁漫天堂资源中文www| 很黄的视频免费| av超薄肉色丝袜交足视频| 99re6热这里在线精品视频| 黄色 视频免费看| 黑人巨大精品欧美一区二区mp4| 一本综合久久免费| 91九色精品人成在线观看| 免费在线观看日本一区| 一边摸一边抽搐一进一小说 | 亚洲五月天丁香| 欧美激情极品国产一区二区三区| avwww免费| 99国产精品一区二区三区| 丰满迷人的少妇在线观看| 亚洲一区高清亚洲精品| 午夜亚洲福利在线播放| 国产精品久久久久久精品古装| 成年女人毛片免费观看观看9 | 久久久精品国产亚洲av高清涩受| 国产不卡av网站在线观看| 国产精品国产高清国产av | 桃红色精品国产亚洲av| 天堂√8在线中文| 成人av一区二区三区在线看| 欧美乱码精品一区二区三区| 熟女少妇亚洲综合色aaa.| 久久久久久久精品吃奶| 99re在线观看精品视频| 亚洲熟女毛片儿| 精品少妇一区二区三区视频日本电影| 欧美 亚洲 国产 日韩一| 窝窝影院91人妻| 亚洲色图综合在线观看| 看片在线看免费视频| 两人在一起打扑克的视频| 久久久久国产一级毛片高清牌| 99热国产这里只有精品6| 亚洲一区二区三区欧美精品| 19禁男女啪啪无遮挡网站| 国产不卡一卡二| 久久久国产欧美日韩av| 极品少妇高潮喷水抽搐| 少妇被粗大的猛进出69影院| 亚洲片人在线观看| 午夜免费观看网址| 国产日韩欧美亚洲二区| 亚洲精品在线观看二区| 久久天堂一区二区三区四区| 自线自在国产av| 纯流量卡能插随身wifi吗| 久久香蕉精品热| 天堂俺去俺来也www色官网| 成年人午夜在线观看视频| 久久久久国产精品人妻aⅴ院 | 妹子高潮喷水视频| 亚洲第一欧美日韩一区二区三区| 一级片免费观看大全| 亚洲九九香蕉| 亚洲熟妇中文字幕五十中出 | 国产成人精品久久二区二区免费| 91老司机精品| 欧美中文综合在线视频| 视频区图区小说| 黄色视频,在线免费观看| 国产又色又爽无遮挡免费看| 水蜜桃什么品种好| 美女扒开内裤让男人捅视频| 夜夜爽天天搞| 免费在线观看完整版高清| 色老头精品视频在线观看| 午夜成年电影在线免费观看| 中国美女看黄片| 亚洲第一青青草原| 动漫黄色视频在线观看| 久久影院123| av超薄肉色丝袜交足视频| 久久精品国产清高在天天线| 婷婷成人精品国产| 亚洲色图综合在线观看| 99久久国产精品久久久| 人妻一区二区av| 丝袜在线中文字幕| 9191精品国产免费久久| 日日爽夜夜爽网站| 亚洲中文日韩欧美视频| 久久国产亚洲av麻豆专区| 搡老乐熟女国产| 99久久综合精品五月天人人| 人妻久久中文字幕网| 亚洲精品国产一区二区精华液| 成人三级做爰电影| 国产精品欧美亚洲77777| 国产男女内射视频| 亚洲avbb在线观看| 多毛熟女@视频| www.精华液| 国产高清videossex| 9191精品国产免费久久| 在线国产一区二区在线| 亚洲av成人一区二区三| av电影中文网址| 欧美精品啪啪一区二区三区| 最新美女视频免费是黄的| 亚洲成av片中文字幕在线观看| 深夜精品福利| 成人影院久久| 91大片在线观看| 亚洲情色 制服丝袜| 亚洲精品久久午夜乱码| 黄色毛片三级朝国网站| 欧美日韩瑟瑟在线播放| 在线观看www视频免费| 成人av一区二区三区在线看| 亚洲成人国产一区在线观看| 成年动漫av网址| 精品人妻在线不人妻| 国产av一区二区精品久久| 国产不卡av网站在线观看| 天堂动漫精品| 色94色欧美一区二区| 日本欧美视频一区| av欧美777| 国产精品免费一区二区三区在线 | 国产精品98久久久久久宅男小说| 12—13女人毛片做爰片一| 嫩草影视91久久| 精品人妻在线不人妻| ponron亚洲| 久热爱精品视频在线9| 国产乱人伦免费视频| 久热爱精品视频在线9| 午夜影院日韩av| 成人特级黄色片久久久久久久| 亚洲专区国产一区二区| 少妇粗大呻吟视频| 三上悠亚av全集在线观看| 亚洲精品国产区一区二| 亚洲全国av大片| av天堂在线播放| 日韩成人在线观看一区二区三区| 欧美av亚洲av综合av国产av| 亚洲国产看品久久| 一进一出抽搐gif免费好疼 | av超薄肉色丝袜交足视频| av福利片在线| 亚洲人成伊人成综合网2020| 国产精品美女特级片免费视频播放器 | 一级a爱视频在线免费观看| 亚洲专区国产一区二区| 国产在线一区二区三区精| 精品久久久久久,| 黄色女人牲交| 女性生殖器流出的白浆| 亚洲黑人精品在线| 欧美国产精品一级二级三级| 国产国语露脸激情在线看| 亚洲avbb在线观看| 国产成人免费无遮挡视频| av不卡在线播放| 91老司机精品| 91在线观看av| 又紧又爽又黄一区二区| 精品久久久久久,| 人妻 亚洲 视频| 免费在线观看完整版高清| 成在线人永久免费视频| 亚洲午夜理论影院| 久久影院123| 看免费av毛片| 久久ye,这里只有精品| 久久精品国产99精品国产亚洲性色 | 香蕉丝袜av| 成人免费观看视频高清| 亚洲精品国产精品久久久不卡| 日韩视频一区二区在线观看| 欧美日韩乱码在线| 久久久久国内视频| 美女 人体艺术 gogo| 精品欧美一区二区三区在线| 高潮久久久久久久久久久不卡| 亚洲国产欧美网| 老司机亚洲免费影院| 欧美日韩黄片免| 波多野结衣av一区二区av| 亚洲第一av免费看| 日本黄色日本黄色录像| 中文字幕人妻丝袜制服| 最新的欧美精品一区二区| 波多野结衣一区麻豆| 最新在线观看一区二区三区| 黄片播放在线免费| 一夜夜www| 大香蕉久久成人网| 日本vs欧美在线观看视频| 国产1区2区3区精品| xxxhd国产人妻xxx| 亚洲欧美一区二区三区久久| 国产精品亚洲一级av第二区| 精品人妻1区二区| 人妻一区二区av| 亚洲精品中文字幕一二三四区| 日本vs欧美在线观看视频| 老司机午夜福利在线观看视频| 亚洲精品久久午夜乱码| 成年版毛片免费区| 久久久久久久国产电影| 两个人免费观看高清视频| 亚洲精品美女久久av网站| 免费女性裸体啪啪无遮挡网站| 国产精品 国内视频| e午夜精品久久久久久久| 制服人妻中文乱码| 欧美人与性动交α欧美软件| 精品一区二区三卡| 精品国产一区二区三区久久久樱花| 男人操女人黄网站| 99久久综合精品五月天人人| 免费在线观看日本一区| 变态另类成人亚洲欧美熟女 | www.自偷自拍.com| 国产不卡av网站在线观看| 大香蕉久久网| 久久精品亚洲av国产电影网| 国产成人啪精品午夜网站| 好看av亚洲va欧美ⅴa在| 不卡一级毛片| 国产又爽黄色视频| 国产一区二区激情短视频| 亚洲七黄色美女视频| 欧美日韩视频精品一区| 国产男女内射视频| 搡老熟女国产l中国老女人| 国产精品综合久久久久久久免费 | 热re99久久精品国产66热6| 在线观看一区二区三区激情| av视频免费观看在线观看| 丝袜美腿诱惑在线| 国产熟女午夜一区二区三区| 成人三级做爰电影| av电影中文网址| 精品一区二区三卡| 国产成人精品在线电影| 国产精华一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 香蕉久久夜色| 人人妻,人人澡人人爽秒播| 亚洲av美国av| 久久天躁狠狠躁夜夜2o2o| 中出人妻视频一区二区| 精品高清国产在线一区| 黑人巨大精品欧美一区二区mp4| 99re在线观看精品视频| 大陆偷拍与自拍| 大香蕉久久成人网| 久久香蕉激情| 9色porny在线观看| 国产乱人伦免费视频| 精品一区二区三卡| 90打野战视频偷拍视频| 后天国语完整版免费观看| 国产日韩一区二区三区精品不卡| 亚洲五月色婷婷综合| 国产精品免费一区二区三区在线 | 久久久久久久午夜电影 | 午夜成年电影在线免费观看| 女性被躁到高潮视频| 韩国精品一区二区三区| 国产区一区二久久| 国产精品免费大片| 亚洲精品在线美女| 性色av乱码一区二区三区2| 亚洲欧美日韩高清在线视频| 啦啦啦免费观看视频1| 欧美日韩一级在线毛片| 色老头精品视频在线观看| 天堂√8在线中文| 少妇的丰满在线观看| 国产精品久久久久成人av| 日韩欧美一区二区三区在线观看 | 成人18禁在线播放| av中文乱码字幕在线| 国产欧美日韩一区二区精品| 欧美日韩精品网址| 久久久国产精品麻豆| 色老头精品视频在线观看| 一本一本久久a久久精品综合妖精| 老司机午夜十八禁免费视频| 九色亚洲精品在线播放| 国产一区二区激情短视频| 亚洲avbb在线观看| 大型黄色视频在线免费观看| 欧美精品人与动牲交sv欧美| 99国产精品99久久久久| 亚洲片人在线观看| 国产欧美日韩综合在线一区二区| 一区二区三区国产精品乱码| 久久狼人影院| 欧美精品人与动牲交sv欧美| 黑人猛操日本美女一级片| 变态另类成人亚洲欧美熟女 | 国产精品亚洲一级av第二区| 久久精品国产亚洲av香蕉五月 | 最新美女视频免费是黄的| 在线观看免费视频日本深夜| 伊人久久大香线蕉亚洲五| 波多野结衣av一区二区av| 久久人妻福利社区极品人妻图片| 亚洲成国产人片在线观看| 亚洲国产精品合色在线| 国产精品久久久av美女十八| 99riav亚洲国产免费| 日日爽夜夜爽网站| 久久久久精品人妻al黑| 免费观看a级毛片全部| 一级毛片精品| 国产高清视频在线播放一区| 成年人免费黄色播放视频| 中文字幕人妻熟女乱码| 久久人人爽av亚洲精品天堂| 午夜成年电影在线免费观看| a级毛片黄视频| 另类亚洲欧美激情| 在线观看午夜福利视频| 黄色a级毛片大全视频| 久久久久久久国产电影| 宅男免费午夜| 男女之事视频高清在线观看| 淫妇啪啪啪对白视频| 一进一出抽搐gif免费好疼 | 国产成人免费无遮挡视频| 好男人电影高清在线观看| 日韩视频一区二区在线观看| 久久草成人影院| 亚洲专区中文字幕在线| 国产有黄有色有爽视频| 老司机福利观看| 国产精品1区2区在线观看. | 午夜影院日韩av| 国产成人精品久久二区二区91| 丰满人妻熟妇乱又伦精品不卡| videos熟女内射| 19禁男女啪啪无遮挡网站| 18在线观看网站| 女性被躁到高潮视频| 叶爱在线成人免费视频播放| 黑人猛操日本美女一级片| 91精品国产国语对白视频| 亚洲全国av大片| 99国产综合亚洲精品| 午夜成年电影在线免费观看| 免费在线观看影片大全网站| 视频区欧美日本亚洲| 老司机午夜福利在线观看视频| 欧美+亚洲+日韩+国产| 国产伦人伦偷精品视频| 99热国产这里只有精品6| 午夜91福利影院| 婷婷成人精品国产| 久久国产精品大桥未久av| 别揉我奶头~嗯~啊~动态视频| netflix在线观看网站| 国产在线观看jvid| 51午夜福利影视在线观看| 黄色丝袜av网址大全| 成年版毛片免费区| 黄频高清免费视频| 欧美大码av| 国产在线一区二区三区精| 亚洲一码二码三码区别大吗| 高清av免费在线| 婷婷成人精品国产| 亚洲人成电影观看| 午夜91福利影院| 日韩精品免费视频一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲午夜理论影院| 人人妻,人人澡人人爽秒播| 日韩有码中文字幕| 亚洲国产精品合色在线| 欧美日本中文国产一区发布| 国产成人欧美在线观看 | 老司机靠b影院| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人 | 国产精品免费视频内射| 大香蕉久久网| 美女扒开内裤让男人捅视频| 悠悠久久av| 美女高潮到喷水免费观看| 国产亚洲欧美精品永久| 久久久久久免费高清国产稀缺| 黄频高清免费视频| 夜夜躁狠狠躁天天躁| 亚洲av成人一区二区三| av视频免费观看在线观看| 国产成人一区二区三区免费视频网站| 久久国产亚洲av麻豆专区| 啦啦啦 在线观看视频| 校园春色视频在线观看| 999久久久国产精品视频| 夜夜爽天天搞| 91老司机精品| 亚洲熟妇中文字幕五十中出 | 91九色精品人成在线观看| 欧美黑人精品巨大| x7x7x7水蜜桃| 黑人巨大精品欧美一区二区蜜桃| 一级,二级,三级黄色视频| 亚洲av成人一区二区三| 91麻豆精品激情在线观看国产 | 国产成人欧美| 一级毛片精品| 久久久久精品人妻al黑| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| 人人妻人人爽人人添夜夜欢视频| av网站在线播放免费| 人人妻,人人澡人人爽秒播| 欧美乱妇无乱码|