Introduction
Dicamba-resistant cotton (Gossypium hirsutum L.) cultivars have been commercialized in the United States, and the adoption of varieties with this trait has been widespread in cotton production, because their use enables management of weed species with evolved resistance to glyphosate and acetolactate synthase (ALS) inhibitors, such as Palmer amaranth (Amaranthus palmeri S. Watson). The dicamba-resistance trait is being incorporated into most glyphosate- and glufosinate-based herbicide-resistant (HR) cultivars (Cahoon et al. Reference Cahoon, York, Jordan, Everman, Seagroves, Culpepper and Eure2015b; Merchant et al. Reference Merchant, Sosnoskie, Culpepper, Steckel, York, Braxton and Ford2013; Meyer et al. Reference Meyer, Norsworthy, Young, Steckel, Bradley, Johnson, Loux, Davis, Kruger, Bararpour, Ikley, Spaunhorst and Butts2015). The increased use of dicamba will increase selection pressure on weed communities and potentially lead to weed species shifts over time (Culpepper Reference Culpepper2006; Shaner Reference Shaner2000). Diversifying cropping systems, integrating cultural and mechanical weed management practices, and reducing overall herbicide use of the same herbicide are required to reduce selection pressure and slow down the evolutionary process leading to herbicide resistance (Boerboom Reference Boerboom1999; Kruger et al. Reference Kruger, Johnson, Weller, Owen, Shaw, Wilcut, Jordan, Wilson, Bernards and Young2009; Neve et al. Reference Neve, Norsworthy, Smith and Zelaya2011; Norsworthy et al. Reference Norsworthy, Ward, Shaw, Llewellyn, Nichols, Webster, Bradley, Frisvold, Powles, Burgos, Witt and Barrett2012; Vencill et al. Reference Vencill, Nichols, Webster, Soteres, Mallory-Smith, Burgos, Johnson and McClelland2012). Furthermore, to steward dicamba-resistant cropping systems and keep them as a viable and effective technology for weed management, weed dynamics must be understood not only within one season but also over multiple seasons.
Introducing new management tools that include both chemical and nonchemical controls can produce weed species shifts (Cordeau et al. Reference Cordeau, Smith, Gallandt, Brown, Salon, DiTommaso and Ryan2017; Culpepper Reference Culpepper2006; de la Fuente et al. Reference de la Fuente, Oreja, Lenardis, Fuentes, Agosti, Barrio, Barberis, Robredo, Gil, Marzetti and Niccia2021 Tuesca et al. Reference Tuesca, Puricelli and Papa2001). The repeated use of a few herbicides over time in weed management can reduce weed diversity (Oreja et al. Reference Oreja, Inman, Jordan, Bardhan and Leon2022), resulting in a few dominant species that are difficult to manage, frequently exhibit herbicide resistance, and ultimately can cause important yield losses (Storkey and Neve Reference Storkey and Neve2018). Herbicide programs play a key role in weed shifts and the pace at which these occur (Menalled et al. Reference Menalled, Gross and Hammond2001). The inclusion of different mechanisms of action (MOAs), rather than a simple herbicide program with a few MOAs applied repeatedly, expands the spectrum of different weed species controlled, and the chance of selecting HR biotypes is reduced (Vencill et al. Reference Vencill, Nichols, Webster, Soteres, Mallory-Smith, Burgos, Johnson and McClelland2012). However, herbicides are not the only factors acting on community assembly; others, such as competition among weeds (Swanton et al. Reference Swanton, Clements and Derksen1993), can have a strong influence, depending on the species. There are no studies that report the weed density variability in dicamba-resistant cotton crops exposed to different herbicide programs that include the repetitive application of this herbicide.
Long term studies are needed to examine weed density variations in response to weed management decisions and determine overall trends. Such studies are useful tools for understanding both weed community changes over time in response to agronomic practices (Clements et al. Reference Clements, Weise and Swanton1994; Hobbs and Humphries Reference Hobbs and Humphries1995) and the selection pressure on HR biotypes. Therefore, the objectives of this research were (1) to determine whether herbicide programs with wider control spectra reduce weed species communities over time more so than simple herbicide programs with narrower control spectra; (2) to ascertain whether the effects on the population trajectories are maintained after changing herbicide programs; and (3) to compare A. palmeri population frequency of glyphosate resistance over 8 yr in response to herbicide programs that included glyphosate, dicamba, and residual herbicides.
Materials and Methods
Two experiments were established in two separate fields in a randomized complete block design with four replications each in 2011 at the Upper Coastal Plain Research Station in Rocky Mount, NC (35.893°N, 77.681°W). Cotton (BollGard II® XtendFlex®; Monsanto, St. Louis, MO) resistant to dicamba, glufosinate, glyphosate, and Lepidopteran insects (events MON-88701-6, MON-88913-6, MON-15985) was planted in conventionally tilled, raised beds at a seeding of 17 seed m−1 of row. Cotton in all years was planted in the second or third week of May. Other than treatments imposed for the experiment, cotton was managed according to North Carolina Cooperative Extension Service recommendations (Edmisten et al. Reference Edmisten, Bullen, Collins, Crozier, Reisig, Thiessen and Hardy2018). At the beginning of the study, fields were naturally infested with A. palmeri, including both glyphosate-susceptible (GS) and ALS-susceptible as well as glyphosate-resistant (GR) and ALS-resistant A. palmeri. The frequency of ALS resistance and glyphosate resistance was approximately 30% and <10%, respectively.
Seven herbicide programs composed of different MOAs were compared (Tables 1 and 2). Herbicides were applied using a CO2-pressurized backpack sprayer equipped with flat-fan nozzles (AIXR 11002, TeeJet® Technologies, Wheaton, IL) calibrated to deliver 140 L ha−1 at 152 kPa. Depending on the treatment, preemergence herbicides were applied immediately after planting (G+R, G+D+R, and G+½D+R), and postemergence herbicides were applied 2 (early postemergence), 4 (medium postemergence), and 6 (late postemergence) wk after planting.
Herbicide active ingredient, trade name, formulation, application rate, and manufacturer.
a Clarity® was used during 2011 to 2015. Engenia® was used during 2016 to 2018.
Preemergence and early (EPOST), medium (MPOST), and late (LPOST) postemergence herbicide treatments applied in Rocky Mount, NC, in 2011–2014. a
a All plots received glyphosate plus dicamba POST for the remainder of the study (2015–2018).
Weed Population Density
The germinable seedbank method was used to determine the weed density of each species (Reinhardt and Leon Reference Reinhardt and Leon2018). Ten soil cores (10.2 cm by 7.6 cm for a total volume of 4,630 ml) were collected from each plot before preemergence herbicide applications, in May of each year, except in 2015, when collection was done in January. Soil cores were placed in flats to a depth of 4 cm with a total surface area of 1,550 cm2, maintained in a climate-controlled greenhouse (25 to 35 C at 80% to 90% relative humidity), and irrigated periodically to promote germination and adequate seedling growth. A seedling count was done 3 wk after establishment, when no new seedlings were observed. Seedlings were identified, and the number of individuals per species was recorded.
Frequencies of HR Biotypes
To determine the frequency of GR A. palmeri individuals in flats from the soil cores each year, plants were treated with glyphosate at 946 g ae ha−1 after seedling count and identification, and the number of surviving plants was determined at 2 wk after application. In May 2019, at the end of the study, plants surviving glyphosate were treated with dicamba at 560 g ae ha−1 to determine whether there were any individuals with tolerance to dicamba after 8 yr of dicamba use. A known dicamba-susceptible A. palmeri population was included as a positive control.
Statistical Analysis
While 17 species were registered throughout the experiment (Table 3), those detected just a few times across all plots, treatments, and years were not considered in the analysis. To evaluate the effect of herbicide programs on seedling population density, data from two fields were pooled, maintaining the independence of each replicate. Seedling number per species was log-transformed before analysis and then subjected to regression procedures using GraphPad Prism 6 (GraphPad Software, San Diego, CA) program testing for linear, quadratic functions (Tables 4 and 5). The PROC Mixed procedure in SAS (v. 9.4; SAS Institute, Cary, NC) was used to analyze the log-transformed frequency of GR A. palmeri glyphosate resistance over time.
Weed species detected in the germinable seedbank during the experiment.
Equations and R-squared (R2) of the models for each species and the treatments glyphosate (G), glyphosate plus dicamba (G+D), glyphosate plus residual herbicides (diuron plus pendimethalin) (G+D+R), glyphosate plus dicamba plus acetochlor (G+D+Ac), glyphosate plus alternating dicamba between years plus residual herbicides (diuron plus pendimethalin) (G+½D+R), and glyphosate plus alternating dicamba between years (G+½D).
Equations and R-squared (R2) of the models for each species and the treatments glyphosate (G), glyphosate plus dicamba (G+D), glyphosate plus residual herbicides (diuron plus pendimethalin) (G+D+R), glyphosate plus dicamba plus acetochlor (G+D+Ac), glyphosate plus alternating dicamba between years plus residual herbicides (diuron plus pendimethalin) (G+½D+R), and glyphosate plus alternating dicamba between years (G+½D).
Results and Discussion
Weed Population Trajectories
Population density varied during the 96 mo that the experiment lasted, showing different patterns among weed species and/or herbicide programs (Figure 1). Some species increased, others decreased, and others maintained their population stability. Furthermore, after 48 mo, subjecting all populations to a single, uniform program based on the postemergence herbicides glyphosate and dicamba did not have the same effect on the population trajectories of all weed species studied.
Seedling population density for (A) Amaranthus palmeri, (B) Mollugo verticillata, (C) Digitaria sanguinalis, (D) Cyperus compressus, and (E) Eclipta prostrata from soil cores in response to herbicide programs: glyphosate (G), glyphosate plus dicamba (G+D), glyphosate plus residual herbicides (diuron plus pendimethalin) (G+D+R), glyphosate plus dicamba plus acetochlor (G+D+Ac), glyphosate plus alternating dicamba between years plus residual herbicides (diuron plus pendimethalin) (G+½D+R) and glyphosate plus alternating dicamba between years (G+½D), throughout the months after experiment initiation. Error bars represent standard error of the mean for each data point; an asterisk (*) indicates regression line slope is different from zero.
Treatments that included the use of dicamba every year (G+D, G+D+R, and G+D+Ac) were the only treatments that reduced the number of A. palmeri seedlings from the beginning of the experiment (Figure 1A; Table 4). The treatment in which only glyphosate postemergence was used during the first 4 yr (G) exhibited the highest number of A. palmeri seedlings, but this was expected, due to the presence of GR individuals in the field from the beginning of the study (Inman et al. Reference Inman, Jordan, York, Jennings, Monks, Everman, Bollman, Fowler, Cole and Soteres2016). In the first 4 yr, the addition of preemergence herbicides (G+R) helped maintain lower A. palmeri densities when compared with glyphosate-only postemergence (G), but the populations were still higher than in those treatments that included dicamba every year. When preemergence herbicides were used with alternating years of glyphosate and glyphosate plus dicamba (G+½D+R), lower A. palmeri densities were observed compared with the program without preemergence herbicides (G+½D) (Figure 1A). After the last 4 yr of glyphosate plus dicamba postemergence applications, a reduction in populations was observed in all treatments. By the end of the experiment, after 96 mo, A. palmeri densities did not differ across treatments (Figure 1A), due to the high efficacy of A. palmeri control with dicamba (Cahoon et al. Reference Cahoon, York, Jordan, Everman, Seagroves, Culpepper and Eure2015b; Tehranchian et al. Reference Tehranchian, Norsworthy, Powles, Bararpour, Bagavathiannan, Barber and Scott2017).
Carpetweed (Mollugo verticillata L.) and large crabgrass [Digitaria sanguinalis (L.) Scop.] populations showed a similar pattern; both were reduced in all treatments in the first 4 yr, with the decline of D. sanguinalis being more obvious. The structure of the herbicide program was not the critical factor affecting population density, as all programs exhibited the same trend. Digitaria sanguinalis seeds have short persistence in the soil seedbank (Oreja et al. Reference Oreja, Batlla and de la Fuente2020); if control is successful and no new seeds enter the soil seedbank by the end of season, populations should decline as observed at 24 and 36 mo after the beginning of the experiment (Figure 1C). Glyphosate applied alone (G) caused a more marked reduction of M. verticillata than other treatments (Figure 1B). This reduction could be due not only to the good control obtained with this herbicide (Culpepper and York Reference Culpepper and York2000; Van Gessel et al. Reference Van Gessel, Ayeni and Majek2001) but also by the absence of the antagonism of dicamba to glyphosate performance when both herbicides are applied together (Meyer et al. Reference Meyer, Peter, Norsworthy and Beffa2020). Diuron and pendimethalin controlled D. sanguinalis well (Cahoon et al. Reference Cahoon, York, Jordan, Everman, Seagroves, Braswell and Jennings2015a); and dicamba, pendimethalin, and acetochlor provided acceptable control of M. verticillata. Once all herbicide programs changed to glyphosate plus dicamba, an increase in seedling number was observed for both M. verticillata and D. sanguinalis (Figure 1B and C). These patterns seem to be not only the result of the change in the herbicide program but also of interspecific competition. Otherwise, the program with glyphosate plus dicamba during the last 4 yr would have maintained the decreasing trend observed during the first 4 yr. For example, the A. palmeri populations increased in most treatments during the first 4 yr (Figure 1A), while D. sanguinalis and M. verticillata populations decreased, and the opposite was observed during the last 4 yr of the study (Figure 1B and C). Interspecific competition strongly influences community assembly (Aschehoug et al. Reference Aschehoug, Brooker, Atwater, Maron and Callaway2016), so it seems that significant decreases in A. palmeri populations during the last 4 yr reduced interference on other weed species allowing those to recover their populations.
In the cases of annual sedge (Cyperus compressus L.) and eclipta [Eclipta prostrata (L.) L.] (Figure 1D and E), recovery after the unification of herbicide programs was minimal, almost null, and the behavior of the populations was more erratic (Figure 1D and E; see low R2 values in Table 5), probably due to a smaller seedbank and lower competitive ability than other weed species exhibiting faster population recovery. In contrast, corn spurry (Spergula arvensis L.) was practically absent in the first 4 yr, but its populations increased rapidly in the second half of the study (Figure 2D). This weed is usually sufficiently controlled with several of the herbicides used in the experiment, such as glyphosate (Ivany Reference Ivany2004), dicamba (Velloso and Dal’Piaz Reference Velloso and Dal’Piaz1982), and pendimethalin (Haar et al. Reference Haar, Fennimore and Lambert2001). Spergula arvensis is a summer annual species that can also behave as a winter annual, depending on location and due to seed dimorphism creating variable dormancy levels (Wagner Reference Wagner1988). Therefore, reduced levels of A. palmeri could free up resources for this plant to grow and increase its reproductive success, as was observed for D. sanguinalis and M. verticillata.
Density of (A) Eleusine indica, (B) Ambrosia artemisiifolia, (C) Chenopodium album, and (D) Spergula arvensis from soil cores with their respective regression lines for the different treatments: glyphosate (G), glyphosate plus dicamba (G+D), glyphosate plus residual herbicides (diuron plus pendimethalin) (G+D+R), glyphosate plus dicamba plus acetochlor (G+D+Ac), glyphosate plus alternating dicamba between years plus residual herbicides (diuron plus pendimethalin) (G+½D+R) and glyphosate plus alternating dicamba between years (G+½D), throughout the months after experiment initiation (Months). Error bars represent standard error of the mean for each data point; an asterisk (*) indicates regression line slope is different from zero.
No herbicide program reduced E. indica populations. Even under the most aggressive programs (G+D+Ac, G+½D+R, and G+½D), populations were more or less stable or tended to increase throughout the 96 mo (Figure 2A). After the first 4 yr, the treatment with glyphosate alone (G) had one of the lowest values of seedling density, but immediately after the herbicide programs were unified, the seedling numbers increased to levels that surpassed the original population (Figure 2A). Glyphosate is an effective herbicide for control of E. indica (Chuah et al. Reference Chuah, Salmijah and Ismail2004), but the antagonism exerted by dicamba (Meyer et al. Reference Meyer, Peter, Norsworthy and Beffa2020) could have reduced glyphosate phytotoxicity in this grass. Furthermore, this is the only species that seems to be unaffected by A. palmeri populations, while the rest of species exhibited low densities when A. palmeri density was high during the first 4 yr (Figure 2A). Despite its low stature, E. indica can maintain high reproductive rates (i.e., >300,000 seeds per plant; Ma et al. Reference Ma, Ma, Wu, Ren, Jiang and Ma2019), even under shaded conditions, by shifting photoassimilate allocation. For example, under 40% to 80% shading, this weed reduced total dry weight and tillers per plant and partitioned more resources to inflorescences and less to roots (Ismail et al. Reference Ismail, Chuah, Salmijah and Teng2003).
For species such as common ragweed (Ambrosia artemisiifolia L.) and common lambsquarters (Chenopodium album L.), initial densities were low, and all herbicide programs were successful at maintaining the populations at low levels (Figure 2B and C). Also, programs that included residual herbicides favored a faster decline in the populations of these species. It is known that pendimethalin can provide adequate levels of preemergence control of C. album (Alebrahim et al. Reference Alebrahim, Majd, Mohassel, Wilcockson, Baghestani, Ghorbani and Kudsk2012), with higher efficacy than acetochlor (Chomas and Kells Reference Chomas and Kells2004; Soltani et al. Reference Soltani, Nurse, Shropshire and Sikkema2013). Residual herbicides can also complement postemergence broad-spectrum herbicides by improving control of A. artemisiifolia (Armel et al. Reference Armel, Wilson, Richardson and Hines2003; Barnes et al. Reference Barnes, Knezevic, Sikkema, Lindquist and Jhala2017).
Frequency of Resistance
At initiation of the experiment, the frequency of GR A. palmeri ranged from 1% to 9% (Figure 3). After 1 yr, the frequency of GR A. palmeri increased in all programs, being greatest in the glyphosate-only program, but after 4 yr, there were no differences in glyphosate resistance frequency among treatments (Figure 3). Shergill et al. (Reference Shergill, Bish, Biggs and Bradley2018) reported evolution of GR giant ragweed (Ambrosia trifida L.) in a susceptible population after 4 yr of continuous use of glyphosate only. In this study, no differences in glyphosate resistance frequency were observed among herbicide programs after 8 yr (Figure 3), despite good control of GR A. palmeri plants using continuous glyphosate plus dicamba (Figure 1A). This was likely the result of GR pollen movement from adjacent plots and surrounding fields. It has been documented that GR A. palmeri pollen disperses up to 300 m under normal field conditions (Sosnoskie et al. Reference Sosnoskie, Webster, Kichler, MacRae, Grey and Culpepper2012).
Frequency of glyphosate resistance of Amaranthus palmeri from soil cores in response to the herbicide treatments: glyphosate (G), glyphosate plus dicamba (G+D), glyphosate plus residual herbicides (diuron plus pendimethalin) (G+D+R), glyphosate plus dicamba plus acetochlor (G+D+Ac), glyphosate plus alternating dicamba between years plus residual herbicides (diuron plus pendimethalin) (G+½D+R) and glyphosate plus alternating dicamba between years (G+½D), throughout the months after experiment initiation (Months). Error bars represent standard error of the mean for each data point.
This study shows the effectiveness of dicamba in decreasing dense populations of GR A. palmeri over the course of several years. Furthermore, these data are consistent with other research that showed a larger increase in the frequency of glyphosate resistance when glyphosate is used alone (Culpepper et al. Reference Culpepper, Grey, Vencill, Kichler, Webster, Brown, York, Davis and Hanna2006; Shergill et al. Reference Shergill, Bish, Biggs and Bradley2018). As observed in the first 4 yr of this study, weed seed production is critical for rapid evolution of GR A. palmeri. To date, research has shown a lack of fitness penalties in GR A. palmeri compared with GS A. palmeri; therefore, there is little delay in the buildup of GR individuals (Jasieniuk et al. Reference Jasieniuk, Brule-Babel and Morrison1996), and as the results of the present study suggest, the frequency of glyphosate resistance in a population is not easily reduced.
Although dicamba has been shown to be an effective tool in managing GR A. palmeri in cotton over 8 continuous years, diligent stewardship of dicamba should be a priority. Residual herbicides should be incorporated into herbicide programs, and postemergence-only programs must be avoided. Also, mixing two different herbicide MOAs would delay herbicide-resistance evolution more than rotating MOAs (Powles et al. Reference Powles, Preston, Bryan and Jutsum1997). Dicamba resistance has been reported in A. palmeri after three generations of exposure to sublethal rates (Tehranchian et al. Reference Tehranchian, Norsworthy, Powles, Bararpour, Bagavathiannan, Barber and Scott2017). However, at the end of our study, there were no survivors from GR A. palmeri populations treated with a field rate (560 g ae ha−1) of dicamba, thus indicating that no tolerance to this herbicide was found among any GR A. palmeri populations. The present study differed from the Tehranchian et al. (Reference Tehranchian, Norsworthy, Powles, Bararpour, Bagavathiannan, Barber and Scott2017) study in its application of full rather than sublethal rates (560 g ae ha−1 vs. 140, 280, and 420 g ae ha−1). In the present research, very few, if any, weed escapes were allowed to grow and contribute to the soil weed seedbank. While we are trying to reduce the incidence of a problematic species (e.g., A. palmeri) in the long term we might end up favoring other species (e.g., E. indica and D. sanguinalis), which ultimately could evolve resistance or increase their weediness. Therefore, special attention must be paid to the rest of the species in the community throughout time. Diverse weed control and production practices are necessary not only to delay the evolution of resistance (Norsworthy et al. Reference Norsworthy, Ward, Shaw, Llewellyn, Nichols, Webster, Bradley, Frisvold, Powles, Burgos, Witt and Barrett2012; Vencill et al. Reference Vencill, Nichols, Webster, Soteres, Mallory-Smith, Burgos, Johnson and McClelland2012), but to maintain weed population levels below those that are economically damaging.
Population trajectories can be very different depending on the herbicide program, weed species, and interactions among them, all of which shape the community over time. The present study illustrates how intensifying herbicide use (number of applications and MOAs) did not result in the reduction of the populations of all weed species present in the field. Thus, species with prolific reproduction can affect weed community assembly and “mask” the overall weed control potential of herbicide programs. The present study shows that intensive herbicide programs targeting dominant weed species can provide excellent control of those species but not necessarily of the entire weed community. In fact, the removal of dominant weed species seemed to have freed resources or reduced competitive interactions, favoring the increase of populations of other weed species. Ignoring those community dynamics can induce researchers to reach erroneous conclusions about the efficacy of weed control strategies and especially of herbicide programs.
Acknowledgments
The first 4 years of the field experiment were partially supported by Monsanto. This research was supported by funding to RGL from USDA-NIFA grants 2017-6505-26807, 2018-70006-28933, 2019-68012-29818, and Hatch Project NC02653. The authors declare no conflicts of interest.
Open access
