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

    基于psbA-trnH分析的何首烏野生居群遺傳多樣性

    2012-09-06 01:51:22白明明孫小芹郭建林李密密杭悅宇
    關(guān)鍵詞:居群何首烏堿基

    白明明,孫小芹,郭建林,李密密,杭悅宇

    〔江蘇省·中國(guó)科學(xué)院植物研究所(南京中山植物園)江蘇省植物遷地保護(hù)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210014〕

    基于psbA-trnH分析的何首烏野生居群遺傳多樣性

    白明明,孫小芹,郭建林,李密密,杭悅宇①

    〔江蘇省·中國(guó)科學(xué)院植物研究所(南京中山植物園)江蘇省植物遷地保護(hù)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210014〕

    對(duì)產(chǎn)自不同省區(qū)的何首烏〔Fallopia multiflora(Thunb.)Harald.〕17個(gè)野生居群85個(gè)單株的psbA-trnH序列進(jìn)行了擴(kuò)增和分析,在此基礎(chǔ)上分析了居群間的遺傳多樣性,并采用NJ法對(duì)85個(gè)單株進(jìn)行了聚類分析。結(jié)果表明:供試85個(gè)單株的psbA-trnH序列長(zhǎng)度為384 bp;其中,變異位點(diǎn)為167 bp,簡(jiǎn)約信息位點(diǎn)為53 bp,分別占序列總長(zhǎng)度的43.5%和13.8%。變異類型主要為堿基缺失和替換;變異位點(diǎn)主要集中在235~281 bp區(qū)域,根據(jù)位點(diǎn)變異情況可將17個(gè)居群大體分為3類。各居群間的遺傳距離為0.000~0.172,其中,貴州居群與其他16個(gè)居群間的遺傳距離為0.167~0.172,而其他16個(gè)居群間的遺傳距離為0.000~0.017。17個(gè)居群間的核苷酸多樣性指數(shù)(Pi)、遺傳分化系數(shù)(Nst)和基因流(Nm)分別為0.028 56、0.918 68和0.04;除貴州居群外其他16個(gè)居群的Pi、Nst和Nm分別為0.015 68、0.837 19和0.10;貴州居群與其周邊省區(qū)(四川、云南、廣西、湖南和湖北)居群的Pi、Nst和Nm分別為0.047 99、0.937 62和0.03。在NJ系統(tǒng)樹(shù)上,17個(gè)居群可聚為4支,且大部分居群的供試單株聚在同一分支中;僅貴州居群?jiǎn)为?dú)聚為一支,與序列分析的劃分結(jié)果基本一致。由研究結(jié)果可見(jiàn):野生何首烏居群總遺傳變異的91.868%來(lái)自居群間,8.132%來(lái)自居群內(nèi),居群間的基因交流較少;除貴州居群外其余16個(gè)居群的整體遺傳多樣性水平偏低,說(shuō)明何首烏居群整體多樣性水平在很大程度上受貴州居群的影響。

    何首烏;居群;psbA-trnH序列;遺傳多樣性;NJ系統(tǒng)樹(shù)

    何首烏〔Fallopia multiflora(Thunb.)Harald.〕為蓼科(Polygonaceae)何首烏屬(FallopiaAdans.)多年生草本植物[1],其塊根和藤莖均為傳統(tǒng)中藥材,具有解毒、消癰、截瘧、潤(rùn)腸通便的功效[2]。何首烏野生資源豐富,主要分布在安徽、福建、廣東、甘肅、廣西、貴州、湖北、河南、湖南、江蘇、江西、四川、山東、上海、陜西、云南和浙江。

    目前,有關(guān)何首烏遺傳多樣性的研究非常少。王凌暉等[3]采用RAPD技術(shù)對(duì)廣西產(chǎn)10個(gè)野生何首烏居群進(jìn)行了遺傳多樣性分析,用40對(duì)隨機(jī)引物共擴(kuò)增出175個(gè)多態(tài)性片段,多態(tài)性百分率達(dá)59.3%;嚴(yán)寒靜等[4]采用PCR直接測(cè)序法測(cè)定何首烏10個(gè)種源核糖體DNA的ITS序列,共找到22個(gè)變異位點(diǎn),其中ITS1區(qū)域有4個(gè)、ITS2區(qū)域有16個(gè)、5.8S區(qū)域有2個(gè),序列間遺傳分化距離0.001 75~0.049 45。對(duì)于野生何首烏來(lái)說(shuō),上述研究涉及的樣品分布范圍并不廣、居群數(shù)也偏少,且所用片段的變異位點(diǎn)不能充分反映何首烏野生種質(zhì)資源的遺傳結(jié)構(gòu)。

    psbA-trnH基因間隔區(qū)位于編碼光合系統(tǒng)Ⅱ反應(yīng)中心DⅠ蛋白的psbA基因和編碼tRNA組氨酸的trnH基因之間,被認(rèn)為是葉綠體基因組中進(jìn)化速率最快的基因間隔區(qū)之一。目前,已將psbA-trnH序列分析用于薯蕷屬(DioscoreaL.)3個(gè)種[5]、石斛屬(DendrobiumSw.)15個(gè)種[6]和連翹〔Forsythia suspensa(Thunb.)Vahl〕3個(gè)居群[7]的遺傳多樣性研究。

    作者對(duì)采自17個(gè)省區(qū)的野生何首烏85個(gè)單株的psbA-trnH序列進(jìn)行了分析,比較分析野生何首烏的遺傳多樣性,探討不同居群的遺傳變異程度,以期為其優(yōu)良品種的選育和種質(zhì)資源的研究提供依據(jù)。

    1 材料和方法

    1.1 材料

    供試的17個(gè)野生何首烏居群的分布地及基本地理信息見(jiàn)表1。憑證標(biāo)本存放于江蘇省·中國(guó)科學(xué)院植物研究所標(biāo)本館(NAS),經(jīng)江蘇省·中國(guó)科學(xué)院植物研究所杭悅宇研究員鑒定。于2010年在每個(gè)居群選取5株樣株,采集幼嫩且無(wú)病蟲(chóng)害的干凈葉片,用變色硅膠干燥后保存、待用。

    實(shí)驗(yàn)用10×PCR buffer(含Mg2+)、TaqDNA聚合酶和dNTPs均購(gòu)自上海博彩生物科技有限公司,引物由上海英駿生物技術(shù)有限公司合成,AxyPrep DNA凝膠回收試劑盒購(gòu)自杭州愛(ài)思進(jìn)生物技術(shù)有限公司。

    1.2 方法

    1.2.1 基因組總DNA的提取參照文獻(xiàn)[8]采用改良CTAB法提取基因組總DNA。取10 mg葉片,在液氮中研磨后加入500 μL提取液于65℃水浴中保溫45 min,期間搖動(dòng)30次;加入500 μL體積比24∶1的氯仿-異戊醇混合液,混合均勻,12 000 r·min-1離心10 min;取上清液,加入2倍體積無(wú)水乙醇,-20℃放置1 h以上;4 000 r·min-1離心10~20 min使DNA沉淀;用700 μL體積分?jǐn)?shù)70%乙醇清洗1次,4 000 r ·min-1離心5~10 min;再用700 μL無(wú)水乙醇清洗1次,4 000 r·min-1離心5~10 min;得到DNA沉淀,晾干,最后加入100 μL水或TE溶解DNA。

    1.2.2psbA-trnH序列的擴(kuò)增及產(chǎn)物檢測(cè)根據(jù)GenBank中發(fā)布的何首烏psbA-trnH序列(GenBank登錄號(hào)為EU554047.1),采用Primer Premier 5.0軟件自行設(shè)計(jì)psbA-trnH序列引物,F(xiàn)HS和FHA引物序列分別為5'-TTCCCGCTAGACCTAGCTGC-3'和5'-ACTGCCTTGATCCACTTGGC-3'。

    擴(kuò)增反應(yīng)在PE-9600型PCR儀(Perkin Elmer公司生產(chǎn))上進(jìn)行。擴(kuò)增反應(yīng)體系總體積20 μL,含30 ng模板DNA、10×PCR buffer(含Mg2+)2 μL、10 mmol·L-1dNTPs 0.3 μL、5 μmol·L-1引物2 μL和5 U·μL-1TaqDNA聚合酶0.2 μL,用滅菌雙蒸水補(bǔ)足至20 μL。擴(kuò)增程序?yàn)?94℃預(yù)變性5 min; 94℃變性40 s,58℃退火40 s,70℃延伸45 s,共30個(gè)循環(huán);最后于72℃延伸5 min。

    表1 供試何首烏野生居群的基本概況1)Table 1 Basic status of wild populations of Fallopia multiflora(Thunb.)Harald.tested1)

    擴(kuò)增產(chǎn)物與6×甘油凝膠上樣液混合后,用質(zhì)量體積分?jǐn)?shù)0.8%瓊脂糖凝膠(含0.5 μg·mL-11×EB)進(jìn)行電泳檢測(cè),電泳時(shí)間約40 min。電泳結(jié)束后用WV-BP330型凝膠掃描分析系統(tǒng)(江蘇捷達(dá)科技發(fā)展有限公司生產(chǎn))觀察擴(kuò)增結(jié)果并拍照。

    1.2.3 擴(kuò)增片段的純化和序列測(cè)定PCR擴(kuò)增產(chǎn)物用AxyPrep DNA凝膠回收試劑盒進(jìn)行純化,純化后使用引物FHS和FHA直接進(jìn)行雙向測(cè)序,測(cè)序由上海華大基因有限公司完成。

    1.3 序列分析

    采用Sequencher軟件進(jìn)行峰圖文件的編輯和拼接,拼接好的序列用MEGA 4.1軟件進(jìn)行多重序列比對(duì)(multiple sequence alignment,MSA),對(duì)序列長(zhǎng)度、保守位點(diǎn)、變異位點(diǎn)和簡(jiǎn)約信息位點(diǎn)等序列特征進(jìn)行分析[9];通過(guò)Kimura 2-parameter(K2P)模型計(jì)算遺傳距離,構(gòu)建距離矩陣[10];應(yīng)用DnaSP 5.10[11]軟件計(jì)算核苷酸多樣性指數(shù)[12]、基因流[13]和居群間遺傳分化系數(shù)[13];以萹蓄(PolygonumaviculareL.) (GenBank登錄號(hào)為FJ395458.1)為外類群,通過(guò)鄰接法(neighbor-joining method,NJ)構(gòu)建基于遺傳距離的聚類圖,構(gòu)建的系統(tǒng)發(fā)育樹(shù)采用1 000次循環(huán)的自舉檢驗(yàn)法(bootstrap test)檢驗(yàn)各分支的置信度,空缺(gap)始終作為缺失(missing)處理[14]。

    2 結(jié)果和分析

    2.1 何首烏居群psbA-trnH序列擴(kuò)增結(jié)果分析

    供試的85個(gè)野生何首烏單株psbA-trnH序列的長(zhǎng)度為243~434 bp,排序后兩端切平獲得的序列長(zhǎng)度為384 bp,當(dāng)空缺(gap)始終作缺失(missing)處理時(shí),保守位點(diǎn)有217 bp,占序列總長(zhǎng)度的56.5%;變異位點(diǎn)有167 bp,占序列總長(zhǎng)度的43.5%,變異類型主要為堿基缺失和替換。其中,簡(jiǎn)約信息位點(diǎn)有53 bp,占序列總長(zhǎng)度的13.8%。

    2.2 何首烏居群psbA-trnH序列的差異分析

    何首烏不同居群psbA-trnH序列變異位點(diǎn)的比較結(jié)果見(jiàn)表2。由表2可見(jiàn):何首烏不同居群間和單株間psbA-trnH序列均存在堿基變異位點(diǎn),主要集中在235~281 bp區(qū)域。福建、廣東、湖北、湖南、江西、四川和浙江居群的所有單株,山東居群的4號(hào)單株,上海居群的2號(hào)和4號(hào)單株,廣西居群的1~4號(hào)單株以及云南居群的2~5號(hào)單株的psbA-trnH序列在235~281 bp位點(diǎn)的堿基序列依次為TATGGCCTCT TTGGTGT(或G)TTTG(或T或A)GTTGGGTGGACT TTTTTTTGATTCAT;貴州居群5個(gè)單株的psbA-trnH序列在235~281 bp位點(diǎn)的堿基序列依次為TCTCCCTCTTT--TGTTTT----GGG-GGAAATTTCTTT--TTCTT(“-”表示堿基缺失)。比較而言,貴州居群psbA-trnH序列235~281 bp區(qū)域間有16個(gè)變異位點(diǎn),與其他11個(gè)居群的差異較大。在251 bp位置,其他11個(gè)居群的堿基分別為T和G,其中,福建、湖北、江西、四川和浙江居群的所有單株,山東居群的4號(hào)單株,上海居群的2號(hào)和4號(hào)單株及云南居群的2~5號(hào)單株的堿基為T;廣東和湖南居群的所有單株以及廣西居群的1~4號(hào)單株的堿基為G。在255 bp位置,福建、湖北、江西和浙江居群的所有單株,山東居群的4號(hào)單株以及上海居群的2號(hào)和4號(hào)單株的堿基為G;廣東和湖南居群的所有單株以及廣西居群的1~4號(hào)單株的堿基為T;四川居群的5個(gè)單株和云南居群的2~5號(hào)單株的堿基為A。

    安徽、甘肅、河南、江蘇和陜西居群所有單株的psbA-trnH序列在235~281 bp位點(diǎn)的堿基均缺失;山東居群的1~3號(hào)和5號(hào)單株,上海居群的1號(hào)、3號(hào)和5號(hào)單株,廣西居群的5號(hào)單株及云南居群的1號(hào)單株在這段序列位點(diǎn)也存在堿基缺失。

    根據(jù)psbA-trnH序列235~281 bp區(qū)域的堿基差異,17個(gè)野生何首烏居群可被劃分為3大組:第1組包括福建、廣東、湖北、湖南、江西、四川、浙江、廣西和云南居群;第2組僅貴州居群,其特征為235~281 bp序列與第1組相比有16個(gè)變異位點(diǎn);第3組包括安徽、甘肅、河南、江蘇、陜西、山東和上海居群,其特征為235~281 bp區(qū)域堿基缺失。另外,第1組又可分為3小組:第1小組包括廣東、湖南和廣西居群,其251 bp位點(diǎn)的堿基為G、255 bp位點(diǎn)的堿基為T;第2小組包括福建、湖北、江西和浙江居群,其251和255 bp 2個(gè)位點(diǎn)的堿基恰好與第1小組完全相反;第3小組包括四川和云南居群,其251和255 bp的堿基分別為T和A。

    表2 何首烏17個(gè)野生居群psbA-trnH序列變異位點(diǎn)的比較Table 2 Comparison of variable sites in psbA-trnH sequence of seventeen wild populations of Fallopia multiflora(Thunb.)Harald.

    續(xù)表2Table 2 (Continued)

    續(xù)表2>Table 2 (Continued)

    2.3 何首烏居群的遺傳多樣性分析

    采用MEGA 4.1軟件中的Kimura 2-parameter模型計(jì)算何首烏不同居群間psbA-trnH序列的遺傳距離并構(gòu)建序列距離矩陣,結(jié)果見(jiàn)表3。何首烏各居群間的遺傳距離為0.000~0.172;其中,貴州居群與其他居群間的遺傳距離均較大,為0.167~0.172;其他居群間的遺傳距離均較小,為0.000~0.017。

    何首烏17個(gè)野生居群的核苷酸多樣性指數(shù)為0.028 56±0.005 43,居群間的遺傳分化系數(shù)達(dá)到0.918 68,基因流為0.04。說(shuō)明17個(gè)何首烏野生居群總遺傳變異的91.868%來(lái)自居群間,8.132%來(lái)自居群內(nèi),居群間的基因交流較少。

    由于貴州居群的psbA-trnH序列與其他居群有較大差異,因此,對(duì)其他16個(gè)居群以及貴州居群與其周邊省區(qū)居群的遺傳多樣性進(jìn)行了分析,結(jié)果顯示:其他16個(gè)居群的核苷酸多樣性指數(shù)為0.015 68± 0.000 77,居群間遺傳分化系數(shù)為0.837 19,基因流為0.10;貴州居群與其周邊省區(qū)(四川、云南、廣西、湖南和湖北)居群的核苷酸多樣性指數(shù)為0.047 99± 0.010 88,居群間遺傳分化系數(shù)為0.937 62,基因流為0.03。說(shuō)明其他16個(gè)何首烏野生居群的整體遺傳多樣性水平偏低;貴州居群與其周邊省區(qū)居群的遺傳分化系數(shù)大于其他16個(gè)居群,基因交流則小于其他16個(gè)居群,說(shuō)明何首烏野生居群的遺傳多樣性整體水平在很大程度上受貴州居群的影響。

    表3 何首烏17個(gè)野生居群間的遺傳距離1)Table 3 Genetic distance among seventeen wild populations of Fallopia multiflora(Thunb.)Harald.1)

    2.4 何首烏居群的聚類分析

    以萹蓄為外類群,依據(jù)遺傳距離、采用鄰接法(NJ)對(duì)何首烏17個(gè)野生居群進(jìn)行聚類分析并構(gòu)建遺傳關(guān)系樹(shù),結(jié)果見(jiàn)圖1。由圖1可看出:17個(gè)何首烏野生居群的85個(gè)單株可聚成4支。浙江、江西、湖北、福建和四川居群的所有單株,山東居群的4號(hào)單株,上海居群的2號(hào)和4號(hào)單株及云南居群的2~5號(hào)單株聚在一起成為第1支,自展支持率達(dá)57%,其中,四川居群的5個(gè)單株和云南居群的2~5號(hào)單株又獨(dú)立成支,與其他居群分離,但自展支持率僅為54%。廣東和湖南居群的所有單株及廣西居群的1~4號(hào)單株聚在一起成為第2支,自展支持率達(dá)86%。第3支包括安徽、甘肅、河南、江蘇和陜西居群的所有單株,云南居群的1號(hào)單株,廣西居群的5號(hào)單株,上海居群的1號(hào)、3號(hào)和5號(hào)單株及山東居群的1~3號(hào)和5號(hào)單株,自展支持率達(dá)99%,其中安徽居群相對(duì)獨(dú)立。第4支是與其他居群遺傳關(guān)系較遠(yuǎn)的貴州居群的5個(gè)單株,自展支持率為99%,與其他居群明顯分離,但與外類群關(guān)系相對(duì)較近。

    通過(guò)比較可見(jiàn):何首烏野生居群的聚類結(jié)果與psbA-trnH序列分析的劃分結(jié)果基本一致。

    3 討論和結(jié)論

    遺傳多樣性是生物在長(zhǎng)期進(jìn)化和發(fā)展過(guò)程中形成的自然屬性,其不僅體現(xiàn)在種群間和種群內(nèi),也體現(xiàn)在不同個(gè)體間,是種群和個(gè)體遺傳變異的總和[15]。特定物種的遺傳多樣性是該物種長(zhǎng)期生存、進(jìn)化和適應(yīng)的結(jié)果[16]。一般來(lái)說(shuō),自交物種的核苷酸多樣性指數(shù)平均值為0.51,異交物種則小于0.1[17],而何首烏野生居群的核苷酸多樣性指數(shù)僅0.028 56。Buso等[18]認(rèn)為遺傳分化系數(shù)大于0.25則表示居群間的分化程度很大,而本研究中何首烏居群間的遺傳分化系數(shù)達(dá)到0.918 68。Slatkin[19-20]將基因流分為高(大于或等于1.000)、中(0.250~0.990)、低(0.000~0.249)3個(gè)等級(jí),本研究中何首烏居群間的基因流為0.04,屬于低等級(jí)。因而,總體上看何首烏野生居群整體遺傳多樣性水平偏低,居群間遺傳分化很大,但居群間的基因流遠(yuǎn)小于標(biāo)準(zhǔn),說(shuō)明遺傳漂變是影響何首烏居群遺傳結(jié)構(gòu)的主導(dǎo)因素[21]。而何首烏居群間的基因流不足以完全抵制由遺傳漂變而引起的居群分化,這也是何首烏居群遺傳多樣性水平偏低和居群間遺傳分化系數(shù)較大的原因。

    圖1 基于psbA-trnH序列分析的何首烏17個(gè)野生居群的NJ系統(tǒng)樹(shù)Fig.1 NJ phylogenetic tree of seventeen wild populations of Fallopia multiflora(Thunb.)Harald.based on psbA-trnH sequence analysis

    研究結(jié)果表明:除貴州居群外,供試的其他16個(gè)居群間的基因流為0.10,遺傳分化系數(shù)為0.837 19;而貴州居群與其周邊省區(qū)(四川、云南、廣西、湖南和湖北)居群間的基因流僅為0.03,遺傳分化系數(shù)卻達(dá)到0.937 62,明顯大于其他16個(gè)居群間。推測(cè)造成這一現(xiàn)象的原因可能是由于云貴高原具有較特殊的喀斯特地貌,由此形成的地理隔離阻礙了基因交流,導(dǎo)致分化較大。因此,何首烏居群的整體多樣性水平在很大程度上受貴州居群的影響。從形態(tài)上看,17個(gè)居群的何首烏外部形態(tài)特征基本一致,但貴州居群植株主、側(cè)葉脈的白斑程度最大,葉竇凹陷程度最小,花序長(zhǎng)度最短,葉長(zhǎng)寬比和宿存花被長(zhǎng)度相對(duì)較小,開(kāi)花期相對(duì)較晚(另文發(fā)表)。貴州居群在形態(tài)上的差異得到了psbA-trnH序列分析及聚類分析(自展支持率99%,與供試的其他16個(gè)居群明顯分離)結(jié)果的支持,證明貴州居群確實(shí)與其他居群存在較大差異,表現(xiàn)出特殊的形態(tài)和遺傳特性。推測(cè)貴州居群可能是何首烏的1個(gè)新變種,但還需多學(xué)科研究驗(yàn)證。

    聚類分析結(jié)果表明:17個(gè)何首烏野生居群中有13個(gè)居群各自的5個(gè)單株聚在一起,僅山東、云南和廣西居群各有1個(gè)單株與同居群大部分單株較遠(yuǎn)而被聚到其他分支中,而上海居群的2號(hào)和4號(hào)單株聚在第1支中,1號(hào)、3號(hào)和5號(hào)單株則聚在第3支中。說(shuō)明各居群不同單株基于psbA-trnH序列的遺傳背景基本一致,僅有少量變異,這一現(xiàn)象可能與這些居群?jiǎn)沃觊g分布的地理環(huán)境生態(tài)因子不同有關(guān)。在進(jìn)化發(fā)展過(guò)程中同種不同個(gè)體所承受的生境選擇存在一定差異,使個(gè)體在各自的自然環(huán)境下發(fā)生了變異。有研究者在黃獨(dú)(Dioscorea bulbiferaL.)遺傳多樣性研究中也觀察到這一現(xiàn)象[22],因而,同一居群不同個(gè)體間存在變異的情況并不少見(jiàn)。

    [1]WU Z Y,RAVEN P H.Flora of China:Vol.5[M].Beijing: Science Press,2003:277-350.

    [2]國(guó)家藥典委員會(huì).中華人民共和國(guó)藥典:2010年版(一部)[M].北京:中國(guó)醫(yī)藥科技出版社,2010:164-165.

    [3]王凌暉,曹福亮,汪貴斌,等.何首烏野生種質(zhì)資源的RAPD指紋圖譜構(gòu)建[J].南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2005,29 (4):37-40.

    [4]嚴(yán)寒靜,房志堅(jiān),余世孝.不同種源何首烏的ITS序列分析及其親緣關(guān)系研究[J].西北植物學(xué)報(bào),2008,28(5):922-927.

    [5]孫華欽,羅科,鄒文俊,等.穿龍薯蕷、黃山藥和盾葉薯蕷psbA-trnH片段序列分析[J].應(yīng)用與環(huán)境生物學(xué)報(bào),2006,12 (6):792-797.

    [6]邵世光,韓麗,馬艷紅,等.楓斗類石斛cpDNApsbA-trnH的序列分析與鑒別[J].藥學(xué)學(xué)報(bào),2009,44(10):1173-1178.

    [7]李永,楊攀,李永華.伏牛山區(qū)連翹遺傳多樣性研究[J].西北植物學(xué)報(bào),2011,31(4):665-670.

    [8]DOYLE J J,DOYLE J L.A rapid DNA isolation procedure for small quantities of fresh leaf tissue[J].Phytochemical Bulletin,1987,19 (1):11-15.

    [9]KUMAR S,TAMURA K,JAKOBSEN I B,et al.MEGA2: molecularevolutionarygeneticsanalysissoftware[J].Bioinformatics,2001,17(12):1244-1245.

    [10]KIMURA M.A simple method for estimating evolutionary rates of basesubstitutionsthroughcomparativestudiesofnucleotide sequences[J].Journal of Molecular Evolution,1980,16(2): 111-120.

    [12]NEI M.Molecular Evolutionary Genetics[M].New York: Columbia University Press,1987.

    [13]LYNCH M,CREASE T J.The analysis of population survey data on DNA sequence variation[J].Molecular Biology and Evolution,1990,7(4):377-394.

    [14]FELSENSTEIN J.Confidence limits on phylogenies:an approach using the bootstrap[J].Evolution,1985,39(4):783-791.

    [15]中國(guó)科學(xué)院生物多樣性委員會(huì).生物多樣性研究的原理與方法[M].北京:中國(guó)科學(xué)技術(shù)出版社,1994.

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

    [17]HAMRICK J L,GODT M J W.Allozyme diversity in plant species[M]∥BROWN A H D,CLEGG M T,KAHLER A L,et al.Plant Population Genetics,Breeding and Genetic Resources.Sunderland: Sinauer Associates,Inc.,1989:43-63.

    [18]BUSO G S C,RANGEL P H,F(xiàn)ERREIRA M E.Analysis of genetic variability in South American wild rice populations(Oryza glumaepatula)with isozymes and RAPD markers[J].Molecular Ecology,1998,7(1):107-117.

    [19]SLATKIN M.Estimating levels of gene flow in natural populations[J].Genetics,1981,99(2):323-335.

    [20]SLATKIN M.Rare alleles as indicators of gene flow[J].Evolution,1985,39(1):53-65.

    [21]WALKER C W,VILC,LANDA A,et al.Genetic variation and population structureinScandinavianwolverine(Gulogulo) populations[J].Molecular Ecology,2001,10(1):53-63.

    [22]鄭玉紅,夏冰,杭悅宇,等.黃獨(dú)遺傳多樣性研究[J].西北植物學(xué)報(bào),2006,26(10):2011-2017.

    (責(zé)任編輯:佟金鳳)

    Genetic diversity of wild populations of Fallopia multiflora based on psbA-trnH analysis

    BAI Ming-ming,SUN Xiao-qin,GUO Jian-lin,LI Mi-mi,HANG Yue-yu①(Jiangsu Province Key Laboratory for PlantEx-situConservation,Institute of Botany,Jiangsu Province and the Chinese Academy of Sciences,Nanjing 210014,China),J.Plant Resour.&Environ.2012,21(2):36-44

    Sequence ofpsbA-trnH of 85 individuals in 17 wild populations ofFallopia multiflora(Thunb.) Harald.from different provinces and regions in China was amplified and analyzed,and on the basis,genetic diversity among populations was analyzed and cluster analysis of 85 individuals was also carried out by NJ method.The results show that the length ofpsbA-trnH sequence of 85 individuals is 384 bp,in which,there are 167 bp variable sites and 53 bp parsimony informative sites,accounting for 43.5%and 13.8%of the total length of sequence,respectively.Variable types are mainly base deletion and substitution.Variable sites mainly concentrate in the region of 235-281 bp.17 populations are almostly divided into three types according to site variation status.The genetic distances among 17 populations are 0.000-0.172,in which,genetic distances between Guizhou population and other 16 populations are 0.167-0.172,and those among other 16 populations are 0.000-0.017.Nucleotide diversity index (Pi),coefficient of gene differentiation(Nst)and gene flow(Nm)among 17 populations are 0.028 56,0.918 68 and 0.04,respectively.Pi,Nst andNm among other 16 populations except Guizhou population are 0.015 68,0.837 19 and 0.10,respectively.AndPi,Nst andNm between Guizhou population and its neighboring populations(Sichuan,Yunnan,Guangxi,Hu’nan and Hubei)are 0.047 99,0.937 62 and 0.03,respectively.On NJ phylogenetic tree,17 populations are clustered into four branches and individuals tested in most populations are clustered in a same branch and only Guizhoupopulation is clustered alone in a branch,which is basically same with the deviation result by sequence analysis.It is suggested that 91.868%of overall genetic variation of 17 wild populations exists among populations and 8.132%within populations,and gene exchange among populations is less.Except Guizhou population,overall genetic diversity level among other 16 populations is low,indicating that the overall diversity level ofF.multiflorapopulations is affected by Guizhou population at a large extent.

    Fallopia multiflora(Thunb.)Harald.;population;psbA-trnH sequence;genetic diversity; NJ phylogenetic tree

    book=2012,ebook=77

    Q946-33;S567.23+9

    A

    1674-7895(2012)02-0036-09

    2011-12-27

    國(guó)家農(nóng)業(yè)部DUS測(cè)試品種信息DNA測(cè)試技術(shù)研究項(xiàng)目(200903008-02-05)

    白明明(1987—),女,山東郯城人,碩士研究生,主要從事藥用植物資源的研究。

    ①通信作者E-mail:hangyueyu@21cn.com

    猜你喜歡
    居群何首烏堿基
    燕山地區(qū)不同居群白頭翁種子萌發(fā)的初步評(píng)價(jià)
    制何首烏引起藥品不良反應(yīng)1例報(bào)道
    基于簡(jiǎn)單重復(fù)序列間擴(kuò)增分子標(biāo)記的金釵石斛遺傳多樣性研究
    應(yīng)用思維進(jìn)階構(gòu)建模型 例談培養(yǎng)學(xué)生創(chuàng)造性思維
    中國(guó)科學(xué)家創(chuàng)建出新型糖基化酶堿基編輯器
    生命“字母表”迎來(lái)4名新成員
    生命“字母表”迎來(lái)4名新成員
    3個(gè)居群野生草地早熟禾耐鹽性比較研究
    香菇何首烏降血脂
    制何首烏中二苯乙烯苷對(duì)光和熱的不穩(wěn)定性
    中成藥(2014年11期)2014-02-28 22:29:49
    美女高潮喷水抽搐中文字幕| 精品无人区乱码1区二区| 制服丝袜大香蕉在线| 日韩欧美免费精品| 天美传媒精品一区二区| 天美传媒精品一区二区| 色精品久久人妻99蜜桃| 久久性视频一级片| 欧美色视频一区免费| 51国产日韩欧美| 老司机午夜十八禁免费视频| 久久久久性生活片| 特大巨黑吊av在线直播| 婷婷丁香在线五月| 天堂影院成人在线观看| 国产精品电影一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲专区国产一区二区| 欧美日韩黄片免| 久久久久久九九精品二区国产| 一a级毛片在线观看| 少妇的逼好多水| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久视频播放| 99国产精品一区二区蜜桃av| 99riav亚洲国产免费| 桃色一区二区三区在线观看| 中文字幕精品亚洲无线码一区| 国产97色在线日韩免费| 亚洲国产中文字幕在线视频| 黄色丝袜av网址大全| 淫妇啪啪啪对白视频| 人人妻人人看人人澡| 久久精品国产亚洲av涩爱 | 亚洲人成电影免费在线| 免费在线观看日本一区| 午夜福利在线观看免费完整高清在 | 天天一区二区日本电影三级| 草草在线视频免费看| 国产91精品成人一区二区三区| 天天添夜夜摸| 精品免费久久久久久久清纯| 嫁个100分男人电影在线观看| 老司机午夜十八禁免费视频| 禁无遮挡网站| 国模一区二区三区四区视频| 久9热在线精品视频| 99热6这里只有精品| 国产v大片淫在线免费观看| av黄色大香蕉| 窝窝影院91人妻| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩精品网址| 人妻夜夜爽99麻豆av| 搡老岳熟女国产| 搡老岳熟女国产| 国产97色在线日韩免费| 很黄的视频免费| 午夜激情欧美在线| 又爽又黄无遮挡网站| 国产在视频线在精品| 91在线观看av| 色av中文字幕| 三级毛片av免费| 亚洲av电影在线进入| 99久久精品热视频| 五月伊人婷婷丁香| 99在线人妻在线中文字幕| 岛国在线免费视频观看| 性色avwww在线观看| 国产精品久久久久久久久免 | 1000部很黄的大片| 国内精品一区二区在线观看| 啦啦啦观看免费观看视频高清| 一本综合久久免费| 亚洲在线观看片| 亚洲国产精品成人综合色| 久久国产精品影院| 三级毛片av免费| 国产私拍福利视频在线观看| 搡老熟女国产l中国老女人| 十八禁人妻一区二区| 无遮挡黄片免费观看| 久久久精品大字幕| 亚洲av第一区精品v没综合| 身体一侧抽搐| 色av中文字幕| 欧美性猛交黑人性爽| 又紧又爽又黄一区二区| 久久久久性生活片| 国产精品免费一区二区三区在线| 91在线观看av| 亚洲成av人片免费观看| av国产免费在线观看| 精品国内亚洲2022精品成人| 女人十人毛片免费观看3o分钟| 亚洲中文字幕日韩| x7x7x7水蜜桃| 国产一区二区在线观看日韩 | 亚洲不卡免费看| 亚洲精品粉嫩美女一区| 91在线观看av| 高清毛片免费观看视频网站| 中文字幕av成人在线电影| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 麻豆一二三区av精品| 夜夜爽天天搞| 亚洲 国产 在线| 欧美极品一区二区三区四区| 国产高清视频在线观看网站| 伊人久久精品亚洲午夜| 一本一本综合久久| 精品久久久久久,| 日本一本二区三区精品| 国内毛片毛片毛片毛片毛片| 不卡一级毛片| 久久久久精品国产欧美久久久| 国产黄色小视频在线观看| a在线观看视频网站| 最好的美女福利视频网| 国产精品电影一区二区三区| 日韩欧美精品免费久久 | 欧美日韩综合久久久久久 | 特大巨黑吊av在线直播| 99国产综合亚洲精品| 精品人妻偷拍中文字幕| 欧美大码av| 日韩av在线大香蕉| 午夜福利18| 看黄色毛片网站| 好男人在线观看高清免费视频| 欧美午夜高清在线| 欧美日韩瑟瑟在线播放| 欧美色欧美亚洲另类二区| 床上黄色一级片| 1000部很黄的大片| 日本精品一区二区三区蜜桃| 一区二区三区国产精品乱码| 亚洲最大成人手机在线| 久久久成人免费电影| 亚洲欧美日韩东京热| 国产亚洲欧美98| 综合色av麻豆| 中文字幕av成人在线电影| 99精品在免费线老司机午夜| 毛片女人毛片| 18禁国产床啪视频网站| 一级黄色大片毛片| 高清在线国产一区| 夜夜爽天天搞| 国产主播在线观看一区二区| 少妇的逼好多水| 成年人黄色毛片网站| 国产欧美日韩一区二区精品| 日本一二三区视频观看| 最后的刺客免费高清国语| 中亚洲国语对白在线视频| 老司机午夜十八禁免费视频| 国产伦人伦偷精品视频| 欧美bdsm另类| 日韩精品青青久久久久久| 午夜免费男女啪啪视频观看 | 高清日韩中文字幕在线| 夜夜看夜夜爽夜夜摸| 九色国产91popny在线| 久久性视频一级片| tocl精华| 特大巨黑吊av在线直播| 蜜桃久久精品国产亚洲av| 蜜桃久久精品国产亚洲av| 日本在线视频免费播放| 99热这里只有精品一区| 久久伊人香网站| 啦啦啦观看免费观看视频高清| 老司机午夜十八禁免费视频| 日日干狠狠操夜夜爽| 一区福利在线观看| 精品不卡国产一区二区三区| 观看免费一级毛片| av女优亚洲男人天堂| 色综合站精品国产| 国产av在哪里看| 少妇熟女aⅴ在线视频| 国产精品久久电影中文字幕| 国产成人av激情在线播放| 亚洲 欧美 日韩 在线 免费| 一a级毛片在线观看| 欧美性猛交╳xxx乱大交人| 色综合亚洲欧美另类图片| 999久久久精品免费观看国产| 免费一级毛片在线播放高清视频| e午夜精品久久久久久久| 搡老妇女老女人老熟妇| 久久久成人免费电影| 老汉色∧v一级毛片| 精品一区二区三区视频在线 | 欧美黄色淫秽网站| 亚洲欧美日韩卡通动漫| 久久中文看片网| 日韩人妻高清精品专区| 熟妇人妻久久中文字幕3abv| 欧美极品一区二区三区四区| a在线观看视频网站| 亚洲av日韩精品久久久久久密| 亚洲成人精品中文字幕电影| 欧美最新免费一区二区三区 | 国产一区二区三区在线臀色熟女| 欧美av亚洲av综合av国产av| h日本视频在线播放| 国产又黄又爽又无遮挡在线| 最近视频中文字幕2019在线8| 12—13女人毛片做爰片一| 99久国产av精品| 午夜免费激情av| 国产一区二区三区在线臀色熟女| 午夜a级毛片| 99久久久亚洲精品蜜臀av| 又黄又粗又硬又大视频| 99久久成人亚洲精品观看| 国产精品美女特级片免费视频播放器| 亚洲av电影在线进入| 搡女人真爽免费视频火全软件 | 精品一区二区三区视频在线 | 亚洲在线自拍视频| 97超级碰碰碰精品色视频在线观看| 久久亚洲真实| 国产成人福利小说| 午夜福利欧美成人| 欧美一区二区精品小视频在线| 少妇熟女aⅴ在线视频| 久久久精品大字幕| 中出人妻视频一区二区| 久久九九热精品免费| 听说在线观看完整版免费高清| av欧美777| 日本免费一区二区三区高清不卡| 久久久久国产精品人妻aⅴ院| 国产色婷婷99| 亚洲av不卡在线观看| av女优亚洲男人天堂| 国产精品久久久久久人妻精品电影| 少妇的逼好多水| 搡女人真爽免费视频火全软件 | 亚洲人成电影免费在线| 日韩av在线大香蕉| 丰满乱子伦码专区| 色综合亚洲欧美另类图片| 国产av在哪里看| 精品国产亚洲在线| 天天一区二区日本电影三级| 一级毛片女人18水好多| 精品久久久久久成人av| 美女 人体艺术 gogo| 99riav亚洲国产免费| 天美传媒精品一区二区| 19禁男女啪啪无遮挡网站| 国产精品香港三级国产av潘金莲| 日韩人妻高清精品专区| 又紧又爽又黄一区二区| 特级一级黄色大片| 窝窝影院91人妻| 日韩欧美国产在线观看| 亚洲精品456在线播放app | 免费av观看视频| 午夜两性在线视频| 91麻豆精品激情在线观看国产| 在线观看美女被高潮喷水网站 | 国产免费av片在线观看野外av| 国产精品日韩av在线免费观看| 欧美成人a在线观看| 51午夜福利影视在线观看| 亚洲性夜色夜夜综合| 哪里可以看免费的av片| 香蕉av资源在线| 一a级毛片在线观看| 亚洲自拍偷在线| 日本撒尿小便嘘嘘汇集6| 十八禁人妻一区二区| 亚洲aⅴ乱码一区二区在线播放| 黄色视频,在线免费观看| 亚洲精品粉嫩美女一区| 欧美黄色片欧美黄色片| 久久久成人免费电影| 男插女下体视频免费在线播放| 欧美极品一区二区三区四区| 国产亚洲精品综合一区在线观看| 欧美区成人在线视频| 日韩成人在线观看一区二区三区| 五月玫瑰六月丁香| 99热6这里只有精品| 亚洲av成人精品一区久久| av福利片在线观看| 亚洲激情在线av| 日韩欧美精品免费久久 | 有码 亚洲区| 99在线人妻在线中文字幕| 99精品欧美一区二区三区四区| 国产91精品成人一区二区三区| 国产蜜桃级精品一区二区三区| 久久亚洲精品不卡| 在线观看一区二区三区| 美女免费视频网站| 搡老岳熟女国产| 一级a爱片免费观看的视频| 亚洲乱码一区二区免费版| 欧美+亚洲+日韩+国产| 欧美+亚洲+日韩+国产| 久久久久久九九精品二区国产| 麻豆久久精品国产亚洲av| 国产极品精品免费视频能看的| 国产精品亚洲一级av第二区| 波野结衣二区三区在线 | 老熟妇仑乱视频hdxx| 欧美一区二区精品小视频在线| 成人av一区二区三区在线看| 日韩欧美 国产精品| 午夜福利高清视频| 亚洲 国产 在线| 又粗又爽又猛毛片免费看| 日韩精品青青久久久久久| 日韩精品青青久久久久久| 欧美又色又爽又黄视频| 婷婷精品国产亚洲av在线| 97碰自拍视频| 精品人妻1区二区| 午夜影院日韩av| 搡老岳熟女国产| 在线观看一区二区三区| 国产伦精品一区二区三区视频9 | 国产精品日韩av在线免费观看| 久久精品综合一区二区三区| 露出奶头的视频| 国产精品女同一区二区软件 | 亚洲avbb在线观看| 国产精品 欧美亚洲| 久久久久精品国产欧美久久久| 亚洲天堂国产精品一区在线| 最新美女视频免费是黄的| 天堂网av新在线| 成人一区二区视频在线观看| 波多野结衣巨乳人妻| 在线观看舔阴道视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品乱码一区二三区的特点| 少妇高潮的动态图| 国产精品亚洲一级av第二区| 999久久久精品免费观看国产| 国产爱豆传媒在线观看| 国产av在哪里看| 色噜噜av男人的天堂激情| 香蕉av资源在线| 国产一区二区在线av高清观看| 亚洲av电影在线进入| 日本a在线网址| 色精品久久人妻99蜜桃| 好男人在线观看高清免费视频| 黄色片一级片一级黄色片| 国产单亲对白刺激| 成人18禁在线播放| 欧美一级a爱片免费观看看| 日日摸夜夜添夜夜添小说| 欧美+日韩+精品| 在线观看66精品国产| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品一区av在线观看| 少妇丰满av| 99精品欧美一区二区三区四区| 九九热线精品视视频播放| 欧美日韩乱码在线| 久久久久国内视频| 一级作爱视频免费观看| 狂野欧美激情性xxxx| 中文在线观看免费www的网站| 国产精品野战在线观看| 天天一区二区日本电影三级| 久久香蕉精品热| 欧美又色又爽又黄视频| 搡女人真爽免费视频火全软件 | 香蕉av资源在线| 国内精品久久久久久久电影| 在线观看66精品国产| 俄罗斯特黄特色一大片| 国产av一区在线观看免费| 成人欧美大片| 又黄又爽又免费观看的视频| 国产午夜福利久久久久久| 男人舔奶头视频| 久久精品国产亚洲av香蕉五月| 九色成人免费人妻av| 欧美激情在线99| 成人性生交大片免费视频hd| 女人被狂操c到高潮| 九色国产91popny在线| 香蕉久久夜色| 国产精品久久久久久久电影 | 最后的刺客免费高清国语| 免费在线观看亚洲国产| 国产精品永久免费网站| 麻豆一二三区av精品| h日本视频在线播放| 欧美+亚洲+日韩+国产| 九九久久精品国产亚洲av麻豆| 色综合站精品国产| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美精品综合久久99| 麻豆久久精品国产亚洲av| 嫩草影院精品99| 婷婷丁香在线五月| 久久国产精品影院| 国产精品一区二区三区四区久久| 成人性生交大片免费视频hd| 18+在线观看网站| 在线观看一区二区三区| 久久久国产精品麻豆| 毛片女人毛片| 亚洲不卡免费看| 无限看片的www在线观看| 老熟妇仑乱视频hdxx| 91麻豆精品激情在线观看国产| 国产伦人伦偷精品视频| 精品熟女少妇八av免费久了| av黄色大香蕉| 欧美日韩瑟瑟在线播放| 波多野结衣高清无吗| 中文字幕人妻丝袜一区二区| 老鸭窝网址在线观看| aaaaa片日本免费| 啦啦啦观看免费观看视频高清| 国产黄片美女视频| 高清日韩中文字幕在线| 亚洲第一欧美日韩一区二区三区| 最好的美女福利视频网| 欧美激情在线99| 国产精品一区二区免费欧美| 不卡一级毛片| 18禁在线播放成人免费| 国产v大片淫在线免费观看| av福利片在线观看| 中文字幕人妻丝袜一区二区| 非洲黑人性xxxx精品又粗又长| 老司机在亚洲福利影院| 日韩欧美精品v在线| 亚洲性夜色夜夜综合| 女警被强在线播放| 高潮久久久久久久久久久不卡| 在线观看免费视频日本深夜| 日本 av在线| 日日干狠狠操夜夜爽| 色视频www国产| 午夜免费男女啪啪视频观看 | 免费观看精品视频网站| 日韩中文字幕欧美一区二区| 男女下面进入的视频免费午夜| 日韩欧美 国产精品| 免费电影在线观看免费观看| 欧美激情在线99| www国产在线视频色| 欧美日韩国产亚洲二区| 波多野结衣高清无吗| 久久精品综合一区二区三区| 国产淫片久久久久久久久 | 国产精品三级大全| av女优亚洲男人天堂| 欧美一级a爱片免费观看看| 日韩中文字幕欧美一区二区| 亚洲激情在线av| 一边摸一边抽搐一进一小说| 精品99又大又爽又粗少妇毛片 | www.色视频.com| 久久天躁狠狠躁夜夜2o2o| 国产精品 欧美亚洲| 国产精品98久久久久久宅男小说| avwww免费| 国产精品 国内视频| 亚洲avbb在线观看| 少妇高潮的动态图| 久久精品夜夜夜夜夜久久蜜豆| 91久久精品国产一区二区成人 | 免费在线观看日本一区| av在线蜜桃| 国产色婷婷99| 亚洲精品在线美女| 久久久久久久久久黄片| 毛片女人毛片| 啦啦啦韩国在线观看视频| 久久久久久人人人人人| 色综合站精品国产| 国产视频内射| 一个人看视频在线观看www免费 | 久久午夜亚洲精品久久| 国产探花在线观看一区二区| 国产极品精品免费视频能看的| 亚洲精品美女久久久久99蜜臀| 午夜精品久久久久久毛片777| 亚洲av中文字字幕乱码综合| 精品人妻1区二区| 性色av乱码一区二区三区2| 精品午夜福利视频在线观看一区| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 午夜老司机福利剧场| 国产色婷婷99| 亚洲精品亚洲一区二区| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 国产精品98久久久久久宅男小说| 国产精品一及| 真人一进一出gif抽搐免费| 国产激情偷乱视频一区二区| 99国产极品粉嫩在线观看| 日本免费a在线| 亚洲精品乱码久久久v下载方式 | 欧美性猛交╳xxx乱大交人| www国产在线视频色| 三级男女做爰猛烈吃奶摸视频| 岛国视频午夜一区免费看| 日韩有码中文字幕| 又黄又粗又硬又大视频| 婷婷六月久久综合丁香| tocl精华| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜| 桃红色精品国产亚洲av| 日日摸夜夜添夜夜添小说| 日韩成人在线观看一区二区三区| 老鸭窝网址在线观看| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av在线| 久久久久久久精品吃奶| 99热这里只有是精品50| 国产av在哪里看| 色精品久久人妻99蜜桃| 美女被艹到高潮喷水动态| 国产真实伦视频高清在线观看 | 亚洲在线观看片| 国产野战对白在线观看| 国产精品 国内视频| 精品国内亚洲2022精品成人| 日本与韩国留学比较| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 美女大奶头视频| 久久九九热精品免费| 国产伦一二天堂av在线观看| 他把我摸到了高潮在线观看| 无人区码免费观看不卡| 日韩精品中文字幕看吧| 波多野结衣高清作品| 久久草成人影院| 欧美日韩黄片免| 精品免费久久久久久久清纯| 日本一本二区三区精品| 在线视频色国产色| 欧美成狂野欧美在线观看| 老司机深夜福利视频在线观看| 97超视频在线观看视频| 欧美3d第一页| 国产99白浆流出| 欧美成人性av电影在线观看| 国产探花在线观看一区二区| 欧美bdsm另类| 热99re8久久精品国产| 在线观看66精品国产| 最近视频中文字幕2019在线8| 一本精品99久久精品77| 国产伦精品一区二区三区视频9 | 久久国产精品影院| 久久亚洲真实| 成年版毛片免费区| 欧美zozozo另类| 又黄又爽又免费观看的视频| 日本免费a在线| 欧美乱色亚洲激情| 免费大片18禁| 欧美日韩国产亚洲二区| 国产欧美日韩精品亚洲av| 色尼玛亚洲综合影院| 日韩成人在线观看一区二区三区| 国产淫片久久久久久久久 | av女优亚洲男人天堂| 一级黄色大片毛片| 亚洲七黄色美女视频| 中文字幕av在线有码专区| 日韩欧美三级三区| 熟女少妇亚洲综合色aaa.| 女人高潮潮喷娇喘18禁视频| 亚洲欧美日韩高清专用| 在线观看66精品国产| 搡老妇女老女人老熟妇| 久久久久久久久中文| 免费高清视频大片| 国内久久婷婷六月综合欲色啪| 午夜影院日韩av| bbb黄色大片| 国产极品精品免费视频能看的| 啪啪无遮挡十八禁网站| 麻豆成人av在线观看| 国产精品嫩草影院av在线观看 | 嫩草影院精品99| 一本综合久久免费| 观看美女的网站| 在线十欧美十亚洲十日本专区| 岛国在线免费视频观看| 成年版毛片免费区| 一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 精品欧美国产一区二区三| 中文亚洲av片在线观看爽| 中亚洲国语对白在线视频| 国产毛片a区久久久久| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www| 午夜福利免费观看在线| 国产激情偷乱视频一区二区|