QTL Mapping Using a High-Density Genetic Map to Identify Candidate Genes Associated With Metribuzin Tolerance in Hexaploid Wheat (Triticum aestivum L.)

Tolerance to metribuzin, a broad-spectrum herbicide, is an important trait for weed control in wheat breeding. However, the genetics of metribuzin tolerance in relation to the underlying quantitative trait loci (QTL) and genes is limited. This study developed F₈ recombinant inbred lines (RILs) from a cross between a highly resistant genotype (Chuan Mai 25) and highly susceptible genotype (Ritchie), which were used for QTL mapping of metribuzin tolerance. Genotyping was done using a diversity arrays technology sequencing (DArTseq) platform, and phenotyping was done in controlled environments. Herbicide tolerance was measured using three traits, visual score (VS), reduction of chlorophyll content (RCC), and mean value of chlorophyll content for metribuzin-treated plants (MCC). A high-density genetic linkage map was constructed using 2,129 DArTseq markers. Inclusive composite interval mapping (ICIM) identified seven QTL, one each on chromosomes 2A, 2D, 3A, 3B, 4A, 5A, and 6A. Three major QTL—Qrcc.uwa.2AS, Qrcc.uwa.5AL, and Qrcc.uwa.6AL—explained 11.39%, 11.06%, and 11.45% of the phenotypic variation, respectively. The 5A QTL was further validated using kompetitive allele-specific PCR (KASP) assays in an F₃ validation population developed from Chuan Mai 25 × Dagger. Blasting the single-nucleotide polymorphisms (SNPs) flanking the QTL in the wheat reference genome RefV1.0 revealed SNP markers within or very close to annotated genes which could be candidate genes responsible for metribuzin tolerance. Most of the candidate genes were related to metabolic detoxification, especially those of P450 pathway and xenobiotic transmembrane transporter activity, which are reportedly key molecules responsible for herbicide tolerance. This study is the first to use specially developed populations to conduct QTL mapping on the metribuzin tolerance trait. The three major QTL and candidate genes identified in this study could facilitate marker-assisted metribuzin breeding in wheat. The QTL could be fine-mapped to locate the genes responsible for metribuzin tolerance, which could be introgressed into elite wheat cultivars.

Introduction

Wheat (Triticum aestivum L.) is a major food crop worldwide, providing approximately 20% of the global daily consumption of calories and proteins (Monneveux et al., 2012; Shiferaw et al., 2013). Weeds are among the many restrictive factors that affect wheat yield, as they compete with the crop for light, water, nutrients, and space (Song et al., 2005). Weeds are generally controlled by herbicide application (Pan et al., 2012). Breeding wheat cultivars with high tolerance to herbicides will improve weed control efficiency and hence crop production.

Metribuzin is an excellent herbicide for controlling a wide range of weeds in dryland farming systems (Kleemann and Gill, 2008; Si et al., 2009). However, most wheat cultivars lack tolerance to this broad-spectrum herbicide. Our previous study identified cultivar Chuan Mai 25 with superior metribuzin-tolerance using a pre-emergence selection method (Bhoite et al., 2018); this cultivar is a valuable resource for increasing wheat tolerance to metribuzin. Genetic studies on metribuzin tolerance are limited. Bhoite et al. (2018) conducted a quantitative trait loci (QTL) mapping study using the International Triticeae Mapping Initiative (ITMI) population (Synthetic W7984 × Opata 85). Pilcher et al. (2017) identified 169 (upregulated) and 127 (downregulated) genes in a metribuzin-tolerant and a susceptible wheat genotype that had significant differential expression in response to metribuzin. Bhoite et al. (2019) further reported that the inheritance of pre-emergent metribuzin tolerance acted as an additive-dominance model in wheat and proposed putative genes through 90k SNP array analysis of seven cultivars with contrasting responses to metribuzin. However, no studies have constructed a genetic linkage map for metribuzin tolerance in wheat using a specially developed mapping population.

The genetic mechanism of metribuzin tolerance in wheat is not fully understood. Different mechanisms for metribuzin tolerance in weeds and crops have been reported. Metribuzin tolerance in weeds is mostly target site-based, targeting the chloroplast psbA gene encoding D1 protein, which inhibits photosynthesis at photosystem II (PSII) (Pan et al., 2012). Specific point mutation in the psbA gene conferring triazine resistance is mainly reported in weeds (Powles and Yu, 2010), while only a few mutations have been found in crops, such as triazine-tolerant canola (Reith and Straus, 1987). Most of the herbicide tolerance in crops is based on non-target site mechanisms. Villarroya et al. (2000) revealed that the inheritance of metribuzin tolerance in durum wheat involves many QTL in the nucleus, but not in the cytoplasm. Kilen and Barrentine (1983) studied the inheritance of metribuzin tolerance in wild soybean and reported that the trait is controlled by alleles at the same locus as the Hm gene. The mechanism of metribuzin tolerance in narrow-leafed lupin is linked to a non-target site, and likely to be based on P450-mediated metribuzin metabolism (Pan et al., 2012). Moreover, Javid et al. (2017) identified major QTL and the underlying cytochrome P450 genes associated with metribuzin tolerance in field pea, suggesting a non-target site mechanism in the crop.

Diversity arrays technology sequencing (DArTseq) is a high-throughput genotyping method used extensively for mapping, genetic diversity, and association mapping studies (Jaccoud et al., 2001; Marone et al., 2012). Gawroński et al. (2016) reported that DArT markers effectively target gene space in the large, complex, and repetitive rye (Secale cereale L.) genome. Xia et al. (2005) revealed DArT as effective markers for high-throughput genotyping of cassava (Manihot esculenta) and its wild relatives. Marone et al. (2012) demonstrated that wheat DArT markers represent the main features of the genome and could serve as an ideal tool for identifying candidate genes. The sequence knowledge of DArT markers can provide functional meaning to those markers facilitating candidate gene identification (Marone et al., 2012; Aitken et al., 2014).

This study aimed to 1) use the specially developed F₈ RIL populations to conduct QTL mapping on metribuzin tolerance; 2) develop a high-density genetic linkage map of hexaploid wheat with DArTseq markers; 3) identify QTL for metribuzin tolerance traits, including visual score (VS), reduction of chlorophyll content (RCC), and mean value of chlorophyll content for metribuzin-treated plants [mean chlorophyll content (MCC)]; and 4) identify candidate genes associated with the QTL for metribuzin tolerance.

Materials and Methods

Plant Material

