Transgenic Cry1Ac cotton does not affect the development and fecundity of Chrysoperla carnea

The development and fecundity of the predator Chrysoperla carnea Stephens (Neuroptera: Chrysopidae) were assessed by feeding Aphis gossypii Glover (Hemiptera: Aphididae) that had been reared on transgenic Bacillus thuringiensis (Bt) cotton SGK321 and a non-Bt cotton control (SY321) for two successive generations. We found no significant differences in the developmental stage duration, stage survival, or egg hatch rate between C. carnea fed A. gossypii reared on the Bt and non-Bt cotton. The fecundity per female over a 25-day observation period was very similar between treatments; for C. carnea fed A. gossypii reared on SGK321 vs. SY321, the amount of eggs laid was not significantly different in both generations. Furthermore, a population dynamics of A. gossypii and lacewing (mainly C. carnea) were highly similar in the SGK321 and SY321 treatments during 2016–2017. These results suggest that Bt cotton does not have a significantly negative or positive effect on C. carnea in terms of development, survival, fecundity, or population dynamics.


Introduction
Transgenic cotton producing insecticidal proteins from the bacterium Bacillus thuringiensis (Bt) plays a significant role in insect pest management around the world. Bt cotton not only successfully controls several insect pests of cotton [e.g., Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae)] but also reduces pesticide use in other crops that previously required protection against target pests, thereby boosting crop yields and generating more income for farmers [1][2][3]. In China, Bt cotton has been used to control cotton bollworm since 1997, and it has already been adopted by 95% of farmers in northern China, where it plays a prominent role in pest control in cotton fields [4]. Transgenic Bt cotton has contributed substantially to the reduction of damage by cotton bollworm. Because of its widespread use, efforts to evaluate its safety have received a great deal of attention worldwide, and its potential for harmful effects a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 on non-target species, especially natural predators [5], has been assessed. Widespread planting of Bt cotton has the potential to change the arthropod community in cotton [6], thus affecting populations of sucking pests, such as cotton aphids [7][8], mirid bugs [9], and others.
Cotton aphid Aphis gossypii (Hemiptera: Aphididae) is an important sucking pest of cotton worldwide that can affect cotton plants both by its direct feeding and through its ability to transmit diseases [10]. Predators in cotton fields, such as lacewings, ladybeetles are very important natural enemies in cotton fields, and can control the aphids effectively [11][12].
Lacewings (Chrysopidae) are important predators in cotton fields due to their high predation rates and ecological plasticity. Some widespread species are common in crop fields worldwide, including in China [13]. Chrysoperla carnea Stephens (Neuroptera: Chrysopidae) is an important natural predator that can be reared in the laboratory and used for pest control in the field [14]. This species preys on various pests in natural and cultivated fields, such as aphids, cotton bollworms, whiteflies, and the eggs and young larvae of Lepidoptera [15][16][17][18][19]. Green lacewings have also been used to assess the potential non-target effects of insecticides [20][21], fungicides [22], and GE plants [23][24][25][26][27][28][29] on the third trophic level in agricultural ecosystems. Overall, lacewings play an important role in biological control and integrated pest management in cotton fields.
The Xinjiang Uygur Autonomous Region (XJR) is the largest commercial cotton production region in China. The area of cotton cultivation in the region has expanded in recent years, reaching 2.15 million hectares in 2016 and accounting for more than 80% of all cotton production in China [30]. In Xinjiang, the semi-arid drought-prone climate characteristics are considerably different from that of other areas of China [31]. Therefore, the population dynamics of arthropods and composition of insect species in cotton fields in the XJR differ widely from that in the Yellow River Region or the Changjiang River Region [13].
Many studies have examined the direct and indirect effects of Bt on non-target species, including predators [5,[32][33][34]. Most studies on the effects of Bt toxins on lacewings have focused on the direct effects using artificial diets, and these studies have demonstrated that green lacewing larvae are not sensitive to the toxins of Cry1Ab, Cry1Ac, or Cry2Aa [27,35]. Because reduced prey quality has been reported from tri-trophic studies with other predators, the effects of the Bt toxin on green lacewings may be indirectly mediated by lower prey quality, but only trace amounts of Cry proteins in aphids fed Bt crops can be tested [36]. Because A. gossypii may acquire Cry Bt toxin upon ingestion of cotton plant sap, it may then transmit the toxin to aphid predators. Initial studies of this interaction focused on the transfer of Bt toxin to predators, and lacewing larvae have been commonly examined in tri-trophic experiments through their predation activity [37][38]. The mite Tetranychus urticae (Koch) and Bt-resistant larvae of H. armigera were used as prey and Bt as the toxin, and no adverse effects on green lacewing were observed [24,34,39]. This study assessed the potential for impacts of Bt cotton on the predator C. carnea when feeding on A. gossypii reared on Bt vs. non-Bt cotton for two generations, and we also monitored their population dynamics under field conditions during 2016-2017. The results are available for identifying the ecological effect of Bt cotton on predatory lacewing.

Cotton cultivation and management
Seeds of SGK321 containing the Cry1Ac gene and the non Bt variety ShiYuan321 (SY321) (the parental cultivar of SGK321) were both provided by the Institute of Plant Protection, Chinese Academy of Agricultural Sciences. The study site was at the Scientific Observation and Experimental Station for Crop Pests in Korla, XJR (85˚48 0 N, 41˚44 0 S). Two cotton varieties were planted on 10 April 2016 and 8 April 2017 in plots using drip irrigation under film technology. The planting density was 180000 plants/hm 2 and plot size was 2000 m 2 for each cotton variety, and three plots were established per variety, the space between the two treatment plots was 2 m, and all the plots were randomly arranged random arrangement in the fields. All cotton plots were not treated with any chemical pesticides for this study, other management practices followed normal agronomic practices for the region.

Development and fecundity in the laboratory
Predator rearing. Chrysoperla carnea individuals were collected in a commercial cotton field in 2012, and adults and larvae were fed 20% honey water and eggs of Sitotroga cerealella (Olivier) [40] at the insectarium of the institute of plant protection, Xinjiang Academy of Agricultural Sciences, under constant conditions for adults oviposition (temperature at 23±1˚C, relative humidity at 65%±5%, and L/D cycle of18/6 h).
This experiment began on June 2017 once the A. gossypii populations of the experimental plots were sufficient for the study. Two generations of C. carnea were reared on Aphis gossypii collected from either the SGK321 or SY321 cotton varieties. Experiments with the first and second predator generations began on 10 June and 8 July, 2017, respectively, and under constant conditions (temperature at 26±1˚C, relative humidity at 65%±5%, and L/D cycle of18/6 h).
Lacewing larval development and survival. C. carnea neonate larvae (<8 h) were collected from containers and reared individually in plastic culture dishes (9 cm in diameter, 1.5 cm high). Neonate larvae were fed A. gossypii on cotton leaves picked from either SGK321 or SY321 cotton plants. Predators in their 1 st , 2 nd and 3 rd instar stages were provided with a minimum of 50, 100, and 300 aphids as prey and leaves with A. gossypii were replaced every 24 h. The experiment had three replicates for each treatment, with 30 first instar larvae in each replication. The stage duration of each larval instar and the pupation period was calculated from observations of the numbers of live or dead lacewings until adult emergence, and observed every 8 hours.
Adult fecundity and hatch rate. Twenty newly emerged C. carnea pairs developed on Bt cotton and non Bt cotton were selected respectively, and each pair was reared in a plastic pot (20 cm in diameter, 30 cm high) with honey water (20%) for food.Cauze (3 cm width, 20 cm length) were hung inside cages for female lacewings to lay eggs on. Gauze and a rubber band were used to seal the pots. The gauze was removed daily to record the number of eggs laid by each female and then replaced with fresh gauze. Put a wet cotton balls(d = 2 cm) for keep moist and changed it daily. The fecundity experiment lasted 25 days, and each day the number of hatched eggs and dead adults were recorded. After obtaining new adults of C. carnea, the experiment was repeated with the second generation of lacewings.

Population dynamics in field plots
According to a previous survey, three species of lacewings were identified in local cotton fields, i.e., C. carnea, C. formosa Brauer and C. sinica Tjeder [41][42], and C. carnea was the dominant species. These plots were investigated once every 7 d from 1 June to 17 August in both years, the average plant height were 30cm in 1 June to 100cm in 17August, and average of leaves were sampled. In each plot, a total of 5 points were randomly chosen in which 20 plants for visual sampling. On each plant, the numbers of A. gossypii and lacewings (larva and adults) were counted [43]. Instars and adults of lacewings were recorded respectively, but combined statistics were performed during analysis.

ELISA for Bt cotton
To confirm that SGK321 expressed the Cry1Ac protein and the control variety did not, one leaf was collected from the upper portion of five randomly selected plants from each plot (cotton variety) in 2017. Leaf sampling was repeated every 10 days from 10 June to 8 September, and A. gossypii or C. carnea can be found in both and control cotton fields. Leaves were stored at -20˚C. ELISA testing was used for the detection of the Bt protein. A Cry1Ac kit (EnviroLogix) was used for the ELISA tests, which followed the manufacturer's instructions. Sample concentrations of Bt toxin (absorption values of 450 nm) were measured against a standard curve of protein standard absorption values of 0, 1, 2, and 4 ng/g.

Data analysis
The datasets on the development time, preoviposition period, pupal period, and number of eggs per female were all log-transformed for normality, and the percentage data sets (survival and hatch rate) were arcsine-transformed. The egg number data were using original data. Larval development was analysed using 90 data points (3×30/replicate) for each treatment. Survival data were analysed based on means of three replicates (30 individuals per replicate). The fecundity and hatch rate datasets were analysed using 20 data points for each treatment; and the fecundity data were adjusted according to the number of dead females for analysis because they could not survive for 25 days. The population dynamics of A. gossypii and lacewings between SGK321 and SY321 were analyzed respectively by t-test for each year. All analyses were performed using SPSS 19.0 software [44].

Cry1Ac protein expression for Bt cotton leaves
Cry1Ac protein of fresh weight (FW) expression in the leaves of SGK321 detected by ELISA was 964.2 ng/g on June 10, and it dropped to 802.7 ng/g on June 20 in 2017. The decrease continued until September, with a minimum content of 35.1 ng/g on September 10 (Fig 1). No Cry1Ac protein was detected in the SY321 leaves.

Effect of Cry1Ac cotton aphids on the developmental duration of C. carnea larvae
The durations of the 1 st , 2 nd and 3 rd instars, the total larval stage for the first and second predator generation in the experiment were no significant difference between fed A. gossypii reared on SGK321 (Bt) and SY321(non-Bt) cotton. Similarly, the preoviposition period the pre-pupal stage, and the pupae for the first and second generation between two treatments were not significant ( Table 1).

Effect of Cry1Ac cotton on the survival-survival of C. carnea
The stage survival of C. carnea fed A. gossypii reared on the two cotton varieties was not different between lacewings fed on aphids reared on Bt vs control cotton, and this finding was true for both the 1 st and 2 nd lacewing generations, suggesting that there was no significant effect at this stage. The survival of the 1 st instar and the pre-pupal stage of C. carnea fed A. gossypii reared on SGK321 (Bt) was lower than C. carnea fed A. gossypii reared on SY321 (control) cotton in both the 1 st and 2 nd lacewing generations, whereas the 2 nd instar and pupal stage's survival and the preoviposition period of C. carnea fed A. gossypii reared on SGK321 (Bt) were longer than those of C. carnea fed A. gossypii reared on SY321 (control) cotton in both the 1 st and 2 nd generations. The survival of the total larval stage of C. carnea fed A. gossypii reared on SGK321 (Bt) was higher than that of C. carnea fed A. gossypii reared on the SY321 (control) in the 1 st generation, although the opposite result was observed in the 2 nd instar ( Table 2),however, there was no significant difference between two varieties, and reversed results between two generation may caused by experimental error.

Field dynamics of Aphis gossypii and lacewing (mainly C. carnea)
The field dynamics of Aphis gossypii and lacewings were investigated in SY321 and SGK321 during 2016-2017. The number of Aphis gossypii and lacewings were similar in SY321 and SGK321 cotton fields in each year, and no significant difference have been found between the transgenic Bt cotton and non-transgenic cotton for each pest and predator group in each year (Fig 2 and Table 3).

Discussion
Safety evaluations of Bt cotton must include assessments of both direct (Cry protein-related) and indirect effect. In this study, we found no influence of Bt cotton on the larval lacewing parameters or female fecundity of C. carnea in either of two successive generations. The developmental times of C. carnea life stages appeared normal for all three larval instars (i.e., the prepupal and pupal stages and preoviposition period), which is similar to the results of a previous study [45]. These results suggest that C. carnea was not negatively affected by Bt cotton.
The approach to assessing the potential impact of GE plants on non-target organisms requires the evaluation of the level of insecticidal protein expression of the plant [46][47]. In this study, we measured the protein content of the Bt toxin in Bt cotton leaves and found that Bt cotton expressed the Bt toxin normally, although the levels of Bt proteins in the tender leaves of SGK321 decreased gradually with increasing age of the Bt cotton, and some results showed Cry-protein concentrations strongly decreased with increasing trophic level to values mostly below the detection limit in predators [48]. Therefore, the Bt toxin may not be actually No effect of Bt cotton on predatory lacewing ingested by A. gossypii,and it is unlikely that the aphids contain significant amounts of Cry protein [36,[49][50] This failure of Bt to be present in aphid diets has been observed with different aphid species reared on Bt maize. Although Bt maize expresses higher levels of Bt toxin than cotton, little or no Bt toxin was detected in aphids fed on Bt corn [51][52][53], and a new study showing only traces of some Cry protein contained in aphids fed Bt maize [36,54]. The above findings can be attributed to the fact that aphids cannot ingest considerable amounts of Bt proteins when fed transgenic Bt crops because they feed solely on the phloem of the plants, and phloem has a very low content of Bt toxin [39]. There is evidence that aphid samples from Bt crops containing relative high amounts of Cry proteins might have been caused by contamination for example thrips or thrips faces [36]. Few studies have focused on the arthropod diversity, communities and population dynamics of Bt cotton fields in the XJR [42,[55][56][57], most investigations over two or three years and finding no negative effects because one-year field research may not provide strong evidence in the ecological assessments of transgenic Bt cotton. Diverse ecological factors change greatly over time and can affect the growth of crops and populations of pest and natural enemies [58]. Investigations over two consecutive years were carried out in this study because the population dynamics of A. gossypii and C. carnea changed considerably over this time period but maintained the same trend both on transgenic Bt cotton and non-Bt cotton in the same year, which indicated that different trends in the population dynamics of pests and predators were greatly affected by annual factors but not by the Bt gene. A combination of bioassays in the laboratory indicated that the transgenic Bt cotton did not produce adverse effects on the predator, and this result was similar to that of previous studies [32][33][34][35].
Many studies have shown that Bt toxins have no direct negative effects on the larvae of green lacewings [23,59] and the larvae of C. sinica are not sensitive to Bt [35]. We collected aphids in cotton fields to better simulate the field situations for predation and found results very similar to those of previous research. Moreover, the population dynamics of A. gossypii and lacewings (mainly C. carnea) fed on Bt cotton and non-Bt cotton during 2016-2017 indicated that transgenic Bt cotton did not affect individual development and fecundity or the population dynamics of C. carnea. Table 3. The T-test on the number of lacewings and aphids in Bt and non-Bt cotton respectively. Table 3. The T-test on the number of lacewings and aphids in Bt and non-Bt cotton respectively.

Years
Date