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

    21種野生百合親緣關(guān)系的ISSR分析

    2016-10-28 01:39:47羅建讓肖菲菲張延龍牛立新
    廣西植物 2016年9期
    關(guān)鍵詞:野百合親緣百合

    侯 珺, 羅建讓, 肖菲菲, 張延龍, 牛立新

    ( 西北農(nóng)林科技大學(xué) 風(fēng)景園林藝術(shù)學(xué)院, 陜西 楊凌 712100 )

    ?

    21種野生百合親緣關(guān)系的ISSR分析

    侯珺, 羅建讓, 肖菲菲, 張延龍*, 牛立新

    ( 西北農(nóng)林科技大學(xué) 風(fēng)景園林藝術(shù)學(xué)院, 陜西 楊凌 712100 )

    該研究選用7條ISSR標(biāo)記引物對(duì)百合屬21種野生百合的親緣關(guān)系進(jìn)行分析,使用POPGEN1.32和MEGA5.1數(shù)據(jù)處理軟件分析數(shù)據(jù)。結(jié)果表明:7條引物擴(kuò)增出149條條帶,其中136條為多態(tài)性條帶,多態(tài)性條帶比率為91.06%,平均有效等位點(diǎn)數(shù)為1.762 4,平均Nei’s基因多樣度為0.421 4,平均Shanon信息指數(shù)0.608 5。遺傳距離變化范圍為0.307 5~0.887 3,紫脊百合和尖被百合的遺傳距離最大,均值為0.887 3,表明21份材料中二者親緣關(guān)系最遠(yuǎn);紫脊百合和宜昌百合遺傳距離最小,均值為0.307 5,表明二者親緣關(guān)系最近。聚類結(jié)果與形態(tài)學(xué)分類大致吻合,21個(gè)供試材料可被分成5個(gè)類群,大理百合、寶興百合、卷丹百合、紫斑百合、乳頭百合、川百合、蘭州百合、山丹百合、綠花百合和毛百合為第Ⅰ類群;紫脊百合、宜昌百合、岷江百合、通江百合和淡黃花百合為第Ⅱ類群;野百合和百合為第Ⅲ類群;青島百合為第Ⅳ類群;第Ⅴ類群包括玫紅百合、有斑百合和尖被百合。毛百合與卷瓣組百合親緣關(guān)系較近,鐘花組百合與卷瓣組百合存在基因交流,說明不能僅僅依靠形態(tài)學(xué)對(duì)野生百合進(jìn)行分類。聚類分析結(jié)果中野百合與百合單獨(dú)聚為一類,這說明是否具葉柄對(duì)野生百合的分類是一個(gè)重要的形態(tài)學(xué)特性。ISSR分子標(biāo)記適合用于百合屬植物的親緣關(guān)系分析。

    百合屬, 野生百合, ISSR, 親緣關(guān)系, 樹狀圖

    百合屬(Lilium)為百合科百合族的一個(gè)多年生屬,有90余種,大多分布在北半球的溫帶和寒帶地區(qū),少數(shù)分布在熱帶高海拔山區(qū),中國(guó)為世界百合起源中心,約46種18個(gè)變種(趙祥云等,1999),分布于27個(gè)省市自治區(qū)(陳心啟等,1980)。百合屬植物具多年生地下鱗莖,花色多樣,藥用及觀賞價(jià)值極高。雖然野生百合種質(zhì)資源豐富,但對(duì)百合屬的屬下類群劃分、起源、進(jìn)化等問題仍爭(zhēng)議頗多?,F(xiàn)國(guó)際最常采用的分類系統(tǒng)是Comber(1949)依據(jù)15個(gè)形態(tài)特征建立的,他將百合屬分為7組。《中國(guó)植物志》將百合屬植物劃分為4個(gè)組(陳心啟等,1980)。Nishikawa et al(1999)最早對(duì)百合屬進(jìn)行分子系統(tǒng)學(xué)研究,發(fā)現(xiàn)國(guó)際分類系統(tǒng)的毛百合組與卷瓣組親緣關(guān)系較近,國(guó)際分類系統(tǒng)中屬于毛百合組的毛百合被劃分在中國(guó)分類系統(tǒng)中的鐘花組,由此可見毛百合的系統(tǒng)位置較復(fù)雜。我國(guó)百合屬分組中屬于百合組的野百合和百合,在Comber的分組中被劃分在具葉柄組(sect.Archelirion)。中國(guó)是百合屬植物的自然分布中心之一,因此研究中國(guó)百合屬植物親緣關(guān)系對(duì)于探討百合屬起源、進(jìn)化、遷移以及建立科學(xué)的分類系統(tǒng)至關(guān)重要。目前國(guó)內(nèi)百合系統(tǒng)研究多以形態(tài)學(xué)(向地英等,2006)、染色體特征(張藝萍等,2014)、花粉形態(tài)(顧欣等,2013)為依據(jù)。利用分子標(biāo)記研究百合親緣關(guān)系多見于栽培品種,原生種的分子標(biāo)記有RAPD,AFLP,SRAP分析,較全面的針對(duì)野生種系統(tǒng)學(xué)的研究較少。

    簡(jiǎn)單重復(fù)序列間擴(kuò)散多態(tài)性(ISSR)是由加拿大蒙特利爾大學(xué)Ziekiewicz et al(1994)基于SSR標(biāo)記創(chuàng)建來的一種新型分子標(biāo)記技術(shù),其基本原理就是在SSR的3′或5′端加錨1~4個(gè)隨機(jī)核苷酸作為引物,對(duì)兩側(cè)具有反向排列SSR的一段DNA序列進(jìn)行擴(kuò)增,然后進(jìn)行電泳、染色,根據(jù)所讀取譜帶的有無(wú)及相對(duì)位置差異,來分析不同樣品間 ISSR標(biāo)記的多態(tài)性。由于 ISSR標(biāo)記技術(shù)結(jié)合了RAPD和SSR技術(shù)的優(yōu)點(diǎn),實(shí)驗(yàn)操作簡(jiǎn)便快捷、重復(fù)性高、穩(wěn)定性好、成本較低、所揭示的遺傳多態(tài)性十分豐富,因此,ISSR技術(shù)已被廣泛應(yīng)用于植物育種、植物遺傳多樣性、系統(tǒng)發(fā)育、品種鑒定、基因定位、遺傳作圖等研究中(Tsumura et al,1996;Qian et al,2011;Pradeep et al,2002;Devarumath et al,2002;Guo et al,2011)。 本研究將此項(xiàng)技術(shù)應(yīng)用于中國(guó)百合屬的親緣關(guān)系分析,首次對(duì)我國(guó)百合屬四個(gè)組21種野生類群進(jìn)行研究,希望為百合屬系統(tǒng)學(xué)研究及資源合理利用提供科學(xué)依據(jù),并對(duì)分類上存在爭(zhēng)議的種的歸屬進(jìn)行確定。

    1 材料與方法

    1.1 材料

    供試的野生百合材料有21個(gè)種或變種,來自陜西、甘肅、云南、四川和山東省(表1),由西北農(nóng)林科技大學(xué)牛立新教授鑒定,種植于西北農(nóng)林科技大學(xué)百合資源溫室內(nèi)。

    1.2 方法

    1.2.1 DNA提取采用Biospin植物基因組DNA提取試劑盒提取21份植物材料的基因組DNA。用1%的瓊脂糖凝膠電泳檢測(cè)DNA的純度與濃度,稀釋至80 ng·μL-1于-20 ℃保存。

    表 1 供試21種野生百合采集地

    表 2 7條ISSR引物的退火溫度及其擴(kuò)增結(jié)果

    Y = (C, T)

    1.2.2 引物篩選及優(yōu)化體系建立基于早期ISSR標(biāo)記在木蘭亞屬植物應(yīng)用中的研究報(bào)道,本研究引物參照加拿大哥倫比亞大學(xué)(UBC)公布的序列和日本Yamagishi et al(2002)在研究百合所采用的3A引物序列,由上海生工生物工程技術(shù)服務(wù)有限公司合成,利用Eppendorf梯度PCR擴(kuò)增儀(Mastercycler pro型)確定引物最佳退火溫度和篩選引物,最終從163條引物中篩選出7條穩(wěn)定且重復(fù)性好的ISSR引物。利用單因素設(shè)計(jì),對(duì)影響PCR反應(yīng)的模板濃度、引物濃度、循環(huán)次數(shù)等進(jìn)行優(yōu)化,建立適合野生百合ISSR-PCR最佳反應(yīng)體系。

    PCR總反應(yīng)體系為20 μL,含10 μL 2 ×TaqPCR Master Mix,RNase-free water 8 μL ,DNA模板80 ng,10 μmol·L-1引物1 μL。PCR擴(kuò)增反應(yīng)程序?yàn)?4 ℃預(yù)變性5 min,94 ℃變性40 s,51.1~58.1 ℃(因引物不同而異)退火45 s,72 ℃ 延伸1.5 min,39個(gè)循環(huán),最后72 ℃延伸8 min,4 ℃ 保存。PCR產(chǎn)物在6%聚丙烯酰胺凝膠上電泳,銀染檢測(cè)。

    1.2.3 數(shù)據(jù)統(tǒng)計(jì)分析根據(jù)擴(kuò)增產(chǎn)物電泳圖譜,對(duì)清晰且易于辨認(rèn)的條帶,在相同遷移位置上有帶標(biāo)為“1”,無(wú)帶標(biāo)為“0”,構(gòu)建出樣品分子標(biāo)記0、1矩陣。應(yīng)用POPGEN 1.32軟件,采用SM相似系數(shù)法計(jì)算21份野生百合材料的多態(tài)性條帶百分比(PPB)、有效等位基因樹(Ne)、Nei’s基因多樣度(H)、Shanon信息指數(shù)(I)、和遺傳距離。應(yīng)用MEGA5.1軟件,采用鄰接法(Neighbor-joining)進(jìn)行聚類分析,得出野生百合親緣關(guān)系聚類圖。

    2 結(jié)果與分析

    2.1 ISSR擴(kuò)增結(jié)果

    選用野生百合四個(gè)組的四種野生百合材料,分別為輪葉組的青島百合、鐘花組的玫紅百合、卷瓣組的卷丹百合和百合組的宜昌百合,篩選ISSR引物,從169條引物中篩選出7條擴(kuò)增效果較好的引物,用于野生百合材料的ISSR擴(kuò)增。多態(tài)性條帶的分子量主要分布在250~750 bp之間,擴(kuò)增結(jié)果及多態(tài)性信息見表2,共得到清晰條帶149條,多態(tài)性條帶136條,多態(tài)性條帶比率為91.06%。經(jīng)重復(fù)試驗(yàn)發(fā)現(xiàn) ISSR標(biāo)記穩(wěn)定性較強(qiáng),多態(tài)性條帶豐富,利用7條引物能有效地將所有樣品進(jìn)行區(qū)分,因而ISSR可用于野生百合種質(zhì)的鑒定。其中引物UBC842對(duì) 21 份百合材料的擴(kuò)增結(jié)果見圖 1 。

    2.2 親緣關(guān)系分析

    21份供試材料的遺傳距離變化范圍為0.307 5~0.887 3。ISSR引物組合平均有效等位基因數(shù)、Nei’s 基因多樣度和 Shannon 信息指數(shù)均值分別為1.762 4,0.421 4和0.608 5(表2)。表2表明供試材料多態(tài)性高,遺傳變異豐富。數(shù)據(jù)統(tǒng)計(jì)結(jié)果顯示,紫脊百合和尖被百合的遺傳距離最大,均值為0.8873,表明21份野生百合材料中二者親緣關(guān)系最遠(yuǎn);紫脊百合和宜昌百合遺傳距離最小,均值為0.3075,表明二者親緣關(guān)系最近(表3)。

    圖 1 引物UBC842對(duì)21個(gè)野生百合樣品擴(kuò)增的凝膠電泳圖譜 1. 有斑百合; 2. 尖被百合; 3. 野百合; 4. 山丹百合; 5. 綠花百合; 6. 通江百合; 7. 川百合; 8. 淡黃花百合; 9. 紫脊百合; 10. 宜昌百合; 11. 岷江百合; 12. 百合; 13. 玫紅百合; 14. 蘭州百合; 15. 青島百合; 16. 乳頭百合; 17. 大理百合; 18. 紫斑百合; 19. 毛百合; 20. 寶興百合; 21. 卷丹百合。Fig. 1 Electrophoresis result of amplification of 21 wild species of Lilium using primer U842 1. L. concolor var. pulchellum; 2. L. lophophorum; 3. L. brownie; 4. L. pumilum; 5. L. fargesii; 6. L. sargentiae; 7. L. davidii; 8. L. sulphureum; 9. L. leucanthum var. centifolium; 10. L. leucanthum; 11. L. regale; 12. L. brownii var. viridulum; 13. L. amoenum; 14. L. davidii var. unicolor; 15. L. tsingtauense; 16. L. papilliferum; 17. L. taliense; 18. L. nepalense; 19. L. dauricum; 20. L. duchartrei; 21. L. lancifolium.

    2102000.53061900.49380.40181800.44670.58190.55591700.43530.44670.48180.48181600.49380.54320.55590.62220.56881500.72300.65000.59520.66420.65000.62221400.62220.49380.50590.43530.51820.55590.50591300.51820.58190.69310.65000.65000.58190.62220.67851200.49380.43530.62220.59520.60860.55590.59520.60860.63601100.67850.50590.47000.58190.58190.56880.54320.53060.54320.54321000.43530.59520.63600.59520.63600.66420.56880.56880.60860.51820.4938900.30750.45830.59520.75380.42400.63600.55590.49380.47000.58190.59520.5952800.51820.59520.40180.63600.73830.58190.62220.59520.63600.58190.70800.63600.6360700.56880.48180.72300.48180.49380.58190.38010.58190.50590.47000.54320.50590.67850.5688600.60860.44670.55590.60860.53060.59520.63600.56880.66420.69310.62220.42400.63600.54320.6785500.54320.54320.50590.56880.65000.56880.60860.59520.50590.54320.51820.55590.50590.54320.58190.5559400.42400.60860.60860.54320.58190.63600.53060.59520.66420.56880.75380.55590.54320.51820.48180.70800.4938300.55590.54320.55590.72300.56880.53060.53060.55590.47000.58190.56880.66420.69310.62220.59520.58190.51820.6500200.78550.75380.65000.72300.78550.65000.88730.66420.63600.65000.58190.62220.72300.63600.80180.80180.60860.67850.6500100.41280.81830.75380.62220.75380.58190.59520.75380.63600.60860.70800.53060.62220.69310.63600.67850.80180.66420.70800.8018123456789101112131415161718192021

    注: 表內(nèi)野生百合順序同圖1。

    Note: see Figure 1 for the sequence ofLilium.

    圖 2 21種野生百合聚類圖Fig. 2 Cluster analysis of 21 wild species of Lilium based on ISSR data using UPGMA

    2.3 聚類分析

    根據(jù)ISSR擴(kuò)增結(jié)果,得到21份野生百合材料的親緣關(guān)系樹狀圖(圖2)。由圖2可見, 21個(gè)供試材料可分成5類:Ⅰ類包括卷瓣組的大理百合、寶興百合、卷丹百合、紫斑百合、乳頭百合、川百合、蘭州百合、山丹百合和綠花百合,鐘花組的毛百合;Ⅱ類包括百合組的紫脊百合、宜昌百合、岷江百合、通江百合和淡黃花百合;Ⅲ類包括百合組的野百合、百合;Ⅳ類有輪葉組的青島百合;Ⅴ類有鐘花組的玫紅百合、有斑百合和尖被百合。

    3 討論與結(jié)論

    本研究采用7個(gè)ISSR標(biāo)記引物對(duì)原產(chǎn)中國(guó)的21個(gè)百合野生種基因組DNA進(jìn)行擴(kuò)增,得到清晰條帶149條,多態(tài)性條帶136條,多態(tài)性比率為91.06%。統(tǒng)計(jì)結(jié)果顯示供試材料間的遺傳距離較大,這說明我國(guó)野生百合種間差異明顯、基因資源豐富。由此可見,ISSR標(biāo)記適合于百合屬植物親緣關(guān)系分析,是一種有效、可靠的分子標(biāo)記。

    21份材料的遺傳距離變化范圍為0.307 5~0.887 3。數(shù)據(jù)統(tǒng)計(jì)結(jié)果顯示,紫脊百合和尖被百合的遺傳距離最大,均值為0.887 3,說明二者親緣關(guān)系最遠(yuǎn);紫脊百合為宜昌百合的變種,紫脊百合和宜昌百合遺傳距離最小,均值為0.307 5,說明二者親緣關(guān)系最近。蘭州百合為川百合變種,二者遺傳距離為0.380 1,百合為野百合變種,二者遺傳距離為0.47,遺傳距離均較小,種間遺傳距離多大于0.47,可見變種與原種親緣關(guān)系比種間親緣關(guān)系近。

    通過對(duì)21種野生百合的聚類分析,可以將供試材料分成5類。其中卷瓣組全部聚為一類,百合組大多數(shù)種聚為一類,輪葉組自成一類,鐘花組大多數(shù)聚為一類,聚類結(jié)果大致上與形態(tài)學(xué)分類一致。聚類結(jié)果中第Ⅰ類中,大多數(shù)野生種花被片反卷程度大,與其他類形態(tài)學(xué)上差異明顯;第Ⅱ類中花呈喇叭狀,花被先端向外反卷,程度小于卷瓣組反卷程度;第Ⅲ類與第Ⅱ類花的性狀一致,區(qū)別在于葉子具葉柄;第Ⅳ類為輪葉組的青島百合,其葉片輪生,與其他類群葉片著生方式大不相同;第Ⅴ類野生百合花呈鐘狀,花瓣先端微彎。由此可見,花被片的反卷以及葉子的性狀都是野生百合分類重要的依據(jù)。

    傳統(tǒng)分類中,毛百合被劃分在鐘花組,本研究中毛百合與卷瓣組所有種聚為一類。根據(jù)國(guó)際野生百合形態(tài)學(xué)分組,屬于中國(guó)野生百合鐘花組(sect.Lophophorum)的毛百合被放入一個(gè)中國(guó)分組中沒有的毛百合組(sect.Daurolirion)中,在Michael & Harris(2002)的百合進(jìn)化樹圖中,鐘花組的毛百合與卷瓣組的卷丹、山丹劃在一個(gè)進(jìn)化樹分支上,Lee(2011)通過DNA測(cè)序研究朝鮮百合屬植物進(jìn)化分支,聚類結(jié)果第6分支將毛百合組與大多數(shù)的卷瓣組百合聚為一類,說明毛百合組與卷瓣組百合之間親緣關(guān)系較近。這與 Dubouzet & Shinoda(1999)研究結(jié)果相同。智利等(2011)在對(duì)中國(guó)23個(gè)野生種通過SRAP分子標(biāo)記對(duì)其進(jìn)行遺傳多樣性研究表明,毛百合與卷瓣組聚為一類。由此可見,鐘花組的毛百合與卷瓣組聚為一類,是符合野生百合形態(tài)學(xué)和分子水平分類的,證實(shí)了Nishikawa et al(1999)的發(fā)現(xiàn),國(guó)際分類系統(tǒng)的毛百合組與卷瓣組親緣關(guān)系較近,鐘花組與卷瓣組之間存在基因交流。毛百合雖未能與同樣花冠的種聚在一起,而與花被反卷的卷瓣組聚在一起,說明毛百合系統(tǒng)位置可能比較復(fù)雜,有待隨后繼續(xù)研究。

    在中國(guó)野生百合分組中,野百合和百合屬于百合組,而在Comber(1949)的野生百合分組中,野百合和百合被放在具葉柄組(sect.Archelirion)中,中國(guó)分組中無(wú)具葉柄組。本研究聚類結(jié)果顯示,野百合與百合單獨(dú)聚為一類,因此支持了國(guó)際野生百合形態(tài)學(xué)上的分組,并且說明是否具葉柄對(duì)野生百合的分類是一個(gè)重要的形態(tài)學(xué)性狀。Nishikawa et al(1999)研究證實(shí)很難依據(jù)形態(tài)特征對(duì)百合屬進(jìn)行分類,因?yàn)橛行┓N雖然親緣關(guān)系較遠(yuǎn),但其形態(tài)特征尤其是花型卻非常相似,本研究種毛百合、野百合、百合的聚類形式與傳統(tǒng)形態(tài)分類結(jié)果不盡相同,一定程度上印證了此觀點(diǎn) 。

    本研究證實(shí)了ISSR分子標(biāo)記技術(shù)能較好地從分子水平揭示野生百合種間的親緣關(guān)系,肯定了形態(tài)學(xué)分類的科學(xué)性,并發(fā)現(xiàn)毛百合與卷瓣組百合的親緣關(guān)系相近,鐘花組與卷瓣組野生百合存在基因交流。我國(guó)野生百合種類豐富、生態(tài)類型多樣,曾為百合育種做出了突出貢獻(xiàn)。如20世紀(jì)初,利用野生百合(岷江百合)的抗逆性,瀕臨滅絕的歐洲百合得以拯救?,F(xiàn)今栽培百合雜種系中的亞洲百合雜種系、麝香百合雜種系、喇叭百合雜種系和白花百合雜種系等均有我國(guó)原生種參與雜交而成。因此,進(jìn)一步加強(qiáng)我國(guó)百合遺傳多樣性與親緣關(guān)系研究,對(duì)于野生資源科學(xué)保護(hù)和有效利用具有重要意義 。

    CHEN XQ,XU JM,LIANG SY, et al, 1980. Flora Reipublicae Popularis Sinicae [M]. Beijing: Science Press: 116. [陳心啟,許介眉,梁松筠, 等, 1980. 中國(guó)植物志 [M]. 北京: 科學(xué)出版社: 116.]

    COMBER HF, 1949. A new classification of the genusLilium[M]. Lily Yearbook, Royale Horticulture Society(London), 15: 86-105.

    DEVARUMATH RM, NANDY S, RANI V, et al, 2002. RAPD, ISSR and RFLP fingerprintings as useful markers to evaluate genetic integrity of micropropagated plants of three diploid and triploid elite tea clones representingCamelliasinensis(China type)andC.assamicassp.Assamica(Assam-India type) [J]. Plant Cell Rep, 21: 166-173.

    DUBOUZET JG, SHINODA K, 1999. Phylogenetic analysis of the internal transcribed spacer region of JapaneseLiliumspecies [J]. Theor Appl Genet, 98: 954-960.

    GU X, ZHANG YL, NIU LX, et al, 2013. Pollen morphology observation of 15 wild lilies from four provinces in Western China [J]. Acta Hortic Sin, 40(7): 1 389-1 398. [顧欣, 張延龍, 牛立新, 等, 2013. 中國(guó)西部四省15種野生百合花粉形態(tài)研究 [J]. 園藝學(xué)報(bào), 40(7): 1 389-1 398.]

    GUO WH, JEONG JH, KIM ZS, et al, 2011. Genetic diversity ofLiliumtsingtauensein China and Korea revealed by ISSR markers and morphological characters [J]. Biochem Syst Ecol, 39(4-6): 352-360.

    LEE CS, KIM SC, YEAU SH,et al, 2011. Major lineages of the genusLilium(Liliaceae)based on nrDNA sequences, with special emphasis on the Korean species [J]. J Plant Biol, 54: 159-171.

    MICHAEL JB, HARRIS H, 2002. The gardener’s guide to growing lilies [M]. Oregon: Timber Press: 1-160.

    NISHIKAWA T, OKAZAKI K, UCHINO T, et al, 1999. A molecular phylogeny ofLiliumin the internal transcribed spacer region in nuclear ribosomal DNA [J]. J Mol Evol, 49: 238-249.

    PRADEEP RM, SARLA N, SIDDIQ EA, et al, 2002. Inter simple sequence repeat(ISSR) polymorphism and its application in plant breeding [J]. Euphytica, 128: 9-17.

    QIAN W, GE S, HONG DY,et al, 2001. Genetic variation within and among populations of a wild riceOryzagranulatafrom China detected by RAPD and ISSR markers [J]. Theor Appl Genet, 102: 440-449.

    TSUMURA Y, OHBA K, STRAUS SH, et al, 1996. Diversity and inheritance of inter-simple sequence repeat polymorphisms in Douglas-fir (Pseudotsuga menziesii) and sugi (Cryptom eria japonica) [J]. Theor Appl Genet, 92: 40-45.

    XIANG DY, ZHANG YL,HAO RJ, et al, 2006. Study on the characters description ofLiliumfrom Qinba Mountain and its adjacent area [J]. Chin Agric Sci Bull,22(10): 97-100. [向地英, 張延龍,郝瑞杰, 等, 2006. 秦巴山區(qū)及毗鄰地區(qū)野生百合性狀描述 [J]. 中國(guó)農(nóng)學(xué)通報(bào),22(10): 97-100.]

    YAMAGISHI M, ABE H , NAKANO M, et al, 2002. PCR-based molecular markers in Asiatic hybrid lily [J]. Sci Hortic-Amsterdam, 96: 225-234.

    ZHANG YP, WANG JH, WU LF, et al, 2014. The karyotype of wild and horticultural varieties of lily [J]. J Fujian Agric For Univ, 43(4): 370-373. [張藝萍, 王繼華, 吳麗芳,等, 2014. 百合屬野生種和園藝品種的染色體核型 [J]. 福建農(nóng)林大學(xué)學(xué)報(bào)·自然科學(xué)版, 43(4): 370-373.]

    ZHAO XY,WANG SD,CHEN XL,et al, 1999.Lilium[M]. Beijing:China Agriculture Press:1. [趙祥云,王樹棟,陳新露,等, 1999. 百合 [M]. 北京:中國(guó)農(nóng)業(yè)出版社:1.]

    ZHI L, TENG ZH, LI XY, et al, 2011. Phylogenetic relationship analysis of 23 wild species ofLiliumby SRAP markers [J]. J Agric Biotechnol, 19(4): 677-684. [智利, 滕中華, 李先源, 等, 2011. 23種野生百合遺傳關(guān)系的SRAP分析 [J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 19(4): 677-684.]ZIEKIEWICZ E, RAFALSKI A, LABUDA D, et al,1994. Genome fingerprinting by simple sequence repeat(SSR)-anchored polymerase chain reaction amplification [J]. Genomics, 20: 176-183.

    Genetic relationship analysis of 21 wild species ofLiliumby ISSR markers

    HOU Jun, LUO Jian-Rang, XIAO Fei-Fei, ZHANG Yan-long*, NIU Li-Xin

    (CollegeofLandscapeArchitectureandArts,NorthwestA&FUniversity, Yangling 712100, China )

    To probe the infrageneric phylogenetic relationship ofLilium, seven inter-simple sequence repeat (ISSR) primers were chosen to analyze genetic diversity of 21 wild species ofLilium. POPGENE1.32 and MEGA5.1 were used for data analysis. Results showed that a total of 149 clear DNA bands were amplified, 136 of which were polymorphic, the proportion was 91.06%. The average value of effective number of alleles was 1.762 4. The average value of Nei’s gene diversity was 0.421 4, and the average value of Shannon's Information index was 0.608 5. The genetic distance ranged from 0.307 5 to 0.887 3. The genetic distance betweenL.leucanthumandL.leucanthumvar.centifoliumwas a maximum of selected 21 parts of wild lily material, average value 0.887 3, which indicated that the genetic relationship between them was the fartherest. Whereas the genetic distance betweenL.leucanthumvar.centifoliumandL.lophophorumwas a minimum of selected 21 parts of wild lily material, average value 0.307 5, which indicated genetic relationship between them was the closest. The clustering results were basically consistent with the morphology classification. Twenty-one wildLiliumspecies were classified into five clustering groups by neighbor joining method.L.taliense,L.duchartrei,L.lancifolium,L.nepalense,L.papilliferum,L.davidii,L.davidiivar.unicolor,L.pumilum,L.fargesiiandL.dauricumwere clustered in the first group. AndL.leucanthumBaker var.centifolium,L.leucanthum,L.nepalense,L.regale,L.sargentiaebelong to the second group.L.brownie,L.browniivar.viridulumwere clustered in the third group.L.tsingtauenseformed a distinct group andL.lophophorum,L.amoenum,L.concolorvar.pulchellumwere clustered in the fifth group.L.dauricumhas close genetic relationship with sect. Sinomartagon, and there were gene exchanges existed between Sect. Lophophorum and Sect. Sinomartagon. This indicated that it was not accurate to classify the wildLiliumall depend on the morphology.L.browniivar. viridulum andL.browniewere clustered into one individual group in cluster analysis results. This result showed that whether it has petioles was an important morphological character in the classification of wild lilies. And all the results indicated that ISSR markers were suited to study the phylogenetic relationship of wild lilies.

    Lilium, wild lilies, ISSR, genetic relationship, tree plot

    10.11931/guihaia.gxzw201504028

    2015-04-19

    2015-05-19

    國(guó)家林業(yè)局重點(diǎn)科研項(xiàng)目(2006-73);陜西省林業(yè)廳項(xiàng)目(陜林計(jì)字 [2011]70號(hào));國(guó)家自然科學(xué)基金(305021110);西北農(nóng)林科技大學(xué)基本科研業(yè)務(wù)費(fèi)專項(xiàng)(Z109021002) [Supported by Key Scientific Research Project of State Forestry Administration (2006-73); Project of Shaanxi Forestry Department ( [2011]70); National Natural Science Foundation of China (305021110); Foundamental Scientific Research Fund of Northwest A & F University (Z109021002)]。

    侯珺(1990-),女,陜西渭南人,碩士,主要從事野生百合遺傳多樣性研究,(E-mail)346144734@qq.com。*

    張延龍,教授,主要從事園林植物資源與育種研究,(E-mail) zhangyanlong@nwsuaf.edu.cn。

    Q949

    A

    1000-3142(2016)09-1032-07

    侯珺, 羅建讓, 肖菲菲, 等. 21種野生百合親緣關(guān)系的ISSR分析 [J]. 廣西植物, 2016, 36(9):1032-1038

    HOU J, LUO JR, XIAO FF,et al. Genetic relationship analysis of 21 wild species ofLiliumby ISSR markers [J]. Guihaia, 2016, 36(9):1032-1038

    猜你喜歡
    野百合親緣百合
    谷子近緣野生種的親緣關(guān)系及其利用研究
    中國(guó)醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    風(fēng)吹百合香
    清心消暑話百合
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    野百合的春天
    賈玲:野百合也會(huì)有春天
    海峽姐妹(2018年4期)2018-05-19 02:12:44
    野百合
    民族音樂(2018年4期)2018-01-24 22:12:47
    小白菜種質(zhì)遺傳多樣性與親緣關(guān)系的SRAP 和SSR分析
    野百合的春天
    国产国语露脸激情在线看| 手机成人av网站| 久久人人爽av亚洲精品天堂| 丰满人妻熟妇乱又伦精品不卡| 级片在线观看| 女性生殖器流出的白浆| 国产精品香港三级国产av潘金莲| 国产又爽黄色视频| 国产精品国产av在线观看| 亚洲精品久久成人aⅴ小说| 欧美午夜高清在线| a级片在线免费高清观看视频| 长腿黑丝高跟| 日韩视频一区二区在线观看| 久久中文字幕一级| 久久久久久久精品吃奶| 在线观看免费日韩欧美大片| 制服诱惑二区| 成人国语在线视频| 色尼玛亚洲综合影院| 精品久久久久久久久久免费视频 | 一区二区三区国产精品乱码| 久久精品亚洲精品国产色婷小说| 亚洲国产精品合色在线| 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三区在线| 国产三级在线视频| 男人舔女人下体高潮全视频| 久久中文看片网| 亚洲熟妇熟女久久| 叶爱在线成人免费视频播放| 男女下面插进去视频免费观看| 嫁个100分男人电影在线观看| tocl精华| 日本撒尿小便嘘嘘汇集6| 可以免费在线观看a视频的电影网站| 夫妻午夜视频| 久久久久久久久免费视频了| 免费一级毛片在线播放高清视频 | 韩国av一区二区三区四区| 免费在线观看视频国产中文字幕亚洲| 国产精品亚洲一级av第二区| 国产精品综合久久久久久久免费 | 美女福利国产在线| 国产伦一二天堂av在线观看| 一二三四在线观看免费中文在| 激情视频va一区二区三区| 久久久久久人人人人人| 亚洲精品一区av在线观看| 午夜福利欧美成人| 在线免费观看的www视频| 欧美激情极品国产一区二区三区| 国产成人精品在线电影| 国产成人精品久久二区二区91| 桃红色精品国产亚洲av| 精品一区二区三区视频在线观看免费 | 亚洲一区二区三区不卡视频| 高清欧美精品videossex| 国产精品1区2区在线观看.| 男人舔女人的私密视频| 国产精品久久久av美女十八| 国产91精品成人一区二区三区| 19禁男女啪啪无遮挡网站| 国内毛片毛片毛片毛片毛片| 母亲3免费完整高清在线观看| 波多野结衣一区麻豆| av片东京热男人的天堂| 亚洲男人的天堂狠狠| 精品国产乱子伦一区二区三区| 咕卡用的链子| 欧美精品亚洲一区二区| 身体一侧抽搐| 巨乳人妻的诱惑在线观看| 国产av一区二区精品久久| 国产一区二区三区综合在线观看| 精品久久久久久久毛片微露脸| 欧美在线一区亚洲| 亚洲黑人精品在线| 麻豆国产av国片精品| 国产成人欧美在线观看| 国产精品一区二区免费欧美| 一本大道久久a久久精品| 国产精品乱码一区二三区的特点 | 午夜免费成人在线视频| 91麻豆精品激情在线观看国产 | 可以在线观看毛片的网站| 国产精品久久久久成人av| 亚洲人成网站在线播放欧美日韩| 久久久久久免费高清国产稀缺| 久久人人精品亚洲av| 亚洲,欧美精品.| 99国产综合亚洲精品| 日韩大码丰满熟妇| 精品久久久久久电影网| 国产日韩一区二区三区精品不卡| 国产亚洲精品综合一区在线观看 | 黄色 视频免费看| 亚洲国产精品sss在线观看 | 精品熟女少妇八av免费久了| 18禁黄网站禁片午夜丰满| 婷婷丁香在线五月| 免费av中文字幕在线| 精品高清国产在线一区| 手机成人av网站| 欧美日韩黄片免| 99久久人妻综合| 男女午夜视频在线观看| 婷婷六月久久综合丁香| 久久精品影院6| av网站免费在线观看视频| 日本一区二区免费在线视频| 久久久久久亚洲精品国产蜜桃av| 成人18禁高潮啪啪吃奶动态图| www日本在线高清视频| 亚洲精品一卡2卡三卡4卡5卡| 国产区一区二久久| 国产97色在线日韩免费| √禁漫天堂资源中文www| 丝袜美足系列| 国产精品免费视频内射| 日韩免费高清中文字幕av| 午夜精品在线福利| 美女 人体艺术 gogo| 成熟少妇高潮喷水视频| 91老司机精品| 欧美性长视频在线观看| 亚洲一区二区三区色噜噜 | 亚洲激情在线av| 男女床上黄色一级片免费看| 欧美激情高清一区二区三区| 午夜福利在线观看吧| 黄色片一级片一级黄色片| 成人18禁在线播放| 黄色毛片三级朝国网站| 老司机在亚洲福利影院| 12—13女人毛片做爰片一| 一级a爱片免费观看的视频| 1024香蕉在线观看| 亚洲av日韩精品久久久久久密| 亚洲国产欧美一区二区综合| 久久性视频一级片| 十分钟在线观看高清视频www| 伦理电影免费视频| 欧美日韩国产mv在线观看视频| 麻豆成人av在线观看| 身体一侧抽搐| 丁香欧美五月| 黄色视频不卡| 电影成人av| 人人妻人人澡人人看| 一级片'在线观看视频| 可以在线观看毛片的网站| 免费在线观看视频国产中文字幕亚洲| 成人特级黄色片久久久久久久| 最好的美女福利视频网| 黄片小视频在线播放| 国产亚洲精品久久久久5区| 天天影视国产精品| 日韩欧美一区视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 日日摸夜夜添夜夜添小说| 亚洲自偷自拍图片 自拍| 久久婷婷成人综合色麻豆| 美女午夜性视频免费| 天堂√8在线中文| 男女之事视频高清在线观看| 青草久久国产| www.精华液| 精品国内亚洲2022精品成人| 水蜜桃什么品种好| 免费看a级黄色片| 一二三四社区在线视频社区8| 男女下面插进去视频免费观看| 女警被强在线播放| 高清在线国产一区| netflix在线观看网站| 欧美日韩精品网址| 神马国产精品三级电影在线观看 | 香蕉丝袜av| 高清毛片免费观看视频网站 | 一本综合久久免费| 老汉色∧v一级毛片| 男女下面进入的视频免费午夜 | 久久午夜综合久久蜜桃| 一本大道久久a久久精品| 每晚都被弄得嗷嗷叫到高潮| 久久中文看片网| 美女高潮到喷水免费观看| 国产一区二区三区视频了| av中文乱码字幕在线| 99热只有精品国产| 国产精品久久久久久人妻精品电影| 久久久久亚洲av毛片大全| 日韩免费av在线播放| 黄色 视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 99在线人妻在线中文字幕| 香蕉丝袜av| 乱人伦中国视频| 另类亚洲欧美激情| а√天堂www在线а√下载| 国产高清国产精品国产三级| 麻豆久久精品国产亚洲av | 国产精品永久免费网站| 在线观看66精品国产| 黄色女人牲交| 欧美日韩视频精品一区| 免费观看精品视频网站| 九色亚洲精品在线播放| 亚洲国产看品久久| 免费人成视频x8x8入口观看| 成人三级做爰电影| 国产亚洲精品综合一区在线观看 | 久久中文字幕人妻熟女| 夜夜爽天天搞| 日本vs欧美在线观看视频| 欧美最黄视频在线播放免费 | 男人舔女人下体高潮全视频| 亚洲性夜色夜夜综合| 老司机午夜十八禁免费视频| 精品第一国产精品| 老司机在亚洲福利影院| 精品久久久久久,| av有码第一页| 成年人黄色毛片网站| 久久久久精品国产欧美久久久| 亚洲一区中文字幕在线| 日韩精品青青久久久久久| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区| 级片在线观看| 亚洲熟女毛片儿| 一进一出好大好爽视频| 国产成人精品无人区| 人妻久久中文字幕网| 国产精品久久久久成人av| 成人国产一区最新在线观看| 久久精品影院6| 久久午夜亚洲精品久久| 人人澡人人妻人| 亚洲精品国产区一区二| 岛国视频午夜一区免费看| 丝袜人妻中文字幕| 国产精品永久免费网站| 国产乱人伦免费视频| 美女午夜性视频免费| 亚洲av熟女| 色综合欧美亚洲国产小说| 一二三四在线观看免费中文在| 村上凉子中文字幕在线| 国产黄a三级三级三级人| 国产av又大| 午夜免费鲁丝| 国产精品亚洲一级av第二区| 他把我摸到了高潮在线观看| 在线观看午夜福利视频| 亚洲第一欧美日韩一区二区三区| 欧美精品啪啪一区二区三区| 免费看a级黄色片| 丝袜美腿诱惑在线| 99久久国产精品久久久| 中文字幕高清在线视频| 精品电影一区二区在线| 国产黄a三级三级三级人| 国产色视频综合| 中出人妻视频一区二区| 国产精品久久久久成人av| 村上凉子中文字幕在线| 日本a在线网址| 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 久久精品国产清高在天天线| 国产精华一区二区三区| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三区在线| 国产三级黄色录像| 两人在一起打扑克的视频| 91九色精品人成在线观看| 成人三级黄色视频| 日本a在线网址| 国产精品久久电影中文字幕| 国产精品99久久99久久久不卡| 久久天堂一区二区三区四区| 免费看十八禁软件| 亚洲在线自拍视频| ponron亚洲| 免费高清在线观看日韩| 国产区一区二久久| 欧美激情久久久久久爽电影 | 国产精品香港三级国产av潘金莲| 一夜夜www| 精品久久久久久成人av| 一级黄色大片毛片| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 中文字幕高清在线视频| 涩涩av久久男人的天堂| 在线看a的网站| 夫妻午夜视频| 夜夜看夜夜爽夜夜摸 | 欧美日韩福利视频一区二区| 国产成人精品久久二区二区免费| 国产成人系列免费观看| 麻豆av在线久日| 久久国产精品人妻蜜桃| 三上悠亚av全集在线观看| 视频区欧美日本亚洲| 性欧美人与动物交配| 美女国产高潮福利片在线看| 午夜福利,免费看| 在线观看日韩欧美| 搡老熟女国产l中国老女人| avwww免费| 真人一进一出gif抽搐免费| 亚洲国产欧美网| 男女下面进入的视频免费午夜 | 手机成人av网站| 亚洲精品一二三| 欧美日韩福利视频一区二区| 中文字幕高清在线视频| 男女下面进入的视频免费午夜 | 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲| 这个男人来自地球电影免费观看| 免费在线观看日本一区| 久久人人97超碰香蕉20202| 久久中文字幕一级| www.熟女人妻精品国产| 身体一侧抽搐| 亚洲九九香蕉| www.熟女人妻精品国产| 久久婷婷成人综合色麻豆| 日本免费a在线| 日本三级黄在线观看| 亚洲精品国产区一区二| 亚洲人成77777在线视频| 美女 人体艺术 gogo| 亚洲一区中文字幕在线| 久久精品国产亚洲av香蕉五月| 五月开心婷婷网| 很黄的视频免费| 手机成人av网站| 91精品三级在线观看| 国产精品久久视频播放| 女人被狂操c到高潮| 热99国产精品久久久久久7| 在线观看66精品国产| 中文字幕人妻丝袜一区二区| 热99re8久久精品国产| 色老头精品视频在线观看| 99精品欧美一区二区三区四区| 99热国产这里只有精品6| 国产成+人综合+亚洲专区| 国产精品久久久人人做人人爽| 如日韩欧美国产精品一区二区三区| 亚洲av片天天在线观看| 国产精品香港三级国产av潘金莲| 香蕉国产在线看| 91九色精品人成在线观看| 国产精品香港三级国产av潘金莲| 香蕉国产在线看| 欧美乱色亚洲激情| 国产精品国产av在线观看| 91国产中文字幕| 亚洲欧美一区二区三区久久| 久久亚洲精品不卡| 后天国语完整版免费观看| 欧美乱色亚洲激情| 国产91精品成人一区二区三区| 黄色片一级片一级黄色片| 亚洲一区中文字幕在线| 在线免费观看的www视频| 18禁黄网站禁片午夜丰满| 国产精品国产高清国产av| 午夜福利,免费看| 日本vs欧美在线观看视频| 欧美精品亚洲一区二区| 五月开心婷婷网| 免费久久久久久久精品成人欧美视频| 精品国产亚洲在线| 欧美日韩亚洲综合一区二区三区_| 男女高潮啪啪啪动态图| 欧美+亚洲+日韩+国产| 18禁美女被吸乳视频| 可以免费在线观看a视频的电影网站| 女同久久另类99精品国产91| 亚洲人成77777在线视频| 法律面前人人平等表现在哪些方面| 丝袜在线中文字幕| 免费观看人在逋| 91老司机精品| 国产av在哪里看| 亚洲五月天丁香| 久久久久久久久中文| 国产蜜桃级精品一区二区三区| 久久精品亚洲av国产电影网| 看免费av毛片| 一级a爱片免费观看的视频| 久久精品亚洲av国产电影网| 精品久久久精品久久久| 黄色女人牲交| 国产精品秋霞免费鲁丝片| 制服人妻中文乱码| 国产成人欧美| 午夜福利免费观看在线| 亚洲男人的天堂狠狠| 欧美+亚洲+日韩+国产| 亚洲一区二区三区欧美精品| 久久 成人 亚洲| 亚洲国产看品久久| 日本wwww免费看| 亚洲欧美日韩另类电影网站| 黄色视频,在线免费观看| 成人影院久久| 欧美日本亚洲视频在线播放| 十八禁人妻一区二区| 搡老熟女国产l中国老女人| 宅男免费午夜| 欧美黄色淫秽网站| 欧美日韩黄片免| 黄频高清免费视频| 精品国产乱子伦一区二区三区| 美女大奶头视频| 久9热在线精品视频| 国产乱人伦免费视频| 女生性感内裤真人,穿戴方法视频| 精品午夜福利视频在线观看一区| 超碰97精品在线观看| 久久人妻熟女aⅴ| 国产成人av激情在线播放| 在线观看一区二区三区激情| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 国产精品日韩av在线免费观看 | 制服诱惑二区| 无限看片的www在线观看| 91精品三级在线观看| 久久久国产欧美日韩av| 一二三四社区在线视频社区8| 人妻久久中文字幕网| 久久人人97超碰香蕉20202| 色精品久久人妻99蜜桃| 亚洲欧美一区二区三区久久| 黄片小视频在线播放| 可以免费在线观看a视频的电影网站| 激情视频va一区二区三区| 色综合站精品国产| 亚洲欧洲精品一区二区精品久久久| 新久久久久国产一级毛片| 成人手机av| 美女国产高潮福利片在线看| 国产一区二区在线av高清观看| 男女高潮啪啪啪动态图| 村上凉子中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 在线十欧美十亚洲十日本专区| 国产欧美日韩精品亚洲av| √禁漫天堂资源中文www| 神马国产精品三级电影在线观看 | 超色免费av| 欧美国产精品va在线观看不卡| 久久精品影院6| 免费av中文字幕在线| www.自偷自拍.com| 99riav亚洲国产免费| 不卡av一区二区三区| 久久久国产精品麻豆| 亚洲在线自拍视频| 亚洲片人在线观看| 亚洲人成伊人成综合网2020| 在线观看www视频免费| 日本免费a在线| 久久 成人 亚洲| 欧美激情 高清一区二区三区| 少妇的丰满在线观看| 97碰自拍视频| 香蕉国产在线看| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 色尼玛亚洲综合影院| 在线观看一区二区三区激情| 国产成人av教育| 午夜两性在线视频| 亚洲黑人精品在线| 国产伦一二天堂av在线观看| a级片在线免费高清观看视频| 亚洲人成电影免费在线| 怎么达到女性高潮| 超碰97精品在线观看| 88av欧美| 午夜精品在线福利| 色尼玛亚洲综合影院| 欧美日韩国产mv在线观看视频| 亚洲专区中文字幕在线| 天堂动漫精品| 中文字幕人妻熟女乱码| 9热在线视频观看99| 国产亚洲精品久久久久5区| 热re99久久国产66热| av视频免费观看在线观看| 亚洲欧美精品综合久久99| 男人舔女人的私密视频| 国产精华一区二区三区| 国产乱人伦免费视频| av网站免费在线观看视频| www日本在线高清视频| 99国产综合亚洲精品| 午夜91福利影院| 五月开心婷婷网| 亚洲精品av麻豆狂野| 黄片小视频在线播放| 久久久久精品国产欧美久久久| 亚洲第一av免费看| 精品一区二区三区四区五区乱码| 女人被躁到高潮嗷嗷叫费观| 日韩大码丰满熟妇| 9热在线视频观看99| 最近最新中文字幕大全免费视频| 韩国av一区二区三区四区| 黄色视频不卡| 99精品在免费线老司机午夜| 99精品欧美一区二区三区四区| 久久久久精品国产欧美久久久| av在线播放免费不卡| 精品欧美一区二区三区在线| 搡老岳熟女国产| 久久久国产一区二区| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 91成人精品电影| av片东京热男人的天堂| 乱人伦中国视频| 97超级碰碰碰精品色视频在线观看| 黄色 视频免费看| 亚洲一码二码三码区别大吗| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 怎么达到女性高潮| 黄色视频不卡| 亚洲精品av麻豆狂野| 19禁男女啪啪无遮挡网站| 亚洲色图综合在线观看| 国产欧美日韩精品亚洲av| 一级毛片女人18水好多| www.999成人在线观看| 成人影院久久| 国产一卡二卡三卡精品| 熟女少妇亚洲综合色aaa.| 搡老岳熟女国产| 午夜久久久在线观看| 日本 av在线| 久久久久久亚洲精品国产蜜桃av| 国产欧美日韩一区二区三区在线| 别揉我奶头~嗯~啊~动态视频| 亚洲成人久久性| 欧美色视频一区免费| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 不卡一级毛片| 后天国语完整版免费观看| 丝袜在线中文字幕| 欧美日韩一级在线毛片| 久久欧美精品欧美久久欧美| 99久久国产精品久久久| 日韩精品青青久久久久久| 黄色片一级片一级黄色片| 日本五十路高清| 亚洲成国产人片在线观看| 巨乳人妻的诱惑在线观看| 久久久国产一区二区| 少妇被粗大的猛进出69影院| 欧美大码av| 88av欧美| 亚洲精品av麻豆狂野| 亚洲中文av在线| 99国产精品99久久久久| 国产一区二区三区在线臀色熟女 | 国产91精品成人一区二区三区| 精品第一国产精品| 99国产精品免费福利视频| 黄色毛片三级朝国网站| 女性生殖器流出的白浆| 精品久久久精品久久久| 日本免费a在线| 女性被躁到高潮视频| 久久精品91无色码中文字幕| 精品国产超薄肉色丝袜足j| 亚洲人成电影免费在线| 99国产极品粉嫩在线观看| 老司机午夜福利在线观看视频| 精品国产一区二区三区四区第35| www.精华液| 在线免费观看的www视频| av国产精品久久久久影院| 亚洲狠狠婷婷综合久久图片| 国产熟女xx| 精品福利永久在线观看| 黑人巨大精品欧美一区二区mp4| 国产色视频综合| 日韩精品青青久久久久久| 少妇的丰满在线观看| 久久久久久人人人人人| 久久久精品国产亚洲av高清涩受| 少妇粗大呻吟视频| 人人妻人人澡人人看| 天天躁狠狠躁夜夜躁狠狠躁| 夫妻午夜视频| 岛国视频午夜一区免费看| 可以在线观看毛片的网站| 国内毛片毛片毛片毛片毛片| a级毛片黄视频| 久久精品亚洲av国产电影网| 久久人人97超碰香蕉20202| 亚洲第一欧美日韩一区二区三区|