Expression of Bt Protein in Transgenic Bt Cotton Plants and Ecological Fitness of These Plants in Different Habitats

Fitness is one of the key parameters to evaluate the effects of transgenic plants on the ecological environment. To evaluate the ecological risk of transgenic Bt cotton plants growing in different habitats, we determined the expression of the exogenous Bt gene and the fitness of transgenic and non-transgenic cotton plants in three habitats (farmland, grassland, and shrub). We observed that the expression of Bt protein in the farmland was significantly higher than that in the natural habitat, and when the growth environment was suitable, the Bt protein expression level showed a downward trend with the advancement of the growth. There were no significant differences in plant height, aboveground biomass, and seed yield between the Bt transgenic and non-transgenic cotton plants at the same growth stage under the same habitat. Nevertheless, in different habitats, the fitness of the same cotton line showed significant differences. In the farmland habitat, the plant height, aboveground biomass, and seed yield of both transgenic cotton and its non-transgenic isoline were significantly higher than that in the other two natural habitats. The results indicate that the expression of Bt protein does not increase the fitness of the parent plants and would not cause the weeding of the recipient cotton plants.

Introduction

With the rapid development of resistance to chemical pesticides in cotton pests, such as Helicoverpa armigera, the production cost of cotton is increasing and consequently the income of farmers is getting lower (Lu et al., 2012). People are eager to find an effective solution to the problem of pesticide resistance in cotton pests. Bacillus thuringiensis (Bt) is a biological insecticide widely used in pest control in agriculture and forestry (Schnepf et al., 1998). In 1996, Bt-Cry1Ac transgenic insect-resistant cotton (Bt cotton) began to be commercialized in the United States (Bollgard I) and Australia (Ingard I), and has been widely planted in the world’s major cotton-producing countries. In 2018, the planting area of Bt cotton in the world reached 24.9 million hectares (ISSAA, 2019). The commercial planting of transgenic cotton effectively reducing the cost of cotton production, increasing the cotton production, and generating huge economic benefits for the farmers (Naseem and Qaim, 2016). With increasing number of genetically modified (GM) crops entering commercial production, determination of their long-term fitness effects in natural habitats has become an important part of safety assessment and monitoring of such crops (Burke and Rieseberg, 2003; Parimala and Muthuchelian, 2010; Arpaia et al., 2014; Sankaranarayanan and Nalayini, 2015; Lu et al., 2016; Jin et al., 2017; Koller et al., 2019).

On the one hand, plants with exogenous genes can be endowed with resistance to stress to improve their fitness, so that they can survive and reproduce better under adverse environmental conditions. On the other hand, without adverse environmental conditions, the expression of exogenous genes may incur a cost of fitness because of the occupancy of energy and material resources. If the available resources of plants are fixed, the energy consumed for the expression of exogenous genes will inevitably lead to a decrease in resources for growth and reproduction, resulting in the lower growth and reproduction capacity of transgenic plants compared to that of non-transgenic plants (Filipecki and Malepszy, 2006; Xia et al., 2009; Yang et al., 2015). Under conditions of low or no selective pressure, GM plants may show the cost of fitness, which to some extent reduces the possibility of exogenous genes escaping from non-GM plants (Zhang et al., 2012). If the fitness cost is small or there is no fitness cost, then the transferred genes may be preserved and spread in non-transgenic or wild weed populations; this means that there is a risk of exogenous genes escaping to non-transgenic plants (Ellstrand, 2003). Therefore, ecological fitness study is an important parameter to evaluate the ecological effect of transgenic plants. The influence of exogenous genes on the fitness of recipient plants is affected by multiple factors, including the type and construction methods of exogenous genes, the characteristics of transgenic plants, genetic background, the existence of selective pressure, and the environmental conditions under which the transgenic plants grow (Labra et al., 2001; Noro et al., 2007; Mahmoodurrahman et al., 2012). With the large-scale and long-term planting of GM crops, the possibility of these crops entering the natural habitat is greatly increased, and it is increasingly possible that exogenous genes in transgenic plants can enter the wild close relatives of their natural habitats through gene drift (Ellstrand and Hoffman., 1990; Snow et al., 2005; Lu and Yang 2009; Ellstrand et al., 2013). Whether transgenic plants entering the natural habitat and wild related species containing exogenous genes can survive or expand their populations is a scientific question that must be answered in the environmental safety evaluation and monitoring of transgenic organisms.

Transgenic Bt cotton is the predominant commercial transgenic crop in China that has been cultivated for the longest period of time (Lu et al., 2012). However, few studies have been conducted on assessing the fitness of GM cotton in some wild or semi-wild environments. Whether transgenic insect-resistant cotton can survive and proliferate for a long time, whether the insertion of exogenous genes affects the suitability of the recipient cotton for poor ecological environment and then causes exogenous gene escape and invasion, and whether there is a risk of weeding have not been thoroughly studied. We performed experiments to evaluate the ecological fitness effects of transgenic cotton Bt under natural conditions facing different environmental. The results will provide new clues and scientific basis for evaluating the impact of Bt expression on plant fitness. After GM crops escape from their field habitat and enter the natural habitat, the expression of exogenous genes in the natural habitat directly determines the performance and ecological consequences of GM crops in natural habitat. Moreover, the expression of exogenous Bt genes under different environmental conditions was studied to assess the ecological risk of Bt-transgenic plants entering the natural habitat after escape.

