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

    不同氮處理下水稻劍葉葉寬的全基因組關(guān)聯(lián)分析

    2017-07-31 23:47:33高易宏燕金香涂政軍冷語佳陳龍黃李超代麗萍張光恒朱麗胡江任德勇郭龍彪錢前王丹英曾大力
    中國農(nóng)業(yè)科學(xué) 2017年14期
    關(guān)鍵詞:劍葉葉寬氮肥

    高易宏,燕金香,涂政軍,冷語佳,陳龍,黃李超,代麗萍,張光恒,朱麗,胡江,任德勇,郭龍彪,錢前,王丹英,曾大力

    ?

    不同氮處理下水稻劍葉葉寬的全基因組關(guān)聯(lián)分析

    高易宏,燕金香,涂政軍,冷語佳,陳龍,黃李超,代麗萍,張光恒,朱麗,胡江,任德勇,郭龍彪,錢前,王丹英,曾大力

    (中國水稻研究所/水稻生物學(xué)國家重點實驗室,杭州 310006)

    【目的】探索水稻劍葉葉寬調(diào)控及其對氮肥響應(yīng)的遺傳機理,為氮高效水稻新品種的培育提供新的種質(zhì)資源和基因標(biāo)記?!痉椒ā恳?34份水稻地方種質(zhì)資源為關(guān)聯(lián)分析材料,通過基因組重測序發(fā)掘獲得了3 356 591個分布于全基因組的高密度SNP位點(single nucleotide polymorphism,SNP)。在大田栽培條件下,以施氮量為主區(qū),品種為裂區(qū)的設(shè)計,設(shè)計低氮(不施氮肥,N0)、正常氮(純氮96 kg·hm-2,N1)和高氮(純氮192 kg·hm-2,N2)3種氮肥處理。于水稻成熟期分別調(diào)查水稻劍葉葉寬在低、中、高3種氮肥處理下的表現(xiàn)及響應(yīng),結(jié)合EMMAX軟件計算群體親緣關(guān)系矩陣和EIGENSOFT軟件分析群體結(jié)構(gòu),采用納入親緣關(guān)系矩陣及群體結(jié)構(gòu)的混合線性模型開展全基因組關(guān)聯(lián)分析?!窘Y(jié)果】水稻劍葉葉寬在N0、N1、N2 3種氮肥處理下均呈正態(tài)分布,并表現(xiàn)豐富的變異。劍葉葉寬受品種差異及氮水平的影響,且與施氮水平呈顯著正相關(guān)。在N0、N1、N2氮處理下共檢測到14個與劍葉葉寬顯著相關(guān)的SNP位點。其中,低氮處理下檢測到的SNP位點的最小等位基因頻率均大于0.46,表明此類SNP在關(guān)聯(lián)群體中廣泛存在;中氮和高氮水平下檢測到的SNP位點的最小等位基因頻率均較小,是一類較為稀有的SNP位點。位于第12染色體上的一個SNP(chr12:15 066 507)位點在正常氮及高氮處理下均被檢測到,在高氮處理下還檢測到的另一顯著性位點,其候選區(qū)間內(nèi)包含一候選基因,該基因與業(yè)已報道的葉寬性狀相關(guān)基因同屬于細(xì)胞色素P450家族。根據(jù)不同氮處理下劍葉葉寬的響應(yīng),鑒定出20與低氮響應(yīng)有關(guān)的SNP位點,8個位點與高氮響應(yīng)有關(guān)。其中與高氮響應(yīng)的顯著性位點中,位于第1染色體顯著性峰候選區(qū)間包含業(yè)已克隆的與氮素利用相關(guān)的基因?!窘Y(jié)論】通過全基因組關(guān)聯(lián)分析共檢測到42個與劍葉葉寬及其在不同氮處理下葉寬響應(yīng)相關(guān)的顯著性關(guān)聯(lián)位點。

    水稻;全基因組關(guān)聯(lián)分析;葉寬;氮肥

    0 引言

    【研究意義】氮是植物生長發(fā)育中至關(guān)重要的營養(yǎng)元素,合理的氮肥施用將顯著的增加作物產(chǎn)量[1],氮營養(yǎng)供應(yīng)不充足則會導(dǎo)致作物減產(chǎn),而過量施用氮肥會導(dǎo)致環(huán)境污染等問題。因而氮素的高效利用是當(dāng)前水稻育種關(guān)注的重要問題,通過選育氮高效吸收利用的水稻品種,從而減少田間氮肥施用以提高水稻產(chǎn)量已成為水稻育種的重要目標(biāo)之一[2]?!厩叭搜芯窟M展】水稻的氮素利用效率(nitrogen use efficiency,NUE)受品種、氮肥用量、氮肥施用方法和土壤含氮量等因素的影響。近年來,通過對水稻氮響應(yīng)相關(guān)突變體的研究和QTL定位等方式克隆與鑒定了多個與NUE相關(guān)的基因及位點。Zhang等[3]鑒定了一個位于第12染色體上有關(guān)氮耐受性的QTL位點,并克隆了。Hu等[4]通過在高世代回交群體分離鑒定了一個硝酸鹽轉(zhuǎn)運基因,該基因通過增強硝酸鹽的攝取、根向莖轉(zhuǎn)運和氮響應(yīng)基因的表達(dá)從而提高氮素利用率。Sun等[5]通過分析RIL群體在不同氮濃度下株高和分蘗數(shù)的變化鑒定了一個水稻主效的NUE數(shù)量性狀位點,實為已報道的直立密穗基因,攜帶的植株表現(xiàn)出對氮吸收和同化的增強,從而增加水稻產(chǎn)量。在長期的進化過程中,植物形成了適應(yīng)不同生存環(huán)境的生理機制和分子機理來響應(yīng)環(huán)境的變化。在應(yīng)對缺氮環(huán)境時,植物提高NUE的主要方式是通過提高氮的吸收、轉(zhuǎn)運、同化和再活化等方式,且與環(huán)境之間存在復(fù)雜的互作及遺傳網(wǎng)絡(luò)關(guān)系[6-9]。Tong等[10]通過在2種氮處理條件下的株高、穗數(shù)、葉綠素含量、干重、濕重和產(chǎn)量變化等指標(biāo)鑒定了31個氮響應(yīng)相關(guān)的QTL位點。Cho等[11]在兩種氮濃度下,通過水稻發(fā)育后期籽粒中氮含量、秸稈中氮含量、收獲指數(shù)等多個指標(biāo)鑒定了多個氮響應(yīng)相關(guān)的QTL位點。同樣,在不同的氮水平下水稻葉片呈現(xiàn)明顯的葉形變化,而水稻葉片的長、寬和厚的改變也直接影響水稻的光合作用[12]。殷春淵等[13]研究表明,氮肥處理對水稻葉片的影響主要表現(xiàn)在倒二葉葉長和劍葉葉寬2個特性上。因此,葉片形態(tài)結(jié)構(gòu)的變化也是水稻對不同水平氮響應(yīng)的指標(biāo)。近年來,全基因組關(guān)聯(lián)分析(genome-wide association study,GWAS)成功的應(yīng)用在一些模式生物研究和農(nóng)作物的育種當(dāng)中[14~18]?!颈狙芯壳腥朦c】在水稻中,利用GWAS分析了包括產(chǎn)量、品質(zhì)、抽穗期、株形等多種重要的農(nóng)藝性狀,同樣也對病蟲害抗性等生物脅迫和鹽堿耐受性等非生物脅迫性狀進行了全基因組關(guān)聯(lián)分析,并鑒定了許多重要的SNP位點[14,16-17,19-22]。但是,與氮素利用相關(guān)的GWAS研究尚未見報道?!緮M解決的關(guān)鍵問題】本研究通過選取134份來自不同地區(qū)的水稻品種作為關(guān)聯(lián)分析的群體,利用在不同氮肥條件下葉寬的變化,開展水稻葉寬及氮響應(yīng)相關(guān)SNP位點分析,旨在為進一步揭示水稻葉片形態(tài)調(diào)控和氮素高效利用的生理機制提供基礎(chǔ)。

    1 材料與方法

    1.1 試驗材料

    選擇保存于中國水稻所(CNRRI),國際水稻所(IRRI)種質(zhì)資源庫及來自多個國家和地區(qū)的134個地方品種,其中包括41份來自國內(nèi)各省農(nóng)科院的地方品種;36份來自20個國家的地方品種;57份國內(nèi)大面積推廣種植品種。

    1.2 試驗設(shè)置及性狀考察

    試驗于2015年在中國水稻研究所實驗基地(浙江,富陽)進行,試驗田土壤為水稻土(土壤含有機質(zhì)36.9 g·kg-1、全氮2.73 g·kg-1、速效鉀104 mg·kg-1、速效磷52 mg·kg-1、堿解氮156 mg·kg-1,pH 6.24)。試驗設(shè)3個氮肥水平:低氮(N0,不另施氮肥)、正常氮(N1,純氮96 kg·hm-2)和高氮(N2,純氮192 kg·hm-2),以施氮量為主區(qū),品種為裂區(qū),3次重復(fù)。主區(qū)面積5.0 m×23.0 m,播種前將其做成3條畦寬1.5 m,溝寬25 cm,長23 m的畦。水稻種于畦上,每品種種植4行,每行6株,株行距為20 cm×20 cm。不同氮處理間作埂隔離,并用塑料薄膜覆蓋埂體,保證單獨排灌。

    試驗用關(guān)聯(lián)群體于2015年5月10日采用點直播栽培,每穴播3—5粒種子,2葉1心時定苗2株/穴。以尿素為氮肥,按基肥﹕分蘗肥﹕穗肥 = 5﹕3﹕2分3次施入,其中基肥(播種前1 d)與分蘗肥(播種后25 d)統(tǒng)一施用,穗肥由于品種間生育期的差異,于孕穗始期分品種單獨施用。6月中旬至9月中旬每隔3d檢查品種的生育期,標(biāo)記孕穗始期相近的品種,將穗肥小心撒于每個標(biāo)記品種葉片上。為防止水分流動導(dǎo)致的串肥,穗肥施用期間(6月中旬至9月中旬)保持畦面濕潤但無水層,僅溝中有水。磷、鉀施用量為P2O575 kg·hm-2、K2O 150 kg·hm-2,其中磷肥為底肥一次施入,鈣、鎂、鉀肥按基、蘗肥各50%分2次施入。其他栽培管理措施同一般大田。

    在抽穗后10 d對劍葉葉寬進行調(diào)查,每品種選取中間的6株無邊際效應(yīng)的植株進行考察。參考王蘭等[23]方法量取每片葉的最寬處為葉寬,每株取主莖分別記錄劍葉的葉寬。

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

    1.3.1 性狀相關(guān)性分析 使用Excel 2010進行數(shù)據(jù)的錄入和整理,并計算其平均值和標(biāo)準(zhǔn)差。使用SAS 9.2進行進行不同氮處理下葉寬變化的相關(guān)分析和方差分析。

    1.3.2 基因型測定 采用CTAB抽提法提取134份品種幼葉基因組DNA[24],檢測合格后用于高通量測序。每品種測序reads數(shù)據(jù)量在5Gb以上,相當(dāng)于10×以上的基因組測序。從reads到SNP的檢測參考Zhou等[17]方法進行。SNP使用ANNOVAR (Version: 2013-08-23)軟件包進行注釋[25]。

    1.3.3 群體結(jié)構(gòu)與親緣關(guān)系分析 根據(jù)SNP信息,使用PHYLIP 3.68軟件包基于鄰接法構(gòu)建品種間系統(tǒng)進化樹[26],使用在線工具iTOL(http://itol.embl.de/)進行進化樹的優(yōu)化。使用EIGENSOFT軟件中的smartpca程序包進行主成分分析(PCA)[27]。

    1.3.4 葉寬性狀的全基因組關(guān)聯(lián)分析 采用基于混合線性模型(mixed linear model,MLM)的EMMAX(Efficient Mixed-Model Association eXpedited)軟件包進行目標(biāo)性狀與SNP標(biāo)記間的關(guān)聯(lián)分析[20, 28],將EMMAX軟件計算的親緣關(guān)系K矩陣及基于EIGENSOFT軟件計算的群體結(jié)構(gòu)Q納入分析模型中,關(guān)聯(lián)分析參數(shù)設(shè)定為:emmax -v,-d 10,-t,-p,-k,-c,-o。在考慮關(guān)聯(lián)分析群體大小和SNP數(shù)量的情況下,對于關(guān)聯(lián)分析的顯著性,采取10-6作為寬松閾值,而進一步的嚴(yán)格閾值采取Bonferroni校驗閾值[29]。

    2 結(jié)果

    2.1 3種氮處理條件下葉寬性狀多樣性及其變化

    在134份材料中,劍葉葉寬表現(xiàn)豐富的變異,正常氮處理下變幅為9.93—29.00 mm。如木邦谷、羅大穗等地方品種表現(xiàn)為株型較大,分蘗少,成熟期葉片少,劍葉葉寬大于20 mm;而地方品種英德大葉拍、合江23等表現(xiàn)為株型小,成熟期葉片多,葉片呈細(xì)長狀,葉寬小于10 mm。由此表明,此關(guān)聯(lián)群體具有充分的表型變異,能夠進行基于葉寬的全基因組關(guān)聯(lián)分析。

    對裂區(qū)設(shè)計整體進行方差分析,變異來源于氮肥使用量、品種及品種與氮肥交互作用的值均小于0.0001,而區(qū)組與氮處理間的交互作用為0.2444,這說明葉寬變異主要受氮處理及品種差異引起,不同品種在不同氮濃度下展現(xiàn)不同的葉寬變化。對134份品種成熟期劍葉寬調(diào)查結(jié)果表明,在低、中、高3種氮處理下,劍葉葉寬均呈正態(tài)分布(圖1-a—圖1-c),從頻率分布來看,中氮和高氮處理下劍葉寬為15—18 mm的品種明顯增多。在反映氮響應(yīng)變化中,正常氮處理相對于低氮處理表現(xiàn)為葉寬變大,而相對于高氮處理變化要?。▓D1-d)。相關(guān)性分析亦表明,氮處理與劍葉葉寬呈顯著正相關(guān)(表1),這表明,葉寬受外界氮水平的影響,且葉寬的變化能夠反映不同氮處理。

    根據(jù)不同氮處理下葉寬相對變化,將葉寬變化分為四種響應(yīng)模式(圖2-e)。有35.2%的品種表現(xiàn)為在高氮條件下葉寬增加而低氮條件下葉寬減小,有18.1%的品種在低氮和高氮條件下劍葉寬均增大,另有36.2%的品種在低氮和高氮條件下葉寬均減小,還有5.7%的品種在低氮條件下劍葉變寬而在高氮條件下反而變窄。由此可知,該關(guān)聯(lián)群體在不同氮處理下葉寬相對變化呈現(xiàn)豐富的變異,可以利用葉寬的相對變化來反映氮素利用效率。

    表1 不同氮濃度下葉寬統(tǒng)計及其相關(guān)性分析

    *、**分別表示差異達(dá)到0.05 和 0.01 的顯著水平

    *,** present significant at 5% and 1% probability levels respectively

    a、b、c分別為低氮、中氮和高氮條件下的葉寬分布;d:3種氮水平下葉寬分布的氣泡圖,氣泡寬度代表株系數(shù),紅色虛線為3種氮條件下葉寬平均值;e:高氮脅迫下葉寬相對變化(N2-N1)和低氮脅迫下葉寬相對變化(N0-N1)的分布

    2.2 關(guān)聯(lián)分析SNP標(biāo)記

    基于134份水稻的測序數(shù)據(jù),結(jié)合使用BWA和GATK軟件進行SNP的提取,使用最小等位基因頻率大于0.03(MAF>0.03)進行SNP的過濾,最終提取到3 356 591個SNP位點,平均每kb有8.6個SNP位點。為進一步分析SNP位點信息,使用ANNOVAR軟件對SNP數(shù)據(jù)注釋(圖2),其中,有2 814個位點注釋為剪接位點相關(guān),此類位點通過影響基因轉(zhuǎn)錄剪接的變化導(dǎo)致基因功能改變;409 387個SNP位點注釋到外顯子區(qū),引起核酸編碼信息的改變,其中59%的SNP為非同義突變,造成編碼信息的改變,從而引起蛋白序列的變化;另外有12 221個位點注釋與終止密碼子相關(guān),終止密碼子的缺失或增加會導(dǎo)致編碼蛋白移碼突變或提前終止。

    圖2 關(guān)聯(lián)分析SNP標(biāo)記分布及其功能注釋

    2.3 群體的遺傳結(jié)構(gòu)分析

    基于134品種重測序的SNP數(shù)據(jù),開展親緣關(guān)系和系統(tǒng)進化分析。結(jié)果表明,此關(guān)聯(lián)群體主要分為秈稻和粳稻兩個亞群。由GAPIT軟件基于SNP計算得到的親緣關(guān)系系數(shù)表明,此134份水稻品種主要分為71份秈稻品種和63份粳稻品種,亞種內(nèi)部分品種的親緣關(guān)系較近,形成小的亞群(圖3-a和圖3-b)?;赟NP數(shù)據(jù)計算品種間遺傳距離和NJ(Neighbor- joining)法構(gòu)建系統(tǒng)發(fā)育樹表明,此關(guān)聯(lián)群體主要為秈粳之間的分歧(圖3-c)。主成分分析(principal component analysis,PCA)同樣也說明了此關(guān)聯(lián)群體秈粳分化的群體結(jié)構(gòu)。此外,基于PCA三點圖可知,秈稻亞種間分化較為集中,粳稻亞種間存在不同程度的分化,較為分散,部分品種間存在小亞群結(jié)構(gòu)(圖3-d)。

    2.4 葉寬性狀的全基因組關(guān)聯(lián)分析

    通過對不同氮處理條件下134份水稻成熟期劍葉葉寬的全基因組關(guān)聯(lián)分析表明,在不同氮處理條件下,與葉寬相關(guān)聯(lián)的SNP位點各不相同(表2)。以-log10(p)大于6為閾值,共檢測到14個與葉寬性狀顯著相關(guān)的SNP位點(圖4)。在低氮和中氮水平下,均檢測到6個顯著的SNP位點與葉寬性狀相關(guān)。而高氮水平只檢測到位于第12染色體上的2個顯著位點,其中一個位點與中氮水平一致,推測該調(diào)控葉寬的位點受氮水平影響較小。

    表2 3種氮處理下劍葉寬的關(guān)聯(lián)位點

    a、b:群體親緣關(guān)系分布圖;c:基于SNP標(biāo)記構(gòu)建的N-J樹;d:主成分分析(PCA)散點分布圖

    在低氮水平下檢測到的顯著性SNP位點最小等位基因頻率均大于0.46,此類SNP在關(guān)聯(lián)群體中分布較為廣泛。而在中氮和高氮水平下檢測到的SNP位點最小等位基因頻率均較小,屬于一類較為稀有的突變位點。在中氮水平下,第12染色體上顯著性SNP位點的值(-lgp)為8.45,高于Bonferroni校正后的閾值水平。對中氮水平下葉寬性狀顯著關(guān)聯(lián)的位點進行候選基因分析時發(fā)現(xiàn),位于第12染色體上的Chr12_14869227顯著位點±200 kb的候選區(qū)間內(nèi),包含一注釋為Cytochrome P450功能相關(guān)蛋白的候選基因,其與業(yè)已報道的葉寬性狀相關(guān)基因同屬于細(xì)胞色素P450家族(圖5)。編碼水稻BR生物合成途徑中的一個關(guān)鍵的BR-6氧化酶,該酶活下降的植株表現(xiàn)為BR缺乏癥狀,而使葉寬受到顯著的影響[30-31]。在高氮水平下,檢測到另一個位于第12染色體上顯著性關(guān)聯(lián)的峰,該區(qū)段包含已知的,編碼一個細(xì)胞壁生物合成和植株生長必須的糖基轉(zhuǎn)移酶家族的纖維素合酶類似蛋白(cellulose synthase-like,CSL),的突變導(dǎo)致水稻葉片變窄、半卷和矮化等表型[32-33]。

    2.5 不同氮水平劍葉葉寬的全基因組關(guān)聯(lián)分析

    進一步使用不同氮水平下葉寬的相對變化值來分析葉寬對氮水平的響應(yīng),共檢測到28個顯著的SNP位點,其中20個與低氮響應(yīng)相關(guān),8個與高氮響應(yīng)相關(guān)(表3、圖6)。在低氮環(huán)境下,葉寬的響應(yīng)位點比高氮環(huán)境下多,這與葉寬變化在低氮環(huán)境下表型變異豐富,而在高氮水平下葉寬變化不明顯有關(guān)。在高氮水平下,位于第7染色體上的顯著性位點Chr7_20639350,其-lg(p)為9.01,超過了Bonferroni檢驗閾值;另外,還檢測到位于第1染色體上的顯著性峰,該區(qū)段存在多個顯著性峰值點(圖7),在此區(qū)段內(nèi)存在一個已知的,該基因與氮素利用相關(guān),表現(xiàn)為葉片明顯早衰,氮的再活化作用被抑制[34-35],推測該基因為此檢測區(qū)間的主效基因。

    左側(cè)為曼哈頓分布點圖,右側(cè)為Q-Q圖;紅色虛線顯示為N1和N2水平下位于第12染色體共同的關(guān)聯(lián)位點。下同

    圖5 第12染色體上在N1水平下的關(guān)聯(lián)顯著位點候選區(qū)間分析

    表3 低氮和高氮水平下劍葉葉寬變化的關(guān)聯(lián)位點

    圖6 劍葉葉寬對低氮和高氮差異響應(yīng)的全基因組關(guān)聯(lián)分析

    圖7 第1染色體上在N2水平的關(guān)聯(lián)顯著位點候選區(qū)間分析

    3 討論

    水稻成熟期劍葉形態(tài)是決定植株光合過程的受光面積和株型結(jié)構(gòu)的關(guān)鍵因素,并進而影響灌漿期光合生物同化量的大小及水稻產(chǎn)量[36-37],因而改善水稻成熟期劍葉形態(tài)成為水稻保持優(yōu)勢株型和高產(chǎn)的重要手段。但水稻葉片形態(tài)是復(fù)雜的受微效多基因控制的數(shù)量性狀,易受光、溫、水、肥等環(huán)境影響[38]。近年來,隨著GWAS在植物領(lǐng)域研究的應(yīng)用,為發(fā)掘重要種質(zhì)資源和快速鑒定新基因提供了新的分析工具[14,16]。

    從GWAS的結(jié)果來看,水稻劍葉葉寬對氮肥脅迫的響應(yīng)更為顯著。在低氮或高氮水平下檢測到的關(guān)聯(lián)位點遠(yuǎn)多于正常氮肥氮肥水平下劍葉葉寬的關(guān)聯(lián)位點,尤其以低氮響應(yīng)為甚,這可能與劍葉葉寬易受氮肥水平的影響以及不同種質(zhì)資源對低氮或高氮的響應(yīng)方式不同有關(guān)。因而,對于易受環(huán)境影響的性狀在開展GWAS分析時,最好能結(jié)合多環(huán)境下的分析或嚴(yán)格控制其他環(huán)境因素的影響。本研究采用3次重復(fù)的裂區(qū)設(shè)計,但沒有作多年/異地的重復(fù),可能會對部分關(guān)聯(lián)位點的精準(zhǔn)選擇產(chǎn)生影響。在進行秈粳型的群體進行關(guān)聯(lián)分析時,將BN (Balding-Nichols)矩陣納入基于MLM模型的關(guān)聯(lián)分析中,進行群體結(jié)構(gòu)的校正,可以有效規(guī)避由于群體結(jié)構(gòu)造成的影響,從而減少假陽性位點的產(chǎn)生。

    氮素的同化和利用對植物至關(guān)重要,前人的研究結(jié)果來看,分別在第1、2、3、5、7、10、11、12染色體上存在已報道的與氮素相關(guān)QTL位點。其中,位于第2染色體的RM5812 標(biāo)記與本研究的低氮脅迫下的檢測到的顯著性位點區(qū)間相近[39]。同樣,在低氮脅迫下檢測到的位于第3、第5染色體上的顯著性位點所在區(qū)間也有QTL位點報道過(http://archive.gramene.org/)。因此,基于顯著性的關(guān)聯(lián)位點,需要進一步對候選區(qū)間內(nèi)的基因進行鑒定、分析和克隆。后期的進一步表達(dá)水平變化或基因敲除研究對關(guān)聯(lián)位點進行功能驗證十分重要,深入展開關(guān)聯(lián)候選基因的功能研究將有助于準(zhǔn)確揭示水稻氮素利用遺傳機制奠定基礎(chǔ)。

    4 結(jié)論

    在不同氮肥水平下對水稻劍葉葉寬進行全基因組關(guān)聯(lián)分析,共檢測到的42個顯著的關(guān)聯(lián)位點中,其中14個與劍葉葉寬相關(guān),28個與氮肥的響應(yīng)相關(guān)。對所檢測到的位點進行候選區(qū)間分析,發(fā)現(xiàn)2個位于候選區(qū)間內(nèi)分別調(diào)控水稻劍葉葉寬及氮響應(yīng)的已知基因;檢測到多條染色體上的顯著位點與已報道與氮素利用相關(guān)QTL區(qū)間相近。此外,基于候選區(qū)間基因的功能注釋,檢測到1個與已知調(diào)控葉寬的同源基因。

    致謝:中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所田志喜研究員和李云海研究員為關(guān)聯(lián)分析提供指導(dǎo)及關(guān)聯(lián)群體材料,在此表示感謝。

    References

    [1] ANDREWS M, RAVEN J A, LEA P J. Do plants need nitrate? The mechanisms by which nitrogen form affects plants., 2013, 163: 174-199.

    [2] ZHANG Q F. Strategies for developing green super rice., 2007, 104(42): 16402-16409.

    [3] ZHANG Y J, TAN L B, ZHU Z F, YUAN L X, XIE D X, SUN C Q. TOND1 confers tolerance to nitrogen deficiency in rice., 2015, 81(3): 367-376.

    [4] HU B, WANG W, OU S, TANG J, LI H, CHE R, ZHANG Z, CHAI X, WANG H, Wang Y, LIANG C, LIU L, PIAO Z, DENG Q, DENG K, XU C, LIANG Y, ZHANG L, LI L, CHU C. Variation in NRT1.1B contributes to nitrate-use divergence between rice subspecies., 2015, 47(7): 834-838.

    [5] SUN H, QIAN Q, WU K, LUO J, WANG S, ZHANG C, MA Y, LIU Q, HUANG X, YUAN Q, HAN R, ZHAO M, DONG G, GUO L, ZHU X, GOU Z, WANG W, WU Y, LIN H, FU X. Heterotrimeric G proteins regulate nitrogen-use efficiency in rice., 2014, 46(2): 652.

    [6] GALLAIS A, HIREL B. An approach to the genetics of nitrogen use efficiency in maize., 2004, 55(396): 295-306.

    [7] HIREL B, LE GOUIS J, NEY B, GALLAIS A. The challenge of improving nitrogen use efficiency in crop plants: towards a more central role for genetic variability and quantitative genetics within integrated approaches., 2007, 58(9): 2369-2387.

    [8] BI Y M, KANT S, CLARKE J, GIDDA S, MING F, XU J, ROCHON A, SHELP B J, HAO L, ZHAO R, MULLEN R T, ZHU T, ROTHSTEIN S J. Increased nitrogen-use efficiency in transgenic rice plants over-expressing a nitrogen-responsive early nodulin gene identified from rice expression profiling., 2009, 32(12):1749-1760.

    [9] XU G, FAN X, MILLER A J. Plant nitrogen assimilation and use efficiency., 2012, 63(63): 153-182.

    [10] TONG H H, MEI H W, YU X Q, XU X Y, LI M S, ZHANG SQ, LUO L J. Identification of related QTLs at late developmental stage in rice (L.) under two nitrogen levels., 2006, 33(5): 458-467.

    [11] CHO Y I, JIANG W, CHIN J H, PIAO Z, CHO Y G, Mccouch S, Koh H J. Identification of QTLs associated with physiological nitrogen use efficiency in rice., 2007, 23(1): 72-79.

    [12] 凌啟鴻. 作物群體質(zhì)量. 上海: 上??萍汲霭嫔? 2000: 42-120.

    LING Q H.. Shanghai: Shanghai Scientific and Technical Publishers, 2000: 42-120. (in Chinese)

    [13] 殷春淵, 王書玉, 薛應(yīng)征, 劉賀梅, 張栩, 孫建權(quán), 王和樂, 胡秀明. 氮肥處理對新稻18號水稻產(chǎn)量及葉片形態(tài)特征的影響. 中國農(nóng)業(yè)科技導(dǎo)報, 2012, 14(3): 101-106.

    YIN Q Y, WANG Y S, XUE Y Z, LIU H M, ZHANG X, SUN J Q, WANG H L, HU X M. Effects of nitrogen treatment on rice yield and leaves conformation characteristics of 'Xindao No.18'., 2012, 14(3): 101-106. (in Chinese)

    [14] ZHAO K, TUNG C W, EIZENGA G C, WRIGHT M H, ALI M L, PRICE A H, NORTON G J, ISLAM M R, REYNOLDS A, MEZEY J, MCCLUNG A M, BUSTAMANTE C D, MCCOUCH S R. Genome-wide association mapping reveals a rich genetic architecture of complex traits in., 2011, 2(1): 1020-1021.

    [15] ATWELL S, HUANG YS, VILHJáLMSSON BJ, WILLEMS G, HORTON M, LI Y, MENG D, PLATT A, TARONE A M, HU T T, JIANG R, MULIYATI N W, ZHANG X, AMER M A, BAXTER I, BRACHI B, CHORY J, DEAN C, DEBIEU M, DE MEAUX J, ECKER JR, FAURE N, KNISKERN JM, JONES JD, MICHAEL T, NEMRI A, ROUX F, SALT DE, TANG C, TODESCO M, TRAW MB, WEIGEL D, MARJORAM P, BOREVITZ JO, BERGELSON J, NORDBORG M. Genome-wide association study of 107 phenotypes ininbred lines., 2010, 465(7298): 627-631.

    [16] HUANG X H, WEI X H, SANG T, ZHAO Q, FENG Q, ZHAO Y, LI C, ZHU C, LU T, ZHANG Z, LI M, FAN D, GUO Y, WANG A, WANG L, DENG L, LI W, LU Y, WENG Q, LIU K, HUANG T, ZHOU T, JING Y, LI W, LIN Z, BUCKLER ES, QIAN Q, ZHANG Q F, LI J, HAN B. Genome-wide association studies of 14 agronomic traits in rice landraces., 2010, 42(11): 961-967.

    [17] ZHOU Z K, JIANG Y, WANG Z, GOU Z, LYU J, LI W, YU Y, SHU L, ZHAO Y, MA Y, FANG C, SHEN Y, LIU T, LI C, LI Q, WU M, WANG M, WU Y, DONG Y, WAN W, WANG X, DING Z, GAO Y, XIANG H, ZHU B, LEE S H, WANG W, TIAN Z X. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean., 2015, 33(4): 408-414.

    [18] WEN W, LI D, LI X, GAO Y, LI W, LI H, LIU J, LIU H, CHEN W, LUO J, YAN J. Metabolome-based genome-wide association study of maize kernel leads to novel biochemical insights., 2015, 5(2): 487-507.

    [19] KUMAR V, SINGH A, MITHRA S V, KRISHNAMURTHY S L, PARIDA S K, JAIN S, TIWARI K K, KUMAR P, RAO A R, SHARMA S K, KHURANA J P, SINGH N K, MOHAPATRA T.Genome-wide association mapping of salinity tolerance in rice ()., 2015, 22(2): 133-145.

    [20] WANG C, YANG Y, YUAN X, XU Q, FENG Y, YU H, WANG Y, WEI X. Genome-wide association study of blast resistance in indica rice., 2014, 14(1): 1-11.

    [21] ZHU D, KANG H, LI Z, LIU M, ZHU X, WANG Y, WANG D, WANG Z, LIU W, WANG G L. A genome-wide association study of field resistance to magnaporthein rice., 2016, 9(1): 44.

    [22] CHEN W, GAO Y, XIE W, GONG L, LU K, WANG W, LI Y, LIU X, ZHANG H, DONG H, ZHANG W, ZHANG L, YU S, WANG G, LIAN X, LUO J. Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism., 2014, 46(7):714-721.

    [23] 王蘭, 黃李超, 代麗萍, 楊窯龍, 徐杰, 冷語佳, 張光恒, 胡江, 朱麗, 高振宇, 董國軍, 郭龍彪, 錢前, 曾大力. 利用日本晴/9311重組自交系群體定位水稻成熟期葉形相關(guān)性狀QTL. 中國水稻科學(xué), 2014, 28(6): 589-597.

    WANG L, HUANG L C, DAI L P, YANG Y L, XUE J, LENG Y J, ZHANG G H, HU J, GAO Z Y, DONG G J, GUO L B, QIAN Q, ZENG D L. QTL analysis for rice leaf morphology at maturity stage using a recombinant inbred line population derived from a cross between nipponbare and 9311., 2014, 28(6): 589-597. (in Chinese)

    [24] Murray M G, Thompson W F. Rapid isolation of high molecular weight plant DNA., 1980, 8(19): 4321-4325.

    [25] Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data., 2010, 38(16): e164.

    [26] Felsenstein J.PHYLIP-phylogeny inference package (version 3.2)., 1989, 5: 164-166.

    [27] Price A L, Patterson N J, Plenge R M, Weinblatt M E, Shadick N A, Reich D. Principal components analysis corrects for stratification in genome-wide association studies., 2006, 38(8): 904-909.

    [28] Kang H M, Sul J H, Service S K, Zaitlen N A, Kong S Y, Freimer N B, Sabatti C, Eskin E. Variance component model to account for sample structure in genome-wide association studies., 2010, 42(4): 348-354.

    [29] Dudbridge F, Gusnanto A. Estimation of significance thresholds for genomewide association scans., 2008, 32(3): 227-234.

    [30] MORI M, NOMURA T, OOKA H, ISHIZAKA M, YOKOTA T, SUGIMOTO K, OKABE K, KAJIWARA H, SATOH K, YAMAMOTO K, HIROCHIKA H, KIKUCHI S. Isolation and characterization of a rice dwarf mutant with a defect in brassinosteroid biosynthesis., 2002, 130(3): 1152-1161.

    [31] HONG Z, UEGUCHI-TANAKA M, SHIMIZU-SATO S, INUKAI Y, FUJIOKA S, SHIMADA Y, TAKATSUTO S, AGETSUMA M, YOSHIDA S, WATANABE Y, UOZU S, KITANO H, ASHIKARI M, MATSUOKA M. Loss-of-function of a rice brassinosteroid biosynthetic enzyme, C-6 oxidase, prevents the organized arrangement and polar elongation of cells in the leaves and stem., 2002, 32(4): 495-508.

    [32] LI M, XIONG G, LI R, CUI J, TANG D, ZHANG B, PAULY M, CHENG Z, ZHOU Y. Rice cellulose synthase-like D4 is essential for normal cell-wall biosynthesis and plant growth., 2009, 60(6):1055-1069.

    [33] HU J, ZHU L, ZENG D, GAO Z, GUO L, FANG Y, ZHANG G, DONG G, YAN M, LIU J, QIAN Q. Identification and characterization of NARROW AND ROLLED LEAF 1, a novel gene regulating leaf morphology and plant architecture in rice., 2010, 73(3): 283-292.

    [34] WADA S, HAYASHIDA Y, IZUMI M, KURUSU T, HANAMATA S, KANNO K, KOJIMA S, YAMAYA T, KUCHITSU K, MAKINO A, ISHIDA H. Autophagy supports biomass production and nitrogen use efficiency at the vegetative stage in rice., 2015, 168(1): 60-73.

    [35] KURUSU T, KOYANO T, HANAMATA S, KUBO T, NOGUCHI Y, YAGI C, NAGATA N, YAMAMOTO T, OHNISHI T, OKAZAKI Y, KITAHATA N, ANDO D, ISHIKAWA M, WADA S, MIYAO A, HIROCHIKA H, SHIMADA H, MAKINO A, SAITO K, ISHIDA H, KINOSHITA T, KURATA N, KUCHITSU K. OsATG7 is required for autophagy-dependent lipid metabolism in rice postmeiotic anther development., 2014, 10(5): 878-888.

    [36] LI Z K, PINSON S R, STANSEL J W, PATERSON A H. Genetic dissection of the source-sink relationship affecting fecundity and yield in rice (L.)., 1998, 4: 419-426.

    [37] YUE B, XUE W Y, LUO L J, XING Y Z. QTL Analysis for flag leaf characteristics and their relationships with yield and yield traits in rice.2006, 33: 824-832.

    [38] ZHU Y, CHANG L, TANG L, JIANG H, ZHANG W, CAO W. Modelling leaf shape dynamics in rice., 2009, 57(1): 73-81.

    [39] Dai G J, Cheng S H, Hua Z T, Zhang M L, Jiang H B, Feng Y, Shen X H, Su Y A, He N, Ma Z B, Ma X Q, Hou S G, Wang Y R. Mapping quantitative trait loci for nitrogen uptake and utilization efficiency in rice (L.) at different nitrogen fertilizer levels., 2015, 14(3): 10404-10414.

    (責(zé)任編輯 李莉)

    Genome-wide association analysis on flag leaf width under different nitrogenlevels in rice

    GAO YiHong, YAN JinXiang, TU ZhengJun, LENG YuJia, CHEN Long, HUANG LiChao, DAI LiPing, ZHANG GuangHeng, ZHU Li, HU Jiang, REN DeYong, GUO LongBiao, QIAN Qian, WANG DanYing, ZENG DaLi

    (China National Rice Research Institute/State Key Laboratory of Rice Biology, HangZhou 310006)

    【Objective】 The objective of this experiment is to study the genetic mechanism of flag leaf width and its response to different nitrogen fertilizer rates, and provide advantageous germplasm resources and genetic markers for the improvement of nitrogen use efficiency in rice breeding. 【Method】Based on re-sequencing of 134 rice landraces, a total of 3 356 591 SNPs distributed on the whole genome were identified. Three different nitrogen levels were assigned as the main plot in the split-plot design, and rice landraces are assigned at random to the subplots within each whole plot. Three nitrogen levels including low nitrogen (no nitrogen fertilizer), normal nitrogen (96 kg·hm-2) and high nitrogen (192 kg·hm-2) were applied under normal field cultivation, respectively. The EMMAX method was used to analyze the genetic relationship and EIGENSOFT was employed to detect the population structure. The mixed linear model was used to detect the potential genome-wide association between single nucleotide polymorphisms (SNPs) and the flag leaf width performance or response under low, medium and high nitrogen treatments. 【Result】The results showed that the flag leaf width displayed normal distribution in N0, N1 and N2 treatments, respectively. The variation of flag leaf width caused by the varietal differences and the different nitrogen levels, and the significant positive correlation between nitrogen fertilizer and flag leaf width were detected at all nitrogen levels. A total of 14 SNPs on five chromosomes presented significant association with the flag leaf width under three different nitrogen levels. The estimated minimum allele frequency was 0.46 at low nitrogen level, which indicate those loci are widely distributed in the association population. While the estimated minimum allele frequencies were lower at normal nitrogen and high nitrogen levels. One SNP located on chromosome 12 was found both at normal nitrogen and high nitrogen levels. Besides, another SNP on chromosome 12 was also detected at high nitrogen. Its flank contained a candidate gene, which belongs to the super-family of cytochrome P450. It is a homologous gene ofwhich was confirmed inregulating the flag leaf width in rice. Based on the variation of leaf width at different nitrogen levels, twenty SNPs and eight SNPs were identified to low nitrogen and high nitrogen response, respectively. Among them, the nitrogen utilization related gene,, was associated with high nitrogen response on chromosome 1. 【Conclusion】In this study, a total of 42 SNPs associated with flag leaf width and response to different nitrogen were identified based on the genome-wide association analysis.

    ; genome-wide association study; leaf width; nitrogen fertilizer

    2017-01-12;接受日期:2017-04-05

    國家自然科學(xué)基金(91435105、31661143006、31371581)、中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程項目

    高易宏,E-mail:gaoyihong93@163.com。燕金香,E-mail:2911669220@qq.com。高易宏和燕金香為同等貢獻作者。通信作者曾大力,E-mail:dalizeng@126.com。通信作者王丹英,E-mail:wangdanying@caas.cn

    猜你喜歡
    劍葉葉寬氮肥
    水稻回交群體劍葉性狀綜合評價及QTL定位
    氮肥供應(yīng)充足 春耕生產(chǎn)有保障
    馬奶子葡萄葉面積評估模型的建立
    回歸方程法測定兩個龍眼品種葉面積研究
    江淮小氮肥 耕耘六十年——紀(jì)念安徽小氮肥誕生六十周年
    抓住機遇 主動作為 努力推進我國氮肥市場穩(wěn)步前行
    基于骨架細(xì)化和坐標(biāo)點計算的水稻劍葉角測量
    基于支持向量機的水稻葉面積測定
    2017春季各地氮肥市場掠影
    洛陽市13種綠化樹種葉面積的回歸測算
    精品午夜福利视频在线观看一区| 变态另类成人亚洲欧美熟女| 国产精品香港三级国产av潘金莲| 91av网站免费观看| 老司机午夜十八禁免费视频| 亚洲午夜理论影院| 久久欧美精品欧美久久欧美| 国产午夜精品久久久久久| 午夜免费激情av| 国产激情久久老熟女| 欧美国产日韩亚洲一区| 91九色精品人成在线观看| 亚洲成人久久性| 深夜精品福利| 制服诱惑二区| 欧美一级a爱片免费观看看 | 国产麻豆成人av免费视频| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲av高清不卡| 日韩成人在线观看一区二区三区| 欧美日韩亚洲综合一区二区三区_| 99久久国产精品久久久| 亚洲av片天天在线观看| 日韩精品青青久久久久久| 九色成人免费人妻av| 亚洲成av人片免费观看| 50天的宝宝边吃奶边哭怎么回事| ponron亚洲| www.999成人在线观看| 一本久久中文字幕| 午夜精品一区二区三区免费看| 一边摸一边做爽爽视频免费| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 国产人伦9x9x在线观看| 亚洲中文字幕一区二区三区有码在线看 | 欧美又色又爽又黄视频| 极品教师在线免费播放| e午夜精品久久久久久久| 在线国产一区二区在线| 此物有八面人人有两片| 琪琪午夜伦伦电影理论片6080| 一本一本综合久久| 波多野结衣巨乳人妻| 欧美一区二区精品小视频在线| 久久精品影院6| 亚洲色图 男人天堂 中文字幕| 精品第一国产精品| 免费无遮挡裸体视频| 一区二区三区高清视频在线| 悠悠久久av| 人人妻人人澡欧美一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久国内视频| 国产亚洲av嫩草精品影院| 黄片大片在线免费观看| 一本一本综合久久| 成年免费大片在线观看| 国产精品久久电影中文字幕| 国产精品电影一区二区三区| 日韩欧美在线乱码| 777久久人妻少妇嫩草av网站| 国产午夜福利久久久久久| 国产午夜福利久久久久久| 村上凉子中文字幕在线| 欧美一区二区精品小视频在线| 亚洲精品在线观看二区| 亚洲精品在线美女| 日韩欧美精品v在线| 欧美日韩乱码在线| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡| 国内揄拍国产精品人妻在线| 男插女下体视频免费在线播放| 黄色a级毛片大全视频| 久久久久久免费高清国产稀缺| 在线免费观看的www视频| 高清毛片免费观看视频网站| 免费观看精品视频网站| 国产精品综合久久久久久久免费| 久久香蕉国产精品| 国产69精品久久久久777片 | 久久天堂一区二区三区四区| 麻豆国产av国片精品| 国产99久久九九免费精品| 亚洲成人久久性| 国产亚洲欧美在线一区二区| 非洲黑人性xxxx精品又粗又长| 91在线观看av| 两个人免费观看高清视频| 中国美女看黄片| 91大片在线观看| 国产一级毛片七仙女欲春2| 五月玫瑰六月丁香| e午夜精品久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 国产av又大| svipshipincom国产片| 日日夜夜操网爽| 国产伦一二天堂av在线观看| 变态另类成人亚洲欧美熟女| www.自偷自拍.com| 国产精品久久久av美女十八| АⅤ资源中文在线天堂| 啪啪无遮挡十八禁网站| 久久久久性生活片| 国产成人av激情在线播放| 亚洲一区中文字幕在线| 大型黄色视频在线免费观看| 一进一出好大好爽视频| 亚洲精品中文字幕一二三四区| 欧美日韩福利视频一区二区| 亚洲欧洲精品一区二区精品久久久| 欧美日本视频| 国产一级毛片七仙女欲春2| 麻豆国产97在线/欧美 | 999久久久国产精品视频| 亚洲人成77777在线视频| 首页视频小说图片口味搜索| 国产亚洲av高清不卡| 免费高清视频大片| 中出人妻视频一区二区| 欧美成人午夜精品| 亚洲狠狠婷婷综合久久图片| 美女扒开内裤让男人捅视频| 一进一出抽搐gif免费好疼| or卡值多少钱| 欧美国产日韩亚洲一区| 身体一侧抽搐| 亚洲精品中文字幕一二三四区| www.熟女人妻精品国产| www.自偷自拍.com| 国产伦人伦偷精品视频| 久久人妻福利社区极品人妻图片| 亚洲成av人片免费观看| 成人永久免费在线观看视频| 亚洲av美国av| 国产69精品久久久久777片 | 18禁国产床啪视频网站| 两人在一起打扑克的视频| 搡老熟女国产l中国老女人| 国内揄拍国产精品人妻在线| 国产精品久久视频播放| 一区二区三区国产精品乱码| 国产高清激情床上av| 青草久久国产| 亚洲中文日韩欧美视频| 丰满的人妻完整版| 国产真实乱freesex| 国产1区2区3区精品| 国产野战对白在线观看| 18美女黄网站色大片免费观看| 亚洲欧美日韩高清在线视频| 一级a爱片免费观看的视频| 一级片免费观看大全| 久久久精品大字幕| 丁香欧美五月| 99精品久久久久人妻精品| 91在线观看av| 99热只有精品国产| 久久久国产成人精品二区| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 亚洲成人免费电影在线观看| 久久精品综合一区二区三区| www.熟女人妻精品国产| 18禁观看日本| 久久亚洲真实| 亚洲av成人一区二区三| 丰满的人妻完整版| 最新在线观看一区二区三区| 身体一侧抽搐| 久久亚洲真实| 国产亚洲精品av在线| 国产高清激情床上av| 麻豆国产av国片精品| 久久精品亚洲精品国产色婷小说| 国产高清视频在线观看网站| 亚洲男人的天堂狠狠| 亚洲中文字幕一区二区三区有码在线看 | 狠狠狠狠99中文字幕| 精品久久久久久久久久久久久| 99久久久亚洲精品蜜臀av| 香蕉丝袜av| 国产人伦9x9x在线观看| 亚洲成人精品中文字幕电影| 一个人免费在线观看电影 | 免费一级毛片在线播放高清视频| 精品一区二区三区视频在线观看免费| 欧美在线黄色| 91成年电影在线观看| 亚洲色图 男人天堂 中文字幕| 在线播放国产精品三级| 欧美不卡视频在线免费观看 | 成年女人毛片免费观看观看9| 亚洲成av人片在线播放无| 伊人久久大香线蕉亚洲五| 一二三四社区在线视频社区8| 久久精品影院6| 久久久久国内视频| 狂野欧美白嫩少妇大欣赏| 成人18禁在线播放| 在线观看免费视频日本深夜| 亚洲激情在线av| 两性夫妻黄色片| 看免费av毛片| 免费av毛片视频| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 在线观看美女被高潮喷水网站 | 在线永久观看黄色视频| 狂野欧美激情性xxxx| 日韩中文字幕欧美一区二区| 国产高清有码在线观看视频 | 午夜免费观看网址| 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 亚洲精品色激情综合| 麻豆av在线久日| 最新在线观看一区二区三区| 男男h啪啪无遮挡| 老熟妇仑乱视频hdxx| 亚洲av成人av| 禁无遮挡网站| 国产高清videossex| 麻豆成人av在线观看| 午夜免费观看网址| 久久久久九九精品影院| 一本精品99久久精品77| 小说图片视频综合网站| 他把我摸到了高潮在线观看| 十八禁人妻一区二区| 精品第一国产精品| 亚洲精品中文字幕在线视频| 88av欧美| tocl精华| 一区二区三区国产精品乱码| 日韩欧美精品v在线| 日日干狠狠操夜夜爽| 色综合欧美亚洲国产小说| 亚洲欧美日韩东京热| ponron亚洲| 免费av毛片视频| 给我免费播放毛片高清在线观看| 怎么达到女性高潮| 人妻夜夜爽99麻豆av| 亚洲国产欧美人成| 免费在线观看成人毛片| 成人永久免费在线观看视频| 啦啦啦免费观看视频1| 亚洲精华国产精华精| 成人18禁高潮啪啪吃奶动态图| 精品国产亚洲在线| 视频区欧美日本亚洲| 久久久久性生活片| 国产免费男女视频| 色噜噜av男人的天堂激情| 国产精品久久视频播放| 日韩大尺度精品在线看网址| 久久欧美精品欧美久久欧美| 亚洲精华国产精华精| 中文字幕久久专区| 嫩草影视91久久| 日日干狠狠操夜夜爽| 男女床上黄色一级片免费看| 久久久久久久午夜电影| 久99久视频精品免费| 黄色a级毛片大全视频| 久久欧美精品欧美久久欧美| 天天添夜夜摸| 国产精品亚洲av一区麻豆| 久久 成人 亚洲| 在线观看66精品国产| 欧美黑人欧美精品刺激| 久久久国产成人免费| 精品久久久久久久人妻蜜臀av| 久久欧美精品欧美久久欧美| 亚洲精品中文字幕一二三四区| www.熟女人妻精品国产| 999久久久精品免费观看国产| 亚洲av日韩精品久久久久久密| av片东京热男人的天堂| 久久久久国产精品人妻aⅴ院| 国产av一区二区精品久久| 91国产中文字幕| 在线国产一区二区在线| 18禁黄网站禁片免费观看直播| 在线观看午夜福利视频| 黄色a级毛片大全视频| 女人高潮潮喷娇喘18禁视频| 久久亚洲真实| 成人一区二区视频在线观看| 亚洲精品一区av在线观看| 毛片女人毛片| 国产午夜福利久久久久久| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 我要搜黄色片| 无限看片的www在线观看| 视频区欧美日本亚洲| 怎么达到女性高潮| 天堂√8在线中文| 一本精品99久久精品77| 国产私拍福利视频在线观看| 日本三级黄在线观看| 精华霜和精华液先用哪个| 日本撒尿小便嘘嘘汇集6| 少妇粗大呻吟视频| av福利片在线| 999久久久国产精品视频| 久久午夜亚洲精品久久| 妹子高潮喷水视频| 国产精华一区二区三区| 亚洲 欧美 日韩 在线 免费| 在线观看免费日韩欧美大片| 日韩中文字幕欧美一区二区| 午夜福利视频1000在线观看| 亚洲精品在线美女| 后天国语完整版免费观看| 欧美成人午夜精品| 午夜福利视频1000在线观看| 动漫黄色视频在线观看| 十八禁人妻一区二区| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 国产伦一二天堂av在线观看| 我要搜黄色片| 99国产精品一区二区三区| 国产精品av久久久久免费| 成人精品一区二区免费| 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 国产精品一区二区三区四区久久| 淫妇啪啪啪对白视频| 国产精品免费一区二区三区在线| 国产久久久一区二区三区| 一a级毛片在线观看| 成人国语在线视频| 国产在线观看jvid| 久久久国产成人精品二区| 人人妻,人人澡人人爽秒播| 99久久国产精品久久久| 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 久久久久亚洲av毛片大全| bbb黄色大片| 熟女电影av网| 制服丝袜大香蕉在线| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 老司机午夜福利在线观看视频| 男人舔女人的私密视频| 亚洲黑人精品在线| 亚洲精品国产一区二区精华液| 亚洲色图av天堂| 十八禁网站免费在线| 少妇的丰满在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲色图av天堂| 国产午夜精品久久久久久| 成熟少妇高潮喷水视频| 国产一区二区三区视频了| 少妇裸体淫交视频免费看高清 | 国产爱豆传媒在线观看 | 白带黄色成豆腐渣| 视频区欧美日本亚洲| 一本精品99久久精品77| 人人妻人人澡欧美一区二区| 午夜免费激情av| 小说图片视频综合网站| 给我免费播放毛片高清在线观看| 12—13女人毛片做爰片一| 国产欧美日韩一区二区精品| 亚洲 国产 在线| 久久这里只有精品19| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| www.自偷自拍.com| 久久久久性生活片| 搡老岳熟女国产| 人妻久久中文字幕网| 日韩精品免费视频一区二区三区| 国产高清videossex| 国产精品一及| 国产99久久九九免费精品| 麻豆国产97在线/欧美 | 久久久久久人人人人人| 18禁黄网站禁片免费观看直播| av欧美777| 欧美日韩福利视频一区二区| 老司机深夜福利视频在线观看| 香蕉av资源在线| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 夜夜爽天天搞| 欧美乱码精品一区二区三区| 黄片小视频在线播放| 成年女人毛片免费观看观看9| 一进一出抽搐gif免费好疼| 18美女黄网站色大片免费观看| 香蕉久久夜色| www.999成人在线观看| 欧洲精品卡2卡3卡4卡5卡区| 最好的美女福利视频网| 制服诱惑二区| 亚洲专区国产一区二区| 日本免费a在线| 午夜日韩欧美国产| 亚洲18禁久久av| 色尼玛亚洲综合影院| 亚洲美女黄片视频| 国产视频内射| 亚洲国产看品久久| 久久国产精品人妻蜜桃| 19禁男女啪啪无遮挡网站| 亚洲中文日韩欧美视频| 亚洲精品中文字幕一二三四区| 免费无遮挡裸体视频| 又大又爽又粗| 国产爱豆传媒在线观看 | 日韩免费av在线播放| 99久久国产精品久久久| 美女免费视频网站| 国内久久婷婷六月综合欲色啪| 久久中文看片网| 国产99久久九九免费精品| 91成年电影在线观看| 无限看片的www在线观看| 日韩av在线大香蕉| 亚洲精品在线美女| 精品无人区乱码1区二区| 国产精品影院久久| 欧美极品一区二区三区四区| 久久天堂一区二区三区四区| 国产精品电影一区二区三区| 我要搜黄色片| 又紧又爽又黄一区二区| 麻豆成人午夜福利视频| 亚洲av美国av| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| 国产精品久久久久久久电影 | 丝袜人妻中文字幕| 国内少妇人妻偷人精品xxx网站 | 久久久久国产精品人妻aⅴ院| 久久精品91无色码中文字幕| 亚洲精品中文字幕在线视频| 桃色一区二区三区在线观看| 男女做爰动态图高潮gif福利片| 男女那种视频在线观看| 91成年电影在线观看| 天天一区二区日本电影三级| 免费在线观看黄色视频的| 男人舔奶头视频| 99久久综合精品五月天人人| 国产精品一及| 中文字幕av在线有码专区| 性色av乱码一区二区三区2| 亚洲欧美精品综合一区二区三区| 最近在线观看免费完整版| av中文乱码字幕在线| 亚洲avbb在线观看| 欧美黑人欧美精品刺激| 韩国av一区二区三区四区| 99久久无色码亚洲精品果冻| 我的老师免费观看完整版| 国产日本99.免费观看| 国产不卡一卡二| 国产精品一区二区免费欧美| 久久久精品欧美日韩精品| 激情在线观看视频在线高清| 制服诱惑二区| 久久国产乱子伦精品免费另类| 色尼玛亚洲综合影院| 久久久水蜜桃国产精品网| 亚洲专区字幕在线| 99riav亚洲国产免费| 亚洲 国产 在线| 91麻豆av在线| 亚洲18禁久久av| 女人爽到高潮嗷嗷叫在线视频| 亚洲精华国产精华精| 久久久久久人人人人人| 欧美成人免费av一区二区三区| 日韩欧美三级三区| 好看av亚洲va欧美ⅴa在| 久久精品91无色码中文字幕| 伊人久久大香线蕉亚洲五| 人妻久久中文字幕网| 日本精品一区二区三区蜜桃| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| 嫩草影视91久久| 熟女少妇亚洲综合色aaa.| 亚洲国产欧洲综合997久久,| 国产免费男女视频| 久久九九热精品免费| 亚洲精品一区av在线观看| 色尼玛亚洲综合影院| 亚洲国产欧美人成| 白带黄色成豆腐渣| 在线观看日韩欧美| 国产精品香港三级国产av潘金莲| 欧美最黄视频在线播放免费| 亚洲av熟女| 999精品在线视频| 国产成人欧美在线观看| 亚洲男人天堂网一区| 亚洲av片天天在线观看| 1024视频免费在线观看| 波多野结衣高清无吗| svipshipincom国产片| 亚洲中文av在线| 曰老女人黄片| 亚洲人成电影免费在线| 99re在线观看精品视频| 成人一区二区视频在线观看| cao死你这个sao货| 18禁观看日本| 国产精品永久免费网站| or卡值多少钱| 国产视频一区二区在线看| 免费在线观看成人毛片| 国产精品久久久久久精品电影| 中文字幕av在线有码专区| 国产精品 国内视频| 久久精品人妻少妇| 高清在线国产一区| 听说在线观看完整版免费高清| 在线观看免费视频日本深夜| 夜夜夜夜夜久久久久| 久久中文看片网| 国产成+人综合+亚洲专区| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人不卡在线观看播放网| 国产片内射在线| 久久久久久大精品| 国产亚洲av高清不卡| av片东京热男人的天堂| 国产精品,欧美在线| 久久精品国产综合久久久| 亚洲中文av在线| 日韩欧美国产在线观看| 国产精品,欧美在线| 国产激情欧美一区二区| 日韩欧美三级三区| 日本三级黄在线观看| 国产欧美日韩一区二区三| 精品第一国产精品| 国产精品一区二区三区四区久久| 成人一区二区视频在线观看| 视频区欧美日本亚洲| 香蕉久久夜色| 亚洲aⅴ乱码一区二区在线播放 | 夜夜爽天天搞| 国产视频一区二区在线看| 好男人在线观看高清免费视频| 午夜福利高清视频| av片东京热男人的天堂| 美女免费视频网站| 亚洲第一电影网av| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合| 久久中文字幕一级| 亚洲一区中文字幕在线| 十八禁人妻一区二区| 日韩欧美 国产精品| 午夜激情福利司机影院| 国产成人系列免费观看| 高清在线国产一区| 老汉色av国产亚洲站长工具| 韩国av一区二区三区四区| 老鸭窝网址在线观看| 亚洲 欧美 日韩 在线 免费| 欧美一级毛片孕妇| 欧美一区二区国产精品久久精品 | 美女黄网站色视频| 窝窝影院91人妻| 国产伦人伦偷精品视频| 国产69精品久久久久777片 | 波多野结衣高清无吗| 蜜桃久久精品国产亚洲av| 老汉色av国产亚洲站长工具| 欧美 亚洲 国产 日韩一| 国产成人影院久久av| 好男人电影高清在线观看| 欧美极品一区二区三区四区| 波多野结衣高清作品| 亚洲国产日韩欧美精品在线观看 | 国产又色又爽无遮挡免费看| 国产高清videossex| 黄色视频不卡| 老司机深夜福利视频在线观看| www.www免费av| 国产精品综合久久久久久久免费| 性欧美人与动物交配| 色综合亚洲欧美另类图片| 丁香六月欧美| 亚洲av熟女| 12—13女人毛片做爰片一| 色综合欧美亚洲国产小说| 久久精品影院6| 欧美精品啪啪一区二区三区| 一本一本综合久久| av中文乱码字幕在线| 成人三级做爰电影| 免费高清视频大片| 亚洲国产精品合色在线| 99在线视频只有这里精品首页| 精品久久久久久,| 久久午夜综合久久蜜桃| 亚洲精品国产精品久久久不卡| 国产午夜精品论理片|