• <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)建及遺傳多樣性分析
    老司机在亚洲福利影院| 国产免费福利视频在线观看| 黄片小视频在线播放| 国产在线观看jvid| 12—13女人毛片做爰片一| 久久久精品区二区三区| 999精品在线视频| 老熟妇仑乱视频hdxx| 99久久人妻综合| 国产日韩欧美在线精品| 青春草视频在线免费观看| 99国产精品免费福利视频| 亚洲第一青青草原| 两个人免费观看高清视频| 午夜福利乱码中文字幕| 亚洲第一青青草原| 久久性视频一级片| 两个人看的免费小视频| 亚洲欧美精品综合一区二区三区| av欧美777| 老司机影院毛片| 伊人亚洲综合成人网| 麻豆国产av国片精品| a在线观看视频网站| 老司机午夜十八禁免费视频| 精品少妇内射三级| www.自偷自拍.com| 久久国产精品大桥未久av| tocl精华| 一区二区三区四区激情视频| 少妇粗大呻吟视频| 日本黄色日本黄色录像| 欧美精品一区二区免费开放| 久久久久久久久免费视频了| av免费在线观看网站| 国产一区有黄有色的免费视频| 美女高潮到喷水免费观看| 国产精品成人在线| 人人澡人人妻人| 一二三四在线观看免费中文在| 丁香六月天网| 亚洲人成77777在线视频| 日本av免费视频播放| www.熟女人妻精品国产| 欧美成人午夜精品| 高清视频免费观看一区二区| 男女边摸边吃奶| 少妇的丰满在线观看| 9191精品国产免费久久| 亚洲 国产 在线| 一二三四社区在线视频社区8| 亚洲精品日韩在线中文字幕| 中国国产av一级| 一边摸一边抽搐一进一出视频| 嫩草影视91久久| 一二三四社区在线视频社区8| 国产精品一区二区在线不卡| 大香蕉久久成人网| 国产男人的电影天堂91| 在线观看免费高清a一片| 免费在线观看视频国产中文字幕亚洲 | 精品人妻熟女毛片av久久网站| 在线精品无人区一区二区三| 国产精品九九99| 成人亚洲精品一区在线观看| 中国美女看黄片| 最新在线观看一区二区三区| 国产精品一区二区在线不卡| 婷婷色av中文字幕| 大香蕉久久网| 成人国语在线视频| 中文字幕精品免费在线观看视频| 国产真人三级小视频在线观看| 国产有黄有色有爽视频| 亚洲综合色网址| 超碰成人久久| 亚洲视频免费观看视频| 久9热在线精品视频| 国产成人av激情在线播放| 国产视频一区二区在线看| 久久精品国产亚洲av香蕉五月 | 成人免费观看视频高清| 99久久精品国产亚洲精品| 丝袜美足系列| 亚洲欧美色中文字幕在线| 久久中文字幕一级| 亚洲avbb在线观看| 国产精品国产三级国产专区5o| a级片在线免费高清观看视频| 视频在线观看一区二区三区| 看免费av毛片| 精品亚洲成a人片在线观看| 国产成人精品在线电影| 国产男人的电影天堂91| 建设人人有责人人尽责人人享有的| 欧美午夜高清在线| 国产熟女午夜一区二区三区| 黄色毛片三级朝国网站| videos熟女内射| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 又黄又粗又硬又大视频| 两个人免费观看高清视频| 男女边摸边吃奶| 欧美大码av| 国产av又大| 1024视频免费在线观看| 免费观看a级毛片全部| 午夜激情久久久久久久| 看免费av毛片| a 毛片基地| 国产欧美日韩一区二区三区在线| 少妇精品久久久久久久| 又大又爽又粗| 亚洲精品av麻豆狂野| 在线亚洲精品国产二区图片欧美| 最近中文字幕2019免费版| 亚洲av男天堂| 一级黄色大片毛片| 在线天堂中文资源库| 欧美日韩一级在线毛片| 中文字幕高清在线视频| 精品一区二区三卡| 中国美女看黄片| 国产一区二区 视频在线| 丰满饥渴人妻一区二区三| 在线精品无人区一区二区三| 日韩大片免费观看网站| 亚洲欧美精品综合一区二区三区| 国产亚洲欧美在线一区二区| 中文字幕色久视频| 在线十欧美十亚洲十日本专区| 在线观看免费视频网站a站| 精品福利永久在线观看| 亚洲 欧美一区二区三区| 18禁观看日本| 国产成人精品无人区| 伦理电影免费视频| 韩国精品一区二区三区| 亚洲精华国产精华精| 国产高清国产精品国产三级| 性色av一级| 久久毛片免费看一区二区三区| 欧美精品人与动牲交sv欧美| 国产福利在线免费观看视频| 母亲3免费完整高清在线观看| 老司机午夜福利在线观看视频 | 永久免费av网站大全| 麻豆国产av国片精品| 国产亚洲欧美精品永久| 丁香六月欧美| 99热全是精品| 国产精品麻豆人妻色哟哟久久| 1024香蕉在线观看| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频 | 国产免费av片在线观看野外av| 国产亚洲精品久久久久5区| 日本一区二区免费在线视频| 国产高清videossex| 动漫黄色视频在线观看| 天堂中文最新版在线下载| 大码成人一级视频| 另类精品久久| 欧美午夜高清在线| av网站免费在线观看视频| 日本撒尿小便嘘嘘汇集6| av天堂久久9| 欧美激情久久久久久爽电影 | 女性被躁到高潮视频| 亚洲国产毛片av蜜桃av| 香蕉丝袜av| 国产区一区二久久| 一区二区三区乱码不卡18| 久久狼人影院| bbb黄色大片| 99国产精品一区二区三区| 国产又爽黄色视频| 大陆偷拍与自拍| a级毛片在线看网站| 欧美成狂野欧美在线观看| 日韩大片免费观看网站| 国产精品一区二区在线不卡| 9热在线视频观看99| 久久久久久久国产电影| 热99久久久久精品小说推荐| 黄频高清免费视频| 国产麻豆69| 午夜福利视频精品| 久久久久久久大尺度免费视频| 欧美午夜高清在线| 亚洲精华国产精华精| 自拍欧美九色日韩亚洲蝌蚪91| 精品卡一卡二卡四卡免费| 丝袜在线中文字幕| 精品乱码久久久久久99久播| 汤姆久久久久久久影院中文字幕| 亚洲avbb在线观看| 欧美乱码精品一区二区三区| 男人舔女人的私密视频| 日本欧美视频一区| 久久久国产精品麻豆| 精品福利永久在线观看| 精品熟女少妇八av免费久了| 国产片内射在线| 高清av免费在线| 国产片内射在线| 久久国产亚洲av麻豆专区| 岛国在线观看网站| 精品亚洲乱码少妇综合久久| 精品乱码久久久久久99久播| 这个男人来自地球电影免费观看| av又黄又爽大尺度在线免费看| 久久99热这里只频精品6学生| 九色亚洲精品在线播放| 一区福利在线观看| 少妇裸体淫交视频免费看高清 | 欧美成人午夜精品| 国产亚洲精品久久久久5区| 国产日韩一区二区三区精品不卡| 国产亚洲精品第一综合不卡| av网站免费在线观看视频| av国产精品久久久久影院| 国产精品麻豆人妻色哟哟久久| 久久精品亚洲av国产电影网| 国产老妇伦熟女老妇高清| 男女床上黄色一级片免费看| 一本—道久久a久久精品蜜桃钙片| 欧美成人午夜精品| 国产成人精品无人区| 久久青草综合色| 无遮挡黄片免费观看| 亚洲av国产av综合av卡| 日日摸夜夜添夜夜添小说| 国产真人三级小视频在线观看| 欧美日韩视频精品一区| 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 中国国产av一级| 波多野结衣av一区二区av| 欧美精品一区二区免费开放| 国产精品99久久99久久久不卡| 日韩三级视频一区二区三区| 久久国产精品人妻蜜桃| h视频一区二区三区| 中文字幕人妻熟女乱码| 老司机在亚洲福利影院| 国产亚洲精品久久久久5区| 我的亚洲天堂| 国产免费av片在线观看野外av| 极品人妻少妇av视频| 亚洲欧美日韩另类电影网站| 久久人妻熟女aⅴ| 亚洲欧美日韩高清在线视频 | 狠狠狠狠99中文字幕| 国产成人av激情在线播放| 国产xxxxx性猛交| 别揉我奶头~嗯~啊~动态视频 | 国产欧美日韩一区二区精品| 亚洲欧美精品综合一区二区三区| 男女之事视频高清在线观看| 精品高清国产在线一区| 久久久久精品人妻al黑| 丰满迷人的少妇在线观看| 午夜福利视频在线观看免费| 国产成人一区二区三区免费视频网站| 久久久精品区二区三区| 在线看a的网站| 国产成人a∨麻豆精品| av视频免费观看在线观看| 久久这里只有精品19| 丝袜人妻中文字幕| av不卡在线播放| 国产一区二区激情短视频 | 亚洲精品粉嫩美女一区| 久久综合国产亚洲精品| 国产欧美日韩一区二区三 | 亚洲人成电影免费在线| 国产老妇伦熟女老妇高清| 久久精品人人爽人人爽视色| 日本av手机在线免费观看| 免费少妇av软件| 欧美日韩福利视频一区二区| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| 窝窝影院91人妻| 他把我摸到了高潮在线观看 | 亚洲精品av麻豆狂野| 中国美女看黄片| 久久香蕉激情| 黄色视频在线播放观看不卡| 人人妻人人澡人人爽人人夜夜| 夫妻午夜视频| 午夜视频精品福利| 性色av乱码一区二区三区2| 中文字幕色久视频| 日韩一区二区三区影片| 女人久久www免费人成看片| 精品国产乱码久久久久久小说| 少妇猛男粗大的猛烈进出视频| 一级黄色大片毛片| 女人久久www免费人成看片| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| h视频一区二区三区| 男女边摸边吃奶| 国产一区二区 视频在线| 国产一级毛片在线| 老司机影院毛片| 久久久久久亚洲精品国产蜜桃av| 在线 av 中文字幕| a级片在线免费高清观看视频| 一区二区三区精品91| 亚洲av电影在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 捣出白浆h1v1| 国产成人一区二区三区免费视频网站| 国产精品熟女久久久久浪| 精品一品国产午夜福利视频| 制服诱惑二区| 一本大道久久a久久精品| 久久狼人影院| 久久久国产成人免费| 99久久综合免费| 新久久久久国产一级毛片| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久| 搡老岳熟女国产| 亚洲伊人久久精品综合| 亚洲av成人不卡在线观看播放网 | 99久久综合免费| 日本黄色日本黄色录像| 男女国产视频网站| 国产av一区二区精品久久| 十分钟在线观看高清视频www| 欧美日韩亚洲高清精品| 汤姆久久久久久久影院中文字幕| 少妇 在线观看| 国产福利在线免费观看视频| 老司机亚洲免费影院| 久久女婷五月综合色啪小说| 精品高清国产在线一区| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图| 99国产极品粉嫩在线观看| 中文字幕人妻熟女乱码| 久久久久久人人人人人| 亚洲美女黄色视频免费看| 777久久人妻少妇嫩草av网站| 国产亚洲精品久久久久5区| 色播在线永久视频| 日韩熟女老妇一区二区性免费视频| 国产日韩一区二区三区精品不卡| 1024香蕉在线观看| 正在播放国产对白刺激| 国产一卡二卡三卡精品| 成人18禁高潮啪啪吃奶动态图| 久久精品aⅴ一区二区三区四区| 成年女人毛片免费观看观看9 | 精品国内亚洲2022精品成人 | 桃花免费在线播放| 欧美+亚洲+日韩+国产| 久久香蕉激情| 999精品在线视频| 在线观看免费日韩欧美大片| 丰满饥渴人妻一区二区三| 视频在线观看一区二区三区| 美女国产高潮福利片在线看| 啪啪无遮挡十八禁网站| 一边摸一边抽搐一进一出视频| 日韩欧美免费精品| av不卡在线播放| 叶爱在线成人免费视频播放| 亚洲精品国产av成人精品| 国产精品自产拍在线观看55亚洲 | 日本猛色少妇xxxxx猛交久久| 国产在线观看jvid| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| 自拍欧美九色日韩亚洲蝌蚪91| 成人国产一区最新在线观看| 国产极品粉嫩免费观看在线| 狠狠精品人妻久久久久久综合| 欧美成狂野欧美在线观看| 一区二区三区激情视频| 国产一区二区三区av在线| 性色av一级| 国产精品偷伦视频观看了| 国产精品国产三级国产专区5o| 亚洲精品日韩在线中文字幕| 丰满饥渴人妻一区二区三| 在线观看免费日韩欧美大片| 国产一级毛片在线| 欧美日韩一级在线毛片| 韩国精品一区二区三区| 国产成人免费观看mmmm| av天堂在线播放| 国产国语露脸激情在线看| 亚洲全国av大片| 中文字幕人妻熟女乱码| 国产深夜福利视频在线观看| 女人精品久久久久毛片| 巨乳人妻的诱惑在线观看| 国产男人的电影天堂91| 精品人妻在线不人妻| 美女扒开内裤让男人捅视频| 国产免费av片在线观看野外av| 男女高潮啪啪啪动态图| a级片在线免费高清观看视频| 国产亚洲午夜精品一区二区久久| 伊人久久大香线蕉亚洲五| 精品福利永久在线观看| 久久人人爽av亚洲精品天堂| 91字幕亚洲| 中国美女看黄片| 国产av国产精品国产| 日韩 亚洲 欧美在线| 精品国产一区二区三区四区第35| 久久人妻熟女aⅴ| 狂野欧美激情性xxxx| 欧美性长视频在线观看| 亚洲欧美色中文字幕在线| 精品免费久久久久久久清纯 | 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 老司机影院毛片| 女人精品久久久久毛片| 日韩免费高清中文字幕av| 少妇精品久久久久久久| 中国美女看黄片| 又黄又粗又硬又大视频| 一个人免费在线观看的高清视频 | 国产成人欧美| 超色免费av| 岛国在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 真人做人爱边吃奶动态| 90打野战视频偷拍视频| 热99国产精品久久久久久7| 亚洲av欧美aⅴ国产| 久久人人爽人人片av| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产av蜜桃| 日韩一卡2卡3卡4卡2021年| 精品亚洲成国产av| av又黄又爽大尺度在线免费看| 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 国产三级黄色录像| 精品国内亚洲2022精品成人 | 在线亚洲精品国产二区图片欧美| 久久这里只有精品19| 少妇精品久久久久久久| 高清在线国产一区| 日韩三级视频一区二区三区| 亚洲精品国产av成人精品| 欧美在线一区亚洲| 亚洲欧美激情在线| 18禁黄网站禁片午夜丰满| 日本91视频免费播放| 国产激情久久老熟女| 亚洲第一青青草原| 窝窝影院91人妻| 俄罗斯特黄特色一大片| 男人爽女人下面视频在线观看| 最新的欧美精品一区二区| 蜜桃国产av成人99| 亚洲av男天堂| 国产av一区二区精品久久| 色视频在线一区二区三区| 国产成人精品在线电影| 一边摸一边抽搐一进一出视频| 亚洲精品久久久久久婷婷小说| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 欧美黑人精品巨大| 黄色a级毛片大全视频| 性高湖久久久久久久久免费观看| 亚洲精品国产av成人精品| 久久久久久久久久久久大奶| 啦啦啦啦在线视频资源| 久久久精品国产亚洲av高清涩受| 青春草视频在线免费观看| 丰满迷人的少妇在线观看| 成人免费观看视频高清| 岛国在线观看网站| 国产一区二区三区av在线| av天堂在线播放| 国产精品免费大片| 深夜精品福利| 人妻人人澡人人爽人人| 亚洲精品日韩在线中文字幕| 亚洲专区中文字幕在线| 久久久国产一区二区| 美女高潮喷水抽搐中文字幕| 咕卡用的链子| 国产精品 国内视频| 久久久久久免费高清国产稀缺| 国产麻豆69| h视频一区二区三区| 亚洲欧美日韩高清在线视频 | 国产成人免费观看mmmm| 国产精品久久久人人做人人爽| 啦啦啦在线免费观看视频4| 免费看美女性在线毛片视频| 在线观看免费视频日本深夜| 国产成人av教育| 搡老岳熟女国产| 久久人妻av系列| 国产一区二区三区在线臀色熟女| 午夜激情av网站| 麻豆国产av国片精品| 日韩成人在线观看一区二区三区| 观看免费一级毛片| 99国产综合亚洲精品| 午夜激情av网站| 成人av在线播放网站| 亚洲成人中文字幕在线播放| 国产成人精品久久二区二区免费| 久久亚洲真实| 亚洲专区中文字幕在线| 亚洲avbb在线观看| 黄色视频不卡| 国产精品一区二区三区四区免费观看 | 亚洲国产精品sss在线观看| 日韩欧美国产在线观看| videosex国产| 国产精品免费视频内射| 操出白浆在线播放| 国内揄拍国产精品人妻在线| 一级毛片精品| 国产伦一二天堂av在线观看| 波多野结衣高清作品| 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 国产99白浆流出| 性欧美人与动物交配| 国产精品,欧美在线| 在线观看舔阴道视频| 日韩高清综合在线| 人人妻人人澡欧美一区二区| 欧美大码av| 亚洲精品国产精品久久久不卡| 岛国在线观看网站| 国产真实乱freesex| 五月伊人婷婷丁香| 午夜精品一区二区三区免费看| 欧美日韩乱码在线| 亚洲中文字幕一区二区三区有码在线看 | 首页视频小说图片口味搜索| 在线免费观看的www视频| 一本综合久久免费| 欧美乱妇无乱码| 美女大奶头视频| 高潮久久久久久久久久久不卡| 黄频高清免费视频| 久久婷婷成人综合色麻豆| 国产区一区二久久| av国产免费在线观看| 日本在线视频免费播放| 国产蜜桃级精品一区二区三区| 在线播放国产精品三级| 久久精品国产清高在天天线| svipshipincom国产片| 日韩欧美一区二区三区在线观看| 亚洲精品国产一区二区精华液| 久久99热这里只有精品18| 国产aⅴ精品一区二区三区波| 欧美黄色片欧美黄色片| 中文字幕久久专区| 亚洲中文av在线| 91在线观看av| 亚洲精品中文字幕一二三四区| 激情在线观看视频在线高清| 欧美黑人巨大hd| 麻豆成人av在线观看| 成年人黄色毛片网站| 午夜视频精品福利| aaaaa片日本免费| 巨乳人妻的诱惑在线观看| 国产久久久一区二区三区| 国产蜜桃级精品一区二区三区| 国产在线精品亚洲第一网站| 亚洲国产欧洲综合997久久,| 99精品在免费线老司机午夜| 亚洲狠狠婷婷综合久久图片| 亚洲va日本ⅴa欧美va伊人久久| 母亲3免费完整高清在线观看| 夜夜看夜夜爽夜夜摸| 亚洲人成电影免费在线| 亚洲美女黄片视频| 欧美国产日韩亚洲一区| 亚洲,欧美精品.| 黄色毛片三级朝国网站| 又大又爽又粗| 搡老岳熟女国产| 国产1区2区3区精品| 亚洲国产精品成人综合色| 国产精品综合久久久久久久免费| 国内久久婷婷六月综合欲色啪| 国产一区二区三区在线臀色熟女| 在线观看www视频免费| 午夜老司机福利片| 国产真人三级小视频在线观看| 中文字幕久久专区| 亚洲精华国产精华精| 亚洲成人久久性| 成人亚洲精品av一区二区| a级毛片a级免费在线| 精品一区二区三区视频在线观看免费| 亚洲av电影不卡..在线观看| or卡值多少钱|