任敏霞,張子恒,曾月霞,李含笑,王利峰,李會(huì)平,2，劉軍俠,2

(1河北農(nóng)業(yè)大學(xué) 林學(xué)院，河北 保定 071000；2河北省林木種質(zhì)資源與森林保護(hù)重點(diǎn)實(shí)驗(yàn)室，河北 保定 071000)

病蟲害作為限制農(nóng)林業(yè)發(fā)展、造成經(jīng)濟(jì)損失的重要因素，已無法通過傳統(tǒng)方法獲得有效抑制[1]，因此利用轉(zhuǎn)基因技術(shù)將特定基因轉(zhuǎn)入植物體使其獲得抗性來抑制病蟲害的發(fā)生成為當(dāng)前研究的熱點(diǎn)[2]。楊樹作為我國造林綠化重要樹種，具有基因組小、離體操作容易的特點(diǎn)，常作為轉(zhuǎn)基因林木模式樹種進(jìn)行研究。目前研究處于前列的為轉(zhuǎn)(Bt)抗蟲基因楊，自1993年將對(duì)鱗翅目昆蟲具有毒性的Bt基因成功轉(zhuǎn)入歐洲黑楊以來，已成功培育出多種抗蟲基因楊[3]，一部分已進(jìn)入田間試驗(yàn)與商業(yè)化種植階段，具有良好抗蟲效果。

隨著研究的深入，轉(zhuǎn)抗蟲基因楊樹生態(tài)安全性問題引起眾多學(xué)者廣泛關(guān)注[4]。生態(tài)系統(tǒng)中，既包括害蟲又包括天敵，害蟲分為主要害蟲與次要害蟲，而天敵在抑制次要害蟲數(shù)量以及調(diào)節(jié)害蟲種群數(shù)量之中有非常大的作用[5,6]。因此近年來相關(guān)研究人員做了一些關(guān)于轉(zhuǎn)抗蟲基因楊對(duì)主要節(jié)肢動(dòng)物種群影響的研究。目前的研究表明：轉(zhuǎn)不同種類型的抗蟲基因楊對(duì)靶標(biāo)害蟲均表現(xiàn)較高抗蟲性，食葉害蟲數(shù)量明顯減少[4,7,8]，且轉(zhuǎn)基因楊樹抗蟲性會(huì)影響到下一代幼蟲，使其死亡率提高，表明轉(zhuǎn)基因楊樹抗蟲具有持續(xù)性[4,9,10]。對(duì)靶標(biāo)害蟲楊扇舟蛾種群生命表的研究中也證明轉(zhuǎn)基因楊樹楊扇舟蛾存活率與對(duì)照林地相比蟲口密度顯著降低、生長發(fā)育遲緩，表明轉(zhuǎn)基因楊可以減少靶標(biāo)害蟲自然種群數(shù)量，甚至影響未來種群的發(fā)展[11-12]。轉(zhuǎn)基因楊在有效抑制靶標(biāo)害蟲的同時(shí)，對(duì)非靶標(biāo)害蟲及天敵昆蟲具有不同的效應(yīng)[7]。一部分研究表明轉(zhuǎn)基因楊在抑制靶標(biāo)害蟲的同時(shí)，對(duì)于非靶標(biāo)害蟲沒有毒殺作用，且寄生于靶標(biāo)害蟲的天敵種類、數(shù)量和寄生率都高于對(duì)照林地[12]，轉(zhuǎn)基因楊使天敵種群數(shù)量增加。而其他研究表明非靶標(biāo)害蟲種類、數(shù)量具有增多趨勢(shì)，轉(zhuǎn)基因楊對(duì)主要捕食性天敵數(shù)量無影響，而寄生性天敵的數(shù)量則明顯減少[13-15]。

目前轉(zhuǎn)基因楊對(duì)主要節(jié)肢動(dòng)物種群影響的研究尤其是對(duì)非靶標(biāo)生物的影響結(jié)論不盡相同，需對(duì)其繼續(xù)進(jìn)行探究，可為其生態(tài)安全性提供科學(xué)依據(jù)并為以后更好地防治提供方法[16]。因此以轉(zhuǎn)基因107楊為研究對(duì)象，通過研究轉(zhuǎn)基因107楊對(duì)主要節(jié)肢動(dòng)物種群數(shù)量及動(dòng)態(tài)的影響，探究其生態(tài)安全性為進(jìn)一步推廣與應(yīng)用提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)林分為轉(zhuǎn)抗蟲基因107楊(BtCry1Ac和API基因)，篩選出對(duì)美國白蛾、楊扇舟蛾等具有高抗蟲的#1、#2、#3、#9株系[17]。

