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

    鉤齒溲疏的SRAP遺傳多樣性分析

    2014-12-23 11:00:07呂俊杰任瑩徐秀榮梁婷婷臧德奎
    山東農(nóng)業(yè)科學(xué) 2014年11期
    關(guān)鍵詞:遺傳多樣性

    呂俊杰+任瑩+徐秀榮+梁婷婷+臧德奎

    摘 要:運(yùn)用SRAP分子標(biāo)記研究了山東省3個(gè)鉤齒溲疏居群的遺傳多樣性,結(jié)果顯示,用15對(duì)引物共檢測(cè)到244個(gè)位點(diǎn),其中多態(tài)位點(diǎn)167個(gè),物種水平多態(tài)位點(diǎn)百分率(PPL) 68.44%,Neis基因多樣性指數(shù)(H) 0.2158,Shannons信息指數(shù)(I)0.3282,數(shù)據(jù)表明鉤齒溲疏有較高的遺傳多樣性水平;居群間遺傳分化系數(shù)(Gst)為0.3685,表明居群間遺傳變異只占36.85%,遠(yuǎn)低于居群內(nèi)遺傳分化;在居群水平上,以嶗山鉤齒溲疏居群的遺傳多樣性最高,徂徠山最低;居群間遺傳一致度(GI)和遺傳距離(GD)變化范圍分別為0.8350~0.8884和0.1184~0.1803,居群間的遺傳距離與地理位置間沒有直接相關(guān)性。

    關(guān)鍵詞:鉤齒溲疏;SRAP;遺傳多樣性

    中圖分類號(hào):S685.990.1 文獻(xiàn)標(biāo)識(shí)號(hào):A 文章編號(hào):1001-4942(2014)11-0018-05

    鉤齒溲疏(Deutzia baroniana)屬于虎耳草科(Saxifragaceae)溲疏屬(Deutzia),落葉灌木,分布于遼寧、河北、山西、陜西、山東、江蘇和河南[1]。鉤齒溲疏花朵潔白,繁密而素凈,花期長(zhǎng),且盛開于初夏,正值少花時(shí)節(jié),宜叢植于山坡、草坪、路旁或林緣,也是巖石園、花籬的常用材料,花枝可作插花觀賞。鉤齒溲疏適應(yīng)性非常強(qiáng),常生于溝谷、林緣或巖石縫中,耐旱、耐瘠薄、抗寒,對(duì)于干旱、少雨的北方城市綠化極為適宜,是北方十分難得的園林綠化、美化材料。此外,鉤齒溲疏根、葉、果還可藥用。因此,鉤齒溲疏具有很高的觀賞和經(jīng)濟(jì)價(jià)值,是一種值得開發(fā)的野生植物資源,值得大力引種并繁育推廣。

    一個(gè)物種之所以能繁殖存活,并且適應(yīng)環(huán)境改變,根本原因在于其遺傳多樣性[2]。居群對(duì)環(huán)境的適應(yīng)能力在一定程度上由居群的遺傳多樣性水平制約著,因此通過研究其遺傳多樣性水平可預(yù)測(cè)這個(gè)居群未來的發(fā)展趨勢(shì)[3]。而目前尚未見從DNA水平上對(duì)鉤齒溲疏的遺傳多樣性進(jìn)行分析的研究和報(bào)道。

    相關(guān)序列擴(kuò)增多態(tài)性(Sequence-related amplified polymorphism,SRAP)是通過設(shè)計(jì)獨(dú)特的引物對(duì)基因ORFs(Open reading frames)的特定區(qū)域進(jìn)行擴(kuò)增并分析的標(biāo)記手段[4]。目前已成功應(yīng)用于作物遺傳多樣性評(píng)價(jià)、指紋圖譜構(gòu)建、重要性狀的標(biāo)記以及相關(guān)基因的克隆等方面[5~8]。本研究對(duì)山東省內(nèi)21份鉤齒溲疏進(jìn)行SRAP遺傳多樣性分析,旨在研究其遺傳多樣性水平,為開發(fā)利用鉤齒溲疏野生資源、在景觀應(yīng)用中選育優(yōu)良品種提供分子水平的依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)樣品于2012~2013年采自山東嶗山、泰山和徂徠山,共21份樣品(表1)。同一居群相鄰植株采樣間距30 m以上,采取新鮮幼嫩葉片放入密封袋中,用變色硅膠干燥。記錄所采樣品的小地名、海拔、經(jīng)緯度并編號(hào)。

    1.2 DNA提取

    選用改良的CTAB法[9]提取樣品總DNA,用0.8%瓊脂糖凝膠電泳檢測(cè)DNA質(zhì)量,用紫外分光光度計(jì)檢測(cè)DNA濃度后用ddH2O稀釋至50 ng/μL,-20℃保存?zhèn)溆谩?/p>

    1.3 SRAP-PCR反應(yīng)體系和程序

    SRAP引物設(shè)計(jì)參照Li等[4],13個(gè)正向引物和19個(gè)反向引物共組成247對(duì)引物組合,由上海生工生物工程技術(shù)服務(wù)有限公司合成。從中篩選出15對(duì)擴(kuò)增條帶數(shù)量多、清晰、穩(wěn)定的引物組合(表2)用于正式擴(kuò)增。

    反應(yīng)體系為25 μL:10×PCR Buffer 2.5 μL,模板DNA 75 ng,Mg2+ 3 mmol/L,dNTP 0.4 mmol/L,上下游引物各0.22 μmol,Taq DNA聚合酶1.5 U,用ddH2O調(diào)整使終體積為25 μL。

    PCR擴(kuò)增程序?yàn)椋?4℃預(yù)變性5 min;94℃變性1 min,35℃退火1 min,72℃延伸1 min,5個(gè)循環(huán);94℃變性1 min,50℃退火1 min,72℃延伸1 min,30個(gè)循環(huán);72℃延伸7 min。

    PCR擴(kuò)增產(chǎn)物采用2%瓊脂糖凝膠電泳分離,溴化乙錠(EB)染色后在北京六一WD-9403CS型紫外儀上采集圖像。

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

    根據(jù)擴(kuò)增條帶結(jié)果,選取易于識(shí)別的條帶,在相同遷移率位置上,有記為1,反之為0,得到原始數(shù)據(jù)矩陣。用POPGENE 1.32統(tǒng)計(jì)各項(xiàng)遺傳參數(shù):多態(tài)性位點(diǎn)數(shù)NPL、多態(tài)位點(diǎn)百分率PPL、等位基因數(shù)Na、有效等位基因數(shù)Ne;Neis基因多樣性指數(shù)H、Shannons信息指數(shù)I;居群間Neis遺傳一致度GI和遺傳距離GD;總遺傳多樣性Ht、群體內(nèi)的遺傳多樣性Hs、群體間遺傳分化系數(shù)Gst=Dst/Ht, Dst為群體間遺傳多樣性;體現(xiàn)基因流強(qiáng)度的群體每代遷移數(shù)(Nm)。應(yīng)用NTSYS pc version 2.10e 軟件根據(jù)Neis遺傳相似系數(shù)利用UPGMA方法進(jìn)行聚類,繪出聚類分析樹狀圖[10]。

    2 結(jié)果與分析

    2.1 SRAP引物擴(kuò)增多態(tài)性分析

    鉤齒溲疏21份樣品用篩選出的15對(duì)引物擴(kuò)增(表3),共擴(kuò)增出244個(gè)位點(diǎn),多態(tài)性位點(diǎn)167個(gè),多態(tài)性位點(diǎn)百分率為68.44%,說明山東省鉤齒溲疏遺傳多樣性較高。每對(duì)引物產(chǎn)生的位點(diǎn)數(shù)為12~24個(gè),多態(tài)性位點(diǎn)數(shù)為7~17個(gè),平均每對(duì)引物產(chǎn)生16.27個(gè)位點(diǎn)和11.13個(gè)多態(tài)性位點(diǎn)。以ME7-EM15的多態(tài)性位點(diǎn)比率最高,為100%,H值和I值也均最高;ME10-EM15產(chǎn)生的多態(tài)位點(diǎn)比率最低,為38.89%,H值和I值也最低。

    2.2 居群遺傳多樣性分析

    在物種水平上,多態(tài)位點(diǎn)百分率(PPL)為68.44%。在居群水平上,多態(tài)位點(diǎn)百分率(PPL)介于32.79%~45.08%,平均為37.57%(表4);以嶗山遺傳多樣性最高,泰山次之,徂徠山最低。由Shannons信息指數(shù)計(jì)算的鉤齒溲疏物種水平的遺傳多樣性為0.3282,各居群的遺傳多樣性變化范圍為0.1776~0.2525,平均為0.2046。由Neis指數(shù)計(jì)算的物種水平的遺傳多樣性為0.2158,各居群基因多樣性介于0.1200~0.1719,平均為0.1383。該指數(shù)計(jì)算的值較Shannons信息指數(shù)低,但兩種方法計(jì)算的三個(gè)居群的遺傳多樣性大小同多態(tài)位點(diǎn)百分率大小順序是一致的,均是嶗山最高,徂徠山最低。數(shù)據(jù)表明山東省鉤齒溲疏物種水平的遺傳多樣性較高,居群水平的遺傳多樣性較低。

    2.3 居群遺傳結(jié)構(gòu)與遺傳一致度分析

    對(duì)鉤齒溲疏3個(gè)自然居群的遺傳結(jié)構(gòu)分析表明,總遺傳多樣性(Ht)為0.2190,其中居群內(nèi)(Hs)為0.1383,居群間(Dst)為0.0807,居群內(nèi)遺傳多樣度明顯高于居群間。居群間遺傳分化系數(shù)(Gst)為0.3685,即居群間遺傳變異僅占36.85%,因此居群內(nèi)遺傳分化遠(yuǎn)高于居群間。在居群遺傳學(xué)研究中,基因流被劃分為高、中、低三水平[11],而鉤齒溲疏各居群間基因流(Nm)大小為0.8566,介于0.25~0.99,屬于中等水平,表明居群間基因流相對(duì)有限。

    統(tǒng)計(jì)鉤齒溲疏各居群間的Neis遺傳距離GD和遺傳一致度GI以進(jìn)一步分析鉤齒溲疏居群的遺傳結(jié)構(gòu)。3個(gè)居群間遺傳距離GD變動(dòng)在0.1184~0.1803,平均為0.1514(表5),以嶗山與泰山的遺傳距離最近,嶗山與徂徠山的最遠(yuǎn)。3個(gè)居群間遺傳一致度GI變化范圍為0.8350~0.8884,平均為0.8598。泰山與嶗山居群的遺傳距離是3個(gè)居群間的最小值,遺傳一致度最高,但在地理位置上距離較遠(yuǎn);同樣,泰山與徂徠山地理距離最近,但遺傳一致度反比泰山與嶗山的低。這表明,鉤齒溲疏居群間的遺傳距離和地理位置并非直接相關(guān)。

    2.4 聚類分析

    UPGMA聚類結(jié)果(圖1)顯示,泰山居群14號(hào)和15號(hào)樣品遺傳相似度最高,為0.91。在相似系數(shù)為0.74時(shí)所有樣品可以分為3大類:第Ⅰ大類包含7個(gè)樣品,全部來自徂徠山;第Ⅱ大類包括來自嶗山的4個(gè)樣品和來自泰山的8個(gè)樣品;第Ⅲ大類只有2個(gè)樣品,均來自嶗山。第Ⅱ大類在相似系數(shù)為0.81處又可分為3個(gè)亞類,第1亞類含1個(gè)樣品,第2亞類含3個(gè)樣品,均來自嶗山;第3亞類含8個(gè)樣品,全部來自泰山。

    由聚類圖可看出,徂徠山的樣品單獨(dú)聚為一個(gè)大類,泰山的樣品聚在一個(gè)亞類里,只有嶗山的樣品分布在2個(gè)大類中,這反映出嶗山居群豐富的遺傳多樣性。同時(shí),徂徠山樣品單獨(dú)聚為第Ⅰ大類,泰山和嶗山的部分樣品聚為第Ⅱ大類,表明了嶗山與泰山的遺傳相似度較高。

    聚類結(jié)果顯示,多數(shù)情況下樣品來源決定樣品間的親緣關(guān)系,地理來源相同的樣品傾向于聚在一起,但供試樣品的聚類結(jié)果又并非嚴(yán)格按照地理來源來劃分。說明山東鉤齒溲疏的遺傳分化并非受地理來源單一因素的影響,是多個(gè)因素綜合作用的結(jié)果。

    3 討論

    在體現(xiàn)居群內(nèi)變異水平方面,多態(tài)位點(diǎn)百分率是重要指標(biāo)之一。在物種水平上鉤齒溲疏的多態(tài)位點(diǎn)百分率為68.44%,表明山東省內(nèi)鉤齒溲疏有較豐富的遺傳多樣性。

    Loveless等[12]曾統(tǒng)計(jì)過不同類型植物的遺傳變異水平和居群分化程度,得出了大范圍連續(xù)分布的異交植物的遺傳變異大部分存在于居群內(nèi)的結(jié)論;而鉤齒溲疏居群間遺傳分化系數(shù)(Gst)為0.3685,表明遺傳分化主要存在于居群內(nèi),同這一結(jié)論類似。鉤齒溲疏各居群間基因流(Nm)大小為0.8566,表明居群間基因流相對(duì)有限,是3個(gè)居群間存在一定程度的遺傳分化的原因之一。

    在居群水平上,徂徠山的遺傳多樣性最低,應(yīng)與徂徠山鉤齒溲疏生境類型變化相對(duì)有限有關(guān)。而嶗山居群的遺傳多樣性最為豐富,可能與嶗山的地理、氣候條件更為優(yōu)越和復(fù)雜有關(guān)。嶗山地處山東半島,東、南、北三面環(huán)海,西面與陸地相接。由于受海洋影響,氣候溫暖濕潤(rùn),南坡和北坡氣候迥異。嶗山植被類型豐富,群落組成上南坡?lián)碛休^多的亞熱帶區(qū)系成分,而北坡則以長(zhǎng)白區(qū)系成分為主[13]。嶗山復(fù)雜多樣的生境類型,高度的空間異質(zhì)性使得鉤齒溲疏不同基因型的植株都可以生長(zhǎng),所以遺傳多樣性處于較高水平。生態(tài)因子可作為控制居群中基因型頻率的一種選擇壓力,在相當(dāng)程度上對(duì)居群的遺傳分化起作用[14]。反過來,遺傳多樣性也是物種適應(yīng)環(huán)境變化的基礎(chǔ)。因此可以推測(cè),嶗山鉤齒溲疏居群較高的遺傳多樣性是對(duì)嶗山優(yōu)越、復(fù)雜的生態(tài)環(huán)境長(zhǎng)期適應(yīng)的結(jié)果。

    在調(diào)查中發(fā)現(xiàn),鉤齒溲疏嶗山居群豐富的遺傳多樣性不僅體現(xiàn)在分子水平上,也反映在形態(tài)上,尤其是花朵的大小,葉片的大小、形狀等。形態(tài)上豐富的變異為開發(fā)利用鉤齒溲疏野生資源、選育優(yōu)良類型提供了基礎(chǔ)。

    參 考 文 獻(xiàn):

    [1] 中國(guó)植物志編委會(huì). 中國(guó)植物志[M]. 北京: 科學(xué)出版社, 1982: 103-104.

    [2] Spiess E B. Genes in population[M]. 2nd ed. New York: John Wiley & Sons, 1989.

    [3] 潘瑩, 趙桂仿. 分子水平的遺傳多樣性及其測(cè)量方法[J]. 西北植物學(xué)報(bào), 1998, 18(4): 645-653.

    [4] Li G, Qurios C F. Sequence related amplified polymorphism(SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica[J]. Theoretical and Applied Genetics, 2001, 103(3): 455-461.

    [5] 張安世, 刑智峰, 劉永英, 等. SRAP分子標(biāo)記及其應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué), 2007, 35(9): 2562-2563.

    [6] 李莉, 彭建營(yíng), 白瑞霞, 等. SRAP與TRAP標(biāo)記及其在園藝植物研究中的應(yīng)用[J]. 西北植物學(xué)報(bào), 2006, 26(8): 1749-1752.

    [7] 張鮮艷, 張飛, 陳發(fā)棣, 等. 12 份不同地理居群野菊的遺傳多樣性分析[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 34(3): 48-54.

    [8] 林萍, 姚小華, 王開良, 等. 油茶長(zhǎng)林系列優(yōu)良無性系的SRAP分子鑒別及遺傳分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2010, 18(2): 272-279.

    [9] Kidwell K K, Osborn T C. Simple plant DNA isolation procedures[M]//Beckman J, Osborn T C. Plant genomes: methods for genetic and physical mapping. Dordrecht: Kluwer Academic Publishers, 1992: 1-13.

    [10] 徐克學(xué). 數(shù)量分類學(xué)[M]. 北京: 科學(xué)出版社,1994.

    [11] Govindaraju D R. Relationship between dispersal ability and levels of gene flow in plants[J]. Oikos, 1988, 52: 31-35.

    [12] Loveless M D, Hamrick J L. Ecological determinants of genetic structure in plant populations[J]. Ann. Rev. Ecol. Syst., 1984, 15(1): 65-95.

    [13] 臧德奎. 山東嶗山森林植物區(qū)系的研究[D].南京: 南京林業(yè)大學(xué),1990.

    [14] 汪小凡, 廖萬金, 宋志平. 小毛茛居群的遺傳分化及其與空間隔離的相關(guān)性[J]. 生物多樣性, 2001, 9(2): 138-144. 山 東 農(nóng) 業(yè) 科 學(xué) 2014,46(11):26~29

    2.3 居群遺傳結(jié)構(gòu)與遺傳一致度分析

    對(duì)鉤齒溲疏3個(gè)自然居群的遺傳結(jié)構(gòu)分析表明,總遺傳多樣性(Ht)為0.2190,其中居群內(nèi)(Hs)為0.1383,居群間(Dst)為0.0807,居群內(nèi)遺傳多樣度明顯高于居群間。居群間遺傳分化系數(shù)(Gst)為0.3685,即居群間遺傳變異僅占36.85%,因此居群內(nèi)遺傳分化遠(yuǎn)高于居群間。在居群遺傳學(xué)研究中,基因流被劃分為高、中、低三水平[11],而鉤齒溲疏各居群間基因流(Nm)大小為0.8566,介于0.25~0.99,屬于中等水平,表明居群間基因流相對(duì)有限。

    統(tǒng)計(jì)鉤齒溲疏各居群間的Neis遺傳距離GD和遺傳一致度GI以進(jìn)一步分析鉤齒溲疏居群的遺傳結(jié)構(gòu)。3個(gè)居群間遺傳距離GD變動(dòng)在0.1184~0.1803,平均為0.1514(表5),以嶗山與泰山的遺傳距離最近,嶗山與徂徠山的最遠(yuǎn)。3個(gè)居群間遺傳一致度GI變化范圍為0.8350~0.8884,平均為0.8598。泰山與嶗山居群的遺傳距離是3個(gè)居群間的最小值,遺傳一致度最高,但在地理位置上距離較遠(yuǎn);同樣,泰山與徂徠山地理距離最近,但遺傳一致度反比泰山與嶗山的低。這表明,鉤齒溲疏居群間的遺傳距離和地理位置并非直接相關(guān)。

    2.4 聚類分析

    UPGMA聚類結(jié)果(圖1)顯示,泰山居群14號(hào)和15號(hào)樣品遺傳相似度最高,為0.91。在相似系數(shù)為0.74時(shí)所有樣品可以分為3大類:第Ⅰ大類包含7個(gè)樣品,全部來自徂徠山;第Ⅱ大類包括來自嶗山的4個(gè)樣品和來自泰山的8個(gè)樣品;第Ⅲ大類只有2個(gè)樣品,均來自嶗山。第Ⅱ大類在相似系數(shù)為0.81處又可分為3個(gè)亞類,第1亞類含1個(gè)樣品,第2亞類含3個(gè)樣品,均來自嶗山;第3亞類含8個(gè)樣品,全部來自泰山。

    由聚類圖可看出,徂徠山的樣品單獨(dú)聚為一個(gè)大類,泰山的樣品聚在一個(gè)亞類里,只有嶗山的樣品分布在2個(gè)大類中,這反映出嶗山居群豐富的遺傳多樣性。同時(shí),徂徠山樣品單獨(dú)聚為第Ⅰ大類,泰山和嶗山的部分樣品聚為第Ⅱ大類,表明了嶗山與泰山的遺傳相似度較高。

    聚類結(jié)果顯示,多數(shù)情況下樣品來源決定樣品間的親緣關(guān)系,地理來源相同的樣品傾向于聚在一起,但供試樣品的聚類結(jié)果又并非嚴(yán)格按照地理來源來劃分。說明山東鉤齒溲疏的遺傳分化并非受地理來源單一因素的影響,是多個(gè)因素綜合作用的結(jié)果。

    3 討論

    在體現(xiàn)居群內(nèi)變異水平方面,多態(tài)位點(diǎn)百分率是重要指標(biāo)之一。在物種水平上鉤齒溲疏的多態(tài)位點(diǎn)百分率為68.44%,表明山東省內(nèi)鉤齒溲疏有較豐富的遺傳多樣性。

    Loveless等[12]曾統(tǒng)計(jì)過不同類型植物的遺傳變異水平和居群分化程度,得出了大范圍連續(xù)分布的異交植物的遺傳變異大部分存在于居群內(nèi)的結(jié)論;而鉤齒溲疏居群間遺傳分化系數(shù)(Gst)為0.3685,表明遺傳分化主要存在于居群內(nèi),同這一結(jié)論類似。鉤齒溲疏各居群間基因流(Nm)大小為0.8566,表明居群間基因流相對(duì)有限,是3個(gè)居群間存在一定程度的遺傳分化的原因之一。

    在居群水平上,徂徠山的遺傳多樣性最低,應(yīng)與徂徠山鉤齒溲疏生境類型變化相對(duì)有限有關(guān)。而嶗山居群的遺傳多樣性最為豐富,可能與嶗山的地理、氣候條件更為優(yōu)越和復(fù)雜有關(guān)。嶗山地處山東半島,東、南、北三面環(huán)海,西面與陸地相接。由于受海洋影響,氣候溫暖濕潤(rùn),南坡和北坡氣候迥異。嶗山植被類型豐富,群落組成上南坡?lián)碛休^多的亞熱帶區(qū)系成分,而北坡則以長(zhǎng)白區(qū)系成分為主[13]。嶗山復(fù)雜多樣的生境類型,高度的空間異質(zhì)性使得鉤齒溲疏不同基因型的植株都可以生長(zhǎng),所以遺傳多樣性處于較高水平。生態(tài)因子可作為控制居群中基因型頻率的一種選擇壓力,在相當(dāng)程度上對(duì)居群的遺傳分化起作用[14]。反過來,遺傳多樣性也是物種適應(yīng)環(huán)境變化的基礎(chǔ)。因此可以推測(cè),嶗山鉤齒溲疏居群較高的遺傳多樣性是對(duì)嶗山優(yōu)越、復(fù)雜的生態(tài)環(huán)境長(zhǎng)期適應(yīng)的結(jié)果。

    在調(diào)查中發(fā)現(xiàn),鉤齒溲疏嶗山居群豐富的遺傳多樣性不僅體現(xiàn)在分子水平上,也反映在形態(tài)上,尤其是花朵的大小,葉片的大小、形狀等。形態(tài)上豐富的變異為開發(fā)利用鉤齒溲疏野生資源、選育優(yōu)良類型提供了基礎(chǔ)。

    參 考 文 獻(xiàn):

    [1] 中國(guó)植物志編委會(huì). 中國(guó)植物志[M]. 北京: 科學(xué)出版社, 1982: 103-104.

    [2] Spiess E B. Genes in population[M]. 2nd ed. New York: John Wiley & Sons, 1989.

    [3] 潘瑩, 趙桂仿. 分子水平的遺傳多樣性及其測(cè)量方法[J]. 西北植物學(xué)報(bào), 1998, 18(4): 645-653.

    [4] Li G, Qurios C F. Sequence related amplified polymorphism(SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica[J]. Theoretical and Applied Genetics, 2001, 103(3): 455-461.

    [5] 張安世, 刑智峰, 劉永英, 等. SRAP分子標(biāo)記及其應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué), 2007, 35(9): 2562-2563.

    [6] 李莉, 彭建營(yíng), 白瑞霞, 等. SRAP與TRAP標(biāo)記及其在園藝植物研究中的應(yīng)用[J]. 西北植物學(xué)報(bào), 2006, 26(8): 1749-1752.

    [7] 張鮮艷, 張飛, 陳發(fā)棣, 等. 12 份不同地理居群野菊的遺傳多樣性分析[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 34(3): 48-54.

    [8] 林萍, 姚小華, 王開良, 等. 油茶長(zhǎng)林系列優(yōu)良無性系的SRAP分子鑒別及遺傳分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2010, 18(2): 272-279.

    [9] Kidwell K K, Osborn T C. Simple plant DNA isolation procedures[M]//Beckman J, Osborn T C. Plant genomes: methods for genetic and physical mapping. Dordrecht: Kluwer Academic Publishers, 1992: 1-13.

    [10] 徐克學(xué). 數(shù)量分類學(xué)[M]. 北京: 科學(xué)出版社,1994.

    [11] Govindaraju D R. Relationship between dispersal ability and levels of gene flow in plants[J]. Oikos, 1988, 52: 31-35.

    [12] Loveless M D, Hamrick J L. Ecological determinants of genetic structure in plant populations[J]. Ann. Rev. Ecol. Syst., 1984, 15(1): 65-95.

    [13] 臧德奎. 山東嶗山森林植物區(qū)系的研究[D].南京: 南京林業(yè)大學(xué),1990.

    [14] 汪小凡, 廖萬金, 宋志平. 小毛茛居群的遺傳分化及其與空間隔離的相關(guān)性[J]. 生物多樣性, 2001, 9(2): 138-144. 山 東 農(nóng) 業(yè) 科 學(xué) 2014,46(11):26~29

    2.3 居群遺傳結(jié)構(gòu)與遺傳一致度分析

    對(duì)鉤齒溲疏3個(gè)自然居群的遺傳結(jié)構(gòu)分析表明,總遺傳多樣性(Ht)為0.2190,其中居群內(nèi)(Hs)為0.1383,居群間(Dst)為0.0807,居群內(nèi)遺傳多樣度明顯高于居群間。居群間遺傳分化系數(shù)(Gst)為0.3685,即居群間遺傳變異僅占36.85%,因此居群內(nèi)遺傳分化遠(yuǎn)高于居群間。在居群遺傳學(xué)研究中,基因流被劃分為高、中、低三水平[11],而鉤齒溲疏各居群間基因流(Nm)大小為0.8566,介于0.25~0.99,屬于中等水平,表明居群間基因流相對(duì)有限。

    統(tǒng)計(jì)鉤齒溲疏各居群間的Neis遺傳距離GD和遺傳一致度GI以進(jìn)一步分析鉤齒溲疏居群的遺傳結(jié)構(gòu)。3個(gè)居群間遺傳距離GD變動(dòng)在0.1184~0.1803,平均為0.1514(表5),以嶗山與泰山的遺傳距離最近,嶗山與徂徠山的最遠(yuǎn)。3個(gè)居群間遺傳一致度GI變化范圍為0.8350~0.8884,平均為0.8598。泰山與嶗山居群的遺傳距離是3個(gè)居群間的最小值,遺傳一致度最高,但在地理位置上距離較遠(yuǎn);同樣,泰山與徂徠山地理距離最近,但遺傳一致度反比泰山與嶗山的低。這表明,鉤齒溲疏居群間的遺傳距離和地理位置并非直接相關(guān)。

    2.4 聚類分析

    UPGMA聚類結(jié)果(圖1)顯示,泰山居群14號(hào)和15號(hào)樣品遺傳相似度最高,為0.91。在相似系數(shù)為0.74時(shí)所有樣品可以分為3大類:第Ⅰ大類包含7個(gè)樣品,全部來自徂徠山;第Ⅱ大類包括來自嶗山的4個(gè)樣品和來自泰山的8個(gè)樣品;第Ⅲ大類只有2個(gè)樣品,均來自嶗山。第Ⅱ大類在相似系數(shù)為0.81處又可分為3個(gè)亞類,第1亞類含1個(gè)樣品,第2亞類含3個(gè)樣品,均來自嶗山;第3亞類含8個(gè)樣品,全部來自泰山。

    由聚類圖可看出,徂徠山的樣品單獨(dú)聚為一個(gè)大類,泰山的樣品聚在一個(gè)亞類里,只有嶗山的樣品分布在2個(gè)大類中,這反映出嶗山居群豐富的遺傳多樣性。同時(shí),徂徠山樣品單獨(dú)聚為第Ⅰ大類,泰山和嶗山的部分樣品聚為第Ⅱ大類,表明了嶗山與泰山的遺傳相似度較高。

    聚類結(jié)果顯示,多數(shù)情況下樣品來源決定樣品間的親緣關(guān)系,地理來源相同的樣品傾向于聚在一起,但供試樣品的聚類結(jié)果又并非嚴(yán)格按照地理來源來劃分。說明山東鉤齒溲疏的遺傳分化并非受地理來源單一因素的影響,是多個(gè)因素綜合作用的結(jié)果。

    3 討論

    在體現(xiàn)居群內(nèi)變異水平方面,多態(tài)位點(diǎn)百分率是重要指標(biāo)之一。在物種水平上鉤齒溲疏的多態(tài)位點(diǎn)百分率為68.44%,表明山東省內(nèi)鉤齒溲疏有較豐富的遺傳多樣性。

    Loveless等[12]曾統(tǒng)計(jì)過不同類型植物的遺傳變異水平和居群分化程度,得出了大范圍連續(xù)分布的異交植物的遺傳變異大部分存在于居群內(nèi)的結(jié)論;而鉤齒溲疏居群間遺傳分化系數(shù)(Gst)為0.3685,表明遺傳分化主要存在于居群內(nèi),同這一結(jié)論類似。鉤齒溲疏各居群間基因流(Nm)大小為0.8566,表明居群間基因流相對(duì)有限,是3個(gè)居群間存在一定程度的遺傳分化的原因之一。

    在居群水平上,徂徠山的遺傳多樣性最低,應(yīng)與徂徠山鉤齒溲疏生境類型變化相對(duì)有限有關(guān)。而嶗山居群的遺傳多樣性最為豐富,可能與嶗山的地理、氣候條件更為優(yōu)越和復(fù)雜有關(guān)。嶗山地處山東半島,東、南、北三面環(huán)海,西面與陸地相接。由于受海洋影響,氣候溫暖濕潤(rùn),南坡和北坡氣候迥異。嶗山植被類型豐富,群落組成上南坡?lián)碛休^多的亞熱帶區(qū)系成分,而北坡則以長(zhǎng)白區(qū)系成分為主[13]。嶗山復(fù)雜多樣的生境類型,高度的空間異質(zhì)性使得鉤齒溲疏不同基因型的植株都可以生長(zhǎng),所以遺傳多樣性處于較高水平。生態(tài)因子可作為控制居群中基因型頻率的一種選擇壓力,在相當(dāng)程度上對(duì)居群的遺傳分化起作用[14]。反過來,遺傳多樣性也是物種適應(yīng)環(huán)境變化的基礎(chǔ)。因此可以推測(cè),嶗山鉤齒溲疏居群較高的遺傳多樣性是對(duì)嶗山優(yōu)越、復(fù)雜的生態(tài)環(huán)境長(zhǎng)期適應(yīng)的結(jié)果。

    在調(diào)查中發(fā)現(xiàn),鉤齒溲疏嶗山居群豐富的遺傳多樣性不僅體現(xiàn)在分子水平上,也反映在形態(tài)上,尤其是花朵的大小,葉片的大小、形狀等。形態(tài)上豐富的變異為開發(fā)利用鉤齒溲疏野生資源、選育優(yōu)良類型提供了基礎(chǔ)。

    參 考 文 獻(xiàn):

    [1] 中國(guó)植物志編委會(huì). 中國(guó)植物志[M]. 北京: 科學(xué)出版社, 1982: 103-104.

    [2] Spiess E B. Genes in population[M]. 2nd ed. New York: John Wiley & Sons, 1989.

    [3] 潘瑩, 趙桂仿. 分子水平的遺傳多樣性及其測(cè)量方法[J]. 西北植物學(xué)報(bào), 1998, 18(4): 645-653.

    [4] Li G, Qurios C F. Sequence related amplified polymorphism(SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica[J]. Theoretical and Applied Genetics, 2001, 103(3): 455-461.

    [5] 張安世, 刑智峰, 劉永英, 等. SRAP分子標(biāo)記及其應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué), 2007, 35(9): 2562-2563.

    [6] 李莉, 彭建營(yíng), 白瑞霞, 等. SRAP與TRAP標(biāo)記及其在園藝植物研究中的應(yīng)用[J]. 西北植物學(xué)報(bào), 2006, 26(8): 1749-1752.

    [7] 張鮮艷, 張飛, 陳發(fā)棣, 等. 12 份不同地理居群野菊的遺傳多樣性分析[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 34(3): 48-54.

    [8] 林萍, 姚小華, 王開良, 等. 油茶長(zhǎng)林系列優(yōu)良無性系的SRAP分子鑒別及遺傳分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2010, 18(2): 272-279.

    [9] Kidwell K K, Osborn T C. Simple plant DNA isolation procedures[M]//Beckman J, Osborn T C. Plant genomes: methods for genetic and physical mapping. Dordrecht: Kluwer Academic Publishers, 1992: 1-13.

    [10] 徐克學(xué). 數(shù)量分類學(xué)[M]. 北京: 科學(xué)出版社,1994.

    [11] Govindaraju D R. Relationship between dispersal ability and levels of gene flow in plants[J]. Oikos, 1988, 52: 31-35.

    [12] Loveless M D, Hamrick J L. Ecological determinants of genetic structure in plant populations[J]. Ann. Rev. Ecol. Syst., 1984, 15(1): 65-95.

    [13] 臧德奎. 山東嶗山森林植物區(qū)系的研究[D].南京: 南京林業(yè)大學(xué),1990.

    [14] 汪小凡, 廖萬金, 宋志平. 小毛茛居群的遺傳分化及其與空間隔離的相關(guān)性[J]. 生物多樣性, 2001, 9(2): 138-144. 山 東 農(nóng) 業(yè) 科 學(xué) 2014,46(11):26~29

    猜你喜歡
    遺傳多樣性
    從葉綠體DNA角度分析云南省砂梨地方品種遺傳多樣性
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    金魚起源及遺傳多樣性研究進(jìn)展
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    金銀花SSR指紋圖譜的構(gòu)建及遺傳多樣性分析
    色视频www国产| 在线观看免费视频日本深夜| 熟女人妻精品中文字幕| 男女之事视频高清在线观看| 亚洲七黄色美女视频| 国产黄a三级三级三级人| 波多野结衣高清作品| 亚洲国产精品专区欧美| 国产精品福利在线免费观看| 免费观看a级毛片全部| 日本色播在线视频| 亚州av有码| 色哟哟·www| 国产成年人精品一区二区| 日日撸夜夜添| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久大av| 美女黄网站色视频| 国产亚洲91精品色在线| 亚洲熟妇中文字幕五十中出| 一级a做视频免费观看| 我要看日韩黄色一级片| 舔av片在线| 最近视频中文字幕2019在线8| 亚洲精品国产成人久久av| 男女国产视频网站| 2018国产大陆天天弄谢| 亚洲av.av天堂| 日本wwww免费看| 免费av毛片视频| 日本午夜av视频| av在线观看视频网站免费| 一边亲一边摸免费视频| 五月天丁香电影| 亚洲不卡免费看| 精品少妇黑人巨大在线播放| 日韩精品有码人妻一区| 精品久久久精品久久久| 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 免费av毛片视频| 好男人视频免费观看在线| 高清毛片免费看| 久久久久久久久久久免费av| 亚洲av男天堂| 亚洲四区av| 国产乱人偷精品视频| 亚洲av二区三区四区| 日韩精品有码人妻一区| 亚洲欧美清纯卡通| 成人美女网站在线观看视频| 久久精品人妻少妇| 超碰av人人做人人爽久久| 国产乱来视频区| kizo精华| 日韩中字成人| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av天美| 国产精品.久久久| 亚洲精品成人久久久久久| av在线天堂中文字幕| 久久精品国产亚洲网站| 伊人久久精品亚洲午夜| 日韩精品有码人妻一区| 亚洲最大成人手机在线| 成人午夜精彩视频在线观看| 国产精品无大码| 亚洲av一区综合| 国产69精品久久久久777片| 尾随美女入室| 日韩精品有码人妻一区| 99久国产av精品国产电影| 真实男女啪啪啪动态图| 国产爱豆传媒在线观看| 中国国产av一级| 建设人人有责人人尽责人人享有的 | 看十八女毛片水多多多| 成人国产麻豆网| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 天堂av国产一区二区熟女人妻| 内射极品少妇av片p| 白带黄色成豆腐渣| 亚洲成色77777| 日本黄大片高清| 男人舔奶头视频| 黄片wwwwww| 少妇丰满av| 99久久中文字幕三级久久日本| eeuss影院久久| 26uuu在线亚洲综合色| 热99在线观看视频| 免费电影在线观看免费观看| 免费av毛片视频| 国产一区亚洲一区在线观看| 五月玫瑰六月丁香| 国产高清国产精品国产三级 | 自拍偷自拍亚洲精品老妇| 美女xxoo啪啪120秒动态图| 亚洲精品中文字幕在线视频 | 最近视频中文字幕2019在线8| 麻豆av噜噜一区二区三区| 中国美白少妇内射xxxbb| 欧美日韩在线观看h| 美女xxoo啪啪120秒动态图| 如何舔出高潮| 欧美三级亚洲精品| 国产精品爽爽va在线观看网站| 麻豆精品久久久久久蜜桃| 伊人久久精品亚洲午夜| 日韩三级伦理在线观看| 久久99热6这里只有精品| 少妇熟女aⅴ在线视频| 久久99热这里只有精品18| 日韩欧美精品免费久久| 免费少妇av软件| 日韩人妻高清精品专区| 两个人视频免费观看高清| 九九爱精品视频在线观看| 免费观看av网站的网址| 亚洲av电影在线观看一区二区三区 | 成人亚洲精品av一区二区| 看十八女毛片水多多多| 国产精品三级大全| 国产91av在线免费观看| 日韩欧美 国产精品| or卡值多少钱| 国产久久久一区二区三区| 精品熟女少妇av免费看| 少妇裸体淫交视频免费看高清| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 中文字幕免费在线视频6| 全区人妻精品视频| 国产精品三级大全| 99久国产av精品| 黄色日韩在线| 综合色av麻豆| 久久精品国产亚洲av天美| 亚洲电影在线观看av| 床上黄色一级片| 国产精品国产三级国产专区5o| 久久精品国产亚洲av涩爱| 搡女人真爽免费视频火全软件| 亚洲va在线va天堂va国产| 男女国产视频网站| 国国产精品蜜臀av免费| 亚洲精品国产成人久久av| 日韩av在线大香蕉| 午夜福利视频1000在线观看| 国产美女午夜福利| 亚洲图色成人| 啦啦啦韩国在线观看视频| 尾随美女入室| .国产精品久久| 中文字幕av成人在线电影| 一级毛片 在线播放| 欧美zozozo另类| 晚上一个人看的免费电影| 亚洲精品影视一区二区三区av| 久久久久久国产a免费观看| 亚洲欧美日韩卡通动漫| 成人性生交大片免费视频hd| 亚洲人与动物交配视频| 简卡轻食公司| 免费少妇av软件| 老司机影院毛片| 日韩三级伦理在线观看| 国内精品宾馆在线| 免费观看精品视频网站| 国产一区二区在线观看日韩| 成人午夜高清在线视频| 久久久久久久大尺度免费视频| 人体艺术视频欧美日本| 欧美性感艳星| 亚洲av成人精品一区久久| 亚洲精品自拍成人| 婷婷色综合www| 午夜视频国产福利| 久久久成人免费电影| 国产人妻一区二区三区在| 伦精品一区二区三区| 99久久中文字幕三级久久日本| 午夜激情福利司机影院| 国产有黄有色有爽视频| 亚洲av一区综合| 在线观看人妻少妇| 国产真实伦视频高清在线观看| 能在线免费看毛片的网站| 秋霞伦理黄片| 国产精品久久久久久精品电影| 日韩 亚洲 欧美在线| www.av在线官网国产| 亚洲av福利一区| 又爽又黄无遮挡网站| 精品国产露脸久久av麻豆 | 啦啦啦啦在线视频资源| 日本-黄色视频高清免费观看| 一级二级三级毛片免费看| 国产在线一区二区三区精| 欧美精品一区二区大全| av国产久精品久网站免费入址| 亚洲丝袜综合中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 国产人妻一区二区三区在| 久久久久国产网址| 国产视频首页在线观看| 国产精品麻豆人妻色哟哟久久 | 麻豆国产97在线/欧美| 日韩欧美 国产精品| 亚洲成人中文字幕在线播放| 激情 狠狠 欧美| 亚洲美女搞黄在线观看| 午夜日本视频在线| 在现免费观看毛片| 波野结衣二区三区在线| 视频中文字幕在线观看| 人人妻人人看人人澡| 国产成人91sexporn| 人妻系列 视频| 高清午夜精品一区二区三区| 一个人观看的视频www高清免费观看| 日韩欧美一区视频在线观看 | 亚洲一区高清亚洲精品| av网站免费在线观看视频 | 人妻一区二区av| 国产精品一及| 成年av动漫网址| 男的添女的下面高潮视频| 少妇猛男粗大的猛烈进出视频 | 男人舔奶头视频| 丰满少妇做爰视频| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 午夜精品国产一区二区电影 | 免费播放大片免费观看视频在线观看| 国产亚洲一区二区精品| 1000部很黄的大片| 好男人视频免费观看在线| 久久亚洲国产成人精品v| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 久久久亚洲精品成人影院| 欧美成人a在线观看| 国产亚洲精品av在线| 日韩三级伦理在线观看| 精品久久久久久久久久久久久| 亚洲精品,欧美精品| 午夜免费激情av| 亚洲国产日韩欧美精品在线观看| 少妇的逼好多水| 99久久人妻综合| 2022亚洲国产成人精品| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载 | 中文乱码字字幕精品一区二区三区 | 国产成人精品婷婷| 高清在线视频一区二区三区| 久久国产乱子免费精品| 非洲黑人性xxxx精品又粗又长| 一级毛片黄色毛片免费观看视频| 少妇猛男粗大的猛烈进出视频 | 国产精品无大码| 国产精品久久久久久精品电影小说 | 免费黄网站久久成人精品| 2022亚洲国产成人精品| 国产黄色视频一区二区在线观看| 两个人的视频大全免费| 美女cb高潮喷水在线观看| 十八禁网站网址无遮挡 | 国产毛片a区久久久久| av一本久久久久| 国产色爽女视频免费观看| 五月伊人婷婷丁香| 国产黄色视频一区二区在线观看| 亚洲第一区二区三区不卡| 啦啦啦中文免费视频观看日本| 精品久久久久久久人妻蜜臀av| 亚洲av男天堂| 内地一区二区视频在线| 男人舔女人下体高潮全视频| 免费看日本二区| ponron亚洲| 91久久精品国产一区二区成人| 免费av观看视频| 国产成人精品婷婷| 黄色欧美视频在线观看| av又黄又爽大尺度在线免费看| 成人亚洲精品av一区二区| 国产黄频视频在线观看| 国产中年淑女户外野战色| 国产成人精品久久久久久| 国精品久久久久久国模美| 午夜免费男女啪啪视频观看| 亚洲熟妇中文字幕五十中出| av黄色大香蕉| 亚洲综合色惰| 极品少妇高潮喷水抽搐| 婷婷六月久久综合丁香| 日本三级黄在线观看| 哪个播放器可以免费观看大片| 国产中年淑女户外野战色| 国产精品一区二区性色av| 美女内射精品一级片tv| 欧美不卡视频在线免费观看| 91久久精品国产一区二区成人| 97超碰精品成人国产| 成人国产麻豆网| 日韩 亚洲 欧美在线| 久久精品熟女亚洲av麻豆精品 | 高清毛片免费看| 女人久久www免费人成看片| 爱豆传媒免费全集在线观看| 国产综合精华液| av黄色大香蕉| av在线蜜桃| 国产高清有码在线观看视频| 国产精品一二三区在线看| 18+在线观看网站| 男人舔女人下体高潮全视频| 久久久久网色| 欧美不卡视频在线免费观看| 国产精品一区二区三区四区久久| 人妻少妇偷人精品九色| 看免费成人av毛片| 天堂√8在线中文| 国产视频内射| 老司机影院毛片| 淫秽高清视频在线观看| 欧美成人精品欧美一级黄| 大陆偷拍与自拍| 一区二区三区四区激情视频| 精品久久久久久久人妻蜜臀av| 日本午夜av视频| 亚洲国产色片| 白带黄色成豆腐渣| 久久久久久久午夜电影| 亚洲av在线观看美女高潮| 22中文网久久字幕| 亚洲精品影视一区二区三区av| 国产成人午夜福利电影在线观看| 九九久久精品国产亚洲av麻豆| 好男人在线观看高清免费视频| 精品少妇黑人巨大在线播放| 2018国产大陆天天弄谢| 午夜久久久久精精品| 黄片wwwwww| 春色校园在线视频观看| 少妇的逼水好多| 精品久久久久久电影网| 97精品久久久久久久久久精品| 国产午夜福利久久久久久| 18禁在线无遮挡免费观看视频| 人妻一区二区av| 亚洲国产欧美人成| 女的被弄到高潮叫床怎么办| 国产精品一区二区性色av| 一区二区三区四区激情视频| 赤兔流量卡办理| 国产av不卡久久| av黄色大香蕉| 最新中文字幕久久久久| 超碰97精品在线观看| 免费观看av网站的网址| 天堂影院成人在线观看| 免费av不卡在线播放| 色哟哟·www| 极品少妇高潮喷水抽搐| 国产不卡一卡二| 国产乱来视频区| 1000部很黄的大片| 免费观看的影片在线观看| 国产探花在线观看一区二区| 午夜精品一区二区三区免费看| 欧美bdsm另类| 美女高潮的动态| 午夜福利在线观看免费完整高清在| 69av精品久久久久久| 高清视频免费观看一区二区 | 亚洲va在线va天堂va国产| 成人高潮视频无遮挡免费网站| 亚洲av电影在线观看一区二区三区 | 免费黄色在线免费观看| 成人一区二区视频在线观看| 美女黄网站色视频| 蜜桃亚洲精品一区二区三区| 大陆偷拍与自拍| 亚洲美女视频黄频| 一夜夜www| 亚洲性久久影院| 麻豆国产97在线/欧美| av黄色大香蕉| 草草在线视频免费看| 麻豆精品久久久久久蜜桃| 国产高潮美女av| 亚洲aⅴ乱码一区二区在线播放| 成年女人在线观看亚洲视频 | 亚洲综合色惰| 狠狠精品人妻久久久久久综合| 人人妻人人澡欧美一区二区| 亚洲综合精品二区| 看十八女毛片水多多多| 亚洲欧美日韩卡通动漫| 色5月婷婷丁香| 国产精品久久久久久精品电影| 亚洲最大成人手机在线| 美女脱内裤让男人舔精品视频| 一级毛片久久久久久久久女| 亚洲精华国产精华液的使用体验| 久久人人爽人人爽人人片va| 狠狠精品人妻久久久久久综合| 男人舔奶头视频| 水蜜桃什么品种好| 边亲边吃奶的免费视频| 一本一本综合久久| 人体艺术视频欧美日本| 亚洲精品久久午夜乱码| 欧美日本视频| 男女国产视频网站| 成年女人在线观看亚洲视频 | 男女下面进入的视频免费午夜| 亚洲国产色片| 欧美日本视频| 亚洲欧美精品专区久久| 女人十人毛片免费观看3o分钟| 亚洲av在线观看美女高潮| ponron亚洲| 丰满乱子伦码专区| 欧美极品一区二区三区四区| 亚洲电影在线观看av| 综合色av麻豆| 蜜桃久久精品国产亚洲av| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 国产v大片淫在线免费观看| 精品国产一区二区三区久久久樱花 | 精品久久久久久成人av| 国产精品国产三级国产av玫瑰| 成人综合一区亚洲| 搡老乐熟女国产| 不卡视频在线观看欧美| 国产精品久久视频播放| 国产高潮美女av| 青春草国产在线视频| 日韩中字成人| 日本一本二区三区精品| 亚洲最大成人av| 看免费成人av毛片| 国产亚洲av嫩草精品影院| 国产乱人视频| 18禁动态无遮挡网站| 国产成人免费观看mmmm| 亚洲av中文字字幕乱码综合| 欧美性猛交╳xxx乱大交人| www.av在线官网国产| 18禁裸乳无遮挡免费网站照片| 我的女老师完整版在线观看| 午夜久久久久精精品| 九色成人免费人妻av| 亚洲无线观看免费| 男女视频在线观看网站免费| 欧美xxⅹ黑人| 99热全是精品| 亚洲av成人精品一二三区| 国产日韩欧美在线精品| 中国国产av一级| 女人久久www免费人成看片| 99久久精品国产国产毛片| 亚洲精品亚洲一区二区| 精品99又大又爽又粗少妇毛片| 免费不卡的大黄色大毛片视频在线观看 | 婷婷六月久久综合丁香| 久久精品熟女亚洲av麻豆精品 | 国产伦精品一区二区三区四那| 欧美日韩视频高清一区二区三区二| 亚洲经典国产精华液单| 中文乱码字字幕精品一区二区三区 | 国产亚洲精品久久久com| 亚洲av中文av极速乱| 欧美极品一区二区三区四区| 亚洲av福利一区| av在线亚洲专区| 一个人观看的视频www高清免费观看| 午夜福利在线观看免费完整高清在| 在线免费观看不下载黄p国产| 黄片wwwwww| 亚洲人成网站高清观看| 激情 狠狠 欧美| 国产白丝娇喘喷水9色精品| 禁无遮挡网站| av播播在线观看一区| 日本三级黄在线观看| 少妇的逼好多水| 一个人看视频在线观看www免费| 国模一区二区三区四区视频| 91午夜精品亚洲一区二区三区| 欧美日韩亚洲高清精品| 精品午夜福利在线看| av在线亚洲专区| 日韩欧美精品免费久久| 日产精品乱码卡一卡2卡三| 亚洲不卡免费看| 日本猛色少妇xxxxx猛交久久| 国产伦在线观看视频一区| av国产久精品久网站免费入址| 精品人妻偷拍中文字幕| 免费黄网站久久成人精品| 国产国拍精品亚洲av在线观看| 搞女人的毛片| 99久国产av精品国产电影| 久久99蜜桃精品久久| 亚洲天堂国产精品一区在线| 我要看日韩黄色一级片| 国产伦在线观看视频一区| 超碰av人人做人人爽久久| 欧美激情久久久久久爽电影| 亚洲精品久久久久久婷婷小说| 热99在线观看视频| 成人亚洲精品av一区二区| eeuss影院久久| 亚洲精品国产av蜜桃| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 日本午夜av视频| 免费黄网站久久成人精品| 日韩制服骚丝袜av| 色尼玛亚洲综合影院| av免费在线看不卡| 天堂网av新在线| 国产av不卡久久| 久久久精品免费免费高清| 国产高清国产精品国产三级 | 免费观看性生交大片5| 亚洲欧美日韩东京热| 中国国产av一级| 亚洲av成人精品一区久久| 亚洲综合色惰| 国产一区有黄有色的免费视频 | 狂野欧美激情性xxxx在线观看| 亚洲精华国产精华液的使用体验| 欧美性感艳星| 色综合色国产| 国产人妻一区二区三区在| 久久久久精品久久久久真实原创| 天堂av国产一区二区熟女人妻| 91精品伊人久久大香线蕉| 亚洲激情五月婷婷啪啪| 极品少妇高潮喷水抽搐| 青青草视频在线视频观看| 高清午夜精品一区二区三区| 亚洲成人一二三区av| 亚洲熟妇中文字幕五十中出| 九草在线视频观看| 国产麻豆成人av免费视频| 国内揄拍国产精品人妻在线| 天堂av国产一区二区熟女人妻| 午夜福利视频精品| 日韩av不卡免费在线播放| 久久久久久久久久久丰满| 乱人视频在线观看| 国产精品久久久久久久电影| 久久草成人影院| 一边亲一边摸免费视频| 亚洲成色77777| 日韩一区二区三区影片| 日韩欧美三级三区| 欧美人与善性xxx| 成年人午夜在线观看视频 | 免费看av在线观看网站| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩东京热| 国产精品一区www在线观看| 热99在线观看视频| 亚洲国产日韩欧美精品在线观看| 18禁裸乳无遮挡免费网站照片| 深夜a级毛片| 可以在线观看毛片的网站| 亚洲高清免费不卡视频| 2021天堂中文幕一二区在线观| 天堂网av新在线| 卡戴珊不雅视频在线播放| 日日啪夜夜撸| 观看免费一级毛片| 亚洲av一区综合| 边亲边吃奶的免费视频| 久久久久久久久久久丰满| 超碰av人人做人人爽久久| 边亲边吃奶的免费视频| 亚洲国产精品国产精品| 听说在线观看完整版免费高清| 国内精品美女久久久久久| 狠狠精品人妻久久久久久综合| 国产高清国产精品国产三级 | 国产精品一区二区在线观看99 | 三级国产精品欧美在线观看| 亚洲av不卡在线观看| 男女国产视频网站| 熟女电影av网| 色视频www国产| 国产成人福利小说| 最近中文字幕2019免费版| 女人被狂操c到高潮| 国产国拍精品亚洲av在线观看| 老司机影院毛片| 国产一级毛片在线| 一级毛片电影观看| 亚洲电影在线观看av| 色视频www国产| 国产欧美另类精品又又久久亚洲欧美| 99久久九九国产精品国产免费| 建设人人有责人人尽责人人享有的 | 人人妻人人看人人澡| 国内精品宾馆在线| 亚洲国产av新网站|