Next Article in Journal
Design and Test of a Crawler-Type Tiger-Nut Combine Harvester
Next Article in Special Issue
Integrated Crop Management in Sustainable Agriculture
Previous Article in Journal
Has Property Rights Reform of China’s Farmland Water Facilities Improved Farmers’ Irrigation Efficiency?—Evidence from a Typical Reform Pilot in China’s Yunnan Province
Previous Article in Special Issue
Adaptation and High Yield Performance of Honglian Type Hybrid Rice in Pakistan with Desirable Agricultural Traits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Transgenic (Bt) Cotton on the Productivity of Various Cotton-Based Cropping Systems in Pakistan

1
Department of Agronomy, Bahauddin Zakariya University, Multan 60800, Pakistan
2
Physiology/Chemistry Section, Central Cotton Research Institute, Multan 60800, Pakistan
3
College of Agriculture, BZU Bahadur Sub-Campus, Layyah 31200, Pakistan
4
College of Agriculture, University of Layyah, Layyah 31200, Pakistan
5
Department of Plant Protection, Faculty of Agriculture, Harran University, Sanlıurfa 63050, Turkey
6
School of Veterinary and Life Sciences, Murdoch University, 90 South Street, Murdoch, WA 6150, Australia
*
Authors to whom correspondence should be addressed.
Agriculture 2023, 13(2), 276; https://doi.org/10.3390/agriculture13020276
Submission received: 26 December 2022 / Revised: 18 January 2023 / Accepted: 20 January 2023 / Published: 23 January 2023
(This article belongs to the Special Issue Integrated Crop Management in Sustainable Agriculture)

Abstract

:
Cotton (Gossypium hirsutum L.) is an important fiber crop in Pakistan with significant economic importance. Transgenic, insect-resistant cotton (carrying a gene from Bacillus thuringiensis (Bt)) was inducted in the cotton-based cropping systems of Pakistan during 2002, and is now sown in >90% of cotton fields in the country. However, concerns are rising that Bt cotton would decrease the productivity of winter crops (sown after cotton), leading to decreased system productivity. This two-year field study determined the impacts of transgenic (Bt) and non-transgenic (non-Bt) cotton genotypes on the productivities of winter crops (i.e., wheat, Egyptian clover, and canola), and the overall productivities of the cropping systems including these crops. Four cotton genotypes (two Bt and two non-Bt) and three winter crops (i.e., wheat, Egyptian clover, and canola) were included in the study. Nutrient availability was assessed after the harvest of cotton and winter crops. Similarly, the yield-related traits of cotton and winter crops were recorded at their harvest. The productivities of the winter crops were converted to net economic returns, and the overall economic returns of the cropping systems with winter crops were computed. The results revealed that Bt and non-Bt cotton genotypes significantly (p < 0.05) altered nutrient availability (N, P, K, B, Zn, and Fe). However, the yield-related attributes of winter crops were not affected by cotton genotypes, whereas the overall profitability of the cropping systems varied among the cotton genotypes. Economic analyses indicated that the Bt cotton–wheat cropping system was the most profitable, with a benefit–cost ratio of 1.55 in the semi-arid region of Pakistan. It is concluded that Bt cotton could be successfully inducted into the existing cropping systems of Pakistan without any decrease to the overall productivity of the cropping system.

1. Introduction

Cotton (Gossypium hirsutum L.) is the most important fiber crop grown in Pakistan. It serves as a foundation for Pakistan’s economy [1]. Pest infestations exert negative effects on cotton productivity, as >160 insect pests infest cotton at various growth stages [2,3]. Subsequently, cotton farmers incur significant amounts of money on pesticides to combat the pest infestation [2,4,5]. Farmers in Pakistan spend ~$300 million every year on pest control, and most of the pesticides (~80%) are sprayed in the cotton crop [6,7]. The extensive use of pesticides causes significant environmental and human health hazards [8,9]. The cultivation of transgenic, insect resistance crops can decrease the use of pesticides [10].
The use of genetically modified (GM) crops is drawing significant attention in Pakistan, due to their insect resistance abilities [11]. Transgenic cotton (commonly known as Bt cotton) is one of the greatest examples of GM crops, which are bollworms-resistant [12]. Bt cotton can reduce the damage caused by specific insects, and enhance crop productivity because of its resistance to bollworms [13]. The reduced use of insecticides and lower pest infestation are two major benefits of GM crops [14,15,16]. Numerous studies have shown that GM crops with Bt Cry proteins are resistant to various lepidopterans [17,18]. Bt cotton could help farmers in reducing the use of pesticides [2,5,19]. Therefore, Bt cotton is considered to be environmentally friendly [20,21].
There are increasing concerns on the potential negative consequences of GM crops on the soil health and productivity of other crops [22,23,24]. Numerous studies have revealed that the cultivation of GM crops increases the levels of Bt toxins in the soil [25,26,27]. Bt toxins are naturally present in the soil, and the recurrent planting of GM crops raises their concentration in the soil and alters the composition and behavior of soil microorganisms [28,29]. Bt toxins produced by GM crops enter the soil and could exert negative impacts on the productivity of other crops [30]. The cultivation of the Bt cotton line ‘GK19’ increased the accumulation of Bt proteins in the soil under salinity stress [31]. Increased levels of Bt toxins might exert negative impacts on agroecosystems [29] and alter the chemical composition of the root zone [27,32]. Similarly, toxins–microbe interactions significantly influence soil properties and nutrient availability. The changing rhizosphere conditions due to the cultivation of GM crops increased soil phosphorus (P) availability [33]. Furthermore, the cultivation of GM crops increases the available N and the oxidative metabolism in soil because of the enhanced activities of urease and dehydrogenase enzymes [34].
The induction of Bt cotton in the cotton-based cropping systems of Pakistan is being criticized due to the associated negative impacts. However, its introduction had favorable consequences on the environment and on farmers’ health [35]. Cases of pesticide poisoning in farmers have decreased in China, India, and Pakistan after the introduction of Bt cotton [36,37]. However, the impacts of Bt cotton cultivation on the productivity of subsequent winter crops and the overall productivity of cropping systems have been less explored in Pakistan. Cotton is followed by wheat crop in the cotton-based cropping systems of Pakistan [38]. However, the recurrent cultivation of the same crops caused significant weed and insect infestation problems. It is suggested to include alternative crops (other than wheat) in the cotton-based cropping system of the country [38]. However, the impacts of Bt cotton on the productivity of alternative crops are still unknown.
This study assessed the impacts of Bt cotton cultivation on nutrient availability, the yield-related traits of different crops (wheat, Egyptian clover, and canola) sown after cotton, and the overall productivity of the cropping systems, including these crops. It was hypothesized that the cultivation of Bt cotton would reduce the yield of winter crops; however, the overall productivity of the cropping systems would not be affected. It was further hypothesized that cotton-based cropping system with Bt genotypes would have higher economic returns compared to those having non-Bt genotypes.

2. Materials and Methods

2.1. Experimental Site

This field study was conducted at CCRI (Central Cotton Research Institute), Multan (30.2° N, 71.43° E, and 122 meters above sea level), Pakistan, during 2016–2017 and 2017–2018. The soil of the experimental site was analyzed to assess the nutrient availability and the physico-chemical characters before and after the experiments. The soil texture was silt–clay–loam. The soil physico-chemical properties are given in Table 1.

2.2. Experimental Treatments

This experiment consisted of two factors, i.e., cotton genotypes (Bt and non-Bt) and winter crops. The Bt genotypes included in the study were ‘CIM-616’ (Bt1) and ‘GH Mubarik’ (Bt2), while the non-Bt genotypes were ‘CIM-620’ (non-Bt1) and ‘N-414’ (non-Bt2). Similarly, wheat (Triticum aestivum L.), Egyptian clover (Trifolium alexandrinum L.), and canola (Brassica napus) were winter crops that were sown after cotton harvest. The cultivation of winter crops resulted in three possible cropping systems, i.e., cotton–wheat, cotton–Egyptian clover, and cotton–canola. These systems may further be classified into to Bt and non-Bt cotton-based cropping systems. The Bt and non-Bt cotton genotypes were sown in May, whereas winter crops were sown in November following the harvest of the cotton crop (Table 2). The experiment was laid out according to two factorial design, where cotton genotypes were kept in main plots (6 m × 10 m), whereas winter crops were randomized in the sub-plots (2 m × 10 m). All treatments had three replications, and the experiment was repeated for two years. The sub-plots were regarded as being experimental units, and each experimental unit had three replications, as described above.