Materials and Methods

Cotton Varieties

The cotton varieties used in this study were transgenic Cry1Ab/c cotton Zhong30 and its non-transgenic isolines, Zhong16, which were provided by the Cotton Research Institute of Chinese Academy of Agricultural Sciences. The transgene has one insertion and used CaMV 35S for transgene expression in Zhong30. All cotton plants were transplanted to a field after cultivation in a greenhouse.

Experimental Design

In 2014 and 2015, three experimental sites were set up in Hengxi, Jiangning District, Nanjing, China, which had three different habitats (farmland, grassland, and shrub habitats, respectively) with significant differences in physicochemical properties of the soil (Table 2). The linear distances among the three habitats did not exceed 500 m. Cotton row spacing was 1 m, plant spacing was 1.5 m. Zhong 16 was planted in the East and Zhong 30 in the West. A total of 200 cotton plants are planted on each experimental site. During the survey, samples were randomly divided into groups, for details, see the survey methods for each index below. The composition and biomass of weeds were analyzed, and the physical and chemical properties of soil were analyzed by soil sampling.

Weed Species and Biomass

Five squares (50 cm × 50 cm) were randomly selected at each experimental site. All the plants in the squares were cut at ground level and packed in nylon mesh bags. Under the condition of direct sunlight, the weeds were dried to a constant mass in the ventilated area, and the mass of weeds was recorded. The composition and biomass of weeds are shown in Table 1.

TABLE 1

Table 1 Weed composition and management practices at the different experimental sites.

Physical and Chemical Properties of the Soil

Five samples were taken from each cotton variety in each experimental site, and the five sampling sites formed an S shape. Clean up the sundries 20 mm deep on the surface and cut off a soil sample 5- to 7-cm thick with a spade. The thickness of soil sample at each point was almost the same. After natural air drying, grinding, and sieving, the samples were sent to Nanjing Institute of Soil Research, Chinese Academy of Sciences for determination of the physical and chemical properties according to Agricultural Standard of People’s Republic of China NY/T 1121.7-2014. The specific conditions at each experimental site are shown in Table 2.

TABLE 2

Table 2 Physicochemical properties of the soil at different experimental fields (mean ± stdev).

Investigation of Insects and Diseases

Five investigation sites were set up in the four corners and the middle position of zhong30 and zhong16 respectively, and 10 cotton plants were investigated in each investigation site. Three times of investigation of insects and diseases were conducted by direct observation method at the seedling stage, flowering stage and boll opening stage of cotton.

The species and number of insects on each cotton plant and the degree to which cotton leaves were eaten by insects were investigated and recorded. The damage degree is divided into 6 levels: 0. no feeding; 1. the damage is needle-like; 2. the damage area is less than 1/4; 3. the damage area is less than 1/2; 4. the damage area is more than 1/2; 5. the leaves were completely eaten. Record the incidence of Fusarium and Verticillium wilts of cotton.

Quantification of Cry1Ab/c Protein by Enzyme-Linked Immunosorbent Assay

Two hundred cotton plants were divided into 3 groups at each experimental site and 10 cotton plants were randomly selected from each group for marking for each cotton variety. Cotton was sampled at the bud, boll, and boll opening stages. A complete leaf from the top of each plant was collected. The leaf was immediately placed in a 10-ml sealed centrifuge tube after being detached and stored in a liquid nitrogen tank. After reaching the laboratory, the samples were stored in a deep-freezer at −70°C. The expression of Cry1Ab/c protein in Zhong30 cotton was quantified using a QualiPlate Kit for Cry1Ab/Cry1Ac (EnviroLogix Inc., Portland, ME, USA). About 20 mg of tissue sample was weighed and its exact mass was recorded. The test was done according to the manufacturer’s instructions. The optical density (OD) was measured at 450 nm using a microplate reader (Infinite M2000; Tecan Group Inc., Männedorf, Switzerland). The calibration was done by generating a standard curve of OD against the protein content using the following concentrations of the Bt protein (Cry1Ab) standard (EnviroLogix Inc.): 0.03125, 0.0625, 0.125, 0.25, 0.5 ng, 1.0, 2.0, and 4.0 ng mL⁻¹. The level of Cry1Ab/c protein in the fresh cotton leaf samples was determined using the standard curve and the dilution ratios of the extract (μg g⁻¹ FW).

Measurement of the Fitness Indices of Vegetative and Reproductive Growth

