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

    青藏高原地區(qū)山生柳遺傳多樣性研究

    2019-09-10 07:22:44吳錦楊茜王久利高慶波張發(fā)起陳世龍
    廣西植物 2019年9期
    關(guān)鍵詞:遺傳多樣性青藏高原

    吳錦 楊茜 王久利 高慶波 張發(fā)起 陳世龍

    摘要:? 利用葉綠體非編碼區(qū)片段研究分布于青藏高原地區(qū)的山生柳居群遺傳多樣性,對未來山生柳生態(tài)環(huán)境和青藏高原地區(qū)物種豐富度的保護(hù)具有指導(dǎo)意義。該研究設(shè)計(jì)并篩選出cpDNA引物5′trnG2G-3′trnG(UUC)和5′rpS12-rpL20,用擴(kuò)增出的片段和對應(yīng)的聯(lián)合片段進(jìn)行后續(xù)的遺傳多樣性分析。結(jié)果表明:通過山生柳的聯(lián)合片段檢測到35種單倍型,單倍型多態(tài)性0.626,核苷酸多態(tài)性0.000 85。中性檢驗(yàn)Tajima’s D(-2.286 70, P<0.01)和Fu’s Fs(-5.298 05, P<0.02)都是顯著負(fù)值,推測山生柳個(gè)體數(shù)近期經(jīng)歷過擴(kuò)張。AMOVA分析顯示,居群內(nèi)和居群間遺傳變異分別為93.70%和6.30%,表明居群內(nèi)的變異是山生柳遺傳變異的主要來源。居群間遺傳分化程度中等偏低(FST=0.063),基因流(Nm)為7.439,說明山生柳各居群的基因交流非常頻繁,不同地理居群間存在一定的基因流動(dòng)。遺傳分化系數(shù)NST(0.075)大于GST(0.068)和基于遺傳距離和單倍型的UPGMA聚類分析,表明山生柳12個(gè)居群分為4組且與居群的地理分布沒有明顯相關(guān)性。山生柳是進(jìn)行有性繁殖還是無性繁殖主要受環(huán)境因素的影響,居群內(nèi)變異是山生柳遺傳變異性的主要來源,居群間基因交流頻繁。

    關(guān)鍵詞: 山生柳, 基因流, 青藏高原, 葉綠體DNA, 遺傳多樣性

    中圖分類號:?; Q943文獻(xiàn)標(biāo)識碼:? A文章編號:? 1000-3142(2019)09-1180-11

    Abstract:? We studied the genetic diversity of Salix oritrepha in the Qinghai-Tibetan Plateau using the non-coding region of chloroplast DNA fragments. The study will have? guiding significance for the protection of ecological environment and animal husbandry of S. oritrepha and species richness in Qinghai-Tibet Plateau in the future. We designed two useful primers 5′trnG2G-3′trnG (UUC) and 5′rpS12-rpL20. The amplified fragments and their combined fragments were used for subsequent genetic diversity analysis of S. oritrepha. The results showed that we detected 35 haplotypes by the combined sequences, haplotype diversity (Hd) was 0.626, nucleotide diversity (Pi) was 0.000 85. Neutrality tests of Tajima’s D (-2.286 70, P<0.01) and Fus Fs (-5.298 05, P<0.02) showed significantly negative values, which meant that S. oritrepha had expansion recently of its individual numbers. Analysis of molecular variance showed that within-population variance was described as 93.70% of the total variance, while among-population variance was 6.30%, which meant that within-population variance was the main resource of genetic variation of S. oritrepha. The fi-xation index (FST) was 0.063 and number of migrants (Nm) was 7.439, which meant that gene flow was very frequent among different geographical populations. An estimation of non-signification NST (0.075) was bigger than gene differentiation coefficient GST (0.068). UPGMA tree based on genetic distance and haplotypes suggested that all populations were divided into four groups and had no significant correlation with the geographical distribution of the population. The mode of reproduction (sexual or asexual) of S. oritrepha mainly depended on environmental factors. The genetic variability mainly came from within-population and it has a frequently gene flow among populations.

    Key words:? Salix oritrepha, gene flow, Qinghai-Tibetan Plateau, chloroplast DNA, genetic diversity

    山生柳(Salix oritrepha)是中國特有物種,落葉矮小灌木,高10~120 cm,藏藥中稱朗瑪,生于海拔3 000~4 700 m的高山地帶,繁衍方式以有性生殖為主,無性繁殖為輔,是青藏高原地區(qū)高寒柳灌叢的建群種和優(yōu)勢種,具有極高的藥用價(jià)值,對生態(tài)環(huán)境的保護(hù)也尤為主要。一方面對山生柳化學(xué)成分的研究,證明其可解熱、消炎止痛、治療婦科疾病和提高機(jī)體免疫力(封士蘭等,2001);另一方面山生柳是深根性灌木,根深長度達(dá)到植株縱向高度的4 ~ 5倍,故可以抵抗干旱和低溫。由于耐高鹽高堿且生長繁殖快速,故能涵養(yǎng)水源、治沙和防沙,使其起到河岸防護(hù)、水土保持的作用(陳文業(yè)等,2008)。常與高山繡線菊(Spiraea alpina)、金露梅(Potentilla fruticosa)、沙棘(Hippophae rhamnoides)、窄葉鮮卑花(Sibiraea angustata)、大苞柳(Salix pseudospissa)等組成高寒灌叢(王芳等,2010),是藏民牲畜的夏季放牧地。

    遺傳多樣性分析是評價(jià)和保護(hù)瀕危物種的主要指標(biāo)(Cires et al., 2011),種內(nèi)遺傳多樣性越豐富,物種對環(huán)境的適應(yīng)能力就越強(qiáng),保護(hù)了物種和生態(tài)系統(tǒng)的多樣化,也能減慢因?yàn)檫m應(yīng)、進(jìn)化所致的滅亡進(jìn)程。近幾十年來,對青藏高原地區(qū)植被的遺傳多樣性的研究越來越深入(Duan et al., 2011; Khan et al., 2014),這些都對保護(hù)青藏高原地區(qū)植被豐富度具有重要的指導(dǎo)意義。研究植物遺傳多樣性的方法中,葉綠體基因組DNA被廣泛使用(Demesure et al., 1996;Rungis et al., 2017;Wu, 2016)。葉綠體DNA的結(jié)構(gòu)緊湊簡單、分子量小,可從分子角度作比較精準(zhǔn)的分析,其比核DNA存在更加穩(wěn)定、進(jìn)化更加緩慢且編碼序列也保守,故廣泛應(yīng)用于植物遺傳多樣性的剖析研究(Shaw et al., 2014)。葉綠體DNA在大多數(shù)被子植物中為母系遺傳,只通過種子進(jìn)行遺傳,不受花粉遺傳的影響(Mogensen, 1996)。

    目前,對高寒山生柳灌叢的群落特征、在畜牧業(yè)中的作用以及在防止沙漠化中起到的作用都有研究(Kang et al., 2015;王緯和王盛,2015;楊路存,2013),但對山生柳遺傳多樣性的研究較少,郭敏等(2012)的研究表明山生柳種群具有較高的遺傳多樣性且隨著海拔梯度的升高而呈現(xiàn)升高趨勢,且遺傳變異主要存在于居群內(nèi),因其所研究的地理分布范圍僅限祁連山,不能準(zhǔn)確分析出山生柳遺傳結(jié)構(gòu)及種群動(dòng)態(tài)。本研究采集了青藏高原大部分地區(qū)的山生柳,設(shè)計(jì)并篩選出兩對特異性針對山生柳的cpDNA引物,用葉綠體非編碼區(qū)的聯(lián)合片段對山生柳進(jìn)行遺傳多樣性研究,將對山生柳生殖系統(tǒng)的研究、種質(zhì)資源的保護(hù)和生態(tài)環(huán)境的保護(hù)都有重要意義。

    1材料與方法

    1.1 材料

    2006—2017年對分布于青藏高原地區(qū)的山生柳進(jìn)行廣泛的野外采集。對每個(gè)野外自然居群選取8~58個(gè)個(gè)體,個(gè)體間距超過100 m,共收集山生柳12個(gè)居群351個(gè)個(gè)體(表1)。采集山生柳生長良好的新鮮幼嫩葉片,迅速用硅膠進(jìn)行干燥,之后放于-20 ℃ 長期保存。所采集的憑證標(biāo)本存放于中國科學(xué)院西北高原生物研究所青藏高原生物標(biāo)本館(HNWP)中。

    1.2 DNA提取與PCR擴(kuò)增

    選取改良的CTAB法從干燥處理的山生柳嫩葉中提取基因組DNA(Porebski et al., 1997)。參照序列選取已上傳到NCBI上的Salix arbutifolia、S. babylonica、S. interior、S. oreinoma、S. purpurea、S. suchowensis和S. tetrasperma等柳屬植物葉綠體基因組序列,引物片段的選取參考之前的研究(Hamilton, 1999; Scarcelli et al., 2011; Shaw et al., 2005),用軟件Primer Premier 5.0(Singh et al., 1998)、BioXM 2.6(黃驥和張紅生,2004)和SnapGene 3.2.1(McKeone et al., 2014)進(jìn)行引物設(shè)計(jì),共設(shè)計(jì)了21對引物用于PCR擴(kuò)增。擴(kuò)增反應(yīng)體系為25 μL:2.5 μL的10×PCR Buffer(含1.5 mmol·L-1 MgCl2)、0.5 μL的10 mmol·L-1 dNTP、正反引物各0.5 μL(10 pmol·L-1)、1個(gè)單位的Taq DNA酶(TaKaRa,大連)、15~50 ng的模板DNA,去離子水補(bǔ)齊到25 μL。PCR反應(yīng)擴(kuò)增程序:94 ℃預(yù)變性3 min;94 ℃變性45 s;60 ℃退火45 s;72 ℃延伸30 s;35個(gè)循環(huán);72 ℃延伸5 min;4 ℃保存。PCR反應(yīng)產(chǎn)物送至生物工程(上海)股份有限公司測序。將獲得21對引物在12個(gè)居群24個(gè)個(gè)體中進(jìn)行PCR擴(kuò)增預(yù)試驗(yàn),篩選出獲得兩對遺傳多態(tài)性較高的引物 [5′trnG2G-3′trnG(UUC)和5′rpS12-rpL20]用于所有個(gè)體的擴(kuò)增測序(表2)。

    1.3 數(shù)據(jù)分析

    運(yùn)用Chromas軟件(http://technelysium.com.au/wp/chromas/)和MEGA 7.0.26(Kumar et al., 2018)對所有測序序列進(jìn)行人工校對,并用MEGA軟件對校對好的所有個(gè)體的擴(kuò)增片段進(jìn)行對位排列。利用DnaSP 6.10.04(Rozas et al., 2017)軟件確定山生柳單倍型,并統(tǒng)計(jì)中性檢驗(yàn)中的Tajima’s D和Fu’s Fs。運(yùn)用SPSS(George & Mallery, 2003)軟件獲得基于遺傳距離和單倍型分布的UPGMA聚類分析圖。利用Arlequin軟件包(Excoffier & Lischer, 2010)中的分子變異分析方法(analysis of molecular variance, AMOVA)檢測居群內(nèi)和居群間的遺傳變異水平,并對單倍型分布的FST進(jìn)行評價(jià)(1 000次重復(fù)置換檢測)。根據(jù)基因流(Nm)與遺傳分化系數(shù)(FST)之間的關(guān)系計(jì)算Nm值和居群平均等位基因數(shù)。失配分布(mismatch distribution)分析也在Arlequin中完成。通過PERMUT軟件計(jì)算居群內(nèi)平均遺傳多樣性(average gene diversity within populations,hS)、總遺傳多樣性(total gene diversity,hT)、居群間遺傳分化系數(shù)(GST)和(NST),并用U-統(tǒng)計(jì)方法對GST和NST進(jìn)行比較(1 000次重復(fù)置換檢測)以檢測山生柳單倍型的分布特點(diǎn)(Caraux & Pinloche, 2005)。利用TCS軟件(Clement et al., 2000)來建立山生柳單倍型之間的TCS簡約網(wǎng)絡(luò)關(guān)系圖。

    2結(jié)果與分析

    2.1 山生柳遺傳多樣性分析

    對351個(gè)個(gè)體的擴(kuò)增產(chǎn)物進(jìn)行單向測序,利用MEGA軟件對位排列且刪除poly結(jié)構(gòu)后,5′trnG2G-3′trnG(UUC)和5′rpS12-rpL20片段矩陣長度分別為481 bp、730 bp,共檢測到36個(gè)變異位點(diǎn),其中17個(gè)為簡約性信息位點(diǎn),變異率為2.97%。從山生柳的聯(lián)合片段中檢測出35個(gè)單倍型(Haplotype, H01~H35),將5′trnG2G-3′trnG(UUC)片段的單倍型序列和5′rpS12-rpL20片段的單倍型序列提交至GeneBank(MK602662~MK602696和MK606373~MK606407)。聯(lián)合片段的單倍型多樣性為0.626 1,單倍型多樣性標(biāo)準(zhǔn)差為0.029(表3)。各居群聯(lián)合片段之間的多態(tài)性表現(xiàn)不同,居群的平均等位基因數(shù)(Na)依次為1.067、1.078、1.278、1.022、1.656、1.089、1.067、1.100、1.056、1.078、1.089、1.033、1.134。遺傳多樣性最豐富的是居群YS,其具有最高的核苷酸多樣性(0.001 45)、平均核苷酸差異(1.722)和最多的單倍型數(shù)目(13)(表4)。

    由各居群的單倍型分布可知,僅H01被全部12個(gè)居群所共享;H04、H05、H06和H07是GZ所擁有的特有單倍型;H11是MQ所擁有的特有單倍型;H17、H18、H19、H20、H21和H32是YS所擁有的特有單倍型;H27、H28、H29和H30是GD所擁有的特有單倍型;H16是MY所擁有的特有單倍型;H15是DR所擁有的特有單倍型;H24、H25和H26是BM所擁有的特有單倍型;H31是NQ所擁有的特有單倍型;H33是QL所擁有的特有單倍型;H34和H35是SN所擁有的特有單倍型;H02,H03,H08,H09,H10,H12,H13,H14,H16,H17,H22,H23是2~10個(gè)居群所共有的單倍型(圖1,表4)。

    2.2 山生柳遺傳結(jié)構(gòu)分析

    分子變異分析(AMOVA)結(jié)果表明,居群間遺傳變異為6.30%,居群內(nèi)遺傳變異為93.70%(P<0.01),說明居群內(nèi)遺傳變異為山生柳遺傳變異的主要來源(表5)。遺傳分化系數(shù)FST=0.062 98,居群間遺傳分化很小。根據(jù)基因流Nm與FST之間的關(guān)系(Kirk & Freeland, 2011),Nm =(1-FST)/2FST,得到Nm為7.439, 基因流很強(qiáng)。用PERMUT和DnaSP軟件計(jì)算得出山生柳居群內(nèi)平均遺傳多樣山生柳各居群間的遺傳距離大部分都較小,其中QM和NQ之間出現(xiàn)了最大遺傳距離0.977,YS和SN出現(xiàn)了最小遺傳距離0.000,大多數(shù)居群之間的遺傳距離差異性不大,總的平均遺傳距離為0.424;居群間的遺傳分化指數(shù)(FST)顯示有18個(gè)負(fù)值,占27.28%。其中QM和NQ之間出現(xiàn)了最大的FST=0.255 0(表6)。

    山生柳12個(gè)居群的遺傳距離的范圍為0.000~1.234,表明山生柳居群內(nèi)遺傳多樣性豐富?;谶z傳距離的UPGMA聚類分析,得到山生柳的12個(gè)居群分為I(YS, SN, MY, DT, QL, BM)、Ⅱ(GZ, MQ, DR)、Ⅲ(GD, NQ)和IV(QM)共四組(圖2)。基于單倍型的UPGMA聚類分析,得到12個(gè)居群分為I(YS, GD, MY)、II(NQ)、III(BM, GZ)和IV(QM, SN, QM, QL, DT, DR)共四組(圖3)。兩者均沒有表現(xiàn)出與12個(gè)居群的地理分布相對應(yīng)的聚類關(guān)系。

    2.3 山生柳種群動(dòng)態(tài)

    中性檢驗(yàn)結(jié)果表明,山生柳各個(gè)引物片段和聯(lián)合片段的Tajima’s D(-2.286 70, P<0.01)和Fu’s Fs(-5.298 05, P<0.02)均為顯著的負(fù)值(表3)。失配分析結(jié)果顯示,在山生柳分布范圍內(nèi)其分布曲線整體呈現(xiàn)單峰分布, 其SSD(sum of square devia-tions)值 (0.003 8, P=0.68)和HRag (Harpending’s raggedness index)值(0.047 9, P=0.72)均不顯著(P>0.05),這與中性檢驗(yàn)的結(jié)果相一致(圖4)。

    用山生柳葉綠體基因所檢測出的35個(gè)單倍型,建立單倍型多樣性的簡約網(wǎng)絡(luò)圖(圖5)?;芈方Y(jié)構(gòu)是趨同演化的結(jié)果,矩形框代表祖先(或根部)單倍型即H01(Saltonstall, 2002)。圖5中略去了單倍型的頻率數(shù)值。圖5中“o”表示理論分析出的,但沒有實(shí)際觀測到的單倍型。除YS(H32)和DG(H27,H29)的單倍型同其它單倍型關(guān)系較遠(yuǎn)外,大部分單倍型之間僅有一到兩步突變,演化關(guān)系呈現(xiàn)星狀輻射結(jié)構(gòu)(楊玉慧等,2004)。

    3討論

    3.1 山生柳的居群遺傳多樣性

    篩選出的山生柳特異性條帶為山生柳無性系鑒定提供可操作性。山生柳的聯(lián)合片段共檢測到36個(gè)變異位點(diǎn)和17個(gè)簡約性的信息位點(diǎn),插入缺失變異位點(diǎn)19個(gè),占總變異的52.78%。插入缺失是山生柳變異位點(diǎn)出現(xiàn)的主要來源,與DNA在進(jìn)行復(fù)制時(shí)存在高頻率的點(diǎn)突變有關(guān)。在殘酷的高山圖 4基于山生柳單倍型的失配分析山地冰雪寒溫帶生存環(huán)境中,山生柳的DNA復(fù)制發(fā)生堿基錯(cuò)配的可能性更大,機(jī)體的自我修復(fù)功能會受到一定阻礙,說明惡劣的生存環(huán)境可能在一定程度上影響DNA復(fù)制過程(Baccarelli & Bollati, 2009; Chinnusamy & Zhu, 2009; Dubrova, 2003)。

    對山生柳居群的聯(lián)合片段進(jìn)行數(shù)據(jù)分析,檢測到35個(gè)單倍型。數(shù)據(jù)分析的結(jié)果總體呈現(xiàn)出較低的遺傳多樣性(Hd=0.626 1, Pi=0.000 85)。山生柳每個(gè)居群間的單倍型多樣性(Hd)表現(xiàn)不同,推測可能是由于在青藏高原地區(qū)的山生柳分布范圍廣,各居群所在的自然生存環(huán)境條件具有差異,經(jīng)過長期不同的生態(tài)環(huán)境和自然選擇,會導(dǎo)致具有一定地理隔離的山生柳居群存在差異的遺傳多樣性。根據(jù)之前對山生柳在祁連山上不同海拔種群的遺傳多樣性水平的研究結(jié)論(郭敏等,2012),推測不同的地理位置、生存環(huán)境和海拔高度都體現(xiàn)山生柳不顯著的分子水平的差異性,但是差異性不會太大。

    3.2 山生柳的居群遺傳結(jié)構(gòu)

    AMOVA分子變異結(jié)果表明,山生柳居群內(nèi)的遺傳變異占總變異的93.70%(P<0.01),居群間的遺傳變異僅為6.30%,揭示了山生柳居群內(nèi)的遺傳變異比居群間的遺傳變異要大,居群內(nèi)的遺傳變異是山生柳遺傳變異的主要來源。推測的原因:(1)35個(gè)單倍型中僅有H01為12個(gè)居群都擁有的共享單倍型,共有208個(gè)個(gè)體擁有,所占比例為59.26%,由于H01在居群中是無規(guī)律分布,因此山生柳居群間的遺傳差異并不明顯;(2)35個(gè)單倍型中有22個(gè)單倍型是各個(gè)居群所獨(dú)有,甚至某一居群會同時(shí)固定多個(gè)特有單倍型,這就說明山生柳的某些居群內(nèi)部會有可辨別的遺傳變異。特有單倍型比例較高的物種往往不能揭示出明顯的譜系地理格局(Geng et al., 2018)。繁育系統(tǒng)和遺傳變異分布之間存在著密切的聯(lián)系(Hogbin & Peakall, 1999),即異交物種的遺傳變異大多分布在群體之內(nèi),群體間的遺傳變異通常占27%以下(趙景奎等,2008),與本研究結(jié)果相符。

    植物居群的遺傳結(jié)構(gòu)受群體進(jìn)化歷史、基因流等影響,群體進(jìn)化歷史一般通過遺傳分化指數(shù)(FST)來反映,其大小可在一定程度上揭示種群間基因流和遺傳漂變的程度(楊艷婷等,2018)。一般來講,當(dāng)FST∈[0,0.05),群體間沒有遺傳分化;當(dāng)FST∈[0.05,0.15),群體分化程度中等;當(dāng)FST∈[0.15,0.25),群體間高度分化;FST∈[0.25,1],群體間分化程度非常高(Wright,1950,1970)。山生柳自然居群的遺傳分化指數(shù)(FST)為0.062 98,表明居群分化程度為中等。山生柳各居群間的遺傳分化指數(shù)(FST)有19個(gè)負(fù)值,占28.79%(正常情況下FST∈[0,1]),表明少部分居群間沒有分化,但總體上反映出的是大部分居群間出現(xiàn)非常普遍的基因交流(Wright, 1950, 1970, 1982),其中QM和NQ之間的FST值最大為0.255 0。山生柳一般為有性生殖,花粉量大隨風(fēng)飄散,但在嚴(yán)酷的生存環(huán)境下會進(jìn)行無性繁殖,本研究所選的山生柳材料地理分布較廣、地理環(huán)境高度復(fù)雜,導(dǎo)致居群間遺傳變異的差異性較大,有的居群幾乎無遺傳分化(QM, DT),有的居群分化程度非常高(QL),但總體還是表現(xiàn)為分化程度中等偏低。

    根據(jù)基因流公式,算得Nm=7.438,當(dāng)Nm=7.438>1,屬于較高水平(Wright,1982),表明不同的居群之間存在頻繁的基因流動(dòng),增大群體間遺傳差異的表現(xiàn)一般為自然選擇以及遺傳漂變,而基因流可以“弱化”群體之間不易察覺的遺傳差異使居群間趨向統(tǒng)一。

    對山生柳單倍型變異的地理結(jié)構(gòu)檢驗(yàn)結(jié)果為GST值(0.068)和NST值(0.075)(NST>GST,P>0.05),表明山生柳不同居群間出現(xiàn)了大量的親緣關(guān)系靠近的單倍型,可以證明山生柳不存在明顯的譜系地理學(xué)的關(guān)系?;谶z傳距離和單倍型的UPGMA聚類分析,得到山生柳的12個(gè)居群各分為四組且與12個(gè)居群地理分布沒有明顯相關(guān)性,進(jìn)一步驗(yàn)證了山生柳居群不具有譜系地理結(jié)構(gòu)。

    3.3 山生柳的種群擴(kuò)張

    基因流分析結(jié)果證實(shí)山生柳的12個(gè)居群存在頻繁的基因交流,山生柳居群間不存在明顯的基因流限制。通過對所有地理居群單倍型進(jìn)行失配分析,所有的山生柳單倍型呈現(xiàn)只有一個(gè)明顯頂峰的單峰型曲線,表示不能拒絕居群擴(kuò)張的假說。中性檢驗(yàn)的結(jié)果表明,Tajima’s D (-2.286 70, P<0.01)和Fu’s Fs(-5.298 05, P<0.02)均為顯著的負(fù)值,這與失配分析的單峰型曲線結(jié)果相一致。TCS網(wǎng)絡(luò)圖顯示H01單倍型為祖先單倍型,由于H01為種群的共享單倍型,故不能確定山生柳種群的擴(kuò)散中心。TCS網(wǎng)絡(luò)圖呈星狀輻射分布,且新衍生出的單倍型與原始單倍型組成一個(gè)“星狀”輻射結(jié)構(gòu),說明該群體曾經(jīng)在遇到瓶頸后迅速擴(kuò)增(Arbogast, 2000)。通過單倍型多樣性與核苷酸多樣性分析得到單倍型多樣性高(Hd>0.5)而核苷酸多樣性低(Pi<0.005),預(yù)示著群體曾經(jīng)歷過瓶頸效應(yīng)后,伴隨了迅速的種群擴(kuò)張與變異的積累。綜上所述,推測山生柳在其分布范圍或者其個(gè)體數(shù)近期經(jīng)歷過擴(kuò)張。

    對青藏高原地區(qū)相同生境的高山繡線菊(Khan et al., 2014)、金露梅(Ma et al., 2014)、沙棘(Bartish et al., 2006)、窄葉鮮卑花(Duan et al., 2011)等灌木的種群動(dòng)態(tài)研究表明其遺傳分組曾經(jīng)歷過快速的種群擴(kuò)張事件。這與山生柳的種群擴(kuò)張結(jié)果一致,推測物種生境和繁育方式的相似性會出現(xiàn)物種為了營養(yǎng)、繁殖、保護(hù)、擴(kuò)散以及對環(huán)境的適應(yīng)出現(xiàn)協(xié)同進(jìn)化的現(xiàn)象(Morgan et al., 2005; Thompson & Burdon, 1992)。有研究表明山生柳種群可能以有性和無性兩種繁殖方式進(jìn)行更新以鞏固其優(yōu)勢地位(李毅等,2002),年齡結(jié)構(gòu)分布圖呈正三角形,種群動(dòng)態(tài)表現(xiàn)為增長型(李毅等,2002),這與本研究中居群近期出現(xiàn)過擴(kuò)張相吻合?;ǚ酆头N子的擴(kuò)散是自然植物種群最主要的基因流方式(李海生和陳桂珠,2004),故山生柳生殖方式的多樣化會導(dǎo)致其基因交流頻繁,預(yù)期未來山生柳居群會處于進(jìn)一步擴(kuò)張狀態(tài)。

    山生柳頻繁的基因交流和經(jīng)歷過近期擴(kuò)張,猜測此結(jié)果可能的原因如下:(1)cpDNA是母系遺傳,估算的基因流大小與花粉傳播無關(guān),與其多樣化的繁育方式有關(guān);(2)山生柳一般以有性生殖為主,其生存環(huán)境是長期積雪(低溫)、干燥風(fēng)大等氣象特點(diǎn),使其種子可以隨處傳播;(3)由于種子的傳播距離遠(yuǎn),故青藏高原地區(qū)和橫斷山脈錯(cuò)綜復(fù)雜地形也只能在很小的程度上阻隔山生柳居群間的基因交流;(4)在嚴(yán)酷環(huán)境壓力下,山生柳種子活力極低,幼苗補(bǔ)充量受限,從而有性生殖受限,此時(shí)主要依托無性繁殖來鞏固其物種數(shù)量的穩(wěn)定(陳文業(yè)等,2008)。

    參考文獻(xiàn):

    ARBOGAST BS, 2000. Phylogeography: The history and formation of species [J]. Integr Comp Biol, 41(1): 134-135.

    BACCARELLI A, BOLLATI V, 2009. Epigenetics and environmental chemicals [J]. Curr Opin Pediatr, 21(2): 243-251.

    BARTISH IV, KADEREIT JW, COMES HP, 2006. Late Quaternary history of Hippophae rhamnoides L. (Elaeagnaceae) inferred from chalcone synthase intron (Chsi) sequences and chloroplast DNA variation [J]. Mol Ecol, 15(13): 4065-4083.

    CARAUX G, PINLOCHE S, 2005. PermutMatrix: A graphical environment to arrange gene expression profiles in optimal linear order [J]. Bioinformatics, 21(1): 1280-1281.

    CHEN WY, CHEN YX, ZHENG HP, et al., 2008. Ecological habits and seedling raising techniques of Salix oritrepha [J]. Prot For Sci Technol, (3): 135-137.? [陳文業(yè), 陳永霞, 鄭華平, 等, 2008. 山生柳的生態(tài)習(xí)性與育苗技術(shù) [J]. 防護(hù)林科技, (3):135-137.]

    CHINNUSAMY V, ZHU JK, 2009. Epigenetic regulation of stress responses in plants [J]. Curr Opin Plant Biol, 12(2): 133-139.

    CIRES E, SAMAIN MS, GOETGHEBEUR P, et al., 2011. Genetic structure in peripheral Western European populations of the endangered species Cochlearia pyrenaica (Brassicaceae) [J]. Plant Syst Evol, 297(1-2): 75-85.

    CLEMENT M, POSADA D, CRANDALL KA, 2000. TCS: A computer program to estimate gene genealogies [J]. Mol Ecol, 9(10): 1657-1659.

    DEMESURE B, COMPS B, PETIT RJ, 1996. Chloroplast DNA phylogeography of the common beech (Fagus sylvatica L.) in Europe [J]. Evolution, 50(6): 2515.

    DUAN YZ, GAO QB, ZHANG FQ, et al., 2011. Phylogeographic analysis of the endemic species Sibiraea angustata reveals a marginal refugium in the Qinghai-Tibet Plateau [J]. Nord J Bot, 29(5): 615-624.

    DUBROVA YE, 2003. Radiation-induced transgenerational instability [J]. Oncogene, 22(45): 7087-7093.

    EXCOFFIER L, LISCHER H, 2010. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows [J]. Mol Ecol Resour, 10(3): 564-567.

    FENG SL, SU C, DUAN SY, et al., 2001. Chemical constituents from Salix oritrepha [J]. Chin J Chin Mat Med, (9): 31-32.? [封士蘭, 蘇策, 段生玉, 等, 2001. 藏藥山生柳的化學(xué)成分研究 [J]. 中國中藥雜志, (9):31-32.]

    GENG JZM, LI Y, JIA LK, et al., 2018. Phylogeography of Saxifraga tangutica Engl. (Saxifragaceae)? [J]. Acta Bot Boreal-Occident Sin, 38(2): 370-380.? [更吉卓瑪, 李彥, 賈留坤, 等, 2018. 唐古特虎耳草譜系地理學(xué)研究 [J]. 西北植物學(xué)報(bào), 38(2):370-380.]

    GEORGE D, MALLERY P, 2003. SPSS for Windows step by step: A simple guide and reference, 16.0 update [J]. Comp Softw, (100): 357.

    GUO M, MA YJ, LI Y, 2012. Analysis of genetic diversity in Salix oritrepha from different altitudes in Gansu Qilian mountains using SSR markers [J]. Acta Pratac Sin, 21(5): 114-121.? [郭敏, 馬彥軍, 李毅, 2012. 祁連山不同海拔梯度山生柳遺傳多樣性的SSR分析 [J]. 草業(yè)學(xué)報(bào), 21(5):114-121.]

    HAMILTON MB, 1999. Four primer pairs for the amplification of chloroplast intergenic regions with intraspecific variation [J]. Mol Ecol, 8(3): 521-523.

    HOGBIN PM, PEAKALL R, 1999. Evaluation of the contribution of genetic research to the management of the endangered plant Zieria prostrata [J]. Conserv Biol, 13(3): 514-522.

    HUANG Y, ZHANG HS, 2004. Development of nucleotide sequence analysis software based on Windows [J]. Chin J Bioinform, 2(1): 13-17. [黃驥, 張紅生, 2004. 基于Windows的核酸序列分析軟件的開發(fā) [J]. 生物信息學(xué), 2(1):13-17.]

    KANG JJ, ZHAO WZ, ZHAO M, et al., 2015. The cutting reproduction technique of Salix oritrepha and its application on degraded grasslands restoration [J]. J Soil Water Conserv, 70(1): 45-53.

    KHAN G, ZHANG FQ, GAO QB, et al., 2014. Molecular phylogeography and intraspecific divergence of Spiraea alpina (Rosaceae) distributed in the Qinghai-Tibetan Plateau and adjacent regions inferred from nrDNA [J]. Biochem Syst Ecol, 57: 278-286.

    KIRK H, FREELAND J R, 2011. Applications and implications of neutral versus non-neutral markers in molecular ecology [J]. Int J Mol Sci, 12(6): 3966-3988.

    KUMAR S, STECHER G, LI M, et al., 2018. MEGA X: Molecular evolutionary genetics analysis across computing platforms [J]. Mol Biol Evol, 35(6): 1547-1549.

    LI HS, CHEN GZ, 2004. Genetic diversity of mangrove plant Sonneratia caseolaris in Hainan Island based on ISSR analysis [J]. Acta Ecol? Sin, 24(8): 1657-1663.? [李海生, 陳桂珠, 2004. 海南島紅樹植物海桑遺傳多樣性的ISSR分析 [J]. 生態(tài)學(xué)報(bào), 24(8):1657-1663.]

    LI Y, HU ZZ, WANG ZT, 2002. Populations distribution patterns of alpine Salix Oritrepha in East Qilian Mountains [J]. Acta Pratac Sin, (3): 48-54.? [李毅, 胡自治, 王志泰, 2002. 東祁連山高寒地區(qū)山生柳種群分布格局研究 [J]. 草業(yè)學(xué)報(bào), (3):48-54.]

    MA YZ, LI ZH, WANG X, et al., 2014. Phylogeography of the genus Dasiphora (Rosaceae) in the Qinghai-Tibetan Plateau: Divergence blurred by expansion [J]. Biol J Linn Soc, 111(4): 777-788.

    MCKEONE R, WIKSTROM M, KIEL C, et al., 2014. Assessing the correlation between mutant rhodopsin stability and the severity of retinitis pigmentosa [J]. Mol Vis, 20: 183-199.

    MOGENSEN HL, 1996. Invited special paper: The hows and whys of cytoplasmic inheritance in seed plants [J]. Am J Bot, 83(3): 383-404.

    MORGAN MT, WILSON WG, KNIGHT TM, 2005. Plant population dynamics, pollinator foraging, and the selection of self-fertilization [J]. Am Nat, 166(2): 169-183.

    POREBSKI S, BAILEY LG, BAUM BR, 1997. Modification of a CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components [J]. Plant Mol Biol Rep, 15(1): 8-15.

    ROZAS J, FERRER-MATA A, SANCHEZ-DELBARRIO JC, et al., 2017. DnaSP 6: DNA sequence polymorphism analysis of large data sets [J]. Mol Biol Evol, 34(12): 3299-3302.

    RUNGIS D, LAIVINS M, GAILITE A, et al., 2017. Genetic analysis of Latvian Salix alba L. and hybrid populations using nuclear and chloroplast DNA markers [J]. iForest, 10(2): 422-429.

    SALTONSTALL K, 2002. Cryptic invasion by a non-native genotype of the common reed, Phragmites australis, into North America [J]. P Natl Acad Sci USA, 99(4): 2445-2449.

    SCARCELLI N, BARNAUD A, EISERHARDT W, et al., 2011. A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in monocotyledons [J]. PLoS ONE, 6(5): e19954.

    SHAW J, LICKEY EB, BECK JT, et al., 2005. The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis [J]. Am J Bot, 92(1): 142.

    SHAW J, SHAFER HL, LEONARD OR, et al., 2014. Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: The tortoise and the hare IV [J]. Am J Bot, 101(11): 1987-2004.

    SINGH VK, MANGALAM AK, DWIVEDI S, et al., 1998. Primer premier: Program for design of degenerate primers from a protein sequence [J]. Biotechniques, 24(2): 318-319.

    THOMPSON JN, BURDON JJ, 1992. Gene-for-gene coevolution between plants and parasites [J]. Nature, 360(6400): 121-125.

    WANG F, CHEN WY, CAI RZM, et al., 2010. Gray system theory based analysis of hard branch cuttage of Salix oritrepha in alpine area [J]. Pratac Sci, (9): 86-90.? [王芳, 陳文業(yè), 才讓卓瑪, 等, 2010. 高寒區(qū)山生柳硬枝扦插灰色系統(tǒng)理論分析研究 [J]. 草業(yè)科學(xué), (9):86-90.]

    WANG W, WANG S, 2015. Analysis of plant community diversity of Salix oritrepha shrub meadow in Maqu County [J]. For Gansu, (6): 28-30.? [王瑋, 王盛, 2015. 瑪曲縣山生柳灌叢草甸植物群落多樣性分析 [J]. 甘肅林業(yè), (6):28-30.]

    WRIGHT S, 1950. Genetic structure of populations [J]. Brit Med J, 2(4669): 36.

    WRIGHT S, 1970. A Second Installment. (Book Reviews: Evolution and the genetics of populations. The theory of gene frequencies) [J]. J Biosoc Sci, 168(2): 722-723.

    WRIGHT S, 1982. The shifting balance theory and macroevolution [J]. Ann Rev Genet, 16(1): 1-20.

    WU ZQ, 2016. The whole chloroplast genome of shrub willows (Salix suchowensis) [J]. Mitochondrial DNA A, 27(3): 1-2.

    YANG LC, ZHONG ZB, LIU HC, et al., 2013. Community characteristics of alpine shrub Salix oritrepha [J]. Bull Bot Res, 33(6): 653-658.? [楊路存, 鐘澤兵, 劉何春, 等, 2013. 高寒山生柳灌叢的群落特征研究 [J]. 植物研究, 33(6):653-658.]

    YANG YH, ZHANG DX, LI YM, et al., 2004. Mitochondrial DNA diversity and preliminary biogeographic inference of the evolutionary history of the black-spotted pond frog Rana nigromaculata populations in China [J]. Acta Zool Sin, 50(2): 193-201.? [楊玉慧, 張德興, 李義明, 等, 2004. 中國黑斑蛙種群的線粒體DNA多樣性和生物地理演化過程的初探 [J]. 動(dòng)物學(xué)報(bào), 50(2):193-201.]

    YANG YT, HOU XY, WEI ZW, et al., 2018. Screening and genetic diversity analysis of chloroplast non-coding region in Leymus chinensis [J]. Acta Pratac Sin, 27(10): 147-157.? [楊艷婷, 侯向陽, 魏臻武, 等, 2018. 羊草葉綠體非編碼區(qū)多態(tài)性標(biāo)記篩選及群體遺傳多樣性 [J]. 草業(yè)學(xué)報(bào), 27(10):147-157.]

    ZHAO JK, XU LA, JIE HF, et al., 2008. RAPD analysis of population genetic diversity of Tamarix chinensis in Yellow River delta [J]. J Nanjing For Univ (Nat Sci Ed), 32(5): 56-60.? [趙景奎, 徐立安, 解荷峰, 等, 2008. 黃河三角洲檉柳群體遺傳多樣性RAPD分析 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 32(5):56-60.]

    猜你喜歡
    遺傳多樣性青藏高原
    破解青藏高原氣候變化的密碼
    軍事文摘(2022年20期)2023-01-10 07:18:22
    青藏高原上的“含羞花”
    給青藏高原的班公湖量體溫
    軍事文摘(2021年20期)2021-11-10 01:58:40
    為了讓青藏高原的天更藍(lán)、水更綠、草原更美
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    山西大豆自然群體遺傳多樣性的研究
    在线观看免费视频网站a站| 波多野结衣一区麻豆| 久热这里只有精品99| 一级黄色大片毛片| 久久国产亚洲av麻豆专区| 国产一区亚洲一区在线观看| 欧美成狂野欧美在线观看| 免费观看人在逋| kizo精华| avwww免费| 91成人精品电影| 午夜视频精品福利| 黄网站色视频无遮挡免费观看| 高清欧美精品videossex| 一区二区三区四区激情视频| 成人国产av品久久久| 视频区图区小说| 这个男人来自地球电影免费观看| 狠狠婷婷综合久久久久久88av| 国产真人三级小视频在线观看| 精品人妻1区二区| 黄色 视频免费看| 黄色视频在线播放观看不卡| 午夜日韩欧美国产| 亚洲中文日韩欧美视频| 黄色怎么调成土黄色| 2018国产大陆天天弄谢| 成人国产av品久久久| 啦啦啦中文免费视频观看日本| 欧美中文综合在线视频| 国产一区二区在线观看av| 日本午夜av视频| 国产成人影院久久av| 免费观看av网站的网址| 日韩熟女老妇一区二区性免费视频| 久久久欧美国产精品| 国产无遮挡羞羞视频在线观看| 国产三级黄色录像| 国产一卡二卡三卡精品| 黄片小视频在线播放| 中文字幕人妻丝袜制服| 国产亚洲av片在线观看秒播厂| 久久久久久免费高清国产稀缺| 老司机在亚洲福利影院| 在线 av 中文字幕| 欧美日韩av久久| 久久人人爽人人片av| 一级a爱视频在线免费观看| www日本在线高清视频| 夜夜骑夜夜射夜夜干| 精品卡一卡二卡四卡免费| 黄色毛片三级朝国网站| 人人妻人人添人人爽欧美一区卜| 两个人看的免费小视频| 久久久久精品国产欧美久久久 | 久久人妻福利社区极品人妻图片 | 校园人妻丝袜中文字幕| xxxhd国产人妻xxx| www.熟女人妻精品国产| 国产精品一区二区免费欧美 | 精品国产乱码久久久久久小说| 久久精品成人免费网站| 老鸭窝网址在线观看| a级毛片在线看网站| 欧美精品一区二区大全| 永久免费av网站大全| 又大又爽又粗| 国产精品偷伦视频观看了| 中国美女看黄片| 美女大奶头黄色视频| 狠狠婷婷综合久久久久久88av| 久久人妻熟女aⅴ| 国产成人精品久久久久久| 最近最新中文字幕大全免费视频 | 国产亚洲av片在线观看秒播厂| 啦啦啦啦在线视频资源| 午夜激情av网站| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品国产一区二区精华液| 欧美中文综合在线视频| 日韩av免费高清视频| 成人三级做爰电影| 国产男人的电影天堂91| 成年女人毛片免费观看观看9 | 欧美xxⅹ黑人| 免费高清在线观看日韩| 大陆偷拍与自拍| 精品人妻在线不人妻| 天堂俺去俺来也www色官网| 亚洲国产精品成人久久小说| 免费久久久久久久精品成人欧美视频| 美女主播在线视频| 国产一区二区三区av在线| 超色免费av| 免费看十八禁软件| 嫁个100分男人电影在线观看 | 五月天丁香电影| 中文乱码字字幕精品一区二区三区| 老司机在亚洲福利影院| 国产片特级美女逼逼视频| 国产黄色视频一区二区在线观看| 欧美成人午夜精品| 国产精品99久久99久久久不卡| 成年人午夜在线观看视频| 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区国产| 免费日韩欧美在线观看| 欧美国产精品va在线观看不卡| a级毛片黄视频| 国产黄色视频一区二区在线观看| 久久国产精品影院| 欧美日韩精品网址| 亚洲欧美精品综合一区二区三区| 777久久人妻少妇嫩草av网站| 超碰成人久久| 亚洲欧美中文字幕日韩二区| 9热在线视频观看99| 亚洲欧美日韩另类电影网站| 国产视频一区二区在线看| 久久久久国产精品人妻一区二区| 国产精品熟女久久久久浪| 国产精品久久久久久精品电影小说| 亚洲欧美一区二区三区国产| 秋霞在线观看毛片| 国产精品一二三区在线看| 三上悠亚av全集在线观看| av天堂在线播放| 国产午夜精品一二区理论片| 亚洲中文av在线| 人人妻人人澡人人爽人人夜夜| 欧美激情高清一区二区三区| 男女边摸边吃奶| 九草在线视频观看| 自线自在国产av| 妹子高潮喷水视频| 蜜桃国产av成人99| 狠狠精品人妻久久久久久综合| 中文字幕制服av| 好男人电影高清在线观看| 亚洲欧美激情在线| 黄网站色视频无遮挡免费观看| 丰满饥渴人妻一区二区三| 999久久久国产精品视频| 晚上一个人看的免费电影| 老汉色∧v一级毛片| 超碰97精品在线观看| 亚洲精品国产区一区二| 丝袜喷水一区| 国产野战对白在线观看| 999精品在线视频| 国产欧美日韩一区二区三区在线| 91字幕亚洲| tube8黄色片| 欧美97在线视频| 午夜福利,免费看| www.av在线官网国产| 天天添夜夜摸| 51午夜福利影视在线观看| 精品福利永久在线观看| a级毛片在线看网站| 免费人妻精品一区二区三区视频| 亚洲美女黄色视频免费看| 99热网站在线观看| kizo精华| 久久狼人影院| 午夜福利乱码中文字幕| 少妇被粗大的猛进出69影院| 国产三级黄色录像| 久久人妻福利社区极品人妻图片 | 国产人伦9x9x在线观看| 日本av手机在线免费观看| 大片免费播放器 马上看| 丝袜脚勾引网站| 日韩一区二区三区影片| 亚洲精品成人av观看孕妇| 亚洲欧美一区二区三区国产| 亚洲五月色婷婷综合| 水蜜桃什么品种好| 99热全是精品| 丰满饥渴人妻一区二区三| 久久人人97超碰香蕉20202| 午夜免费成人在线视频| 男女午夜视频在线观看| 国产在线一区二区三区精| 国产日韩欧美亚洲二区| av在线老鸭窝| 国产又色又爽无遮挡免| 人人妻人人爽人人添夜夜欢视频| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩精品亚洲av| 岛国毛片在线播放| 国产精品一区二区在线观看99| 日本a在线网址| 1024视频免费在线观看| 女警被强在线播放| 最新在线观看一区二区三区 | 久久人人爽av亚洲精品天堂| 精品少妇久久久久久888优播| 久久久久久久精品精品| 午夜日韩欧美国产| a 毛片基地| 婷婷成人精品国产| 国产成人系列免费观看| 黄片小视频在线播放| 天天添夜夜摸| 超色免费av| e午夜精品久久久久久久| 国产片内射在线| 日日摸夜夜添夜夜爱| 女人爽到高潮嗷嗷叫在线视频| 久久久久国产精品人妻一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产野战对白在线观看| 另类亚洲欧美激情| 成人亚洲欧美一区二区av| 亚洲精品自拍成人| 欧美日韩一级在线毛片| xxx大片免费视频| 亚洲精品av麻豆狂野| 男人爽女人下面视频在线观看| 97人妻天天添夜夜摸| www.自偷自拍.com| 热99久久久久精品小说推荐| videosex国产| 久久人妻熟女aⅴ| 亚洲色图 男人天堂 中文字幕| 操美女的视频在线观看| 亚洲情色 制服丝袜| 又黄又粗又硬又大视频| 久久精品久久久久久久性| 欧美+亚洲+日韩+国产| 国产亚洲精品第一综合不卡| 久久久精品免费免费高清| 国产成人精品在线电影| 久久精品国产亚洲av高清一级| 老鸭窝网址在线观看| 男人操女人黄网站| 亚洲国产欧美在线一区| av天堂久久9| 亚洲欧美清纯卡通| 国产精品久久久久久精品古装| av在线app专区| 日韩欧美一区视频在线观看| 99re6热这里在线精品视频| 国产亚洲精品久久久久5区| 久久性视频一级片| 七月丁香在线播放| 50天的宝宝边吃奶边哭怎么回事| 一级毛片我不卡| 亚洲欧美中文字幕日韩二区| 亚洲av综合色区一区| 十八禁高潮呻吟视频| 亚洲成人国产一区在线观看 | 精品一区二区三区av网在线观看 | 999精品在线视频| videos熟女内射| 亚洲,欧美,日韩| 欧美精品一区二区大全| 爱豆传媒免费全集在线观看| 久久99热这里只频精品6学生| 久热爱精品视频在线9| 亚洲图色成人| 午夜福利,免费看| 国产成人精品无人区| 最新在线观看一区二区三区 | 亚洲中文字幕日韩| 在线观看一区二区三区激情| 国产一区二区三区综合在线观看| 国产免费福利视频在线观看| 丰满饥渴人妻一区二区三| 欧美亚洲日本最大视频资源| 美女中出高潮动态图| 九草在线视频观看| 超色免费av| 欧美中文综合在线视频| 热99国产精品久久久久久7| 亚洲图色成人| 五月开心婷婷网| 青草久久国产| 久久久国产一区二区| 婷婷色综合www| 国产无遮挡羞羞视频在线观看| 免费少妇av软件| 女人久久www免费人成看片| 中国美女看黄片| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜爱| 免费一级毛片在线播放高清视频 | 久久99热这里只频精品6学生| 亚洲欧美精品综合一区二区三区| 欧美av亚洲av综合av国产av| 十八禁网站网址无遮挡| 亚洲欧美精品自产自拍| 精品一区二区三区四区五区乱码 | 黄色毛片三级朝国网站| 2021少妇久久久久久久久久久| 18禁观看日本| 午夜免费观看性视频| 久久久国产精品麻豆| cao死你这个sao货| 首页视频小说图片口味搜索 | 97在线人人人人妻| 欧美黑人欧美精品刺激| 欧美在线一区亚洲| 日韩制服骚丝袜av| 欧美久久黑人一区二区| 欧美性长视频在线观看| 中文字幕亚洲精品专区| 精品一品国产午夜福利视频| 捣出白浆h1v1| 国产精品 国内视频| 久久久精品国产亚洲av高清涩受| 国产人伦9x9x在线观看| 欧美日韩国产mv在线观看视频| 国产日韩欧美视频二区| 国产高清视频在线播放一区 | 欧美 日韩 精品 国产| 精品人妻在线不人妻| 51午夜福利影视在线观看| 黄色片一级片一级黄色片| 男女高潮啪啪啪动态图| 午夜福利影视在线免费观看| 妹子高潮喷水视频| 亚洲欧美色中文字幕在线| 国产一区二区三区综合在线观看| 老鸭窝网址在线观看| 视频区图区小说| 国精品久久久久久国模美| 久久人人爽人人片av| 午夜福利乱码中文字幕| 你懂的网址亚洲精品在线观看| 成人国语在线视频| 超碰成人久久| 亚洲国产看品久久| 亚洲一区二区三区欧美精品| 久久久久精品国产欧美久久久 | 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 欧美 日韩 精品 国产| 99国产精品免费福利视频| 超色免费av| 国产人伦9x9x在线观看| a 毛片基地| 亚洲国产毛片av蜜桃av| 老司机影院毛片| 亚洲免费av在线视频| 高清欧美精品videossex| 天天影视国产精品| 国产熟女欧美一区二区| 女警被强在线播放| netflix在线观看网站| 精品亚洲成国产av| 国产有黄有色有爽视频| 精品久久蜜臀av无| 国产精品一区二区精品视频观看| 国产一区有黄有色的免费视频| 午夜福利免费观看在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av片天天在线观看| 亚洲精品中文字幕在线视频| 人妻一区二区av| 国产精品偷伦视频观看了| 亚洲欧美日韩高清在线视频 | 国产黄频视频在线观看| 91精品国产国语对白视频| 久久久国产欧美日韩av| 国产亚洲精品久久久久5区| 两人在一起打扑克的视频| 日日摸夜夜添夜夜爱| 日韩制服丝袜自拍偷拍| 丝袜美腿诱惑在线| 婷婷丁香在线五月| 天天躁日日躁夜夜躁夜夜| 国产成人系列免费观看| 天堂中文最新版在线下载| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 考比视频在线观看| 亚洲一区二区三区欧美精品| 欧美大码av| 亚洲欧美激情在线| 91精品三级在线观看| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看 | 国产亚洲欧美在线一区二区| 99国产综合亚洲精品| 日本一区二区免费在线视频| 午夜91福利影院| 欧美 亚洲 国产 日韩一| 精品国产超薄肉色丝袜足j| 亚洲精品一二三| 午夜老司机福利片| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠躁躁| 蜜桃在线观看..| 最黄视频免费看| 国产精品 欧美亚洲| 色婷婷久久久亚洲欧美| 另类亚洲欧美激情| 亚洲黑人精品在线| 大香蕉久久成人网| 免费久久久久久久精品成人欧美视频| 超碰成人久久| 久热这里只有精品99| 男女之事视频高清在线观看 | 又大又爽又粗| 亚洲欧洲国产日韩| 制服人妻中文乱码| 国产免费一区二区三区四区乱码| 国产精品一区二区精品视频观看| 亚洲第一青青草原| 一级,二级,三级黄色视频| 日本vs欧美在线观看视频| 亚洲情色 制服丝袜| 狠狠精品人妻久久久久久综合| 中文字幕人妻丝袜一区二区| 欧美在线一区亚洲| 18禁裸乳无遮挡动漫免费视频| 又粗又硬又长又爽又黄的视频| 99久久99久久久精品蜜桃| 亚洲av综合色区一区| 欧美黄色片欧美黄色片| 亚洲激情五月婷婷啪啪| 91精品伊人久久大香线蕉| 国产无遮挡羞羞视频在线观看| kizo精华| 天天操日日干夜夜撸| 色视频在线一区二区三区| 久久久久国产精品人妻一区二区| 另类亚洲欧美激情| 在线观看免费午夜福利视频| 日日夜夜操网爽| 999精品在线视频| 免费女性裸体啪啪无遮挡网站| 午夜激情久久久久久久| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 只有这里有精品99| 黄片播放在线免费| 精品亚洲成国产av| 一级毛片黄色毛片免费观看视频| 欧美黄色淫秽网站| 亚洲三区欧美一区| 一二三四社区在线视频社区8| 亚洲欧美一区二区三区国产| 五月天丁香电影| 亚洲三区欧美一区| 日本黄色日本黄色录像| 婷婷色麻豆天堂久久| 国产人伦9x9x在线观看| 日韩一区二区三区影片| 99精国产麻豆久久婷婷| 99热网站在线观看| 人妻人人澡人人爽人人| 精品亚洲成国产av| 欧美激情高清一区二区三区| 亚洲免费av在线视频| 国产精品一二三区在线看| 男女之事视频高清在线观看 | 久久精品久久久久久噜噜老黄| 免费日韩欧美在线观看| 国产野战对白在线观看| 亚洲国产精品一区二区三区在线| 日本av手机在线免费观看| 成年av动漫网址| 久久亚洲精品不卡| 在线看a的网站| 一区二区三区四区激情视频| 成年av动漫网址| 久久中文字幕一级| 欧美日韩视频高清一区二区三区二| 午夜两性在线视频| 国产午夜精品一二区理论片| 美女脱内裤让男人舔精品视频| 麻豆av在线久日| 女警被强在线播放| 啦啦啦视频在线资源免费观看| 亚洲国产欧美日韩在线播放| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| e午夜精品久久久久久久| 热99国产精品久久久久久7| 国产欧美日韩精品亚洲av| 欧美人与性动交α欧美精品济南到| 国产1区2区3区精品| 伊人久久大香线蕉亚洲五| 日韩,欧美,国产一区二区三区| 日本午夜av视频| 亚洲 欧美一区二区三区| 午夜福利乱码中文字幕| 视频在线观看一区二区三区| 美女扒开内裤让男人捅视频| 9热在线视频观看99| 国产伦理片在线播放av一区| 国产成人影院久久av| 激情五月婷婷亚洲| 亚洲欧美中文字幕日韩二区| 男女无遮挡免费网站观看| 伊人亚洲综合成人网| 丰满迷人的少妇在线观看| 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 国产黄色免费在线视频| 一本综合久久免费| av欧美777| 最新的欧美精品一区二区| 国产人伦9x9x在线观看| 在线观看www视频免费| 夫妻性生交免费视频一级片| 精品国产乱码久久久久久小说| 欧美日韩国产mv在线观看视频| 丝袜在线中文字幕| 日韩熟女老妇一区二区性免费视频| 男女下面插进去视频免费观看| 亚洲综合色网址| 国产无遮挡羞羞视频在线观看| 成人免费观看视频高清| 亚洲伊人久久精品综合| 丰满迷人的少妇在线观看| 亚洲国产毛片av蜜桃av| 亚洲av日韩在线播放| 少妇人妻 视频| 天天躁日日躁夜夜躁夜夜| 91成人精品电影| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合色网址| 精品国产超薄肉色丝袜足j| 尾随美女入室| 又粗又硬又长又爽又黄的视频| 女人爽到高潮嗷嗷叫在线视频| 国产主播在线观看一区二区 | 日韩大片免费观看网站| 色播在线永久视频| 国产亚洲av片在线观看秒播厂| 国产有黄有色有爽视频| 久久久国产一区二区| 国产老妇伦熟女老妇高清| 国产一卡二卡三卡精品| 男人操女人黄网站| netflix在线观看网站| 午夜福利乱码中文字幕| 人成视频在线观看免费观看| 精品少妇久久久久久888优播| 久久国产精品男人的天堂亚洲| 欧美精品人与动牲交sv欧美| 欧美大码av| 黄色 视频免费看| 老司机影院毛片| 首页视频小说图片口味搜索 | 欧美国产精品va在线观看不卡| 秋霞在线观看毛片| 日本黄色日本黄色录像| 精品国产一区二区久久| 欧美精品啪啪一区二区三区 | 女性生殖器流出的白浆| 久久这里只有精品19| 国产老妇伦熟女老妇高清| 亚洲av电影在线进入| 国产男女内射视频| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 2018国产大陆天天弄谢| 黄网站色视频无遮挡免费观看| 大陆偷拍与自拍| 亚洲国产中文字幕在线视频| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠久久av| 欧美国产精品一级二级三级| 无遮挡黄片免费观看| 这个男人来自地球电影免费观看| 国产免费视频播放在线视频| 美女午夜性视频免费| 亚洲激情五月婷婷啪啪| 伦理电影免费视频| 久久狼人影院| 免费人妻精品一区二区三区视频| 丁香六月欧美| 考比视频在线观看| 欧美日韩成人在线一区二区| 久久99热这里只频精品6学生| 深夜精品福利| 国产精品亚洲av一区麻豆| 黄网站色视频无遮挡免费观看| 中文字幕人妻丝袜制服| 99九九在线精品视频| 午夜免费男女啪啪视频观看| 最黄视频免费看| 欧美精品啪啪一区二区三区 | 精品国产乱码久久久久久男人| 黑人猛操日本美女一级片| 丁香六月天网| av福利片在线| 亚洲自偷自拍图片 自拍| 亚洲,欧美精品.| 91麻豆av在线| 91九色精品人成在线观看| 黄色一级大片看看| 一级片'在线观看视频| 老司机靠b影院| 宅男免费午夜| 免费在线观看黄色视频的| 国产1区2区3区精品| 在线观看免费视频网站a站| 久久免费观看电影| 欧美日韩av久久| 亚洲av日韩精品久久久久久密 | 美女国产高潮福利片在线看| 男女床上黄色一级片免费看| 91老司机精品| 涩涩av久久男人的天堂| 汤姆久久久久久久影院中文字幕| 美国免费a级毛片| a级片在线免费高清观看视频| 亚洲国产av新网站| 久久久精品区二区三区|