CN106093369B - Bt水稻对非靶标害虫褐飞虱的安全性评价方法 - Google Patents
Bt水稻对非靶标害虫褐飞虱的安全性评价方法 Download PDFInfo
- Publication number
- CN106093369B CN106093369B CN201610492939.3A CN201610492939A CN106093369B CN 106093369 B CN106093369 B CN 106093369B CN 201610492939 A CN201610492939 A CN 201610492939A CN 106093369 B CN106093369 B CN 106093369B
- Authority
- CN
- China
- Prior art keywords
- age
- mrow
- paddy rice
- msub
- plant hopper
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 153
- 235000009566 rice Nutrition 0.000 title claims abstract description 153
- 241001498622 Cixius wagneri Species 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 23
- 241000238631 Hexapoda Species 0.000 title claims abstract description 22
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 153
- 230000004083 survival effect Effects 0.000 claims abstract description 35
- 230000035558 fertility Effects 0.000 claims abstract description 28
- 238000011161 development Methods 0.000 claims abstract description 15
- 230000009261 transgenic effect Effects 0.000 claims description 55
- 241000607479 Yersinia pestis Species 0.000 claims description 15
- 235000013601 eggs Nutrition 0.000 claims description 12
- 238000009395 breeding Methods 0.000 claims description 8
- 230000001488 breeding effect Effects 0.000 claims description 8
- 239000011521 glass Substances 0.000 claims description 8
- 230000001850 reproductive effect Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 101150065438 cry1Ab gene Proteins 0.000 claims description 5
- 230000035755 proliferation Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000010835 comparative analysis Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 abstract 1
- 235000013305 food Nutrition 0.000 description 14
- 230000018109 developmental process Effects 0.000 description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 108010000912 Egg Proteins Proteins 0.000 description 7
- 102000002322 Egg Proteins Human genes 0.000 description 7
- 210000004681 ovum Anatomy 0.000 description 7
- 241000819999 Nymphes Species 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 3
- 241001556089 Nilaparvata lugens Species 0.000 description 3
- 244000037671 genetically modified crops Species 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000255777 Lepidoptera Species 0.000 description 2
- 241000700605 Viruses Species 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 230000001418 larval effect Effects 0.000 description 2
- 230000017448 oviposition Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000426497 Chilo suppressalis Species 0.000 description 1
- 241001414720 Cicadellidae Species 0.000 description 1
- 241000008892 Cnaphalocrocis patnalis Species 0.000 description 1
- 241001466044 Delphacidae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 241001249129 Scirpophaga incertulas Species 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000008144 egg development Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/5005—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
- G01N33/5008—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics
- G01N33/5082—Supracellular entities, e.g. tissue, organisms
- G01N33/5085—Supracellular entities, e.g. tissue, organisms of invertebrates
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Immunology (AREA)
- Hematology (AREA)
- Cell Biology (AREA)
- Chemical & Material Sciences (AREA)
- Urology & Nephrology (AREA)
- Molecular Biology (AREA)
- Toxicology (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Tropical Medicine & Parasitology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Microbiology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Pretreatment Of Seeds And Plants (AREA)
- Catching Or Destruction (AREA)
Abstract
本发明涉及一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法。该方法通过年龄‑龄期两性生命表对褐飞虱的生态适合度进行研究,从而阐明Bt水稻对褐飞虱的生态安全性风险;具体包括五个部分:(1)发育历期和寿命;(2)年龄‑龄期存活率;(3)生殖力和种群参数;(4)寿命期望和生殖价值;(5)年龄分布和种群预测;通过五个方面的比较分析综合评价Bt水稻对非靶标害虫褐飞虱的生态风险。本发明在实验室内利用两性种群生命表建立了***评价Bt水稻非靶标效应的评价方法,对完善Bt水稻对非靶标生物的影响研究以及确定Bt水稻种植的生态安全性具有重要意义。
Description
技术领域
本发明属于转基因作物生态安全风险评价技术领域,尤其是一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法。
背景技术
水稻是我国的重要粮食作物,在粮食生产占据着至关重要的战略地位(夏如兵,2009),同时水稻也是虫害最严重的粮食作物(张文辉和刘光杰,2001)。近年来由于全球气候变暖,二化螟、三化螟、大螟和稻纵卷叶螟等鳞翅目害虫对水稻的为害程度日益严重,给农民带来了严重的经济损失(韩永强等,2008;雷惠质和江建云,1996)。利用基因工程技术培育的转Bt抗虫水稻能够有效地防治鳞翅目害虫,其表达的Bt毒蛋白对靶标害虫有极高的抗性(蔡万伦等,2008;肖汉祥等,2011;徐雪亮等,2013)。然而,靶标害虫生态位的缺失,是否会引起非靶标害虫(叶蝉类,飞虱类等)种群密度的上升,从而由次要害虫转变为主要害虫(Lu et al,2010)。截止目前,全球已有20多个国家种植转基因作物(James,2015)。为确保转基因作物的安全使用,针对非靶标生物的生态安全风险评估是必要的(Yu et al,2011)。
褐飞虱[Nilaparvata lugens()],属半翅目飞虱科,是亚洲水稻上的一种次要的单食性害虫,主要通过刺吸式口器直接吸取稻株韧皮部的植物汁液为害水稻,并且能传播水稻病毒病(王荣富等,1998)。由于转Bt抗虫水稻表达的Bt蛋白的主要靶标对象是鳞翅目昆虫,褐飞虱属于非靶标害虫。生命表一直是研究昆虫生态学的重要工具。然而传统雌性年龄生命表由于忽略龄期重叠、不考虑雄性个体及不包含成虫前期死亡率,从而使得到的种群参数存在偏差,不能确切反映昆虫种群的实际动态(Chi and Liu,1985;Huang andChi,2012;吴坤君等,1994)。为了改进和补充传统生命表的不足、不完整,Chi和Liu(1985)提出了两性生命表理论,并在种群生态学研究中得到广泛应用(Li et al,2015;Tuan etal,2016;Wang et al,2016;Yin et al,2013)。
发明内容
本发明的目的在于提供一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法,该方法通过构建褐飞虱在Bt水稻和对照非转基因水稻上的实验种群年龄-龄期两性生命表,从种群生态学角度***地研究Bt水稻对非靶标害虫褐飞虱生长发育、寿命以及繁殖的影响,旨在提供一种***、准确、简便、科学的转基因水稻对非靶标害虫褐飞虱的安全性评价方法。
为实现上述目的,本发明的技术方案是:一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法,包括如下步骤,
S1:计算并记录褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命以及生殖力的参数;
S2:计算褐飞虱在Bt水稻及其对照非转基因水稻上的年龄-龄期存活率;
S3:计算褐飞虱在Bt水稻及其对照非转基因水稻上的生殖力和种群动态参数;
S4:计算褐飞虱在Bt水稻及其对照非转基因水稻上的寿命期望和生殖价值;
S5:分析褐飞虱在Bt水稻及其对照非转基因水稻上的稳定年龄-龄期分布和种群预测。
在本发明一实施例中,所述Bt水稻选择转cry1Ab基因抗虫水稻mf-MH3301-1,所述对照非转基因水稻为对照亲本非转基因水稻MH3301。
在本发明一实施例中,所述步骤S1中是通过Bootstrap方法来计算褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命以及生殖力的参数。
在本发明一实施例中,所述S2-S5采用年龄-龄期两性生命表方法,具体即通过年龄-龄期两性生命表计算机程序Twosex-MS Chart对步骤S1记录参数进行分析获得年龄-龄期存活率、生殖力、种群动态参数、寿命期望、生殖价值、稳定年龄-龄期分布和种群预测。
在本发明一实施例中,所述褐飞虱在Bt水稻及其对照非转基因水稻上的种群动态参数计算公式如下:
λ=er;GRR=Σmx
其中,n是褐飞虱全龄期数,x和i为年龄,j为龄期;lx指整个种群存活到年龄x的概率;Sxj指新孵化个体存活到x年龄j龄期的概率;mx指种群在年龄x的平均产卵量;fxj指雌成虫在年龄x的平均产卵量;R0为净增殖率;r为内禀增长率;λ为周限增长率;T为世代平均周期;GRR为总繁殖力GRR;Rx为年龄特征增殖率。
在本发明一实施例中,褐飞虱在玻璃试管中单头饲养,该玻璃试管直径3cm高度20cm。
相较于现有技术,本发明具有以下有益效果:本发明通过Bootstrap(n=100000)方法来计算褐飞虱的发育历期、寿命、生殖力、种群数量,以及种群动态参数等的标准误(SE),并通过程序Twosex-MS Chart的paired bootstrap test功能来比较取食两种水稻的褐飞虱发育历期、生殖、寿命、种群参数等的差异。
附图说明
图1为本发明方法流程图。
图2是褐飞虱在Bt水稻(A)和非转基因水稻(B)上的年龄-龄期存活率曲线(sxj)。
图3是在Bt水稻(A)和非转基因水稻(B)上褐飞虱的种群年龄特征存活率(lx)、年龄-龄期特征生殖力(fxj)、种群年龄特征生殖力(mx)、种群年龄特征生殖值(lxmx)。
图4是褐飞虱在Bt水稻(A)和非转基因水稻(B)上的寿命期望(exj)。
图5是褐飞虱在Bt水稻(A)和非转基因水稻(B)上的生殖价值(vxj)。
图6是褐飞虱在Bt水稻(A)和非转基因水稻(B)上的稳定年龄分布。
图7是褐飞虱在Bt水稻(A)和非转基因水稻(B)上的种群大小预测。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
如图1所示,本发明的一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法,包括如下步骤,
S1:计算并记录褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命以及生殖力的参数;
S2:计算褐飞虱在Bt水稻及其对照非转基因水稻上的年龄-龄期存活率;
S3:计算褐飞虱在Bt水稻及其对照非转基因水稻上的生殖力和种群动态参数;
S4:计算褐飞虱在Bt水稻及其对照非转基因水稻上的寿命期望和生殖价值;
S5:分析褐飞虱在Bt水稻及其对照非转基因水稻上的稳定年龄-龄期分布和种群预测。
所述Bt水稻选择转cry1Ab基因抗虫水稻mf-MH3301-1,所述对照非转基因水稻为对照亲本非转基因水稻MH3301。
所述步骤S1中是通过Bootstrap方法来计算褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命和生殖力的参数。
所述S2-S5采用年龄-龄期两性生命表方法,具体即通过年龄-龄期两性生命表计算机程序Twosex-MS Chart对步骤S1记录参数进行分析获得年龄-龄期存活率、生殖力、种群动态参数、寿命期望、生殖价值、稳定年龄-龄期分布和种群预测。
所述褐飞虱在Bt水稻及其对照非转基因水稻上的种群动态参数计算公式如下:
λ=er;GRR=Σmx
其中,n是褐飞虱全龄期数,x和i为年龄,j为龄期;lx指整个种群存活到年龄x的概率;Sxj指新孵化个体存活到x年龄j龄期的概率;mx指种群在年龄x的平均产卵量;fxj指雌成虫在年龄x的平均产卵量;R0为净增殖率;r为内禀增长率;λ为周限增长率;T为世代平均周期;GRR为总繁殖力GRR;Rx为年龄特征增殖率。
褐飞虱在玻璃试管中单头饲养,该玻璃试管直径3cm高度20cm。
以下为本发明的具体实施例。
本发明的供试水稻为转cry1Ab基因抗虫水稻mf-MH3301-1和对照亲本非转基因水稻MH3301,水稻品种均由福建省农科院农业遗传工程重点实验室提供。
褐飞虱采自福建农林大学植物病毒研究所,后在福建农林大学应用生态研究所气候箱(温度28士1℃,RH 80%士l%,光照L:D=16:8)中分别用转cry1Ab基因抗虫水稻及其对照非转基因水稻连续饲养多年至今。
分别取45日龄两种水稻苗10株分别放入10个大玻璃试管中水培,每个玻璃试管接入1头怀卵雌褐飞虱,用100目的纱布密封,在人工气候箱(温度28士1℃,RH 80%士l%,光照L:D=16/8)产卵24h后,将褐飞虱移走。每日定时观察并记录卵的孵化情况。选取刚孵化的若虫分别用对应的转基因和非转基因水稻苗单管(直径3cm高度20cm)饲养,每个处理69头褐飞虱,供试水稻苗每隔10d更换一次,水培液每隔2d更换1次,每天定时记录褐飞虱蜕皮、羽化、死亡等情况。等褐飞虱羽化12h内取对应水稻品种上的褐飞虱成虫配对,分别记录其每日产卵量并更换产卵的水稻苗,直至雌虫死亡。当褐飞虱成虫全部死亡后整个实验结束。
褐飞虱所有个体的发育历期、存活情况、雌虫日产卵量等原始数据通过年龄-龄期两性生命表计算机程序Twosex-MSChart(Chi,2015)来进行分析,从而获得许多生命表参数。如:年龄-龄期特征存活率Sxj,年龄特征存活率lx,年龄特征生殖力mx,fxj年龄-龄期特征生殖力,年龄-龄期特征寿命期望exj,年龄-龄期特征生殖价值vxj(x=年,j=龄期),内禀增长率r,周限增长率λ,净增殖率R0,世代平均周期T,总繁殖力GRR,产卵前期APOP,总产卵前期TPOP。
其中,lx指整个种群(包括雌、雄褐飞虱及未发育至成虫死亡的所有个体)存活到年龄x的概率,sxj指新孵化个体存活到x年龄j龄期的概率。fxj指雌成虫在年龄x的平均产卵量,mx指种群在年龄x的平均产卵量。exj指种群个体到x年龄j龄期的预期寿命,vxj指种群个体在x年龄j龄期对下一代的平均贡献值。产卵前期指从褐飞虱羽化到第一次产卵的时间,总产卵前期指从卵发育到成虫开始产卵的时间。
种群动态参数计算公式如下: λ=er;GRR=Σmx;
其中,n是褐飞虱全龄期数,x和i为年龄,j为龄期;lx指整个种群存活到年龄x的概率;Sxj指新孵化个体存活到x年龄j龄期的概率;mx指种群在年龄x的平均产卵量;fxj指雌成虫在年龄x的平均产卵量;R0为净增殖率;r为内禀增长率;λ为周限增长率;T为世代平均周期;GRR为总繁殖力GRR;Rx为年龄特征增殖率。
(1)褐飞虱在Bt水稻及其对照上的发育历期、寿命以及生殖力等参数
由表1所示,除一龄和五龄外,褐飞虱各龄期的发育时间在Bt水稻mf-MH3301-1和非转基因对照MH3301没有显著差异(P>0.05)。其中褐飞虱一龄在Bt水稻mf-MH3301-1和非转基因对照MH3301的发育时间分别为2.43天和2.77天,差异达到极显著水平(P<0.01)。褐飞虱五龄在Bt水稻mf-MH3301-1和非转基因对照MH3301的发育时间分别为3.82天和4.20天,达到显著差异(P<0.05)。褐飞虱在Bt水稻mf-MH3301-1和非转基因对照MH3301上的雌虫寿命、雄虫寿命和总体寿命均没有显著差异(P>0.05)。雌褐飞虱在Bt水稻mf-MH3301-1上的生殖力显著大于非转基因对照MH3301(P<0.01);而产卵前期和总产卵前期却没有显著差异(P>0.05)。
表1
表中数据是平均值±标准误,同一行数据后不同字母分别表示经pairedbootstrap test检验差异显著。
(2)褐飞虱在Bt水稻和非转基因对照上的年龄-龄期存活率
年龄-龄期存活率sxj指新孵化个体存活到x年龄j龄期的概率。由于昆虫个体之间以及雌雄之间发育速率的差异,因此存在明显的重叠现象。年龄-龄期存活率曲线(sxj)不仅描述了不同龄期褐飞虱在两种水稻上随时间变化的存活情况,并且能够反映龄期分化和龄期重叠现象,所以年龄-龄期存活率曲线(sxj)更准确的描述了取食两种水稻的褐飞虱种群发育以及存活率的变化情况。从卵期到成虫之前,褐飞虱在两种水稻上的存活率曲线基本相似。然而在成虫阶段,取食Bt水稻mf-MH3301-1的褐飞虱雄虫的存活率明显低于雌虫;取食非转基因对照MH3301的褐飞虱雌雄虫的存活率并无显著差异。当第五龄若虫在26d存活率降为零时,褐飞虱成虫的存活率达到最大值(图2,其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301)。
(3)褐飞虱在转基因抗虫水稻及其对照上的生殖力和种群动态参数
年龄特征存活率lx指整个种群(包括雌、雄褐飞虱及未发育至成虫死亡的所有个体)存活到年龄x的概率。与取食非转基因水稻MH3301相比,随着年龄增加取食Bt水稻mf-MH3301-1的褐飞虱年龄特征存活率(lx)下降速度明显比较缓慢。年龄-龄期特征生殖力fxj指雌成虫在年龄x的平均产卵量,年龄特征生殖力mx指种群在年龄x的平均产卵量。取食Bt水稻mf-MH3301-1的褐飞虱年龄-龄期特征生殖力(fxj)和种群年龄特征生殖力(mx)显著高于对照MH3301,因此取食Bt水稻mf-MH3301-1的褐飞虱的种群年龄特征生殖值(lxmx)也显著高于对照MH3301。当褐飞虱在Bt水稻上取食时,第43d达到最大生殖力13.5;然而在非转基因水稻上取食时,第49d才达到最大生殖力10,延缓了6d(图3,其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301)。
由于总生殖力GRR=∑mx,净增殖率因此取食mf-MH3301-1的褐飞虱的总生殖力和净增殖率显著高于对照MH3301,是其2倍多(P<0.05)。同理,取食mf-MH3301-1的褐飞虱的内禀增长率r和周限增长率λ也显著大于对照MH3301(P<0.05)。而褐飞虱在两种水稻上的世代平均周期T并无显著差异(P>0.05)(表2)。
表2
注:表中数据是平均值±标准误,同一行数据后不同字母分别表示经pairedbootstrap test检验差异显著。
(4)褐飞虱在两种水稻上的寿命期望和生殖价值
年龄-龄期寿命期望exj指种群个体到x年龄j龄期的预期寿命。考虑到性别以及个体之间发育速率的差异,利用褐飞虱在两种水稻上的寿命期望曲线,可以看到不同性别以及不同龄期寿命期望的变化。卵期初始零点的寿命期望(e0,1)即为种群寿命的值,这和表1中的寿命相同。随着时间的增加,褐飞虱在两种水稻的寿命期望均逐渐降低。但是,在Bt水稻上褐飞虱雌虫的寿命期望高于雄虫;而在非转基因对照上却恰恰相反,褐飞虱雄虫的寿命期望高于雌虫(图4,其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301)。
年龄-龄期生殖价值vxj指种群个体在x年龄j龄期对下一代的平均贡献值。卵期在零点对应的生殖价值(v0,1)即为褐飞虱种群在两种水稻上的周限增长率。随着龄期的增加,褐飞虱个体的对后代的生殖价值越来越高,在雌虫阶段达到最大值。在Bt水稻和非转基因对照上,最大生殖价值分别为(v37,7=50.89)和(v35,7=36.05)。其中,褐飞虱在第五龄和雌成虫阶段,在Bt水稻上取食褐飞虱的生殖价值高于非转基因对照(图5,其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301)。
(5)褐飞虱在两种水稻上的稳定年龄-龄期分布和种群预测
如图6所示(其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301),根据年龄-龄期两性生命表数据分别求得在理想条件下(食物、空间不受限制,无天敌干扰)褐飞虱在Bt水稻mf-MH3301-1和非转基因水稻MH3301上的稳定年龄-龄期分布。其中,褐飞虱在Bt水稻上的最终稳定年龄分布是:卵期=60.42%,若虫期=34.00%成虫期=5.58%。褐飞虱在非转基因水稻上的最终稳定年龄分布是:卵期=52.26%,若虫期=39.48%,成虫期=8.26%。
如图7所示(其中A Bt水稻mf-MH-3301-1,B非转基因水稻MH3301),褐飞虱在Bt水稻mf-MH3301-1和非转基因水稻MH3301上的种群增长速率预测。通过年龄-龄期两性生命表构建的种群预测很容易描述出龄期结构的变化趋势。由于年龄-龄期两性生命表可以描述龄期分化,褐飞虱每个龄期都可以显示出来。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (5)
1.一种Bt水稻对非靶标害虫褐飞虱的安全性评价方法,其特征在于:包括如下步骤,
S1:计算并记录褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命以及生殖力的参数;
S2:计算褐飞虱在Bt水稻及其对照非转基因水稻上的年龄-龄期存活率;
S3:计算褐飞虱在Bt水稻及其对照非转基因水稻上的生殖力和种群动态参数;
S4:计算褐飞虱在Bt水稻及其对照非转基因水稻上的寿命期望和生殖价值;
S5:分析褐飞虱在Bt水稻及其对照非转基因水稻上的稳定年龄-龄期分布和种群预测;
所述褐飞虱在Bt水稻及其对照非转基因水稻上的种群动态参数计算公式如下:
<mrow>
<msub>
<mi>l</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>s</mi>
<mrow>
<mi>x</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>;</mo>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>s</mi>
<mrow>
<mi>x</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>x</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>s</mi>
<mrow>
<mi>x</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>l</mi>
<mi>x</mi>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
<mo>;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msub>
<mi>l</mi>
<mi>x</mi>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>;</mo>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<msup>
<mi>e</mi>
<mi>r</mi>
</msup>
<mo>;</mo>
<mi>T</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>lnR</mi>
<mn>0</mn>
</msub>
</mrow>
<mi>r</mi>
</mfrac>
<mo>;</mo>
<mi>G</mi>
<mi>R</mi>
<mi>R</mi>
<mo>=</mo>
<mo>&Sigma;</mo>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
<mo>;</mo>
<mi>R</mi>
<mi>x</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>x</mi>
</munderover>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mrow>
其中,n是褐飞虱全龄期数,x和i为年龄,j为龄期;lx指整个种群存活到年龄x的概率;Sxj指新孵化个体存活到x年龄j龄期的概率;mx指种群在年龄x的平均产卵量;fxj指雌成虫在年龄x的平均产卵量;R0为净增殖率;r为内禀增长率;λ为周限增长率;T为世代平均周期;GRR为总繁殖力GRR;Rx为年龄特征增殖率。
2.根据权利要求1所述的Bt水稻对非靶标害虫褐飞虱的安全性评价方法,其特征在于:所述Bt水稻选择转cry1Ab基因抗虫水稻mf-MH3301-1,所述对照非转基因水稻为对照亲本非转基因水稻MH3301。
3.根据权利要求1所述的Bt水稻对非靶标害虫褐飞虱的安全性评价方法,其特征在于:所述步骤S1中是通过Bootstrap方法来计算褐飞虱在Bt水稻及其对照非转基因水稻上包括发育历期、寿命以及生殖力的参数。
4.根据权利要求1所述的Bt水稻对非靶标害虫褐飞虱的安全性评价方法,其特征在于:所述S2-S5采用年龄-龄期两性生命表方法,具体即通过年龄-龄期两性生命表计算机程序Twosex-MS Chart对步骤S1记录参数进行分析获得年龄-龄期存活率、生殖力、种群动态参数、寿命期望、生殖价值、稳定年龄-龄期分布和种群预测。
5.根据权利要求1所述的Bt水稻对非靶标害虫褐飞虱的安全性评价方法,其特征在于:褐飞虱在玻璃试管中单头饲养,该玻璃试管直径3cm高度20cm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492939.3A CN106093369B (zh) | 2016-06-29 | 2016-06-29 | Bt水稻对非靶标害虫褐飞虱的安全性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492939.3A CN106093369B (zh) | 2016-06-29 | 2016-06-29 | Bt水稻对非靶标害虫褐飞虱的安全性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093369A CN106093369A (zh) | 2016-11-09 |
CN106093369B true CN106093369B (zh) | 2017-10-20 |
Family
ID=57213716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610492939.3A Expired - Fee Related CN106093369B (zh) | 2016-06-29 | 2016-06-29 | Bt水稻对非靶标害虫褐飞虱的安全性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093369B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156187A (zh) * | 2010-12-10 | 2011-08-17 | 中国农业大学 | 一种检测转抗虫基因玉米对非靶标生物是否安全的方法 |
CN102577897A (zh) * | 2012-02-16 | 2012-07-18 | 扬州大学 | 基于农药介入后转Bt基因水稻对褐飞虱生态安全性评价的方法 |
CN102986462A (zh) * | 2012-12-03 | 2013-03-27 | 南京大学 | 一种评价转基因棉花对非靶标生物影响的检测方法 |
CN103502458A (zh) * | 2011-02-17 | 2014-01-08 | 宾夕法尼亚大学托管会 | 用于改变组织特异性和改善aav9-介导的基因转移的组合物和方法 |
CN103875608A (zh) * | 2014-03-24 | 2014-06-25 | 华中农业大学 | 转基因抗虫水稻对捕食性天敌黑肩绿盲蝽安全性评价方法 |
CN103931563A (zh) * | 2014-04-15 | 2014-07-23 | 华中农业大学 | 转基因抗虫水稻对捕食性天敌草间钻头蛛安全性评价方法 |
CN103931564A (zh) * | 2014-04-15 | 2014-07-23 | 华中农业大学 | 转基因抗虫水稻对寄生性天敌稻虱缨小蜂安全性评价方法 |
CN104133054A (zh) * | 2014-07-04 | 2014-11-05 | 中国农业科学院植物保护研究所 | 一种评价转基因棉花对大草蛉影响的生物测定方法 |
-
2016
- 2016-06-29 CN CN201610492939.3A patent/CN106093369B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156187A (zh) * | 2010-12-10 | 2011-08-17 | 中国农业大学 | 一种检测转抗虫基因玉米对非靶标生物是否安全的方法 |
CN103502458A (zh) * | 2011-02-17 | 2014-01-08 | 宾夕法尼亚大学托管会 | 用于改变组织特异性和改善aav9-介导的基因转移的组合物和方法 |
CN102577897A (zh) * | 2012-02-16 | 2012-07-18 | 扬州大学 | 基于农药介入后转Bt基因水稻对褐飞虱生态安全性评价的方法 |
CN102986462A (zh) * | 2012-12-03 | 2013-03-27 | 南京大学 | 一种评价转基因棉花对非靶标生物影响的检测方法 |
CN103875608A (zh) * | 2014-03-24 | 2014-06-25 | 华中农业大学 | 转基因抗虫水稻对捕食性天敌黑肩绿盲蝽安全性评价方法 |
CN103931563A (zh) * | 2014-04-15 | 2014-07-23 | 华中农业大学 | 转基因抗虫水稻对捕食性天敌草间钻头蛛安全性评价方法 |
CN103931564A (zh) * | 2014-04-15 | 2014-07-23 | 华中农业大学 | 转基因抗虫水稻对寄生性天敌稻虱缨小蜂安全性评价方法 |
CN104133054A (zh) * | 2014-07-04 | 2014-11-05 | 中国农业科学院植物保护研究所 | 一种评价转基因棉花对大草蛉影响的生物测定方法 |
Non-Patent Citations (1)
Title |
---|
农药对转Bt基因水稻生理生化及褐飞虱种群增长的影响;石朝鹏;《中国博士学位论文全文数据库》;20140415;第33-46页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106093369A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103931564B (zh) | 转基因抗虫水稻对寄生性天敌稻虱缨小蜂安全性评价方法 | |
CN103875608B (zh) | 转基因抗虫水稻对捕食性天敌黑肩绿盲蝽安全性评价方法 | |
CN102172231A (zh) | 烟蚜越冬保种繁育方法 | |
CN102640733B (zh) | 一种调控烟蚜茧蜂滞育的方法 | |
CN102939943A (zh) | 烟蚜的人工繁殖方法 | |
CN104026078B (zh) | 一种散白蚁室内集约化饲养方法 | |
CN104322452A (zh) | 转Bt基因抗虫水稻对捕食性天敌青翅蚁形隐翅甲安全性评价方法 | |
CN103392595B (zh) | 一种马铃薯脱毒苗试管薯的培育方法 | |
CN104255669B (zh) | 红肩瓢虫滞育诱导方法 | |
CN106093369B (zh) | Bt水稻对非靶标害虫褐飞虱的安全性评价方法 | |
CN103766324A (zh) | 冬枣园利用红颈常室茧蜂生物防治绿盲蝽的方法 | |
CN108935341A (zh) | 一种菜粉蝶绒茧蜂的饲养方法 | |
CN112166912A (zh) | 一种草坪病虫害防治方法 | |
CN104054634B (zh) | 一种云实蛀虫的人工繁殖技术 | |
CN106577381A (zh) | 一种吉富罗非鱼抗无乳链球菌病优良品系的选育方法 | |
CN103125451A (zh) | 一种扩繁伞群追寄蝇的方法 | |
CN103931563B (zh) | 转基因抗虫水稻对捕食性天敌草间钻头蛛安全性评价方法 | |
CN104542397A (zh) | 一种罗非鱼多元亲本杂交制种方法 | |
CN102106319A (zh) | 一种利用梨小食心虫卵室内繁殖暗黑赤眼蜂的方法 | |
CN112106733B (zh) | 一种稻飞虱室内生物测定方法 | |
CN102499183A (zh) | 一种暗黑鳃金龟室内饲养方法 | |
CN104152533B (zh) | 一种烟草小苗期快速鉴定黑胫病抗性的方法 | |
CN102257925B (zh) | 秦岭冷杉保护方法 | |
CN112219626A (zh) | 一种有害生物的绿色防控方法 | |
CN208016721U (zh) | 一种尺蠖类幼虫大量饲养的装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171020 |