One hundred cotton plants were divided into 3 groups in each cotton variety at each experimental site and 10 cotton plants were randomly selected and marked from each group. The height of 30 randomly marked cotton from soil to terminal tip was measured (accurate to centimeter) using a measuring tape. To determine the aboveground biomass, five unmarked cotton plants were randomly selected from two kinds of cotton at each site. The plant was cut off from a position nearest to the surface, dried under natural shade to a constant mass, and then weighed using a balance (PB602-N, Mettler Toledo).

To detect the effect on the reproductive growth, seeds of marking cotton were harvested. The seeds were dried and cotton was removed from around their short nap with a portable cotton clothes sorting mill (mj57-110d, Xinxing Cotton Machinery Factory, Dong-guang county, Hebei province, China). All the cottonseeds were counted and weighed using an electronic balance (tp-214, Denver Instruments Co. Ltd.), accurate to 0.01 g.

Data Analysis

The expression of the exogenous Cry1Ab/c protein in the Zhong30 cotton plants might be affected by the growth stage and growth environment. Therefore, Tukey’s HSD test with the alpha level at 0.05 was used to determine the difference of exogenous Cry1Ab/c protein expression in different growth stages of cotton in different experimental sites.

In this study, two-factor (factor 1 is cotton variety, factor 2 is experimental site) was used to compare the difference of fitness index between zhong16 and zhong30 cotton. Test factor 1, the relative fitness was calculated by dividing the fitness index of zhong30 cotton by the fitness index of zhong16 cotton, and the relative fitness of each index is added to take the mean value to calculate the comprehensive fitness index. Test factor 2, with farmland ecosystem experimental site I as control, each fitness index of cotton in other experimental sites was divided by corresponding fitness index of experimental site I cotton, and the relative fitness of each index is added to take the mean value to calculate the comprehensive fitness index. Tukey’s HSD test was performed to determine the significant differences of the fitness.

All statistical analyses were performed using SPSS v. 20.0 for Windows (IBM Corp., Armonk, NY, USA).

Results

Insects and Diseases in Cotton Fields

Based on the investigation of cotton diseases and insect pests, the results showed that because the experimental area was located in the non-cotton area and there was no cotton planting in the surrounding kilometers, there were almost no common cotton diseases in the experimental area. Only a small number of neutral insects, such as Araneida, Coccinellidae, and Locusta migratoria manilensis (Meyen) were found at the three experimental sites. Cotton pests such as Helicoverpa armigera, Aphis gossypii Glover and Apolygus lucorμm (Meyer-Dür.) have hardly been found. The target insects of transgenic cotton with Bt gene are Lepidoptera insects, so it can be determined that our experimental sites are in a state of low insect pressure.

Bt Expression in Transgenic Lines Under Different Habitats

In the two-year survey, the trend of change in the content of Bt protein in the leaves of transgenic Bt cotton plants was different at different growth periods at the same experimental site (Table 3). The protein content in the plants at experimental site I and experimental site II decreased gradually with the development of the plant, but the protein content in the plants at experimental site III did not show the same trend. During the three growth stages of the leaves, the content of Bt protein at experimental site I was significantly higher than at the other experimental sites (p< 0.05). The content of Bt protein at experimental site II was higher than that at experimental site III. There were significant differences in the content at the bud stage and flowering and boll-forming stages (p< 0.05), but no significant differences were observed in the content at the boll-opening stage (p> 0.05) (Table 3).

TABLE 3

Table 3 Content of Bt protein in the leaves of Zhong30 plants at different experiment sites during the different growth stages (μg/g fresh mass) (mean ± stdev).

Plant Height of Transgenic and Conventional Cotton Under the Three Habitats

As shown in Figure 1, there were no significant differences in the height between Zhong30 and Zhong16 plants under the same biotope at the same stage in 2014 and 2015 (p>0.05). The height of the same kind of cotton plants was significantly different at the different experimental sites during all the survey periods (p<0.0001) (Figure 1). The height of cotton plants at experimental site I was significantly higher than the height at the other experimental sites with the natural system (p<0.0001). Interestingly, significant difference in heights was observed between the two experimental sites with the natural habitat; the height at experimental site II was significantly higher than that at experimental site III (p<0.0001).

FIGURE 1

Figure 1 Plant height (mean ± SEM) of Zhong30 and Zhong16 cotton plants grown under field (I), grassland (II), and shrub (III) habitat at two stages in 2014 and 2015: (A) Flowering and boll-forming stage in 2014, (B) boll-opening stage in 2014, (C) flowering and boll-forming stage in 2015, (D) boll-opening stage in 2015. Data followed by the same lowercase letters are not significantly different at 0.05 level.

Aboveground Biomass of Transgenic and Conventional Cotton Under Three Habitats