2.3. Crop Husbandry

Cotton and winter crops were planted following the recommendations in production technology provided by the local agriculture extension department (https://www.agripunjab.gov.pk/, accessed on 12 January 2015). The recommendations followed in the study are given in Table 2. Irrigation was applied according to the moisture needs of the crops. All crop protection measures were taken to protect crops from insect and disease infestation. All crops were harvested when they reached physiological maturity.

2.4. Data Collection

2.4.1. Soil Properties

Particle size distribution was determined using the hydrometer method [39]. Soil EC was recorded using a digital EC meter following the standard procedures detailed by Dellavalle [40]. A digital pH meter was used to measure the soil pH from a saturated soil paste [40].

2.4.2. Nutrient Availability

Soil nitrogen (N) availability was measured spectrophotometrically with a segmented-flow system. The phosphorus (P) was determined using the vanadomolybdate method, potassium (K) through flame photometry, and zinc (Zn) and iron (Fe) through atomic absorption spectrophotometry [41]. Soil organic matter was measured using a loss-on-ignition protocol, as introduced by Hoogsteen et al. [42].

2.4.3. Weed Infestation

Weed infestation was evaluated 45 days after the sowing of cotton and winter crops. The weeds present in a l m2 quadrat were counted from each experimental unit at three different places. The weeds were identified, grouped into narrow and broad-leaved, and their densities were computed. The density of each experimental unit was averaged from different locations within a replication.

2.5. Morphological and Yield-Related Traits

2.5.1. Cotton

The number of sympodial and monopodial branches were counted from 10 randomly selected plants in each experimental unit and averaged. The weights of 10 opened bolls from a single plant were measured using a sensitive balance from 10 randomly selected plants in each experimental unit, and averaged. Three manual pickings were performed from each experimental unit to record the seed cotton yield. The seed cotton yields of three pickings were added and converted to seed cotton yield per hectare, using a unitary method. The cotton stalks were harvested and left in the field for two weeks. Afterwards, the stalks were weighed to record the total biomass (biological yield), and they were expressed in kg ha−1. The harvest index was estimated by dividing the seed cotton yield to biological yield, and this was expressed as a percentage.

2.5.2. Wheat

The number of productive (spike-bearing) tillers present in a 1 m2 area were counted. The lengths of 10 randomly selected spikes were recorded from four central rows in each treatment and averaged. The number of grains were counted from 10 randomly selected spikes and averaged. The weight of 1000 grains from three random samples in each treatment was measured using a sensitive balance. The grain yield per plot was measured on a sensitive balance when the seed reached the required moisture level, i.e., 12%, and converted to t ha−1.

2.5.3. Canola

The number of siliques per plant were counted from three randomly selected plants in each experimental unit. Ten random siliques were opened, and the number of grains in them were counted and averaged. The weight (g) of 1000 seeds from randomly sampled seeds per plot was measured on a sensitive balance. The mature crop was harvested, sundried, and threshed manually to record the seed yield, which was converted into kg ha−1. The biological yield was recorded by weighing the total above ground biomass harvested from four central rows of each experimental plot, and this was converted into kg ha−1.

2.5.4. Egyptian Clover

All plants within the experimental unit were harvested during each cut, and weighed to record the fresh forage yield. A pre-weighed amount of fresh forage was oven-dried, and fresh forage yield was converted into dry forage yield using a unitary method. The fresh and dry forage yields were converted into t/ha. Crude protein was measured via a near-infrared spectroscopy system [43].

2.6. Economic Analysis

The profit abilities of different cotton genotypes via winter crops interactions were computed following CIMMYT [44]. The input costs and outputs obtained in monetary terms were calculated. The costs regarding land rent, irrigation, labor, seeds, fertilizers, pesticides, sowing, harvesting, etc., were computed. The existing market prices of the produce were used to compute the gross income. These expenses were deducted from the gross income to obtain the net income. The benefit–cost ratio (BCR) was computed by dividing the net economic returns with the expenses incurred.

2.7. Statistical Analysis

The collected data for the nutrient availability, weed density, and yield-related parameters of different crops were checked for normality using the Shapiro-Wilk normality test [45]. The parameters having non-normal distributions were normalized using the Arcsine transformation technique to meet the normality assumption of the Analysis of Variance (ANOVA). The differences among the years were tested, which were significant; therefore, the data of both years were analyzed, presented, and interpreted separately. A two-way ANOVA was used to test the significance among the treatments, and the means were compared using a least significant difference post hoc test at 95% probability, where ANOVA denoted significant differences [46]. The interactive effects of cotton genotypes and winter crops were significant for most of the studied traits. Therefore, the interactive effects were presented and interpreted. A one-way ANOVA was used to analyze the data on the yield-related traits of cotton. All statistical computations were performed on SPSS statistical software, version 21.0 [47].

3. Results

3.1. Nutrient Availability

Soil nutrient availability and organic matter content were significantly altered by cotton genotypes via winter crops interaction during both years (Table 3). Wheat sown after the Bt cotton genotype ‘GH-Mubarik’ had the highest available N during each year, while canola sown after the non-Bt genotypes ‘N-414’ and ‘CIM-620’ resulted in the lowest available N during both years (Table 3). Egyptian clover sown after the non-Bt cotton genotype ‘CIM-620’ during the first year, and the Bt genotype ‘CIM-616’ during the second year, recorded the highest P, whereas the lowest values were recorded for wheat sown after the non-Bt genotype ‘N-414’ during both years (Table 3). The Egyptian clover sown after the non-Bt genotype ‘N-414’ resulted in the highest available K, which was statistically similar to the wheat cultivation after the Bt genotype ‘GH-Mubarik’. The lowest available K was noted for canola sown after the Bt cotton genotypes during the first year, and the non-Bt genotype ‘N-414’ during the second year (Table 3). Egyptian clover following the non-Bt genotypes resulted in the highest available Zn during both years, while the lowest Zn was recorded for canola sown after the Bt genotype ‘GH-Mubarik’ during both years (Table 3). The interactive effects of wheat and the non-Bt genotype ‘CIM-620’ resulted in the highest available Fe during each year, while the lowest Fe was recorded for Egyptian clover sown after the non-Bt genotype ‘CIM-620’ during first year, and canola sown after the non-Bt genotype ‘N-414’ during the second year (Table 3). The highest organic matter content was recorded for wheat sown after the Bt genotypes during both years, while canola sown after the non-Bt genotypes resulted in the lowest value of soil organic matter during both years (Table 3).

3.2. Weed Density

The densities of broadleaved, narrow-leaved, and total weeds were significantly affected by the interactive effect of cotton genotypes and winter crops. Overall, Bt genotypes recorded lesser weed infestation compared with the non-Bt genotypes included in the current study (Figure 1).
Similarly, the wheat crop had the highest density of narrow-leaved, broadleaved, and total weeds during both years (Figure 2). Wheat sown after the non-Bt genotype ‘CIM-620’ had the highest density of broad-leaved weeds, whereas the lowest density was observed in Egyptian clover and canola crops sown after both Bt genotypes (Table 4).
Similarly, the lowest density of narrow-leaved weeds was recorded for Egyptian clover sown after the Bt genotype ‘GH-Mubarik’ and non-Bt genotype ‘CIM-620’ during both years, whereas the highest density was noted for wheat sown after the Bt genotype ‘CIM-616’ and the non-Bt genotype ‘CIM-620’ during both years (Table 4). Likewise, the highest density of total weeds was recorded in wheat crop sown after the non-Bt genotype ‘CIM-620’, whereas the lowest density was recorded in Egyptian clover sown after the Bt genotypes during both years (Table 4).

3.3. Yield-Related Attributes of Cotton

Various yield-related traits of cotton were significantly affected by genotypes, except for the number of monopodial branches and the harvest index (Table 5).
The Bt genotype ‘CIM-616’ and non-Bt genotype ‘CIM-620’ recorded the highest number of sympodial branches during the first year, whereas both Bt genotypes recorded the highest number of sympodial branches during the second year. The non-Bt genotypes had the lowest number of sympodial branches during the second year (Table 5). The Bt genotype ‘CIM-616’ recorded the highest seed cotton yield during the first year, while it was not affected by genotypes during the second year (Table 5).

3.4. Yield-Related Attributes of Wheat

Yield-related traits of wheat crop were significantly impacted by different cotton genotypes with some exceptions, i.e., grain yield during the first year, and harvest index during both years. Overall, the wheat crop planted after the non-Bt genotypes recorded higher values for the number of productive tillers, number of grains per spike, 1000-grain weight, and grain and biological yields (Table 6). Wheat sown after the Bt genotypes recorded lower values for the number of productive tillers, number of grains per spike, 1000-grain weight, and grain, and biological yields, compared to non-Bt genotypes during both years (Table 6).

3.5. Yield-Related Attributes of Canola

Yield-related traits of canola were not affected by different cotton genotypes except for seed yield during the first year, where the crop sown after non-Bt genotypes recorded higher values compared to the Bt genotypes (Table 7).

3.6. Yield-Related Attributes of Egyptian Clover

The yield-related traits of Egyptian clover were significantly affected by different cotton genotypes during both years (Table 8). The crop sown after non-Bt genotypes recorded higher values for the total forage yield, dry matter content, and crude protein during both years of the study, compared to the Bt genotypes (Table 8).

3.7. Economic Returns/System Productivity

Economic analysis showed that the Bt genotypes–wheat cropping system resulted in the highest net benefits, whereas the non-Bt genotypes–canola cropping system recorded the lowest net benefits (Table 9). Similarly, the Bt genotypes–wheat cropping system resulted in the highest benefit–cost ratio (BCR), while the non-Bt genotypes–canola cropping system resulted in the lowest BCR (Table 9).

4. Discussion

The results of the current study indicated that nutrient availability, yield-related attributes of winter crops, and overall system productivity were significantly affected by the cotton genotypes. Bt cotton genotypes improved the nutrient availability and system productivity, which are directly linked to the higher fertilizer input and the better yield of the Bt genotypes. Increased fertilizer use after the introduction of Bt genotypes have increased crop yields because of better pest control [48]. Furthermore, farmers started applying more fertilizers after the induction of Bt genotypes in the existing cropping systems of Pakistan [49]. Cotton yields and related benefits may also be affected by numerous factors, such as changes to irrigation systems, crop production technologies, agronomic practices, farmer training, or weather fluctuations, etc. [48]. Earlier studies have reported that nutrient availability may vary across Bt and non-Bt genotypes because of the differences in their nutrient requirements and absorption [33]. Therefore, the results of the current study regarding nutrient availability are in agreement with the earlier studies.
The highest P, N, Zn, Fe, and organic matter contents were recorded from the soil cultivated with winter crops after the Bt genotypes. Wheat crop sown after Bt genotypes resulted in higher values of available N, Fe, and organic matter contents, whereas Egyptian clover following Bt genotypes resulted in higher P, K, and Zn contents. Nutrient uptake is dependent on plants and their genetic makeup. Several factors affect the nutrient uptake capacity of plants [50]. These factors include root surface area, and the type and quantity of root exudates released in the rhizosphere and microbial communities [50]. Moreover, the plant characteristics and interactive effect between roots and soil microorganisms also play significant roles in nutrient uptake [51]. The quantity of nutrients available to plants depends mainly upon their availability in the root zone [52]. Genetically modified (GM) crops can disrupt soil nutrient cycles due to changes in the root zone [53]. The quantity and quality of root exudates affects microbial activity, which alters the solubility of mineral or fixed P, and P availability [2,21,27]. The availability of soil nutrients is significantly affected by the cultivation of Bt cotton. Growing Bt cotton also decreased the available N and K, while increasing Zn and P [29,54]. However, our study indicated that the available N was increased in the soil cultivated with Bt genotypes. The Bt genotypes received higher amounts of nutrients, which can be linked to the increased nutrient availability. Similarly, the availabilities of K, Zn, and P also increased in the treatments with Bt genotypes. Higher root biomass-mediated exudation is responsible for the enhanced availability of Zn and Fe in Bt cotton cultivated soil, compared to non-Bt cotton [55,56].
Weed infestation exerts significant negative impacts on crop yields [57,58]. Weed infestation decreases cotton yield by 0.26 to 66%, depending on the weed species and their densities [59]. The cultural practices used in the existing cropping systems of Pakistan encourage the growth of several weed species [60]. The recurrent cultivation of Bt genotypes may result in the proliferation of specific weed species. Different genotypes significantly vary in their weed competitive ability [61]. Several earlier studies have reported that cotton genotypes significantly differ in their competitive ability with weeds [62,63]. The weed competitive abilities of these cultivars were linked with their potential to establish a crop canopy. The cultivars which developed a dense canopy in a shorter period were more competitive. However, these studies did not include any Bt genotypes. Low weed density was recorded in the soil cultivated with Bt genotypes in the current study. Weed competition of Bt genotypes in the current study could be linked with its quicker canopy development, compared to non-Bt cotton genotypes.
Wheat and Egyptian clover sown after non-Bt genotypes had better yield-related traits. However, the yield-related traits of canola were not affected by cotton genotypes. The lower yields of winter crops in the fields cultivated with Bt genotypes can be linked with the increased levels of Bt toxins in the soil and the higher nutrient consumption by cotton plants. Moreover, the improvement in the yield-related traits of winter crops in the fields cultivated with non-Bt genotypes can be linked to the low Bt toxin levels in these soils. It has been observed that toxins produced in the aerial parts and roots of Bt cotton may cause soil pollution upon their release [34,64,65]. Bt toxins released from the plants become absorbed or bound to the soil particles, and then they become safe from degradation from other microorganisms that are present in the soil [66]. The recurrent cultivation of Bt cotton in the same field increases the level of Bt toxins in the soil, which can change the activity and composition of the soil microbes and the soil biochemical nature [29,67,68,69].
The Bt cotton–wheat cropping system revealed the highest net income and benefit–cost ratio (BCR). The highest productivity of this system was due to a higher production of Bt cotton and wheat. The Bt genotypes produced the highest yield due to a lower rate of insect infestation, compared with non-Bt genotypes. However, wheat yield was higher in the fields cultivated with non-Bt genotypes, and a lesser number of sprays in Bt cotton for pest management decreased the input costs. This eventually reduced expenses, which resulted in a higher net income and BCR than with non-Bt cotton.
The results confirmed the hypothesis of the study, where Bt genotypes exerted negative impacts on the yield-related traits of winter crops (with some exceptions) and improved the overall system productivity. Therefore, Bt genotypes can be included in the cotton-based cropping systems without any decrease in the productivity and economic returns.

5. Conclusions

The results of the current study indicated that nutrient availability, weed infestation, the yield-related traits of winter crops, and the system productivity of various cotton-based cropping systems were significantly affected by Bt and non-Bt cotton genotypes. Overall, Bt genotypes had higher yields than non-Bt genotypes. The soil cultivated with Bt genotypes resulted in higher N, P, and Zn availabilities. The yield-related traits of winter crops were negatively affected by the Bt genotypes. Economic analysis indicated that Bt cotton, followed by wheat, resulted in the highest economic returns and benefit–cost ratios. Therefore, Bt cotton can be successfully inducted in the cotton–wheat cropping systems of semi-arid regions in Pakistan in order to obtain higher economic benefits.

Author Contributions

Conceptualization, S.U.-A., S.F. and M.H.; data curation, M.W.R.M. and F.A.; formal analysis, M.W.R.M.; investigation, M.W.R.M. and A.-u.-R.; methodology, S.U.-A., A.-u.-R., S.F. and M.H.; project administration, M.H.; resources, F.A. and M.H.; software, S.U.-A. and A.-u.-R.; supervision, M.H.; validation, S.F.; visualization, S.U.-A.; writing—original draft, M.W.R.M.; writing—review and editing, F.A., S.U.-A., S.F. and M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

All data are within the manuscript.

Acknowledgments

This manuscript has been prepared from the Ph.D. thesis of the first author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. GOP Economic Survey of Pakistan; Economic Advisory Wing: Islamabad, Pakistan, 2021.
  2. Kouser, S.; Spielman, D.J.; Qaim, M. Transgenic Cotton and Farmers’ Health in Pakistan. PLoS ONE 2019, 14, e0222617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Naeem-Ullah, U.; Ramzan, M.; Bokhari, S.H.M.; Saleem, A.; Qayyum, M.A.; Iqbal, N.; Habib ur Rahman, M.; Fahad, S.; Saeed, S. Insect Pests of Cotton Crop and Management Under Climate Change Scenarios. In Environment, Climate, Plant and Vegetation Growth; Springer International Publishing: Cham, Switzerland, 2020; pp. 367–396. [Google Scholar]
  4. Zehr, U.B. Cotton: Biotechnological Advances; Springer Science & Business Media: Cham, Switzerland, 2010; Volume 65, ISBN 3642047963. [Google Scholar]
  5. Kouser, S.; Qaim, M. Valuing Financial, Health, and Environmental Benefits of Bt Cotton in Pakistan. Agric. Econ. 2013, 44, 323–335. [Google Scholar] [CrossRef]
  6. Arshad, M.; Suhail, A.; Asghar, M.; Tayyib, M.; Hafeez, F. Factors Influencing the Adoption of Bt Cotton in the Punjab, Pakistan. J. Agric. Soc. Sci. 2007, 11, 19. [Google Scholar]
  7. Khan, M.; Mahmood, H.Z.; Damalas, C.A. Pesticide Use and Risk Perceptions among Farmers in the Cotton Belt of Punjab, Pakistan. Crop Prot. 2015, 67, 184–190. [Google Scholar] [CrossRef]
  8. Tooker, J.F.; Pearsons, K.A. Newer Characters, Same Story: Neonicotinoid Insecticides Disrupt Food Webs through Direct and Indirect Effects. Curr. Opin. Insect Sci. 2021, 46, 50–56. [Google Scholar] [CrossRef]
  9. Ziółkowska, E.; Topping, C.J.; Bednarska, A.J.; Laskowski, R. Supporting Non-Target Arthropods in Agroecosystems: Modelling Effects of Insecticides and Landscape Structure on Carabids in Agricultural Landscapes. Sci. Total Environ. 2021, 774, 145746. [Google Scholar] [CrossRef]
  10. Peshin, R.; Hansra, B.S.; Singh, K.; Nanda, R.; Sharma, R.; Yangsdon, S.; Kumar, R. Long-Term Impact of Bt Cotton: An Empirical Evidence from North India. J. Clean. Prod. 2021, 312, 127575. [Google Scholar] [CrossRef]
  11. Rana, M.A. When Seed Becomes Capital: Commercialization of Bt Cotton in Pakistan. J. Agrar. Chang. 2021, 21, 702–719. [Google Scholar] [CrossRef]
  12. Cheema, H.M.N.; Khan, A.A.; Noor, K. Bt Cotton in Pakistan. In Genetically Modified Crops in Asia Pacific; CSIRO Publishing: Clayton, Australia, 2021; p. 91. [Google Scholar]
  13. Lv, N.; Liu, Y.; Guo, T.; Liang, P.; Li, R.; Liang, P.; Gao, X. The Influence of Bt Cotton Cultivation on the Structure and Functions of the Soil Bacterial Community by Soil Metagenomics. Ecotoxicol. Environ. Saf. 2022, 236, 113452. [Google Scholar] [CrossRef]
  14. Smyth, S.J.; Kerr, W.A.; Phillips, P.W.B. Global Economic, Environmental and Health Benefits from GM Crop Adoption. Glob. Food Sec. 2015, 7, 24–29. [Google Scholar] [CrossRef]
  15. Hutchison, W.D.; Burkness, E.C.; Mitchell, P.D.; Moon, R.D.; Leslie, T.W.; Fleischer, S.J.; Abrahamson, M.; Hamilton, K.L.; Steffey, K.L.; Gray, M.E.; et al. Areawide Suppression of European Corn Borer with Bt Maize Reaps Savings to Non-Bt Maize Growers. Science 2010, 330, 222–225. [Google Scholar] [CrossRef] [Green Version]
  16. Kumar, K.; Gambhir, G.; Dass, A.; Tripathi, A.K.; Singh, A.; Jha, A.K.; Yadava, P.; Choudhary, M.; Rakshit, S. Genetically Modified Crops: Current Status and Future Prospects. Planta 2020, 251, 91. [Google Scholar] [CrossRef]
  17. Di Lelio, I.; Barra, E.; Coppola, M.; Corrado, G.; Rao, R.; Caccia, S. Transgenic Plants Expressing Immunosuppressive DsRNA Improve Entomopathogen Efficacy against Spodoptera littoralis Larvae. J. Pest Sci. 2022, 95, 1413–1428. [Google Scholar] [CrossRef]
  18. Katta, S.; Talakayala, A.; Reddy, M.K.; Addepally, U.; Garladinne, M. Development of Transgenic Cotton (Narasimha) Using Triple Gene Cry2Ab-Cry1F-Cry1Ac Construct Conferring Resistance to Lepidopteran Pest. J. Biosci. 2020, 45, 31. [Google Scholar] [CrossRef]
  19. Klümper, W.; Qaim, M. A Meta-Analysis of the Impacts of Genetically Modified Crops. PLoS ONE 2014, 9, e111629. [Google Scholar] [CrossRef] [Green Version]
  20. Kumar, S.; Chandra, A.; Pandey, K.C. Bacillus thuringiensis (Bt) Transgenic Crop: An Environment Friendly Insect-Pest Management Strategy. J. Environ. Biol. 2008, 29, 641–653. [Google Scholar]
  21. Tokel, D.; Genc, B.N.; Ozyigit, I.I. Economic Impacts of Bt (Bacillus thuringiensis) Cotton. J. Nat. Fibers 2022, 19, 4622–4639. [Google Scholar] [CrossRef]
  22. Halford, N.G.; Shewry, P.R. Genetically Modified Crops: Methodology, Benefits, Regulation and Public Concerns. Br. Med. Bull. 2000, 56, 62–73. [Google Scholar] [CrossRef] [Green Version]
  23. Kuzma, J.; Grieger, K. Community-Led Governance for Gene-Edited Crops. Science 2020, 370, 916–918. [Google Scholar] [CrossRef]
  24. Sendhil, R.; Nyika, J.; Yadav, S.; Mackolil, J.; Prashat, G.R.; Workie, E.; Ragupathy, R.; Ramasundaram, P. Genetically Modified Foods: Bibliometric Analysis on Consumer Perception and Preference. GM Crops Food 2022, 13, 65–85. [Google Scholar] [CrossRef]
  25. Liu, J.; Liang, Y.; Hu, T.; Zeng, H.; Gao, R.; Wang, L.; Xiao, Y. Environmental Fate of Bt Proteins in Soil: Transport, Adsorption/Desorption and Degradation. Ecotoxicol. Environ. Saf. 2021, 226, 112805. [Google Scholar] [CrossRef] [PubMed]
  26. Zhou, X.; Li, H.; Liu, D.; Hao, J.; Liu, H.; Lu, X. Effects of Toxin from Bacillus thuringiensis (Bt) on Sorption of Pb (II) in Red and Black Soils: Equilibrium and Kinetics Aspects. J. Hazard. Mater. 2018, 360, 172–181. [Google Scholar] [CrossRef] [PubMed]
  27. Mandal, A.; Sarkar, B.; Owens, G.; Thakur, J.K.; Manna, M.C.; Niazi, N.K.; Jayaraman, S.; Patra, A.K. Impact of Genetically Modified Crops on Rhizosphere Microorganisms and Processes: A Review Focusing on Bt Cotton. Appl. Soil Ecol. 2020, 148, 103492. [Google Scholar] [CrossRef]
  28. Stotzky, G. Persistence and Biological Activity in Soil of the Insecticidal Proteins from Bacillus thuringiensis, Especially from Transgenic Plants. Plant Soil 2005, 266, 77–89. [Google Scholar] [CrossRef]
  29. Sarkar, B.; Patra, A.K.; Purakayastha, T.J.; Megharaj, M. Assessment of Biological and Biochemical Indicators in Soil under Transgenic Bt and Non-Bt Cotton Crop in a Sub-Tropical Environment. Environ. Monit. Assess. 2009, 156, 595–604. [Google Scholar] [CrossRef]
  30. Fleming, D.; Musser, F.; Reisig, D.; Greene, J.; Taylor, S.; Parajulee, M.; Lorenz, G.; Catchot, A.; Gore, J.; Kerns, D.; et al. Effects of Transgenic Bacillus thuringiensis Cotton on Insecticide Use, Heliothine Counts, Plant Damage, and Cotton Yield: A Meta-Analysis, 1996–2015. PLoS ONE 2018, 13, e0200131. [Google Scholar] [CrossRef] [Green Version]
  31. Luo, J.; Zhang, S.; Zhu, X.; Lu, L.; Wang, C.; LI, C.; Cui, J.; Zhou, Z. Effects of Soil Salinity on Rhizosphere Soil Microbes in Transgenic Bt Cotton Fields. J. Integr. Agric. 2017, 16, 1624–1633. [Google Scholar] [CrossRef] [Green Version]
  32. Dunfield, K.E.; Germida, J.J. Impact of Genetically Modified Crops on Soil- and Plant-Associated Microbial Communities. J. Environ. Qual. 2004, 33, 806. [Google Scholar] [CrossRef]
  33. Sarkar, B.; Patra, A.K.; Purakayastha, T.J. Transgenic Bt-Cotton Affects Enzyme Activity and Nutrient Availability in a Sub-Tropical Inceptisol. J. Agron. Crop Sci. 2008, 194, 289–296. [Google Scholar] [CrossRef]
  34. Sun, C.X.; Chen, L.J.; Wu, Z.J.; Zhou, L.K.; Shimizu, H. Soil Persistence of Bacillus thuringiensis (Bt) Toxin from Transgenic Bt Cotton Tissues and Its Effect on Soil Enzyme Activities. Biol. Fertil. Soils 2007, 43, 617–620. [Google Scholar] [CrossRef]
  35. Noman, A.; Bashir, R.; Aqeel, M.; Anwer, S.; Iftikhar, W.; Zainab, M.; Zafar, S.; Khan, S.; Islam, W.; Adnan, M. Success of Transgenic Cotton (Gossypium hirsutum L.): Fiction or Reality? Cogent Food Agric. 2016, 2, 1207844. [Google Scholar] [CrossRef]
  36. Flachs, A. Transgenic cotton: High hopes and farming reality. Nat. Plants 2017, 3, 16212. [Google Scholar] [CrossRef]
  37. Smyth, S.J. The Human Health Benefits from GM Crops. Plant Biotechnol. J. 2020, 18, 887–888. [Google Scholar] [CrossRef]
  38. Matloob, A.; Aslam, F.; Rehman, H.U.; Khaliq, A.; Ahmad, S.; Yasmeen, A.; Hussain, N. Cotton-Based Cropping Systems and Their Impacts on Production. In Cotton Production and Uses; Springer: Singapore, 2020; pp. 283–310. ISBN 9789811514722. [Google Scholar]
  39. Kroetsch, D.; Wang, C. Particle Size Distribution. Soil Sampl. Methods Anal. 2008, 2, 713–725. [Google Scholar]
  40. Dellavalle, N.B. Determination of Soil-Paste PH and Conductivity of Saturation Extract. In Reference Methods for Soil Analysis; Soil and Plant Analysis Council, Inc.: Athens, GA, USA, 1992; pp. 40–43. [Google Scholar]
  41. Cunniff, P.; AOAC International. Official Methods of Analysis of AOAC International, 16th ed.; AOAC International: Gaithersburg, MD, USA, 1997. [Google Scholar]
  42. Hoogsteen, M.J.J.; Lantinga, E.A.; Bakker, E.J.; Groot, J.C.J.; Tittonell, P.A. Estimating Soil Organic Carbon through Loss on Ignition: Effects of Ignition Conditions and Structural Water Loss. Eur. J. Soil Sci. 2015, 66, 320–328. [Google Scholar] [CrossRef]
  43. Jafari, A.; Connolly, V.; Frolich, A.; Walsh, E.J. A Note on Estimation of Quality Parameters in Perennial Ryegrass by near Infrared Reflectance Spectroscopy. Ir. J. Agric. Food Res. 2003, 42, 293–299. [Google Scholar]
  44. CIMMYT. From Agronomic Data to Farmer Recommendations: An Economics Workbook; CIMMYT: Heroica Veracruz, Mexico, 1988; ISBN 9686127194. [Google Scholar]
  45. Shapiro, S.S.; Wilk, M.B. An Analysis of Variance Test for Normality (Complete Samples). Biometrika 1965, 52, 591–611. [Google Scholar] [CrossRef]
  46. Steel, R.; Torrei, J.; Dickey, D. Principles and Procedures of Statistics a Biometrical Approach; McGraw-Hill College: New York, NY, USA, 1997. [Google Scholar]
  47. IBM, Inc. SPSS Statistics for Windows (Version 20); IBM SPSS Inc.: Armonk, NY, USA, 2012; pp. 1–8. [Google Scholar]
  48. Kranthi, K.R.; Stone, G.D. Long-Term Impacts of Bt Cotton in India. Nat. Plants 2020, 6, 188–196. [Google Scholar] [CrossRef]
  49. Ahmad, I.; Zhou, G.; Zhu, G.; Ahmad, Z.; Song, X.; Hao, G.; Jamal, Y.; Ibrahim, M.E.H. Response of Leaf Characteristics of BT Cotton Plants to Ratio of Nitrogen, Phosphorus and Potassium. Pak. J. Bot. 2021, 53, 873–881. [Google Scholar] [CrossRef]
  50. Marschner, P. Marschner’s Mineral Nutrition of Higher Plants, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2011; ISBN 9780123849052. [Google Scholar]
  51. Jones, M.L.M.; Wallace, H.L.; Norris, D.; Brittain, S.A.; Haria, S.; Jones, R.E.; Rhind, P.M.; Reynolds, B.R.; Emmett, B.A. Changes in Vegetation and Soil Characteristics in Coastal Sand Dunes along a Gradient of Atmospheric Nitrogen Deposition. Plant Biol. 2004, 6, 598–605. [Google Scholar] [CrossRef] [Green Version]
  52. Marschner, P.; Fu, Q.; Rengel, Z. Manganese Availability and Microbial Populations in the Rhizosphere of Wheat Genotypes Differing in Tolerance to Mn Deficiency. J. Plant Nutr. Soil Sci. 2003, 166, 712–718. [Google Scholar] [CrossRef]
  53. Hu, H.; Xie, M.; Yu, Y.; Zhang, Q. Transgenic Bt Cotton Tissues Have No Apparent Impact on Soil Microorganisms. Plant Soil Environ. 2013, 59, 366–371. [Google Scholar] [CrossRef] [Green Version]
  54. O’Callaghan, M.; Glare, T.R.; Burgess, E.P.J.; Malone, L.A. Effects of Plants Genetically Modified for Insect Resistance on Nontarget Organisms. Annu. Rev. Entomol. 2005, 50, 271–292. [Google Scholar] [CrossRef] [Green Version]
  55. Beura, K.; Rakshit, A. Effect of Bt Cotton on Nutrient Dynamics under Varied Soil Type. Ital. J. Agron. 2011, 6, e35. [Google Scholar] [CrossRef] [Green Version]
  56. Kumari, S.; Manjhi, B.K.; Beura, K.S.; Rakshit, A. Decomposition Bt Cotton Residues Affecting Soil Microbial Activity under Varied Soils. Int. J. Agric. Environ. Biotechnol. 2015, 8, 359. [Google Scholar] [CrossRef]
  57. Shahzad, M.; Farooq, M.; Hussain, M. Weed Spectrum in Different Wheat-Based Cropping Systems under Conservation and Conventional Tillage Practices in Punjab, Pakistan. Soil Tillage Res. 2016, 163, 71–79. [Google Scholar] [CrossRef]
  58. Shahzad, M.; Hussain, M.; Jabran, K.; Farooq, M.; Farooq, S.; Gašparovič, K.; Barboricova, M.; Aljuaid, B.S.; El-Shehawi, A.M.; Zuan, A.T.K. The Impact of Different Crop Rotations by Weed Management Strategies’ Interactions on Weed Infestation and Productivity of Wheat (Triticum aestivum L.). Agronomy 2021, 11, 2088. [Google Scholar] [CrossRef]
  59. Tariq, M.; Abdullah, K.; Ahmad, S.; Abbas, G.; Rahman, M.H.; Khan, M.A. Weed Management in Cotton. In Cotton Production and Uses; Springer: Singapore, 2020; pp. 145–161. [Google Scholar]
  60. Farkas, A. Soil Management and Tillage Possibilities in Weed Control. Herbologia 2006, 7, 9–23. [Google Scholar]
  61. Mahajan, G.; Chauhan, B.S. The Role of Cultivars in Managing Weeds in Dry-Seeded Rice Production Systems. Crop Prot. 2013, 49, 52–57. [Google Scholar] [CrossRef]
  62. Rezakhanlou, A.; Mirshekari, B.; Zand, E.; Farahvash, F.; Baghestani, M.A. Evaluation of Competitiveness of Cotton Varieties to Cocklebur (Xanthium srumarium L.). J. Food Agric. Environ. 2013, 11, 308–311. [Google Scholar]
  63. Chandler, J.M.; Meredith, W.R. Yields of Three Cotton (Gossypium hirsutum) Cultivars as Influenced by Spurred Anoda (Anoda cristata) Competition. Weed Sci. 1983, 31, 303–307. [Google Scholar] [CrossRef]
  64. Tabashnik, B.E.; Carrière, Y.; Dennehy, T.J.; Morin, S.; Sisterson, M.S.; Roush, R.T.; Shelton, A.M.; Zhao, J.-Z. Insect Resistance to Transgenic Bt Crops: Lessons from the Laboratory and Field. J. Econ. Entomol. 2003, 96, 1031–1038. [Google Scholar] [CrossRef]
  65. Heckel, D.G. How Do Toxins from Bacillus thuringiensis Kill Insects? An Evolutionary Perspective. Arch. Insect Biochem. Physiol. 2020.
  66. Tapp, H.; Stotzky, G. Insecticidal Activity of the Toxins from Bacillus thuringiensis Subspecies Kurstaki and Tenebrionis Adsorbed and Bound on Pure and Soil Clays. Appl. Environ. Microbiol. 1995, 61, 1786–1790. [Google Scholar] [CrossRef] [Green Version]
  67. Rui, Y.-K.; Yi, G.-X.; Zhao, J.; Wang, B.-M.; Li, Z.-H.; Zhai, Z.-X.; He, Z.-P.; Li, Q.X. Changes of Bt Toxin in the Rhizosphere of Transgenic Bt Cotton and Its Influence on Soil Functional Bacteria. World J. Microbiol. Biotechnol. 2005, 21, 1279–1284. [Google Scholar] [CrossRef]
  68. Wei, X.-D.; Zou, H.-L.; Chu, L.-M.; Liao, B.; Ye, C.-M.; Lan, C.-Y. Field Released Transgenic Papaya Affects Microbial Communities and Enzyme Activities in Soil. Plant Soil 2006, 285, 347–358. [Google Scholar] [CrossRef]
  69. Fang, M.; Motavalli, P.P.; Kremer, R.J.; Nelson, K.A. Assessing Changes in Soil Microbial Communities and Carbon Mineralization in Bt and Non-Bt Corn Residue-Amended Soils. Appl. Soil Ecol. 2007, 37, 150–160. [Google Scholar] [CrossRef]
Figure 1. The densities (±standard errors of means, n = 3) of broadleaved, narrow-leaved, and total weed species recorded in different cotton genotypes. Here, Bt and NBt stand for Bt and non-Bt genotypes.
Figure 1. The densities (±standard errors of means, n = 3) of broadleaved, narrow-leaved, and total weed species recorded in different cotton genotypes. Here, Bt and NBt stand for Bt and non-Bt genotypes.
Agriculture 13 00276 g001
Figure 2. The impacts of different winter crops on weed density (±standard errors of means, n = 3) in different winter crops included in the study.
Figure 2. The impacts of different winter crops on weed density (±standard errors of means, n = 3) in different winter crops included in the study.
Agriculture 13 00276 g002
Table 1. Physicochemical characteristics of the soil in the experimental area before the initiation of the experiments.
Table 1. Physicochemical characteristics of the soil in the experimental area before the initiation of the experiments.
Soil PropertiesUnit2016–20172017–2018
Organic matter content%0.590.56
Total nitrogen (N)kg ha−122.1222.23
Available phosphorus (P)kg ha−118.0218.08
Available potassium (K)kg ha−1245.15249.15
pH 8.178.19
ECdS m−14.965.00
Silt%54.1554.00
Sand%25.7526.10
Clay%20.1019.90
Table 2. The production practices used for the cultivation of cotton and winter crops used in study.
Table 2. The production practices used for the cultivation of cotton and winter crops used in study.
Crops NameGenotype NamePlanting Time *Seed Rate (kg ha−1)Fertilizer
NPK (kg ha−1)
R × R
(cm)
P × P
(cm)
Harvesting
Time
CottonGH Mubarik and CIM-616 (Bt)
CIM-620 and CIM-554
(non-Bt)
08 and 10 May25250-175-125 (Bt)
200-145-100 (non-Bt)
7520Last picking in October
WheatGalaxy-201313 and 16 November125130-100-6225 21 and 23 April
CanolaHyola-42012 and 13 November590-60-50304-56 and 10 April
Egyptian cloverAnmol berseem9 and 11 November2522-115-0 Last cutting in April
Different dates in planting and harvesting time column indicate the dates in first and second years of the experiment, R × R = row to row spacing, P × P = plant to plant spacing. *, the first and second dates denote the planting time of the respective crop during 1st and 2nd year of the study, respectively.
Table 3. The influences of cotton genotypes from winter crops interaction on nutrient availability and soil organic matter contents in the soil after the harvest of winter crops.
Table 3. The influences of cotton genotypes from winter crops interaction on nutrient availability and soil organic matter contents in the soil after the harvest of winter crops.
Treatments2016–20172017–2018
WheatEgyptian CloverCanolaWheatEgyptian CloverCanola
Available nitrogen (kg ha−1)
CIM-616 (Bt1)0.17 ± 0.001 a–c0.15 ± 0.003 c–e0.14 ± 0.001 de0.18 ± 0.003 ab0.16 ± 0.005 b–d0.16 ± 0.001 b–d
GH-Mubarik (Bt2)0.19 ± 0.003 a0.14 ± 0.002 de0.16 ± 0.001 b–d0.19 ± 0.002 a0.16 ± 0.004 b–d0.17 ± 0.003 a–c
CIM-620 (NBt1)0.18 ± 0.002 ab0.15 ± 0.001 c–e0.14 ± 0.004 de0.18 ± 0.001 ab0.14 ± 0.003 cd0.14 ± 0.002 d
N-414 (NBt2)0.16 ± 0.004 b–d0.15 ± 0.002 c–e0.13 ± 0.002 e0.18 ± 0.002 ab0.16 ± 0.002 b–d0.16 ± 0.002 b–d
LSD (p ≤ 0.05)0.0200.020
Available phosphorous (kg ha−1)
CIM-616 (Bt1)19.36 ± 0.02 a–e19.50 ± 0.03 a–c19.38 ± 0.02 a–e19.40 ± 0.01a–c19.59 ± 0.07 a19.28 ± 0.04 b–d
GH-Mubarik (Bt2)19.20 ± 0.04 de19.52 ± 0.02 ab19.24 ± 0.05 de19.10 ± 0.02 de19.42 ± 0.04 ab19.14 ± 0.06 de
CIM-620 (NBt1)19.28 ± 0.06 c–e19.58 ± 0.07 a19.30 ± 0.04 b–e19.18 ± 0.06 de19.48 ± 0.02 ab19.20 ± 0.05 c–e
N-414 (NBt2)19.18 ± 0.04 e19.40 ± 0.03 a–d19.40 ± 0.03 a–d19.08 ± 0.05 e19.30 ± 0.03 a–d19.30 ± 0.04 a–d
LSD (p ≤ 0.05)0.100.12
Available potassium (kg ha−1)
CIM-616 (Bt1)394 ± 6.1 b–d400 ± 4.3 ab394 ± 3.3 e402 ± 2.2 a–c404 ± 2.2 a–c406 ± 3.2 a–c
GH-Mubarik (Bt2)388 ± 5.3 de396 ± 6.1 a–c394 ± 3.2 e408 ± 6.1 a402 ± 3.1 a–c404 ± 2.6 a–c
CIM-620 (NBt1)390 ± 3.4 c–e398 ± 3.3 ab386 ± 8.3 de400 ± 3.4 bc402 ± 3.4 a–c402 ± 2.7 a–c
N-414 (NBt2)392 ± 2.4 b–e402 ± 1.2 a390 ± 4.5 c–e400 ± 3.3 bc406 ± 4.3 ab400 ± 2.3 c
LSD (p ≤ 0.05)7.586.30
Available zinc (kg ha−1)
CIM-616 (Bt1)1.46 ± 0.01 d1.60 ± 0.03 b1.44 ± 0.03 de1.56 ± 0.04 cd1.60 ± 0.02 bc1.48 ± 0.02 ef
GH-Mubarik (Bt2)1.46 ± 0.02 d1.58 ± 0.02 bc1.40 ± 0.02 e1.50 ± 0.04 de1.62 ± 0.02 b1.42 ± 0.03 f
CIM-620 (NBt1)1.58 ± 0.02 bc1.66 ± 0.01 a1.44 ± 0.04 c1.60 ± 0.03 bc1.68 ± 0.04 a1.56 ± 0.04 cd
N-414 (NBt2)1.58 ± 0.01 bc1.68 ± 0.02 a1.56 ± 0.02 bc1.58 ± 0.02 bc1.68 ± 0.05 a1.58 ± 0.02 bc
LSD (p ≤ 0.05)0.040.06
Available iron (kg ha−1)
CIM-616 (Bt1)7.62 ± 0.12 d–f7.42 ± 0.10 fg7.72 ± 0.09 c–e7.68 ± 0.11 de7.78 ± 0.13 c–e7.78 ± 0.10 c–e
GH-Mubarik (Bt2)7.82 ± 0.14 b–e7.68 ± 0.11 de7.84 ± 0.11 b–d7.90 ± 0.12 bc8.04 ± 0.17 ab7.80 ± 0.14 c–e
CIM-620 (NBt1)8.14 ± 0.19 a7.32 ± 0.09 g7.94 ± 0.11 a–c8.20 ± 0.16 a7.96 ± 0.11 bc7.80 ± 0.13 c–e
N-414 (NBt2)8.04 ± 0.11 ab7.56 ± 0.14 ef7.60 ± 0.09 ef7.96 ± 0.10 bc7.84 ± 0.12 b–d7.62 ± 0.12 e
LSD (p ≤ 0.05)0.240.20
Soil organic matter (%)
CIM-616 (Bt1)0.59 ± 0.02 a0.53 ± 0.04 c–e0.51 ± 0.02 de0.62 ± 0.01 ab0.62 ± 0.01 ab0.59 ± 0.01 cd
GH-Mubarik (Bt2)0.58 ± 0.03 ab0.53 ± 0.04 c–e0.51 ± 0.02 de0.60 ± 0.02 bc0.63 ± 0.01 a0.60 ± 0.01 bc
CIM-620 (NBt1)0.57 ± 0.04 a–c0.52 ± 0.03 de0.49 ± 0.01 ef0.59 ± 0.02 cd0.59 ± 0.01 cd0.57 ± 0.01 d
N-414 (NBt2)0.58 ± 0.02 ab0.54 ± 0.04 b–d0.45 ± 0.02 f0.58 ± 0.01 cd0.58 ± 0.02 cd0.60 ± 0.01 bc
LSD (p ≤ 0.05)0.040.03
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 4. Interactive effects of different cotton genotypes and winter crops on the densities of broadleaved, narrow-leaved, and total weeds recorded in winter crops.
Table 4. Interactive effects of different cotton genotypes and winter crops on the densities of broadleaved, narrow-leaved, and total weeds recorded in winter crops.
Treatments2016–20172017–2018
WheatEgyptian CloverCanolaWheatEgyptian CloverCanola
Broadleaved weeds density (m−2)
CIM-616 (Bt1)81.3 ± 3.1 b61.0 ± 3.1 f62.3 ± 2.8 f83.0 ± 3.4 cd66.7 ± 3.9 g65.0 ± 3.1 g
GH-Mubarik (Bt2)68.7 ± 4.3 e60.3 ± 2.4 f64.0 ± 2.9 f74.3 ± 2.6 ef69.3 ± 5.2 fg72.7 ± 3.4 ef
CIM-620 (NBt1)90.3 ± 3.4 a77.7 ± 2.8 bc75.3 ± 2.4 cd97.0 ± 6.3 a88.7 ± 2.8 b83.7 ± 5.1 bc
N-414 (NBt2)76.0 ± 2.7 cd72.7 ± 3.3 de69.7 ± 1.4 e82.0 ± 4.7 cd81.0 ± 3.0 cd78.0 ± 3.1 de
LSD (p ≤ 0.05)4.425.57
Narrow-leaved weeds density (m−2)
CIM-616 (Bt1)57.0 ± 2.2 a16.0 ± 3.6 e27.3 ± 2.1 c64.0 ± 4.6 a24.0 ± 3.4 e35.3 ± 4.4 c
GH-Mubarik (Bt2)48.0 ± 3.3 b12.7 ± 3.4 ef20.3 ± 3.4d55.7 ± 3.3 b22.0 ±3.3 e28.0 ± 3.2 d
CIM-620 (NBt1)56.0 ± 4.5 a11.7 ± 2.6 f24.3 ± 2.9 c66.3 ± 4.1 a21.7 ± 4.5 e34.3 ± 4.0 c
N-414 (NBt2)46.7 ± 2.6 b20.0 ± 2.9 d14.3 ± 1.7 ef55.7 ± 4.8 b28.0 ± 3.7 d25.0 ± 4.2 de
LSD (p ≤ 0.05)3.473.53
Total weeds density (m−2)
CIM-616 (Bt1)138 ± 6 b77.0 ± 8 h89.7 ± 3 f147 ± 10 b90.7 ± 5 h100 ± 5 g
GH-Mubarik (Bt2)116 ± 5 d73.0 ± 11 i84.3 ± 2 g130 ± 11 d91.3 ± 6 h100 ± 6 g
CIM-620 (NBt1)146 ± 10 a89.3 ± 4 f99.7 ± 8 e163 ± 19 a110 ± 8 f118 ± 11 e
N-414 (NBt2)122 ± 4 c92.7 ± 5 f84.0 ± 4 g137 ± 6 c109 ± 7 f103 ± 5 g
LSD (p ≤ 0.05)3.845.43
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 5. The influences of different Bt and non-Bt genotypes on yield-related attributes of cotton crop.
Table 5. The influences of different Bt and non-Bt genotypes on yield-related attributes of cotton crop.
Treatments201620172016201720162017
Monopodial Branches (Plant−1)Sympodial Branches (Plant−1)Boll Weight (g)
CIM-616 (Bt1)1.78 ± 0.21.78 ± 0.323.6 ± 0.9 b26.0 ± 1.1 a3.2 ± 0.1 a3.2 ± 0.04 a
GH-Mubarik (Bt2)1.67 ± 0.31.67 ± 0.422.3 ± 0.8 b25.6 ± 1.2 a3.1 ± 0.1 ab3.1 ± 0.02 b
CIM-620 (NBt1)1.89 ± 0.31.89 ± 0.225.0 ± 1.4 a23.0 ± 0.8 b3.0 ± 0.05 bc3.0 ± 0.08 b
N-414 (NBt2)1.67 ± 0.41.67 ± 0.322.0 ± 1.1 bc24.0 ± 1.2 b2.9 ± 0.04 c2.9 ± 005 c
LSD (p ≤ 0.05)NSNS1.491.290.150.07
Seed cotton yield (kg ha−1)Harvest index (%)
CIM-616 (Bt1)2892 ± 141 a2832 ± 22132.7 ± 2.1232.4 ± 1.9
GH-Mubarik (Bt2)2685 ± 123 b2635 ± 21329.5 ± 2.2129.4 ± 3.1
CIM-620 (NBt1)2645 ± 129 b2563 ± 30333.9 ± 2.3931.2 ± 2.2
N-414 (NBt2)2613 ± 147 b2570 ± 30932.3 ± 2.5331.1 ± 2.6
LSD (p ≤ 0.05)150.14NSNSNS
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 6. Yield-related parameters of wheat crop sown after the harvest of transgenic and non-transgenic cotton genotypes.
Table 6. Yield-related parameters of wheat crop sown after the harvest of transgenic and non-transgenic cotton genotypes.
Treatments2016–20172017–20182016–20172017–20182016–20172017–2018
Productive Tillers (m−2)Grains (Spike−1)1000-Grain Weight (g)
CIM-616 (Bt1)189 ± 15 ab191 ± 17 ab55.7 ± 1.9 b56.0 ± 2.0 b36.2 ± 1.8 c36.8 ± 1.4 c
GH-Mubarik (Bt2)181 ± 11 b176 ± 14 b53.3 ± 1.7 b53.9 ± 1.8 c37.9 ± 1.7 bc37.9 ± 1.6 bc
CIM-620 (NBt1)197 ± 10 a201 ± 14 a59.5 ± 1.2 a58.8 ± 1.4 a40.2 ± 1.5 a40.7 ± 1.3 a
N-414 (NBt2)202 ± 12 a202 ± 16 a59.1 ± 1.4 a58.1 ± 1.6 a39.8 ± 1.6 ab39.6 ± 2.3 ab
LSD (p ≤ 0.05)14.414.42.41.72.01.8
Grain yield (t ha−1)Biological yield (t ha−1)Harvest index (%)
CIM-616 (Bt1)5.82 ± 0.85.95 ± 0.2 b17.6 ± 0.5 bc15.7 ± 0.4 bc33.1 ± 1.237.8 ± 2.1
GH-Mubarik (Bt2)5.98 ± 0.75.92 ± 0.1 b17.1 ± 0.6 c15.3 ± 0.5 c34.9 ± 1.638.8 ± 2.2
CIM-620 (NBt1)6.30 ± 0.76.26 ± 0.2 a18.2 ± 0.6 ab16.4 ± 0.6 ab34.6 ± 1.838.2 ± 2.1
N-414 (NBt2)6.21 ± 0.86.31 ± 0.2 a18.7 ± 0.5 a16.9 ± 0.5 a33.2 ± 2.037.4 ± 2.4
LSD (p ≤ 0.05)NS0.250.860.86NSNS
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 7. Yield-related traits of canola sown after different transgenic and non-transgenic cotton genotypes.
Table 7. Yield-related traits of canola sown after different transgenic and non-transgenic cotton genotypes.
Treatments2016–20172017–20182016–20172017–20182016–20172017–2018
Siliques (Plant−1)Seeds (Silique−1)1000-Seed Weight (g)
CIM-616 (Bt1)106 ± 22105 ± 7 ab26.7 ± 2.326.9 ± 3.72.77 ± 0.32.73 ± 0.4
GH-Mubarik (Bt2)103 ± 12102 ± 9 b24.0 ± 3.425.0 ± 4.02.90 ± 0.42.85 ± 0.3
CIM-620 (NBt1)109 ± 6109 ± 7 ab25.2 ± 3.625.8 ± 3.12.87 ± 0.62.90 ± 0.2
N-414 (NBt2)112 ± 16113 ± 8 a26.3 ± 4.127.0 ± 3.32.83 ± 0.52.93 ± 0.5
LSD (p ≤ 0.05)NS8.4NSNSNSNS
Biological yield (kg ha−1)Seed yield (kg ha−1)Harvest index (%)
CIM-616 (Bt1)4800 ± 3435271 ± 2341650 ± 212 b1797 ± 15834.4 ± 2.234.1 ± 3.4
GH-Mubarik (Bt2)5132 ± 4125070 ± 2671700 ± 223 b1833 ± 12333.2 ± 3.136.2 ± 2.1
CIM-620 (NBt1)5233 ± 3455345 ± 3121950 ± 201 a1850 ± 11237.4 ± 2.634.7 ± 3.3
N-414 (NBt2)4876 ± 3215478 ± 4341900 ± 198 a1900 ± 12139.1 ± 2.834.7 ± 3.1
LSD (p ≤ 0.05)NSNS197.7NSNSNS
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 8. The influence of various cotton varieties on yield-related traits of Egyptian clover.
Table 8. The influence of various cotton varieties on yield-related traits of Egyptian clover.
TreatmentsFresh Forage Yield (t ha−1)Dry forage Yield (t ha−1)Crude Protein (%)
2016–20172017–20182016–20172017–20182016–20172017–2018
CIM-616 (Bt1)28.3 ± 1.21 b30.7 ± 1.98 b2.91±0.11 b3.62 ± 0.09 b21.0 ± 1.2 b20.6 ± 1.7 b
GH-Mubarik (Bt2)28.2 ± 1.26 b32.0 ± 2.02 b2.97 ± 0.17 b3.72 ± 0.16 ab20.2 ± 1.6 b20.3 ± 2.4 b
CIM-620 (NBt1)34.1 ± 2.34 a34.8 ± 1.12 a3.50 ± 0.21 a3.87 ± 0.11 a24.0 ± 2.1 a22.3 ± 1.6 ab
N-414 (NBt2)32.3 ± 2.31 a33.2 ± 1.63 ab3.35 ± 0.18 a3.76 ± 0.12 ab23.7 ± 1.9 a23.6 ± 1.2 a
LSD (p ≤ 0.05)1.962.630.180.172.442.52
Means followed by different letters significantly differ (p ≤ 0.05) from each other. The values of different traits are means of three replications ± standard errors of means. NS = non-significant; Bt = transgenic genotypes; NBt = conventional or non-transgenic genotypes.
Table 9. The impacts of different transgenic and non-transgenic cotton genotypes on system productivity of various cotton-based cropping systems.
Table 9. The impacts of different transgenic and non-transgenic cotton genotypes on system productivity of various cotton-based cropping systems.
Treatments2016–20172017–2018
TEGINIBCRTEGINIBCR
Bt1 × Wheat1563.592607.751044.161.671563.592531.58967.991.62
Bt2 × Wheat1563.592516.55952.961.611563.592423.59860.001.55
NBt1 × Wheat1629.852572.23942.381.581629.852466.93837.091.51
NBt2 × Wheat1629.852564.07934.221.571629.852491.97862.121.53
Bt1 × Canola1500.761979.11478.351.321500.762009.13508.371.34
Bt2 × Canola1500.761910.59409.831.271500.761923.63422.871.28
NBt1 × Canola1567.021966.87399.851.261567.021904.60337.591.22
NBt2 × Canola1567.021925.86358.851.231567.021926.66359.651.23
Bt1 × Egyptian clover1621.021989.23368.211.231621.022016.24395.221.24
Bt2 × Egyptian clover1621.021893.59272.571.171621.021957.23336.211.21
NBt1 × Egyptian clover1687.282010.18322.901.191687.281987.73300.451.18
NBt2 × Egyptian clover1687.281953.95266.681.161687.281954.70267.421.16
BCR = benefit–cost ratio; Bt1 = CIM-616; Bt2 = GH-Mubarik; NBt1 = CIM-620; NBt2 = N-414; Bt = transgenic cotton; NBt = non-transgenic cotton; TE= total expenditure; GI= gross income; NI= net income, the values of TE, GI, and NI are in US$.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marral, M.W.R.; Ahmad, F.; Ul-Allah, S.; Atique-ur-Rehman; Farooq, S.; Hussain, M. Influence of Transgenic (Bt) Cotton on the Productivity of Various Cotton-Based Cropping Systems in Pakistan. Agriculture 2023, 13, 276. https://doi.org/10.3390/agriculture13020276

AMA Style

Marral MWR, Ahmad F, Ul-Allah S, Atique-ur-Rehman, Farooq S, Hussain M. Influence of Transgenic (Bt) Cotton on the Productivity of Various Cotton-Based Cropping Systems in Pakistan. Agriculture. 2023; 13(2):276. https://doi.org/10.3390/agriculture13020276

Chicago/Turabian Style

Marral, Muhammad Waseem Riaz, Fiaz Ahmad, Sami Ul-Allah, Atique-ur-Rehman, Shahid Farooq, and Mubshar Hussain. 2023. "Influence of Transgenic (Bt) Cotton on the Productivity of Various Cotton-Based Cropping Systems in Pakistan" Agriculture 13, no. 2: 276. https://doi.org/10.3390/agriculture13020276

APA Style

Marral, M. W. R., Ahmad, F., Ul-Allah, S., Atique-ur-Rehman, Farooq, S., & Hussain, M. (2023). Influence of Transgenic (Bt) Cotton on the Productivity of Various Cotton-Based Cropping Systems in Pakistan. Agriculture, 13(2), 276. https://doi.org/10.3390/agriculture13020276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop