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

    水分管理對水稻籽粒硒積累及根際土壤細菌群落多樣性的影響*

    2021-11-15 05:23:00王瑞昕張思宇周鑫斌
    土壤學報 2021年6期
    關(guān)鍵詞:菌門根際籽粒

    王瑞昕,楊 靜,方 正,張思宇,周鑫斌

    水分管理對水稻籽粒硒積累及根際土壤細菌群落多樣性的影響*

    王瑞昕,楊 靜,方 正,張思宇,周鑫斌??

    (西南大學資源環(huán)境學院,重慶 400716)

    采用盆栽試驗方法,研究不同水分管理方式對水稻根際土壤硒組分、籽粒硒積累以及根際土壤細菌群落多樣性的影響。結(jié)果表明:在水稻的各生育期,好氧和干濕交替較淹水灌溉一定程度上提高了土壤pH和氧化還原電位(Eh),土壤水溶態(tài)和可交換態(tài)硒含量增加,從而提高了土壤硒的有效性。水稻成熟后,不同部位的含硒量由高到低依次為根(1.38~2.22 mg·kg–1)、葉(0.55~0.85 mg·kg–1)、莖(0.53~0.61 mg·kg–1)和籽粒(0.15~0.53 mg·kg–1)。水稻籽粒含硒量以干濕交替灌溉最高,淹水灌溉最低,二者含硒量差異達顯著水平。干濕交替灌溉的水稻產(chǎn)量顯著高于常規(guī)淹水灌溉,且較淹水灌溉提高了7.83%,較好氧灌溉提高13.51%。水稻根際土壤優(yōu)勢菌為變形菌門、綠彎菌門、擬桿菌門、酸桿菌門、Patescibacteria和芽單胞菌門,變形菌門是不同水分管理方式下水稻根際土壤中豐度最高的細菌,水分管理措施顯著影響其豐富度,干濕交替和好氧灌溉中變形菌門的豐富度明顯高于淹水灌溉。從綱水平看,Gammaproteobacteria的豐度與土壤有效硒含量呈正相關(guān),故Gammaproteobacteria豐度的增加可能是土壤硒生物有效性增加的另一個重要原因。綜上,干濕交替灌溉不但能提高水稻產(chǎn)量和稻米硒含量,且較正常淹水管理節(jié)約用水,在水稻生產(chǎn)中,是一種優(yōu)先推薦的水分管理方法。

    水稻;硒;干濕交替;根際土壤;細菌多樣性

    硒是人體必需的微量營養(yǎng)元素,人體硒攝入不足會導致一系列健康問題,諸如免疫力低下、生育能力降低和增加癌癥患病風險[1]。通過主食稻米來增加人體硒攝入是一條非常理想的途徑,然而75%的水稻樣品不能夠滿足人類日常硒需求[2]。葉面噴硒可增加水稻籽粒硒含量,但長期噴硒不僅加大生產(chǎn)成本,還可能造成局部環(huán)境污染,過量易對人體健康造成毒害[3]。如能通過合理的農(nóng)藝措施來調(diào)控土壤硒的生物有效性,增加水稻可食部分的硒含量,從而增加人類硒攝入,是硒生物強化最理想的策略,是從源頭上調(diào)控硒營養(yǎng)的根本出路。

    土壤硒生物地球化學循環(huán)由化學轉(zhuǎn)化過程和生物轉(zhuǎn)化過程共同參與。土壤氧化還原電位(Eh)和酸堿度(pH)對土壤硒形態(tài)的化學轉(zhuǎn)化起著主要作用[4],而水分管理措施很大程度上影響著稻田土壤Eh和pH,最終控制著稻田硒的生物有效性[5]。在好氧的中性至堿性土壤環(huán)境中,硒主要以硒酸鹽形態(tài)存在;在厭氧的土壤環(huán)境中,主要以亞硒酸鹽形態(tài)存在。在酸性條件下,元素硒和硒化物則是主要的存在形態(tài)[6]。不同形態(tài)硒的生物有效性不同,硒酸鹽和亞硒酸鹽是水溶性的,其中硒酸鹽不能被土壤吸附,而亞硒酸鹽可被土壤固相中鐵錳氧化物或水化物牢牢吸附,導致溶解性降低[7],所以硒酸鹽較亞硒酸鹽具有更高的植物可利用性。因此,科研人員試圖發(fā)現(xiàn)一個經(jīng)濟有效的稻田水分管理方法來增加水稻籽粒硒含量。研究表明,旱作較淹水栽培的土壤孔隙水中硒濃度更高,進而增加了水稻籽粒硒含量[8]。另一方面,硒的微生物轉(zhuǎn)化是將硒從一種形態(tài)變?yōu)榱硪环N形態(tài)的驅(qū)動力,近年來的研究表明,微生物能夠還原Se(Ⅳ)、Se(Ⅵ)→Se(0)[9],該過程可能降低了硒的生物有效性;而微生物對硒的氧化作用Se(0)→Se(Ⅳ)、Se(Ⅵ)報道較少,該過程有望增加硒的生物有效性。然而,到目前為止,稻田水分管理方式和根域相對微氧環(huán)境對土壤硒生物有效性及根際微生物對稻田硒循環(huán)的重要作用機制尚未闡明。

    近年來,各地研制并推廣了多種水分管理模式(如干濕交替、好氧灌溉),不僅能夠減少灌溉用水量,而且是保證水稻高產(chǎn)高效的重要措施[10]。水分管理很大程度上影響著稻田土壤的氧化還原電位和酸堿度,進而控制著稻田硒的化學轉(zhuǎn)化以及微生物轉(zhuǎn)化過程,最終影響水稻籽粒硒積累。基于以上背景,研究不同水分管理模式下水稻根際土壤中硒的形態(tài)轉(zhuǎn)化規(guī)律、水稻籽粒硒積累、水稻根際土壤細菌群落多樣性三者間的關(guān)系,揭示水分管理影響水稻吸收積累硒的作用機制,為從源頭上、生產(chǎn)上有效調(diào)控硒由稻田向農(nóng)產(chǎn)品傳遞以及滿足人類硒營養(yǎng)健康具有重要的理論和現(xiàn)實意義。

    1 材料與方法

    1.1 供試材料

    供試水稻品種為晶兩優(yōu)534。供試土壤采自重慶市北碚區(qū)西南大學試驗農(nóng)場,采用多點取樣法采集20~30 cm土層,土壤類型為普通紫色濕潤雛形土(中性紫色土),其基本理化性質(zhì)如下:pH 6.8,有機質(zhì)10.14 g·kg–1,堿解氮73.56 mg·kg–1,有效磷15.32 mg·kg–1,速效鉀171.3 mg·kg–1,全硒0.27 mg·kg–1。供試土壤經(jīng)自然風干后磨碎過篩,充分混勻用于盆栽試驗。

    1.2 試驗設(shè)計

    本試驗在西南大學資源環(huán)境學院植物營養(yǎng)學植物培養(yǎng)室內(nèi)進行(白天光照時間14 h,光強度260~350 μmol·m–2·s–1,白天溫度28℃,夜晚溫度20℃,相對濕度60%~70%)。盆栽試驗所用的培養(yǎng)缽為聚乙烯盆(直徑24 cm,高30 cm),每盆裝土9 kg。2019年6月5日施基肥,肥料用量分別為N 300 mg·kg–1,P2O5100 mg·kg–1,K2O 300 mg·kg–1,分別以尿素、過磷酸鈣和硫酸鉀作為肥源,同時所有處理添加亞硒酸鈉1.00 mg·kg–1與土壤充分混勻,加水平衡一周后進行水稻移栽,移栽60 d后進行追肥,追施尿素0.3 g·kg–1,硫酸鉀0.05 g·kg–1,每盆根際袋內(nèi)移栽2株,水稻移栽后即進行不同水分管理處理。本試驗設(shè)置3種水分管理處理,(1)淹水灌溉F(長期淹水,保持水層2 cm至作物收獲前2周);(2)好氧灌溉A(長期保持土壤濕潤,無明顯積水);(3)干濕交替灌溉AFA(淹水與灌溉交替進行,24 d為一個周期),每個處理重復3次。

    1.3 樣品采集與分析

    自水稻移栽20 d后,每隔12天測定一次土壤pH、Eh:使用pH、Eh計測定[11]。

    分別在水稻的分蘗期(40 d)、抽穗期(60 d)和成熟期(110 d)取根際袋內(nèi)的土壤,風干后過0.149 mm篩待用。測定土壤中不同硒形態(tài)的含量,分級提取依據(jù)瞿建國等[12]的逐級連續(xù)浸提方法。提取之后,各形態(tài)提取液需在95℃、6 mol·L–1鹽酸介質(zhì)中還原15 min后用原子熒光光譜儀(AFS,AFS-920,北京吉天儀器有限公司)測定。

    成熟期測定水稻根、莖、葉、籽粒中的硒含量:各部位進行分別粉碎,樣品用HNO3︰HClO4(4︰1)混合消煮,消化同Zhang等[13]的方法。消化待測液采用原子熒光光譜儀(AFS,AFS-920,北京吉天儀器有限公司)進行測定??瞻讟悠吠瑫r進行消化,硒回收率為95%~105%。

    收獲后測定水稻株高、有效穗數(shù)以及產(chǎn)量。

    1.4 土壤微生物多樣性分析

    水稻收獲前,取根際袋中的鮮土用于土壤微生物多樣性測定。1)DNA 抽提和PCR擴增:水稻根際土壤總DNA的抽提使用E.Z.N.A.?soil試劑盒(Omega Bio-tek,Norcross,GA,美國),該試劑盒可快速(60 min內(nèi)可處理1 g土壤)可靠地從各種土壤樣本中分離出高質(zhì)量的基因組DNA。用1%瓊脂糖凝膠電泳檢測DNA提取質(zhì)量后,用338F(5′-ACTCCTACGGGAGGCAGCAG-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)引物對V3~V4可變區(qū)進行PCR擴增,擴增程序為:95℃預變性3 min,27個循環(huán)(95℃變性30 s,55℃退火30 s,72℃延伸30 s),最后72℃延伸 10 min(PCR儀:ABI GeneAmp?9700型)。2)Illumina Miseq 測序:使用2%瓊脂糖凝膠回收PCR產(chǎn)物,純化后洗脫檢測。根據(jù)Illumina MiSeq 平臺(Illumina,San Diego,美國)標準操作規(guī)程構(gòu)建文庫,利用Illumina公司的Miseq PE300平臺進行測序(上海美吉生物醫(yī)藥科技有限公司)。Miseq測序得到的原始數(shù)據(jù)經(jīng)過拼接質(zhì)控后得到優(yōu)化數(shù)據(jù),區(qū)分樣本后進行OTU(Operational taxonomic unit,操作分類單元)聚類分析和物種分類學分析,可以對OTU進行多種多樣性指數(shù)分析,以及在各個分類水平上進行群落結(jié)構(gòu)的統(tǒng)計分析,進而完成樣本群落組成和系統(tǒng)發(fā)育信息的多元分析和差異顯著性檢驗等一系列深入的統(tǒng)計學和可視化分析。

    1.5 數(shù)據(jù)處理

    采用Microsoft Excel 2007軟件整理數(shù)據(jù)后,用SPSS 18.0 軟件進行數(shù)據(jù)的方差分析,用鄧肯(Duncan)檢驗法進行差異顯著性比較。應用Mothur軟件根據(jù)97%序列的相似度,將序列歸為多個OTU并生成稀釋曲線,計算覆蓋率(Coverage)、超(Chao)、香農(nóng)(Shannon)、香農(nóng)埃文(Shannoneven)指數(shù)。有關(guān)指數(shù)的算法如下[14-16]:

    式中,、chao、shannon、shannon分別表示覆蓋率、超、香農(nóng)、香農(nóng)埃文指數(shù)。1為僅含有一條序列的OTU數(shù)目;2為僅含有兩條序列的OTU數(shù)目;n為第個OTU所含的序列數(shù);obs為實際觀測到的OTU數(shù);為抽樣中出現(xiàn)的總序列數(shù)目;max為在物種豐富度相同時能夠達到的最大香農(nóng)指數(shù);為群落物種豐富度指數(shù);通常使用2、e等作為底數(shù)。

    2 結(jié) 果

    2.1 水分管理對土壤pH和Eh的影響

    水稻移栽40 d后,每隔12天測定一次水稻土壤的pH、Eh。如圖1所示,好氧和干濕交替灌溉較淹水灌溉一定程度上提高了土壤pH,且隨著處理時間的延長,好氧灌溉的土壤pH呈現(xiàn)小幅度的升高后降低的趨勢,總體變化幅度不大。在整個水稻生育期,淹水(6.50~6.73)和干濕交替灌溉(6.68~6.83)的土壤pH呈現(xiàn)弱酸性,好氧灌溉(7.19~7.28)的土壤則接近中性。

    對于土壤Eh而言,淹水和干濕交替灌溉在不同生育時期均顯著低于好氧灌溉(<0.05)。淹水灌溉的土壤處于還原狀態(tài),氧化還原電位處于–203~–261 mV,在移栽88 d之前變化幅度不大,88 d之后土壤Eh上升,這是由于在水稻收獲前曬田,使得土壤含水量下降所致。干濕交替的土壤氧化還原電位隨著水分管理方式的變化而改變,好氧灌溉的水稻土壤處于氧化狀態(tài),土壤Eh的變化范圍為163~253 mV,隨著處理時間的延長,土壤Eh有所下降。

    2.2 水分管理對土壤硒組分含量變化的影響

    通過逐級連續(xù)浸提的方法測定了水稻根際土壤中不同形態(tài)硒含量,結(jié)果表明,不同形態(tài)硒含量由高到低依次為F5(殘渣態(tài))、F3(酸溶態(tài))、F4(有機物結(jié)合態(tài))、F2(可交換態(tài))和F1(水溶態(tài))。其中,殘渣態(tài)占土壤全硒含量的52.11%~59.50%;酸溶態(tài)硒占土壤全硒含量的13.36%~16.12%;有機物結(jié)合態(tài)硒(包括碳酸鹽以及鐵錳氧化物結(jié)合態(tài))占土壤全硒含量的11.73%~13.63%;可交換態(tài)硒占根際土壤全硒含量的10.01%~14.76%;水溶態(tài)硒含量最低,占土壤全硒含量的1.65%~7.28%(圖2)。

    在水稻分蘗期、抽穗期、成熟期,水稻根際土壤水溶態(tài)硒、可交換態(tài)硒含量呈現(xiàn)A>AFA>F處理的趨勢,酸溶態(tài)硒、有機物結(jié)合態(tài)硒、殘渣態(tài)硒含量呈現(xiàn)F>AFA>A的趨勢。即好氧和干濕交替處理下,土壤有效態(tài)硒含量高于淹水處理。

    2.3 水分管理對水稻各器官硒含量變化的影響

    從圖3中可以看出,不同水分管理條件下水稻各器官的硒含量由高到低依次為根(1.38~2.22 mg·kg–1)、葉(0.55~0.85 mg·kg–1)、莖(0.53~0.61 mg·kg–1)、籽粒(0.15~0.53 mg·kg–1)。其中,水稻根中硒含量為水稻葉硒含量的2.26倍~2.59倍,是水稻莖的2.56倍~3.63倍,是水稻籽粒的4.21倍~9.01倍。A處理的水稻根中硒濃度最低,分別較F、AFA處理低38.33%、60.62%(<0.05)。三種不同水分管理條件對成熟期水稻莖中硒濃度的影響不顯著,好氧灌溉下的水稻葉硒含量顯著低于其他兩種處理。干濕交替灌溉的水稻籽粒硒含量最高,且顯著高于淹水灌溉,好氧灌溉的水稻籽粒硒含量較淹水高出0.16 mg·kg–1。

    2.4 水分管理對水稻各器官硒轉(zhuǎn)移系數(shù)的影響

    圖4呈現(xiàn)出不同水分管理條件對水稻硒轉(zhuǎn)移系數(shù)的影響,TCS/R、TCL/S、TCG/L分別表示水稻莖/根、葉/莖、籽粒/葉中硒濃度的比值,反映出水稻硒在各器官中的轉(zhuǎn)移狀況。在水稻成熟期,好氧灌溉的TCS/R最高(0.39),顯著高于淹水(0.29)及干濕交替(0.28)處理。好氧灌溉下處理水稻硒由莖向葉的轉(zhuǎn)移系數(shù)顯著低于其他處理,淹水灌溉下水稻硒由葉到籽粒的轉(zhuǎn)移系數(shù)顯著低于其他兩處理(圖4),且分別較干濕交替、好氧灌溉低70.55%、67.94%。

    2.5 水分管理對水稻產(chǎn)量的影響

    水稻收獲時期,三個處理之間的株高差異不顯著。干濕交替灌溉的水稻有效穗數(shù)最多,且顯著高于好氧灌溉,分別較好氧、淹水灌溉高出38.10%、19.05%。干濕交替灌溉產(chǎn)量顯著高于好氧灌溉,且較常規(guī)淹水灌溉提高了7.83%,較好氧灌溉提高13.51%(表1)。

    表1 不同水分管理下的水稻產(chǎn)量

    注:表中數(shù)據(jù)為平均值±標準誤差。同列不同小寫字母表示處理間差異顯著(<0.05)。下同。Note:The data in the table are of mean± standard deviation. Different lowercase letters in the same column indicate significant difference between treatments at<0.05. The same below.

    2.6 水分管理對水稻根際土壤細菌α多樣性的影響

    稀釋曲線可用來比較測序數(shù)據(jù)量不同的樣本中物種的豐富度、均一性或多樣性,也可用來說明樣本的測序數(shù)據(jù)量是否合理。如圖5所示,隨著取樣數(shù)的增加,稀釋曲線逐漸趨于平穩(wěn),OTU數(shù)目逐漸飽和,說明本試驗取樣合理,樣本數(shù)足夠大且能夠較為全面地反映出水稻根際土壤細菌群落狀況。

    經(jīng)過Mothur軟件分析,在OTU相似水平為0.97的情況下得到能夠反映水稻根際土壤細菌α多樣性的指數(shù),包括覆蓋率、超、香農(nóng)、香農(nóng)埃文指數(shù)。不同水分管理措施下,土壤細菌群落覆蓋率指數(shù)之間無顯著差異,即水分管理對細菌群落的覆蓋度不產(chǎn)生明顯影響。土壤菌群豐富度表現(xiàn)為F(4 556)>AFA(4 140)>A(3 761),三者之間存在顯著性差異,淹水處理下的根際土壤細菌豐富度最高。對于土壤細菌多樣性而言,淹水處理的香農(nóng)指數(shù)顯著高于其他兩處理,分別較A、AFA處理高出9.20%、3.14%。淹水處理的土壤菌群均勻度(0.847)最好,顯著高于好氧處理(0.803),干濕交替處理處于中間水平,但與淹水處理之間差異未達顯著水平(表2)。對于土壤細菌群落覆蓋度、多樣性、豐富度和均勻度而言,均呈現(xiàn)F>AFA>A的規(guī)律,即土壤由淹水變?yōu)闈駶櫊顟B(tài)時,水稻根際土壤細菌α多樣性呈下降趨勢。

    表2 不同水分管理條件下水稻根際土壤細菌多樣性指數(shù)

    2.7 水分管理對水稻根際土壤細菌群落結(jié)構(gòu)多樣性的影響

    不同水分管理對細菌在門水平上的物種組成影響如圖6所示,好氧灌溉下的根際土壤優(yōu)勢菌為變形菌門(41.27%~48.18%)、綠彎菌門(8.49%~10.84%)、擬桿菌門(8.23%~10.81%)、酸桿菌門(8.09%~8.20%)、Patescibacteria(5.72%~8.96%)、芽單胞菌門(6.20%~8.46%)6種細菌,占總菌數(shù)的88.71%~90.35%。干濕交替灌溉下的根際土壤優(yōu)勢菌為變形菌門(39.54%~46.41%)、綠彎菌門

    (10.45%~14.76%)、擬桿菌門(8.10%~9.63%)、酸桿菌(8.26%~10.01%)、Patescibacteria(6.15%~6.73%)、芽單胞菌門(3.46%~5.12%)、放線菌門(2.83%~4.05%)、厚壁菌門(2.36%~3.54%)8種細菌,占總菌數(shù)的89.97%~91.73%。淹水灌溉下的根際土壤優(yōu)勢菌為變形菌門(35.02%~43.73%)、綠彎菌門(9.98%~15.21%)、擬桿菌門(10.01%~10.77%)、酸桿菌門(6.41%~9.62%)、Patescibacteria(4.91%~7.59%)、芽單胞菌門(2.54%~3.60%)、放線菌門(2.77%~4.80%)、厚壁菌門(2.80%~3.81%)8種細菌,占總菌數(shù)的87.00%~88.03%。

    如圖7所示,對A、AFA、F三種不同水分管理方式下根際土壤樣本進行優(yōu)勢屬的組間差異顯著性檢驗分析,結(jié)果表明在屬水平,A處理下豐度最高,AFA處理次之,F(xiàn)處理最低,三者間存在顯著差異(<0.05)。的豐度在AFA處理下最高,F(xiàn)處理下最低。A處理下,豐度顯著高于其他兩處理。豐度在A處理下最高,F(xiàn)處理下豐度最低,且差異達極顯著水平(<0.01)。在A處理下豐度最低,豐度在F處理下顯著低于其他處理,豐度在F處理時最低,且與其他兩處理之間存在極顯著差異。

    水稻根際土壤豐度前5位的細菌菌綱與不同硒形態(tài)的環(huán)境因子關(guān)聯(lián)分析如圖8所示,圖中不同菌綱以實線表示,虛線表示土壤Eh、pH以及不同土壤硒形態(tài),環(huán)境因子箭頭間的夾角代表正(銳角)、負相關(guān)性(鈍角)。結(jié)果表明,土壤Eh、pH與土壤中有效態(tài)硒(水溶態(tài)、可交換態(tài)硒)含量呈顯著正相關(guān),且對于物種數(shù)據(jù)的影響程度(解釋量)較大。土壤中豐度最高的Gammaproteobacteria與土壤Eh、pH以及水溶態(tài)、可交換態(tài)、酸溶態(tài)硒濃度呈現(xiàn)顯著正相關(guān),與有機物結(jié)合態(tài)、殘渣態(tài)硒含量呈負相關(guān)。Deltaproteobacteria、Anaerolineae、Alphaproteobacteria、Bacteroidia豐度與土壤有機物結(jié)合態(tài)、殘渣態(tài)硒含量呈正相關(guān),與土壤Eh、pH以及土壤水溶態(tài)、可交換態(tài)、酸溶態(tài)硒含量呈負相關(guān)。

    3 討 論

    3.1 水分管理方式對水稻根際土壤硒形態(tài)和籽粒硒積累的影響

    土壤Eh和pH是影響稻田土壤硒有效性的關(guān)鍵因子[4]。在好氧和干濕交替灌溉下,土壤Eh和pH要高于淹水灌溉(圖1)。化學純培養(yǎng)體系得出,在酸性及還原條件下(7.5 < pE+ pH < 15,pE表示電子活度),亞硒酸鹽是土壤中硒主要的存在形態(tài),土壤中的鐵氧化物和氫氧化物對亞硒酸鹽有很強的吸附作用,降低其在土壤溶液中的濃度進而影響亞硒酸鹽的生物有效性[17];在堿性和氧化環(huán)境中(pE + pH >15),硒酸鹽是土壤中硒的主要存在形態(tài)[4],且不易被土壤膠體吸附,所以,硒酸鹽較亞硒酸鹽具有更高的植物有效性[7];在強還原條件下,硒還可能轉(zhuǎn)化為Se(-Ⅱ,0),這兩種硒形態(tài)不溶于水,所以對植物而言是無效的??茖W家在含硒量較高的淹水沉積物中發(fā)現(xiàn),當亞鐵氧化成三價鐵之后,Se(-Ⅱ,0)可部分轉(zhuǎn)化為Se(Ⅳ),當氧化還原電位超過200 mV時,亞硒酸鹽可緩慢氧化為硒酸鹽[4,18]。因此,土壤中硒的存在形態(tài)很大程度上取決于土壤氧化還原條件以及酸堿度,不同形態(tài)的硒可在微生物、土壤腐殖化、氧化還原作用下進行轉(zhuǎn)化,最終在土壤中維持一個相對平衡的狀態(tài)。

    本試驗中,不同形態(tài)硒含量由高到低依次為:殘渣態(tài)(52.11%~59.50%)、酸溶態(tài)(13.36%~16.12%)、有機物結(jié)合態(tài)(包括碳酸鹽以及鐵錳氧化物結(jié)合態(tài))(11.73%~13.63%)、可交換態(tài)(10.01%~14.76%)和水溶態(tài)(1.65%~7.28%),這與先前的研究[19-20]一致。好氧和干濕交替處理下土壤有效態(tài)硒(水溶態(tài)和可交換態(tài)硒)含量高于淹水處理,而淹水處理的酸溶態(tài)硒、有機物結(jié)合態(tài)硒、殘渣態(tài)硒含量則升高(圖2),說明淹水條件下土壤硒的生物有效性較低,這與先前的研究[8,21-22]一致。

    在淹水條件下,稻田土壤長期處于還原狀態(tài),土壤中高價態(tài)硒可被還原為Se(0,Ⅱ-),可溶性和可交換態(tài)硒逐漸轉(zhuǎn)化為結(jié)合態(tài)、殘渣態(tài)硒[4]。另一方面,長期的淹水處理使水稻形成了發(fā)達的通氣組織[23],在通氣組織的作用下,水稻根際處于一種相對氧化的狀態(tài),根際土壤溶液中的Fe2+、Mn2+被固定于根表形成鐵膜,將土壤中的硒更多地固定在水稻根表鐵膜中,降低其有效性[24]。綜上所述,淹水灌溉相較干濕交替、好氧灌溉而言,土壤硒的生物有效性較低(圖8)。

    就水稻籽粒硒含量而言,以干濕交替處理最高,表現(xiàn)為AFA>A>F的趨勢(圖3),這說明,土壤中水溶態(tài)硒和可交換態(tài)硒含量越高,則土壤硒生物有效性越高,相應地提高了水稻籽粒中硒含量。本研究結(jié)果與Li等[8]的相同,土壤中添加亞硒酸鈉和硒酸鈉后,好氧灌溉處理的籽粒硒含量分別為淹水處理的2倍和25倍,可見好氧栽培可使水稻籽粒中積累更多硒。此外,本研究發(fā)現(xiàn),三種水分管理中,干濕交替灌溉的水稻產(chǎn)量最高(表1),這與張慎鳳[25]的研究結(jié)果一致,且干濕交替灌溉可顯著提高稻米出糙率、精米率、透明度、淀粉黏滯譜的最高黏度和崩解值等稻米品質(zhì)指標。連同本研究結(jié)果,干濕交替灌溉可顯著提高稻米硒含量(圖3),在水稻生產(chǎn)中,可通過控制土壤水的落干程度,可以獲得高產(chǎn)、優(yōu)質(zhì)、節(jié)水的效果。但特別需要注意的是,干濕交替灌溉必須采取輕干濕交替灌溉(自然落干至土壤水勢-15 kPa)方可顯著提高產(chǎn)量,如果是重干濕交替灌溉(自然落干至土壤水勢-30 kPa)則降低了產(chǎn)量和品質(zhì)。

    3.2 不同水分管理方式下水稻根際土壤細菌區(qū)系和活性與水稻硒吸收的關(guān)系

    土壤硒微生物轉(zhuǎn)化也是硒生物地球化學循環(huán)的一個重要過程。本研究的高通量測序結(jié)果表明,水稻根際土壤優(yōu)勢菌為變形菌門、綠彎菌門、擬桿菌門、酸桿菌門、Patescibacteria、芽單胞菌門等,這些菌群也常被認為是稻田土壤中常見的優(yōu)勢微生物[26],不同水分管理顯著影響著上述菌群的豐度(圖6)。變形菌門是不同水分管理方式下水稻根際土壤中豐度最高的細菌,有研究表明,變形菌門中Alphaproteobacteria是一類具有多種生理特性的有機體,它與木質(zhì)素、芳香化合物的降解有關(guān)[27]。Betaproteobacteria與土壤凈碳礦化率呈正相關(guān),同時對土壤氮素循環(huán)起到關(guān)鍵作用[26]。Gammaproteobacteria能參與土壤固氮、溶磷過程,對植物根際鐵載體的形成有促進作用[28]。Deltaproteobacteria(由、、、、、組成)被證明在水稻生態(tài)系統(tǒng)中參與還原Fe、S化合物[29-30]。就屬水平而言,干濕交替和好氧灌溉中的細菌屬、、、、、等顯著高于淹水灌溉(圖7),而目前尚未見報道上述菌屬是否參與土壤硒的生物轉(zhuǎn)化以及Fe或S的氧化還原反應,有待進一步驗證。

    變形菌門是本試驗土壤中豐度最高的菌門,且Gammaproteobacteria豐度與土壤pH、Eh以及有效硒(水溶態(tài)、可交換態(tài)硒)含量之間呈顯著正相關(guān)(圖8),即干濕交替以及好氧灌溉下Gammaproteobacteria豐度的增加在一定意義上提高了土壤硒的生物有效性。有證據(jù)證明,Gammaproteobacteria可產(chǎn)生的根際鐵載體能與土壤中的金屬離子結(jié)合形成可溶性的金屬螯合物從而提高硒的生物有效性[31]。也有研究表明,添加高產(chǎn)鐵載體根際菌后,土壤中水溶態(tài)硒含量提高了2.50倍~7.85倍[32],即變形菌所產(chǎn)生的根際鐵載體通過螯合作用,使鐵氧化物釋放出硒,這也可能是干濕交替和好氧灌溉較淹水灌溉的土壤有效硒含量較高的原因。與變形菌門相比,其他菌門不直接參與鐵和硫的氧化還原過程,如綠彎菌門是一類生長在極端環(huán)境中,可參與土壤碳、氮循環(huán)的菌群[33]。擬桿菌門主要作用是水解和分解大分子物質(zhì)如纖維素和半纖維素[34]。芽單胞菌門、酸桿菌門被認為是寡營養(yǎng)細菌(Oligotrophic bacteria)中的一員,在貧營養(yǎng)條件下能夠分解復雜的土壤基質(zhì)[35]。本研究也發(fā)現(xiàn),不同水分管理條件對以上菌門(除變形菌門外)無顯著的影響(圖6)。

    4 結(jié) 論

    在整個水稻生育期內(nèi),干濕交替和好氧灌溉的土壤pH、Eh較淹水灌溉高,使得土壤有效態(tài)硒含量(水溶態(tài)硒和可交換態(tài)硒之和)高于淹水灌溉,從而增加了土壤硒的化學有效性。干濕交替的水分管理方式,既能達到節(jié)水的目的,而且可以顯著增加水稻籽粒硒含量。通過對根際土壤細菌多樣性進行分析,發(fā)現(xiàn)變形菌門是不同水分管理方式下水稻根際土壤中豐度最高的細菌,水分管理措施對其豐富度產(chǎn)生影響,干濕交替和好氧灌溉中變形菌門的豐度明顯高于淹水處理,且Gammaproteobacteria豐度的增加可能是土壤硒生物有效性增加的另一個重要原因。

    [1] Combs G F. Selenium in global food systems[J]. British Journal of Nutrition,2001,85(5):517—547.

    [2] Williams P N,Lombi E,Sun G X,et al. Selenium characterization in the global rice supply chain[J]. Environmental Science & Technology,2009,43(15):6024—6030.

    [3] Zhou X B,Shi W M,Yang L Z. Effect of foliar application of selenite on selenium accumulation and distribution in rice[J]. Acta Pedologica Sinica,2007,44(1):73—78. [周鑫斌,施衛(wèi)明,楊林章. 葉面噴硒對水稻籽粒硒富集及分布的影響[J]. 土壤學報,2007,44(1):73—78.]

    [4] Masscheleyn P H,Delaune R D,Patrick W H. Transformations of selenium as affected by sediment oxidation-reduction potential and pH[J]. Environmental Science & Technology,1990,24(1):91—96.

    [5] Lyubenova L,Sabodash X,Schr?der P,et al. Selenium species in the roots and shoots of chickpea plants treated with different concentrations of sodium selenite[J]. Environmental Science and Pollution Research,2015,22(21):16978—16986.

    [6] Shahid M,Niazi N K,Khalid S,et al. A critical review of selenium biogeochemical behavior in soil-plant system with an inference to human health[J]. Environmental Pollution,2018,234:915—934.

    [7] Neal R H,Sposito G. Selenate adsorption on alluvial soils[J]. Soil Science Society of America Journal,1989,53(1):70—74.

    [8] Li H F,Lombi E,Stroud J L,et al. Selenium speciation in soil and rice:Influence of water management and Se fertilization[J]. Journal of Agricultural and Food Chemistry,2010,58(22):11837—11843.

    [9] Bao Q L,Xiao K Q,Chen Z,et al. Methane production and methanogenic archaeal communities in two types of paddy soil amended with different amounts of rice straw[J]. FEMS Microbiology Ecology,2014,88(2):372—385.

    [10] Mostafazadehfard F B,Jafari F,Mousavi S,et al. Effects of irrigation water management on yield and water use efficiency of rice in cracked paddy soils[J]. Australian Journal of Crop Science,2010,4(3):136—141.

    [11] Lu R K. Analytical methods for soil and agro-chemistry [M]. China Agricultural Science and Technology Press,2000:125—233. [魯如坤. 土壤農(nóng)業(yè)化學分析方法[M]. 北京:中國農(nóng)業(yè)科學技術(shù)出版社,2000:125—233.]

    [12] Qu J G,Xu B X,Gong S C. Sequential extraction techniques for determination of selenium speciation in soils and sediments[J]. Environmental Chemistry,1997,16(3):277—283. [瞿建國,徐伯興,龔書椿. 連續(xù)浸提技術(shù)測定土壤和沉積物中硒的形態(tài)[J]. 環(huán)境化學,1997,16(3):277—283.]

    [13] Zhang L H,Shi W M,Wang X C. Difference in selenite absorption between high- and low-selenium rice cultivars and its mechanism[J]. Plant and Soil,2006,282(1/2):183—193.

    [14] Chao A N,Yang M C K. Stopping rules and estimation for recapture debugging with unequal failure rates[J]. Biometrika,1993,80(1):193—201.

    [15] Faith D P. Conservation evaluation and phylogenetic diversity[J]. Biological Conservation,1992,61(1):1—10.

    [16] Hill M O. Diversity and evenness:A unifying notation and its consequences[J]. Ecology,1973,54(2):427—432.

    [17] Barrow N J,Whelan B R. Testing a mechanistic model. VII. The effects of pH and of electrolyte on the reaction of selenite and selenate with a soil[J]. Journal of Soil Science,1989,40(1):17—28.

    [18] Jayaweera G R,Biggar J W. Role of redox potential in chemical transformations of selenium in soils[J]. Soil Science Society of America Journal,1996,60(4):1056—1063.

    [19] Zhang D,Dong T,Ye J,et al. Selenium accumulation in wheat(L.)as affected by coapplication of either selenite or selenate with phosphorus[J]. Soil Science and Plant Nutrition,2017,63(1):37—44.

    [20] Bajaj M,Eiche E,Neumann T,et al. Hazardous concentrations of selenium in soil and groundwater in North-West India[J]. Journal of Hazardous Materials,2011,189(3):640—646.

    [21] Funwie A V. Effect of soil properties on availability and mobility of selenium[D]. Faculty of Bioscience Engineering,University of Gent,Belgium,2012:1—65.

    [22] Peng Q,Wang M K,Cui Z W,et al. Assessment of bioavailability of selenium in different plant-soil systems by diffusive gradients in thin-films(DGT)[J]. Environmental Pollution,2017,225:637—643.

    [23] Armstrong W. The oxidising activity of roots in waterlogged soils[J]. Physiologia Plantarum,1967,20(4):920—926.

    [24] Shiba H,Daimon H. Histological observation of secondary aerenchyma formed immediately after flooding inand[J]. Plant and Soil,2003,255(1):209—215.

    [25] Zhang S F. Effect of alternate wetting and drying on the growth and development,grain yield and quality of rice[D]. Yangzhou:Yangzhou University,2009. [張慎鳳. 干濕交替灌溉對水稻生長發(fā)育、產(chǎn)量與品質(zhì)的影響[D].揚州:揚州大學,2009.]

    [26] Huang Q,Wang J L,Wang C,et al. The 19-years inorganic fertilization increased bacterial diversity and altered bacterial community composition and potential functions in a paddy soil[J]. Applied Soil Ecology,2019,144:60—67.

    [27] Kersters K,de Vos P,Gillis M,et al. Introduction to the proteobacteria[M]//The Prokaryotes. New York,NY:Springer New York,2006:3—37.

    [28] Weller D M. Pseudomonas biocontrol agents of soilborne pathogens:Looking back over 30 years[J]. Phytopathology,2007,97(2):250—256.

    [29] Muyzer G,Stams A J M. The ecology and biotechnology of sulphate-reducing bacteria[J]. Nature Reviews Microbiology,2008,6(6):441—454.

    [30] Ding L J,Su J Q,Xu H J,et al. Long-term nitrogen fertilization of paddy soil shifts iron-reducing microbial community revealed by RNA-13C-acetate probing coupled with pyrosequencing[J]. The ISME Journal,2015,9(3):721—734.

    [31] Karag?z K,Ate? F,Karag?z H,et al. Characterization of plant growth-promoting traits of bacteria isolated from the rhizosphere of grapevine grown in alkaline and acidic soils[J]. European Journal of Soil Biology,2012,50:144—150.

    [32] Long Y C,Chen X,Zhou S Q. Isolation,identification and assessment on selenium biofortification of siderophore-producing rhizobacteria[J]. Current Biotechnology,2017(5):402—408. [龍云川,陳軒,周少奇. 高產(chǎn)鐵載體根際菌的篩選鑒定及硒活化特性評價[J]. 生物技術(shù)進展,2017(5):402—408.]

    [33] Kragelund C,Caterina L,Borger A,et al. Identity,abundance and ecophysiology of filamentous Chloroflexi species present in activated sludge treatment plants[J]. FEMS Microbiology Ecology,2007,59(3):671—682.

    [34] Wegner C E,Liesack W. Microbial community dynamics during the early stages of plant polymer breakdown in paddy soil[J]. Environmental Microbiology,2016,18(9):2825—2842.

    [35] Fierer N,Jackson R B. The diversity and biogeography of soil bacterial communities[J]. Proceedings of the National Academy of Sciences of the United States of America,2006,103(3):626—631.

    Effects of Water Management on Selenium Accumulation in Rice Grains and Bacterial Community Diversity in Rhizosphere Soil

    WANG Ruixin, YANG Jing, FANG Zheng, ZHANG Siyu, ZHOU Xinbin??

    (College of Resources and Environment, Southwest University, Chongqing 400716, China)

    【Objective】For this paper a pot experiment was conducted to evaluate effects of water management varying in pattern on composition of soil selenium, bacterial community diversity in rhizosphere soil and Se accumulation in grains of rice, and further to explore mechanism of water management affecting the absorption and accumulation of selenium in rice, so as to provide considerable theoretical and practical support to effective regulation of selenium transfer from rice field to agricultural products, to meet the demand of human for selenium nutrient and health. 【Method】The experiment was designed to have three treatments in water management treatments, that is, (I) flooded irrigation (F); (II) aerobic irrigation (A); and (III) alternation of flooded and aerobic irrigation (AFA), and three replicates for each treatment. Starting from D20 (20 days after rice transplanting), soil pH and Eh were measured every 12 days. The selenium in the rhizosphere soil was fractionated at the tillering stage, heading stage and maturity stage of the crop. Content of selenium in rice roots, stems, leaves and grains as well as yield of the crop were measured at maturity. Before rice harvest, the fresh soil in the rhizosphere bag was retrieved and placed in a bag, which was then sealed and stored at low temperature (–80℃) for high-throughput sequencing analysis. 【Result】Results show that at all the rice growth stages, Treatments A and AFA were higher than Treatment F in soil pH and Eh to a certain extent and in content of soluble and exchangeable selenium as well, thereby in availability of soil selenium. In terms of selenium content, the organs of a rice plant exhibited an order of root (1.38–2.22 mg·kg–1) > leaf (0.55–0.85 mg·kg–1) > stem (0.53–0.61 mg·kg–1) > grain (0.15– 0.53 mg·kg–1). Treatment AFA was the highest in selenium content in rice grains, while Treatment F was the lowest. The difference between the two was obvious. Treatment AFA was significantly or 7.83% higher than Treatment F and 13.51% higher than Treatment A in grain yield. The dominant species of bacteria in the rhizosphere soil of rice were Proteobacteria, Chloroflexi, Bacteroidetes, Acidobacteria, Patescibacteria and Gemmatimonadetes, with Proteobacteria being the most abundant one in rice rhizosphere soil in all the treatments. Obviously water management had a significant impact on richness of the bacteria. Treatments AFA and A were significantly higher than Treatment F in richness of Proteobacteria. On the level of class, abundance of Gammaproteobacteria was positively related to content of available Se in the soil, so the increase in the former may be deemed as an important reason for higher bioavailability of selenium in the soil. 【Conclusion】 Comprehensive analysis shows that alternation of flooded and aerobic irrigation can not only improve yield and selenium content of rice grains, but also save water for irrigation as compared to normal flooding irrigation. Therefore it can be recommended as preferential water management method in rice production.

    Rice; Selenium; Alternation of flooded and aerobic irrigation; Rhizosphere soil; Bacterial diversity

    S143.7+1

    A

    10.11766/trxb202004090050

    王瑞昕,楊靜,方正,張思宇,周鑫斌. 水分管理對水稻籽粒硒積累及根際土壤細菌群落多樣性的影響[J]. 土壤學報,2021,58(6):1574–1584.

    WANG Ruixin,YANG Jing,F(xiàn)ANG Zheng,ZHANG Siyu,ZHOU Xinbin. Effects of Water Management on Selenium Accumulation in Rice Grains and Bacterial Community Diversity in Rhizosphere Soil[J]. Acta Pedologica Sinica,2021,58(6):1574–1584.

    *國家自然科學基金項目(31672238)和重慶市基礎(chǔ)科學與前沿技術(shù)研究項目(cstc2017jcyjAX0425)資助 Supported by the National Natural Science Foundation of China(No. 31672238)and the Research Program of Basic Research and Frontier Technology of Chongqing City in China(No. cstc2017jcyjAX0425)

    Corresponding author,E-mail:zxbissas@163.com

    王瑞昕(1996—),女,山西長治人,碩士研究生,主要從事植物硒營養(yǎng)生理研究。E-mail:17835424992@163.com

    2020–04–09;

    2020–06–28;

    2020–09–22

    (責任編輯:陳榮府)

    猜你喜歡
    菌門根際籽粒
    籽粒莧的飼用價值和高產(chǎn)栽培技術(shù)
    根際微生物對植物與土壤交互調(diào)控的研究進展
    籽粒莧的特性和種植技術(shù)
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    饑餓與重攝食對河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進展
    妊娠期糖尿病腸道菌群變化及臨床價值分析
    玉米機械脫粒籽粒含水量與破碎率的相關(guān)研究
    黃花蒿葉水提物對三七根際尖孢鐮刀菌生長的抑制作用
    商麥1619 籽粒灌漿的特性
    久久99热6这里只有精品| 亚洲成色77777| 亚洲精品乱久久久久久| 色视频www国产| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产日韩一区二区| 中文资源天堂在线| 国产色爽女视频免费观看| 插阴视频在线观看视频| 各种免费的搞黄视频| 肉色欧美久久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 国产精品一区二区三区四区免费观看| 国产精品一区二区在线不卡| 成人综合一区亚洲| 十分钟在线观看高清视频www | 美女内射精品一级片tv| 欧美高清成人免费视频www| 又粗又硬又长又爽又黄的视频| 久久久久久久久久久丰满| 深夜a级毛片| 少妇 在线观看| 国产女主播在线喷水免费视频网站| 18禁在线无遮挡免费观看视频| 色5月婷婷丁香| 一级毛片电影观看| 夜夜爽夜夜爽视频| 女人精品久久久久毛片| 亚洲美女搞黄在线观看| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 欧美日韩在线观看h| 成年人午夜在线观看视频| 久久久久久久久久人人人人人人| 丰满人妻一区二区三区视频av| 国产爽快片一区二区三区| 国产精品无大码| 国产精品一二三区在线看| 午夜福利网站1000一区二区三区| 好男人视频免费观看在线| 国产精品99久久久久久久久| 国产高清国产精品国产三级| 午夜免费观看性视频| 国产男人的电影天堂91| 乱码一卡2卡4卡精品| 少妇人妻 视频| 成年av动漫网址| 又大又黄又爽视频免费| 香蕉精品网在线| 免费观看a级毛片全部| 亚洲av综合色区一区| 黄片无遮挡物在线观看| 韩国av在线不卡| freevideosex欧美| 99热全是精品| 亚洲真实伦在线观看| 成人午夜精彩视频在线观看| 蜜臀久久99精品久久宅男| 亚洲人与动物交配视频| 成人毛片60女人毛片免费| 国产一区亚洲一区在线观看| 人人澡人人妻人| 99久久精品国产国产毛片| 久久99热6这里只有精品| 97在线人人人人妻| 另类亚洲欧美激情| 欧美人与善性xxx| 天堂中文最新版在线下载| 男女免费视频国产| 婷婷色综合大香蕉| 少妇的逼好多水| 啦啦啦中文免费视频观看日本| 国产亚洲精品久久久com| 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 美女脱内裤让男人舔精品视频| 日韩欧美精品免费久久| 黄色日韩在线| 午夜91福利影院| 精品熟女少妇av免费看| 亚洲av日韩在线播放| av网站免费在线观看视频| 欧美激情国产日韩精品一区| 久久久久久久久久久丰满| 亚洲国产精品成人久久小说| 我要看日韩黄色一级片| 五月开心婷婷网| 日本欧美国产在线视频| 黄色怎么调成土黄色| av天堂久久9| 在线播放无遮挡| 啦啦啦中文免费视频观看日本| 人人妻人人爽人人添夜夜欢视频 | 色94色欧美一区二区| 日韩一区二区视频免费看| 久久精品国产亚洲网站| 国产成人免费无遮挡视频| 日韩精品有码人妻一区| 午夜免费观看性视频| 久久精品久久精品一区二区三区| 老司机影院毛片| 欧美成人午夜免费资源| 国产精品国产三级国产av玫瑰| 一级毛片aaaaaa免费看小| 国产精品久久久久久av不卡| 肉色欧美久久久久久久蜜桃| 欧美xxxx性猛交bbbb| 亚洲精品日韩av片在线观看| 中文字幕人妻熟人妻熟丝袜美| 免费久久久久久久精品成人欧美视频 | 成年人免费黄色播放视频 | 国产亚洲精品久久久com| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 亚洲精品成人av观看孕妇| 赤兔流量卡办理| 久久久亚洲精品成人影院| 大片电影免费在线观看免费| 亚洲av免费高清在线观看| 在线观看一区二区三区激情| 三级国产精品片| 欧美日韩国产mv在线观看视频| 亚洲内射少妇av| 国产免费视频播放在线视频| 日韩精品有码人妻一区| 一级爰片在线观看| 校园人妻丝袜中文字幕| 中文字幕人妻丝袜制服| 欧美3d第一页| 爱豆传媒免费全集在线观看| 亚洲,欧美,日韩| 热99国产精品久久久久久7| 精华霜和精华液先用哪个| 成人综合一区亚洲| 国产免费视频播放在线视频| 成人毛片a级毛片在线播放| 亚洲美女搞黄在线观看| 男女国产视频网站| 最近中文字幕2019免费版| 99热这里只有是精品在线观看| 国产精品国产三级国产av玫瑰| 99re6热这里在线精品视频| 国产精品嫩草影院av在线观看| 国产 一区精品| 久久久久久人妻| 欧美另类一区| 精品国产露脸久久av麻豆| 欧美3d第一页| 免费久久久久久久精品成人欧美视频 | 亚洲色图综合在线观看| 亚洲av欧美aⅴ国产| 亚洲图色成人| 一区二区三区免费毛片| 欧美日韩一区二区视频在线观看视频在线| 国产极品天堂在线| 18禁裸乳无遮挡动漫免费视频| 久久精品久久精品一区二区三区| 亚洲精品日本国产第一区| 国产精品一区www在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 夜夜看夜夜爽夜夜摸| 最新的欧美精品一区二区| 九九久久精品国产亚洲av麻豆| 国产精品欧美亚洲77777| 国产成人免费观看mmmm| 美女xxoo啪啪120秒动态图| 国产成人免费无遮挡视频| 欧美成人精品欧美一级黄| 伦精品一区二区三区| 亚洲精品亚洲一区二区| 欧美另类一区| 熟女人妻精品中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 日韩精品有码人妻一区| 一边亲一边摸免费视频| 日韩欧美精品免费久久| 伊人久久国产一区二区| 久久鲁丝午夜福利片| 国产成人精品无人区| 日本色播在线视频| 亚洲欧洲日产国产| 中文字幕av电影在线播放| 精品国产国语对白av| 下体分泌物呈黄色| 国产91av在线免费观看| 在线观看人妻少妇| 国产日韩欧美在线精品| 欧美激情极品国产一区二区三区 | 人人澡人人妻人| 人妻系列 视频| 纯流量卡能插随身wifi吗| 久久这里有精品视频免费| 99视频精品全部免费 在线| 日本vs欧美在线观看视频 | 国产男女超爽视频在线观看| 久久精品夜色国产| 日本欧美视频一区| 久久久久久久久久久免费av| 夫妻性生交免费视频一级片| 大陆偷拍与自拍| 欧美日本中文国产一区发布| 国产有黄有色有爽视频| 美女主播在线视频| 色视频在线一区二区三区| 国产精品久久久久久精品电影小说| 青春草视频在线免费观看| 我要看日韩黄色一级片| 99久久人妻综合| 高清不卡的av网站| 国产极品粉嫩免费观看在线 | av免费在线看不卡| 黑丝袜美女国产一区| 国产日韩欧美亚洲二区| 国产日韩欧美亚洲二区| 精品国产国语对白av| 人妻系列 视频| 人人妻人人添人人爽欧美一区卜| 国产一区二区三区av在线| 欧美+日韩+精品| 18禁动态无遮挡网站| 一区二区三区四区激情视频| 日韩中字成人| 性高湖久久久久久久久免费观看| 亚洲美女搞黄在线观看| 性色avwww在线观看| 十分钟在线观看高清视频www | 中文天堂在线官网| 少妇被粗大的猛进出69影院 | 一本大道久久a久久精品| 国产成人精品婷婷| 亚洲综合精品二区| 免费看av在线观看网站| 久久精品熟女亚洲av麻豆精品| 国产美女午夜福利| 秋霞在线观看毛片| 极品教师在线视频| 免费大片黄手机在线观看| 日日摸夜夜添夜夜添av毛片| 日韩欧美 国产精品| 最新中文字幕久久久久| 亚洲欧美日韩卡通动漫| 国产无遮挡羞羞视频在线观看| 国产精品伦人一区二区| 日日摸夜夜添夜夜添av毛片| 这个男人来自地球电影免费观看 | 美女cb高潮喷水在线观看| 欧美丝袜亚洲另类| 久久毛片免费看一区二区三区| 久久免费观看电影| 亚洲av福利一区| 高清欧美精品videossex| 高清在线视频一区二区三区| 黄色配什么色好看| 国产真实伦视频高清在线观看| 久久久久精品久久久久真实原创| 国产精品三级大全| 日本欧美国产在线视频| 国产亚洲av片在线观看秒播厂| 国产精品麻豆人妻色哟哟久久| 夫妻性生交免费视频一级片| 欧美bdsm另类| 久久6这里有精品| 午夜福利影视在线免费观看| 一级黄片播放器| 亚洲国产精品国产精品| 国产毛片在线视频| 日韩av在线免费看完整版不卡| 亚洲一级一片aⅴ在线观看| 久久99热这里只频精品6学生| 一级毛片我不卡| 成人二区视频| 色5月婷婷丁香| 国产精品国产三级专区第一集| 欧美国产精品一级二级三级 | 人妻 亚洲 视频| 人人妻人人澡人人爽人人夜夜| 在线观看免费日韩欧美大片 | 欧美日韩精品成人综合77777| 国产精品一区二区在线不卡| 欧美日韩国产mv在线观看视频| 日本91视频免费播放| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产欧美在线一区| 国产 精品1| 日本欧美视频一区| 久久久久久久久久人人人人人人| 国产熟女午夜一区二区三区 | 午夜福利网站1000一区二区三区| 日日撸夜夜添| 一级毛片我不卡| 国产美女午夜福利| av专区在线播放| 国内少妇人妻偷人精品xxx网站| 最新中文字幕久久久久| 美女cb高潮喷水在线观看| 久久久精品免费免费高清| 深夜a级毛片| av不卡在线播放| 久热这里只有精品99| 3wmmmm亚洲av在线观看| 在线 av 中文字幕| 亚洲精品国产av蜜桃| 秋霞伦理黄片| 亚洲精品视频女| 观看美女的网站| 97在线人人人人妻| 精品人妻偷拍中文字幕| 午夜免费鲁丝| 久久精品国产鲁丝片午夜精品| 国产av码专区亚洲av| 精品午夜福利在线看| 精品人妻熟女av久视频| 97超碰精品成人国产| 成人黄色视频免费在线看| av天堂久久9| 最近2019中文字幕mv第一页| 日韩电影二区| 男女啪啪激烈高潮av片| 视频中文字幕在线观看| 国产高清三级在线| 在线观看av片永久免费下载| 免费看日本二区| 又黄又爽又刺激的免费视频.| 欧美精品一区二区大全| 久久综合国产亚洲精品| 热99国产精品久久久久久7| 嘟嘟电影网在线观看| 少妇高潮的动态图| 日韩大片免费观看网站| 国产亚洲一区二区精品| 久久久久久久久久成人| 久久精品久久久久久噜噜老黄| 黄色怎么调成土黄色| av线在线观看网站| 国产精品一区二区性色av| 国产高清国产精品国产三级| 久久精品夜色国产| 少妇被粗大的猛进出69影院 | 老女人水多毛片| 日本91视频免费播放| 久久久久久久久久成人| 国内揄拍国产精品人妻在线| 国产在视频线精品| 亚洲无线观看免费| 少妇人妻久久综合中文| 亚洲国产色片| 精品人妻熟女毛片av久久网站| 久久久久久久久大av| 免费在线观看成人毛片| 十八禁网站网址无遮挡 | 性色av一级| 成人亚洲精品一区在线观看| 国产在线免费精品| 亚洲怡红院男人天堂| 简卡轻食公司| 人妻 亚洲 视频| 蜜臀久久99精品久久宅男| 免费观看无遮挡的男女| 日韩 亚洲 欧美在线| 欧美+日韩+精品| 少妇人妻 视频| 在现免费观看毛片| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片 | 人人妻人人澡人人看| 我要看黄色一级片免费的| 18+在线观看网站| 99久久中文字幕三级久久日本| 妹子高潮喷水视频| 国产黄片美女视频| 日日撸夜夜添| 国产乱来视频区| 日韩av在线免费看完整版不卡| 自拍欧美九色日韩亚洲蝌蚪91 | 在线看a的网站| 一级a做视频免费观看| 一级av片app| 最近中文字幕2019免费版| 夫妻午夜视频| 老司机影院成人| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 免费人妻精品一区二区三区视频| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久国产精品麻豆| 亚洲一级一片aⅴ在线观看| 欧美xxxx性猛交bbbb| 国产高清有码在线观看视频| 久久久久久久久久久丰满| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久av不卡| 一级爰片在线观看| 亚洲精品,欧美精品| 成人无遮挡网站| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| 精品少妇内射三级| 久热这里只有精品99| 国产伦精品一区二区三区视频9| 2021少妇久久久久久久久久久| 看免费成人av毛片| 日韩av免费高清视频| 日韩一区二区三区影片| 三级经典国产精品| 大陆偷拍与自拍| 女性生殖器流出的白浆| 亚洲自偷自拍三级| 超碰97精品在线观看| 99视频精品全部免费 在线| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 狂野欧美激情性bbbbbb| www.av在线官网国产| 久久久久人妻精品一区果冻| 精品亚洲乱码少妇综合久久| 亚洲精品aⅴ在线观看| 18禁在线播放成人免费| 亚洲不卡免费看| 亚洲精品乱码久久久久久按摩| 美女福利国产在线| 成人免费观看视频高清| 国产成人免费观看mmmm| 国产精品久久久久久精品古装| 国产精品一区二区性色av| 美女内射精品一级片tv| 人人妻人人爽人人添夜夜欢视频 | 噜噜噜噜噜久久久久久91| 欧美成人午夜免费资源| 欧美精品一区二区免费开放| 乱系列少妇在线播放| 日本欧美视频一区| 看非洲黑人一级黄片| 亚洲精品456在线播放app| av福利片在线| 日日啪夜夜爽| 五月玫瑰六月丁香| 成年女人在线观看亚洲视频| 国产亚洲最大av| 欧美老熟妇乱子伦牲交| 国产高清有码在线观看视频| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 国产熟女欧美一区二区| 亚洲国产精品成人久久小说| 国产黄色免费在线视频| 久热久热在线精品观看| 精品一区二区免费观看| 欧美97在线视频| 欧美另类一区| 国产探花极品一区二区| 欧美日韩视频高清一区二区三区二| 成年女人在线观看亚洲视频| 亚洲综合色惰| 新久久久久国产一级毛片| 精华霜和精华液先用哪个| 51国产日韩欧美| 老司机影院毛片| 少妇裸体淫交视频免费看高清| 三上悠亚av全集在线观看 | 夫妻午夜视频| 国产精品久久久久久久久免| 色5月婷婷丁香| 久久人人爽人人爽人人片va| 嘟嘟电影网在线观看| 日韩人妻高清精品专区| 春色校园在线视频观看| 看十八女毛片水多多多| 国产综合精华液| 日韩欧美一区视频在线观看 | 亚洲图色成人| 伦理电影免费视频| 插阴视频在线观看视频| 九九爱精品视频在线观看| 亚洲精品日韩在线中文字幕| 曰老女人黄片| 99精国产麻豆久久婷婷| 视频中文字幕在线观看| h视频一区二区三区| 亚洲综合色惰| 国产欧美日韩一区二区三区在线 | 精品酒店卫生间| 亚洲,一卡二卡三卡| 中国美白少妇内射xxxbb| 久久ye,这里只有精品| 国产一区二区在线观看日韩| 少妇人妻精品综合一区二区| 国产免费一级a男人的天堂| 亚洲四区av| 久久久久久久久久人人人人人人| 国产免费视频播放在线视频| 人妻制服诱惑在线中文字幕| 日韩强制内射视频| 日韩精品免费视频一区二区三区 | 99久久综合免费| 国产伦在线观看视频一区| 欧美少妇被猛烈插入视频| 又爽又黄a免费视频| 春色校园在线视频观看| 色网站视频免费| av在线app专区| 草草在线视频免费看| 99热这里只有是精品在线观看| 91午夜精品亚洲一区二区三区| 我要看黄色一级片免费的| 久久国产精品大桥未久av | 99久久精品一区二区三区| 国产亚洲91精品色在线| 欧美精品高潮呻吟av久久| 日韩大片免费观看网站| 十八禁高潮呻吟视频 | 性色av一级| 日本与韩国留学比较| 久久精品久久久久久久性| 亚洲真实伦在线观看| 三级国产精品片| 亚洲精品视频女| 水蜜桃什么品种好| 欧美精品一区二区大全| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产av蜜桃| 如日韩欧美国产精品一区二区三区 | av视频免费观看在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | av免费在线看不卡| av黄色大香蕉| 国产精品秋霞免费鲁丝片| 狠狠精品人妻久久久久久综合| 五月伊人婷婷丁香| 免费看av在线观看网站| 日本av免费视频播放| 久久99热6这里只有精品| 国产黄频视频在线观看| 性色av一级| 在线观看免费日韩欧美大片 | 韩国高清视频一区二区三区| 国产成人精品久久久久久| 六月丁香七月| 免费少妇av软件| 国产欧美另类精品又又久久亚洲欧美| 久热久热在线精品观看| 婷婷色麻豆天堂久久| 亚洲精品第二区| 蜜臀久久99精品久久宅男| 午夜激情福利司机影院| 久久免费观看电影| 少妇人妻精品综合一区二区| 精品酒店卫生间| 久久精品久久精品一区二区三区| 亚洲久久久国产精品| 精品少妇久久久久久888优播| 亚洲久久久国产精品| 午夜免费鲁丝| 色5月婷婷丁香| 日本爱情动作片www.在线观看| 久久亚洲国产成人精品v| 成人漫画全彩无遮挡| 国产精品人妻久久久久久| 少妇的逼好多水| 成年女人在线观看亚洲视频| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 欧美精品国产亚洲| 十八禁高潮呻吟视频 | 午夜精品国产一区二区电影| 国产一区二区三区av在线| 热re99久久精品国产66热6| 亚洲熟女精品中文字幕| 极品教师在线视频| 国产欧美日韩综合在线一区二区 | 亚洲欧美日韩东京热| 午夜福利视频精品| 欧美最新免费一区二区三区| 亚洲国产最新在线播放| 大陆偷拍与自拍| av不卡在线播放| 成人毛片60女人毛片免费| 曰老女人黄片| 男女边摸边吃奶| 久久久久国产精品人妻一区二区| 午夜日本视频在线| 久久精品国产亚洲网站| 午夜日本视频在线| 青春草国产在线视频| 欧美精品一区二区免费开放| 丝袜在线中文字幕| 日韩伦理黄色片| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 一级毛片aaaaaa免费看小| 久久影院123| 人人妻人人看人人澡| 我要看日韩黄色一级片| 欧美激情国产日韩精品一区| 成年人免费黄色播放视频 | 建设人人有责人人尽责人人享有的| 久久99精品国语久久久| 一级av片app| 一区二区av电影网| 在线观看一区二区三区激情| 日本猛色少妇xxxxx猛交久久| 国产精品国产av在线观看| 午夜免费鲁丝| av天堂中文字幕网| 亚洲va在线va天堂va国产| 天堂中文最新版在线下载| 国产精品一区www在线观看| 婷婷色麻豆天堂久久| 91精品国产九色| 精品久久久久久电影网| 在线观看一区二区三区激情| 97在线视频观看|