As shown in Figure 2, there were no significant differences in the aboveground biomass between Zhong30 and Zhong16 under the same biotope in 2014 and 2015 (p>0.05). The aboveground biomass of the same cotton line was significantly higher for those grown under the field habitat (experimental site I) than for those grown under the natural conditions (experimental site II and III) by about 4.0 to 10.1-fold in 2014 and 2015. The biomass of the cotton plants grown under the grassland habitat (experimental site II) was 1.7 to 2.0-fold higher than of those grown under the shrub habitat (experimental site III).

FIGURE 2

Figure 2 Aboveground biomass (mean ± SEM) of Zhong30 and Zhong16 cotton plants grown under field (I), grassland (II), and shrub (III) habitat in 2014 (A) and 2015 (B). Data followed by the same lowercase letters are not significantly different at 0.05 level.

Reproductive Growth of Transgenic and Conventional Cotton Under Different Biotopes

In 2014, the seed mass per plant of Bt cotton was less than that of non-Bt cotton, but the difference was not significant at experimental site I (Table 4). Only one Zhong 16 cotton produced seed in experimental site II and experimental site III, which was not enough for statistical analysis. In 2015, the seed mass per plant of Zhong 30 was greater than that of Zhong 16, but the difference was not significant (p=0.0827). There were only two Zhong 30 cotton plants producing bolls at experimental site III, and no seeds were found at experimental site II, and therefore, no statistical analysis could be done.

TABLE 4

Table 4 Production of seeds per plant of different cotton varieties at different experimental sites in 2014 and 2015 (g) (mean ± stdev).

Comparison of the Fitness of Transgenic and Conventional Cotton Plants Under Different Biotopes

In 2014, the total fitness of Zhong30 relative to that of Zhong16 at experimental site II was 1.07, indicating that Zhong30 showed fitness benefits in grassland habitat; the total fitness of Zhong30 relative to that of Zhong16 at experimental site I and experimental site III was 0.89 and 0.97, respectively, indicating that Zhong30 showed a certain fitness cost in field and shrub habitats. In 2015, the total fitness of Zhong30 relative to that of Zhong16 at experimental sites I, II, and III was 1.03, 1.02, and 1.05, respectively (Figure 3). Both in 2014 and 2015, compared with cotton in experimental site I, cotton in experimental site II and III showed significant fitness cost (Figure 4). The results showed that cotton growth environment rather than expression of Bt protein had more significant effect on cotton fitness.

FIGURE 3

Figure 3 The overall fitness of transgenic Bt cotton Zhong30 relative to that of non-transgenic cotton Zhong16 grown under field (I), grassland (II), and shrub (III) habitat in 2014 (A) and 2015 (B). Data followed by the same lowercase letters are not significantly different at 0.05 level.

FIGURE 4

Figure 4 The overall fitness of cotton in grassland (II), and shrub (III) habitat relative to that of cotton in field (I) habitat in 2014 and 2015. Data followed by the same lowercase letters are not significantly different at 0.05 level.

Discussion

There is an important relationship between the fitness of transgenic crops and the plant growth environment (Burke and Rieseberg, 2003). The purpose of this experiment was to simulate the changes of Bt protein expression and fitness of transgenic cotton after escaping to the natural ecosystem, so the selected experimental site and its surroundings had no cotton planting for many years. The cotton plants were disease-free and under the condition of no pest pressure or low pest pressure, which could better simulate the natural ecosystem such as grassland and shrub habitat. Experimental site I belongs to field habitat. The content of soil organic matter is high, and the content of mineral elements needed for the growth of all kinds of plants is closest to the level of ordinary farmland. In the process of cotton growth, weeding, fertilization and irrigation are carried out according to the way of farmland management, but no pesticides are applied. Experimental site II belongs to grassland habitat, the content of organic matter is high, but the content of other mineral elements is low, and the stress of weeds in experimental site II is very serious. Experimental site III belongs to shrub habitat. There are few weeds, but the soil fertility is very low, soil erosion is serious, and cotton is vulnerable to drought stress. After transplanting, cottons in experimental sites II and three are allowed to grow freely without any artificial management.

The properties of soil were reported to affect the expression of Bt protein; drought or saline soil conditions could affect the growth and survival of plants, thereby, affecting the expression of Bt protein (Rochester, 2006). Soil water deficit resulted in the decline of the insecticidal protein expression in bolls of Sikang 3, a hybrid cultivar cotton transgenically expressing a Bt gene (Zhang et al., 2017). Under high and low temperature, drought, and flooding stress habitats, the expression of the exogenous Bt protein in transgenic maize was significantly lower than that under normal habitats (Trtikova et al., 2015). At the same time, the expression of exogenous Bt protein in transgenic rice in habitats stressed with NaCl, flooding (Luo et al., 2008), low, medium, and high salinity (Luo et al., 2017; Fu et al., 2018) showed a decreasing trend compared to the expression in normal habitats. In our study, we compared the content of Bt protein in transgenic cotton plants at different growth stages under different environments. The results for the two years showed that the content of Bt protein in the leaves of Bt transgenic cotton plants in field habitat was significantly higher than that in the other habitats. This may be because under adverse conditions, transgenic plants will produce more substances to resist adverse survival factors (for example, under drought conditions, transgenic plants will produce more drought-resistant and high-temperature resistant proteins), thus affecting the expression of foreign genes. The abovementioned results showed that the external environment exerts a great effect on the expression of the exogenous Bt protein.

