• <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)展
    山西大豆自然群體遺傳多樣性的研究
    人体艺术视频欧美日本| 国产美女午夜福利| 乱系列少妇在线播放| 久久精品国产亚洲av天美| 中文字幕久久专区| 特大巨黑吊av在线直播| 联通29元200g的流量卡| 老司机影院成人| av卡一久久| 各种免费的搞黄视频| 男女免费视频国产| 插逼视频在线观看| 少妇高潮的动态图| 人人妻人人看人人澡| 国产大屁股一区二区在线视频| 国产一区有黄有色的免费视频| 成人黄色视频免费在线看| 免费人妻精品一区二区三区视频| 日韩av不卡免费在线播放| 亚洲精品aⅴ在线观看| 三级国产精品欧美在线观看| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 亚洲国产精品一区三区| 人妻系列 视频| 精品一品国产午夜福利视频| 国产高潮美女av| 久久久久视频综合| 亚洲国产精品成人久久小说| 秋霞伦理黄片| 亚洲第一区二区三区不卡| 色婷婷av一区二区三区视频| 美女cb高潮喷水在线观看| 亚洲图色成人| 人体艺术视频欧美日本| 国产亚洲欧美精品永久| 国产 精品1| 国产色婷婷99| 久久热精品热| 我要看日韩黄色一级片| 热re99久久精品国产66热6| 久久久欧美国产精品| av网站免费在线观看视频| 五月玫瑰六月丁香| 国产精品一区www在线观看| 国产成人a区在线观看| 亚洲第一区二区三区不卡| 97超视频在线观看视频| 在线免费十八禁| 又黄又爽又刺激的免费视频.| 老司机影院毛片| 亚洲欧洲国产日韩| 插阴视频在线观看视频| tube8黄色片| 99视频精品全部免费 在线| 亚洲av.av天堂| 99久久综合免费| 亚洲精品中文字幕在线视频 | 国产无遮挡羞羞视频在线观看| 欧美高清成人免费视频www| 午夜福利在线观看免费完整高清在| 亚洲综合色惰| 丰满少妇做爰视频| 狂野欧美激情性bbbbbb| 一区二区三区四区激情视频| 精品久久久久久久末码| 在线观看国产h片| 国产精品爽爽va在线观看网站| 国产综合精华液| 精品久久久久久久末码| 麻豆成人午夜福利视频| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久久性| 国产爽快片一区二区三区| 精品亚洲成国产av| 夜夜爽夜夜爽视频| 男女无遮挡免费网站观看| 一个人看视频在线观看www免费| 看免费成人av毛片| 看非洲黑人一级黄片| 免费观看无遮挡的男女| 纯流量卡能插随身wifi吗| 日韩国内少妇激情av| 观看av在线不卡| 搡老乐熟女国产| 亚洲欧美日韩另类电影网站 | 欧美日韩视频精品一区| 亚洲高清免费不卡视频| 亚洲第一av免费看| 亚洲精品久久久久久婷婷小说| 18禁在线播放成人免费| 97在线人人人人妻| 日韩欧美 国产精品| 国产精品三级大全| 免费看av在线观看网站| 免费av不卡在线播放| 亚洲国产精品专区欧美| 国产精品一区www在线观看| 国产又色又爽无遮挡免| 一区在线观看完整版| 国产高清有码在线观看视频| 熟妇人妻不卡中文字幕| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 麻豆成人午夜福利视频| 不卡视频在线观看欧美| 久久99热这里只有精品18| 一区二区三区乱码不卡18| 精品久久久噜噜| 交换朋友夫妻互换小说| 国产精品福利在线免费观看| 亚洲精品色激情综合| 一本久久精品| 美女中出高潮动态图| 日韩一本色道免费dvd| 国产高清国产精品国产三级 | 亚洲伊人久久精品综合| 亚洲国产色片| 国产毛片在线视频| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| h视频一区二区三区| 男人舔奶头视频| 免费大片18禁| 国产日韩欧美亚洲二区| 久久精品人妻少妇| 亚洲成人手机| 99久久精品一区二区三区| 国产黄色免费在线视频| 一区二区三区免费毛片| 人妻 亚洲 视频| 国产黄频视频在线观看| 成人无遮挡网站| 国产免费又黄又爽又色| 日本与韩国留学比较| 亚洲色图综合在线观看| 欧美一级a爱片免费观看看| 亚洲精品自拍成人| 天天躁日日操中文字幕| 99久久精品一区二区三区| 欧美精品一区二区免费开放| 亚洲综合精品二区| 欧美老熟妇乱子伦牲交| 久久久久精品久久久久真实原创| 久久久久网色| 五月开心婷婷网| 99久久中文字幕三级久久日本| 下体分泌物呈黄色| 久久久久国产网址| 久久人人爽人人爽人人片va| 99热全是精品| 亚洲成色77777| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩另类电影网站 | 成人18禁高潮啪啪吃奶动态图 | 99久久综合免费| 欧美97在线视频| 汤姆久久久久久久影院中文字幕| 亚洲成人手机| 国产黄色视频一区二区在线观看| 国国产精品蜜臀av免费| 我要看日韩黄色一级片| 国产精品久久久久久精品电影小说 | 在线精品无人区一区二区三 | 亚洲国产精品专区欧美| 欧美性感艳星| 国产无遮挡羞羞视频在线观看| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在| 视频区图区小说| 国产av国产精品国产| 美女视频免费永久观看网站| 91精品一卡2卡3卡4卡| 男女边摸边吃奶| 国国产精品蜜臀av免费| 国产 精品1| 精品久久久噜噜| 永久网站在线| 国产精品秋霞免费鲁丝片| 91精品伊人久久大香线蕉| 午夜老司机福利剧场| 热99国产精品久久久久久7| 18禁动态无遮挡网站| 国产伦理片在线播放av一区| 边亲边吃奶的免费视频| xxx大片免费视频| 亚洲av免费高清在线观看| 久久精品久久久久久久性| 国产高清三级在线| 成人亚洲精品一区在线观看 | 91在线精品国自产拍蜜月| 在线免费观看不下载黄p国产| 免费黄频网站在线观看国产| 五月开心婷婷网| 91aial.com中文字幕在线观看| 国产一区亚洲一区在线观看| 美女内射精品一级片tv| 91精品国产九色| 赤兔流量卡办理| 欧美日韩综合久久久久久| 少妇人妻 视频| 卡戴珊不雅视频在线播放| 国产 精品1| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线| 国产午夜精品久久久久久一区二区三区| 久久久久视频综合| 久久青草综合色| 午夜激情福利司机影院| 亚洲无线观看免费| 日韩一本色道免费dvd| av又黄又爽大尺度在线免费看| 国产成人a∨麻豆精品| 久久精品国产亚洲av天美| 免费av不卡在线播放| 亚洲激情五月婷婷啪啪| 日韩中文字幕视频在线看片 | 十八禁网站网址无遮挡 | 久久国产乱子免费精品| 青青草视频在线视频观看| 插阴视频在线观看视频| 老熟女久久久| 国产欧美亚洲国产| 只有这里有精品99| 国产精品国产三级国产av玫瑰| 成人18禁高潮啪啪吃奶动态图 | 色网站视频免费| 国产精品av视频在线免费观看| 舔av片在线| 日韩 亚洲 欧美在线| 国产成人freesex在线| 国产有黄有色有爽视频| 少妇人妻一区二区三区视频| 日本av手机在线免费观看| 男的添女的下面高潮视频| 男女免费视频国产| 日韩制服骚丝袜av| 啦啦啦中文免费视频观看日本| 成年免费大片在线观看| 亚洲,欧美,日韩| 街头女战士在线观看网站| 久久久国产一区二区| 亚洲四区av| 国产成人精品一,二区| 男女免费视频国产| 国产成人精品婷婷| 高清欧美精品videossex| 亚洲不卡免费看| 国产伦理片在线播放av一区| 日产精品乱码卡一卡2卡三| 舔av片在线| 久热久热在线精品观看| 久久久久久久久久成人| 只有这里有精品99| 少妇裸体淫交视频免费看高清| 婷婷色av中文字幕| 久久久久久久久久成人| 五月伊人婷婷丁香| 亚洲色图综合在线观看| 91精品国产九色| 三级经典国产精品| 菩萨蛮人人尽说江南好唐韦庄| 99久久综合免费| 极品少妇高潮喷水抽搐| 性高湖久久久久久久久免费观看| 3wmmmm亚洲av在线观看| av卡一久久| 日日啪夜夜爽| 欧美日韩综合久久久久久| xxx大片免费视频| 国产人妻一区二区三区在| 久久久精品免费免费高清| 久久久亚洲精品成人影院| av一本久久久久| 美女高潮的动态| 看免费成人av毛片| 亚洲av成人精品一二三区| 国产成人免费无遮挡视频| 99久久人妻综合| 免费观看av网站的网址| 777米奇影视久久| 狂野欧美白嫩少妇大欣赏| 久久久欧美国产精品| 亚洲精品aⅴ在线观看| 多毛熟女@视频| 麻豆精品久久久久久蜜桃| 99re6热这里在线精品视频| 日本色播在线视频| 久久人人爽av亚洲精品天堂 | 91久久精品电影网| 久久ye,这里只有精品| 天堂俺去俺来也www色官网| 成人一区二区视频在线观看| 国产爱豆传媒在线观看| 黑人猛操日本美女一级片| 国产精品99久久99久久久不卡 | 国产亚洲5aaaaa淫片| 午夜免费鲁丝| 女的被弄到高潮叫床怎么办| 国产熟女欧美一区二区| 精品午夜福利在线看| 中文乱码字字幕精品一区二区三区| 国产精品国产av在线观看| 国产黄色免费在线视频| 天美传媒精品一区二区| 一级av片app| 人妻一区二区av| 欧美日韩视频精品一区| 99九九线精品视频在线观看视频| 51国产日韩欧美| 18+在线观看网站| 五月伊人婷婷丁香| 六月丁香七月| 国产欧美日韩一区二区三区在线 | 久久久久网色| 中文字幕久久专区| 成年美女黄网站色视频大全免费 | 看免费成人av毛片| 五月伊人婷婷丁香| 国语对白做爰xxxⅹ性视频网站| 亚洲第一av免费看| 亚洲激情五月婷婷啪啪| av国产精品久久久久影院| 久久久久久人妻| 国产成人精品久久久久久| 国产亚洲欧美精品永久| 亚洲精品第二区| 色5月婷婷丁香| 男女下面进入的视频免费午夜| 一本久久精品| 青春草视频在线免费观看| 国产黄频视频在线观看| 91久久精品国产一区二区成人| 精品亚洲成a人片在线观看 | 亚洲国产色片| 国产精品不卡视频一区二区| 免费播放大片免费观看视频在线观看| 亚洲国产欧美人成| 妹子高潮喷水视频| 亚洲欧美日韩东京热| www.色视频.com| 久久国内精品自在自线图片| 国产又色又爽无遮挡免| 精品人妻偷拍中文字幕| 高清av免费在线| 国产 精品1| 亚洲国产日韩一区二区| 精品国产三级普通话版| 国产老妇伦熟女老妇高清| 久久毛片免费看一区二区三区| 亚洲精品乱码久久久久久按摩| 在线观看美女被高潮喷水网站| 内地一区二区视频在线| 国产精品久久久久久久电影| 亚洲内射少妇av| 久久久久久久精品精品| 精品人妻熟女av久视频| 日本欧美视频一区| 成人毛片60女人毛片免费| 丰满迷人的少妇在线观看| 特大巨黑吊av在线直播| 蜜桃亚洲精品一区二区三区| 日韩视频在线欧美| 看十八女毛片水多多多| 少妇人妻一区二区三区视频| 欧美日韩国产mv在线观看视频 | 亚洲精品日韩av片在线观看| 亚洲综合精品二区| h视频一区二区三区| 国产一区二区三区av在线| 欧美xxxx黑人xx丫x性爽| 中文字幕人妻熟人妻熟丝袜美| 欧美三级亚洲精品| 搡老乐熟女国产| 亚洲av在线观看美女高潮| 国产av精品麻豆| 亚洲内射少妇av| 国产爽快片一区二区三区| 色网站视频免费| 亚洲欧洲日产国产| 日韩一区二区三区影片| 蜜桃在线观看..| 亚洲av欧美aⅴ国产| 99九九线精品视频在线观看视频| 久热久热在线精品观看| 三级国产精品欧美在线观看| 人人妻人人添人人爽欧美一区卜 | 美女脱内裤让男人舔精品视频| 久久ye,这里只有精品| 亚洲精品乱码久久久v下载方式| 久久国产精品大桥未久av | 亚洲国产欧美在线一区| 日韩一区二区三区影片| 精品熟女少妇av免费看| 最近中文字幕高清免费大全6| 精品久久久噜噜| 免费高清在线观看视频在线观看| 美女内射精品一级片tv| 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区综合在线观看 | 小蜜桃在线观看免费完整版高清| 亚洲欧美一区二区三区黑人 | 美女高潮的动态| 日日啪夜夜爽| 十八禁网站网址无遮挡 | 欧美激情国产日韩精品一区| www.av在线官网国产| 人妻夜夜爽99麻豆av| 日韩伦理黄色片| 在线免费十八禁| 99国产精品免费福利视频| 日韩大片免费观看网站| 欧美 日韩 精品 国产| 男人舔奶头视频| 午夜福利视频精品| 国产av一区二区精品久久 | 亚洲精品国产成人久久av| 丝瓜视频免费看黄片| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| 深夜a级毛片| 最近的中文字幕免费完整| 国产黄片美女视频| 伦理电影免费视频| 亚洲第一av免费看| 精品国产乱码久久久久久小说| 国内精品宾馆在线| 国产高清三级在线| 久久99精品国语久久久| 丰满少妇做爰视频| 少妇精品久久久久久久| 亚洲精品久久久久久婷婷小说| 女性被躁到高潮视频| 欧美激情极品国产一区二区三区 | 亚洲av福利一区| 欧美成人a在线观看| 国产又色又爽无遮挡免| 成人国产麻豆网| 建设人人有责人人尽责人人享有的 | 免费观看的影片在线观看| 国产精品伦人一区二区| 国产探花极品一区二区| 高清av免费在线| 热re99久久精品国产66热6| 亚洲成人av在线免费| 亚洲精品国产色婷婷电影| 国产高清三级在线| 亚洲一区二区三区欧美精品| 能在线免费看毛片的网站| 亚洲,欧美,日韩| 五月玫瑰六月丁香| 精品人妻一区二区三区麻豆| 日韩av不卡免费在线播放| 在线观看免费高清a一片| 久久久精品免费免费高清| 日韩制服骚丝袜av| 一级av片app| 亚洲av不卡在线观看| videos熟女内射| 亚洲美女黄色视频免费看| 亚洲人与动物交配视频| 午夜福利在线观看免费完整高清在| 卡戴珊不雅视频在线播放| 久久精品熟女亚洲av麻豆精品| 久久久久久久久大av| 亚洲不卡免费看| 亚洲精品日韩av片在线观看| 日韩一区二区三区影片| 国产免费福利视频在线观看| 91精品伊人久久大香线蕉| 赤兔流量卡办理| 亚洲综合色惰| 亚洲精品乱码久久久v下载方式| 国产精品无大码| 97精品久久久久久久久久精品| 免费观看在线日韩| 久久 成人 亚洲| 十八禁网站网址无遮挡 | 久久影院123| 伊人久久国产一区二区| 黄片wwwwww| 啦啦啦啦在线视频资源| 边亲边吃奶的免费视频| 国产在视频线精品| 久久久久久久久久久免费av| 国产精品精品国产色婷婷| 日韩欧美一区视频在线观看 | 亚洲av福利一区| 国产高清三级在线| 美女国产视频在线观看| 晚上一个人看的免费电影| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 久久人妻熟女aⅴ| 插阴视频在线观看视频| 三级经典国产精品| 国精品久久久久久国模美| 欧美日本视频| 夫妻午夜视频| 男女无遮挡免费网站观看| 精品久久久久久电影网| 91精品国产九色| 亚洲美女黄色视频免费看| 国产美女午夜福利| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 欧美+日韩+精品| av一本久久久久| 中文字幕久久专区| 亚洲av福利一区| 国产爽快片一区二区三区| 少妇人妻一区二区三区视频| 99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 精品一区二区三区视频在线| 国产精品.久久久| 一二三四中文在线观看免费高清| 内射极品少妇av片p| 又粗又硬又长又爽又黄的视频| 亚洲国产毛片av蜜桃av| 一级毛片 在线播放| 国产久久久一区二区三区| 3wmmmm亚洲av在线观看| 性色av一级| 精品人妻偷拍中文字幕| 女性被躁到高潮视频| 一级毛片黄色毛片免费观看视频| 国语对白做爰xxxⅹ性视频网站| 五月开心婷婷网| 午夜福利视频精品| 国产白丝娇喘喷水9色精品| 男男h啪啪无遮挡| 夫妻性生交免费视频一级片| 少妇猛男粗大的猛烈进出视频| 91精品一卡2卡3卡4卡| 激情 狠狠 欧美| 亚洲精品一二三| 丰满迷人的少妇在线观看| 国产又色又爽无遮挡免| 深爱激情五月婷婷| 免费大片黄手机在线观看| 久久99热这里只频精品6学生| 日韩成人伦理影院| 欧美性感艳星| 在线精品无人区一区二区三 | 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| av专区在线播放| 一本色道久久久久久精品综合| 最近最新中文字幕大全电影3| 亚洲国产精品专区欧美| 99久久人妻综合| 亚洲国产精品专区欧美| 中文字幕人妻熟人妻熟丝袜美| 搡女人真爽免费视频火全软件| 久久久久久人妻| 日韩精品有码人妻一区| 免费黄网站久久成人精品| 国内揄拍国产精品人妻在线| 国产一级毛片在线| 黄片无遮挡物在线观看| 三级国产精品片| 黄片无遮挡物在线观看| 一级黄片播放器| 日韩一本色道免费dvd| 韩国av在线不卡| 男女边摸边吃奶| 性色av一级| 99热国产这里只有精品6| 亚洲不卡免费看| 97超视频在线观看视频| 联通29元200g的流量卡| 黑人猛操日本美女一级片| 在线播放无遮挡| 久久青草综合色| 老司机影院毛片| 最近手机中文字幕大全| 在线看a的网站| 国产黄色视频一区二区在线观看| 2022亚洲国产成人精品| 欧美3d第一页| 日韩精品有码人妻一区| 国产 一区 欧美 日韩| 国产精品人妻久久久影院| h视频一区二区三区| 日韩中文字幕视频在线看片 | 国产高潮美女av| 国产精品国产三级国产专区5o| 精华霜和精华液先用哪个| 国产乱人偷精品视频| 日韩中文字幕视频在线看片 | 夫妻午夜视频| 亚洲性久久影院| 男男h啪啪无遮挡| 熟女电影av网| 丰满少妇做爰视频| 我的老师免费观看完整版| 日韩 亚洲 欧美在线| 一区二区三区精品91| 狂野欧美白嫩少妇大欣赏| 久久国产乱子免费精品| 伊人久久国产一区二区| xxx大片免费视频| 日韩在线高清观看一区二区三区| 国精品久久久久久国模美| 亚洲第一区二区三区不卡| 亚洲欧美中文字幕日韩二区| 国产午夜精品久久久久久一区二区三区| 国产黄色视频一区二区在线观看| 亚洲熟女精品中文字幕| 九九爱精品视频在线观看| 欧美精品国产亚洲|