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

    寧麥9號與揚麥158株高及其構(gòu)成因素的遺傳解析

    2020-05-21 00:43:42張平平馬鴻翔
    作物學(xué)報 2020年6期

    姜 朋 何 漪 張 旭 吳 磊 張平平 馬鴻翔

    寧麥9號與揚麥158株高及其構(gòu)成因素的遺傳解析

    姜 朋 何 漪 張 旭 吳 磊 張平平 馬鴻翔*

    江蘇省農(nóng)業(yè)科學(xué)院 / 江蘇省農(nóng)業(yè)生物學(xué)重點實驗室 / 江蘇省現(xiàn)代作物生產(chǎn)協(xié)同創(chuàng)新中心, 江蘇南京 210014

    寧麥9號與揚麥158是我國長江中下游麥區(qū)的主栽品種和骨干親本, 長江中下游麥區(qū)近3年來審定品種中80%都是其衍生后代, 研究其性狀的遺傳具重要意義。以寧麥9號與揚麥158為親本構(gòu)建的包含282個家系的重組自交系群體為材料, 利用Illumina 90k芯片對群體進(jìn)行基因型分析, 建立高密度遺傳圖譜。連續(xù)3個生長季對株高及節(jié)間長度、穗長等株高構(gòu)成因素進(jìn)行測定, 結(jié)合遺傳圖譜對株高及相關(guān)性狀進(jìn)行QTL定位, 獲得14個控制株高及其構(gòu)成因素的穩(wěn)定表達(dá)位點。通過進(jìn)一步位置比對, 聚焦到6個染色體區(qū)段, 初步明確了各節(jié)間對株高的遺傳調(diào)控機(jī)制。同時, 將6個染色體區(qū)段中同源性較低的連鎖標(biāo)記轉(zhuǎn)化為適用于高通量篩選的KASP標(biāo)記, 利用101份區(qū)域試驗材料進(jìn)行標(biāo)記效應(yīng)驗證, 結(jié)果顯示聚合與兩個位點具有較高的選擇效率, 繼續(xù)聚合后, 中選材料顯著減少, 可能降低選擇效率; 對與兩個一因多效位點的選擇建議以降低株高的等位變異為主;可作為穗下節(jié)間(D1)的選擇標(biāo)記對株高展開優(yōu)化選擇。期望以上結(jié)果能為長江中下游麥區(qū)的小麥株高遺傳改良提供幫助。

    小麥; 寧麥9號; 揚麥158; 株高; KASP標(biāo)記

    株高是小麥重要的農(nóng)藝性狀, 影響植株的形態(tài)結(jié)構(gòu), 并與田間群體產(chǎn)量密切相關(guān)。小麥矮稈基因的利用是綠色革命的重要組成部分, 對現(xiàn)代小麥育種具有深遠(yuǎn)影響[1]。經(jīng)典遺傳學(xué)研究表明, 小麥株高是一個復(fù)雜性狀, 由多個基因控制, 存在主效基因, 也有微效位點。迄今為止, 已有25個基因被命名[2-4]。位于 4B 和 4D 染色體上的和基因, 以及2D染色體上的基因在世界范圍內(nèi)廣為應(yīng)用, 其相關(guān)分子標(biāo)記已成功開發(fā)[5-6]。此外, 在小麥21條染色體上均檢測到影響株高的QTL[7-11]。McCartney等[12]利用RL4452בAC’ Domain構(gòu)建的DH群體對包括株高在內(nèi)的多個農(nóng)藝性狀進(jìn)行QTL分析, 在 2D、 4B、4D、5B、7A和7B染色體檢測到6個株高QTL, 其中和的定位區(qū)間與和位置重合。Griffiths等[13]利用4個DH群體進(jìn)行了株高meta-QTL分析, 在除3D、4A、5D及第7同源群外的15條染色體上找到16個meta-QTL。Liu等[9]在5個生長階段對小麥株高進(jìn)行了跟蹤調(diào)查, 利用條件與非條件QTL作圖方法進(jìn)行分析, 檢測到8個條件QTL與9個非條件QTL。其中在2個時期被重復(fù)檢測到, 表型貢獻(xiàn)率達(dá)13.42%~16.13%, 此研究表明株高基因的表達(dá)具有一定的時空特異性。小麥株高是主穗長與各節(jié)間長的總和, 但控制主穗長及各莖節(jié)間長的遺傳位點對小麥株高的遺傳貢獻(xiàn)并不一致。Cui等[14]利用條件QTL與非條件QTL相結(jié)合的分析方法在單個QTL水平揭示了株高與穗長及各節(jié)間長之間的關(guān)系, 其中倒三節(jié)對株高影響最大。Zhang等[15]測定了不同時期的各個節(jié)間長度, 明確了5個穩(wěn)定的株高QTL的表達(dá)模式及貢獻(xiàn)率。

    遺傳連鎖圖為遺傳定位、標(biāo)記開發(fā)及候選基因發(fā)掘等提供了強(qiáng)有力的工具。近年來發(fā)展起來的SNP標(biāo)記具有遺傳穩(wěn)定、數(shù)量多、分布廣等特點, 并且適于高通量檢測, 基于其開發(fā)的9k、90k、660k等基因型芯片集合了成千上萬個SNP標(biāo)記, 不僅加快了小麥遺傳圖譜構(gòu)建速度, 也大大提高了小麥遺傳圖譜的密度[16-18]。中國春小麥參考基因組序列的公布, 有助于小麥遺傳圖譜質(zhì)量的提高, 并可直接基于目標(biāo)區(qū)段進(jìn)行候選基因預(yù)測[19]。

    長江中下游麥區(qū)是我國第二大麥區(qū), 同時也是我國最大的弱筋小麥產(chǎn)區(qū), 是我國小麥生產(chǎn)的重要組成部分。寧麥9號與揚麥158分別是由江蘇省農(nóng)業(yè)科學(xué)院與江蘇里下河地區(qū)農(nóng)業(yè)科學(xué)研究所育成的高產(chǎn)優(yōu)質(zhì)抗病小麥品種, 都曾經(jīng)是生產(chǎn)上的主栽品種, 具有較大的推廣面積, 而且目前仍是重要的骨干親本, 以此為基礎(chǔ)育成了數(shù)十個小麥品種[20-21]。寧麥9號的衍生品種寧麥13、揚輻麥4號等, 揚麥158的衍生品種揚麥20等均已成為當(dāng)前生產(chǎn)中的主栽品種及常用親本。在生產(chǎn)中, 寧麥9號株高一般在80~ 85 cm, 揚麥158約為90~95 cm, 呈現(xiàn)較大差異, 明確其株高遺傳機(jī)制對小麥育種工作具有重要的指導(dǎo)意義。本研究以來源于寧麥9號/揚麥158的282份重組自交系為材料, 結(jié)合高密度遺傳圖譜, 對株高及其構(gòu)成開展QTL定位研究, 以期為長江中下游麥區(qū)育種工作中的株高選擇提供幫助。

    1 材料與方法

    1.1 試驗材料

    以寧麥9號×揚麥158構(gòu)建重組自交系(RIL)群體(F2:8), 包括282個家系。寧麥9號是揚麥6號與日本西風(fēng)小麥的雜交后代, 前者具有骨干親本南大2419與江東門的遺傳背景; 揚麥158來源于St1472/506與揚麥4號雜交組合, St1472/506是黃淮麥區(qū)的常用親本, 揚麥4號含有南大2419、勝利麥與阿夫等多個骨干親本的遺傳背景。寧麥9號與揚麥158株高及各節(jié)間長度、穗長等性狀差異顯著(圖1)。通過Ellis等[5]和Asplund等[6]報道的分子標(biāo)記確定寧麥9號與揚麥158均為/變異類型, 均不含位點(圖2), 其株高差異可能來源于其他位點。

    圖1 寧麥9號與揚麥158的田間表現(xiàn)

    圖2 寧麥9號與揚麥158株高相關(guān)分子標(biāo)記檢測

    a: 寧麥9號; b: 揚麥158; M: marker。

    a: Ningmai 9; b: Yangmai 158; M: marker.

    1.2 田間試驗及數(shù)據(jù)處理

    2016—2017、2017—2018和2018—2019連續(xù)3個生長季將RIL群體及其親本種植于江蘇省農(nóng)業(yè)科學(xué)院六合基地, 為方便描述, 以收獲年份2017、2018和2019分別表示3個環(huán)境。采用隨機(jī)區(qū)組設(shè)計, 單行種植, 每行60粒, 行長1.6 m, 行距0.25 m, 2次重復(fù), 常規(guī)田間管理。

    于小麥乳熟期從每個家系中取10個單莖, 從莖基部開始向上利用直尺測量每個節(jié)間-倒五節(jié)(The fifth internode from the top, D5)、倒四節(jié)(The fourth internode from the top, D4)、倒三節(jié)(The third internode from the top, D3)、倒二節(jié)(The second internode from the top, D2)、倒一節(jié)(The first internode from the top, D1)長度及穗長, 株高為各節(jié)間與穗長之和, 并計算平均值用于進(jìn)一步統(tǒng)計分析。另同時測量101份2018國家及省區(qū)域試驗材料的株高, 用于后續(xù)驗證。

    采用Microsoft Excel 2016進(jìn)行表型初步統(tǒng)計與相關(guān)分析, 采用SPSS 19.0進(jìn)行方差分析及檢測。依據(jù)He等[22]的方法, 按2=δ2/(δ2+δ2/+δ/)計算遺傳力, 其中δ2表示基因型方差,δ2表示基因型與環(huán)境互作方差,δ2為誤差,代表環(huán)境數(shù)目,代表重復(fù)數(shù)。

    1.3 遺傳圖譜與QTL分析

    采用Illumina 90k芯片獲取基因型。遺傳圖譜覆蓋21條染色體, 包含41個連鎖群, 2285個bin標(biāo)記, 總長為3022 cM[23]。采用QTL IciMapping 4.1軟件[24]的完備區(qū)間作圖法(inclusive composite interval mapping, ICIM)進(jìn)行QTL定位[25], 設(shè)walking step為0.1 cM, LOD閾值為2.5。

    1.4 KASP標(biāo)記的開發(fā)

    根據(jù)SNP位點和側(cè)翼序列設(shè)計PCR擴(kuò)增引物, 開發(fā)KASP分子標(biāo)記。設(shè)計每個標(biāo)記2條SNP特異性引物(F1/F2)和一條通用引物(R), F1尾部添加能夠與FAM熒光結(jié)合的特異性序列, F2尾部添加能夠與HEX熒光結(jié)合的特異性序列。利用Polymarker (http://www.polymarker.info/)設(shè)計KASP引物, 由生工生物工程(上海)股份有限公司合成。

    從RILs群體中隨機(jī)挑選46份材料進(jìn)行幼嫩葉片取樣, 采用CTAB法[25]提取基因組DNA。KASP反應(yīng)總體系為5 μL, 包含2×KASP Master Mix 2.5 μL、KASP Assay Mix (引物混合工作液) 0.07 μL、濃度為20 ng μL–1的模板DNA 2.43 μL。KASP反應(yīng)程序第一步為94℃, 15 min; 第二步為94℃, 20 s, 61~55℃, 1 min, 每個循環(huán)降低0.6℃, 共進(jìn)行10個循環(huán); 第三步為94℃, 20 s, 55℃, 1 min, 共進(jìn)行26個循環(huán), 在LGC公司為Hydrocycler-16水浴PCR儀中進(jìn)行PCR。通過KASP熒光分析儀(LGC公司型號為PHERAstar plus)掃描分析PCR結(jié)果。

    2 結(jié)果與分析

    2.1 株高及其構(gòu)成因素的表型分析

    在連續(xù)3年的調(diào)查中(表1), 揚麥158的各節(jié)間長度、穗長及株高均高于寧麥9號, 其中D5與D4差異較小, 一般為1~2 cm, D3、D2與D1差異較大, 約為2~3 cm, 穗長相差約1 cm, 最終造成株高相差約10 cm。在RIL群體中, 各節(jié)間長度、穗長及株高均呈現(xiàn)較大變異, D5變異系數(shù)最大, 約20%, 其余節(jié)間、穗長及株高變異系數(shù)在10%左右, 最大株高與最小株高差異均超過30 cm, 節(jié)間差距從D5的5 cm左右遞增至D1的近20 cm。從D5至D1, 性狀遺傳力從0.43逐步增大至0.83, 穗長與株高遺傳力均超過0.7。

    表1 株高及其構(gòu)成因素的表型統(tǒng)計

    D5: the fifth internode from the top; D4: the fourth internode from the top; D3: the third internode from the top; D2: the second internode from the top; D1: the first internode from the top.

    株高、穗長及各節(jié)間長度在不同基因型與不同年份間差異極顯著, 同時基因型與年份的互作對其也有極顯著影響(表2)。各節(jié)間長度間均呈極顯著正相關(guān), 而與穗長相關(guān)系數(shù)普遍較低。株高與各節(jié)間長度、穗長均呈極顯著正相關(guān), 而同一節(jié)間在不同環(huán)境間亦呈現(xiàn)出良好的相關(guān)性(表3)。

    表2 株高及其構(gòu)成因素的方差分析(F值)

    *和**分別表示0.05和0.01顯著水平。

    *and**indicate significant difference at the 0.05 and 0.01 probability levels, respectively. Abbreviations are the same as those given in Table 1.

    表3 株高及其構(gòu)成因素的相關(guān)分析

    *和**分別表示0.05和0.01顯著水平; 對角線上括號內(nèi)數(shù)字為同一性狀不同環(huán)境間的相關(guān)系數(shù), 2017為2017年度與2018年度間相關(guān)系數(shù), 2018為2018年度與2019年度間相關(guān)系數(shù), 2019為2019年度與2017年度間相關(guān)系數(shù)。

    *and**indicate significant difference at the 0.05 and 0.01 probability levels, respectively. The figures in brackets on the diagonal line indicate correlation coefficients between the same traits in different environments. The correlation coefficients between 2017 and 2018 are put in the line of 2017, and those between 2018 and 2019, 2019 and 2017 are put in lines of 2018 and 2019, respectively. Abbreviations are the same as those given in Table 1.

    2.2 株高及其構(gòu)成因素的QTL分析

    綜合3個環(huán)境表型數(shù)據(jù), 利用完備區(qū)間作圖法共鑒定到96個QTL, 分布在18條染色體上, 表型貢獻(xiàn)率為1.72%~12.32%。其中, D5至D1分別檢測到7、14、17、13和13個QTL, 與穗長相關(guān)的QTL有11個, 控制最終株高的QTL有21個。進(jìn)一步對這些QTL比較分析, 最終在2A、2D、5A、5D、7A等5條染色體上獲得14個多環(huán)境表達(dá)的穩(wěn)定QTL (圖2和表4)。其中, 控制D5、D4及D3長度的QTL均為1個, 控制D2與D1均為3個, 而控制株高的為5個, 沒有檢測到與穗長相關(guān)的穩(wěn)定QTL。

    通過位置比對發(fā)現(xiàn),、、、、和遺傳位置接近或重合, 很可能為1個一因多效QTL, 同時控制多個節(jié)間的長度。、與,與, 均是相同的情況。為方便描述, 將與分別簡寫為、與。除此3個QTL外,對D1具有調(diào)節(jié)作用,同時對D4與D3產(chǎn)生影響。以株高及各節(jié)間長度減少定義其加性效應(yīng), 這3個位點對應(yīng)多個性狀的加性效應(yīng)表現(xiàn)一致,的加性效應(yīng)來源于揚麥158, 而與的加性效應(yīng)來源于寧麥9號。此外,加性效應(yīng)來源于寧麥9號,與加性效應(yīng)來源于揚麥158。

    2.3 KASP標(biāo)記開發(fā)

    根據(jù)定位結(jié)果, 我們最終獲得6個株高及其構(gòu)成因素控制區(qū)段(圖3), 每個區(qū)段內(nèi)含有多個SNP標(biāo)記, 選擇其中序列同源性較低的標(biāo)記進(jìn)行KASP標(biāo)記開發(fā), 最終將IAAV880 ()、IACX9152 ()、BobWhite_c1796_701 ()、BS00098207_51 ()、Kukri_c7786_81 ()與wsnp_JD_c6050_ 7214383 ()等成功轉(zhuǎn)化(表5和圖4)。

    圖3 株高及其構(gòu)成因素的QTL定位

    Fig. 3 QTL mapping for plant height and its components

    圖中數(shù)字與表4中QTL的序號一致。

    The numbers in the figure are consistent with the sequence numbers of QTL in Table 4.

    2.4 KASP標(biāo)記的驗證

    為檢測已開發(fā)KASP標(biāo)記的功效, 利用101份2018國家及省區(qū)域試驗材料進(jìn)行驗證。結(jié)果顯示, 除與外, 其他4個標(biāo)記的加性效應(yīng)與定位結(jié)果一致(表6)。對于與位點, 僅有不到10%的檢測材料中含有寧麥9號等位變異, 這可能影響統(tǒng)計檢測結(jié)果的準(zhǔn)確性; 而主要調(diào)控D1長度, 對最終株高影響不顯著。在其他3個位點中,對株高有顯著影響, 可產(chǎn)生近2 cm的株高差異,、分別產(chǎn)生0.59 cm、1.32 cm的差異, 但未達(dá)到顯著水平。

    表5 株高及其構(gòu)成因素的KASP標(biāo)記序列

    F1、F2為SNP特異性引物, F1尾部添加能夠與FAM熒光結(jié)合的特異性序列, F2尾部添加能夠與HEX熒光結(jié)合的特異性序列; R為通用引物。

    F1 and F2 are the specific primers for the SNP, specific sequence with FAM binding fluorescence is added to the end of F1, and that with HEX binding fluorescence is added to the end of F2; R is common primer.

    圖4 株高及其構(gòu)成因素的KASP標(biāo)記開發(fā)

    Fig, 4 Development of KASP markers for plant height and its components

    A表示寧麥9號等位變異類型, B表示揚麥158等位變異類型。

    A indicates the allele of Ningmai 9 and B indicates the allele of Yangmai 158.

    表6 區(qū)域試驗材料中不同等位變異的t-檢測

    *表示0.05顯著水平; 括號中數(shù)字表示攜帶相應(yīng)等位變異的材料數(shù)量。

    *indicates significant difference at the 0.05 probability level. The figures in brackets indicate the number of the materials with corresponding alleles.

    按照差異大小逐步分析幾個位點的聚合效應(yīng)發(fā)現(xiàn), 當(dāng)聚合與位點時, 株高可相差4 cm, 進(jìn)一步聚合后, 差異超過10 cm, 但材料數(shù)量急劇減少, 從降低株高的角度出發(fā), 結(jié)合與較單一分子標(biāo)記選擇更加有效, 當(dāng)聚合更多標(biāo)記后, 株高降低不顯著, 而入選材料有所減少, 選擇效率降低。

    表7 QTL聚合的效應(yīng)分析

    “+”表示增加株高的等位變異組合, “?”表示降低株高的等位變異組合;*表示0.05顯著水平; 括號中數(shù)字表示攜帶相應(yīng)等位變異的材料數(shù)量。

    “+” indicates the combination of the alleles which increase plant height, and “?” indicates the combination of the alleles which reduce plant height;*indicates significant difference at the 0.05 probability level; the figures in brackets indicate the number of the materials with corres-ponding alleles.

    3 討論

    寧麥9號與揚麥158分別于1997年和1993年通過品種審定, 育成后均在生產(chǎn)上大面積種植, 為長江中下游麥區(qū)乃至全國具有重大影響力的品種。揚麥158適應(yīng)性廣, 且具有良好的農(nóng)藝性狀和較好的赤霉病抗性[27], 寧麥9號具有突出的弱筋小麥品質(zhì), 且多花多粒、結(jié)實性高、豐產(chǎn)性好, 還是我國小麥栽培品種赤霉病抗性QTL-的主要供體[28]。這2個品種已經(jīng)育成超過20年, 在生產(chǎn)中已被后代品種替代, 但它們目前仍被作為雜交親本應(yīng)用于小麥育種。當(dāng)前小麥生產(chǎn)中大面積推廣的寧麥13、揚輻麥4號、揚麥20等均為寧麥9號或揚麥158的衍生后代。據(jù)統(tǒng)計, 近3年江蘇省淮南麥區(qū)與國家長江中下游麥區(qū)通過審定的小麥品種中, 超過80%是寧麥9號或揚麥158的衍生后代(http://www.jsseed. cn/, http://www.moa.gov.cn/)。此外, 作為我國小麥品種中的主要來源, 寧麥9號已被引入其他麥區(qū)作為赤霉病抗性親本加以利用[29]; 輪選16、輪選66及輪選166等國家黃淮南片麥區(qū)的審定品種均含有揚麥158的遺傳背景。鑒于此, 針對其重要農(nóng)藝性狀開展遺傳研究對更好地在育種中利用這2個材料具實踐意義。

    寧麥9號與揚麥158株高差異較大。本研究顯示揚麥158各節(jié)間長度與穗長均大于寧麥9號, 同時遺傳因素在株高及其構(gòu)成因素中占主導(dǎo)作用, 與前人研究結(jié)果一致[30-31]。我們通過分子標(biāo)記檢測確定寧麥9號與揚麥158均為/變異類型, 均不含位點。對遺傳定位結(jié)果比對分析, 最終將影響株高的QTL聚焦到6個染色體區(qū)段。通過物理位置比對,位于2D染色體短臂中部, 與位于2D染色體短臂末端的[32]相距較遠(yuǎn), 而與Wang等[33]利用條件與非條件QTL作圖檢測到的接近。梁子英等[34]在5A染色體–區(qū)段檢測到株高、D4與D5的控制位點, 其位置與基本重合, 而在2019環(huán)境中與D4、D3也表現(xiàn)相關(guān)。位于5A染色體長臂, 與的關(guān)聯(lián)標(biāo)記位置接近, 并且來源于日本小麥品種, 而亦來源于含有日本西風(fēng)小麥血統(tǒng)的寧麥9號, 推測可能為同一位點[32]。位于5D染色體, 其連鎖標(biāo)記、均與相近, 但由于5D染色體標(biāo)記較少, 覆蓋度不足, 是否為同一位點有待進(jìn)一步驗證[32]。前人在2A染色體上鑒定到多個株高控制位點[10,33-34], 但與均相距較遠(yuǎn)。Peng等[35]在7A染色體短臂末端–區(qū)間定位到, 而位于7A染色體短臂中部, 可能為不同的位點。

    在小麥育種中, 一般將株高較低的基因型作為優(yōu)勢等位變異。與在101份區(qū)域試驗材料中呈現(xiàn)較強(qiáng)的選擇效應(yīng), 多數(shù)材料在位點處呈現(xiàn)來源于揚麥158的優(yōu)勢等位變異,與所有節(jié)間長度都有關(guān), 可能為一個非常重要的選擇位點, 并且很可能經(jīng)歷過早代選擇, 使其在高代材料中分布頻率較低; 針對此位點, 今后可在低代材料中進(jìn)行驗證, 并同時構(gòu)建次級分離群體, 為進(jìn)一步精細(xì)定位奠定基礎(chǔ)。通過比較基因組學(xué)分析發(fā)現(xiàn),附近區(qū)域為許多小麥重要農(nóng)藝、抗性、品質(zhì)性狀控制位點的熱點聚集區(qū), 包括開花期、葉枯病、黃色素含量等[36-39], 對其中某些性狀的定向選擇促使某種基因型被大量保留, 最終導(dǎo)致呈現(xiàn)大量的非優(yōu)勢等位變異。

    各節(jié)間長度與不同的農(nóng)藝性狀具有一定的相關(guān)性[40-41], 特別是植株抗倒性, 研究表明矮稈、D5較短、D1較長的植株有利于協(xié)調(diào)高產(chǎn)與抗倒性[42]。本研究顯示同時調(diào)控D1與株高,同時調(diào)控D5、D1與株高, 且其對應(yīng)各性狀的加性效應(yīng)來源一致, 無法通過這2個位點實現(xiàn)D1與D5、株高的差異化選擇。前述在多數(shù)區(qū)域試驗材料中呈現(xiàn)來源于揚麥158的降低株高、縮短各節(jié)間長度的等位變異,也是多數(shù)保留了來源于寧麥9號的降低株高的等位變異, 表明育種家在株高選擇中仍以整株株高為主要選擇目標(biāo), 未考慮各節(jié)間的差異性選擇, 在選擇方向上對這2個一因多效位點的選擇以降低株高的等位變異為主。僅控制D1長度, 且對整體株高影響不大(表4和表6), 可作為D1的選擇標(biāo)記。

    4 結(jié)論

    共定位到6個株高控制區(qū)段, 并成功轉(zhuǎn)化為適用于大規(guī)模篩選的KASP標(biāo)記, 經(jīng)過初步驗證, 聚合、標(biāo)記對整體株高具有較高的選擇效率; 繼續(xù)聚合后, 中選材料顯著減少, 可能降低選擇效率, 而對需要在早代育種材料進(jìn)一步驗證, 并且對與一因多效位點的選擇建議以降低株高的等位變異為主;可作為D1的選擇標(biāo)記對株高進(jìn)行優(yōu)化選擇。

    [1] Hedden P. The genes of the green revolution., 2003, 19: 5–9.

    [2] Mo Y, Vanzetti L S, Hale I, Spagnolo E J, Guidobaldi F, Al-Oboudi J, Odle N, Pearce S, Helguera M, Dubcovsky J. Identification and characterization of, a locus on chromosome arm 6AS affecting wheat plant height, heading time, and spike development., 2018, 131: 2021–2035.

    [3] Tian X, Wen W, Xie L, Fu L, Xu D, Fu C, Wang D, Chen X, Xia X, Chen Q, He Z, Cao S. Molecular mapping of reduced plant height genein bread wheat., 2017, https://doi.org/10.3389/fpls.2017.01379.

    [4] McIntosh R A, Dubcovsky J, Rogers W J, Morris C, Xia X C. Catalogue of gene symbols for wheat: 2017 supplement. https://shigen.nig.ac.jp/wheat/komugi/genes/macgene/supplement2017. pdf. Accessed 17 Feb 2018.

    [5] Ellis M, Spielmeyer W, Gale K, Rebetzke G, Richards R. “Perfect” markers for theanddwarfing genes in wheat., 2002, 105: 1038–1042.

    [6] Asplund L, Leino M W, Hagenblad J. Allelic variation at thelocus in a 19th century wheat collection., 2012, 2012: 385610.

    [7] B?rner A, Schumann E, Fürste A, C?ster H, Leithold B, R?der S, Weber E. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (L.)., 2002, 105: 921–936.

    [8] Peng J, Yefim R, Tzion F, RiDer M S, Youchun L, Eviatar N, Abraham K. Domestication quantitative trait loci in, the progenitor of wheat., 2003, 100: 2489–2494.

    [9] Liu G, Xu S B, Ni Z F, Xie C J, Qin D D, Li J, Lu L H, Zhang J P, Peng H R, Sun Q X. Molecular dissection of plant height QTLs using recombinant inbred lines from hybrids between common wheat (L.) and spelt wheat (L.)., 2011, 56: 1897–1903.

    [10] Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Snape J. Meta-QTL analysis of the genetic control of crop height in elite European winter wheat germplasm., 2012, 29: 159–171.

    [11] Tobias W, Langer S M, Longin C F H. Genetic control of plant height in European winter wheat cultivars., 2015, 128: 865–874.

    [12] Mccartney C A, Somers D J, Humphreys D G, Lukow O, Ames N, Noll J, Cloutier S, McCallum B D. Mapping quantitative trait loci controlling agronomic traits in the spring wheat cross Rl4452x ‘AC Domain’., 2005, 48: 870.

    [13] Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Herry L, Faure S, Laurie D, Bilham L, Snape J. Meta-QTL analysis of the genetic control of ear emergence in elite European winter wheat germplasm., 2009, 119: 383–395.

    [14] Cui F, Li J, Ding A, Zhao C, Wang L, Wang X, Li S, Bao Y, Li X, Feng D, Kong L, Wang H. Conditional QTL mapping for plant height with respect to the length of the spike and internode in two mapping populations of wheat., 2011, 122: 1517–1536.

    [15] Zhang N, Fan X, Cui F, Zhao C, Zhang W, Zhao X, Yang L, Pan R, Chen M, Han J. Characterization of the temporal and spatial expression of wheat (L.) plant height at the QTL level and their influence on yield-related traits., 2017, 130: 1235–1252.

    [16] Cavanagh C R, Shiaoman C, Shichen W, Bevan Emma H, Stuart S, Seifollah K, Kerrie F, Cyrille S, Brown-Guedira G L, Alina A. Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars., 2013, 110: 8057–8062.

    [17] Wang S, Wong D, Forrest K, Allen A, Chao S, Huang B E, Maccaferri M, Salvi S, Milner S G, Cattivelli L. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array., 2014, 12: 787–796.

    [18] Cui F, Zhang N, Fan X, Zhang W, Zhao C, Yang L, Pan R, Chen M, Han J, Zhao X. Utilization of a wheat660k SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number., 2017, 7: 3788.

    [19] Consortium T I W G S. Shifting the limits in wheat research and breeding using a fully annotated reference genome., 2018, 361: eaar7191.

    [20] 張曉, 張伯橋, 江偉, 呂國鋒, 張曉祥, 李曼, 高德榮. 揚麥系列品種品質(zhì)性狀相關(guān)基因的分子檢測. 中國農(nóng)業(yè)科學(xué), 2015, 48: 3779–3793. Zhang X, Zhang B Q, Jiang W, Lyu G F, Zhang X X, Li M, Gao D R. Molecular detection for quality traits-related genes in Yangmai series wheat cultivars., 2015, 48: 3779–3793 (in Chinese with English abstract).

    [21] Jiang P, Zhang P P, Zhang X, Ma H X. Genetic contribution of Ningmai 9 wheat to its derivatives evaluated by using SNP markers., 2016, 2016: 3602986.

    [22] He X, Lillemo M, Shi J, Wu J, Bj?rnstad ?, Belova T, Dreisigacker S, Duveiller E, Singh P. QTL characterization ofhead blight resistance in CIMMYT bread wheat line Soru#1., 2016, 11: e0158052.

    [23] Jiang P, Zhang X, Wu L, He Y, Zhuang W, Cheng X, Ge W, Ma H, Kong L. A novel QTL on chromosome 5AL of Yangmai 158 increases resistance tohead blight in wheat., 2020, 69: 249–258.

    [24] Meng L, Li H H, Zhang L Y, Wang J K. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations., 2015, 3: 269–283.

    [25] Li H, Ye G, Wang J. A modified algorithm for the improvement of composite interval mapping., 2007, 175: 361–374.

    [26] Saghai-Maroof M A, Soliman K M, Jorgensen R A, Allard R W. Ribosomal DNA spacer-length polymorphisms in barley: mendelian inheritance, chromosomal location, and population dynamics., 1984, 81: 8014–8018.

    [27] 程順和, 張勇, 張伯橋, 高德榮, 吳宏亞, 陸成彬, 呂國鋒, 王朝順. 小麥抗赤霉病育種2條技術(shù)路線的探討. 揚州大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版), 2003, 24: 59–62. Cheng S H, Zhang Y, Zhang B Q, Gao D R, Wu H Y, Lu C B, Lyu G F, Wang C S. Discussion of two ways of breeding scab resistance in wheat.(Agric & Life Sci), 2003, 24: 59–62 (in Chinese with English abstract).

    [28] 朱展望, 徐登安, 程順和, 高春保, 夏先春, 郝元峰, 何中虎. 中國小麥品種抗赤霉病基因的鑒定與溯源. 作物學(xué)報, 2018, 44: 473–482. Zhu Z W, Xu D A, Cheng S H, Gao C B, Xia X C, Hao Y F, He Z H. Characterization ofhead blight resistance geneand its putative ancestor in Chinese wheat germplasm., 2018, 44: 473?482 (in Chinese with English abstract).

    [29] 張宏軍, 宿振起, 柏貴華, 張旭, 馬鴻翔, 李騰, 鄧云, 買春艷, 于立強(qiáng), 劉宏偉, 楊麗, 李洪杰, 周陽. 利用基因功能標(biāo)記選擇提高黃淮冬麥區(qū)小麥品種對赤霉病的抗性. 作物學(xué)報, 2018, 44: 505–511. Zhang H J, Su Z Q, Bai G H, Zhang X, Ma H X, Li T, Deng Y, Mai C Y, Yu L Q, Liu H W. Improvement of resistance of wheat cultivars tohead blight in the Yellow-Huai rivers valley winter wheat zone with functional marker selection ofgene., 2018, 44: 505?511 (in Chinese with English abstract).

    [30] 陳廣鳳, 陳建省, 田紀(jì)春. 小麥株高相關(guān)性狀與SNP標(biāo)記全基因組關(guān)聯(lián)分析. 作物學(xué)報, 2015, 41: 1500–1509. Chen G F, Chen J S, Tian J C. Genome-wide association analysis between SNP markers and plant height related traits in wheat., 2015, 41: 1500–1509 (in Chinese with English abstract).

    [31] 武炳瑾, 馮潔, 崔紫霞, 張傳量, 孫道杰. 利用90k基因芯片進(jìn)行小麥株高QTL分析. 麥類作物學(xué)報, 2017, 5: 578?584. Wu B J, Feng J, Cui Z X, Zhang C L, Sun D J. QTL analysis of plant height by using 90k chip technology., 2017, 5: 578–584 (in Chinese with English abstract).

    [32] Ellis M H, Rebetzke G J, Azanza F, Richards R A, Spielmeyer W. Molecular mapping of gibberellin-responsive dwarfing genes in bread wheat., 2005, 111: 423–430.

    [33] Wang Z H, Wu X S, Ren Q, Chang X P, Li R Z, Jing R L. QTL mapping for developmental behavior of plant height in wheat (L.)., 2010, 174: 447–458.

    [34] 梁子英, 李美霞, 王竹林, 沈瑋囡, 奚亞軍, 孫風(fēng)麗, 劉曙東. 小麥株高相關(guān)性狀的QTL分析. 西北農(nóng)業(yè)學(xué)報, 2014, 23: 64–72. Liang Z Y, Li M X, Wang Z L, Shen W N, Xi Y J, Sun F L, Liu S D. Mapping quantitative trait loci for plant height related characteristics in wheat.,2014, 23: 64–72 (in Chinese with English abstract).

    [35] Peng Z S, Li X, Yang Z J, Liao M L. A new reduced height gene found in the tetraploid semi-dwarf wheat landrace Aiganfanmai., 2011, 10: 2349.

    [36] Roncallo P F, Cervigni G L, Jensen C, Miranda R, Carrera A D, Helguera M, Echenique V. QTL analysis of main and epistatic effects for flour color traits in durum wheat., 2012, 185: 77–92.

    [37] Kamran A, Iqbal M, Spaner D. Flowering time in wheat (L.): a key factor for global adaptability., 2014, 197: 1–26.

    [38] Botwright Acu?a T L, Rebetzke G J, He X, Maynol E, Wade L J. Mapping quantitative trait loci associated with root penetration ability of wheat in contrasting environments., 2014, 34: 631–642.

    [39] Dreisigacker S, Wang X, Martinez Cisneros B A, Jing R, Singh P K. Adult-plant resistance toblotch in hexaploid spring wheat., 2015, 128: 2317–2329.

    [40] Kelbert A J, Spaner D, Briggs K G, King J R. The association of culm anatomy with lodging susceptibility in modern spring wheat genotypes., 2004, 136: 211–221.

    [41] 朱新開, 郭文善, 李春燕, 封超年, 彭永欣. 小麥株高及其構(gòu)成指數(shù)與產(chǎn)量及品質(zhì)的相關(guān)性. 麥類作物學(xué)報, 2009, 29: 1034?1038.Zhu X K, Guo W S, Li C Y, Feng C N, Peng Y X. Relationship of plant height component indexes with grain yield and quality in wheat., 2009, 29: 1034–1038 (in Chinese with English abstract).

    [42] 朱新開, 王祥菊, 郭凱泉, 郭文善, 封超年, 彭永欣. 小麥倒伏的莖稈特征及對產(chǎn)量與品質(zhì)的影響. 麥類作物學(xué)報, 2006, 26: 87–92. Zhu X K, Wang X J, Guo K Q, Guo W S, Feng C N, Peng Y X. Stem characteristics of wheat with stem lodging and effects of lodging on grain yield and quality., 2006, 26: 87–92 (in Chinese with English abstract).

    Genetic analysis of plant height and its components for wheat (L.) cultivars Ningmai 9 and Yangmai 158

    JIANG Peng, HE Yi, ZHANG Xu, WU Lei, ZHANG Ping-Ping, and MA Hong-Xiang*

    Jiangsu Academy of Agricultural Sciences / Jiangsu Provincial Key Laboratory for Agrobiology / Jiangsu Collaborative Innovation Center for Modern Crop Production, Nanjing 210014, Jiangsu, China

    Ningmai 9 and Yangmai 158 are the main wheat cultivars and core parents in the middle and lower reaches of the Yangtze River in China. In the past three years, 80% of the released varieties in the middle and lower reaches of the Yangtze River had the background of Ningmai 9 or Yangmai 158. To make better use of these two parents, the genetic mechanism of their traits need to be further clarified. A high-density genetic map was constructed by Illumina 90k chip using 282 recombinant inbred lines (RILs) from the cross between Ningmai 9 and Yangmai 158. In this study, the traits including plant height, internode length, and spike length were determined in three consecutive growing seasons, and 14 stable QTLs were obtained by QTL mapping. By further position alignment, we focused on six chromosome intervals, which preliminarily revealed the genetic regulatory mechanism of the internode on plant height. KASP markers suitable for high-throughput analysis were developed based on the low-homology markers in the six chromosome intervals, and they were further validated in 101 wheat accessions. The polymerization ofandhad high selection efficiency which might be decreased if further intruding. It suggests that the selection ofandshould mainly focus on the alleles reducing plant height, andcould be used in marker-assisted selection for internode length below spike (D1). The results in this study may provide assistance for wheat height genetic improvement in the middle and lower reaches of the Yangtze River.

    wheat (L.); Ningmai 9; Yangmai 158; plant height; KASP marker

    10.3724/SP.J.1006.2020.91063

    本研究由國家重點研發(fā)計劃項目(2017YFD0100801), 國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(CARS-3)和國家自然科學(xué)基金項目(31671690)資助。

    This study was supported by the National Key Project for the Research and Development of China (2017YFD0100801), the China Agriculture Research System (CARS-3), and the National Natural Science Foundation of China (31671690).

    馬鴻翔, E-mail: hxma@jaas.ac.cn

    E-mail: hmjp2005@163.com

    2019-10-15;

    2020-01-15;

    2020-02-17.

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

    av国产免费在线观看| 99国产精品99久久久久| 麻豆国产av国片精品| 亚洲国产欧洲综合997久久,| 2021天堂中文幕一二区在线观| 亚洲精品在线美女| 天天一区二区日本电影三级| 国产av麻豆久久久久久久| 亚洲人成网站高清观看| 欧美丝袜亚洲另类 | 俺也久久电影网| 免费看美女性在线毛片视频| 88av欧美| 少妇被粗大的猛进出69影院| 亚洲国产精品999在线| 免费在线观看影片大全网站| 无限看片的www在线观看| 老司机深夜福利视频在线观看| 色尼玛亚洲综合影院| 国产成人影院久久av| 黄色视频不卡| 免费电影在线观看免费观看| 国产男靠女视频免费网站| 国内少妇人妻偷人精品xxx网站 | 看片在线看免费视频| 国产片内射在线| 成人三级做爰电影| 中文字幕精品亚洲无线码一区| 男人的好看免费观看在线视频 | 国产成人av激情在线播放| 天堂影院成人在线观看| 亚洲av电影不卡..在线观看| 午夜激情福利司机影院| 精品熟女少妇八av免费久了| 国产亚洲av嫩草精品影院| 欧美成人一区二区免费高清观看 | 男女之事视频高清在线观看| 91成年电影在线观看| 深夜精品福利| 欧美在线黄色| 女人高潮潮喷娇喘18禁视频| 一进一出抽搐gif免费好疼| 成年版毛片免费区| 日韩欧美在线乱码| 国产aⅴ精品一区二区三区波| 免费在线观看黄色视频的| 禁无遮挡网站| 久久精品91无色码中文字幕| 最近最新免费中文字幕在线| 亚洲一码二码三码区别大吗| 亚洲国产欧美人成| av国产免费在线观看| 波多野结衣高清无吗| 精品国产亚洲在线| 精品少妇一区二区三区视频日本电影| 国内精品一区二区在线观看| 波多野结衣巨乳人妻| 久久国产精品人妻蜜桃| 国内久久婷婷六月综合欲色啪| 法律面前人人平等表现在哪些方面| 成人亚洲精品av一区二区| 可以免费在线观看a视频的电影网站| 免费在线观看黄色视频的| 久久九九热精品免费| 午夜精品在线福利| 99riav亚洲国产免费| 亚洲精品久久国产高清桃花| 亚洲五月天丁香| 岛国在线免费视频观看| 亚洲av第一区精品v没综合| 亚洲av电影在线进入| 天天添夜夜摸| 叶爱在线成人免费视频播放| 美女黄网站色视频| 欧美av亚洲av综合av国产av| 我的老师免费观看完整版| 亚洲专区国产一区二区| 小说图片视频综合网站| 国产精品99久久99久久久不卡| www日本在线高清视频| 国产乱人伦免费视频| 一二三四在线观看免费中文在| 国产精品亚洲美女久久久| 国产午夜精品论理片| 看片在线看免费视频| 国产激情偷乱视频一区二区| 黄色毛片三级朝国网站| 国产激情偷乱视频一区二区| 夜夜躁狠狠躁天天躁| 精品高清国产在线一区| 欧美高清成人免费视频www| 亚洲精品av麻豆狂野| 波多野结衣巨乳人妻| 欧美乱码精品一区二区三区| 色综合婷婷激情| 男插女下体视频免费在线播放| 精品国产乱码久久久久久男人| 国产午夜福利久久久久久| 在线观看舔阴道视频| 国产成人aa在线观看| 黄片大片在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品av麻豆狂野| 搡老熟女国产l中国老女人| 欧美性长视频在线观看| 欧美精品亚洲一区二区| 久久久久久免费高清国产稀缺| av免费在线观看网站| 不卡一级毛片| 91字幕亚洲| 欧美日本亚洲视频在线播放| 日本免费一区二区三区高清不卡| av中文乱码字幕在线| 色噜噜av男人的天堂激情| 亚洲国产看品久久| 亚洲,欧美精品.| 五月伊人婷婷丁香| 亚洲精品在线观看二区| 成人三级黄色视频| 亚洲自拍偷在线| 国产蜜桃级精品一区二区三区| 欧美在线黄色| 精品国产美女av久久久久小说| 天堂动漫精品| 91成年电影在线观看| 亚洲七黄色美女视频| 嫁个100分男人电影在线观看| 一进一出抽搐动态| 成人国语在线视频| 老司机福利观看| 午夜久久久久精精品| 国产成人av激情在线播放| 韩国av一区二区三区四区| 国产亚洲精品综合一区在线观看 | 久久精品综合一区二区三区| 国产精品av久久久久免费| 中亚洲国语对白在线视频| 一本久久中文字幕| 操出白浆在线播放| 国产视频内射| 久久久久久亚洲精品国产蜜桃av| 国产v大片淫在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 特大巨黑吊av在线直播| 成人18禁在线播放| 午夜久久久久精精品| 一级毛片精品| 中文资源天堂在线| 欧美精品啪啪一区二区三区| 国产精品99久久99久久久不卡| 国产99白浆流出| 午夜影院日韩av| 欧美日本亚洲视频在线播放| 一a级毛片在线观看| 免费av毛片视频| 欧美黑人欧美精品刺激| 两个人免费观看高清视频| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美在线一区二区| 色综合站精品国产| 亚洲av电影在线进入| 丝袜人妻中文字幕| 国产熟女午夜一区二区三区| xxxwww97欧美| √禁漫天堂资源中文www| 男女那种视频在线观看| 麻豆一二三区av精品| 曰老女人黄片| 窝窝影院91人妻| 国产亚洲精品一区二区www| av福利片在线观看| 亚洲国产高清在线一区二区三| 动漫黄色视频在线观看| 97碰自拍视频| 久久亚洲精品不卡| 99在线人妻在线中文字幕| 久久中文字幕人妻熟女| а√天堂www在线а√下载| 香蕉国产在线看| 亚洲人成伊人成综合网2020| 欧美成人免费av一区二区三区| 国产aⅴ精品一区二区三区波| 日韩欧美在线乱码| 婷婷亚洲欧美| 女同久久另类99精品国产91| 一进一出抽搐gif免费好疼| 女人高潮潮喷娇喘18禁视频| 亚洲激情在线av| 九色成人免费人妻av| 少妇裸体淫交视频免费看高清 | 免费在线观看亚洲国产| 成人国语在线视频| 欧美久久黑人一区二区| 色综合婷婷激情| 老熟妇乱子伦视频在线观看| 悠悠久久av| 国产又色又爽无遮挡免费看| 女生性感内裤真人,穿戴方法视频| a级毛片a级免费在线| 欧美成人一区二区免费高清观看 | av超薄肉色丝袜交足视频| 亚洲五月天丁香| 午夜精品久久久久久毛片777| 精品一区二区三区四区五区乱码| 99久久综合精品五月天人人| 久久人妻av系列| 啦啦啦观看免费观看视频高清| 欧美乱色亚洲激情| 好男人电影高清在线观看| 在线观看一区二区三区| 两个人免费观看高清视频| 90打野战视频偷拍视频| 母亲3免费完整高清在线观看| 欧美成人一区二区免费高清观看 | 亚洲国产欧美人成| 日韩av在线大香蕉| 精品欧美国产一区二区三| 可以在线观看的亚洲视频| 欧美黄色淫秽网站| 欧美中文综合在线视频| 丁香六月欧美| 亚洲va日本ⅴa欧美va伊人久久| 在线观看www视频免费| 久久亚洲精品不卡| 精品免费久久久久久久清纯| 女警被强在线播放| 国产成人啪精品午夜网站| 曰老女人黄片| 国产男靠女视频免费网站| 国产精品av视频在线免费观看| 深夜精品福利| 欧美在线一区亚洲| xxxwww97欧美| 精品国内亚洲2022精品成人| 手机成人av网站| 老司机靠b影院| 国产免费av片在线观看野外av| 国内毛片毛片毛片毛片毛片| 18美女黄网站色大片免费观看| 黄色片一级片一级黄色片| 日本三级黄在线观看| 亚洲国产欧美网| 久久国产精品影院| 国产三级在线视频| 天堂√8在线中文| 欧美性猛交╳xxx乱大交人| 禁无遮挡网站| 蜜桃久久精品国产亚洲av| 日韩国内少妇激情av| 欧美一区二区精品小视频在线| 亚洲av成人不卡在线观看播放网| 久久天堂一区二区三区四区| 深夜精品福利| 丝袜人妻中文字幕| 特大巨黑吊av在线直播| 特级一级黄色大片| 又粗又爽又猛毛片免费看| 久久久水蜜桃国产精品网| 男人的好看免费观看在线视频 | 欧美又色又爽又黄视频| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 久久精品国产99精品国产亚洲性色| 欧美最黄视频在线播放免费| 久热爱精品视频在线9| 99re在线观看精品视频| 视频区欧美日本亚洲| 好男人在线观看高清免费视频| 精品久久久久久久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美+亚洲+日韩+国产| 黄色视频不卡| 禁无遮挡网站| 久久热在线av| 夜夜躁狠狠躁天天躁| 精品国产亚洲在线| 在线免费观看的www视频| 国产精品av久久久久免费| 亚洲成人久久性| 日韩免费av在线播放| 亚洲激情在线av| 国产亚洲欧美在线一区二区| 一个人观看的视频www高清免费观看 | 午夜福利视频1000在线观看| 中出人妻视频一区二区| 成人永久免费在线观看视频| 国产精品久久久久久亚洲av鲁大| 黄色视频不卡| 最近最新免费中文字幕在线| 久久久久精品国产欧美久久久| 制服人妻中文乱码| 黄色成人免费大全| 国产高清激情床上av| 国产伦人伦偷精品视频| 精品久久蜜臀av无| 国内精品一区二区在线观看| 亚洲美女视频黄频| 亚洲精华国产精华精| 成熟少妇高潮喷水视频| 日本在线视频免费播放| 怎么达到女性高潮| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线| 久久精品91无色码中文字幕| 国产亚洲欧美98| 国产亚洲精品一区二区www| 美女大奶头视频| 91av网站免费观看| 又爽又黄无遮挡网站| 天天一区二区日本电影三级| 国产精品香港三级国产av潘金莲| 欧美色视频一区免费| 欧美绝顶高潮抽搐喷水| 五月玫瑰六月丁香| 久久久久性生活片| 色老头精品视频在线观看| 国产av一区二区精品久久| 制服诱惑二区| 成人永久免费在线观看视频| 欧美在线黄色| 亚洲av第一区精品v没综合| 亚洲欧美精品综合久久99| www日本在线高清视频| 久久这里只有精品中国| 在线观看免费日韩欧美大片| 黑人操中国人逼视频| 97人妻精品一区二区三区麻豆| 亚洲av电影在线进入| 欧美色欧美亚洲另类二区| 国产精品亚洲美女久久久| 日韩成人在线观看一区二区三区| 制服丝袜大香蕉在线| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲美女久久久| 亚洲精品一区av在线观看| 变态另类丝袜制服| 成年人黄色毛片网站| 精品第一国产精品| 欧美成人午夜精品| 国产精品久久久久久精品电影| cao死你这个sao货| 高清在线国产一区| 老司机在亚洲福利影院| 一进一出好大好爽视频| 最近在线观看免费完整版| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 精品国产美女av久久久久小说| 久久午夜亚洲精品久久| 一级毛片女人18水好多| 午夜视频精品福利| 欧美日韩中文字幕国产精品一区二区三区| 亚洲色图 男人天堂 中文字幕| 国产成人av激情在线播放| 国产视频一区二区在线看| 十八禁网站免费在线| 十八禁人妻一区二区| 国产久久久一区二区三区| 亚洲欧美精品综合一区二区三区| www.999成人在线观看| 无人区码免费观看不卡| 国产伦一二天堂av在线观看| 日韩精品免费视频一区二区三区| 亚洲国产欧美一区二区综合| 精品第一国产精品| 国产三级中文精品| 亚洲精品中文字幕一二三四区| √禁漫天堂资源中文www| 国产99白浆流出| 国产精品免费一区二区三区在线| 欧美日韩福利视频一区二区| 日韩大尺度精品在线看网址| 97碰自拍视频| 国内毛片毛片毛片毛片毛片| 亚洲成a人片在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器 | 日韩欧美精品v在线| a级毛片在线看网站| 可以在线观看毛片的网站| 女人高潮潮喷娇喘18禁视频| 成人午夜高清在线视频| 老司机福利观看| av福利片在线观看| 久久精品综合一区二区三区| 精品熟女少妇八av免费久了| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av日韩精品久久久久久密| 在线看三级毛片| 国产精品99久久99久久久不卡| 国产欧美日韩精品亚洲av| 桃红色精品国产亚洲av| 精品高清国产在线一区| 日韩欧美国产一区二区入口| 黄频高清免费视频| 中文在线观看免费www的网站 | a在线观看视频网站| 9191精品国产免费久久| 99re在线观看精品视频| 久久久久国内视频| 99国产精品99久久久久| 色噜噜av男人的天堂激情| 青草久久国产| 欧美乱色亚洲激情| 免费电影在线观看免费观看| 国产69精品久久久久777片 | 久久亚洲精品不卡| 中国美女看黄片| netflix在线观看网站| 国产av一区二区精品久久| 国产日本99.免费观看| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 法律面前人人平等表现在哪些方面| 老司机靠b影院| 岛国在线免费视频观看| 久久久久久久久免费视频了| 国产成人系列免费观看| 欧美成狂野欧美在线观看| 黑人操中国人逼视频| 美女黄网站色视频| 首页视频小说图片口味搜索| 久久精品91蜜桃| 一夜夜www| 国产在线精品亚洲第一网站| 亚洲真实伦在线观看| 十八禁人妻一区二区| 中文资源天堂在线| 亚洲精品色激情综合| 国产成人欧美在线观看| 亚洲国产高清在线一区二区三| 叶爱在线成人免费视频播放| 成人18禁高潮啪啪吃奶动态图| 2021天堂中文幕一二区在线观| a在线观看视频网站| 亚洲自偷自拍图片 自拍| 日本精品一区二区三区蜜桃| av片东京热男人的天堂| 国产精品乱码一区二三区的特点| 亚洲美女视频黄频| 非洲黑人性xxxx精品又粗又长| 国语自产精品视频在线第100页| 男人舔女人的私密视频| 人妻久久中文字幕网| 亚洲欧美精品综合一区二区三区| 中文字幕熟女人妻在线| 日韩欧美国产在线观看| av天堂在线播放| 亚洲精华国产精华精| 精品高清国产在线一区| 国产一区二区在线观看日韩 | 精品国内亚洲2022精品成人| 午夜影院日韩av| 制服诱惑二区| 国产精品一区二区三区四区免费观看 | 亚洲全国av大片| 夜夜爽天天搞| 一区二区三区国产精品乱码| 免费av毛片视频| xxxwww97欧美| 亚洲熟女毛片儿| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频| 久久亚洲精品不卡| 首页视频小说图片口味搜索| 亚洲成人久久性| 日本 欧美在线| 久久精品国产综合久久久| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看| 国产爱豆传媒在线观看 | 久久久久国产精品人妻aⅴ院| 哪里可以看免费的av片| 亚洲精品国产一区二区精华液| 久久欧美精品欧美久久欧美| 两性夫妻黄色片| 激情在线观看视频在线高清| 一个人免费在线观看的高清视频| 正在播放国产对白刺激| 嫩草影视91久久| 免费看十八禁软件| 免费在线观看视频国产中文字幕亚洲| 波多野结衣高清作品| 亚洲国产精品久久男人天堂| 97人妻精品一区二区三区麻豆| 一本久久中文字幕| 久久久精品大字幕| 黄色成人免费大全| 好男人电影高清在线观看| 99久久精品热视频| 中文字幕熟女人妻在线| 十八禁人妻一区二区| 国产麻豆成人av免费视频| 精品久久久久久成人av| 在线观看免费视频日本深夜| 久久久久久大精品| 正在播放国产对白刺激| 男人舔女人下体高潮全视频| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| 麻豆av在线久日| 欧美精品啪啪一区二区三区| 久久久久久久午夜电影| 脱女人内裤的视频| 免费av毛片视频| 国产单亲对白刺激| 天天躁夜夜躁狠狠躁躁| 久久性视频一级片| 国产精品一区二区精品视频观看| 听说在线观看完整版免费高清| 丁香六月欧美| 深夜精品福利| 最近最新中文字幕大全电影3| 婷婷精品国产亚洲av在线| 男人的好看免费观看在线视频 | 国产一区二区在线av高清观看| 成人亚洲精品av一区二区| 亚洲美女黄片视频| 国产真实乱freesex| 精品不卡国产一区二区三区| 一夜夜www| 国产亚洲精品一区二区www| 老汉色av国产亚洲站长工具| 观看免费一级毛片| 他把我摸到了高潮在线观看| 中出人妻视频一区二区| 成人国产一区最新在线观看| 欧美成狂野欧美在线观看| 婷婷精品国产亚洲av| 日本a在线网址| 女人爽到高潮嗷嗷叫在线视频| 九九热线精品视视频播放| 亚洲成人免费电影在线观看| 色播亚洲综合网| 超碰成人久久| 91av网站免费观看| 亚洲国产欧美网| 一个人免费在线观看的高清视频| 日本三级黄在线观看| 美女大奶头视频| 久久久久国产一级毛片高清牌| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 亚洲中文字幕一区二区三区有码在线看 | 黄色 视频免费看| 久久精品91无色码中文字幕| 亚洲av五月六月丁香网| 亚洲av第一区精品v没综合| 欧美 亚洲 国产 日韩一| 国产精品自产拍在线观看55亚洲| 真人做人爱边吃奶动态| 九色成人免费人妻av| 国产精品电影一区二区三区| 哪里可以看免费的av片| 99热这里只有是精品50| 久久中文看片网| 天堂影院成人在线观看| 五月伊人婷婷丁香| 91麻豆av在线| 亚洲真实伦在线观看| 免费在线观看亚洲国产| 久久久久亚洲av毛片大全| 午夜精品在线福利| 国产精品久久久久久久电影 | 国产精品98久久久久久宅男小说| 一夜夜www| 国产精品乱码一区二三区的特点| 我要搜黄色片| 99国产精品一区二区三区| 亚洲欧美一区二区三区黑人| 欧美 亚洲 国产 日韩一| 亚洲av第一区精品v没综合| 欧美精品啪啪一区二区三区| 欧美成人性av电影在线观看| 熟妇人妻久久中文字幕3abv| 国产伦人伦偷精品视频| 欧美日韩国产亚洲二区| 叶爱在线成人免费视频播放| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 国产亚洲精品久久久久5区| 亚洲最大成人中文| 亚洲国产精品成人综合色| av免费在线观看网站| 99riav亚洲国产免费| 日本黄色视频三级网站网址| 在线观看免费日韩欧美大片| 不卡一级毛片| 国产欧美日韩一区二区精品| 国内少妇人妻偷人精品xxx网站 | bbb黄色大片| 日本五十路高清| 91麻豆av在线| 久久这里只有精品中国| 国产精品永久免费网站| 亚洲av美国av| 国产激情偷乱视频一区二区| 九九热线精品视视频播放| 欧美乱码精品一区二区三区| 性色av乱码一区二区三区2| av福利片在线观看| 男插女下体视频免费在线播放| 精品国产亚洲在线| 后天国语完整版免费观看| ponron亚洲| 欧美色视频一区免费| 欧美日韩黄片免| 国产欧美日韩一区二区精品| 午夜福利在线在线| 免费看日本二区| 黄色女人牲交| 高潮久久久久久久久久久不卡| 欧美一区二区国产精品久久精品 |