• <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)展
    999精品在线视频| 欧美zozozo另类| 免费在线观看黄色视频的| 日本在线视频免费播放| 老熟妇乱子伦视频在线观看| 男人舔女人下体高潮全视频| 最近最新中文字幕大全免费视频| 精品久久蜜臀av无| 99国产精品一区二区三区| 老司机靠b影院| 宅男免费午夜| 在线观看午夜福利视频| 91大片在线观看| 手机成人av网站| 日韩成人在线观看一区二区三区| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 91麻豆精品激情在线观看国产| 欧美在线一区亚洲| www.精华液| 国产熟女xx| 此物有八面人人有两片| 美女高潮到喷水免费观看| 给我免费播放毛片高清在线观看| 精品一区二区三区av网在线观看| 成人特级黄色片久久久久久久| 久久精品aⅴ一区二区三区四区| 一夜夜www| xxxwww97欧美| 一本大道久久a久久精品| 神马国产精品三级电影在线观看 | 午夜精品久久久久久毛片777| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 黄片小视频在线播放| 又大又爽又粗| 中文字幕久久专区| 日本免费一区二区三区高清不卡| 一个人免费在线观看的高清视频| 久久久久久亚洲精品国产蜜桃av| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 十八禁人妻一区二区| 久久伊人香网站| 国产黄片美女视频| 国产不卡一卡二| 亚洲精品中文字幕一二三四区| 成人国语在线视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲真实伦在线观看| 90打野战视频偷拍视频| 很黄的视频免费| 波多野结衣巨乳人妻| 国产爱豆传媒在线观看 | 给我免费播放毛片高清在线观看| 欧美最黄视频在线播放免费| 最好的美女福利视频网| 成人午夜高清在线视频 | 88av欧美| 欧美日韩福利视频一区二区| 在线观看舔阴道视频| 亚洲真实伦在线观看| 99国产综合亚洲精品| 日韩一卡2卡3卡4卡2021年| 亚洲av成人一区二区三| 精品久久久久久,| 免费看美女性在线毛片视频| 美女扒开内裤让男人捅视频| 成人国产综合亚洲| 一区二区日韩欧美中文字幕| 亚洲成av人片免费观看| 真人做人爱边吃奶动态| 成人三级黄色视频| 男女之事视频高清在线观看| 高清在线国产一区| 成人av一区二区三区在线看| 日本免费a在线| 成年女人毛片免费观看观看9| 亚洲自偷自拍图片 自拍| 日韩成人在线观看一区二区三区| 免费看十八禁软件| 看片在线看免费视频| www.自偷自拍.com| 亚洲欧美一区二区三区黑人| www.熟女人妻精品国产| 丝袜在线中文字幕| 黑丝袜美女国产一区| 在线av久久热| 亚洲欧美精品综合一区二区三区| 岛国视频午夜一区免费看| 91国产中文字幕| 成熟少妇高潮喷水视频| 18禁黄网站禁片午夜丰满| 久久久久久国产a免费观看| 国产极品粉嫩免费观看在线| 看片在线看免费视频| 丁香欧美五月| 69av精品久久久久久| 亚洲avbb在线观看| 欧美日韩黄片免| 久久久水蜜桃国产精品网| 日本一区二区免费在线视频| 久久久久久久久中文| 女人高潮潮喷娇喘18禁视频| 亚洲av第一区精品v没综合| 午夜福利成人在线免费观看| 999久久久精品免费观看国产| 免费高清视频大片| aaaaa片日本免费| 搞女人的毛片| 国产91精品成人一区二区三区| 丝袜人妻中文字幕| 日日夜夜操网爽| www.自偷自拍.com| 日韩欧美免费精品| 高清在线国产一区| 男人舔奶头视频| 精品免费久久久久久久清纯| 欧美大码av| 夜夜看夜夜爽夜夜摸| 久久久久久久久久黄片| 欧美日韩一级在线毛片| 欧美乱色亚洲激情| 两个人免费观看高清视频| 国产色视频综合| 久久国产亚洲av麻豆专区| 人人妻人人澡欧美一区二区| 一夜夜www| svipshipincom国产片| 欧美三级亚洲精品| 久久中文字幕人妻熟女| 国产高清videossex| 999久久久精品免费观看国产| 视频区欧美日本亚洲| 夜夜夜夜夜久久久久| 欧美激情久久久久久爽电影| 一二三四社区在线视频社区8| 久久这里只有精品19| 国产精品乱码一区二三区的特点| www.999成人在线观看| 久久这里只有精品19| 国产精品乱码一区二三区的特点| 亚洲精品美女久久久久99蜜臀| 三级毛片av免费| 搞女人的毛片| 99在线人妻在线中文字幕| 99久久精品国产亚洲精品| 在线十欧美十亚洲十日本专区| 免费av毛片视频| 亚洲色图 男人天堂 中文字幕| 人人妻人人看人人澡| 久久国产乱子伦精品免费另类| 日本 欧美在线| 久久青草综合色| 国产三级在线视频| 男人的好看免费观看在线视频 | 一级黄色大片毛片| 听说在线观看完整版免费高清| 欧美日韩亚洲国产一区二区在线观看| 妹子高潮喷水视频| 欧美在线一区亚洲| 国产成人精品久久二区二区免费| 国产在线观看jvid| 后天国语完整版免费观看| 91在线观看av| 日本精品一区二区三区蜜桃| 午夜福利成人在线免费观看| 18禁观看日本| 欧美成人一区二区免费高清观看 | 午夜久久久久精精品| 亚洲成人久久性| 俄罗斯特黄特色一大片| 久久伊人香网站| 国产爱豆传媒在线观看 | 校园春色视频在线观看| 韩国精品一区二区三区| 欧美黑人巨大hd| 亚洲精品av麻豆狂野| 亚洲国产精品999在线| 国产精品 国内视频| 窝窝影院91人妻| www日本黄色视频网| 亚洲五月婷婷丁香| 亚洲av成人一区二区三| 亚洲成av片中文字幕在线观看| 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区mp4| 亚洲精品粉嫩美女一区| 国产精品永久免费网站| 午夜福利高清视频| 1024视频免费在线观看| 国产伦一二天堂av在线观看| 国产精品一区二区三区四区久久 | 一进一出抽搐gif免费好疼| 国产伦人伦偷精品视频| 中文字幕久久专区| 日韩欧美一区二区三区在线观看| 操出白浆在线播放| 日韩中文字幕欧美一区二区| 亚洲七黄色美女视频| xxx96com| 久久精品夜夜夜夜夜久久蜜豆 | 久久草成人影院| 久久香蕉精品热| av在线天堂中文字幕| 国产人伦9x9x在线观看| 亚洲欧美日韩无卡精品| 欧美成狂野欧美在线观看| 免费在线观看视频国产中文字幕亚洲| 免费高清视频大片| √禁漫天堂资源中文www| 91字幕亚洲| 色播亚洲综合网| 国产激情偷乱视频一区二区| 亚洲第一青青草原| 性欧美人与动物交配| 国产99白浆流出| 18禁美女被吸乳视频| 欧美最黄视频在线播放免费| 久久精品国产清高在天天线| 久久久国产成人精品二区| 一区二区三区激情视频| 69av精品久久久久久| 免费一级毛片在线播放高清视频| 欧美激情 高清一区二区三区| 国产区一区二久久| 老熟妇仑乱视频hdxx| 69av精品久久久久久| 看片在线看免费视频| 草草在线视频免费看| x7x7x7水蜜桃| 一级黄色大片毛片| 人成视频在线观看免费观看| 国产精品电影一区二区三区| 免费在线观看完整版高清| 亚洲一区高清亚洲精品| 91字幕亚洲| 久久精品成人免费网站| 黑人操中国人逼视频| 久久精品国产清高在天天线| 亚洲国产中文字幕在线视频| 成人国产一区最新在线观看| 国产色视频综合| 岛国在线观看网站| 欧美激情高清一区二区三区| 国产99久久九九免费精品| 午夜福利欧美成人| 国产视频内射| 午夜福利高清视频| 免费高清在线观看日韩| 精品一区二区三区av网在线观看| 美女午夜性视频免费| 日本一区二区免费在线视频| 老汉色∧v一级毛片| 欧美黑人巨大hd| 91九色精品人成在线观看| 亚洲久久久国产精品| 午夜久久久在线观看| 男人舔奶头视频| 国产精华一区二区三区| www日本在线高清视频| 女生性感内裤真人,穿戴方法视频| 黄色丝袜av网址大全| 在线视频色国产色| 99久久精品国产亚洲精品| 国产99白浆流出| 久久亚洲真实| 日本熟妇午夜| 叶爱在线成人免费视频播放| 午夜免费鲁丝| 国产不卡一卡二| 在线看三级毛片| 色哟哟哟哟哟哟| 免费搜索国产男女视频| 一本大道久久a久久精品| 国产午夜精品久久久久久| 亚洲第一电影网av| 国产在线观看jvid| 久久久国产成人精品二区| 日韩免费av在线播放| www.www免费av| 日韩国内少妇激情av| 成人免费观看视频高清| 一级毛片精品| 免费高清视频大片| 久久热在线av| 久久天堂一区二区三区四区| 国产一区在线观看成人免费| 男人的好看免费观看在线视频 | av在线播放免费不卡| 日韩一卡2卡3卡4卡2021年| 12—13女人毛片做爰片一| 日本 欧美在线| 一级毛片高清免费大全| 亚洲精品美女久久久久99蜜臀| 国产av一区在线观看免费| 日韩成人在线观看一区二区三区| 国产区一区二久久| 国产亚洲精品第一综合不卡| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 麻豆成人午夜福利视频| 亚洲一区二区三区色噜噜| 亚洲五月天丁香| 无限看片的www在线观看| 午夜福利免费观看在线| 看免费av毛片| 久久久水蜜桃国产精品网| 久久国产精品人妻蜜桃| 精品国产美女av久久久久小说| www国产在线视频色| 99国产极品粉嫩在线观看| 好男人在线观看高清免费视频 | 人人妻人人澡人人看| 精品国产亚洲在线| 欧美乱妇无乱码| 一二三四在线观看免费中文在| 免费女性裸体啪啪无遮挡网站| 9191精品国产免费久久| 黄色视频,在线免费观看| 国内精品久久久久精免费| 他把我摸到了高潮在线观看| 级片在线观看| 欧美午夜高清在线| 亚洲欧洲精品一区二区精品久久久| 18禁黄网站禁片免费观看直播| 精品一区二区三区四区五区乱码| 伊人久久大香线蕉亚洲五| 国产真实乱freesex| 久久精品91蜜桃| 国产成人精品久久二区二区91| 亚洲在线自拍视频| 2021天堂中文幕一二区在线观 | 99久久无色码亚洲精品果冻| 亚洲人成网站高清观看| 一区福利在线观看| 欧美色欧美亚洲另类二区| 一边摸一边做爽爽视频免费| 黄频高清免费视频| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜一区二区| 欧美一级毛片孕妇| 97碰自拍视频| 午夜精品在线福利| 国产精品亚洲美女久久久| 亚洲av成人av| 亚洲第一电影网av| 欧美在线黄色| 国内精品久久久久久久电影| 亚洲人成网站在线播放欧美日韩| 久久久久久久精品吃奶| 欧美乱色亚洲激情| 淫秽高清视频在线观看| 国产精品永久免费网站| 美女午夜性视频免费| 日本一区二区免费在线视频| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久人妻精品电影| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 久久九九热精品免费| 欧美精品亚洲一区二区| 国产精品自产拍在线观看55亚洲| 高潮久久久久久久久久久不卡| 国产精品久久久人人做人人爽| 欧美丝袜亚洲另类 | 99国产精品一区二区蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 美女大奶头视频| 一本久久中文字幕| 亚洲一区二区三区不卡视频| 亚洲男人天堂网一区| 校园春色视频在线观看| av中文乱码字幕在线| 精品欧美国产一区二区三| 777久久人妻少妇嫩草av网站| 亚洲九九香蕉| 一级毛片精品| 欧美午夜高清在线| av片东京热男人的天堂| 丝袜在线中文字幕| 久久天躁狠狠躁夜夜2o2o| 婷婷亚洲欧美| 天天添夜夜摸| 一级毛片精品| 天堂√8在线中文| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品第一综合不卡| 午夜精品久久久久久毛片777| 国产精品乱码一区二三区的特点| 亚洲狠狠婷婷综合久久图片| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av香蕉五月| 精品久久蜜臀av无| 此物有八面人人有两片| 久久性视频一级片| 曰老女人黄片| 天天添夜夜摸| 精品久久久久久久人妻蜜臀av| www.999成人在线观看| 午夜久久久在线观看| 色在线成人网| 在线观看日韩欧美| 亚洲专区国产一区二区| 丝袜美腿诱惑在线| 国产黄色小视频在线观看| 欧美成人免费av一区二区三区| 欧美亚洲日本最大视频资源| 久久久久九九精品影院| 亚洲一区二区三区色噜噜| 欧美日韩福利视频一区二区| 国产野战对白在线观看| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 精品高清国产在线一区| 少妇裸体淫交视频免费看高清 | 亚洲精品一卡2卡三卡4卡5卡| 日韩高清综合在线| 老司机靠b影院| 午夜激情福利司机影院| 国产在线观看jvid| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 男男h啪啪无遮挡| 欧美性猛交╳xxx乱大交人| 亚洲av成人一区二区三| 亚洲aⅴ乱码一区二区在线播放 | 亚洲,欧美精品.| 夜夜躁狠狠躁天天躁| 国产亚洲欧美在线一区二区| 国产精品亚洲美女久久久| 欧美在线一区亚洲| 国产午夜福利久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播| 成人国语在线视频| 色综合亚洲欧美另类图片| 亚洲欧美日韩高清在线视频| 国产精品精品国产色婷婷| 欧美绝顶高潮抽搐喷水| 中国美女看黄片| 麻豆久久精品国产亚洲av| x7x7x7水蜜桃| 中文字幕人成人乱码亚洲影| 最近最新中文字幕大全电影3 | 两个人看的免费小视频| 久久久久久久午夜电影| 久久精品亚洲精品国产色婷小说| 日本免费a在线| 亚洲精品一卡2卡三卡4卡5卡| 嫁个100分男人电影在线观看| 国产精品 国内视频| 婷婷精品国产亚洲av| 日韩大码丰满熟妇| 亚洲成人久久性| 国产区一区二久久| 久久99热这里只有精品18| 啦啦啦观看免费观看视频高清| 亚洲狠狠婷婷综合久久图片| 国语自产精品视频在线第100页| a在线观看视频网站| 12—13女人毛片做爰片一| 亚洲免费av在线视频| 日韩欧美一区视频在线观看| 日日夜夜操网爽| 免费女性裸体啪啪无遮挡网站| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片午夜丰满| 美女大奶头视频| 国产aⅴ精品一区二区三区波| 欧美激情极品国产一区二区三区| 日本熟妇午夜| 欧美乱妇无乱码| 日本五十路高清| 成年女人毛片免费观看观看9| 香蕉av资源在线| 波多野结衣巨乳人妻| 国产成人精品久久二区二区91| 欧美三级亚洲精品| 欧美黄色淫秽网站| 香蕉av资源在线| 一级毛片高清免费大全| 色播亚洲综合网| 欧美亚洲日本最大视频资源| 久久精品国产99精品国产亚洲性色| 男女做爰动态图高潮gif福利片| 久久青草综合色| 日本免费a在线| 中文亚洲av片在线观看爽| 国产成人一区二区三区免费视频网站| 亚洲片人在线观看| 免费观看精品视频网站| 亚洲欧美一区二区三区黑人| 999久久久精品免费观看国产| 国产视频内射| 亚洲成人国产一区在线观看| 亚洲国产日韩欧美精品在线观看 | 999久久久国产精品视频| www国产在线视频色| 黄片小视频在线播放| 亚洲午夜理论影院| 日本 欧美在线| 日韩欧美 国产精品| 老司机午夜福利在线观看视频| 久久精品国产清高在天天线| 91在线观看av| a在线观看视频网站| 99国产极品粉嫩在线观看| 日本在线视频免费播放| 欧美日韩中文字幕国产精品一区二区三区| 女性被躁到高潮视频| 国产亚洲av嫩草精品影院| 国产野战对白在线观看| 久久国产精品人妻蜜桃| 久久人妻福利社区极品人妻图片| 女同久久另类99精品国产91| 97超级碰碰碰精品色视频在线观看| 色综合亚洲欧美另类图片| 看黄色毛片网站| 老鸭窝网址在线观看| 亚洲全国av大片| 欧美中文综合在线视频| 欧美激情极品国产一区二区三区| 亚洲av五月六月丁香网| 欧美日韩福利视频一区二区| 日本 av在线| 窝窝影院91人妻| 亚洲黑人精品在线| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| 91成年电影在线观看| 国语自产精品视频在线第100页| 一进一出抽搐动态| 国产精品久久久人人做人人爽| 久久久久久久久中文| 麻豆av在线久日| 亚洲精品在线观看二区| 观看免费一级毛片| 黄频高清免费视频| 亚洲成人久久性| 亚洲精品国产一区二区精华液| 女生性感内裤真人,穿戴方法视频| 白带黄色成豆腐渣| 中文字幕精品亚洲无线码一区 | 特大巨黑吊av在线直播 | 亚洲最大成人中文| 精品熟女少妇八av免费久了| 国内揄拍国产精品人妻在线 | www.999成人在线观看| 少妇熟女aⅴ在线视频| 国产精华一区二区三区| 美女免费视频网站| 啦啦啦观看免费观看视频高清| 国产成人精品无人区| 成人av一区二区三区在线看| 午夜成年电影在线免费观看| 亚洲人成77777在线视频| 麻豆成人av在线观看| 两人在一起打扑克的视频| 免费电影在线观看免费观看| 香蕉久久夜色| 国产亚洲av嫩草精品影院| 啦啦啦免费观看视频1| 久久 成人 亚洲| 一级片免费观看大全| 神马国产精品三级电影在线观看 | 老汉色∧v一级毛片| 熟女少妇亚洲综合色aaa.| 90打野战视频偷拍视频| 亚洲国产高清在线一区二区三 | 欧美+亚洲+日韩+国产| 国内久久婷婷六月综合欲色啪| or卡值多少钱| 黄网站色视频无遮挡免费观看| 成年女人毛片免费观看观看9| 日韩三级视频一区二区三区| 成年版毛片免费区| xxxwww97欧美| 久久精品影院6| 婷婷精品国产亚洲av| 国产区一区二久久| 国产av不卡久久| 国产精品1区2区在线观看.| 免费在线观看亚洲国产| 亚洲,欧美精品.| 久久久久国产精品人妻aⅴ院| 三级毛片av免费| 久久久久久久久久黄片| 婷婷亚洲欧美| xxxwww97欧美| 久久精品影院6| 在线观看66精品国产| e午夜精品久久久久久久| 天堂√8在线中文| 亚洲精品一卡2卡三卡4卡5卡| 成年免费大片在线观看| 日韩三级视频一区二区三区| 精品少妇一区二区三区视频日本电影| 欧美中文日本在线观看视频| 美女高潮到喷水免费观看| 亚洲成人免费电影在线观看| 成熟少妇高潮喷水视频| 亚洲狠狠婷婷综合久久图片| 亚洲一区二区三区色噜噜| 欧美国产日韩亚洲一区| 久久人妻福利社区极品人妻图片| 午夜视频精品福利| 国产又黄又爽又无遮挡在线| 国产成人精品久久二区二区91| 久久久久亚洲av毛片大全| 日韩大尺度精品在线看网址|