1.2 試驗(yàn)林概況

試驗(yàn)林地位于河北省灤南縣國有林地苗圃場(chǎng)，屬暖溫帶半濕潤季風(fēng)性氣候，年平均氣溫10.6 ℃，常年降水量658 mm。試驗(yàn)林采用隨機(jī)區(qū)組設(shè)計(jì)的方法，劃分為4個(gè)區(qū)組，每個(gè)區(qū)組種植有#1、#2、#3、#9株系與對(duì)照107楊，每個(gè)株系為1個(gè)小區(qū)，每個(gè)小區(qū)種植30株樹，株距2 m，行間距5 m。試驗(yàn)林地貌氣候一致，管理方式一致，未施用任何化學(xué)農(nóng)藥。

1.3 調(diào)查方法

2019年5-10月，每20 d調(diào)查1次，共調(diào)查8次。采取隨機(jī)選取樣方的方法，每個(gè)區(qū)組隨機(jī)選取5株樣樹，樹冠：每株樣樹(東、西、南、北)4個(gè)方位上下各選取50 cm樣枝，觀察并記錄其上主要節(jié)肢動(dòng)物種類和數(shù)量；主干：調(diào)查距地表3 m內(nèi)主干主要節(jié)肢動(dòng)物種類與數(shù)量并記錄；地面：選取的樣地內(nèi)5個(gè)1 m×1 m的樣方進(jìn)行調(diào)查，記錄主要節(jié)肢動(dòng)物種類與數(shù)量。對(duì)無法識(shí)別的昆蟲放入裝有酒精的指形管中，帶回實(shí)驗(yàn)室并借助工具進(jìn)行鑒定[18-19]。

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

利用SPSS 21.0中單因素(One-Way ANOVA)方差分析對(duì)主要節(jié)肢動(dòng)物種群數(shù)量進(jìn)行顯著差異性分析(P<0.05)，相關(guān)數(shù)據(jù)整理與繪圖在Microsoft Excel 2010中完成。

2 結(jié)果與分析

2.1 轉(zhuǎn)抗蟲基因107楊節(jié)肢動(dòng)物群落組成

轉(zhuǎn)抗蟲基因107楊節(jié)肢動(dòng)物群落組成見表1。

表1 轉(zhuǎn)抗蟲基因107楊節(jié)肢動(dòng)物群落組成Table 1 Composition of arthropod community of transgenic insect-resistant gene 107 poplar

由表1可見，共獲得節(jié)肢動(dòng)物14 299頭，分屬11個(gè)目，105種。物種數(shù)以鱗翅目、半翅目、蜘蛛目、鞘翅目、雙翅目為優(yōu)勢(shì)類群，占比為17.1%、20.0%、19.0%、16.2%。個(gè)體數(shù)量半翅目占優(yōu)勢(shì)地位，占比為29.3%，其次為鱗翅目、膜翅目與蜘蛛目，占比分別為25.4%、18.0%、17.2%。半翅目個(gè)體數(shù)量較多可能是由于轉(zhuǎn)基因107楊對(duì)鱗翅目靶標(biāo)害蟲的控制作用，靶標(biāo)害蟲數(shù)量減少，導(dǎo)致半翅目個(gè)體數(shù)量增多。

2.2 主要靶標(biāo)害蟲個(gè)體數(shù)量及動(dòng)態(tài)變化

2.2.1 主要靶標(biāo)害蟲個(gè)體數(shù)量 主要靶標(biāo)害蟲個(gè)體數(shù)量方差分析見表2。

表2 主要靶標(biāo)害蟲個(gè)體數(shù)量方差分析Table 2 Variance analysis of main target pest in individual number

由表2可見，轉(zhuǎn)基因107美國白蛾、楊扇舟蛾與楊白潛葉蛾個(gè)體數(shù)量顯著低于對(duì)照組(P<0.05)，楊小舟蛾個(gè)體數(shù)量相對(duì)較少，與對(duì)照組差異不顯著(P>0.05)，各株系間主要靶標(biāo)害蟲個(gè)體數(shù)量不存在顯著性差異，#2株系害蟲數(shù)量相對(duì)較少。上述結(jié)果表明轉(zhuǎn)基因107楊對(duì)主要靶標(biāo)害蟲具顯著抗蟲效果。

