• <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在线不卡| 国产色婷婷99| 亚洲第一区二区三区不卡| 男人舔女人下体高潮全视频| 亚洲aⅴ乱码一区二区在线播放| 51国产日韩欧美| 午夜精品一区二区三区免费看| 人人妻人人澡人人爽人人夜夜 | 麻豆精品久久久久久蜜桃| 国产大屁股一区二区在线视频| 国产一区二区三区综合在线观看 | 寂寞人妻少妇视频99o| 欧美bdsm另类| 国产成人a∨麻豆精品| 久热久热在线精品观看| 看免费成人av毛片| 国产精品福利在线免费观看| 免费观看性生交大片5| 久久精品国产亚洲网站| 可以在线观看毛片的网站| 亚洲国产高清在线一区二区三| 五月伊人婷婷丁香| 欧美成人一区二区免费高清观看| videos熟女内射| 国产成人精品福利久久| 日韩精品有码人妻一区| 成人午夜高清在线视频| 内地一区二区视频在线| 波野结衣二区三区在线| 日本wwww免费看| 亚洲欧美日韩卡通动漫| av线在线观看网站| 女人久久www免费人成看片| 黄片无遮挡物在线观看| 免费观看性生交大片5| 亚洲成人一二三区av| 日日摸夜夜添夜夜添av毛片| 日韩欧美精品v在线| 久久精品综合一区二区三区| 日韩一区二区三区影片| av在线亚洲专区| 国产极品天堂在线| 亚洲精品久久午夜乱码| 狠狠精品人妻久久久久久综合| 十八禁网站网址无遮挡 | 精品亚洲乱码少妇综合久久| 天堂中文最新版在线下载 | 国产精品爽爽va在线观看网站| 久久久久精品性色| 中文字幕av在线有码专区| 日韩精品青青久久久久久| 欧美另类一区| 国产精品国产三级国产专区5o| 美女被艹到高潮喷水动态| 身体一侧抽搐| 精品久久久久久久久av| 51国产日韩欧美| 少妇的逼水好多| 午夜日本视频在线| 亚洲av一区综合| 亚洲成人一二三区av| 亚洲,欧美,日韩| 免费看a级黄色片| 中国国产av一级| 人妻系列 视频| 少妇的逼水好多| 亚洲av中文av极速乱| 亚洲内射少妇av| 国产精品美女特级片免费视频播放器| 色5月婷婷丁香| 国产三级在线视频| 69av精品久久久久久| 亚洲av不卡在线观看| 特级一级黄色大片| av播播在线观看一区| 亚洲av电影在线观看一区二区三区 | 美女内射精品一级片tv| 国产成人精品一,二区| 久久久久久久久久人人人人人人| 五月玫瑰六月丁香| 国产乱来视频区| 97超碰精品成人国产| 亚洲av国产av综合av卡| 嫩草影院新地址| 嘟嘟电影网在线观看| 欧美+日韩+精品| av免费观看日本| 亚洲精品久久午夜乱码| 国产 一区 欧美 日韩| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 久久久久性生活片| 欧美不卡视频在线免费观看| 亚洲欧美成人综合另类久久久| 国产精品久久久久久av不卡| 亚洲在线自拍视频| 国产黄色免费在线视频| 免费大片18禁| 极品少妇高潮喷水抽搐| 国产一级毛片七仙女欲春2| 午夜激情久久久久久久| 久久精品久久久久久久性| 亚洲欧美日韩东京热| 免费播放大片免费观看视频在线观看| 日韩欧美精品免费久久| 中文字幕人妻熟人妻熟丝袜美| 精品国产露脸久久av麻豆 | 国产在线一区二区三区精| 日本免费在线观看一区| 淫秽高清视频在线观看| 丝瓜视频免费看黄片| 全区人妻精品视频| 精品一区在线观看国产| 亚洲不卡免费看| 三级男女做爰猛烈吃奶摸视频| 亚洲av男天堂| 国产麻豆成人av免费视频| 波多野结衣巨乳人妻| 国产精品国产三级专区第一集| 综合色丁香网| 国产大屁股一区二区在线视频| 我的老师免费观看完整版| 国产精品三级大全| 久热久热在线精品观看| 精品不卡国产一区二区三区| 国产亚洲最大av| 美女高潮的动态| 国内精品一区二区在线观看| 亚洲av成人av| 晚上一个人看的免费电影| 国产精品一区二区三区四区免费观看| 国产成年人精品一区二区| 亚洲国产精品成人久久小说| 国产精品99久久久久久久久| av免费在线看不卡| 亚洲成人久久爱视频| 久久久久久久久久成人| xxx大片免费视频| 秋霞在线观看毛片| 麻豆成人av视频| 国产不卡一卡二| 午夜老司机福利剧场| 成人一区二区视频在线观看| 26uuu在线亚洲综合色| 欧美xxxx黑人xx丫x性爽| 真实男女啪啪啪动态图| 亚洲自拍偷在线| 卡戴珊不雅视频在线播放| 国产综合懂色| 少妇人妻一区二区三区视频| 亚洲精品日本国产第一区| 亚洲精华国产精华液的使用体验| 日韩制服骚丝袜av| 午夜精品在线福利| 精品酒店卫生间| 99久久中文字幕三级久久日本| 国产成人aa在线观看| 欧美激情国产日韩精品一区| 久久精品夜夜夜夜夜久久蜜豆| 日本欧美国产在线视频| 亚洲最大成人中文| 一级a做视频免费观看| 国产精品国产三级国产专区5o| 欧美不卡视频在线免费观看| 国产亚洲最大av| 亚洲精品国产av蜜桃| 国产熟女欧美一区二区| 国产黄色小视频在线观看| 久久精品国产鲁丝片午夜精品| 午夜福利网站1000一区二区三区| 国产精品久久久久久av不卡| 国产精品麻豆人妻色哟哟久久 | 中文字幕制服av| av.在线天堂| 丰满乱子伦码专区| 亚洲国产精品国产精品| 国产黄频视频在线观看| 在线 av 中文字幕| 国产成人一区二区在线| 精品午夜福利在线看| 国产69精品久久久久777片| 久久久久久久久久人人人人人人| 国产综合精华液| 天堂影院成人在线观看| 国产精品福利在线免费观看| 午夜福利视频1000在线观看| 精品酒店卫生间| 中文精品一卡2卡3卡4更新| 日韩一区二区三区影片| 精品一区二区三区人妻视频| 91精品伊人久久大香线蕉| 国产日韩欧美在线精品| 国产一区二区在线观看日韩| 毛片女人毛片| 岛国毛片在线播放| 国产精品一及| 国产乱来视频区| 色网站视频免费| 色综合站精品国产| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 麻豆成人午夜福利视频| 最近中文字幕2019免费版| 亚洲精品国产成人久久av| 中文字幕免费在线视频6| 国产女主播在线喷水免费视频网站 | 大又大粗又爽又黄少妇毛片口| 欧美日韩在线观看h| 亚洲综合精品二区| 伦理电影大哥的女人| 精品午夜福利在线看| 国产精品久久久久久久久免| 精品一区二区三区人妻视频| 一区二区三区免费毛片| 久久鲁丝午夜福利片| 又粗又硬又长又爽又黄的视频| 免费看av在线观看网站| 国产免费一级a男人的天堂| 校园人妻丝袜中文字幕| 80岁老熟妇乱子伦牲交| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 欧美精品一区二区大全| 欧美日韩国产mv在线观看视频 | 欧美一级a爱片免费观看看| 黄片无遮挡物在线观看| 亚洲人与动物交配视频| 身体一侧抽搐| 精品久久久精品久久久| 国产一区二区三区综合在线观看 | 国产精品一及| 国产91av在线免费观看| av在线观看视频网站免费| 最后的刺客免费高清国语| 亚洲内射少妇av| 中文字幕免费在线视频6| 国产永久视频网站| 插阴视频在线观看视频| 伦理电影大哥的女人| 91久久精品国产一区二区三区| 久久精品国产亚洲av天美| 日韩一本色道免费dvd| av国产久精品久网站免费入址| 精品久久久久久久人妻蜜臀av| 亚洲国产av新网站| 中文字幕制服av| 亚洲精品亚洲一区二区| 亚洲国产色片| 永久网站在线| 美女内射精品一级片tv| 九草在线视频观看| 欧美性感艳星| 中文字幕av在线有码专区| 免费电影在线观看免费观看| 别揉我奶头 嗯啊视频| 亚洲乱码一区二区免费版| 亚洲精品日本国产第一区| 国产精品.久久久| 久久精品综合一区二区三区| 联通29元200g的流量卡| 日日撸夜夜添| 国产单亲对白刺激| 自拍偷自拍亚洲精品老妇| 在线 av 中文字幕| 精品熟女少妇av免费看| 性色avwww在线观看| 一个人看视频在线观看www免费| 高清av免费在线| 伊人久久国产一区二区| 又大又黄又爽视频免费| 婷婷色av中文字幕| 自拍偷自拍亚洲精品老妇| 精品一区二区三区视频在线| 国产午夜福利久久久久久| 国产av国产精品国产| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久 | 国产午夜福利久久久久久| 免费观看性生交大片5| 少妇的逼水好多| 国产成年人精品一区二区| 亚洲美女搞黄在线观看| av免费观看日本| 秋霞在线观看毛片| 我的老师免费观看完整版| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单| 亚洲国产精品成人综合色| 赤兔流量卡办理| 精品午夜福利在线看| 99久久人妻综合| 自拍偷自拍亚洲精品老妇| 国产人妻一区二区三区在| 国产高清不卡午夜福利| 91久久精品国产一区二区成人| 欧美高清成人免费视频www| 日韩欧美三级三区| 免费黄频网站在线观看国产| 搡女人真爽免费视频火全软件| 人妻系列 视频| 国产高潮美女av| av国产久精品久网站免费入址| 天天一区二区日本电影三级| 国产黄片视频在线免费观看| 伊人久久精品亚洲午夜| 青青草视频在线视频观看| 赤兔流量卡办理| 一级片'在线观看视频| 久久久久久久久久久丰满| av国产免费在线观看| 久久久久性生活片| 国产极品天堂在线| 麻豆乱淫一区二区| 亚洲真实伦在线观看| 亚洲成人久久爱视频| 亚洲欧美一区二区三区黑人 | 亚洲激情五月婷婷啪啪| 国产麻豆成人av免费视频| 三级国产精品欧美在线观看| 在现免费观看毛片| 免费观看在线日韩| 国产在视频线精品| 国产av码专区亚洲av| 性插视频无遮挡在线免费观看| 麻豆成人av视频| 激情 狠狠 欧美| 精品国内亚洲2022精品成人| 床上黄色一级片| 日本一本二区三区精品| 亚洲18禁久久av| 亚洲国产精品国产精品| 美女国产视频在线观看| 内地一区二区视频在线| 精品久久久久久电影网| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| 色播亚洲综合网| 日韩国内少妇激情av| 青春草亚洲视频在线观看| 久久精品熟女亚洲av麻豆精品 | 直男gayav资源| av国产免费在线观看| 国产单亲对白刺激| 欧美日韩精品成人综合77777| 国产单亲对白刺激| 美女主播在线视频| 欧美 日韩 精品 国产| 国产精品国产三级专区第一集| 久久精品夜色国产| 国产乱人视频| 国产中年淑女户外野战色| av网站免费在线观看视频 | 欧美最新免费一区二区三区| 欧美日韩国产mv在线观看视频 | 免费大片黄手机在线观看| 嫩草影院精品99| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 国产精品人妻久久久久久| 激情五月婷婷亚洲| 日日摸夜夜添夜夜添av毛片| 一个人免费在线观看电影| 欧美人与善性xxx| 欧美日韩视频高清一区二区三区二| 国产午夜精品论理片| 噜噜噜噜噜久久久久久91| 人人妻人人看人人澡| 一级毛片电影观看| 亚洲欧美成人精品一区二区| av在线天堂中文字幕| 尾随美女入室| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 欧美日韩视频高清一区二区三区二| 麻豆成人av视频| 免费大片18禁| 免费看光身美女| 一夜夜www| 成人高潮视频无遮挡免费网站| 啦啦啦中文免费视频观看日本| 日韩,欧美,国产一区二区三区| 一区二区三区四区激情视频| 网址你懂的国产日韩在线| 亚洲国产精品成人综合色| 网址你懂的国产日韩在线| 国产精品久久视频播放| 国产熟女欧美一区二区| 精品一区二区三区视频在线| 亚洲精品日韩在线中文字幕| 麻豆av噜噜一区二区三区| 国产黄频视频在线观看| 国产精品蜜桃在线观看| 国产精品av视频在线免费观看| 免费黄网站久久成人精品| 亚洲熟女精品中文字幕| 精品久久久久久成人av| 99久久中文字幕三级久久日本| 国产色爽女视频免费观看| 亚洲av福利一区| 亚洲国产精品国产精品| 777米奇影视久久| 一级毛片aaaaaa免费看小| 欧美另类一区| 搞女人的毛片| 最近中文字幕高清免费大全6| 国产精品一区二区三区四区免费观看| 国模一区二区三区四区视频| 在线a可以看的网站| 97热精品久久久久久| 肉色欧美久久久久久久蜜桃 | 国产成人精品一,二区| 国产av不卡久久| 午夜精品一区二区三区免费看| 亚洲av成人av| 在线播放无遮挡| 超碰av人人做人人爽久久| 国产一区二区三区综合在线观看 | 乱人视频在线观看| 欧美激情国产日韩精品一区| 日韩av在线大香蕉| av.在线天堂| 午夜爱爱视频在线播放| 街头女战士在线观看网站| 18禁动态无遮挡网站| 少妇丰满av| av黄色大香蕉| 日韩av免费高清视频| 免费看不卡的av| 亚洲av一区综合| 色综合色国产| 毛片女人毛片| 2022亚洲国产成人精品| 日本黄色片子视频| 在线观看一区二区三区| 国产精品麻豆人妻色哟哟久久 | 国内精品一区二区在线观看| 一区二区三区乱码不卡18| 国产成人精品久久久久久| 能在线免费观看的黄片| 黄片无遮挡物在线观看| 亚洲图色成人| 九色成人免费人妻av| 国产又色又爽无遮挡免| 国精品久久久久久国模美| 精品不卡国产一区二区三区| 久99久视频精品免费| ponron亚洲| 国产中年淑女户外野战色| 哪个播放器可以免费观看大片| 国产一级毛片在线| 亚洲av男天堂| 亚洲丝袜综合中文字幕| 天美传媒精品一区二区| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 两个人的视频大全免费| 日韩 亚洲 欧美在线| 午夜福利高清视频| 久热久热在线精品观看| 美女黄网站色视频| 免费看不卡的av| 久久97久久精品| 99视频精品全部免费 在线| 国产精品久久久久久精品电影小说 | 人人妻人人看人人澡| 亚洲成人精品中文字幕电影| 午夜精品国产一区二区电影 | 成人漫画全彩无遮挡| 色综合色国产| 在线 av 中文字幕| 黄片无遮挡物在线观看| 伊人久久国产一区二区| 国产精品日韩av在线免费观看| 熟妇人妻久久中文字幕3abv| 一本久久精品| 你懂的网址亚洲精品在线观看| 日韩 亚洲 欧美在线| 97人妻精品一区二区三区麻豆| 又爽又黄a免费视频| 夫妻午夜视频| 国产毛片a区久久久久| 2021少妇久久久久久久久久久| 婷婷色麻豆天堂久久| 久久久久久久国产电影| 午夜激情久久久久久久| 内射极品少妇av片p| 欧美激情在线99| 国产成人精品久久久久久| 国产精品国产三级国产专区5o| 精品国产三级普通话版| 成人一区二区视频在线观看| 69人妻影院| 精品熟女少妇av免费看| 韩国高清视频一区二区三区| 亚洲av男天堂| 91精品一卡2卡3卡4卡| 天天躁日日操中文字幕| 在线观看av片永久免费下载| 中文字幕免费在线视频6| 国产亚洲91精品色在线| 日韩国内少妇激情av| 男插女下体视频免费在线播放| 神马国产精品三级电影在线观看| 别揉我奶头 嗯啊视频| 色哟哟·www| 看免费成人av毛片| 亚洲精品日本国产第一区| 日韩欧美精品v在线| 在线观看人妻少妇| 99九九线精品视频在线观看视频| 人体艺术视频欧美日本| 久久久久性生活片| 亚洲av男天堂| 日韩成人伦理影院| 18禁裸乳无遮挡免费网站照片| 神马国产精品三级电影在线观看| 国产在线一区二区三区精| 国产一区二区亚洲精品在线观看| 少妇的逼水好多| 91精品伊人久久大香线蕉| 亚洲精品aⅴ在线观看| 三级毛片av免费| 国产精品av视频在线免费观看| 午夜免费观看性视频| 乱人视频在线观看| 高清午夜精品一区二区三区| 麻豆精品久久久久久蜜桃| 美女国产视频在线观看| 国产亚洲精品av在线| 最后的刺客免费高清国语| 日产精品乱码卡一卡2卡三| 亚洲精品一区蜜桃| 日本猛色少妇xxxxx猛交久久| 日本熟妇午夜| 午夜精品在线福利| 亚洲精品成人久久久久久| 中文字幕久久专区| 国精品久久久久久国模美| 边亲边吃奶的免费视频| 精品久久久久久成人av| 色综合站精品国产| 三级男女做爰猛烈吃奶摸视频| 99热全是精品| 国产 亚洲一区二区三区 | 精品一区二区三卡| 两个人视频免费观看高清| 边亲边吃奶的免费视频| 欧美高清性xxxxhd video| 丰满乱子伦码专区| 色网站视频免费| 国产乱来视频区| 亚洲一级一片aⅴ在线观看| 看黄色毛片网站| 国模一区二区三区四区视频| 自拍偷自拍亚洲精品老妇| 又黄又爽又刺激的免费视频.| 一个人免费在线观看电影| 三级国产精品欧美在线观看| 午夜精品国产一区二区电影 | 在线观看免费高清a一片| 中文字幕av成人在线电影| 久久久久久久久久久丰满| 乱人视频在线观看| 精品熟女少妇av免费看| 欧美日韩亚洲高清精品| 国产精品三级大全| av天堂中文字幕网| 久久精品久久久久久久性| 五月玫瑰六月丁香| 国产毛片a区久久久久| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 亚洲无线观看免费| 亚洲四区av| 美女大奶头视频| 国产淫语在线视频| 六月丁香七月| 日本与韩国留学比较| 观看免费一级毛片| 一区二区三区四区激情视频| 中文字幕av成人在线电影| 天堂影院成人在线观看| 欧美一区二区亚洲| 亚洲第一区二区三区不卡| 91精品国产九色| 男女下面进入的视频免费午夜| 可以在线观看毛片的网站| 毛片一级片免费看久久久久| 免费高清在线观看视频在线观看| 国产成人精品一,二区| 一级毛片久久久久久久久女| 1000部很黄的大片| 亚洲精品一区蜜桃| 少妇裸体淫交视频免费看高清| 午夜福利网站1000一区二区三区| 日本三级黄在线观看| 久久综合国产亚洲精品| 午夜福利网站1000一区二区三区| 别揉我奶头 嗯啊视频| 在线观看美女被高潮喷水网站| 边亲边吃奶的免费视频| 伊人久久国产一区二区| 青青草视频在线视频观看| 两个人的视频大全免费| 色尼玛亚洲综合影院| 偷拍熟女少妇极品色| 欧美精品一区二区大全| 又爽又黄a免费视频| 一夜夜www| 亚洲精品自拍成人| 国产一区有黄有色的免费视频 | 久久这里有精品视频免费| 亚洲aⅴ乱码一区二区在线播放|