Seeds of four cultivars of Triticum aestivum L.—Chuan Mai 25, Ritchie, Dagger, and Spear—were obtained from the Australian Winter Cereals Collection. The responses of the four cultivars to metribuzin have been reported or screened among 946 genotypes, as described by Bhoite et al. (2017). Chuan Mai 25 is among the most tolerant to metribuzin, whereas Ritchie, Dagger, and Spear are susceptible to metribuzin. The four cultivars were selected to develop one mapping population (Chuan Mai 25 × Ritchie) and two validation populations (Chuan Mai 25 × Dagger and Chuan Mai 25 × Spear).

Population Development

One hundred and thirty-seven F₈ recombinant inbred lines (RILs), developed from a cross between the metribuzin tolerant cultivar Chuan Mai 25 and susceptible cultivar Ritchie (Bhoite et al., 2017), were used for QTL mapping of metribuzin tolerance. Following the QTL mapping, crosses were made between Chuan Mai 25 and two susceptible genotypes (Spear and Dagger) to develop populations to validate the phenotypic effects of the identified QTL in different genetic backgrounds. Seventy lines were randomly selected from each F₃ validation population for the validation experiment.

Plant Growth and Treatment

The F₈ RILs were grown in seedling trays with 5 × 6 cells (50 × 50 × 60 mm) filled with homogenous sand. One seed per cell was sown for each of the RILs in a greenhouse at The University of Western Australia (31°57′S, 115°47′E) where they were watered every 2 days. Three biological replicates (in different trays) with six plants in a row (within a tray) for each replicate were used for the control and treatment, respectively. The trays were sprayed with water to 100% field capacity one day before the herbicide application. The pre-emergent treatment of metribuzin (400 g ai ha⁻¹) was as described in Bhoite et al. (2018). Specifically, metribuzin was sprayed in two passes, with the cabinet sprayer calibrated to deliver a spray volume of 106.19 L ha⁻¹, 50 cm above the target, with a forward speed of 3.6 km h⁻¹ and a pressure of 200 kPa in two flat-fan nozzles (TeeJet XR11001 flat-fan, Spraying Systems Co, Wheaton, IL, USA) (Figure S1, doi: 10.6084/m9.figshare.12570113). After being sprayed, the seedlings were transferred to a phytotron set at 25/15°C (day/night), 14-h/10-h photoperiod, and 800 µmol m⁻² s⁻¹ intensity. The plants were watered every two days. About 15 days after treatment, samples were selected for the phenotypic evaluation. The same spray and glasshouse settings were used to evaluate the effects of metribuzin in the two validation populations.

Phenotyping

