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

    內(nèi)蒙古荒漠草原植物遺傳多樣性對模擬增溫處理的響應(yīng)

    2016-12-28 07:23:26李春瑞田青松杜建材王忠武
    生態(tài)學(xué)報(bào) 2016年21期
    關(guān)鍵詞:花針茅生活型細(xì)葉

    曹 路,李春瑞,田青松,杜建材,王忠武,韓 冰,*

    1 中國農(nóng)業(yè)科學(xué)院草原研究所,呼和浩特 010018 2 內(nèi)蒙古農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院,呼和浩特 010018 3 內(nèi)蒙古農(nóng)業(yè)大學(xué)生態(tài)與環(huán)境學(xué)院, 呼和浩特 010018

    內(nèi)蒙古荒漠草原植物遺傳多樣性對模擬增溫處理的響應(yīng)

    曹 路1,2,李春瑞1,2,田青松1,杜建材1,王忠武3,韓 冰1,2,*

    1 中國農(nóng)業(yè)科學(xué)院草原研究所,呼和浩特 010018 2 內(nèi)蒙古農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院,呼和浩特 010018 3 內(nèi)蒙古農(nóng)業(yè)大學(xué)生態(tài)與環(huán)境學(xué)院, 呼和浩特 010018

    為探究全球變暖對溫帶荒漠草原地上種群的遺傳影響,對已經(jīng)接受模擬增溫處理6年的短花針茅草原4種不同生活型植物,即半灌木、多年生禾草、多年生雜類草和一年生植物,應(yīng)用AFLP分子標(biāo)記方法研究了其遺傳多樣性和遺傳結(jié)構(gòu)。結(jié)果顯示,對照處理與增溫處理下的木地膚、短花針茅、細(xì)葉蔥、豬毛菜4種植物的多態(tài)位點(diǎn)百分率(PPB)分別為11.32%,11.32%;40.83%,39.91%;14.29%,13.10%;19.85%,19.12%。Nei′s基因多樣性指數(shù)(He)分別為0.0274,0.0259;0.0812,0.0899;0.0131,0.0084;0.0506,0.0456。Shannon′s 信息指數(shù)值(I)分別為0.0447,0.0430;0.1354,0.1466;0.0267,0.0182;0.0811,0.0733。分子方差分析(AMOVA)顯示4種植物的變異主要來源于實(shí)驗(yàn)處理內(nèi)部,木地膚為85.03%,短花針茅為66.35%,細(xì)葉蔥為70.00%,豬毛菜為66.52%;增溫處理間的變異分別占-2.81%,-5.47%,-3.60%,2.53%(P>0.05)。4種植物增溫處理與變異程度之間在統(tǒng)計(jì)學(xué)上并無相關(guān)性。研究表明雖然短時(shí)間的模擬增溫并不足以使4種生活型植物種群遺傳多樣性和遺傳結(jié)構(gòu)發(fā)生顯著變化,但相對于3種多年生植物,一年生植物豬毛菜更容易受到增溫影響。多年生和一年生植物對增溫具有不同的遺傳響應(yīng)。

    木地膚;短花針茅;細(xì)葉蔥;豬毛菜;遺傳多樣性;擴(kuò)增片段長度多態(tài)性(AFLP);模擬增溫

    氣候是決定地球上植被類型及物種分布最主要的因素,而植被類型和物種分布是地球氣候最鮮明的反映和標(biāo)志。目前,人們普遍認(rèn)同環(huán)境變化必然引起植物群落結(jié)構(gòu)和功能的改變,如果這種變化持續(xù)下去,影響將是深刻而持久的[1- 2]。已有的研究結(jié)果表明,增溫使發(fā)草(Deschampsiacaespitosa)、黑褐苔草(Carexalrofusca)、羊草(Leymuschinensis)的生物量有所增加;矮篙草草甸建群種和主要伴生種春季物候期提前,秋季物候期推遲;川西北高寒草甸植物發(fā)草(Deschampsiacaespitosa)的植物光合速率增加;而對植物氣孔導(dǎo)度變化沒有規(guī)律,葉片中鈉、鉀、磷的含量顯著下降[3- 4]。李英年[5]對矮嵩草(Kobresiahumilis)草甸的模擬增溫研究表明,短期增溫條件下群落的總生物量變化不明顯,增溫導(dǎo)致禾草生物量顯著增加,而雜類草生物量顯著減少。田云錄[6]利用開放式主動增溫系統(tǒng)模擬增溫使冬小麥(Triticumaestivum)地上生物量、產(chǎn)量、有效穗數(shù)和千粒重顯著提高。在賀剛[7]對煙苗的研究中,增溫增加了煙苗莖粗、葉片數(shù)、干物質(zhì)含量。而在模擬增溫條件下有關(guān)植物種群遺傳多樣性和遺傳結(jié)構(gòu)的研究還相對較少。

    在全球氣候變化的背景下,中國的氣候與環(huán)境也發(fā)生了顯著變化[8]。近50年的觀測發(fā)現(xiàn),內(nèi)蒙古自治區(qū)的氣候也有明顯變暖的趨勢。根據(jù)馬瑞芳的研究[9],內(nèi)蒙古荒漠草原氣溫呈現(xiàn)波動式升高的規(guī)律,平均溫度從20世紀(jì)50年代的8.1℃增加到90年代的9.0℃。

    為了更準(zhǔn)確地預(yù)測未來氣候變化對生態(tài)系統(tǒng)的影響,由國家自然科學(xué)基金委資助,在我國溫帶典型草原進(jìn)行了“全球變暖背景下白天和夜間增溫對溫帶草原群落結(jié)構(gòu)和生態(tài)系統(tǒng)功能的特異性影響”的增溫實(shí)驗(yàn)。本研究所在的“內(nèi)蒙古農(nóng)牧科學(xué)院四子王旗實(shí)驗(yàn)基地”是設(shè)置在該增溫樣帶上荒漠草原的一個觀測站點(diǎn),地處內(nèi)蒙古烏蘭察布市四子王旗王府一隊(duì)。在模擬全球氣候變暖的條件下,以紅外線輻射器照射提高溫度為手段[10],在野外自然條件中進(jìn)行生態(tài)系統(tǒng)控制性增溫實(shí)驗(yàn)。經(jīng)測定經(jīng)過增溫處理1年后樣地內(nèi)的0cm,7.5cm,15cm,30cm 和 50cm 土壤平均溫度與不增溫樣地相比,平均增加了1.32℃,0.92℃,0.88℃,0.80℃,0.74℃[11]。目前在該實(shí)驗(yàn)區(qū)內(nèi)已經(jīng)取得了大量的實(shí)驗(yàn)成果,模擬增溫使荒漠草原上的短花針茅(Stipabreviflora)、冷蒿(Artemisiafrigida)、無芒隱子草(Cleistogenessongorica)、木地膚(Kochiaprostrata)、阿氏旋花(Convolvulusammannii)、細(xì)葉蔥(Alliumtenuissimum)春季返青提前,總物候持續(xù)天數(shù)延長[11];但未能顯著促進(jìn)土壤呼吸[12];沒有明顯改變植物群落的組成[13];模擬增溫和氮素添加降低荒漠草原多年生植物氮素回收效率[14]。本研究選擇了4種生活型植物,分別屬于灌木、半灌木(Shrubs、halfshrubs)的木地膚、多年生禾草(Perennialgrasses)的短花針茅、多年生雜類草(Perennialforbs)的細(xì)葉蔥以及一、二年生植物(Annualsandbiennials)中的豬毛菜(Salsolacollina)。運(yùn)用擴(kuò)增片段長度多態(tài)性(amplified fragment length polymorphism,AFLP)方法,研究溫度增加對荒漠草原4種生活型植物種群結(jié)構(gòu)和遺傳多樣性的影響,揭示從2006年5月3日開始增溫至2012年7月14日,上述4種植物種群遺傳多樣性對模擬增溫的響應(yīng)趨勢,為預(yù)測野生牧草的產(chǎn)量、品質(zhì)變化及維持生態(tài)系統(tǒng)穩(wěn)定性,探究荒漠草原不同生活型的野生牧草對環(huán)境變化的適應(yīng)潛能提供實(shí)驗(yàn)依據(jù)。

    1 材料與方法

    1.1 樣地概況及增溫設(shè)置

    實(shí)驗(yàn)樣地設(shè)立在內(nèi)蒙古四子王旗境內(nèi)的“內(nèi)蒙古農(nóng)牧科學(xué)院四子王旗實(shí)驗(yàn)基地”。地理位置為東經(jīng)111°53′41 7",北緯41°46′43 6",海拔高度1456m,年均降水量約為248mm,年平均氣溫2—5℃,土壤為淡栗鈣土。該區(qū)是處于內(nèi)蒙古中部、陰山山脈以北烏蘭察布高原地區(qū)的荒漠草原,是內(nèi)蒙古草原的重要組成部分,也是草原區(qū)向荒漠區(qū)過渡的旱生性最強(qiáng)的草原生態(tài)系統(tǒng)。荒漠草原以其獨(dú)特的種類組成、結(jié)構(gòu)和功能,以及地理分布規(guī)律和獨(dú)特的自然景觀,成為亞洲中部的一種特殊的、十分重要的生態(tài)系統(tǒng)類型,具有強(qiáng)烈的大陸性特點(diǎn)?;哪菰遣菰凶詈瞪念愋?建群種由強(qiáng)旱生叢生禾草組成,經(jīng)常生有大量的強(qiáng)旱生的小半灌木、并在群落中形成穩(wěn)定的優(yōu)勢層片。試驗(yàn)區(qū)屬于以短花針茅為建群種的荒漠草原的地帶性植被,草地植物群落類型為短花針茅+冷篙+無芒隱子草,植被草層低矮,平均高度為8cm,且植被較稀疏,蓋度為12%—18%,種類組成較貧乏,植物群落主要由20多種植物組成,建群種為短花針茅,優(yōu)勢種為冷篙、無芒隱子草。研究區(qū)域內(nèi)植物類型劃分見表1。

    該樣地在野外自然條件下進(jìn)行生態(tài)系統(tǒng)控制增溫實(shí)驗(yàn),設(shè)置增溫處理小區(qū)和對照處理小區(qū),每處理6次重復(fù),共12個小區(qū)。樣地內(nèi)增溫處理小區(qū)和對照小區(qū)為隨機(jī)分布,每個小區(qū)面積為3m×4m=12m2;在6個增溫小區(qū)內(nèi)各安裝一個紅外線輻射器(Infrared Radiator)進(jìn)行加熱,該裝置通過懸掛于樣地中央上方、可以散發(fā)紅外線輻射的燈管來實(shí)現(xiàn)模擬氣候變暖。紅外線輻射器由美國伯利恒Kalglo電子公司(Kalglo Electronics Bethlehem, PA,USA)生產(chǎn),在對照小區(qū)同樣安裝用于懸掛燈管的裝置,但不安裝燈管,用以消除增溫裝置對樣地內(nèi)植物生長環(huán)境的影響。從2006年5月3日開始加熱。樣地示意圖見圖1。

    表1 樣地內(nèi)植物生活型功能群組成

    圖1 樣地示意圖Fig.1 Schematic plots 灰色陰影小區(qū)(W1—W6)為增溫處理(Warming)小區(qū),白色小區(qū)(C1—C6)為空白對照(Control)小區(qū)

    1.2 樣品采集與DNA提取

    2012年7月,在植物生長旺盛期,分別采集增溫小區(qū)和對照小區(qū)的木地膚、短花針茅、細(xì)葉蔥、豬毛菜(因W6、C6小區(qū)中木地膚和豬毛菜數(shù)量極少,數(shù)據(jù)分析過程中予以剔除)。在每個小區(qū)每種植物分單株采樣,選擇健康無病斑的植物將嫩葉剪下,用無菌水沖洗干凈后置放于2ml凍存管中,立即投放在液氮罐中備用。具體采樣情況見表2。

    采用CTAB法提取基因組DNA[15]。用瓊脂糖凝膠電泳檢測DNA完整性,用紫外分光光度計(jì)檢測DNA純度和濃度。條帶單一明亮,OD(A260/280)值在1.8—2.0之間,表明所提取的DNA純度高,可用于AFLP的操作。將DNA濃度稀釋至50ng/μL,于-20℃保存?zhèn)溆谩?/p>

    1.3 PCR擴(kuò)增與數(shù)據(jù)統(tǒng)計(jì)

    AFLP擴(kuò)增參照Vos[16]等的方法。EcoR I和Mse I酶切后連接接頭,進(jìn)行預(yù)擴(kuò)增,反應(yīng)程序:94℃預(yù)變性3min; (94℃, 40s; 60℃, 40s; 72℃, 60s) ×35循環(huán);72℃,10min。每種植物從32對選擇性擴(kuò)增引物中篩選出5對信號強(qiáng)度一致、條帶分布均勻、清晰并具有多態(tài)性的引物進(jìn)行后續(xù)實(shí)驗(yàn)。擴(kuò)增產(chǎn)物用8%的聚丙烯酰胺凝膠電泳檢測。銀染后室溫干燥過夜,掃描圖像。AFLP條帶大小從100—1500 bp,按照電泳結(jié)果中,同一遷移率上條帶的有無進(jìn)行統(tǒng)計(jì),有條帶的記為“1”,無條帶的記為“0”,獲得01矩陣。

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

    根據(jù)AFLP獲得的01矩陣,分別利用POPGEN 1.31[17]軟件計(jì)算4種生活型植物處理間的位點(diǎn)數(shù),多態(tài)位點(diǎn)數(shù),多態(tài)位點(diǎn)百分率(PPB),Nei′s基因多樣性指數(shù)(He),Shannon′s信息指數(shù)(I),遺傳一致度和遺傳距離。利用NTSYS 2.1[18]軟件進(jìn)行非加權(quán)算術(shù)平均法(UPGMA)聚類。并利用Arlequin3.11[19]軟件進(jìn)行分子方差分析及等級剖分分析。

    表2 4種生活型植物的采樣情況

    2 結(jié)果與分析

    2.1 遺傳多樣性分析

    經(jīng)過POPGEN 1.32軟件的計(jì)算,木地膚5對引物擴(kuò)增出106個位點(diǎn),其中對照小區(qū)多態(tài)性位點(diǎn)12個,多態(tài)位點(diǎn)百分率11.32%;增溫小區(qū)多態(tài)性位點(diǎn)12個,多態(tài)位點(diǎn)百分率11.32%。短花針茅5對引物擴(kuò)增出218個位點(diǎn),其中對照小區(qū)多態(tài)性位點(diǎn)89個,多態(tài)位點(diǎn)比例40.83%;增溫小區(qū)多態(tài)性位點(diǎn)87個,多態(tài)位點(diǎn)百分率39.91%。細(xì)葉蔥5對引物擴(kuò)增出84個位點(diǎn),其中對照小區(qū)多態(tài)性位點(diǎn)12個,多態(tài)位點(diǎn)百分率14.29%;增溫小區(qū)多態(tài)性位點(diǎn)11個,多態(tài)位點(diǎn)百分率13.10%。豬毛菜5對引物擴(kuò)增出136個位點(diǎn),其中對照小區(qū)多態(tài)性位點(diǎn)27個,多態(tài)位點(diǎn)百分率19.85%;增溫小區(qū)多態(tài)性位點(diǎn)26個,多態(tài)位點(diǎn)百分率19.12%。Nei′s基因多樣性指數(shù)與Shannon′s 信息指數(shù)普遍較低,4種植物對照與增溫處理下的Nei′s基因多樣性指數(shù)(He)分別為0.0274,0.0259(木地膚);0.0812,0.0899(短花針茅);0.0131,0.0084(細(xì)葉蔥);0.0506,0.04569(豬毛菜)。Shannon′s 信息指數(shù)值(I)分別為0.0447,0.0430(木地膚);0.1354,0.14669(短花針茅);0.0267,0.0182(細(xì)葉蔥);0.0811,0.07339(豬毛菜)。4種植物在增溫處理和對照間的遺傳差異較小(表3)。增溫處理并未對其遺傳多樣性產(chǎn)生明顯影響。

    2.2 遺傳分化

    遺傳一致度和遺傳距離代表了不同處理間的植物的親緣關(guān)系。遺傳一致度越大、遺傳距離越小,則親緣關(guān)系較近。由表4可見木地膚、短花針茅、細(xì)葉蔥和豬毛菜增溫與對照之間的的遺傳一致度分別為0.9999,0.9963,0.9994,0.9929;遺傳距離分別為0.0001,0.0037,0.0006,0.0071。其中豬毛菜的遺傳一致度最小為0.9929,遺傳距離最遠(yuǎn)為0.0071。

    根據(jù)遺傳一致度,利用NTSYS 2.10軟件進(jìn)行非加權(quán)組平均法(unweighted pair-group method with arithmetic means, UPGMA)聚類。圖2為4種植物的不同小區(qū)間的聚類圖,其中木地膚增溫處理小區(qū)W4、W1、W3分別依次和對照小區(qū)的C5、C4、C3形成3個一級聚類,依次形成的二、三級及高級聚類都沒有將增溫和對照小區(qū)聚為大類;短花針茅雖然C1和C2、W1和W2、C3和C6這6個小區(qū)分別形成3個相同處理的一級聚類,但是在更高級的聚類中沒有相同規(guī)律;細(xì)葉蔥的C1和W4形成一級聚類,W2、W3和C3形成一級聚類,然后依次與C4、W6和C5形成二、三、四級聚類;豬毛菜W3和W5,C4和W4,W1和W2分別形成一級聚類,進(jìn)一步的聚類也沒有規(guī)律??梢?種植物并沒有按照增溫與對照處理的聚為不同的兩支,聚類皆沒有規(guī)律可循。

    表3 4種植物的遺傳多樣性分析

    表 4 4種植物處理間的遺傳一致度和遺傳距離

    圖2 4種生活型植物的UPGMA聚類圖Fig.2 Dendrogram by UPGMA analysis of four life forms plantsA: 木地膚10個小區(qū)的UPGMA聚類圖;B: 短花針茅12個小區(qū)的UPGMA聚類圖;C: 細(xì)葉蔥12個小區(qū)的UPGMA聚類圖;D: 豬毛菜10個小區(qū)的UPGMA聚類圖

    2.3 分子方差分析

    對4種植物進(jìn)行AMOVA分析,結(jié)果(表5)顯示4種植物的變異主要來源于種群內(nèi)。其中,木地膚種群內(nèi)變異最多85.03%,細(xì)葉蔥70.00%,豬毛菜次之66.52%,短花針茅最少66.35%。而4種植物來源于處理間的變異都較少分別為-2.81%(木地膚)、-5.47%(短花針茅)、-3.60%(細(xì)葉蔥)、2.53%(豬毛菜)。其中木地膚、短花針茅、細(xì)葉蔥的處理間的變異百分率都為負(fù)值,說明增溫處理與其變異并沒有相關(guān)性。而豬毛菜來源于處理間的變異占到了總變異的2.53%,雖然P值仍然大于0.05,但表明在模擬增溫處理6年之后發(fā)生了一定程度上的遺傳結(jié)構(gòu)的變化。豬毛菜作為一年生植物,表現(xiàn)出與其他生活型植物不同的響應(yīng)。

    表5 4種植物的AMOVA分析

    通過等級剖分(nested analysis of molecular variance)分析4種植物種群內(nèi)、處理內(nèi)種群間和處理間的遺傳分化系數(shù)(表6)。處理間的遺傳分化都較低,其中木地膚、短花針茅、細(xì)葉蔥都為負(fù)值(-0.02807、-0.05468、-0.03602),只有豬毛菜數(shù)為正值(0.02527)。說明木地膚、短花針茅、細(xì)葉蔥3種植物的遺傳分化與增溫處理之間并無相關(guān)性。而豬毛菜在增溫處理之下發(fā)生了一定程度遺傳變異。

    表6 4種植物的遺傳分化系數(shù)

    3 討論與結(jié)論

    氣候預(yù)測顯示,未來我國氣溫仍將升高[3]。在野外自然條件下的生態(tài)系統(tǒng)增溫實(shí)驗(yàn)成為了研究全球變暖與陸地生態(tài)系統(tǒng)的關(guān)鍵[20]。本研究所使用的紅外線輻射器(Infrared radiator)可以同時(shí)影響顯熱、潛熱和土壤熱通量[21],因而有效的模擬了全球變暖的情況。

    在本研究中4種生活型植物木地膚、短花針茅、細(xì)葉蔥、豬毛菜的多態(tài)位點(diǎn)百分率分別為12.26%、41.28、14.29、20.59,Nei′s基因多樣性指數(shù)分別為0.0268、0.0863、0.0109、0.0522和Shannon′s 信息指數(shù)分別為0.0445、0.1428、0.0229、0.0834,都較低,遺傳多樣性并不豐富。對4種生活型植物的聚類圖分析也顯示增溫處理組與對照組并未聚為兩個類群。首先,這與樣地設(shè)置有關(guān),在以往有關(guān)遺傳多樣性的研究中大都是以不同品種或是地理間隔較遠(yuǎn)的植物為研究目標(biāo)[22- 24],如在對大針茅(Stipagrandis)[25]和克氏針茅(Stipakrylovii)[26]遺傳多樣性的研究中,研究者在同一地點(diǎn)不同樣本之間的采樣距離也會間隔50m以上,因?yàn)橐酝际茄芯孔匀环N群遺傳差異,由于研究目的不同,地理種群間距離較遠(yuǎn),種群間有較高的遺傳多樣性(克氏針茅Nei′s基因多樣性指數(shù)0. 2989,Shannon′s 信息指數(shù)0. 4014;大針茅Nei′s基因多樣性指數(shù)0.3486,Shannon′s 信息指數(shù)0.2314)而在本研究中12個小區(qū)間相隔只有3m,由于間隔近,不同小區(qū)間植物的遺傳背景較一致,有的不同小區(qū)內(nèi)的植株甚至可能來源于同一個母株的種子。這是種群間遺傳多樣性較低的直接原因。其次,由于增溫處理的時(shí)間有限,6a的時(shí)間很有可能并不足以使植物種群遺傳發(fā)生明顯分化。在模擬增溫對青藏高原草甸生態(tài)系統(tǒng)影響的研究中發(fā)現(xiàn)[27],在增溫實(shí)驗(yàn)的前兩年出現(xiàn)植物初級生產(chǎn)量大幅增加,生物多樣性增加的現(xiàn)象,而在實(shí)驗(yàn)的第3年之后初級生產(chǎn)量開始下降,并且生物多樣性也在降低。這說明在植物響應(yīng)模擬增溫的過程,其反應(yīng)隨處理時(shí)間的長短會產(chǎn)生不同的結(jié)果,推測短花針茅草原各種生活型植物在遺傳上對增溫的響應(yīng)還沒有出現(xiàn)。另外多年生植物在樣地內(nèi)并沒有適宜的條件允許實(shí)生苗建植,即使增溫處理使種子產(chǎn)生了某些變化,這些種子也未必能在小區(qū)內(nèi)形成植株,繼而無法被采集到相關(guān)遺傳改變的信息;而且多年生植物對這種溫和的增溫處理在新芽和新生葉片生長過程中可能未受到影響。

    分子方差分析是1992年由Excoffier[19]發(fā)展出來的一種進(jìn)行遺傳變異等級剖分的方法,并提出了與F-statistics類似的PHI-statistics等方法來有效地度量亞群體的分化。在本研究中,4種生活型植物當(dāng)中的多年生植物處理間變異占總變異的百分率(Total variance)以及遺傳分化系數(shù)(PHI-statistics)都為負(fù)值。這種現(xiàn)象在張麗[28]對中華白蛉(Phlebotomuschinensis)的研究中也曾經(jīng)出現(xiàn)過,負(fù)值說明處理間不存在遺傳分化。而一年生植物豬毛菜則表現(xiàn)出了不同,其處理間的遺傳分化系數(shù)為0.0253,處理間的變異百分率也達(dá)到了2.53%,雖然P值仍然大于0.05,但表現(xiàn)出一定程度的遺傳分化趨勢。這4種生活型植物都是荒漠草原的典型植物,從形態(tài)結(jié)構(gòu)上來看,他們都具有典型的旱生結(jié)構(gòu),可以有效的反射陽光,控制水分的蒸發(fā)[29- 32]。以往的研究結(jié)果顯示一年生植物對增溫處理有較敏感的響應(yīng),常兆豐[33]對民勤荒漠區(qū)不同生活型植物響應(yīng)氣候變暖的研究中發(fā)現(xiàn),草本植物的春季物候和秋季物候響應(yīng)氣候變暖的敏感程度都大于灌木。在楊健[34]的研究中發(fā)現(xiàn),增溫處理對多年生植物加拿大一枝黃花(Hemerocalliscitrina)地下部分沒有影響,但趨于降低一年生植物鬼針草(Bidenspilosa)對根系生物量的分配。本研究中作為一年生植物的豬毛菜很可能是因?yàn)榫邆淞烁痰纳钪芷?才能夠在相同的時(shí)間內(nèi)較早的對模擬增溫作出了遺傳層面上的應(yīng)答。推測每次含有一些變異的種子飄落后都會在小區(qū)內(nèi)參與下一年群落的形成,因而增溫處理對其影響在每個世代(每年)中都有一些積累,雖然這種積累非常的微弱。這種變異在一定程度上反映了增溫對植物影響的趨勢。

    而王迪[35]等認(rèn)為,當(dāng)植物遇到環(huán)境壓力時(shí),即使DNA序列水平上沒有發(fā)生變化,也會發(fā)生表觀遺傳學(xué)上的變化,這種變化一部分是可以遺傳的,一部分會因?yàn)閴毫Φ娜コ谕淮?xì)胞中發(fā)生回復(fù)。近年來表觀遺傳學(xué)已經(jīng)越來越受到關(guān)注[36]。和傳統(tǒng)意義上的遺傳學(xué)相比,所謂的表觀遺傳學(xué),就是不涉及DNA序列改變,然而最終可能影響表型變化的可遺傳的現(xiàn)象[37]。在作者其他的研究中發(fā)現(xiàn),短花針茅雖然在增溫處理與未增溫處理間DNA水平?jīng)]有發(fā)生改變,但是在轉(zhuǎn)錄組水平上發(fā)現(xiàn)了一些基因的差異表達(dá),表明增溫影響了基因的表達(dá),這其中不能否定有表觀遺傳學(xué)水平的變化,模擬增溫使表觀遺傳的修飾程度會發(fā)生改變。自然選擇是塑造和維持植物自然遺傳多樣性與種群遺傳結(jié)構(gòu)的主要動力[38]。隨著全球氣候變暖的趨勢,植物的種群[39]、群落[40- 41]、生態(tài)系統(tǒng)[42-45]乃至整個生物圈[46]必將發(fā)生巨大改變。而本研究發(fā)現(xiàn)不同的生活型植物面對這一變化產(chǎn)生不同的結(jié)果,一年生植物會容易受到增溫影響。隨著時(shí)間的推移,研究的深入,在之后的研究中將會加入對增溫去除組的關(guān)注,并進(jìn)行表觀遺傳學(xué)的研究,將對模擬增溫下不同生活型植物的適應(yīng)機(jī)制會有更加透徹的了解。

    [1] Dwire K A, Kauffman J B, Brookshire E N J, Baham J E. Plant biomass and species composition along an environmental gradient in montane riparian meadows. Oecologia, 2004, 139(2): 309- 317.

    [2] Alemayehu F R, Frenck G, van der Linden L, Mikkelsen T N, J?rgensen R B. Can barley (HordeumvulgareL. s.l.) adapt to fast climate changes? A controlled selection experiment. Genetic Resources and Crop Evolution, 2014, 61(1): 151- 161.

    [3] 李娜, 王根緒, 楊燕, 高永恒, 柳林安, 劉光生. 短期增溫對青藏高原高寒草甸植物群落結(jié)構(gòu)和生物量的影響. 生態(tài)學(xué)報(bào), 2011, 31(4): 895- 905.

    [4] 黃文華, 王樹彥, 韓冰, 焦志軍, 韓國棟. 草地生態(tài)系統(tǒng)對模擬大氣增溫的響應(yīng). 草業(yè)科學(xué), 2014, 31(11): 2069- 2076.

    [5] 李英年, 趙亮, 趙新全, 周華坤. 5年模擬增溫后矮嵩草草甸群落結(jié)構(gòu)及生產(chǎn)量的變化. 草地學(xué)報(bào), 2004, 12(3): 236- 239.

    [6] 田云錄, 鄭建初, 張彬, 陳金, 董文軍, 楊飛, 張衛(wèi)建. 麥田開放式晝夜不同增溫系統(tǒng)的設(shè)計(jì)及增溫效果. 中國農(nóng)業(yè)科學(xué), 2010, 43(18): 3724- 3731.

    [7] 賀剛, 周輝, 廖衡斌, 劉喜慶. 烤煙漂浮育苗中增溫對煙苗的影響. 作物研究, 2014, 28(7): 798- 800.

    [8] 丁一匯, 任國玉, 石廣玉, 宮鵬, 鄭循華, 翟盤茂, 張德二, 趙宗慈, 王紹武, 王會軍, 羅勇, 陳德亮, 高學(xué)杰, 戴曉蘇. 氣候變化國家評估報(bào)告(Ⅰ): 中國氣候變化的歷史和未來趨勢. 氣候變化研究進(jìn)展, 2006, 2(1): 3- 8.

    [9] 馬瑞芳. 內(nèi)蒙古草原區(qū)近50年氣候變化及其對草地生產(chǎn)力的影響[D]. 北京: 中國農(nóng)業(yè)科學(xué)院研究生院, 2007.

    [10] Shaver G R, Johnson L C, Cades D H, Murray G, Laundre J A, Rastetter E B, Nadelhoffer K J, Giblin A E. Biomass and CO2flux in wet sedge tundra: responses to nutrients, temperature, and light. Ecological Monographs, 1998, 68(1): 75- 97.

    [11] 珊丹, 韓國棟, 趙萌莉, 王珍, 韓雄, 高福光. 控制性增溫和施氮對荒漠草原土壤呼吸的影響. 干旱區(qū)資源與環(huán)境, 2009, 23(9): 106- 112.

    [12] 張宇, 紅梅. 內(nèi)蒙古荒漠草原土壤呼吸對模擬增溫和氮素添加的響應(yīng). 草地學(xué)報(bào), 2014, 22(6): 1227- 1231.

    [13] 王晨晨, 王珍, 張新杰, 李寅龍, 劉軍利, 韓國棟. 增溫對荒漠草原植物群落組成及物種多樣性的影響. 生態(tài)環(huán)境學(xué)報(bào), 2014, 23(1): 43- 49.

    [14] 李元恒, 韓國棟, 王珍, 王正文, 趙萌莉, 王薩仁娜. 增溫和氮素添加降低荒漠草原多年生植物氮素回收效率. 生態(tài)學(xué)報(bào), 2015, 35(18):. (未出版刊物)

    [15] 李榮華, 夏巖石, 劉順枝, 孫莉麗, 郭培國, 繆紳裕, 陳健輝. 改進(jìn)的CTAB提取植物DNA方法. 實(shí)驗(yàn)室研究與探索, 2009, 28(9): 14- 16.

    [16] Vos P, Hogers R, Bleeker M, Reijans, M, van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M. AFLP: a new technique for DNA fingerprinting. Nucleic Acids Research, 1995, 23(21): 4407- 4414.

    [17] Yeh F C, Yang R C, Boyle T B J. POPGENE version 1. 31. Microsoft Windows-based freeware for population genetic analysis. Edmonton: University of Alberta, 1997.

    [18] Rohlf F J. NTSYS-PC: Numerical Taxonomy and Multivariate Analysis System, Version 2. 1. New York: Exeter Software, Setauket, 2000.

    [19] Excoffier L, Smouse P E, Quattro J M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics, 1992, 131(2): 479- 491.

    [20] Shen K P, Harte J. Ecosystem climate manipulations // Methods in Ecosystem Science. New York: Springer-Verlag, 2000: 353- 369.

    [21] 牛書麗, 韓興國, 馬克平, 萬師強(qiáng). 全球變暖與陸地生態(tài)系統(tǒng)研究中的野外增溫裝置. 植物生態(tài)學(xué)報(bào), 2007, 31(2): 262- 271.

    [22] 李娜, 楊鷹, 萬潔, 韋雷飛, 鐘金城, 武志娟, 姜雪鷗, 張利亞. 攀西黑山羊群體的AFLP遺傳多樣性研究. 家畜生態(tài)學(xué)報(bào), 2013, 34(2): 22- 26.

    [23] 趙冰, 徐曼, 司國臣, 李厚華, 張延龍. 秦嶺秀雅杜鵑野生種群遺傳多樣性和遺傳分化的AFLP分析. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(11): 2983- 2990.

    [24] 張慶, 牛建明, 董建軍. 內(nèi)蒙古地區(qū)短花針茅(Stipabreviflora)種群遺傳多樣性. 生態(tài)學(xué)報(bào), 2008, 28(7): 3447- 3455.

    [25] 珊丹, 趙萌莉, 韓冰, 韓國棟. 不同放牧壓力下大針茅種群的遺傳多樣性. 生態(tài)學(xué)報(bào), 2006, 26(10): 3175- 3183.

    [26] 王金龍, 高玉葆, 趙念席, 劉景玲, 林楓, 宋濤. 內(nèi)蒙古中東部草原區(qū)克氏針茅形態(tài)特征和RAPD遺傳分化的相關(guān)性研究. 植物研究, 2006, 26(6): 708- 714.

    [27] 趙新全, 陳世龍, 曹廣民, 師生波, 徐世曉, 王啟蘭. 研究方向: 青藏高原高寒草甸生態(tài)系統(tǒng)與全球氣候變化的相互作用機(jī)理研究. 科技和產(chǎn)業(yè), 2003, 3(8): 51- 59.

    [28] 張麗. 我國中華白蛉的鑒別和分子群體遺傳結(jié)構(gòu)研究[D]. 上海: 第二軍醫(yī)大學(xué), 2010.

    [29] 康薩如拉, 牛建明, 張慶, 陳麗萍. 短花針茅葉片解剖結(jié)構(gòu)及與氣候因子的關(guān)系. 草業(yè)學(xué)報(bào), 2013, 22(1): 77- 86.

    [30] 高松, 蘇培璽, 嚴(yán)巧娣, 丁松爽, 張嶺梅. C4荒漠植物豬毛菜與木本豬毛菜的葉片解剖結(jié)構(gòu)及光合生理特征. 植物生態(tài)學(xué)報(bào), 2009, 33(2): 347- 354.

    [31] 劉濤, 李柱, 安沙舟, 許幗英. 干旱脅迫對木地膚幼苗生理生化特性的影響. 干旱區(qū)研究, 2008, 25(2): 231- 235.

    [32] 栗利元, 張未芳, 鄭聯(lián)壽. 細(xì)葉韭生活習(xí)性的研究. 黑龍江農(nóng)業(yè)科學(xué), 2011, (2): 63- 65.

    [33] 常兆豐, 王強(qiáng)強(qiáng), 韓福貴, 仲生年. 民勤荒漠區(qū)不同生活型植物物候響應(yīng)氣候變暖的差異. 生態(tài)學(xué)雜志, 2011, 30(9): 1921- 1929.

    [34] 楊健. 增溫和氮沉降對加拿大一枝黃花入侵性的影響[D]. 成都: 成都理工大學(xué), 2014.

    [35] 王迪, 傅彬英, 張立軍. 植物表觀遺傳變化與環(huán)境壓力研究進(jìn)展. 分子植物育種, 2008, 6(3): 569- 573.

    [36] Bender J. Plant epigenetics. Current Biology, 2002, 12(12): R412-R414.

    [37] Bird A. Perceptions of epigenetics. Nature, 2007, 447(7143): 396- 398.

    [38] Linhart Y B, Grant M C. Evolutionary significance of local genetic differentiation in plants. Annual Review of Ecology and Systematics, 1996, 27: 237- 277.

    [39] Klanderud K. Climate change effects on species interactions in an alpine plant community. Journal of Ecology, 2005, 93(1): 127- 137.

    [40] Klein J A, Harte J, Zhao X Q. Experimental warming causes large and rapid species loss, dampened by simulated grazing, on the Tibetan Plateau. Ecology Letters, 2004, 7(12): 1170- 1179.

    [41] Harte J, Shaw R. Shifting dominance within a montane vegetation community: results of a climate-warming experiment. Science, 1995, 267(5199): 876- 880.

    [42] Tilman D, Knops J, Wedin D, Reich P, Ritchie M, Siemann E. The influence of functional diversity and composition on ecosystem processes. Science, 1997, 277(5330): 1300- 1302.

    [43] Zhou G S, Wang Y H, Zhang X S. Advance in the responses of vegetation/ecosystems in China to global change. Bulletin of the Chinese Academy of Sciences, 1999, 13(3): 158- 165.

    [44] Aerts R, Cornelissen J H C, Dorrepaal E. Plant performance in a warmer world: general responses of plants from cold, northern biomes and the importance of winter and spring events. Plant Ecology, 2006, 182(1/2): 65- 77.

    [45] Morgan J A, LeCain D R, Pendall E, Blumenthal D M, Kimball B A, Carrillo Y, Williams D G, Heisler-White J, Dijkstra F A, West M. C4 grasses prosper as carbon dioxide eliminates desiccation in warmed semi-arid grassland. Nature, 2011, 476(7359): 202- 205.

    [46] Yin H J, Lai T, Cheng X Y, Jiang X M, Liu Q. Warming effect on growth and physiology of seedlings ofBetulaalbo-sinensisandAbiesfaxonianaunder two contrasting light conditions in subalpine coniferous forests of western Sichuan, China. Frontiers of Forestry in China, 2009, 4(4): 432- 442.

    The response of genetic diversity in desert steppe plants to simulated warming in Inner Mongolia, China

    CAO Lu1,2,LI Chunrui1,2,TIAN Qingsong1,DU Jiancai1,WANG Zhongwu3,HAN Bing1,2,*

    1InstituteofGrasslandResearch,ChineseAcademyofAgriculturalSciences,Hohhot010018,China2CollegeofLifeScience,InnerMongoliaAgriculturalUniversity,Hohhot010018,China3CollegeofEcologyandEnvironmentalScience,InnerMongoliaAgriculturalUniversity,Hohhot010018,China

    Climate is the most important factor that determines vegetation types and the distribution of species, and, accordingly, these features are two of the most prominent indicators of climate change. Currently, scientists generally agree that climate change will inevitably lead to changes in plant community structure and function, and if this change continues, the effects will be profound and enduring. Existing research shows that simulated warming causes an increase in biomass inDeschampsiacaespitosa,Carexalrofusca, andLeymuschinensis; as climates change, these constructive species and their main companion species within aKobresiahumilismeadow have experienced an advancement in their spring phenology and a delay in their autumn phenophase. In addition, the photosynthetic rate ofDeschampsiacaespitosain a Northwest Sichuan alpine meadow has increased. Although the change in stomatal conductance was irregular, a significant decrease has also been observed in the sodium, potassium, and phosphorus content of leaves. However, relatively little research has been conducted on the effects of simulated warming on the genetic structure and diversity of plant populations. Against the background of global climate change, the temperatures of the Inner Mongolian desert steppe have become unevenly elevated, with average annual temperatures increasing from 8.1℃ in the 1950s to 9.0℃ in the 1990s. The present study explores the effects of global warming on the genetics of wild forage plant populations with different life forms in desert steppe habitat, in an effort to elucidate their potential to adapt to environmental change. The study site was located in the desert steppe in Siziwangqi territory of Inner Mongolia. Suspension infrared radiators were used to create a controlled warming experiment under otherwise natural field conditions. Warming began on May 3, 2006. The average soil temperature at depths of 0, 7.5, 15, 30, and 50 cm increased by 1.32, 0.92, 0.88, 0.80, and 0.74℃, respectively, after warming for 1 year compared with the average in plots not exposed to warming. For this study, changes in genetic diversity and structure were analyzed in four plant populations: (1) small half shrubs, represented byKochiaprostrata; (2) perennial grasses, represented byStipabreviflora; (3) perennial forbs, represented byAlliumtenuissimum; and (4) annuals and biennials, represented bySalsolacollina). The study was conducted under simulated warming pressure and genetic analysis was performed using amplified fragment length polymorphism. The percentages of polymorphic loci inK.prostrata,S.breviflora,A.tenuissimum, andS.collinaunder non-warming were 11.32%, 40.83%, 14.29%, and 19.85%, whereas those under simulated warming were 11.32%, 39.91%, 13.10%, and 19.12%, respectively. The genetic diversity of the four populations measured under control and simulated warming, measured by the Shannon′s information index, were as follows:K.prostrata(0.0274, 0.0259),S.breviflora(0.0812, 0.0899),Alliumtenuissimum(0.0131, 0.0084), andS.collina(0.0506, 0.0456). These findings exhibited the same distributional pattern as that of Nei′s genetic diversity index forK.prostrata(0.0447, 0.0430),S.breviflora(0.1354, 0.1466),Alliumtenuissimum(0.0267, 0.0182), andS.collina(0.0811, 0.0733). Cluster analysis of these four species showed that known of the species reacted significantly to the warming process, and that inter-individual clusters were not significantly different. The results of an analysis of molecular variance (AMOVA) indicated that the main source of variation among the four life form populations was within-population variation:K.prostrata(85.03%),S.breviflora(66.35%),A.tenuissimum(70.00%), andS.collina(66.52%). The among-groups variation was not significant and accounted for the following percentages of variation:K.prostrata(-2.81%),S.breviflora(-5.47%),Alliumtenuissimum(-3.60%), andS.collina(2.53%). No statistically significant correlation was found between simulated warming and genetic differentiation. This study shows that a short period of simulated warming was not sufficient to create a significant change in genetic diversity and structure for the four life form populations studied here; however, compared with the three types of perennials studied, the annual plantS.collina, is more susceptible to the effects of warming. Perennials and annuals have different genetic responses to warming. This study provides experimental evidence that can reveal the potential adaptation of plants to environmental change for different life forms of wild forage plants of the desert steppe, and will help researchers to predict forage yield and changes in forage quality.

    Kochiaprostrata;Stipabreviflora;Alliumtenuissimum;Salsolacollina; genetic diversity; AFLP; simulated warming

    內(nèi)蒙古自治區(qū)自然科學(xué)基金資助項(xiàng)目(2015MS0305);中國科學(xué)院西部之光“人才培養(yǎng)”項(xiàng)目;中國農(nóng)業(yè)科學(xué)院草原研究所科技創(chuàng)新工程資助項(xiàng)目

    2015- 04- 22;

    日期:2016- 03- 03

    10.5846/stxb201504220826

    *通訊作者Corresponding author.E-mail: hb_nmg@163.com

    曹路,李春瑞,田青松,杜建材,王忠武,韓冰.內(nèi)蒙古荒漠草原植物遺傳多樣性對模擬增溫處理的響應(yīng).生態(tài)學(xué)報(bào),2016,36(21):6909- 6918.

    Cao L,Li C R,Tian Q S,Du J C,Wang Z W,Han B.The response of genetic diversity in desert steppe plants to simulated warming in Inner Mongolia, China.Acta Ecologica Sinica,2016,36(21):6909- 6918.

    猜你喜歡
    花針茅生活型細(xì)葉
    明堂山主要藥用植物資源調(diào)查分析
    短花針茅葉片熱值及還原糖含量對載畜率的響應(yīng)
    放牧強(qiáng)度對短花針茅植物種群空間分布的影響
    不同載畜率對荒漠草原建群種短花針茅種群時(shí)空變異性的影響
    放牧強(qiáng)度對荒漠草原建群種短花針茅空間異質(zhì)性的影響
    紫蕓英
    論數(shù)學(xué)教學(xué)中有效運(yùn)用數(shù)學(xué)模型
    成才之路(2017年21期)2017-07-27 20:52:40
    寬闊水國家級自然保護(hù)區(qū)珙桐群落特征及演替趨勢分析
    廣西細(xì)葉云南松森林資源動態(tài)分析
    關(guān)于初中數(shù)學(xué)作業(yè)布置的幾點(diǎn)思考
    国产爱豆传媒在线观看 | 香蕉丝袜av| 久久精品亚洲精品国产色婷小说| 熟女少妇亚洲综合色aaa.| 日本a在线网址| 国产精品久久久久久精品电影 | 成人av一区二区三区在线看| 制服人妻中文乱码| 亚洲欧美日韩无卡精品| 亚洲精品av麻豆狂野| 99久久国产精品久久久| a在线观看视频网站| 亚洲人成网站在线播放欧美日韩| 97人妻精品一区二区三区麻豆 | 亚洲男人天堂网一区| 精品国产超薄肉色丝袜足j| 日日夜夜操网爽| 嫁个100分男人电影在线观看| 国内精品久久久久精免费| 成人三级做爰电影| 满18在线观看网站| 男女做爰动态图高潮gif福利片| 香蕉av资源在线| 婷婷六月久久综合丁香| 男女下面进入的视频免费午夜 | av有码第一页| 久久久久久亚洲精品国产蜜桃av| 亚洲av第一区精品v没综合| 在线免费观看的www视频| 亚洲国产看品久久| 久热这里只有精品99| 亚洲精品国产精品久久久不卡| 久久精品aⅴ一区二区三区四区| 久久精品人妻少妇| 看免费av毛片| 少妇被粗大的猛进出69影院| 在线观看66精品国产| 亚洲精品中文字幕在线视频| 国产一区在线观看成人免费| 亚洲精品色激情综合| 亚洲七黄色美女视频| 久久久久九九精品影院| 国产国语露脸激情在线看| 在线播放国产精品三级| 成人永久免费在线观看视频| 亚洲专区国产一区二区| 真人做人爱边吃奶动态| 成人一区二区视频在线观看| 最近最新中文字幕大全电影3 | 99在线视频只有这里精品首页| 免费在线观看日本一区| 日本五十路高清| 国语自产精品视频在线第100页| 91老司机精品| 午夜久久久在线观看| 国产伦人伦偷精品视频| 欧美人与性动交α欧美精品济南到| 免费在线观看影片大全网站| 欧美精品亚洲一区二区| 欧美乱码精品一区二区三区| 欧美成人免费av一区二区三区| 老司机福利观看| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 精品无人区乱码1区二区| 亚洲成人免费电影在线观看| 欧美日韩一级在线毛片| 亚洲中文av在线| 亚洲五月婷婷丁香| 国产国语露脸激情在线看| 免费在线观看成人毛片| 夜夜夜夜夜久久久久| 中文在线观看免费www的网站 | 麻豆国产av国片精品| 免费在线观看成人毛片| 中文字幕最新亚洲高清| 日日摸夜夜添夜夜添小说| 在线十欧美十亚洲十日本专区| 成人av一区二区三区在线看| 久久精品夜夜夜夜夜久久蜜豆 | 国产在线精品亚洲第一网站| 丁香六月欧美| 亚洲午夜精品一区,二区,三区| 精品国产亚洲在线| 成人三级黄色视频| 免费av毛片视频| 亚洲人成77777在线视频| 国产精品九九99| 国产成人精品久久二区二区91| 在线观看免费午夜福利视频| 精品日产1卡2卡| 91九色精品人成在线观看| 999久久久精品免费观看国产| 国产男靠女视频免费网站| 99精品欧美一区二区三区四区| 免费搜索国产男女视频| x7x7x7水蜜桃| 亚洲欧洲精品一区二区精品久久久| netflix在线观看网站| 黄网站色视频无遮挡免费观看| 一a级毛片在线观看| 亚洲av第一区精品v没综合| 国产精品久久视频播放| 成人国语在线视频| 精品国内亚洲2022精品成人| 在线观看www视频免费| 久久天堂一区二区三区四区| 最近最新中文字幕大全免费视频| 在线看三级毛片| 国产又爽黄色视频| 成在线人永久免费视频| 精品熟女少妇八av免费久了| 午夜久久久在线观看| 欧美国产精品va在线观看不卡| 欧美+亚洲+日韩+国产| 日韩欧美一区二区三区在线观看| 日韩视频一区二区在线观看| 丝袜人妻中文字幕| 亚洲成国产人片在线观看| 欧美一级a爱片免费观看看 | 久久久国产精品麻豆| 天堂影院成人在线观看| 欧美最黄视频在线播放免费| 国产精品98久久久久久宅男小说| 一区二区三区高清视频在线| 日本 av在线| 国产成年人精品一区二区| 制服诱惑二区| 日韩欧美国产在线观看| 欧美大码av| 欧美日本视频| 国产一级毛片七仙女欲春2 | 亚洲男人天堂网一区| 波多野结衣高清作品| 国产精品一区二区三区四区久久 | 亚洲国产精品久久男人天堂| 亚洲中文日韩欧美视频| 亚洲成人免费电影在线观看| 男男h啪啪无遮挡| 亚洲va日本ⅴa欧美va伊人久久| 日本三级黄在线观看| 欧美大码av| av片东京热男人的天堂| 男女之事视频高清在线观看| 在线天堂中文资源库| 久久国产精品男人的天堂亚洲| 自线自在国产av| 最好的美女福利视频网| 国内揄拍国产精品人妻在线 | 少妇熟女aⅴ在线视频| 男人操女人黄网站| 午夜久久久久精精品| 精品一区二区三区四区五区乱码| 亚洲欧美精品综合久久99| 夜夜躁狠狠躁天天躁| 黄色成人免费大全| 黄色女人牲交| 国产精品久久久久久人妻精品电影| 色精品久久人妻99蜜桃| 亚洲五月天丁香| 嫩草影院精品99| 欧美成狂野欧美在线观看| 国产精品 欧美亚洲| 免费观看精品视频网站| 精品国产国语对白av| 亚洲午夜精品一区,二区,三区| 精品人妻1区二区| 亚洲 欧美 日韩 在线 免费| 亚洲真实伦在线观看| 久久久久久人人人人人| 国产亚洲精品第一综合不卡| 十八禁网站免费在线| 国产亚洲精品久久久久5区| 99精品久久久久人妻精品| 成熟少妇高潮喷水视频| 中文字幕高清在线视频| 久久久久久人人人人人| 国产一区二区激情短视频| 给我免费播放毛片高清在线观看| 三级毛片av免费| 日本成人三级电影网站| 成年人黄色毛片网站| 变态另类丝袜制服| 俄罗斯特黄特色一大片| 91国产中文字幕| 十分钟在线观看高清视频www| 日本三级黄在线观看| 十八禁人妻一区二区| 人成视频在线观看免费观看| www.999成人在线观看| 午夜免费鲁丝| 欧美最黄视频在线播放免费| 欧美日本视频| 国产精品国产高清国产av| 一区二区三区高清视频在线| 9191精品国产免费久久| 身体一侧抽搐| 在线视频色国产色| 一二三四在线观看免费中文在| 久久久久久人人人人人| 国产精品亚洲一级av第二区| 身体一侧抽搐| 欧美黑人精品巨大| 亚洲激情在线av| 免费高清在线观看日韩| 人人妻,人人澡人人爽秒播| 99久久久亚洲精品蜜臀av| 亚洲国产欧美网| 亚洲午夜精品一区,二区,三区| 欧美乱色亚洲激情| 中文字幕精品免费在线观看视频| 在线观看www视频免费| 嫩草影视91久久| 侵犯人妻中文字幕一二三四区| 搞女人的毛片| 妹子高潮喷水视频| 欧美又色又爽又黄视频| 女性生殖器流出的白浆| 久久这里只有精品19| 9191精品国产免费久久| 韩国精品一区二区三区| 国产午夜精品久久久久久| 别揉我奶头~嗯~啊~动态视频| 中文在线观看免费www的网站 | 日韩欧美三级三区| 69av精品久久久久久| 欧美中文日本在线观看视频| 国产1区2区3区精品| 亚洲真实伦在线观看| 久久中文字幕人妻熟女| 在线观看免费午夜福利视频| 看免费av毛片| 又黄又爽又免费观看的视频| av在线播放免费不卡| 国产精品av久久久久免费| 日本a在线网址| 亚洲第一欧美日韩一区二区三区| 国产免费男女视频| 亚洲真实伦在线观看| 中文字幕另类日韩欧美亚洲嫩草| 熟女少妇亚洲综合色aaa.| 最近最新中文字幕大全电影3 | 久久国产精品男人的天堂亚洲| 亚洲中文字幕一区二区三区有码在线看 | av电影中文网址| 国产男靠女视频免费网站| 国产精品爽爽va在线观看网站 | 99精品在免费线老司机午夜| 欧美日韩黄片免| 男人舔女人的私密视频| 88av欧美| 久久久国产成人免费| 成人精品一区二区免费| 久久欧美精品欧美久久欧美| 啦啦啦 在线观看视频| av福利片在线| 亚洲男人的天堂狠狠| 91老司机精品| 午夜两性在线视频| 国产成年人精品一区二区| 国产亚洲欧美在线一区二区| 麻豆成人午夜福利视频| 热re99久久国产66热| 少妇粗大呻吟视频| 女性生殖器流出的白浆| 日本撒尿小便嘘嘘汇集6| 一区二区三区国产精品乱码| 免费人成视频x8x8入口观看| 国产精品精品国产色婷婷| 成人手机av| 黄色片一级片一级黄色片| 在线十欧美十亚洲十日本专区| 在线观看一区二区三区| 51午夜福利影视在线观看| 91大片在线观看| 婷婷精品国产亚洲av在线| 最近最新免费中文字幕在线| 欧美中文综合在线视频| 18禁观看日本| 99久久无色码亚洲精品果冻| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美精品综合一区二区三区| 男人舔奶头视频| 国产精品日韩av在线免费观看| 成人国语在线视频| 精品国产亚洲在线| 亚洲九九香蕉| 色综合欧美亚洲国产小说| 婷婷亚洲欧美| 亚洲国产看品久久| 嫩草影院精品99| 黑丝袜美女国产一区| 男女午夜视频在线观看| 国产成年人精品一区二区| a在线观看视频网站| 国产99久久九九免费精品| 成人国产一区最新在线观看| 性色av乱码一区二区三区2| 欧美国产日韩亚洲一区| 国内少妇人妻偷人精品xxx网站 | 日韩大尺度精品在线看网址| 成人av一区二区三区在线看| www.自偷自拍.com| 午夜亚洲福利在线播放| 免费在线观看亚洲国产| 久久久国产成人免费| 亚洲精品国产一区二区精华液| 亚洲av成人一区二区三| 中文字幕人成人乱码亚洲影| a级毛片在线看网站| 欧美日韩亚洲国产一区二区在线观看| 黄色a级毛片大全视频| 欧美丝袜亚洲另类 | 精品无人区乱码1区二区| 黄色丝袜av网址大全| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av香蕉五月| 久久久久九九精品影院| 国产区一区二久久| 国产精品亚洲av一区麻豆| 免费在线观看视频国产中文字幕亚洲| 亚洲专区中文字幕在线| 亚洲真实伦在线观看| 在线观看66精品国产| 久久精品夜夜夜夜夜久久蜜豆 | 日本黄色视频三级网站网址| 他把我摸到了高潮在线观看| 18禁美女被吸乳视频| 啦啦啦免费观看视频1| 日韩三级视频一区二区三区| 亚洲成av人片免费观看| 精品一区二区三区av网在线观看| 老司机靠b影院| 午夜福利免费观看在线| 欧美另类亚洲清纯唯美| 一区二区三区高清视频在线| 国产精品久久久久久人妻精品电影| 日日夜夜操网爽| 久久狼人影院| 男女之事视频高清在线观看| 亚洲人成网站高清观看| 午夜两性在线视频| 老汉色av国产亚洲站长工具| 欧美色欧美亚洲另类二区| 亚洲五月色婷婷综合| 国产黄a三级三级三级人| 久久久精品国产亚洲av高清涩受| 亚洲国产精品合色在线| 搡老岳熟女国产| 91国产中文字幕| 国产伦一二天堂av在线观看| 亚洲第一青青草原| 很黄的视频免费| 久热爱精品视频在线9| 亚洲精品粉嫩美女一区| 色综合欧美亚洲国产小说| 国产午夜福利久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 精品久久久久久久久久免费视频| 欧美又色又爽又黄视频| 日韩大码丰满熟妇| 热re99久久国产66热| 久久99热这里只有精品18| 人人妻人人澡人人看| 婷婷精品国产亚洲av| 日本精品一区二区三区蜜桃| 亚洲成av人片免费观看| 亚洲av熟女| xxx96com| 波多野结衣av一区二区av| 男人舔奶头视频| 国产亚洲精品久久久久久毛片| 亚洲国产欧洲综合997久久, | 国产午夜福利久久久久久| 日本免费a在线| xxx96com| 悠悠久久av| 久久精品国产清高在天天线| 午夜精品久久久久久毛片777| 午夜激情av网站| 精品免费久久久久久久清纯| 久久 成人 亚洲| av欧美777| 一进一出抽搐动态| 欧美成人午夜精品| 久久精品国产亚洲av高清一级| 香蕉久久夜色| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 91在线观看av| 香蕉av资源在线| 国产亚洲欧美98| 精品少妇一区二区三区视频日本电影| 亚洲精品美女久久久久99蜜臀| 婷婷丁香在线五月| 在线观看66精品国产| 免费高清在线观看日韩| 亚洲欧美日韩高清在线视频| 午夜福利18| 日韩av在线大香蕉| 50天的宝宝边吃奶边哭怎么回事| 天天躁狠狠躁夜夜躁狠狠躁| 禁无遮挡网站| 天天一区二区日本电影三级| 巨乳人妻的诱惑在线观看| 18美女黄网站色大片免费观看| 久久久久久久久中文| 91成人精品电影| 日韩 欧美 亚洲 中文字幕| 午夜老司机福利片| 少妇熟女aⅴ在线视频| 黄色视频不卡| 亚洲天堂国产精品一区在线| 精品国产一区二区三区四区第35| 黄频高清免费视频| 大型av网站在线播放| 中文在线观看免费www的网站 | 欧美日韩精品网址| 亚洲成av片中文字幕在线观看| 亚洲男人的天堂狠狠| av天堂在线播放| 欧美成人一区二区免费高清观看 | 精品久久久久久久末码| 精品国产乱子伦一区二区三区| 久久人人精品亚洲av| 一级作爱视频免费观看| 美女国产高潮福利片在线看| 天天添夜夜摸| 久久青草综合色| 日韩欧美一区二区三区在线观看| 成人国语在线视频| 亚洲精品在线美女| 色精品久久人妻99蜜桃| 精品高清国产在线一区| 变态另类丝袜制服| 国产在线精品亚洲第一网站| 啦啦啦免费观看视频1| 老司机午夜十八禁免费视频| 丰满人妻熟妇乱又伦精品不卡| svipshipincom国产片| 一级黄色大片毛片| 黑丝袜美女国产一区| 91大片在线观看| 亚洲第一av免费看| 一进一出抽搐动态| 亚洲精品国产区一区二| av超薄肉色丝袜交足视频| 亚洲欧美激情综合另类| 人妻丰满熟妇av一区二区三区| 精华霜和精华液先用哪个| 非洲黑人性xxxx精品又粗又长| 亚洲 欧美 日韩 在线 免费| 人成视频在线观看免费观看| 啦啦啦观看免费观看视频高清| av欧美777| 欧美黄色淫秽网站| 久久精品人妻少妇| 国产熟女午夜一区二区三区| 老熟妇乱子伦视频在线观看| 久久精品国产亚洲av高清一级| 不卡av一区二区三区| 国产精品精品国产色婷婷| 99国产精品99久久久久| 久久久久久国产a免费观看| 精品一区二区三区av网在线观看| 久久久久久久久久黄片| 在线播放国产精品三级| 日韩精品免费视频一区二区三区| 亚洲欧美激情综合另类| 麻豆久久精品国产亚洲av| 久久久久久久精品吃奶| 国产亚洲精品一区二区www| 亚洲国产欧洲综合997久久, | 免费av毛片视频| 热99re8久久精品国产| 大型黄色视频在线免费观看| 免费看十八禁软件| 亚洲国产日韩欧美精品在线观看 | 久久亚洲真实| 国产精品亚洲美女久久久| 亚洲av成人av| 中文字幕最新亚洲高清| av在线播放免费不卡| www国产在线视频色| 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| 国语自产精品视频在线第100页| 两个人视频免费观看高清| 亚洲最大成人中文| 在线免费观看的www视频| 啪啪无遮挡十八禁网站| 在线国产一区二区在线| 国产又色又爽无遮挡免费看| 欧美中文日本在线观看视频| 麻豆久久精品国产亚洲av| 十八禁人妻一区二区| 一边摸一边做爽爽视频免费| a级毛片在线看网站| 国语自产精品视频在线第100页| а√天堂www在线а√下载| 午夜日韩欧美国产| 少妇被粗大的猛进出69影院| 熟女少妇亚洲综合色aaa.| 欧美黑人精品巨大| 国产片内射在线| 丰满的人妻完整版| 最近最新中文字幕大全电影3 | 久久久精品欧美日韩精品| 成人永久免费在线观看视频| 男女之事视频高清在线观看| 国产成+人综合+亚洲专区| 欧美性猛交黑人性爽| 亚洲欧美日韩高清在线视频| 久热这里只有精品99| 成年女人毛片免费观看观看9| 久久久国产欧美日韩av| 操出白浆在线播放| 一进一出抽搐gif免费好疼| 青草久久国产| 夜夜看夜夜爽夜夜摸| 久久天躁狠狠躁夜夜2o2o| 亚洲va日本ⅴa欧美va伊人久久| 露出奶头的视频| 看片在线看免费视频| 久久久精品欧美日韩精品| 日日爽夜夜爽网站| 免费高清视频大片| 午夜免费鲁丝| 国产欧美日韩精品亚洲av| 国内精品久久久久精免费| 免费在线观看完整版高清| 高清在线国产一区| 欧美日本亚洲视频在线播放| 亚洲欧美激情综合另类| 精品少妇一区二区三区视频日本电影| 亚洲午夜精品一区,二区,三区| 90打野战视频偷拍视频| 我的亚洲天堂| 欧美黄色片欧美黄色片| 精品久久久久久久久久免费视频| 男女床上黄色一级片免费看| 久久久久久久久免费视频了| 日韩欧美一区视频在线观看| 午夜福利成人在线免费观看| 午夜福利一区二区在线看| 久久热在线av| 在线天堂中文资源库| 日日摸夜夜添夜夜添小说| 欧美性猛交黑人性爽| 麻豆一二三区av精品| 久久久久久久午夜电影| 久久精品国产综合久久久| 叶爱在线成人免费视频播放| 欧美精品啪啪一区二区三区| 亚洲国产欧美一区二区综合| 黄色成人免费大全| 欧美黑人巨大hd| 色老头精品视频在线观看| 制服丝袜大香蕉在线| 久久久久国产精品人妻aⅴ院| 99久久无色码亚洲精品果冻| 一级黄色大片毛片| 香蕉国产在线看| 国产爱豆传媒在线观看 | 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 99久久国产精品久久久| 午夜福利在线在线| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 俺也久久电影网| 麻豆久久精品国产亚洲av| 日韩有码中文字幕| 人人妻人人澡欧美一区二区| 老司机福利观看| 宅男免费午夜| 日韩大尺度精品在线看网址| 99在线人妻在线中文字幕| 国产99久久九九免费精品| 色综合站精品国产| 91在线观看av| 免费看美女性在线毛片视频| 欧美 亚洲 国产 日韩一| 69av精品久久久久久| 在线观看免费视频日本深夜| 麻豆国产av国片精品| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 巨乳人妻的诱惑在线观看| 老司机福利观看| 亚洲激情在线av| 99热这里只有精品一区 | 国产成人一区二区三区免费视频网站| 成人国产一区最新在线观看| 色尼玛亚洲综合影院| 少妇被粗大的猛进出69影院| 看免费av毛片| 国产一区二区三区在线臀色熟女| 欧美性长视频在线观看| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| 1024手机看黄色片| 神马国产精品三级电影在线观看 | 成人国语在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 久久欧美精品欧美久久欧美| 免费在线观看成人毛片| 人人妻人人澡欧美一区二区| 亚洲精品久久成人aⅴ小说| 欧美成人一区二区免费高清观看 | 18禁裸乳无遮挡免费网站照片 | 日韩三级视频一区二区三区| 欧美最黄视频在线播放免费| 香蕉丝袜av| 夜夜看夜夜爽夜夜摸|