Cry1C rice doesn’t affect the ecological fitness of rice brown planthopper, Nilaparvata lugens either under RDV stress or not

The potential risks of Bt rice on non-target arthropods (NTAs) should be evaluated and defined before commercial production. Recently, effects of Bt rice on NTAs under abiotic and biotic stress conditions attracted much attention. Here we reported the effects of Bt rice T1C-19 (Cry1C rice) on the non-target herbivore, Nilaparvata lugens (rice brown planthopper, BPH) with or without RDV (rice dwarf virus) infection conditions. BPH showed no feeding and oviposition preference between Bt rice T1C-19 and its non-Bt parental rice Minghui 63 (MH63), as well as between RDV-infected and RDV-free rice plants. Meanwhile, rice type, RDV infection status, and their interaction had little impacts on the survival, development and fecundity of BPH. By comparison with non-Bt control, Bt rice T1C-19 with or without RDV infection had no significant effects on the life-table parameters of BPH including rm, R0, T, DT and λ. Thus, it could be concluded that Bt rice T1C-19 doesn’t affect the ecological fitness of BPH either under RDV stress or not.

Recently, effects of Bt crops on indirect multi-trophic interactions received much attention from researchers as complicated interactions are common in real environment. It was reported that reduced level of Bt-sensitive lepidopteran caterpillars-induced terpenoids in Bt cotton benefited the non-target herbivore Aphis gossypii (Hemiptera: Aphididae) 13 . Changes in plant terpenoid levels as a response to Heliothis virescens damage might also influence the host quality of aphids for the parasitoid Lysiphlebus testaceipes (Hymenoptera: Aphidiidae). Fewer mummies of L. testaceipes were formed on H. virescens infested non-Bt cotton than on Bt cotton 14 . While caterpillar infested Bt rice attracted less BPH due to a relative lower level of plant volatiles released 15 . Plant pathogens including bacteria, fungi, virus would also interact with plant internal metabolism, impair plant growth and reproduction, and thus influence the behavior and performance of vector or non-vector insects. The non-vector insect Bemisia tabaci preferred to settle and oviposit on Tomato spotted wilt virus (TSWV)-infected Datura stramonium due to the alteration of plant volatiles 16 . Rice black streak dwarf virus (RBSDV) infection induced the increment of amino acids and soluble sugar in rice plants, and these changed metabolites improved the ecological fitness of the non-vector BPH 17 . There were also a few reports on combined effects of Bt rice and pathogen attacks on vector or non-vector insects. Southern rice black-streaked dwarf virus (SRBSDV)-infected T1C-19 or T2A-1 had little effects on the ecological adaptability of the non-vector insect BPH 18 . Main biological parameters of vector insect GRLHs such as egg hatching rate, nymph survival and female fecundity were also not affected by RDV (rice dwarf virus)-infected T1C-19 and T2A-1 19 .
Rice dwarf virus (RDV), a member of the genus Phytoreovirus in the family Reoviridae, is mainly transmitted by GRLHs in a persistent-propagative manner 19 . GRLHs acquire the virus by feeding on RDV-infected rice plants for a few minutes to several days. RDV could proliferates in the nearly 70% of GRLH individuals, which became viruliferous after a latent period of 2-3 weeks 20 . Viruliferous GRLHs can transmit RDV to rice plants throughout the lifetime. Infection with RDV would cause significant phenotype changes in rice plants such as dwarfism, white chlorotic specks on leaves, and incomplete panicle exsertion [21][22][23] . Here, the combined influences of Bt rice and RDV infection on an important insect BPH were investigated, which can transmit plant viruses such as rice ragged stunt virus (RRSV) and rice grassy stunt virus (RGSV), but not RDV 22,[24][25][26] . The feeding and oviposition preference of BPH were investigated between Bt rice T1C-19 and its non-Bt parent rice MH63, as well as those between RDV-infected and RDV-free rice plants. Life-table parameters were used to evaluate the ecological fitness of BPH on both RDV-infected and RDV-free T1C-19 or MH63 rice plants. The present study would fulfill the multi-factor interactions between Bt rice and other biotic factors on NTAs, which can help to predict the sustainability of Bt rice in the system of rice pest integrated management.

