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

    煙稈生物炭還田對植煙黃棕壤養(yǎng)分、細菌群落結構與功能的影響

    2023-07-22 14:56:13戴華鑫梁太波李彩斌羅貞寶關羅浩鐘帥翟振邊文杰黃衛(wèi)張艷玲
    南方農(nóng)業(yè)學報 2023年2期
    關鍵詞:生物炭土壤養(yǎng)分

    戴華鑫 梁太波 李彩斌 羅貞寶 關羅浩 鐘帥 翟振 邊文杰 黃衛(wèi) 張艷玲

    摘要:【目的】明確煙稈生物炭對植煙黃棕壤養(yǎng)分、細菌群落結構與功能的影響,為提高煙稈的資源化利用提供理論依據(jù)。【方法】基于大田試驗,比較不同用量煙稈生物炭(0、5、15、20和40 t/ha)施用4年后煙田土壤養(yǎng)分指標及細菌豐度變化;采用16S rDNA序列高通量測序技術和PICRUSt基因預測方法,分析耕層土壤細菌群落結構、多樣性及代謝功能變化?!窘Y果】隨著生物炭施用量的增加,土壤pH、有機質、總氮、有效磷、速效鉀含量及細菌豐度呈增加趨勢。主坐標分析(PCoA)和聚類分析結果表明,施用生物炭改變了細菌群落組成。生物炭用量在15 t/ha時,土壤細菌OTU豐富度及Alpha多樣性指數(shù)最高;在門水平上,變形菌門(Proteobacteria)、放線菌門(Actinobacteria)相對豐度較低,而芽單胞菌門(Gemmatimonadetes)、綠彎菌門(Chloroflexi)、浮霉菌門(Planctomycetes)和硝化螺旋菌門(Nitrospirae)相對豐度較高;在屬水平上,芽單胞菌屬(Gemmatimonas)、黃色土源菌屬(Flavisolibacter)和Roseiflexus相對豐度較高。冗余分析(RDA)結果表明,生物炭用量與土壤擬桿菌門顯著正相關(P<0.05,下同),與酸桿菌門顯著負相關。PICRUSt功能預測表明,施用生物炭增加了與土壤碳代謝相關的細菌基因豐度。【結論】貴州植煙黃棕壤施用煙稈生物炭后,土壤養(yǎng)分、烤煙農(nóng)藝性狀和經(jīng)濟性狀指標均有提升,土壤細菌群落豐富度及多樣性顯著提高,土壤細菌碳代謝能力增強,對改善煙田土壤微生態(tài)具有積極作用。以15 t/ha煙稈生物炭作用效果較好,可在生產(chǎn)上推廣應用。

    關鍵詞:生物炭;植煙黃棕壤;土壤養(yǎng)分;細菌群落;功能預測

    中圖分類號:S572.01? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼:A 文章編號:2095-1191(2023)02-0476-12

    Abstract:【Objective】In order to explore the effects of tobacco-stalk biochar (BC) on soil nutrients the structure and functioning of bacterial communities in tobacco-planting yellow-brown fields,and provide theoretical basis for the utilization of tobacco stalk. 【Method】A field experiment with different biochar dosages(0,5,15,20,40 t/ha) was set up to compare the changes of soil nutrient indexes and bacterial abundance of tobacco field after four years application. The 16S rDNA high-throughput sequencing and PICRUSt gene prediction analysis were used to decipher soil community composition, diversity and metabolic function. 【Result】The results showed that soil pH,organic matter,total nitrogen,available phosphorus,available potassium and the abundance of soil bacteria significantly increased with the dosage of biochar application. The results of principal coordinate analysis (PCoA) and cluster analysis show that soil bacterial community composition was significantly altered by the biochar application. When the amount of tobacco-stalk biochar was 15 t/ha,bacteria OTU richness and Alpha diversity indexes were higher than other treatments. At the phylum level,the relative abundance of Proteobacteria and Actinobacteria was lower under 15 t/ha than other treatments,but those of Gemmatimona-detes,Chloroflexi,Planctomycetes and Nitrospirae were high. At the genus level,the relative abundance of Gemmatimonas,F(xiàn)lavisolibacter,Roseiflexus was higher under 15 t/ha than other treatments. Redundancy analysis(RDA) analysis showed that biochar application dosage was significantly positively correlated with the relative abundance of Bacteroidetes(P<0.05, the same below),while negatively correlated with Acidobacteria. PICRUSt functional prediction analysis suggested that biochar application increased the functional gene abundance related to soil carbon metabolism in bacterial community. 【Conclusion】After applying tobacco-stalk biochar to the yellow-brown tobacco-planting soil in Guizhou, soil nu-trients, agronomic and economic indicators of flue-cured tobacco are improved, and the richness and diversity of soil bacterial communities are significantly increased. The carbon metabolism ability of soil bacteria is enhanced, which has a positive effect on improving the soil microecology of tobacco fields. In sum, the effect of using 15 t/ha tobacco-stalk biochar isrelatively better and can be popularized in tobacco production.

    Key words: biochar; tobacco-planting yellow-brown soil; soil nutrients; bacterial community; function prediction

    Foundation items: Key Scientific and Technological Project of China National Tobacco Corporation (110201902004);Science and Technology Project of China Tobacco Guangdong Industrial Co., Ltd.(QJ/GY G03.009/02);Youth Talent Support Project of Zhengzhou Tobacco Research Institute(112020CR0250)

    0 引言

    【研究意義】生物炭是農(nóng)業(yè)廢棄有機物等經(jīng)高溫低氧或限氧環(huán)境熱解后產(chǎn)生的一種富碳物質(Blanco-Canqui,2017)。生物炭具有豐度的微孔結構、較大的比表面積、較高的穩(wěn)定性和吸附能力,在改良土壤、固碳減排、修復環(huán)境等方面具有很大的應用前景(Oni et al.,2019)。煙草是我國主要的經(jīng)濟作物之一,我國煙區(qū)每年生產(chǎn)近450萬t煙草秸稈(李文豪,2019),主要處理方式是直接燃燒和秸稈還田。秸稈燃燒在造成環(huán)境污染的同時也浪費了生物質資源,直接還田則可能引起病蟲害爆發(fā)進而造成煙草產(chǎn)質量下滑。有研究報道,合理施用生物炭可提高土壤保水性能(Wang et al.,2020;曾鳳玲等,2022),增加土壤碳氮儲量和有效養(yǎng)分含量(黃雁飛等,2020;張珂珂等,2022),修復土壤重金屬污染(白珊等,2021),改良作物根系發(fā)育狀況(陳懿等,2017),有利于作物生長和品質提升(李青山等,2021;李松偉等,2021;王志丹等,2022)。貴州是我國第二大煙葉產(chǎn)區(qū),黃棕壤是當?shù)乜緹煼N植的主要土壤類型之一。因此,研究煙稈生物炭還田對植煙黃棕壤養(yǎng)分、細菌群落結構與功能的影響,對改善當?shù)刂矡熗寥佬再|、煙稈的合理利用及生態(tài)保護均具有重要意義。【前人研究進展】土壤微生物是土壤碳氮循環(huán)的“推進器”和土壤養(yǎng)分有效性的“轉換器”,由于土壤微生物對環(huán)境變化的敏感性,其自身種群多樣性及優(yōu)勢菌屬的變化可在一定程度上指示土壤肥力和質量狀況的高低(Luo et al.,2016;桂意云等,2022)。王光飛等(2017)通過盆栽試驗研究秸稈生物炭(0~2%)對種植辣椒土壤性狀的影響,結果發(fā)現(xiàn)土壤pH、有機質、有效磷和速效鉀含量隨著生物炭用量的增加逐漸上升,土壤細菌及4種功能菌數(shù)量也逐步增加,微生物代謝能力、多樣性及均勻度均與生物炭用量呈倒U型曲線關系,其在1.33%用量時達最高值。況帥等(2021)在田間植煙紅壤中添加油菜秸稈生物炭(0~9 t/ha),發(fā)現(xiàn)土壤養(yǎng)分、細菌OTU數(shù)目、菌群豐度和多樣性等明顯增加,土壤中主要優(yōu)勢菌群變形菌門(Proteobacteria)、放線菌門(Actinobacteria)與酸桿菌門(Acidobacteria)與土壤有效磷和速效鉀含量呈顯著正相關。姚麗茹等(2022)在半干旱農(nóng)業(yè)地區(qū)通過長期田間定位試驗研究施用小麥秸稈生物炭(0~20 t/ha)后作物生長及土壤細菌群落的變化,結果顯示隨著生物炭施用量的增加,土壤pH、有機碳、全氮、小麥生物量和產(chǎn)量呈上升趨勢,變形菌門(Proteobacteria)和浮霉菌門(Planctomycetes)相對豐度增加,土壤細菌群落變化主要受生物炭施用量和采樣時期的影響,以10~20 t/ha生物炭用量較合適。然而,Zhang等(2021)在東北黑土單次施用10~50 t/ha玉米秸稈生物炭3年后,發(fā)現(xiàn)土壤細菌香農(nóng)指數(shù)顯著降低。由此看出,生物炭對土壤細菌群落結構的影響并不一致,這可能與土壤類型、利用方式及生物炭種類、用量等因素有關。在烤煙生產(chǎn)中,葉協(xié)鋒等(2015)研究指出,適量添加生物炭可改善烤后煙葉單料煙評吸質量,以常規(guī)施肥配施600 kg/ha生物炭用量的處理效果最好,而過量的生物炭可能會對烤后煙葉品質造成負面影響;鄧建強等(2018)研究指出,適量生物炭可降低土壤可能致病菌屬的相對豐度,但高劑量施用將增大煙葉生產(chǎn)感病風險?!颈狙芯壳腥朦c】前人關于生物炭的研究多集中于施用當年單一階段的變化,而關于生物炭施用多年后土壤養(yǎng)分和細菌群落變化的研究相對較少,尤其是針對貴州植煙黃棕壤施用生物炭后土壤養(yǎng)分和細菌群落變化的研究尚無報道?!緮M解決的關鍵問題】通過大田試驗,重點研究煙稈生物炭施用4年后對植煙土壤養(yǎng)分、細菌群落結構及功能的影響,旨在闡明生物炭對貴州典型植煙黃棕壤微生態(tài)的調控效應,為提高烤煙秸稈的資源化利用提供理論依據(jù)。

    1 材料與方法

    1. 1 試驗材料

    生物炭(煙稈)材料來自貴州金葉豐農(nóng)業(yè)科技有限公司。煙稈在380 ℃下裂解2 h制備成煙稈生物炭,其基本性質為:pH 9.18、碳48.3%、氫3.93%、氧31.4%、氮1.50%、全磷2.38%、全鉀2.44%、速效鉀161.9 mg/kg、比表面積1.47 m2/g、總孔隙體積0.0075 cm3/g、平均孔直徑20.31 nm。供試煙草品種為云煙87。土壤類型為黃棕壤,2018年3月下旬拔除前茬煙稈和雜草后采集基礎土樣,檢測土壤理化指標:pH 5.40、有機質21.0 g/kg、總氮0.96 g/kg、全磷1.00 g/kg、全鉀1.71%、堿解氮102.90 mg/kg、有效磷28.30 mg/kg、速效鉀196.40 mg/kg。

    1. 2 試驗方法

    試驗于2018—2021年在貴州威寧黑石鎮(zhèn)科技園(東經(jīng)103°59′55″、北緯26°45′14″)進行。2018年5月,將煙稈生物炭一次性施入大田,考慮生物炭實際成本及田間承載狀況,用量分別為0 t/ha(CK)、5 t/ha(T1)、15 t/ha(T2)、20 t/ha(T3)和40 t/ha(T4),之后不再施入生物炭;生物炭在起壟前均勻撒施在土壤中,隨后將其與耕層土壤迅速旋耕混勻。施肥用量參照當?shù)爻R?guī)優(yōu)質烤煙生產(chǎn),每公頃施用酒糟有機肥1800.0 kg、復合肥(9-13-22)525.0 kg、提苗肥(15-8-7)37.5 kg和追肥(13-0-26)300.0 kg。隨機區(qū)組設計,3次重復,共15個小區(qū),株距0.55 m,行距1.10 m,小區(qū)面積74.8 m2(長11.0 m、寬6.8 m),試驗地面積1122.0 m2。基肥在移栽前結合有機肥窩施,提苗肥于移栽時結合定根水施用(肥料融入水中澆灌)、移栽25 d后穴施追肥,其他各項田間生產(chǎn)管理措施統(tǒng)一按當?shù)丶夹g規(guī)范進行。

    連續(xù)4年種植烤煙,2021年8月在圓頂期測量煙株農(nóng)藝性狀,9月大田煙葉采收結束后,在每小區(qū)選取8株煙株,采集每2株煙株中心點處的耕層(0~20 cm)土壤,均勻混合成1個樣品。每個樣品分成2份,1份放入冰盒運回實驗室,用于土壤細菌豐度及多樣性檢測;另1份裝于布袋中,帶回實驗室自然風干,過篩后用于土壤養(yǎng)分指標測定。采收煙葉統(tǒng)一裝于烤房中部,按當?shù)爻R?guī)烘烤工藝烘烤和分級收購,統(tǒng)計煙葉經(jīng)濟性狀。

    1. 3 測定指標及方法

    1. 3. 1 土壤養(yǎng)分指標測定 參照《土壤農(nóng)業(yè)化學分析方法》(魯如坤,2000)測定土壤養(yǎng)分指標。土壤pH采用水浸提電位法(水土比2.5∶1)測定;土壤有機質含量采用重鉻酸鉀外加熱法測定;土壤總氮和堿解氮含量分別采用半微量凱氏定氮法和堿解擴散法測定;土壤有效磷及速效鉀含量分別采用用氟化銨—鹽酸浸提—鉬銻抗比色法及醋酸銨浸提—火焰光度法測定。

    1. 3. 2 土壤細菌DNA提取與擴增 采用土壤DNA提取試劑盒(PowerSoil DNA Isolation Kit,MOBIO,USA)對土壤樣本提取基因組DNA,提取后的DNA利用瓊脂糖凝膠電泳進行檢驗,檢驗合格后,將DNA樣本放入-20 ℃條件保存,后續(xù)用于qPCR檢測及二代高通量測序。qPCR檢測和二代高通量測序中均采用16S rRNA基因V4~V5區(qū)域通用引物,引物序列為515F(5'-GTGCCAGCMGCCGCGGTAA-3')和926R(5'-CCGTCAATTCMTTTGAGTTT-3'),由微基生物科技(上海)有限公司完成測序及分析。

    1. 3. 3 土壤細菌豐度的qPCR檢測 運用染料法熒光定量PCR檢測土壤樣本中細菌16S rRNA基因V4~V5區(qū)的拷貝數(shù)。采用16S rRNA基因V4~V5區(qū)通用引物擴增基因組DNA,采用瓊脂糖凝膠回收試劑盒(Axygen,USA)回收PCR產(chǎn)物,隨后進行TA克隆,篩選陽性克隆進行測序,并采用質粒提取試劑盒(Axygen,USA)進行質粒提取。采用Qubit 3.0(Invitrogen,USA)測定質粒DNA濃度,根據(jù)質粒濃度和質粒DNA相對分子量計算標準質粒拷貝數(shù)。將構建好的標準質粒進行10倍梯度稀釋,作為模板,構建標準曲線。qPCR反應總體系20.0 μL:2×SYBR Premix Ex TaqTM10.0 μL,上、下游引物(10 μmol/L)各0.4 μL,ROX 0.4 μL,DNA模板3.0 μL,ddH2O補足至20.0 μL。每個反應3個平行,陰性對照采用ddH2O作為模板。PCR擴增程序:95 ℃預變性30 s;95 ℃ 10 s,55 ℃ 30 s,72 ℃ 30 s,進行40次循環(huán)。熒光數(shù)據(jù)采集在延伸時進行。實時定量PCR在ABI 7900HT型熒光定量PCR儀中進行。結果以每克樣本中16S V4~V5區(qū)基因拷貝數(shù)統(tǒng)計。

    1. 3. 4 土壤細菌多樣性的二代高通量測序 采用二代高通量測序技術檢測土壤樣本中細菌多樣性。采用兩步PCR擴增方法構建文庫(戴華鑫等,2017),第一步采用含有部分測序引物和16S rRNA基因V4~V5通用引物的融合引物進行PCR擴增目的片段;以回收的PCR產(chǎn)物為模板,進行第二步擴增,添加Barcode和測序接頭序列,并用FTC-3000TM Real-Time PCR儀進行熒光定量,均一化混勻后完成文庫構建,文庫質控合格后,采用Illumina NovaSeq 6000進行測序。

    1. 3. 5 PICRUSt2功能預測 基于OTU的代表序列和豐度信息,采用PICRUSt2進行功能和代謝途徑預測分析,通過比對KEGG數(shù)據(jù)庫獲得不同代謝通路的豐度信息(Langille et al.,2013)。

    1. 4 數(shù)據(jù)分析

    二代高通量原始測序結果經(jīng)Trimmomatic初步質控后,采用Flash進行序列拼接,采用UPARSE按照序列97%相似性進行OTU(Operational taxonomic unit)聚類,同時去掉singletons序列,每個OTU得到1條代表序列,采用mothur的classify.seqs命令與SILVA 128數(shù)據(jù)庫比對,進行物種信息注釋,閾值設置0.6。根據(jù)注釋結果,分別在門、屬分類水平上統(tǒng)計微生物相對豐度,并采用R語言ggplot2進行可視化(Gao et al.,2020)。利用R語言vegan包中的estimateR函數(shù)計算Chao和ACE豐富度指數(shù),diversity函數(shù)計算Shannon多樣性指數(shù),picante包的pd函數(shù)計算PD_whole_tree,rda函數(shù)進行冗余分析(RDA),采用ggplot2程序包對結果進行可視化。Beta多樣性分析,采用R語言ape程序包進行weighted.unifrac分析,dendextend包進行聚類樹分析,基于OTU數(shù)據(jù),采用bray curtis算法進行主坐標分析(PCoA),聚類樹形圖和PCoA圖采用ggplot2進行可視化。采用vegan包中的adonis函數(shù),對bray curtis結果進行adonis分析。ggpubr包進行Kruskal-Wallis秩和檢驗,對樣本間屬水平微生物豐度進行差異分析(Wang et al.,2019)。采用SPSS 24.0的單因素方差分析和Duncans新復極差法進行差異顯著性檢驗(P<0.05)。

    2 結果與分析

    2. 1 不同生物炭用量對土壤養(yǎng)分的影響

    由表1可知,隨著生物炭用量增加,土壤pH、有機質、總氮、有效磷和速效鉀含量均呈增加趨勢,堿解氮含量無顯著變化(P>0.05,下同)。T4處理的土壤有機質、速效鉀含量分別為49.53 g/kg、691.00 mg/kg,顯著高于其他處理(P<0.05,下同);土壤pH、總氮、有效磷含量高于CK和T1處理。由此可見,施用煙稈生物炭對土壤養(yǎng)分具有明顯的提升作用。

    2. 2 不同生物炭用量對烤煙農(nóng)藝性狀和經(jīng)濟性狀的影響

    由表2可知,施用生物炭后烤煙的莖圍、上部葉最大葉長、中部葉最大葉長顯著高于CK,株高、上部葉最大葉寬、中部葉最大葉寬、葉數(shù)等農(nóng)藝性狀與在各處理間均無顯著差異。從烤煙經(jīng)濟性狀來看,T1、T2處理的煙葉產(chǎn)量和產(chǎn)值顯著高于其他處理,不同處理間烤煙均價、上等煙和中等煙比例無顯著差異(表3)。

    2. 3 不同生物炭用量對土壤細菌豐度的影響

    由圖1可看出,土壤細菌豐度隨著生物炭用量的增加呈增加趨勢。與CK相比,增施5、15、20 t/ha生物炭后,T1、T2和T3處理的土壤細菌豐度增加不顯著;生物炭施用量增加至40 t/ha時,T4處理的土壤細菌豐度為1.82×108 CFU/g,顯著高于其他處理。

    2. 4 不同生物炭用量對土壤細菌Alpha多樣性的影響

    由表4可知,樣品序列覆蓋率均達0.99以上,說明文庫中包含細菌群落的絕大多數(shù)細菌類群,基本能反映群落結構組成。與其他處理相比,CK的土壤細菌OTUs和Alpha多樣性指數(shù)處于較低水平;隨著生物炭用量的增加,細菌OTUs和Alpha多樣性指數(shù)(Chao、ACE、Shannon、PD_whole_tree)均呈先升高后降低的變化趨勢,T2處理的土壤細菌群落Alpha多樣性指數(shù)均為最高。

    2. 5 不同生物炭用量對土壤細菌群落結構的影響

    由圖2可看出,土壤優(yōu)勢細菌相對豐度高于1.00%的細菌門有變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria)、芽單胞菌門(Gemmatimonadetes)、綠彎菌門(Chloroflexi)、擬桿菌門(Bacteroidetes)、浮霉菌門(Planctomycetes)、厚壁菌門(Firmicutes)和硝化螺旋菌門(Nitrospirae)。疣微菌門(Verrucomicrobia)和裝甲菌門(Armatimo-nadetes)的相對豐度在0.50%~1.00%,其他細菌門的相對豐度均低于0.50%。與其他處理相比,生物炭用量為15 t/ha時(T2處理)土壤變形菌門、放線菌門相對豐度較低,分別為25.70%和22.58%;芽單胞菌門、綠彎菌門、浮霉菌門、硝化螺旋菌門相對豐度較高,分別為12.35%、8.13%、3.34%、1.69%和0.63%。隨著生物炭用量的增加(CK~T4處理),酸桿菌門相對豐度呈明顯降低趨勢,分別為17.70%、14.39%、13.10%、12.54%和10.48%,擬桿菌門相對豐度呈明顯增加趨勢,分別為3.46%、5.14%、7.79%、8.39%和11.06%。

    由圖3可看出,生物炭用量對不同細菌屬相對豐度的影響存在差異。隨著生物炭用量的增加,多數(shù)細菌群落屬水平相對豐度總體呈下降趨勢。與其他處理相比,T4處理的Gemmatirosa、Bryobacter、水恒桿菌屬(Mizugakiibacter)、Gaiella、慢生根瘤菌屬(Bradyrhizobium)、念珠菌固體桿菌屬(Candidatus_Solibacter)、Reyranella、Jatrophihabitans、Dyella、馬杜拉放線菌屬(Actinomadura)、放線異壁酸菌屬(Actinoallomurus)、類芽孢桿菌屬(Paenibacillus)等細菌豐度較低;T2處理的芽單胞菌屬(Gemmatimonas)、黃色土源菌(Flavisolibacter)、Roseiflexus等細菌豐度相對較高;T3處理的放線菌屬(Actinophytocola)、Reyranella、類芽孢桿菌屬(Paenibacillus)等細菌豐度相對較高,即總體呈先升高后降低的拋物線趨勢;Acidibacter的相對豐度則呈先降低后升高的變化趨勢。

    基于OTU層次的非加權組平均法聚類分析結果表明,包括CK在內(nèi)的不同樣品可分為三大類;CK與T1處理距離相對較近,T2處理與T3處理距離相對較近,T4處理距離其他處理相對較遠(圖4)。該結果與試驗處理設置相一致,同時表明,與CK相比,生物炭施用4年后仍然改變著煙田耕層土壤的細菌群落結構。

    基于OTU水平和Bray-curtis距離的PCoA分析結果(圖5)表明,不同生物炭用量處理的土壤細菌群落組成出現(xiàn)明顯區(qū)分,前兩軸共解釋75.03%的群落變異,adonis結果也說明不同處理間存在極顯著差異(P<0.001)。CK樣品沿第1軸與T2、T3和T4處理的3組樣品明顯分開,T2、T3和T4處理樣品沿第2軸分開。

    2. 6 土壤細菌與環(huán)境因子的相關性

    細菌群落門水平與土壤環(huán)境因子的RDA分析結果(圖6)顯示,第1軸的解釋變異量為80.22%,第2軸的解釋變異量為14.68%,兩軸共解釋總變異的94.90%。在環(huán)境因子方面,生物炭與pH、有效磷、有機質、總氮、速效鉀呈銳角,呈正相關關系,與土壤堿解氮呈鈍角,呈負相關關系。細菌群落與土壤理化性質方面,擬桿菌門(Bacteroidetes)與土壤有機質、pH、生物炭、總氮、速效鉀呈正相關關系,與堿解氮呈負相關關系;酸桿菌門(Acidobacteria)、厚壁菌門(Firmicutes)均與土壤總氮、速效鉀、有效磷、pH、生物炭呈負相關關系,與堿解氮呈正相關關系;綠彎菌門(Chloroflexi)與土壤堿解氮呈負相關關系;變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、芽單胞菌門(Gemmatimonadetes)受生物炭影響相對較小。因此,生物炭用量與土壤擬桿菌門呈正相關關系,與酸桿菌門呈負相關關系。

    2. 7 土壤細菌群落PICRUSt功能預測分析

    將預測得到的功能基因家族(圖7)與KEGG數(shù)據(jù)庫比對發(fā)現(xiàn),對于第一等級,所有菌群的基因序列注釋到的預測功能可分為6類,其中前4類分別是:新陳代謝(Metabolism)、環(huán)境信息加工(Environmental Information Processing)、遺傳信息加工(Genetic Information Processing)和細胞過程(Cellular Processes)。對于第二等級,參與新陳代謝的菌群中,以參與碳水化合物代謝(Carbohydrate metabolism)、氨基酸代謝(Amino acid metabolism)、能量代謝(Energy meta-bolism)、輔助因子及維他命代謝(Metabolism of cofactors and vitamins)的為最大類群;參與遺傳信息加工的菌群中,參與遺傳信息翻譯(Translation)、復制和修復機制(Replication and repair)的類群具有明顯優(yōu)勢;參與環(huán)境信息加工的菌群中,參與膜運轉類群(Membrane transport)具有明顯優(yōu)勢,其次是參與信號轉導(Signal transduction)的類群;參與細胞過程的菌群中,以參與原核生物細胞群落代謝(Cellular community-prokaryotes)的為最大類群。

    2. 8 代謝相關基因分析

    圖8為所測土壤樣品KEGG直系同源基因簇(KO)豐度熱圖,橫坐標為不同生物炭用量處理,縱坐標表示代謝相關基因編碼酶的豐度。分析發(fā)現(xiàn),氨基酸代謝、碳水化合物代謝和能量代謝通路上的相關基因表達差異較明顯;可大致聚為CK和生物炭兩類,不同生物炭用量間也存在差異,其中T4與其他3個生物炭處理差異相對較大。以細菌群落多樣性指數(shù)較高的T2處理與CK相比,土壤代謝相關基因豐度升高的占比較大,包括與能源物質代謝密切相關的轉酮醇酶(K00615),與碳代謝相關的脫氫酶(K00058、K00382、K01754、K01687),與脂肪酸合成代謝相關的?;d體蛋白(K02078)、3-氧?;埘;d體蛋白]合酶(K00648、K09458)及?;D移酶(K00655),作用于組氨酸的氨基轉移酶(K00817)及支鏈氨基酸的氨基轉移酶(K00826),與色氨酸代謝相關酶(K00128),參與核酸錯配修復的酶(K03657),參與氨基酸合成相關酶(K01733)及ABC轉運蛋白相關酶(K02066、K02067)。

    3 討論

    通常認為,生物炭本身呈堿性,其高pH主要源于組分中的灰分包含大量K、Ca、Na等金屬礦物,可有效降低土壤交換性H+含量,因此對于增加酸性土壤pH具有明顯作用(Van Zwieten et al.,2010)。研究表明,施用生物炭后土壤總氮、有機質、有效磷和速效鉀含量顯著上升,且隨著生物炭用量的增加效果較明顯(張廣雨等,2020);煙稈生物炭施用于田間后可提高烤煙根際土壤pH及有機碳含量(Jones et al.,2012)。本研究同樣發(fā)現(xiàn),隨著生物炭用量逐漸增加,土壤pH、有機質、總氮、有效磷、速效鉀含量呈增加趨勢,可能與生物炭和肥料的互補或協(xié)同作用有關。生物炭巨大的內(nèi)表面積、較高的孔隙結構、豐富的官能團能吸附多種離子,延緩土壤中肥料養(yǎng)分的釋放,提高土壤保肥性能和碳氮比(李彩斌等,2018),調節(jié)土壤肥力平衡,改變土壤微生物定殖棲息地環(huán)境(張璐等,2019)。此外,土壤pH是影響微生物生長代謝的主要因子之一;當pH從酸性向中性范圍變化時,土壤微生物量會隨著pH的升高而增加(Aciego et al.,2008)。本研究中,施入生物炭后土壤細菌豐度顯著增加,可能與生物炭的多孔結構、土壤pH增加及微生態(tài)環(huán)境的改善等多重因素疊加有關(Oni et al.,2019)。生物炭中部分易分解炭可作為微生物生長的基質,施用生物炭可刺激土壤微生物的繁殖和活力,并發(fā)生起爆效應(Chen et al.,2017);適量生物炭可提高土壤細菌群落的豐度和多樣性,但超過一定用量后細菌OTUs數(shù)量則會下降,致病菌豐度增加(鄧建強等,2018),其原因可能與高用量生物炭對土壤性狀產(chǎn)生的負面影響有關(王光飛等,2017)。本研究發(fā)現(xiàn),施用生物炭處理的土壤細菌群落Alpha多樣性指數(shù)相對較高,尤其以15 t/ha處理的土壤細菌群落較豐富。因此,對于貴州煙區(qū)黃棕壤土壤細菌群落組成來說,結合煙稈生物炭自身特性、土壤類型及試驗地條件等多種因素,施用該用量的煙稈生物炭較理想。

    酸桿菌門是土壤生態(tài)系統(tǒng)中一個較豐富、難以培養(yǎng)的細菌群落,在關鍵碳、氮和硫等生物地質化學循環(huán)中扮演重要生態(tài)角色(Kalam et al.,2020)。大多數(shù)酸桿菌更喜歡在酸性條件下生長,較低的pH條件下酸桿菌豐度通常較高(Lladó et al.,2018)。本研究中,由于生物炭施用提高了土壤pH,可能是酸桿菌豐度隨著生物用量增加而降低的原因之一。研究表明,在土壤微生物區(qū)系中,擬桿菌門能分泌多種碳水化合物活性酶,降解來源于動植物殘體的多聚糖及纖維素(Larsbrink and McKee,2020)。施用生物炭可顯著提高土壤放線菌門和擬桿菌門的多樣性和生物量,改善土壤微生物群落結構(Liu et al.,2020)。在生物炭改良土壤的根際相關群落中,擬桿菌門的相對豐度較高(Kolton et al.,2011)。Liao等(2019)利用13CO2標記,發(fā)現(xiàn)在豆科間作系統(tǒng)中,施用生物炭不僅改變了根際細菌群落結構和多樣性,還提高了擬桿菌成員(芽孢桿菌、梭狀芽孢桿菌等)的相對豐度。本研究發(fā)現(xiàn),生物炭用量與土壤擬桿菌門顯著正相關,在15或20 t/ha用量時,土壤放線菌屬、芽單胞菌屬、類芽孢桿菌屬相對豐度明顯增加。放線菌是產(chǎn)生抗生素的主要細菌之一,土壤中放線菌數(shù)量及其比例對土壤微生物生態(tài)平衡調整有重要作用(張璐等,2019);芽孢桿菌屬、假單胞菌屬為植物根際促生菌,其主要功能是分泌植物生長調節(jié)劑而促進植物生長(戴華鑫等,2020)。因此,從細菌群落多樣性指數(shù)及細菌屬功能來看,生物炭施用量在15或20 t/ha時效果較好。

    農(nóng)田土壤中存在大量具有特定生態(tài)功能的活躍微生物,這些微生物與新陳代謝活動緊密相關(LeBrun and Kang,2018);短期生物炭添加后即可改變土壤細菌群落組成(Wong et al.,2019),其對土壤碳庫的固持能力可維持3年以上(Wang et al.,2021)。轉酮醇酶是磷酸戊糖途徑非氧化分支中的關鍵限速酶,對于能量轉導和產(chǎn)生用于核酸合成的核糖至關重要,在微生物碳代謝中扮演重要角色(Wilkinson and Dalby,2019);?;d體蛋白是脂肪酸合成中的關鍵蛋白質,3-氧?;埘;d體蛋白]合酶是一種長鏈脂肪酸合成酶,烯酰CoA水合酶是脂肪酸合成和分解代謝的關鍵分子(Khurana et al.,2010)。本研究中,上述與脂肪酸合成、能源物質及碳代謝有關酶編碼基因豐度高于對照,說明生物炭施入土壤4年后土壤細菌碳源利用能力仍處于較高水平。支鏈氨基酸是一類具有多種生理和代謝作用的必需氨基酸(呂娜娜等,2018)。研究表明,氨基酸能促進土壤微生物活動,加快有機物礦化和土壤養(yǎng)分稀釋,增加土壤微生物數(shù)量(張貴云等,2019);生物炭處理顯著提升了與能量生產(chǎn)和轉換、氨基酸轉運代謝、核苷酸及碳水化合物轉運代謝等相關功能的基因表達量(Gao et al.,2021)。ABC轉運基因的功能不僅局限于轉運,還包括通道調控和分子量變換等(Gong and Wang,2022),且與根系分泌物的釋放量有關(Cox et al.,2019)。本研究發(fā)現(xiàn),細菌豐度隨著生物炭施入量增加而增加,且生物炭處理后的土壤養(yǎng)分相對較高,土壤細菌類群功能代謝及相關酶基因豐度改變可能是其原因之一。

    4 結論

    貴州植煙黃棕壤施用煙稈生物炭后,土壤養(yǎng)分、烤煙農(nóng)藝性狀和經(jīng)濟性狀指標均有提升,土壤細菌群落豐富度及多樣性顯著提高,土壤細菌碳代謝能力增強,對改善煙田土壤微生態(tài)環(huán)境具有積極作用。其中以15 t/ha煙稈生物炭作用效果較好,可在生產(chǎn)上應用推廣。

    參考文獻:

    白珊,倪幸,楊瑗羽,方先芝,柳丹,葉正錢. 2021. 不同原材料生物炭對土壤重金屬Cd、Zn的鈍化作用[J]. 江蘇農(nóng)業(yè)學報,37(5):1199-1205. [Bai S,Ni X,Yang Y Y,F(xiàn)ang X Z,Liu D,Ye Z Q. 2021. Immobilization of soil cadmium and zinc by different raw material derived biochars[J]. Jiangsu Journal of Agricultural Sciences,37(5):1199-1205.] doi:10.3969/j.issn.1000-4440.2021.05.015.

    陳懿,陳偉,高維常,程建中,林葉春,潘文杰. 2017. 煙稈生物炭對烤煙根系生長的影響及其作用機理[J]. 煙草科技,50(6):26-32. [Chen Y,Chen W,Gao W C,Cheng J Z,Lin Y C,Pan W J. 2017. Effects of tobacco stalk biochar on root growth of flue-cured tobacco and its action me-chanism[J]. Tobacco Science & Technology,50(6):26-32.] doi:10.16135/j.issn1002-0861.2016.0268.

    戴華鑫,陳麗燕,陳彥春,胡利偉,王建偉,梁太波,張仕祥,陳釗,張艷玲,尹啟生. 2017. 豫中南煙區(qū)不同質地土壤理化性質、酶活性及微生物群落分析[J]. 煙草科技,50(9):7-14. [Dai H X,Chen L Y,Chen Y C,Hu L W,Wang J W,Liang T B,Zhang S X,Chen Z,Zhang Y L,Yin Q S. 2017. Effects of soil texture on physicochemical properties,enzyme activities and microbial communities of soils in central and south Henan tobacco-planting areas[J]. Tobacco Science & Technology,50(9):7-14.] doi:10.16135/j.issn1002-0861.2017.0010.

    戴華鑫,張艷玲,段衛(wèi)東,陳小龍,蘇新宏,任應斌,劉文濤,李亮,毛家偉,張翔. 2020. 豆?jié){灌根對烤煙生長及土壤細菌群落的影響[J]. 煙草科技,53(4):1-10. [Dai H X,Zhang Y L,Duan W D,Chen X L,Su X H,Ren Y B,Liu W T,Li L,Mao J W,Zhang X. 2020. Effects of soybean milk root irrigation on growth of flue-cured tobacco and soil bacterial communities[J]. Tobacco Science & Technology,53(4):1-10.] doi:10.16135/j.issn1002-0861. 2019.0463.

    鄧建強,譚軍,施河麗,樊俊,向必坤,王瑞. 2018. 生物炭對土地整治區(qū)土壤微生物調控效應[J]. 中國煙草學報,24(3):46-52. [Deng J Q,Tan J,Shi H L,F(xiàn)an J,Xiang B K,Wang R. 2018. Control effect of biochar on soil microorganism in land consolidation region[J]. Acta Taba-caria Sinica,24(3):46-52.] doi:10.16472/j.chinatobacco.2017.291.

    桂意云,李海碧,韋金菊,毛蓮英,張榮華,區(qū)惠平,祝開,趙培方,周會,劉昔輝. 2022. 生物炭對旱坡地宿根甘蔗土壤養(yǎng)分、酶活性及微生物多樣性的影響[J]. 南方農(nóng)業(yè)學報,53(3):776-784. [Gui Y Y,Li H B,Wei J J,Mao L Y,Zhang R H,Ou H P,Zhu K,Zhao P F,Zhou H,Liu X H. 2022. Effects of biochar on soil nutrients,enzyme activities and microbial diversity of ratoon sugarcane in dry slope land[J]. Journal of Southern Agriculture,53(3):776-784.] doi:10.3969/j.issn.2095-1191.2022.03.019.

    黃雁飛,陳桂芬,熊柳梅,劉斌,劉永賢,黃玉溢,唐其展. 2020. 不同秸稈生物炭對水稻生長及土壤養(yǎng)分的影響[J]. 南方農(nóng)業(yè)學報,51(9):2113-2119. [Huang Y F,Chen G F,Xiong L M,Liu B,Liu Y X,Huang Y Y,Tang Q Z. 2020. Effects of different straw biochars on rice growth and soil nutrients[J]. Journal of Southern Agriculture,51(9):2113-2119.] doi:10.3969/j.issn.2095-1191. 2020.09.008.

    況帥,段焰,劉芮,程昌新,胡志明,劉浩,何曉健,王松峰,宋文靜,叢萍. 2021. 油菜秸稈生物炭對植煙紅壤養(yǎng)分及細菌群落多樣性的影響[J]. 中國煙草科學,42(1):20-26. [Kuang S,Duan Y,Liu R,Cheng C X,Hu Z M,Liu H,He X J,Wang S F,Song W J,Cong P. 2021. Effects of straw biochar on soil nutrients and diversity of bacterial communities in tobacco-planting red soil[J]. Chinese Tobacco Science,42(1):20-26.] doi:10.13496/j.issn.1007-5119.2021.01.004.

    李彩斌,張久權,陳雪,張繼光,翟欣,凌愛芬. 2018. 生物炭施用對土壤健康的影響及其對烤煙生產(chǎn)的潛在風險[J]. 中國煙草科學,39(6):91-97. [Li C B,Zhang J Q,Chen X,Zhang J G,Zhai X,Ling A F. 2018. Effect of biochar application on soil health and its potential risks to flue-cured tobacco production[J]. Chinese Tobacco Science,39(6):91-97.] doi:10.13496/j.issn.1007-5119.2018.06.013.

    李青山,王德權,杜傳印,王毅,高凱,王剛,王慎強. 2021. 有機無機肥與生物炭配施對烤煙生長發(fā)育和煙葉質量的影響[J]. 土壤通報,52(6):1393-1401. [Li Q S,Wang D Q,Du C Y,Wang Y,Gao K,Wang G,Wang S Q. 2021. Effect of combined application of organic and inorganic fertilizers and biochar on the growth and development of flue-cured tobacco and leaf quality[J]. Chinese Journal of Soil Science,52(6):1393-1401.] doi:10.19336/j.cnki.trtb.2021012601.

    李松偉,王發(fā)展,陳彪,謝俊明,史久長. 2021. 生物炭配施海藻肥對連作植煙土壤理化特性及烤煙生長的影響[J]. 河南農(nóng)業(yè)大學學報,55(5):852-861. [Li S W,Wang F Z,Chen B,Xie J M,Shi J C. 2021. Effects of biochar combined with seaweed fertilizer on the physicochemical characteristics of continuous tobacco planting soil and flue-cured tobacco growth[J]. Journal of Henan Agricultural University,55(5):852-861.] doi:1016445/j.cnki. 1000-2340.20210804.001.

    李文豪. 2019. 煙草秸稈有機肥的研制與土壤改良應用[D]. 武漢:湖北大學. [Li W H. 2019. Development and application of organic tobacco straw fertilizer for soil improvement[D]. Wuhan:Hubei University.] doi:10.27130/d.cnki. ghubu.2019.000468.

    魯如坤. 2000. 土壤農(nóng)業(yè)化學分析方法[M]. 北京:中國農(nóng)業(yè)科技出版社. [Lu R K. 2000. Analytical methods of soil agricultural chemistry[M]. Beijing:China Agricultural Science and Technology Press.]

    呂娜娜,沈宗專,王東升,劉紅軍,薛超,李榮,沈啟榮. 2018. 施用氨基酸有機肥對黃瓜產(chǎn)量及土壤生物學性狀的影響[J]. 南京農(nóng)業(yè)大學學報,41(3):456-464. [Lü N N,Shen Z Z,Wang D S,Liu H J,Xue C,Li R,Shen Q R. 2018. Effects of amino acid organic fertilizer on cucumber yield and soil biological characters[J]. Journal of Nanjing Agricultural University,41(3):456-464.] doi:10.7685/jnau.201709023.

    王光飛,馬艷,郭德杰,曹云,羅德旭,趙建鋒,孫玉東. 2017. 不同用量秸稈生物炭對辣椒疫病防控效果及土壤性狀的影響[J]. 土壤學報,54(1):204-215.[Wang G F,Ma Y,Guo D J,Cao Y,Luo D X,Zhao J F,Sun Y D. 2017. Application-rate-dependent effects of straw biochar on control of phytophthora blight of chilli pepper and soil properties[J]. Acta Pedologica Sinica,54(1):204-215.] doi:10.11766/trxb20160414002

    王志丹,代曉華,劉吉利,于祥,吳娜. 2022. 生物炭對寧夏干旱地區(qū)玉米生長及產(chǎn)量的影響[J]. 甘肅農(nóng)業(yè)大學學報,57(1):74-82. [Wang Z D,Dai X H,Liu J L,Yu X,Wu N. 2022. Effect of biochar on growth and yield of maize in arid area of Ningxia[J]. Journal of Gansu Agricultural University,57(1):74-82.] doi:10.13432/j.cnki.jgsau.2022. 01.009.

    姚麗茹,李偉,朱員正,普布倉決,韓娟. 2022. 施用生物炭對麥田土壤細菌群落多樣性和冬小麥生長的影響[J/OL]. 環(huán)境科學,https://kns.cnki.net/kcms/detail/11.1895.X.20220919. 1141.032.html. [Yao L R,Li W,Zhu Y Z,Pubu C J,Han J. 2022. Effects of biochar application on soil bacterial community diversity and winter wheat growth in wheat fields[J/OL]. Environmental Science,https://kns.cnki.net/kcms/detail/11.1895.X.20220919. 1141.032.html.] doi:10.13227/j.hjkx.202207125.

    葉協(xié)鋒,李志鵬,于曉娜,牛瑜德,李先鋒,管賽賽,宗勝杰,孟琦. 2015. 生物炭用量對植煙土壤碳庫及烤后煙葉質量的影響[J]. 中國煙草學報,21(5):33-41. [Ye X F,Li Z P,Yu X N,Niu Y D,Li X F,Guan S S,Zong S J,Meng Q. 2015. Effect of biochar application rate on quality of flue-cured tobacco leaves and carbon pool in tobacco growing soil[J]. Acta Tabacaria Sinica,21(5):33-41.] doi:10.16472/j.chinatobacco.2014.546.

    曾鳳鈴,鄒玉霞,張衛(wèi)華,運劍葦,張偉,邵清清. 2022. 生物炭用量與粒徑對紫色土水分入滲及再分布特征的影響[J]. 西南大學學報(自然科學版),44(12):136-144. [Zeng F L,Zou Y X,Zhang W H,Yun J W,Zhang W,Shao Q Q. 2022. Effects of rateand particle size of biochar on water infiltration and redistribution characteristics of purple soil[J]. Journal of Southwest University(Natural Science Edition),44(12):136-144.] doi:10.13718/j.cnki.xdzk. 2022.12.014.

    張廣雨,胡志明,褚德朋,戴永平,尹紹靜,嚴明雄,段焰,孫康,Naila Ilyas,尤祥偉,李義強. 2020. 生物炭對根際土壤微生態(tài)的調控及對煙草青枯病的防控作用[J]. 中國煙草學報,26(6):81-88. [Zhang G Y,Hu Z M,Chu D P,Dai Y P,Yin S J,Yan M X,Duan Y,Sun K,Naila I,You X W,Li Y Q. 2020. Regulation of biochar on rhizosphere soil microecology and its control effect on tobacco bacterial wilt[J]. Acta Tabacaria Sinica,26(6):81-88.] doi:10.16472/j.chinatobacco.2020.135.

    張貴云,呂貝貝,張麗萍,劉珍,范巧蘭,魏明峰,姚眾,袁嘉瑋,柴躍進. 2019. 黃土高原旱地麥田26年免耕覆蓋對土壤肥力及原核微生物群落多樣性的影響[J]. 中國生態(tài)農(nóng)業(yè)學報,27(3):358-368. [Zhang G Y,Lü B B,Zhang L P,Liu Z,F(xiàn)an Q L,Wei M F,Yao Z,Yuan J W,Chai Y J. 2019. Effect of long-term no-tillage with stubble on soil fertility and diversity of prokaryotic microbiome in dryland wheat soils on the Loess Plateau,China[J]. Chinese Journal of Eco-Agriculture,27(3):358-368.] doi:10.13930/j.cnki.cjea.180604.

    張珂珂,宋曉,郭斗斗,黃紹敏,岳克,郭騰飛,張水清,岳艷軍. 2022. 生物炭對潮土土壤肥力特征和氮肥利用效率的影響[J]. 河南農(nóng)業(yè)科學,51(12):73-80. [Zhang K K,Song X,Guo D D,Huang S M,Yue K,Guo T F,Zhang S Q,Yue Y J. 2022. Effects of biochar on fertility and nitrogen utilization efficiency in fluvo-aquic soil[J]. Journal of Henan Agricultural Sciences,51(12):73-80.] doi:10.15933/j.cnki.1004-3268.2022.12.009.

    張璐,閻海濤,任天寶,李帥,楊永鋒,彭桂新,于建春,劉國順. 2019. 有機物料對植煙土壤養(yǎng)分、酶活性和微生物群落功能多樣性的影響[J]. 中國煙草學報,25(2):55-62. [Zhang L,Yan H T,Ren T B,Li S,Yang Y F,Peng G X,Yu J C,Liu G S. 2019. Effects of organic matter on nu-trient,enzyme activity and functional diversity of microbial community in tobacco planting soil[J]. Acta Tabacaria Sinica,25(2):55-62.] doi:10.16472/j.chinatobacco. 2018.195.

    Aciego Pietri J C,Brookes P C. 2008. Relationships between soil pH and microbial properties in a UK arable soil[J]. Soil Biology and Biochemistry,40(7):1856-1861.doi:10.1016/j.soilbio.2008.03.020.

    Blanco-Canqui H. 2017. Biochar and soil physical properties[J]. Soil Science Society of America Journal,81(4):? ? ? 687-711. doi:10.2136/sssaj2017.01.0017.

    Chen J H,Li S H,Liang C F,Xu Q F,Li Y C,Qin H,F(xiàn)uhrmann J J. 2017. Response of microbial community structure and function to short-term biochar amendment in an intensively managed bamboo(Phyllostachys praecox)plantation soil:Effect of particlesize and addition rate[J]. Science of the Total Environment,574:24-33. doi:10.1016/j.scitotenv.2016.08.190.

    Cox D E,Dyer S,Weir R,Cheseto X,Sturrock M,Coyne D,Torto B,Maule A G,Dalzell J J. 2019. ABC transporter genes ABC-C6 and ABC-G33 alter plant-microbe-parasite interactions in the rhizosphere[J]. Scientific Reports,9(1):19899. doi:10.1038/s41598-019-56493-w.

    Gao M L,Zhang Z,Dong Y M,Song Z G,Dai H X. 2020. Responses of bacterial communities in wheat rhizospheres in different soils to di-n-butyl and di(2-ethylhexyl)phtha-late contamination[J]. Geoderma,362:114126.doi:10.1016/j.geoderma.2019.114126.

    Gao M Y,Yang J F,Liu C M,Gu B W,Han M,Li J W,Li N,Liu N,An N,Dai J,Liu X H,Han X R. 2021. Effects of long-term biochar and biochar-based fertilizer application on brown earth soil bacterial communities[J]. Agriculture,Ecosystems & Environment,309:107285.doi:10.1016/j.agee.2020.107285.

    Gong X P,Wang S H. 2022. New insights into evolution of the ABC transporter family in Mesostigma viride,a unicellular charophyte algae[J]. Current Issues in Molecular Biology,44(4):1646-1660. doi:10.3390/cimb44040112.

    Jones D L,Rousk J,Edwards-Jones G,DeLuca T H,Murphy D V. 2012. Biochar-mediated changes in soil quality and plant growth in a three year field trial[J]. Soil Biology and Biochemistry,45:113-124. doi:10.1016/j.soilbio.2011. 10.012.

    Kalam S,Basu A,Ahmad I,Sayyed R Z,El-Enshasy H A,Dailin D J,Suriani N L. 2020. Recent understanding of soil acidobacteria and their ecological signi?cance:A critical review[J]. Frontiers in Microbiology,11:580024. doi:10. 3389/fmicb.2020.580024.

    Khurana P,Gokhale R S,Mohanty D. 2010. Genome scale prediction of substrate specificity for acyl adenylate superfamily of enzymes based on active site residue profiles[J]. BMC Bioinformatics,11(1):57. doi:10.1186/1471-2105-11-57.

    Kolton M,Harel Y M,Pasternak Z,Graber E R,Elad Y,Cytryn E. 2011. Impact of biochar application to soil on the root-associated bacterial community structure of fully developed greenhouse pepper plants[J]. Applied and Environmental Microbiology,77(14):4924-4930. doi:10.1128/AEM.00148-11.

    Langille M G I,Zaneveld J,Caporaso J G,McDonald D,Knights D,Reyes J A,Clemente J C,Burkepile D E,Thurber R L V,Knight R,Beiko R G,Huttenhower C. 2013. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences[J]. Nature Biotechnology,31(9):814-821. doi:10.1038/nbt.2676.

    Larsbrink J,McKee L S. 2020. Bacteroidetes bacteria in the soil:Glycan acquisition,enzyme secretion,and gliding motility[J]. Advances in Applied Microbiology,110:63-98. doi:10.1016/bs.aambs.2019.11.001.

    LeBrun E S,Kang S. 2018. A comparison of computationally predicted functional metagenomes and microarray analysis for microbial P cycle genes in a unique basalt-soil forest[J]. F1000Research,7:179. doi:10.12688/f1000research. 13841.1.

    Liao H K,Li Y Y,Yao H Y. 2019. Biochar amendment stimulates utilization of plant-derived carbon by soil bacteria in an intercropping system[J]. Frontiers in Microbiology,10:1361. doi:10.3389/fmicb.2019.01361.

    Liu W,Li Y L,F(xiàn)eng Y,Qiao J C,Zhao H W,Xie J X,F(xiàn)ang Y Y,Shen S G,Liang S X. 2020. The effectiveness of nanobiochar for reducing phytotoxicity and improving soil remediation in cadmium-contaminated soil[J]. Scientific Reports,10(1):858. doi:10.1038/s41598-020-57954-3.

    Lladó S,López-Mondéjar R,Baldrian P. 2018. Drivers of microbial community structure in forest soils[J]. Applied Microbiology and Biotechnology,102(10):4331-4338. doi:10.1007/s00253-018-8950-4.

    Luo X S,F(xiàn)u X Q,Yang Y,Cai P,Peng S B,Chen W L,Huang Q Y. 2016. Microbial communities play important roles in modulating paddy soil fertility[J]. Scientific Reports,6(1):20326. doi:10.1038/srep20326.

    Oni B A,Oziegbe O,Olawole O O. 2019. Significance of biochar application to the environment and economy[J]. Annals of Agricultural Sciences,64(2):222-236. doi:10.1016/j.aoas.2019.12.006.

    Van Zwieten L,Kimber S,Morris S,Chan K Y,Downie A,Rust J,Joseph S,Cowie A. 2010. Effects of biochar from slow pyrolysis of papermill waste on agronomic performance and soil fertility[J]. Plant Soil,327:235-246. doi:10.1007/s11104-009-0050-x.

    Wang C,Chen D,Shen J L,Yuan Q,F(xiàn)an F L,Wei W X,Li Y,Wu J S. 2021.Biochar alters soil microbial communities and potential functions 3-4 years after amendment in a double rice cropping system[J]. Agriculture,Ecosystems & Environment,311:107291. doi:10.1016/j.agee.2020.107291.

    Wang H,Liu S R,Schindlbacher A,Wang J X,Yang Y J,Song Z C,You Y M,Shi Z M,Li Z Y,Chen L,Ming A G,Lu L H,Cai D X. 2019. Experimental warming reduced topsoil carbon content and increased soil bacterial diversity in a subtropical planted forest[J]. Soil Biology and Biochemistry,133:155-164. doi:10.1016/j.soilbio.2019.03.004.

    Wang L W,OConnor D,Rinklebe J,Ok Y S,Tsang D C W,Shen Z T,Hou D Y. 2020. Biochar aging:Mechanisms,physicochemical changes,assessment,and implications for field applications[J]. Environmental Science & Technology,54(23):14797-14814. doi:10.1021/acs.est.0c04033.

    Wilkinson H C,Dalby P A. 2019. Novel insights into transketolase activation by cofactor binding identifies two native species subpopulations[J]. Scientific Reports,9(1):16116. doi:10.1038/s41598-019-52647-y.

    Wong J T F,Chen X W,Deng W J,Chai Y M,Ng C W W,Wong M H. 2019. Effects of biochar on bacterial communities in a newly established landfill cover topsoil[J]. Journal of Environmental Management,236:667-673. doi:10.1016/j.jenvman.2019.02.010.

    Zhang H J,Wang S J,Zhang J X,Tian C J,Luo S S. 2021. Biochar application enhances microbial interactions in mega-aggregates of farmland black soil[J]. Soil and Tilla-ge Research,213:105145. doi:10.1016/j.still.2021.105145.

    (責任編輯 王 暉)

    猜你喜歡
    生物炭土壤養(yǎng)分
    我國秸稈生物炭資源化的難點與對策研究
    不同施肥模式對油茶植株營養(yǎng)生長和土壤養(yǎng)分的影響
    廣靈縣平川區(qū)土壤養(yǎng)分變化及施肥建議
    生物炭的制備與表征比較研究
    人間(2016年27期)2016-11-11 17:45:25
    生物炭的應用研究進展
    稻蟹共作模式下稻蟹產(chǎn)出與土壤理化性質的研究
    促進生物炭工業(yè)生產(chǎn)與土地利用 推動21世紀農(nóng)業(yè)“黑色革命”
    孝義市不同種植方式耕作土壤養(yǎng)分狀況研究
    茶園土壤養(yǎng)分狀況與分布
    基于生物炭的生態(tài)浮床設計
    亚洲男人的天堂狠狠| 法律面前人人平等表现在哪些方面| 大型黄色视频在线免费观看| 亚洲av成人一区二区三| 国产激情偷乱视频一区二区| 怎么达到女性高潮| 色婷婷久久久亚洲欧美| 熟妇人妻久久中文字幕3abv| 老熟妇乱子伦视频在线观看| 99久久久亚洲精品蜜臀av| 少妇 在线观看| 丰满人妻熟妇乱又伦精品不卡| 久99久视频精品免费| 成人三级做爰电影| 国产精品永久免费网站| 看免费av毛片| 国产亚洲精品一区二区www| 精品国内亚洲2022精品成人| 欧美一级a爱片免费观看看 | 免费高清视频大片| 91麻豆精品激情在线观看国产| 两个人免费观看高清视频| 亚洲片人在线观看| 久久精品影院6| 国产精品久久视频播放| 色av中文字幕| 国产免费男女视频| 午夜免费成人在线视频| 19禁男女啪啪无遮挡网站| 国产av又大| 日本一本二区三区精品| 久久国产精品影院| 一区福利在线观看| 亚洲精品中文字幕一二三四区| 欧美黄色淫秽网站| 午夜日韩欧美国产| 久久精品人妻少妇| 日韩三级视频一区二区三区| 亚洲色图 男人天堂 中文字幕| 欧美久久黑人一区二区| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产| 18美女黄网站色大片免费观看| www.精华液| 美女国产高潮福利片在线看| 欧美在线黄色| 日日干狠狠操夜夜爽| 国产成人欧美在线观看| 91九色精品人成在线观看| 国产精品影院久久| 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 啦啦啦观看免费观看视频高清| 亚洲精品在线美女| 黄片小视频在线播放| www.www免费av| 级片在线观看| 侵犯人妻中文字幕一二三四区| 两个人免费观看高清视频| 亚洲精华国产精华精| 成人18禁在线播放| 熟女电影av网| 岛国视频午夜一区免费看| 在线观看免费日韩欧美大片| 制服人妻中文乱码| 熟女电影av网| 亚洲精品色激情综合| 这个男人来自地球电影免费观看| 国产又黄又爽又无遮挡在线| 欧美成人午夜精品| 色综合亚洲欧美另类图片| 免费看十八禁软件| 制服人妻中文乱码| 国产精品99久久99久久久不卡| 久久婷婷成人综合色麻豆| 在线十欧美十亚洲十日本专区| 亚洲熟妇熟女久久| 亚洲成av片中文字幕在线观看| 亚洲精品色激情综合| 美女高潮到喷水免费观看| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 欧美中文日本在线观看视频| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 国产av一区在线观看免费| 精品少妇一区二区三区视频日本电影| 国产蜜桃级精品一区二区三区| 免费看十八禁软件| 两个人看的免费小视频| 亚洲中文字幕一区二区三区有码在线看 | 淫秽高清视频在线观看| 免费在线观看亚洲国产| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 亚洲最大成人中文| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡免费网站照片 | 一a级毛片在线观看| 老鸭窝网址在线观看| 级片在线观看| 亚洲精品国产区一区二| 国产精品九九99| 很黄的视频免费| 中文字幕精品免费在线观看视频| 中文字幕人成人乱码亚洲影| 国产精品 国内视频| 91大片在线观看| 精品一区二区三区视频在线观看免费| 1024香蕉在线观看| 大型av网站在线播放| aaaaa片日本免费| 日日爽夜夜爽网站| 天天一区二区日本电影三级| 韩国精品一区二区三区| 国产一卡二卡三卡精品| 黑丝袜美女国产一区| 大型av网站在线播放| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美精品综合一区二区三区| 深夜精品福利| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 91老司机精品| 午夜福利高清视频| av在线天堂中文字幕| 色av中文字幕| 亚洲中文字幕日韩| 国产一区二区在线av高清观看| 久久久精品国产亚洲av高清涩受| 成人永久免费在线观看视频| 久久狼人影院| 亚洲国产精品成人综合色| 成人国语在线视频| 国产1区2区3区精品| 黄色 视频免费看| 久久天堂一区二区三区四区| 欧美日韩乱码在线| 亚洲男人的天堂狠狠| 国产亚洲精品久久久久久毛片| 亚洲熟妇熟女久久| 国产亚洲欧美98| 色老头精品视频在线观看| 欧美午夜高清在线| 久久香蕉激情| 女人爽到高潮嗷嗷叫在线视频| 日韩 欧美 亚洲 中文字幕| 桃红色精品国产亚洲av| 日日夜夜操网爽| 亚洲三区欧美一区| 两个人免费观看高清视频| 久久久精品欧美日韩精品| 久久香蕉激情| 国产成人av教育| 日本一本二区三区精品| 在线十欧美十亚洲十日本专区| 久久久久久免费高清国产稀缺| 99国产极品粉嫩在线观看| 啦啦啦 在线观看视频| cao死你这个sao货| 成人国语在线视频| 搡老熟女国产l中国老女人| 91字幕亚洲| 男女之事视频高清在线观看| 亚洲天堂国产精品一区在线| 12—13女人毛片做爰片一| 国产精品自产拍在线观看55亚洲| 老司机福利观看| 亚洲三区欧美一区| 久久久精品欧美日韩精品| 国产成人啪精品午夜网站| 不卡一级毛片| 极品教师在线免费播放| 午夜成年电影在线免费观看| 色哟哟哟哟哟哟| 亚洲性夜色夜夜综合| 国产精品 欧美亚洲| 每晚都被弄得嗷嗷叫到高潮| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 亚洲av熟女| 久久久久久久久免费视频了| 国产精品久久久久久亚洲av鲁大| 欧美zozozo另类| 亚洲人成网站高清观看| av中文乱码字幕在线| 国产精品 国内视频| 国产精品野战在线观看| www日本黄色视频网| 欧美乱妇无乱码| 又黄又粗又硬又大视频| 香蕉av资源在线| 黄色片一级片一级黄色片| 亚洲国产精品合色在线| 亚洲av成人av| 老司机午夜十八禁免费视频| 亚洲 欧美一区二区三区| av天堂在线播放| 90打野战视频偷拍视频| 在线观看免费午夜福利视频| 亚洲av中文字字幕乱码综合 | 国产成人精品无人区| 午夜激情福利司机影院| svipshipincom国产片| 国产免费男女视频| 国产av在哪里看| 免费高清在线观看日韩| 亚洲欧美激情综合另类| 一区二区三区国产精品乱码| 丝袜在线中文字幕| 成人亚洲精品一区在线观看| 一级作爱视频免费观看| 欧美激情久久久久久爽电影| 久久精品影院6| 色av中文字幕| 精品电影一区二区在线| 夜夜爽天天搞| 看黄色毛片网站| 国产免费av片在线观看野外av| 久久这里只有精品19| av免费在线观看网站| 久久久久久久午夜电影| 欧洲精品卡2卡3卡4卡5卡区| 大型黄色视频在线免费观看| 精品久久久久久久末码| 热re99久久国产66热| 一进一出好大好爽视频| 老司机午夜十八禁免费视频| 日韩中文字幕欧美一区二区| 久久中文字幕一级| 黄色a级毛片大全视频| 亚洲国产精品999在线| 精品欧美国产一区二区三| 不卡av一区二区三区| 变态另类丝袜制服| 日本三级黄在线观看| 亚洲av成人不卡在线观看播放网| 色精品久久人妻99蜜桃| 欧美日本视频| 欧美黄色片欧美黄色片| 色综合欧美亚洲国产小说| 免费在线观看影片大全网站| a级毛片在线看网站| 9191精品国产免费久久| 亚洲成人国产一区在线观看| 久久精品国产99精品国产亚洲性色| 久久香蕉精品热| 丝袜美腿诱惑在线| 看免费av毛片| 国产免费av片在线观看野外av| 桃色一区二区三区在线观看| 精品国产亚洲在线| 狠狠狠狠99中文字幕| 亚洲精品av麻豆狂野| 国产精品久久久久久精品电影 | 在线永久观看黄色视频| 国产成人欧美在线观看| 国内少妇人妻偷人精品xxx网站 | 国产高清有码在线观看视频 | 精品久久久久久久久久免费视频| 老熟妇乱子伦视频在线观看| 一本大道久久a久久精品| 亚洲成人国产一区在线观看| 一本一本综合久久| 一级毛片精品| www.自偷自拍.com| 国产精品久久久久久精品电影 | 亚洲国产欧美网| 色综合亚洲欧美另类图片| 久久午夜亚洲精品久久| 热re99久久国产66热| 午夜福利欧美成人| 日韩大码丰满熟妇| 色综合婷婷激情| 香蕉av资源在线| 一级a爱片免费观看的视频| 香蕉国产在线看| 亚洲精品av麻豆狂野| 午夜a级毛片| 中文字幕人妻熟女乱码| 真人做人爱边吃奶动态| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站 | 午夜激情av网站| 亚洲精品美女久久久久99蜜臀| 可以在线观看的亚洲视频| 国产亚洲精品综合一区在线观看 | 久久精品国产亚洲av香蕉五月| 午夜福利视频1000在线观看| 丁香六月欧美| 搡老岳熟女国产| 又大又爽又粗| 国产激情欧美一区二区| 午夜影院日韩av| 国产精品 国内视频| 啦啦啦韩国在线观看视频| 国产男靠女视频免费网站| 一个人观看的视频www高清免费观看 | 岛国在线观看网站| 日韩欧美三级三区| 老司机午夜十八禁免费视频| 午夜激情av网站| 哪里可以看免费的av片| 国产成人欧美在线观看| 中文亚洲av片在线观看爽| 黄色视频不卡| 精品日产1卡2卡| 亚洲成人免费电影在线观看| 国产一区二区在线av高清观看| 99久久精品国产亚洲精品| 一本综合久久免费| 久久精品aⅴ一区二区三区四区| 国产精华一区二区三区| 精品久久久久久久人妻蜜臀av| 美女大奶头视频| 久久久久九九精品影院| 久久国产亚洲av麻豆专区| 欧美在线黄色| 国产精品久久久久久人妻精品电影| 韩国av一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 午夜久久久久精精品| 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| 亚洲成国产人片在线观看| 淫秽高清视频在线观看| 老司机福利观看| 不卡一级毛片| av有码第一页| 在线av久久热| 手机成人av网站| 亚洲熟妇熟女久久| 一级片免费观看大全| 午夜激情福利司机影院| 女性被躁到高潮视频| 国产亚洲精品第一综合不卡| 免费在线观看影片大全网站| 丝袜美腿诱惑在线| 久久伊人香网站| 国产一级毛片七仙女欲春2 | 午夜久久久在线观看| 又紧又爽又黄一区二区| 黄色视频,在线免费观看| 中文在线观看免费www的网站 | 亚洲色图av天堂| 国产精品一区二区三区四区久久 | 久热爱精品视频在线9| 51午夜福利影视在线观看| 人妻丰满熟妇av一区二区三区| 在线视频色国产色| 欧美日韩一级在线毛片| 免费一级毛片在线播放高清视频| 婷婷六月久久综合丁香| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 桃红色精品国产亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 天天躁夜夜躁狠狠躁躁| 18美女黄网站色大片免费观看| 免费在线观看日本一区| 精品国产美女av久久久久小说| 久久久国产成人精品二区| 成人av一区二区三区在线看| 国产成人精品久久二区二区91| 欧美一级a爱片免费观看看 | 黄色女人牲交| 中文字幕精品亚洲无线码一区 | 国产成年人精品一区二区| 婷婷丁香在线五月| 最好的美女福利视频网| 午夜精品在线福利| 亚洲专区国产一区二区| 欧美激情高清一区二区三区| 怎么达到女性高潮| 日本黄色视频三级网站网址| 国产精品精品国产色婷婷| 亚洲 欧美一区二区三区| 久久久久久久久久黄片| 欧美日韩一级在线毛片| 亚洲熟女毛片儿| 亚洲五月天丁香| 美国免费a级毛片| 嫩草影院精品99| 国产亚洲欧美精品永久| 国产精品综合久久久久久久免费| 国产av在哪里看| 亚洲成av人片免费观看| 欧美精品啪啪一区二区三区| 久久香蕉国产精品| 国产精品一区二区免费欧美| 国产一区在线观看成人免费| 午夜老司机福利片| 悠悠久久av| 青草久久国产| 亚洲av成人一区二区三| 久久精品aⅴ一区二区三区四区| 成年女人毛片免费观看观看9| 每晚都被弄得嗷嗷叫到高潮| 男人的好看免费观看在线视频 | 成人免费观看视频高清| 免费在线观看亚洲国产| 黄色 视频免费看| 欧美乱妇无乱码| 熟女电影av网| 男女视频在线观看网站免费 | 国语自产精品视频在线第100页| 深夜精品福利| 日韩欧美在线二视频| 亚洲国产看品久久| 夜夜夜夜夜久久久久| 亚洲全国av大片| 欧美日韩福利视频一区二区| 丝袜在线中文字幕| 一本久久中文字幕| 免费无遮挡裸体视频| 久久精品国产亚洲av高清一级| 国产激情偷乱视频一区二区| 亚洲国产欧洲综合997久久, | 久热这里只有精品99| 性欧美人与动物交配| 亚洲精品在线观看二区| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 国产私拍福利视频在线观看| 三级毛片av免费| 国产乱人伦免费视频| 亚洲无线在线观看| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av在线| 久久午夜亚洲精品久久| 亚洲专区国产一区二区| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 国产一卡二卡三卡精品| 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 日韩高清综合在线| 成人国产综合亚洲| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 波多野结衣av一区二区av| 日本免费a在线| 日本三级黄在线观看| 亚洲美女黄片视频| 久久香蕉激情| 黄色 视频免费看| 久9热在线精品视频| 午夜两性在线视频| 99国产综合亚洲精品| 99精品在免费线老司机午夜| 一区二区日韩欧美中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇乱子伦视频在线观看| 久久久久久久精品吃奶| 国产又爽黄色视频| 久久精品91无色码中文字幕| 日韩欧美 国产精品| 午夜激情福利司机影院| 男女下面进入的视频免费午夜 | 亚洲国产欧美日韩在线播放| 日本撒尿小便嘘嘘汇集6| 午夜福利欧美成人| 50天的宝宝边吃奶边哭怎么回事| or卡值多少钱| 美女国产高潮福利片在线看| 给我免费播放毛片高清在线观看| 午夜免费成人在线视频| 久久伊人香网站| 长腿黑丝高跟| 老司机深夜福利视频在线观看| 18美女黄网站色大片免费观看| 12—13女人毛片做爰片一| 亚洲欧美激情综合另类| 日韩欧美一区视频在线观看| 人妻丰满熟妇av一区二区三区| 黄色视频,在线免费观看| 国产乱人伦免费视频| 后天国语完整版免费观看| 国产三级在线视频| 正在播放国产对白刺激| 亚洲中文av在线| 欧美国产精品va在线观看不卡| 一本久久中文字幕| 国产av在哪里看| 色尼玛亚洲综合影院| 国产黄a三级三级三级人| 国产精品亚洲美女久久久| 视频在线观看一区二区三区| 97碰自拍视频| 国产精品 欧美亚洲| 亚洲熟妇熟女久久| 在线十欧美十亚洲十日本专区| 亚洲五月婷婷丁香| 久久久久亚洲av毛片大全| 国产成人啪精品午夜网站| 99在线人妻在线中文字幕| 免费在线观看视频国产中文字幕亚洲| 日韩大码丰满熟妇| 日韩精品青青久久久久久| 熟女少妇亚洲综合色aaa.| 免费高清视频大片| 欧美国产精品va在线观看不卡| 国产亚洲av嫩草精品影院| 国产精品久久久久久亚洲av鲁大| 欧美国产日韩亚洲一区| 欧美日韩中文字幕国产精品一区二区三区| www日本在线高清视频| 青草久久国产| 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| 日本五十路高清| 久久久久久亚洲精品国产蜜桃av| 久久午夜综合久久蜜桃| 亚洲国产欧美日韩在线播放| 国产精华一区二区三区| 亚洲国产高清在线一区二区三 | 高清在线国产一区| 宅男免费午夜| 国产在线精品亚洲第一网站| 亚洲av美国av| 日日摸夜夜添夜夜添小说| 久久香蕉国产精品| 久久精品国产综合久久久| 香蕉丝袜av| 亚洲国产精品合色在线| 久久久久国内视频| 欧美日韩黄片免| 女警被强在线播放| 女人高潮潮喷娇喘18禁视频| 成年免费大片在线观看| 精华霜和精华液先用哪个| 国产精品永久免费网站| 久久精品影院6| 国产激情偷乱视频一区二区| 国产一区二区在线av高清观看| 在线国产一区二区在线| 午夜福利免费观看在线| 真人做人爱边吃奶动态| 成人亚洲精品一区在线观看| 9191精品国产免费久久| 90打野战视频偷拍视频| 波多野结衣巨乳人妻| 99国产综合亚洲精品| 又黄又粗又硬又大视频| 欧美成人性av电影在线观看| 美女国产高潮福利片在线看| 久久亚洲真实| 国产视频一区二区在线看| 俺也久久电影网| 12—13女人毛片做爰片一| 午夜激情福利司机影院| 亚洲成人精品中文字幕电影| 成人三级做爰电影| 男女做爰动态图高潮gif福利片| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 级片在线观看| 俄罗斯特黄特色一大片| 不卡一级毛片| 99热这里只有精品一区 | 欧美+亚洲+日韩+国产| 亚洲国产精品久久男人天堂| 日韩视频一区二区在线观看| 在线观看www视频免费| 可以免费在线观看a视频的电影网站| 国产欧美日韩精品亚洲av| 免费无遮挡裸体视频| 精品无人区乱码1区二区| 亚洲成人国产一区在线观看| 欧美 亚洲 国产 日韩一| 国产真实乱freesex| 国产av一区二区精品久久| netflix在线观看网站| 男人的好看免费观看在线视频 | 老汉色av国产亚洲站长工具| 欧美激情久久久久久爽电影| 亚洲人成电影免费在线| 国产精品久久久久久精品电影 | 丝袜人妻中文字幕| 波多野结衣av一区二区av| 欧美最黄视频在线播放免费| 亚洲五月天丁香| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 久久天躁狠狠躁夜夜2o2o| 精品人妻1区二区| 夜夜躁狠狠躁天天躁| 母亲3免费完整高清在线观看| 精品电影一区二区在线| 欧美国产日韩亚洲一区| 日韩欧美免费精品| 久久久久久九九精品二区国产 | 日本黄色视频三级网站网址| 老司机福利观看| 亚洲人成77777在线视频| 久久久久久免费高清国产稀缺| 亚洲精品国产一区二区精华液| 欧美激情极品国产一区二区三区| 黑人巨大精品欧美一区二区mp4| 精品久久蜜臀av无| 国产一区二区在线av高清观看| 欧美成人午夜精品| 国产v大片淫在线免费观看| 欧美性猛交黑人性爽| 在线观看www视频免费| 午夜久久久在线观看| 亚洲狠狠婷婷综合久久图片| 精品国产亚洲在线| 亚洲精品国产一区二区精华液| 极品教师在线免费播放|