Cotton plants were reported to be most resistant to insects in the early stages of their growth as the concentration of Bt toxin declined during the growing season (Sivasupramaniam et al., 2008; Brévault et al., 2013). Cot102 insect-resistant cotton exhibited the best resistance to H. armigera in the early growth stage; with the development of insect-resistant plants, the resistance to phytophagous insects showed a decreasing trend (Llewellyn et al., 2007). The concentration of Bt protein in the terminal leaves declined significantly during the growing season, with the mean toxin concentration being 9.2-times higher for Cry1Ac and 2.9-times higher for Cry2Ab in the young (squaring) cotton than in the old (fruiting) cotton (Carrière et al., 2019). In 2014 and 2015, with the development of Bt transgenic cotton plants, the concentration of Bt protein in field and grassland habitats decreased gradually, while the concentration of Bt protein in shrub habitats remained at a low level with no uniform trend. The soil of cotton in shrub habitat is barren and vulnerable to drought stress, which may be the reason why the Bt protein content of Bt transgenic cotton in shrub habitat is the lowest among the three habitats and does not show significant change in cotton growing season.

The existence of exogenous genes did not change the fitness effect of transgenic cotton plants in 12 natural habitats in Australia, and therefore, the ecological risk of Bt transgenic cotton was inferred to be low (Eastick and Hearnden, 2006). Under low insect pressure, the seed yield of transgenic insect-resistant rice was lower than that of non-transgenic rice, but the difference was not significant (Yang et al., 2011). The growth and reproductive capacity of Cry1Ab/c transgenic rice (Su et al., 2013) was similar to those of their parents under a semi-wild growth environment (high and low insect pressure), which indicated that Cry1Ab/c transgenic rice had lower ecological risk under the semi-wild growth environment. The comparative study on Cry1Ab/c transgenic rice, HH1, and its parent rice, MH63, in Nanchang (having a natural ecosystem) showed that the reproductive and competitive abilities of HH1 and MH63 rice were similar in the natural ecosystem with high insect pressure. (Liu et al., 2015). In this study, in the agroecosystem (field habitat (I)), the growth process of cotton was carefully managed by humans, and there was sufficient water and fertilizer. The growth and reproduction indices (height, aboveground biomass, and seed production per plant) of cotton plants at this site were significantly higher than those for cotton plants of the same kind but grown at other sites. Both the natural ecosystems (grassland (II) and shrub (III) habitats) are subjected to severe environmental stress and almost no seeds were produced in the Bt transgenic and non-transgenic cotton plants. There was no significant difference in the relative fitness between Zhong30 and Zhong16 grown in the same habitat. According to the results of our research, there is no significant difference in the fitness between the Bt transgenic cotton and conventional cotton plants grown under insect-free or low insect pressure conditions in the field ecosystem and in the natural environment under drought and weed stress conditions.

However, the insertion of exogenous genes results in the fitness cost of the host reproductive indicators. In saline-alkaline soils, the insect-resistant transgenic Cry1C* rice T1C-19 showed a strong reproductive capacity, and significantly reduced the loss of yield caused by insects, therefore leading to a higher yield than that of its non-transgenic counterparts MH63 grown (Fu and Liu, 2020). The vegetative and reproductive growth abilities showed a significantly higher fitness cost for the Cry1Ab/c transgenic rice, Huahui1, than that for the parental rice, Minghui63, grown under the same normal farmland and saline-alkaline soil conditions (Fu et al., 2018). Under the same drought conditions, the yield of GM insect-resistant cotton, RCHB 708, was significantly higher than that of non-GM cotton, LRA 5166, and GM cotton showed a fitness advantage (Sankaranarayanan and Nalayini, 2015). Under low insect pressure, Bt/CpT1 transgenic rice showed potential fitness cost in the yield traits (yield per plant, number of grains per plant, number of effective panicles, etc.) compared to that for the parent rice (Xia et al., 2010). The yield (Kg/100 hill) of Bt-SY63 Cry1Ab/c transgenic rice was lower than that of its parent rice, SY63, under low insect pressure with the application of chemical insecticides for two years (Wang et al., 2010). The 1000-grain mass of the insect-resistant rice HH1 was significantly lower than that of its parent rice MH63 under low insect pressure (Xia et al., 2011). In this experiment, the growth of cotton plants was very unsatisfactory due to the environmental stress at the experimental sites II and III. The yield of both the GM and non-GM cotton seeds was low and even zero for some plants.

Conclusion

