• <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ì)公園
    久久精品国产亚洲av涩爱| 女人十人毛片免费观看3o分钟| 亚洲18禁久久av| 又粗又硬又长又爽又黄的视频| 亚洲av免费在线观看| 日韩三级伦理在线观看| 日本猛色少妇xxxxx猛交久久| 九九热线精品视视频播放| 观看美女的网站| 国产亚洲av嫩草精品影院| 老司机福利观看| 国产一区有黄有色的免费视频 | 最新中文字幕久久久久| 寂寞人妻少妇视频99o| 建设人人有责人人尽责人人享有的| 热re99久久精品国产66热6| 大香蕉97超碰在线| 少妇 在线观看| 飞空精品影院首页| 亚洲第一av免费看| 少妇被粗大的猛进出69影院 | 午夜免费观看性视频| 国产色婷婷99| 蜜桃国产av成人99| 一边摸一边做爽爽视频免费| 狠狠精品人妻久久久久久综合| 午夜影院在线不卡| 日韩熟女老妇一区二区性免费视频| 亚洲,欧美,日韩| 欧美激情极品国产一区二区三区 | 一区二区三区乱码不卡18| 欧美日韩av久久| 欧美日韩亚洲高清精品| 天天影视国产精品| 91精品三级在线观看| 九九在线视频观看精品| 国产一区二区三区av在线| 欧美bdsm另类| √禁漫天堂资源中文www| 亚洲婷婷狠狠爱综合网| 国产69精品久久久久777片| 国产亚洲av片在线观看秒播厂| 国产精品久久久久久久久免| 欧美成人精品欧美一级黄| 99久久人妻综合| 国产精品不卡视频一区二区| 国产成人精品福利久久| 亚洲精品日韩在线中文字幕| 国产精品人妻久久久影院| xxxhd国产人妻xxx| 精品一区二区三区视频在线| 丰满饥渴人妻一区二区三| 涩涩av久久男人的天堂| 欧美变态另类bdsm刘玥| 伦精品一区二区三区| 国产精品久久久久久久电影| 男女免费视频国产| 天堂中文最新版在线下载| 欧美精品一区二区大全| 国产一区二区在线观看日韩| 国产国拍精品亚洲av在线观看| 久久久久人妻精品一区果冻| 99久久综合免费| 嫩草影院入口| 热re99久久精品国产66热6| 成年女人在线观看亚洲视频| 国产白丝娇喘喷水9色精品| 中文欧美无线码| 日韩av免费高清视频| 咕卡用的链子| 在线亚洲精品国产二区图片欧美| 国产成人午夜福利电影在线观看| 久久人妻熟女aⅴ| 婷婷色综合大香蕉| 久久久国产一区二区| 免费观看性生交大片5| 少妇 在线观看| 精品人妻一区二区三区麻豆| 尾随美女入室| 天天躁夜夜躁狠狠久久av| 日韩大片免费观看网站| 欧美日韩国产mv在线观看视频| 纯流量卡能插随身wifi吗| 国产精品国产av在线观看| 在线天堂最新版资源| 日韩精品有码人妻一区| 国产精品国产三级国产av玫瑰| 免费av中文字幕在线| 国产在视频线精品| 免费日韩欧美在线观看| 狠狠精品人妻久久久久久综合| 久久久精品免费免费高清| 免费观看在线日韩| 欧美人与性动交α欧美精品济南到 | 国产1区2区3区精品| 欧美国产精品va在线观看不卡| 韩国精品一区二区三区 | 成人影院久久| 在线观看免费高清a一片| 色婷婷久久久亚洲欧美| 丝瓜视频免费看黄片| 夜夜爽夜夜爽视频| av有码第一页| 国产精品偷伦视频观看了| 大片电影免费在线观看免费| 99九九在线精品视频| 一级黄片播放器| 久久久久久人人人人人| 99久久人妻综合| 在线观看免费视频网站a站| 人人妻人人爽人人添夜夜欢视频| 看免费成人av毛片| 免费人妻精品一区二区三区视频| av播播在线观看一区| 欧美人与性动交α欧美精品济南到 | 五月天丁香电影| 狂野欧美激情性bbbbbb| 晚上一个人看的免费电影| 国产免费又黄又爽又色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av黄色大香蕉| 岛国毛片在线播放| 国产熟女午夜一区二区三区| 青青草视频在线视频观看| 日日爽夜夜爽网站| 国产精品一区www在线观看| 老女人水多毛片| 免费高清在线观看日韩| 久久精品久久精品一区二区三区| 日本欧美国产在线视频| 丝瓜视频免费看黄片| 十八禁网站网址无遮挡| 22中文网久久字幕| 女人精品久久久久毛片| 黄色怎么调成土黄色| 亚洲国产色片| 在线亚洲精品国产二区图片欧美| kizo精华| 亚洲成人av在线免费| 男人操女人黄网站| 视频中文字幕在线观看| 欧美丝袜亚洲另类| 国产午夜精品一二区理论片| 看十八女毛片水多多多| 国产精品无大码| 91久久精品国产一区二区三区| av不卡在线播放| 久久国产精品男人的天堂亚洲 | 亚洲综合色惰| 美女xxoo啪啪120秒动态图| 一级黄片播放器| 国产毛片在线视频| 午夜av观看不卡| 久久精品久久精品一区二区三区| 亚洲精品中文字幕在线视频| 久久久精品免费免费高清| 国产成人精品在线电影| 成人亚洲精品一区在线观看| 久久99一区二区三区| 日韩中字成人| 国产日韩一区二区三区精品不卡| 国产精品一区www在线观看| 丝袜喷水一区| 国产成人精品在线电影| 人人妻人人澡人人看| 国精品久久久久久国模美| 激情视频va一区二区三区| 国产av一区二区精品久久| 成人毛片a级毛片在线播放| 国产男人的电影天堂91| 国产在视频线精品| 涩涩av久久男人的天堂| 在线天堂最新版资源| 色婷婷av一区二区三区视频| 午夜老司机福利剧场| 十八禁高潮呻吟视频| 国产免费现黄频在线看| 免费在线观看黄色视频的| 久久 成人 亚洲| 久久久a久久爽久久v久久| 日韩av不卡免费在线播放| 久久97久久精品| 成人综合一区亚洲| 亚洲国产av新网站| 51国产日韩欧美| 九草在线视频观看| 久久ye,这里只有精品| 久久久久久人妻| 久久99一区二区三区| 黄色怎么调成土黄色| 亚洲欧美色中文字幕在线| 欧美日韩成人在线一区二区| 国产精品嫩草影院av在线观看| 欧美激情 高清一区二区三区| 精品一区二区三区四区五区乱码 | 99国产综合亚洲精品| www.熟女人妻精品国产 | 18禁裸乳无遮挡动漫免费视频| 免费黄色在线免费观看| 永久网站在线| 男男h啪啪无遮挡| 秋霞伦理黄片| av免费在线看不卡| 久久久久久久大尺度免费视频| 男女免费视频国产| videossex国产| 久久毛片免费看一区二区三区| 视频区图区小说| 国产黄色视频一区二区在线观看| 国产熟女欧美一区二区| 国精品久久久久久国模美| 男女无遮挡免费网站观看| 国产一区二区在线观看日韩| a级毛色黄片| 美女脱内裤让男人舔精品视频| 国产亚洲最大av| 久久精品国产鲁丝片午夜精品| 丝袜人妻中文字幕| 欧美国产精品va在线观看不卡| 久久久久精品人妻al黑| 有码 亚洲区| 高清毛片免费看| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 国产av一区二区精品久久| 日韩制服骚丝袜av| 国产一级毛片在线| 丝袜人妻中文字幕| 欧美日韩av久久| 国产成人精品在线电影| 亚洲精品久久午夜乱码| 国产成人午夜福利电影在线观看| 久久久国产一区二区| 久久久久久人妻| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 亚洲成av片中文字幕在线观看 | 人人澡人人妻人| 国产精品久久久av美女十八| 五月伊人婷婷丁香| 亚洲五月色婷婷综合| 美女脱内裤让男人舔精品视频| 国产精品女同一区二区软件| 国产又爽黄色视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇熟女欧美另类| 国产成人精品无人区| 毛片一级片免费看久久久久| 久久精品国产亚洲av涩爱| 一边亲一边摸免费视频| 欧美精品高潮呻吟av久久| 日本午夜av视频| 黄色视频在线播放观看不卡| 亚洲av日韩在线播放| 国精品久久久久久国模美| 少妇的逼水好多| 最近最新中文字幕免费大全7| 欧美3d第一页| 2022亚洲国产成人精品| av不卡在线播放| 女性生殖器流出的白浆| av国产精品久久久久影院| 秋霞伦理黄片| 2018国产大陆天天弄谢| 国产精品久久久av美女十八| 七月丁香在线播放| 久久久久久人妻| 大香蕉久久成人网| 一级黄片播放器| 日本-黄色视频高清免费观看| 国精品久久久久久国模美| 国产不卡av网站在线观看| 中国三级夫妇交换| 免费观看在线日韩| 一区二区三区乱码不卡18| 亚洲国产日韩一区二区| 校园人妻丝袜中文字幕| 高清视频免费观看一区二区| 中文字幕最新亚洲高清| 80岁老熟妇乱子伦牲交| 精品久久国产蜜桃| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 欧美 日韩 精品 国产| 亚洲精品日韩在线中文字幕| 狂野欧美激情性bbbbbb| 国产永久视频网站| 人妻少妇偷人精品九色| 亚洲国产精品一区二区三区在线| 亚洲综合精品二区| 波野结衣二区三区在线| 亚洲成人手机| 亚洲经典国产精华液单| 免费观看无遮挡的男女| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 亚洲三级黄色毛片| xxx大片免费视频| 久久久久久人人人人人| 99香蕉大伊视频| 欧美成人午夜精品| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说| 看十八女毛片水多多多| 美女福利国产在线| www.av在线官网国产| 久久综合国产亚洲精品| 欧美+日韩+精品| 91久久精品国产一区二区三区| 欧美97在线视频| 十分钟在线观看高清视频www| 三上悠亚av全集在线观看| 中文欧美无线码| 亚洲精品成人av观看孕妇| 国产女主播在线喷水免费视频网站| 午夜福利影视在线免费观看| 全区人妻精品视频| 春色校园在线视频观看| 日本免费在线观看一区| 人妻人人澡人人爽人人| 国产精品久久久久久久电影| 中文字幕精品免费在线观看视频 | 亚洲国产色片| 亚洲精品aⅴ在线观看| 中文字幕精品免费在线观看视频 | 久久精品夜色国产| 男女边摸边吃奶| 女性被躁到高潮视频| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| 国产不卡av网站在线观看| √禁漫天堂资源中文www| 国产xxxxx性猛交| av片东京热男人的天堂| 91成人精品电影| 精品一区二区免费观看| 精品亚洲成a人片在线观看| 人人妻人人澡人人看| 女的被弄到高潮叫床怎么办| 精品国产乱码久久久久久小说| 国产一级毛片在线| 欧美bdsm另类| 久久精品国产亚洲av天美| 少妇被粗大猛烈的视频| 亚洲精品国产av蜜桃| 日本色播在线视频| 天天躁夜夜躁狠狠躁躁| 精品久久久久久电影网| 伊人久久国产一区二区| 99热网站在线观看| 丁香六月天网| 黑人高潮一二区| av在线app专区| 国产成人91sexporn| 黑人巨大精品欧美一区二区蜜桃 | 男人添女人高潮全过程视频| 啦啦啦在线观看免费高清www| 在线观看一区二区三区激情| 啦啦啦在线观看免费高清www| 巨乳人妻的诱惑在线观看| 免费不卡的大黄色大毛片视频在线观看| 新久久久久国产一级毛片| 亚洲综合色惰| 天天影视国产精品| 精品久久国产蜜桃| 97超碰精品成人国产| 一级毛片电影观看| 日韩成人伦理影院| 中国国产av一级| tube8黄色片| 26uuu在线亚洲综合色| 国内精品宾馆在线| 欧美少妇被猛烈插入视频| 熟女人妻精品中文字幕| 国产伦理片在线播放av一区| 极品人妻少妇av视频| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀 | 久久人人爽av亚洲精品天堂| 丁香六月天网| 日本色播在线视频| 亚洲欧美中文字幕日韩二区| 国产白丝娇喘喷水9色精品| 狂野欧美激情性xxxx在线观看| 丁香六月天网| 国产一区二区在线观看日韩| 少妇 在线观看| 亚洲人成77777在线视频| 在线观看国产h片| 国产免费又黄又爽又色| 日韩伦理黄色片| 欧美 日韩 精品 国产| 精品国产乱码久久久久久小说| 内地一区二区视频在线| 亚洲一级一片aⅴ在线观看| 午夜影院在线不卡| 欧美少妇被猛烈插入视频| 国产精品人妻久久久久久| 黑人欧美特级aaaaaa片| 黄色配什么色好看| 最后的刺客免费高清国语| 中文字幕人妻丝袜制服| 男人爽女人下面视频在线观看| 精品人妻偷拍中文字幕| 天堂俺去俺来也www色官网| 久久99蜜桃精品久久| 亚洲精品456在线播放app| 日韩伦理黄色片| 国产乱人偷精品视频| 99热网站在线观看| 日韩电影二区| 两个人看的免费小视频| 性高湖久久久久久久久免费观看| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品一区二区三区在线| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 晚上一个人看的免费电影| 草草在线视频免费看| 国产亚洲精品第一综合不卡 | 人妻一区二区av| 国产成人av激情在线播放| 在线天堂中文资源库| 国产精品久久久久久av不卡| 在线天堂最新版资源| 国产老妇伦熟女老妇高清| 91午夜精品亚洲一区二区三区| 国产69精品久久久久777片| 少妇被粗大猛烈的视频| 高清在线视频一区二区三区| 一边亲一边摸免费视频| 国产成人一区二区在线| 亚洲高清免费不卡视频| 纵有疾风起免费观看全集完整版| 国产亚洲一区二区精品| 欧美bdsm另类| 日韩制服丝袜自拍偷拍| 亚洲国产欧美在线一区| 久久亚洲国产成人精品v| 国精品久久久久久国模美| 大陆偷拍与自拍| 国产精品嫩草影院av在线观看| 91午夜精品亚洲一区二区三区| 啦啦啦啦在线视频资源| 亚洲伊人久久精品综合| 精品国产一区二区三区四区第35| 成人18禁高潮啪啪吃奶动态图| 免费观看av网站的网址| 最后的刺客免费高清国语| 国产女主播在线喷水免费视频网站| 亚洲五月色婷婷综合| 哪个播放器可以免费观看大片| 免费黄网站久久成人精品| 18禁裸乳无遮挡动漫免费视频| 新久久久久国产一级毛片| 国产一区二区在线观看av| 日韩制服丝袜自拍偷拍| 国产av码专区亚洲av| 亚洲国产精品专区欧美| 国产精品国产三级国产专区5o| 久久99精品国语久久久| 三级国产精品片| 99re6热这里在线精品视频| 亚洲 欧美一区二区三区| 亚洲精品美女久久av网站| 亚洲成国产人片在线观看| 飞空精品影院首页| 国产一区二区在线观看日韩| 国产深夜福利视频在线观看| 国产国语露脸激情在线看| 国产深夜福利视频在线观看| 国产国语露脸激情在线看| 中文字幕免费在线视频6| 日本午夜av视频| 18禁动态无遮挡网站| 国产欧美日韩综合在线一区二区| 亚洲国产欧美日韩在线播放| 成人手机av| 久久99热6这里只有精品| 成人毛片a级毛片在线播放| 一本久久精品| 欧美最新免费一区二区三区| 91精品国产国语对白视频| 国产在线视频一区二区| 午夜免费鲁丝| 精品少妇黑人巨大在线播放| 精品第一国产精品| 亚洲,欧美精品.| 久久热在线av| 97超碰精品成人国产| 搡老乐熟女国产| 久久av网站| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 成人影院久久| 日韩成人伦理影院| 9色porny在线观看| 中国美白少妇内射xxxbb| 国产日韩一区二区三区精品不卡| 尾随美女入室| 综合色丁香网| www.av在线官网国产| 国产免费视频播放在线视频| 成人亚洲精品一区在线观看| xxxhd国产人妻xxx| 国产一区二区激情短视频 | 各种免费的搞黄视频| 两个人免费观看高清视频| 亚洲精品av麻豆狂野| 欧美xxⅹ黑人| 美女视频免费永久观看网站| 日韩三级伦理在线观看| 欧美人与善性xxx| 高清欧美精品videossex| 大香蕉久久成人网| 午夜福利视频在线观看免费| 国产精品女同一区二区软件| 免费观看av网站的网址| 热re99久久精品国产66热6| 色吧在线观看| 久久av网站| 国产在视频线精品| 91aial.com中文字幕在线观看| 精品亚洲成a人片在线观看| 国产高清国产精品国产三级| 性高湖久久久久久久久免费观看| 午夜激情av网站| 又大又黄又爽视频免费| 80岁老熟妇乱子伦牲交| 亚洲婷婷狠狠爱综合网| 亚洲av国产av综合av卡| 精品亚洲成a人片在线观看| 中国国产av一级| 九九在线视频观看精品| 国产69精品久久久久777片| 中国三级夫妇交换| 亚洲在久久综合| 老司机影院成人| 婷婷色综合大香蕉| 亚洲精品久久成人aⅴ小说| 一二三四中文在线观看免费高清| 国产精品 国内视频| 日韩制服丝袜自拍偷拍| 国产片内射在线| 大陆偷拍与自拍| 国产精品一二三区在线看| 一本—道久久a久久精品蜜桃钙片| 亚洲精华国产精华液的使用体验| 人妻人人澡人人爽人人| 两性夫妻黄色片 | 日本wwww免费看| 欧美变态另类bdsm刘玥| 咕卡用的链子| 如日韩欧美国产精品一区二区三区| 天天影视国产精品| 亚洲,一卡二卡三卡| 欧美精品一区二区大全| 亚洲图色成人| 中国国产av一级| 久久人妻熟女aⅴ| 国产精品蜜桃在线观看| 多毛熟女@视频| 久久国产精品大桥未久av| 国产精品一二三区在线看| 各种免费的搞黄视频| 欧美日本中文国产一区发布| 97人妻天天添夜夜摸| 少妇人妻 视频| av视频免费观看在线观看| 丰满少妇做爰视频| 街头女战士在线观看网站| 国产xxxxx性猛交| 国产免费视频播放在线视频| 婷婷色麻豆天堂久久| av在线老鸭窝| 看免费av毛片| 成人毛片a级毛片在线播放| 久久精品国产亚洲av天美| 亚洲av电影在线观看一区二区三区| 免费女性裸体啪啪无遮挡网站| 如何舔出高潮| 丰满乱子伦码专区| 国产白丝娇喘喷水9色精品| 麻豆乱淫一区二区| 777米奇影视久久| 国产成人精品福利久久| 久久人人爽人人片av| 色吧在线观看| 黄色配什么色好看| 亚洲欧美日韩另类电影网站| 国产精品国产三级国产av玫瑰| 欧美日本中文国产一区发布| 最近的中文字幕免费完整| 成人国语在线视频| 黄色一级大片看看| 欧美成人午夜精品| 欧美精品国产亚洲| 亚洲人成77777在线视频| 热re99久久精品国产66热6| 九色成人免费人妻av| 欧美日韩综合久久久久久| 国产成人欧美| 97在线视频观看| 免费女性裸体啪啪无遮挡网站| 久久久久久久国产电影| 青春草亚洲视频在线观看| 亚洲国产精品一区二区三区在线| 欧美成人午夜免费资源| 久久久久人妻精品一区果冻| 亚洲av日韩在线播放| 99久久人妻综合| 日本av手机在线免费观看| 亚洲精品乱码久久久久久按摩| 精品亚洲成国产av| 一二三四中文在线观看免费高清|