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

    唐古特虎耳草譜系地理學(xué)研究

    2018-04-08 03:09:53更吉卓瑪賈留坤夏銘澤高慶波
    西北植物學(xué)報 2018年2期
    關(guān)鍵詞:物種

    更吉卓瑪,李 彥,賈留坤,夏銘澤,高慶波

    (1 中國科學(xué)院西北高原生物研究所 中國科學(xué)院高山植物適應(yīng)與進(jìn)化重點(diǎn)實驗室,西寧 810001;2 中國科學(xué)院大學(xué),北京 100039;3 中國科學(xué)院西北高原生物研究所 青海省作物分子育種重點(diǎn)實驗室,西寧 810001)

    由于復(fù)雜的地理拓?fù)浣Y(jié)構(gòu)、微環(huán)境以及更廣泛的生態(tài)位,高山地區(qū)往往具有較高的生物多樣性,如高加索山、安第斯山、阿爾卑斯山等[1-3]。青藏高原及其周邊的高山地區(qū)包含了3個全球生物多樣性熱點(diǎn)地區(qū),即東部的橫斷山區(qū)、南部的喜馬拉雅山系及西部的阿爾泰和天山山脈[4]。一般認(rèn)為,與造山運(yùn)動相關(guān)的地理隔離和異域分化是導(dǎo)致青藏高原地區(qū)植物多樣性較高的原因之一[5-6]??焖佥椛浞只谇嗖馗咴貐^(qū)的很多植物類群中被檢測到,如馬先蒿屬(PedicularisL.)[7]、紅景天屬(RhodiolaL.)[8]、風(fēng)毛菊屬(SaussureaDC.)[9]、虎耳草屬(SaxifragaL.)[10-11]、槖吾屬-垂頭菊屬-蟹甲草屬復(fù)合群(LigulariaCass.-CremanthodiumBenth.-ParasenecioW. W. Sm. & J. Small)[12]等。然而,與較高分類階元的分化模式不同,青藏高原地區(qū)植物種內(nèi)居群間的分化可能更多的受到第四紀(jì)冰期氣候波動的影響。更新世晚期以來,青藏高原經(jīng)歷了至少4次規(guī)模較大的冰期,其中聶聶雄拉冰期(第二冰期,約720~500 ka)強(qiáng)度最大,冰川擴(kuò)張最廣,此后冰期強(qiáng)度逐漸減弱[13-14]。而全球末次盛冰期時青藏高原冰川的擴(kuò)張程度非常有限,這與歐洲和北美大陸末次冰期形成統(tǒng)一的冰蓋不同[15]。第四紀(jì)冰期與間冰期交替所引起的氣候波動對于現(xiàn)存生物類群的地理分布格局和遺傳結(jié)構(gòu)產(chǎn)生了極大的影響[16-17]。由于復(fù)雜的地理地貌特征、未被統(tǒng)一的大冰蓋覆蓋以及東亞季風(fēng)氣候和印度季風(fēng)氣候的影響,青藏高原地區(qū)的植物可能具有更復(fù)雜的第四紀(jì)冰期進(jìn)化歷史[18-19]。一方面,某些植物,如祁連圓柏(Juniperusprzewalskii)[20]、青海云杉(Piceacrassifolia)[21]等,在第四紀(jì)冰期時由高原臺面退縮到海拔較低的東南部邊緣避難所內(nèi),間冰期或冰期后再從邊緣避難所回遷到高原臺面,從而形成現(xiàn)有的地理分布格局。另一方面,由于第四紀(jì)冰期時青藏高原未被統(tǒng)一的大冰蓋覆蓋,僅有不同程度的山地冰川發(fā)育,高原臺面可能存在適宜的局部微環(huán)境而成為某些耐寒植物的避難所。對銀露梅(Potentillaglabra)[22]、露蕊烏頭(Aconitumgymnandrum)[23]和西川紅景天(Rhodiolaalsia)[24]等的研究揭示,在第四紀(jì)冰期時高原臺面的居群并沒有退縮到高原邊緣避難所,而是在高原臺面上的一個或幾個較大的避難所中保留下來,間冰期或冰期后再由避難所向附近擴(kuò)散。此外,對西藏圓柏復(fù)合群(Juniperustibeticacomplex)[25]、西藏沙棘(Hippophaetibetana)[26]、菊葉紅景天(Rhodiolachrysanthemifolia)[27]等的研究發(fā)現(xiàn),這些物種第四紀(jì)冰期時在高原臺面存在多個微型避難所,這些微型避難所在物種的分布范圍內(nèi)呈均勻分布。然而,具有不同分布范圍、生境和生活史特征的植物可能具有不同的第四紀(jì)冰期進(jìn)化歷史。對青藏高原地區(qū)更多的植物進(jìn)行譜系地理學(xué)研究,不但可為高山植物演化歷史與第四紀(jì)冰期的關(guān)系提供依據(jù),同時對青藏高原生物多樣性保護(hù)具有重要啟示。

    虎耳草屬是虎耳草科(Saxifragaceae)最大的屬,約有450~500種,主要分布在歐洲和亞洲的高山地區(qū),某些物種環(huán)北極分布,另有一些物種擴(kuò)散到南美洲的落基山脈、安第斯山脈和火地島[10-11]。中國產(chǎn)虎耳草屬植物約220種,主要分布在青藏高原-喜馬拉雅地區(qū)[28]?;⒍輰僦参锉粡V泛地用于譜系地理學(xué)和系統(tǒng)發(fā)育學(xué)研究,來揭示北極和高山地區(qū)植物的分化式樣[10-11, 29-32]。然而,青藏高原作為虎耳草屬物種多樣性中心之一,該地區(qū)虎耳草屬植物種內(nèi)分化與第四紀(jì)冰期的關(guān)系尚不明確。唐古特虎耳草(S.tanguticaEngl.)為多年生草本植物,生于海拔2 900~4 600 m的高山草甸、灌叢、石隙,分布于青海、甘肅、四川、西藏等地,是青藏高原地區(qū)高山草甸的重要組成成分,對維持高寒草甸生態(tài)系統(tǒng)的穩(wěn)定具有重要意義。本研究選取葉綠體DNA(cpDNA)非編碼區(qū)trnL-trnF和rpl16片段對青藏高原地區(qū)的唐古特虎耳草進(jìn)行譜系地理學(xué)研究,目的在于揭示該物種居群內(nèi)與居群間的遺傳分布式樣,討論唐古特虎耳草對于第四紀(jì)氣候波動的響應(yīng),探討其現(xiàn)有地理分布格局的歷史成因。

    1 材料和方法

    1.1 實驗材料

    研究材料于2007年至2014年采自青海、西藏、四川、甘肅,共采集唐古特虎耳草18個居群209個個體(表1,圖1)。每個居群隨機(jī)采集2~24個個體,個體之間相隔至少10 m以上。采集生長良好植株的新鮮嫩葉后立即用硅膠干燥,帶回實驗室存放于-20 ℃冰箱備用。憑證標(biāo)本由R. J. Gornall教授和高慶波副研究員共同鑒定,館藏于中國科學(xué)院西北高原植物研究所青藏高原生物標(biāo)本館(HNWP)。

    表1 唐古特虎耳草18個居群的采集信息

    餅狀圖表示每個居群中單倍型的頻率圖1 唐古特虎耳草18個居群的采樣圖以及葉綠體單倍型分布圖Pie charts show the proportions of haplotypes within each populationFig.1 Map of the 18 sampled populations of S. tangutica and the distribution of cpDNA haplotypes in the species

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

    本研究采用改良的CTAB法從硅膠干燥的葉片中提取總DNA[33]。采用通用引物對所有個體的cpDNAtrnL-trnF和rpl16片段進(jìn)行PCR擴(kuò)增[34-35]。PCR擴(kuò)增反應(yīng)體系為50 μL:包含5 μL的10 × PCR Buffer(含1.5 mmol/L MgCl2),2 μL的10 mmol/L dNTPs,正反引物各1.0 μL(5 pmol/L),TaqDNA聚合酶(TaKaRa,大連)0.4 μL,1 μL的DNA模板,用雙蒸水補(bǔ)足到50 μL。PCR擴(kuò)增反應(yīng)程序:94 ℃預(yù)變性10 min;94 ℃變性1 min,58 ℃退火50 s,72 ℃延伸1 min,32循環(huán);72 ℃延伸10 min。PCR產(chǎn)物用1%的瓊脂糖凝膠電泳進(jìn)行檢測,然后送測序公司進(jìn)行雙向測序。

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

    利用Chromas 2.6.2軟件(http://www.technelysium.com.au)對測序峰圖進(jìn)行肉眼核對。利用MEGA 7.0.26軟件對序列進(jìn)行對位排列并統(tǒng)計序列的堿基組成[36],發(fā)現(xiàn)trnL-trnF序列上存在一段約20 bp的poly A結(jié)構(gòu),為了確保堿基變異的可靠性,poly A及兩側(cè)約40 bp的堿基被刪除。同樣,對rpl16序列中出現(xiàn)的大于11 bp的poly A和poly T結(jié)構(gòu)也進(jìn)行了類似的刪除處理。利用DnaSP 5.10軟件統(tǒng)計cpDNA聯(lián)合序列的變異位點(diǎn)并確定單倍型[37]??傔z傳多樣性(HT)、居群內(nèi)平均遺傳多樣性(Hs)、居群間遺傳分化系數(shù)GST和NST值均利用PERMUT軟件進(jìn)行計算[38]。其中GST僅考慮單倍型的頻率,而NST兼顧了單倍型頻率和單倍型之間的相似性。使用U-統(tǒng)計方法對GST和NST進(jìn)行比較(1 000次重復(fù)的置換檢驗)以檢測單倍型變異的地理分布模式。NST大于GST且顯著,表明居群在分布范圍內(nèi)具有顯著的譜系地理結(jié)構(gòu)。由于PERMUT軟件要求每個居群至少要有3個個體,所以,居群P13(2個個體)在進(jìn)行PERMUT分析時未包含在內(nèi)。

    運(yùn)用Arlequin 3.5.2計算每個居群的單倍型多樣性(He)和核苷酸多樣性(π)[39]。利用Arlequin軟件包中的分子變異分析(AMOVA)檢測居群內(nèi)和居群間的遺傳變異水平,并對單倍型分布的FST進(jìn)行評價(1 000次置換檢驗)[40]。Tajima’sD、Fu & Li’sD*兩種無限突變位點(diǎn)模型的中性檢驗及歧點(diǎn)分布分析都在DnaSP 5.10程序中完成[41-42]。

    以最大簡約法為原則,利用NETWORK 4.6.1.3軟件構(gòu)建單倍型中央連接網(wǎng)狀圖[43-44]。在運(yùn)行該軟件時,對于多堿基的插入、缺失或突變位點(diǎn),被認(rèn)為是通過一步突變形成。在單倍型Network圖中,原始的單倍型多存在于網(wǎng)絡(luò)圖的內(nèi)部,而近代起源的單倍型發(fā)生在網(wǎng)絡(luò)圖的邊緣位置[45-46]。

    2 結(jié)果與分析

    2.1 唐古特虎耳草的單倍型分布

    對18個居群209個個體的trnL-trnF和rpl16片段進(jìn)行測序,MEGA軟件對位排列后的矩陣長度為1 626 bp,通過DnaSP軟件鑒定出74個單倍型,有25個簡約信息位點(diǎn)(表2)。在檢測出的74個單倍型中,有53個單倍型是居群特有單倍型,單倍型H5出現(xiàn)的頻率最高,有46個個體擁有此單倍型。表3顯示在18個居群中,16個居群擁有特有單倍型。居群的遺傳多樣性(He)范圍為0.490 9~1.000 0,平均為0.853 6;居群的核苷酸多樣性(π)范圍為0.007 506~0.046 481,平均為0.022 145。

    2.2 居群遺傳結(jié)構(gòu)特征

    通過PERMUT程序計算得出唐古特虎耳草居群內(nèi)平均遺傳多樣性HS值(0.843)、總遺傳多樣性HT值(0.933)、居群間遺傳分化GST(0.097)和NST值(0.109)。使用U-統(tǒng)計方法對唐古特虎耳草單倍型變異的地理結(jié)構(gòu)進(jìn)行檢驗后發(fā)現(xiàn)NST大于GST,但不顯著(P> 0.05),且居群間的遺傳分化水平很低,表明唐古特虎耳草在現(xiàn)有分布區(qū)域內(nèi)沒有明顯的譜系地理結(jié)構(gòu)。分子變異分析(AMOVA)的結(jié)果表明,居群間的遺傳變異僅為8.15%(FST= 0.081 45;P< 0.01),而居群內(nèi)的遺傳變異為91.85%(表4),進(jìn)一步揭示了唐古特虎耳草遺傳變異主要存在于居群內(nèi)。

    2.3 基于cpDNA聯(lián)合序列的居群歷史動態(tài)分析

    對唐古特虎耳草cpDNA聯(lián)合序列所有個體進(jìn)行中性檢驗,結(jié)果表明,Tajima’sD(-2.045 07,P< 0.05)和Fu & Li’D*(-3.629 27,P< 0.05)均為顯著的負(fù)值,推測其分布范圍或者個體數(shù)近期經(jīng)歷過擴(kuò)張。歧點(diǎn)分布分析結(jié)果顯示,在唐古特虎耳草整個分布范圍內(nèi)觀測到的歧點(diǎn)分布呈單峰曲線(圖2),其SSD值(0.007 2)和Raggednes index值(0.007 6)均不顯著(P> 0.05),說明該物種在整個分布范圍內(nèi)經(jīng)歷過近期擴(kuò)張現(xiàn)象,這與中性檢驗結(jié)果一致。

    2.4 單倍型之間的系統(tǒng)發(fā)育關(guān)系

    最大簡約法(maximum parsimony, MP)、最大似然法(maximum likelihood, ML)、貝葉斯法(Bayesian inference, BI)構(gòu)建的單倍型系統(tǒng)發(fā)育樹均呈現(xiàn)“刷子”結(jié)構(gòu),不能很好地揭示單倍型之間的系統(tǒng)發(fā)育關(guān)系,也表明了唐古特虎耳草單倍型的快速輻射分化(結(jié)果未顯示)。NETWORK軟件可更好揭示遺傳分化較淺序列之間的系統(tǒng)發(fā)育關(guān)系?;谧畲笏迫环?gòu)建74個單倍型的中央鏈接網(wǎng)狀圖,結(jié)果顯示(圖3),居群共享單倍型,如H1、H2、H3、H5、H16等位于網(wǎng)狀圖的中央位置,而居群特有單倍型一般位于網(wǎng)狀圖的邊緣。根據(jù)溯祖理論,共享單倍型為較古老的單倍型,而居群特有單倍型可能是通過近期輻射分化衍生而來的較年輕的單倍型。

    表2 唐古特虎耳草74個單倍型trnL-trnF和rpl16片段簡約信息位點(diǎn)

    續(xù)表2Continued Table 2

    單倍型Hap簡約信息位點(diǎn)ParsimonyinformativesitetrnL-trnFrpl16230355414456611792813928966988100010901124119213341352135913601414145114571465147915511562H53-------------------------H54--------------?A----------H55---------?C--?A------------H56--------?T--------?GC-----?TH57-----------?T-?A------?A---?TH58-------------------------H59-------------------------H60----------------?A--------H61----------?A--------------H62-------------------------H63---------?C--?A------------H64--------?T---------------?TH65-C-----C-?C---------------H66----C------?T-----------?GTH67--------?T--?T-------?AA-?A-?TH68-------------------------H69-C----A-------?A----------H70--------------?A----------H71-------------------------H72A----T------------------?TH73----------------?A--------H74----------?A--------------

    注:-代表與單倍型H1相同的堿基

    Note: - stands for the same base ashaplotypeH1

    表3 唐古特虎耳草18個居群的單倍型組成、遺傳多樣性和核苷酸多樣性

    表4 唐古特虎耳草18個居群的分子遺傳變異(AMOVA)分析結(jié)果

    圖2 唐古特虎耳草209個個體葉綠體DNA2個片段聯(lián)合數(shù)據(jù)的歧點(diǎn)分布分析Fig.2 Mismatch distribution for 209 individuals of S. tangutica based on two cpDNA fragments

    3 討 論

    3.1 唐古特虎耳草的居群遺傳結(jié)構(gòu)

    遺傳變異分析表明唐古特虎耳草cpDNA片段總遺傳多樣性水平較高,HT= 0.933。類似的結(jié)果也在青藏高原很多高山植物中被檢測到,如偏花報春(Primulasecundiflora,HT= 0.966)[47]、西藏沙棘(Hippophaetibetana,HT= 0.956)[26]、綿參(Eriophytonwallichii,HT= 0.979)[48]、西川紅景天(Rhodiolaalsia,HT= 0.950)[24]、菊葉紅景天(Rhodiolachrysanthemifolia,HT= 0.937)[27]等。然而,這些物種大多數(shù)都具有低的居群內(nèi)平均遺傳多樣性。在本研究中,唐古特虎耳草不僅具有較高的總遺傳多樣性,居群內(nèi)的平均遺傳多樣性也很高(HS= 0.843),表明唐古特虎耳草較高的遺傳變異主要由居群內(nèi)遺傳變異所致。這一結(jié)果也得到了分子變異分析(AMOVA)的證實。AMOVA分析表明,唐古特虎耳草居群內(nèi)遺傳變異占總變異的91.85%,而居群間遺傳變異僅為8.15%,進(jìn)一步揭示了唐古特虎耳草遺傳變異主要存在于居群內(nèi)。如此高的居群內(nèi)遺傳變異在青藏高原的植物中并不多見,推測造成該遺傳結(jié)構(gòu)的原因如下:一方面,由于檢測到的74個單倍型中只有少數(shù)共享單倍型,但其個體數(shù)所占比例很高(67.46%),并且這些共享單倍型在居群中隨機(jī)分布,因此降低了居群間的遺傳差異;另一方面,居群特有單倍型比例高(53/74),某些居群固定多個特有單倍型,從而增加了居群內(nèi)的遺傳變異。特有單倍型比例較高的物種往往不能揭示出明顯的譜系地理格局[22, 24, 27]。 此外,對虎耳草屬山羊臭組(Saxifragasect.CiliataeHaw.)系統(tǒng)發(fā)育和物種形成的研究表明,唐古特虎耳草所在的唐古拉亞組(subsect.HirculoideaeEngl. & Irmsch.)具有較快的分化速率,該亞組約110個物種是在最近的2.12百萬年(Ma)分化而來[10]。由此推測,唐古特虎耳草具有更短的居群分化歷史。事實上,基于NETWORK構(gòu)建的單倍型系統(tǒng)發(fā)育關(guān)系表明,大量的特有單倍型位于NETWORK圖的邊緣,揭示了唐古特虎耳草單倍型的近期快速分化。在第四紀(jì)冰期前,唐古特虎耳草可能廣泛分布于青藏高原及其周邊地區(qū),并廣泛固定某些古老單倍型(如H5)。隨著第四紀(jì)冰期氣候的波動,原來連續(xù)的居群片段化,形成了相互隔離的居群,并經(jīng)歷了異域分化,產(chǎn)生大量特有單倍型。在此過程中,瓶頸效應(yīng)和遺傳漂變等作用使得特有單倍型在居群中固定下來,抹去了居群原有的部分遺傳結(jié)構(gòu)。由于單獨(dú)進(jìn)化時間較短,特有單倍型之間的差異大多表現(xiàn)在單個堿基的變異。

    圓的大小表示單倍型頻率,實心方塊代表缺失單倍型圖3 唐古特虎耳草cpDNA聯(lián)合片段74種單倍型的中央連接網(wǎng)狀圖Circle size is proportional to haplotype frequencies, parallelograms represent missing haplotypesFig.3 Median-joining network of 74 haplotypes of S. tangutica based on two cpDNA fragments

    3.2 冰期避難所及居群進(jìn)化歷史的推測

    第四紀(jì)循環(huán)的氣候波動對生物區(qū)系的分布和遺傳結(jié)構(gòu)產(chǎn)生了巨大的影響[49]。植物對第四紀(jì)冰期的反應(yīng)可以根據(jù)現(xiàn)有的遺傳結(jié)構(gòu)進(jìn)行追溯。于海彬等[50]綜述了36種高山植物的譜系地理格局,總結(jié)出高山植物面對青藏高原復(fù)雜的地理結(jié)構(gòu)、微氣候和生境條件所呈現(xiàn)出的譜系歷史。一部分高山植物在第四紀(jì)冰期時退縮到青藏高原邊緣避難所,冰期后回遷到高原臺面。符合這種“退縮-回遷”模式的物種,其遺傳結(jié)構(gòu)表現(xiàn)為遺傳多樣性較高的居群和古老單倍型往往集中在高原邊緣避難所內(nèi),居群遺傳多樣性由高原邊緣向高原臺面遞減。此外,一些耐寒植物,在冰期時沒有完全退縮到低海拔的邊緣地區(qū),而是在高原臺面形成了一個或幾個相互隔離的避難所,冰期后發(fā)生局部擴(kuò)張。符合這種“高原臺面避難所-局部擴(kuò)張”模式的物種,其遺傳結(jié)構(gòu)表現(xiàn)為特有單倍型或者遺傳多樣性較高的居群分布在高原臺面上的一個或幾個相互隔離的區(qū)域。另外,某些物種第四紀(jì)冰期時在高原臺面形成多個微型避難所,間冰期或冰期后發(fā)生小范圍的擴(kuò)張,其遺傳多樣性較高的居群和特有單倍型在物種的分布范圍內(nèi)均勻分布。

    本研究沒有揭示出唐古特虎耳草遺傳多樣性較高的某一地理區(qū)域,也沒有發(fā)現(xiàn)遺傳多樣性由高原邊緣向高原臺面遞減。相反,本研究所涉及的唐古特虎耳草居群普遍具有較高的遺傳多樣性,特有單倍型在分布范圍內(nèi)均勻分布。我們推測,由于青藏高原地形的高度異質(zhì)性以及復(fù)雜的局部微氣候和生境條件,在第四紀(jì)冰期時,唐古特虎耳草在高原臺面上未被冰川覆蓋的多個微型避難所中存活下來。唐古特虎耳草居群對第四紀(jì)冰期氣候波動的反應(yīng)可能更多地表現(xiàn)為垂直海拔高度的遷移,而非大規(guī)模的水平遷移,這與銀露梅(Potentillaglabra)[22]和西川紅景天(Rhodiolachrysanthemifolia)[27]的進(jìn)化歷史相似。居群在垂直遷移的過程中發(fā)生了瓶頸效應(yīng)和奠基者效應(yīng),居群隨機(jī)固定少量古老單倍型。隨后的異域分化和遺傳漂變使得較年輕的衍生特有單倍型在居群中固定下來,形成了唐古特虎耳草現(xiàn)有的居群遺傳結(jié)構(gòu)。居群歷史動態(tài)檢驗結(jié)果表明,Tajima’sD和Fu & Li’D*均為顯著的負(fù)值,歧點(diǎn)分布分析呈單峰曲線,說明唐古特虎耳草經(jīng)歷過近期分布范圍或者個體數(shù)的擴(kuò)張[51-53]。物種經(jīng)歷分布范圍擴(kuò)張事件,其原分布中心應(yīng)具有較高的遺傳多樣性和特有單倍型,并且原分布中心的居群包含擴(kuò)散區(qū)域居群所具有的單倍型。但在本研究中,大量的特有單倍型在唐古特虎耳草的整個分布范圍內(nèi)均勻分布,遺傳多樣性和核苷酸多樣性較高的居群也呈現(xiàn)均勻分布,并未檢測到分布范圍擴(kuò)張的信號。由此推測,該物種可能經(jīng)歷過個體數(shù)快速擴(kuò)張事件,從而導(dǎo)致了中性檢驗的負(fù)值和錯配分布的單峰曲線。

    綜上所述,唐古特虎耳草具有較高的遺傳變異,且主要由居群內(nèi)遺傳變異所致,大量特有單倍型是近期輻射分化的結(jié)果。唐古特虎耳草的進(jìn)化歷史與西藏圓柏復(fù)合群(Juniperustibeticacomplex)[25]、西藏沙棘(Hippophaetibetana)[26]、菊葉紅景天(Rhodiolachrysanthemifolia)[27]相似,第四紀(jì)冰期時在高原臺面存在多個微型避難所。第四紀(jì)冰期氣候波動使得原來連續(xù)的居群片段化,形成了相互隔離的居群,并發(fā)生異域快速分化,形成大量特有單倍型。唐古特虎耳草居群對第四紀(jì)冰期氣候波動的反應(yīng)可能更多地表現(xiàn)為垂直海拔高度的遷移,而非大規(guī)模的水平遷移。此研究為青藏高原地區(qū)譜系地理學(xué)研究提供了一個新的案例,為了解第四紀(jì)冰期高山植物的避難所和現(xiàn)有分布格局的成因提供了參考依據(jù)。

    參考文獻(xiàn):

    [1]MYERS N, MITTERMEIER R A, MITTERMEIER C G,etal. Biodiversity hotspots for conservation priorities [J].Nature, 2000,403(6 772): 853-858.

    [2]HOORN C, MOSBRUGGER V, MULCH A,etal. Biodiversity from mountain building [J].NatureGeoscience, 2013,6(3): 154.

    [3]HUGHES C E. The tropical Andean plant diversity powerhouse [J].NewPhytologist, 2016,210(4): 1 152-1 154.

    [5]WEN J, ZHANG J-Q, NIE Z-L,etal. Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau [J].FrontiersinGenetics, 2014,5: 4.

    [6]XING Y, REE R H. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2017,114(17): E3444-E3451.

    [7]EATON D A R, FENSTER C B, HEREFORD J,etal. Floral diversity and community structure inPedicularis(Orobanchaceae) [J].Ecology, 2012,93(8): S182-S194.

    [8]ZHANG J Q, MENG S Y, WEN J,etal. Phylogenetic relationships and character evolution ofRhodiola(Crassulaceae) based on nuclear ribosomal ITS and plastidtrnL-FandpsbA-trnHsequences [J].SystematicBotany, 2014,39(2): 441-451.

    [9]WANG Y J, SUSANNA A, VON RAAB-STRAUBE E,etal. Island-like radiation ofSaussurea(Asteraceae: Cardueae) trigged by uplifts of the Qinghai-Tibetan Plateau [J].BotanicalJournaloftheLinneanSociety, 2009,97(4): 893-903.

    [10]GAO Q B, LI Y H, GORNALL R J,etal. Phylogeny and speciation inSaxifragasect.Ciliatae(Saxifragaceae): evidence frompsbA-trnH,trnL-Fand ITS sequences [J].Taxon, 2015,64(4): 703-713.

    [11]EBERSBACH J, MUELLNER-RIEHL A N, MICHALAK I,etal. In and out of the Qinghai-Tibet Plateau: divergence time estimation and historical biogeography of the large arctic-alpine genusSaxifragaL [J].JournalofBiogeography, 2017,44(4): 900-910.

    [12]LIU J Q, WANG Y J, WANG A L,etal. Radiation and diversification within theLigularia-Cremanthodium-Paraseneciocomplex (Asteraceae) trigged by uplift of the Qinghai-Tibetan Plateau [J].MolecularPhylogeneticsandEvolution, 2006,38(1): 31-49.

    [13]ZHENG B, XU Q, SHEN Y. The relationship between climate change and Quaternary glacial cycles on the Qinghai-Tibetan Plateau: review and speculation [J].QuaternaryInternational, 2002,97-98(1): 93-101.

    [14]OWEN L A, BENN D I. Equilibrium-line altitudes of the Last Glacial Maximum for the Himalaya and Tibet: an assessment and evaluation of results [J].QuaternaryInternational, 2005,138-139(3): 55-78.

    [15]LEHMKUHL F, OWEN L A. Late Quaternary glaciation of Tibet and the bordering mountains: a review [J].Boreas, 2005,34(2): 87-100.

    [16]HEWITT G M. Post-glacial re-colonization of European biota [J].BiologicalJournaloftheLinneanSociety, 1999,68(1/2): 87-112.

    [17]HEWITT G M. Speciation, hybrid zones and phylogeography-or seeing genes in space and time [J].MolecularEcology, 2001,10(3): 537-549.

    [18]ZHANG D, FENGQUAN L, JIANMIN B. Eco-environmental effects of the Qinghai-Tibet Plateau uplift during the Quaternary in China [J].EnvironmentalGeology, 2000,39(12): 1 352-1 358.

    [19]HEWITT G M. Genetic consequences of climatic oscillations in the Quaternary [J].PhilosophicalTransactionsoftheRoyalSocietyofLondonB:BiologicalSciences, 2004,359(1 442): 183-195.

    [20]ZHANG Q, CHIANG T Y, GEORGE M,etal. Phylogeography of the Qinghai-Tibetan Plateau endemicJuniperusprzewalskii(Cupressaceae) inferred from chloroplast DNA sequence variation [J].MolecularEcology, 2005,14(11): 3 513-3 524.

    [21]MENG L, YANG R, ABBOTT R J,etal. Mitochondrial and chloroplast phylogeography ofPiceacrassifoliaKom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands [J].MolecularEcology, 2007,16(19): 4 128-4 137.

    [22]WANG L Y, IKEDA H, LIU T L,etal. Repeated range expansion and glacial endurance ofPotentillaglabra(Rosaceae) in the Qinghai-Tibetan Plateau [J].JournalofIntegrativePlantBiology, 2009,51(7): 698-706.

    [23]WANG L, ABBOTT R J, ZHENG W,etal. History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau:Aconitumgymnandrum(Ranunculaceae) [J].MolecularEcology, 2009,18(4): 709-721.

    [24]GAO Q, ZHANG D, DUAN Y,etal. Intraspecific divergences ofRhodiolaalsia(Crassulaceae) based on plastid DNA and internal transcribed spacer fragments [J].BotanicalJournaloftheLinneanSociety, 2012,168(2): 204-215.

    [25]OPGENOORTH L, VENDRAMIN G G, MAO K,etal. Tree endurance on the Tibetan Plateau marks the world’s highest known tree line of the Last Glacial Maximum [J].NewPhytologist, 2010,185(1): 332-342.

    [26]WANG H, QIONG L, SUN K,etal. Phylogeographic structure ofHippophaetibetana(Elaeagnaceae) highlights the highest microrefugia and the rapid uplift of the Qinghai-Tibetan Plateau [J].MolecularEcology, 2010,19(14): 2 964-2 979.

    [27]GAO Q B, ZHANG F Q, XING R,etal. Phylogeographic study revealed microrefugia for an endemic species on the Qinghai-Tibetan Plateau:Rhodiolachrysanthemifolia(Crassulaceae) [J].PlantSystematicsandEvolution, 2016,302(9):1 179-1 193.

    [28]WU C Y, RAVEN P H. Flora of China [M]. Beijing: Science Press; St. Louis: Missouri Botanical Garden Press, 2001,8: 208-344.

    [29]ABBOTT R J, COMES H P. Evolution in the Arctic: a phylogeographic analysis of the circumarctic plant,Saxifragaoppositifolia(Purple saxifrage) [J].NewPhytologist, 2003,161(1): 211-224.

    [30]DECHAINE E G, ANDERSON S A, MCNEW J M,etal. On the evolutionary and biogeographic history ofSaxifragasect.Trachyphyllum(Gaud.) Koch (Saxifragaceae Juss.) [J].PLoSOne, 2013,8(7): e69814.

    [31]EBERSBACH J, SCHNITZLER J, FAVRE A,etal. Evolutionary radiations in the species-rich mountain genusSaxifragaL [J].BMCEvolutionaryBiology, 2017,17(1): 119.

    [32]GAO Q B, LI Y, GENGJI Z M,etal. Population genetic differentiation and taxonomy of three closely related species ofSaxifraga(Saxifragaceae) from southern Tibet and the Hengduan Mountains [J].FrontiersinPlantScience, 2017,8: 1 325.

    [33]DOYLE J J, DOYLE J L. A rapid DNA isolation procedure for small quantities of fresh leaf material [J].PhytochemicalBulletin,BotanicalSocietyofAmerica, 1987,19(1): 11-15.

    [34]TABERLET P, GIELLY L, PAUTOU G,etal. Universal primers for amplification of three non-coding regions of chloroplast [J].PlantMolecularBiology, 1991,17(5): 1 105-1 109.

    [35]SCARCELLI N, BARNAUD A, EISERHARDT W,etal. A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in Monocotyledons [J].PLoSOne, 2011,6(5): e19954.

    [36]KUMAR S, STECHER G, TAMURA K. MEGA 7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets [J].MolecularBiologyandEvolution, 2016,33(7): 1 870-1 874.

    [37]LIBRADO P, ROZAS J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [J].Bioinformatics, 2009,25(11): 1 451-1 452.

    [38]PONS O, PETIT R J. Measuring and testing genetic differentiation with ordered versus unordered alleles [J].Genetics, 1996,144(3): 1 237-1 245.

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

    [40]WEIR B S, COCKERHAM C C. Estimating F-Statistics for the analysis of population structure [J].Evolution, 1984,38(6): 1 358-1 370.

    [41]TAJIMA F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J].Genetics, 1989,123(3): 585-595.

    [42]FU Y X, LI W H. Statistical tests of neutrality of mutations [J].Genetics, 1993,133(3): 693-709.

    [43]BANDELT H J, FORSTER P, R?HL A. Median-joining networks for inferring intraspecific phylogenies [J].MolecularBiologyandEvolution, 1999,16(1): 37-48.

    [44]POLZIN T, DANESHMAND S V. On Steiner trees and minimum spanning trees in hypergrahs [J].OperationsResearchLetters, 2003,31(1): 12-20.

    [45]GOLDING G B. The detection of deleterious selection using ancestors inferred from a phylogenetic history [J].GeneticsResearch, 1987,49(1): 71-82.

    [46]CRANDALL K A, TEMPLETON A R. Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction [J].Genetics, 1993,134(3): 959-969.

    [47]WANG F Y, GONG X, HU C M,etal. Phylogeography of an alpine speciesPrimulasecundiflorainferred from the chloroplast DNA sequence variation [J].JournalofSystematicsandEvolution, 2008,46(1): 13-22.

    [48]王曉雄,樂霽培,孫航,等.青藏高原高山流石灘特有植物綿參的譜系地理學(xué)研究 [J]. 植物分類與資源學(xué)報, 2011,33(6): 605-614.

    WANG X X, YUE J P, SUN H,etal. Phylogeographical study onEriophytonwallichii(Labiatae) from Alpine scree of Qinghai-Tibetan Plateau [J].PlantDiversity&Resources, 2011,33(6): 605-614.

    [49]DYNESIUS M, JANSSON R. Evolutionary consequences of changes in species' geographical distributions driven by Milankovitch climate oscillations [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2000,97(16): 9 115-9 120.

    [50]于海彬, 張鐿鋰. 青藏高原及其周邊地區(qū)高山植物譜系地理學(xué)研究進(jìn)展 [J]. 西北植物學(xué)報, 2013,33(6): 1 268-1 278.

    YU H B, ZHANG Y L. Advances in phylogeography of alpine plants in the Tibetan Plateau and adjacent regions [J].ActaBotanicaBoreali-OccidentaliaSinica, 2013,33(6): 1 268-1 278.

    [51]RAY N, CURRAT M, EXCOFFIER L. Intra-deme molecular diversity in spatially expanding populations [J].MolecularBiologyandEvolution, 2003,20(1): 76-86.

    [52]EXCOFFIER L. Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model [J].MolecularEcology, 2004,13(4): 853-864.

    [53]HARPENDING H C, BATZER M A, GURVEN M,etal. Genetic traces of ancient demography [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 1998,95(4): 1 961-1 967.

    猜你喜歡
    物種
    物種大偵探
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    生日禮物種草合集
    物種大滅絕
    麗水發(fā)現(xiàn)新物種
    誰在“摧毀”澳大利亞——可怕的物種入侵
    一億年后,地球上可能出現(xiàn)哪些新物種
    回首2018,這些新物種值得關(guān)注
    国产视频一区二区在线看| 一进一出好大好爽视频| 国产高清视频在线播放一区| 国产爱豆传媒在线观看 | 一个人免费在线观看电影 | 嫁个100分男人电影在线观看| 人成视频在线观看免费观看| 黄色视频,在线免费观看| 国产一区二区激情短视频| 久久久国产精品麻豆| 露出奶头的视频| 成人午夜高清在线视频| 久久久久九九精品影院| 亚洲一卡2卡3卡4卡5卡精品中文| 国产视频内射| 日本三级黄在线观看| 人妻丰满熟妇av一区二区三区| 窝窝影院91人妻| 亚洲欧美日韩东京热| 三级男女做爰猛烈吃奶摸视频| 小说图片视频综合网站| 特级一级黄色大片| 怎么达到女性高潮| 啦啦啦韩国在线观看视频| 淫妇啪啪啪对白视频| 精品一区二区三区av网在线观看| 91大片在线观看| 亚洲av成人一区二区三| 国产精品乱码一区二三区的特点| av福利片在线观看| 夜夜看夜夜爽夜夜摸| 成人特级黄色片久久久久久久| 国产高清视频在线播放一区| 国产高清有码在线观看视频 | 最新在线观看一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久| 亚洲人成伊人成综合网2020| 最新在线观看一区二区三区| 午夜成年电影在线免费观看| 国产精品一区二区免费欧美| 国产探花在线观看一区二区| 精品国内亚洲2022精品成人| 亚洲成av人片在线播放无| 又黄又爽又免费观看的视频| 欧美三级亚洲精品| 久久久久亚洲av毛片大全| av天堂在线播放| av欧美777| 国产精品,欧美在线| 欧美日韩中文字幕国产精品一区二区三区| 精品国产乱码久久久久久男人| 99在线视频只有这里精品首页| 91老司机精品| 国产精品综合久久久久久久免费| 国产一区二区三区在线臀色熟女| 久久久精品国产亚洲av高清涩受| 国产亚洲精品综合一区在线观看 | 成熟少妇高潮喷水视频| 又黄又爽又免费观看的视频| 久久精品国产清高在天天线| 别揉我奶头~嗯~啊~动态视频| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 亚洲av美国av| 久久性视频一级片| 免费一级毛片在线播放高清视频| 村上凉子中文字幕在线| av视频在线观看入口| 91麻豆av在线| 国产精品98久久久久久宅男小说| 亚洲18禁久久av| 国产亚洲av嫩草精品影院| 久久热在线av| 成人特级黄色片久久久久久久| 欧美黑人欧美精品刺激| 一进一出抽搐动态| 欧美性长视频在线观看| 男女做爰动态图高潮gif福利片| 99久久无色码亚洲精品果冻| xxx96com| 国产三级在线视频| 久久精品国产亚洲av高清一级| 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 91老司机精品| 级片在线观看| 美女扒开内裤让男人捅视频| 亚洲第一电影网av| 别揉我奶头~嗯~啊~动态视频| 一进一出抽搐动态| 成人三级做爰电影| 精品久久久久久久久久免费视频| 中文资源天堂在线| 在线观看免费视频日本深夜| 男人舔奶头视频| 国产主播在线观看一区二区| 国产午夜精品久久久久久| 香蕉丝袜av| 麻豆av在线久日| 久99久视频精品免费| xxxwww97欧美| 超碰成人久久| 国产欧美日韩一区二区三| 久久久久国内视频| 免费一级毛片在线播放高清视频| 日本五十路高清| 我要搜黄色片| 日韩欧美在线二视频| 中文亚洲av片在线观看爽| 亚洲精品久久成人aⅴ小说| 男人舔女人的私密视频| 免费看a级黄色片| 手机成人av网站| 大型黄色视频在线免费观看| 在线观看一区二区三区| 国产三级黄色录像| 午夜福利在线观看吧| 国产成人av教育| 精品久久久久久久久久久久久| 亚洲av成人一区二区三| 久久久久免费精品人妻一区二区| 亚洲精品一卡2卡三卡4卡5卡| av中文乱码字幕在线| 免费电影在线观看免费观看| 日本a在线网址| 久99久视频精品免费| 在线观看免费午夜福利视频| 国产亚洲精品av在线| 国产真人三级小视频在线观看| 1024视频免费在线观看| 亚洲中文日韩欧美视频| 少妇被粗大的猛进出69影院| 99久久国产精品久久久| 久热爱精品视频在线9| 深夜精品福利| 久久久国产成人精品二区| 婷婷亚洲欧美| 亚洲精品粉嫩美女一区| 人妻久久中文字幕网| 久久国产精品影院| or卡值多少钱| 国产av一区在线观看免费| 欧美日韩黄片免| 国产午夜精品久久久久久| 十八禁网站免费在线| 日韩欧美在线乱码| 亚洲自拍偷在线| 久久久久国产精品人妻aⅴ院| 亚洲一区高清亚洲精品| 99在线视频只有这里精品首页| 可以免费在线观看a视频的电影网站| 国产探花在线观看一区二区| 在线a可以看的网站| 黄色成人免费大全| 一级a爱片免费观看的视频| 久久久国产精品麻豆| 两个人视频免费观看高清| 最近视频中文字幕2019在线8| 久9热在线精品视频| 国产久久久一区二区三区| 亚洲av中文字字幕乱码综合| 国内毛片毛片毛片毛片毛片| 99久久精品国产亚洲精品| 岛国在线观看网站| 免费在线观看完整版高清| 欧美日韩福利视频一区二区| 三级男女做爰猛烈吃奶摸视频| 19禁男女啪啪无遮挡网站| 免费看日本二区| 夜夜看夜夜爽夜夜摸| 露出奶头的视频| 日韩欧美在线二视频| 日韩欧美 国产精品| 制服丝袜大香蕉在线| 国产高清videossex| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| www.精华液| 12—13女人毛片做爰片一| 免费观看人在逋| 久久人妻av系列| 日韩欧美国产一区二区入口| 十八禁人妻一区二区| 久久久久久久久中文| 国产精品香港三级国产av潘金莲| 亚洲成av人片免费观看| 免费在线观看视频国产中文字幕亚洲| 国产三级黄色录像| 国产aⅴ精品一区二区三区波| 此物有八面人人有两片| 国产成人影院久久av| 亚洲人成网站高清观看| 亚洲自拍偷在线| 国产精华一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲激情在线av| 精品一区二区三区四区五区乱码| 亚洲熟女毛片儿| videosex国产| 伦理电影免费视频| 啦啦啦免费观看视频1| 两性夫妻黄色片| 欧美绝顶高潮抽搐喷水| 国产三级黄色录像| 国产在线精品亚洲第一网站| 2021天堂中文幕一二区在线观| 国产精品久久视频播放| 欧美高清成人免费视频www| 在线永久观看黄色视频| 成年免费大片在线观看| 老司机靠b影院| 叶爱在线成人免费视频播放| 国产aⅴ精品一区二区三区波| 女人高潮潮喷娇喘18禁视频| 精品熟女少妇八av免费久了| 欧美日韩精品网址| 黄频高清免费视频| 亚洲av片天天在线观看| 天堂动漫精品| www.精华液| 日韩欧美国产一区二区入口| 99在线视频只有这里精品首页| 又黄又爽又免费观看的视频| 亚洲国产欧美人成| 国产高清videossex| 国产亚洲精品第一综合不卡| 麻豆成人午夜福利视频| 欧美绝顶高潮抽搐喷水| 久久久久久大精品| 日本在线视频免费播放| 午夜免费成人在线视频| 亚洲成人久久爱视频| 变态另类成人亚洲欧美熟女| 亚洲av第一区精品v没综合| 精品国内亚洲2022精品成人| 亚洲,欧美精品.| 桃色一区二区三区在线观看| 999久久久国产精品视频| 色播亚洲综合网| 日韩欧美在线乱码| 亚洲美女视频黄频| 精品一区二区三区四区五区乱码| 身体一侧抽搐| 99riav亚洲国产免费| 老司机深夜福利视频在线观看| 亚洲av成人av| 级片在线观看| 国产亚洲精品第一综合不卡| 少妇裸体淫交视频免费看高清 | 一级作爱视频免费观看| 老鸭窝网址在线观看| 99riav亚洲国产免费| 亚洲国产精品久久男人天堂| 国产精品 欧美亚洲| 韩国av一区二区三区四区| 亚洲精品中文字幕在线视频| 色噜噜av男人的天堂激情| 欧美日韩乱码在线| 国产熟女午夜一区二区三区| 亚洲avbb在线观看| 国产精品乱码一区二三区的特点| 一区福利在线观看| 欧美日韩乱码在线| 一级作爱视频免费观看| aaaaa片日本免费| 又爽又黄无遮挡网站| 999久久久精品免费观看国产| 亚洲aⅴ乱码一区二区在线播放 | 一个人免费在线观看电影 | 日韩精品青青久久久久久| 国产精华一区二区三区| 黄色成人免费大全| 亚洲国产欧美一区二区综合| 老熟妇乱子伦视频在线观看| 91av网站免费观看| 啦啦啦观看免费观看视频高清| 亚洲片人在线观看| 欧美一级毛片孕妇| 精品人妻1区二区| 最近在线观看免费完整版| 国产精品1区2区在线观看.| 少妇裸体淫交视频免费看高清 | 在线免费观看的www视频| 俺也久久电影网| 黑人操中国人逼视频| 欧美日韩中文字幕国产精品一区二区三区| 一进一出好大好爽视频| 舔av片在线| 亚洲美女黄片视频| 久久香蕉国产精品| 亚洲欧美日韩高清专用| 高潮久久久久久久久久久不卡| 成人av一区二区三区在线看| 色综合站精品国产| 一级作爱视频免费观看| 人人妻人人看人人澡| 亚洲成人中文字幕在线播放| 色精品久久人妻99蜜桃| 精品久久久久久久久久久久久| 精品国产亚洲在线| 久久久久久久午夜电影| 国产亚洲欧美98| 激情在线观看视频在线高清| 久久久久久国产a免费观看| 99热这里只有是精品50| 国产精品免费一区二区三区在线| 亚洲天堂国产精品一区在线| 国产私拍福利视频在线观看| 法律面前人人平等表现在哪些方面| 日韩成人在线观看一区二区三区| 亚洲欧美激情综合另类| 亚洲精品在线观看二区| 99久久精品国产亚洲精品| 日本免费一区二区三区高清不卡| 99精品久久久久人妻精品| 欧美三级亚洲精品| 制服诱惑二区| 身体一侧抽搐| 青草久久国产| 999久久久精品免费观看国产| 久久精品91无色码中文字幕| 国产精品免费视频内射| 亚洲黑人精品在线| 亚洲色图av天堂| 欧美成人性av电影在线观看| 国产精品乱码一区二三区的特点| 最近最新中文字幕大全免费视频| 色老头精品视频在线观看| 51午夜福利影视在线观看| 日韩欧美在线乱码| 99国产综合亚洲精品| 国产精品野战在线观看| 麻豆国产97在线/欧美 | 成人国产一区最新在线观看| 国产成人一区二区三区免费视频网站| 神马国产精品三级电影在线观看 | 丁香欧美五月| 国产单亲对白刺激| 一级a爱片免费观看的视频| 欧美黄色淫秽网站| 亚洲成人精品中文字幕电影| 精品久久蜜臀av无| 国产成人系列免费观看| 成人亚洲精品av一区二区| 精品人妻1区二区| 精品一区二区三区四区五区乱码| 2021天堂中文幕一二区在线观| 老司机午夜十八禁免费视频| 亚洲真实伦在线观看| 黄色片一级片一级黄色片| 欧美一区二区精品小视频在线| 国产熟女午夜一区二区三区| 熟女电影av网| 制服丝袜大香蕉在线| 高潮久久久久久久久久久不卡| 丰满人妻熟妇乱又伦精品不卡| 无限看片的www在线观看| 男人舔女人下体高潮全视频| 18禁裸乳无遮挡免费网站照片| 成人一区二区视频在线观看| 国产欧美日韩一区二区三| 18禁黄网站禁片午夜丰满| 欧美成人午夜精品| 亚洲人成77777在线视频| 亚洲国产中文字幕在线视频| 久久久久久免费高清国产稀缺| 亚洲18禁久久av| 高清在线国产一区| 99国产极品粉嫩在线观看| 青草久久国产| 757午夜福利合集在线观看| 波多野结衣高清作品| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 亚洲一区中文字幕在线| 露出奶头的视频| 欧美最黄视频在线播放免费| 免费av毛片视频| 中出人妻视频一区二区| 三级国产精品欧美在线观看 | 午夜福利免费观看在线| 91老司机精品| av欧美777| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 91成年电影在线观看| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 亚洲熟妇中文字幕五十中出| 蜜桃久久精品国产亚洲av| 成人三级黄色视频| 亚洲国产欧美一区二区综合| 国产成+人综合+亚洲专区| 亚洲熟妇熟女久久| 国产aⅴ精品一区二区三区波| 亚洲美女视频黄频| www国产在线视频色| 99国产极品粉嫩在线观看| 99在线视频只有这里精品首页| 亚洲欧美日韩东京热| 久久精品综合一区二区三区| 国产av一区在线观看免费| www日本黄色视频网| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 午夜免费成人在线视频| 中国美女看黄片| 在线观看舔阴道视频| 99精品欧美一区二区三区四区| 狂野欧美激情性xxxx| 免费在线观看影片大全网站| 麻豆av在线久日| 欧美不卡视频在线免费观看 | 国产v大片淫在线免费观看| 国产亚洲精品一区二区www| 99热只有精品国产| 精品午夜福利视频在线观看一区| 少妇熟女aⅴ在线视频| 亚洲中文av在线| 欧美国产日韩亚洲一区| 999久久久国产精品视频| 亚洲av熟女| 国产激情偷乱视频一区二区| 亚洲激情在线av| 女人被狂操c到高潮| 一进一出抽搐gif免费好疼| 久久中文字幕一级| 可以在线观看毛片的网站| 午夜a级毛片| 操出白浆在线播放| 成年版毛片免费区| 俄罗斯特黄特色一大片| 岛国视频午夜一区免费看| 无限看片的www在线观看| 亚洲精品美女久久av网站| 丰满人妻熟妇乱又伦精品不卡| 丁香欧美五月| 99久久精品国产亚洲精品| 天天躁夜夜躁狠狠躁躁| 午夜福利视频1000在线观看| 亚洲专区国产一区二区| 久久久久九九精品影院| 欧美+亚洲+日韩+国产| 少妇的丰满在线观看| 久久精品国产亚洲av香蕉五月| 精品久久久久久久久久免费视频| 一个人免费在线观看的高清视频| 日韩精品中文字幕看吧| 禁无遮挡网站| 国产av在哪里看| 亚洲精品国产一区二区精华液| 美女黄网站色视频| 欧美高清成人免费视频www| 精品免费久久久久久久清纯| 国产亚洲精品久久久久久毛片| 久久久国产成人免费| 国产精品久久久久久亚洲av鲁大| 欧美另类亚洲清纯唯美| www.999成人在线观看| 91在线观看av| 中亚洲国语对白在线视频| 动漫黄色视频在线观看| 亚洲成a人片在线一区二区| 日本 av在线| 精品国产乱子伦一区二区三区| 最近视频中文字幕2019在线8| 婷婷亚洲欧美| 久久久国产成人免费| 久久人人精品亚洲av| 成人18禁在线播放| 久久久久久久午夜电影| 成人国产综合亚洲| 欧美成狂野欧美在线观看| av中文乱码字幕在线| 法律面前人人平等表现在哪些方面| 国产成年人精品一区二区| 9191精品国产免费久久| 欧美黑人巨大hd| 少妇粗大呻吟视频| 91麻豆av在线| 亚洲国产欧洲综合997久久,| 韩国av一区二区三区四区| 最近最新免费中文字幕在线| 日本免费一区二区三区高清不卡| 色av中文字幕| 99在线人妻在线中文字幕| 搞女人的毛片| 国产亚洲av嫩草精品影院| av中文乱码字幕在线| 国内久久婷婷六月综合欲色啪| 波多野结衣高清作品| 香蕉久久夜色| 在线观看午夜福利视频| 人人妻,人人澡人人爽秒播| 午夜成年电影在线免费观看| 国产三级黄色录像| 亚洲av五月六月丁香网| 在线观看免费午夜福利视频| 国产精品香港三级国产av潘金莲| 99精品在免费线老司机午夜| 亚洲七黄色美女视频| 妹子高潮喷水视频| 国产区一区二久久| 亚洲激情在线av| 亚洲国产欧美网| 欧美黑人巨大hd| 国产爱豆传媒在线观看 | 亚洲自拍偷在线| 国产又黄又爽又无遮挡在线| 国产午夜精品久久久久久| 久久久久国产精品人妻aⅴ院| 一本一本综合久久| av国产免费在线观看| 夜夜躁狠狠躁天天躁| 国产精品av久久久久免费| 麻豆成人av在线观看| 精品无人区乱码1区二区| avwww免费| 一级毛片精品| 老鸭窝网址在线观看| 男人舔女人的私密视频| 国产精品 欧美亚洲| videosex国产| 黑人巨大精品欧美一区二区mp4| 国内精品久久久久久久电影| 男女做爰动态图高潮gif福利片| 一级毛片精品| 亚洲精华国产精华精| 亚洲av成人不卡在线观看播放网| 久久久久精品国产欧美久久久| 视频区欧美日本亚洲| 欧美不卡视频在线免费观看 | 亚洲国产看品久久| 国产av麻豆久久久久久久| 亚洲国产精品成人综合色| 午夜视频精品福利| 美女扒开内裤让男人捅视频| 国产精品日韩av在线免费观看| 又黄又爽又免费观看的视频| 可以在线观看的亚洲视频| av在线播放免费不卡| 国产99久久九九免费精品| 亚洲av第一区精品v没综合| 我的老师免费观看完整版| 久久精品国产99精品国产亚洲性色| 亚洲18禁久久av| 成人高潮视频无遮挡免费网站| 欧美人与性动交α欧美精品济南到| 亚洲成人中文字幕在线播放| 成人欧美大片| 99久久无色码亚洲精品果冻| 久久欧美精品欧美久久欧美| 最近最新中文字幕大全电影3| 午夜亚洲福利在线播放| 悠悠久久av| 成人三级做爰电影| 国产三级黄色录像| 小说图片视频综合网站| av片东京热男人的天堂| 丰满人妻一区二区三区视频av | 欧美极品一区二区三区四区| 久久香蕉国产精品| 在线a可以看的网站| 欧美丝袜亚洲另类 | 禁无遮挡网站| 精华霜和精华液先用哪个| 亚洲 欧美一区二区三区| 老鸭窝网址在线观看| 国产av在哪里看| 夜夜躁狠狠躁天天躁| 男女之事视频高清在线观看| 成年版毛片免费区| 日韩高清综合在线| 老熟妇乱子伦视频在线观看| 精华霜和精华液先用哪个| 亚洲avbb在线观看| www日本黄色视频网| 黄色丝袜av网址大全| 亚洲片人在线观看| 成人一区二区视频在线观看| 性色av乱码一区二区三区2| 国产黄片美女视频| 一本一本综合久久| 91字幕亚洲| 老司机午夜十八禁免费视频| 99久久精品热视频| 看免费av毛片| 一个人免费在线观看电影 | 一进一出抽搐gif免费好疼| 国产日本99.免费观看| 国产精品久久久久久精品电影| 丁香六月欧美| 亚洲av第一区精品v没综合| 又紧又爽又黄一区二区| 免费看十八禁软件| 97碰自拍视频| 亚洲第一电影网av| 欧美一区二区国产精品久久精品 | 色综合亚洲欧美另类图片| 麻豆一二三区av精品| av视频在线观看入口| 我的老师免费观看完整版| 精品欧美一区二区三区在线| 丰满人妻一区二区三区视频av | 搡老妇女老女人老熟妇| 国产精品99久久99久久久不卡| av福利片在线| av国产免费在线观看| 亚洲在线自拍视频|