2.2.2 主要靶標(biāo)害蟲個(gè)體數(shù)量動(dòng)態(tài)變化 美國白蛾、楊扇舟蛾、楊小舟蛾及楊白潛葉蛾個(gè)體數(shù)量動(dòng)態(tài)變化見圖1-圖4。

圖1 美國白蛾個(gè)體數(shù)量動(dòng)態(tài)變化Figure 1 The individual number dynamics of Hyphantria cunea

圖2 楊扇舟蛾個(gè)體數(shù)量動(dòng)態(tài)變化Figure 2 The individual number dynamics of Clostera anachoreta

圖3 楊小舟蛾個(gè)體數(shù)量動(dòng)態(tài)變化Figure 3 The individual number dynamics ofMicromelalopha troglodyta Graeser

圖4 楊白潛葉蛾個(gè)體數(shù)量動(dòng)態(tài)變化Figure 4 The individual number dynamics of Leucopterasusinella Herrich-Schfer

由圖1可見，各株系美國白蛾個(gè)體數(shù)量相對(duì)低于對(duì)照組，對(duì)照組美國白蛾個(gè)體數(shù)量在5月8日-5月29日為個(gè)體數(shù)量急劇增加的時(shí)期，5月29日達(dá)到發(fā)生最高峰，之后急劇下降，又上升，最終處于平穩(wěn)狀態(tài)。#1、#2、#3、#9株系美國白蛾個(gè)體數(shù)量相差不大，處于相對(duì)穩(wěn)定偏低狀態(tài)。由圖2可見，對(duì)照組楊扇舟蛾個(gè)體數(shù)量相對(duì)高于各株系，8月22日-9月31日為其個(gè)體數(shù)量發(fā)生高峰期，#3株系楊扇舟蛾個(gè)體數(shù)量同樣在8月22日-9月31日為發(fā)生高峰期，個(gè)體數(shù)量相對(duì)較多。其他各株系楊扇舟蛾個(gè)體數(shù)量較少一直處于相對(duì)穩(wěn)定狀態(tài)。由圖3可見，對(duì)照組楊小舟蛾個(gè)體數(shù)量相對(duì)高于各株系，對(duì)照組楊小舟蛾在8月上旬為其發(fā)生高峰期，8月下旬至9月初個(gè)體數(shù)量減少，9月下旬至10月初個(gè)體數(shù)量又增多，之后個(gè)體數(shù)量減少趨于平穩(wěn)。#2、#3、#9株系在7月31日個(gè)體數(shù)量開始增多，8月22日達(dá)到發(fā)生最高峰，之后逐漸減少趨于平穩(wěn)，#1株系楊小舟蛾個(gè)體數(shù)量相對(duì)較少，一直處于較低水平。由圖4可見，整體上各株系楊白潛葉蛾個(gè)體數(shù)量低于對(duì)照。對(duì)照組楊白潛葉蛾個(gè)體數(shù)量在6月21日-7月11日相對(duì)較多，7月11日達(dá)到發(fā)生最高峰。各株系楊白潛葉蛾個(gè)體數(shù)量在5月29日-6月21日為發(fā)生高峰期，之后呈逐漸下降趨勢(shì)。

2.3 主要非靶標(biāo)害蟲個(gè)體數(shù)量及動(dòng)態(tài)變化

2.3.1 主要非靶標(biāo)害蟲個(gè)體數(shù)量 主要非靶標(biāo)害蟲個(gè)體數(shù)量方差分析見表3。

由表3可見，主要非靶標(biāo)害蟲葉蟬sp.、皺大球蚧及大青葉蟬個(gè)體數(shù)量與對(duì)照組差異不顯著(P>0.05)，白楊毛蚜個(gè)體數(shù)量顯著低于對(duì)照組(P<0.05)，各株系間葉蟬個(gè)體數(shù)量存在一定差異，#3個(gè)體數(shù)量顯著低于#1、#9，各株系間白楊毛蚜、皺大球蚧及大青葉蟬個(gè)體數(shù)量差異不顯著。表明轉(zhuǎn)基因107楊對(duì)主要的非靶標(biāo)害蟲未產(chǎn)生顯著不利影響。

