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

    基于ISSR數(shù)據(jù)探討節(jié)莖曲柄蘚(Campylopus umbellatus (Arn.) Paris)遺傳多樣性和居群遺傳結(jié)構(gòu)

    2019-06-10 05:13單壹壹申琳李丹丹于晶
    關(guān)鍵詞:嵊泗列島遺傳變異

    單壹壹 申琳 李丹丹 于晶

    摘 要: 選用簡(jiǎn)單序列重復(fù)區(qū)間擴(kuò)增多態(tài)性(ISSR)分子標(biāo)記,對(duì)分布于浙江清涼峰國(guó)家級(jí)自然保護(hù)區(qū)順溪塢和涼源、浙江臺(tái)州臨海天臺(tái)山、浙江舟山嵊泗列島、江蘇蘇州天平山、廣西上思十萬(wàn)大山、臺(tái)灣宜蘭縣馬告生態(tài)公園和云南西雙版納納板河流域國(guó)家自然保護(hù)區(qū)8個(gè)自然居群的40份節(jié)莖曲柄蘚(Campylopus umbellatus (Arn.) Par.)的遺傳多樣性進(jìn)行了研究.應(yīng)用篩選出的12條引物,共擴(kuò)增出172條清晰、重復(fù)性高的條帶.使用POPGENE 3.2分析發(fā)現(xiàn):其中158條條帶的多態(tài)性位點(diǎn)百分率(PPB)為91.86%,根井正利基因多樣性指數(shù)(H)為0.2922,香農(nóng)-維納多樣性指數(shù)(I)為0.4459,說(shuō)明節(jié)莖曲柄蘚居群遺傳多樣性水平較高.8個(gè)居群的遺傳分化系數(shù)(Gst)為0.6346,居群間的基因流(Nm)為0.2879,63.33%的遺傳變異存在于居群間,36.67%的遺傳變異存在于居群內(nèi),說(shuō)明節(jié)莖曲柄蘚居群間的遺傳分化明顯.基于ISSR數(shù)據(jù)的聚類(lèi)分析表明:以遺傳距離65為分組閾值,8個(gè)居群可被分為6組,節(jié)莖曲柄蘚居群間的遺傳分化主要由地理因素造成,居群內(nèi)的遺傳分化可能與生境的異質(zhì)性有關(guān).

    關(guān)鍵詞: 節(jié)莖曲柄蘚; 簡(jiǎn)單序列重復(fù)區(qū)間擴(kuò)增多態(tài)性(ISSR); 遺傳多樣性

    中圖分類(lèi)號(hào): Q 948文獻(xiàn)標(biāo)志碼: A文章編號(hào): 1000-5137(2019)02-0207-10

    0 引 言

    在分子水平上,對(duì)于苔蘚植物的遺傳研究一般運(yùn)用分子標(biāo)記法[1].目前,DNA分子標(biāo)記技術(shù)已發(fā)展出10余種,常用手段包括隨機(jī)擴(kuò)增多態(tài)性DNA標(biāo)記(RAPD)、簡(jiǎn)單序列重復(fù)區(qū)間擴(kuò)增多態(tài)性(ISSR)、擴(kuò)增片段長(zhǎng)度多態(tài)性(AFLP)、相關(guān)序列擴(kuò)增多態(tài)性(SRAP)等[1-2].ISSR是由加拿大ZIETKIEWICZ等[3]提出的,以錨定的微衛(wèi)星為引物,擴(kuò)增間隔較小的兩個(gè)簡(jiǎn)單序列長(zhǎng)度多態(tài)性(SSR)間的序列,再通過(guò)凝膠電泳條帶分析樣本DNA多態(tài)性的方法[2].其優(yōu)點(diǎn)是無(wú)需預(yù)先知道基因組的任何信息就可以進(jìn)行標(biāo)記,可以在生長(zhǎng)周期的任何階段進(jìn)行檢測(cè),而且引物設(shè)計(jì)容易,成本低,DNA用量少[1-2,4],被認(rèn)為是一種非常理想的分子標(biāo)記方法.ISSR分子標(biāo)記技術(shù)已被廣泛應(yīng)用于許多被子植物[2,5-7]的遺傳多樣性分析中,但對(duì)于苔蘚植物,如:黃絲瓜蘚(Pohlia nutans (Hedw.) Lindb.)[8]、鼠尾蘚(Myuroclada maximowiczii (G.G.Borshch.) Steere & W.B.Schofield)[9]、小羽蘚屬(Haplocladium (Müll.Hal.) Müll.Hal.)[10]和中華蓑蘚(Macromitrium cavaleriei Cardot & Thér.)[11]的遺傳多樣性的研究,常使用的是RAPD標(biāo)記法.汪琛穎等[12]對(duì)狹邊大葉蘚(Rhodobryum ontariense (Kindb.) Paris)[4]和真蘚科(Bryaceae Schwgr.),王瑩瑩[13]對(duì)浙江千島湖大灰蘚(Hypnum plumiforme Wilson)和東亞小金發(fā)蘚(Pogonatum inflexum (Lindb.) Sande Lac.),買(mǎi)買(mǎi)提明·蘇來(lái)曼等[14]對(duì)水蘚科(Fontinalaceae Schimp.),魏青永等[15]對(duì)卷葉鳳尾蘚(Fissidens dubius P.Beauv.)使用ISSR標(biāo)記法進(jìn)行了遺傳多樣性的研究.

    曲尾蘚科(Dicranaceae)的曲柄蘚屬(Campylopus Brid.)是泛熱帶分布的類(lèi)群[16].自BRIDEL[16]基于Campylopus flexosus (Hedw.) Brid.創(chuàng)立曲柄蘚屬以來(lái),該屬植物先后被報(bào)道多達(dá)1480 種,目前接受的種名有160種左右[17].

    中國(guó)曲柄蘚屬植物目前被接受的有20余種[18-21].其中,節(jié)莖曲柄蘚(Campylopus umbellatus (Arn.) Paris)分布于中國(guó)十幾個(gè)省.該種植物體形態(tài)、大小以及不同部位葉片特征變化很大,毛尖有無(wú)、毛尖長(zhǎng)度、莖長(zhǎng)度以及葉片和中肋寬度的比值等都會(huì)有強(qiáng)烈變異[21].其典型特點(diǎn)是:植物體中等至大型,粗壯,墨綠色或黃綠色,常具透明長(zhǎng)毛尖,葉基部寬,具明顯角細(xì)胞,不突出或略突出,中肋占葉片基部1/3~1/2,中肋背部有2~4個(gè)細(xì)胞高的櫛片,葉片中上部細(xì)胞為紡錘形至短蟲(chóng)形.BROTHERUS[22]把葉片無(wú)透明毛尖和透明毛尖短的類(lèi)型歸類(lèi)為C.sinensis (Müll.Hal.) J.-P.Frahm;DUBY[23]把葉片有毛尖、透明毛尖短、中肋上部櫛片低的類(lèi)型歸類(lèi)為C.nigrescens (Mitt.) Duby.DIXION[23]根據(jù)該種植物體與C.coreensis Cardot相似,因其葉片稍長(zhǎng)、葉寬、中肋背面櫛片稍高等,將其歸類(lèi)為一個(gè)變種C.coreensis Card.var.amoyensis Dixon & Thér;高謙[21]認(rèn)為以上各類(lèi)型均屬本種特點(diǎn)范疇,故都擬為本種異名,并且認(rèn)為本種并非中美洲和南美洲產(chǎn)的C.richardii Brid..FRAHM[14]在對(duì)東亞曲柄蘚屬植物進(jìn)行修訂時(shí),列出了該種的12個(gè)異名,其中包括了C.coreensis Cardot.由密蘇里植物園主辦的Tropics(http://www.tropicos.org/)網(wǎng)站收錄的該種異名達(dá)38個(gè).

    FRAHM[16]認(rèn)為C.umbellatus (Arn.) Paris植物的大小在其分布范圍內(nèi),從內(nèi)部到邊緣呈梯度變化,在東南亞和中國(guó)南部高達(dá)幾厘米,但在北部?jī)H1 cm左右.為了解這些形態(tài)變異是否有遺傳基礎(chǔ),本文作者基于ISSR技術(shù),對(duì)該種的8個(gè)地理居群進(jìn)行了遺傳多樣性分析.

    1 實(shí)驗(yàn)材料與方法

    1.1 實(shí)驗(yàn)材料

    本實(shí)驗(yàn)對(duì)象為采自浙江、江蘇、廣西、云南和臺(tái)灣的節(jié)莖曲柄蘚8個(gè)居群40份樣品,具體采集信息見(jiàn)表1.為確保采集的樣品是不同克隆的植株,每個(gè)居群采樣間距最小為5 m [14].樣品采集后裝于牛皮紙標(biāo)本袋中,自然風(fēng)干后冷凍保存?zhèn)溆?實(shí)驗(yàn)所涉及的憑證標(biāo)本保存于上海師范大學(xué)(SHNU)苔蘚植物標(biāo)本館.

    1.2 試劑

    植物基因組DNA提取試劑盒(天根生化科技(北京)有限公司,DP305);Es Taq MasterMix(康為世紀(jì)生物科技有限公司,CW0690);所有引物均由上海生物工程有限公司合成并提供.

    1.3 DNA提取及聚合酶鏈?zhǔn)椒磻?yīng)(PCR)擴(kuò)增

    選取40份節(jié)莖曲柄蘚的新鮮植物體莖尖,無(wú)菌水沖洗干凈并用變色硅膠粒干燥后,分別用試劑盒提取基因組DNA,制成DNA樣品,-20℃保存?zhèn)溆?從加拿大哥倫比亞大學(xué)(UBC)公布的100條ISSR通用引物中,參考文獻(xiàn)[2]選取退火溫度為50~55 ℃之間的56條引物進(jìn)行篩選,篩出的12條引物見(jiàn)表2;采用康為公司提供的Es Taq MasterMix 50 μL反應(yīng)體系進(jìn)行PCR擴(kuò)增,擴(kuò)增程序參考汪琛穎等[4]對(duì)狹邊大葉蘚遺傳多樣性ISSR分析的擴(kuò)增程序;PCR擴(kuò)增產(chǎn)物在含溴化乙錠(EB)核酸染料的質(zhì)量濃度為20 g·mL-1的瓊脂糖凝膠電泳中分離,電壓為80 V;應(yīng)用全自動(dòng)數(shù)碼凝膠成像分析系統(tǒng)(Tanon 2500)對(duì)電泳結(jié)果進(jìn)行觀察并拍照記錄.

    1.4 數(shù)據(jù)記錄與分析

    將實(shí)驗(yàn)得到的凝膠電泳圖中的每一個(gè)條帶均視為一個(gè)分子標(biāo)記,統(tǒng)計(jì)電泳圖譜是否有處于同一位置的DNA帶,有的記為“1”,無(wú)帶的記為“0”,得到一個(gè)0/1二元數(shù)據(jù)矩陣[4].采用POPGENE 3.2分析軟件對(duì)節(jié)莖曲柄蘚所有居群進(jìn)行遺傳多樣性參數(shù)分析[4],分別計(jì)算多態(tài)性位點(diǎn)百分率(PPB)、根井正利基因多樣性指數(shù)(H)、香農(nóng)-維納多樣性指數(shù)(I)等遺傳多樣性參數(shù)以及各居群間的遺傳分化系數(shù)和基因流等遺傳分化指數(shù),并分別計(jì)算12條引物的PPB.運(yùn)用PCORD 4.0軟件,根據(jù)40份樣品的遺傳距離構(gòu)建反映各居群遺傳關(guān)系的聚類(lèi)圖,以研究各居群間遺傳變異與地理距離的相關(guān)性[2,4].

    2 結(jié)果與分析

    2.1 節(jié)莖曲柄蘚ISSR擴(kuò)增結(jié)果

    利用篩選出擴(kuò)增條帶清晰、重復(fù)性好的12條引物,對(duì)8個(gè)居群40份節(jié)莖曲柄蘚樣品進(jìn)行PCR擴(kuò)增,共得到172個(gè)擴(kuò)增位點(diǎn),其中158個(gè)為多態(tài)性位點(diǎn),PPB達(dá)91.86%,表明8個(gè)居群的40份節(jié)莖曲柄蘚樣品在分子水平上具有豐富的多態(tài)性,遺傳基礎(chǔ)較廣.12個(gè)引物擴(kuò)增出的條帶大小多集中在300~2500 bp之間.由表2可知,平均每個(gè)引物擴(kuò)增出位點(diǎn)14.33個(gè),多態(tài)性位點(diǎn)13.17個(gè).擴(kuò)增位點(diǎn)最多的是UBC841,共擴(kuò)增出19個(gè)位點(diǎn);引物UBC808,UBC880次之,分別為18個(gè);最少的是UBC825和UBC827,只有9個(gè).引物多態(tài)性位點(diǎn)百分率最高的是UBC825,UBC827和UBC842,可達(dá)100%.引物UBC826的PPB最低,但也達(dá)81.82%,說(shuō)明引物UBC825,UBC827和UBC842對(duì)于節(jié)莖曲柄蘚這8個(gè)居群遺傳多樣性的研究效果最好.本實(shí)驗(yàn)篩選出的12條ISSR引物,堿基A含量都非常豐富,其中5條含有豐富的堿基C,另外8條富含堿基G.

    2.2 居群的遺傳多樣性分析

    多樣性指數(shù)通常用于表示居群內(nèi)和居群間的遺傳變異程度[11].基于ISSR標(biāo)記的節(jié)莖曲柄蘚的遺傳多樣性指數(shù)如表3所示.由表3可知:8個(gè)居群的PPB變化范圍為9.30%~36.05%,PPB最高的居群是云南西雙版納納板河居群,為36.05%;廣西上思十萬(wàn)大山居群和浙江省清涼峰保護(hù)區(qū)的順溪塢居群次之,都為35.47%;接著依次是浙江舟山嵊泗列島居群、臺(tái)灣宜蘭縣馬告生態(tài)公園居群、浙江臺(tái)州臨海天臺(tái)山居群和江蘇蘇州天平山居群,分別為34.30%,30.23%,29.07%和13.95%;多態(tài)性位點(diǎn)百分率最低的是浙江清涼峰保護(hù)區(qū)涼源居群,僅為9.30%;各居群內(nèi)平均PPB為27.98%,物種水平上居群總的PPB為91.86%.

    I的變化范圍為0.0576~0.2210(表3),表明不同地理來(lái)源的節(jié)莖曲柄蘚居群內(nèi)的遺傳多樣性存在一定差異,其中I最高的是云南西雙版納納板河居群,接著依次是十萬(wàn)大山居群、嵊泗列島居群、順溪塢居群、馬告生態(tài)公園居群、天臺(tái)山居群和天平山居群,最低的是涼源居群.各個(gè)居群內(nèi)平均I為0.1635,物種水平上居群總的I為0.4459.

    節(jié)莖曲柄蘚居群的H的變化范圍為0.0400~0.1526(表3),其中H最高的是云南西雙版納納板河居群,接著依次是十萬(wàn)大山居群、嵊泗列島居群、順溪塢居群、天臺(tái)山居群、馬告生態(tài)公園居群和天平山居群,最低的是涼源居群.各個(gè)居群內(nèi)H的平均值為0.1114,物種水平上居群總的H為0.2922.

    由表3還可知,各居群的PPB,I和H比較一致,云南西雙版納納板河、廣西上思十萬(wàn)大山、浙江舟山嵊泗列島和浙江清涼峰順溪塢居群的遺傳多樣性指數(shù)比浙江臺(tái)州臨海天臺(tái)山、臺(tái)灣宜蘭縣馬告生態(tài)公園、江蘇蘇州天平山和浙江清涼峰涼源居群的高.僅在個(gè)別居群之間,3個(gè)參數(shù)略有差異,如I和H反映嵊泗列島居群的遺傳多樣性高于順溪塢居群,而PPB反映的結(jié)果則相反;PPB和I反映馬告生態(tài)公園居群的遺傳多樣性高于天臺(tái)山居群,而H反映的結(jié)果則相反.

    2.3 居群的遺傳分化分析

    H是衡量居群遺傳分化的最常用指標(biāo),常用來(lái)計(jì)算總的遺傳變異中,居群間變異所占的比例[10].表4顯示:由I估算的總遺傳變異中,63.33%的遺傳變異存在于居群間,36.67%存在于居群內(nèi);由H估算的遺傳分化系數(shù)表明總遺傳變異中,61.88%的遺傳變異存在于居群間,38.12%存在于居群內(nèi).I和H估算的結(jié)果一致,均表明節(jié)莖曲柄蘚的遺傳變異大部分存在于居群間.同時(shí),居群的基因流較小,僅為0.2879,說(shuō)明居群間的基因流動(dòng)水平低,居群之間的遺傳隔離明顯,不同的居群通過(guò)遺傳變異以適應(yīng)不同的環(huán)境,但是基因交流較少.

    2.4 居群的遺傳距離和聚類(lèi)分析

    基于ISSR標(biāo)記的8個(gè)節(jié)莖曲柄蘚居群的遺傳距離和遺傳相似系數(shù)見(jiàn)表5,節(jié)莖曲柄蘚8個(gè)居群間遺傳相似系數(shù)在0.6408~0.8766之間變化,平均值為0.7517,說(shuō)明節(jié)莖曲柄蘚8個(gè)居群間的遺傳相似性較小,居群間的遺傳分化明顯,這與2.3節(jié)的結(jié)果一致;其中浙江清涼峰保護(hù)區(qū)順溪塢居群和浙江臺(tái)州臨海天臺(tái)山居群遺傳相似系數(shù)最高,為0.8766,其次是順溪塢居群和涼源居群,遺傳相似系數(shù)為0.8587,江蘇蘇州天平山居群和浙江清涼峰保護(hù)區(qū)涼源居群的遺傳相似系數(shù)相對(duì)較小,為0.6408;8個(gè)居群的遺傳距離在0.1317~0.4451之間變化,平均值為0.2895.

    基于ISSR標(biāo)記,應(yīng)用POPGENE 3.2軟件,根據(jù)遺傳距離作出8個(gè)居群40個(gè)樣品的遺傳關(guān)系聚類(lèi)圖(圖1).以遺傳距離65為分組閾值,由圖1可知,這8個(gè)居群可以分為6組,分別是浙江清涼峰保護(hù)區(qū)的順溪塢居群、涼源居群和浙江臺(tái)州臨海天臺(tái)山居群聚集成的一組,剩下的5個(gè)居群(臺(tái)灣宜蘭馬告生態(tài)公園居群、浙江舟山嵊泗列島居群、江蘇蘇州天平山居群、云南納板河居群、廣西上思十萬(wàn)大山居群)分別自成一組.這8個(gè)居群主要分成兩個(gè)大的分支:分支①由順溪塢居群、涼源居群、臺(tái)州天臺(tái)山居群、臺(tái)灣宜蘭居群及舟山嵊泗列島居群組成;分支②由江蘇天平山居群、云南納板河居群和廣西十萬(wàn)大山居群聚集形成.分支①又由2個(gè)小的分支組成:同屬于浙江省的順溪塢居群、涼源居群和臺(tái)州天臺(tái)山居群聚為一個(gè)小分支;同屬于島嶼氣候的舟山嵊泗列島居群和臺(tái)灣宜蘭居群聚集為另一個(gè)小分支.分支②中江蘇天平山居群先和云南納板河居群聚集再與廣西十萬(wàn)大山居群聚集.研究結(jié)果顯示:節(jié)莖曲柄蘚的遺傳變異與地理位置及生境氣候相關(guān).

    基于節(jié)莖曲柄蘚40個(gè)樣品,8個(gè)居群的形態(tài)特征數(shù)據(jù),即ISSR電泳位點(diǎn)和8個(gè)居群分布的地理信息數(shù)據(jù),應(yīng)用R語(yǔ)言程序計(jì)算得到該種40個(gè)樣品的形態(tài)特征與ISSR數(shù)據(jù)相關(guān)性分析圖(圖2).由圖2可知,該種40個(gè)樣品的形態(tài)分化與其地理分布的相關(guān)性系數(shù)R為0.607,在P<0.05的水平上達(dá)到顯著相關(guān).

    3 討 論

    3.1 節(jié)莖曲柄蘚的遺傳多樣性

    在種水平上,應(yīng)用12個(gè)引物對(duì)節(jié)莖曲柄蘚8個(gè)居群40個(gè)樣品進(jìn)行了ISSR分析,共擴(kuò)增出172條清晰、重復(fù)性高的條帶,POPGENE 3.2分析結(jié)果表明:有158個(gè)為多態(tài)性位點(diǎn),PPB為91.86%;I為0.4459,居群內(nèi)平均值為0.1635;H為0.2922,居群內(nèi)平均值為0.1114.郭水良等[11]基于RAPD數(shù)據(jù)探討了中華蓑蘚4個(gè)居群的遺傳變異基礎(chǔ),中華蓑蘚的PPB為79.96%,物種水平的I為0.5126,4個(gè)居群的平均值為0.4055,總的H為0.3307,4個(gè)居群的平均值為0.2755,認(rèn)為中華蓑蘚的遺傳變異主要存在于居群之內(nèi).劉麗等[9]也應(yīng)用RAPD分子標(biāo)記對(duì)鼠尾蘚8個(gè)居群的遺傳多樣性進(jìn)行了研究,得到其PPB為84.60%,I為0.4877,H為0.3326,居群內(nèi)平均值為0.2228;李倩影[10]研究了小羽蘚屬遺傳多樣性,在104個(gè)檢測(cè)位點(diǎn)中發(fā)現(xiàn)94個(gè)多態(tài)性位點(diǎn),PPB為91.35%,總居群的H為0.2682,I為0.4108,該種居群表現(xiàn)出較高的遺傳多樣性,總的遺傳居群間遺傳分化系數(shù)為0.5591,居群內(nèi)遺傳分化系數(shù)為0.1182,居群間的變異程度遠(yuǎn)遠(yuǎn)大于居群內(nèi)的變異程度;汪琛穎等[4]對(duì)9個(gè)狹邊大葉蘚居群的遺傳多樣性進(jìn)行了ISSR分析,PPB為37.75%,H為0.1530,I為0.2228,遺傳分化指數(shù)為0.2997,遺傳變異主要存在于居群內(nèi),認(rèn)為該種的居群遺傳距離和地理距離沒(méi)有相關(guān)性.比較發(fā)現(xiàn),節(jié)莖曲柄蘚的遺傳多樣性與鼠尾蘚和小羽蘚屬相近,而高于狹邊大葉蘚,低于中華蓑蘚.可能是因?yàn)樵谝巴猓苌僖?jiàn)有節(jié)莖曲柄蘚的孢子體,說(shuō)明該種有性繁殖能力弱,配子體主要通過(guò)營(yíng)養(yǎng)繁殖擴(kuò)展種群,這種繁殖方式使其遺傳多樣性主要存在于居群之間.

    3.2 節(jié)莖曲柄蘚的遺傳結(jié)構(gòu)與分化

    從聚類(lèi)分析圖(圖1)可看出:節(jié)莖曲柄蘚的40個(gè)樣本8個(gè)居群以遺傳距離65為分組閾值可以分為6組,同屬于浙江省的順溪塢居群、涼源居群以及臺(tái)州臨海天臺(tái)山居群聚集為一組,說(shuō)明這3個(gè)居群的遺傳相似性高,遺傳距離近,這個(gè)結(jié)論與其相近的地理位置關(guān)系是吻合的;剩下的5個(gè)居群(臺(tái)灣宜蘭居群、浙江舟山嵊泗列島居群、江蘇天平山居群、云南納板河居群、廣西十萬(wàn)大山居群)各自分別形成5組,這5個(gè)居群的遺傳相似性低,遺傳距離遠(yuǎn),這與其所處的空間距離較遠(yuǎn)吻合;聚類(lèi)圖中分支②是由江蘇天平山居群、云南納板河居群、廣西十萬(wàn)大山居群聚集形成,云南納板河居群和廣西十萬(wàn)大山居群所處的地理位置毗鄰且處于同一緯度,而這兩個(gè)居群與江蘇天平山居群的氣候特點(diǎn)相似,說(shuō)明空間距離相近的居群,或者生長(zhǎng)環(huán)境相近的居群遺傳相似度高.分支①中:雖然舟山嵊泗列島也屬于浙江省,但其同臺(tái)灣宜蘭居群先聚集在一起,最終才與浙江省另外3個(gè)居群形成的分支聚集在一起,這是由于這兩個(gè)地區(qū)都屬于島嶼氣候.說(shuō)明節(jié)莖曲柄蘚的遺傳變異受地理生境和環(huán)境氣候影響大,這8個(gè)居群的基因流動(dòng)水平較低,居群間遺傳分化明顯,這與居群所處環(huán)境的異質(zhì)性有關(guān).同時(shí),形態(tài)特征與ISSR數(shù)據(jù)相關(guān)性分析圖(圖2)反映出該種40個(gè)樣品的形態(tài)分化與其地理分布的相關(guān)性指數(shù)為0.607,在P<0.05的水平上達(dá)到顯著相關(guān),也證明了節(jié)莖曲柄蘚8個(gè)居群間的遺傳分化與所處環(huán)境的異質(zhì)性有關(guān).

    WRIGHT[25]認(rèn)為,居群間基因流大于1時(shí),則能發(fā)揮其均質(zhì)化作用,反之若基因流小于1時(shí),則表明基因流是造成遺傳分化的主要因素,不同的基因型在不同的微生境上的適合度不同,導(dǎo)致具有相同基因型的個(gè)體聚集在相似的微生境上,從而產(chǎn)生居群間的遺傳分化.本研究中,節(jié)莖曲柄蘚居群的基因流為0.2879,遺傳分化系數(shù)為0.6346,說(shuō)明節(jié)莖曲柄蘚居群間雖然有一定的基因流動(dòng),但是流動(dòng)水平較低,且主要存在于生境相似的居群間,而在居群內(nèi)部,由于居群生境連續(xù)性差,像江蘇的天平山、清涼峰保護(hù)區(qū)的涼源,發(fā)生大量的遺傳漂變[4],從而導(dǎo)致該種植物遺傳變異的產(chǎn)生.

    3.3 島嶼生境的節(jié)莖曲柄蘚居群的遺傳多樣性分析

    由于島嶼的特殊環(huán)境,關(guān)于其苔蘚植物的遺傳多樣性研究較少.沿海島嶼,由于長(zhǎng)期處于隔離狀態(tài)下,并且相當(dāng)一部分苔蘚植物又很少產(chǎn)生孢子,無(wú)法基于孢子進(jìn)行遠(yuǎn)距離傳播,因而產(chǎn)生基因流減弱,苔蘚植物遺傳漂變等[26].WANG等[27]研究發(fā)現(xiàn):千島湖區(qū)破碎化的生境中,以無(wú)性繁殖為主的大灰蘚(Hypnum plumiforme Wilson),其島嶼居群的遺傳多樣性顯著高于湖邊居群,而經(jīng)常以有性繁殖為主的東亞小金發(fā)蘚(Pogonatum inflexum (Lindb.) Lac.),其島嶼群落和湖邊群落的遺傳分化程度不高.

    本研究中采自舟山嵊泗列島的節(jié)莖曲柄蘚居群,其位于長(zhǎng)江口與杭州灣的交匯處,在浙江省最東部、舟山群島最北部,處于華南褶皺系的次級(jí)構(gòu)造單元新昌—定?!B(niǎo)山斷隆的東北部邊緣,是浙東天臺(tái)山脈向東北延伸沉陷入海的外露部分[26].從表5可知,舟山嵊泗列島居群與臺(tái)州天臺(tái)山居群的遺傳相似性系數(shù)最高,其次是與同樣位于浙江的順溪塢和涼源居群.聚類(lèi)分析圖(圖1)也顯示:舟山嵊泗列島節(jié)莖曲柄蘚居群最后與浙江順溪塢和涼源居群聚類(lèi)在一起,說(shuō)明節(jié)莖曲柄蘚居群遺傳距離和地理距離有相關(guān)性.多態(tài)性位點(diǎn)百分率、根井正利基因多樣性指數(shù)和香濃-維納多樣性指數(shù)(表3)分析表明:舟山嵊泗列島節(jié)莖曲柄蘚居群的遺傳多樣性又高于浙江臺(tái)州天臺(tái)山居群.這與WANG等[27]對(duì)千島湖區(qū)大灰蘚居群遺傳多樣性的研究結(jié)果相符,說(shuō)明島嶼生境的長(zhǎng)期隔離使苔蘚植物發(fā)生遺傳漂變,這也印證了HUTSEMKERS等[28]的觀點(diǎn):海洋島嶼不是孢子植物遺傳多樣性的“洼地”,也不是進(jìn)化上的死胡同.

    本研究中采自舟山嵊泗列島的節(jié)莖曲柄蘚居群與采自臺(tái)灣的居群聚類(lèi)在一起,這是由島嶼環(huán)境的特殊性導(dǎo)致的.嵊泗列島屬亞熱帶海洋性季風(fēng)氣候,冬無(wú)畏寒,夏無(wú)酷暑,常年溫和濕潤(rùn),四季分明[29];臺(tái)灣是位于亞洲東部、太平洋西北側(cè)的島嶼,北回歸線將其劃為南北兩個(gè)氣候區(qū),北部屬亞熱帶季風(fēng)氣候,南部則為熱帶季風(fēng)氣候,因其四面環(huán)海,氣候又受海洋性季風(fēng)調(diào)節(jié).本研究所用樣品采自臺(tái)灣最北端的宜蘭縣馬告生態(tài)公園明池森林游樂(lè)區(qū),氣候特點(diǎn)與嵊泗列島相似,且經(jīng)度相近,故節(jié)莖曲柄蘚在這兩個(gè)地區(qū)的居群遺傳距離相近,在聚類(lèi)圖上聚集在一起.

    參考文獻(xiàn):

    [1] 王慧,徐杰.分子生物學(xué)在苔蘚植物分類(lèi)學(xué)中的應(yīng)用 [J].科學(xué)時(shí)代,2013(3):1-3.

    WANG H,XU J.Application of molecular biology in taxonomy of bryophytes [J].Science Times,2013(3):1-3.

    [2] 劉德好.中國(guó)主要稻區(qū)稗屬植物分類(lèi)與多樣性研究 [D].上海:上海師范大學(xué),2014.

    LIU D H.On classification and diversity of Echinochloa Beauv.in main rice regions of China [D].Shanghai:Shanghai Normal University,2014.

    [3] ZIETKIEWICZ E,RAFALSKI A,LABUDA D.Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification [J].Genomics,1994,20(2):176-183.

    [4] 汪琛穎,趙建成.狹邊大葉蘚(Rhodobryum ontariense)遺傳多樣性的ISSR分析 [J].貴州師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,28(4):33-40.

    WANG C Y,ZHAO J C.Genetic diversity in a medicinal plant species:Rhodobryum ontariense (Bryaceae),detected by ISSR[J].Journal of Guizhou Normal University (Natural Sciences),2010,28(4):33-40.

    [5] 吳雪琴.南嶺地區(qū)觀光木居種群遺傳多樣性的ISSR分析 [D].長(zhǎng)沙:中南林業(yè)科技大學(xué),2013.

    WU X Q.ISSR analysis on genetic diversity of Tsoongiodendron odorum chun populations from the Nanling Mountains [D].Changsha:Central South University of Forestry and Technology,2013.

    [6] 王仕玉,郭鳳根,張應(yīng)華,等.18份滇產(chǎn)巖白菜資源的ISSR分析 [J].中國(guó)農(nóng)學(xué)通報(bào),2012,28(7):114-118.

    WANG S Y,GUO F G,ZHANG Y H,et al.ISSR analysis of 18 Bergenia purpurascens resources in Yunnan [J].Chinese Agricultural Science Bulletin,2012,28(7):114-118.

    [7] 吳祝華,施季森,席夢(mèng)利,等.百合屬野生種及品種親緣關(guān)系的 ISSR 分析 [J].江蘇農(nóng)業(yè)科學(xué),2012,40(1):19-22.

    WU Z H,SHI J S,XI M L,et al.ISSR analysis of the relationship between wild species and varieties of Lilium [J].Jiangsu Agricultural Science,2012,40(1):19-22.

    [8] SKOTNICKI M,BARGAGLI R,NINHAM J.Genetic diversity in the moss Pohlia nutans on geothermal ground of Mount Rittmann,Victoria Land,Antarctica [J].Polar Biology,2002,25(10):771-777.

    [9] 劉麗,朱永青,王幼芳.鼠尾蘚不同居群間形態(tài)及RAPD分析 [J].云南植物研究,2006,28(6):570-574.

    LIU L,ZHU Y Q,WANG Y F.Morphological and RAPD analyses of different Myuroclada maximowiczii (Brachytheciaceae) populations [J].Acta Botanica Yunnanica,2006,28(6):570-574.

    [10] 李倩影.小羽蘚屬(Haplocladium)植物遺傳多樣性及其與環(huán)境因子關(guān)系分析 [D].上海:上海師范大學(xué),2010.

    LI Q Y.Genetic diodiversity of moss genus Haplocladium and its relationship with the environmental factors [D].Shanghai:Shanghai Normal University,2010.

    [11] 郭水良,鄭園園,婁玉霞,等.基于RAPD 數(shù)據(jù)探討中華蓑蘚(Macromitrium cavaleriei Card.& Ther.) 形態(tài)變異的遺傳基礎(chǔ) [J].植物科學(xué)學(xué)報(bào),2011,29(3):312-318.

    GUO S L,ZHENG Y Y,LOU Y X,et al.Genetic background of morphological variation in Macromitrium cavaleriei Card.& Ther.) based on RAPD [J].Plant Science Journal,2011,29(3):312-318.

    [12] 汪琛穎,趙建成.真蘚科植物ISSR-PCR反應(yīng)體系的優(yōu)化及ISSR指紋圖譜的初步構(gòu)建 [J].安徽農(nóng)業(yè)科學(xué),2011,39(27):16490-16493.

    WANG C Y,ZHAO J C.Optimization of ISSR-PCR reaction system and preliminary construction of ISSR fingerprinting of some species in Bryaceae [J].Journal of Anhui Agricultural Sciences,2011,39(27):16490-16493.

    [13] 王瑩瑩.浙江千島湖生境片斷化對(duì)苔蘚植物物種及遺傳多樣性的影響 [D].上海:華東師范大學(xué),2011.

    WANG Y Y.Effects of habitat fragmentation on bryophytes richness and genetic diversity in Thousand-Island Lake Region,Zhejiang [D].Shanghai:East China Normal University,2011.

    [14] 買(mǎi)買(mǎi)提明·蘇來(lái)曼,徐紅紅,瑪爾孜亞·阿不力米提,等.水蘚科植物遺傳多樣性研究 [J].新疆大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,30(4):390-395.

    MAMTIMIN S,XU H H,MARZIYA A,et al.Genetic diversity of fontinalaceae by ISSR analysis [J].Journal of Xinjiang University (Natural Science),2013,30(4):390-395.

    [15] 魏青永,郭水良,曹同,等.基于 ISSR 數(shù)據(jù)探討卷葉鳳尾蘚 (Fissidens dubius P.Beauv.) 遺傳多樣性 [J].植物科學(xué)學(xué)報(bào),2016,34(2):238-245.

    WEI Q Y,GUO S L,CAO T,et al.Genetic diversity of Fissidens dubius P.Beauv.based on ISSR data [J].Plant Science Journal,2016,34(2):238-245.

    [16] BRIDEL S E.Muscologiae Recentiorum Supplementum [M].Gotha:C.G.Ettinger,1819,4:220.

    [17] FRAHM J P.A revision of the East-Asian species of Campylopus [J].Journal of the Hattori Botanical Laboratory,1992,71:133-164.

    [18] CROSBY M R,MAGILL R E,ALLEN B,et al.A Checklist of the Mosses [M].St.Louis:Missouri Botanical Garden,2000:55-246.

    [19] FRAHM J P.The Campylopus species of the world-a picture book:photographs of 59 species,subspecies and varieties [J].Archive for Bryology,2012(Special 9):1-63.

    [20] 賈渝,何思.中國(guó)生物物種名錄:第1卷植物,苔蘚植物 [M].北京:科學(xué)出版社,2013.

    JIA Y,HE S.Species Catalogue of China:Volume 1 Plants,Bryophytes [M].Beijing:Science Press,2013.

    [21] 高謙.中國(guó)苔蘚植物志:第1卷 [M].北京:科學(xué)出版社,1994.

    GAO C.Flora Bryophytarum Sinicorum:Volume 1 [M].Beijing:Science Press,1994.

    [22] BROTHERUS V F.Musci novi sinenses collecti a Dr.Henr.Handel-Mazzetti.II [J].Sitzungsberichte,Akademie der Wissenschaften in Wien.Mathematisch-naturwissenschaftliche Klasse.Abteilung 1,1924,133:559-584.

    [23] DUBY J .Choix de cryptogames exotiques nouvelles ou mal counues:I.mousses [J].Mémoires de la Société de Physique et d'Histoire Naturelle de Genève,1868,19:291-304.

    [24] DIXON H N.Mosses of Hong Kong with other Chinese mosses [J].Hong Kong Naturalist,1933,2(Suppl.):1-31.

    [25] WRIGHT S.Isolation by distance under diverse systems of mating [J].Genetics,1946,31(1):39-59.

    [26] 于晶,孫鵬雁,郭水良.島嶼苔蘚植物多樣性和生物地理學(xué)研究的意義、概況和方向 [J].上海師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,45(1):70-75.

    YU J,SUN P Y,GUO S L.Studies on diversity and biogeography of island bryophytes:significances,progress and dircetion [J].Journal of Shanghai Normal University (Natural Sciences),2016,45(1):70-75.

    [27] WANG Y Y,ZHU Y Q,WANG Y F.Differences in spatial genetic structure and diversity in two mosses with different dispersal strategies in a fragmented landscape [J].Journal of Bryology,2012,34(1):9-16.

    [28] HUTSEMKERS V,SZVNYI P,SHWA A J,et al.Oceanic islands are not sinks of biodiversity in spore-producing plants [J].Proceedings of the National Academy of Sciences,2011,108(47):18989-18994.

    [29] 任淑華,蔡克勤.嵊泗列島旅游氣候資源開(kāi)發(fā)利用 [J].海洋開(kāi)發(fā)與管理,2008,25(5):112-114.

    REN S H,CAI K Q.Development and utilization of tourism climate resources in the Shengsi Islands [J].Ocean Development and Management,2008,25(5):112-114.

    (責(zé)任編輯:顧浩然)

    猜你喜歡
    嵊泗列島遺傳變異
    先導(dǎo)編輯技術(shù)可編輯近90%的人類(lèi)遺傳變異
    基于改進(jìn)遺傳變異算子的海島算法
    NOWHERE LAND
    日本列島封土石室墓的展開(kāi)與連云港石室墓
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    長(zhǎng)山列島國(guó)家地質(zhì)公園
    国产午夜精品论理片| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 在线观看免费高清a一片| 午夜福利在线观看免费完整高清在| 欧美日韩亚洲高清精品| 69人妻影院| 久久精品久久久久久噜噜老黄| 欧美最新免费一区二区三区| 日日摸夜夜添夜夜爱| 我的老师免费观看完整版| 极品少妇高潮喷水抽搐| 好男人视频免费观看在线| 精品久久久久久电影网| 大话2 男鬼变身卡| 人人妻人人看人人澡| 三级毛片av免费| 国产黄片视频在线免费观看| 91aial.com中文字幕在线观看| 在线a可以看的网站| 成人综合一区亚洲| 午夜爱爱视频在线播放| 有码 亚洲区| av在线蜜桃| 一个人看视频在线观看www免费| 国产不卡一卡二| 精品亚洲乱码少妇综合久久| 亚洲av电影在线观看一区二区三区 | 特级一级黄色大片| 小蜜桃在线观看免费完整版高清| 男人爽女人下面视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 成年女人在线观看亚洲视频 | 国产黄色小视频在线观看| 国产亚洲5aaaaa淫片| 成人综合一区亚洲| 久久精品夜色国产| 欧美一区二区亚洲| 国产精品av视频在线免费观看| 一个人看的www免费观看视频| 成人午夜精彩视频在线观看| 插逼视频在线观看| 黄色欧美视频在线观看| 三级国产精品片| 嫩草影院新地址| 777米奇影视久久| 免费看光身美女| 久久精品国产鲁丝片午夜精品| 精品久久久精品久久久| 久久人人爽人人爽人人片va| 免费高清在线观看视频在线观看| 一级毛片黄色毛片免费观看视频| 国产综合懂色| 日韩 亚洲 欧美在线| 国产一区亚洲一区在线观看| 精品人妻一区二区三区麻豆| 午夜亚洲福利在线播放| 亚洲欧洲国产日韩| 男女啪啪激烈高潮av片| 欧美97在线视频| 人妻少妇偷人精品九色| 免费看美女性在线毛片视频| 亚洲欧洲国产日韩| 蜜臀久久99精品久久宅男| 午夜福利网站1000一区二区三区| 欧美激情国产日韩精品一区| 日本三级黄在线观看| 国产麻豆成人av免费视频| 日日摸夜夜添夜夜爱| 久久久久久久久久久免费av| 又爽又黄无遮挡网站| 亚洲国产成人一精品久久久| 免费看日本二区| 国产黄色视频一区二区在线观看| 男插女下体视频免费在线播放| 精品久久国产蜜桃| 男女下面进入的视频免费午夜| 亚洲电影在线观看av| 日韩在线高清观看一区二区三区| av在线老鸭窝| 国产69精品久久久久777片| 99久国产av精品国产电影| 搡老乐熟女国产| 欧美一区二区亚洲| 亚洲国产欧美在线一区| 美女主播在线视频| 纵有疾风起免费观看全集完整版 | 国产一级毛片在线| 美女内射精品一级片tv| 亚洲国产欧美在线一区| 又大又黄又爽视频免费| 国内精品美女久久久久久| 干丝袜人妻中文字幕| 国产在视频线精品| 亚洲精品影视一区二区三区av| 久久久精品94久久精品| 国产 亚洲一区二区三区 | 在线天堂最新版资源| 成人亚洲精品一区在线观看 | www.av在线官网国产| 亚洲人成网站在线观看播放| 蜜桃久久精品国产亚洲av| 非洲黑人性xxxx精品又粗又长| 国产黄色视频一区二区在线观看| 成人鲁丝片一二三区免费| av专区在线播放| 国产 一区 欧美 日韩| 嘟嘟电影网在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品影视一区二区三区av| 视频中文字幕在线观看| 男女边吃奶边做爰视频| 热99在线观看视频| 国产真实伦视频高清在线观看| 男女边摸边吃奶| 成人亚洲欧美一区二区av| 色综合色国产| 久久久久免费精品人妻一区二区| 国产精品人妻久久久久久| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品一区二区| 国产在线男女| 国产激情偷乱视频一区二区| 别揉我奶头 嗯啊视频| 晚上一个人看的免费电影| 麻豆av噜噜一区二区三区| 成人av在线播放网站| av国产久精品久网站免费入址| 日韩强制内射视频| 色哟哟·www| 亚洲熟妇中文字幕五十中出| 男人舔奶头视频| 久久久久久九九精品二区国产| 免费av不卡在线播放| 少妇高潮的动态图| 国产成人freesex在线| 18+在线观看网站| 秋霞伦理黄片| 日日啪夜夜爽| 中文字幕免费在线视频6| 天堂俺去俺来也www色官网 | 欧美xxxx性猛交bbbb| 久久精品久久久久久久性| 在现免费观看毛片| 五月天丁香电影| 午夜免费激情av| 亚洲aⅴ乱码一区二区在线播放| 视频中文字幕在线观看| 亚洲精品乱久久久久久| 亚洲国产精品专区欧美| 国产黄色免费在线视频| 国产免费福利视频在线观看| 亚洲国产精品成人久久小说| 国产伦精品一区二区三区视频9| 亚洲精品视频女| 国产乱人视频| 欧美潮喷喷水| 欧美丝袜亚洲另类| 简卡轻食公司| 国内精品宾馆在线| 一级毛片我不卡| 国产精品嫩草影院av在线观看| 午夜福利在线观看吧| 国产精品久久久久久av不卡| 亚洲国产精品成人综合色| 男插女下体视频免费在线播放| 一区二区三区乱码不卡18| 日韩精品青青久久久久久| 国产精品国产三级专区第一集| 国产真实伦视频高清在线观看| 久久精品国产亚洲av涩爱| 99re6热这里在线精品视频| 伊人久久国产一区二区| 99久久人妻综合| 直男gayav资源| 欧美日韩亚洲高清精品| 晚上一个人看的免费电影| 99久国产av精品| 免费看日本二区| 国产乱来视频区| 国产精品一区www在线观看| 最近的中文字幕免费完整| 国产乱来视频区| 99热全是精品| 亚洲图色成人| 免费黄频网站在线观看国产| 一本一本综合久久| 久久99蜜桃精品久久| 国内精品美女久久久久久| 人体艺术视频欧美日本| 欧美xxⅹ黑人| 久久久久久九九精品二区国产| 哪个播放器可以免费观看大片| 国产亚洲午夜精品一区二区久久 | 成人亚洲精品一区在线观看 | 免费看光身美女| 日韩在线高清观看一区二区三区| 中国美白少妇内射xxxbb| 国产成人一区二区在线| av在线亚洲专区| 久久精品国产亚洲av涩爱| 免费av观看视频| 国产精品不卡视频一区二区| 日韩成人伦理影院| 色综合色国产| 日本免费在线观看一区| 久久久午夜欧美精品| 亚洲av中文av极速乱| 欧美bdsm另类| 九色成人免费人妻av| 搡老乐熟女国产| 国产免费视频播放在线视频 | 天天一区二区日本电影三级| 美女被艹到高潮喷水动态| 一个人免费在线观看电影| 有码 亚洲区| 777米奇影视久久| 国产精品三级大全| 九色成人免费人妻av| 久久这里只有精品中国| 中文天堂在线官网| 亚洲成人精品中文字幕电影| 色综合亚洲欧美另类图片| 久久久久久久午夜电影| 亚洲怡红院男人天堂| 你懂的网址亚洲精品在线观看| 亚洲精品乱码久久久久久按摩| av在线观看视频网站免费| 亚洲真实伦在线观看| 精品国产三级普通话版| 免费不卡的大黄色大毛片视频在线观看 | 久久热精品热| 你懂的网址亚洲精品在线观看| av在线蜜桃| 国产爱豆传媒在线观看| av播播在线观看一区| 久久国产乱子免费精品| 又大又黄又爽视频免费| 我的女老师完整版在线观看| 久久人人爽人人爽人人片va| 一区二区三区免费毛片| 国产成人一区二区在线| 一级av片app| 成年人午夜在线观看视频 | 三级国产精品片| 亚洲图色成人| 午夜精品一区二区三区免费看| 久久久欧美国产精品| 成人高潮视频无遮挡免费网站| 大陆偷拍与自拍| freevideosex欧美| www.色视频.com| 草草在线视频免费看| 男女边摸边吃奶| 欧美精品一区二区大全| 黄片无遮挡物在线观看| eeuss影院久久| 精品人妻视频免费看| 国产亚洲5aaaaa淫片| 色综合亚洲欧美另类图片| 成人特级av手机在线观看| 成人午夜高清在线视频| 精品久久久精品久久久| 一级黄片播放器| 色尼玛亚洲综合影院| 成人午夜精彩视频在线观看| 欧美不卡视频在线免费观看| 日本黄大片高清| 黄片无遮挡物在线观看| 色综合亚洲欧美另类图片| 国产av不卡久久| 国产av国产精品国产| 精品一区二区三区视频在线| 中文字幕制服av| 波多野结衣巨乳人妻| 看免费成人av毛片| 日韩视频在线欧美| 全区人妻精品视频| 亚洲欧洲国产日韩| xxx大片免费视频| 午夜亚洲福利在线播放| 午夜精品在线福利| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| 亚洲精品日韩在线中文字幕| 久久久精品欧美日韩精品| 亚洲av福利一区| 97在线视频观看| av卡一久久| 简卡轻食公司| 欧美成人午夜免费资源| 中国美白少妇内射xxxbb| 丰满乱子伦码专区| 久久久久久久久久久免费av| 1000部很黄的大片| 国产黄频视频在线观看| 免费av毛片视频| 成人漫画全彩无遮挡| 国产精品嫩草影院av在线观看| 欧美 日韩 精品 国产| 99久久精品一区二区三区| 天堂中文最新版在线下载 | 有码 亚洲区| 日本av手机在线免费观看| 国产精品美女特级片免费视频播放器| 亚洲人与动物交配视频| 国产大屁股一区二区在线视频| 亚洲av中文字字幕乱码综合| 人体艺术视频欧美日本| 国产一区二区三区av在线| 免费不卡的大黄色大毛片视频在线观看 | 日本爱情动作片www.在线观看| 五月天丁香电影| 韩国av在线不卡| 欧美极品一区二区三区四区| 国产成人午夜福利电影在线观看| 九九在线视频观看精品| 久久99精品国语久久久| 自拍偷自拍亚洲精品老妇| 国产一区二区在线观看日韩| 亚洲伊人久久精品综合| 777米奇影视久久| 欧美日韩国产mv在线观看视频 | 国产有黄有色有爽视频| 久久久久久伊人网av| 欧美xxxx黑人xx丫x性爽| 免费观看的影片在线观看| 亚洲美女搞黄在线观看| 内射极品少妇av片p| 男女啪啪激烈高潮av片| 成人高潮视频无遮挡免费网站| 永久网站在线| 日韩欧美国产在线观看| 亚洲最大成人av| 日韩亚洲欧美综合| 直男gayav资源| 不卡视频在线观看欧美| 久久久色成人| 亚洲av电影不卡..在线观看| 尾随美女入室| 自拍偷自拍亚洲精品老妇| 亚洲精品色激情综合| 精品久久久久久久末码| 色视频www国产| 亚洲不卡免费看| 欧美不卡视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 韩国av在线不卡| 久久久久久久亚洲中文字幕| 国产白丝娇喘喷水9色精品| 亚洲自拍偷在线| 激情五月婷婷亚洲| 99久久精品国产国产毛片| videos熟女内射| av福利片在线观看| 街头女战士在线观看网站| 嫩草影院精品99| 一级黄片播放器| 亚洲av中文字字幕乱码综合| 久久久久久久久久人人人人人人| 国产一区二区在线观看日韩| 偷拍熟女少妇极品色| 久久久久性生活片| 日韩在线高清观看一区二区三区| 在线观看av片永久免费下载| 一级片'在线观看视频| 插阴视频在线观看视频| 精品久久久久久久久久久久久| 亚洲色图av天堂| 成人国产麻豆网| 精品久久久久久久末码| 99热这里只有精品一区| 国产精品一区二区在线观看99 | 亚洲欧美成人综合另类久久久| 纵有疾风起免费观看全集完整版 | 乱人视频在线观看| 韩国av在线不卡| 久久亚洲国产成人精品v| 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 中文乱码字字幕精品一区二区三区 | 天堂网av新在线| 神马国产精品三级电影在线观看| 欧美日韩国产mv在线观看视频 | 亚洲国产高清在线一区二区三| 免费观看a级毛片全部| 一个人免费在线观看电影| 国产v大片淫在线免费观看| 尾随美女入室| av播播在线观看一区| 欧美zozozo另类| 别揉我奶头 嗯啊视频| 欧美区成人在线视频| 国产国拍精品亚洲av在线观看| 国产精品一二三区在线看| 久久久久久久久中文| 2021天堂中文幕一二区在线观| 99九九线精品视频在线观看视频| 婷婷色av中文字幕| 女人久久www免费人成看片| 亚洲伊人久久精品综合| 美女cb高潮喷水在线观看| 日韩欧美国产在线观看| 你懂的网址亚洲精品在线观看| 九色成人免费人妻av| 精品久久国产蜜桃| av国产久精品久网站免费入址| 一级爰片在线观看| 97人妻精品一区二区三区麻豆| 一区二区三区免费毛片| 成人午夜精彩视频在线观看| 婷婷色综合大香蕉| 日日干狠狠操夜夜爽| 乱码一卡2卡4卡精品| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久久末码| 亚洲在线自拍视频| 久久久午夜欧美精品| 亚洲国产精品成人综合色| 亚洲国产精品国产精品| 日本黄大片高清| 国产黄色小视频在线观看| 99久久人妻综合| 国产一级毛片七仙女欲春2| 一级毛片电影观看| 美女高潮的动态| 国产黄片美女视频| 亚洲av福利一区| 麻豆乱淫一区二区| 中文字幕免费在线视频6| 黄色一级大片看看| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 三级国产精品片| 免费大片黄手机在线观看| 中文字幕av成人在线电影| 亚洲精品成人久久久久久| 成人午夜高清在线视频| 又爽又黄a免费视频| ponron亚洲| 午夜亚洲福利在线播放| 国产成人精品一,二区| 国产亚洲5aaaaa淫片| 午夜老司机福利剧场| 亚洲精华国产精华液的使用体验| 最新中文字幕久久久久| 2021少妇久久久久久久久久久| 成人午夜精彩视频在线观看| 日本wwww免费看| 午夜激情欧美在线| 免费无遮挡裸体视频| 国产一区二区亚洲精品在线观看| 18禁动态无遮挡网站| 欧美精品一区二区大全| 亚洲性久久影院| 国产一区二区在线观看日韩| 国产老妇伦熟女老妇高清| 特级一级黄色大片| 中文精品一卡2卡3卡4更新| 国产69精品久久久久777片| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 精品久久久久久久末码| 一个人免费在线观看电影| 久久精品国产亚洲av天美| 欧美三级亚洲精品| 亚洲av一区综合| 国产精品美女特级片免费视频播放器| 国产午夜精品一二区理论片| 2021天堂中文幕一二区在线观| 在线观看免费高清a一片| 精品国内亚洲2022精品成人| 国产精品无大码| 精品久久久久久久末码| 天堂av国产一区二区熟女人妻| 2021天堂中文幕一二区在线观| 午夜福利视频精品| 国产日韩欧美在线精品| 日日啪夜夜撸| 少妇猛男粗大的猛烈进出视频 | 成人午夜高清在线视频| 免费黄频网站在线观看国产| 亚洲精品视频女| 女人十人毛片免费观看3o分钟| 国产成人a区在线观看| 久久99精品国语久久久| 91精品伊人久久大香线蕉| 亚洲丝袜综合中文字幕| 3wmmmm亚洲av在线观看| 听说在线观看完整版免费高清| 大香蕉97超碰在线| 人妻一区二区av| 国产成人免费观看mmmm| 成人性生交大片免费视频hd| 99视频精品全部免费 在线| 国产av国产精品国产| 伦精品一区二区三区| a级毛片免费高清观看在线播放| 国产免费又黄又爽又色| 亚洲欧美成人综合另类久久久| 天堂俺去俺来也www色官网 | 亚洲av日韩在线播放| 日日摸夜夜添夜夜爱| 乱码一卡2卡4卡精品| 亚洲天堂国产精品一区在线| av网站免费在线观看视频 | 能在线免费看毛片的网站| 亚洲精品一二三| 97超视频在线观看视频| 91精品国产九色| 中文字幕av成人在线电影| 六月丁香七月| 久久鲁丝午夜福利片| 亚洲精品久久久久久婷婷小说| 亚洲aⅴ乱码一区二区在线播放| a级一级毛片免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲一区二区精品| 大话2 男鬼变身卡| 亚洲欧美日韩卡通动漫| 国产探花在线观看一区二区| 婷婷色综合大香蕉| 日本欧美国产在线视频| 黄色一级大片看看| 国产一区二区亚洲精品在线观看| 99热这里只有精品一区| 国内少妇人妻偷人精品xxx网站| 日日干狠狠操夜夜爽| 国产欧美日韩精品一区二区| 亚洲av成人精品一二三区| 久久久久性生活片| 九九在线视频观看精品| 啦啦啦中文免费视频观看日本| 免费av观看视频| 高清av免费在线| av线在线观看网站| 91aial.com中文字幕在线观看| 国产精品美女特级片免费视频播放器| 夫妻性生交免费视频一级片| 欧美性猛交╳xxx乱大交人| 丰满少妇做爰视频| av线在线观看网站| 婷婷色综合大香蕉| 99热全是精品| 成人毛片a级毛片在线播放| 亚洲国产最新在线播放| 亚洲色图av天堂| 偷拍熟女少妇极品色| 99久久精品一区二区三区| 亚洲一级一片aⅴ在线观看| 久久久久久久国产电影| 国产精品无大码| 亚洲欧美一区二区三区国产| 欧美三级亚洲精品| 亚洲高清免费不卡视频| 欧美人与善性xxx| 狂野欧美激情性xxxx在线观看| 久久鲁丝午夜福利片| 日韩大片免费观看网站| 日本三级黄在线观看| 91aial.com中文字幕在线观看| 亚洲成人久久爱视频| 网址你懂的国产日韩在线| 一个人免费在线观看电影| 免费看不卡的av| 国产亚洲午夜精品一区二区久久 | 久久精品夜色国产| 男的添女的下面高潮视频| 久久久色成人| 欧美xxⅹ黑人| 亚洲四区av| 精品国产露脸久久av麻豆 | 精品人妻一区二区三区麻豆| 在线观看一区二区三区| 丝瓜视频免费看黄片| 精品国产三级普通话版| 内地一区二区视频在线| 少妇丰满av| 精品99又大又爽又粗少妇毛片| 日韩制服骚丝袜av| 在线观看免费高清a一片| 国产高清有码在线观看视频| 最近中文字幕2019免费版| 在线 av 中文字幕| 婷婷色麻豆天堂久久| 狠狠精品人妻久久久久久综合| 一级黄片播放器| 久久久亚洲精品成人影院| 国产人妻一区二区三区在| 精品人妻偷拍中文字幕| 水蜜桃什么品种好| 麻豆久久精品国产亚洲av| 国产精品美女特级片免费视频播放器| av网站免费在线观看视频 | 最近的中文字幕免费完整| 精品人妻偷拍中文字幕| 日韩 亚洲 欧美在线| 日韩精品青青久久久久久| 亚洲精品第二区| av在线观看视频网站免费| 欧美xxxx性猛交bbbb| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| ponron亚洲| 韩国av在线不卡| 日韩中字成人| 日韩欧美国产在线观看| 国产av不卡久久| 国产又色又爽无遮挡免| 国产高清国产精品国产三级 | 国产精品一及| 蜜桃亚洲精品一区二区三区| 97人妻精品一区二区三区麻豆| 国产v大片淫在线免费观看|