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

    利用剩余雜合體衍生的近等基因系精細(xì)定位水稻粒長微效QTL qGL1.1

    2020-03-31 14:07:08李盼盼朱玉君郭梁莊杰云樊葉楊
    中國水稻科學(xué) 2020年2期
    關(guān)鍵詞:粒長株系區(qū)間

    李盼盼 朱玉君 郭梁 莊杰云 樊葉楊,*

    利用剩余雜合體衍生的近等基因系精細(xì)定位水稻粒長微效QTL

    李盼盼1朱玉君1郭梁2莊杰云1樊葉楊1,*

    (1中國水稻研究所 水稻生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室/國家水稻改良中心, 杭州 310006;2袁隆平農(nóng)業(yè)高科技股份有限公司, 長沙 410001;*通信聯(lián)系人, E-mail: fanyeyang@caas.cn)

    【】本研究旨在對前期在水稻第1染色體長臂521.8 kb的區(qū)間內(nèi)定位到的進(jìn)行精細(xì)定位。從和所在區(qū)間分別呈雜合的2個BC2F9單株配組衍生的F4群體中,篩選到Wn28826-RM1231區(qū)間內(nèi)雜合片段呈梯系排列的3個單株,構(gòu)建了3套F5:6近等基因系。2017年種植于浙江杭州,考查千粒重、粒長和粒寬。利用SAS軟件的GLM程序進(jìn)行雙因素方差分析,對的效應(yīng)進(jìn)行了驗(yàn)證。在此基礎(chǔ)上,篩選出雜合片段更小且呈交迭排列的6個剩余雜合體,發(fā)展了6套F8:9近等基因系,2018年種植于海南陵水。對每套近等基因系中雙親基因型株系的表型差異進(jìn)行雙因素方差分析。在2個試驗(yàn)中對粒長和千粒重均呈極顯著差異,效應(yīng)方向一致且大小穩(wěn)定。密陽46等位基因能分別增加粒長0.027 mm和提高千粒重0.17 g,貢獻(xiàn)率分別達(dá)到27.12%和19.09%。鑒于在前后試驗(yàn)中對粒長影響最為顯著,而對粒寬作用不顯著,故將重新命名為。通過比較各套近等基因系的分離區(qū)間的基因組位置,最終將定位于Wn29077和Wn29154之間約76.8 kb的區(qū)間內(nèi)。

    水稻;近等基因系;微效作用;粒長;精細(xì)定位

    水稻(L.)是保障糧食安全的重要作物之一,其單產(chǎn)的提高對糧食安全發(fā)揮了重要的作用。千粒重是水稻產(chǎn)量構(gòu)成因子的三要素之一,主要由籽粒大小和充實(shí)度控制,而籽粒大小一般以粒長、粒寬和粒厚等來衡量,此外粒長和粒寬還與稻米外觀品質(zhì)息息相關(guān)[1]。因此,分離和克隆籽粒大小相關(guān)的基因可為水稻高產(chǎn)優(yōu)質(zhì)育種奠定重要的理論基礎(chǔ)并提供可用的基因資源。

    籽粒大小屬于典型的數(shù)量性狀,由大量的數(shù)量性狀座位(Quantitative Trait Loci,QTL)控制。迄今已有17個QTL獲得克隆,它們分布于除水稻第1、10、11和12染色體外的8條染色體上。其中,/、、/、、/、//、、、和等10個QTL主要控制粒長和千粒重[2-16];、、和等4個QTL主要控制粒寬和千粒重[17-21];對粒長和粒寬的效應(yīng)大小相近且方向相同,故最終控制千粒重[22];/和雖然也對粒長和粒寬兼具效應(yīng),但其效應(yīng)方向相反,故對千粒重并不表現(xiàn)出顯著效應(yīng)[23-25]。這些QTL主要涉及植物激素、泛素-蛋白酶體通路和G蛋白信號等分子途徑[26-27]。這些研究豐富了我們對水稻籽粒大小的分子調(diào)控機(jī)制的了解,但其調(diào)控通路網(wǎng)絡(luò)尚不清楚。

    主效QTL和微效QTL均是解析復(fù)雜性狀遺傳機(jī)制的重要組成部分[28-29]。主效QTL因其效應(yīng)較大,能在不同遺傳背景和環(huán)境中穩(wěn)定檢測到,故而是精細(xì)定位和克隆的首選;相對地,微效QTL效應(yīng)較小,更易受遺傳背景和環(huán)境的干擾,故而在初級定位中較難穩(wěn)定檢測到[30],目前尚未有籽粒大小微效QTL克隆的報道。因此,籽粒大小微效QTL的發(fā)掘、精細(xì)定位和克隆,對全面揭示籽粒大小調(diào)控機(jī)制具有重要意義。

    在前期研究中,我們應(yīng)用珍汕97///珍汕97//珍汕97/密陽46衍生的高代回交群體,在水稻第1染色體長臂將千粒重QTL分解為和,其中,被界定于RM11554-RM11569約521.8 kb的區(qū)間內(nèi),密陽46等位基因顯著增加粒長,提高粒重[31]。本研究針對,新構(gòu)建了分離區(qū)間呈梯系排列的3套近等基因系(Near Isogenic Lines, NILs),對其效應(yīng)進(jìn)行驗(yàn)證,并進(jìn)一步發(fā)展了分離區(qū)間更為精細(xì)且交迭排列的6套近等基因系,開展精細(xì)定位。

    1 材料與方法

    1.1 水稻材料

    本研究所用定位群體為9套近等基因系,其中3套F5:6近等基因系用于效應(yīng)驗(yàn)證,6套F8:9近等基因系用于精細(xì)定位,材料構(gòu)建過程(圖1)如下。

    利用珍汕97///珍汕97//珍汕97/密陽46組合衍生的1個BC2F9群體,從中挑選在區(qū)間(RM11437-RM11615)和區(qū)間(RM11615-RM11800)[32]分別呈雜合的2個單株,配組衍生1個F4群體。經(jīng)分離區(qū)間標(biāo)記檢測,篩選到Wn28826-RM1231區(qū)間內(nèi)雜合片段呈梯系排列的3個單株,自交形成3個NIL-F2型的F5群體,分別從每個群體中挑選分離區(qū)間呈母本純合型(即珍汕97純合型)和父本純合型(即密陽46純合型)的單株,自交后構(gòu)建了3套F5:6近等基因系,命名為Y1、Y2和Y3。Y1包含珍汕97純合型株系30個和密陽46純合型株系29個,Y2和Y3分別包含珍汕97純合型和密陽46純合型株系各30個。3套近等基因系的分離區(qū)間分別為Wn28826-RM1231、Wn28990-RM1231和Wn29154-RM1231(圖2-A)。

    Fig. 1. Procedure for developing the nine sets of near isogenic lines (NILs).

    同時,為進(jìn)一步篩選雜合區(qū)間更為精細(xì)的梯系剩余雜合體材料,從衍生出Y1的F5群體中挑選整個分離區(qū)間(Wn28826-RM1231)呈雜合的2個單株,經(jīng)過“剩余雜合體-NIL-F2型群體-剩余雜合體-NIL-F2型群體”2次標(biāo)記篩選,獲得了雜合片段更小且呈交疊排列的6個F7單株,自交形成6個NIL-F2型F8群體。從中挑選珍汕97純合型和密陽46純合型單株,自交后構(gòu)建了6套F8:9近等基因系,分別命名為LP1、LP2、LP3、LP4、LP5和LP6,每套近等基因系分別包含珍汕97純合型株系和密陽46純合型株系各30個。LP3和LP5來源于同一個F5單株,LP1、LP2、LP4和LP6則來源于另一個F5單株。6套近等基因系的分離區(qū)間分別為Wn28826-Wn28893、Wn28826-Wn28990、Wn28990-Wn29048、Wn29048-Wn29125、Wn29125-RM1231和Wn29154-RM1231(圖2-B)。

    1.2 田間試驗(yàn)及性狀考查

    3套F5:6近等基因系于2017年5月種植在浙江省杭州市中國水稻研究所試驗(yàn)基地(以下簡稱杭州試驗(yàn))。6套F8:9近等基因系于2018年12月種植在海南省陵水市中國水稻研究所試驗(yàn)基地(簡稱陵水試驗(yàn))。采用完全隨機(jī)區(qū)組設(shè)計,2次重復(fù),每株系種植1行,每行8株,行株距16.7 cm×26.7 cm,正常田間管理。每株系混收中間4株。根據(jù)Zhang等[31]方法挑選飽谷600粒,采用萬深SC-G自動種子考種分析及千粒重儀測定千粒重、粒長和粒寬。

    1.3 標(biāo)記檢測

    采用Zheng等[33]的方法提取微量DNA,其中NIL-F2型群體分單株提取DNA,近等基因系每株系取中間6株的葉片混提DNA。在RM11554-RM11569區(qū)間共開發(fā)了8對特異分子標(biāo)記,包括2對InDel、2對CAPS和4對dCAPS標(biāo)記(表1)。按常規(guī)方法進(jìn)行PCR擴(kuò)增[34],采用2%瓊脂糖凝膠電泳或6%~10%的非變性聚丙烯酰胺凝膠電泳檢測PCR擴(kuò)增產(chǎn)物和酶切產(chǎn)物。

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

    采用Microsoft Excel計算2個試驗(yàn)9套近等基因系千粒重、粒長和粒寬的基本統(tǒng)計數(shù)據(jù),包括平均值、標(biāo)準(zhǔn)差、變異系數(shù)、變異范圍、偏斜度和峰度。應(yīng)用SAS軟件的一般線性模型(Proc GLM)[35]對每套近等基因系中珍汕97純合型株系和密陽46純合型株系的表型差異進(jìn)行雙因素方差分析。當(dāng)檢測到顯著的(<0.05)基因型差異時,則計算加性效應(yīng)和貢獻(xiàn)率[36]。

    A–用于驗(yàn)證的3套F5:6近等基因系; B–用于精細(xì)定位的6套F8:9近等基因系。

    Fig. 2. Genotypic compositions of the nine NILs sets in the target region.

    2 結(jié)果與分析

    2.1 性狀分布與表型變異

    在千粒重、粒長和粒寬這3個性狀上,3套F5:6近等基因系在杭州試驗(yàn)以及6套F8:9近等基因系在陵水試驗(yàn)的性狀均呈連續(xù)分布,且偏斜度和峰度較小,表現(xiàn)出數(shù)量性狀的特點(diǎn)(表2)。

    從平均值以及變異范圍來看,各套近等基因系的粒長均值和變異范圍在兩地試驗(yàn)中相差不大,但各套近等基因系的千粒重和粒寬均值和變異范圍在杭州試驗(yàn)中均比陵水試驗(yàn)中小,表明粒長受地域環(huán)境的影響比千粒重和粒寬小。

    比較同一世代中各套近等基因系在同一個性狀上的平均值,發(fā)現(xiàn)杭州試驗(yàn)中Y1和Y3同一性狀的平均值較為接近,而Y2均值最大;陵水試驗(yàn)中6套近等基因系的平均值呈線性連續(xù)分布,且LP4在3個性狀上的平均值都高于其他5套近等基因系。表明同一世代各套近等基因系的遺傳背景還存在控制這3個性狀的微效QTL差異。

    表1 本研究新開發(fā)的DNA標(biāo)記

    表2 9套近等基因系千粒重、粒長和粒寬的表型變異

    按2種基因型劃分,比較各套近等基因系同一性狀雙親基因型分布,發(fā)現(xiàn)千粒重和粒長存在明顯的群體間差異(圖3)。Y1和LP4中2種基因型株系在千粒重和粒長性狀上,表現(xiàn)出珍汕97純合型株系趨向于低值區(qū)而密陽46純合型株系趨向于高值區(qū)的特點(diǎn),特別是粒長,呈現(xiàn)出孟德爾式分離。表明Y1和LP4的分離區(qū)間可能存在控制粒長和千粒重的QTL。

    圖3 9套近等基因系千粒重和粒長分布

    Fig. 3. Distributions of 1000-grain weight and grain length in nine NILs sets.

    表3 9套近等基因系千粒重、粒長和粒寬的QTL檢測結(jié)果

    ?密陽46等位基因取代珍汕97等位基因所產(chǎn)生的遺傳效應(yīng)。R?效應(yīng)對表型方差的貢獻(xiàn)率。

    ,Additive effect of replacing a Zhenshan 97 allele with a Milyang 46 allele.R,Proportion of phenotypic variance explained by the QTL effect.

    2.2 QTL效應(yīng)驗(yàn)證

    前期將定位于RM11554-RM11569約521.8 kb的區(qū)間,對千粒重和粒長表現(xiàn)顯著效應(yīng),其加性效應(yīng)大小分別為0.06 g和0.017 mm,增效等位基因均來自于密陽46[31]。杭州試驗(yàn)中3套近等基因系千粒重、粒長和粒寬的雙因素方差分析表明。在分離區(qū)間為Wn28826-RM1231的Y1中,3個性狀均在雙親純合基因型株系間檢測到極顯著差異(<0.01),加性效應(yīng)分別為0.13 g、0.023 mm和0.008 mm,增效等位基因均來自密陽46;在分離區(qū)間為Wn28990-RM1231的Y2中,僅粒長在雙親純合基因型株系間分別檢測到顯著差異(<0.05),加性效應(yīng)大小為0.015 mm;而在分離區(qū)間為Wn29154-RM1231的Y3中,3個性狀在雙親純合基因型株系間的差異均未達(dá)到顯著水平(表3)。比較本研究中QTL對粒長的效應(yīng)大小和方向,發(fā)現(xiàn)其與效應(yīng)大小相近,方向一致,表明Y1和Y2很可能存在的分離,而Y3則不存在的分離。比較3套近等基因系的分離區(qū)間(圖3-A),將初步界定于Wn28944-Wn29154之間約210.2 kb的范圍內(nèi)。

    2.3 精細(xì)定位

    Wn29048-Wn29125區(qū)間呈分離的LP4中,千粒重和粒長雙親純合基因型株系間均檢測到極顯著差異(<0.0001)。加性效應(yīng)大小分別為0.17 g和0.027 mm,增效等位基因均來自密陽46(圖4),貢獻(xiàn)率達(dá)到19.09%和27.12%,粒寬則未在雙親純合基因型株系間檢測到顯著差異,表明LP4存在的分離(表3)。

    圖4 攜帶ZS97純合型和MY46純合型qGL1.1等位基因的近等基因系的籽粒(標(biāo)尺:10 mm)

    Fig. 4. Grains of NILZS97and NILMY46for(Scale bar: 10 mm).

    分離區(qū)間為Wn28990-Wn29048的LP3中也檢測到千粒重和粒長的顯著差異,但值僅為0.0394,其加性效應(yīng)大小分別為0.08 g和0.012 mm,但其增效等位基因來自珍汕97,貢獻(xiàn)率為4.30%和5.25%。鑒于LP3中檢測到的千粒重和粒長QTL效應(yīng)大小僅為效應(yīng)值的一半左右,更重要的是效應(yīng)方向相反。因此,該QTL很可能不是,而是一個新的控制千粒重和粒長的QTL。

    分離區(qū)間為Wn29154-RM1231的LP6中僅檢測到粒長的顯著差異(< 0.05),其加性效應(yīng)大小和方向與LP3類似:加性效應(yīng)大小為0.013 mm,珍汕97增效等位基因增加粒長,貢獻(xiàn)率為6.03%。鑒于LP3與LP6的分離區(qū)間相互獨(dú)立,我們推測LP6可能還存在另外一個控制粒長的QTL。

    而分離區(qū)間為Wn28826-Wn28893、Wn28826-Wn28990和Wn29125-RM1231的LP1、LP2和LP5中則均未檢測到3個性狀在雙親純合基因型株系間的顯著差異。有意思的是,LP5是除LP4外唯一1套粒長性狀上密陽46純合型均值大于珍汕97純合型均值的近等基因系,一個密陽46等位基因取代珍汕97等位基因所產(chǎn)生的粒長遺傳效應(yīng)為(8.105–8.083) / 2 = 0.011 mm,該值正好約等于LP4和LP6的粒長加性效應(yīng)值之和。因此我們推測,LP5粒長的表現(xiàn)很可能是和LP6所檢測到的QTL共同作用的結(jié)果。

    綜上所述,最終將定位于LP4和LP5兩套近等基因系共有分離標(biāo)記Wn29125及其兩側(cè)重組區(qū)間,即Wn29077-Wn29154,大小約76.8 kb。鑒于在前期試驗(yàn)和本研究中對粒長表現(xiàn)最為顯著,且對粒寬作用不顯著,故將重新命名為。

    3 討論

    微效QTL由于效應(yīng)較小,很容易被遺傳背景和環(huán)境的效應(yīng)掩蓋,所以較難在初定位中穩(wěn)定檢測到。而構(gòu)建遺傳背景高度純合的近等基因系,可以有效消除遺傳背景的干擾,進(jìn)而對效應(yīng)較小的微效QTL 進(jìn)行驗(yàn)證和分解, 乃至精細(xì)定位和克隆[37]。近年來,本實(shí)驗(yàn)室利用剩余雜合體構(gòu)建近等基因系,通過多輪“梯系剩余雜合體篩選-近等基因系構(gòu)建- QTL分析”,在水稻第1染色體長臂定位和分解出6個微效QTL,其中、、和,分別界定在120.4 kb、77.5 kb、35.2 Mb和125.5 kb的染色體區(qū)間內(nèi)[31,38-40]。本研究將()定位至76.8 kb的范圍內(nèi),進(jìn)一步佐證了該方法和策略在微效QTL精細(xì)定位中的有效性和高效性。

    在F5:6近等基因系Y1中,粒長和千粒重的加性效應(yīng)分別為0.023 mm和0.13 g,值和貢獻(xiàn)率分別為0.0018和0.0017、8.13%和8.09%;而在F8:9近等基因系LP4中,粒長和千粒重的加性效應(yīng)分別為0.027 mm和0.17 g,值均小于0.0001,貢獻(xiàn)率分別為27.12%和19.09%。這表明隨著自交的連續(xù),遺傳背景更加同質(zhì),在加性效應(yīng)基本不變的情況下,顯著性水平和貢獻(xiàn)率均得到了大幅提升。此外,不同世代的近等基因系雖然分別種植于浙江杭州和海南陵水,但表現(xiàn)一致,表明在不同環(huán)境中能穩(wěn)定發(fā)揮作用,在高產(chǎn)育種中具有重要的應(yīng)用潛力和價值。

    LP5中雙親純合基因型株系間盡管沒有檢測到粒長的顯著差異,但從粒長的加性效應(yīng)大小來看,恰好近似于LP4和LP6的粒長加性效應(yīng)值之和,效應(yīng)的方向也與效應(yīng)較大的LP4中檢測到的QTL,即的效應(yīng)方向一致?;诖?,我們推測LP5同時包含了以及LP6中檢測到的另一個微效粒長QTL。因此,我們將定位于LP4和LP5共有分離標(biāo)記Wn29125及其兩側(cè)重組區(qū)間,即Wn29077-Wn29154區(qū)間內(nèi)。

    經(jīng)Rice Genome Annotation Project Database (http://rice.plantbiology.msu.edu/)檢索,在Wn29077-Wn29154區(qū)間76.8 kb的范圍內(nèi)預(yù)測有13個注釋基因。其中5個編碼已知功能域蛋白,8個是編碼未知功能蛋白。5個注釋基因分別是LOC_Os01g 50680(編碼OsSub2-假定同源類枯草桿菌蛋白酶)、LOC_Os01g50690(編碼包含G蛋白-β亞基WD重復(fù)結(jié)構(gòu)域的蛋白),LOC_Os01g50700(編碼脫水家族蛋白),LOC_Os01g50720(MYB家族轉(zhuǎn)錄因子)和LOC_Os01g50750(編碼C3H4型結(jié)構(gòu)域的鋅指蛋白)。8個未知功能的注釋基因包括4個編碼逆轉(zhuǎn)座子:LOC_Os01g50640、LOC_Os01g50650、LOC_ Os01g50660和LOC_Os01g50670,4個編碼未知功能的表達(dá)蛋白或假設(shè)蛋白:LOC_Os01 g50730、LOC_Os01g50740、LOC_Os01g50630和LOC_ Os01g50710。在已克隆的籽粒大小相關(guān)QTL中,編碼的G 蛋白γ亞基功能缺失時,促進(jìn)了細(xì)胞的伸長,進(jìn)而產(chǎn)生長粒的表型,此外,G蛋白的α亞基和β亞基也分別報道與籽粒大小有關(guān)[27]。因此,編碼包含G蛋白-β亞基WD重復(fù)結(jié)構(gòu)域蛋白的LOC_Os01g50690可能是的候選基因。

    LP3和LP6也均檢測到粒長的顯著效應(yīng),與相比,不但效應(yīng)較小,而且方向相反。考慮到這兩套近等基因系的分離區(qū)間相互獨(dú)立,我們推測在的上下游還各存在1個與相斥的粒長QTL,其效應(yīng)尚需進(jìn)一步驗(yàn)證。

    [1] Xing Y Z, Zhang Q F. Genetic and molecular bases of rice yield[J]., 2010, 61: 421-442.

    [2] Duan P G, Ni S, Wang J M, Zhang B L, Xu R, Wang Y X, Chen H Q, Zhu X D, Li Y H.Regulation ofby OsmiR396 controls grain size and yield in rice[J]., 2016, 2: 15203.

    [3] Hu J, Wang Y X, Fang Y X, Zeng L J, Xu J, Yu H P, Shi Z Y, Pan J J, Zhang D, Kang S J, Zhu L, Dong G J, Guo L B, Zeng D L, Zhang G H, Xie L H, Xiong G S, Li J Y, Qian Q. A rare allele ofenhances grain size and grain yield in rice[J]., 2015, 8: 1455-1465.

    [4] Yu J P, Xiong H Y, Zhu X Y,Zhang H L, Li H H, Miao J L, Wang W S, Tang Z S, Zhang Z Y, Yao G X, Zhang Q, Pan Y H, Wang X, Rashid M A R, Li J J, Gao Y M, Li Z K, Yang W C, Fu X D, Li Z C.contributing to rice grain length and yield was mined by Ho-LAMap[J]., 2017, 15: 28.

    [5] Yu J P, Miao J L, Zhang Z Y, Xiong H Y, Zhu X Y, Sun X M, Pan Y H, Liang Y T, Zhang Q, Abdul Rehman R M, Li J J, Zhang H L, Li Z C. Alternative splicing ofcontrols grain length and yield inrice[J]., 2018, 16: 1667-1678.

    [6] Liu Q, Han R X, Wu K, Zhang J Q, Ye Y F, Wang S S, Chen J F, Pan Y J, Li Q, Xu X P, Zhou J W, Tao D Y, Wu Y J, Fu X D. G-protein βγ subunits determine grain size through interaction with MADS-domain transcription factors in rice[J]., 2018, 9: 852.

    [7] Fan C C, Xing Y Z, Mao H L, Lu T T, Han B, Xu C G, Li X H, Zhang Q F., a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein[J]., 2006, 112: 1164-1171.

    [8] Qi P, Lin Y S, Song X J, Shen J B, Huang W, Shan J X, Zhu M Z, Jiang L W, Gao P J, Lin H X. The novel quantitative trait locuscontrols rice grain size and yield by regulating Cyclin-T1;3[J]., 2012, 22: 1666-1680.

    [9] Zhang X J, Wang J F, Huang J, Lan H X, Wang C L, Yin C F, Wu Y Y, Tang H J, Qian Q, Li J Y, Zhang H S. Rare allele ofassociated with grain length causes extra-large grain and a significant yield increase in rice[J]., 2012, 109(52): 21534-21539.

    [10] Hu Z J, Lu S J, Wang M J, He H H, Sun L, Wang H R, Liu X H, Jiang L, Sun J L, Xin X Y, Kong W, Chu C C, Xue H W, Yang J S, Luo X J, Liu J X. A novel QTLencodes the GSK3/SHAGGY-like kinase OsGSK5/OsSK41 that interacts with OsARF4 to negatively regulate grain size and weight in rice[J]., 2018, 11:736-749.

    [11] Xia D, Zhou H, Liu R j, Dan W H, Li P B, Wu B, Chen J X, Wang L Q, Gao G J, Zhang Q L, He Y Q., a novel QTL encoding a GSK3/SHAGGY-like kinase, epistatically interacts withto produce extra-long grains in rice[J]., 2018, 11: 754-756.

    [12] Ying J Z, Ma M, Bai C, Huang X H, Liu J L, Fan Y Y, Song X J., a major QTL that negatively modulates grain length and weight in rice[J]., 2018, 11: 750-753.

    [13] Wu W G, Liu X Y, Wang M H, Meyer R S, Luo X J, Ndjiondjop M N, Tan L B, Zhang J W, Wu J Z, Cai H W, Sun C Q, Wang X K, Wing R A, Zhu Z F. A single-nucleotide polymorphism causes smaller grain size and loss of seed shattering during African rice domestication[J]., 2017, 3: 17064.

    [14] Wang A H, Hou Q Q, Si L Z,Huang X H, Luo J H, Lu D F, Zhu J J, Shangguan Y Y, Miao J S, Xie Y F, Wang Y C, Zhao Q, Feng Q, Zhou C C, Li Y, Fan D L, Lu Y Q, Tian Q L, Wang Z X, Han B.The PLATZ transcription factor GL6 affects grain length and number in rice[J]., 2019, 180 : 2077-2090.

    [15] Ishimaru K, Hirotsu N, Madoka Y, Murakami N, Hara N, Onodera H, Kashiwagi T, Ujiie K, Shimizu B I, Onishi A, Miyagawa H, Katoh E. Loss of function of the IAA-glucose hydrolase geneenhances rice grain weight and increases yield[J]., 2013, 45(6): 707-711.

    [16] Si L Z, Chen J Y, Huang X H, Gong H, Luo J H, Hou Q Q, Zhou T Y, Lu T T, Zhu J J, Shangguan Y Y, Chen E W, Gong C X, Zhao Q, Jing Y F, Zhao Y, Li Y, Cui L L, Fan D L, Lu Y Q, Weng Q J, Wang Y C, Zhan Q L, Liu K Y, Wei X H, An K, An G, Han B.controls grain size in cultivated rice[J]., 2016, 48(4): 447-456.

    [17] Song X J, Huang W, Shi M, Zhu M Z, Lin H X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase[J]., 2007, 39(5): 623-630.

    [18] Li Y B, Fan C C, Xing Y Z, Jiang Y H, Luo L J, Sun L, Shao D, Xu C J, Li X H, Xiao J H, He Y Q, Zhang Q F. Natural variation inplays an important role in regulating grain size and yield in rice[J].2011, 43(12): 1266-1269.

    [19] Duan P G, Xu J S, Zeng D L, Zhang B L, Geng M F, Zhang G Z, Huang K, Huang L J, Xu R, Ge S, Qian Q, Li Y H. Natural variation in the promoter ofcontributes to grain size diversity in rice[J]., 2017, 10: 685-694.

    [20] Liu J F, Chen J, Zheng X M, Wu F Q, Lin Q B, Heng Y Q, Tian P, Cheng Z J, Yu X W, Zhou K N, Zhang X, Guo X P, Wang J L, Wang H Y, Wan J M.acts in the brassinosteroid signaling pathway to regulate grain width and weight in rice[J].2017, 3: 17043.

    [21] Wang S K, Wu K, Yuan Q B, Liu X Y, Liu Z B, Lin X Y, Zeng R Z, Zhu H T, Dong G J, Qian Q, Zhang G Q, Fu X D. Control of grain size, shape and quality byin rice[J].2012, 44(8): 950-954.

    [22] Song X J, Kuroha T, Ayano M, Furuta T, Nagai K, Komeda N, Segami S, Miura K, Ogawa D, Kamura T, Suzuki T, Higashiyama T, Yamasaki M, Mori H, Inukai Y, Wu J Z, Kitano H, Sakakibara H, Jacobsen S E, Ashikari M. Rare allele of a previously unidentified histone H4 acetyltransferase enhances grain weight, yield, and plant biomass in rice[J]., 2015, 112(1): 76-81.

    [23] Wang Y X, Xiong G S, Hu J, Jiang L, Yu H, Xu J, Fang Y X, Zeng L J, Xu E B, Xu J, Ye W J, Meng X B, Liu R F, Chen H Q, Jing Y H, Wang Y H, Zhu X D, Li J Y, Qian Q. Copy number variation at thelocus contributes to grain size diversity in rice[J].2015, 47(8): 944-948.

    [24] Wang S K, Li S, Liu Q, Wu K, Zhang J Q, Wang S S, Wang Y, Chen X B, Zhang Y, Gao C X, Wang F, Huang H X, Fu X D. Theregulatory module determines grain shape and simultaneously improves rice yield and grain quality[J]., 2015, 47(8): 949-954.

    [25] Zhao D S, Li Q F, Zhang C Q, Zhang C, Yang Q Q, Pan L X, Ren X Y, Lu J, Gu M H, Liu Q Q.acts as a transcriptional activator to regulate rice grain shape and appearance quality[J]., 2018, 9: 1240.

    [26] Li N, Xu R, Duan P G, Li Y H. Control of grain size in rice[J]., 2018, 31: 237-251.

    [27] 劉喜,牟昌鈴,周春雷,程治軍,江玲,萬建民. 水稻粒型基因克隆和調(diào)控機(jī)制研究進(jìn)展[J]. 中國水稻科學(xué),2018, 32(1): 1-11.

    Liu X, Mou C L, Zhou C L, Cheng Z J, Jiang L, Wang J M. Research progress on cloning and regulation mechanism of rice grain shape genes[J]., 2018, 32(1): 1-11. (in Chinese with English abstract)

    [28] Mackay T F C, Stone E A, Ayroles J F. The genetics of quantitative traits: challenges and prospects[J].2009, 10: 565-577.

    [29] Kumar J, Gupta D S, Gupta S, Dubey S, Gupta P, Kumar S. Quantitative trait loci from identification to exploitation for crop improvement[J].2017, 36: 1187-1213.

    [30] Yamamoto T, Yonemaru J, Yano M. Towards the understanding of complex traits in rice: substantially or superficially[J]., 2009, 16: 141-154.

    [31] Zhang H W, Fan Y Y, Zhu Y J, Chen J Y, Yu S B, Zhuang J Y. Dissection of theregion into two tightly-linked minor QTLs having stable effects for grain weight in rice[J]., 2016, 17: 98.

    [32] Guo L, Wang K, Chen J Y, Huang D R, Fan Y Y, Zhuang J Y. Dissection of two quantitative trait loci for grain weight linked in repulsion on the long arm of chromosome 1 of rice (L.) [J]2013, 1: 70-76.

    [33] Zheng K L, Huang N, Bennett J, Khush G S. PCR-based marker-assisted selection in rice breeding//IRRI Discussion Paper Series No. 12. Los Banos: International Rice Research Institute, 1995.

    [34] Chen X, Temnykh S, Xu Y, Cho Y G, McCouch S R. Development of a microsatellite framework map providing genome-wide coverage in rice(L.) [J].1997, 95: 553-567.

    [35] SAS Institute Inc. SAS/STAT User’s Guide[M]. Cary, NC: SAS Institute, 1999.

    [36] Dai W M, Zhang K Q, Wu J R, Wang L, Duan B W, Zheng K L, Cai R, Zhuang J Y. Validating a segment on the short arm of chromosome 6 responsible for genetic variation in the hull silicon content and yield traits of rice[J].2008, 160: 317-324.

    [37] Abiola O, Angel J M, Avner P, Bachmanov A A, Belknap J K, Bennett B, Blankenhorn E P, Blizard D A, Bolivar V, Brockmann G A, Buck K J, Bureau J F, Casley W L, Chesler E J, Cheverud J M, Churchill G A, Cook M, Crabbe J C, Crusio W E, Darvasi A, Haan G D, Dermant P, Doerge R W, Elliot R W, Farber C R, Flaherty L, Flint J, Gershenfeld H, Gibson J P, Gu J, Gu W, Himmelbauer H, Hitzemann R, Hsu H C, Hunter K, Iraqi F F, Jansen R C, Johnson T E, Jones B C, Kempermann G, Lammert F, Lu L, Manly K F, Matthews D B, Medrano J F, Mehrabian M, Mittlemann G, Mock B A, Mogil J S, Montagutelli X, Morahan G, Mountz J D, Nagase H, Nowakowski R S, O'Hara B F, Osadchuk A V, Paigen B, Palmer A A, Peirce J L, Pomp D, Rosemann M, Rosen G D, Schalkwyk L C, Seltzer Z, Settle S, Shimomura K, Shou S, Sikela J M, Siracusa L D, Spearow J L, Teuscher C, Threadgill D W, Toth L A, Toye A A, Vadasz C, Van Zant G, Wakeland E, Williams R W, Zhang H G, Zou F; Complex Trait Consortium. The nature and identification of quantitative trait loci: A community’s view[J].2003, 4(11): 911-916.

    [38] Wang L L, Chen Y Y, Guo L, Zhang H W, Fan Y Y, Zhuang J Y. Dissection ofto three QTLs for grain weight and grain size in rice (L.) [J]., 2015, 202: 119-127.

    [39] Dong Q, Zhang Z H, Wang L L,Zhu Y J, Fan Y Y, Mou T M, Ma L Y, Zhuang J Y. Dissection and fine-mapping of two QTL for grain size linked in a 460-kb region on chromosome 1 of rice[J]., 2018, 11: 44.

    [40] Wang W H, Wang L L, Zhu Y J, Fan Y Y, Zhuang J Y. Fine-mapping of, a quantitative trait locus for 1000-grain weight in rice[J]., 2019, 26(4): 220-228.

    Fine Mapping of, a Minor QTL for Grain Length, Using Near Isogenic Lines Derived from Residual Heterozygotes in Rice

    LI Panpan1, ZHU Yujun1, GUO Liang2, ZHUANG Jieyun1, FAN Yeyang1,*

    (1State Key Laboratory of Rice Biology and Chinese National Center for Rice Improvement, China National Rice Research Institute, Hangzhou 310006, China;2Yuan Longping High-tech Agriculture Co, Changsha 410001, China;*Corresponding author, E-mail: fanyeyang@caas.cn)

    The objective of this study is to fine-mapwhich was previously mapped in a 521.8-kb region on the long arm of rice chromosome 1.In this study, two BC2F9plants carryingand,respectively, was crossed and produced an F4population. Three plants carrying sequential heterozygous segments in the interval Wn28826-RM1231 were selected and selfed to develop three sets of near isogenic lines (NILs) comprising the two homozygous genotypes. The NILs in F5:6were planted in Hangzhou in 2017 and the 1000-grain weight, grain length and grain width were measured. The effect ofwas validated by two-way analysis of variance (ANOVA) using SAS program. Then six new sequential residual heterozygotes carrying smaller heterozygous segments overlapped in the interval Wn28826-RM1231 were identified. Six sets of pairwise NILs in F8:9were planted in Lingshui in 2018 and the three traits were measured. Phenotypic differences between the two homozygous genotypic groups were analyzed using ANOVA.The allelic direction ofremained consistent and the genetic effects were stable in the two trials. The Milyang 46 allele could significantly increase grain length and 1000-grain weight by 0.027 mm and 0.17 g with the contribution up to 27.12% and 19.09%, respectively.Considering thatstably showed significant effect to grain length but not to grain width in both previous and this study,was renamed asAs a result,was delimited into a 76.8-kb region flanked by Wn29077 and Wn29154 based on comparing of the segregating regions among the six sets of NILs.

    rice; near isogenic line; minor effect; grain length; fine mapping

    Q943.2; S511.03

    A

    1001-7216(2020)02-0125-10

    10.16819/j.1001-7216.2020.9125

    2019-11-21;

    2019-12-28。

    國家重點(diǎn)研發(fā)計劃資助項(xiàng)目(2017YFD0100305);中央公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2017RG001-2)。

    猜你喜歡
    粒長株系區(qū)間
    解兩類含參數(shù)的復(fù)合不等式有解與恒成立問題
    你學(xué)會“區(qū)間測速”了嗎
    過表達(dá)NtMYB4a基因增強(qiáng)煙草抗旱能力
    水稻粒長遺傳及其功能基因研究進(jìn)展
    嫦娥5號返回式試驗(yàn)衛(wèi)星小麥育種材料研究進(jìn)展情況
    秈稻粒長與稻米品質(zhì)的相關(guān)性及其育種應(yīng)用
    日本晴/R1126水稻重組自交系群體粒形性狀QTL定位
    衢州椪柑變異株系—黃皮椪柑相關(guān)特性研究
    浙江柑橘(2016年1期)2016-03-11 20:12:31
    區(qū)間對象族的可鎮(zhèn)定性分析
    小麥粒型相關(guān)性狀的QTL定位分析
    一级黄片播放器| 日本黄色视频三级网站网址| 亚洲精品在线观看二区| 天堂√8在线中文| 久久久久国内视频| 欧美高清性xxxxhd video| 欧美绝顶高潮抽搐喷水| 国产av在哪里看| 国产高清视频在线观看网站| 在现免费观看毛片| 国产大屁股一区二区在线视频| 成人国产综合亚洲| 亚洲片人在线观看| 欧美激情在线99| 性欧美人与动物交配| 波多野结衣巨乳人妻| 国产一区二区三区视频了| www.色视频.com| 欧美zozozo另类| 性欧美人与动物交配| 日日摸夜夜添夜夜添av毛片 | 亚洲专区中文字幕在线| 亚洲av.av天堂| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲欧美98| 真人做人爱边吃奶动态| 波多野结衣高清无吗| 长腿黑丝高跟| 18禁黄网站禁片午夜丰满| 久久久精品大字幕| ponron亚洲| 波多野结衣高清作品| 日韩高清综合在线| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 日本精品一区二区三区蜜桃| 性色avwww在线观看| 级片在线观看| av天堂中文字幕网| 久久精品夜夜夜夜夜久久蜜豆| 观看美女的网站| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 亚洲成人精品中文字幕电影| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 看十八女毛片水多多多| 亚洲av.av天堂| 午夜老司机福利剧场| 日韩欧美国产在线观看| 精品国产亚洲在线| 久久精品国产清高在天天线| 18美女黄网站色大片免费观看| 亚洲久久久久久中文字幕| 欧美精品啪啪一区二区三区| 美女大奶头视频| 俄罗斯特黄特色一大片| 老司机福利观看| 伦理电影大哥的女人| 国产不卡一卡二| 久久6这里有精品| 99久国产av精品| 国产免费男女视频| 精品国产三级普通话版| 午夜免费激情av| 中文字幕人妻熟人妻熟丝袜美| 国产精品98久久久久久宅男小说| 久久久久久久久大av| 97热精品久久久久久| 国产精品永久免费网站| 动漫黄色视频在线观看| 久久久久久国产a免费观看| 国产成人av教育| 国产在线男女| 亚洲精品日韩av片在线观看| 精品久久久久久久末码| 亚洲在线观看片| 色哟哟哟哟哟哟| 亚洲18禁久久av| 简卡轻食公司| 亚洲一区二区三区不卡视频| 国产免费av片在线观看野外av| 亚洲五月天丁香| 国产精品永久免费网站| 欧美一区二区精品小视频在线| 亚洲av日韩精品久久久久久密| 热99re8久久精品国产| 麻豆一二三区av精品| 欧美+亚洲+日韩+国产| 国产三级黄色录像| 国产成人欧美在线观看| 亚洲va日本ⅴa欧美va伊人久久| 成人三级黄色视频| 亚洲精品久久国产高清桃花| 黄色丝袜av网址大全| 精品不卡国产一区二区三区| 国产成人福利小说| 欧美3d第一页| 国产真实乱freesex| 精品一区二区三区av网在线观看| 日韩免费av在线播放| 成人特级av手机在线观看| 一级黄片播放器| 欧美激情国产日韩精品一区| 精品一区二区免费观看| 久久久国产成人免费| 欧美最黄视频在线播放免费| 午夜影院日韩av| 淫秽高清视频在线观看| 99久国产av精品| 免费av不卡在线播放| 亚洲av.av天堂| 亚洲国产色片| 看黄色毛片网站| 自拍偷自拍亚洲精品老妇| 91在线精品国自产拍蜜月| 午夜福利成人在线免费观看| 两个人的视频大全免费| 国产高清有码在线观看视频| 十八禁国产超污无遮挡网站| 欧美又色又爽又黄视频| 国模一区二区三区四区视频| 天堂网av新在线| 亚洲色图av天堂| 全区人妻精品视频| 男人狂女人下面高潮的视频| 国产单亲对白刺激| 夜夜躁狠狠躁天天躁| 一个人看的www免费观看视频| 成人av一区二区三区在线看| 91字幕亚洲| 亚洲国产精品sss在线观看| 九色成人免费人妻av| 老鸭窝网址在线观看| 久久久国产成人精品二区| 免费人成在线观看视频色| 国产av麻豆久久久久久久| 好男人电影高清在线观看| 亚洲成人精品中文字幕电影| 床上黄色一级片| 午夜亚洲福利在线播放| 欧美日韩福利视频一区二区| 一级av片app| 亚洲av中文字字幕乱码综合| 国产毛片a区久久久久| 国产白丝娇喘喷水9色精品| 波多野结衣高清无吗| 久久人人爽人人爽人人片va | 色精品久久人妻99蜜桃| 美女 人体艺术 gogo| 级片在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲乱码一区二区免费版| 天堂√8在线中文| 久久6这里有精品| 日韩有码中文字幕| 精品久久久久久久久久久久久| av在线老鸭窝| av福利片在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲专区中文字幕在线| 国产成年人精品一区二区| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片午夜丰满| 亚洲人成网站高清观看| 国产熟女xx| 老熟妇仑乱视频hdxx| 日本精品一区二区三区蜜桃| 中亚洲国语对白在线视频| 欧美极品一区二区三区四区| 国内精品久久久久精免费| 国产精品国产高清国产av| 大型黄色视频在线免费观看| 婷婷六月久久综合丁香| 国产黄色小视频在线观看| 丝袜美腿在线中文| 成年人黄色毛片网站| 美女 人体艺术 gogo| 亚洲欧美精品综合久久99| 亚洲va日本ⅴa欧美va伊人久久| 我的女老师完整版在线观看| 国产精品亚洲av一区麻豆| 国产精品1区2区在线观看.| 波多野结衣高清作品| 成人鲁丝片一二三区免费| av天堂中文字幕网| 亚洲av一区综合| 高清在线国产一区| 国产三级在线视频| 性插视频无遮挡在线免费观看| 精品一区二区三区av网在线观看| www日本黄色视频网| 免费观看的影片在线观看| 18+在线观看网站| 草草在线视频免费看| 一级作爱视频免费观看| 日韩国内少妇激情av| 欧美黑人欧美精品刺激| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 国产精品精品国产色婷婷| 欧美国产日韩亚洲一区| 国产精品久久久久久亚洲av鲁大| 成年人黄色毛片网站| 久久久久久久久久成人| 久久精品国产清高在天天线| 国产精品亚洲av一区麻豆| 精品一区二区三区视频在线观看免费| 我要搜黄色片| 成人精品一区二区免费| 在线免费观看的www视频| xxxwww97欧美| 人妻夜夜爽99麻豆av| 一个人免费在线观看电影| 国产大屁股一区二区在线视频| 午夜免费成人在线视频| 麻豆成人午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 久久精品人妻少妇| 日本一二三区视频观看| 国产黄片美女视频| 国产精品精品国产色婷婷| 最近最新中文字幕大全电影3| 91久久精品电影网| 99riav亚洲国产免费| 亚洲午夜理论影院| 免费高清视频大片| 在线播放国产精品三级| 97人妻精品一区二区三区麻豆| www.999成人在线观看| 国产淫片久久久久久久久 | 国产精品av视频在线免费观看| 一级作爱视频免费观看| 一区二区三区激情视频| 欧美成人性av电影在线观看| 欧美潮喷喷水| 在线免费观看的www视频| 亚洲第一欧美日韩一区二区三区| 婷婷精品国产亚洲av在线| 久久久久国内视频| 网址你懂的国产日韩在线| 嫩草影院精品99| 国产成人啪精品午夜网站| 日本免费一区二区三区高清不卡| 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久午夜电影| 免费在线观看日本一区| 久久亚洲精品不卡| 欧美成人免费av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 性色avwww在线观看| 国产探花极品一区二区| 天堂√8在线中文| 亚洲国产精品久久男人天堂| 日韩成人在线观看一区二区三区| 亚洲欧美激情综合另类| 人妻夜夜爽99麻豆av| 一个人观看的视频www高清免费观看| 色综合婷婷激情| 亚洲av美国av| 校园春色视频在线观看| 亚洲自拍偷在线| 69人妻影院| 精品日产1卡2卡| 日韩欧美国产一区二区入口| 久久久久久久精品吃奶| www.熟女人妻精品国产| 18禁黄网站禁片免费观看直播| bbb黄色大片| 一级黄色大片毛片| 亚洲精华国产精华精| 亚洲,欧美,日韩| 国产一区二区亚洲精品在线观看| 男女那种视频在线观看| 宅男免费午夜| 中文字幕免费在线视频6| 国产三级中文精品| 国产欧美日韩一区二区精品| 国产视频一区二区在线看| 一个人免费在线观看电影| 麻豆一二三区av精品| 欧美黄色片欧美黄色片| 最近最新中文字幕大全电影3| 国产不卡一卡二| 97超级碰碰碰精品色视频在线观看| 成年女人看的毛片在线观看| 色播亚洲综合网| 欧美激情国产日韩精品一区| 超碰av人人做人人爽久久| 天堂av国产一区二区熟女人妻| 久久中文看片网| 免费黄网站久久成人精品 | 亚洲第一电影网av| 婷婷六月久久综合丁香| 亚洲精品乱码久久久v下载方式| 亚洲欧美日韩无卡精品| 国产视频内射| 国产三级在线视频| 夜夜夜夜夜久久久久| 精品人妻偷拍中文字幕| 亚洲成人久久爱视频| 国产午夜福利久久久久久| 99久久无色码亚洲精品果冻| 老熟妇乱子伦视频在线观看| 91在线精品国自产拍蜜月| 在线播放无遮挡| 女生性感内裤真人,穿戴方法视频| 久久久久久久亚洲中文字幕 | 人人妻人人看人人澡| 久久精品影院6| av在线老鸭窝| 亚洲自偷自拍三级| 午夜精品一区二区三区免费看| 亚洲综合色惰| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 国产v大片淫在线免费观看| 最好的美女福利视频网| 国产69精品久久久久777片| 日本免费一区二区三区高清不卡| 一进一出抽搐动态| 日韩免费av在线播放| 最近在线观看免费完整版| 亚洲电影在线观看av| 亚洲精华国产精华精| 成人国产综合亚洲| 欧美黑人巨大hd| 国内精品一区二区在线观看| 欧美在线一区亚洲| 国产精品国产高清国产av| 国产黄色小视频在线观看| 好看av亚洲va欧美ⅴa在| 国产精品爽爽va在线观看网站| 免费一级毛片在线播放高清视频| 色综合站精品国产| 91在线观看av| 制服丝袜大香蕉在线| 欧美日韩国产亚洲二区| 一个人免费在线观看电影| 91狼人影院| 听说在线观看完整版免费高清| 国产精品野战在线观看| 色哟哟·www| 可以在线观看毛片的网站| 美女高潮喷水抽搐中文字幕| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 啦啦啦观看免费观看视频高清| 欧美黑人欧美精品刺激| 精品国内亚洲2022精品成人| 国产一区二区亚洲精品在线观看| 18+在线观看网站| 90打野战视频偷拍视频| 99久久99久久久精品蜜桃| 国产亚洲欧美在线一区二区| 精品久久久久久久人妻蜜臀av| 亚洲内射少妇av| 99精品在免费线老司机午夜| 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 亚洲av第一区精品v没综合| 日韩国内少妇激情av| 国产伦人伦偷精品视频| 在线十欧美十亚洲十日本专区| 嫩草影院新地址| 欧美乱色亚洲激情| 69av精品久久久久久| 伦理电影大哥的女人| 一个人观看的视频www高清免费观看| 国产精品一区二区三区四区免费观看 | 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区| 久久久久久九九精品二区国产| 精品熟女少妇八av免费久了| 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 午夜激情福利司机影院| 国产真实伦视频高清在线观看 | 欧美日韩瑟瑟在线播放| 中亚洲国语对白在线视频| 三级男女做爰猛烈吃奶摸视频| 天天躁日日操中文字幕| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 在线观看免费视频日本深夜| 国产亚洲av嫩草精品影院| 色哟哟哟哟哟哟| 熟女人妻精品中文字幕| 波多野结衣高清无吗| 久久99热6这里只有精品| 国产午夜精品久久久久久一区二区三区 | 丰满的人妻完整版| 老女人水多毛片| 丰满人妻一区二区三区视频av| 精品久久久久久久久久免费视频| 久久人人爽人人爽人人片va | 男女床上黄色一级片免费看| 怎么达到女性高潮| 久久久久久久午夜电影| 欧美日本视频| 51午夜福利影视在线观看| 伦理电影大哥的女人| 国产精品久久视频播放| 精品一区二区三区人妻视频| 精品人妻熟女av久视频| 亚洲自拍偷在线| 色av中文字幕| 午夜福利在线观看吧| 成人特级黄色片久久久久久久| 亚洲欧美清纯卡通| 日本免费一区二区三区高清不卡| 夜夜夜夜夜久久久久| 18禁黄网站禁片免费观看直播| 亚洲男人的天堂狠狠| 久久精品久久久久久噜噜老黄 | 午夜两性在线视频| 精品欧美国产一区二区三| 国产高潮美女av| 国产成人aa在线观看| 最好的美女福利视频网| 蜜桃久久精品国产亚洲av| 国产91精品成人一区二区三区| 夜夜躁狠狠躁天天躁| 久久精品综合一区二区三区| 色5月婷婷丁香| 亚洲欧美日韩卡通动漫| 熟女电影av网| 别揉我奶头~嗯~啊~动态视频| 日韩欧美国产一区二区入口| 亚洲人成伊人成综合网2020| 3wmmmm亚洲av在线观看| 九色国产91popny在线| 欧美黑人欧美精品刺激| 国产亚洲精品综合一区在线观看| 久久6这里有精品| 小说图片视频综合网站| 国产人妻一区二区三区在| 男人舔奶头视频| 亚洲第一欧美日韩一区二区三区| 亚洲激情在线av| 亚洲三级黄色毛片| 国内毛片毛片毛片毛片毛片| 波野结衣二区三区在线| 一本一本综合久久| 自拍偷自拍亚洲精品老妇| 中文字幕免费在线视频6| 高清日韩中文字幕在线| 午夜视频国产福利| 搡老妇女老女人老熟妇| 国产私拍福利视频在线观看| 又爽又黄a免费视频| 天堂√8在线中文| 欧美日本视频| 精品久久久久久久末码| 国产爱豆传媒在线观看| 成人鲁丝片一二三区免费| 男女做爰动态图高潮gif福利片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 婷婷丁香在线五月| 精品乱码久久久久久99久播| 真人做人爱边吃奶动态| 亚洲欧美日韩卡通动漫| 国内精品美女久久久久久| 99在线视频只有这里精品首页| 此物有八面人人有两片| 欧美高清性xxxxhd video| 在线免费观看的www视频| 午夜视频国产福利| 欧美在线一区亚洲| 我的老师免费观看完整版| 亚洲成a人片在线一区二区| 一区福利在线观看| 亚洲av成人精品一区久久| 噜噜噜噜噜久久久久久91| 婷婷亚洲欧美| 精品一区二区三区视频在线| 少妇人妻一区二区三区视频| 精品人妻一区二区三区麻豆 | 日韩精品青青久久久久久| 又爽又黄无遮挡网站| 免费av毛片视频| 国产伦精品一区二区三区四那| 国产精品女同一区二区软件 | 亚洲,欧美精品.| 人人妻人人看人人澡| 日韩av在线大香蕉| 日韩精品中文字幕看吧| 狠狠狠狠99中文字幕| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 香蕉av资源在线| 久久精品国产亚洲av天美| 成人三级黄色视频| 亚洲成人中文字幕在线播放| 中国美女看黄片| 两个人的视频大全免费| 中出人妻视频一区二区| 国产免费av片在线观看野外av| 精品久久久久久久久久久久久| 国产主播在线观看一区二区| 免费人成视频x8x8入口观看| 亚洲一区高清亚洲精品| 一区二区三区高清视频在线| 久久九九热精品免费| 精品人妻偷拍中文字幕| 亚洲成a人片在线一区二区| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 丁香欧美五月| 特大巨黑吊av在线直播| 久久精品国产自在天天线| 久久精品国产亚洲av天美| 亚洲av美国av| 免费一级毛片在线播放高清视频| 国产av不卡久久| 亚洲在线观看片| 美女黄网站色视频| 国产成人a区在线观看| 美女黄网站色视频| 午夜日韩欧美国产| 天天躁日日操中文字幕| 欧美+亚洲+日韩+国产| 日韩欧美在线二视频| 国产国拍精品亚洲av在线观看| 搡老熟女国产l中国老女人| 最后的刺客免费高清国语| 成人国产一区最新在线观看| 亚洲人成网站在线播| 精品久久国产蜜桃| 一个人看的www免费观看视频| 一进一出抽搐动态| 精品一区二区免费观看| 91在线观看av| 禁无遮挡网站| 18禁黄网站禁片免费观看直播| 国内精品久久久久精免费| 亚洲七黄色美女视频| av在线蜜桃| 色综合欧美亚洲国产小说| 少妇的逼好多水| 欧美不卡视频在线免费观看| bbb黄色大片| 黄色女人牲交| 韩国av一区二区三区四区| 男人的好看免费观看在线视频| 精品人妻1区二区| 观看免费一级毛片| 51午夜福利影视在线观看| 亚洲成人久久性| 国产精品久久电影中文字幕| a级毛片a级免费在线| 国产精品98久久久久久宅男小说| 国内精品久久久久久久电影| 国产精品,欧美在线| 国产精品亚洲美女久久久| 久久午夜亚洲精品久久| 国产不卡一卡二| 精品久久久久久久人妻蜜臀av| 久久午夜福利片| 人人妻人人澡欧美一区二区| 亚洲在线自拍视频| 久久久久久久精品吃奶| 97碰自拍视频| 九色国产91popny在线| 51午夜福利影视在线观看| 一进一出抽搐gif免费好疼| 久久精品国产亚洲av天美| 我要看日韩黄色一级片| 国产极品精品免费视频能看的| 精品午夜福利视频在线观看一区| 国产一区二区亚洲精品在线观看| 久久国产精品影院| 美女高潮的动态| 哪里可以看免费的av片| 国产精品1区2区在线观看.| 两人在一起打扑克的视频| 精品午夜福利在线看| 成人高潮视频无遮挡免费网站| 国产成人aa在线观看| 亚洲激情在线av| 99国产精品一区二区三区| 久久久国产成人免费| 久久婷婷人人爽人人干人人爱| 亚洲无线观看免费| 久9热在线精品视频| 男女下面进入的视频免费午夜| 久久久久久久久中文| 免费观看的影片在线观看| 少妇被粗大猛烈的视频| 久久久久久久久中文| 免费观看的影片在线观看| 欧美bdsm另类| 亚洲中文日韩欧美视频| 毛片女人毛片| 日本 av在线| 欧美高清成人免费视频www| 午夜精品在线福利| 国产精品av视频在线免费观看| 色5月婷婷丁香| 成人av在线播放网站| 国产伦精品一区二区三区视频9| 亚洲成av人片在线播放无| 在现免费观看毛片| 99久久精品热视频| 久久精品国产亚洲av涩爱 | 国产一区二区三区视频了| 欧美高清性xxxxhd video| 欧美性感艳星| 久久久久久九九精品二区国产| 中文在线观看免费www的网站| 欧美成人性av电影在线观看|