表3 主要非靶標(biāo)害蟲個(gè)體數(shù)量方差分析Table 3 Variance analysis of main non-target pest in individual number

2.3.2 主要非靶標(biāo)害蟲個(gè)體數(shù)量動(dòng)態(tài)變化 葉蟬sp.、皺大球蚧、白楊毛蚜及大青葉蟬個(gè)體數(shù)量動(dòng)態(tài)變化見圖5-圖8。

圖5 葉蟬sp.個(gè)體數(shù)量動(dòng)態(tài)變化Figure 5 The individual number dynamics ofCicadellidae sp.

圖6 皺大球蚧個(gè)體數(shù)量動(dòng)態(tài)變化Figure 6 The individual number dynamics ofEulecanium kuwanai Kanda

圖7 白楊毛蚜個(gè)體數(shù)量動(dòng)態(tài)變化Figure 7 The individual number dynamics ofChaitophorus populeti

圖8 大青葉蟬個(gè)體數(shù)量動(dòng)態(tài)變化

由圖5可見，#9株系與對(duì)照組葉蟬sp.個(gè)體數(shù)量5月29日為發(fā)生最高峰，5月29日之后，對(duì)照與各株系葉蟬sp.個(gè)體數(shù)量均略有增加，相對(duì)平穩(wěn)。由圖6可見，#1株系、#9株系與對(duì)照皺大球蚧個(gè)體數(shù)量在6月21日為發(fā)生最高峰，9月13日為#1、#9株系與對(duì)照個(gè)體數(shù)量發(fā)生小高峰期。#2株系、#3株系個(gè)體數(shù)量相對(duì)較少，一直處于較低水平狀態(tài)。由圖7可見，各株系白楊毛蚜個(gè)體數(shù)量低于對(duì)照組，各株系及對(duì)照組白楊毛蚜個(gè)體數(shù)量表現(xiàn)為相同變化趨勢(shì)，在5月8日-9月13日個(gè)體數(shù)量一直較少，10月9日以后白楊毛蚜個(gè)體數(shù)量急劇增加，達(dá)到發(fā)生最高峰。由圖8可見，對(duì)照、#3大青葉蟬個(gè)體數(shù)量在7月11日以后逐漸增加，8月上旬達(dá)到發(fā)生最高峰，之后數(shù)量逐漸減少，#9株系在7月11日為發(fā)生最高峰，之后個(gè)體數(shù)量下降，8月上旬個(gè)體數(shù)量又增多，為大青葉蟬個(gè)體數(shù)量發(fā)生小高峰期，#1、#2個(gè)體數(shù)量較少，變化幅度較小，相對(duì)穩(wěn)定。

2.4 主要天敵個(gè)體數(shù)量及動(dòng)態(tài)變化

2.4.1 主要天敵個(gè)體數(shù)量 主要天敵個(gè)體數(shù)量方差分析見表4。

由表4可見，轉(zhuǎn)基因107楊主要天敵異色瓢蟲、狼蛛sp.、圓蛛sp.、跳蛛sp.個(gè)體數(shù)量與對(duì)照組無顯著性差異(P>0.05)，部分株系主要天敵個(gè)體數(shù)量表現(xiàn)為一定程度的增多，表明轉(zhuǎn)基因107楊對(duì)主要天敵未造成不利影響。

表4 主要天敵個(gè)體數(shù)量方差分析Table 4 Variance analysis of main natural enemiesin individual number

2.4.2 主要天敵個(gè)體數(shù)量動(dòng)態(tài)變化 異色瓢蟲、狼蛛sp.、圓蛛sp.及跳蛛sp.個(gè)體數(shù)量動(dòng)態(tài)變化見圖9-圖12。

圖9 異色瓢蟲個(gè)體數(shù)量動(dòng)態(tài)變化

圖10 狼蛛sp.個(gè)體數(shù)量動(dòng)態(tài)變化

圖11 圓蛛sp.個(gè)體數(shù)量動(dòng)態(tài)變化

圖12 跳蛛sp.個(gè)體數(shù)量動(dòng)態(tài)變化

