• <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定位分析
    日本一区二区免费在线视频| 久久人妻福利社区极品人妻图片 | 亚洲欧美精品综合一区二区三区| 一级毛片电影观看| 精品国产一区二区久久| 久久精品国产亚洲av高清一级| 男女之事视频高清在线观看 | 别揉我奶头~嗯~啊~动态视频 | 满18在线观看网站| 久久天躁狠狠躁夜夜2o2o | 美女福利国产在线| 黄片播放在线免费| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品久久午夜乱码| 精品第一国产精品| 久久影院123| 久久热在线av| 青春草亚洲视频在线观看| 国产亚洲午夜精品一区二区久久| 亚洲综合色网址| 欧美日韩亚洲国产一区二区在线观看 | 大码成人一级视频| 熟女少妇亚洲综合色aaa.| 亚洲综合色网址| 精品少妇久久久久久888优播| 久久青草综合色| 亚洲国产最新在线播放| 国产精品久久久久成人av| 精品少妇内射三级| 曰老女人黄片| 日韩一区二区三区影片| 国产亚洲一区二区精品| 久久久久久免费高清国产稀缺| 亚洲国产最新在线播放| 又黄又粗又硬又大视频| 伦理电影免费视频| 亚洲欧洲日产国产| 亚洲情色 制服丝袜| 久久青草综合色| 亚洲欧美一区二区三区久久| 亚洲精品国产av成人精品| 大片免费播放器 马上看| 人人妻人人添人人爽欧美一区卜| 国产免费视频播放在线视频| 最近手机中文字幕大全| 久久久久久久精品精品| 免费在线观看日本一区| 精品一区二区三区av网在线观看 | 一边摸一边抽搐一进一出视频| 手机成人av网站| 婷婷成人精品国产| 涩涩av久久男人的天堂| 国产日韩欧美亚洲二区| av在线老鸭窝| 精品一品国产午夜福利视频| 午夜免费成人在线视频| 国产爽快片一区二区三区| 悠悠久久av| 久久精品国产a三级三级三级| 欧美av亚洲av综合av国产av| 韩国精品一区二区三区| 日本五十路高清| 日韩电影二区| 一本一本久久a久久精品综合妖精| av欧美777| 日本黄色日本黄色录像| 成人亚洲精品一区在线观看| 亚洲九九香蕉| 亚洲激情五月婷婷啪啪| 天堂中文最新版在线下载| 在线看a的网站| 久久性视频一级片| 国产成人免费观看mmmm| 女警被强在线播放| 国产亚洲午夜精品一区二区久久| 黑人猛操日本美女一级片| 成人影院久久| 日韩免费高清中文字幕av| 亚洲一区中文字幕在线| 老司机在亚洲福利影院| 在线 av 中文字幕| 啦啦啦啦在线视频资源| 啦啦啦视频在线资源免费观看| netflix在线观看网站| 婷婷色综合www| 黄色视频在线播放观看不卡| 建设人人有责人人尽责人人享有的| 在线亚洲精品国产二区图片欧美| 精品人妻1区二区| 免费在线观看完整版高清| 美女高潮到喷水免费观看| 国产成人精品久久久久久| 亚洲精品久久久久久婷婷小说| 别揉我奶头~嗯~啊~动态视频 | 少妇的丰满在线观看| 欧美少妇被猛烈插入视频| 日本av手机在线免费观看| 国产免费视频播放在线视频| 人妻人人澡人人爽人人| 国产在视频线精品| 亚洲三区欧美一区| 日本五十路高清| 亚洲国产精品999| 成年美女黄网站色视频大全免费| 久久中文字幕一级| 制服诱惑二区| 狂野欧美激情性xxxx| 国产激情久久老熟女| 亚洲精品国产色婷婷电影| 日日摸夜夜添夜夜爱| 超碰成人久久| 久久人人97超碰香蕉20202| 亚洲一区中文字幕在线| 久久国产精品影院| 欧美 日韩 精品 国产| 久久人人97超碰香蕉20202| 新久久久久国产一级毛片| 国产精品久久久人人做人人爽| 五月开心婷婷网| 亚洲五月婷婷丁香| 最新在线观看一区二区三区 | 国产人伦9x9x在线观看| 久久热在线av| 午夜福利影视在线免费观看| 女性被躁到高潮视频| 亚洲成人手机| 侵犯人妻中文字幕一二三四区| a 毛片基地| 亚洲国产精品国产精品| avwww免费| 日韩熟女老妇一区二区性免费视频| 国产亚洲精品久久久久5区| 男女之事视频高清在线观看 | 少妇人妻久久综合中文| 一本—道久久a久久精品蜜桃钙片| 日本av免费视频播放| 操出白浆在线播放| 中国国产av一级| 热99国产精品久久久久久7| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产欧美在线一区| 黄片播放在线免费| 欧美日韩精品网址| 亚洲伊人色综图| 91精品三级在线观看| 精品人妻在线不人妻| 国产成人精品在线电影| 日韩一本色道免费dvd| 日韩av在线免费看完整版不卡| 精品国产乱码久久久久久小说| 狂野欧美激情性xxxx| 国产av国产精品国产| 国产精品久久久久久精品电影小说| 久久久精品区二区三区| 成人国产一区最新在线观看 | 成人国语在线视频| 久久热在线av| 97人妻天天添夜夜摸| 国产亚洲一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品一二三| 国产福利在线免费观看视频| 成年人免费黄色播放视频| 午夜福利在线免费观看网站| 51午夜福利影视在线观看| 久9热在线精品视频| 大片电影免费在线观看免费| 99久久综合免费| 亚洲国产欧美日韩在线播放| 天天影视国产精品| 久久人妻福利社区极品人妻图片 | 精品少妇一区二区三区视频日本电影| 午夜影院在线不卡| 欧美日韩综合久久久久久| 免费高清在线观看视频在线观看| 国产精品人妻久久久影院| 久久久久久久精品精品| 亚洲欧美日韩另类电影网站| 精品人妻1区二区| 欧美激情极品国产一区二区三区| 国产精品国产三级专区第一集| 成人国产av品久久久| 欧美黑人精品巨大| 丝袜美腿诱惑在线| 国产一区二区在线观看av| 欧美精品av麻豆av| 好男人电影高清在线观看| 丝袜喷水一区| 大陆偷拍与自拍| 操出白浆在线播放| av在线播放精品| 每晚都被弄得嗷嗷叫到高潮| 精品人妻在线不人妻| 成人免费观看视频高清| 母亲3免费完整高清在线观看| 国产一区二区在线观看av| 脱女人内裤的视频| 悠悠久久av| 五月开心婷婷网| 日本wwww免费看| 亚洲欧美一区二区三区久久| 亚洲中文字幕日韩| 宅男免费午夜| 色精品久久人妻99蜜桃| 久久久久视频综合| 久久精品亚洲熟妇少妇任你| 亚洲少妇的诱惑av| 19禁男女啪啪无遮挡网站| 电影成人av| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av成人精品| 老司机靠b影院| 亚洲久久久国产精品| 男的添女的下面高潮视频| 菩萨蛮人人尽说江南好唐韦庄| 色视频在线一区二区三区| 又黄又粗又硬又大视频| 午夜激情av网站| 两个人免费观看高清视频| 久久99一区二区三区| bbb黄色大片| 欧美少妇被猛烈插入视频| 国产精品免费大片| 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 久久亚洲国产成人精品v| 亚洲国产欧美一区二区综合| 亚洲欧洲日产国产| 久久精品久久久久久久性| 欧美黑人精品巨大| 一级,二级,三级黄色视频| 国产精品三级大全| 日韩中文字幕欧美一区二区 | 国产精品免费大片| 亚洲欧美一区二区三区国产| 亚洲国产看品久久| 高清黄色对白视频在线免费看| 日韩大片免费观看网站| 国产免费一区二区三区四区乱码| 少妇粗大呻吟视频| 色网站视频免费| 777米奇影视久久| 成年人午夜在线观看视频| 亚洲成人免费av在线播放| 免费高清在线观看日韩| 视频区图区小说| 亚洲中文日韩欧美视频| 国产高清视频在线播放一区 | 久久久国产精品麻豆| 观看av在线不卡| 色播在线永久视频| 亚洲成人免费av在线播放| 大陆偷拍与自拍| 大片电影免费在线观看免费| 不卡av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片黄色毛片免费观看视频| 国产成人免费观看mmmm| 国产在视频线精品| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 国产精品欧美亚洲77777| 水蜜桃什么品种好| 亚洲国产精品999| 欧美日韩福利视频一区二区| 亚洲,一卡二卡三卡| 美女视频免费永久观看网站| 男女床上黄色一级片免费看| 男女边摸边吃奶| 伦理电影免费视频| 亚洲精品国产区一区二| 另类亚洲欧美激情| 国产一区二区三区av在线| 搡老乐熟女国产| 午夜精品国产一区二区电影| 日韩熟女老妇一区二区性免费视频| 日韩制服骚丝袜av| 色婷婷久久久亚洲欧美| 搡老岳熟女国产| 汤姆久久久久久久影院中文字幕| av天堂久久9| 久久人人97超碰香蕉20202| 色94色欧美一区二区| 久久这里只有精品19| 91九色精品人成在线观看| 国产一区有黄有色的免费视频| 国产亚洲精品久久久久5区| 国产一区亚洲一区在线观看| 国产老妇伦熟女老妇高清| 亚洲精品久久成人aⅴ小说| 亚洲av欧美aⅴ国产| 欧美日韩亚洲国产一区二区在线观看 | 日本一区二区免费在线视频| 男男h啪啪无遮挡| 精品一区二区三卡| 久久99精品国语久久久| 午夜免费鲁丝| videos熟女内射| 日韩一本色道免费dvd| 中文字幕最新亚洲高清| 韩国高清视频一区二区三区| 欧美黑人精品巨大| 日本av免费视频播放| 免费在线观看日本一区| 嫩草影视91久久| 男的添女的下面高潮视频| 成人影院久久| 伊人久久大香线蕉亚洲五| 免费女性裸体啪啪无遮挡网站| 国产男人的电影天堂91| 国产精品99久久99久久久不卡| 两个人免费观看高清视频| 丝袜人妻中文字幕| www.熟女人妻精品国产| 十分钟在线观看高清视频www| 亚洲色图综合在线观看| 亚洲图色成人| 亚洲一码二码三码区别大吗| 免费av中文字幕在线| 国产黄色视频一区二区在线观看| 黄片小视频在线播放| 别揉我奶头~嗯~啊~动态视频 | 男男h啪啪无遮挡| 国产在线观看jvid| 欧美在线一区亚洲| av天堂久久9| 91字幕亚洲| 老鸭窝网址在线观看| 色综合欧美亚洲国产小说| 9热在线视频观看99| 日韩 欧美 亚洲 中文字幕| 美女大奶头黄色视频| 男女高潮啪啪啪动态图| 亚洲精品在线美女| www.av在线官网国产| 日本一区二区免费在线视频| 99久久综合免费| 美女中出高潮动态图| 少妇裸体淫交视频免费看高清 | 韩国高清视频一区二区三区| 日本黄色日本黄色录像| 少妇人妻 视频| 国产在视频线精品| 午夜福利视频在线观看免费| 国产免费福利视频在线观看| 男女国产视频网站| 欧美+亚洲+日韩+国产| 欧美日韩综合久久久久久| 久久久精品国产亚洲av高清涩受| 黄色毛片三级朝国网站| 女人被躁到高潮嗷嗷叫费观| 亚洲成色77777| 国产激情久久老熟女| 亚洲欧美精品综合一区二区三区| 宅男免费午夜| 国产亚洲精品久久久久5区| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 在线天堂中文资源库| 国产麻豆69| 97精品久久久久久久久久精品| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 99久久综合免费| 国产91精品成人一区二区三区 | 搡老岳熟女国产| 亚洲精品自拍成人| 一本综合久久免费| 亚洲成av片中文字幕在线观看| 麻豆av在线久日| 亚洲国产欧美网| 熟女av电影| 中文字幕色久视频| 黄网站色视频无遮挡免费观看| 国产精品一区二区在线观看99| 亚洲免费av在线视频| 九色亚洲精品在线播放| 亚洲第一青青草原| 婷婷丁香在线五月| 亚洲欧美激情在线| 色网站视频免费| 久久久久精品国产欧美久久久 | 丁香六月欧美| 午夜福利视频在线观看免费| 大型av网站在线播放| 久久ye,这里只有精品| 老汉色∧v一级毛片| 国产精品九九99| 欧美另类一区| 人人妻,人人澡人人爽秒播 | 极品人妻少妇av视频| 国产精品久久久久久精品电影小说| 18禁国产床啪视频网站| 日本色播在线视频| 美国免费a级毛片| 一本一本久久a久久精品综合妖精| 中文精品一卡2卡3卡4更新| 亚洲欧美激情在线| 欧美激情极品国产一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 亚洲伊人色综图| 日本av免费视频播放| 曰老女人黄片| 国产福利在线免费观看视频| 黄色a级毛片大全视频| 久久精品国产a三级三级三级| 国产三级黄色录像| 美女扒开内裤让男人捅视频| 99热全是精品| 肉色欧美久久久久久久蜜桃| av天堂久久9| 一个人免费看片子| 国产熟女午夜一区二区三区| 国语对白做爰xxxⅹ性视频网站| 一区二区日韩欧美中文字幕| 男人操女人黄网站| 黄色视频在线播放观看不卡| 国产欧美亚洲国产| 欧美中文综合在线视频| 美女视频免费永久观看网站| 久久人妻熟女aⅴ| 亚洲精品av麻豆狂野| 精品国产一区二区三区四区第35| 如日韩欧美国产精品一区二区三区| 丰满少妇做爰视频| 99国产精品免费福利视频| 成人国语在线视频| 999久久久国产精品视频| 国产成人精品久久二区二区91| 午夜免费鲁丝| 精品人妻在线不人妻| 午夜老司机福利片| 久久久精品区二区三区| 国产熟女欧美一区二区| 久久ye,这里只有精品| 久久精品国产a三级三级三级| 男女边摸边吃奶| 亚洲自偷自拍图片 自拍| 热99久久久久精品小说推荐| 亚洲国产av新网站| 天堂俺去俺来也www色官网| 欧美激情高清一区二区三区| 久久综合国产亚洲精品| 在线观看人妻少妇| 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| av有码第一页| 欧美久久黑人一区二区| 天天影视国产精品| 黄色片一级片一级黄色片| 黄色视频在线播放观看不卡| 成人免费观看视频高清| 国产片内射在线| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 日韩av免费高清视频| 脱女人内裤的视频| 婷婷成人精品国产| 天天躁夜夜躁狠狠久久av| 亚洲成人手机| 精品少妇黑人巨大在线播放| 少妇的丰满在线观看| 建设人人有责人人尽责人人享有的| 久久久精品免费免费高清| 亚洲国产av影院在线观看| 又紧又爽又黄一区二区| 丝袜美腿诱惑在线| 国产成人a∨麻豆精品| a级片在线免费高清观看视频| 99热网站在线观看| 三上悠亚av全集在线观看| 亚洲国产欧美一区二区综合| 男男h啪啪无遮挡| 一二三四社区在线视频社区8| 大片电影免费在线观看免费| 国产一级毛片在线| 一区二区三区精品91| 亚洲一区中文字幕在线| 久久女婷五月综合色啪小说| 国产一区二区 视频在线| 久久国产精品男人的天堂亚洲| 久久精品久久久久久噜噜老黄| 一级黄片播放器| 久久精品国产亚洲av高清一级| 脱女人内裤的视频| 久久av网站| 两个人看的免费小视频| 涩涩av久久男人的天堂| 亚洲av日韩在线播放| 亚洲色图综合在线观看| 亚洲成人免费av在线播放| 观看av在线不卡| 超碰97精品在线观看| 日韩人妻精品一区2区三区| 午夜福利影视在线免费观看| 大香蕉久久网| 午夜免费鲁丝| 亚洲久久久国产精品| 多毛熟女@视频| 成年动漫av网址| 中文精品一卡2卡3卡4更新| 自线自在国产av| 日本一区二区免费在线视频| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 午夜日韩欧美国产| 99国产精品免费福利视频| 欧美少妇被猛烈插入视频| 亚洲国产毛片av蜜桃av| 日本一区二区免费在线视频| 久久精品久久久久久久性| 久久久亚洲精品成人影院| av不卡在线播放| 肉色欧美久久久久久久蜜桃| 91精品国产国语对白视频| 亚洲av美国av| 亚洲国产精品一区三区| 日韩视频在线欧美| 精品欧美一区二区三区在线| 精品人妻一区二区三区麻豆| 日本一区二区免费在线视频| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 精品国产国语对白av| 日韩,欧美,国产一区二区三区| 欧美少妇被猛烈插入视频| 欧美97在线视频| 国产男人的电影天堂91| videos熟女内射| 十分钟在线观看高清视频www| 午夜视频精品福利| 亚洲一区二区三区欧美精品| 亚洲午夜精品一区,二区,三区| 久久人妻福利社区极品人妻图片 | 男女床上黄色一级片免费看| 欧美另类一区| 免费看不卡的av| 青春草亚洲视频在线观看| 一区二区日韩欧美中文字幕| 各种免费的搞黄视频| 黑人巨大精品欧美一区二区蜜桃| 成人黄色视频免费在线看| 麻豆乱淫一区二区| 一区二区三区精品91| 亚洲精品成人av观看孕妇| 久久99一区二区三区| 一区二区三区乱码不卡18| 在线精品无人区一区二区三| 欧美少妇被猛烈插入视频| 日韩 欧美 亚洲 中文字幕| 极品人妻少妇av视频| 国产高清视频在线播放一区 | 日韩大码丰满熟妇| 一边摸一边抽搐一进一出视频| 老司机影院毛片| 国产免费一区二区三区四区乱码| 欧美人与善性xxx| 丰满迷人的少妇在线观看| 一区二区日韩欧美中文字幕| 最黄视频免费看| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美一区二区综合| 高清欧美精品videossex| 日韩熟女老妇一区二区性免费视频| 精品第一国产精品| 午夜福利,免费看| 人人澡人人妻人| 国产午夜精品一二区理论片| 国产成人欧美在线观看 | 国产日韩欧美在线精品| xxx大片免费视频| 欧美激情 高清一区二区三区| 青青草视频在线视频观看| 中文字幕人妻丝袜一区二区| 高清不卡的av网站| 久久中文字幕一级| 亚洲欧美一区二区三区国产| 免费高清在线观看视频在线观看| 亚洲精品中文字幕在线视频| 亚洲熟女精品中文字幕| 18禁观看日本| 亚洲免费av在线视频| 丝袜喷水一区| 男女高潮啪啪啪动态图| 免费在线观看完整版高清| 国产99久久九九免费精品| 国产一区二区激情短视频 | 我要看黄色一级片免费的| 亚洲国产精品一区二区三区在线| 亚洲精品日本国产第一区| 麻豆av在线久日| 婷婷成人精品国产| 国产成人一区二区三区免费视频网站 | 亚洲国产精品国产精品| 大码成人一级视频| 午夜福利,免费看| 亚洲 欧美一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 满18在线观看网站| 亚洲熟女毛片儿| 丝袜人妻中文字幕| 欧美成人午夜精品| 国产黄频视频在线观看| 国产视频一区二区在线看| 男女边摸边吃奶| 麻豆乱淫一区二区| 午夜av观看不卡| 亚洲国产精品国产精品| 女人精品久久久久毛片| 国产又爽黄色视频| 天堂8中文在线网| 色网站视频免费| 97人妻天天添夜夜摸| 婷婷色麻豆天堂久久|