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

    林窗對凋落物分解過程中細菌群落結(jié)構(gòu)和多樣性的影響

    2015-12-07 06:43:19張明錦張健紀托未劉華李勛張艷楊萬勤陳良華
    生態(tài)環(huán)境學報 2015年8期
    關(guān)鍵詞:林窗類群均勻度

    張明錦,張健,紀托未,劉華,李勛,張艷,楊萬勤,陳良華*

    1. 長江上游林業(yè)生態(tài)工程重點實驗室//生態(tài)林業(yè)研究所,四川農(nóng)業(yè)大學,四川 成都 611130;2. 高山森林生態(tài)系統(tǒng)定位研究站//華西雨屏區(qū)生態(tài)環(huán)境監(jiān)測站,四川 成都 611130;3. 長江上游生態(tài)安全協(xié)同創(chuàng)新中心,四川 成都 611130

    林窗對凋落物分解過程中細菌群落結(jié)構(gòu)和多樣性的影響

    張明錦1,2,3,張健1,2,3,紀托未1,2,3,劉華1,2,3,李勛1,2,3,張艷1,2,3,楊萬勤1,2,3,陳良華1,2,3*

    1. 長江上游林業(yè)生態(tài)工程重點實驗室//生態(tài)林業(yè)研究所,四川農(nóng)業(yè)大學,四川 成都 611130;2. 高山森林生態(tài)系統(tǒng)定位研究站//華西雨屏區(qū)生態(tài)環(huán)境監(jiān)測站,四川 成都 611130;3. 長江上游生態(tài)安全協(xié)同創(chuàng)新中心,四川 成都 611130

    基于細菌16S rDNA的PCR-DGGE方法,以林下(US)為對照,研究了長江上游低山丘陵區(qū)人為采伐形成的馬尾松(Pinus massoniana)人工林7種不同大小林窗(G1:100 m2、G2:225 m2、G3:400 m2、G4:625 m2、G5:900 m2、G6:1 225 m2、G7:1 600 m2)在冬夏兩季對凋落物分解過程中細菌群落結(jié)構(gòu)和多樣性的影響。結(jié)果表明,凋落物中的細菌群落結(jié)構(gòu)和多樣性受到林窗大小的顯著影響,這種差異性受到季節(jié)的影響。有些菌綱只出現(xiàn)在特定大小的林窗,例如,綠彎菌綱(Chloroflexi)和放線菌綱(Actinobacteria)都只出現(xiàn)在US和G1林窗,芽孢桿菌綱(Bacilli)只出現(xiàn)在大中型(G4~G7)林窗,黃桿菌綱(Flavobacteria)只出現(xiàn)在中型(G4和G5)林窗,浮霉菌綱(Planctomycetacia)只存在于G5、G7和US。冬夏兩季細菌群落結(jié)構(gòu)差異較大,冬季的細菌類群高于夏季。在冬季,細菌群落的多樣性與林窗大小、凋落物含水率、日平均溫度、凋落物的全碳、可溶性有機碳和全氮含量顯著相關(guān);在夏季,細菌群落的多樣性僅與日平均溫度和全碳含量顯著相關(guān)??傮w來看,變形菌門是凋落物分解過程中的優(yōu)勢類群,有利于凋落物分解過程中的碳氮循環(huán)。林窗的形成改變了凋落物中的細菌群落結(jié)構(gòu)和組成,減小了凋落物分解過程中細菌群落的豐富度指數(shù)、多樣性指數(shù)和均勻度指數(shù),但增大了優(yōu)勢度指數(shù)。冬季的細菌群落的豐富度、多樣性和均勻度高于夏季,優(yōu)勢度則低于夏季。

    林窗;凋落物分解;細菌群落;多樣性;PCR-DGGE

    微生物群落是凋落物分解過程中重要的生物因素,在陸地生態(tài)系統(tǒng)中的生物地球化學循環(huán)中扮演著重要的角色。已有研究表明,不同微生物群落的生長和發(fā)展可以改變土壤有機質(zhì)的周轉(zhuǎn)和命運(Creamer et al.,2015)。微生物的活動推動著土壤中碳和能量的流動,有研究表明土壤微生物具有很高的多樣性,但相關(guān)的研究較少(Prosser,2002)。微生物群落主要包括真菌和細菌,其豐富度(Bekku et al.,2004)和組成(Lipson et al.,2002;Schadt et al.,2003;Lipson,2007)具有較大的時空變化。這些變化與溫度(Wardle,1998)、水分(Michelsen et al.,2004)和養(yǎng)分有效性(Lipson et al.,2000)等非生物因素變化以及植物養(yǎng)分競爭(Bardgett et al.,2002)、植物與微生物互利共生(Wardle et al.,2004)等生物因素密切相關(guān)。同樣地,微生物群落的變化對凋落物分解、土壤養(yǎng)分循環(huán)的速度和養(yǎng)分的有效性也有影響(Corre et al.,2002)??梢姡⑸锶郝涞淖兓瘡娏矣绊懼鷳B(tài)系統(tǒng)的功能和生產(chǎn)力(van der Heijden et al.,2008)。微生物群落的變化通??梢詺w因于水分和溫度的變化,以及碳和氮的可用性等方面(Zhang et al.,1995)。

    林窗是森林生態(tài)系統(tǒng)中普遍存在的干擾方式之一,也是森林更新和演替的重要驅(qū)動力之一(Splechtna et al.,2005)。林窗的存在會直接影響太陽輻射量和到達地面的降水量。這些變量的差異影響著林窗和周圍林下的空氣溫度、相對濕度、土壤溫度和土壤含水量。與林下相比,林窗內(nèi)具有更高的太陽輻射量,能夠直接導致空氣溫度、土壤溫度的增加和相對濕度下降(Clinton,2003;Abd Latif et al.,2010)。林窗內(nèi)的非生物和生物條件的變化主要取決于林窗的大小和林窗內(nèi)的位置(Nagel et al.,2010)。一些研究發(fā)現(xiàn)較大的林窗由于太陽輻射的增加,因而具有較高的土壤溫度(Ritter,2005)。但也有研究發(fā)現(xiàn)較小的林窗中植被的保溫效應好,因而較更大的林窗具有更高的土壤溫度(Sariyildiz,2008)。土壤水分與林窗大小之間的關(guān)系也并不明確,一些研究發(fā)現(xiàn)為負相關(guān)關(guān)系(Prescott et al.,2003),而有的研究表明兩者沒有關(guān)系(Ritter et al.,2005;Sariyildiz,2008),這些研究的不同結(jié)果可能是因為土壤含水率受較多因素的影響。較大的林窗往往減少了植物的蒸騰作用并且有較高的凈降水量,但同時具有較高的徑流和地面蒸發(fā)量(Gray et al.,2002)??梢?,林窗的形成以及不同大小的林窗改變了森林生態(tài)系統(tǒng)中的微氣候環(huán)境,而微生物對周圍環(huán)境條件的變化非常敏感。雖然目前關(guān)于林窗對氣候因子的研究已取得了大量的成果,但關(guān)于林窗對森林凋落物分解過程中細菌群落結(jié)構(gòu)的影響的研究還很少。凋落葉分解過程中細菌群落結(jié)構(gòu)的變化不僅受到溫度、濕度等環(huán)境因子的影響,而且還可能會受到凋落物基質(zhì)養(yǎng)分有效性的影響。因此,有必要針對林窗對凋落物分解過程中細菌群落結(jié)構(gòu)的影響進行深入研究,以進一步認識林窗中微生物對有機質(zhì)分解的作用以及相關(guān)養(yǎng)分和碳循環(huán)的影響。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究區(qū)域位于長江上游四川省宜賓市高縣境內(nèi) ( 28°35′45″N~28°36′19″N , 104°33′37″E~104°34′22″E),地處岷江、金沙江交匯段,海拔高度410~560 m,是典型的低山丘陵地區(qū);屬于亞熱帶季風性濕潤氣候,雨量充沛,雨熱同季;最高氣溫40.3 ℃,最低氣溫-3.7 ℃,年均氣溫18.1 ℃;年均降水量1021.8 mm,相對濕度81%。試驗地位于來復鎮(zhèn)毛巔坳工區(qū)經(jīng)營小班,地處川云公路中段,總面積約為33.4 hm2。土壤為山地黃壤,成土母巖為黃色砂巖,土壤發(fā)育完善,層次過渡明顯,土壤養(yǎng)分狀況差,缺乏氮和磷。林地為42年生成熟馬尾松人工純林,坡度平緩,坡向S、SE,郁閉度0.7左右。林下灌木主要有枹櫟(Quercus serrata)、黃荊(Vitex negundo)、黃牛奶(Symplocos laurina)、梨葉懸鉤子(Rubus pirifolius)、野桐(Mallotus japonicus)、展毛野牡丹(Melastoma normale)等;草本主要有淡竹葉(Lophatherum gracile)、金星蕨(Parathelypteris glanduligera)、馬唐(Digitaria sanguinalis)、芒(Miscanthus sinensis)、芒萁(Dicranopteris dichotoma)、烏蕨(Stenoloma chusanum)、皺葉狗尾草(Setaria plicata)等。

    1.2 樣地設置

    2011年8月在試驗區(qū)域?qū)拥氐闹参锒鄻有?、土壤理化性質(zhì)、林分結(jié)構(gòu)等作初步調(diào)查。2011年10月在踏查基礎上在實驗區(qū)選擇立地條件相近的馬尾松(Pinus massoniana)人工純林進行采伐,形成7種面積不等的近正方形林窗各3個,面積大小分別為G1:100 m2、G2:225 m2、G3:400 m2、G4:625 m2、G5:900 m2、G6:1225 m2、G7:1600 m2,未砍伐的人工林作為對照-US(林下)。林窗斑塊主要位于山體中部和中上坡位,坡向基本為東,坡度為20°~30°,海拔400 m左右,林窗間距大于10 m。將林窗內(nèi)砍伐的樹枝、樹干等殘體移出,但未進行掘根處理。樣地基本情況見表1。

    1.3 動態(tài)采樣

    在2013年10月下旬凋落物高峰期,在研究區(qū)域收集人工林下新鮮的馬尾松凋落葉,帶回實驗室風干。實驗處理采用凋落物袋分解法,凋落物袋大小約為20 cm×20 cm,孔徑0.5 mm。將風干的馬尾松凋落葉按每袋10 g裝入凋落物分解袋并封口,在7種不同大小林窗的林窗中心和林窗邊緣2個位置和林下對照放置凋落物分解袋,每種林窗大小和林下對照都設置了3個重復樣地,計劃進行7次采樣,每次采樣在每個放置點采 3袋,所以共計(7×2×3+3)×3×7=945袋。于2013年11月中旬,在 7個大小梯度的林窗中心和林窗邊緣及林下對照,去除樣地表面凋落物層,將凋落物樣袋放置于土壤表面,讓凋落葉自然分解。在7個林窗的中心、邊緣及林下的凋落物袋中放入紐扣式溫度記錄器(Button DS1923-F5,Maxim/Dallas Semiconductor,Sunnyvale,CA,USA),每2 h記錄一次溫度數(shù)據(jù)。分別在冬季(2014年2月中旬)和夏季(2014年8月中旬)進行采樣,在每個設置的樣點各采集3袋凋落物分解袋。將樣品裝入已滅菌的封口聚乙烯袋,并標記,迅速放入冰盒中帶回實驗室。除去樣品中的土壤、根系等雜質(zhì)后,剪碎過2 mm篩,一部分裝入已滅菌的封口聚乙烯袋于 4 ℃保存用于測定細菌群落結(jié)構(gòu),另一部分烘干用于凋落物含水率、全碳、全氮、全磷和可溶性有機碳的測定。

    1.4 測定方法

    凋落物 C、N、P含量的測定參照《中華人民共和國林業(yè)行業(yè)標準 LY/T─1999》。全碳用重鉻酸鉀氧化外加熱法測定。樣品用硫酸-高氯酸消煮,用凱式定氮儀(浙江托普儀器有限公司KDN系列)測定全N,用鉬銻抗比色法測定全P。DOC含量的測定步驟如下(王春陽等,2011):將烘干的植物樣品粉碎過2 mm篩,精確稱取植物樣品0.05 g于20 mL離心管中,加入10 mL蒸餾水,振蕩0.5 h,4000 r·min-1、4 ℃離心 15 min,取上清液在-0.08 MPa下經(jīng)孔徑為0.45 μm濾膜(上海市新亞凈化器件廠生產(chǎn),使用前首先要用蒸餾水淋洗濾膜,以淋洗濾膜可能殘留的有機組分)過濾,即獲得植物凋落物DOC提取液,4 ℃低溫保存待測液,用TOC分析儀(multi N/C 2100,AnalytikJena)測定待測液中DOC濃度。

    凋落物中細菌基因組DNA的提取采用北京艾德萊生物科技有限公司(Aidlab Biotechnologies Co.,Ltd.,China)的SoilPure試劑盒,按操作說明提取并純化,純化后的DNA在1%瓊脂糖凝膠上進行電泳檢測。采用細菌通用引物對總基因組 16S rDNA的V3區(qū)進行PCR擴增,其引物序列為:F338:5′GC-ACTCCTACGGGAGGCAGCAG3′;R518:5′ATTACCGCGGCTGCTGG3′;F338的5′GC夾子結(jié)構(gòu)為:5′CGCCCGCCGCGCCCCGCGCCCGG CCCGCCGCCCCCGCCCC3′。PCR反應體系為:25 ng模板、2.5 pmol正反引物,10 μg Bovine Serum Albumin(TaKaRa)、12.5 μL Premix Taq(TaKaRa)、無菌純水補足到25 μL。反應程序為:95 ℃預變性3 min,94 ℃變性30 s,55 ℃退火45 s,72 ℃延伸1 min,30個循環(huán),循環(huán)結(jié)束后72 ℃延伸5 min。用1.5%瓊脂糖凝膠電泳檢測PCR擴增產(chǎn)物。

    PCR 產(chǎn)物檢測后,通過 CBS Denaturing Gradient Gel Electrophoresis System(CBS Scientific Company,Inc)電泳儀進行 DGGE電泳分離。電泳緩沖液為 1×TAE,膠濃度為 8%,變性梯度為40%~70%。PCR產(chǎn)物和6×LoadingBuffer按照5∶1混合均勻后點樣。電泳條件設置:電壓為100 V,溫度為60 ℃,時間為930 min,電泳結(jié)束后進行銀染,并用 GS-800TM Calibrated Densitometer(Bio-Rad Laboratories,Hercules,USA)進行掃描成像。使用Quantity One(Bio-Rad,Hercules,USA)軟件對 DGGE圖譜進行細菌群落組成結(jié)構(gòu)和多樣性的分析。選擇電泳膠上的優(yōu)勢條帶進行切膠回收,回收DNA用不含GC夾子的引物進行PCR擴增,擴增后進行純化。純化產(chǎn)物連接到pMD19-T載體(TaKaRa),再轉(zhuǎn)入大腸桿菌感受態(tài)細胞DH5α在培養(yǎng)皿中進行培養(yǎng),挑陽性單克隆送華大基因研究院進行測序,以確定優(yōu)勢條帶所對應的細菌種屬。

    1.5 數(shù)據(jù)處理與統(tǒng)計分析

    DGGE主要條帶的測序結(jié)果去除載體序列后,在GenBank中進行BLAST分析,找出相似性大于95%的相似序列,以確定條帶所對應的細菌類群。用Quantity One軟件分析出各泳道條帶的相對強度后,結(jié)合Blast比對的測序結(jié)果,計算細菌類群在各樣品中的相對比例。用Quantity One軟件對DGGE圖譜進行分析,通過計算各個樣品的Patrick豐富度指數(shù)即條帶數(shù)(S)、Shannon多樣性指數(shù)(H)、Pielou均勻度指數(shù)(EH)、Simpson優(yōu)勢度指數(shù)(C)來評價樣品中細菌群落的多樣性。其中,計算公式如下:

    式中,S表示每個泳道檢測到的不同條帶的總和,Pi表示每個泳道中各個條帶的強度占總強度的百分含量。

    采用SPSS 17.0作統(tǒng)計分析,用單因素方差分析方法(one-way ANOVA)和Duncan(D)多重比較法檢驗不同季節(jié)凋落物含水率在不同大小林窗間的差異。用配對樣本T檢驗(Paired Sample T test)檢驗不同大小林窗下不同林窗位置處凋落物含水率在冬夏兩季之間的差異。細菌多樣性指數(shù)與環(huán)境因子的相關(guān)關(guān)系采用Pearson法進行分析。數(shù)據(jù)整理和計算采用Microsoft Excel 2007軟件,作圖采用Origin 9.0。

    表1 馬尾松人工林不同面積林窗樣地基本概況Table 1 General characteristics in the sampling plots in different size forest gaps of Pinus massoniana plantations

    2 結(jié)果與分析

    2.1 日平均溫度和凋落物含水率

    從采樣當月的日平均溫度來看(表2),在冬季,林窗中心的最低和最高日平均溫度分別是 9.27 ℃(G2)和10.15 ℃(G7),林窗邊緣是8.78 ℃(G5)和10.67 ℃(G6)。在G1~G5林窗,林窗中心的日平均溫度高于林窗邊緣,但在G6和G7,林窗中心的日平均溫度低于林窗邊緣。在夏季,林窗中心的最低和最高日平均溫度分別是 23.81 ℃(G2)和28.10 ℃(G7),林窗邊緣是 22.97 ℃(G2)和26.13 ℃(G3)。林窗中心的日平均溫度都高于林窗邊緣。

    表2 采樣當月的日平均溫度Table 2 Average daily temperature in sampling month ℃

    從凋落物含水率來看(圖 1所示),在冬季,無論在林窗中心還是林窗邊緣,大型林窗(G6和G7)中的凋落物含水率都只有15.43%~19.78%,顯著低于林下和 G1~G4林窗。在夏季,凋落物的含水率隨林窗大小無顯著變化。在G5~G7和G3林窗中心和邊緣,以及G1林窗中心,凋落物在冬季的含水率顯著低于夏季,其余無顯著差異。

    圖1 凋落物含水率Fig. 1 Water content of litter

    2.2 細菌群落組成

    如圖2所示,從已知的細菌類群來看,在冬季,凋落物中的細菌主要歸屬于變形菌門( Proteobacteria) 中 的 α-變 形 菌 綱( Alphaproteobacteria )、 β-變 形 菌 綱( Betaproteobacteria ) 、 γ-變 形 菌 綱(Gammaproteobacteria),酸桿菌門(Acidobacteria)的酸桿菌綱(Acidobacteria),異常球菌-棲熱菌門( Deinococcus-Thermus) 的 異 常 球 菌 綱(Deinococci),擬桿菌門(Bacteroidetes)的鞘脂桿 菌 綱 ( Sphingobacteria) 、 黃 桿 菌 綱(Flavobacteria),厚壁菌門(Firmicutes)的芽孢桿菌綱(Bacilli),綠彎菌門(Chloroflexi)的綠彎菌綱(Chloroflexi),放線菌門(Actinobacteria)的放線菌綱(Actinobacteria),共7大門10大綱。在林窗中心(圖 2a),已知細菌類群只占了16.99%~40.95%。在林窗邊緣(圖2b),已知細菌類群也只有29.12%~65.56%,多數(shù)類群為未知類群。從已知的細菌類群來看,無論是在林窗中心還是林窗邊緣,除了G3林窗是以γ-變形菌綱為優(yōu)勢類群外,其余林窗中心和邊緣以及林下都是以α-變形菌綱為優(yōu)勢類群。在林窗中心,β-變形菌綱只出現(xiàn)在G2~G5林窗;在林窗邊緣,β-變形菌綱也只出現(xiàn)在G2、G3和G5林窗。黃桿菌綱只出現(xiàn)在G4林窗中心和邊緣,以及G5林窗中心。芽孢桿菌綱只出現(xiàn)在G4、G6、G7林窗中心和邊緣,以及G5林窗中心。只有G1林窗中心和邊緣以及US出現(xiàn)了綠彎菌綱和放線菌綱。

    圖2 冬季不同林窗內(nèi)凋落物中不同細菌類群的相對含量Fig. 2 Percentage of relative content of different bacterial groups of litter in different forest gap in winter

    如圖3所未,從已知的細菌類群來看,夏季凋落物中的細菌主要歸屬于變形菌門(Proteobacteria)的 α-變形菌綱(Alphaproteobacteria)、β-變形菌綱( Betaproteobacteria ) 、 γ- 變 形 菌 綱(Gammaproteobacteria),擬桿菌門(Bacteroidetes)的鞘脂桿菌綱(Sphingobacteria),放線菌門(Actinobacteria)的放線菌綱(Actinobacteria),浮霉 菌 門 ( Planctomycetes) 的 浮 霉 菌 綱(Planctomycetacia),酸桿菌門(Acidobacteria)的酸桿菌綱(Acidobacteria),共5大門7大綱。林窗中心的已知類群占細菌群落的11.14%~57.82%(圖3a),林窗邊緣的已知類群為28.37%~56.59%(圖3b),多數(shù)類群仍為未知類群。從已知的細菌類群來看,在US和G5林窗中心,放線菌綱為優(yōu)勢類群;在G1、G3和G6林窗中心,α-變形菌綱為優(yōu)勢類群;在G2林窗中心,β-變形菌綱占優(yōu)勢類群;在G4和G7林窗中心,γ-變形菌綱所占比例最高。在林窗邊緣,除了G5和G6林窗的優(yōu)勢類群為鞘脂桿菌綱外,其余林窗邊緣的優(yōu)勢類群都是α-變形菌綱。酸桿菌綱只在G1和G3林窗邊緣以及G4林窗中心和邊緣出現(xiàn),浮霉菌綱只出現(xiàn)在 G5和 G7林窗中心以及US。

    圖3 季不同林窗內(nèi)凋落物中不同細菌類群的相對含量Fig. 3 Percentage of relative content of different bacterial groups of litter in different forest gap in summer

    表3 細菌群落多樣性指數(shù)Table 3 bacterial community diversity index

    2.3 細菌群落多樣性

    如表3所示,在冬季,無論在林窗中心還是林窗邊緣,細菌群落的Patrick豐富度指數(shù)和Shannon多樣性指數(shù)除了G1高于林下外,其余林窗都低于林下,且隨著林窗面積的增大,兩者均先減小后增大,其中G4林窗的最小。Simpson優(yōu)勢度指數(shù)與前兩者相反,除了G1林窗略低于林下外,其余林窗都高于林下,且隨著林窗面積的增大先增大后減小,在 G4林窗達到最大。在林窗中心,除了 G1和G6以外,其余林窗的Pielou均勻度指數(shù)都低于林下;在林窗邊緣,除了G1、G6和G7以外,其余林窗也低于林下。在夏季,在林窗中心,除了G4和G6的Patrick豐富度指數(shù)高于林下外,其余林窗都低于林下,Shannon多樣性指數(shù)除了G4高于林下外,其余林窗都低于林下。林窗中心的Pielou均勻度指數(shù)都低于林下,但Simpson優(yōu)勢度指數(shù)都高于林窗。在林窗邊緣,只有G4和G5的Patrick豐富度指數(shù)低于林下。除了 G2和 G3林窗的Shannon多樣性指數(shù)高于林下外,其余林窗都低于林下。Simpson優(yōu)勢度指數(shù)與此相反,除了G2和G3低于林下外,其余林窗都高于林下。在林窗邊緣,除了G3林窗外,其余林窗的Pielou均勻度指數(shù)都低于林下??傮w來看,林窗的形成減小了凋落物分解過程中細菌群落的 Patrick豐富度指數(shù)、Shannon多樣性指數(shù)和Pielou均勻度指數(shù),但增大了Simpson優(yōu)勢度指數(shù)。從季節(jié)變化來看,夏季的patrick豐富度指數(shù)、Shannon多樣性指數(shù)和Pielou均勻度指數(shù)低于冬季,但Simpson優(yōu)勢度指數(shù)高于冬季。

    2.4 細菌群落多樣性與環(huán)境因子和凋落物化學質(zhì)量的相關(guān)分析

    從相關(guān)性分析結(jié)果來看(表4),在冬、夏兩季,環(huán)境因子和凋落物化學質(zhì)量對細菌群落結(jié)構(gòu)多樣性的影響差異較大。在冬季,豐富度和多樣性指數(shù)與林窗大小顯著相關(guān),豐富度指數(shù)還與WC、DOC、N顯著相關(guān),多樣性指數(shù)與DOC、N呈顯著正相關(guān),均勻度指數(shù)與WC呈顯著負相關(guān),與AT、TC呈顯著正相關(guān),優(yōu)勢度指數(shù)與N呈顯著負相關(guān)。在夏季,細菌群落結(jié)構(gòu)的豐富度、多樣性、均勻度和優(yōu)勢度指數(shù)與林窗大小的相關(guān)性都不顯著,但都與AT顯著相關(guān),均勻度指數(shù)還與TC顯著相關(guān)。

    表4 細菌群落多樣性指數(shù)的相關(guān)性分析Table 4 Correlation analysis of bacterial community diversity index

    3 討論與結(jié)論

    3.1 討論

    林窗大小是林窗的重要特征之一,不同大小的林窗干擾可能產(chǎn)生不同的森林微環(huán)境(光照、溫度及水分等)(Galhidy et al.,2006),從而影響到土壤養(yǎng)分有效性以及土壤微生物群落等(Muscolo et al.,2007)。在冬季,細菌群落的豐富度指數(shù)和多樣性指數(shù)與林窗大小的相關(guān)性都極顯著(表4)。一些菌綱的細菌類群只出現(xiàn)在特定大小的林窗內(nèi)。這可能與冬季凋落物的含水率隨林窗大小有顯著變化有關(guān)。細菌群落結(jié)構(gòu)不僅在不同大小林窗間有差異,在冬、夏兩季間的差異也很明顯(圖2和圖3)。從已知類群來看,在冬季凋落物分解中,出現(xiàn)了 7大菌門中的 10個菌綱的細菌類群,但在夏季,只出現(xiàn)了5大門中7個菌綱的細菌類群。冬季的細菌類群高于夏季。一方面,可能是因為不同細菌類群的生態(tài)適應性有所差異(蒲玲玲,2006;鄭賀云等,

    2012)。研究區(qū)域冬夏兩季溫差很大,冬冷夏熱(表2),同時夏季的雨水較多,凋落物的含水率明顯高于冬季(圖1)。另一方面,可能與凋落物養(yǎng)分有效性的改變有關(guān)(Marschner et al.,2011),在本研究中,馬尾松凋落葉中DOC、N、P的平均含量分別從冬季的9.65、10.60、0.45 g·kg-1減少到只有4.40、2.00、0.42 g·kg-1,DOC/TC從0.03減小到0.01,TC/TN從34.75增大到313.99。有研究發(fā)現(xiàn)微生物底物有效性較差,所能支持的微生物類群及數(shù)量相對較少(Taylor et al.,1989)。林窗大小的改變和凋落物分解的不同階段導致的溫度、水分、養(yǎng)分有效性的差異使細菌群落發(fā)生著變化,細菌群落多樣性和結(jié)構(gòu)組成的變化直接關(guān)系到微生物類群的功能和作用(Torsvik et al.,2002;Wellington et al.,2003)。

    從已知類群來看,在冬夏兩季,變形菌門、放線菌門、擬桿菌門和酸桿菌門均存在,這說明在調(diào)落物分解過程中,這4類細菌對環(huán)境的適應能力和凋落物的分解能力較強。在冬季,變形菌門是凋落物分解過程中的優(yōu)勢菌類群。在夏季,變形菌門、擬桿菌門和放線菌門分別在不同大小林窗的中心和邊緣具有優(yōu)勢地位。細菌種類不同其功能也有所差異,細菌群落組成的變化將直接影響到各種細菌類群不同功能的發(fā)揮(Taylor et al.,1989),以及相互間的互作關(guān)系,從而對各種營養(yǎng)元素的循環(huán)、有機物的分解利用、植物對養(yǎng)分的吸收和生長發(fā)育產(chǎn)生影響。關(guān)于變形菌門的研究發(fā)現(xiàn)變形菌門是土壤中較為豐富的一個細菌門(Janssen,2006;楊鑫等,2008),對全球的碳、氮循環(huán)具有極為重要的作用(Kersters et al.,2006)。變形菌門在凋落物分解過程中的優(yōu)勢地位有利于凋落物分解過程中的碳氮循環(huán)。馬尾松分解過程中的擬桿菌主要以黃桿菌綱和鞘脂桿菌綱為主,該類菌群也是調(diào)落物分解過程中的重要細菌類群,這類細菌對凋落物中纖維素的降解起著重要作用(趙野逸,2013)。放線菌能促使凋落物迅速腐爛,在氮素循環(huán)中起著重要作用。在冬季,不同大小林窗內(nèi)凋落物中的優(yōu)勢細菌類群都屬于變形菌門。但到了夏季,不同大小林窗內(nèi)的優(yōu)勢細菌類群出現(xiàn)了差異,且不同細菌類群的相對含量也有很大變化。這說明隨著馬尾松凋落葉分解,不同大小的林窗對凋落物中細菌群落結(jié)構(gòu)有著顯著影響,進而可能會影響凋落物的分解過程。

    從細菌群落多樣性的相關(guān)性分析可以看出,環(huán)境因子和凋落物質(zhì)量與細菌群落結(jié)構(gòu)多樣性的相關(guān)性在冬夏兩季之間差異很大(表4)。在冬季,細菌多樣性不僅與WC和AT顯著相關(guān),同時與TC、DOC和N也有著很強的相關(guān)性。前人的研究也表明土壤溫度(Rinnan et al.,2007;Frey et al.,2008;Feng et al.,2009)、碳氮的可利用性(Hogberg et al.,2007;Ushio et al.,2008)影響著土壤微生物群落。而在夏季,除了均勻度指數(shù)與TC有顯著的相關(guān)性,其余指數(shù)只與AT有顯著的相關(guān)性。WC和AT的變化受林窗大小的影響,而TC、DOC和N含量作為凋落物基質(zhì)的基本特征,隨凋落物分解過程不斷發(fā)生改變。細菌主要通過對凋落物的分解而獲取碳源和氮源,隨著凋落物的分解,凋落物的可溶性有機碳(DOC)、全氮(N)等分解逐漸釋放,含量不斷減少,養(yǎng)分有效性降低,進而導致可供細菌利用的能源物質(zhì)減少,對細菌群落結(jié)構(gòu)和多樣性影響也逐漸減?。ㄖ軙詰c等,2011)

    3.2 結(jié)論

    綜上所述,人工林窗的形成對凋落物中細菌群落的多樣性有著顯著影響。林窗的形成減小了凋落物分解過程中細菌群落的豐富度指數(shù)、多樣性指數(shù)和均勻度指數(shù),但增大了優(yōu)勢度指數(shù)。從季節(jié)變化來看,夏季的豐富度指數(shù)、多樣性指數(shù)和均勻度指數(shù)都低于冬季,但優(yōu)勢度指數(shù)高于冬季。由于冬、夏兩季氣候特征的差異以及凋落物分解引起的自身有機質(zhì)和養(yǎng)分含量的變化導致了冬、夏兩季凋落物分解過程中細菌群落結(jié)構(gòu)和多樣性的差異。相同氣候下不同大小林窗的形成使得在林窗內(nèi)形成了不同的微氣候環(huán)境,微氣候環(huán)境的差異能在一定程度上改變凋落物分解過程中的細菌群落結(jié)構(gòu)和多樣性,進而影響到森林生態(tài)系統(tǒng)中的物質(zhì)和能量循環(huán)。雖然林窗的形成導致了細菌菌落結(jié)構(gòu)的改變,但由于對很多細菌類群的生物學和生態(tài)學特性的研究還相對缺乏,因此林窗與細菌群落之間的關(guān)系有待進一步研究。

    ABD LATIF Z, BLACKBUMl G A. 2010. The effects of gap size on some microclimate variables during late summer and autumn in a temperate broadleaved deciduous forest [J]. International Journal of Biometeorology, 54(2): 119-129.

    BARDGETT R D, STREETER T C, COLE L, et al. 2002. Linkages between soil biota, nitrogen availability, and plant nitrogen uptake in a mountain ecosystem in the Scottish Highlands [J]. Applied Soil Ecology, 19(2): 121-134.

    BEKKU Y S, NAKATSUBO T, KUME A, et al. 2004. Soil microbial biomass, respiration rate, and temperature dependence on a successional glacier foreland in Ny-?lesund, Svalbard [J]. Arctic Antarctic and Alpine Research, 36: 395-399.

    CLINTON B D. 2003. Light, temperature and soil moisture response to elevation, evergreen understory, and small canopy gaps in the southern Appalachians [J]. Forest Ecology and Management, 186(1-3): 243-255.

    CORRE M, SCHNABEL R, STOUT W. 2002. Spatial and seasonal variation of gross nitrogen transformations and microbial biomass in a Northeastern US grassland [J]. Soil Biology and Biochemistry, 34(4): 445-457.

    CREAMER C A, DE MENEZES A B, KRULL E S, et al. 2015. Microbial community structure mediates response of soil C decomposition to litter addition and warming [J]. Soil Biology and Biochemistry, 80: 175-188.

    FENG X, SIMPSON M J. 2009. Temperature and substrate controls on microbial phospholipid fatty acid composition during incubation of grassland soils contrasting in organic matter quality [J]. Soil Biology and Biochemistry, 41(4): 804-812.

    FREY S D, DRIJBER R, SMITH H, et al. 2008. Microbial biomass, functional capacity, and community structure after 12 years of soil warming [J]. Soil Biology and Biochemistry, 40(11): 2904-2907.

    GALHIDY L, MIHOK B, HAGYO A, et al. 2006. Effects of gap size and associated changes in light and soil moisture on the understorey vegetation of Hungarian beech forest [J]. Plant Ecology, 183(1): 133-145.

    GRAY A, SPIES T, EASTER M. 2002. Microclimatic and soil moisture responses to gap formation in coastal Douglas-fir forests [J]. Canadian Journal of Forest Research, 32(2): 332-343.

    HOGBERG M N, HOGBERG P, MYROLD D D. 2007. Is microbial community composition in boreal forest soils determined by pH, C-to-N ratio, the trees, or all three? [J]. Oecologia, 150(4): 590-601.

    JANSSEN P H. 2006. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes [J]. Applied and Environmental Microbiology, 72(3): 1719-1728.

    KERSTERS K, DE V P, GILLIS M, et al. 2006. Introduction to the Proteobacteria [M]. New York: Springer.

    LIPSON D A. 2007. Relationships between temperature responses and bacterial community structure along seasonal and altitudinal gradients [J]. Fems Microbiology Ecology, 59(2): 418-427.

    LIPSON D A, SCHADT C W, SCHMIDT S K. 2002. Changes in soil microbial community structure and function in an alpine dry meadow following spring snow melt [J]. Microbial Ecology, 43(3): 307-314.

    LIPSON D A, SCHMIDT S K, MONSON R K. 2000. Carbon availability and temperature control the post-snowmelt decline in alpine soil microbial biomass [J]. Soil Biology and Biochemistry, 32(4): 441-448.

    MARSCHNER P, UMAR S, BAUMANN K. 2011. The microbial community composition changes rapidly in the early stages of decomposition of wheat residue [J]. Soil Biology and Biochemistry, 43(2): 445-451.

    MICHELSEN A, ANDERSSON M, JENSEN M, et al. 2004. Carbon stocks, soil respiration and microbial biomass in fire-prone tropical grassland, woodland and forest ecosystems [J]. Soil Biology and Biochemistry, 36(11): 1707-1717.

    MUSCOLO A, SIDARI M, MERCURIO R. 2007. Influence of gap size on organic matter decomposition, microbial biomass and nutrient cycle in Calabrian pine (Pinus laricio, Poiret) stands [J]. Forest Ecology and Management, 242(2): 412-418.

    NAGEL T A, SVOBODA M, RUGANI T, et al. 2010. Gap regeneration and replacement patterns in an old-growth Fagus-Abies forest of Bosnia-Herzegovina [J]. Plant Ecology, 208(2): 307-318.

    PRESCOTT C, HOPE G, BLEVINS L. 2003. Effect of gap size on litter decomposition and soil nitrate concentrations in a high-elevation spruce-fir forest [J]. Canadian Journal of Forest Research, 33(11): 2210-2220.

    PROSSER J I. 2002. Molecular and functional diversity in soil microorganisms [J]. Plant and Soil, 244(1/2): 9-17.

    RINNAN R, MICHELSEN A, B??TH E, et al. 2007. Fifteen years of climate change manipulations alter soil microbial communities in a subarctic heath ecosystem [J]. Global ChangeBiology, 13(1): 28-39.

    RITTER E. 2005. Litter decomposition and nitrogen mineralization in newly formed gaps in a Danish beech (Fagus sylvatica) forest [J]. Soil Biology and Biochemistry, 37(7): 1237-1247.

    RITTER E, STARR M, VESTERDAL L. 2005. Losses of nitrate from gaps of different sizes in a managed beech (Fagus sylvatica) forest [J]. Canadian Journal of Forest Research, 35: 308-319.

    SARIYILDIZ T. 2008. Effects of Gap-Size Classes on Long-Term Litter Decomposition Rates of Beech, Oak and Chestnut Species at High Elevations in Northeast Turkey [J]. Ecosystems, 11(6): 841-853.

    SCHADT C W, MARTIN A P, LIPSON D A, et al. 2003. Seasonal dynamics of previously unknown fungal lineages in tundra soils [J]. Science, 301(5638): 1359-1361.

    SPLECHTNA B E, GRATZER G, BLACK B A. 2005. Disturbance history of a European old-growth mixed-species forest——a spatial dendroecological analysis [J]. Journal of vegetation science, 16(5): 511-522.

    TAYLOR B R, PARKINSON D, PARSONS W F J. 1989. Nitrogen and lignin content as predictors of litter decay rates: a microcosm test [J]. Ecology, 70(1): 97-104.

    TORSVIK V, ?VRE?S L. 2002. Microbial diversity and function in soil: from genes to ecosystems [J]. Current Opinion in Microbiology, 3(5): 240-245.

    USHIO M, WAGAI R, BALSER T C, et al. 2008. Variations in the soil microbial community composition of a tropical montane forest ecosystem: Does tree species matter? [J]. Soil Biology and Biochemistry, 40(10): 2699-2702.

    VAN DER HEIJDEN M G A, BARDGETT R D, VAN STRAALEN N M. 2008. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems [J]. Ecology Letters, 11(3): 296-310.

    WARDLE D A. 1998. Controls of temporal variability of the soil microbial biomass: a global-scale synthesis [J]. Soil Biology and Biochemistry, 30(13): 1627-1637.

    WARDLE D A, BARDGETT R D, KLIRONOMOS J N, et al. 2004. Ecological linkages between aboveground and belowground biota [J]. Science, 304(5677): 1629-1633.

    WELLINGTON E M H, BERRY A, KRSEK M. 2003. Resolving functional diversity in relation to microbial community structure in soil: exploiting genomics and stable isotope probing [J]. Current Opinion in Microbiology, 3(6): 295-301.

    ZHANG Q, ZAK J. 1995. Effects of gap size on litter decomposition and microbial activity in a subtropical forest [J]. Ecology, 76(7): 2196-2204.

    蒲玲玲. 2006. 青藏高原冰川與凍土微生物多樣性的研究[D]. 蘭州: 蘭州大學碩士學位論文.

    王春陽, 周建斌, 王祥, 等. 2011. 黃土高原區(qū)不同植物凋落物可溶性有機碳的含量及生物降解特性[J]. 環(huán)境科學, 32(4): 1139-1145.

    楊鑫, 曹靖, 董茂星, 等. 2008. 外來樹種日本落葉松對森林土壤質(zhì)量及細菌多樣性的影響[J]. 應用生態(tài)學報, 19(10): 2109-2116.

    趙野逸. 2013. 雪被對高山森林凋落物分解過程中微生物群落的影響[D].四川: 四川農(nóng)業(yè)大學碩士學位論文.

    鄭賀云, 黎志坤, 李超, 等. 2012. 新疆阿克蘇地區(qū)鹽堿地細菌類群多樣性及優(yōu)勢菌群分析[J]. 微生物學通報, 39(7): 1031-1043.

    周曉慶, 吳福忠, 楊萬勤, 等. 2011. 高山森林凋落物分解過程中的微生物生物量動態(tài)[J]. 生態(tài)學報, 31(14): 4144-4152.

    Influence of Forest Gap on Bacterial Community Structure and Diversity during Litter Decomposition

    ZHANG Mingjin1,2,3, ZHANG Jian1,2,3, JI Tuowei1,2,3, LIU Hua1,2,3, LI Xun1,2,3, ZHANG Yan1,2,3, YANG Wanqin1,2,3, CHEN Lianghua1,2,3*
    1. Key laboratory of Ecological Forestry Engineering of Sichuan Province//Institute of Ecology & Forest, Sichuan Agricultural University, Chengdu 611130, China; 2. Long-term Research Station of Alpine Forest Ecosystems//Monitoring Station for Eco-environments in the Rainy Zone of SW China, Chengdu 611130, China; 3. Collaborative Innovation Center of Ecological Security in the Upper Reaches of Yangtze River, Chengdu 611130, China

    The research on the bacterial community structure and diversity in litter in winter and summer was based on polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) for 16S rDNA in seven gap sizes of a Pinus massoniana plantation (G1: 100 m2, G2: 225 m2, G3: 400 m2, G4: 625 m2, G5: 900 m2, G6: 1 225 m2, G7: 1 600 m2) in a hilly area of the upper Yangtze River. Results indicated that the bacterial community structure and diversity in litter were influenced by gap size, and this difference was also affected by the seasons. Some bacteria classes only appeared in forest gap of certain size, for example, Chloroflexi and Actinobacteria only appeared in US and G1, Bacilli only appeared in large and medium-sized (G4 to G7) forest gap, Flavobacteria only appeared in medium-sized (G4 and G5) forest gap, Planctomycetacia only existed in G5, G7 and US. The difference of bacteria community structure between winter and summer was very big. The bacteria groups of winter were more than summer. In winter, the diversity of bacterial community showed very close correlation with gap size, litter water content, daily average temperature, total carbon, dissolved organic carbon and total nitrogen content of litter; in summer, it was only related to daily average temperature and total carbon content of litter. Overall, Proteobacteria was the dominant in the process of litter decomposition, was conducive to the cycle of carbon and nitrogen in the process of litter decomposition. The formation of forest gap changed bacterial community structure and composition in litter, reduced the richness index, diversity index and evenness index of bacterial community in the process of litter decomposition, but increased the dominance index. Richness, diversity and evenness of bacterial community in winter were higher than those in summer, but dominance was lower than that in summer.

    forest gap; litter decomposition; bacterial community; diversity; PCR-DGGE

    10.16258/j.cnki.1674-5906.2015.08.005

    Q948;X171.1

    A

    1674-5906(2015)08-1287-08

    張明錦,張健,紀托未,劉華,李勛,張艷,楊萬勤,陳良華. 林窗對凋落物分解過程中細菌群落結(jié)構(gòu)和多樣性的影響[J]. 生態(tài)環(huán)境學報, 2015, 24(8): 1287-1294.

    ZHANG Mingjin, ZHANG Jian, JI Tuowei, LIU Hua, LI Xun, ZHANG Yan, YANG Wanqin, CHEN Lianghua. Influence of Forest Gap on Bacterial Community Structure and Diversity during Litter Decomposition [J]. Ecology and Environmental Sciences, 2015, 24(8): 1287-1294.

    國家自然科學基金項目(31370628);國家科技支撐計劃項目(2011BAC09B05);四川省科技支撐計劃項目(12ZC0017);四川省科技廳應用基礎項目(2012JY0047);四川省教育廳科技創(chuàng)新團隊計劃項目(11TD006)

    張明錦(1990年生),碩士研究生,主要從事退化人工林生態(tài)恢復與重建研究。E-mail: 958192283@qq.com *通信作者:陳良華(1982年生),副研究員,研究生導師,主要從事植物生理生態(tài),生態(tài)環(huán)境修復等研究。E-mail: chenlh@sicau.edu.cn

    2015-06-01

    猜你喜歡
    林窗類群均勻度
    低播量下雜交稻產(chǎn)量形成對種植均勻度的響應
    作物研究(2023年2期)2023-05-28 13:44:14
    木荷林窗微環(huán)境特征及幼苗更新研究
    亞熱帶闊葉林林窗干擾特征與分布模式研究
    均勻度控制不佳可致肉種雞晚產(chǎn)
    薏苡種質(zhì)資源ISSR分子標記篩選及親緣關(guān)系分析
    天然林林窗與自然更新研究進展*
    黑石頂自然保護區(qū)土壤動物功能類群對季節(jié)動態(tài)的響應
    錦綸長絲染色均勻度判色新方法
    2011年春夏季黃海和東海微型浮游動物類群組成及其攝食的研究
    杉木人工林林窗特征因子與閩楠幼林生長關(guān)系的研究
    侵犯人妻中文字幕一二三四区| 亚洲精品中文字幕在线视频| 久久精品国产99精品国产亚洲性色 | 国产精品一区二区免费欧美| 麻豆av在线久日| 国产欧美日韩一区二区三| 亚洲 国产 在线| 飞空精品影院首页| 正在播放国产对白刺激| 90打野战视频偷拍视频| 一进一出好大好爽视频| 嫩草影视91久久| 9191精品国产免费久久| 久久香蕉激情| 少妇被粗大的猛进出69影院| 久久精品熟女亚洲av麻豆精品| 国产高清视频在线播放一区| 岛国在线观看网站| 老司机午夜福利在线观看视频 | 汤姆久久久久久久影院中文字幕| 新久久久久国产一级毛片| 动漫黄色视频在线观看| 麻豆成人av在线观看| 在线观看66精品国产| 亚洲欧美日韩高清在线视频 | 日韩制服丝袜自拍偷拍| 欧美性长视频在线观看| 丰满迷人的少妇在线观看| 操美女的视频在线观看| 国产熟女午夜一区二区三区| 如日韩欧美国产精品一区二区三区| 国产精品九九99| 丝瓜视频免费看黄片| 99国产精品一区二区三区| 久久 成人 亚洲| 成年女人毛片免费观看观看9 | 日韩免费av在线播放| 精品福利观看| 亚洲午夜理论影院| 国产日韩一区二区三区精品不卡| 别揉我奶头~嗯~啊~动态视频| 亚洲男人天堂网一区| 久久中文看片网| 国产视频一区二区在线看| 久久 成人 亚洲| 久久久国产成人免费| 黄频高清免费视频| 午夜福利,免费看| 91成人精品电影| 午夜两性在线视频| 欧美乱码精品一区二区三区| 亚洲欧美日韩高清在线视频 | 国产高清视频在线播放一区| 老司机亚洲免费影院| 高潮久久久久久久久久久不卡| e午夜精品久久久久久久| 两性夫妻黄色片| 黄色片一级片一级黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 1024视频免费在线观看| 一个人免费在线观看的高清视频| 日韩 欧美 亚洲 中文字幕| 999精品在线视频| 日本黄色视频三级网站网址 | 国产成人免费观看mmmm| 成人黄色视频免费在线看| 国产精品自产拍在线观看55亚洲 | 久久久久视频综合| 亚洲熟女精品中文字幕| 日韩有码中文字幕| 亚洲伊人久久精品综合| 日本av免费视频播放| 香蕉丝袜av| 免费在线观看完整版高清| 午夜福利一区二区在线看| 丝袜美腿诱惑在线| 成年人黄色毛片网站| 色婷婷av一区二区三区视频| 精品久久久久久电影网| 可以免费在线观看a视频的电影网站| 精品久久蜜臀av无| 欧美成人免费av一区二区三区 | 欧美变态另类bdsm刘玥| 99国产极品粉嫩在线观看| 国产高清视频在线播放一区| 少妇粗大呻吟视频| 国产日韩欧美亚洲二区| 亚洲第一av免费看| 欧美变态另类bdsm刘玥| 免费看a级黄色片| 日韩成人在线观看一区二区三区| 国产不卡一卡二| 少妇精品久久久久久久| 免费看a级黄色片| 久久久国产欧美日韩av| 少妇被粗大的猛进出69影院| 日本黄色日本黄色录像| 一二三四社区在线视频社区8| 涩涩av久久男人的天堂| 飞空精品影院首页| 不卡一级毛片| 涩涩av久久男人的天堂| 两个人免费观看高清视频| 十八禁高潮呻吟视频| 久久久久久久国产电影| 欧美成人午夜精品| 九色亚洲精品在线播放| 欧美精品亚洲一区二区| 曰老女人黄片| www日本在线高清视频| 黄色a级毛片大全视频| 欧美久久黑人一区二区| 国产成人免费无遮挡视频| 伦理电影免费视频| 热99re8久久精品国产| 波多野结衣一区麻豆| 久久精品91无色码中文字幕| 18禁美女被吸乳视频| 日韩视频一区二区在线观看| 日本黄色日本黄色录像| 欧美激情 高清一区二区三区| 国产av一区二区精品久久| 巨乳人妻的诱惑在线观看| 97在线人人人人妻| 91字幕亚洲| 国产成人欧美| 天堂俺去俺来也www色官网| 侵犯人妻中文字幕一二三四区| 亚洲一区中文字幕在线| 欧美国产精品va在线观看不卡| 精品久久久精品久久久| 亚洲av成人一区二区三| 日日夜夜操网爽| 国产精品久久久久久精品电影小说| 另类精品久久| 男人舔女人的私密视频| 日韩欧美一区视频在线观看| 中文字幕人妻丝袜制服| 国产精品.久久久| 亚洲精品美女久久久久99蜜臀| 精品乱码久久久久久99久播| 男女床上黄色一级片免费看| 高清视频免费观看一区二区| 美女高潮喷水抽搐中文字幕| 精品国产一区二区三区久久久樱花| 久久精品成人免费网站| 18禁黄网站禁片午夜丰满| 99久久精品国产亚洲精品| 欧美成狂野欧美在线观看| 国产精品偷伦视频观看了| a级毛片在线看网站| 国产麻豆69| 夫妻午夜视频| 99精品在免费线老司机午夜| 亚洲色图综合在线观看| 我要看黄色一级片免费的| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 国产又色又爽无遮挡免费看| 免费在线观看完整版高清| 99国产精品免费福利视频| 啦啦啦免费观看视频1| videos熟女内射| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| 多毛熟女@视频| 成人三级做爰电影| 日韩大片免费观看网站| 人妻久久中文字幕网| 色尼玛亚洲综合影院| 国产aⅴ精品一区二区三区波| 久久午夜综合久久蜜桃| 美国免费a级毛片| 丝袜在线中文字幕| 久久精品熟女亚洲av麻豆精品| 日本欧美视频一区| 国产精品国产av在线观看| 大陆偷拍与自拍| 国产成人精品无人区| 久久天躁狠狠躁夜夜2o2o| 免费在线观看黄色视频的| 色综合欧美亚洲国产小说| 亚洲精品久久午夜乱码| 中文字幕高清在线视频| 脱女人内裤的视频| 99国产精品99久久久久| 在线十欧美十亚洲十日本专区| 精品国产超薄肉色丝袜足j| 巨乳人妻的诱惑在线观看| 亚洲伊人色综图| 男女无遮挡免费网站观看| 蜜桃在线观看..| 19禁男女啪啪无遮挡网站| 啦啦啦 在线观看视频| 久久中文字幕一级| 亚洲七黄色美女视频| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 色尼玛亚洲综合影院| 搡老熟女国产l中国老女人| 成人精品一区二区免费| 亚洲 欧美一区二区三区| 亚洲av片天天在线观看| 午夜日韩欧美国产| netflix在线观看网站| 国产男女内射视频| 最新在线观看一区二区三区| 在线观看免费午夜福利视频| 亚洲精品国产一区二区精华液| netflix在线观看网站| 亚洲国产欧美日韩在线播放| 午夜视频精品福利| 免费看a级黄色片| 国产主播在线观看一区二区| 啦啦啦 在线观看视频| 丝袜人妻中文字幕| 操美女的视频在线观看| 青青草视频在线视频观看| 精品一区二区三区av网在线观看 | 国产精品98久久久久久宅男小说| 午夜免费成人在线视频| 欧美老熟妇乱子伦牲交| 国产精品自产拍在线观看55亚洲 | 国产精品香港三级国产av潘金莲| 精品人妻在线不人妻| 亚洲性夜色夜夜综合| 国产av又大| 黑人猛操日本美女一级片| 9191精品国产免费久久| 国产精品国产高清国产av | 久久精品亚洲熟妇少妇任你| 精品一区二区三区视频在线观看免费 | 91九色精品人成在线观看| 亚洲五月色婷婷综合| 免费黄频网站在线观看国产| 老汉色av国产亚洲站长工具| 黄频高清免费视频| 嫁个100分男人电影在线观看| 成人国语在线视频| 我的亚洲天堂| 我要看黄色一级片免费的| 午夜精品国产一区二区电影| 人人妻人人爽人人添夜夜欢视频| 久久精品91无色码中文字幕| 国产单亲对白刺激| 国产在视频线精品| 久久国产精品影院| 日韩制服丝袜自拍偷拍| 亚洲国产看品久久| av有码第一页| 99riav亚洲国产免费| 久久久久国产一级毛片高清牌| 久久人人爽av亚洲精品天堂| av不卡在线播放| 另类精品久久| 免费少妇av软件| 亚洲中文av在线| 视频在线观看一区二区三区| 色婷婷av一区二区三区视频| 最新美女视频免费是黄的| 久久精品亚洲精品国产色婷小说| 亚洲欧美激情在线| 久久精品aⅴ一区二区三区四区| 免费观看av网站的网址| a在线观看视频网站| 午夜老司机福利片| 中文字幕最新亚洲高清| 午夜福利视频精品| 久热爱精品视频在线9| 搡老乐熟女国产| 亚洲国产欧美网| 久热爱精品视频在线9| 精品人妻1区二区| 国产一卡二卡三卡精品| 久久中文字幕人妻熟女| 国产精品成人在线| 成年女人毛片免费观看观看9 | kizo精华| 国产精品av久久久久免费| 成年人黄色毛片网站| 国产精品欧美亚洲77777| 国产成+人综合+亚洲专区| 日韩制服丝袜自拍偷拍| 99国产精品免费福利视频| 久久精品成人免费网站| 久久ye,这里只有精品| 热99国产精品久久久久久7| 国产日韩欧美视频二区| 国产日韩欧美在线精品| 女性生殖器流出的白浆| 欧美日韩一级在线毛片| 亚洲精品av麻豆狂野| 欧美亚洲日本最大视频资源| 天堂动漫精品| 亚洲一区二区三区欧美精品| 国产精品久久电影中文字幕 | 国产成人av教育| 久久 成人 亚洲| 久久婷婷成人综合色麻豆| 蜜桃在线观看..| 亚洲性夜色夜夜综合| 成人黄色视频免费在线看| 欧美大码av| 亚洲欧美精品综合一区二区三区| 久久精品aⅴ一区二区三区四区| 国产精品久久久久久人妻精品电影 | 精品久久蜜臀av无| 一区二区三区国产精品乱码| 一区二区三区激情视频| 精品人妻在线不人妻| 亚洲av国产av综合av卡| 免费av中文字幕在线| 免费在线观看视频国产中文字幕亚洲| 在线观看人妻少妇| 母亲3免费完整高清在线观看| 美女扒开内裤让男人捅视频| 两性午夜刺激爽爽歪歪视频在线观看 | 精品久久久久久久毛片微露脸| 国产有黄有色有爽视频| 亚洲国产欧美一区二区综合| 国产精品美女特级片免费视频播放器 | 99久久99久久久精品蜜桃| 老熟妇仑乱视频hdxx| 欧美激情极品国产一区二区三区| 久久久国产成人免费| 天天操日日干夜夜撸| 首页视频小说图片口味搜索| 亚洲av成人一区二区三| 亚洲全国av大片| 亚洲精品国产精品久久久不卡| 男人操女人黄网站| 最近最新中文字幕大全电影3 | 日韩欧美一区视频在线观看| 久久天躁狠狠躁夜夜2o2o| 熟女少妇亚洲综合色aaa.| 夜夜骑夜夜射夜夜干| 亚洲精品久久成人aⅴ小说| 久久精品亚洲av国产电影网| xxxhd国产人妻xxx| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看 | 国产成人免费无遮挡视频| 精品一区二区三区视频在线观看免费 | 国产男靠女视频免费网站| 国产一区二区激情短视频| 免费观看av网站的网址| aaaaa片日本免费| 一级毛片精品| 久久久久久久大尺度免费视频| 国产精品九九99| 男人操女人黄网站| 亚洲 欧美一区二区三区| 久久精品国产a三级三级三级| 国内毛片毛片毛片毛片毛片| 国产成人免费观看mmmm| 色老头精品视频在线观看| 亚洲第一av免费看| 大陆偷拍与自拍| 亚洲国产av新网站| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 亚洲性夜色夜夜综合| 首页视频小说图片口味搜索| 国产精品1区2区在线观看. | 我要看黄色一级片免费的| 日韩欧美一区视频在线观看| 亚洲天堂av无毛| 高清欧美精品videossex| 在线永久观看黄色视频| 欧美人与性动交α欧美精品济南到| 91九色精品人成在线观看| 最新在线观看一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美亚洲二区| 极品教师在线免费播放| 色婷婷av一区二区三区视频| av电影中文网址| 国产又色又爽无遮挡免费看| 亚洲精品自拍成人| 国产在线精品亚洲第一网站| 大型黄色视频在线免费观看| 中文字幕人妻丝袜制服| 婷婷丁香在线五月| 97人妻天天添夜夜摸| 亚洲视频免费观看视频| 久久九九热精品免费| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 美女高潮喷水抽搐中文字幕| 精品久久久久久电影网| 国产精品99久久99久久久不卡| 久久久久精品人妻al黑| 国产成人欧美在线观看 | 国产1区2区3区精品| 在线观看人妻少妇| 久久精品成人免费网站| 国精品久久久久久国模美| 乱人伦中国视频| 日韩欧美免费精品| 久久影院123| 国产av精品麻豆| a在线观看视频网站| 国产精品 欧美亚洲| 女同久久另类99精品国产91| 五月开心婷婷网| 人妻久久中文字幕网| 欧美成人午夜精品| 国产午夜精品久久久久久| 最新在线观看一区二区三区| 午夜日韩欧美国产| 免费少妇av软件| 免费看十八禁软件| 日日夜夜操网爽| 亚洲美女黄片视频| a级毛片在线看网站| 我的亚洲天堂| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| 欧美国产精品一级二级三级| 另类精品久久| 黄色视频,在线免费观看| 国产片内射在线| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 十八禁网站网址无遮挡| 亚洲人成电影免费在线| netflix在线观看网站| 亚洲国产欧美日韩在线播放| 国产一区有黄有色的免费视频| 久久久久久久精品吃奶| 亚洲欧洲精品一区二区精品久久久| 成人特级黄色片久久久久久久 | 男女之事视频高清在线观看| 国产精品欧美亚洲77777| 露出奶头的视频| 岛国毛片在线播放| 亚洲综合色网址| 久热爱精品视频在线9| 久久久久精品人妻al黑| 大片电影免费在线观看免费| 纯流量卡能插随身wifi吗| 桃红色精品国产亚洲av| 九色亚洲精品在线播放| 女性生殖器流出的白浆| 99精品在免费线老司机午夜| videos熟女内射| 久久午夜亚洲精品久久| 成年人午夜在线观看视频| 极品教师在线免费播放| 黑人欧美特级aaaaaa片| 在线观看人妻少妇| 日本av免费视频播放| 免费日韩欧美在线观看| 捣出白浆h1v1| 夫妻午夜视频| videos熟女内射| 国产成人一区二区三区免费视频网站| 精品国产超薄肉色丝袜足j| 亚洲国产精品一区二区三区在线| 黄色怎么调成土黄色| 久久亚洲精品不卡| 免费不卡黄色视频| 国产视频一区二区在线看| 90打野战视频偷拍视频| www.精华液| 精品亚洲成国产av| 在线观看免费午夜福利视频| 一二三四在线观看免费中文在| 99精国产麻豆久久婷婷| 亚洲欧美日韩高清在线视频 | 久久久精品国产亚洲av高清涩受| 一级a爱视频在线免费观看| 成在线人永久免费视频| 91老司机精品| 久久精品熟女亚洲av麻豆精品| 亚洲国产av新网站| 亚洲三区欧美一区| 香蕉丝袜av| 亚洲中文日韩欧美视频| 午夜福利,免费看| av又黄又爽大尺度在线免费看| 天天躁夜夜躁狠狠躁躁| 国产成人影院久久av| 欧美人与性动交α欧美软件| 日韩三级视频一区二区三区| 久久久久精品国产欧美久久久| 亚洲少妇的诱惑av| 后天国语完整版免费观看| 黄色视频在线播放观看不卡| 亚洲性夜色夜夜综合| 丁香六月欧美| 美女高潮喷水抽搐中文字幕| 亚洲少妇的诱惑av| 亚洲色图 男人天堂 中文字幕| 成年人午夜在线观看视频| 欧美日韩中文字幕国产精品一区二区三区 | 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 纯流量卡能插随身wifi吗| 免费看十八禁软件| 欧美日韩国产mv在线观看视频| 18禁国产床啪视频网站| 中文字幕色久视频| 99久久人妻综合| 成人黄色视频免费在线看| 一级黄色大片毛片| 亚洲av第一区精品v没综合| 色综合婷婷激情| 久久久水蜜桃国产精品网| 久久久国产一区二区| 欧美+亚洲+日韩+国产| 超色免费av| 老熟妇乱子伦视频在线观看| 午夜福利在线免费观看网站| 精品免费久久久久久久清纯 | 日本wwww免费看| 国产精品久久久久久精品古装| 欧美精品人与动牲交sv欧美| av不卡在线播放| 欧美激情高清一区二区三区| 亚洲成人手机| 99在线人妻在线中文字幕 | 日本撒尿小便嘘嘘汇集6| 不卡一级毛片| 亚洲国产欧美一区二区综合| 黑人操中国人逼视频| 精品少妇黑人巨大在线播放| 夫妻午夜视频| 国产一区有黄有色的免费视频| 精品国产乱码久久久久久小说| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 中文字幕高清在线视频| 久久中文看片网| 亚洲国产中文字幕在线视频| 精品亚洲成a人片在线观看| tocl精华| 一边摸一边做爽爽视频免费| 最近最新中文字幕大全免费视频| 亚洲av片天天在线观看| 91精品三级在线观看| av一本久久久久| 日韩 欧美 亚洲 中文字幕| 日韩一卡2卡3卡4卡2021年| 欧美 日韩 精品 国产| a在线观看视频网站| 欧美日韩av久久| 日本a在线网址| 丰满少妇做爰视频| 欧美日韩中文字幕国产精品一区二区三区 | 巨乳人妻的诱惑在线观看| 久久久精品国产亚洲av高清涩受| 国产一区二区三区视频了| 国产成人影院久久av| 国产精品一区二区精品视频观看| 一夜夜www| 香蕉久久夜色| videosex国产| 99精品在免费线老司机午夜| 亚洲精品美女久久av网站| 50天的宝宝边吃奶边哭怎么回事| 亚洲,欧美精品.| 制服人妻中文乱码| 人妻一区二区av| 精品一区二区三区av网在线观看 | 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 成人亚洲精品一区在线观看| 国精品久久久久久国模美| 99在线人妻在线中文字幕 | 99精品在免费线老司机午夜| 国产午夜精品久久久久久| 午夜福利视频精品| 露出奶头的视频| 天天躁日日躁夜夜躁夜夜| 国内毛片毛片毛片毛片毛片| 黑人操中国人逼视频| 午夜福利在线免费观看网站| 欧美乱码精品一区二区三区| 日韩欧美免费精品| 女性被躁到高潮视频| 欧美在线一区亚洲| 亚洲精品一卡2卡三卡4卡5卡| 麻豆av在线久日| 黄色成人免费大全| 高清毛片免费观看视频网站 | 亚洲成人免费av在线播放| 亚洲国产av新网站| 丝瓜视频免费看黄片| 极品人妻少妇av视频| 新久久久久国产一级毛片| 国产97色在线日韩免费| 丝袜喷水一区| 老司机福利观看| 久久精品国产亚洲av香蕉五月 | 久久久久视频综合| 国产高清激情床上av| 制服人妻中文乱码| 考比视频在线观看| av天堂久久9| 国产精品九九99| 大片免费播放器 马上看| 国产精品av久久久久免费| 99国产精品一区二区三区| cao死你这个sao货| 女同久久另类99精品国产91| 丝袜人妻中文字幕| 欧美亚洲日本最大视频资源| 99精品久久久久人妻精品| 午夜免费成人在线视频| 欧美成人免费av一区二区三区 | 一边摸一边抽搐一进一出视频| 麻豆av在线久日| 久久久久久免费高清国产稀缺| av福利片在线|