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

    大連地區(qū)島嶼與大陸玉竹種群遺傳多樣性的ISSR分析

    2017-11-10 01:47:17郝久程穆曉紅張恒慶
    植物研究 2017年5期
    關(guān)鍵詞:玉竹海島種群

    郝久程 賈 鑫 穆曉紅 張恒慶

    (遼寧師范大學(xué)生命科學(xué)學(xué)院,大連 116081)

    大連地區(qū)島嶼與大陸玉竹種群遺傳多樣性的ISSR分析

    郝久程 賈 鑫 穆曉紅 張恒慶*

    (遼寧師范大學(xué)生命科學(xué)學(xué)院,大連 116081)

    選取大連地區(qū)大陸與海島共有植物玉竹為研究對(duì)象,采用ISSR分子標(biāo)記技術(shù)對(duì)來自5個(gè)海島和4個(gè)大陸種群的262個(gè)玉竹個(gè)體進(jìn)行遺傳多樣性的比較和分析。從10個(gè)篩選出的ISSR引物擴(kuò)增得到120個(gè)位點(diǎn)信息,其中多態(tài)性條帶百分率為91.67%,Nei’s基因多樣性指數(shù)(h)為0.346 0,Shannon信息指數(shù)(I)為0.510 8。其遺傳分化系數(shù)(Gst)為0.117 4,基因流(Nm)為3.758 5。研究結(jié)果表明玉竹天然種群的遺傳多樣性較為豐富,種群間基因交流較為頻繁,遺傳距離與地理距離具有一定的相關(guān)性。通過海島與大陸種群遺傳多樣性的比對(duì)發(fā)現(xiàn),海島種群的遺傳多樣性略高于大陸種群,表明在孤立的生境和更為復(fù)雜的選擇壓力下,海島玉竹種群可能會(huì)積累更多的遺傳變異從而形成較高的遺傳多樣性水平。本文研究結(jié)果將為進(jìn)一步探討隔離生境中天然植物種群遺傳進(jìn)化規(guī)律提供證據(jù)。

    ISSR-PCR;玉竹;遺傳多樣性;遺傳結(jié)構(gòu)

    玉竹(Polygonatumodoratum)是百合科(Liliaceae)黃精屬(Polygonatum)多年生草本植物。野生玉竹生于涼爽、濕潤(rùn)、無積水、土壤肥沃的林下或山野陰坡地,宜溫暖濕潤(rùn)氣候,喜陰濕環(huán)境,較耐寒[1]。

    大連位于遼寧省遼東半島南端,氣候適宜,植物資源豐富。玉竹在我國(guó)分布較廣,屬于歐亞大陸溫帶地區(qū)廣布種,在大連地區(qū)的分布范圍也比較廣泛,在采樣過程中發(fā)現(xiàn)北至莊河仙人洞國(guó)家級(jí)自然保護(hù)區(qū)南至旅順蛇島均有玉竹種群的分布。國(guó)內(nèi)外關(guān)于玉竹方面的研究很多,主要集中于易混種鑒定[2]、化學(xué)成分[3]、藥理作用[4]、栽培技術(shù)及應(yīng)用開發(fā)等方面。

    簡(jiǎn)單重復(fù)序列間擴(kuò)增即ISSR(inter-simple sequence repeat)是自加拿大蒙特利爾大學(xué)的Zietkiewicz等[5]于1994年創(chuàng)建的,它是一種建立在PCR反應(yīng)基礎(chǔ)上的新型的分子標(biāo)記,具有操作簡(jiǎn)單、多態(tài)性豐富、試驗(yàn)重復(fù)性好的特點(diǎn)。國(guó)外關(guān)于玉竹遺傳多樣性方面的研究還未見報(bào)道。國(guó)內(nèi)有對(duì)玉竹遺傳多樣性方面的研究,但全部集中于大陸間的不同省份之間或不同栽培品種間的玉竹遺傳多樣性比較,路放[6]利用ISSR分子標(biāo)記技術(shù)對(duì)吉林地區(qū)8個(gè)品系的栽培玉竹品種進(jìn)行分析,發(fā)現(xiàn)不同品系玉竹種間存在較高的多態(tài)性,遺傳多樣性較為豐富;卜靜等[7]通過對(duì)我國(guó)5個(gè)省份共11個(gè)種群的野生玉竹和1個(gè)種群的栽培玉竹進(jìn)行了ISSR分析,發(fā)現(xiàn)不同產(chǎn)地野生玉竹種間遺傳多樣性較為豐富,對(duì)環(huán)境變化的適應(yīng)能力強(qiáng),野生玉竹的遺傳多樣性高于栽培居群的遺傳多樣性。對(duì)于玉竹在海島隔離生境的遺傳多樣性研究還沒有報(bào)道,也沒有相關(guān)學(xué)者對(duì)遼寧省地區(qū)內(nèi)野生玉竹遺傳多樣性進(jìn)行研究。雖然我國(guó)玉竹分布廣泛,但是隨著人們不斷的采挖野生玉竹根莖進(jìn)行栽培,野生玉竹的資源也在減少,遺傳多樣性會(huì)逐漸喪失。

    本實(shí)驗(yàn)采用ISSR分子標(biāo)記技術(shù)對(duì)遼寧省大連市的6個(gè)縣市(區(qū))共9個(gè)玉竹天然種群(其中5個(gè)為海島種群,4個(gè)為大陸種群)的遺傳多樣性進(jìn)行了比較和分析,獲取9個(gè)種群的遺傳學(xué)數(shù)據(jù),在物種水平和種群水平上揭示玉竹種群遺傳多樣性的高低,揭示玉竹種群的遺傳結(jié)構(gòu)并探討其形成的原因。同時(shí),比較研究海島和大陸玉竹種群的遺傳結(jié)構(gòu)、遺傳多樣性水平、遺傳分化的特點(diǎn),其結(jié)果將為探討隔離生境中天然植物種群遺傳進(jìn)化規(guī)律提供證據(jù)。

    1 實(shí)驗(yàn)材料

    實(shí)驗(yàn)材料于2016年5~6月分別采自于大連地區(qū)9個(gè)不同的地點(diǎn),分單株采集,樣品間距離不少于20 m,采樣時(shí)記錄周邊環(huán)境,折斷莖部整株采集,采集后帶回實(shí)驗(yàn)室進(jìn)行二次鑒定,確保采集樣品的準(zhǔn)確性,然后將材料的嫩葉采集后用封口帶單株封裝,保存于實(shí)驗(yàn)室-20℃冰柜中備用,詳細(xì)的采集信息見表1。

    表1玉竹樣品數(shù)量及采集地

    Table1ThesamplesizeandlocationofP.odoratumpopulation

    種群編號(hào)Code樣品數(shù)No.ofsamples緯度Latitude(N)經(jīng)度Longitude(E)地理位置LocationSD3038°57.083′120°58.937′旅順,蛇島(海島)Lvshun,Snakeisland(Island)XS2638°55.513′121°29.283′甘井子,西山(大陸)Ganjingzi,Westmount(Land)XRD3439°59.020′122°57.450′莊河,仙人洞(大陸)Zhuanghe,Xianrendong(Land)GLD2939°38.752′123°02.032′莊河,蛤蜊島(海島)Zhuanghe,Clamsisland(Island)WJD2839°26.753′123°05.243′莊河,王家島(海島)Zhuanghe,Wangjiaisland(Island)SCD3039°31.175′122°59.561′莊河,石城島(海島)Zhuanghe,Shichengisland(Island)XHS2539°18.276′121°54.449′金州,小黑山(大陸)Jinzhou,Xiaoheishan(Land)XPD2938°50.640′121°28.920′高新區(qū),小平島(大陸)High?techzone,Xiaopingisland(Island)DCS3139°15.051′122°34.925′長(zhǎng)海,大長(zhǎng)山島(海島)Changhai,Dachangshanisland(Island)

    2 實(shí)驗(yàn)方法

    2.1 DNA的提取

    玉竹DNA的提取采取張恒慶等[8]改良的CTAB法。粗提DNA樣品進(jìn)一步純化處理后,進(jìn)行瓊脂糖凝膠電泳,用美國(guó)冷泉港P41G型凝膠成像系統(tǒng)照相。采用Gel Pro1.1軟件對(duì)比樣品條帶與DL15000分子量標(biāo)準(zhǔn)的含量,確定DNA的含量和分子量,用雙蒸水稀釋成統(tǒng)一濃度后作為PCR反應(yīng)的模板。

    2.2 引物的篩選

    ISSR引物為上海生工生物工程公司合成,稀釋成10 μmol·L-1備用。在合成的27個(gè)引物中進(jìn)行了篩選,選取擴(kuò)增條帶穩(wěn)定、清晰、數(shù)量適中的10個(gè)引物作為本研究的引物,對(duì)全部樣品進(jìn)行擴(kuò)增。具體引物名稱及退火溫度見表2。

    表2引物名稱、擴(kuò)增位點(diǎn)數(shù)及最佳退火溫度

    Table2NumberofpolymorphiclociandoptimalannealingtemperatureofISSRprimer

    引物Primer序列Sequence擴(kuò)增位點(diǎn)數(shù)Locus退火溫度Annealingtemperature(℃)UBC817(CA)8A1152.0UBC818(CA)8G1451.3UBC841(GA)8YC1252.5UBC843(CT)8RA1150.3UBC844(CT)8RC1252.0UBC845(CT)8RG948.5UBC853(TC)8RT1152.5UBC855(AC)8YT1353.0UBC856(AC)8YA1253.0UBC857(AC)8YG1550.8

    2.3 ISSR反應(yīng)體系及PCR擴(kuò)增

    本實(shí)驗(yàn)參考陳玉秀等[9]的優(yōu)化結(jié)果,并通過各組成成分的梯度實(shí)驗(yàn),最終確定反應(yīng)體系如下:總體積25 μL,包括10 ng模板DNA,0.4 μmol·L-1ISSR引物,10×PCR buffer 2.5 μL,0.2 mmol·L-1dNTP混合物,2.0 mmol·L-1MgCl2,1.0 UTaqDNA聚合酶,最后加ddH2O補(bǔ)齊至25 μL。

    使用英國(guó)TECHNE TC-512型PCR循環(huán)儀進(jìn)行PCR擴(kuò)增,玉竹的ISSR-PCR反應(yīng)程序?yàn)?94℃ 5 min進(jìn)行預(yù)變性,35個(gè)循環(huán),條件為94℃變性45 s,48.5~53.0℃退火45 s(退火溫度因引物不同),72℃延伸90 s,最后在72℃保持7 min使擴(kuò)增的片段充分延伸。

    2.4 擴(kuò)增產(chǎn)物檢測(cè)及數(shù)據(jù)統(tǒng)計(jì)分析方法

    PCR反應(yīng)結(jié)束后,采用水平板瓊脂糖凝膠電泳來檢測(cè)擴(kuò)增結(jié)果。具體操作為:1×TBE電泳緩沖液,凝膠濃度為1.2%,以DL2000作為分子量標(biāo)準(zhǔn),80 V穩(wěn)壓電泳2.5 h。電泳結(jié)束后用凝膠成像系統(tǒng)照相,記錄電泳結(jié)果。

    根據(jù)電泳結(jié)果,選取多個(gè)位點(diǎn)記錄不同物種的同一位點(diǎn)上的條帶分布情況。記錄時(shí)只記錄那些電泳條帶清晰,并能在重復(fù)實(shí)驗(yàn)中穩(wěn)定出現(xiàn)的位點(diǎn),有帶的記為“1”,無帶的記為“0”,不確定的記為“.”,最后形成數(shù)據(jù)庫(kù),經(jīng)POPGENE1.32軟件[10]進(jìn)行遺傳多樣性分析,得出多態(tài)性條帶百分率(PPB)、Nei’s基因多樣性指數(shù)(h)、Shannon信息指數(shù)(I)、群體總基因多樣性(Ht)、群體內(nèi)基因多樣性(Hs)、遺傳分化系數(shù)(Gst)和基因流(Nm)等,同時(shí)利用軟件對(duì)玉竹種群進(jìn)行聚類分析。

    3 實(shí)驗(yàn)結(jié)果

    3.1 玉竹ISSR-PCR擴(kuò)增結(jié)果及遺傳多樣性分析

    本實(shí)驗(yàn)用篩選出的10個(gè)引物對(duì)大連地區(qū)的9個(gè)玉竹天然種群共262個(gè)個(gè)體進(jìn)行了PCR擴(kuò)增,選出的10個(gè)引物共檢測(cè)到120個(gè)位點(diǎn),每個(gè)引物檢測(cè)到的位點(diǎn)數(shù)介于9~15個(gè),擴(kuò)增的DNA片段長(zhǎng)度在150~2 000 bp,平均每個(gè)引物檢測(cè)到12個(gè)位點(diǎn)。在檢測(cè)到的120個(gè)位點(diǎn)中多態(tài)位點(diǎn)為110個(gè),其種內(nèi)多態(tài)位點(diǎn)比率為91.67%,9個(gè)種群內(nèi)多態(tài)位點(diǎn)比率在69.17%~87.50%,其中多態(tài)位點(diǎn)比率最高的是SCD(石城島)種群,最低的是XHS(小黑山)種群(表3),引物UBC845的部分?jǐn)U增結(jié)果見圖1。

    表3在種群和物種水平上玉竹遺傳多樣性ISSR分析結(jié)果

    Table3GeneticdiversityinP.odoratumdeterminedbyISSRmarkersatpopulationandspecieslevel

    種群名稱Population擴(kuò)增位點(diǎn)數(shù)Locus多態(tài)位點(diǎn)比率PercentageofpolymorphismlociNei’s指數(shù)Nei’sgenediversity(h)Shannon指數(shù)Shannon’sinformationindex(I)WJD9680.000.3163(0.1853)0.4625(0.2583)XHS8369.170.2524(0.1987)0.3746(0.2809)XS8974.170.2771(0.1972)0.4085(0.2765)SCD10587.500.3308(0.1673)0.4880(0.2270)SD9579.170.3011(0.1814)0.4448(0.2546)XPD10184.170.3215(0.1748)0.4737(0.2402)GLD9982.500.3103(0.1767)0.4585(0.2463)DCS9780.830.3120(0.1816)0.4589(0.2525)XRD10184.170.3206(0.1754)0.4718(0.2427)Specieslevel11091.670.3460(0.1532)0.5108(0.2039)

    圖1 部分玉竹種群ISSR電泳圖譜 M. DL2000Fig.1 ISSR electrophorogram of P.odoratum populations M. DL2000

    3.2 玉竹種群間的遺傳分化

    根據(jù)POPGENE1.32軟件計(jì)算,9個(gè)玉竹種群的遺傳分化計(jì)算結(jié)果為玉竹種群的Ht為0.345 2,Hs為0.304 7,遺傳分化系數(shù)Gst=0.117 4,說明有11.74%的遺傳變異發(fā)生在種群間,玉竹的遺傳變異主要發(fā)生在種群內(nèi),占總的遺傳變異的88.26%?;蛄?Nm)較大為3.758 5,說明種群間基因交流比較頻繁。

    3.3 玉竹的種群遺傳距離和聚類分析

    遺傳距離(D)和遺傳一致度(I)是衡量玉竹種群間的遺傳差異和遺傳相似性的重要指標(biāo),本研究運(yùn)用POPGENE1.32軟件進(jìn)行計(jì)算,得到9個(gè)玉竹種群間的遺傳距離(D)和遺傳一致度(I),如表4所示,我們可以看到,GLD與XRD的遺傳一致度最高,為0.978 2;XHS與SD的遺傳一致度最低,為0.887 0。

    根據(jù)所得到的Nei’s遺傳距離(D)的數(shù)據(jù),利用POPGENE1.32軟件,采用非加權(quán)組平均法(unweighted pair-group method with arithmetic means,UPGMA),對(duì)9個(gè)玉竹種群的進(jìn)行聚類分析,從而得到種群間的遺傳關(guān)系聚類圖(圖2)。從圖2中可以看出,9個(gè)玉竹種群被分為4個(gè)主要集群,其中SD和XHS種群為獨(dú)立的兩個(gè)集群,XS和XPD聚為一個(gè)集群,其他種群構(gòu)成一個(gè)集群。

    圖2 玉竹種群的UPGMA聚類圖Fig.2 Dendrogram obtained by UPGMA cluster analysis based on Nei’s genetic distance(D) among P.odoratum populations

    表4玉竹種群的遺傳距離(斜線下)與遺傳一致度(斜線上)

    Table4Geneticdistance(undertheobliqueline)andgeneticidentity(ontheobliqueline)amongP.odoratumpopulations

    WJDXHSXSSCDSDXPDGLDDCSXRDWJD????0.91340.93290.97310.91410.94980.95880.95330.9549XHS0.0906????0.93170.91680.88700.93270.91910.92470.9327XS0.06950.0707????0.93990.88760.96950.94560.95560.9474SCD0.02720.08690.0620????0.93770.95870.97030.95410.9662SD0.08980.11990.11920.0643????0.93030.92460.91680.9343XPD0.05150.06970.03100.04210.0723????0.96680.96390.9718GLD0.04210.08440.05590.03020.07840.0337????0.97020.9782DCS0.04780.07820.04540.04700.08690.03680.0303????0.9664XRD0.04620.06960.05400.03430.06790.02860.02200.0342????

    4 討論

    4.1 玉竹種群遺傳多樣性水平

    遺傳變異的程度直接表達(dá)了遺傳多樣性水平的高低,高水平的遺傳多樣性是維持物種長(zhǎng)期生存和進(jìn)化的基礎(chǔ)條件。遺傳多樣性越高或遺傳變異越豐富,植物適應(yīng)生存環(huán)境變化的能力就越強(qiáng),而遺傳多樣性的降低將會(huì)導(dǎo)致植物的適應(yīng)能力下降,甚至?xí)袦缃^的風(fēng)險(xiǎn)。本研究利用ISSR分子標(biāo)記揭示了大連地區(qū)9個(gè)野生玉竹種群的在物種水平上的多態(tài)性條帶比率(PPB)為91.67%,Nei’s基因多樣性指數(shù)(h)為0.346 0,Shannon信息指數(shù)(I)為0.510 8。說明了玉竹種群的遺傳多樣性較為豐富,這與卜靜等[6]通過對(duì)我國(guó)5個(gè)省份共11個(gè)種群的野生玉竹和1個(gè)種群的栽培玉竹進(jìn)行ISSR分析(PPB=94.13%,h=0.432 1,I=0.619 7)所得到的結(jié)論基本一致。

    本研究得到的玉竹遺傳多樣性較高的原因可能與玉竹自身的繁育系統(tǒng)有關(guān),植物的繁育系統(tǒng)是決定植物遺傳多樣性的最主要的因素[11]。玉竹的繁殖方式為有性繁殖和無性繁殖兩種方式并存。有性繁殖的遺傳信息來源于親本雙方,這樣更有利于基因的重組和突變,從而加大遺傳變異幾率,人們也普遍接受這一觀點(diǎn),但近年來研究[12]表明無性繁殖的種群也可以保持較高水平的遺傳多樣性。因而玉竹的兩種繁殖方式使得其遺傳多樣性水平較高。通過結(jié)果我們也可以看到XHS種群的多態(tài)性條帶比率(PPB)最低,為69.17%。這可能與強(qiáng)烈的人為干擾有關(guān),XHS種群采集地位于森林公園的人行道旁,經(jīng)常會(huì)有人在其種群分布范圍內(nèi)活動(dòng),強(qiáng)烈的人為干擾可能會(huì)導(dǎo)致其遺傳多樣性下降。人為干擾等外部因素通過某種間接方式(如影響植物的交配系統(tǒng)、瓶頸效應(yīng)等)使等位基因數(shù)目或頻率發(fā)生改變,最終使植物群體的遺傳多樣性水平和遺傳結(jié)構(gòu)發(fā)生變化[13]。

    4.2 玉竹種群的遺傳結(jié)構(gòu)

    本研究得到玉竹種群總體的遺傳分化系數(shù)為0.117 4,基因流為3.758 5。表明了玉竹種群間基因交流較為頻繁,種群間的分化程度較低?;蛄髦饕ㄟ^花粉擴(kuò)散和種子傳播兩種形式來實(shí)現(xiàn),玉竹種群之所以基因流較大,這可能與玉竹種群在大連地區(qū)分布較廣以及大陸種群及個(gè)別海島種群間地理距離較近有關(guān),XS與XPD種群地理距離較近,XRD和GLD之間,GLD與SCD、WJD之間的地理距離也較近,各島嶼間的人員流動(dòng)增大了基因交流的可能性,島嶼一般在玉竹的開花期風(fēng)比較大,加速了花粉的傳播,間接的促進(jìn)了基因的交流,島嶼的片段化生境并沒有完全阻止基因流。同時(shí)玉竹又是歐亞地區(qū)廣布種,在大連地區(qū)分布也較廣,也是導(dǎo)致玉竹具有較大基因流的原因。

    通過對(duì)遺傳距離和UPGMA聚類分析的對(duì)比可以發(fā)現(xiàn),XRD和GLD的遺傳距離為0.022 0,說明二者的親緣關(guān)系最近,而這兩個(gè)種群的地理距離較近,可以進(jìn)行基因交流,GLD雖然是一個(gè)島嶼,但其離岸距離很近,且有一條大壩與大陸相連,海水的隔離并不會(huì)對(duì)基因交流產(chǎn)生影響。SD與XHS的遺傳距離最大,為0.119 9,說明了二者的親緣關(guān)系最遠(yuǎn),有較大的遺傳背景差異,SD種群位于渤海海域,其它島嶼種群均來自黃海海域,較遠(yuǎn)的地理距離以及海水的阻隔阻止了與其他種群的基因交流,從而導(dǎo)致SD與其他種群的遺傳差異較大。

    UPGMA聚類圖可以直接反映出各種群的親緣關(guān)系,聚類結(jié)果與地理分布上具有一定的相關(guān)性,DCS種群與XRD、GLD的地理距離較遠(yuǎn),但也聚為一類,表明地理距離在玉竹天然種群的遺傳分化中可能并非唯一的主導(dǎo)因素,可能還受到周邊生境及人為干擾等影響,具體原因還有待于更多島嶼種群并結(jié)合玉竹的生殖生物學(xué)與生態(tài)學(xué)方面的研究作出科學(xué)的評(píng)價(jià)。

    4.3 玉竹海島種群與大陸種群遺傳多樣性的比較

    一般來說,一個(gè)物種的海島種群遺傳多樣性由于奠基者效應(yīng)或自然選擇等原因會(huì)低于大陸種群的遺傳多樣性[14]。分析本研究所得結(jié)果,發(fā)現(xiàn)在玉竹的9個(gè)種群中,從總體上看海島種群的遺傳多樣性略高于大陸種群。李麗等[15]利用SSR標(biāo)記技術(shù)對(duì)來自海南、廣東、廣西的19個(gè)美麗雞血藤野生種群的遺傳多樣性進(jìn)行比較,發(fā)現(xiàn)海南種質(zhì)較內(nèi)陸種質(zhì)擁有更高的遺傳變異,推測(cè)可能由于美麗雞血藤在島嶼生存面臨較之大陸更為復(fù)雜的環(huán)境,而進(jìn)化了較高的遺傳多樣性。分析玉竹形成海島種群的遺傳多樣性略高于大陸種群的原因,一方面可能與種群規(guī)模和人為干擾有關(guān),大陸種群中XRD種群采集地建有國(guó)家級(jí)自然保護(hù)區(qū),XPD種群位于密林深處,很少受到人為干擾,玉竹的棲息地也保存的較好,使得這兩個(gè)種群的遺傳多樣性相對(duì)較高;XHS種群受到人為干擾程度較大,在野外調(diào)查時(shí)還發(fā)現(xiàn)有人直接采挖玉竹根莖,對(duì)種群的破壞較大;XS種群規(guī)模較小。而海島種群雖然會(huì)受到地理隔離和生境片段化的影響,但是其種群規(guī)模較大,同時(shí)海島種群受到的人為干擾程度適中。導(dǎo)致大陸種群的遺傳多樣性差異較大且略低于海島種群。另一方面,海島種群也面臨著海島上更為復(fù)雜多樣的氣候及生存環(huán)境,多變的環(huán)境和孤立的生境使玉竹面臨更為復(fù)雜的選擇壓力,可能會(huì)影響玉竹進(jìn)化,導(dǎo)致海島玉竹種群可能會(huì)積累更多的遺傳變異從而形成較高的遺傳多樣性水平,用以適應(yīng)復(fù)雜的生態(tài)環(huán)境。

    1.中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì).中國(guó)植物志:第15卷[M].北京:科學(xué)出版社,1978.

    Editorial Committee Flora of China Chinese Academy of Sciences.Flora of China:vol.15[M].Beijing:Science Press,1978.

    2.潘清平,周日寶,陳玉秀,等.玉竹不同品種的ISSR分子鑒定[J].中國(guó)現(xiàn)代中藥,2008,10(10):28-30.

    Pan Q P,Zhou R B,Chen Y X,et al.Identification ofPolygonatumodoratumby ISSR[J].Modern Chinese Medicine,2008,10(10):28-30.

    3.竺平暉,陳愛萍.GC-MS法對(duì)湖南產(chǎn)玉竹揮發(fā)油成分的分析研究[J].中草藥,2010,41(8):1264-1265.

    Zhu P H,Chen A P.Analysis of volatile oils of PolygonatumOdoratumfrom Hunan by GC/MS[J].Chinese Traditional and Herbal Drugs,2010,41(8):1264-1265.

    4.Choi S B,Park S.A steroidal glycoside fromPolygonatumodoratum(Mill.) Druce.improves insulin resistance but does not alter insulin secretion in 90% pancreatectomized rats[J].Bioscience,Biotechnology,and Biochemistry,2002,66(10):2036-2043.

    5.Zietkiewicz E,Rafalski A,Labuda D.Genome fingerprinting by simple sequence repeat(SSR)-anchored polymerase chain reaction amplification[J].Genomics,1994,20(2):176-183.

    6.路放.玉竹種質(zhì)資源評(píng)價(jià)研究[D].長(zhǎng)春:吉林農(nóng)業(yè)大學(xué),2014.

    Lu F.Study on assessment of germplasm resources ofPolygonatumodoratum[D].Changchun:Jilin Agricultural University,2014.

    7.卜靜,王冬梅,李登武.不同產(chǎn)地野生玉竹種質(zhì)資源多樣性與親緣關(guān)系的ISSR分析[J].中草藥,2012,43(9):1824-1828.

    Bu J,Wang D M,Li D W.ISSR analysis on diversity and genetic relationship in germplasm resources of wildPo1ygonatumodoratumfrom different habitats[J].Chinese Traditional and Herbal Drugs,2012,43(9):1824-1828.

    8.張恒慶,安利佳,祖元?jiǎng)?涼水國(guó)家自然保護(hù)區(qū)天然紅松林遺傳變異的RAPD分析[J].植物研究,2000,20(2):201-206.

    Zhang H Q,An L J,Zu Y G.RAPD analysis on genetic variation ofPinuskoraiensisin Liang Shui national natural reserve[J].Bulletin of Botanical Research,2000,20(2):201-206.

    9.陳玉秀,周日寶,潘清平,等.玉竹ISSR反應(yīng)體系的優(yōu)化[J].湖南中醫(yī)藥大學(xué)學(xué)報(bào),2008,28(2):37-39.

    Chen Y X,Zhou R B,Pan Q P,et al.Optimization for reaction system ofPolygonatumodoratumISSR[J].Journal of TCM University of Hunan,2008,28(2):37-39.

    10.Yeh F C,Yang R C,Boyle T B,et al.POPGENE Version 1.32 Microsoft Windows-based freeware for populations genetic analysis.University of Alberta,Edmonton[M].Edmonton:Molecular Biology and Biotechnology Center,University of Alberta,1999.

    11.王洪新,胡志昂.植物的繁育系統(tǒng)、遺傳結(jié)構(gòu)和遺傳多樣性保護(hù)[J].生物多樣性,1996,4(2):92-96.

    Wang H X,Hu Z A.Plant breeding system,genetic structure and conservation of genetic diversity[J].Chinese Biodiversity,1996,4(2):92-96.

    12.Burd On R D,Aimers H J.Risk management for clonal forestry withPinusradiate-analysis and review.1:Strategic issues and risk spread[J].New Zealand Journal of forestry Science,2003,33(2):156-180.

    13.文亞峰,韓文軍,吳順.植物遺傳多樣性及其影響因素[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2010,30(12):80-87.

    Wen Y F,Han W J,Wu S.Plant genetic diversity and its influencing factors[J].Journal of Central South University of Forestry & Technology,2010,30(12):80-87.

    14.Muller F,Voccia M,B? A,et al.Genetic diversity and gene flow in a Caribbean treePterocarpusofficinalisJacq:a study based on chloroplast and nuclear microsatellites[J].Genetica,2009,135(2):185-198.

    15.李麗,徐立,李志英,等.不同地理來源美麗雞血藤遺傳多樣性的SSR分析[J].廣東農(nóng)業(yè)科學(xué),2013,40(22):156-160.

    Li L,Xu L,LI Z Y,et al.Genetic diversity ofCalleryaspeciosafrom different geographical origins by SSR markers[J].Guangdong Agricultural Sciences,2013,40(22):156-160.

    introduction:HAO Jiu-Cheng(1992—),male,Master,mainly engaged in population genetics.

    date:2017-03-04

    GeneticDiversitybetweenIslandandMainlandNaturalPopulationsofPolygonatumodoratuminDalianAreabyISSR

    HAO Jiu-Cheng JIA Xin MU Xiao-Hong ZHANG Heng-Qing*

    (College of Life Sciences,Liaoning Normal University,Dalian 116081)

    WithPolygonatumodoratum, a plant on both the mainland and islands in Dalian, we studied the genetic diversity of 262P.odoratumfrom 5 islands and 4 mainland populations by ISSR molecular marker technology. Information of 120 loci was obtained by amplification of 10 ISSR primers that were screened out, where the percentage of polymorphic bands was 91.67%, the gene diversity index of Nei’s was 0.346 0, and Shannon information index was 0.510 8. Its genetic differentiation coefficient was 0.117 4, and its gene flow was 3.758 5. The genetic diversity of natural populations ofP.odoratumwas relatively rich with relatively frequent gene exchange among different populations, and its genetic distance had certain correlation with its geographic distance. By comparison between genetic diversity of populations from the mainland and islands, the genetic diversity of populations from islands was slightly higher than that from the mainland, indicating that under a more isolated ecological environment and more complicated selection pressure, populations ofP.odoratumon islands might accumulate more genetic variation to form high levels of genetic diversity. Our results would provide evidence for further discussion on laws of genetic evolution of populations of natural plants in isolated ecological environments.

    ISS-PCR;Polygonatumodoratum;genetic diversity;genetic structure

    郝久程(1992—),男,碩士研究生,主要從事種群遺傳學(xué)方面的研究。

    * 通信作者:E-mail:hengqingzhang@163.com

    2017-03-04

    * Corresponding author:E-mail:hengqingzhang@163.com

    Q949.71+8.23

    A

    10.7525/j.issn.1673-5102.2017.05.010

    猜你喜歡
    玉竹海島種群
    邢氏水蕨成功繁衍并建立種群 等
    山西省發(fā)現(xiàn)刺五加種群分布
    養(yǎng)陰潤(rùn)燥話玉竹
    冰與火共存的海島
    奧秘(2020年5期)2020-06-30 10:12:10
    編竹簍
    在海島度假
    玉竹的化學(xué)成分研究
    玉竹
    崗更湖鯉魚的種群特征
    神奇的海島
    精品少妇久久久久久888优播| 午夜91福利影院| 亚洲一区中文字幕在线| 欧美av亚洲av综合av国产av| 成人免费观看视频高清| 亚洲色图 男人天堂 中文字幕| www.熟女人妻精品国产| 国产成+人综合+亚洲专区| 国产精品久久久久久精品古装| 久久精品aⅴ一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品免费一区二区三区在线 | 丝袜人妻中文字幕| 成人18禁高潮啪啪吃奶动态图| av不卡在线播放| 亚洲少妇的诱惑av| 国产男女超爽视频在线观看| av片东京热男人的天堂| 久久精品亚洲精品国产色婷小说| 日本欧美视频一区| 久久精品国产亚洲av高清一级| 日本撒尿小便嘘嘘汇集6| 91精品国产国语对白视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲中文av在线| 如日韩欧美国产精品一区二区三区| 久久久久精品国产欧美久久久| 欧美日韩av久久| 久久精品熟女亚洲av麻豆精品| 天堂中文最新版在线下载| 菩萨蛮人人尽说江南好唐韦庄| 我的亚洲天堂| 欧美日韩亚洲高清精品| 国产成人精品在线电影| 国产97色在线日韩免费| 美国免费a级毛片| 亚洲伊人色综图| 国产亚洲欧美在线一区二区| 亚洲伊人久久精品综合| 99国产精品一区二区三区| 国产男靠女视频免费网站| 90打野战视频偷拍视频| 日韩中文字幕视频在线看片| 岛国在线观看网站| 久久久欧美国产精品| 亚洲精华国产精华精| 亚洲国产中文字幕在线视频| 黄网站色视频无遮挡免费观看| 另类亚洲欧美激情| 后天国语完整版免费观看| 青青草视频在线视频观看| 亚洲精品在线美女| 国产精品.久久久| 精品国产乱码久久久久久小说| 脱女人内裤的视频| 午夜免费成人在线视频| 国产精品av久久久久免费| 免费不卡黄色视频| 我要看黄色一级片免费的| 91成人精品电影| 香蕉丝袜av| 一本—道久久a久久精品蜜桃钙片| 香蕉久久夜色| 成年动漫av网址| 天天操日日干夜夜撸| 一级片'在线观看视频| 岛国在线观看网站| 亚洲国产毛片av蜜桃av| 19禁男女啪啪无遮挡网站| 最黄视频免费看| xxxhd国产人妻xxx| av天堂在线播放| 大陆偷拍与自拍| 国产成人免费观看mmmm| 国产亚洲精品久久久久5区| 亚洲国产av影院在线观看| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 久久精品人人爽人人爽视色| 久久这里只有精品19| 精品卡一卡二卡四卡免费| 人人妻,人人澡人人爽秒播| 一个人免费在线观看的高清视频| 91麻豆av在线| 亚洲欧美日韩另类电影网站| 捣出白浆h1v1| 日韩有码中文字幕| 中文字幕精品免费在线观看视频| 久热这里只有精品99| 18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 超碰97精品在线观看| 99re6热这里在线精品视频| 别揉我奶头~嗯~啊~动态视频| 国精品久久久久久国模美| 男人舔女人的私密视频| 男女之事视频高清在线观看| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 久久中文看片网| 欧美乱妇无乱码| 国产99久久九九免费精品| 亚洲专区中文字幕在线| 国产免费福利视频在线观看| 国产一区二区 视频在线| 在线观看www视频免费| 亚洲成av片中文字幕在线观看| 成人黄色视频免费在线看| 日日摸夜夜添夜夜添小说| 老熟女久久久| 啪啪无遮挡十八禁网站| 亚洲 欧美一区二区三区| 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色 | 丁香六月天网| 久久精品91无色码中文字幕| 中文字幕av电影在线播放| 日韩欧美国产一区二区入口| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲日产国产| 欧美日韩亚洲高清精品| 国产伦人伦偷精品视频| 最近最新中文字幕大全电影3 | 日韩视频在线欧美| 视频在线观看一区二区三区| 精品高清国产在线一区| 久久亚洲精品不卡| 一区二区三区精品91| 亚洲精品av麻豆狂野| 极品少妇高潮喷水抽搐| 这个男人来自地球电影免费观看| 黄色怎么调成土黄色| 亚洲全国av大片| 香蕉丝袜av| 亚洲国产中文字幕在线视频| 亚洲专区中文字幕在线| 男女免费视频国产| 搡老乐熟女国产| 亚洲国产欧美日韩在线播放| 亚洲av第一区精品v没综合| √禁漫天堂资源中文www| 日本黄色视频三级网站网址 | 国产欧美日韩综合在线一区二区| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 777米奇影视久久| 一二三四在线观看免费中文在| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱久久久久久| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 久久午夜综合久久蜜桃| 国产精品影院久久| 桃红色精品国产亚洲av| 老司机深夜福利视频在线观看| 国产av一区二区精品久久| 青草久久国产| 国产精品熟女久久久久浪| 久久毛片免费看一区二区三区| 一区在线观看完整版| 搡老熟女国产l中国老女人| 日韩精品免费视频一区二区三区| 中文欧美无线码| 99久久99久久久精品蜜桃| av欧美777| 制服诱惑二区| 欧美日韩av久久| 国产欧美日韩精品亚洲av| 99久久精品国产亚洲精品| 水蜜桃什么品种好| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 亚洲欧美精品综合一区二区三区| 国产男女超爽视频在线观看| 欧美久久黑人一区二区| 色在线成人网| 在线观看免费午夜福利视频| 精品福利观看| 久久天堂一区二区三区四区| 这个男人来自地球电影免费观看| 王馨瑶露胸无遮挡在线观看| 精品乱码久久久久久99久播| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 国产精品一区二区在线观看99| 欧美另类亚洲清纯唯美| 老熟妇仑乱视频hdxx| 国产免费福利视频在线观看| 精品少妇久久久久久888优播| 少妇精品久久久久久久| 99精品久久久久人妻精品| 久久精品亚洲熟妇少妇任你| 亚洲,欧美精品.| 亚洲人成伊人成综合网2020| 精品福利观看| 国产精品九九99| 欧美一级毛片孕妇| 热re99久久精品国产66热6| 亚洲九九香蕉| 国产高清国产精品国产三级| 亚洲国产毛片av蜜桃av| 一进一出抽搐动态| 日本黄色视频三级网站网址 | 精品欧美一区二区三区在线| 999久久久国产精品视频| 国产成人精品无人区| 精品少妇久久久久久888优播| 国产在线观看jvid| 免费观看a级毛片全部| 不卡av一区二区三区| 一二三四社区在线视频社区8| 国产成人欧美在线观看 | 日本精品一区二区三区蜜桃| 国产男女内射视频| 狠狠精品人妻久久久久久综合| av超薄肉色丝袜交足视频| 国产精品亚洲一级av第二区| 美女午夜性视频免费| 日本av免费视频播放| 少妇猛男粗大的猛烈进出视频| av在线播放免费不卡| 麻豆国产av国片精品| 国产精品熟女久久久久浪| www.999成人在线观看| 精品欧美一区二区三区在线| 老熟妇仑乱视频hdxx| 亚洲国产欧美一区二区综合| 亚洲专区字幕在线| 色94色欧美一区二区| 黑人猛操日本美女一级片| 日韩制服丝袜自拍偷拍| 亚洲五月婷婷丁香| 午夜激情av网站| 国产精品一区二区在线观看99| www.熟女人妻精品国产| 午夜福利欧美成人| 老司机午夜十八禁免费视频| 亚洲成人免费电影在线观看| 母亲3免费完整高清在线观看| 少妇猛男粗大的猛烈进出视频| 国产精品久久久人人做人人爽| 性少妇av在线| 好男人电影高清在线观看| 久久精品成人免费网站| 在线十欧美十亚洲十日本专区| 国产欧美亚洲国产| 欧美久久黑人一区二区| 他把我摸到了高潮在线观看 | 国产欧美日韩一区二区三| 国产区一区二久久| 久久久精品区二区三区| 一边摸一边抽搐一进一小说 | 亚洲国产毛片av蜜桃av| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 亚洲欧美一区二区三区黑人| 国产免费视频播放在线视频| 久久中文看片网| 777米奇影视久久| 亚洲 国产 在线| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区久久| 亚洲av日韩精品久久久久久密| 国产高清videossex| 国产一区有黄有色的免费视频| 精品国产乱码久久久久久小说| 女同久久另类99精品国产91| 交换朋友夫妻互换小说| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 精品卡一卡二卡四卡免费| 一本一本久久a久久精品综合妖精| 精品午夜福利视频在线观看一区 | 制服诱惑二区| 九色亚洲精品在线播放| 国产成人精品无人区| 精品卡一卡二卡四卡免费| 中文字幕人妻丝袜一区二区| 国产精品久久电影中文字幕 | 在线永久观看黄色视频| 亚洲精品国产一区二区精华液| 十八禁高潮呻吟视频| 一区在线观看完整版| 精品少妇一区二区三区视频日本电影| 久久人妻熟女aⅴ| 中文字幕高清在线视频| 岛国毛片在线播放| xxxhd国产人妻xxx| www.999成人在线观看| 男女高潮啪啪啪动态图| 两个人免费观看高清视频| 亚洲精品乱久久久久久| 午夜福利视频在线观看免费| 极品人妻少妇av视频| 999精品在线视频| 麻豆av在线久日| 不卡av一区二区三区| 国产成人精品在线电影| 亚洲av片天天在线观看| 性高湖久久久久久久久免费观看| 久久青草综合色| 欧美午夜高清在线| 色尼玛亚洲综合影院| 日韩成人在线观看一区二区三区| 免费不卡黄色视频| 精品亚洲成a人片在线观看| 69精品国产乱码久久久| 欧美亚洲日本最大视频资源| 成年人黄色毛片网站| 国产精品1区2区在线观看. | 别揉我奶头~嗯~啊~动态视频| 国产视频一区二区在线看| 在线看a的网站| 精品久久久精品久久久| 亚洲性夜色夜夜综合| 国产1区2区3区精品| 深夜精品福利| 91老司机精品| 亚洲久久久国产精品| 搡老乐熟女国产| 亚洲精品久久成人aⅴ小说| 亚洲精品乱久久久久久| 在线观看一区二区三区激情| 欧美一级毛片孕妇| 国产又色又爽无遮挡免费看| 国产男女内射视频| 久热这里只有精品99| 一级片'在线观看视频| 国产一区有黄有色的免费视频| 日本vs欧美在线观看视频| 狠狠精品人妻久久久久久综合| 国产黄色免费在线视频| 美女福利国产在线| 国产精品.久久久| 亚洲色图综合在线观看| 国产精品九九99| 日本av手机在线免费观看| 午夜福利,免费看| 中文字幕最新亚洲高清| 亚洲精品一二三| √禁漫天堂资源中文www| 久久狼人影院| 色尼玛亚洲综合影院| 久久精品熟女亚洲av麻豆精品| 精品欧美一区二区三区在线| 精品视频人人做人人爽| 成人18禁在线播放| 久久久欧美国产精品| 成人18禁在线播放| 国产精品免费一区二区三区在线 | 色综合婷婷激情| videosex国产| 亚洲国产欧美网| 国产精品.久久久| 一级毛片精品| 露出奶头的视频| 国产99久久九九免费精品| 一二三四社区在线视频社区8| 91老司机精品| 亚洲精品国产精品久久久不卡| av国产精品久久久久影院| 亚洲精品美女久久av网站| 精品国内亚洲2022精品成人 | 免费观看av网站的网址| 一本综合久久免费| 亚洲欧美一区二区三区黑人| 9热在线视频观看99| 亚洲黑人精品在线| 一级黄色大片毛片| 一边摸一边抽搐一进一小说 | 亚洲男人天堂网一区| 9色porny在线观看| 国产午夜精品久久久久久| 丁香六月欧美| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 亚洲精品中文字幕一二三四区 | 91国产中文字幕| 黄网站色视频无遮挡免费观看| 国产精品成人在线| 日韩免费av在线播放| 日韩视频一区二区在线观看| 欧美乱码精品一区二区三区| 久久精品国产亚洲av高清一级| 在线 av 中文字幕| 色精品久久人妻99蜜桃| 动漫黄色视频在线观看| 精品熟女少妇八av免费久了| 午夜免费成人在线视频| 久久久久视频综合| 精品人妻在线不人妻| 男女高潮啪啪啪动态图| 国产精品二区激情视频| 久久 成人 亚洲| 国产精品.久久久| 99久久99久久久精品蜜桃| 国内毛片毛片毛片毛片毛片| 在线 av 中文字幕| 啪啪无遮挡十八禁网站| 黄网站色视频无遮挡免费观看| 国产色视频综合| 美女视频免费永久观看网站| 99久久国产精品久久久| 免费黄频网站在线观看国产| 国产激情久久老熟女| 久久精品国产亚洲av香蕉五月 | 九色亚洲精品在线播放| 国产色视频综合| 成人av一区二区三区在线看| 99热国产这里只有精品6| 嫩草影视91久久| 成人国产av品久久久| 在线av久久热| 日韩 欧美 亚洲 中文字幕| 欧美国产精品一级二级三级| 日本精品一区二区三区蜜桃| 热re99久久国产66热| 欧美日韩亚洲综合一区二区三区_| 欧美人与性动交α欧美软件| 一级,二级,三级黄色视频| 一区二区三区激情视频| 免费在线观看黄色视频的| 亚洲伊人色综图| 可以免费在线观看a视频的电影网站| 久久狼人影院| 黄片播放在线免费| 午夜福利在线免费观看网站| 岛国在线观看网站| 大香蕉久久成人网| 狠狠精品人妻久久久久久综合| 日本vs欧美在线观看视频| 人妻久久中文字幕网| 黄网站色视频无遮挡免费观看| 男人操女人黄网站| av天堂在线播放| 91麻豆av在线| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产免费现黄频在线看| 欧美人与性动交α欧美软件| 免费看十八禁软件| 老熟妇仑乱视频hdxx| 操美女的视频在线观看| 在线观看免费视频网站a站| 搡老乐熟女国产| 国产午夜精品久久久久久| 国产亚洲精品一区二区www | 久久久精品94久久精品| www.999成人在线观看| 亚洲成人手机| 国产成人精品无人区| 女同久久另类99精品国产91| 操出白浆在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧洲精品一区二区精品久久久| 99国产精品免费福利视频| 精品卡一卡二卡四卡免费| 最近最新中文字幕大全免费视频| 成人特级黄色片久久久久久久 | 日韩中文字幕欧美一区二区| 午夜老司机福利片| 亚洲第一av免费看| 久久人妻av系列| 国产午夜精品久久久久久| 精品久久久精品久久久| 午夜福利在线观看吧| 老汉色∧v一级毛片| 性少妇av在线| 国产精品久久久人人做人人爽| 精品少妇久久久久久888优播| 大陆偷拍与自拍| 欧美日韩精品网址| 午夜精品久久久久久毛片777| 亚洲伊人色综图| 国产精品免费大片| 久久久水蜜桃国产精品网| 咕卡用的链子| 久久久久国内视频| 丝瓜视频免费看黄片| 久久久久视频综合| 欧美黄色片欧美黄色片| 在线观看舔阴道视频| 亚洲黑人精品在线| 国产高清视频在线播放一区| av又黄又爽大尺度在线免费看| 亚洲精品一二三| 亚洲中文日韩欧美视频| 老鸭窝网址在线观看| 亚洲国产欧美日韩在线播放| 老司机午夜福利在线观看视频 | 高清av免费在线| 中文字幕人妻丝袜一区二区| 亚洲专区字幕在线| 在线观看免费日韩欧美大片| 国产男靠女视频免费网站| 黑人猛操日本美女一级片| 国产精品九九99| 在线 av 中文字幕| 69精品国产乱码久久久| 久久国产精品男人的天堂亚洲| 国产精品一区二区在线不卡| 91成年电影在线观看| 欧美日韩精品网址| 日本一区二区免费在线视频| 一区福利在线观看| 一级毛片精品| 啪啪无遮挡十八禁网站| 狠狠精品人妻久久久久久综合| 91麻豆av在线| 国产欧美日韩精品亚洲av| 国产精品久久久久成人av| 久久性视频一级片| 精品亚洲成a人片在线观看| 国产人伦9x9x在线观看| 亚洲av电影在线进入| 国产精品1区2区在线观看. | √禁漫天堂资源中文www| 欧美激情极品国产一区二区三区| 国产精品国产高清国产av | 国产免费视频播放在线视频| 人成视频在线观看免费观看| 欧美 日韩 精品 国产| 婷婷丁香在线五月| 十八禁网站免费在线| 国产精品久久久久成人av| 两性夫妻黄色片| 亚洲国产毛片av蜜桃av| 高清av免费在线| kizo精华| 日韩免费高清中文字幕av| 国产精品偷伦视频观看了| 色在线成人网| 人人妻人人添人人爽欧美一区卜| 亚洲七黄色美女视频| 亚洲国产欧美日韩在线播放| 2018国产大陆天天弄谢| 国产精品.久久久| 97人妻天天添夜夜摸| 视频区图区小说| 两个人免费观看高清视频| 精品高清国产在线一区| 久久精品成人免费网站| 国产亚洲精品第一综合不卡| 欧美久久黑人一区二区| 一夜夜www| 12—13女人毛片做爰片一| 91麻豆av在线| 国产三级黄色录像| 精品福利观看| 老司机午夜十八禁免费视频| 亚洲第一青青草原| 两性午夜刺激爽爽歪歪视频在线观看 | a级片在线免费高清观看视频| 99久久99久久久精品蜜桃| 欧美亚洲 丝袜 人妻 在线| 免费黄频网站在线观看国产| 精品国产乱子伦一区二区三区| 精品国产乱码久久久久久男人| 国产亚洲一区二区精品| 老司机午夜福利在线观看视频 | 最近最新中文字幕大全电影3 | 国产色视频综合| 日本欧美视频一区| 麻豆av在线久日| 精品人妻1区二区| 最近最新中文字幕大全电影3 | 亚洲七黄色美女视频| 欧美日韩亚洲高清精品| 色尼玛亚洲综合影院| 精品国产一区二区三区久久久樱花| 国产日韩欧美视频二区| 久久久久视频综合| 啦啦啦 在线观看视频| 九色亚洲精品在线播放| 视频区图区小说| 99re在线观看精品视频| 正在播放国产对白刺激| 国产欧美亚洲国产| 中文字幕精品免费在线观看视频| 久久精品91无色码中文字幕| 精品亚洲成国产av| 精品国产乱子伦一区二区三区| 午夜老司机福利片| 久久中文看片网| 757午夜福利合集在线观看| 欧美日韩福利视频一区二区| 国产精品国产av在线观看| 一本大道久久a久久精品| 欧美日韩一级在线毛片| 精品人妻在线不人妻| 一本久久精品| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费视频网站a站| 亚洲欧美一区二区三区久久| 亚洲熟女精品中文字幕| 国产av一区二区精品久久| 日韩免费高清中文字幕av| 欧美日韩一级在线毛片| 宅男免费午夜| 岛国在线观看网站| 老司机午夜十八禁免费视频| 99国产精品免费福利视频| 国产精品99久久99久久久不卡| av片东京热男人的天堂| 99在线人妻在线中文字幕 | 午夜福利,免费看| 国产精品国产高清国产av | 国产成+人综合+亚洲专区| 99精品久久久久人妻精品| 91麻豆精品激情在线观看国产 | 天天躁日日躁夜夜躁夜夜| 久久香蕉激情| 精品一品国产午夜福利视频| 丰满少妇做爰视频|