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

    裸果木種群遺傳多樣性及其與土壤因子的關(guān)聯(lián)性研究

    2017-10-16 01:30:07徐振朋宛濤蔡萍張曉明伊衛(wèi)東宛詣超
    關(guān)鍵詞:果木樣地甘肅

    徐振朋,宛濤*,蔡萍,張曉明,伊衛(wèi)東,宛詣超

    裸果木種群遺傳多樣性及其與土壤因子的關(guān)聯(lián)性研究

    徐振朋1,宛濤1*,蔡萍2,張曉明1,伊衛(wèi)東3,宛詣超1

    1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)草原與資源環(huán)境學(xué)院,內(nèi)蒙古 呼和浩特 010019;2. 錫林郭勒職業(yè)學(xué)院,內(nèi)蒙古 錫林浩特 026000;3. 內(nèi)蒙古農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,內(nèi)蒙古 呼和浩特 010019

    為深入了解裸果木(Gymnocarpos przewalskii)遺傳基礎(chǔ)、斑塊狀分布的種群遺傳結(jié)構(gòu)及遺傳多樣性與土壤因子間的相關(guān)性,對(duì)分布于西北荒漠區(qū)的12個(gè)裸果木種群遺傳多樣性及土壤養(yǎng)分及進(jìn)行分析,為裸果木遺傳多樣性的地理變異趨勢(shì)預(yù)測(cè)及種質(zhì)資源的科學(xué)管理和保育提供依據(jù)。結(jié)果表明,裸果木居群遺傳多樣性較高,其中,Nei’s基因多樣性指數(shù)(He)變化范圍為0.235~0.269,Shannom信息指數(shù)(I)變化范圍為0.351~0.405,居群內(nèi)具有豐富的遺傳變異,基因分化系數(shù)(Gst)為0.279。基于遺傳一致度對(duì)種群進(jìn)行PCA分析,12個(gè)種群裸果木可分為4大類:甘肅肅北縣種群為分布的核心區(qū),單獨(dú)為一類;以河西走廊中西部為中心分化2類,甘肅瓜州及甘肅柳園鎮(zhèn)為河西走廊西部中心種群,與向西南和向西延伸的甘肅阿克塞縣和新疆伊州區(qū)種群親緣關(guān)系較近,聚為一類;甘肅金塔縣、甘肅高臺(tái)縣、甘肅清泉鄉(xiāng)、甘肅肅州區(qū)為中心的河西走廊東部分布聚為一類;向東和向東南延伸的內(nèi)蒙古阿右旗、寧夏沙坡頭區(qū)域和甘肅民勤縣聚為一類。依據(jù)裸果木分布區(qū)土壤養(yǎng)分含量進(jìn)行主成分聚類,將種群劃分為3大類:寧夏沙坡頭樣地中有機(jī)質(zhì)、磷、氮和水分含量較高,被單獨(dú)劃為一類;甘肅肅北縣等種群因土壤養(yǎng)分含量處于中等水平,聚為一類;內(nèi)蒙古阿右旗等種群因鉀元素含量高,被劃分為一類。將遺傳特征參數(shù)與土壤養(yǎng)分含量進(jìn)行相關(guān)性分析,土壤中氮含量及全鹽含量對(duì)種群遺傳多樣性水平有一定影響,處于貧瘠土壤上種群機(jī)體能夠更加有效地利用土壤中的限制性養(yǎng)分。

    裸果木;土壤因子;遺傳多樣性;遺傳結(jié)構(gòu)

    裸果木(Gymnocarpos przewalskii)為石竹科(Caryophyllaceae)裸果木屬(Gymnocarpos)多年生小半灌木,為國(guó)家Ⅰ級(jí)保護(hù)植物,起源于第三紀(jì),古地中海殘遺成分(傅立國(guó),1992;劉媖心,1985)。該植物主要分布于甘肅河西走廊、內(nèi)蒙古西部、寧夏中西部和新疆東部(馬松梅等,2010)。生境較為單一,生長(zhǎng)在海拔為800~2500 m的荒漠區(qū)的干河床、干河道、山前洪積扇及礫石戈壁灘等地。盡管裸果木種群生長(zhǎng)稀疏、生物量相對(duì)較低,但它對(duì)維持脆弱的荒漠生態(tài)環(huán)境,遏制荒漠化,維護(hù)當(dāng)?shù)厣鷳B(tài)系統(tǒng)穩(wěn)定具有不可替代的作用(劉媖心,1985)。然而,因群落結(jié)構(gòu)簡(jiǎn)單、生態(tài)系統(tǒng)脆弱、自然地理環(huán)境變遷及人類活動(dòng)干擾等影響,近年來(lái)裸果木種群數(shù)量銳減,種群斑塊狀分布格局趨于明顯(馬松梅等,2010)。目前對(duì)于裸果木的研究主要集中于種群生態(tài)學(xué)(柴永青等,2010;王立龍等,2015)、分類學(xué)(Petrusson et al.,1996;Oxelman et al.,2002)、繁殖生物學(xué)(汪之波等,2004;劉生龍等,1995;李新蓉等,2016)等方面。遺傳多樣性是生物體攜帶遺傳信息的總和,是生物多樣性的重要組成部分,是維持生物繁殖活力及與環(huán)境相互適應(yīng)的基礎(chǔ)。物種適應(yīng)環(huán)境產(chǎn)生的遺傳變異體現(xiàn)出它們應(yīng)對(duì)不同環(huán)境的能力,在應(yīng)對(duì)長(zhǎng)期以來(lái)環(huán)境壓力的影響,其主要是通過改變遺傳基礎(chǔ)以適應(yīng)環(huán)境(Gienapp et al.,2008)。因此,資源異質(zhì)性是種群遺傳多樣性研究的基礎(chǔ)(Yeaman et al.,2006;Huang et al.,2016)。在全球變化的大背景下,研究環(huán)境因子對(duì)瀕危物種遺傳多樣性的影響具有重要的生態(tài)學(xué)與保護(hù)生物學(xué)意義(Taberlet et al.,2012)。種群或物種遺傳多樣性越高,對(duì)環(huán)境變化的適應(yīng)能力越強(qiáng),越容易擴(kuò)展其分布范圍(Barrière et al.,2005)。種群遺傳結(jié)構(gòu)的研究不僅有助于闡明自然條件下生物變異與環(huán)境的關(guān)系,而且能夠解析物種親緣關(guān)系及演替進(jìn)程。探討物種的遺傳多樣性及其與生態(tài)因子間的相關(guān)性,評(píng)價(jià)其遺傳變異及其產(chǎn)生的相應(yīng)格局,將有助于更好地保護(hù)和利用種質(zhì)資源。為了深入了解裸果木遺傳基礎(chǔ)、斑塊狀分布的種群遺傳結(jié)構(gòu)、遺傳多樣性與土壤因子間的相關(guān)性,以分布中國(guó)西部干旱區(qū)的裸果木為研究對(duì)象,對(duì)其種群的遺傳變異與分化進(jìn)行研究,分析裸果木種群遺傳變異與土壤因子的關(guān)系,探討影響裸果木種群遺傳變異的主要生態(tài)因子,為裸果木遺傳多樣性的地理變異趨勢(shì)預(yù)測(cè)及其種質(zhì)資源的科學(xué)管理和保育提供依據(jù)。

    1 材料與方法

    1.1 樣品采集

    2015年5月,依據(jù)裸果木分布特點(diǎn)對(duì)分布于內(nèi)蒙古、寧夏、新疆及甘肅河西走廊區(qū)域 12個(gè)裸果木種群土壤樣品及葉片進(jìn)行取樣,每個(gè)種群采集20個(gè)植株的健康葉片,每株植株間隔距離10 m以上,樣地分布信息見表1和圖1。

    在每個(gè)樣地內(nèi)選取50 m×50 m樣方3個(gè),每個(gè)樣方采集5個(gè)10~20 cm土層土樣,等量混合均勻,室內(nèi)自然風(fēng)干,去除粗枝、石塊等,過2目篩后進(jìn)行土壤理化性質(zhì)分析。

    1.2 樣品分析方法

    1.2.1 土壤樣品分析

    土壤pH采用電位法測(cè)定,有機(jī)質(zhì)采用硫酸-重鉻酸鉀法測(cè)定;全氮采用半微量凱氏法測(cè)定;全磷采用鉬藍(lán)比色法測(cè)定;全鉀采用火焰光度法測(cè)定;速效氮采用擴(kuò)散皿法測(cè)定;速效磷采用碳酸氫鈉法測(cè)定;速效鉀采用乙酸銨提取法測(cè)定(鮑士旦,1999)。

    表1 試驗(yàn)區(qū)裸果木種群樣地信息Table 1 The test site and sample plot’s information of the population of Gymnocarpos przewalskii

    1.2.2 遺傳多樣性的ISSR分析

    采用改進(jìn)CTAB法進(jìn)行DNA提取(李金璐等,2013);利用紫外分光光度法檢測(cè)DNA濃度,并用1.0%的瓊脂糖凝膠電泳檢測(cè)其完整性,-20 ℃保存?zhèn)溆?。依?jù)正交實(shí)驗(yàn)建立裸果木ISSR-PCR反應(yīng)體系,在 25 μL反應(yīng)體系中 dNTP濃度為 0.30 mmol?L-1,引物為 0.50 μmol?L-1,Mg2+為 2.50 mmol·L-1,Taq DNA 酶為 1.00 U?μL-1,DNA 質(zhì)量濃度為40.00 ng?μL-1時(shí)效果最佳。對(duì)100個(gè)ISSR引物進(jìn)行篩選,從中選取擴(kuò)增條帶清晰、重復(fù)性較好的引物,所選引物條件如表 2所示。采用UBC(GA)8C引物進(jìn)行PCR擴(kuò)增條件篩選,初步確定 ISSR-PCR的最佳擴(kuò)增程序?yàn)椋?4 ℃預(yù)變性 3 min;94 ℃變性30 s,54 ℃復(fù)性45 s,72 ℃延伸1.5 min,35個(gè)循環(huán);72 ℃延伸10 min,4 ℃保存。

    圖1 裸果木種群調(diào)查樣點(diǎn)空間分布示意圖Fig. 1 Observed spatial clustering of sample plots in the survey region of Gymnocarpos przewalskii

    表2 用于檢測(cè)裸果木居群遺傳分化的引物Table 2 The primers in analyzing genetic variations of Gymnocarpos przewalskii populations

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

    運(yùn)用SPSS 18.0軟件對(duì)不同種群土壤因子進(jìn)行主成分分析,POPGEN 32軟件對(duì)種群間和種群內(nèi)Nei’s基因多樣性指數(shù)(He)、總的遺傳多樣性(Ht)、居群內(nèi)的遺傳多樣性(Hs)、居群內(nèi)的遺傳多樣性(Hs)、Shannom信息指數(shù)(I)、等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、多態(tài)位點(diǎn)(PPL)等遺傳參數(shù)進(jìn)行分析(李明等,2013);基因分化系數(shù)(Gst)由公式Gst=1-Hs/Ht算得;居群間基因流(Nm)由公式Nm=(1-Gst)/4Gst算得;采用NTSYS-pc 2.1軟件中的遺傳一致度參數(shù)構(gòu)建種群 PCA主成分關(guān)系圖,將遺傳參數(shù)與土壤因子進(jìn)行相關(guān)分析。

    2 結(jié)果分析

    2.1 裸果木種群樣地土壤因子分析

    不同裸果木種群樣地土壤養(yǎng)分含量見表 3,土壤含水量在 3.93%~8.77%之間;有機(jī)質(zhì)變化范圍為2.9~6.3 g?kg-1;全氮變化范圍為 0.13~0.43 g?kg-1;速效氮變化范圍為26.00~51.07 mg?kg-1;全磷含量變化范圍為 0.31~0.86 g?kg-1,速效磷變化范圍為3.07~5.07 mg?kg-1。其中,寧夏沙坡頭種群樣地中的有機(jī)質(zhì)、土壤含水量、磷元素及氮元素均表現(xiàn)為最高。全鉀變化范圍為13.6~20.17 g?kg-1,速效鉀變化范圍為 144.67~261.33 mg?kg-1;甘肅高臺(tái)縣種群土壤鉀元素含量較高。土壤含鹽量變化范圍為2.17~4.11 g?kg-1,其中,甘肅阿克塞縣及甘肅肅北縣種群土壤全鹽含量較高。土壤pH變化范圍為6.33~8.42,寧夏沙坡頭區(qū)、甘肅瓜州種群的裸果木種群土壤為弱酸性,其余大部分區(qū)域土壤為堿性或弱堿性,表明堿性或弱堿性土壤更適合裸果木生長(zhǎng)。

    對(duì)不同裸果木種群樣地土壤養(yǎng)分含量進(jìn)行主成分分析(PCA),提取了前3個(gè)主成分,累計(jì)貢獻(xiàn)率為88.70%。PC1主要解釋了除全鉀和速效鉀以外磷、氮及土壤含水量等含量;PC2解釋全鉀及速效鉀等含量;PC3主要解釋了全鹽含量及pH等。同時(shí),繪制土壤養(yǎng)分三維聚類圖,如圖2所示,甘肅肅北縣、甘肅瓜州縣、甘肅金塔縣和甘肅阿克塞縣等樣地聚為一類,為第1類;內(nèi)蒙古阿右旗、新疆伊州區(qū)、甘肅民勤縣、甘肅省高臺(tái)縣、甘肅肅州區(qū)、甘肅清泉鄉(xiāng)、甘肅柳園鎮(zhèn)等樣地聚為一類,為第2類;寧夏沙坡頭樣地單獨(dú)為第3類。對(duì)土壤主成分因子進(jìn)行分析可知,在第1類樣地土壤中,有機(jī)質(zhì)、磷、氮和水分含量處于中等水平;在第3類樣地中,即寧夏沙坡頭樣地中,有機(jī)質(zhì)、磷、氮和水分含量最高的一類。在第2類樣地中,分異特征為鉀和速效鉀高被劃分一類,而第一、三主成分的得分范圍介于第1類和第3類樣地之間。

    表3 不同裸果木種群土壤的理化性質(zhì)Table 3 The soil physicochemical properties in different distribution areas of Gymnocarpos przewalskii

    圖2 不同種群土壤因子主成分得分聚類圖Fig. 2 Scores plots from PCA soils from different originsPC 1為第1主成分;PC 2為第2主成分;PC 3為第3主成分PC 1 The first PC; PC 2 The second PC; PC 3 the third PC

    2.2 裸果木種群遺傳特征分析

    依據(jù)PCR反應(yīng)體系對(duì)100條ISSR引物進(jìn)行篩選,最終篩選出 15條擴(kuò)增條帶清晰、重復(fù)性好的引物,用于不同種群裸果木的遺傳多樣性分析。對(duì)12個(gè)種群的240個(gè)裸果木植株進(jìn)行擴(kuò)增,總計(jì)擴(kuò)增出條帶數(shù)為 143,具有多態(tài)位條帶數(shù) 127,引物多態(tài)位平均比例為88.81%。不同種群裸果木遺傳參數(shù)如表4所示,平均多態(tài)位點(diǎn)百分率(PPL)為67.37%,較高的有甘肅阿克塞縣種群(81.12%),其次分別為甘肅肅北縣居群、寧夏沙坡頭、甘肅柳園鎮(zhèn)、甘肅清泉鄉(xiāng)及甘肅金塔縣等種群,最低為內(nèi)蒙古阿右旗種群(61.54%)。居群中有效等位基因數(shù)為(Ne)1.403~1.462,較高的為甘肅柳園鎮(zhèn)種群,最低為甘肅金塔縣種群。

    居群間He變化范圍為0.235~0.269,平均值為0.248;居群間I變化范圍為0.351~0.405,平均值為0.369。甘肅省阿克塞與甘肅肅北縣種群的 He及 I值均高于其他種群,具有豐富的遺傳多樣性??偟倪z傳多樣性(Ht)為 0.345,而居群內(nèi)的遺傳多樣性(Hs)為0.248,群體間遺傳分化系數(shù)Gst為0.279,即居群間的遺傳變異占總變異的27.90%,表明裸果木遺傳多樣性主要存在于居群內(nèi),而居群間也存在一定的分化。居群間基因流(Nm)為 1.292,表明居群之間存在著遺傳交換。

    圖3 ISSR標(biāo)記對(duì)不同種群裸果木PCA二維圖Fig. 3 Two-dimensional PCA plot based on ISSR from Gymnocarpos przewalskii

    表4 基于ISSR標(biāo)記的不同種群裸果木遺傳多態(tài)性參數(shù)Table 4 Genetic diversity from different populations based on ISSR of Gymnocarpos przewalskii

    對(duì)12個(gè)種群裸果木進(jìn)行PCA分析,前3個(gè)主成分累積解釋了總變異的93.01%?;谶z傳一致度對(duì)種群進(jìn)行主成分分類,結(jié)果如圖3所示,遺傳距離與地理分布有明顯關(guān)系,以河西走廊中西部為中心分化為2類,甘肅肅北縣種群為裸果木分布的核心區(qū),單獨(dú)為一類。甘肅金塔縣、甘肅高臺(tái)縣、甘肅清泉鄉(xiāng)、甘肅肅州區(qū)為中心的河西走廊東部分布親緣關(guān)系較近,聚為一類;東和向東南延伸的內(nèi)蒙古阿右旗、寧夏沙坡頭區(qū)域和甘肅民勤縣親緣關(guān)系較近,聚為一類;甘肅瓜州及甘肅柳園鎮(zhèn)為核心的河西走廊西部中心種群,與向西南和向西延伸的甘肅阿克塞縣和新疆伊州區(qū)種群親緣關(guān)系較近,聚為一類。

    2.3 遺傳特征與土壤因子相關(guān)性分析

    對(duì)不同種群裸果木 ISSR標(biāo)記的各遺傳多樣性參數(shù)與土壤因子進(jìn)行相關(guān)分析(表5),結(jié)果表明,ISSR標(biāo)記的各遺傳多樣性參數(shù)與其所處的土壤因子有一定的相關(guān)性。土壤因子中氮元素及速效氮與Ne、He、I呈顯著負(fù)相關(guān);土壤中鹽分含量與遺傳參數(shù)相關(guān)性較強(qiáng),與Na、Ne、PPL呈顯著正相關(guān),與He、I呈極顯著正相關(guān),其余各土壤因子與遺傳參數(shù)之間無(wú)顯著相關(guān)性。

    3 討論

    3.1 裸果木總體遺傳多樣性

    Nybom et al.(2004)認(rèn)為期望雜合度(He)和Shannon信息指數(shù)(I)遺傳信息更為重要,其注重物種在居群水平的遺傳多樣性,而 Hamrick et al.(1996)認(rèn)為多態(tài)位點(diǎn)百分率(P)和期望雜合度(He)遺傳信息比較重要,其注重物種水平的多樣性。本研究結(jié)果表明,不同種群裸果木平均多態(tài)位點(diǎn)百分率(PPL)為67.37%,居群間He平均值為0.248,I平均值為0.369。裸果木總的遺傳多樣性(Ht)為0.345,而居群內(nèi)的遺傳多樣性(Hs)為0.248,綜合分析,裸果木具有較高的遺傳多樣性,維持著較高的遺傳變異。裸果木具有較高遺傳水平的原因可能與其本身的生物學(xué)特性有關(guān),一方面裸果木以有性繁殖為主,世代周期較長(zhǎng),為異花授粉植物具蟲媒及風(fēng)媒傳粉;另一方面裸果木為古地中海孑遺種,在中國(guó)西北旱區(qū)“避難所”中第三紀(jì)殘遺植物得以保留;末次間冰期裸果木分布范圍較廣,但受冰期的影響數(shù)量驟減趨于瀕危(吳征鎰等,1983),受第四紀(jì)冰川的影響當(dāng)時(shí)大部分植物消失,然而裸果木在強(qiáng)烈的氣候變遷過程中仍保存至今,能夠開拓適應(yīng)地來(lái)應(yīng)對(duì)環(huán)境變化,遺留下的種群通過積累豐富的遺傳基礎(chǔ)而適應(yīng)環(huán)境的變化,同四合木(Tetraenamongolica)及長(zhǎng)葉紅砂(Reaumuria trigyna)(張穎娟等,2008;Ge et al.,2003)等孑遺植物相同,表現(xiàn)出較高的遺傳多樣性水平。因此,可推斷出裸果木遺傳多樣性水平較高與其生物學(xué)特性、演化歷史和分布區(qū)生態(tài)因子的特殊性分不開。Ma et al.(2012)采用cpDNA標(biāo)記對(duì)裸果木的遺傳結(jié)構(gòu)進(jìn)行研究表明,中國(guó)西北地區(qū)的干旱和沙漠的擴(kuò)展造成裸果木片段化。本研究結(jié)果表明裸果木基因流為 1.292,基因流較大的可能原因是裸果木繁育系統(tǒng)屬兼性異交類型,以異交為主,自交為輔,異交導(dǎo)致裸果木的基因流增大。同時(shí),裸果木兼有風(fēng)媒和蟲媒傳粉,盛花期裸果木花期開花量較大(李新蓉等,2016),傳粉昆蟲較多,使得基因在較多個(gè)體間交流。亦有研究表明,高原隆升形成的主要山系阻礙了植物群體間的基因流,促進(jìn)了種群間分化(Liu et al.,2013;Li et al.,2011)。本研究中,裸果木種群在西北的分布主要沿著天山山脈延伸分布于新疆部分區(qū)域,而在東部的分布沿著祁連山分布于河西走廊東部及陰山山脈的內(nèi)蒙古西部。甘肅肅北縣及甘肅阿克塞縣為裸果木分布中心分布區(qū),由中心分布向兩邊擴(kuò)展,使得種群間的基因流變小,遺傳水平亦隨著分布范圍擴(kuò)大而變低;其次,種群沿著祁連山、阿爾金山北側(cè)分布,種群間不存在山系間的隔離,種群間的地理隔離因子不強(qiáng)烈,可能存在較大的基因交流。

    3.2 土壤因子對(duì)裸果木遺傳特征的影響

    Levins(1969)研究了自然種群的遺傳結(jié)構(gòu)與微生境土壤因子的關(guān)系,建立了環(huán)境幅度模型,主張分子遺傳變異與環(huán)境變異相關(guān),高水平的遺傳變異是植物響應(yīng)顯著異質(zhì)環(huán)境的適應(yīng)策略。裸果木種群分布地土壤有機(jī)質(zhì)含量較低,地表多礫石,一些表層土中具有龜裂化、漆皮化和碳酸鹽表聚集等現(xiàn)象。分布區(qū)土壤類型為灰漠土、灰棕漠土及棕漠土,土壤為中性或弱堿性。主成分分析顯示,裸果木植物個(gè)體在聚類特征上呈現(xiàn)從東到西分布的特征,與水分、土壤養(yǎng)分變化的分布狀況相吻合。研究發(fā)現(xiàn),在裸果木種群分布中,土壤鹽分含量作用較大,這與王孝安(1997)對(duì)荒漠戈壁古老植被土壤養(yǎng)分的研究結(jié)果一致。古南大路及古地中海起源孑遺植物,在新的地區(qū)的生存和繁衍需要有“預(yù)適性”過程,強(qiáng)耐鹽性島嶼植物可演化為旱生植物,這也從側(cè)面證明極耐旱植物裸果木在一定程度上對(duì)鹽分土壤具有適應(yīng)性。本研究也發(fā)現(xiàn)生長(zhǎng)在甘肅省肅北縣、甘肅瓜州縣、甘肅金塔縣和甘肅省阿克塞等樣地一類土壤類型的裸果木多生長(zhǎng)在干河床、沖積扇等地,大部分為裸果木種群分布的核心區(qū)。這些分布區(qū)的土壤含水量較低,氮含量低,鹽分含量高,土壤較貧瘠,但是卻具有較高的遺傳多樣性。與荒漠中的叉毛蓬(Petrosimonia sibirica)(王祎玲等,2006)、紅砂(Reaumuria songarica)等種群存在相似的適應(yīng)機(jī)制,表明荒漠區(qū)的裸果木處于貧瘠土壤上種群機(jī)體能夠更加有效地利用土壤中的限制性養(yǎng)分。大部分物種遺傳多樣性與土壤中的總氮量呈顯著正相關(guān)(李丹等,2001),而本研究裸果木種群的遺傳多樣性與樣地土壤的總氮及速效氮含量存在的顯著負(fù)相關(guān)性表明,相對(duì)于他未能顯現(xiàn)顯著相關(guān)性的土壤因子,裸果木對(duì)土壤中的氮元素具有一定的選擇適應(yīng)性。在土壤養(yǎng)分供應(yīng)不足但細(xì)胞內(nèi)滲透壓調(diào)節(jié)較高的環(huán)境下,在長(zhǎng)期貧瘠土壤環(huán)境下,裸果木為了生存,不斷地增強(qiáng)應(yīng)對(duì)環(huán)境的能力,積累了豐富的遺傳變異,以更好地維持著亞種群的穩(wěn)定性和正常的代謝功能。

    表5 裸果木遺傳參數(shù)與土壤因子相關(guān)性分析Table 5 Correlation analysis between parameters of genetic diversity and soil for Gymnocarpos przewalskii

    4 結(jié)論

    物種的遺傳多樣性水平在一定程度上體現(xiàn)著這個(gè)物種適應(yīng)環(huán)境的能力,制約著物種適應(yīng)性進(jìn)化水平,同時(shí)物種的遺傳多樣性水平可為其種群現(xiàn)狀和保護(hù)價(jià)值的評(píng)估提供重要信息。遺傳多樣性水平是瀕危植物保護(hù)的核心問題,裸果木種群遺傳多樣性研究表明,在自然選擇中能夠產(chǎn)生較為豐富的遺傳變異,表明裸果木在遺傳多樣性上沒有表現(xiàn)出瀕危特性,導(dǎo)致其瀕危的原因可能是其本身的繁殖特性引起的,今后對(duì)該種的研究需加強(qiáng)其生物學(xué)特點(diǎn)、生殖生態(tài)、傳粉生態(tài)等方面的研究,以揭示其瀕危的機(jī)理。研究發(fā)現(xiàn),甘肅阿克塞縣及甘肅肅北縣分布有Nei’s遺傳多樣性指數(shù)(He)和Shannom信息指數(shù)(I)較高的裸果木種群,且這些地區(qū)也表現(xiàn)出較高的多態(tài)位點(diǎn)百分率(PPL),表明這些居群具有較為豐富的遺傳多樣性,也反映出這些地區(qū)的裸果木種群適應(yīng)力較強(qiáng)。

    實(shí)地踏查發(fā)現(xiàn),甘肅阿克塞及甘肅肅北縣是裸果木分布的中心地帶,推測(cè)其冰期避難所可能位于河西走廊的中部地區(qū),其種群在西北的分布主要沿著天山山脈延伸分布于新疆部分區(qū)域,在東部的分布沿著祁連山分布于河西走廊東部及陰山山脈的內(nèi)蒙古西部。裸果木分布由中心向兩邊擴(kuò)展,使得種群間的基因流變小,導(dǎo)致遺傳水平隨著分布范圍擴(kuò)大變低。

    土壤養(yǎng)分分析表明,堿性或弱堿性土壤更適合裸果木生長(zhǎng)。土壤氮含量及全鹽含量對(duì)種群遺傳特征有一定影響,氮素含量少的瘠地表現(xiàn)出高的遺傳多樣性,表明處于貧瘠土壤上裸果木種群能夠更加有效地利用土壤中的限制性養(yǎng)分。生態(tài)系統(tǒng)作為一個(gè)整體,各環(huán)境因子之間相互聯(lián)系、共同作用,本文僅分析了裸果木不同種群間部分土壤因子與遺傳多樣性的關(guān)系,有關(guān)多因子的綜合作用還有待進(jìn)一步研究。

    BARRIèRE A, FéLIX M A. 2005. High local genetic diversity and low out crossing rate in Caenorhabditis elegans natural populations [J]. Current Biology, 15(13): 1176-1184.

    GE X J, YU Y, ZHAO N X, et al. 2003. Genetic variation in the endangered Inner Mongolia endemic shrub Tetraena mongolica Maxim(Zygophyllaceae) [J]. Biological Conservation, 111(3): 427-434.

    GIENAPP P, TEPLITSKY C, ALHO J S, et al. 2008. Climate change and evolution:disentangling environmental and genetic responses [J].Molecular ecology, 17(1): 167-178.

    HAMRICK J L, GODT M J W. 1996. Effects of life history traits on genetic diversity in plant species [J]. Philosophical Transactions of the Royal Society of London B:Biological Sciences, 351(1345): 1291-1298.

    HUANG Y, TRAN I, AGRAWAL A F. 2016. Does genetic variation maintained by environmental heterogeneity facilitate adaptation to novel selection? [J]. The American Naturalist, 188(1): 27-37.

    LEVINS R. 1969. Some demographic and genetic consequences of environmental heterogeneity for biological control [J]. Bulletin of the Entomological society of America, 15(3): 237-240.

    LI Y, ZHAI S N, QIU Y X, et al. 2011. Glacial survival east and west of the‘Mekong-Salween Divide’ in the Himalaya-Hengduan Mountains region as revealed by AFLPs and cpDNA sequence variation in Sinopodophyllum hexandrum (Berberidaceae) [J]. Molecular Phylogenetics and Evolution, 59(2): 412-424.

    LIU J, MOLLER M, PROVAN J, et al. 2013. Geological and ecological factors drive cryptic speciation of yews in a biodiversity hotspot [J].New Phytologist, 199(4): 1093-1108.

    MA S M, ZhANG M L, SANDERSON S C. 2012. Phylogeography of the rare Gymnocarpos przewalskii (Caryophyllaceae): indications of multiple glacial refugia in north-western China [J]. Australian Journal of Botany, 60(1): 20-31.

    NYBOM H. 2004. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants [J]. Molecular Ecology, 13(5): 1143-1155.

    OXELMAN B, AHLGREN B, THULIN M. 2002. Circumscription and phylogenetic relationships of Gymnocarpos (Caryophyllaceae Paronychioideae) [J]. Edinburgh Journal of Botany, 59(2): 221-237.

    PETRUSSON L, THULIN M. 1996. Taxonomy and biogeography of Gymnocarpos (Caryophyllaceae) [J]. Edinburgh Journal of Botany,53(1): 1-26.

    TABERLET P, ZIMMERMANN N E, ENGLISCH T, et al. 2012. Genetic diversity in widespread species is not congruent with species richness in alpine plant communities [J]. Ecology Letters, 15(12): 1439-1448.

    YEAMAN S, JARVIS A. 2006. Regional heterogeneity and gene flow maintain variance in a quantitative trait within populations of lodgepole pine [J]. Proceedings of the Royal Society of London B: Biological Sciences, 273(1594): 1587-1593.

    鮑士旦. 1999. 土壤農(nóng)化分析[M]. 北京: 中國(guó)農(nóng)業(yè)出版社.

    柴永青, 曹致中, 蔡卓山, 等. 2010. 肅北地區(qū)稀有植物裸果木種群的空間分布格局[J]. 草業(yè)學(xué)報(bào), 19(5): 239-249.

    傅立國(guó). 1992. 中國(guó)植物紅皮書[M]. 北京: 科學(xué)出版社: 68-371.

    李丹, 彭少麟. 2001. 三個(gè)不同海拔梯度馬尾松種群的遺傳多樣性及其與生態(tài)因子的相關(guān)性[J]. 生態(tài)學(xué)報(bào), 21(3): 415-421.

    李金璐, 王碩, 于婧, 等. 2013. 一種改良的植物 DNA 提取方法[J]. 植物學(xué)報(bào), 48(1): 72-78.

    李明, 王樹香, 高寶嘉. 油松天然次生林居群遺傳多樣性及與產(chǎn)地地理氣候因子的關(guān)聯(lián)分析[J]. 生態(tài)學(xué)報(bào), 2013, 33(12): 3602-3610.

    李新蓉, 唐欣, 付文潔. 2016. 荒漠孑遺植物裸果木的花部特征及繁育系統(tǒng)[J]. 生態(tài)學(xué)雜志, 35(10): 2592-2598.

    劉生龍, 郭志中, 王理德, 等. 1995. 裸果木引種栽培試驗(yàn)研究[J]. 中國(guó)沙漠, 15(1): 88-91.

    劉媖心. 1985. 中國(guó)沙漠植物志[M]. 北京: 科學(xué)出版社.

    馬松梅, 張明理, 張宏祥, 等. 2010. 利用最大熵模型和規(guī)則集遺傳算法模型預(yù)測(cè)孑遺植物裸果木的潛在地理分布及格局[J]. 植物生態(tài)學(xué)報(bào), 34(11): 1327-1335.

    汪之波, 高清祥, 孫繼周. 2004. 稀有植物裸果木的組織培養(yǎng)及植株再生[J]. 西北植物學(xué)報(bào), 24(7): 1319-1321.

    王立龍, 王亮, 張麗芳, 等. 2015. 不同生境下瀕危植物裸果木種群結(jié)構(gòu)及動(dòng)態(tài)特征[J]. 植物生態(tài)學(xué)報(bào), 39(10): 980-989.

    王孝安. 1997. 甘南瑪曲植物群落的多元分析與環(huán)境解釋[J]. 生態(tài)學(xué)報(bào),17(1): 61-65.

    王祎玲, 畢潤(rùn)成, 閆桂琴, 等. 2006. 新疆阜康綠洲荒漠過渡帶叉毛蓬種群遺傳多樣性及環(huán)境適應(yīng)性研究[J]. 西北植物學(xué)報(bào), 26(6):1133-1141.

    吳征鎰, 王荷生. 1983. 中國(guó)自然地理——植物地理(上冊(cè))[M]. 北京:科學(xué)出版社: 1-125.

    張穎娟, 王玉山. 2008. 瀕危小灌木長(zhǎng)葉紅砂種群的遺傳多樣性[J]. 林業(yè)科學(xué), 44(12): 43-47.

    Abstract: To study the genetic basis Gymnocarpos przewalskii, population distribution of plaque structure and correlation between genetic diversity and soil factors. In this study, the analysis of soil nutrient and population genetic diversity of 12 Gymnocarpos przewalskii population in the northwest desert area provided scientific guidance for predicting the geographical variation of genetic diversity of Gymnocarpos przewalski, and the scientific management and conservation of germplasm resources provide evidence.The results showed that the genetic diversity of Gymnocarpos przewalskii was higher, the Nei's gene diversity (Ne) varied from 0.2356 to 0.2697, and the Shannon Information Index (I) varied from 0.3509 to 0.4057. Genetic variation was abundant while the genetic differentiation coefficient (Gst) was 0.2790. PCA was used to analysis the Genetic consistency of population, and 12 populations of Gymnocarpos przewalskii can be divided into four classes. The core area of the population distribution of Subei town of Gansu can be classified as one class. To the central and western part of the Hexi Corridor as the center of differentiation grouped in to 2, Guazhou of Gansu and Liuyuan town of Gansu were the Hexi Corridor West Center population, with the southwest and westward extension of Akesai county Gansu and Yizhou district of Xinjiang population had the closer relationship were clustered as a class. The population was divided into three groups according to the principal component analysis (PCA) of soil nutrition content in the distribution area. Shapotou area of Ningxia which contained higher organic matter, phosphorus, nitrogen and water content were classified as a class. Subei country of Gansu and others which had moderate level of soil nutrient were classified as a class. A Youqi of Inner Mongolia and other populations which had high levels of potassium were classified as a class. The correlation between genetic characteristics and soil nutrient content was analyzed. The genetic diversity of the population was affected by soil nitrogen content and total salt content. The population in the barren soil could make effectively use of the restricted nutrient in the soil.

    Key words: Gymnocarpos przewalskii; soil factors; genetic diversity; genetic structure

    Study on the Relevance between Genetic Diversity of Gymnocarpos przewalskii Populations and Soil Factors

    XU Zhenpeng1, WAN Tao1*, CAI Ping2, ZHANG Xiaoming1, YI Weidong3, WAN Yichao1

    1. College of Grassland, Resources and Environment, Inner Mongolia Agricultural University, Hohhot 010019, China;2. Xilingol Vocational College, Xilinhot 026000, China; 3. College of Agricultural, Inner Mongolia Agricultural University, Hohhot 010019, China

    10.16258/j.cnki.1674-5906.2017.09.003

    Q 948; X176

    A

    1674-5906(2017)09-1473-07

    徐振朋, 宛濤, 蔡萍, 張曉明, 伊衛(wèi)東, 宛詣超. 2017. 裸果木種群遺傳多樣性及其與土壤因子的關(guān)聯(lián)性研究[J]. 生態(tài)環(huán)境學(xué)報(bào), 26(9): 1473-1479.

    XU Zhenpeng, WAN Tao, CAI Ping, ZHANG Xiaoming, YI Weidong, WAN Yichao. 2017. Study on the relevance between genetic diversity of gymnocarpos przewalskii populations and soil factors [J]. Ecology and Environmental Sciences, 26(9): 1473-1479.

    內(nèi)蒙古自然科學(xué)基金項(xiàng)目(2016MS0324);內(nèi)蒙古自治區(qū)研究生科研創(chuàng)新項(xiàng)目(B20151012919)

    徐振朋(1988年生),男,博士研究生,研究方向牧草種質(zhì)資源。E-mail: 947464263@qq.com*通信作者:宛濤(1959年生),男,教授,博士生導(dǎo)師,主要從事牧草種質(zhì)資源研究。E-mail: wantao425@sohu.com

    2017-08-16

    猜你喜歡
    果木樣地甘肅
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    昆明市主要綠化樹種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    洄流
    行走甘肅
    甘肅卷
    大美甘肅,我們來(lái)了
    甘肅卷
    角果木屬的化學(xué)成分及生物活性研究
    二氧化硅改良頂果木木材表面性質(zhì)研究
    69精品国产乱码久久久| 亚洲成国产人片在线观看| 国产人伦9x9x在线观看| 97人妻天天添夜夜摸| 国产精品麻豆人妻色哟哟久久| 国产一卡二卡三卡精品 | 欧美久久黑人一区二区| 一级a爱视频在线免费观看| 老熟女久久久| 黄色视频不卡| 制服诱惑二区| 久久久久久久久免费视频了| 一边摸一边抽搐一进一出视频| 亚洲精品一区蜜桃| 免费在线观看完整版高清| 亚洲精品日韩在线中文字幕| 丰满乱子伦码专区| 777久久人妻少妇嫩草av网站| 少妇精品久久久久久久| 国产免费福利视频在线观看| 国产成人欧美在线观看 | 亚洲图色成人| 国产精品av久久久久免费| 又大又黄又爽视频免费| 午夜免费鲁丝| 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 亚洲精品久久久久久婷婷小说| 亚洲精品国产一区二区精华液| 一本久久精品| 国产精品久久久av美女十八| 纵有疾风起免费观看全集完整版| 欧美另类一区| 啦啦啦啦在线视频资源| 久久鲁丝午夜福利片| 国产精品成人在线| 久久综合国产亚洲精品| 激情视频va一区二区三区| 1024香蕉在线观看| 国产精品无大码| 久久午夜综合久久蜜桃| 中文乱码字字幕精品一区二区三区| 少妇精品久久久久久久| 电影成人av| 高清av免费在线| 亚洲欧美激情在线| 免费在线观看黄色视频的| 亚洲欧美成人综合另类久久久| 免费久久久久久久精品成人欧美视频| 中文乱码字字幕精品一区二区三区| av女优亚洲男人天堂| av在线app专区| 精品一区二区三卡| 日韩伦理黄色片| 男女国产视频网站| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 各种免费的搞黄视频| 午夜福利视频精品| 又黄又粗又硬又大视频| 久热爱精品视频在线9| 国产爽快片一区二区三区| 成年动漫av网址| 亚洲av电影在线进入| 国产极品天堂在线| 最黄视频免费看| 国产精品欧美亚洲77777| 黄色毛片三级朝国网站| 久久久久国产精品人妻一区二区| 晚上一个人看的免费电影| 国产亚洲av片在线观看秒播厂| 久久精品熟女亚洲av麻豆精品| 国产毛片在线视频| 99国产综合亚洲精品| 日韩一本色道免费dvd| 国产精品一国产av| 国产97色在线日韩免费| 男女高潮啪啪啪动态图| 热re99久久国产66热| av女优亚洲男人天堂| 中文字幕制服av| 久热这里只有精品99| 十八禁网站网址无遮挡| 亚洲第一区二区三区不卡| av在线观看视频网站免费| 精品一区二区免费观看| 欧美中文综合在线视频| 免费观看性生交大片5| 老司机影院成人| 国产精品久久久久久精品古装| 午夜福利视频在线观看免费| a级片在线免费高清观看视频| 波野结衣二区三区在线| 超色免费av| 人人澡人人妻人| 欧美成人精品欧美一级黄| 国产在线免费精品| 精品久久蜜臀av无| 麻豆精品久久久久久蜜桃| 久久av网站| 欧美精品一区二区免费开放| 99久久99久久久精品蜜桃| 乱人伦中国视频| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 夜夜骑夜夜射夜夜干| 亚洲第一av免费看| 中文欧美无线码| 午夜免费鲁丝| 无遮挡黄片免费观看| 男女下面插进去视频免费观看| www.精华液| 亚洲五月色婷婷综合| 亚洲第一青青草原| 国产熟女午夜一区二区三区| 人妻人人澡人人爽人人| 18禁国产床啪视频网站| 99热国产这里只有精品6| 毛片一级片免费看久久久久| 男女午夜视频在线观看| 亚洲精品自拍成人| 中文字幕av电影在线播放| 亚洲国产精品国产精品| 制服丝袜香蕉在线| 男女边吃奶边做爰视频| 一区二区三区激情视频| 国产av精品麻豆| 777米奇影视久久| 国产精品99久久99久久久不卡 | 青青草视频在线视频观看| 晚上一个人看的免费电影| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 老司机在亚洲福利影院| 成年人免费黄色播放视频| 久久久久精品久久久久真实原创| 最近手机中文字幕大全| 在线观看免费视频网站a站| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 久久国产精品男人的天堂亚洲| 亚洲激情五月婷婷啪啪| 好男人视频免费观看在线| 搡老岳熟女国产| 男女边摸边吃奶| 日韩成人av中文字幕在线观看| 久久久精品区二区三区| 国产精品av久久久久免费| 免费久久久久久久精品成人欧美视频| 亚洲精品国产av蜜桃| 免费观看人在逋| 色综合欧美亚洲国产小说| 麻豆乱淫一区二区| 国产av码专区亚洲av| 高清av免费在线| 性色av一级| 精品酒店卫生间| 中文字幕色久视频| 国产精品久久久久久精品古装| 免费高清在线观看日韩| 成人18禁高潮啪啪吃奶动态图| 国产国语露脸激情在线看| 最近手机中文字幕大全| 精品人妻熟女毛片av久久网站| 日本wwww免费看| 国产精品熟女久久久久浪| 久久天堂一区二区三区四区| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 国产探花极品一区二区| 少妇的丰满在线观看| 亚洲欧美一区二区三区国产| 不卡视频在线观看欧美| 国产男女内射视频| 久久人人97超碰香蕉20202| 制服丝袜香蕉在线| 亚洲精品av麻豆狂野| 无遮挡黄片免费观看| 国产黄色免费在线视频| 国产成人a∨麻豆精品| 国产精品久久久av美女十八| 亚洲天堂av无毛| 热re99久久国产66热| 青春草视频在线免费观看| 精品免费久久久久久久清纯 | 男男h啪啪无遮挡| 国产日韩欧美在线精品| 在线观看免费午夜福利视频| 最近2019中文字幕mv第一页| 国产在线免费精品| av天堂久久9| 日韩制服丝袜自拍偷拍| 美女大奶头黄色视频| 亚洲 欧美一区二区三区| 国产亚洲av片在线观看秒播厂| 18在线观看网站| 男男h啪啪无遮挡| 日韩精品免费视频一区二区三区| 男人操女人黄网站| 国产片特级美女逼逼视频| 99re6热这里在线精品视频| 少妇的丰满在线观看| 性色av一级| 中文字幕另类日韩欧美亚洲嫩草| 国产99久久九九免费精品| 无遮挡黄片免费观看| 只有这里有精品99| 日本欧美国产在线视频| 岛国毛片在线播放| 久久久久人妻精品一区果冻| 久久国产精品男人的天堂亚洲| 亚洲熟女精品中文字幕| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 青春草亚洲视频在线观看| 免费黄网站久久成人精品| 中文天堂在线官网| 天天躁日日躁夜夜躁夜夜| 国产又色又爽无遮挡免| 色综合欧美亚洲国产小说| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 久久97久久精品| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 欧美在线一区亚洲| av国产久精品久网站免费入址| 亚洲精品一二三| 日本vs欧美在线观看视频| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说| 久久精品久久精品一区二区三区| 宅男免费午夜| 国产精品人妻久久久影院| 妹子高潮喷水视频| 午夜福利网站1000一区二区三区| 午夜91福利影院| 蜜桃在线观看..| 色吧在线观看| 国产欧美日韩一区二区三区在线| 午夜福利网站1000一区二区三区| 大片免费播放器 马上看| 国产一卡二卡三卡精品 | 三上悠亚av全集在线观看| 成人毛片60女人毛片免费| 久久 成人 亚洲| 亚洲人成网站在线观看播放| 欧美日韩亚洲国产一区二区在线观看 | 国产精品秋霞免费鲁丝片| 国产精品免费视频内射| 亚洲七黄色美女视频| 久久 成人 亚洲| 亚洲视频免费观看视频| 各种免费的搞黄视频| 国产精品蜜桃在线观看| 欧美日韩成人在线一区二区| 国产片特级美女逼逼视频| 久久精品人人爽人人爽视色| 久久97久久精品| 久热这里只有精品99| 操美女的视频在线观看| 免费日韩欧美在线观看| 欧美黄色片欧美黄色片| 操出白浆在线播放| 国产精品久久久久成人av| 精品一品国产午夜福利视频| 日韩av在线免费看完整版不卡| 男女边吃奶边做爰视频| 亚洲欧美激情在线| 国产伦理片在线播放av一区| 国产黄频视频在线观看| 在线观看国产h片| 亚洲三区欧美一区| 国产探花极品一区二区| 99热国产这里只有精品6| 高清不卡的av网站| 亚洲国产精品国产精品| 亚洲第一区二区三区不卡| 各种免费的搞黄视频| 一级a爱视频在线免费观看| 在线天堂最新版资源| 9色porny在线观看| 熟女av电影| 大片电影免费在线观看免费| 国产成人系列免费观看| 各种免费的搞黄视频| 伊人久久大香线蕉亚洲五| 久久精品亚洲av国产电影网| 秋霞在线观看毛片| 国产精品久久久久久人妻精品电影 | 美女中出高潮动态图| 人人妻人人澡人人看| 性色av一级| 天天躁夜夜躁狠狠久久av| 人妻人人澡人人爽人人| 最近中文字幕2019免费版| 日日爽夜夜爽网站| 国产一区二区三区综合在线观看| 欧美日韩视频高清一区二区三区二| 午夜福利一区二区在线看| av又黄又爽大尺度在线免费看| 黑人巨大精品欧美一区二区蜜桃| 最黄视频免费看| 性色av一级| 日韩av不卡免费在线播放| 中文字幕精品免费在线观看视频| 国产精品av久久久久免费| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看| xxxhd国产人妻xxx| 欧美日本中文国产一区发布| 青青草视频在线视频观看| 国产精品国产av在线观看| 中国三级夫妇交换| 少妇被粗大的猛进出69影院| 成人三级做爰电影| 日韩电影二区| 大码成人一级视频| 国产免费又黄又爽又色| 日日撸夜夜添| 久久久国产欧美日韩av| 亚洲av男天堂| 免费黄频网站在线观看国产| 亚洲一级一片aⅴ在线观看| 韩国精品一区二区三区| 日韩精品有码人妻一区| 熟妇人妻不卡中文字幕| 又大又黄又爽视频免费| 日韩制服丝袜自拍偷拍| 日韩一卡2卡3卡4卡2021年| 晚上一个人看的免费电影| 黄片无遮挡物在线观看| 热re99久久精品国产66热6| 亚洲一区二区三区欧美精品| a 毛片基地| 欧美最新免费一区二区三区| 亚洲国产最新在线播放| 90打野战视频偷拍视频| 热re99久久国产66热| 精品酒店卫生间| 亚洲 欧美一区二区三区| 高清视频免费观看一区二区| 综合色丁香网| 国产欧美亚洲国产| 高清不卡的av网站| 午夜福利免费观看在线| 女人爽到高潮嗷嗷叫在线视频| 夫妻性生交免费视频一级片| 午夜日韩欧美国产| 久久久国产精品麻豆| 国产成人免费无遮挡视频| 19禁男女啪啪无遮挡网站| 欧美日韩视频高清一区二区三区二| 宅男免费午夜| 色视频在线一区二区三区| 亚洲视频免费观看视频| 久久久久久久久免费视频了| 亚洲精品美女久久av网站| 亚洲精品久久午夜乱码| 亚洲中文av在线| 黄色怎么调成土黄色| 欧美黄色片欧美黄色片| 久久精品久久精品一区二区三区| 成人免费观看视频高清| 午夜91福利影院| 欧美日韩综合久久久久久| av卡一久久| 亚洲av电影在线进入| 好男人视频免费观看在线| 免费高清在线观看日韩| 精品久久久精品久久久| 你懂的网址亚洲精品在线观看| 精品少妇内射三级| 99精国产麻豆久久婷婷| 超碰成人久久| 永久免费av网站大全| 嫩草影院入口| 亚洲第一青青草原| 777米奇影视久久| 国产欧美日韩一区二区三区在线| 51午夜福利影视在线观看| 9色porny在线观看| 青春草视频在线免费观看| 成人国产麻豆网| 精品少妇一区二区三区视频日本电影 | 91精品三级在线观看| av又黄又爽大尺度在线免费看| 亚洲成人一二三区av| 大话2 男鬼变身卡| 国产av国产精品国产| 成人午夜精彩视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久热爱精品视频在线9| 国产精品一区二区在线观看99| 多毛熟女@视频| www.自偷自拍.com| 国产黄色免费在线视频| 国产成人a∨麻豆精品| 国产午夜精品一二区理论片| 国产精品秋霞免费鲁丝片| 久久久精品区二区三区| 欧美xxⅹ黑人| 国产成人免费无遮挡视频| 久久天堂一区二区三区四区| 亚洲成国产人片在线观看| 另类精品久久| 老鸭窝网址在线观看| 一区二区三区乱码不卡18| 亚洲在久久综合| 成人国语在线视频| 婷婷成人精品国产| 免费女性裸体啪啪无遮挡网站| 久久97久久精品| 午夜免费男女啪啪视频观看| 99久久精品国产亚洲精品| 激情视频va一区二区三区| 在线天堂最新版资源| 国精品久久久久久国模美| 赤兔流量卡办理| 操出白浆在线播放| 一级黄片播放器| 你懂的网址亚洲精品在线观看| 天堂8中文在线网| 欧美97在线视频| 99精品久久久久人妻精品| 一级黄片播放器| 一级片免费观看大全| 日韩欧美精品免费久久| 久久99一区二区三区| 国产黄色视频一区二区在线观看| 国产一区二区三区av在线| 亚洲专区中文字幕在线 | 久久久精品免费免费高清| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 一区二区三区乱码不卡18| 国产成人精品久久二区二区91 | 天堂中文最新版在线下载| a 毛片基地| 久久性视频一级片| av线在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 新久久久久国产一级毛片| 久久久久视频综合| 国产不卡av网站在线观看| 看免费av毛片| 一级毛片黄色毛片免费观看视频| 国产黄色视频一区二区在线观看| 亚洲精品日本国产第一区| 欧美黑人欧美精品刺激| 王馨瑶露胸无遮挡在线观看| 国产女主播在线喷水免费视频网站| 亚洲欧美中文字幕日韩二区| 中文欧美无线码| 嫩草影视91久久| 国产精品人妻久久久影院| 国产午夜精品一二区理论片| 国产99久久九九免费精品| 最近中文字幕2019免费版| 曰老女人黄片| 国产精品国产三级专区第一集| 999精品在线视频| 日韩不卡一区二区三区视频在线| 在线观看一区二区三区激情| 国产精品一区二区在线不卡| 性高湖久久久久久久久免费观看| 在现免费观看毛片| 成人手机av| 熟妇人妻不卡中文字幕| 2018国产大陆天天弄谢| 一本一本久久a久久精品综合妖精| 色视频在线一区二区三区| 色网站视频免费| 老司机在亚洲福利影院| 国产精品一国产av| 午夜福利视频在线观看免费| 搡老乐熟女国产| 国产成人精品福利久久| 我要看黄色一级片免费的| 伊人亚洲综合成人网| 亚洲精品在线美女| 亚洲美女视频黄频| 亚洲精品一区蜜桃| 一级毛片电影观看| 19禁男女啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 伊人亚洲综合成人网| 亚洲精品在线美女| 午夜老司机福利片| 秋霞伦理黄片| 中文字幕av电影在线播放| 亚洲第一区二区三区不卡| 韩国av在线不卡| 亚洲欧洲国产日韩| 午夜免费观看性视频| 黄色怎么调成土黄色| 精品国产乱码久久久久久男人| 免费黄色在线免费观看| 亚洲美女视频黄频| 久久久精品免费免费高清| 国产av国产精品国产| 成年美女黄网站色视频大全免费| 久久久精品区二区三区| 色吧在线观看| 亚洲第一青青草原| 美国免费a级毛片| 婷婷色麻豆天堂久久| www.自偷自拍.com| 免费av中文字幕在线| 精品酒店卫生间| a级毛片黄视频| 国产 一区精品| av在线老鸭窝| 91老司机精品| 51午夜福利影视在线观看| 在线观看www视频免费| 日日撸夜夜添| 欧美老熟妇乱子伦牲交| 高清不卡的av网站| 伊人亚洲综合成人网| 在线亚洲精品国产二区图片欧美| 自拍欧美九色日韩亚洲蝌蚪91| 在线精品无人区一区二区三| 中文字幕精品免费在线观看视频| 黄片无遮挡物在线观看| 国产精品麻豆人妻色哟哟久久| 9色porny在线观看| 中文字幕人妻丝袜制服| 亚洲成人一二三区av| 各种免费的搞黄视频| 深夜精品福利| 久久精品亚洲熟妇少妇任你| 亚洲色图 男人天堂 中文字幕| tube8黄色片| 国产一区二区 视频在线| 欧美日韩亚洲综合一区二区三区_| 国产日韩欧美视频二区| 少妇人妻久久综合中文| 宅男免费午夜| 又黄又粗又硬又大视频| 国产 一区精品| 久久国产精品男人的天堂亚洲| 久久久久久久精品精品| 满18在线观看网站| 午夜精品国产一区二区电影| 多毛熟女@视频| 国产精品一区二区在线观看99| 少妇人妻精品综合一区二区| 日本vs欧美在线观看视频| 久久亚洲国产成人精品v| 亚洲综合精品二区| 91国产中文字幕| 国产一区有黄有色的免费视频| 91老司机精品| 国产成人精品无人区| 999久久久国产精品视频| 久久久国产一区二区| 美女高潮到喷水免费观看| 男人爽女人下面视频在线观看| 日韩伦理黄色片| 80岁老熟妇乱子伦牲交| 99久久99久久久精品蜜桃| 欧美日韩成人在线一区二区| av在线老鸭窝| 国产99久久九九免费精品| 美女大奶头黄色视频| 乱人伦中国视频| 丰满迷人的少妇在线观看| 免费观看人在逋| 黄色 视频免费看| 亚洲国产中文字幕在线视频| 色94色欧美一区二区| 最近的中文字幕免费完整| av片东京热男人的天堂| 中文字幕人妻丝袜一区二区 | 这个男人来自地球电影免费观看 | 亚洲精品国产一区二区精华液| 黄色视频不卡| 免费黄色在线免费观看| 国产精品亚洲av一区麻豆 | 咕卡用的链子| 亚洲熟女精品中文字幕| 一区二区日韩欧美中文字幕| av国产久精品久网站免费入址| 久久性视频一级片| 亚洲国产看品久久| 在线观看免费高清a一片| 18在线观看网站| 国产精品久久久av美女十八| 久久 成人 亚洲| 18在线观看网站| 免费少妇av软件| 中文字幕高清在线视频| 18在线观看网站| 亚洲成人免费av在线播放| 啦啦啦中文免费视频观看日本| 久久久亚洲精品成人影院| 日日啪夜夜爽| 亚洲精品久久午夜乱码| 久久免费观看电影| 波多野结衣av一区二区av| av国产精品久久久久影院| 美女大奶头黄色视频| 丁香六月天网| 国产av国产精品国产| 婷婷色麻豆天堂久久| 天天影视国产精品| 午夜福利视频精品| av在线观看视频网站免费| 美女午夜性视频免费| 精品少妇内射三级| 国产一区二区三区av在线| 人妻 亚洲 视频| 免费在线观看视频国产中文字幕亚洲 |