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

    標(biāo)記引物數(shù)量對(duì)白樺遺傳多樣性的影響

    2016-05-11 01:43:46趙雙菁劉瑩瑩魏敏靜楊傳平魏志剛
    防護(hù)林科技 2016年4期
    關(guān)鍵詞:遺傳多樣性白樺

    趙雙菁,劉瑩瑩,魏敏靜,楊傳平,魏志剛

    〔林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室(東北林業(yè)大學(xué)),黑龍江 哈爾濱 150040〕

    ?

    標(biāo)記引物數(shù)量對(duì)白樺遺傳多樣性的影響

    趙雙菁,劉瑩瑩,魏敏靜,楊傳平,魏志剛

    〔林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室(東北林業(yè)大學(xué)),黑龍江 哈爾濱 150040〕

    摘要以帽兒山實(shí)驗(yàn)林場(chǎng)15個(gè)種源的白樺種源試驗(yàn)林300個(gè)樣本為試驗(yàn)材料,通過(guò)不同數(shù)量的SRAP引物對(duì)白樺遺傳多樣性研究結(jié)果的影響做出分析,并提出標(biāo)記引物確定的計(jì)算公式。結(jié)果表明,SRAP引物數(shù)量不影響白樺各種源間遺傳參數(shù)的變化趨勢(shì),但隨著引物數(shù)量的增加,區(qū)分不同種源間遺傳差異的能力隨之提高,當(dāng)SRAP引物數(shù)量為15對(duì)時(shí),可探測(cè)到的多態(tài)位點(diǎn)百分率達(dá)98.5%以上。

    關(guān)鍵詞白樺;標(biāo)記引物數(shù)量;遺傳多樣性

    遺傳多樣性是生物多樣性的基礎(chǔ)和核心,是生態(tài)系統(tǒng)多樣性和物種多樣性的基礎(chǔ)方面,保護(hù)生物多樣性最終是要保護(hù)其遺傳多樣性[1]。林木遺傳多樣性的研究開(kāi)始于形態(tài)學(xué)為基礎(chǔ)的種源試驗(yàn),如杉木(Cunninghamialanceolata)、桉樹(shù)(Eucalyptusrobusta)、紅松(Pinuskoraiensis)、樟子松(Pinussylvestrisvar.mongolica)和長(zhǎng)白落葉松(Larixolgensis)等樹(shù)種的種源試驗(yàn)研究[2,3]。同工酶標(biāo)記發(fā)明后,被先后運(yùn)用到歐洲赤松(Pinussylvestris)、云杉(Piceaasperata)、紅松、顫楊(Populustremuloide)、花旗松(Pseudotsugamenziesii)和歐洲山毛櫸(Fagussylvatica)等重要林木遺傳多樣性研究中[4-6]。隨后,DNA標(biāo)記的出現(xiàn)及其技術(shù)自身的優(yōu)點(diǎn),被迅速運(yùn)用到林木遺傳多樣性研究中并極大地推動(dòng)了該領(lǐng)域的研究進(jìn)程,如歐洲赤松、栲樹(shù)(Castanopsisfargesii)、挪威云杉(Piceaabies)、山毛櫸、歐洲刺槐(Robiniapseudoacacia)、云杉、海岸松(Pinuspinaster)、雪松(Cedrusdeodara)、紅松、木荷(Schimasuperba)、長(zhǎng)葉榧(Torreyajackii)和白皮松(Pinusbungeana)等一些主要用材和珍稀瀕危樹(shù)種均采用DNA標(biāo)記進(jìn)行了研究,是林木遺傳多樣性研究中應(yīng)用最廣泛的技術(shù)[7]。目前為止,雖然大規(guī)模測(cè)序技術(shù)的不斷發(fā)展和完善,相繼有楊樹(shù)、桉樹(shù)和歐洲赤松等少數(shù)幾個(gè)樹(shù)種完成了全基因測(cè)序,這為基于SNP技術(shù)為基礎(chǔ)的林木遺傳多樣性及輔助選擇育種提供了可能,然而對(duì)于多數(shù)植物與樹(shù)木而言仍然缺乏全基因組序列,因此,目前基于基因組非編碼區(qū)變異變?yōu)榛A(chǔ)的DNA標(biāo)記技術(shù)仍是多數(shù)植物遺傳多樣性研究中的主要方法與技術(shù)。

    群體遺傳學(xué)的研究基礎(chǔ)是DNA序列變異。同源DNA序列的遺傳分化程度是衡量群體遺傳結(jié)構(gòu)的主要指標(biāo),其分化式樣則是理解群體遺傳結(jié)構(gòu)產(chǎn)生和維持的進(jìn)化內(nèi)在驅(qū)動(dòng)力諸如遺傳突變、重組、基因轉(zhuǎn)換的前提[8]。因此,在遺傳多樣性研究中,要盡可能地揭示出該物種全基因組區(qū)域的DNA序列的變異,這樣才能全面反映出其遺傳多樣性和群體遺傳結(jié)構(gòu),為后期資源的保存和利用奠定基礎(chǔ)。DNA分子標(biāo)記,不同的標(biāo)記及引物反映的基因組區(qū)域的不同DNA序列變異,因此如果標(biāo)記數(shù)量較少時(shí),可能無(wú)法全部覆蓋基因組上產(chǎn)生變異的所有區(qū)域,因此不可能揭示該變異的全部信息[9],而數(shù)量過(guò)多,不能增加信息量還會(huì)造成額外的實(shí)驗(yàn)成本。因此,對(duì)于遺傳多樣性研究中,究竟采用多少標(biāo)記足以揭示其群體的遺傳結(jié)構(gòu)是一個(gè)值得重視的問(wèn)題。如Julio等在研究薔薇的遺傳多樣性時(shí)發(fā)現(xiàn),對(duì)于雜合度的估算,必須有10個(gè)以上的個(gè)體和20個(gè)以上的標(biāo)記[10]。而Reyazu等在黃麻(Corchorusspecies)遺傳多樣性研究中發(fā)現(xiàn),15對(duì)與41對(duì)SSRs引物揭示的信息量相同[11]。白樺(Betulaplatyphylla)是一種分布范圍廣、適用性強(qiáng)、用途廣泛的闊葉速生樹(shù)種,然而,有關(guān)白樺遺傳多樣性研究報(bào)道中 ,同樣也沒(méi)有重視到標(biāo)記數(shù)量對(duì)遺傳參數(shù)結(jié)果與遺傳結(jié)果影響的問(wèn)題。

    SRAP ( Sequence-related amplified polymorphism,序列相關(guān)擴(kuò)增多態(tài)性) 標(biāo)記建立以來(lái),由于其操作簡(jiǎn)便快速、成本低、可信度高、易于測(cè)序等特點(diǎn)備受分子生物學(xué)家的青睞[9],被迅速應(yīng)用到蓮藕(Nelumbonucif)[10]、桃(Amygdaluspersica)[11]花生(Arachishypogaea)[12]、油菜(Brassicanapus)[13]、黃瓜(Cucumissativus)[14]、甘薯(Dioscoreaesculenta)[15]和西瓜(Citrulluslanatus)[16]等植物的種質(zhì)資源鑒定評(píng)價(jià)、遺傳圖譜構(gòu)建(包括轉(zhuǎn)錄圖譜)、重要性狀標(biāo)記乃至基因分離克隆等方面得到了應(yīng)用。本文通過(guò)對(duì)不同SRAP引物對(duì)白樺遺傳多樣性研究結(jié)果的影響的研究,提出SRAP標(biāo)記研究白樺遺傳多樣性引物數(shù)量的標(biāo)記,在此基礎(chǔ)上提出標(biāo)記引物確定的計(jì)算公式。

    1材料與方法

    1.1試驗(yàn)材料來(lái)源

    試驗(yàn)材料來(lái)自東北林業(yè)大學(xué)帽兒山實(shí)驗(yàn)林場(chǎng)白樺種源試驗(yàn)林,共15個(gè)種源,分別為:新疆、涼水、汪清、露水河、小北湖、輝南、桓仁、草河口、綽爾、寧夏、莫爾道嘎、清源、東方紅、帽兒山、烏伊嶺。每個(gè)種源隨機(jī)選取20個(gè)個(gè)體,共計(jì)300個(gè)樣本。2012年6月采取當(dāng)年生嫩葉,用密封袋封好,做好標(biāo)記,放在冰盒中帶回實(shí)驗(yàn)室,置于-80 ℃冰箱中保存。

    1.2試驗(yàn)方法

    1.2.1白樺葉片基因組DNA提取白樺葉片基因組DNA提取方法見(jiàn)參考文獻(xiàn)[12]。

    1.2.2反應(yīng)體系及反應(yīng)程序SRAP該反應(yīng)體系為20 μL:0.21×106μgL-1DNA,1.5 μL;25 mmolL-1Mg2+,1.4 μL;5 UμL-1Taq酶,0.25 μL;2.5 mmolL-1dNTPs,2 μL;10 μmL-1引物,0.35 μL。PCR反應(yīng)程序?yàn)椋?4 ℃預(yù)變性5 min,94 ℃變性1 min,35 ℃復(fù)性1 min,72 ℃延伸1 min,5個(gè)循環(huán);94 ℃變性1 min,50 ℃復(fù)性1 min,72 ℃延伸1 min,30個(gè)循環(huán),72 ℃延伸7 min[13,14]。

    1.2.3引物篩選本研究中所用SRAP引物根據(jù)文獻(xiàn)[15]中的引物序列,在上海生工合成。從15個(gè)種源中選取地理位置最近的2個(gè)種源(莫爾道嘎和桓仁),每個(gè)種源選5個(gè)個(gè)體用于引物的篩選。以條帶清晰并且群體與個(gè)體間至少有一條差異條帶的引物為選擇標(biāo)準(zhǔn),從合成的30對(duì)引物中,篩選出17對(duì)條帶清晰、多態(tài)性高的引物用于實(shí)驗(yàn)分析。引物序列如表1所示。

    表1 SRAP引物序列

    1.2.4SRAP引物數(shù)量按5、7、9、11、13、15、17對(duì)SRAP引物對(duì)15個(gè)白樺種源進(jìn)行擴(kuò)增和分析。

    1.3數(shù)據(jù)統(tǒng)計(jì)及分析

    電泳圖譜中的每一條帶均代表了引物與模板DNA互補(bǔ)的一對(duì)結(jié)合位點(diǎn),可記為一個(gè)分子標(biāo)記,根據(jù)分子量標(biāo)準(zhǔn)對(duì)照反應(yīng)產(chǎn)物在膠上的對(duì)應(yīng)位置,估計(jì)擴(kuò)增產(chǎn)物的分子量大小,有帶的記為“1”,無(wú)帶的記為“0”[16]。所得結(jié)果為一二元數(shù)據(jù)矩陣,利用Popgene32軟件計(jì)算各項(xiàng)遺傳參數(shù)。

    2結(jié)果與分析

    2.1不同標(biāo)記引物數(shù)量的白樺群體遺傳多樣性

    按5、7、9、11、13、15、17對(duì)SRAP引物對(duì)15個(gè)白樺種源、300個(gè)個(gè)體進(jìn)行SRAP擴(kuò)增并對(duì)各遺傳參數(shù)統(tǒng)計(jì)分析,結(jié)果見(jiàn)表2。表2表明,引物數(shù)量不同時(shí),得到的白樺群體遺傳多樣性各項(xiàng)參數(shù)的大小也不同。如5對(duì)特異性引物(me2em3、me2em4、me2em5、me2em6、me2em8)共檢測(cè)到46個(gè)位點(diǎn),其中多態(tài)位點(diǎn)有45,多態(tài)位點(diǎn)百分率達(dá)97.83%,多態(tài)位點(diǎn)數(shù)量為21~33、多態(tài)位點(diǎn)百分率在45.65%~71.74%,平均等位基因數(shù)在1.456 5~1.717 4,有效等位基因數(shù)在1.200 8~1.333 9,Nei’s指數(shù)在0.123 1~0.202 6,Shannon指數(shù)在0.192 6~0.314 1,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,最小為草河口種源。白樺種源總基因多樣性為0.234 3,種群內(nèi)基因多樣性為0.168 2,種群間的遺傳多樣性為0.066 1,種源間分化系數(shù)為0.283 0,說(shuō)明71.70%的遺傳變異存在于種源內(nèi),28.30%存在于種源間,基因流Nm為 1.272 8。

    7對(duì)SRAP引物(me2em3、me2em4、me2em5、me2em6、me2em8、me2em11、me3em3)共檢測(cè)到62個(gè)位點(diǎn),其中多態(tài)位點(diǎn)有61,多態(tài)位點(diǎn)百分率為98.38%,多態(tài)位點(diǎn)數(shù)量在26~61之間、多態(tài)位點(diǎn)百分率在41.94%~98.38%,平均等位基因數(shù)在1.419 4~1.709 7,有效等位基因數(shù)在1.205 8~1.295 3、Nei’s指數(shù)在0.125 2~0.214 1、Shannon指數(shù)在0.193 1~0.320 6,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大與最小種源分別為涼水和淖爾種源。白樺種源總基因多樣性為0.243 7,種群內(nèi)和種群間基因多樣性為0.166 1和0.077 6,種源間分化系數(shù)為0.318 3,說(shuō)明68.16%的遺傳變異存在于種源內(nèi),31.84%存在于種源間,基因流Nm為1.070 6。

    9對(duì)SRAP引物(me2em3、me2em4、me2em5、me2em6、me2em8、me2em11、me3em3、me3em4、me4em9)檢測(cè)到白樺各種源平均等位基因數(shù)為1.369 9~1.726 0,有效等位基因數(shù)1.176 3~1.373 9,Nei’s指數(shù)在0.107 6~0.216 4,Shannon指數(shù)在0.166 8~0.325 3,多態(tài)位點(diǎn)數(shù)量為27~72、多態(tài)位點(diǎn)百分率達(dá)36.99%~98.63%,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,而最小為淖爾種源。白樺種源總基因多樣性為0.248 8,種群內(nèi)基因多樣性Hs為0.163 8,種群間的遺傳多樣性為0.077 6,種源間分化系數(shù)為0.341 9,說(shuō)明65.84%的遺傳變異存在于種源內(nèi),34.16%存在于種源間,基因流Nm為0.962 5。

    11對(duì)特異性引物(me2em3、me2em4、me2em5、me2em6、me2em8、me2em11、me3em3、me3em4、me4em9、me4em11、me5em4)檢測(cè)到白樺各種源平均等位基因數(shù)為1.348 3~1.752 8,有效等位基因數(shù)為1.165 8~1.367 2,Nei’s指數(shù)為0.101 3~0.223 4,Shannon指數(shù)為0.157 1~0.344 2,多態(tài)位點(diǎn)數(shù)量為31~87,多態(tài)位點(diǎn)百分率為34.8%~75.68%,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,而最小為淖爾種源。白樺種源總基因多樣性為0.248 1,種群內(nèi)基因多樣性為0.162 9,種群間的遺傳多樣性為0.085 2,種源間分化系數(shù)為0.343 5,說(shuō)明65.66%的遺傳變異存在于種源內(nèi),34.34%存在于種源間,基因流Nm為0.955 7。

    13對(duì)特異性引物(me2em3、me2em4、me2em5、me2em6、me2em8、me2em11、me3em3、me3em4、me4em9、me4em11、me5em4、me5em11、me6em12) 檢測(cè)到白樺各種源平均等位基因數(shù)為1.403 8~1.769 2,有效等位基因數(shù)1.191 3~1.379 0,Nei’s指數(shù)在0.118 5~0.229 1,Shannon指數(shù)0.184 0~0.352 0,多態(tài)位點(diǎn)數(shù)量為42~80,多態(tài)位點(diǎn)百分率達(dá)40.38%~76.92%,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,而最小為淖爾種源。白樺種源總基因多樣性為0.255 6,種群內(nèi)基因多樣性為0.168 4,種群間的遺傳多樣性為0.087 2,種源間分化系數(shù)為0.341 1,說(shuō)明65.88%的遺傳變異存在于種源內(nèi),34.12%存在于種源間,但基因流Nm為0.966 0。

    15對(duì)特異性引物(me2em3、me2em4、me2em5、me2em6、me2em8、me2em11、me3em3、me3em4、me4em9、me4em11、me5em4、me5em11、me6em12)檢測(cè)到白樺各種源平均等位基因數(shù)為1.408 3~1.791 7,有效等位基因數(shù)為1.211 6~1.396 4,Nei’s指數(shù)為0.125 3~0.239 3,Shannon指數(shù)為0.190 9~0.367 0,多態(tài)位點(diǎn)數(shù)量為42~80,多態(tài)位點(diǎn)百分率為40.38%~76.92%,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,而最小為淖爾種源。白樺種源總基因多樣性為0.262 3,種群內(nèi)基因多樣性Hs為0.173 0,種群間的遺傳多樣性為0.089 3,種源間分化系數(shù)為0.340 6,說(shuō)明65.96%的遺傳變異存在于種源內(nèi),34.05%存在于種源間,但基因流Nm為0.967 9。

    利用篩選出的17對(duì)特異性SRAP引物檢測(cè)到白樺各種源平均等位基因數(shù)為1.363 0~1.829 6,有效等位基因數(shù)為1.188 1~1.415 3,Nei’s指數(shù)為0.111 4~0.248 3,Shannon指數(shù)為0.169 7~0.379 9,多態(tài)位點(diǎn)數(shù)量為49~112,多態(tài)位點(diǎn)百分率為36.3%~82.96%,各種源上述遺傳參數(shù)的變化趨勢(shì)基本一致,最大均為涼水種源,最小為清源種源。白樺種源總基因多樣性為0.2591,種群內(nèi)基因多樣性為0.170 7,種群間的遺傳多樣性為0.088 4,種源間分化系數(shù)為0.340 9,說(shuō)明65.88%的遺傳變異存在于種源內(nèi),34.12%存在于種源間,基因流Nm為0.966 5。

    表2 不同引物數(shù)量遺傳參數(shù)比較

    注:a為引物數(shù)量為5對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,b為引物數(shù)量7對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,c為引物數(shù)量9對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,d為引物數(shù)量11對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,e為引物數(shù)量13對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,f為引物數(shù)量15對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果,g為引物數(shù)量17對(duì)時(shí)白樺種源的聚類(lèi)結(jié)果

    2.2白樺種源不同樣本數(shù)量時(shí)群體間親緣關(guān)系

    為了進(jìn)一步分析引物數(shù)量對(duì)于群體遺傳結(jié)構(gòu)的分析結(jié)果的影響,對(duì)不同標(biāo)記數(shù)量分析了白樺15個(gè)種源間的遺傳距離并進(jìn)行了聚類(lèi)分析。結(jié)果表明,引物數(shù)量不同,不僅對(duì)不同種源間的遺傳距離的估算結(jié)果產(chǎn)生影響,而且得到的白樺種源間進(jìn)化關(guān)系存在一定的差異(圖1)。不同引物數(shù)量,反映出遺傳距離最近與最遠(yuǎn)的種源也存在差異,如當(dāng)引物數(shù)為7、13、15和17對(duì)時(shí),遺傳距離最小的2個(gè)種源均為莫爾道嘎和桓仁;而引物數(shù)為5、9和11對(duì)時(shí),遺傳距離最小的2個(gè)種源為清源和寧夏。引物數(shù)量為9、11、13、15和17對(duì)時(shí),遺傳距離最大的2個(gè)種源均為涼水和草河口,而引物數(shù)量為7對(duì)時(shí),涼水和東方紅種源的遺傳距離最大;引物數(shù)量為5對(duì)時(shí),涼水和烏伊嶺種源的遺傳距離最大。從上述不同引物數(shù)量白樺種源聚類(lèi)結(jié)果來(lái)看,在相同的遺傳距離處(如0.11處),不同引物數(shù)量導(dǎo)致白樺種源聚類(lèi)類(lèi)群數(shù)量各不相同,如當(dāng)引物數(shù)量為5對(duì)時(shí),15個(gè)種源聚成3類(lèi);當(dāng)引物數(shù)量為5對(duì)時(shí),15個(gè)種源聚成3大類(lèi)群;引物數(shù)量為7和9對(duì)時(shí),聚成4類(lèi)。同時(shí),上述不同聚類(lèi)結(jié)果中,包含的種源也各不相同,總體來(lái)說(shuō),隨著引物數(shù)的增多,聚類(lèi)數(shù)量與每個(gè)類(lèi)別中包括的種源趨于相同。如當(dāng)引物數(shù)量為5對(duì)時(shí), 15個(gè)種源可聚成3個(gè)類(lèi)群:涼水種源聚為一類(lèi),烏伊嶺聚為一類(lèi),其余種源聚成一大類(lèi)(圖1 a)。當(dāng)引物數(shù)量為7對(duì)時(shí),可聚成4個(gè)類(lèi)群:涼水種源聚為一類(lèi),輝南聚為一類(lèi),汪清、新疆、露水河聚成一類(lèi),其余種源聚成一大類(lèi)(圖1 b)。當(dāng)引物數(shù)量為9對(duì)時(shí),可聚成4個(gè)類(lèi)群:涼水種源聚為一類(lèi),輝南聚為一類(lèi),汪清、新疆、露水河聚成一類(lèi),其余種源聚成一大類(lèi)(圖1 c)。當(dāng)引物數(shù)量為11對(duì)時(shí),可聚成5個(gè)類(lèi)群:涼水種源單獨(dú)聚為一類(lèi);汪清、新疆、露水河種源聚為一類(lèi);小北湖、桓仁、莫爾道嘎、輝南聚為一類(lèi);草河口、淖爾、寧夏、清源種源聚為一類(lèi),東方紅、烏伊嶺、帽兒山種源聚為一類(lèi)(圖1d)。當(dāng)引物數(shù)量為13對(duì)時(shí),可聚成5個(gè)類(lèi)群:涼水、汪清種源聚為一類(lèi),新疆、露水河種源聚為一類(lèi),小北湖、輝南、莫爾道嘎和桓仁種源聚為一類(lèi),草河口、淖爾、寧夏、清源種源聚為一類(lèi)(圖1 e)。當(dāng)引物數(shù)量為15對(duì)時(shí),可聚成6個(gè)類(lèi)群:涼水、汪清種源聚為一類(lèi);新疆、露水河種源聚為一類(lèi),桓仁為一類(lèi)、東方紅、烏伊嶺、帽兒山聚為一類(lèi),小北湖、輝南、莫爾道嘎種源聚為一類(lèi),草河口、淖爾、寧夏、清源種源聚為一類(lèi)(圖1 f)。當(dāng)引物數(shù)量為17對(duì)時(shí),可聚成5個(gè)類(lèi)群:涼水、汪清種源聚為一類(lèi);新疆、露水河種源聚為一類(lèi);小北湖、輝南、桓仁、莫爾道嘎種源聚為一類(lèi);草河口、淖爾、寧夏、清源種源聚為一類(lèi);東方紅、烏伊嶺、帽兒山聚為一類(lèi)(圖1 g)。

    2.3引物數(shù)量和白樺遺傳參數(shù)的相關(guān)分析

    為了進(jìn)一步弄清SRAP引物數(shù)量對(duì)白樺主要遺傳參數(shù)影響的大小,對(duì)引物數(shù)量和各遺傳參數(shù)進(jìn)行了相關(guān)性分析(表3),結(jié)果表明,引物數(shù)量白樺群體的有效等位基因數(shù)量Ne、Nei’s基因多樣性指數(shù)H、Shannon’s信息多樣性指數(shù)I、多態(tài)位點(diǎn)、多態(tài)位點(diǎn)百分率和群體基因多樣性間呈極顯著相關(guān),而群體間分化系數(shù)與基因流與引物數(shù)量之間顯著相關(guān),而與其他遺傳參數(shù)之間無(wú)顯著相關(guān)。

    表3 SRAP引物數(shù)量與主要遺傳參數(shù)的相關(guān)性

    注:表格內(nèi)上面數(shù)字為相關(guān)系數(shù),下面為顯著性水平

    2.4白樺遺傳多樣性引物數(shù)量的確定

    在群體遺傳多樣性研究中,檢測(cè)多態(tài)位點(diǎn)數(shù)越高越有利于分析和研究群體遺傳多樣性和遺傳結(jié)構(gòu),因此對(duì)白樺種源不同樣本數(shù)量與多態(tài)位點(diǎn)百分率進(jìn)行擬合分析(圖2),在95%的置信區(qū)間內(nèi),兩者之間極顯著線性回歸,其決定系數(shù)為93.9%,表明兩者之間存在較好的相關(guān)性。

    我們?cè)诎讟暹z傳多樣性研究時(shí),可在綜合考慮成本、技術(shù)難度等基礎(chǔ)上,根據(jù)所需多態(tài)位點(diǎn)百分率大小,計(jì)算出所需樣本的數(shù)量。如當(dāng)白樺遺傳多樣性研究時(shí),以最低所需多態(tài)位點(diǎn)百分率98.5%為標(biāo)準(zhǔn),則每個(gè)種源所需引物數(shù)量為14、16對(duì)時(shí),實(shí)際操作中可設(shè)為15對(duì)標(biāo)記。

    3結(jié)論

    3.1SRAP引物數(shù)量對(duì)白樺種源遺傳參數(shù)影響不同。其中,引物數(shù)量與白樺群體的有效等位基因數(shù)量、Nei’s基因多樣性指數(shù)、Shannon’s信息多樣性指數(shù)I、多態(tài)位點(diǎn)、多態(tài)位點(diǎn)百分率和群體基因多樣性間呈極顯著相關(guān),而群體間分化系數(shù)與基因流與引物數(shù)量之間顯著相關(guān),而與其他遺傳參數(shù)之間無(wú)顯著相關(guān)。不同標(biāo)記數(shù)量與存在顯著差異的遺傳參數(shù)之間均呈正相關(guān),然而,除5個(gè)SARP標(biāo)記數(shù)時(shí),15個(gè)白樺種源內(nèi)多態(tài)位百分率最低為草河口種源,而7個(gè)標(biāo)記以上時(shí),分析發(fā)現(xiàn)多態(tài)位百分率最低均為淖爾種源。

    除此之外,不同標(biāo)記數(shù)量反映出的白樺不同種源間各遺傳參數(shù)的變化趨勢(shì)基本相似,表明不同引物數(shù)量并不影響白樺各種源間遺傳參數(shù)的變化趨勢(shì)。

    3.2SRAP不同引物數(shù)量對(duì)白樺15個(gè)種源間進(jìn)化關(guān)系的分析結(jié)果有較大影響。當(dāng)引物數(shù)量分別為5和7時(shí),反映出的遺傳距離最小與最大的種源均不同,當(dāng)引物數(shù)量為9和11對(duì)時(shí),反映出的遺傳距離最小種源不同,但遺傳距離最大的2個(gè)種源相同,均為涼水和草河口;引物數(shù)量超過(guò)13對(duì)時(shí),種源間遺傳距離最小與最大的2對(duì)種源分別為莫爾道嘎和桓仁、涼水和草河口。而且,隨著SRAP引物數(shù)量增加,在遺傳距離為0.11處,分別可將15個(gè)種源聚類(lèi)成3、4、4、5、6、5和5個(gè)類(lèi)別,表明隨著引物數(shù)量的增加,區(qū)分不同種源間遺傳差異的能力總體上也隨之提高。

    3.3SRAP引物數(shù)量與其多態(tài)位點(diǎn)百分率呈線性相關(guān),且決定系數(shù)達(dá)93.9%。并以此為公式計(jì)算出開(kāi)展白樺遺傳多樣性研究時(shí),如設(shè)定要探測(cè)到的多態(tài)位點(diǎn)百分率達(dá)98.5%以上時(shí),白樺遺傳多樣性研究中SRAP引物數(shù)量設(shè)定為15個(gè)。

    參考文獻(xiàn):

    [1] 顧萬(wàn)春,李文英.我國(guó)林木種質(zhì)資源共享現(xiàn)狀及建議[J].世界林業(yè)研究,2007,20:66-69

    [2] 楊書(shū)文,楊傳平,王秋玉,等.東北內(nèi)蒙主要用材樹(shù)種的種源選擇[J].東北林業(yè)大學(xué)學(xué)報(bào),1993(4):1-7

    [3] 沈熙環(huán).林木育種成果亟待鞏固發(fā)展[J].中國(guó)林業(yè),2002(6):29-30

    [4] J.Buiteveld et al.Genetic diversity and differentiation in European beech (Fagus sylvatica L.) stands varying in management history[J].Forest Ecology and Management,2007,247: 98-106

    [5] 羅建勛,顧萬(wàn)春,陳少瑜.云杉天然群體遺傳多樣性的等位酶變異[J].植物生態(tài)學(xué)報(bào),2006(1):165-173

    [6] 張春曉,李?lèi)?沈熙環(huán).林木同工酶遺傳多樣性研究進(jìn)展[J].北京林業(yè)大學(xué)學(xué)報(bào),1998(3):61-69

    [7] A Roy, A Bandyopadhyay, A K Mahapatra, et al. Evaluation of genetic diversity in jute (Corc-horus species) using STMS, ISSR and RAPD markers[J]. Plant Breeding, 2006, 125(3): 292- 297.

    [8] 王云生,黃宏文,王瑛.植物分子群體遺傳學(xué)研究動(dòng)態(tài)[J].遺傳,2007(10):1191-1198

    [9] Navinder Sainil, Neelu Jainl, Sunita Jain, et a1. Assessment of genetic diversity within and among Basmati and non-Basmati rice varieties using AFLP, ISSR and SSR markers[J].Euphytica, 2004,140(1):133-146

    [10] Norma Julio, Ana Sobral, Juan Rondan Duen, et a1. RAPD and ISSR markers indicate diminished gene flow due to recent fragmentation of Polylepis australis woodlands in central Arge-ntina[J]. Biochemical Systematics and Ecology,2008,36(5-6),329-335

    [11] Reyazul R, Mir.Sachin Rustgi, Shailendra Sharma,et a1.A preliminary genetic analysis of fibre traits and the use of new genomic SSRs for genetic diversity in jute[J]. Euphytica, 2008,161(3):413-427

    [12] 魏志剛,楊傳平,張力杰,等. 一種白樺葉片基因組DNA提取方法[P].黑龍江:CN103013985A,2013-04-03

    [13] 楊傳平,潘華,魏志剛.白樺ISSR~PCR反應(yīng)體系的優(yōu)化[J].東北林業(yè)大報(bào),2003,3(6):1-3

    [14] 連蓮.白樺ISSR、AFLP遺傳圖譜構(gòu)建[D].哈爾濱:東北林業(yè)大學(xué),2008:26

    [15] Lie G F. Quiros Sequence-related amplified polymorphism(SRAP), a new marker system base on a simple PCR reaction:its application to mapping and gene tagging in Brassiea[J]. Theoretical and Applied Genetics, 2001,103(8):455-461

    [16] 李雪峰.興安落葉松種群、無(wú)性系遺傳多樣性的研究[D].哈爾濱:東北林業(yè)大學(xué),2009:22-23

    Effects of Number of Labeling Primers on Genetic Diversity ofBetulaplatyphylla

    Zhao Shuangjing, Liu Yingying, Wei Minjing, Yang Chuanping, Wei Zhigang

    〔State Key Laboratory of Forest Genetics and Breeding (Northeast Forestry University), Harbin 150040,China〕

    AbstractTaking 300 samples from 15 provenance in experimental forest farm for Betula platyphylla in Maoershan Experimental Forest Farm as experimental materials,the effects of different number of SRAP on primers genetic diversity of Betula platyphylla were analyzed. The formula for labeling primers was proposed. Result shows that the number of SRAP primer does not affect the trends for genetic parameters among provenances of Betula platyphylla. The ability of distinguishing genetic differences among different provenances increase with the increase of the number for labeling primers. While SRAP primer number being 15 pairs,polymorphic percentage which can be detected reach above 98.5%.

    Key wordsBetula platyphylla;number of labeling primers;genetic diversity

    中圖分類(lèi)號(hào):S792.153

    文獻(xiàn)標(biāo)識(shí)碼:A

    doi:10.13601/j.issn.1005-5215.2016.04.002

    作者簡(jiǎn)介:趙雙菁,碩士,Email:zshuangjing2006@163.com通訊作者:魏志剛,教授,博士,Email:zhigangwei1973@163.com

    基金項(xiàng)目:863課題白樺、桉樹(shù)等分子育種與品質(zhì)創(chuàng)制(2011AA100202)

    收稿日期:2016-03-21

    文章編號(hào):1005-5215(2016)04-0006-05

    猜你喜歡
    遺傳多樣性白樺
    《白樺》素養(yǎng)提升
    《白樺》素養(yǎng)提升
    白樺生北國(guó)
    文苑(2019年20期)2019-11-20 02:12:33
    白樺生北國(guó)
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    俄羅斯見(jiàn)聞:浴血白樺
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    亚洲熟女毛片儿| 色综合欧美亚洲国产小说| 一级作爱视频免费观看| 99riav亚洲国产免费| 夜夜爽天天搞| 久久九九热精品免费| 亚洲人成伊人成综合网2020| 国产亚洲欧美精品永久| 成人三级做爰电影| 成熟少妇高潮喷水视频| 精品日产1卡2卡| 91国产中文字幕| 亚洲欧美一区二区三区黑人| 欧美色视频一区免费| av网站在线播放免费| 欧美日韩国产mv在线观看视频| 在线十欧美十亚洲十日本专区| 999久久久精品免费观看国产| 亚洲五月色婷婷综合| 成人特级黄色片久久久久久久| 日韩一卡2卡3卡4卡2021年| 男人的好看免费观看在线视频 | 91精品三级在线观看| 欧美黑人欧美精品刺激| 成人国产一区最新在线观看| 99re在线观看精品视频| 亚洲成av片中文字幕在线观看| 久久精品国产亚洲av高清一级| 欧美人与性动交α欧美精品济南到| 午夜亚洲福利在线播放| 叶爱在线成人免费视频播放| 91成年电影在线观看| 午夜福利,免费看| 亚洲av成人av| 久久中文字幕一级| 日韩人妻精品一区2区三区| 如日韩欧美国产精品一区二区三区| 欧美一级毛片孕妇| 看黄色毛片网站| 久久久久久免费高清国产稀缺| 国产99白浆流出| 精品电影一区二区在线| 久久欧美精品欧美久久欧美| 亚洲精品久久成人aⅴ小说| 波多野结衣高清无吗| 亚洲成人精品中文字幕电影 | 亚洲熟妇中文字幕五十中出 | 亚洲中文日韩欧美视频| 亚洲午夜理论影院| 久久香蕉精品热| 国产一区在线观看成人免费| 日本黄色日本黄色录像| av在线天堂中文字幕 | 精品久久久久久久久久免费视频 | 国产精品久久久av美女十八| 国产91精品成人一区二区三区| 啦啦啦 在线观看视频| 国产蜜桃级精品一区二区三区| 亚洲美女黄片视频| 日韩三级视频一区二区三区| 国产精品久久久久久人妻精品电影| 欧美色视频一区免费| 色综合欧美亚洲国产小说| 99国产精品99久久久久| 中文亚洲av片在线观看爽| 久久人人精品亚洲av| 亚洲成av片中文字幕在线观看| 一夜夜www| 母亲3免费完整高清在线观看| 超色免费av| 后天国语完整版免费观看| 国产区一区二久久| 99久久精品国产亚洲精品| 国产真实乱freesex| 性色avwww在线观看| 成年女人永久免费观看视频| 久久久久久久久大av| 深夜a级毛片| 久久欧美精品欧美久久欧美| 身体一侧抽搐| 日本熟妇午夜| 波多野结衣巨乳人妻| 日本一本二区三区精品| 91午夜精品亚洲一区二区三区 | bbb黄色大片| 色视频www国产| 尤物成人国产欧美一区二区三区| 美女被艹到高潮喷水动态| 国产伦在线观看视频一区| 亚洲av二区三区四区| 老熟妇仑乱视频hdxx| ponron亚洲| 国产野战对白在线观看| 丰满的人妻完整版| 成年女人永久免费观看视频| 啪啪无遮挡十八禁网站| 最近在线观看免费完整版| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 亚洲第一欧美日韩一区二区三区| 久久久精品大字幕| 国产熟女xx| 亚洲三级黄色毛片| 久久久久免费精品人妻一区二区| 国产白丝娇喘喷水9色精品| 女同久久另类99精品国产91| 亚洲 欧美 日韩 在线 免费| 每晚都被弄得嗷嗷叫到高潮| 色在线成人网| 十八禁网站免费在线| 俄罗斯特黄特色一大片| 丝袜美腿在线中文| 美女cb高潮喷水在线观看| 亚洲人成电影免费在线| 可以在线观看的亚洲视频| 特大巨黑吊av在线直播| 国产久久久一区二区三区| 看免费av毛片| 伦理电影大哥的女人| 午夜福利在线观看免费完整高清在 | 亚洲av电影不卡..在线观看| 天堂av国产一区二区熟女人妻| 久久香蕉精品热| 一进一出好大好爽视频| 毛片女人毛片| 天堂网av新在线| 一进一出抽搐gif免费好疼| 亚洲无线在线观看| 桃色一区二区三区在线观看| 国产伦在线观看视频一区| 一区福利在线观看| 亚洲五月天丁香| 观看免费一级毛片| 三级国产精品欧美在线观看| 国产亚洲欧美在线一区二区| 国产视频内射| 亚洲va日本ⅴa欧美va伊人久久| 亚洲美女视频黄频| 午夜免费激情av| 免费在线观看日本一区| 此物有八面人人有两片| а√天堂www在线а√下载| 婷婷丁香在线五月| 久久国产精品影院| 亚洲国产欧洲综合997久久,| aaaaa片日本免费| 久久精品国产亚洲av涩爱 | 国产在线男女| 91av网一区二区| 成人无遮挡网站| 亚洲欧美日韩高清在线视频| 亚洲 国产 在线| 欧美性猛交╳xxx乱大交人| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久成人av| 啦啦啦观看免费观看视频高清| 欧美+日韩+精品| 精品久久久久久久久亚洲 | 搡老妇女老女人老熟妇| 中亚洲国语对白在线视频| 51国产日韩欧美| 免费黄网站久久成人精品 | 亚州av有码| 69人妻影院| 国产人妻一区二区三区在| 一级av片app| 久久久久久久久中文| 白带黄色成豆腐渣| 日韩欧美国产一区二区入口| 欧美黑人欧美精品刺激| 最近最新免费中文字幕在线| 女人被狂操c到高潮| 久久精品国产亚洲av涩爱 | 我要看日韩黄色一级片| 欧美xxxx性猛交bbbb| 国产成人啪精品午夜网站| 免费看光身美女| 国产乱人视频| 在线免费观看的www视频| 18禁在线播放成人免费| 午夜老司机福利剧场| 又爽又黄无遮挡网站| 脱女人内裤的视频| 婷婷亚洲欧美| 欧美+亚洲+日韩+国产| 国产精品免费一区二区三区在线| 亚洲五月天丁香| 天堂影院成人在线观看| 国语自产精品视频在线第100页| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 国产精华一区二区三区| 欧美国产日韩亚洲一区| 精品福利观看| 国产一区二区在线观看日韩| 99国产综合亚洲精品| 嫩草影院入口| 精品无人区乱码1区二区| 国产精品一区二区免费欧美| 一个人看视频在线观看www免费| 在线观看av片永久免费下载| 亚洲综合色惰| 亚洲精品影视一区二区三区av| 中出人妻视频一区二区| 午夜日韩欧美国产| 国产爱豆传媒在线观看| 亚洲av成人av| 亚洲成av人片免费观看| 久久久久久久亚洲中文字幕 | 国产高清有码在线观看视频| 久9热在线精品视频| 香蕉av资源在线| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 久久人妻av系列| 国内精品美女久久久久久| 国产主播在线观看一区二区| 日本a在线网址| x7x7x7水蜜桃| 亚洲成av人片在线播放无| 久久香蕉精品热| 九九久久精品国产亚洲av麻豆| 最近最新免费中文字幕在线| 别揉我奶头 嗯啊视频| 日本成人三级电影网站| 久久久久久久久大av| 免费观看的影片在线观看| 久久久精品大字幕| 国产一区二区三区视频了| 丰满人妻一区二区三区视频av| 日韩高清综合在线| 国产色爽女视频免费观看| 永久网站在线| 99热只有精品国产| 午夜老司机福利剧场| 蜜桃久久精品国产亚洲av| 好看av亚洲va欧美ⅴa在| 美女免费视频网站| 亚洲欧美日韩卡通动漫| 国产成人欧美在线观看| 在线国产一区二区在线| 又粗又爽又猛毛片免费看| 久久人人爽人人爽人人片va | 精品人妻一区二区三区麻豆 | 女人十人毛片免费观看3o分钟| 国产成人av教育| 中文字幕久久专区| ponron亚洲| 久久久久久久午夜电影| 毛片女人毛片| 真人做人爱边吃奶动态| 亚洲 国产 在线| 中文字幕熟女人妻在线| 在线免费观看不下载黄p国产 | 亚洲天堂国产精品一区在线| 精品不卡国产一区二区三区| 91麻豆精品激情在线观看国产| 最近最新免费中文字幕在线| 亚洲精品乱码久久久v下载方式| 黄色配什么色好看| 午夜免费成人在线视频| 久久久精品欧美日韩精品| 一级黄片播放器| 亚洲中文日韩欧美视频| 自拍偷自拍亚洲精品老妇| 免费看美女性在线毛片视频| 欧美黄色淫秽网站| 757午夜福利合集在线观看| 久久这里只有精品中国| 制服丝袜大香蕉在线| 最新在线观看一区二区三区| 搞女人的毛片| 黄色女人牲交| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 毛片一级片免费看久久久久 | 精品不卡国产一区二区三区| 观看美女的网站| 国产伦精品一区二区三区四那| 91麻豆精品激情在线观看国产| 亚洲第一欧美日韩一区二区三区| 91午夜精品亚洲一区二区三区 | 亚洲欧美日韩高清在线视频| 久久99热这里只有精品18| 一级av片app| 亚洲精品亚洲一区二区| 国产精品嫩草影院av在线观看 | 成人av一区二区三区在线看| 日韩欧美精品免费久久 | 十八禁人妻一区二区| 黄色一级大片看看| av在线天堂中文字幕| 成人一区二区视频在线观看| www.色视频.com| 欧美成人免费av一区二区三区| 亚洲av成人不卡在线观看播放网| 久久久色成人| av福利片在线观看| 岛国在线免费视频观看| 日韩欧美在线乱码| 国产国拍精品亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 精品久久久久久,| 在线免费观看不下载黄p国产 | 久久99热这里只有精品18| 午夜免费激情av| 国产精品国产高清国产av| 免费电影在线观看免费观看| 亚洲黑人精品在线| 国产精品野战在线观看| 2021天堂中文幕一二区在线观| 欧美在线一区亚洲| 精品久久久久久久久久久久久| 国产亚洲精品久久久久久毛片| 国产亚洲欧美在线一区二区| 成熟少妇高潮喷水视频| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 日韩中字成人| or卡值多少钱| www.熟女人妻精品国产| 日本成人三级电影网站| 丰满的人妻完整版| 午夜福利视频1000在线观看| 国产色婷婷99| 久久热精品热| 免费观看人在逋| 亚洲av不卡在线观看| 欧美激情久久久久久爽电影| 五月伊人婷婷丁香| 日韩欧美精品免费久久 | 久久这里只有精品中国| 欧美成人a在线观看| 欧美潮喷喷水| 身体一侧抽搐| 美女xxoo啪啪120秒动态图 | 日本精品一区二区三区蜜桃| 中文在线观看免费www的网站| 在线观看一区二区三区| 国产v大片淫在线免费观看| 国产野战对白在线观看| 国产日本99.免费观看| 午夜福利高清视频| 久9热在线精品视频| 身体一侧抽搐| 久久午夜亚洲精品久久| 黄色配什么色好看| 亚洲七黄色美女视频| 搞女人的毛片| 国产黄色小视频在线观看| 18+在线观看网站| 成人特级av手机在线观看| 免费一级毛片在线播放高清视频| 床上黄色一级片| 国产色婷婷99| 一二三四社区在线视频社区8| 亚洲在线自拍视频| 久久精品国产清高在天天线| 精品人妻一区二区三区麻豆 | 亚洲电影在线观看av| 18美女黄网站色大片免费观看| 99久久九九国产精品国产免费| 精品免费久久久久久久清纯| 久久久久免费精品人妻一区二区| 香蕉av资源在线| 国产成人av教育| 伊人久久精品亚洲午夜| www.999成人在线观看| 男人狂女人下面高潮的视频| 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 级片在线观看| 国产午夜福利久久久久久| 在线观看美女被高潮喷水网站 | 身体一侧抽搐| 两个人视频免费观看高清| 亚洲三级黄色毛片| 日韩欧美精品v在线| 日韩欧美精品免费久久 | 国产精品日韩av在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 永久网站在线| 亚洲av一区综合| 日本一本二区三区精品| 亚洲成人中文字幕在线播放| 亚洲五月婷婷丁香| netflix在线观看网站| 精品午夜福利在线看| 成年女人毛片免费观看观看9| 精品一区二区三区人妻视频| 在线观看美女被高潮喷水网站 | 亚洲精品亚洲一区二区| 99riav亚洲国产免费| 欧美精品啪啪一区二区三区| 九色成人免费人妻av| 日本在线视频免费播放| 日本与韩国留学比较| 日韩欧美精品免费久久 | 男女之事视频高清在线观看| 亚洲最大成人中文| 国产高清有码在线观看视频| 国产乱人伦免费视频| 国模一区二区三区四区视频| 精品不卡国产一区二区三区| 青草久久国产| 日韩欧美 国产精品| 热99re8久久精品国产| 日韩精品青青久久久久久| 1024手机看黄色片| 亚洲专区国产一区二区| 久久久久久久精品吃奶| 老司机深夜福利视频在线观看| 女人十人毛片免费观看3o分钟| 国产黄片美女视频| 中国美女看黄片| 久久香蕉精品热| 国产精品嫩草影院av在线观看 | eeuss影院久久| 日本三级黄在线观看| 日韩欧美一区二区三区在线观看| 国产在视频线在精品| 亚洲自偷自拍三级| 脱女人内裤的视频| 亚洲最大成人av| 久久精品国产亚洲av涩爱 | 成人av一区二区三区在线看| 日韩 亚洲 欧美在线| 午夜免费成人在线视频| 如何舔出高潮| 亚洲综合色惰| 中文字幕久久专区| or卡值多少钱| 757午夜福利合集在线观看| 亚洲五月天丁香| 精品一区二区三区视频在线观看免费| 可以在线观看的亚洲视频| av中文乱码字幕在线| 首页视频小说图片口味搜索| 午夜福利欧美成人| 欧美国产日韩亚洲一区| 五月玫瑰六月丁香| 欧美最黄视频在线播放免费| 欧美一区二区国产精品久久精品| 亚洲,欧美,日韩| 夜夜爽天天搞| 亚洲中文字幕一区二区三区有码在线看| 好男人在线观看高清免费视频| 99久久久亚洲精品蜜臀av| 国产日本99.免费观看| 婷婷色综合大香蕉| 中文在线观看免费www的网站| АⅤ资源中文在线天堂| 日日摸夜夜添夜夜添av毛片 | 亚洲av不卡在线观看| 精品久久久久久久久亚洲 | 国产色爽女视频免费观看| 日韩成人在线观看一区二区三区| av专区在线播放| 91在线观看av| 久久人人爽人人爽人人片va | 国产亚洲精品久久久久久毛片| 啪啪无遮挡十八禁网站| 熟妇人妻久久中文字幕3abv| 90打野战视频偷拍视频| 成人一区二区视频在线观看| 欧美黑人欧美精品刺激| 国产高清三级在线| 国产91精品成人一区二区三区| 蜜桃久久精品国产亚洲av| 五月伊人婷婷丁香| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 老鸭窝网址在线观看| 婷婷亚洲欧美| 99久久无色码亚洲精品果冻| 亚洲欧美日韩东京热| 欧美日韩乱码在线| 国产精品99久久久久久久久| 久久中文看片网| 黄色丝袜av网址大全| 在线观看一区二区三区| 久久草成人影院| 自拍偷自拍亚洲精品老妇| 国产中年淑女户外野战色| 日韩人妻高清精品专区| 亚洲欧美日韩高清专用| 精品不卡国产一区二区三区| 午夜福利在线观看吧| 亚洲经典国产精华液单 | 日日夜夜操网爽| 亚洲自偷自拍三级| 亚洲人成网站高清观看| 91麻豆精品激情在线观看国产| 欧美黄色淫秽网站| ponron亚洲| 久久久国产成人免费| 国产精品一区二区性色av| 国产精品嫩草影院av在线观看 | 国产精品不卡视频一区二区 | 亚洲av成人精品一区久久| www.色视频.com| 直男gayav资源| 搡女人真爽免费视频火全软件 | 亚洲美女黄片视频| 桃色一区二区三区在线观看| av中文乱码字幕在线| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 99热这里只有是精品在线观看 | 亚洲一区高清亚洲精品| 亚洲美女黄片视频| 国产成人欧美在线观看| av在线观看视频网站免费| 一本一本综合久久| 又紧又爽又黄一区二区| 国产乱人视频| 一级黄片播放器| 91久久精品电影网| 51国产日韩欧美| 国产av在哪里看| 国产单亲对白刺激| 亚洲成人中文字幕在线播放| 校园春色视频在线观看| 3wmmmm亚洲av在线观看| 中出人妻视频一区二区| 搡女人真爽免费视频火全软件 | 一个人观看的视频www高清免费观看| 悠悠久久av| 在线看三级毛片| 91麻豆精品激情在线观看国产| 国产亚洲精品综合一区在线观看| 97热精品久久久久久| 午夜福利在线在线| 一本精品99久久精品77| 九九热线精品视视频播放| 亚洲av二区三区四区| 国产成人av教育| 自拍偷自拍亚洲精品老妇| 国产 一区 欧美 日韩| 美女xxoo啪啪120秒动态图 | aaaaa片日本免费| 国产欧美日韩一区二区三| 亚洲午夜理论影院| 午夜免费激情av| 亚洲自拍偷在线| 精品久久久久久久久久免费视频| 51午夜福利影视在线观看| 麻豆久久精品国产亚洲av| 一级av片app| 免费高清视频大片| 黄色一级大片看看| 一进一出抽搐动态| 亚洲美女视频黄频| 国产亚洲精品综合一区在线观看| a在线观看视频网站| 色在线成人网| 中文字幕久久专区| 啪啪无遮挡十八禁网站| 久久欧美精品欧美久久欧美| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 高清在线国产一区| a级毛片a级免费在线| 亚洲精品一卡2卡三卡4卡5卡| 很黄的视频免费| 国产精品免费一区二区三区在线| 欧美乱色亚洲激情| 最近在线观看免费完整版| 在线观看一区二区三区| 色综合婷婷激情| 午夜免费激情av| 久久国产乱子免费精品| 窝窝影院91人妻| 亚洲国产精品成人综合色| 国产爱豆传媒在线观看| 黄色丝袜av网址大全| 国内毛片毛片毛片毛片毛片| 国产精品久久久久久久电影| 日本撒尿小便嘘嘘汇集6| 亚洲av免费在线观看| 熟妇人妻久久中文字幕3abv| 国产淫片久久久久久久久 | 一级av片app| 精品久久久久久久久av| 国内精品一区二区在线观看| 三级毛片av免费| 免费av毛片视频| 91九色精品人成在线观看| 真人一进一出gif抽搐免费| 国产主播在线观看一区二区| 韩国av一区二区三区四区| 日韩亚洲欧美综合| 国产主播在线观看一区二区| 在线观看舔阴道视频| 日韩亚洲欧美综合| 国产伦精品一区二区三区视频9| 中文资源天堂在线| 国产欧美日韩一区二区精品| 久久婷婷人人爽人人干人人爱| 国产午夜精品论理片| 午夜精品久久久久久毛片777| 久久久成人免费电影| 熟妇人妻久久中文字幕3abv| 亚洲人成电影免费在线| 91麻豆精品激情在线观看国产| 免费看光身美女| 少妇人妻精品综合一区二区 | 欧洲精品卡2卡3卡4卡5卡区| 九九热线精品视视频播放| 成年版毛片免费区| 成年免费大片在线观看| 天堂动漫精品| 免费av不卡在线播放| 男人狂女人下面高潮的视频| 久久久久久久午夜电影| 动漫黄色视频在线观看|