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

    轉(zhuǎn)基因水稻秸稈還田對土壤硝化反硝化微生物群落的影響*

    2018-01-04 01:43:34王沛譞宋亞娜
    關鍵詞:水稻深度

    王沛譞, 徐 焱, 宋亞娜

    ?

    轉(zhuǎn)基因水稻秸稈還田對土壤硝化反硝化微生物群落的影響*

    王沛譞1,3, 徐 焱2, 宋亞娜1**

    (1. 福建省農(nóng)業(yè)科學院生物技術研究所 福州 350003; 2. 福州時代中學 福州 350007; 3. 福建師范大學附中福州 350007)

    轉(zhuǎn)基因作物可能通過根系分泌物和植株殘體組成的改變及外源基因的轉(zhuǎn)移釋放令土壤微生物群落產(chǎn)生變化, 影響土壤微生物的生態(tài)功能。氨氧化細菌和反硝化細菌是驅(qū)動土壤硝化和反硝化過程的關鍵微生物, 其群落結(jié)構的變化直接關系土壤氮素的轉(zhuǎn)化與利用。本研究利用熒光定量PCR和PCR-DGGE技術分析了轉(zhuǎn)雙價抗蟲基因水稻‘Kf8’秸稈還田降解過程中, 土壤氨氧化細菌和反硝化細菌群落豐度與組成的變化, 探討轉(zhuǎn)基因水稻是否存在影響稻田土壤氮素轉(zhuǎn)化與N2O排放的可能。結(jié)果顯示: 無論是氨氧化細菌基因還是反硝化細菌基因, 其豐度在轉(zhuǎn)基因水稻‘Kf8’與非轉(zhuǎn)基因水稻‘Mh86’的秸稈還田土壤中都沒有顯著差異; 轉(zhuǎn)基因水稻‘Kf8’和非轉(zhuǎn)基因水稻‘Mh86’秸稈還田降解過程中0~10 cm土層中的基因豐度均顯著高于10~20 cm及20~30 cm土層(<0.05); 各深度土層中的基因豐度均存在隨秸稈還田時間延長而增加的趨勢。水稻秸稈還田降解過程中, 轉(zhuǎn)基因水稻‘Kf8’的土壤氨氧化細菌和反硝化細菌的群落多樣性指數(shù)及組成, 均與非轉(zhuǎn)基因水稻‘Mh86’沒有顯著差異。相關分析結(jié)果表明土壤氨氧化細菌和反硝化細菌群落組成均與水稻秸稈還田時間存在顯著相關性(=0.002), 反硝化細菌群落組成還與土層深度顯著相關(=0.024)。本研究表明轉(zhuǎn)抗蟲基因水稻秸稈還田對稻田土壤硝化和反硝化關鍵微生物群落不會產(chǎn)生明顯影響。就土壤微生物群落而言, 轉(zhuǎn)抗蟲基因水稻秸稈還田不存在影響土壤氮素轉(zhuǎn)化與N2O排放的可能。

    轉(zhuǎn)基因水稻; 秸稈還田; 氨氧化細菌;基因; 反硝化細菌;基因; 群落結(jié)構

    秸稈還田能增加土壤有機質(zhì), 改良土壤結(jié)構, 促進微生物活力和作物根系的發(fā)育, 在杜絕了秸稈焚燒所造成的大氣污染的同時具有增肥增產(chǎn)作用。轉(zhuǎn)基因作物可能通過外源基因改變作物根系分泌物、根系及植株殘體的組成及其降解過程, 從而在生長過程中及秸稈還田時影響土壤微生物群落組成及其功能[1]。

    硝化和反硝化作用是由微生物驅(qū)動的土壤氮素轉(zhuǎn)化過程, 硝化和反硝化關鍵微生物群落結(jié)構直接關系到土壤氮素轉(zhuǎn)化與利用、氮肥有效性以及溫室氣體N2O排放等環(huán)境問題。由含有氨單加氧酶(amo)的氨氧化細菌驅(qū)動的氨氧化過程是硝化作用的速率限制性步驟[2], 氨單加氧酶基因()是檢測環(huán)境中氨氧化細菌的標記基因[3-4]。反硝化作用是在硝酸鹽還原酶、亞硝酸鹽還原酶、氧化氮還原酶和氧化亞氮還原酶4種酶連續(xù)催化下完成的, 這4種酶的編碼基因分別為、、和。目前亞硝酸鹽還原酶基因()常用于反硝化細菌的分子標記[5-6]。

    轉(zhuǎn)抗蟲基因水稻具有良好的抗蟲和增產(chǎn)效果, 其外源基因及蛋白隨水稻秸稈降解進入土壤后, 若改變了土壤氨氧化細菌或反硝化細菌的群落結(jié)構, 勢必影響稻田土壤的硝化和反硝化過程, 從而關系到稻田氮肥利用及溫室氣體N2O的排放等問題。

    目前, 已有研究表明轉(zhuǎn)Bt基因水稻秸稈的降解對土壤細菌和真菌群落組成沒有顯著影響[7-8]。但有關轉(zhuǎn)基因水稻秸稈降解對土壤中與氮素轉(zhuǎn)化密切相關的硝化反硝化微生物群落的影響尚少見報道。本研究從硝化反硝化微生物功能基因著手, 開展了轉(zhuǎn)基因水稻秸稈還田土壤中硝化反硝化微生物群落的研究。

    我們以轉(zhuǎn)雙價抗蟲基因水稻為對象, 利用分子生物學技術, 以和為標記基因, 對水稻秸稈還田降解過程中的土壤氨氧化細菌和反硝化細菌群落結(jié)構進行研究, 探討轉(zhuǎn)基因水稻是否存在影響稻田土壤氮素轉(zhuǎn)化與N2O排放的可能。

    1 材料與方法

    1.1 供試水稻材料

    供試水稻材料為轉(zhuǎn)基因恢復系‘Kf8’及其對應的非轉(zhuǎn)基因親本恢復系‘Mh86’?!甂f8’是以三系‘Mh86’為受體, 通過農(nóng)桿菌介導法將雙T-抗蟲PCDMARUBAC載體中(蘇云金芽胞桿菌殺蟲晶體蛋白基因)和(修飾的豇豆胰蛋白酶抑制劑基因)兩個抗蟲基因?qū)胨? 獲得的能穩(wěn)定遺傳和表達的轉(zhuǎn)基因恢復系。轉(zhuǎn)基因恢復系‘Kf8’已于2006年獲得農(nóng)業(yè)部轉(zhuǎn)基因安全委員會批準進行生產(chǎn)性試驗。

    1.2 田間試驗設計

    本研究在連續(xù)進行了6年的田間定位試驗基礎上進行。定位試驗在具有防護措施的聯(lián)動網(wǎng)室中進行, 試驗地點位于福州市郊的吳鳳村(26o01′N, 119o20′E)。該試驗基地屬暖濕亞熱帶季風氣候, 年無霜期長達326 d, 年均氣溫19.6 ℃, 平均濕度77%, 年均降水量1 342.5 mm。供試土壤為砂壤, 土壤有機質(zhì)含量為23.32 g×kg-1, pH 6.15, 全氮1.65 g×kg-1, 全磷(P2O5) 0.76 g×kg-1, 全鉀(K2O) 18.28 g×kg-1, 速效氮154.53 mg×kg-1, 速效磷31.34 mg×kg-1, 速效鉀70.52 mg×kg-1。6年定位試驗均為中稻栽培, 5月播種, 6月插秧, 10月收獲。2個水稻材料, 每個材料3次重復, 共6個小區(qū)隨機排列, 每小區(qū)種植水稻350株, 株行距為20 cm×20 cm, 單本插。定位試驗除不使用殺蟲劑外其余水肥措施和田間管理等均按當?shù)爻R?guī)進行。本研究于定位試驗第6年水稻收獲后(2013年10月), 將水稻全部秸稈粉碎原位填埋于各田間小區(qū)內(nèi)自然降解還田, 填埋深度約30 cm。

    1.3 土壤樣品采集

    分別于秸稈還田后50 d、100 d、150 d和200 d采集0~10 cm、10~20 cm和20~30 cm深度土層的土壤。采用五點法取樣, 每個小區(qū)選取5點, 用土鉆分別采集不同深度土層土壤, 五點采集的土樣混合為一個樣品。田間采集的土樣放入4 ℃冰盒中保存, 帶回實驗室后將混合后的土壤去除根系、雜草、土壤動物和石塊等雜質(zhì)后混勻, 裝入50 mL塑料離心管中,-70 ℃冷凍保存用于土壤微生物群落分析。

    1.4 土壤微生物總DNA提取

    稱取0.5 g解凍后的土壤樣品, 采用FastDNAòSPIN Kit For Soil (QBIOgene)試劑盒進行土壤微生物總DNA的提取, DNA樣品-20 ℃冰箱保存待用。

    1.5 熒光定量PCR

    利用氨氧化細菌基因及反硝化細菌的基因的特異引物(表1), 在ABI PRISM7500 Real-Time PCR System擴增儀上進行了絕對熒光定量PCR分析。反應體系為25 μL, 其中2 μL稀釋5倍的DNA模板加反應液23 μL, 反應液包括12.5 μL SYBR Premix Ex TaqTM(2 x), 0.5 μL ROX Reference Dye II(50 x), 前、后引物各1 μL(5 pmol×L-1, 英駿生物工程公司, 上海)和8 μL超純水。標準曲線的建立方法為: 對土壤DNA用目的基因的特異引物進行PCR擴增, 將PCR擴增產(chǎn)物克隆到pMD-18(TaKaRa, 大連)載體中, 然后對插入正確片段的克隆進行質(zhì)粒提取, 用Nanodrop ND-1000 UV-Vis測定質(zhì)粒濃度, 計算質(zhì)粒DNA拷貝數(shù), 按10倍梯度稀釋質(zhì)粒DNA制作標準曲線[9]。本試驗中和基因的標準曲線范圍分別為9.09′108~9.09′104copies×μL-1和3.67′109~3.67′105copies×μL-1。和基因標準曲線的2均達0.99, 擴增效率分別為1.01和0.99。

    表1 PCR引物及反應條件

    *: 分別在amoA1F和nirS4F引物序列前加一個GC夾序列(CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGG)用于DGGE的PCR擴增。*: A GC clamp (CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGG) was attached to the sequence of amoA1F or nirS4F for PCR of DGGE.

    1.6 聚合酶鏈反應-變性梯度凝膠電泳(PCR- DGGE)

    分別利用氨氧化細菌基因和反硝化細菌基因的帶GC夾特異引物進行PCR擴增[6,10](表1)。反應體系50mL, 其中2mL稀釋5倍的DNA模板加反應液48mL, 反應液包括1mL Taq DNA聚合酶(2.5 U×mL-1), 4mL dNTPs(2.5 mmol×L-1), 5mL 10×PCR- buffer(10倍的PCR緩沖液), 前、后引物各4mL(5 pmol×mL-1)和30mL超純水。PCR產(chǎn)物經(jīng)瓊脂糖電泳EB(溴化乙腚)染色檢查后, 用40mL的PCR產(chǎn)物進行變性梯度凝膠電泳, 采用梯度為35%~50%(基因)和40%~60%(基因)的8%聚丙烯酰胺凝膠[化學變性劑為100%尿素7 mol×L-1和40%(/)的去離子甲酰胺]在1倍TAE緩沖液中150 V 60 ℃下電泳5 h。電泳后用10 mL SyBR green I(1倍TAE稀釋10 000倍)核酸染料染色45 min, 然后用FluoChem SP(Alpha Innotech)成像系統(tǒng)拍照。

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

    采用Quantity One (Bio-Rad)軟件對DGGE圖譜中DNA條帶的位置和亮度進行數(shù)字化, 計算微生物群落多樣性指數(shù)Shannon-Weiner index (′)=-?()× ln(/) (為每個帶的亮度,為某一樣品所有帶亮度的和); 利用Canoco 4.5生態(tài)學軟件進行微生物群落結(jié)構的主成分分析(PCA); 采用SPSS 13.0軟件進行數(shù)據(jù)方差分析。

    2 結(jié)果與分析

    2.1 土壤硝化反硝化微生物豐度變化

    利用熒光定量PCR技術分析了土壤硝化作用關鍵微生物氨氧化細菌基因與反硝化細菌基因豐度, 結(jié)果如圖1所示。從圖1A可見, 在水稻秸稈還田過程中稻田0~10 cm土壤中氨氧化細菌豐度較高, 4個取樣時間內(nèi)的基因數(shù)量都可以達106copies?g-1(干土)以上; 而10~20 cm及20~30 cm土層中的基因數(shù)量只到達104~105copies?g-1(干土); 前者顯著(<0.05)高于后者。在水稻秸稈還田的各個取樣時間內(nèi), 轉(zhuǎn)基因水稻‘Kf8’的氨氧化細菌基因豐度與其對應的常規(guī)水稻‘Mh86’沒有顯著差異, 0~10 cm、10~20 cm和20~30 cm 3個深度土層中‘Kf8’秸稈還田土壤的基因數(shù)量與‘Mh86’均處于同一數(shù)量級上, 且差異不顯著。說明‘Kf8’與‘Mh86’兩個品種水稻秸稈降解過程中土壤氨氧化細菌基因豐度變化是一致的。

    水稻秸稈還田過程中土壤反硝化細菌基因豐度變化如圖1B所示, 結(jié)果表明秸稈還田的各個取樣時間內(nèi), 轉(zhuǎn)基因水稻‘Kf8’與非轉(zhuǎn)基因水稻‘Mh86’在各個深度土層中的反硝化細菌基因豐度都沒有顯著差異。隨秸稈還田時間的延長各深度土層的基因豐度存在逐漸增加的趨勢, ‘Kf8’處理的0~10 cm、10~20 cm和20~30 cm土層中基因數(shù)量在秸稈還田200 d時分別是50 d時的1.58倍、1.91倍和1.41倍, ‘Mh86’處理分別是1.70倍、1.77倍和1.26倍。說明隨著‘Kf8’與‘Mh86’兩品種水稻秸稈的降解土壤反硝化細菌基因豐度都在增高。

    圖1 轉(zhuǎn)基因水稻‘Kf8’和常規(guī)水稻‘Mh86’秸稈還田后不同時間不同深度土壤氨氧化細菌amoA基因(A)和反硝化細菌nirS基因(B)豐度的變化

    同一深度土層不同字母表示差異顯著(0.05)。Values of the same depth of soil in the same time followed by different letters are significantly different (< 0.05).

    2.2 稻田土壤硝化反硝化微生物群落多樣性變化

    通過對稻田土壤中基因和基因的PCR擴增及其產(chǎn)物的變性梯度凝膠電泳(DGGE), 分析了稻田土壤中硝化和反硝化細菌的群落多樣性。圖2A顯示在水稻秸稈還田降解的200 d內(nèi)土壤氨氧化細菌群落多樣性指數(shù)的變化。在各取樣時間內(nèi)及各深度土層中氨氧化細菌群落多樣性指數(shù)在‘Kf8’與‘Mh86’間均沒有顯著差異; 同一品種水稻的氨氧化細菌群落多樣性指數(shù)在不同取樣時間及不同深度土層中也都沒有明顯差異; 秸稈還田降解的200 d內(nèi)2個品種水稻的土壤氨氧化細菌群落多樣性指數(shù)基本維持在2.0左右。

    在水稻秸稈還田降解的200 d內(nèi), 各深度土層的反硝化細菌群落多樣性指數(shù)在‘Kf8’和‘Mh86’間也均沒有差異(圖2B); 同一秸稈還田時間內(nèi), 各個品種水稻的土壤反硝化細菌群落多樣性指數(shù)在不同深度土層間同樣沒有明顯差異(圖2B)。秸稈還田降解的200 d內(nèi)2個品種水稻的土壤反硝化細菌群落多樣性指數(shù)基本維持在2.5左右。

    圖2 轉(zhuǎn)基因水稻‘Kf8’和常規(guī)水稻‘Mh86’秸稈還田后不同時間不同深度土壤氨氧化細菌(A)和反硝化細菌(B)多樣性指數(shù)的變化

    同一深度土層不同字母表示差異顯著(0.05)。Values of the same depth of soil in the same time followed by different letters are significantly different (< 0.05).

    2.3 稻田土壤硝化反硝化微生物群落組成的變化

    基于DGGE圖譜利用Canoco 4.5生態(tài)學軟件進行了秸稈還田過程中土壤硝化和反硝化微生物群落組成的主成分分析(PCA), 并將微生物群落組成與水稻品種、取樣時間和土層深度3個因素進行了相關性分析。

    圖3a顯示, 2個品種水稻各深度土層4個取樣時間的氨氧化細菌群落組成在變異較大的主成分1(PC1)方向上存在分布差異, 其中秸稈還田100 d和150 d的分布較為接近, 并分別與50 d、200 d的分布差異較大。同一取樣時間內(nèi)水稻‘Kf8’與水稻‘Mh86’間以及不同深度土層間的氨氧化細菌群落組成的分布差異均表現(xiàn)在變異較小的主成分2(PC2)方向上。同時, 相關性分析結(jié)果表明土壤氨氧化細菌的群落組成僅與取樣時間即秸稈還田時間長短存在顯著相關(=0.002), 與水稻品種和土層深度都沒有顯著相關性。

    土壤反硝化細菌群落組成的變化趨勢與土壤氨氧化細菌一致(圖3b), 50 d、100 d、150 d和 200 d的不同水稻品種及不同深度土層的反硝化細菌群落組成均在PC1方向上分布差異明顯, 同一取樣時間內(nèi)同一品種水稻的不同深度土層的反硝化細菌群落組成在PC2方向上也存在較明顯的分布差異。相關性分析結(jié)果表明, 土壤反硝化細菌的群落組成與取樣時間(=0.002)和土層深度(=0.024)均存在顯著相關, 但與水稻品種沒有明顯相關性。

    圖3 轉(zhuǎn)基因水稻‘Kf8’和常規(guī)水稻‘Mh86’秸稈還田后不同時間不同深度土壤氨氧化細菌(a)和反硝化細菌(b)群落組成的主成分分析

    3 討論與結(jié)論

    目前, 關于轉(zhuǎn)基因作物對環(huán)境微生物影響的研究結(jié)果依據(jù)作物種類、目標基因、微生物種類等因素而存在差異, 且研究多集中于對土壤總微生物如細菌或真菌的群落方面。如: 在轉(zhuǎn)抗病基因木瓜的研究中發(fā)現(xiàn), 轉(zhuǎn)基因品種的種植能夠改變土壤真菌或細菌的群落組成[11]; 而土壤微生物群落結(jié)構在轉(zhuǎn)Bt基因抗蟲玉米和非轉(zhuǎn)基因玉米間沒有差異[12]。本研究所在定位試驗的前期研究證實了轉(zhuǎn)雙價抗蟲基因水稻定位種植4年后對土壤細菌和真菌的群落組成與豐度都沒有顯著影響[13]。還有研究表明轉(zhuǎn)Bt基因抗蟲水稻的秸稈或根系還田降解過程中對土壤細菌或真菌的群落組成沒有影響[7-8]。

    近年來對轉(zhuǎn)基因作物土壤中功能微生物群落變化的研究逐漸增多, 如對驅(qū)動土壤氮素轉(zhuǎn)化關鍵微生物群落結(jié)構的研究。本研究定位試驗前期工作已證實在一定時期內(nèi)種植轉(zhuǎn)雙價抗蟲基因水稻不會影響土壤氨氧化細菌的群落組成與豐度[14]。本研究進一步分析了轉(zhuǎn)雙價抗蟲基因水稻秸稈還田降解過程中對土壤硝化和反硝化微生物群落組成與豐度的影響。結(jié)果表明轉(zhuǎn)雙價抗蟲基因水稻秸稈還田降解對土壤氨氧化細菌和反硝化細菌群落豐度與結(jié)構都沒有明顯影響。土壤氨氧化細菌和反硝化細菌的群落結(jié)構僅與秸稈還田降解的時間長短或土層深度相關, 而與轉(zhuǎn)基因與非轉(zhuǎn)基因的水稻品種沒有相關性。另外, 由于氨氧化是個好氧的微生物轉(zhuǎn)化過程, 所以秸稈還田降解過程中通氣性較好的表層土壤中含有更豐富的氨氧化細菌。

    轉(zhuǎn)基因技術研發(fā)的目的是將一些有利于作物生長發(fā)育的外源基因?qū)肽繕俗魑? 使其具有抗蟲、抗病、抗逆、改善品質(zhì)及增產(chǎn)的優(yōu)勢。但外源基因的轉(zhuǎn)入有可能引起作物根系分泌物、根系及植株殘體組成的改變, 外源基因及其表達蛋白也可以通過根系分泌物或植株殘體降解直接進入到土壤中, 這些過程都有可能使種植轉(zhuǎn)基因作物的土壤的微生物群落及其生態(tài)功能發(fā)生變化, 從而影響到土壤環(huán)境的生態(tài)可持續(xù)性。

    轉(zhuǎn)Bt基因作物中的殺蟲晶體蛋白主要通過根系分泌物、花粉飄落和秸稈還田釋放到環(huán)境中[15]。進入土壤的Bt蛋白能夠迅速被土壤活性顆粒吸附, 1~3 h就能達到吸附平衡[16]。吸附態(tài)Bt蛋白的殺蟲活性可保持數(shù)周或數(shù)月, 轉(zhuǎn)Bt基因玉米的盆栽和大田試驗證明 Bt毒素的殺蟲活性能保持180 d以上[17]; 轉(zhuǎn)Bt基因玉米、棉花和馬鈴薯秸稈分解釋放的Bt毒素殺蟲活性都可保持數(shù)周以上[18-20]。Bt毒蛋白在土壤環(huán)境中失活或消除的主要途徑為進入昆蟲體內(nèi)后被降解、被陽光中紫外線照射分解及微生物的降解和最終礦化作用, 其中土壤微生物是影響B(tài)t毒蛋白在土壤中存留累積的最關鍵因子[21]。

    本研究結(jié)果初步揭示轉(zhuǎn)雙價抗蟲基因水稻的秸稈還田降解對土壤氮素轉(zhuǎn)化關鍵微生物群落豐度與組成沒有明顯影響, 不存在影響土壤氮素轉(zhuǎn)化與N2O排放的可能。但由于土壤N2O的產(chǎn)生有多種途徑, 由多個酶催化進行, 且不同途徑的不同酶催化反應間存在相互制約或促進的復雜關系[22-24]。因此, 還需要通過對參與土壤N2O產(chǎn)生的關鍵微生物功能基因多樣性及其作用進行系統(tǒng)的分析研究, 才能深入揭示轉(zhuǎn)基因水稻對稻田溫室氣體N2O排放的影響及機制。

    [1] MASOERO F, MOSCHINI M, ROSSI F, et al. Nutritive value, mycotoxin contamination and in vitro rumen fermentation of normal and genetically modified corn () grown in northern Italy[J]. Maydica, 1999, 44(3): 205–209

    [2] CAPONE D G. Methane, nitrogen oxides and halomethanes[M]// Rodgers J E, Whitman W B. Microbial Consumption of Greenhouse Gases. Washington, DC: American Society for Microbiology, 1991: 225–275

    [3] STEPHEN J R, McCAIG A E, SMITH Z, et al. Molecular diversity of soil and marine 16S rRNA gene sequences related to β-subgroup ammonia-oxidizing bacteria[J]. Applied and Environmental Microbiology, 1996, 62(11): 4147–4154

    [4] WOESE C R, WEISBURG W G, HAHN C M, et al. The phylogeny of purple bacteria: The gamma subdivision[J]. Systematic and Applied Microbiology, 1985, 6(1): 25–33

    [5] ZUMFT W G. Cell biology and molecular basis of denitrification[J]. Microbiology and Molecular Biology Reviews, 1997, 61(4): 533–616

    [6] 宋亞娜, 林智敏, 林艷. 氮肥對稻田土壤反硝化細菌群落結(jié)構和豐度的影響[J]. 中國生態(tài)農(nóng)業(yè)學報, 2012, 20(1): 7–12 SONG Y N, LIN Z M, LIN Y. Response of denitrifying bacteria community structure and abundance to nitrogen in paddy fields[J]. Chinese Journal of Eco-Agriculture, 2012, 20(1): 7–12

    [7] LU H H, WU W X, CHEN Y X, et al. Decomposition oftransgenic rice residues and response of soil microbial community in rapeseed-rice cropping system[J]. Plant and Soil, 2010, 336(1/2): 279–290

    [8] WU W X, YE Q F, MIN H, et al. Bt-transgenic rice straw affects the culturable microbiota and dehydrogenase and phosphatase activities in a flooded paddy soil[J]. Soil Biology and Biochemistry, 2004, 36(2): 289–295

    [9] 宋亞娜, 陳在杰, 林智敏. 水稻生育期內(nèi)紅壤稻田氨氧化微生物數(shù)量和硝化勢的變化[J]. 中國生態(tài)農(nóng)業(yè)學報, 2010, 18(5): 954–958 SONG Y N, CHEN Z J, LIN Z M. Abundance of ammonia-oxidizer and potential nitrification rate of quaternary red-clay paddy soil during rice growth[J]. Chinese Journal of Eco-Agriculture, 2010, 18(5): 954–958

    [10] 袁飛, 冉煒, 胡江, 等. 變性梯度凝膠電泳法研究我國不同土壤氨氧化細菌群落組成及活性[J]. 生態(tài)學報, 2005, 25(6): 1318–1324 YUAN F, RAN W, HU J, et al. Ammonia-oxidizing bacteria communities and their influence on the nitrification potential of Chinese soils measured by denaturing gradient gel electrophoresis (DGGE)[J]. Acta Ecologica Sinica, 2005, 25(6): 1318–1324

    [11] WEI X D, ZOU H L, CHU L M, et al. Field released transgenic papaya effect on soil microbial communities and enzyme activities[J]. Journal of Environmental Sciences, 2006, 18(4): 734–740

    [12] GRIFFITHS B S, CAUL S, THOMPSON J, et al. A comparison of soil microbial community structure, protozoa and nematodes in field plots of conventional and genetically modified maize expressing theCryIAb toxin[J]. Plant and Soil, 2005, 275(1/2): 135–146

    [13] SONG Y N, SU J, CHEN R, et al. Diversity of microbial community in a paddy soil withtransgenic rice[J]. Pedosphere, 2014, 24(3): 349–358

    [14] 宋亞娜, 蘇軍, 林艷, 等. 轉(zhuǎn)1Ac基因水稻對土壤氨氧化細菌群落組成和豐度的影響[J]. 生物安全學報, 2012, 21(1): 67–73 SONG Y N, SU J, LIN Y, et al. Effect of1Actransgenic rice on community composition and abundance of ammonia-oxidizing bacteria in paddy soil[J]. Journal of Biosafety, 2012, 21(1): 67–73

    [15] 王建武, 駱世明, 馮遠嬌, 等. 轉(zhuǎn)Bt基因作物Bt毒素在土壤中的環(huán)境去向及其生態(tài)效應[J]. 生態(tài)學報, 2003, 23(4): 797–804 WANG J W, LUO S M, FENG Y J, et al. Environmental fate and ecological effects of Bt toxin from transgenic Bt crops in soil[J]. Acta Ecologica Sinica, 2003, 23(4): 797–804

    [16] ZHOU X Y, HUANG Q Y, CHEN S W, et al. Adsorption of the insecticidal protein ofon montmorillonite, kaolinite, silica, goethite and red soil[J]. Applied Clay Science, 2005, 30(2): 87–93

    [17] SAXENA D, FLORES S, STOTZKY G. Bt toxin is released in root exudates from 12 transgenic corn hybrids representing three transformation events[J]. Soil Biology and Biochemistry, 2002, 34(1): 133–137

    [18] SIMS S R, HOLDEN L R. Insect bioassay for determining soil degradation ofsubspCryIA (b) protein in corn tissue[J]. Environmental Entomology, 1996, 25(3): 659–664

    [19] DONEGAN K, PALM C J, FIELAND V J, et al. Changes in levels, species and DNA fingerprints of soil microorganisms associated with cotton expressing thevarendotoxin[J]. Applied Soil Ecology, 1995, 2(2): 111–124

    [20] DONEGAN K K, SCHALLER D L, STONE J K, et al. Microbial populations, fungal species diversity and plant pathogen levels in field plots of potato plants expressing thevarendotoxin[J]. Transgenic Research, 1996, 5(1): 25–35

    [21] 白耀宇, 蔣明星, 程家安, 等. 轉(zhuǎn)Bt基因作物Bt毒蛋白在土壤中的安全性研究[J]. 應用生態(tài)學報, 2003, 14(11): 2062–2066 BAI Y Y, JIANG M X, CHENG J A, et al. Advances in safetystudies of soil Bt toxin proteins released from transgenic Bt crops[J]. Chinese Journal of Applied Ecology, 2003, 14(11): 2062–2066

    [22] FIRESTONE M K, DAVIDSON E A. Microbial basis of NO and N2O production and consumption in soils[M]//Andreae M O, Schimel D S. Exchange of Trace Gases Between Terrestrial Ecosystems and the Atmosphere. New York: John Wiley, 1989: 7–21

    [23] ZHU X, BURGER M, DOANE T A, et al. Ammonia oxidation pathways and nitrifier denitrification are significant sources of N2O and NO under low oxygen availability[J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(16): 6328–6333

    [24] 朱永官, 王曉輝, 楊小茹, 等. 農(nóng)田土壤N2O產(chǎn)生的關鍵微生物過程及減排措施[J]. 環(huán)境科學, 2014, 35(2): 792–800ZHU Y G, WANG X H, YANG X R, et al. Key microbial processes in nitrous oxide emissions of agricultural soil and mitigation strategies[J]. Environmental Science, 2014, 35(2): 792–800

    Effect of transgenic rice straw return to soil on nitrification and denitrification microbial community*

    WANG Peixuan1,3, XU Yan2, SONG Yana1**

    (1. Institute of Biological Technology, Fujian Academy of Agricultural Sciences, Fuzhou 350003, China; 2. Fuzhou Shidai Middle School, Fuzhou 350007, China; 3. Affiliated High School of Fujian Normal University, Fuzhou 350007, China)

    The cultivation of genetically modified plants may have the potential to alter microbial community structure and function in soils through root exudates and plant residues. Ammonia-oxidizing bacteria or denitrifying bacteria are the key microbes for nitrification or denitrification. A change in community structure of ammonia-oxidizing bacteria or denitrifying bacteria can affect the conversion and utilization of nitrogen in soil. The purpose of this study was to explore the possibility of transgenic rice to induce change in nitrogen transformation and N2O emission in paddy soils. In the study, the abundance and composition of ammonia-oxidizing bacteria or denitrifying bacteria in paddy soils under straw return to soil undertransgenic gene rice ‘Kf8’ or non-transgenic rice ‘Mh86’ were analyzed by real-time PCR and denaturing gradient gel electrophoresis (DGGE) based ongene orgene. The results showed that there were no differences in the abundance ofgene orgene in the soil with returned transgenic rice ‘Kf8’ straw and non-transgenic rice ‘Mh86’ straw. The abundance ofgene in the 0–10 cm soil layer was significantly (< 0.05) higher than that in the 10–20 cm and 20–30 cm soil layers under degraded transgenic rice ‘Kf8’ straw or non-transgenic rice ‘Mh86’ straw. The abundance ofgene in the soil increased with the time of returned straw to soil for either transgenic rice ‘Kf8’ or no-transgenic rice ‘Mh86’. At the same time, the Shannon-Weiner index and composition of ammonia-oxidizing bacteria or denitrifying bacteria in the soil under degraded transgenic rice ‘Kf8’ straw were similar to those under degraded non-transgenic rice ‘Mh86’ straw. The composition of ammonia-oxidizing bacteria in the soil was significantly correlated with the time of rice straw return (= 0.002). Also the abundance of denitrifying bacteria was significantly correlated with the time of rice straw return (= 0.002) and depth of soil (= 0.024). The findings demonstrated that there was no significant effect of returnedtransgenic rice straw to soil on key microbial communities for nitrification or denitrification in soil. In terms of soil microbial community, there was no signiifanct effect of returnedtransgenic rice straw to soil on nitrogen transformation and N2O emission in paddy soils.

    Transgenic rice; Straw return; Ammonia-oxidizing bacteria;gene; Denitrifying bacteria;gene; Community structure

    Corresponding author, E-mail: syana@sina.com

    Feb. 13, 2017; accepted May 18, 2017

    10.13930/j.cnki.cjea.170794

    S154.3

    A

    1671-3990(2018)01-0008-08

    通信作者:宋亞娜, 主要從事土壤微生物分子生態(tài)學研究。E-mail: syana@sina.com 王沛譞, 主要從事土壤微生物分子生態(tài)學研究。E-mail: 1092546064@qq.com

    2017-02-13

    2017-05-18

    *This study was supported by the Natural Science Foundation of Fujian Province (2015J01105), the Public Welfare Foundation of Fujian Province (2015R1019-12) and the 30thYouth Innovation Competition Project in Fujian Province.

    *福建省自然科學基金項目(2015J01105)、福建省公益類科研院所專項(2015R1019-12)和第30屆福建省青少年創(chuàng)新大賽項目資助

    王沛譞, 徐焱, 宋亞娜. 轉(zhuǎn)基因水稻秸稈還田對土壤硝化反硝化微生物群落的影響[J]. 中國生態(tài)農(nóng)業(yè)學報, 2018, 26(1): 8-15

    Wang P X, Xu Y, Song Y N. Effect of transgenic rice straw return to soil on nitrification and denitrification microbial community[J]. Chinese Journal of Eco-Agriculture, 2018, 26(1): 8-15

    猜你喜歡
    水稻深度
    什么是海水稻
    有了這種合成酶 水稻可以耐鹽了
    水稻種植60天就能收獲啦
    軍事文摘(2021年22期)2021-11-26 00:43:51
    油菜可以像水稻一樣實現(xiàn)機插
    深度理解一元一次方程
    一季水稻
    文苑(2020年6期)2020-06-22 08:41:52
    水稻花
    文苑(2019年22期)2019-12-07 05:29:00
    深度觀察
    深度觀察
    深度觀察
    国产一区二区在线观看日韩| 久久人人精品亚洲av| 国产精品福利在线免费观看| 国产av不卡久久| 精品免费久久久久久久清纯| 色视频www国产| 亚洲人成网站高清观看| 国产精品乱码一区二三区的特点| 久久国产乱子免费精品| 免费大片18禁| 久久99热6这里只有精品| 欧美一级a爱片免费观看看| 大香蕉久久网| 草草在线视频免费看| 欧美性猛交╳xxx乱大交人| 久久久色成人| 久久久精品大字幕| 免费观看在线日韩| 精品一区二区免费观看| 精品久久久久久久末码| 欧美另类亚洲清纯唯美| 黄色配什么色好看| 3wmmmm亚洲av在线观看| 女人十人毛片免费观看3o分钟| av女优亚洲男人天堂| 欧美国产日韩亚洲一区| 看免费成人av毛片| 五月伊人婷婷丁香| 国产精品久久久久久av不卡| 亚洲图色成人| 午夜精品一区二区三区免费看| 成人精品一区二区免费| 国产午夜精品论理片| 99久久精品一区二区三区| 99热精品在线国产| 成人一区二区视频在线观看| 精品久久久噜噜| 精品人妻一区二区三区麻豆 | 99九九线精品视频在线观看视频| 久久精品国产清高在天天线| 美女免费视频网站| 无遮挡黄片免费观看| 久久久成人免费电影| 国产精品一及| 亚洲自偷自拍三级| 亚洲av.av天堂| 香蕉av资源在线| 亚洲内射少妇av| 精品久久国产蜜桃| 听说在线观看完整版免费高清| 麻豆国产97在线/欧美| 亚洲无线在线观看| 国产精品精品国产色婷婷| 国产精品精品国产色婷婷| 亚洲av熟女| 俺也久久电影网| 国产真实乱freesex| 男插女下体视频免费在线播放| 最后的刺客免费高清国语| 亚洲在线观看片| 91精品国产九色| 女人被狂操c到高潮| 午夜精品国产一区二区电影 | 小蜜桃在线观看免费完整版高清| 国产激情偷乱视频一区二区| 老熟妇乱子伦视频在线观看| 免费av不卡在线播放| 亚洲精品亚洲一区二区| 国产一区二区三区av在线 | 色吧在线观看| 激情 狠狠 欧美| 少妇裸体淫交视频免费看高清| 18禁在线播放成人免费| 午夜老司机福利剧场| 亚洲人与动物交配视频| 狂野欧美白嫩少妇大欣赏| 干丝袜人妻中文字幕| 亚洲精品乱码久久久v下载方式| 中文字幕人妻熟人妻熟丝袜美| 特大巨黑吊av在线直播| 级片在线观看| 99热精品在线国产| 伊人久久精品亚洲午夜| 免费看光身美女| 免费一级毛片在线播放高清视频| 人妻丰满熟妇av一区二区三区| 91在线观看av| 91在线观看av| 三级经典国产精品| 亚洲激情五月婷婷啪啪| 日韩中字成人| 日本与韩国留学比较| 又粗又爽又猛毛片免费看| 久久精品久久久久久噜噜老黄 | 国产伦一二天堂av在线观看| 变态另类成人亚洲欧美熟女| 日韩精品青青久久久久久| 午夜a级毛片| 久久久精品欧美日韩精品| 嫩草影院新地址| 亚洲av熟女| 国产精品嫩草影院av在线观看| 国产真实乱freesex| 18禁裸乳无遮挡免费网站照片| 成人漫画全彩无遮挡| 国国产精品蜜臀av免费| 亚洲aⅴ乱码一区二区在线播放| 午夜精品在线福利| 精品人妻一区二区三区麻豆 | 婷婷精品国产亚洲av| 三级国产精品欧美在线观看| 精华霜和精华液先用哪个| 久久精品影院6| 噜噜噜噜噜久久久久久91| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美中文字幕日韩二区| 亚洲精品456在线播放app| 成人一区二区视频在线观看| 国产免费男女视频| 小说图片视频综合网站| 12—13女人毛片做爰片一| 搡老岳熟女国产| 波多野结衣高清作品| 狂野欧美激情性xxxx在线观看| 两个人视频免费观看高清| 亚洲国产色片| 成年女人毛片免费观看观看9| 乱码一卡2卡4卡精品| 欧美潮喷喷水| 丝袜喷水一区| 麻豆一二三区av精品| 可以在线观看毛片的网站| h日本视频在线播放| 午夜福利18| 亚洲精品色激情综合| 黄色视频,在线免费观看| 春色校园在线视频观看| 好男人在线观看高清免费视频| 亚洲欧美中文字幕日韩二区| 午夜精品在线福利| 别揉我奶头~嗯~啊~动态视频| 麻豆乱淫一区二区| 国产精品伦人一区二区| 国产av麻豆久久久久久久| 天天一区二区日本电影三级| 2021天堂中文幕一二区在线观| 亚洲熟妇中文字幕五十中出| 国模一区二区三区四区视频| 午夜福利成人在线免费观看| 搡老熟女国产l中国老女人| 国产视频一区二区在线看| 美女免费视频网站| 国产精品三级大全| 给我免费播放毛片高清在线观看| 日本三级黄在线观看| 一个人看视频在线观看www免费| 悠悠久久av| 国产成人91sexporn| 久久韩国三级中文字幕| 又爽又黄无遮挡网站| 看黄色毛片网站| 国产精品一区二区性色av| 2021天堂中文幕一二区在线观| 免费高清视频大片| 最近在线观看免费完整版| 国产成人a区在线观看| 成人亚洲欧美一区二区av| 大型黄色视频在线免费观看| 九色成人免费人妻av| 亚洲国产精品成人久久小说 | videossex国产| 简卡轻食公司| 欧美极品一区二区三区四区| 精品久久久久久久久久免费视频| 黄色一级大片看看| 日韩精品青青久久久久久| 色噜噜av男人的天堂激情| 99久久中文字幕三级久久日本| 欧美日韩乱码在线| 国产成人aa在线观看| 十八禁国产超污无遮挡网站| 在线国产一区二区在线| 白带黄色成豆腐渣| 成人毛片a级毛片在线播放| 在线观看一区二区三区| 一夜夜www| 国产日本99.免费观看| 日本一本二区三区精品| 欧美最新免费一区二区三区| 无遮挡黄片免费观看| 国产一区二区在线av高清观看| 性插视频无遮挡在线免费观看| 麻豆国产av国片精品| 人妻久久中文字幕网| 人妻丰满熟妇av一区二区三区| 精品不卡国产一区二区三区| 国产高潮美女av| 亚洲美女视频黄频| 舔av片在线| 最近在线观看免费完整版| 久久久色成人| 在线播放国产精品三级| 亚洲国产精品合色在线| 91午夜精品亚洲一区二区三区| 国产一区亚洲一区在线观看| 在线观看午夜福利视频| 欧美成人免费av一区二区三区| 人妻丰满熟妇av一区二区三区| 精品不卡国产一区二区三区| 精品久久久久久久久亚洲| 有码 亚洲区| 亚洲欧美日韩高清专用| 欧美成人精品欧美一级黄| 亚洲国产日韩欧美精品在线观看| 国产探花极品一区二区| 亚洲精品亚洲一区二区| 亚洲最大成人av| 少妇人妻精品综合一区二区 | 最近手机中文字幕大全| 亚洲av电影不卡..在线观看| 国产亚洲精品av在线| 国产乱人视频| 九九久久精品国产亚洲av麻豆| 亚洲婷婷狠狠爱综合网| 亚洲国产精品久久男人天堂| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 我的女老师完整版在线观看| 国产成人a∨麻豆精品| 成年女人毛片免费观看观看9| 亚洲欧美成人精品一区二区| 成人av一区二区三区在线看| av免费在线看不卡| 免费观看人在逋| 性欧美人与动物交配| 亚洲欧美日韩东京热| 亚洲av中文av极速乱| 亚洲熟妇中文字幕五十中出| 别揉我奶头~嗯~啊~动态视频| 中出人妻视频一区二区| 特级一级黄色大片| 久久国产乱子免费精品| 美女xxoo啪啪120秒动态图| 亚洲精品国产av成人精品 | 男女啪啪激烈高潮av片| 九九久久精品国产亚洲av麻豆| 日韩成人av中文字幕在线观看 | .国产精品久久| 久久久久久伊人网av| 午夜福利在线观看吧| 成年av动漫网址| 久久久久久久久久成人| 国内少妇人妻偷人精品xxx网站| 变态另类丝袜制服| 97在线视频观看| 99热6这里只有精品| 日日摸夜夜添夜夜添av毛片| 亚洲真实伦在线观看| 久久久精品94久久精品| 麻豆成人午夜福利视频| 淫妇啪啪啪对白视频| 午夜福利在线观看免费完整高清在 | 欧美xxxx黑人xx丫x性爽| 干丝袜人妻中文字幕| 欧美绝顶高潮抽搐喷水| 国产乱人视频| 日韩欧美免费精品| 女人被狂操c到高潮| 国产成人一区二区在线| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| av在线播放精品| 中文字幕人妻熟人妻熟丝袜美| 国内少妇人妻偷人精品xxx网站| 免费看美女性在线毛片视频| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 一级毛片电影观看 | 午夜激情福利司机影院| 91在线观看av| 全区人妻精品视频| 麻豆国产97在线/欧美| 色视频www国产| 亚洲av五月六月丁香网| 黑人高潮一二区| 毛片女人毛片| 日韩一区二区视频免费看| 免费无遮挡裸体视频| 香蕉av资源在线| 在线免费十八禁| 青春草视频在线免费观看| 蜜臀久久99精品久久宅男| 欧美zozozo另类| 欧美国产日韩亚洲一区| 日本-黄色视频高清免费观看| 99在线人妻在线中文字幕| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| a级一级毛片免费在线观看| 国产乱人视频| 91久久精品国产一区二区三区| 在线免费十八禁| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区 | 天天躁夜夜躁狠狠久久av| 99在线人妻在线中文字幕| 日本欧美国产在线视频| 欧美人与善性xxx| 91麻豆精品激情在线观看国产| 日韩 亚洲 欧美在线| 精品久久久久久久人妻蜜臀av| 国产真实伦视频高清在线观看| 六月丁香七月| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 联通29元200g的流量卡| 婷婷精品国产亚洲av在线| 一级毛片久久久久久久久女| 日韩欧美国产在线观看| 国产白丝娇喘喷水9色精品| 久久亚洲精品不卡| 亚洲va在线va天堂va国产| 联通29元200g的流量卡| 久久精品91蜜桃| 黄色视频,在线免费观看| 中文字幕免费在线视频6| 小说图片视频综合网站| 日韩三级伦理在线观看| 国产69精品久久久久777片| 亚洲欧美日韩东京热| 婷婷六月久久综合丁香| 丝袜美腿在线中文| 成人亚洲精品av一区二区| 亚洲一区二区三区色噜噜| 我要看日韩黄色一级片| 国产视频内射| 白带黄色成豆腐渣| 亚洲无线在线观看| 久久精品国产亚洲av天美| 午夜福利18| 日韩 亚洲 欧美在线| 亚洲电影在线观看av| 丝袜美腿在线中文| av在线亚洲专区| 亚洲熟妇熟女久久| 久久久久久久久久久丰满| av在线观看视频网站免费| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av| 俄罗斯特黄特色一大片| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| 免费观看的影片在线观看| 一个人看视频在线观看www免费| 一级毛片电影观看 | 亚洲自拍偷在线| 一级毛片我不卡| 一级毛片aaaaaa免费看小| 大型黄色视频在线免费观看| 精品一区二区三区人妻视频| 久久婷婷人人爽人人干人人爱| а√天堂www在线а√下载| aaaaa片日本免费| 国产伦精品一区二区三区四那| 麻豆久久精品国产亚洲av| 欧美性感艳星| 秋霞在线观看毛片| 成人一区二区视频在线观看| 久久精品人妻少妇| 亚洲欧美成人综合另类久久久 | 性插视频无遮挡在线免费观看| 亚洲不卡免费看| 精品福利观看| 美女大奶头视频| 久久久精品94久久精品| 99久国产av精品| 久久久久久久久久成人| 日韩av不卡免费在线播放| 国产一区亚洲一区在线观看| 色视频www国产| 欧美最黄视频在线播放免费| 亚洲人成网站在线播放欧美日韩| 久久精品国产清高在天天线| 欧美zozozo另类| 欧美色视频一区免费| 国产久久久一区二区三区| 日韩在线高清观看一区二区三区| 麻豆久久精品国产亚洲av| 午夜福利18| 岛国在线免费视频观看| 亚洲不卡免费看| 亚洲专区国产一区二区| 床上黄色一级片| 一级黄片播放器| 天堂网av新在线| 亚洲中文字幕一区二区三区有码在线看| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 99精品在免费线老司机午夜| 亚洲精品一区av在线观看| 欧美日韩乱码在线| 成人美女网站在线观看视频| 亚洲精品一区av在线观看| 最近中文字幕高清免费大全6| 国产久久久一区二区三区| 91久久精品国产一区二区成人| av黄色大香蕉| 亚洲激情五月婷婷啪啪| 亚洲高清免费不卡视频| 青春草视频在线免费观看| 亚洲激情五月婷婷啪啪| 两个人视频免费观看高清| 国产精品一区二区免费欧美| 亚洲经典国产精华液单| 久久久成人免费电影| 久久午夜亚洲精品久久| 尤物成人国产欧美一区二区三区| 三级国产精品欧美在线观看| 精品乱码久久久久久99久播| 在现免费观看毛片| 免费av观看视频| 寂寞人妻少妇视频99o| 国产高清不卡午夜福利| 色噜噜av男人的天堂激情| 日本a在线网址| 女的被弄到高潮叫床怎么办| 中文亚洲av片在线观看爽| 黑人高潮一二区| 男人的好看免费观看在线视频| 日韩成人伦理影院| 中文字幕av成人在线电影| 亚洲av美国av| 久久99热6这里只有精品| 日本免费一区二区三区高清不卡| avwww免费| 美女cb高潮喷水在线观看| av黄色大香蕉| 在线观看免费视频日本深夜| 精品一区二区免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品久久国产高清桃花| 国产又黄又爽又无遮挡在线| 97碰自拍视频| 亚洲精品久久国产高清桃花| 久久久久免费精品人妻一区二区| 精品乱码久久久久久99久播| 中文字幕免费在线视频6| 神马国产精品三级电影在线观看| 日韩高清综合在线| 国产成人一区二区在线| 亚洲av免费在线观看| 国产精品久久电影中文字幕| 午夜福利高清视频| 嫩草影院入口| 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| 我要看日韩黄色一级片| 高清毛片免费看| 高清午夜精品一区二区三区 | 成人国产麻豆网| 国产中年淑女户外野战色| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 此物有八面人人有两片| 哪里可以看免费的av片| 久久久国产成人免费| 免费观看精品视频网站| 日韩制服骚丝袜av| 国产精品嫩草影院av在线观看| 国产色婷婷99| 欧美bdsm另类| 国产精华一区二区三区| 亚洲欧美成人精品一区二区| 99热只有精品国产| 国产视频内射| 亚洲一区二区三区色噜噜| 99riav亚洲国产免费| 日韩欧美在线乱码| 国产高清激情床上av| 可以在线观看的亚洲视频| 成年版毛片免费区| 日韩av在线大香蕉| 日日摸夜夜添夜夜添小说| 免费av毛片视频| 天美传媒精品一区二区| 99热6这里只有精品| 国产精品久久久久久久电影| 免费观看人在逋| 亚洲av电影不卡..在线观看| 亚洲美女黄片视频| 国产男人的电影天堂91| 男女之事视频高清在线观看| 色5月婷婷丁香| 一a级毛片在线观看| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 亚洲一级一片aⅴ在线观看| 国产高清三级在线| 一级a爱片免费观看的视频| 国产亚洲欧美98| 欧美日韩综合久久久久久| 欧美xxxx黑人xx丫x性爽| 免费人成在线观看视频色| 最好的美女福利视频网| 韩国av在线不卡| 国产精品国产高清国产av| 国国产精品蜜臀av免费| 十八禁国产超污无遮挡网站| 神马国产精品三级电影在线观看| 最近在线观看免费完整版| 内射极品少妇av片p| 亚洲自偷自拍三级| 久久久精品大字幕| 久久人人精品亚洲av| 国产精品久久久久久精品电影| 久久久国产成人免费| 我的老师免费观看完整版| 免费观看人在逋| 成人美女网站在线观看视频| 精品人妻视频免费看| 性色avwww在线观看| 国产毛片a区久久久久| 综合色av麻豆| 一夜夜www| 欧美色欧美亚洲另类二区| 变态另类成人亚洲欧美熟女| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 最好的美女福利视频网| 国产69精品久久久久777片| ponron亚洲| 51国产日韩欧美| 久久精品91蜜桃| 成人av在线播放网站| 国产精品久久久久久久电影| 午夜福利视频1000在线观看| 国产精品不卡视频一区二区| 丰满的人妻完整版| 99热6这里只有精品| 亚洲无线观看免费| 一级毛片久久久久久久久女| 欧美潮喷喷水| 色综合站精品国产| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩无卡精品| 国产片特级美女逼逼视频| 国产成人freesex在线 | 亚洲av一区综合| 秋霞在线观看毛片| 日本成人三级电影网站| 女人十人毛片免费观看3o分钟| 国产欧美日韩精品一区二区| 亚洲熟妇中文字幕五十中出| 久久人妻av系列| 深夜a级毛片| 国模一区二区三区四区视频| 网址你懂的国产日韩在线| 国产成人福利小说| 看片在线看免费视频| 村上凉子中文字幕在线| 97热精品久久久久久| 中国美白少妇内射xxxbb| 欧美bdsm另类| 欧美激情在线99| 国产精品电影一区二区三区| 长腿黑丝高跟| 精品99又大又爽又粗少妇毛片| 亚洲电影在线观看av| 最后的刺客免费高清国语| 国产精品99久久久久久久久| av在线天堂中文字幕| 日韩一区二区视频免费看| 长腿黑丝高跟| 亚洲精品乱码久久久v下载方式| 91久久精品电影网| 欧美xxxx黑人xx丫x性爽| 特级一级黄色大片| 女人十人毛片免费观看3o分钟| 国产又黄又爽又无遮挡在线| 久久精品夜夜夜夜夜久久蜜豆| av女优亚洲男人天堂| 精品福利观看| 久久久久九九精品影院| 国产成人一区二区在线| 国产成人91sexporn| 久99久视频精品免费| 九九在线视频观看精品| 夜夜爽天天搞| 中国美白少妇内射xxxbb| 亚洲美女黄片视频| 成人亚洲精品av一区二区| 久久人人精品亚洲av| 亚洲精品456在线播放app| 国产爱豆传媒在线观看| 乱系列少妇在线播放| 美女被艹到高潮喷水动态| 亚洲av中文av极速乱| 最近手机中文字幕大全| 免费av不卡在线播放| 欧美在线一区亚洲| 看免费成人av毛片| 亚洲最大成人av| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 啦啦啦观看免费观看视频高清| 精品人妻偷拍中文字幕| 欧美一区二区亚洲| 午夜福利成人在线免费观看| 成人特级黄色片久久久久久久| 国产欧美日韩一区二区精品|