The expression of exogenous genes in transgenic plants is crucial to the suitability of such plants. We detected the expression of exogenous Bt gene in Bt transgenic cotton plants grown under different environments, and found that the expression levels decreased significantly under the natural environment in the presence of different environmental factors. Therefore, we speculate that under the stress imposed by environmental factors, if the exogenous Bt gene escapes into a natural environment with low insect pressure, the expression behavior of the exogenous Bt gene would be affected by environmental factors, such as drought and weed stress; the fitness benefit of the exogenous Bt gene would not be observed in such cases, and no competitive advantage would be achieved with respect to the species growing in the natural environment. In the present study, there were no significant differences in plant height, aboveground biomass, and seed yield between the Bt cotton and non-Bt cotton plants grown in the natural ecosystem under low insect pressure. Although cotton plants can grow in natural ecosystem, they have weak reproductive capacity and very low seed yield. If the exogenous Bt gene enters the natural ecosystem, it will not increase the fitness of the parental donors, and therefore, the exogenous Bt gene will not cause the weeding of the recipient cotton.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Author Contributions

LL and RG are co-first authors of the article. LL analyzed the data and wrote the manuscript. RG and QQ carried out the experiments, JF helped analyze the data, and BL designed and helped to calibrate the manuscript.

Funding

This work was supported by the “Basic Scientific Research Program in National Nonprofit Scientific Research Institutes (GYZX190103)”; “National Major Special Projects for Genetically Modified Organisms (2016ZX08012-005)”; “Jiangsu Planned Projects for Postdoctoral Research Funds (2018K167C)”.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was funded by grants from the Basic Scientific Research Program in National Nonprofit Scientific Research Institutes (GYZX190103), the National Major Special Projects for Genetically Modified Organisms (2016ZX08012-005), and the Jiangsu Planned Projects for Postdoctoral Research Funds (2018K167C).The authors thank Editage for revising the language and style.

References

Arpaia, S., MesséAn, ,. A., Birch, N. A., Hokannen, H., HäRtel, S., Loon, J. V., Lovei, G., et al. (2014). Assessing and monitoring impacts of genetically modified plants on agro-ecosystems: the approach of AMIGA project. Entomol Exp. Appl. 2, 288–298. doi: 10.4081/entomologia.2014.154

CrossRef Full Text | Google Scholar

Brévault, T., Heuberger, S., Zhang, M., Ellers-Kirk, C., Ni, X., Masson, L., et al. (2013). Potential shortfall of pyramided transgenic cotton for insect resistance managementt. Proc. Natl. Acad. Sci. USA. 110, 5806–5811. doi: 10.1073/pnas.1216719110

CrossRef Full Text | Google Scholar

Burke, J. M., Rieseberg, L. H. (2003). Fitness effects of transgenic disease resistance in sunflowers. Science. 300, 1250–1250. doi: 10.1126/science.1084960

CrossRef Full Text | Google Scholar

Carrière, Y., Degain, B., Unnithan, G. C., Harpold, V. S., Li, X., Tabashnik, B. E. (2019). Seasonal declines in Cry1Ac and Cry2Ab concentration in maturing cotton favor faster evolution of resistance to pyramided Bt cotton in Helicoverpa zea (Lepidoptera: Noctuidae). J. Econ Entomol. 112, 2907–2914. doi: 10.1093/jee/toz236

PubMed Abstract | CrossRef Full Text | Google Scholar

Eastick, R. J., Hearnden, M. N. (2006). Potential for weediness of Bt cotton in northern Australia. Weed Sci. 54, 1142–1151. doi: 10.1614/WS-06-077R.1

CrossRef Full Text | Google Scholar

Ellstrand, N. C., Hoffman, C. A. (1990). Hybridization as an avenue of escape for engineered genes-strategies for risk reduction. BioSciene. 40, 438–442. doi: 10.2307/1311390

CrossRef Full Text | Google Scholar

Ellstrand, N. C., Meirmans, P. G., Rong, J., Bartsch, D., Ghosh, A., De Jong, T. J., et al. (2013). Introgression of crop alleles into wild or weedy populations. Annu. Rev. Ecol. Evol. S. 44, 325–345. doi: 10.1146/annurev-ecolsys-110512-135840

CrossRef Full Text | Google Scholar

Ellstrand, N. C. (2003). Dangerous Liaisons?: When cultivated plants mate with their wild relatives. Plant Sci. 167, 187–188. doi: 10.1016/j.plantsci.2004.02.020

CrossRef Full Text | Google Scholar

Filipecki, M., Malepszy, S. (2006). Unintended consequences of plant transformation: a molecular insight. J. Appl. Genet. 47, 277–286. doi: 10.1007/BF03194637

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, J., Liu, B. (2020). Enhanced yield performance of transgenic Cry1C* rice in saline-alkaline soil. GM Crops Food. 11, 97–112. doi: 10.1080/21645698.2019.1709383.

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, J., Song, X., Liu, B., Shi, Y., Shen, W., Fang, Z., et al. (2018). Fitness cost of transgenic Cry1Ab/c rice under saline-alkaline soil condition. Front. Plant Sci. 9, 1552. doi: 10.3389/fpls.2018.01552