The chlorophyll content of the RILs was determined using a portable Minolta SPAD-502 chlorophyll meter (Spectrum Technologies, Inc., Plainfield, IL, U.S.) to evaluate the phytotoxic effects of metribuzin on wheat seedlings, with three replicate reads for each measurement. The average read from a fully unfolded leaf represented the final result, which is linearly correlated with plant chlorophyll concentration (Hamblin et al., 2014; Bhoite et al., 2018). The effects of metribuzin were assessed using three parameters: RCC, relative to the control (high RCC values representing greater susceptibility and low RCC values representing greater tolerance); MCC, for metribuzin-treated plants; and VS from 0 (100% senescence/phytotoxicity) to 10 (no senescence/no sign of phytotoxicity) to estimate leaf senescence (for further details, see Bhoite et al. (2018). The control plants remained healthy, with a VS of 10 for all plants.

Phenotypic Data Analysis

A phenotypic data analysis was conducted using GenStat statistical software 17^th edition (VSN International, 2014). Analysis of variance (ANOVA) was performed based on the fixed effects model of Y_ij = μ + g_j + Ɛ_ij, where Y_ij is the observed phenotypic mean, μ is the population mean, g_j is the effect due to the j^th genotype, and Ɛ_ij is the random error. Heritability was measured using the formula of ${\text{h}}^2={\delta }_g^2/({\delta }_g^2+{\delta }_e^2)$, where ${\delta }_g^2$ and ${\delta }_e^2$ are estimated genotypic and error variances as ${\delta }_g^2=\frac{MSg-MSe}{r}$ and ${\delta }_e^2=\frac{MSe}{r}$, respectively, where MSg is the mean square of the RILs, MSe is the residual error, and r is the number of replicates (Nyquist and Baker, 1991; Onyemaobi et al., 2018).

Genotyping and QTL Mapping

Genomic DNA was extracted from the leaves of 15-day-old seedlings of individual plants from parental lines of Chuan Mai 25, Ritchie, Dagger, and Spear, and each of the 137 F₈ mapping population and 140 F₃ validation populations (70 lines from Chuan Mai 25 × Dagger and 70 lines from Chuan Mai 25 × Spear) using the cetyl trimethyl ammonium bromide (CTAB) method with some modifications.

Total genomic DNA quality and quantity were checked on NanoDrop™ 2000 Spectrophotometer (Thermo Fisher Scientific Australia). DArT sequencing was conducted by Diversity Arrays Technology Pty Ltd. (Canberra, Australia) following the protocol described on the company’s website (https://www.diversityarrays.com/). The SNPs within the available genomic fragments were used for the linkage map construction. The DArT sequencing generated 28,668 polymorphic SNP markers. The raw data were sorted according to the marker alleles for each RIL, which were compared to that of their parental lines, Chuan Mai 25 and Ritchie. Those with the Chuan Mai 25 allele were coded as A, Ritchie as B, and heterozygous loci were coded as H. The data were then inputted into IciMapping 4.0 software, with redundant markers deleted using the “bin” function, based on the anchor information to generate the map file. After performing grouping, ordering, and rippling, a linkage map was generated.

QTL analysis was conducted for the three trait values using the inclusive composite interval mapping-additive and dominance (ICIM-ADD) method of IciMapping 4.0. The logarithm of odds (LOD) > 2.5 was selected to declare significant QTL for the VS, RCC, and MCC traits. The R² rate represents the percentage of variance explained by each QTL in proportion to the total phenotypic variance (Bearzoti and Vencovsky, 1998). The QTL were classified as major when their R² value was more than 10%. The graphical representations of QTL on linkage groups were drawn using MapChart 2.32 software (Voorrips, 2002).

QTL Validation

Kompetitive allele-specific PCR (KASP) markers developed from the flanking markers of the identified QTL were used to genotype the segregating lines in the validation populations. The KASP reagents were obtained from Geneworks Pty Ltd, Australia. In order to detect the markers, PCR was performed in a total volume of 10 μl, containing 15-ng genomic DNA as a template, 5-μl KASP 2 × reaction mix, 0.14 μl, primer mix, and DNase/RNase-free water. The PCR reactions were amplified in an Eppendorf Mastercycler programmed as 94°C for 15 min, 10 cycles of 94°C for 20 s and 61°C (with a decrement of −0.6°C per cycle until 55°C) for 60 s, and 26 cycles of 94°C for 20 s and 50°C for 60 s.

The KASP genotyping assays were read by an Applied Biosystems 7500 Fast Real-Time PCR System. Two F₃ cross populations (Chuan Mai 25 × Dagger and Chuan Mai 25 × Spear) were genotyped and phenotyped to validate a major QTL on chromosome 5A, using marker DArT1101715. The mean phenotypic performances of the lines based on the two types of allele combinations (AA and aa) were compared using the student’s t-test (P < 0.05).

Identification of Putative Candidate Genes

To identify putative coding gene regions, the flanking DArT marker sequences of the identified QTL were used to perform a BLASTN search (expect threshold = 10) against the wheat genome database Refv1.0 using the website blast tool (https://urgi.versailles.inra.fr/blast/?dbgroup=wheat_iwgsc_refseq_v1_chromosomes&program=blastn). Genes within the QTL interval with known functions for metabolic detoxification and stress resistance, especially those having SNP variation within the genes, were considered potential candidate genes for metribuzin tolerance. Gene IDs beginning with Traes were obtained from JBrowse. Stress-relevant genes were selected from the other two protein databases (https://www.uniprot.org; http://www.ebi.ac.uk/interpro/entry/IPR026509) to obtain more information, including those with protein domain, family, molecular, and biological functions of the candidate genes.

Results

Phenotypic Variation of Wheat Under Metribuzin Stress

The phenotypic data analysis revealed a normal frequency distribution among the 137 F₈ RILs of Chuan Mai 25 × Ritchie for the metribuzin response (Figure 1). A wide range of phenotypic variation was observed in the RILs for the three traits in the control and metribuzin treatment (Table S1, doi: 10.6084/m9.figshare.12570113). After being treated with metribuzin, the MCC, RCC, and VS values among the RILs ranged from 0.43–26.33, 5.70–41.47, and 0.00–10.00, respectively. Broad sense heritability for chlorophyll content was 83% in the metribuzin treatment and 60% under normal growing conditions, indicating genetic factors strongly influence chlorophyll content variation in the mapping population (Table 1).

FIGURE 1

Figure 1 Phenotypic distribution of (A) visual score, (B) reduction of chlorophyll content, and (C) mean value of chlorophyll content for metribuzin-treated plants in 137 RILs of the Chuan Mai 25 × Ritchie mapping population.

TABLE 1

Table 1 Analysis of variance for metribuzin tolerance and associated traits and their heritability estimates in the 137 RILs of Chuan Mai 25 × Ritchie mapping population.

The Chuan Mai 25 × Ritchie Map

A total of 2,129 polymorphic SNP markers were used to generate the genetic linkage map, which were distributed over 21 linkage groups and assigned to the 21 wheat chromosomes based on the anchor DArTseq markers. The constructed genetic map covered 13,599 cM of the genome, with an average distance of 6.39 cM between adjacent markers. Linkage analysis assigned 38.99%, 36.92%, and 24.10% of the markers to the A, B, and D genomes, respectively. The DArTseq markers were not uniformly distributed on different chromosomes; chromosome 4D had the fewest markers (32), and chromosome 1A had the most markers (152). Moreover, chromosome 4D had the largest marker gap (15.85 cM), and chromosome 3A had the smallest gap (3.69 cM).

QTL for Metribuzin Tolerance

The QTL for metribuzin tolerance was mapped by integrating appropriate phenotypic and genotypic data (Tables S2, S3, doi: 10.6084/m9.figshare.12570113). In Chuan Mai 25 × Ritchie, seven QTL for metribuzin tolerance were detected for VS, RCC, and MCC, with one each on chromosomes 2A, 2D, 3A, 3B, 4A, 5A, and 6A (Table 2 and Figure 2). The LOD scores ranged from 2.60–5.49, explaining 8.62–11.45% of the phenotypic variation in RCC (Table 1). The seven QTL were mostly located on long chromosome arms, apart from QRCC.UWA.2AS, QMCC.UWA.2DS, and QRCC.UWA.4AS that were located on the short chromosome arms of 2A, 2D, and 4A at 203.8 cM, 144.0 cM, and 107.3 cM, respectively. The QTL on the long chromosome arm of 6A (QVS.UWA.6AL, QRCC.UWA.6AL, and QMCC.UWA.6AL) were detected at 429 cM for all three traits (VS, RCC, and MCC), explaining up to 11.45% of the phenotypic variation, with a LOD score of 5.49. Another major QTL (QRCC.UWA.5AL) contributed 11.06% of the phenotypic variation, with a LOD score of 5.27, which was mapped on the long chromosome arm of 5A at 164 cM for RCC. The other two QTL (QRCC.UWA.3AL and QRCC.UWA.3BL) were detected on the chromosomes of 3AL and 3BL at 486.4 cM and 612 cM, respectively. Among the seven identified QTL for metribuzin tolerance, those on chromosomes 2AS, 2DS, 3BL, 4AS, 5AL, and 6AL gained favorable alleles from Chuan Mai 25, while that on 3AL gained its favorable allele from Ritchie (Table 1).

TABLE 2

Table 2 Quantitative trait loci (QTL) for metribuzin tolerance identified in the Chuan Mai 25 × Ritchie mapping population.

FIGURE 2

Figure 2 Locations of detected QTL for metribuzin tolerance in Chuan Mai 25 × Ritchie RILs for visual score (VS), reduction in chlorophyll content (RCC), and mean value of chlorophyll content for metribuzin-treated plants (MCC). QTL are indicated by solid bars; blue bars indicate QTL for RCC, red bars indicate QTL for VS and pink bars indicate QTL for MCC. Map distances are indicated on the left in Kosambi centimorgans and DArT markers are indicated on the right of each chromosome.

Validation of the Major QTL for Metribuzin Tolerance Under Herbicide Stress

A major QTL (QRCC.UWA.5AL) explained 11.06% of the phenotypic variation, with a high LOD score of 5.27 detected for RCC, and was contributed by the resistant genotype Chuan Mai 25. This QTL was selected for validation in 70 randomly selected lines from the two F₃ populations (Chuan Mai 25 × Dagger and Chuan Mai 25 × Spear). DArT1101715, located just 0.29 cM from the peak QTL location, was used to validate the 5A QTL.

The validation lines with different alleles for each DArT marker were separated into allele group 1 (susceptible allele) and allele group 2 (resistant allele). DArT1101715, targeting QRCC.UWA.5AL (Figure 3), had a 21.2% higher average VS and 27.1% higher chlorophyll content for the Chuan Mai 25 × Dagger progenies with homozygous alleles from Chuan Mai 25 (allele group 2) than Ritchie (allele group 1). However, no polymorphism occurred for the marker in the Chuan Mai 25 × Spear population. As DArT1101715 is 0.29 cM from the QTL peak, monomorphism of this marker did not necessarily mean no variance at the QRCC.UWA.5AL locus, suggesting that more markers need to be developed for this locus.

FIGURE 3

Figure 3 Genotyping plot of KASP assay for DArT1101715 on Chuan Mai 25 × Dagger. X-axis: Allele1, reported by FAM fluorescence; Y-axis, Allele 2, reported by HEX fluorescence; Red dots, homozygous allele group 1 (susceptible); Blue dots: homozygous allele group 2 (resistant); Green dots, heterozygous alleles; Black crosses: no call.

Putative Candidate Genes

From the seven QTL for metribuzin tolerance in Chuan Mai 25 × Ritchie, the flanking markers of each QTL were blasted in the URGI-JBrowse and Uniprot databases, which identified several putative candidate genes related to metabolic detoxification, including cytochrome P450 pathway, xenobiotic transmembrane transporter, and defense responses (Table 3).

TABLE 3

Table 3 Putative candidate genes for the identified QTL conferring metribuzin tolerance in wheat.

For Qrcc.uwa.2AS and Qmcc.uwa.2AS, the blastN search of the markers in Refv1.0 identified three potential candidate genes with known functions in stress tolerance, including a gene coding cytochrome P450 which regulates oxidoreductase activity. Qrcc.uwa.3AL was also co-located with three candidate genes related to the defense response, negative regulation of translation, and protein detoxification. One candidate gene was identified from Qrcc.uwa.3BL, which is involved in ion transport, regulation of jasmonic acid (JA) biosynthetic, and calcium-mediated signaling. Two important candidate genes were identified from Qrcc.uwa.4AS with known functions of F-box-like domain superfamily, transmembrane protein, and fatty acid biosynthesis. One candidate gene was identified from Qrcc.uwa.5AL, which codes auxin canalization and plant pleckstrin homology-like region, and plays an important role in plant growth and resistance. One potential gene was identified from Qvs.uwa.6AL, Qrcc.uwa.6AL, and Qmcc.uwa.6AL, which was located in the cytoplasm with a phenylalanine ammonia-lyase function involved in trans-cinnamate biosynthesis.

Discussion

This study is the first to use specially developed population lines to conduct QTL mapping on metribuzin tolerance in wheat. Seven QTL for metribuzin tolerance were identified from the F₈ RILs derived from a highly metribuzin-resistant genotype (Chuan Mai 25) and a highly susceptible genotype (Ritchie). The major 5A QTL was further validated in an F₃ validation population, developed from Chuan Mai 25 × Dagger crosses, using KASP assays. Most of the identified QTL were related to metabolic detoxification, especially the P450 pathway and xenobiotic transmembrane transporter activity.

RCC offers a rapid and non-destructive method for evaluating senescence in wheat seedlings in the metribuzin treatment as chlorophyll content and senescence are highly correlated (Bhoite et al., 2018). Metribuzin mainly inhibits photosynthesis by competing with plastoquinone at the plastoquinone-binding site on the D1 protein in the PSII complex (Pan et al., 2012). A series of reactions occur in plants in response to herbicides, which reduce chlorophyll content due to the production of active oxygen species (Xu et al., 2015). Higher RCC values in this study indicated that part of the PSII system was damaged, as reported under various stresses in plants (Gill et al., 2015; Xu et al., 2015; Wang et al., 2016).

One major QTL was located on the short arm of wheat chromosome 2A. This genomic region harbored cytochrome P450 genes, which are widely reported as responsible for tolerance to different stresses, including herbicides. Cytochrome P450 enzymes can accelerate metabolism and facilitate the detoxification of certain herbicide molecules (Zimmerlin and Durst, 1992). Bhoite et al. (2019) also reported that the candidate gene on chromosome 2A has the important biological function of metabolic detoxification, which is consistent with our findings. An increase in P450-dependent metabolism was first reported in Lolium rigidum herbicide-resistant biotypes from Australia and Alopecurus myosuroides biotypes from Europe (Holt et al., 1993). Siminszky (2006) reported plant cytochrome P450-mediated herbicide metabolism and the pivotal role of plant P450s in herbicide metabolism. Furthermore, enhanced detoxification in herbicide-resistant weeds was associated with elevated levels of P450 activity (Siminszky, 2006). Cytochrome P450 enzymes facilitate herbicide metabolism in many crops, including wheat, maize, and sorghum (Werck-Reichhart et al., 2000). Pan et al. (2012) indicated that the mechanism of metribuzin tolerance in lupin is attributed to P450-mediated metabolism. Yu and Powles (2014) revealed that herbicide resistance is often due to metabolic resistance and that cytochrome P450 monooxygenase is usually implicated in herbicide metabolic resistance. Future studies on metabolic pathways could identify the specific genes related to herbicide metabolic resistance (Fang et al., 2018; Yang et al., 2018a; Yang et al., 2018b).

A major QTL (Qrcc.uwa.3AL) contributing to variations in metribuzin tolerance was located on wheat chromosome 3A, a region known to carry major genes related to tolerance and adaptability (Kulwal et al., 2005). The candidate gene encodes a ribosome-inactivating protein related to the defense response. The results are consistent with those of Yang et al. (2016), who reported a positive correlation between the level of ribosome-inactivating protein (RIP) and resistance to cucumber mosaic virus in Momordica charantia. RIPs can induce apoptosis in a wide variety of cells (Narayanan et al., 2005), which plays a role in the regulation of cell death induced by the herbicide. Qrcc.uwa.3AL was also related to antiporter activity and xenobiotic transmembrane transporter activity to realize the protein detoxification by multidrug and toxic compound extrusion (MATE) protein. Transporter proteins from the MATE family are vital for metabolite transport in plants (He et al., 2010).

The genomic region of Qrcc.uwa.3BL harbors genes encoding ion transport proteins which operate in the ion channel activity and are related to JA regulation. JA is a plant hormone synthesized in the chloroplast membrane, which plays an important role in the abiotic stress responses (Kazan, 2015). It is also involved in a series of stress-related processes, such as stomatal closure and the detoxification activity of antioxidant enzymes, which could increase plant tolerance to various abiotic stresses (Zhou et al., 2015; Ahmad et al., 2016). Hu et al. (2019) indicated that SlF3HL is a positive regulator of chilling stress tolerance in tomato, possibly by regulating JA biosynthesis and signaling. Chapman et al. (2018) revealed that JA contributed to soybean’s tolerance to soybean aphid. Interestingly, JA triggered stress tolerance to the herbicide imazapic in tobacco (Kaya and Doganlar, 2016). JA enhanced plant tolerance to stresses in wheat by increasing antioxidant activity (Qiu et al., 2014). Thus, the 3BL QTL related to ion transport and JA pathway might play a vital role in metribuzin tolerance in wheat.

For the 4AS QTL, the candidate gene TraesCS4A01G452300 belongs to the F-box-like domain superfamily. Gupta et al. (2015) indicated that F-box genes played a crucial role in the response to biotic and abiotic stresses in chickpea. The wheat F-box protein gene TaFBA1 was related to plant tolerance to heat stress (Li et al., 2018). Various abiotic stresses induced the expression of TaFBA1, including drought, NaCl, and abscisic acid (ABA) (Zhou et al., 2014). Peng et al. (2014) demonstrated that F-box proteins conferred salinity tolerance in Arabidopsis. Moreover, the F-box protein COI1 plays an vital role in the jasmonate-signaling pathway and regulates primary carbohydrate metabolism in tobacco (Wang et al., 2014). F-box proteins also enhance polyphenol production and UV tolerance in Arabidopsis (Zhang et al., 2015). Another candidate gene TraesCS4A01G007400 in the QTL region coding the important 3-ketoacyl-CoA synthase (KCS) in the fatty acid biosynthesis. KCS is a key enzyme for the synthesis of long-chain fatty acids and wax, which serve as the first line of defense against pathogens, phytophagous insects, and environmental stresses, such as drought, UV damage, and frost (Todd et al., 1999; Joubès et al., 2008). Todd et al. (1999) demonstrated that KCS1 plays a vital role in the resistance of Arabidopsis thaliana to low humidity stress at a young age. Azachi et al. (2002) cloned a KCS gene from Dunaliella salina, which has outstanding salt tolerance, that might play a role in balancing adverse external osmotic pressure. Fan et al. (2018) found that stress conditions dramatically increased the expression of the MaKCS gene in Mychonastes afer.

The genomic region of major QTL Qrcc.uwa.5AL harbored genes coding auxin canalization and related to the plant pleckstrin homology-like region. Lee et al. (2014) demonstrated the canalization-based vein formation in a growing leaf and cell division along the shortest axis contributed to the relaxation of stress, as reported by Alim et al. (2012). Auxin acts as a general coordinator of plant growth and development (Leyser, 2018). Therefore, auxin canalization may be crucial for herbicide transportation and cell metabolic detoxification.

The candidate gene for the major 6AL QTL was related to coding aromatic amino acid lyase involved in trans-cinnamate biosynthesis. The aromatic amino acid has been related to the detoxification of glyphosate (Samsel and Seneff, 2013), and regulation of salt stress in wheat (Jiang et al., 2017). One flanking DArT marker of the QTL overlapped the mRNA CJ661337 (https://www.ncbi.nlm.nih.gov/nuccore/93059904), which is involved in the stress response (Mochida et al., 2006), suggesting the important role this locus may play in response to metribuzin stress.

In comparison to the study done by Bhoite et al. (2018), who reported metribuzin tolerance QTL on chromosomes 1A, 2D, and 4A, we detected six new loci in this study. Blasting the QTL flanking markers with the wheat reference genome RefV1.0 revealed that the 2D QTL overlapped in the two studies, as a flanking marker “DArT1219191” of QMCC.UWA.2DS fell into the interval of Qsns.uwa.2DS. However, no candidate genes were suggested for this QTL as the physical position of “DArT3064546”, another flanking marker of QMCC.UWA.2DS, was too far away, making the physical interval of the QTL too large (617,964,302 bp) for identifying any possible candidate genes. We suggest that the difference between the Chinese Spring (reference genome) and the parental cultivars used in this mapping study might be due to chromosomal fragment translocations. Although both studies reported QTL on chromosome 4A, their physical positions were far from each other, suggesting that chromosome 4A might harbor many genes responsible for metribuzin tolerance. Bhoite et al. (2019) later reported putative genes for metribuzin tolerance on chromosomes 2A, 2D, 3B, 4A, 4B, 7A, 7B, and 7D. We also identified QTL with candidate genes on chromosomes 2A, 3B, and 4A, but at different physical positions. Shi et al. (2020) investigated genomic regions underlying the tolerance to two other herbicides (carfentrazone-ethyl and tribenuron-methyl) and detected QTL on chromosomes 1B, 2D, 5B, 5D, 6D, and 7D in wheat. No overlapping was found with the QTL identified in this study, suggesting different regulatory mechanisms for the tolerance to different herbicides.

Controlling weeds effectively is vital for wheat production, especially in areas with a Mediterranean climate, such as Australia. Herbicides are the main management method for controlling weeds in crops. Thus, understanding the mechanisms of herbicide tolerance and breeding varieties with high herbicide resistance is critical for future economic food production. This study revealed that chlorophyll traits can be used to select for herbicide tolerance in wheat. The identified SNP markers and candidate genes will be useful for marker-assisted selection of metribuzin tolerance, which could be an efficient way to control weeds and further improve wheat production.

Conclusion

A high-density genetic linkage map was constructed using 2,129 DArTseq markers. ICIM identified seven QTL, one each on chromosomes 2A, 2D, 3A, 3B, 4A, 5A, and 6A. One major QTL for metribuzin tolerance on chromosome 5A was further validated using KASP assays in the F₃ validation population developed from a Chuanmai 25 × Dagger cross. Three major QTL (Qrcc.uwa.2AS, Qrcc.uwa.5AL, and Qrcc.uwa.6AL) explained a total of 33.90% of the phenotypic variation, and the identified candidate genes could facilitate marker-assisted selection of metribuzin tolerance in wheat. Moreover, the QTL could be fine-mapped to locate the causal genes responsible for metribuzin tolerance for sustainable wheat production.

Data Availability Statement

QTL data is available in figshare data repository portal, with a doi: 10.6084/m9.figshare.12570113.

Author Contributions

LX and HL contributed equally to this work. LX, HL, PS, and GY designed and conceived the study. LX, HL, GL, JW, and RB performed the experiments. LX, HL, and AK analyzed the data. LX, HL, and WZ wrote the manuscript.

Funding

This work was supported by the Global Innovation Linkage program (GIL53853) from the Australian Department of Industry, Innovation and Science, and the China Scholarship Council (201808330058).

Conflict of Interest

Author AK was employed by company Diversity Arrays Technology Pty. Ltd.

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

The reviewer HA declared a past co-authorship with several of the authors HL and GY to the handling editor.

Acknowledgments

The authors would like to thank Mechelle Owen at the UWA Australian Herbicide Resistance Initiative and UWA glasshouse staff Robert Creasy and Bill Piasini for their help with the experimental setup. The authors also thank Christine Davies of Tweak Editing for correcting English language errors.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.573439/full#supplementary-material

References

Ahmad, P., Rasool, S., Gul, A., Sheikh, S. A., Akram, N. A., Ashraf, M., et al. (2016). Jasmonates: multifunctional roles in stress tolerance. Front. Plant Sci. 7, 813. doi: 10.3389/fpls.2016.00813

PubMed Abstract | CrossRef Full Text | Google Scholar

Aitken, K. S., McNeil, M. D., Berkman, P. J., Hermann, S., Kilian, A., Bundock, P. C., et al. (2014). Comparative mapping in the Poaceae family reveals translocations in the complex polyploid genome of sugarcane. BMC Plant Biol. 14, 190. doi: 10.1186/s12870-014-0190-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Alim, K., Hamant, O., Boudaoud, A. (2012). Regulatory role of cell division rules on tissue growth heterogeneity. Front. Plant Sci. 3, 174. doi: 10.3389/fpls.2012.00174

PubMed Abstract | CrossRef Full Text | Google Scholar

Azachi, M., Sadka, A., Fisher, M., Goldshlag, P., Gokhman, I., Zamir, A. (2002). Salt induction of fatty acid elongase and membrane lipid modifications in the extreme halotolerant alga Dunaliella salina. Plant Physiol. 129, 1320–1329. doi: 10.1104/pp.001909

CrossRef Full Text | Google Scholar

Bearzoti, E., Vencovsky, R. (1998). Estimation of the proportion of genetic variance explained by molecular markers. Genet. Mol. Biol. 21, 4. doi: 10.1590/S1415-47571998000400025

CrossRef Full Text | Google Scholar

Bhoite, R. N., Si, P., Stefanova, K. T., Siddique, K. H. M., Yan, G. (2017). Identification of new metribuzin-tolerant wheat (Triticum spp.) genotypes. Crop Pasture Sci. 68, 401–408. doi: 10.1071/cp17017

CrossRef Full Text | Google Scholar

Bhoite, R., Onyemaobi, I., Si, P., Siddique, K. H. M., Yan, G. (2018). Identification and validation of QTL and their associated genes for pre-emergent metribuzin tolerance in hexaploid wheat (Triticum aestivum L.). BMC Genet. 19, 102. doi: 10.1186/s12863-018-0690-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhoite, R., Si, P., Liu, H., Xu, L., Siddique, K. H. M., Yan, G. J. (2019). Inheritance of pre-emergent metribuzin tolerance and putative gene discovery through high-throughput SNP array in wheat (Triticum aestivum L.). BMC Plant Biol. 19, 457. doi: 10.1186/s12870-019-2070-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chapman, K. M., Marchi-Werle, L., Hunt, T. E., Heng-Moss, T. M., Louis, J. (2018). Abscisic and jasmonic acids contribute to soybean tolerance to the soybean aphid (Aphis glycines Matsumura). Sci. Rep. 8, 15148. doi: 10.1038/s41598-018-33477-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Yuan, C., Jin, Y., Hu, G. R., Li, F. L. (2018). Characterization of 3-ketoacyl-coA synthase in a nervonic acid producing oleaginous microalgae Mychonastes afer. Algal Res. 31, 225–231. doi: 10.1016/j.algal.2018.02.017

CrossRef Full Text | Google Scholar

Fang, Y. M., Hou, Z. N., Zhang, X. D., Yang, D. F., Liang, Z. S. (2018). Diverse specialized metabolism and their responses to lactalbumin hydrolysate in hairy root cultures of Salvia miltiorrhiza Bunge and Salvia castanea Diels f. tomentosa Stib. Biochem. Eng. J. 131, 58–69. doi: 10.1016/j.bej.2017.12.007

CrossRef Full Text | Google Scholar

Gawroński, P., Pawełkowicz, M., Tofil, K., Uszyński, G., Sharifova, S., Ahluwalia, S., et al. (2016). DarT markers effectively target gene space in the rye genome. Front. Plant Sci. 7, 1600. doi: 10.3389/fpls.2016.01600

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, R. A., Zang, L. L., Ali, B., Farooq, M. A., Cui, P., Yang, S., et al. (2015). Chromium-induced physio-chemical and ultrastructural changes in four cultivars of Brassica napus L. Chemosphere 120, 154–164. doi: 10.1016/j.chemosphere.2014.06.029

CrossRef Full Text | Google Scholar

Gupta, S., Garg, V., Kant, C., Bhatia, S. (2015). Genome-wide survey and expression analysis of F-box genes in chickpea. BMC Genomics 16, 67. doi: 10.1186/s12864-015-1293-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamblin, J., Stefanova, K., Angessa, T. T. (2014). Variation in chlorophyll content per unit leaf area in spring wheat and implications for selection in segregating material. PloS One 9, e92529. doi: 10.1371/journal.pone.0092529

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Szewczyk, P., Karyakin, A., Evin, M., Hong, W. X., Zhang, Q. H., et al. (2010). Structure of a cation-bound multidrug and toxic compound extrusion transporter. Nature 467, 991–994. doi: 10.2210/pdb3mku/pdb

CrossRef Full Text | Google Scholar

Holt, J. S., Powles, S. B., Holtum, J. A. M. (1993). Mechanisms and agronomic aspects of herbicide resistance. Annu. Rev. Plant Physiol. Plant Mol. Biol. 44, 203–229. doi: 10.2134/jae1991.0063

CrossRef Full Text | Google Scholar

Hu, T., Wang, Y., Wang, Q., Dang, N., Wang, L., Liu, C., et al. (2019). The tomato 2-oxoglutarate-dependent dioxygenase gene SlF3HL is critical for chilling stress tolerance. Hortic. Res. 6, 45. doi: 10.1038/s41438-019-0127-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaccoud, D., Peng, K., Feinstein, D., Kilian, A. (2001). Diversity arrays: a solid state technology for sequence information independent genotyping. Nucleic Acids Res. 29, E25. doi: 10.1093/nar/29.4.e25

PubMed Abstract | CrossRef Full Text | Google Scholar

Javid, M., Noy, D., Sudheesh, S., Forster, J., Kaur, S. (2017). Identification of QTLs associated with metribuzin tolerance in field pea (Pisum sativum L.). Euphytica 213, 91. doi: 10.1007/s10681-017-1878-4

CrossRef Full Text | Google Scholar

Jiang, Q. Y., Niu, F. Y., Sun, X. J., Hu, Z., Li, X. H., Ma, Y. Z., et al. (2017). RNA-seq analysis of unintended effects in transgenic wheat overexpressing the transcription factor GmDREB1. Crop J. 5, 207–218. doi: 10.1016/j.cj.2016.12.001

CrossRef Full Text | Google Scholar

Joubès, J., Raffaele, S., Bourdenx, B., Garcia, C., Laroche-Traineau, J., Moreau, P., et al. (2008). The VLCFA elongase gene family in Arabidopsis thaliana: phylogenetic analysis, 3D modelling and expression profiling. Plant Mol. Biol. 67, 547–566. doi: 10.1007/s11103-008-9339-z

CrossRef Full Text | Google Scholar

Kaya, A., Doganlar, Z. B. (2016). Exogenous jasmonic acid induces stress tolerance in tobacco (Nicotiana tabacum) exposed to imazapic. Ecotoxicol. Environ. Saf. 124, 470–479. doi: 10.1016/j.ecoenv.2015.11.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazan, K. (2015). Diverse roles of jasmonates and ethylene in abiotic stress tolerance. Trends Plant Sci. 20, 219–229. doi: 10.1016/j.tplants.2015.02.001

CrossRef Full Text | Google Scholar

Kilen, T. C., Barrentime, W. L. (1983). Linkage relationships in soybean between genes controlling reactions to phytophthora rot and metribuzin. Crop Sci. 23, 894–896. doi: 10.2135/cropsci1983.0011183x002300050019x

CrossRef Full Text | Google Scholar

Kleemann, S. G. L., Gill, G. S. (2008). Applications of metribuzin for the control of rigid brome (Bromus rigidus) in no-till barley crops of southern Australia. Weed Technol. 22, 34–37. doi: 10.1614/wt-07-017.1

CrossRef Full Text | Google Scholar

Kulwal, P. L., Kumar, N., Gaur, A., Khurana, P., Khurana, J. P., Tyagi, A. K., et al. (2005). Mapping of a major QTL for pre-harvest sprouting tolerance on chromosome 3A in bread wheat. Theor. Appl. Genet. 111, 1052–1059. doi: 10.1007/s00122-005-0021-4

CrossRef Full Text | Google Scholar

Lee, S. W., Feugier, F. G., Morishita, Y. (2014). Canalization-based vein formation in a growing leaf. J. Theor. Biol. 353, 104–120. doi: 10.1016/j.jtbi.2014.03.005

CrossRef Full Text | Google Scholar

Leyser, O. (2018). Auxin signaling. Plant Physiol. 176, 465–479. doi: 10.1104/pp.17.00765

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q. X., Wang, W. Q., Wang, W. L., Zhang, G. Q., Liu, Y., Wang, Y., et al. (2018). Wheat F-box protein gene TaFBA1 is involved in plant tolerance to heat stress. Front. Plant Sci. 9, 521. doi: 10.3389/fpls.2018.00521

PubMed Abstract | CrossRef Full Text | Google Scholar

Marone, D., Panio, G., Ficco, D. B., Russo, M. A., De Vita, P., Papa, R., et al. (2012). Characterization of wheat DArT markers: genetic and functional features. Mol. Genet. Genomics 287, 741–753. doi: 10.1007/s00438-012-0714-8

CrossRef Full Text | Google Scholar

Mochida, K., Kawaura, K., Shimosaka, E., Kawakami, N., Shin, I. T., Kohara, Y., et al. (2006). Tissue expression map of a large number of expressed sequence tags and its application to in silico screening of stress response genes in common wheat. Mol. Genet. Genomics 276 (3), 304–312. doi: 10.1007/s00438-006-0120-1

CrossRef Full Text | Google Scholar

Monneveux, P., Jing, R., Misra, S. C. (2012). Phenotyping for drought adaptation in wheat using physiological traits. Front. Physiol. 3, 429. doi: 10.3389/fphys.2012.00429

PubMed Abstract | CrossRef Full Text | Google Scholar

Narayanan, S., Surendranath, K., Bora, N., Surolia, A., Karande, A. A. (2005). Ribosome inactivating proteins and apoptosis. FEBS Lett. 579, 1324–1331. doi: 10.1016/j.febslet.2005.01.038

CrossRef Full Text | Google Scholar

Nyquist, W. E., Baker, R. J. (1991). Estimation of heritability and prediction of selection response in plant populations. Crit. Rev. Plant Sci. 10, 235–322. doi: 10.1080/07352689109382313

CrossRef Full Text | Google Scholar

Onyemaobi, I., Ayalew, H., Liu, H., Siddique, K. H. M., Yan, G. (2018). Identification and validation of a major chromosome region for high grain number per spike under meiotic stage water stress in wheat (Triticum aestivum L.). PloS One 13, e0194075. doi: 10.1371/journal.pone.0194075

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, G., Si, P., Yu, Q., Tu, J. M., Powles, S. (2012). Non-target site mechanism of metribuzin tolerance in induced tolerant mutants of narrow-leafed lupin (Lupinus angustifolius L.). Crop Pasture Sci. 63, 452–458. doi: 10.1071/CP12065

CrossRef Full Text | Google Scholar

Peng, J., Li, Z., Wen, X., Li, W., Shi, H., Yang, L., et al. (2014). Salt-induced stabilization of EIN3/EIL1 confers salinity tolerance by deterring ROS accumulation in Arabidopsis. PloS Genet. 10, e1004664. doi: 10.1371/journal.pgen.1004664

PubMed Abstract | CrossRef Full Text | Google Scholar

Pilcher, W., Zandkamiri, H., Arceneaux, K., Harrison, S., Baisakh, N. (2017). Genome-wide microarray analysis leads to identification of genes in response to herbicide, metribuzin in wheat leaves. PloS One 12, e0189639. doi: 10.1371/journal.pone.0189639

PubMed Abstract | CrossRef Full Text | Google Scholar

Powles, S. B., Yu, Q. (2010). Evolution in action: plants resistant to herbicides. Annu. Rev. Plant Biol. 61, 317–347. doi: 10.1146/annurev-arplant-042809-112119

CrossRef Full Text | Google Scholar

Qiu, Z., Guo, J., Zhu, A., Zhang, L., Zhang, M. (2014). Exogenous jasmonic acid can enhance tolerance of wheat seedlings to salt stress. Ecotoxicol. Environ. Saf. 104, 202–208. doi: 10.1016/j.ecoenv.2014.03.014

CrossRef Full Text | Google Scholar

Reith, M. E., Straus, N. A. (1987). Nucleotide sequence of the chloroplast gene responsible for triazine resistance in canola. Theor. Appl. Genet. 73, 357–363. doi: 10.1007/bf00262501

CrossRef Full Text | Google Scholar

Samsel, A., Seneff, S. (2013). Glyphosate’s suppression of cytochrome P450 enzymes and amino acid biosynthesis by the gut microbiome: pathways to modern diseases. Entropy 15, 1416–1463. doi: 10.3390/e15041416

CrossRef Full Text | Google Scholar

Shi, C. N., Zheng, Y. T., Geng, J. Y., Liu, C. Y., Pei, H., Ren, Y., et al. (2020). Identification of herbicide resistance loci using a genome-wide association study and linkage mapping in Chinese common wheat. Crop J. 8, 666–675. doi: 10.1016/j.cj.2020.02.004

CrossRef Full Text | Google Scholar

Shiferaw, B., Smale, M., Braun, H.-J., Duveiller, E., Reynolds, M., Muricho, G. (2013). Crops that feed the world 10. Past successes and future challenges to the role played by wheat in global food security. Food Secur. 5, 291–317. doi: 10.1007/s12571-013-0263-y

CrossRef Full Text | Google Scholar

Si, P., Buirchell, B., Sweetingham, M. W. (2009). Improved metribuzin tolerance in narrow-leafed lupin (Lupinus angustifolius L.) by induced mutation and field selection. Field Crops Res. 113, 282–286. doi: 10.1016/j.fcr.2009.06.003

CrossRef Full Text | Google Scholar

Siminszky, B. (2006). Plant cytochrome P450-mediated herbicide metabolism. Phytochem. Rev. 5, 445–458. doi: 10.1007/s11101-006-9011-7

CrossRef Full Text | Google Scholar

Song, W. J., Zhou, W. J., Jin, Z. L., Cao, D. D., Joel, D. M., Takeuchi, Y., et al. (2005). Germination response of Orobanche seeds subjected to conditioning temperature, water potential and growth regulator treatments. Weed Res. 45, 467–476. doi: 10.1111/j.1365-3180.2005.00477.x

CrossRef Full Text | Google Scholar

Todd, J., Post-Beittenmiller, D., Jaworski, J. G. (1999). KCS1 encodes a fatty acid elongase 3-ketoacyl-CoA synthase affecting wax biosynthesis in Arabidopsis thaliana. Plant J. 17, 119–130. doi: 10.1046/j.1365-313x.1999.00352.x

CrossRef Full Text | Google Scholar

Villarroya, M., Escorial, M., Garcia-Baudin, J., Chueca, M. (2000). Inheritance of tolerance to metribuzin in durum wheat. Weed Res. 40, 293–300. doi: 10.1046/j.1365-3180.2000.00188.x

CrossRef Full Text | Google Scholar

Voorrips, R. E. (2002). MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 93 (1), 77–78. doi: 10.1093/jhered/93.1.77

CrossRef Full Text | Google Scholar

VSN International (2014). Genstat for windows 17th edition (Hemel Hempstead, UK: VSN International).

Google Scholar

Wang, W., Liu, G., Niu, H., Timko, M. P., Zhang, H. (2014). The F-box protein COI1 functions upstream of MYB305 to regulate primary carbohydrate metabolism in tobacco (Nicotiana tabacum L. cv. TN90). J. Exp. Bot. 65, 2147–2160. doi: 10.1093/jxb/eru084

CrossRef Full Text | Google Scholar

Wang, J., Lv, M., Islam, F., Gill, R. A., Yang, C., Ali, B., et al. (2016). Salicylic acid mediates antioxidant defense system and ABA pathway related gene expression in Oryza sativa against quinclorac toxicity. Ecotoxicol. Environ. Saf. 133, 146–156. doi: 10.1016/j.ecoenv.2016.07.002

CrossRef Full Text | Google Scholar

Werck-Reichhart, D., Hehn, A., Didierjean, L. (2000). Cytochromes P450 for engineering herbicide tolerance. Trends Plant Sci. 5, 116–123. doi: 10.1016/s1360-1385(00)01567-3

CrossRef Full Text | Google Scholar

Xia, L., Peng, K., Yang, S., Wenzl, P., de Vicente, M. C., Fregene, M., et al. (2005). DArT for high-throughput genotyping of Cassava (Manihot esculenta) and its wild relatives. Theor. Appl. Genet. 110, 1092–1098. doi: 10.1007/s00122-005-1937-4

CrossRef Full Text | Google Scholar

Xu, L., Zhang, W. F., Ali, B., Islam, F., Zhu, J. W., Zhou, W. J. (2015). Synergism of herbicide toxicity by 5-aminolevulinic acid is related to physiological and ultra-structural disorders in crickweed (Malachium aquaticum L.). Pestic. Biochem. Physiol. 125, 53–61. doi: 10.1016/j.pestbp.2015.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, T., Meng, Y., Chen, L. J., Lin, H. H., Xi, D. H. (2016). The roles of alpha-momorcharin and jasmonic acid in modulating the response of Momordica charantia to Cucumber Mosaic Virus. Front. Microbiol. 7, 1796. doi: 10.3389/fmicb.2016.01796

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Fang, Y., Xia, P., Zhang, X., Liang, Z. (2018a). Diverse responses of tanshinone biosynthesis to biotic and abiotic elicitors in hairy root cultures of Salvia miltiorrhiza and Salvia castanea Diels f.tomentosa. Gene 643, 61–67. doi: 10.1016/j.gene.2017.11.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Huang, Z., Jin, W., Xia, P., Jia, Q., Yang, Z., et al. (2018b). DNA methylation: a new regulator of phenolic acids biosynthesis in Salvia miltiorrhiza. Ind. Crop Prod. 124, 402–411. doi: 10.1016/j.indcrop.2018.07.046

CrossRef Full Text | Google Scholar

Yu, Q., Powles, S. (2014). Metabolism-based herbicide resistance and cross-resistance in crop weeds: a threat to herbicide sustainability and global crop production. Plant Physiol. 166, 1106–1118. doi: 10.1104/pp.114.242750

CrossRef Full Text | Google Scholar

Zhang, X. B., Gou, M. Y., Guo, C. R., Yang, H. J., Liu, C. J. (2015). Down-regulation of the kelch domain-containing F-box protein 1 in Arabidopsis enhances the production of (poly)phenols and tolerance to ultraviolet radiation. Plant Physiol. 167, 337–350. doi: 10.1104/pp.114.249136

CrossRef Full Text | Google Scholar

Zhou, S., Sun, X., Yin, S., Kong, X., Zhou, S., Xu, Y., et al. (2014). The role of the F-box gene TaFBA1 from wheat (Triticum aestivum L.) in drought tolerance. Plant Physiol. Biochem. 84, 213–223. doi: 10.1016/j.plaphy.2014.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., Wu, Y., Yang, Y., Du, M., Zhang, X., Guo, Y., et al. (2015). An arabidopsis plasma membrane proton ATPase modulates JA signaling and is exploited by the Pseudomonas syringae effector protein AvrB for stomatal invasion. Plant Cell 27, 2032–2041. doi: 10.1105/tpc.15.00466

CrossRef Full Text | Google Scholar

Zimmerlin, A., Durst, F. (1992). Aryl hydroxylation of the herbicide diclofop by a wheat cytochrome-P-450 monooxygenase–substrate specificity and physiological activity. Plant Physiol. 100, 874–881. doi: 10.1104/pp.100.2.874

CrossRef Full Text | Google Scholar