由圖9可見，異色瓢蟲各株系與對(duì)照組個(gè)體數(shù)量相差不大，#1、#3、#9株系與對(duì)照組個(gè)體數(shù)量在6月21日達(dá)到發(fā)生最高峰，9月13日為發(fā)生小高峰，#2株系在5月29日與9月13日個(gè)體數(shù)量相對(duì)較多，為發(fā)生高峰期。由圖10可見，各株系與對(duì)照狼蛛sp.個(gè)體數(shù)量相差不大，對(duì)照組、#1、#2、#3狼蛛sp.個(gè)體數(shù)量在7月11日達(dá)到發(fā)生最高峰，#9株系7月31日達(dá)到發(fā)生最高峰，7月31日之后，各株系與對(duì)照個(gè)體數(shù)量呈逐漸下降趨勢(shì)。由圖11可見，5月8日-5月29日對(duì)照組與#2株系圓蛛sp.個(gè)體數(shù)量相對(duì)較多達(dá)到峰值，之后又逐漸下降，趨于平穩(wěn)。其余各株系圓蛛sp.個(gè)體數(shù)量增幅不大，總體呈偏低、均衡的狀態(tài)。由圖12可見，各株系與對(duì)照組跳蛛sp.個(gè)體數(shù)量變化幅度較大，整體呈現(xiàn)為上升、下降、上升再下降的趨勢(shì)。5-6月份為各株系及對(duì)照個(gè)體數(shù)量發(fā)生高峰期，7、8月份跳蛛sp.個(gè)體數(shù)量下降，相對(duì)較少，之后個(gè)體數(shù)量又表現(xiàn)為上升趨勢(shì)，不同月份間跳蛛個(gè)體數(shù)量波動(dòng)較大。

2.5 主要中性昆蟲個(gè)體數(shù)量及動(dòng)態(tài)變化

2.5.1 主要中性昆蟲個(gè)體數(shù)量 主要中性昆蟲個(gè)體數(shù)量見表5。

表5 主要中性昆蟲個(gè)體數(shù)量方差分析Table 5 Variance analysis of main neutral enemiesin individual number

由表5可見，轉(zhuǎn)基因107楊家蠅、家蟻個(gè)體數(shù)量與對(duì)照組無顯著性差異(P>0.05)，各株系家蠅個(gè)體數(shù)量相對(duì)于對(duì)照組表現(xiàn)為一定程度的增多，表明轉(zhuǎn)基因107楊未對(duì)主要中性昆蟲產(chǎn)生顯著不利影響。

2.5.2 主要中性昆蟲個(gè)體數(shù)量動(dòng)態(tài)變化 家蠅、家蟻個(gè)體數(shù)量動(dòng)態(tài)變化如圖13、圖14所示。

圖13 家蠅個(gè)體數(shù)量動(dòng)態(tài)變化Figure 13 The individual number dynamicsOf Musca domestica

圖14 家蟻個(gè)體數(shù)量動(dòng)態(tài)變化Figure 14 The individual number dynamicsOf Monamorium pharaonis

由圖13可見，各株系與對(duì)照組家蠅個(gè)體數(shù)量變化趨勢(shì)相同，5月8日-5月29日個(gè)體數(shù)量較多，保持相對(duì)較高水平，為發(fā)生高峰期，其余月份各株系與對(duì)照個(gè)體數(shù)量隨月份的變化幅度較大，整體呈現(xiàn)為下降、升高再下降的趨勢(shì)。由圖14可見，對(duì)照組與各株系間家蟻個(gè)體數(shù)量相差不大，但總體呈現(xiàn)的動(dòng)態(tài)趨勢(shì)不穩(wěn)定，CK、#3、#9在5月29日-6月21日為發(fā)生高峰期，6月21日達(dá)到最高峰，#1株系在7月11日-7月31日為發(fā)生高峰期，在7月31日達(dá)到個(gè)體數(shù)量發(fā)生最高峰，不同月份之間家蟻個(gè)體數(shù)量波動(dòng)較大。

3 結(jié)論與討論