PubMed Abstract | CrossRef Full Text | Google Scholar

ISAAA. (2019). Global status of commercialized biotech/GM crops in 2018. ISAAA Brief No. 54. (Ithaca, NY: ISAAA). doi: 10.13523/j.cb.20190801

CrossRef Full Text | Google Scholar

Jin, X., Chen, Y., Liu, P., Li, C., Cai, X., Rong, J., et al. (2017). Introgression from cultivated rice alters genetic structures of wild relative populations: implications for in situ conservation. Aob Plants. 10, 1–13. doi: 10.1093/aobpla/plx055

CrossRef Full Text | Google Scholar

Koller, T., Brunner, S., Herren, G., Sanchezmartin, J., Hurni, S., Keller, B. (2019). Field grown transgenic Pm3e wheat lines show powdery mildew resistance and no fitness costs associated with high transgene expression. Transgenic Res. 28, 9–20. doi: 10.1007/s11248-018-0099-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Labra, M., Savini, C., Bracale, M., Pelucchi, N., Colombo, L., Bardini, M., et al. (2001). Genomic changes in transgenic rice (Oryza sativa L.) plants produced by infecting calli with Agrobacterium tumefaciens. Plant Cell Rep. 20, 325–330. doi: 10.1007/s002990100329

CrossRef Full Text | Google Scholar

Liu, Y., Ge, F., Liang, Y., Wu, G., Li, J. (2015). Characterization of competitive interactions in the coexistence of Bt-transgenic and conventional rice. BMC Biotechnol. 15, 27–27. doi: 10.1186/s12896-015-0141-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Llewellyn, D. J., Mares, C. L., Fitt, G. P. (2007). Field performance and seasonal changes in the efficacy against Helicoverpa armigera (Hübner) of transgenic cotton expressing the insecticidal protein vip3A. Agr For. Entomol. 9, 93–101. doi: 10.1111/j.1461-9563.2007.00332.x

CrossRef Full Text | Google Scholar

Lu, B., Yang, C. (2009). Gene flow from genetically modified rice to its wild relatives: Assessing potential ecological consequences. Biotechnol. Adv. 27, 1083–1091. doi: 10.1016/j.biotechadv.2009.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Y., Wu, K., Jiang, Y., Guo, Y., Desneux, N. (2012). Widespread adoption of bt cotton and insecticide decrease promotes biocontrol services. Nature 487, 362–365. doi: 10.1038/nature11153

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, B., Yang, X., Ellstrand, N. C. (2016). Fitness correlates of crop transgene flow into weedy populations: A case study of weedy rice in China and other examples. Evol. Appl. 9, 857–870. doi: 10.1111/eva.12377

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Z., Dong, H., Li, W., Ming, Z., Zhu, Y. (2008). Individual and combined effects of salinity and waterlogging on Cry1Ac expression and insecticidal efficacy of Bt cotton. Crop Protect. 27, 1485–1490. doi: 10.1016/j.cropro.2008.06.006

CrossRef Full Text | Google Scholar

Luo, J., Zhang, S., Peng, J., Zhu, X., Lv, L., Wang, C., et al. (2017). Effects of Soil Salinity on the Expression of Bt Toxin (Cry1Ac) and the Control Efficiency of Helicoverpa armigera in Field-Grown Transgenic Bt Cotton. PloS One 12, e170379. doi: 10.1371/journal.pone.0170379

CrossRef Full Text | Google Scholar

Mahmoodurrahman, R. Q. A., Batool, F., Azam, S., Shahid, A. A., Husnain, T. (2012). Transgene copy number and phenotypic variations in transgenic Basmati rice. J. Anim Plant Sci. 22, 1004–1013.

Google Scholar

Naseem, A., Qaim, M. (2016). Genetically modified crops and agricultural development. Am. J. Agr. Econ. 98, 1561–1564. doi: 10.1093/ajae/aaw047

CrossRef Full Text | Google Scholar

Noro, Y., Takanoshimizu, T., Syono, K., Kishima, Y., Sano, Y. (2007). Genetic variations in rice in vitro cultures at the EPSPs-RPS20 region. Theor. Appl. Genet. 114, 705–711. doi: 10.1007/s00122-006-0470-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Parimala, P., Muthuchelian, K. (2010). Physiological response of non-Bt and Bt cotton to short-term drought stress. Photosynthetica. 48, 630–634. doi: 10.1007/s11099-010-0081-9

CrossRef Full Text | Google Scholar

Rochester, I. J. (2006). Effect of Genotype, Edaphic, Environmental Conditions, and Agronomic Practices on Cry1Ac Protein Expression in Transgenic Cotton. J. Cotton Sci. 10, 252–262.

Google Scholar

