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

    不同海拔梯度下小葉錦雞兒的居群遺傳多樣性響

    2015-02-16 06:21:46黃文達趙學勇李玉霖羅亞勇王少昆潘成臣中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所奈曼沙漠化研究站甘肅蘭州730000甘肅省寒區(qū)旱區(qū)逆境生理與生態(tài)重點實驗室甘肅蘭州730000
    草業(yè)科學 2015年4期
    關(guān)鍵詞:香農(nóng)錦雞兒居群

    黃文達,趙學勇,李玉霖,羅亞勇,王少昆,潘成臣(.中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所奈曼沙漠化研究站, 甘肅 蘭州 730000; 2. 甘肅省寒區(qū)旱區(qū)逆境生理與生態(tài)重點實驗室,甘肅 蘭州 730000)

    植物生產(chǎn)層

    不同海拔梯度下小葉錦雞兒的居群遺傳多樣性響

    黃文達1,2,趙學勇1,李玉霖1,羅亞勇1,王少昆1,潘成臣1
    (1.中國科學院寒區(qū)旱區(qū)環(huán)境與工程研究所奈曼沙漠化研究站, 甘肅 蘭州 730000; 2. 甘肅省寒區(qū)旱區(qū)逆境生理與生態(tài)重點實驗室,甘肅 蘭州 730000)

    以小葉錦雞兒(Caraganamicrophylla)為研究對象,利用ISSR分子標記技術(shù),對不同海拔梯度下小葉錦雞兒6個居群遺傳多樣性進行分析。研究發(fā)現(xiàn),8條ISSR擴增引物得到142條清晰的DNA條帶,其中多態(tài)性條帶為126條,多態(tài)位點百分率為88.73%。用POPGENE軟件計算不同海拔梯度下各居群間和居群內(nèi)的遺傳參數(shù),結(jié)果表明, Nei’s基因多樣性為0.292 9,香農(nóng)信息指數(shù)為0.441 7,表明在物種水平上小葉錦雞兒居群具有較高的遺傳多樣性;低海拔區(qū)域居群的所有遺傳多樣性高于高海拔區(qū)域居群,表明海拔高度增加阻礙了居群間的基因交流,導(dǎo)致小葉錦雞兒居群遺傳多樣性水平降低。

    小葉錦雞兒;海拔梯度;ISSR

    小葉錦雞兒(Caraganamicrophylla)屬豆科錦雞兒屬多年生灌木,在我國多分布于海拔介于200~1 500 m的東北及內(nèi)蒙古、河北、山西、陜西等地的固定和半固定沙地[1]。具有較強的抗風蝕和固氮能力,在植被恢復(fù)過程中具有重要作用[2]。ISSR(Inter Simple Sequence Repeat)分子標記通常用于估算植物種群的遺傳多樣性水平和遺傳結(jié)構(gòu),尤其對基因組序列還不清楚的物種特別有用[3-4]。截至目前,關(guān)于小葉錦雞兒遺傳多樣性的研究主要集中在以下幾個方面: 1)小葉錦雞兒作為錦雞兒屬的一種植物,研究其種間遺傳差異[5-9]; 2)種質(zhì)資源遺傳多樣性[10]; 3)形態(tài)學性狀遺傳多樣性[11];4)內(nèi)蒙古草原區(qū)小葉錦雞兒遺傳多樣性和遺傳分化[12-13]; 5)SSR-PCR體系優(yōu)化及應(yīng)用[14]。但是針對不同海拔梯度下小葉錦雞兒居群的遺傳多樣性研究還未見報道。為此,利用ISSR分子標記技術(shù)研究了小葉錦雞兒居群在不同海拔梯度下的遺傳差異,以期揭示小葉錦雞兒居群間和居群內(nèi)的遺傳多樣性水平,分析其遺傳多樣性與海拔梯度之間的相互關(guān)系,闡明該物種的生態(tài)學適應(yīng)過程,從而為干旱-半干旱地區(qū)生態(tài)系統(tǒng)的持續(xù)穩(wěn)定恢復(fù)提供遺傳學理論依據(jù)。

    1 材料與方法

    1.1 材料

    采樣包括6個居群,67個個體,每一居群的樣本量為10~12,居群內(nèi)個體與個體之間至少間隔30 m(表1)。樣本具有代表性,在每一個呈島嶼狀不連續(xù)分布的居群中,取樣個體代表該地區(qū)所有分布的個體。每個居群的海拔和經(jīng)緯度在采集地中心使用GIS儀(Garmin)測定。

    采集每個個體的幼嫩葉片10 g以上,用變色硅膠迅速干燥后,裝入密封塑料袋帶回實驗室,放于-20 ℃冰箱中保存。

    表1 小葉錦雞兒6個居群的個體數(shù)、生境類型及地理分布Table 1 Location, number of plants and habitats of 6 populations ofCaraganamicrophylla

    1.2 總DNA提取

    采用常規(guī)試劑盒AxyPrep Genomic DNA Mini Kits對所有小葉錦雞兒個體進行總DNA提取。

    1.3 PCR擴增、純化及序列測定

    從100條ISSR擴增引物(加拿大哥倫比亞大學UBC Primer Set No.9, http://www.biotech.ubc.ca/services/NAPS/Primer_Set/Primers.pdf)中篩選出8條條帶清晰、重復(fù)性和穩(wěn)定性好且多態(tài)性條帶相對較多的引物用于全部DNA樣品的PCR擴增。PCR擴增程序為94 ℃預(yù)變性3 min,94 ℃變性30 s,退火30 s,72 ℃延伸1 min,共30個循環(huán),72 ℃延伸5 min,4 ℃保存。擴增反應(yīng)結(jié)束后,以100 bp Ladder marker做參照,用1.8%的瓊脂糖凝膠電泳(含EB 0.5 μg·mol-1)分離,穩(wěn)壓100 V,電泳時間2 h, 凝膠成像系統(tǒng)照相,觀察擴增結(jié)果。

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

    將擴增產(chǎn)物中有條帶的記做“1”,無條帶的記做“0”,只記載清晰易辨的擴增帶,數(shù)據(jù)輸入數(shù)據(jù)矩陣,應(yīng)用POPGENE 1.32軟件,假定這些小葉錦雞兒居群在這些ISSR標記位點處于Hardy-Weinberg平衡狀態(tài),對全部居群和不同海拔梯度下的群體分別計算下列遺傳參數(shù):多態(tài)位點百分率;觀測等位基因數(shù);有效等位基因數(shù);Nei’s基因多樣性;香農(nóng)信息指數(shù)(1972)[15];總基因多樣性;群體內(nèi)基因多樣性;基因分化系數(shù);基因流;遺傳一致度和Nei’s遺傳距離[15]。應(yīng)用SPSS 17.0 軟件分析居群內(nèi)的遺傳多樣性水平與海拔梯度之間的相關(guān)關(guān)系。應(yīng)用ARLEQUIN軟件包分析低海拔區(qū)域、高海拔區(qū)域和區(qū)域內(nèi)居群內(nèi)的遺傳變異[16]。應(yīng)用NTSYS pc 2.02軟件構(gòu)建小葉錦雞兒居群的聚類圖[17]。

    2 結(jié)果

    2.1 多態(tài)位點比率

    通過ISSR擴增,從100條引物中篩選出8條擴增條帶多、信號強、背景清晰且重復(fù)性好的引物用于ISSR-PCR擴增反應(yīng)。經(jīng)統(tǒng)計,這8條擴增引物得到142條清晰的DNA條帶,其中多態(tài)性條帶為126條,所占的百分比是88.73%,擴增的DNA片段長度介于150~2 000 bp。8條引物擴增的DNA條帶數(shù)平均為17.75條。引物857所產(chǎn)生的條帶數(shù)和多態(tài)性條帶最多,分別為22條和21條。引物826所產(chǎn)生的條帶數(shù)和多態(tài)性條帶最少,分別為12條和10條。引物864的多態(tài)百分率最高,達100.00%。引物842的多態(tài)百分率最低,為80.00%(表2)。

    2.2 遺傳多樣性和遺傳分化

    不同海拔區(qū)域6個小葉錦雞兒居群的遺傳多樣性結(jié)果(表3)表明,8條引物在6個居群所產(chǎn)生的多態(tài)位點數(shù)為8~61條;多態(tài)位點百分率為5.63%~42.96%;觀測等位基因數(shù)為1.056 3~1.429 6;有效等位基因數(shù)為1.036 1~1.305 0;Nei’s基因多樣性為0.021 4~0.173 3;香農(nóng)信息指數(shù)為0.031 8~0.253 3。其中居群2的所有遺傳多樣性參數(shù)(多態(tài)位點百分率、觀察等位基因數(shù)、有效等位基因數(shù)、Nei’s基因多樣性和香農(nóng)信息指數(shù))都明顯高于其他居群,居群1的所有遺傳多樣性參數(shù)均低于其他居群。在物種水平上,多態(tài)位點百分率為88.73%,觀察等位基因數(shù)為1.887 3,有效等位基因數(shù)為1.496 7,Nei’s基因多樣性為0.292 9,香農(nóng)信息指數(shù)為0.441 7,該結(jié)果揭示了小葉錦雞兒居群具有豐富的遺傳多樣性。在海拔梯度上,低海拔區(qū)域小葉錦雞兒居群的所有遺傳多樣性參數(shù)都高于高海拔區(qū)域居群,其中重要的遺傳多樣性參數(shù)Nei’s基因多樣性和香農(nóng)信息指數(shù)差異最為明顯。低海拔區(qū)域Nei’s基因多樣性和香農(nóng)信息指數(shù)均明顯高于多年生植物的平均值0.25和0.35,分別為0.280 6和0.415 8。而高海拔區(qū)域Nei’s基因多樣性和香農(nóng)信息指數(shù)均低于多年生植物的平均值,分別為0.193 0和0.288 6。

    表2 用于小葉錦雞兒ISSR擴增的引物及其退火溫度和多態(tài)位點百分率Table 2 Primers sequence, melting temperature and percentage of polymorphism in ISSR analysis ofCaraganamicrophylla

    表3 小葉錦雞兒居群遺傳多樣性Table 3 Genetic diversity indices ofCaraganamicrophyllapopulations

    本試驗研究居群為隨機交配群體,因此,可以假定為Hardy-Weinberg平衡群體。根據(jù)Hardy- Weinberg平衡定律,采用Nei’s基因多樣性指數(shù)對6個小葉錦雞兒居群間基因多樣性和遺傳分化進行比較(表4)。

    試驗結(jié)果表明,海拔梯度對小葉錦雞兒居群的遺傳變異產(chǎn)生了一定程度的影響。低海拔區(qū)域小葉錦雞兒居群大部分遺傳分化參數(shù)都高于高海拔區(qū)域居群。在低海拔區(qū)域,小葉錦雞兒居群的總基因多樣性為0.278 3;群體內(nèi)基因多樣性為0.112 2;在同一區(qū)域中不同居群間的基因分化系數(shù)為0.596 6,即小葉錦雞兒在不同居群間產(chǎn)生了一定程度的遺傳分化,居群間的遺傳變異占總變異的59.66%;居群間的基因流為0.338 0(表4)。高海拔區(qū)域小葉錦雞兒居群的基因流為0.391 0,高于低海拔區(qū)域居群。在居群水平上,小葉錦雞兒居群的總基因多樣性為0.287 6,高于區(qū)域水平上;群體內(nèi)基因多樣性和基因分化系數(shù)分別為0.098 5和0.582 1,都低于低水熱區(qū)域居群(表4)。AMOVA 分析結(jié)果顯示(表5),小葉錦雞兒居群的主要遺傳變異來自于居群間,居群間的遺傳差異百分率達到72.79%(P<0.001),居群內(nèi)遺傳差異百分率為27.21%。

    為了進一步說明小葉錦雞兒各居群間的分化程度,用Popgen 1.32軟件計算了居群間的遺傳距離(表6),發(fā)現(xiàn)居群5和居群4之間的遺傳相似度最大,為0.844 1,遺傳距離最小,為0.169 5。居群1和居群4之間的遺傳相似度最小,為0.646 6,遺傳距離最大,為0.436 0。從UPGMA聚類圖可以看出(圖1),高海拔區(qū)域的小葉錦雞兒居群4、5和6聚為一支,低海拔區(qū)域的居群各聚一支。相關(guān)性分析表明,小葉錦雞兒居群的海拔高度與觀察等位基因的相關(guān)系數(shù)為-0.027,與有效等位基因的相關(guān)系數(shù)為-0.014,與Nei’s基因多樣性的相關(guān)系數(shù)為-0.003,與香農(nóng)指數(shù)的相關(guān)系數(shù)為-0.001,這些遺傳多樣性參數(shù)與海拔均呈負相關(guān),但均未達到顯著水平。

    表4 海拔梯度對小葉錦雞兒居群遺傳分化的影響Table 4 Effects of altitude gradient on the genetic differentiation ofCaraganamicrophyllapopulations

    表5 小葉錦雞兒居群分子變異分析Table 5 Analysis of molecular variance (AMOVA) forCaraganamicrophyllapopulations

    表6 小葉錦雞兒各居群遺傳相似度(對角線以上)和遺傳距離(對角線以下)Table 6 Nei’s genetic identity (above diagonal) and genetic distance (below diagonal) betweenCaraganamicrophyllapopulations

    圖1 基于Nei’s (1972)遺傳距離的小葉錦雞兒居群UPGMA聚類圖Fig.1 UPGMA dendrogram ofCaraganamicrophyllapopulations based on Nei’s (1972) genetic distances from ISSR data

    3 討論與結(jié)論

    物種的遺傳多樣性與其對環(huán)境變化的適應(yīng)性緊密相關(guān),估測遺傳差異是物種保護和生態(tài)恢復(fù)的重要步驟,通過估測遺傳多樣性指數(shù)(如多態(tài)位點百分率等)來衡量物種的遺傳豐富度。本研究結(jié)果顯示,在物種水平上,小葉錦雞兒居群的多態(tài)位點百分率為88.73%,這個值與3種錦雞兒的平均多態(tài)位點百分率(87%)相近[8],說明不同海拔區(qū)域的小葉錦雞兒居群具有較高水平的遺傳多樣性,這與以往的研究結(jié)果一致[5,7,10]。相對于多態(tài)位點百分數(shù)來說,Nei’s基因多樣性和香農(nóng)信息指數(shù)能更好地反映條帶的豐富度和均勻程度,而豐富度和均勻度是衡量多樣性的兩個重要指標,Nei’s基因多樣性還與物種的生物學特性緊密相關(guān)。本研究發(fā)現(xiàn),物種水平上小葉錦雞兒居群的Nei’s基因多樣性是0.292 9,香農(nóng)信息指數(shù)為0.441 7,結(jié)果與楊明博等[18]利用ISSR分子標記研究中間錦雞兒(C.davazamcii)的值相近(Nei’s基因多樣性為0.296 3,香農(nóng)信息指數(shù)為0.454 2),而明顯高于郭宏宇等[12]利用RAPD分子標記對內(nèi)蒙古草原區(qū)小葉錦雞兒的研究結(jié)果(Nei’s基因多樣性為0.184 5,香農(nóng)信息指數(shù)為0.394 2)。這可能是分子標記技術(shù)的不同所導(dǎo)致的。Hamrick和Godt[19]、Nybom和Bartish[20]認為植物的遺傳多樣性水平雖依賴物種的生活史特征,例如繁殖系統(tǒng)、種子擴散機制、地理分布和生活型,但主要受環(huán)境因子的影響[21-22]。小葉錦雞兒是靠蟲媒傳粉的多年生灌木,具有廣泛的生態(tài)幅,除了能異花授粉外,還可以進行自花授粉,具有極強的生命力與競爭力,這些生物學特征都可以使小葉錦雞兒獲得高水平的遺傳多樣性[23-25];從而使得小葉錦雞兒成為中國北方的廣布種和優(yōu)勢種。Ayala和Keger[26]指出,經(jīng)過長期的進化過程,物種對環(huán)境變化的適應(yīng)性主要依靠其種內(nèi)的遺傳變異水平。

    在海拔梯度上,低海拔區(qū)域小葉錦雞兒居群的所有遺傳多樣性參數(shù)都高于高海拔區(qū)域居群,這表明海拔增加阻礙了小葉錦雞兒居群間的基因交流,原因可能是高海拔地區(qū)溫度低,降水豐富,在一定程度上限制了花粉及果實的散布,導(dǎo)致遺傳多樣性水平降低,這與馬丹煒等[27]對巖生植物金發(fā)草(Pogonatherumpaniceum)的遺傳多樣性研究一致。低海拔區(qū)域小葉錦雞兒居群總基因多樣性、群體內(nèi)基因多樣性和基因分化系數(shù)都高于高海拔區(qū)域居群,而居群間的基因流低于高海拔區(qū)域居群, 這說明低海拔區(qū)域小葉錦雞兒居群的遺傳分化程度高于高海拔區(qū)域居群。這主要是因為低海拔區(qū)域為沙漠化草地,生境的片段化導(dǎo)致小葉錦雞兒居群間產(chǎn)生較大的遺傳漂變,這與Huang等[13]對科爾沁沙地小葉錦雞兒居群的研究結(jié)果一致。也有研究證明,生境片段化會影響植物居群間的遺傳多樣性[28-33]。在海拔水平上,低海拔區(qū)域和高海拔區(qū)域小葉錦雞兒居群間的遺傳變異都占總變異的50%以上,其群體內(nèi)分化系數(shù)分別為0.596 6和0.561 2, 這說明海拔梯度對小葉錦雞兒居群的遺傳分化有一定影響。顯著性檢測也顯示小葉錦雞兒居群間的變異所占比例較大,且對總體遺傳變異有顯著影響(P<0.001),說明各地理居群間有一定程度的分化,這與郭宏宇等[12]對小葉錦雞兒和楊明博等[18]對中間錦雞兒遺傳分化的研究不一致,這可能是由物種差異和采樣地的地理差異引起的。本研究發(fā)現(xiàn),居群4和居群5之間的遺傳相似度最大,遺傳距離最小, 這是由于居群4和居群5都處于高海拔區(qū)域,居群間的地理距離較小,居群間會有基因流動。有研究表明,小尺度范圍內(nèi)物種的遺傳多樣性指數(shù)和地理距離之間沒有相關(guān)性[34-35]。本研究結(jié)果發(fā)現(xiàn),居群1和居群4之間的遺傳相似度最小,遺傳距離最大,這主要因為生境差異,居群1來自于丘間低地,其微生境條件比流動沙地好[36-44],而且居群1和居群4之間存在一定的地理距離。根據(jù)UPGMA聚類,小葉錦雞兒居群的遺傳多樣性呈現(xiàn)出一定的分布規(guī)律,聚類結(jié)果將高海拔區(qū)域居群聚在一起。相關(guān)性分析表明,小葉錦雞兒居群的遺傳多樣性參數(shù)與海拔呈負相關(guān),說明海拔高度的變化對小葉錦雞兒的遺傳分化產(chǎn)生了一定的影響,但結(jié)果未達到顯著水平,這與王贊和高洪文[10]對檸條錦雞兒(C.korshinskii)表型性狀與地理因子的相關(guān)性研究結(jié)果一致。

    概而言之,不同海拔梯度影響了小葉錦雞兒居群的遺傳多樣性,而且低海拔區(qū)域小葉錦雞兒居群的遺傳多樣性參數(shù)高于高海拔區(qū)域居群。隨著海拔高度的增加,小葉錦雞兒居群產(chǎn)生了一定的遺傳分化,這表明海拔增加阻礙了居群間的基因交流,導(dǎo)致其遺傳多樣性水平降低[27]。

    [1] 趙一之.小葉、中間和檸條三種錦雞兒的分布式樣及其生態(tài)適應(yīng)[J].生態(tài)學報,2005,25(12):3411-3414.

    [2] 富象乾.內(nèi)蒙古植物志[M].第二版.呼和浩特:內(nèi)蒙古人民出版社,1993:639-648.

    [3] Meloni M,Perini D,F(xiàn)iligheddu R,Binelli G.Genetic variation in five Mediterranean populations ofJuniperusphoeniceaas revealed by inter-simple sequence repeat (ISSR) markers[J].Annals of Botany,2006,97(2):299-304.

    [4] 錢迎倩,馬克平.生物多樣性研究的原理與方法[M].北京:中國科學技術(shù)出版社,1994:13-16.

    [5] 郭強,時永杰,魏臻武,楊占花,盧娟,賈艷瓊.河西走廊14種錦雞兒遺傳多樣性SSR分析[J].草地學報,2008,16(3):227-233.

    [6] 楊九艷,楊劼,楊明博,鞠愛華,渠弼. 鄂爾多斯高原錦雞兒屬藥用植物的ISSR分析[J].中草藥,2006,37(10):1562-1566.

    [7] 盛紅梅,陳拓,安黎哲,鄭曉玲,蒲玲玲,劉亞潔,孫學剛.錦雞兒屬植物的遺傳多樣性及其種間關(guān)系[J].中國沙漠,2005,25(5):697-701.

    [8] 宋俊雙,王贊,高洪文.三種錦雞兒遺傳多樣性ISSR分析(簡報)[J].草地學報,2006,14(4):384-386.

    [9] Ma C C,Gao Y B,Wang J L,Guo H Y.Interspecific transition amongCaraganamicrophylla,C.davazamciiandC.korshinskiialong geographic gradient.Ⅰ. Ecological and RAPD evidence[J].Acta Botanica Sinica,2003,45(10):1218-1227.

    [10] 王贊,高洪文.錦雞兒屬植物種質(zhì)資源遺傳多樣性研究進展[J].植物遺傳資源學報,2008,9(3):397-400.

    [11] 陸景偉,王贊,張新全,裴玉紅,高洪文.小葉錦雞兒天然居群形態(tài)學性狀遺傳多樣性分析[J].湖北農(nóng)業(yè)科學,2008,47(3):328-331.

    [12] 郭宏宇,高玉葆,馬成倉,吳建波,王銀華.內(nèi)蒙古高原小葉錦雞兒遺傳分化的研究[J].西北植物學報,2008,28(1):0072-0077.

    [13] Huang W D,Zhao X Y,Zhao X,Luo Y Q,F(xiàn)eng J,Su N.Genetic variation within the sand-fixation speciesCaraganamicrophylla(Leguminosae) in Horqin sandy land detected by inter-simple sequence repeats analysis[J].Biochemical Systematics and Ecology,2013,51(11):343-348.

    [14] 韓永增,王贊,高洪文.小葉錦雞兒SSR-PCR體系優(yōu)化及應(yīng)用[J].草業(yè)科學,2011,28(3):399-413.

    [15] Yeh F C,Yang R C,Boyle T.POPGENE,Miceosoft Windows-based Freeware for Population Genetic Analysis[M].Edmonton,Canada:Molecular Biology and Biotechnology Center,University of Alberta,1999.

    [16] Schneider S,Roessli D,Excoffier L.Arlequin Ver.2.000:A software for population genetics data analysis.Geneva:Genetics and Biometry Laboratory,University of Geneva,2000.

    [17] Rohl F J.NTSYSpc:Numerical Taxonomy and Multivariate Analysis System,Version 2.02..New York:Exeter Software,Setauket,1998.

    [18] 楊明博,楊劼,楊九艷,臧春鑫.鄂爾多斯高原不同降雨量地區(qū)中間錦雞兒(CaraganadavazamciiSancz)種群的ISSR遺傳分析[J].植物生理學通訊,2006,42(3):411-416.

    [19] Hamrick J L,Godt M J.Allozyme diversity in plant species[A].Brown A D H,Clegg T,Kahler A L,Weir B S.(Eds.),Plant Population Genetics,Breeding and Genetic Resources[M].Sunderland,Massachusetts:Sinauer Associates,1990:43-63.

    [20] Nybom H,Bartish I V.Effects of life history traits and sampling strategies on genetic diversity estimates obtained with RAPD markers in plants[J].Perspectives in Plant Ecology,Evolution and Systematics,2000,3(2):93-114.

    [21] Loveless M D,Hamrick J L.Ecological determinants of genetic structure in plant populations[J].Annual Review of Ecology and Systematics,1984,15(1):65-95.

    [22] Hamrick J L,Godt M J W.Effects of life history traits on genetic diversity in plant species[J].Philosophical Transactions:Biological Sciences,1996,351(9):1291-1298.

    [23] Babbel G R,Selander R K.Genetic variability in edaphically restricted and widespread plant species[J].Evolution,1974,28(4):619-630.

    [24] Pearse D E,Crandall K A.Beyond FST:Analysis of population genetic data for conservation[J].Conservation Genetics,2004,5(5):585-602.

    [25] Ge S,Wang K Q,Hong D Y,Zhang W H,Zu Y G.Comparisons of genetic diversity in the endangeredAdenophoralobophyllaand its widespread congener,A.potaninii[J].Conservation Biology,1999,13(3):509-513.

    [26] Ayala F J,Keger J A.Modern Genetics[M].Second edit Menlo Park,USA:Benjami-Cumings Publishing Company Inc.,1984.

    [27] 馬丹煒,王勝華,羅通,王文國,莊國慶,陳放.環(huán)境因子對鹽生植物金發(fā)草遺傳多樣性的影響[J].中山大學學報(自然科學版),2006,45(2):73-77.

    [28] 李靜,葉萬輝,葛學軍. 生境片段化對植物的遺傳影響[J].中山大學學報(自然科學版),2005,44(2):193-199.

    [29] 黃文達,趙學勇,左小安,李玉強,連杰.不同生境梯度下差不嘎蒿(Artemisiahalodendron)種群ITS基因的序列分析[J].中國沙漠,2013,33(2):554-559.

    [30] Huang W D,Zhao X Y,Zhao X,Zhao H L,Wang S K,Lian J.A combined approach using ISSR and ITS analysis for the characterization ofArtermisiahalodendronfrom Horqin Sandy Land,Northern China[J].Biochemical Systematics and Ecology,2011,39(4):346-351.

    [31] 張穎娟,高瑞霞,李青豐.不同干擾生境中荒漠小灌木紅砂種群遺傳多樣性研究[J].干旱區(qū)資源與環(huán)境,2008,22(3):147-151.

    [32] 陳克霞,王嶸,陳小勇.天然片段生境中山姜(Alpiniajaponica)種群遺傳結(jié)構(gòu)[J].生態(tài)學報,2008,28(6):2480-2485.

    [33] 黃文達,趙學勇,趙昕,趙哈林,連杰,王少昆.分子標記在種群遺傳學研究中的應(yīng)用[J].草業(yè)科學,2010,27(11):115-120.

    [34] 府宇雷,錢吉,馬玉虹,李軍,鄭師章.不同尺度下野大豆種群的遺傳分化[J].生態(tài)學報,2002,22(2):176-184.

    [35] 李珊,蔡宇良,錢增強,趙桂仿.云南金錢槭形態(tài)變異與遺傳變異的相關(guān)性研究[J].生態(tài)學報,2005,24(5):925-931.

    [36] Zuo X A,Zhao H L,Zhao X Y.Plant distribution at the mobile dune scale and its relevance to soil properties and topographic features[J].Environmental Geology,2008,54(5):1111-1120.

    [37] Zuo X A,Zhao H L,Zhao X Y.Spatial pattern and heterogeneity of soil properties in sand dunes under grazing and restoration in Horqin sandy land,Northern China[J].Soil and Tillage Research,2008,99:202-212.

    [38] Burke I C,Lauenroth W K,Riggle R.Spatial variability in soil properties in the shortgrass steppe:the relative importance of topography,grazing,microsite,and plant species in controlling spatial patterns[J].Ecosystems,1999,2(5):422-438.

    [39] Sebastiá M T. Role of topography and soils in grassland structuring at the landscape and community scales[J].Basic and Applied Ecology,2004,5(4):331-346.

    [40] Li Y Q,Zhao H L,Li Y L,Cui J Y,Zuo X A.Soil nitrogen mineralization and nitrification in different habitats,Horqin sandy land[J].Journal of Desert Research,2009,29(3):438-444.

    [41] Li Y Q,Zhao H L,Yi X Y,Zuo X A,Chen Y P.Dynamics of carbon and nitrogen storages in plant-soil system during desertification process in Horqin sandy land[J].Environmental Science,2006,27(4):635-640.

    [42] Holmgren M,Scheffer M.El Nin? as a window of opportunity for the restoration of degraded arid ecosystems[J].Ecosystems,2001,4(2):151-159.

    [43] Su Y Z,Zhao H L,Li Y L,Cui J Y.Influencing mechanisms of several shrubs on soil chemical properties in semi-arid Horqin sandy land,China[J].Arid Land Research and Management,2004,18(3):251-263.

    [44] Li S G,Harazono Y,Zhao H L,He L,Chang X L.Micrometeorological changes following establishment of artificially establishedArtemisiavegetationon desertified sandy land in the Horqin sandy land,China and their implication in regional environmental change[J].Journal of Arid Environments,2002,52(1):101-119.

    (責任編輯 王芳)

    Genetic diversity analysis ofCaraganamicrophyllapopulation in different altitude gradients

    HUANG Wen-da1,2, ZHAO Xue-yong1, LI Yu-lin1, LUO Ya-yong1, WANG Shao-kun1, PAN Cheng-chen1
    (1.Naiman Desertification Research Station, Cold and Arid Region Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000, China;2.Key Laboratory of Stress Physiology and Ecology in Cold and Arid Regions, Lanzhou 730000;China)

    Caraganamicrophyllais the most dominant and constructive shrub species. We evaluated the level of genetic variation within and among populations sampled from two different altitude gradients using inter-simple sequence repeat polymorphism (ISSR) molecular markers. The results showed that eight ISSR primers generated 142 bands, of which 126 (88.73%) were polymorphic. Based on Nei’s Gst value calculated by POPGENE, the gene diversity (h) and Shannon’s information index (I) was 0.292 9 and 0.441 7, respectively indicating considerable genetic variations at the species level. At the altitude gradients, the genetic diversity index in populations from low altitude region was higher than from high altitude region, shows that the gene exchange between populations was blocked by the increase of altitude.

    Caraganamicrophylla; altitude gradients; ISSR

    HUANG Wen-da E-mail:huangwenda2008@163.com

    10.11829j.issn.1001-0629.2014-0379

    2014-08-13 接受日期:2014-01-19

    國家自然科學基金項目(41201561);中科院寒旱所人才基金項目(Y351151001);甘肅省自然科學基金(145RJYA269)

    黃文達(1980-),女,甘肅平川人,助理研究員,博士,主要從事干旱區(qū)植物分子生態(tài)學研究。E-mail:huangwenda2008@163.com

    Q943;S542+.903

    A

    1001-0629(2015)04-0552-08*

    黃文達,趙學勇,李玉霖,羅亞勇,王少昆,潘成臣.不同海拔梯度下小葉錦雞兒的居群遺傳多樣性[J].草業(yè)科學,2015,32(4):552-559.

    HUANG Wen-da,ZHAO Xue-yong,LI Yu-lin,LUO Ya-yong,WANG Shao-kun,PAN Cheng-chen.Genetic diversity analysis ofCaraganamicrophyllapopulation in different altitude gradients[J].Pratacultural Science,2015,32(4):552-559.

    猜你喜歡
    香農(nóng)錦雞兒居群
    燕山地區(qū)不同居群白頭翁種子萌發(fā)的初步評價
    大衛(wèi),不可以
    基于簡單重復(fù)序列間擴增分子標記的金釵石斛遺傳多樣性研究
    8種野生錦雞兒在烏魯木齊的引種試驗
    防護林科技(2020年5期)2020-07-21 01:53:32
    17種錦雞兒屬植物葉片解剖結(jié)構(gòu)及抗旱性分析
    3個居群野生草地早熟禾耐鹽性比較研究
    草地學報(2018年5期)2018-11-07 02:25:00
    不同溫度對變色錦雞兒種子發(fā)芽特征的影響
    校園恩仇錄:小混混和易拉罐女王的故事
    艾米麗的呼嚕
    錦雞兒屬植物種子性狀和幼苗生長特征比較研究
    天堂av国产一区二区熟女人妻| 日韩欧美精品v在线| 亚洲成人中文字幕在线播放| 国产精品久久久久久亚洲av鲁大| 一区二区三区高清视频在线| 精品一区二区三区视频在线观看免费| 婷婷亚洲欧美| 在线观看免费视频日本深夜| 国产极品精品免费视频能看的| 日韩免费av在线播放| 亚洲 国产 在线| 日韩成人在线观看一区二区三区| 久久亚洲精品不卡| 黄片小视频在线播放| 国产三级黄色录像| 免费看美女性在线毛片视频| 夜夜爽天天搞| 婷婷亚洲欧美| 国产高清视频在线播放一区| 一二三四社区在线视频社区8| 亚洲国产欧美一区二区综合| 中文字幕最新亚洲高清| 午夜两性在线视频| 十八禁网站免费在线| 国模一区二区三区四区视频 | 国产视频一区二区在线看| 欧美在线黄色| 99久久无色码亚洲精品果冻| 婷婷丁香在线五月| 天天一区二区日本电影三级| 亚洲自拍偷在线| 噜噜噜噜噜久久久久久91| bbb黄色大片| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久 | 欧美日韩综合久久久久久 | 真人做人爱边吃奶动态| 搞女人的毛片| 变态另类成人亚洲欧美熟女| 欧美+亚洲+日韩+国产| 99久久精品热视频| 我要搜黄色片| 国产精品美女特级片免费视频播放器 | 日韩欧美在线二视频| 亚洲av电影不卡..在线观看| 国产亚洲精品久久久com| 这个男人来自地球电影免费观看| 日韩欧美免费精品| 又大又爽又粗| 非洲黑人性xxxx精品又粗又长| 欧美成人免费av一区二区三区| 99热这里只有是精品50| 好看av亚洲va欧美ⅴa在| 黄色日韩在线| 日韩精品青青久久久久久| 欧美在线黄色| 一二三四社区在线视频社区8| 亚洲av第一区精品v没综合| 禁无遮挡网站| 欧美成人免费av一区二区三区| 国产精品电影一区二区三区| 日韩国内少妇激情av| 久久国产精品影院| 一级黄色大片毛片| 美女黄网站色视频| 宅男免费午夜| 亚洲中文av在线| 精品国产乱码久久久久久男人| 亚洲精品久久国产高清桃花| 日韩av在线大香蕉| 亚洲欧美日韩卡通动漫| 一个人免费在线观看电影 | 精品午夜福利视频在线观看一区| 最新美女视频免费是黄的| 真实男女啪啪啪动态图| 精品久久久久久久人妻蜜臀av| 中文资源天堂在线| 欧美日韩精品网址| 人妻夜夜爽99麻豆av| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美三级三区| 免费看美女性在线毛片视频| 欧美乱妇无乱码| 欧美国产日韩亚洲一区| 免费在线观看影片大全网站| 久久久久久大精品| 美女高潮的动态| 久久久久免费精品人妻一区二区| 久久亚洲精品不卡| 久久精品91无色码中文字幕| 亚洲欧洲精品一区二区精品久久久| 日本成人三级电影网站| 亚洲精品国产精品久久久不卡| 最近在线观看免费完整版| 日本一本二区三区精品| 日韩欧美 国产精品| 一二三四社区在线视频社区8| 后天国语完整版免费观看| 欧美中文日本在线观看视频| 精品久久蜜臀av无| 在线免费观看不下载黄p国产 | 十八禁人妻一区二区| bbb黄色大片| av黄色大香蕉| 超碰成人久久| 日韩免费av在线播放| 一进一出好大好爽视频| 色综合站精品国产| 国产亚洲精品久久久com| 国产欧美日韩一区二区精品| 国产久久久一区二区三区| 亚洲自偷自拍图片 自拍| 亚洲成人精品中文字幕电影| 嫩草影院入口| 99久国产av精品| 国产av麻豆久久久久久久| 丰满人妻熟妇乱又伦精品不卡| 欧美3d第一页| 99久国产av精品| 色av中文字幕| 91在线精品国自产拍蜜月 | 99精品欧美一区二区三区四区| 久久天堂一区二区三区四区| 亚洲熟妇熟女久久| 欧美激情在线99| 狂野欧美白嫩少妇大欣赏| 国产久久久一区二区三区| 欧美日本视频| 夜夜看夜夜爽夜夜摸| 日韩成人在线观看一区二区三区| 91麻豆精品激情在线观看国产| 国产成人福利小说| 亚洲精华国产精华精| 精品无人区乱码1区二区| 中文在线观看免费www的网站| 超碰成人久久| 男女午夜视频在线观看| 久久久久久国产a免费观看| 一夜夜www| 国产精品久久久久久人妻精品电影| 桃红色精品国产亚洲av| 国产成人av教育| 制服人妻中文乱码| 国产毛片a区久久久久| 免费看a级黄色片| 18禁国产床啪视频网站| 色精品久久人妻99蜜桃| 国产v大片淫在线免费观看| 中文字幕久久专区| 久久国产精品影院| 男人舔女人下体高潮全视频| 日韩欧美三级三区| 九九热线精品视视频播放| 色播亚洲综合网| 久久这里只有精品19| 亚洲第一欧美日韩一区二区三区| 美女被艹到高潮喷水动态| 成人国产综合亚洲| 日韩欧美精品v在线| 叶爱在线成人免费视频播放| cao死你这个sao货| 午夜激情欧美在线| 久久热在线av| 国产精品久久久久久久电影 | 国产爱豆传媒在线观看| 日韩成人在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 国产成人aa在线观看| 黑人欧美特级aaaaaa片| 脱女人内裤的视频| 亚洲 国产 在线| 淫秽高清视频在线观看| 99在线视频只有这里精品首页| 国产综合懂色| 亚洲专区国产一区二区| 国产真人三级小视频在线观看| 午夜亚洲福利在线播放| 国产单亲对白刺激| 波多野结衣巨乳人妻| 国产精品久久久久久久电影 | 1024香蕉在线观看| 国产黄色小视频在线观看| 热99re8久久精品国产| 午夜久久久久精精品| 成人永久免费在线观看视频| 久久久久久九九精品二区国产| 精品国产美女av久久久久小说| 亚洲av第一区精品v没综合| bbb黄色大片| 淫妇啪啪啪对白视频| 亚洲精品在线美女| 长腿黑丝高跟| 精品一区二区三区视频在线 | 欧美绝顶高潮抽搐喷水| 在线视频色国产色| 欧美日韩中文字幕国产精品一区二区三区| 精品免费久久久久久久清纯| 亚洲人成电影免费在线| 日本黄色片子视频| 99久久精品一区二区三区| 婷婷亚洲欧美| 日本五十路高清| 99在线人妻在线中文字幕| 1024手机看黄色片| 黑人欧美特级aaaaaa片| 熟女电影av网| 久久天堂一区二区三区四区| 一边摸一边抽搐一进一小说| 亚洲国产精品久久男人天堂| 757午夜福利合集在线观看| 色在线成人网| 天天添夜夜摸| xxxwww97欧美| 中文资源天堂在线| 亚洲自拍偷在线| 亚洲av第一区精品v没综合| 少妇的丰满在线观看| 久久精品亚洲精品国产色婷小说| 国产蜜桃级精品一区二区三区| 午夜激情欧美在线| 亚洲一区二区三区色噜噜| 国产黄色小视频在线观看| 精品一区二区三区av网在线观看| 少妇丰满av| 国产成人aa在线观看| 黄色日韩在线| 国产高清有码在线观看视频| 看片在线看免费视频| 国产不卡一卡二| 欧美日韩福利视频一区二区| 一级黄色大片毛片| av视频在线观看入口| 国产精品野战在线观看| 日本精品一区二区三区蜜桃| 日本在线视频免费播放| 五月伊人婷婷丁香| aaaaa片日本免费| 中文亚洲av片在线观看爽| 国产精品女同一区二区软件 | 噜噜噜噜噜久久久久久91| 午夜久久久久精精品| 精品久久久久久,| av片东京热男人的天堂| 久久精品91蜜桃| 亚洲欧美一区二区三区黑人| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| 亚洲欧美激情综合另类| 国产高潮美女av| 国产精品 国内视频| 99久久99久久久精品蜜桃| 精品乱码久久久久久99久播| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品久久久久久毛片| 午夜亚洲福利在线播放| 黑人操中国人逼视频| 色吧在线观看| 欧美另类亚洲清纯唯美| 欧美av亚洲av综合av国产av| x7x7x7水蜜桃| 日韩欧美 国产精品| 天天躁日日操中文字幕| 免费看光身美女| 欧美黑人欧美精品刺激| 99热6这里只有精品| 在线国产一区二区在线| 999久久久国产精品视频| 在线观看一区二区三区| 久久久久国内视频| 国产综合懂色| 亚洲中文字幕日韩| 亚洲精品美女久久av网站| 美女 人体艺术 gogo| 欧美乱色亚洲激情| 美女午夜性视频免费| 中国美女看黄片| 精品一区二区三区av网在线观看| 好男人电影高清在线观看| 在线观看美女被高潮喷水网站 | 午夜免费成人在线视频| 51午夜福利影视在线观看| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 精品久久久久久久毛片微露脸| 真人一进一出gif抽搐免费| 亚洲国产精品999在线| 琪琪午夜伦伦电影理论片6080| 久久九九热精品免费| 国产人伦9x9x在线观看| 久久精品人妻少妇| 日本免费a在线| 丰满人妻熟妇乱又伦精品不卡| 国产午夜福利久久久久久| 国产亚洲欧美在线一区二区| 精品无人区乱码1区二区| 757午夜福利合集在线观看| 亚洲,欧美精品.| 午夜激情福利司机影院| 国产成+人综合+亚洲专区| 国产伦精品一区二区三区视频9 | 日韩欧美三级三区| 亚洲人与动物交配视频| 超碰成人久久| 国产免费av片在线观看野外av| 成人av在线播放网站| 日韩中文字幕欧美一区二区| 18禁裸乳无遮挡免费网站照片| 欧美不卡视频在线免费观看| 色播亚洲综合网| 免费在线观看视频国产中文字幕亚洲| 欧美性猛交╳xxx乱大交人| 一级毛片精品| 女生性感内裤真人,穿戴方法视频| 黑人欧美特级aaaaaa片| 深夜精品福利| 九色国产91popny在线| 久久精品夜夜夜夜夜久久蜜豆| 国产69精品久久久久777片 | www.熟女人妻精品国产| 九九久久精品国产亚洲av麻豆 | 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区精品| 悠悠久久av| 免费观看的影片在线观看| av天堂中文字幕网| 最近最新免费中文字幕在线| 亚洲专区国产一区二区| 久久久久国内视频| 亚洲中文字幕日韩| 成年人黄色毛片网站| 小说图片视频综合网站| a在线观看视频网站| 免费无遮挡裸体视频| 天天添夜夜摸| 成人永久免费在线观看视频| 18禁美女被吸乳视频| 国产精品99久久久久久久久| 性色av乱码一区二区三区2| 日韩欧美在线乱码| 国产精品av久久久久免费| 国产成人精品久久二区二区免费| 露出奶头的视频| 一区二区三区国产精品乱码| 舔av片在线| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 亚洲av中文字字幕乱码综合| 中文字幕熟女人妻在线| 香蕉久久夜色| 日日干狠狠操夜夜爽| 亚洲va日本ⅴa欧美va伊人久久| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 制服丝袜大香蕉在线| 国产成+人综合+亚洲专区| 极品教师在线免费播放| 国产伦在线观看视频一区| 国产精品自产拍在线观看55亚洲| 亚洲国产日韩欧美精品在线观看 | 国产97色在线日韩免费| 无遮挡黄片免费观看| 国产亚洲欧美98| 日本与韩国留学比较| 国产亚洲欧美98| 午夜福利18| 国产单亲对白刺激| 一a级毛片在线观看| 国产成人欧美在线观看| 午夜a级毛片| 欧美3d第一页| 一a级毛片在线观看| 亚洲国产中文字幕在线视频| 国产蜜桃级精品一区二区三区| 国内久久婷婷六月综合欲色啪| 久久这里只有精品中国| 国产 一区 欧美 日韩| 97超视频在线观看视频| 国产探花在线观看一区二区| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 国产一区二区在线av高清观看| 精品国内亚洲2022精品成人| 综合色av麻豆| 一个人看视频在线观看www免费 | 美女扒开内裤让男人捅视频| 在线永久观看黄色视频| 中文在线观看免费www的网站| 国产私拍福利视频在线观看| 亚洲在线自拍视频| 亚洲国产欧美人成| 欧美日韩国产亚洲二区| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 精品久久久久久久末码| 国产亚洲精品久久久久久毛片| 国产精华一区二区三区| 国产爱豆传媒在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲自拍偷在线| 精品久久蜜臀av无| 黄色女人牲交| 国产精品一区二区免费欧美| 19禁男女啪啪无遮挡网站| 国产激情偷乱视频一区二区| 免费看美女性在线毛片视频| 亚洲av五月六月丁香网| 免费观看的影片在线观看| 日韩欧美在线二视频| 久久精品国产清高在天天线| 一级毛片女人18水好多| 天天添夜夜摸| 国产毛片a区久久久久| 色精品久久人妻99蜜桃| 999精品在线视频| 国产美女午夜福利| 欧美av亚洲av综合av国产av| 午夜成年电影在线免费观看| 亚洲,欧美精品.| 精品国内亚洲2022精品成人| av在线天堂中文字幕| 在线十欧美十亚洲十日本专区| 男女视频在线观看网站免费| 亚洲 欧美 日韩 在线 免费| 国产成人av教育| 色在线成人网| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 看黄色毛片网站| 在线观看美女被高潮喷水网站 | 精品久久蜜臀av无| 他把我摸到了高潮在线观看| 成人一区二区视频在线观看| 香蕉丝袜av| 婷婷亚洲欧美| 国产亚洲av嫩草精品影院| 欧美一区二区精品小视频在线| 国产成人aa在线观看| 91字幕亚洲| 欧美黄色淫秽网站| 国产伦精品一区二区三区视频9 | 两性午夜刺激爽爽歪歪视频在线观看| 免费人成视频x8x8入口观看| 国产亚洲精品久久久久久毛片| 久久精品91无色码中文字幕| 免费看光身美女| 国内毛片毛片毛片毛片毛片| 免费人成视频x8x8入口观看| 日韩欧美在线乱码| 亚洲av日韩精品久久久久久密| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 免费看日本二区| 人人妻,人人澡人人爽秒播| 国产麻豆成人av免费视频| 老汉色∧v一级毛片| 毛片女人毛片| 精品国产亚洲在线| 亚洲五月婷婷丁香| 变态另类成人亚洲欧美熟女| 久久天躁狠狠躁夜夜2o2o| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 可以在线观看毛片的网站| 久久中文看片网| 国产成人欧美在线观看| 欧美+亚洲+日韩+国产| 久9热在线精品视频| 久久精品91蜜桃| 丝袜人妻中文字幕| 久久久久国产精品人妻aⅴ院| 亚洲九九香蕉| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 怎么达到女性高潮| 精品乱码久久久久久99久播| 高清在线国产一区| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片 | 97超视频在线观看视频| 日韩欧美国产在线观看| 久久99热这里只有精品18| 欧洲精品卡2卡3卡4卡5卡区| 国产精品av视频在线免费观看| 黄色丝袜av网址大全| 亚洲精品国产精品久久久不卡| 搞女人的毛片| 亚洲五月婷婷丁香| 日本与韩国留学比较| www.熟女人妻精品国产| 久久久国产精品麻豆| 欧美不卡视频在线免费观看| av女优亚洲男人天堂 | 日本一本二区三区精品| 听说在线观看完整版免费高清| h日本视频在线播放| 成人av在线播放网站| 91av网一区二区| 一本精品99久久精品77| 天天添夜夜摸| 在线观看免费视频日本深夜| 99热6这里只有精品| 国产真人三级小视频在线观看| 亚洲人成网站高清观看| 少妇裸体淫交视频免费看高清| 久久香蕉精品热| 老司机午夜十八禁免费视频| 九九久久精品国产亚洲av麻豆 | 亚洲av美国av| 香蕉久久夜色| 美女cb高潮喷水在线观看 | 国产高清videossex| 国产精品日韩av在线免费观看| 亚洲最大成人中文| 99在线视频只有这里精品首页| 在线视频色国产色| 欧美一级a爱片免费观看看| 偷拍熟女少妇极品色| 熟女电影av网| 欧美丝袜亚洲另类 | 99国产精品99久久久久| 母亲3免费完整高清在线观看| 悠悠久久av| 色av中文字幕| 国产精品九九99| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 美女大奶头视频| 淫妇啪啪啪对白视频| 午夜亚洲福利在线播放| 不卡一级毛片| 中出人妻视频一区二区| 长腿黑丝高跟| 欧美不卡视频在线免费观看| 女警被强在线播放| 国产伦精品一区二区三区视频9 | 在线观看66精品国产| 中文字幕精品亚洲无线码一区| av欧美777| 波多野结衣巨乳人妻| 午夜福利在线观看吧| 亚洲精品在线观看二区| 国产又黄又爽又无遮挡在线| 亚洲av电影在线进入| 成年版毛片免费区| 美女 人体艺术 gogo| 一进一出好大好爽视频| 国产午夜精品论理片| 亚洲成av人片免费观看| 在线视频色国产色| 脱女人内裤的视频| 老鸭窝网址在线观看| 久久欧美精品欧美久久欧美| 一本精品99久久精品77| 欧美日韩乱码在线| 12—13女人毛片做爰片一| 精品99又大又爽又粗少妇毛片 | 精品福利观看| 99热6这里只有精品| 亚洲激情在线av| 中文资源天堂在线| 香蕉国产在线看| 日韩 欧美 亚洲 中文字幕| 国产激情欧美一区二区| 色精品久久人妻99蜜桃| 男女午夜视频在线观看| 亚洲欧美日韩高清在线视频| 久久久久久久精品吃奶| 此物有八面人人有两片| av黄色大香蕉| 国产黄色小视频在线观看| 成人av一区二区三区在线看| svipshipincom国产片| 黄色视频,在线免费观看| 九九在线视频观看精品| 国产单亲对白刺激| 国产午夜精品久久久久久| 一二三四在线观看免费中文在| 国产视频内射| 免费高清视频大片| 欧美不卡视频在线免费观看| 亚洲人成伊人成综合网2020| 久久久久久久久久黄片| 国产毛片a区久久久久| 久久这里只有精品19| 欧美黄色片欧美黄色片| 婷婷亚洲欧美| 日本撒尿小便嘘嘘汇集6| 亚洲18禁久久av| 亚洲色图av天堂| 国产日本99.免费观看| 国产高清视频在线观看网站| 国产精品一区二区三区四区免费观看 | 国产97色在线日韩免费| 国产欧美日韩精品亚洲av| 国产精品98久久久久久宅男小说| 国产亚洲欧美98| 久久精品影院6| 变态另类成人亚洲欧美熟女| 首页视频小说图片口味搜索| 19禁男女啪啪无遮挡网站| 欧美日韩综合久久久久久 | 免费在线观看日本一区| 91在线精品国自产拍蜜月 | 国产精品av久久久久免费| 老汉色∧v一级毛片| 色尼玛亚洲综合影院| 欧美日韩乱码在线| 男人舔女人的私密视频| xxx96com| 男女视频在线观看网站免费| 亚洲美女黄片视频| 丁香六月欧美| a在线观看视频网站|