Results
Feeding and oviposition preference of BPH. Feeding and oviposition preference of BPH were investigated in homemade plastic cages in pairs independently (MH63 vs T1C-19, MH63-RDV vs T1C-19-RDV, MH63 vs MH63-RDV, T1C-19 vs T1C-19-RDV). For feeding preference, the number of BPH female adults settled on each rice type was calculated within 72 h post-inoculation (PI) in each test. Number of BPH on T1C-19 did not differ from that on MH63 either under RDV-free or RDV infection status (Fig. 1a, P > 0.05). Meanwhile, number of BPH on RDV-infected T1C-19 or MH63 showed no significant difference with that on RDV-free rice plants either (Fig. 1b, P > 0.05). Thus, rice type and RDV infection status had no impact on BPH feeding preference within 72 h PI. For oviposition preference, the number of egg mass and total number of eggs per female adult were counted on each rice plants at 72 h PI. No significant difference was seen in the two parameters between BPH on different rice types both under RDV-free and RDV infection status (Fig. 2, P > 0.05). Similarly, www.nature.com/scientificreports/ the fecundity of BPH on RDV-infected rice plants also had no significant difference with those on RDV-free rice plants ( Fig. 2, P > 0.05). Therefore, rice type and RDV infection had no effects on BPH oviposition preference.
Combined effects of Bt rice and RDV infection on biological parameters of BPH. Survival, development and fecundity of BPH. The survivorship curves of BPH nymphs showed no significant difference among the four treatments ( Fig. 3, χ 2 = 5.734; df = 3; P = 0.125). Rice type, RDV infection, and their interaction (rice type × RDV infection status) had no impact on BPH survival rate (rice type, χ 2 = 1.519; df = 1; P = 0.218; RDV infection, χ 2 = 1.332; df = 1; P = 0.249; rice type × RDV infection, χ 2 = 5.300; df = 3; P = 0.151). Biological parameters including egg hatching rate and egg duration, nymph duration, adult longevity and body weight, as well as female fecundity were listed in Table 1. Rice type, RDV infection status and their interaction had no significant effects on egg hatching rate, adult longevity and body weight, and also the pre-oviposition period, oviposition period and fecundity (total number of eggs laid per female) of female adults (Table 1, P > 0.05). The egg duration was significantly affected by rice type (P < 0.001), and the interaction of rice type and RDV infection (P = 0.001), but not by RDV infection status (P > 0.05). Egg duration of BPH was significantly longer on T1C-19-RDV rice plants than on other rice plants (T1C-19, MH63, MH63-RDV) ( Table 1). Male nymph duration was significantly affected by rice type (P = 0.018), RDV infection status (P = 0.017) and their interaction (P = 0.013). But female nymph duration was only affected by rice type (P = 0.005). Nymph duration on T1C-19 was significantly longer than that on MH63 under RDV-free condition, but not RDV infection status (Table 1). Rice type, RDV infection and their interaction had no significant impacts on the daily fecundity of BPH either ( Life-table parameters of BPH. Life-table parameters of BPH including the intrinsic rate of natural increase (r m ), the net reproductive rate (R 0 ), the mean generation time (T), the doubling time (DT) and the finite capacity of increase (λ) were not affected by rice type or RDV infection status ( Table 2, Supplementary Table S1, P > 0.05).  www.nature.com/scientificreports/

Discussion
The effects of Bt rice on NTAs including herbivores, predators, parasitoids and detritivores have been studied thoroughly in past 20 years 27 . Recently, much attention has been paid to effects of Bt crops on indirect multitrophic interactions, whereas few attentions were paid to the combined effects of Bt rice with other biotic or abiotic factors [13][14][15]18,19,[28][29][30][31] . Interactions between plants and pathogens on insects are complex. In field system, vector and non-vector insects usually shared the same host plants and compete for plant resources 32 . It was reported that RDV-infected rice plants were more attractive to non-viruliferous vector insect GRLHs for feeding 33 . However, no such effect was found on the non-vector insect BPH in the current study. Meanwhile, BPH showed no feeding and oviposition preference between Bt rice T1C-19 and its non-Bt parental rice MH63 under RDV-free or RDV infection status (Fig. 1a), which is similar to the previous study on GRLHs 19 . In general, our results indicated  www.nature.com/scientificreports/ that rice type and RDV infection status had no effect on feeding and oviposition preference of BPH (Figs. 1, 2). Consistently, BPH showed no feeding preference between undamaged Huahui 1 (Cry1Ab/Cry1Ac rice) and its non-Bt parental rice MH63 in both laboratory and field-cage experiments 15 . And rough rice grains of T1C-19 did not change the host preference of the non-target stored-product pest, Rhyzopertha dominica (Coleoptera: Bostrichidae) and its parasitoid wasp Anisopteromalus calandrae (Hymenoptera: Pteromalidae) 34 . The effects of rice type, RDV infection status and their interaction on biological parameters of BPH were also investigated. Although difference in some biological parameters were found between BPH fed on Bt rice T1C-19 and its non-Bt parental rice MH63, the survival, female fecundity and life-table parameters of BPH including r m , R 0 , T, DT and λ showed no significant difference between T1C-19 and MH63 rice plants ( Fig. 3; Tables 1, 2). It indicated that T1C-19 had no effect on population growth of BPH, which was in consistent with previous study 35 . RDV infection status also showed no impacts on BPH survival, development and female fecundity except for male nymph duration, although it significantly changed the phenotypes of rice plants. As a result, life-table parameters showed no significant changes between BPH fed on rice plants with or without RDV infection (Tables 1, 2). No significant effects of interaction between rice type and RDV infection status were found on BPH survival, development and female fecundity except for egg duration and male nymph duration. Based on our results, we concluded that the ecological fitness of BPH was not affected by rice type, RDV infection status and their interaction, as important biological parameters such as egg hatching rate, nymph survival rate and female fecundity showed no significant changes. Similarly, T1C-19 and T2A-1 infected with SRBSDV had no effect on survival, female fecundity and egg hatchability of the non-vector insect BPH 18 . However, there were also reports indicated that biological parameters of non-vector insects would be affected by virus infection. Ecological fitness of non-vector BPH and white-backed rice planthopper (WBPH) Sogatella furcifera (Hemiptera: Delphacidae) were improved on rice plants infected with rice black streak dwarf virus (RBSDV), which is transmitted by small BPH, Laodelphax striatellus (Hemiptera: Delphacidae) 17,36 . The nymph survival rate and female longevity of the non-vector insect BPH were significantly decreased when fed on SRBSDV-infected Y-liangyou1 rice plants 37 . Therefore, the impacts on biological parameters of vector or non-vector insects depend on the interactions between rice plants and viruses. Different viruses, or different rice lines, would result in different impacts.
In summary, our study showed that BPH had no feeding and oviposition preference between Bt rice T1C-19 and its non-Bt parental rice MH63 under RDV-free or RDV infection status. Similarly, Bt rice T1C-19 with or without RDV infection had negligible impact on biological parameters and life-table parameters of BPH at most cases, compared with the non-Bt control rice MH63. Thus, it is clear that Bt rice T1C-19 doesn't affect the ecological fitness of BPH either under RDV stress or not. Moreover, life-table parameters are recommended to be used in the assessment of the effects of Bt crops on NTAs, as it can provide a comprehensive view of the biological characteristics of a given population under defined environmental conditions 38 . It would give us a more intuitive conclusion of the effects on the population growth, even when significant difference exists in some biological parameters.

Materials and methods
Insects. BPH and GRLH colonies were first collected from the paddy fields at the experiment farm of Zhejiang University, Hangzhou, China and then reared continuously on the susceptible rice 'Taichung Native1′ (TN1) seedlings in 80-mesh nylon cage (50 cm length × 50 cm width × 50 cm height) in a climate chamber at 26 ± 1 °C, 75 ± 5% relative humidity and a photoperiod of 14 L:10 D. RDV acquisition of GRLHs was conducted as described by Wang et al. 19 . The 2nd or 3rd instar GRLH nymphs were fed with RDV-infected TN1 rice plants (preserved in our laboratory) after 2 h' starving. These nymphs were transferred to RDV-free TN1 plants 48 h later to go through the latent period. Then the GRLH nymphs were individually released to RDV-free TN1 seedlings in plastic tubes (diameter 2.5 cm, height 25 cm, one seedling in each tube) and labeled. Seedlings were replaced every 2 days. All these seedlings were collected and transplanted into greenhouse, marked with number corresponding to GRLH nymphs individually. Viruliferous GRLHs were selected as RDV symptoms appear on the relative rice seedlings and collected approximately 15 days later 19 . These viruliferous GRLHs were then kept on the RDV-infected TN1 plants under the same conditions as described above until experiment.

Plants.
Bt rice line T1C-19 (Cry1C rice) and its non-Bt parental rice Minghui 63 (MH63) were used in this study 39 . To obtain RDV-infected rice plants, 4th to 5th instar viruliferous GRLH nymphs were inoculated on rice  41 using the jackknife method. Means ± standard error (SE) in the same rows followed by the same letters are not significantly different at P < 0.05 level. b r m : the intrinsic rate of natural increase, R 0 : the net reproductive rate; T: the mean generation time; DT: the doubling time; λ: the finite capacity of increase. Feeding and oviposition preference of BPH. Feeding and oviposition preference experiments were conducted using plastic pots (diameter 20 cm) and foam board with two holes spacing 10 cm apart. One seedling of 45-day-old was transplanted into each hole and cultured in Kimura B nutrient solution. The plots were covered by a transparent polyethylene-plastic cylinder (diameter 18 cm, height 50 cm) with two side-windows of nylon mesh in the middle. The top of cylinder was also covered by nylon mesh for ventilation 19  Statistical analysis. Analysis of data on feeding and oviposition preference, biological parameters, the survivorship curves of BPH nymphs and the daily fecundity of female BPH were performed using SPSS 20.0 software. The impacts of rice type or RDV infection on feeding and oviposition preference of BPH were analyzed by paired-sample t-test. Rice type and RDV infection on the survivorship curves of BPH nymphs were analyzed using the Kaplan-Meier method with Wilcoxon test. Interaction of rice type and RDV infection on survival rate of BPH nymphs were determined using Cox multivariate model. Data on biological parameters including development duration of nymphs, adult longevity and body weight, pre-oviposition period, oviposition period and female fecundity of BPH were analyzed using general linear models (GLMs) followed by Tukey's multiple range test. The daily fecundity of female BPH was analyzed using a mixed-model, repeated-measures analysis of variance (ANOVA) 35 . The life-table parameters, including the intrinsic rate of natural increase (r m ), the net reproductive rate (R 0 ), the mean generation time (T), the doubling time (DT) and the finite capacity of increase (λ) were calculated and analyzed using SAS software (version 9.1) based on the Jackknife program 40,41 .