Sankaranarayanan, K., Nalayini, P. (2015). Performance and behaviour of Bt cotton hybrids under sub-optimal rainfall situation. Arch. Agronom. Soil Sci. 61, 1179–1197. doi: 10.1080/03650340.2014.986112

CrossRef Full Text | Google Scholar

Schnepf, E., Crickmore, N., Van Rie, J., Lereclus, D., Baum, J., Feitelson, J., et al. (1998). Bacillus thuringiensis and its pesticidal crystal proteins. Microbiol. Mol. Biol. Rev. 62, 775–806. doi: 10.1128/MMBR.62.3.775-806.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Sivasupramaniam, S., Moar, W. J., Ruschke, L., Osborn, J. A., Jiang, C., Sebaugh, J. L., et al. (2008). Toxicity and Characterization of Cotton Expressing Bacillus thuringiensis Cry1Ac and Cry2Ab2 Proteins for Control of Lepidopteran Pests. J. Econ. Entomol. 101, 546–554. doi: 10.1603/0022-0493(2008)101[546:TACOCE]2.0.CO;2

PubMed Abstract | CrossRef Full Text | Google Scholar

Snow, A. A., Andow, D. A., Gepts, P., Hallerman, E. M., Power, A., Wolfenbarger, J. M. T. L. (2005). Genetically engineered organisms and the environment: current status and recommendations. Ecoll. Appl. 15, 377–404. doi: 10.1890/04-0539

CrossRef Full Text | Google Scholar

Su, J., Song, Y. N., Yao, Y. X., Chen, J. M., Wu, M. J., Li, G., et al. (2013). Fitness of transgenic insecticidal rice under different growing condition. Mol. Plant Breed. 11, 469–476. doi: 10.3969/mpb.011.000469

CrossRef Full Text | Google Scholar

Trtikova, M., Wikmark, O. G., Zemp, N., Widmer, A., Hilbeck, A. (2015). Transgene Expression and Bt Protein Content in Transgenic Bt Maize (MON810) under Optimal and Stressful Environmental Conditions. PloS One 10, e123011. doi: 10.1371/journal.pone.0123011

CrossRef Full Text | Google Scholar

Wang, Y., Zhang, G., Du, J., Liu, B., Wang, M. (2010). Influence of transgenic hybrid rice expressing a fused gene derived from Cry1Ab and Cry1Ac on primary insect pests and rice yield. Crop Protect. 29, 128–133. doi: 10.1016/j.cropro.2009.10.004

CrossRef Full Text | Google Scholar

Xia, H., Lu, B., Su, J., Chen, R., Rong, J., Song, Z., et al. (2009). Normal expression of insect-resistant transgene in progeny of common wild rice crossed with genetically modified rice: its implication in ecological biosafety assessment. Theor. Appl. Genet. 119, 635–644. doi: 10.1007/s00122-009-1075-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, H., Chen, L., Wang, F., Lu, B. (2010). Yield benefit and underlying cost of insect-resistance transgenic rice: Implication in breeding and deploying transgenic crops. Field Crops Res. 118, 215–220. doi: 10.1016/j.fcr.2010.05.008

CrossRef Full Text | Google Scholar

Xia, H., Lu, B., Xu, K., Wang, W., Yang, X., Yang, C., et al. (2011). Enhanced yield performance of Bt rice under target-insect attacks: implications for field insect management. Transgenic Res. 20, 655–664. doi: 10.1007/s11248-010-9449-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Xia, H., Wang, W., Wang, F., Su, J., Snow, A. A., et al. (2011). Transgenes for insect resistance reduce herbivory and enhance fecundity in advanced generations of crop-weed hybrids of rice. Evol. Appl. 4, 672–684. doi: 10.1111/j.1752-4571.2011.00190.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Cai, X. X., Wang, F., Lu, B. R. (2015). Efficacy of insect-resistance Bt/CpTI transgenes in F₅-F₇ generations of rice crop-weed hybrid progeny: implications for assessing ecological impact of transgene flow. Chin. Sci. Bullet. 60, 1563–1571. doi: 10.1007/s11434-015-0885-x

CrossRef Full Text | Google Scholar

Zhang, F., Lei, S., Liu, Y., Guo, L., Yin, Q., Liu, J., et al. (2012). Differences in ecological fitness between bt transgenic rice and conventional rice under insect-infestation pressures. Chin. J. Appl. Environ. Biol. 18 (1), 35–41. doi: 10.3724/SP.J.1145.2012.00035

CrossRef Full Text | Google Scholar

Zhang, X., Wang, J., Peng, S., Li, Y., Tian, X., Wang, G., et al. (2017). Effects of Soil Water Deficit on Insecticidal Protein Expression in Boll Shells of Transgenic Bt Cotton and the Mechanism. Front. Plant Sci. 8:2107. doi: 10.3389/fpls.2017.02107

PubMed Abstract | CrossRef Full Text | Google Scholar