Glyphosate’s impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching
© Doğramacı et al.; licensee BioMed Central. 2015
Received: 6 June 2014
Accepted: 11 May 2015
Published: 19 May 2015
Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.
Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.
Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.
After commercialization of glyphosate [N-(phosphonomethyl)glycine] in 1974, it has become one of the most common broad spectrum herbicides used in agricultural systems . However, the evolution of glyphosate-resistant weeds has increased the cost of weed management, reduced yield potential, and mandated the need for new weed management strategies . No evidence of naturally occurring glyphosate-resistant plants was reported before the introduction of glyphosate-resistant transgenic crops, although several reports did indicate that biotypes of field bindweed (Convolvulus arvensis) and bermudagrass (Cynodon dactylon) were somewhat tolerant to glyphosate [3, 4]. Laticiferous perennial weeds such as hemp dogbane (Apocynum cannabinum), common milkweed (Asclepias syriaca), and leafy spurge (Euphorbia esula) have also been reported as glyphosate tolerant weeds [5–7]. Regardless, control of invasive perennial weeds in non-cultivated ecosystems of North America, including leafy spurge, generally requires integrated pest management programs that incorporate the use of glyphosate .
Reproduction of leafy spurge occurs through both seeds and vegetative propagules; however, persistence is mainly attributed to vegetative reproduction from an abundance of underground adventitious buds (UABs) . Dormancy in UABs of leafy spurge also contributes to escape from conventional control measures and long-term management often requires follow up application with herbicides. Because glyphosate is a non-selective herbicide and impacts desired plants as well, in long-term management programs glyphosate is often used at lower than recommended rates or in combination with other herbicides. In fact, recommended rates of glyphosate for pastures, rangelands, roadsides, riverbanks and recreational areas (~1 kg ha−1) destroy the aboveground shoots of leafy spurge plants; however, it does little or no damage to UABs, and leafy spurge regenerates from these UABs. However, greater concentrations of glyphosate (~2-6 kg ha−1) generally cause sub-lethal effects in underground root system of leafy spurge; these higher rates still do not destroy the plants completely, but induce uncontrolled growth patterns (e.g., tillering and branching) from UABs in following generations, similar to the glyphosate-induced tillering observed in quackgrass (Agropyron repense) and common milkweed [10, 11]. Further, the effects of glyphosate on leafy spurge  were proposed to involve glyphosate’s effect on release of apical dominance in lateral buds similar to that observed in sorghum (Sorghum bicolor), soybean (Glycine max), and pea (Pisum sativum), which have been linked to changes in auxin levels and transport [13, 14].
Translocation of glyphosate to the apical meristems, root meristems, and underground reproductive organs of perennial plants  inhibits nuclear encoded and chloroplast localized  5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) of the shikimate biosynthetic pathway [17, 18]. Glyphosate’s known mode of action occurs by inhibiting the EPSPS-catalyzed conversion of shikimate-3-phosphate (S3P) to 5-enolpyruvylshikimate 3-phosphate (EPSP) by forming an EPSPS–S3P–glyphosate complex that competes for phosphoenolpyruvate binding [19–21]. These are critical steps for the production of chorismate, which is a precursor of aromatic amino acids (phenylalanine, tyrosine and tryptophan), auxin, and many other secondary products essential for plant growth and development. Thus, the glyphosate-induced inhibition of EPSPS reduces chorismate production leading to increased accumulation of shikimate . However, beyond the known mode of action, the impact of glyphosate-treatment on transcriptional regulation has mainly focused on leaf tissues exposed to glyphosate in annual [23–25] and perennial grass species [26, 27]. Sub-lethal concentrations of foliar applied glyphosate also results in new shoot growth from UABs that produce a phenotype with reduced main stem elongation and increased branching from axillary bud outgrowth . Thus, investigating glyphosate’s impact on molecular processes and phytohormones affecting vegetative regrowth from UABs could provide new insights for novel weed management solutions.
The impact of foliar applied glyphosate has already provided some insights into the involvement of molecular processes associated with auxin, ethylene, and gibberellic acid (GA) biosynthesis and signaling pathways in crown buds of leafy spurge . The most significant changes observed in these UABs were for ENT-COPALYL DIPHOSPHATE SYNTHETASE 1, which is involved in a committed step for GA biosynthesis, and auxin transporters including PIN-FORMED (PINs), PIN-LIKES (PILS) and ATP-BINDING CASSETTE (ABC) TRANSPORTERS. Further, most of the significantly affected processes were associated with the endoplasmic reticulum, suggesting that this organelle plays some role in cellular physiology in response to glyphosate-treatment in leafy spurge. However, since Doğramacı et al.  only examined a select set of transcripts from UABs, it may not reflect the effects of foliar applied glyphosate on global gene expression or actual changes in phytohormone and shikimate levels. Thus, to gain further insights into glyphosate’s impact on molecular processes associated with phenotypes showing reduced stem elongation and increased branching due to axillary bud outgrowth, our objectives were 1) using RNAseq to determine global transcriptome abundance in aerial tissues derived from crown buds of glyphosate-treated and -untreated leafy spurge, 2) obtain shikimate and phytohormone concentration profiles from the same samples, and 3) elucidate relationships between transcriptome, shikimate, and phytohormone profiles to identify molecular processes that are altered by glyphosate.
Vegetative growth and phenotypic changes in response to glyphosate
Transcriptome changes in aerial tissues
Transcriptome analyses (see Additional file 2) revealed that 12,918 transcripts had significant changes in abundance (FDR ≤ 0.05) and 6,239 of these transcripts had significant changes ≥ 2-fold in either direction; 3,632 decreased and 2,607 increased abundance in glyphosate-treated samples relative to controls. Of the 12,918 differentially expressed genes, 9,283 transcripts had similarity (BlastX hits at E values < 10E-3) to genes in the Arabidopsis database (www.arabidopsis.org). Therefore, leafy spurge transcripts discussed throughout this manuscript are based on the most homologous match to Arabidopsis (see Additional file 2). Quantitative real-time PCR (qRT-PCR) based on a select set of genes (see Additional file 3) was used to validate the reliability of the RNAseq data, which revealed a similar trend in expression patterns - with the exception of CKX1. Additionally, Pearson’s (P-value = 5.1E-05) or Spearman’s (P-value = 0.012) correlation coefficients confirmed significant correlation between RNAseq and qRT-PCR data and validated use of the global transcript dataset.
Shikimate abundance and chorismate biosynthesis
Changes in abundance of transcripts involved in chorismate, auxin and strigolactone biosynthesis/signaling
Gene ID Abv.
EMB1144, EMBRYO DEFECTIVE 1144
SHIKIMATE KINASE 1
TRYPTOPHAN SYNTHASE BETA-SUBUNIT 2
TRYPTOPHAN SYNTHASE BETA TYPE 2
CYTOCHROME P450, FAMILY 79
CYTOCHROME P450, FAMILY 83
TRYPTOPHAN AMINOTRANSFERASE ARABIDOPSIS1
AUXIN RESPONSE FACTOR 2
AUXIN RESPONSE FACTOR 5
AUXIN RESPONSE FACTOR 6
AUXIN RESPONSE FACTOR 16
AUXIN-RESPONSIVE FAMILY PROTEIN
INDOLE-3-ACETIC ACID INDUCIBLE 19
INDOLE-3-ACETIC ACID INDUCIBLE 26
INDOLE-3-ACETIC ACID INDUCIBLE 29
LATERAL ROOT PRIMORDIUM 1
SMALL AUXIN UPREGULATED RNA 32
SMALL AUXIN UPREGULATED RNA 55
SMALL AUXIN UPREGULATED RNA 59
LIKE AUXIN RESISTANT 2
LIKE AUXIN RESISTANT 3
ATP-BINDING CASSETTE B19
IAA-ALANINE RESISTANT 3
IAA-LEUCINE RESISTANT (ILR)-LIKE GENE 6
MORE AXILLARY BRANCHES 1/CYTOCHROME P450
MORE AXILLARY BRANCHES 3
MORE AXILLARY BRANCHES 2
Glyphosate-treatment impacts molecular mechanisms regulating the shoot apical meristem
Sub-network enrichment analysis (SNEA) identified numerous over-represented ontologies among transcripts with increased or decreased abundance (see Additional file 4) in aerial shoots derived from crown buds of glyphosate-treated plants. An overview of the SNEA data indicated that hubs associated with increased transcript abundance are involved with numerous processes (e.g., cell cycle, cell differentiation and proliferation, cytokinesis, DNA replication, organ formation), proteins (e.g., ASYMMETRIC LEAVES1/2, CLAVATA3, CUP-SHAPED COTYLEDON1/2, E2F TRANSCRIPTION FACTOR, KNOTTED-LIKE FROM ARABIDOPSIS THALIANA1, SHOOT MERISTEMLESS, WUSCHEL), and small molecules (i.e., brassinolide, cytokinin) associated with shoot apical meristem (SAM) maintenance. However, with the exception of the process ‘stem growth’ (see Additional file 4), most hubs associated with decreased transcript abundance are involved in stress responses (e.g., apoptosis, photoinhibition, protein degradation, ROS generation, and virulence and xenobiotic clearance).
Phytohormone profiling and abundance of associated transcripts
Many of the aforementioned processes involved in SAM maintenance are regulated by phytohormone signaling and cross-talk [31, 32]. Indeed, previous studies  confirmed that pathways associated with phytohormone biosynthesis or signaling are impacted in crown buds of glyphosate-treated plants. To obtain a better understanding for the phenotypes observed in this study (Fig. 1), which could not be attributed to shikimate levels (Fig. 2) or transcripts involved in chorismate biosynthesis (Table 1), we quantified the abundance of phytohormones (see Additional file 5) in the vegetative shoots derived from crown buds of glyphosate-treated and control plants and their correlation to global transcriptome abundance.
Auxin abundance and associated transcripts
Tryptophan-dependent IAA biosynthesis is proposed to diverge into four pathways in plants [33, 34], which include the indole-3-acetamide (IAM), the indole-3-pyruvic acid (IPA), the tryptamine (TAM), and the indole-3-acetaldoxime (IAOX) pathways (Fig. 3). In subsequent-generations of aerial shoots obtained from crown buds of glyphosate-treated plants, transcripts homologous to Arabidopsis genes involved in the IPA pathway had increased abundance (TAA1, YUC4, YUC9; Table 1), whereas transcripts involved in the IAOX (CYP79B2, CYP83B1; Table 1) pathway had decreased abundance and no significant changes were observed for transcripts in the TAM and IAM pathways. Although overexpressing CYP79B2 in Arabidopsis results in increased levels of free auxin , in this study, auxin levels were increased (Fig. 3) even though transcript abundance for a leafy spurge homolog of CYP79B2 (of the IAOX pathway) was significantly decreased (−4.59 log2-fold; Table 1). Thus, the results presented in this study indicate that auxin is likely synthesized through the IPA pathway in subsequent generations of aerial tissues derived from plants exposed to glyphosate.
The increased abundance of auxin and/or its conjugates in response to glyphosate correlates well with the significant changes observed for transcripts homologous to Arabidopsis genes involved in auxin responses, which have either increased (ARF5, ARF6, GH3.1, IAA19, IAA26, LRP1) or decreased (ARF2, ARF16, ARFP, IAA29, SAUR32, SAUR55, SAUR59) abundance (Table 1). Other homologs of interest included several auxin transporters (ABCB19, LAX2, LAX3, PID, PIN1), which have a significant increase in transcript abundance. Although strigolactone abundance was not measured in this study, several homologs of transcripts involved in strigolactone biosynthesis (MAX1 and MAX3) were increased; whereas, abundance of transcript involved in downstream strigolactone signaling (MAX2) was decreased. Strigolactones can inhibit branching by depleting auxin transporters such as PIN1 [31, 32]. Thus, assuming these transcripts perform similar function in leafy spurge, the increased abundance of a transcript homologous to Arabidopsis PIN1 (Table 1), along with the increased branching observed in response to glyphosate treatment (Fig. 1) suggest that strigolactone biosynthesis or signaling are negatively impacted.
Cytokinins and cell division
Changes in abundance of transcripts involved in cytokinin biosynthesis/response and cell cycle
Gene ID Abv.
ADENOSINE KINASE 2
ADENINE PHOSPHORIBOSYL TRANSFERASE 2
ADENINE PHOSPHORIBOSYL TRANSFERASE 4
LONELY GUY 1
LONELY GUY 5
LONELY GUY 5
LONELY GUY 7
CYTOKININ OXIDASE 3
UDP-GLUCOSYL TRANSFERASE 73C1
URIDINE DIPHOSPHATE GLYCOSYLTRANSFERASE 74E2
HISTIDINE KINASE 5
HISTIDINE-CONTAINING PHOSPHOTRANSMITTER 1
ARABIDOPSIS RESPONSE REGULATOR 3
CELL DIVISION CONTROL 2
CYCLIN-DEPENDENT KINASE B1;2
CYCLIN-DEPENDENT KINASE B2;2
CYCLIN B 1;2
CYCLIN-DEPENDENT PROTEIN KINASE 3;2
GA abundance and associated transcripts
Changes in abundance of transcripts involved in GA biosynthesis, catabolism and signaling
Gene ID Abv.
GIBBERELLIN 3-OXIDASE 1
ENT-KAURENE OXIDASE 1
GIBBERELLIN 2-OXIDASE 1
GIBBERELLIN 2-OXIDASE 2
GIBBERELLIN 2-OXIDASE 4
GIBBERELLIN 2-OXIDASE 4
GA-STIMULATED ARABIDOPSIS 6
GA-STIMULATED ARABIDOPSIS 9
GA-STIMULATED ARABIDOPSIS 11
ABA abundance and associated transcripts
Changes in abundance of transcripts involved in ABA biosynthesis, catabolism and signaling
Gene ID Abv.
ALDEHYDE OXIDASE 4
ABA DEFICIENT 2
ABA DEFICIENT 2
ABA DEFICIENT 2
NINE-CIS-EPOXYCAROTENOID DIOXYGENASE 3
CYTOCHROME P450, FAMILY 707
UDP-GLUCOSYL TRANSFERASE 71B1
UDP-GLUCOSYL TRANSFERASE 71B1
ABSCISIC ACID RESPONSIVE ELEMENT-BINDING FAC 1
CALCIUM-DEPENDENT PROTEIN KINASE 1A
DEHYDRATION-RESPONSIVE ELEMENT-BINDING 1D
HOMEOBOX PROTEIN 6
PHOSPHOLIPASE D ALPHA 1
REGULATORY COMPONENTS OF ABA RECEPTOR 10
Other phytohormones and growth related processes
Changes in abundance of transcripts involved in ethylene and JA biosynthesis, metabolism and signaling
Gene ID Abv.
ACC SYNTHASE 1
ACC SYNTHASE 6
ACC SYNTHASE 10
ACC SYNTHASE 11
ETHYLENE RESPONSE FACTOR 1
ETHYLENE RESPONSIVE ELEMENT BINDING FAC 6
ERF DOMAIN PROTEIN 12
ETHYLENE-RESPONSIVE ELEMENT BINDING FAC 13
ETHYLENE-RESPONSIVE ELEMENT BINDING FAC 105
ETHYLENE RESPONSE FACTOR 110
RELATED TO AP2 3
ALLENE OXIDE SYNTHASE
CYTOCHROME P450, FAMILY 94
JASMONATE RESISTANT 1
WRKY DNA-BINDING PROTEIN 40
WRKY DNA-BINDING PROTEIN 50
WRKY DNA-BINDING PROTEIN 70
ARABIDOPSIS THALIANA SULFOTRANSFERASE 2A
Changes in transcript abundance of ABC Transporters and Glutathione Transferases (GSTs)
Gene ID Abv.
ATP-BINDING CASSETTE G2
ATP-BINDING CASSETTE G5
ATP-BINDING CASSETTE G6
ATP-BINDING CASSETTE G6
ATP-BINDING CASSETTE G8
ATP-BINDING CASSETTE G14
ATP-BINDING CASSETTE G25
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
ATP-BINDING CASSETTE G40
GLUTATHIONE S-TRANSFERASE 30
GLUTATHIONE S-TRANSFERASE F11
GLUTATHIONE TRANSFERASE LAMBDA 2
GLUTATHIONE S-TRANSFERASE TAU 6
GLUTATHIONE S-TRANSFERASE TAU 25
GLUTATHIONE S-TRANSFERASE TAU 25
GLUTATHIONE S-TRANSFERASE TAU 25
GLUTATHIONE S-TRANSFERASE THETA 1
GLUTATHIONE S-TRANSFERASE ZETA 1
Direct exposure of leafy spurge aerial tissues to sub-lethal concentrations of glyphosate causes an indirect effect that impacts molecular processes in UABs prior to growth induction by decapitation . In this study, reduced main stem elongation and increased branching was observed in four subsequent generations derived from UABs present in the same root system of treated plants (Fig. 1, and Additional file 1), and the effects of glyphosate on phenotype did not diminish in these generations. However, after four generations of vegetative growth no more UABs were left in the root system of these plants (grown in cone-tainers) to observe if the plants would ultimately outgrow the effects of glyphosate. Parallel studies performed under field conditions revealed that higher rates of foliar-applied glyphosate (3 or 6 kg ha−1) affects vegetative growth from UABs of leafy spurge and causes similar phenotypical alterations for at least two consecutive years after the glyphosate-treatment (Doğramaci et al. unpublished data). Sub-lethal rates of glyphosate have also been reported to induce similar phenotypes in annual species such as corn (Zea mays), cotton (Gossypium hirsutum), sorghum, and wheat (Triticum aestivum) [13, 41–43]. To gain insights on factors involved in these glyphosate-induced phenotypes, here we determined the effects of foliar applied glyphosate on the abundance of metabolites (shikimate and hormones) and their correlation to global transcript abundance in subsequent generation of aerial tissues.
Glyphosate-treatment caused differential accumulation of shikimate
Little information is available on the abundance of shikimate in underground buds or in subsequent generations derived from glyphosate-treated perennials, with the exception of glyphosate-induced accumulation of shikimate in sprouted tubers of the perennial weed purple nutsedge (Cyperus rotundus) . Overall, changes in shikimate abundance in leafy spurge tissues after glyphosate-treatment (Fig. 2) were relatively low compared to those reported in other species [45, 46], but it is important to note that the rate used in this study had only a sub-lethal effect. Minor changes in abundance of transcripts involved in the chorismate biosynthetic pathway in aerial tissues (Table 1) are similar to that previously observed in crown buds of leafy spurge  or other weedy species such as palmer amaranth and horseweed following glyphosate-treatment [46, 47]. The increased accumulation of shikimate in aerial tissues of leafy spurge directly treated with glyphosate or in crown buds that were indirectly exposed (Fig. 2) is not surprising, given that significant shikimate accumulation is detectable several days prior to visual injury in response to glyphosate-treatment in other species . Regardless of their seasonal dormancy phase (i.e., para-, endo-, or eco-dormancy), UABs of leafy spurge are active and rapidly respond to environmental changes; in fact, during the shifts between seasonal phases there is some limited growth and development  suggesting UABs may serve as a sink tissue. We hypothesize that the increased abundance of shikimate in crown buds is due to either 1) translocation of glyphosate from foliar tissues into the root system and UABs to locally inhibit EPSPS, and/or 2) translocation and compartmentalization of foliar shikimate into the root system or UABs. However, there is more support for the first option based on previous reports indicating that glyphosate or its metabolites are translocated into the roots and UABs of leafy spurge , similar to that reported in other weedy species [15, 50].
Although crown buds of foliar glyphosate-treated leafy spurge initially had increased shikimate abundance prior to decapitation, in subsequent generations of aerial shoots (six weeks post-decapitation) shikimate abundance decreased and was similar to that observed in untreated controls (Fig. 2). These results suggest that duration of time after growth induction either allowed for dissipation of shikimate originally observed in crown buds, or EPSPS is no longer inhibited. However, the reduction in shikimate levels did not reverse the glyphosate-induced phenotype (Fig. 1, and Additional file 1). Studies in Arabidopsis revealed that shikimate can be consumed as a secondary metabolite as an ester . Further, silencing of a gene (HYDROXYCINNAMOYL-COA SHIKIMATE/QUINATE HYDROXYCINNAMOYLTRANSFERASE) downstream of the chorismate pathway (i.e., phenylpropanoid metabolism) that involves shikimate recycling causes stunting  similar to that observed in this study, which appears to inhibit polar auxin transport . The results of this study (Fig. 2) are consistent with abundance of shikimate initially increasing after glyphosate-treatment in other plant systems and then decreasing over time [48, 53-55]. For example, abundance of shikimate generally peaks between 4 or 7 days after glyphosate-treatment in crops such as corn, soybean (Glycine max), sunflower (Helianthus annuus), and wheat [48, 53], and in several weedy species .
Even though glyphosate is known to be metabolized in plants and presence of glyphosate oxidoreductase-like enzymes (involved in the degradation pathway of glyphosate in soil) have been speculated, to date no corresponding enzymes or genes have been identified [55, 56]. Therefore, we are not able to determine the extent of glyphosate metabolism based on transcript abundance in leafy spurge. Previous studies  have shown that 14C-glyphosate levels remain similar in leaves and roots of field bindweed, Canada thistle (Cirsium arvense), and tall morning glory (Ipomoea purpurea) 3–4 weeks post-application, with only a limited accumulation of the main glyphosate metabolite aminomethylphosphonic acid (AMPA). However, in glyphosate-resistant soybean, the highest levels of glyphosate and AMPA were found in treated leaves within 1 day after glyphosate treatment and gradually decreased within 3 weeks . Further, while exogenous application of AMPA did not affect shikimate levels it had phytotoxic effects and induced some injuries (e.g., reduced chlorophyll content, shoot fresh weight) in glyphosate resistant and non-resistant soybean ; but, stunted or bushy phenotypes were not reported after AMPA application.
Because growth-induction of leafy spurge crown buds results in decreased abundance of starch and sucrose and increased abundance of fructose , it appears that storage reserves in the root system are remobilized into free sugars to provide an energy source for new vegetative shoot growth. It is plausible that decapitation-induced growth causes a similar process to occur where root sequestered glyphosate or its metabolites, if any, are remobilized into UABs. Remobilization of these substances could cause subsequent-generations of vegetative growth to continually be exposed to residual toxicity and impact SAM maintenance and axillary bud outgrowth. Glyphosate’s effects on meristematic tissues over extended periods in white birch (Betula papyrifera) has also been speculated to involve sugar flow into the root system in the fall and a reverse flow of sugar to developing leaves in the spring . This hypothesis is consistent with glyphosate being mobilized to sink tissues such as growing meristems via the phloem . Alternatively, it is possible that glyphosate has an irreversible effect on the main SAM of UABs prior to decapitation-induced growth. The effects of glyphosate on a select set of molecular processes in crown buds of foliar glyphosate-treated plants , and increased abundance of shikimate in crown buds 7 days post-treatment (Fig. 2), would support this hypothesis.
Auxin, cytokinin, strigolactone cross-talk and shoot branching
Shoot branching is known to involve the key hormones auxin, cytokinin, and strigolactones [31, 60]. The increased abundance of transcripts homologous to Arabidopsis genes involved in the IPA biosynthetic pathway leading to auxin (Fig. 3), along with the increased levels of auxin in subsequent generations of aerial tissues from glyphosate-treated plants suggest that tryptophan-dependent auxin biosynthesis is not inhibited at this stage of development. These results also fit well with the abundance of shikimate observed in the same aerial tissues (Fig. 2). Because increased abundance of shikimate is a marker for inhibition of EPSPS, the increased abundance of auxin and similar levels of shikimate in both glyphosate treated and control plants suggest that EPSPS is not inhibited in subsequent generations of aerial shoots (Fig. 2). These results may suggest either lack of glyphosate or inhibition of its transport into plastids in these shoots, therefore not affecting the chorismate pathway or its downstream auxin biosynthesis pathway. However, we cannot rule out the possibility that the increase in auxin could also be due to an increase in the number of growing meristems.
In this study, the bushy phenotypes derived from crown buds of foliar glyphosate-treated plants do not appear to have a dominant/leading apical meristem to suppress outgrowth of axillary buds (Fig. 1), suggesting that apical dominance is lacking in the system and/or the polar auxin transport stream is affected. Previous studies have shown that exogenous application of auxin represses stem elongation of shoots derived from crown buds of leafy spurge [57, 61], which may have some bearing on the reduced main stem elongation and increased branching in response to glyphosate (Fig. 1). Indeed, auxin has been proposed to regulate shoot branching through at least two processes [31, 32]. First, the canalization theory involves the effects of apical dominance on abundance and transport of auxin in stems, which is proposed to inhibit auxin transport from axillary buds to repress axillary bud outgrowth. The second proposed mechanism involves auxin’s impact on acropetal movement of cytokinins and strigolactone into the buds to regulate bud outgrowth [62–66]. However, whereas cytokinins induce bud outgrowth, strigolactone signaling can affect PIN1 polarization in the plasma membrane to inhibit auxin transport .
PINs encode auxin efflux carriers that mediate tissue-specific, cell-to-cell polar auxin transport  critical for developmental patterning and differential growth responses . Strigolactone biosynthesis or signaling mutants also have increased axillary bud outgrowth, resulting in a bushy and dwarfed phenotype in Arabidopsis ; increased branching is also correlated with high levels of PIN1 and auxin in the main stem [70, 71]. However, in this study, the increased transcript abundance for homologs of Arabidopsis MAX1 and MAX3 (Table 1), involved in the strigolactone biosynthesis , do not correlate well with the glyphosate-induced branching phenotype. Indeed, the minimal decrease in transcript abundance for a homolog of MAX2, involved in downstream strigolactone signaling through regulation of PIN1  could account for the increased transcript abundance for PIN1 (Table 1). Assuming that the leafy spurge homologs perform similar functions as their Arabidopsis counterparts, these results suggest that PIN1 could be abundant and affecting auxin transport within the polar auxin transport stream of leafy spurge plants in response to glyphosate, in a manner similar to that previously reported .
The increased abundance of bioactive cytokinins (Fig. 4) in the aerial shoots derived from crown buds of foliar glyphosate-treated plants is consistent with axillary bud outgrowth (Fig. 1). Ironically, although bioactive auxin levels were increased, auxin’s regulatory role on cytokinin levels and acropetal movement appears to be negated in glyphosate-treated plants. Possible explanations for the disparity of the auxin-cytokinin interaction include the local biosynthesis of cytokinin in aerial shoots, particularly since cytokinins are known to be synthesized in most plant tissues . Indeed, the increased abundance of transcripts (Table 2) homologous to Arabidopsis genes involved in cytokinin biosynthesis (AK2, APT2, APT4, IPT3, LOG5), and decreased abundance of transcripts involved in cytokinin metabolism (CKX3, UGT73C1, UGT74E2) in subsequent generations of aerial shoots in response to glyphosate is consistent with the abundance of cytokinins (Fig. 4). Additionally, cytokinins are also known to impact polar auxin transport through modulation of auxin efflux carrier activity . Interestingly, numerous transcripts homologous to Arabidopsis genes involved in auxin transport had a significant increase in abundance (ABCB19, LAX2, LAX3, PID, PIN1; Table 1). Some of these auxin transporter family members (e.g., PINs, and ABC TRANSPORTERS) also had significant changes in transcript abundance in crown buds of foliar glyphosate-treated plants prior to growth induction by decapitation .
Interestingly, sugar signaling, not auxin, has been identified as the main regulator of initial axillary bud outgrowth in pea . Although it is difficult to determine if a similar mechanism involving sugar signaling impacts axillary bud outgrowth in aerial shoots of leafy spurge in response to glyphosate-treatment, Gene Set Enrichment Analyses (GSEA, Additional file 7) did identify both decreased and increased abundance of transcripts associated with carbohydrate metabolism and transport. Although limiting sugars through translocation to the shoot tips of axillary buds is proposed to maintain apical dominance, long-term signaling involving auxin, cytokinin and strigolactone are still thought to play a significant role in branching .
Glyphosate induces transcript abundance of large gene families
Large gene families associated with detoxification processes and resistance to multiple herbicides [72–74] had significant changes in transcript abundance in aerial shoots derived from crown buds of glyphosate-treated plants. For example, members of the P450s (46 of 95 increased), UGTs (13 of 51 increased), GSTs (6 of 10 increased), and ABC TRANSPORTERS (8 of 15 increased) had significant changes in abundance (see Additional file 2). In other plant systems, P450s together with GSTs and UGTs are involved in herbicide biochemical modification through metabolism, whereas ABC TRANSPORTERS are involved in compartmentalizing herbicides and their metabolites; this four-step process is associated with herbicide resistance in some weeds [24, 72, 75-78]. As previously discussed, several homologs of Arabidopsis P450s, UGTs and ABC TRASPORTERs differentially expressed in this study are linked to hormone biosynthesis or signaling (Tables 1–5), and it is possible that other members of these families (see Additional file 2) could be playing some role in reducing the phytotoxic effects of residual glyphosate. Assuming these leafy spurge homologs perform similar functions, significant changes in abundance of these gene families suggest that detoxification processes could play some direct or indirect role in reducing the phytotoxic effects of glyphosate. Alternatively, because glyphosate can rapidly disrupt plastids, chloroplasts, and the endoplasmic reticulum in plants [79, 80], and some leafy spurge homologs of P450s (see Additional file 2) are known to be localized in these organelles in other systems , it is also possible that residual glyphosate impacts cellular organelles in UABs and could help explain the differential abundance of P450s observed in this study.
Plants treated with glyphosate have reduced levels of bioactive GA
Even though bioactive GAs were not detectable in any samples (Fig. 5), the results indicate that precursors of bioactive GAs were increased, but the abundance of bioactive GAs may be reduced in subsequent generation of aerial tissues derived from crown buds of glyphosate-treated plants. Further, the decreased abundance of a transcript (GA3OX1) involved in synthesis of bioactive GAs and increased abundance of transcripts (GA2OX2, GA2OX4) involved in catabolism of bioactive GAs (Table 3, and Fig. 5), support the decreased abundance of bioactive GAs observed in glyphosate-treated plants. Overall, these results indicate that downstream GA biosynthesis and signaling are negatively impacted, which could contribute to inhibition of stem elongation of tissues derived from crown buds of foliar glyphosate-treated leafy spurge. This hypothesis is supported by the fact that certain commercial chemicals used to stunt growth (e.g., phosphon D, cycocel, ancymidol and paclobutrazol) do so by inhibiting GA biosynthesis .
Translocation of foliar applied glyphosate and its sequestration into the root system of leafy spurge  caused irreversible effects to the main SAM of UABs prior to growth-induction by decapitation of the directly treated aerial tissues. This theory is supported by the increased abundance of shikimate in crown buds post foliar glyphosate-treatment and prior to growth-induction (Fig. 2). After growth-induction by decapitation, this irreversible effect, in turn, causes the main stem of subsequent generations of aerial tissue to lose apical dominance, thus impacting axillary bud outgrowth. Alternatively, the potential remobilization of root-sequestered glyphosate would be consistent with apical meristems of growing shoots acting as sink tissues [15, 83], which is supported by the fact that SNEA (see Additional file 4) highlighted numerous molecular processes associated with SAM maintenance containing transcripts with increased abundance (6 weeks post-decapitation). Collectively, these two hypotheses provide insights into potential mechanisms resulting in glyphosate-induced phenotypes that have reduced main stem elongation and increased branching due to axillary bud outgrowth (Fig. 1).
The results of this study are consistent with the known involvement of phytohormone cross-talk to regulate branching. The increased abundance of bioactive auxin and cytokinin (Fig. 2) and PIN1 (Table 1) in subsequent generations of aerial shoots (6 weeks post-decapitation) in response to glyphosate-treatment is consistent with axillary bud outgrowth, assuming auxin is being exported from buds into the polar auxin transport stream [70, 71]. Additionally, the increased abundance of cytokinins in subsequent generations of aerial tissues derived from UABs of glyphosate-treated plants would also be consistent with the increased abundance of transcripts involved in cytokinin biosynthesis (Fig. 4) or associated with cell division processes (Table 2). The increased abundance for homologs of MAX1 and MAX3 suggest that strigolactone biosynthesis may be activated in response to the residual effects of glyphosate. However, in our system, the increased abundance of a PIN1 homolog and branching in response to glyphosate are counterintuitive to increased strigolactone signaling, which can inhibit branching through regulation of auxin transporters such as PIN1 [31, 32, 69].
The differential abundance of transcripts for P450s, known to be localized in plastids, endoplasmic reticulum, and mitochondria, support glyphosate’s impact on disruption of these intracellular organelles, as demonstrated in other plant systems [79, 80]. Further, the differential expression of numerous P450s, GSTs, UGTs, and ABC TRANSPORTERS in subsequent generations of aerial shoots in response to glyphosate (Additional file 2) could suggest the occurrence of detoxification processes through metabolism and/or compartmentalization in these tissues. The fact that homologs of ABC TRANSPORTERS, PINs, and PILS also had differential abundance in crown buds of leafy spurge in response to glyphosate  add further support for molecular processes associated with membrane transport and/or sequestration being impacted by glyphosate-treatment prior to growth-inducing decapitation. Although further research will be needed to determine the exact mechanism leading to these glyphosate-induced phenotypes, this research provides insights for molecular processes associated with stunted growth and increased branching from underground adventitious buds of leafy spurge.
Material and Methods
Plant material, glyphosate-treatments, and vegetative growth
Leafy spurge plants were propagated from a uniform biotype (1984-ND001) in cone-tainers and maintained in a greenhouse as described by Anderson and Davis . After conducting dose–response experiments using concentration of 0 to 6.72 kg ha−1 glyphosate as active ingredient, aerial tissues of four-month old greenhouse grown leafy spurge plants were treated with 0 or 2.24 kg ha−1 glyphosate plus 0.25 % surfactant as described by Doğramacı et al. . Plants were returned to growth-conducive conditions in a greenhouse after foliar glyphosate-treatment. Each experiment included four biological replicates; aerial tissues of eight plants from each replicate were decapitated 7 days post-application to determine the vegetative growth rate from crown buds as previously described by Foley et al. . After six weeks of vegetative growth from crown buds of foliar glyphosate-treated and control plants, shoot heights were measured and analyzed using the generalized linear mixed model (PROC GLIMMIX) procedure of SAS 9.2. Further, these same cone-tainers were used to determine the longevity of glyphosate’s effects on vegetative growth from UABs of glyphosate-treated and control plants; every six weeks aboveground shoots were decapitated for three additional re-growth cycle. However, only the first generation of aerial tissues derived from UABs of glyphosate-treated and control plants were harvested and flash frozen in liquid nitrogen and stored in −80 °C for further experiments i.e., shikimate and phytohormone quantification, and transcript analyses by RNAseq and qRT-PCR. Additionally, aerial tissues of leafy spurge plants, which were directly exposed to the treatments and their crown buds were collected for shikimate quantification 7 days after the treatments.
Quantification of shikimate levels
Quantification of shikimate in aerial tissues directly exposed to the treatments (7 days post-treatment), crown buds of foliar treated plants (7 days post-treatment), and aerial tissues derived from crown buds of foliar treated plants (six weeks after growth-inducing decapitation) were accomplished spectrophotometrically [45, 86]. In brief, 0.5 gram of frozen tissue was extracted in 2 mL 0.25 N HCl at room temperature. 250 μL of this extract was added into one ml 0.25 % (w/v) periodic acid/0.25 % (w/v) meta-periodate and incubated at 37 °C for 15 minutes. 1 mL of 0.6 M NaOH/0.22 M Na2SO3 was added into the sample and read by a spectrophotometer (Beckman, DU 7400); 510 μM shikimate (Sigma) was used as the standard. These assays included four biological and two technical replicates; technical replicates were averaged prior to calculating mean of biological replicates and standard error.
Quantification of phytohormones (ABA and ABA metabolites, auxins, cytokinins, and gibberellins) in the eight aerial tissue samples (2 treatments [controls and 2.24 kg ha−1 glyphosate] × 4 biological replicates) derived from crown buds of foliar glyphosate-treated plants was conducted at the NRCC-PBI (Saskatoon, SK, Canada) by UPLC-ESI-MS/MS. Deuterated forms of each of the hormones were used as internal standards [87, 88]. Extraction and purification of samples and quantification of phytohormones was performed following a procedure described in Chiwocha et al. [89, 90]. Profiling data obtained from four biological replicates were used to calculate the mean and standard error. The list and abbreviations of hormone profiles are provided in Additional file 5.
RNA extraction, cDNA library preparation, and transcript analyses by RNAseq
Total RNA was purified from the eight aerial tissue samples derived from crown buds of foliar glyphosate-treated plants according to the pine tree RNA extraction protocol . The RNA was treated with amplification grade DNase1 (Invitrogen), quantified by the Qubit® 2.0 Fluorometer, and quality was confirmed by agarose gel electrophoresis. The RNAseq libraries were prepared with the TruSeq Stranded mRNA LT Sample Preparation Kit (Illumina, Cat. N° RS–122–2101) starting from 1 μg of total RNA. The pool of barcoded RNAseq libraries was quantified by qRT-PCR using the Library Quantification kit (Kapa Biosystems, Cat. N° KK4824). The size range of the final cDNA libraries was determined on an Agilent bioanalyzer DNA7500 DNA chip (Agilent Technologies). The cDNA libraries were sequenced on one lane for 151 cycles from each end of the cDNA fragments on a HiSeq2500 using a TruSeq SBS sequencing kit version 1 (Illumina). The sequence images were transformed with the Real Time Analysis 184.108.40.206 Illumina software to bcl files, which were demultiplexed to fastq files with CASAVA version 1.8.2. The quality-scores line in fastq files processed with Casava1.8.2 use an ASCII offset of 33 known as Sanger scores.
Transcriptome assembly, annotation and expression analyses
The Trinity package version trinityrnaseq-r2013-02-25  was utilized to assemble the transcriptome from the RNASeq data. Prior to assembly, adaptor sequences were trimmed and so were any low quality bases (<30 phred score) from both ends of the sequence. Sequences of 50 nucleotide and higher were utilized for the Trinity assembly. The assembly was performed as per guidelines in the published protocols , and on the Trinity website (http://trinityrnaseq.sourceforge.net/#running_trinity). The abundance of every gene (component) in each of the eight samples was determined by performing alignment and abundance estimation using the RSEM package  (version 1.2.3). Only those genes that had at least one sample with the transcripts per million (TPM) value of ≥ 0.5 were used as input for the next step. Differential gene expression was performed at the level of genes with the EBSeq package  (version 1.0) using the median normalization approach and 10 iterations of the algorithm used by EBSeq. Genes were called differentially expressed if they had an FDR value of ≤0.05.
Annotation of the assembled transcriptome and predicted proteome was performed with the Trinotate package (version trinotate_r20130826). Transdecoder, a Perl script packaged with Trinity, was used to predict proteins, yielding the most likely longest-ORF peptide candidates file best_candidates.eclipsed_orfs_removed.pep. Analyses were performed as detailed in the Trinotate manual (http://trinotate.sourceforge.net/), with the following noted exceptions. BLASTx of assembled transcripts and BLASTp most likely longest-ORF peptide candidates against UniProt-SwissProt database were run with the options “-evalue 10E-3 –max_target_seqs 5” . For more reliable functional annotation, BLASTx of Trinity.fasta and BLASTp of best_candidates.eclipsed_orfs_removed.pep against Arabidopsis peptides were performed (Athaliana_167_protein.fa from Phytozome v9.0). BLAST output files were parsed to 1 top hit per query with an in-house Perl script prior to insertion in the Trinotate SQLite database (http://trinotate.sourceforge.net/#LoadingSQLite). After extraction of the Trinotate annotation report from the database (http://trinotate.sourceforge.net/#OutputReport), Arabidopsis annotations were added using an in-house Perl script. The expression data is deposited at Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) as GEO dataset query GSE56509.
qRT-PCR was utilized to confirm transcript profiles obtained from RNAseq for selected genes. Sequences from a leafy spurge EST database  were used for designing primer pairs using the Primer-Select of Lasergene 8 software program (DNASTAR, Inc., Madison, WI). These primer pairs were used to quantify abundance of selected transcripts (see Additional file 3). Total RNA samples were reverse transcribed into complementary DNA (cDNA) as described in Doğramacı et al. . In brief, 5 μg of total RNA was treated with DNase1 amplification grade (Invitrogen), and reverse transcribed using Super Script First-Strand Synthesis Kit3 (Invitrogen) in a 20 μl volume according to the manufacturer’s instructions. After cDNA synthesis, each reaction was diluted to 800 μl total volume and stored at −80 °C. For qRT-PCR reactions, 1 μl total cDNA was added to a 10 μl PCR reaction mixture containing 5 μl of LightCycler® 480 SYBR Green I Master and 0.5 μl of each primer-set. Transcript abundance was measured from three technical and four biological replicates using a LightCycler 480 II (Roche). All transcript values were normalized using the reference genes ARF2, PTB, SAND, and ORE9 identified by Chao et al. . QbasePLUS version 2.4 software (Biogazelle, Ghent, Belgium) was used to normalize expression values. Values from four biological replicates were averaged and data from control samples (0 kg ha−1) were used for baseline expression. Gene abbreviations and descriptions of all putative homologous leafy spurge genes included throughout this report were obtained from an Arabidopsis website (www.arabidopsis.org) and are presented in (see Additional file 2). Pathways in Figs. 2–6 were developed based on the PlantCyc website and/or available literature on hormone biosynthesis/signaling, and more details can be found at http://pmn.plantcyc.org/PLANT/class-tree?object=Pathways.
SNEA and GSEA of transcriptome data
The entire transcriptome dataset (i.e., transcripts per million values; see Additional file 2) was analyzed using Ariadne Pathway Studio 9.0 Software-Resnet Plant Version 2.1 (Ariadne Genomics Inc., Rockville, MD, USA) to obtain SNEA and GSEA using Mann–Whitney U test algorithm with a 0.05 enrichment p value cutoff. Central hubs of sub-networks were obtained by SNEA (Additional file 4) using the customized advanced parameters to identify the most important biological processes, transcription factors, binding partners, miRNA targets, and protein modification targets for the transcriptome data. GSEA was used to identify predefined sets of genes over-represented between treatments (Additional file 7) based on gene ontology (http://www.geneontology.org/) or AraCyc metabolic pathways (http://pmn.plantcyc.org/).
The authors thank Brant B. Bigger, Leonard W. Cook and Cheryl A. Huckle for their technical assistance, and Dr. Suzanne Abrams and Dr. L. Irina Zaharia for collaboration on phytohormone profiling. We also thank Dr. Dale L. Shaner for providing the protocol for quantifying shikimate abundance.
Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer.
- Duke SO, Powels SB. Glyphosate: a once-in-a-century herbicide. Pest Manag Sci. 2008;64:319–25.View ArticlePubMedGoogle Scholar
- Duke SO, Powles SB. Glyphosate-resistant crops and weeds: Now and in the future. AgBioforum. 2009;12:346–57.Google Scholar
- Bryson CT, Wills GD. Susceptibility of bermudagrass (Cynodon dactylon) biotypes to several herbicides. Weed Sci. 1985;33:848–52.Google Scholar
- Westwood JH, Tominaga T, Weller SC. Characterization and breakdown of self-incompatibility in field bindweed (Convolvulus arvensis L.). J Hered. 1997;88:459–65.View ArticleGoogle Scholar
- Gottrup O, O'Sullivan PA, Schraa RJ, VandenBorn WH. Uptake, translocation, metabolism and selectivity of glyphosate in Canada thistle and leafy spurge. Weed Res. 1976;16:197–201.View ArticleGoogle Scholar
- Wyrill III JB, Burnside OC. Absorption, translocation, and metabolism of 2,4-D and glyphosate in common milkweed and hemp dogbane. Weed Sci. 1976;24:557–66.Google Scholar
- Cramer GL, Burnside OC. Control of common milkweed (Asclepias syriaca). Weed Sci. 1981;29:636–40.Google Scholar
- Lym RG. Leafy spurge (Euphorbia esula) control with glyphosate plus 2,4-D. J Range Manage. 2000;53:68–72.View ArticleGoogle Scholar
- Anderson JV, Horvath DP, Chao WS, Foley ME: Bud dormancy in perennial plants: A mechanism for survival. In: Lubzens E, Cerda J and Clark M, editors. Dormancy and Resistance in Harsh Environments. Topics in Current Genetics vol 21. Springer-Verlag Berlin Heidelberg; 2010. p. 69–90.Google Scholar
- Coupland D, Caseley JC. Reduction of silica and increase in tillering induced in Agropyron repens by glyphosate. J Exp Bot. 1975;26:138–44.View ArticleGoogle Scholar
- Waldecker MA, Wyse DL. Chemical and physical effects of the accumulation of glyphosate in common milkweed (Asclepias syriaca) root buds. Weed Sci. 1985;33:605–11.Google Scholar
- Maxwell BD, Foley ME, Fay PK. The influence of glyphosate on bud dormancy in leafy spurge (Euphorbia esula). Weed Sci. 1987;35:6–10.Google Scholar
- Baur JR. Reduction of glyphosate induced tillering in sorghum (Sorghum bicolor) by several chemicals. Weed Sci. 1979;27:69–73.Google Scholar
- Lee TT. Release of lateral buds from apical dominance by glyphosate in soybean and pea seedlings. J Plant Growth Regul. 1984;3:227–35.View ArticleGoogle Scholar
- Shaner DL. Role of translocation as a mechanism of resistance to glyphosate. Weed Sci. 2009;57:118–23.View ArticleGoogle Scholar
- Della-Cioppa G, Bauer SC, Klein BK, Shah DM, Fraley RT, Kishore GM. Translocation of the precursor of 5-enolpyruvulshikimate-3-phosphate synthase into chloroplasts of higher plants in vitro. Proc Natl Acad Sci U S A. 1986;83:6873–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Amrhein N, Deus B, Gehrke P, Steinrucken HC. The site of the inhibition of the shikimate pathway by glyphosate. II. Interference of glyphosate with chorismate formation in vivo and in vitro. Plant Physiol. 1980;66:830–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Steinrucken HC, Amrhein N. The herbicide glyphosate is a potent inhibitor of 5-enolpyruvylshikimic acid-3-phospahte synthase. Biochem Biophys Res Comm. 1980;94:1207–12.View ArticlePubMedGoogle Scholar
- Boocock MR, Coggins JR. Kinetics of 5-enolpyruvylshikimate-3-phosphate synthase inhibition by glyphosate. FEBS Lett. 1983;154:127–33.View ArticlePubMedGoogle Scholar
- Gruys KJ, Walker MC, Sikorski JA. Substrate synergism and the steady-state kinetic reaction mechanism for EPSP synthase from Escherichia coli. Biochemistry. 1992;31:5534–44.View ArticlePubMedGoogle Scholar
- Franz JE, Mao MK, Sikorski JA. Glyphosate: a unique global herbicide. Am Chem Soc Monogr. 1997;189:521–615.Google Scholar
- Herrmann KM, Weaver LM. The shikimate pathway. Annu Rev Plant Physiol Plant Mol Biol. 1999;50:473–503.View ArticlePubMedGoogle Scholar
- Zhu J, Patzoldt WL, Shealy RT, Vodkin LO, Clough SJ, Tranel PJ. Transcriptome response to glyphosate in sensitive and resistant soybean. J Agric Food Chem. 2008;56:6355–63.View ArticlePubMedGoogle Scholar
- Peng Y, Abercrombie LL, Yuan JS, Riggins CW, Sammons RD, Tranel PJ, et al. Characterization of the horseweed (Conyza canadensis) transcriptome using GS-FLX 454 pyrosequencing and its application for expression analysis of candidate non-target herbicide resistance genes. Pest Manag Sci. 2010;66:1053–62.View ArticlePubMedGoogle Scholar
- Yuan JS, Abercrombie LLG, Cao Y, Halfhill MD, Zhou X, Peng Y, et al. Functional genomics analysis of horseweed (Conyza canadensis) with special reference to the evolution of non target-site glyphosate resistance. Weed Sci. 2010;58:109–17.View ArticleGoogle Scholar
- Cebeci O, Budak H. Global expression patterns of three Festuca species exposed to different doses of glyphosate using the Affymetrix GeneChip wheat genome array. Comp Funct Genomics. 2009;2009:505701.View ArticlePubMed CentralGoogle Scholar
- Unver T, Bakar M, Shearman RC, Budak H. Genome-wide profiling and analysis of Festuca arundinacea miRNAs and transcriptomes in response to foliar glyphosate application. Mol Genet Genom. 2010;283:397–413.View ArticleGoogle Scholar
- Doğramacı M, Anderson JV, Chao WS, Foley ME. Foliar application of glyphosate affects molecular mechanisms in underground adventitious buds of leafy spurge (Euphorbia esula) and alters their vegetative growth patterns. Weed Sci. 2014;62:217–29.View ArticleGoogle Scholar
- Lang GA. Dormancy: A new universal terminology. Hort Sci. 1987;22:817–20.Google Scholar
- Horvath DP. Role of mature leaves in inhibition of root bud growth in Euphorbia esula L. Weed Sci. 1999;47:544–50.Google Scholar
- Vanstraelen M, Benková E. Hormonal interactions in the regulation of plant development. Annu Rev Cell Dev Biol. 2012;28:463–87.View ArticlePubMedGoogle Scholar
- Shinohara N, Taylor C, Leyser O. Strigolactone can promote or inhibit shoot branching by triggering rapid depletion of the auxin efflux protein PIN1 from the plasma membrane. PLoS Biol. 2013;11, e1001474.View ArticlePubMed CentralPubMedGoogle Scholar
- Mashiguchi K, Tanaka K, Sakai T, Sugawara S, Kawaide H, Natsume M, et al. The main auxin biosynthesis pathway in Arabidopsis. Proc Natl Acad Sci U S A. 2011;108:18512–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Mano Y, Nemoto K. The pathway of auxin biosynthesis in plants. J Exp Bot. 2012;63:2853–72.View ArticlePubMedGoogle Scholar
- Zhao Y, Hull AK, Gupta NR, Goss KA, Alonso J, Ecker JR, et al. Trp-dependent auxin biosynthesis in Arabidopsis: involvement of cytochrome P450s CYP79B2 and CYP79B3. Genes Develop. 2002;16:3100–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Werner T, Motyka V, Strnad M, Schmulling T. Regulation of plant growth by cytokinin. Proc Natl Acad Sci U S A. 2001;98:10487–92.View ArticlePubMed CentralPubMedGoogle Scholar
- Jones B, Gunneras SA, Petersson SV, Tarkowski P, Graham N, May S, et al. Cytokinin regulation of auxin synthesis in Arabidopsis involves a homeostatic feedback loop regulated via auxin and cytokinin signal transduction. Plant Cell. 2010;22:2956–69.View ArticlePubMed CentralPubMedGoogle Scholar
- Hirsch S, Oldroyd GED. GRAS-domain transcription factors that regulate plant development. Plant Signal Behav. 2009;4:698–700.View ArticlePubMed CentralPubMedGoogle Scholar
- Yamashino T, Kitayama M, Mizuno T. Transcription of ST2A encoding a sulfotransferase family protein that is involved in jasmonic acid metabolism is controlled according to the circadian clock- and PIF4/PIF5-mediated external coincidence mechanism in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2013;77:2454–60.View ArticlePubMedGoogle Scholar
- Doğramacı M, Foley ME, Chao WS, Christoffers MJ, Anderson JV. Induction of endodormancy in crown buds of leafy spurge (Euphorbia esula L.) implicates a role for ethylene and cross-talk between photoperiod and temperature. Plant Mol Biol. 2013;81:577–93.View ArticlePubMedGoogle Scholar
- Baur JR, Bovey RW, Veech JA. Growth responses in sorghum and wheat induced by glyphosate. Weed Sci. 1977;25:238–40.Google Scholar
- Roider CA, Griffin JL, Harrison SA, Jones CA. Wheat response to simulated glyphosate drift. Weed Technol. 2007;21:1010–5.View ArticleGoogle Scholar
- Kutman BY, Kutman UB, Cakmak I. Foliar nickel application alleviates detrimental effects of glyphosate drift on yield and seed quality of wheat. J Agr Food Chem. 2013;61:8364–72.View ArticleGoogle Scholar
- Wang CY. Effect of glyphosate on aromatic amino acid metabolism in purple nutsedge (Cyperus rotundus). Weed Technol. 2001;15:628–35.View ArticleGoogle Scholar
- Shaner DL, Nadler-Hassar T, Henry WB, Koger CH. A rapid in vivo shikimate accumulation assay with excised leaf discs. Weed Sci. 2005;53:769–74.View ArticleGoogle Scholar
- Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL, et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri. Proc Natl Acad Sci U S A. 2010;107:1029–34.View ArticlePubMed CentralPubMedGoogle Scholar
- Nol N, Tsikou D, Eid M, Livieratos IC, Giannopolitis CN. Shikimate leaf disc assay for early detection of glyphosate resistance in Conyza canadensis and relative transcript levels of EPSPS and ABC transporter genes. Weed Res. 2012;52:233–41.View ArticleGoogle Scholar
- Henry WB, Shaner DL, West MS. Shikimate accumulation in sunflower, wheat, and proso millet after glyphosate application. Weed Sci. 2007;55:1–5.View ArticleGoogle Scholar
- Anderson JV, Gesch RW, Jia Y, Chao WS, Horvath DP. Seasonal shifts in dormancy status, carbohydrate metabolism, and related gene expression in crown buds of leafy spurge. Plant Cell Environ. 2005;28:1567–78.View ArticleGoogle Scholar
- Sandberg CL, Meggit WF, Penner D. Absorption, translocation and metabolism of 14C labeled glyphosate in several weed species. Weed Res. 1980;20:195–200.View ArticleGoogle Scholar
- Hoffmann L, Besseau S, Geoffroy P, Ritzenthaler C, Meyer D, Lapierre C, et al. Silencing of hydroxycinnamoyl-coenzyme a shikimate/quinate hydroxycinnamoyltransferase affects phenylpropanoid biosynthesis. Plant Cell. 2004;16:1446–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Li X, Bonawitz ND, Weng JK, Chapple C. The growth reduction associated with repressed lignin biosynthesis in Arabidopsis thaliana is independent of flavonoids. Plant Cell. 2010;22:1620–32.View ArticlePubMed CentralPubMedGoogle Scholar
- Henry WB. Koger CH. Accumulation of shikimate in corn and soybean exposed to various rates of glyphosate. Crop Manage: Shaner DL; 2005.Google Scholar
- Mueller TC, Tellis A, Beeler JE, Sharma SD, Singh M. Shikimate accumulation in nine weedy species following glyphosate application. Weed Res. 2008;48:455–60.View ArticleGoogle Scholar
- Duke SO. Glyphosate degradation in glyphosate-resistant and -susceptible crops and weeds. J Agr Food Chem. 2010;59:5835–41.View ArticleGoogle Scholar
- Reddy KN, Rimando AM, Duke SO. Aminomethylphosphonic acid, a metabolite of glyphosate, causes injury in glyphosate-treated, glyphosate-resistant soybean. J Agr Food Chem. 2004;52:5139–43.View ArticleGoogle Scholar
- Chao WS, Serpe MD, Anderson JA, Gesch RW, Horvath DP. Sugars, hormones, and environment affects the dormancy status in underground adventitious buds of leafy spurge (Euphorbia esula). Weed Sci. 2006;54:59–68.View ArticleGoogle Scholar
- Stasiak MA, Hofstra G, Fletcher RA. Physiological changes induced in birch seedlings by sublethal applications of glyphosate. Can J for Res. 1992;22:812–7.View ArticleGoogle Scholar
- Geiger DR, Bestman HD. Self-limitation of herbicide mobility by phytotoxic action. Weed Sci. 1990;38:324–9.Google Scholar
- Mason MG, Ross JJ, Babst BA, Wienclaw BN, Beveridge CA. Sugar demand, not auxin, is the initial regulator of apical dominance. Proc Natl Acad Sci U S A. 2014;111:6092–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Horvath DP. The role of specific plant organs and polar auxin transport in correlative inhibition of leafy spurge (Euphorbia esula L.) root buds. Can J Bot. 1998;76:1227–31.Google Scholar
- Stirnberg P, Furner IJ, Ottoline Leyser HM. MAX2 participates in an SCF complex which acts locally at the node to suppress shoot branching. Plant J. 2007;50:80–94.View ArticlePubMedGoogle Scholar
- Gomez-Roldan V, Fermas S, Brewer PB, Puech-Pages V, Dun EA, Pillot JP, et al. Strigolactone inhibition of shoot branching. Nature. 2008;455:189–94.View ArticlePubMedGoogle Scholar
- Umehara M, Hanada A, Yoshida S, Akiyama K, Arite T, Takeda-Kamiya N, et al. Inhibition of shoot branching by new terpenoid plant hormones. Nature. 2008;455:195–200.View ArticlePubMedGoogle Scholar
- Hayward A, Stirnberg P, Beveridge C, Leyser O. Interactions between auxin and strigolactone in shoot branching control. Plant Physiol. 2009;151:400–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Hamiaux C, Drummond RSM, Janssen BJ, Ledger SE, Cooney JM, Newcomb RD, et al. DAD2 is an alpha_beta hydrolase likely to be involved in the perception of the plant branching hormone, strigolactone. Curr Biol. 2012;22:2032–6.View ArticlePubMedGoogle Scholar
- Peer WA. From perception to attenuation: auxin signalling and responses. Curr Opin Plant Biol. 2013;16:561–8.View ArticlePubMedGoogle Scholar
- Blilou I, Xu J, Wildwater M, Willemsen V, Paponov I, Friml J, et al. The PIN auxin efflux facilitator network controls growth and patterning in Arabidopsis roots. Nature. 2005;433:39–44.View ArticlePubMedGoogle Scholar
- Domagalska MA, Leyser O. Signal integration in the control of shoot branching. Nat Rev Mol Cell Biol. 2011;12:211–21.View ArticlePubMedGoogle Scholar
- Bennett T, Sieberer T, Willett B, Booker J, Luschnig C, Leyser O. The Arabidopsis MAX pathway controls shoot branching by regulating auxin transport. Curr Biol. 2006;16:553–63.View ArticlePubMedGoogle Scholar
- Prusinkiewicz P, Crawford S, Smith RS, Ljung K, Bennett T, Ongaro V, et al. Control of bud activation by an auxin transport switch. Proc Natl Acad Sci U S A. 2009;106:17431–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Cummins I, Wortley DJ, Sabbadin F, He Z, Coxon CR, Straker HE, et al. Key role for a glutathione transferase in multiple-herbicide resistance in grass weeds. Proc Natl Acad Sci U S A. 2013;110:5812–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Cummins I, Moss S, Cole DJ, Edwards R. Glutathione trasferases in herbicide-resistant and herbicide-susceptible black-grass (Alopercurus myosuroides). Pestic Sci. 1997;51:244–50.View ArticleGoogle Scholar
- Conte SS, Lloyd AM. Exploring multiple drug and herbicide resistance in plant - spotlight on transporter proteins. Plant Sci. 2011;180:196–203.View ArticlePubMedGoogle Scholar
- Werck-Reichhart D, Hehn A, Didierjean L. Cytochromes P450 for engineering herbicide tolerance. Trends Plant Sci. 2000;5:116–23.View ArticlePubMedGoogle Scholar
- Brazier M, Cole DJ, Edwards R. O-Glucosyltransferase activities toward phenolic natural products and xenobiotics in wheat and herbicide-resistant and herbicide-susceptible black-grass (Alopecurus myosuroides). Phytochemistry. 2002;59:149–56.View ArticlePubMedGoogle Scholar
- Yanhui C, Xiaoyuan Y, Kun H, Meihua L, Jigang L, Zhaofeng G, et al. The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2007;60:107–24.View ArticleGoogle Scholar
- Busi R, Vila-Aiub MM, Powles SB. Genetic control of cytochrome P450 metabolism-based herbicide resistance mechanism in Lolium rigidum. Heredity. 2011;106:817–24.View ArticlePubMed CentralPubMedGoogle Scholar
- Campbell WF, Evans JO, Reed SC. Effect of glyphosate on chloroplast ultrastructure of quackgrass mesophyll cells. Weed Sci. 1976;24:22–5.Google Scholar
- Mollenhauer C, Smart CC, Amrhein N. Glyphosate toxicity in the shoot apical region of the tomato plant I. Plastid swelling is the initial ultrastructural feature following in vivo inhibition of 5-enolpyruvylshikimic acid 3-phosphate synthase. Pestic Biochem and Phys. 1987;29:55–65.View ArticleGoogle Scholar
- Neve EPA, Ingelman-Sundberg M. Cytochrome P450 proteins: retention and distribution from the endoplasmic reticulum. Curr Opin Drug Discov Devel. 2010;13:78–85.PubMedGoogle Scholar
- Sponsel VM, Hedden P. Gibberellin biosynthesis and inactivation. In: Davies PJ, editor. Plant Hormones: Biosynthesis, Signal Transduction, Action! The Netherlands: Kluwer Academic Publishers; 2004. p. 63–94.Google Scholar
- Ge X, d'Avignon DA, Ackerman JJH, Douglas Sammons R. Rapid vacuolar sequestration: The horseweed glyphosate resistance mechanism. Pest Manag Sci. 2010;66:345–8.PubMed CentralPubMedGoogle Scholar
- Anderson JV, Davis DG. Abiotic stress alters transcript profiles and activity of glutathione S-transferase, glutathione peroxidase, and glutathione reductase in Euphorbia esula. Physiol Plantarum. 2004;120:421–33.View ArticleGoogle Scholar
- Foley ME, Anderson JV, Horvath DP. The effects of temperature, photoperiod, and vernalization on regrowth and flowering competence in Euphorbia esula (Euphorbiaceae) crown buds. Botany. 2009;87:986–92.View ArticleGoogle Scholar
- Cromartie TH, Polge ND. An improved assay for shikimic acid and its use as a monitor for the activity of sufosate. Proc Weed Sci Soc Am. 2000;40:291.Google Scholar
- Abrams SR, Nelson K, Ambrose SJ. Deuterated abscisic acid analogs for mass spectrometry and metabolism studies. J Labelled Compd Radiopharm. 2003;46:273–83.View ArticleGoogle Scholar
- Zaharia LI, Galka MM, Ambrose SJ, Abrams SR. Preparation of deuterated abscisic acid metabolites for use in mass spectrometry and feeding studies. J Labelled Compd Radiopharm. 2005;48:435–45.View ArticleGoogle Scholar
- Chiwocha SDS, Abrams SR, Ambrose SJ, Cutler AJ, Loewen M, Ross ARS, et al. A method for profiling classes of plant hormones and their metabolites using liquid chromatography-electrospray ionization tandem mass spectrometry: an analysis of hormone regulation of thermodormancy of lettuce (Lactuca sativa L.) seeds. Plant J. 2003;35:405–17.View ArticlePubMedGoogle Scholar
- Chiwocha SDS, Cutler AJ, Abrams SR, Ambrose SJ, Yang J, Ross ARS, et al. The etr1-2 mutation in Arabidopsis thaliana affects the abscisic acid, auxin, cytokinin and gibberellin metabolic pathways during maintenance of seed dormancy, moist-chilling and germination. Plant J. 2005;42:35–48.View ArticlePubMedGoogle Scholar
- Chang S, Puryear J, Cairney J. A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993;11:113–6.View ArticleGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nature Prot. 2013;8:1494–512.View ArticleGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. MC Bioinformatics. 2011;12:323.View ArticleGoogle Scholar
- Leng N, Dawson JA, Thompson JA, Ruotti V, Rissman AI, Smits BMG, et al. EBSeq: An empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics. 2013;29:1035–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Anderson JV, Horvath DP, Chao WS, Foley ME, Hernandez AG, Thimmapuram J, et al. Characterization of an EST database for the perennial weed leafy spurge: An important resource for weed biology research. Weed Sci. 2007;55:193–203.View ArticleGoogle Scholar
- Chao WS, Doğramacı M, Foley ME, Horvath DP, Anderson JV. Selection and validation of endogenous reference genes for qRT-PCR analysis in leafy spurge (Euphorbia esula). PLoS One. 2012;7, e42839.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.