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

    馬鈴薯種質(zhì)的遺傳特性與抗旱性的關(guān)系

    2016-08-02 09:38:13白江平楊宏羽李維婷高慧娟張俊蓮張金林王旺田
    草業(yè)科學(xué) 2016年3期
    關(guān)鍵詞:遺傳多樣性抗旱性

    婁 艷,白江平,楊宏羽,李維婷,高慧娟,張俊蓮,王 蒂,張金林,王旺田

    (1.甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室/甘肅省干旱生境作物學(xué)重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730070;2.甘肅農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,甘肅 蘭州 730070; 3.甘肅農(nóng)業(yè)大學(xué)園藝學(xué)院,甘肅 蘭州 730070;4.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020)

    ?

    馬鈴薯種質(zhì)的遺傳特性與抗旱性的關(guān)系

    婁 艷1,2,白江平1,2,楊宏羽1,3,李維婷1,2,高慧娟1,2,張俊蓮1,3,王 蒂1,2,張金林4,王旺田2

    (1.甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室/甘肅省干旱生境作物學(xué)重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730070;2.甘肅農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,甘肅 蘭州 730070; 3.甘肅農(nóng)業(yè)大學(xué)園藝學(xué)院,甘肅 蘭州 730070;4.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020)

    摘要:為探究馬鈴薯(Solanum tuberosum)不同遺傳特性與其對(duì)水分脅迫響應(yīng)之間的相關(guān)關(guān)系,本研究通過(guò)SSR分析了供試材料的遺傳多樣性,并通過(guò)PEG-6000和不澆水脅迫處理馬鈴薯移栽苗的方法,研究了干旱脅迫下不同馬鈴薯品種中過(guò)氧化物酶(POD、CAT)活性及可溶性蛋白含量的變化。結(jié)果顯示,在遺傳相似系數(shù)(GS)為0.700水平上,44份馬鈴薯材料可被聚為5類,大部分品種聚在一起,說(shuō)明供試品種的遺傳基礎(chǔ)相對(duì)狹窄。通過(guò)GGE-Biplot綜合分析POD、CAT及可溶性蛋白質(zhì)變化量,將44份供試材料分為6組。比對(duì)分析遺傳多樣性分析與生理結(jié)果,發(fā)現(xiàn)二者之間并不具有明顯的一致性,即品種遺傳特性的差異并不能直接反映品種的抗旱性能差異,而兩種脅迫方式在研究品種抗旱特性上具有較高的一致性。該研究結(jié)果可以為今后抗旱試驗(yàn)的進(jìn)行提供參考。

    關(guān)鍵詞:遺傳多樣性;PEG脅迫;雙標(biāo)圖;抗旱性

    自1995年Veilleux[1]首次使用SSR標(biāo)記研究馬鈴薯(Solanumtuberosum)以來(lái),SSR標(biāo)記因其多態(tài)性高、可重復(fù)性好、操作簡(jiǎn)單等優(yōu)點(diǎn)而被廣泛應(yīng)用于馬鈴薯遺傳關(guān)系分析研究[2-4]、品種鑒定[5]、重點(diǎn)種質(zhì)資源的收集和評(píng)價(jià)[6]等方面。近年來(lái),中國(guó)學(xué)者在應(yīng)用SSR對(duì)種質(zhì)材料遺傳多樣性分析[7-10]、作物抗旱性相關(guān)基因的SSR分子標(biāo)記[11]、水分脅迫下抗旱相關(guān)性狀的QTL分析[12]及對(duì)控制水分脅迫應(yīng)答蛋白的相關(guān)基因的SSR標(biāo)記[13]等方面也做了大量的研究。

    干旱是影響植物正常生長(zhǎng)發(fā)育的主要環(huán)境脅迫因素之一[14]。植物在長(zhǎng)期進(jìn)化過(guò)程中,演化出一系列形態(tài)學(xué)的和生理生化方面的適應(yīng)性機(jī)制和策略,能有效抵制脅迫逆境如調(diào)節(jié)保護(hù)酶系統(tǒng)[15-16],改變過(guò)氧化物酶(POD)、過(guò)氧化氫酶(CAT) 的活性等,在逆境脅迫下動(dòng)員酶和非酶防御系統(tǒng)來(lái)清除活性氧,保護(hù)植物細(xì)胞免受氧化損害。自1979 年,Heyser和Nabors[17]首次以PEG作為誘導(dǎo)劑和篩選劑篩選出抗旱的煙草(Nicotianatabacum)細(xì)胞系以來(lái),高分子量的聚乙二醇(PEG)(分子量≥6 000)常作為一種非穿透滲透劑被廣泛用于植物研究領(lǐng)域[18-21]。然而,關(guān)于不同馬鈴薯品種在PEG-6000溶液處理下進(jìn)行模擬干旱脅迫與自然干旱水分脅迫處理方面的研究很少。同時(shí),眾所周知,作物在水分脅迫條件下的抗性因品種遺傳特性的不同而異,而就能否從遺傳多樣性的角度來(lái)解釋不同種質(zhì)間脅迫響應(yīng)差異方面的研究還未見報(bào)道。

    本研究通過(guò)SSR分子標(biāo)記對(duì)44份馬鈴薯材料的遺傳多樣性進(jìn)行分析,并利用GGE-Biplot對(duì)44份馬鈴薯材料在水分逆境脅迫下的響應(yīng)能力進(jìn)行綜合評(píng)價(jià)分析,嘗試從遺傳多樣性的角度解釋不同種質(zhì)間脅迫響應(yīng)差異,并探索兩種水分脅迫方式在研究作物抗旱性能上是否具有一致性,以期為馬鈴薯抗旱性材料的篩選和鑒定,品種特性的改良等方面提供理論依據(jù)。

    1材料與方法

    1.1供試材料

    試驗(yàn)于2014年6-9月在甘肅農(nóng)業(yè)大學(xué)作物遺傳育種和干旱生境作物學(xué)重點(diǎn)實(shí)驗(yàn)室進(jìn)行,本研究選用42份從國(guó)際馬鈴薯研究中心(CIP)引進(jìn)的馬鈴薯材料(表1)和兩份本地品種,即抗逆性較強(qiáng)的栽培種隴薯3號(hào)(L3)和對(duì)逆境脅迫較敏感的栽培種大西洋(DXY)脫毒試管苗為試驗(yàn)材料。所有品種均由甘肅農(nóng)業(yè)大學(xué)馬鈴薯品種改良重點(diǎn)實(shí)驗(yàn)室提供。

    1.2試驗(yàn)方法

    1.2.1馬鈴薯材料試管苗的培養(yǎng)以同一時(shí)期轉(zhuǎn)接的長(zhǎng)勢(shì)相對(duì)一致的脫毒試管苗為基礎(chǔ)苗,將其剪成至少帶一片葉的小莖段,插入裝有固體培養(yǎng)基(MS)50 mL的三角瓶中,每瓶扦插6~8個(gè)莖段。置于平均溫度為(22±2) ℃,光照為3 000 lx,光周期為16 h·d-1的培養(yǎng)室進(jìn)行培養(yǎng)。

    1.2.2基因組DNA的提取與檢測(cè)及試管苗的移栽1)剪取培養(yǎng)30 d左右的各供試材料試管苗新鮮莖葉,用液氮研磨并參考 CTAB[22]法提取總基因組DNA。采用1%瓊脂糖凝膠電泳檢測(cè)DNA質(zhì)量,用紫外分光光度計(jì)法檢測(cè)其質(zhì)量和濃度。用ddH2O將各樣品稀釋到25 ng·μL-1,置于-20 ℃冰箱保存。2)將經(jīng)過(guò)煉苗的試管苗洗去培養(yǎng)基,剪去根部(保留一部分),移栽入裝有蛭石(基質(zhì))的72孔育苗盤中(皿數(shù)6×12,口徑40 mm×40 mm, 底部20 mm×20 mm,皿深45 mm),種植深3~4 cm,每個(gè)品種每個(gè)處理各種12株,移栽在不同的穴盤上,共36株。扦插后,澆透水,將育苗盤擺放在培養(yǎng)室,蓋上塑料薄膜保濕5~7 d。每隔1 d澆1次水(1周澆3次水),每隔1周澆1次營(yíng)養(yǎng)液(即不加蔗糖的液體MS培養(yǎng)基),營(yíng)養(yǎng)液pH值在5.8左右。移栽苗在育苗盤生長(zhǎng)20 d后進(jìn)行試驗(yàn)處理。PEG處理組用10% PEG(-0.120 MPa)澆灌(第1次10%PEG10mL·穴-1,其余幾次均為10%PEG5mL·穴-1);自然脅迫處理組不澆水;對(duì)照(CK)正常澆水,水量和處理組PEG的量一致。處理6 d后,最早發(fā)現(xiàn)2號(hào)材料PEG處理組的植株出現(xiàn)萎蔫癥狀,拍照、采樣后分裝并冷凍于-80 ℃冰箱備用。

    1.2.3SSR條帶檢測(cè)與相關(guān)生理指標(biāo)的測(cè)定參照已發(fā)表的24對(duì)SSR引物[5,8,23]并由上海生工生物工程技術(shù)服務(wù)有限公司合成(23對(duì)表2已列出),進(jìn)行SSR-PCR擴(kuò)增,PCR產(chǎn)物用8.0%非變性聚丙烯酰胺凝膠電泳分離檢測(cè),銀染顯色并統(tǒng)計(jì)多態(tài)性條帶。

    表1 42份引進(jìn)馬鈴薯種質(zhì)材料編號(hào)及親本材料

    過(guò)氧化氫酶(CAT)的活性測(cè)定采用紫外吸收法[24];過(guò)氧化物酶(POD)活性測(cè)定采用愈創(chuàng)木酚法[24],可溶性蛋白質(zhì)的提取及含量測(cè)定采用考馬斯亮藍(lán)G-250法[24]。

    1.2.4條帶統(tǒng)計(jì)與數(shù)據(jù)分析采用人工讀膠法統(tǒng)計(jì)電泳譜帶,在同一位置上清晰的條帶記為1,無(wú)條帶記為0,擴(kuò)增不出的記為-9,建立0,1數(shù)據(jù)庫(kù)并保存到Excel中。利用 NTSYS-pc 2.10e軟件[25],通過(guò)UPGMA方法進(jìn)行聚類分析。通過(guò)Popgene1.32軟件[26]計(jì)算每個(gè)位點(diǎn)的等位變異數(shù)目,有效等位基因數(shù)目[27],Shannon’s信息指數(shù)[28]及等位基因頻率;計(jì)算多態(tài)性信息含量(Polymorphism Information Content):PIC=1-∑(Pi2)[29]。其中Pi指在第i位點(diǎn)處檢測(cè)到的等位基因頻率。

    用Microsoft Excel 2007對(duì)生理數(shù)據(jù)進(jìn)行輸入、整理,并計(jì)算出對(duì)應(yīng)品種脅迫前后酶活性的變化量,相對(duì)變化量=(脅迫后植株體內(nèi)酶含量-處理后對(duì)照組中酶含量)/對(duì)照中酶含量×100%;通過(guò)SPSS Statistics 18.0 軟件進(jìn)行SPSS方差分析(Analysis of Variance ANOVA)(P≤0.05)和聚類分析,多重比較采用最小顯著差數(shù)法(LSDα)和新復(fù)極差法(Duncan法),用來(lái)檢測(cè)不同品種脅迫前后酶活性變化量差異。檢測(cè)到顯著差異P≤0.05后,使用GGE-Biplot軟件[23,30]對(duì)兩向數(shù)據(jù)進(jìn)行主成分分析。

    2結(jié)果與分析

    2.1SSR標(biāo)記多態(tài)性分析

    本研究用24對(duì)SSR引物對(duì)44份馬鈴薯供試材料進(jìn)行PCR,選出其中多態(tài)性較好的23對(duì)SSR引物(表2)。結(jié)果顯示,23對(duì)SSR引物在44份供試材料中共檢測(cè)到66個(gè)SSR標(biāo)記,平均每個(gè)等位位點(diǎn)的等位變異數(shù)為2.869 6,變異范圍為2~5;平均等位變異數(shù)為2.463 1,平均有效等位變異所占的比重為72.88%;23對(duì)SSR引物的Shannon’s信息指數(shù)的變化范圍為0.248 9~1.334 8,平均為0.776 2;多態(tài)性信息含量 (PIC)變化范圍為0.127 1~0.727 1,平均為0.465 0,表明供試材料的遺傳多樣性較高(表2)。

    2.2遺傳相似系數(shù)和聚類分析

    為了分析供試材料的遺傳相似系數(shù)(GS),并以更直觀的方式揭示材料間的親緣關(guān)系。本研究利用NTSYS-pc2.10e軟件分析了供試馬鈴薯材料的遺傳相似系數(shù)(GS),并根據(jù)UPGMA方法進(jìn)行聚類分析。結(jié)果顯示,44份供試馬鈴薯材料的遺傳相似系數(shù)變異范圍為0.484 9~0.909 1。其中31和36號(hào)材料之間的遺傳相似系數(shù)最大(GS=0.909 1),表明這兩份材料的遺傳背景相似;其次是材料8與9 和28與29(GS=0.893 9),所以不建議用材料31與36、8與9和28與29進(jìn)行親本雜交組合配置;遺傳相似系數(shù)最小的是材料4與14(GS=0.484 9),說(shuō)明材料4與14之間的親緣關(guān)系較遠(yuǎn),遺傳差異性較大。UPGMA方法聚類分析發(fā)現(xiàn)(圖1),在遺傳相似系數(shù)(GS)為0.700水平上,44份馬鈴薯材料可被聚為5類,分別記為A、B、C、D和E,每類分別包含有6、32、4、1和1份馬鈴薯材料。B類中所包含的32份馬鈴薯材料,占所有供試材料的72.7%。

    表2 23對(duì)SSR引物在44份馬鈴薯材料中的多態(tài)性

    注:SSR引物選自已發(fā)表的論文研究:a[5],b[23], c[8]。

    Note:SSR primer sequences were selected from previous studies: a[5], b[23], c[8].

    2.3水分脅迫處理相關(guān)數(shù)據(jù)分析

    為了分析兩種處理下各不同變量之間的相關(guān)性,判斷分析兩種水分脅迫方式之間是否有一定的一致性。采用SPSS R型聚類(變量聚類),對(duì)PEG處理下可溶性蛋白相對(duì)于對(duì)照的變幅(PEG-P或PEG protein)、 PEG處理下POD相對(duì)于對(duì)照的變幅(PEG-POD)、PEG處理下CAT相對(duì)于對(duì)照的變幅(PEG-CAT)、自然脅迫下可溶性蛋白相對(duì)于對(duì)照的變幅(ZR-P或ZR protein)、自然脅迫下POD相對(duì)于對(duì)照的變幅(ZR-POD)、自然脅迫下CAT相對(duì)于對(duì)照的變幅(ZR-CAT)共6個(gè)變量進(jìn)行相關(guān)性分析(表3),結(jié)果顯示,ZR-POD與PEG-POD之間顯著相關(guān),相關(guān)系數(shù)最大,為0.965,ZR-CAT與PEG-CAT之間顯著相關(guān),相關(guān)系數(shù)為0.933,ZR protein與PEG protein間顯著相關(guān),相關(guān)系數(shù)為0.888,均在P<0.01水平上達(dá)到極顯著。不同脅迫處理下,同種酶酶活性變量間相關(guān)性較高,說(shuō)明不論是PEG模擬水分脅迫還是自然水分脅迫,對(duì)馬鈴薯相關(guān)生理生化指標(biāo)的影響具有較高的一致性。材料1因在前期試驗(yàn)中死亡,故后續(xù)生理測(cè)定試驗(yàn)中只測(cè)定了43份材料。

    為了對(duì)供試馬鈴薯材料進(jìn)行生理上群體劃分并解釋各處理變量之間的關(guān)系,本研究利用GGE-Biplot對(duì)43個(gè)品種不同處理下不同酶活性變化量進(jìn)行了分析(圖2)。圖2中的多邊形由連接同一方向上距雙標(biāo)圖原點(diǎn)最遠(yuǎn)的品種的酶活性變化量位點(diǎn)構(gòu)成,其余的品種都被包括在多邊形內(nèi)部。由原點(diǎn)出發(fā)并且與多邊形各邊垂直的射線將整個(gè)雙標(biāo)圖分為6個(gè)大小不同的區(qū)域,并由此將不同品種分為6個(gè)組,即a、b、c、d、e和f。

    圖1 供試馬鈴薯材料聚類圖

    變量Variable變量VariablePEGproteinZRproteinPEG-PODZR-PODPEG-CATZR-CATPEGprotein1.0001.000ZRprotein0.888**1.000PEG-POD0.4590.4761.000ZR-POD0.5010.5230.965**1.000PEG-CAT0.5590.5500.5380.5171.000ZR-CAT0.4960.5300.5210.4910.933**1.000

    注:**代表在0,01水平上顯著相關(guān)。

    Note: ** show significant difference at 0.01 level.

    圖2 GGE-Biplot雙標(biāo)圖分析不同馬鈴薯品種與不同脅迫處理方式及脅迫前后不同酶活性變化間的關(guān)系

    處理變量可以被分為兩個(gè)組:PEGP(PEG protein)、ZRP(ZR protein)、PEG-CAT及ZR-CAT為一組,PEG-POD和ZR-POD為一組,而且第一主成分(PC1)與第二主成分(PC2)一起可以解釋84.1%的變量,說(shuō)明雙標(biāo)圖足以解釋不同處理下不同酶活性變化量與不同品種間的差異。不論是在PEG處理下,還是在自然水分脅迫下,可溶性蛋白含量及CAT活性在品種38中變化量最大, 而POD在品種35中的變化量最大,說(shuō)明相比于其它品種,35和38的抗氧化系統(tǒng)對(duì)干旱脅迫的響應(yīng)更為靈敏。為了分析不同供試材料在脅迫前后酶活性變化量的大小及同一品種在不同酶活性測(cè)定下,其變化量的穩(wěn)定性,本研究利用GGE-Biplot對(duì)43份供試馬鈴薯材料在水分脅迫誘導(dǎo)前后酶活性的變化量進(jìn)行排序,并對(duì)其變化量的穩(wěn)定性進(jìn)行分析。圖中帶箭頭的直線是平均處理軸,它同時(shí)經(jīng)過(guò)雙標(biāo)圖原點(diǎn)和平均處理值(紅色圓圈),其位置取決于各處理坐標(biāo)的平均值[31],箭頭所指的方向是酶活性變化量在所有處理下的近似平均變化量,箭頭方向代表變化量越大。通過(guò)雙標(biāo)圖原點(diǎn)并與平均處理軸垂直的雙箭頭直線代表各品種酶活性變化量的穩(wěn)定性,箭頭所指的方向?yàn)樽兓环€(wěn)定,并且越偏離平均處理軸,說(shuō)明變化量越不穩(wěn)定。如圖3所示,品種23、13、20、25、43、35、18和38的變化量均高于平均值,而且品種38的變化量最大。但是,13、20、35、38這4個(gè)品種的位點(diǎn)與平均處理軸間的距離較遠(yuǎn),所以它們的變化易受不同酶活性測(cè)定的影響,穩(wěn)定性較差。為進(jìn)一步分析不同處理間及不同酶活性測(cè)定間的關(guān)系,利用GGE-Biplot進(jìn)行區(qū)別力與代表性的分析。如圖4所示,連接各處理變量的位點(diǎn)與雙標(biāo)圖原點(diǎn)的直線稱為處理向量。向量的不同長(zhǎng)度反映的是不同處理變量對(duì)酶活性變化量的區(qū)別能力,向量之間的夾角的余弦值近似代表兩個(gè)處理變量之間的相關(guān)性系數(shù)[30]。除PEG-POD與ZR-POD外,各處理的向量長(zhǎng)度比較接近,說(shuō)明同一酶活性變化量在不同處理變量下的變化趨勢(shì)較接近。PEG-P與ZR-P、PEG-CAT與ZR-CAT、PEG-POD與ZR-POD之間的夾角均很小,說(shuō)明PEG脅迫處理和自然脅迫處理之間具有較高的相關(guān)性,這與通過(guò)SPSS R型聚類對(duì)各處理變量間的相關(guān)性分析結(jié)果一致。而兩種脅迫處理下POD測(cè)定與可溶性蛋白、CAT測(cè)定間的夾角較大,雖然小于90°,但仍能說(shuō)明,即使在相同的脅迫處理下POD活性變化量與可溶性蛋白、CAT活性變化量之間的變化差異較為明顯,尤其與可溶性蛋白含量的變化更為明顯。

    3討論與結(jié)論

    通過(guò)對(duì)種質(zhì)材料遺傳相似系數(shù)的研究分析發(fā)現(xiàn),供試材料遺傳背景的相似性與其遺傳相似系數(shù)大小之間具有一定的相關(guān)性。如本研究中材料28與29具有相同的父本和母本,材料33與34則只具有相同的父本,育成品種親緣關(guān)系亦比較近,這與材料28與29(GS=0.893 9)和33與34(GS=0.878 8)較高的遺傳相似系數(shù)是一致的,這也與劉建霞等[7-9]的相關(guān)研究中遺傳背景相似的品種遺傳差異較小的結(jié)論是一致的。然而,31和36號(hào)材料的親本雖然完全不同,但二者之間的遺傳相似系數(shù)最大(GS=0.909 1),仍表明這兩份材料的親緣關(guān)系是較近的,這可能是由于材料本身遺傳特性的差異或發(fā)生遺傳變異引起的,也有可能是由于SSR分子標(biāo)記試驗(yàn)所用的引物較少,未能充分表現(xiàn)材料本身的遺傳特異性等原因引起的。

    圖3 利用GGE-Biplot分析43份馬鈴薯材料在兩種脅迫處理下酶活性變化量的排序

    圖4 利用GGE-Biplot雙標(biāo)圖分析各處理之間的相關(guān)性

    馬鈴薯根系較淺,對(duì)水分虧缺非常敏感[32]。水分虧缺會(huì)對(duì)植物的生長(zhǎng)狀況、形態(tài)結(jié)構(gòu)與生理生化產(chǎn)生顯著影響[33-34]。本研究發(fā)現(xiàn),受到水分脅迫后,不同馬鈴薯植株體內(nèi)POD、CAT活性及可溶性蛋白含量較對(duì)照均有所上升,這與大量研究報(bào)道[35-37]一致。且?guī)追N酶活性變化在兩種不同的脅迫處理間具有較高的相關(guān)性,這與抗艷紅等[34]、張立軍等[38]的PEG模擬干旱試驗(yàn)測(cè)得的生理指標(biāo)與室外旱池的實(shí)際抗旱指標(biāo)有較好的一致關(guān)系的結(jié)論相符合。

    綜合比較SSR分子標(biāo)記聚類分析與生理分析結(jié)果,發(fā)現(xiàn)兩種聚類分析結(jié)果間并不具有較高一致性。例如SSR聚類分析B類中的材料31和36間因具有較高的遺傳相似系數(shù)(GS=0.909 1)而被聚在一起,與此同時(shí),二者在受到水分脅迫后在生理生化反應(yīng)也具有一定的一致性,而被聚為同一組(e)中。然而,對(duì)于大部分品種來(lái)說(shuō),遺傳親緣關(guān)系的遠(yuǎn)近與其對(duì)生理脅迫的響應(yīng)間并不具有高度的一致性。如材料28、29具有相同的父本[392 820.1=(C93.154)]和母本(CHIEFTAIN),遺傳背景相似而聚在一起,然而,二者在對(duì)水分脅迫的響應(yīng)并不一致,分別屬于不同的組,即d和e。同樣,材料8、9間因具有較高的遺傳相似系數(shù)(GS=0.893 9)而被聚在一起,但二者在受到水分脅迫后其生理生化反應(yīng)并不具有較高一致性,二者分別屬于不同的組,即b和c。綜合多個(gè)品種遺特性及對(duì)生理生化響應(yīng)的差異,發(fā)現(xiàn)二者之間并不具有簡(jiǎn)單的相關(guān)性。其主要原因可能為:1)材料多樣性分析的出發(fā)點(diǎn)不一樣。SSR分子標(biāo)記通過(guò)分析,且廣泛分布在真核生物基因組上的,由于其重復(fù)次數(shù)不同或重復(fù)程度不完全而形成每個(gè)座位的多態(tài)性,對(duì)物種種質(zhì)資源進(jìn)行遺傳多樣性分析[39]。然而,GGE-Biplot聚類分析是通過(guò)檢測(cè)植物受到水分逆境脅迫后,其體內(nèi)抗氧化物酶活性的變化情況,對(duì)供試材料的水分逆境脅迫響應(yīng)快慢進(jìn)行的分析。2)兩種試驗(yàn)方法采用了不同的統(tǒng)計(jì)分析軟件,數(shù)據(jù)算法不一樣,如SSR分子標(biāo)記使用NTSYS-pc2.10e,采用UPGMA法構(gòu)建聚類圖,而水分脅迫后,酶活性變化的分析采用SPSS R型聚類,因此也使得聚類結(jié)果有差異。3)SSR分子標(biāo)記試驗(yàn)所用的引物較少,因此擴(kuò)增獲得的多態(tài)性條帶數(shù)目也有限,未能充分表現(xiàn)材料本身的遺傳特異性。另一方面,植物對(duì)水分脅迫逆境的響應(yīng)機(jī)制是由多種因素共同參與的復(fù)雜過(guò)程[7,40],單純的對(duì)水分脅迫逆境下對(duì)抗氧化物酶活性的測(cè)定并不能充分表現(xiàn)品種的生理特性。

    本研究通過(guò)SSR分子標(biāo)記技術(shù)將44份馬鈴薯材料在遺傳相似系數(shù)(GS)為0.700水平上被聚為5類,且大部分品種被聚在一個(gè)類群,說(shuō)明供試品種遺傳基礎(chǔ)還是相對(duì)狹窄。遺傳相似系數(shù)較大的品種往往被聚在一起,說(shuō)明它們的親緣關(guān)系較近。通過(guò)水分脅迫試驗(yàn)發(fā)現(xiàn),馬鈴薯移栽苗受到PEG-6000澆灌和自然不澆水兩種脅迫處理后其體內(nèi)的POD、CAT及可溶性蛋白含量相對(duì)于對(duì)照均上升,且同一種酶在兩種脅迫處理下的變化量具有較高的相關(guān)性,說(shuō)明兩種水分脅迫方式在鑒定馬鈴薯抗旱特性上具有較高的一致性。綜合分析3種酶在兩種脅迫處理下,不同品種中的變化量,將43份材料聚為大小不一的6組。綜合比較兩種聚類結(jié)果發(fā)現(xiàn),二者并不具有一致的相關(guān)性,說(shuō)明不能從遺傳多樣性的角度很好地解釋不同種質(zhì)間對(duì)脅迫響應(yīng)差異,其原因還有待于進(jìn)一步研究。

    參考文獻(xiàn)(References)

    [1]Veilleux R E,Shen L Y,Paz M M.Analysis of the genetic composition of anther-derived potato by randomly amplified polymorphic DNA and simple sequence repeat.Genome,1995,38(6):1153-1162.

    [2]Provan J,Powell W,Waugh R.Microsatellite analysis of relationships within cultivated potato (Solanumtuberosum).Theoretical and Applied Genetics,1996,92(8):1078-1084.

    [3]Rocha E A,Paiva L V,de Carvalho H H,Guimar?es C T.Molecular characterization and genetic diversity of potato cultivars using SSR and RAPD markers.Crop Breeding and Applied Biotechnology,2010,10(3):204-210.

    [4]de Galarreta J I R,Barandalla L,Rios D J,Lopez R,Ritter E.Genetic relationships among local potato cultivars from Spain using SSR markers.Genetic Resources and Crop Evolution,2011,58(3):383-395.

    [5]Ghislain M,Spooner D M,Rodríguez F,Villamn F,Nunez J,Vásquez C,Waugh R,Bonierbale M.Selection of highly informative and user-friendly microsatellites (SSRs) for genotyping of cultivated potato.Theoretical and Applied Genetics,2004,108(5):881-890.

    [6]Spooner D M,Núez J,Trujillo G,María del R H,Guzmán F,Ghislain M.Extensive simple sequence repeat genotyping of potato landraces supports a major reevaluation of their gene pool structure and classification.Proceedings of the National Academy of Sciences of USA,2007,104(49):19398-19403.

    [7]劉建霞,雷海英,溫日宇,郭耀東,王潤(rùn)梅,周鳳.山西省馬鈴薯主栽品種遺傳多樣性的SSR分析.華北農(nóng)學(xué)報(bào),2013,27(6):72-77.

    Liu J X,Lei H Y,Wen R Y,Guo Y D,Wang R M,Zhou F.SSR analysis on genetic diversity of leading potato cultivars in Shanxi Province.Acta Agriculturae Boreali-sinica,2013,27(6):72-77.(in Chinese)

    [8]葉景秀,張鳳軍,張永成.青海省20個(gè)主要馬鈴薯審定品種的SSR標(biāo)記遺傳分析.種子,2013,32(6):1-4.

    Ye J X,Zhang F J,Zhang Y C.Genetic diversity analysis of 20 potato (SolanumtuberosumL.) genetypes in Qinghai Province with SSR markers.Seeds,2013,32(6):1-4.(in Chinese)

    [9]段艷鳳,劉杰,卞春松,段紹光,徐建飛,金黎平.中國(guó)88個(gè)馬鈴薯審定品種SSR指紋圖譜構(gòu)建與遺傳多樣性分析.作物學(xué)報(bào),2009,35(8):1451-1457.

    Duan Y F,Liu J,Bian C S,Duan S G,Xu J F,Jin L P.Construction of fingerprinting and analysis of genetic diversity with SSR markers for eighty-eight approved potato cultivars (SolanumtuberosumL.) in China.Acta Agronomica Sinica,2009,35(8):1451-1457.(in Chinese)

    [10]陳立強(qiáng),師尚禮.42份紫花苜蓿種質(zhì)資源遺傳多樣性的SSR分析.草業(yè)科學(xué),2015,32(3):372-381.

    Chen L Q,Shi S L.Genetic diversity among 42 alfalfa accession revealed by SSR marker.Pratacultural Science,2015,32(3):372-381.(in Chinese)

    [11]黃素華.玉米自交系多樣性研究以及玉米苗期耐旱相關(guān)基因的克隆、分析.北京:中國(guó)農(nóng)業(yè)科學(xué)院博士后學(xué)位論文,2005.

    Huang S H.Study on diversity of maize (ZeamaysL.) inbred lines and cloning of related genes for drought resistance(tolerance)from corn seedling.Post Doctoral Thesis.Beijing:Chinese Academy of Agricultural Sciences,2005.(in Chinese)

    [12]王陽(yáng),劉成,王天宇,石云素,宋燕春,黎裕.干旱脅迫和正常灌溉條件下玉米產(chǎn)量性狀的QTL分析.植物遺傳資源學(xué)報(bào),2007,8(2):179-183.

    Wang Y,Liu C,Wang T Y,Shi Y S,Song Y C,Li Y.QTL analysis of yield components in maize under different water regimes.Journal of Plant Genetic Resources,2007,8(2):179-183.(in Chinese)

    [13]張潔.冬小麥幼苗期水分脅迫應(yīng)答蛋白及其分子標(biāo)記的研究.楊凌:西北農(nóng)林科技大學(xué)碩士學(xué)位論文,2007.

    Zhang J.Molecular marker and response protein of winter wheat during seedling stage in the water stress condition.Master Thesis.Yangling:Northwest Agriculture and Forestry University,2007.(in Chinese)

    [14]胡國(guó)霞,馬蓮菊,陳強(qiáng),趙貴林,褚妍,任菲,李雪梅.植物抗氧化系統(tǒng)對(duì)水分脅迫及復(fù)水響應(yīng)研究進(jìn)展.安徽農(nóng)業(yè)科學(xué),2011,39(3):1278-1280.

    Hu G X,Ma L J,Chen Q,Zhao G L,Chu Y,Ren F,Li X M.Research progress of plant antioxidant system to water stress and rewatering. Journal of Anhui Agricultural Sciences,2011,39(3):1278-1280.(in Chinese)

    [15]趙海超,抗艷紅,龔學(xué)臣,尹江,田再民.干旱脅迫對(duì)不同馬鈴薯品種苗期生理生化指標(biāo)的影響.作物雜志,2013(6):63-69.

    Zhao H C,Kang Y H,Gong X C,Yin J,Tian Z M.Effects of drought stress on physiological and biochemical indices of different potato variety.Crops,2013(6):63-69.(in Chinese)

    [16]何淼,李文鶴,卓麗環(huán).野菊幼苗對(duì)自然干旱脅迫的生理響應(yīng).草業(yè)科學(xué),2011,28(8):1456-1460.

    He M,Li W H,Zhuo L H.Physiological response ofDendrathemaindicumseedlings to nature drought stress.Pratacultural Science,2011,28(8):1456-1460.(in Chinese)

    [17]Heyser J W,Nabors M W.Osmotic adjustment of tobacco cells and plants to penetrating and non-penetrating solutes.Plant Physiology,1979,63:77-77.

    [18]Carpita N,Sabularse D,Montezinos D,Delmer D P.Determination of the pore size of cell walls of living plant cells.Science,1979,205:1144-1147.

    [19]王瑾,劉桂茹,楊學(xué)舉.PEG脅迫下不同抗旱性小麥品種幼苗形態(tài)及主要理化特性的比較.河北農(nóng)業(yè)大學(xué)學(xué)報(bào),2005,28(5):6-10.

    Wang J,Liu G R,Yang X J.Comparison of the morphology and biochemical characteristics of wheat seedlings between drought-resistant and drought-sensitive varieties under PEG water-stress.Journal of Agricultural University of Hebei,2005,28(5):6-10.(in Chinese)

    [20]李建武,王蒂.水分脅迫對(duì)馬鈴薯試管苗抗氧化酶活性的影響.北方園藝,2008(1):7-9.

    Li J W,Wang D.The effect of water stress on the activity of antioxidant enzymes ofSolanumtuberosumL.Northern Horticulture,2008 (1):7-9.(in Chinese)

    [21]秦玉芝,陳玨,劉明月,熊興耀,何長(zhǎng)征,周玄,邢錚.聚乙二醇模擬干旱對(duì)馬鈴薯幼苗生長(zhǎng)與細(xì)胞膜透性的影響.湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,37(6):627-631.

    Qin Y Z,Chen Y,Liu M Y,Xiong X Y,He C Z,Zhou X,Xing Z.Effect of PEG-induced water stress on morphologic characters and membrane permeability of potato.Journal of Hunan Agricultural University(Natural Sciences),2011,37(6):627-631.(in Chinese)

    [22]Jones A S,Walker R T.Isolation and analysis of the deoxyribonucleic acid ofMycoplasmamycoidesvar.capri. Nature,1963,198:588-589.

    [23]Feingold S,Lloyd J,Norero N,Bonierbale M,Lorenzen J.Mapping and characterization of new EST-derived microsatellites for potato (SolanumtuberosumL.).Theoretical and Applied Genetics,2005,111(3):456-466.

    [24]鄒奇.植物生理實(shí)驗(yàn)指導(dǎo).北京:中國(guó)農(nóng)業(yè)出版社,2000,163-165.

    Zou Q.Experiment Instruction of Plant Physiology.Beijing:China Agricultural Press,2000:163-165.

    [25]Rolf F.Ntsys-Pc:Numerical Taxonomy and Multivariate Analysis System Version 2.1.Setauket,NY,USA:Exeter Publishing,2000.

    [26]Yeh F C,Yang R C,Boyle T.PopGene Version 131:Microsoft Window-based freeware for population genetic analysis.University of Alberta and Centre for International Forestry Research,1999:11-23.

    [27]Kimura M,Crow J F.The number of alleles that can be maintained in a finite population.Genetics,1964,49(4):725.

    [28]Lewontin R C.The Apportionment of Human Diversity.Evolutionary Biology.US:Springer,1995:381-398.

    [29]Nei M.Analysis of gene diversity in subdivided populations.Proceedings of the National Academy of Sciences of USA,1973,70(12):3321-3323.

    [30]嚴(yán)威凱.雙標(biāo)圖分析在農(nóng)作物品種多點(diǎn)試驗(yàn)中的應(yīng).作物學(xué)報(bào),2010,36(11):1805-1819.

    Yan W K.Optimal use of biplots in analysis of multi-location variety test data.Acta gronomica Sinica,2010,36(11):1805-1819.(in Chinese)

    [31]白江平,崔同霞,張金文,張俊蓮,王蒂.馬鈴薯塊SGAs合成途徑關(guān)鍵基因表達(dá)量的GGE雙標(biāo)圖分析.植物研究,2014,34(2):238-244.

    Bai J P,Cui T X,Zhang J W,Zhang J L,Wang D.GGE-Biplot analysis on gene expression of seven key genes mediating biosynthesis of steroidal glycoalkaloid in potato.Bulletin of Botanical Research,2014,34(2):238-244.(in Chinese)

    [32]Gopal J,Iwama K.In vitro screening of potato against water-stress mediated through sorbitol and polyethylene glycol.Plant Cell Reports,2007,26(5):693-700.

    [33]焦志麗,李勇,呂典秋,王晶英.不同程度干旱脅迫對(duì)馬鈴薯幼苗生長(zhǎng)和生理特性的影響.中國(guó)馬鈴薯,2012,25(6):329-333.

    Jiao Z L,Li Y,Lyu D Q,Wang J Y.Effects of different drought treatments on growth indicators and physiological characters of potato seedlings.Chinese Potato Journal,2012,25(6):329-333.(in Chinese)

    [34]抗艷紅,龔學(xué)臣,田再民,趙海超,齊海明.聚乙二醇處理馬鈴薯脫毒試管苗的生理反應(yīng).江蘇農(nóng)業(yè)科學(xué),2011(2):162-164.

    Kang Y H,Gong X C,Tian Z M,Zhao H C,Qi H M.Physiological response of potato virus-free test-tube plantlets under peg-induced water stress.Jiangsu Agricultural Sciences,2011(2):162-164.(in Chinese)

    [35]李建武,王蒂,司懷軍,彭曉莉.水分脅迫下馬鈴薯試管苗的生理響應(yīng).甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào),2005(3):319-323.

    Li J W,Wang D,Si H J,Peng X L.Physiological response of test-tuber plantlets of potato under water stress.Journal of Gansu Agricultural University,2005(3):319-323.(in Chinese)

    [36]Devi R,Kaur N,Gupta A K.Potential of antioxidant enzymes in depicting drought tolerance of wheat (TriticumaestivumL.).Indian Journal of Biochemistry & Biophysics,2012,49(4):257-265.

    [37]Mirzaee M,Moieni A,Ghanati F.Effects of Drought stress on the lipid peroxidation and antioxidant enzyme activities in two canola(BrassicanapusL.)Cultivars.Journal of Agricultural Science and Technology,2013,15(3):593-602.

    [38]張立軍,樊金娟,阮燕曄,關(guān)義新.聚乙二醇在植物滲透脅迫生理研究中的應(yīng)用.植物生理學(xué)通訊,2004,40(3):361-364.

    Zhang L J,Fan J J,Ruan Y Y,Guan Y X.Application of polyethylene glycol in the study of plant osmotic stress physiology.Plant Physiology Communications,2004,40(3):361-364.(in Chinese)

    [39]楊文柱,焦燕.SSR分子標(biāo)記技術(shù)在生物遺傳學(xué)領(lǐng)域的應(yīng)用.安徽農(nóng)業(yè)科學(xué),2012,40(2):640-642.

    Yang W Z,Jiao Y.Application of SSR molecular marker technique in the field of biogenetics.Journal of Anhui Agricultural Sciences,2012,40(2):640-642.(in Chinese)

    [40]張斌,秦嶺,王瑩瑩,范仲學(xué),畢玉平.作物抗旱分子機(jī)制研究進(jìn)展.山東農(nóng)業(yè)科學(xué),2012(12):22-26.

    Zhang B,Qin L,Wang Y Y,Fan Z X,Bi Y P.Research progress of crop drought-resistant mechanism.Shandong Agricultural Sciences,2012(12):22-26.(in Chinese)

    (責(zé)任編輯張瑾)

    DOI:10.11829/j.issn.1001-0629.2015-0182

    *收稿日期:2015-03-31接受日期:2015-08-20

    基金項(xiàng)目:科技部國(guó)際科技合作項(xiàng)目(2014DFG31570);甘肅科技基金(1308RJZA131、1308RJIA005);蘭州科技研究項(xiàng)目(2013-4-156)

    通信作者:王蒂(1955-),男,陜西延安人,教授,博導(dǎo),博士,研究方向?yàn)樽魑镞z傳育種研究。E-mail:wangd@gsau.edu.cn

    中圖分類號(hào):S532.03

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1001-0629(2016)3-0431-11* 1

    Corresponding author:Wang DiE-mail:wangd@gsau.edu.cn

    Analyze the genetic variation and drought resistance of the potato germplasm

    Lou Yan1,2, Bai Jiang-ping1,2, Yang Hong-yu1,3, Li Wei-ting1,2, Gao Hui-juan1,2,Zhang Jun-lian1,3, Wang Di1,2, Zhang Jin-lin4, Wang Wang-tian2

    (1.Gansu Key Laboratory of Crop Improvement and Germplasm Enhancement/Gansu Provincial Key Laboratory of Aridland Crop Science, Lanzhou 730070, China;2.College of Agronomy, Gansu Agricultural University, Lanzhou 730070, China;3.College of Horticulture, Gansu Agricultural University, Lanzhou 730070, China;4.College of pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730000, China)

    Abstract:To explore the relationship between potato (Solanum tuberosum) different genetic traits and its response to water stress, we analyzed the genetic diversity of test materials by the methods of SSR, studied the content variation of peroxidase (POD and CAT) activity and soluble protein of different potato varieties under water stresses by processing potato transplanting seedling for PEG-6000 simulated water stress and natural water stress method. The results showed that the genetic similarity coefficient (GS) at 0.700 level, 44 potato materials can be divided into five types, most breed clustered together, explaining the genetic basis of testing varieties is still relatively narrow. Through a GGE-biplot comprehensive analysis of soluble protein, POD and CAT variation, 44 selected materials were divided into 6 groups. Comparative analysis of genetic diversity analysis and physiological results found that there were no obvious consistency between them, namely the species differences in genetic traits can’t react their drought resistance performance difference directly. But two stress modes with high consistency on studying the characteristics of varieties of drought resistance. This research provides the reference for further drought resistance research.

    Key words:genetic diversity; PEG stress; GGE-biplot; drought resistance

    婁艷,白江平,楊宏羽,李維婷,高慧娟,張俊蓮,王蒂,張金林,王旺田.馬鈴薯種質(zhì)的遺傳特性與抗旱性的關(guān)系.草業(yè)科學(xué),2016,33(3):431-441.

    Lou Y,Bai J P,Yang H Y,Li W T,Gao H J,Zhang J L,Wang D,Zhang J L,Wang W T.Analyze the genetic variation and drought resistance of the potato germplasm.Pratacultural Science,2016,33(3):431-441.

    第一作者:婁艷(1989-),女,甘肅臨夏人,在讀碩士生,研究方向?yàn)轳R鈴薯遺傳育種。E-mail:lou_yan0931@163.com

    共同第一作者:白江平(1978-),男,甘肅天水人,副教授,碩導(dǎo),博士,研究方向?yàn)樽魑镞z傳育種。E-mail:baijp@gsau.edu.cn

    猜你喜歡
    遺傳多樣性抗旱性
    山東農(nóng)業(yè)大學(xué)揭示了PLATZ4轉(zhuǎn)錄因子調(diào)控植物抗旱性的新機(jī)制
    蔬菜(2023年5期)2023-07-29 02:52:35
    不同光質(zhì)對(duì)黃瓜幼苗抗旱性的影響
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    基于稱重法的煙草伸根期抗旱性鑒定
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    樹木抗旱性及抗旱造林技術(shù)研究綜述
    山西大豆自然群體遺傳多樣性的研究
    国产色婷婷99| 成人18禁在线播放| 黑人欧美特级aaaaaa片| 激情在线观看视频在线高清| 国产一区二区在线观看日韩 | 在线免费观看不下载黄p国产 | 亚洲不卡免费看| 久久精品亚洲精品国产色婷小说| 91久久精品电影网| 成年女人永久免费观看视频| 热99re8久久精品国产| 亚洲在线观看片| 国产乱人伦免费视频| 欧美日本亚洲视频在线播放| 观看美女的网站| 免费在线观看亚洲国产| 9191精品国产免费久久| 大型黄色视频在线免费观看| 有码 亚洲区| 免费观看的影片在线观看| svipshipincom国产片| 国内精品一区二区在线观看| 2021天堂中文幕一二区在线观| 国产一区二区激情短视频| 成年人黄色毛片网站| 一本精品99久久精品77| 精品一区二区三区视频在线观看免费| 十八禁网站免费在线| 国产亚洲av嫩草精品影院| 级片在线观看| 欧美三级亚洲精品| 国产成年人精品一区二区| 日韩高清综合在线| 国产又黄又爽又无遮挡在线| 日韩中文字幕欧美一区二区| www.999成人在线观看| 国产亚洲av嫩草精品影院| 久久人妻av系列| 欧美区成人在线视频| 黄色女人牲交| 国产亚洲精品久久久com| 亚洲男人的天堂狠狠| 天天一区二区日本电影三级| 国产一区二区在线观看日韩 | 9191精品国产免费久久| 欧美又色又爽又黄视频| 日本免费a在线| e午夜精品久久久久久久| 午夜免费男女啪啪视频观看 | 精品午夜福利视频在线观看一区| h日本视频在线播放| 男女下面进入的视频免费午夜| 少妇裸体淫交视频免费看高清| 国产日本99.免费观看| 亚洲国产欧美网| 亚洲专区中文字幕在线| 欧美一区二区亚洲| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看 | 免费大片18禁| 校园春色视频在线观看| 熟女电影av网| 女生性感内裤真人,穿戴方法视频| 在线观看免费午夜福利视频| 午夜免费男女啪啪视频观看 | 久久伊人香网站| 网址你懂的国产日韩在线| av专区在线播放| 人妻夜夜爽99麻豆av| www.999成人在线观看| 99在线视频只有这里精品首页| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 99热只有精品国产| 最后的刺客免费高清国语| 国产av一区在线观看免费| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 欧美乱妇无乱码| 国内精品久久久久精免费| 亚洲成人久久性| 哪里可以看免费的av片| 欧美黄色淫秽网站| 变态另类成人亚洲欧美熟女| 免费电影在线观看免费观看| 午夜激情欧美在线| 老司机深夜福利视频在线观看| 97人妻精品一区二区三区麻豆| 一夜夜www| 老司机福利观看| 精品福利观看| 19禁男女啪啪无遮挡网站| 在线观看日韩欧美| 成人永久免费在线观看视频| 成人三级黄色视频| 亚洲精品影视一区二区三区av| 国产av一区在线观看免费| 一区二区三区激情视频| 舔av片在线| 亚洲一区高清亚洲精品| 又紧又爽又黄一区二区| 亚洲无线在线观看| 特大巨黑吊av在线直播| 精品国产亚洲在线| 国产精品自产拍在线观看55亚洲| 日本在线视频免费播放| 成人av一区二区三区在线看| 香蕉丝袜av| 99久久综合精品五月天人人| 99视频精品全部免费 在线| 深夜精品福利| 亚洲人成电影免费在线| 欧美绝顶高潮抽搐喷水| 亚洲av美国av| xxxwww97欧美| 波野结衣二区三区在线 | 亚洲av不卡在线观看| h日本视频在线播放| 在线免费观看的www视频| 最近在线观看免费完整版| 中文字幕高清在线视频| 97人妻精品一区二区三区麻豆| 免费搜索国产男女视频| 99热精品在线国产| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| av专区在线播放| 国产三级在线视频| 日韩av在线大香蕉| 日日摸夜夜添夜夜添小说| 日本五十路高清| 最近最新中文字幕大全免费视频| 色尼玛亚洲综合影院| 精品国产亚洲在线| 免费观看人在逋| 欧美另类亚洲清纯唯美| 精品一区二区三区人妻视频| 亚洲av一区综合| 日本三级黄在线观看| 午夜福利在线观看免费完整高清在 | 麻豆国产97在线/欧美| 99热6这里只有精品| 国产精品乱码一区二三区的特点| 国产精品,欧美在线| 亚洲人成电影免费在线| 免费无遮挡裸体视频| 18禁国产床啪视频网站| 亚洲成a人片在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 岛国视频午夜一区免费看| 亚洲人成电影免费在线| 国产精品99久久99久久久不卡| 精品乱码久久久久久99久播| 国产亚洲欧美98| 婷婷亚洲欧美| 国产精品三级大全| 精品久久久久久成人av| 精品电影一区二区在线| 嫩草影院入口| 99久久精品一区二区三区| 成人三级黄色视频| 亚洲国产色片| 99国产精品一区二区三区| 免费av毛片视频| 1024手机看黄色片| 国产精品亚洲美女久久久| 黄色女人牲交| 国产视频内射| av福利片在线观看| 噜噜噜噜噜久久久久久91| 特大巨黑吊av在线直播| 老熟妇仑乱视频hdxx| 国产伦人伦偷精品视频| 国产探花在线观看一区二区| 男人舔女人下体高潮全视频| 亚洲成人久久爱视频| 欧美乱色亚洲激情| www.999成人在线观看| 美女高潮的动态| 又紧又爽又黄一区二区| 天天添夜夜摸| 国产成人av激情在线播放| 人人妻人人澡欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 免费观看精品视频网站| 国产精品香港三级国产av潘金莲| 亚洲人与动物交配视频| 小说图片视频综合网站| 少妇裸体淫交视频免费看高清| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 搡老岳熟女国产| 性色avwww在线观看| a级一级毛片免费在线观看| 亚洲精品成人久久久久久| 亚洲最大成人手机在线| 色老头精品视频在线观看| 深夜精品福利| 国产久久久一区二区三区| 欧美一区二区精品小视频在线| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 久久精品国产亚洲av香蕉五月| 熟女人妻精品中文字幕| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 精品人妻一区二区三区麻豆 | 亚洲av免费高清在线观看| 一本综合久久免费| 国产三级黄色录像| 真实男女啪啪啪动态图| 久久精品影院6| www国产在线视频色| 在线天堂最新版资源| 深爱激情五月婷婷| 高潮久久久久久久久久久不卡| 亚洲自拍偷在线| 日本一二三区视频观看| 国产精品乱码一区二三区的特点| 午夜福利免费观看在线| 99久久成人亚洲精品观看| 亚洲专区中文字幕在线| 给我免费播放毛片高清在线观看| 老熟妇乱子伦视频在线观看| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 91字幕亚洲| 亚洲精品影视一区二区三区av| 在线国产一区二区在线| 97超视频在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 国产精品久久久久久人妻精品电影| 美女被艹到高潮喷水动态| 国产不卡一卡二| 99精品久久久久人妻精品| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 国产欧美日韩精品一区二区| 久久中文看片网| 成人鲁丝片一二三区免费| 一区二区三区激情视频| 69av精品久久久久久| 国产激情欧美一区二区| 久久婷婷人人爽人人干人人爱| 午夜激情欧美在线| 日韩欧美 国产精品| 亚洲五月天丁香| 国产亚洲精品一区二区www| 日韩欧美三级三区| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区| 无人区码免费观看不卡| 身体一侧抽搐| 99riav亚洲国产免费| 国产熟女xx| 国产欧美日韩精品亚洲av| 国产精品久久电影中文字幕| 一级a爱片免费观看的视频| 国产精品美女特级片免费视频播放器| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| 九九热线精品视视频播放| 男女那种视频在线观看| 他把我摸到了高潮在线观看| 内地一区二区视频在线| 欧美在线黄色| 性欧美人与动物交配| ponron亚洲| 国模一区二区三区四区视频| 搞女人的毛片| 国产精品自产拍在线观看55亚洲| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 亚洲人成网站在线播| 欧美又色又爽又黄视频| 国产麻豆成人av免费视频| 国产精品亚洲一级av第二区| 69av精品久久久久久| 在线播放无遮挡| 亚洲国产欧洲综合997久久,| 黑人欧美特级aaaaaa片| 黄片小视频在线播放| 国产精品久久久久久久久免 | 中文字幕人成人乱码亚洲影| 国产69精品久久久久777片| 老熟妇仑乱视频hdxx| 午夜免费激情av| 欧美精品啪啪一区二区三区| 制服丝袜大香蕉在线| 国产激情欧美一区二区| 少妇的逼好多水| 日韩欧美国产在线观看| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 在线天堂最新版资源| 99久久无色码亚洲精品果冻| 国产伦一二天堂av在线观看| 国产一区二区三区在线臀色熟女| 免费一级毛片在线播放高清视频| 老汉色av国产亚洲站长工具| 怎么达到女性高潮| 免费高清视频大片| 国产色婷婷99| 精华霜和精华液先用哪个| 中国美女看黄片| 伊人久久精品亚洲午夜| 一夜夜www| 欧美日本视频| 亚洲无线在线观看| 色尼玛亚洲综合影院| 国产高清视频在线播放一区| 久久久色成人| 757午夜福利合集在线观看| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久久久免 | 国产精品av视频在线免费观看| 女人被狂操c到高潮| 国产精品久久久久久亚洲av鲁大| 99精品久久久久人妻精品| 久久精品91蜜桃| 国产精品 国内视频| 搞女人的毛片| 51午夜福利影视在线观看| 人人妻人人看人人澡| 亚洲熟妇熟女久久| 成人欧美大片| 床上黄色一级片| 99热这里只有是精品50| 精品乱码久久久久久99久播| 黑人欧美特级aaaaaa片| 国产高清视频在线观看网站| 亚洲在线自拍视频| 国产精品一区二区三区四区久久| 两个人看的免费小视频| 岛国视频午夜一区免费看| 精品人妻一区二区三区麻豆 | 99久久精品国产亚洲精品| av在线蜜桃| 午夜老司机福利剧场| 老熟妇乱子伦视频在线观看| 嫩草影院入口| 在线a可以看的网站| 欧美黑人巨大hd| 亚洲国产精品合色在线| 亚洲黑人精品在线| 免费av不卡在线播放| 夜夜爽天天搞| 黄色成人免费大全| 日本 欧美在线| 老熟妇乱子伦视频在线观看| 91久久精品国产一区二区成人 | 久久久国产精品麻豆| 国产精品电影一区二区三区| 国产一区二区亚洲精品在线观看| 99久久精品一区二区三区| 99在线人妻在线中文字幕| 欧美一级a爱片免费观看看| 国产成年人精品一区二区| 男女那种视频在线观看| 亚洲av成人av| 99久久精品一区二区三区| 国产黄a三级三级三级人| 久久久精品大字幕| 亚洲精品美女久久久久99蜜臀| 亚洲第一电影网av| 精品国产美女av久久久久小说| 日本三级黄在线观看| 搡女人真爽免费视频火全软件 | 亚洲精品一区av在线观看| 18禁黄网站禁片午夜丰满| 人妻夜夜爽99麻豆av| 日韩亚洲欧美综合| 91久久精品电影网| 中文字幕精品亚洲无线码一区| 亚洲精品亚洲一区二区| 亚洲成a人片在线一区二区| 国产av在哪里看| 很黄的视频免费| av天堂在线播放| 国产中年淑女户外野战色| 亚洲av成人不卡在线观看播放网| 国产日本99.免费观看| av天堂在线播放| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 国产激情欧美一区二区| 露出奶头的视频| 一级毛片高清免费大全| 欧美性感艳星| 又黄又爽又免费观看的视频| 精品久久久久久久末码| 国产精品综合久久久久久久免费| 在线观看美女被高潮喷水网站 | 亚洲精品一卡2卡三卡4卡5卡| 免费人成视频x8x8入口观看| av天堂中文字幕网| 俺也久久电影网| 国产精品久久久久久久电影 | 神马国产精品三级电影在线观看| 美女 人体艺术 gogo| 精品人妻1区二区| 最新在线观看一区二区三区| eeuss影院久久| 成人亚洲精品av一区二区| 国内揄拍国产精品人妻在线| 最新美女视频免费是黄的| 亚洲精品色激情综合| 观看免费一级毛片| 亚洲精品粉嫩美女一区| 免费av观看视频| 亚洲av一区综合| 欧美精品啪啪一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产视频一区二区在线看| 久久九九热精品免费| 日韩欧美 国产精品| 美女 人体艺术 gogo| 日本与韩国留学比较| 女人被狂操c到高潮| 3wmmmm亚洲av在线观看| 国产在视频线在精品| 国产91精品成人一区二区三区| 亚洲人成网站高清观看| 久久久久久国产a免费观看| netflix在线观看网站| 在线观看美女被高潮喷水网站 | 精品国产超薄肉色丝袜足j| 国产淫片久久久久久久久 | 啦啦啦免费观看视频1| 精品久久久久久成人av| 久久久久久久久久黄片| 窝窝影院91人妻| 午夜福利免费观看在线| 国产精品乱码一区二三区的特点| 欧美xxxx黑人xx丫x性爽| 成年女人看的毛片在线观看| bbb黄色大片| 精品人妻1区二区| 国产v大片淫在线免费观看| 男人舔女人下体高潮全视频| xxx96com| a级毛片a级免费在线| 欧美日韩国产亚洲二区| 国产精品久久久久久人妻精品电影| 精品不卡国产一区二区三区| 麻豆成人午夜福利视频| 午夜福利18| 亚洲不卡免费看| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美网| 亚洲欧美一区二区三区黑人| 亚洲人成网站在线播放欧美日韩| 欧美日韩一级在线毛片| 99热这里只有精品一区| 九九久久精品国产亚洲av麻豆| 真实男女啪啪啪动态图| 色综合站精品国产| 国产精品乱码一区二三区的特点| a在线观看视频网站| 久久久久久人人人人人| www国产在线视频色| 黄色女人牲交| 国产av麻豆久久久久久久| or卡值多少钱| 国产黄片美女视频| 毛片女人毛片| 婷婷精品国产亚洲av在线| 久久草成人影院| 久久婷婷人人爽人人干人人爱| 亚洲在线观看片| 精品久久久久久久毛片微露脸| 成人高潮视频无遮挡免费网站| 日韩高清综合在线| 男人舔奶头视频| 精品日产1卡2卡| 色吧在线观看| 99久久99久久久精品蜜桃| 中文字幕人妻熟人妻熟丝袜美 | 国产午夜精品论理片| 夜夜爽天天搞| 12—13女人毛片做爰片一| 色视频www国产| 久久这里只有精品中国| 国产精品电影一区二区三区| 精品免费久久久久久久清纯| 日本黄色视频三级网站网址| 51国产日韩欧美| eeuss影院久久| 久久久久久久亚洲中文字幕 | 欧美色欧美亚洲另类二区| 男女下面进入的视频免费午夜| 精品久久久久久久久久免费视频| 亚洲av电影在线进入| 亚洲av五月六月丁香网| av女优亚洲男人天堂| 色在线成人网| 九色国产91popny在线| 欧美性猛交黑人性爽| 国产精品乱码一区二三区的特点| 毛片女人毛片| 色在线成人网| 国产高清三级在线| 听说在线观看完整版免费高清| 国产精品98久久久久久宅男小说| av女优亚洲男人天堂| 在线视频色国产色| 久久久久久久亚洲中文字幕 | 国产精品亚洲av一区麻豆| 国产伦精品一区二区三区视频9 | 久久久久久久久久黄片| 久久久久免费精品人妻一区二区| 最后的刺客免费高清国语| 嫩草影视91久久| 欧美中文综合在线视频| 国产国拍精品亚洲av在线观看 | 欧美日韩中文字幕国产精品一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 99riav亚洲国产免费| 波多野结衣高清无吗| 国产久久久一区二区三区| 国产精品三级大全| 国产精品一区二区三区四区免费观看 | 搞女人的毛片| 欧美乱妇无乱码| 亚洲成人精品中文字幕电影| 两个人看的免费小视频| 亚洲av一区综合| 日本黄色视频三级网站网址| 欧美丝袜亚洲另类 | 无遮挡黄片免费观看| 最新美女视频免费是黄的| 亚洲性夜色夜夜综合| 欧美区成人在线视频| 国内精品久久久久久久电影| www.色视频.com| 看黄色毛片网站| 精品久久久久久久末码| 99国产极品粉嫩在线观看| 国产精品 欧美亚洲| 久久久国产精品麻豆| 黄色成人免费大全| 免费在线观看成人毛片| 18禁国产床啪视频网站| 国产伦一二天堂av在线观看| 免费观看的影片在线观看| 日韩成人在线观看一区二区三区| 国产成人av激情在线播放| 亚洲欧美一区二区三区黑人| 两个人的视频大全免费| 精品99又大又爽又粗少妇毛片 | 美女黄网站色视频| 亚洲无线观看免费| 成年人黄色毛片网站| 精华霜和精华液先用哪个| 看免费av毛片| 亚洲欧美日韩卡通动漫| 日韩人妻高清精品专区| 黄色丝袜av网址大全| 欧美性猛交黑人性爽| 又紧又爽又黄一区二区| 日本黄色片子视频| 亚洲国产欧洲综合997久久,| 岛国在线观看网站| 欧美日韩一级在线毛片| 黄色女人牲交| 久久久久免费精品人妻一区二区| 亚洲最大成人手机在线| 国产不卡一卡二| 成人精品一区二区免费| 日韩高清综合在线| 欧美av亚洲av综合av国产av| 欧美最黄视频在线播放免费| 69av精品久久久久久| 欧美+亚洲+日韩+国产| 色尼玛亚洲综合影院| 露出奶头的视频| 亚洲中文字幕日韩| 老司机在亚洲福利影院| 亚洲内射少妇av| 国产单亲对白刺激| 久久午夜亚洲精品久久| 亚洲午夜理论影院| 久久精品国产自在天天线| 欧美在线一区亚洲| 亚洲男人的天堂狠狠| 亚洲一区二区三区色噜噜| av在线天堂中文字幕| 午夜久久久久精精品| 每晚都被弄得嗷嗷叫到高潮| 欧美激情在线99| 日本三级黄在线观看| av天堂在线播放| 青草久久国产| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久人人做人人爽| e午夜精品久久久久久久| 国产免费一级a男人的天堂| 午夜福利在线观看吧| 色综合婷婷激情| 色吧在线观看| 亚洲avbb在线观看| 欧美av亚洲av综合av国产av| 麻豆成人午夜福利视频| 久久精品91蜜桃| 夜夜爽天天搞| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 欧美另类亚洲清纯唯美| 最好的美女福利视频网| 18禁黄网站禁片午夜丰满| 成年女人看的毛片在线观看| 嫩草影院精品99| 在线播放国产精品三级| 精华霜和精华液先用哪个|