Introduction
Halauxifen-methyl (HM) is categorized in the 6-aryl-picolinic acid subclass of synthetic auxin herbicides (Group 4) and was commercialized in 2015 (Epp et al. Reference Epp, Alexander, Balko, Buysse, Brewster, Bryan, Daeuble, Fields, Gast, Green, Irvine, Lo, Lowe, Renga and Richburg2016; Schmitzer et al. Reference Schmitzer, Balko, Daeuble, Epp, Satchivi, Siddall, Weimer and Yerkes2015). Like other synthetic auxins, HM mimics indole-3-acetic acid, which regulates numerous aspects of plant growth and development (Demeulenaere and Beeckman Reference Demeulenaere and Beeckman2014; McSteen Reference McSteen2010). Halauxifen-methyl provides selective postemergence control of dicot weeds in cereal crops (Dzikowski et al. Reference Dzikowski, Becker, Larelle, Kamerichs and Gast2016; Epp et al. Reference Epp, Alexander, Balko, Buysse, Brewster, Bryan, Daeuble, Fields, Gast, Green, Irvine, Lo, Lowe, Renga and Richburg2016). While HM was initially developed as a tank-mix partner for postemergence weed control in cereal crops at relatively low application rates (5 to 7.5 g ha−1), it can also be used as a burndown treatment before soybean [Glycine max (L.) Merr.] planting at even lower rates (1 to 2 g ha−1) due to its short soil half-life (10 to 25 d) and minimal soil residual activity (Epp et al. Reference Epp, Alexander, Balko, Buysse, Brewster, Bryan, Daeuble, Fields, Gast, Green, Irvine, Lo, Lowe, Renga and Richburg2016, Reference Epp, Schmitzer and Crouse2018).
Allohexaploid bread wheat (Triticum aestivum L.) (2n = 6x = 42; AABBDD) is tolerant to HM via rapid detoxification of the biologically active form, halauxifen acid (HA) (Dzikowski et al. Reference Dzikowski, Becker, Larelle, Kamerichs and Gast2016) following phase I and II metabolic reactions. Phase I reactions typically involve oxidation, reduction, or hydrolysis of phytotoxic parent molecules; most reactions are catalyzed by cytochrome P450 enzymes (CYPs) or bioactivation of proherbicide esters by carboxyl-esterases (CXEs) (Gershater and Edwards Reference Gershater and Edwards2007; Kreuz et al. Reference Kreuz, Tommasini and Martinoia1996). Phase II reactions involve conjugation of phytotoxic parent molecules or phase I metabolites with endogenous compounds, such as sugars or reduced glutathione (Kreuz et al. Reference Kreuz, Tommasini and Martinoia1996; Riechers et al. Reference Riechers, Kreuz and Zhang2010). Following de-esterification of HM to HA, HA is O-demethylated and subsequently conjugated with glucose (Figure 1). Although HM metabolism in wheat has been characterized (Dzikowski et al. Reference Dzikowski, Becker, Larelle, Kamerichs and Gast2016), the genes encoding HA-detoxifying enzymes have not been identified. Identification of such genes would be valuable for plant transformation to endow HM tolerance in HM-sensitive plants, such as soybean (Zobiole and Kalsing Reference Zobiole and Kalsing2017), cotton (Gossypium hirsutum L.), and other dicots. Additionally, these enzymes could potentially detoxify multiple herbicides, making them valuable in breeding programs aimed at improving herbicide tolerance in dicot crops (Dimaano and Iwakami Reference Dimaano and Iwakami2021; Gharabli et al. Reference Gharabli, Della Gala and Welner2023).
Various genetic resources derived from the allohexaploid wheat cultivar ‘Chinese Spring’ (CS) are available for evaluating phenotypic responses via alteration of the base chromosome number or substituting chromosomes, including alien substitution and nullisomic-tetrasomic (NT) lines (Jiang et al. Reference Jiang, Friebe and Gill1994; Law et al. Reference Law, Snape and Worland1987). Alien substitution lines have one pair of endogenous homoeologous chromosomes replaced with the homoeologous chromosomes from another species, denoted as the “alien” genome (Jiang et al. Reference Jiang, Friebe and Gill1994). For NT lines, one pair of homoeologues are removed (nullisomic) and additional copies of another pair of homoeologues are added (tetrasomic) (Law et al. Reference Law, Snape and Worland1987; Sears Reference Sears1954, Reference Sears1966). These lines were utilized in previous HM phenotyping experiments with the goals of identifying lines with increased HM sensitivity and providing evidence concerning which chromosomes harbor genes necessary to confer HM tolerance in wheat (Obenland and Riechers Reference Obenland and Riechers2020).
Increased sensitivity to HM was clearly evident in wheat lines that lack chromosome 5A, relative to CS (Obenland and Riechers Reference Obenland and Riechers2020), but HM tolerance in wheat lines possessing 5A was similar to that of CS. In our prior research, we hypothesized that increased HM sensitivity in these wheat lines is due to reduced HA detoxification rates, which likely resulted from losing genes on chromosome 5A encoding HA-detoxifying enzymes (Obenland and Riechers Reference Obenland and Riechers2020). To directly address this hypothesis, the objective of this study was to measure and compare HA levels in CS and wheat lines utilized in our previous research via a combination of excised-leaf assays and liquid chromatography–mass spectrometry (LC-MS) analysis.
Materials and Methods
Chemicals, Plant Materials, and Growth Conditions
Analytical-grade standards of HM and HA (95.9% and 99.5% pure, respectively) were provided by Corteva Agriscience (Indianapolis, IN, USA). All other analytical-grade reagents were purchased through Fisher Scientific (Thermo-Fisher, Hanover Park, IL, USA) or Sigma Chemical (Millipore Sigma, St Louis, MO, USA). Seeds for CS, alien substitution lines (5A, 5B, and 5D), and NT lines were acquired from the Kansas State Wheat Genetics Resource Center; seed of Sears’ goatgrass (Aegilops searsii M. Feldman & M. Kislev) (AS; PI 599163), a diploid wheat relative (2n = 2x = 14; S genome), was obtained from the National Small Grains Collection (Aberdeen, ID, USA) of the U.S. Department of Agriculture–Agricultural Research Service. To promote uniform germination, seeds were subjected to a cold treatment by placing them on water-soaked filter paper in Petri dishes stored at 5 C for 3 d. Seedlings were cultivated in a growth chamber (Conviron Gen1000, Controlled Environments, Winnipeg, Canada) until seedlings produced two to three leaves (Zadoks stage 12-13) under conditions of 28/22 C day/night and a 16/8-h photoperiod. The LED lights in the growth chamber provided 550 μmol m−2 s−1 photon flux at the plant canopy level.
Excised-Leaf Assay, Tissue Harvest and Extraction, and LC-MS Analysis
The excised-leaf assay and metabolite extraction protocols utilized in our experiments are based on a protocol modified from previous studies (Concepcion et al. Reference Concepcion, Kaundun, Morris, Hutchings, Strom, Lygin and Riechers2021; Lygin et al. Reference Lygin, Kaundun, Morris, Mcindoe, Hamilton and Riechers2018; Ma et al. Reference Ma, Skelton and Riechers2015). For each experiment, the two oldest wheat leaves were cut at the collar, the cut surface of the leaf was washed with deionized water, transferred to 15-ml plastic tubes containing 400 µl of 0.1 M Tris-Cl (pH 6.0), and preincubated in the growth chamber for 1 h. After preincubation, leaves designated for HM treatment were transferred to tubes containing 400 µl of 300 μM HM (in 0.1 M Tris-Cl, pH 6.0) and untreated leaves were transferred to tubes containing only 400 µl of 0.1 M Tris-Cl (pH 6.0). Approximately 1.0 cm of tissue above the cut surface of the excised leaf was exposed to the uptake solution in these assays. At 2 h after treatment (HAT), leaves designated for the 2-HAT time point were rinsed with deionized water, and fresh weights were recorded (approx. 200 mg of fresh weight per sample). Samples designated for later time points (8, 12, and 24 HAT) were rinsed with deionized water and transferred to tubes containing 400 µl of quarter-strength Murashige and Skoog solution (MP Biomedicals, Solon, OH, USA), which provides mineral salts, vitamins, and organic compounds to support plant growth (Murashige and Skoog Reference Murashige and Skoog1962). Following each harvest time point, leaves were frozen in liquid nitrogen, and samples were stored at −80 C until metabolite extraction. Based on preliminary experiments quantifying HM remaining in the uptake solution at 2 HAT via high-performance liquid chromatography with photodiode array detection (HPLC-PDA) using methods (Supplementary Material 1) modified from a previous study (Concepcion et al. Reference Concepcion, Kaundun, Morris, Hutchings, Strom, Lygin and Riechers2021), at least 19% of HM supplied was absorbed by excised wheat leaves (Supplementary Material 2). However, AS absorbed a greater amount of HM than CS, which may be related to anatomical differences between these species. Substitution line 5D and the nullisomic 5D-tetrasomic 5A (N5D-T5A) line also absorbed more HM than CS (Supplementary Material 2), which positively correlates with a relatively higher degree of HM tolerance (Obenland and Riechers Reference Obenland and Riechers2020) and a putative faster rate of HM metabolism.
Frozen tissue was ground in liquid nitrogen using a mortar and pestle, and compounds were extracted in 1 ml of 90% (v/v) methanol. The supernatant was removed after centrifugation at 12,000 × g, and a second extraction was performed with the remaining plant material by adding 1 ml of 90% (v/v) methanol. After another centrifugation at 12,000 × g, the pellet was discarded, and the second supernatant was combined with the first supernatant, resulting in a final volume of 2 ml. All samples were spiked with 4-chloro-DL-PHE as an internal standard before extraction. Samples were dried and concentrated under nitrogen gas and reconstituted with 250 μl of water:acetonitrile (1:1, v/v) containing 0.1% formic acid. Quality-control (QC) samples were then prepared by combining aliquots of each sample for injection throughout each experimental run (Dunn et al. Reference Dunn, Broadhurst, Begley, Zelena, Francis-Mcintyre, Anderson, Brown, Knowles, Halsall, Haselden, Nicholls, Wilson, Kell and Goodacre2011). An HA standard was also prepared. All samples were stored at −80 C until further analysis.
Samples were submitted to the Roy J. Carver Metabolomics Facility, University of Illinois at Urbana-Champaign, for LC-MS analysis using previously described methods (Concepcion et al. Reference Concepcion, Kaundun, Morris, Hutchings, Strom, Lygin and Riechers2021; Elolimy et al. Reference Elolimy, Alharthi, Zeineldin, Parys, Helmbrecht and Loor2019). Samples were analyzed with a Q-Exactive MS system (Thermo, Bremen, Germany), and separation was conducted with a Dionex Ultimate 3000 series HPLC equipped with a Phenomenex C18 column (4.6 by100 mm, 2.6-μm particle size). Mobile phases consisted of A (H2O with 0.1% formic acid [v/v]) and B (acetonitrile with 0.1% formic acid [v/v]). The flow rate was 0.25 ml min−1 with a linear gradient starting at 100% A for 3 min. The gradient then transitioned to 100% B (20 to 30 min) and returned to 100% A (31 to 36 min). Twenty microliters of each sample was injected, and the autosampler temperature was 15 C. Mass spectra were then acquired under both positive (sheath gas flow rate: 45; auxiliary gas flow rate: 11; sweep gas flow rate: 2; spray voltage: 3.5 kV; capillary temp: 250 C; auxiliary gas heater temp: 415 C) and negative electrospray ionization (sheath gas flow rate: 45; auxiliary gas flow rate: 11; sweep gas flow rate: 2; spray voltage: −2.5 kV; capillary temp: 250 C; auxiliary gas heater temp: 415 C). The full-scan mass-spectrum resolution was set to 70,000 with a scan range of m/z 67 to 1,000, and the automatic gain control target was 1 × 106 with a maximum injection time of 200 ms. Chromatographic analysis was conducted in a randomized sequence, including QC samples injected initially to equilibrate the analytical platform and after every 10 test samples to evaluate the stability of the experimental procedure (Dunn et al. Reference Dunn, Broadhurst, Begley, Zelena, Francis-Mcintyre, Anderson, Brown, Knowles, Halsall, Haselden, Nicholls, Wilson, Kell and Goodacre2011; Godzien et al. Reference Godzien, Alonso-Herranz, Barbas and Armitage2015; Sangster et al. Reference Sangster, Major, Plumb, Wilson and Wilson2006; Wehrens et al. Reference Wehrens, Hageman, van Eeuwijk, Kooke, Flood, Wijnker, Keurentjes, Lommen, van Eekelen, Hall, Mumm and de Vos2016). Raw data files obtained in full-MS mode (samples, procedural blank, and QC) and data obtained in full-MS followed by data-dependent MS2 were processed using MS-DIAL v. 4.9221218 following previously described parameters (Concepcion et al. Reference Concepcion, Kaundun, Morris, Hutchings, Strom, Lygin and Riechers2021; Tsugawa et al. Reference Tsugawa, Cajka, Kind, Ma, Higgins, Ikeda, Kanazawa, Vandergheynst, Fiehn and Arita2015). The peak area of each metabolite feature was normalized based on QC samples and peak area of the internal standard (Supplementary Material 2). Normalized peak areas for HA were then exported into a .csv file for further statistical analysis.
Experimental Setup and Statistical Analysis
The first experiment included CS, AS, and three alien substitution lines derived from CS wheat but with endogenous chromosomes 5A, 5B, or 5D replaced with the corresponding homoeologous chromosome (5S) from AS (Obenland and Riechers Reference Obenland and Riechers2020). Based on the results of the first experiment, NT lines with different numbers of chromosomes 5A and 5D (i.e., 5A or 5D were either omitted [nullisomic] or doubled [tetrasomic]) were used for the second experiment, which included CS, nullisomic 5A-tetrasomic 5D (N5A-T5D), and N5D-T5A. Both experiments were designed to meet or exceed minimum requirements for replicating and reporting metabolite profiling results (Sumner et al. Reference Sumner, Amberg, Barrett, Beale, Beger, Daykin, Fan, Fiehn, Goodacre, Griffin, Hankemeier, Hardy, Harnly, Higashi and Kopka2007). The first experiment was independently conducted twice using either three or four biological replicates per treatment per time point (seven total biological replicates). The second experiment was conducted once and utilized five biological replicates per treatment per time point. Chinese Spring served as a positive control in both experiments, as it is HM tolerant, while AS served as a negative control due to its sensitivity to HM (Obenland and Riechers Reference Obenland and Riechers2020). Comparisons of HA normalized peak areas among populations and time points were performed with the lme4 package (Bates et al. Reference Bates, Mächler, Bolker and Walker2015) in R (v. 4.2.0) using RStudio (v. 2023.03.0). Both experiments used a mixed-effects model in which population, time point, and their interactions were treated as fixed effects and replicates were random effects. Homogeneity of variance was confirmed with the Levene’s test. Data were subjected to ANOVA, and treatment means were separated with Fisher’s LSD (α = 0.05).
Results and Discussion
LC-MS Analysis of HA Levels in CS, AS, Alien Substitution, and NT Lines
Peak abundances of HA in the first experiment with CS and AS exhibited anticipated results as positive and negative controls, respectively; the amounts of HA in CS remained relatively low throughout the time course, while HA peak abundances in AS increased during the time course and were greater than HA levels in CS and all substitution lines at 8, 12, and 24 HAT (Figure 2). At these same time points, HA levels in AS were 8.5-, 7.4-, and 5.6-fold higher, respectively, than HA levels in CS. Among substitution lines, only 5A accumulated HA relative to CS at 8, 12, and 24 HAT (Figure 2), with HA peak abundances 3.8-, 3.2- and 2.4-fold higher than HA levels in CS, respectively. In contrast, the level of HA in 5D was not different from HA levels in CS at any time point.
By comparison, HA levels in the 5B line (Figure 2) varied throughout the time course relative to CS, 5A, and 5D, which were more consistent between 8 and 24 HAT. The peak abundance of HA in 5B was not different from HA levels in CS or 5D at 8, 12, and 24 HAT, but HA abundance in 5B increased to levels also not different from HA levels in 5A at 12 and 24 HAT (Figure 2). This finding indicates the rate of HA metabolism in 5B at 12 and 24 HAT is slower than that in CS and 5D compared with its rate of metabolism at 8 HAT or, alternatively, that rates of HA formation from HM are higher in 5B at 12 and 24 HAT compared with 8 HAT. These metabolism results are in accord with previous phenotypic experiments in the greenhouse, where the 5B line displayed intermediate HM sensitivity relative to the highly sensitive 5A line and relatively tolerant 5D line (Obenland and Riechers Reference Obenland and Riechers2020). The variation in HA peak abundance in the 5B line during the time course may indicate competing CXE versus CYP enzyme activities (further described later), but additional research is needed to investigate this hypothesis. However, because HA peak abundances in lines 5A and 5D consistently differed at 8, 12, and 24 HAT, we focused on NT lines with varying numbers of 5A and 5D chromosomes (e.g., N5A-T5D and N5D-T5A) for subsequent LC-MS analysis.
In the second experiment, HA peak abundance in CS remained relatively low and did not significantly change throughout the time course (Figure 3), which was also noted in the previous experiment (Figure 2). Only line N5A-T5D displayed higher amounts of HA relative to CS at 12 and 24 HAT (approximately 3-fold; Figure 3). HA peak abundances in N5D-T5A did not differ from HA levels in CS throughout the time course (Figure 3).
These results indicate that species and lines lacking chromosome 5A (AS, N5A-T5D, and the 5A substitution line) have a relatively slow rate of HA metabolism. By contrast, lines possessing chromosome 5A (CS, N5D-T5A, and the 5B and 5D substitution lines) maintain HA metabolism rates similar to rates in CS. These findings support the hypothesis that genes encoding HA-detoxifying enzymes are located on wheat chromosome 5A, and as a result, lines lacking 5A consistently accumulate greater amounts of phytotoxic HA relative to lines possessing 5A through 24 HAT. Additionally, these results corroborate findings from previous greenhouse experiments in which 5A, N5A-T5D, and AS were relatively sensitive to HM, while CS, N5D-T5A, 5B, and 5D maintained tolerance to HM (Obenland and Riechers Reference Obenland and Riechers2020).
As HM is absorbed by the cuticle and underlying leaf tissue, it is de-esterified to HA (Figure 1) at an unknown rate. HA is a phytotoxic, transient, phase I compound that is subsequently O-demethylated, then conjugated with glucose (Dzikowski et al. Reference Dzikowski, Becker, Larelle, Kamerichs and Gast2016; Figure 1). O-demethylation and glucose conjugation reactions are likely catalyzed by CYPs and UDP-dependent glucosyltransferases (UGTs), respectively, which commonly mediate synthetic auxin herbicide–detoxification reactions in tolerant grasses (Frear Reference Frear1995; Mithila et al. Reference Mithila, Hall, Johnson, Kelley and Riechers2011; Sterling and Hall Reference Sterling and Hall1997; Zhang and Yang Reference Zhang and Yang2021). The presence of HA in all species and lines investigated (Figures 2 and 3) implies the activity of CXEs that convert HM to HA (Gershater and Edwards Reference Gershater and Edwards2007). However, the activity of CXEs, CYPs, and UGTs and resulting metabolic flux could vary, which could shed more light on our current results. Given differential absorption of HM among some species and lines (Supplementary Material 2), future investigation of CXEs is warranted. For example, the abundance of HA in the 5B or 5D substitution lines during the time course (Figure 2) could be relatively low due to loss of an active CXE that converts HM to HA. Conducting a follow-up metabolism study and supplying HA instead of HM in our assay may help elucidate the possible role of differential CXE activity among wheat lines in future research.
Additional biokinetic analyses should also measure abundances of the O-demethylated metabolite and glucose conjugate of HA (Figure 1) to comprehensively understand altered HA detoxification rates or patterns in these substitution and NT lines. Such investigations would require analytical metabolite standards or radiolabeled HM, which were not available for the current study. Because these initial metabolites could be present but occur in low abundance—particularly the O-demethylated HA metabolite if it is rapidly turned over to phase II conjugates—a whole-plant assay may be more informative, especially if quantifying HM absorption, HA translocation, or measuring HA metabolites beyond 24 h is desired.
The discovery of HM sensitivity in AS was fortuitous for our prior research (Obenland and Riechers Reference Obenland and Riechers2020). Detection of a relatively low HA metabolism rate (Figure 2) is notable because grasses are generally tolerant to synthetic auxins, with the exception of quinoline carboxylic acids (Grossmann Reference Grossmann2010; Sterling and Hall Reference Sterling and Hall1997) and a different 6-aryl-picolinic acid herbicide, florpyrauxifen-benzyl (Epp et al. Reference Epp, Alexander, Balko, Buysse, Brewster, Bryan, Daeuble, Fields, Gast, Green, Irvine, Lo, Lowe, Renga and Richburg2016; Lim et al. Reference Lim, Kim, Noh, Lim, Yook, Kim, Yi and Kim2021; Miller et al. Reference Miller, Norsworthy and Scott2018; Takano et al. Reference Takano, Greenwalt, Ouse, Zielinski and Schmitzer2023). Interestingly, the wheat safener cloquintocet-mexyl (3.75 g ai ha−1) reduced injury in AS even at an elevated rate of HM (60 g ae ha−1 plus adjuvants; Obenland and Riechers Reference Obenland and Riechers2020), which indicates enhancement of herbicide tolerance (presumably via induction of metabolic enzymes) by safeners is not limited to cultivated cereal crops (Riechers et al. Reference Riechers, Kreuz and Zhang2010).
Our results inform future experiments, including RNA-seq or genome-wide association studies, aimed at identifying potential candidate genes on chromosome 5A. However, these experiments often yield numerous (>50) potential candidate genes (Baek et al. Reference Baek, Goodrich, Brown, James, Moose, Lambert and Riechers2019). This is especially true for wheat, which has a relatively large (approximately 16 GB) allohexaploid genome (IWGSC Reference Epp, Alexander, Balko, Buysse, Brewster, Bryan, Daeuble, Fields, Gast, Green, Irvine, Lo, Lowe, Renga and Richburg2018), meaning a significant gene might be identified along with its homoeologues. This concern is exacerbated since CYPs and UGTs are encoded by large multigene families per diploid genome (Casey and Dolan Reference Casey and Dolan2023; Gharabli et al. Reference Gharabli, Della Gala and Welner2023). Previous studies utilizing various wheat species, alien substitution lines, or NT lines have identified differential homoeologous chromosome contributions to herbicide tolerance, metabolism, and gene expression patterns (Riechers et al. Reference Riechers, Yang, Irzyk, Jones and Fuerst1996, Reference Riechers, Irzyk, Jones and Fuerst1997; Xu et al. Reference Xu, Lagudah, Moose and Riechers2002; Yu et al. Reference Yu, Guo, Cui, Chen and Li2023).
Potential candidate genes encoding HA-detoxifying enzymes include CYPs and UGTs (Frear Reference Frear1995; Mithila et al. Reference Mithila, Hall, Johnson, Kelley and Riechers2011; Riechers et al. Reference Riechers, Kreuz and Zhang2010; Zhang and Yang Reference Zhang and Yang2021). CYPs have been implicated in detoxification of numerous herbicides in the Poaceae, which has the highest number of herbicide-detoxifying plant CYPs (Brazier-Hicks et al. Reference Brazier-Hicks, Franco-Ortega, Watson, Rougemont, Cohn, Dale, Hawkes, Goldberg-Cavalleri, Onkokesung and Edwards2022; Casey and Dolan Reference Casey and Dolan2023; Dimaano and Iwakami Reference Dimaano and Iwakami2021; Han et al. Reference Han, Yu, Beffa, González, Maiwald, Wang and Powles2020; Pan et al. Reference Pan, Guo, Wang, Shi, Yang, Zhou, Yu and Bai2022; Yu et al. Reference Yu, Guo, Cui, Chen and Li2023; Zheng et al. Reference Zheng, Yu, Sun, Zhang, Zhang, Tang, Lin and Shen2022). Members of the plant CYP81A subfamily are the most common, but CYP72A and CYP71C members have also been identified (Supplementary Material 1). Compared with CYPs, genes encoding herbicide-detoxifying UGTs in plants are generally not as well characterized. UGT activities in wheat have been directly studied with natural product and xenobiotic substrates (Brazier et al. Reference Brazier, Cole and Edwards2002), and the activities of herbicide-detoxifying UGTs have been inferred from metabolism studies in wheat in which glycosylated metabolites of isoproturon (Lu et al. Reference Lu, Zhang and Yang2015) and florasulam (DeBoer et al. Reference DeBoer, Thornburgh and Ehr2006) were detected.
Regulatory elements and factors, such as promoters, untranslated regions (UTRs), and transcription factors (TFs), of candidate genes may also contribute to the prominent role of genes on wheat chromosome 5A in controlling HA detoxification compared with its homoeologous chromosomes 5B and 5D. The role of TFs in herbicide detoxification is often inferred, as RNA-seq has identified TFs regulated by herbicides and safeners (Baek et al. Reference Baek, Goodrich, Brown, James, Moose, Lambert and Riechers2019; Brazier-Hicks et al. Reference Brazier-Hicks, Howell, Cohn, Hawkes, Hall, McIndoe and Edwards2020; Chen et al. Reference Chen, Su, Yang, Peng, Wu, Zhao, Lin, Bai and Wang2023; Duhoux et al. Reference Duhoux, Carrère, Gouzy, Bonin and Délye2015; Gaines et al. Reference Gaines, Lorentz, Figge, Herrmann, Maiwald, Ott, Han, Busi, Yu, Powles and Beffa2014; Zhang et al. Reference Zhang, Gao, Fang, Wang, Chen, Li and Dong2022). Homoeologue expression bias is commonly observed in allopolyploid plants (Grover et al. Reference Grover, Gallagher, Szadkowski, Yoo, Flagel and Wendel2012), partly attributed to variations in cis- and trans-regulatory elements between homoeologues (He et al. Reference He, Wang, Rutter, Jordan, Ren, Taagen, DeWitt, Sehgal, Sukumaran, Dreisigacker, Reynolds, Halder, Sehgal, Liu and Chen2022; Kremling et al. Reference Kremling, Chen, Su, Lepak, Romay, Swarts, Lu, Lorant, Bradbury and Buckler2018). We therefore hypothesize potential candidate genes on chromosome 5A and their 5B and 5D homoeologues vary in promoter/UTR sequences and/or TF binding that regulates gene expression, resulting in relatively higher expression and activity of 5A homoeologues. However, gene copy number variation or differences in encoded protein sequences among homoeologous genes should not be precluded without further experimentation.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/wsc.2024.24
Data availability
Data from LC-MS analyses are available on MassIVE (https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp; reference numbers are MSV000092323, MSV000092320, and MSV000092256).
Acknowledgments
We thank Norbert Satchivi and Carla Yerkes (Corteva Agrisciences) for providing technical advice and standards for experimentation, Harold Bockelman of the USDA-ARS and Jon Raupp of the Kansas State Wheat Genetics Resource Center for assistance in obtaining wheat seeds, and Michael La Frano and Alexander Ulanov (University of Illinois) for providing LC-MS services.
Funding statement
This research was partially funded by Corteva Agrisciences.
Competing interests
The authors declare no competing interests.