通過對(duì)轉(zhuǎn)基因107楊主要節(jié)肢動(dòng)物種群種類、數(shù)量及動(dòng)態(tài)進(jìn)行調(diào)查研究，得出結(jié)論如下：轉(zhuǎn)基因107楊各株系美國白蛾、楊扇舟蛾等靶標(biāo)害蟲個(gè)體數(shù)量明顯減少低于對(duì)照楊，轉(zhuǎn)基因107楊對(duì)主要靶標(biāo)害蟲具有良好抑制效果；主要非靶標(biāo)害蟲葉蟬 sp.、皺大球蚧、大青葉蟬個(gè)體數(shù)量與對(duì)照無明顯差異，白楊毛蚜個(gè)體數(shù)量低于對(duì)照，轉(zhuǎn)基因107楊對(duì)非靶標(biāo)害蟲未產(chǎn)生不利影響；轉(zhuǎn)基因107楊對(duì)主要天敵與中性昆蟲未造成顯著不利影響，部分株系異色瓢蟲、狼蛛sp.、跳蛛sp.及家蠅個(gè)體數(shù)量相對(duì)于對(duì)照表現(xiàn)為一定程度的增多；大部分主要節(jié)肢動(dòng)物種群發(fā)生高峰期集中在5-6月、8-9月，但楊白潛葉蛾7月為發(fā)生高峰期，白楊毛蚜在10月為發(fā)生高峰期。

轉(zhuǎn)基因楊樹為防治蟲害、減少化學(xué)農(nóng)藥使用量方面做出巨大貢獻(xiàn)的同時(shí)[20]其生態(tài)安全性受到廣泛關(guān)注，本研究以轉(zhuǎn)基因107楊為研究對(duì)象，通過研究其對(duì)主要節(jié)肢動(dòng)物種群的影響為轉(zhuǎn)基因楊生態(tài)安全性研究提供依據(jù)。結(jié)果表明：轉(zhuǎn)基因107楊靶標(biāo)害蟲美國白蛾、楊扇舟蛾等個(gè)體數(shù)量明顯減少，轉(zhuǎn)基因107楊對(duì)主要靶標(biāo)害蟲具有顯著的抑制效果，這與張益文等的研究結(jié)果相似[21]。昆蟲之間是由食物網(wǎng)相連接的，轉(zhuǎn)基因楊在抑制靶標(biāo)害蟲的同時(shí)對(duì)其他種群也有可能產(chǎn)生不利影響，因此研究過程中應(yīng)綜合考慮，且試驗(yàn)林地處于開放環(huán)境中，在分析時(shí)也應(yīng)考慮種群的遷移和遷出。研究表明轉(zhuǎn)基因107楊在抑制靶標(biāo)害蟲的同時(shí)，對(duì)非靶標(biāo)害蟲未產(chǎn)生顯著不利影響，這與Zhang等的研究結(jié)果相似[22]。天敵在抑制害蟲方面起到非常重要的作用，利用轉(zhuǎn)抗蟲基因楊樹進(jìn)行抗蟲時(shí)不能對(duì)天敵產(chǎn)生不利影響[23]，因此在對(duì)轉(zhuǎn)抗蟲基因楊樹生態(tài)安全性進(jìn)行評(píng)價(jià)時(shí)，需對(duì)其是否對(duì)天敵產(chǎn)生影響進(jìn)行研究，本研究結(jié)果為轉(zhuǎn)基因107楊對(duì)主要天敵未產(chǎn)生不利影響，部分天敵及中性昆蟲個(gè)體數(shù)量表現(xiàn)為一定程度的增多，表明轉(zhuǎn)基因107楊未產(chǎn)生毒殺作用，未對(duì)天敵及中性昆蟲個(gè)體數(shù)量產(chǎn)生不利影響，這與姚麗等的研究結(jié)果相似[24]。部分株系昆蟲數(shù)量波動(dòng)較大，起伏不定可能是由于采用隨機(jī)取樣調(diào)查的方法、調(diào)查選取的時(shí)間、氣候因素如風(fēng)雨等及其它各種環(huán)境因素的變化對(duì)節(jié)肢動(dòng)物種群的數(shù)量產(chǎn)生影響，導(dǎo)致節(jié)肢動(dòng)物種群產(chǎn)生波動(dòng)[25]。

轉(zhuǎn)基因楊作為多年生林木，研究其對(duì)主要節(jié)肢動(dòng)物種群的影響需對(duì)其進(jìn)行高頻率的調(diào)查與長期調(diào)監(jiān)測(cè)，為科學(xué)評(píng)價(jià)轉(zhuǎn)基因楊的生態(tài)安全性提供依據(jù)。